Source author record

Martin Kunz

Martin Kunz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

67works
16topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

67 published item(s)

preprint2022arXiv

A Cautionary Tale: Dark Energy in Single-Field, Slow-Roll Inflationary Models

The current epoch of accelerated cosmic expansion is postulated to be driven by dark energy, which in the standard model takes the form of a cosmological constant with equation of state parameter $w=-1$. We propose an innovative perspective over the nature of dark energy by drawing a parallel with inflation, which we assume to be driven by a single scalar field, the inflaton. The inflaton was not a cosmological constant, as indicated by the fact that inflation ended and by the Planck satellite's constraint of $n_s\neq 1$ at $8σ$ confidence. Therefore, it is interesting to verify whether its equation of state parameter was measurably different from $-1$. We analyze this question for a class of single-field slow-roll inflationary models, where the hierarchy of Hubble slow-roll parameters is truncated at different orders. Based on the latest Planck and BICEP2/Keck data, we obtain a $68\%$ upper bound of $1+w < 0.0014$ for the three-parameter model, which gives the best description to the data. This provides a cautionary tale for drawing conclusions about the nature of today's dark energy based upon the non-detection of a deviation from $w=-1$ with current and upcoming cosmological surveys.

preprint2022arXiv

An Estimator for the lensing potential from galaxy number counts

We derive an estimator for the lensing potential from galaxy number counts which contains a linear and a quadratic term. We show that this estimator has a much larger signal-to-noise ratio than the corresponding estimator from intensity mapping. We show that this is due to the additional lensing term in the number count angular power spectrum which is present already at linear order. We estimate the signal-to-noise ratio for future photometric surveys. We find that particularly at high redshifts, $z\gtrsim 1.5$, the signal to noise ratio can become of order 30. We therefore claim that number counts in photometric surveys are an excellent means to measure tomographic lensing spectra.

preprint2022arXiv

The Hydrogen Intensity and Real-time Analysis eXperiment: 256-Element Array Status and Overview

The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a radio interferometer array currently in development, with an initial 256-element array to be deployed at the South African Radio Astronomy Observatory (SARAO) Square Kilometer Array (SKA) site in South Africa. Each of the 6m, $f/0.23$ dishes will be instrumented with dual-polarisation feeds operating over a frequency range of 400-800 MHz. Through intensity mapping of the 21 cm emission line of neutral hydrogen, HIRAX will provide a cosmological survey of the distribution of large-scale structure over the redshift range of $0.775 < z < 2.55$ over $\sim$15,000 square degrees of the southern sky. The statistical power of such a survey is sufficient to produce $\sim$7 percent constraints on the dark energy equation of state parameter when combined with measurements from the Planck satellite. Additionally, HIRAX will provide a highly competitive platform for radio transient and HI absorber science while enabling a multitude of cross-correlation studies. In this paper, we describe the science goals of the experiment, overview of the design and status of the sub-components of the telescope system, and describe the expected performance of the initial 256-element array as well as the planned future expansion to the final, 1024-element array.

preprint2021arXiv

A new estimator for gravitational lensing using galaxy and intensity mapping surveys

We introduce the Galaxy Intensity Mapping cross-COrrelation estimator (GIMCO), which is a new tomographic estimator for the gravitational lensing potential, based on a combination of intensity mapping (IM) and galaxy number counts. The estimator can be written schematically as IM$(z_f)\times$galaxy$(z_b)$ $-$ galaxy$(z_f)\times$IM$(z_b)$ for a pair of distinct redshifts $(z_f,z_b)$; this combination allows to greatly reduce the contamination by density-density correlations, thus isolating the lensing signal. As an estimator constructed only from cross-correlations, it is additionally less susceptible to systematic effects. We show that the new estimator strongly suppresses cosmic variance and consequently improves the signal-to-noise ratio (SNR) for the detection of lensing, especially on linear scales and intermediate redshifts. %This makes it particularly valuable for future studies of dark energy and modified gravity. For cosmic variance dominated surveys, the SNR of our estimator is a factor 30 larger than the SNR obtained from the correlation of galaxy number counts only. Shot noise and interferometer noise reduce the SNR. For the specific example of the Dark Energy Survey (DES) cross-correlated with the Hydrogen Intensity mapping and Real time Analysis eXperiment (HIRAX), the SNR is around 4, whereas for Euclid cross-correlated with HIRAX it reaches 52. This corresponds to an improvement of a factor 4-5 compared to the SNR from DES alone. For Euclid cross-correlated with HIRAX the improvement with respect to Euclid alone strongly depends on the redshift. We find that the improvement is particularly important for redshifts below 1.6, where it reaches a factor of 5. This makes our estimator especially valuable to test dark energy and modified gravity, that are expected to leave an impact at low and intermediate redshifts.

preprint2020arXiv

Dark degeneracy I: Dynamical or interacting dark energy?

We revisit the dark degeneracy that arises from the Einstein equations relating geometry to the total cosmic substratum but not resolving its individual components separately. We establish the explicit conditions for the dark degeneracy in the fluid description of the dark sector. At the background level, this degeneracy can be formally understood in terms of a unified dark sector Equation of State (EoS) that depends both on the dynamical nature of the dark energy (DE) as well as on its interaction with the pressureless dark matter. For linear perturbations, the degeneracy arises for specified DE pressure perturbations (or sound speed, equivalently) and DE anisotropic stress. Specializing to the degeneracy between non-interacting dynamical DE and interacting vacuum DE models, we perform a parameter estimation analysis for a range of dynamical DE parametrizations, where for illustration we explicitly break the degeneracy at the linear level by adopting a luminal sound speed for both scenarios. We conduct this analysis using cosmological background data alone and in combination with Planck measurements of the cosmic microwave background radiation. We find that although the overall phenomenology between the dynamical DE and interacting approaches is similar, there are some intriguing differences. In particular, there is an ambiguity in the strength of constraints on $Ω_{m0}$ and $σ_8$, which are considerably weakened for interacting vacuum DE, indicating that the dark degeneracy can change the significance of tensions in cosmological parameters inferred from different data sets.

preprint2020arXiv

Distinguishing cosmologies using the turn-around radius near galaxy clusters

