Phantom dark energy can produce amplified cosmic acceleration at late times, thus increasing the value of $H_0$ favored by CMB data and releasing the tension with local measurements of $H_0$. We show that the best fit value of $H_0$ in the context of the CMB power spectrum is degenerate with a constant equation of state parameter $w$, in accordance with the approximate effective linear equation $H_0 + 30.93\; w - 36.47 = 0$ ($H_0$ in $km \; sec^{-1} \; Mpc^{-1}$). This equation is derived by assuming that both $\Omega_{0 \rm m}h^2$ and $d_A=\int_0^{z_{rec}}\frac{dz}{H(z)}$ remain constant (for invariant CMB spectrum) and equal to their best fit Planck/$\Lambda$CDM values as $H_0$, $\Omega_{0 \rm m}$ and $w$ vary. For $w=-1$, this linear degeneracy equation leads to the best fit $H_0=67.4 \; km \; sec^{-1} \; Mpc^{-1}$ as expected. For $w=-1.22$ the corresponding predicted CMB best fit Hubble constant is $H_0=74 \; km \; sec^{-1} \; Mpc^{-1}$ which is identical with the value obtained by local distance ladder measurements while the best fit matter density parameter is predicted to decrease since $\Omega_{0 \rm m}h^2$ is fixed. We verify the above $H_0-w$ degeneracy equation by fitting a $w$CDM model with fixed values of $w$ to the Planck TT spectrum showing also that the quality of fit ($\chi^2$) is similar to that of $\Lambda$CDM. However, when including SnIa, BAO or growth data the quality of fit becomes worse than $\Lambda$CDM when $w< -1$. Finally, we generalize the $H_0-w(z)$ degeneracy equation for $w(z)=w_0+w_1\; z/(1+z)$ and identify analytically the full $w_0-w_1$ parameter region that leads to a best fit $H_0=74\; km \; sec^{-1} \; Mpc^{-1}$ in the context of the Planck CMB spectrum. This exploitation of $H_0-w(z)$ degeneracy can lead to immediate identification of all parameter values of a given $w(z)$ parametrization that can potentially resolve the $H_0$ tension.

A redshift tomography of the Pantheon type Ia supernovae (SnIa) data focusing on the best fit value of the absolute magnitude $M$ and/or Hubble constant $H_0$ in the context of $\Lambda$CDM indicates a local variation ($z\lesssim 0.2$) at $2\sigma$ level, with respect to the best fit of the full dataset. If this variation is physical, it can be interpreted either as a locally higher value of $H_0$, corresponding to a local matter underdensity $\delta \rho_0/\rho_0 \simeq -0.10 \pm 0.04$ or as a time variation of Newton's constant which implies an evolving Chandrasekhar mass and thus an evolving absolute magnitude $M$ of SnIa. The local void scenario would predict an anisotropy in the best fit value of $H_0$ since it is unlikely that we are located at the center of a local spherical underdensity. Using a hemisphere comparison method we find an anisotropy level consistent with simulated isotropic datasets. We show however, that the anisotropic sky distribution of the Pantheon SnIa data induces a preferred range of directions even in simulated Pantheon data obtained in the context of isotropic $\Lambda$CDM. We thus construct a more isotropically distributed subset of the Pantheon SnIa and show that the preferred range of directions disappears. Using this subset we again find no evidence for anisotropy using either the hemisphere comparison method or the dipole fit method. In the context of the modified gravity scenario, we allow for an evolving normalized Newton's constant consistent with General Relativity (GR) at early and late times $\mu(z)=1+g_a z^2/(1+z)^2-g_a z^4/(1+z)^4$ and fit for $g_a$ assuming $L\sim G_{\rm{eff}}^b$. For $b=-3/2$ indicated by previous studies we find $g_a=-0.47 \pm 0.36$ which is more than $1.5\sigma$ away from the GR value of $g_a=0$. This weak hint for weaker gravity at low $z$ is consistent with similar evidence from growth and weak lensing data.

The $E_G$ statistic is a powerful probe for detecting deviations from GR by combining weak lensing (WL), real-space clustering and redshift space distortion (RSD) measurements thus probing both the lensing and the growth effective Newton constants ($G_L$ and $G_{eff}$). We construct an up to date compilation of $E_G$ statistic data including both redshift and scale dependence ($E_G(R,z)$). We combine this $E_G$ data compilation with an up to date compilation of $f\sigma_8$ data from RSD observations to identify the current level of tension between the Planck/$\Lambda$CDM standard model based on general relativity and a general model independent redshift evolution parametrization of $G_L$ and $G_{eff}$. Each $f\sigma_8$ datapoint considered has been published separately in the context of independent analyses of distinct galaxy samples. However, there are correlations among the datapoints considered due to overlap of the analyzed galaxy samples. Due to these correlations the derived levels of tension of the best fit parameters with Planck/$\Lambda$CDM are somewhat overestimated but this is the price to pay for maximizing the information encoded in the compilation considered. We find that the level of tension increases from about $3.5\sigma$ for the $f\sigma_8$ data compilation alone to about $6\sigma$ when the $E_G$ data are also included in the analysis. The direction of the tension is the same as implied by the $f\sigma_8$ RSD growth data alone (lower $\Omega_m$ and/or weaker effective Newton constant at low redshifts for both the lensing and the growth effective Newton constants ($G_L$ and $G_{eff}$)). These results further amplify the hints for weakening modified gravity discussed in other recent analyses.

The existence of the cosmological particle horizon as the maximum measurable length $l_{max}$ in the universe leads to a generalization of the quantum uncertainty principle (GUP) to the form $\Delta x \Delta p \geq \frac{\hbar}{2}\frac{1}{1-\alpha\Delta x^2} $, where $\alpha\equiv l_{max}^{-2}$. The effects of this GUP on simple quantum mechanical systems has been shown recently by one of the authors\cite{Perivolaropoulos:2017rgq} to be extremely small (beyond current measurements) due to the extremely large scale of the current particle horizon. This is not the case in the Early Universe during the quantum generation of the inflationary primordial fluctuation spectrum. We estimate the effects of such GUP on the primordial fluctuation spectrum and on the corresponding spectral index. We generalize the field commutation (GFC) relation to $[\varphi(k),\pi_{\varphi}(k')]$=$i\delta(k-k')\frac{1}{1-\mu\varphi^2(k)}$, where $\mu\sim \alpha^2\equiv l_{max}^{-4}$ is a GFC parameter, $\varphi$ denotes a scalar field and $\pi_{\varphi}$ denotes its canonical conjugate momentum. We obtain the predicted primordial perturbation spectrum as $P_S(k)=P_S^{(0)}(k)\left(1+\frac{\bar{\mu}}{k}\right)$ where $\bar{\mu}\equiv\mu V_* \simeq \sqrt{\alpha}= l_{max}^{-1}$ (here $V_*\simeq l_{max}^3$ is the volume corresponding to $l_{max}$) and $P_S^{(0)}(k)$ is the standard primordial spectrum obtained in the context of the Heisenberg uncertainty principle ($\mu=0$). We show that the predicted scalar spectral index is $n_s=1-\lambda-\frac{\bar{\mu}}{k}$ where $\lambda$ is a slow-roll parameter. Using observational constraints on the scale dependence of the spectral index $n_s$ we show that the $2\sigma$ range of $\alpha$ corresponds to $l_{max}\gtrsim 10^{26} m $ which is of the same order as the current particle horizon.

Dynamical observational probes of the growth of density perturbations indicate that gravity may be getting weaker at low redshifts $z$. This evidence is at about $2-3\sigma$ level and comes mainly from weak lensing data that measure the parameter $S_8=\sigma_8 \sqrt{\Omega_{0m}/0.3}$ and redshift space distortion data that measure the growth rate times the amplitude of the linear power spectrum parameter $f\sigma_8 (z)$. The measured $f\sigma_8$ appears to be lower than the prediction of General Relativity (GR) in the context of the standard $\Lambda$CDM model as defined by the Planck best fit parameter values. This is the well known $f\sigma_8$ tension of $\Lambda$CDM, which constitutes one of the two main large scale challenges of the model along with the $H_0$ tension. We review the observational evidence that leads to the $f\sigma_8$ tension and discuss some theoretical implications. If this tension is not a systematic effect it may be an early hint of modified gravity with an evolving effective Newton's constant $G_{eff}$ and gravitational slip parameter $\eta$. We discuss such best fit parametrizations of $G_{eff}(z)$ and point out that they can not be reproduced by simple scalar-tensor and $f(R)$ modified gravity theories because these theories generically predict stronger gravity than General Relativity (GR) at low $z$ in the context of a $\Lambda$CDM background $H(z)$. Finally, we show weak evidence for an evolving reduced absolute magnitude of the SnIa of the Pantheon dataset at low redshifts ($z<0.1$) which may also be explained by a reduced strength of gravity and may help resolve the $H_0$ tension.

General Relativity (GR) is consistent with a wide range of experiments/observations from millimeter scales up to galactic scales and beyond. However, there are reasons to believe that GR may need to be modified because it includes singularities (it is an incomplete theory) and also it requires fine-tuning to explain the accelerating expansion of the universe through the cosmological constant. Thus, it is important to check various experiments and observations beyond the above range of scales for possible hints of deviations from the predictions of GR. If such hints are found it is important to understand which classes of modified gravity theories are consistent with them. The goal of this review is to summarize recent progress on these issues. On sub millimeter scales we show an analysis of the data of the Washington experiment (Kapner et al. (2007)) searching for modifications of Newton's Law on sub-millimeter scales and demonstrate that a spatially oscillating signal is hidden in this dataset. We show that even though this signal cannot be explained in the context of standard modified theories (viable scalar tensor and $f(R)$ theories), it is a rather generic prediction of nonlocal gravity theories. On cosmological scales we review recent analyses of Redshift Space Distortion data which measure the growth rate of cosmological perturbations at various redshifts and show that these data are in some tension with the $\Lambda$CDM parameter values indicated by Planck/2015 CMB data at about 3$\sigma$ level. This tension can be reduced by allowing for an evolution of the effective Newton constant that determines the growth rate of cosmological perturbations. We conclude that even though this tension between the data and the predictions of GR could be due to systematic/statistical uncertainties of the data, it could also constitute early hints of a new gravitational theory.

A cosmological observable measured in a range of redshifts can be used as a probe of a set of cosmological parameters. Given the cosmological observable and the cosmological parameter, there is an optimum range of redshifts where the observable can constrain the parameter in the most effective manner. For other redshift ranges the observable values may be degenerate with respect to the cosmological parameter values and thus inefficient in constraining the given parameter. These are blind redshift ranges. We determine the optimum and the blind redshift ranges of cosmological observables with respect to the cosmological parameters: matter density parameter $\Omega_m$, equation of state parameter $w$ and a modified gravity parameter $g_a$ which parametrizes the evolution of an effective Newton's constant. We consider the observables: growth rate of matter density perturbations expressed through $f(z)$ and $f\sigma_8$, the distance modulus $\mu(z)$, Baryon Acoustic Oscillation observables $D_V(z) \times \frac{r_s^{fid}}{r_s}$, $H \times \frac{r_s}{r_s^{fid}}$ and $D_A \times \frac{r_s^{fid}}{r_s}$, $H(z)$ measurements and the gravitational wave luminosity distance. We introduce a new statistic $S_P^O(z)\equiv \frac{\Delta O}{\Delta P}(z) \cdot V_{eff}^{1/2}$, including the effective survey volume $V_{eff}$, as a measure of the constraining power of a given observable $O$ with respect to a cosmological parameter $P$ as a function of redshift $z$. We find blind redshift spots $z_b$ ($S_P^O(z_b)\simeq 0$) and optimal redshift spots $z_s$ ($S_P^O(z_s)\simeq max$) for these observables with respect to the parameters $\Omega_m$, $w$ and $g_a$. For $O=f\sigma_8$ and $P=(\Omega_{m},w,g_a)$ we find blind spots at $z_b\simeq(1,2,2.7)$ respectively and optimal (sweet) spots at $z_s=(0.5,0.8,1.2)$. Thus probing higher redshifts may be less effective than probing lower redshifts with higher accuracy.

Using the Navarro-Frenk-White (NFW) dark matter density profile we reconstruct an effective field theory model for gravity at large distances from a central object by demanding that the vacuum solution has the same gravitational properties as the NFW density profile has in the context of General Relativity (GR). The dimensionally reduced reconstructed action for gravity leads to a vacuum metric that includes a modified Rindler acceleration term in addition to the Schwarzschild and cosmological constant terms. The new term is free from infrared curvature singularities and leads to a much better fit of observed galaxy velocity rotation curves than the corresponding simple Rindler term of the Grumiller metric, at the expense of one additional parameter. When the new parameter is set to zero the new metric term reduces to a Rindler constant acceleration term. We use galactic velocity rotation data to find the best fit values of the parameters of the reconstructed geometric potential and discuss possible cosmological implications.

We investigate analytically and numerically the orbits of spinning particles around black holes in the post Newtonian limit and in the presence of cosmic expansion. We show that orbits that are circular in the absence of spin, get deformed when the orbiting particle has spin. We show that the origin of this deformation is twofold: a. the background expansion rate which induces an attractive (repulsive) interaction due to the cosmic background fluid when the expansion is decelerating (accelerating) and b. a spin-orbit interaction which can be attractive or repulsive depending on the relative orientation between spin and orbital angular momentum and on the expansion rate.

A recent analysis by one of the authors\cite{Perivolaropoulos:2018cgr} has pointed out that Derrick's theorem can be evaded in curved space. Here we extend that analysis by demonstrating the existence of a static metastable solution in a wide class of metrics that include a Schwarzschild-Rindler-AntideSitter spacetime (Grumiller metric) defined as $ds^2= f(r) dt^2 - f(r)^{-1} dr^2 - r^2 (d\theta^2 +\sin^2\theta d\phi^2)$ with $f(r)=1-\frac{2Gm}{r}+2br-\frac{\Lambda}{3} r^2$ ($\Lambda<0\; b<0$). This metric emerges generically as a spherically symmetric vacuum solution in a class of scalar-tensor theories\cite{Grumiller:2010bz} as well as in Weyl conformal gravity\cite{Mannheim:1988dj}. It also emerges in General Relativity (GR) in the presence of a cosmological constant and a proper spherically symmetric perfect fluid. We demonstrate that this metric supports a static spherically symmetric metastable soliton scalar field solution that corresponds to a spherical domain wall. We derive the static solution numerically and identify a range of parameters $m, b, \Lambda$ of the metric for which the spherical wall is metastable. Our result is supported by both a minimization of the scalar field energy functional with proper boundary conditions and by a numerical simulation of the scalar field evolution. The metastable solution is very well approximated as $\phi(r) = Tanh\left[q (r-r_0)\right]$ where $r_0$ is the radius of the metastable wall that depends on the parameters of the metric and $q$ determines the width of the wall. We also find the gravitational effects of the thin spherical wall solution and its backreaction on the background metric that allows its formation. We show that this backreaction does not hinder the metastability of the solution even though it can change the range of parameters that correspond to metastability.

Recent analyses \cite{Nesseris:2017vor,Kazantzidis:2018rnb} have indicated that an effective Newton's constant $G_{\rm eff}(z)$ decreasing with redshift may relieve the observed tension between the Planck15 best fit $\Lambda$CDM cosmological background ({\it i.e.} Planck15/$\Lambda$CDM) and the corresponding $\Lambda$CDM background favored by growth $f\sigma_8$ and weak lensing data. We investigate the consistency of such a decreasing $G_{\rm eff}(z)$ with some viable scalar-tensor models and $f(R)$ theories. We stress that $f(R)$ theories generically can not lead to a decreasing $G_{\rm eff}(z)$ for any cosmological background. For scalar-tensor models we deduce that in the context of a $\Lambda$CDM cosmological background, a decreasing $G_{\rm eff}(z)$ is not consistent with a large Brans-Dicke parameter $\omega_{BD,0}$ today. This inconsistency remains and amplifies in the presence of a phantom dark energy equation of state parameter ($w < -1$). However it can be avoided for $w >-1$. We also find that any modified gravity model with the required decreasing $G_{\rm eff}(z)$ and $G_{{\rm eff},0}=G$, would have a characteristic signature in its growth index $\gamma$ with $0.61\lesssim \gamma_0\lesssim 0.69$ and large slopes $\gamma_0'$, $0.16\lesssim \gamma_0'\lesssim 0.4$, which is a characteristic signature of a decreasing (with $z$) $G_{\rm eff}(z)<G$ on small redshifts. This is a substantial departure today from the quasi-static behaviour in $\Lambda$CDM with $(\gamma_0,\gamma_0')\approx (0.55,-0.02)$.

