We discuss the production of gravity waves from the fragmentation of a supersymmetric condensate in the early universe. Supersymmetry predicts the existence of flat directions in the potential. At the end of inflation, the scalar fields develop large time-dependent vacuum expectation values along these flat directions. Under some general conditions, the scalar condensates undergo a fragmentation into non-topological solitons, Q-balls. We study this process numerically and confirm the recent analytical calculations showing that it can produce gravity waves observable by Advanced Laser Interferometer Gravitational-Wave Observatory (LIGO), Laser Interferometer Space Antenna (LISA), and Big Bang Observer (BBO). The fragmentation can generate gravity waves with an amplitude as large as Omega_{GW}~10^{-11} and with a peak frequency ranging from mHz to 10 Hz, depending on the parameters. The discovery of such a relic gravitational background radiation can open a new window on the physics at the high scales, even if supersymmetry is broken well above the electroweak scale.

In string gas cosmology, extra dimensions are stabilised by a gas of strings. In the matter-dominated era, competition between matter pushing the extra dimensions to expand and the string gas pulling them back can lead to oscillations of the extra dimensions and acceleration in the visible dimensions. We fit this model to supernova data, taking into account the Big Bang Nucleosynthesis constraint on the energy density of the string gas. The fit to the Union set of supernova data is acceptable, but the fit to the ESSENCE data is poor.

It is plausible that the scalar density perturbations are created by a relatively low scale model of inflation which predicts the CMB anisotropy and excites Standard Model baryon and cold dark matter, but negligible gravity waves. Nevertheless a significantly large tensor perturbations can be observed if there exists a prior phase of high scale inflation separated by a matter or radiation dominated epoch. In this paper we provide a simple example when the gravity waves generated at high scales trickle through the horizon of the second phase of inflation and leave a distinct imprint in the spectrum of the tensor modes. For a high scale inflation occurring at $H\sim 10^{13}$ GeV while the second phase of inflation happening at $H\sim 1$ GeV, the largest tensor to scalar ratio is bounded by $r_{\rm observed}\leq0.8$.

The generation of non-gaussianity is studied in a three fluid curvaton model. By utilizing second order perturbation theory we derive general formulae for the large scale temperature fluctuation and non-gaussianity parameter, $f_{NL}$, that includes the possibility of a non-adiabatic final state. In the adiabatic limit we recover previously known results. The results are applied to a three fluid curvaton model where the curvaton decays into radiation and matter. We find that the amount of non-gaussianity decreases as the final state of the system becomes more adiabatic and that the generated non-gaussianity in the scenario is small, $|f_{NL}| \sim \mathcal{O}(1)$.

Despite the fact that modified theories of gravity, in particular the f(R) gravity models have attracted much attention in the last years, the problem of the energy localization in the framework of these models has not been addressed. In the present work the concept of energy-momentum complexes is presented in this context. We generalize the Landau-Lifshitz prescription of calculating the energy-momentum complex to the framework of f(R) gravity. As an important special case, we explicitly calculate the energy-momentum complex for the Schwarzschild-de Sitter metric for a general f(R) theory as well as for a number of specific, popular choices of f(R).

We consider general metric $f(R)$ theories of gravity by solving the field equations in the presence of a spherical static mass distribution by analytical perturbative means. Expanding the field equations systematically in $\cO(G)$, we solve the resulting set of equations and show that $f(R)$ theories which attempt to solve the dark energy problem very generally lead to $\gamma_{PPN}=1/2$ in the solar system. This excludes a large class of theories as possible explanations of dark energy. We also present the first order correction to $\gamma_{PPN}$ and show that it cannot have a significant effect.

A three fluid system describing the decay of the curvaton is studied by numerical and analytical means. We place constraints on the allowed interaction strengths between the fluids and initial curvaton density by requiring that the curvaton decays before nucleosynthesis while nucleosynthesis, radiation-matter equality and decoupling occur at correct temperatures. We find that with a continuous, time-independent interaction, a small initial curvaton density is naturally preferred along with a low reheating temperature. Allowing for a time-dependent interaction, this constraint can be relaxed. In both cases, a purely adiabatic final state can be generated, but not without fine-tuning. Unlike in the two fluid system, the time-dependent interactions are found to have a small effect on the curvature perturbation itself due to the different nature of the system. The presence of non-gaussianity in the model is discussed.

