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.

A fast, fragment-based hybrid many-body interaction model is used to optimize the structures of five small-molecule organic crystals (with fixed experimental lattice parameters) and predict their lattice energies with accuracies of ∼2-4 kJ/mol compared to experiment. This model treats individual molecules in the central unit cell and their short-range two-body interactions quantum mechanically, while long-range electrostatics and many-body induction are treated with a classical polarizable force field. For the hydrogen bonded ice, formamide, and acetamide crystals, MP2 calculations extrapolated to the complete-basisset limit provide good accuracy. However, MP2 exhibits difficulties for crystals such as benzene and imidazole, where π-stacking dispersion interactions are important, and post-MP2 corrections determined from small-basis-set CCSD(T) calculations are required to achieve chemical accuracy. Using these techniques, accurate crystal lattice energy predictions for small-molecule organic crystals are feasible with currently available computing power.

Elementary processes associated with ionization of liquid water provide a framework for understanding radiation-matter interactions in chemistry and biology. Although numerous studies have been conducted on the dynamics of the hydrated electron, its partner arising from ionization of liquid water, H2O+, remains elusive. We used tunable femtosecond soft x-ray pulses from an x-ray free electron laser to reveal the dynamics of the valence hole created by strong-field ionization and to track the primary proton transfer reaction giving rise to the formation of OH. The isolated resonance associated with the valence hole (H2O+/OH) enabled straightforward detection. Molecular dynamics simulations revealed that the x-ray spectra are sensitive to structural dynamics at the ionization site. We found signatures of hydrated-electron dynamics in the x-ray spectrum.

The equation-of-motion coupled-cluster (EOM-CC) methods provide a robust description of electronically excited states and their properties. Here, we present a formalism for two-photon absorption (2PA) cross sections for the equation-of-motion for excitation energies CC with single and double substitutions (EOM-CC for electronically excited states with single and double substitutions) wave functions. Rather than the response theory formulation, we employ the expectation-value approach which is commonly used within EOM-CC, configuration interaction, and algebraic diagrammatic construction frameworks. In addition to canonical implementation, we also exploit resolution-of-the-identity (RI) and Cholesky decomposition (CD) for the electron-repulsion integrals to reduce memory requirements and to increase parallel efficiency. The new methods are benchmarked against the CCSD and CC3 response theories for several small molecules. We found that the expectation-value 2PA cross sections are within 5% from the quadratic response CCSD values. The RI and CD approximations lead to small errors relative to the canonical implementation (less than 4%) while affording computational savings. RI/CD successfully address the well-known issue of large basis set requirements for 2PA cross sections calculations. The capabilities of the new code are illustrated by calculations of the 2PA cross sections for model chromophores of the photoactive yellow and green fluorescent proteins.

Significant advances in fragment-based electronic structure methods have created a real alternative to force-field and density functional techniques in condensed-phase problems such as molecular crystals. This perspective article highlights some of the important challenges in modeling molecular crystals and discusses techniques for addressing them. First, we survey recent developments in fragment-based methods for molecular crystals. Second, we use examples from our own recent research on a fragment-based QM/MM method, the hybrid many-body interaction (HMBI) model, to analyze the physical requirements for a practical and effective molecular crystal model chemistry. We demonstrate that it is possible to predict molecular crystal lattice energies to within a couple kJ mol(-1) and lattice parameters to within a few percent in small-molecule crystals. Fragment methods provide a systematically improvable approach to making predictions in the condensed phase, which is critical to making robust predictions regarding the subtle energy differences found in molecular crystals.