{"id":"edc06214-a699-4e41-9fec-0d444ac17a85","arxiv_id":"2501.15242","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A large dark-quark chemical potential in a PNJL toy model makes the dark-QCD chiral transition first-order and shifts its gravitational-wave signal into the BBO band for T_c between 1 and 100 GeV.","lead":"Dark quarks with a large chemical potential can turn the dark-QCD chiral transition into a first-order one, releasing gravitational waves in a band that the proposed BBO detector could see. The paper maps this transition in a PNJL model and finds that a higher chemical potential lengthens the transition and boosts the signal.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The first-order transition and BBO-detectable GW signal rest on treating the Polyakov loop's l and l* as independent real fields; if the physical l*=conj(l) constraint removes the barrier, the central claim collapses. This needs a constrained re-computation.","rationale":"The reader's weakest assumption and my most load-bearing concern coincide: the phase diagram that produces the first-order chiral transition is computed in an enlarged field space where l and l* are independent real fields. The paper acknowledges the physical constraint but bypasses it because the constrained potential is complex. This is a structural vulnerability rather than a numerical quibble: the transition barrier, and therefore the entire gravitational-wave prediction, arises from the Polyakov-loop sector in this unphysical parametrization. The proposed test directly probes whether the barrier survives on the physical submanifold, which would settle whether the central claim is real or an artifact of the ansatz. I do not treat the paper as dishonest; the issue is an acknowledged approximation inherited from earlier PNJL literature, but its consequences here are unusually large because the entire signal depends on it. The Z_sigma positivity patch is also concerning, but it affects the strength and duration of the transition rather than its existence, so I regard it as secondary. The reader's CONDITIONAL verdict already captures this uncertainty, so no verdict change is needed.","tokens_in":28191,"tokens_out":5243,"duration_ms":54515,"concrete_test":"Recompute the phase diagram of Section III.B with the Polyakov loop constrained to its physical form: L = diag(e^{i theta1}, e^{i theta2}, e^{i theta3}) with theta1+theta2+theta3 = 0, l = (e^{i theta1}+e^{i theta2}+e^{i theta3})/3, l* = conj(l), and minimize the real part of the effective potential over (theta1, theta2, sigma) for the same GS*Lambda^2 and mu/T0 values as Fig. 2. If no first-order jump in sigma survives at any mu/T0 up to the edge of the phase diagram, or if the critical point moves outside the claimed range, the central claim fails. A cheaper variant: at a claimed first-order point, e.g., GS*Lambda^2 = 2.2 and mu/T0 = 1.65, compute the effective potential along the straight-line path between the two minima restricted to the physical submanifold; if no barrier exists there, the three-dimensional-field-space barrier is an artifact.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—a first-order dark chiral transition driven by chemical potential, with BBO-detectable GWs—rests on the phase diagram of Section III.B. That diagram is obtained by treating l and l* as independent real fields (just before Eq. 3.14), after the paper itself notes they should be complex conjugates and that the effective potential then becomes complex. This is not a harmless parametrization: it enlarges the configuration space from the physical two real degrees of freedom of a complex Polyakov loop to a three-dimensional field space (l, l*, sigma). The claimed first-order transition is identified through a global minimum with <l> != <l*> (Fig. 1, right). Such a minimum lies outside the physical subspace, so the barrier between the two vacua in sigma may be an artifact of the extra unphysical direction. The paper's own check that the NJL model without the Polyakov loop has no first-order transition shows that the barrier is generated by the Polyakov-loop sector—exactly the sector treated with this ansatz. If the constrained problem has no true barrier, the first-order transition, the phase diagram, and every gravitational-wave number in Section V collapse. The Z_sigma positivity patch in Section IV.A is a second ad hoc element, but it affects the tunneling action and hence beta/H and alpha, not the existence of the transition; the l/l* ansatz is the more load-bearing assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript proposes that a large dark-quark chemical potential, mu ~ T, can change the chiral transition in a three-flavor dark-QCD from a crossover into a first-order phase transition. The authors work in the PNJL model in the chiral limit with a vanishing KMT term, compute the effective potential V(l,l*,sigma,T,mu), plot the phase diagram for several values of G_S, and then reduce the tunneling problem to a single field sigma to compute the nucleation temperature, alpha, and beta/H. These parameters are fed into standard sound-wave and turbulence formulas to predict gravitational-wave spectra. The paper also sketches an Affleck-Dine mechanism for generating the dark-quark asymmetry. The central claim is that for 1 GeV lesssim T_c lesssim 100 GeV the gravitational-wave signal from the dark chiral phase transition may reach the sensitivity of BBO.","tokens_in":28472,"tokens_out":9046,"duration_ms":90794,"significance":"If the first-order transition is real, this is a plausible new route to observable gravitational waves from a dark chiral transition and usefully complements existing studies of dark Yang-Mills deconfinement transitions. The paper is transparent about its free parameters, checks the known mu=0 crossover, and points to a GitHub implementation of the tunneling solver. Its quantitative gravitational-wave predictions are useful order-of-magnitude targets. At present, however, the existence and location of the first-order region rest on a technically uncontrolled treatment of the Polyakov-loop degrees of freedom, so the significance is conditional on the constrained calculation requested below.","major_comments":[{"comment":"The paper explicitly replaces the physical constraint l* = conj(l) by treating l and l* as two independent real fields, and Fig. 1 (right panel) identifies the global minimum with <l> different from <l*>. This minimum lies outside the physical configuration space. This is not a harmless parametrization: the NJL-only check at the end of Section III.B.1 shows that without the Polyakov-loop coupling there is no first-order transition, so the barrier is generated by the Polyakov-loop sector, exactly the sector whose treatment is relaxed. Please repeat the phase-diagram calculation with l and l* constrained to be complex conjugates, minimizing over the real and imaginary parts of l, and report whether the first-order region in Fig. 2 survives. If the complex effective potential makes this minimization ill-defined, state explicitly which physical approximation the independent-field extremum represents. Every gravitational-wave prediction in Sections IV and V inherits this uncertainty.","section":"Section III.B.1, before Eq. (3.14)"},{"comment":"The paper notes that Z(sigma) is negative in some region and then forces positive definiteness by replacing dI/dq^2 with dI/dq^2 + integral d^3k/(2 pi)^3 1/(8 E_k^5), citing Ref. [14]. This is an ad hoc subtraction, not derived in the present manuscript, and it enters the bounce equation (4.2), and hence the values of T_n, alpha, and beta/H in Table I. Please either derive this term from a systematic renormalization condition or quantify the sensitivity of the gravitational-wave spectra to the subtraction, including the option of no subtraction.","section":"Section IV.A, after Eq. (4.19)"},{"comment":"The note states that the altered tunneling-potential method gives S3/T approximately 0.85 times the exact numerical solution, while the main text and Section VI describe the two approaches as showing good consistency. A 15% difference in the action is not negligible for the computed nucleation temperature, which is set by S_E/T_n ~ 140, or for beta/H, which is obtained from the temperature derivative of the action. Please reconcile the two methods or quote the resulting systematic uncertainty in the phase-transition parameters and in the final gravitational-wave amplitudes.","section":"Section IV.B.1, note 4"},{"comment":"The authors acknowledge that different choices of the Polyakov-loop potential may give different results, but no test with an alternative form is provided. Since the first-order transition is generated by the interaction between the chiral condensate and the Polyakov-loop sector, this model uncertainty is directly load-bearing for the phase diagram. Please add a comparison with at least one alternative Polyakov-loop potential, for example the logarithmic form or the quasigluon-based model of Ref. [50], to show that the first-order region and the qualitative gravitational-wave conclusions are robust.","section":"Section III.A.2, after Eq. (3.11)"}],"minor_comments":[{"comment":"The text refers to the phase diagram as the right panel of Fig. 1, but the phase diagram appears to be in Fig. 2; the figure cross-references should be corrected.","section":"Section III.B.1"},{"comment":"There are numerous typos, including 'meta stable', 'psuedo', 'unphysic', 'zhelare', 'Afflect-Dine', and 'bonce' near Eq. (4.2); a careful proofread is needed.","section":"Throughout"},{"comment":"The notation ebeta could be misread as an exponential of beta; consider introducing an explicit symbol such as beta/H at first use and using it consistently.","section":"Eq. (4.20) and Table I"},{"comment":"The conversion from dark-quark number density to mu/T uses the degenerate-fermion formula although the relevant regime has mu/T ~ O(1); the paper should state explicitly that this is only an order-of-magnitude estimate.","section":"Section II, Eq. (2.11)"},{"comment":"The claimed check that the single-field reduction of the tunneling problem is accurate to one percent relies on the unpublished Ref. [53]; please include the comparison in the paper or cite a public version.","section":"Section IV, after Eq. (4.1)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a toy-model phenomenological study; its novelty relative to Refs. [31,32,39] is the transposition of the lepton-asymmetry mechanism to the dark sector, where no BBN bound applies. If the constrained Polyakov-loop calculation confirms the first-order transition, I expect the paper can become a solid contribution. I would ask the editor to insist on that constrained calculation before acceptance, because as it stands the central claim is not established. The authors are otherwise transparent about their assumptions and provide useful numerical checks."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Let me cut to the main point. This is a serious and mostly clean study of a dark chiral phase transition driven by a large quark chemical potential, with a gravitational-wave forecast. The new content is not the mechanism—large chemical potentials making chiral transitions first-order in PNJL models is established visible-QCD lore (refs 21–29)—but the application to a dark sector free of baryon-asymmetry constraints, plus the demonstration that mu can dramatically improve GW detectability. The AD toy model for generating the asymmetry is a nice addition, even if it is clearly schematic.\n\nWhat the paper does well: the numerical setup reproduces the known crossover at mu=0, the phase diagram is internally consistent, and the authors are honest about several limitations, including the comment about l and l*. The GW analysis is standard, and the observation that increasing mu at fixed T_c raises alpha and lowers beta/H is coherent and potentially interesting.\n\nNow the soft spots, in order of weight. First, the load-bearing issue is the treatment of l and l* as independent real fields (Sec. III.B.1). The authors acknowledge the physical constraint but then enlarge the field space to three real dimensions. The first-order transition is identified at a minimum with <l> != <l*>, which lies outside the physical subspace. The stress test is right: the NJL check without the Polyakov loop shows no FOPT, so the barrier is coming from the sector that is being treated with this ansatz. If the constrained problem removes the barrier, the central claim collapses. This needs a constrained re-computation. Second, the Z_sigma positive-definiteness patch (Sec. IV.A) is ad hoc and affects the tunneling action, so the quantitative alpha and beta/H numbers carry model uncertainty. Third, the abstract's claim that mu increases latent heat lacks the fixed-T_c caveat; with fixed T_0, latent heat actually decreases. Finally, no PNJL code or phase-diagram data are shipped, which limits independent checking—the linked VacuumTunneling package is a separate tool.\n\nOverall: the paper is a plausible exploratory study, not a settled prediction. It deserves a serious referee, because the parameter space is new and the technical concern is addressable rather than obviously fatal. I would recommend conditional acceptance with the constrained l/l* computation and code/data release as requirements.","headline":"A plausible PNJL study of a dark chiral transition driven by large chemical potential, but the quantitative GW signal rests on an under-validated Polyakov-loop ansatz.","tokens_in":29032,"tokens_out":3810,"would_cite":true,"duration_ms":35892,"reading_group":"yes","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 claims that a large dark-quark chemical potential turns the dark-QCD chiral transition first-order and could make its gravitational waves detectable by BBO.","keywords":["dark QCD","chiral phase transition","chemical potential","Polyakov loop","gravitational waves","first-order phase transition","Affleck-Dine mechanism","PNJL model"],"falsifier":"Run a sign-problem-free lattice simulation of three-flavor SU(3) dark QCD at imaginary chemical potential (or with a conserved isospin charge) and analytically continue to real $\\mu$; if no first-order chiral transition appears for $\\mu/T \\approx 1.4$--$1.8$ at the couplings studied, the central claim fails. More directly, re-minimize the same PNJL potential on the constrained subspace $l^* = l^\\dagger$ and check whether a barrier between the two vacua survives.","tokens_in":27903,"feed_emoji":"🌊","tokens_out":5093,"duration_ms":45589,"temperature":0.7,"pith_summary":"This paper claims that a hidden QCD-like sector—dark-QCD—with a large chemical potential for its quarks undergoes a first-order chiral phase transition in the early universe, even without the instanton-induced KMT interaction that earlier work required. The authors show that the chemical potential makes the Polyakov loop complex, tilts the partial-confinement interaction, and creates a barrier between the chirally symmetric and broken vacua. The first-order transition is the chiral one; the confinement transition itself remains a cross-over. Because a large chemical potential also lengthens the transition's duration and boosts its latent-heat release, the gravitational-wave signal from a transition at temperatures between 1 GeV and 100 GeV could fall within the sensitivity of a future BBO-class space detector.","feed_headline":"Dark quark density can turn dark-QCD transition first-order","feed_subtitle":"If right, the gravitational waves from a 1-100 GeV transition could fall in BBO's sensitivity band.","key_machinery":"The central object is the PNJL effective potential $V_{\\mathrm{PNJL}}(l, l^*, \\sigma, T, \\mu)$, which couples the chiral condensate $\\sigma$ to the traced Polyakov loop $l$ and its conjugate $l^*$. Its thermal part is written with a generalized chemical potential shifted by the background gauge field, producing the complex-valued loop that breaks the degeneracy between quark and antiquark contributions. The first-order nature of the transition comes from a fine balance between the zero-temperature and finite-temperature pieces of the potential near $T_c$, a balance that a large $\\mu$ weakens, thereby lowering $\\beta/H$ and increasing the released latent heat. The bounce calculation is modified by the field-dependent renormalization factor $Z_\\sigma$, which the authors derive at one loop in the Polyakov-loop background and then include in a modified bounce equation for $\\sigma$.","core_discovery":"Working in the three-flavor PNJL model in the chiral limit with vanishing KMT term, the authors find that a dark-quark chemical potential $\\mu$ of order of the temperature turns the dark chiral phase transition into a first-order one. The traced Polyakov loop $l$ and its conjugate pick up distinct vacuum values ($l \\neq l^*$), so the transition connects a deconfined, chirally symmetric vacuum to a partially confined, chirally broken vacuum through a barrier. In the phase diagram in the $T$-$\\mu$ plane, first-order transitions occur above a critical chemical potential that depends on the four-quark coupling $G_S$, with the region of first-order behavior shrinking as $G_S$ decreases and disappearing below $G_S\\Lambda^2 \\approx 1.8$. The authors further compute the tunneling action for the composite field $\\sigma$, including its loop-generated wave-function renormalization $Z_\\sigma$ in the Polyakov-loop background, and use the resulting phase-transition parameters to predict the gravitational-wave spectrum, concluding that for $T_c$ between 1 and 100 GeV the signal may reach BBO.","pith_inferences":["If the sign-problem artifact is not the real reason for $l \\neq l^*$, the quantitative gravitational-wave predictions (peak frequency and amplitude) could shift, though the qualitative mechanism of $\\mu$ enhancing the first-order character may still survive in a constrained treatment.","The same mechanism could be expected to work for larger color number $N$, where the latent heat release grows, potentially giving even stronger signals; the paper notes this but does not compute it.","A null detection by BBO in the $10^{-3}$--$1$ Hz band would translate into an upper bound on the combination of dark-quark chemical potential and transition temperature, provided the dark sector thermalizes with the visible one.","The role of $\\mu$ in relaxing the fine-tuned cancellation between $V_{\\mathrm{zero}}$ and $V_T$ suggests that other strongly coupled dark sectors with conserved charges may generically host stronger first-order transitions than previously estimated."],"forward_implications":["A dark-QCD sector with a large quark chemical potential can have a first-order chiral phase transition without relying on a large KMT instanton term.","For transition temperatures $T_c$ between 1 GeV and 100 GeV, the produced gravitational waves peak in the intermediate-frequency band and may be detectable by BBO.","Increasing $\\mu$ prolongs the duration of the phase transition (decreases $\\beta/H$) and raises the released latent heat for fixed $T_c$, both of which raise the gravitational-wave amplitude.","The confinement-deconfinement transition in this setup is a cross-over; the first-order transition is purely chiral, so a detected signal would point to the chiral sector rather than confinement.","The toy Affleck-Dine model shows that a large dark-quark asymmetry can be generated in the early universe, making the scenario physically realizable."],"supporting_citations":[{"why":"Supplies the PNJL model with Polyakov loop and chemical potential whose effective potential is used to plot the phase diagram.","marker":"[25]"},{"why":"Establishes the chiral effective model with the Polyakov loop that the dark-QCD calculation extends.","marker":"[40]"},{"why":"Provides the Polyakov-loop modeling framework, including the treatment of l and l* as independent real fields and the fermion-sign-problem rationale.","marker":"[44]"},{"why":"Earlier finite-density PNJL study that treats l and l* independently, the precedent the paper follows for the three-dimensional field space.","marker":"[29]"},{"why":"Dark chiral phase transition with KMT term; the paper shows the same first-order behaviour can be achieved by chemical potential instead.","marker":"[13]"},{"why":"Gives the wave-function renormalization factor calculation in a Polyakov-loop background that the authors adapt to include μ.","marker":"[14]"},{"why":"Provides the sound-wave gravitational wave spectrum formula used to convert phase-transition parameters into a predicted signal.","marker":"[67]"},{"why":"Defines the BBO detector sensitivity curve against which the predicted spectra are compared.","marker":"[78]"}],"fun_headline_variants":["Dark quark chemical potential flips chiral transition to first-order","Heavy dark quark density triggers first-order chiral transition","Gravitational waves from dark chiral transition may hit BBO","Dark QCD at high mu yields first-order transition and GW signal","Chemical potential makes dark chiral transition first-order"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation minimizes the effective potential with $l$ and $l^*$ treated as two independent real fields, even though they should be complex conjugates; if the true constrained problem has only conjugate minima, the barrier that makes the transition first-order could disappear.","fun_headline_variants_meta":{"raw":{"variants":["Dark quark chemical potential flips chiral transition to first-order","Heavy dark quark density triggers first-order chiral transition","Gravitational waves from dark chiral transition may hit BBO","Dark QCD at high mu yields first-order transition and GW signal","Chemical potential makes dark chiral transition first-order"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000195,"raw_usage":{"total_tokens":1383,"prompt_tokens":994,"completion_tokens":389,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":610,"completion_tokens_details":{"reasoning_tokens":309}},"tokens_in":610,"tokens_out":389,"duration_ms":3934,"temperature":1.0,"reasoning_tokens":309,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T14:29:44.816321+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a sign-problem-free lattice simulation of three-flavor SU(3) dark QCD at imaginary chemical potential (or with a conserved isospin charge) and analytically continue to real $\\mu$; if no first-order chiral transition appears for $\\mu/T \\approx 1.4$--$1.8$ at the couplings studied, the central claim fails. More directly, re-minimize the same PNJL potential on the constrained subspace $l^* = l^\\dagger$ and check whether a barrier between the two vacua survives.","supporting_citations":[{"cited_title":"Gravitational waves from composite dark sectors","cited_arxiv_id":null,"evidence_quote":"Supplies the PNJL model with Polyakov loop and chemical potential whose effective potential is used to plot the phase diagram."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Polyakov-loop modeling framework, including the treatment of l and l* as independent real fields and the fermion-sign-problem rationale."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier finite-density PNJL study that treats l and l* independently, the precedent the paper follows for the three-dimensional field space."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the sound-wave gravitational wave spectrum formula used to convert phase-transition parameters into a predicted signal."},{"cited_title":"The stochastic gravitational wave back- ground from turbulence and magnetic fields generated by a first-order phase transition.JCAP, 12:024, 2009","cited_arxiv_id":null,"evidence_quote":"Defines the BBO detector sensitivity curve against which the predicted spectra are compared."}],"review_version":1}