First-Principles Modeling of Equilibration Dynamics of Hyperthermal Products of Surface Reactions Using Scalable Neural Network Potential

Qidong Lin,Bin Jiang
2023-04-21
Abstract:Equilibration dynamics of hot oxygen atoms following O2 dissociation on Pd(100) and Pd(111) surfaces are investigated by molecular dynamics simulations based on a scalable neural network potential enabling first-principles description of O2 and O interacting with variable Pd supercells. We find that to accurately describe the equilibration dynamics after dissociation, the simulation cell length necessarily exceeds twice the maximum distance of equilibrated oxygen adsorbates. By analyzing hundreds of trajectories with appropriate initial sampling, the measured distance distribution of equilibrated atom pairs on Pd(111) is well reproduced. However, our results on Pd(100) suggest that the ballistic motion of hot atoms predicted previously is a rare event under ideal conditions, while initial molecular orientation and surface thermal fluctuation could significantly affect the overall post-dissociation dynamics. On both surfaces, dissociated oxygen atoms remain primarily locate their nascent positions and then randomly cross bridge sites nearby.
Materials Science,Other Condensed Matter
What problem does this paper attempt to address?
The paper aims to study the equilibrium kinetics of hyperthermal products following the dissociation of oxygen molecules (O2) on palladium (Pd) surfaces. Specifically, the authors utilized scalable neural network potential energy surfaces (PES) to conduct molecular dynamics simulations, describing the interactions between O2 and its dissociation products with Pd superlattices of different sizes from first principles. The main objectives of the paper are: 1. **Validate the simulation unit size**: Determine that the length of the simulation unit must be more than twice the maximum distance of the final stable position of the dissociated oxygen atoms to accurately describe the post-dissociation kinetic process. 2. **Explore the rarity of spherical motion**: Discover that the previously assumed spherical motion of thermal oxygen atoms on the Pd(100) surface is actually a rare event under ideal initial conditions. 3. **Analyze kinetic behavior**: Through extensive trajectory simulations, reveal that dissociated oxygen atoms tend to first locate near their initial positions and subsequently undergo random walk-like diffusion, consistent with experimental results. In summary, the paper addresses the high computational cost and limited predictive accuracy issues present in traditional quantum mechanics/metal embedding (QM/Me) methods through an efficient and scalable neural network approach, providing new insights into the kinetics of hyperthermal products in surface reactions.