Based on tensions between the early and late time cosmology, we proposed a double valued cosmological constant which could undergo a phase transition in its history. It is named "double-$\Lambda$ Cold Dark Matter": $\Lambda\Lambda$CDM. An occurred phase transition results in (micro-) structures for the dark sector with a proper (local) interaction. In this paper, inspired by the physics of critical phenomena, we study a simplified model such that the cosmological constant has two values before a transition scale factor, $a_t$, and afterwards it becomes single-valued. We consider both the background and perturbation data sets including CMB, BAO distances and R19 data point. $\Lambda\Lambda$CDM has its maximum likelihood for $a_t= 0.916^{+0.055}_{-0.0076}$ and $H_0= 72.8\pm 1.6$. This result shows no inconsistency between early and late time measurements of Hubble parameter in $\Lambda\Lambda$CDM model. In comparison to $\Lambda$CDM, our model has better fit to data such that $\Delta\chi^2=-11$ and even if we take care of two additional degrees of freedom we do have better AIC quantity $\Delta$AIC$=-7$. We conclude that a phase transition in the behavior of dark energy can address $H_0$ tension successfully and may be responsible for the other cosmological tensions.

A dark energy model (DE) is proposed based on Ginzburg-Landau theory of phase transition (GLT). This model, GLTofDE, surprisingly provides a framework to study not only temporal tensions in cosmology e.g. $H_0$ tension but also spatial anomalies of CMB e.g. the hemispherical power asymmetry and quadrupole-octopole alignment. In the mean field (or Landau) approximation of GLTofDE, there is a spontaneously symmetry breaking exactly like the Higgs potential. We modeled this transition, phenomenologically, and showed that GLTofDE can resolve both the $H_0$ tension and Lyman-$\alpha$ anomaly in a non-trivial way. According to $\chi^2$-analysis the transition happens at $z_t=0.738\pm0.028$ while $H_0=71.89\pm0.93$ km/s/Mpc and $\Omega^{like}_{k}=-0.225\pm0.049$ which are consistent with the latest $H(z)$ reconstructions. In addition, the GLTofDE proposes a framework to address the CMB anomalies when it is considered beyond the mean field approximation. In this regime existence of a long wavelength mode is a typical consequence which is named the Goldstone mode in the case of continuous symmetries. This mode, which is an automatic byproduct in GLTofDE, makes cosmological constant, direction dependent. This means one side of the sky should be colder than the other side in agreement with what has been already observed in CMB. In addition between initial stochastic pattern and the final state with one long wavelength mode, we can observe smaller patches or protrusions of the biggest remaining patch in the simulation. Our simulations show these protrusions are few in numbers and will be evolved according to Alan-Cahn mechanism. These protrusions can give an additional effect on CMB which is the existence of aligned quadrupole-octopole mode and its direction should be orthogonal to the dipole direction. We conclude that GLTofDE is a fertile framework both theoretically and phenomenologically.

We propose a cosmological model, \"u$\Lambda$CDM, based on {\it \"uber-gravity}, which is a canonical ensemble average of many theories of gravity. In this model, we have a sharp transition from (a purely) $\Lambda$CDM era to a phase in which the Ricci scalar is a constant. This transition occurs when the Ricci scalar reaches a critical scale or alternatively at a transition redshift, $z_\oplus$. We use the observations of baryonic acoustic oscillations (BAO) and Supernovae Ia (SNe), as well as the cosmic microwave background (CMB) data to constrain \"u$\Lambda$CDM. This yields $H_0=70.6_{-1.3}^{+1.1}$ km/s/Mpc, $\Omega_m=0.2861 \pm{0.0092} $ and $z_\oplus=0.537_{-0.375}^{+0.277}$, providing a marginally better fit with a Akaike information criterion of 0.8. %After adding in the data from BAO in the Lyman-$\alpha$ forest in quasar spectra, we get $\Delta\chi^2=-2.14$, where the best fit occurs at $H_0=71.2_{-1.5}^{+1.8}$ km/s/Mpc, $\Omega_m=0.280\pm 0.013$ and $z_{\oplus}= 0.41^{+0.31}_{-0.29}$. Therefore, \"u$\Lambda$CDM can ease the $H_0$-tension, albeit marginally, with one additional free parameter. We also provide a preliminary study of the linear perturbation theory in \"u$\Lambda$CDM which points to interesting potential {\it smoking guns} in the observations of large scales structure at $z < z_{\oplus}$.

