Motivated by a scenario of magnetogenesis in which a homogeneous magnetic field is generated during inflation, we study the magnetohydrodynamic evolution of the primordial plasma motions for two kinds of initial conditions -- (i) a spatially homogeneous imposed field with an unlimited correlation length, and (ii) a zero flux scale-invariant statistically homogeneous magnetic field. In both cases, we apply, for a short initial time interval, monochromatic forcing at a certain wave number so that the correlation length is finite, but much smaller than the typical length scale of turbulence. In particular, we investigate the decay of non-helical and helical hydromagnetic turbulence. We show that, in the presence of a homogeneous magnetic field, the decay of helical and nonhelical small-scale fields can occur rapidly. This is a special property of a system with a perfectly homogeneous magnetic field, which is sometimes considered as a local approximation to a slowly varying background field. It can never change and acts as an imposed magnetic field. This is in a sharp contrast to the case of a statistically homogeneous magnetic field, where we recover familiar decay properties: a much slower decay of magnetic energy and a faster growth of the correlation length, especially in the case with magnetic helicity. The result suggests that a homogeneous magnetic field, if generated during inflation, should persist under the influence of small-scale fields and could be the origin of the large-scale magnetic field in the universe.

This paper focuses on two aspects of the statistics of cosmological observables that are important for the next stages of precision cosmology. First, we note that the theory of reduced angular $N$-point spectra has only been developed in detail up to the trispectrum case and in a fashion that makes it difficult to go beyond. To fill this gap, here we present a constructive approach that provides a systematic description of reduced angular $N$-point spectra and their covariance matrices, for arbitrary $N$. Second, we focus on the common practice in the literature on cosmological observables, which consists in simply discarding a part of the expression, namely, the terms containing fields evaluated at the observer position. We point out that this is not justified beyond linear order in perturbation theory, as these terms contribute to all the multipoles of the corresponding spectra and with a magnitude that is of the same order as the rest of the non-linear corrections. We consider the possibility that the reason for neglecting these terms is a conceptual discomfort when using ensemble averages, which originates in an apparent tension between the ergodic hypothesis and the privileged position of the observer on the light-cone. We clarify this subtle issue by performing a careful derivation of the relation between the theoretical statistical predictions and the observational estimators for all $N$. We conclude that there is no inconsistency whatsoever in ensemble-averaging fields at and near the observer position, thus clearing the way for consistent and robust high-precision calculations.

In this paper we derive and test the flat sky approximation for galaxy number counts. We show that, while for the lensing term it reduces to the Limber approximation, for the standard density and redshift space distortion it is different and very accurate already at low $\ell$ while the corresponding Limber approximation completely fails. At equal redshift the accuracy of the standard terms is around 0.2% at low redshifts and 0.5% for redshift $z=5$, even to low $\ell$. At unequal redshifts the precision is less impressive and can only be trusted for very small redshift differences, $\Delta z<\Delta z_0 \simeq 3.6\times10^{-4}(1+z)^{2.14}$, but the lensing terms dominate for $\Delta z>\Delta z_1 \simeq 0.33(r(z)H(z))/(z+1)$. The Limber approximation achieves an accuracy of 0.5% above $\ell\simeq 40$ for the pure lensing term and above $\ell\simeq 80$ for the lensing-density cross-correlation. Besides being very accurate, the flat sky approximation is also very fast and can therefore be useful for data analysis and forecasts with MCMC methods.

In this paper we study gravitational lensing of gravitational wave events. The main point of the present work is to introduce a semi-analytic approach so that each ingredient can be varied and tested individually. Our analytic model for the source population is motivated by a numerical study and we compare semi-analytical and numerical results. We determine the expected magnification for events seen at a given observed luminosity distance. We find that while the probability of significant magnification of the observed LIGO-Virgo events is very small, the probability distribution of the magnification has a significant tail to high $\mu$ such that e.g. the variance of the magnification is very large and even diverges in the geometric optics approximation. For the 10 binary black hole mergers observed by LIGO-Virgo in the O1+O2 observation campaigns, the probability that one of them has been magnified with magnification of 5 or bigger is $\mathcal{P}_{\text{obs}}(>5)\sim 0.01$ while the probability of magnification by 50 or bigger is $\mathcal{P}_{\text{obs}}(>50)\sim 0.005$.

It is well known that helical magnetic fields undergo a so-called inverse cascade by which their correlation length grows due to the conservation of magnetic helicity in classical ideal magnetohydrodynamics (MHD). At high energies above approximately $10$ MeV, however, classical MHD is necessarily extended to chiral MHD and then the conserved quantity is $\langle\mathcal{H}\rangle + 2 \langle\mu_5\rangle / \lambda$ with $\langle\mathcal{H}\rangle$ being the mean magnetic helicity and $\langle\mu_5\rangle$ being the mean chiral chemical potential of charged fermions. Here, $\lambda$ is a (phenomenological) chiral feedback parameter. In this paper, we study the evolution of the chiral MHD system with the initial condition of nonzero $\langle\mathcal{H}\rangle$ and vanishing $\mu_5$. We present analytic derivations for the time evolution of $\langle\mathcal{H}\rangle$ and $\langle\mu_5\rangle$ that we compare to a series of laminar and turbulent three-dimensional direct numerical simulations. We find that the late-time evolution of $\langle\mathcal{H}\rangle$ depends on the magnetic and kinetic Reynolds numbers ${\rm Re}_{_\mathrm{M}}$ and ${\rm Re}_{_\mathrm{K}}$. For a high ${\rm Re}_{_\mathrm{M}}$ and ${\rm Re}_{_\mathrm{K}}$ where turbulence occurs, $\langle\mathcal{H}\rangle$ eventually evolves in the same way as in classical ideal MHD where the inverse correlation length of the helical magnetic field scales with time $t$ as $k_\mathrm{p} \propto t^{-2/3}$. For a low Reynolds numbers where the velocity field is negligible, the scaling is changed to $k_\mathrm{p} \propto t^{-1/2}\mathrm{ln}\left(t/t_\mathrm{log}\right)$. After being rapidly generated, $\langle\mu_5\rangle$ always decays together with $k_\mathrm{p}$, i.e. $\langle\mu_5\rangle \approx k_\mathrm{p}$, with a time evolution that depends on whether the system is in the limit of low or high Reynolds numbers.

The subject of this paper are shape and size correlations of galaxies due to weak gravitational lensing and due to direct tidal interaction of elliptical galaxies with gravitational fields sourced by the cosmic large-scale structure. Setting up a linear intrinsic alignment model for elliptical galaxies which parameterises the reaction of the galaxy to an external tidal shear field through the velocity dispersion, we predict intrinsic correlations and cross-correlations with weak lensing for both shapes and sizes, juxtaposing both types of spectra with lensing. We quantify the observability of the intrinsic shape and size correlations and estimate with the Fisher-formalism how well the alignment parameter can be determined from the Euclid weak lensing survey. Specifically, we find a contamination of the weak lensing convergence spectra with an intrinsic size correlation amounting to up to 10% over a wide multipole range $\ell=100\ldots 300$, with a corresponding cross-correlation exhibiting a sign change, similar to the cross-correlation between weak lensing shear and intrinsic shapes. A determination of the alignment parameter yields a precision of a few percent forecasted for Euclid, and we show that all shape and many size correlations should be measurable with Euclid.

We study the importance of gravitational lensing in the modelling of the number counts of galaxies. We confirm previous results for photometric surveys, showing that lensing cannot be neglected in a survey like LSST since it would infer a significant shift of cosmological parameters. For a spectroscopic survey like SKA2, we find that neglecting lensing in the monopole, quadrupole and hexadecapole of the correlation function also induces an important shift of parameters. For $\Lambda$CDM parameters, the shift is moderate, of the order of $0.6\sigma$ or less. However, for a model-independent analysis, that measures the growth rate of structure in each redshift bin, neglecting lensing introduces a shift of up to $2.3\sigma$ at high redshift. Since the growth rate is directly used to test the theory of gravity, such a strong shift would wrongly be interpreted as the breakdown of General Relativity. This shows the importance of including lensing in the analysis of future surveys. On the other hand, for a survey like DESI, we find that lensing is not important, mainly due to the value of the magnification bias parameter of DESI, $s(z)$, which strongly reduces the lensing contribution at high redshift. We also propose a way of improving the analysis of spectroscopic surveys, by including the cross-correlations between different redshift bins (which is neglected in spectroscopic surveys) from the spectroscopic survey or from a different photometric sample. We show that including the cross-correlations in the SKA2 analysis does not improve the constraints. On the other hand replacing the cross-correlations from SKA2 by cross-correlations measured with LSST improves the constraints by 10 to 20 %. Interestingly, for $\Lambda$CDM parameters, we find that LSST and SKA2 are highly complementary, since they are affected differently by degeneracies between parameters.

We investigate for the first time weak lensing in a high-resolution relativistic N-body simulation, which includes all relevant effects from general relativity. We integrate the photon geodesics backwards from the observer to the emitters. Our analysis is fully relativistic and non-perturbative for the scalar part of the gravitational potential and first-order in the vector part, frame dragging. We solve the Sachs optical equations and study in detail the weak-lensing convergence, ellipticity and rotation. We present the angular power spectra and one-point probability distribution functions for the weak-lensing variables, which we find are broadly in agreement with comparable Newtonian simulations. Our geometric approach, however, is more robust and flexible, and can therefore be applied consistently to non-standard cosmologies and modified theories of gravity.

In this short paper we determine the effects of structure on the cosmological consistency relation which is valid in a perfect Friedmann Universe. We show that within $\Lambda$CDM the consistency relation is violated by about 1.5% for redshifts $z\simeq 2$ and this violation raises up to 2% at $z\simeq 5$ and 3% at $z\simeq 10$ after which it settles at about 2.7% for $z>10$. This effect of cosmic structure on the distance redshift relation is also very sensitive to the determination of the dark energy equation of state via cosmic distances. It actually leads to an apparent unphysical behavior of $w(z)$ which even diverges at $z\sim 2.5$ and which we also discuss here.

The background photon temperature $\bar T$ is one of the fundamental cosmological parameters. Despite its significance, $\bar T$ has never been allowed to vary in the data analysis, owing to the precise measurement of the comic microwave background (CMB) temperature by COBE FIRAS. However, even in future CMB experiments, $\bar T$ will remain unknown due to the unknown monopole contribution $\Theta_0$ at our position to the observed (angle-averaged) temperature $\langle T\rangle^{\rm obs}$. By fixing $\bar T\equiv\langle T\rangle^{\rm obs}$, the standard analysis underestimates the error bars on cosmological parameters, and the best-fit parameters obtained in the analysis are biased in proportion to the unknown amplitude of $\Theta_0$. Using the Fisher formalism, we find that these systematic errors are smaller than the error bars from the $Planck$ satellite. However, with $\bar T\equiv\langle T\rangle^{\rm obs}$, these systematic errors will always be present and irreducible, and future cosmological surveys might misinterpret the measurements.

Future galaxy clustering surveys will probe small scales where non-linearities become important. Since the number of modes accessible on intermediate to small scales is very high, having a precise model at these scales is important especially in the context of discriminating alternative cosmological models from the standard one. In the mildly non-linear regime, such models typically differ from each other, and galaxy clustering data will become very precise on these scales in the near future. As the observable quantity is the angular power spectrum in redshift space, it is important to study the effects of non-linear density and redshift space distortion (RSD) in the angular power spectrum. We compute non-linear contributions to the angular power spectrum using a flat-sky approximation that we introduce in this work, and compare the results of different perturbative approaches with $N$-body simulations. We find that the TNS perturbative approach is significantly closer to the $N$-body result than Eulerian or Lagrangian 1-loop approximations, effective field theory of large scale structure or a halofit-inspired model. However, none of these prescriptions is accurate enough to model the angular power spectrum well into the non-linear regime. In addition, for narrow redshift bins, $\Delta z \lesssim 0.01$, the angular power spectrum acquires non-linear contributions on all scales, right down to $\ell=2$, and is hence not a reliable tool at this time. To overcome this problem, we need to model non-linear RSD terms, for example as TNS does, but for a matter power spectrum that remains reasonably accurate well into the deeply non-linear regime, such as halofit.

An important part of cosmological model fitting relies on correlating distance indicators of objects (for example type Ia supernovae) with their redshift, often illustrated on a Hubble diagram. Comparing the observed correlation with a homogeneous model is one of the key pieces of evidence for dark energy. The presence of cosmic structures introduces a bias and scatter, mainly due to gravitational lensing and peculiar velocities, but also due to smaller non-linear relativistic contributions which are more difficult to account for. For the first time we perform ray tracing onto halos in a relativistic N-body simulation. Our simulation is the largest that takes into account all leading-order corrections from general relativity in the evolution of structure, and we present a novel methodology for working out the non-linear projection of that structure onto the observer's past light cone. We show that the mean of the bias in the Hubble diagram is indeed as small as expected from perturbation theory. However, the distribution of sources is significantly skewed with a very long tail of highly magnified objects and we illustrate that the bias of cosmological parameters strongly depends on the function of distance which we consider.

We investigate the weak lensing corrections to the CMB polarization anisotropies. We concentrate on the effect of rotation and show that the rotation of polarisation is a true physical effect which has to be taken into account at second order in perturbation theory. We clarify inconsistencies on the treatment of this rotation in the recent literature. We also show that at first order in perturbation theory there is no rotation of polarisation also for vector and tensor modes.

Even though we know that physical observations are frame independent, the frame dependence of cosmological perturbations is relatively subtle and has led to confusion in the past. In this paper we show that while the (unobservable) matter power spectrum is frame dependent, the observable number counts are not. We shall also determine how the frame dependence of the power spectrum depends on scale.

Recently Moradinezhad Dizgah & Durrer have shown that the $E_g$ statistics, useful to test theories of modified gravity, is plagued by additional scale and bias dependent lensing contributions. In this work we develop and illustrate a method to remove these lensing terms by using in addition to the galaxy clustering data also shear data and the correlations of shear and galaxy clustering. We introduce a truly observable statistics termed $\tilde E_g$ which conserves the properties of scale and bias independence on linear scales. The method discussed here is best adapted to photometric surveys. It is found that the corrections to the original $E_g$ statistics are small for the present DES data, but for future surveys of the quality of Euclid they are very substantial.

We develop a theoretical framework to describe the cosmological observables on the past light cone such as the luminosity distance, weak lensing, galaxy clustering, and the cosmic microwave background anisotropies. We consider that all the cosmological observables include not only the background quantity, but also the perturbation quantity, and they are subject to cosmic variance, which sets the fundamental limits on the cosmological information that can be derived from such observables, even in an idealized survey with an infinite number of observations. To quantify the maximum cosmological information content, we apply the Fisher information matrix formalism and spherical harmonic analysis to cosmological observations, in which the angular and the radial positions of the observables on the light cone carry different information. We discuss the maximum cosmological information that can be derived from five different observables: (1) type Ia supernovae, (2) cosmic microwave background anisotropies, (3) weak gravitational lensing, (4) local baryon density, and (5) galaxy clustering. We compare our results with the cosmic variance obtained in the standard approaches, which treat the light cone volume as a cubic box of simultaneity. We discuss implications of our formalism and ways to overcome the fundamental limit.

