Pith. sign in

REVIEW 3 major objections 5 minor 45 references

Multicritical behavior of the fidelity susceptibility for the 2D quantum transverse-field $XY$ model

T0 review · 3 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The crossover exponent at the 2D XY multicritical point is $\phi = 1.0(2)$, so the Ising-branch phase boundary rises linearly from $(h,\eta)=(2,0)$.

desk verdict A careful ED fidelity-susceptibility study that gives a plausible but not rigorous first estimate of the crossover exponent φ=1.0(2) for the 2D quantum XY model. read the letter →

arxiv 1908.01248 v1 pith:PDI3XWRD submitted 2019-08-04 cond-mat.stat-mech

Yoshihiro Nishiyama (Okayama University) This is my paper · ORCID
classification cond-mat.stat-mech
keywords two-dimensionalquantumXYmodelfidelitysusceptibilitycrossovercriticalexponentmulticriticalpointexactdiagonalizationfinite-sizescalingIsinguniversalityclassphasetransition
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper sets out to determine how the Ising-universality branch of the two-dimensional quantum transverse-field $XY$ model ends at the multicritical point $(h,\eta)=(2,0)$, where the transverse field and the $XY$-plane anisotropy meet. Using the fidelity susceptibility $\chi_F$ as an order-parameter-free probe and exact diagonalization of clusters up to $6\times 6$, it argues that the data collapse onto the crossover scaling form only when the crossover exponent is $\phi=1.0(2)$. If correct, this means the phase boundary obeys $h_c(\eta)-h_c(0)\sim \eta^{1/\phi}\sim \eta$ near the multicritical point, rising at least monotonically with anisotropy. That resolves a discrepancy: the earlier exact-diagonalization data suggested monotonic increase, whereas the spin-anisotropic spherical model predicted a reentrant boundary. A reader should care because the crossover exponent controls the shape of the phase boundary at a point where quantum Monte Carlo is hampered by slowing down, and the paper offers a comparatively cheap exact-diagonalization route to it.

What carries the argument

The load-bearing object is the fidelity susceptibility $\chi_F=(1/N)\,\partial^2_{\Delta h}F|_{\Delta h=0}$ for the ground-state overlap $F=|\langle h|h+\Delta h\rangle|$, which peaks at a quantum phase transition without assuming an order parameter. The argument is carried by the crossover scaling ansatz $\chi_F = L^{\dot x}\,g\big((h-h_c(\eta))L^{1/\dot\nu},\, \eta L^{\phi/\dot\nu}\big)$, with multicritical indices $\dot x=3$ and $\dot\nu=1/2$ fixed by the known $z=2$, $\nu=1/2$ endpoint singularity. Holding the second argument constant at $\eta L^{2\phi}=10.8$ reduces the data to a one-parameter collapse in $\phi$, and the quality of that collapse as $\phi$ is varied is what selects $\phi=1.0(2)$.

What would settle it

Compute the fidelity susceptibility on larger square clusters, for example $L=8$ or $L=10$, at the same scaled anisotropy value $\eta L^{2\phi}=10.8$; if the $L=5$ and $L=6$ overlap does not persist, or a different $\phi$ is needed for the collapse, the $\phi=1.0(2)$ estimate fails. Equivalently, measure the critical field directly for very small $\eta$ and test whether $h_c(\eta)-2$ grows linearly with $\eta$; a sublinear or reentrant approach would contradict the claim.

Watch

Extended reading notes

Core claim

On its own terms, the paper's central claim is that the multicritical point at $(h,\eta)=(2,0)$ has crossover exponent $\phi=1.0(2)$. The evidence is a set of crossover scaling plots of $\chi_F$: with the scaled anisotropy fixed at $\eta L^{2\phi}=10.8$, the $L=4,5,6$ curves overlap for $\phi=1$; postulating $\phi=1.25$ or $\phi=0.75$ produces visible scatter on one side of the peak. Interpreting $h_c-h_c(0)\sim \eta^{1/\phi}$ gives a phase boundary that increases at least monotonically with $\eta$, in agreement with the earlier exact-diagonalization energy-gap data and against the reentrant scenario. As a consistency check, the same machinery at $\eta=1$ recovers $h_c=3.06(2)$ and $\nu=0.614(8)$, compatible with prior estimates.

