The standard cosmological model determined from the accurate cosmic microwave background measurements made by the Planck satellite implies a value of the Hubble constant $H_0$ that is $4.4$ standard deviations lower than the one determined from Type Ia supernovae. The Planck best fit model also predicts higher values of the matter density fraction $\Omega_m$ and clustering amplitude $S_8$ compared to those obtained from the Dark Energy Survey Year 1 data. Here we show that accounting for the enhanced recombination rate due to additional small-scale inhomogeneities in the baryon density may solve both the $H_0$ and the $S_8-\Omega_m$ tensions. The additional baryon inhomogeneities can be induced by primordial magnetic fields present in the plasma prior to recombination. The required field strength to solve the Hubble tension is just what is needed to explain the existence of galactic, cluster, and extragalactic magnetic fields without relying on dynamo amplification. Our results show clear evidence for this effect and motivate further detailed studies of primordial magnetic fields, setting several well-defined targets for future observations.

Astrophysics spans an enormous range of questions on scales from individual planets to the entire cosmos. To address the richness of 21st century astrophysics requires a corresponding richness of telescopes spanning all bands and all messengers. Much scientific benefit comes from having the multi-wavelength capability available at the same time. Most of these bands,or measurement sensitivities, require space-based missions. Historically, NASA has addressed this need for breadth with a small number of flagship-class missions and a larger number of Explorer missions. While the Explorer program continues to flourish, there is a large gap between Explorers and strategic missions. A fortunate combination of new astrophysics technologies with new, high capacity, low dollar-per-kg to orbit launchers, and new satellite buses allow for cheaper missions with capabilities approaching strategic mission levels. NASA has recognized these developments by calling for Probe-class mission ideas for mission studies, spanning most of the electromagnetic spectrum from GeV gamma-rays to the far infrared, and the new messengers of neutrinos and ultra-high energy cosmic rays. The key insight from the Probes exercise is that order-of-magnitude advances in science performance metrics are possible across the board for initial total cost estimates in the range 500M-1B dollars.

We numerically study the decay of cosmic global string loops due to radiation of Goldstone bosons and massive scalar ($\chi$) particles. The length of loops we study range from 200-1000 times the width of the string core. We find that the lifetime of a loop is $\approx 1.4 L$. The energy spectrum of Goldstone boson radiation has a $k^{-1}$ fall off, where $k$ is the wavenumber, and a sharp peak at $k \approx m_\chi/2$, where $m_\chi$ is the mass of $\chi$. The latter is a new feature and implies a peak at high energies (MeV-GeV) in the cosmological distribution of QCD axions.

The Probe of Inflation and Cosmic Origins (PICO) is a proposed probe-scale space mission consisting of an imaging polarimeter operating in frequency bands between 20 and 800 GHz. We describe the science achievable by PICO, which has sensitivity equivalent to more than 3300 Planck missions, the technical implementation, the schedule and cost.

We provide an overview of the science case, instrument configuration and project plan for the next-generation ground-based cosmic microwave background experiment CMB-S4, for consideration by the 2020 Decadal Survey.

The expansion history of the Universe reconstructed from a combination of recent data indicates a preference for a changing Dark Energy (DE) density. Moreover, the DE density appears to be increasing with cosmic time, with its equation of state being below -1 on average, and possibly crossing the so-called phantom divide. Scalar-tensor theories, in which the scalar field mediates a force between matter particles, offer a natural framework in which the effective DE equation of state can be less than -1 and cross the phantom barrier. We consider the generalized Brans-Dicke (GBD) class of scalar-tensor theories and reconstruct their Lagrangian given the effective DE density extracted from recent data. Then, given the reconstructed Lagrangian, we solve for the linear perturbations and investigate the characteristic signatures of these reconstructed GBD in the cosmological observables, such as the cosmic microwave background (CMB) anisotropy, the galaxy number counts, and their cross-correlations. In particular, we demonstrate that the Integrated Sachs-Wolfe (ISW) effect probed by the cross-correlation of CMB with the matter distribution can rule out scalar-tensor theories as the explanation of the observed DE dynamics independently from the laboratory and solar system fifth force constraints.

The primary scientific target of the CMB polarization experiments that are currently being built and proposed is the detection of primordial tensor perturbations. As a byproduct, these instruments will significantly improve constraints on cosmic birefringence, or the rotation of the CMB polarization plane. If convincingly detected, cosmic birefringence would be a dramatic manifestation of physics beyond the standard models of particle physics and cosmology. We forecast the bounds on the cosmic polarization rotation (CPR) from the upcoming ground-based Simons Observatory (SO) and the space-based LiteBIRD experiments, as well as a "fourth generation" ground-based CMB experiment like CMB-S4 and the mid-cost space mission PICO. We examine the detectability of both a stochastic anisotropic rotation field and an isotropic rotation by a constant angle. CPR induces new correlations of CMB observables, including spectra of parity-odd type in the case of isotropic CPR, and mode-coupling correlations in the anisotropic rotation case. We find that LiteBIRD and SO will reduce the 1$\sigma$ bound on the isotropic CPR from the current value of 30 arcmin to 1.5 and 0.6 arcmin, respectively, while CMB-S4-like and PICO will reduce it to $\sim 0.1$ arcmin. The bounds on the amplitude of a scale-invariant CPR spectrum will be reduced by 1, 2 and 3 orders of magnitude by LiteBIRD, SO and CMB-S4-like/PICO, respectively. We discuss implications of the forecasted CPR bounds for pseudoscalar fields, primordial magnetic fields (PMF), and violations of Lorentz invariance. We find that CMB-S4-like and PICO can reduce the 1$\sigma$ bound on the amplitude of the scale-invariant PMF from 1 nG to 0.1 nG, while also probing the magnetic field of the Milky Way. They will also significantly improve bounds on the axion-photon coupling, placing stringent constraints on the string theory axions.

We present the science case, reference design, and project plan for the Stage-4 ground-based cosmic microwave background experiment CMB-S4.

Constraints on the tension and the abundance of cosmic strings depend crucially on the rate at which they decay into particles and gravitational radiation. We study the decay of cosmic string loops in the Abelian-Higgs model by performing field theory simulations of loop formation and evolution. We find that our set of string loops emit particle radiation primarily due to kink collisions, and that their decay time due to these losses is proportional to $L^p$ with $p \approx 2$ where $L$ is the loop length. In contrast, the decay time to gravitational radiation scales in proportion to $L$, and we conclude that particle emission is the primary energy loss mechanism for loops smaller than a critical length scale, while gravitational losses dominate for larger loops.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

While Type Ia Supernovae (SNe Ia) are one of the most mature cosmological probes, the next era promises to be extremely exciting in the number of different ways SNe Ia are used to measure various cosmological parameters. Here we review the experiments in the 2020s that will yield orders of magnitudes more SNe Ia, and the new understandings and capabilities to constrain systematic uncertainties at a level to match these statistics. We then discuss five different cosmological probes with SNe Ia: the conventional Hubble diagram for measuring dark energy properties, the distance ladder for measuring the Hubble constant, peculiar velocities and weak lensing for measuring sigma8 and strong-lens measurements of H0 and other cosmological parameters. For each of these probes, we discuss the experiments that will provide the best measurements and also the SN Ia-related systematics that affect each one.

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

