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 →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
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.
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
- 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.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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)
- [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.
- [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.'
- [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)
- [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.
- [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.
- [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.
- [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.
- [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
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
free parameters (2)
- Finite-size scaling coefficients (a and b)
- Superconducting correlation decay exponents alpha_s and alpha_d =
alpha_s ~ 4.9, alpha_d ~ 4.4
assumptions (5)
- domain assumption The fPEPS variational ansatz with bond dimension D=12 spans the relevant ground-state manifold at 1/8 doping.
- domain assumption The thermodynamic-limit energy is obtained by a second-order polynomial in (L1L2)^-1/2.
- domain assumption Open boundary conditions do not qualitatively change the stripe orientation.
- domain assumption E0=-0.467775 at zero doping from Sandvik (Ref 41) is accurate.
- standard math Variational principle: a lower variational energy implies a state closer to the true ground state.
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
Reference graph
Works this paper leans on
-
[1]
J. G. Bednorz and K. A. M¨ uller, Z. Phys. B: Condens. Matter 64, 189 (1986)
1986
-
[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)
1987
-
[3]
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]
Damascelli, Z
A. Damascelli, Z. Hussain, and Z.-X. Shen, Rev. Mod. Phys. 75, 473 (2003)
2003
- [5]
-
[6]
P. A. Lee, N. Nagaosa, and X.-G. Wen, Rev. Mod. Phys. 78, 17 (2006)
2006
-
[7]
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)
work page 2000
-
[8]
J. Hubbard, Proceedings of the Royal Society of Lon- don A: Mathematical, Physical and Engineering Sciences 276, 238 (1963), ISSN 0080-4630
work page 1963
Show all 46 references
-
[9]
P. W. ANDERSON, Science 235, 1196 (1987), ISSN 0036-8075
1987
-
[10]
F. C. Zhang and T. M. Rice, Phys. Rev. B 37, 3759 (1988)
1988
-
[11]
Han, Q.-H
J. Han, Q.-H. Wang, and D.-H. Lee, International Jour- nal of Modern Physics B 15, 1117 (2001)
2001
-
[12]
Sherman and M
A. Sherman and M. Schreiber, European Physical Jour- nal B 32, 203 (2003)
2003
-
[13]
Capello, M
M. Capello, M. Raczkowski, and D. Poilblanc, Phys. Rev. B 77, 224502 (2008)
2008
-
[14]
C. S. Hellberg and E. Manousakis, Phys. Rev. Lett. 83, 132 (1999)
1999
-
[15]
S. R. White and D. J. Scalapino, Phys. Rev. Lett. 80, 1272 (1998)
1998
-
[16]
J. M. Tranquada, J. D. Axe, N. Ichikawa, A. R. Mood- enbaugh, Y. Nakamura, and S. Uchida, Phys. Rev. Lett. 78, 338 (1997)
1997
-
[17]
W.-J. Hu, F. Becca, and S. Sorella, Phys. Rev. B 85, 081110(R) (2012)
2012
-
[18]
Corboz, S
P. Corboz, S. R. White, G. Vidal, and M. Troyer, Phys. Rev. B 84, 041108(R) (2011)
2011
-
[19]
Corboz, T
P. Corboz, T. M. Rice, and M. Troyer, Phys. Rev. Lett. 113, 046402 (2014)
2014
-
[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)
2003
-
[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)
2003
-
[22]
S. R. White and D. J. Scalapino, Phys. Rev. B 61, 6320 (2000)
2000
-
[23]
S. R. White and D. J. Scalapino, Phys. Rev. B 70, 220506(R) (2004)
2004
-
[24]
S. R. White and D. J. Scalapino, Phys. Rev. B 79, 220504(R) (2009)
2009
-
[25]
Schollw¨ ock, Annals of Physics 326, 96 (2011)
U. Schollw¨ ock, Annals of Physics 326, 96 (2011)
2011
-
[26]
Verstraete, V
F. Verstraete, V. Murg, and J. Cirac, Advances in Physics 57, 143 (2008)
2008
-
[27]
H. C. Jiang, Z. Y. Weng, and T. Xiang, Phys. Rev. Lett. 101, 090603 (2008)
2008
- [28]
-
[29]
Sfondrini, J
A. Sfondrini, J. Cerrillo, N. Schuch, and J. I. Cirac, Ph ys. Rev. B 81, 214426 (2010)
2010
-
[30]
Verstraete, M
F. Verstraete, M. M. Wolf, D. Perez-Garcia, and J. I. Cirac, Phys. Rev. Lett. 96, 220601 (2006)
2006
-
[31]
Z.-C. Gu, F. Verstraete, and X.-G. Wen, arXiv:1004.256 3 (2010)
2010
-
[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)
2013
-
[33]
Corboz, R
P. Corboz, R. Or´ us, B. Bauer, and G. Vidal, Phys. Rev. B 81, 165104 (2010)
2010
-
[34]
C. V. Kraus, N. Schuch, F. Verstraete, and J. I. Cirac, Phys. Rev. A 81, 052338 (2010)
2010
-
[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)
2017
-
[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)
2018
-
[37]
S.-J. Dong, C. Wang, Y. Han, G.-c. Guo, and L. He, Phys. Rev. B 99, 195153 (2019)
2019
-
[38]
Normand and A
B. Normand and A. P. Kampf, Phys. Rev. B 64, 024521 (2001)
2001
-
[39]
Barthel, C
T. Barthel, C. Pineda, and J. Eisert, Phys. Rev. A 80, 042333 (2009)
2009
-
[40]
W. M. C. Foulkes, L. Mitas, R. J. Needs, and G. Ra- jagopal, Rev. Mod. Phys. 73, 33 (2001)
2001
-
[41]
A. W. Sandvik and G. Vidal, Phys. Rev. Lett. 99, 220602 (2007)
2007
-
[42]
A. W. Sandvik, Phys. Rev. B 56, 11678 (1997)
1997
-
[43]
See Supplemental Material for more detailed results, which includes Ref. [18]
-
[44]
J. L. CARDY, ed., 1 - Introduction to Theory of Finite- Size Scaling, vol. 2 of Current Physics–Sources and Com- ments (Elsevier, 1988)
1988
-
[45]
A. P. Kampf, D. J. Scalapino, and S. R. White, Phys. Rev. B 64, 052509 (2001)
2001
-
[46]
Chou and T.-K
C.-P. Chou and T.-K. Lee, Phys. Rev. B 81, 060503(R) (2010)
2010
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.