Load-bearing premise

The estimate rests on assuming that the data for $L=4,5,6$ at the chosen fixed values of $\eta L^{2\phi}$ are already close enough to the infinite-size scaling curve that small-size deviations can be neglected; the paper supports this only by the observed overlap of the $L=5$ and $L=6$ curves, not by a systematic check of corrections.

Editorial extensions

If this is right

  • Near $(h,\eta)=(2,0)$, the Ising-branch phase boundary behaves as $h_c(\eta)-2\sim \eta$, so it leaves the multicritical point linearly, at least monotonically.
  • The $\phi\approx 1$ result backs the earlier exact-diagonalization reading of a monotonic $h_c(\eta)$ and rules out the reentrant phase boundary predicted by the spin-anisotropic spherical model.
  • At the multicritical point itself the fidelity susceptibility diverges with system size as $L^3$, a stronger growth than the $L^{1.259}$ at the Ising transition, making the endpoint easy to locate through $\chi_F$ peaks.
  • Because $\chi_F$ is order-parameter free, the same crossover analysis can be applied on the $XX$-symmetric side of the phase diagram without choosing an order parameter.
  • At $\eta=1$ the method reproduces $h_c=3.06(2)$ and $\nu=0.614(8)$ with modest system sizes, indicating that fidelity-susceptibility scaling gives unbiased criticality estimates even with small clusters.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • My inference: if $\phi=1$ holds to higher precision, the linear boundary $h_c(\eta)-2\sim\eta$ may be exact rather than approximate. A direct small-$\eta$ measurement of the critical field could reveal whether a symmetry protects the linear shape.
  • My inference: the same fixed-$\eta L^{2\phi}$ crossover analysis could be run with quantum Monte Carlo tuned to the $z=2$ endpoint by adjusting the imaginary-time aspect ratio. That would push the test to larger $L$ and check whether the $L=5,6$ data are already asymptotic.
  • My inference: the paper's $\phi=1.0(2)$ suggests that models where reentrant boundaries have been reported, such as spin-$S=1$ $XY$ chains or two-band Hubbard models, deserve the same fidelity-susceptibility crossover analysis. The goal would be to see whether those curved boundaries are genuine or artifacts of the spherical-model approximation.
  • My inference: for experiments, the $L^3$ growth of $\chi_F$ at the multicritical point implies that fidelity-like response functions would show a sharper signal at $(h,\eta)=(2,0)$ than along the Ising branch. This may help locate the endpoint in magnetic or cold-atom systems.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper reports exact-diagonalization calculations of the fidelity susceptibility χ_F for the 2D quantum transverse-field XY model on clusters up to 6×6. After validating the approach at the Ising limit η=1, where the estimated h_c=3.06(2) and ν=0.614(8) are consistent with earlier numerical work, the author performs a crossover-scaling analysis around the multicritical point (h,η)=(2,0). Using the known multicritical exponents x˙=3 and ν˙=1/2 and the scaling ansatz χ_F=L^3 g((h-h_c(η))L^2, ηL^{2φ}), the author tests φ=1, 1.25, and 0.75 by visual data collapse for L=4,5,6 and concludes φ=1.0(2), which implies a linear leading growth h_c(η)-h_c(0)∼η and thus favors a monotonically increasing phase boundary over the reentrant scenario of the spherical model.

Significance. If the estimate φ=1.0(2) is correct, it resolves a discrepancy between earlier exact-diagonalization data and the spin-anisotropic spherical model and anchors the multicritical scaling of the 2D quantum XY model. The paper also demonstrates that the fidelity susceptibility is a useful ED probe for both Ising and XX critical behavior. The Ising-limit results are credible, and the use of independently established multicritical indices avoids circularity in the main scaling analysis. However, the central crossover exponent is inferred from visual collapse of only three lattice sizes at large fixed values of the crossover variable, so the quantitative support for the quoted uncertainty is presently weak and needs to be strengthened before the claim can be accepted as stated.

