Development and Application of a ReaxFF Reactive Force Field for Cerium Oxide/Water Interfaces

The Open DAC 2023 Dataset and Challenges for Sorbent Discovery in Direct Air Capture

Direct air capture (DAC) of CO2 with porous adsorbents such as metal-organic frameworks (MOFs) has the potential to aid large-scale decarbonization. Previous screening of MOFs for DAC relied on empirical force fields and ignored adsorbed H2O and MOF deformation. We performed quantum chemistry calculations overcoming these restrictions for thousands of MOFs. The resulting data enable efficient descriptions using machine learning.

Long-range proton and hydroxide ion transfer dynamics at the water/CeO2 interface in the nanosecond regime: reactive molecular dynamics simulations and kinetic analysis

The structural properties, dynamical behaviors, and ion transport phenomena at the interface between water and cerium oxide are investigated by reactive molecular dynamics (MD) simulations employing neural network potentials (NNPs). The NNPs are trained to reproduce density functional theory (DFT) results, and DFT-based MD (DFT-MD) simulations with enhanced sampling techniques and refinement schemes are employed to efficiently and systematically acquire training data that include diverse hydrogen-bonding configurations caused by proton hopping events. The water interfaces with two low-index surfaces of (111) and (110) are explored with these NNPs, and the structure and long-range proton and hydroxide ion transfer dynamics are examined with unprecedented system sizes and long simulation times. Various types of proton hopping events at the interface are categorized and analyzed in detail. Furthermore, in order to decipher the proton and hydroxide ion transport phenomena along the surface, a counting analysis based on the semi-Markov process is formulated and applied to the MD trajectories to obtain reaction rates by considering the transport as stochastic jump processes. Through this model, the coupling between hopping events, vibrational motions, and hydrogen bond networks at the interface are quantitatively examined, and the high activity and ion transport phenomena at the water/CeO2 interface are unequivocally revealed in the nanosecond regime.

Impact of Atomic Defects on Ceria Surfaces on Chemical Mechanical Polishing of Silica Glass Surfaces

This study investigates the impact of atomic defects, such as oxygen vacancies and Ce3+ ions, on cerium oxide (ceria) surfaces during chemical mechanical polishing (CMP) for silica glass finishing. Using density functional theory (DFT) and reactive molecular dynamics simulations, the interaction of orthosilicic molecules and silica glass with dry and wet ceria surfaces is explored. Defects alter the surface reactivity, leading to the dissociation of orthosilicic acid on oxygen vacancies, forming a strong Si-O-Ce bond. Hydroxylated surfaces exhibit easier oxygen vacancy formation and thermodynamically favored substitution of hydroxyl groups with orthosilicic acid. A new ReaxFF library for silica/ceria interfaces with defects is validated using DFT outcomes. Reactive MD simulations demonstrate that ceria surfaces with 30% Ce3+ ions on (111) planes exhibit higher polishing efficiency, attributed to increased Si-O-Ce bond formation. The simultaneous presence of oxygen vacancies and various acidic and basic sites on ceria surfaces enhances the polishing efficiency, involving acid-base reactions with silica. Defective surfaces show superior efficiency by removing silicate chains, contrasting with nondefective surfaces removing isolated orthosilicate units. This study provides insights into optimizing CMP processes for high-precision glass industry surface finishing.

Disclosing gate-opening/closing events inside a flexible metal-organic framework loaded with CO2 by reactive and essential dynamics

We have combined reactive molecular dynamics simulations with principal component analysis to provide a clearer view of the interactions and motion of the CO2 molecules inside a metal-organic framework and the movements of the MOF components that regulate storage, adsorption, and diffusion of the guest species. The tens-of-nanometer size of the simulated model, the capability of the reactive force field tuned to reproduce the inorganic-organic material confidently, and the unconventional use of essential dynamics have effectively disclosed the gate-opening/closing phenomenon, possible coordinations of CO2 at the metal centers, all the diffusion steps inside the MOF channels, the primary motions of the linkers, and the effects of their concerted rearrangements on local CO2 relocations.

New Atomistic Insights on the Chemical Mechanical Polishing of Silica Glass with Ceria Nanoparticles

