{"id":"434d9bd1-b756-4e67-9191-228a7e4187d5","arxiv_id":"2412.10111","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"In the unstable branch of quadratic gravity, singularity-free bouncing cosmologies exist but their second-order action variation has negative eigenvalues, indicating instability and a possible transition to a Big Bang.","lead":"This paper studies universes in a modified theory of gravity with extra curvature terms and finds a family of solutions that never hit a Big Bang or Big Crunch singularity. It then argues, using a stability calculation, that these smooth universes are fragile and can turn into a Big Bang universe when the theory's parameters are nudged slightly.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The instability claim rests on pointwise Hessian eigenvalues in (a, ȧ, ä); for this higher-derivative constrained system that is not a valid dynamical stability criterion, and no nearby trajectory is actually shown to cross a = 0.","rationale":"To support the abstract's claim, it would suffice to show either (i) that the Hessian criterion is a correct and sufficient stability test for this constrained higher-derivative system, or (ii) that nearby solutions actually undergo a transition to a singular Big Bang. Neither is established. The Hessian claim is the least secure because quadratic gravity is known to have an unbounded-below Hamiltonian (Ostrogradsky ghost), so a pointwise negative eigenvalue in (a, ȧ, ä) is expected and cannot distinguish a generic pathology from instability of the chosen trajectory. The numerical evidence is also thin: Figure 10 is one solution, its parameters are merely said to be close to a transition, and no code or data are provided. The reader's CONDITIONAL verdict already flags the Hessian criterion; my concern agrees with that and sharpens it. I do not propose moving the verdict, because the current conditional status correctly asks for the missing dynamical check; if that check fails, the paper should be revised to a statement about ghost pathologies rather than solution instability.","tokens_in":9558,"tokens_out":9123,"duration_ms":102479,"concrete_test":"With the Figure 10 parameters (ω = 0, ρ0 = 1.66, Λ = 0.5, α = 4, β = 19, γ = 10, a(1) = 1, ȧ(1) = 0.8, ä(1) = 0.7), integrate the full fourth-order equation (15) for ρ0 in a small interval around 1.66 (for example 1.60 to 1.72 in steps of 0.01) and record whether a(t) crosses zero, and also integrate a few trajectories with initial (ȧ(1), ä(1)) perturbed by 10^-6 to 10^-2. In parallel, linearize Eq. (15) about the reported nonsingular solution and integrate the resulting fourth-order linear ODE for δa with smooth compact initial data. The central claim is corroborated if infinitesimal parameter or initial-data variations (or an exponentially growing linearized mode) drive a to zero; if nearby solutions remain nonsingular and bounded, the Hessian criterion in Section V is not a valid instability test.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section V computes the matrix (21) as the second-order variation of the action in the minisuperspace variables a, ȧ, ä and concludes instability from the negative eigenvalue λ1, whose eigenvector is identified with expansion. The load-bearing premise is that pointwise negativity of this Hessian proves dynamical instability. That premise is unproved and suspect. For the fourth-order Lagrangian (9) the second variation is an integral operator: δȧ and δä are time derivatives of δa, not independent directions, so diagonalizing ∂²L/∂x_i∂x_j at a fixed time does not give a Lyapunov spectrum. Negative curvature of the action is the usual Ostrogradsky/ghost symptom of quadratic gravity, not a statement that the specific nonsingular solution is repelled by the equations of motion. The text also says the eigenvector near tc has ȧ ≈ 0 and ä > 0, but that identifies an instantaneous direction in (a, ȧ, ä) space, not a growing perturbation. Finally, the paper asserts that the family behaves qualitatively the same without a parameter scan, and the abstract's claim that a Big Bang emerges under slight parameter changes is never demonstrated by integrating nearby parameters or perturbed initial data.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies homogeneous and isotropic FLRW solutions in quadratic gravity with action (1), focusing on the regime 3α < β where the de Sitter branch is unstable. It presents numerical solutions showing that this regime admits singularity-free bouncing universes (Figure 9, Figure 10) in addition to collapsing and singular solutions. The central claim is that these nonsingular solutions are unstable: the authors compute the second variation of the action with respect to (a, ȧ, ä), forming the matrix (21), and from its time-dependent eigenvalues (Figure 11) and the associated eigenvector (Figure 12) conclude that a Big Bang can emerge from the singularity-free regime when parameters are slightly modified. The paper also observes that on FLRW backgrounds the quadratic gravity action reduces to a Starobinsky-like f(R) form (19), since the Weyl tensor vanishes.","tokens_in":9869,"tokens_out":3388,"duration_ms":38967,"significance":"If the stability conclusion were rigorously established, the paper would make a useful contribution to the debate on singularity avoidance in higher-derivative gravity: it would show that the nonsingular classical solutions in the 3α < β branch are dynamically unstable, so that they cannot serve as a complete resolution of the initial singularity. The reduction to an f(R) form on FLRW backgrounds and the numerical identification of a large family of nonsingular solutions are interesting and potentially valuable. However, the paper's central claim rests on a stability criterion that is not justified for this constrained higher-derivative system, and the evidence is limited to one representative numerical solution with no error estimates or parameter scan. The paper is clearly written in its descriptive parts, but the stability analysis as presented does not meet the burden of proof for the strong statement in the abstract.","major_comments":[{"comment":"The Hessian-positivity criterion used to conclude instability is not established for this higher-derivative, diffeomorphism-invariant system. The second variation of the action is an integral functional in which δȧ and δä are time derivatives of δa, not independent directions at a fixed time. Diagonalizing the pointwise matrix ∂²L/∂x_i∂x_j in (a, ȧ, ä) space does not yield a Lyapunov spectrum, and a negative eigenvalue of this matrix does not by itself imply a growing solution of the linearized equations of motion. The manuscript must either prove that this criterion is equivalent to a standard dynamical stability analysis for this constrained fourth-order system, or replace it with a direct analysis of the linearized perturbation equations (e.g., solving the perturbed equations of motion with appropriate boundary conditions). Without this, the central instability claim in the abstract is not supported.","section":"Section V, Eq. (21) and Figures 10–12"},{"comment":"The stability verdict is based on a single numerically integrated solution, with no code, data, or error estimates provided. The text asserts that the family of nonsingular solutions behaves qualitatively the same, but no parameter scan or second example is shown. Moreover, the abstract claims that a Big Bang emerges when parameters are slightly modified, yet the paper never integrates nearby parameters or perturbed initial data to demonstrate that the scale factor actually crosses a = 0. Providing a small parameter scan around the Figure 10 solution and explicitly integrating a nearby trajectory that becomes singular would be necessary to substantiate that load-bearing claim.","section":"Section V, 'the behavior of the family ... is qualitatively the same' and Figure 10"},{"comment":"The identification of the eigenvector direction with 'expansion' is not meaningful without a valid dynamical interpretation of the Hessian eigenvectors. Even if the Hessian criterion were accepted, the statement that v1y ≈ 0 and v1z > 0 corresponds to expansion around the critical time conflates an instantaneous direction in the finite-dimensional configuration space (a, ȧ, ä) with the time evolution of a physical perturbation. The manuscript should clarify how the eigenvector components relate to a growing mode of the actual equations of motion, or avoid this interpretive claim.","section":"Section V, eigenvector analysis around t_c"}],"minor_comments":[{"comment":"The manuscript contains several typos and grammatical errors that should be corrected, e.g., 'homegeneous' (Section IV), 'wievpoint' and 'o the fact' (Section VI), 'the of the vanishing' (Section V), and 'under homogeneous perturbations are stable or not' (Section V).","section":"Throughout"},{"comment":"The matrix (21) is introduced without derivation, and the prefactor '24' as well as the layout of the matrix entries are not fully clear as typeset. The authors should present the explicit second-order variation of the action, including the treatment of boundary terms, so that the reader can verify the quadratic form.","section":"Section V, Eq. (21)"},{"comment":"The abstract states that 'the complete analysis shows that a Big Bang can emerge from a singularity-free regime,' while Section V concludes that the analysis 'suggests the existence of an instability.' The wording should be aligned; the stronger claim is not supported by the current evidence.","section":"Abstract and Section VI"},{"comment":"The relation ˙E_b = ˙a b E_a appears dimensionally odd; if this is a typo, it should be fixed, and if not, the notation should be defined more carefully.","section":"Section III, Eq. (17)"}],"recommendation":"major_revision","confidential_remarks":"The central instability claim rests on a stability criterion that is, at best, unproven for this constrained higher-derivative system and, at worst, invalid. If the authors can replace the Hessian-eigenvalue argument with a proper linearized stability analysis (or rigorously connect it to the equations of motion), the paper could become a useful contribution. As it stands, the abstract overstates what is demonstrated. I would recommend major revision rather than immediate rejection because the descriptive results (the f(R) reduction and the numerical identification of nonsingular solutions) are of some value, and the stability question is in principle addressable within the manuscript's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nRead the Asorey-Ezquerro-Pardina paper on singularity-free solutions in quadratic gravity. My take: the paper has a genuinely new piece of analysis—the explicit second-variation eigenvector study of the nonsingular FLRW branch—but the load-bearing stability conclusion is not supported as it stands. The Hessian matrix in (a, ȧ, ä) is diagonalized pointwise, and its negativity is read as dynamical instability. That is suspect for a constrained fourth-order system: δȧ and δä are time derivatives of δa, not independent directions, so the second variation is an integral operator. Pointwise negativity of that matrix is essentially the Ostrogradsky/ghost symptom, not a statement that this particular solution is repelled by the equations of motion.\n\nWhat the paper does well: it is clear about the reduction of quadratic gravity to a Starobinsky-like f(R) model on FLRW, the gauge-fixing issue, and the relation to the ghost. The numerical survey of solutions—big crunch, nonsingular bounces—is useful and reproduces known branches. The citations to [24] for the ghost interpretation and [22] for de Sitter stability are appropriate. The genuinely new element is the eigenvalue analysis, which could be meaningful if paired with a proper stability criterion.\n\nThe soft spots are real but addressable. There is no code or data, only one representative solution, and no demonstration that nearby parameters or perturbed initial data actually drive the scale factor through zero. The abstract's claim that a Big Bang emerges from a singularity-free regime under slight parameter changes is not shown by integrating nearby trajectories. The analysis is also limited to homogeneous perturbations. So the central instability verdict is conditional at best.\n\nIf I were refereeing, I would ask for: (i) a justification of the Hessian criterion (or a better stability notion), (ii) a parameter scan or a second example, and (iii) the actual perturbed evolution. The paper is not a throwaway—the question matters and the work is honestly presented. It deserves serious peer review, not desk rejection.\n\nOverall: worth engaging, but the instability claim needs much stronger support.","headline":"The instability claim rests on a questionable pointwise Hessian criterion, but the paper contains a genuinely new eigenvalue analysis and deserves a serious referee.","tokens_in":10357,"tokens_out":2241,"would_cite":false,"duration_ms":23878,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83F05","83D05"],"pacs":["04.50.Kd","98.80.Qc","04.20.Dw"],"model":"deepseek-v4-flash","headline":"Quadratic gravity's singularity-free universes are unstable and can tip into a Big Bang.","keywords":["quadratic gravity","singularity-free cosmology","Big Bang emergence","higher-derivative gravity","FLRW cosmology","de Sitter stability","ghost instability","minisuperspace stability"],"falsifier":"Compute the full linearized perturbation spectrum, including inhomogeneous scalar, vector, and tensor modes, around the nonsingular solution of Figure 10: if every physical gauge-invariant mode decays for the same parameters, the claimed instability is not established, while a growing inhomogeneous mode would make it stronger than claimed. A simpler check is to integrate nearby initial conditions around the same background and look for the first trajectory whose scale factor reaches zero.","tokens_in":9369,"feed_emoji":"💥","tokens_out":11723,"duration_ms":104606,"temperature":0.7,"pith_summary":"The paper studies homogeneous, isotropic universes in quadratic gravity -- Einstein-Hilbert action plus $R^2$ and $R_{\\mu\\nu}R^{\\mu\\nu}$ terms -- and claims that, in the parameter regime $3\\alpha < \\beta$, the theory admits a family of solutions with no Big Bang or Big Crunch singularity: the universe contracts to a finite minimum size and re-expands. It then computes the second variation of the action for such a solution and finds that the associated quadratic form is not positive definite. Near the minimal scale factor, the largest-instability eigenvector points in the direction of expansion ($\\dot a \\approx 0$, $\\ddot a > 0$). The paper concludes that these singularity-free solutions are unstable: a slight change of parameters makes a Big Bang emerge from a nonsingular regime. If correct, this means the singularity-avoiding power of quadratic gravity is not robust under homogeneous perturbations.","feed_headline":"Nonsingular bounces in quadratic gravity are unstable","feed_subtitle":"In the 3α<β regime, a slight parameter shift turns a bounce into a Big Bang.","key_machinery":"The load-bearing object is the Hessian matrix (21), the second variation of the gravitational action on the minisuperspace, i.e. the reduced configuration space parametrized by $(a,\\dot a,\\ddot a)$ for a flat FLRW metric. Its entries mix the parameters $\\gamma$, $\\Lambda$, $\\alpha$, $\\beta$, the matter density $\\rho_0$ and equation-of-state $\\omega$, together with the scale factor and its time derivatives. The stability verdict is read from the time-dependent eigenvalues and eigenvectors of this matrix evaluated on a background nonsingular solution: positive definiteness would mean stability, while the observed non-positive eigenvalues -- and the dominant eigenvector pointing to expansion at the bounce -- signal instability.","core_discovery":"The central claim is that the nonsingular FLRW solutions of quadratic gravity found in the unstable de Sitter branch $3\\alpha<\\beta$ are dynamically unstable. The evidence is the second-order variation of the reduced action: the Hessian matrix (21) in the variables $(a,\\dot a,\\ddot a)$ has a spectrum that is not positive definite, with the dominant unstable eigenvalue growing without bound at the time $t_c\\approx2.12$ when the scale factor reaches its near-zero minimum $a(t_c)\\approx0.028$. At that moment the associated eigenvector has $\\dot a\\approx0$ and $\\ddot a>0$, exactly the signature of expansion. The paper therefore asserts that perturbations drive a transition from a contraction-then-expansion nonsingular universe into a singular Big Bang solution, and that this instability is intrinsic to the singularity-free branch rather than an artifact of initial conditions.","pith_inferences":["The paper only perturbs within the homogeneous, isotropic ansatz; whether inhomogeneous scalar, vector, or tensor modes make the nonsingular branch more unstable, or instead stabilize it, remains untested.","Because the Hessian diverges as the scale factor approaches zero, the analysis suggests a testable gradient: nonsingular solutions with smaller minimal $a(t)$ should be progressively more unstable, a prediction that could be checked by scanning the family of solutions numerically.","If the instability survives a gauge-invariant treatment, bouncing and ekpyrotic scenarios built on this branch of quadratic gravity would need an external stabilizing mechanism or would not persist as long-lived cosmological histories."],"forward_implications":["In the regime $3\\alpha<\\beta$, the nonsingular cosmological solutions sit at a dynamical threshold: small parameter changes can turn them into solutions with a Big Bang.","The same analysis implies that the bounce-like phase is transient rather than periodically recurring, because the fastest-growing mode near the minimal scale factor drives expansion away from the contracting phase.","The instability is visible already within the homogeneous, isotropic ansatz, so restoring isotropy is not what protects the nonsingular branch; the destabilizing direction is present in the reduced action itself.","Since the condition $3\\alpha=\\beta$ reduces the FLRW action to an $f(R)$ model of the $R^2$ type, the nonsingular family belongs specifically to the sector where the scalaron mass is imaginary, linking the instability to the ghost-like content of the theory."],"supporting_citations":[{"why":"Establishes quadratic gravity as a renormalizable higher-derivative theory, the framework the paper works in.","marker":"[2]"},{"why":"Supplies the particle spectrum of the theory, including the scalaron and the spin-2 ghost whose masses depend on $\\alpha$ and $\\beta$.","marker":"[3]"},{"why":"Shows the FLRW reduction of quadratic gravity to an $f(R)$ model when $3\\alpha=\\beta$.","marker":"[21]"},{"why":"Provides the result that de Sitter space is stable only for $3\\alpha>\\beta$, defining the branch where nonsingular solutions appear.","marker":"[22]"},{"why":"Analyzes stability of isotropic cosmological singularities in higher-order gravity, the singular counterpart the paper contrasts with its nonsingular solutions.","marker":"[23]"},{"why":"Interprets singularity avoidance in terms of the ghost particle's repulsive energy, the mechanism the paper invokes for the nonsingular branch.","marker":"[24]"}],"fun_headline_variants":["Nonsingular bounces unstable: Big Bang emerges","Quadratic gravity bounce solutions are unstable","Perturbations send nonsingular bounces to Big Bang","Bounce to Big Bang: nonsingular cosmologies unstable","Small tweak makes nonsingular bounce go Big Bang"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The instability verdict assumes that positivity of the Hessian matrix (21) in the reduced variables $(a,\\dot a,\\ddot a)$ is the correct stability criterion for this constrained, higher-derivative gravitational system, and that the single numerically integrated solution of Figure 10 is representative of the whole nonsingular family.","fun_headline_variants_meta":{"raw":{"variants":["Nonsingular bounces unstable: Big Bang emerges","Quadratic gravity bounce solutions are unstable","Perturbations send nonsingular bounces to Big Bang","Bounce to Big Bang: nonsingular cosmologies unstable","Small tweak makes nonsingular bounce go Big Bang"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001218,"raw_usage":{"total_tokens":4951,"prompt_tokens":825,"completion_tokens":4126,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":441,"completion_tokens_details":{"reasoning_tokens":4046}},"tokens_in":441,"tokens_out":4126,"duration_ms":29702,"temperature":1.0,"reasoning_tokens":4046,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T16:20:10.711708+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the full linearized perturbation spectrum, including inhomogeneous scalar, vector, and tensor modes, around the nonsingular solution of Figure 10: if every physical gauge-invariant mode decays for the same parameters, the claimed instability is not established, while a growing inhomogeneous mode would make it stronger than claimed. A simpler check is to integrate nearby initial conditions around the same background and look for the first trajectory whose scale factor reaches zero.","supporting_citations":[{"cited_title":"De Felice and S","cited_arxiv_id":null,"evidence_quote":"Shows the FLRW reduction of quadratic gravity to an $f(R)$ model when $3\\alpha=\\beta$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the result that de Sitter space is stable only for $3\\alpha>\\beta$, defining the branch where nonsingular solutions appear."},{"cited_title":"Middleton and J","cited_arxiv_id":null,"evidence_quote":"Analyzes stability of isotropic cosmological singularities in higher-order gravity, the singular counterpart the paper contrasts with its nonsingular solutions."},{"cited_title":"Kuntz and R","cited_arxiv_id":null,"evidence_quote":"Interprets singularity avoidance in terms of the ghost particle's repulsive energy, the mechanism the paper invokes for the nonsingular branch."}],"review_version":1}