The present work deals with holographic dark energy models with Hubble horizon as the infra-red cut-off. The interaction rate between dark energy and dark matter has been reconstructed with three different choices of the interaction term. It is shown that the coupling parameter of the interaction term should evolve with redshift to allow the successful transition from decelerated to accelerated phase of expansion. Constraints on the model parameters are obtained from Markov Chain Monte Carlo (MCMC) analysis using the supernova distance modulus data and observational measurements of the Hubble parameter. Results show that the model with the coupling parameter increasing with redshift (z) or equivalently decreasing with the evolution, are ruled out. On the other hand, coupling parameters, increasing or slowly varying with the evolution, are consistent with the observed evolution scenario. A Bayesian evidence calculation has been carried out for statistical selection of the reconstructed models. Though the kinematical parameters are well behaved for these models, the physical variables which determine the nature of the components in the matter sector, are not at all realistic. We have concluded that the existence of spatial curvature is essential for this particular type of dark energy models.

We consider the ghost-condensate model of dark energy with a generic potential term. The inclusion of the potential is shown to give greater freedom in realising the phantom regime. The self-consistency of the analysis is demonstrated using WMAP7+BAO+Hubble data.

In this paper we combine partial-wave (`modal') methods with a wavelet analysis of the CMB bispectrum. Our implementation exploits the advantages of both approaches to produce robust, reliable and efficient estimators which can constrain the amplitude of arbitrary primordial bispectra. This will be particularly important for upcoming surveys such as \emph{Planck}. A key advantage is the computational efficiency of calculating the inverse covariance matrix in wavelet space, producing an error bar which is close to optimal. We verify the efficacy and robustness of the method by applying it to WMAP7 data, finding $\fnllocal=38.4\pm 23.6$ and $\fnlequil=-119.2\pm123.6$.

The clustering of galaxies observed in future redshift surveys will provide a wealth of cosmological information. Matching the signal at different redshifts constrains the dark energy driving the acceleration of the expansion of the Universe. In tandem with these geometrical constraints, redshift-space distortions (RSD) depend on the build up of large-scale structure. As pointed out by many authors measurements of these effects are intrinsically coupled. We investigate this link, and argue that it strongly depends on the cosmological assumptions adopted when analysing data. Using representative assumptions for the parameters of the Euclid survey in order to provide a baseline future experiment, we show how the derived constraints change due to different model assumptions. We argue that even the assumption of a Friedman-Robertson-Walker (FRW) space-time is sufficient to reduce the importance of the coupling to a significant degree. Taking this idea further, we consider how the data would actually be analysed and argue that we should not expect to be able to simultaneously constrain multiple deviations from the standard $\Lambda$CDM model. We therefore consider different possible ways in which the Universe could deviate from the $\Lambda$CDM model, and show how the coupling between geometrical constraints and structure growth affects the measurement of such deviations.

We present an analysis of supernova light curves simulated for the upcoming Dark Energy Survey (DES) supernova search. The simulations employ a code suite that generates and fits realistic light curves in order to obtain distance modulus/redshift pairs that are passed to a cosmology fitter. We investigated several different survey strategies including field selection, supernova selection biases, and photometric redshift measurements. Using the results of this study, we chose a 30 square degree search area in the griz filter set. We forecast 1) that this survey will provide a homogeneous sample of up to 4000 Type Ia supernovae in the redshift range 0.05<z<1.2, and 2) that the increased red efficiency of the DES camera will significantly improve high-redshift color measurements. The redshift of each supernova with an identified host galaxy will be obtained from spectroscopic observations of the host. A supernova spectrum will be obtained for a subset of the sample, which will be utilized for control studies. In addition, we have investigated the use of combined photometric redshifts taking into account data from both the host and supernova. We have investigated and estimated the likely contamination from core-collapse supernovae based on photometric identification, and have found that a Type Ia supernova sample purity of up to 98% is obtainable given specific assumptions. Furthermore, we present systematic uncertainties due to sample purity, photometric calibration, dust extinction priors, filter-centroid shifts, and inter-calibration. We conclude by estimating the uncertainty on the cosmological parameters that will be measured from the DES supernova data.

We examine how dark energy constraints from current observational data depend on the analysis methods used: the analysis of Type Ia supernovae (SNe Ia), and that of galaxy clustering data. We generalize the flux-averaging analysis method of SNe Ia to allow correlated errors of SNe Ia, in order to reduce the systematic bias due to weak lensing of SNe Ia. We find that flux-averaging leads to larger errors on dark energy and cosmological parameters if only SN Ia data are used. When SN Ia data (the latest compilation by the SNLS team) are combined with WMAP 7 year results (in terms of our Gaussian fits to the probability distributions of the CMB shift parameters), the latest Hubble constant (H_0) measurement using the Hubble Space Telescope (HST), and gamma ray burst (GRB) data, flux-averaging of SNe Ia increases the concordance with other data, and leads to significantly tighter constraints on the dark energy density at z=1, and the cosmic curvature \Omega_k. The galaxy clustering measurements of H(z=0.35)r_s(z_d) and r_s(z_d)/D_A(z=0.35) (where H(z) is the Hubble parameter, D_A(z) is the angular diameter distance, and r_s(z_d) is the sound horizon at the drag epoch) by Chuang & Wang (2011) are consistent with SN Ia data, given the same pirors (CMB+H_0+GRB), and lead to significantly improved dark energy constraints when combined. Current data are fully consistent with a cosmological constant and a flat universe.

We demonstrate a methodology for optimizing the ability of future dark energy surveys to answer model selection questions, such as `Is acceleration due to a cosmological constant or a dynamical dark energy model?'. Model selection Figures of Merit are defined, exploiting the Bayes factor, and surveys optimized over their design parameter space via a Monte Carlo method. As a specific example we apply our methods to generic multi-fibre baryon acoustic oscillation spectroscopic surveys, comparable to that proposed for SuMIRe PFS, and present implementations based on the Savage-Dickey Density Ratio that are both accurate and practical for use in optimization. It is shown that whilst the optimal surveys using model selection agree with those found using the Dark Energy Task Force (DETF) Figure of Merit, they provide better informed flexibility of survey configuration and an absolute scale for performance; for example, we find survey configurations with close to optimal model selection performance despite their corresponding DETF Figure of Merit being at only 50% of its maximum. This Bayes factor approach allows us to interpret the survey configurations that will be good enough for the task at hand, vital especially when wanting to add extra science goals and in dealing with time restrictions or multiple probes within the same project.

COrE (Cosmic Origins Explorer) is a fourth-generation full-sky, microwave-band satellite recently proposed to ESA within Cosmic Vision 2015-2025. COrE will provide maps of the microwave sky in polarization and temperature in 15 frequency bands, ranging from 45 GHz to 795 GHz, with an angular resolution ranging from 23 arcmin (45 GHz) and 1.3 arcmin (795 GHz) and sensitivities roughly 10 to 30 times better than PLANCK (depending on the frequency channel). The COrE mission will lead to breakthrough science in a wide range of areas, ranging from primordial cosmology to galactic and extragalactic science. COrE is designed to detect the primordial gravitational waves generated during the epoch of cosmic inflation at more than $3\sigma $ for $r=(T/S)>=10^{-3}$. It will also measure the CMB gravitational lensing deflection power spectrum to the cosmic variance limit on all linear scales, allowing us to probe absolute neutrino masses better than laboratory experiments and down to plausible values suggested by the neutrino oscillation data. COrE will also search for primordial non-Gaussianity with significant improvements over Planck in its ability to constrain the shape (and amplitude) of non-Gaussianity. In the areas of galactic and extragalactic science, in its highest frequency channels COrE will provide maps of the galactic polarized dust emission allowing us to map the galactic magnetic field in areas of diffuse emission not otherwise accessible to probe the initial conditions for star formation. COrE will also map the galactic synchrotron emission thirty times better than PLANCK. This White Paper reviews the COrE science program, our simulations on foreground subtraction, and the proposed instrumental configuration.

The proposed CMBPol mission will be able to detect the imprint of topological defects on the cosmic microwave background (CMB) provided the contribution is sufficiently strong. We quantify the detection threshold for cosmic strings and for textures, and analyse the satellite's ability to distinguish between these different types of defects. We also assess the level of danger of misidentification of a defect signature as from the wrong defect type or as an effect of primordial gravitational waves. A 0.002 fractional contribution of cosmic strings to the CMB temperature spectrum at multipole ten, and similarly a 0.001 fractional contribution of textures, can be detected and correctly identified at the 3{\sigma} level. We also confirm that a tensor contribution of r = 0.0018 can be detected at over 3{\sigma}, in agreement with the CMBpol mission concept study. These results are supported by a model selection analysis.

We explore the possibility that a new physics interaction can provide an explanation for the knee just above $10^6$ GeV in the cosmic ray spectrum. We model the new physics modifications to the total proton-proton cross section with an incoherent term that allows for missing energy above the scale of new physics. We add the constraint that the new physics must also be consistent with published $pp$ cross section measurements, using cosmic ray observations, an order of magnitude and more above the knee. We find that the rise in cross section required at energies above the knee is radical. The increase in cross section suggests that it may be more appropriate to treat the scattering process in the black disc limit at such high energies. In this case there may be no clean separation between the standard model and new physics contributions to the total cross section. We model the missing energy in this limit and find a good fit to the Tibet III cosmic ray flux data. We comment on testing the new physics proposal for the cosmic ray knee at the Large Hadron Collider.

A space-based galaxy redshift survey would have enormous power in constraining dark energy and testing general relativity, provided that its parameters are suitably optimized. We study viable space-based galaxy redshift surveys, exploring the dependence of the Dark Energy Task Force (DETF) figure-of-merit (FoM) on redshift accuracy, redshift range, survey area, target selection, and forecast method. Fitting formulae are provided for convenience. We also consider the dependence on the information used: the full galaxy power spectrum P(k, P(k) marginalized over its shape, or just the Baryon Acoustic Oscillations (BAO). We find that the inclusion of growth rate information (extracted using redshift space distortion and galaxy clustering amplitude measurements) leads to a factor of ~ 3 improvement in the FoM, assuming general relativity is not modified. This inclusion partially compensates for the loss of information when only the BAO are used to give geometrical constraints, rather than using the full P(k) as a standard ruler. We find that a space-based galaxy redshift survey covering ~20,000 deg^2 over 0.5< z < 2 with \sigma_z/(1+z)<= 0.001 exploits a redshift range that is only easily accessible from space, extends to sufficiently low redshifts to allow both a vast 3-D map of the universe using a single tracer population, and overlaps with ground-based surveys to enable robust modeling of systematic effects. We argue that these parameters are close to their optimal values given current instrumental and practical constraints.

We present a Gaussianity analysis of the WMAP 5-year Cosmic Microwave Background (CMB) temperature anisotropy data maps. We use several third order estimators based on the spherical Mexican hat wavelet. We impose constraints on the local non-linear coupling parameter fnl using well motivated non-Gaussian simulations. We analyse the WMAP maps at resolution of 6.9 arcmin for the Q, V, and W frequency bands. We use the KQ75 mask recommended by the WMAP team which masks out 28% of the sky. The wavelet coefficients are evaluated at 10 different scales from 6.9 to 150 arcmin. With these coefficients we compute the third order estimators which are used to perform a chi-squared analysis. The chi-squared statistic is used to test the Gaussianity of the WMAP data as well as to constrain the fnl parameter. Our results indicate that the WMAP data are compatible with the Gaussian simulations, and the fnl parameter is constrained to -8 < fnl < +111 at 95% CL for the combined V+W map. This value has been corrected for the presence of undetected point sources, which add a positive contribution of Delta_fnl = 3+-5 in the V+W map. Our results are very similar to those obtained by Komatsu et al (2008) using the bispectrum.

We use Bayesian model selection tools to forecast the Planck satellite's ability to distinguish between different models for the reionization history of the Universe, using the large angular scale signal in the cosmic microwave background polarization spectrum. We find that Planck is not expected to be able to distinguish between an instantaneous reionization model and a two-parameter smooth reionization model, except for extreme values of the additional reionization parameter. If it cannot, then it will be unable to distinguish between different two-parameter models either. However, Bayesian model averaging will be needed to obtain unbiased estimates of the optical depth to reionization. We also generalize our results to a hypothetical future cosmic variance limited microwave anisotropy survey, where the outlook is more optimistic.

While observations indicate that the predominant source of cosmic inhomogeneities are adiabatic perturbations, there are a variety of candidates to provide auxiliary trace effects, including inflation-generated primordial tensors and cosmic defects which both produce B-mode cosmic microwave background (CMB) polarization. We investigate whether future experiments may suffer confusion as to the true origin of such effects, focusing on the ability of Planck to distinguish tensors from cosmic strings, and show that there is no significant degeneracy.

Although cosmic microwave background (CMB) anisotropy data alone cannot constrain simultaneously the spatial curvature and the equation of state of dark energy, CMB data provide a valuable addition to other experimental results. However computing a full CMB power spectrum with a Boltzmann code is quite slow; for instance if we want to work with many dark energy and/or modified gravity models, or would like to optimize experiments where many different configurations need to be tested, it is possible to adopt a quicker and more efficient approach. In this paper we consider the compression of the projected Planck CMB data into four parameters, R (scaled distance to last scattering surface), l_a (angular scale of sound horizon at last scattering), Omega_b h^2 (baryon density fraction) and n_s (powerlaw index of primordial matter power spectrum), all of which can be computed quickly. We show that, although this compression loses information compared to the full likelihood, such information loss becomes negligible when more data is added. We also demonstrate that the method can be used for scalar field dark energy independently of the parametrisation of the equation of state, and discuss how this method should be used for other kinds of dark energy models.

Current observational bounds on dark energy depend on assumptions about the curvature of the universe. We present a simple and efficient method for incorporating constraints from CMB anisotropy data, and use it to derive constraints on cosmic curvature and dark energy density as a free function of cosmic time using current data. We show that there are two CMB shift parameters, R=sqrt{\Omega_m H_0^2} r(z_{CMB}) (scaled distance to recombination) and l_a=\pi r(z_{CMB})/r_s(z_{CMB})(angular scale of the sound horizon at recombination), with measured values that are nearly uncorrelated with each other. Allowing nonzero cosmic curvature, the three-year WMAP data give R =1.71 +/- 0.03, l_a =302.5 +/- 1.2, and \Omega_b h^2 = 0.02173 +/- 0.00082, independent of the dark energy model. The corresponding bounds for a flat universe are R =1.70 +/- 0.03, l_a =302.2 +/- 1.2, and \Omega_b h^2 = 0.022 +/- 0.00082. We give the covariance matrix of (R, l_a, \Omega_b h^2) from the three-year WMAP data. We find that (R, l_a, \Omega_b h^2) provide an efficient and intuitive summary of CMB data as far as dark energy constraints are concerned. Using current CMB, SN Ia, and BAO data, we find that dark energy density is consistent with a constant in cosmic time, with marginal deviations from a cosmological constant that may reflect current systematic uncertainties or true evolution in dark energy. A flat universe is allowed by current data: \Omega_k=-0.006_{-0.012}^{+0.013}_{-0.025}^{+0.025} for w_X(z)=const., and \Omega_k=-0.002_{-0.018}^{+0.018}_{-0.032}^{+0.041} for w_X(z)=w_0+w_a(1-a)(68% and 95% C.L.). The bounds on cosmic curvature are less stringent if dark energy density is allowed to be a free function of cosmic time, and are also dependent on the assumption about the early time property of dark energy.

We use model selection forecasting to assess the ability of the Planck satellite to make a positive detection of spectral index running. We simulate Planck data for a range of assumed cosmological parameter values, and carry out a three-way Bayesian model comparison of a Harrison-Zel'dovich model, a power-law model, and a model including running. We find that Planck will be able to strongly support running only if its true value satisfies |dn/d ln k| > 0.02.

Modified gravity provides a possible explanation for the currently observed cosmic accelaration. In this paper, we study general classes of modified gravity models. The Einstein-Hilbert action is modified by using general functions of the Ricci and the Gauss-Bonnet scalars, both in the metric and in the Palatini formalisms. We do not use an explicit form for the functions, but a general form with a valid Taylor expansion up to second order about redshift zero in the Riemann-scalars. The coefficients of this expansion are then reconstructed via the cosmic expansion history measured using current cosmological observations. These are the quantities of interest for theoretical considerations relating to ghosts and instabilities. We find that current data provide interesting constraints on the coefficients. The next-generation dark energy surveys should shrink the allowed parameter space for modifed gravity models quite dramatically.

Probing the geometry of the universe is one of the most important endevours in cosmology. Current observational data from the Cosmic Microwave Background anisotropy (CMB), galaxy surveys and type Ia supernovae (SNe Ia) strongly constrain the curvature of the universe to be close to zero for a universe dominated by a cosmological constant or dark energy with a constant equation of state. Here we investigate the role of cosmic priors on deriving these tight bounds on geometry, by considering a landscape motivated scenario with an oscillating curvature term. We perform a likelihood analysis of current data under such a model of non-trivial geometry and find that the uncertainties on curvature, and correspondingly on parameters of the matter and dark energy sectors, are larger. Future dark energy experiments together with CMB data from experiments like Planck could dramatically improve our ability to constrain cosmic curvature under such models enabling us to probe possible imprints of quantum gravity.

We examine the choice of scale at which constraints on inflationary observables are presented. We describe an implementation of the hierarchy of inflationary consistency equations which ensures that they remain enforced on different scales, and then seek to optimize the scale for presentation of constraints on marginalized inflationary parameters from WMAP3 data. For models with spectral index running, we find a strong variation of the constraints through the range of observational scales available, and optimize by finding the scale which decorrelates constraints on the spectral index n_S and the running. This scale is k=0.017 Mpc^{-1}, and gives a reduction by a factor of more than four in the allowed parameter area in the n_S-r plane (r being the tensor-to-scalar ratio) relative to k=0.002 Mpc^{-1}. These optimized constraints are similar to those obtained in the no-running case. We also extend the analysis to a larger compilation of data, finding essentially the same conclusions.

In astro-ph/0702542v2, Linder and Miquel seek to criticize the use of Bayesian model selection for data analysis and for survey forecasting and design. Their discussion is based on three serious misunderstandings of the conceptual underpinnings and application of model-level Bayesian inference, which invalidate all their main conclusions. Their paper includes numerous further inaccuracies, including an erroneous calculation of the Bayesian Information Criterion. Here we seek to set the record straight.

We compute the Bayesian evidences for one- and two-parameter models of evolving dark energy, and compare them to the evidence for a cosmological constant, using current data from Type Ia supernova, baryon acoustic oscillations, and the cosmic microwave background. We use only distance information, ignoring dark energy perturbations. We find that, under various priors on the dark energy parameters, LambdaCDM is currently favoured as compared to the dark energy models. We consider the parameter constraints that arise under Bayesian model averaging, and discuss the implication of our results for future dark energy projects seeking to detect dark energy evolution. The model selection approach complements and extends the figure-of-merit approach of the Dark Energy Task Force in assessing future experiments, and suggests a significantly-modified interpretation of that statistic.

Model selection aims to determine which theoretical models are most plausible given some data, without necessarily asking about the preferred values of the model parameters. A common model selection question is to ask when new data require introduction of an additional parameter, describing a newly-discovered physical effect. We review several model selection statistics, and then focus on use of the Bayesian evidence, which implements the usual Bayesian analysis framework at the level of models rather than parameters. We describe our CosmoNest code, which is the first computationally-efficient implementation of Bayesian model selection in a cosmological context. We apply it to recent WMAP satellite data, examining the need for a perturbation spectral index differing from the scale-invariant (Harrison-Zel'dovich) case.

We use the three-year WMAP observations to determine the normalization of the matter power spectrum in inflationary cosmologies. In this context, the quantity of interest is not the normalization marginalized over all parameters, but rather the normalization as a function of the inflationary parameters n and r with marginalization over the remaining cosmological parameters. We compute this normalization and provide an accurate fitting function. The statistical uncertainty in the normalization is 3 percent, roughly half that achieved by COBE. We use the k-l relation for the standard cosmological model to identify the pivot scale for the WMAP normalization. We also quote the inflationary energy scale corresponding to the WMAP normalization.

Type Ia supernova (SN Ia), galaxy clustering, and cosmic microwave background anisotropy (CMB) data provide complementary constraints on the nature of the dark energy in the universe. We find that the three-year Wilkinson Microwave Anisotropy Probe (WMAP) observations give a CMB shift parameter of R = (\Omega_m H_0^2)^{1/2} \int_0^{z_{CMB}} dz'/H(z')= 1.70 \pm 0.03. Using this new measured value of the CMB shift parameter, together with the baryon acoustic oscillation (BAO) measurement from the Sloan Digital Sky Survey (SDSS), and SN Ia data from the HST/GOODS program and the first year Supernova Legacy Survey, we derive model-independent constraints on the dark energy density rho_X(z) and the cosmic expansion rate H(z). We also derive constraints on the dark energy equation of state w_X(z)=w_0+w'z (with cutoff at z=2) and w_X(a)=w_0+(1-a)w_a. We find that current data provide slightly tighter constraints on rho_X(z) and H(z) as free functions in redshift, and roughly a factor of two improvement in constraining w_X(z). A cosmological constant remains consistent with data, however, uncertainties remain large for model-independent constraints of dark energy. Significant increase in the number of observed SNe Ia between redshifts of 1 and 2, complemented by improved BAO and weak lensing cosmography measurements (as expected from the JEDI mission concept for the Joint Dark Energy Mission), will be required to dramatically tighten model-independent dark energy constraints.

The recent WMAP3 results have placed measurements of the spectral index n_S in an interesting position. While parameter estimation techniques indicate that the Harrison-Zel'dovich spectrum n_S=1 is strongly excluded (in the absence of tensor perturbations), Bayesian model selection techniques reveal that the case against n_S=1 is not yet conclusive. In this paper, we forecast the ability of the Planck satellite mission to use Bayesian model selection to convincingly exclude (or favour) the Harrison-Zel'dovich model.

We present a Bayesian model selection analysis of WMAP3 data using our code CosmoNest. We focus on the density perturbation spectral index $n_S$ and the tensor-to-scalar ratio $r$, which define the plane of slow-roll inflationary models. We find that while the Bayesian evidence supports the conclusion that $n_S \neq 1$, the data are not yet powerful enough to do so at a strong or decisive level. If tensors are assumed absent, the current odds are approximately 8 to 1 in favour of $n_S \neq 1$ under our assumptions, when WMAP3 data is used together with external data sets. WMAP3 data on its own is unable to distinguish between the two models. Further, inclusion of $r$ as a parameter weakens the conclusion against the Harrison-Zel'dovich case (n_S = 1, r=0), albeit in a prior-dependent way. In appendices we describe the CosmoNest code in detail, noting its ability to supply posterior samples as well as to accurately compute the Bayesian evidence. We make a first public release of CosmoNest, now available at http://www.cosmonest.org.

A key science goal of upcoming dark energy surveys is to seek time evolution of the dark energy. This problem is one of {\em model selection}, where the aim is to differentiate between cosmological models with different numbers of parameters. However, the power of these surveys is traditionally assessed by estimating their ability to constrain parameters, which is a different statistical problem. In this paper we use Bayesian model selection techniques, specifically forecasting of the Bayes factors, to compare the abilities of different proposed surveys in discovering dark energy evolution. We consider six experiments -- supernova luminosity measurements by the Supernova Legacy Survey, SNAP, JEDI, and ALPACA, and baryon acoustic oscillation measurements by WFMOS and JEDI -- and use Bayes factor plots to compare their statistical constraining power. The concept of Bayes factor forecasting has much broader applicability than dark energy surveys.

Several analyses of the microwave sky maps from the Wilkinson Microwave Anisotropy Probe (WMAP) have drawn attention to alignments amongst the low-order multipoles. Amongst the various possible explanations, an effect of cosmic topology has been invoked by several authors. We focus on an alignment of the first four multipoles (\ell = 2 to 5) found by Land and Magueijo (2005), and investigate the distribution of their alignment statistic for a set of simulated cosmic microwave background maps for cosmologies with slab-like topology. We find that this topology does offer a modest increase in the probability of the observed value, but that even for the smallest topology considered the probability of the observed value remains below one percent.

In order to reconstruct the initial conditions of the universe it is important to devise a method that can efficiently constrain the shape of the power spectrum of primordial matter density fluctuations in a model-independent way from data. In an earlier paper we proposed a method based on the wavelet expansion of the primordial power spectrum. The advantage of this method is that the orthogonality and multiresolution properties of wavelet basis functions enable information regarding the shape of $P_{\rm in}(k)$ to be encoded in a small number of non-zero coefficients. Any deviation from scale-invariance can then be easily picked out. Here we apply this method to simulated data to demonstrate that it can accurately reconstruct an input $P_{\rm in}(k)$, and present a prescription for how this method should be used on future data.

The abundance of new cosmological data becoming available means that a wider range of cosmological models are testable than ever before. However, an important distinction must be made between parameter fitting and model selection. While parameter fitting simply determines how well a model fits the data, model selection statistics, such as the Bayesian Evidence, are now necessary to choose between these different models, and in particular to assess the need for new parameters. We implement a new evidence algorithm known as nested sampling, which combines accuracy, generality of application and computational feasibility, and apply it to some cosmological datasets and models. We find that a five-parameter model with Harrison-Zel'dovich initial spectrum is currently preferred.

We study the statistical properties of the 1st year WMAP data on different scales using the spherical mexican hat wavelet transform. Consistent with the results of Vielva et al. (2003) we find a deviation from Gaussianity in the form of kurtosis of wavelet coefficients on $3-4^\circ$ scales in the southern Galactic hemisphere. This paper extends the work of Vielva et al. as follows. We find that the non-Gaussian signal shows up more strongly in the form of a larger than expected number of cold pixels and also in the form of scale-scale correlations amongst wavelet coefficients. We establish the robustness of the non-Gaussian signal under more wide-ranging assumptions regarding the Galactic mask applied to the data and the noise statistics. This signal is unlikely to be due to the usual quadratic term parametrized by the non-linearity parameter $f_{NL}$. We use the skewness of the spherical mexican hat wavelet coefficients to constrain $f_{NL}$ with the 1st year WMAP data. Our results constrain $f_{NL}$ to be $50\pm 80$ at 68% confidence, and less than 280 at 99% confidence.

A primordial cosmological magnetic field induces and supports vorticity or Alfven waves, which in turn generate cosmic microwave background (CMB) anisotropies. A homogeneous primordial magnetic field with fixed direction induces correlations between the $a_{l-1,m}$ and $a_{l+1,m}$ multipole coefficients of the CMB temperature anisotropy field. We discuss the constraints that can be placed on the strength of such a primordial magnetic field using CMB anisotropy data from the WMAP experiment. We place 3 $\sigma$ upper limits on the strength of the magnetic field of $B < 15$ nG for vector perturbation spectral index $n=-5$ and $B<1.7$ nG for $n=-7$.

We reconstruct the dark energy density $\rho_X(z)$ as a free function from current type Ia supernova (SN Ia) data (Tonry et al. 2003; Barris et al. 2003; Knop et al. 2003), together with the Cosmic Microwave Background (CMB) shift parameter from CMB data (WMAP, CBI, and ACBAR), and the large scale structure (LSS) growth factor from 2dF galaxy survey data. We parametrize $\rho_X(z)$ as a continuous function, given by interpolating its amplitudes at equally spaced $z$ values in the redshift range covered by SN Ia data, and a constant at larger $z$ (where $\rho_X(z)$ is only weakly constrained by CMB data). We assume a flat universe, and use the Markov Chain Monte Carlo (MCMC) technique in our analysis. We find that the dark energy density $\rho_X(z)$ is constant for $0 \la z \la 0.5$ and increases with redshift $z$ for $0.5 \la z \la 1$ at 68.3% confidence level, but is consistent with a constant at 95% confidence level. For comparison, we also give constraints on a constant equation of state for the dark energy. Flux-averaging of SN Ia data is required to yield cosmological parameter constraints that are free of the bias induced by weak gravitational lensing \citep{Wang00b}. We set up a consistent framework for flux-averaging analysis of SN Ia data, based on \cite{Wang00b}. We find that flux-averaging of SN Ia data leads to slightly lower $\Omega_m$ and smaller time-variation in $\rho_X(z)$. This suggests that a significant increase in the number of SNe Ia from deep SN surveys on a dedicated telescope \citep{Wang00a} is needed to place a robust constraint on the time-dependence of the dark energy density.

In order to constrain and possibly detect unusual physics during inflation, we allow the power spectrum of primordial matter density fluctuations, P_{in}(k), to be an arbitrary function in the estimation of cosmological parameters from data. The multi-resolution and good localization properties of orthogonal wavelets make them suitable for detecting features in P_{in}(k). We expand P_{in}(k) directly in wavelet basis functions. The likelihood of the data is thus a function of the wavelet coefficients of P_{in}(k), as well as the H_0, \Omega_b h^2, $\Omega_c h^2 and the \tau_{ri}, in a flat $\Lambda$CDM cosmology. We derive constraints on these parameters from CMB anisotropy data (WMAP, CBI, and ACBAR) and large scale structure (LSS) data (2dFGRS and PSCZ) using the Markov Chain Monte Carlo (MCMC) technique. The direct wavelet expansion method is different and complimentary to the wavelet band power method of Mukherjee & Wang (2003a,b), and results from the two methods are consistent. In addition, as we demonstrate, the direct wavelet expansion method has the advantage that once the wavelet coefficients have been constrained, the reconstruction of P_{in}(k) can be effectively denoised, i.e., P_{in}(k) can be reconstructed using only the coefficients that, say, deviate from zero at greater than 1\sigma. In doing so the essential properties of P_{in}(k) are retained. The reconstruction also suffers much less from the correlated errors of binning methods. The shape of the primordial power spectrum, as reconstructed in detail here, reveals an interesting new feature at 0.001 \la k/{Mpc}^{-1} \la 0.005. It will be interesting to see if this feature is confirmed by future data. The reconstructed and denoised P_{in}(k) is favored over the scale-invariant and power-law forms at \ga 1\sigma.

Reconstructing the shape of the primordial power spectrum in a model independent way from cosmological data is a useful consistency check on what is usually assumed regarding early universe physics. It is also our primary window to unknown physics during the inflationary era. Using a power-law form for the primordial power spectrum $P_{in}(k)$ and constraining the scalar spectral index and its running, \cite{Peiris03} found that the first year WMAP data seem to indicate a preferred scale in $P_{in}(k)$. We use two complementary methods: the wavelet band powers method of \cite{pia1}, and the top-hat binning method of \cite{Wang99} to reconstruct $P_{in}(k)$ as a free function from CMB data alone (WMAP, CBI, and ACBAR), or from CMB data together with large scale structure data (2dFGRS and PCSZ). The shape of the reconstructed $P_{in}(k)$ is consistent with scale-invariance, although it allows some indication of a preferred scale at $k \sim 0.01$Mpc$^{-1}$. While consistent with the possible evidence for a running of the scalar spectral index found by the WMAP team, our results highlight the need of more stringent and independent constraints on cosmological parameters (the Hubble constant in particular) in order to more definitively constrain deviations of $P_{in}(k)$ from scale-invariance without making assumptions about the inflationary model.

Likelihood analyses of the COBE-DMR sky maps are used to determine the normalization of the inverse-power-law-potential scalar field dark energy model. Predictions of the DMR-normalized model are compared to various observations to constrain the allowed range of model parameters. Although the derived constraints are restrictive, evolving dark energy density scalar field models remain an observationally-viable alternative to the constant cosmological constant model.

We constrain Galactic foreground contamination of the Python V cosmic microwave background anisotropy data by cross correlating it with foreground contaminant emission templates. To model foreground emission we use 100 and 12 $\mu$m dust emission templates and two point source templates based on the PMN survey. The analysis takes account of inter-modulation correlations in 8 modulations of the data that are sensitive to a large range of angular scales and also densely sample a large area of sky. As a consequence the analysis here is highly constraining. We find little evidence for foreground contamination in an analysis of the whole data set. However, there is indication that foregrounds are present in the data from the larger-angular-scale modulations of those Python V fields that overlap the region scanned earlier by the UCSB South Pole 1994 experiment. This is an independent consistency cross-check of findings from the South Pole 1994 data.

Measuring the primordial matter power spectrum is our primary means of probing unknown physics in the very early universe. We allow the primordial power spectrum to be an arbitrary function, and parametrize it in terms of wavelet band powers. Current cosmological data correspond to 11 such wavelet bands. We derive constraints on these band powers as well as $H_0$, $\Omega_b h^2$ and $\Omega_m h^2$ from current Cosmic Microwave Background Anisotropy (CMB) data using the Markov Chain Monte Carlo (MCMC) technique. Our results indicate a feature in the primordial power spectrum at $ 0.008 \la k/(h {Mpc}^{-1}) \la 0.1$. MAP and Planck data should allow us to put tighter constraints on the primordial power spectrum.

We use joint likelihood analyses of combinations of fifteen cosmic microwave background (CMB) anisotropy data sets from the DMR, UCSB South Pole 1994, Python I--III, ARGO, MAX 4 and 5, White Dish, OVRO, and SuZIE experiments to constrain cosmogonies. We consider open and spatially-flat-Lambda cold dark matter cosmogonies, with nonrelativistic-mass density parameter Omega_0 in the range 0.1--1, baryonic-mass density parameter Omega_B in the range (0.005--0.029) h^{-2}, and age of the universe t_0 in the range (10--20) Gyr. Marginalizing over all parameters but Omega_0, the data favor Omega_0 \simeq 0.9--1 (0.4--0.6) flat-Lambda (open) models. The range in deduced Omega_0 values is partially a consequence of the different combinations of smaller-angular-scale CMB anisotropy data sets used in the analyses, but more significantly a consequence of whether the DMR quadrupole moment is accounted for or ignored in the analysis. For both flat-Lambda and open models, after marginalizing over all other parameters, a lower Omega_B h^2 \simeq 0.005--0.009 is favored. This is also marginally at odds with estimates from more recent CMB anisotropy data and some estimates from standard nucleosynthesis theory and observed light element abundances. For both sets of models a younger universe with t_0 \simeq 12--15 Gyr is favored, consistent with other recent non-CMB indicators. We emphasize that since we consider only a small number of data sets, these results are tentative. More importantly, the analyses here do not rule out the currently favored flat-Lambda model with Omega_0 \sim 0.3, nor the larger Omega_B h^2 values favored by some other data.

We use Owens Valley Radio Observatory (OVRO) cosmic microwave background (CMB) anisotropy data to constrain cosmological parameters. We account for the OVRO beamwidth and calibration uncertainties, as well as the uncertainty induced by the removal of non-CMB foreground contamination. We consider open and spatially-flat-Lambda cold dark matter cosmogonies, with nonrelativistic-mass density parameter Omega_0 in the range 0.1--1, baryonic-mass density parameter Omega_B in the range (0.005--0.029)h^{-2}, and age of the universe t_0 in the range (10--20) Gyr. Marginalizing over all parameters but Omega_0, the OVRO data favors an open (spatially-flat-Lambda) model with Omega_0 ~= 0.33 (0.1). At the 2 sigma confidence level model normalizations deduced from the OVRO data are mostly consistent with those deduced from the DMR, UCSB South Pole 1994, Python I-III, ARGO, MAX 4 and 5, White Dish, and SuZIE data sets.

We study Galactic emission foreground contamination of the OVRO and UCSB South Pole 1994 cosmic microwave background (CMB) anisotropy data by cross-correlating with templates of infrared dust emission and new high resolution $H_{\alpha}$ data. $H_{\alpha}$ data provide rough upper limits on the level of free-free emission in the data sets. The cross-correlation analysis does not contradict a two component foreground emission hypothesis, with the two dust-correlated components being free-free emission and spinning dust emission.

Recently correlation analyses between Galactic dust emission templates and a number of CMB data sets have led to differing claims on the origin of the Galactic contamination at low frequencies. de Oliviera-Costa {\em et al} (1999) have presented work based on Tenerife data supporting the spinning dust hypothesis. Since the frequency coverage of these data is ideal to discriminate spectrally between spinning dust and free-free emission, we used the latest version of the Tenerife data, which have lower systematic uncertainty, to study the correlation in more detail. We found however that the evidence in favor of spinning dust originates from a small region at low Galactic latitude where the significance of the correlation itself is low and is compromised by systematic effects in the Galactic plane signal. The rest of the region was found to be uncorrelated. Regions that correlate with higher significance tend to have a steeper spectrum, as is expected for free-free emission. Averaging over all correlated regions yields dust-correlation coefficients of $180\pm47$ and $123\pm16$ $\mu$K /MJy sr$^{-1}$ at 10 and 15 GHz respectively. These numbers however have large systematic uncertainties that we have identified and care should be taken when comparing with results from other experiments. We do find evidence for synchrotron emission with spectral index steepening from radio to microwave frequencies, but we cannot make conclusive claims about the origin of the dust-correlated component based on the spectral index estimates. Data with higher sensitivity are required to decide about the significance of the dust-correlation at high Galactic latitudes and other Galactic templates, in particular H$_\alpha$ maps, will be necessary for constraining its origin.

We investigate the detection of non-Gaussianity in the 4-year COBE data reported by Pando, Valls-Gabaud & Fang (1998), using a technique based on the discrete wavelet transform. Their analysis was performed on the two DMR faces centred on the North and South Galactic poles respectively, using the Daubechies 4 wavelet basis. We show that these results depend critically on the orientation of the data, and so should be treated with caution. For two distinct orientations of the data, we calculate unbiased estimates of the skewness, kurtosis and scale-scale correlation of the corresponding wavelet coefficients in all of the available scale domains of the transform. We obtain several detections of non-Gaussianity in the DMR-DSMB map at greater than the 99 per cent confidence level, but most of these occur on pixel-pixel scales and are therefore not cosmological in origin. Indeed, after removing all multipoles beyond $\ell = 40$ from the COBE maps, only one robust detection remains. Moreover, using Monte-Carlo simulations, we find that the probability of obtaining such a detection by chance is 0.59. We repeat the analysis for the 53+90 GHz coadded COBE map. In this case, after removing $\ell > 40$ multipoles, two non-Gaussian detections at the 99 per cent level remain. Nevertheless, again using Monte-Carlo simulations, we find that the probability of obtaining two such detections by chance is 0.28. Thus, we conclude the wavelet technique does {\em not} yield strong evidence for non-Gaussianity of cosmological origin in the 4-year COBE data.

We use data from the Tenerife 10, 15 and 33 GHz beamswitching experiments along with the COBE 53 and 90 GHz data to separate the cosmic microwave background (CMB) signal from the Galactic signal and create two maps at high Galactic latitude. The new multi-MEM technique is used to obtain the best reconstruction of the two channels. The two maps are presented and known features are identified within each. We find that the Galactic contribution to both the 15 and 33 GHz Tenerife data is small enough to be ignored when compared to the errors in the data and the magnitude of the CMB signal.

We study the effect of Galactic foregrounds with spatially varying spectral indices on analysis of simulated data from the Planck Satellite. We also briefly mention the effect the extra-galactic point sources have on the data analysis and summarise the most recent constraints on Galactic emission at GHz frequencies.