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

Magnetization process and ordering of the S=12 pyrochlore Heisenberg antiferromagnet in a magnetic field

Imre Hagymási, Robin Schäfer, Roderich Moessner, David J. Luitz

Published Aug 31, 2022• 7 citations
arXivHTMLPDF

Abstract

We study the $S=\frac{1}{2}$ pyrochlore Heisenberg antiferromagnet in a magnetic field. Using large-scale density-matrix renormalization group calculations for clusters with up to 128 spins, we find indications of a finite triplet gap, causing a threshold field to nonzero magnetization in the magnetization curve. We obtain a robust saturation field consistent with a magnon crystal, although the corresponding $5/6$ magnetization plateau is very slim and possibly unstable. Most remarkably, there is a pronounced and apparently robust $1/2$ magnetization plateau where the ground state breaks the rotational symmetry of the lattice, exhibiting oppositely polarized spins on alternating kagome and triangular planes. Reminiscent of the kagome ice plateau of the pyrochlore Ising antiferromagnet known as spin ice, it arises via a much more subtle ``quantum order by disorder'' mechanism.

DOI
10.1103/physrevb.106.l060411

Related papers

A Tale of Two Plateaus: Competing Orders in Spin-1 and Spin-$\tfrac{3}{2}$ Pyrochlore Magnets

Imre Hagymási

We use large-scale density-matrix renormalization group simulations with bond dimensions up to $20\ 000$ to determine the magnetization curves of spin-1 and spin-$\tfrac{3}{2}$ pyrochlore Heisenberg antiferromagnets. Both models exhibit a robust half-magnetization plateau, and we find that the same 16-site state (quadrupled unit cell) is selected in both cases on the largest 64-site cubic cluster we consider for the plateau state. This contrasts sharply with the effective quantum dimer model prediction which favors the ``R'' state, and demonstrates the breakdown of the perturbative mechanism at the Heisenberg point. These results provide a nonperturbative characterization of field-induced phases in pyrochlore magnets and predictive guidance for spin-1 and spin-$\tfrac{3}{2}$ materials.

arXiv:2602.16661•Feb 18, 2026•0 citations

Effect of competing orders on the half-magnetization plateaus in spin-1 and spin-3/2 pyrochlore magnets

Anonymous

We use large-scale density-matrix renormalization group simulations with bond dimensions up to 20000 to determine the magnetization curves of spin-1 and spin-3/2 pyrochlore Heisenberg antiferromagnets. Both models exhibit a robust half-magnetization plateau. On the largest cubic cluster for which controlled energy extrapolations are feasible, containing 64 sites, we find that a plateau state with a 16-site magnetic unit cell (quadrupled structural unit cell) consistently has lower energy than competing ordering patterns for both spin lengths. This contrasts sharply with the prediction of the effective quantum dimer model and suggests a breakdown of the perturbative mechanism at the Heisenberg point. These results provide a nonperturbative characterization of field-induced phases in pyrochlore magnets and predictive guidance for spin-1 and spin-3/2 materials.

Physical review. B./Physical review. B•Aug 20, 2026•0 citations

Quantum magnetism of the spin-1 kagome-lattice antiferromagnet

Katsuhiro Morita

We investigate the spin-1 kagome-lattice Heisenberg antiferromagnet using large-scale Lanczos diagonalization and the finite-temperature Lanczos method. The zero-temperature magnetization process exhibits plateaus at m = 0 , 1 / 3 , 7 / 9 , and 8 / 9 , where m is the normalized magnetization. The m = 0 plateau is identified as a trimer valence-bond-crystal state, while the high-field plateaus at m = 7 / 9 and 8 / 9 are identified as magnon crystals. In particular, the m = 8 / 9 plateau corresponds to the exact localized-magnon crystal state. A smoothed zero-temperature magnetization curve constructed using the Gaussian-kernel smoothing method indicates magnetization jumps at the lower-field edge of the m = 1 / 3 plateau and at the upper-field edges of the m = 7 / 9 and 8 / 9 plateaus. At finite temperatures, the specific heat exhibits a double-peak structure with peaks around T / J ≃ 0.1 and T / J ≃ 1.1 , and the low-temperature peak may be related to trimer valence-bond-crystal ordering. The finite-temperature magnetization curves show that the m = 1 / 3 plateau remains visible at low temperatures, whereas the high-field plateaus are rapidly smeared out by thermal effects. These results provide benchmark data for thermodynamic and high-field magnetization measurements in candidate spin-1 kagome-lattice materials.

arXiv:2607.03086•Materials Today Quantum•Sep 12, 2026•0 citations

Possible Inversion Symmetry Breaking in the S = 1 / 2 Pyrochlore Heisenberg Magnet

