Past and future isotropization in $f(R)$ dark energy models

We perform the dynamical system analysis for the homogeneous and anisotropic Bianchi I cosmology in the context of Hu-Sawicki $f(R)$ gravity with parameters $\{n,C_{1}\}=\{1,1\}$. Deriving the corresponding four-dimensional dynamical system from the field equations allows us to find the fixed points of the system in both the vacuum and matter cases alongside their stability. This approach also provides analytical expressions for both the cosmic expansion and shear of each fixed point. We show that, for the chosen parameters, Hu-Sawicki gravity allows for isotropization both in the past and in the future for both vacuum and matter. For the matter case, the analysis was subdivided into the dust, radiation and cosmological constant cases, and we show that all three cases exhibit past and future isotropization. The obtained results also revealed interesting aspects of the phase space such as the existence and phenomenology of lines and two-dimensional sheets of fixed points.