Abstract: We have studied the molecular energy terms of the hole configurations (h(u)(+))(m), m=2,3,4,5 of C-60(m+) cations and the electronic configurations (t(1u))(n) n=2,3,4, as well as (t(1u))(n-1)t(1g) of the C-60(n-) anions. The lowest terms (within an energy span of 0.03 eV) for C-60(2+) are three triplets T-3(1g), (3)G(g), T-3(2g) and for C-60(3+) are three quartets T-4(1u), (4)G(u), T-4(2u), which favor Jahn-Teller distortions. For the ground state of C-60(2-) we find a triplet T-3(1g) in agreement with Hund's rules. Our method takes into account intramolecular direct and exchange multipolar Coulomb interactions.

Abstract: A microscopic theory, which describes the orientational dynamics of C60 molecules in the face-centered cubic phase of C60-fullerite, is formulated or the case of a complex molecular structure. Interaction centers which comprise atoms, double bonds, and single bonds as molecular constituents contribute to the intermolecular potential. Orientation dependent physical properties are described in terms of symmetry-adapted rotator functions. It is found that a same set of rotator functions is sufficient even in the case of a complex molecular structure. Phase transition temperatures are discussed for various models of molecular structure.

Abstract: In present work, we have performed the Grand Canonical Monte Carlo (GCMC) simulations to quantify CO2 capture inside porous silica at high operating temperatures of 673.15 K and 873.15 K; and over a operating pressure range of 500 kPa – 4000 kPa that are methane steam reforming process parameters. Related chemical potential values at these thermodynamic conditions are obtained from the bulk phase simulations in the Canonical ensemble in conjunction with Widom’s insertion technique, where the CO2 has been accurately represented by TraPPE force field. Present structure of the porous silica is a single slit pore geometry of various heights (H = 20 Å, 31.6 Å, 63.2 Å and 126.5 Å), dimensions in which possible vapour-liquid equilibria for generic square well fluids has been reported in literature. Estimation of the pore-fluid interactions show a higher interaction between silica pore and adsorbed CO2 compared to the reported pore-fluid interactions between homogeneous carbon slit pore and adsorbed CO2; thus resulting in an enhancement of adsorption inside silica pores of H = 20 Å and H = 126.5 Å, which are respectively 3.5 times and 1.5 times higher than that in homogeneous carbon slit pores of same dimensions and at 673.15 K and 500 kPa. Estimated local density plots indicate the presence of structured layers due to more molecular packing, which confirms possible liquid-like and vapour-like phase coexistence of the supercritical bulk phase CO2 under confinement.

Abstract: MFI zeolite nanosheets with thickness of 2 and 8 nm were synthesized, transformed into bifunctional catalysts by loading with platinum and tested in n-decane isomerization and hydrocracking. Detailed analysis of skeletal isomers and hydrocracked products revealed that the MFI nanosheets display transition-state shape-selectivity similar to bulk MFI zeolite crystals. The suppressed formation of bulky skeletal isomers and C5 cracking products are observed both in the nanosheets and the bulk crystals grown in three dimensions. This is typical for restricted transition-state shape-selectivity, characteristic for the MFI type pores. It is a first clear example of transition-state shape-selectivity inside a zeolitic nanosheet. Owing to the short diffusion path across the sheets, expression of diffusion-based discrimination of reaction products in the MFI nanosheets was limited. The 2-methylnonane formation among monobranched C10 isomers and 2,7-dimethyloctane among dibranched C10 isomers, which in MFI zeolite are favored by product diffusion, was much less favored on the nanosheets compared to the reference bulk ZSM-5 material.

Abstract: We have studied the dependence of the binding energy of a cubane dimer on the mutual orientation of and the distance between the composing monomers employing the second-order Møller-Plesset perturbation scheme (MP2) with the cc-pVDZ molecular basis set. We have found that the MP2 contribution from the molecular correlations is responsible for the bound state of the cubane dimer, whereas the Hartree-Fock contribution remains anti-bonding at all intermolecular distances. Starting with two molecules in the standard orientation and centers of mass at (0,0,0) and (0,0,d), respectively, the maximal binding energy is found at d = 5.125 Å and one of the monomers rotated by 45° about the z-axis. This configuration implies that the hydrogen atoms belonging to different monomers tend to repel each other. The results are in agreement with experimental data on the optimal packing of cubane molecules in the solid state.

