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

Renormalization group analysis of dipolar Heisenberg model on square lattice

Ahmet Keleş, Erhai Zhao

Published Jun 6, 2018• 19 citations
arXivHTMLPDF

Abstract

We present a detailed functional renormalization group analysis of spin-$\frac{1}{2}$ dipolar Heisenberg model on square lattice. This model is similar to the well-known ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ model and describes the pseudospin degrees of freedom of polar molecules confined in deep optical lattice with long-range anisotropic dipole-dipole interactions. Previous study of this model based on tensor network ansatz indicates a paramagnetic ground state for certain dipole tilting angles which can be tuned in experiments to control the exchange couplings. The tensor ansatz formulated on a small cluster unit cell is inadequate to describe the spiral order, and therefore the phase diagram at high azimuthal tilting angles remains undetermined. Here, we obtain the full phase diagram of the model from numerical pseudofermion functional renormalization group calculations. We show that an extended quantum paramagnetic phase is realized between the N\'eel and stripe/spiral phases. In this region, the spin susceptibility flows smoothly down to the lowest numerical renormalization group scales with no sign of divergence or breakdown of the flow, in sharp contrast to the flow towards the long-range-ordered phases. Our results provide further evidence that the dipolar Heisenberg model is a fertile ground for quantum spin liquids.

DOI
10.1103/physrevb.97.245105

Related papers

Absence of Long-Range Order in a Triangular Spin System with Dipolar Interactions

Ahmet Keleş, Erhai Zhao

The antiferromagnetic Heisenberg model on the triangular lattice is perhaps the best known example of frustrated magnets, but it orders at low temperatures. Recent density matrix renormalization group (DMRG) calculations find that the next nearest neighbor interaction J_{2} enhances the frustration, and it leads to a spin liquid for J_{2}/J_{1}∈(0.08,0.15). In addition, a DMRG study of a dipolar Heisenberg model with longer range interactions gives evidence for a spin liquid at a small dipole tilting angle θ∈[0,10°). In both cases, the putative spin liquid region appears to be small. Here, we show that for the triangular lattice dipolar Heisenberg model, a robust quantum paramagnetic phase exists in a surprisingly wide region, θ∈[0,54°), for dipoles tilted along the lattice diagonal direction. We obtain the phase diagram of the model by functional renormalization group (RG), which treats all magnetic instabilities on equal footing. The quantum paramagnetic phase is characterized by a smooth continuous flow of vertex functions and spin susceptibility down to the lowest RG scale, in contrast to the apparent breakdown of RG flow in phases with stripe or spiral order. Our finding points to a promising direction to search for quantum spin liquids in ultracold dipolar molecules.

arXiv:1801.00842•Physical Review Letters•May 3, 2018•42 citations

Frustrated Magnetism of Dipolar Molecules on a Square Optical Lattice: Prediction of a Quantum Paramagnetic Ground State

Haiyuan Zou, Erhai Zhao, W. Vincent Liu

Motivated by the experimental realization of quantum spin models of polar molecule KRb in optical lattices, we analyze the spin 1/2 dipolar Heisenberg model with competing anisotropic, long-range exchange interactions. We show that, by tilting the orientation of dipoles using an external electric field, the dipolar spin system on square lattice comes close to a maximally frustrated region similar, but not identical, to that of the J_{1}-J_{2} model. This provides a simple yet powerful route to potentially realize a quantum spin liquid without the need for a triangular or kagome lattice. The ground state phase diagrams obtained from Schwinger-boson and spin-wave theories consistently show a spin disordered region between the Néel, stripe, and spiral phase. The existence of a finite quantum paramagnetic region is further confirmed by an unbiased variational ansatz based on tensor network states and a tensor renormalization group.

arXiv:1702.08517•Physical Review Letters•Jul 31, 2017•18 citations

Evidence of valence bond solid and quantum paramagnetism in the $J_1$--$J_2$--$J_3$ Heisenberg model from pseudofermion functional renormalization group

GÜN GÜNAL, Erhai Zhao, Ahmet Keleş

