Abstract: The evaluation of nucleation rates from molecular dynamics trajectories is hampered by the slow nucleation time scale and impact of finite size effects. Here, we show that accurate nucleation rates can be obtained in a very general fashion relying only on the free energy barrier, transition state theory, and a simple dynamical correction for diffusive recrossing. In this setup, the time scale problem is overcome by using enhanced sampling methods, in casu metadynamics, whereas the impact of finite size effects can be naturally circumvented by reconstructing the free energy surface from an appropriate ensemble. Approximations from classical nucleation theory are avoided. We demonstrate the accuracy of the approach by calculating macroscopic rates of droplet nucleation from argon vapor, spanning 16 orders of magnitude and in excellent agreement with literature results, all from simulations of very small (512 atom) systems.

Abstract: In order to increase the accuracy of the linearized augmented plane wave (LAPW) method, we present a new approach where the plane wave basis function is augmented by two different atomic radial components constructed at two different linearization energies corresponding to two different electron bands (or energy windows). We demonstrate that this case can be reduced to the standard treatment within the LAPW paradigm where the usual basis set is enriched by the basis functions of the tight binding type, which go to zero with zero derivative at the sphere boundary. We show that the task is closely related with the problem of extended core states which is currently solved by applying the LAPW method with local orbitals (LAPW+LO). In comparison with LAPW+LO, the number of supplemented basis functions in our approach is doubled, which opens up a new channel for the extension of the LAPW and LAPW+LO basis sets. The appearance of new supplemented basis functions absent in the LAPW+LO treatment is closely related with the existence of the ul-component in the canonical LAPW method. We discuss properties of additional tight binding basis functions and apply the extended basis set for computation of electron energy bands of lanthanum (face and body centered structures) and hexagonal close packed lattice of cadmium. We demonstrate that the new treatment gives lower total energies in comparison with both canonical LAPW and LAPW+LO, with the energy difference more pronounced for intermediate and poor LAPW basis sets.

Abstract: Superexchange interactions in alkali fullerides AC(60) are derived for C-60 molecular ions separated by interstitial alkali-metal ions. We use a multiconfiguration approach which comprises the lowest molecular orbital states of the C-60 molecule and the excited s and d states of the alkali-metal atom A. Interactions are described by the valence bond (Heitler-London) method for a complex (C-60 – A – C-60) – with two valence electrons. The electronic charge transfer between the alkali-metal atom and a neighboring C-60 molecule is not complete. The occupation probability of excited d and s states of the alkali atom is not negligible. In correspondence with the relative positions of the C-60 molecules and A atoms in the polymer crystal, we consider 180degrees and 90degrees (angle) superexchange pathways. For the former case the ground state is found to be a spin singlet separated from a triplet at similar to20 K. For T < 20 K there appear strong spin correlations for the 180degrees superexchange pathway. The results are related to spin lattice relaxation experiments on CsC60 in the polymerized and in the quenched cubic phase. (C) 2005 American Institute of Physics.

Abstract: In this work, motivated by the fabrication of an AlSb monolayer, we have focused on the electronic, mechanical and optical properties of AlSb and InSb monolayers with double-layer honeycomb structures, employing the density functional theory approach. The phonon band structure and cohesive energy confirm the stability of the XSb (X = Al and In) monolayers. The mechanical properties reveal that the XSb monolayers have a brittle nature. Using the GGA + SOC (HSE + SOC) functionals, the bandgap of the AlSb monolayer is predicted to be direct, while InSb has a metallic character using both functionals. We find that XSb (X = Al, In) two-dimensional bodies can absorb ultraviolet light. The present findings suggest several applications of AlSb and InSb monolayers in novel optical and electronic usages.

