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:2112.09594 · Physical review. B./Physical review. B

Combined unitary and symmetric group approach applied to low-dimensional Heisenberg spin systems

Werner Dobrautz, Vamshi M. Katukuri, Nikolay A. Bogdanov, Daniel Kats, Giovanni Li Manni, Ali Alavi

Published May 18, 2022• 27 citations
arXivHTMLPDF

Abstract

A novel combined unitary and symmetric group approach is used to study the spin-$\frac{1}{2}$ Heisenberg model and related Fermionic systems in a total spin-adapted representation, using a linearly-parameterised Ansatz for the many-body wave function. We show that a more compact ground-state wave function representation---indicated by a larger leading ground-state coefficient---is obtained when combining the symmetric group ${\mathcal{S}}_{n}$, in the form of permutations of the underlying lattice site ordering, with the cumulative spin coupling based on the unitary group $\mathrm{U}(n)$. In one-dimensional systems the observed compression of the wave function is reminiscent of block-spin renormalization group approaches, and allows us to study larger lattices (here taken up to 80 sites) with the spin-adapted full configuration interaction quantum Monte Carlo method, which benefits from the sparsity of the Hamiltonian matrix and the corresponding sampled eigenstates that emerge from the reordering. We find that in an optimal lattice ordering the configuration state function with highest weight already captures with high accuracy the spin-spin correlation function of the exact ground-state wave function. This feature is found for more general lattice models, such as the Hubbard model, and ab initio quantum chemical models, exemplified by one-dimensional hydrogen chains. We also provide numerical evidence that the optimal lattice ordering for the unitary group approach is not generally equivalent to the optimal ordering obtained for methods based on matrix-product states, such as the density-matrix renormalization group approach.

DOI
10.1103/physrevb.105.195123

Related papers

Resolution of Ground and Excited Electronic States within the Unitary Group Approach

Giovanni Li Manni, Daniel Kats, Niklas Liebermann

In this work ground and excited electronic states of Heisenberg cluster models, in the form of con- figuration interaction many-body wave functions, are characterized within the spin-adapted Graph- ical Unitary Group Approach framework, and relying on a novel combined unitary and symmetric group approach. Finite-size cluster models of well-defined point-group symmetry and of general local-spin Slocal ≥ 21 are presented, including J1–J2 triangular and tetrahedral clusters, which are often used to describe magnetic interactions in biological and bio-mimetic polynuclear transition metal clusters with unique catalytic activity, such as nitrogen fixation and photosynthesis. We show that a unique block-diagonal structure of the underlying Hamiltonian matrix in the spin-adapted basis emerges when an optimal lattice site ordering is chosen that reflects the internal symmetries of the model investigated. The block-diagonal structure is bound to the commutation relations between cumulative spin operators and the Hamiltonian operator, that in turn depend on the geometry of the cluster investigated. The many-body basis transformation, in the form of the orbital/site reorder- ing, exposes such commutation relations. These commutation relations represent a rigorous and formal demonstration of the block-diagonal structure in Hamiltonian matrices and the compression of the corresponding spin-adapted many-body wave functions. As a direct consequence of the block- diagonal structure of the Hamiltonian matrix it is possible to selectively optimize electronic excited states without the overhead of calculating the lower-energy states by simply relying on the initial ansatz for the targeted wave function. Additionally, more compact many-body wave functions are obtained. In extreme cases, electronic states are precisely described by a single configuration state function, despite the curse of dimensionality of the corresponding Hilbert space. These findings are crucial in the electronic structure theory framework, for they offer a conceptual route towards wave functions of reduced multi-reference character, that can be optimized more easily by approximated eigensolvers and are of more facile physical interpretation. They open the way to study larger ab initio and model Hamiltonians of increasingly larger number of correlated electrons, while keeping the computational costs at their lowest.

•Oct 12, 2022•0 citations

Density-matrix algorithms for quantum renormalization groups

Steven R. White