We study stellar configurations and the space-time around them in metric $f(R)$ theories of gravity. In particular, we focus on the polytropic model of the Sun in the $f(R)=R-\mu^4/R$ model. We show how the stellar configuration in the $f(R)$ theory can, by appropriate initial conditions, be selected to be equal to that described by the Lane-Emden -equation and how a simple scaling relation exists between the solutions. We also derive the correct solution analytically near the center of the star in $f(R)$ theory. Previous analytical and numerical results are confirmed, indicating that the space-time around the Sun is incompatible with Solar System constraints on the properties of gravity. Numerical work shows that stellar configurations, with a regular metric at the center, lead to $\gamma_{PPN}\simeq1/2$ outside the star ie. the Schwarzschild-de Sitter -space-time is not the correct vacuum solution for such configurations. Conversely, by selecting the Schwarzschild-de Sitter -metric as the outside solution, we find that the stellar configuration is unchanged but the metric is irregular at the center. The possibility of constructing a $f(R)$ theory compatible with the Solar System experiments and possible new constraints arising from the radius-mass -relation of stellar objects is discussed.

The so-called shift parameter is related to the position of the first acoustic peak in the power spectrum anisotropies of the cosmic microwave background (CMB) anisotropies. It is an often used quantity in simple tests of dark energy models. However, the shift parameter is not directly measurable from the cosmic microwave background, and its value is usually derived from the data assuming a spatially flat cosmology with dark matter and a cosmological constant. Our aim in this paper is to evaluate the effectiveness of the shift parameter as a constraint on dark energy models, and the potential pitfalls in using it as a test of non-standard models of dark energy. By comparing to full CMB fits, we show that combining the shift parameter with the position of the first acoustic peak in the CMB power spectrum improves the accuracy of the test considerably.

The large-scale dynamics of a two-fluid system with a time dependent interaction is studied analytically and numerically. We show how a rapid transition can significantly suppress the large-scale curvature perturbation and present approximative formulae for estimating the effect. By comparing to numerical results, we study the applicability of the approximation and find good agreement with exact calculations.

Static spherically symmetric perfect fluid solutions are studied in metric $f(R)$ theories of gravity. We show that pressure and density do not uniquely determine $f(R)$ ie. given a matter distribution and an equation state, one cannot determine the functional form of $f(R)$. However, we also show that matching the outside Schwarzschild-de Sitter-metric to the metric inside the mass distribution leads to additional constraints that severely limit the allowed fluid configurations.

Spherically symmetric static empty space solutions are studied in f(R) theories of gravity. We reduce the set of modified Einstein's equations to a single equation and show how one can construct exact solutions in different f(R) models. In particular, we show that for a large class models, including e.g. the f(R)=R-\mu^4/R model, the Schwarzschild-de Sitter metric is an exact solution of the field equations. The significance of these solutions is discussed in light of solar system constraints on $f(R)$ theories of gravity.

We investigate f(R) theories of gravity within the Palatini approach and show how one can determine the expansion history, H(a), for an arbitrary choice of f(R). As an example, we consider cosmological constraints on such theories arising from the supernova type Ia, large scale structure formation and cosmic microwave background observations. We find that best fit to the data is a non-null leading order correction to the Einstein gravity, but the current data exhibits no significant preference over the concordance LCDM model. Our results show that the often considered 1/R models are not compatible with the data. The results demonstrate that the background expansion alone can act as a good discriminator between modified gravity models when multiple data sets are used.

Assuming only a homogeneous and isotropic universe and using both the 'Gold' Supernova Type Ia sample of Riess et al. and the results from the Supernova Legacy Survey, we calculate the Bayesian evidence of a range of different parameterizations of the deceleration parameter. We consider both spatially flat and curved models. Our results show that although there is strong evidence in the data for an accelerating universe, there is little evidence that the deceleration parameter varies with redshift.