Reactive molecular dynamics simulations have been used to simulate the chemical mechanical polishing (CMP) process of silica glass surfaces with the ceria (111) and (100) surfaces, which are predominantly found in ceria nanoparticles. Since it is known that an alteration layer is formed at the glass surface as a consequence of the chemical interactions with the slurry solutions used for polishing, we have created several glass surface models with different degrees of hydroxylation and porosity for investigating their morphology and chemistry after the interaction with acidic, neutral, and basic water solutions and the ceria surfaces. Both the chemical and mechanical effects under different pressure and temperature conditions have been studied and clarified. According to the simulation results, we have found that the silica slab with a higher degree of hydroxylation (thicker alteration layer) is more reactive, suggesting that proper chemical treatment is fundamental to augment the polishing efficiency. The reactivity between the silica and ceria (111) surfaces is higher at neutral pH since more OH groups present at the two surfaces increased the Si-O-Ce bonds formed at the interface. Usually, an outermost tetrahedral silicate unit connected to the rest of the silicate network through a single bond was removed during the polishing simulations. We observed that higher pressure and temperature accelerated the removal of more SiO₄ units. However, excessively high pressure was found to be detrimental since the heterogeneous detachment of SiO₄ units led to rougher surfaces and breakage of the Si-O-Si bond, even in the bulk of the glass. Despite the lower concentration of Ce ions at the surface resulting in the lower amount of Si-O-Ce formed, the (100) ceria surface was intrinsically more reactive than (111). The different atomic-scale mechanisms of silica removal at the two ceria surfaces were described and discussed.

Hydrophobic Nanoconfinement Enhances CO2 Conversion to H2CO3

Understanding the formation of H₂CO₃ in water from CO₂ is important in environmental and industrial processes. Although numerous investigations have studied this reaction, the conversion of CO₂ to H₂CO₃ in nanopores, and how it differs from that in bulk water, has not been understood. We use ReaxFF metadynamics molecular simulations to demonstrate striking differences in the free energy of CO₂ conversion to H₂CO₃ in bulk and nanoconfined aqueous environments. We find that nanoconfinement not only reduces the energy barrier but also reverses the reaction from endothermic in bulk water to exothermic in nanoconfined water. Also, charged intermediates are observed more often under nanoconfinement than in bulk water. Stronger solvation and more favorable proton transfer with increasing nanoconfinement enhance the thermodynamics and kinetics of the reaction. Our results provide a detailed mechanistic understanding of an important step in the carbonation process, which depends intricately on confinement, surface chemistry, and CO₂ concentration.

Toward a Consistent Prediction of Defect Chemistry in CeO2

Polarizable shell-model potentials are widely used for atomic-scale modeling of charged defects in solids using the Mott-Littleton approach and hybrid Quantum Mechanical/Molecular Mechanical (QM/MM) embedded-cluster techniques. However, at the pure MM level of theory, the calculated defect energetics may not satisfy the requirement of quantitative predictions and are limited to only certain charged states. Here, we proposed a novel interatomic potential development scheme that unifies the predictions of all relevant charged defects in CeO₂ based on the Mott-Littleton approach and QM/MM electronic-structure calculations. The predicted formation energies of oxygen vacancies accompanied by different excess electron localization patterns at the MM level of theory reach the accuracy of density functional theory (DFT) calculations using hybrid functionals. The new potential also accurately reproduces a wide range of physical properties of CeO₂, showing excellent agreement with experimental and other computational studies. These findings provide opportunities for accurate large-scale modeling of the partial reduction and nonstoichiometry in CeO₂, as well as a prototype for developing robust interatomic potentials for other defective crystals.

H2O2adsorption and dissociation on various CeO2(111) surface models: a first-principles study

Periodic density functional theory (DFT) calculations using the hybrid PBE0 functional and atom-centered Gaussian functions as basis sets were carried out to investigate the absorption and the first steps involved in the decomposition of hydrogen peroxide (H₂O₂) on three different models of the ceria (111) surface. One of the models is a clean surface, and the others are defective and partially hydroxylated ceria surfaces. On the clean surface, we found that the minimum energy path of hydrogen peroxide decomposition involves a three-step process, i.e., adsorption, deprotonation, and formation of the peroxide anion, stabilized through its interaction with the surface at a Ce (IV) site, with activation barriers of less than about 0.5 eV. The subsequent formation of superoxide anions and molecular oxygen species is attributed to electron transfer from the reactants to the Ce (IV) ions underneath. On the defective surface, H₂O₂dissociation is an energetically downhill reaction thermodynamically driven by the healing of the O vacancies, after the reduction and decomposition of H₂O₂into oxygen and water. On the hydroxylated surface, H₂O₂is first adsorbed by forming a favorable H-bond and then undergoes heterolytic dissociation, forming two hydroxyl groups at two vicinal Ce sites.