{"id":"ffd92cc9-ba0a-434f-8451-f36f16adf3f9","arxiv_id":"2412.12973","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A large-volume approximant for the constrained fermionic path integral is derived and tested in the chiral Gross-Neveu model, showing the chiral condensate appears as the edge of a flat disk in the constraint potential.","lead":"Physicists construct the constraint effective potential for fermionic condensates, using a saddle-point expansion to handle the Grassmann-valued constraint. The method locates the chiral condensate as the edge of a flat region in the potential, simplifying the standard double-extrapolation approach.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The η=0 quadratic expansion in Eq. (27) is validated only on minimizing configurations; its failure on off-constraint configurations could shift the constrained potential.","rationale":"The reader's weakest assumption identifies exactly the spot where the argument is least secure: the quadratic expansion of u around η=0, Eq. (27), is applied to all bosonic configurations in Eq. (29), but the paper validates it numerically only on the minimizing configurations of Eq. (53) in Section VD. My analysis sharpens why this matters: for a fixed configuration with M−φ of order one, the exact η-integral has its saddle at η of order V, where the omitted η³/V³ terms are not small, and the Gaussian suppression factor in Eq. (29) does not by itself guarantee that such configurations are negligible at the correct exponential rate. The moment comparison in Sec. IIIC1 shows agreement of low moments up to finite-volume corrections, but pointwise control of the log-density—needed for the central claim Ω=Γ and for locating the disk edge—requires checking the full ensemble or an analytic bound on higher cumulants. This is a genuine gap, but it is not a demonstrated failure: the constrained weight may well render off-constraint configurations irrelevant, and the large-Nf test plausibly samples the dominant sector. The paper's independent support includes exact normalization, the O(1/V) moment matching, the Ward-identity consistency check in Eq. (54), and the numerical flatting of Ω with volume, all of which make the central construction credible. The concern therefore does not overturn the reader's CONDITIONAL verdict; it strengthens the reason for conditioning on an ensemble-level test of the η=0 expansion. No ad hominem or theatrical framing is intended: the issue is a technical assumption about the validity of a per-configuration large-volume expansion, and it is explicitly localizable to Eqs. (26)–(29) and Section VD.","tokens_in":21321,"tokens_out":17894,"duration_ms":198256,"concrete_test":"For a small lattice (e.g., Ns=20 and 28) in the χGN model, generate an ensemble of ρ configurations sampled with the constrained weight exp[−Nf Seff,φ[ρ]] (plus fluctuations around the minimizing configurations of Eq. (53)). On each sampled configuration, compute numerically the exact η-integral ∫ dη0 dη1 exp[−V u(η)] using the full determinant det(Q − iη·σ/V), and compare its logarithm with the Gaussian result Eq. (28)/(29). Weight the relative error by the constrained weight and monitor its behavior across Ns=20,28,40,60; if the weighted error vanishes as V→∞ (preferably as 1/V or faster) the expansion is validated on the full ensemble. As a cheaper diagnostic, evaluate the third-order coefficient κabc of Eq. (38) on these configurations and check that V·κ·(Δη)³, with the saddle shift Δη∼Vχ⁻¹(M−φ), is subleading for configurations contributing non-negligible weight.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that Eq. (29) is a valid finite-volume representation of the constrained fermionic partition function rests on truncating u(η) to quadratic order around η=0. For a fixed bosonic configuration, the exact η-integral is controlled by the saddle satisfying ⟨M(η̄)⟩=φ, giving η̄≈−iVχ⁻¹(M(0)−φ). When M(0)−φ is O(1), η̄ is O(V), and the discarded O(η³/V³) terms contribute O(V) to the exponent; the quadratic expansion is then not a controlled approximation for that configuration. The Gaussian weight exp[−V/2(φ−M)χ⁻¹(φ−M)] suppresses such configurations only if their true contribution is exponentially suppressed at a matching rate; if the exact suppression exponent differs, or if the entropy of such configurations is large, the resulting Ω(φ) can differ. Section VD explicitly checks Re u only on the global minimizing configurations of Eq. (53), which satisfy M≈φ in the disk, and the paper's own wording concedes this validates the approximation 'at the very least for these configurations.' The moment matching in Sec. IIIC1 (corrections of order 1/V or smaller) does not by itself control pointwise convergence of the log-density for off-constraint field configurations, which is precisely what the central equality Ω=Γ requires. This is the load-bearing gap in the finite-volume justification of the method.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs a constraint effective potential for fermionic order parameters. Starting from a Dirac-delta constraint on the fermionic bilinear (Eq. 9), the authors introduce a Fourier representation of the delta, carry out the fermionic integral, and then perform a quadratic expansion of the resulting log-determinant functional u(eta) around eta = 0, justified as a large-volume expansion. This yields the approximate constrained partition function of Eq. (29), a standard path integral over bosonic fields with a Gaussian weight that suppresses configurations for which the observable M differs from the constraint phi. The paper shows that this approximation is normalized, matches the first two moments of the exact constrained distribution up to O(1/V) corrections, and argues that the resulting potential Omega(phi) coincides with the effective potential Gamma(phi) in the thermodynamic limit. The method is tested in the large-N_f chiral Gross-Neveu model, where the constrained condensate follows phi, the potential flattens for |phi| below the symmetry-breaking value, and the edge of the flat disk reproduces the standard value phi-bar. Inhomogeneous, spin-wave-like configurations are shown to dominate inside the disk.","tokens_in":21636,"tokens_out":8190,"duration_ms":81137,"significance":"If valid, the construction is significant: it is the first constraint effective potential for a continuous fermionic symmetry, and it offers a practical alternative to the double extrapolation procedure used to locate spontaneous symmetry breaking, with potential applications to QCD. The paper contains several genuine analytic checks: the approximate measure is normalized, the exactness of the constraint is tested through the constrained condensate, and the moment matching in Section IIIC1 is a nontrivial consistency condition. The numerical demonstration in the large-N_f chiral Gross-Neveu model is clean and uses established methodology. The main weakness is that the central finite-volume validity claim rests on an expansion whose control is demonstrated only on minimizing configurations, not on the full ensemble of bosonic fields that the approximate path integral must describe. The contribution is therefore promising, but the advertised generality is not yet fully supported.","major_comments":[{"comment":"The paper's central claim is that Eq. (29) is a valid finite-volume representation of the exact constrained path integral, obtained by expanding u(eta) quadratically around eta = 0. For a fixed bosonic configuration, the exact eta-integral in Eq. (24) is controlled by the saddle satisfying <M(eta-bar)> = phi; when M(0) - phi is O(1), eta-bar is O(V), and the discarded O(eta^3/V^3) terms in Eq. (26) contribute O(V) to the exponent, so the quadratic expansion is not controlled for that configuration. The Gaussian weight in Eq. (29) suppresses such configurations only if the exact suppression exponent equals the approximate one at leading order. Section VD checks Re u only on the minimizing configurations of Eq. (53), on which M is close to phi, and the paper itself concedes that this validates the approximation \"at the very least for these configurations.\" Since Eq. (29) is proposed before the bosonic integral is performed, this leaves a load-bearing gap for off-constraint bosonic configurations, and hence for finite-N_f applications such as QCD. Please either provide a uniform large-volume bound on the remainder of the quadratic expansion over the field configurations that contribute to Eq. (29), or test the Gaussian approximation against the exact one-dimensional eta-integral on a sample of non-minimizing rho configurations.","section":"Sec. IIIB, Eqs. (24)-(29)"},{"comment":"The moment matching in Table I shows that the first two moments of Z_approx^phi/Z agree with the exact distribution up to O(1/V) corrections. This does not by itself control the pointwise convergence of the log-density, i.e. of Omega(phi) = -(1/V) log Z_phi, which is a large-deviation functional. The equality Omega = Gamma in the thermodynamic limit asserted in Eq. (23) requires a large-deviation or relative-entropy statement for the approximated measure, not merely moment convergence. The numerical evidence at large-N_f mitigates this concern for the chiGN test, but the analytic claim in Section III is stronger than what is proved. Please either supply such a bound or explicitly state that the thermodynamic-limit equality is established only under the additional assumption that the Gaussian approximation of the eta-integral is uniform in the bosonic configurations.","section":"Sec. IIIC1 and Eq. (23)"},{"comment":"The numerical demonstration is entirely in the large-N_f limit, where the rho path integral collapses to the global minimizers of Eq. (53). Consequently, the numerical tests do not exercise the full ensemble of bosonic configurations that Eq. (29) must control in a finite-N_f simulation; the off-constraint-configuration error identified above could in principle shift the finite-volume potential and the edge location in Monte Carlo applications. A direct test of the approximation on configurations sampled with the weight of Eq. (29), or a finite-N_f simulation, would address this concern. As the manuscript stands, the claim of general applicability to QCD rests on uniformity of the eta = 0 expansion that is neither proved nor numerically demonstrated on the full configuration ensemble.","section":"Secs. IV and V"}],"minor_comments":[{"comment":"The exchange of the eta-integral and the epsilon to 0 limit with the path integral is formal; Appendix A justifies the finite-dimensional Grassmann case, but the main text would benefit from a sentence noting that on a finite lattice the relevant integrals are finite-dimensional and the exchange is permitted.","section":"Sec. III, Eq. (13)"},{"comment":"The statement that definiteness of Gamma implies existence of a purely imaginary solution of Eq. (17) is too quick; footnote 2 already indicates non-uniqueness at large |eta|, so the branch selection should be stated as part of the definition of the Legendre transform.","section":"Sec. IIIA, after Eq. (17)"},{"comment":"The axes for eta_0 and eta_1 are not labeled with the lattice-scale normalization, and only Re u is shown; since the phase of the integrand (Im u) matters for the Gaussian approximation away from the minimum, a statement about Im u would be useful.","section":"Fig. 5"},{"comment":"The continuum extrapolation yields a nonzero value beta = 0.000503 at fixed finite volume; the statement that this deviation vanishes in the infinite-volume limit is an expectation based on Sec. VA, not a demonstrated extrapolation, and should be presented as such.","section":"Sec. VE, Fig. 7"},{"comment":"The exact constrained partition function is a distribution rather than a positive measure, as Appendix A shows; the probabilistic language in Section IIB should be understood as heuristic, and a short caveat at the first use of \"probability density\" would prevent confusion.","section":"Appendix A and Sec. IIB"}],"recommendation":"major_revision","confidential_remarks":"The paper is inventive and the chiral Gross-Neveu test is well executed, but the advertised finite-volume representation is not fully controlled. The main gap is the uniformity of the eta = 0 quadratic expansion over bosonic configurations; the required fix may be local, e.g., adding a numerical comparison of the exact and Gaussian eta-integrals on off-constraint configurations and softening the generality claim. If that is done, the paper could be acceptable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this one if you care about how we will eventually get the chiral condensate in QCD without the painful m→0, V→∞ double limit. The paper constructs a constraint effective potential for fermionic order parameters with continuous symmetry, using a large-volume expansion of the Grassmann-valued constraint. The final expression, Eq. (29), is a standard path integral with a Gaussian weight that biases bosonic configurations toward the desired condensate, and it can be simulated. That is genuinely new; Refs [10,11] used the same Fourier representation but evaluated the characteristic function by exact diagonalization, which does not scale. The derivation is mostly transparent, and the moment-matching argument in Sec IIIC1 shows the approximated distribution reproduces the exact one up to 1/V corrections. The numerical test in the chiral Gross-Neveu model is suggestive: the constrained condensate tracks the constraint well inside the disk, the potential flattens with volume, and the edge of the disk extrapolates to the standard phibar. The relation between the flattening region and inhomogeneous spin-wave configurations is clearly demonstrated and fits the O(2) bosonic picture.\n\nThe soft spots are not fatal but they are real. The expansion of u(η) around η=0 is not the true saddle point in finite volume, and the paper justifies it by 1/V scaling plus a check of Re u only on minimizing configurations. The stress-test worry is legitimate: for off-constraint bosonic configurations where M(0)−φ is O(1), the saddle sits at η~V and the discarded O(η^3/V^3) terms contribute O(V) to the exponent. The Gaussian weight suppresses those configurations, but the paper does not prove the suppression rate matches the exact integrand; the moment matching does not control pointwise log-density convergence, which is what Ω=Γ needs. For the large-Nf test this gap is harmless because the path integral collapses to minimizers. For a Monte Carlo application to QCD it matters. The authors should either prove the suppression argument or test the approximation on off-minimizer configurations, e.g. by reweighting or comparing Eq. (29) against the exact η integral for a small system.\n\nSecond, the numerical demonstration normalizes everything by phibar from the standard method, so the edge-of-disk prediction is not fully independent. That is acceptable for a first test, but the paper should present at least one absolute, unnormalized determination of phibar from the edge. Third, the extrapolations in Figs 2 and 7 have no error bars; the fits are clean, but error propagation is missing. Fourth, the novelty discussion relative to Refs [10,11] is fair; it could be sharper about what exactly was impossible there.\n\nVerdict: this deserves a serious referee. The method is important if it survives; the central idea is sound and the demonstration is credible. I would recommend conditional acceptance after the authors tighten the off-constraint issue and add error bars. I'd take it to reading group, and I'd cite it if I worked on lattice fermions.","headline":"A genuinely new constrained-path-integral method for fermionic condensates, demonstrated in the chiral Gross-Neveu model; the quadratic-expansion gap on off-constraint configurations needs tightening, but the paper deserves peer review.","tokens_in":22118,"tokens_out":2761,"would_cite":true,"duration_ms":26832,"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":"The paper constructs the first constraint effective potential for fermionic order parameters and tests it in the chiral Gross-Neveu model.","keywords":["constraint effective potential","fermionic condensate","spontaneous symmetry breaking","chiral Gross-Neveu model","lattice field theory","saddle-point approximation","Grassmann variables","chiral condensate"],"falsifier":"Compute the neglected cubic and higher terms in the auxiliary field on configurations sampled from the unconstrained ensemble of the same theory and check whether their contribution to $Z_\\phi$ is suppressed as $V\\to\\infty$ for $|\\phi|<\\bar\\phi$; if the suppression fails, the approximated path integral in Eq.~(29) would not reproduce the exact constrained distribution.","tokens_in":21152,"feed_emoji":"⚛️","tokens_out":6207,"duration_ms":49897,"temperature":0.7,"pith_summary":"The paper constructs, for the first time, a constraint effective potential for order parameters built from fermion fields. Such a potential would give physicists a direct finite-volume handle on spontaneous symmetry breaking: the physically realized condensate is read off from the edge of the flat disk of the potential rather than from a double extrapolation in volume and explicit symmetry breaking. Because the constraint must fix a Grassmann-valued composite, the exact constrained path integral is converted into a practical form by expanding the integrand around zero auxiliary field and keeping terms up to quadratic order. The method is tested in the chiral Gross-Neveu model, where the edge of the flat region locates the chiral condensate consistently with the standard approach.","feed_headline":"Read condensates off the edge of a flat potential","feed_subtitle":"New constraint method replaces the double extrapolation by one volume limit, tested in the Gross-Neveu model.","key_machinery":"The carrying mechanism is the saddle-point, large-volume expansion of the Grassmann-valued delta constraint. Representing the Dirac delta by a Fourier integral over an auxiliary field $\\eta$ and expanding the logarithm of the shifted fermion determinant to quadratic order in $\\eta/V$ converts the constraint into a Gaussian weight: every bosonic configuration gets the extra factor $\\exp[-V(\\phi-M)\\chi^{-1}(\\phi-M)/2]$, where $M$ is the fermion-bilinear observable and $\\chi$ its fluctuations, a susceptibility matrix. This turns the constrained path integral into an ordinary path integral over bosonic fields with a modified action, Eq.~(29), which is the object a simulation can sample.","core_discovery":"The central claim is that the approximated constrained path integral in Eq.~(29) is a valid finite-volume representation of the exact constrained partition function for fermionic condensates. It is normalized, reproduces the exact moments up to $\\mathcal{O}(1/V)$ corrections, and yields $\\Omega(\\phi) = \\Gamma(\\phi)$ in the thermodynamic limit, so that the edge of the flat disk in $\\Omega$ gives the physically realized condensate. The authors demonstrate this in the chiral Gross-Neveu model in the large-$N_f$ limit, where the constrained condensate tracks the constraint value, the potential flattens inside the disk as the volume grows, and the valley edge extrapolates to the condensate obtained by the standard source method.","pith_inferences":["If Eq.~(29) survives a check on the full ensemble, it could be implemented in QCD simulations using standard determinant and trace estimates, without full diagonalization of the Dirac operator.","The comparison in Section~IIID suggests the new constraint carries more fluctuation information than binning the observable $M$ as in a density-of-states calculation, so it may resolve finite-volume chiral observables more sharply.","A direct numerical test would be to compare the exact third moment $\\langle(\\bar\\psi\\psi)^3\\rangle$ with the third moment of the approximated distribution; agreement to $\\mathcal{O}(1/V)$ would confirm the expansion on the full ensemble rather than only on energy-minimizing configurations."],"forward_implications":["In the chiral Gross-Neveu model the edge of the flat disk in the constraint effective potential gives the chiral condensate, replacing the double limit by a single infinite-volume extrapolation.","The flat region is populated by inhomogeneous, spin-wave-like condensates whose winding behavior changes at a cusp in the finite-volume potential.","The method applies to any theory with a bilinear fermion action and a real, positive fermion determinant, including QCD.","The approximated constrained distribution matches the first two moments of the exact constrained distribution exactly and the higher moments up to $\\mathcal{O}(1/V)$ corrections."],"supporting_citations":[{"why":"Origin of the constraint effective potential as a constrained path integral with a fixed order parameter.","marker":"[6]"},{"why":"Establishes the constraint effective potential and its equality to the quantum effective potential in the thermodynamic limit, the result being extended here to fermions.","marker":"[7]"},{"why":"Constrained scalar O(2) model showing the flat disk is populated by inhomogeneous configurations, the pattern reproduced for fermionic condensates.","marker":"[8]"},{"why":"Prior constrained Gross-Neveu study with discrete symmetry, providing the fermionic-constraint setup and comparison to a bosonic constraint.","marker":"[9]"},{"why":"Characteristic-function approach to constrained fermion path integrals in the Schwinger model, from which the present construction differs in using separate constraints in all order-parameter directions.","marker":"[10]"},{"why":"Probability-distribution-function method for the chiral condensate, an earlier constrained approach the paper compares against.","marker":"[11]"},{"why":"Supplies the superanalysis definition of delta distributions with Grassmann arguments used to make the constraint meaningful.","marker":"[25]"},{"why":"Textbook saddle-point approximation underlying the large-volume expansion of the auxiliary-field integral.","marker":"[26]"}],"fun_headline_variants":["First constraint effective potential for fermionic condensates","Fermionic condensates: single volume limit instead of double extrapolation","Saddle-point constraint method tested in Gross-Neveu model","Edge of flat potential yields condensate via new constraint approach","New avenue to study chiral symmetry breaking via fermionic constraints"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The quadratic expansion of the exponent in the auxiliary field must be accurate on every bosonic configuration that dominates the constrained path integral in the disk region, while the paper checks it directly only on the minimizing configurations.","fun_headline_variants_meta":{"raw":{"variants":["First constraint effective potential for fermionic condensates","Fermionic condensates: single volume limit instead of double extrapolation","Saddle-point constraint method tested in Gross-Neveu model","Edge of flat potential yields condensate via new constraint approach","New avenue to study chiral symmetry breaking via fermionic constraints"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000189,"raw_usage":{"total_tokens":1244,"prompt_tokens":764,"completion_tokens":480,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":380,"completion_tokens_details":{"reasoning_tokens":397}},"tokens_in":380,"tokens_out":480,"duration_ms":4513,"temperature":1.0,"reasoning_tokens":397,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T13:32:03.148186+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the neglected cubic and higher terms in the auxiliary field on configurations sampled from the unconstrained ensemble of the same theory and check whether their contribution to $Z_\\phi$ is suppressed as $V\\to\\infty$ for $|\\phi|<\\bar\\phi$; if the suppression fails, the approximated path integral in Eq.~(29) would not reproduce the exact constrained distribution.","supporting_citations":[{"cited_title":"Fukuda and E","cited_arxiv_id":null,"evidence_quote":"Prior constrained Gross-Neveu study with discrete symmetry, providing the fermionic-constraint setup and comparison to a bosonic constraint."},{"cited_title":"O’Raifeartaigh, A","cited_arxiv_id":null,"evidence_quote":"Characteristic-function approach to constrained fermion path integrals in the Schwinger model, from which the present construction differs in using separate constraints in all order-parameter directions."},{"cited_title":"Spontaneous symmetry breaking via inhomogeneities and the differential surface tension","cited_arxiv_id":"2109.03668","evidence_quote":"Probability-distribution-function method for the chiral condensate, an earlier constrained approach the paper compares against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the superanalysis definition of delta distributions with Grassmann arguments used to make the constraint meaningful."},{"cited_title":"Maxwell Construction for Scalar Field Theories with Spontaneous Symmetry Breaking","cited_arxiv_id":"1211.0921","evidence_quote":"Textbook saddle-point approximation underlying the large-volume expansion of the auxiliary-field integral."}],"review_version":1}