We construct an updated extended compilation of distinct (but possibly correlated) $f\sigma_8(z)$ Redshift Space Distortion (RSD) data published between 2006 and 2018. It consists of 63 datapoints and is significantly larger than previously used similar datasets. After fiducial model correction we obtain the best fit $\Omega_{0m}-\sigma_8$ $\Lambda$CDM parameters and show that they are at a $5\sigma$ tension with the corresponding Planck15/$\Lambda$CDM values. Introducing a nontrivial covariance matrix correlating randomly $20\%$ of the RSD datapoints has no significant effect on the above tension level. We show that the tension disappears (becomes less than $1\sigma$) when a subsample of the 20 most recently published data is used. A partial cause for this reduced tension is the fact that more recent data tend to probe higher redshifts (with higher errorbars) where there is degeneracy among different models due to matter domination. Allowing for a nontrivial evolution of the effective Newton's constant as $G_{\textrm{eff}}(z)/G_{\textrm{N}} = 1 + g_a \left(\frac{z}{1+z}\right)^2 - g_a \left(\frac{z}{1+z}\right)^4$ ($g_a$ is a parameter) and fixing a \plcdm background we find $g_a=-0.91\pm 0.17$ from the full $f\sigma_8$ dataset while the 20 earliest and 20 latest datapoints imply $g_a=-1.28^{+0.28}_{-0.26}$ and $g_a=-0.43^{+0.46}_{-0.41}$ respectively. Thus, the more recent $f\sigma_8$ data appear to favor GR in contrast to earlier data. Finally, we show that the parametrization $f\sigma_8(z)=\lambda \sigma_8 \Omega(z)^\gamma /(1+z)^\beta$ provides an excellent fit to the solution of the growth equation for both GR ($g_a=0$) and modified gravity ($g_a\neq 0$).

It is well known that global topological defects induce a repulsive gravitational potential for test particles. 'What is the gravitational potential induced by black holes with a cosmological constant (Schwarzschild-de Sitter (S-dS) metric) on finite thickness global topological defects?'. This is the main question addressed in the present analysis. We also discuss the validity of Derrick's theorem when scalar fields are embedded in non-trivial gravitational backgrounds. In the context of the above question, we consider three global defect configurations: a finite thickness spherical domain wall with a central S-dS black hole, a global string loop with a S-dS black hole in the center and a global monopole near a S-dS black hole. Using an analytical model and numerical simulations of the evolving spherical wall we show that the spherical wall experiences a repelling gravitational potential due to the mass of the central black hole. This potential is further amplified by the presence of a cosmological constant. For initial domain wall radius larger than a critical value, the repulsive potential dominates over the wall tension and the wall expands towards the cosmological horizon of the S-dS metric where it develops ghost instabilities. For smaller initial radius, tension dominates and the wall contracts towards the black hole horizon where it also develops ghost instabilities. We also show, using the same analytical model and energetic arguments that a global monopole is gravitationally attracted by a black hole while a cosmological constant induces a repulsive gravitational potential as in the case of test particles. Finally we show that a global string loop with finite thickness experiences gravitational repulsion due to the cosmological constant which dominates over its tension for a radius larger than a critical radius leading to an expanding rather than contracting loop.

A recent analysis by one of the authors\cite{Perivolaropoulos:2016ucs} has indicated the presence of a $2\sigma$ signal of spatially oscillating new force residuals in the torsion balance data of the Washington experiment. We extend that study and analyse the data of the Stanford Optically Levitated Microsphere Experiment (SOLME) \cite{Rider:2016xaq} (kindly provided by the authors of \cite{Rider:2016xaq}) searching for sub-mm spatially oscillating new force signals. We find a statistically significant oscillating signal for a force residual of the form $F(z)=\alpha \; cos(\frac{2\pi}{\lambda}\; z +c)$ where $z$ is the distance between the macroscopic interacting masses (levitated microsphere and cantilever). The best fit parameter values are $\alpha=(1.1 \pm 0.4)\times 10^{-17}N$, $\lambda=(35.2\pm 0.6)\mu m$. Monte Carlo simulation of the SOLME data under the assumption of zero force residuals has indicated that the statistical significance of this signal is at about $2\sigma$ level. Private communication with members of the SOLME group has indicated that this previously unnoticed signal is most probably due to a systematic effect (diffraction of non-Gaussian tails of the laser beam). Thus it can only be useful as an upper bound of the amplitude of new spatially oscillating forces on sub-mm scales. Assuming gravitational origin emerging from a fundamental modification of the Newtonian potential of the form $V_{eff}(r)=-G\frac{M}{r}(1+\alpha_O \cos(\frac{2\pi}{\lambda}\; r+\theta))\equiv V_{N}(r)+V_{osc}(r)$, we evaluate the source integral and obtain a bound $\alpha_O < 10^7$ for $\lambda \simeq 35 \mu m$. Thus, we get no useful constraints on the modified gravity parameter $\alpha_O$. However, the constraints on the general phenomenological parameter $\alpha$ can be useful in constraining other fifth force models related to dark energy (chameleon oscillating potentials etc).