Motivated by recent tensor network studies reporting a quantum spin liquid (QSL) within the quantum paramagnetic (QP) regime, we investigate numerically the zero-temperature phase diagram of the $J_1$--$J_2$--$J_3$ Heisenberg model on the square lattice using the pseudofermion functional renormalization group (pf-fRG). Thanks to improved frequency resolution and larger lattice truncations, we are able to explore a significantly larger parameter space than previous pf-fRG work to map out the phase boundary of an extended QP region. We further resolve possible bond order within the QP regime by measuring the dimer response to small deformations in $J_1$, i.e. bond anisotropy that breaks the point-group symmetry of the Hamiltonian. The functional renormalization group flow of the dimer susceptibility suggests that the QP window splits into two: a region where it grows continuously down to the lowest energy scale to signal valence bond solid (VBS) order, and another region where it reaches a peak at a finite scale and is subsequently suppressed to indicate a QSL.

arXiv:2610.26135•Oct 19, 2026•0 citations

Quantum phase diagram of spin-1 J1−J2 Heisenberg model on the square lattice: An infinite projected entangled-pair state and density matrix renormalization group study

Reza Haghshenas, Wang-Wei Lan, Shou-Shu Gong +1 more

We study the spin-1 Heisenberg model on the square lattice with the antiferromagnetic nearest-neighbor ${J}_{1}$ and the next-nearest-neighbor ${J}_{2}$ couplings by using the infinite projected entangled-pair state (iPEPS) ansatz and the density matrix renormalization group (DMRG) calculation. The iPEPS simulation, which studies the model directly in the thermodynamic limit, finds a crossing of the ground state from the N\'eel magnetic state to the stripe magnetic state at ${J}_{2}/{J}_{1}\ensuremath{\simeq}0.549$, showing a direct phase transition. In the finite-size DMRG calculation on the cylinder geometry up to the cylinder width ${L}_{y}=10$, we find that around the same critical point the N\'eel and the stripe orders are strongly suppressed, which implies the absence of an intermediate phase. Both calculations identify that the stripe order comes with a first-order transition at ${J}_{2}/{J}_{1}\ensuremath{\simeq}0.549$. Our results indicate that unlike the spin-1/2 ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ square model, quantum fluctuations in the spin-1 model may not be strong enough to stabilize an intermediate nonmagnetic phase.

arXiv:1802.00874•Physical review. B./Physical review. B•May 31, 2018•21 citations

J1−J2 frustrated two-dimensional Heisenberg model: Random phase approximation and functional renormalization group

Johannes Reuther, P. Wölfle

We study the ground-state properties of the two-dimensional spin-1/2 ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model on a square lattice, within diagrammatic approximations using an auxiliary-fermion formulation with exact projection. In a first approximation, we assume a phenomenological width of the pseudofermion spectral function to calculate the magnetization, susceptibilities, and the spin-correlation length within random-phase approximation, demonstrating the appearance of a paramagnetic phase between the N\'eel-ordered and Collinear-ordered phases, at sufficiently large pseudofermion damping. Second we use a functional renormalization-group formulation. We find that the conventional truncation scheme omitting three-particle and higher-order vertices is not sufficient. We therefore include self-energy renormalizations in the single-scale propagator as recently proposed by Katanin, to preserve Ward identities in a better way. We find N\'eel order at $g={J}_{2}/{J}_{1}\ensuremath{\lesssim}{g}_{\text{c}1}\ensuremath{\approx}0.4\dots{}0.45$ and Collinear order at $g\ensuremath{\gtrsim}{g}_{\text{c}2}\ensuremath{\approx}0.66\dots{}0.68$, which is in good agreement with results obtained by numerical studies. In the intervening quantum paramagnetic phase, we find enhanced columnar dimer and plaquette fluctuations of equal strength.

arXiv:0912.0860•Physical Review B•Apr 8, 2010•161 citations

J_1-J_2 frustrated 2d Heisenberg model: RPA and Functional RG

Johannes Reuther, P. Wölfle

