Pith. sign in

REVIEW 3 major objections 5 minor 46 references

Stable diagonal stripes in the t-J model at $\bar{n}_h$=1/8 doping from fPEPS calculations

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

Pith's one-line read Using fermionic projected entangled pair states, this paper finds that the ground state of the t-J model at 1/8 hole doping and J/t=0.4 has stable diagonal stripes, and that long-range superconductivity is suppressed.

desk verdict A credible fPEPS claim of stable diagonal stripes at 1/8 doping with the lowest energy so far, but the orientation claim needs a same-method comparison that isn't in the paper. read the letter →

arxiv 1908.10000 v1 pith:PC27CSK2 submitted 2019-08-27 cond-mat.str-el

classification cond-mat.str-el
keywords t-Jmodel1/8holedopingdiagonalstripesfermionicprojectedentangledpairstatesfPEPSstripeordersuperconductivitysuppressioncorrelations
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 addresses a long-standing dispute over the ground state of the two-dimensional t-J model at 1/8 hole doping, a parameter region relevant to stripe order in copper-oxide superconductors. Using fermionic projected entangled pair states (fPEPS) with gradient optimization and Monte Carlo sampling, it reports a thermodynamic-limit hole energy of -1.6186 at J/t=0.4, lower than previous variational, DMRG, and iPEPS estimates. The central finding is that the ground state has stable stripes running along the lattice diagonal, with hole-density period 4, staggered-magnetization period 8, and a hole filling of rho_l=0.5 per stripe; earlier vertical-stripe results are attributed to boundary conditions that favor the vertical orientation. The same calculation shows that both s-wave and d-wave superconducting pair correlations decay quickly, with power-law exponents of about 4.9 and 4.4, so long-range superconductivity is suppressed in this stripe state. If correct, the paper establishes that the true ground state at this doping is diagonal stripes rather than uniform or vertical-stripe order.

What carries the argument

The load-bearing machinery is the fermionic projected entangled pair state (fPEPS) tensor network, a systematically improvable variational ansatz for interacting fermions on a lattice. The wavefunction is first prepared by imaginary-time evolution with a simple update and then refined by stochastic gradient descent, with energies and gradients estimated by Monte Carlo sampling; U(1) particle-number symmetry is enforced and open boundary conditions are used. With bond dimension D=12 and environment truncation Dc=48, energies on lattices from 4x4 to 12x12 are extrapolated to the thermodynamic limit by a second-order polynomial in 1/sqrt(L1 L2). Because every tensor is optimized independently with no imposed periodicity, the method claims to be unbiased with respect to stripe orientation; the 4x4 result matches exact diagonalization to about 1e-4, validating the larger D=12 runs.

What would settle it

A direct same-method test: on an 8x8 or 12x12 lattice with boundary conditions that do not favor either orientation, optimize fPEPS wavefunctions initialized respectively with diagonal stripes, vertical stripes, and a uniform state, all at the same D=12 and Dc=48; if the vertical-stripe or uniform state has lower energy, or if increasing D to 16 changes the ordering, the central claim fails. Equivalently, a DMRG calculation with periodic boundary conditions along the diagonal direction that finds vertical stripes at equal or lower energy would contradict the orientation claim.

Watch

Extended reading notes

Core claim

The discovery is that the ground state of the t-J model at n_h=1/8, t=1, J/t=0.4, is a diagonal stripe phase rather than the vertical stripe phase found in earlier DMRG and iPEPS calculations. In 12x12 and other open-boundary lattices, the hole density forms site-centered stripes with period 4 along the diagonal direction, while the staggered magnetization has period 8 with a pi phase shift across the domain wall; the stripe hole filling is rho_l=0.5. The authors emphasize that their simulation is unbiased: all L1 x L2 tensors are independent, random initial states converge to the same pattern, and the thermodynamic-limit hole energy (-1.6186) is lower than previously reported values. They also find that s-wave and d-wave pair correlations decay as $r^{{-alpha}}$ with $\alpha$ about 4.9 and 4.4 respectively, so superconductivity does not survive to long range in this phase. The difference in stripe orientation from earlier work is attributed to the different boundary conditions used in those calculations.