Recently, the idea of taking ensemble average over gravity models has been introduced. Based on this idea, we study the ensemble average over (effectively) all the gravity models (constructed from Ricci scalar) dubbing the name \"uber-gravity which is a {\it{fixed point}} in the model space. The \"uber-gravity has interesting universal properties, independent from the choice of basis: $i)$ it mimics Einstein-Hilbert gravity for high-curvature regime, $ii)$ it predicts stronger gravitational force for an intermediate-curvature regime, $iii)$ surprisingly, for low-curvature regime, i.e. $R<R_0$ where $R$ is Ricci scalar and $R_0$ is a given scale, the Lagrangian vanishes automatically and $iiii)$ there is a sharp transition between low- and intermediate-curvature regimes at $R=R_0$. We show that the \"uber-gravity response is robust to all values of vacuum energy, $\rho_{vac}$ when there is no other matter. So as a toy model, \"uber-gravity, gives a way to think about the hierarchy problems e.g. the cosmological constant problem. Due to the transition at $R=R_0$ there is a chance for \"uber-gravity to bypass Weinberg's no-go theorem. The cosmology of this model is also promising because of its non-trivial predictions for small curvature scales in comparison to $\Lambda$CDM model.

We compute a general expression for the contribution of vector perturbations to the redshift-space distortion of galaxy surveys. We show that they contribute to the same multipoles of the correlation function as scalar perturbations and should thus in principle be taken into account in data analysis. We derive constraints for next-generation surveys on the amplitude of two sources of vector perturbations, namely non-linear clustering and topological defects. While topological defects leave a very small imprint on redshift-space distortions, we show that the multipoles of the correlation function are sensitive to vorticity induced by non-linear clustering. Therefore future redshift surveys such as DESI or the SKA should be capable of measuring such vector modes, especially with the hexadecapole which appears to be the most sensitive to the presence of vorticity.

We put forward the idea that all the theoretically consistent models of gravity have contributions to the observed gravity interaction. In this formulation, each model comes with its own Euclidean path-integral weight where general relativity (GR) has automatically the maximum weight in high-curvature regions. We employ this idea in the framework of Lovelock models and show that in four dimensions the result is a specific form of the $f(R,G)$ model. This specific $f(R,G)$ satisfies the stability conditions and possesses self-accelerating solutions. Our model is consistent with the local tests of gravity since its behavior is the same as in GR for the high-curvature regime. In the low-curvature regime the gravitational force is weaker than in GR, which can be interpreted as the existence of a repulsive fifth force for very large scales. Interestingly, there is an intermediate-curvature regime where the gravitational force is stronger in our model compared to GR. The different behavior of our model in comparison with GR in both low- and intermediate-curvature regimes makes it observationally distinguishable from $\Lambda$CDM.

We discuss a scale-free model of bigravity, in which the mass parameter of the standard bigravity potential is promoted to a dynamical scalar field. This modification retains the ghost-free bigravity structure, in particular it remains free of the Boulware-Deser ghost. We investigate the theory's interaction structure, focusing on its consistent scaling limits and strong coupling scales. Furthermore we explore the model's quadratic action, both around generic background configurations and paying special attention to cosmological backgrounds and to the associated background evolution. Finally we consider the possibility of realizing a phase of late-time acceleration as well as a quasi-de Sitter inflationary stage at early times, when the promoted "mass scalar" becomes the inflaton.

Weyl bi-connection model manifests a natural framework to automatically produce the Galileon structure. It is shown that this framework can explain scalar Galileon, vector Galileon as well as their interactions by generalizing the Weyl non-metricity. So it can be interpreted as a geometrical realization for Galileons. The non-metricity part enjoys a $U(1)$ gauge invariance. The result is interestingly non-trivial since the Galileon structure appears spontaneously and not by demanding the absence of the Ostrogradsky ghost. This fact suggests a possible deeper conceptual relation between Weyl bi-connection model and the absence of Ostrogradsky ghost.

