Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2207.02435 · The Journal of Chemical Physics

A comparison between the one- and two-step spin–orbit coupling approaches based on the ab initio density matrix renormalization group

Huanchen Zhai, Garnet Kin‐Lic Chan

Published Sep 29, 2022• 19 citations
arXivHTMLPDF

Abstract

The efficient and reliable treatment of both spin-orbit coupling (SOC) and electron correlation is essential for understanding f-element chemistry. We analyze two approaches to the problem: the one-step approach, where both effects are treated simultaneously, and the two-step state interaction approach. We report an implementation of the ab initio density matrix renormalization group with a one-step treatment of the SOC effect, which can be compared to prior two-step treatments on an equal footing. Using a dysprosium octahedral complex and bridged dimer as benchmark systems, we identify characteristics of problems where the one-step approach is beneficial for obtaining the low-energy spectrum.

DOI
10.1063/5.0107805

Related papers

Spin orbit coupling for molecular ab initio density matrix renormalization group calculations: Application to g-tensors

Michael Roemelt

Spin Orbit Coupling (SOC) is introduced to molecular ab initio density matrix renormalization group (DMRG) calculations. In the presented scheme, one first approximates the electronic ground state and a number of excited states of the Born-Oppenheimer (BO) Hamiltonian with the aid of the DMRG algorithm. Owing to the spin-adaptation of the algorithm, the total spin S is a good quantum number for these states. After the non-relativistic DMRG calculation is finished, all magnetic sublevels of the calculated states are constructed explicitly, and the SOC operator is expanded in the resulting basis. To this end, spin orbit coupled energies and wavefunctions are obtained as eigenvalues and eigenfunctions of the full Hamiltonian matrix which is composed of the SOC operator matrix and the BO Hamiltonian matrix. This treatment corresponds to a quasi-degenerate perturbation theory approach and can be regarded as the molecular equivalent to atomic Russell-Saunders coupling. For the evaluation of SOC matrix elements, the full Breit-Pauli SOC Hamiltonian is approximated by the widely used spin-orbit mean field operator. This operator allows for an efficient use of the second quantized triplet replacement operators that are readily generated during the non-relativistic DMRG algorithm, together with the Wigner-Eckart theorem. With a set of spin-orbit coupled wavefunctions at hand, the molecular g-tensors are calculated following the scheme proposed by Gerloch and McMeeking. It interprets the effective molecular g-values as the slope of the energy difference between the lowest Kramers pair with respect to the strength of the applied magnetic field. Test calculations on a chemically relevant Mo complex demonstrate the capabilities of the presented method.

The Journal of Chemical Physics•Jul 28, 2015•29 citations

A Cost-Effective Treatment of Spin–Orbit Couplings in the State-Averaged Driven Similarity Renormalization Group Second-Order Perturbation Theory

Meng Wang, Chenyang Li

We present an economical approach to treat spin-orbit coupling (SOC) in the state-averaged driven similarity renormalization group second-order perturbation theory (SA-DSRG-PT2). The electron correlation is first introduced by forming the SA-DSRG-PT2 dressed spin-free Hamiltonian. This Hamiltonian is then augmented with the Breit-Pauli Hamiltonian and diagonalized using spin-pure reference states to obtain the SOC-corrected energy spectrum. The spin-orbit mean-field approximation is also assumed to reduce the cost associated with the two-electron spin-orbit integrals. The resulting method is termed BP1-SA-DSRG-PT2c, and it possesses the same computational scaling as the non-relativistic counterpart, where only the one- and two-body density cumulants are required to obtain the vertical transition energy. The accuracy of BP1-SA-DSRG-PT2c is assessed on a few atoms and small molecules, including main-group diatomic molecules, transition-metal atoms, and actinide dioxide cations. Numerical results suggest that BP1-SA-DSRG-PT2c performs comparably to other internally contracted multireference perturbation theories with SOC treated using the state interaction scheme.

Molecules•May 7, 2025•4 citations

