Heterogeneous CPU + GPU Algorithm for Variational Two-Electron Reduced-Density Matrix-Driven Complete Active-Space Self-Consistent Field Theory

Cavity Quantum Electrodynamics Complete Active Space Configuration Interaction Theory

Polariton chemistry has attracted great attention as a potential route to modify chemical structure, properties, and reactivity through strong interactions among molecular electronic, vibrational, or rovibrational degrees of freedom. A rigorous theoretical treatment of molecular polaritons requires the treatment of matter and photon degrees of freedom on equal quantum mechanical footing. In the limit of molecular electronic strong or ultrastrong coupling to one or a few molecules, it is desirable to treat the molecular electronic degrees of freedom using the tools of ab initio quantum chemistry, yielding an approach we refer to as ab initio cavity quantum electrodynamics, where the photon degrees of freedom are treated at the level of cavity quantum electrodynamics. Here, we present an approach called Cavity Quantum Electrodynamics Complete Active Space Configuration Interaction theory to provide ground- and excited-state polaritonic surfaces with a balanced description of strong correlation effects among electronic and photonic degrees of freedom. This method provides a platform for ab initio cavity quantum electrodynamics when both strong electron correlation and strong light-matter coupling are important and is an important step toward computational approaches that yield multiple polaritonic potential energy surfaces and couplings that can be leveraged for ab initio molecular dynamics simulations of polariton chemistry.

Dynamic Correlation on the Adaptive Sampling Configuration Interaction (ASCI) Reference Function: ASCI-DSRG-MRPT2

A balanced description of static and dynamic electron correlations is at the heart of quantum chemical methods. To obtain accurate results in strongly correlated systems using wave-function-based methods, a large active space is necessary to ensure correct descriptions of static correlations. Correcting the results for dynamic correlations is also necessary. In this work, we present implementations of second-order perturbation theory for dynamic correlations based on the adaptive sampling configuration interaction self-consistent field (ASCI-SCF) method. In particular, we implemented spin-free driven similarity renormalization group second-order multireference perturbation theory (DSRG-MRPT2). The extrapolation of the ASCI + PT2 energy based on the relaxed Hamiltonian in DSRG-MRPT2 gives a reasonable approximation of DSRG-MRPT2 based on CASSCF. We demonstrate the application of the ASCI-DSRG-MRPT2 method in evaluations of the spin-state energy gaps in iron porphyrins, polyacenes, and periacenes along with the reaction energies of methane oxidation by FeO+ and electrocyclic ring formation in cethrene.

Data-Driven Refinement of Electronic Energies from Two-Electron Reduced-Density-Matrix Theory

The exponential computational cost of describing strongly correlated electrons can be mitigated by adopting a reduced-density matrix (RDM)-based description of the electronic structure. While variational two-electron RDM (v2RDM) methods can enable large-scale calculations on such systems, the quality of the solution is limited by the fact that only a subset of known necessary N-representability constraints can be applied to the 2RDM in practical calculations. Here, we demonstrate that violations of partial three-particle (T1 and T2) N-representability conditions, which can be evaluated with knowledge of only the 2RDM, can serve as physics-based features in a machine-learning (ML) protocol for improving energies from v2RDM calculations that consider only two-particle (PQG) conditions. Proof-of-principle calculations demonstrate that the model yields substantially improved energies relative to reference values from configuration-interaction-based calculations.

Distributed memory, GPU accelerated Fock construction for hybrid, Gaussian basis density functional theory

With the growing reliance of modern supercomputers on accelerator-based architecture such a graphics processing units (GPUs), the development and optimization of electronic structure methods to exploit these massively parallel resources has become a recent priority. While significant strides have been made in the development GPU accelerated, distributed memory algorithms for many modern electronic structure methods, the primary focus of GPU development for Gaussian basis atomic orbital methods has been for shared memory systems with only a handful of examples pursing massive parallelism. In the present work, we present a set of distributed memory algorithms for the evaluation of the Coulomb and exact exchange matrices for hybrid Kohn-Sham DFT with Gaussian basis sets via direct density-fitted (DF-J-Engine) and seminumerical (sn-K) methods, respectively. The absolute performance and strong scalability of the developed methods are demonstrated on systems ranging from a few hundred to over one thousand atoms using up to 128 NVIDIA A100 GPUs on the Perlmutter supercomputer.