Load-bearing premise

The load-bearing premise is that the diagonal-stripe pattern found at bond dimension D=12 on finite open-boundary lattices is the variational ground state, so that no lower-energy vertical-stripe, uniform, or other state was missed; the paper does not report a same-method energy comparison of diagonal versus vertical stripes.

Editorial extensions

If this is right

  • At J/t=0.4 and 1/8 doping, the t-J model ground state is charge-ordered rather than uniform: the hole density has period 4 and the staggered spin period 8, robust across lattice sizes and aspect ratios.
  • The stripe hole filling rho_l=0.5 matches the earlier stripe picture, but the orientation is diagonal; vertical stripes found by DMRG and iPEPS may be artifacts of the anisotropic boundary conditions those methods use.
  • Long-range superconductivity does not coexist with this stripe state: s- and d-wave pair correlations decay with exponents near 4.9 and 4.4, well above the threshold for true long-range order.
  • The reported thermodynamic-limit hole energy of -1.6186 is the lowest among the methods compared, so if the variational energies are converged, the diagonal-stripe state is the best available estimate of the true ground state at this doping.

Reading between the lines

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

  • An implication the authors leave implicit: if diagonal stripes are the true ground state, the widely used vertical-stripe results from DMRG and iPEPS at the same parameters are probably influenced by the boundary conditions, and a direct same-method comparison of both orientations on equivalent lattices would be a decisive check.
  • A testable extension would be to increase the bond dimension D beyond 12 and confirm that the diagonal-stripe hole energy continues to decrease and remains below a constrained vertical-stripe state; if the gap closes, the orientation claim weakens.
  • The paper's suppression of pairing suggests a sharp experimental signature: in a clean two-dimensional t-J-like system at 1/8 doping, one should see diagonal charge modulations and no bulk superconducting phase, which could distinguish this strong-coupling stripe scenario from a uniform resonating-valence-bond picture.
Share X Bluesky LinkedIn Reddit HN

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. This manuscript applies a fermionic projected entangled pair states (fPEPS) method with stochastic gradient descent and Monte Carlo sampling to the two-dimensional t-J model at J/t=0.4 and hole doping n_h=1/8 on finite open-boundary lattices up to 12x12. The authors report a hole energy E_h=-1.6186 after extrapolation to the thermodynamic limit, and they find stripe order along diagonal directions with a hole-density period of 4 and a staggered-magnetization period of 8. They further report that s-wave and d-wave superconducting pair correlations decay with exponents fit to alpha_s≈4.9 and alpha_d≈4.4, from which they conclude that long-range superconductivity is suppressed. The paper compares its hole energy with previous VQMC, iPEPS, and DMRG results and argues that the diagonal orientation is the true ground-state stripe orientation.

Significance. If the diagonal-stripe result is correct, it would be a significant revision of the prevailing vertical-stripe picture for the t-J model at 1/8 doping and would sharpen the debate about stripe order and superconductivity. The paper has genuine strengths: the 4x4 fPEPS energy agrees with exact diagonalization to about 1e-4, the decision to use a larger environment truncation (Dc=6D) for correlation functions is a thoughtful accuracy measure, and the fPEPS optimizations use independent tensors on each site, which is less constrained than the supercell ansatze used in earlier iPEPS calculations. However, the central orientation claim and the superconducting-suppression claim are not yet supported by the evidence actually reported, because no same-method energy comparison against vertical stripes is given and no error estimates accompany the extrapolations and fits.