Outside galaxy clusters the competition between the inwards gravitational attraction and the outwards expansion of the Universe leads to a special radius of velocity cancellation, which is called the turn-around radius. Measurements of the turn-around radius hold promises of constraining cosmological parameters, and possibly even properties of gravity. Such a measurement is, however, complicated by the fact that the surroundings of galaxy clusters are not spherical, but instead are a complicated collection of filaments, sheets and voids. In this paper we use the results of numerically simulated universes to quantify realistic error-bars of the measurement of the turn-around radius. We find that for a $Λ$CDM cosmology these error-bars are typically of the order of $20\%$. We numerically simulate three different implementations of dark energy models and of a scalar dark sector interaction to address whether the turn-around radius can be used to constrain non-trivial cosmologies, and we find that only rather extreme models can be distinguished from a $Λ$CDM universe due to the large error-bars arising from the non-trivial cluster environments.

preprint2019arXiv

A Flexible Framework for Anomaly Detection via Dimensionality Reduction

Anomaly detection is challenging, especially for large datasets in high dimensions. Here we explore a general anomaly detection framework based on dimensionality reduction and unsupervised clustering. We release DRAMA, a general python package that implements the general framework with a wide range of built-in options. We test DRAMA on a wide variety of simulated and real datasets, in up to 3000 dimensions, and find it robust and highly competitive with commonly-used anomaly detection algorithms, especially in high dimensions. The flexibility of the DRAMA framework allows for significant optimization once some examples of anomalies are available, making it ideal for online anomaly detection, active learning and highly unbalanced datasets.

preprint2019arXiv

Non-linear contributions to angular power spectra

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

preprint2019arXiv

Nucleation of Dislocations in 3.9 nm Nanocrystals at High Pressure

As circuitry approaches single nanometer length scales, it is important to predict the stability of metals at these scales. The behavior of metals at larger scales can be predicted based on the behavior of dislocations, but it is unclear if dislocations can form and be sustained at single nanometer dimensions. Here, we report the formation of dislocations within individual 3.9 nm Au nanocrystals under nonhydrostatic pressure in a diamond anvil cell. We used a combination of x-ray diffraction, optical absorbance spectroscopy, and molecular dynamics simulation to characterize the defects that are formed, which were found to be surface-nucleated partial dislocations. These results indicate that dislocations are still active at single nanometer length scales and can lead to permanent plasticity.

preprint2016arXiv

Constraints on dark-matter properties from large-scale structure

We use large-scale cosmological observations to place constraints on the dark-matter pressure, sound speed and viscosity, and infer a limit on the mass of warm-dark-matter particles. Measurements of the cosmic microwave background (CMB) anisotropies constrain the equation of state and sound speed of the dark matter at last scattering at the per mille level. Since the redshifting of collisionless particles universally implies that these quantities scale like $a^{-2}$ absent shell crossing, we infer that today $w_{\rm (DM)}< 10^{-10.0}$, $c_{\rm s,(DM)}^2 < 10^{-10.7}$ and $c_{\rm vis, (DM)}^{2} < 10^{-10.3}$ at the $99\%$ confidence level. This very general bound can be translated to model-dependent constraints on dark-matter models: for warm dark matter these constraints imply $m> 70$ eV, assuming it decoupled while relativistic around the same time as the neutrinos; for a cold relic, we show that $m>100$ eV. We separately constrain the properties of the DM fluid on linear scales at late times, and find upper bounds $c_{\rm s, (DM)}^2<10^{-5.9}$, $c_{\rm vis, (DM)}^{2} < 10^{-5.7}$, with no detection of non-dust properties for the DM.

preprint2016arXiv

Cosmic Infrared Background anisotropies as a window into primordial non-Gaussianity

The angular power spectrum of the cosmic infrared background (CIB) is a sensitive probe of the local primordial bispectrum. CIB measurements are integrated over a large volume so that the scale dependent bias from the primordial non-Gaussianity leaves a strong signal in the CIB power spectrum. Although galactic dust dominates over the non-Gaussian CIB signal, it is possible to mitigate the dust contamination with enough frequency channels, especially if high frequencies such as the Planck 857 GHz channel are available. We show that, in this case, measurements of the cosmic microwave background from future space missions should be able to probe the local bispectrum shape down to an amplitude |f_nl| < 1.

preprint2016arXiv

gevolution: a cosmological N-body code based on General Relativity

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

preprint2016arXiv

Lensing convergence and the neutrino mass scale in galaxy redshift surveys

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

preprint2016arXiv

Non-local gravity and comparison with observational datasets. II. Updated results and Bayesian model comparison with $Λ$CDM

We present a comprehensive and updated comparison with cosmological observations of two non-local modifications of gravity previously introduced by our group, the so called RR and RT models. We implement the background evolution and the cosmological perturbations of the models in a modified Boltzmann code, using CLASS. We then test the non-local models against the {\em Planck} 2015 TT, TE, EE and Cosmic Microwave Background (CMB) lensing data, isotropic and anisotropic Baryonic Acoustic Oscillations (BAO) data, JLA supernovae, $H_0$ measurements and growth rate data, and we perform Bayesian parameter estimation. We then compare the RR, RT and $Λ$CDM models, using the Savage-Dickey method. We find that the RT model and $Λ$CDM perform equally well, while the performance of the RR model with respect to $Λ$CDM depends on whether or not we include a prior on $H_0$ based on local measurements.

preprint2015arXiv

Bayesian Inference for Radio Observations

New telescopes like the Square Kilometre Array (SKA) will push into a new sensitivity regime and expose systematics, such as direction-dependent effects, that could previously be ignored. Current methods for handling such systematics rely on alternating best estimates of instrumental calibration and models of the underlying sky, which can lead to inadequate uncertainty estimates and biased results because any correlations between parameters are ignored. These deconvolution algorithms produce a single image that is assumed to be a true representation of the sky, when in fact it is just one realization of an infinite ensemble of images compatible with the noise in the data. In contrast, here we report a Bayesian formalism that simultaneously infers both systematics and science. Our technique, Bayesian Inference for Radio Observations (BIRO), determines all parameters directly from the raw data, bypassing image-making entirely, by sampling from the joint posterior probability distribution. This enables it to derive both correlations and accurate uncertainties, making use of the flexible software MEQTREES to model the sky and telescope simultaneously. We demonstrate BIRO with two simulated sets of Westerbork Synthesis Radio Telescope data sets. In the first, we perform joint estimates of 103 scientific (flux densities of sources) and instrumental (pointing errors, beamwidth and noise) parameters. In the second example, we perform source separation with BIRO. Using the Bayesian evidence, we can accurately select between a single point source, two point sources and an extended Gaussian source, allowing for 'super-resolution' on scales much smaller than the synthesized beam.

preprint2015arXiv

Can dark energy viscosity be detected with the Euclid survey?

