arxiv: 2604.06054 · v1 · submitted 2026-04-07 · ✦ hep-ph

Recognition: 2 theorem links

· Lean Theorem

Quarkyonic Meson Matter for Finite Isospin Density

Larry McLerran

Authors on Pith no claims yet

Pith reviewed 2026-05-10 19:18 UTC · model grok-4.3

classification ✦ hep-ph

keywords finite isospin densityquarkyonic matterlarge Nclinear sigma modelmeson condensateFermi seaQCD phases

0 comments

The pith

For isospin chemical potentials between the QCD scale and a color-number scaled multiple of it, meson matter forms a quarkyonic state with quarks bound into mesons filling a Fermi sea.

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

The paper studies QCD at finite isospin density using large number of colors. It employs a linear sigma model to describe the system at low chemical potential, where mesons form a Bose condensate. In the intermediate range from the QCD scale up to the square root of Nc times that scale, the author argues the matter is quarkyonic: quarks remain bound into mesons of size set by the strong scale yet occupy a filled Fermi sea, possibly with surface Bose condensation or finite-width Cooper pairs. This description bridges low-density confined matter and higher-density regimes without requiring deconfinement. A sympathetic reader would care because the argument identifies a new phase of dense QCD matter whose properties could influence equations of state or observable signals in systems with controlled isospin imbalance.

Core claim

For Λ_QCD ≤ μ_I ≤ √N_c Λ_QCD, meson matter is quarkyonic, with quarks bound into mesons on a size scale of order Λ_QCD corresponding to a filled Fermi sea of quarks, with possible Bose condensation at the Fermi surface and/or Cooper pairs with finite width of the surface of order Λ_QCD.

What carries the argument

The quarkyonic regime for meson matter, identified by matching the isospin chemical potential against the QCD scale and the large-Nc scaled scale inside the linear sigma model description of meson dynamics.

If this is right

At much lower μ_I the system is a simple Bose condensate of mesons.
At much higher μ_I the matter stays confined but consists of small-size mesons and Cooper pairs treatable by weak-coupling methods.
In the intermediate window the filled Fermi sea of mesons can support surface condensates or Cooper pairing of finite width.
The transition into and out of the quarkyonic window is controlled by the relative size of μ_I and √N_c Λ_QCD.

Where Pith is reading between the lines

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

The same intermediate-density logic might apply to other chemical potentials or to systems with both baryon and isospin density.
Observables such as the speed of sound or specific heat could show distinct behavior inside the window compared with pure hadronic or pure quark matter.
The picture suggests that quarkyonic behavior is a general feature of confined dense matter rather than a phenomenon tied exclusively to baryon number.

Load-bearing premise

The linear sigma model remains a valid description of meson dynamics across the low-to-intermediate density crossover, and large-Nc scaling arguments suffice to identify the quarkyonic regime without needing explicit dynamical calculations.

What would settle it

An explicit computation or lattice simulation at μ_I of order Λ_QCD that shows either the absence of a filled Fermi sea or that the typical meson size deviates substantially from the QCD scale.

Figures

Figures reproduced from arXiv: 2604.06054 by Larry McLerran.