Efficient evaluation of spin–orbit couplings in single-molecule magnets using DMRG within the DMET framework

Xinyu Sun, Haibo Ma

The accurate characterization of complex electronic structures in single-molecule magnets (SMMs) is crucial for guiding their rational design. However, conventional multireference wavefunction theory faces significant computational bottlenecks for these systems, stemming from their intricate ligand environments and the necessity for large active spaces to capture strong electron correlation. To overcome these limitations, we introduce an efficient quantum embedding framework that synergistically combines density matrix embedding theory with ab initio density matrix renormalization group Relativistic spin-orbit coupling (SOC) effect is incorporated through a state-interaction procedure using the eXact two-component atomic mean-field spin-orbit (X2CAMF-SO) operator. Benchmarking on representative single-ion magnets confirms that our DMET-embedded DMRG-SI approach delivers results highly consistent with all-electron calculations, while offering a clear path for systematic improvement via active space expansion. Furthermore, by investigating SOC effects in a chromium dimer, we demonstrate the method's potential for high-accuracy and scalable ab initio studies of larger-sized polynuclear SMMs.

The Journal of Chemical Physics•Oct 23, 2025•3 citations

Quantum chemistry using the density matrix renormalization group II

Alexander O. Mitrushenkov, Roberto Linguerri, Paolo Palmieri +1 more

We have compared different strategies for ab initio quantum chemistry density matrix renormalization group treatments. The two starting orbital blocks include all valence and active orbitals of the reference complete active space self consistent field wave function. To generate the remaining blocks, we propose following the order of the contributions to the correlation energy: a posteriori using approximate occupation numbers or a priori, using a Møller–Plesset type of arguments, by explicit evaluation of second-order interactions. We have compared two different schemes for orbital localization to identify the important and less important orbital interactions and simplify the generation of the orbital blocks. To truncate the expansion we have compared two approaches, keeping constant the number m of components or the threshold λ to fix the residue of the expansion at each step. The extrapolation of the energies is found to provide accurate estimates of the full configuration interaction energy, making the expansion independent on the actual values of the two parameters m and λ. We propose to generate the factors for the two blocks from ground and excited eigenvectors of the Hamiltonian matrix.

The Journal of Chemical Physics•Aug 14, 2003•72 citations

A state interaction spin-orbit coupling density matrix renormalization group method

Elvira R. Sayfutyarova, Garnet Kin‐Lic Chan

We describe a state interaction spin-orbit (SISO) coupling method using density matrix renormalization group (DMRG) wavefunctions and the spin-orbit mean-field (SOMF) operator. We implement our DMRG-SISO scheme using a spin-adapted algorithm that computes transition density matrices between arbitrary matrix product states. To demonstrate the potential of the DMRG-SISO scheme we present accurate benchmark calculations for the zero-field splitting of the copper and gold atoms, comparing to earlier complete active space self-consistent-field and second-order complete active space perturbation theory results in the same basis. We also compute the effects of spin-orbit coupling on the spin-ladder of the iron-sulfur dimer complex [Fe2S2(SCH3)4](3-), determining the splitting of the lowest quartet and sextet states. We find that the magnitude of the zero-field splitting for the higher quartet and sextet states approaches a significant fraction of the Heisenberg exchange parameter.

The Journal of Chemical Physics•Jun 15, 2016•49 citations

Ab initio density matrix renormalization group study of magnetic coupling in dinuclear iron and chromium complexes

Travis V. Harris, Yuki Kurashige, Takeshi Yanai +1 more