The Probe of Inflation and Cosmic Origins (PICO) is an imaging polarimeter that will scan the sky for 5 years in 21 frequency bands spread between 21 and 799 GHz. It will produce full-sky surveys of intensity and polarization with a final combined-map noise level of 0.87 $\mu$K arcmin for the required specifications, equivalent to 3300 Planck missions, and with our current best-estimate would have a noise level of 0.61 $\mu$K arcmin (6400 Planck missions). PICO will either determine the energy scale of inflation by detecting the tensor to scalar ratio at a level $r=5\times 10^{-4}~(5\sigma)$, or will rule out with more than $5\sigma$ all inflation models for which the characteristic scale in the potential is the Planck scale. With LSST's data it could rule out all models of slow-roll inflation. PICO will detect the sum of neutrino masses at $>4\sigma$, constrain the effective number of light particle species with $\Delta N_{\rm eff}<0.06~(2\sigma)$, and elucidate processes affecting the evolution of cosmic structures by measuring the optical depth to reionization with errors limited by cosmic variance and by constraining the evolution of the amplitude of linear fluctuations $\sigma_{8}(z)$ with sub-percent accuracy. Cross-correlating PICO's map of the thermal Sunyaev-Zeldovich effect with LSST's gold sample of galaxies will precisely trace the evolution of thermal pressure with $z$. PICO's maps of the Milky Way will be used to determine the make up of galactic dust and the role of magnetic fields in star formation efficiency. With 21 full sky legacy maps in intensity and polarization, which cannot be obtained in any other way, the mission will enrich many areas of astrophysics. PICO is the only single-platform instrument with the combination of sensitivity, angular resolution, frequency bands, and control of systematic effects that can deliver this compelling, timely, and broad science.

We present a major upgrade of MGCAMB, a patch for the Einstein-Boltzmann solver CAMB used for phenomenological tests of general relativity against cosmological datasets. This new version is compatible with the latest CosmoMC code and includes a consistent implementation of massive neutrinos and dynamical dark energy. The code has been restructured to make it easier to modify to fit the custom needs of specific users. To illustrate the capabilities of the code, we present joint constraints on the modified growth, massive neutrinos and the dark energy equation of state from the latest cosmological observations, including the recent galaxy counts and weak lensing measurements from the Dark Energy Survey, and find a good consistency with the {\Lambda}CDM model.

Ongoing and upcoming cosmological surveys will significantly improve our ability to probe the equation of state of dark energy, $w_{\rm DE}$, and the phenomenology of Large Scale Structure. They will allow us to constrain deviations from the $\Lambda$CDM predictions for the relations between the matter density contrast and the weak lensing and the Newtonian potential, described by the functions $\Sigma$ and $\mu$, respectively. In this work, we derive the theoretical prior for the joint covariance of $w_{\rm DE}$, $\Sigma$ and $\mu$, expected in general scalar-tensor theories with second order equations of motion (Horndeski gravity), focusing on their time-dependence at certain representative scales. We employ Monte-Carlo methods to generate large ensembles of statistically independent Horndeski models, focusing on those that are physically viable and in broad agreement with local tests of gravity, the observed cosmic expansion history and the measurement of the speed of gravitational waves from a binary neutron star merger. We identify several interesting features and trends in the distribution functions of $w_{\rm DE}$, $\Sigma$ and $\mu$, as well as in their covariances; we confirm the high degree of correlation between $\Sigma$ and $\mu$ in scalar-tensor theories. The derived prior covariance matrices will allow us to reconstruct jointly $w_{\rm DE}$, $\Sigma$ and $\mu$ in a non-parametric way.

We reconstruct evolution of the dark energy (DE) density using a nonparametric Bayesian approach from a combination of latest observational data. We caution against parameterizing DE in terms of its equation of state as it can be singular in modified gravity models, and using it introduces a bias preventing negative effective DE densities. We find a $3.7\sigma$ preference for an evolving effective DE density with interesting features. For example, it oscillates around the $\Lambda$CDM prediction at $z\lesssim0.7$, and could be negative at $z\gtrsim2.3$; dark energy can be pressure-less at multiple redshifts, and a short period of cosmic deceleration is allowed at $0.1 \lesssim z\lesssim 0.2$. We perform the reconstruction for several choices of the prior, as well as a evidence-weighted reconstruction. We find that some of the dynamical features, such as the oscillatory behaviour of the DE density, are supported by the Bayesian evidence, which is a first detection of a dynamical DE with a positive Bayesian evidence. The evidence-weighted reconstruction prefers a dynamical DE at a $(2.5\pm0.06)\sigma$ significance level.

Was the primordial universe magnetized? The answer to this question would help explain the origin of micro-Gauss strength magnetic fields observed in galaxies. It is also of fundamental importance in developing a complete theory of the early universe. While there can be other signatures of cosmological magnetic fields, a signature in the cosmic microwave background (CMB) would prove their primordial origin. The B-mode polarization of CMB is particularly promising in this regard because there are relatively few other sources of B-modes, and because the vortical modes sourced by the primordial magnetic field (PMF) survive diffusion damping up to a small fraction of the Silk length. At present, the Planck temperature and polarization spectra combined with the B-mode spectrum measured by the South Pole Telescope (SPT) constrain the PMF strength to be no more than $\sim 1$ nano-Gauss (nG). Because of the quartic scaling of the CMB anisotropy spectra with the PMF strength, this bound will not change by much even with the significantly better measurements of the B-mode spectrum by the Stage III and Stage IV CMB experiments. On the other hand, tightening the bound well below the $1$ nG threshold would rule out the purely primordial origin (requiring no dynamo action) of galactic fields. Considering Faraday rotation, which converts some of the E-modes into B-modes and scales linearly with the field strength, will help to achieve this goal. As we demonstrate, the upcoming experiments, such as SPT-3G and the Simons Observatory, will be sensitive to fields of $\sim 0.5$ nG strength thanks to the mode-coupling signature induced by Faraday rotation. A future Stage IV ground based experiment or a Space Probe will be capable of probing fields below $0.1$ nG, and would detect a scale-invariant PMF of $0.2$ nG strength without de-lensing or subtracting the galactic rotation measure.

Phenomenological functions $\Sigma$ and $\mu$, also known as $G_{\rm light}/G$ and $G_{\rm matter}/G$, are commonly used to parameterize modifications of the growth of large-scale structure in alternative theories of gravity. We study the values these functions can take in Horndeski theories, i.e. the class of scalar-tensor theories with second order equations of motion. We restrict our attention to models that are in a broad agreement with tests of gravity and the observed cosmic expansion history. In particular, we require the speed of gravity to be equal to the speed of light today, as required by the recent detection of gravitational waves and electromagnetic emission from a binary neutron star merger. We examine the correlations between the values of $\Sigma$ and $\mu$ analytically within the quasi-static approximation, and numerically, by sampling the space of allowed solutions. We confirm that the conjecture made in [Pogosian:2016pwr], that $(\Sigma-1)(\mu -1) \ge 0$ in viable Horndeski theories, holds very well. Along with that, we check the validity of the quasi-static approximation within different corners of Horndeski theory. Our results show that, even with the tight bound on the present day speed of gravitational waves, there is room within Horndeski theories for non-trivial signatures of modified gravity at the level of linear perturbations.

A flat Friedman-Roberson-Walker universe dominated by a cosmological constant ($\Lambda$) and cold dark matter (CDM) has been the working model preferred by cosmologists since the discovery of cosmic acceleration. However, tensions of various degrees of significance are known to be present among existing datasets within the $\Lambda$CDM framework. In particular, the Lyman-$\alpha$ forest measurement of the Baryon Acoustic Oscillations (BAO) by the Baryon Oscillation Spectroscopic Survey (BOSS) prefers a smaller value of the matter density fraction $\Omega_{\rm M}$ compared to the value preferred by cosmic microwave background (CMB). Also, the recently measured value of the Hubble constant, $H_0=73.24\pm1.74 \ {\rm km}\ {\rm s}^{-1} \ {\rm Mpc}^{-1}$, is $3.4\sigma$ higher than $66.93\pm0.62 \ {\rm km}\ {\rm s}^{-1} \ {\rm Mpc}^{-1}$ inferred from the Planck CMB data. In this work, we investigate if these tensions can be interpreted as evidence for a non-constant dynamical dark energy (DE). Using the Kullback-Leibler (KL) divergence to quantify the tension between datasets, we find that the tensions are relieved by an evolving DE, with the dynamical DE model preferred at a $3.5\sigma$ significance level based on the improvement in the fit alone. While, at present, the Bayesian evidence for the dynamical DE is insufficient to favour it over $\Lambda$CDM, we show that, if the current best fit DE happened to be the true model, it would be decisively detected by the upcoming DESI survey.