major comments (3)
  1. [Sec. 2.3, Eq. (10), Figs. 6–8] The estimate φ=1.0(2) rests entirely on visual data collapse for L=4, 5, and 6, with the crossover variable held at u=ηL^{2φ}=10.8 (φ=1), 26.5 (φ=1.25), and 4.4 (φ=0.75). These values of u are of order 10, not asymptotically small, so corrections to the crossover scaling form are uncontrolled; the only evidence that the asymptotic regime has been reached is the statement in Sec. 2.3 that the L=5 and L=6 data 'almost overlap each other, entering at the crossover-scaling regime.' No quantitative collapse residual, no scan in u, and no correction-to-scaling term L^{-ω}g_1(·,u) are provided. As a consequence, the quoted uncertainty 0.2 is not justified; testing only three discrete values of φ cannot bracket a continuous error bar.
  2. [Sec. 2.3, horizontal axis of Figs. 6–8] The estimates h_c(η) used to form the horizontal axis are never reported, and their statistical and systematic uncertainties are not propagated into the collapse. The text states that they were obtained by the same scheme as Sec. 2.1 but using ν=0.63002 [43], whereas Sec. 2.1 uses the author's own ν=0.614(8) for the same Ising branch. This inconsistency needs to be explained, and the sensitivity of the collapse to h_c(η) must be assessed; otherwise the scatter in Figs. 7 and 8 cannot be attributed to φ rather than to inaccurate h_c(η).
  3. [Sec. 2.3, Figs. 6–8] The crossover scaling plots display only negative values of (h-h_c(η))L^2, i.e., only the ordered-side branch; neither the peak region nor the disordered side is shown. In contrast, the η=1 scaling plot in Fig. 5 covers both sides of the transition. A collapse limited to one tail of the scaling function is a much weaker test of Eq. (10), and the paper should either show the full-range data or justify the truncation.
minor comments (5)
  1. [Eq. (2)] The definition is missing a minus sign; since F≈1-(1/2)Nχ_F(Δh)^2 for small Δh, ∂^2 F/∂(Δh)^2|_{Δh=0} is negative, whereas χ_F is used as a positive quantity throughout the paper.
  2. [Fig. 6 caption] The parenthetical '(second argument of the crossover scaling function g (9))' should refer to Eq. (10), not Eq. (9).
  3. [Sec. 2.3 and Sec. 3] The sentence 'our result φ=1.0(2) strongly suggests a linear increase of hc(η) with η' is acceptable for the leading small-η behavior, but the earlier phrasing 'the phase boundary increases, at least, monotonically with η' (Sec. 2.3) does not follow from the scaling analysis and should be softened to a statement about the initial linear growth.
  4. [General] The numerical values of h_c(η,L), the peak values of χ_F, and the collapse curves are not tabulated; making these data available (e.g., as supplementary material) would materially aid verification.
  5. [General] There are a number of typos: 'preceeding' for 'preceding' (Secs. 1 and 2), 'filed' for 'field' in the Hamiltonian paragraph, and 'wherea s' in the introduction; these should be corrected in a revision.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the crossover exponent is estimated from an external scaling ansatz, not reduced to fitted inputs.

full rationale

