Second-order N-electron Valence State Perturbation Research Articles

Chelate Complexes of 3d Transition Metal Ions─A Challenge for Electronic-Structure Methods?

Different electronic-structure methods were assessed for their ability to predict two important properties of the industrially relevant chelating agent nitrilotriacetic acid (NTA): its selectivity with respect to six different first-row transition metal ions and the spin-state energetics of its complex with Fe(III). The investigated methods encompassed density functional theory (DFT), the random phase approximation (RPA), coupled cluster (CC) theory, and the auxiliary-field quantum Monte Carlo (AFQMC) method, as well as the complete active space self-consistent field (CASSCF) method and the respective on-top methods: second-order N-electron valence state perturbation theory (NEVPT2) and multiconfiguration pair-density functional theory (MC-PDFT). Different strategies for selecting active spaces were explored, and the density matrix renormalization group (DMRG) approach was used to solve the largest active spaces. Despite somewhat ambiguous multi-reference diagnostics, most methods gave relatively good agreement with experimental data for the chemical reactions connected to the selectivity, which only involved transition-metal complexes in their high-spin state. CC methods yielded the highest accuracy followed by range-separated DFT and AFQMC. We discussed in detail that even higher accuracies can be obtained with NEVPT2, under the prerequisite that consistent active spaces along the entire chemical reaction can be selected, which was not the case for reactions involving Fe(III). A bigger challenge for electronic-structure methods was the prediction of the spin-state energetics, which additionally involved lower spin states that exhibited larger multi-reference diagnostics. Conceptually different, typically accurate methods ranging from CC theory via DMRG-NEVPT2 in combination with large active spaces to AFQMC agreed well that the high-spin state is energetically significantly favored over the other spin states. This was in contrast to most DFT functionals and RPA which yielded a smaller stabilization and some common DFT functionals and MC-PDFT even predicting the low-spin state to be energetically most favorable.

Quasi-degenerate extension of local N-electron valence state perturbation theory with pair-natural orbital method based on localized virtual molecular orbitals.

Chemical phenomena involving near-degenerate electronic states, such as conical intersections or avoided crossing, can be properly described using quasi-degenerate perturbation theory. This study proposed a highly scalable quasi-degenerate second-order N-electron valence state perturbation theory (QD-NEVPT2) using the local pair-natural orbital (PNO) method. Our recent study showed an efficient implementation of the PNO-based state-specific NEVPT2 method using orthonormal localized virtual molecular orbitals (LVMOs) as an intermediate local basis. This study derived the state-coupling (or off-diagonal) terms to implement QD-NEVPT2 in an alternative manner to enhance efficiency based on the internally contracted basis and PNO overlap matrices between different references. To facilitate further acceleration, a local resolution-of-the-identity (RI) three-index integral generation algorithm was developed using LMOs and LVMOs. Although the NEVPT2 theory is considered to be less susceptible to the intruder-state problem (ISP), this study revealed that it can easily suffer from ISP when calculating high-lying excited states. We ameliorated this instability using the imaginary level shift technique. The PNO-QD-NEVPT2 calculations were performed on small organic molecules for the 30 lowest-lying states, as well as photoisomerization involving the conical intersection of 1,1-dimethyldibenzo[b,f] silepin with a cis-stilbene skeleton. These calculations revealed that the PNO-QD-NEVPT2 method yielded negligible errors compared to the canonical QD-NEVPT2 results. Furthermore, we tested its applicability to a large photoisomerization system using the green fluorescent protein model and the ten-state calculation of the large transition metal complex, showcasing that off-diagonal elements can be evaluated at a relatively low cost.

Spinless formulation of linearized adiabatic connection approximation and its comparison with the second order N-electron valence state perturbation theory.

The adiabatic connection (AC) approximation, along with its linearized variant AC0, was introduced as a method of obtaining dynamic correlation energy. When using a complete active space self-consistent field (CASSCF) wave function as a reference, the AC0 approximation is considered one of the most efficient multi-reference perturbation theories. It only involves the use of 1st- and 2nd-order reduced density matrices. However, some numerical results have indicated that the excitation energies predicted by AC0 are not as reliable as those from the second-order N-electron valence state perturbation theory (NEVPT2). In this study, we develop a spinless formulation of AC0 based on the Dyall Hamiltonian and provide a detailed comparison between AC0 and NEVPT2 approaches. We demonstrate the components within the correlation energy expressions that are common to both methods and those unique to either AC0 or NEVPT2. We investigate the role of the terms exclusive to NEVPT2 and explore the possibility of enhancing AC0's performance in this regard.

