Active Orbital Space Research Articles

Reduced Scaling Correlated Natural Transition Orbitals for Multilevel Coupled Cluster Calculations.

Multilevel coupled cluster theory offers reduced scaling computation of intensive properties in systems that are too large for standard coupled cluster calculations. A significant benefit of the multilevel coupled cluster framework is the possibility of calculating intensive properties that are not tightly localized if an appropriate set of active orbitals is used. Correlated natural transition orbitals (CNTOs) are tailored to describe excitation processes. For multilevel coupled cluster singles and doubles (MLCCSD) and singles and perturbative doubles (MLCC2) calculations, the construction of CNTOs generally becomes the computational bottleneck. Here, we demonstrate how CNTOs can be obtained with operations, eliminating the -scaling steps involved in the original approach. This reduction in scaling moves the bottleneck of MLCC2 and MLCCSD calculations from the active orbital space preparation to the MLCC2 and MLCCSD equations with -scaling.

Optical Properties of Defects in Solids via Quantum Embedding with Good Active Space Orbitals

Optical Properties of Defects in Solids via Quantum Embedding with Good Active Space Orbitals

The trials and triumphs of modelling X-ray absorption spectra of transition metal phthalocyanines.

This study explores the electronic structure of Co, Fe and Mn phthalocyanines (TMPcs) as well as their perfluorinated counterparts through a series of electronic structure calculations utilizing multireference methods and by simulating their metal L-edge and ligand (nitrogen and fluorine) K-edge X-ray absorption spectra (XAS) in an angle-resolved manner. Simulations targeting different ground-state symmetries, where relevant, have been conducted to observe changes in the N K-edge lineshape. The applicability of the quasi-degenerate formulation of n-electron valence state perturbation theory (QD-NEVPT2) for L-edge X-ray absorption spectroscopy (XAS) is evaluated, alongside the use of a restricted active space (RAS) formalism to describe the final-state multiplets generated by L-shell X-ray processes. Our findings provide valuable insights into the electronic properties of TMPcs, in particular with respect to the effect of fluorination, and demonstrate the broad applicability of various formulations of NEVPT2 in spectral simulations. Moreover, this study highlights the utility of manual truncation of the configuration spaces in order to allow for large active orbital spaces in aforementioned calculations.

PyBEST: Improved functionality and enhanced performance

PyBEST: Improved functionality and enhanced performance

Automated, Consistent, and Even-Handed Selection of Active Orbital Spaces for Quantum Embedding.

A widely used strategy to reduce the computational cost of quantum-chemical calculations is to partition the system into an active subsystem, which is the focus of the computational efforts, and an environment that is treated at a lower computational level. The system partitioning is mostly based on localized molecular orbitals. When reaction paths or energy differences are to be calculated, it is crucial to keep the orbital space consistent for all structures. Inconsistencies in orbital space can lead to unpredictable errors on the potential energy surface. While successful strategies to ensure this consistency have been established for organic and even metal-organic systems, these methods often fail for metal clusters or nanoparticles with a high density of near-degenerate and delocalized molecular orbitals. However, such systems are highly relevant for catalysis. Accurate yet feasible quantum-mechanical ab initio calculations are therefore highly desired. In this work, we present an approach based on the subsystem projected atomic orbital decomposition algorithm that allows us to ensure automated and consistent partitioning even for systems with delocalized and near-degenerate molecular orbitals and demonstrate the validity of this method for the binding energies of small molecules on transition-metal clusters.

Generalization of the Tensor Product Selected CI Method for Molecular Excited States.

