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

Spin and thermal transport and critical phenomena in three-dimensional antiferromagnets

Kazushi Aoyama

Published Dec 7, 2022• 4 citations
arXivHTMLPDF

Abstract

We investigate spin and thermal transport near the N\'eel transition temperature ${T}_{N}$ in three dimensions by numerically analyzing the classical antiferromagnetic $XXZ$ model on the cubic lattice, where in the model, the anisotropy of the exchange interaction $\mathrm{\ensuremath{\Delta}}={J}_{z}/{J}_{x}$ plays a role to control the universality class of the transition. It is found by means of the hybrid Monte Carlo and spin-dynamics simulations that in the $XY$ and Heisenberg cases of $\mathrm{\ensuremath{\Delta}}\ensuremath{\le}1$, the longitudinal spin conductivity ${\ensuremath{\sigma}}_{\ensuremath{\mu}\ensuremath{\mu}}^{s}$ exhibits a divergent enhancement on cooling toward ${T}_{N}$ while not in the Ising case of $\mathrm{\ensuremath{\Delta}}>1$. In all three cases, the temperature dependence of the thermal conductivity ${\ensuremath{\kappa}}_{\ensuremath{\mu}\ensuremath{\mu}}$ is featureless at ${T}_{N}$, being consistent with experimental results. The divergent enhancement of ${\ensuremath{\sigma}}_{\ensuremath{\mu}\ensuremath{\mu}}^{s}$ toward ${T}_{N}$ is attributed to the spin-current relaxation time which gets longer toward ${T}_{N}$, showing a power-law divergence characteristic of critical phenomena. It is also found that in contrast to the $XY$ case where the divergence in ${\ensuremath{\sigma}}_{\ensuremath{\mu}\ensuremath{\mu}}^{s}$ is rapidly suppressed below ${T}_{N}, {\ensuremath{\sigma}}_{\ensuremath{\mu}\ensuremath{\mu}}^{s}$ likely remains divergent even below ${T}_{N}$ in the Heisenberg case, which might experimentally be observed in the ideally isotropic antiferromagnet ${\mathrm{RbMnF}}_{3}$.

DOI
10.1103/physrevb.106.224407

Related papers

Effects of magnetic anisotropy on spin and thermal transport in classical antiferromagnets on the square lattice

Kazushi Aoyama, Hikaru Kawamura

Transport properties of the classical antiferromagnetic XXZ model on the square lattice have been theoretically investigated, putting emphasis on how the occurrence of a phase transition is reflected in spin and thermal transports. As is well known, the anisotropy of the exchange interaction $\mathrm{\ensuremath{\Delta}}\ensuremath{\equiv}{J}_{z}/{J}_{x}$ plays a role to control the universality class of the transition of the model, i.e., either a second-order transition at ${T}_{N}$ into a magnetically ordered state or the Kosterlitz-Thouless (KT) transition at ${T}_{\text{KT}}$, which respectively occur for the Ising-type ($\mathrm{\ensuremath{\Delta}}>1$) and XY-type ($\mathrm{\ensuremath{\Delta}}<1$) anisotropies, while for the isotropic Heisenberg case of $\mathrm{\ensuremath{\Delta}}=1$, a phase transition does not occur at any finite temperature. It is found by means of the hybrid Monte Carlo and spin-dynamics simulations that the spin current probes the difference in the ordering properties, while the thermal current does not. For the XY-type anisotropy, the longitudinal spin-current conductivity ${\ensuremath{\sigma}}_{xx}^{s}$ ($={\ensuremath{\sigma}}_{yy}^{s}$) exhibits a divergence at ${T}_{\text{KT}}$ of the exponential form ${\ensuremath{\sigma}}_{xx}^{s}\ensuremath{\propto}exp\left[B/\sqrt{T/{T}_{\text{KT}}\ensuremath{-}1}\phantom{\rule{0.16em}{0ex}}\right]$ with $B=O(1)$, while for the Ising-type anisotropy, the temperature dependence of ${\ensuremath{\sigma}}_{xx}^{s}$ is almost monotonic without showing a clear anomaly at ${T}_{N}$ and such a monotonic behavior is also the case in the Heisenberg-type spin system. The significant enhancement of ${\ensuremath{\sigma}}_{xx}^{s}$ at ${T}_{\text{KT}}$ is found to be due to the exponential rapid growth of the spin-current relaxation time toward ${T}_{\text{KT}}$, which can be understood as a manifestation of the topological nature of a vortex whose lifetime is expected to get longer toward ${T}_{\text{KT}}$. Possible experimental platforms for the spin-transport phenomena associated with the KT topological transition are discussed.