Optical Properties of Neutral F Centers in Bulk MgO with Density Matrix Embedding.

The optical spectra of neutral oxygen vacancies (F0 centers) in the bulk MgO lattice are investigated using density matrix embedding theory. The impurity Hamiltonian is solved with the complete active space self-consistent field and second-order n-electron valence state perturbation theory (NEVPT2-DMET) multireference methods. To estimate defect-localized vertical excitation energies at the nonembedding and thermodynamic limits, a double extrapolation scheme is employed. The extrapolated NEVPT2-DMET vertical excitation energy value of 5.24 eV agrees well with the experimental absorption maxima at 5.03 eV, whereas the excitation energy value of 2.89 eV at the relaxed triplet defect-localized state geometry overestimates the experimental emission at 2.4 eV by only nearly 0.5 eV, indicating the involvement of the triplet-singlet decay pathway.

Local Excitations of a Charged Nitrogen Vacancy in Diamond with Multireference Density Matrix Embedding Theory

We investigate the negatively charged nitrogen-vacancy center in diamond using periodic density matrix embedding theory (pDMET). To describe the strongly correlated excited states of this system, the complete active space self-consistent field (CASSCF) followed by n-electron valence state second-order perturbation theory (NEVPT2) was used as the impurity solver. Since the NEVPT2-DMET energies show a linear dependence on the inverse of the size of the embedding subspace, we performed an extrapolation of the excitation energies to the nonembedding limit using a linear regression. The extrapolated NEVPT2-DMET first triplet-triplet excitation energy is 2.31 eV and that for the optically inactive singlet-singlet transition is 1.02 eV, both in agreement with the experimentally observed vertical excitation energies of ∼2.18 eV and ∼1.26 eV, respectively. This is the first application of pDMET to a charged periodic system and the first investigation of the NV- defect using NEVPT2 for periodic supercell models.

Local N-electron valence state perturbation theory using pair-natural orbitals based on localized virtual molecular orbitals.

Second-order N-electron valence state perturbation theory (NEVPT2) is an exactly size-consistent and intruder-state-free multi-reference theory. To accelerate the NEVPT2 computation, Guo and Neese combined it with the local pair-natural orbital (PNO) method using the projected atomic orbitals (PAOs) as the underlying local basis [Guo et al., J. Chem. Phys. 144, 094111 (2016)]. In this paper, we report the further development of the PNO-NEVPT2 method using the orthonormal and non-redundant localized virtual molecular orbitals (LVMOs) instead of PAOs. The LVMOs were previously considered to perform relatively poor compared to PAOs because the resulting orbital domains were unacceptably large. Our prior work, however, showed that this drawback can be remedied by re-forming the domain construction scheme using differential overlap integrals [Saitow et al., J. Chem. Phys. 157, 084101 (2022)]. In this work, we develop further refinements to enhance the feasibility of using LVMOs. We first developed a two-level semi-local approach for screening out so-called weak-pairs to select or truncate the pairs for PNO constructions more flexibly. As a refinement specific to the Pipek-Mezey localization for LVMOs, we introduced an iterative scheme to truncate the Givens rotations using varying thresholds. We assessed the LVMO-based PNO-NEVPT2 method through benchmark calculations for linear phenyl alkanes, which demonstrate that it performs comparably well relative to the PAO-based approach. In addition, we evaluated the Co-C bond dissociation energies for the cobalamin derivatives composed of 200 or more atoms, which confirms that the LVMO-based method can recover more than 99.85% of the canonical NEVPT2 correlation energy.

Tetracoordinate Co(II) complexes with semi-coordination as stable single-ion magnets for deposition on graphene.

