We present a detailed investigation of the Rastall gravity extension of the standard $\Lambda$CDM model. We review the model for two simultaneous modifications of different nature in the Friedmann equation due to the Rastall gravity: the new contributions of the material (actual) sources (considered as effective source) and the altered evolution of the material sources. We discuss the role/behavior of these modifications with regard to some low redshift tensions, including the $H_0$ tension, prevailing within the standard $\Lambda$CDM. We constrain the model at the level of linear perturbations, and obtain the first constraints through a robust and accurate analysis using the latest full Planck CMB data, with and without including BAO data. We find that the Rastall parameter $\epsilon$ (null for general relativity) is consistent with zero at 68% CL (with a tendency towards positive values, $-0.0001 < \epsilon < 0.0007$ (CMB+BAO) at 68% CL), which in turn implies no significant statistical evidence for deviation from general relativity, and also a precision of $\mathcal{O}(10^{-4})$ for the coefficient $-1/2$ of the term $g_{\mu\nu}R$ in the Einstein field equations of general relativity (guaranteeing the local energy-momentum conservation). We explore the consequences led by the Rastall gravity on the cosmological parameters in the light of the observational analyses. It turns out that the effective source dynamically screens the usual vacuum energy at high redshift, but this mechanism barely works due to the opposition by the altered evolution of CDM. Consequently, two simultaneous modifications of different nature in the Friedmann equation act against each other, and do not help to considerably relax the so called low redshift tensions. Our results may offer a guide for the research community that studies the Rastall gravity in various aspects of gravitation and cosmology.

We present an explicit detailed theoretical and observational investigation of an anisotropic massive Brans-Dicke (BD) gravity extension of the standard $\Lambda$CDM model, wherein the extension is characterized by two additional degrees of freedom; the BD parameter, $\omega$, and the present day density parameter corresponding to the shear scalar, $\Omega_{\sigma^2,0}$. The BD parameter, determining the deviation from general relativity (GR), by alone characterizes both the dynamics of the effective dark energy (DE) and the redshift dependence of the shear scalar. These two affect each other depending on $\omega$, namely, the shear scalar contributes to the dynamics of the effective DE, and its anisotropic stress --which does not exist in scalar field models of DE within GR-- controls the dynamics of the shear scalar deviating from the usual $\propto(1+z)^6$ form in GR. We mainly confine the current work to non-negative $\omega$ values as it is the right sign --theoretically and observationally-- for investigating the model as a correction to the $\Lambda$CDM. By considering the current cosmological observations, we find that $\omega\gtrsim 250$, $\Omega_{\sigma^2,0}\lesssim 10^{-23}$ and the contribution of the anisotropy of the effective DE to this value is insignificant. We conclude that the simplest anisotropic massive BD gravity extension of the standard $\Lambda$CDM model exhibits no significant deviations from it all the way to the Big Bang Nucleosynthesis. We also point out the interesting features of the model in the case of negative $\omega$ values; for instance, the constraints on $\Omega_{\sigma^2,0}$ could be relaxed considerably, the values of $\omega\sim-1$ (relevant to string theories) predict dramatically different dynamics for the expansion anisotropy.

In this paper, we introduce a scale-independent energy-momentum squared gravity (EMSG) that allows different gravitational couplings for different types of sources, which may lead to scenarios with many interesting applications/implications in cosmology. In the present study, to begin with, we study a modification of the $\Lambda$ cold dark matter ($\Lambda$CDM) model, where photons and baryons couple to the spacetime as in general relativity, while the cold dark matter and relativistic relics (neutrinos and any other relativistic relics) couple to the spacetime in accordance with EMSG. This scenario induces pseudo nonminimal interactions on these components, leading to modification at both the background and perturbative levels. A consequence of this scenario is that the dimensionless free parameter of the theory may induce direct changes on the effective number of the relativistic species, without the need to introduce new extra species. In order to quantify the observational consequences of the cosmological scenario, we use the cosmic microwave background Planck data (temperature, polarization, and lensing power spectrum) and baryonic acoustic oscillations data. We find that the free model parameter is too small to induce statistically significant corrections on the $\Lambda$CDM model due to EMSG. We deduce that the model presented here is quite rich with promising cosmological applications/implications that deserve further investigations.

