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.

Anisotropy is a fundamental property of particle interactions. It occupies a central role in cold and ultra-cold molecular processes, where long-range forces have been found to significantly depend on orientation in ultra-cold polar molecule collisions 1,2 . Recent experiments have demonstrated the emergence of quantum phenomena such as scattering resonances in the cold collisions regime due to quantization of the intermolecular degrees of freedom 3-8 . Although these states have been shown to be sensitive to interaction details, the effect of anisotropy on quantum resonances has eluded experimental observation so far. Here, we directly measure the anisotropy in atom-molecule interactions via quantum resonances by changing the quantum state of the internal molecular rotor. We observe that a quantum scattering resonance at a collision energy of appears in the Penning ionization of molecular hydrogen with metastable helium only if the molecule is rotationally excited. We use state of the art ab initio and multichannel quantum molecular dynamics calculations to show that the anisotropy contributes to the effective interaction only for molecules in the first excited rotational state, whereas rotationally ground state interacts purely isotropically with metastable helium. Control over the quantum state of the internal molecular rotation allows us to switch the anisotropy on or off and thus disentangle the isotropic and anisotropic parts of the interaction. These quantum phenomena provide a challenging benchmark for even the most advanced theoretical descriptions, highlighting the advantage of using cold collisions to advance the microscopic understanding of particle interactions.

The role of internal molecular degrees of freedom, such as rotation, has scarcely been explored experimentally in low-energy collisions despite their significance to cold and ultracold chemistry. Particularly important to astrochemistry is the case of the most abundant molecule in interstellar space, hydrogen, for which two spin isomers have been detected, one of which exists in its rotational ground state whereas the other is rotationally excited. Here we demonstrate that quantization of molecular rotation plays a key role in cold reaction dynamics, where rotationally excited ortho-hydrogen reacts faster due to a stronger long-range attraction. We observe rotational state-dependent non-Arrhenius universal scaling laws in chemi-ionization reactions of para-H2 and ortho-H2 by He(2(3)P2), spanning three orders of magnitude in temperature. Different scaling laws serve as a sensitive gauge that enables us to directly determine the exact nature of the long-range intermolecular interactions. Our results show that the quantum state of the molecular rotor determines whether or not anisotropic long-range interactions dominate cold collisions.

State-of-the-art ab initio techniques have been applied to compute the potential energy curves for the electronic states in the A(1)Σ(u)(+), c(3)Π(u), and a(3)Σ(u)(+) manifold of the strontium dimer, the spin-orbit and nonadiabatic coupling matrix elements between the states in the manifold, and the electric transition dipole moment from the ground X(1)Σ(g)(+) to the nonrelativistic and relativistic states in the A+c+a manifold. The potential energy curves and transition moments were obtained with the linear response (equation of motion) coupled cluster method limited to single, double, and linear triple excitations for the potentials and limited to single and double excitations for the transition moments. The spin-orbit and nonadiabatic coupling matrix elements were computed with the multireference configuration interaction method limited to single and double excitations. Our results for the nonrelativistic and relativistic (spin-orbit coupled) potentials deviate substantially from recent ab initio calculations. The potential energy curve for the spectroscopically active (1)0(u)(+) state is in quantitative agreement with the empirical potential fitted to high-resolution Fourier transform spectra [A. Stein, H. Knöckel, and E. Tiemann, Eur. Phys. J. D 64, 227 (2011)]. The computed ab initio points were fitted to physically sound analytical expressions, and used in converged coupled channel calculations of the rovibrational energy levels in the A+c+a manifold and line strengths for the A(1)Σ(u)(+)←X(1)Σ(g (+) transitions. Positions and lifetimes of quasi-bound Feshbach resonances lying above the (1)S(0) + (3)P(1) dissociation limit were also obtained. Our results reproduce (semi)quantitatively the experimental data observed thus far. Predictions for on-going and future experiments are also reported.

Weakly bound molecules have physical properties without atomic analogues, even as the bond length approaches dissociation. In particular, the internal symmetries of homonuclear diatomic molecules result in formation of two-body superradiant and subradiant excited states. While superradiance [1][2][3] has been demonstrated in a variety of systems, subradiance [4][5][6] is more elusive due to the inherently weak interaction with the environment. Here we characterize the properties of deeply subradiant molecular states with intrinsic quality factors exceeding 10 13 via precise optical spectroscopy with the longest molecule-light coherent interaction times to date. We find that two competing effects limit the lifetimes of the subradiant molecules, with different asymptotic behaviors. The first is radiative decay via weak magnetic-dipole and electric-quadrupole interactions. We prove that its rate increases quadratically with the bond length, confirming quantum mechanical predictions. The second is nonradiative decay through weak gyroscopic predissociation, with a rate proportional to the vibrational mode spacing and sensitive to short-range physics. This work bridges the gap between atomic and molecular metrology based on lattice-clock techniques [7], yielding new understanding of long-range interatomic interactions and placing ultracold molecules at the forefront of precision measurements.Simple molecules provide a wealth of opportunities for precision measurements. Their richer internal structure compared to atoms enables experiments that push the boundaries in determinations of the electric dipole moment of the electron [8], the electron-to-proton mass ratio and its variations [9,10], and parity violation [11]. Diatomic molecules are moving to the forefront of manybody science [12] and quantum chemistry [13], providing glimpses into fundamental laws [14]. However, this attractive complexity of molecular structure has historically posed difficulties for manipulation and modeling [15]. This work removes many of these barriers by employing techniques of optical lattice atomic clocks [16,17] to control the quantum states of weakly bound homonuclear diatomic strontium molecules, in particular by using state-insensitive optical lattices [18] for molecular transitions with three types of optical transition moments. We observe strongly forbidden optical transitions in this asymptotic diatomic system, an ideal regime for studying the breakdown of the ubiquitous dipole approximation where the size of the quantum particle is a significant fraction of the resonant wavelength. We explain these observations with a state-of-the-art ab initio molecular model [19] and asymptotic scaling laws. The results prove that the quantum mechanical effect of subradiance can be exploited for precision spectroscopy, and demonstrate the promise of combining precise state control, coherent manipulation, and accurate ab initio calculations with recently available ultracold molecular systems.We create Sr 2 molecules by photoassociation [20] from ...