Active-space two-electron reduced-density-matrix method: Complete active-space calculations without diagonalization of the N-electron Hamiltonian

Effect of Oriented External Electric Fields on the Electronic Properties of Linear Acenes: A Thermally Assisted Occupation DFT Study

Recently, oriented external electric fields (OEEFs) have earned much attention due to the possibility of tuning the properties of electronic systems. From a theoretical perspective, one can resort to electronic structure calculations to understand how the direction and strength of OEEFs affect the properties of electronic systems. However, for multi-reference (MR) systems, calculations employing the popular Kohn-Sham density functional theory with the traditional semilocal and hybrid exchange-correlation energy functionals can yield erroneous results. Owing to its decent compromise between accuracy and efficiency for MR systems at the nanoscale (i.e., MR nanosystems), in this study, thermally assisted occupation density functional theory (TAO-DFT) is adopted to explore the electronic properties of n-acenes (n = 2-10), containing n linearly fused benzene rings, in OEEFs, where the OEEFs of various electric field strengths are applied along the long axes of n-acenes. According to our TAO-DFT calculations, the ground states of n-acenes in OEEFs are singlets for all the cases examined. The effect of OEEFs is shown to be significant on the vertical ionization potentials and vertical electron affinities of ground-state n-acenes with odd-number fused benzene rings. Moreover, the MR character of ground-state n-acenes in OEEFs increases with the increase in the acene length and/or the electric field strength.

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.

Electronic Properties of Graphene Nano-Parallelograms: A Thermally Assisted Occupation DFT Computational Study

In this computational study, we investigate the electronic properties of zigzag graphene nano-parallelograms (GNPs), which are parallelogram-shaped graphene nanoribbons of various widths and lengths, using thermally assisted occupation density functional theory (TAO-DFT). Our calculations revealed a monotonic decrease in the singlet-triplet energy gap as the GNP length increased. The GNPs possessed singlet ground states for all the cases examined. With the increase of GNP length, the vertical ionization potential and fundamental gap decreased monotonically, while the vertical electron affinity increased monotonically. Some of the GNPs studied were found to possess fundamental gaps in the range of 1-3 eV, lying in the ideal region relevant to solar energy applications. Besides, as the GNP length increased, the symmetrized von Neumann entropy increased monotonically, denoting an increase in the degree of the multi-reference character associated with the ground state GNPs. The occupation numbers and real-space representation of active orbitals indicated that there was a transition from the nonradical nature of the shorter GNPs to the increasing polyradical nature of the longer GNPs. In addition, the edge/corner localization of the active orbitals was found for the wider and longer GNPs.

Single Reference Treatment of Strongly Correlated H4 and H10 Isomers with Richardson-Gaudin States

Richardson-Gaudin (RG) states are employed as a variational wave function ansatz for strongly correlated isomers of H4 and H10. In each case, a single RG state describes the seniority-zero sector quite well. Simple natural orbital functionals offer a cheap and reasonable approximation of the outstanding weak correlation in the seniority-zero sector, while systematic improvement is achieved by performing a configuration interaction in terms of RG states.

Cyclo[n]carbons and catenanes from different perspectives: disentangling the molecular thread

All-carbon atomic rings, cyclo[n]carbons, have recently attracted vivid attention of experimentalists and theoreticians. Among them, cyclo[18]carbon is the most studied system. In this paper, we summarize and review various properties of cyclo[n]carbons, emphasising the aspects of their aromaticity/antiaromaticity. In the first part, the trends in bonding patterns and selected aromaticity indices with the increasing size of the rings are discussed. In the second part we explore the properties of catenane models based on interlocked cyclo[18]carbon rings from different perspectives and investigate their behaviour under the action of external force using computational experiments.

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.

Characterizing Excited States of a Copper-Based Molecular Qubit Candidate with Correlated Electronic Structure Methods

Molecular spins have a variety of potential advantages as qubits for quantum computation, such as tunability and well-understood design pathways through organometallic synthesis. Organometallic and heavy-metal-based molecular spin qubits can also exhibit rich electronic structures due to ligand field interactions and electron correlation. These features make consistent and reliable modeling of these species a considerable challenge for contemporary electronic structure techniques. Here, we elucidate the electronic structure of a Cu(II) complex analogous to a recently proposed room-temperature molecular spin qubit. Using active space methods to describe the electron correlation, we show the nuanced interaction between the metal d orbitals and ligand σ and π orbitals makes these systems challenging to model, both in terms of the delocalized spin density and the excited state ordering. We show that predicting the correct spin delocalization requires special consideration of the Cu d orbitals and that the excited state spectrum for the Cu(III) complex also requires the explicit inclusion of the π orbitals in the active space. These interactions are rather common in molecular spin qubit motifs and may play an important role in spin-decoherence processes. Our results may lend insight into future studies of the orbital interactions and electron delocalization of similar complexes.

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.