We construct an updated and extended compilation of growth rate data consisting of 34 points and including corrections for model dependence. In order to maximize the independence of the datapoints we also construct a subsample of this compilation (`Gold' growth dataset) which consists of 18 datapoints. We test the consistency of this dataset with the best fit Planck15/$\Lambda$CDM parameters in the context of General Relativity (GR) using the evolution equation for the growth factor $\delta(a)$ with a $w$CDM background. We find tension at the $\sim 3 \sigma$ level between the best fit parameters $w$ (the dark energy equation of state), $\Omega_{0m}$ (the matter density parameter) and $\sigma_8$ (the matter power spectrum normalization on scales $8h^{-1}$Mpc) and the corresponding Planck15/$\Lambda$CDM parameters. We show that the tension disappears if we allow for evolution of the effective Newton's constant, parametrized as $G_{eff}(a)/G_N = 1 + g_a(1-a)^n-g_a(1-a)^{2n}$ with $n\ge2$ where $g_a$, $n$ are parameters of the model, $a$ is the scale factor and $z = 1/a-1$ is the redshift. This parametrization satisfies three criteria: a. $G_{eff} > 0$, b. Consistency with Big Bang Nucleosynthesis ($G_{eff}(a\ll 1)/G_N=1$), c. Consistency with solar system tests ($G_{eff}(a=1)/G_N=1$ and $G_{eff}'(a=1)/G_N=0$). We show that the best fit form of $G_{eff}(z)$ obtained from the growth data corresponds to weakening gravity at recent redshifts (decreasing function of $z$) and we demonstrate that this behavior is not consistent with any scalar-tensor Lagrangian with a real scalar field. Finally, we use MGCAMB to find the best fit $G_{eff}(z)$ obtained from the Planck CMB power spectrum on large angular scales and show that it is a mildly increasing function of $z$, in $3\sigma$ tension with the corresponding decreasing best fit $G_{eff}(z)$ obtained from the growth data.

We investigate the viability of sub-millimeter wavelength oscillating deviations from the Newtonian potential at both the theoretical and the experimental/observational level. At the theoretical level such deviations are generic predictions in a wide range of extensions of General Relativity (GR) including $f(R)$ theories, massive Brans-Dicke theories, compactified extra dimension models and nonlocal extensions of GR. However, the range of parameters associated with such oscillating deviations is usually connected with instabilities. An exception emerges in nonlocal gravity theories where oscillating deviations from Newtonian potential occur naturally on sub-millimeter scales without instabilities. As an example of a model with unstable Newtonian oscillations we review an $f(R)$ expansion around General Relativity of the form $f(R)=R+\frac{1}{6 m^2} R^2$ with $m^2<0$ pointing out possible stabilization mechanisms. As an example of a model with stable Newtonian oscillations we discuss nonlocal gravity theories. If such oscillations are realized in Nature on sub-millimeter scales, a signature is expected in torsion balance experiments testing Newton's law. We search for such a signature in the torsion balance data of the Washington experiment (combined torque residuals) testing Newton's law at sub-millimeter scales. We show that an oscillating residual ansatz with spatial wavelength $\lambda \simeq 0.1mm$ provides a better fit to the data compared to the residual Newtonian constant ansatz by $\Delta \chi^2 = -15$. Similar improved fits, also occur in about $10\%$ of Monte Carlo realizations of Newtonian data. Thus, the significance level of this improved fit is at a level of not more than $2\sigma$. The energy scale corresponding to this best fit wavelength is identical to the dark energy length scale $\lambda_{de} \equiv\sqrt[4]{\hbar c/\rho_{ de}}\approx 0.1mm$.

The cosmological particle horizon is the maximum measurable length in the Universe. The existence of such a maximum observable length scale implies a modification of the quantum uncertainty principle. Thus due to non-locality of quantum mechanics, the global properties of the Universe could produce a signature on the behaviour of local quantum systems. A Generalized Uncertainty Principle (GUP) that is consistent with the existence of such a maximum observable length scale $l_{max}$ is $\Delta x \Delta p \geq \frac{\hbar}{2}\;\frac{1}{1-\alpha \Delta x^2}$ where $\alpha = l_{max}^{-2}\simeq (H_0/c)^2$ ($H_0$ is the Hubble parameter and $c$ is the speed of light). In addition to the existence of a maximum measurable length $l_{max}=\frac{1}{\sqrt \alpha}$, this form of GUP implies also the existence of a minimum measurable momentum $p_{min}=\frac{3 \sqrt{3}}{4}\hbar \sqrt{\alpha}$. Using appropriate representation of the position and momentum quantum operators we show that the spectrum of the one dimensional harmonic oscillator becomes $\bar{\mathcal{E}}_n=2n+1+\lambda_n \bar{\alpha}$ where $\bar{\mathcal{E}}_n\equiv 2E_n/\hbar \omega$ is the dimensionless properly normalized $n^{th}$ energy level, $\bar{\alpha}$ is a dimensionless parameter with $\bar{\alpha}\equiv \alpha \hbar/m \omega$ and $\lambda_n\sim n^2$ for $n\gg 1$ (we show the full form of $\lambda_n$ in the text). For a typical vibrating diatomic molecule and $l_{max}=c/H_0$ we find $\bar{\alpha}\sim 10^{-77}$ and therefore for such a system, this effect is beyond reach of current experiments. However, this effect could be more important in the early universe and could produce signatures in the primordial perturbation spectrum induced by quantum fluctuations of the inflaton field.

Sudden singularities occur in FRW spacetimes when the scale factor remains finite and different from zero while some of its derivatives diverge. After proper rescaling, the scale factor close to such a singularity at $t=0$ takes the form $a(t)=1+ c |t|^\eta$ (where $c$ and $\eta$ are parameters and $\eta\geq 0$). We investigate analytically and numerically the geodesics of free and gravitationally bound particles through such sudden singularities. We find that even though free particle geodesics go through sudden singularities for all $\eta\geq 0$, bound systems get dissociated for a wide range of the parameter $c$. For $\eta <1$ bound particles receive a diverging impulse at the singularity and get dissociated for all positive values of the parameter $c$. For $\eta > 1$ (Sudden Future Singularities (SFS)) bound systems get a finite impulse that depends on the value of $c$ and get dissociated for values of $c$ larger than a critical value $c_{cr}(\eta,\omega_0)>0$ that increases with the value of $\eta$ and the rescaled angular velocity $\omega_0$ of the bound system. We obtain an approximate equation for the analytical estimate of $c_{cr}(\eta,\omega_0)$. We also obtain its accurate form by numerical derivation of the bound system orbits through the singularities. Bound system orbits through Big Brake singularities ($c<0$, $1<\eta<2$) are also derived numerically and are found to get disrupted at the singularity. However, they remain bound for all values of the parameter $c$ considered.

Dark energy equation of state $w(z)$ parametrizations with two parameters and given monotonicity are generically either convex or concave functions. This makes them suitable for fitting either freezing or thawing quintessence models but not both simultaneously. Fitting a dataset based on a freezing model with an unsuitable (concave when increasing) $w(z)$ parametrization (like CPL) can lead to significant misleading features like crossing of the phantom divide line, incorrect $w(z=0)$, incorrect slope \etc that are not present in the underlying cosmological model. To demonstrate this fact we generate scattered cosmological data both at the level of $w(z)$ and the luminosity distance $D_L(z)$ based on either thawing or freezing quintessence models and fit them using parametrizations of convex and of concave type. We then compare statistically significant features of the best fit $w(z)$ with actual features of the underlying model. We thus verify that the use of unsuitable parametrizations can lead to misleading conclusions. In order to avoid these problems it is important to either use both convex and concave parametrizations and select the one with the best $\chi^2$ or use principal component analysis thus splitting the redshift range into independent bins. In the latter case however, significant information about the slope of $w(z)$ at high redshifts is lost. Finally, we propose a new family of parametrizations (nCPL) $w(z)=w_0+w_a (\frac{z}{1+z})^n$ which generalizes the CPL and interpolates between thawing and freezing parametrizations as the parameter $n$ increases to values larger than 1.

We investigate the geodesics of a Schwarzschild spacetime embedded in an isotropic expanding cosmological background (McVittie metric). We focus on bound particle geodesics in a background including matter and phantom dark energy with constant dark energy equation of state parameter $w<-1$ involving a future Big Rip singularity at a time $t_{\ast}$. Such geodesics have been previously studied in the Newtonian approximation and found to lead to dissociation of bound systems at a time $t_{rip}<t_{\ast}$ which for fixed background $w$, depends on a single dimensionless parameter $\bar{\omega}_0$ related to the angular momentum and depending on the mass and the size of the bound system. We extend this analysis to large massive bound systems where the Newtonian approximation is not appropriate and we compare the derived dissociation time with the corresponding time in the context of the Newtonian approximation. By identifying the time when the effective potential minimum disappears due to the repulsive force of dark energy we find that the dissociation time of bound systems occurs earlier than the prediction of the Newtonian approximation. However, the effect is negligible for all existing cosmological bound systems and it would become important only in hypothetical bound extremely massive ($10^{20} M_\odot$) and large ($100Mpc$) bound systems. We verify this result by explicit solution of the geodesic equations. This result is due to an interplay between the repulsive phantom dark energy effects and the existence of the well known innermost stable orbits of Schwarzschild spacetimes.

All quintessence potentials that are either monotonic with negative interval or have a minimum at negative values of the potential, generically predict a future collapse of the scale factor to a "doomsday" singularity. We show that this doomsday is generically avoided in models with a proper non-minimal coupling of the quintessence scalar field to the curvature scalar $R$. For simplicity we consider linear quintessence potential $V=-s\phi$ and linear non-minimal coupling $F=1-\lambda \phi$. However our result is generic and is due to the fact that the non-minimal coupling modifies the effective potential that determines the dynamics of the scalar field. Thus for each positive value of the parameter $s$ we find a critical value $\lambda_{crit}(s)$ such that for $\lambda>\lambda_{crit}(s)$ the negative potential energy does not dominate the universe and the cosmic doomsday Big Crunch singularity is avoided because the scalar field eventually rolls up its potential. We find that $\lambda_{crit}(s)$ increases approximately linearly with $s$. For $\lambda>\lambda_{crit}(s)$ the potential energy of the scalar field becomes positive and it eventually dominates while the dark energy equation of state parameter tends to $w=-1$ leading to a deSitter Universe.

We consider Lagrangians in 3+1 dimensions admitting topological defects where there is an additional coupling between the defect scalar field $\Phi$ and the gauge field kinetic term (eg $B(\vert \Phi \vert^2) F_{\mu \nu}F^{\mu \nu}$). Such a {\it dilatonic} coupling in the context of a static defect, induces a spatially dependent effective gauge charge and effective mass for the scalar field which leads to modified properties of the defect core. In particular, the scale of the core gets modified while the stability properties of the corresponding embedded defects are also affected. These modifications are illustrated for gauged (Nielsen-Olesen) vortices and for gauged ('t Hooft-Polyakov) monopoles. The corresponding dilatonic global defects are also studied in the presence of an external gauge field.

A wide range of large scale observations hint towards possible modifications on the standard cosmological model which is based on a homogeneous and isotropic universe with a small cosmological constant and matter. These observations, also known as "cosmic anomalies" include unexpected Cosmic Microwave Background perturbations on large angular scales, large dipolar peculiar velocity flows of galaxies ("bulk flows"), the measurement of inhomogenous values of the fine structure constant on cosmological scales ("alpha dipole") and other effects. The presence of the observational anomalies could either be a large statistical fluctuation in the context of {\lcdm} or it could indicate a non-trivial departure from the cosmological principle on Hubble scales. Such a departure is very much constrained by cosmological observations for matter. For dark energy however there are no significant observational constraints for Hubble scale inhomogeneities. In this brief review I discuss some of the theoretical models that can naturally lead to inhomogeneous dark energy, their observational constraints and their potential to explain the large scale cosmic anomalies.

We investigate the cosmological predictions of several $f(T)$ models, with up to two parameters, at both the background and the perturbation levels. Using current cosmological observations (geometric supernovae type Ia, cosmic microwave background and baryonic acoustic oscillation and dynamical growth data) we impose constraints on the distortion parameter, which quantifies the deviation of these models from the concordance $\Lambda$ cosmology at the background level. In addition we constrain the growth index $\gamma$ predicted in the context of these models using the latest perturbation growth data in the context of three parametrizations for $\gamma$. The evolution of the best fit effective Newton constant, which incorporates the $f(T)$-gravity effects, is also obtained along with the corresponding $1\sigma$ error regions. We show that all the viable parameter sectors of the $f(T)$ gravity models considered practically reduce these models to $\Lambda$CDM. Thus, the degrees of freedom that open up to $\Lambda$CDM in the context of $f(T)$ gravity models are not utilized by the cosmological data leading to an overall disfavor of these models.

We demonstrate that the stability of the semilocal vortex can be significantly improved by the presence of a dilatonic coupling of the form e^\frac{q | \Phi |^2}{\eta^2} F_{\mu \nu}F^{\mu \nu} with q>0 where \eta is the scale of symmetry breaking that gives rise to the vortex. For q=0 we obtain the usual embedded (semilocal) Nielsen-Olesen vortex. We find the stability region of the parameter \beta = (\frac{m_\Phi}{m_A})^2 (m_\Phi and m_A are the masses of the scalar and gauge fields respectively). We show that the stability region of \beta is 0<\beta<\beta_{max}(q) where \beta_{max}(q=0)=1 (as expected) and \beta_{max}(q) is an increasing function of q. This result may have significant implications for the stability of the electroweak vortex in the presence of a dilatonic coupling (dilatonic electroweak vortex).

We demonstrate that a wide range of viable f(R) parameterizations (including the Hu & Sawicki and the Starobinsky models) can be expressed as perturbations deviating from the LCDM Lagrangian. We constrain the deviation parameter b using a combination of geometrical and dynamical observational probes. In particular, we perform a joint likelihood analysis of the recent Supernovae Type Ia data, the Cosmic Microwave Background shift parameters, the Baryonic Acoustic Oscillations and the growth rate data provided by the various galaxy surveys. This analysis provides constraints for the following parameters: the matter density Omega_{m0}, the deviation from LCDM parameter b and the growth index gamma(z). We parametrize the growth index gamma(z) in three manners (constant, Taylor expansion around z=0, and Taylor expansion around the scale factor). We point out the numerical difficulty for solving the generalized f(R) Friedman equation at high redshifts due to stiffness of the resulting ordinary differential equation. We resolve this problem by constructing an efficient analytical perturbative method in the deviation parameter b. We demonstrate that this method is highly accurate, by comparing the resulting analytical expressions for the Hubble parameter, with the numerical solutions at low and intermediate redshifts. Surprisingly, despite of its perturbative nature, the accuracy of the method persists even for values of b that are of O(1).

We use a global pixel based estimator to identify the axis of the residual Maximum Temperature Asymmetry (MTA) (after the dipole subtraction) of the WMAP 7 year Internal Linear Combination (ILC) CMB temperature sky map. The estimator is based on considering the temperature differences between opposite pixels in the sky at various angular resolutions (4 degrees-15 degrees and selecting the axis that maximizes this difference. We consider three large scale Healpix resolutions (N_{side}=16 (3.7 degrees), N_{side}=8 (7.3 degrees) and N_{side}=4 (14.7 degrees)). We compare the direction and magnitude of this asymmetry with three other cosmic asymmetry axes (\alpha dipole, Dark Energy Dipole and Dark Flow) and find that the four asymmetry axes are abnormally close to each other. We compare the observed MTA axis with the corresponding MTA axes of 10^4 Gaussian isotropic simulated ILC maps (based on LCDM). The fraction of simulated ILC maps that reproduces the observed magnitude of the MTA asymmetry and alignment with the observed \alpha dipole is in the range of 0.1%-0.5%$ (depending on the resolution chosen for the CMB map). The corresponding magnitude+alignment probabilities with the other two asymmetry axes (Dark Energy Dipole and Dark Flow) are at the level of about 1%. We propose Extended Topological Quintessence as a physical model qualitatively consistent with this coincidence of directions.

We present a detailed analysis (including redshift tomography) of the cosmic dipoles in the Keck+VLT quasar absorber and in the Union2 SnIa samples. We show that the fine structure constant cosmic dipole obtained through the Keck+VLT quasar absorber sample at $4.1\sigma$ level is anomalously aligned with the corresponding dark energy dipole obtained through the Union2 sample at $2\sigma$ level. The angular separation between the two dipole directions is $11.3^\circ \pm 11.8^\circ$. We use Monte Carlo simulations to find the probability of obtaining the observed dipole magnitudes with the observed alignment, in the context of an isotropic cosmological model with no correlation between dark energy and fine structure constant $\alpha$. We find that this probability is less than one part in $10^6$. We propose a simple physical model (extended topological quintessence) which naturally predicts a spherical inhomogeneous distribution for both dark energy density and fine structure constant values. The model is based on the existence of a recently formed giant global monopole with Hubble scale core which also couples non-minimally to electromagnetism. Aligned dipole anisotropies would naturally emerge for an off-centre observer for both the fine structure constant and for dark energy density. This model smoothly reduces to \lcdm for proper limits of its parameters. Two predictions of this model are (a) a correlation between the existence of strong cosmic electromagnetic fields and the value of $\alpha$ and (b) the existence of a dark flow on Hubble scales due to the repulsive gravity of the global defect core (`Great Repulser') aligned with the dark energy and $\alpha$ dipoles. The direction of the dark flow is predicted to be towards the spatial region of lower accelerating expansion. Existing data about the dark flow are consistent with this prediction.

We investigate the cosmic evolution of the linear bias in the framework of a flat FLRW spacetime. We consider metric perturbations in the Newtonian gauge, including Hubble scale effects. Making the following assumptions, (i) scale independent current epoch bias $b_0$, (ii) equal accelerations between tracers and matter, (iii) unimportant halo merging effects (which is quite accurate for $z<3$), we analytically derive the scale dependent bias evolution. The identified scale dependence is only due to Hubble scale evolution GR effects, while other scale dependence contributions are ignored. We find that up to galaxy cluster scales the fluctuations of the metric do not introduce a significant scale dependence in the linear bias. Our bias evolution model is then used to derive a connection between the matter growth index $\gamma$ and the observable value of the tracer power spectrum normalization $\sigma_8(z)$. We show how this connection can be used as an observational test of General Relativity on extragalactic scales.

A global monopole (or other topological defect) formed during a recent phase transition with core size comparable to the present Hubble scale, could induce the observed accelerating expansion of the universe. In such a model, topological considerations trap the scalar field close to a local maximum of its potential in a cosmologically large region of space. We perform detailed numerical simulations of such an inhomogeneous dark energy system (topological quintessence) minimally coupled to gravity, in a flat background of initially homogeneous matter. We find that when the energy density of the field in the monopole core starts dominating the background density, the spacetime in the core starts to accelerate its expansion in accordance to a \Lambda CDM model with an effective inhomogeneous spherical dark energy density parameter \Omega_\Lambda(r). The matter density profile is found to respond to the global monopole profile via an anti-correlation (matter underdensity in the monopole core). Away from the monopole core, the spacetime is effectively Einstein-deSitter (\Omega_\Lambda(r_{out}) -> 0) while at the center \Omega_\Lambda(r ~ 0) is maximum. We fit the numerically obtained expansion rate at the monopole core to the Union2 data and show that the quality of fit is almost identical to that of \Lambda CDM. Finally, we discuss potential observational signatures of this class of inhomogeneous dark energy models.

The consistency level of LCDM with geometrical data probes has been increasing with time during the last decade. Despite of these successes, there are some puzzling conflicts between LCDM predictions and dynamical data probes (bulk flows, alignment and magnitude of low CMB multipoles, alignment of quasar optical polarization vectors, cluster halo profiles). Most of these puzzles are related to the existence of preferred anisotropy axes which appear to be unlikely close to each other. A few models that predict the existence of preferred cosmological axes are briefly discussed.

We consider on-center and off-center observers in an inhomogeneous, spherically symmetric, isocurvature (flat) concentration of dark energy with typical size of a few Gpc. Such a concentration could be produced e.g. by a recently formed global monopole with core size that approaches the Hubble scale. In this case we would have what may be called `topological quintessence' in analogy with the well-known topological inflation. We show that the minimum comoving radius r_{0min} of such a dark energy inhomogeneity that is consistent with the Union2 Type Ia supernovae (SnIa) data at the 3\sigma level is r_{0min}\simeq 1.8 Gpc. As expected, the best-fit fractional dark energy density at the center, \Omega_X,in, approaches the corresponding LCDM value \Omega_X,in =0.73 for large enough values of the inhomogeneity radius r_0 (r_0 > 4Gpc). Using the Union2 data, we show that the maximum allowed shift r_{obs-max} of the observer from the center of the inhomogeneity is about 0.7 r_0 which respects the Copernican principle. The model naturally predicts the existence of a preferred axis and alignment of the low CMB multipoles. However, the constraints on r_{obs-max} coming from the magnitude of the CMB dipole remain a severe challenge to the Copernican principle and lead to r_{obs-max}< 110 Mpc even for an inhomogeneity radius as large as r_0=7 Gpc.

We review, compare and extend recent studies searching for evidence for a preferred cosmological axis. We start from the Union2 SnIa dataset and use the hemisphere comparison method to search for a preferred axis in the data. We find that the hemisphere of maximum accelerating expansion rate is in the direction $(l,b)=({309^\circ}^{+23^\circ}_{-3^\circ}, {18^\circ}^{+11^\circ}_{-10^\circ})$ ($\omm=0.19$) while the hemisphere of minimum acceleration is in the opposite direction $(l,b)=({129^\circ}^{+23^\circ}_{-3^\circ},{-18^\circ}^{+10^\circ}_{-11^\circ})$ ($\omm=0.30$). The level of anisotropy is described by the normalized difference of the best fit values of $\omm$ between the two hemispheres in the context of \lcdm fits. We find a maximum anisotropy level in the Union2 data of $\frac{\Delta \ommax}{\bomm}=0.43\pm 0.06$. Such a level does not necessarily correspond to statistically significant anisotropy because it is reproduced by about $30%$ of simulated isotropic data mimicking the best fit Union2 dataset. However, when combined with the axes directions of other cosmological observations (bulk velocity flow axis, three axes of CMB low multipole moments and quasar optical polarization alignment axis), the statistical evidence for a cosmological anisotropy increases dramatically. We estimate the probability that the above independent six axes directions would be so close in the sky to be less than $1%$. Thus either the relative coincidence of these six axes is a very large statistical fluctuation or there is an underlying physical or systematic reason that leads to their correlation.

We propose a scale dependent analytic approximation to the exact linear growth of density perturbations in Scalar-Tensor (ST) cosmologies. In particular, we show that on large subhorizon scales, in the Newtonian gauge, the usual scale independent subhorizon growth equation does not describe the growth of perturbations accurately, as a result of scale-dependent relativistic corrections to the Poisson equation. A comparison with exact linear numerical analysis indicates that our approximation is a significant improvement over the standard subhorizon scale independent result on large subhorizon scales. A comparison with the corresponding results in the Synchronous gauge demonstrates the validity and consistency of our analysis.

The LCDM cosmological model assumes the existence of a small cosmological constant in order to explain the observed accelerating cosmic expansion. Despite the dramatic improvement of the quality of cosmological data during the last decade it remains the simplest model that fits remarkably well (almost) all cosmological observations. In this talk I review the increasingly successful fits provided by LCDM on recent geometric probe data of the cosmic expansion. I also briefly discuss some emerging shortcomings of the model in attempting to fit specific classes of data (eg cosmic velocity dipole flows and cluster halo profiles). Finally, I summarize recent results on the theoretically predicted matter overdensity ($\delta_m=\frac{\delta \rho_m}{\rho_m}$) evolution (a dynamical probe of the cosmic expansion), emphasizing its scale and gauge dependence on large cosmological scales in the context of general relativity. A new scale dependent parametrization which describes accurately the growth rate of perturbations even on scales larger than 100h^{-1}Mpc is shown to be a straightforward generalization of the well known scale independent parametrization f(a)=\omms(a)^\gamma valid on smaller cosmological scales.

We solve analytically and numerically the generalized Einstein equations in scalar-tensor cosmologies to obtain the evolution of dark energy and matter linear perturbations. We compare our results with the corresponding results for minimally coupled quintessence perturbations. Our results for natural (O(1)) values of parameters in the Lagrangian which lead to a background expansion similar to LCDM are summarized as follows: 1. Scalar-Tensor dark energy density perturbations are amplified by a factor of about 10^6 compared to minimally coupled quintessence perturbations on scales less than about 100 h^{-1} Mpc. On these scales dark energy perturbations constitute a fraction of about 10% compared to matter density perturbations. 2. Scalar-Tensor dark energy density perturbations are anti-correlated with matter linear perturbations on sub-Hubble scales. Thus clusters of galaxies are predicted to overlap with voids of dark energy. 3. This anti-correlation of matter with negative pressure perturbations induces a mild amplification of matter perturbations by about 10% on sub-Hubble scales. 4. The evolution of scalar field perturbations on sub-Hubble scales, is scale independent and therefore it corresponds to a vanishing effective speed of sound (c_{s\Phi}=0). It also involves large oscillations at early times induced by the amplified effective mass of the field. This amplification is due to the non-minimal coupling of the field to the Ricci curvature scalar and (thus) to matter. No such oscillations are present in minimally coupled quintessence perturbations which are suppressed on sub-Hubble scales (c_{s\Phi}=1). Implications of our results include modified predictions for cluster halo profiles observed properties. These profiles are known to be in some tension with the predictions of LCDM.

Previous solar system constraints of the Brans-Dicke (BD) parameter $\omega$ have either ignored the effects of the scalar field potential (mass terms) or assumed a highly massive scalar field. Here, we interpolate between the above two assumptions and derive the solar system constraints on the BD parameter $\omega$ for {\it any} field mass. We show that for $\omega=O(1)$ the solar system constraints relax for a field mass $m \gsim 20 \times m_{AU}= 20\times 10^{-27}GeV$.

We rank the six latest Type Ia supernova (SnIa) datasets (Constitution (C), Union (U), ESSENCE (Davis) (E), Gold06 (G), SNLS 1yr (S) and SDSS-II (D)) in the context of the Chevalier-Polarski-Linder (CPL) parametrization $w(a)=w_0+w_1 (1-a)$, according to their Figure of Merit (FoM), their consistency with the cosmological constant ($\Lambda$CDM), their consistency with standard rulers (Cosmic Microwave Background (CMB) and Baryon Acoustic Oscillations (BAO)) and their mutual consistency. We find a significant improvement of the FoM (defined as the inverse area of the 95.4% parameter contour) with the number of SnIa of these datasets ((C) highest FoM, (U), (G), (D), (E), (S) lowest FoM). Standard rulers (CMB+BAO) have a better FoM by about a factor of 3, compared to the highest FoM SnIa dataset (C). We also find that the ranking sequence based on consistency with $\Lambda$CDM is identical with the corresponding ranking based on consistency with standard rulers ((S) most consistent, (D), (C), (E), (U), (G) least consistent). The ranking sequence of the datasets however changes when we consider the consistency with an expansion history corresponding to evolving dark energy $(w_0,w_1)=(-1.4,2)$ crossing the phantom divide line $w=-1$ (it is practically reversed to (G), (U), (E), (S), (D), (C)). The SALT2 and MLCS2k2 fitters are also compared and some peculiar features of the SDSS-II dataset when standardized with the MLCS2k2 fitter are pointed out. Finally, we construct a statistic to estimate the internal consistency of a collection of SnIa datasets. We find that even though there is good consistency among most samples taken from the above datasets, this consistency decreases significantly when the Gold06 (G) dataset is included in the sample.

We generalize the small scale dynamics of the universe by taking into account models with an equation of state which evolves with time, and provide a complete formulation of the cluster virialization attempting to address the nonlinear regime of structure formation. In the context of the current dark energy models, we find that galaxy clusters appear to form at z\sim 1-2, in agreement with previous studies. Also, we investigate thoroughly the evolution of spherical matter perturbations, as the latter decouple from the background expansion and start to ``turn around'' and finally collapse. Within this framework, we find that the concentration parameter depends on the choice of the considered dark energy (homogeneous or clustered). In particular, if the distribution of the dark energy is clustered then we produce more concentrated structures with respect to the homogeneous dark energy. Finally, comparing the predicted concentration parameter with the observed concentration parameter, measured for four massive galaxy clusters, we find that the scenario which contains a pure homogeneous dark energy is unable to reproduce the data. The situation becomes somewhat better in the case of an inhomogeneous (clustered) dark energy.

We demonstrate the scale dependence of the growth function of cosmological perturbations in dark energy models based on General Relativity. This scale dependence is more prominent on cosmological scales of $100h^{-1}Mpc$ or larger. We derive a new scale dependent parametrization which generalizes the well known Newtonian approximation result $f_0(a)\equiv \frac{d\ln \delta_0}{d\ln a}=\om(a)^\gamma$ ($\gamma ={6/11}$ for \lcdm) which is a good approximation on scales less than $50h^{-1}Mpc$. Our generalized parametrization is of the form $f(a)=\frac{f_0(a)}{1+\xi(a,k)}$ where $\xi(a,k)=\frac{3 H_{0}^{2} \omm}{a k^2}$. We demonstrate that this parametrization fits the exact result of a full general relativistic evaluation of the growth function up to horizon scales for both \lcdm and dynamical dark energy. In contrast, the scale independent parametrization does not provide a good fit on scales beyond 5% of the horizon scale ($k\simeq 0.01 h^{-1}Mpc$).

It has recently been pointed out by Kowalski et. al. (arxiv:0804.4142) that there is `an unexpected brightness of the SnIa data at z>1'. We quantify this statement by constructing a new statistic which is applicable directly on the Type Ia Supernova (SnIa) distance moduli. This statistic is designed to pick up systematic brightness trends of SnIa datapoints with respect to a best fit cosmological model at high redshifts. It is based on binning the normalized differences between the SnIa distance moduli and the corresponding best fit values in the context of a specific cosmological model (eg LCDM). We then focus on the highest redshift bin and extend its size towards lower redshifts until the Binned Normalized Difference (BND) changes sign (crosses 0) at a redshift z_c (bin size N_c). The bin size N_c of this crossing (the statistical variable) is then compared with the corresponding crossing bin size N_{mc} for Monte Carlo data realizations based on the best fit model. We find that the crossing bin size N_c obtained from the Union08 and Gold06 data with respect to the best fit LCDM model is anomalously large compared to N_{mc} of the corresponding Monte Carlo datasets obtained from the best fit LCDM in each case. In particular, only 2.2% of the Monte Carlo LCDM datasets are consistent with the Gold06 value of N_c while the corresponding probability for the Union08 value of N_c is 5.3%. Thus, according to this statistic, the probability that the high redshift brightness bias of the Union08 and Gold06 datasets is realized in the context of a (w_0,w_1)=(-1,0) model (LCDM cosmology) is less than 6%. The corresponding realization probability in the context of a (w_0,w_1)=(-1.4,2) model is more than 30% for both the Union08 and the Gold06 datasets.

We demonstrate numerically that an oscillation mode in 1+1 dimensions (eg a breather or an oscillon) can decay into a kink-antikink pair by a sudden distortion of the evolution potential which occurs within a certain time or space domain. In particular, we consider the transition of a sine-Gordon potential into a \Phi^4 potential. The breather field configuration is assumed to initially evolve in a sine-Gordon potential with velocity $v$ and oscillation frequency $\omega$. We then consider two types of numerical experiments: a. An abrupt transition of the potential to a $\Phi^4$ form at t_0=0 over the whole 1-dimensional lattice and b. The impact of the breather on a region x>x_0=0 where the potential has the \Phi^4 form which is different from the sine-Gordon form valid at x<x_0=0. We find that in both cases there is a region of parameters (v,\omega) such that the breather decays to a kink-antikink pair. This region of parameters for kink-antikink formation is qualitatively similar with the parameter region where the energy of the breather exceeds the energy of the kink-antikink pair in the \Phi^4 potential. We demonstrate that the same mechanism for soliton formation is realized when using a gaussian oscillator (oscillon) instead of a breather. We briefly discuss the implications of our results for realistic experiments as well as their extension to soliton formation in two and three space dimensions.

The LCDM cosmological model is a well defined, simple and predictive model which is consistent with the majority of current cosmological observations. Despite of these successes there are specific cosmological observations which differ from the predictions of LCDM at a level of 2\sigma or higher. These observations include the following: 1. Large Scale Velocity Flows (LCDM predicts significantly smaller amplitude and scale of flows than what observations indicate), 2. Brightness of Type Ia Supernovae (SnIa) at High Redshift z (LCDM predicts fainter SnIa at High z), 3. Emptiness of Voids (LCDM predicts more dwarf or irregular galaxies in voids than observed), 4. Profiles of Cluster Haloes (LCDM predicts shallow low concentration and density profiles in contrast to observations which indicate denser high concentration cluster haloes) 5. Profiles of Galaxy Haloes (LCDM predicts halo mass profiles with cuspy cores and low outer density while lensing and dynamical observations indicate a central core of constant density and a flattish high dark mass density outer profile), 6. Sizable Population of Disk Galaxies (LCDM predicts a smaller fraction of disk galaxies due to recent mergers expected to disrupt cold rotationally supported disks). Even though the origin of some of the above challenges may be astrophysical or related to dark matter properties, it should be stressed that even on galactic and cluster scales, the effects of dark energy on the equilibrium and stability of astrophysical systems are not negligible and they may play a key role in the resolution of the above puzzles. Here, I briefly review these six challenges of LCDM and discuss the possible dark energy properties required for their resolution.

Several correlations among Gamma Ray Burst (GRB) observables with available redshifts have been recently identified. Proper evaluation and calibration of these correlations may facilitate the use of GRBs as standard candles constraining the expansion history of the universe up to redshifts of $z>6$. Here we use the 69 GRB dataset recently compiled by Schaefer (astro-ph/0612285) and we test the calibration of five of the above correlations ($1:\epkk-E_\gamma$, $2:\epkk-L$, $3:\tlag-L$, $4:V-L$, $5:\trt-L$) with respect to two potential sources of systematics: Evolution with redshift and cosmological model used in the calibration. In examining the model dependence we assume flat \lcdm and vary $\omm$. Our approach avoids the circularity problem of previous studies since we do not fix $\omm$ to find the correlation parameters. Instead we simultaneously minimize $\chi^2$ with respect to both the log-linear correlation parameters $a$, $b$ and the cosmological parameter $\omm$. We find no statistically significant evidence for redshift dependence of $a$ and $b$ in any of the correlation relations tested. We also find that one of the five correlation relations tested ($\epkk-E_\gamma$) has a significantly lower intrinsic dispersion compared to the other correlations. For this correlation relation, the maximum likelihood method favors the existence of a cosmological constant while the other four correlation favor a flat matter dominated universe $\omm \simeq 1$. Finally, a cross-correlation analysis between the GRBs and SnIa data for various values of $\omm$ has shown that the $E_{peak}-E_\gamma$ relation traces well the SnIa redshift regime. However, even the tightest correlation relation ($E_{peak}-E_\gamma$) provides much weaker constraints on $\omm$ than current SnIa data.

We verify numerically that in the context of general relativity (GR), flat models which have the same $\Omega_{\rm m}$ and CMB shift parameter $R$ but different $H(a)$ and $w(a)$ also have very similar (within less than 8%) growth of perturbations even though the dark energy density evolution is quite different. This provides a direct connection between geometrical and dynamical tests of dark energy and may be used as a cosmological test of general relativity.

We consider a universe with a compact extra dimension and a cosmological constant emerging from a suitable ultraviolet cutoff on the zero point energy of the vacuum. We derive the Casimir force between parallel conducting plates as a function of the following scales: plate separation, radius of the extra dimension and cutoff energy scale. We find that there are critical values of these scales where the Casimir force between the plates changes sign. For the cutoff energy scale required to reproduce the observed value of the cosmological constant, we find that the Casimir force changes sign and becomes repulsive for plate separations less than a critical separation $d_0=0.6mm$, assuming a zero radius of the extra dimension (no extra dimension). This prediction contradicts Casimir experiments which indicate an attractive force down to plate separations of $100nm$. For a non-zero extra dimension radius, the critical separation $d_0$ gets even larger than $0.6mm$ and remains inconsistent with Casimir force experiments. We conclude that with or without the presence of a compact extra dimension, vacuum energy with any suitable cutoff can not play the role of the cosmological constant.

We compare the dark energy model constraints obtained by using recent standard ruler data (Baryon Acoustic Oscillations (BAO) at z=0.2 and z=0.35 and Cosmic Microwave Background (CMB) shift parameters R and l_a) with the corresponding constraints obtained by using recent Type Ia Supernovae (SnIa) standard candle data (ESSENCE+SNLS+HST from Davis et. al.). We find that, even though both classes of data are consistent with LCDM at the 2\sigma level, there is a systematic difference between the two classes of data. In particular, we find that for practically all values of the parameters (\Omega_0m,\Omega_b) in the 2\sigma range of the the 3-year WMAP data (WMAP3) best fit, LCDM is significantly more consistent with the SnIa data than with the CMB+BAO data. For example for (\Omega_0m,\Omega_b)=(0.24,0.042) corresponding to the best fit values of WMAP3, the dark energy equation of state parametrization w(z)=w_0 + w_1 (z/(1+z)) best fit is at a 0.5\sigma distance from LCDM (w_0=-1,w_1=0) using the SnIa data and 1.7\sigma away from LCDM using the CMB+BAO data. There is a similar trend in the earlier data (SNLS vs CMB+BAO at z=0.35). This trend is such that the standard ruler CMB+BAO data show a mild preference for crossing of the phantom divide line w=-1, while the recent SnIa data favor LCDM. Despite of this mild difference in trends, we find no statistically significant evidence for violation of the cosmic distance duality relation \eta \equiv d_L(z)/(d_A(z) (1+z)^2)=1. For example, using a prior of \Omega_0m=0.24, we find \eta=0.95 \pm 0.025 in the redshift range 0<z<2, which is consistent with distance duality at the 2\sigma level.

We have compiled a dataset consisting of 22 datapoints at a redshift range (0.15,3.8) which can be used to constrain the linear perturbation growth rate f=\frac{d\ln\delta}{d\ln a}. Five of these data-points constrain directly the growth rate f through either redshift distortions or change of the power spectrum with redshift. The rest of the datapoints constrain f indirectly through the rms mass fluctuation \sigma_8(z) inferred from Ly-\alpha at various redshifts. Our analysis tests the consistency of the LCDM model and leads to a constraint of the Wang-Steinhardt growth index \gamma (defined from f=\Omega_m^\gamma) as \gamma=0.67^{+0.20}_{-0.17}. This result is clearly consistent at $1\sigma$ with the value \gamma={6/11}=0.55 predicted by LCDM. A first order expansion of the index \gamma in redshift space leads to similar results.We also apply our analysis on a new null test of LCDM which is similar to the one recently proposed by Chiba and Nakamura (arXiv:0708.3877) but does not involve derivatives of the expansion rate $H(z)$. This also leads to the fact that LCDM provides an excellent fit to the current linear growth data.

We consider a braneworld model with bulk-brane energy exchange. This allows for crossing of the w=-1 phantom divide line without introducing phantom energy with quantum instabilities. We use the latest SnIa data included in the Gold06 dataset to provide an estimate of the preferred parameter values of this braneworld model. We use three fitting approaches which provide best fit parameter values and hint towards a bulk energy component that behaves like relativistic matter which is propagating in the bulk and is moving at a speed v along the fifth dimension, while the bulk-brane energy exchange component corresponds to negative pressure and signifies energy flowing from the bulk into the brane. We find that the best fit effective equation of state parameter $w_{eff}$ marginally crosses the phantom divide line w=-1. Thus, we have demonstrated both the ability of this class of braneworld models to provide crossing of the phantom divide and also that cosmological data hint towards natural values for the model parameters.

We consider scalar-tensor theories and reconstruct their potential U(\Phi) and coupling F(\Phi) by demanding a background LCDM cosmology. In particular we impose a background cosmic history H(z) provided by the usual flat LCDM parameterization through the radiation (w_{eff}=1/3), matter (w_{eff}=0) and deSitter (w_{eff}=-1) eras. The cosmological dynamical system which is constrained to obey the LCDM cosmic history presents five critical points in each era, one of which corresponding to the standard General Relativity (GR). In the cases that differ from GR, the reconstructed coupling and potential are of the form F(\Phi)\sim \Phi^2 and U(\Phi)\sim F(\Phi)^m where m is a constant. This class of scalar tensor theories is also theoretically motivated by a completely independent approach: imposing maximal Noether symmetry on the scalar-tensor Lagrangian. This approach provides independently: i) the form of the coupling and the potential as F(\Phi)\sim \Phi^2 and U(\Phi)\sim F(\Phi)^m, ii) a conserved charge related to the potential and the coupling and iii) allows the derivation of exact solutions by first integrals of motion.