The paper's central estimate phi = 1.0(2) is obtained by testing three trial values of phi in the crossover scaling variable eta L^(2 phi) and comparing the resulting data-collapse quality of chi_F L^-3 versus (h - h_c(eta)) L^2 (Sec. 2.3, Figs. 6-8). The exponents x_dot = 3 and nu_dot = 1/2 that set these axes are taken from external literature (Refs. [31, 41, 42, 9]), and h_c(eta) is fixed by an independent finite-size extrapolation using nu = 0.63002 from Hasenbusch [43]. Thus the estimated phi is not an input that the collapse is forced to reproduce; it is the variable being discriminated. The Ising-limit parameters nu = 0.614(8) and h_c = 3.06(2) are fitted in Sec. 2.1 and used only for calibration against known results, not to construct the multicritical conclusion. There are no load-bearing self-citations: the cited prior works (Henkel [27], Wald and Henkel [29], Riedel-Wegner [39], Pfeuty et al. [40], etc.) are independent of the present author. Weaknesses the paper itself acknowledges - small system sizes L = 4-6, possible systematic errors, and reliance on visual overlap rather than a quantitative collapse metric - are validity concerns, not circularity. No equation in the paper defines phi in terms of the data collapse, nor is any quantity called a prediction obtained from the same data that fixed it by construction.

Assumptions & free parameters 4 free parameters · 4 assumptions · 0 invented entities

The main free parameter is the crossover exponent phi, which is the central output of the paper rather than an input. The multicritical exponents (x_dot = 3, nu_dot = 1/2) and the Ising nu = 0.63002 are assumed from prior literature and are load-bearing for the scaling collapse. No new entities are introduced. The key additional assumption is that the small-system data are in the asymptotic scaling regime.

free parameters (4)
  • crossover exponent phi = 1.0(2)
    Inferred from the visual quality of the crossover scaling collapse in Figs. 6-8 for phi = 1, 1.25, and 0.75; no quantitative collapse measure is provided.
  • alpha_F/nu at eta=1 = 1.259(38)
    Least-squares fit to three system sizes L = 3, 4, 5 plus a systematic error estimate; used only to validate the method against prior results, not for the main phi estimate.
  • hc at eta=1 = 3.06(2)
    Extrapolated from hc(L) using nu = 0.614; used as a calibration check.
  • nu at eta=1 = 0.614(8)
    Derived from alpha_F/nu via the scaling relation alpha_F/nu = 2/nu - 2; used only for the Ising-limit scaling plot.
assumptions (4)
  • domain assumption At the multicritical point (h, eta) = (2, 0), the fidelity susceptibility and correlation length scale with exponents x_dot = 3 and nu_dot = 1/2, respectively (Eq. (10) and surrounding text).
    These exponents are taken from Refs. [31,41,42] and from the scaling relation using z = 2 and alpha = 1/2; the central crossover collapse assumes them.
  • domain assumption For fixed eta > 0, the phase boundary hc(eta) belongs to the Ising universality class with correlation-length exponent nu ~ 0.63002 [43]; this value is used to extrapolate hc(eta).
    Used in Sec. 2.3 to determine hc; if wrong, the x-axis of the crossover scaling plots shifts.
  • standard math The fidelity susceptibility obeys the finite-size scaling form chi_F = L^x f((h - hc)L^{1/nu}) for Ising transitions and the crossover extension Eq. (10).
    Standard finite-size scaling and crossover scaling theory from Refs. [12,39,40]; accepted background.
  • ad hoc to paper The data for L = 5 and 6 with eta chosen such that eta L^(2 phi) = 10.8 are already in the asymptotic crossover-scaling regime, so corrections to scaling can be ignored.
    Stated in Sec. 2.3: 'the data for L = 5 and 6 almost overlap each other, entering at the crossover-scaling regime.' This is not quantitatively justified.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Multicritical behavior of the fidelity susceptibility for the 2D quantum transverse-field $XY$ model." pith.science (2026). https://pith.science/paper/PDI3XWRD

@misc{pith2026190801248,
  author       = {Pith},
  title        = {Pith review of: Multicritical behavior of the fidelity susceptibility for the 2D quantum transverse-field $XY$ model},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/PDI3XWRD}},
  note         = {Machine review of arXiv:1908.01248}
}
abstract

