Showing papers in "Journal of Cosmology and Astroparticle Physics in 2016"

Science with the space-based interferometer eLISA. II: Gravitational waves from cosmological phase transitions

Abstract: We investigate the potential for the eLISA space-based interferometer to detect the stochastic gravitational wave background produced by strong first-order cosmological phase transitions. We discuss the resulting contributions from bubble collisions, magnetohydrodynamic turbulence, and sound waves to the stochastic background, and estimate the total corresponding signal predicted in gravitational waves. The projected sensitivity of eLISA to cosmological phase transitions is computed in a model-independent way for various detector designs and configurations. By applying these results to several specific models, we demonstrate that eLISA is able to probe many well-motivated scenarios beyond the Standard Model of particle physics predicting strong first-order cosmological phase transitions in the early Universe.

A 6% measurement of the Hubble parameter at z∼0.45: direct evidence of the epoch of cosmic re-acceleration

Abstract: Deriving the expansion history of the Universe is a major goal of modern cosmology. To date, the most accurate measurements have been obtained with Type Ia Supernovae (SNe) and Baryon Acoustic Oscillations (BAO), providing evidence for the existence of a transition epoch at which the expansion rate changes from decelerated to accelerated. However, these results have been obtained within the framework of specific cosmological models that must be implicitly or explicitly assumed in the measurement. It is therefore crucial to obtain measurements of the accelerated expansion of the Universe independently of assumptions on cosmological models. Here we exploit the unprecedented statistics provided by the Baryon Oscillation Spectroscopic Survey (BOSS, [1-3]) Data Release 9 to provide new constraints on the Hubble parameter H(z) using the cosmic chronometers approach. We extract a sample of more than 130000 of the most massive and passively evolving galaxies, obtaining five new cosmology-independent H(z) measurements in the redshift range 0.3 < z < 0.5, with an accuracy of ~11–16% incorporating both statistical and systematic errors. Once combined, these measurements yield a 6% accuracy constraint of H(z = 0.4293) = 91.8 ± 5.3 km/s/Mpc. The new data are crucial to provide the first cosmology-independent determination of the transition redshift at high statistical significance, measuring zt = 0.4 ± 0.1, and to significantly disfavor the null hypothesis of no transition between decelerated and accelerated expansion at 99.9% confidence level. This analysis highlights the wide potential of the cosmic chronometers approach: it permits to derive constraints on the expansion history of the Universe with results competitive with standard probes, and most importantly, being the estimates independent of the cosmological model, it can constrain cosmologies beyond—and including—the ΛCDM model.

The trouble with $H_0$