Abstract: Atomic scale simulations have been shown to be a powerful tool for elucidating the growth mechanisms of carbon nanotubes. The growth picture is however not entirely clear yet due to the gap between current simulations and real experiments. We here simulate for the first time the nucleation and subsequent growth of single-wall carbon nanotubes (SWNTs) from oxygen-containing hydrocarbon feedstocks using the hybrid Molecular Dynamics/Monte Carlo technique. The underlying nucleation mechanisms of Ni-catalysed SWNT growth are discussed in detail. Specifically, we find that as a function of the feedstock, different carbon fractions may emerge as the main growth species, due to a competition between the feedstock decomposition, its rehydroxylation and its contribution to etching of the growing SWNT. This study provides a further understanding of the feedstock effects in SWNT growth in comparison with available experimental evidence as well as with<italic>ab initio</italic>and other simulation data, thereby reducing the simulation–experiment gap.

Abstract: Biomembranes are under constant attack of free radicals that may lead to lipid oxidation in conditions of oxidative stress. The products generated during lipid oxidation are responsible for structural and dynamical changes which may jeopardize the membrane function. For instance, the local rearrangements of oxidized lipid molecules may induce membrane rupture. In this study, we investigated the effects of mechanical stress on oxidized phospholipid bilayers (PLBs). Model bilayers were stretched until pore formation (or poration) using nonequilibrium molecular dynamics simulations. We studied single-component homogeneous membranes composed of lipid oxidation products, as well as two-component heterogeneous membranes with coexisting native and oxidized domains. In homogeneous membranes, the oxidation products with —OH and —OOH groups reduced the areal strain required for pore formation, whereas the oxidation product with ]O group behaved similarly to the native membrane. In heterogeneous membranes composed of oxidized and non-oxidized domains, we tested the hypothesis according to which poration may be facilitated at the domain interface region. However, results were inconclusive due to their large statistical variance and sensitivity to simulation setup parameters. We pointed out important technical issues that need to be considered in future simulations of mechanically-induced poration of heterogeneous membranes. This research is of interest for photodynamic therapy and plasma medicine, because ruptured and intact plasma membranes are experimentally considered hallmarks of necrotic and apoptotic cell death.

Abstract: Molecular dynamics simulations using improved TersoffBrenner potential parameters were performed to investigate Cl+ etching of a {2×1} reconstructed Si(100) surface. Steady-state Si etching accompanying the Cl coverage of the surface is observed. Furthermore, a steady-state chlorinated reaction layer is formed. The thickness of this reaction layer is found to increase with increasing energy. The stoichiometry of SiClx species in the reaction layer is found to be SiCl:SiCl2:SiCl3 = 1.0:0.14:0.008 at 50 eV. These results are in excellent agreement with available experimental data. While elemental Si products are created by physical sputtering, most SiClx (0<x<4) etch products are produced by chemical-enhanced physical sputtering.

Abstract: Plasma-surface interactions are in general highly complex due to the interplay of many concurrent processes. Molecular dynamics simulations provide insight in some of these processes, subject to the accessible time and length scales, and the availability of suitable force fields. In this introductory tutorial-style review, we aim to describe the current capabilities and limitations of molecular dynamics simulations in this field, restricting ourselves to low-temperature nonthermal plasmas. Attention is paid to the simulation of the various fundamental processes occurring, including sputtering, etching, implantation, and deposition, as well as to what extent the basic plasma components can be accounted for, including ground state and excited species, electric fields, ions, photons, and electrons. A number of examples is provided, giving an bird’s eye overview of the current state of the field.

Abstract: Molecular dynamics method was employed to investigate the effects of the reaction layer formed near the surface region on CF(3)(+) etching of Si at different temperatures. The simulation results show that the coverages of F and C are sensitive to the surface temperature. With increasing temperature, the physical etching is enhanced, while the chemical etching is weakened. It is found that with increasing surface temperature, the etching rate of Si increases. As to the etching products, the yields of SiF and SiF(2) increase with temperature, whereas the yield of SiF(3) is not sensitive to the surface temperature. And the increase of the etching yield is mainly due to the increased desorption of SiF and SiF(2). The comparison shows that the reactive layer plays an important part in the subsequeat impacting, which enhances the etching rate of Si and weakens the chemical etching intensity.