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 study the cosmological constant ($\Lambda$) in the standard $\Lambda$CDM model by introducing the \textit{graduated dark energy} (gDE) characterised by a minimal dynamical deviation from the null inertial mass density of the $\Lambda$ in the form $\rho_{\rm inert}\propto \rho^{\lambda}<0$ with $\lambda<1$ being a ratio of two odd integers, for which its energy density $\rho$ dynamically takes negative values in the finite past. For large negative values of $\lambda$, it creates a phenomenological model described by a smooth function that approximately describes the $\Lambda$ spontaneously switching sign in the late universe to become positive today. We confront the model with the latest combined observational data sets of PLK+BAO+SN+$H$. It is striking that the data predict bimodal posterior probability distributions for the parameters of the model along with large negative $\lambda$ values; the new maximum significantly excludes the $\Lambda$, and the old maximum contains the $\Lambda$. The improvement in the goodness of fit for the $\Lambda$ reaches highly significant levels, $\Delta\chi_{\rm min}^2=6.4$ for the new maxima, while it remains at insignificant levels, $\Delta\chi_{\rm min}^2\lesssim0.02$, for the old maxima. We show that, in contrast to the old maxima, which do not distinguish from the $\Lambda$, the new maxima agree with the model-independent $H_0$ measurements, high-precision Ly-$\alpha$ data, and model-independent $Omh^2$ diagnostic estimates. Our results provide strong hints of a spontaneous sign switch in the cosmological constant and lead us to conjecture that the universe has transitioned from AdS vacua to dS vacua, at a redshift $z\approx 2.32$ and triggered the late-time acceleration, and suggests looking for such mechanisms in string theory constructions.

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.

We consider higher-dimensional massive Brans-Dicke theory with Ricci-flat internal space. The background model is perturbed by a massive gravitating source which is pressureless in the external (our space) but has an arbitrary equation-of-state parameter $\Omega$ in the internal space. We obtain the exact solution of the system of linearized equations for the perturbations of the metric coefficients and scalar field. For a massless scalar field, relying on the fine-tuning between the Brans-Dicke parameter $\omega$ and $\Omega$, we demonstrate that (i) the model does not contradict gravitational tests relevant to the parameterized post-Newtonian parameter $\gamma$, and (ii) the scalar field is not ghost in the case of nonzero $|\Omega|\sim O(1)$ along with the natural value $|\omega|\sim O(1)$. In the general case of a massive scalar field, the metric coefficients acquire the Yukawa correction terms, where the Yukawa mass scale $m$ is defined by the mass of the scalar field. For the natural value $\omega\sim O(1)$, the inverse-square-law experiments impose the following restriction on the lower bound of the mass: $m\gtrsim 10^{-11}\,$GeV. The experimental constraints on $\gamma$ requires that $\Omega$ must be extremely close to $-1/2$.