A class of supergravity inspired quintessence models is studied by comparing to cosmological data. The set of considered models includes several previously studied quintessential potentials, as well as the $\Lambda$CDM model. We find that even though the commonly studied supergravity inspired quintessence models fit the data better than the $\Lambda$CDM model, they are a relatively poor fit when compared to the best fit model in the studied class. Our results suggest a low energy scale, less than ${\cal{M}}\sim 1$ TeV, for the effective supergravity potential.

We argue that string theory has all the ingredients to provide us with candidates for the cold dark matter and explain the current acceleration of our Universe. In any generic string compactification the dilaton plays an important role as it couples to the Standard Model and other heavy non-relativistic degrees of freedom such as the string winding modes and wrapped branes, we collectively call them stringy cold dark matter. These couplings are non-universal which results in an interesting dynamics for a rolling dilaton. Initially, its potential can track radiation and matter while beginning to dominate the dynamics recently, triggering a phase of acceleration. This scenario can be realized as long as the dilaton also couples strongly to some heavy modes. We furnish examples of such modes. We provide analytical and numerical results and compare them with the current supernovae result. This favors certain stringy candidates.

Recent studies by a number of independent collaborations, have correlated the CMB temperatures measured by the WMAP satellite with different galaxy surveys that trace the matter distribution with light from the whole range of the electromagnetic spectrum: radio, far-infrared, optical and X-ray surveys. The new data systematically finds positive correlations, indicating a rapid slow down in the growth of structure in the universe. Individual cross-correlation measurements are of low significance, but we show that combining data at different redshifts introduces important new constraints. Contrary to what happens at low redshifts, for a fixed $\Omm$, the higher the dark energy contend, $\Ol$, the lower the ISW cross-correlation amplitude. At 68% confidence level, the data finds new independent evidence of dark energy: $\Ol =0.42-1.22$ . It also confirms, to higher significance, the presence of a large dark matter component: $\Omm =0.18-0.34$, exceeding the density of baryonic matter, but far from the critical value. Combining these new constraints with the prior of a flat universe, or the prior of an accelerating universe provides strong new evidence for a dark cosmos. Combination with supernova data yields $\Ol = 0.71 \pm 0.13$, $\Omm = 0.29 \pm 0.04$. If we also assume a flat universe, we find $\Ol = 0.70 \pm 0.05$ and $w = -1.02 \pm 0.17$ for a constant dark energy equation of state.

Modified theories of gravity have recently been studied by several authors as possibly viable alternatives to the cosmological concordance model. Such theories attempt to explain the accelerating expansion of the universe by changing the theory of gravity, instead of introducing dark energy. In particular, a class of models based on higher order curvature invariants, so-called $f(R)$ gravity models, has drawn attention. In this letter we show that within this framework, the expansion history of the universe does not uniquely determine the form of the gravitational action and it can be radically different from the standard Einstein-Hilbert action. We demonstrate that for any barotropic fluid, there always exists a class of $f(R)$ models that will have exactly the same expansion history as that arising from the Einstein-Hilbert action. We explicitly show how one can extend the Einstein-Hilbert action by constructing a $f(R)$ theory that is equivalent on the classical level. Due to the classical equivalence between $f(R)$ theories and Einstein-Hilbert gravity with an extra scalar field, one can also hence construct equivalent scalar-tensor theories with standard expansion.

Motivated by the string landscape we examine scenarios for which inflation is a two-step process, with a comparatively short inflationary epoch near the string scale and a longer period at a much lower energy (like the TeV scale). We quantify the number of $e$-foldings of inflation which are required to yield successful inflation within this picture. The constraints are very sensitive to the equation of state during the epoch between the two inflationary periods, as the extra-horizon modes can come back inside the horizon and become reprocessed. We find that the number of $e$-foldings during the first inflationary epoch can be as small as 12, but only if the inter-inflationary period is dominated by a network of cosmic strings (such as might be produced if the initial inflationary period is due to the brane-antibrane mechanism). In this case a further 20 $e$-foldings of inflation would be required at lower energies to solve the late universe's flatness and horizon problems.