We consider f(R) modified gravity theories in the metric variation formalism and attempt to reconstruct the function f(R) by demanding a background LCDM cosmology. In particular we impose the following requirements: a. A background cosmic history H(z) provided by the usual flat LCDM parametrization though the radiation (w_eff=1/3), matter (w_eff=0) and deSitter (w_eff=-1) eras. b. Matter and radiation dominate during the `matter' and `radiation' eras respectively i.e. \Omega_m =1 when w_eff=0 and \Omega_r=1 when w_eff=1/3. We have found that the cosmological dynamical system constrained to obey the LCDM cosmic history has four critical points in each era which correspondingly lead to four forms of f(R). One of them is the usual general relativistic form f(R)=R-2\Lambda. The other three forms in each era, reproduce the LCDM cosmic history but they do not satisfy requirement b. stated above. Only one of these forms (different from general relativity) is found to be an attractor of the dynamical cosmological evolution. It has (\Omega_DE=1, \Omega_r=0, \Omega_m=0) throughout the evolution. Its phase space trajectory is numerically obtained.

The Gold06 SnIa dataset recently released in astro-ph/0611572 consists of five distinct subsets defined by the group or instrument that discovered and analyzed the corresponding data. These subsets are: the SNLS subset (47 SnIa), the HST subset (30 SnIa), the HZSST subset (41 SnIa), the SCP subset (26 SnIa) and the Low Redshift (LR) subset (38 SnIa). These subsets sum up to the 182 SnIa of the Gold06 dataset. We use Monte-Carlo simulations to study the statistical consistency of each one of the above subsets with the full Gold06 dataset. In particular, we compare the best fit $w(z)$ parameters (w_0,w_1) obtained by subtracting each one of the above subsets from the Gold06 dataset (subset truncation), with the corresponding best fit parameters (w^r_0,w^r_1) obtained by subtracting the same number of randomly selected SnIa from the same redshift range of the Gold06 dataset (random truncation). We find that the probability for (w^r_0,w^r_1)=(w_0,w_1) is large for the Gold06 minus SCP (Gold06-SCP) truncation but is less than 5% for the Gold06-SNLS, Gold06-HZSST and Gold06-HST truncations. This result implies that the Gold06 dataset is not statistically homogeneous. By comparing the values of the best fit (w_0,w_1) for each subset truncation we find that the tension among subsets is such that the SNLS and HST subsets are statistically consistent with each other and `pull' towards LCDM (w_0=-1,w_1=0) while the HZSST subset is statistically distinct and strongly `pulls' towards a varying w(z) crossing the line $w=-1$ from below (w_0<-1,w_1>0). We also isolate six SnIa that are mostly responsible for this behavior of the HZSST subset.