In a recent paper [JCTC, 2020, 16, 6098], we introduced a new approach for accurately approximating full CI ground states in large electronic active-spaces called Tensor Product Selected CI (TPSCI). In TPSCI, a large orbital active space is first partitioned into disjoint sets (clusters) for which the exact, local many-body eigenstates are obtained. Tensor products of these locally correlated many-body states are taken as the basis for the full, global Hilbert space. By folding correlation into the basis states themselves, the low-energy eigenstates become increasingly sparse, creating a more compact selected CI expansion. While we demonstrated that this approach can improve accuracy for a variety of systems, there is even greater potential for applications to excited states, particularly those which have some excited-state character. In this paper, we report on the accuracy of TPSCI for excited states, including a far more efficient implementation in the Julia programming language. In traditional SCI methods that use a Slater determinant basis, accurate excitation energies are obtained only after a linear extrapolation and at a large computational cost. We find that TPSCI with perturbative corrections provides accurate excitation energies for several excited states of various polycyclic aromatic hydrocarbons with respect to the extrapolated result (i.e., near exact result). Further, we use TPSCI to report highly accurate estimates of the lowest 31 eigenstates for a tetracene tetramer system with an active space of 40 electrons in 40 orbitals, giving direct access to the initial bright states and the resulting 18 doubly excited (biexcitonic) states.

Complete Active Space Methods for NISQ Devices: The Importance of Canonical Orbital Optimization for Accuracy and Noise Resilience.

To avoid the scaling of the number of qubits with the size of the basis set, one can divide the molecular space into active and inactive regions, which is also known as complete active space methods. However, selecting the active space alone is not enough to accurately describe quantum mechanical effects such as correlation. This study emphasizes the importance of optimizing the active space orbitals to describe correlation and improve the basis-dependent Hartree-Fock energies. We will explore classical and quantum computation methods for orbital optimization and compare the chemically inspired ansatz, UCCSD, with the classical full CI approach for describing the active space in both weakly and strongly correlated molecules. Finally, we will investigate the practical implementation of a quantum CASSCF, where hardware-efficient circuits must be used and noise can interfere with accuracy and convergence. Additionally, we will examine the impact of using canonical and noncanonical active orbitals on the convergence of the quantum CASSCF routine in the presence of noise.

Corresponding Active Orbital Spaces along Chemical Reaction Paths.

The accuracy of reaction energy profiles calculated with multiconfigurational electronic structure methods and corrected by multireference perturbation theory depends crucially on consistent active orbital spaces selected along the reaction path. However, it has been challenging to choose molecular orbitals that can be considered corresponding in different molecular structures. Here, we demonstrate how active orbital spaces can be selected consistently along reaction coordinates in a fully automatized way. The approach requires no structure interpolation between reactants and products. Instead, it emerges from a synergy of the Direct Orbital Selection orbital mapping ansatz combined with our fully automated active space selection algorithm autoCAS. We demonstrate our algorithm for the potential energy profile of the homolytic carbon-carbon bond dissociation and rotation around the double bond of 1-pentene in the electronic ground state. However, our algorithm also applies to electronically excited Born-Oppenheimer surfaces.

Selected configuration interaction wave functions in phaseless auxiliary field quantum Monte Carlo

We present efficient algorithms for using selected configuration interaction (sCI) trial wave functions in phaseless auxiliary field quantum Monte Carlo (ph-AFQMC). These advances, geared toward optimizing computational performance for longer configuration interaction expansions, allow us to use up to a million configurations in the trial state for ph-AFQMC. In one example, we found the cost of ph-AFQMC per sample to increase only by a factor of about 3 for a calculation with 104 configurations compared to that with a single one, demonstrating the tiny computational overhead due to a longer expansion. This favorable scaling allows us to study the systematic convergence of the phaseless bias in auxiliary field quantum Monte Carlo calculations with an increasing number of configurations and provides a means to gauge the accuracy of ph-AFQMC with other trial states. We also show how the scalability issues of sCI trial states for large system sizes could be mitigated by restricting them to a moderately sized orbital active space and leveraging the near-cancellation of out of active space phaseless errors.

Machine Learning-Assisted Selection of Active Spaces for Strongly Correlated Transition Metal Systems.