We study the ground state properties of the two-dimensional spin-1/2 J_1-J_2-Heisenberg model on a square lattice, within diagrammatic approximations using an auxiliary fermion formulation with exact projection. In a first approximation we assume a phenomenological width of the pseudofermion spectral function to calculate the magnetization, susceptibilities and the spin correlation length within RPA, demonstrating the appearance of a paramagnetic phase between the Neel ordered and Collinear ordered phases, at sufficiently large pseudo fermion damping. Secondly we use a Functional Renormalization Group formulation. We find that the conventional truncation scheme omitting three-particle and higher order vertices is not sufficient. We therefore include self-energy renormalizations in the single-scale propagator as recently proposed by Katanin, to preserve Ward identities in a better way. We find Neel order at g = J_2/J_1 0.66 ... 0.68, which is in good agreement with results obtained by numerical studies. In the intervening quantum paramagnetic phase we find enhanced columnar dimer and plaquette fluctuations of equal strength.

•Dec 4, 2009•0 citations
We present a detailed functional renormalization group analysis of spin-12\frac{1}{2}21​ dipolar Heisenberg model on square lattice. This model is similar to the well-known J1\ensuremath-J2{J}_{1}\text{\ensuremath{-}}{J}_{2}J1​\ensuremath-J2​ model and describes the pseudospin degrees of freedom of polar molecules confined in deep optical lattice with long-range anisotropic dipole-dipole interactions. Previous study of this model based on tensor network ansatz indicates a paramagnetic ground state for certain dipole tilting angles which can be tuned in experiments to control the exchange couplings. The tensor ansatz formulated on a small cluster unit cell is inadequate to describe the spiral order, and therefore the phase diagram at high azimuthal tilting angles remains undetermined. Here, we obtain the full phase diagram of the model from numerical pseudofermion functional renormalization group calculations. We show that an extended quantum paramagnetic phase is realized between the N'eel and stripe/spiral phases. In this region, the spin susceptibility flows smoothly down to the lowest numerical renormalization group scales with no sign of divergence or breakdown of the flow, in sharp contrast to the flow towards the long-range-ordered phases. Our results provide further evidence that the dipolar Heisenberg model is a fertile ground for quantum spin liquids.
The antiferromagnetic Heisenberg model on the triangular lattice is perhaps the best known example of frustrated magnets, but it orders at low temperatures. Recent density matrix renormalization group (DMRG) calculations find that the next nearest neighbor interaction J_{2} enhances the frustration, and it leads to a spin liquid for J_{2}/J_{1}∈(0.08,0.15). In addition, a DMRG study of a dipolar Heisenberg model with longer range interactions gives evidence for a spin liquid at a small dipole tilting angle θ∈[0,10°). In both cases, the putative spin liquid region appears to be small. Here, we show that for the triangular lattice dipolar Heisenberg model, a robust quantum paramagnetic phase exists in a surprisingly wide region, θ∈[0,54°), for dipoles tilted along the lattice diagonal direction. We obtain the phase diagram of the model by functional renormalization group (RG), which treats all magnetic instabilities on equal footing. The quantum paramagnetic phase is characterized by a smooth continuous flow of vertex functions and spin susceptibility down to the lowest RG scale, in contrast to the apparent breakdown of RG flow in phases with stripe or spiral order. Our finding points to a promising direction to search for quantum spin liquids in ultracold dipolar molecules.
Motivated by the experimental realization of quantum spin models of polar molecule KRb in optical lattices, we analyze the spin 1/2 dipolar Heisenberg model with competing anisotropic, long-range exchange interactions. We show that, by tilting the orientation of dipoles using an external electric field, the dipolar spin system on square lattice comes close to a maximally frustrated region similar, but not identical, to that of the J_{1}-J_{2} model. This provides a simple yet powerful route to potentially realize a quantum spin liquid without the need for a triangular or kagome lattice. The ground state phase diagrams obtained from Schwinger-boson and spin-wave theories consistently show a spin disordered region between the Néel, stripe, and spiral phase. The existence of a finite quantum paramagnetic region is further confirmed by an unbiased variational ansatz based on tensor network states and a tensor renormalization group.
Evidence of valence bond solid and quantum paramagnetism in the J1J_1J1​--J2J_2J2​--J3J_3J3​ Heisenberg model from pseudofermion functional renormalization group
Motivated by recent tensor network studies reporting a quantum spin liquid (QSL) within the quantum paramagnetic (QP) regime, we investigate numerically the zero-temperature phase diagram of the J1J_1J1​--J2J_2J2​--J3J_3J3​ Heisenberg model on the square lattice using the pseudofermion functional renormalization group (pf-fRG). Thanks to improved frequency resolution and larger lattice truncations, we are able to explore a significantly larger parameter space than previous pf-fRG work to map out the phase boundary of an extended QP region. We further resolve possible bond order within the QP regime by measuring the dimer response to small deformations in J1J_1J1​, i.e. bond anisotropy that breaks the point-group symmetry of the Hamiltonian. The functional renormalization group flow of the dimer susceptibility suggests that the QP window splits into two: a region where it grows continuously down to the lowest energy scale to signal valence bond solid (VBS) order, and another region where it reaches a peak at a finite scale and is subsequently suppressed to indicate a QSL.
We study the spin-1 Heisenberg model on the square lattice with the antiferromagnetic nearest-neighbor J1{J}_{1}J1​ and the next-nearest-neighbor J2{J}_{2}J2​ couplings by using the infinite projected entangled-pair state (iPEPS) ansatz and the density matrix renormalization group (DMRG) calculation. The iPEPS simulation, which studies the model directly in the thermodynamic limit, finds a crossing of the ground state from the N\'eel magnetic state to the stripe magnetic state at ${J}_{2}/{J}_{1}\ensuremath{\simeq}0.549$, showing a direct phase transition. In the finite-size DMRG calculation on the cylinder geometry up to the cylinder width Ly=10{L}_{y}=10Ly​=10, we find that around the same critical point the N\'eel and the stripe orders are strongly suppressed, which implies the absence of an intermediate phase. Both calculations identify that the stripe order comes with a first-order transition at ${J}_{2}/{J}_{1}\ensuremath{\simeq}0.549$. Our results indicate that unlike the spin-1/2 ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ square model, quantum fluctuations in the spin-1 model may not be strong enough to stabilize an intermediate nonmagnetic phase.
We study the ground-state properties of the two-dimensional spin-1/2 ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model on a square lattice, within diagrammatic approximations using an auxiliary-fermion formulation with exact projection. In a first approximation, we assume a phenomenological width of the pseudofermion spectral function to calculate the magnetization, susceptibilities, and the spin-correlation length within random-phase approximation, demonstrating the appearance of a paramagnetic phase between the N\'eel-ordered and Collinear-ordered phases, at sufficiently large pseudofermion damping. Second we use a functional renormalization-group formulation. We find that the conventional truncation scheme omitting three-particle and higher-order vertices is not sufficient. We therefore include self-energy renormalizations in the single-scale propagator as recently proposed by Katanin, to preserve Ward identities in a better way. We find N\'eel order at $g={J}_{2}/{J}_{1}\ensuremath{\lesssim}{g}_{\text{c}1}\ensuremath{\approx}0.4\dots{}0.45$ and Collinear order at $g\ensuremath{\gtrsim}{g}_{\text{c}2}\ensuremath{\approx}0.66\dots{}0.68$, which is in good agreement with results obtained by numerical studies. In the intervening quantum paramagnetic phase, we find enhanced columnar dimer and plaquette fluctuations of equal strength.
We study the ground state properties of the two-dimensional spin-1/2 J_1-J_2-Heisenberg model on a square lattice, within diagrammatic approximations using an auxiliary fermion formulation with exact projection. In a first approximation we assume a phenomenological width of the pseudofermion spectral function to calculate the magnetization, susceptibilities and the spin correlation length within RPA, demonstrating the appearance of a paramagnetic phase between the Neel ordered and Collinear ordered phases, at sufficiently large pseudo fermion damping. Secondly we use a Functional Renormalization Group formulation. We find that the conventional truncation scheme omitting three-particle and higher order vertices is not sufficient. We therefore include self-energy renormalizations in the single-scale propagator as recently proposed by Katanin, to preserve Ward identities in a better way. We find Neel order at g = J_2/J_1 0.66 ... 0.68, which is in good agreement with results obtained by numerical studies. In the intervening quantum paramagnetic phase we find enhanced columnar dimer and plaquette fluctuations of equal strength.