{"id":"c9ab8477-aad3-4291-b411-329c082c8e80","arxiv_id":"2608.12706","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A dark QCD crossover around 5.5 times the SM QCD scale imprints a frequency-shifted distortion on scalar induced gravitational wave spectra computed from rescaled effective degrees of freedom.","lead":"This paper computes how gravitational waves induced by primordial density fluctuations would be altered if a hidden 'dark QCD' sector confines at about 5.5 times the ordinary QCD scale. It finds a shifted spectral distortion that could, in principle, be searched for in pulsar timing array data.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The dark QCD feature location assumes the twin sector shares the SM temperature; if the twin QCD bath is colder, the claimed k/k_QCD~5.5 feature shifts and can overlap the SM feature, so the signature location is not fixed by the confinement scale alone.","rationale":"The reader's weakest_assumption listed two related issues: the rescaled dark QCD equation of state and the implicit common-temperature assumption. I agree that both deserve attention, but I identify the common-temperature assumption as the single most load-bearing concern because it controls the location of the predicted feature, which is the paper's headline signature. The rescaling uncertainty alters the shape and amplitude of the distortion but does not change the existence of a crossover-induced feature; a different crossover shape would still produce a dip in c_s^2, just with modified depth and width. In contrast, if the dark sector is colder than the SM, the feature location in physical frequency changes by the factor 1/r, potentially moving it on top of the SM QCD feature or far away from it. The paper's own benchmark motivation, asymmetric twin baryon dark matter, typically requires a suppressed twin radiation density, which often implies T_d < T, so the r=1 choice is not self-evident. The paper explicitly says it does not model the twin leptonic sector or dark photon dynamics (Section 2.2), and therefore the temperature ratio is left unspecified. Because the reader already issued a CONDITIONAL verdict, this concern does not change the verdict; it sharpens the condition that must be stated and checked before the '5.5' prediction is presented as robust. A two-temperature recomputation is straightforward and would settle whether the feature location is genuinely tied to the confinement scale ratio or to the unknown temperature ratio.","tokens_in":24360,"tokens_out":7870,"duration_ms":83867,"concrete_test":"Repeat the calculation for a dark sector with temperature ratio r = T_d/T in {0.2, 0.5, 1.0}: replace g*_d(T) by g*_d(rT) when constructing the total g* and g*_s, recompute w(eta) and c_s^2(eta), and then recompute the SIGW spectra for the benchmark k_* values in Table 2. If the trough position and the relative-difference peak in Fig. 5 shift by more than O(1) in k/k_QCD, the common-temperature assumption is load-bearing; if the feature remains near k/k_QCD ~ 5.5 for all r, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim is that the dark QCD imprint appears at k/k_QCD ~ 5.5 (Conclusion, Section 5). This follows only if the dark QCD bath temperature equals the SM bath temperature, an assumption the paper makes implicitly by constructing g*(T) with the dark sector contribution evaluated at the same T (Section 2.2). In the cited asymmetric twin baryon scenarios, the twin sector need not share the SM temperature: constraints on Delta Neff typically require the twin electromagnetic/leptonic sector to be colder or diluted, and the twin QCD bath may decouple with a temperature ratio r = T_d/T != 1. If T_d = r T, the dark crossover occurs when the SM temperature is T ~ (5.5/r) T_QCD, shifting the feature to k/k_QCD ~ 5.5/r (up to g* corrections). For r ~ 1/5.5, the dark feature coincides with the SM feature and the claimed frequency-shifted distortion disappears. The paper gives no argument that r=1 is enforced in the asymmetric twin baryon benchmark; it explicitly omits the twin leptonic sector where the temperature ratio would be set. Because the headline signal is a frequency-shifted feature, this missing parameter controls the central prediction. This is not a disagreement with the general mechanism that a hidden crossover will imprint on SIGWs, but the specific '5.5' signature is conditional on r=1.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the energy-density spectrum of scalar induced gravitational waves (SIGWs) produced by a monochromatic primordial curvature perturbation, in a thermal history that includes not only the Standard Model QCD crossover but also a dark QCD sector with confinement scale Λ_dQCD ≃ 5.5 Λ_QCD, motivated by asymmetric twin baryon dark matter. The authors construct effective energy and entropy degrees of freedom for the SM plus a three-flavor or six-flavor dark QCD sector, using a temperature rescaling of SM hadron resonance gas, lattice, and perturbative QCD results. They then solve the first-order scalar perturbation equation and the second-order tensor equation numerically, with w(η) and c_s^2(η) derived from g_*(T) and g_*s(T). The main findings are that the realistic SM thermal history modifies the radiation-dominated SIGW spectrum in the trough region, and that adding the dark QCD sector produces an additional, frequency-shifted distortion when the scalar mode reenters the horizon near the dark confinement scale, with the dark feature located near k/k_QCD ~ 5.5.","tokens_in":24737,"tokens_out":15497,"duration_ms":166325,"significance":"If the assumptions hold, the paper provides a concrete template for using SIGWs to probe hidden confining sectors, extending the known SM QCD imprint to a well-motivated twin Higgs benchmark. Its strengths are the careful numerical treatment of the scalar transfer function with a time-dependent equation of state, the detailed reconstruction of g_* and g_*s from established SM QCD inputs, and the transparent comparison with the analytic radiation-dominated limit. The central quantitative claim, however, depends on an unstated and unargued assumption that the dark QCD bath has the same temperature as the SM bath, and on the assumption that the dark crossover is an exact rescaled copy of the SM one. These conditions control the location and shape of the predicted feature and need to be made explicit and tested before the result can be regarded as a robust 'characteristic signature.'","major_comments":[{"comment":"The dark-sector g_* and g_*s are constructed by evaluating SM QCD results at temperatures rescaled by a factor of 5.5 relative to the SM bath temperature T. This implicitly assumes T_d = T, i.e., that the dark QCD and SM baths share a common temperature. The paper gives no argument for this in the asymmetric twin baryon benchmark: it explicitly omits the twin leptonic sector and the dark photon, and it does not specify the decoupling or reheating dynamics that would set the temperature ratio. If T_d = r T, the dark crossover occurs when the SM temperature is T ≈ (5.5/r) T_QCD, shifting the feature to k/k_QCD ≈ 5.5/r up to g_* corrections. For r ≈ 1/5.5 the dark feature would coincide with the SM QCD feature, and the claimed frequency-shifted distortion would disappear. Because the load-bearing quantitative claim in the Conclusion is the location k/k_QCD ≈ 5.5, this missing parameter controls the central prediction. I request either an explicit model-based derivation of r = 1 or a scan over r, including the r < 1 regionmotivated by asymmetric reheating or entropy injection, with the resulting shift of the feature shown.","section":"§2.2 (SM + dQCD3 construction) and §5 (Conclusion)"},{"comment":"The dark QCD equation of state is modeled by applying a single temperature rescaling of 5.5 to the SM hadron resonance gas, lattice trace anomaly, and perturbative QCD expressions. The spectral distortion computed in §4.2 is therefore conditional on the dark crossover being an exact rescaled copy of the SM crossover. A dark gauge theory with a first-order transition, a different N_f behavior, or a different relation between quark masses and the confinement scale would give a different distortion shape and amplitude. While the paper does state this as an approximation, the abstract and conclusion present the result as 'a characteristic signature of the dark QCD sector.' I recommend adding an explicit limitation and, if feasible, a simple robustness check (for example, a sharper crossover or a different N_f) to show which features of the spectrum are generic to any hidden confining sector and which are specific to the assumed benchmark.","section":"§2.2, §4.2 and Abstract/Conclusion"}],"minor_comments":[{"comment":"The caption of Table 2 says the cases are 'used to study the SM QCD crossover effect'; this should read 'dark QCD crossover effect' since the table describes SM + dQCD3 benchmarks.","section":"Table 2 caption"},{"comment":"Cases d–f select the peak scale k_* by hand so that it lies on the dark crossover region. A figure scanning k_* over a wider range would make the claim that a distinct feature appears near k/k_QCD ~ 5.5 visually supported rather than inherited from the input benchmark Λ_dQCD = 5.5 Λ_QCD and the chosen k_* values.","section":"§4.2 and Fig. 5"},{"comment":"The dark hadron masses in the HRG are rescaled by 5.5 while the dark quark masses are taken to be 3 times the SM quark masses; for pseudo-Goldstone bosons these two scalings are not mutually consistent, and the text should clarify that the 5.5 rescaling of the full hadron spectrum is a simplification.","section":"§2.2, HRG rescaling"},{"comment":"The lower panel of Fig. 5 plots an absolute relative difference whose denominator can be very small near the trough; the large values in that panel should be discussed so that the reader does not interpret them as a large absolute change in the spectral amplitude.","section":"§4.2, lower panel of Fig. 5"}],"recommendation":"major_revision","confidential_remarks":"The paper is a competent extension of the Abe–Tada–Ueda SM QCD SIGW analysis to a twin QCD sector, and the numerical implementation appears sound. The main issue is not the SIGW formalism but the model input: the central prediction of a feature at k/k_QCD ~ 5.5 assumes a common SM/dark bath temperature without argument, and the dark-sector EoS is a rescaled copy of the SM QCD EoS. If the authors add an explicit treatment or scan of the temperature ratio r = T_d/T, soften the 'characteristic signature' claim accordingly, and clarify the HRG mass-scalings, I think the paper can be made acceptable for JCAP."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my take on 2608.12706. It's a competent extension of the known SIGW-through-QCD-crossover computation to a benchmark dark QCD sector. The genuinely new pieces are the effective g*(T) and g_s(T) curves for three- and six-flavor dark QCD with Lambda_d = 5.5 Lambda_QCD, and the demonstration that a dark crossover leaves a frequency-shifted dip in the monochromatic SIGW spectrum. The machinery is standard and the numerical work is described carefully; the appendix on the SM g* reconstruction is solid and useful.\n\nWhere I'd push back: the characteristic feature at k/k_QCD ~ 5.5 is not derived from dark QCD dynamics; it is built in by choosing Lambda_d/Lambda_QCD = 5.5 and scaling the SM EoS. That is fine as a benchmark, but the paper should be clearer that the 'signature' is conditional on that benchmark. More importantly, the calculation implicitly assumes the dark QCD bath shares the SM temperature T. The stress-test note is right: if T_d = r T with r < 1, the dark crossover occurs at a different SM temperature and the feature location shifts to k/k_QCD ~ 5.5/r; for r ~ 1/5.5 it overlaps the SM QCD feature entirely. The paper never states the equal-temperature assumption or discusses how Delta Neff constraints or twin-lepton dynamics fix r. Since the headline is a frequency shift, this is the softest spot in the paper. Also, the thermodynamic uncertainty band for g* is constructed but not propagated to the SIGW spectra, so we don't know whether the claimed spectral differences are robust within their own EoS uncertainty. No code is provided, which is a minor reproducibility annoyance.\n\nNone of this kills the paper. The mechanism — a hidden confining crossover imprints on SIGWs — is plausible and the proof-of-principle is legitimate. The authors just need to state the temperature assumption, show how the feature depends on the temperature ratio, and propagate their EoS uncertainties. A JCAP referee should ask for those changes, not reject the paper.\n\nI'd send it to a serious referee; for my own work I probably wouldn't cite it, and I'm not sure I'd schedule it in the reading group.","headline":"A careful but conditional SIGW probe of a dark QCD crossover; the headline frequency shift depends on an unstated equal-temperature assumption.","tokens_in":25251,"tokens_out":5013,"would_cite":false,"duration_ms":45714,"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":"A dark QCD crossover imprints a frequency-shifted distortion on scalar-induced gravitational waves, offering a way to probe hidden confining sectors.","keywords":["scalar induced gravitational waves","dark QCD","QCD crossover","twin Higgs","asymmetric twin baryon dark matter","equation of state","sound speed","pulsar timing arrays"],"falsifier":"A measurement of the SIGW spectrum that resolves the SM QCD trough near k/k_QCD ≈ 1 but shows no enhanced dip near k/k_QCD ≈ 5.5 would rule out the benchmark dark QCD signal; likewise, a lattice computation of a dark gauge theory with a first-order or differently shaped transition would change the predicted distortion and could be compared directly with the smooth-rescaling prediction.","tokens_in":1884,"feed_emoji":"🌊","tokens_out":2599,"duration_ms":75605,"temperature":0.7,"pith_summary":"This paper argues that scalar-induced gravitational waves (SIGWs) can act as a cosmological probe of a hidden, strongly coupled dark QCD sector. It constructs the effective energy and entropy degrees of freedom for the Standard Model plus a dark QCD sector with a confinement scale about 5.5 times that of SM QCD, builds the resulting equation of state and sound speed, and solves for the second-order gravitational waves sourced by a monochromatic primordial curvature perturbation. The central result is that the dark QCD crossover produces its own spectral distortion, shifted to higher frequency than the SM QCD imprint, so the two features would appear at k/k_QCD ≈ 1 and k/k_QCD ≈ 5.5 respectively. If real, this gives gravitational-wave observatories a way to detect strong dynamics that is otherwise almost invisible to direct experiments.","feed_headline":"Dark QCD shifts a gravitational-wave signature to higher frequencies","feed_subtitle":"The dark confinement feature should appear at about 5.5 times the QCD scale, a higher frequency than the visible imprint.","key_machinery":"The load-bearing object is the SIGW kernel built from the first-order scalar transfer function T_phi(k, eta), solved numerically with the time-dependent equation of state w(eta) and sound speed $c_s^{2}$(eta) derived from the effective degrees of freedom g_*(T) and g_{*s}(T). The dark sector input is obtained by a simple temperature rescaling by 5.5 of the SM hadron resonance gas, lattice trace anomaly, and perturbative QCD pressure, with quark masses tripled. The paper uses the diagnostic that the trough in the gravitational wave energy density spectrum sits at approximately $\\sqrt$(2) c_s(eta_cancel), where eta_cancel is fixed by k_* eta_cancel approximately equal to 4, which ties the spectral feature to the sound-speed dip of the crossover.","core_discovery":"The paper establishes that a dark QCD crossover — a smooth confinement transition in a hidden gauge sector at roughly 5.5 times the SM QCD scale — leaves a characteristic, frequency-shifted imprint on the scalar-induced gravitational wave spectrum. Treating the dark QCD thermodynamics as a temperature-rescaled copy of the SM QCD thermodynamics, with hadron masses scaled by 5.5 and quark masses by 3, the paper computes the equation of state parameter w(T) and the sound speed $c_s^{2}$(T), solves the first-order scalar perturbation equation with these time-dependent coefficients, and feeds the result into the second-order tensor source. For a monochromatic primordial power spectrum, the exact zero of the radiation-dominated spectrum is lifted to a finite trough and the resonant logarithmic singularity becomes a finite cusp. When the scalar mode reenters near the dark QCD crossover, the depth of that trough differs measurably from the SM-only prediction, with the dark feature sitting near k/k_QCD ≈ 5.5, at higher frequency than the SM feature near k/k_QCD ≈ 1.","pith_inferences":["The same rescaling logic could be turned around: detecting a shifted trough would measure the ratio of the dark confinement scale to the visible one, while its absence would push the dark confinement scale away from the 5.5 benchmark or require the dark sector to be colder than the SM bath.","If the dark transition is first order rather than a crossover, the evolution of w and c_s^2 would be sharper and the spectral distortion, especially the shape of the trough and cusp, would differ from the smooth-rescaling prediction; this makes the SIGW shape a potential discriminator of transition order.","Because the paper omits the twin leptonic sector and dark photon on the grounds of Delta N_eff constraints, a more complete thermal history including those states could shift g_* slightly and move the dark feature's location; quantifying that shift is a natural follow-up.","The method should transfer to other hidden confining sectors with different numbers of flavors and quark masses, so the same pipeline could map out a family of dark confinement signatures in the SIGW band."],"forward_implications":["The realistic SM thermal history changes the SIGW spectrum relative to the idealized radiation-dominated case: the exact zero becomes a finite trough and the logarithmic resonance becomes a smooth cusp.","A dark QCD crossover at 5.5 times the SM QCD scale adds a second, higher-frequency distortion; in units of the SM QCD scale, the SM feature appears near k/k_QCD ≈ 1 and the dark feature near k/k_QCD ≈ 5.5.","For monochromatic primordial spectra whose peak scale reenters near the dark crossover, the trough value of the gravitational wave energy density fraction is significantly larger than in the SM-only history, while for modes reentering well below the dark crossover the spectra nearly coincide.","The six-flavor dark QCD benchmark gives a spectrum qualitatively similar to the minimal three-flavor case in the temperature range most relevant for the dark crossover.","With suitable template fits, pulsar timing array data could in principle search for both the SM QCD and the shifted dark QCD features against smooth astrophysical backgrounds."],"supporting_citations":[{"why":"Supplies the reconstructed SM g_* and g_{*s} over a wide temperature range and the matching procedure (hadron resonance gas, lattice, perturbative QCD) that the dark sector reuses by rescaling.","marker":"[1]"},{"why":"Established that the SM QCD crossover imprints on SIGWs and provides the reference peak scale that the paper rescales by 5.5 to probe the dark crossover.","marker":"[5]"},{"why":"Provides the lattice QCD trace-anomaly fitting functions for 2+1+1 flavors that are temperature-rescaled to model the dark QCD crossovers.","marker":"[33]"},{"why":"Motivates asymmetric twin baryon dark matter with an enhanced dark confinement scale and the observed dark-matter to baryon coincidence.","marker":"[13]"},{"why":"Sets the benchmark dark proton mass of about 5.5 GeV, giving the dark confinement scale approximately 5.5 times the SM QCD scale.","marker":"[15]"},{"why":"Introduces the twin Higgs framework that motivates the existence of a dark QCD sector.","marker":"[16]"},{"why":"Gives the analytic radiation-dominated SIGW spectrum whose zero, resonant peak, and cutoff are the baseline that the crossover distorts.","marker":"[37]"},{"why":"Supplies the second-order SIGW formalism and scalar transfer function used in the numerical solution.","marker":"[36]"}],"fun_headline_variants":["Dark QCD leaves a higher-frequency mark on gravitational waves","Dark QCD shifts gravity-wave imprint to 5.5x QCD frequency","Dark QCD crossover leaves a higher-frequency ripple in gravity waves","Gravitational waves can probe dark QCD via a shifted frequency peak","Dark QCD's gravitational-wave imprint appears at 5.5x the QCD scale"],"cache_read_input_tokens":27264,"weakest_assumption_plain":"The calculation assumes the dark QCD crossover is a rescaled copy of the SM QCD crossover, with temperatures shifted by 5.5 and quark masses tripled, and that the dark sector bath shares the same temperature as the visible one.","fun_headline_variants_meta":{"raw":{"variants":["Dark QCD leaves a higher-frequency mark on gravitational waves","Dark QCD shifts gravity-wave imprint to 5.5x QCD frequency","Dark QCD crossover leaves a higher-frequency ripple in gravity waves","Gravitational waves can probe dark QCD via a shifted frequency peak","Dark QCD's gravitational-wave imprint appears at 5.5x the QCD scale"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000906,"raw_usage":{"total_tokens":3924,"prompt_tokens":1002,"completion_tokens":2922,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":618,"completion_tokens_details":{"reasoning_tokens":2825}},"tokens_in":618,"tokens_out":2922,"duration_ms":21934,"temperature":1.0,"reasoning_tokens":2825,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T04:08:19.587272+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A measurement of the SIGW spectrum that resolves the SM QCD trough near k/k_QCD ≈ 1 but shows no enhanced dip near k/k_QCD ≈ 5.5 would rule out the benchmark dark QCD signal; likewise, a lattice computation of a dark gauge theory with a first-order or differently shaped transition would change the predicted distortion and could be compared directly with the smooth-rescaling prediction.","supporting_citations":[],"review_version":1}