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: Plasma-treated liquids have great potential for biomedical applications. However, insight into the underlying mechanisms and the exact chemistry is still scarce. In this study, we present the combination of a 0D chemical kinetics and a 2D fluid dynamics model to investigate the plasma treatment of a buffered water solution with the kINPen (R) plasma jet. Using this model, we calculated the gas and liquid flow profiles and the transport and chemistry of all species in the gas and the liquid phase. Moreover, we evaluated the stability of the reactive oxygen and nitrogen species after plasma treatment. We found that of all species, only H2O2, HNO2/NO2-, and HNO3/NO3- are stable in the buffered solution after plasma treatment. This is because both their production and loss processes in the liquid phase are dependent on short-lived radicals (e.g. OH, NO, and NO2). Apart from some discrepancy in the absolute values of the concentrations, which can be explained by the model, all general trends and observations in our model are in qualitative agreement with experimental data and literature.

Abstract: Single-chain elasticity of polyethylene at theta point up to 90% of stretching with respect to its contour length is computed by Monte Carlo simulation of an atomistic model in continuous space. The elasticity law together with the free-energy and the internal energy variations with stretching are found to be very well represented by the wormlike chain model up to 65% of the chain elongation, provided the persistence length is treated as a temperature-dependent parameter. Beyond this value of elongation simple ideal chain models are not able to describe the Monte Carlo data in a thermodynamic consistent way. This study reinforces the use of the wormlike chain model to interpret experimental data on the elasticity of synthetic polymers in the finite extensibility regime, provided the chain is not yet in its fully stretched regime. Specific solvent effects on the elasticity law and the partition between energetic and entropic contributions to single chain elasticity are investigated. (C) 2002 American Institute of Physics.

Abstract: Gold nanocluster properties exhibit unique size-dependence. In this contribution, we employ reactive molecular dynamics simulations to calculate the size- and temperature-dependent surface energies, strain energies and atomic displacements for icosahedral, cuboctahedral, truncated octahedral and decahedral Au-nanoclusters. The calculations demonstrate that the surface energy decreases with increasing cluster size at 0 K but increases with size at higher temperatures. The calculated melting curves as a function of cluster size demonstrate the Gibbs-Thomson effect. Atomic displacements and strain are found to strongly depend on the cluster size and both are found to increase with increasing cluster size. These results are of importance for understanding the size-and temperature-dependent surface processes on gold nanoclusters.

Abstract: Two-dimensional materials composed of transition metal carbides and nitrides (MXenes) are poised to revolutionize energy conversion and storage. In this work, we used density functional theory (DFT) to investigate the adsorption of Mg and Na adatoms on five M2CS2 monolayers (where M = Mo, Nb, Ti, V, and Zr) for battery applications. We assessed the stability of the adatom (i.e. Na and Mg)-monolayer systems by calculating adsorption and formation energies, as well as voltages as a function of surface coverage. For instance, we found that Mo2CS2 cannot support a full layer of Na nor even a single Mg atom. Na and Mg exhibit the strongest binding on Zr2CS2, followed by Ti2CS2, Nb2CS2 and V2CS2. Using the nudged elastic band method (NEB), we computed promising diffusion barriers for both dilute and nearly full ion surface coverage cases. In the dilute ion adsorption case, a single Mg and Na atom on Ti2CS2 experience similar to 0.47 eV and similar to 0.10 eV diffusion barriers between the lowest energy sites, respectively. For a nearly full surface coverage, a Na ion moving on Ti2CS2 experiences a similar to 0.33 eV energy barrier, implying a concentration-dependent diffusion barrier. Our molecular dynamics results indicate that the three (one) layers (layer) of the Mg (Na) ion on both surfaces of Ti2CS2 remain stable at T = 300 K. While, according to voltage calculations, Zr2CS2 can store Na up to three atomic layers, our MD simulations predict that the outermost layers detach from the Zr2CS2 monolayer due to the weak interaction between Na ions and the monolayer. This suggests that MD simulations are essential to confirm the stability of an ion-electrode system – an insight that is mostly absent in previous studies.

Abstract: Modeling the oxidation process of silicon nanowires through reactive force field based molecular dynamics simulations suggests that the formation of Si epoxide defects occurs both at the Si/SiOx interface and at the nanowire surface, whereas for flat surfaces, this defect is experimentally observed to occur only at the interface as a result of stress. In this paper, we argue that the increasing curvature stabilizes the defect at the nanowire surface, as suggested by our density functional theory calculations. The latter can have important consequences for the opto-electronic properties of thin silicon nanowires, since the epoxide induces an electronic state within the band gap. Removing the epoxide defect by hydrogenation is expected to be possible but becomes increasingly difficult with a reduction of the diameter of the nanowires.

