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

Classical dynamics of the antiferromagnetic Heisenberg S = 1 2 spin ladder

David Dahlbom, Jinu Thomas, Steven Johnston, Kipton M. Barros, Cristian D. Batista

Published Sep 3, 2024• 8 citations
arXivHTMLPDF

Abstract

We employ a classical limit grounded in SU(4) coherent states to investigate the temperature-dependent dynamical spin structure factor of the $S=1/2$ ladder consisting of weakly coupled dimers. By comparing the outcomes of this classical approximation with density matrix renormalization group and exact diagonalization calculations in finite-size ladders, we demonstrate that the classical dynamics offers an accurate approximation across the entire temperature range when the interdimer coupling is weak and a good approximation in the high-temperature regime even when the interdimer coupling is strong. This agreement is achieved after appropriately rescaling the temperature axis and renormalizing expectation values to satisfy a quantum sum rule, following D. Dahlbom et al. [Phys. Rev. B 109, 014427 (2024)]. We anticipate the method will be particularly effective when applied to 2D and 3D lattices composed of weakly coupled dimers, situations that remain challenging for alternative numerical methods.

DOI
10.1103/physrevb.110.104403

Related papers

Spin dynamics and transport in gapped one-dimensional Heisenberg antiferromagnets at nonzero temperatures

Kedar Damle, Subir Sachdev A

We discuss the theory of the nonzero temperature $(T)$ spin dynamics and transport in one-dimensional Heisenberg antiferromagnets with a gap $\ensuremath{\Delta}$. For $T\ensuremath{\ll}\ensuremath{\Delta}$, we develop a semiclassical picture of thermally excited particles. Multiple inelastic collisions between the particles are crucial, and are described by a two-particle $\mathcal{S}$ matrix which is shown to have a superuniversal form at low momenta. This is established by computations on the $O(3)$ $\ensuremath{\sigma}$ model, and strong- and weak-coupling expansions (the latter using a Majorana fermion representation) for the two-leg $S=1/2$ Heisenberg antiferromagnetic ladder. As an aside, we note that the strong-coupling calculation reveals an $S=1$, two-particle bound state which leads to the presence of a second peak in the $T=0$ inelastic neutron-scattering (INS) cross section for a range of values of momentum transfer. We obtain exact, or numerically exact, universal expressions for the thermal broadening of the quasiparticle peak in the INS cross section, the spin diffusivity, and for the field dependence of the NMR relaxation rate ${1/T}_{1}$ of the effective semiclassical model; these are expected to be asymptotically exact for the quantum antiferromagnets in the limit $T\ensuremath{\ll}\ensuremath{\Delta}$. The results for ${1/T}_{1}$ are compared with the experimental findings of Takigawa et al. [Phys. Rev. Lett. 76, 2173 (1996)] and the agreement is quite good. In the regime $\ensuremath{\Delta}<T<$(a typical microscopic exchange) and we argue that a complementary description in terms of semiclassical waves applies, and give some exact results for the thermodynamics and dynamics.

arXiv:cond-mat/9711014•Physical review. B, Condensed matter•Apr 1, 1998•175 citations

Understanding temperature-dependent SU(3) spin dynamics in the S = 1 antiferromagnet Ba2FeSi2O7

Seung-Hwan Do, Hao Zhang, David Dahlbom +9 more

Abstract Quantum magnets admit more than one classical limit and N-level systems with strong single-ion anisotropy are expected to be described by a classical approximation based on SU(N) coherent states. Here we test this hypothesis by modeling finite temperature inelastic neutron scattering (INS) data of the effective spin-one antiferromagnet Ba2FeSi2O7. The measured dynamic structure factor is calculated with a generalized Landau-Lifshitz dynamics for SU(3) spins. Unlike the traditional classical limit based on SU(2) coherent states, the results obtained with classical SU(3) spins are in good agreement with the measured temperature-dependent spectrum. The SU(3) approach developed here provides a general framework to understand the broad class of materials comprising weakly coupled antiferromagnetic dimers, trimers, or tetramers, and magnets with strong single-ion anisotropy.

arXiv:2205.11770•npj Quantum Materials•Jan 12, 2023•18 citations