We study non-Gaussianity in two distinct models of preheating: instant and tachyonic. In instant preheating non-Gaussianity is sourced by the local terms generated through the coupled perturbations of the two scalar fields. We find that the non-Gaussianity parameter is given by $f_{NL}^{\phi}\sim 2g < O(1)$, where $g$ is a coupling constant, so that instant preheating is unlikely to be constrained by WMAP or Planck. In the case of tachyonic preheating non-Gaussianity arises solely from the instability of the tachyon matter and is found to be large. We find that for single field inflation the present WMAP data implies a bound $V_{0}^{1/4}/M_{P}\leq 10^{-4}$ on the scale of tachyonic instability. We argue that the tachyonic preheating limits are useful also for string-motivated inflationary models.

We consider a two-field model for inflation where the second order metric perturbations can be amplified by a parametric resonance during preheating. We demonstrate that there can arise a considerable enhancement of non-Gaussianity sourced by the local terms generated through the coupled perturbations. We argue that the non-Gaussianity parameter could be as large as f_{NL} ~ 50. Our results may provide a useful test of preheating in future CMB experiments.

We demonstrate that string motivated inflation ending via tachyonic instability leaves a detectable imprint on the cosmic microwave background (CMB) radiation by virtue of the excitation of non-Gaussian gravitational fluctuations. The present WMAP bound on non-Gaussianity is shown to constrain the string scale by $M_S/M_P\leq 10^{-4}$ for string coupling $g_{s}<0.1$, hence improving the existing bounds. If tachyon fluctuations during inflation are not negligible, we find the stringent constraint $g_s\sim 10^{-9}$ for $M_{S}/M_{P}<10^{-3}$. This case may soon be ruled out by the forthcoming CMB non-Gaussinianity bounds.

We study the one-loop effective action for gravity in a cosmological setup to determine possible cosmological effects of quantum corrections to Einstein theory. By considering the effect of the universal non-local terms in a toy model, we show that they can play an important role in the very early universe. We find that during inflation, the non-local terms are significant, leading to deviations from the standard inflationary expansion.

Spatially flat loitering universe models have recently been shown to arise in the context of brane world scenarios. Such models allow more time for structure formation to take place at high redshifts, easing, e.g., the tension between the observed and predicted evolution of the quasar population with redshift. While having the desirable effect of boosting the growth of structures, we show that in such models the position of the first peak in the power spectrum of the cosmic microwave background anisotropies severely constrains the amount of loitering at high redshifts.

We confront Cardassian models with recent observational data. These models can be viewed either as purely phenomenological modifications of the Friedmann equation, or as arising from cosmic fluids with non-standard properties. In the first case, we find that the models are consistent with the data for a wide range of parameters but no significant preference over the Lambda-CDM model is found. In the latter case we find that the fact that the sound speed is non-zero in these models makes them inconsistent with the galaxy power spectrum from the Sloan Digital Sky Survey.

An alternative to dark energy as an explanation for the present phase of accelerated expansion of the Universe is that the Friedmann equation is modified, e.g. by extra dimensional gravity, on large scales. We explore a natural parametrization of a general modified Friedmann equation, and find that the present supernova type Ia and cosmic microwave background data prefer a correction of the form 1/H over a cosmological constant. We also explore the constraints that can be expected in the future, and find that there are good prospects for distinguishing this model from the standard cosmological constant to very high significance if one combines supernova data with a precise measurement of the matter density.

The Integrated Sachs-Wolfe effect is studied in non-standard cosmologies. By considering flat universes with a non-fluctuating dark energy component, it is shown how the quadrupole power can be suppressed by atypical evolution of the scale factor. For example, a brief period of non-standard evolution at a high redshift can suppress the quadrupole significantly. The effect on the overall normalization of the CMB power spectrum is also discussed. Non-standard cosmologies can affect the overall normalization significantly and enhance the primordial fluctuations. The possibility of constraining such non-standard models with CMB and independent measures of $\sigma_8$, is considered.

The growth of large scale structure is studied in a universe containing both cold dark matter (CDM) and generalized Chaplygin gas (GCg). GCg is assumed to contribute only to the background evolution of the universe while the CDM component collapses and forms structures. We present some new analytical as well as numerical results for linear and non-linear growth in such model. The model passes the standard cosmological distance test without the need of a cosmological constant (LCDM). But we find that the scenario is severely constrained by current observations of large scale structure. Any small deviations of the GCg parameters away from the standard Lambda dominated cosmology (LCDM) produces substantial suppression for the growth of structures.