We study the cosmological consequences of a recently proposed nonlocal modification of general relativity, obtained by adding a term $m^2R\,\Box^{-2}R$ to the Einstein-Hilbert action. The model has the same number of parameters as $\Lambda$CDM, with $m$ replacing $\Omega_{\Lambda}$, and is very predictive. At the background level, after fixing $m$ so as to reproduce the observed value of $\Omega_M$, we get a pure prediction for the equation of state of dark energy as a function of redshift, $w_{\rm DE}(z)$, with $w_{\rm DE}(0)$ in the range $[-1.165,-1.135]$ as $\Omega_M$ varies over the broad range $\Omega_M\in [0.20,0.36]$. We find that the cosmological perturbations are well-behaved, and the model fully fixes the dark energy perturbations as a function of redshift $z$ and wavenumber $k$. The nonlocal model provides a good fit to supernova data and predicts deviations from General Relativity in structure formation and in weak lensing at the level of 3-4%, therefore consistent with existing data but readily detectable by future surveys. For the logarithmic growth factor we obtain $\gamma\simeq 0.53$, to be compared with $\gamma\simeq 0.55$ in $\Lambda$CDM. For the Newtonian potential on subhorizon scales our results are well fitted by $\Psi(a;k)=[1+\mu_s a^s]\Psi_{\rm GR}(a;k)$ with a scale-independent $\mu_s\simeq 0.09$ and $s\simeq 2$, while the anisotropic stress is negligibly small.

Gravitational vector degrees of freedom typically arise in many examples of modified gravity models. We start to systematically explore their role in these scenarios, studying the effects of coupling gravitational vector and scalar degrees of freedom. We focus on set-ups that enjoy a Galilean symmetry in the scalar sector and an Abelian gauge symmetry in the vector sector. These symmetries, together with the requirement that the equations of motion contain at most two space-time derivatives, only allow for a small number of operators in the Lagrangian for the gravitational fields. We investigate the role of gravitational vector fields for two broad classes of phenomena that characterize modified gravity scenarios. The first is self-acceleration: we analyze in general terms the behavior of vector fluctuations around self-accelerating solutions, and show that vanishing kinetic terms of vector fluctuations lead to instabilities on cosmological backgrounds. The second phenomenon is the screening of long range fifth forces by means of Vainshtein mechanism. We show that if gravitational vector fields are appropriately coupled to a spherically symmetric source, they can play an important role for defining the features of the background solution and the scale of the Vainshtein radius. Our general results can be applied to any concrete model of modified gravity, whose low-energy vector and scalar degrees of freedom satisfy the symmetry requirements that we impose.

What is the right way to interpret a massive graviton? We generalize the kinematical framework of general relativity to multiple connections. The average of the connections is itself a connection and plays the role of the canonical connection in standard General Relativity. At the level of dynamics, the simplest choice of the Einstein-Hilbert action is indistinguishable from the single-connection case. However, inspired by Weyl geometry, we show how one can construct massive gravity to all orders in perturbation theory compatible with the de Rham-Gabadadze-Tolley ghost-free model. We conclude that the mass of the graviton can be interpreted as a geometrical property of spacetime arising from two connections. Furthermore in the multi-connection framework there is no ambiguity in the definition of physical metric and consequently coupling to matter.

We study linear perturbations around time dependent spherically symmetric solutions in the Lambda_3 massive gravity theory, which self-accelerate in the vacuum. We find that the dynamics of the scalar perturbations depend on the coordinate choice for the background solutions. For particular choices of coordinates there is a symmetry enhancement, leaving no propagating scalar degrees of freedom at linear order in perturbations. In contrast, any other coordinate choice propagates a single scalar mode. We find that the Hamiltonian of this scalar mode is unbounded from below for all self-accelerating solutions, signalling an instability.

In an open Friedmann-Robertson-Walker (FRW) space background, we study the classical and quantum cosmological models in the framework of the recently proposed nonlinear massive gravity theory. Although the constraints which are present in this theory prevent it from admitting the flat and closed FRW models as its cosmological solutions, for the open FRW universe, it is not the case. We have shown that, either in the absence of matter or in the presence of a perfect fluid, the classical field equations of such a theory adopt physical solutions for the open FRW model, in which the mass term shows itself as a cosmological constant. These classical solutions consist of two distinguishable branches: One is a contacting universe which tends to a future singularity with zero size, while another is an expanding universe having a past singularity from which it begins its evolution. A classically forbidden region separates these two branches from each other. We then employ the familiar canonical quantization procedure in the given cosmological setting to find the cosmological wave functions. We use the resulting wave function to investigate the possibility of the avoidance of classical singularities due to quantum effects. It is shown that the quantum expectation values of the scale factor, although they have either contracting or expanding phases like their classical counterparts, are not disconnected from each other. Indeed, the classically forbidden region may be replaced by a bouncing period in which the scale factor bounces from the contraction to its expansion eras. Using the Bohmian approach of quantum mechanics, we also compute the Bohmian trajectory and the quantum potential related to the system, which their analysis shows are the direct effects of the mass term on the dynamics of the universe.

