Ab Initio Molecular Dynamics Simulations Research Articles

Intermolecular Correlations in Liquid Lactic Acid Based on ab initio Molecular Dynamics Simulations Combined with High-Energy X-ray Diffraction Measurements

Intermolecular Correlations in Liquid Lactic Acid Based on ab initio Molecular Dynamics Simulations Combined with High-Energy X-ray Diffraction Measurements

Exploring Li-Ion Transport Properties of Li3TiCl6: A Machine Learning Molecular Dynamics Study

Abstract We performed large-scale molecular dynamics simulations based on a machine-learning force field (MLFF) to investigate the Li-ion transport mechanism in cation-disordered Li3TiCl6 cathode at six different temperatures, ranging from 25°C to 100°C. In this work, the deep neural network method and data generated by ab-initio molecular dynamics (AIMD) simulations were deployed to build a high-fidelity MLFF. Radial distribution functions, Li-ion mean square displacements (MSD), diffusion coefficients, ionic conductivity, activation energy, and crystallographic direction-dependent migration barriers were calculated and compared with corresponding AIMD and experimental data to benchmark the accuracy of the MLFF. Our MSD analysis captured both the self and distinct parts of Li-ion dynamics. The latter reveals that the Li-ions are involved in anti-correlation motion that was rarely reported for solid-state materials. The obtained trajectory from molecular dynamics infers that the Li-ion transportation is mainly through interstitial hopping, which was confirmed by intra- and inter-layer Li-ion displacement with respect to simulation time. Ionic conductivity (1.06 mS/cm) and activation energy (0.29eV) calculated by our simulation are highly comparable with that of experimental values. Overall, the MLFF-based simulation is a promising method to explain the intricate electrochemical properties of the Li3TiCl6 cathode with remarkably reduced computational time

Highly Efficient Acidic Electrosynthesis of Hydrogen Peroxide at Industrial-Level Current Densities Promoted by Alkali Metal Cations.

Acidic H2O2 synthesis through electrocatalytic 2e- oxygen reduction presents a sustainable alternative to the energy-intensive anthraquinone oxidation technology. Nevertheless, acidic H2O2 electrosynthesis suffers from low H2O2 Faradaic efficiencies primarily due to the competing reactions of 4e- oxygen reduction to H2O and hydrogen evolution in environments with high H+ concentrations. Here, we demonstrate the significant effect of alkali metal cations, acting as competing ions with H+, in promoting acidic H2O2 electrosynthesis at industrial-level currents, resulting in an effective current densities of 50‒421 mA cm‒2 with 84‒100% Faradaic efficiency and a production rate of 856‒7842 μmol cm-2 h-1 that far exceeds the performance observed in pure acidic electrolytes or low-current electrolysis. Finite-element simulations indicate that high interfacial pH near the electrode surface formed at high currents is crucial for activating the promotional effect of K+. In situ attenuated total reflection Fourier transform infrared spectroscopy and ab initio molecular dynamics simulations reveal the central role of alkali metal cations in stabilizing the key *OOH intermediate to suppress 4e- oxygen reduction through interacting with coordinated H2O.

Highly Efficient Acidic Electrosynthesis of Hydrogen Peroxide at Industrial‐Level Current Densities Promoted by Alkali Metal Cations

Acidic H2O2 synthesis through electrocatalytic 2e– oxygen reduction presents a sustainable alternative to the energy‐intensive anthraquinone oxidation technology. Nevertheless, acidic H2O2 electrosynthesis suffers from low H2O2 Faradaic efficiencies primarily due to the competing reactions of 4e– oxygen reduction to H2O and hydrogen evolution in environments with high H+ concentrations. Here, we demonstrate the significant effect of alkali metal cations, acting as competing ions with H+, in promoting acidic H2O2 electrosynthesis at industrial‐level currents, resulting in an effective current densities of 50‒421 mA cm‒2 with 84‒100% Faradaic efficiency and a production rate of 856‒7842 μmol cm‐2 h‐1 that far exceeds the performance observed in pure acidic electrolytes or low‐current electrolysis. Finite‐element simulations indicate that high interfacial pH near the electrode surface formed at high currents is crucial for activating the promotional effect of K+. In situ attenuated total reflection Fourier transform infrared spectroscopy and ab initio molecular dynamics simulations reveal the central role of alkali metal cations in stabilizing the key *OOH intermediate to suppress 4e– oxygen reduction through interacting with coordinated H2O.

