Orbital Surface Hopping with an Electron Thermostat Yields Accurate Dynamics and Detailed Balance

Abstract

In mixed quantum-classical simulations of molecule-metal surface interactions, the discretization of the metallic electronic continuum typically results in a closed-system representation that fails to capture the open-system nature of the true physical process. This approximation can introduce significant artifacts, including deviations in the dynamical evolution and a violation of the principle of detailed balance. To address this fundamental challenge, we introduce an electronic thermostat into our previously developed orbital surface hopping (OSH) framework, generalizing the method to efficiently handle many discrete electronic states. We first outline the derivation of electronic thermostat orbital surface hopping, where the amplitude of the electronic thermostat is well justified. We then demonstrate that this method can reproduce accurate dynamics and detailed balance in long time, whereas without electronic thermostat the detailed balance is violated. Thus, this method offers a reliable tool for studying nonadiabatic dynamics near metal surfaces.

References in corpus (9)

• Scattering theory of current-induced forces in mesoscopic systems
• Current-induced atomic dynamics, instabilities, and Raman signals: Quasi-classical Langevin equation approach
• Born-Oppenheimer Dynamics, Electronic Friction, and the Inclusion of Electron-Electron Interactions
• Surface hopping from the perspective of quantum-classical Liouville dynamics
• Nuclear dynamics at molecule-metal interfaces: A pseudoparticle perspective
• HierarchicalEOM.jl: An efficient Julia framework for hierarchical equations of motion in open quantum systems
• Efficient implementation and performance analysis of the independent electron surface hopping method for dynamics at metal surfaces
• QuantumToolbox.jl: An efficient Julia framework for simulating open quantum systems
• Fewest switches surface hopping with decoherence in the Marcus inverted regime: correct rates but wrong thermal populations