Pith. sign in

REVIEW 3 major objections 4 minor 5 cited by

The paper establishes that in three-color strong-coupling QCD with Kogut–Susskind quarks, the first-order chiral and nuclear transitions terminate at the same critical quark mass—m_c^χ=2.0545(34) and m_c^n=2.075(23) at Nτ=8—and that the fir

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

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

T0 review · deepseek-v4-flash

2026-08-03 07:14 UTC pith:YB6TILOO

load-bearing objection A credible first TRG computation of (3+1)d SU(3) strong-coupling QCD at finite density; the endpoint masses are plausible, but the power-law extrapolation and the 1024^4 convergence check need tightening before the quoted numbers are treated as definitive. the 3 major comments →

arxiv 2601.20690 v2 pith:YB6TILOO submitted 2026-01-28 hep-lat

Tensor renormalization group study of cold and dense QCD in the strong coupling limit

Yuto Sugimoto , Shinichiro Akiyama , Yoshinobu Kuramashi This is my paper
classification hep-lat
keywords strong coupling QCDtensor renormalization groupchiral transitionnuclear transitioncritical endpointfinite densityKogut-Susskind quarksfirst-order phase transition
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

Building on a tensor-network treatment of dense two-color QCD, the paper extends the method to three-color strong-coupling QCD with Kogut–Susskind quarks and uses it to map the cold, dense phase diagram. At temporal size Nτ=8 it finds that both the chiral condensate and the quark number density jump discontinuously at the same chemical potential for quark masses between 1.90 and 2.025, and that the two jumps vanish at the same critical mass within errors: m_c^χ=2.0545(34) and m_c^n=2.075(23). It also reports that at m=2.07 a first-order transition survives on a 1024^4 lattice, which it interprets as the thermodynamic limit at zero temperature, consistent with mean-field analysis. Because the finite-density region is inaccessible to ordinary Monte Carlo, the result suggests tensor renormalization group methods can probe cold dense QCD directly.

Core claim

On Nτ=8, the free energy develops a kink in the chemical potential for each quark mass m∈[1.90,2.025], and the size of the kink—measured as the gap in the quark number density and the chiral condensate across the transition—shrinks as m grows. Fitting the six gap values to Δ=A(m_c−m)^p yields m_c^χ=2.0545(34) for the chiral condensate and m_c^n=2.075(23) for the density, in agreement with each other, and for every m the transition chemical potentials from the two observables coincide. On a 1024^4 lattice at m=2.07, the kink sharpens as volume increases and both observables jump near μ≈1.591, so the paper concludes the transition is first order in the thermodynamic limit at zero temperature.

What carries the argument

The calculation is built on a Grassmann tensor network: after integrating out the SU(3) link variables and quarks at each site, the partition function becomes a network of sparse site tensors with twelve Grassmann indices per leg. Contraction is performed with the anisotropic tensor renormalization group, using bond dimension D=55 (and D3d=120 for the finite-temperature Nτ=8 computation), accelerated on GPUs. The paper locates first-order transitions by the kink in the thermodynamic potential and extracts the discontinuity sizes from local polynomial fits bracketing the kink; the critical masses are obtained by fitting those discontinuities to a power law in m.

Load-bearing premise

The endpoint values rest on the assumption that the discontinuities vanish as a simple power law in quark mass across six fitted points in m∈[1.90,2.025], with the extrapolation assumed insensitive to the chosen fit form and range; the 1024^4 zero-temperature claim separately assumes bond dimension D=55 suffices, a convergence check performed only at Nτ=8.

What would settle it

Run the Nτ=8 calculation at m=2.04 and 2.05 with bond dimension D≈70 and compare the measured gaps to the power-law predictions: the fitted curves predict specific nonzero values there, so a zero gap or a clear deviation would invalidate the endpoint extrapolation. Separately, re-running the 1024^4 contraction at m=2.07 with larger D would show whether the kink near μ≈1.591 is a truncation artifact.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • At Nτ=8 the chiral and nuclear transitions are degenerate: both observables jump at the same chemical potential for every mass studied, and the endpoint masses agree within errors.
  • The critical quark mass m_c≈2.05–2.08 sits between the dual Monte Carlo estimate (≈1.7) and the mean-field prediction (≈2.4), giving a numerical benchmark for strong-coupling QCD.
  • The persistence of a first-order transition at m=2.07 on a 1024^4 lattice supports the mean-field expectation that the finite-mass transition is first order rather than a crossover at zero temperature.
  • The method sidesteps the complex-action/sign problem that blocks direct Monte Carlo simulation in this region, so similar tensor-network contractions can be applied to nearby cold-dense models.