Stable, while Still Active? A DFT Study of Cu, Ag, and Au Single Atoms at the C3N4/TiO2 Interface.

Hybrid DFT calculations are employed to compare the adsorption and stabilization of Cu, Ag, and Au atoms on graphitic C3N4 and on the heterojunction formed by g- C3N4 and TiO2. While Cu and Ag can be strongly chemisorbed in form of cations on g- C3N4, Au is only weakly physisorbed. On g- C3N4/TiO2, all coinage metal adatoms can be strongly chemisorbed, but, while Cu and Ag forms cations, Au form an Au- species. Ab Initio Molecular Dynamics simulations confirm that the metal adatoms on g-C3N4 are highly mobile at room temperature, while they remain confined in the interfacial spacing between C3N4 and TiO2 on the heterojunction, being both stably bound and reachable for the reactants in a catalytic cycle. Doping g-C3N4/TiO2 with metal single atoms permits thus to generate catalytic systems with tunable charge and chemical properties and improved stability with respect to bare C3N4. Moreover, the changes in the electronic structure of g-C3N4/TiO2 induced by the presence of the metal single atoms are beneficial also for photocatalytic applications.

Anion-Dependent Strength Scale of Interactions in Ionic Liquids from X-ray Photoelectron Spectroscopy, Ab Initio Molecular Dynamics, and Density Functional Theory.

Using a combination of experiments and calculations, we have gained new insights into the nature of anion-cation interactions in ionic liquids (ILs). An X-ray photoelectron spectroscopy (XPS)-derived anion-dependent electrostatic interaction strength scale, determined using XPS core-level binding energies for IL cations, is presented here for 39 different anions, with at least 18 new anions included. Linear correlations of experimental XPS core-level binding energies for IL cations with (a) calculated core binding energies (ab initio molecular dynamics (AIMD) simulations were used to generate high-quality model IL structures followed by single-point density functional theory (DFT) to obtain calculated core binding energies), (b) experimental XPS core-level binding energies for IL anions, and (c) other anion-dependent interaction strength scales led to three main conclusions. First, the effect of different anions on the cation can be related to ground-state interactions. Second, the variations of anion-dependent interactions with the identity of the anion are best rationalized in terms of electrostatic interactions and not occupied valence state/unoccupied valence state interactions or polarizability-driven interactions. Therefore, the XPS-derived anion-dependent interaction strength scale can be explained using a simple electrostatic model based on electrostatic site potentials. Third, anion-probe interactions, irrespective of the identity of the probe, are primarily electrostatic, meaning that our electrostatic interaction strength scale captures some inherent, intrinsic property of anions independent of the probe used to measure the interaction strength scale.

Li-decorated C48B12 and Li12C48B12-impregnated MOF-5 for hydrogen storage: A multi-scale simulation study

In this study, we employed density functional theory (DFT) to investigate the adsorption sites and quantities of Li atoms on the C48B12 surface. Additionally, we examined the hydrogen adsorption properties of the formed Li12C48B12 through ab initio molecular dynamics (AIMD) simulations. Our AIMD simulations confirmed the stability of adsorbed Li atoms on the C48B12 surface, computed adsorbed H2 quantities, and predicted gravimetric densities of H2 adsorption on Li12C48B12 at 77 K and 298 K. Grand Canonical Monte Carlo simulations explored the H2 gravimetric densities of pure MOF-5, C48B12-impregnated MOF-5, and Li12C48B12-impregnated MOF-5 at 77 K and 298 K within 100 bar. Li doping on the C48B12 surface did not significantly enhance the gravimetric densities of H2 in MOF-5 but enhanced the stability of H2 adsorption in MOF-5.

