Combining Internally Contracted States and Matrix Product States To Perform Multireference Perturbation Theory

Multireference Correlated Oscillator Strengths from Adiabatic Connection Approaches Based on Extended Random Phase Approximation

We show that accurate oscillator strengths can be obtained from adiabatic connection (AC) approaches based on the extended random phase approximation (ERPA) combined with multireference (complete active space, CAS) wave functions. The oscillator strengths calculated using the perturbation-corrected ERPA transition density matrices, proposed in this work, and the excitation energies calculated with recently introduced AC correlation energy methods, AC0 and AC0D, compete with accuracy in the perturbational CASPT2 approach and require less computational effort. AC0 and AC0D methods scale more favorably with the number of active orbitals than multiconfigurational perturbation approaches like CASPT2 and NEVPT2 thanks to their dependence on reduced density matrices up to the order of 2. Importantly, the newly developed approach for computing correlated transition dipole moments does not entail any additional costs, as all intermediate quantities become available when AC0 energies are being computed. We also test the performance of the recently proposed AC method corrected for the negative-transition contributions to the correlation energy, AC0D, for triplet excitation energies. Similarly, as for the singlet excitations, the correction improves the performance of the AC0 method, particularly for the low-lying excited states.

Renormalized-Residue-Based Multireference Configuration Interaction Method for Strongly Correlated Systems

The implementation of multireference configuration interaction (MRCI) methods in quantum systems with large active spaces is hindered by the expansion of configuration bases or the intricate handling of reduced density matrices (RDMs). In this work, we present a spin-adapted renormalized-residue-based MRCI (RR-MRCI) approach that leverages renormalized residues to effectively capture the entanglement between active and inactive orbitals. This approach is reinforced by a novel efficient algorithm, which also facilitates an efficient deployment of spin-adapted matrix product state MRCI (MPS-MRCI). The RR-MRCI framework possesses several advantages: (1) It considers the orbital entanglement and utilizes highly compressed MPS structure, improving computational accuracy and efficiency compared with internally contracted (ic) MRCI. (2) Utilizing small-sized buffer environments of a few external orbitals as probes based on quantum information theory, it enhances computational efficiency over MPS-MRCI and offers potential application to large molecular systems. (3) The RR framework can be implemented in conjunction with ic-MRCI, eliminating the need for high-rank RDMs, by using distinct renormalized residues. We evaluated this method across nine diverse molecular systems, including Cu2O22+ with an active space of (24e,24o) and two complexes of lanthanide and actinide with active space (38e,36o), demonstrating the method's versatility and efficacy.

Block2: A comprehensive open source framework to develop and apply state-of-the-art DMRG algorithms in electronic structure and beyond

block2 is an open source framework to implement and perform density matrix renormalization group and matrix product state algorithms. Out-of-the-box it supports the eigenstate, time-dependent, response, and finite-temperature algorithms. In addition, it carries special optimizations for ab initio electronic structure Hamiltonians and implements many quantum chemistry extensions to the density matrix renormalization group, such as dynamical correlation theories. The code is designed with an emphasis on flexibility, extensibility, and efficiency and to support integration with external numerical packages. Here, we explain the design principles and currently supported features and present numerical examples in a range of applications.

Intruder-free cumulant-truncated driven similarity renormalization group second-order multireference perturbation theory

Accurate multireference electronic structure calculations are important for constructing potential energy surfaces. Still, even in the case of low-scaling methods, their routine use is limited by the steep growth of the computational and storage costs as the active space grows. This is primarily due to the occurrence of three- and higher-body density matrices or, equivalently, their cumulants. This work examines the effect of various cumulant truncation schemes on the accuracy of the driven similarity renormalization group second-order multireference perturbation theory. We test four different levels of three-body reduced density cumulant truncations that set different classes of cumulant elements to zero. Our test cases include the singlet-triplet gap of CH2, the potential energy curves of the XΣg+1 and AΣu+3 states of N2, and the singlet-triplet splittings of oligoacenes. Our results show that both relative and absolute errors introduced by these cumulant truncations can be as small as 0.5 kcal mol-1 or less. At the same time, the amount of memory required is reduced from O(NA6) to O(NA5), where NA is the number of active orbitals. No additional regularization is needed to prevent the intruder state problem in the cumulant-truncated second-order driven similarity renormalization group multireference perturbation theory methods.