We use a low redshift expansion of the cosmological equations of extended (scalar-tensor) quintessence to divide the observable Hubble history parameter space in four sectors: A forbidden sector I where the scalar field of the theory becomes imaginary (the kinetic term becomes negative), a forbidden sector II where the scalar field rolls up (instead of down) its potential, an allowed `freezing' quintessence sector III where the scalar field is currently decelerating down its potential towards freezing and an allowed `thawing' sector IV where the scalar field is currently accelerating down its potential. The dividing lines between the sectors depend sensitively on the time derivatives of the Newton's constant G over powers of the Hubble parameter. For minimally coupled quintessence which appears as a special case for a constant G our results are consistent with previous studies. Observable parameter \chi^2 contours based on current data (SNLS dataset) are also constructed on top of the sectors, for a prior of \Omega_m=0.24. By demanding that the observed 2\sigma \chi^2 parameter contours do not lie entirely in the forbidden sectors we derive stringent constraints on the current second time derivative of Newton's constant G. In particular we find {\ddot G}/G >-1.91 H_0^2=-2 10^{-20}h^2 yrs^{-2} at the 2\sigma level which is complementary to solar system tests which constrain only the first derivative of G as |{\dot G}/G|<10^{-14}yrs^{-1} at 1\sigma.

If the dark energy equation of state parameter w(z) crosses the phantom divide line w=-1 (or equivalently if the expression d(H^2(z))/dz - 3\Omega_m H_0^2 (1+z)^2 changes sign) at recent redshifts, then there are two possible cosmological implications: Either the dark energy consists of multiple components with at least one non-canonical phantom component or general relativity needs to be extended to a more general theory on cosmological scales. The former possibility requires the existence of a phantom component which has been shown to suffer from serious theoretical problems and instabilities. Therefore, the later possibility is the simplest realistic theoretical framework in which such a crossing can be realized. After providing a pedagogical description of various dark energy observational probes, we use a set of such probes (including the Gold SnIa sample, the first year SNLS dataset, the 3-year WMAP CMB shift parameter, the SDSS baryon acoustic oscillations peak (BAO), the X-ray gas mass fraction in clusters and the linear growth rate of perturbations at z=0.15 as obtained from the 2dF galaxy redshift survey) to investigate the priors required for cosmological observations to favor crossing of the phantom divide. We find that a low \Omega_m prior (0.2<\Omega_m <0.25) leads, for most observational probes (except of the SNLS data), to an increased probability (mild trend) for phantom divide crossing. An interesting degeneracy of the ISW effect in the CMB perturbation spectrum is also pointed out.

