Configuration State Functions Research Articles

Unified Implementation of Relativistic Wave Function Methods: 4C-iCIPT2 as a Showcase.

In parallel to the unified construction of relativistic Hamiltonians based solely on physical arguments (J. Chem. Phys. 2024, 160, 084111), a unified implementation of relativistic wave function methods is achieved here via programming techniques (e.g., template metaprogramming and polymorphism in C++). That is, once the code for constructing the Hamiltonian matrix is made ready, all the rest can be generated automatically from existing templates used for the nonrelativistic counterparts. This is facilitated by decomposing a second-quantized relativistic Hamiltonian into diagrams that are topologically the same as those required for computing the basic coupling coefficients between spin-free configuration state functions (CSF). Moreover, both time reversal and binary double point group symmetries can readily be incorporated into molecular integrals and Hamiltonian matrix elements. The latter can first be evaluated in the space of (randomly selected) spin-dependent determinants and then transformed to that of spin-dependent CSFs, thanks to simple relations in between. As a showcase, we consider here the no-pair four-component relativistic iterative configuration interaction with selection and perturbation correction (4C-iCIPT2), which is a natural extension of the spin-free iCIPT2 (J. Chem. Theory Comput. 2021, 17, 949), and can provide near-exact numerical results within the manifold of positive energy states (PES), as demonstrated by numerical examples.

Second order Rayleigh–Schrödinger perturbation theory for the Grasp2018 package: Core–core correlations