Measuring Electron Correlation: The Impact of Symmetry and Orbital Transformations

In this perspective, the various measures of electron correlation used in wave function theory, density functional theory and quantum information theory are briefly reviewed. We then focus on a more traditional metric based on dominant weights in the full configuration solution and discuss its behavior with respect to the choice of the N-electron and the one-electron basis. The impact of symmetry is discussed, and we emphasize that the distinction among determinants, configuration state functions and configurations as reference functions is useful because the latter incorporate spin-coupling into the reference and should thus reduce the complexity of the wave function expansion. The corresponding notions of single determinant, single spin-coupling and single configuration wave functions are discussed and the effect of orbital rotations on the multireference character is reviewed by analyzing a simple model system. In molecular systems, the extent of correlation effects should be limited by finite system size and in most cases the appropriate choices of one-electron and N-electron bases should be able to incorporate these into a low-complexity reference function, often a single configurational one.

TAO-DFT with the Polarizable Continuum Model

For the ground-state properties of gas-phase nanomolecules with multi-reference character, thermally assisted occupation (TAO) density functional theory (DFT) has recently been found to outperform the widely used Kohn-Sham DFT when traditional exchange-correlation energy functionals are employed. Aiming to explore solvation effects on the ground-state properties of nanomolecules with multi-reference character at a minimal computational cost, we combined TAO-DFT with the PCM (polarizable continuum model). In order to show its usefulness, TAO-DFT-based PCM (TAO-PCM) was used to predict the electronic properties of linear acenes in three different solvents (toluene, chlorobenzene, and water). According to TAO-PCM, in the presence of these solvents, the smaller acenes should have nonradical character, and the larger ones should have increasing polyradical character, revealing striking similarities to the past findings in the gas phase.

Quantum Many-Body Theory from a Solution of the N-Representability Problem

Here we present a many-body theory based on a solution of the N-representability problem in which the ground-state two-particle reduced density matrix (2-RDM) is determined directly without the many-particle wave function. We derive an equation that re-expresses physical constraints on higher-order RDMs to generate direct constraints on the 2-RDM, which are required for its derivation from an N-particle density matrix, known as N-representability conditions. The approach produces a complete hierarchy of 2-RDM constraints that do not depend explicitly upon the higher RDMs or the wave function. By using the two-particle part of a unitary decomposition of higher-order constraint matrices, we can solve the energy minimization by semidefinite programming in a form where the low-rank structure of these matrices can be potentially exploited. We illustrate by computing the ground-state electronic energy and properties of the H_{8} ring.

Electronic Properties of Hexagonal Graphene Quantum Rings from TAO-DFT

The reliable prediction of electronic properties associated with graphene nanosystems can be challenging for conventional electronic structure methods, such as Kohn-Sham (KS) density functional theory (DFT), due to the presence of strong static correlation effects in these systems. To address this challenge, TAO (thermally assisted occupation) DFT has been recently proposed. In the present study, we employ TAO-DFT to predict the electronic properties of n-HGQRs (i.e., the hexagonal graphene quantum rings consisting of n aromatic rings fused together at each side). From TAO-DFT, the ground states of n-HGQRs are singlets for all the cases investigated (n = 3-15). As the system size increases, there should be a transition from the nonradical to polyradical nature of ground-state n-HGQR. The latter should be intimately related to the localization of active TAO-orbitals at the inner and outer edges of n-HGQR, which increases with increasing system size.

Efficient Computation of Two-Electron Reduced Density Matrices via Selected Configuration Interaction

We create an approach to efficiently calculate two-electron reduced density matrices (2-RDMs) using selected configuration interaction wavefunctions. This is demonstrated using the specific example of Monte Carlo configuration interaction (MCCI). The computation of the 2-RDMs is accelerated by using ideas from fast implementations of full configuration interaction (FCI) and recent advances in implementing the Slater–Condon rules using hardware bitwise operations. This method enables a comparison of MCCI and truncated CI 2-RDMs with FCI values for a range of molecules, which includes stretched bonds and excited states. The accuracy in energies, wavefunctions, and 2-RDMs is seen to exhibit a similar behavior. We find that MCCI can reach sufficient accuracy of the 2-RDM using significantly fewer configurations than truncated CI, particularly for systems with strong multireference character.