Constraining the Dark Energy (DE) equation of state, w, is one of the primary science goals of ongoing and future cosmological surveys. In practice, with imperfect data and incomplete redshift coverage, this requires making assumptions about the evolution of w with redshift z. These assumptions can be manifested in a choice of a specific parametric form, which can potentially bias the outcome, or else one can reconstruct w(z) non-parametrically, by specifying a prior covariance matrix that correlates values of w at different redshifts. In this work, we derive the theoretical prior covariance for the effective DE equation of state predicted by general scalar-tensor theories with second order equations of motion (Horndeski theories). This is achieved by generating a large ensemble of possible scalar-tensor theories using a Monte Carlo methodology, including the application of physical viability conditions. We also separately consider the special sub-case of the minimally coupled scalar field, or quintessence. The prior shows a preference for tracking behaviors in the most general case. Given the covariance matrix, theoretical priors on parameters of any specific parametrization of w(z) can also be readily derived by projection.

Phenomenological functions $\Sigma$ and $\mu$ (also known as $G_{\rm light}/G$ and $G_{\rm matter}/G$) are commonly used to parameterize possible modifications of the Poisson equation relating the matter density contrast to the lensing and the Newtonian potentials, respectively. They will be well constrained by future surveys of large scale structure. But what would the implications of measuring particular values of these functions be for modified gravity theories? We ask this question in the context of general Horndeski class of single field scalar-tensor theories with second order equations of motion. We find several consistency conditions that make it possible to rule out broad classes of theories based on measurements of $\Sigma$ and $\mu$ that are independent of their parametric forms. For instance, a measurement of $\Sigma \ne 1$ would rule out all models with a canonical form of kinetic energy, while finding $\Sigma-1$ and $\mu-1$ to be of opposite sign would strongly disfavour the entire class of Horndeski models. We separately examine the large and the small scale limits, the possibility of scale-dependence, and the consistency with bounds on the speed of gravitational waves. We identify sub-classes of Horndeski theories that can be ruled out based on the measured difference between $\Sigma$ and $\mu$.

A primordial magnetic field (PMF) present before recombination can leave specific signatures on the cosmic microwave background (CMB) fluctuations. Of particular importance is its contribution to the B-mode polarization power spectrum. Indeed, vortical modes sourced by the PMF can dominate the B-mode power spectrum on small scales, as they survive damping up to a small fraction of the Silk length. Therefore, measurements of the B-mode polarization at high-$\ell$ , such as the one recently performed by the South Pole Telescope (SPT), have the potential to provide stringent constraints on the PMF. We use the publicly released SPT B-mode polarization spectrum, along with the temperature and polarization data from the Planck satellite, to derive constraints on the magnitude, the spectral index and the energy scale at which the PMF was generated. We find that, while Planck data constrains the magnetic amplitude to $B_{1 \, \text{Mpc}} < 3.3$ nG at 95\% confidence level (CL), the SPT measurement improves the constraint to $B_{1 \, \text{Mpc}} < 1.5$ nG. The magnetic spectral index, $n_B$, and the time of the generation of the PMF are unconstrained. For a nearly scale-invariant PMF, predicted by simplest inflationary magnetogenesis models, the bound from Planck+SPT is $B_{1 \, \text{Mpc}} < 1.2$ nG at 95% CL. For PMF with $n_B=2$, expected for fields generated in post-inflationary phase transitions, the 95% CL bound is $B_{1 \, \text{Mpc}} < 0.002$ nG, corresponding to the magnetic fraction of the radiation density $\Omega_{B\gamma} < 10^{-3}$ or the effective field $B_{\rm eff} < 100$ nG. The patches for the Boltzmann code CAMB and the Markov Chain Monte Carlo engine CosmoMC, incorporating the PMF effects on CMB, are made publicly available.

This book lays out the scientific goals to be addressed by the next-generation ground-based cosmic microwave background experiment, CMB-S4, envisioned to consist of dedicated telescopes at the South Pole, the high Chilean Atacama plateau and possibly a northern hemisphere site, all equipped with new superconducting cameras. CMB-S4 will dramatically advance cosmological studies by crossing critical thresholds in the search for the B-mode polarization signature of primordial gravitational waves, in the determination of the number and masses of the neutrinos, in the search for evidence of new light relics, in constraining the nature of dark energy, and in testing general relativity on large scales.

We present a science forecast for the eBOSS survey, part of the SDSS-IV project, which is a spectroscopic survey using multiple tracers of large-scale structure, including luminous red galaxies (LRGs), emission line galaxies (ELGs) and quasars (both as a direct probe of structure and through the Ly-$\alpha$ forest). Focusing on discrete tracers, we forecast the expected accuracy of the baryonic acoustic oscillation (BAO), the redshift-space distortion (RSD) measurements, the $f_{\rm NL}$ parameter quantifying the primordial non-Gaussianity, the dark energy and modified gravity parameters. We also use the line-of-sight clustering in the Ly-$\alpha$ forest to constrain the total neutrino mass. We find that eBOSS LRGs ($0.6<z<1.0$) (combined with the BOSS LRGs at $z>0.6$), ELGs ($0.6<z<1.2$) and Clustering Quasars (CQs) ($0.6<z<2.2$) can achieve a precision of 1%, 2.2% and 1.6% precisions, respectively, for spherically averaged BAO distance measurements. Using the same samples, the constraint on $f\sigma_8$ is expected to be 2.5%, 3.3% and 2.8% respectively. For primordial non-Gaussianity, eBOSS alone can reach an accuracy of $\sigma(f_{\rm NL})\sim10-15$, depending on the external measurement of the galaxy bias and our ability to model large-scale systematic errors. eBOSS can at most improve the dark energy Figure of Merit (FoM) by a factor of $3$ for the Chevallier-Polarski-Linder (CPL) parametrisation, and can well constrain three eigenmodes for the general equation-of-state parameter (Abridged).

We constrain anisotropic cosmic birefringence using four-point correlations of even-parity $E$-mode and odd-parity $B$-mode polarization in the cosmic microwave background measurements made by the POLARization of the Background Radiation (POLARBEAR) experiment in its first season of observations. We find that the anisotropic cosmic birefringence signal from any parity-violating processes is consistent with zero. The Faraday rotation from anisotropic cosmic birefringence can be compared with the equivalent quantity generated by primordial magnetic fields if they existed. The POLARBEAR nondetection translates into a 95% confidence level (C.L.) upper limit of 93 nanogauss (nG) on the amplitude of an equivalent primordial magnetic field inclusive of systematic uncertainties. This four-point correlation constraint on Faraday rotation is about 15 times tighter than the upper limit of 1380 nG inferred from constraining the contribution of Faraday rotation to two-point correlations of $B$-modes measured by Planck in 2015. Metric perturbations sourced by primordial magnetic fields would also contribute to the $B$-mode power spectrum. Using the POLARBEAR measurements of the $B$-mode power spectrum (two-point correlation), we set a 95% C.L. upper limit of 3.9 nG on primordial magnetic fields assuming a flat prior on the field amplitude. This limit is comparable to what was found in the Planck 2015 two-point correlation analysis with both temperature and polarization. We perform a set of systematic error tests and find no evidence for contamination. This work marks the first time that anisotropic cosmic birefringence or primordial magnetic fields have been constrained from the ground at subdegree scales.

An interaction between the vacuum energy and dark matter is an intriguing possibility which may offer a way of solving the cosmological constant problem. Adopting a general prescription for momentum exchange between the two dark components, we reconstruct $\alpha(a)$, the temporal evolution of the coupling strength between dark matter and vacuum energy, in a nonparametric Bayesian approach using combined observational data sets from the cosmic microwave background, supernovae and large scale structure. An evolving interaction between the vacuum energy and dark matter removes some of the tensions between different data sets. However, it is not preferred over $\Lambda$CDM in the Bayesian sense, as improvement in the fit is not sufficient to compensate for the increase in the volume of the parameter space.