Recent work has demonstrated that it is important to constrain the dynamics of cosmological perturbations, in addition to the evolution of the background, if we want to distinguish among different models of the dark sector. Especially the anisotropic stress of the (possibly effective) dark energy fluid has been shown to be an important discriminator between modified gravity and dark energy models. In this paper we use approximate analytical solutions of the perturbation equations in the presence of viscosity to study how the anisotropic stress affects the weak lensing and galaxy power spectrum. We then forecast how sensitive the photometric and spectroscopic Euclid surveys will be to both the speed of sound and the viscosity of our effective dark energy fluid when using weak lensing tomography and the galaxy power spectrum. We find that Euclid alone can only constrain models with very small speed of sound and viscosity, while it will need the help of other observables in order to give interesting constraints on models with a sound speed close to one. This conclusion is also supported by the expected Bayes factor between models.

preprint2015arXiv

Cosmology and fundamental physics with the Euclid satellite

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

preprint2015arXiv

Does small scale structure significantly affect cosmological dynamics?

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

preprint2015arXiv

General relativity and cosmic structure formation

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

preprint2015arXiv

Semi-analytical Study on the Generic Degeneracy for Galaxy Clustering Measurements

From the galaxy power spectrum in redshift space, we derive semi-analytical results on the generic degeneracy of galaxy clustering measurements. Defining the observables $\bar{A}= Gbσ_8$ and $\bar{R} = Gfσ_8$, (being $G$ the growth function, $b$ the bias, $f$ the growth rate, and $σ_8$ the amplitude of the power spectrum), we perform a Fisher matrix formalism to forecast the expected precision of these quantities for a Euclid-like survey. Among the results we found that galaxy surveys have generically a slightly negative correlation between $\bar{A}$ and $\bar{R}$, and they can always measure $\bar{R}$ about 3.7 to 4.7 times better than $\bar{A}$.

preprint2015arXiv

Testing foundations of modern cosmology with SKA all-sky surveys

Continuum and HI surveys with the Square Kilometre Array (SKA) will allow us to probe some of the most fundamental assumptions of modern cosmology, including the Cosmological Principle. SKA all-sky surveys will map an enormous slice of space-time and reveal cosmology at superhorizon scales and redshifts of order unity. We illustrate the potential of these surveys and discuss the prospects to measure the cosmic radio dipole at high fidelity. We outline several potentially transformational tests of cosmology to be carried out by means of SKA all-sky surveys.

preprint2015arXiv

Using dark energy to suppress power at small scales

The latest Planck results reconfirm the existence of a slight but chronic tension between the best-fit Cosmic Microwave Background (CMB) and low-redshift observables: power seems to be consistently lacking in the late universe across a range of observables (e.g.~weak lensing, cluster counts). We propose a two-parameter model for dark energy where the dark energy is sufficiently like dark matter at large scales to keep the CMB unchanged but where it does not cluster at small scales, preventing concordance collapse and erasing power. We thus exploit the generic scale-dependence of dark energy instead of the more usual time-dependence to address the tension in the data. The combination of CMB, distance and weak lensing data somewhat prefer our model to $Λ$CDM, at $Δχ^2=2.4$. Moreover, this improved solution has $σ_8=0.79 \pm 0.02$, consistent with the value implied by cluster counts.

preprint2014arXiv

Anisotropic stress as signature of non-standard propagation of gravitational waves

We make precise the heretofore ambiguous statement that anisotropic stress is a sign of a modification of gravity. We show that in cosmological solutions of very general classes of models extending gravity --- all scalar-tensor theories (Horndeski), Einstein-Aether models and bimetric massive gravity --- a direct correspondence exists between perfect fluids apparently carrying anisotropic stress and a modification in the propagation of gravitational waves. Since the anisotropic stress can be measured in a model-independent manner, a comparison of the behavior of gravitational waves from cosmological sources with large-scale-structure formation could in principle lead to new constraints on the theory of gravity.

preprint2014arXiv

Can Self-Ordering Scalar Fields explain the BICEP2 B-mode signal?

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

preprint2014arXiv

Can topological defects mimic the BICEP2 B-mode signal?

We show that the B-mode polarization signal detected at low multipoles by BICEP2 cannot be entirely due to topological defects. This would be incompatible with the high-multipole B-mode polarization data and also with existing temperature anisotropy data. Adding cosmic strings to a model with tensors, we find that B-modes on their own provide a comparable limit on the defects to that already coming from Planck satellite temperature data. We note that strings at this limit give a modest improvement to the best-fit of the B-mode data, at a somewhat lower tensor-to-scalar ratio of $r \simeq 0.15$.

preprint2014arXiv

Constraining topological defects with temperature and polarization anisotropies

We analyse the possible contribution of topological defects to cosmic microwave anisotropies, both temperature and polarisation. We allow for the presence of both inflationary scalars and tensors, and of polarised dust foregrounds that may contribute to or dominate the B-mode polarisation signal. We confirm and quantify our previous statements that topological defects on their own are a poor fit to the B-mode signal. However, adding topological defects to a models with a tensor component or a dust component improves the fit around $\ell=200$. Fitting simultaneously to both temperature and polarisation data, we find that textures fit almost as well as tensors ($Δχ^2 = 2.0$), while Abelian Higgs strings are ruled out as the sole source of the B-mode signal at low $\ell$. The 95% confidence upper limits on models combining defects and dust are $Gμ< 2.7\times 10^{-7}$ (Abelian Higgs strings), $Gμ< 9.8\times 10^{-7}$ (semilocal strings) and $Gμ< 7.3\times 10^{-7}$ (textures), a small reduction on the $Planck$ bounds. The most economical fit overall is obtained by the standard $Λ$CDM model with a polarised dust component.

preprint2014arXiv

Cosmic Microwave Background temperature and polarization anisotropies from the large-N limit of global defects

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

preprint2014arXiv

Cosmological perturbations and structure formation in nonlocal infrared modifications of general relativity

