The Chronus Quantum (ChronusQ) software package is an open source (under the GNU General Public License v2) software infrastructure which targets the solution of challenging problems that arise in ab initio electronic structure theory. Special emphasis is placed on the consistent treatment of time dependence and spin in the electronic wave function, as well as the inclusion of relativistic effects in said treatments. In addition, ChronusQ provides support for the inclusion of uniform finite magnetic fields as external perturbations through the use of gauge‐including atomic orbitals. ChronusQ is a parallel electronic structure code written in modern C++ which utilizes both message passing implementation and shared memory (OpenMP) parallelism. In addition to the examination of the current state of code base itself, a discussion regarding ongoing developments and developer contributions will also be provided. This article is categorized under: Software > Quantum Chemistry Electronic Structure Theory > Ab Initio Electronic Structure Methods Electronic Structure Theory > Density Functional Theory

The quantum mechanical treatment of both electrons and nuclei is crucial in nonadiabatic dynamical processes such as proton-coupled electron transfer. The nuclear−electronic orbital (NEO) method provides an elegant framework for including nuclear quantum effects beyond the Born–Oppenheimer approximation. To enable the study of nonequilibrium properties, we derive and implement a real-time NEO (RT-NEO) approach based on time-dependent Hatree-Fock or density functional theory, in which the electronic and nuclear degrees of freedom are propagated in a time-dependent variational framework. Nuclear and electronic spectral features can be resolved from the time-dependent dipole moment computed using the RT-NEO method. The test cases show the dynamical interplay between the quantum nuclei and the electrons through vibronic coupling. Moreover, vibrational excitation in the RT-NEO approach is demonstrated by applying a resonant driving field, and electronic excitation is demonstrated by simulating excited state intramolecular proton transfer. This work shows that the RT-NEO approach is a promising tool to study nonadiabatic quantum dynamical processes within a time-dependent variational description for the coupled electronic and nuclear degrees of freedom.

The recent development of the Ehrenfest dynamics approach in the nuclear-electronic orbital (NEO) framework provides a promising way to simulate coupled nuclear-electronic dynamics. Our previous study showed that the NEO-Ehrenfest approach with a semiclassical traveling proton basis method yields accurate predictions of molecular vibrational frequencies. In this work, we provide a more thorough analysis of the semiclassical traveling proton basis method to elucidate its validity and convergence behavior. We also conduct NEO-Ehrenfest dynamics simulations to study an excited state intramolecular proton transfer process. These simulations reveal that nuclear quantum effects influence the predictions of proton transfer reaction rates and kinetic isotope effects due to the intrinsic delocalized nature of the quantum nuclear wave function. This work illustrates the importance of nuclear quantum effects in coupled nuclear-electronic dynamical processes and shows that the NEO-Ehrenfest approach can be a powerful tool for providing insights and predictions for these processes.

The block-localized wave function method is useful to provide insights on chemical bonding and intermolecular interactions through energy decomposition analysis. The method relies on block localization of molecular orbitals (MOs) by constraining the orbitals to basis functions within given blocks. Here, a generalized block-localized orbital (GBLO) method is described to allow both physically localized and delocalized MOs to be constrained in orbital-block definitions. Consequently, GBLO optimization can be conveniently tailored by imposing specific constraints. The GBLO method is illustrated by three examples: (1) constrained polarization response orbitals through dipole and quadrupole perturbation in a water dimer complex, (2) the ground and first excited-state potential energy curves of ethene about its C–C bond rotation, and (3) excitation energies of double electron excited states. Multistate density functional theory is used to determine the energies of the adiabatic ground and excited states using a minimal active space (MAS) comprising specifically charge-constrained and excited determinant configurations that are variationally optimized by the GBLO method. We find that the GBLO expansion that includes delocalized MOs in configurational blocks significantly reduces computational errors in comparison with physical block localization, and the computed ground- and excited-state energies are in good accordance with experiments and results obtained from multireference configuration interaction calculations.

The recently developed real-time nuclear–electronic orbital (RT-NEO) approach provides an elegant framework for treating electrons and selected nuclei, typically protons, quantum mechanically in nonequilibrium dynamical processes. However, the RT-NEO approach neglects the motion of the other nuclei, preventing a complete description of the coupled nuclear–electronic dynamics and spectroscopy. In this work, the dynamical interactions between the other nuclei and the electron–proton subsystem are described with the mixed quantum–classical Ehrenfest dynamics method. The NEO-Ehrenfest approach propagates the electrons and quantum protons in a time-dependent variational framework, while the remaining nuclei move classically on the corresponding average electron–proton vibronic surface. This approach includes the non-Born–Oppenheimer effects between the electrons and the quantum protons with RT-NEO and between the classical nuclei and the electron–proton subsystem with Ehrenfest dynamics. Spectral features for vibrational modes involving both quantum and classical nuclei are resolved from the time-dependent dipole moments. This work shows that the NEO-Ehrenfest method is a powerful tool to study dynamical processes with coupled electronic and nuclear degrees of freedom.