We study a new model of Energy-Momentum Squared Gravity (EMSG), called Energy-Momentum Log Gravity (EMLG), constructed by the addition of the term $f(T_{\mu\nu}T^{\mu\nu})=\alpha \ln(\lambda\,T_{\mu\nu}T^{\mu\nu})$, envisaged as a correction, to the Einstein-Hilbert action with cosmological constant $\Lambda$. The choice of this modification is made as a specific way of including new terms in the right-hand side of the Einstein field equations, resulting in constant effective inertial mass density and, importantly, leading to an explicit exact solution of the matter energy density in terms of redshift. We look for viable cosmologies, in particular, an extension of the standard $\Lambda$CDM model. EMLG provides an effective dynamical dark energy passing below zero at large redshifts, accommodating a mechanism for screening $\Lambda$ in this region, in line with suggestions for alleviating some of the tensions that arise between observational data sets within the standard $\Lambda$CDM model. We present a detailed theoretical investigation of the model and then constrain the free parameter $\alpha'$, a normalisation of $\alpha$, using the latest observational data. The data does not rule out the $\Lambda$CDM limit of our model ($\alpha'= 0$), but prefers slightly negative values of the EMLG model parameter ($\alpha'= -0.032\pm 0.043$), which leads to the screening of $\Lambda$. We also discuss how EMLG relaxes the persistent tension that appears in the measurements of $H_0$ within the standard $\Lambda$CDM model.

We consider the simplest anisotropic generalization, as a correction, to the standard $\Lambda$CDM model, by replacing the spatially flat Robertson-Walker metric by the Bianchi type-I metric, which brings in a new term $\Omega_{\sigma 0}a^{-6}$ (mimicking the stiff fluid) in the average expansion rate $H(a)$ of the Universe. From Hubble and Pantheon data, relevant to the late Universe ($z\lesssim 2.4$), we obtain the constraint $\Omega_{\sigma0}\lesssim10^{-3}$, in line with the model-independent constraints. When the baryonic acoustic oscillations and cosmic microwave background (CMB) data are included, the constraint improves by 12 orders of magnitude, i.e., $\Omega_{\sigma0}\lesssim10^{-15}$. We find that this constraint could alter neither the matter-radiation equality redshift nor the peak of the matter perturbations. Demanding that the expansion anisotropy has no significant effect on the standard big bang nucleosynthesis (BBN), we find the constraint $\Omega_{\sigma0}\lesssim10^{-23}$. We show explicitly that the constraint from BBN renders the expansion anisotropy irrelevant to make a significant change in the CMB quadrupole temperature, whereas the constraint from the cosmological data in our model provides the temperature change up to $\sim11\, \rm mK$, though it is much beyond the CMB quadrupole temperature.

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.

We investigate cosmological perturbations for nonlinear $f(R)$ models within the cosmic screening approach. Matter is considered both in the form of a set of discrete point-like massive bodies and in the form of a continuous pressureless perfect fluid. We perform full relativistic analysis of the first-order theory of scalar perturbations for arbitrary nonlinear $f(R)$ models and demonstrate that scalar potentials $\Phi(t,\mathbf{r})$ and $\Psi(t,\mathbf{r})$ are determined by a system of only two master equations. Our equations are applicable at all spatial scales as long as the approximation $\delta R/\bar R \ll 1$ (which is usually assumed in studies devoted to cosmological perturbations in $f(R)$ models) works.

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 consider 5D brane world models with broken global 4D Poincar\a'{e} invariance (4D part of the spacetime metric is not conformal to the Minkowski spacetime). The bulk is filled with the negative cosmological constant and may contain a perfect fluid. In the case of empty bulk (the perfect fluid is absent), it is shown that one brane solution always has either a physical or a coordinate singularity in the bulk. We cut off these singularities in the case of compact two brane model and obtain regular exact solutions for both 4D Poincar\a'{e} broken and restored invariance. When the perfect fluid is present in the bulk, we get the master equation for the metric coefficients in the case of arbitrary bulk perfect fluid equation of state (EoS) parameters. In two particular cases of EoS, we obtain the analytic solutions for thin and thick branes. First one generalizes the well known Randall-Sundrum model with one brane to the case of the bulk anisotropic perfect fluid. In the second solution, the 4D Poincar\a'{e} invariance is restored. Here, the spacetime goes asymptotically to the anti-de Sitter one far from the thick brane.

We consider Kaluza-Klein models where internal spaces are compact flat or curved Einstein spaces. This background is perturbed by a compact gravitating body with the dust-like equation of state (EoS) in the external/our space and an arbitrary EoS parameter $\Omega$ in the internal space. Without imposing any restrictions on the form of the perturbed metric and the distribution of the perturbed energy densities, we perform the general analysis of the Einstein and conservation equations in the weak-field limit. All conclusions follow from this analysis. For example, we demonstrate that the perturbed model is static and perturbed metric preserves the block-diagonal form. In a particular case $\Omega=-1/2$, the found solution corresponds to the weak-field limit of the black strings/branes. The black strings/branes are compact gravitating objects which have the topology (four-dimensional Schwarzschild spacetime)$\times$ ($d$-dimensional internal space) with $d\geq 1$. We present the arguments in favour of these objects. First, they satisfy the gravitational tests for the parameterized post-Newtonian parameter $\gamma$ at the same level of accuracy as General Relativity. Second, they are preferable from the thermodynamical point of view. Third, averaging over the Universe, they do not destroy the stabilization of the internal space. These are the astrophysical and cosmological aspects of the black strings/branes.

We consider the non-minimal model of gravity in $Y(R) F^2$-form. We investigate a particular case of the model, for which the higher order derivatives are eliminated but the scalar curvature $R$ is kept to be dynamical via the constraint $Y_RF_{mn}F^{mn} =-\frac{2}{\kappa^2}$. The effective fluid obtained can be represented by interacting electromagnetic field and vacuum depending on $Y(R)$, namely, the energy density of the vacuum tracks $R$ while energy density of the conventional electromagnetic field is dynamically scaled with the factor $\frac{Y(R)}{2}$. We give exact solutions for anisotropic inflation by assuming the volume scale factor of the Universe exhibits a power-law expansion. The directional scale factors do not necessarily exhibit power-law expansion, which would give rise to a constant expansion anisotropy, but expand non-trivially and give rise to a non-monotonically evolving expansion anisotropy that eventually converges to a non-zero constant. Relying on this fact, we discuss the anisotropic e-fold during the inflation by considering observed scale invariance in CMB and demanding the Universe to undergo the same amount of e-folds in all directions. We calculate the residual expansion anisotropy at the end of inflation, though as a result of non-monotonic behaviour of expansion anisotropy all the axes of the Universe undergo the same of amount of e-folds by the end of inflation. We also discuss the generation of the modified electromagnetic field during the first few e-folds of the inflation and its persistence against to the vacuum till end of inflation.

We show that if the $\alpha$-attractor model is realized by the spontaneous breaking of the scale symmetry, then the stability and the dynamics of the vector field that gauges the scale symmetry can severely constrain the $\alpha$-parameter as $5/6 < \alpha < 1$ restricting the inflationary predictions in a very tiny region in the $n_s - r$ plane that are in great agreement with the latest Planck data. Although the different values of $\alpha$ do not make a tangible difference for $n_s$ and $r$, they provide radically different scenarios for the post-inflationary dynamics which determines the standard BBN processes and the large scale isotropy of the universe.

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 study the late-time evolution of the Universe where dark energy (DE) is presented by a barotropic fluid on top of cold dark matter (CDM). We also take into account the radiation content of the Universe. Here by the late stage of the evolution we refer to the epoch where CDM is already clustered into inhomogeneously distributed discrete structures (galaxies, groups and clusters of galaxies). Under this condition the mechanical approach is an adequate tool to study the Universe deep inside the cell of uniformity. More precisely, we study scalar perturbations of the FLRW metric due to inhomogeneities of CDM as well as fluctuations of radiation and DE. For an arbitrary equation of state for DE we obtain a system of equations for the scalar perturbations within the mechanical approach. First, in the case of a constant DE equation of state parameter $w$, we demonstrate that our method singles out the cosmological constant as the only viable dark energy candidate. Then, we apply our approach to variable equation of state parameters in the form of three different linear parametrizations of $w$, e.g., the Chevallier-Polarski-Linder perfect fluid model. We conclude that all these models are incompatible with the theory of scalar perturbations in the late Universe.

We introduce a new parametrization for the dark energy, led by the same idea to the linear expansion of the equation of state in scale factor $a$ and in redshift $z$, which diverges neither in the past nor future and contains the same number of degrees of freedom with the former two. We present constraints of the cosmological parameters using the most updated baryon acoustic oscillation (BAO) measurements along with cosmic microwave background (CMB) data and a recent reanalysis of Type Ia supernova (SN) data. This new parametrization allowed us to carry out successive observational analyses by decreasing its degrees of freedom systematically until ending up with a dynamical dark energy model that has the same number of parameters with $\Lambda$CDM. We found that the dark energy source with a dynamical equation of state parameter equal $-2/3$ at the early universe and $-1$ today fits the data slightly better than $\Lambda$.

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.

The parametrizations $q=q_0+q_1 z$ and $q=q_0+q_1 (1-a/a_0)$ (Chevallier-Polarski-Linder parametrization) of deceleration parameter, which are linear in cosmic redshift $z$ and scale factor $a$, have been frequently utilized in the literature to study kinematics of Universe. In this paper, we follow a strategy that leads to these two well known parametrizations of deceleration parameter as well as an additional new parametrization $q=q_0+q_1(1-t/t_0)$, which is linear in cosmic time $t$. We study the features of this linearly time-varying deceleration parameter in contrast with the other two linear parametrizations. We investigate in detail the kinematics of the Universe by confronting the three models with the latest observational data. We further study the dynamics of the Universe by considering the linearly time-varying deceleration parameter model in comparison with the standard $\Lambda$CDM model. We also discuss future of the Universe in the context of the models under consideration.

The field equations of Brans-Dicke gravity coupled to a mass-varying vector field are derived. Anisotropic cosmological solutions with a locally rotationally symmetric Bianchi type I metric and time-dependent scalar and electric vector fields are studied. A particular class of exact solutions for which all the variable parameters have a power-law time dependence is given. The universe expands with a constant expansion anisotropy within this class of solutions. We show that the accelerating expansion is driven by the scalar field and the electric vector field can be interpreted as an anisotropic dark matter source.

In this paper, we consider a simple form of expansion history of Universe referred to as the hybrid expansion law - a product of power-law and exponential type of functions. The ansatz by construction mimics the power-law and de Sitter cosmologies as special cases but also provides an elegant description of the transition from deceleration to cosmic acceleration. We point out the Brans-Dicke realization of the cosmic history under consideration. We construct potentials for quintessence, phantom and tachyon fields, which can give rise to the hybrid expansion law in general relativity. We investigate observational constraints on the model with hybrid expansion law applied to late time acceleration as well as to early universe a la nucleosynthesis.

A class of cosmological solutions of higher dimensional Einstein field equations with the energy-momentum tensor of a homogeneous, isotropic fluid as the source are considered with an anisotropic metric that includes the direct sum of a 3-dimensional (physical, flat) external space metric and an n-dimensional (compact, flat) internal space metric. A simple kinematical constraint is postulated that correlates the expansion rates of the external and internal spaces in terms of a real parameter {\lambda}. A specific solution for which both the external and internal spaces expand at different rates is given analytically for n=3. Assuming that the internal dimensions were at Planck length scales when the external space starts with a Big Bang (t=0), they expand only 1.49 times and stay at Planck length scales even in the present age of the universe (13.7 Gyr). The effective four dimensional universe would exhibit a behavior consistent with our current understanding of the observed universe. It would start in a stiff fluid dominated phase and evolve through radiation dominated and pressureless matter dominated phases, eventually going into a de Sitter phase at late times.

We present cosmological solutions for (1+3+n)-dimensional steady state universe in dilaton gravity with an arbitrary dilaton coupling constant w and exponential dilaton self-interaction potentials in the string frame. We focus particularly on the class in which the 3-space expands with a time varying deceleration parameter. We discuss the number of the internal dimensions and the value of the dilaton coupling constant to determine the cases that are consistent with the observed universe and the primordial nucleosynthesis. The 3-space starts with a decelerated expansion rate and evolves into accelerated expansion phase subject to the values of w and n, but ends with a Big Rip in all cases. We discuss the cosmological evolution in further detail for the cases w=1 and w=1/2 that permit exact solutions. We also comment on how the universe would be conceived by an observer in four dimensions who is unaware of the internal dimensions and thinks that the conventional general relativity is valid at cosmological scales.

We compare the cosmological kinematics obtained via our law of linearly varying deceleration parameter (LVDP) with the kinematics obtained in the {\Lambda}CDM model. We show that the LVDP model is almost indistinguishable from the {\Lambda}CDM model up to the near future of our universe as far as the current observations are concerned, though their predictions differ tremendously into the far future.