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

The Journal of Chemical Physics

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

Xinyu Sun, Haibo Ma

Published Oct 23, 2025• 3 citations
Publisher

Abstract

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.

DOI
10.1063/5.0297017

Related papers

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

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

High-performance ab initio density matrix renormalization group method: Applicability to large-scale multireference problems for metal compounds

Yuki Kurashige, Takeshi Yanai

This article presents an efficient and parallelized implementation of the density matrix renormalization group (DMRG) algorithm for quantum chemistry calculations. The DMRG method as a large-scale multireference electronic structure model is by nature particularly efficient for one-dimensionally correlated systems, while the present development is oriented toward applications for polynuclear transition metal compounds, in which the macroscopic one-dimensional structure of electron correlation is absent. A straightforward extension of the DMRG algorithm is proposed with further improvements and aggressive optimizations to allow its application with large multireference active space, which is often demanded for metal compound calculations. Special efficiency is achieved by making better use of sparsity and symmetry in the operator and wave function representations. By accomplishing computationally intensive DMRG calculations, the authors have found that a large number of renormalized basis states are required to represent high entanglement of the electron correlation for metal compound applications, and it is crucial to adopt auxiliary perturbative correction to the projected density matrix during the DMRG sweep optimization in order to attain proper convergence to the solution. Potential energy curve calculations for the Cr(2) molecule near the known equilibrium precisely predicted the full configuration interaction energies with a correlation space of 24 electrons in 30 orbitals [denoted by (24e,30o)]. The energies are demonstrated to be accurate to 0.6mE(h) (the error from the extrapolated best value) when as many as 10,000 renormalized basis states are employed for the left and right DMRG block representations. The relative energy curves for [Cu(2)O(2)](2+) along the isomerization coordinate were obtained from DMRG and other correlated calculations, for which a fairly large orbital space (32e,62o) is modeled as a full correlation space. The DMRG prediction nearly overlaps with the energy curve from the coupled cluster with singles, doubles, and perturbative triple [CCSD(T)] calculations, while the multireference complete active space self-consistent field (CASSCF) calculations with the small reference configuration (8e,8o) are found to overestimate the biradical character of the electronic state of [Cu(2)O(2)](2+) according to the one-electron density matrix analysis.

The Journal of Chemical Physics•Jun 19, 2009•193 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 Theory for Multireference Quantum Chemistry

Takeshi Yanai

A major challenge in molecular electronic structure theory arises from chemical systems where the mean-field picture based upon a single electronic configuration fails to properly describe mechanics of many electrons. Such complex electronic states should be accounted for using a correlated superposition of multiple determinants and are referred to as multireference in quantum chemistry. Application domain of traditional multireference methods is severely limited to systems containing a small number of correlated electrons because of their exponential complexity. In recent years, much attention has been drawn to the density matrix renormalization group (DMRG) theory as a promising alternative. It has been introduced in ab initio quantum chemistry calculations and shown to be a highly-scalable multireference method that can handle full quantum degrees of freedom for the strongly-correlated wavefunction with unprecedented large configuration space. In this article, we provide an explicative introduction to the underlying theory and formalism built into the DMRG method and give a brief overview of our methodological extensions that combine DMRG and active-space correlation models.

Molecular Science•Jan 1, 2014•0 citations

Calculation of NMR Shieldings for Strongly Correlated Molecules Using the Density Matrix Renormalization Group Self-Consistent-Field Method with Large Active Spaces

Jiaao Su, Wei Huang, Yinxuan Song +1 more

Accurate prediction of nuclear magnetic resonance (NMR) shielding constants for strongly correlated systems remains challenging. In this context, we present a density matrix renormalization group self-consistent-field (DMRG-SCF) implementation for NMR shielding calculations with large active spaces. Within a gauge-including atomic orbital DMRG-SCF (GIAO-DMRG-SCF) framework, we formulate and solve the coupled-perturbed DMRG-SCF (CP-DMRG-SCF) equations using local-site effective Hamiltonians. This approach enables the direct evaluation of the wave function's response to an external magnetic field. In the limit of small active spaces, the computed NMR shieldings are in agreement with conventional complete active space self-consistent-field (CASSCF) results, while the present approach naturally extends to larger active spaces inaccessible to standard multiconfigurational methods. Applications to a series of molecular systems, particularly strongly correlated transition metal complexes, demonstrate that the proposed implementation can be both reliable and accurate.

