Calculation of spin–orbit couplings using RASCI spinless one-particle density matrices: Theory and applications

Coupled-cluster treatment of complex open-shell systems: the case of single-molecule magnets

We investigate the reliability of two cost-effective coupled-cluster methods for computing spin-state energetics and spin-related properties of a set of open-shell transition-metal complexes. Specifically, we employ the second-order approximate coupled-cluster singles and doubles (CC2) method and projection-based embedding that combines equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) with density functional theory (DFT). The performance of CC2 and EOM-CCSD-in-DFT is assessed against EOM-CCSD. The chosen test set includes two hexaaqua transition-metal complexes containing Fe(II) and Fe(III), and a large Co(II)-based single-molecule magnet with a non-aufbau ground state. We find that CC2 describes the excited states more accurately, reproducing EOM-CCSD excitation energies within 0.05 eV. However, EOM-CCSD-in-DFT excels in describing transition orbital angular momenta and spin-orbit couplings. Moreover, for the Co(II) molecular magnet, using EOM-CCSD-in-DFT eigenstates and spin-orbit couplings, we compute spin-reversal energy barriers, as well as temperature-dependent and field-dependent magnetizations and magnetic susceptibilities that closely match experimental values within spectroscopic accuracy. These results underscore the efficiency of CC2 in computing state energies of multi-configurational, open-shell systems and highlight the utility of the more cost-efficient EOM-CCSD-in-DFT for computing spin-orbit couplings and magnetic properties of complex and large molecular magnets.

Unveiling Hidden Shake-Up Features in the Uranyl M4-Edge Spectrum

The M4,5-edge high energy resolution X-ray absorption near-edge structure (HR-XANES) spectra of actinyls offer valuable insights into the electronic structure and bonding properties of heavy-element complexes. To conduct a comprehensive spectral analysis, it is essential to employ computational methods that accurately account for relativistic effects and electron correlation. In this work, we utilize variational relativistic multireference configurational interaction methods to compute and analyze the X-ray M4-edge absorption spectrum of uranyl. By employing these advanced computational techniques, we achieve excellent agreement between the calculated spectral features and experimental observations. Moreover, the calculations unveil significant shake-up features, which arise from the intricate interplay between strongly correlated 3d core-electron and ligand excitations. This research provides important theoretical insights into the spectral characteristics of heavy-element complexes. Furthermore, it establishes the foundation for utilizing M4,5-edge spectroscopy as a means to investigate the chemical activities of such complexes. By leveraging this technique, we can gain a deeper understanding of the bonding behavior and reactivity of heavy-element compounds.

Comparison of Variational and Perturbative Spin-Orbit Coupling within Two-Component CASSCF

The modeling of spin-orbit coupling (SOC) remains a challenge in computational chemistry due to the high computational cost. With the rising popularity of spin-driven processes and f-block metals in chemistry and materials science, it is incumbent on the community to develop accurate multiconfigurational SOC methods that scale to large systems and understand the limits of different treatments of SOC. Herein, we introduce an implementation of perturbative SOC in scalar-relativistic two-component CASSCF (srX2C-CASSCF-SO). Perspectives on the limitations and accuracy of srX2C-CASSCF-SO are presented via benchmark calculations.

Triplet-Singlet Emission of d-Block Metal Complexes Characterized by Spin-Orbit Natural Transition Orbitals

Spin-orbit natural transition orbital (SO-NTO) methodology, recently developed in our group for complete and restricted active space (CAS/RAS) wavefunction calculations, is applied to analyze triplet-to-singlet emission in transition metal complexes. The lowest-energy (longest-wavelength) spin-forbidden transitionT 1 → S 0 ${{{\rm T}}_{1}\to {{\rm S}}_{0}}$ is studied for for [Ir(pbt)2(acac)] and [Re(CO)4(pbt)] and the complexes [W(CO)4(bpy)] and [Mo(CO)4(bpy)]. For the latter complexes, spin-forbidden transitions from higher spin-triplet levels are additionally analyzed. SO-NTOs are compared with spin-free NTOs for the transitions under consideration. The major assignment of a spin-forbidden transition is obtained from the spin-free NTO analysis, while the source of intensity of the electronic transition is revealed by the SO-NTOs.