We study the cosmological consequences of a recently proposed nonlocal modification of general relativity, obtained by adding a term $m^2R\,\Box^{-2}R$ to the Einstein-Hilbert action. The model has the same number of parameters as $Λ$CDM, with $m$ replacing $Ω_Λ$, and is very predictive. At the background level, after fixing $m$ so as to reproduce the observed value of $Ω_M$, we get a pure prediction for the equation of state of dark energy as a function of redshift, $w_{\rm DE}(z)$, with $w_{\rm DE}(0)$ in the range $[-1.165,-1.135]$ as $Ω_M$ varies over the broad range $Ω_M\in [0.20,0.36]$. We find that the cosmological perturbations are well-behaved, and the model fully fixes the dark energy perturbations as a function of redshift $z$ and wavenumber $k$. The nonlocal model provides a good fit to supernova data and predicts deviations from General Relativity in structure formation and in weak lensing at the level of 3-4%, therefore consistent with existing data but readily detectable by future surveys. For the logarithmic growth factor we obtain $γ\simeq 0.53$, to be compared with $γ\simeq 0.55$ in $Λ$CDM. For the Newtonian potential on subhorizon scales our results are well fitted by $Ψ(a;k)=[1+μ_s a^s]Ψ_{\rm GR}(a;k)$ with a scale-independent $μ_s\simeq 0.09$ and $s\simeq 2$, while the anisotropic stress is negligibly small.

preprint2014arXiv

Distance-redshift relation in plane symmetric universes

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

preprint2014arXiv

Extending BEAMS to incorporate correlated systematic uncertainties