Journal of Chemical Theory and Computation•Sep 2, 2026•0 citations
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 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.
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.
This article presents an efficient and parallelized implementation of the density matrix renormalization group (DMRG) algorithm for quantum chemistry calculations. The DMRG method as a large-scale multireference electronic structure model is by nature particularly efficient for one-dimensionally correlated systems, while the present development is oriented toward applications for polynuclear transition metal compounds, in which the macroscopic one-dimensional structure of electron correlation is absent. A straightforward extension of the DMRG algorithm is proposed with further improvements and aggressive optimizations to allow its application with large multireference active space, which is often demanded for metal compound calculations. Special efficiency is achieved by making better use of sparsity and symmetry in the operator and wave function representations. By accomplishing computationally intensive DMRG calculations, the authors have found that a large number of renormalized basis states are required to represent high entanglement of the electron correlation for metal compound applications, and it is crucial to adopt auxiliary perturbative correction to the projected density matrix during the DMRG sweep optimization in order to attain proper convergence to the solution. Potential energy curve calculations for the Cr(2) molecule near the known equilibrium precisely predicted the full configuration interaction energies with a correlation space of 24 electrons in 30 orbitals [denoted by (24e,30o)]. The energies are demonstrated to be accurate to 0.6mE(h) (the error from the extrapolated best value) when as many as 10,000 renormalized basis states are employed for the left and right DMRG block representations. The relative energy curves for [Cu(2)O(2)](2+) along the isomerization coordinate were obtained from DMRG and other correlated calculations, for which a fairly large orbital space (32e,62o) is modeled as a full correlation space. The DMRG prediction nearly overlaps with the energy curve from the coupled cluster with singles, doubles, and perturbative triple [CCSD(T)] calculations, while the multireference complete active space self-consistent field (CASSCF) calculations with the small reference configuration (8e,8o) are found to overestimate the biradical character of the electronic state of [Cu(2)O(2)](2+) according to the one-electron density matrix analysis.
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.
A major challenge in molecular electronic structure theory arises from chemical systems where the mean-field picture based upon a single electronic configuration fails to properly describe mechanics of many electrons. Such complex electronic states should be accounted for using a correlated superposition of multiple determinants and are referred to as multireference in quantum chemistry. Application domain of traditional multireference methods is severely limited to systems containing a small number of correlated electrons because of their exponential complexity. In recent years, much attention has been drawn to the density matrix renormalization group (DMRG) theory as a promising alternative. It has been introduced in ab initio quantum chemistry calculations and shown to be a highly-scalable multireference method that can handle full quantum degrees of freedom for the strongly-correlated wavefunction with unprecedented large configuration space. In this article, we provide an explicative introduction to the underlying theory and formalism built into the DMRG method and give a brief overview of our methodological extensions that combine DMRG and active-space correlation models.
Accurate prediction of nuclear magnetic resonance (NMR) shielding constants for strongly correlated systems remains challenging. In this context, we present a density matrix renormalization group self-consistent-field (DMRG-SCF) implementation for NMR shielding calculations with large active spaces. Within a gauge-including atomic orbital DMRG-SCF (GIAO-DMRG-SCF) framework, we formulate and solve the coupled-perturbed DMRG-SCF (CP-DMRG-SCF) equations using local-site effective Hamiltonians. This approach enables the direct evaluation of the wave function's response to an external magnetic field. In the limit of small active spaces, the computed NMR shieldings are in agreement with conventional complete active space self-consistent-field (CASSCF) results, while the present approach naturally extends to larger active spaces inaccessible to standard multiconfigurational methods. Applications to a series of molecular systems, particularly strongly correlated transition metal complexes, demonstrate that the proposed implementation can be both reliable and accurate.