We present a theoretical and experimental study of two tetracoordinate Co(II)-based complexes with semi-coordination interactions, i.e., non-covalent interactions involving the central atom. We argue that such interactions enhance the thermal and structural stability of the compounds, making them appropriate for deposition on substrates, as demonstrated by their successful deposition on graphene. DC magnetometry and high-frequency electron spin resonance (HF-ESR) experiments revealed an axial magnetic anisotropy and weak intermolecular antiferromagnetic coupling in both compounds, supported by theoretical predictions from complete active space self-consistent field calculations complemented by N-electron valence state second-order perturbation theory (CASSCF-NEVPT2), and broken-symmetry density functional theory (BS-DFT). AC magnetometry demonstrated that the compounds are field-induced single-ion magnets (SIMs) at applied static magnetic fields, with slow relaxation of magnetization governed by a combination of quantum tunneling, Orbach, and direct relaxation mechanisms. The structural stability under ambient conditions and after deposition was confirmed by X-ray photoelectron spectroscopy (XPS) and Raman spectroscopy. Theoretical modeling by DFT of different configurations of these systems on graphene revealed n-type doping of graphene originating from electron transfer from the deposited molecules, confirmed by electrical transport measurements and Raman spectroscopy.

Large-Scale Benchmarking of Multireference Vertical-Excitation Calculations via Automated Active-Space Selection.

We have calculated state-averaged complete-active-space self-consistent-field (SA-CASSCF), multiconfiguration pair-density functional theory (MC-PDFT), hybrid MC-PDFT (HMC-PDFT), and n-electron valence state second-order perturbation theory (NEVPT2) excitation energies with the approximate pair coefficient (APC) automated active-space selection scheme for the QUESTDB benchmark database of 542 vertical excitation energies. We eliminated poor active spaces (20–40% of calculations) by applying a threshold to the SA-CASSCF absolute error. With the remaining calculations, we find that NEVPT2 performance is significantly impacted by the size of the basis set the wave functions are converged in, regardless of the quality of their description, which is a problem absent in MC-PDFT. Additionally, we find that HMC-PDFT is a significant improvement over MC-PDFT with the translated PBE (tPBE) density functional and that it performs about as well as NEVPT2 and second-order coupled cluster on a set of 373 excitations in the QUESTDB database. We optimized the percentage of SA-CASSCF energy to include in HMC-PDFT when using the tPBE on-top functional, and we find the 25% value used in tPBE0 to be optimal. This work is by far the largest benchmarking of MC-PDFT and HMC-PDFT to date, and the data produced in this work are useful as a validation of HMC-PDFT and of the APC active-space selection scheme. We have made all the wave functions produced in this work (orbitals and CI vectors) available to the public and encourage the community to utilize this data as a tool in the development of further multireference model chemistries.

Assessing the Performances of CASPT2 and NEVPT2 for Vertical Excitation Energies.

Methods able to simultaneously account for both static and dynamic electron correlations have often been employed, not only to model photochemical events but also to provide reference values for vertical transition energies, hence allowing benchmarking of lower-order models. In this category, both the complete-active-space second-order perturbation theory (CASPT2) and the N-electron valence state second-order perturbation theory (NEVPT2) are certainly popular, the latter presenting the advantage of not requiring the application of the empirical ionization-potential-electron-affinity (IPEA) and level shifts. However, the actual accuracy of these multiconfigurational approaches is not settled yet. In this context, to assess the performances of these approaches, the present work relies on highly accurate (±0.03 eV) aug-cc-pVTZ vertical transition energies for 284 excited states of diverse character (174 singlet, 110 triplet, 206 valence, 78 Rydberg, 78 n → π*, 119 π → π*, and 9 double excitations) determined in 35 small- to medium-sized organic molecules containing from three to six non-hydrogen atoms. The CASPT2 calculations are performed with and without IPEA shift and compared to the partially contracted (PC) and strongly contracted (SC) variants of NEVPT2. We find that both CASPT2 with IPEA shift and PC-NEVPT2 provide fairly reliable vertical transition energy estimates, with slight overestimations and mean absolute errors of 0.11 and 0.13 eV, respectively. These values are found to be rather uniform for the various subgroups of transitions. The present work completes our previous benchmarks focused on single-reference wave function methods ( J. Chem. Theory Comput. 2018, 14, 4360; J. Chem. Theory Comput. 2020, 16, 1711), hence allowing for a fair comparison between various families of electronic structure methods. In particular, we show that ADC(2), CCSD, and CASPT2 deliver similar accuracies for excited states with a dominant single-excitation character.