Abstract: Strain engineering is an effective way to tune the intrinsic properties of a material. Here, we show by using first-principles calculations that both uniaxial and biaxial tensile strain applied to monolayer TiS3 are able to significantly modify its intrinsic mobility. From the elastic modulus and the phonon dispersion relation we determine the tensile strain range where structure dynamical stability of the monolayer is guaranteed. Within this region, we find more than one order of enhancement of the acoustic phonon limited mobility at 300 K (100 K), i.e. from 1.71 x 10(4) (5.13 x 10(4)) cm(2) V-1 s(-1) to 5.53 x 10(6) (1.66 x 10(6)) cm(2) V-1 s(-1). The degree of anisotropy in both mobility and effective mass can be tuned by using tensile strain. Furthermore, we can either increase or decrease the band gap of TiS3 monolayer by applying strain along different crystal directions. This property allows us to use TiS3 not only in electronic but also in optical applications.

Abstract: We report on a first-principles study of the structural and electronic properties of a Sigma3 (112) grain boundary model in CuInSe2. The study focuses on a coherent, stoichiometry preserving, cation–Se terminated grain boundary, addressing the properties of the grain boundary as such, as well as the effect
of well known defects in CuInSe2. We show that in spite of its apparent simplicity, such a grain boundary exhibits a very rich phenomenology, providing an explanation for several of the experimentally observed properties of grain boundaries in CuInSe2 thin films. In particular, we show that the combined effect of Cu vacancies and cation antisites can result in the observed Cu depletion with no In enrichment at the grain boundaries. Furthermore, Cu vacancies are unlikely to produce a hole barrier at the grain boundaries, but Na may indeed have such an effect. We find that Na-on-Cu defects will tend to form abundantly at
the grain boundaries, and can provide a mechanism for the carrier depletion and/or type inversion experimentally reported.

Abstract: Copper deposited by thermal evaporation onto pristine and oxygen plasma treated carbon nanotubes (CNTs) diffuse over the CNT surface, coalescing and forming crystalline islands. The nucleation sites of the islands are preferentially defects, and more homogeneous island dispersion was observed at the CNT oxygen functionalized surface. The presence of weakly bound oxygen atoms at the CNT surface induces the formation of CuO bonds at the Cu/CNT interface, as described through density functional calculations. Exposure to air allows further oxidation to facetted crystalline Cu2O. Oxygen plasma pre-treatment represents a promising route for homogenous disperse Cu2O nanoparticle decoration of CNTs.

Abstract: Nickel foams were used as catalysts to dissociate acetylene and deposit carbon atoms. Graphite nanofibers with distinct structures were developed at 550degreesC with nickel foams pretreated with hydrogen. HREM observations showed that the graphite layers of the nanofibers were aligned at a certain angle to the fiber axis. It is suggested that hydrogen treatment and metal catalysts have a tremendous impact on the yields and microstructures of the graphite nanofibers. The growth mechanism of these fish-bone graphite nanofibers is also discussed.

Abstract: Nitrogen-doped graphene was recently synthesized and was reported to be a catalyst for hydrogen dissociative adsorption under a perpendicular applied electric field (F). In this work, the diffusion of H atoms on N-doped graphene, in the presence and absence of an applied perpendicular electric field, is studied using density functional theory. We demonstrate that the applied field can significantly facilitate the binding of hydrogen molecules on N-doped graphene through dissociative adsorption and diffusion on the surface. By removing the applied field the absorbed H atoms can be released efficiently. Our theoretical calculation indicates that N-doped graphene is a promising hydrogen storage material with reversible hydrogen adsorption/desorption where the applied electric field can act as a switch for the uptake/release processes.

Abstract: ZnM2O4 (M = Co, Rh, Ir) spinels are considered as a class of potential p-type transparent conducting oxides (TCOs). We report the formation energy of acceptor-like defects using first principles calculations with an advanced hybrid exchange-correlation functional (HSE06) within density functional theory (DFT). Due to the discrepancies between the theoretically obtained band gaps with this hybrid functional and the – scattered – experimental results, we also perform GW calculations to support the validity of the description of these spinels with the HSE06 functional. The considered defects are the cation vacancy and antisite defects, which are supposed to be the leading source of disorder in the spinel structures. We also discuss the band alignments in these spinels. The calculated formation energies indicate that the antisite defects ZnM (Zn replacing M, M = Co, Rh, Ir) and VZn act as shallow acceptors in ZnCo2O4, ZnRh2O4 and ZnIr2O4, which explains the experimentally observed p-type conductivity in those systems. Moreover, our systematic study indicates that the ZnIr antisite defect has the lowest formation energy in the group and it corroborates the highest p-type conductivity reported for ZnIr2O4 among the group of ZnM2O4 spinels. To gain further insight into factors affecting the p-type conductivity, we have also investigated the formation of localized small polarons by calculating the self-trapping energy of the holes.