{"id":"61182c04-17de-416d-b714-1cdf635ae5ba","arxiv_id":"2501.11048","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A thesis that re-presents the author's prior published work, concluding that f(T,B), f(T,T_G), f(T,B,T_G,B_G), and scalar-field variants admit stable dark-energy attractors only for carefully hand-picked parameter ranges.","lead":"This doctoral thesis builds dynamical-system models of modified teleparallel gravity, adding scalar fields and boundary terms to torsion-based gravity and checking which phases of cosmic expansion appear as stable critical points. A generalist would read it as a consolidated test of whether boundary-term modifications of teleparallel gravity can replace the cosmological constant as an explanation for cosmic acceleration.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Constant-λ closure in the f(T,B) and f(T,B,T_G,B_G) analyses restricts the phase space to exact power-law solutions, so the claimed radiation→matter→de Sitter attractor sequence is not established for the full dynamics.","rationale":"The reader's weakest assumption—that λ=Ḧ/H^3 is treated as constant—is exactly the most load-bearing technical weakness in the thesis. The central claim is that modified teleparallel models provide a viable framework for both early and late cosmic acceleration, which requires a dynamical trajectory that passes through radiation, matter, and dark-energy domination. The constant-λ ansatz is not a harmless parameter choice because λ is not a free coupling but a derived kinematic quantity: its value changes as the universe transitions between epochs, and the dynamical equations used to classify critical points either omit this evolution or replace it with a fixed number. This affects the validity of the critical points themselves and, more seriously, the stability classifications that support the 'late-time attractor' claims. The test I propose directly addresses whether the concern lands: promoting λ to a dynamical variable or integrating the original second-order equations would reveal whether λ(N) is approximately constant along relevant trajectories. If it is not, the phase portraits and tables in Chapters 4 and 5 describe a different, restricted system, and the viability conclusion would need to be downgraded or explicitly qualified. Because the reader's own conditional verdict already hinges on this same assumption, my stress-test does not change the verdict: the thesis remains CONDITIONAL pending this check. I found no separate internal inconsistency that would warrant rejection beyond this closure issue, and no independent evidence (such as machine-checked proofs or released code) that would override it.","tokens_in":67776,"tokens_out":3817,"duration_ms":46061,"concrete_test":"Promote λ to a dynamical variable and recompute the stability analysis of §4.3. Specifically, differentiate the f(T,B) Friedmann constraint or the second-order field equations (4.2) for the logarithmic model to obtain dλ/dN as a function of X,Y,Z,W, and integrate the enlarged autonomous system together with the original equations using the parameter values and initial conditions of Figs. 4.4–4.5. Then compare λ(N) with the assumed constant values 8, 9/2, and 0 along the radiation→matter→de Sitter trajectory, and check whether the de Sitter critical point remains a global attractor of the enlarged system. If λ(N) changes by more than a few percent over the transition, or if the attractor structure changes, the constant-λ critical-point analysis does not support the claimed viability.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central viability claim requires a single dynamical system capable of describing radiation-, matter-, and dark-energy-dominated epochs in sequence. The f(T,B) analysis in §4.3 achieves closure by declaring λ=Ḧ/H^3 constant 'throughout the analysis' (the paragraph before Eq. 4.10). For a power-law scale factor a∝t^p, λ=2/p^2, so the quoted values λ=8 (radiation), λ=9/2 (matter), and λ=0 (de Sitter) are exact only for isolated power-law solutions. A real transition between these epochs has a time-dependent λ; no equation for its evolution is supplied in §4.3 or in the f(T,B,T_G,B_G) chapter that builds on the same assumption. Consequently, Eqs. (4.10) and the subsequent critical-point/stability tables describe a restricted system in which λ is frozen, not the actual phase flow of the cosmological model. The f(T,T_G) mixed-power-law model partially avoids this by solving λ from the variables, but the summary claim treats all models as one viable framework, and the constant-λ construction is load-bearing for the f(T,B) and f(T,B,T_G,B_G) results. The stability classifications of the de Sitter points C4/P4, which rely on zero-eigenvalue arguments plus numerical phase portraits, are particularly exposed: without a correct closure for λ, these points need not be attractors of the full system.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This PhD-thesis manuscript compiles dynamical-system analyses of several modified teleparallel gravity theories: teleparallel-Horndeski scalar-tensor models with power-law kinetic couplings (Ch. 2), f(T,φ) gravity with power-law torsion coupling and exponential/power-law potentials (Ch. 3), f(T,B) and f(T,T_G) gravity with boundary/Gauss-Bonnet terms (Ch. 4), the combined f(T,B,T_G,B_G) formalism (Ch. 5), and nonminimal boundary-coupled quintessence with observational Hubble and supernova comparisons (Ch. 6). The recurring method is to introduce dimensionless variables, form an autonomous system, locate critical points, compute eigenvalues, and identify which points can represent radiation-, matter-, and dark-energy-dominated epochs. The abstract concludes that these models provide a viable framework for early and late-time cosmic acceleration.","tokens_in":68163,"tokens_out":5669,"duration_ms":63733,"significance":"If the dynamical-system results were fully sound, the thesis would give a broad and useful phenomenology of teleparallel gravity modifications, with explicit stability tables, phase portraits, and some comparison to H(z) and distance-modulus data. Strengths of the manuscript include the unusually explicit algebraic presentation of the critical points and eigenvalues, the large number of models treated in a unified notation, and the inclusion of observational evolution plots in Chapter 6. However, the central viability claim is weakened by two load-bearing technical issues: the constant-λ closure in the f(T,B) and f(T,B,T_G,B_G) analyses, and the incomplete stability classification of non-hyperbolic de Sitter points. The present-day values Ω_m≈0.3 and Ω_DE≈0.7 are obtained from hand-picked initial conditions and parameter values rather than from a fit, so they are illustrative consistency checks rather than falsifiable predictions. The thesis is therefore a valuable collection of model studies, but the advertised single-framework explanation of the full radiation→matter→dark-energy sequence is not yet established.","major_comments":[{"comment":"The paragraph preceding Eq. (4.10) states that λ=Ḧ/H^3 is 'treated as a constant throughout the analysis.' For a power-law scale factor a∝t^p one has λ=2/p^2, so the values λ=8 (radiation) and λ=9/2 (matter) correspond exactly to isolated power-law solutions. During a physical transition between these epochs, λ is time-dependent, and no evolution equation for λ is supplied for the f(T,B) system or for the f(T,B,T_G,B_G) system in Chapter 5 that builds on the same assumption. Consequently, the critical points and stability tables in §4.3 describe a restricted system with λ frozen, not the actual phase flow of the cosmological model. The f(T,T_G) mixed-power-law model in §4.5.1 is less affected because λ is solved from the variables in Eq. (4.28), but the summary claim treats all models as one viable framework, and for the f(T,B) and f(T,B,T_G,B_G) analyses the constant-λ closure is load-bearing.","section":"§4.3, Eqs. (4.10) and surrounding text"},{"comment":"The de Sitter critical points C4 and P4 have eigenvalues {0,0,0,-4}, so linear stability theory is inconclusive. The text states that the central-manifold condition is not satisfied and therefore the method is not applied, but then concludes from 2D phase portraits that the point is an attractor. A two-dimensional projection of a four-dimensional flow cannot establish asymptotic stability; at most it suggests attraction in the projected directions. Without a valid center-manifold reduction, a Lyapunov function, or a rigorous normally-hyperbolic argument, the late-time attractor status of these points is unproven, and this status is central to the claimed dark-energy epoch.","section":"§4.3.1, point C4, and §4.3.2, point P4 (Tables 4.2, 4.6 and Figs. 4.3, 4.6)"},{"comment":"The presented present-day values Ω_m≈0.3, Ω_DE≈0.7, q0≈-0.61, and the transition redshifts are obtained by evolving the autonomous system from a specific hand-picked initial condition (e.g., x0=10^-8.89, y0=10^-2.89, u0=10^-5.96, ρ0=10^-0.9) with chosen parameters β=-0.2, σ=-0.30, λ=-0.2. No parameter estimation, likelihood, or goodness-of-fit is performed, so statements that these values 'agree with Planck' or are 'compatible with ΛCDM' are consistency checks rather than tests of the model. This does not invalidate the dynamical-system analysis, but it should be reframed in the abstract and conclusions as an illustrative compatibility exercise, not as evidence that the model is observationally preferred.","section":"§3.3.1 and §3.3.2, Figs. 3.3-3.4 and 3.7-3.8"},{"comment":"The thesis repeatedly infers a viable cosmic history from the coexistence of a stable dark-energy point and saddle radiation/matter points in the same parameter range. Strictly, one must also exhibit a trajectory (heteroclinic orbit or explicit numerical evolution) that connects the radiation saddle to the matter saddle and then to the dark-energy attractor. The 2D phase portraits and the single evolution plots from selected initial conditions do not establish this sequence for the full higher-dimensional systems. For the central claim of a single model describing radiation→matter→dark-energy in sequence, at least one representative heteroclinic chain or a full-dimension numerical trajectory with stated initial conditions is needed.","section":"Chapters 2-5 (general method)"}],"minor_comments":[{"comment":"The sentence 'Teleparallel gravity, is an alternative to General Relativity, explains gravitation through torsion' contains a misplaced comma; it should read 'Teleparallel gravity, an alternative to General Relativity, explains gravitation through torsion.'","section":"Abstract"},{"comment":"The tetrad for flat FLRW spacetime is written as (1,a(t),a(t),a(t)), which is ambiguous; it should be written as the diagonal tetrad diag(1,a(t),a(t),a(t)) to avoid confusion with a four-vector.","section":"Eq. (1.15)"},{"comment":"There are missing closing parentheses in references to 'action equation Eq. (2.11' and similar; please proofread the equation cross-references.","section":"§2.3.1 and §2.3.2"},{"comment":"The cosmological model name is written inconsistently as FLRW, FLRW, and 'FLR W'; one standard form should be used throughout.","section":"Throughout"},{"comment":"The subsection heading 'Radiaiton-dominated Critical points' contains a typo: it should be 'Radiation-dominated Critical points.'","section":"§4.3.2 heading"},{"comment":"Some phase-portrait captions list parameter values such as τ, ζ, δ without defining them in the caption or immediately preceding text; please clarify each symbol when it first appears in a figure.","section":"Figure captions, Chapter 2"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a compilation of previously published articles (EPJC 82, 680; EPJC 83, 809; EPJC 84, 1088; Annals of Physics 460, 169563). If the journal expects original unpublished content, the degree of overlap and the novelty relative to these papers must be clarified. The constant-λ issue is not a cosmetic defect: it affects the core viability claim for the f(T,B) and f(T,B,T_G,B_G) systems. The author should either provide a dynamical closure for λ, restrict the claims appropriately, or present the results as conditional on λ being constant."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Key takeaway: this is a PhD thesis assembled from four already-published papers, and it doesn't pretend otherwise. The useful part is the worked machinery: explicit critical-point tables, full eigenvalue lists, and phase portraits for the f(T,B), f(T,T_G), f(T,B,T_G,B_G), and scalar-torsion models. Anyone learning dynamical systems in modified teleparallel gravity will find it a decent reference. It is not a new research claim on top of the papers it compiles.\n\nThe soft spot is load-bearing. In Chapter 4, right before Eq. (4.10), the text defines λ = Ḧ/H^3 and says it is treated as a constant throughout. For power-law scale factors, λ = 2/p^2, so the quoted values λ=8 (radiation), 9/2 (matter), 0 (de Sitter) correspond to isolated power-law solutions. A sequence of epochs means λ evolves; no evolution equation for λ is supplied. So the critical-point and stability tables describe a reduced, frozen-λ system, not the phase flow of the full cosmological model. The abstract's blanket claim that these models are a viable framework for early and late acceleration is stronger than what the analysis supports. The f(T,T_G) mixed power-law section does derive λ from the variables, which is better, but the summary lumps all models together.\n\nSecond, the non-hyperbolic de Sitter points C4/P4 are called stable attractors even though the text reports that the center-manifold condition is not satisfied. The classification then rests on numerical phase portraits. That may be true, but without a Lyapunov function or a valid center-manifold reduction it isn't established. Third and minor: the H(z)/SNe Ia comparison in Chapter 6 is visually fine but statistically thin—no chi-square, no error budget, no evidence of a fit rather than a tune.\n\nTo the author's credit, the algebra is explicit and checkable, the provenance of each chapter is stated, and the weak spots I just named are not hidden. As a thesis the compilation is defensible. As a new journal submission it is double publication with an unsupported central claim. I would not send it to referees; if it were a thesis defense I'd pass conditionally and require the λ-dynamics issue to be addressed before the unified viability statement is kept.","headline":"Competent compilation of already-published dynamical-system papers, but the central viability claim rests on freezing λ=Ḧ/H^3 to isolated power-law trajectories.","tokens_in":68678,"tokens_out":3807,"would_cite":false,"duration_ms":43720,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83D05","83F05","37C75"],"pacs":["04.50.Kd","95.36.+x","98.80.-k"],"model":"deepseek-v4-flash","headline":"This thesis claims that modified teleparallel gravity theories can reproduce the observed cosmic sequence—radiation, matter, and accelerating dark-energy eras—as stable critical points of a single dynamical system, with late-time…","keywords":["modified teleparallel gravity","dynamical system analysis","cosmic acceleration","dark energy","boundary terms","teleparallel Gauss-Bonnet","scalar field cosmology","cosmological phase space"],"falsifier":"Compute $\\lambda = \\ddot{H}/H^3$ along the best-fit $H(z)$ curves the thesis reports for the $f(T,B)$ and $f(T,T_G)$ models; if $\\lambda$ moves appreciably across the radiation-to-dark-energy epochs instead of holding near the values used to find the critical points, such as $8$ for radiation and $9/2$ for matter, the fixed-point analysis does not govern the models' real trajectories.","tokens_in":67526,"feed_emoji":"🌌","tokens_out":7719,"duration_ms":83420,"temperature":0.7,"pith_summary":"This thesis sets out to show that modified teleparallel gravity—a family of theories that describes gravity through torsion rather than curvature—can account for the full sequence of cosmic phases, from radiation and matter domination to the present accelerated expansion, without invoking a cosmological constant. The author constructs autonomous dynamical systems for several such theories, including scalar-field-coupled models and extensions with boundary and Gauss-Bonnet terms, and identifies critical points whose stability matches the expected epochs of cosmic history. For the models examined, stable late-time attractors with equation-of-state parameter near minus one emerge, and the evolution of density parameters and the deceleration parameter comes close to standard Lambda-CDM behaviour. The thesis concludes that these modified teleparallel models form a viable framework for late-time cosmic acceleration.","feed_headline":"Modified teleparallel gravity can explain cosmic acceleration","feed_subtitle":"Dynamical-system analysis yields stable matter and dark-energy phases with present-day parameters near Lambda-CDM values.","key_machinery":"The engine of the analysis is the autonomous dynamical system: the Friedmann and Klein-Gordon equations of each theory are rewritten in dimensionless phase-space variables, and the fixed points of the resulting ordinary differential equations are classified by the eigenvalues of the Jacobian matrix, with center-manifold theory used for non-hyperbolic points. A constant dimensionless parameter $\\lambda = \\ddot{H}/H^3$ closes the system in the boundary-term chapters, and exponential and power-law potentials fix the scalar-field sector. The distinctive load-bearing objects are the teleparallel boundary term $B$ and the teleparallel Gauss-Bonnet term $T_G$, which enter the Lagrangians as new couplings and reshape the phase space; their presence is what distinguishes these models from ordinary $f(T)$ gravity and allows the de Sitter and scaling solutions the thesis highlights.","core_discovery":"The central claim is that, across the studied families of modified teleparallel gravity—$f(T,\\phi)$, $f(T,B)$, $f(T,T_G)$, $f(T,B,T_G,B_G)$, and the general teleparallel scalar-tensor formalism—the cosmological field equations can be rewritten as autonomous dynamical systems whose fixed points reproduce radiation-, matter-, and dark-energy-dominated eras in a single phase space. The radiation and matter critical points are saddles, so the Universe can pass through them, while the dark-energy critical points are stable late-time attractors; de Sitter solutions appear as limiting cases. The thesis reports that for representative parameter choices the models yield present-day density parameters near $\\Omega_m \\approx 0.3$ and $\\Omega_{\\mathrm{DE}} \\approx 0.7$, a transition from deceleration to acceleration at redshift $z \\approx 0.6$, and Hubble and distance-modulus curves compatible with observational fits. If correct, this means the observed cosmic acceleration can be explained by the geometry of torsion with added boundary couplings, without a cosmological constant.","pith_inferences":["Inference: The viability claim is not yet tested against the full cosmological data set; a natural next step is to confront these models with CMB, baryon-acoustic-oscillation, and large-scale-structure likelihoods, which the thesis does not do.","Inference: If the constant-$\\lambda$ assumption is relaxed, the fixed-point structure of the boundary-term chapters can change qualitatively; integrating the full field equations numerically along the same trajectories would show which stability classifications survive.","Inference: Because the $f(T,B)$ power-law models are described as useful for addressing the $H_0$ tension, the same dynamical-system machinery could be used to check whether the stable attractors shift the inferred Hubble constant, a comparison the thesis leaves implicit.","Inference: The same phase-space method could be applied to teleparallel models with direct dark-matter and dark-energy interactions, which the thesis does not investigate."],"forward_implications":["If the central claim holds, modified teleparallel gravity can generate the full radiation-to-matter-to-accelerating-dark-energy sequence from a single Lagrangian, with no cosmological constant input.","The stable dark-energy critical points mean the models converge to an accelerating late-time state from a broad set of initial conditions, making the acceleration phase an attractor rather than a fine-tuned choice.","The boundary and Gauss-Bonnet terms alter the phase-space structure and produce de Sitter and scaling solutions, offering a dynamical route to the coincidence problem within the scope of the paper.","The scalar-field models with exponential and power-law potentials give present-day matter and dark-energy densities, transition redshifts, and equation-of-state values close to those inferred from current observations.","The fits to Hubble and Supernovae Ia data in the boundary-coupling chapter indicate that the nonminimal boundary coupling with several potentials can match the observed distance-modulus curve, supporting the models' observational viability."],"supporting_citations":[{"why":"The 1998 Type Ia supernova results that established late-time cosmic acceleration, the observational problem the thesis addresses.","marker":"[1,2]"},{"why":"Supplies the CMB-based values of present density parameters and Hubble constant used for comparison with the models' outcomes.","marker":"[20]"},{"why":"Introduces $f(T,B)$ gravity and its field equations, the theory analysed in chapter 4.","marker":"[37]"},{"why":"Cited for teleparallel Gauss-Bonnet gravity and for scaling solutions that alleviate the coincidence problem.","marker":"[48]"},{"why":"Provides the $f(T,\\phi)$ dynamical-system construction and critical-point comparison used in chapters 2 and 3.","marker":"[60]"},{"why":"Gives the general teleparallel scalar-tensor action on which chapter 2's models are built.","marker":"[63]"},{"why":"The source for the normally-hyperbolic critical-point stability conditions applied to non-hyperbolic fixed points.","marker":"[75]"},{"why":"Supports the power-law $f(T,B)$ form by referring to observational tests with Supernovae Ia data.","marker":"[83]"},{"why":"The observational constraint on the present deceleration parameter used to validate the reported values of $q_0$.","marker":"[90]"},{"why":"The origin of the constant $\\lambda=\\ddot{H}/H^3$ parametrisation that closes the autonomous systems in the boundary-term chapters.","marker":"[93]"}],"fun_headline_variants":["Teleparallel gravity explains cosmic acceleration without Λ","Torsion-based gravity matches dark-energy observations","Modified teleparallel gravity yields stable late-time attractors","No cosmological constant: teleparallel gravity fits cosmic data","Teleparallel gravity: one phase space from matter to dark energy"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the dimensionless quantity $\\lambda = \\ddot{H}/H^3$ stays constant along the cosmological trajectories analysed in the boundary-term chapters, so the dynamical systems are autonomous; if the real Universe's $\\lambda$ varies with time, the critical points and stability classifications derived under this assumption need not describe the actual evolution.","fun_headline_variants_meta":{"raw":{"variants":["Teleparallel gravity explains cosmic acceleration without Λ","Torsion-based gravity matches dark-energy observations","Modified teleparallel gravity yields stable late-time attractors","No cosmological constant: teleparallel gravity fits cosmic data","Teleparallel gravity: one phase space from matter to dark energy"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001049,"raw_usage":{"total_tokens":4435,"prompt_tokens":1000,"completion_tokens":3435,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":616,"completion_tokens_details":{"reasoning_tokens":3360}},"tokens_in":616,"tokens_out":3435,"duration_ms":24004,"temperature":1.0,"reasoning_tokens":3360,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T18:41:03.290159+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute $\\lambda = \\ddot{H}/H^3$ along the best-fit $H(z)$ curves the thesis reports for the $f(T,B)$ and $f(T,T_G)$ models; if $\\lambda$ moves appreciably across the radiation-to-dark-energy epochs instead of holding near the values used to find the critical points, such as $8$ for radiation and $9/2$ for matter, the fixed-point analysis does not govern the models' real trajectories.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The source for the normally-hyperbolic critical-point stability conditions applied to non-hyperbolic fixed points."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supports the power-law $f(T,B)$ form by referring to observational tests with Supernovae Ia data."},{"cited_title":"Capozziello, R","cited_arxiv_id":null,"evidence_quote":"The observational constraint on the present deceleration parameter used to validate the reported values of $q_0$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The origin of the constant $\\lambda=\\ddot{H}/H^3$ parametrisation that closes the autonomous systems in the boundary-term chapters."}],"review_version":1}