The Journal of Chemical Physics

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

Michael Roemelt

Published Jul 28, 2015• 29 citations

Abstract

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.

DOI: 10.1063/1.4927432

Related papers

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 PhysicsJun 15, 201649 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 PhysicsOct 23, 20253 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.

MoleculesMay 7, 20254 citations

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

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.

arXiv:2207.02435The Journal of Chemical PhysicsSep 29, 202219 citations

Spin-adapted density matrix renormalization group algorithms for quantum chemistry

Sandeep Sharma, Garnet Kin-Lic Chan

We extend the spin-adapted density matrix renormalization group (DMRG) algorithm of McCulloch and Gulacsi [Europhys. Lett. 57, 852 (2002)] to quantum chemical Hamiltonians. This involves using a quasi-density matrix, to ensure that the renormalized DMRG states are eigenfunctions of Ŝ(2), and the Wigner-Eckart theorem, to reduce overall storage and computational costs. We argue that the spin-adapted DMRG algorithm is most advantageous for low spin states. Consequently, we also implement a singlet-embedding strategy due to Tatsuaki [Phys. Rev. E 61, 3199 (2000)] where we target high spin states as a component of a larger fictitious singlet system. Finally, we present an efficient algorithm to calculate one- and two-body reduced density matrices from the spin-adapted wavefunctions. We evaluate our developments with benchmark calculations on transition metal system active space models. These include the Fe(2)S(2), [Fe(2)S(2)(SCH(3))(4)](2-), and Cr(2) systems. In the case of Fe(2)S(2), the spin-ladder spacing is on the microHartree scale, and here we show that we can target such very closely spaced states. In [Fe(2)S(2)(SCH(3))(4)](2-), we calculate particle and spin correlation functions, to examine the role of sulfur bridging orbitals in the electronic structure. In Cr(2) we demonstrate that spin-adaptation with the Wigner-Eckart theorem and using singlet embedding can yield up to an order of magnitude increase in computational efficiency. Overall, these calculations demonstrate the potential of using spin-adaptation to extend the range of DMRG calculations in complex transition metal problems.

arXiv:1408.5039The Journal of Chemical PhysicsMar 28, 2012295 citations

Direct density matrix renormalization group approaches for strong correlation effects in quantum chemistry

Philipp Bertone Snajberk

An implementation of the single-site density matrix renormalization group (DMRG) approach for the quantum-chemical application to calculate approximate as well as exact electronic ground-state energies and wave functions of atoms and molecules is presented.The algorithm is formulated in terms of the variational and sizeconsistent class of wave functions, matrix product states, and their operator analogues, matrix product operators.A procedure to construct quantum-chemical Hamiltonian operators with long-range Coulomb interactions as matrix product operators is developed.The many-electron Hamiltonian operators are represented as sums and products of one-electron operators.For the efficient matrix product operator construction of quantum-chemical Hamiltonians, different basis spaces including discrete position and momentum space as well as molecular and atomic orbitals are investigated.Discrete position space yields the best matrix product operator size scaling, but introduces artificial long-range couplings in the kinetic energy term of the Hamiltonian which are a consequence of the DMRG algorithm in terms of matrix product states.The huge amount of discrete space points needed for an accurate description poses a problem because of the electron-electron interaction.Therefore, a momentum space or Fourier series expansion of the electron-electron interaction is shown which makes it possible to compress its matrix product operator representation.DMRG applied to quantum-chemical Hamiltonians in matrix product operator form can serve as means to directly (without mean-field treatment in advance) investigate quantum-chemical problems and find the analytically exact ground-state solution of small atoms and molecules.The next basis space analyzed in detail is the atomic orbital basis.The matrix product operator construction and the direct application of the DMRG method in this basis is examined.A physical compression scheme based on a singular value decomposition of the electron-electron interaction in the atomic orbital basis is demonstrated.Furthermore, DMRG calculations in the atomic orbital basis for certain systems where mean-field theory yields qualitatively wrong results and where higher-order approximative correction methods like truncated configuration interaction and coupled cluster cannot fully capture strong electron correlation effects are compared to results obtained in the molecular orbital basis.The comparison shows that atomic orbitals describe the physical situation well and thus lead to a tremendous reduction of the matrix product state dimension and the computational effort.i

Jan 1, 20171 citation