New Density Matrix Renormalization Group Approaches for Strongly Correlated Systems Coupled with Large Environments

Thanks to the high compression of the matrix product state (MPS) form of the wave function and the efficient site-by-site iterative sweeping optimization algorithm, the density matrix normalization group (DMRG) and its time-dependent variant (TD-DMRG) have been established as powerful computational tools in accurately simulating the electronic structure and quantum dynamics of strongly correlated molecules with a large number (10^1-2) of quantum degrees of freedom (active orbitals or vibrational modes). However, the quantitative characterization of the quantum many-body behaviors of realistic strongly correlated systems requires a further consideration of the interaction between the embedded active subsystem and the remaining correlated environment, e.g., a larger number (10^2-3) of external orbitals in electronic structure or infinite condensed-phase phononic modes in nucleus dynamics. To this end, we introduced three new post-DMRG and TD-DMRG approaches, namely (1) DMRG2sCI-MRCI and DMRG2sCI-ENPT by the reconstruction of selected configuration interaction (sCI) type of compact reference function from DMRG coefficients and the use of externally contracted MRCI (multireference configuration interaction) and Epstein-Nesbet perturbation theory (ENPT), without recourse to the expensive high order n-electron reduced density matrices (n-RDMs). (2) DMRG combined with RR-MRCI (renormalized residue-based MRCI), which improves the computational accuracy and efficiency of internally contracted (ic) MRCI by renormalizing the contracted bases with small-sized buffer environment(s) of a few external orbitals as probes based on quantum information theory. (3) HM (hierarchical mapping)-TD-DMRG in which a large environment is reduced to a small number of renormalized environmental modes (which accounts for the most vital system-environment interactions) through stepwise mapping transformation. These advances extend the efficacy of highly accurate DMRG/TD-DMRG computations to the quantitative characterization of the electronic structure and quantum dynamics in realistic strongly correlated systems coupled with large environments and are reviewed in this paper.

Kylin 1.0: An ab-initio density matrix renormalization group quantum chemistry program

The accurate evaluation of electron correlations is highly necessary for the proper descriptions of the electronic structures in strongly correlated molecules, ranging from bond-dissociating molecules, polyradicals, to large conjugated molecules and transition metal complexes. For this purpose, in this paper, a new ab-initio quantum chemistry program Kylin 1.0 for electron correlation calculations at various quantum many-body levels, including configuration interaction (CI), perturbation theory (PT), and density matrix renormalization group (DMRG), is presented. Furthermore, fundamental quantum chemistry methods such as Hartree-Fock self-consistent field (HF-SCF) and the complete active space SCF (CASSCF) are also implemented. The Kylin 1.0 program possesses the following features: (1) a matrix product operator (MPO) formulation-based efficient DMRG implementation for describing static electron correlation within a large active space composed of more than 100 orbitals, supporting both U 1 n × U 1 S z $$ \mathrm{U}{(1)}_{\mathrm{n}}\times \mathrm{U}{(1)}_{{\mathrm{S}}_{\mathrm{z}}} $$ and U 1 n × SU 2 S $$ \mathrm{U}{(1)}_{\mathrm{n}}\times \mathrm{SU}{(2)}_{\mathrm{S}} $$ symmetries; (2) an efficient second-order DMRG-self-consistent field (SCF) implementation; (3) an externally contracted multi-reference CI (MRCI) and Epstein-Nesbet PT with DMRG reference wave functions for including the remaining dynamic electron correlation outside the large active spaces. In this paper, we introduce the capabilities and numerical benchmark examples of the Kylin 1.0 program.

Accurate Property Prediction by Second Order Perturbation Theory: The REMP and OO-REMP Hybrids