Modified gravity theories often contain a scalar field of gravitational strength which interacts with matter. We examine constraints on the range and the coupling strength of a scalar gravitational degree of freedom using a subset of current data that can be safely analyzed within the linear perturbation theory. Using a model-independent implementation of scalar-tensor theories in MGCAMB in terms of two functions of the scale factor describing the mass and the coupling of the scalar degree of freedom, we derive constraints on the $f(R)$, generalized chameleon, Symmetron and Dilaton models. Since most of the large scale structure data available today is from relatively low redshifts, only a limited range of observed scales is in the linear regime, leading to relatively weak constraints. We then perform a forecast for a future large scale structure survey, such as Large Synoptic Survey Telescope (LSST), which will map a significant volume at higher redshifts, and show that it will produce much stronger constraints on scalar interactions in specific models. We also perform a principal component analysis and find that future surveys should be able to provide tight constraints on several eigenmodes of the scalar mass evolution.

Interactions of different types of topological defects can play an important role in the aftermath of a phase transition. We study interactions of fundamental magnetic monopoles and stable domain walls in a Grand Unified theory in which $SU(5) \times Z_2$ symmetry is spontaneously broken to $SU(3)\times SU(2)\times U(1)/Z_6$. We find that there are only two distinct outcomes depending on the relative orientation of the monopole and the wall in internal space. In one case, the monopole passes through the wall, while in the other it unwinds on hitting the wall.

This article, written for Scolarpedia, provides a brief introduction into the subject of cosmic strings, together with a review of their main properties, cosmological evolution and observational signatures.

We study properties of domain walls in the symmetron model, in which the scalar gravitational degree of freedom decouples from matter in regions of high density, and exhibits a spontaneously broken $Z_2$ symmetry at low densities. The non-minimal coupling of the scalar to matter leads to a host of interesting properties of the domain walls that are not present in minimally coupled theories. We estimate the cosmological energy fraction in domain walls and find that this leads to an upper bound on the redshift of the symmetry breaking. We also show that a spherical symmetron wall can remain stable if it is "pinned" on matter halos and derive a criterion for the stability. In addition, we present results of numerical simulations of representative interactions between domain walls and matter over-densities.

Scaling networks of cosmic defects, such as strings and textures, actively generate scalar, vector and tensor metric perturbations throughout the history of the universe. In particular, {\em vector} modes sourced by defects are an efficient source of the CMB B-mode polarization. We use the recently released BICEP2 and POLARBEAR B-mode polarization spectra to constrain properties of a wide range of different types of cosmic strings networks. We find that in order for strings to provide a satisfactory fit on their own, the effective inter-string distance needs to be extremely large -- spectra that fit the data best are more representative of global strings and textures. When a local string contribution is considered together with the inflationary B-mode spectrum, the fit is improved. We discuss implications of these results for theories that predict cosmic defects.

At linear order in cosmological perturbations, departures from the growth in the cosmological standard model can be quantified in terms of two functions of redshift $z$ and Fourier number $k$. Previous studies have performed principal component forecasts for several choices of these two functions, based on expected capabilities of upcoming large structure surveys. It is typically found that there will be many well-constrained degrees of freedom. However, not all and, probably most, of these degrees of freedom were physical if the parametrization had allowed for an arbitrary $k$-dependence. In this paper, we restrict the $k$-dependence to that allowed in local theories of gravity under the quasi-static approximation, i.e. ratios of polynomials in $k$, and identify the best constrained features in the ($z$,$k$)-dependence of the commonly considered functions $\mu$ and $\gamma$ as measured by an LSST-like weak lensing survey. We estimate the uncertainty in the measurements of the eigenmodes of modified growth. We find that imposing the theoretical prior on $k$-dependence reduces the number of degrees of freedom and the covariance between parameters. On the other hand, imaging surveys like LSST are not as sensitive to the $z$-dependence as they are to the $k$-dependence of the modified growth functions. This trade off provides us with, more or less, the same number of well-constrained eigenmodes (with respect to our prior) as found before, but now these modes are physical.

Bounds on the amplitude of a scale-invariant stochastic primordial magnetic field (PMF) can be significantly improved by measurements of the Faraday Rotation (FR) of CMB polarization. The mode-coupling correlations induced by FR make it possible to extract it from cross-correlations of the B-mode polarization with the E-mode and the temperature anisotropy. In this paper, we construct an estimator of the rotation measure that appropriately combines measurements of the FR from multiple frequency channels. We study the dependence of the signal-to-noise in the PMF detection on the resolution and the noise of the detectors, as well as the removal of the weak lensing contribution and the galactic FR. We show that a recently proposed space-based experiment PRISM can detect magnetic fields of 0.1 nano-Gauss strength at a 2 sigma level. Higher detection levels can be achieved by reducing the detector noise and improving the resolution or increasing the number of channels in the 30-70 GHz frequency range.

Faraday Rotation (FR) of CMB polarization, as measured through mode-coupling correlations of E and B modes, can be a promising probe of a stochastic primordial magnetic field (PMF). While the existence of a PMF is still hypothetical, there will certainly be a contribution to CMB FR from the magnetic field of the Milky Way. We use existing estimates of the Milky Way rotation measure (RM) to forecast its detectability with upcoming and future CMB experiments. We find that the galactic RM will not be seen in polarization measurements by Planck, but that it will need to be accounted for by CMB experiments capable of detecting the weak lensing contribution to the B-mode. We then discuss prospects for constraining the PMF in the presence of FR due to the galaxy under various assumptions that include partial de-lensing and partial subtraction of the galactic FR. We find that a realistic future sub-orbital experiment, covering a patch of the sky near the galactic poles, can detect a scale-invariant PMF of 0.1 nano-Gauss at better than 95% confidence level, while a dedicated space-based experiment can detect even smaller fields.

We perform a principal component analysis to assess ability of future observations to measure departures from General Relativity in predictions of the Poisson and anisotropy equations on linear scales. In particular, we focus on how the measurements of redshift-space distortions (RSD) observed from spectroscopic galaxy redshift surveys will improve the constraints when combined with lensing tomographic surveys. Assuming a Euclid-like galaxy imaging and redshift survey, we find that adding the 3D information decreases the statistical uncertainty by a factor between 3 and 10 compared to the case when only observables from lensing tomographic surveys are used. We also find that the number of well-constrained modes increases by a factor between 3 and 7. Our study indicates the importance of joint galaxy imaging and redshift surveys such as SuMIRe and Euclid to give more stringent tests of the {\Lambda}CDM model and to distinguish between various modified gravity and dark energy models.

The next generation of large scale surveys will not only measure cosmological parameters within the framework of General Relativity, but will also allow for precision tests of the framework itself. At the order of linear perturbations, departures from the growth in the LCDM model can be quantified in terms of two functions of time and Fourier number k. We argue that in local theories of gravity, in the quasi-static approximation, these functions must be ratios of polynomials in k, with the numerator of one function being equal to the denominator of the other. Moreover, the polynomials are even and of second degree in practically all viable models considered today. This means that, without significant loss of generality, one can use data to constraint only five functions of a single variable, instead of two functions of two variables. Furthermore, since the five functions are expected to be slowly varying, one can fit them to data in a non-parametric way with the aid of an explicit smoothness prior. We discuss practical application of this parametrization to forecasts and fits.

Primordial magnetic fields (PMF) can create polarization $B$-modes in the cosmic microwave background (CMB) through Faraday rotation (FR), leading to non-trivial 2-point and 4-point correlators of the CMB temperature and polarization. We discuss the detectability of primordial magnetic fields using different correlators and evaluate their relative merits. We have fully accounted for the contamination by weak lensing, which contributes to the variance, but whose contribution to the 4-point correlations is orthogonal to that of FR. We show that a Planck-like experiment can detect scale-invariant PMF of nG strength using the FR diagnostic at 90GHz, while realistic future experiments at the same frequency can detect 10^{-10} G. Utilizing multiple frequencies will improve on these prospects, making FR of CMB a powerful probe of scale-invariant PMF.