Abstract: Inspired by the recent successful growth of Ti2C and Ti3C2 monolayers, here, we investigate the structural, electronic, and mechanical properties of functionalized Ti2C and Ti3C2 monolayers by means of density functional theory calculations. The results reveal that monolayers of Ti2C and Ti3C2 are dynamically stable metals. Phonon band dispersion calculations demonstrate that two-surface functionalization of Ti2C and Ti(3)C(2)via chalcogenides (S, Se, and Te), halides (F, Cl, Br, and I), and oxygen atoms results in dynamically stable novel functionalized monolayer materials. Electronic band dispersions and density of states calculations reveal that all functionalized monolayer structures preserve the metallic nature of both Ti2C and Ti3C2 except Ti2C-O-2, which possesses the behavior of an indirect semiconductor via full-surface oxygen passivation. In addition, it is shown that although halide passivated Ti3C2 structures are still metallic, there exist multiple Dirac-like cones around the Fermi energy level, which indicates that semi-metallic behavior can be obtained upon external effects by tuning the energy of the Dirac cones. In addition, the computed linear-elastic parameters prove that functionalization is a powerful tool in tuning the mechanical properties of stiff monolayers of bare Ti2C and Ti3C2. Our study discloses that the electronic and structural properties of Ti2C and Ti3C2 MXene monolayers are suitable for surface modification, which is highly desirable for material property engineering and device integration.

Abstract: Research progress on single layer group III monochalcogenides has been increasing rapidly owing to their interesting physics. Herein, we investigate the dynamically stable single layer forms of XBi (X = Ge, Si or Sn) using density functional theory calculations. Phonon band dispersion calculations and ab initio molecular dynamics simulations reveal the dynamical and thermal stability of the considered monolayers. Raman spectra calculations indicate the existence of 5 Raman active phonon modes, 3 of which are prominent and can be observed in possible Raman measurements. The electronic band structures of the XBi single layers were investigated with and without the effects of spin-orbit coupling (SOC). Our results show that XBi single layers show semiconducting properties with narrow band gap values without SOC. However, only single layer SiBi is an indirect band gap semiconductor, while GeBi and SnBi exhibit metallic behaviors when adding spin-orbit coupling effects. In addition, the calculated linear elastic parameters indicate the soft nature of the predicted monolayers. Moreover, our predictions for the thermoelectric properties of single layer XBi reveal that SiBi is a good thermoelectric material with increasing temperature. Overall, it is proposed that single layer XBi structures can be alternative, stable 2D single layers with varying electronic and thermoelectric properties.

Abstract: While spectroscopic data on small hydrocarbons in interstellar media in combination with crossed molecular beam (CMB) experiments have provided a wealth of information on astrochemically relevant species, much of the underlying mechanistic pathways of their formation remain elusive. Therefore, in this work, the chemical reaction mechanisms of C(³P_J) + C₆H₆and C⁺(²P) + C₆H₆systems using the quantum mechanical molecular dynamics (QMMD) technique at the PBE0-D3(BJ) level of theory is investigated, mimicking a CMB experiment. Both the dynamics of the reactions as well as the electronic structure for the purpose of the reaction network are evaluated. The method is validated for the first reaction by comparison to the available experimental data. The reaction scheme for the C(³P_J) + C₆H₆system covers the literature data,<italic>e.g.</italic>the major products are the 1,2-didehydrocycloheptatrienyl radical (C₇H₅) and benzocyclopropenyl radical (C₆H₅–CH), and it reveals the existence of less common pathways for the first time. The chemistry of the C⁺(²P_J) + C₆H₆system is found to be much richer, and we have found that this is because of more exothermic reactions in this system in comparison to those in the C(³P_J) + C₆H₆system. Moreover, using the QMMD simulation, a number of reaction paths have been revealed that produce three distinct classes of reaction products with different ring sizes. All in all, at all the collision energies and orientations, the major product is the heptagon molecular ion for the ionic system. It is also revealed that the collision orientation has a dominant effect on the reaction products in both systems, while the collision energy mostly affects the charged system. These simulations both prove the applicability of this approach to simulate crossed molecular beams, and provide fundamental information on reactions relevant for the interstellar medium.

