Spin Purification in Full-CI Quantum Monte Carlo via a First-Order Penalty Approach.

In this article, we demonstrate that a first-order spin penalty scheme can be efficiently applied to the Slater determinant based Full-CI Quantum Monte Carlo (FCIQMC) algorithm, as a practical route toward spin purification. Two crucial applications are presented to demonstrate the validity and robustness of this scheme: the 1Δg ← 3Σg vertical excitation in O2 and key spin gaps in a [Mn3(IV)O4] cluster. In the absence of a robust spin adaptation/purification technique, both applications would be unattainable by Slater determinant based ground state methods, with any starting wave function collapsing into the higher-spin ground state during the optimization. This strategy can be coupled to other algorithms that use the Slater determinant based FCIQMC algorithm as configuration interaction eigensolver, including the Stochastic Generalized Active Space, the similarity-transformed FCIQMC, the tailored-CC, and second-order perturbation theory approaches. Moreover, in contrast to the GUGA-FCIQMC technique, this strategy features both spin projection and total spin adaptation, making it appealing when solving anisotropic Hamiltonians. It also provides spin-resolved reduced density matrices, important for the investigation of spin-dependent properties in polynuclear transition metal clusters, such as the hyperfine-coupling constants.

Quenched Lewis Acidity: Studies on the Medium Dependent Fluorescence of Zinc(II) Complexes.

Three new zinc(II) coordination units [Zn(1-3)] based on planar-directing tetradentate Schiff base-like ligands H2 (1-3) were synthesized. Their solid-state structures were investigated by single crystal X-ray diffraction, showing the tendency to overcome the square-planar coordination sphere by axial ligation. Affinity in solution towards axial ligation has been tested by extended spectroscopic studies, both in the absorption and emission mode. The electronic spectrum of the pyridine complex [Zn(1)(py)] has been characterized by MC-PDFT to validate the results of extended TD-DFT studies. Green emission of non-emissive solutions of [Zn(1-3)] in chloroform could be switched on in the presence of potent Lewis-bases. While interpretation in terms of an equilibrium of stacked/non-fluorescent and destacked/fluorescent species is in line with precedents from literature, the sensitivity of [Zn(1-3)] was greatly reduced. Results of a computation-based structure search allow to trace the hidden Lewis acidity of [Zn(1-3)] to a new stacking motif, resulting in a strongly enhanced stability of the dimers.

Modeling magnetic interactions in high-valent trinuclear [Mn3(IV)O4]4+ complexes through highly compressed multi-configurational wave functions.

In this work we apply a quantum chemical framework, recently designed in our laboratories, to rationalize the low-energy electronic spectrum and the magnetic properties of an homo-valent trinuclear [Mn3(IV)O4]4+ model of the oxygen-evolving center in photosystem II. The method is based on chemically motivated molecular orbital unitary transformations, and the optimization of spin-adapted many-body wave functions, both for ground- and excited-states, in the transformed MO basis. In this basis, the configuration interaction Hamiltonian matrix of exchange-coupled multi-center clusters is extremely sparse and characterized by a unique block diagonal structure. This property leads to highly compressed wave functions (oligo- or single-reference) and crucially enables state-specific optimizations. This work is the first showing that compression and selective targeting of ground- and excited-states wave functions is possible for systems with three magnetic centers that are not exactly half-filled, and that potentially exhibit frustrated spin interactions. The reduced multi-reference character of the wave function greatly simplifies the interpretation of the ground- and excited-state electronic structures, and provides a route for the direct rationalization of magnetic interactions in these compounds, often considered a challenge in polynuclear transition-metal chemistry. In this study, strong electron correlation effects have explicitly been described by conventional and stochastic multiconfigurational methodologies, while dynamic correlation effects have been accounted for by multiconfigurational second order perturbation theory, CASPT2. Ab initio results for the [Mn3(IV)O4]4+ system have been mapped to a three-site Heisenberg model with two magnetic coupling constants. The magnetic coupling constants and the temperature dependence of the effective magnetic moment predicted by the ab initio calculations are in good agreement with the available experimental data, and confirm the antiferromagnetic interaction among the three magnetic centers, while providing a simple and rigorous description of the noncollinearity of the local spins, that characterize most of the low-energy states for this system.