The two-dimensional quantum $XY$ model with a transverse magnetic field was investigated with the exact diagonalization method. Upon turning on the magnetic field $h$ and the $XY$-plane anisotropy $\eta$, there appear a variety of phase boundaries, which meet at the multicritical point $(h,\eta)=(2,0)$. We devote ourselves to the Ising-universality branch, placing an emphasis on the multicritical behavior. As a probe to detect the underlying phase transitions, we adopt the fidelity susceptibility $\chi_F$. The fidelity susceptibility does not rely on any presumptions as to the order parameter involved. We made a finite-size-scaling analysis of $\chi_F$ for $\eta=1$ (Ising limit), where a number of preceding results are available. Thereby, similar analyses with $\eta$ scaled were carried out around the multicritical point. We found that the $\chi_F$ data are described by the crossover scaling theory. A comparison with the preceding studies of the multicriticality is made.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

45 extracted references · 45 canonical work pages

  1. [43]

    Hasenbusch, Phys

    M. Hasenbusch, Phys. Rev. B 82 (2010) 174433

    work page 2010

  2. [1]

    Uhlmann, Rep

    A. Uhlmann, Rep. Math. Phys. 9 (1976) 273

    work page 1976

  3. [2]

    Jozsa, J

    R. Jozsa, J. Mod. Opt. 41 (1994) 2315

    work page 1994

  4. [3]

    Peres, Phys

    A. Peres, Phys. Rev. A 30 (1984) 1610

    work page 1984

  5. [4]

    Gorin, T

    T. Gorin, T. Prosen, T. H. Seligman, and M. ˇZnidariˇ c, Phys. Rep.435 (2006) 33

    work page 2006

  6. [5]

    H. T. Quan, Z. Song, X. F. Liu, P. Zanardi, and C. P. Sun, Phys . Rev. Lett. 96 (2006) 140604

    work page 2006

  7. [6]

    Zanardi and N

    P. Zanardi and N. Paunkovi´ c, Phys. Rev. E 74 (2006) 031123

    work page 2006

  8. [7]

    Zhou, and J

    H.-Q. Zhou, and J. P. Barjaktarevi˜ c, J. Phys. A: Math. The or. 41 (2008) 412001

    work page 2008

