{"id":"3b5da7ec-e29c-40d3-b08f-9cec8624558e","arxiv_id":"hep-ph/0406159","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A two-colour Nambu-Jona-Lasinio model whose quark-quark and quark-antiquark couplings are tied by colour-current Fierz transformation reproduces the chiral and diquark condensates, baryon density and pion mass seen in two-colour lattice QCD.","lead":"A two-colour, two-flavour quark model with a single contact coupling reproduces lattice-QCD data for the transition into a diquark-condensed phase at finite density. The paper benchmarks how far a simple quasiparticle model can track the thermodynamics of a tractable gauge theory.","discovery_kind":"extension","skeptic_critique":{"model":"selfhosted/deepseek-v4.1-flash","headline":"The 'very good agreement' claim (abstract:s4) overlays a single model curve at m0 = 4.5 MeV on lattice points from ref. [27] that span several quark masses; the paper never states which lattice mass matches the model, so the comparison cannot separate agreement from coincidence.","rationale":"The paper is a careful mean-field study with a clearly defined model, parameters fixed against hadronic observables rather than the lattice data, and acknowledged artifacts (the |Delta| turnover near Lambda, the density deficit at large mu). The Fierz coefficients G = H = 3/2 G_c are correctly derived in the appendix, and the mu_c = m_pi/2 result is a genuine consistency check of the Pauli-Gursey structure, honestly discussed in §5. The reader's CONDITIONAL verdict is appropriate. My concern sharpens the reader's first point: the 'very good agreement' claim is not merely unquantified, it is a comparison of a single-mass model curve against multi-mass lattice data, and the paper does not establish the mapping. This is fixable by replotting at matched masses. The reader's second point (mu_c = m_pi/2 partly guaranteed) is correct but the paper acknowledges it in §5, so it is a framing issue rather than a correctness issue. The reader's third point (private communication [57]) is a real but minor concern: the justification for dropping coloured qqbar modes rests on an unpublished communication, but the truncation is a standard NJL approximation and the empirical comparison does not depend on it. I would keep the verdict at CONDITIONAL.","tokens_in":17074,"tokens_out":17289,"duration_ms":329948,"concrete_test":"Recompute the model curves of Figs. 2, 5 and 6 at each of the lattice quark masses used in ref. [27] (or state explicitly the mapping between the model's m0 = 4.5 MeV and the lattice masses) and overlay each curve on the corresponding lattice points. If the model curve at the matched mass deviates from the lattice points by more than the lattice statistical error, the 'very good agreement' claim in abstract:s4 is overstated and should be qualified. A secondary check: verify from ref. [27] that the positive-determinant sector of SU(2) adjoint QCD has the same Goldstone spectrum (five Goldstones, SU(4) to Sp(4)) and the same m_pi/f_pi ratio as fundamental two-flavour QCD at the simulated parameters.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central empirical claim is abstract:s4, 'Very good agreement is found with lattice QCD (N_c=2) results,' supported by Figs. 2, 5 and 6, which overlay model curves on data from ref. [27]. Two features of that comparison are load-bearing and unexamined. First, the model is computed at a single current quark mass, m0 = 4.5 MeV (Table 1), while the lattice points in Figs. 2, 5 and 6 are explicitly labelled as corresponding to 'different values for the quark masses' (open circles, squares, diamonds). The paper never states which lattice quark mass matches m0 = 4.5 MeV, nor does it show the model curve at the other masses. If the lattice points at different masses do not collapse onto a single curve, then a single model curve cannot be in 'very good agreement' with all of them; it can at best pass through the middle of a spread. Second, the comparison is to SU(2) adjoint staggered fermions, not fundamental quarks; the bridge is one sentence in §4 ('the positive determinant sector behaves like a two-flavour theory') with no demonstration that the Goldstone spectrum and condensate physics coincide at the quantitative level required. Neither point is acknowledged in the text. A related missing support is the private communication [57], cited in the conclusions to justify dropping the coloured qqbar modes that would otherwise break G = H; that truncation is what enforces the Pauli-Gursey structure, so its justification should not rest on an unpublished note.","agreement_with_reader":"partial"},"referee_report":{"model":"selfhosted/deepseek-v4.1-flash","summary":"The paper studies a two-flavour, two-colour Nambu–Jona-Lasinio model at finite temperature and baryon chemical potential, with the quark-antiquark and diquark couplings tied by the Fierz transformation of the colour-current interaction (G = H = 3G_c/2, Eq. (4)). The three model parameters (Λ = 0.78 GeV, G = 10.3 GeV^-2, m0 = 4.5 MeV) are fixed by f_π, the quark condensate and GMOR (Table 1). Minimizing the thermodynamic potential (Eq. (22)) with the gap equations (Eq. (23)), the authors find a second-order transition at µ_c = m_π/2, a nonzero diquark condensate above µ_c, a correlated drop of the chiral condensate, and a linear rise m_π = 2µ. They compare the T = 0 condensates, baryon density and pion mass to lattice data from ref. [27] (Figs. 2, 5, 6) and to leading-order chiral effective field theory, and discuss the chiral limit (§5).","tokens_in":17373,"tokens_out":15455,"duration_ms":264990,"significance":"The paper provides a quantitatively successful mean-field description of N_c = 2 QCD thermodynamics, with parameters fixed by hadronic observables rather than tuned to the lattice. The G = H relation from the Fierz transform is a principled constraint that ties the chiral and diquark channels, and the resulting Pauli-Gursey symmetry is a genuine feature of two-colour QCD. The comparison to lattice data and chiral EFT is valuable, and the chiral-limit analysis (§5) is instructive. The paper is a useful benchmark for NJL-type quasiparticle models.","major_comments":[{"comment":"The model is computed at a single current quark mass m0 = 4.5 MeV (Table 1), while the lattice points in these figures are explicitly labelled as corresponding to different quark masses. The paper does not state which lattice mass matches m0 = 4.5 MeV, nor does it show the model curve at the other masses, nor does it give a quantitative measure of agreement. If the lattice points at different masses do not collapse onto a single curve, a single model curve cannot be in 'very good agreement' with all of them. Please state the mapping and/or show the model at the lattice masses.","section":"Figs. 2, 5, 6; abstract:s4","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The 'very good agreement' is demonstrated only at T = 0 (Figs. 2, 5, 6). The finite-temperature phase diagram (Fig. 4) and the finite-T condensates (Figs. 3, 10) are not compared to lattice data, although lattice studies of the N_c = 2 phase diagram are cited (refs. [13, 14, 37–40]). The abstract should be qualified to T = 0, or the comparison extended.","section":"Abstract:s4","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The mapping from adjoint staggered fermions to two-flavour fundamental QCD rests on one sentence ('the positive determinant sector behaves like a two-flavour theory'). Since this is a cited published result, I do not ask for a re-proof, but the paper should state explicitly what is being compared and what residual systematic remains.","section":"§4, ref. [27]","authors_acknowledge":null,"authors_acknowledge_p":null}],"minor_comments":[{"comment":"The sharp cutoff θ(Λ²−p²) multiplies the thermal terms as well as the vacuum term; please clarify or justify this choice, since it can affect the finite-T results.","section":"Eq. (22)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"f_π = 89.6 MeV vs the empirical 92.4 MeV; please comment on the 3% discrepancy and its effect on the comparison.","section":"Table 1","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The justification for dropping colour-triplet qqbar modes rests on a private communication; please provide a published reference or a quantitative estimate.","section":"Conclusions, ref. [57]","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"'Very good agreement' is qualitative; please give a quantitative measure (e.g., χ² or a table of values) to support the claim.","section":"Abstract:s4","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The dashed line m_π = 2µ is the chiral EFT prediction; please clarify the range of validity and the expected deviations.","section":"Fig. 6","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The mixing of σ, Δ, Δ* at µ > µ_c is discussed; please give the explicit form of the mass matrix elements or the eigenmodes.","section":"§4.1","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"Please state explicitly whether the transition at µ = 0 in the exact chiral limit is first order, since the potential depends only on R² = σ² + |Δ|².","section":"§5","authors_acknowledge":null,"authors_acknowledge_p":null}],"recommendation":"minor_revision"},"author_rebuttal":null,"desk_editor":{"model":"selfhosted/deepseek-v4.1-flash","letter":"Ratti and Weise take the two-colour NJL model with the quark-antiquark and diquark couplings tied by the Fierz transform of the colour-current interaction (G = H = 3G_c/2, eq. 4) and run it against the Hands et al. lattice data. The genuinely new content is the quantitative overlay: chiral and diquark condensates, baryon density and pion mass compared point-by-point with N_c = 2 lattice results, plus a mass-spectrum section that solves the mixed sigma-diquark-antidiquark matrix above mu_c and identifies the massless U(1)_B Goldstone mode. That last piece is a real addition, not a rehash of Kogut et al. The parameter fixing is done the right way round: f_pi, the quark condensate and GMOR, none of them tuned to the lattice observables being compared. The known artifacts, the diquark condensate turning over near the cutoff and the baryon density falling short at large mu, are stated in the text rather than buried. That counts for something. The soft spots are reporting gaps, not structural errors. First, the 'very good agreement' of abstract:s4 rests on a single model curve at m0 = 4.5 MeV laid over lattice points that the captions themselves say correspond to different quark masses. The paper never says which lattice mass matches 4.5 MeV, and there is no error band anywhere. As printed you cannot separate agreement from a curve passing through the middle of a spread. Second, the bridge from SU(2) adjoint staggered fermions to a two-flavour theory is one sentence in section 4; that is the load-bearing identification for the whole comparison and it deserves a paragraph. Third, the truncation that drops the coloured quark-antiquark modes, and thereby enforces Pauli-Gursey, is justified by a private communication [57]. That is the one place where an unpublished note carries real weight. On the reader's circularity worry: mu_c = m_pi/2 does follow from the imposed G = H plus a model-computed m_pi, so it is a consistency check of the symmetry pattern rather than an independent prediction. The paper says 'confirm', which is fair; I would not call this a flaw, just a framing point. Who should read it: anyone modelling colour superconductivity with NJL, and anyone using N_c = 2 lattice data as a benchmark. It deserves referees, mainly to force the mass-matching and error-band issues into the open.","headline":"A careful, honestly bounded mean-field study that puts a Pauli-Gursey-symmetric two-colour NJL model against the N_c = 2 lattice data; the physics is sound, but the headline 'very good agreement' is looser than the figures can support as printed.","tokens_in":18074,"tokens_out":1415,"duration_ms":21268},"pith_extraction":{"msc":[],"pacs":["12.38.Aw","11.30.Rd","12.38.Mh","21.65.+f"],"model":"selfhosted/deepseek-v4.1-flash","headline":"A two-colour quark model with one coupling reproduces the lattice thermodynamics of dense matter, including a transition to a diquark condensate at half the pion mass.","keywords":["two-colour QCD","Nambu-Jona-Lasinio model","diquark condensation","colour superconductivity","finite chemical potential","chiral symmetry restoration","Pauli-Gürsey symmetry","lattice QCD comparison"],"falsifier":"A two-colour lattice calculation at $μ > 2μ_c$ that finds the scaled baryon density departing from the model curve in the same direction and by the same amount as the mean-field result would support the picture, while one that finds the density following the mean-field curve where the model currently underestimates it would show that the discrepancy is not, as claimed, a mean-field artifact; likewise, a lattice measurement of the lightest scalar mass above $μ_c$ that does not vanish would break the baryon-number Goldstone identification.","tokens_in":18335,"feed_emoji":"⚛️","tokens_out":3195,"duration_ms":53848,"temperature":0.7,"pith_summary":"The paper asks whether the thermodynamics of QCD with two colours can be captured by a Nambu–Jona-Lasinio model: a quark model in which gluons are integrated out and replaced by a local four-fermion contact interaction. The distinguishing move is that the strength of the quark–antiquark channel and the strength of the diquark channel are not tuned independently; a Fierz rearrangement of a single colour-current interaction fixes them equal, $G = H = 3G_c/2$, which is exactly the Pauli-Gürsey symmetry relating pions and scalar diquarks in the two-colour theory. Minimizing the mean-field thermodynamic potential then yields a second-order transition at $μ_c = m_π/2$, above which scalar diquarks Bose-condense and the diquark condensate grows while the chiral condensate shrinks. The authors overlay their curves on lattice data for the two condensates, the baryon density and the pion mass, and report good agreement below about twice the critical chemical potential. The payoff, if the claim holds, is that the low-energy quasiparticle physics of a dense gauge theory is set by symmetry and a single coupling rather than by the details of gluon exchange.","feed_headline":"Two-colour QCD's diquark transition falls out of one coupling","feed_subtitle":"A single-strength quark model reproduces the lattice condensates, the baryon density and a pion mass rising as 2μ.","key_machinery":"The load-bearing object is the identification $G = H = 3G_c/2$ obtained by Fierz-transforming the colour-current interaction $-G_c Σ_a J^a_μ J^{μa}$ into colour-singlet quark–antiquark and scalar-diquark channels. This single relation is what makes the model predictive rather than fitted: it forces the pion and the scalar diquark to be related by the Pauli-Gürsey symmetry, so the diquark condensate appears at exactly the chemical potential where a half pion mass can be supplied, $μ_c = m_π/2$. The machinery that turns this into numbers is the Nambu–Gorkov propagator for the bispinor of quarks and charge-conjugated quarks, the mean-field thermodynamic potential of eq. (22) with a sharp three-momentum cutoff $Λ = 0.78$ GeV, and a second-order expansion of the effective action whose mass matrix in the $(π, σ, Δ, Δ^*)$ basis yields the mixed eigenmodes once $|Δ| ≠ 0$.","core_discovery":"The central claim is that a two-flavour, two-colour NJL model whose quark–antiquark and diquark couplings are tied together by the Fierz transform of the colour-current interaction reproduces the thermodynamics that lattice simulations of $N_c = 2$ QCD observe. Concretely: a second-order phase transition at $μ_c = m_π/2$; a nonzero scalar-diquark condensate for $μ > μ_c$ that tracks the fall of the chiral condensate; a baryon charge density in agreement with lattice data for $μ < 2μ_c$; and a pion mass that rises linearly as $m_π = 2μ$ in the dense phase, matching both the lattice points and the leading-order chiral effective Lagrangian prediction. The authors also work out the bosonic spectrum: below $μ_c$ the diquark and antidiquark masses split as $m_π ∓ 2μ$ according to their baryon number, and above $μ_c$ the $σ$, diquark and antidiquark fields mix, leaving one exactly massless mode — the Goldstone boson of spontaneously broken baryon number — alongside a massive mode following $2μ√(1 + 3(m_π/2μ)^4)$. In the chiral limit the model exhibits the Pauli-Gürsey degeneracy along $σ^2 + |Δ|^2 = $ const, so that the ordering of the two limits $m_0 → 0$ and $μ → 0$ matters.","pith_inferences":["If the single-coupling relation is what buys the agreement, then the same Fierz bookkeeping should predict the ratio of diquark to chiral observables in any theory where the colour-current interaction dominates — a check that does not require a lattice, only a well-measured meson and diquark spectrum.","The authors note that $|Δ|$ turns over near $μ ∼ Λ$ and call that an artifact of the cutoff; a natural test is to repeat the calculation with a smooth or larger cutoff and see whether the turnover moves with $Λ$ or with a physical scale.","The underestimate of the baryon density beyond $μ ≈ 2μ_c$ is attributed to missing quasiparticle correlations; a beyond-mean-field treatment should therefore improve the density first and the condensates later, which is an ordering that a future calculation could confirm or refute.","Because the transition is second order and driven by a conserved charge, the same structure — a Goldstone mode of broken baryon number appearing at half the meson mass — plausibly recurs in any theory where diquarks are the lightest charged excitations, which makes it a template rather than a two-colour curiosity."],"forward_implications":["The critical chemical potential for the onset of dense matter in two-colour QCD is fixed by the pion mass alone, $μ_c = m_π/2$, with no independent diquark coupling to fit.","In the dense phase the pion mass grows as $m_π = 2μ$, so a measurement of the pion mass at finite density is a direct probe of the condensate's BEC-type nature.","Baryon-number superconductivity in this theory has a definite signature: one scalar mode becomes exactly massless above $μ_c$ while its partner is pushed up to $2μ√(1 + 3(m_π/2μ)^4)$.","The chiral and diquark condensates are not independent order parameters in the chiral limit; they rotate into one another along a fixed radius, so a small quark mass or a small chemical potential is enough to pick the diquark direction.","The model's agreement with lattice data up to $μ ≈ 2μ_c$ and its failure beyond that range locates where mean-field quasiparticles stop being an adequate description of dense two-colour matter."],"supporting_citations":[{"why":"Supplies the two-colour lattice data for the chiral and diquark condensates, the baryon density and the pion mass against which the model curves are overlaid.","marker":"[27]"},{"why":"Provides the leading-order chiral effective Lagrangian predictions — $μ_c = m_π/2$, $m_π = 2μ$, the $√(1-(m_π/2μ)^4)$ condensate form and the massive-mode formula — used as the analytic baseline.","marker":"[11]"},{"why":"Supplies the general argument that the transition to finite baryon density occurs at the lowest energy per quark, namely half the lightest diquark mass.","marker":"[12]"},{"why":"Background review of lattice QCD thermodynamics cited for the overall state of the field.","marker":"[1]"},{"why":"Original Nambu–Jona-Lasinio construction on which the four-fermion quasiparticle model is built.","marker":"[42]"},{"why":"Source for the form of the second-order effective action and the fluctuation matrix used to extract the meson and diquark masses.","marker":"[56]"},{"why":"Earlier lattice study that reported mixing between the $σ$, diquark and antidiquark modes once the diquark condensate is nonzero, which the mass-matrix analysis here reproduces.","marker":"[39]"},{"why":"Mean-field lattice study of the $(T, μ)$ phase diagram of two-colour QCD, the earlier framework this work compares against.","marker":"[13]"}],"fun_headline_variants":["One Fierz-tied coupling reproduces two-colour QCD diquark thermodynamics","Diquark condensate switches on at half the pion mass in two-colour QCD","Tied quark and diquark couplings match N_c=2 lattice thermodynamics","Broken baryon number leaves one massless mode above μ_c","NJL model tracks lattice diquark condensate, baryon density and pion mass"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The argument rests on replacing the local colour gauge dynamics by a global-symmetry contact interaction with one coupling and a sharp momentum cutoff, and on treating the auxiliary fields in mean field; if that replacement misses the relevant physics at the densities probed, the agreement with lattice data is a coincidence of the fit rather than a derivation.","fun_headline_variants_meta":{"raw":{"variants":["One Fierz-tied coupling reproduces two-colour QCD diquark thermodynamics","Diquark condensate switches on at half the pion mass in two-colour QCD","Tied quark and diquark couplings match N_c=2 lattice thermodynamics","Broken baryon number leaves one massless mode above μ_c","NJL model tracks lattice diquark condensate, baryon density and pion mass"]},"model":"selfhosted/deepseek-v4.1-flash","effort":"low","attempts":1,"cost_usd":0.0,"raw_usage":{"total_tokens":1358,"prompt_tokens":1009,"completion_tokens":349,"prompt_tokens_details":null,"completion_tokens_details":{"reasoning_tokens":256}},"tokens_in":1009,"tokens_out":349,"duration_ms":5309,"image_pages":[],"temperature":1.0,"reasoning_budget":50,"reasoning_tokens":256,"transport_provider":"selfhosted","cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-11-13T08:40:37.259776+00:00","model_set":{"reader":"selfhosted/deepseek-v4.1-flash"},"pipeline_version":"pith-pipeline@v1.5.0-archive-selfhost","input_kind":null,"input_chars":null,"input_total_chars":null,"passes_missing":null,"published_standing":null,"published_confidence":null,"objection_audit":null,"run_record":null,"basis":{"kind":"not_recorded","source":null,"chars_given":null,"chars_total":null},"receipt":null,"falsifier":"A two-colour lattice calculation at $μ > 2μ_c$ that finds the scaled baryon density departing from the model curve in the same direction and by the same amount as the mean-field result would support the picture, while one that finds the density following the mean-field curve where the model currently underestimates it would show that the discrepancy is not, as claimed, a mean-field artifact; likewise, a lattice measurement of the lightest scalar mass above $μ_c$ that does not vanish would break the baryon-number Goldstone identification.","supporting_citations":[{"cited_title":"Eur. Ph ys. J. C22, 451","cited_arxiv_id":null,"evidence_quote":"Supplies the two-colour lattice data for the chiral and diquark condensates, the baryon density and the pion mass against which the model curves are overlaid."},{"cited_title":"Nucl. Phys. B582, 477","cited_arxiv_id":null,"evidence_quote":"Provides the leading-order chiral effective Lagrangian predictions — $μ_c = m_π/2$, $m_π = 2μ$, the $√(1-(m_π/2μ)^4)$ condensate form and the massive-mode formula — used as the analytic baseline."},{"cited_title":"a nd J. J. M. Verbaarschot, Phys. Rev. D58, 096007","cited_arxiv_id":null,"evidence_quote":"Supplies the general argument that the transition to finite baryon density occurs at the lowest energy per quark, namely half the lightest diquark mass."},{"cited_title":"Lecture Notes in Physics (Springer) 583, 209","cited_arxiv_id":null,"evidence_quote":"Background review of lattice QCD thermodynamics cited for the overall state of the field."},{"cited_title":"Phys. Rev. 122, 345","cited_arxiv_id":null,"evidence_quote":"Original Nambu–Jona-Lasinio construction on which the four-fermion quasiparticle model is built."},{"cited_title":"Abnormal number of Nambu-Goldstone bosons in the color-asymmetric 2SC phase of an NJL-type model","cited_arxiv_id":"hep-ph/0403151","evidence_quote":"Source for the form of the second-order effective action and the fluctuation matrix used to extract the meson and diquark masses."},{"cited_title":"Phys. Rev. D68, 054507","cited_arxiv_id":null,"evidence_quote":"Earlier lattice study that reported mixing between the $σ$, diquark and antidiquark modes once the diquark condensate is nonzero, which the mass-matrix analysis here reproduces."},{"cited_title":"Phys. Lett. B169, 421","cited_arxiv_id":null,"evidence_quote":"Mean-field lattice study of the $(T, μ)$ phase diagram of two-colour QCD, the earlier framework this work compares against."}],"lane_accuracy":{"status":"not yet measured","lane":{"model":"selfhosted/deepseek-v4.1-flash","input_kind":null},"url":"https://pith.science/accuracy"},"review_version":1}