Revisiting thermodynamics in (LiF, NaF, KF, CrF2)–CrF3 by first-principles calculations and CALPHAD modeling

The thermodynamic description of the (LiF, NaF, KF, CrF2)–CrF3 systems has been revisited, aiming for a better understanding of the effects of Cr on the FLiNaK molten salt. First-principles calculations based on density functional theory (DFT) were performed to determine the electronic and structural properties of each compound, including the formation enthalpy, volume, and bulk modulus. DFT-based phonon calculations were carried out to determine the thermodynamic properties of compounds, for example, enthalpy, entropy, and heat capacity as functions of temperature. Phonon-based thermodynamic properties show a good agreement with experimental data of binary compounds LiF, NaF, KF, CrF3, and CrF2, establishing a solid foundation to determine thermodynamic properties of ternary compounds as well as to verify results estimated by the Neumann-Kopp rule. Additionally, DFT-based ab initio molecular dynamics (AIMD) simulations were employed to predict the mixing enthalpies of liquid salts. Using DFT-based results and experimental data in the literature, the (LiF, NaF, KF, CrF2)–CrF3 system has been remodeled in terms of the CALculation of PHAse Diagrams (CALPHAD) approach using the modified quasichemical model with quadruplet approximation (MQMQA) for liquid. Calculated phase stability in the present work shows an excellent agreement with experiments, indicating the effectiveness of combining DFT-based total energy, phonon, and AIMD calculations, and CALPHAD modeling to provide the thermodynamic description in complex molten salt systems.

Neutron scattering and neural-network quantum molecular dynamics investigation of the vibrations of ammonia along the solid-to-liquid transition

Vibrational spectroscopy allows us to understand complex physical and chemical interactions of molecular crystals and liquids such as ammonia, which has recently emerged as a strong hydrogen fuel candidate to support a sustainable society. We report inelastic neutron scattering measurement of vibrational properties of ammonia along the solid-to-liquid phase transition with high enough resolution for direct comparisons to ab-initio simulations. Theoretical analysis reveals the essential role of nuclear quantum effects (NQEs) for correctly describing the intermolecular spectrum as well as high energy intramolecular N-H stretching modes. This is achieved by training neural network models using ab-initio path-integral molecular dynamics (PIMD) simulations, thereby encompassing large spatiotemporal trajectories required to resolve low energy dynamics while retaining NQEs. Our results not only establish the role of NQEs in ammonia but also provide general computational frameworks to study complex molecular systems with NQEs.

The development of thermodynamically consistent and physics-informed equation-of-state model through machine learning

Ab initio molecular dynamics (AIMD) simulations have become an important tool used in the construction of equations of state (EOS) tables for warm dense matter. Due to computational costs, only a limited number of system state conditions can be simulated, and the remaining EOS surface must be interpolated for use in radiation-hydrodynamic simulations of experiments. In this work, we develop a thermodynamically consistent EOS model that utilizes a physics-informed machine learning approach to implicitly learn the underlying Helmholtz free-energy from AIMD generated energies and pressures. The model, referred to as PIML-EOS, was trained and tested on warm dense polystyrene producing a fit within a 1% relative error for both energy and pressure and is shown to satisfy both the Maxwell and Gibbs–Duhem relations. In addition, we provide a path toward obtaining thermodynamic quantities, such as the total entropy and chemical potential (containing both ionic and electronic contributions), which are not available from current AIMD simulations.

In Situ Species Analysis of a Lithium-Ion Battery Electrolyte Containing LiTFSI and Propylene Carbonate.

In this study, we analyzed the species in a model electrolyte consisting of a lithium salt, lithium bis(trifluoromethane sulfone)imide (LiTFSI), and a widely used neutral solvent propylene carbonate (PC) with excess infrared (IR) spectroscopy, ab initio molecular dynamics simulations (AIMD), and quantum chemical calculations. Complexing species including the charged ones [Li+(PC)4, TFSI-, TFSI-(PC), TFSI-(PC)2, and Li(TFSI)2-] are identified in the electrolyte. Quantum chemical calculations show strong Li+···O(PC) interaction, which suggests that Li+ would transport in the mode of solvation-carriage. However, the interaction energy of each hydrogen bond in TFSI-(PC) is very weak, suggesting that TFSI- would transport in hopping mode. In addition, the concentration dependences of the relative population of the species were also derived, providing a scenario for the dissolving process of the salt in PC. These in-depth studies provide physical insights into the structural and interactive properties of the electrolyte of lithium-ion batteries.