Abstract: We perform a comprehensive cosmological study of the H0 tension between the direct local measurement and the model-dependent value inferred from the Cosmic Microwave Background. With the recent measurement of H0 this tension has raised to more than 3 σ. We consider changes in the early time physics without modifying the late time cosmology. We also reconstruct the late time expansion history in a model independent way with minimal assumptions using distance measurements from Baryon Acoustic Oscillations and Type Ia Supernovae, finding that at z < 0.6 the recovered shape of the expansion history is less than 5% different than that of a standard ΛCDM model. These probes also provide a model insensitive constraint on the low-redshift standard ruler, measuring directly the combination rsh where H0 = h × 100 Mpc−1km/s and rs is the sound horizon at radiation drag (the standard ruler), traditionally constrained by CMB observations. Thus rs and H0 provide absolute scales for distance measurements (anchors) at opposite ends of the observable Universe. We calibrate the cosmic distance ladder and obtain a model-independent determination of the standard ruler for acoustic scale, rs. The tension in H0 reflects a mismatch between our determination of rs and its standard, CMB-inferred value. Without including high-l Planck CMB polarization data (i.e., only considering the ``recommended baseline" low-l polarisation and temperature and the high l temperature data), a modification of the early-time physics to include a component of dark radiation with an effective number of species around 0.4 would reconcile the CMB-inferred constraints, and the local H0 and standard ruler determinations. The inclusion of the ``preliminary" high-l Planck CMB polarisation data disfavours this solution.

Degenerate higher derivative theories beyond Horndeski: evading the Ostrogradski instability

Abstract: Theories with higher order time derivatives generically suffer from ghost-like instabilities, known as Ostrogradski instabilities. This fate can be avoided by considering ``degenerate'' Lagrangians, whose kinetic matrix cannot be inverted, thus leading to constraints between canonical variables and a reduced number of physical degrees of freedom. In this work, we derive in a systematic way the degeneracy conditions for scalar-tensor theories that depend quadratically on second order derivatives of a scalar field. We thus obtain a classification of all degenerate theories within this class of scalar-tensor theories. The quartic Horndeski Lagrangian and its extension beyond Horndeski belong to these degenerate cases. We also identify new families of scalar-tensor theories with the property that they are degenerate despite the nondegeneracy of the purely scalar part of their Lagrangian.

Physics reach of the XENON1T dark matter experiment

Abstract: The XENON1T experiment is currently in the commissioning phase at the Laboratori Nazionali del Gran Sasso, Italy. In this article we study the experiment's expected sensitivity to the spin-independent WIMP-nucleon interaction cross section, based on Monte Carlo predictions of the electronic and nuclear recoil backgrounds. The total electronic recoil background in 1 tonne fiducial volume and (1, 12) keV electronic recoil equivalent energy region, before applying any selection to discriminate between electronic and nuclear recoils, is (1.80 ± 0.15) · 10(−)(4) (kg·day·keV)(−)(1), mainly due to the decay of (222)Rn daughters inside the xenon target. The nuclear recoil background in the corresponding nuclear recoil equivalent energy region (4, 50) keV, is composed of (0.6 ± 0.1) (t·y)(−)(1) from radiogenic neutrons, (1.8 ± 0.3) · 10(−)(2) (t·y)(−)(1) from coherent scattering of neutrinos, and less than 0.01 (t·y)(−)(1) from muon-induced neutrons. The sensitivity of XENON1T is calculated with the Profile Likelihood Ratio method, after converting the deposited energy of electronic and nuclear recoils into the scintillation and ionization signals seen in the detector. We take into account the systematic uncertainties on the photon and electron emission model, and on the estimation of the backgrounds, treated as nuisance parameters. The main contribution comes from the relative scintillation efficiency Script L(eff), which affects both the signal from WIMPs and the nuclear recoil backgrounds. After a 2 y measurement in 1 t fiducial volume, the sensitivity reaches a minimum cross section of 1.6 · 10(−)(47) cm(2) at m(χ) = 50 GeV/c(2).

DARWIN: towards the ultimate dark matter detector

Abstract: DARk matter WImp search with liquid xenoN (DARWIN(2)) will be an experiment for the direct detection of dark matter using a multi-ton liquid xenon time projection chamber at its core. Its primary g ...

Science with the space-based interferometer LISA. IV: probing inflation with gravitational waves

Abstract: We investigate the potential for the LISA space-based interferometer to detect the stochastic gravitational wave background produced from different mechanisms during inflation. Focusing on well-motivated scenarios, we study the resulting contributions from particle production during inflation, inflationary spectator fields with varying speed of sound, effective field theories of inflation with specific patterns of symmetry breaking and models leading to the formation of primordial black holes. The projected sensitivities of LISA are used in a model-independent way for various detector designs and configurations. We demonstrate that LISA is able to probe these well-motivated inflationary scenarios beyond the irreducible vacuum tensor modes expected from any inflationary background.

Breaking a dark degeneracy with gravitational waves

Abstract: We identify a scalar-tensor model embedded in the Horndeski action whose cosmological background and linear scalar fluctuations are degenerate with the concordance cosmology. The model admits a self-accelerated background expansion at late times that is stable against perturbations with a sound speed attributed to the new field that is equal to the speed of light. While degenerate in scalar fluctuations, self-acceleration of the model implies a present cosmological tensor mode propagation at 95 % of the speed of light with a damping of the wave amplitude that is 5 % less efficient than in general relativity. We show that these discrepancies are endemic to self-accelerated Horndeski theories with degenerate large-scale structure and are tested with measurements of gravitational waves emitted by events at cosmological distances. Hence, gravitational-wave cosmology breaks the dark degeneracy in observations of the large-scale structure between two fundamentally different explanations of cosmic acceleration—a cosmological constant and a scalar-tensor modification of gravity. The gravitational wave event GW150914 recently detected with the aLIGO instruments and its potential association with a weak short gamma-ray burst observed with the Fermi GBM experiment may have provided this crucial measurement.

Extended scalar-tensor theories of gravity

Abstract: We study new consistent scalar-tensor theories of gravity recently introduced by Langlois and Noui with potentially interesting cosmological applications. We derive the conditions for the existence of a primary constraint that prevents the propagation of an additional dangerous mode associated with higher order equations of motion. We then classify the most general, consistent scalar-tensor theories that are at most quadratic in the second derivatives of the scalar field. In addition, we investigate the possible connection between these theories and (beyond) Horndeski through conformal and disformal transformations. Finally, we point out that these theories can be associated with new operators in the effective field theory of dark energy, which might open up new possibilities to test dark energy models in future surveys.

Limits to dark matter annihilation cross-section from a combined analysis of MAGIC and Fermi-LAT observations of dwarf satellite galaxies

Abstract: We present the first joint analysis of gamma-ray data from the MAGIC Cherenkov telescopes and the Fermi Large Area Telescope (LAT) to search for gamma-ray signals from dark matter annihilation in dwarf satellite galaxies. We combine 158 hours of Segue 1 observations with MAGIC with 6-year observations of 15 dwarf satellite galaxies by the Fermi-LAT. We obtain limits on the annihilation cross-section for dark matter particle masses between 10 GeV and 100 TeV – the widest mass range ever explored by a single gamma-ray analysis. These limits improve on previously published Fermi-LAT and MAGIC results by up to a factor of two at certain masses. Our new inclusive analysis approach is completely generic and can be used to perform a global, sensitivity-optimized dark matter search by combining data from present and future gamma-ray and neutrino detectors.

Relic neutrino decoupling with flavour oscillations revisited

Abstract: We study the decoupling process of neutrinos in the early universe in the presence of three-flavour oscillations. The evolution of the neutrino spectra is found by solving the corresponding momentum-dependent kinetic equations for the neutrino density matrix, including for the first time the proper collision integrals for both diagonal and off-diagonal elements. This improved calculation modifies the evolution of the off-diagonal elements of the neutrino density matrix and changes the deviation from equilibrium of the frozen neutrino spectra. However, it does not vary the contribution of neutrinos to the cosmological energy density in the form of radiation, usually expressed in terms of the effective number of neutrinos, Neff. We find a value of Neff = 3.045, in agreement with previous theoretical calculations and consistent with the latest analysis of Planck data. This result does not depend on the ordering of neutrino masses. We also consider the effect of non-standard neutrino-electron interactions (NSI), predicted in many theoretical models where neutrinos acquire mass. For two sets of NSI parameters allowed by present data, we find that Neff can be reduced down to 3.040 or enhanced up to 3.059.

Cosmic ray nuclei, antiprotons and gamma rays in the galaxy: a new diffusion model

Abstract: We model the transport of cosmic ray nuclei in the galaxy by means of a new numerical code. At variance with previous numerical models, ours accounts for a generic spatial distribution of the diffusion coefficient. We found that in the case of radially uniform diffusion, the main secondary/primary ratios (B/C, N/O and sub-Fe/Fe) and the modulated antiproton spectrum consistently match the available observations. Convection and re-acceleration do not seem to be required in the energy range that we consider: 1≤E≤103 GeV/nucleon. We generalize these results, accounting for radial dependence of the diffusion coefficient, which is assumed to trace that of the cosmic ray sources. While this does not affect the prediction of secondary/primary ratios, the simulated longitude profile of the diffuse γ-ray emission is significantly different from that for the uniform case and may agree with EGRET measurements without requiring invoking ad hoc assumptions on the galactic gas density distribution.

CRPropa 3-a public astrophysical simulation framework for propagating extraterrestrial ultra-high energy particles

Abstract: We present the simulation framework CRPropa version 3 designed for efficient development of astrophysical predictions for ultra-high energy particles. Users can assemble modules of the most relevant propagation effects in galactic and extragalactic space, include their own physics modules with new features, and receive on output primary and secondary cosmic messengers including nuclei, neutrinos and photons. In extension to the propagation physics contained in a previous CRPropa version, the new version facilitates high-performance computing and comprises new physical features such as an interface for galactic propagation using lensing techniques, an improved photonuclear interaction calculation, and propagation in time dependent environments to take into account cosmic evolution effects in anisotropy studies and variable sources. First applications using highlighted features are presented as well.

Gravitational waves at interferometer scales and primordial black holes in axion inflation

Abstract: We study the prospects of detection at terrestrial and space interferometers, as well as at pulsar timing array experiments, of a stochastic gravitational wave background which can be produced in models of axion inflation. This potential signal, and the development of these experiments, open a new window on inflation on scales much smaller than those currently probed with Cosmic Microwave Background and Large Scale Structure measurements. The sourced signal generated in axion inflation is an ideal candidate for such searches, since it naturally grows at small scales, and it has specific properties (chirality and non-gaussianity) that can distinguish it from an astrophysical background. We study under which conditions such a signal can be produced at an observable level, without the simultaneous overproduction of scalar perturbations in excess of what is allowed by the primordial black hole limits. We also explore the possibility that scalar perturbations generated in a modified version of this model may provide a distribution of primordial black holes compatible with the current bounds, that can act as a seeds of the present black holes in the universe.

A fresh look at linear cosmological constraints on a decaying dark matter component

Abstract: We consider a cosmological model in which a fraction fdcdm of the Dark Matter (DM) is allowed to decay in an invisible relativistic component, and compute the resulting constraints on both the decay width (or inverse lifetime) Γdcdm and fdcdm from purely gravitational arguments. We report a full derivation of the Boltzmann hierarchy, correcting a mistake in previous literature, and compute the impact of the decay—as a function of the lifetime—on the CMB and matter power spectra. From CMB only, we obtain that no more than 3.8% of the DM could have decayed in the time between recombination and today (all bounds quoted at 95% CL). We also comment on the important application of this bound to the case where primordial black holes constitute DM, a scenario notoriously difficult to constrain. For lifetimes longer than the age of the Universe, the bounds can be cast as fdcdmΓdcdm < 6.3×10−3 Gyr−1. For the first time, we also checked that degeneracies with massive neutrinos are broken when information from the large scale structure is used. Even secondary effects like CMB lensing suffice to this purpose. Decaying DM models have been invoked to solve a possible tension between low redshift astronomical measurements of σ8 and Ωm and the ones inferred by Planck. We reassess this claim finding that with the most recent BAO, HST and σ8 data extracted from the CFHT survey, the tension is only slightly reduced despite the two additional free parameters. Nonetheless, the existing tension explains why the bound on fdcdmΓdcdm loosens to fdcdmΓdcdm < 15.9×10−3 Gyr−1 when including such additional data. The bound however improves to fdcdmΓdcdm < 5.9 ×10−3 Gyr−1 if only data consistent with the CMB are included. This highlights the importance of establishing whether the tension is due to real physical effects or unaccounted systematics, for settling the reach of achievable constraints on decaying DM.

Hamiltonian analysis of higher derivative scalar-tensor theories

Abstract: We perform a Hamiltonian analysis of a large class of scalar-tensor Lagrangians which depend quadratically on the second derivatives of a scalar field. By resorting to a convenient choice of dynamical variables, we show that the Hamiltonian can be written in a very simple form, where the Hamiltonian and the momentum constraints are easily identified. In the case of degenerate Lagrangians, which include the Horndeski and beyond Horndeski quartic Lagrangians, our analysis confirms that the dimension of the physical phase space is reduced by the primary and secondary constraints due to the degeneracy, thus leading to the elimination of the dangerous Ostrogradsky ghost. We also present the Hamiltonian formulation for nondegenerate theories and find that they contain four degrees of freedom, including a ghost, as expected. We finally discuss the status of the unitary gauge from the Hamiltonian perspective.

Scale-dependent gravitational waves from a rolling axion

Abstract: We consider a model in which a pseudo-scalar field σ rolls for some e-folds during inflation, sourcing one helicity of a gauge field. These fields are only gravitationally coupled to the inflaton, and therefore produce scalar and tensor primordial perturbations only through gravitational interactions. These sourced signals are localized on modes that exit the horizon while the roll of σ is significant. We focus our study on cases in which the model can simultaneously produce (i) a large gravitational wave signal, resulting in observable B-modes of the CMB polarizations, and (ii) sufficiently small scalar perturbations, so to be in agreement with the current limits from temperature anisotropies. Different choices of parameters can instead lead to a localized and visible departure from gaussianity in the scalar sector, either at CMB or LSS scales.

Healthy degenerate theories with higher derivatives

Abstract: In the context of classical mechanics, we study the conditions under which higher-order derivative theories can evade the so-called Ostrogradsky instability. More precisely, we consider general Lagrangians with second order time derivatives, of the form L(a, a, a; qi, qi) with a = 1,⋯,n and i = 1,⋯,m. For n = 1, assuming that the qi's form a nondegenerate subsystem, we confirm that the degeneracy of the kinetic matrix eliminates the Ostrogradsky instability. The degeneracy implies, in the Hamiltonian formulation of the theory, the existence of a primary constraint, which generates a secondary constraint, thus eliminating the Ostrogradsky ghost. For n > 1, we show that, in addition to the degeneracy of the kinetic matrix, one needs to impose extra conditions to ensure the presence of a sufficient number of secondary constraints that can eliminate all the Ostrogradsky ghosts. When these conditions that ensure the disappearance of the Ostrogradsky instability are satisfied, we show that the Euler-Lagrange equations, which involve a priori higher order derivatives, can be reduced to a second order system.

Stellar equilibrium configurations of compact stars in f(R,T) theory of gravity

Abstract: In this article we study the hydrostatic equilibrium configuration of neutron stars and strange stars, whose fluid pressure is computed from the equations of state p=ωρ5/3 and p=028(ρ−4ℬ), respectively, with ω and ℬ being constants and ρ the energy density of the fluid We start by deriving the hydrostatic equilibrium equation for the f(R,T) theory of gravity, with R and T standing for the Ricci scalar and trace of the energy-momentum tensor, respectively Such an equation is a generalization of the one obtained from general relativity, and the latter can be retrieved for a certain limit of the theory For the f(R,T)=R+2λ T functional form, with λ being a constant, we find that some physical properties of the stars, such as pressure, energy density, mass and radius, are affected when λ is changed We show that for a fixed central star energy density, the mass of neutron and strange stars can increase with λ Concerning the star radius, it increases for neutron stars and it decreases for strange stars with the increment of λ Thus, in f(R,T) theory of gravity we can push the maximum mass above the observational limits This implies that the equation of state cannot be eliminated if the maximum mass within General Relativity lies below the limit given by observed pulsars

Production Regimes for Self-Interacting Dark Matter

Abstract: In the context of Self-Interacting Dark Matter as a solution for the small-scale structure problems, we consider the possibility that Dark Matter could have been produced without being in thermal equilibrium with the Standard Model bath. We discuss one by one the following various dark matter production regimes of this kind: freeze-in, reannihilation and dark freeze-out. We exemplify how these mechanisms work in the context of the particularly simple Hidden Vector Dark Matter model. In contrast to scenarios where there is thermal equilibrium with the Standard Model bath, we find two regimes which can easily satisfy all the laboratory and cosmological constraints. These are dark freeze-out with 3-to-2 annihilations and freeze-in via a light mediator. In the first regime, different temperatures in the visible and the Dark Matter sectors allow us to avoid the constraints coming from cosmic structure formation as well as the use of non-perturbative couplings to reproduce the observed relic density. For the second regime, different couplings are responsible for Dark Matter relic density and self-interactions, permitting to surpass BBN, X-ray, CMB and direct detection constraints.

Generalized Proca action for an Abelian vector field

Abstract: We revisit the most general theory for a massive vector field with derivative self-interactions, extending previous works on the subject to account for terms having trivial total derivative interactions for the longitudinal mode. In the flat spacetime (Minkowski) case, we obtain all the possible terms containing products of up to five first-order derivatives of the vector field, and provide a conjecture about higher-order terms. Rendering the metric dynamical, we covariantize the results and add all possible terms implying curvature.

Robust forecasts on fundamental physics from the foreground-obscured, gravitationally-lensed CMB polarization

Abstract: Recent results from the BICEP, Keck Array and Planck Collaborations demonstrate that Galactic foregrounds are an unavoidable obstacle in the search for evidence of inflationary gravitational waves in the cosmic microwave background (CMB) polarization. Beyond the foregrounds, the effect of lensing by intervening large-scale structure further obscures all but the strongest inflationary signals permitted by current data. With a plethora of ongoing and upcoming experiments aiming to measure these signatures, careful and self-consistent consideration of experiments' foreground- and lensing-removal capabilities is critical in obtaining credible forecasts of their performance. We investigate the capabilities of instruments such as Advanced ACTPol, BICEP3 and Keck Array, CLASS, EBEX10K, PIPER, Simons Array, SPT-3G and SPIDER, and projects as COrE+, LiteBIRD-ext, PIXIE and Stage IV, to clean contamination due to polarized synchrotron and dust from raw multi-frequency data, and remove lensing from the resulting co-added CMB maps (either using iterative CMB-only techniques or through cross-correlation with external data). Incorporating these effects, we present forecasts for the constraining power of these experiments in terms of inflationary physics, the neutrino sector, and dark energy parameters. Made publicly available through an online interface, this tool enables the next generation of CMB experiments to foreground-proof their designs, optimize their frequency coverage to maximize scientific output, and determine where cross-experimental collaboration would be most beneficial. We find that analyzing data from ground, balloon and space instruments in complementary combinations can significantly improve component separation performance, delensing, and cosmological constraints over individual datasets. In particular, we find that a combination of post-2020 ground- and space-based experiments could achieve constraints such as σ(r)~1.3×10(−)(4), σ(n(t))~0.03, σ( n(s) )~1.8×10(−)(3), σ(α(s))~1.7×10(−)(3), σ( M(ν) )~31 meV, σ( w )~0.09, σ( w(0) )~ 0.25, 0σ( wa )~ 0.5, σ( N(eff) )~0.024 and σ( Ω(k) )~1.5×10(−)(3), after component separation and iterative delensing.

Lyman-alpha forests cool warm dark matter

Abstract: The free-streaming of keV-scale particles impacts structure growth on scales that are probed by the Lyman-alpha forest of distant quasars. Using an unprecedentedly large sample of medium-resolution QSO spectra from the ninth data release of SDSS, along with a state-of-the-art set of hydrodynamical simulations to model the Lyman-alpha forest in the non-linear regime, we issue one of the tightest bounds to date, from Ly-α data alone, on pure dark matter particles: mX>4.09 keV (95% CL) for early decoupled thermal relics such as a hypothetical gravitino, and correspondingly ms>24.4 keV (95% CL) for a non-resonantly produced right-handed neutrino. This limit depends on the value on ns, and Planck measures a higher value of ns than SDSS-III/BOSS. Our bounds thus change slightly when Ly-α data are combined with CMB data from Planck 2016. The limits shift to mX>2.96 keV (95% CL) and ms>16.0 keV (95% CL). Thanks to SDSS-III data featuring smaller uncertainties and covering a larger redshift range than SDSS-I data, our bounds confirm the most stringent results established by previous works and are further at odds with a purely non-resonantly produced sterile neutrino as dark matter.

Toward (finally!) ruling out Z and Higgs mediated dark matter models

Abstract: In recent years, direct detection, indirect detection, and collider experiments have placed increasingly stringent constraints on particle dark matter, exploring much of the parameter space associated with the WIMP paradigm. In this paper, we focus on the subset of WIMP models in which the dark matter annihilates in the early universe through couplings to either the Standard Model Z or the Standard Model Higgs boson. Considering fermionic, scalar, and vector dark matter candidates within a model-independent context, we find that the overwhelming majority of these dark matter candidates are already ruled out by existing experiments. In the case of Z mediated dark matter, the only scenarios that are not currently excluded are those in which the dark matter is a fermion with an axial coupling and with a mass either within a few GeV of the Z resonance (mDM ≃ mZ/2) or greater than 200 GeV, or with a vector coupling and with mDM > 6 TeV . Several Higgs mediated scenarios are currently viable if the mass of the dark matter is near the Higgs pole (mDM ≃ mH/2). Otherwise, the only scenarios that are not excluded are those in which the dark matter is a scalar (vector) heaviermore » than 400 GeV (1160 GeV) with a Higgs portal coupling, or a fermion with a pseudoscalar (CP violating) coupling to the Standard Model Higgs boson. Furthermore, with the exception of dark matter with a purely pseudoscalar coupling to the Higgs, it is anticipated that planned direct detection experiments will probe nearly the entire range of models considered in this study.« less

Improved limits on dark matter annihilation in the Sun with the 79-string IceCube detector and implications for supersymmetry

Abstract: We present an improved event-level likelihood formalism for including neutrino telescope data in global fits to new physics. We derive limits on spin-dependent dark matter-proton scattering by employing the new formalism in a re-analysis of data from the 79-string IceCube search for dark matter annihilation in the Sun, including explicit energy information for each event. The new analysis excludes a number of models in the weak-scale minimal supersymmetric standard model (MSSM) for the first time. This work is accompanied by the public release of the 79-string IceCube data, as well as an associated computer code for applying the new likelihood to arbitrary dark matter models.

Precision measurement of the local bias of dark matter halos

Abstract: We present accurate measurements of the linear, quadratic, and cubic local bias of dark matter halos, using curved "separate universe" N-body simulations which effectively incorporate an infinite-wavelength overdensity. This can be seen as an exact implementation of the peak-background split argument. We compare the results with the linear and quadratic bias measured from the halo-matter power spectrum and bispectrum, and find good agreement. On the other hand, the standard peak-background split applied to the Sheth & Tormen (1999) and Tinker et al. (2008) halo mass functions matches the measured linear bias parameter only at the level of 10%. The prediction from the excursion set-peaks approach performs much better, which can be attributed to the stochastic moving barrier employed in the excursion set-peaks prediction. We also provide convenient fitting formulas for the nonlinear bias parameters b2(b1) and b3(b1), which work well over a range of redshifts.

Phases of New Physics in the CMB

Abstract: Fluctuations in the cosmic neutrino background are known to produce a phase shift in the acoustic peaks of the cosmic microwave background. It is through the sensitivity to this effect that the recent CMB data has provided a robust detection of free-streaming neutrinos. In this paper, we revisit the phase shift of the CMB anisotropy spectrum as a probe of new physics. The phase shift is particularly interesting because its physical origin is strongly constrained by the analytic properties of the Green's function of the gravitational potential. For adiabatic fluctuations, a phase shift requires modes that propagate faster than the speed of fluctuations in the photon-baryon plasma. This possibility is realized by free-streaming relativistic particles, such as neutrinos or other forms of dark radiation. Alternatively, a phase shift can arise from isocurvature fluctuations. We present simple models to illustrate each of these effects. We then provide observational constraints from the Planck temperature and polarization data on additional forms of radiation. We also forecast the capabilities of future CMB Stage IV experiments. Whenever possible, we give analytic interpretations of our results.

Beginning inflation in an inhomogeneous universe

Abstract: Using numerical solutions of the full Einstein field equations coupled to a scalar inflaton field in 3+1 dimensions, we study the conditions under which a universe that is initially expanding, highly inhomogeneous and dominated by gradient energy can transition to an inflationary period. If the initial scalar field variations are contained within a sufficiently flat region of the inflaton potential, and the universe is spatially flat or open on average, inflation will occur following the dilution of the gradient and kinetic energy due to expansion. This is the case even when the scale of the inhomogeneities is comparable to the initial Hubble length, and overdense regions collapse and form black holes, because underdense regions continue expanding, allowing inflation to eventually begin. This establishes that inflation can arise from highly inhomogeneous initial conditions and solve the horizon and flatness problems, at least as long as the variations in the scalar field do not include values that exceed the inflationary plateau.

Matter power spectrum and the challenge of percent accuracy

Abstract: Future galaxy surveys require one percent precision in the theoretical knowledge of the power spectrum over a large range including very nonlinear scales. While this level of accuracy is easily obtained in the linear regime with perturbation theory, it represents a serious challenge for small scales where numerical simulations are required. In this paper we quantify the precision of present-day N -body methods, identifying main potential error sources from the set-up of initial conditions to the measurement of the final power spectrum. We directly compare three widely used N -body codes, Ramses, Pkdgrav3, and Gadget3 which represent three main discretisation techniques: the particle-mesh method, the tree method, and a hybrid combination of the two. For standard run parameters, the codes agree to within one percent at k≤1 hMpc −1 and to within three percent at k≤10 hMpc −1. We also consider the bispectrum and show that the reduced bispectra agree at the sub-percent level for k≤2 hMpc −1 . In a second step, we quantify potential errors due to initial conditions, box size, and resolution using an extended suite of simulations performed with our fastest code Pkdgrav3. We demonstrate that the simulation box size should not be smaller than L=0.5 h −1 Gpc to avoid systematic finite-volume effects (while much larger boxes are required to beat down the statistical sample variance). Furthermore, a maximum particle mass of M p =10 9 h −1 M ⊙ is required to conservatively obtain one percent precision of the matter power spectrum. As a consequence, numerical simulations covering large survey volumes of upcoming missions such as DES, LSST, and Euclid will need more than a trillion particles to reproduce clustering properties at the targeted accuracy.

Evidence for dark matter interactions in cosmological precision data

Abstract: We study a two-parameter extension of the cosmological standard model ΛCDM in which cold dark matter interacts with a new form of dark radiation. The two parameters correspond to the energy density in the dark radiation fluid ΔNfluid and the interaction strength between dark matter and dark radiation. The interactions give rise to a very weak ``dark matter drag'' which damps the growth of matter density perturbations throughout radiation domination, allowing to reconcile the tension between predictions of large scale structure from the CMB and direct measurements of σ8. We perform a precision fit to Planck CMB data, BAO, large scale structure, and direct measurements of the expansion rate of the universe today. Our model lowers the χ-squared relative to ΛCDM by about 12, corresponding to a preference for non-zero dark matter drag by more than 3σ. Particle physics models which naturally produce a dark matter drag of the required form include the recently proposed non-Abelian dark matter model in which the dark radiation corresponds to massless dark gluons.