If Newton's constant G evolves on cosmological timescales as predicted by extended gravity theories then Type Ia supernovae (SnIa) can not be treated as standard candles. The magnitude-redshift datasets however can still be useful. They can be used to simultaneously fit for both H(z) and G(z) (so that local G(z) constraints are also satisfied) in the context of appropriate parametrizations. Here we demonstrate how can this analysis be done by applying it to the Gold SnIa dataset. We compare the derived effective equation of state parameter w(z) at best fit with the corresponding result obtained by neglecting the evolution G(z). We show that even though the results clearly differ from each other, in both cases the best fit w(z) crosses the phantom divide w=-1. We then attempt to reconstruct a scalar tensor theory that predicts the derived best fit forms of H(z) and G(z). Since the best fit G(z) fixes the scalar tensor potential evolution F(z), there is no ambiguity in the reconstruction and the potential U(z) can be derived uniquely. The particular reconstructed scalar tensor theory however, involves a change of sign of the kinetic term $\Phi'(z)^2$ as in the minimally coupled case.

This is a pedagogical review of the recent observational data obtained from type Ia supernova surveys that support the accelerating expansion of the universe. The methods for the analysis of the data are reviewed and some of the theoretical implications obtained from their analysis are discussed.

We have performed a comparative analysis of three recent and reliable SnIa datasets available in the literature: the Full Gold (FG) dataset (157 data points $0<z<1.75$), a Truncated Gold (TG) dataset (140 data points $0<z<1$) and the most recent Supernova Legacy Survey (SNLS) dataset (115 data points $0<z<1$). We have analyzed and compared the likelihood of cosmological constant and dynamical dark energy parametrizations allowing for crossing of the phantom divide line (PDL). We find that even though the constraints obtained using the three datasets are consistent with each other at the 95% confidence level, the latest (SNLS) dataset shows distinct trends which are not shared by the Gold datasets. We find that the best fit dynamical $w(z)$ obtained from the SNLS dataset does not cross the PDL $w=-1$ and remains above and close to the $w=-1$ line for the whole redshift range $0<z<1$ showing no evidence for phantom behavior. The LCDM parameter values ($w_0=-1$, $w_1=0$) almost coincide with the best fit parameters of the dynamical $w(z)$ parametrizations. In contrast, the best fit dynamical $w(z)$ obtained from the Gold datasets (FG and TG) clearly crosses the PDL and departs significantly from the PDL $w=-1$ line while the LCDM parameter values are about $2\sigma$ away from the best fit $w(z)$. In addition, the $(\Omega_{0m},\Omega_\Lambda)$ parameters in a LCDM parametrization without a flat prior, fit by the SNLS dataset, favor the minimal flat LCDM concordance model. The corresponding fit with the Gold datasets mildly favors a closed universe and the flat LCDM parameter values are $1\sigma - 2\sigma$ away from the best fit $(\Omega_{0m},\Omega_\Lambda)$.

We use the SnIa Gold dataset to compare LCDM with 10 representative parametrizations of the recent Hubble expansion history $H(z)$. For the comparison we use two statistical tests; the usual $\chi_{min}^2$ which is insensitive to the parametrization number of parameters, and a statistic we call the p-test which depends on both the value of $\chi_{min}^2$ and the number $n$ of the parametrization parameters. The p-test measures the confidence level to which the parameter values corresponding to LCDM are excluded from the viewpoint of the parametrization tested. For example, for a linear equation of state parametrization $w(z)=w_0 + w_1 z$ the LCDM parameter values ($w_0=-1$, $w_1=0$) are excluded at 75% confidence level. We use a flat prior and $\Omega_{0m}=0.3$. All parametrizations tested are consistent with the Gold dataset at their best fit. According to both statistical tests, the worst fits among the 10 parametrizations, correspond to the Chaplygin gas, the brane world and the Cardassian parametrizations. The best fit is achieved by oscillating parametrizations which can exclude the parameter values corresponding to LCDM at 85% confidence level. Even though this level of significance does not provide a statistically significant exclusion of LCDM (it is less than $2\sigma$) and does not by itself constitute conclusive evidence for oscillations in the cosmological expansion, when combined with similar independent recent evidence for oscillations coming from the CMB and matter power spectra it becomes an issue worth of further investigation.

There is accumulating observational evidence (based on SnIa data) that the dark energy equation of state parameter $w$ may be evolving with time and crossing the phantom divide barrier $w=-1$ at recent times. The confirmation of these indications by future data would indicate that minimally coupled quintessence can not reproduce the observed expansion rate $H(z)$ for any scalar field potential. Here we explicitly demonstrate that scalar tensor theories of gravity (extended quintessence) can predict crossing of the phantom divide barrier. We reconstruct phenomenologically viable scalar-tensor potentials $F(\Phi)$ and $U(\Phi)$ that can reproduce a polynomial best fit expansion rate $H(z)$ and the corresponding dark energy equation of state parameter $w(z)$ crossing the $w=-1$ line. The form of the reconstructed scalar tensor potentials is severely constrained but is not uniquely determined. This is due to observational uncertainties in the form of $H(z)$ and also because a single observed function $H(z)$ does not suffice for the reconstruction of the two potential functions $F(\Phi)$ and $U(\Phi)$.