Reducing the Quantum Many-Electron Problem to Two Electrons with Machine Learning

An outstanding challenge in chemical computation is the many-electron problem where computational methodologies scale prohibitively with system size. The energy of any molecule can be expressed as a weighted sum of the energies of two-electron wave functions that are computable from only a two-electron calculation. Despite the physical elegance of this extended "aufbau" principle, the determination of the distribution of weights─geminal occupations─for general molecular systems has remained elusive. Here we introduce a new paradigm for electronic structure where approximate geminal-occupation distributions are "learned" via a convolutional neural network. We show that the neural network learns the N-representability conditions, constraints on the distribution for it to represent an N-electron system. By training on hydrocarbon isomers with only 2-7 carbon atoms, we are able to predict the energies for isomers of octane as well as hydrocarbons with 8-15 carbons. The present work demonstrates that machine learning can be used to reduce the many-electron problem to an effective two-electron problem, opening new opportunities for accurately predicting electronic structure.

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.

Generation and Aerobic Oxidative Catalysis of a Cu(II) Superoxo Complex Supported by a Redox-Active Ligand

Cu systems feature prominently in aerobic oxidative catalysis in both biology and synthetic chemistry. Metal ligand cooperativity is a common theme in both areas as exemplified by galactose oxidase and by aminoxyl radicals in alcohol oxidations. This has motivated investigations into the aerobic chemistry of Cu and specifically the isolation and study of Cu-superoxo species that are invoked as key catalytic intermediates. While several examples of complexes that model biologically relevant Cu(II) superoxo intermediates have been reported, they are not typically competent aerobic catalysts. Here, we report a new Cu complex of the redox-active ligand ^tBu,TolDHP (2,5-bis((2-t-butylhydrazono)(p-tolyl)methyl)-pyrrole) that activates O₂ to generate a catalytically active Cu(II)-superoxo complex via ligand-based electron transfer. Characterization using ultraviolet (UV)-visible spectroscopy, Raman isotope labeling studies, and Cu extended X-ray absorption fine structure (EXAFS) analysis confirms the assignment of an end-on κ¹ superoxo complex. This Cu-O₂ complex engages in a range of aerobic catalytic oxidations with substrates including alcohols and aldehydes. These results demonstrate that bioinspired Cu systems can not only model important bioinorganic intermediates but can also mediate and provide mechanistic insight into aerobic oxidative transformations.

Bright, Modular, and Switchable Near-Infrared II Emission from Compact Tetrathiafulvalene-Based Diradicaloid Complexes

Near-infrared (NIR)-emitting molecules are promising candidates for biological sensing and imaging applications; however, many NIR dyes are large conjugated systems which frequently have issues with stability, solubility, and tunability. Here, we report a novel class of compact and tunable fluorescent diradicaloid complexes which are air-, water-, light-, and temperature-stable. These properties arise from a compressed π manifold which promotes an intense ligand-centered π-π transition in the NIR II (1000-1700 nm) region and which subsequently emits at ∼1200 nm. This emission is among the brightest known for monomolecular lumiphores with deep NIR II (>1100 nm) emission, nearly an order of magnitude brighter than the commercially available NIR II dye IR 26. Furthermore, this fluorescence is electrochemically sensitive, with efficient switching upon addition of redox agents. The brightness, stability, and modularity of this system distinguish it as a promising candidate for the development of new technologies built around NIR emission.

Interplay of Electronic and Geometric Structure Tunes Organic Biradical Character in Bimetallic Tetrathiafulvalene Tetrathiolate Complexes

The synthesis and design of organic biradicals with tunable singlet-triplet gaps have become the subject of significant research interest, owing to their possible photochemical applications and use in the development of molecular switches and conductors. Recently, tetrathiafulvalene tetrathiolate (TTFtt) has been demonstrated to exhibit such organic biradical character in doubly ionized bimetallic complexes. In this article we use high-level ab initio calculations to interrogate the electronic structure of a series of TTFtt-bridged metal complexes, resolving the factors governing their biradical character and singlet-triplet gaps. We show that the degree of biradical character correlates with a readily measured experimental predictor, the central TTFtt C-C bond length, and that it may be described by a one-parameter model, providing valuable insight for the future rational design of TTFtt based biradical compounds and materials.

