Determining the range of validity of Migdal's approximation for electron-phonon (e-ph) coupled systems is a long-standing problem. Many attempts to answer this question employ the Holstein Hamiltonian, where the electron density couples linearly to local lattice displacements. When these displacements are large, however, nonlinear corrections to the interaction must also be included, which can significantly alter the physical picture obtained from this model. Using determinant quantum Monte Carlo and the self-consistent Migdal approximation, we compared superconducting and charge-density-wave correlations in the Holstein model with and without second-order nonlinear interactions. We find a disagreement between the two cases, even for relatively small values of the e-ph coupling strength, and, importantly, that this can occur in the same parameter regions where Migdal's approximation holds. Our results demonstrate that questions regarding the validity of Migdal's approximation go hand in hand with questions of the validity of a linear e-ph interaction.

We describe a new open-source Python-based package for high accuracy correlated electron calculations using quantum Monte Carlo (QMC) in real space: PyQMC. PyQMC implements modern versions of QMC algorithms in an accessible format, enabling algorithmic development and easy implementation of complex workflows. Tight integration with the PySCF environment allows for a simple comparison between QMC calculations and other many-body wave function techniques, as well as access to high accuracy trial wave functions.

No abstract

We demonstrate the applicability of extended Lagrangian Born-Oppenheimer quantumbased molecular dynamics (XL-BOMD) to model electron transfer reactions occurring on solidliquid interfaces. Specifically, we consider the reduction of O 2 as catalyzed at the interface of an N-doped graphene sheet and H 2 O at fuel cell cathodes. This system is a good testbed for nextgeneration computational chemistry methods since the electrochemical functionalities strongly depend on atomic-scale quantum mechanics. As opposed to prior iterations of first principles molecular dynamics, XL-BOMD only requires a full self-consistent-charge relaxation during the initial time step. The electronic ground state and total energy are stabilized thereafter through nuclear and electronic equations of motion assisted by an inner-product kernel updated with low-rank approximations. A species charge analysis reveals that the kernel-based XL-BOMD simulation can capture an electron transfer between the PGM-free catalyst and a solvated O 2 molecule mediated by H 2 O, which results in the molecular dissociation of O 2 .

No abstract