Eur. Phys. J. C (2025) 85: 1425
https://doi.org/10.1140/epjc/s10052-025-15151-4

Regular Article - Theoretical Physics

A unified dynamical systems framework for cosmology in f(Q) gravity: generic features across the connection branches

¹ Mathematics Division, Department of Basic Sciences and Social Sciences, North Eastern Hill University, 793022, Shillong, Meghalaya, India
² Visiting Associate, Inter University Centre for Astronomy and Astrophysics, 411 007, Pune, India
³ Department of Mathematics, St. Anthony’s College, 793001, Shillong, Meghalaya, India
⁴ High Energy Physics Theory Group, Department of Physics, Faculty of Science, Chulalongkorn University, 10330, Bangkok, Thailand
⁵ Center for Space Research, North-West University, 2520, Potchefstroom, South Africa
⁶ School of Science, Jiangsu University of Science and Technology, 212100, Zhenjiang, China
⁷ The Institute for Fundamental Study, Naresuan University, 65000, Phitsanulok, Thailand

^a saikat.c@chula.ac.th, saikat.chakraborty@nwu.ac.za

Received: 10 September 2025
Accepted: 3 December 2025
Published online: 14 December 2025

Abstract

We develop a unified dynamical systems framework for spatially flat FLRW cosmology in f(Q) gravity, covering all three connection branches using a single set of Hubble-normalised variables without fixing the function f(Q) a priori. This model-independent and connection-agnostic approach enables direct comparison across connection choices and uncovers structural features of the cosmological dynamics not visible in connection-specific formulations. While the existing works narrowly focus only on fixed point analysis, in our work we put special efforts to identify the invariant submanifolds, model-independent trajectories and physically viable regions of phase space across connections. For a broad class of viable f(Q) models, we establish the generic existence of de Sitter attractors and matter-dominated fixed points in the two branches of connection, offering a robust route to late-time acceleration without fine-tuning. We further identify an invariant submanifold that yields $\mathrm{\Lambda }$CDM-like background evolution despite underlying dynamics distinct from General Relativity, providing a geometric origin for cosmic acceleration distinguishable only at the perturbation level. We also derive a first integral on this submanifold, allowing analytic reconstruction of the dynamical connection and uncovering hidden conservation laws. Another key feature we found is that while trivial connections exhibit strong parameter dependence, the nontrivial branches often feature parameter-independent behaviours. We also study the variation of the effective gravitational coupling, $k_{\text{eff}}$, across branches and show how this can be constrained using astrophysical observations, which bridges theoretical viability with observational consistency in a novel way. Applying our framework to the illustrative model $f(Q)=\alpha Q+\beta {(-Q)}^n$, we find late-time acceleration and $\mathrm{\Lambda }$CDM-like behaviour without vacuum energy. Finally, we propose a general route for extending dynamical systems analysis to broader classes of f(Q) models using the $m_i$-hierarchy method. This framework enables closure for models previously inaccessible to standard approaches and offers a diagnostic tool for identifying structurally viable cosmologies within modified gravity theories.