The water/ceria(111) interface: Computational overview and new structures

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.

Predicting Frequency from the External Chemical Environment: OH Vibrations on Hydrated and Hydroxylated Surfaces

Robust correlation curves are essential to decipher structural information from IR-vibrational spectra. However, for surface-adsorbed water and hydroxides, few such correlations have been presented in the literature. In this paper, OH vibrational frequencies are correlated against 12 structural descriptors representing the quantum mechanical or geometrical environment, focusing on those external to the vibrating molecule. A nonbiased fitting procedure based on Gaussian process regression (GPR) was used alongside simple analytical functional forms. The training data consist of 217 structurally unique OH groups from 38 water/metal oxide interface systems for MgO, CaO and CeO₂, all optimized at the DFT level, and the fully anharmonic and uncoupled OH vibrational signatures were calculated. Among our results, we find the following: (i) The intermolecular R(H···O) hydrogen bond distance is particularly strong, indicating the primary cause of the frequency shift. (ii) Similarly, the electric field along the H-bond vector is also a good descriptor. (iii) Highly detailed machine learning descriptors (ACSF, SOAP) are less intuitive but were found to be more capable descriptors. (iv) Combinations of geometric and QM descriptors give the best predictions, supplying complementary information.

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.

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

Ceria (CeO₂) is a well-known catalytic oxide with many environmental, energy production, and industrial applications, most of them involving water as a reactant, byproduct, solvent, or simple spectator. In this work, we parameterized a Ce/O/H ReaxFF for the study of ceria and ceria/water interfaces. The parameters were fitted to an ab initio training set obtained at the DFT/PBE0 level, including the structures, cohesive energies, and elastic properties of the crystalline phases Ce, CeO₂, and Ce₂O₃; the O-defective structures and energies of vacancy formation on CeO₂ bulk and CeO₂ (111) surface, as well as the absorption and reaction energies of H₂ and H₂O molecules on CeO₂ (111). The new potential reproduced reasonably well all the fitted properties as well as the relative stabilities of the different ceria surfaces, the oxygen vacancies formation, and the energies and structures of associative and dissociative water molecules on them. Molecular dynamics simulations of the liquid water on the CeO₂ (111) and CeO₂ (100) surfaces were carried out to study the coverage and the mechanism of water dissociation. After equilibration, on average, 35% of surface sites of CeO₂ (111) are hydroxylated whereas 15% of them are saturated with molecular water associatively adsorbed. As for the CeO₂ (100) surface, we observed that water preferentially dissociates covering 90% of the available surface sites in excellent agreement with recent experimental findings.

Design strategies for ceria nanomaterials: untangling key mechanistic concepts

The morphologies of ceria nanocrystals play an essential role in determining their redox and catalytic performances in many applications, yet the effects of synthesis variables on the formation of ceria nanoparticles of different morphologies and their related growth mechanisms have not been systematised. The design of these morphologies is underpinned by a range of fundamental parameters, including crystallography, optical mineralogy, the stabilities of exposed crystallographic planes, CeO_2-x stoichiometry, phase equilibria, thermodynamics, defect equilibria, and the crystal growth mechanisms. These features are formalised and the key analytical methods used for analysing defects, particularly the critical oxygen vacancies, are surveyed, with the aim of providing a source of design parameters for the synthesis of nanocrystals, specifically CeO_2-x. However, the most important aspect in the design of CeO_2-x nanocrystals is an understanding of the roles of the main variables used for synthesis. While there is a substantial body of data on CeO_2-x morphologies fabricated using low cerium concentrations ([Ce]) under different experimental conditions, the present work fully maps the effects of the relevant variables on the resultant CeO_2-x morphologies in terms of the commonly used raw materials [Ce] (and [NO₃^-] in Ce(NO₃)₃·6H₂O) as feedstock, [NaOH] as precipitating agent, temperature, and time (as well as the complementary vapour pressure). Through the combination of consideration of the published literature and the generation of key experimental data to fill in the gaps, a complete mechanistic description of the development of the main CeO_2-x morphologies is illustrated. Further, the mechanisms of the conversion of nanochains into the two variants of nanorods, square and hexagonal, have been elucidated through crystallographic reasoning. Other key conclusions for the crystal growth process are the critical roles of (1) the formation of Ce(OH)₄ crystallite nanochains as the precursors of nanorods and (2) the disassembly of the nanorods into Ce(OH)₄ crystallites and NO₃^--assisted reassembly into nanocubes (and nanospheres) as an unrecognised intermediate stage of crystal growth.