Show all 45 references
  1. [8]

    You and Y.-L

    W.-L. You and Y.-L. Dong, Phys. Rev. B 84 (2011) 174426

    2011

  2. [9]

    Rossini and E

    D. Rossini and E. Vicari, Phys. Rev. E 98 (2018) 062137

    2018

  3. [10]

    V. R. Vieira, J. Phys: Conference Series 213 (2010) 012005

    2010

  4. [11]

    S.-J. Gu, Int. J. Mod. Phys. B 24 (2010) 4371

    2010

  5. [12]

    A. F. Albuquerque, F. Alet, C. Sire, and S. Capponi, Phys. Rev. B 81 (2010) 064418

    2010

  6. [13]

    Schwandt, F

    D. Schwandt, F. Alet, and S. Capponi, Phys. Rev. Lett. 103 (2009) 170501

    2009

  7. [14]

    De Grandi, A

    C. De Grandi, A. Polkovnikov, and A. W. Sandvik, Phys. Rev . B 84 (2011) 224303

    2011

  8. [15]

    W ang, Y.-H

    L. W ang, Y.-H. Liu, J. Imriˇ ska, P. N. Ma, and M. Troyer, Ph ys. Rev. X 5 (2015) 031007

    2015

  9. [16]

    Zhang, X

    J. Zhang, X. Peng, N. Rajendran, and D. Suter, Phys. Rev. L ett. 100 (2008) 100501

    2008

  10. [17]

    Kolodrubetz, V

    M. Kolodrubetz, V. Gritsev, and A. Polkovnikov, Phys. Re v. B 88 (2013) 064304

    2013

  11. [18]

    Gu and W

    S.-J. Gu and W. C. Yu, Europhys. Lett. 108 (2014) 20002

    2014

  12. [19]

    Q. Luo, J. Zhao, and X. W ang, Phys. Rev. E 98 (2018) 022106

    2018

  13. [20]

    Mukherjee, A

    V. Mukherjee, A. Polkovnikov, and A. Dutta, Phys. Rev. B 83 (2011) 075118

    2011

  14. [21]

    Maziero, H

    J. Maziero, H. C. Guzman, L. C. C´ eleri, M. S. Sarandy, and R. M. Serra, Phys. Rev. A 82 (2010) 012106

    2010

  15. [22]

    Sun, Y.-Y

    Z.-Y. Sun, Y.-Y. W u, J. Xu, H.-L. Huang, B.-F. Zhan, B. W an g, and C.-B. Duanpra, Phys. Rev. A 89 (2014) 022101

    2014

  16. [23]

    Karpat, B

    G. Karpat, B. C ¸ akmak, and F. F. Fanchini, Phys. Rev. B 90 (2014) 104431

    2014

  17. [24]

    Katsura, Phys

    S. Katsura, Phys. Rev. 127 (1962) 1508

    1962

  18. [25]

    Barouch, B

    E. Barouch, B. M. McCoy, and M. Dresden, Phys. Rev. A 2 (1970) 1075

    1970

  19. [26]

    Suzuki, Progress Theor

    M. Suzuki, Progress Theor. Phys. 46 (1971) 1337

    1971

  20. [27]

    Henkel, J

    M. Henkel, J. Phys. A: Mathematical and Theoretical 17 (1984) L795

    1984

  21. [28]

    Jalal, R

    S. Jalal, R. Khare, and S. Lal, arXiv:1610.09845

    arXiv

  22. [29]

    W ald and M

    S. W ald and M. Henkel, J. Stat. Mech.: Theory and Experime nt (2015) P07006

    2015

  23. [30]

    V. Z. Kashurnikov, N. V. Prokof’ev, B. V. Scistunov, and M . Troyer, Phys. Rev. B 59 (1999) 1162

    1999

  24. [31]

    V. Zapf, M. Jaime, and C.D. Batista, Rev. Mod. Phys. 86 (2014) 563

    2014

  25. [32]

    M. S. L. du Croo de Jongh and J.M.J. van Leeuwen, Phys. Rev. B 57 (1998) 8494

    1998

  26. [33]

    Yu, H.-M

    W.-C. Yu, H.-M. Kwok, J. Cao and S.-J. Gu, Phys. Rev. E 80 (2009) 021108

    2009

  27. [34]

    Montakhab and A

    A. Montakhab and A. Asadian, Phys. Rev. A 82 (2010) 062313

    2010

  28. [35]

    Huang and F.-L

    C.-Y. Huang and F.-L. Lin, Phys. Rev. A 81 (2010) 032304

    2010

  29. [36]

    Braiorr-Orrs, M

    B. Braiorr-Orrs, M. W eyrauch, and M. V. Rakov, Quantum In formation and Compu- tation 16 (2016) 0885,

    2016

  30. [37]

    Amico, R

    L. Amico, R. Fazio, A. Osterloh, and V. Vedral, Rev. Mod. P hys. 80 (2008) 517

    2008

  31. [38]

    Horodecki, P

    R. Horodecki, P. Horodecki, M. Horodecki, and K. Horodec ki, Rev. Mod. Phys. 81 (2009) 865

    2009

  32. [39]

    Riedel and F

    E.K. Riedel and F. W egner, Z. Phys. 225 (1969) 195

    1969

  33. [40]

    Pfeuty, D

    P. Pfeuty, D. Jasnow, and M. E. Fisher, Phys. Rev. B 10 (1974) 2088

    1974

  34. [41]

    Adamski, J

    M. Adamski, J. J¸ edrzejewski, and T. Krokhmalskii, arXi v:1502.05268

    arXiv

  35. [42]

    Hoeger, G

    C. Hoeger, G. V. Gehlen, and V. Rittenberg, J. Phys. A: Mat h. Gen. 18 (1985) 1813

    1985

  36. [44]

    Hofstetter and M

    W. Hofstetter and M. Henkel, J. Phys. A: Math. Gen. 29 (1996) 1359

    1996

  37. [45]

    De Franco, L

    C. De Franco, L. F. Tocchio and F. Becca, Phys. Rev. B 98 (2018) 075117

    2018

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.