{"id":"186985cc-4ca7-4fac-830e-3ce24f664a8a","arxiv_id":"2607.23965","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"An explicit three-scalar model realizes a thermal stasis attractor with a finite lifetime, whose fast and slow trajectories make the duration of stasis strongly initial-condition dependent.","lead":"This paper builds an explicit particle-physics model in which a cosmological 'stasis' attractor—one that holds matter and radiation abundances fixed—operates only for a finite time. It shows that whether the universe ever reaches stasis can depend sensitively on initial conditions, and that it can spend many e-folds being pulled toward stasis without ever arriving.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Attractor dynamics check out, but 'satisfying all constraints' rests on an unspecified φ→SM decay operator not present in Eq. (3.27); absent that operator the model fails BBN/ΔN_eff.","rationale":"I agree with the reader's primary weakest assumption: the visible-sector decay operator is missing from the explicit model, and the claim to satisfy observational constraints is conditional on an unspecified extension. The reader's secondary concern about the GR correction coefficient b_X^(1) is not load-bearing, because the paper explicitly demonstrates in Section VI and Fig. 6 that the corresponding bound is vastly subleading; even O(10^6) coefficient would not move the constraint into the stasis wedge. I verified the central linear-algebra claims: the fixed point and eigenvalues for q=-2 are correct, and the logistic trajectory follows from Eq. (3.24). The analytic attractor dynamics therefore survive scrutiny. However, the abstract's broad statement that the model 'satisfies all relevant phenomenological and cosmological constraints' is not supported by the model as written; it is a promise about an operator that is never specified. This does not falsify the attractor mechanism, but it does mean the full central claim is conditional. The reader's CONDITIONAL verdict already captures this, so I recommend no change.","tokens_in":47643,"tokens_out":9268,"duration_ms":100724,"concrete_test":"Add an explicit Z2-breaking decay operator to Eq. (3.27), e.g., y φ ψ̄ψ with a small Yukawa y, and for the Fig. 7 benchmark parameters (gχ=1.25e-5, gφ=1.0, g_G=9.21e-10, μ satisfying Eq. (4.2)) scan y/Λ over the window t_end < t_ϕ < t_BBN. Require that Eq. (4.87) holds with ΔN_eff<0.16, that Γ_φ<H during the stasis epoch so φ is not depleted before the expiration date, and that the reheating temperature from φ decay exceeds T_BBN ≈ 10 MeV. Re-run the Boltzmann evolution with this operator included and check whether any choice of y/Λ simultaneously satisfies all three conditions. If no such choice exists, the assumption in Section IV.G is unrealizable and the 'satisfying all constraints' claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The core attractor calculation (Eqs. 3.13–3.25) is internally consistent: for q=-2 the fixed point is Ω_M=2/5 with eigenvalues (-1/5,-1), and the logistic fastest trajectory is a genuine feature of the two-dimensional flow. The load-bearing gap is the abstract's claim that this model 'satisfies all relevant phenomenological and cosmological constraints.' That claim depends entirely on Section IV.G, which assumes—without adding any term to the Lagrangian in Eq. (3.27)—that φ couples to visible-sector fields through 'highly suppressed operators' and decays after stasis to reheat the SM before BBN. The exact Z2 symmetries imposed in Section III.E forbid any such coupling; a coupling that permits φ→SM decay requires an explicit Z2-breaking extension. If that operator is absent, φ is stable and the late universe remains dominated by dark radiation plus φ, violating BBN and ΔN_eff bounds. If it is present, its strength and back-reaction during the stasis epoch must be checked; a decay rate that is not sufficiently suppressed will deplete φ before the expiration date and change the pump. Thus the 'satisfying all constraints' portion of the central claim is not a property of the model as defined; it is a conditional on an unspecified extension. The uncomputed GR coefficient b_X^(1) in Section IV.B is a lesser concern: the paper shows the associated bound is subleading, and Fig. 6 confirms the contour lies far outside the stasis wedge.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes a three-real-scalar model (ϕ, X, χ) that realizes a thermal stasis attractor. The core dynamics are reduced to a two-dimensional flow in the matter abundance Ω_M and a 'coldness' variable Ξ. For an annihilation cross-section σv ∝ |p_CM|^{-2} (q = -2), the fixed point is Ω_M = 2/5 with Jacobian eigenvalues (-1/5, -1), so the fixed point is attracting. The authors identify a 'fastest' trajectory that is exactly logistic, derive approximate analytic trajectories for the approach to stasis, and show that the validity of the q = -2 form is restricted to a finite temperature window [T_min, T_max], giving an 'expiration date' (and a 'manufacture date'). They then impose a long list of internal consistency conditions and observational bounds, identify a 'stasis wedge' in the (τ, ϱ_M) plane, and show numerically that the duration of stasis is sensitive to initial conditions, ranging from about 2 to 10 e-folds in the benchmark cases of Fig. 7.","tokens_in":48149,"tokens_out":4415,"duration_ms":51873,"significance":"If the model is taken as stated, the paper makes a useful conceptual point: a stasis attractor need not produce a stasis epoch; the same attractor can pull the system for many e-folds without ever reaching stasis, and the duration of the realized stasis epoch can depend sharply on initial conditions. The analytic core is a strength: the eigenvalue computation, the logistic solution along the fastest trajectory, and the explicit initial-condition conditions (Eqs. (5.6), (5.20), (5.21)) are transparent and checkable. The paper is also honest about the extent to which it builds on Ref. [10], and the constraint analysis is unusually explicit, with Tables I and II summarizing the conditions that define the allowed region. The main weakness is that the claim of 'satisfying all relevant phenomenological and cosmological constraints' is conditional on a φ→SM decay mechanism that is not part of the defined Lagrangian; absent that mechanism, the late universe in the model is dark-sector dominated and fails BBN/ΔN_eff bounds. This is a fixable but load-bearing gap.","major_comments":[{"comment":"The abstract and Introduction claim the model 'satisfies all relevant phenomenological and cosmological constraints,' but the only mechanism that transfers energy to the visible sector is introduced conditionally in Section IV.G: 'If φ couples to the fields of the visible sector via highly suppressed operators...'. The Lagrangian in Eq. (3.27) contains no such operator, and the two Z2 symmetries imposed in Section III.E forbid any renormalizable φ-SM coupling (φ is odd under the first Z2). Without an explicit Z2-breaking extension, φ is stable, no visible-sector reheating occurs, and the post-stasis universe is dominated by dark radiation plus the stable ϕ gas. The condition is therefore a property of an unspecified extension, not of the model as defined. The manuscript should either add an explicit φ→SM operator to the Lagrangian and re-derive the stasis-pump and late-time constraints w","section":"Section IV.G and Eq. (3.27)"},{"comment":"Even if one accepts that a highly suppressed φ→SM decay operator is present, the analysis does not check the back-reaction of that decay during the stasis epoch. Equations (4.83) through (4.87) assume that φ decays can be neglected until t_ϕ and that ρ_M thereafter redshifts as matter until decay reheats the SM. But a nonzero Γ_φ depletes ρ_M during stasis, modifying the pump equation (3.13) through an additional loss term in dρ_M/dt and an additional source term in dρ_γ/dt. The benchmark durations N_s quoted in Section VI.B and the BBN bound in Eq. (4.87) assume Φ decays only after t_end. The authors should demonstrate that for the benchmark parameters the φ decay rate satisfies Γ_φ ≪ H throughout the stasis window, and that the φ abundance at t_end is consistent with the value used in Eq. (4.87).","section":"Section IV.G (decay back-reaction)"}],"minor_comments":[{"comment":"The coefficient b_X^(1)(p^2) is not computed, and the text says 'deriving it is beyond the scope of this paper.' The order-of-magnitude assumption is probably harmless because Fig. 6 shows the corresponding contour is very subleading, but this insensitivity should be stated quantitatively in the text: for example, the contour in Fig. 6 would need to move by many orders of magnitude before affecting the stasis wedge.","section":"Section IV.B, Eq. (4.18)"},{"comment":"The Monte-Carlo coefficient ε ≈ 0.669 is quoted without describing the integration method, the number of samples, or the statistical uncertainty. Since ε enters the constraint contours (4.52)–(4.55) and hence the boundary of the 'stasis wedge' in Fig. 6, the numerical procedure should be documented, or the text should clearly state that only order-of-magnitude accuracy is intended.","section":"Section IV.C, Eq. (4.48)"},{"comment":"The arrival criterion depends on the arbitrary cutoff δ. The text uses δ = 0.001 in Fig. 1 but should state the value used for the N_s values in Fig. 7 and confirm that the reported e-fold counts are not strongly δ-dependent. A one-sentence sensitivity test would settle this.","section":"Section III.D, Eq. (3.22)"},{"comment":"The footnote about 'nine horizontal lines' is not appropriate for a journal article and should be removed.","section":"Eq. (4.51), footnote"}],"recommendation":"major_revision","confidential_remarks":"I agree with the reader's and skeptic's assessment. The attractor calculation itself is sound and the central dynamical result is interesting, but the advertised claim of satisfying all phenomenological and cosmological constraints is not yet a property of the model as defined: it depends on an unspecified Z2-breaking φ→SM decay operator. This is a significant but potentially fixable issue, so I recommend major revision rather than rejection. I do not see a circularity problem: the model is deliberately constructed to have q = -2, and the new results about initial-condition sensitivity and expiration dates are derived within that construction."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The headline: this paper is worth reading for anyone working on early-universe stasis, but its central 'satisfies all relevant constraints' claim is not a property of the model as written. It depends on a decay operator for φ into SM states that never appears in the Lagrangian of Eq. (3.27), and the Z2 symmetries imposed there forbid it. The stress-test note is right about that.\n\nWhat is actually new: the authors take their earlier thermal-stasis mechanism (Ref [10]), re-derive the attractor cleanly, and then show that the attractor has a finite temperature window, i.e., a shelf life. The logistic fastest trajectory (Eq. 3.25), the analytic reachability condition Eq. (5.21), and the 'stasis wedge' constraint map are new and internally consistent. The linearized dynamics for q=-2 are simple: fixed point Ω_M=2/5, eigenvalues (-1/5,-1). The figures and analytic approximations agree, and the paper is transparent about which pieces come from earlier work.\n\nSoft spots, in order of size. The missing φ→SM decay operator is the load-bearing one. Without it, φ is stable, and the late universe is dark radiation plus φ, which violates BBN and ΔN_eff. The authors state the assumption explicitly in Section IV.G, but the abstract and several summary statements don't carry the caveat. This is addressable—add a tiny Z2-breaking coupling to SM fermions, check back-reaction—but it isn't done here, so the constraint-satisfaction claim is conditional. Second, the uncomputed GR correction b_X^(1) in IV.B is a real gap, but the authors show the resulting bound is subleading and Fig. 6 confirms the contour sits far outside the stasis wedge; I don't think it changes the results. Third, the Monte-Carlo coefficient ε≈0.669 in Eq. (4.48) is quoted without code, but the analytic upper bound preceding it gives the same order of magnitude, so I treat that as minor.\n\nThere are also a few places where the text overuses phrases like 'wide range of behaviors,' but the actual results support the main claim about initial-condition sensitivity: N_s ranges from ~10.5 e-folds if you start in stasis to ~2 e-folds for typical starts. That is a concrete, useful observation.\n\nWho this is for: BSM cosmologists interested in pre-BBN dynamics and anyone working on stasis or attractors in cosmology. It deserves a serious referee: the core dynamics are checkable and the flaw is a missing ingredient, not a math error. I would send it to peer review with a request to fix the decay-operator issue, either by adding it to the model or by softening the abstract.","headline":"A solid, transparent derivation of the thermal stasis attractor's finite shelf life, undermined by an unmodeled φ→SM decay operator that the 'satisfies all constraints' claim depends on.","tokens_in":48560,"tokens_out":8269,"would_cite":true,"duration_ms":59595,"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 claims that a three-scalar model with a finite-temperature stasis attractor produces stasis epochs whose length—and even occurrence—depends sharply on initial conditions, ranging from about 10 e-folds to none.","keywords":["thermal stasis","cosmological attractor","early universe","scalar field model","matter-radiation stasis","annihilation pump","expiration date","big-bang nucleosynthesis constraints"],"falsifier":"Compute the one-loop coefficient b_X^(1)(p^2) of the H-dependent correction to the mediator propagator in the FRW background; if it is not O(1) but instead large enough that the GR term dominates within the temperature window for the paper's benchmark couplings, the q = −2 scaling and the attractor fail. Equivalently, for a specified ϕ→SM decay operator, compute the BBN dark-radiation abundance: if ΔN_eff exceeds about 0.16 for any consistent decay time, the model is excluded.","tokens_in":47559,"feed_emoji":"⏳","tokens_out":7442,"duration_ms":79103,"temperature":0.7,"pith_summary":"Cosmological stasis is a proposed early-universe epoch in which matter and radiation abundances stay fixed even though the universe keeps expanding. This paper tries to establish that stasis can come from an explicit, self-consistent particle-physics model: three real scalar fields with an s-channel annihilation pump that makes the annihilation rate grow as the particles cool. The resulting attractor pulls the system toward a fixed 2/5 matter, 3/5 radiation split, but only while the gas temperature lies in a finite window; below the lower end the attractor expires. The paper's central quantitative claim is that the number of e-folds of stasis depends sensitively on initial conditions: about 10.5 e-folds when the system starts at the fixed point, about 6.3 along the fastest trajectory, roughly 2 e-folds when the matter abundance starts at unity, and zero for many other starts—though those universes still spend time under the attractor's influence. A sympathetic reader cares because this is a concrete existence proof that near-stasis and partial-stasis cosmologies are viable and observationally distinguishable from standard pre-BBN history.","feed_headline":"Stasis can last 10 e-folds, or never happen","feed_subtitle":"A new three-scalar model gives stasis an expiration date, so its duration swings with initial conditions.","key_machinery":"The central object is the 'coldness', Ξ ≡ T^q ρ_M / m^{q+4}, which is the quantity that stays constant during stasis; for q = −2, Ξ = ρ_M / (m^4 τ^2). Together with Ω_M, it forms a two-variable dynamical system whose Jacobian at the fixed point has eigenvalues (−1/5, −1), creating fast and slow approach directions. The microscopic engine is the s-channel annihilation process ϕϕ → X → χχ, whose one-loop mediator propagator develops a term linear in the center-of-mass momentum, making the swept-volume rate scale as |p_CM|^{-2} (q = −2) only within a temperature window; the lower end of that window is the attractor's expiration date.","core_discovery":"The central claim is that the thermal stasis mechanism can be realized in a three-scalar model with an s-channel annihilation pump, and that this realization satisfies the relevant internal consistency and observational constraints. For q = −2, the fixed point is Ω_M = 2/5, with Jacobian eigenvalues (−1/5, −1), so the fixed point is a genuine attractor; the fastest approach is exactly logistic, ds/dN = s(1 − s). Because the q = −2 scaling of the annihilation cross-section only holds for T_min ≲ T ≲ T_max, the attractor has a finite lifetime, with a 'manufacture date' and an 'expiration date'. Consequently, the duration of stasis depends on initial conditions: near the fixed point or along th","pith_inferences":["The authors leave implicit that the same initial-condition sensitivity means any attempt to use stasis to solve pre-BBN puzzles must specify initial conditions, not just particle-physics parameters; otherwise e-fold counts are not predictive.","The expiration-date structure suggests a model-building principle: any stasis realized through a momentum-dependent cross-section has a finite temperature shelf life, so 'manufacture' and 'expiration' temperatures are tunable design parameters that could place stasis in observationally interesting windows.","The paper sets aside density-perturbation growth; an immediate testable extension would be to compute whether near-stasis epochs generate enhanced halo formation, which would modify the annihilation pump and could tighten or shift the allowed parameter space.","The fastest trajectory's exact logistic solution, which the paper notes also appears in another known stasis realization, invites a test of whether a universal fastest-path structure governs stasis attractors more generally."],"forward_implications":["In a broad class of initial conditions, the universe spends a significant number of e-folds under the attractor yet never reaches stasis; such 'near-stasis' cosmologies still differ from standard matter/radiation domination before BBN.","The stasis epoch has an upper bound set by the temperature window: even from the best initial condition, this model gives about 10.5 e-folds, providing a benchmark for observational searches.","The fastest approach to stasis is exactly solvable as a logistic equation, so along that trajectory the number of e-folds to stasis can be predicted analytically.","The model satisfies its consistency constraints—kinetic equilibrium, negligible 4→2 annihilation, no Bose-Einstein condensate, negligible mediator relics, thermodynamic limit, GR corrections, and BBN/dark-radiation bounds—in large regions of parameter space, with the BBN constraint not excluding any shown region of the (τ, ϱ_M) plane.","Late-time dark-radiation bounds force reheating via ϕ decay to occur after stasis ends; otherwise too much dark radiation survives to BBN."],"fun_headline_variants":["Stasis gets an expiration date in three-scalar model","How long stasis lasts hinges on initial conditions","Thermal stasis: finite pull, not eternal","Three scalars give stasis a variable lifetime","Stasis duration swings with where you start"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that ϕ decays to visible-sector states through couplings that are assumed but not part of the model's Lagrangian (and, secondarily, that the uncalculated GR correction to the mediator propagator is O(1)); if either fails, the universe is not reheated before BBN or the q = −2 window closes, and the claim of satisfying all constraints collapses.","fun_headline_variants_meta":{"raw":{"variants":["Stasis gets an expiration date in three-scalar model","How long stasis lasts hinges on initial conditions","Thermal stasis: finite pull, not eternal","Three scalars give stasis a variable lifetime","Stasis duration swings with where you start"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000403,"raw_usage":{"total_tokens":1920,"prompt_tokens":711,"completion_tokens":1209,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":455,"completion_tokens_details":{"reasoning_tokens":1135}},"tokens_in":455,"tokens_out":1209,"duration_ms":9048,"temperature":1.0,"reasoning_tokens":1135,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-31T23:23:30.134451+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the one-loop coefficient b_X^(1)(p^2) of the H-dependent correction to the mediator propagator in the FRW background; if it is not O(1) but instead large enough that the GR term dominates within the temperature window for the paper's benchmark couplings, the q = −2 scaling and the attractor fail. Equivalently, for a specified ϕ→SM decay operator, compute the BBN dark-radiation abundance: if ΔN_eff exceeds about 0.16 for any consistent decay time, the model is excluded.","supporting_citations":[],"review_version":1}