Dynamical analysis of coupled curvature-matter scenario in viable f(R) dark energy models at de Sitter phase

We explore the interaction between dark matter and curvature-driven dark energy within viable $f(R)$ gravity models, employing the phase-space analysis approach of linear stability theory. By incorporating an interacting term, denoted as $\mathcal{Q}=\alpha H \tilde{\rho}_{\rm m}\left(\frac{\kappa^2 }{3H^2}\rho_{\rm curv} + 1 \right)$, into the continuity equations of both sectors, we examine dynamics of two $f(R)$ gravity models that adhere to local gravity constraints and fulfill cosmological viability criteria. In de Sitter phase, our investigation reveals modifications of critical points compared to the conventional form, attributed to the introduced interaction. Through a comprehensive phase-space analysis, we illustrate the trajectories near critical points and outline the constraints on the phase space based on cosmological and cosmographic parameters. Furthermore, we can also observe stable behavior of this interacting system during the de Sitter phase of the universe.