We employ the effective field theory approach for multi-field inflation which is a generalization of Weinberg's work. In this method the first correction terms in addition to standard terms in the Lagrangian have been considered. These terms contain up to the fourth derivative of the fields including the scalar field and the metric. The results show the possible shapes of the interaction terms resulting eventually in non-Gaussianity in a general formalism. In addition generally the speed of sound is different but almost unity. Since in this method the adiabatic mode is not discriminated initially so we define the adiabatic as well as entropy modes for a specific two-field model. It has been shown that the non-Gaussianity of the adiabatic mode and the entropy mode are correlated in shape and amplitude. It is shown that even for speed close to unity large non-Gaussianities are possible in multi-field case. The amount of the non-Gaussianity depends on the curvature of the classical path in the phase-space in the Hubble unit such that it is large for the large curvature. In addition it is emphasized that the time derivative of adiabatic and entropy perturbations do not transform due to the shift symmetry as well as the original perturbations. Though two specific combinations of them are invariant under such a symmetry and these combinations should be employed to construct an effective field theory of multi-field inflation.

We investigate the effects of a special kind of dynamical deformation between the momenta of the scalar field of the Brans-Dicke theory and the scale factor of the FRW metric. This special choice of deformation includes linearly a deformation parameter. We trace the deformation footprints in the cosmological equations of motion when the BD coupling parameter goes to infinity. One class of the solutions gives a constant scale factor in the late time that confirms the previous result obtained via another approach in the literature. This effect can be interpreted as a quantum gravity footprint in the coarse grained explanation. The another class of the solutions removes the big bang singularity, and the accelerating expansion region has an infinite temporal range which overcomes the horizon problem. After this epoch, there is a graceful exiting by which the universe enters in the radiation dominated era.

We study the effects of noncommutativity, in the form of a Lie-algebraically deformed Poisson commutation relations, on the evolution of a Bianchi type I cosmological model with a positive cosmological constant. The phase space variables turn out to correspond to the scale factors of this model in $x$, $y$ and $z$ directions. According to the conditions that the structure constants (deformation parameters) should satisfy, we argue that there are two types of noncommutative phase space with Lie-algebraic structure. The exact classical solutions in commutative and type I noncommutative cases are presented. In the framework of this type of deformed phase space, we investigate the possibility of building a Bianchi I model with cyclic scale factors in which the size of the universe in each direction experiences an endless sequence of contractions and re-expansions. We also obtain some approximate solutions for the type II noncommutative structure by numerical methods and show that the cyclic behavior is repeated as well. These results are compared with the standard commutative case, and similarities and differences of these solutions are discussed.

We investigate brane inflation driven by two stacks of mobile branes in a throat. The stack closest to the bottom of the throat annihilates first with antibranes, resulting in particle production and a change of the equation of state parameter w. We calculate analytically some observable signatures of the collision; related decays are common in multi-field inflation, providing the motivation for this case study. The discontinuity in w enters the matching conditions relating perturbations in the remaining degree of freedom before and after the collision, affecting the power-spectrum of curvature perturbations. We find an oscillatory modulation of the power-spectrum for scales within the horizon at the time of the collision, and a slightly redder spectrum on super-horizon scales. We comment on implications for staggered inflation.

We consider a class of generalized FRW type metrics in the context of higher dimensional Einstein gravity in which the extra dimensions are allowed to have different scale factors. It is shown that noncommutativity between the momenta conjugate to the internal space scale factors controls the power-law behavior of the scale factors in the extra dimensions, taming it to an oscillatory behavior. Hence noncommutativity among the internal momenta of the mini-super-\emph{phase}-space can be used to explain stabilization of the compactification volume of the internal space in a higher dimensional gravity theory.