Where Pith is reading between the lines

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

  • The near-1/2 fitted exponents (p_n=0.47(10), p_χ=0.514(20)) suggest the endpoint may be mean-field-like, but the paper does not claim universality; checking the same exponent at larger Nτ would test that.
  • A single new data point at m=2.04 with larger bond dimension would test the power-law extrapolation directly, since the fitted curves make specific predictions there.
  • The paper's own outlook notes there is no symmetry forcing the chiral and nuclear transitions to share a chemical potential; turning on the gauge coupling is therefore a sharp test of whether the degeneracy seen here is special to the strong-coupling limit.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 4 minor

Summary. The paper applies the Grassmann anisotropic tensor renormalization group to (3+1)-dimensional strong-coupling QCD with Kogut--Susskind quarks at finite chemical potential. At Nτ=8 it computes the thermodynamic potential, quark number density, and chiral condensate, and locates first-order transition points for quark masses m∈[1.90,2.025]. The discontinuities in ⟨n⟩ and ⟨χ̄χ⟩ are fit to power laws, Eqs. (12)-(13), yielding critical masses m_c^n(Nτ=8)=2.075(23) and m_c^χ(Nτ=8)=2.0545(34), which are claimed to be consistent. The paper also reports a first-order transition at m=2.07 on a 1024^4 lattice, interpreted as essentially the thermodynamic limit at zero temperature, and compares the endpoints with dual Monte Carlo and mean-field results.

Significance. If the results hold, this is a significant technical advance: it extends TRG methods to a (3+1)-dimensional non-Abelian gauge theory with dynamical fermions at finite density, a regime where standard Monte Carlo suffers from the complex-action problem. The paper is self-contained in its Grassmann tensor formulation, provides explicit convergence checks (δ_3df<10^-5, δf=O(10^-6)), and demonstrates consistency between two independent observables—the chiral condensate and the quark number density—at the same transition chemical potential. This consistency is not imposed by construction and is a genuine strength. The comparison with dual-MC and mean-field estimates is useful for positioning the result in the literature. The main weakness is that the central endpoint extraction relies on a single six-point power-law fit whose systematic stability is not established.

major comments (3)
  1. [Sec. III A, Eqs. (12)-(13), Table I] The endpoint masses are the central result, but they rest entirely on six gap values at m=1.900,...,2.025 fitted to Δ=A(m_c−m)^p with p free. Over this interval the gaps change by only about a factor of two (Δ⟨n⟩: 0.633→0.350; Δ⟨χχ⟩: 0.692→0.308), and the fitted endpoints lie only Δm≈0.03–0.05 above the heaviest fitted mass. The paper reports no goodness-of-fit, no comparison with a fixed mean-field exponent p=1/2, no correction-to-scaling term, and no scan over the fit range. Because the agreement of m_c^χ and m_c^n is the paper's main claim, a shared or opposing bias in the assumed power-law form could change the conclusion. Please add fit-form and range systematics and quote the resulting systematic error.
  2. [Sec. III A, Table I and gap extraction] The gap magnitudes Δ⟨n⟩ and Δ⟨χχ⟩ are obtained by fitting four points with a quadratic in each phase for the number density and seven points with a third-order polynomial for the chiral condensate. The quoted errors appear to reflect only the location uncertainty |μ+−μ−| and not the uncertainty from the fitting window, polynomial order, or finite-difference step. Since these gaps are the input to the power-law fits of Eqs. (12)-(13), the parameter errors on A, m_c, and p are likely underestimated. Please propagate a systematic from the local fitting procedure, or justify that the chosen windows are uniquely determined by the data.
  3. [Sec. III B, Figs. 8-9] The zero-temperature run at m=2.07 on a 1024^4 lattice uses bond dimension D=55, but the convergence checks in Sec. III A (δ_3df, δf) are performed only on a 32^3×8 lattice at Nτ=8. The claim that the first-order jump seen at 1024^4 is 'essentially in the thermodynamic limit at zero temperature' requires a check that D=55 is sufficient at this volume and temporal extent. At minimum, show D-scan data at larger Nτ (e.g., Nτ=16 or 32) demonstrating that the jump and its location are stable, or state explicitly the computational restrictions preventing such a check.