Excited States of Crystalline Point Defects with Multireference Density Matrix Embedding Theory

Accurate and affordable methods to characterize the electronic structure of solids are important for targeted materials design. Embedding-based methods provide an appealing balance in the trade-off between cost and accuracy─particularly when studying localized phenomena. Here, we use the density matrix embedding theory (DMET) algorithm to study the electronic excitations in solid-state defects with a restricted open-shell Hartree-Fock (ROHF) bath and multireference impurity solvers, specifically, complete active space self-consistent field (CASSCF) and n-electron valence state second-order perturbation theory (NEVPT2). We apply the method to investigate the electronic excitations in an oxygen vacancy (OV) on a MgO(100) surface and find absolute deviations within 0.05 eV between DMET using the CASSCF/NEVPT2 solver, denoted as CAS-DMET/NEVPT2-DMET, and the nonembedded CASSCF/NEVPT2 approach. Next, we establish the practicality of DMET by extending it to larger supercells for the OV defect and a neutral silicon vacancy in diamond where the use of nonembedded CASSCF/NEVPT2 is extremely expensive.

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.

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.

Locating conical intersections using the quasidegenerate partially and strongly contracted NEVPT2 methods

Conical intersections (CIs) play an important role in photochemistry because they sometimes govern the non-radiative decay process. However, accurate characterizations of CIs are not always straightforward. In this study, analytic gradients and interstate coupling vectors for the quasidegenerate partially and strongly contracted n-electron valence state second-order perturbation theory (QD-PC-NEVPT2 and QD-SC-NEVPT2) were developed and applied to locating CIs of benzene. The pilot application demonstrates that the results of both the methods are similar to the extended multistate complete active space second-order perturbation theory (XMS-CASPT2), while the lack of dynamic electron correlation resulted in a deviation of 0.7 eV.

Improvement of Ab Initio Ligand Field Theory by Means of Multistate Perturbation Theory.

Over the last few years, ab initio ligand field theory (AILFT) has evolved into an important tool for the extraction of ligand field models from ab initio calculations. The inclusion of dynamic correlation on top of complete active space self-consistent field (CASSCF) reference functions, which is important for accurate results, was so far realized at the level of second-order N-electron valence state perturbation theory (NEVPT2). In this work, we introduce two alternative methods for the inclusion of dynamic correlation into AILFT calculations, the second-order dynamic correlation dressed complete active space method (DCD-CAS(2)) and the Hermitian quasi-degenerate NEVPT2 (HQD-NEVPT2). These methods belong to the class of multistate perturbation theory approaches, which allow for the mixing of CASSCF states under the effect of dynamic correlation (state-mixing). The two new versions of AILFT were tested for a diverse set of transition-metal complexes. It was found that the multistate methods have, compared to NEVPT2, an AILFT fit with smaller root mean square deviations (rmsds) between ab initio and AILFT energies. A comparison of AILFT excitation energies with the experiment shows that for some systems, the agreement gets better at the multistate level because of the smaller rmsds. However, for some systems, the agreement gets worse, which could be attributed to a cancellation of errors at the NEVPT2 level that is partly removed at the multistate level. An investigation of trends in the extracted ligand field parameters shows that at the multistate level, the ligand field splitting Δ gets larger, whereas the Racah parameters B and C get smaller and larger, respectively. An investigation of the reasons for the observed improvement for octahedral CrIII halide complexes shows that the possibility of state-mixing relaxes constraints that are present at the NEVPT2 level and that keep Δ and B from following their individual preferences.

Spectroscopic and Computational Comparisons of Thiolate-Ligated Ferric Nonheme Complexes to Cysteine Dioxygenase: Second-Sphere Effects on Substrate (Analogue) Positioning.