New supernova surveys such as the Dark Energy Survey, Pan-STARRS and the LSST will produce an unprecedented number of photometric supernova candidates, most with no spectroscopic data. Avoiding biases in cosmological parameters due to the resulting inevitable contamination from non-Ia supernovae can be achieved with the BEAMS formalism, allowing for fully photometric supernova cosmology studies. Here we extend BEAMS to deal with the case in which the supernovae are correlated by systematic uncertainties. The analytical form of the full BEAMS posterior requires evaluating 2^N terms, where N is the number of supernova candidates. This `exponential catastrophe' is computationally unfeasible even for N of order 100. We circumvent the exponential catastrophe by marginalising numerically instead of analytically over the possible supernova types: we augment the cosmological parameters with nuisance parameters describing the covariance matrix and the types of all the supernovae, τ_i, that we include in our MCMC analysis. We show that this method deals well even with large, unknown systematic uncertainties without a major increase in computational time, whereas ignoring the correlations can lead to significant biases and incorrect credible contours. We then compare the numerical marginalisation technique with a perturbative expansion of the posterior based on the insight that future surveys will have exquisite light curves and hence the probability that a given candidate is a Type Ia will be close to unity or zero, for most objects. Although this perturbative approach changes computation of the posterior from a 2^N problem into an N^2 or N^3 one, we show that it leads to biases in general through a small number of misclassifications, implying that numerical marginalisation is superior.

preprint2014arXiv

Fitting BICEP2 with defects, primordial gravitational waves and dust

In this work we discuss the possibility of cosmic defects being responsible for the B-mode signal measured by the BICEP2 collaboration. We also allow for the presence of other cosmological sources of B-modes such as inflationary gravitational waves and polarized dust foregrounds, which might contribute to or dominate the signal. On the one hand, we find that defects alone give a poor fit to the data points. On the other, we find that defects help to improve the fit at higher multipoles when they are considered alongside inflationary gravitational waves or polarized dust. Finally, we derive new defect constraints from models combining defects and dust. This proceeding is based on previous works [1,2].

preprint2014arXiv

General Relativistic N-body simulations in the weak field limit

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

preprint2014arXiv

Intrinsic uncertainty on the nature of dark energy

We argue that there is an intrinsic noise on measurements of the equation of state parameter $w=p/ρ$ from large-scale structure around us. The presence of the large-scale structure leads to an ambiguity in the definition of the background universe and thus there is a maximal precision with which we can determine the equation of state of dark energy. To study the uncertainty due to local structure, we model density perturbations stemming from a standard inflationary power spectrum by means of the exact Lemaître-Tolman-Bondi solution of Einstein's equation, and show that the usual distribution of matter inhomogeneities in a $Λ$CDM cosmology causes a variation of $w$ -- as inferred from distance measures -- of several percent. As we observe only one universe, or equivalently because of the cosmic variance, this uncertainty is systematic in nature.

preprint2014arXiv

N-body methods for relativistic cosmology

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

preprint2014arXiv

Non-local gravity and comparison with observational datasets

We study the cosmological predictions of two recently proposed non-local modifications of General Relativity. Both models have the same number of parameters as $Λ$CDM, with a mass parameter $m$ replacing the cosmological constant. We implement the cosmological perturbations of the non-local models into a modification of the CLASS Boltzmann code, and we make a full comparison to CMB, BAO and supernova data. We find that the non-local models fit these datasets as well as $Λ$CDM. For both non-local models parameter estimation using Planck+JLA+BAO data gives a value of $H_0$ higher than in $Λ$CDM, and in better agreement with the values obtained from local measurements.

preprint2014arXiv

PRISM (Polarized Radiation Imaging and Spectroscopy Mission): An Extended White Paper

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to ESA in May 2013 as a large-class mission for investigating within the framework of the ESA Cosmic Vision program a set of important scientific questions that require high resolution, high sensitivity, full-sky observations of the sky emission at wavelengths ranging from millimeter-wave to the far-infrared. PRISM's main objective is to explore the distant universe, probing cosmic history from very early times until now as well as the structures, distribution of matter, and velocity flows throughout our Hubble volume. PRISM will survey the full sky in a large number of frequency bands in both intensity and polarization and will measure the absolute spectrum of sky emission more than three orders of magnitude better than COBE FIRAS. The aim of this Extended White Paper is to provide a more detailed overview of the highlights of the new science that will be made possible by PRISM

preprint2014arXiv

The traces of anisotropic dark energy in light of Planck

We study a dark energy model with non-zero anisotropic stress, either linked to the dark energy density or to the dark matter density. We compute approximate solutions that allow to characterise the behaviour of the dark energy model and to assess the stability of the perturbations. We also determine the current limits on such an anisotropic stress from the cosmic microwave background data by the Planck satellite, and derive the corresponding constraints on the modified growth parameters like the growth index, the effective Newton's constant and the gravitational slip.

preprint2014arXiv

Towards the Future of Supernova Cosmology

For future surveys, spectroscopic follow-up for all supernovae will be extremely difficult. However, one can use light curve fitters, to obtain the probability that an object is a Type Ia. One may consider applying a probability cut to the data, but we show that the resulting non-Ia contamination can lead to biases in the estimation of cosmological parameters. A different method, which allows the use of the full dataset and results in unbiased cosmological parameter estimation, is Bayesian Estimation Applied to Multiple Species (BEAMS). BEAMS is a Bayesian approach to the problem which includes the uncertainty in the types in the evaluation of the posterior. Here we outline the theory of BEAMS and demonstrate its effectiveness using both simulated datasets and SDSS-II data. We also show that it is possible to use BEAMS if the data are correlated, by introducing a numerical marginalisation over the types of the objects. This is largely a pedagogical introduction to BEAMS with references to the main BEAMS papers.

preprint2013arXiv

Cosmology with Photometrically-Classified Type Ia Supernovae from the SDSS-II Supernova Survey

We present the cosmological analysis of 752 photometrically-classified Type Ia Supernovae (SNe Ia) obtained from the full Sloan Digital Sky Survey II (SDSS-II) Supernova (SN) Survey, supplemented with host-galaxy spectroscopy from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). Our photometric-classification method is based on the SN typing technique of Sako et al. (2011), aided by host galaxy redshifts (0.05<z<0.55). SNANA simulations of our methodology estimate that we have a SN Ia typing efficiency of 70.8%, with only 3.9% contamination from core-collapse (non-Ia) SNe. We demonstrate that this level of contamination has no effect on our cosmological constraints. We quantify and correct for our selection effects (e.g., Malmquist bias) using simulations. When fitting to a flat LambdaCDM cosmological model, we find that our photometric sample alone gives omega_m=0.24+0.07-0.05 (statistical errors only). If we relax the constraint on flatness, then our sample provides competitive joint statistical constraints on omega_m and omega_lambda, comparable to those derived from the spectroscopically-confirmed three-year Supernova Legacy Survey (SNLS3). Using only our data, the statistics-only result favors an accelerating universe at 99.96% confidence. Assuming a constant wCDM cosmological model, and combining with H0, CMB and LRG data, we obtain w=-0.96+0.10-0.10, omega_m=0.29+0.02-0.02 and omega_k=0.00+0.03-0.02 (statistical errors only), which is competitive with similar spectroscopically confirmed SNe Ia analyses. Overall this comparison is re-assuring, considering the lower redshift leverage of the SDSS-II SN sample (z<0.55) and the lack of spectroscopic confirmation used herein. These results demonstrate the potential of photometrically-classified SNe Ia samples in improving cosmological constraints.

preprint2013arXiv

Crystal nucleation and near-epitaxial growth in nacre

Nacre is a layered, iridescent lining found inside many mollusk shells, with a unique brick-and-mortar periodic structure at the sub-micron scale, and remarkable resistance to fracture. Despite extensive studies, it remains unclear how nacre forms. Here we present 20-nm, 2°-resolution Polarization-dependent Imaging Contrast (PIC) images of shells from 15 mollusk shell species, mapping nacre tablets and their orientation patterns, showing where new crystal orientations appear and how they propagate across organic sheets as nacre grows. In all shells we found stacks of co-oriented aragonite (CaCO3) tablets arranged into vertical columns or staggered diagonally. Only near the nacre-prismatic boundary are disordered crystals nucleated, as spherulitic aragonite. Overgrowing nacre tablet crystals are most frequently co-oriented with the underlying spherulitic aragonite or with another tablet, connected by mineral bridges. Therefore aragonite crystal growth in nacre is epitaxial or near-epitaxial, with abrupt or gradual changes in orientation, with c-axes within 20°. Based on these data, we propose that there is one mineral bridge per tablet, and that "bridge-tilting" is a possible mechanism to introduce small, gradual or abrupt changes in the orientation of crystals within a stack of tablets as nacre grows.

preprint2013arXiv

Host Galaxy Spectra and Consequences for SN Typing From The SDSS SN Survey

We present the spectroscopy from 5254 galaxies that hosted supernovae (SNe) or other transient events in the Sloan Digital Sky Survey II (SDSS-II). Obtained during SDSS-I, SDSS-II, and the Baryon Oscillation Spectroscopic Survey (BOSS), this sample represents the largest systematic, unbiased, magnitude limited spectroscopic survey of supernova (SN) host galaxies. Using the host galaxy redshifts, we test the impact of photometric SN classification based on SDSS imaging data with and without using spectroscopic redshifts of the host galaxies. Following our suggested scheme, there are a total of 1166 photometrically classified SNe Ia when using a flat redshift prior and 1126 SNe Ia when the host spectroscopic redshift is assumed. For 1024 (87.8%) candidates classified as likely SNe Ia without redshift information, we find that the classification is unchanged when adding the host galaxy redshift. Using photometry from SDSS imaging data and the host galaxy spectra, we also report host galaxy properties for use in future nalysis of SN astrophysics. Finally, we investigate the differences in the interpretation of the light curve properties with and without knowledge of the redshift. When using the SALT2 light curve fitter, we find a 21% increase in the number of fits that converge when using the spectroscopic redshift. Without host galaxy redshifts, we find that SALT2 light curve fits are systematically biased towards lower photometric redshift estimates and redder colors in the limit of low signal-to-noise data. The general improvements in performance of the light curve fitter and the increased diversity of the host galaxy sample highlights the importance of host galaxy spectroscopy for current photometric SN surveys such as the Dark Energy Survey and future surveys such as the Large Synoptic Survey Telescope.

preprint2013arXiv

Model-independent constraints on the cosmological anisotropic stress

The effective anisotropic stress or gravitational slip $η=-Φ/Ψ$ is a key variable in the characterisation of the physical origin of the dark energy, as it allows to test for a non-minimal coupling of the dark sector to gravity in the Jordan frame. It is however important to use a fully model-independent approach when measuring $η$ to avoid introducing a theoretical bias into the results. In this paper we forecast the precision with which future large surveys can determine $η$ in a way that only relies on directly observable quantities. In particular, we do not assume anything concerning the initial spectrum of perturbations, nor on its evolution outside the observed redshift range, nor on the galaxy bias. We first leave $η$ free to vary in space and time and then we model it as suggested in Horndeski models of dark energy. Among our results, we find that a future large scale lensing and clustering survey can constrain $η$ to within 10% if $k$-independent, and to within 60% or better at $k=0.1 h/$Mpc if it is restricted to follow the Horndeski model.

preprint2013arXiv

Probing Dark Energy through Scale Dependence

We consider the consequences of having no prior knowledge of the true dark energy model for the interpretation of cosmological observations. The magnitude of redshift-space distortions and weak-lensing shear is determined by the metric on the geodesics of which galaxies and light propagate. We show that, given precise enough observations, we can use these data to completely reconstruct the metric on our past lightcone and therefore to measure the scale- and time-dependence of the anisotropic stress and the evolution of the gravitational potentials in a model-independent manner. Since both dark matter and dark energy affect the visible sector only through the gravitational field they produce, they are inseparable without a model for dark energy: galaxy bias cannot be measured and therefore the distribution of dark matter determined; the peculiar velocity of dark matter can be identified with that of the galaxies only when the equivalence principle holds. Given these limitations, we show how one can nonetheless build tests for classes of dark energy models which depend on making measurements at multiple scales at a particular redshift. They are null tests on the model-independent observables, do not require modeling evolution in time and do not require any parametrization of the free functions of these models, such as the sound speed. We show how one can rule out or constrain the whole class of the most-general scalar-tensor theories even without assuming the quasi-static limit.

preprint2012arXiv

BEAMS: separating the wheat from the chaff in supernova analysis

We introduce Bayesian Estimation Applied to Multiple Species (BEAMS), an algorithm designed to deal with parameter estimation when using contaminated data. We present the algorithm and demonstrate how it works with the help of a Gaussian simulation. We then apply it to supernova data from the Sloan Digital Sky Survey (SDSS), showing how the resulting confidence contours of the cosmological parameters shrink significantly.

preprint2012arXiv

Consistent perturbations in an imperfect fluid

We present a new prescription for analysing cosmological perturbations in a more-general class of scalar-field dark-energy models where the energy-momentum tensor has an imperfect-fluid form. This class includes Brans-Dicke models, f(R) gravity, theories with kinetic gravity braiding and generalised galileons. We employ the intuitive language of fluids, allowing us to explicitly maintain a dependence on physical and potentially measurable properties. We demonstrate that hydrodynamics is not always a valid description for describing cosmological perturbations in general scalar-field theories and present a consistent alternative that nonetheless utilises the fluid language. We apply this approach explicitly to a worked example: k-essence non-minimally coupled to gravity. This is the simplest case which captures the essential new features of these imperfect-fluid models. We demonstrate the generic existence of a new scale separating regimes where the fluid is perfect and imperfect. We obtain the equations for the evolution of dark-energy density perturbations in both these regimes. The model also features two other known scales: the Compton scale related to the breaking of shift symmetry and the Jeans scale which we show is determined by the speed of propagation of small scalar-field perturbations, i.e. causality, as opposed to the frequently used definition of the ratio of the pressure and energy-density perturbations.

preprint2012arXiv

Nonlinear cosmological consistency relations and effective matter stresses

We propose a fully nonlinear framework to construct consistency relations for testing generic cosmological scenarios using the evolution of large scale structure. It is based on the covariant approach in combination with a frame that is purely given by the metric, the normal frame. As an example, we apply this framework to the LambdaCDM model, by extending the usual first order conditions on the metric potentials to second order, where the two potentials start to differ from each other. We argue that working in the normal frame is not only a practical choice but also helps with the physical interpretation of nonlinear dynamics. In this frame, effective pressures and anisotropic stresses appear at second order in perturbation theory, even for "pressureless" dust. We quantify their effect and compare them, for illustration, to the pressure of a generic clustering dark energy uid and the anisotropic stress in the DGP model. Besides, we also discuss the effect of a mismatch of the potentials on the determination of galaxy bias.

preprint2012arXiv

Observables and unobservables in dark energy cosmologies

The aim of this paper is to answer the following two questions: (1) Given cosmological observations of the expansion history and linear perturbations in a range of redshifts and scales as precise as is required, which of the properties of dark energy could actually be reconstructed without imposing any parameterization? (2) Are these observables sufficient to rule out not just a particular dark energy model, but the entire general class of viable models comprising a single scalar field? This paper bears both good and bad news. On one hand, we find that the goal of reconstructing dark energy models is fundamentally limited by the unobservability of the present values of the matter density Omega_m0, the perturbation normalization sigma_8 as well as the present matter power spectrum. On the other, we find that, under certain conditions, cosmological observations can nonetheless rule out the entire class of the most general single scalar-field models, i.e. those based on the Horndeski Lagrangian.

preprint2012arXiv

Particle motion in weak relativistic gravitational fields

We derive the geodesic equation of motion in the presence of weak gravitational fields produced by relativistic sources such as cosmic strings, decomposed into scalar, vector and tensor parts. We find that the vector (gravito-magnetic) force is an important contributor, and for non-relativistic particles we recover the well-known result for the impulse from a moving straight string. Our results can be straightforwardly incorporated into N-body simulations to allow for the presence of cosmic defects or other sources of weak gravitational fields.

preprint2012arXiv

The phenomenological approach to modeling the dark energy

In this mini-review we discuss first why we should investigate cosmological models beyond LCDM. We then show how to describe dark energy or modified gravity models in a fluid language with the help of one background and two perturbation quantities. We review a range of dark energy models and study how they fit into the phenomenological framework, including generalizations like phantom crossing, sound speeds different from c and non-zero anisotropic stress, and how these effective quantities are linked to the underlying physical models. We also discuss the limits of what can be measured with cosmological data, and some challenges for the framework.

preprint2011arXiv

Anisotropic stress and stability in modified gravity models

The existence of anisotropic stress of a purely geometrical origin seems to be a characteristic of higher order gravity models, and has been suggested as a probe to test these models observationally, for example in weak lensing experiments. In this paper, we seek to find a class of higher order gravity models of f(R,G) type that would give us a zero anisotropic stress and study the consequences for the viability of the actual model. For the special case of a de Sitter background, we identify a subclass of models with the desired property. We also find a direct link between anisotropic stress and the stability of the model as well as the presence of extra degrees of freedom, which seems to be a general feature of higher order gravity models. Particularly, setting the anisotropic stress equal to zero for a de Sitter background leads to a singularity that makes it impossible to reach the de Sitter evolution.

preprint2011arXiv

Observation of new intrinsic properties of VO2

The single crystal VO2, exihibiting a first-order metal-insulator transition (MIT) at 67.2 degrees C and an insulator-insulator transition (IIT) at ~49.7 degrees C, is grown. From synchrotron-based x-ray microdiffraction analysis, the IIT shows structural phase transition (SPT) of monoclinic M2 to M1 phases while the MIT displays M1 to rutile R phases. The IIT exhibits percolative SPT while the MIT shows abrupt transition width of < 0.02 degrees C, supporting Mott's prediction. The MIT occurs non-percolatively with a sharp boundary between R and M1 phases. The MIT onset temperature shows significant variation.

preprint2011arXiv

Parameter Estimation with BEAMS in the presence of biases and correlations

The original formulation of BEAMS - Bayesian Estimation Applied to Multiple Species - showed how to use a dataset contaminated by points of multiple underlying types to perform unbiased parameter estimation. An example is cosmological parameter estimation from a photometric supernova sample contaminated by unknown Type Ibc and II supernovae. Where other methods require data cuts to increase purity, BEAMS uses all of the data points in conjunction with their probabilities of being each type. Here we extend the BEAMS formalism to allow for correlations between the data and the type probabilities of the objects as can occur in realistic cases. We show with simple simulations that this extension can be crucial, providing a 50% reduction in parameter estimation variance when such correlations do exist. We then go on to perform tests to quantify the importance of the type probabilities, one of which illustrates the effect of biasing the probabilities in various ways. Finally, a general presentation of the selection bias problem is given, and discussed in the context of future photometric supernova surveys and BEAMS, which lead to specific recommendations for future supernova surveys.

preprint2011arXiv

Photometric Supernova Cosmology with BEAMS and SDSS-II

Supernova cosmology without spectroscopic confirmation is an exciting new frontier which we address here with the Bayesian Estimation Applied to Multiple Species (BEAMS) algorithm and the full three years of data from the Sloan Digital Sky Survey II Supernova Survey (SDSS-II SN). BEAMS is a Bayesian framework for using data from multiple species in statistical inference when one has the probability that each data point belongs to a given species, corresponding in this context to different types of supernovae with their probabilities derived from their multi-band lightcurves. We run the BEAMS algorithm on both Gaussian and more realistic SNANA simulations with of order 10^4 supernovae, testing the algorithm against various pitfalls one might expect in the new and somewhat uncharted territory of photometric supernova cosmology. We compare the performance of BEAMS to that of both mock spectroscopic surveys and photometric samples which have been cut using typical selection criteria. The latter typically are either biased due to contamination or have significantly larger contours in the cosmological parameters due to small data-sets. We then apply BEAMS to the 792 SDSS-II photometric supernovae with host spectroscopic redshifts. In this case, BEAMS reduces the area of the (Ω_m,Ω_Λ) contours by a factor of three relative to the case where only spectroscopically confirmed data are used (297 supernovae). In the case of flatness, the constraints obtained on the matter density applying BEAMS to the photometric SDSS-II data are Ω_m(BEAMS)=0.194\pm0.07. This illustrates the potential power of BEAMS for future large photometric supernova surveys such as LSST.

preprint2010arXiv

A dark energy view of inflation

Traditionally, inflationary models are analyzed in terms of parameters such as the scalar spectral index ns and the tensor to scalar ratio r, while dark energy models are studied in terms of the equation of state parameter w. Motivated by the fact that both deal with periods of accelerated expansion, we study the evolution of w during inflation, in order to derive observational constraints on its value during an earlier epoch likely dominated by a dynamic form of dark energy. We find that the cosmic microwave background and large-scale structure data is consistent with w_inflation=-1 and provides an upper limit of 1+w <~ 0.02. Nonetheless, an exact de Sitter expansion with a constant w=-1 is disfavored since this would result in ns=1.

preprint2010arXiv

CMB power spectra from cosmic strings: predictions for the Planck satellite and beyond

We present a significant improvement over our previous calculations of the cosmic string contribution to cosmic microwave background (CMB) power spectra, with particular focus on sub-WMAP angular scales. These smaller scales are relevant for the now-operational Planck satellite and additional sub-orbital CMB projects that have even finer resolutions. We employ larger Abelian Higgs string simulations than before and we additionally model and extrapolate the statistical measures from our simulations to smaller length scales. We then use an efficient means of including the extrapolations into our Einstein-Boltzmann calculations in order to yield accurate results over the multipole range 2 < l < 4000. Our results suggest that power-law behaviour cuts in for l > 3000 in the case of the temperature power spectrum, which then allows cautious extrapolation to even smaller scales. We find that a string contribution to the temperature power spectrum making up 10% of power at l=10 would be larger than the Silk-damped primary adiabatic contribution for l > 3500. Astrophysical contributions such as the Sunyaev-Zeldovich effect also become important at these scales and will reduce the sensitivity to strings, but these are potentially distinguishable by their frequency-dependence.

preprint2010arXiv

Detecting and distinguishing topological defects in future data from the CMBPol satellite

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σ level. We also confirm that a tensor contribution of r = 0.0018 can be detected at over 3σ, in agreement with the CMBpol mission concept study. These results are supported by a model selection analysis.

preprint2010arXiv

Fingerprinting Dark Energy II: weak lensing and galaxy clustering tests

The characterization of dark energy is a central task of cosmology. To go beyond a cosmological constant, we need to introduce at least an equation of state and a sound speed and consider observational tests that involve perturbations. If dark energy is not completely homogeneous on observable scales then the Poisson equation is modified and dark matter clustering is directly affected. One can then search for observational effects of dark energy clustering using dark matter as a probe. In this paper we exploit an analytical approximate solution of the perturbation equations in a general dark energy cosmology to analyze the performance of next-decade large scale surveys in constraining equation of state and sound speed. We find that tomographic weak lensing and galaxy redshift surveys can constrain the sound speed of the dark energy only if the latter is small, of the order of $c_{s}\lesssim0.01$ (in units of $c$). For larger sound speeds the error grows to 100% and more. We conclude that large scale structure observations contain very little information about the perturbations in canonical scalar field models with a sound speed of unity. Nevertheless, they are able to detect the presence of "cold" dark energy, i.e. a dark energy with non-relativistic speed of sound.

preprint2010arXiv

Results from the Supernova Photometric Classification Challenge

We report results from the Supernova Photometric Classification Challenge (SNPCC), a publicly released mix of simulated supernovae (SNe), with types (Ia, Ibc, and II) selected in proportion to their expected rate. The simulation was realized in the griz filters of the Dark Energy Survey (DES) with realistic observing conditions (sky noise, point-spread function and atmospheric transparency) based on years of recorded conditions at the DES site. Simulations of non-Ia type SNe are based on spectroscopically confirmed light curves that include unpublished non-Ia samples donated from the Carnegie Supernova Project (CSP), the Supernova Legacy Survey (SNLS), and the Sloan Digital Sky Survey-II (SDSS-II). A spectroscopically confirmed subset was provided for training. We challenged scientists to run their classification algorithms and report a type and photo-z for each SN. Participants from 10 groups contributed 13 entries for the sample that included a host-galaxy photo-z for each SN, and 9 entries for the sample that had no redshift information. Several different classification strategies resulted in similar performance, and for all entries the performance was significantly better for the training subset than for the unconfirmed sample. For the spectroscopically unconfirmed subset, the entry with the highest average figure of merit for classifying SNe~Ia has an efficiency of 0.96 and an SN~Ia purity of 0.79. As a public resource for the future development of photometric SN classification and photo-z estimators, we have released updated simulations with improvements based on our experience from the SNPCC, added samples corresponding to the Large Synoptic Survey Telescope (LSST) and the SDSS, and provided the answer keys so that developers can evaluate their own analysis.

preprint2010arXiv

Statistical Classification Techniques for Photometric Supernova Typing

Future photometric supernova surveys will produce vastly more candidates than can be followed up spectroscopically, highlighting the need for effective classification methods based on lightcurves alone. Here we introduce boosting and kernel density estimation techniques which have minimal astrophysical input, and compare their performance on 20,000 simulated Dark Energy Survey lightcurves. We demonstrate that these methods are comparable to the best template fitting methods currently used, and in particular do not require the redshift of the host galaxy or candidate. However both methods require a training sample that is representative of the full population, so typical spectroscopic supernova subsamples will lead to poor performance. To enable the full potential of such blind methods, we recommend that representative training samples should be used and so specific attention should be given to their creation in the design phase of future photometric surveys.

preprint2010arXiv

The local B-polarization of the CMB: a very sensitive probe of cosmic defects

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

preprint2009arXiv

Optimizing baryon acoustic oscillation surveys II: curvature, redshifts, and external datasets

We extend our study of the optimization of large baryon acoustic oscillation (BAO) surveys to return the best constraints on the dark energy, building on Paper I of this series (Parkinson et al. 2007). The survey galaxies are assumed to be pre-selected active, star-forming galaxies observed by their line emission with a constant number density across the redshift bin. Star-forming galaxies have a redshift desert in the region 1.6 < z < 2, and so this redshift range was excluded from the analysis. We use the Seo & Eisenstein (2007) fitting formula for the accuracies of the BAO measurements, using only the information for the oscillatory part of the power spectrum as distance and expansion rate rulers. We go beyond our earlier analysis by examining the effect of including curvature on the optimal survey configuration and updating the expected `prior' constraints from Planck and SDSS. We once again find that the optimal survey strategy involves minimizing the exposure time and maximizing the survey area (within the instrumental constraints), and that all time should be spent observing in the low-redshift range (z<1.6) rather than beyond the redshift desert, z>2. We find that when assuming a flat universe the optimal survey makes measurements in the redshift range 0.1 < z <0.7, but that including curvature as a nuisance parameter requires us to push the maximum redshift to 1.35, to remove the degeneracy between curvature and evolving dark energy. The inclusion of expected other data sets (such as WiggleZ, BOSS and a stage III SN-Ia survey) removes the necessity of measurements below redshift 0.9, and pushes the maximum redshift up to 1.5. We discuss considerations in determining the best survey strategy in light of uncertainty in the true underlying cosmological model.