The prediction of molecular properties such as equilibrium structures or vibrationalwavenumbers is a routine task in computational chemistry. If very high accuracy is required, however, the use of computationally demanding ab initio wavefunction methods is mandatory. We present property calculations utilizing the REMP and OO-REMP hybrid perturbation theories showing that with the latter approach, very accurate results are obtained at second order in perturbation theory. Specifically, equilibrium structures and harmonic vibrational wavenumbers as well as dipole moments of closed and open shell molecules were calculated and compared to the best available experimental results or very accurate calculations.OO-REMP is capable of predicting bond lengths of small closed and open shell molecules with an accuracy of 0.2 pm and 0.5 pm, respectively, often within the range of experimental uncertainty. Equilibrium harmonic vibrational wavenumbers are predicted with an accuracy better than 20 cm−1 . Dipole moments of small closed and open shell molecules are reproduced with a relative error of less than 3 %. Across all investigated properties it turns out that a 20 %:80 % MP:RE mixing ratio consistently provides the best results. This is in line with our previous findings featuring closed and open shell reaction energies.

Low-Lying Excited States of Natural Carotenoids Viewed by Ab Initio Methods

Low-lying excited states of carotenoids (the optically dark 2A_g^- and bright 1B_u⁺) are deeply involved in energy transfer processes in photosynthetic antennas, such as light harvesting and non-photochemical quenching. Though any ab initio modeling of these phenomena has to rely on precise energies of the carotenoid electronic states, the accurate evaluation of these states remains a challenging problem due to their different natures. The paper aims to study the accuracy of the excitation energies of the low-lying excited states of certain open- and closed-chain carotenoids obtained by a state-of-the-art multireference approach for electronic structure calculation. Here, density matrix renormalization group SCF (DMRGSCF) and a perturbative approach based on driven similarity renormalization group second-order multireference perturbation theory (DSRG-MRPT2) were used to treat the static and dynamic correlation, respectively. Nuclear geometries of the electronic states were optimized with DFT-based approaches. It is demonstrated that spin-flip TD-DFT can replace multiconfigurational methods for the geometry optimization of the 2A_g^- state but not for the calculation of the excitation energy. Adiabatic excitation energies to the 1B_u⁺ state were shown to be within a margin of 1000 cm^-1 with an appropriate flow parameter value. Adiabatic excitation energies to the 2A_g^- state for the open-chain carotenoids lie within a range of experimental values (taking into account the broad range of experimental estimates); for the closed-chain ones, the error does not exceed 2000 cm^-1. Ab initio stationary (1A_g^- → 1B_u⁺) and transient (2A_g^- → 1B_u⁺) absorption spectra were modeled for violaxanthin and lycopene, and these spectra showed good agreement with the experimental ones both in terms of the vibronic structure and the transition energies.

Post-Density Matrix Renormalization Group Methods for Describing Dynamic Electron Correlation with Large Active Spaces

The ab initio density matrix renormalization group (DMRG) method has been well-established and has become one of the most accurate numerical methods for the precise electronic structure solution of large active spaces. In the past few years, to capture the missing dynamic correlation, various post-DMRG approaches have been proposed through the combination of DMRG and multireference quantum chemical methods or density functional theory. With this in mind, this work provides a brief overview of ab initio DMRG principles and the new developments within post-DMRG methods. For clarity, post-DMRG methods are classified into two main categories depending on whether high-order n-electron reduced density matrices are used, and their merits and disadvantages are properly discussed. Finally, we conclude by discussing unsolved bottlenecks and giving development perspectives of post-DMRG approaches, which are expected to yield quantitative descriptions of complex electronic structures in large strongly correlated molecules and materials.

An efficient implementation of the NEVPT2 and CASPT2 methods avoiding higher-order density matrices

A factorization of the matrix elements of the Dyall Hamiltonian in N-electron valence state perturbation theory allowing their evaluation with a computational effort comparable to the one needed for the construction of the third-order reduced density matrix at the most is presented. Thus, the computational bottleneck arising from explicit evaluation of the fourth-order density matrix is avoided. It is also shown that the residual terms arising in the case of an approximate complete active space configuration interaction solution and containing even the fifth-order density matrix for two excitation classes can be evaluated with little additional effort by choosing again a favorable factorization of the corresponding matrix elements. An analogous argument is also provided for avoiding the fourth-order density matrix in complete active space second-order perturbation theory. Practical calculations indicate that such an approach leads to a considerable gain in computational efficiency without any compromise in numerical accuracy or stability.