major comments (3)
  1. [Fig. 1, Table I, and the paragraphs comparing with iPEPS] The central claim that the ground state has diagonal stripes rather than vertical stripes is not supported by a same-method variational comparison. The paper reports fPEPS optimizations that converge to diagonal stripes from random initial states, which demonstrates reproducibility of a single local minimum, not that the diagonal orientation is globally lower in energy. The energy comparisons in Table I are against DMRG and iPEPS results that use different boundary conditions and different ansatz restrictions, so they cannot isolate the stripe orientation. The text explicitly acknowledges that iPEPS calculations found diagonal stripes to have somewhat higher energy than vertical stripes, and it proposes a DMRG calculation with diagonal periodic boundary conditions as a possible resolution without performing it. A direct fPEPS calculation initialized with vertical-stripe states, or a constrained vertical-stripe ansatz with the same D and Dc, is needed to establish that diagonal stripes are lower in energy.
  2. [Fig. 1 and E∞=-0.6701 extrapolation] The thermodynamic limit extrapolation is presented without error bars, fit coefficients, or a documented fit procedure. The reported hole energy E_h=-1.6186 depends on a second-order polynomial fit in sqrt(L1L2) over lattices from 4x4 to 12x12, combined with the zero-doping energy E0 from Ref. 41, but the paper does not state the statistical errors on the individual energies, the sensitivity of E∞ to the fitting form, or the bond-dimension convergence data for the largest systems beyond the assertion that Dc=48 is converged at D=12. Because the competition between stripe and uniform phases is known to be very close in energy, an estimate of the uncertainty in E∞ is essential for claiming the 'most competitive ground state hole energy.'
  3. [Fig. 4 and Eqs. (2)-(3)] The conclusion that long-range superconductivity is suppressed is based on power-law fits with exponents alpha_s≈4.9 and alpha_d≈4.4 extracted from pair correlation functions over distances up to r~10 on a 12x12 lattice with a single reference point r_i=(6,2). No fit residuals, error bars, or alternative decay forms are reported, and the available distance range is short relative to the system size and likely affected by the stripe structure and boundaries. To support the claim, the authors should show log-log fits for multiple reference points and system sizes, or otherwise provide a more controlled finite-size analysis of the correlation decay.
minor comments (5)
  1. [Eq. (1)] In Eq. (1), the phrase 'where ⟨i,j⟩ are the nearest-neighbor sites' should be 'where ⟨i,j⟩ runs over nearest-neighbor pairs of sites' to avoid the impression that each site itself is a pair.
  2. [Abstract and Fig. 3] The abstract describes 'width of 4 unit cells' for the stripes, while the text and Fig. 3 report a hole-density period of 4 and a staggered-magnetization period of 8; 'width' is ambiguous and should be replaced with a precise statement about the charge and spin periodicities.
  3. [Table I] The full-update iPEPS row in Table I corresponds to n_h=0.120, not exactly 1/8; this caveat should appear in the table caption so that the reader does not compare the energies as if they were at identical doping.
  4. [Fig. 3] Figure 3 does not fully label its vertical axes; the green curve should be labeled as the average hole density per site and the red curve as the staggered magnetization (-1)^{i-1}⟨S_i^z⟩, with units indicated.
  5. [Paragraph on robustness against random initial states] The statement that 'we always obtain the same stripe ground states for randomly chosen initial states' would be more informative if the number of random initializations, the range of final energies, and the criterion for identifying the same state were reported.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the diagonal-stripe ground-state claim is an unbiased variational output benchmarked against exact and external results.

full rationale

No circularity is present in the derivation chain. The central claim—that the fPEPS ground state of the t-J model at n_h=1/8 and J/t=0.4 has stable diagonal stripes with stripe filling rho_l=0.5 and suppressed long-range superconductivity—is the output of an unbiased variational optimization, not an input. The paper states that no constraints are imposed on the tensors: 'All L1 x L2 tensors are independent and free to change during the optimization. We always obtain the same stripe ground states for randomly chosen initial states.' The 4x4 fPEPS energy is checked against exact diagonalization (difference about 1e-4), and the thermodynamic-limit hole energy is compared with independent DMRG, iPEPS, and vQMC values. The zero-doping energy E0 used to define E_hole is taken from Sandvik (Ref. 41), an external source, so the hole-energy definition is not self-referential. The only self-citations are to the authors' own fPEPS method papers for optimization details and convergence criteria, but these are not load-bearing in a circular manner: the exact 4x4 agreement and external benchmarks independently support the method. The absence of a same-method vertical-stripe energy comparison is a real limitation, and the paper explicitly acknowledges it by proposing future DMRG calculations with diagonal periodic boundary conditions; however, that is a completeness or correctness risk, not circularity, because the diagonal-stripe result is not derived from an assumption of vertical stripes or from any parameter fitted to the target energy. The stripe filling rho_l = W * n_h = 0.5 is an arithmetic consequence of the observed period-4 hole stripes, not a self-definitional prediction.