**Figure 1.** Figure 1: FIG. 1: Evolution of [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗

read the original abstract

QCD at finite isospin density is considered for a large number of colors $N_c$. A linear sigma model is used to model the meson content of the theory at low density. At isospin chemical potential $\mu_I << \Lambda_{QCD}$, this matter forms a Bose condensate. For $\mu_I >> \sqrt{N_c} \Lambda_{QCD}$, unlike QCD remains confined, but the degrees of freedom of the system are mesons and Cooper pairs bound on size scales small compared to the QCD size scale determined by the superfluid gap. For most purposes this matter may be analyzed using weak coupling methods. For $ \Lambda_{QCD} \le \mu_I \le \sqrt{Nc} \Lambda_{QCD}$, we argue that meson matter is quarkyonic, with quarks bound into mesons on a size scale of order $\Lambda_{QCD}$ corresponding to a filled Fermi sea of quarks, with possible Bose condensation at the Fermi surface and/or Cooper pairs with finite width of the surface of order $\Lambda_{QCD}$.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

McLerran sketches a quarkyonic meson phase at intermediate isospin density via large-Nc scaling, but the identification stays at the level of argument without explicit derivations.

read the letter

McLerran takes the quarkyonic-matter framework he has used for baryon density and applies it to finite isospin chemical potential in large-Nc QCD. The paper lays out three regimes: a low-density Bose condensate of mesons described by the linear sigma model when mu_I is well below Lambda_QCD; an intermediate window where the matter is called quarkyonic, with quarks bound into mesons on the QCD scale but forming a filled Fermi sea; and a high-density regime above sqrt(Nc) Lambda_QCD where the system remains confined yet can be treated with weak-coupling methods for mesons and Cooper pairs. The separation of scales is the clearest part of the note and gives a compact picture that could be tried in neutron-star or heavy-ion modeling with isospin asymmetry. The linear sigma model is used only at low density, and the high-density side leans on standard large-Nc counting for confinement and pairing. The intermediate phase is the new claim, presented as a plausible matching between the two limits. The main limitation is that the quarkyonic identification rests on scaling arguments and the phrase “we argue,” without gap equations, effective potentials solved at finite mu_I, or any computation of fermionic propagators. The sigma model itself contains only bosonic fields, so the filled quark Fermi sea is an external interpretation rather than a result extracted from the model’s equations. No numerical checks or binding-energy calculations are shown. This is a short conceptual extension aimed at people already working on large-Nc dense matter. Readers who want a quick phase sketch for isospin density will get something to think about; those looking for quantitative predictions or model solutions will not find them here. The work is coherent on its own terms and comes from someone who has done this kind of counting before, so it deserves a serious referee even if the next version will need more explicit calculations to be useful.

Referee Report

3 major / 2 minor

Summary. The paper considers QCD at finite isospin chemical potential μ_I in the large-N_c limit. It employs a linear sigma model to describe the low-density regime (μ_I ≪ Λ_QCD), where a Bose condensate of mesons forms. For μ_I ≫ √N_c Λ_QCD it argues that the system remains confined but is dominated by mesons and Cooper pairs on scales much smaller than 1/Λ_QCD, amenable to weak-coupling analysis. The central claim is that in the intermediate window Λ_QCD ≲ μ_I ≲ √N_c Λ_QCD the matter is “quarkyonic”: quarks are bound into mesons of size ∼1/Λ_QCD while forming a filled Fermi sea, possibly with surface Bose condensation or finite-width Cooper pairing.

Significance. If substantiated, the work extends the quarkyonic-matter concept from baryon density to isospin density and supplies a concrete large-N_c phase diagram that interpolates between a low-density mesonic condensate and a high-density weakly coupled confined phase. The use of the linear sigma model plus standard large-N_c counting is a strength, but the absence of explicit dynamical calculations limits immediate falsifiability.

major comments (3)

[Abstract / intermediate-density argument] Abstract and the paragraph introducing the intermediate-density regime: the identification of a filled quark Fermi sea with p_F ∼ μ_I is asserted via large-N_c scaling alone, without an explicit derivation of an effective quark chemical potential, a fermionic propagator, or a binding-energy calculation inside the linear sigma model. Because the model contains only bosonic fields (σ, π), the mapping from meson dynamics to a quark Fermi surface must be shown rather than inferred from dimensional analysis.
[Intermediate-density regime] The claim that meson size remains O(1/Λ_QCD) while μ_I reaches √N_c Λ_QCD is load-bearing for the quarkyonic identification, yet no gap equation, effective potential minimization, or numerical check at finite μ_I is presented to confirm that the condensate scale does not shrink or that the meson dispersion remains gapped on the Λ_QCD scale.
[Phase-structure discussion] The transition points Λ_QCD and √N_c Λ_QCD are stated without a quantitative matching condition between the low-density linear-sigma-model description and the high-density weak-coupling regime; an explicit overlap region or continuity argument for the order parameters would be required to establish the window.

minor comments (2)

Notation: √N_c versus √Nc appears inconsistently; standardize to √N_c throughout.
The phrase “for most purposes this matter may be analyzed using weak coupling methods” in the high-density regime would benefit from a brief statement of the relevant expansion parameter (e.g., g(μ_I) or 1/N_c).

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments. We address each major comment point by point below, providing clarifications based on the large-N_c framework of the manuscript and indicating revisions that will be incorporated to strengthen the arguments.

read point-by-point responses

Referee: Abstract / intermediate-density argument: Abstract and the paragraph introducing the intermediate-density regime: the identification of a filled quark Fermi sea with p_F ∼ μ_I is asserted via large-N_c scaling alone, without an explicit derivation of an effective quark chemical potential, a fermionic propagator, or a binding-energy calculation inside the linear sigma model. Because the model contains only bosonic fields (σ, π), the mapping from meson dynamics to a quark Fermi surface must be shown rather than inferred from dimensional analysis.

Authors: The linear sigma model is the effective bosonic theory obtained from the underlying quark-gluon dynamics in the large-N_c limit, where the σ and π fields represent quark bilinears. The isospin chemical potential μ_I enters through the covariant derivative on the meson fields, which in the quark picture corresponds to opposite chemical potentials for up and down quarks. This induces a filled Fermi sea with p_F ∼ μ_I while the binding remains on the Λ_QCD scale set by the model parameters. Although the manuscript relies on this standard large-N_c correspondence rather than a new microscopic derivation, we agree that an explicit mapping paragraph will improve clarity. We will revise the introduction to derive the effective quark chemical potential from the large-N_c counting and the structure of the linear sigma model. revision: yes
Referee: [Intermediate-density regime] The claim that meson size remains O(1/Λ_QCD) while μ_I reaches √N_c Λ_QCD is load-bearing for the quarkyonic identification, yet no gap equation, effective potential minimization, or numerical check at finite μ_I is presented to confirm that the condensate scale does not shrink or that the meson dispersion remains gapped on the Λ_QCD scale.

Authors: We acknowledge that an explicit minimization of the effective potential at finite μ_I would make the scale stability more transparent. In the linear sigma model the vacuum parameters are fixed by zero-density phenomenology, and the chemical-potential term does not drive the condensate scale below Λ_QCD for μ_I ≲ √N_c Λ_QCD because of the 1/N_c suppression of quantum corrections and the form of the potential. The meson dispersion remains gapped throughout this window. To address the comment we will add a short subsection (or appendix) that minimizes the effective potential at finite μ_I and confirms that both the vev and the meson masses stay O(Λ_QCD). revision: yes
Referee: [Phase-structure discussion] The transition points Λ_QCD and √N_c Λ_QCD are stated without a quantitative matching condition between the low-density linear-sigma-model description and the high-density weak-coupling regime; an explicit overlap region or continuity argument for the order parameters would be required to establish the window.

Authors: The boundaries are determined by where each description ceases to be valid: the linear sigma model applies for dilute mesonic matter (μ_I ≪ Λ_QCD) while weak-coupling methods become reliable once the effective coupling is small (μ_I ≫ √N_c Λ_QCD). We agree that an explicit continuity argument for the order parameters would strengthen the case for a well-defined intermediate window. We will expand the phase-structure section to include a qualitative matching discussion, showing that the meson condensate and superfluid gap vary smoothly across the overlap region on the basis of large-N_c scaling. revision: yes

Circularity Check

0 steps flagged

No circularity: regime identification rests on standard large-Nc scaling applied to an effective bosonic model without self-referential reduction or fitted predictions.

full rationale

The paper's central argument proceeds by dividing the μ_I range into three regimes and invoking the linear sigma model only for the low-density Bose condensate while using large-Nc scaling to characterize the intermediate and high-density regimes. No equations appear that define a quantity in terms of itself, rename a fit as a prediction, or reduce the quarkyonic claim to a self-citation chain. The identification of a filled quark Fermi sea is presented as an interpretive consequence of scaling rather than a derived output from the model's equations of motion; this is an external physical argument, not a circular construction. The derivation chain is therefore self-contained against external benchmarks and receives the default non-circularity finding.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 1 invented entities

The central claim rests on the large-Nc limit simplifying the theory and on the linear sigma model capturing meson physics at low density; the quarkyonic phase itself is introduced via scaling without independent dynamical evidence.

axioms (2)

domain assumption The large-Nc limit of QCD applies and controls the phase structure at finite isospin density
Invoked to justify weak-coupling analysis at high mu_I and to define the sqrt(Nc) Lambda_QCD scale separating regimes.
domain assumption The linear sigma model provides an accurate effective description of meson dynamics for mu_I much less than Lambda_QCD
Used to establish the Bose-condensate phase at low density and to anchor the low-density side of the intermediate regime.

invented entities (1)

Quarkyonic meson matter no independent evidence
purpose: To label the intermediate-density phase in which quarks remain confined into mesons yet occupy a filled Fermi sea on the QCD scale
Postulated on the basis of scaling arguments connecting low- and high-density regimes; no independent falsifiable signature is supplied in the abstract.

pith-pipeline@v0.9.0 · 5480 in / 1744 out tokens · 59011 ms · 2026-05-10T19:18:44.920606+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

A linear sigma model is used to model the meson content... For Λ_QCD ≤ μ_I ≤ √Nc Λ_QCD, we argue that meson matter is quarkyonic... filled Fermi sea of quarks
IndisputableMonolith/Foundation/ArithmeticFromLogic.lean embed_injective unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

ρI_Q(k) = 1/(2Nc) ∫ KM(k−p/2) ρI_M(p) d³p (Kojo Filling Criteria)

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

16 extracted references · 16 canonical work pages

[1]

D. T. Son and M. A. Stephanov, Phys. Rev. Lett.86, 592 (2001), arXiv:hep- ph/0005225

work page arXiv 2001
[2]

J.B.KogutandD.K.Sinclair,Phys.Rev.D66,014508(2002),arXiv:hep-lat/0201017. 18

work page arXiv 2002
[3]

J. B. Kogut and D. K. Sinclair, Nucl. Phys. B Proc. Suppl.119, 556 (2003), arXiv:hep- lat/0209054

work page arXiv 2003
[4]

B. B. Brandt, G. Endrodi, and S. Schmalzbauer, EPJ Web Conf.175, 07020 (2018), arXiv:1709.10487 [hep-lat]

work page arXiv 2018
[5]

B. B. Brandt, F. Cuteri, and G. Endrodi, JHEP07, 055 (2023), arXiv:2212.14016 [hep-lat]

work page arXiv 2023
[6]

Abbott, W

R. Abbott, W. Detmold, F. Romero-López, Z. Davoudi, M. Illa, A. Parreño, R. J. Perry, P. E. Shanahan, and M. L. Wagman (NPLQCD), Phys. Rev. D108, 114506 (2023), arXiv:2307.15014 [hep-lat]

work page arXiv 2023
[7]

Abbott, W

R. Abbott, W. Detmold, M. Illa, A. Parreño, R. J. Perry, F. Romero-López, P. E. Shanahan, and M. L. Wagman (NPLQCD), Phys. Rev. Lett.134, 011903 (2025), arXiv:2406.09273 [hep-lat]

work page arXiv 2025
[8]

Ayala, A

A. Ayala, A. Bandyopadhyay, R. L. S. Farias, L. A. Hernández, and J. L. Hernández, Phys. Rev. D107, 074027 (2023), arXiv:2301.13633 [hep-ph]

work page arXiv 2023
[9]

M. G. Alford, K. Rajagopal, and F. Wilczek, Phys. Lett. B422, 247 (1998), arXiv:hep- ph/9711395

work page arXiv 1998
[10]

Instantons in QCD

T. Schäfer and E. V. Shuryak, Rev. Mod. Phys.70, 323 (1998), arXiv:hep-ph/9610451

work page Pith review arXiv 1998
[11]

Fujimoto, Phys

Y. Fujimoto, Phys. Rev. D109, 054035 (2024), arXiv:2312.11443 [hep-ph]

work page arXiv 2024
[12]

D. T. Son, Physical Review D59(1999), 10.1103/physrevd.59.094019

work page doi:10.1103/physrevd.59.094019 1999
[13]

Fujimoto, T

Y. Fujimoto, T. Kojo, and L. D. McLerran, Phys. Rev. Lett.132, 112701 (2024), arXiv:2306.04304 [nucl-th]

work page arXiv 2024
[14]

McLerran and R

L. McLerran and R. D. Pisarski, Nucl. Phys. A796, 83 (2007), arXiv:0706.2191 [hep- ph]

work page arXiv 2007
[15]

Kojo, Phys

T. Kojo, Phys. Rev. D104, 074005 (2021), arXiv:2106.06687 [nucl-th]

work page arXiv 2021
[16]

McLerran and S

L. McLerran and S. Reddy, Phys. Rev. Lett.122, 122701 (2019), arXiv:1811.12503 [nucl-th]

work page arXiv 2019