Analytic Energy Gradients for the Driven Similarity Renormalization Group Multireference Second-Order Perturbation Theory

We derive analytic energy gradients of the driven similarity renormalization group (DSRG) multireference second-order perturbation theory (MRPT2) using the method of Lagrange multipliers. In the Lagrangian, we impose constraints for a complete-active-space self-consistent-field reference wave function and the semicanonical orthonormal molecular orbitals. Solving for the associated Lagrange multipliers is found to share the same asymptotic scaling of a single DSRG-MRPT2 energy computation. A pilot implementation of the DSRG-MRPT2 analytic gradients is used to optimize the geometry of the singlet and triplet states of p-benzyne. The equilibrium bond lengths and angles are similar to those computed via other MRPT2s and Mukherjee's multireference coupled cluster theory. An approximate DSRG-MRPT2 method that neglects the contributions of the three-body density cumulant is found to introduce negligible errors in the geometry of p-benzyne, lending itself to a promising low-cost approach for molecular geometry optimizations using large active spaces.

Multireference Algebraic Diagrammatic Construction Theory for Excited States: Extended Second-Order Implementation and Benchmark

We present an implementation and benchmark of new approximations in multireference algebraic diagrammatic construction theory for simulations of neutral electronic excitations and UV/vis spectra of strongly correlated molecular systems (MR-ADC). Following our work on the first-order MR-ADC approximation [J. Chem. Phys. 2018, 149, 204113], we report the strict and extended second-order MR-ADC methods (MR-ADC(2) and MR-ADC(2)-X) that combine the description of static and dynamic electron correlation in the ground and excited electronic states without relying on state-averaged reference wave functions. We present an extensive benchmark of the new MR-ADC methods for excited states in several small molecules, including the carbon dimer, ethylene, and butadiene. Our results demonstrate that, for weakly correlated electronic states, the MR-ADC(2) and MR-ADC(2)-X methods outperform the third-order single-reference ADC approximation and are competitive with the results from equation-of-motion coupled cluster theory. For states with multireference character, the performance of the MR-ADC methods is similar to that of an N-electron valence perturbation theory. In contrast to conventional multireference perturbation theories, the MR-ADC methods have many attractive features, such as a straightforward and efficient calculation of excited-state properties and a direct access to excitations outside of the frontier (active) orbitals.

Modern multireference methods and their application in transition metal chemistry

Transition metal chemistry is a challenging playground for quantum chemical methods owing to the simultaneous presence of static and dynamic electron correlation effects in many systems. Wavefunction based multireference (MR) methods constitute a physically sound and systematically improvable Ansatz to deal with this complexity but suffer from some conceptual difficulties and high computational costs. The latter problem partially arises from the unfavorable scaling of the Full Configuration Interaction (Full-CI) problem which in the majority of MR methods is solved for a subset of the molecular orbital space, the so-called active space. In the last years multiple methods such as modern variants of selected CI, Full-CI Quantum Monte Carlo (FCIQMC) and the density matrix renormalization group (DMRG) have been developed that solve the Full-CI problem approximately for a fraction of the computational cost required by conventional techniques thereby significantly extending the range of applicability of modern MR methods. This perspective review outlines recent advancements in the field of MR electronic structure methods together with the resulting chances and challenges for theoretical studies in the field of transition metal chemistry. In light of its emerging importance a special focus is put on the selection of adequate active spaces and the concomitant development of numerous selection aides in recent years.

Toward an Accurate Ab Initio Description of Low-Lying Singlet Excited States of Polyenes