preprint2009arXiv

Unified dark energy and dark matter from a scalar field different from quintessence

We explore unification of dark matter and dark energy in a theory containing a scalar field of non-Lagrangian type, obtained by direct insertion of a kinetic term into the energy-momentum tensor. This scalar is different from quintessence, having an equation of state between -1 and 0 and a zero sound speed in its rest frame. We solve the equations of motion for an exponential potential via a rewriting as an autonomous system, and demonstrate the observational viability of the scenario, for sufficiently small exponential potential parameter λ, by comparison to a compilation of kinematical cosmological data.

preprint2008arXiv

Testing backreaction effects with observations

In order to quantitatively test the ability of averaged inhomogeneous cosmologies to correctly describe observations of the large scale properties of the Universe, we introduce a smoothed template metric corresponding to a constant spatial curvature model at any time, but with an evolving curvature parameter. This metric is used to compute quantities along an approximate effective lightcone of the averaged model of the Universe. As opposed to the standard Friedmann model, we parameterize this template metric by exact scaling properties of an averaged inhomogeneous cosmology, and we also motivate this form of the metric by results on a geometrical smoothing of inhomogeneous cosmological hypersurfaces. We test our hypothesis for the template metric against supernova data and the position of the CMB peaks, and infer the goodness--of--fit and parameter uncertainties. We find that averaged inhomogeneous models can reproduce the observations without requiring an additional Dark Energy component (though a volume acceleration is still needed), and that current data do not disfavour our main assumption on the effective lightcone structure. We also show that the experimental uncertainties on the angular diameter distance and the Hubble parameter from Baryon Acoustic Oscillations measurements -- forseen in future surveys like the proposed EUCLID satellite project -- are sufficiently small to distinguish between a FLRW template geometry and the template geometry with consistently evolving curvature. (Abridged)