The antiferromagnetic cross-coupled spin ladder: Quantum fidelity and tensor networks approach

Xi-Hao Chen, Sam Young Cho, Huan-Qiang Zhou +1 more

We investigate the phase diagram of the cross-coupled Heisenberg spin ladder with antiferromagnetic couplings. For this model, the results for the existence of the columnar dimer phase, which was predicted on the basis of weak coupling field theory renormalization group arguments, have been conflicting. The numerical work on this model has been based on various approaches, including exact diagonalization, series expansions and density-matrix renormalization group calculations. Using the recently-developed tensor network states and groundstate fidelity approach for quantum spin ladders, we find no evidence for the existence of the columnar dimer phase. We also provide an argument based on the symmetry of the Hamiltonian, which suggests that the phase diagram for antiferromagnetic couplings consists of a single line separating the rung-singlet and the Haldane phases.

arXiv:1508.03941•Journal of the Korean Physical Society•May 1, 2016•11 citations

Dimerized phase in the cross-coupled antiferromagnetic spin ladder

Gergely Barcza, Örs Legeza, R. M. Noack +1 more

We revisit the phase diagram of the frustrated $s=1/2$ spin ladder with antiferromagnetic rung and diagonal couplings. In particular, we reexamine the evidence for the columnar dimer phase, which has been predicted from analytic treatment of the model and has been claimed to be found in numerical calculations. By performing density-matrix renormalization-group (DMRG) calculations on long ladders, keeping a large number of states within the DMRG, and carrying out careful scaling in number of states kept and system size, we calculate the behavior of the dimer order parameter as a function of frustration strength. We find no positive numerical evidence for a finite dimer order parameter in the thermodynamic limit anywhere in the parameter regime in which the columnar dimer phase is expected to appear. Conservative error estimates in the scaling of our data place stringent limits on the maximum possible value of the dimer order parameter and the maximum possible extent of the parameter regime in which a sufficiently weak dimer phase could still exist.

arXiv:1104.3990•Physical Review B•Aug 20, 2012•24 citations

Statics and dynamics of weakly coupled antiferromagnetic spin- 1 2 ladders in a magnetic field

Pierre Bouillot, Corinna Kollath, Andreas Martin Läuchli +9 more

We investigate weakly coupled spin-$1/2$ ladders in a magnetic field. The work is motivated by recent experiments on the compound ${(}_{\mathrm{C}}{\mathrm{H}}_{12}{\mathrm{N}}_{)}{\mathrm{CuBr}}_{4}$ (BPCB). We use a combination of numerical and analytical methods, in particular, the density-matrix renormalization group (DMRG) technique, to explore the phase diagram and the excitation spectra of such a system. We give detailed results on the temperature dependence of the magnetization and the specific heat, and the magnetic-field dependence of the nuclear-magnetic-resonance relaxation rate of single ladders. For coupled ladders, treating the weak interladder coupling within a mean-field or quantum Monte Carlo approach, we compute the transition temperature of triplet condensation and its corresponding antiferromagnetic order parameter. Existing experimental measurements are discussed and compared to our theoretical results. Furthermore, we compute, using time-dependent DMRG, the dynamical correlations of a single spin ladder. Our results allow to describe directly the inelastic neutron scattering cross section up to high energies. We focus on the evolution of the spectra with the magnetic field and compare their behavior for different couplings. The characteristic features of the spectra are interpreted using different analytical approaches such as the mapping onto a spin chain, a Luttinger liquid or onto a $t$-$J$ model. For values of parameters for which such measurements exist, we compare our results to inelastic neutron scattering experiments on the compound BPCB and find excellent agreement. We make additional predictions for the high-energy part of the spectrum that are potentially testable in future experiments.

arXiv:1009.0840•Physical Review B•Feb 9, 2011•139 citations

Finite-temperature dynamics of the spin- 12 bond alternating Heisenberg antiferromagnetic chain

H.-J. Mikeska, C. Luckmann