Beginnings of exciton condensation in coronene analog of graphene double layer

Exciton condensation, a Bose-Einstein condensation of excitons into a single quantum state, has recently been achieved in low-dimensional materials including twin layers of graphene and van der Waals heterostructures. Here, we computationally examine the beginnings of exciton condensation in a double layer composed of coronene, a seven-benzene-ring patch of graphene. As a function of interlayer separation, we compute the exciton population in a single coherent quantum state, showing that the population peaks around 1.8 at distances near 2 Å. Visualization reveals interlayer excitons at the separation distance of the condensate. We determine the exciton population as a function of the twist angle between two coronene layers to reveal the magic angles at which the condensation peaks. As with previous recent calculations showing some exciton condensation in hexacene double layers and benzene stacks, the present two-electron reduced-density-matrix calculations with coronene provide computational evidence for the ability to realize exciton condensation in molecular-scale analogs of extended systems such as the graphene double layer.

iCISCF: An Iterative Configuration Interaction-Based Multiconfigurational Self-Consistent Field Theory for Large Active Spaces

An iterative configuration interaction (iCI)-based multiconfigurational self-consistent field (SCF) theory, iCISCF, is proposed to handle systems that require large active spaces. The success of iCISCF stems from three ingredients: (1) efficient selection of individual configuration state functions spanning the active space while maintaining full spin symmetry; (2) the use of Jacobi rotation for optimization of the active orbitals in conjunction with a quasi-Newton algorithm for the core/active-virtual and core-active orbital rotations; (3) a second-order perturbative treatment of the residual space left over by the selection procedure (i.e., iCISCF(2)). Several examples that go beyond the capability of CASSCF are taken as showcases to reveal the efficacy of iCISCF and iCISCF(2), facilitated by iCAS for imposed automatic selection and localization of active orbitals.

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.

The Static-Dynamic-Static Family of Methods for Strongly Correlated Electrons: Methodology and Benchmarking

A series of methods (SDSCI, SDSPT2, iCI, iCIPT2, iCISCF(2), iVI, and iCAS) is introduced to accurately describe strongly correlated systems of electrons. Born from the (restricted) static-dynamic-static (SDS) framework for designing many-electron wave functions, SDSCI is a minimal multireference (MR) configuration interaction (CI) approach that constructs and diagonalizes a [Formula: see text] matrix for [Formula: see text] states, regardless of the numbers of orbitals and electrons to be correlated. If the full molecular Hamiltonian H in the QHQ block (which describes couplings between functions of the first-order interaction space Q) of the SDSCI CI matrix is replaced with a zeroth-order Hamiltonian [Formula: see text] before the diagonalization is taken, we obtain SDSPT2, a CI-like second-order perturbation theory (PT2). Unlike most variants of MRPT2, SDSPT2 treats single and multiple states in the same way and is particularly advantageous in the presence of near degeneracy. On the other hand, if the SDSCI procedure is repeated until convergence, we will have iterative CI (iCI), which can converge quickly from the above to the exact solutions (full CI) even when starting with a poor guess. When further combined with the selection of important configurations followed by a PT2 treatment of dynamic correlation, iCI becomes iCIPT2, which is a near-exact theory for medium-sized systems. The microiterations of iCI for relaxing the coefficients of contracted many-electron functions can be generalized to an iterative vector interaction (iVI) approach for finding exterior or interior roots of a given matrix, in which the dimension of the search subspace is fixed by either the number of target roots or the user-specified energy window. Naturally, iCIPT2 can be employed as the active space solver of the complete active space (CAS) self-consistent field, leading to iCISCF(2), which can further be combined with iCAS for automated selection of active orbitals and assurance of the same CAS for all states and all geometries. The methods are calibrated by taking the Thiel set of benchmark systems as examples. Results for the corresponding cations, a new set of benchmark systems, are also reported.

Exciton Condensation in Molecular-Scale van der Waals Stacks

Recent experiments have realized the Bose-Einstein condensation of excitons, known as exciton condensation, in extended systems such as bilayer graphene and van der Waals heterostructures. Here we computationally demonstrate the beginnings of exciton condensation in multilayer, molecular-scale van der Waals stacks composed of benzene subunits. The populations of excitons, which are computed from the largest eigenvalue of the particle-hole reduced density matrix (RDM) through advanced variational RDM calculations, are shown to increase with the length of the stack. The large eigenvalue indicates a nonclassical long-range ordering of the excitons that can support the frictionless flow of energy. Moreover, we use chemical substitutions and geometric modifications to tune the extent of the condensation. Results suggest exciton condensation in a potentially large family of molecular systems with applications to energy-efficient transport.