We examine the accuracy of the quasi-static approximation and a parametric form commonly employed for evolving linear perturbations in f(R) gravity theories and placing cosmological constraints. In particular, we analyze the nature and the importance of the near horizon effects that are often neglected. We find that for viable models, with a small present value of the scalaron field, such corrections are entirely negligible with no impact on observables. We also find that the one-parameter form, commonly used to represent the modified equations for linear perturbations in f(R), leads to theoretical systematic errors that are relatively small and, therefore, is adequate for placing constraint on f(R) models.

Cosmic Microwave Background (CMB) polarization B-modes induced by Faraday Rotation (FR) can provide a distinctive signature of primordial magnetic fields because of their characteristic frequency dependence and because they are only weakly damped on small scales. FR also leads to mode-coupling correlations between the E and B type polarization, and between the temperature and the B-mode. These additional correlations can further help distinguish magnetic fields from other sources of B-modes. We review the FR induced CMB signatures and present the constraints on primordial magnetism that can be expected from upcoming CMB experiments. Our results suggest that FR of CMB will be a promising probe of primordial magnetic fields.

We apply a new non-parametric Bayesian method for reconstructing the evolution history of the equation-of-state $w$ of dark energy, based on applying a correlated prior for $w(z)$, to a collection of cosmological data. We combine the latest supernova (SNLS 3-year or Union2.1), cosmic microwave background, redshift space distortion and the baryonic acoustic oscillation measurements (including BOSS, WiggleZ and 6dF) and find that the cosmological constant appears consistent with current data, but that a dynamical dark energy model which evolves from $w<-1$ at $z\sim0.25$ to $w > -1$ at higher redshift is mildly favored. Estimates of the Bayesian evidences show little preference between the cosmological constant model and the dynamical model for a range of correlated prior choices. Looking towards future data, we find that the best fit models for current data could be well distinguished from the $\Lambda$CDM model by observations such as Planck and Euclid-like surveys.

The next generation of weak lensing surveys will trace the evolution of matter perturbations and gravitational potentials from the matter dominated epoch until today. Along with constraining the dynamics of dark energy, they will probe the relations between matter overdensities, local curvature, and the Newtonian potential. We work with two functions of time and scale to account for any modifications of these relations in the linear regime from those in the LCDM model. We perform a Principal Component Analysis (PCA) to find the eigenmodes and eigenvalues of these functions for surveys like DES and LSST. This paper builds on and significantly extends the PCA analysis of Zhao et al. (2009) in several ways. In particular, we consider the impact of some of the systematic effects expected in weak lensing surveys. We also present the PCA in terms of other choices of the two functions needed to parameterize modified growth on linear scales, and discuss their merits. We analyze the degeneracy between the modified growth functions and other cosmological parameters, paying special attention to the effective equation of state w(z). Finally, we demonstrate the utility of the PCA as an efficient data compression stage which enables one to easily derive constraints on parameters of specific models without recalculating Fisher matrices from scratch.

We develop an efficient, non-parametric Bayesian method for reconstructing the time evolution of the dark energy equation of state w(z) from observational data. Of particular importance is the choice of prior, which must be chosen carefully to minimise variance and bias in the reconstruction. Using a principal component analysis, we show how a correlated prior can be used to create a smooth reconstruction and also avoid bias in the mean behaviour of w(z). We test our method using Wiener reconstructions based on Fisher matrix projections, and also against more realistic MCMC analyses of simulated data sets for Planck and a future space-based dark energy mission. While the accuracy of our reconstruction depends on the smoothness of the assumed w(z), the relative error for typical dark energy models is <10% out to redshift z=1.5.

Faraday rotation induced B-modes can provide a distinctive signature of primordial magnetic fields because of their characteristic frequency dependence and because they are only weakly damped on small scales, allowing them to dominate B-modes from other sources. By numerically solving the full CMB radiative transport equations, we study the B-mode power spectrum induced by stochastic magnetic fields that have significant power on scales smaller than the thickness of the last scattering surface. Constraints on the magnetic field energy density and inertial scale are derived from WMAP 7-year data, and are stronger than the big bang nucleosynthesis (BBN) bound for a range of parameters. Observations of the CMB polarization at smaller angular scales are crucial to provide tighter constraints or a detection.

We study signatures of cosmic superstring networks containing strings of multiple tensions and Y-junctions, on the cosmic microwave background (CMB) temperature and polarisation spectra. Focusing on the crucial role of the string coupling constant $g_s$, we show that the number density and energy density of the scaling network are dominated by different types of string in the $g_s \sim 1$ and $g_s \ll 1$ limits. This can lead to an observable shift in the position of the B-mode peak --- a distinct signal leading to a direct constraint on $g_s$. We forecast the joint bounds on $g_s$ and the fundamental string tension $\mu_F$ from upcoming and future CMB polarisation experiments, as well as the signal to noise in detecting the difference between B-mode signals in the limiting cases of large and small $g_s$. We show that such a detectable shift is within reach of planned experiments.

We introduce a patch to the commonly used public codes CAMB and CosmoMC that allows the user to implement a general modification of the equations describing the growth of cosmological perturbations, while preserving the covariant conservation of the energy-momentum. This patch replaces the previously publicly released code MGCAMB, while also extending it in several ways. The new version removes the limitation of late-time-only modifications to the perturbed Einstein equations, and includes several parametrization introduced in the literature. To demonstrate the use of the patch, we obtain joint constraints on the neutrino mass and parameters of a scalar-tensor gravity model from CMB, SNe and ISW data as measured from the correlation of CMB with large scale structure.

We review the existence, formation and properties of cosmic strings in string theory, the wide variety of observational techniques that are being employed to detect them, and the constraints that current observations impose on string theory models.

The next generation of weak lensing surveys will trace the growth of large scale perturbations through a sequence of epochs, offering an opportunity to test General Relativity (GR) on cosmological scales. We review in detail the parametrization used in MGCAMB to describe the modified growth expected in alternative theories of gravity and generalized dark energy models. We highlight its advantages and examine several theoretical aspects. In particular, we show that the same set of equations can be consistently used on super-horizon and sub-horizon linear scales. We also emphasize the sensitivity of data to scale-dependent features in the growth pattern, and propose using Principal Component Analysis to converge on a practical set of parameters which is most likely to detect departures from GR. The connection with other parametrizations is also discussed.

We study the cosmic microwave background temperature and polarisation spectra sourced by multi-tension cosmic superstring networks. First we obtain solutions for the characteristic length scales and velocities associated with the evolution of a network of F-D strings, allowing for the formation of junctions between strings of different tensions. We find two distinct regimes describing the resulting scaling distributions for the relative densities of the different types of strings, depending on the magnitude of the fundamental string coupling g_s. In one of them, corresponding to the value of the coupling being of order unity, the network's stress-energy power spectrum is dominated by populous light F and D strings, while the other regime, at smaller values of g_s, has the spectrum dominated by rare heavy D strings. These regimes are seen in the CMB anisotropies associated with the network. We focus on the dependence of the shape of the B-mode polarisation spectrum on g_s and show that measuring the peak position of the B-mode spectrum can point to a particular value of the string coupling. Finally, we assess how this result, along with pulsar bounds on the production of gravitational waves from strings, can be used to constrain a combination of g_s and the fundamental string tension mu_F. Since CMB and pulsar bounds constrain different combinations of the string tensions and densities, they result in distinct shapes of bounding contours in the (mu_F, g_s) parameter plane, thus providing complementary constraints on the properties of cosmic superstrings.

We explore the complementarity of weak lensing and galaxy peculiar velocity measurements to better constrain modifications to General Relativity. We find no evidence for deviations from GR on cosmological scales from a combination of peculiar velocity measurements (for Luminous Red Galaxies in the Sloan Digital Sky Survey) with weak lensing measurements (from the CFHT Legacy Survey). We provide a Fisher error forecast for a Euclid-like space-based survey including both lensing and peculiar velocity measurements, and show that the expected constraints on modified gravity will be at least an order of magnitude better than with present data, i.e. we will obtain 5% errors on the modified gravity parametrization described here. We also present a model--independent method for constraining modified gravity parameters using tomographic peculiar velocity information, and apply this methodology to the present dataset.