Parallel spectroscopic and computational studies of iron(III) cysteine dioxygenase (CDO) and synthetic models are presented. The synthetic complexes utilize the ligand tris(4,5-diphenyl-1-methylimidazol-2-yl)phosphine (Ph2TIP), which mimics the facial three-histidine triad of CDO and other thiol dioxygenases. In addition to the previously reported [FeII(CysOEt)(Ph2TIP)]BPh4 (1; CysOEt is the ethyl ester of anionic l-cysteine), the formation and crystallographic characterization of [FeII(2-MTS)(Ph2TIP)]BPh4 (2) is reported, where the methyl 2-thiosalicylate anion (2-MTS) resembles the substrate of 3-mercaptopropionate dioxygenase (MDO). One-electron chemical oxidation of 1 and 2 yields ferric species that bind cyanide and azide anions, which have been used as spectroscopic probes of O2 binding in prior studies of FeIII-CDO. The six-coordinate FeIII-CN and FeIII-N3 adducts are examined with UV-vis absorption, electron paramagnetic resonance (EPR), and resonance Raman (rRaman) spectroscopies. In addition, UV-vis and rRaman studies of cysteine- and cyanide-bound FeIII-CDO are reported for both the wild-type (WT) enzyme and C93G variant, which lacks the Cys-Tyr cross-link that is present in the second coordination sphere of the WT active site. Density functional theory (DFT) and ab initio calculations are employed to provide geometric and electronic structure descriptions of the synthetic and enzymatic FeIII adducts. In particular, it is shown that the complete active space self-consistent field (CASSCF) method, in tandem with n-electron valence state second-order perturbation theory (NEVPT2), is capable of elucidating the structural basis of subtle shifts in EPR g values for low-spin FeIII species.

Analytic first-order derivatives of partially contracted n-electron valence state second-order perturbation theory (PC-NEVPT2)

A balanced treatment of dynamic and static electron correlation is important in computational chemistry, and multireference perturbation theory (MRPT) is able to do this at a reasonable computational cost. In this paper, analytic first-order derivatives, specifically gradients and dipole moments, are developed for a particular MRPT method, state-specific partially contracted n-electron valence state second-order perturbation theory (PC-NEVPT2). Only one linear equation needs to be solved for the derivative calculation if the Z-vector method is employed, which facilitates the practical application of this approach. A comparison of the calculated results with experimental geometrical parameters of O3 indicates excellent agreement although the calculated results for O3 - are slightly outside the experimental error bars. The 0-0 transition energies of various methylpyrimidines and trans-polyacetylene are calculated by performing geometry optimizations and seminumerical second-order geometrical derivative calculations. In particular, the deviations of 0-0 transition energies of trans-polyacetylene from experimental values are consistently less than 0.1 eV with PC-NEVPT2, indicating the reliability of the method. These results demonstrate the importance of adding dynamic electron correlation on top of methods dominated by static electron correlation and of developing analytic derivatives for highly accurate methods.

UV-Vis Absorption Spectroscopy of Polonium(IV) Chloride Complexes: An Electronic Structure Theory Study.

More than a hundred years after its discovery, the chemistry of the polonium radioelement is still largely unknown. However, it is quite clear that the properties of this heavy element ( Z = 84) may be affected by relativistic effects, in particular scalar relativistic effects and the so-called spin-orbit coupling (SOC). In this Article, we revisit the interpretation of UV-vis absorption spectra of polonium(IV) complexes in HCl medium, reported decades ago. From the data, two complexes were hypothesized, complex I with a maximum of absorption at 344 nm (at low HCl concentration) and complex II with a maximum at 418 nm (the only visible peak for HCl concentrations above 0.5 M). By monitoring the absorbance at 344 and 418 nm as a function of both the HCl concentration and the pH, complex I was formulated as [Po(OH)Cl x]3- x while complex II was formulated as [PoCl2+ x]2- x. In this work, we study the ground-state geometries of the [Po(OH)Cl x]3- x and [PoCl2+ x]2- x complexes for x = 4-2, i.e. for the most probable complexes, with density functional theory (DFT), demonstrating that solvation can remarkably change the geometries of such systems. The electronic excitation energies are then computed with time-dependent DFT (TD-DFT), second-order N-electron valence state perturbation theory (NEVPT2), and contracted spin-orbit configuration interaction (c-SOCI), showing (i) that the SOC must be at play to obtain excitation energies in the right energy domain and (ii) that the quantum chemical calculations point toward x = 4, i.e., toward the [Po(OH)Cl4]- and [PoCl6]2- complexes.

