Abstract: The electronic properties of graphene nano-flakes (GNFs) with different edge passivation are investigated by using density functional theory. Passivation with F and H atoms is considered: C-Nc X-Nx (X = F or H). We studied GNFs with 10 < N-c < 56 and limit ourselves to the lowest energy configurations. We found that: (i) the energy difference Delta between the highest occupied molecular orbital and the lowest unoccupied molecular orbital decreases with N-c, (ii) topological defects (pentagon and heptagon) break the symmetry of the GNFs and enhance the electric polarization, (iii) the mutual interaction of bilayer GNFs can be understood by dipole-dipole interaction which were found sensitive to the relative orientation of the GNFs, (iv) the permanent dipoles depend on the edge terminated atom, while the energy gap is independent of it, and (v) the presence of heptagon and pentagon defects in the GNFs results in the largest difference between the energy of the spin-up and spin-down electrons which is larger for the H-passivated GNFs as compared to F-passivated GNFs. Our study shows clearly the effect of geometry, size, termination, and bilayer on the electronic properties of small GNFs. This study reveals important features of graphene nano-flakes which can be detected using Raman spectroscopy. (C) 2014 AIP Publishing LLC.

Abstract: Uniform-acceptance force-bias Monte Carlo (fbMC) methods have been shown to be a powerful technique to access longer timescales in atomistic simulations allowing, for example, phase transitions and growth. Recently, a new fbMC method, the time-stamped force-bias Monte Carlo (tfMC) method, was derived with inclusion of an estimated effective timescale; this timescale, however, does not seem able to explain some of the successes the method. In this contribution, we therefore explicitly quantify the effective timescale tfMC is able to access for a variety of systems, namely a simple single-particle, one-dimensional model system, the Lennard-Jones liquid, an adatom on the Cu(100) surface, a silicon crystal with point defects and a highly defected graphene sheet, in order to gain new insights into the mechanisms by which tfMC operates. It is found that considerable boosts, up to three orders of magnitude compared to molecular dynamics, can be achieved for solid state systems by lowering of the apparent activation barrier of occurring processes, while not requiring any system-specific input or modifications of the method. We furthermore address the pitfalls of using the method as a replacement or complement of molecular dynamics simulations, its ability to explicitly describe correct dynamics and reaction mechanisms, and the association of timescales to MC simulations in general.

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: 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 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: Atomistic simulation methods for the quantification of free energies are in wide use. These methods operate by sampling the probability density of a system along a small set of suitable collective variables (CVs), which is, in turn, expressed in the form of a free energy surface (FES). This definition of the FES can capture the relative stability of metastable states but not that of the transition state because the barrier height is not invariant to the choice of CVs. Free energy barriers therefore cannot be consistently computed from the FES. Here, we present a simple approach to calculate the gauge correction necessary to eliminate this inconsistency. Using our procedure, the standard FES as well as its gauge-corrected counterpart can be obtained by reweighing the same simulated trajectory at little additional cost. We apply the method to a number of systems—a particle solvated in a Lennard-Jones fluid, a Diels–Alder reaction, and crystallization of liquid sodium—to demonstrate its ability to produce consistent free energy barriers that correctly capture the kinetics of chemical or physical transformations, and discuss the additional demands it puts on the chosen CVs. Because the FES can be converged at relatively short (sub-ns) time scales, a free energy-based description of reaction kinetics is a particularly attractive option to study chemical processes at more expensive quantum mechanical levels of theory.

Abstract: We calculate bubble nucleation rates in a Lennard-Jones fluid through explicit molecular dynamics simulations. Our approach-based on a recent free energy method (dubbed reweighted Jarzynski sampling), transition state theory, and a simple recrossing correction-allows us to probe a fairly wide range of rates in several superheated and cavitation regimes in a consistent manner. Rate predictions from this approach bridge disparate independent literature studies on the same model system. As such, we find that rate predictions based on classical nucleation theory, direct brute force molecular dynamics simulations, and seeding are consistent with our approach and one another. Published rates derived from forward flux sampling simulations are, however, found to be outliers. This study serves two purposes: First, we validate the reliability of common modeling techniques and extrapolation approaches on a paradigmatic problem in materials science and chemical physics. Second, we further test our highly generic recipe for rate calculations, and establish its applicability to nucleation processes.

Abstract: Binding of the SARS-CoV-2 S-glycoprotein to cell receptors is vital for the entry of the virus into cells and subsequent infection. ACE2 is the main cell receptor for SARS-CoV-2, which can attach to the C-terminal receptor-binding domain (RBD) of the SARS-CoV-2 S-glycoprotein. The GRP78 receptor plays an anchoring role, which attaches to the RBD and increases the chance of other RBDs binding to ACE2. Although high levels of reactive oxygen and nitrogen species (RONS) are produced during viral infections, it is not clear how they affect the RBD structure and its binding to ACE2 and GRP78. In this research, we apply molecular dynamics simulations to study the effect of oxidation of the highly reactive cysteine (Cys) amino acids of the RBD on its binding to ACE2 and GRP78. The interaction energy of both ACE2 and GRP78 with the whole RBD, as well as with the RBD main regions, is compared in both the native and oxidized RBDs. Our results show that the interaction energy between the oxidized RBD and ACE2 is strengthened by 155 kJ/mol, increasing the binding of the RBD to ACE2 after oxidation. In addition, the interaction energy between the RBD and GRP78 is slightly increased by 8 kJ/mol after oxidation, but this difference is not significant. Overall, these findings highlight the role of RONS in the binding of the SARS-CoV-2 S-glycoprotein to host cell receptors and suggest an alternative mechanism by which RONS could modulate the entrance of viral particles into the cells.

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.