Electronic Properties of Carbon Nanobelts Predicted by Thermally-Assisted-Occupation DFT

Accurate prediction of properties of large-scale multi-reference (MR) electronic systems remains difficult for traditional computational methods (e.g., the Hartree-Fock theory and Kohn-Sham density functional theory (DFT)). Recently, thermally-assisted-occupation (TAO)-DFT has been demonstrated to offer reliable description of electronic properties of various large-scale MR electronic systems. Consequently, in this work, TAO-DFT is used to unlock the electronic properties associated with C-Belt[n] (i.e., the carbon nanobelts containing n fused 12-membered carbon rings). Our calculations show that for all the system sizes reported (n = 4-24), C-Belt[n] have singlet ground states. In general, the larger the size of C-Belt[n], the more pronounced the MR character of ground-state C-Belt[n], as evident from the symmetrized von Neumann entropy and the occupation numbers of active TAO-orbitals. Furthermore, the active TAO-orbitals are delocalized along the circumference of C-Belt[n], as evident from the visualization of active TAO-orbitals.

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.

iCAS: Imposed Automatic Selection and Localization of Complete Active Spaces

It is shown that in the spirit of "from fragments to molecule" for localizing molecular orbitals [J. Chem. Theory Comput. 2011, 7, 3643], a prechosen set of occupied/virtual valence/core atomic/fragmental orbitals can be transformed to an equivalent set of localized occupied/virtual pre-localized molecular orbitals (pre-LMO), which can then be taken as probes to select the same number of maximally matching localized occupied/virtual Hartree-Fock (HF) or restricted open-shell HF (ROHF) molecular orbitals as the initial local orbitals spanning the desired complete active space (CAS). In each cycle of the self-consistent field (SCF) calculation, the CASSCF orbitals can be localized by means of the noniterative "top-down least-change" algorithm for localizing ROHF orbitals [J. Chem. Phys. 2017, 146, 104104] such that the maximum matching between the orbitals of two adjacent iterations can readily be monitored, leading finally to converged localized CASSCF orbitals that overlap most the guess orbitals. Such an approach is to be dubbed as "imposed CASSCF" (iCASSCF or simply iCAS in short) for good reasons: (1) it has been assumed that only those electronic states that have largest projections onto the active space defined by the prechosen atomic/fragmental orbitals are to be targeted. This is certainly an imposed constraint but has wide applications in organic and transition metal chemistry where valence (or core) atomic/fragmental orbitals can readily be identified. (2) The selection of both initial and optimized local active orbitals is imposed from the very beginning by the pre-LMOs (which span the same space as the prechosen atomic/fragmental orbitals). Apart from the (imposed) automation and localization, iCAS has two additional merits: (a) the guess orbitals are guaranteed to be the same for all geometries, for the pre-LMOs do not change in character with geometry and (b) the use of localized orbitals facilitates the SCF convergence, particularly for large active spaces. Both organic molecules and transition-metal complexes are taken as showcases to reveal the efficacy of iCAS.

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.

Correlation-driven phenomena in periodic molecular systems from variational two-electron reduced density matrix theory

Correlation-driven phenomena in molecular periodic systems are challenging to predict computationally not only because such systems are periodically infinite but also because they are typically strongly correlated. Here, we generalize the variational two-electron reduced density matrix (2-RDM) theory to compute the energies and properties of strongly correlated periodic systems. The 2-RDM of the unit cell is directly computed subject to necessary N-representability conditions such that the unit-cell 2-RDM represents at least one N-electron density matrix. Two canonical but non-trivial systems, periodic metallic hydrogen chains and periodic acenes, are treated to demonstrate the methodology. We show that while single-reference correlation theories do not capture the strong (static) correlation effects in either of these molecular systems, the periodic variational 2-RDM theory predicts the Mott metal-to-insulator transition in the hydrogen chains and the length-dependent polyradical formation in acenes. For both hydrogen chains and acenes, the periodic calculations are compared with previous non-periodic calculations with the results showing a significant change in energies and increase in the electron correlation from the periodic boundary conditions. The 2-RDM theory, which allows for much larger active spaces than are traditionally possible, is applicable to studying correlation-driven phenomena in general periodic molecular solids and materials.

