Search

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 O 2 and key spin gaps in a [Mn 3 (IV) O 4 ] 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.

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) 4 S 4 ] 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.

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 Fe 2+ and Fe 3+ oxidation states.