minor comments (4)
  1. [Sec. III B] The sentence 'calculating μ dependence of the chiral condensate and the quark number density using Eqs. (13) and (11)' is incorrect: Eq. (13) is the power-law fit for Δ⟨χχ⟩, not the finite-difference formula for the chiral condensate. Please cite the correct equation or number the chiral condensate derivative.
  2. [Sec. III A, Fig. 3 inset] The volume dependence is shown only for f(m=1.9, μ). Please state explicitly which volumes are included and whether Ns=32 is sufficient for all quark masses studied, since the table reports results at a single volume.
  3. [Sec. III A, Table I] The error definition is unclear: the text says the quoted error corresponds to |μ+−μ−|, but Table I lists errors of 0.0005 for μ_c^n and 0.001 for μ_c^χ. Clarify how these errors are obtained and why they differ.
  4. [Sec. III A, Figs. 5 and 7] The figures would be more informative with residuals or reduced χ² values for the fits. This would also help the reader judge the adequacy of the power-law ansatz.

Circularity Check

0 steps flagged

No significant circularity; the agreement of m_c^chi and m_c^n is an output of separate power-law extrapolations, not a fitted input.

full rationale

The derivation is self-contained in the sense required by the circularity standard. The chiral condensate and quark number density are both numerical derivatives of the same ln Z (Eq. (11) and the analogous m-derivative), so they can share the location of a first-order kink; nevertheless, the two critical masses are extracted from separately measured discontinuities Delta< n >(m) and Delta< chichi >(m), fitted to independent power-law forms Eqs. (12)-(13) with free endpoints and free exponents. Nothing in these fits enforces m_c^n = m_c^chi; their consistency is an output, not an input. The power-law ansatz with free exponent is a fit-form/statistical-modeling choice and is a legitimate extrapolation-reliability concern (six points, endpoints only ~0.03-0.05 above the heaviest fitted mass), but it is not an instance of fitting a parameter to a subset and then 'predicting' that subset, nor of defining X in terms of Y. The comparisons to dual-formulation Monte Carlo and mean-field results are external benchmarks rather than self-citations. Self-citations that occur [26,28,31,58-60] are to methodological prior work (ATRG/GTRG algorithms and the QC2D tensor construction) and are not used as unverified authority to force the phase structure; the paper presents its own tensor formulation and numerical checks. Hence no load-bearing circular step is found.

Axiom & Free-Parameter Ledger

8 free parameters · 7 axioms · 0 invented entities

The central numerical results rest on fit parameters A, m_c, p that are fitted to data rather than predicted, and on TRG truncation and strong-coupling model assumptions. No new particles, forces, or theoretical entities are introduced.

free parameters (8)
  • A_n (amplitude of power-law fit to quark number density jump) = 1.43(16)
    Fitted to six discontinuity values in Eq. (12); affects m_c^n estimate.
  • m_c^n (critical quark mass from number density) = 2.075(23)
    Central result; not directly observed, obtained by extrapolation of fitted power law.
  • p_n (exponent of power-law fit) = 0.47(10)
    Free exponent in Eq. (12); no theory value imposed.
  • A_χ (amplitude of power-law fit to chiral condensate jump) = 1.81(5)
    Fitted to discontinuity values in Eq. (13).
  • m_c^χ (critical quark mass from chiral condensate) = 2.0545(34)
    Central result; extracted from power-law extrapolation.
  • p_χ (exponent of power-law fit) = 0.514(20)
    Free exponent in Eq. (13).
  • TRG bond dimensions (D, D3d) = D=55, D3d=120
    Chosen by hand; convergence checked at Nτ=8, not at 1024^4.
  • Finite-difference steps (Δμ, Δm) = 0.001
    Used for numerical derivatives; near first-order transitions finite steps can mask discontinuities, as the authors note.