Galván I, Alavi A, Aleotti F, Aquilante F, Autschbach J, Avagliano D, Baiardi A, Bao JJ, Battaglia S, Birnoschi L, Blanco-González A, Bokarev SI, Broer R, Cacciari R, Calio PB, Carlson RK, Carvalho Couto R, Cerdán L, Chibotaru LF, Chilton NF, Church JR, Conti I, Coriani S, Cuéllar-Zuquin J, Daoud RE, Dattani N, Decleva P, de Graaf C, Delcey M, De Vico L, Dobrautz W, Dong SS, Feng R, Ferré N, Filatov(Gulak) M, Gagliardi L, Garavelli M, González L, Guan Y, Guo M, Hennefarth MR, Hermes MR, Hoyer CE, Huix-Rotllant M, Jaiswal VK, Kaiser A, Kaliakin DS, Khamesian M, King DS, Kochetov V, Krośnicki M, Kumaar AA, Larsson ED, Lehtola S, Lepetit MB, Lischka H, López Ríos P, Lundberg M, Ma D, Mai S, Marquetand P, Merritt ICD, Montorsi F, Mörchen M, Nenov A, Nguyen VHA, Nishimoto Y, Oakley MS, Olivucci M, Oppel M, Padula D, Pandharkar R, Phung QM, Plasser F, Raggi G, Rebolini E, Reiher M, Rivalta I, Roca-Sanjuán D, Romig T, Safari AA, Sánchez-Mansilla A, Sand AM, Schapiro I, Scott TR, Segarra-Martí J, Segatta F, Sergentu DC, Sharma P, Shepard R, Shu Y, Staab JK, Straatsma TP, Sørensen LK, Tenorio BNC, Truhlar DG, Ungur L, Vacher M, Veryazov V, Voß TA, Weser O, Wu D, Yang X, Yarkony D, Zhou C, Zobel JP, Lindh R. 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.

State-Interaction Approach for Evaluating g-Tensors within EOM-CC and RAS-CI Frameworks: Theory and Benchmarks

Among various techniques designed for studying open-shell species, electron paramagnetic resonance (EPR) spectroscopy plays an important role. The key quantity measured by EPR is the g-tensor, describing the coupling between an external magnetic field and molecular electronic spin. One theoretical framework for quantum chemistry calculations of g-tensors is based on response theory, which involves substantial developments that are specific to the underlying electronic structure models. A simplified and easier-to-implement approach is based on the state-interaction scheme, in which perturbation is included by considering a small number of states. We describe and benchmark the state-interaction approach using equation-of-motion coupled-cluster and restricted-active-space configuration interaction wave functions. The analysis confirms that this approach can deliver accurate results and highlights caveats of applying it, such as a choice of the reference state, convergence with respect to the number of states used in calculations, etc. The analysis also contributes toward a better understanding of challenges in calculations of higher-order properties using approximate wave functions.

State Interaction Linear Response Time-Dependent Density Functional Theory with Perturbative Spin-Orbit Coupling: Benchmark and Perspectives

Spin-orbit coupling (SOC) is an important driving force in photochemistry. In this work, we develop a perturbative spin-orbit coupling method within the linear response time-dependent density function theory framework (TDDFT-SO). A full state interaction scheme, including singlet-triplet and triplet-triplet coupling, is introduced to describe not only the coupling between the ground and excited states, but also between excited states with all couplings between spin microstates. In addition, expressions to compute spectral oscillator strengths are presented. Scalar relativity is included variationally using the second-order Douglas-Kroll-Hess Hamiltonian, and the TDDFT-SO method is validated against variational SOC relativistic methods for atomic, diatomic, and transition metal complexes to determine the range of applicability and potential limitations. To demonstrate the robustness of TDDFT-SO for large-scale chemical systems, the UV-Vis spectrum of Au₂₅(SR)₁₈ ^- is computed and compared to experiment. Perspectives on the limitation, accuracy, and capability of perturbative TDDFT-SO are presented via analyses of benchmark calculations. Additionally, an open-source Python software package (PyTDDFT-SO) is developed and released to interface with the Gaussian 16 quantum chemistry software package to perform this calculation.

Spin-Orbit Couplings of Open-Shell Systems with Restricted Active Space Configuration Interaction

In this work we perform electronic structure calculations to unravel the origin of spin-orbit couplings (SOCs) in open-shell molecules. For that, we select systems displaying di or polyradical character, e.g., trimethylene, and analyze the changes in the magnitude of SOC constants along molecular distortions of ethylene and in the presence of intermolecular interactions between open and closed-shell moieties in the O₂-C₂H₄ system. Calculations were performed by using nonrelativistic wave functions obtained with the restricted active space configuration interaction (RASCI) method, in conjunction with a recent implementation for the calculation of SOC based on the spin-orbit mean field approximation. Our results demonstrate the suitability of RASCI in the calculation of SOCs of open-shell systems, while providing a deep understanding of the relationship between couplings and the nature of the electronic states. Moreover, we introduce a new definition of the SOC constant for the study of molecular aggregates.