Active space quantum chemical methods could provide very accurate description of strongly correlated electronic systems, which is of tremendous value for natural sciences. The proper choice of the active space is crucial but a nontrivial task. In this article, we present a neural network-based approach for automatic selection of active spaces, focused on transition metal systems. The training set has been formed from artificial systems composed of one transition metal and various ligands, on which we have performed the density matrix renormalization group and calculated the single-site entropy. On the selected set of systems, ranging from small benchmark molecules up to larger challenging systems involving two metallic centers, we demonstrate that our machine learning models could predict the active space orbitals with reasonable accuracy. We also tested the transferability on out-of-the-model systems, including bimetallic complexes and complexes with ligands, which were not involved in the training set. Also, we tested the correctness of the automatically selected active spaces on a Fe(II)-porphyrin model, where we studied the lowest states at the DMRG level and compared the energy difference between spin states or the energy difference between conformations of ferrocene with recent studies.

Active orbital preservation for multiconfigurational self-consistent field.

We introduce Active Orbital Preservation for Multiconfigurational Self-Consistent Field (AOP-MCSCF), an automated approach to improving the consistency of active space orbitals over multiple molecular configurations. Our approach is based on maximum overlap with a reference set of active space orbitals taken from a single geometry of a chromophore in the gas phase and can be used to automatically preserve the appropriate orbitals of the chromophore across multiple thermally sampled configurations, even when the chromophore is solvated by quantum-mechanically treated water molecules. In particular, using the singular value decomposition of a Molecular Orbital (MO) overlap matrix between the system and reference, we rotate the MOs of the system to align with the reference active space orbitals and use the resulting rotated orbitals as an initial guess to a MCSCF calculation. We demonstrate the approach on aqueous p-hydroxybenzylidene-imidazolinone (HBI) and find that AOP-MCSCF converges to the "correct" orbitals for over 90% of 3000 thermally sampled configurations. In addition, we compute the linear absorption spectrum and find excellent agreement with new experimental measurements up to 5.4eV (230nm). We show that electrostatic contributions to the solvation energy of HBI largely explain the observed state-dependent solvatochromism.

Computational Investigations of the Lithium Superoxide Dimer Rearrangement on Noisy Quantum Devices.

Quantum chemistry studies of biradical systems are challenging due to the required multiconfigurational nature of the wavefunction. In this work, Variational Quantum Eigensolver (VQE) is used to compute the energy profile for the lithium superoxide dimer rearrangement, involving biradical species, on quantum simulators and devices. Considering that current quantum devices can only handle limited number of qubits, we present guidelines for selecting an appropriate active space to perform computations on chemical systems that require many qubits. We show that with VQE performed with a quantum simulator reproduces results obtained with full-configuration interaction (Full CI) for the chosen active space. However, results deviate from exact values by about 39 mHa for calculations on a quantum device. This deviation can be improved to about 4 mHa using the readout mitigation approach and can be further improved to 2 mHa, approaching chemical accuracy, using the state tomography technique to purify the calculated quantum state.

Tailored coupled cluster theory in varying correlation regimes.

The tailored coupled cluster (TCC) approach is a promising ansatz that preserves the simplicity of single-reference coupled cluster theory while incorporating a multi-reference wave function through amplitudes obtained from a preceding multi-configurational calculation. Here, we present a detailed analysis of the TCC wave function based on model systems, which require an accurate description of both static and dynamic correlation. We investigate the reliability of the TCC approach with respect to the exact wave function. In addition to the error in the electronic energy and standard coupled cluster diagnostics, we exploit the overlap of TCC and full configuration interaction wave functions as a quality measure. We critically review issues, such as the required size of the active space, size-consistency, symmetry breaking in the wave function, and the dependence of TCC on the reference wave function. We observe that possible errors caused by symmetry breaking can be mitigated by employing the determinant with the largest weight in the active space as reference for the TCC calculation. We find the TCC model to be promising in calculations with active orbital spaces which include all orbitals with a large single-orbital entropy, even if the active spaces become very large and then may require modern active-space approaches that are not restricted to comparatively small numbers of orbitals. Furthermore, utilizing large active spaces can improve on the TCC wave function approximation and reduce the size-consistency error because the presence of highly excited determinants affects the accuracy of the coefficients of low-excited determinants in the active space.