arXiv:1908.06630•Physical review. B./Physical review. B•Oct 9, 2019•18 citations

Critical entropy of quantum Heisenberg magnets on simple-cubic lattices

Stefan Weßel

We analyze the temperature dependence of the entropy of the spin-1/2 Heisenberg model on the three-dimensional simple-cubic lattice, for both the case of antiferromagnetic and ferromagnetic nearest-neighbor exchange interactions. Using optimized extended ensemble quantum Monte Carlo simulations, we extract the entropy at the critical temperature for magnetic order from a finite-size scaling analysis. For the antiferromagnetic case, the critical entropy density equals $0.341(5){k}_{B}$, whereas for the ferromagnet, a larger value of $0.401(5){k}_{B}$ is obtained. We compare our simulation results to estimates put forward recently in studies assessing means of realizing the antiferromagnetic N\'eel state in ultracold fermion gases in optical lattices.

arXiv:0912.3990•Physical Review B•Feb 10, 2010•32 citations

Critical Temperature and the Transition from Quantum to Classical Order Parameter Fluctuations in the Three-Dimensional Heisenberg Antiferromagnet

Anders W. Sandvik

We present results of extensive quantum Monte Carlo simulations of the three-dimensional $S\phantom{\rule{0ex}{0ex}}=\phantom{\rule{0ex}{0ex}}1/2$ Heisenberg antiferromagnet. Finite-size scaling of the spin stiffness and the sublattice magnetization gives the critical temperature ${T}_{\mathrm{c}}/J\phantom{\rule{0ex}{0ex}}=\phantom{\rule{0ex}{0ex}}0.946\ifmmode\pm\else\textpm\fi{}0.001$. The critical behavior is consistent with the classical 3D Heisenberg universality class, as expected. We discuss the general nature of the transition from quantum mechanical to classical (thermal) order parameter fluctuations at a continuous ${T}_{\mathrm{c}}>0$ phase transition.

arXiv:cond-mat/9804234•Physical Review Letters•Jun 8, 1998•108 citations

Dynamical critical phenomena in three-dimensional Heisenberg spin glasses

Mitsuhiro Matsumoto, Koji Hukushima, Hajime Takayama

Spin-glass (SG) and chiral-glass (CG) orderings in three-dimensional (3D) Heisenberg spin glass with and without magnetic anisotropy are studied by using large-scale off-equilibrium Monte Carlo simulations. A characteristic time of relaxation, which diverges at a transition temperature in the thermodynamic limit, is obtained as a function of the temperature and the system size. Based on the finite-size-scaling analysis for the relaxation time, it is found that in the isotropic Heisenberg spin glass, the CG phase transition occurs at a finite temperature, while the SG transition occurs at a lower temperature, which is compatible with zero. Our results of the anisotropic case support the chirality scenario for the phase transitions in the 3D Heisenberg spin glasses.

arXiv:cond-mat/0204225•Physical review. B, Condensed matter•Sep 4, 2002•22 citations

Weak First-Order Transition and Pseudoscaling Behavior in the Universality Class of the O(N) Ising Model

А. О. Сорокин

Using Monte Carlo and renormalization group methods, we investigate systems with critical behavior described by two order parameters: continuous (vector) and discrete (scalar). We consider two models of classical three-dimensional Heisenberg magnets with different numbers of spin components N = 1,…,4: the model on a cubic lattice with an additional competing antiferromagnetic exchange interaction in a layer and the model on a body-centered lattice with two competing antiferromagnetic interactions. In both models, we observe a first-order transition for all values of N. In the case where competing exchanges are approximately equal, the first order of a transition is close to the second order, and pseudoscaling behavior is observed with critical exponents differing from those of the O(N) model. In the case N = 2, the critical exponents are consistent with the well-known indices of the class of magnets with a noncollinear spin ordering. We also give a possible explanation of the observed pseudoscaling in the framework of the renormalization group analysis.

Theoretical and Mathematical Physics•Aug 1, 2019•5 citations

Universality in the three-dimensional random bond quantum Heisenberg antiferromagnet

Ulvi Kanbur, Erol Vatansever, Hamza Polat

The three-dimensional quenched random bond diluted $({J}_{1}\ensuremath{-}{J}_{2})$ quantum Heisenberg antiferromagnet is studied on a simple-cubic lattice. Using extensive stochastic series expansion quantum Monte Carlo simulations, we perform very long runs for an $L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L$ lattice up to $L=48$. By employing a standard finite-size scaling method, the numerical values of the N\'eel temperature are determined with high precision as a function of the coupling ratio $r={J}_{2}/{J}_{1}$. Based on the estimated critical exponents, we find that the critical behavior of the considered model belongs to the pure classical three-dimensional $O(3)$ Heisenberg universality class.