Assumptions & free parameters 2 free parameters · 5 assumptions · 0 invented entities

No new physical entities are introduced. The central inputs are the t-J model, the fPEPS method, and an external zero-doping energy. The main free parameters are finite-size and correlation-function fits; the main unproved load-bearing assumptions are D-convergence and boundary-condition neutrality.

free parameters (2)
  • Finite-size scaling coefficients (a and b)
    A second-order polynomial in (L1L2)^-1/2 is fit to five finite-size energies to obtain the thermodynamic-limit energy; the fitted coefficients are not reported.
  • Superconducting correlation decay exponents alpha_s and alpha_d = alpha_s ~ 4.9, alpha_d ~ 4.4
    Power-law fits P_s,d(r) ~ r^{-alpha} to pair correlations on a 12x12 lattice are used to argue that long-range superconductivity is suppressed.
assumptions (5)
  • domain assumption The fPEPS variational ansatz with bond dimension D=12 spans the relevant ground-state manifold at 1/8 doping.
    D is fixed at 12; convergence is checked via Dc and a 4x4 exact comparison, but no systematic D sweep is shown at large lattice sizes.
  • domain assumption The thermodynamic-limit energy is obtained by a second-order polynomial in (L1L2)^-1/2.
    The functional form is a standard finite-size scaling choice, but its validity for these stripe states is not demonstrated in detail.
  • domain assumption Open boundary conditions do not qualitatively change the stripe orientation.
    The paper attributes DMRG vertical stripes to its boundary conditions, but it does not test its own open-boundary geometry against other boundary conditions.
  • domain assumption E0=-0.467775 at zero doping from Sandvik (Ref 41) is accurate.
    This external value is used to define the hole energy and is taken without re-evaluation in this paper.
  • standard math Variational principle: a lower variational energy implies a state closer to the true ground state.
    This is the basis for comparing hole energies across different numerical methods.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Stable diagonal stripes in the t-J model at $\bar{n}_h$=1/8 doping from fPEPS calculations." pith.science (2026). https://pith.science/paper/PC27CSK2

@misc{pith2026190810000,
  author       = {Pith},
  title        = {Pith review of: Stable diagonal stripes in the t-J model at $\barn_h$=1/8 doping from fPEPS calculations},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/PC27CSK2}},
  note         = {Machine review of arXiv:1908.10000}
}
abstract

We investigate the 2D t-J model at a hole doping of $\bar{n}_h$=1/8 using recently developed high accuracy fermionic projected entangled pair states(fPEPS) method. By applying stochastic gradient descent method combined with Monte Carlo sampling technique, we obtain the ground state hole energy $E_{\rm hole}$=-1.6186 for $J/t$=0.4. We show that the ground state has stable diagonal stripes instead of vertical stripes with width of 4 unit cells, and stripe filling $\rho_l$=0.5. We further show that the long range superconductivity order is suppressed at this point.

Figures

Figures reproduced from arXiv: 1908.10000 by the authors.