preprint2007arXiv

The dark degeneracy: On the number and nature of dark components

We use that gravity probes only the total energy momentum tensor to show how this leads to a degeneracy for generalised dark energy models. Because of this degeneracy, Omega_m cannot be measured. We demonstrate this explicitely by showing that the CMB and supernova data is compatible with very large and very small values of Omega_m for a specific family of dark energy models. We also show that for the same reason interacting dark energy is always equivalent to a family of non-interacting models. We argue that it is better to face this degeneracy and to parametrise the actual observables.

preprint2004arXiv

Model-independent dark energy test with sigma_8 using results from the Wilkinson Microwave Anisotropy Probe

By combining the recent WMAP measurements of the cosmic microwave background anisotropies and the results of the recent luminosity distance measurements to type-Ia supernovae, we find that the normalization of the matter power spectrum on cluster scales, sigma_8, can be used to discriminate between dynamical models of dark energy (quintessence models) and a conventional cosmological constant model (LCDM).

preprint2002arXiv

The role of point-like topological excitations at criticality: from vortices to global monopoles

We determine the detailed thermodynamic behavior of vortices in the O(2) scalar model in 2D and of global monopoles in the O(3) model in 3D. We construct new numerical techniques, based on cluster decomposition algorithms, to analyze the point defect configurations. We find that these criteria produce results for the Kosterlitz-Thouless temperature in agreement with a topological transition between a polarizable insulator and a conductor, at which free topological charges appear in the system. For global monopoles we find no pair unbinding transition. Instead a transition to a dense state where pairs are no longer distinguishable occurs at T<Tc, without leading to long range disorder. We produce both extensive numerical evidence of this behavior as well as a semi-analytic treatment of the partition function for defects. General expectations for N=D>3 are drawn, based on the observed behavior.

preprint1996arXiv

Microwave Background Anisotropies Induced by Global Scalar Fields: The Large N Limit

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