Resolution of Low-Energy States in Spin-Exchange Transition-Metal Clusters: Case Study of Singlet States in [Fe(III)4S4] Cubanes.

Polynuclear transition-metal (PNTM) clusters owe their catalytic activity to numerous energetically low-lying spin states and stable oxidation states. The characterization of their electronic structure represents one of the greatest challenges of modern chemistry. We propose a theoretical framework that enables the resolution of targeted electronic states with ease and apply it to two [Fe(III)4S4] cubanes. Through direct access to their many-body wave functions, we identify important correlation mechanisms and their interplay with the geometrical distortions observed in these clusters, which are core properties in understanding their catalytic activity. The simulated magnetic coupling constants predicted by our strategy allow us to make qualitative connections between spin interactions and geometrical distortions, demonstrating its predictive power. Moreover, despite its simplicity, the strategy provides magnetic coupling constants in good agreement with the available experimental ones. The complexes are intrinsically frustrated anti-ferromagnets, and the obtained spin structures together with the geometrical distortions represent two possible ways to release spin frustration (spin-driven Jahn-Teller distortion). Our paradigm provides a simple, yet rigorous, route to uncover the electronic structure of PNTM clusters and may be applied to a wide variety of such clusters.

Multiconfiguration Pair-Density Functional Theory: A New Way To Treat Strongly Correlated Systems.

The electronic energy of a system provides the Born-Oppenheimer potential energy for internuclear motion and thus determines molecular structure and spectra, bond energies, conformational energies, reaction barrier heights, and vibrational frequencies. The development of more efficient and more accurate ways to calculate the electronic energy of systems with inherently multiconfigurational electronic structure is essential for many applications, including transition metal and actinide chemistry, systems with partially broken bonds, many transition states, and most electronically excited states. Inherently multiconfigurational systems are called strongly correlated systems or multireference systems, where the latter name refers to the need for using more than one ("multiple") configuration state function to provide a good zero-order reference wave function. This Account describes multiconfiguration pair-density functional theory (MC-PDFT), which was developed as a way to combine the advantages of wave function theory (WFT) and density functional theory (DFT) to provide a better treatment of strongly correlated systems. First we review background material: the widely used Kohn-Sham DFT (which uses only a single Slater determinant as reference wave function), multiconfiguration WFT methods that treat inherently multiconfigurational systems based on an active space, and previous attempts to combine multiconfiguration WFT with DFT. Then we review the formulation of MC-PDFT. It is a generalization of Kohn-Sham DFT in that the electron kinetic energy and classical electrostatic energy are calculated from a reference wave function, while the rest of the energy is obtained from a density functional. However, there are two main differences with respent to Kohn-Sham DFT: (i) The reference wave function is multiconfigurational rather than being a single Slater determinant. (ii) The density functional is a function of the total density and the on-top pair density rather than being a function of the spin-up and spin-down densities. In work carried out so far, the multiconfigurational wave function is a multiconfiguration self-consistent-field wave function. The new formulation has the advantage that the reference wave function has the correct spatial and spin symmetry and can describe bond dissociation (of both single and multiple bonds) and electronic excitations in a formally and physically correct way. We then review the formulation of density functionals in terms of the on-top pair density. Finally we review successful applications of the theory to bond energies and bond dissociation potential energy curves of main-group and transition metal bonds, to barrier heights (including pericyclic reactions), to proton affinities, to the hydrogen bond energy of water dimer, to ground- and excited-state charge transfer, to valence and Rydberg excitations of molecules, and to singlet-triplet splittings of radicals. We find that that MC-PDFT can give accurate results not only with complete-active-space multiconfiguration wave functions but also with generalized-active-space multiconfiguration wave functions, which are practical for larger numbers of active electrons and active orbitals than are complete-active-space wave functions. The separated-pair approximation, which is a special case of generalized active space self-consistent-field theory, is especially promising. MC-PDFT, because it requires much less computer time and storage than pure WFT methods, has the potential to open larger and more complex strongly correlated systems to accurate simulation.