Localized Active Space Pair-Density Functional Theory

Accurate quantum chemical methods for the prediction of spin-state energy gaps for strongly correlated systems are computationally expensive and scale poorly with the size of the system. This makes calculations for many experimentally interesting molecules impractical even with abundant computational resources. Previous work has shown that the localized active space (LAS) self-consistent field (SCF) method can be an efficient way to obtain multiconfiguration SCF wave functions of comparable quality to the corresponding complete active space (CAS) ones. To obtain quantitative results, a post-SCF method is needed to estimate the complete correlation energy. One such method is multiconfiguration pair-density functional theory (PDFT), which calculates the energy based on the density and on-top pair density obtained from a multiconfiguration wave function. In this work, we introduce localized-active-space PDFT, which uses a LAS wave function for subsequent PDFT calculations. The method is tested by computing spin-state energies and gaps in conjugated organic molecules and a bimetallic compound and comparing to the corresponding CAS-PDFT values.

Spin-Flip Pair-Density Functional Theory: A Practical Approach To Treat Static and Dynamical Correlations in Large Molecules

We present a practical approach to treat static and dynamical correlation accurately in large multiconfigurational systems. The static correlation is taken into account by using the spin-flip approach, which is well-known for capturing static correlation accurately at low-computational expense. Unlike previous approaches to add dynamical correlation to spin-flip models which use perturbation theory or coupled-cluster theory, we explore the ability to use the on-top pair-density functional theory approaches recently developed by Gagliardi and co-workers (J. Comput. Theor. Chem., 2014, 10, 3669). External relaxations are performed in the spin-flip calculations through a restricted active space framework for which a truncation scheme for the orbitals used in the external excitation is presented. The performance of the approach is demonstrated by computing energy gaps between ground and excited states for diradicals, triradicals, and linear polyacene chains ranging from naphthalene to dodecacene. Accurate results are obtained using the new approach for these challenging open-shell molecular systems.

Energetics of Formation of Cyclacenes from 2,3-Didehydroacenes and Implications for Astrochemistry

The carriers of the diffuse interstellar bands (DIBs) are still largely unknown although polycyclic aromatic hydrocarbons, carbon chains, and fullerenes are likely candidates. A recent analysis of the properties of n-acenes of general formula C4n+2 H2n+4 suggested that these could be potential carriers of some DIBs. Dehydrogenation reactions of n-acenes after absorption of an interstellar UV photon may result in dehydroacenes. Here the reaction energies and barriers for formation of n-cyclacenes from 2,3-didehydroacenes (n-DDA) by intramolecular Diels-Alder reaction to dihydro-etheno-cyclacenes (n-DEC) followed by ejection of ethyne by retro-Diels-Alder reactions are analyzed using thermally assisted occupation density functional theory (TAO-DFT) for n=10-20. It is found that the barriers for each of the steps depend on the ring strain of the underlying n-cyclacene, and that the ring strain of n-DEC is about 75 % of that of the corresponding n-cyclacene. In each case, ethyne extrusion is the step with the highest energy barrier, but these barriers are smaller than CH bond dissociation energies, suggesting that formation of cyclacenes is an energetically conceivable fate of n-acenes after multiple absorption of UV photons.

Quantum Information and Algorithms for Correlated Quantum Matter

Discoveries in quantum materials, which are characterized by the strongly quantum-mechanical nature of electrons and atoms, have revealed exotic properties that arise from correlations. It is the promise of quantum materials for quantum information science superimposed with the potential of new computational quantum algorithms to discover new quantum materials that inspires this Review. We anticipate that quantum materials to be discovered and developed in the next years will transform the areas of quantum information processing including communication, storage, and computing. Simultaneously, efforts toward developing new quantum algorithmic approaches for quantum simulation and advanced calculation methods for many-body quantum systems enable major advances toward functional quantum materials and their deployment. The advent of quantum computing brings new possibilities for eliminating the exponential complexity that has stymied simulation of correlated quantum systems on high-performance classical computers. Here, we review new algorithms and computational approaches to predict and understand the behavior of correlated quantum matter. The strongly interdisciplinary nature of the topics covered necessitates a common language to integrate ideas from these fields. We aim to provide this common language while weaving together fields across electronic structure theory, quantum electrodynamics, algorithm design, and open quantum systems. Our Review is timely in presenting the state-of-the-art in the field toward algorithms with nonexponential complexity for correlated quantum matter with applications in grand-challenge problems. Looking to the future, at the intersection of quantum information science and algorithms for correlated quantum matter, we envision seminal advances in predicting many-body quantum states and describing excitonic quantum matter and large-scale entangled states, a better understanding of high-temperature superconductivity, and quantifying open quantum system dynamics.