Two-Component Multireference Restricted Active Space Configuration Interaction for the Computation of L-Edge X-ray Absorption Spectra

X-ray absorption spectroscopy is a powerful probe of local electronic and nuclear structures, providing insights into chemical processes. The theoretical prediction and interpretation of metal L-edge X-ray absorption spectra are complicated by both relativistic effects, including spin-orbit coupling and the multiconfigurational nature of the states involved. This work details an exact two-component multireference restricted active space (RAS) configuration interaction scheme that uses an exact two-component state-averaged complete active space self-consistent-field method, which includes the spin-orbit coupling in a variational manner, for the accurate description of the electronic structure before using a RAS configuration interaction method to describe the core excited states of the X-ray spectrum. Benchmark calculations are presented for a series of iron-containing complexes, with results showing key features of the spectrum being reproduced, including ligand-to-metal charge transfer and shake-up excitations.

Spin-Orbit Natural Transition Orbitals and Spin-Forbidden Transitions

Natural transition orbitals (NTOs) are in widespread use for visualizing and analyzing electronic transitions. The present work introduces the analysis of formally spin-forbidden transitions with the help of complex-valued spin-orbit (SO) NTOs. The analysis specifically focuses on the components in such transitions that cause their intensity to be nonzero because of SO coupling. Transition properties such as transition dipole moments are partitioned into SO-NTO hole-particle pairs, such that contributions to the intensity from specific occupied and unoccupied orbitals are obtained. The method has been implemented within the restricted active space (RAS) self-consistent field wave function theory framework, with SO coupling treated by RAS state interaction. SO-NTOs have a broad range of potential applications, which is illustrated by the T₂-S₁ state mixing in pyrazine, spin-forbidden versus spin-allowed 4f-5d transitions in the Tb³⁺ ion, and the phosphorescence of tris(2-phenylpyridine) iridium [Ir(ppy)₃].

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Evaluation of two-particle properties within finite-temperature self-consistent one-particle Green's function methods: Theory and application to GW and GF2

One-particle Green's function methods can model molecular and solid spectra at zero or non-zero temperatures. One-particle Green's functions directly provide electronic energies and one-particle properties, such as dipole moment. However, the evaluation of two-particle properties, such as ⟨S²⟩ and ⟨N²⟩, can be challenging because they require a solution of the computationally expensive Bethe-Salpeter equation to find two-particle Green's functions. We demonstrate that the solution of the Bethe-Salpeter equation can be completely avoided. Applying the thermodynamic Hellmann-Feynman theorem to self-consistent one-particle Green's function methods, we derive expressions for two-particle density matrices in a general case and provide explicit expressions for GF2 and GW methods. Such density matrices can be decomposed into an antisymmetrized product of correlated one-electron density matrices and the two-particle electronic cumulant of the density matrix. Cumulant expressions reveal a deviation from ensemble representability for GW, explaining its known deficiencies. We analyze the temperature dependence of ⟨S²⟩ and ⟨N²⟩ for a set of small closed-shell systems. Interestingly, both GF2 and GW show a non-zero spin contamination and a non-zero fluctuation of the number of particles for closed-shell systems at the zero-temperature limit.

Equation-of-Motion Coupled-Cluster Protocol for Calculating Magnetic Properties: Theory and Applications to Single-Molecule Magnets

We present a new computational protocol for computing macroscopic magnetic properties of transition-metal complexes using the equation-of-motion coupled-cluster (EOM-CC) framework. The approach follows a two-step state-interaction scheme: we first compute zero-order states using nonrelativistic EOM-CC and then use these states to evaluate matrix elements of the spin-orbit and Zeeman operators. Diagonalization of the resulting Hamiltonian yields spin-orbit- and field-perturbed eigenstates. Temperature- and field-dependent magnetization and susceptibility are computed by numerical differentiation of the partition function. To compare with powder-sample experiments, these quantities are numerically averaged over field orientations. We applied this protocol to several single-molecule magnets (SMMs) with Fe(II) and Fe(III) in trigonal pyramidal, linear, and trigonal bipyramidal coordination environments. We described the underlying electronic structure by the electron-attachment (EOM-EA) and spin-flip (EOM-SF) variants of EOM-CC. The computed energy barriers for spin inversion, and macroscopic magnetization and susceptibility agree well with experimental data. Trends in magnetic anisotropy and spin-reversal energy barriers are explained in terms of a molecular orbital picture rigorously distilled from spinless transition density matrices between many-body states. The results illustrate excellent performances of EOM-CC in describing magnetic behavior of mononuclear transition-metal SMMs.