pith. sign in

arxiv: 2601.20328 · v3 · pith:ZKKDTKVZnew · submitted 2026-01-28 · ❄️ cond-mat.str-el

Ground-State Phase Diagram of (1/2,1/2,1) Mixed Diamond Chains with Single-Site Anisotropy

Kazuo Hida This is my paper

Pith reviewed 2026-05-25 06:52 UTC · model grok-4.3

classification ❄️ cond-mat.str-el
keywords mixed diamond chainsingle-site anisotropyground-state phase diagramTomonaga-Luttinger liquidferrimagnetic phasesNéel orderspin chainsexact diagonalization
0
0 comments X

The pith

Mixed diamond chains with (1/2,1/2,1) spins exhibit anisotropy inversion where easy-plane and easy-axis behaviors swap their phase types.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper determines the ground-state phases of mixed diamond chains with vertex spins of magnitude 1/2 and apical spins of 1/2 and 1, including single-site anisotropy D acting on the S=1 site. Numerical exact diagonalization and DMRG calculations, supplemented by analytical approximations in limiting cases, map out the phases as a function of anisotropy and exchange parameters. The central result is the identification of a Néel ordered phase, a nonmagnetic Tomonaga-Luttinger liquid phase, and quantized and partial ferrimagnetic phases, plus a region of anisotropy inversion. In this region the Ising-like Néel phase occurs for easy-plane anisotropy D greater than zero while the XY-like Tomonaga-Luttinger liquid occurs for easy-axis anisotropy D less than zero. A sympathetic reader would care because the inversion reverses the usual link between anisotropy sign and ordering type in these one-dimensional frustrated magnets.

Core claim

The ground-state phase diagram of the (1/2,1/2,1) mixed diamond chain with single-site anisotropy D consists of a Néel ordered phase, a nonmagnetic Tomonaga-Luttinger liquid phase, and quantized and partial ferrimagnetic phases. A region with anisotropy inversion is found where the Ising-like Néel phase is realized for the easy-plane anisotropy D >0 and the XY-like Tomonaga-Luttinger liquid phase is realized for the easy-axis anisotropy D <0 on the S=1 sites.

What carries the argument

The (1/2,1/2,1) mixed diamond chain with vertex-apical exchange couplings modified by parameters delta and lambda, plus single-site anisotropy D on the S=1 apical spin; this lattice structure and its Hamiltonian enable the phase identification through finite-size numerics and limiting-case analysis.

If this is right

  • The Néel phase is Ising-like for D>0 inside the inversion region.
  • The Tomonaga-Luttinger liquid is XY-like for D<0 inside the inversion region.
  • Quantized ferrimagnetic phases appear at specific values of magnetization per unit cell.
  • Partial ferrimagnetic phases occupy additional regions of the diagram.
  • Phase boundaries obtained numerically match analytical results in various limiting cases.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The inversion region may be located by scanning the sign of D while monitoring the character of correlations in the Néel and liquid phases.
  • Materials approximating this mixed-spin diamond geometry could be tested for swapped anisotropy responses in magnetization or specific heat.
  • The same numerical protocol could be applied to nearby values of the apical spin magnitudes to check whether the inversion persists.

Load-bearing premise

Finite-size exact diagonalization and DMRG calculations together with analytical approximations correctly identify the thermodynamic-limit phases and the location of the anisotropy inversion boundary.

What would settle it

DMRG results on significantly larger system sizes that show the anisotropy inversion region shrinking to zero width or vanishing entirely.

Figures

Figures reproduced from arXiv: 2601.20328 by Kazuo Hida.