DFT and AIMD insights into heterogeneous dissociation of 2-chlorothiophenol on CuO(111) surface: Impact of H2O and OH

Copper oxides are vital catalysts in facilitating the formation of polychlorinated thianthrenes/dibenzothiophenes (PCTA/DTs) through heterogeneous reactions in high-temperature industrial processes. Chlorothiophenols (CTPs) are the most crucial precursors for PCTA/DT formation. The initial step in this process is the metal-catalyzed production of chlorothiophenoxy radicals (CTPRs) from CTPs via dissociation reactions. This work combines density functional theory (DFT) calculations with ab initio molecular dynamics (AIMD) simulations to explore the formation mechanism of the adsorbed 2-CTPR from 2-CTP, with the assistance of CuO(111). Our study demonstrates that flat adsorption configurations of 2-CTP on the CuO(111) surface are more stable than vertical configurations. The CuO(111) surface acts as a strong catalyst, facilitating the dissociation of 2-CTP into the adsorbed 2-CTPR. Surface oxygen vacancies enhance the adsorption of 2-CTP on the CuO(111) surface, while moderately suppressing the dissociation of 2-CTP. More importantly, water molecules and surface hydroxyl groups actively promote the dissociation of 2-CTP. Specifically, water directly participates in the reaction through “water bridge”, enabling a barrier-free process. This research provides molecular-level insights into the heterogeneous generation of dioxins with the catalysis of metal oxides in fly ash from static and dynamic aspects, providing novel approaches for reducing dioxin emissions and establishing dioxin control strategies.

High-Performance Dendrite-Free Lithium Metal Anode Based on Metal-Organic Framework Glass.

The performance of lithium metal batteries is severely hampered by uncontrollable dendrite growth and volume change within the anode. This work tackles these obstacles by introducing a novel strategy: applying an isotropic and internal grain-boundary-free layer, specifically, a metal-organic framework (MOF) glass layer with nano-porosity onto the electrochemical plated lithium metal anode. Both ab initio and classical molecular dynamics simulations indicate that the MOF glass layer makes the lithium transport smooth and uniform via its internal monolithic and interfacial advantages. This MOF glass layer with the fast and more uniform lithium diffusion in the monolithic interior and its interface enables dendrite-free lithium plating and stripping through surface confinement effect and interfacial effect. When employed in symmetric batteries, the achieved Li metal anode can operate over 300 hours at 1mA cm-2. The full batteries matched with LiFePO4 exhibit high capacity (148 mAh g-1), excellent rate performance (61 mAh g-1 at 5 C), and outstanding cycling stability (with capacity retention of ∼90% after 1000 cycles). The full batteries matched with high-voltage LiCoO2 also show superior performances. Therefore, the strategy of utilizing a MOF glass layer enables the development of high-performance lithium metal anodes. This article is protected by copyright. All rights reserved.

Finding Natural, Dense, and Stable Frustrated Lewis Pairs on Wurtzite Crystal Surfaces for Small‐Molecule Activation

AbstractThe surface frustrated Lewis pairs (SFLPs) open up new opportunities for substituting noble metals in the activation and conversion of stable molecules. However, the applications of SFLPs on a larger scale are impeded by the complex construction process, low surface density, and sensitivity to the reaction environment. Herein, wurtzite‐structured crystals such as GaN, ZnO, and AlP are found for developing natural, dense, and stable SFLPs. It is revealed that the SFLPs can naturally exist on the (100) and (110) surfaces of wurtzite‐structured crystals. All the surface cations and anions serve as the Lewis acid and Lewis base in SFLPs, respectively, contributing to the surface density of SFLPs as high as 7.26×1014 cm−2. Ab initio molecular dynamics simulations indicate that the SFLPs can keep stable under high temperatures and the reaction atmospheres of CO and H2O. Moreover, outstanding performance for activating the given small molecules is achieved on these natural SFLPs, which originates from the optimal orbital overlap between SFLPs and small molecules. Overall, these findings not only provide a simple method to obtain dense and stable SFLPs but also unfold the nature of SFLPs toward the facile activation of small molecules.