Accelerating the density-functional tight-binding method using graphical processing units

Acceleration of the density-functional tight-binding (DFTB) method on single and multiple graphical processing units (GPUs) was accomplished using the MAGMA linear algebra library. Two major computational bottlenecks of DFTB ground-state calculations were addressed in our implementation: the Hamiltonian matrix diagonalization and the density matrix construction. The code was implemented and benchmarked on two different computer systems: (1) the SUMMIT IBM Power9 supercomputer at the Oak Ridge National Laboratory Leadership Computing Facility with 1-6 NVIDIA Volta V100 GPUs per computer node and (2) an in-house Intel Xeon computer with 1-2 NVIDIA Tesla P100 GPUs. The performance and parallel scalability were measured for three molecular models of 1-, 2-, and 3-dimensional chemical systems, represented by carbon nanotubes, covalent organic frameworks, and water clusters.

Challenges for Variational Reduced-Density-Matrix Theory: Total Angular Momentum Constraints

The variational two-electron reduced density matrix (v2RDM) method is generalized for the description of total angular momentum (J) and projection of total angular momentum (M_J) states in atomic systems described by nonrelativistic Hamiltonians, and it is shown that the approach exhibits serious deficiencies. Under ensemble N-representability constraints, v2RDM theory fails to retain the appropriate degeneracies among various J states for fixed spin (S) and orbital angular momentum (L), and for fixed L, S, and J, the manifold of M_J states is not necessarily degenerate. Moreover, a substantial energy error is observed for a system for which the two-electron reduced density matrix is exactly ensemble N-representable; in this case, the error stems from violations in pure-state N-representability conditions. Unfortunately, such violations do not appear to be good indicators of the reliability of energies from v2RDM theory in general. Several states are identified for which energy errors are near zero and yet pure-state conditions are clearly violated.

Fermi.jl: A Modern Design for Quantum Chemistry

Approximating molecular wave functions involves heavy numerical effort; therefore, codes for such tasks are written completely or partially in efficient languages such as C, C++, and Fortran. While these tools are dominant throughout quantum chemistry packages, the efficient development of new methods is often hindered by the complexity associated with code development. In order to ameliorate this scenario, some software packages take a dual approach where a simpler, higher-level language, such as Python, substitutes the traditional ones wherever performance is not critical. Julia is a novel, dynamically typed, programming language that aims to solve this two-language problem. It gained attention because of its modern and intuitive design, while still being highly optimized to compete with "low-level" languages. Recently, some chemistry-related projects have emerged exploring the capabilities of Julia. Herein, we introduce the quantum chemistry package Fermi.jl, which contains the first implementations of post-Hartree-Fock methods written in Julia. Its design makes use of many Julia core features, including multiple dispatch, metaprogramming, and interactive usage. Fermi.jl is a modular package, where new methods and implementations can be easily added to the existing code. Furthermore, it is designed to maximize code reusability by relying on general functions with specialized methods for particular cases. The feasibility of the project is explored through evaluating the performance of popular ab initio methods. It is our hope that this project motivates the usage of Julia within the community and brings new contributions into Fermi.jl.

Challenges for variational reduced-density-matrix theory with three-particle N-representability conditions

The direct variational optimization of the two-electron reduced density matrix (2RDM) can provide a reference-independent description of the electronic structure of many-electron systems that naturally capture strong or nondynamic correlation effects. Such variational 2RDM approaches can often provide a highly accurate description of strong electron correlation, provided that the 2RDMs satisfy at least partial three-particle N-representability conditions (e.g., the T2 condition). However, recent benchmark calculations on hydrogen clusters [N. H. Stair and F. A. Evangelista, J. Chem. Phys. 153, 104108 (2020)] suggest that even the T2 condition leads to unacceptably inaccurate results in the case of two- and three-dimensional clusters. We demonstrate that these failures persist under the application of full three-particle N-representability conditions (3POS). A variety of correlation metrics are explored in order to identify regimes under which 3POS calculations become unreliable, and we find that the relative squared magnitudes of the cumulant three- and two-particle reduced density matrices correlate reasonably well with the energy error in these systems. However, calculations on other molecular systems reveal that this metric is not a universal indicator for the reliability of the reduced-density-matrix theory with 3POS conditions.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

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.