The low-lying excited states of carotenoids play a crucial role in many important biophysical processes such as photosynthesis. Most of these excited states are strongly correlated, which makes them both challenging for a qualitative ab initio description and an engaging model system for trying out emerging multireference methods. Among these methods, driven similarity renormalization group (DSRG) and its perturbative version (DSRG-MRPT2) are especially attractive in terms of both accuracy and moderate numerical complexity. In this paper, we applied density matrix renormalization group (DMRG) followed by DSRG-MRPT2 for the calculation of vertical and adiabatic excitation energies into the 2A_g^-, 1B_u^-, and 1B_u⁺ electronic states of polyenes containing from 8 to 13 conjugating double bonds acting as a model for natural carotenoids. It was shown that the DSRG flow parameter should be adjusted to ensure both the energy convergence with respect to it and the agreement with the experimental data. With the increased flow parameter, the proposed combination of methods provides a reasonable agreement with the experiment. The deviations of the adiabatic excitation energies are less than 1000 cm^-1 for the 2A_g^- and less than 3000 cm^-1 for the excited states of the B_u symmetry, which in terms of accuracy significantly outperforms the N-electron valence state perturbation theory. At the same time, DSRG-MRPT2 is shown to be robust with respect to variation of quality of the DMRG reference wave function such as the orbital optimization or the number of electronic states in the averaging.

Approximations of density matrices in N-electron valence state second-order perturbation theory (NEVPT2). II. The full rank NEVPT2 (FR-NEVPT2) formulation

In Paper I, the performances of pre-screening (PS), extended PS (EPS), and cumulant (CU) approximations to the fourth-order density matrix were examined in the context of second-order N-electron valence state perturbation theory (NEVPT2). It has been found that the CU, PS, and even EPS approximations with loose thresholds may introduce intruder states. In the present work, the origin of these "false intruder" states introduced by approximated density matrices is discussed. Canonical NEVPT2 implementations employ a rank reduction trick. By analyzing its residual error, we find that the omission of the rank reduction leads to a more stable multireference perturbation theory for incomplete active space reference wave functions. Such a full rank (FR)-NEVPT2 formulation is equivalent to the conventional NEVPT2 method for the complete active space self-consistent field/complete active space configuration interaction reference wave function. A major drawback of the FR-NEVPT2 formulation is the necessity of the fifth-order density matrix. To avoid the construction of the high-order density matrices, the combination of the FR-NEVPT2 with the CU approximation is studied. However, we find that the CU approximation remains problematic as it still introduces intruder states. The question of how to robustly and efficiently perform internally contracted multireference perturbation theories with approximate densities remains a challenging field of investigation.

Approximations of density matrices in N-electron valence state second-order perturbation theory (NEVPT2). I. Revisiting the NEVPT2 construction

Over the last decade, the second-order N-electron valence state perturbation theory (NEVPT2) has developed into a widely used multireference perturbation method. To apply NEVPT2 to systems with large active spaces, the computational bottleneck is the construction of the fourth-order reduced density matrix. Both its generation and storage become quickly problematic beyond the usual maximum active space of about 15 active orbitals. To reduce the computational cost of handling fourth-order density matrices, the cumulant approximation (CU) has been proposed in several studies. A more conventional strategy to address the higher-order density matrices is the pre-screening approximation (PS), which is the default one in the ORCA program package since 2010. In the present work, the performance of the CU, PS, and extended PS (EPS) approximations for the fourth-order density matrices is compared. Following a pedagogical introduction to NEVPT2, contraction schemes, as well as the approximations to density matrices, and the intruder state problem are discussed. The CU approximation, while potentially leading to large computational savings, virtually always leads to intruder states. With the PS approximation, the computational savings are more modest. However, in conjunction with conservative cutoffs, it produces stable results. The EPS approximation to the fourth-order density matrices can reproduce very accurate NEVPT2 results without any intruder states. However, its computational cost is not much lower than that of the canonical algorithm. Moreover, we found that a good indicator of intrude states problems in any approximation to high order density matrices is the eigenspectra of the Koopmans matrices.

Transcorrelated density matrix renormalization group

We introduce the transcorrelated Density Matrix Renormalization Group (tcDMRG) theory for the efficient approximation of the energy for strongly correlated systems. tcDMRG encodes the wave function as a product of a fixed Jastrow or Gutzwiller correlator and a matrix product state. The latter is optimized by applying the imaginary-time variant of time-dependent (TD) DMRG to the non-Hermitian transcorrelated Hamiltonian. We demonstrate the efficiency of tcDMRG with the example of the two-dimensional Fermi-Hubbard Hamiltonian, a notoriously difficult target for the DMRG algorithm, for different sizes, occupation numbers, and interaction strengths. We demonstrate fast energy convergence of tcDMRG, which indicates that tcDMRG could increase the efficiency of standard DMRG beyond quasi-monodimensional systems and provides a generally powerful approach toward the dynamic correlation problem of DMRG.