Imre Hagymási, Robin Schäfer, Roderich Moessner +1 more

We address the ground-state properties of the long-standing and much-studied three-dimensional quantum spin liquid candidate, the S=1/2 pyrochlore Heisenberg antiferromagnet. By using SU(2) density-matrix renormalization group (DMRG), we are able to access cluster sizes of up to 128 spins. Our most striking finding is a robust spontaneous inversion symmetry breaking, reflected in an energy density difference between the two sublattices of tetrahedra, familiar as a starting point of earlier perturbative treatments. We also determine the ground-state energy, E_{0}/N_{sites}=-0.490(6)J, by combining extrapolations of DMRG with those of a numerical linked cluster expansion. These findings suggest a scenario in which a finite-temperature spin liquid regime gives way to a symmetry-broken state at low temperatures.

arXiv:2010.03563•Physical Review Letters•Mar 18, 2021•54 citations

Enhanced symmetry-breaking tendencies in the S=1 pyrochlore antiferromagnet

Imre Hagymási, Vincent Noculak, Johannes Reuther

We investigate the ground-state properties of the nearest-neighbor $S=1$ pyrochlore Heisenberg antiferromagnet using two complementary numerical methods, the density-matrix renormalization group (DMRG) and pseudofermion functional renormalization group (PFFRG). Within DMRG, we are able to reliably study clusters with up to 48 spins by keeping 20 000 SU(2) states. The investigated 32-site and 48-site clusters both show indications of a robust ${C}_{3}$ rotation symmetry breaking of the ground-state spin correlations and the 48-site cluster additionally features inversion symmetry breaking. Our PFFRG analysis of various symmetry-breaking perturbations corroborates the findings of either ${C}_{3}$ or a combined ${C}_{3}$/inversion symmetry breaking. Moreover, in both methods the symmetry-breaking tendencies appear to be more pronounced than in the $S=1/2$ system.

Physical review. B./Physical review. B•Dec 19, 2022•16 citations

Theory of the $1/9$ Magnetization Plateau of the Spin-$1/2$ Kagome Heisenberg Antiferromagnet

Pratyay Ghosh, Frédéric Mila

We investigate the $1/9$ magnetization plateau of the spin-$1/2$ kagome Heisenberg antiferromagnet using infinite projected entangled-pair states and identify a spin-density-wave state with a $27$-site unit cell and a lower variational energy than previous approaches. In contrast to higher-magnetization plateaus, which are well described by localized magnons on nonoverlapping hexagons and preserve a sixfold rotational symmetry, the $1/9$ plateau only keeps a threefold rotational symmmetry and consists of a tiling by three types of $9$-site clusters related by $\pm2π/3$ rotations. These conclusions are supported by a complementary cluster mean-field theory which correctly preserves the symmetry of the higher plateaus but leads to the same spontaneous symmetry breaking for the $1/9$ plateau with a qualitatively similar pattern of nine different local magnetizations.