Understanding the Mechanism Stabilizing Intermediate Spin States in Fe(II)-Porphyrin.

Spin fluctuations in Fe(II)-porphyrins are at the heart of heme-proteins functionality. Despite significant progress in porphyrin chemistry, the mechanisms that rule spin state stabilization remain elusive. Here, it is demonstrated by using multiconfigurational quantum chemical approaches, including the novel Stochastic-CASSCF method, that electron delocalization between the metal center and the π system of the macrocycle differentially stabilizes the triplet spin states over the quintet. This delocalization takes place via charge-transfer excitations, involving the π system of the macrocycle and the out-of-plane iron d orbitals, key linking orbitals between metal and macrocycle. Through a correlated breathing mechanism the 3d electrons can make transitions toward the π orbitals of the macrocycle. This guarantees a strong coupling between the on-site radial correlation on the metal and electron delocalization. Opposite-spin 3d electrons of the triplet can effectively reduce electron repulsion in this manner. Constraining the out-of-plane orbitals from breathing hinders delocalization and reverses the spin ordering. The breathing mechanism is made effective by strong electron correlation effects in the π system of the macrocycle. Reducing the correlation treatment on the macrocycle to second-order only also reverses the spin ordering. High order (beyond second-order) correlation on the macrocycle reduces the energetic cost of the additional electron to a sufficient extent to stabilize the triplet state. Our results find a qualitative analogue in six resonance structures involving the metal center in the Fe2+ and Fe3+ oxidation states.

Molcas 8: New capabilities for multiconfigurational quantum chemical calculations across the periodic table.

In this report, we summarize and describe the recent unique updates and additions to the Molcas quantum chemistry program suite as contained in release version 8. These updates include natural and spin orbitals for studies of magnetic properties, local and linear scaling methods for the Douglas-Kroll-Hess transformation, the generalized active space concept in MCSCF methods, a combination of multiconfigurational wave functions with density functional theory in the MC-PDFT method, additional methods for computation of magnetic properties, methods for diabatization, analytical gradients of state average complete active space SCF in association with density fitting, methods for constrained fragment optimization, large-scale parallel multireference configuration interaction including analytic gradients via the interface to the Columbus package, and approximations of the CASPT2 method to be used for computations of large systems. In addition, the report includes the description of a computational machinery for nonlinear optical spectroscopy through an interface to the QM/MM package Cobramm. Further, a module to run molecular dynamics simulations is added, two surface hopping algorithms are included to enable nonadiabatic calculations, and the DQ method for diabatization is added. Finally, we report on the subject of improvements with respects to alternative file options and parallelization.

Comment on "Fe2: As simple as a Herculean labour. Neutral (Fe2), cationic (Fe2(+)), and anionic (Fe2(-)) species" [J. Chem. Phys. 142, 244304 (2015)].

A recent paper on Fe2 [A. Kalemos, J. Chem. Phys. 142, 244304 (2015)] critiqued our previous work on the system [Hoyer et al., J. Chem. Phys. 141, 204309 (2014)]. In this comment, we explain the nature of our previously reported potential energy curve for Fe2 and we discuss our computed properties for Fe2. Additionally, we fix a labeling error that was present in our previous work, although this error is unrelated to the main point of discussion.

A two-coordinate manganese(0) complex with an unsupported Mn-Mg bond: allowing access to low coordinate homo- and heterobimetallic compounds.

This study details the synthesis and characterization of an unprecedented two-coordinate, high-spin manganese(0) complex that incorporates an unsupported Mn-Mg bond, viz. L(†)MnMg((Mes)Nacnac) (L(†) = -N(Ar(†))(SiPr(i)3), Ar(†) = C6H2{C(H)Ph2}2Pr(i)-2,6,4; (Mes)Nacnac = [(MesNCMe)2CH](-); Mes = mesityl). This compound has been utilized as an "inorganic Grignard reagent" in the preparation of the first two-coordinate manganese(I) dimer, L(†)MnMnL* (L* = -N(Ar*)(SiMe3), Ar* = C6H2{C(H)Ph2}2Me-2,6,4), and the related mixed valence, bis(amido)-hetereobimetallic complex, Mn(II)(µ-L(†))(µ-L*)Cr(0). It is also shown to act as a two-electron reducing agent in reactions with unsaturated substrates.