We study quintessence and phantom field theory models based on linear-negative potentials of the form $V(\phi)=s \phi$. We investigate the predicted redshift dependence of the equation of state parameter $w(z$ for a wide range of slopes $s$ in both quintessence and phantom models. We use the gold dataset of 157 SnIa and place constraints on the allowed range of slopes $s$. We find $s=0\pm 1.6$ for quintessence and $s=\pm 0.7\pm 1$ for phantom models (the range is at the $2\sigma$ level and the units of $s$ are in $\sqrt{3}M_p H_0^2\simeq 10^{-38}eV^3$ where $M_p$ is the Planck mass). In both cases the best fit is very close to $s\simeq 0$ corresponding to a cosmological constant. We also show that specific model independent parametrizations of $w(z)$ which allow crossing of the phantom divide line $w=-1$ (hereafter PDL) provide significantly better fits to the data. Unfortunately such crossings are not allowed in any phantom or quintessence single field model minimally coupled to gravity. Mixed models (coupled phantom-quintessence fields) can in principle lead to a $w(z)$ crossing the PDL but a preliminary investigation indicates that this does not happen for natural initial conditions.

The cosmic string theory for cosmological fluctuations is a good example of healthy scientific progress in cosmology. It is a well defined physically motivated model that has been tested by cosmological observations and has been ruled out as a primary source of primordial fluctuations. Until about fifteen years ago, the cosmic string theory of cosmological perturbations provided one of the two physically motivated candidate theories for the generation of primordial perturbations. The cosmological data that appeared during the last decade have been compared with the well defined predictions of the theory and have ruled out cosmic strings as a primary source of primordial cosmological perturbations. Since cosmic strings are predicted to form after inflation in a wide range of microphysical theories including (supersymmetric and fundamental string theories) their observational bounds may serve a source of serious constraints for these theories. This is a pedagogical review of the historical development, the main predictions of the cosmic string theory and the constraints that have been imposed on it by cosmological observations. Recent lensing events that could be attributed to lighter cosmic strings are also discussed.

We study analytically and numerically the evolution of bound systems in universes with accelerating expansion where the acceleration either increases with time towards a Big Rip singularity (phantom cosmologies) or decreases with time (quintessence). We confirm the finding of Caldwell et. al. [arXiv:astro-ph/0302506] (hereafter CKW) that bound structures get dissociated in phantom cosmologies but we demonstrate that this happens earlier than anticipated by CKW. In particular we find that the `rip time' when a bound system gets unbounded is not the time when the repulsive phantom energy gravitational potential due to the average $(\rho + 3p)$ balances the attractive gravitational potential of the mass M of the system. Instead, the `rip time' is the time when the minimum of the time dependent effective potential (including the centrifugal term) disappears. For the Milky Way galaxy this happens approximately 180Myrs before the Big Rip singularity instead of approximately 60Myrs indicated by CKW for a phantom cosmology with w=-1.5. A numerical reconstruction of the dissociating bound orbits is presented.

We study the expansion history of the universe up to a redshift of z=1.75 using the 194 recently published SnIa data by Tonry et. al. and Barris et. al. In particular we find the best fit forms of several cosmological models and $H(z)$ ansatze, determine the best fit values of their parameters and rank them according to increasing value of $\chi_{min}^2$ (the minimum value of $\chi^2$ for each $H(z)$ ansatz). We use a prior of $\Omega_{0m} = 0.3$ and assume flat geometry of the universe. No prior assumptions are made about validity of energy conditions. The fitted models are fourteen and include SCDM, LCDM, dark energy with constant equation of state parameter $w$ (quiessence), third order polynomial for $H(1+z)$, Chaplygin gas, Cardassian model, $w(z)=w_0 + w_1 z$, $w(z)=w_0 + z w_1/(1+z)$, an oscillating ansatz for $H(z)$ etc. All these models with the exception of SCDM are consistent with the present data. However, the quality of the fit differs significantly among them and so do the predicted forms of $w(z)$ and $H(z)$ at best fit. The worst fit among the data-consistent models considered corresponds to the simplest model LCDM ($\chi_{min}^2 = 198.7$ for $\Omega_{0m} = 0.34$) while the best fit is achieved by the three parameter oscillating ansatz ($\chi_{min}^2 = 194.1$). Most of the best fit ansatze have an equation of state parameter $w(z)$ that varies between $w(z) \simeq -1$ for $z<0.5$ to $w(z) > 0$ for $z>1$. This implies that the sign of the pressure of the dark energy may be alternating as the redshift increases. The goodness of fit of the oscillating $H(z)$ ansatz lends further support to this possibility.

We study analytically and numerically the core structure of global vortices forming on topologically deformed brane-worlds with a single toroidally compact extra dimension. It is shown that for an extra dimension size larger than the scale of symmetry breaking the magnitude of the complex scalar field at the vortex center can dynamically remain non-zero. Singlevaluedness and regularity are not violated. Instead, the winding escapes to the extra dimension at the vortex center. As the extra dimension size decreases the field magnitude at the core dynamically decreases also and in the limit of zero extra dimension size we reobtain the familiar global vortex solution. Extensions to other types of defects and gauged symmetries are also discussed.

We consider a Brans-Dicke scalar field stabilized by a general power law potential with power index $n$ at a finite equilibrium value. Redshifting matter induces oscillations of the scalar field around its equilibrium due to the scalar field coupling to the trace of the energy momentum tensor. If the stabilizing potential is sufficiently steep these high frequency oscillations are consistent with observational and experimental constraints for arbitrary value of the Brans-Dicke parameter $\omega$. We study analytically and numerically the equation of state of these high frequency oscillations in terms of the parameters $\omega$ and $n$ and find the corresponding evolution of the universe scale factor. We find that the equation of state parameter can be negative and less than -1 but it is not related to the evolution of the scale factor in the usual way. Nevertheless, accelerating expansion is found for a certain parameter range. Our analysis applies also to oscillations of the size of extra dimensions (the radion field) around an equilibrium value. This duality between self-coupled Brans-Dicke and radion dynamics is applicable for $\omega= -1 + 1/D$ where D is the number of extra dimensions.

We show that the redshift of pressureless matter density due to the expansion of the universe generically induces small oscillations in the stabilized radius of extra dimensions (the radion field). The frequency of these oscillations is proportional to the mass of the radion and can have interesting cosmological consequences. For very low radion masses $m_b$ ($m_b\sim10-100 H_0\simeq10^{-32} eV$) these low frequency oscillations lead to oscillations in the expansion rate of the universe. The occurrence of acceleration periods could naturally lead to a resolution of the coincidence problem, without need of dark energy. Even though this scenario for low radion mass is consistent with several observational tests it has difficulty to meet fifth force constraints. If viewed as an effective Brans-Dicke theory it predicts $\omega=-1+\frac{1}{D}$ ($D$ is the number of extra dimensions), while experiments on scales larger than $1mm$ imply $\omega>2500$. By deriving the generalized Newtonian potential corresponding to a massive toroidally compact radion we demonstrate that Newtonian gravity is modified only on scales smaller than $m_b^{-1}$. Thus, these constraints do not apply for $m_b>10^{-3} eV$ (high frequency oscillations) corresponding to scales less than the current experiments ($0.3mm$). Even though these high frequency oscillations can not resolve the coincidence problem they provide a natural mechanism for dark matter generation. This type of dark matter has many similarities with the axion.

The linearity and quietness of the Local ($< 10 Mpc$) Hubble Flow (LHF) in view of the very clumpy local universe is a long standing puzzle in standard and in open CDM cosmogony. The question addressed in this paper is whether the antigravity component of the recently discovered dark energy can cool the velocity flow enough to provide a solution to this puzzle. We calculate the growth of matter fluctuations in a flat universe containing a fraction $\Omega_X(t_0)$ of dark energy obeying the time independent equation of state $p_X = w \rho_X$. We find that dark energy can indeed cool the LHF. However the dark energy parameter values required to make the predicted velocity dispersion consistent with the observed value $v_{rms}\simeq 40km/sec$ have been ruled out by other observational tests constraining the dark energy parameters $w$ and $\Omega_X$. Therefore despite the claims of recent qualitative studies dark energy with time independent equation of state can not by itself explain the quietness and linearity of the Local Hubble Flow.

We investigate the effects of negative pressure induced by dark energy (cosmological constant or quintessence) on the dynamics at various astrophysical scales. Negative pressure induces a repulsive term (antigravity) in Newton's law which dominates on large scales. Assuming a value of the cosmological constant consistent with the recent SnIa data we determine the critical scale $r_c$ beyond which antigravity dominates the dynamics ($r_c \sim 1Mpc $) and discuss some of the dynamical effects implied. We show that dynamically induced mass estimates on the scale of the Local Group and beyond are significantly modified due to negative pressure. We also briefly discuss possible dynamical tests (eg effects on local Hubble flow) that can be applied on relatively small scales (a few $Mpc$) to determine the density and equation of state of dark energy.

We use the publicly available code CMBFAST as modified by L. Pogosian and T. Vachaspati to simulate the effects of wiggly cosmic strings on the Cosmic Microwave Background (CMB). Using the modified CMBFAST code which takes into account vector modes and models wiggly cosmic strings by the one scale model, we go beyond the angular power spectrum to construct CMB temperature maps with resolution of a few degrees. The statistics of these maps are then studied using conventional and recently proposed statistical tests optimized for the detection of hidden temperature discontinuities induced by the Gott-Kaiser-Stebbins effect. We show however, that these realistic maps can not be distinguished in a statistically significant way from purely Gaussian maps with identical power spectrum.

Newton's law gets modified in the presence of a cosmological constant by a small repulsive term (antigarvity) that is proportional to the distance. Assuming a value of the cosmological constant consistent with the recent SnIa data ($\Lambda \simeq 10^{-52} m^{-2}$) we investigate the significance of this term on various astrophysical scales. We find that on galactic scales or smaller (less than a few tens of kpc) the dynamical effects of the vacuum energy are negligible by several orders of magnitude. On scales of 1Mpc or larger however we find that vacuum energy can significantly affect the dynamics. For example we show that the velocity data in the Local Group of galaxies correspond to galactic masses increased by 35% in the presence of vacuum energy. The effect is even more important on larger low density systems like clusters of galaxies or superclusters.

We present a detailed analytical and numerical study of a novel type of static, superconducting, classically stable string texture in a renormalizable topologically trivial massive U(1) gauge model with one charged and one neutral scalar. An upper bound on the mass of the charged scalar as well as on the current that the string can carry are established. A preliminary unsuccesful search for stable solutions corresponding to large superconducting loops is also reported.

We apply new statistical tests on the four year 53 GHz DMR data and show that there is significant probability for the existence of coherent temperature discontinuities hidden in the CMB maps around the Galactic Poles. Comparing Monte Carlo simulations with the DMR maps we find that the probability for the existence of a coherent discontinuity with amplitude $0.2 \times ({{\delta T} \over T})_{rms}$ superposed on the gaussian fluctuations is more than double than the corresponding probability when no discontinuity is present. This result is consistent with the existence of a horizon size long cosmic string with mass per unit length $\mu$ given by $G\mu v_s \gamma_s \simeq 10^{-6}/\pi$ where $v_s$ is the string velocity and $\gamma_s$ is the corresponding Lorentz factor.

Using analytical methods and Monte Carlo simulations, we analyze new statistics designed to detect isolated step-like discontinuities which are coherent over large areas of Cosmic Microwave Background (CMB) pixel maps. Such coherent temperature discontinuities are predicted by the Kaiser-Stebbins effect to form due to long cosmic strings present in our present horizon. The background of the coherent step-like seed is assumed to be a scale invariant Gaussian random field which could have been produced by a superposition of seeds on smaller scales and/or by inflationary quantum fluctuations. We find that the proposed statistics can detect the presense of a coherent discontinuity at a sensitivity level almost an order of magnitude better compared to more conventional statistics like the skewness or the kurtosis.

We demonstrate the existence of global monopole and vortex configurations whose core exhibits a phase structure. We determine the critical values of parameters for which the transition from the symmetric to the non-symmetric phase occurs and discuss the novel dynamics implied by the non-symmetric cores for defect interactions. We model phase transitions in the core of global embedded topological defects by identifying the relevant parameters with the vacuum expectation value of a dynamical scalar field. Finally, we argue that superheavy defects that undergo a core phase transition in the very early universe provide a novel realization for topological inflation.

We review recent work on a new class of topological defects which possess a nonsymmetric core. They arise in scalar field theories with global symmetries, U(1) for domain walls and SU(2) for vortices, which are explicitly broken to $Z_2$ and U(1) respectively. Both of the latter symmetries are spontaneously broken. For a particular range of parameters both types of defect solutions are shown to become unstable and decay to the well known stable walls and vortices with symmetric cores.

A new statistic is briefly reviewed, designed to detect isolated coherent step-like discontinuities produced by cosmic strings present at late times. As a background I superpose a scale invariant Gaussian random field which could have been produced by a superposition of seeds on all scales and/or by inflationary quantum fluctuations. The statistical variable considered is the Sample Mean Difference (SMD) between large neighbouring sectors of CMB maps, separated by lines in two dimensional maps and points in one dimensional maps. I find that the SMD statistics can detect at the $1\sigma$ level the presense of a long string with $G\mu (v_s \gamma_s)={1\over {8\pi}}({{\delta T}\over T})_{rms} \simeq 0.5 \times 10^{-7}$ while more conventional statistics like the skewness or the kurtosis require a value of $G\mu$ almost an order of magnitude larger for detectability at a comparable level.

Using analytical methods and Monte Carlo simulations, we analyze a new statistic designed to detect isolated step-like discontinuities which are coherent over large areas of Cosmic Microwave Background (CMB) pixel maps. Such coherent temperature discontinuities are predicted by the Kaiser-Stebbins effect to form due to long cosmic strings present in our present horizon. The background of the coherent step-like seed is assumed to be a scale invariant Gaussian random field which could have been produced by a superposition of seeds on smaller scales and/or by inflationary quantum fluctuations. The effects of uncorrelated Gaussian random noise are also studied. The statistical variable considered is the Sample Mean Difference (SMD) between large neighbouring sectors of CMB maps, separated by a straight line in two dimensional maps and a point in one dimensional maps. We find that including noise, the SMD statistics can detect at the $1 \sigma$ to $2 \sigma$ level the presense of a long string with $G\mu (v_s \gamma_s)= 1/(8\pi) ((\delta T)/T)_{rms} \simeq 0.5 \times 10^{-7}$ while more conventional statistics like the skewness or the kurtosis require a value of $G\mu$ almost an order of magnitude larger for detectability at a comparable level.

We demonstrate that field theories involving explicit breaking of continous symmetries, incorporate two generic classes of topological defects each of which is stable for a particular range of parameters. The first class includes defects of the usual type where the symmetry gets restored in the core and vacuum energy gets trapped there. We show however that these defect solutions become unstable for certain ranges of parameters and decay not to the vacuum but to another type of stable defect where the symmetry in not restored in the core. In the wall case, initially spherical, bubble-like configurations are simulated numerically and shown to evolve generically towards a planar collapse. In the string case, the decay of the symmetric core vortex resembles the decay of a semilocal string to a skyrmion with the important difference that while the skyrmion is unstable and decays to the vacuum, the resulting non-symmetric vortex is topologically stable.

Using Monte Carlo simulations of perturbations induced by cosmic strings on the microwave background, we demonstrate the scale invariance of string fluctuation patterns. By comparing string-induced fluctuation patterns with gaussian random phase ones, we show that the non-gaussian signatures of the string patterns are detectable by tests based on the moments of the distributions only for angular scales smaller than a few arcminutes and for maps based on the gradient of temperature fluctuations. However, we find that tests of the gaussianity of the moments fail when we include a reasonable amount of instrumental noise in a pattern. Signal to noise ratios of $3.3$ or greater completely suppress a string pattern's non-gaussian features even at the highest resolutions.

Using a toy model Lagrangian we investigate the formation of vortices in first order phase transitions. The evolution and interactions of vacuum bubbles are also studied using both analytical approximations and a numerical simulation of scalar field dynamics. A long lived bubble wall bound state is discovered and its existence is justified by using a simplified potential for the bubble wall interaction. The conditions that need to be satisfied for vortex formation by bubble collisions are also studied with particular emphasis placed on geometrical considerations. These conditions are then implemented in a Monte Carlo simulation for the study of the probability of defect formation. It is shown that the probability of vortex formation by collision of relativistically expanding bubbles gets reduced by about 10% due to the above mentioned geometric effects.

This is a pedagogical introduction to the basics of the cosmic string theory and also a review of the recent progress made with respect to the macrophysical predictions of the theory. Topics covered include, string formation and evolution, large scale structure formation, generation of peculiar velocity flows, cosmic microwave background (CMB) fluctuations, lensing, gravitational waves and constraints from the msec pulsar. Particular emphasis is placed on the signatures predicted on the CMB and the corresponding non-gaussian features.

By using a simple analytical model based on counting random multiple impulses inflicted on photons by a network of cosmic strings we show how to construct the general q-point temperature correlation function of the Cosmic Microwave Background radiation. Our analysis is sensible specially for large angular scales where the Kaiser-Stebbins effect is dominant. Then we concentrate our study on the four-point function and in particular on its zero-lag limit, namely, the excess kurtosis parameter, for which we obtain a predicted value of $\sim 10^{-2}$. In addition, we estimate the cosmic variance for the kurtosis due to a Gaussian fluctuation field, showing its dependence on the primordial spectral index of density fluctuations $n$ and finding agreement with previous published results for the particular case of a flat Harrison-Zel'dovich spectrum. Our value for the kurtosis compares well with previous analyses but falls below the threshold imposed by the cosmic variance when commonly accepted parameters from string simulations are considered. In particular the non-Gaussian signal is found to be inversely proportional to the scaling number of defects, as could be expected by the central limit theorem.

I show the existence of a new type of vortex solution which is non-static but stationary and carries angular momentum. This {\it spinning vortex} can be embedded in models with trivial vacuum topology like a model with $SU(2)_{global}\times U(1)_{local} \rightarrow U(1)_{global}$ symmetry breaking. The stability properties of the embedded spinning vortex are also studied in detail and it is shown that stability improves drastically as angular momentum increases. The implications of this result for vortices embedded in the standard electroweak model are under study.

I review recent progress towards a detailed understanding of the power spectrum of cosmic string perturbations induced on the Microwave Background (CMB). The use of a simple analytic model allows to include important effects that have not been included in previous studies. These effects are potential fluctuations on the last scattering surface (the Sachs-Wolfe effect) and Doppler CMB perturbations. These previously neglected fluctuations are shown to dominate over post-recombination fluctuations (Gott-Kaiser-Stebbins effect) on angular scales less than about $2^\circ$ assuming that no reionization occurs. The effective power spectrum index $n_{eff}$ is calculated on COBE scales and is shown to be somewhat larger than 1 ($1.35 \lsim n_{eff} \lsim 1.5$). After normalizing the model on COBE data and string simulations, I derive the predicted values of $({{\delta T}\over T})_{rms}$ for several ongoing CMB experiments on medium and small angular scales. A comparison is also made with current observations and predictions of models based on inflation.

The spinning vortex is a stationary generalization of the Nielsen-Olesen vortex involving a linear time dependence of the Goldstone boson. Here we show that this vortex can be embedded in models with $SU(2)_{global}\times U(1)_{local}$ symmetry. We also map the stability sector in parameter space and show that for sufficiently large spinning velocities the vortex is stable for {\it any} value of the Higgs self coupling parameter $\beta$. This is a significant improvement of stability compared to the semilocal vortex introduced by Vachaspati and Achucarro which is stable only for $\beta<1$. This result may have significant implications for electroweak vortices.

Using a simple analytic model based on seed superposition we obtain the spectrum of microwave background perturbations induced by cosmic strings on all angular scales larger than about 2 armin. We assume standard recombination in an Einstein de Sitter universe with $h=1/2$ and study the fluctuation spectrum along a great circle in the sky. Doppler and potential perturbations on the last scattering surface (Sachs-Wolfe effect) are shown to dominate over post-recombination perturbations on scales smaller than about 2 degrees. Using a filter function corresponding to the COBE experiment we obtain an effective power spectrum index $n_{eff}\simeq 1.35$ in good agreement with the recently announced second year COBE data showing $n_{eff}\simeq 1.5$. The only free parameter ($G\mu$) of the string model is fixed by normalizing on the COBE detection leading to $G\mu=1.6$. Other parameters (e.g. the {\it rms} string velocity) are fixed by comparing with string simulations. Using these parameter values we compare the {\it rms} temperature fluctuations $({{\Delta T}\over T})_{\it rms}$ predicted for ongoing experiments (Tenerife, SP91, MAX, OVRO etc) with detections and with the corresponding predictions of inflationary models.

ABSRACT: We review recent progress on testing the hypothesis of the existence of cosmic string perturbations in microwave background maps. Using an analytical model for the string network we show that the predicted amplitude and spectrum of MBR fluctuations are consistent with the COBE data for a reasonable value of the single free parameter of the string model (the mass per unit length of the string $\mu$). To {\it distinguish} the predictions of cosmic strings from those of Gaussian models it is necessary to apply specific {\it statistical tests} which are sensitive to non-random phases of the MBR temperature maps. We discuss two such statistical tests: First, the probability distribution and the moments of fluctuations ${{\delta T}\over T}$ and second the corresponding statistical quantities for the {\it gradient} of fluctuations along a fixed direction. We show that the second statistic can detect cosmic string specific signatures on MBR temperature maps with resolution of a few arcminutes or smaller.

We investigate analytically and numerically the asymptotic behavior of the Nielsen-Olesen vortex solutions and show that they approach their asymptotic values exponentially but with exponents that differ from the ones quoted in the literature. In particular, it is shown that for values of the Higgs self-coupling that are larger than a critical value, both the Higgs field and the gauge field approach exponentially their asymptotic values with {\it equal} exponents that are {\it independent} of the Higgs self-coupling.

Using an analytical model for the string network we show that the kurtosis of cosmic microwave background (CMB) temperature gradient maps is a good statistic to distinguish between the cosmic string model and inflationary models of structure formation. The difference between the stringy and inflationary value for the kurtosis is inversely proportional to the angular resolution and to the number of strings per Hubble volume of the strings' scaling solution. If strings are indeed responsible for CMB anisotropies then experiments with resolutions of a couple of arcminutes or smaller could determine it using this statistic.

We propose a new test for distinguishing observationally cosmological models based on seed-like primordial perturbations (like cosmic strings or textures), from models based on Gaussian fluctuations. We investigate analytically the {\it Fourier space} statistical properties of temperature or density fluctuation patterns generated by seed-like objects and compare these properties with those of Gaussian fluctuations generated during inflation. We show that the proposed statistical test can easily identify temperature fluctuations produced by a superposition of a small number of seeds per horizon scale for {\it any} observational angular resolution and {\it any} seed geometry. However, due to the Central Limit Theorem, the distinction becomes more difficult as the number of seeds in the fluctuation pattern increases.

We show analytically and numerically the existence of double vortex solutions in two-Higgs systems. These solutions are generalizations of the \no vortices and exist for all values of the parameters in the Lagrangians considered. We derive analytically the asymptotic behavior of the solutions and confirm it numerically by solving the field equations. Finally, we show that these solutions can be embedded in realistic theories like the two-doublet extension of the standard model.

Using an analytic model we show that the predictions of the cosmic string model for the peculiar velocities and the microwave background (MBR) anisotropy depend on similar combinations of string evolution parameters. Normalizing from the COBE detection of MBR anisotropy, and for certain reasonable values of string network evolution parameters, we find that the magnitude of predicted velocity flows is in good agreement with observations on small scales but is inconsistent with observations on large scales ($ > 50 h^{-1}$Mpc). The Cosmic Mach Number obtained from the cosmic string scenario is found to depend on a single network evolution parameter and is consistent with observations on scales $5h^{-1}$Mpc to $20h^{-1}$Mpc.

We use the analytical model recently introduced in Ref. \cite{lp92}, to investigate the statistics of temperature fluctuations on the cosmic microwave background (CMB), induced by topological defects. The cases of cosmic strings and textures are studied. We derive analytically the characteristic function of the probability distribution for ${\delta T}\over T$ and use it to obtain the lowest twelve moments including the skewness and the kurtosis. The distribution function is also obtained and it is compared with the Gaussian distribution thus identifying long non-Gaussian tails. We show that for both cosmic strings and textures all odd moments (including skewness) vanish while the relative deviation from the Gaussian for even moments increases with the order of the moment. The non-Gaussian signatures of textures, derived from the distribution function and the moments, are found to be much more prominent than the corresponding signatures for strings. We discuss the physical origin of this result.

We construct a simple analytical model to study the effects of cosmic strings on the microwave background radiation. Our model is based on counting random multiple impulses inflicted on photon trajectories by the string network between the time of recombination and today. We construct the temperature auto-correlation function and use it to obtain the effective power spectrum index n, the rms-quadrupole-normalized amplitude $Q_{rms-PS}$ and the rms temperature variation smoothed on small angular scales. For the values of the scaling solution parameters obtained in Refs.\cite{bb90},\cite{as90} we obtain $n=1.14 \pm 0.5$, $Q_{rms-PS}=(4.5\pm 1.5) G\mu$ and $({{\Delta T}\over T})_{rms}=5.5 G\mu$. Demanding consistency of these results with the COBE data leads to $G\mu=(1.7 \pm 0.7)\times 10^{-6}$ (where $\mu$ is the string mass per unit length), in good agreement with direct normalizations of $\mu$ from observations.