Abstract: Using first-principles calculations, we study the effect of embedding various atoms into the nanopore sites of both C6N6 and C6N8 monolayers. Our results indicate that the embedded atoms significantly affect the electronic and magnetic properties of C6N6 and C6N8 monolayers and lead to extraordinary and multifarious electronic properties, such as metallic, half-metallic, spin-glass semiconductor and dilute-magnetic semiconductor behaviour. Our results reveal that the H atom concentration dramatically affects the C6N6 monolayer. On increasing the H coverage, the impurity states also increase due to H atoms around the Fermi-level. C6N6 shows metallic character when the H atom concentration reaches 6.25%. Moreover, the effect of charge on the electronic properties of both Cr@C6N6 and C@C6N8 is also studied. Cr@C6N6 is a ferromagnetic metal with a magnetic moment of 2.40 mu(B), and when 0.2 electrons are added and removed, it remains a ferromagnetic metal with a magnetic moment of 2.57 and 2.77 mu(B), respectively. Interestingly, one can observe a semi-metal, in which the VBM and CBM in both spin channels touch each other near the Fermi-level. C@C6N8 is a semiconductor with a nontrivial band gap. When 0.2 electrons are removed, it remains metallic, and under excess electronic charge, it exhibits half-metallic behaviour.

Abstract: Using first-principles calculations the effect of topological defects, vacancies, Stone-Wales and anti-site and substitution of atoms, on the structure and electronic properties of monolayer C3N are investigated. Vacancy defects introduce localized states near the Fermi level and a local magnetic moment. While pristine C3N is an indirect semiconductor with a 0.4 eV band gap, with substitution of O, S and Si atoms for C, it remains a semiconductor with a band gap in the range 0.25-0.75 eV, while it turns into a metal with H, Cl, B, P, Li, Na, K, Be and Mg substitution. With F substitution, it becomes a dilute-magnetic semiconductor, while with Ca substitution it is a ferromagnetic-metal. When replacing the N host atom, C3N turns into: a metal (H, O, S, C, Si, P, Li and Be), ferromagnetic-metal (Mg), half-metal (Ca) and spin-glass semiconductor (Na and K). Moreover, the effects of charging and strain on the electronic properties of Na atom substitution in C3N are investigated. We found that the magnetic moment decreases or increases depending on the type and size of strain (tensile or compression). Our study shows how the band gap and magnetism in monolayer C3N can be tuned by introducing defects and atom substitution. The so engineered C3N can be a good candidate for future low dimensional devices.

Abstract: Defects are inevitably present in materials, and their existence in a material strongly affects its fundamental physical properties. We have systematically investigated the effects of surface adsorption, substitutional impurities, defect engineering, an electric field and strain engineering on the structural, electronic and magnetic properties of antimonene nanosheets, using spin-polarized density functional calculations based on first-principles. The adsorption or substitution of atoms can locally modify the atomic and electronic structures as well as induce a variety of electronic behaviors including metal, half-metal, ferromagnetic metal, dilute magnetic semiconductor and spin-glass semiconductor. Our calculations show that the presence of typical defects (vacancies and Stone-Wales defect) in antimonene affects the geometrical symmetry as well as the band gap in the electronic band structure and induces magnetism to antimonene. Moreover, by applying an external electric field and strain (uniaxial and biaxial), the electronic structure of antimonene can be easily modified. The calculation results presented in this paper provide a fundamental insight into the tunable nature of the electronic properties of antimonene, supporting its promise for use in future applications.