Ab Initio Study of Low-Lying Excited States of Carotenoid-Derived Polyenes

Knowledge about excited states of carotenoids is essential for understanding photophysical processes underlying photosynthesis. However, due to the presence of a large number of optically dark states, experimental study of the excited-state manifold is limited to a significant extent. In this paper, we apply high-level ab initio quantum chemical methods to study the low-lying excited states of polyenes containing from 8 to 13 conjugated double bonds, which serve as a model for natural carotenoids. Vertical and adiabatic excitation energies from the ground 1A_g^- state to the excited 2A_g^-, 1B_u⁺, and 1B_u^- states were evaluated by means of density matrix renormalization group (DMRG) with NEVPT2 perturbative correction. The energies of all excited states are highly sensitive to nuclear geometry, especially the 2A_g^- state. Thus, the 2A_g^- and 1B_u⁺ states interchange their relative positions upon geometry relaxation, while the vertical excitation energy to the 2A_g^- state is rather high. At the same time, the 1B_u^- state energy is shown to be higher than other studied excited states at any geometry. With relaxed geometries of the excited states, absorption and transient absorption spectra were calculated within the Franck-Condon approximation bridging the gap between experimental spectroscopic data and computational results.

Reduced scaling extended multi-state CASPT2 (XMS-CASPT2) using supporting subspaces and tensor hyper-contraction

We present a reduced scaling formulation of the extended multi-state CASPT2 (XMS-CASPT2) method, which is based on our recently developed state-specific CASPT2 (SS-CASPT2) formulation using supporting subspaces and tensor hyper-contraction. By using these two techniques, the off-diagonal elements of the effective Hamiltonian can be computed with only O(N³) operations and O(N²) memory, where N is the number of basis functions. This limits the overall computational scaling to O(N⁴) operations and O(N²) memory. Thus, excited states can now be obtained at the same reduced (relative to previous algorithms) scaling we achieved for SS-CASPT2. In addition, we also investigate how the energy denominators can be factorized with the Laplace quadrature when some of the denominators are negative, which is critical for excited state calculations. An efficient implementation of the method has been developed using graphical processing units while also exploiting spatial sparsity in tensor operations. We benchmark the accuracy of the new method by comparison to non-THC formulated XMS-CASPT2 for the excited states of various molecules. In our tests, the THC approximation introduces negligible errors (≈0.01 eV) compared to the non-THC reference method. Scaling behavior and computational timings are presented to demonstrate performance. The new method is also interfaced with quantum mechanics/molecular mechanics (QM/MM). In an example study of green fluorescent protein, we show how the XMS-CASPT2 potential energy surfaces and excitation energies are affected by increasing the size of the QM region up to 278 QM atoms with more than 2300 basis functions.

Near-Linear Scaling in DMRG-Based Tailored Coupled Clusters: An Implementation of DLPNO-TCCSD and DLPNO-TCCSD(T)

We present a new implementation of density matrix renormalization group based tailored coupled clusters method (TCCSD), which employs the domain-based local pair natural orbital approach (DLPNO). Compared to the previous local pair natural orbital (LPNO) version of the method, the new implementation is more accurate, offers more favorable scaling, and provides more consistent behavior across the variety of systems. On top of the singles and doubles, we include the perturbative triples correction (T), which is able to retrieve even more dynamic correlation. The methods were tested on three systems: tetramethyleneethane, oxo-Mn(Salen), and iron(II)-porphyrin model. The first two were revisited to assess the performance with respect to LPNO-TCCSD. For oxo-Mn(Salen), we retrieved between 99.8 and 99.9% of the total canonical correlation energy which is an improvement of 0.2% over the LPNO version in less than 63% of the total LPNO runtime. Similar results were obtained for iron(II)-porphyrin. When the perturbative triples correction was employed, irrespective of the active space size or system, the obtained energy differences between two spin states were within the chemical accuracy of 1 kcal/mol using the default DLPNO settings.