We test General Relativity (GR) using current cosmological data: the cosmic microwave background (CMB) from WMAP5 (Komatsu et al. 2009), the integrated Sachs-Wolfe (ISW) effect from the cross-correlation of the CMB with six galaxy catalogs (Giannantonio et al. 2008), a compilation of supernovae Type Ia (SNe) including the latest SDSS SNe (Kessler et al. 2009), and part of the weak lensing (WL) data from CFHTLS (Fu et al. 2008, Kilbinger et al. 2009) that probe linear and mildly non-linear scales. We first test a model where the effective Newton's constant, mu, and the ratio of the two gravitational potentials, eta, transit from the GR value to another constant at late times; in this case, we find that standard GR is fully consistent with the combined data. The strongest constraint comes from the ISW effect which would arise from this gravitational transition; the observed ISW signal imposes a tight constraint on a combination of mu and eta that characterizes the lensing potential. Next, we consider four pixels in time and space for each function mu and eta, and perform a Principal Component Analysis (PCA) finding that seven of the resulting eight eigenmodes are consistent with GR within the errors. Only one eigenmode shows a 2-sigma deviation from the GR prediction, which is likely to be due to a systematic effect. However, the detection of such a deviation demonstrates the power of our time- and scale-dependent PCA methodology when combining observations of structure formation and expansion history to test GR.

We study the utility of wavelets for detecting the redshift evolution of the dark energy equation of state w(z) from the combination of supernovae, CMB and BAO data. We show that local features in w, such as bumps, can be detected efficiently using wavelets. To demonstrate, we first generate a mock supernovae (SNe) data sample for a SNAP-like survey with a bump feature in w(z) hidden in, then successfully discover it by performing a blind wavelet analysis. We also apply our method to analyze the recently released "Constitution" SNe data, combined with WMAP and BAO from SDSS, and find weak hints of dark energy dynamics. Namely, we find that models with w(z) < -1 for 0.2 < z < 0.5, and w(z)> -1 for 0.5 < z <1, are mildly favored at 95% confidence level. This is in good agreement with several recent studies using other methods, such as redshift binning with principal component analysis (PCA) (e.g. Zhao and Zhang, arXiv:0908.1568)

Cosmic string networks generate cosmological perturbations actively throughout the history of the universe. Thus, the string sourced anisotropy of the cosmic microwave background is not affected by Silk damping as much as the anisotropy seeded by inflation. The spectrum of perturbations generated by strings does not match the observed CMB spectrum on large angular scales (l<1000) and is bounded to contribute no more than 10% of the total power on those scales. However, when this bound is marginally saturated, the anisotropy created by cosmic strings on small angular scales l>2000 will dominate over that created by the primary inflationary perturbations. This range of angular scales in the CMB is presently being measured by a number of experiments; their results will test this prediction of cosmic string networks soon.

Neutral hydrogen is ubiquitous, absorbing and emitting 21 cm radiation throughout much of the Universe's history. Active sources of perturbations, such as cosmic strings, would generate simultaneous perturbations in the distribution of neutral hydrogen and in the Cosmic Microwave Background (CMB) radiation from recombination. Moving strings would create wakes leading to 21 cm brightness fluctuations, while also perturbing CMB light via the Gott-Kaiser-Stebbins effect. This would lead to spatial correlations between the 21 cm and CMB anisotropies. Passive sources, like inflationary perturbations, predict no cross correlations prior to the onset of reionization. Thus, observation of any cross correlation between CMB and 21 cm radiation from dark ages would constitute evidence for new physics. We calculate the cosmic string induced correlations between CMB and 21 cm and evaluate their observability.

We study the evolution of linear perturbations in metric f(R) models of gravity and identify a potentially observable characteristic scale-dependent pattern in the behavior of cosmological structures. While at the background level viable f(R) models must closely mimic LCDM, the differences in their prediction for the growth of large scale structures can be sufficiently large to be seen with future weak lensing surveys. While working in the Jordan frame, we perform an analytical study of the growth of structures in the Einstein frame, demonstrating the equivalence of the dynamics in the two frames. We also provide a physical interpretation of the results in terms of the dynamics of an effective dark energy fluid with a non-zero shear. We find that the growth of structure in f(R) is enhanced, but that there are no small scale instabilities associated with the additional attractive "fifth force". We then briefly consider some recently proposed observational tests of modified gravity and their utility for detecting the f(R) pattern of structure growth.

Future weak lensing surveys will map the evolution of matter perturbations and gravitational potentials, yielding a new test of general relativity on cosmic scales. They will probe the relations between matter overdensities, local curvature, and the Newtonian potential. These relations can be modified in alternative gravity theories or by the effects of massive neutrinos or exotic dark energy fluids. We introduce two functions of time and scale which account for any such modifications in the linear regime. We use a principal component analysis to find the eigenmodes of these functions that cosmological data will constrain. The number of constrained modes gives a model-independent forecast of how many parameters describing deviations from general relativity could be constrained, along with $w(z)$. The modes' scale and time dependence tell us which theoretical models will be better tested.

We use a principal component approach to contrast different kinds of probes of dark energy, and to emphasize how an array of probes can work together to constrain an arbitrary equation of state history w(z). We pay particular attention to the role of the priors in assessing the information content of experiments and propose using an explicit prior on the degree of smoothness of w(z) that is independent of the binning scheme. We also show how a figure of merit based on the mean squared error probes the number of new modes constrained by a data set, and use it to examine how informative various experiments will be in constraining the evolution of dark energy.

In alternative theories of gravity, designed to produce cosmic acceleration at the current epoch, the growth of large scale structure can be modified. We study the potential of upcoming and future tomographic surveys such as DES and LSST, with the aid of CMB and supernovae data, to detect departures from the growth of cosmic structure expected within General Relativity. We employ parametric forms to quantify the potential time- and scale-dependent variation of the effective gravitational constant, and the differences between the two Newtonian potentials. We then apply the Fisher matrix technique to forecast the errors on the modified growth parameters from galaxy clustering, weak lensing, CMB, and their cross-correlations across multiple photometric redshift bins. We find that even with conservative assumptions about the data, DES will produce non-trivial constraints on modified growth, and that LSST will do significantly better.

We summarize the utility of precise cosmic microwave background (CMB) polarization measurements as probes of the physics of inflation. We focus on the prospects for using CMB measurements to differentiate various inflationary mechanisms. In particular, a detection of primordial B-mode polarization would demonstrate that inflation occurred at a very high energy scale, and that the inflaton traversed a super-Planckian distance in field space. We explain how such a detection or constraint would illuminate aspects of physics at the Planck scale. Moreover, CMB measurements can constrain the scale-dependence and non-Gaussianity of the primordial fluctuations and limit the possibility of a significant isocurvature contribution. Each such limit provides crucial information on the underlying inflationary dynamics. Finally, we quantify these considerations by presenting forecasts for the sensitivities of a future satellite experiment to the inflationary parameters.

How did the universe evolve? The fine angular scale (l>1000) temperature and polarization anisotropies in the CMB are a Rosetta stone for understanding the evolution of the universe. Through detailed measurements one may address everything from the physics of the birth of the universe to the history of star formation and the process by which galaxies formed. One may in addition track the evolution of the dark energy and discover the net neutrino mass. We are at the dawn of a new era in which hundreds of square degrees of sky can be mapped with arcminute resolution and sensitivities measured in microKelvin. Acquiring these data requires the use of special purpose telescopes such as the Atacama Cosmology Telescope (ACT), located in Chile, and the South Pole Telescope (SPT). These new telescopes are outfitted with a new generation of custom mm-wave kilo-pixel arrays. Additional instruments are in the planning stages.

Detecting the parity-odd, or B-mode, polarization pattern in the cosmic microwave background radiation due to primordial gravity waves is considered to be the final observational key to confirming the inflationary paradigm. The search for viable models of inflation from particle physics and string theory has (re)discovered another source for B-modes: cosmic strings. Strings naturally generate as much vector mode perturbation as they do scalar, producing B-mode polarization with a spectrum distinct from that expected from inflation itself. In a large set of models, B-modes arising from cosmic strings are more prominent than those expected from primordial gravity waves. In light of this, we study the physical underpinnings of string-sourced B-modes and the model dependence of the amplitude and shape of the $C_l^{BB}$ power spectrum. Observational detection of a string-sourced B-mode spectrum would be a direct probe of post-inflationary physics near the GUT scale. Conversely, non-detection would put an upper limit on a possible cosmic string tension of G\mu < 10^{-7} within the next three years.