Abstract: Perovskite oxides are active room-temperature bifunctional oxygen electrocatalysts in alkaline media, capable of performing the oxygen reduction reaction (ORR) and oxygen evolution reaction (OER) with lower combined overpotentials relative to their precious metal counterparts. However, their semiconducting nature necessitates the use of activated carbons as conductive supports to generate applicably relevant current densities. In efforts to advance the performance and theory of oxide electrocatalysts, the chemical and physical properties of the oxide material often take precedence over contributions from the conductive additive. In this work, we find that carbon plays an important synergistic role in improving the performance of La1-xSrxCoO3- (0 x 1) electrocatalysts through the activation of O-2 and spillover of radical oxygen intermediates, HO2- and O-2(-), which is further reduced through chemical decomposition of HO2- on the perovskite surface. Through a combination of thin-film rotating disk electrochemical characterization of the hydrogen peroxide intermediate reactions (hydrogen peroxide reduction reaction (HPRR), hydrogen peroxide oxidation reaction (HPOR)) and oxygen reduction reaction (ORR), surface chemical analysis, HR-TEM, and microkinetic modeling on La1-xSrxCoO3- (0 x 1)/carbon (with nitrogen and non-nitrogen doped carbons) composite electrocatalysts, we deconvolute the mechanistic aspects and contributions to reactivity of the oxide and carbon support.

Abstract: A number of recent computational material design studies based on density functional theory (DFT) calculations have put forward a new class of materials with electrically switchable chemical characteristics that can be exploited in the development of tunable gas storage and electrocatalytic applications. We find systematic flaws in almost every computational study of gas adsorption on polarized or charged surfaces, stemming from an improper and unreproducible treatment of periodicity, leading to very large errors of up to 3 eV in some cases. Two simple corrective procedures that lead to consistent results are proposed, constituting a crucial course correction to the research in the field.

Abstract: The interaction between cold atmospheric pressure plasma and liquids is receiving increasing attention for various applications. In particular, the use of plasma-treated liquids (PTL) for biomedical applications is of growing importance, in particular for sterilization and cancer treatment. However, insight into the
underlying mechanisms of plasma–liquid interactions is still scarce. Here, we present a 2D fluid dynamics model for the interaction between a plasma jet and liquid water. Our results indicate that the formed reactive species originate from either the gas phase (with further solvation) or are formed at the liquid interface. A clear increase in the aqueous density of H2O2, HNO2/NO2- and NO3-
is observed as a function of time, while the densities of O3, HO2/O2- and ONOOH/ONOO- are found to quickly reach a maximum due to chemical reactions in solution. The trends observed in our model correlate well with experimental observations from the literature.

Abstract: In this study, we investigated the diffusion of H-atoms to the subsurface and their further diffusion into the bulk of a Ni(111) crystal by means of density functional theory calculations in the context of thermal and plasma-assisted catalysis. The H-atoms at the surface can originate from the dissociative adsorption of H2 or CH4 molecules, determining the surface H-coverage. When a threshold H-coverage is passed, corresponding to 1.00 ML for the crystalline Ni(111) surface, the surface-bound H-atoms start to diffuse to the subsurface. A similar threshold coverage is observed for the interstitial H-coverage. Once the interstitial sites are filled up with a coverage above 1.00 ML of H, dissolution of interstitial H-atoms to the layer below the interstitial sites will be initiated. Hence, by applying a high pressure or inducing a reactive plasma and high temperature, increasing the H-flux to the surface, a large amount of hydrogen can diffuse in a crystalline metal like Ni and can be absorbed. The formation of metal hydride may modify the entire reaction kinetics of the system. Equivalently, the H-atoms in the bulk can easily go back to the surface and release a large amount of heat. In a plasma process, H-atoms are formed in the plasma, and therefore the energy barrier for dissociative adsorption is dismissed, thus allowing achievement of the threshold coverage without applying a high pressure as in a thermal process. As a result, depending on the crystal plane and type of metal, a large number of H-atoms can be dissolved (absorbed) in the metal catalyst, explaining the high efficiency of plasma-assisted catalytic reactions. Here, the mechanism of H-dissolution is established as a chemical identifier for the investigation of the reaction kinetics of a chemical process.