A formulation of numerical real-space renormalization groups for quantum many-body problems is presented and several algorithms utilizing this formulation are outlined. The methods are presented and demonstrated using S=1/2 and S=1 Heisenberg chains as test cases. The key idea of the formulation is that rather than keep the lowest-lying eigenstates of the Hamiltonian in forming a new effective Hamiltonian of a block of sites, one should keep the most significant eigenstates of the block density matrix, obtained from diagonalizing the Hamiltonian of a larger section of the lattice which includes the block. This approach is much more accurate than the standard approach; for example, energies for the S=1 Heisenberg chain can be obtained to an accuracy of at least ${10}^{\mathrm{\ensuremath{-}}9}$. The method can be applied to almost any one-dimensional quantum lattice system, and can provide a wide variety of static properties.

Physical review. B, Condensed matter•Oct 1, 1993•3,244 citations

Applications of Density Matrix Renormalization Group to Problems in Magnetism

G A Gehring, Robert J. Bursill, Tao Xiang

The application of real space renormalization group methods to quantum lattice models has become a topic of great interest following the development of the density matrix renormalization group by White.This method has been used to find the ground and low-lying excited state energies and wave functions of quantum spin models in which the form of the ground state is not clear, for instance because the interactions are frustrated.It has also been applied to fermion problems where the tendency for localization due to the strong Coulomb repulsion is opposed by the lowering of the kinetic energy which occurs as a result of electron transfer.The approach is particularly suitable for one-, or quasi-one-dimensional problems.The method involves truncating the Hilbert space in a systematic and optimised manner.Results for the ground state energy are thus variational bounds.The results for low-lying energies and correlation functions for one-dimensional systems have unprecedented accuracy and the method has become the method of choice for solving one-dimensional quantum spin problems.We review the method and results obtained for the spin-1 chain with biquadratic exchange as well as the spin-1/2 model with competing nearest and next nearest neighbour exchange will be described.More recently, the density matrix renormalization group has been applied to reformulate the coupling constant renormalization group approach which is appropriate for the study of critical properties.This approach has been applied to the anisotropic spin-1/2 Heisenberg chain.Finally, we discuss recent work which has borne promising applications in two dimensions -the Ising model and the two-dimensional Hubbard model.

Acta Physica Polonica A•Jan 1, 1997•18 citations

Renormalization-group method for weakly coupled quantum chains: Application to the spin- 12 Heisenberg model

S. Moukouri

The Kato-Bloch perturbation formalism is used to present a density-matrix renormalization-group method for strongly anisotropic two-dimensional systems. This method is used to study Heisenberg chains weakly coupled by the transverse couplings ${J}_{\ensuremath{\perp}}$ and ${J}_{d}$ (along the diagonals). An extensive comparison of the renormalization group and quantum Monte Carlo results for parameters where the simulations by the latter method are possible shows a very good agreement between the two methods. It is found, by analyzing ground state energies and spin-spin correlation functions, that there is a transition between two ordered magnetic states. When ${J}_{d}∕{J}_{\ensuremath{\perp}}\ensuremath{\lesssim}0.5$, the ground state displays a N\'eel order. When ${J}_{d}∕{J}_{\ensuremath{\perp}}\ensuremath{\gtrsim}0.5$, a collinear magnetic ground state in which interchain spin correlations are ferromagnetic becomes stable. In the vicinity of the transition point, ${J}_{d}∕{J}_{\ensuremath{\perp}}\ensuremath{\approx}0.5$, the ground state is disordered. But, the nature of this disordered ground state is unclear. While the numerical data seem to show that the chains are disconnected, the possibility of a genuine disordered two-dimensional state, hidden by finite size effects, cannot be excluded.

arXiv:cond-mat/0305608•Physical Review B•Jul 7, 2004•25 citations

Incorporation of density-matrix wave functions in Monte Carlo simulations: Application to the frustrated Heisenberg model

M. S. L. du Croo de Jongh, J. M. J. van Leeuwen, Wim van Saarloos