Near-Exact CASSCF-Level Geometry Optimization with a Large Active Space using Adaptive Sampling Configuration Interaction Self-Consistent Field Corrected with Second-Order Perturbation Theory (ASCI-SCF-PT2)

An accurate description of electron correlation is one of the most challenging problems in quantum chemistry. The exact electron correlation can be obtained by means of full configuration interaction (FCI). A simple strategy for approximating FCI at a reduced computational cost is selected CI (SCI), which diagonalizes the Hamiltonian within only the chosen configuration space. Recovery of the contributions of the remaining configurations is possible with second-order perturbation theory. Here, we apply adaptive sampling configuration interaction (ASCI) combined with molecular orbital optimizations (ASCI-SCF) corrected with second-order perturbation theory (ASCI-SCF-PT2) for geometry optimization by implementing the analytical nuclear gradient algorithm for ASCI-PT2 with the Z-vector (Lagrangian) formalism. We demonstrate that for phenalenyl radicals and anthracene, optimized geometries and the number of unpaired electrons can be obtained at nearly the CASSCF accuracy by incorporating PT2 corrections and extrapolating them. We demonstrate the current algorithm's utility for optimizing the equilibrium geometries and electronic structures of six-ring-fused polycyclic aromatic hydrocarbons and 4-periacene.

Pushing the Limits of Acene Chemistry: The Recent Surge of Large Acenes

Acenes, consisting of linearly fused benzene rings, are an important fundamental class of organic compounds with various applications. Hexacene is the largest acene that was synthesized and isolated in the 20th century. The next largest member of the acene family, heptacene, was observed in 2007 and since then significant progress in preparing acenes has been reported. Significantly larger acenes, up to undecacene, could be studied by means of low-temperature matrix isolation spectroscopy with in situ photolytic generation, and up to dodecacene by means of on-surface synthesis employing innovative precursors and highly defined crystalline metal surfaces under ultrahigh vacuum conditions. The review summarizes recent experimental and theoretical advances in the area of acenes that give a significantly deeper insight into the fundamental properties and nature of the electronic structure of this fascinating class of organic compounds.

TAO-DFT-Based Ab Initio Molecular Dynamics

Recently, AIMD (ab initio molecular dynamics) has been extensively employed to explore the dynamical information of electronic systems. However, it remains extremely challenging to reliably predict the properties of nanosystems with a radical nature using conventional electronic structure methods (e.g., Kohn-Sham density functional theory) due to the presence of static correlation. To address this challenge, we combine the recently formulated TAO-DFT (thermally-assisted-occupation density functional theory) with AIMD. The resulting TAO-AIMD method is employed to investigate the instantaneous/average radical nature and infrared spectra of n-acenes containing n linearly fused benzene rings (n = 2-8) at 300 K. According to the TAO-AIMD simulations, on average, the smaller n-acenes (up to n = 5) possess a nonradical nature, and the larger n-acenes (n = 6-8) possess an increasing radical nature, showing remarkable similarities to the ground-state counterparts at 0 K. Besides, the infrared spectra of n-acenes obtained with the TAO-AIMD simulations are in qualitative agreement with the existing experimental data.

TAO-DFT investigation of electronic properties of linear and cyclic carbon chains

It has been challenging to adequately investigate the properties of nanosystems with radical nature using conventional electronic structure methods. We address this challenge by calculating the electronic properties of linear carbon chains (l-CC[n]) and cyclic carbon chains (c-CC[n]) with n = 10-100 carbon atoms, using thermally-assisted-occupation density functional theory (TAO-DFT). For all the cases investigated, l-CC[n]/c-CC[n] are ground-state singlets, and c-CC[n] are energetically more stable than l-CC[n]. The electronic properties of l-CC[n]/c-CC[n] reveal certain oscillation patterns for smaller n, followed by monotonic changes for larger n. For the smaller carbon chains, odd-numbered l-CC[n] are more stable than the adjacent even-numbered ones; c-CC[[Formula: see text]]/c-CC[4m] are more/less stable than the adjacent odd-numbered ones, where m are positive integers. As n increases, l-CC[n]/c-CC[n] possess increasing polyradical nature in their ground states, where the active orbitals are delocalized over the entire length of l-CC[n] or the whole circumference of c-CC[n].