axioms (7)
  • domain assumption Strong-coupling limit: gauge action absent (1/g^2=0)
    The entire calculation is for strong-coupling lattice QCD, not full QCD; the authors state that going to finite couplings is future work (Sec. IV).
  • domain assumption Kogut-Susskind staggered action with one flavor and N_c=3
    Used throughout Eq. (1); results may depend on this particular lattice fermion action.
  • standard math Grassmann Gaussian integration and fermionic sign counting in Appendix A
    The explicit coefficient tensor assumes correctness of Grassmann algebra identities and the sign factor R; there is no independent check in the paper.
  • domain assumption TRG truncation with finite bond dimensions approximates the exact tensor contraction
    Convergence is checked for Nτ=8 in Eqs. (9)-(10) and Figs. 1-2, but is assumed for the 1024^4 computation in Sec. III B without an explicit check.
  • ad hoc to paper Power-law scaling of discontinuities near the endpoint (Eqs. (12)-(13))
    No RG derivation of the functional form or of the free exponent p is given; the fit assumes all data lie in the scaling regime.
  • domain assumption Derivative approximations with Δμ=Δm=0.001 faithfully represent thermodynamic quantities away from the transition
    The authors note finite Δμ can mask discontinuities (Sec. III A) and use a kink-fitting procedure instead, but the same step sizes underlie the reported values.
  • domain assumption Equivalence of 1024^4 lattice to the thermodynamic limit at zero temperature
    The paper argues the kink develops rapidly with volume, but no systematic finite-size scaling is shown.

pith-pipeline@v1.3.0-alltime-deepseek · 15945 in / 14718 out tokens · 150928 ms · 2026-08-03T07:14:44.360886+00:00 · methodology

0 comments
read the original abstract

We study the phase structure of the (3+1)-dimensional cold and dense QCD with the Kogut--Susskind quark in the strong coupling limit using the tensor renormalization group method. The chiral and nuclear transitions are investigated by calculating the chiral condensate and the quark number density as a function of the chemical potential. For a fixed temporal extent $N_\tau=8$, we determine the critical quark masses $m_c^{\chi}$ and $m_c^{n}$ for the chiral condensate and the quark number density, respectively, at which the first-order phase transition terminates with the vanishing discontinuity in thermodynamic quantities. We find that both quantities at the same quark mass exhibit a discontinuity at the same chemical potential, and the resulting critical quark masses are consistent with each other. We also compare our results for the critical quark masses with those obtained from the Monte Carlo simulation in the dual formulation and from the mean-field analysis. We further confirm the first-order phase transition at finite quark mass on a $1024^4$ lattice, which is essentially in the thermodynamic limit at zero temperature, as expected from the mean-field analysis.

Figures

Figures reproduced from arXiv: 2601.20690 by Shinichiro Akiyama, Yoshinobu Kuramashi, Yuto Sugimoto.