We present results for the dynamic structure factor of the $S=1∕2$ bond alternating Heisenberg chain over a large range of frequencies and temperatures. Data are obtained from a numerical evaluation of thermal averages based on the calculation of all eigenvalues and eigenfunctions for chains of up to 20 spins. Interpretation is guided by the exact temperature dependence in the noninteracting dimer limit which remains qualitatively valid up to an interdimer exchange $\ensuremath{\lambda}\ensuremath{\approx}0.5$. The temperature induced central peak around zero frequency is clearly identified and aspects of the crossover to spin diffusion in its variation from low to high temperatures are discussed. The one-magnon peak acquires an asymmetric shape with increasing temperature. The two-magnon peak is dominated by the $S=1$ bound state which remains well defined up to temperatures of the order of $J$. The variation with temperature and wave vector of the integrated intensity for one-magnon and two-magnon scattering and of the central peak are discussed.

arXiv:cond-mat/0602387•Physical Review B•May 19, 2006•19 citations
We employ a classical limit grounded in SU(4) coherent states to investigate the temperature-dependent dynamical spin structure factor of the S=1/2S=1/2S=1/2 ladder consisting of weakly coupled dimers. By comparing the outcomes of this classical approximation with density matrix renormalization group and exact diagonalization calculations in finite-size ladders, we demonstrate that the classical dynamics offers an accurate approximation across the entire temperature range when the interdimer coupling is weak and a good approximation in the high-temperature regime even when the interdimer coupling is strong. This agreement is achieved after appropriately rescaling the temperature axis and renormalizing expectation values to satisfy a quantum sum rule, following D. Dahlbom et al. [Phys. Rev. B 109, 014427 (2024)]. We anticipate the method will be particularly effective when applied to 2D and 3D lattices composed of weakly coupled dimers, situations that remain challenging for alternative numerical methods.
We discuss the theory of the nonzero temperature (T)(T)(T) spin dynamics and transport in one-dimensional Heisenberg antiferromagnets with a gap $\ensuremath{\Delta}$. For $T\ensuremath{\ll}\ensuremath{\Delta}$, we develop a semiclassical picture of thermally excited particles. Multiple inelastic collisions between the particles are crucial, and are described by a two-particle S\mathcal{S}S matrix which is shown to have a superuniversal form at low momenta. This is established by computations on the O(3)O(3)O(3) $\ensuremath{\sigma}$ model, and strong- and weak-coupling expansions (the latter using a Majorana fermion representation) for the two-leg S=1/2S=1/2S=1/2 Heisenberg antiferromagnetic ladder. As an aside, we note that the strong-coupling calculation reveals an S=1S=1S=1, two-particle bound state which leads to the presence of a second peak in the T=0T=0T=0 inelastic neutron-scattering (INS) cross section for a range of values of momentum transfer. We obtain exact, or numerically exact, universal expressions for the thermal broadening of the quasiparticle peak in the INS cross section, the spin diffusivity, and for the field dependence of the NMR relaxation rate 1/T1{1/T}_{1}1/T1​ of the effective semiclassical model; these are expected to be asymptotically exact for the quantum antiferromagnets in the limit $T\ensuremath{\ll}\ensuremath{\Delta}$. The results for 1/T1{1/T}_{1}1/T1​ are compared with the experimental findings of Takigawa et al. [Phys. Rev. Lett. 76, 2173 (1996)] and the agreement is quite good. In the regime $\ensuremath{\Delta}<T<$(a typical microscopic exchange) and we argue that a complementary description in terms of semiclassical waves applies, and give some exact results for the thermodynamics and dynamics.
Abstract Quantum magnets admit more than one classical limit and N-level systems with strong single-ion anisotropy are expected to be described by a classical approximation based on SU(N) coherent states. Here we test this hypothesis by modeling finite temperature inelastic neutron scattering (INS) data of the effective spin-one antiferromagnet Ba2FeSi2O7. The measured dynamic structure factor is calculated with a generalized Landau-Lifshitz dynamics for SU(3) spins. Unlike the traditional classical limit based on SU(2) coherent states, the results obtained with classical SU(3) spins are in good agreement with the measured temperature-dependent spectrum. The SU(3) approach developed here provides a general framework to understand the broad class of materials comprising weakly coupled antiferromagnetic dimers, trimers, or tetramers, and magnets with strong single-ion anisotropy.
We investigate the phase diagram of the cross-coupled Heisenberg spin ladder with antiferromagnetic couplings. For this model, the results for the existence of the columnar dimer phase, which was predicted on the basis of weak coupling field theory renormalization group arguments, have been conflicting. The numerical work on this model has been based on various approaches, including exact diagonalization, series expansions and density-matrix renormalization group calculations. Using the recently-developed tensor network states and groundstate fidelity approach for quantum spin ladders, we find no evidence for the existence of the columnar dimer phase. We also provide an argument based on the symmetry of the Hamiltonian, which suggests that the phase diagram for antiferromagnetic couplings consists of a single line separating the rung-singlet and the Haldane phases.
We revisit the phase diagram of the frustrated s=1/2s=1/2s=1/2 spin ladder with antiferromagnetic rung and diagonal couplings. In particular, we reexamine the evidence for the columnar dimer phase, which has been predicted from analytic treatment of the model and has been claimed to be found in numerical calculations. By performing density-matrix renormalization-group (DMRG) calculations on long ladders, keeping a large number of states within the DMRG, and carrying out careful scaling in number of states kept and system size, we calculate the behavior of the dimer order parameter as a function of frustration strength. We find no positive numerical evidence for a finite dimer order parameter in the thermodynamic limit anywhere in the parameter regime in which the columnar dimer phase is expected to appear. Conservative error estimates in the scaling of our data place stringent limits on the maximum possible value of the dimer order parameter and the maximum possible extent of the parameter regime in which a sufficiently weak dimer phase could still exist.
We investigate weakly coupled spin-1/21/21/2 ladders in a magnetic field. The work is motivated by recent experiments on the compound (CH12N)CuBr4{(}_{\mathrm{C}}{\mathrm{H}}_{12}{\mathrm{N}}_{)}{\mathrm{CuBr}}_{4}(C​H12​N)​CuBr4​ (BPCB). We use a combination of numerical and analytical methods, in particular, the density-matrix renormalization group (DMRG) technique, to explore the phase diagram and the excitation spectra of such a system. We give detailed results on the temperature dependence of the magnetization and the specific heat, and the magnetic-field dependence of the nuclear-magnetic-resonance relaxation rate of single ladders. For coupled ladders, treating the weak interladder coupling within a mean-field or quantum Monte Carlo approach, we compute the transition temperature of triplet condensation and its corresponding antiferromagnetic order parameter. Existing experimental measurements are discussed and compared to our theoretical results. Furthermore, we compute, using time-dependent DMRG, the dynamical correlations of a single spin ladder. Our results allow to describe directly the inelastic neutron scattering cross section up to high energies. We focus on the evolution of the spectra with the magnetic field and compare their behavior for different couplings. The characteristic features of the spectra are interpreted using different analytical approaches such as the mapping onto a spin chain, a Luttinger liquid or onto a ttt-JJJ model. For values of parameters for which such measurements exist, we compare our results to inelastic neutron scattering experiments on the compound BPCB and find excellent agreement. We make additional predictions for the high-energy part of the spectrum that are potentially testable in future experiments.
We present results for the dynamic structure factor of the S=1∕2S=1∕2S=1∕2 bond alternating Heisenberg chain over a large range of frequencies and temperatures. Data are obtained from a numerical evaluation of thermal averages based on the calculation of all eigenvalues and eigenfunctions for chains of up to 20 spins. Interpretation is guided by the exact temperature dependence in the noninteracting dimer limit which remains qualitatively valid up to an interdimer exchange $\ensuremath{\lambda}\ensuremath{\approx}0.5$. The temperature induced central peak around zero frequency is clearly identified and aspects of the crossover to spin diffusion in its variation from low to high temperatures are discussed. The one-magnon peak acquires an asymmetric shape with increasing temperature. The two-magnon peak is dominated by the S=1S=1S=1 bound state which remains well defined up to temperatures of the order of JJJ. The variation with temperature and wave vector of the integrated intensity for one-magnon and two-magnon scattering and of the central peak are discussed.