The applicability of ab initio multireference wavefunction-based methods to the study of magnetic complexes has been restricted by the quickly rising active-space requirements of oligonuclear systems and dinuclear complexes with S > 1 spin centers. Ab initio density matrix renormalization group (DMRG) methods built upon an efficient parameterization of the correlation network enable the use of much larger active spaces, and therefore may offer a way forward. Here, we apply DMRG-CASSCF to the dinuclear complexes [Fe2OCl6](2-) and [Cr2O(NH3)10](4+). After developing the methodology through systematic basis set and DMRG M testing, we explore the effects of extended active spaces that are beyond the limit of conventional methods. We find that DMRG-CASSCF with active spaces including the metal d orbitals, occupied bridging-ligand orbitals, and their virtual double shells already capture a major portion of the dynamic correlation effects, accurately reproducing the experimental magnetic coupling constant (J) of [Fe2OCl6](2-) with (16e,26o), and considerably improving the smaller active space results for [Cr2O(NH3)10](4+) with (12e,32o). For comparison, we perform conventional MRCI+Q calculations and find the J values to be consistent with those from DMRG-CASSCF. In contrast to previous studies, the higher spin states of the two systems show similar deviations from the Heisenberg spectrum, regardless of the computational method.

The Journal of Chemical Physics•Feb 3, 2014•45 citations
The efficient and reliable treatment of both spin-orbit coupling (SOC) and electron correlation is essential for understanding f-element chemistry. We analyze two approaches to the problem: the one-step approach, where both effects are treated simultaneously, and the two-step state interaction approach. We report an implementation of the ab initio density matrix renormalization group with a one-step treatment of the SOC effect, which can be compared to prior two-step treatments on an equal footing. Using a dysprosium octahedral complex and bridged dimer as benchmark systems, we identify characteristics of problems where the one-step approach is beneficial for obtaining the low-energy spectrum.
Spin Orbit Coupling (SOC) is introduced to molecular ab initio density matrix renormalization group (DMRG) calculations. In the presented scheme, one first approximates the electronic ground state and a number of excited states of the Born-Oppenheimer (BO) Hamiltonian with the aid of the DMRG algorithm. Owing to the spin-adaptation of the algorithm, the total spin S is a good quantum number for these states. After the non-relativistic DMRG calculation is finished, all magnetic sublevels of the calculated states are constructed explicitly, and the SOC operator is expanded in the resulting basis. To this end, spin orbit coupled energies and wavefunctions are obtained as eigenvalues and eigenfunctions of the full Hamiltonian matrix which is composed of the SOC operator matrix and the BO Hamiltonian matrix. This treatment corresponds to a quasi-degenerate perturbation theory approach and can be regarded as the molecular equivalent to atomic Russell-Saunders coupling. For the evaluation of SOC matrix elements, the full Breit-Pauli SOC Hamiltonian is approximated by the widely used spin-orbit mean field operator. This operator allows for an efficient use of the second quantized triplet replacement operators that are readily generated during the non-relativistic DMRG algorithm, together with the Wigner-Eckart theorem. With a set of spin-orbit coupled wavefunctions at hand, the molecular g-tensors are calculated following the scheme proposed by Gerloch and McMeeking. It interprets the effective molecular g-values as the slope of the energy difference between the lowest Kramers pair with respect to the strength of the applied magnetic field. Test calculations on a chemically relevant Mo complex demonstrate the capabilities of the presented method.
We present an economical approach to treat spin-orbit coupling (SOC) in the state-averaged driven similarity renormalization group second-order perturbation theory (SA-DSRG-PT2). The electron correlation is first introduced by forming the SA-DSRG-PT2 dressed spin-free Hamiltonian. This Hamiltonian is then augmented with the Breit-Pauli Hamiltonian and diagonalized using spin-pure reference states to obtain the SOC-corrected energy spectrum. The spin-orbit mean-field approximation is also assumed to reduce the cost associated with the two-electron spin-orbit integrals. The resulting method is termed BP1-SA-DSRG-PT2c, and it possesses the same computational scaling as the non-relativistic counterpart, where only the one- and two-body density cumulants are required to obtain the vertical transition energy. The accuracy of BP1-SA-DSRG-PT2c is assessed on a few atoms and small molecules, including main-group diatomic molecules, transition-metal atoms, and actinide dioxide cations. Numerical results suggest that BP1-SA-DSRG-PT2c performs comparably to other internally contracted multireference perturbation theories with SOC treated using the state interaction scheme.
The accurate characterization of complex electronic structures in single-molecule magnets (SMMs) is crucial for guiding their rational design. However, conventional multireference wavefunction theory faces significant computational bottlenecks for these systems, stemming from their intricate ligand environments and the necessity for large active spaces to capture strong electron correlation. To overcome these limitations, we introduce an efficient quantum embedding framework that synergistically combines density matrix embedding theory with ab initio density matrix renormalization group Relativistic spin-orbit coupling (SOC) effect is incorporated through a state-interaction procedure using the eXact two-component atomic mean-field spin-orbit (X2CAMF-SO) operator. Benchmarking on representative single-ion magnets confirms that our DMET-embedded DMRG-SI approach delivers results highly consistent with all-electron calculations, while offering a clear path for systematic improvement via active space expansion. Furthermore, by investigating SOC effects in a chromium dimer, we demonstrate the method's potential for high-accuracy and scalable ab initio studies of larger-sized polynuclear SMMs.
We have compared different strategies for ab initio quantum chemistry density matrix renormalization group treatments. The two starting orbital blocks include all valence and active orbitals of the reference complete active space self consistent field wave function. To generate the remaining blocks, we propose following the order of the contributions to the correlation energy: a posteriori using approximate occupation numbers or a priori, using a Møller–Plesset type of arguments, by explicit evaluation of second-order interactions. We have compared two different schemes for orbital localization to identify the important and less important orbital interactions and simplify the generation of the orbital blocks. To truncate the expansion we have compared two approaches, keeping constant the number m of components or the threshold λ to fix the residue of the expansion at each step. The extrapolation of the energies is found to provide accurate estimates of the full configuration interaction energy, making the expansion independent on the actual values of the two parameters m and λ. We propose to generate the factors for the two blocks from ground and excited eigenvectors of the Hamiltonian matrix.
We describe a state interaction spin-orbit (SISO) coupling method using density matrix renormalization group (DMRG) wavefunctions and the spin-orbit mean-field (SOMF) operator. We implement our DMRG-SISO scheme using a spin-adapted algorithm that computes transition density matrices between arbitrary matrix product states. To demonstrate the potential of the DMRG-SISO scheme we present accurate benchmark calculations for the zero-field splitting of the copper and gold atoms, comparing to earlier complete active space self-consistent-field and second-order complete active space perturbation theory results in the same basis. We also compute the effects of spin-orbit coupling on the spin-ladder of the iron-sulfur dimer complex [Fe2S2(SCH3)4](3-), determining the splitting of the lowest quartet and sextet states. We find that the magnitude of the zero-field splitting for the higher quartet and sextet states approaches a significant fraction of the Heisenberg exchange parameter.
The applicability of ab initio multireference wavefunction-based methods to the study of magnetic complexes has been restricted by the quickly rising active-space requirements of oligonuclear systems and dinuclear complexes with S > 1 spin centers. Ab initio density matrix renormalization group (DMRG) methods built upon an efficient parameterization of the correlation network enable the use of much larger active spaces, and therefore may offer a way forward. Here, we apply DMRG-CASSCF to the dinuclear complexes [Fe2OCl6](2-) and [Cr2O(NH3)10](4+). After developing the methodology through systematic basis set and DMRG M testing, we explore the effects of extended active spaces that are beyond the limit of conventional methods. We find that DMRG-CASSCF with active spaces including the metal d orbitals, occupied bridging-ligand orbitals, and their virtual double shells already capture a major portion of the dynamic correlation effects, accurately reproducing the experimental magnetic coupling constant (J) of [Fe2OCl6](2-) with (16e,26o), and considerably improving the smaller active space results for [Cr2O(NH3)10](4+) with (12e,32o). For comparison, we perform conventional MRCI+Q calculations and find the J values to be consistent with those from DMRG-CASSCF. In contrast to previous studies, the higher spin states of the two systems show similar deviations from the Heisenberg spectrum, regardless of the computational method.