Quadratic energy-momentum squared gravity: constraints from big bang nucleosynthesis

In this work, we extend the standard cosmological model within the quadratic energy-momentum squared gravity (qEMSG) framework, introducing a nonminimal interaction between the usual material field ($T_{\mu\nu}$) and its accompanying partner field (qEMSF, $T_{\mu\nu}^{\rm qEMSF}$), defined by $f(\mathbf{T}^2)=\alpha\mathbf{T}^2$ with $\mathbf{T^2}=T_{\mu\nu}T^{\mu\nu}$. Adopting an analytical approach within the qEMSG framework, we present a comprehensive exploration of Big Bang Nucleosynthesis (BBN) dynamics. Our analysis selects the radiation-dominated universe solution compatible with the standard cosmological model limit as $\alpha\rightarrow0$ and reveals that qEMSF interaction model can modify the radiation energy density's evolution, potentially altering neutron-proton interconversion rates and consequently affecting $^4$He abundance in various ways. By explicitly defining modifications to the predicted primordial $^4$He mass fraction, $Y_{\rm p}$, we establish the most stringent cosmological constraints on the parameter $\alpha$ based on recent measurements of $Y_{\rm p}$: $(-8.81\leq\alpha\leq8.14)\times10^{-27}\,\mathrm{eV}^{-4}$ (68% CL) from Aver et al.'s primordial $^4$He abundance measurements, aligning with $\alpha=0$. Additionally, $(3.48\leq\alpha\leq4.43)\,\times 10^{-27}\rm{eV}^{-4}$ (68% CL) from Fields et al.'s estimates, utilizing the Planck-CMB estimated baryon density within the standard cosmological model framework, diverges from $\alpha=0$, thereby lending support to the qEMSF interaction model. The study also highlights the bidirectional nature of energy-momentum/entropy transfer in qEMSF interaction model, depending on the sign of $\alpha$. The implications of qEMSF in the presence of additional relativistic relics are also explored, showcasing the model's potential to accommodate deviations from standard cosmology and the Standard Model of particle physics.