CO Intermediate‐Assisted Dynamic Cu Sintering During Electrocatalytic CO2 Reduction on Cu−N−C Catalysts

AbstractThe electrochemical CO2 reduction reaction (eCO2RR) to multicarbon products has been widely recognized for Cu‐based catalysts. However, the structural changes in Cu‐based catalysts during the eCO2RR pose challenges to achieving an in‐depth understanding of the structure–activity relationship, thereby limiting catalyst development. Herein, we employ constant‐potential density functional theory calculations to investigate the sintering process of Cu single atoms of Cu−N−C single‐atom catalysts into clusters under eCO2RR conditions. Systematic constant‐potential ab initio molecular dynamics simulations revealed that the leaching of Cu−(CO)x moieties and subsequent agglomeration into clusters can be facilitated by synergistic adsorption of H and eCO2RR intermediates (e.g., CO). Increasing the Cu2+ concentration or the applied potential can efficiently suppress Cu sintering. Both microkinetic simulations and experimental results further confirm that sintered Cu clusters play a crucial role in generating C2 products. These findings provide significant insights into the dynamic evolution of Cu‐based catalysts and the origin of their activity toward C2 products during the eCO2RR.

CAIMD Simulations Guided Design of Atomic Praseodymium Doping In–Bi Nanofibers for High‐Energy‐Efficiency CO2 Electrolysis to Formate in Ultra‐Wide Potential Window

AbstractThe electrochemical CO2 reduction reaction (ECO2RR) has emerged as a promising technology for achieving carbon neutralization. Even though considerable efforts are dedicated to gain deep insight into the understanding of ECO2RR on a mechanism level through density functional theory (DFT) studies, effects of solvent molecules and temperature have long been neglected by conventional DFT calculations as a consequence of limitations in current technologies and computational power of supercomputers. Under this context, the lack of comprehensive understanding over the energy changes in the reaction derived from the only concern on free energy changes between reaction intermediates have arouse an urgent call for exploring feasible calculation options toward generalized theoretical study. Here, a systematic mechanism study is provided toward ECO2RR via constrained ab initio molecular dynamics (cAIMD) simulations, in which the effects of solvent water molecules and temperature are taken into consideration to guide the synthesis of single‐atom alloy (SAA) catalyst. Consequently, the resultant Pr0.05/InBi achieves a maximum Faradic efficiency (FE) of 96.4% and an energy efficacy (EE) of 59.41% for formate. This work offers a novel approach to the design and screening of SAA catalysts, presenting foreseeable future in accelerating the industrial application process of ECO2RR.

An investigation of self-interstitial diffusion in α-zirconium by an on-the-fly machine learning force field

The on-the-fly machine learning force field approach, based on the Gaussian approximation potential and Bayesian error estimation, was used to study the diffusion of self-interstitial atoms in α-zirconium. Ab initio molecular dynamics simulations of lattice vibration and interstitial diffusion at different temperatures were employed to develop the force field. The radial and angular descriptors of the potential were further optimized to achieve better agreement with first-principles results. Subsequent long-term diffusion simulations were performed to assess the diffusion behavior based on the obtained force field. Tracer diffusion coefficients and diffusion anisotropy were studied at temperatures of 600–1200 K, and the Bayesian errors were estimated throughout the diffusion simulations. The mean and maximum estimated Bayesian errors of atomic force were approximately twice as large as those observed during the learning period. The basal diffusion was greatly favored compared to the interstitial diffusion along the c-axis, consistent with previous simulations based on first-principles results and classical potentials. The accuracy and applicability of the current on-the-fly machine learning approach were critically evaluated.

