Structure, dynamics and stability of water/scCO2/mineral interfaces from ab initio molecular dynamics simulations.
Ontology highlight
ABSTRACT:
SUBMITTER: Lee MS
PROVIDER: S-EPMC4600984 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Tin-halide perovskites (THPs) have emerged as promising lead-free perovskites for photovoltaics and photocatalysis applications but still fall short in terms of stability and efficiency with respect to their lead-based counterpart. A detailed understanding of the degradation mechanism of THPs in a water environment is missing. This Letter presents ab initio molecular dynamics (AIMD) simulations to unravel atomistic details of THP/water interfaces comparing methylammonium tin iodide, MASnI3, with the lead-based MAPbI3. Our results reveal facile solvation of surface tin-iodine bonds in MASnI3, while MAPbI3 remains more robust to degradation despite a larger amount of adsorbed water molecules. Additional AIMD simulations on dimethylammonium tin bromide, DMASnBr3, investigate the origins of their unprecedented water stability. Our results indicate the presence of amorphous surface layers of hydrated zero-dimensional SnBr3 complexes which may protect the inner structure from degradation and explain their success as photocatalysts. We believe that the atomistic details of the mechanisms affecting THP (in-)stability may inspire new strategies to stabilize THPs.
| S-EPMC8935372 | biostudies-literature
Project description:Unraveling the atomistic structures of electric double layers (EDL) at electrified interfaces is of paramount importance for understanding the mechanisms of electrocatalytic reactions and rationally designing electrode materials with better performance. Despite numerous efforts dedicated in the past, a molecular level understanding of the EDL is still lacking. Combining the state-of-the-art ab initio molecular dynamics (AIMD) and recently developed computational standard hydrogen electrode (cSHE) method, it is possible to realistically simulate the EDL under well-defined electrochemical conditions. In this work, we report extensive AIMD calculation of the electrified Pt(111)-Had/water interfaces at the saturation coverage of adsorbed hydrogen (Had) corresponding to the typical hydrogen evolution reaction conditions. We calculate the electrode potentials of a series of EDL models with various surface charge densities using the cSHE method and further obtain the Helmholtz capacitance that agrees with experiment. Furthermore, the AIMD simulations allow for detailed structural analyses of the electrified interfaces, such as the distribution of adsorbate Had and the structures of interface water and counterions, which can in turn explain the computed dielectric property of interface water. Our calculation provides valuable molecular insight into the electrified interfaces and a solid basis for understanding a variety of electrochemical processes occurring inside the EDL.
| S-EPMC8395682 | biostudies-literature
Project description:The stability of magnetic states is essential for potential spintronic applications. Here we report on the thermal stability of magnetic states of monovacancy graphene using ab initio molecular dynamics simulations. At room temperature, thermal fluctuations of the graphene lattice induce a rapid magnetic switching between two states with a high and low magnetic moment, indicating that due to the instability of the atomic structure of the vacancy, the associated magnetic moment is thermodynamically unstable. Lowering the temperature can significantly reduce the rate of the switching process and enhance the resident time on the high magnetic state. It stabilizes in the high magnetic state at as low as 30?K. Analyzing the atomic trajectories and the instant electronic structures confirms that these two magnetic states in MD simulations correspond to the magnetic and nonmagnetic states reported in the literatures. Such fluctuations of local magnetic moments are associated with the vertical displacement of the carbon atoms with the unsaturated dangling bond. This study reveals the dynamical correlation between atomic movement and the magnetic switching, and a comprehensive picture of vacancy magnetism in graphene. It has implications in graphene based spintronic devices.
| S-EPMC6345904 | biostudies-other
Project description:We have used ab initio molecular dynamics (AIMD) to characterize water properties using two meta-generalized gradient approximation (meta-GGA) functionals, M06-L-D3 and B97M-rV, and compared their performance against a standard GGA corrected for dispersion, revPBE-D3, at ambient conditions (298 K, and 1 g cm-3 or 1 atm). Simulations of the equilibrium density, radial distribution functions, self-diffusivity, the infrared spectrum, liquid dipole moments, and characterizations of the hydrogen bond network show that all three functionals have overcome the problem of the early AIMD simulations that erroneously found ambient water to be highly structured, but they differ substantially among themselves in agreement with experiment on this range of water properties. We show directly using water cluster data up through the pentamer that revPBE-D3 benefits from a cancellation of its intrinsic functional error by running classical trajectories, whereas the meta-GGA functionals are demonstrably more accurate and would require the simulation of nuclear quantum effects to realize better agreement with all cluster and condensed phase properties.
| S-EPMC6092720 | biostudies-literature
Project description:Simulations of Density Functional Theory-based ab initio molecular dynamics (AIMD) have been performed for a series of aqueous lithium bis(trifluoromethylsulfonyl)imide (LiTFSI) solutions with concentrations ranging from salt-in-water to water-in-salt systems. Analysis of the structure of electrolytes has revealed a preference of Li+ cations to interact with water molecules. In concentrated LiTFSI solutions, water molecules form small associates. The total number of hydrogen bonds (HBs) in the system decreases with salt concentrations, with bonds to water acceptors being only partially replaced by interactions with TFSI anions. Infrared (IR) spectra in the region of the O-H stretching frequency calculated from AIMD trajectories are in good agreement with experimental data. Statistics of oscillations of individual O-H bonds have shown correlations between vibrational frequencies and the structure of HBs formed by water. The changes in the IR spectrum have been related to the varying contributions of different local environments of the water molecules. The abundances of the three spectral components calculated from the simulations agree well with the decomposition of the experimental IR spectra reported in the literature.
| S-EPMC10839825 | biostudies-literature
Project description:We investigated excess electron solvation dynamics in (NH3)n- ammonia clusters in the n = 8-32 size range by performing finite temperature molecular dynamics simulations. In particular, we focused on three possible scenarios. The first case is designed to model electron attachment to small neutral ammonia clusters (n ≤ ∼10) that form hydrogen-bonded chains. The excess electron is bound to the clusters via dipole bound states, and persists with a VDE of ∼160 meV at 100 K for the n = 8 cluster. The coupled nuclear and electronic relaxation is fast (within ∼100 fs) and takes place predominantly by intermolecular librational motions and the intramolecular umbrella mode. The second scenario illustrates the mechanism of excess electron attachment to cold compact neutral clusters of medium size (8 ≤ n ≤ 32). The neutral clusters show increasing tendency with size to bind the excess electron on the surface of the clusters in weakly bound, diffuse, and highly delocalized states. Anionic relaxation trajectories launched from these initial states provide no indication for excess electron stabilization for sizes n < 24. Excess electrons are likely to autodetach from these clusters. The two largest investigated clusters (n = 24 and 32) can accommodate the excess electron in electronic states that are mainly localized on the surface, but they may be partly embedded in the cluster. In the last 500 fs of the simulated trajectories, the VDE of the solvated electron fluctuates around ∼200 meV for n = 24 and ∼500 meV for n = 32, consistent with the values extrapolated from the experimentally observed linear VDE-n-1/3 trend. In the third case, we prepared neutral ammonia cluster configurations, including an n = 48 cluster, that contain possible electron localization sites within the interior of the cluster. Excess electrons added to these clusters localize in cavities with high VDEs up to 1.9 eV for n = 48. The computed VDE values for larger clusters are considerably higher than the experimentally observed photoelectric threshold energy for the solvated electron in bulk ammonia (∼1.4 eV). Molecular dynamics simulations launched from these initial cavity states strongly indicate, however, that these cavity structures exist only for ∼200 fs. During the relaxation the electron leaves the cavity and becomes delocalized, while the cluster loses its initial compactness.
| S-EPMC7458421 | biostudies-literature
Project description:An accurate and efficient ab initio molecular dynamics (AIMD) simulation of liquid water was made possible using the fragment-based approach (J. F. Liu, X. He and J. Z. H. Zhang, Phys. Chem. Chem. Phys., 2017, 19, 11931-11936). In this study, we advance the AIMD simulations using the fragment-based coupled cluster (CC) theory, more accurately revealing the structural and dynamical properties of liquid water under ambient conditions. The results show that the double-donor hydrogen-bond configurations in liquid water are nearly in balance with the single-donor configurations, with a slight bias towards the former. Our observation is in contrast to the traditional tetrahedral water structure. The hydrogen-bond switching dynamics in liquid water are very fast, with a hydrogen-bond life time of around 0.78 picoseconds, determined using AIMD simulation at the CCD/aug-cc-pVDZ level. This time scale is remarkably shorter than the ∼3.0 picoseconds that is commonly obtained from traditional nonpolarized force fields and density functional theory (DFT) based first-principles simulations. Additionally, the obtained radial distribution functions, triplet oxygen angular distribution, diffusion coefficient, and the dipole moment of the water molecule are uniformly in good agreement with the experimental observations. The current high-level AIMD simulation sheds light on the understanding of the structural and dynamical properties of liquid water.
| S-EPMC5885775 | biostudies-literature
Project description:Today's fertilizers rely heavily on mining phosphorus (P) rocks. These rocks are known to become exhausted in near future, and therefore effective P use is crucial to avoid food shortage. A substantial amount of P from fertilizers gets adsorbed onto soil minerals to become unavailable to plants. Understanding P interaction with these minerals would help efforts that improve P efficiency. To this end, we performed a molecular level analysis of the interaction of common organic P compounds (glycerolphosphate (GP) and inositol hexaphosphate (IHP)) with the abundant soil mineral (goethite) in presence of water. Molecular dynamics simulations are performed for goethite-IHP/GP-water complexes using the multiscale quantum mechanics/molecular mechanics method. Results show that GP forms monodentate (M) and bidentate mononuclear (B) motifs with B being more stable than M. IHP interacts through multiple phosphate groups with the 3M motif being most stable. The order of goethite-IHP/GP interaction energies is GP M < GP B < IHP M < IHP 3M. Water is important in these interactions as multiple proton transfers occur and hydrogen bonds are formed between goethite-IHP/GP complexes and water. We also present theoretically calculated infrared spectra which match reasonably well with frequencies reported in literature.
| S-EPMC7795625 | biostudies-literature
Project description:Controlling friction and wear at silica-diamond interfaces is crucial for their relevant applications in tribology such as micro-electromechanical systems and atomic force microscopes. However, the tribological performance on diamond surfaces is highly affected by the working environment where atmospheric gases are present. In this work, we investigate the effects of adsorbed oxygen on the friction and wear of diamond surfaces sliding against silica by massive ab initio molecular dynamics simulations. Different surface orientations, O-coverages, and tribological conditions are considered. The results suggest that diamond surfaces with full oxygen passivation are very effective in preventing surface adhesion, and as a result present extremely low friction and wear. At low oxygen coverage, Si-O-C bond formation was observed as well as atomistic wear initiated from C-C bond breaking at extreme pressure. The analysis of electronic structures of the configurations resulting from key tribochemical reactions clarifies the mechanisms of friction reduction and atomistic wear. Overall, our accurate in silico experiments shed light on the influence of adsorbed oxygen on the tribological properties and wear mechanisms of diamond against silica.
| S-EPMC10520978 | biostudies-literature
Project description:Understanding the adsorption of isolated metal cations from water on to mineral surfaces is critical for toxic waste retention and cleanup in the environment. Heterogeneous nucleation of metal oxyhydroxides and other minerals on material surfaces is key to crystal growth and dissolution. The link connecting these two areas, namely cation dimerization and polymerization, is far less understood. In this work we apply ab initio molecular dynamics calculations to examine the coordination structure of hydroxide-bridged Cu(II) dimers, and the free energy changes associated with Cu(II) dimerization on silica surfaces. The dimer dissociation pathway involves sequential breaking of two Cu2+-OH− bonds, yielding three local minima in the free energy profiles associated with 0-2 OH− bridges between the metal cations, and requires the design of a (to our knowledge) novel reaction coordinate for the simulations. Cu(II) adsorbed on silica surfaces are found to exhibit stronger tendency towards dimerization than when residing in water. Cluster-plus-implicit-solvent methods yield incorrect trends if OH− hydration is not correctly depicted. The predicted free energy landscapes are consistent with fast equilibrium times (seconds) among adsorbed structures, and favor Cu2+ dimer formation on silica surfaces over monomer adsorption. Heterogeneous nucleation of metal oxyhydroxides on mineral surfaces is key to environmentally relevant crystal growth and dissolution, but the cation dimerization steps remain largely unexplored. Here, the authors use ab initio molecular dynamics calculations to examine the coordination structure of hydroxide-bridged Cu(II) dimers, and the free energy changes associated with Cu(II) dimerization on silica surfaces.
| S-EPMC9814296 | biostudies-literature