Portably parallel construction of a configuration-interaction wave function from a matrix-product state using the Charm++ framework

The construction of configuration-interaction (CI) expansions from a matrix product state (MPS) involves numerous matrix operations and the skillful sampling of important configurations in a large Hilbert space. In this work, we present an efficient procedure for constructing CI expansions from MPS employing the parallel object-oriented Charm++ programming framework, upon which automatic load-balancing and object migrating facilities can be employed. This procedure was employed in the MPS-to-CI utility (Moritz et al., J. Chem. Phys. 2007, 126, 224109), the sampling-reconstructed complete active-space algorithm (SR-CAS, Boguslawski et al., J. Chem. Phys. 2011, 134, 224101), and the entanglement-driven genetic algorithm (EDGA, Luo et al., J. Chem. Theory Comput. 2017, 13, 4699). It enhances productivity and allows the sampling programs to evolve to their population-expansion versions, for example, EDGA with population expansion (PE-EDGA). Further, examples of 1,2-dioxetanone and firefly dioxetanone anion (FDO^- ) molecules demonstrated the following: (a) parallel efficiencies can be persistently improved by simply by increasing the proportions of the asynchronous executions and (b) a sampled CAS-type CI wave function of a bi-radical-state FDO^- molecule utilizing the full valence (30e,26o) active space can be constructed within a few hours with using thousands of cores.

A New Mixing of Nonlocal Exchange and Nonlocal Correlation with Multiconfiguration Pair-Density Functional Theory

We propose a hybrid multiconfiguration pair-density functional theory (HMC-PDFT) that is a weighted average of complete-active-space self-consistent-field (CASSCF) and multiconfiguration pair-density functional theory (MC-PDFT) energies with a semiempirical parameter to control the fraction of CASSCF energy. We also explore a more general two-parameter hybrid method with a scaled correlation energy that allows us to compare to the recently proposed λ-MC-PDFT method. We scan the parameter space for the scaled-correlation method using test sets consisting of electronic excitation energies and diatomic bond energies, and we find no significant improvement by introducing the scaling parameter. We find that unscaled HMC-PDFT offers significantly improved accuracy over both CASSCF and the original MC-PDFT for a wide range of systems, and we present as an example of this approach "tPBE0", the "translated" MC-PDFT generalization of the popular PBE0 hybrid Kohn-Sham density functional.

Efficient local energy evaluation for multi-Slater wave functions in orbital space quantum Monte Carlo

We present an algorithm for calculating the local energy of a multi-Slater wave function in orbital space quantum Monte Carlo (QMC). Recent developments in selected configuration interaction methods have led to increased interest in using multi-Slater trial wave functions in various QMC methods. For an ab initio Hamiltonian, our algorithm has a cost scaling of O(n⁵ + n_c), as opposed to the O(n⁴n_c) scaling of existing orbital space algorithms, where n is the system size and n_c is the number of configurations in the wave function. We present our method using variational Monte Carlo calculations with the Jastrow multi-Slater wave function, although the formalism should be applicable for auxiliary field QMC. We apply it to polyacetylene and demonstrate the possibility of using a much larger number of configurations than possible using existing methods.

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].

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.

Entangled Electrons Drive a Non-superexchange Mechanism in a Cobalt Quinoid Dimer Complex

A central theme in chemistry is the understanding of the mechanisms that drive chemical transformations. A well-known, highly cited mechanism in organometallic chemistry is the superexchange mechanism in which unpaired electrons on two or more metal centers interact through an electron pair of the bridging ligand. We use a combination of novel synthesis and computation to show that such interactions may in fact occur by a more direct mechanism than superexchange that is based on direct quantum entanglement of the two metal centers. Specifically, we synthesize and experimentally characterize a novel cobalt dimer complex with benzoquinoid bridging ligands and investigate its electronic structure with the variational two-electron reduced density matrix method using large active spaces. The result draws novel connections between inorganic mechanisms and quantum entanglement, thereby opening new possibilities for the design of strongly correlated organometallic compounds whose magnetic and spin properties have applications in superconductors, energy storage, thermoelectrics, and spintronics.