Figure 1
Figure 1. FIG. 1: The ground state energies of t-J model with [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2: The ground state hole density and spin moment on [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3: The average hole density [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: FIG. 4: The pair correlation function for the d-wave and the s [PITH_FULL_IMAGE:figures/full_fig_p004_4.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

46 extracted references · 33 canonical work pages

  1. [1]

    J. G. Bednorz and K. A. M¨ uller, Z. Phys. B: Condens. Matter 64, 189 (1986)

  2. [2]

    M. K. Wu, J. R. Ashburn, C. J. Torng, P. H. Hor, R. L. Meng, L. Gao, Z. J. Huang, Y. Q. Wang, and C. W. Chu, Phys. Rev. Lett. 58, 908 (1987)

  3. [3]

    +” for s- wave and “−

    The staggered magnetization (red line) shows a period of 8, whereas the hole density shows a period of 4. The site- centered nature of the hole stripes is evident from the hole density ⟨nh i,j ⟩ (green line). These hole and spin pattern with period of 4 and 8, which are robust for different size of systems, can be used to explain why the stripe order canno...

  4. [4]

    Damascelli, Z

    A. Damascelli, Z. Hussain, and Z.-X. Shen, Rev. Mod. Phys. 75, 473 (2003)

  5. [5]

    Valla, A

    T. Valla, A. V. Fedorov, J. Lee, J. C. Davis, and G. D. Gu, Science 314, 1914 (2006), ISSN 0036-8075

  6. [6]

    P. A. Lee, N. Nagaosa, and X.-G. Wen, Rev. Mod. Phys. 78, 17 (2006)

  7. [7]

    Ichikawa, S

    N. Ichikawa, S. Uchida, J. M. Tranquada, T. Niem¨ oller, P. M. Gehring, S.-H. Lee, and J. R. Schneider, Phys. Rev. Lett. 85, 1738 (2000)

  8. [8]

    Hubbard, Proceedings of the Royal Society of Lon- don A: Mathematical, Physical and Engineering Sciences 276, 238 (1963), ISSN 0080-4630

    J. Hubbard, Proceedings of the Royal Society of Lon- don A: Mathematical, Physical and Engineering Sciences 276, 238 (1963), ISSN 0080-4630

Show all 46 references
  1. [9]

    P. W. ANDERSON, Science 235, 1196 (1987), ISSN 0036-8075

  2. [10]

    F. C. Zhang and T. M. Rice, Phys. Rev. B 37, 3759 (1988)

  3. [11]

    Han, Q.-H

    J. Han, Q.-H. Wang, and D.-H. Lee, International Jour- nal of Modern Physics B 15, 1117 (2001)

  4. [12]

    Sherman and M

    A. Sherman and M. Schreiber, European Physical Jour- nal B 32, 203 (2003)

  5. [13]

    Capello, M

    M. Capello, M. Raczkowski, and D. Poilblanc, Phys. Rev. B 77, 224502 (2008)

  6. [14]

    C. S. Hellberg and E. Manousakis, Phys. Rev. Lett. 83, 132 (1999)

  7. [15]

    S. R. White and D. J. Scalapino, Phys. Rev. Lett. 80, 1272 (1998)

  8. [16]

    J. M. Tranquada, J. D. Axe, N. Ichikawa, A. R. Mood- enbaugh, Y. Nakamura, and S. Uchida, Phys. Rev. Lett. 78, 338 (1997)

  9. [17]

    W.-J. Hu, F. Becca, and S. Sorella, Phys. Rev. B 85, 081110(R) (2012)

  10. [18]

    Corboz, S

    P. Corboz, S. R. White, G. Vidal, and M. Troyer, Phys. Rev. B 84, 041108(R) (2011)

  11. [19]

    Corboz, T

    P. Corboz, T. M. Rice, and M. Troyer, Phys. Rev. Lett. 113, 046402 (2014)

  12. [20]

    S. A. Kivelson, I. P. Bindloss, E. Fradkin, V. Oganesyan , J. M. Tranquada, A. Kapitulnik, and C. Howald, Rev. Mod. Phys. 75, 1201 (2003)

  13. [21]

    E. W. Carlson, S. A. Kivelson, D. Orgad, and V. J. Emery, in Physics of Conventional and Unconventional Superconductors, edited by K. H. Bennemann and J. B. Ketterson (Springer, Berlin, 2003)

  14. [22]

    S. R. White and D. J. Scalapino, Phys. Rev. B 61, 6320 (2000)

  15. [23]

    S. R. White and D. J. Scalapino, Phys. Rev. B 70, 220506(R) (2004)

  16. [24]

    S. R. White and D. J. Scalapino, Phys. Rev. B 79, 220504(R) (2009)

  17. [25]

    Schollw¨ ock, Annals of Physics 326, 96 (2011)

    U. Schollw¨ ock, Annals of Physics 326, 96 (2011)

  18. [26]

    Verstraete, V

    F. Verstraete, V. Murg, and J. Cirac, Advances in Physics 57, 143 (2008)

  19. [27]

    H. C. Jiang, Z. Y. Weng, and T. Xiang, Phys. Rev. Lett. 101, 090603 (2008)

  20. [28]

    Verstraete and J

    F. Verstraete and J. I. Cirac, arXiv:cond-mat/0407066 (2004)

  21. [29]

    Sfondrini, J

    A. Sfondrini, J. Cerrillo, N. Schuch, and J. I. Cirac, Ph ys. Rev. B 81, 214426 (2010)

  22. [30]

    Verstraete, M

    F. Verstraete, M. M. Wolf, D. Perez-Garcia, and J. I. Cirac, Phys. Rev. Lett. 96, 220601 (2006)

  23. [31]

    Z.-C. Gu, F. Verstraete, and X.-G. Wen, arXiv:1004.256 3 (2010)

  24. [32]

    Gu, H.-C

    Z.-C. Gu, H.-C. Jiang, D. N. Sheng, H. Yao, L. Balents, and X.-G. Wen, Phys. Rev. B 88, 155112 (2013)

  25. [33]

    Corboz, R

    P. Corboz, R. Or´ us, B. Bauer, and G. Vidal, Phys. Rev. B 81, 165104 (2010)

  26. [34]

    C. V. Kraus, N. Schuch, F. Verstraete, and J. I. Cirac, Phys. Rev. A 81, 052338 (2010)

  27. [35]

    Liu, S.-J

    W.-Y. Liu, S.-J. Dong, Y.-J. Han, G.-C. Guo, and L. He, Phys. Rev. B 95, 195154 (2017)

  28. [36]

    L. He, H. An, C. Yang, F. Wang, J. Chen, C. Wang, W. Liang, S. Dong, Q. Sun, W. Han, et al., IEEE Trans- actions on Parallel and Distributed Systems 29, 2838 (2018)

  29. [37]

    S.-J. Dong, C. Wang, Y. Han, G.-c. Guo, and L. He, Phys. Rev. B 99, 195153 (2019)

  30. [38]

    Normand and A

    B. Normand and A. P. Kampf, Phys. Rev. B 64, 024521 (2001)

  31. [39]

    Barthel, C

    T. Barthel, C. Pineda, and J. Eisert, Phys. Rev. A 80, 042333 (2009)

  32. [40]

    W. M. C. Foulkes, L. Mitas, R. J. Needs, and G. Ra- jagopal, Rev. Mod. Phys. 73, 33 (2001)

  33. [41]

    A. W. Sandvik and G. Vidal, Phys. Rev. Lett. 99, 220602 (2007)

  34. [42]

    A. W. Sandvik, Phys. Rev. B 56, 11678 (1997)

  35. [43]

    See Supplemental Material for more detailed results, which includes Ref. [18]

  36. [44]

    J. L. CARDY, ed., 1 - Introduction to Theory of Finite- Size Scaling, vol. 2 of Current Physics–Sources and Com- ments (Elsevier, 1988)

  37. [45]

    A. P. Kampf, D. J. Scalapino, and S. R. White, Phys. Rev. B 64, 052509 (2001)

  38. [46]

    Chou and T.-K

    C.-P. Chou and T.-K. Lee, Phys. Rev. B 81, 060503(R) (2010)

Pith tools

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