Deviations from the predictions of general relativity due to energy-momentum squared gravity (EMSG) are expected to become pronounced in the high density cores of neutron stars. We derive the hydrostatic equilibrium equations in EMSG and solve them numerically to obtain the neutron star mass-radius relations for four different realistic equations of state. We use the existing observational measurements of the masses and radii of neutron stars to constrain the free parameter, $\alpha ,$ that characterizes the coupling between matter and spacetime in EMSG. We show that $-10^{-38}\,\mathrm{cm^{3}/erg}<\alpha <+10^{-37}\,\mathrm{cm^{3}/erg}$. Under this constraint, we discuss what contributions EMSG can provide to the physics of neutron stars, in particular, their relevance to the so called \textit{hyperon puzzle} in neutron stars. We also discuss how EMSG alters the dynamics of the early universe from the predictions of the standard cosmological model. We show that EMSG leaves the standard cosmology safely unaltered back to $t\sim 10^{-4}$ seconds at which the energy density of the universe is $\sim 10^{34}\,\mathrm{erg\,cm^{-3}}$.

We propose a modified theory of gravitation constructed by the addition of the term $f(T_{\mu\nu}T^{\mu\nu})$ to the Einstein-Hilbert action, and elaborate a particular case $f(T_{\mu\nu}T^{\mu\nu})=\alpha(T_{\mu\nu}T^{\mu\nu})^{\eta}$, where $\alpha$ and $\eta$ are real constants, dubbed as energy-momentum powered gravity (EMPG). We search for viable cosmologies arising from EMPG especially in the context of the late-time accelerated expansion of the Universe. We investigate the ranges of the EMPG parameters $(\alpha,\eta)$ on theoretical as well as observational grounds leading to the late-time acceleration of the Universe with pressureless matter only, while keeping the successes of standard general relativity at early times. We find that $\eta=0$ corresponds to the $\Lambda$CDM model, whereas $\eta\neq 0$ leads to a $w$CDM-type model. However, the underlying physics of the EMPG model is entirely different in the sense that the energy in the EMPG Universe is sourced by pressureless matter only. Moreover, the energy of the pressureless matter is not conserved, namely, in general it does not dilute as $\rho\propto a^{-3}$ with the expansion of the Universe. Finally, we constrain the parameters of an EMPG-based cosmology with a recent compilation of 28 Hubble parameter measurements, and find that this model describes an evolution of the Universe similar to that in the $\Lambda$CDM model. We briefly discuss that EMPG can be unified with Starobinsky gravity to describe the complete history of the Universe including the inflationary era.

We present a mini review of the Stueckelberg mechanism, which was proposed to make the abelian gauge theories massive as an alternative to Higgs mechanism, within the framework of Minkowski as well as curved spacetimes. The higher the scale the tighter the bounds on the photon mass, which might be gained via the Stueckelberg mechanism, may be signalling that even an extremely small mass of the photon which cannot be measured directly could have far reaching effects in cosmology. We present a cosmological model where Stueckelberg fields, which consist of both scalar and vector fields, are non-minimally coupled to gravity and the universe could go through a decelerating expansion phase sandwiched by two different accelerated expansion phases. We discuss also the possible anisotropic extensions of the model.

We investigate a cosmological model in which the Stueckelberg fields are non-minimally coupled to the scalar curvature in a gauge invariant manner. We present not only a solution that can be considered in the context of the late time acceleration of the universe but also a solution compatible with the inflationary cosmology. Distinct behaviors of the scalar and vector fields together with the real valued mass gained by the Stueckelberg mechanism lead the universe to go through the two different accelerated expansion phases with a decelerated expansion phase between them. On the other hand, in the solutions we present, if the mass is null then the universe is either static or exhibits a simple power law expansion due to the vector field potential.

We consider the single field chaotic m^2\phi^2 inflationary model with a period of preheating, where the inflaton decays to another scalar field \chi\ in the parametric resonance regime. In a recent work, one of us has shown that the \chi\ modes circulating in the loops during preheating notably modify the <\zeta\zeta> correlation function. We first rederive this result using a different gauge condition hence reconfirm that superhorizon \zeta\ modes are affected by the loops in preheating. Further, we examine how \chi\ loops give rise to non-gaussianity and affect the tensor perturbations. For that, all cubic and some higher order interactions involving two \chi\ fields are determined and their contribution to the non-gaussianity parameter f_{NL} and the tensor power spectrum are calculated at one loop. Our estimates for these corrections show that while a large amount of non-gaussianity can be produced during reheating, the tensor power spectrum receive moderate corrections. We observe that the loop quantum effects increase with more \chi\ fields circulating in the loops indicating that the perturbation theory might be broken down. These findings demonstrate that the loop corrections during reheating are significant and they must be taken into account for precision inflationary cosmology.

When Brans-Dicke Theory is formulated in terms of the Jordan scalar field \phi, dark energy is related to the mass of this field. We show that if \phi is taken to be a complex scalar field then an exact solution of the vacuum equations shows that Friedmann equation possesses a term, proportional to the inverse sixth power of the scale factor, as well as a constant term. Possible interpretations and phenomenological implications of this result are discussed.