Constraints on Horndeski Gravity with Phantom Crossing

Gravity models in which the dark energy equation of state crosses $w=-1$, also known as the phantom divide, have received extensive interest due to recent analyses favouring this behaviour. We introduce a new subclass of Horndeski scalar-tensor models capable of generating phantom crossing, whilst remaining minimally coupled to matter: the Asymptotic Cubic Galileon (ACG) models. We show that ACG models can jointly fit the expansion history inferred from observations of the Planck cosmic microwave background, baryon acoustic oscillation measurements from the Dark Energy Spectroscopic Instrument, and distance-ladder supernovae measurements from the Dark Energy Survey. We then demonstrate that perturbative observables, including the galaxy-ISW cross-correlation and void force profile, provide powerful constraints that confine viable and testable ACG models to a well-defined region of the broader Horndeski landscape. Model comparison metrics, including $χ^{2}$ and Bayesian evidence, favour both ACG and $w_{0}w_{a}$CDM models over $Λ$CDM, with ACG providing a fit of comparable quality to $w_{0}w_{a}$CDM. Crucially, ACG models ground the observationally preferred $w_{0}w_{a}$CDM behaviour in a robust Lagrangian formulation. This enables interpretation beyond mere phenomenological fits, and motivates further tests of these models on nonlinear scales.