Ab Initio Wave Function-Based Determination of Element Specific Shifts for the Efficient Calculation of X-ray Absorption Spectra of Main Group Elements and First Row Transition Metals.

In this study, a detailed calibration of the performance of modern ab initio wave function methods in the domain of X-ray absorption spectroscopy (XAS) is presented. It has been known for some time that for a given level of approximation, for example, using time-dependent density functional theory (TD-DFT) in conjunction with a given basis set, there are systematic deviations of the calculated transition energies from their experimental values that depend on the functional, the basis set, and the chosen treatment of scalar relativistic effects. This necessitates a linear correlation for a given element/functional/basis set combination to be established before chemical applications can be performed. This is a laborious undertaking since it involves sourcing trustworthy experimental data, lengthy geometry optimizations, and detailed comparisons between theory and experiment. In this work, reference values for the element-specific shifts of all the first-row transition metal atoms and the main group elements C, N, O, F, Si, P, S, and Cl have been computed by using a protocol that is based on the complete active space configuration interaction in conjunction with second-order N-electron valence state perturbation theory (CASCI/NEVPT2). It is shown that by extrapolating the results to the basis set limit of the method and taking into account scalar relativistic effects via the second-order Douglas-Kroll-Hess (DKH2) corrections, the predicted element shifts are on average less than 0.75 eV across all the absorption edges and a very good agreement between theory and experiment in all the studied XAS cases is observed. The transferability of these errors to molecular systems is thoroughly investigated. The constructed CASCI/NEVPT2 database of element shifts is used to evaluate the performance and to automatically calibrate prior to comparison with the experiment two commonly used methods in X-ray spectroscopy, namely, DFT/Restricted open shell configuration interaction singles (DFT/ROCIS) and TD-DFT methods.

Development of Xi'an-CI package – applying the hole–particle symmetry in multi-reference electronic correlation calculations

ABSTRACTThis mini-review introduces our works on the Xi'an-CI (configuration interaction) package using graphical unitary group approach (GUGA). Taking advantage of the hole-particle symmetry in GUGA, the Galfand states used to span the CI space are classified into CI subspaces according to the number of holes and particles, and the coupling coefficients used to calculate Hamiltonian matrix elements could be factorised into the segment factors in the hole, active and external spaces. An efficient multi-reference CI with single and double excitations (MRCISD) algorithm is thus developed that reduces the storage requirement and increases the number of correlated electrons significantly. The hole–particle symmetry also gives rise to a doubly contracted MRCISD approach. Moreover, the internally contracted Gelfand states are defined within the CI subspace arising from the hole–particle symmetry, which makes the implementation of internally contracted MRCISD in the framework of GUGA possible. In addition to MRCISD, the development of multi-reference second-order perturbation theory (MRPT2) also benefits from the hole–particle symmetry. A configuration-based MRPT2 algorithm is proposed and extended to the multi-state n-electron valence-state second-order perturbation theory.

Further development of SDSPT2 for strongly correlated electrons

ABSTRACTThe static-dynamic-static multi-state multi-reference second-order perturbation theory (SDS-MS-MRPT2 or SDSPT2 in short), which employs a multi-partitioned (state-dependent) Møller–Plesset-like diagonal operator for the zeroth-order Hamiltonian and individual configuration state functions for the perturbers, is reconsidered by making use of the state-universal Dyall Hamiltonian for the zeroth-order Hamiltonian and the union of spin-adapted, internally contracted basis functions generated from all the reference functions for the perturbers. The graphical unitary group approach along with particle–hole correspondence is employed in the evaluation of Hamiltonian matrix elements. This particular variant of SDSPT2 can also be regarded as an extension of the spin-adapted, partially contracted variant of NEVPT2 (n-electron valence state second-order perturbation theory). The performance of SDSPT2 is examined by taking C2, LiF, and Fe(C5H5)2 as examples. The results show that the low-lying states of the systems can be obtained very accurately.