arXiv:2008.01464•Physical review. B./Physical review. B•Aug 14, 2020•6 citations
We investigate spin and thermal transport near the N'eel transition temperature TN{T}_{N}TN​ in three dimensions by numerically analyzing the classical antiferromagnetic XXZXXZXXZ model on the cubic lattice, where in the model, the anisotropy of the exchange interaction \ensuremathΔ=Jz/Jx\mathrm{\ensuremath{\Delta}}={J}_{z}/{J}_{x}\ensuremathΔ=Jz​/Jx​ plays a role to control the universality class of the transition. It is found by means of the hybrid Monte Carlo and spin-dynamics simulations that in the XYXYXY and Heisenberg cases of \ensuremathΔ\ensuremath≤1\mathrm{\ensuremath{\Delta}}\ensuremath{\le}1\ensuremathΔ\ensuremath≤1, the longitudinal spin conductivity \ensuremathσ\ensuremathμ\ensuremathμs{\ensuremath{\sigma}}_{\ensuremath{\mu}\ensuremath{\mu}}^{s}\ensuremathσ\ensuremathμ\ensuremathμs​ exhibits a divergent enhancement on cooling toward TN{T}_{N}TN​ while not in the Ising case of \ensuremathΔ>1\mathrm{\ensuremath{\Delta}}>1\ensuremathΔ>1. In all three cases, the temperature dependence of the thermal conductivity \ensuremathκ\ensuremathμ\ensuremathμ{\ensuremath{\kappa}}_{\ensuremath{\mu}\ensuremath{\mu}}\ensuremathκ\ensuremathμ\ensuremathμ​ is featureless at TN{T}_{N}TN​, being consistent with experimental results. The divergent enhancement of \ensuremathσ\ensuremathμ\ensuremathμs{\ensuremath{\sigma}}_{\ensuremath{\mu}\ensuremath{\mu}}^{s}\ensuremathσ\ensuremathμ\ensuremathμs​ toward TN{T}_{N}TN​ is attributed to the spin-current relaxation time which gets longer toward TN{T}_{N}TN​, showing a power-law divergence characteristic of critical phenomena. It is also found that in contrast to the XYXYXY case where the divergence in \ensuremathσ\ensuremathμ\ensuremathμs{\ensuremath{\sigma}}_{\ensuremath{\mu}\ensuremath{\mu}}^{s}\ensuremathσ\ensuremathμ\ensuremathμs​ is rapidly suppressed below TN,\ensuremathσ\ensuremathμ\ensuremathμs{T}_{N}, {\ensuremath{\sigma}}_{\ensuremath{\mu}\ensuremath{\mu}}^{s}TN​,\ensuremathσ\ensuremathμ\ensuremathμs​ likely remains divergent even below TN{T}_{N}TN​ in the Heisenberg case, which might experimentally be observed in the ideally isotropic antiferromagnet RbMnF3{\mathrm{RbMnF}}_{3}RbMnF3​.
Transport properties of the classical antiferromagnetic XXZ model on the square lattice have been theoretically investigated, putting emphasis on how the occurrence of a phase transition is reflected in spin and thermal transports. As is well known, the anisotropy of the exchange interaction $\mathrm{\ensuremath{\Delta}}\ensuremath{\equiv}{J}_{z}/{J}_{x}$ plays a role to control the universality class of the transition of the model, i.e., either a second-order transition at TN{T}_{N}TN​ into a magnetically ordered state or the Kosterlitz-Thouless (KT) transition at TKT{T}_{\text{KT}}TKT​, which respectively occur for the Ising-type ($\mathrm{\ensuremath{\Delta}}>1$) and XY-type ($\mathrm{\ensuremath{\Delta}}<1$) anisotropies, while for the isotropic Heisenberg case of $\mathrm{\ensuremath{\Delta}}=1$, a phase transition does not occur at any finite temperature. It is found by means of the hybrid Monte Carlo and spin-dynamics simulations that the spin current probes the difference in the ordering properties, while the thermal current does not. For the XY-type anisotropy, the longitudinal spin-current conductivity ${\ensuremath{\sigma}}_{xx}^{s}$ ($={\ensuremath{\sigma}}_{yy}^{s}$) exhibits a divergence at TKT{T}_{\text{KT}}TKT​ of the exponential form ${\ensuremath{\sigma}}_{xx}^{s}\ensuremath{\propto}exp\left[B/\sqrt{T/{T}_{\text{KT}}\ensuremath{-}1}\phantom{\rule{0.16em}{0ex}}\right]$ with B=O(1)B=O(1)B=O(1), while for the Ising-type anisotropy, the temperature dependence of ${\ensuremath{\sigma}}_{xx}^{s}$ is almost monotonic without showing a clear anomaly at TN{T}_{N}TN​ and such a monotonic behavior is also the case in the Heisenberg-type spin system. The significant enhancement of ${\ensuremath{\sigma}}_{xx}^{s}$ at TKT{T}_{\text{KT}}TKT​ is found to be due to the exponential rapid growth of the spin-current relaxation time toward TKT{T}_{\text{KT}}TKT​, which can be understood as a manifestation of the topological nature of a vortex whose lifetime is expected to get longer toward TKT{T}_{\text{KT}}TKT​. Possible experimental platforms for the spin-transport phenomena associated with the KT topological transition are discussed.
We analyze the temperature dependence of the entropy of the spin-1/2 Heisenberg model on the three-dimensional simple-cubic lattice, for both the case of antiferromagnetic and ferromagnetic nearest-neighbor exchange interactions. Using optimized extended ensemble quantum Monte Carlo simulations, we extract the entropy at the critical temperature for magnetic order from a finite-size scaling analysis. For the antiferromagnetic case, the critical entropy density equals 0.341(5)kB0.341(5){k}_{B}0.341(5)kB​, whereas for the ferromagnet, a larger value of 0.401(5)kB0.401(5){k}_{B}0.401(5)kB​ is obtained. We compare our simulation results to estimates put forward recently in studies assessing means of realizing the antiferromagnetic N\'eel state in ultracold fermion gases in optical lattices.
We present results of extensive quantum Monte Carlo simulations of the three-dimensional S=1/2S\phantom{\rule{0ex}{0ex}}=\phantom{\rule{0ex}{0ex}}1/2S=1/2 Heisenberg antiferromagnet. Finite-size scaling of the spin stiffness and the sublattice magnetization gives the critical temperature ${T}_{\mathrm{c}}/J\phantom{\rule{0ex}{0ex}}=\phantom{\rule{0ex}{0ex}}0.946\ifmmode\pm\else\textpm\fi{}0.001$. The critical behavior is consistent with the classical 3D Heisenberg universality class, as expected. We discuss the general nature of the transition from quantum mechanical to classical (thermal) order parameter fluctuations at a continuous Tc>0{T}_{\mathrm{c}}>0Tc​>0 phase transition.
Spin-glass (SG) and chiral-glass (CG) orderings in three-dimensional (3D) Heisenberg spin glass with and without magnetic anisotropy are studied by using large-scale off-equilibrium Monte Carlo simulations. A characteristic time of relaxation, which diverges at a transition temperature in the thermodynamic limit, is obtained as a function of the temperature and the system size. Based on the finite-size-scaling analysis for the relaxation time, it is found that in the isotropic Heisenberg spin glass, the CG phase transition occurs at a finite temperature, while the SG transition occurs at a lower temperature, which is compatible with zero. Our results of the anisotropic case support the chirality scenario for the phase transitions in the 3D Heisenberg spin glasses.
Using Monte Carlo and renormalization group methods, we investigate systems with critical behavior described by two order parameters: continuous (vector) and discrete (scalar). We consider two models of classical three-dimensional Heisenberg magnets with different numbers of spin components N = 1,…,4: the model on a cubic lattice with an additional competing antiferromagnetic exchange interaction in a layer and the model on a body-centered lattice with two competing antiferromagnetic interactions. In both models, we observe a first-order transition for all values of N. In the case where competing exchanges are approximately equal, the first order of a transition is close to the second order, and pseudoscaling behavior is observed with critical exponents differing from those of the O(N) model. In the case N = 2, the critical exponents are consistent with the well-known indices of the class of magnets with a noncollinear spin ordering. We also give a possible explanation of the observed pseudoscaling in the framework of the renormalization group analysis.
The three-dimensional quenched random bond diluted $({J}_{1}\ensuremath{-}{J}_{2})$ quantum Heisenberg antiferromagnet is studied on a simple-cubic lattice. Using extensive stochastic series expansion quantum Monte Carlo simulations, we perform very long runs for an $L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L$ lattice up to L=48L=48L=48. By employing a standard finite-size scaling method, the numerical values of the N\'eel temperature are determined with high precision as a function of the coupling ratio r=J2/J1r={J}_{2}/{J}_{1}r=J2​/J1​. Based on the estimated critical exponents, we find that the critical behavior of the considered model belongs to the pure classical three-dimensional O(3)O(3)O(3) Heisenberg universality class.