arXiv:2610.34439•Oct 25, 2026•0 citations
We study the S=12S=\frac{1}{2}S=21​ pyrochlore Heisenberg antiferromagnet in a magnetic field. Using large-scale density-matrix renormalization group calculations for clusters with up to 128 spins, we find indications of a finite triplet gap, causing a threshold field to nonzero magnetization in the magnetization curve. We obtain a robust saturation field consistent with a magnon crystal, although the corresponding 5/65/65/6 magnetization plateau is very slim and possibly unstable. Most remarkably, there is a pronounced and apparently robust 1/21/21/2 magnetization plateau where the ground state breaks the rotational symmetry of the lattice, exhibiting oppositely polarized spins on alternating kagome and triangular planes. Reminiscent of the kagome ice plateau of the pyrochlore Ising antiferromagnet known as spin ice, it arises via a much more subtle ``quantum order by disorder'' mechanism.
A Tale of Two Plateaus: Competing Orders in Spin-1 and Spin-32\tfrac{3}{2}23​ Pyrochlore Magnets
We use large-scale density-matrix renormalization group simulations with bond dimensions up to 20 00020\ 00020 000 to determine the magnetization curves of spin-1 and spin-32\tfrac{3}{2}23​ pyrochlore Heisenberg antiferromagnets. Both models exhibit a robust half-magnetization plateau, and we find that the same 16-site state (quadrupled unit cell) is selected in both cases on the largest 64-site cubic cluster we consider for the plateau state. This contrasts sharply with the effective quantum dimer model prediction which favors the ``R'' state, and demonstrates the breakdown of the perturbative mechanism at the Heisenberg point. These results provide a nonperturbative characterization of field-induced phases in pyrochlore magnets and predictive guidance for spin-1 and spin-32\tfrac{3}{2}23​ materials.
We use large-scale density-matrix renormalization group simulations with bond dimensions up to 20000 to determine the magnetization curves of spin-1 and spin-3/2 pyrochlore Heisenberg antiferromagnets. Both models exhibit a robust half-magnetization plateau. On the largest cubic cluster for which controlled energy extrapolations are feasible, containing 64 sites, we find that a plateau state with a 16-site magnetic unit cell (quadrupled structural unit cell) consistently has lower energy than competing ordering patterns for both spin lengths. This contrasts sharply with the prediction of the effective quantum dimer model and suggests a breakdown of the perturbative mechanism at the Heisenberg point. These results provide a nonperturbative characterization of field-induced phases in pyrochlore magnets and predictive guidance for spin-1 and spin-3/2 materials.
We investigate the spin-1 kagome-lattice Heisenberg antiferromagnet using large-scale Lanczos diagonalization and the finite-temperature Lanczos method. The zero-temperature magnetization process exhibits plateaus at m = 0 , 1 / 3 , 7 / 9 , and 8 / 9 , where m is the normalized magnetization. The m = 0 plateau is identified as a trimer valence-bond-crystal state, while the high-field plateaus at m = 7 / 9 and 8 / 9 are identified as magnon crystals. In particular, the m = 8 / 9 plateau corresponds to the exact localized-magnon crystal state. A smoothed zero-temperature magnetization curve constructed using the Gaussian-kernel smoothing method indicates magnetization jumps at the lower-field edge of the m = 1 / 3 plateau and at the upper-field edges of the m = 7 / 9 and 8 / 9 plateaus. At finite temperatures, the specific heat exhibits a double-peak structure with peaks around T / J ≃ 0.1 and T / J ≃ 1.1 , and the low-temperature peak may be related to trimer valence-bond-crystal ordering. The finite-temperature magnetization curves show that the m = 1 / 3 plateau remains visible at low temperatures, whereas the high-field plateaus are rapidly smeared out by thermal effects. These results provide benchmark data for thermodynamic and high-field magnetization measurements in candidate spin-1 kagome-lattice materials.
We address the ground-state properties of the long-standing and much-studied three-dimensional quantum spin liquid candidate, the S=1/2 pyrochlore Heisenberg antiferromagnet. By using SU(2) density-matrix renormalization group (DMRG), we are able to access cluster sizes of up to 128 spins. Our most striking finding is a robust spontaneous inversion symmetry breaking, reflected in an energy density difference between the two sublattices of tetrahedra, familiar as a starting point of earlier perturbative treatments. We also determine the ground-state energy, E_{0}/N_{sites}=-0.490(6)J, by combining extrapolations of DMRG with those of a numerical linked cluster expansion. These findings suggest a scenario in which a finite-temperature spin liquid regime gives way to a symmetry-broken state at low temperatures.
We investigate the ground-state properties of the nearest-neighbor S=1S=1S=1 pyrochlore Heisenberg antiferromagnet using two complementary numerical methods, the density-matrix renormalization group (DMRG) and pseudofermion functional renormalization group (PFFRG). Within DMRG, we are able to reliably study clusters with up to 48 spins by keeping 20 000 SU(2) states. The investigated 32-site and 48-site clusters both show indications of a robust C3{C}_{3}C3​ rotation symmetry breaking of the ground-state spin correlations and the 48-site cluster additionally features inversion symmetry breaking. Our PFFRG analysis of various symmetry-breaking perturbations corroborates the findings of either C3{C}_{3}C3​ or a combined C3{C}_{3}C3​/inversion symmetry breaking. Moreover, in both methods the symmetry-breaking tendencies appear to be more pronounced than in the S=1/2S=1/2S=1/2 system.
Theory of the 1/91/91/9 Magnetization Plateau of the Spin-1/21/21/2 Kagome Heisenberg Antiferromagnet
We investigate the 1/91/91/9 magnetization plateau of the spin-1/21/21/2 kagome Heisenberg antiferromagnet using infinite projected entangled-pair states and identify a spin-density-wave state with a 272727-site unit cell and a lower variational energy than previous approaches. In contrast to higher-magnetization plateaus, which are well described by localized magnons on nonoverlapping hexagons and preserve a sixfold rotational symmetry, the 1/91/91/9 plateau only keeps a threefold rotational symmmetry and consists of a tiling by three types of 999-site clusters related by ±2π/3\pm2π/3±2π/3 rotations. These conclusions are supported by a complementary cluster mean-field theory which correctly preserves the symmetry of the higher plateaus but leads to the same spontaneous symmetry breaking for the 1/91/91/9 plateau with a qualitatively similar pattern of nine different local magnetizations.