Comparison of Methods for Active Orbital Selection in Multiconfigurational Calculations.

Several methods of constructing the active orbital space for multiconfigurational wave functions are compared on typical moderately strongly or strongly correlated ground-state molecules. The relative merits of these methods and problems inherent in multiconfigurational calculations are discussed. Strong correlation in the ground electronic state is found typically in larger conjugated and in antiaromatic systems, transition states which involve bond breaking or formation, and transition metal complexes. Our examples include polyenes, polyacenes, the reactant, product and transition state of the Bergman cyclization, and two transition metal complexes: Hieber's anion [(CO)3FeNO]- and ferrocene. For the systems investigated, the simplest and oldest selection method, based on the fractional occupancy of unrestricted Hartree-Fock natural orbitals (the UNO criterion), yields the same active space as much more expensive approximate full CI methods. A disadvantage of this method used to be the difficulty of finding broken spin symmetry UHF solutions. However, our analytical method, accurate to fourth order in the orbital rotation angles (Tóth and Pulay J. Chem. Phys. 2016, 145, 164102.), has solved this problem. Two further advantages of the UNO criterion are that, unlike most other methods, it measures not only the energetic proximity to the Fermi level but also the magnitude of the exchange interaction with strongly occupied orbitals and therefore allows the estimation of the correlation strength for orbital selection in Restricted Active Space methods.

Equation-of-Motion MLCCSD and CCSD-in-HF Oscillator Strengths and Their Application to Core Excitations.

We present an implementation of equation-of-motion oscillator strengths for the multilevel CCSD (MLCCSD) model where CCS is used as the lower level method (CCS/CCSD). In this model, the double excitations of the cluster operator are restricted to an active orbital space, whereas the single excitations are unrestricted. Calculated nitrogen K-edge spectra of adenosine, adenosine triphosphate (ATP), and an ATP-water system are used to demonstrate the performance of the model. Projected atomic orbitals (PAOs) are used to partition the virtual space into active and inactive orbital sets. Cholesky decomposition of the Hartree–Fock density is used to partition the occupied orbitals. This Cholesky-PAO partitioning is cheap, scaling as , and is suitable for the calculation of core excitations, which are localized in character. By restricting the single excitations of the cluster operator to the active space, as well as the double excitations, the CCSD-in-HF model is obtained. A comparison of the two models—MLCCSD and CCSD-in-HF—is presented for the core excitation spectra of the adenosine and ATP systems.

Selected Configuration Interaction in a Basis of Cluster State Tensor Products.

Selected configuration interaction (SCI) methods are currently enjoying a resurgence due to several recent developments which improve either the overall computational efficiency or the compactness of the resulting SCI vector. These recent advances have made it possible to get full CI (FCI) quality results for much larger orbital active spaces compared to conventional approaches. However, due to the starting assumption that the FCI vector has only a small number of significant Slater determinants, SCI becomes intractable for systems with strong correlation. This paper introduces a method for developing SCI algorithms in a way which exploits local molecular structure to significantly reduce the number of SCI variables. The proposed method is defined by first grouping the orbitals into clusters over which we can define many-particle cluster states. We then directly perform the SCI algorithm in a basis of tensor products of cluster states instead of Slater determinants. While the approach is general for arbitrarily defined cluster states, we find significantly improved performance by defining cluster states through a Tucker decomposition of the global (and sparse) SCI vector. To demonstrate the potential of this method, called tensor product selected configuration interaction (TPSCI), we present numerical results for a diverse set of examples: (1) modified Hubbard model with different inter- and intracluster hopping terms, (2) less obviously clusterable cases of bond breaking in N2 and F2, and (3) ground state energies of large planar π-conjugated systems with active spaces of up to 42 electrons in 42 orbitals. These numerical results show that TPSCI can be used to significantly reduce the number of SCI variables in the variational space, thus paving a path for extending these deterministic and variational SCI approaches to a wider range of physical systems.