Psi4 1.4: Open-source software for high-throughput quantum chemistry

PSI4 is a free and open-source ab initio electronic structure program providing implementations of Hartree-Fock, density functional theory, many-body perturbation theory, configuration interaction, density cumulant theory, symmetry-adapted perturbation theory, and coupled-cluster theory. Most of the methods are quite efficient, thanks to density fitting and multi-core parallelism. The program is a hybrid of C++ and Python, and calculations may be run with very simple text files or using the Python API, facilitating post-processing and complex workflows; method developers also have access to most of PSI4's core functionalities via Python. Job specification may be passed using The Molecular Sciences Software Institute (MolSSI) QCSCHEMA data format, facilitating interoperability. A rewrite of our top-level computation driver, and concomitant adoption of the MolSSI QCARCHIVE INFRASTRUCTURE project, makes the latest version of PSI4 well suited to distributed computation of large numbers of independent tasks. The project has fostered the development of independent software components that may be reused in other quantum chemistry programs.

Global Hybrid Multiconfiguration Pair-Density Functional Theory

A global hybrid extension of multiconfiguration pair-density functional theory (MC-PDFT) is developed. Using a linear decomposition of the electron-electron repulsion term, a fraction λ of the nonlocal exchange interaction, obtained from variational two-electron reduced-density matrix (v2RDM)-driven complete active-space self-consistent field (CASSCF) theory, is combined with its local counterpart, obtained from an on-top pair-density functional. The resulting scheme (called λ-MC-PDFT) inherits the benefits of MC-PDFT (e.g., its simplicity and the resolution of the symmetry dilemma) and, when combined with the v2RDM approach to CASSCF, requires only polynomially scaling computational effort. As a result, λ-MC-PDFT can efficiently describe static and dynamical correlation effects in strongly correlated systems. The efficacy of the approach is assessed for several challenging multiconfigurational problems, including the dissociation of molecular nitrogen, the double dissociation of a water molecule, and the 1,3-dipolar cycloadditions of ozone to ethylene and ozone to acetylene in the O3ADD6 benchmark set.

Elucidating the multi-configurational character of the firefly dioxetanone anion and its prototypes in the biradical region using full valence active spaces

We analyzed the near-degenerate states of the firefly dioxetanone anion (FDO-) and its prototypes, especially in the biradical region, using multi-configurational approaches. The importance of utilizing full valence active spaces by means of density-matrix renormalization group self-consistent field (DMRG-SCF) calculations was described. Our results revealed that the neglect of some valence orbitals can affect the quantitative accuracy in later multi-reference calculations or the qualitative conclusion when optimizing conical intersections. Using all of the relevant valence orbitals of FDO-, we confirmed that there were two conical intersections, as reported in previous work, and that the intersecting states were changed when the active space was enlarged. Beyond these, we found that there were strong interactions between states in the biradical regions, in which the changes in entanglements can be used to visualize the interacting state evolution.

Non-equilibrium steady state conductivity in cyclo[18]carbon and its boron nitride analogue

A ring-shaped carbon allotrope was recently synthesized for the first time, reinvigorating theoretical interest in this class of molecules.

An adiabatic connection for doubly-occupied configuration interaction wave functions

An adiabatic connection (AC) is developed as an electron correlation correction for doubly occupied configuration interaction (DOCI) wave functions. Following the work of Pernal [Phys. Rev. Lett. 120, 013001 (2018)], the working equations of the approach, termed AC-DOCI, are rooted in the extended random phase approximation (ERPA) and require knowledge of only the ground-state two-electron reduced density matrix (2RDM) from the DOCI. As such, the AC is naturally suited to modeling electron correlation in variational 2RDM (v2RDM)-based approximations to the DOCI. The v2RDM-driven AC-DOCI is applied to the dissociation of molecular nitrogen and the double dissociation of water; the approach yields energies that are similar in quality to those from second-order multireference perturbation theory near equilibrium, but the quality of the AC-DOCI energy degrades at stretched geometries. The exact adiabatic connection path suggests the assumption that the one-electron reduced-density matrix is constant along the AC path is invalid at stretched geometries, but this deficiency alone cannot explain the observed behavior. Rather, it appears that the ERPA's single-particle-transition ansatz cannot, in general, provide good approximations to the 2RDM along the AC path. The AC-DOCI is also applied to a set of 45 reaction energies; for these systems, the approach has an average accuracy that is comparable to that of single-reference second-order many-body perturbation theory.