{"id":"f3c01e4e-6f37-4eb8-8f98-445da7eb9f86","arxiv_id":"2602.21705","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"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.","lead":"This paper maps the full phase diagram of a two-dimensional toy fermion theory using a sign-problem-free tensor-network method, finding where a symmetry-breaking 'Aoki' phase exists and how it ends at strong coupling. It also identifies the universal behavior at each phase boundary and distinguishes a topological insulator from a trivial one.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Admitted anomaly in finite-entanglement scaling exponent leaves the c=1/2 and c=1 assignments in need of independent confirmation.","rationale":"The reader's weakest assumption correctly identifies the finite-entanglement scaling relation as the most load-bearing point. The paper's own Appendix C admits a large unexplained deviation in κ, so this is an explicit, documented soft spot rather than a manufactured one. All central-charge assignments in Fig. 5 use Eq. (III.16), so if the anomaly reflects a failure of that scaling law, the phase diagram's universal classes collapse. The concern does not warrant rejection: the direct c fits are internally consistent, the strong-coupling termination of the Aoki phase is independently corroborated by HOTRG in Appendix B, and the overall phase structure is consistent with the Hamiltonian-formalism results of Ref. [23] where applicable. But because the quantitative content of the phase diagram is the c values, and those values lack an independent determination, the CONDITIONAL verdict is appropriate. No change to the reader's verdict is needed.","tokens_in":20431,"tokens_out":15900,"duration_ms":155326,"concrete_test":"Recompute the central charge at (M,g^2)=(0.1,0.77181) and (0.1,0.92161) without relying on Eq. (III.16): fit the pseudoscalar condensate data of Fig. 8 to π ∝ |g^2 − g_c^2|^β with g_c fixed; the 2D Ising universality required for c=1/2 predicts β=1/8. If β is inconsistent with 1/8 (beyond the quoted errors), the c=1/2 assignments are not supported. As a complementary direct test, run CTMRG at D=192 and 256 at the same couplings and check whether the S_D versus log ξ_D slope stays fixed at c/6 and whether κ approaches the predicted 2.03; a persistent κ≈1.48 or a drifting slope would confirm uncontrolled corrections to the finite-entanglement scaling ansatz.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim—that the Aoki phase is bounded by c=1/2 Ising lines and that the topological/trivial boundary has c=1—rests entirely on Eq. (III.16), S_D ≃ (c/6) log ξ_D, applied to the CTMRG environments at each purported critical point. Appendix C reports that the associated finite-entanglement scaling ξ_D ∼ D^κ gives κ ≈ 1.48 at the c=1/2 points, far from the prediction κ ≈ 2.03 in Eq. (C.2) for c=1/2, and the authors state this deviation is 'significantly larger' than earlier reports and leave its origin to future work. If this discrepancy signals that the finite-entanglement scaling ansatz has uncontrolled corrections in this Grassmann CTMRG implementation, then the slopes in Figs. 10 and 12 need not equal c/6, and the extracted values c=0.498(3), c=0.500(4), c=1.01(3) would not establish the claimed universal classes. No independent check of c—for example from critical exponents, finite-size scaling, or a different tensor-network algorithm—is provided. The numerical phase diagram itself is plausible and partially cross-checked (HOTRG confirms the strong-coupling termination), but the universal-class labels are load-bearing and rely on a scaling law whose validity the paper itself calls into question.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a Grassmann corner transfer matrix renormalization group (CTMRG) algorithm for two-dimensional lattice fermions and applies it to the single-flavor Gross–Neveu–Wilson (GNW) model. The partition function is represented as a Grassmann tensor network and contracted by CTMRG, with the pseudoscalar condensate used to identify the Aoki phase, the entanglement entropy used to extract central charges of phase boundaries, and the entanglement spectrum used to identify a topological insulating phase. The authors report a phase diagram with an Aoki phase bounded by c=1/2 critical lines, topological and trivial phases separated by c=1 critical lines, and no Aoki phase at strong coupling, with a triple point around (M,g^2)=(0.812,0.89). The numerical machinery is benchmarked against free Wilson fermions and cross-checked partially with HOTRG.","tokens_in":20678,"tokens_out":2468,"duration_ms":24279,"significance":"If the central-charge assignments are reliable, this is an important first complete Lagrangian-formulation phase diagram of the N_f=1 GNW model, with nontrivial implications for lattice QCD with odd flavor numbers and for tensor-network methods for fermions. The paper ships a new Grassmann CTMRG implementation and demonstrates its superior accuracy relative to TRG/BTRG/HOTRG in benchmark tests (Fig. 4), plus an independent HOTRG cross-check of the strong-coupling termination of the Aoki phase (Appendix B). The entanglement-spectrum doubling inside the lobe is a clear, falsifiable signature. However, the universal-class labels in Fig. 5 rest on a finite-entanglement scaling relation whose validity the paper itself calls into question, and the pseudoscalar-condensate definition involves an h->0 limit whose numerical implementation is not documented.","major_comments":[{"comment":"The central-charge assignments c=1/2 and c=1, which are load-bearing for the phase diagram in Fig. 5, rely entirely on Eq. (III.16), S_D approx (c/6) log xi_D, applied at each claimed critical point. Appendix C reports that the finite-entanglement scaling xi_D ~ D^kappa requires kappa approx 1.48 for the best data collapse at the c=1/2 points, while Eq. (C.2) predicts kappa approx 2.03 for c=1/2. The authors state that this deviation is 'significantly larger' than previously reported and leave its origin to future work. If Eq. (III.16) has uncontrolled corrections in this Grassmann CTMRG implementation, the slopes in Figs. 10 and 12 do not necessarily equal c/6, and the extracted values c=0.498(3), c=0.500(4), c=1.01(3) do not establish the claimed universality classes. An independent confirmation, e.g. from critical exponents, finite-size scaling, or a different tensor-network algorithm","section":"Appendix C and Sec. IV.B.2"},{"comment":"The pseudoscalar condensate is defined by Eq. (III.3) with a double limit: first the thermodynamic limit, then h->0. The impurity tensor I_n is introduced, but the numerical sections (Figs. 6-8 and 15) never state the values of h used, whether results are extrapolated in h, or how the h->0 limit is implemented in the CTMRG contraction. This is not a mere presentation issue because the magnitude of the condensate, including its vanishing at strong coupling, is a central claim. Without a specified h-extrapolation procedure, the reader cannot assess systematic errors in the order parameter or in the location of the Aoki-phase boundaries extracted from it.","section":"Sec. III.A, Eqs. (III.3)-(III.7)"},{"comment":"The triple-point estimate (M,g^2) approx (0.812,0.89) is based on the difference Delta M between correlation-length peaks as a function of g^2, but no extrapolation to D->infinity or a criterion for 'vanishingly small' Delta M is given. The text says the two transition points are close and a reliable finite-entanglement analysis is left for future work, yet the triple point is used in the schematic phase diagram (Fig. 14) and in the summary. An uncertainty estimate for the triple-point location would be needed to support this part of the phase diagram.","section":"Sec. IV.B.4, Fig. 17"}],"minor_comments":[{"comment":"The heat map legend is not labeled; it is unclear whether it represents the absolute value of the pseudoscalar condensate on a linear or logarithmic scale. Adding a color-bar label and a scale would improve readability.","section":"Fig. 5"},{"comment":"The reduced density matrix rho_D is defined graphically; the text would benefit from a brief verbal description of how the four corner matrices are contracted and traced to yield a D x D density matrix.","section":"Sec. III.B.1, Eq. (III.15)"},{"comment":"The caption calls Fig. 13(c) 'the SPT phase', while the main text says it is inside the lobe and later identifies it as a topological insulator. This terminology is inconsistent; the figure label should match the phase name used in the text.","section":"Sec. IV.B.3, Fig. 13"},{"comment":"The formula kappa = 6/[c(sqrt(12/c)+1)] is quoted from MPS literature. It would be useful to clarify whether this expression is expected to hold exactly for CTMRG of a two-dimensional classical system or only approximately, given that Appendix C itself finds a substantially different kappa.","section":"Appendix C, Eq. (C.2)"},{"comment":"Several spots have missing spaces or typographical issues, e.g. 'N f ' in the introduction and 'g2 = 0.9' in Sec. IV.B.4. A careful proofread is recommended.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is technically substantial and the numerical method appears sound in its benchmarks. The main obstacle is that the central universal-class labels depend on a finite-entanglement scaling relation whose own internal consistency check (Appendix C) shows a large unexplained deviation in kappa. This is not a minor caveat but a load-bearing issue. I would support publication if the authors either provide an independent determination of c (e.g. from the critical exponent nu via data collapse, or from a different algorithm) or substantially reduce the kappa discrepancy. The h->0 implementation gap should also be fixed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis paper reports the first path-integral (Lagrangian) phase diagram of the single-flavor Gross-Neveu-Wilson model, computed with Grassmann CTMRG. The qualitative picture—an Aoki phase bounded by c=1/2 lines, a two-lobe topological-insulator region separated from the trivial phase by c=1 lines, and a strong-coupling termination of the Aoki phase—is new, and the method is benchmarked against free Wilson fermions and cross-checked with an independent HOTRG calculation. That cross-check gives the qualitative phase diagram some weight.\n\nThe genuinely new results are the complete Lagrangian phase diagram and the strong-coupling termination. The two-lobe structure was anticipated by large-Nf and by the Hamiltonian MPS work, so the novelty is in the path-integral confirmation, not in the concept. The CTMRG implementation looks like a sound extension of existing Grassmann tensor methods, and the quoted central charges come with statistical uncertainties.\n\nThe main soft spot is the central charge extraction. It rests on the finite-entanglement scaling relation S_D = (c/6) log xi_D, and the paper itself reports in Appendix C that the best data collapse at the c=1/2 points requires kappa ~ 1.48, whereas Eq. (C.2) for c=1/2 gives kappa ~ 2.03. The authors say the deviation is 'significantly larger' than earlier reports and leave it to future work. If the scaling ansatz is not under control in this implementation, the slopes in Figs. 10 and 12 need not equal c/6, and the extracted values c=0.498(3), c=0.500(4), c=1.01(3) do not establish the claimed universal classes. No independent check—critical exponents, finite-size scaling, another algorithm—is offered. This is not a minor technicality; it is the load-bearing assumption for the universality labels.\n\nSmaller issues: the h->0 limit in Eq. (III.3) is not described (how small h is, how the extrapolation is done); the SPT phase is identified only through doubly degenerate entanglement spectra, not through a topological invariant like the Zak phase, which is mentioned but not computed; and the tension with the 't Hooft anomaly argument in Ref. [15] is acknowledged but not resolved. No code or data are shipped, so independent reproduction is not possible from the text alone.\n\nNone of this is fatal. The qualitative phase diagram is credible, and the HOTRG check argues against a methodological artifact. But the universal-class labels are conditional until the kappa discrepancy is understood or c is obtained independently. This paper deserves a serious referee. I would send it to review, asking the authors to address the scaling problem head-on, give full details of the h->0 procedure, and ideally compute a topological invariant or report a second independent determination of c. It is a useful contribution, and the questions it leaves open are the right ones to ask.","headline":"First Lagrangian phase diagram of the Nf=1 Gross-Neveu-Wilson model with a credible qualitative picture but an unresolved scaling discrepancy that makes the central-charge labels conditional.","tokens_in":21255,"tokens_out":3699,"would_cite":true,"duration_ms":33058,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper determines the full (M, g²) phase diagram of the single-flavor Gross–Neveu–Wilson model by contracting its Grassmann tensor-network path integral, finding an Aoki phase bounded by c=1/2 critical lines that terminates at finite st","keywords":["Gross-Neveu model","Wilson fermions","Aoki phase","topological insulator","Grassmann tensor network","corner transfer matrix renormalization group","entanglement entropy","central charge"],"falsifier":"At M=0, compute the pseudoscalar condensate for g² > 0.9 with bond dimension much larger than 208 under periodic boundary conditions using an independent contraction scheme, and check whether it extrapolates to zero; also test data collapse at the claimed critical points using the predicted κ = 6/(c(√(12/c)+1)). A nonzero extrapolated condensate or a collapse requiring κ far from the predicted value would overturn the strong-coupling termination or the central-charge labels.","tokens_in":20243,"feed_emoji":"","tokens_out":5864,"duration_ms":51854,"temperature":0.7,"pith_summary":"The paper aims to establish the complete phase diagram of the single-flavor Gross–Neveu model with Wilson fermions in the Lagrangian path-integral formulation, using a Grassmann corner transfer matrix renormalization group (CTMRG) to contract the two-dimensional Grassmann tensor network without a sign problem. It claims that the parity-broken Aoki phase, detected by a non-zero pseudoscalar condensate, is enclosed by critical lines of central charge c=1/2 (two-dimensional Ising universality), while the boundary between the topological insulator and trivial phases carries c=1. A further claim is that the Aoki phase does not persist in the strong-coupling regime, terminating at triple points near (M, g²) ≈ (±0.812, 0.89), in contrast to the large-N_f prediction. If these results hold, they provide a sign-problem-free determination of the parity-broken phase for odd-flavor Wilson fermions and a numerically consistent route toward the continuum limit.","feed_headline":"Aoki phase ends at finite coupling in Gross-Neveu-Wilson model","feed_subtitle":"Tensor-network contraction maps the parity-broken Aoki phase, the topological insulator, and a c=1/2 Ising boundary — no sign problem needed","key_machinery":"The central mechanism is the Grassmann corner transfer matrix renormalization group (CTMRG): the lattice path integral is written as a uniform two-dimensional Grassmann tensor network with local bond dimension 4, and the infinite environment is approximated by corner and edge tensors truncated to bond dimension D, updated with Grassmann projectors derived from singular value decomposition. The universality classes are read off from the finite-entanglement scaling relation S_D ≈ (c/6) log ξ_D, where ξ_D is the effective correlation length obtained from the row-to-row or column-to-column transfer matrices. The pseudoscalar condensate is evaluated by inserting a local impurity Grassmann tensor,","core_discovery":"The central numerical result is the (M, g²) phase diagram of the single-flavor Gross–Neveu–Wilson model obtained by approximately contracting the two-dimensional Grassmann tensor network representation of the path integral. The Aoki phase, detected via the pseudoscalar condensate computed with an impurity tensor, is bounded by critical lines with central charge c=1/2, consistent with the two-dimensional Ising universality class. The topological insulator and trivial phases are separated by critical lines with c=1. The Aoki phase terminates at triple points near (M, g²) ≈ (±0.812, 0.89), contrary to the large-N_f phase diagram, and the topological insulating lobes are identified by a fully do","pith_inferences":["If the strong-coupling termination is confirmed, the single-flavor Aoki phase is a weak-to-intermediate coupling phenomenon, and the large-N_f phase diagram is misleading at N_f=1; a testable extension is to measure the pseudoscalar condensate at even larger D or with alternative boundary conditions to settle whether parity is broken at all g² at M=0.","The discrepancy between the fitted κ ≈ 1.48 and the predicted κ ≈ 2.03 for c=1/2 suggests the finite-entanglement scaling of this Grassmann CTMRG is not identical to the standard MPS form; if so, the quantitative c values carry an unquantified systematic error, and a cross-check using an independent method such as direct transfer-matrix spectra at fixed large D would be valuable.","The same Grassmann CTMRG pipeline could be run for N_f=2 to test whether the Aoki phase survives and whether the two-lobe phase remains topological, connecting the parity-broken phase to flavor dependence in Wilson-fermion theories.","A direct computation of the Zak phase or another topological invariant inside the two lobes would turn the entanglement-spectrum signature into a quantitative identification of the symmetry-protected topological phase and could be done with the same converged environments."],"forward_implications":["If the phase diagram is correct, the continuum limit of the single-flavor theory is approached through c=1/2 Ising critical lines (Aoki boundaries) and c=1 lines (topological/trivial boundary), and the Aoki phase is confined to finite g².","The impurity-tensor measurement of the pseudoscalar condensate provides a sign-problem-free order parameter for spontaneous Z₂ parity breaking in an odd-flavor Wilson-fermion theory.","The doubly degenerate entanglement spectrum inside the two lobes provides a practical diagnostic for the topological insulator phase without computing a topological invariant.","The triple-point location (M, g²) ≈ (±0.812, 0.89) is a concrete prediction that can be sharpened by higher-bond-dimension simulations.","The qualitative agreement with Hamiltonian-formalism results, with differences attributed to temporal doublers in the Lagrangian formulation, suggests the phase structure is robust across formalisms."],"fun_headline_variants":["Aoki phase ends at finite coupling in Gross-Neveu-Wilson","Gross-Neveu-Wilson: Aoki phase terminates at strong coupling","Tensor network maps Gross-Neveu phase diagram without sign problem","Ising critical lines bound Aoki phase in Gross-Neveu model","Single-flavor Gross-Neveu: Aoki phase fades at finite coupling"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The central-charge labels rest on the finite-entanglement scaling formula S_D ≈ (c/6) log ξ_D; Appendix C reports a fitted exponent κ ≈ 1.48 that deviates from the predicted 2.03, so if that scaling form is not valid in this Grassmann CTMRG implementation, the c=1/2 and c=1 phase-boundary assignments would not be established.","fun_headline_variants_meta":{"raw":{"variants":["Aoki phase ends at finite coupling in Gross-Neveu-Wilson","Gross-Neveu-Wilson: Aoki phase terminates at strong coupling","Tensor network maps Gross-Neveu phase diagram without sign problem","Ising critical lines bound Aoki phase in Gross-Neveu model","Single-flavor Gross-Neveu: Aoki phase fades at finite coupling"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000464,"raw_usage":{"total_tokens":2167,"prompt_tokens":770,"completion_tokens":1397,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":514,"completion_tokens_details":{"reasoning_tokens":1299}},"tokens_in":514,"tokens_out":1397,"duration_ms":9383,"temperature":1.0,"reasoning_tokens":1299,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T20:56:40.444371+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"At M=0, compute the pseudoscalar condensate for g² > 0.9 with bond dimension much larger than 208 under periodic boundary conditions using an independent contraction scheme, and check whether it extrapolates to zero; also test data collapse at the claimed critical points using the predicted κ = 6/(c(√(12/c)+1)). A nonzero extrapolated condensate or a collapse requiring κ far from the predicted value would overturn the strong-coupling termination or the central-charge labels.","supporting_citations":[],"review_version":1}