Large-scale modes in the temperature anisotropy power spectrum C_l measured by the Wilkinson Microwave Anisotropy Probe (WMAP), seem to have lower amplitudes (C_2, C_3 and C_4) than that expected in the so called concordance LCDM model. In particular, the quadrupole C_2 is reported to have a smaller value than allowed by cosmic variance. This has been interpreted as a possible indication of new physics. In this paper we re-analyse the WMAP data using the 2-point angular correlation and its higher-order moments. This method, which requires a full covariance analysis, is more direct and provides better sampling of the largest modes than the standard harmonic decomposition. We show that the WMAP data is in good agreement (~ 30% probability) with a LCDM model when the WMAP data is considered as a particular realization drawn from a set of realistic LCDM simulations with the corresponding covariance. This is also true for the higher-order moments, shown here up to 6th order, which are consistent with the Gaussian hypothesis. The sky mask plays a major role in assessing the significance of these agreements. We recover the best fit model for the low-order multipoles based on the 2-point correlation with different assumptions for the covariance. Assuming that the observations are a fair sample of the true model, we find C_2 = 123 +/- 233, C_3= 217 +/- 241 and C_4 = 212 +/- 162 (in mu K^2). The errors increase by about a factor of 5 if we assume the \lcdm model. If we exclude the Galactic plane |b|<30 from our analysis, we recover very similar values within the errors (ie C_2=172, C_3= 89, C_4=129). This indicates that the Galactic plane is not responsible for the lack of large-scale power in the WMAP data.

We study the growth of large scale structure in two recently proposed non-standard cosmological models: the brane induced gravity model of Dvali, Gabadadze and Porrati (DGP) and the Cardassian models of Freese and Lewis. A general formalism for calculating the growth of fluctuations in models with a non-standard Friedman equation and a normal continuity equation of energy density is discussed. Both linear and non-linear growth are studied, together with their observational signatures on higher order statistics and abundance of collapsed objects. In general, models which show similar cosmic acceleration at z ~ 1, can produce quite different normalization for large scale density fluctuations, ie sigma_8, cluster abundance or higher order statistics, such as the normalized skewness S_3, which is independent of the linear normalization. For example, for a flat universe with Omega_M ~ 0.22, DGP and standard Cardassian cosmologies predict about 2 and 3 times more clusters respectively than the standard Lambda model at z=1.5. When normalized to CMB fluctuations the sigma_8 amplitude turns out to be lower by a few tens of a percent. We also find that, for a limited red-shift range, the linear growth rate can be faster in some models (eg modified polytropic Cardassian with q>1) than in the Einstein-deSitter universe. The value of the skewness S_3 is found to have up to ~10 percent variations (up or down) from model to model.

We consider different types of Q-balls as self-interacting dark matter. For the Q-balls to act as the dark matter of the universe they should not evaporate, which requires them to carry very large charges; depending on the type, the minimum charge could be as high as Q \sim 10^{33} or the Q-ball coupling to ordinary matter as small as \sim 10^{-35}. The cross-section-to-mass ratio needed for self-interacting dark matter implies a mass scale of m \sim O(1) MeV for the quanta that the Q-balls consist of, which is very difficult to achieve in the MSSM.

We present numerical simulations of fragmentation of the Affleck-Dine condensate in two spatial dimensions. We argue analytically that the final state should consist of both Q-balls and anti-Q-balls in a state of maximum entropy, with most of the balls small and relativistic. Such a behaviour is found in simulations on a 100x100 lattice with cosmologically realistic parameter values. During fragmentation process, we observe filament-like texture in the spatial distribution of charge. The total charge in Q-balls is found to be almost equal to the charge in anti-Q-balls and typically orders of magnitude larger than charge asymmetry. Analytical considerations indicate that, apart from geometrical factors, the results of the simulated two dimensional case should apply also to the fully realistic three dimensional case.