Oxidative stretching of metal-metal bonds to their limits.

Oxidation of quadruply bonded Cr2(dpa)4, Mo2(dpa)4, MoW(dpa)4, and W2(dpa)4 (dpa = 2,2'-dipyridylamido) with 2 equiv of silver(I) triflate or ferrocenium triflate results in the formation of the two-electron-oxidized products [Cr2(dpa)4](2+) (1), [Mo2(dpa)4](2+) (2), [MoW(dpa)4](2+) (3), and [W2(dpa)4](2+) (4). Additional two-electron oxidation and oxygen atom transfer by m-chloroperoxybenzoic acid results in the formation of the corresponding metal-oxo compounds [Mo2O(dpa)4](2+) (5), [WMoO(dpa)4](2+) (6), and [W2O(dpa)4](2+) (7), which feature an unusual linear M···M≡O structure. Crystallographic studies of the two-electron-oxidized products 2, 3, and 4, which have the appropriate number of orbitals and electrons to form metal-metal triple bonds, show bond distances much longer (by >0.5 Å) than those in established triply bonded compounds, but these compounds are nonetheless diamagnetic. In contrast, the Cr-Cr bond is completely severed in 1, and the resulting two isolated Cr(3+) magnetic centers couple antiferromagnetically with J/kB= -108(3) K [-75(2) cm(-1)], as determined by modeling of the temperature dependence of the magnetic susceptibility. Density functional theory (DFT) and multiconfigurational methods (CASSCF/CASPT2) provide support for "stretched" and weak metal-metal triple bonds in 2, 3, and 4. The metal-metal distances in the metal-oxo compounds 5, 6, and 7 are elongated beyond the single-bond covalent radii of the metal atoms. DFT and CASSCF/CASPT2 calculations suggest that the metal atoms have minimal interaction; the electronic structure of these complexes is used to rationalize their multielectron redox reactivity.

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 Slocal>12 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.

Magnetic Interactions in a [Co(II)3Er(III)(OR)4] Model Cubane through Forefront Multiconfigurational Methods.

Strong electron correlation effects are one of the major challenges in modern quantum chemistry. Polynuclear transition metal clusters are peculiar examples of systems featuring such forms of electron correlation. Multireference strategies, often based on but not limited to the concept of complete active space, are adopted to accurately account for strong electron correlation and to resolve their complex electronic structures. However, transition metal clusters already containing four magnetic centers with multiple unpaired electrons make conventional active space based strategies prohibitively expensive, due to their unfavorable scaling with the size of the active space. In this work, forefront techniques, such as density matrix renormalization group (DMRG), full configuration interaction quantum Monte Carlo (FCIQMC), and multiconfiguration pair-density functional theory (MCPDFT), are employed to overcome the computational limitation of conventional multireference approaches and to accurately investigate the magnetic interactions taking place in a [Co(II)3Er(III)(OR)4] (chemical formula [Co(II)3Er(III)(hmp)4(µ2-OAc)2(OH)3(H2O)], hmp = 2-(hydroxymethyl)-pyridine) model cubane water oxidation catalyst. Complete active spaces with up to 56 electrons in 56 orbitals have been constructed for the seven energetically lowest different spin states. Relative energies, local spin, and spin-spin correlation values are reported and provide crucial insights on the spin interactions for this model system, pivotal in the rationalization of the catalytic activity of this system in the water-splitting reaction. A ferromagnetic ground state is found with a very small, ∼50 cm-1, highest-to-lowest spin gap. Moreover, for the energetically lowest states, S = 3-6, the three Co(II) sites exhibit parallel aligned spins, and for the lower states, S = 0-2, two Co(II) sites retain strong parallel spin alignment.

The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry.

The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.

Assessing metal-metal multiple bonds in Cr-Cr, Mo-Mo, and W-W compounds and a hypothetical U-U compound: a quantum chemical study comparing DFT and multireference methods.

To gain insights into the trends in metal-metal multiple bonding among the Group 6 elements, density functional theory has been employed in combination with multiconfigurational methods (CASSCF and CASPT2) to investigate a selection of bimetallic, multiply bonded compounds. For the compound [Ar-MM-Ar] (Ar=2,6-(C(6)H(5))(2)-C(6)H(3), M=Cr, Mo, W) the effect of the Ar ligand on the M(2) core has been compared with the analogous [Ph-MM-Ph] (Ph=phenyl, M=Cr, Mo, W) compounds. A set of [M(2)(dpa)(4)] (dpa=2,2'-dipyridylamide, M=Cr, Mo, W, U) compounds has also been investigated. All of the compounds studied here show important multiconfigurational behavior. For the Mo(2) and W(2) compounds, the σ(2)π(4)δ(2) configuration dominates the ground-state wavefunction, contributing at least 75%. The Cr(2) compounds show a more nuanced electronic structure, with many configurations contributing to the ground state. For the Cr, Mo, and W compounds the electronic absorption spectra have been studied, combining density functional theory and multireference methods to make absorption feature assignments. In all cases, the main features observed in the visible spectra may be assigned as charge-transfer bands. For all compounds investigated the Mayer bond order (MBO) and the effective bond order (EBO) were calculated by density functional theory and CASSCF methods, respectively. The MBO and EBO values share a similar trend toward higher values at shorter normalized metal-metal bond lengths.

Computational insights into uranium complexes supported by redox-active α-diimine ligands.

The electronic structures of two uranium compounds supported by redox-active α-diimine ligands, ((Mes)DAB(Me))(2)U(THF) (1) and Cp(2)U((Mes)DAB(Me)) (2) ((Mes)DAB(Me) = [ArN═C(Me)C(Me)═NAr]; Ar = 2,4,6-trimethylphenyl (Mes)), have been investigated using both density functional theory and multiconfigurational self-consistent field methods. Results from these studies have established that both uranium centers are tetravalent, that the ligands are reduced by two electrons, and that the ground states of these molecules are triplets. Energetically low-lying singlet states are accessible, and some transitions to these states are visible in the electronic absorption spectrum.

FCIQMC-Tailored Distinguishable Cluster Approach: Open-Shell Systems.

A recently proposed tailored approach based on the distinguishable cluster method and the stochastic FCI solver, FCIQMC [J. Chem. Theory Comput. 2020, 16, 5621], is extended to open-shell molecular systems. The method is employed to calculate spin gaps of various Fe(II) complexes, including a Fe(II) porphyrin model system. Both distinguishable cluster and fully relaxed CASSCF natural orbitals were used in this work as reference for the subsequent tailored distinguishable cluster calculations. The distinguishable cluster natural orbitals occupation numbers were also used as an aid to the selection of the active space. The effect of the active space sizes and of the explicit correlation correction (F12) onto the predicted spin gaps is investigated. The tailored distinguishable cluster with singles and doubles yields consistently more accurate results compared to the tailored coupled cluster with singles and doubles.

Stochastic Generalized Active Space Self-Consistent Field: Theory and Application.

An algorithm to perform stochastic generalized active space calculations, Stochastic-GAS, is presented, that uses the Slater determinant based FCIQMC algorithm as configuration interaction eigensolver. Stochastic-GAS allows the construction and stochastic optimization of preselected truncated configuration interaction wave functions, either to reduce the computational costs of large active space wave function optimizations, or to probe the role of specific electron correlation pathways. As for the conventional GAS procedure, the preselection of the truncated wave function is based on the selection of multiple active subspaces while imposing restrictions on the interspace excitations. Both local and cumulative minimum and maximum occupation number constraints are supported by Stochastic-GAS. The occupation number constraints are efficiently encoded in precomputed probability distributions, using the precomputed heat bath algorithm, which removes nearly all runtime overhead of GAS. This strategy effectively allows the FCIQMC dynamics to a priori exclude electronic configurations that are not allowed by GAS restrictions. Stochastic-GAS reduced density matrices are stochastically sampled, allowing orbital relaxations via Stochastic-GASSCF, and direct evaluation of properties that can be extracted from density matrices, such as the spin expectation value. Three test case applications have been chosen to demonstrate the flexibility of Stochastic-GAS: (a) the Stochastic-GASSCF [5·(6, 6)] optimization of a stack of five benzene molecules, that shows the applicability of Stochastic-GAS toward fragment-based chemical systems; (b) an uncontracted stochastic MRCISD calculation that correlates 96 electrons and 159 molecular orbitals, and uses a large (32, 34) active space reference wave function for an Fe(II)-porphyrin model system, showing how GAS can be applied to systematically recover dynamic electron correlation, and how in the specific case of the Fe(II)-porphyrin dynamic correlation further differentially stabilizes the 3Eg over the 5A1g spin state; (c) the study of an Fe4S4 cluster's spin-ladder energetics via highly truncated stochastic-GAS [4·(5, 5)] wave functions, where we show how GAS can be applied to understand the competing spin-exchange and charge-transfer correlating mechanisms in stabilizing different spin-states.

Strong correlation treated via effective hamiltonians and perturbation theory.

We propose a new approach to determine a suitable zeroth-order wavefunction for multiconfigurational perturbation theory. The same ansatz as in complete active space (CAS) wavefunction optimization is used but it is split in two parts, a principal space (A) and a much larger extended space (B). Löwdin's partitioning technique is employed to map the initial eigenvalue problem to a dimensionality equal to that of (A) only. Combined with a simplified expression for the (B) portion of the wavefunction, we are able to drastically reduce the storage and computational demands of the wavefunction optimization. This scheme is used to produce reference wavefunctions and energies for subsequent second-order perturbation theory (PT2) corrections. Releasing the constraint of computing the exact CAS energy and wavefunction prior to the PT2 treatment introduces a nonstandard paradigm for multiconfigurational methods. Based on the results of test calculations, we argue that principal parts with only few percents of the total number of CAS configurations could provide final multiconfigurational PT2 energies of the same accuracy as in the standard paradigm. In the future, algorithmic improvements for this scheme will bring into reach active spaces much beyond the present limit of CAS-based methods, therefore allowing for accurate studies of systems featuring strong correlation.

The generalized active space concept in multiconfigurational self-consistent field methods.

A multiconfigurational self-consistent field method based on the concept of generalized active space (GAS) is presented. GAS wave functions are obtained by defining an arbitrary number of active spaces with arbitrary occupation constraints. By a suitable choice of the GAS spaces, numerous ineffective configurations present in a large complete active space (CAS) can be removed, while keeping the important ones in the CI space. As a consequence, the GAS self-consistent field approach retains the accuracy of the CAS self-consistent field (CASSCF) ansatz and, at the same time, can deal with larger active spaces, which would be unaffordable at the CASSCF level. Test calculations on the Gd atom, Gd(2) molecule, and oxoMn(salen) complex are presented. They show that GAS wave functions achieve the same accuracy as CAS wave functions on systems that would be prohibitive at the CAS level.

Spin-Pure Stochastic-CASSCF via GUGA-FCIQMC Applied to Iron-Sulfur Clusters.

In this work, we demonstrate how to efficiently compute the one- and two-body reduced density matrices within the spin-adapted full configuration interaction quantum Monte Carlo (FCIQMC) method, which is based on the graphical unitary group approach (GUGA). This allows us to use GUGA-FCIQMC as a spin-pure configuration interaction (CI) eigensolver within the complete active space self-consistent field (CASSCF) procedure and hence to stochastically treat active spaces far larger than conventional CI solvers while variationally relaxing orbitals for specific spin-pure states. We apply the method to investigate the spin ladder in iron-sulfur dimer and tetramer model systems. We demonstrate the importance of the orbital relaxation by comparing the Heisenberg model magnetic coupling parameters from the CASSCF procedure to those from a CI-only (CASCI) procedure based on restricted open-shell Hartree-Fock orbitals. We show that the orbital relaxation differentially stabilizes the lower-spin states, thus enlarging the coupling parameters with respect to the values predicted by ignoring orbital relaxation effects. Moreover, we find that, while CASCI results are well fit by a simple bilinear Heisenberg Hamiltonian, the CASSCF eigenvalues exhibit deviations that necessitate the inclusion of biquadratic terms in the model Hamiltonian.