We consider predictions for structure formation from modifications to general relativity in which the Einstein-Hilbert action is replaced by a general function of the Ricci scalar. We work without fixing a gauge, as well as in explicit popular coordinate choices, appropriate for the modification of existing cosmological code. We present the framework in a comprehensive and practical form that can be directly compared to standard perturbation analyses. By considering the full evolution equations, we resolve perceived instabilities previously suggested, and instead find a suppression of perturbations. This result presents significant challenges for agreement with current cosmological structure formation observations. The findings apply to a broad range of forms of f(R) for which the modification becomes important at low curvatures, disfavoring them in comparison with the LCDM scenario. As such, these results provide a powerful method to rule out a wide class of modified gravity models aimed at providing an alternative explanation to the dark energy problem.

Anthropic probability distributions for the cosmological constant and for the sum of neutrino masses are updated using the WMAP-3 data release. The new distribution for Lambda is in a better agreement with observation than the earlier one. The typicality of the observed value, defined as the combined probability of all values less likely than the observed, is no less than 22%. We discuss the dependence of our results on the simplifying assumptions used in deriving the distribution for Lambda and show that the agreement of the anthropic prediction with the data is rather robust. The distribution for the sum of the neutrino masses is peaked at 1 eV, suggesting degenerate masses, but a hierarchical mass pattern is still marginally allowed at a 2 sigma level.

Recent detections of the Integrated Sachs-Wolfe effect through the correlation of the cosmic microwave background temperature anisotropy with traces of large scale structure provided independent evidence for the expansion of the universe being dominated by something other than matter. Even with perfect data, statistical errors will limit the accuracy of such measurements to worse than 10%. On the other hand, the extraordinary sensitivity of the ISW effect to the details of structure formation should help to make up for the lack of precision. In these conference proceedings I discuss the extent to which future ISW measurements can help in testing the physics responsible for the observed cosmic acceleration.

We find the constraints from WMAP and SDSS data on the fraction of cosmological fluctuations sourced by local cosmic strings using a Markov Chain Monte Carlo (MCMC) analysis. In addition to varying the usual 6 cosmological parameters and the string tension ($\mu$), we also varied the amount of small-scale structure on the strings. Our results indicate that cosmic strings can account for up to 7 (14)% of the total power of the microwave anisotropy at 68 (95)% confidence level. The corresponding bound on the string mass per unit length, within our string model, is $G\mu < 1.8 (2.7) \times 10^{-7}$ at 68 (95)% c.l., where this constraint has been altered from what appears below following the correction of errors in our cosmic string code outlined in a recent erratum, astro-ph/0604141. We also calculate the B-type polarization spectra sourced by cosmic strings and discuss the prospects of their detection.

We investigate the effect of wiggly cosmic strings on the cosmic microwave background radiation anisotropy and matter power spectrum by modifying the string network model used by Albrecht et al.. We employ the wiggly equation of state for strings and the one-scale model for the cosmological evolution of certain network characteristics. For the same choice of simulation parameters we compare the results with and without including wiggliness in the model and find that wiggliness together with the accompanying low string velocities lead to a significant peak in the microwave background anisotropy and to an enhancement in the matter power spectrum. For the cosmologies we have investigated (standard CDM, and, CDM plus cosmological constant), and within the limitations of our modeling of the string network, the anisotropy is in reasonable agreement with current observations but the COBE normalized amplitude of density perturbations is lower than what the data suggests. In the case of a cosmological constant and CDM model, a bias factor of about 2 is required.

Overall, brane inflation is compatible with the recent analysis of the WMAP data. Here we explore the constraints of WMAP and 2dFGRS data on the various brane inflationary scenarios. Brane inflation naturally ends with the production of cosmic strings, which may provide a way to distinguish these models observationally. We argue that currently available data cannot exclude a non-negligible contribution from cosmic strings definitively. We perform a partial statistical analysis of mixed models that include a sub-dominant contribution from cosmic strings. Although the data favor models without cosmic strings, we conclude that they cannot definitively rule out a cosmic-string-induced contribution of $\sim 10 %$ to the observed temperature, polarization and galaxy density fluctuations. These results imply that $G\mu \lesssim 3.5\times 10^{-7}(\lambda/0.25)\sqrt{B/0.1}$, where $\lambda$ is a dimensionless parameter related to the interstring distance, and $B$ measures the importance of perturbations induced by cosmic strings. We argue that, conservatively, the data available currently still permit $B\lesssim 0.1$. Precision measurements sensitive to the B-mode polarization produced by vector density perturbation modes driven by the string network could provide evidence for these models. Accurate determinations of $n_s(k)$, the scalar fluctuation index, could also distinguish among various brane inflation models.

In a recent publication, we used the data from WMAP and SDSS to constrain the primordial perturbations and to predict the B-mode polarization sourced by cosmic string networks. We have been alerted by A. Slosar to the existence of errors in the code we used to calculate the Cosmic Microwave Background anisotropies from cosmic strings. Correcting the errors leads to a significant increase in the vector mode contribution to the CMB temperature and polarization anisotropies as well as an overall renormalization of the various string spectra. In these notes we explain the nature of the errors and discuss their implications for previously published constraints on cosmic strings based on this code. The chief change in our results is that our derived limit for the cosmic string tension is strengthened: G\mu < 1.8 (2.7) * 10^{-7} at 68 (95)% confidence. We also note that the newly-enhanced vector mode contribution produces a greatly-increased amplitude for B-mode polarization in the CMB which could exceed the B-mode power produced by the lensing of primordial E-mode polarization into B-mode polarization.

We present an analysis of the constraining power of future measurements of the Integrated Sachs-Wolfe (ISW) effect on models of the equation of state of dark energy as a function of redshift, w(z). To achieve this, we employ a new parameterization of w, which utilizes the mean value of w(z) (<w>) as an explicit parameter. This helps to separate the information contained in the estimation of the distance to the last scattering surface (from the CMB) from the information contained in the ISW effect. We then use Fisher analysis to forecast the expected uncertainties in the measured parameters from future ISW observations for two models of dark energy with very different time evolution properties. For example, we demonstrate that the cross-correlation of Planck CMB data and LSST galaxy catalogs will provide competitive constraints on w(z), compared to a SNAP-like SNe project, for models of dark energy with a rapidly changing equation of state (e.g. 'Kink' models). Our work confirms that, while SNe measurements are more suitable for constraining variations in w(z) at low redshift, the ISW effect can provide important independent constraints on w(z) at high z.

CMB power spectra and the SNIa luminosity-redshift curves have difficulty distinguishing between models with the same average value of the dark energy equation of state. We propose using the cross-correlation of the CMB temperature anisotropy with large scale structure to help break this degeneracy.

It is argued that small values of the neutrino masses may be due to anthropic selection effects. If this is the case, then the combined mass of the three neutrino species is expected to be ~1eV, neutrinos causing a non-negligible suppression of galaxy formation.

Current data exclude cosmic strings as the primary source of primordial density fluctuations. However, in a wide class of inflationary models, strings can form at later stages of inflation and have potentially detectable observational signatures. We study the constraints from WMAP and SDSS data on the fraction of primordial fluctuations sourced by local cosmic strings. The Bayesian analysis presented in this brief report is restricted to the minimal number of parameters. Yet it is useful for two reasons. It confirms the results of Pogosian et al (2003) using an alternative statistical method. Secondly, it justifies the more costly multi-parameter analysis. Already, varying only three parameters -- the spectral index and the amplitudes of the adiabatic and string contributions -- we find that the upper bound on the cosmic string contribution is of order 10%. We expect that the full multi-parameter study, currently underway, will likely loosen this bound.