We combine the density-matrix renormalization-group (DMRG) technique with Green function Monte Carlo (GFMC) simulations using a special representation of the DMRG wave function. As a test case we apply the method to the two-dimensional frustrated Heisenberg antiferromagnet. By supplementing the branching in GFMC simulations with stochastic reconfiguration we get a stable simulation with a small variance also in the region where the fluctuations due to the minus sign problem are maximal. The sensitivity of the results to the choice of the guiding wave function is extensively investigated. In agreement with earlier calculations it is found from the DMRG wave function that for small ratios of the next-nearest-to-nearest neighbor coupling strength the system orders as a N\'eel-type antiferromagnet and for large ratios as a columnar antiferromagnet. The spin stiffness suggests an intermediate regime without magnetic long-range order. The energy curve indicates that the columnar phase is separated from the intermediate phase by a first-order transition. The combination of the DMRG and GFMC techniques allows us to substantiate this picture by calculating also the spin correlations in the system. We observe a pattern of spin correlations in the intermediate regime which is in between dimerlike and plaquette-type ordering, states that have recently been suggested. It is a state with strong dimerization in one direction and weaker dimerization in the perpendicular direction and thus it lacks the square symmetry of the plaquette state.

arXiv:cond-mat/0002116•Physical review. B, Condensed matter•Dec 1, 2000•45 citations

Strongly Correlated Models in One Dimension from a Density Matrix Renormalization Group Perspective

Stefan Rommer

The density matrix renormalization group (DMRG) is a numerical method for studying low dimensional strongly correlated models. The study of these models is one of the most active areas of condensed matter physics today and has been spurred not only by the discovery of high temperature superconductivity but also by experiments on new compounds. In this thesis we study fundamentals of the density matrix renormalization group method as well as applications of the method to Heisenberg spin chains. The introductory text is meant as a brief overview of the DMRG as well as Heisenberg spin systems and provides an introduction to the papers. In papers I and II we investigate the density matrix renormalization group with the aim to better understand the underlying structure of the method. We find that if the renormalization converges to a fixed point in the thermodynamic limit, the resulting basis states can be written as a matrix product form. We also show an explicit relation of the DMRG to variational calculations using the fact that the matrix product states can be rederived through a variational ansatz without any reference to the DMRG calculation. It is then shown how the matrix product states can be used to define a set of Bloch waves for describing elementary excitations of the Heisenberg spin-1 chain. The low lying energy spectrum is calculated. In papers III and IV we investigate the thermodynamics of impurities in a spin-1/2 chain. For this purpose we have shown that a recently developed transfer matrix DMRG method can be used to calculate impurity properties. Both impurity corrections to thermodynamic properties and local properties at finite temperature can be obtained. The method gives quicker and more accurate results than quantum Monte Carlo. The numerical results show good agreement with field theory calculations for a number of impurity configurations and the cross-over function for the magnetic susceptibility from high- to low-temperature behavior shows the expected data collapse. Our results also give quantitative predictions for experiments on doped spin-1/2 chains.

