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 aim to observe a chemical reaction in real time using gas-phase X-ray diffraction. In our initial experiment at the Linac Coherent Light Source (LCLS), we investigated the model system 1,3-cyclohexadiene (CHD) at very low vapor pressures. This reaction serves as a benchmark for numerous transformations in organic synthesis and natural product biology. Excitation of CHD by an ultraviolet optical pulse initiates an electrocyclic reaction that transforms the closed ring system into the open-chain structure of 1,3,5-hexatriene. We describe technical points of the experimental method and present first results. We also outline an approach to analyze the data involving nonlinear least-square optimization routines that match the experimental observations with predicted diffraction patterns calculated from trajectories for nonadiabatic vibronic wave packets.

We describe a gas-phase x-ray scattering experiment capable of capturing molecular motions with atomic spatial resolution and femtosecond time resolution. X-ray free electron lasers can deliver intense x-ray pulses of ultrashort duration, making them suitable to study ultrafast chemical reaction dynamics in an ultraviolet pump, x-ray probe scheme. A cell diffractometer balances sample flow with gas density and laser focusing conditions to provide adequate scattering vector resolution with high signal intensity and near-uniform excitation probability. Images from a pixel-array x-ray detector, spatially and electronically calibrated, allow for detection of scattering intensity changes below 1%. First experiments on the ring-opening reaction of 1,3-cyclohexadiene to form 1, 3, 5-hexatriene show a rapid initial reaction on an 80 fs time scale.

The accurate description of large molecular systems in complex environments remains an ongoing challenge for the field of computational chemistry. This problem is even more pronounced for photoinduced processes, as multiple excited electronic states and their corresponding nonadiabatic couplings must be taken into account. Multiscale approaches such as hybrid quantum mechanics/molecular mechanics (QM/MM) offer a balanced compromise between accuracy and computational burden. Here, we introduce an open-source software package (INAQS) for nonadiabatic QM/MM simulations that bridges the sampling capabilities of the GROMACS MD package and the excited-state infrastructure of the Q-CHEM electronic structure software. The interface is simple and can be adapted easily to other MD codes. The code supports a variety of different trajectory-based molecular dynamics, ranging from Born–Oppenheimer to surface hopping dynamics. To illustrate the power of this combination, we simulate electronic absorption spectra, free-energy surfaces along a reaction coordinate, and the excited-state dynamics of 1,3-cyclohexadiene in solution.

With the notable exception of some illustrative two-degree-of-freedom models whose surprising classical dynamics has been worked out in detail, theories of roaming have largely bypassed the issue of when and why the counterintuitive phenomenon of roaming occurs. We propose that a useful way to begin to address these issues is to look for the geodesic (most efficient) pathways through the potential surfaces of candidate systems. Although roaming manifests itself in an unusual behavior at asymptotic geometries, we found in the case of formaldehyde dissociation that it was the pathways traversing the parts of the potential surface corresponding to highly vibrationally excited reactants that were the most revealing. An examination of the geodesics for roaming pathways in this region finds that they are much less tightly defined than the geodesics in that same region that lead directly to dissociation (whether into closed-shell products or into radical products). Thus, the broader set of options available to the roaming channel gives it an entropic advantage over more conventional reaction channels. These observations suggest that what leads to roaming in other systems may be less the presence of a localized "roaming transition state," than the existence of an entire region of the potential surface conducive to multiple equivalent pathways.

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 present an efficient set of methods for propagating excited-state dynamics involving a large number of electronic states based on a 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 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.

No abstract