Figure 1
Figure 1. Figure 1: Structure of the diamond chain investigated in this work. 3. Limiting Cases 3.1 λ = 0 The system is unfrustrated. In the isotropic case D = 0, the ground state is the QF phase with spontaneous magnetization msp = 1 per unit cell according to the Lieb-Mattis (LM) theorem.8, 9) In the presence of easy￾axis anisotropy D < 0, this state remains as the QF phase with mz sp = 1. This phase is called the LM1 phase… view at source ↗
Figure 2
Figure 2. Figure 2: Ground-state phase diagrams based on the NED data with D = 1 for (a)L = 4, (b)L = 6, and (c) L = 8. The number indicated for each phase is Mz sp of the ground state. two N´eel ordered states that are related by spin inver￾sion, and the other is their antisymmetric superposition. On the other hand, the spin inversion does not change the spontaneously dimerized state. The latter behavior is not found within … view at source ↗
Figure 4
Figure 4. Figure 4: This suggests that this narrow N´eel phase sur [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
Figure 3
Figure 3. Figure 3: Maximum value of spontaneous magnetization mz sp:max within the PF phase plotted against 1/L for D = 1. 0 0.1 0.2 0 0.2 0.4 0.6 1+δmin λmin 1/L D=1 [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 5
Figure 5. Figure 5: for D = −1. Again, around the center of the phase diagram, there are several ferrimagnetic phases that are expected to form a PF phase in the thermody￾namic limit. The N´eel phase in the shaded region of [PITH_FULL_IMAGE:figures/full_fig_p004_5.png] view at source ↗
read the original abstract

The ground-state phases of mixed diamond chains with ($S, \tau^{(1)}, \tau^{(2)})=(1/2,1/2,1)$, where $S$ is the magnitude of vertex spins, and $\tau^{(1)}$ and $\tau^{(2)}$ are those of apical spins, are investigated with the single-site anisotropy $D$ on the $\tau^{(2)}$-site. The two apical spins in each unit cell are coupled by an exchange coupling $\lambda$. The vertex spins are coupled with the top and bottom apical spins by exchange couplings $1+\delta$ and $1-\delta$, respectively. The ground-state phase diagram is determined using the numerical exact diagonalization and DMRG method in addition to the analytical approximations in various limiting cases. The phase diagram consists of a N\'eel ordered phase, a nonmagnetic Tomonaga-Luttinger liquid phase, and quantized and partial ferrimagnetic phases. A region with anisotropy inversion is found where the Ising-like N\'eel phase is realized for the easy-plane anisotropy $D >0$ and the XY-like Tomonaga-Luttinger liquid phase is realized for the easy-axis anisotropy $D <0$ on the $S=1$ sites.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The manuscript determines the ground-state phase diagram of the (1/2,1/2,1) mixed diamond chain with single-site anisotropy D acting on the S=1 sites. The Hamiltonian includes tunable inter-apical coupling λ and bond alternation δ. Using exact diagonalization, DMRG, and limiting-case analytic approximations, the authors map out Néel, Tomonaga-Luttinger liquid (TLL), quantized ferrimagnetic, and partial ferrimagnetic phases, and report an “anisotropy inversion” region in which the Néel phase (Ising-like) appears for D>0 while the TLL phase (XY-like) appears for D<0.

Significance. If the reported inversion survives the thermodynamic limit, the result is noteworthy because it inverts the conventional expectation that easy-axis anisotropy stabilizes Ising order and easy-plane anisotropy stabilizes XY order. The work supplies a concrete numerical example in a frustrated mixed-spin geometry and supplies analytic limits that can serve as benchmarks for future studies.

major comments (1)
  1. [Results and Discussion (phase-boundary figures)] The location and even the existence of the anisotropy-inversion boundary are extracted from finite-L ED and DMRG data. No explicit finite-size extrapolation (gap scaling, order-parameter scaling, or Luttinger-parameter flow versus 1/L) is shown for the inversion line itself. Because distinguishing a gapped Néel phase from a gapless TLL and locating the D-value at which their characters invert both require controlled L→∞ extrapolation, the central claim remains vulnerable to residual finite-size drift.
minor comments (2)
  1. [Abstract and §2] The abstract and introduction should state the system sizes used for the DMRG runs and the maximum bond dimension retained, so that readers can immediately assess the numerical effort.
  2. [Model Hamiltonian] Notation for the two apical spins (τ¹, τ²) and the vertex spin S should be introduced once in the Hamiltonian definition and used consistently thereafter.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and constructive feedback. The major comment raises a valid point about finite-size effects on the anisotropy-inversion boundary. We address it below.

read point-by-point responses

  1. Referee: [Results and Discussion (phase-boundary figures)] The location and even the existence of the anisotropy-inversion boundary are extracted from finite-L ED and DMRG data. No explicit finite-size extrapolation (gap scaling, order-parameter scaling, or Luttinger-parameter flow versus 1/L) is shown for the inversion line itself. Because distinguishing a gapped Néel phase from a gapless TLL and locating the D-value at which their characters invert both require controlled L→∞ extrapolation, the central claim remains vulnerable to residual finite-size drift.

    Authors: We agree that the absence of explicit finite-size scaling for the inversion line leaves the central claim open to finite-size drift. While our ED (L≤20) and DMRG (L≤100) data show stable boundaries across the sizes examined, no gap or order-parameter extrapolations versus 1/L were performed specifically for the inversion point. In the revised manuscript we will add (i) DMRG gap scaling for several cuts across the inversion region and (ii) order-parameter scaling to extrapolate the location of the Néel–TLL boundary in the thermodynamic limit. revision: yes

Circularity Check

0 steps flagged

No circularity; phase diagram obtained by direct numerical computation

full rationale

The work determines the ground-state phases via exact diagonalization, DMRG, and limiting-case analytical approximations applied to the given Hamiltonian. No load-bearing step reduces by construction to a fitted quantity extracted from the same data, nor to a self-citation chain, nor to an ansatz smuggled via prior work. The anisotropy-inversion region is reported as an output of those computations, not an input. This is the standard non-circular outcome for a direct numerical survey.

Axiom & Free-Parameter Ledger

3 free parameters · 2 axioms · 0 invented entities

The model is defined by a standard quantum Heisenberg Hamiltonian plus a single-ion anisotropy term; no new entities are introduced. The couplings λ, δ, and D are varied parameters of the Hamiltonian rather than fitted outputs.

free parameters (3)
  • λ
    Exchange coupling between the two apical spins in each unit cell; scanned as a model parameter.
  • δ
    Asymmetry parameter controlling the difference between top and bottom vertex-apical couplings.
  • D
    Single-site anisotropy strength on the τ^(2) site; scanned across positive and negative values.
axioms (2)
  • standard math The system is described by a local spin Hamiltonian with nearest-neighbor exchange and on-site anisotropy terms.
    Invoked in the definition of the mixed diamond chain model.
  • domain assumption Ground states can be classified into Néel, Tomonaga-Luttinger liquid, and ferrimagnetic phases using standard order parameters and correlation functions.
    Used to label the phases found by numerics.

pith-pipeline@v0.9.0 · 5758 in / 1450 out tokens · 30423 ms · 2026-05-25T06:52:19.388977+00:00 · methodology

Review history (2 revisions) →

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

21 extracted references · 21 canonical work pages · 1 internal anchor

  1. [1]

    Introduction Various exotic quantum phases emerge in low- dimensional frustrated quantum magnets as a result of the interplay of quantum fluctuation and frustration. 1, 2) The conventional diamond chain3, 4) is one of the simplest examples of such systems that can be analyzed thanks to the infinite number of local conservation laws ana- lytically. Extensive...

    work page

  2. [2]

    Ground-State Phase Diagram of (1/2,1/2,1) Mixed Diamond Chains with Single-Site Anisotropy

    Model We investigate the ground-state phases of (1 / 2, 1/ 2, 1) mixed diamond chains with the single-site anisotropy de- scribed by the following Hamiltonian: H = L∑ l=1 [ (1 + δ)Sl(τ (1) l + τ (1) l− 1) + (1 − δ)Sl(τ (2) l + τ (2) l− 1) + λτ (1) l τ (2) l + Dτ (2)z2 l ] , (1) where Sl, τ (1) l and τ (2) l are spin operators with magni- tudes Sl = τ(1) l...

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  3. [3]

    In the isotropic case D = 0, the ground state is the QF phase with spontaneous magnetization msp = 1 per unit cell according to the Lieb-Mattis (LM) theorem

    Limiting Cases 3.1 λ = 0 The system is unfrustrated. In the isotropic case D = 0, the ground state is the QF phase with spontaneous magnetization msp = 1 per unit cell according to the Lieb-Mattis (LM) theorem. 8, 9) In the presence of easy- axis anisotropy D < 0, this state remains as the QF phase with mz sp = 1. This phase is called the LM1 phase. On th...

    work page

  4. [4]

    F and AF stand for ferromagnet and antiferro- magnet, respectively

    6(1 + 2 D/ 5λ ). F and AF stand for ferromagnet and antiferro- magnet, respectively. alized for easy-plane anisotropy D > 0 on the S = 1 sites and the XY-like TLL phase is realized for easy- axis anisotropy D < 0 on the S = 1 sites in the range δc > δ > 0. This implies that the effect of D is inverted in this regime. This is also confirmed numerically in th...

    work page

  5. [5]

    2 for D = 1

    Numerical Results 4.1 Case of easy plane anisotropy : D = 1 We have carried out the numerical exact diagonaliza- tion (NED) calculation for L = 4, 6, and 8 to obtain the ground-state phase diagram of Fig. 2 for D = 1. Around the center of the phase diagram, there are sev- eral ferrimagnetic phases. In the thermodynamic limit, these phases are expected to ...

    work page

  6. [6]

    The phase diagram consists of a N´ eel ordered phase, a nonmagnetic TLL phase, QF and PF phases

    Summary and Discussion The ground-state phase diagram of mixed diamond chains with the Hamiltonian (1) is determined using the NED and the DMRG method in addition to the analyt- ical approximations in various limiting cases. The phase diagram consists of a N´ eel ordered phase, a nonmagnetic TLL phase, QF and PF phases. A region with anisotropy inversion ...

    work page

  7. [7]

    Introduction to Frustrated Magnetism: Materials, Experi- ments, Theory , ed. C. Lacroix, P. Mendels, and F. Mila (Springer Series in Solid-State Sciences, Springer, Heide lberg, 2011)

    work page 2011

  8. [8]

    Frustrated Spin Systems , ed. H. T. Diep, (W orld Scientific, Singapore, 2013) 2nd ed

    work page 2013

  9. [9]

    Takano, K

    K. Takano, K. Kubo, and H. Sakamoto, J. Phys.: Condens. Matter 8, 6405 (1996)

    work page 1996

  10. [10]

    Hida and K

    K. Hida and K. Takano, J. Phys. Soc. Jpn. 86, 033707 (2017)

    work page 2017

  11. [11]

    Kikuchi, Y

    H. Kikuchi, Y. Fujii, M. Chiba, S. Mitsudo, T. Idehara, T. Tonegawa, K. Okamoto, T. Sakai, T. Kuwai, and H. Ohta, Phys. Rev. Lett. 94, 227201 (2005)

    work page 2005

  12. [12]

    Kikuchi, Y

    H. Kikuchi, Y. Fujii, M. Chiba, S. Mitsudo, T. Idehara, T. Tone- gawa, K. Okamoto, T. Sakai, T. Kuwai T, K. Kindo, A. Mat- suo, W. Higemoto, K. Nishiyama, M. Horovi´ c, and C. Bertheir, Prog. Theor. Phys. Suppl. 159, 1 (2005)

    work page 2005

  13. [13]

    K. Hida, J. Phys. Soc. Jpn. 90, 054701 (2021)

    work page 2021

  14. [14]

    K. Hida, J. Phys. Soc. Jpn. 93, 044703 (2024)

    work page 2024

  15. [15]

    Lieb and D

    E. Lieb and D. Mattis, J. Math. Phys. 3, 749 (1962)

    work page 1962

  16. [16]

    C. N. Yang and C. P. Yang : Phys. Rev. 151, 258 (1966)

    work page 1966

  17. [17]

    Mikeska and A

    H-J. Mikeska and A. K. Kolezhuk, in Quantum Magnetism , ed. U. Schollw¨ ock, J. Richter, D. J. J. Farnell, R. F. Bishop (Lecture Notes in Physics, Springer, Berlin Heidelberg, 2004).Chap .1.2

    work page 2004

  18. [18]

    Okamoto, Proceedings of the 12th Asia Pacific Physics Conference (APPC12), JPS Conf

    K. Okamoto, Proceedings of the 12th Asia Pacific Physics Conference (APPC12), JPS Conf. Proc. 012031 (2014)

    work page 2014

  19. [19]

    K Okamoto and Y Ichikawa, J. Phys. Chem. Solids 63 1575 (2002)

    work page 2002

  20. [20]

    Tokuno and K

    A. Tokuno and K. Okamoto, J. Phys. Soc. Jpn. Suppl. 74 157 (2005)

    work page 2005

  21. [21]

    Okamoto, Prog

    K. Okamoto, Prog. Theor. Phys. Suppl. 145 208 (2002). 5

    work page 2002