In order to study electron-transfer mediated chemical processes on a metal surface, one requires not one but two potential energy surfaces (one ground state and one excited state) as in Marcus theory. In this letter, we report that a novel, dynamically weighted, state-averaged constrained CASSCF(2,2) (DW-SA-cCASSCF(2,2)) can produce such surfaces for the Anderson impurity model. Both ground and excited state potentials are smooth, they incorporate states with a charge transfer character, and the accuracy of the ground state surface can be verified for some model problems by renormalization group theory. Future development of gradients and nonadiabatic derivative couplings should allow for the study of nonadiabatic dynamics for molecules near metal surfaces.

We present an efficient set of methods for propagating excited-state dynamics involving a large number of electronic states based on a configuration interaction singles (CIS) electronic state overlap scheme. Specifically, (i) following Head-Gordon et al we implement an exact evaluation of the overlap of singly-excited electronic states at different nuclear geometries using a biorthogonal basis and (ii) we employ a unified protocol 1 for choosing the correct phase for each adiabat at each geometry. For many-electron systems, the combination of these techniques significantly reduces the computational cost of integrating the electronic Schrodinger equation and imposes minimal overhead on top of the underlying electronic structure calculation. As a demonstration, we calculate the electronic excited-state dynamics for a hydrogen molecule scattering off a silver metal cluster, focusing on high-lying excited states where many electrons can be excited collectively and crossings are plentiful. Interestingly, we find that the high-lying, plasmon-like collective excitation spectrum changes with nuclear dynamics, highlighting the need to simulate non-adiabatic nuclear dynamics and plasmonic excitations simultaneously. In the future, the combination of methods presented here should help theorists build a mechanistic understanding of plasmon-assisted charge transfer and excitation energy relaxation processes near a nanoparticle or metal surface.

We derive and implement the necessary equations for solving a dynamically weighted, state-averaged constrained CASSCF(2,2) wavefunction describing a molecule on a metal surface. We show that a partial constraint is far more robust than a full constraint. We further calculate the system-bath electronic couplings that arise because, near a metal, there is a continuum (rather than discrete) number of electronic states. This approach should be very useful for simulating heterogeneous electron transfer going forward.

We test a set of multiconfigurational wavefunction approaches for calculating the ground state electron population for a two-site Anderson model representing a molecule on a metal surface. In particular, we compare (i) a Hartree Fock like wavefunction where frontier orbitals are allowed to be nonorthogonal versus (ii) a fully non-orthogonal configuration interaction wavefunction based on constrained Hartree–Fock states. We test both the strong and weak metal-molecule hybridization (Γ) limits as well as the strong and weak electron–electron repulsion (U) limits. We obtain accurate results as compared with exact numerical renormalization group theory, recovering charge transfer states where appropriate. The current framework should open a path to run molecular non-adiabatic dynamics on metal surfaces.

We derive and implement the necessary equations for solving a dynamically weighted, state-averaged constrained CASSCF(2,2) wave function describing a molecule on a metal surface, where we constrain the overlap between two active orbitals and the impurity atomic orbitals to be a finite number. We show that a partial constraint is far more robust than a full constraint. We further calculate the system-bath electronic couplings that arise because, near a metal, there is a continuum (rather than discrete) number of electronic states. This approach should be very useful for simulating heterogeneous electron transfer and electrochemical dynamics going forward.