This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange–correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear–electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an “open teamware” model and an increasingly modular design.

We present a current, up-to-date review of the surface hopping methodology for solving nonadiabatic problems, 25 years after Tully published the fewest switches surface hopping algorithm. After reviewing the original motivation for and failures of the algorithm, we give a detailed examination of modern advances, focusing on both theoretical and practical issues. We highlight how one can partially derive surface hopping from the Schrödinger equation in the adiabatic basis, how one can change basis within the surface hopping algorithm, and how one should understand and apply the notions of decoherence and wavepacket bifurcation. The question of time reversibility and detailed balance is also examined at length. Recent applications to photoexcited conjugated polymers are discussed briefly.

We study several recent mean-field semiclassical dynamics methods, focusing on the ability to recover detailed balance for long time (equilibrium) populations. We focus especially on Miller and Cotton's [J. Phys. Chem. A 117, 7190 (2013)] suggestion to include both zero point electronic energy and windowing on top of Ehrenfest dynamics. We investigate three regimes: harmonic surfaces with weak electronic coupling, harmonic surfaces with strong electronic coupling, and anharmonic surfaces with weak electronic coupling. In most cases, recent additions to Ehrenfest dynamics are a strong improvement upon mean-fieldtheory. However, for methods that include zero point electronic energy, we show that anharmonic potential energy surfaces often lead to numerical instabilities, as caused by negative populations and forces. We also show that, though the effect of negative forces can appear hidden in harmonic systems, the resulting equilibrium limits do remain dependent on any windowing and zero point energy parameters.

We present a preliminary extension of the fewest switches surface hopping (FSSH) algorithm to the case of complex Hamiltonians as appropriate for modeling the dynamics of photoexcited molecules in magnetic fields. We make ansätze for the direction of momentum rescaling and we account for Berry's phase effects through "magnetic" forces as applicable in the adiabatic limit. Because Berry's phase is a nonlocal, topological characteristic of a set of entangled potential energy surfaces, we find that Tully's local FSSH algorithm can only partially capture the correct physics.

Although the quantum classical Liouville equation (QCLE) arises by cutting off the exact equation of motion for a coupled nuclear-electronic system at order 1 (1 = 0 ), we show that the QCLE does include Berry's phase effects and Berry's forces (which are proportional to a higher order, = 1 ). Thus, the fundamental equation underlying mixed quantum-classical dynamics does not need a correction for Berry's phase effects and is valid for the case of complex Hamiltonians. Furthermore, we also show that, even though Tully's surface hopping model ignores Berry's phase, Berry's phase effects are included automatically within Ehrenfest dynamics. These findings should be of great importance if we seek to model coupled nuclearelectronic dynamics for systems with spin-orbit coupling, where the complex nature of the Hamiltonian is