A multireference coupled-electron pair approximation combined with complete-active space perturbation theory in local pair-natural orbital framework

The Complete-Active Space Second-order Perturbation Theory (CASPT2) has been one of the most widely-used methods for reliably calculating electronic structures of multireference systems. Because of its lowest level treatment of dynamic correlation, it has a high computational feasibility; however, its accuracy in some cases falls short of needs. Here, as a simple yet higher-order alternative, we introduce a hybrid theory of the CASPT2 and a multireference variant of the Coupled-Electron Pair Approximation (CEPA), which is a class of high level correlation theory. A central feature of our theory (CEPT2) is to use the two underlying theories for describing different divisions of correlation components based on the full internal contraction framework. The external components, which usually give a major contribution to the dynamic correlation, are intensively described using the CEPA Ansatz, while the rests are treated at the CASPT2 level. Furthermore, to drastically reduce the computational demands, we have incorporated the pair-natural orbital (PNO) method into our multireference implementations. This development, thus, requires highly complex derivations and coding, while it has been largely facilitated with an automatic expression and code generation technique. To highlight the accuracy of the CEPT2 approach and to assess the errors caused by the PNO truncation, benchmark calculations are shown on small- to medium-size molecules, illustrating the high accuracy of the present CEPT2 model. By tightening the truncation thresholds, the PNO-CEPT2 energy converges toward the canonical counterpart and is more accurate than that of PNO-CASPT2 as long as the same truncation thresholds are used.

An ab initio multireference study of reductive eliminations from organoferrates(iii) in the gas-phase: it is all about the spin state

The reductive elimination reaction from organoferrates(iii) of the composition [FeR₃R′]⁻ is studied by state-of-the-art multireference electronic structure calculations.

Toward the efficient local tailored coupled cluster approximation and the peculiar case of oxo-Mn(Salen)

We introduce a new implementation of the coupled cluster method with single and double excitations tailored by the matrix product state wave functions (DMRG-TCCSD), which employs the local pair natural orbital (LPNO) approach. By exploiting locality in the coupled cluster stage of the calculation, we were able to remove some of the limitations that hindered the application of the canonical version of the method to larger systems and/or with larger basis sets. We assessed the accuracy of the approximation using two systems: tetramethyleneethane (TME) and oxo-Mn(Salen). Using the default cut-off parameters, we were able to recover over 99.7% and 99.8% of the canonical correlation energy for the triplet and singlet state of TME, respectively. In the case of oxo-Mn(Salen), we found that the amount of retrieved canonical correlation energy depends on the size of the complete active space (CAS)-we retrieved over 99.6% for the larger 27 orbital CAS and over 99.8% for the smaller 22 orbital CAS. The use of LPNO-TCCSD allowed us to perform these calculations up to quadruple-ζ basis set, amounting to 1178 basis functions. Moreover, we examined dependence of the ground state of oxo-Mn(Salen) on the CAS composition. We found that the inclusion of 4d_xy orbital plays an important role in stabilizing the singlet state at the DMRG-CASSCF level via double-shell effect. However, by including dynamic correlation, the ground state was found to be triplet regardless of the size of the basis set or the composition of CAS, which is in agreement with previous findings by canonical DMRG-TCCSD in smaller basis.

Density matrix renormalization group pair-density functional theory (DMRG-PDFT): singlet-triplet gaps in polyacenes and polyacetylenes

The density matrix renormalization group (DMRG) is a powerful method to treat static correlation.

The density matrix renormalization group (DMRG) is a powerful method to treat static correlation. Here we present an inexpensive way to calculate correlation energy starting from a DMRG wave function using pair-density functional theory (PDFT). We applied this new approach, called DMRG-PDFT, to study singlet–triplet gaps in polyacenes and polyacetylenes that require active spaces larger than the feasibility limit of the conventional complete active-space self-consistent field (CASSCF) method. The results match reasonably well with the most reliable literature values and have only a moderate dependence on the compression of the initial DMRG wave function. Furthermore, DMRG-PDFT is significantly less expensive than other commonly applied ways of adding additional correlation to DMRG, such as DMRG followed by multireference perturbation theory or multireference configuration interaction.