•Jan 1, 1998•0 citations
A novel combined unitary and symmetric group approach is used to study the spin-12\frac{1}{2}21​ Heisenberg model and related Fermionic systems in a total spin-adapted representation, using a linearly-parameterised Ansatz for the many-body wave function. We show that a more compact ground-state wave function representation---indicated by a larger leading ground-state coefficient---is obtained when combining the symmetric group Sn{\mathcal{S}}_{n}Sn​, in the form of permutations of the underlying lattice site ordering, with the cumulative spin coupling based on the unitary group U(n)\mathrm{U}(n)U(n). In one-dimensional systems the observed compression of the wave function is reminiscent of block-spin renormalization group approaches, and allows us to study larger lattices (here taken up to 80 sites) with the spin-adapted full configuration interaction quantum Monte Carlo method, which benefits from the sparsity of the Hamiltonian matrix and the corresponding sampled eigenstates that emerge from the reordering. We find that in an optimal lattice ordering the configuration state function with highest weight already captures with high accuracy the spin-spin correlation function of the exact ground-state wave function. This feature is found for more general lattice models, such as the Hubbard model, and ab initio quantum chemical models, exemplified by one-dimensional hydrogen chains. We also provide numerical evidence that the optimal lattice ordering for the unitary group approach is not generally equivalent to the optimal ordering obtained for methods based on matrix-product states, such as the density-matrix renormalization group approach.
In this work ground and excited electronic states of Heisenberg cluster models, in the form of con- figuration interaction many-body wave functions, are characterized within the spin-adapted Graph- ical Unitary Group Approach framework, and relying on a novel combined unitary and symmetric group approach. Finite-size cluster models of well-defined point-group symmetry and of general local-spin Slocal ≥ 21 are presented, including J1–J2 triangular and tetrahedral clusters, which are often used to describe magnetic interactions in biological and bio-mimetic polynuclear transition metal clusters with unique catalytic activity, such as nitrogen fixation and photosynthesis. We show that a unique block-diagonal structure of the underlying Hamiltonian matrix in the spin-adapted basis emerges when an optimal lattice site ordering is chosen that reflects the internal symmetries of the model investigated. The block-diagonal structure is bound to the commutation relations between cumulative spin operators and the Hamiltonian operator, that in turn depend on the geometry of the cluster investigated. The many-body basis transformation, in the form of the orbital/site reorder- ing, exposes such commutation relations. These commutation relations represent a rigorous and formal demonstration of the block-diagonal structure in Hamiltonian matrices and the compression of the corresponding spin-adapted many-body wave functions. As a direct consequence of the block- diagonal structure of the Hamiltonian matrix it is possible to selectively optimize electronic excited states without the overhead of calculating the lower-energy states by simply relying on the initial ansatz for the targeted wave function. Additionally, more compact many-body wave functions are obtained. In extreme cases, electronic states are precisely described by a single configuration state function, despite the curse of dimensionality of the corresponding Hilbert space. These findings are crucial in the electronic structure theory framework, for they offer a conceptual route towards wave functions of reduced multi-reference character, that can be optimized more easily by approximated eigensolvers and are of more facile physical interpretation. They open the way to study larger ab initio and model Hamiltonians of increasingly larger number of correlated electrons, while keeping the computational costs at their lowest.
A formulation of numerical real-space renormalization groups for quantum many-body problems is presented and several algorithms utilizing this formulation are outlined. The methods are presented and demonstrated using S=1/2 and S=1 Heisenberg chains as test cases. The key idea of the formulation is that rather than keep the lowest-lying eigenstates of the Hamiltonian in forming a new effective Hamiltonian of a block of sites, one should keep the most significant eigenstates of the block density matrix, obtained from diagonalizing the Hamiltonian of a larger section of the lattice which includes the block. This approach is much more accurate than the standard approach; for example, energies for the S=1 Heisenberg chain can be obtained to an accuracy of at least ${10}^{\mathrm{\ensuremath{-}}9}$. The method can be applied to almost any one-dimensional quantum lattice system, and can provide a wide variety of static properties.
The application of real space renormalization group methods to quantum lattice models has become a topic of great interest following the development of the density matrix renormalization group by White.This method has been used to find the ground and low-lying excited state energies and wave functions of quantum spin models in which the form of the ground state is not clear, for instance because the interactions are frustrated.It has also been applied to fermion problems where the tendency for localization due to the strong Coulomb repulsion is opposed by the lowering of the kinetic energy which occurs as a result of electron transfer.The approach is particularly suitable for one-, or quasi-one-dimensional problems.The method involves truncating the Hilbert space in a systematic and optimised manner.Results for the ground state energy are thus variational bounds.The results for low-lying energies and correlation functions for one-dimensional systems have unprecedented accuracy and the method has become the method of choice for solving one-dimensional quantum spin problems.We review the method and results obtained for the spin-1 chain with biquadratic exchange as well as the spin-1/2 model with competing nearest and next nearest neighbour exchange will be described.More recently, the density matrix renormalization group has been applied to reformulate the coupling constant renormalization group approach which is appropriate for the study of critical properties.This approach has been applied to the anisotropic spin-1/2 Heisenberg chain.Finally, we discuss recent work which has borne promising applications in two dimensions -the Ising model and the two-dimensional Hubbard model.
The Kato-Bloch perturbation formalism is used to present a density-matrix renormalization-group method for strongly anisotropic two-dimensional systems. This method is used to study Heisenberg chains weakly coupled by the transverse couplings ${J}_{\ensuremath{\perp}}$ and Jd{J}_{d}Jd​ (along the diagonals). An extensive comparison of the renormalization group and quantum Monte Carlo results for parameters where the simulations by the latter method are possible shows a very good agreement between the two methods. It is found, by analyzing ground state energies and spin-spin correlation functions, that there is a transition between two ordered magnetic states. When ${J}_{d}∕{J}_{\ensuremath{\perp}}\ensuremath{\lesssim}0.5$, the ground state displays a N\'eel order. When ${J}_{d}∕{J}_{\ensuremath{\perp}}\ensuremath{\gtrsim}0.5$, a collinear magnetic ground state in which interchain spin correlations are ferromagnetic becomes stable. In the vicinity of the transition point, ${J}_{d}∕{J}_{\ensuremath{\perp}}\ensuremath{\approx}0.5$, the ground state is disordered. But, the nature of this disordered ground state is unclear. While the numerical data seem to show that the chains are disconnected, the possibility of a genuine disordered two-dimensional state, hidden by finite size effects, cannot be excluded.
We combine the density-matrix renormalization-group (DMRG) technique with Green function Monte Carlo (GFMC) simulations using a special representation of the DMRG wave function. As a test case we apply the method to the two-dimensional frustrated Heisenberg antiferromagnet. By supplementing the branching in GFMC simulations with stochastic reconfiguration we get a stable simulation with a small variance also in the region where the fluctuations due to the minus sign problem are maximal. The sensitivity of the results to the choice of the guiding wave function is extensively investigated. In agreement with earlier calculations it is found from the DMRG wave function that for small ratios of the next-nearest-to-nearest neighbor coupling strength the system orders as a N\'eel-type antiferromagnet and for large ratios as a columnar antiferromagnet. The spin stiffness suggests an intermediate regime without magnetic long-range order. The energy curve indicates that the columnar phase is separated from the intermediate phase by a first-order transition. The combination of the DMRG and GFMC techniques allows us to substantiate this picture by calculating also the spin correlations in the system. We observe a pattern of spin correlations in the intermediate regime which is in between dimerlike and plaquette-type ordering, states that have recently been suggested. It is a state with strong dimerization in one direction and weaker dimerization in the perpendicular direction and thus it lacks the square symmetry of the plaquette state.
The density matrix renormalization group (DMRG) is a numerical method for studying low dimensional strongly correlated models. The study of these models is one of the most active areas of condensed matter physics today and has been spurred not only by the discovery of high temperature superconductivity but also by experiments on new compounds. In this thesis we study fundamentals of the density matrix renormalization group method as well as applications of the method to Heisenberg spin chains. The introductory text is meant as a brief overview of the DMRG as well as Heisenberg spin systems and provides an introduction to the papers. In papers I and II we investigate the density matrix renormalization group with the aim to better understand the underlying structure of the method. We find that if the renormalization converges to a fixed point in the thermodynamic limit, the resulting basis states can be written as a matrix product form. We also show an explicit relation of the DMRG to variational calculations using the fact that the matrix product states can be rederived through a variational ansatz without any reference to the DMRG calculation. It is then shown how the matrix product states can be used to define a set of Bloch waves for describing elementary excitations of the Heisenberg spin-1 chain. The low lying energy spectrum is calculated. In papers III and IV we investigate the thermodynamics of impurities in a spin-1/2 chain. For this purpose we have shown that a recently developed transfer matrix DMRG method can be used to calculate impurity properties. Both impurity corrections to thermodynamic properties and local properties at finite temperature can be obtained. The method gives quicker and more accurate results than quantum Monte Carlo. The numerical results show good agreement with field theory calculations for a number of impurity configurations and the cross-over function for the magnetic susceptibility from high- to low-temperature behavior shows the expected data collapse. Our results also give quantitative predictions for experiments on doped spin-1/2 chains.