Density functional theory-based screening of Ti4C3O2-loaded single atoms for efficient selective catalytic oxidation of formaldehyde

Indoor formaldehyde (HCHO) pollution poses a major risk to human health. Low-temperature catalytic oxidation is an effective method for HCHO removal. The high activity and selectivity of single atomic catalysts provide a possibility for the development of efficient non-precious metal catalysts. In this study, the most stable single-atom catalyst Ti-Ti4C3O2 was screened by density functional theory among many single atomic catalysts with two-dimensional (2D) monolayer Ti4C3O2 as the support. The computational results show that Ti-Ti4C3O2 is highly selective to HCHO and O2 in complex environments. The HCHO oxidation reaction pathways are proposed based on the Eley-Rideal (E-R) and Langmuir-Hinshelwood (L-H) mechanisms. According to the reaction energy and energy span models, the E-R mechanism has a lower maximum energy barrier and higher catalytic efficiency than the L-H mechanism. In addition, the stability of the Ti-Ti4C3O2 structure and active center was verified by diffusion energy barrier and ab initio molecular dynamics simulations. The above results indicate that Ti-Ti4C3O2 is a promising non-precious metal catalyst. The present study provides detailed theoretical insights into the catalytic oxidation of HCHO by Ti-Ti4C3O2, as well as an idea for the development of efficient non-precious metal catalysts based on 2D materials.

Novel arsenborane as a sensing material for dichlorofluoromethane and dichlorodifluoromethane vapours – a first-principles study

Presently, two-dimensional (2D) materials are cutting-edge research to detect air pollutants in the environment. In this research study, we used arsenborane as a sensing element towards halomethanes viz. dichlorofluoromethane (Freon 21) and dichlorodifluoromethane (Freon 12). At first, the structural firmness of bare arsenborane is confirmed based on formation energy, phonon band spectrum, and ab initio molecular dynamics simulation. Also, the electronic property of bare stable arsenborane is investigated by plotting bands and projected-density-of-states (PDOS) maps. The energy gap of arsenborane is computed to be 0.927eV which confirms semiconductor nature. Significantly, the interaction behaviour of Freon 21 and Freon 12 on arsenborane is studied with the influence of relative band gap changes, Mulliken-charge analysis, and adsorption energy. In the present work, the calculated value of adsorption energy is observed in the scope of -0.083eV to -0.618eV (physisorption regime). Thus, based on the outcome we suggest that the new 2D material arsenborane is a good base material to detect Freon 21 and Freon 12 vapours in the air environment.

Insight into the effects of gas molecules-adsorbed on 2D-FeS2: A DFT study

Lithium-sulfur (Li-S) batteries are especially competitive in the energy sector due to their excellent performances, like preferable energy density and economic benefits. Studying the adsorption of gas molecules on electrode materials has potential engineering significance for Li-S batteries since they have a highly osmotic potential, which causes unavoidable damage to batteries. In this work, the adsorption phenomenon of common gas molecules (H2O, N2, H2, CO2, and O2) on the two-dimensional pyrite (2D-FeS2) cathode material surface, as well as the effects on the electronic and electrochemical properties, were investigated by the first-principles calculations. The adsorption capabilities were estimated by adsorption energy and Mulliken population analysis. Simulation results demonstrated that whole adsorption energies were less than -1.0 eV and larger than -0.6 eV, which shows a physisorption nature. Among them, the O-S bond of O2/2D-FeS2 has the strongest strength. Electronic structure calculations suggested that 2D-FeS2 maintained good conductivity after gas molecules were adsorbed, achieving efficient transfer between electron, lithium, and sulfur intermediates. Additionally, ab initio molecular dynamics (AIMD) simulations showed that Li+ exhibits excellent diffusion performance and low activation energy at different temperatures. 2D-FeS2 still has a stable electrochemical working window (1.87 ∼ 2.47 V), while the theoretical open current voltage is damaged by gas molecule adsorption. Consequently, this work theoretically reveals the effect of gas molecules on the cathode materials for Li-S batteries, which has guide meaning for engineering.