Figure 1
Figure 1. Figure 1: FIG. 1. Relative error of thermodynamic potential in terms of [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Relative error of thermodynamic potential in terms of [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4 [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Fit of ∆ [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6 [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Fit of ∆ [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8 [PITH_FULL_IMAGE:figures/full_fig_p010_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9 [PITH_FULL_IMAGE:figures/full_fig_p011_9.png] view at source ↗

discussion (0)

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

Forward citations

Cited by 5 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Deconfinement from Thermal Tensor Networks: Universal CFT signature in (2+1)-dimensional $\mathbb{Z}_N$ lattice gauge theory

    hep-th 2026-02 conditional novelty 7.0

    Tensor-network contraction of finite-temperature Z_N gauge theory yields central charges and scaling dimensions consistent with Svetitsky–Yaffe universality for N=2,3,5, including a U(1)-symmetric BKT phase for N=5, a...

  2. Phase diagram of the single-flavor Gross--Neveu--Wilson model from the Grassmann corner transfer matrix renormalization group

    hep-lat 2026-02 conditional novelty 6.0

    For one-flavor Gross-Neveu-Wilson fermions, the Aoki phase is bounded by c=1/2 Ising critical lines and terminates at strong coupling, while c=1 lines separate topological and trivial insulators.

  3. Multi-particle states investigation with tensor renormalization group method

    hep-lat 2026-06 unverdicted novelty 5.0

    A TRG-based spectroscopy scheme identifies multi-particle states in the 1+1d Ising model and extracts consistent two-particle scattering phase shifts via Lüscher's formula.

  4. Lattice field theories with a sign problem

    hep-lat 2026-04 unverdicted novelty 2.0

    A review of holomorphic extensions, dual variables, tensor renormalization group, and machine learning approaches for controlling the sign problem in lattice field theories.

  5. Lattice field theories with a sign problem

    hep-lat 2026-04 unverdicted novelty 1.0

    Reviews approaches such as Lefschetz thimbles, complex Langevin dynamics, dual variables, tensor renormalization group, and machine learning to control the sign problem in lattice field theories.

Reference graph

Works this paper leans on

64 extracted references · 51 linked inside Pith · cited by 4 Pith papers

  1. [1]

    B. C. Barrois,Superconducting Quark Matter,Nucl. Phys. B129(1977) 390–396

    1977

  2. [2]

    Asakawa and K

    M. Asakawa and K. Yazaki,Chiral Restoration at Finite Density and Temperature,Nucl. Phys. A504(1989) 668–684

    1989

  3. [3]

    M. G. Alford, K. Rajagopal and F. Wilczek,QCD at finite baryon density: Nucleon droplets and color superconductivity, Phys. Lett. B422(1998) 247–256, [hep-ph/9711395]

    Pith/arXiv arXiv 1998

  4. [4]

    R. Rapp, T. Sch¨ afer, E. V. Shuryak and M. Velkovsky,Diquark Bose condensates in high density matter and instantons, Phys. Rev. Lett.81(1998) 53–56, [hep-ph/9711396]

    Pith/arXiv arXiv 1998

  5. [5]

    M. G. Alford, K. Rajagopal and F. Wilczek,Color flavor locking and chiral symmetry breaking in high density QCD, Nucl. Phys. B537(1999) 443–458, [hep-ph/9804403]

    Pith/arXiv arXiv 1999

  6. [6]

    Karsch and K

    F. Karsch and K. H. Mutter,STRONG COUPLING QCD AT FINITE BARYON NUMBER DENSITY,Nucl. Phys. B 313(1989) 541–559

    1989

  7. [7]

    Delgado Mercado, C

    Y. Delgado Mercado, C. Gattringer and A. Schmidt,Dual Lattice Simulation of the Abelian Gauge-Higgs Model at Finite Density: An Exploratory Proof of Concept Study,Phys. Rev. Lett.111(2013) 141601, [1307.6120]

    Pith/arXiv arXiv 2013

  8. [8]

    Wolff,Baryons in Lattice QCD at Strong Coupling,Phys

    U. Wolff,Baryons in Lattice QCD at Strong Coupling,Phys. Lett. B153(1985) 92–96

    1985

  9. [9]

    Fromm and P

    M. Fromm and P. de Forcrand,Revisiting strong coupling QCD at finite temperature and baryon density,PoS LA TTICE2008(2008) 191, [0811.1931]

    Pith/arXiv arXiv 2008

  10. [10]

    Miura, T

    K. Miura, T. Z. Nakano, A. Ohnishi and N. Kawamoto,Phase diagram evolution at finite coupling in strong coupling lattice QCD,Phys. Rev. D80(2009) 074034, [0907.4245]

    Pith/arXiv arXiv 2009

  11. [11]

    de Forcrand and M

    P. de Forcrand and M. Fromm,Nuclear Physics from lattice QCD at strong coupling,Phys. Rev. Lett.104(2010) 112005, [0907.1915]

    Pith/arXiv arXiv 2010

  12. [12]

    de Forcrand, J

    P. de Forcrand, J. Langelage, O. Philipsen and W. Unger,Lattice QCD Phase Diagram In and Away from the Strong Coupling Limit,Phys. Rev. Lett.113(2014) 152002, [1406.4397]

    Pith/arXiv arXiv 2014

  13. [13]

    Kawamoto and J

    N. Kawamoto and J. Smit,Effective Lagrangian and Dynamical Symmetry Breaking in Strongly Coupled Lattice QCD, Nucl. Phys. B192(1981) 100

    1981

  14. [14]

    Bilic, F

    N. Bilic, F. Karsch and K. Redlich,Flavor dependence of the chiral phase transition in strong coupling QCD,Phys. Rev. D45(1992) 3228–3236

    1992

  15. [15]

    Nishida,Phase structures of strong coupling lattice QCD with finite baryon and isospin density,Phys

    Y. Nishida,Phase structures of strong coupling lattice QCD with finite baryon and isospin density,Phys. Rev. D69 (2004) 094501, [hep-ph/0312371]

    Pith/arXiv arXiv 2004

  16. [16]

    T. Z. Nakano, K. Miura and A. Ohnishi,Effective Potential in the Strong-coupling Lattice QCD with Next-to-Next-to-Leading Order Effects,Prog. Theor. Phys.123(2010) 825–851, [0911.3453]

    Pith/arXiv arXiv 2010

  17. [17]

    J. Kim, P. Pattanaik and W. Unger,Nuclear liquid-gas transition in the strong coupling regime of lattice QCD,Phys. Rev. D107(2023) 094514, [2303.01467]

    Pith/arXiv arXiv 2023

  18. [18]

    Tsutsui, Y

    S. Tsutsui, Y. Asano, Y. Ito, H. Matsufuru, Y. Namekawa, J. Nishimura et al.,Color superconductivity in a small box: a complex Langevin study,PoSLA TTICE2021(2022) 533, [2111.15095]

    Pith/arXiv arXiv 2022

  19. [19]

    Levin and C

    M. Levin and C. P. Nave,Tensor renormalization group approach to two-dimensional classical lattice models,Phys. Rev. Lett.99(2007) 120601, [cond-mat/0611687]

    Pith/arXiv arXiv 2007

  20. [20]

    Z. Y. Xie, J. Chen, M. P. Qin, J. W. Zhu, L. P. Yang and T. Xiang,Coarse-graining renormalization by higher-order singular value decomposition,Phys. Rev. B86(Jul, 2012) 045139, [1201.1144]

    Pith/arXiv arXiv 2012

  21. [21]

    Shimizu and Y

    Y. Shimizu and Y. Kuramashi,Grassmann tensor renormalization group approach to one-flavor lattice Schwinger model, Phys. Rev.D90(2014) 014508, [1403.0642]

    Pith/arXiv arXiv 2014

  22. [22]

    Evenbly and G

    G. Evenbly and G. Vidal,Tensor network renormalization,Phys. Rev. Lett.115(Oct, 2015) 180405

    2015

  23. [23]

    Sakai, S

    R. Sakai, S. Takeda and Y. Yoshimura,Higher order tensor renormalization group for relativistic fermion systems,PTEP 2017(2017) 063B07, [1705.07764]. 3 One can impose stricter constraints. However, the present framework is sufficient for this model. 15

    Pith/arXiv arXiv 2017

  24. [24]

    Yang, Z.-C

    S. Yang, Z.-C. Gu and X.-G. Wen,Loop optimization for tensor network renormalization,Phys. Rev. Lett.118(Mar,

  25. [25]

    Hauru, C

    M. Hauru, C. Delcamp and S. Mizera,Renormalization of tensor networks using graph independent local truncations, Phys. Rev.B97(2018) 045111, [1709.07460]

    Pith/arXiv arXiv 2018

  26. [26]

    Adachi, T

    D. Adachi, T. Okubo and S. Todo,Anisotropic Tensor Renormalization Group,Phys. Rev. B102(2020) 054432, [1906.02007]

    Pith/arXiv arXiv 2020

  27. [27]

    Kadoh and K

    D. Kadoh and K. Nakayama,Renormalization group on a triad network,1912.02414

    Pith/arXiv arXiv 1912

  28. [28]

    Akiyama, Y

    S. Akiyama, Y. Kuramashi, T. Yamashita and Y. Yoshimura,Restoration of chiral symmetry in cold and dense Nambu–Jona-Lasinio model with tensor renormalization group,JHEP01(2021) 121, [2009.11583]

    Pith/arXiv arXiv 2021

  29. [29]

    Adachi, T

    D. Adachi, T. Okubo and S. Todo,Bond-weighted tensor renormalization group,Phys. Rev. B105(Feb, 2022) L060402, [2011.01679]

    Pith/arXiv arXiv 2022

  30. [30]

    Akiyama,Bond-weighting method for the Grassmann tensor renormalization group,JHEP11(2022) 030, [2208.03227]

    S. Akiyama,Bond-weighting method for the Grassmann tensor renormalization group,JHEP11(2022) 030, [2208.03227]

    Pith/arXiv arXiv 2022

  31. [31]

    Sugimoto, S

    Y. Sugimoto, S. Akiyama and Y. Kuramashi,Phase structure of (3+1)-dimensional dense two-color QCD atT= 0in the strong coupling limit with tensor renormalization group,2509.23637

    arXiv

  32. [32]

    Bloch and R

    J. Bloch and R. Lohmayer,Grassmann higher-order tensor renormalization group approach for two-dimensional strong-coupling QCD,Nucl. Phys. B986(2023) 116032, [2206.00545]

    Pith/arXiv arXiv 2023

  33. [33]

    Asaduzzaman, S

    M. Asaduzzaman, S. Catterall, Y. Meurice, R. Sakai and G. C. Toga,Tensor network representation of non-abelian gauge theory coupled to reduced staggered fermions,JHEP05(2024) 195, [2312.16167]

    Pith/arXiv arXiv 2024

  34. [34]

    K. H. Pai, S. Akiyama and S. Todo,Grassmann tensor renormalization group approach to (1+1)-dimensional two-color lattice QCD at finite density,JHEP03(2025) 027, [2410.09485]

    Pith/arXiv arXiv 2025

  35. [35]

    K. H. Pai, S. Akiyama and S. Todo,Two-color lattice QCD in (1 + 1) dimensions with Grassmann tensor renormalization group,PoSLA TTICE2024(2025) 364, [2501.18918]

    Pith/arXiv arXiv 2025

  36. [36]

    Shimizu and Y

    Y. Shimizu and Y. Kuramashi,Critical behavior of the lattice Schwinger model with a topological term atθ=πusing the Grassmann tensor renormalization group,Phys. Rev.D90(2014) 074503, [1408.0897]

    Pith/arXiv arXiv 2014

  37. [37]

    Kawauchi and S

    H. Kawauchi and S. Takeda,Tensor renormalization group analysis of CP(N-1) model,Phys. Rev.D93(2016) 114503, [1603.09455]

    Pith/arXiv arXiv 2016

  38. [38]

    Kawauchi and S

    H. Kawauchi and S. Takeda,Loop-TNR analysis of CP(1) model with theta term,EPJ Web Conf.175(2018) 11015, [1710.09804]

    Pith/arXiv arXiv 2018

  39. [39]

    L.-P. Yang, Y. Liu, H. Zou, Z. Xie and Y. Meurice,Fine structure of the entanglement entropy in the O(2) model,Phys. Rev. E93(2016) 012138, [1507.01471]

    Pith/arXiv arXiv 2016

  40. [40]

    Shimizu and Y

    Y. Shimizu and Y. Kuramashi,Berezinskii-Kosterlitz-Thouless transition in lattice Schwinger model with one flavor of Wilson fermion,Phys. Rev.D97(2018) 034502, [1712.07808]

    Pith/arXiv arXiv 2018

  41. [41]

    Takeda and Y

    S. Takeda and Y. Yoshimura,Grassmann tensor renormalization group for the one-flavor lattice Gross-Neveu model with finite chemical potential,PTEP2015(2015) 043B01, [1412.7855]

    Pith/arXiv arXiv 2015

  42. [42]

    Kadoh, Y

    D. Kadoh, Y. Kuramashi, Y. Nakamura, R. Sakai, S. Takeda and Y. Yoshimura,Tensor network formulation for two-dimensional latticeN= 1 Wess-Zumino model,JHEP03(2018) 141, [1801.04183]

    Pith/arXiv arXiv 2018

  43. [43]

    Kadoh, Y

    D. Kadoh, Y. Kuramashi, Y. Nakamura, R. Sakai, S. Takeda and Y. Yoshimura,Investigation of complexϕ 4 theory at finite density in two dimensions using TRG,JHEP02(2020) 161, [1912.13092]

    Pith/arXiv arXiv 2020

  44. [44]

    Takeda,Tensor network approach to real-time path integral,PoSLA TTICE2019(2019) 033, [1908.00126]

    S. Takeda,Tensor network approach to real-time path integral,PoSLA TTICE2019(2019) 033, [1908.00126]

    Pith/arXiv arXiv 2019

  45. [45]

    N. Butt, S. Catterall, Y. Meurice, R. Sakai and J. Unmuth-Yockey,Tensor network formulation of the massless Schwinger model with staggered fermions,Phys. Rev. D101(2020) 094509, [1911.01285]

    Pith/arXiv arXiv 2020

  46. [46]

    Kuramashi and Y

    Y. Kuramashi and Y. Yoshimura,Tensor renormalization group study of two-dimensional U(1) lattice gauge theory with aθterm,JHEP04(2020) 089, [1911.06480]

    Pith/arXiv arXiv 2020

  47. [47]

    Akiyama, D

    S. Akiyama, D. Kadoh, Y. Kuramashi, T. Yamashita and Y. Yoshimura,Tensor renormalization group approach to four-dimensional complexϕ 4 theory at finite density,JHEP09(2020) 177, [2005.04645]

    Pith/arXiv arXiv 2020

  48. [48]

    Nakayama, L

    K. Nakayama, L. Funcke, K. Jansen, Y.-J. Kao and S. K¨ uhn,Phase structure of the CP(1) model in the presence of a topologicalθ-term,Phys. Rev. D105(2022) 054507, [2107.14220]

    Pith/arXiv arXiv 2022

  49. [49]

    Bloch, R

    J. Bloch, R. G. Jha, R. Lohmayer and M. Meister,Tensor renormalization group study of the three-dimensional O(2) model,Phys. Rev. D104(2021) 094517, [2105.08066]

    Pith/arXiv arXiv 2021

  50. [50]

    Akiyama and Y

    S. Akiyama and Y. Kuramashi,Critical endpoint of (3+1)-dimensional finite densityZ 3 gauge-Higgs model with tensor renormalization group,JHEP10(2023) 077, [2304.07934]

    Pith/arXiv arXiv 2023

  51. [51]

    Akiyama and Y

    S. Akiyama and Y. Kuramashi,Tensor renormalization group study of (1 + 1)-dimensional U(1) gauge-Higgs model atθ =πwith L¨ uscher’s admissibility condition,JHEP09(2024) 086, [2407.10409]

    Pith/arXiv arXiv 2024

  52. [52]

    Luo and Y

    X. Luo and Y. Kuramashi,Tensor renormalization group approach to (1+1)-dimensional SU(2) principal chiral model at finite density,Phys. Rev. D107(2023) 094509, [2208.13991]

    Pith/arXiv arXiv 2023

  53. [53]

    Hite and Y

    M. Hite and Y. Meurice,Quantum real-time evolution using tensor renormalization group methods,Phys. Rev. D111 (2025) 034502, [2411.05301]

    Pith/arXiv arXiv 2025

  54. [54]

    Kanno, S

    H. Kanno, S. Akiyama, K. Murakami and S. Takeda,Grassmann tensor renormalization group for the massive Schwinger 16 model with aθterm using staggered fermions,JHEP11(2025) 036, [2412.08959]

    arXiv 2025

  55. [55]

    Luo and Y

    X. Luo and Y. Kuramashi,Quantum phase transition of (1+1)-dimensional O(3) nonlinear sigma model at finite density with tensor renormalization group,JHEP11(2024) 144, [2406.08865]

    Pith/arXiv arXiv 2024

  56. [56]

    Luo and Y

    X. Luo and Y. Kuramashi,Critical endpoints of three-dimensional finite density SU(3) spin model with tensor renormalization group,JHEP07(2025) 036, [2503.05144]

    Pith/arXiv arXiv 2025

  57. [57]

    Aizawa, S

    H. Aizawa, S. Takeda and Y. Yoshimura,Phase structure analysis of 2d lattice CP(1) model withθterm using tensor renormalization group method,2510.06624

    arXiv

  58. [58]

    Akiyama and D

    S. Akiyama and D. Kadoh,More about the Grassmann tensor renormalization group,JHEP10(2021) 188, [2005.07570]

    Pith/arXiv arXiv 2021

  59. [59]

    Oba,Cost Reduction of Swapping Bonds Part in Anisotropic Tensor Renormalization Group,PTEP2020(2020) 013B02, [1908.07295]

    H. Oba,Cost Reduction of Swapping Bonds Part in Anisotropic Tensor Renormalization Group,PTEP2020(2020) 013B02, [1908.07295]

    Pith/arXiv arXiv 2020

  60. [60]

    Sugimoto and S

    Y. Sugimoto and S. Sasaki,Triad representation for the anisotropic tensor renormalization group in four dimensions, Phys. Rev. D112(2025) 094514, [2507.21909]

    arXiv 2025

  61. [61]

    Akiyama and Y

    S. Akiyama and Y. Kuramashi,Tensor renormalization group approach to (1+1)-dimensional Hubbard model,Phys. Rev. D104(2021) 014504, [2105.00372]

    Pith/arXiv arXiv 2021

  62. [62]

    Akiyama, Y

    S. Akiyama, Y. Kuramashi and T. Yamashita,Metal-insulator transition in (2+1)-dimensional Hubbard model with tensor renormalization group,PTEP2022(9, 2021) 023, [2109.14149]

    Pith/arXiv arXiv 2021

  63. [63]

    A. Ueda, S. De Meyer, A. Naravane, V. Vanthilt and F. Verstraete,Global Tensor Network Renormalization for 2D Quantum systems: A new window to probe universal data from thermal transitions,2508.05406

    Pith/arXiv arXiv

  64. [64]

    R. B. Lehoucq, D. C. Sorensen and C. Yang,ARPACK Users’ Guide. Society for Industrial and Applied Mathematics, 1998, 10.1137/1.9780898719628. 17

    doi:10.1137/1.9780898719628 1998