{"id":"ad632205-d9e7-43d6-8696-ffaff9221781","arxiv_id":"2507.12968","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A five-phase bouncing cosmology yields a broken power-law gravitational wave spectrum whose amplitude bound automatically keeps the bounce energy below the Planck scale.","lead":"This paper adds an extra early contraction phase to a bouncing cosmology and derives the gravitational wave spectrum it produces, a broken power law. The result suggests that any bouncing model passing current cosmic radiation bounds stays below the Planck energy scale, making it potentially testable in future detectors.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The ΔNeff bound does not force w1≤0; blue-tilted low-amplitude spectra satisfy it, so the paper's proof that all viable NMBC models are sub-Planckian is incomplete.","rationale":"The reader's weakest_assumption focused on the high-frequency branch equality and phase-0 leakage, but the reader's rationale explicitly noted that the argument forcing w1 ≤ 0 is logically unsound and that a proof for w1 > 0 is absent. That is the most load-bearing issue because the abstract's universal claim ('all NMBC models satisfying the bound') hinges on excluding w1 > 0. The paper does not provide that exclusion: an upper bound on a blue-tilted spectrum can always be satisfied by lowering the overall amplitude, so no tilt restriction follows. This is an internal logical gap rather than a disagreement with external consensus, and it directly affects the core conclusion. The concrete test would settle whether the missing w1 > 0 case actually permits trans-Planckian bounce scales; until that test is run, the central claim is unproven, though not necessarily false. The reader's conditional verdict already requires fixing this type of issue, so no change to the verdict is needed; the concern reinforces the conditional status.","tokens_in":15589,"tokens_out":11131,"duration_ms":119829,"concrete_test":"Numerically solve Eq. (19) against Eq. (24) for the w1 > 0 branch, using the highest constrained frequency f = 10^9 Hz (the rightmost pivot in Fig. 2), with C^{(1)}(ν1) from Eq. (20) and the stated cosmological constants. For w1 = 1/9, 1/3, 1, set Ω_GW h²(10^9 Hz) = 1.7×10⁻⁶ and solve for ρ_{s↓}^{1/4}. If any result exceeds 0.79 m_pl, the central claim is false and the verdict should be REJECT; if all results are below, the missing w1 > 0 proof can be supplied and the conditional acceptance can stand. Also verify that the w1 = 0 case reproduces 0.79 m_pl as a sanity check.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that every NMBC model obeying Ω_GW h²(f) < 1.7×10⁻⁶ (Eq. 24) has ρ_{s↓}^{1/4} < 0.79 m_pl. The paper's derivation passes through the assertion that the ΔNeff bound 'requires a red or scale-invariant high-frequency SGWB', i.e. n_T^{(0)}(f ≥ f⋆) ≤ 0, which with Eq. (21) gives w1 ≤ 0. This inference is a logical non-sequitur. Eq. (24) is an upper limit on the SGWB amplitude, not a condition on its tilt. For any w1 > 0 (blue tilt, n_T > 0), the spectrum Ω ∝ f^{n_T} Teq is maximized at the highest constrained frequency; by choosing the amplitude factor A (proportional to ρ_{s↓}^{(4−n_T)/4}, Eq. (19)) small enough, one satisfies Eq. (24) at all frequencies. No independent normalization fixes A, so models with w1 > 0 are not excluded. The subsequent 0.79 m_pl bound is derived only on the w1 ≤ 0 branch (Fig. 2), and the paper contains no equivalent calculation for w1 > 0. Consequently, the 'all NMBC models' statement in the abstract and summary is not established. This is a correctness risk in the central claim, not merely a matter of presentation.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper introduces the next-to-minimal bouncing cosmology (NMBC), a five-phase bouncing model obtained by adding an early contraction phase (phase 0) to the minimal bouncing cosmology (MBC). The authors extend their earlier matrix-propagation method to derive a closed-form, broken power-law stochastic gravitational-wave background (SGWB), with the high-frequency branch identical to the MBC result and the low-frequency branch modified by phase 0. They then claim that the current ΔNeff bound Ω_GW h²(f) < 1.7×10⁻⁶ forces w1 ≤ 0, and hence that all viable NMBC models have bounce energy scale ρ_{s↓}^{1/4} < 0.79 m_pl, automatically avoiding the trans-Planckian problem. Four illustrative NMBC spectra are presented as examples.","tokens_in":15902,"tokens_out":8005,"duration_ms":86920,"significance":"The paper provides a genuinely useful analytic tool for SGWB predictions in multi-phase bouncing cosmologies: the matrix method yields closed-form spectra, the two internal consistency checks in Eqs. (31)-(32) are sensible, and the four examples are explicit, falsifiable, and tied to concrete detector sensitivities. If the sub-Planckian claim were fully established, it would be an important result connecting gravitational-wave observatories to the bounce energy scale. However, the central inference from the ΔNeff bound to w1 ≤ 0 is not logically valid, so the headline claim is not currently supported; the strength of the paper lies in the analytic machinery rather than in the demonstrated universality of the sub-Planckian conclusion.","major_comments":[{"comment":"The inference that the ΔNeff upper limit Ω_GW h²(f) < 1.7×10⁻⁶ requires a red or scale-invariant high-frequency SGWB, and hence w1 ≤ 0, is a non-sequitur. An upper limit on the amplitude does not restrict the tilt: for any w1 > 0 (n_T > 0), the spectrum is blue, but the amplitude factor in Eq. (19) contains (ρ_{s↓}^{1/4})^{4-n_T}, so the amplitude can be made arbitrarily small by choosing ρ_{s↓}^{1/4} small enough. No independent normalization fixes this amplitude, and the bound at any single frequency (or even the integrated BBN bound) can be satisfied by a sufficiently low bounce scale. Consequently, the subsequent derivation of the 0.79 m_pl bound applies only to the w1 ≤ 0 branch shown in Fig. 2, not to all NMBC models. The abstract's claim that 'all NMBC models satisfying the current ΔNeff bound automatically avoid the trans-Planckian problem' is therefore not established. Please either extend the calculation to w1 > 0, demonstrating (if true) that the amplitude constraint alone still forces sub-Planckian ρ_{s↓}^{1/4}, or explicitly restrict the claim to the w1 ≤ 0 branch. A concrete test would be to evaluate Eq. (19) at f = 10⁻⁷ Hz with, say, w1 = 1/3 and choose ρ_{s↓}^{1/4} so that Ω_GW h² = 10⁻⁷; this directly contradicts the asserted exclusion of w1 > 0.","section":"Constraint on ρ_{s↓}^{1/4} from ΔNeff, Eqs. (24)-(25)"},{"comment":"The four illustrative examples are intended to validate the analytic framework, but the parameters given for Examples 2 and 3 are internally inconsistent. In the main text, Eq. (36) states ρ_{s↓}^{1/4} = 0.06 m_pl for Example 2, while SM Sec. VII gives 0.06×10⁻⁷ m_pl, a discrepancy of seven orders of magnitude. For Example 3, Eq. (37) gives w0 = 1.3×10⁵, while SM Sec. VII gives w0 = 4×10⁵/3; moreover, neither value is consistent with the stated low-frequency tilt n_T^{(0)}(f < f⋆) = 4−10⁻⁵ through Eq. (30) for the ν0 ≤ 1/2 branch, which would require w0 ≈ 1. These inconsistencies undermine the claim that the examples validate the analytic formulas, even though they are illustrative rather than central to the sub-Planckian argument.","section":"Cosmological Applications and Supplemental Material Sec. VII"}],"minor_comments":[{"comment":"The initial Bunch-Davies vacuum (A0, B0) for phase 0 is written with exp(−i(ν̃1 π/2 + π/4)), but the phase should presumably involve ν̃0, not ν̃1.","section":"Supplemental Material, Eq. (S25)"},{"comment":"Equation (24) is stated for f ≳ 10⁻⁷ Hz, but the text then evaluates it at f = 7.75×10⁻¹⁷ Hz. Please clarify which experimental bound is being used at that frequency and whether the ΔNeff bound actually applies there.","section":"Constraint on ρ_{s↓}^{1/4}, after Eq. (24)"},{"comment":"The sentence 'which is belong to the branch ν1 > 1/2' is grammatically incorrect; please rephrase.","section":"Constraint on ρ_{s↓}^{1/4}, paragraph after Eq. (25)"},{"comment":"The reduction to MBC is verified only for the full kernel in the SM, while the simplified long-phase-1 formulas in Eqs. (28)-(30) do not apply in the η0↓ = ηs↓ limit. Please state this explicitly in the main text so that readers do not attempt to use the simplified formulas for the consistency check.","section":"Relation between MBC and NMBC, Eqs. (31)-(32)"}],"recommendation":"major_revision","confidential_remarks":"The paper fits the journal's scope as an analytic cosmology/GW paper. The reliance on the author's prior MBC results (Refs. [7,49]) is acceptable given that those are derivations, but the present manuscript does not reproduce the core kernel derivation and the central claim rests on a misreading of the ΔNeff bound. The parameter discrepancies in the examples should also be fixed before resubmission."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read it. The genuinely new result is the closed-form broken power-law SGWB for the five-phase NMBC, obtained by extending the matrix method from your MBC paper. The algebra is clean, the low-frequency kernel simplifies to a power law, and the two consistency limits (ν0=ν1 and η0↓=ηs↓) check out. The pivot matching factor ε is a good internal test. The four examples are explicitly illustrative, not fits, so I'm not holding that against the paper.\n\nThe problem is the headline claim. The abstract says all NMBC models satisfying the ΔNeff bound automatically avoid the trans-Planckian problem, with ρ_s↓^{1/4}<0.79 m_pl. That does not follow. Eq. (24) bounds the amplitude of Ω_GW h² at the relevant frequencies; it says nothing about the tilt. A blue-tilted spectrum (w1>0, n_T>0) can be made to satisfy Eq. (24) by choosing a small overall amplitude A, which is proportional to ρ_s↓^{(4-n_T)/4}. The paper's step from the bound to w1≤0 (Eq. (25)) is a non-sequitur. The 0.79 m_pl intersection is computed only on the w1≤0 branch (Fig. 2). So the 'all NMBC models' statement is not established. That is a load-bearing flaw, not a presentation issue.\n\nThere are also a couple of SM typos a referee will trip on: Example 2 lists ρ_s↓^{1/4}=0.06×10^{-7} m_pl while the main text says 0.06 m_pl, and Eq. (S25) uses ν1 in the phase-0 vacuum phase where it should be ν0. Minor, but sloppy.\n\nThe high-frequency branch assumes phase 0 doesn't affect modes exiting during phase 1; that's plausible, and the paper is honest about the pivot discontinuity, but the trans-Planckian conclusion leans on that assumption plus the normalization from Ref. [49]. The self-citation doesn't bother me because that prior result is a parameter-free derivation.\n\nWho this is for: anyone constructing SGWB templates from multi-phase bouncing cosmologies. The broken power-law formula and the matrix method will be useful. But the abstract should be revised to claim sub-Planckian for the red/scale-invariant branch, or the blue branch needs a real proof.\n\nRecommendation: send it to peer review. The framework is worth a referee's time. The overclaim is fixable, and the derivations deserve scrutiny from someone who can check the SM kernel.","headline":"Solid analytic SGWB template for five-phase bounces, but the abstract's 'all NMBC models' overclaims: ΔNeff bounds amplitude, not tilt, so the sub-Planckian conclusion is unproven for blue-tilted spectra.","tokens_in":16429,"tokens_out":3557,"would_cite":true,"duration_ms":37067,"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":"This paper derives a closed-form broken power-law gravitational-wave background for a five-phase bouncing cosmology and shows that the ΔNeff bound automatically forces the bounce energy scale below the Planck mass.","keywords":["bouncing cosmology","stochastic gravitational-wave background","primordial gravitational waves","trans-Planckian problem","broken power law","matrix propagation method","ΔNeff bound","pivot frequency"],"falsifier":"Numerically integrate the tensor perturbation equation across the five phases for the paper's four example parameter sets without imposing $k\\eta_{s\\downarrow}\\ll k\\eta_{0\\downarrow}\\ll1$ or the $k\\eta_{0\\downarrow}\\ll1$ approximation, and compare the exact spectrum with Eqs. (22) and (28); any material departure on the $f\\ge f_\\star$ branch would shift the $0.79\\,m_{\\rm pl}$ intersection and weaken the automatic sub-Planckian conclusion.","tokens_in":15341,"feed_emoji":"🌌","tokens_out":19990,"duration_ms":182751,"temperature":0.7,"pith_summary":"Bouncing cosmologies replace the initial singularity with a finite minimum scale, but the minimal four-phase version predicts a featureless power-law gravitational-wave background that is hard to observe and often requires super-Planckian bounce energies. This paper adds one early contraction phase, making a next-to-minimal bouncing cosmology (NMBC) whose gravitational-wave spectrum is a broken power law: high frequencies match the minimal model, while low frequencies carry a new tilt set by the extra phase. Using a matrix-propagation method built from an inequality algebra, the author derives a closed-form expression for the spectrum and shows that the CMB/BBN bound $\\Omega_{\\rm GW}h^2<1.7\\times10^{-6}$ forces the bounce energy scale below $\\rho_{s\\downarrow}^{1/4}=0.79\\,m_{\\rm pl}$. The result matters because it makes bouncing cosmologies with no trans-Planckian physics potentially detectable across CMB, pulsar-timing, interferometer, and laboratory gravitational-wave experiments.","feed_headline":"Bounce models passing ΔNeff stay below Planck scale","feed_subtitle":"An added phase yields a broken power-law gravitational-wave signal testable across existing and planned detectors.","key_machinery":"The carrying object is the PGW propagation kernel $N^{(0)}_{22}(\\{\\tilde\\nu_i\\},\\{\\eta_{i\\downarrow/\\uparrow}\\})$, the (2,2) entry of $X^{(0)\\dagger}X^{(0)}$, where $X^{(0)}$ is a product of phase transformation matrices $T_i$ and boundary-matching matrices $M_{i\\uparrow/\\downarrow}$ that propagate Bunch-Davies vacuum amplitudes across the five phases. An inequality-based algebra—$\\chi_i=0$ or $\\bar\\chi_i=0$ depending on whether $\\nu_i>1/2$ or $\\nu_i\\le1/2$—collapses the matching matrices into two branches and lets the kernel be evaluated in closed form. The extra phase 0 enters only for $f<f_\\star$ through the factor $T_1^{-1}M_{0\\downarrow}T_0$, producing the broken power law; the pivot frequency $f_\\star=(2\\pi a_0\\eta_{0\\downarrow})^{-1}$ is set by the conformal time at the phase-0→1 transition.","core_discovery":"The central claim is that adding a fifth, early contraction phase (phase 0) to the minimal bouncing cosmology produces a stochastic gravitational-wave background with an analytically computed broken power law, with pivot frequency $f_\\star=(2\\pi a_0\\eta_{0\\downarrow})^{-1}$. The amplitude matrix $X^{(0)}$ equals the MBC matrix for $f>f_\\star$ and acquires a phase-0 factor $T_1^{-1}M_{0\\downarrow}T_0$ for $f<f_\\star$; in the long-phase-1 limit the low-frequency spectrum reduces to a closed form with tilt $n_T^{(0)}=3\\mp2\\tilde\\nu_0$ and coefficient $C^{(0)}$, while the high-frequency branch is exactly the earlier minimal-bouncing result. Imposing the $\\Delta N_{\\rm eff}$ bound $\\Omega_{\\rm GW}h^2\\lesssim1.7\\times10^{-6}$ at the five pivot frequencies $7.75\\times10^{-17}$, $10^{-8}$, $30$, $10^7$, and $10^9$ Hz forces $w_1\\le0$; at $w_1=0$ all constraint curves intersect at $\\rho_{s\\downarrow}^{1/4}=0.79\\,m_{\\rm pl}$, so every viable NMBC model bounces below the Planck scale. Four worked parameter sets reproduce features such as the PTA best fit and lie within projected sensitivities of CMB, space, ground, and laboratory detectors.","pith_inferences":["If the matrix algebra generalizes as the paper suggests, closed-form spectra could be written for any multi-phase history—kination, reheating, or dark phase transitions—turning parameter scans into analytic exercises rather than numerical grids.","The case of zero high-frequency tilt acts as a selection rule for model builders: targeting a red or scale-invariant high-frequency tilt keeps the bounce sub-Planckian, while blue tilts push the required bounce energy past the Planck scale.","A future broken power-law detection would effectively measure two pre-bounce conformal-time scales and two equation-of-state parameters, turning the framework into an observational probe of the contraction phase.","The small discontinuity at the pivot frequency is a direct diagnostic of the approximation used; precise numerical spectra near the pivot could test whether phase-0 leakage is genuinely negligible for the sub-Planckian bound."],"forward_implications":["Every NMBC model satisfying the current ΔNeff bound has a red or scale-invariant high-frequency tilt, so the bounce energy scale is automatically sub-Planckian (≤ 0.79 Planck mass), removing the trans-Planckian problem without fine-tuning.","The high-frequency branch inherits the full minimal bouncing cosmology spectrum and its constraints, so all existing MBC bounds continue to apply unchanged to those modes.","The broken power law with pivot frequency set by the phase-0 transition opens broad observational windows; the four example spectra sit within projected CMB, PTA, LVK, and superconducting circuit/cavity sensitivities.","The closed-form kernel makes parameter inversion direct: specified low- and high-frequency tilts, pivot frequency, and amplitude uniquely determine the equation-of-state parameters, the bounce scale, and the cutoff frequency.","The matrix-propagation algebra extends to other multi-phase early-universe histories, giving closed-form SGWB predictions that complement numerical approaches."],"supporting_citations":[{"why":"Supplies the matrix-representation method, the phase transformation and boundary-matching matrices, and the normalization C^(1)(ν1) that the NMBC high-frequency branch inherits.","marker":"[49]"},{"why":"Derives the minimal bouncing cosmology SGWB power law to which the NMBC spectrum reduces at f ≥ f⋆.","marker":"[7]"},{"why":"Sets the CMB extra-radiation bound Ω_GW h² < 1.7×10⁻⁶ used to derive w1 ≤ 0 and the 0.79 m_pl intersection.","marker":"[50]"},{"why":"Together with [50], provides the BBN/CMB ΔNeff constraint quoted in Eq. (24).","marker":"[51]"},{"why":"Provides current and projected sensitivity curves (CMB, PTA, LVK) against which the four example spectra are tested.","marker":"[8]"},{"why":"Supplies laboratory detector (superconducting circuit and cavity) sensitivities and the high-frequency trans-Planckian considerations.","marker":"[9]"}],"fun_headline_variants":["Extra bounce phase creates broken gravitational-wave spectrum","Next-to-minimal bounce predicts broken GW spectrum","Broken GW power law from extra bounce phase","Next-to-minimal bounce yields testable GW background"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the new early contraction phase changes only low-frequency modes, leaving the high-frequency spectrum exactly equal to the minimal bouncing result; if phase-0 effects leak above the pivot frequency $f_\\star$, or the deep-bounce approximation $k\\eta_{s\\downarrow}\\ll1$ fails near $f_\\star$, the automatic sub-Planckian conclusion changes.","fun_headline_variants_meta":{"raw":{"variants":["Extra bounce phase creates broken gravitational-wave spectrum","Next-to-minimal bounce predicts broken GW spectrum","Broken GW power law from extra bounce phase","Next-to-minimal bounce yields testable GW background"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000643,"raw_usage":{"total_tokens":3014,"prompt_tokens":1059,"completion_tokens":1955,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":675,"completion_tokens_details":{"reasoning_tokens":1896}},"tokens_in":675,"tokens_out":1955,"duration_ms":16595,"temperature":1.0,"reasoning_tokens":1896,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T16:35:02.279345+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically integrate the tensor perturbation equation across the five phases for the paper's four example parameter sets without imposing $k\\eta_{s\\downarrow}\\ll k\\eta_{0\\downarrow}\\ll1$ or the $k\\eta_{0\\downarrow}\\ll1$ approximation, and compare the exact spectrum with Eqs. (22) and (28); any material departure on the $f\\ge f_\\star$ branch would shift the $0.79\\,m_{\\rm pl}$ intersection and weaken the automatic sub-Planckian conclusion.","supporting_citations":[{"cited_title":"Interpreting Pulsar Timing Array data of Gravitational Waves with Ekpyrosis-Bouncing Cosmology","cited_arxiv_id":"2408.06582","evidence_quote":"Together with [50], provides the BBN/CMB ΔNeff constraint quoted in Eq. (24)."},{"cited_title":"Dual Inflation and Bounce Cosmologies Interpretation of Pulsar Timing Array Data","cited_arxiv_id":"2405.15889","evidence_quote":"Supplies laboratory detector (superconducting circuit and cavity) sensitivities and the high-frequency trans-Planckian considerations."}],"review_version":1}