We compute the redshift-dependent angular bispectrum of galaxy number counts at tree-level, including nonlinear clustering bias and estimating numerically for the first time the effect of redshift space distortions (RSD). We show that for narrow redshift bins the amplitude of nonlinear RSD is comparable with the matter density perturbations. While our numerical results only include terms relevant on sub-horizon scales, the formalism can readily be extended to the full tree-level bispectrum. Our approach does not rely on the flat-sky approximation and it can be easily generalized to different sources by including the appropriate bias expansion. We test the accuracy of Limber approximation for different z-bins. We highlight the subtle but relevant differences in the angular bispectrum of galaxy number counts with respect to CMB, due to the different scale dependence of perturbations. Our formalism can also be directly applied to the angular HI intensity mapping bispectrum.

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic microwave background (CMB) have focused on temperature and polarization anisotropies. CMB spectral distortions - tiny departures of the CMB energy spectrum from that of a perfect blackbody - provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, which highlight the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe - its thermal history - thereby providing additional insight into processes within the cosmological standard model (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space could be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets, some of which are already observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. The absence of these signals would pose a huge theoretical challenge, immediately pointing to new physics.

We discuss a model which can generate scale-invariant helical magnetic fields on large scales ($\lesssim 1$Mpc) in the primordial universe. It is also shown that the electric conductivity becomes significant and terminates magnetogenesis even before reheating is completed. By solving the electromagnetic dynamics taking conductivity into account, we find that magnetic fields with amplitude $B\simeq 10^{-15}{\rm G}$ at present can be generated without encountering a backreaction or strong coupling problem.

During single clock inflation, hidden fields (i.e. fields coupled to the inflaton only gravitationally) in their adiabatic vacua can ordinarily only affect observables through virtual effects. After renormalizing background quantities (fixed by observations at some pivot scale), all that remains are logarithmic runnings in correlation functions that are both Planck and slow roll suppressed. In this paper we show how a large number of hidden fields can partially compensate this suppression and generate a potentially observable running in the tensor two point function, consistently inferable courtesy of a large $N$ resummation. We detour to address certain subtleties regarding loop corrections during inflation, extending the analysis of [1]. Our main result is that one can extract bounds on the hidden field content of the universe from bounds on violations of the consistency relation between the tensor spectral index and the tensor to scalar ratio, were primordial tensors ever detected. Such bounds are more competitive than the naive bound inferred from requiring inflation to occur below the strong coupling scale of gravity if deviations from the consistency relation can be bounded to within the sub-percent level. We discuss how one can meaningfully constrain the parameter space of various phenomenological scenarios and constructions that address naturalness with a large number of species (such as `N-naturalness') with CMB observations up to cosmic variance limits, and possibly future 21cm and gravitational wave observations.

A persistent theme in the study of dark energy is the question of whether it really exists or not. It is often claimed hat we are mis-calculating the cosmological model by neglecting the effects associated with averaging over large-scale structures. In the Newtonian approximation this is clear: there is no effect. Within the full relativistic picture this remains an important open question however, owing to the complex mathematics involved. We study this issue using particle numerical simulations which account for all relevant relativistic effects without any problems from shell crossing. In this context we show for the first time that the backreaction from structure can differ by many orders of magnitude depending upon the slicing of spacetime one chooses to average over. In the worst case, where smoothing is carried out in synchronous spatial surfaces, the corrections can reach ten percent and more. However, when smoothing on the constant time hypersurface of the Newtonian gauge, backreaction contributions remain 3-5 orders of magnitude smaller.

The geometric optics approximation traditionally used to study the propagation of gravitational waves on a curved background, breaks down in the vicinity of compact and extended astrophysical objects, where wave-like effects like diffusion and generation of polarization occur. We provide a framework to study the generation of polarization of a stochastic background of gravitational waves propagating in an inhomogeneous universe. The framework is general and can be applied to both cosmological and astrophysical gravitational wave backgrounds in any frequency range. We derive an order of magnitude estimate of the amount of polarization generated for cosmological and astrophysical backgrounds, in the frequency range covered by present and planned gravitational wave experiments. For an astrophysical background in the PTA and LISA band, the amount of polarization generated is suppressed by a factor $10^{-4}$ $(10^{-5})$ with respect to anisotropies. For a cosmological background we get an additional $10^{-2}$ suppression. We speculate on using our approach to map the distribution of (unresolvable) structures in the Universe.

Clustering of a perfect fluid does not lead to the generation of vorticity. It is the collisionless nature of dark matter, inducing velocity dispersion and shell crossing, which is at the origin of cosmological vorticity generation. In this paper we investigate the generation of vorticity during the formation of cosmological large scale structure using the public relativistic N-body code $\textit{gevolution}$. We test several methods to compute the vorticity power spectrum and we study its convergence with respect to the mass and grid resolution of our simulations. We determine the power spectrum, the spectral index on large-scales, the amplitude of the peak position and their time evolution. We also compare the vorticity extracted from our simulations with the vector perturbations of the metric. Our results are accompanied by resolution studies and compared with previous studies in the literature.

We investigate the statistical properties of isotropic, stochastic, Gaussian distributed, helical magnetic fields characterized by different shapes of the energy spectra at large length scales and study the associated realizability condition. We discuss smoothed magnetic fields that are commonly used when the primordial magnetic field is constrained by observational data. We are particularly interested in scale-invariant magnetic fields that can be generated during the inflationary stage by quantum fluctuations. We determine the correlation length of such magnetic fields and relate it to the infrared cutoff of perturbations produced during inflation. We show that this scale determines the observational signatures of the inflationary magnetic fields on the cosmic microwave background. At smaller scales, the scale-invariant spectrum changes with time. It becomes a steeper weak-turbulence spectrum at progressively larger scales. We show numerically that the critical length scale where this happens is the turbulent-diffusive scale, which increases with the square root of time.

In this paper, we study lensing of 21cm intensity mapping (IM). Like in the cosmic microwave background (CMB), there is no first order lensing in intensity mapping. The first effects in the power spectrum are therefore of second and third order. Despite this, lensing of the CMB power spectrum is an important effect that needs to be taken into account, which motivates the study of the impact of lensing on the IM power spectrum. We derive a general formula up to third order in perturbation theory including all the terms with two derivatives of the gravitational potential, i.e. the dominant terms on sub-Hubble scales. We then show that in intensity mapping there is a new lensing term which is not present in the CMB. We obtain that the signal-to-noise of 21 cm lensing for futuristic surveys like SKA2 is about 10. We find that surveys probing only large scales, lmax < 700, can safely neglect the lensing of the intensity mapping power spectrum, but that otherwise this effect should be included.

We investigate the weak lensing corrections to the CMB temperature and polarization anisotropies. We consider all the effects beyond the leading order: post-Born corrections, LSS corrections and, for the polarization anisotropies, the correction due to the rotation of the polarization direction between the emission at the source and the detection at the observer. We show that the full next-to-leading order correction to the B-mode polarization is not negligible on small scales and is dominated by the contribution from the rotation, this is a new effect not taken in account in previous works. Considering vanishing primordial gravitational waves, the B-mode correction due to rotation is comparable to cosmic variance for $\ell \gtrsim 3500$, in contrast to all other spectra where the corrections are always below that threshold for a single multipole. Moreover, the sum of all the effects is larger than cosmic variance at high multipoles, showing that higher-order lensing corrections to B-mode polarization are in principle detectable.

We study the impact on the galaxy correlation function of the presence of a vector component in the tracers' peculiar velocities, in the case in which statistical isotropy is violated. We present a general framework - based on the bipolar spherical harmonics expansion - to study this effect in a model independent way, without any hypothesis on the origin or the properties of these vector modes. We construct six new observables, that can be directly measured in galaxy catalogs in addition to the standard monopole, quadrupole and hexadecapole, and we show that they completely describe any deviations from isotropy. We then perform a Fisher analysis in order to quantify the constraining power of future galaxy surveys. As an example, we show that the SKA2 would be able to detect anisotropic rotational velocities with amplitudes as low as 1% of that of the vorticity generated during shell-crossing in standard dark matter scenarios.

We present a public version of the code COFFE (COrrelation Function Full-sky Estimator) available at https://github.com/JCGoran/coffe. The code computes the galaxy two-point correlation function and its multipoles in linear perturbation theory, including all relativistic and wide angle corrections. COFFE also calculates the covariance matrix for two physically relevant estimators of the correlation function multipoles. We illustrate the usefulness of our code by a simple but relevant example: a forecast of the detectability of the lensing signal in the multipoles of the two-point function. In particular, we show that lensing should be detectable in the multipoles of the two-point function, with a signal-to-noise larger than 10, in future surveys like Euclid or the SKA.

This paper revisits the controversy concerning whether gravitational lensing effects make a significant difference to estimation of distance to the Cosmic Microwave Background last scattering surface in cosmology. A recent paper by Kaiser and Peacock supports a previous paper by Weinberg stating that such affects average to zero because of energy conservation. In this note, problems are pointed out in the Kaiser and Peacock analysis related to their choice of endpoint of integration, and to the `wrinkly surface' argument.

We investigate the corrections which relativistic light-cone computations induce on the correlation of the tangential shear with galaxy number counts, also known as galaxy-galaxy lensing. The standard-approach to galaxy-galaxy lensing treats the number density of sources in a foreground bin as observable, whereas it is in reality unobservable due to the presence of relativistic corrections. We find that already in the redshift range covered by the DES first year data, these currently neglected relativistic terms lead to a systematic correction of up to 50% in the density-shear correlation function for the highest redshift bins. This correction is dominated by the the fact that a redshift bin of number counts does not only lens sources in a background bin, but is itself again lensed by all masses between the observer and the counted source population. Relativistic corrections are currently ignored in the standard galaxy-galaxy analyses, and the additional lensing of a counted source populations is only included in the error budget (via the covariance matrix). At increasingly higher redshifts and larger scales, these relativistic and lensing corrections become however increasingly more important, and we here argue that it is then more efficient, and also cleaner, to account for these corrections in the density-shear correlations.

We derive an exact expression for the correlation function in redshift shells including all the relativistic contributions. This expression, which does not rely on the distant-observer or flat-sky approximation, is valid at all scales and includes both local relativistic corrections and integrated contributions, like gravitational lensing. We present two methods to calculate this correlation function, one which makes use of the angular power spectrum C_ell(z1,z2) and a second method which evades the costly calculations of the angular power spectra. The correlation function is then used to define the power spectrum as its Fourier transform. In this work theoretical aspects of this procedure are presented, together with quantitative examples. In particular, we show that gravitational lensing modifies the multipoles of the correlation function and of the power spectrum by a few percent at redshift z=1 and by up to 30% and more at z=2. We also point out that large-scale relativistic effects and wide-angle corrections generate contributions of the same order of magnitude and have consequently to be treated in conjunction. These corrections are particularly important at small redshift, z=0.1, where they can reach 10%. This means in particular that a flat-sky treatment of relativistic effects, using for example the power spectrum, is not consistent.

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

We consider a scale-invariant helical magnetic field generated during inflation. We show that, if the mean magnetic helicity density of such a field is measured, it can be used to determine a lower bound on the duration of inflation. Upper bounds can be used to derive constraints on the minimal duration of inflation if one assumes that the magnetic field generated during inflation is helical. Using three-dimensional simulations, we show that an initially scale-invariant field develops, which is similar both with and without magnetic helicity. In the fully helical case, however, the magnetic field appears to have a more pronounced folded structure.

Some of the dark matter in the Universe is made up of massive neutrinos. Their impact on the formation of large scale structure can be used to determine their absolute mass scale from cosmology, but to this end accurate numerical simulations have to be developed. Due to their relativistic nature, neutrinos pose additional challenges when one tries to include them in N-body simulations that are traditionally based on Newtonian physics. Here we present the first numerical study of massive neutrinos that uses a fully relativistic approach. Our N-body code, gevolution, is based on a weak-field formulation of general relativity that naturally provides a self-consistent framework for relativistic particle species. This allows us to model neutrinos from first principles, without invoking any ad-hoc recipes. Our simulation suite comprises some of the largest neutrino simulations performed to date. We study the effect of massive neutrinos on the nonlinear power spectra and the halo mass function, focusing on the interesting mass range between 0.06 eV and 0.3 eV and including a case for an inverted mass hierarchy.

Theoretical descriptions of observable quantities in cosmological perturbation theory should be independent of coordinate systems. This statement is often referred to as gauge-invariance of observable quantities, and the sanity of their theoretical description is verified by checking its gauge-invariance. We argue that cosmological observables are invariant scalars under diffeomorphisms and as a consequence their theoretical description is gauge-invariant, only at linear order in perturbations. Beyond linear order, they are usually not gauge-invariant, and we provide the general law for the gauge-transformation that the perturbation part of an observable does obey. We apply this finding to derive the second-order expression for the observational light-cone average in cosmology and demonstrate that our expression is indeed invariant under diffeomorphisms.

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.

In this Letter we study the impact on cosmological parameter estimation, from present and future surveys, due to lensing corrections on CMB temperature and polarization anisotropies beyond leading order. In particular, we show how post-Born corrections, LSS effects and the correction due to the change in the polarization direction between the emission at the source and the detection at the observer, are non-negligible in the determination of the polarization spectra. They have to be taken into account for an accurate estimation of cosmological parameters sensitive to or even based on these spectra. We study in detail the impact of higher order lensing on the determination of the tensor-to-scalar ratio $r$ and on the estimation of the effective number of relativistic species $N_\text{eff}$. We find that neglecting higher order lensing terms can lead to misinterpreting these corrections as a primordial tensor-to-scalar ratio of about $\mathcal{O}(10^{-3})$. Furthermore, it leads to a shift of the parameter $N_\text{eff}$ by nearly two sigma considering the level of accuracy aimed by future S4 surveys.

We compute the generation of vorticity from velocity dispersion in the dark matter fluid. For dark matter at zero temperature Helmholtz's theorem dictates that no vorticity is generated and we therefore allow the dark matter fluid to have a non-vanishing velocity dispersion. This implies a modification to the usual hydrodynamical system (continuity and Euler equations): we have to consider the Boltzmann hierarchy up to the second moment. This means that the Euler equation is modified with a source term that describes the effect of non-zero velocity dispersion. We write an equation for the Eulerian vorticity in Lagrangian coordinates and show that it has a growing mode already at second order in perturbation theory. We compute the power spectrum of the vorticity and the rotational velocity at second order in perturbation theory.

We review the number counts to second order concentrating on the terms which dominate on sub horizon scales. We re-derive the result for these terms and compare it with the different versions found in the literature. We generalize our derivation to higher order terms, especially the third order number counts which are needed to compute the 1-loop contribution to the power spectrum.

High-precision constraints on primordial non-Gaussianity (PNG) will significantly improve our understanding of the physics of the early universe. Among all the subtleties in using large scale structure observables to constrain PNG, accounting for relativistic corrections to the clustering statistics is particularly important for the upcoming galaxy surveys covering progressively larger fraction of the sky. We focus on relativistic projection effects due to the fact that we observe the galaxies through the light that reaches the telescope on perturbed geodesics. These projection effects can give rise to an effective $f_{\rm NL}$ that can be misinterpreted as the primordial non-Gaussianity signal and hence is a systematic to be carefully computed and accounted for in modelling of the bispectrum. We develop the technique to properly account for relativistic effects in terms of purely observable quantities, namely angles and redshifts. We give some examples by applying this approach to a subset of the contributions to the tree-level bispectrum of the observed galaxy number counts calculated within perturbation theory and estimate the corresponding non-Gaussianity parameter, $f_{\rm NL}$, for the local, equilateral and orthogonal shapes. For the local shape, we also compute the local non-Gaussianity resulting from terms obtained using the consistency relation for observed number counts. Our goal here is not to give a precise estimate of $f_{\rm NL}$ for each shape but rather we aim to provide a scheme to compute the non-Gaussian contamination due to relativistic projection effects. For the terms considered in this work, we obtain contamination of $f_{\rm NL}^{\rm loc} \sim {\mathcal O}(1)$.

We investigate the Alcock Paczynski (AP) test applied to the Baryon Acoustic Oscillation (BAO) feature in the galaxy correlation function. By using a general formalism that includes relativistic effects, we quantify the importance of the linear redshift space distortions and gravitational lensing corrections to the galaxy number density fluctuation. We show that redshift space distortions significantly affect the shape of the correlation function, both in radial and transverse directions, causing different values of galaxy bias to induce offsets up to 1% in the AP test. On the other hand, we find that the lensing correction around the BAO scale modifies the amplitude but not the shape of the correlation function and therefore does not introduce any systematic effect. Furthermore, we investigate in details how the AP test is sensitive to redshift binning: a window function in transverse direction suppresses correlations and shifts the peak position toward smaller angular scales. We determine the correction that should be applied in order to account for this effect, when performing the test with data from three future planned galaxy redshift surveys: Euclid, the Dark Energy Spectroscopic Instrument (DESI) and the Square Kilometer Array (SKA).

We study the impact of the often neglected lensing contribution to galaxy number counts on the $E_g$ statistics which is used to constrain deviations from GR. This contribution affects both the galaxy-galaxy and the convergence-galaxy spectra, while it is larger for the latter. At higher redshifts probed by upcoming surveys, neglecting this term induces a substantial error in the spectra and therefore on the $E_g$ statistics which is constructed from the combination of the two. Moreover, including it, renders the $E_g$ statistics scale and bias-dependent and hence puts into question its very objective.

We investigate the weak lensing corrections to the cosmic microwave background temperature anisotropies considering effects beyond the Born approximation. To this aim, we use the small deflection angle approximation, to connect the lensed and unlensed power spectra, via expressions for the deflection angles up to third order in the gravitational potential. While the small deflection angle approximation has the drawback to be reliable only for multipoles $\ell\lesssim 2500$, it allows us to consistently take into account the non-Gaussian nature of cosmological perturbation theory beyond the linear level. The contribution to the lensed temperature power spectrum coming from the non-Gaussian nature of the deflection angle at higher order is a new effect which has not been taken into account in the literature so far. It turns out to be the leading contribution among the post-Born lensing corrections. On the other hand, the effect is smaller than corrections coming from non-linearities in the matter power spectrum, and its imprint on CMB lensing is too small to be seen in present experiments.

We demonstrate the importance of including the lensing contribution in galaxy clustering analyses with large galaxy redshift surveys. It is well known that radial cross-correlations between different redshift bins of galaxy surveys are dominated by lensing. But we show here that also neglecting lensing in the autocorrelations within one bin severely biases cosmological parameter estimation with redshift surveys. It leads to significant shifts for several cosmological parameters, most notably the scalar spectral index and the neutrino mass scale. Especially the latter parameter is one of the main targets of future galaxy surveys.

We present a new N-body code, gevolution, for the evolution of large scale structure in the Universe. Our code is based on a weak field expansion of General Relativity and calculates all six metric degrees of freedom in Poisson gauge. N-body particles are evolved by solving the geodesic equation which we write in terms of a canonical momentum such that it remains valid also for relativistic particles. We validate the code by considering the Schwarzschild solution and, in the Newtonian limit, by comparing with the Newtonian N-body codes Gadget-2 and RAMSES. We then proceed with a simulation of large scale structure in a Universe with massive neutrinos where we study the gravitational slip induced by the neutrino shear stress. The code can be extended to include different kinds of dark energy or modified gravity models and going beyond the usually adopted quasi-static approximation. Our code is publicly available.

We derive the contribution to relativistic galaxy number count fluctuations from vector and tensor perturbations within linear perturbation theory. Our result is consistent with the the relativistic corrections to number counts due to scalar perturbation, where the Bardeen potentials are replaced with line-of-sight projection of vector and tensor quantities. Since vector and tensor perturbations do not lead to density fluctuations the standard density term in the number counts is absent. We apply our results to vector perturbations which are induced from scalar perturbations at second order and give numerical estimates of their contributions to the power spectrum of relativistic galaxy number counts.

Euclid is a European Space Agency medium class mission selected for launch in 2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to understand the origin of the accelerated expansion of the universe. Euclid will explore the expansion history of the universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

We modified the CLASS code in order to include relativistic galaxy number counts in spatially curved geometries; we present the formalism and study the effect of relativistic corrections on spatial curvature. The new version of the code is now publicly available. Using a Fisher matrix analysis, we investigate how measurements of the spatial curvature parameter $\Omega_K$ with future galaxy surveys are affected by relativistic effects, which influence observations of the large scale galaxy distribution. These effects include contributions from cosmic magnification, Doppler terms and terms involving the gravitational potential. As an application, we consider angle and redshift dependent power spectra, which are especially well suited for model independent cosmological constraints. We compute our results for a representative deep, wide and spectroscopic survey, and our results show the impact of relativistic corrections on the spatial curvature parameter estimation. We show that constraints on the curvature parameter may be strongly biased if, in particular, cosmic magnification is not included in the analysis. Other relativistic effects turn out to be subdominant in the studied configuration. We analyze how the shift in the estimated best-fit value for the curvature and other cosmological parameters depends on the magnification bias parameter, and find that significant biases are to be expected if this term is not properly considered in the analysis.

We introduce a new formalism to study perturbations of Hassan-Rosen bigravity theory, around general backgrounds for the two dynamical metrics. In particular, we derive the general expression for the mass term of the perturbations and we explicitly compute it for cosmological settings. We study tensor perturbations in a specific branch of bigravity using this formalism. We show that the tensor sector is affected by a late-time instability, which sets in when the mass matrix is no longer positive definite.

We estimate the sensitivity of future galaxy surveys to cosmological parameters, using the redshift dependent angular power spectra of galaxy number counts, $C_\ell(z_1,z_2)$, calculated with all relativistic corrections at first order in perturbation theory. We pay special attention to the redshift dependence of the non-linearity scale and present Fisher matrix forecasts for Euclid-like and DES-like galaxy surveys. We compare the standard $P(k)$ analysis with the new $C_\ell(z_1,z_2)$ method. We show that for surveys with photometric redshifts the new analysis performs significantly better than the $P(k)$ analysis. For spectroscopic redshifts, however, the large number of redshift bins which would be needed to fully profit from the redshift information, is severely limited by shot noise. We also identify surveys which can measure the lensing contribution and we study the monopole, $C_0(z_1,z_2)$.

We present some accurate and efficient computations of large scale structure observables, obtained with a modified version of the CLASS code which is made publicly available. This code includes all relativistic corrections and computes both the power spectrum Cl(z1,z2) and the corresponding correlation function xi(theta,z1,z2) in linear perturbation theory. For Gaussian initial perturbations, these quantities contain the full information encoded in the large scale matter distribution at the level of linear perturbation theory. We illustrate the usefulness of our code for cosmological parameter estimation through a few simple examples.

We compute the angular power spectra of the E-type and B-type lensing potentials for gravitational waves from inflation and for tensor perturbations induced by scalar perturbations. We derive the tensor-lensed CMB power spectra for both cases. We also apply our formalism to determine the linear lensing potential for a Bianchi I spacetime with small anisotropy.

We discuss the dominant terms of the relativistic galaxy number counts to second order in cosmological perturbation theory on sub-Hubble scales and on intermediate to large redshifts. In particular, we determine their contribution to the bispectrum. In addition to the terms already known from Newtonian second order perturbation theory, we find that there are a series of additional `lensing-like' terms which contribute to the bispectrum. We derive analytical expressions for the full leading order bispectrum and we evaluate it numerically for different configurations, indicating how they can be measured with upcoming surveys. In particular, the new `lensing-like' terms are not negligible for wide redshift bins and even dominate the bispectrum at well separated redshifts. This offers us the possibility to measure them in future surveys.

We investigate how the lensing potential can be measured tomographically with future galaxy surveys using their number counts. Such a measurement is an independent test of the standard $\Lambda$CDM framework and can be used to discern modified theories of gravity. We perform a Fisher matrix forecast based on galaxy angular-redshift power spectra, assuming specifications consistent with future photometric Euclid-like surveys and spectroscopic SKA-like surveys. For the Euclid-like survey we derive a fitting formula for the magnification bias. Our analysis suggests that the cross correlation between different redshift bins is very sensitive to the lensing potential such that the survey can measure the amplitude of the lensing potential at the same level of precision as other standard $\Lambda$CDM cosmological parameters.

In this paper we study the generation of primordial perturbations in a cosmological setting of bigravity during inflation. We consider a model of bigravity which can reproduce the $\Lambda$CDM background and large scale structure and a simple model of inflation with a single scalar field and a quadratic potential. Reheating is implemented with a toy-model in which the energy density of the inflaton is entirely dissipated into radiation. We present analytic and numerical results for the evolution of primordial perturbations in this cosmological setting. We find that the amplitude of tensor perturbations generated during inflation is sufficiently suppressed to avoid the effects of the tensor instability discovered in Refs.[1,2] which develops during the cosmological evolution in the physical sector. We argue that from a pure analysis of the tensor perturbations this bigravity model is compatible with present observations. However, we derive rather stringent limits on inflation from the vector and scalar sectors.

Numerical simulations are a versatile tool providing insight into the complicated process of structure formation in cosmology. This process is mainly governed by gravity, which is the dominant force on large scales. To date, a century after the formulation of general relativity, numerical codes for structure formation still employ Newton's law of gravitation. This approximation relies on the two assumptions that gravitational fields are weak and that they are only sourced by non-relativistic matter. While the former appears well justified on cosmological scales, the latter imposes restrictions on the nature of the "dark" components of the Universe (dark matter and dark energy) which are, however, poorly understood. Here we present the first simulations of cosmic structure formation using equations consistently derived from general relativity. We study in detail the small relativistic effects for a standard {\Lambda}CDM cosmology which cannot be obtained within a purely Newtonian framework. Our particle-mesh N-body code computes all six degrees of freedom of the metric and consistently solves the geodesic equation for particles, taking into account the relativistic potentials and the frame-dragging force. This conceptually clean approach is very general and can be applied to various settings where the Newtonian approximation fails or becomes inaccurate, ranging from simulations of models with dynamical dark energy or warm/hot dark matter to core collapse supernova explosions.

We show that at second order, ensemble averages of observables and directional averages do not commute due to gravitational lensing -- observing the same thing in many directions over the sky is not the same as taking an ensemble average. In principle this non-commutativity is significant for a variety of quantities that we often use as observables and can lead to a bias in parameter estimation. We derive the relation between the ensemble average and the directional average of an observable, at second order in perturbation theory. We discuss the relevance of these two types of averages for making predictions of cosmological observables, focusing on observables related to distances and magnitudes. In particular, we show that the ensemble average of the distance in a given observed direction is increased by gravitational lensing, whereas the directional average of the distance is decreased. For a generic observable, there exists a particular function of the observable that is not affected by second-order lensing perturbations. We also show that standard areas have an advantage over standard rulers, and we discuss the subtleties involved in averaging in the case of supernova observations.

It has recently been shown that second-order corrections to the background distance-redshift relation can build up significantly at large redshifts, due to an aggregation of gravitational lensing events. This shifts the expectation value of the distance to the CMB by 1%. In this paper we show that this shift is already properly accounted for in standard CMB analyses. We clarify the role that the area distance to the CMB plays in the presence of second-order lensing corrections.

Three relativistic particles in addition to the photon are detected in the cosmic microwave background (CMB). In the standard model of cosmology, these are interpreted as the three neutrino species. However, at the time of CMB-decoupling, neutrinos are not only relativistic but they are also freestreaming. Here, we investigate, whether the CMB is sensitive to this defining feature of neutrinos, or whether the CMB-data allow to replace neutrinos with a relativistic fluid. We show that free streaming particles are preferred over a relativistic perfect fluid with $\Delta\chi^2\simeq 21$. We also study the possibility to replace the neutrinos by a viscous fluid and find that a relativistic viscous fluid with either the standard values $c_{\rm eff}^2=c_{\rm vis}^2=1/3$ or best fit values for $c_{\rm eff}^2$ and $c_{\rm vis}^2$ has $\Delta\chi^2=20$ and thus cannot provide a good fit to present CMB data either.

We determine the number counts to second order in cosmological perturbation theory in the Poisson gauge and allowing for anisotropic stress. The calculation is performed using an innovative approach based on the recently proposed "geodesic light-cone" gauge. This allows us to determine the number counts in a purely geometric way, without using Einstein's equation. The result is valid for general dark energy models and (most) modified gravity models. We then evaluate numerically the relevant contributions to the number counts bispectrum. In particular we consider the terms involving the density, redshift space distortion and lensing.

This review describes the discovery of the cosmic microwave background radiation in 1965 and its impact on cosmology in the 50 years that followed. This discovery has established the Big Bang model of the Universe and the analysis of its fluctuations has confirmed the idea of inflation and led to the present era of precision cosmology. I discuss the evolution of cosmological perturbations and their imprint on the CMB as temperature fluctuations and polarization. I also show how a phase of inflationary expansion generates fluctuations in the spacetime curvature and primordial gravitational waves. In addition I present findings of CMB experiments, from the earliest to the most recent ones. The accuracy of these experiments has helped us to estimate the parameters of the cosmological model with unprecedented precision so that in the future we shall be able to test not only cosmological models but General Relativity itself on cosmological scales.

We investigate the cosmological dependence and the constraining power of large-scale galaxy correlations, including all redshift-distortions, wide-angle, lensing and gravitational potential effects on linear scales. We analyze the cosmological information present in the lensing convergence and in the gravitational potential terms describing the so-called "relativistic effects," and we find that, while smaller than the information contained in intrinsic galaxy clustering, it is not negligible. We investigate how neglecting them does bias cosmological measurements performed by future spectroscopic and photometric large-scale surveys such as SKA and Euclid. We perform a Fisher analysis using the CLASS code, modified to include scale-dependent galaxy bias and redshift-dependent magnification and evolution bias. Our results show that neglecting relativistic terms introduces an error in the forecasted precision in measuring cosmological parameters of the order of a few tens of percent, in particular when measuring the matter content of the Universe and primordial non-Gaussianity parameters. Therefore, we argue that radial correlations and integrated relativistic terms need to be taken into account when forecasting the constraining power of future large-scale number counts of galaxy surveys.

In this paper we study gravitational wave perturbations in a cosmological setting of bigravity which can reproduce the {\Lambda}CDM background and large scale structure. We show that in general gravitational wave perturbations are unstable and only for very fine tuned initial conditions such a cosmology is viable. We quantify this fine tuning. We argue that similar fine tuning is also required in the scalar sector in order to prevent the tensor instability to be induced by second order scalar perturbations. Finally, we show that due to this power law instability, models of bigravity can lead to a large tensor to scalar ratio even for low scale inflation.

The large-scale homogeneity and isotropy of the universe is generally thought to imply a well defined background cosmological model. It may not. Smoothing over structure adds in an extra contribution, transferring power from small scales up to large. Second-order perturbation theory implies that the effect is small, but suggests that formally the perturbation series may not converge. The amplitude of the effect is actually determined by the ratio of the Hubble scales at matter-radiation equality and today - which are entirely unrelated. This implies that a universe with significantly lower temperature today could have significant backreaction from more power on small scales, and so provides the ideal testing ground for understanding backreaction. We investigate this using two different N-body numerical simulations - a 3D Newtonian and a 1D simulation which includes all relevant relativistic effects. We show that while perturbation theory predicts an increasing backreaction as more initial small-scale power is added, in fact the virialisation of structure saturates the backreaction effect at the same level independently of the equality scale. This implies that backreaction is a small effect independently of initial conditions. Nevertheless, it may still contribute at the percent level to certain cosmological observables and therefore it cannot be neglected in precision cosmology.

The success of precision cosmology depends not only on accurate observations, but also on the theoretical model - which must be understood to at least the same level of precision. Subtle relativistic effects can lead to biased measurements if they are neglected. One such effect gives a systematic shift in the distance-redshift relation away from its background value, due to the non- linear relativistic conservation of total photon flux. We also show directly how this shift follows from a fully relativistic analysis of the geodesic deviation equation. We derive the expectation value of the shift using second-order perturbations about a concordance background, and show that the distance to last scattering is increased by 1%. We argue that neglecting this shift could lead to a significant bias in the background cosmological parameters, because it alters the meaning of the background model. A naive adjustment of CMB parameter estimation if this shift is really a correction to the background would raise the H0 value inferred from the CMB by 5%, potentially removing the tension with local measurements of H0. Other CMB parameters which depend on the distance would also be shifted by ~1sigma when combined with local H0 data. While our estimations rely on a simplistic analysis, they nevertheless illustrate that accurately defining the background model in terms of the expectation values of observables is critical when we aim to determine the model parameters at the sub-percent level.

We present a framework for general relativistic N-body simulations in the regime of weak gravitational fields. In this approach, Einstein's equations are expanded in terms of metric perturbations about a Friedmann-Lema\^itre background, which are assumed to remain small. The metric perturbations themselves are only kept to linear order, but we keep their first spatial derivatives to second order and treat their second spatial derivatives as well as sources of stress-energy fully non-perturbatively. The evolution of matter is modelled by an N-body ensemble which can consist of free-streaming nonrelativistic (e.g. cold dark matter) or relativistic particle species (e.g. cosmic neutrinos), but the framework is fully general and also allows for other sources of stress-energy, in particular additional relativistic sources like modified-gravity models or topological defects. We compare our method with the traditional Newtonian approach and argue that relativistic methods are conceptually more robust and flexible, at the cost of a moderate increase of numerical difficulty. However, for a LambdaCDM cosmology, where nonrelativistic matter is the only source of perturbations, the relativistic corrections are expected to be small. We quantify this statement by extracting post-Newtonian estimates from Newtonian N-body simulations.

Local measurements of the Hubble expansion rate are affected by structures like galaxy clusters or voids. Here we present a fully relativistic treatment of this effect, studying how clustering modifies the mean distance (modulus)-redshift relation and its dispersion in a standard $\Lambda$CDM universe. The best estimates of the local expansion rate stem from supernova observations at small redshifts (0.01<z<0.1). It is interesting to compare these local measurements with global fits to data from cosmic microwave background anisotropies. In particular, we argue that cosmic variance (i.e. the effects of the local structure) is of the same order of magnitude as the current observational errors and must be taken into account in local measurements of the Hubble expansion rate.

If the B-mode signal in the CMB polarization seen by the BICEP2 experiment is confirmed, it has dramatic implications for models of inflation. The result is also in tension with Planck limits on standard inflationary models. It is therefore important to investigate whether this signal can arise from alternative sources. If so, this could lessen the pressure on inflationary models and the tension with Planck data. We investigate whether vector and tensor modes from primordial magnetic fields can explain the signal. We find that in principle, magnetic fields generated during inflation can indeed produce the required B-mode, for a suitable range of energy scales of inflation. In this case, the primordial gravitational wave amplitude is negligible, so that there is no tension with Planck and no problems posed for current inflationary models. However, the simplest magnetic model is in tension with Planck limits on non-Gaussianity in the trispectrum. It may be possible to fine-tune the magnetogenesis model so that this non-Gaussianity is suppressed. Alternatively, a weaker magnetic field can pass the non-Gaussianity constraints and allow the primordial tensor mode to be reduced to r~0.09, thus removing the tension with Planck data and alleviating the problems with simple inflationary models.

Distance measurements are usually thought to probe the background metric of the universe, but in reality the presence of perturbations will lead to deviations from the result expected in an exactly homogeneous and isotropic universe. At least in principle the presence of perturbations could even explain the observed distance-redshift relation without the need for dark energy. In this paper we re-investigate a toy model where perturbations are plane symmetric, and for which exact solutions are known in the fluid limit. However, if perturbations are large, shell-crossing occurs and the fluid approximation breaks down. This prevents the study of the most interesting cases. Here we use a general-relativistic N-body simulation that does not suffer from this problem and which allows us to go beyond previous works. We show that even for very large plane-symmetric perturbations we are not able to mimic the observed distance-redshift relation. We also discuss how the synchronous comoving gauge breaks down when shell-crossing occurs, while metric perturbations in the longitudinal gauge remain small. For this reason the longitudinal (Newtonian) gauge appears superior for relativistic N-body simulations of large-scale structure formation.

We show that self-ordering scalar fields (SOSF), i.e. non-topological cosmic defects arising after a global phase transition, cannot explain the B-mode signal recently announced by BICEP2. We compute the full $C_\ell^{B}$ angular power spectrum of B-modes due to the vector and tensor perturbations of SOSF, modeled in the large-N limit of a spontaneous broken global O(N) symmetry. We conclude that the low-$\ell$ multipoles detected by BICEP2 cannot be due mainly to SOSF, since they have the wrong spectrum at low multipoles. As a byproduct we derive the first cosmological constraints on this model, showing that the BICEP2 B-mode polarization data admits at most a 2-3% contribution from SOSF in the temperature anisotropies, similar to (but somewhat tighter than) the recently studied case of cosmic strings.

We determine the best-fit values and confidence limits for dynamical dark energy parameters together with other cosmological parameters on the basis of different datasets which include WMAP9 or Planck-2013 results on CMB anisotropy, BAO distance ratios from recent galaxy surveys, magnitude-redshift relations for distant SNe Ia from SNLS3 and Union2.1 samples and the HST determination of the Hubble constant. We use a Markov Chain Monte Carlo routine to map out the likelihood in the multi-dimensional parameter space. We show that the most precise determination of cosmological parameters with the narrowest confidence limits is obtained for the Planck{+}HST{+}BAO{+}SNLS3 dataset. The best-fit values and 2$\sigma$ confidence limits for cosmological parameters in this case are $\Omega_{de}=0.718\pm0.022$, $w_0=-1.15^{+0.14}_{-0.16}$, $c_a^2=-1.15^{+0.02}_{-0.46}$, $\Omega_bh^2=0.0220\pm0.0005$, $\Omega_{cdm}h^2=0.121\pm0.004$, $h=0.713\pm0.027$, $n_s=0.958^{+0.014}_{-0.010}$, $A_s=(2.215^{+0.093}_{-0.101})\cdot10^{-9}$, $\tau_{rei}=0.093^{+0.022}_{-0.028}$. For this dataset, the $\Lambda$CDM model is just outside the 2$\sigma$ confidence region, while for the dataset WMAP9{+}HST{+}BAO{+}SNLS3 the $\Lambda$CDM model is only 1$\sigma$ away from the best fit. The tension in the determination of some cosmological parameters on the basis of two CMB datasets WMAP9 and Planck-2013 is highlighted.

A theory of massive gravity depends on a non-dynamical 'reference metric' f_{\mu\nu} which is often taken to be the flat Minkowski metric. In this paper we examine the theory of perturbations on a background with metric g_{\mu\nu} which does not coincide with the reference metric f_{\mu\nu}. We derive the mass term for general perturbations on this background and show that it generically is not of the form of the Fierz-Pauli mass term. We explicitly compute it for some cosmological situations.

We determine the full C_l spectra and correlation functions of the temperature and polarization anisotropies in the CMB, generated by a source modeled by the large N limit of spontaneously broken global O(N)-theories. We point out a problem in the standard approach of treating the radiation-matter transition by interpolating the eigenvectors of the unequal-time correlators of the source energy-momentum tensor. This affects the CMB predictions from all type of cosmic defects. We propose a method to overcome this difficulty, and find that in the large-N global model that we study, differences in the final CMB power spectra amplitudes reach up to 25%, when compared to implementations of the eigenvector interpolation technique. We discuss as well how to optimally search for the contribution in the CMB from active sources such as cosmic defects, in experiments like Planck, COrE and PRISM.

We develop a formalism for General Relativistic N-body simulations in the weak field regime, suitable for cosmological applications. The problem is kept tractable by retaining the metric perturbations to first order, the first derivatives to second order and second derivatives to all orders, thus taking into account the most important nonlinear effects of Einstein gravity. It is also expected that any significant "backreaction" should appear at this order. We show that the simulation scheme is feasible in practice by implementing it for a plane-symmetric situation and running two test cases, one with only cold dark matter, and one which also includes a cosmological constant. For these plane-symmetric situations, the deviations from the usual Newtonian N-body simulations remain small and, apart from a non-trivial correction to the background, can be accurately estimated within the Newtonian framework. The correction to the background scale factor, which is a genuine "backreaction" effect, can be robustly obtained with our algorithm. Our numerical approach is also naturally suited for the inclusion of extra relativistic fields and thus for dark energy or modified gravity simulations.

We argue that the result presented in "Origin of cosmic magnetic fields" by L. Campanelli [arXiv:1304.6534] is unphysical.

We study the Horndeski vector-tensor theory that leads to second order equations of motion and contains a non-minimally coupled abelian gauge vector field. This theory is remarkably simple and consists of only 2 terms for the vector field, namely: the standard Maxwell kinetic term and a coupling to the dual Riemann tensor. Furthermore, the vector sector respects the U(1) gauge symmetry and the theory contains only one free parameter, M^2, that controls the strength of the non-minimal coupling. We explore the theory in a de Sitter spacetime and study the presence of instabilities and show that it corresponds to an attractor solution in the presence of the vector field. We also investigate the cosmological evolution and stability of perturbations in a general FLRW spacetime. We find that a sufficient condition for the absence of ghosts is M^2>0. Moreover, we study further constraints coming from imposing the absence of Laplacian instabilities. Finally, we study the stability of the theory in static and spherically symmetric backgrounds (in particular, Schwarzschild and Reissner-Nordstr\"om-de Sitter). We find that the theory, quite generally, do have ghosts or Laplacian instabilities in regions of spacetime where the non-minimal interaction dominates over the Maxwell term. We also calculate the propagation speed in these spacetimes and show that superluminality is a quite generic phenomenon in this theory.

The cosmic microwave background (CMB) temperature anisotropies from primordial magnetic fields are studied. In addition to the known passive and the compensated mode we discuss an inflationary magnetic mode in the curvature perturbation, present when magnetic fields are generated during inflation. This mode is absent if the generation mechanism is causal, e.g. a phase transition. We compute and discuss the effect of this mode on the observed CMB anisotropy spectrum, in comparison with the passive and compensated ones. We find that it dominates the CMB anisotropy, and consequently leads to stronger constraints on the amplitude $B_{1\,{\rm Mpc}}$ and spectral index $n_B$ of the magnetic field than what usually found in CMB analyses from the compensated mode. This happens in particular for spectral indexes $n_B>-3$: the inflationary magnetic mode is always scale invariant, therefore through this mode even a magnetic field with a spectrum which is not scale invariant can leave a detectable signal in the CMB at large scales.

We analyze the possibility to distinguish between quintessence and phantom scalar field models of dark energy using observations of luminosity distance moduli of SNe Ia, CMB anisotropies and polarization, matter density perturbations and baryon acoustic oscillations. Among the present observations only Planck data on CMB anisotropy and SDSS DR9 data on baryon acoustic oscillations may be able to decide between quintessence or phantom scalar field models, however for each model a set of best-fit parameters exists, which matches all data with similar goodness of fit. We compare the relative differences of best-fit model predictions with observational uncertainties for each type of data and we show that the accuracy of SNe Ia luminosity distance data is far from the one necessary to distinguish these types of dark energy models, while the CMB data (WMAP, ACT, SPT and especially Planck) are close to being able to reliably distinguish them. Also an improvement of the large-scale structure data (future releses of SDSS BOSS and e.g. Euclid or BigBOSS) will enable us to surely decide between quintessence and phantom dark energy.

We review the possible mechanisms for the generation of cosmological magnetic fields, discuss their evolution in an expanding Universe filled with the cosmic plasma and provide a critical review of the literature on the subject. We put special emphasis on the prospects for observational tests of the proposed cosmological magnetogenesis scenarios using radio and gamma-ray astronomy and ultra high energy cosmic rays. We argue that primordial magnetic fields are observationally testable. They lead to magnetic fields in the intergalactic medium with magnetic field strength and correlation length in a well defined range. We also state the unsolved questions in this fascinating open problem of cosmology and propose future observations to address them.

We show that a longitudinal gauge degree of freedom for a vector field is equivalent to a Pais-Uhlenbeck scalar field. With the help of this equivalence, we can determine natural interactions of this field with scalars and fermions. Since the theory has a global U(1) symmetry, we have the usual conserved current of the charged fields, thanks to which the dynamics of the scalar field is not modified by the interactions. We use this fact to consistently quantize the theory even in the presence of interactions. We argue that such a degree of freedom can only be excited by gravitational effects like the inflationary era of the early universe and may play the role of dark energy in the form of an effective cosmological constant whose value is linked to the inflation scale.

We consider a model of inflation which has recently been proposed in the literature and where inflation is induced by corrections to the energy density coming from the non-commutativity of spacetime. We show that the very rapid inflationary expansion typical of this model is responsible for a burst of particle production which ends inflation and leads to a radiation-dominated phase. We analytically estimate the energy density of these particles and we confront the results with more precise numerical calculations. We estimate the number of inflationary e-folds before the back-reaction of the radiation energy density overcomes the non-commutative effects and terminate inflation naturally.

The dynamics of expansion and large scale structure formation of the Universe are analyzed for models with dark energy in the form of a phantom scalar field which initially mimics a $\Lambda$-term and evolves slowly to the Big Rip singularity. The discussed model of dark energy has three parameters -- the density and the equation of state parameter at the current epoch, $\Omega_{de}$ and $w_0$, and the asymptotic value of the equation of state parameter at $a\rightarrow\infty$, $c_a^2$. Their best-fit values are determined jointly with all other cosmological parameters by the MCMC method using observational data on CMB anisotropies and polarization, SNe Ia luminosity distances, BAO measurements and more. Similar computations are carried out for $\Lambda$CDM and a quintessence scalar field model of dark energy. It is shown that the current data slightly prefer the phantom model, but the differences in the maximum likelihoods are not statistically significant. It is also shown that the phantom dark energy with monotonically increasing density in future will cause the decay of large scale linear matter density perturbations due to the gravitational domination of dark energy perturbations long before the Big Rip singularity.

We argue that from observations alone, only the transverse power spectrum $C_\ell(z_1,z_2)$ and the corresponding correlation function $\xi(\theta,z_1,z_2)$ can be measured and that these contain the full three dimensional information. We determine the two point galaxy correlation function at linear order in perturbation theory. Redshift space distortions are taken into account for arbitrary angular and redshift separations. We discuss the shape of the longitudinal and the transversal correlation functions which are very different from each other and from their real space counterpart. We then go on and suggest how to measure both, the Hubble parameter, $H(z)$, and the angular diameter distance, $D_A(z)$, separately from these correlation functions and perform an Alcock-Paczynski test.

We compute the vector and tensor contributions to the luminosity distance fluctuations in first order perturbation theory and we expand them in spherical harmonics. This work presents the formalism with a first application to a stochastic background of primordial gravitational waves.

We compute the contribution to the scalar metric perturbations from large-scale magnetic fields which are generated during inflation. We show that apart from the usual passive and compensated modes, the magnetic fields also contribute to the constant mode from inflation. This is different from the causal (post inflationary) generation of magnetic fields where such a mode is absent and it might lead to significant, non-Gaussian CMB anisotropies.

In this paper we compute the quantity which is truly measured in a large galaxy survey. We take into account the effects coming from the fact that we actually observe galaxy redshifts and sky positions and not true spatial positions. Our calculations are done within linear perturbation theory for both the metric and the observer velocities but they can be used for non-linear matter power spectra. We shall see that the complications due to the fact that we only observe on our background lightcone and that we do not truly know the distance of the observed galaxy, but only its redshift is not only an additional difficulty, but even more a new opportunity for future galaxy surveys.

We examine the mode functions of the electromagnetic field on spherically symmetric backgrounds with special attention to the subclass which allows for the foliation as open Friedmann-Lemaitre (FL) spacetime. It is well-known that in certain scalar field theories on open FL background there can exist so-called supercurvature modes, their existence depending on parameters of the theory. Looking at specific open universe models, such as open inflation and the Milne universe, we find that no supercurvature modes are present in the spectrum of the electromagnetic field. This excludes the possibility for superadiabatic evolution of cosmological magnetic fields within these models without relying on new physics or breaking the conformal invariance of electromagnetism.

The generation of helical magnetic fields during single field inflation due to an axial coupling of the electromagnetic field to the inflaton is discussed. We find that such a coupling always leads to a blue spectrum of magnetic fields during slow roll inflation. Though the helical magnetic fields further evolve during the inverse cascade in the radiation era after inflation, we conclude that the magnetic fields generated by such an axial coupling can not lead to observed field strength on cosmologically relevant scales.

This paper has been withdrawn by the authors due to a mistake in the calculations leading to an erroneous conclusion.

We study the distance-redshift relation in a universe filled with 'walls' of pressure-less dust separated by under dense regions. We show that as long as the density contrast of the walls is small, or the diameter of the under dense regions is much smaller than the Hubble scale, the distance-redshift relation remains close to what is obtained in a Friedmann universe. However, when arbitrary density contrasts are allowed, every prescribed distance-redshift relation can be reproduced with such models.

We investigate the statistics of the cosmic microwave background using the Kolmogorov-Smirnov test. We show that, when we correctly de-correlate the data, the partition function of the Kolmogorov stochasticity parameter is compatible with the Kolmogorov distribution and, contrary to previous claims, the CMB data are compatible with Gaussian fluctuations with the correlation function given by standard Lambda-CDM. We then use the Kolmogorov-Smirnov test to derive upper bounds on residual point source power in the CMB, and indicate the promise of this statistics for further datasets, especially Planck, to search for deviations from Gaussianity and for detecting point sources and Galactic foregrounds.

We derive an expression for the luminosity distance in a perturbed Friedmann universe. We define the correlation function and the power spectrum of the luminosity distance fluctuations and express them in terms of the initial spectrum of the Bardeen potential. We present semi-analytical results for the case of a pure CDM universe. We argue that the luminosity distance power spectrum represents a new observational tool which can be used to determine cosmological parameters. In addition, our results shed some light into the debate whether second order small scale fluctuations can mimic an accelerating universe.

We provide a simple transfer function that determines the effect of an early matter dominated era on the gravitational wave background and show that a large class of compactifications of superstring theory might be tested by observations of the gravitational wave background from inflation. For large enough reheating temperatures $> 10^9 \GeV$ the test applies to all models containing at least one scalar with mass $< 10^{12}\GeV$ that acquires a large initial oscillation amplitude after inflation and has only gravitational interaction strength, i.e., a field with the typical properties of a modulus.

In this paper I discuss what we truly know about dark energy. I shall argue that up to date our single indication for the existence of dark energy comes from distance measurements and their relation to redshift. Supernovae, CMB anisotropies and observations of baryon acoustic oscillations, they all simply tell us that the observed distance to a given redshift is larger than the one expected from a Friedmann Lemaitre universe with matter only and the locally measured Hubble parameter.

We constrain the parameters of dynamical dark energy in the form of a classical or tachyonic scalar field with barotropic equation of state jointly with other cosmological ones using the combined datasets which include the CMB power spectra from WMAP7, the baryon acoustic oscillations in the space distribution of galaxies from SDSS DR7, the power spectrum of luminous red galaxies from SDSS DR7 and the light curves of SN Ia from 2 different compilations: Union2 (SALT2 light curve fitting) and SDSS (SALT2 and MLCS2k2 light curve fittings). It has been found that the initial value of dark energy equation of state parameter is constrained very weakly by most of the data while the rest of main cosmological parameters are well constrained: their likelihoods and posteriors are similar, have the forms close to Gaussian (or half-Gaussian) and their confidential ranges are narrow. The most reliable determinations of the best fitting value and $1\sigma$ confidence range for the initial value of dark energy equation of state parameter were obtained from the combined datasets including SN Ia data from the full SDSS compilation with MLCS2k2 fitting of light curves. In all such cases the best fitting value of this parameter is lower than the value of corresponding parameter for current epoch. Such dark energy loses its repulsive properties and in future the expansion of the Universe will change into contraction. We also perform an error forecast for the Planck mock data and show that they narrow essentially the confidential ranges of cosmological parameters values, moreover, their combination with SN SDSS compilation with MLCS2k2 light curve fitting may exclude the fields with initial equation of state parameter $>-0.1$ at 2$\sigma$ confidential level.

The fitting formula for the location of the first acoustic peak in the matter power spectrum is revised. We discuss the physics that leads to baryon acoustic oscillations: the recombination history, the tight coupling approximation and the velocity overshoot effect. A new fitting formula is proposed, which is in accordance within 5% with numerical results for a suitable range of cosmological parameters, whereas previous results yield deviations of up to 20%. The crucial improvement turns out to be the accuracy of the recombination history.

We study a consistent infrared and ultraviolet regularization scheme for the cosmological perturbations. The infrared divergences are cured by assuming that the Universe undergoes a transition between a non-singular pre-inflationary, radiation-dominated phase and a slow-roll inflationary evolution. The ultraviolet divergences are eliminated via adiabatic subtraction. A consistent regularization of the field fluctuations through this transition is obtained by performing a mode matching for both the gauge invariant Mukhanov variable and its adiabatic expansion. We show that these quantities do not generate ultraviolet divergences other than the standard ones, when evolving through the matching time. We also show how the de Witt-Schwinger expansion, which can be used to construct the counter-terms regularizing the ultraviolet divergences, ceases to be valid well before horizon exit of the scales of interest. Thus, such counter-terms should not be used beyond the time of the horizon exit so it is unlikely that the observed power spectrum is modified by adiabatic subtraction as claimed in the literature. However, the infrared regularization might have an impact on the observed spectrum, and we briefly discuss this possibility.

n this paper we study the effect of biasing on the power spectrum at large scales. We show that even though non-linear biasing does introduce a white noise contribution on large scales, the $P(k)\propto k^n$ behavior of the matter power spectrum on large scales may still be visible and above the white noise for about one decade. We show, that the Kaiser biasing scheme which leads to linear bias of the correlation function on {\em large} scales, also generates a linear bias of the {\rm power spectrum} on rather small scales. This is a consequence of the divergence on small scales of the pure Harrison-Zeldovich spectrum. However, biasing becomes k-dependent when we damp the underlying power spectrum on small scales. We also discuss the effect of biasing on the baryon acoustic oscillations.

We study the generation of helical magnetic fields during single field inflation induced by an axial coupling of the electromagnetic field to the inflaton. During slow roll inflation, we find that such a coupling always leads to a blue spectrum with $B^2(k) \propto k$, as long as the theory is treated perturbatively. The magnetic energy density at the end of inflation is found to be typically too small to backreact on the background dynamics of the inflaton. We also show that a short deviation from slow roll does not result in strong modifications to the shape of the spectrum. We calculate the evolution of the correlation length and the field amplitude during the inverse cascade and viscous damping of the helical magnetic field in the radiation era after inflation. We conclude that except for low scale inflation with very strong coupling, the magnetic fields generated by such an axial coupling in single field slow roll inflation with perturbative coupling to the inflaton are too weak to provide the seeds for the observed fields in galaxies and clusters.

We study a homogeneous and nearly-isotropic Universe permeated by a homogeneous magnetic field. Together with an isotropic fluid, the homogeneous magnetic field, which is the primary source of anisotropy, leads to a plane-symmetric Bianchi I model of the Universe. However, when free-streaming relativistic particles are present, they generate an anisotropic pressure which counteracts the one from the magnetic field such that the Universe becomes isotropized. We show that due to this effect, the CMB temperature anisotropy from a homogeneous magnetic field is significantly suppressed if the the neutrino masses are smaller than 0.3 eV.

We present a new and especially powerful signature of cosmic strings and other topological or non-topological defects in the polarization of the cosmic microwave background (CMB). We show that even if defects contribute 1% or less in the CMB temperature anisotropy spectrum, their signature in the local $\tilde{B}$-polarization correlation function at angular scales of tens of arc minutes is much larger than that due to gravitational waves from inflation, even if the latter contribute with a ratio as big as $r\simeq 0.1$ to the temperature anisotropies. We show that when going from non-local to local $\tilde{B}$-polarization, the ratio of the defect signal-to-noise with respect to the inflationary value increases by about an order of magnitude. Proposed B-polarization experiments, with a good sensitivity on arcminute scales, may either detect a contribution from topological defects produced after inflation or place stringent limits on them. Even Planck should be able to improve present constraints on defect models by at least an order of magnitude, to the level of $\ep <10^{-7}$. A future full-sky experiment like CMBpol, with polarization sensitivities of the order of $1\mu$K-arcmin, will be able to constrain the defect parameter $\ep=Gv^2$ to a few $\times10^{-9}$, depending on the defect model.

If the cosmological QCD phase transition is strongly first order and lasts sufficiently long, it generates a background of gravitational waves which may be detected via pulsar timing experiments. We estimate the amplitude and the spectral shape of such a background and we discuss its detectability prospects.

We analyse CMB data in a manner which is as model-independent as possible. We encode the effects of late-time cosmology into a single parameter which determines the distance to the last scattering surface. We exclude low multipoles $\ell<40$ from the analysis. We consider the WMAP5 and ACBAR data. We obtain the cosmological parameters $100\omega_b =2.13\pm 0.05$, $\omega_c=0.124\pm 0.007$, $n_s=0.93\pm 0.02$ and $\theta_A=0.593\pm 0.001$ degrees (68% C.L.). The last number is the angular scale subtended by the sound horizon at decoupling. There is a systematic shift in the parameters as more low $\ell$ data is omitted, towards smaller values of $\omega_b$ and $n_s$ and larger values of $\omega_c$. The scale $\theta_A$ remains stable and very well determined.

I discuss the generation of a stochastic background of gravitational waves during a first order phase transition. I present simple general arguments which explain the main features of the gravitational wave spectrum like the $k^3$ power law growth on large scales and a estimate for the peak amplitude. In the second part I concentrate on the electroweak phase transition and argue that the nucleosynthesis bound on its gravitational wave background seriously limits seed magnetic fields which may have been generated during this transition.

We analytically derive the spectrum of gravitational waves due to magneto-hydrodynamical turbulence generated by bubble collisions in a first-order phase transition. In contrast to previous studies, we take into account the fact that turbulence and magnetic fields act as sources of gravitational waves for many Hubble times after the phase transition is completed. This modifies the gravitational wave spectrum at large scales. We also model the initial stirring phase preceding the Kolmogorov cascade, while earlier works assume that the Kolmogorov spectrum sets in instantaneously. The continuity in time of the source is relevant for a correct determination of the peak position of the gravitational wave spectrum. We discuss how the results depend on assumptions about the unequal-time correlation of the source and motivate a realistic choice for it. Our treatment gives a similar peak frequency as previous analyses but the amplitude of the signal is reduced due to the use of a more realistic power spectrum for the magneto-hydrodynamical turbulence. For a strongly first-order electroweak phase transition, the signal is observable with the space interferometer LISA.

Gravitational waves were copiously produced in the early Universe whenever the processes taking place were sufficiently violent. The spectra of several of these gravitational wave backgrounds on subhorizon scales have been extensively studied in the literature. In this paper we analyze the shape and amplitude of the gravitational wave spectrum on scales which are superhorizon at the time of production. Such gravitational waves are expected from the self ordering of randomly oriented scalar fields which can be present during a thermal phase transition or during preheating after hybrid inflation. We find that, if the gravitational wave source acts only during a small fraction of the Hubble time, the gravitational wave spectrum at frequencies lower than the expansion rate at the time of production behaves as $\Omega_{\rm GW}(f) \propto f^3$ with an amplitude much too small to be observable by gravitational wave observatories like LIGO, LISA or BBO. On the other hand, if the source is active for a much longer time, until a given mode which is initially superhorizon ($k\eta_* \ll 1$), enters the horizon, for $k\eta \gtrsim 1$, we find that the gravitational wave energy density is frequency independent, i.e. scale invariant. Moreover, its amplitude for a GUT scale scenario turns out to be within the range and sensitivity of BBO and marginally detectable by LIGO and LISA. This new gravitational wave background can compete with the one generated during inflation, and distinguishing both may require extra information.

We discuss the impact of adiabatic renormalization on the power spectrum of scalar and tensor perturbations from inflation. We show that adiabatic regularization is ambiguous as it leads to very different results, for different adiabatic subtraction schemes, both in the range $v\equiv k/(aH) \gsim 0.1$ and in the infrared regime. All these schemes agree in the far ultraviolet, $v\gg 1$. Therefore, we argue that in the far infrared regime, $v\ll 1$, the adiabatic expansion is no longer valid, and the unrenormalized spectra are the physical, measurable quantities. These findings cast some doubt on the validity of the adiabatic subtraction at horizon exit, $v=1$, to determine the perturbation spectra from inflation which has recently advocated in the literature.

In recent work by two of us, [Durrer & Ruser, PRL 99, 071601 (2007); Ruser & Durrer PRD 76, 104014 (2007)], graviton production due to a moving spacetime boundary (braneworld) in a five dimensional bulk has been considered. In the same way as the presence of a conducting plate modifies the electromagnetic vacuum, the presence of a brane modifies the graviton vacuum. As the brane moves, the time dependence of the vacuum leads to particle creation via the so called 'dynamical Casimir effect'. In our previous work a term in the boundary condition which is linear in the brane velocity has been neglected. In this work we develop a new approach which overcomes this approximation. We show that the previous results are not modified if the brane velocity is low.

Gravitational wave production induces a strong constraint on the amplitude of a primordial magnetic field. It has been shown that the nucleosynthesis bound for a stochastic gravitational wave background implies that causally generated fields cannot have enough power on large scales to provide the seeds necessary for the observed magnetic fields in galaxies and clusters, even by the most optimistic dynamo amplification. Magnetic fields generated at inflation can have high enough amplitude only if their spectrum is very red. Here we show that helicity, which leads to an inverse cascade, can mitigate these limits. In particular, we find that helical fields generated at the QCD phase transition or at inflation with red spectrum are possible seeds for the dynamo. Helical fields generated at the electroweak phase transition are instead excluded as seeds at large scales. We also calculate the spectrum of gravitational waves generated by helical magnetic fields.

We investigate models where structure formation is initiated by scaling seeds: We consider rapidly expanding relativistic shells of energy and show that they can fit current CMB and large scale structure data if they expand with super-luminal velocities. These acausally expanding shells provide a viable alternative to inflation for cosmological structure formation with the same minimal number of parameters to characterize the initial fluctuations. Causally expanding shells alone cannot fit present data. Hybrid models where causal shells and inflation are mixed also provide good fits.

In this paper we discuss some general aspects of the gravitational wave background arising from post-inflationary short-lasting cosmological events such as phase transitions. We concentrate on the physics which determines the shape and the peak frequency of the gravitational wave spectrum. We then apply our general findings to the case of bubble collisions during a first order phase transition and compare different results in the recent literature.

We discuss the creation of massless particles in a Universe, which transits from a radiation-dominated era to any other expansion law. We calculate in detail the generation of gravitons during the transition to a matter dominated era. We show that the resulting gravitons generated in the standard radiation/matter transition are negligible. We use our result to constrain one or more previous matter dominated era, or any other expansion law, which may have taken place in the early Universe. We also derive a general formula for the modification of a generic initial graviton spectrum by an early matter dominated era.

If our Universe is a 3+1 brane in a warped 4+1 dimensional bulk so that its expansion can be understood as the motion of the brane in the bulk, the time dependence of the boundary conditions for arbitrary bulk fields can lead to particle creation via the dynamical Casimir effect. In this talk I report results for the simplest such scenario, when the only particle in the bulk is the graviton and the bulk is the 5 dimensional anti-de Sitter spacetime.

Gravitational wave production from bubble collisions was calculated in the early nineties using numerical simulations. In this paper, we present an alternative analytic estimate, relying on a different treatment of stochasticity. In our approach, we provide a model for the bubble velocity power spectrum, suitable for both detonations and deflagrations. From this, we derive the anisotropic stress and analytically solve the gravitational wave equation. We provide analytical formulae for the peak frequency and the shape of the spectrum which we compare with numerical estimates. In contrast to the previous analysis, we do not work in the envelope approximation. This paper focuses on a particular source of gravitational waves from phase transitions. In a companion article, we will add together the different sources of gravitational wave signals from phase transitions: bubble collisions, turbulence and magnetic fields and discuss the prospects for probing the electroweak phase transition at LISA.

The energy momentum tensor of a magnetic field always contains a spin-2 component in its anisotropic stress and therefore generates gravitational waves. It has been argued in the literature (Caprini & Durrer \cite{CD}) that this gravitational wave production can be very strong and that back-reaction cannot be neglected. On the other hand, a gravitational wave background does affect the evolution of magnetic fields. It has also been argued (Tsagas et al. \cite{Tsagas:2001ak},\cite{Tsagas:2005ki}) that this can lead to very strong amplification of a primordial magnetic field. In this paper we revisit these claims and study back reaction to second order.

We consider different renormalizable models of Lorentz invariance violation. We show that the limits on birefringence of the propagation of cosmic microwave background photons from the five year data of the Wilkinson microwave anisotropy probe (WMAP) can be translated into a limit of Lorentz symmetry violation. The obtained limits on Lorentz invariance violation are stronger than other published limits. We also cast them in terms of limits on a birefringent effective photo "mass" and on a polarization dependence of the speed of light.

Explanations of the late-time cosmic acceleration within the framework of general relativity are plagued by difficulties. General relativistic models are mostly based on a dark energy field with fine-tuned, unnatural properties. There is a great variety of models, but all share one feature in common -- an inability to account for the gravitational properties of the vacuum energy, and a failure to solve the so-called coincidence problem. Two broad alternatives to dark energy have emerged as candidate models: these typically address only the coincidence problem and not the vacuum energy problem. The first is based on general relativity and attempts to describe the acceleration as an effect of inhomogeneity in the universe. If this alternative could be shown to work, then it would provide a dramatic resolution of the coincidence problem; however, a convincing demonstration of viability has not yet emerged. The second alternative is based on infra-red modifications to general relativity, leading to a weakening of gravity on the largest scales and thus to acceleration. Most examples investigated so far are scalar-tensor or brane-world models, and we focus on the simplest candidates of each type: $f(R)$ models and DGP models respectively. Both of these provide a new angle on the problem, but they also face serious difficulties. However, investigation of these models does lead to valuable insights into the properties of gravity and structure formation, and it also leads to new strategies for testing the validity of General Relativity itself on cosmological scales.

A primordial gravitational wave background with positive(blue) spectral index is expected in several non-standard inflationary cosmologies where the stress-energy tensor violates the null energy condition. Here we show that a sizable amount of blue gravitational waves is compatible with current cosmological and astrophysical data. So far most of the works on parameter estimation from cosmic microwave background data have assumed a negative or negligible spectral index. The present limits on cosmological parameters, especially on the scalar spectral index, widen up considerably when one allows also for blue tilts of the tensor spectrum. Since the amplitude of the CMB B-mode polarization is larger in these models, future data from Planck are likely to provide crucial measurements.

Observations provide increasingly strong evidence that the universe is accelerating. This revolutionary advance in cosmological observations confronts theoretical cosmology with a tremendous challenge, which it has so far failed to meet. Explanations of cosmic acceleration within the framework of general relativity are plagued by difficulties. General relativistic models are nearly all based on a dark energy field with fine-tuned, unnatural properties. There is a great variety of models, but all share one feature in common -- an inability to account for the gravitational properties of the vacuum energy. Speculative ideas from string theory may hold some promise, but it is fair to say that no convincing model has yet been proposed. An alternative to dark energy is that gravity itself may behave differently from general relativity on the largest scales, in such a way as to produce acceleration. The alternative approach of modified gravity (or dark gravity) provides a new angle on the problem, but also faces serious difficulties, including in all known cases severe fine-tuning and the problem of explaining why the vacuum energy does not gravitate. The lack of an adequate theoretical framework for the late-time acceleration of the universe represents a deep crisis for theory -- but also an exciting challenge for theorists. It seems likely that an entirely new paradigm is required to resolve this crisis.

We consider a two-brane system in a five-dimensional anti-de Sitter spacetime. We study particle creation due to the motion of the physical brane which first approaches the second static brane (contraction) and then recedes from it(expansion). The spectrum and the energy density of the generated gravitons are calculated. We show that the massless gravitons have a blue spectrum and that their energy density satisfies the nucleosynthesis bound with very mild constraints on the parameters. We also show that the Kaluza-Klein modes cannot provide the dark matter in an anti-de-Sitter braneworld. However, for natural choices of parameters, backreaction from the Kaluza-Klein gravitons may well become important. The main findings of this work have been published in the form of a Letter [R. Durrer and M. Ruser, Phys. Rev. Lett. 99, 071601 (2007), arXiv:0704.0756].

In braneworld cosmology the expanding Universe is realized as a brane moving through a warped higher-dimensional spacetime. Like a moving mirror causes the creation of photons out of vacuum fluctuations, a moving brane leads to graviton production. We show that, very generically, Kaluza-Klein (KK) particles scale like stiff matter with the expansion of the Universe and can therefore not represent the dark matter in a warped braneworld. We present results for the production of massless and KK gravitons for bouncing branes in five-dimensional anti de Sitter space. We find that for a realistic bounce the back reaction from the generated gravitons will be most likely relevant. This letter summarizes the main results and conclusions from numerical simulations which are presented in detail in a long paper [M.Ruser and R. Durrer, Phys. Rev. D 76, 104014 (2007), arXiv:0704.0790]

It has been shown that generalized Einstein-Aether theories may lead to significant modifications to the non-relativistic limit of the Einstein equations. In this paper we study the effect of a general class of such theories on the Solar System. We consider corrections to the gravitational potential in negative and positive powers of distance from the source. Using measurements of the perihelion shift of Mercury and time delay of radar signals to Cassini, we place constraints on these corrections. We find that a subclass of generalized Einstein-Aether theories are compatible with these constraints.

We discuss some properties of Lorentz invariant theories which allow for superluminal motion. We show that, if signals are always sent forward in time, closed curves along which signals propagate can be formed. This leads to problems with causality and with the second law of thermodynamics. Only if one singles out one frame with respect to which all signals travel forward in time, the formation of 'closed signal curves' can be prevented. However, the price to pay is that in some reference frames perturbations propagate towards the past or towards the future, depending on the direction of emission.

We investigate microlensing in the case where the lens is considered as an extended object. We use a multipolar expansion of the lens potential and show that the time-varying nature of the quadrupole contribution allows to separate it from the mass and spin contributions and leads to specific modulations of the amplification signal. As example we study the case of binary system lenses in our galaxy. The modulation is observable if the rotation period of the system is smaller than the time over which the amplification is significant and if the impact parameter of the passing light ray is sufficiently close to the Einstein radius so that the amplification is large. Observations of this modulation can reveal important information on the quadrupole and thus on the gravitational radiation emitted by the binary lens. Even if not observed directly, because of their importance the quadrupole modulation has to be included in the error budget for high magnification ($\mu\leq7$) microlensing events.

We show that there are physically relevant situations where gravitational waves do not inherit the frequency spectrum of their source but its wavenumber spectrum.

The power spectrum of a homogeneous and isotropic stochastic variable, characterized by a finite correlation length, does in general not vanish on scales larger than the correlation scale. If the variable is a divergence free vector field, we demonstrate that its power spectrum is blue on large scales. Accounting for this fact, we compute the gravitational waves induced by an incompressible turbulent fluid and by a causal magnetic field present in the early universe. The gravitational wave power spectra show common features: they are both blue on large scales, and peak at the correlation scale. However, the magnetic field can be treated as a coherent source and it is active for a long time. This results in a very effective conversion of magnetic energy in gravitational wave energy at horizon crossing. Turbulence instead acts as a source for gravitational waves over a time interval much shorter than a Hubble time, and the conversion into gravitational wave energy is much less effective. We also derive a strong constraint on the amplitude of a primordial magnetic field when the correlation length is much smaller than the horizon.

I describe the imprint of primordial magnetic fields on the CMB. I show that these are observable only if the field amplitude is of the order of $B\gsim 10^{-9}G$ on Mpc scale. I further argue that such fields are strongly constrained by the stochastic background of gravity waves which they produce. Primordial magnetic fields, which are strong enough to be seen in the CMB, are compatible with the nucleosynthesis bound, only if their spectrum is close to scale invariant, or maybe if helical magnetic fields provoke an inverse cascade. For helical fields, the CMB signature is especially interesting. It contains parity violating T-B and E--B correlations.

We demonstrate that if k-essence can solve the coincidence problem and play the role of dark energy in the universe, the fluctuations of the field have to propagate superluminally at some stage. We argue that this implies that successful k-essence models violate causality. It is not possible to define a time ordered succession of events in a Lorentz invariant way. Therefore, k-essence cannot arise as low energy effective field theory of a causal, consistent high energy theory.

We show that the dipole of the luminosity distance is a useful observational tool which allows us to determine the Hubble parameter as a function of the redshift H(z). We determine the number of supernovae needed to achieve a given precision for H(z) and to distinguish between different models for dark energy. We analyse a sample of nearby supernovae and find a dipole consistent with the cosmic microwave background at a significance of more than 2 sigma.

We consider super-gravity models in which the lightest supersymmetric particle (LSP) is a stable gravitino. The next-to-lightest supersymmetric particle (NLSP) freezes out with its thermal relic density and then decays after $(10^5-10^{10})$ sec, injecting high-energy photons into the cosmic plasma. These photons heat up the electron plasma which then thermalizes with the cosmic microwave background (CMB) via Compton scattering, bremsstrahlung and double-Compton scattering. Contrary to previous studies which assume instantaneous energy injection, we solve the full kinetic equation for the photon number density with a source term describing the decay of the NLSP. This source term is based on the requirement that the injected energy be almost instantaneously redistributed by Compton scattering, hence leading to a time-dependent chemical potential. We investigate the case of a stau NLSP and determine the constraints on the gravitino and stau masses from observations of the CMB spectrum by assuming that all gravitino LSPs come from stau NLSP decays. Unlike the analytical approximations, we find that there may be a stau mass below which the constraint from the CMB spectrum vanishes.

We study braneworlds in a five dimensional bulk, where cosmological expansion is mimicked by motion through AdS$_5$. We show that the five dimensional graviton reduces to the four dimensional one in the late time approximation of such braneworlds. Inserting a fixed regulator brane far from the physical brane, we investigate quantum graviton production due to the motion of the brane. We show that the massive Kaluza-Klein modes decouple completely from the massless mode and they are not generated at all in the limit where the regulator brane position goes to infinity. In the low energy limit, the massless four dimensional graviton obeys the usual 4d equation and is therefore also not generated in a radiation-dominated universe.

This course is an introduction to the physics of braneworlds. We concentrate on braneworlds with only one extra-dimension and discuss their gravity. We derive the gravitational equations on the brane from the bulk Einstein equation and explore some limits in which they reduce to 4-dimensional Einstein gravity. We indicate how cosmological perturbations from braneworlds are probably very different from usual cosmological perturbations and give some examples of the preliminary results in this active field of research. For completeness, we also present an introduction to 4-dimensional cosmological perturbation theory and, especially its application to the anisotropies of the cosmic microwave background.

We show that scalar as well as vector and tensor metric perturbations in the Randall-Sundrum II braneworld allow normalizable tachyonic modes, i.e., possible instabilities. These instabilities require nonvanishing initial anisotropic stresses on the brane. We show with a specific example that within the Randall-Sundrum II model, even though the tachyonic modes are excited, no instability develops. We argue, however, that in the cosmological context instabilities might in principle be present. We conjecture that the tachyonic modes are due to the singularity of the orbifold construction. We illustrate this with a simple but explicit toy model.

In their recent paper ``Faraday rotation of the cosmic microwave background polarization by a stochastic magnetic field'', Kosowsky et al. Phys.Rev. D71, 043006 (2005) have commented about our paper [C.Caprini and R.Durrer, Phys. Rev. D65, 023517 (2002)], in which we derived very strong limits on the amplitude of a primordial magnetic field from gravitational wave production. They argue that our limits are erroneous. In this short comment we defend our result.

In the context of string theory we argue that higher dimensional Dp-branes unwind and evaporate so that we are left with D3-branes embedded in a (9+1)-dimensional bulk. One of these D3-branes plays the role of our Universe. Within this picture, the evaporation of the higher dimensional Dp-branes provides the entropy of our Universe.

We study the effect of a possible helicity component of a primordial magnetic field on the tensor part of the cosmic microwave background temperature anisotropies and polarization. We give analytical approximations for the tensor contributions induced by helicity, discussing their amplitude and spectral index in dependence of the power spectrum of the primordial magnetic field. We find that an helical magnetic field creates a parity odd component of gravity waves inducing parity odd polarization signals. However, only if the magnetic field is close to scale invariant and if its helical part is close to maximal, the effect is sufficiently large to be observable. We also discuss the implications of causality on the magnetic field spectrum.

This is a review on cosmological perturbation theory. After an introduction, it presents the problem of gauge transformation. Gauge invariant variables are introduced and the Einstein and conservation equations are written in terms of these variables. Some examples, especially perfect fluids and scalar fields are presented in detail. The generation of perturbations during inflation is studied. Lightlike geodesics and their relevance for CMB anisotropies are briefly discussed. Perturbation theory in braneworlds is also introduced.

In this letter we calculate the emission of gravity waves by the binary pulsar in the framework of five dimensional braneworlds. We consider only spacetimes with one compact extra-dimension. We show that the presence of additional degrees of freedom, especially the 'gravi-scalar' leads to a modification of Einstein's quadrupole formula. We compute the induced change for the binary pulsar PSR 1913+16 where it amounts to about 20% which is by far excluded by present experimental data.

We introduce the concepts of adiabatic (curvature) and isocurvature (entropy) cosmological perturbations and present their relevance for parameter estimation from cosmic microwave background anisotropies data. We emphasize that, while present-day data are in excellent agreement with pure adiabaticity, subdominant isocurvature contributions cannot be ruled out. We discuss model independent constraints on the isocurvature contribution. Finally, we argue that the Planck satellite will be able to do precision cosmology even if the assumption of adiabaticity is relaxed.

We discuss the implications of causality on a primordial magnetic field. We show that the residual field on large scales is much more suppressed than usually assumed and that a helical component is even more reduced. Due to this strong suppression, even maximal primordial fields generated at the electroweak phase transition can just marginally seed the fields in clusters, but they cannot leave any detectable imprint on the cosmic microwave background.

The vector perturbations induced on the brane by gravitational waves propagating in the bulk are studied in a cosmological framework. Cosmic expansion arises from the brane motion in a non-compact five-dimensional anti-de Sitter spacetime. By solving the vector perturbation equations in the bulk, for generic initial conditions, we find that they give rise to growing modes on the brane in the Friedmann-Lema\^{\i}tre era. Among these modes, we exhibit a class of normalizable perturbations, which are exponentially growing with respect to conformal time on the brane. The presence of these modes is, at least, strongly constrained by the current observations of the cosmic microwave background (CMB). We estimate the anisotropies they induce in the CMB, and derive quantitative constraints on the allowed amplitude of their primordial spectrum. Our results provide stringent constraints for all braneworld models with bulk inflation.

Threshold biasing of a Gaussian random field gives a linear amplification of the reduced two point correlation function at large distances. We show that for standard cosmological models this does not translate into a linear amplification of the power spectrum (PS) at small k. For standard CDM type models this means that the``turn-over'' at small k of the original PS disappears in the PS of the biased field for the physically relevant range of the threshold parameter. In real space this difference is manifest in the asymptotic behaviour of the normalised mass variance in spheres of radius R, which changes from the ``super-homogeneous'' behaviour to a Poisson-like behaviour.This qualitative change results from the intrinsic stochasticity of the threshold sampling. While our quantitative results are specific to the simplest threshold biasing model, we argue that our qualitative conclusions should be valid generically for any biasing mechanism involving a scale-dependent amplification of the correlation function. One implication is that the real-space correlation function will be a better instrument to probe for the underlying Harrison Zeldovich spectrum in the distribution of visible matter, as the characteristic asymptotic negative power-law \xi (r) \sim -r^{-4} tail is undistorted by biasing.

We investigate in detail the question whether a non-vanishing cosmological constant is required by present-day cosmic microwave background and large scale structure data when general isocurvature initial conditions are allowed for. We also discuss differences between the usual Bayesian and the frequentist approaches in data analysis. We show that the COBE-normalized matter power spectrum is dominated by the adiabatic mode and therefore breaks the degeneracy between initial conditions which is present in the cosmic microwave background anisotropies. We find that in a flat universe the Bayesian analysis requires \Omega_\Lambda \neq 0 to more than 3 \sigma, while in the frequentist approach \Omega_\Lambda = 0 is still within 3 \sigma for a value of h < 0.48. Both conclusions hold regardless of initial conditions.

The locations and amplitudes of three acoustic peaks and two dips in the last releases of the Boomerang, MAXIMA and DASI measurements of the cosmic microwave background (CMB) anisotropy power spectra as well as their statistical confidence levels are determined in a model-independent way. It is shown that the Boomerang-2001 data (Netterfield et al. 2001) fixes the location and amplitude of the first acoustic peak at more than 3\sigma confidence level. The next two peaks and dips are determined at a confidence level above 1\sigma but below 2\sigma. The locations and amplitudes of the first three peaks and two dips are 212+/-17, 5426+/-1218\mu K^2, 544+/-56, 2266+/-607\mu K^2, 843+/-35, 2077+/-876\mu K^2, 413+/-50, 1960+/-503\mu K^2, 746+/-89, 1605+/-650\mu K^2 respectively (1\sigma errors include statistical and systematic errors). The MAXIMA and DASI experiments give similar values for the extrema which they determine. The determined cosmological parameters from the CMB acoustic extrema data show good agreement with other determinations, especially with the baryon content as deduced from standard nucleosynthesis constraints. These data supplemented by the constraints from direct measurements of some cosmological parameters and data on large scale structure lead to a best-fit model which agrees with practically all the used experimental data within 1\sigma. The best-fit parameters are: \Omega_{\Lambda}=0.64^{+0.14}_{-0.27}, \Omega_{m}= 0.36^{+0.21}_{-0.11}, \Omega_b=0.047^{+0.093}_{-0.024}, n_s=1.0^{+0.59}_{-0.17}, h=0.65^{+0.35}_{-0.27} and \tau_c=0.15^{+0.95}_{-0.15} (plus/minus values show 1\sigma upper/lower limits obtained by marginalization over all other model parameters). The best-fit values of \Omega_{\nu} and T/S are close to zero, their 1\sigma upper limits are 0.17 and 1.7 respectively.

We present some phenomenological aspects of the pre-big-bang cosmological model inspired by the duality properties of string theory. In particular, assuming the spatial sections of the homogeneous background geometry to be isotropic, we discuss the quantum production of perturbations of the background fields (gravitons, dilatons, moduli fields), as well as the production of particles which do not contribute to the background, which we call ``seeds''. As such we consider the cases of electromagnetic and axionic seeds. We also discuss their possible observational consequences, for example, we study whether they can provide the origin of primordial galactic magnetic fields, and whether they can generate the initial fluctuations leading to the formation of large-scale structure and the measured cosmic microwave background anisotropies. We finally analyze axion and photon production in four dimensional anisotropic pre-big-bang cosmological models.

At low energy, the four-dimensional effective action of the ekpyrotic model of the universe is equivalent to a slightly modified version of the pre big bang model. We discuss cosmological perturbations in these models. In particular we address the issue of matching the perturbations from a collapsing to an expanding phase in full generality. We show that, generically, one obtains $n=0$ for the spectrum of scalar perturbations in the original pre big model (with vanishing potential). When an exponential potential for the dilaton is included, a scale invariant spectrum ($n=1$) of adiabatic scalar perturbations is produced under very generic matching conditions, both in a modified pre big bang and ekpyrotic scenario. We also derive general results valid for power law scale factors matched to a radiation dominated era.

We derive the gauge invariant perturbation equations for a 5-dimensional bulk spacetime in the presence of a brane. The equations are derived in full generality, without specifying a particular energy content of the bulk or the brane. We do not assume Z_2 symmetry, and show that the degree of freedom associated with brane motion plays a crucial role. Our formalism may also be used in the Z_2 symmetric case where it simplifies considerably.

In this talk I summarize the topics addressed during the conference. We have discussed the most relevant cosmological observations of the last say 10 years and their implications for our understanding of the Universe. My finding throughout this meeting has been that it is amazing how our knowledge of cosmological parameters has improved during the last years. It comes as a surprise to many, that the simplest inflationary model of adiabatic perturbations agrees well with present data. But it is frustrating that we have not made any progress in identifying the dark matter, in contrary, we have supplemented it by what we call the 'dark energy', a gravitationally and esthetically repulsive component which is even more mysterious than the dark matter. We have no understanding about the origin of these two most abundant energy components of the present Universe.

In this note I try to clarify the problem of perturbations in the ekpyrotic universe. I write down the most general matching conditions and specify the choices taken by the two debating sides. I also bring up the problem of surface stresses which always have to be present when a transition from a collapsing to an expanding phase is made.

Topological defects are ubiquitous in physics. Whenever a symmetry breaking phase transition occurs, topological defects may form. The best known examples are vortex lines in type II super conductors or in liquid Helium, and declination lines in liquid crystals. In an adiabatically expanding universe which cools down from a very hot initial state, it is quite natural to postulate that topological defects may have emerged during a phase transition in the early universe and that they may have played the role of initial inhomogeneities seeding the formation of cosmic structure. This basic idea goes back to Kibble (1976). In this report we summarize the progress made in the investigation of Kibble's idea during the last 25 years. Our understanding of the formation and evolution of topological defects is reported almost completely in the beautiful book by Vilenkin & Shellard or the excellent Review by Hindmarsh & Kibble, and we shall hence be rather short on that topic. Nevertheless, in order to be self contained, we have included a short chapter on spontaneous symmetry breaking and defect formation. Our main topic is however the calculation of structure formation with defects, results which are not included in the above references.

We compute the gravity waves induced by anisotropic stresses of stochastic primordial magnetic fields. The nucleosynthesis bound on gravity waves is then used to derive a limit on the magnetic field amplitude as function of the spectral index. The obtained limits are extraordinarily strong: If the primordial magnetic field is produced by a causal process, leading to a spectral index $n\ge 2$ on super horizon scales, galactic magnetic fields produced at the electroweak phase transition or earlier have to be weaker than $B_\la \le 10^{-27}$Gauss! If they are induced during an inflationary phase (reheating temperature $T\sim 10^{15}$GeV) with a spectral index $n\sim 0$, the magnetic field has to be weaker than $B_\la \le 10^{-39}$Gauss! Only very red magnetic field spectra, $n\sim -3$ are not strongly constrained. We also find that a considerable amount of the magnetic field energy is converted into gravity waves. The gravity wave limit derived in this work rules out most of the proposed processes for primordial seeds for the large scale coherent magnetic fields observed in galaxies and clusters.

This is a review of the theory of CMB anisotropies, an updated version of a course given at the troisieme cycle de la Suisse Romande. An introduction to gauge invariant cosmological perturbation theory is given and the theory CMB anisotropies is develiped in this context. Simple analytical approximations for the acoustic peak positions for adiabatic and isocurvature perturbations are derived. Silk damping is discussed by an analytic approximation. A short description of the present status of observations and parameter estimation followed by a critical discussion terminate the review. The full system of differential equations for CMB anisotropies and polarization needed in a numerical treatment is also developed and given in an appendix.

The physics of the origin and evolution of CMB anisotropies is described. I explain the idea and status of cosmic parameter estimation and follow it up with critical comments on its dependence on model assumptions and initial conditions.

Recently high quality data of the cosmic microwave background anisotropies have been published. In this work we study to which extent the cosmological parameters determined by using this data depend on assumptions about the initial conditions. We show that for generic initial conditions, not only the best fit values are very different but, and this is our main result, the allowed parameter range enlarges dramatically.

We derive dynamical equations to describe a single 3-brane containing fluid matter and a scalar field coupling to the dilaton and the gravitational field in a five dimensional bulk. First, we show that a scalar field or an arbitrary fluid on the brane cannot evolve to cancel the cosmological constant in the bulk. Then we show that the Randall-Sundrum model is unstable under small deviations from the fine-tuning between the brane tension and the bulk cosmological constant and even under homogeneous gravitational perturbations. Implications for brane world cosmologies are discussed.

We present an alternative scenario for cosmic structure formation where initial fluctuations are due to Kalb-Ramond axions produced during a pre-big bang phase of inflation. We investigate whether this scenario, where the fluctuations are induced by seeds and therefore are of isocurvature nature, can be brought in agreement with present observations by a suitable choice of cosmological parameters. We also discuss several observational signatures which can distinguish axion seeds from standard inflationary models. We finally discuss the gravitational wave background induced in this model and we show that it may be well within the range of future observations.

We use observational data on the large scale structure (LSS) of the Universe measured over a wide range of scales from sub-galactic up to horizon scale and on the cosmic microwave background anisotropies to determine cosmological parameters within the class of adiabatic inflationary models. We show that a mixed dark matter model with cosmological constant ($\Lambda$MDM model) and parameters $\Omega_m=0.37^{+0.25}_{-0.15}$, $\Omega_{\Lambda}=0.69^{+0.15}_{-0.20}$, $\Omega_{\nu}=0.03^{+0.07}_{-0.03}$, $N_{\nu}=1$, $\Omega_b=0.037^{+0.033}_{-0.018}$, $n_s=1.02^{+0.09}_{-0.10}$, $h=0.71^{+0.22}_{-0.19}$, $b_{cl}=2.4^{+0.7}_{-0.7}$ (1$\sigma$ confidence limits) matches observational data on LSS, the nucleosynthesis constraint, direct measurements of Hubble constant, the high redshift supernova type Ia results and the recent measurements of the location and amplitude of the first acoustic peak in the CMB anisotropy power spectrum. The best model is $\Lambda$ dominated (65% of the total energy density) and has slightly positive curvature, $\Omega=1.06$. The clustered matter consists in 8% massive neutrinos, 10% baryons and 82% cold dark matter (CDM). The upper 2$\sigma$ limit on the neutrino content can be expressed in the form $\Omega_{\nu}h^2/N_{\nu}^{0.64}\le0.042$ or, via the neutrino mass, $m_{\nu}\le4.0$eV. The upper 1(2)$\sigma$ limit for the contribution of a tensor mode to the COBE DMR data is T/S$<1(1.5)$. Furthermore, it is shown that the LSS observations together with the Boomerang (+MAXIMA-1) data on the first acoustic peak rule out zero-$\Lambda$ models at more than $2\sigma$ confidence limit.

During the last years it has become clear that global O(N) defects and U(1) cosmic strings do not lead to the pronounced first acoustic peak in the power spectrum of anisotropies of the cosmic microwave background which has recently been observed to high accuracy. Inflationary models cannot easily accommodate the low second peak indicated by the data. Here we construct causal scaling seed models which reproduce the first and second peak. Future, more precise CMB anisotropy and polarization experiments will however be able to distinguish them from the ordinary adiabatic models.

We present an alternative scenario for cosmic structure formation where initial fluctuations are due to Kalb-Ramond axions produced during a pre-big bang phase of inflation. We investigate whether this scenario, where the fluctuations are induced by seeds and therefore are of isocurvature nature, can be brought in agreement with present observations by a suitable choice of cosmological parameters. We also discuss several observational signatures which can distinguish axion seeds from standard inflationary models.

Observational data on the large scale structure (LSS) of the Universe are used to establish an upper limit for the amplitude of the tensor mode marginalized over all other cosmological parameters within the class of adiabatic inflationary models. It is shown that the upper 1$\sigma$ limit for the contribution of a tensor mode to the COBE DMR data is T/S$<1$.

A primordial stochastic magnetic field will induce temperature and polarization fluctuations in the cosmic microwave background. We outline a calculation of the resulting fluctuation power spectra and present numerical results.

The observational data on the large scale structure (LSS) of the Universe are used to establish the upper limit on the neutrino content marginalized over all other cosmological parameters within the class of adiabatic inflationary models. It is shown that the upper 2$\sigma$ limit on the neutrino content can be expressed in the form $\Omega_{\nu}h^2/N_{\nu}^{0.64}\le0.042$ or, via the neutrino mass, $m_{\nu}\le4.0$eV.

The observational data on the large scale structure (LSS) of the Universe are used to determine cosmological parameters within the class of adiabatic inflationary models. We show that a mixed dark matter model with cosmological constant ($\Lambda$MDM model) and parameters $\Omega_m=0.37^{+0.25}_{-0.15}$, $\Omega_{\Lambda}=0.69^{+0.15}_{-0.20}$, $\Omega_{\nu}=0.03^{+0.07}_{-0.03}$, $N_{\nu}=1$, $\Omega_b=0.037^{+0.033}_{-0.018}$, $n_s=1.02^{+0.09}_{-0.10}$, $h=0.71^{+0.22}_{-0.19}$, $b_{cl}=2.4^{+0.7}_{-0.7}$ (1$\sigma$ confidence limits) matches observational data on LSS, the nucleosynthesis constraint, direct measurements of the Hubble constant, the high redshift supernova type Ia results and the recent measurements of the location and amplitude of the first acoustic peak in the CMB anisotropy power spectrum. The best model is $\Lambda$ dominated (65% of the total energy density) and has slightly positive curvature, $\Omega=1.06$. The clustered matter consists in 8% massive neutrinos, 10% baryons and 82% cold dark matter (CDM). It is shown that the LSS observations together with the Boomerang (+MAXIMA-1) data on the first acoustic peak rule out zero-$\Lambda$ models at more than $2\sigma$ confidence level.

We compute the skewness of the matter distribution arising from non-linear evolution and from non-Gaussian initial perturbations. We apply our result to a very generic class of models with non-Gaussian initial conditions and we estimate analytically the ratio between the skewness due to non-linear clustering and the part due to the intrinsic non-Gaussianity of the models. We finally extend our estimates to higher moments.

Structure formation with topological defects is described. The main differences from inflationary models are highlighted. The results are compared with recent observations. It is concluded that all the defect models studied so far are in disagreement with recent observations of CMB anisotropies. Furthermore, present observations do not support 'decoherence', a generic feature of structure formation from topological defects.

We compute the energy spectra for massless Kalb-Ramond axions in four-dimensional anisotropic string cosmological models. We show that, when integrated over directions, the four-dimensional anisotropic model leads to infra-red divergent spectra similar to the one found in the isotropic case.

The possibility of determining cosmological parameters on the basis of a wide set of observational data including the Abell-ACO cluster power spectrum and mass function, peculiar velocities of galaxies, the distribution of Ly-$\alpha$ clouds and CMB temperature fluctuations is analyzed. Using a $\chi^2$ minimization method, assuming $\Omega_{\Lambda}+\Omega_m =1$ and no contribution from gravity waves, we found that a tilted $\Lambda$MDM model with one sort of massive neutrinos and the parameters $n\approx 1.12$, $\Omega_m\approx 0.4$ ($\Omega_{\Lambda}\approx 0.6$), $\Omega_{cdm}\approx 0.3$, $\Omega_{\nu}\approx 0.06$, $\Omega_b\approx 0.04$ and $h\approx 0.7$ matches observational data best. $\Omega_{\nu}$ is higher for more species of massive neutrinos, $\sim 0.1$ for two and $\sim 0.13$ for three species. $\Omega_m$ raises by $\sim 0.08$ and $\sim 0.15$ respectively. Varying only a subset of parameters and fixing the others shows also that the observational data set used here rules out pure CDM models with $h\ge 0.5$, scale invariant primordial power spectrum, zero cosmological constant and spatial curvature at a very high confidence level, $>99.99%$. The corresponding class of MDM models are ruled out at $\sim 95%$ C.L. It is notable also that this data set determines the amplitude of scalar fluctuations approximately at the same level as COBE four-year data. It indicates that a possible tensor component in the COBE data cannot be very substantial.

The possibility of determining cosmological parameters on the basis of a wide set of observational data including the Abell-ACO cluster power spectrum and mass function, peculiar velocities of galaxies, the distribution of Ly-$\alpha$ clouds and CMB temperature fluctuations is analyzed. Using a $\chi^2$ minimization method, assuming $\Omega_{\Lambda}+\Omega_{\rm{matter}} =1 $ and no contribution from gravity waves, we found that a tilted $\Lambda$MDM model with one sort of massive neutrinos and the parameters $n=1.12\pm 0.10$, $\Omega_m=0.41\pm 0.11$ ($\Omega_{\Lambda}=0.59\pm0.11$), $\Omega_{cdm}=0.31\pm 0.15$, $\Omega_{\nu}=0.059\pm 0.028$, $\Omega_b=0.039\pm 0.014$ and $h=0.70\pm 0.12$ matches observational data best. The 1$\sigma$ (68.3%) confidence limits on each cosmological parameter, which are obtained by marginalizing over the other parameters, are $0.82\le n\le1.39$, $0.19\le\Omega_m\le 1$ ($0\le\Omega_{\Lambda}\le 0.81$), $0\le\Omega_{\nu}\le 0.17$, $0.021\le \Omega_b\le 0.13$ and $0.38\le h\le 0.85$ $1.5\le b_{cl}\le 3.5$. Here $b_{cl}$ is the cluster bias parameter. The best-fit parameters for 31 models which are inside of $1\sigma$ range of the best model are presented. It is shown also that observational data set used here rules out the class of CDM models with $h\ge 0.5$, scale invariant primordial power spectrum, zero cosmological constant and spatial curvature at very high confidence level, $>99.99%$. The corresponding class of MDM models are ruled out at $\sim 95%$ C.L.

In this letter we show that in a Gaussian random field the correlation length, the typical size of correlated structures, does not change with biasing. We interpret the amplification of the correlation functions of subsets identified by different thresholds being due to the increasing sparseness of peaks over threshold. This clarifies an long-standing misconception in the literature. We also argue that this effect does not explain the observed increase of the amplitude of the correlation function xi(r) when galaxies of brighter luminosity or galaxy clusters of increasing richness are considered.

In this letter we investigate whether the isotropy problem is naturally solved in inflationary cosmologies inspired by string theory, so called pre-big-bang cosmologies. We find that, in contrast to what happens in the more common 'potential inflation' models, initial anisotropies do not decay during pre-big-bang inflation.

A recently proposed mechanism for large-scale structure in string cosmology --based on massless axionic seeds-- is further analyzed and extended to the acoustic-peak region. Existence, structure, and normalization of the peaks turn out to depend crucially on the overall evolution of extra dimensions during the pre-big bang phase: conversely, precise cosmic microwave background anisotropy data in the acoustic-peak region will provide, within the next decade, a window on string-theory's extra dimensions before their eventual compactification.

We derive an expression for the angular power spectrum of cosmic microwave background anisotropies due to gravity waves generated by a stochastic magnetic field and compare the result with current observations; we take into account the non-linear nature of the stress energy tensor of the magnetic field. For almost scale invariant spectra, the amplitude of the magnetic field at galactic scales is constrained to be of order 10^{-9} Gauss. If we assume that the magnetic field is damped below the Alfven damping scale, we find that its amplitude at 0.1 h^{-1}Mpc, B_\lambda, is constrained to be B_\lambda<7.9 x10^{-6} e^{3n} Gauss, for n<-3/2, and B_\lambda<9.5x10^{-8} e^{0.37n} Gauss, for n>-3/2, where n is the spectral index of the magnetic field and H_0=100h km s^{-1}Mpc^{-1} is the Hubble constant today.

An accurate analytic approximation of the transfer function for the power spectra of primordial density perturbations in mixed dark matter models is presented. The fitting formula in a matter-dominated Universe ($\Omega_0=\Omega_M=1$) is a function of wavenumber $k$, redshift $z$ and four cosmological parameters: the density of massive neutrinos, $\Omega_{\nu}$, the number of massive neutrino species, $N_{\nu}$, the baryon density, $\Omega_{b}$ and the dimensionless Hubble constant, $h$. Our formula is accurate in a broad range of parameters: $k\le 100 h/Mpc$, $z\le 30$, $\Omega_{\nu}\le 0.5$, $N_{\nu}\le 3$, $\Omega_{b}\le 0.3$, $0.3\le h\le 0.7$. The deviation of the variance of density fluctuations calculated with our formula from numerical results obtained with CMBfast is less than 6% for the entire range of parameters. It increases with $\Omega_bh^{2}$ and is less than $\le 3%$ for $\Omega_bh^{2}\le 0.05$. The performance of the analytic approximation of MDM power spectra proposed here is compared with other approximations found in the literature. Our approximation turns out to be closest to numerical results in the parameter space considered here.

We investigate the global texture model of structure formation in cosmogonies with non-zero cosmological constant for different values of the Hubble parameter. We find that the absence of significant acoustic peaks and little power on large scales are robust predictions of these models. However, from a careful comparison with experiments we conclude that at present we cannot safely reject the model on the grounds of published CMB anisotropy data. If bias is close to one on large scales, galaxy correlation data rules out the models. New, very stringent constraints come from peculiar velocities. Investigating the large-N limit, we argue that our main conclusions apply to all global O(N) models of structure formation.

We consider harmonic maps from Minkowski space into the three sphere. We are especially interested in solutions which are asymptotically constant, i.e. converge to the same value in all directions of spatial infinity. Physical 3-space can then be compactified and can be identified topologically (but not metrically!) with a three sphere. Therefore, at fixed time, the winding of the map is defined. We investigate whether static solutions with non-trivial winding number exist. The answer which we can proof here is only partial: We show that within a certain family of maps no static solutions with non-zero winding number exist. We discuss the existing static solutions in our family of maps. An extension to other maps or a proof that our family of maps is sufficiently general remains an open problem.

We derive a simple closed-form expression, relating $\vs(r)$ -- the mean relative velocity of pairs of galaxies at fixed separation $r$ -- to the two-point correlation function of mass density fluctuations, $\xi(r)$. We compare our analytic model for $\vs(r)$ with N-body simulations, and find excellent agreement in the entire dynamical range probed by the simulations ($0.1 \lsim \xi \lsim 1000$). Our results can be used to estimate the cosmological density parameter, $\Om$, directly from redshift-distance surveys, like Mark III.

We investigate the global texture model of structure formation in cosmogonies with non-zero cosmological constant for different values of the Hubble parameter. We find that the absence of significant acoustic peaks and little power on large scales are robust predictions of these models. However, from a careful comparison with data we conclude that at present we cannot safely reject the model on the grounds of present CMB data. Exclusion by means of galaxy correlation data requires assumptions on biasing and statistics. New, very stringent constraints come from peculiar velocities. Investigating the large-N limit, we argue that our main conclusions apply to all global O(N) models of structure formation.

In this letter we present an idea which reconciles a homogeneous and isotropic Friedmann universe with a fractal distribution of galaxies. We use two observational facts: The flat rotation curves of galaxies and the (still debated) fractal distribution of galaxies with fractal dimension D=2. Our idea can also be interpreted as a redefinition of the notion of bias.

We investigate microwave background anisotropies in the presence of primordial magnetic fields. We show that a homogeneous field with fixed direction can amplify vector perturbations. We calculate the correlations of $\delta T/T$ explicitly and show that a large scale coherent field induces correlations between $a_{\ell-1,m}$ and $a_{\ell+1,m}$. We discuss constraints on amplitude and spectrum of a primordial magnetic field imposed by observations of CMB anisotropies.

A primordial stochastic background of very weakly coupled massless (pseudo-)scalars can seed CMB anisotropy, when large-scale fluctuations of their stress-tensor re-enter the horizon during the matter-dominated era. A general relation between multipole coefficients of the CMB anisotropy and the seed's energy spectrum is derived. Magnitude and tilt of the observed anisotropies can be reproduced for the nearly scale-invariant axion spectra that are predicted in a particularly symmetric class of string cosmology backgrounds.

Anisotropies in the cosmic microwave background radiation due to gravity waves are investigated. An initial spectrum of gravity waves may have been induced during an epoch of inflation. We study the propagation of such a spectrum in a Friedmann background of hot and cold dark matter, radiation and (possibly) a cosmological constant. We finally calculate its imprint as anisotropies on the cosmic microwave background. We also take into account that massless particles can source gravity waves by their anisotropic stresses. We consider general mixed dark matter models with and without cosmological constant. For a given, scale invariant input spectrum of gravity we determine the dependence of the resulting spectrum of CMB anisotropies on the different model parameters.

Pre-big bang cosmology predicts tiny first-order dilaton and metric perturbations at very large scales. Here we discuss the possibility that other -- more copiously generated -- perturbations may act, at second order, as scalar seeds of large-scale structure and CMB anisotropies. We study, in particular, the cases of electromagnetic and axionic seeds. We compute the stochastic fluctuations of their energy-momentum tensor and determine the resulting contributions to the multipole expansion of the temperature anisotropy. In the axion case it is possible to obtain a flat or slightly tilted blue spectrum that fits present data consistently, both for massless and for massive (but very light) axions.

We compute cosmic microwave background angular power spectra for scaling seed models of structure formation. A generic parameterization of the energy momentum tensor of the seeds is employed. We concentrate on two regions of parameter space inspired by global topological defects: O(4) texture models and the large-N limit of O(N) models. We use $\chi^{2}$ fitting to compare these models to recent flat-band power measurements of the cosmic microwave background. Only scalar perturbations are considered.

In this work we present a partially new method to analyze fluctuations which are induced by causal scaling seeds. We show that the power spectra due to this kind of seed perturbations are determined by five analytic functions, which we determine numerically for a special example. We put forward the view that, even if recent work disfavors the models with cosmic strings and global O(4) texture, causal scaling seed perturbations merit a more thorough and general analysis, which we initiate in this paper.

We discuss various questions which arise when one considers the central projection of three dimensional fractal sets (galaxy catalogs) onto the celestial globe. The issues are related to how fractal such projections look. First we show that the lacunarity in the projection can be arbitrarily small. Further characteristics of the projected set---in particular scaling---depend sensitively on how the apparent sizes of galaxies are taken into account. Finally, we discuss the influence of opacity of galaxies. Combining these ideas, seemingly contradictory statements about lacunarity and apparent projections can be reconciled.

We study microwave background anisotropies induced by scaling seed perturbations in a universe dominated by cold dark matter. Using a gauge invariant linear perturbation analysis, we solve the perturbation equations on super-horizon scales, for CMB anisotropies triggered by generic gravitational seeds. We find that perturbations induced by seeds --- under very mild restrictions --- are nearly isocurvature. Thus, compensation, which is mainly the consequence of physically sensible initial conditions, is very generic. We then restrict our study to the case of scaling sources, motivated by global scalar fields. We parameterize the energy momentum tensor of the source by ``seed functions'' and calculate the Sachs-Wolfe and acoustic contributions to the CMB anisotropies. We discuss the dependence of the anisotropy spectrum on the parameters of the model considered. Even within the restricted class of models investigated in this work, we find a surprising variety of results for the position and height of the first acoustic peak as well as for the overall amplitude. In particular, for certain choices of parameters, the spectrum resembles very much the well known adiabatic inflationary spectrum, whereas for others, the position of the first acoustic peak is significantly shifted towards smaller angular scales.

An introduction to topological defects in cosmology is given. We discuss their possible relevance for structure formation. Especial emphasis is given on the signature of topological defects in the spectrum of anisotropies in the cosmic microwave background. We present simple analytic estimates for the CMB spectrum on large and intermediate scales and compare them with the corresponding approximations for models where initial perturbations are generated during an inflationary epoch.

We present an analysis of CMB anisotropies induced by global scalar fields in the large N limit. In this limit, the CMB anisotropy spectrum can be determined without cumbersome 3D simulations. We determine the source functions and their unequal time correlation functions and show that they are quite similar to the corresponding functions in the texture model. This leads us to the conclusion that the large N limit provides a 'cheap approximation' to the texture model of structure formation.

The analysis of anisotropies in the cosmic microwave background (CMB) has become an extremely valuable tool for cosmology. We even have hopes that planned CMB anisotropy experiments may revolutionize cosmology. Together with determinations of the CMB spectrum, they represent the first cosmological precision measurements. This is illustrated in the talk by Anthony Lasenby. The value of CMB anisotropies lies to a big part in the simplicity of the theoretical analysis. Fluctuations in the CMB can be determined almost fully within linear cosmological perturbations theory and are not severely influenced by complicated nonlinear physics. In this contribution the different physical processes causing or influencing anisotropies in the CMB are discussed. The geometry perturbations at and after last scattering, the acoustic oscillations in the baryon-photon-plasma prior to recombination, and the diffusion damping during the process of recombination. The perturbations due to the fluctuating gravitational field, the so called Sachs-Wolfe contribution, is described in a very general form using the Weyl tensor of the perturbed geometry.

We review recent work aimed at showing how global topological defects influence the shape of the angular power spectrum of the CMB radiation on small scales. While Sachs-Wolfe fluctuations give the dominant contribution on angular scales larger than about a few degrees, on intermediate scales the main r\^ole is played by coherent oscillations in the baryon radiation plasma before recombination. In standard cosmological models these oscillations lead to the `Doppler peaks' in the angular power spectrum. Inflation-based cold dark matter models predict the location of the first peak to be at $\ell\sim 220/\sqrt{\Omega_0}$, with a height which is a few times the level of anisotropies at large scales. Here we focus on perturbations induced by global textures. We find that the height of the first peak is reduced and is shifted to $\ell\sim 350$.

The aim of this paper is to show, that the 'oscillating universe' is a viable alternative to inflation. We remind that this model provides a natural solution to the flatness or entropy and to the horizon problem of standard cosmology. We study the evolution of density perturbations and determine the power spectrum in a closed universe. The results lead to constraints of how a previous cycle might have looked like. We argue that most of the radiation entropy of the present universe may have originated from gravitational entropy produced in a previous cycle. We show that measurements of the power spectrum on very large scales could in principle decide whether our universe is closed, flat or open.

The Doppler peaks (Sacharov peaks) in the angular power spectrum of the cosmic microwave background anisotropies, are mainly due to coherent oscillations in the baryon radiation plasma before recombination. Here we present a calculation of the Doppler peaks for perturbations induced by global textures and cold dark matter. We find that the height of the first Doppler peak is smaller than in standard cold dark matter models, and that its position is shifted to $\ell\sim 350$. We believe that our analysis can be easily extended to other types of global topological defects and general global scalar fields.

We investigate cosmological structure formation seeded by topological defects which may form during a phase transition in the early universe. First we derive a partially new, local and gauge invariant system of perturbation equations to treat microwave background and dark matter fluctuations induced by topological defects or any other type of seeds. We then show that this system is well suited for numerical analysis of structure formation by applying it to seeds induced by fluctuations of a global scalar field. Our numerical results are complementary to previous investigations since we use substantially different methods. The resulting microwave background fluctuations are compatible with older simulations. We also obtain a scale invariant spectrum of fluctuations with about the same amplitude. However, our dark matter results yield a smaller bias parameter compatible with $b\sim 2$ on a scale of $20 Mpc$ in contrast to previous work which yielded to large bias factors. Our conclusions are thus more positive. According to the aspects analyzed in this work, global topological defect induced fluctuations yield viable scenarios of structure formation and do better than standard CDM on large scales.

We discuss gravitational effects of global scalar fields and, especially, of global topological defects. We first give an introduction to the dynamics of global fields and the formation of defects. Next we investigate the induced gravitational fields, first in a flat background and then in the expanding universe. In flat space, we explicitly calculate the gravitational fields of exact global monopole and global texture solutions and discuss the motion of photons and massive particles in these geometries. We also show that slowly moving particles and the energy of photons are not affected in static scalar field configurations with vanishing potential energy. In expanding space, we explore the possibility that global topological defects from a phase transition in the very early universe may have seeded inhomogeneities in the energy distribution which yielded the observed large scale structure in the Universe, the sheets of galaxies, clusters, voids ... . We outline numerical simulations which have been performed to tackle this problem and briefly discuss their results.

The topic of this letter is structure formation with topological defects. We first present a partially new, fully local and gauge invariant system of perturbation equations to treat microwave background and dark matter fluctuations induced by topological defects (or any other type of seeds). We show that this treatment is extremly well suited for linear numerical analysis of structure formation by applying it to the texture scenario. Our numerical results cover a larger dynamical range than previous investigations and are complementary since we use substantially different methods.

We show that there are inflationary models for which perturbations in the energy momentum tensor, which are of second order in the scalar field, cannot be neglected. We first specify the conditions under which the usual first order perturbations are absent. We then analyze classically, the growth and decay of our new type of perturbations for one mode of fluctuations $\de\phi_k$ in the scalar field. We generalize this analysis, considering the contribution from the whole spectrum of $\de\phi$ to a given wavelength of geometrical perturbations. Finally, we discuss the evolution of the perturbations during the subsequent radiation dominated era and discuss the resulting spectrum of density fluctuations. In the case of a massless scalar field we find a spectral index $n=4$. For massive scalar fields we obtain $n=0$ but the resulting amplitude of fluctuations for inflation around a GUT scale are by far too high. Hence, 'conventional' inflationary models must not be influenced by this new type of perturbations, in order to lead to acceptable matter and CMB fluctuations.

In this letter a new formula for light deflection is derived using only physically observable concepts. The general result is specialized to cosmological perturbation theory and expressed in terms of gauge--invariant perturbation variables. The resulting scalar, vector and tensor equations are supplemented by simple examples for illustration. The gravity wave example may be of more than academic interest and even represent a new way to detect gravitational waves.

After an introduction to the problem of cosmological structure formation, we develop gauge invariant cosmological perturbation theory. We derive the first order perturbation equations of Einstein's equations and energy momentum ``conservation''. Furthermore, the perturbations of Liouville's equation for collisionless particles and Boltzmann's equation for Compton scattering are worked out. We fully discuss the propagation of photons in a perturbed Friedmann universe, calculating the Sachs--Wolfe effect and light deflection. The perturbation equations are extended to accommodate also perturbations induced by seeds. With these general results we discuss some of the main aspects of the texture model for the formation of large scale structure in the Universe (galaxies, clusters, sheets, voids). In this model, perturbations in the dark matter are induced by texture seeds. The gravitational effects of a spherically symmetric collapsing texture on dark matter, baryonic matter and photons are calculated in first order perturbation theory. We study the characteristic signature of the microwave background fluctuations induced in this scenario and compare it with the COBE observations.

The cosmic microwave background (CMB) anisotropies have turned out to represent one of the most stringent 'bottle necks' for scenarios of large scale structure formation. As a possibility to relax this constraint, it has been proposed that early reionization can damp CMB fluctuations on small scales due to photon diffusion in the ionized plasma. As an example, I investigate the recently proposed scenario with cold dark matter (CDM) and texture seeds. There, an analysis of CMB anisotropies shows that early reionization is a crucial ingredient for this scenario. Without damping, the small scale anisotropies would dominate and exceed observed limits. In this paper I present analytical and numerical results for the amount of damping due to early reionization for CMB perturbations induced by a collapsing texture. Furthermore, the spectral distortion of the CMB due to Compton scattering of the hotter plasma electrons is calculated. Next I discuss the physical processes which lead to a system of coupled ordinary differential equations for the degree of ionization, the electron temperature and the evolution of the ionizing radiation.

The cosmic microwave anisotropies in a scenario of large scale structure formation with cold dark matter and texture are discussed and compared with recent observational results of the COBE satellite. A couple of important statistical parameters are determined. The fluctuations are slightly non gaussian. The quadrupole anisotropy is $1.5\pm 1.2\times 10^{-5}$ and the fluctuations on a angular scale of 10 degrees are $ (3.8\pm 2.6)\times 10^{-5}$. The COBE are within about one standard deviation of the typical texture + CDM model discussed in this paper. Furthermore, we calculate fluctuations on intermediate scales (about 2 degrees) with the result $\De T/T(\theta \sim 2^o) = 3.9\pm 0.8)\times 10^{-5}$. Collapsing textures are modeled by spherically symmetric field configurations. This leads to uncertainties of about a factor of~2.