Analytic gradients for spin multiplets in natural orbital functional theory

Analytic energy gradients with respect to nuclear motion are derived for non-singlet compounds in the natural orbital functional theory. We exploit the formulation for multiplets in order to obtain a simple formula valid for any many-electron system in its ground mixed state with a total spin S and all possible spin projection S_z values. We demonstrate that the analytic gradients can be obtained without resorting to linear response theory or involving iterative procedures. A single evaluation is required, so integral derivatives can be computed on-the-fly along the calculation, thus improving the effectiveness of screening by the Schwarz inequality. The results for small- and medium-sized molecules with many spin multiplicities are shown. Our results are compared with the experimental data and accurate theoretical equilibrium geometries.

Reduced Density Matrix-Driven Complete Active Apace Self-Consistent Field Corrected for Dynamic Correlation from the Adiabatic Connection

The recently proposed multireference adiabatic connection (AC) formalism [Pernal, Phys. Rev. Lett. 120, 013001 (2018)] is applied to recover dynamic electron correlation effects lacking in variational two-electron reduced density matrix (v2RDM)-driven complete active space self-consistent field theory (CASSCF). The AC approach is validated by computing potential energy curves for the dissociation of molecular nitrogen and the symmetric double dissociation of H₂O while enforcing two sets of approximate N-representability conditions in the underlying v2RDM-driven CASSCF calculations (either two-particle or two-particle plus partial three-particle conditions). The AC yields smaller absolute errors than second-order N-electron perturbation theory (NEVPT2) at all molecular geometries for both sets of the N-representability conditions considered. The efficacy of the approach for thermochemistry is also assessed for a set of 31 small-molecule reactions. When imposing partial three-particle N-representability conditions, mean and maximum unsigned errors in reaction energies from the AC are superior to those from NEVPT2.

TAO-DFT Study on the Electronic Properties of Diamond-Shaped Graphene Nanoflakes

At the nanoscale, it has been rather troublesome to properly explore the properties associated with electronic systems exhibiting a radical nature using traditional electronic structure methods. Graphene nanoflakes, which are graphene nanostructures of different shapes and sizes, are typical examples. Recently, TAO-DFT (i.e., thermally-assisted-occupation density functional theory) has been formulated to tackle such challenging problems. As a result, we adopt TAO-DFT to explore the electronic properties associated with diamond-shaped graphene nanoflakes with n = 2–15 benzenoid rings fused together at each side, designated as n-pyrenes (as they could be expanded from pyrene). For all the n values considered, n-pyrenes are ground-state singlets. With increasing the size of n-pyrene, the singlet-triplet energy gap, vertical ionization potential, and fundamental gap monotonically decrease, while the vertical electron affinity and symmetrized von Neumann entropy (which is a quantitative measure of radical nature) monotonically increase. When n increases, there is a smooth transition from the nonradical character of the smaller n-pyrenes to the increasing polyradical nature of the larger n-pyrenes. Furthermore, the latter is shown to be related to the increasing concentration of active orbitals on the zigzag edges of the larger n-pyrenes.

CASSCF with Extremely Large Active Spaces Using the Adaptive Sampling Configuration Interaction Method

The complete active space self-consistent field (CASSCF) method is the principal approach employed for studying strongly correlated systems. However, exact CASSCF can only be performed on small active spaces of ∼20 electrons in ∼20 orbitals due to exponential growth in the computational cost. We show that employing the Adaptive Sampling Configuration Interaction (ASCI) method as an approximate Full CI solver in the active space allows CASSCF-like calculations within chemical accuracy (<1 kcal/mol for relative energies) in active spaces with more than ∼50 active electrons in ∼50 active orbitals, significantly increasing the sizes of systems amenable to accurate multiconfigurational treatment. The main challenge with using any selected CI-based approximate CASSCF is the orbital optimization problem; they tend to exhibit large numbers of local minima in orbital space due to their lack of invariance to active-active rotations (in addition to the local minima that exist in exact CASSCF). We highlight methods that can avoid spurious local extrema as a practical solution to the orbital optimization problem. We employ ASCI-SCF to demonstrate a lack of polyradical character in moderately sized periacenes with up to 52 correlated electrons and compare against heat-bath CI on an iron porphyrin system with more than 40 correlated electrons.