Measurements of the CMB temperature anisotropy and the temperature-polarization cross correlation by WMAP suggest a reionization redshift of $z \sim 17 \pm 5$. On the other hand, observations of high redshift galaxies indicate a presence of a significant fraction of neutral hydrogen at redshift $z\sim 6-7$. We show that cosmic strings with tensions well within, but not far from, current observation bounds could cause early star formation at a level sufficient to explain the high reionization redshift.

It is argued that the observed vacuum energy density and the small values of the neutrino masses could be due to anthropic selection effects. Until now, these two quantities have been treated separately from each other and, in particular, anthropic predictions for the vacuum energy were made under the assumption of zero neutrino masses. Here we consider two cases. In the first, we calculate predictions for the vacuum energy for a fixed (generally non-zero) value of the neutrino mass. In the second we allow both quantities to vary from one part of the universe to another. We find that the anthropic predictions for the vacuum energy density are in a better agreement with observations when one allows for non-zero neutrino masses. We also find that the individual distributions for the vacuum energy and the neutrino masses are reasonably robust and do not change drastically when one adds the other variable.

One of the most tantalizing results from the WMAP experiment is the suggestion that the power at large scales is anomalously low when compared to the prediction of the ``standard'' Lambda-CDM model. The same anomaly, although with somewhat larger uncertainty, was also previously noted in the COBE data. In this work we discuss possible alternate models that give better fits on large scales and apply a model-comparison technique to select amongst them. We find that models with a cut off in the power spectrum at large scales are indeed preferred by data, but only by a factor of 3.6, at most, in the likelihood ratio, corresponding to about ``1.5 sigma'' if interpreted in the traditional manner. Using the same technique, we have also examined the possibility of a systematic error in the measurement (or prediction) of the large-scale power. Ignoring other evidence that the large-scale modes are properly measured and predicted, we find this possibility somewhat more likely, with roughly a 2.5 sigma evidence.

A broad class of dark energy models, which have been proposed in attempts at solving the cosmological constant problems, predict a late time variation of the equation of state with redshift. The variation occurs as a scalar field picks up speed on its way to negative values of the potential. The negative potential energy eventually turns the expansion into contraction and the local universe undergoes a big crunch. In this paper we show that cross-correlations of the CMB anisotropy and matter distribution, in combination with other cosmological data, can be used to forecast the imminence of such cosmic doomsday.

We examine the possibility that thermal, rather than quantum, fluctuations are responsible for seeding the structure of our universe. We find that while the thermalization condition leads to nearly Gaussian statistics, a Harrisson-Zeldovich spectrum for the primordial fluctuations can only be achieved in very special circumstances. These depend on whether the universe gets hotter or colder in time, while the modes are leaving the horizon. In the latter case we find a no-go theorem which can only be avoided if the fundamental degrees of freedom are not particle-like, such as in string gases near the Hagedorn phase transition. The former case is less forbidding, and we suggest two potentially successful ``warming universe'' scenarios. One makes use of the Phoenix universe, the other of ``phantom'' matter.

We have studied ways in which helical primordial magnetic fields could be constrained by measurements of the cosmic microwave background radiation (CMBR). If there were helical flows in the primordial plasma at the time of recombination, they would produce parity violating temperature-polarization correlations ($C_l^{TB}$ and $C_l^{EB}$). However, the magnitude of helical flows induced by helical magnetic fields is unobservably small. We discuss an alternate scheme for extracting the helicity of a stochastically homogeneous and isotropic primordial magnetic field using Faraday rotation measure maps of the CMBR and the power spectrum of B-type polarization ($C_l^{BB}$).

P and CP violation in cosmology can be manifested as large-scale helical velocity flows in the ambient plasma and as primordial helical magnetic fields. We show that kinetic helicity at last scattering leads to temperature-polarization correlations ($C_l^{TB}$ and $C_l^{EB}$) in the cosmic microwave background radiation (CMBR) and calculate the magnitude of the effect. Helical primordial magnetic fields, expected from cosmic events such as electroweak baryogenesis, can lead to helical velocity flows and hence to non-vanishing correlations of the temperature and B-type polarization. However we show that the magnitude of the induced helical flow is unobservably small because the helical component of a magnetic field is almost force-free. We discuss an alternate scheme for extracting the helicity of a stochastically homogeneous and isotropic primordial magnetic field using observations of the CMBR. The scheme involves constructing Faraday rotation measure maps of the CMBR and thus determining the sum of the helical and non-helical components of the primordial magnetic field. The power spectrum of B-type polarization fluctuations, on the other hand, are sensitive only to the non-helical component of the primordial magnetic field. The primordial magnetic helicity can then be derived by combining these two sets of observations.

We briefly review certain aspects of cosmic microwave background anisotropies as generated in passive and active models of structure formation. We then focus on cosmic strings based models and discuss their status in the light of current high-resolution observations from the BOOMERanG, MAXIMA and DASI collaborations. Upcoming megapixel experiments will have the potential to look for non-Gaussian features in the CMB temperature maps with unprecedented accuracy. We therefore devote the last part of this review to treat the non-Gaussianity of the microwave background and present a method for computation of the bispectrum from simulated string realizations.

We study dependence of the CMB temperature anisotropy spectrum on the value of the fine structure constant $\alpha$ and the equation of state of the dark energy component of the total density of the universe. We find that bounds imposed on the variation of $\alpha$ from the analysis of currently available CMB data sets can be significantly relaxed if one also allows for a change in the equation of state.

We propose a new method for a numerical computation of the angular bispectrum of the CMB anisotropies arising from active models such as cosmic topological defects, using a modified Boltzmann code based on CMBFAST. The method does not use CMB sky maps and requires moderate computational power. As a first implementation, we apply our method to a recently proposed model of simulated cosmic strings and estimate the observability of the non-Gaussian bispectrum signal. A comparison with the cosmic variance of the bispectrum estimator shows that the bispectrum for the simulated string model we used is not observable.

Recent measurements of the cosmic microwave background (CMB) anisotropies by BOOMERANG and MAXIMA collaborations have tightened the observational constraints on theories of structure formation. They disagree with the predictions of conventional topological defect models. Considering the fact that topological defects are predicted by the majority of realistic particle physics models, the exact nature of the constraints imposed by the recent data on the population and the properties of the defects must be fully understood. We show that the predictions of current cosmic string models can be brought into a closer agreement with the observations by choosing a closed universe with $\Omega=1.3$ and by including the effects of the small-scale structure and radiation products of the strings. These alone, however, are not sufficient for obtaining a good fit to the measured shape of the angular power spectrum. To fit the data cosmic strings would either have to be correlated on large (perhaps superhorizon) scales or would have to possess a higher degree of coherence, i.e. be more ``time-correlated'.

We consider the Grand Unified SU(5) model with a small or vanishing cubic term in the adjoint scalar field in the potential. This gives the model an approximate or exact Z$_2$ symmetry whose breaking leads to domain walls. The simplest domain wall has the structure of a kink across which the Higgs field changes sign ($\Phi \to -\Phi$) and inside which the full SU(5) is restored. The kink is shown to be perturbatively unstable for all parameters. We then construct a domain wall solution that is lighter than the kink and show it to be perturbatively stable for a range of parameters. The symmetry in the core of this domain wall is smaller than that outside. The interactions of the domain wall with magnetic monopole is discussed and it is shown that magnetic monopoles with certain internal space orientations relative to the wall pass through the domain wall. Magnetic monopoles in other relative internal space orientations are likely to be swept away on collision with the domain walls, suggesting a scenario where the domain walls might act like optical polarization filters, allowing certain monopole ``polarizations'' to pass through but not others. As SU(5) domain walls will also be formed at small values of the cubic coupling, this leads to a very complicated picture of the evolution of defects after the Grand Unified phase transition.

We study the interaction of magnetic monopoles and domain walls in a model with SU(5)\times Z_2 symmetry by numerically evolving the field equations. We find that the monopoles unwind and dissipate their magnetic energy on collision with domain walls within which the full SU(5) symmetry is restored.