Time-dependent optimized coupled-cluster method for multielectron dynamics. II. A coupled electron-pair approximation.

We report the implementation of a cost-effective approximation method within the framework of the time-dependent optimized coupled-cluster (TD-OCC) method [T. Sato et al., J. Chem. Phys. 148, 051101 (2018)] for real-time simulations of intense laser-driven multielectron dynamics. The method, designated as TD-OCEPA0, is a time-dependent extension of the simplest version of the coupled-electron pair approximation with optimized orbitals [U. Bozkaya and C. D. Sherrill, J. Chem. Phys. 139, 054104 (2013)]. It is size extensive, gauge invariant, and computationally much more efficient than the TD-OCC method with double excitations. We employed this method to simulate the electron dynamics in Ne and Ar atoms exposed to intense near infrared laser pulses with various intensities. The computed results, including high-harmonic generation spectra and ionization yields, are compared with those of various other methods ranging from uncorrelated time-dependent Hartree-Fock to fully correlated (within the active orbital space) time-dependent complete-active-space self-consistent field (TD-CASSCF). The TD-OCEPA0 results show good agreement with TD-CASSCF ones for moderate laser intensities. For higher intensities, however, TD-OCEPA0 tends to overestimate the correlation effect, as occasionally observed for CEPA0 in the ground-state correlation energy calculations.

Multilevel CC2 and CCSD Methods with Correlated Natural Transition Orbitals.

In the multilevel coupled cluster approach, an active orbital space is treated at a higher level of coupled cluster theory than the remaining inactive orbitals. We introduce the multilevel CC2 method where CC2 is used for the active orbital space. Furthermore, we present a simplified formulation of the multilevel CCSD method where CCSD is used for the active space. The simplification lies in the evaluation of the CC2 amplitudes in the inactive space; these CC2 amplitudes have previously been determined iteratively. We use correlated natural transition orbitals to determine the active orbital spaces. The convergence of the multilevel CC2 and multilevel CCSD valence excitation energies is established with proof-of-concept calculations. The methods are also applied to two larger systems: p-nitroaniline in water and amoxicillin. The calculations on the p-nitroaniline-water system illustrate the usefulness of multilevel coupled cluster methods for molecules in solution and for charge transfer excitations.

Effects of core space and excitation levels on ground-state correlation and photoionization dynamics of Be and Ne.

We explore the effects of correlation on the ground-state energies and on photoionization dynamics in atomic Be and Ne. We apply the time-dependent restricted-active-space self-consistent-field method for several excitation schemes and active orbital spaces with and without a dynamic core to address the effects systematically at different levels of approximation. For the ground-state many-electron wave functions, we compare the correlation energies with entropic measures of entanglement. A larger magnitude of the correlation energy does not always correspond to a larger value of the considered entanglement measures. To evaluate the impact of correlation in a process involving continua, we consider photoionization by attosecond pulses. The photoelectron spectra may be significantly affected by including a dynamical core.

Constructing Molecular π-Orbital Active Spaces for Multireference Calculations of Conjugated Systems.

Molecules with conjugated π systems often feature strong electron correlation and therefore require multireference methods for a reliable computational description. A key prerequisite for the successful application of such methods is the choice of a suitable active space. Herein the automated π-orbital space (PiOS) method for selecting active spaces for multireference calculations of conjugated π systems is presented. This approach allows the construction of small but effective active spaces based on Hückel theory. To demonstrate its performance, π → π* excitations for benzene, octatetraene, and free-base porphine are computed. In addition, this technique can be combined with the automated atomic valence active space method to compute excitations in complex systems with multiple conjugated fragments. This combined approach was used to generate two-dimensional potential energy surfaces for multiple electronic states associated with photoinduced electron-coupled double proton transfer in the blue-light-using flavin photoreceptor protein. These types of methods for the automated selection of active space orbitals are important for ensuring consistency and reproducibility of multireference approaches for a wide range of chemical and biological systems.