We introduce a dynamical Smart Monte Carlo algorithm and assess its applicability for simulating the rare-event dynamics of adsorbate diffusion. Using the dynamical Smart Monte Carlo method, we simulate the self-diffusion of an adatom in the Cu/Cu(001) and Rh/Rh(111) systems and we compare the simulated diffusion coefficients to values arising from molecular dynamics and transition-state theory. We find that the accuracy of Smart Monte Carlo is sensitive to details of the potential-energy surface. For Cu/Cu(001), the agreement between dynamical Smart Monte Carlo, molecular dynamics, and transition-state theory is excellent. A similar comparison for the Rh/Rh(111) systems shows discrepancies between these three techniques. We find that the origins of the discrepancies in the Rh/Rh(111) system are transition-state recrossings, for small simulation time steps, and low escape rates of the adatom from the binding sites, at large time steps. We examine the sampling and dynamics in trajectories using a smaller time step for motion perpendicular to the surface than that for parallel motion. These studies show that low Smart Monte Carlo escape rates in the Rh/Rh(111) system can be correlated to excessive sampling, beyond the configurational space of the potential-energy minimum, at large time steps. Recrossings can be understood to arise from the absence of velocity correlations in the low-friction, transition-state region and can be minimized through the use of a large time step for parallel motion. With the appropriate choice of simulation time steps it is possible to improve the agreement between dynamical Smart Monte Carlo and more rigorous dynamical techniques.
All Science Journal Classification (ASJC) codes
- Physics and Astronomy(all)
- Physical and Theoretical Chemistry