GRASP package is based on the relativistic configuration interaction in which accurate calculations, accounting for valence, valence–valence, core–valence, core and core–core electron correlations, often rely on massive CSF expansions. This paper presents a further development of the method based on the second-order perturbation theory for finding the most important CSFs that have the greatest influence on the core–valence, core and core–core correlations. This method is based on a combination of the relativistic configuration interaction method and the stationary second-order Rayleigh–Schrödinger many-body perturbation theory in an irreducible tensorial form [G. Gaigalas, P. Rynkun, and L. Kitovienė, Second-order Rayleigh–Schrödinger perturbation theory for the GRASP2018 package: Core–valence correlations, Lith. J. Phys. 64(1), 20–39 (2024), https://doi.org/10.3952/physics.2024.64.1.3, and G. Gaigalas, P. Rynkun, and L. Kitovienė, Second-order Rayleigh–Schrödinger perturbation theory for the GRASP2018 package: Core correlations, Lith. J. Phys. 64(2), 73–81 (2024), https://doi.org/10.3952/physics.2024.64.2.1]. In this extension, the perturbation theory takes into account electron core–valence, core and core–core correlations, where an atom or ion has any number of valence electrons, for calculation of energy spectra and other properties. Meanwhile, the rest of the correlations are taken into account in a traditional way. This allows a significant reduction of the space of the configuration state function for complex atoms and ions. We also demonstrate how this method works for calculations of the energy structure and E1 transition properties of Fe XV ion.

Ab Initio Manganese Kα and Kβ Energy Eigenvalues, Shake-Off Probabilities, Auger Rates, with Convergence Tests.

This work presents ab initio calculations for the Kα spectrum of manganese (Z = 25, [Ar]3d54s2), a highly complex system due to the five open orbitals in the 3d shell. The spectrum is composed of the canonical diagram line [1s]→[2p] and shake-off satellite lines [1snl]→[2pnl] (nl∈{2s,2p,3s,3p,3d,4s}), where square brackets denote a hole state. The multiconfiguration Dirac-Hartree-Fock method with the active set approach provides the initial and final atomic wavefunctions. Results are presented as energy eigenvalue spectra for the diagram and satellite transitions. The calculated wavefunctions include over one hundred million configuration state functions and over 280,000 independent transition energies for the seven sets of spectra considered. Shake-off probabilities and Auger transition rates determine satellite intensities. The number of configuration state functions ensures highly-converged wavefunctions. Several measures of convergence demonstrate convergence in the calculated parameters. We obtain convergence of the transition energies in all eight transitions to within 0.06 eV and shake-off probabilities to within 4.5%.

Restricted Open-Shell Hartree-Fock Method for a General Configuration State Function Featuring Arbitrarily Complex Spin-Couplings.

In this work, we present a general spin restricted open-shell Hartree-Fock (ROHF) implementation that is able to generate self-consistent field (SCF) wave functions for an arbitrary configuration state function (CSF). These CSFs can contain an arbitrary number of unpaired electrons in arbitrary spin-couplings. The resulting method is named CSF-ROHF. We demonstrate that starting from the ROHF energy expression, for example, the one given by Edwards and Zerner, it is possible to obtain the values of the ROHF vector-coupling coefficients by setting up an open-shell for each group of consecutive parallel-coupled spins dictated by the unique spin-coupling pattern of any given CSF. To achieve this important and nontrivial goal, we employ the machinery of the iterative configuration expansion configuration interaction (ICE-CI) method, which is able to tackle general CI problems on the basis of spin-adapted CSFs. This development allows for the efficient generation of SCF spin-eigenfunctions for systems with complex spin-coupling patterns, such as polymetallic chains and metal clusters, while maintaining SCF scaling with system size (quadratic or less, depending on the specific algorithm and approximations chosen).

Natural Orbitals and Targeted Non-Orthogonal Orbital Sets for Atomic Hyperfine Structure Multiconfiguration Calculations

Hyperfine structure constants have many applications, but are often hard to calculate accurately due to large and canceling contributions from different terms of the hyperfine interaction operator, and also from different closed and spherically symmetric core subshells that break up due to electron correlation effects. In multiconfiguration calculations, the wave functions are expanded in terms of configuration state functions (CSFs) built from sets of one-electron orbitals. The orbital sets are typically enlarged within the layer-by-layer approach. The calculations are energy-driven, and orbitals in each new layer of correlation orbitals are spatially localized in regions where the weighted total energy decreases the most, overlapping and breaking up different closed core subshells in an irregular pattern. As a result, hyperfine structure constants, computed as expectation values of the hyperfine operators, often show irregular or oscillating convergence patterns. Large orbital sets, and associated large CSF expansions, are needed to obtain converged values of the hyperfine structure constants. We analyze the situation for the states of the {2s22p3,2s22p23p,2s22p24p} odd and {2s22p23s,2s2p4,2s22p24s,2s22p23d} even configurations in N I, and show that the convergence with respect to the increasing sets of orbitals is radically improved by introducing separately optimized orbital sets targeted for describing the spin- and orbital-polarization effects of the 1s and 2s core subshells that are merged with, and orthogonalized against, the ordinary energy-optimized orbitals. In the layer-by-layer approach, the spectroscopic orbitals are kept frozen from the initial calculation and are not allowed to relax in response to the introduced layers of correlation orbitals. To compensate for this lack of variational freedom, the orbitals are transformed to natural orbitals prior to the final calculation based on single and double substitutions from an increased multireference set. The use of natural orbitals has an important impact on the states of the 2s22p23s configuration, bringing the corresponding hyperfine interaction constants in closer agreement with experiment. Relying on recent progress in methodology, the multiconfiguration calculations are based on configuration state function generators, cutting down the time for spin-angular integration by factors of up to 50, compared to ordinary calculations.

Second-order Rayleigh–Schrödinger perturbation theory for the GRASP2018 package: Core–valence correlations

The General Relativistic Atomic Structure package [GRASP2018, C. Froese Fischer, G. Gaigalas, P. Jönsson, and J. Bieroń, Comput. Phys. Commun. (2019), DOI: 10.1016/j.cpc.2018.10.032] is based on multiconfiguration Dirac– Hartree–Fock and relativistic configuration interaction (RCI) methods for energy structure calculations. The atomic state function used in the program is built from the set of configuration state functions (CSFs). The valence–valence, core–valence and core–core correlations are explicitly included through expansions over CSFs in RCI. We present a combination of RCI and the stationary second-order Rayleigh–Schrödinger many-body perturbation theory in an irreducible tensorial form to account for electron core–valence correlations when an atom or ion has any number of valence electrons. This newly developed method, which offers two ways of use, allows a significant reduction of the CSF space for complex atoms and ions. We also demonstrate how the method and program works for the energy structure calculation of Cl III ion.

Spin-Adapted Externally Contracted Multireference Configuration Interaction Method Based on Selected Reference Configurations.

As one kind of approximation of the full configuration interaction solution, the selected configuration interaction (sCI) methods have been shown to be valuable for large active spaces. However, the inclusion of dynamic correlation beyond large active spaces is necessary for more quantitative results. Since the sCI wave function can provide a compact reference for multireference methods, previously, we proposed an externally contracted multireference configuration interaction method using the sCI reference reconstructed from the density matrix renormalization group wave function [J. Chem. Theory Comput. 2018, 14, 4747-4755]. The DMRG2sCI-EC-MRCI method is promising for dealing with more than 30 active orbitals and large basis sets. However, it suffers from two drawbacks: spin contamination and low efficiency when using Slater determinant bases. To solve these problems, in this work, we adopt configuration state function bases and introduce a new algorithm based on the hybrid of tree structure for convenient configuration space management and the graphical unitary group approach for efficient matrix element calculation. The test calculation of naphthalene shows that the spin-adapted version could achieve a speed-up of 6.0 compared with the previous version based on the Slater determinant. Examples of dinuclear copper(II) compound as well as Ln(III) and An(III) complexes show that the sCI-EC-MRCI can give quantitatively accurate results by including dynamic correlation over sCI for systems with large active spaces and basis sets.

Energies and radiative properties of 1s22l2l′ and 1s(2l)22l′ levels in tungsten ion, W70+

Multiconfiguration Dirac–Hartree–Fock (MCDHF) calculations of energy levels, oscillator strengths and transition rates of open K-shell transitions in beryllium-like tungsten ion, W[Formula: see text], are performed using the fully relativistic GRASP2018 code. For only principal quantum number [Formula: see text], energy levels of 1[Formula: see text] 2[Formula: see text], 1[Formula: see text]2[Formula: see text], 1[Formula: see text]2[Formula: see text], 1[Formula: see text]2p, 1[Formula: see text] and 1[Formula: see text] at different angular momenta and parity ([Formula: see text] are presented. Excitations of up to four electrons from the ground level, 1[Formula: see text] 2[Formula: see text], to excited levels are allowed. Correlations up to 5l orbitals are included in the configuration state function (CSF), [Formula: see text], p, d, f and g. The fine-structure calculations of E1, E2, M1 and M2 oscillator strengths and radiative rates between intercombination transitions in W[Formula: see text] ions are evaluated. The accuracy of these results is estimated using different ways, such as the comparison with previously published results and calculating the uncertainty of the transition data. The uncertainty of the transition probabilities showed low values for most of transitions.

Predicting the X-ray Absorption Spectrum of Ozone with Single Configuration State Functions.

X-ray absorption spectra (XAS) of biradicaloid species are often thought to represent a challenge to theoretical methods. This has led to the testing of recently developed multireference techniques on the XAS of ozone, but reproduction of the experimental spectral profile has proven difficult. We utilize a minimal model consisting of a single configuration state function (CSF) per excited state to model core-level excitations of ozone, with the orbitals of each CSF optimized using the restricted open-shell Kohn-Sham (ROKS) method. This protocol leads to semiquantitative agreement with experimental XAS. In fact, we find that low-lying core-hole excited states in biradicaloids can be approximated with individual CSFs, despite the presence of multireference character in the ground state. We also report that the 1s → π* and 1s → σ* transitions have quite distinct widths for O3. This reveals the importance of sampling over a representative range of geometries from the vibrational ground state for properly assessing the accuracy of electronic structure methods against experiments instead of the popular procedure of uniformly broadening stick spectra at the equilibrium geometry.

PyBEST: Improved functionality and enhanced performance

PyBEST: Improved functionality and enhanced performance

Modular Approach to Selected Configuration Interaction in an Arbitrary Spin Basis: Implementation and Comparison of Approaches.

A modular selected configuration interaction (SCI) code has been developed that is based on the existing Monte-Carlo configuration interaction code (MCCI). The modularity allows various selection protocols to be implemented with ease and allows for fair comparison between wave functions built via different criteria. We have initially implemented adaptations of existing SCI theories, which are based on either energy- or coefficient-driven selection schemes. These codes have been implemented not only in the basis of Slater determinants (SDs) but also in the basis of configuration state functions (CSFs) and extended to state-averaged regimes. This allows one to take advantage of the reduced dimensionality of the wave function in the CSF basis and also the guarantee of pure spin states. All SCI methods were found to be able to predict potential energy surfaces to high accuracy, producing compact wave functions, when compared to full configuration interaction (FCI) for a variety of bond-breaking potential energy surfaces. The compactness of the error-controlled adaptive configuration interaction approach, particularly in the CSF basis, was apparent with nonparallelity errors within chemical accuracy while containing as little as 0.02% of the FCI CSF space. The size-to-accuracy was also extended to FCI spaces approaching one billion configurations.

Local Oxidation States in {FeNO}6-8 Porphyrins: Insights from DMRG/CASSCF-CASPT2 Calculations.

A first DMRG/CASSCF-CASPT2 study of a series of paradigmatic {FeNO}6, {FeNO}7, and {FeNO}8 heme-nitrosyl complexes has led to substantial new insight as well as uncovered key shortcomings of the DFT approach. By virtue of its balanced treatment of static and dynamic correlation, the calculations have provided some of the most authoritative information available to date on the energetics of low- versus high-spin states of different classes of heme-nitrosyl complexes. Thus, the calculations indicate low doublet-quartet gaps of 1-4 kcal/mol for {FeNO}7 complexes and high singlet-triplet gaps of ≳20 kcal/mol for both {FeNO}6 and {FeNO}8 complexes. In contrast, DFT calculations yield widely divergent spin state gaps as a function of the exchange-correlation functional. DMRG-CASSCF calculations also help calibrate DFT spin densities for {FeNO}7 complexes, pointing to those obtained from classic pure functionals as the most accurate. The general picture appears to be that nearly all the spin density of Fe[P](NO) is localized on the Fe, while the axial ligand imidazole (ImH) in Fe[P](NO)(ImH) pushes a part of the spin density onto the NO moiety. An analysis of the DMRG-CASSCF wave function in terms of localized orbitals and of the resulting configuration state functions in terms of resonance forms with varying NO(π*) occupancies has allowed us to address the longstanding question of local oxidation states in heme-nitrosyl complexes. The analysis indicates NO(neutral) resonance forms [i.e., Fe(II)-NO0 and Fe(III)-NO0] as the major contributors to both {FeNO}6 and {FeNO}7 complexes. This finding is at variance with the common formulation of {FeNO}6 hemes as Fe(II)-NO+ species but is consonant with an Fe L-edge XAS analysis by Solomon and co-workers. For the {FeNO}8 complex {Fe[P](NO)}-, our analysis suggests a resonance hybrid description: Fe(I)-NO0 ↔ Fe(II)-NO-, in agreement with earlier DFT studies. Vibrational analyses of the compounds studied indicate an imperfect but fair correlation between the NO stretching frequency and NO(π*) occupancy, highlighting the usefulness of vibrational data as a preliminary indicator of the NO oxidation state.

The efficient conversion between linear wave function expansions and nonlinear graphically contracted function expansions

Several algorithms are examined for the transformation of a graphically contracted function (GCF) to a configuration state function (CSF) vector, denoted , and for the reverse transformation from CSF vector to GCF, denoted . The initial implementation of these algorithms is applied to a sequence of full-CI wave functions of similar molecules. The most efficient transformation is based on a depth-first search (DFS) and scales as where is an average facet count for the GCF wave function, is the dimension of the vector x, and is the dimension of the GCF basis. Two different transformations are proposed, one based on a pairwise recursive merge algorithm and the other based on a least-squares fitting algorithm. The recursive merge is a reliable algorithm, but our initial implementation does not display the expected scaling behaviour with respect to . The iterative least-squares fitting algorithm is based on an efficient DFS procedure to compute the error ε and its gradient g. However, the resulting iterative least-squares optimisation is slowly convergent, requiring excessive iterations to converge to a target error tolerance. For most of the calculations in this case study, the recursive merge procedure is found to be more efficient than the least-squares optimisation procedure.

Chemical Bonding in the C2 Molecule

Bonding in the C2 molecule is investigated with CAS(8,8) wave functions using canonical MOs. In a subsequent step, orthogonal atomic orbitals are constructed by localizing the CASSCF MOs on the two carbon atoms with an orthogonal transformation. This orbital transformation causes an orthogonal transformation of the configuration state functions (CSF) spanning the function space of the singlet ground state of C2. Instead of CSFs built from canonical MOs, one obtains CSFs of orthogonal deformed atomic orbitals (AO). This approach resembles the orthogonal valence bond (OVB) methods’ CSFs, which are very different from the conventional VB, based on non-orthogonal AOs. To become used to the different argumentation, the bonding situations in ethane (single bond), ethene (double bond), and the nitrogen molecule (triple bond) are also studied. The complex bonding situation in C2 is caused by the possibility to excite an electron with a spin flip from the doubly occupied 2s AO into the 2p subshell, and the resulting high-spin 5Su state of the carbon atom allows for a better reduction of the Pauli repulsion. However, the electron structure around the equilibrium distance does not allow one to say that C2 in its ground state has a double, or triple, or even a quadruple bond.

Fine-Tuning of Atomic Energies in Relativistic Multiconfiguration Calculations

Ab initio calculations sometimes do not reproduce the experimentally observed energy separations at a high enough accuracy. Fine-tuning of diagonal elements of the Hamiltonian matrix is a process which seeks to ensure that calculated energy separations of the states that mix are in agreement with experiment. The process gives more accurate measures of the mixing than can be obtained in ab initio calculations. Fine-tuning requires the Hamiltonian matrix to be diagonally dominant, which is generally not the case for calculations based on jj-coupled configuration state functions. We show that this problem can be circumvented by a method that transforms the Hamiltonian in jj-coupling to a Hamiltonian in LSJ-coupling for which fine-tuning applies. The fine-tuned matrix is then transformed back to a Hamiltonian in jj-coupling. The implementation of the method into the General Relativistic Atomic Structure Package is described and test runs to validate the program operations are reported. The new method is applied to the computation of the 2s21S0−2s2p1,3P1 transitions in C III and to the computation of Rydberg transitions in B I, for which the 2s2p22S1/2 perturber enters the 2s2ns2S1/2 series. Improved convergence patterns and results are found compared with ab initio calculations.

Measuring Electron Correlation: The Impact of Symmetry and Orbital Transformations.

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

State-Specific Configuration Interaction for Excited States.

We introduce and benchmark a systematically improvable route for excited-state calculations, labeled state-specific configuration interaction (ΔCI), which is a particular realization of multiconfigurational self-consistent field and multireference configuration interaction. Starting with a reference built from optimized configuration state functions, separate CI calculations are performed for each targeted state (hence, state-specific orbitals and determinants). Accounting for single and double excitations produces the ΔCISD model, which can be improved with second-order Epstein-Nesbet perturbation theory (ΔCISD+EN2) or a posteriori Davidson corrections (ΔCISD+Q). These models were gauged against a vast and diverse set of 294 reference excitation energies. We have found that ΔCI is significantly more accurate than standard ground-state-based CI, whereas close performances were found between ΔCISD and EOM-CC2 and between ΔCISD+EN2 and EOM-CCSD. For larger systems, ΔCISD+Q delivers more accurate results than EOM-CC2 and EOM-CCSD. The ΔCI route can handle challenging multireference problems, singly and doubly excited states, from closed- and open-shell species, with overall comparable accuracy and thus represents a promising alternative to more established methodologies. In its current form, however, it is reliable only for relatively low-lying excited states.

Resolution of Electronic States in Heisenberg Cluster Models within the Unitary Group Approach.

In this work ground and excited electronic states of Heisenberg cluster models, in the form of configuration interaction many-body wave functions, are characterized within the spin-adapted Graphical Unitary Group Approach framework, and relying on a novel combined unitary and symmetric group approach. Finite-size cluster models of well-defined point-group symmetry and of general local-spin are presented, including J1-J2 triangular and tetrahedral clusters, which are often used to describe magnetic interactions in biological and biomimetic polynuclear transition metal clusters with unique catalytic activity, such as nitrogen fixation and photosynthesis. We show that a unique block-diagonal structure of the underlying Hamiltonian matrix in the spin-adapted basis emerges when an optimal lattice site ordering is chosen that reflects the internal symmetries of the model investigated. The block-diagonal structure is bound to the commutation relations between cumulative spin operators and the Hamiltonian operator, that in turn depend on the geometry of the cluster investigated. The many-body basis transformation, in the form of the orbital/site reordering, exposes such commutation relations. These commutation relations represent a rigorous and formal demonstration of the block-diagonal structure in Hamiltonian matrices and the compression of the corresponding spin-adapted many-body wave functions. As a direct consequence of the block-diagonal structure of the Hamiltonian matrix, it is possible to selectively optimize electronic excited states without the overhead of calculating the lower-energy states by simply relying on the initial ansatz for the targeted wave function. Additionally, more compact many-body wave functions are obtained. In extreme cases, electronic states are precisely described by a single configuration state function, despite the curse of dimensionality of the corresponding Hilbert space. These findings are crucial in the electronic structure theory framework, for they offer a conceptual route toward wave functions of reduced multireference character, that can be optimized more easily by approximated eigensolvers and are of more facile physical interpretation. They open the way to study larger ab initio and model Hamiltonians of increasingly larger number of correlated electrons, while keeping the computational costs at their lowest. In particular, these elements will expand the potential of electronic structure methods in understanding magnetic interactions in exchange-coupled polynuclear transition metal clusters.

Independently Optimized Orbital Sets in GRASP—The Case of Hyperfine Structure in Li I

In multiconfiguration Dirac–Hartree–Fock (MCDHF) calculations, there is a strong coupling between the localization of the orbital set and the configuration state function (CSF) expansion used to determine it. Furthermore, it is well known that an orbital set resulting from calculations, including CSFs describing core–core correlation and other effects, which aims to lower the weighted energies of a number of targeted states as much as possible, may be inadequate for building CSFs that account for correlation effects that are energetically unimportant but decisive for computed properties, e.g., hyperfine structures or transition rates. This inadequacy can be traced in irregular or oscillating convergence patterns of the computed properties as functions of the increasing orbital set. In order to alleviate the above problems, we propose a procedure in which the orbital set is obtained by merging several separately optimized, and mutually non-orthogonal, orbital sets. This computational strategy preserves the advantages of capturing electron correlation on the total energy through the variational MCDHF method and allows to target efficiently the correlation effects on the considered property. The orbital sets that are merged are successively orthogonalized against each other to retain orthonormality. The merged orbital set is used to build CSFs that efficiently lower the energy and also adequately account for the correlation effects that are important for the property. We apply the procedure to compute the hyperfine structure constants for the 1s22s2S1/2 and 1s22p2P1/2,3/2o states in 7Li and show that it leads to considerably improved convergence patterns with respect to the increasing orbital set compared to standard calculations based on a single orbital set, energy-optimized in the variational procedure. The perspectives of the new procedure are discussed in a broader context in the summary.

The rate constant for the automerization of vinyl radical: A theoretical approach using the Complete Active Space Self Consistent Field method

The Complete Active Space Self Consistent Field method (CASSCF) was utilized to investigate the rate constant for the automerization of vinyl radical. Using an enlarged active space including all valence electrons, thus producing a total number of 213,444 Configuration State Functions (CSFs) the rate constant for the automerization is of the order of 109 s−1 or 1010 s−1 in case of using the vibrational energy splitting tunneling method with the corresponding energy barrier height, involving the zero point vibrational energy corrections, varying between 4.69 and 4.71 kcal/mol (1639.69–1647.34 cm−1). Furthermore depending on the basis set, the ground state vibrational proton tunneling splitting, ΔEsplit, varies between 0.27 and 0.29 cm−1.