Search

The steepness of the central density profiles of dark matter (DM) in low-mass galaxy halos (e.g. dwarf galaxies) is a powerful probe of the nature of DM. We propose a novel scheme to probe the inner profiles of galaxy subhalos using stellar streams. We show that the present day morphological and dynamical properties of accreted globular cluster (GC) streams - those produced from tidal stripping of GCs that initially evolved within satellite galaxies and later merged with the Milky Way (MW) - are sensitive to the central DM density profile and mass of their parent satellites. GCs that accrete within cuspy CDM subhalos produce streams that are physically wider and dynamically hotter than streams that accrete inside cored subhalos. A first comparison of MW streams "GD-1" and "Jhelum" (likely of accreted GC origin) with our simulations indicates a preference for cored subhalos. If these results hold up in future data, the implication is that either the DM cusps were erased by baryonic feedback, or their subhalos naturally possessed cored density profiles implying DM models beyond CDM. Moreover, accreted GC streams are highly structured and exhibit complex morphological features (e.g., parallel structures and "spurs"). This implies that the accretion scenario can naturally explain the recently observed peculiarities in some of the MW streams. We also propose a novel mechanism for forming "gaps" in streams when the remnant of the parent subhalo later passes through the stream. This encounter can last a longer time (and have more of an impact) than the random encounters with DM subhalos previously considered, because the GC stream and its parent subhalo are on similar orbits with small relative velocities. Current and future surveys of the MW halo will uncover numerous faint stellar streams and provide the data needed to substantiate our preliminary tests with this new probe of DM., Comment: 24 pages, 15 figures, 1 table, accepted for publication in MNRAS

In standard stellar evolution, stars with masses ranging from approximately $150$ to $240 M_\odot$ are expected to evolve to a pair instability supernova with no black hole (BH) remnant. This evolutionary behavior leads to a predicted gap in the black hole mass function from approximately $50$ to $140 M_\odot$. Yet the LIGO and Virgo Collaborations[1] recently discovered black holes of masses $66 M_\odot$ and $85 M_\odot$ in the gravitational wave event GW190521. We propose a new method to populate the BH mass gap. If an energy source is added throughout the star in addition to nuclear fusion, it is possible for the altered evolution to avoid the complete destruction of a pair instability supernova, and instead a BH remnant is left behind. An example of an extra energy source is dark matter annihilation within the star, but our results hold more generally. We show this phenomenon by exploring the effect of adding an energy source independent of temperature and density to a $180 M_\odot$ star, using the MESA one-dimensional stellar evolution software. If $\sim50$\% of the star's energy is due to this new source, the star is capable of avoiding the pair instability entirely and evolving towards a core-collapse supernova and ultimately a BH remnant with mass $\sim 120 M_\odot$., 20 pages, 9 figures, 2 tables

Recently, the DAMA/LIBRA collaboration released updated results from their search for the annual modulation signal from Dark Matter (DM) scattering in the detector. Besides approximately doubling the exposure of the DAMA/LIBRA data set, the updated photomultiplier tubes of the experiment allow a lower recoil energy threshold of 1\,keV electron equivalent compared to the previous threshold of 2 keV electron equivalent. We study the compatibility of the observed modulation signal with DM scattering. Due to a conspiracy of multiple effects, the new data at low recoil energies is very powerful for testing the DM hypothesis. We find that canonical (isospin conserving) spin-independent DM-nucleon interactions are no longer a good fit to the observed modulation signal in the standard halo model. The canonical spin-independent case is disfavored by the new data, with best fit points of a DM mass of $\sim 8\,$GeV, disfavored by $5.2\,\sigma$, or a mass of $\sim 54\,$GeV, disfavored by $2.5\,\sigma$. Allowing for isospin violating spin independent interactions, we find a region with a good fit to the data with suppressed effective couplings to iodine for DM masses of $\sim 10\,$GeV. We also consider spin-dependent DM-nucleon interactions, which yield good fits for similar DM masses of $\sim 10\,$GeV or $\sim 45\,$GeV, Comment: 8 pages, 3 figures. v2: 10 pages, 5 figures. Updated to published DAMA/LIBRA-phase2 data, added results for spin-dependent WIMP-nucleon scattering, and additional references. v3: matches published version

Paleo-detectors are a proposed experimental technique to search for dark matter (DM). In lieu of the conventional approach of operating a tonne-scale real-time detector to search for DM-induced nuclear recoils, paleo-detectors take advantage of small samples of naturally occurring rocks on Earth that have been deep underground ($\gtrsim 5$ km), accumulating nuclear damage tracks from recoiling nuclei for $\mathcal{O}(1)$ Gyr. Modern microscopy techniques promise the capability to read out nuclear damage tracks with nanometer resolution in macroscopic samples. Thanks to their $\mathcal{O}(1)$ Gyr integration times, paleo-detectors could constitute nuclear recoil detectors with keV recoil energy thresholds and 100 kilotonne-yr exposures. This combination would allow paleo-detectors to probe DM-nucleon cross sections orders of magnitude below existing upper limits from conventional direct detection experiments. In this article, we use improved background modeling and a new spectral analysis technique to update the sensitivity forecast for paleo-detectors. We demonstrate the robustness of the sensitivity forecast to the (lack of) ancillary measurements of the age of the samples and the parameters controlling the backgrounds, systematic mismodeling of the spectral shape of the backgrounds, and the radiopurity of the mineral samples. Specifically, we demonstrate that even if the uranium concentration in paleo-detector samples is $10^{-8}$ (per weight), many orders of magnitude larger than what we expect in the most radiopure samples obtained from ultra basic rock or marine evaporite deposits, paleo-detectors could still probe DM-nucleon cross sections below current limits. For DM masses $\lesssim 10$ GeV/$c^2$, the sensitivity of paleo-detectors could still reach down all the way to the conventional neutrino floor in a Xe-based direct detection experiment., Comment: Invited contribution to Instruments "Special Issue Innovative Experimental Techniques for Direct Dark Matter Detection)". 30 pages, 5 figures, 1 table. Code available at https://github.com/sbaum90/paleoSpec and https://github.com/sbaum90/paleoSens

It is commonly assumed that the energy density of the Universe was dominated by radiation between reheating after inflation and the onset of matter domination 54,000 years later. While the abundance of light elements indicates that the Universe was radiation dominated during Big Bang Nucleosynthesis (BBN), there is scant evidence that the Universe was radiation dominated prior to BBN. It is therefore possible that the cosmological history was more complicated, with deviations from the standard radiation domination during the earliest epochs. Indeed, several interesting proposals regarding various topics such as the generation of dark matter, matter-antimatter asymmetry, gravitational waves, primordial black holes, or microhalos during a nonstandard expansion phase have been recently made. In this paper, we review various possible causes and consequences of deviations from radiation domination in the early Universe - taking place either before or after BBN - and the constraints on them, as they have been discussed in the literature during the recent years., 67 pages, 18 figures. v2: Discussion and references added. Accepted for publication in The Open Journal of Astrophysics

In the currently favored cosmological paradigm galaxies form hierarchically through the accretion of numerous satellite galaxies. Since the satellites are much less massive than the host halo, they occupy a small fraction of the volume in action space defined by the potential of the host halo. Since actions are conserved when the potential of the host halo changes adiabatically, stars from an accreted satellite are expected to remain clustered in action space as the host halo evolves. In this paper, we identify accreted satellites in three Milky Way like disk galaxies from the cosmological baryonic FIRE-2 simulations by tracking satellite galaxies through simulation snapshots. We then try to recover these satellites by applying the cluster analysis algorithm Enlink to the orbital actions of accreted star particles in the present-day snapshot. We define several metrics to quantify the success of the clustering algorithm and use these metrics to identify well-recovered and poorly-recovered satellites. We plot these satellites in the infall time-progenitor mass (or stellar mass) space, and determine the boundaries between the well-recovered and poorly-recovered satellites in these two spaces with classification tree method. The groups found by Enlink are more likely to correspond to a real satellite if they have high significance, a quantity assigned by Enlink. Since cosmological simulations predict that most stellar halos have a population of insitu stars, we test the ability of Enlink to recover satellites when the sample is contaminated by 10-50% of insitu star particles, and show that most of the satellites well-recovered by Enlink in the absence of insitu stars, stay well-recovered even with 50% contamination. We thus expect that, in the future, cluster analysis in action space will be useful in upcoming data sets (e.g. Gaia) for identifying accreted satellites in the Milky Way., Comment: 16 pages, 9 figures, fixed some typos and added two references

Chain Inflation is an alternative to slow roll inflation in which the universe undergoes a series of transitions between different vacua. The density perturbations (studied in this paper) are seeded by the probabilistic nature of tunneling rather than quantum fluctuations of the inflaton. We find the scalar power spectrum of chain inflation and show that it is fully consistent with a $\Lambda$CDM cosmology. In agreement with some of the previous literature (and disagreement with others), we show that $10^4$ phase transitions per e-fold are required in order to agree with the amplitude of Cosmic Microwave Background anisotropies within the observed range of scales. Interestingly, the amplitude of perturbations constrains chain inflation to a regime of highly unstable de Sitter spaces, which may be favorable from a quantum gravity perspective since the Swampland Conjecture on Trans-Planckian Censorship is automatically satisfied. We provide new analytic estimates for the bounce action and the tunneling rate in periodic potentials which replace the thin-wall approximation in the regime of fast tunneling. Finally, we study model implications and derive an upper limit of $\sim 10^{10}$ GeV on the axion decay constant in viable chain inflation with axions., Comment: 19 pages, 9 figures

In light of the improved sensitivities of cosmological observations, we examine the status of quasi-degenerate neutrino mass scenarios. Within the simplest extension of the standard cosmological model with massive neutrinos, we find that quasi-degenerate neutrinos are severely constrained by present cosmological data and neutrino oscillation experiments. % % We find that Planck 2018 observations of cosmic microwave background (CMB) anisotropies disfavour quasi-degenerate neutrino masses at $2.4$ Gaussian $\sigma$'s, while adding Baryon acoustic oscillations (BAO) data brings the rejection to 5.9$\sigma$'s. % The highest statistical significance with which one would be able to rule out quasi-degeneracy would arise if the sum of neutrino masses is $\Sigma m_\nu = 60$ \meV (the minimum allowed by neutrino oscillation experiments); % indeed a sensitivity of 15 meV, as expected from a combination of future cosmological probes, would further improve the rejection level up to 17$\sigma$. % We discuss the robustness of these projections with respect to assumptions on the underlying cosmological model, and also compare them with bounds from $\beta$ decay endpoint and neutrinoless double beta decay studies., Comment: 19 pages, 6 figures, 1 table

Upcoming cosmic microwave background (CMB) data can be used to explore harmonic 3-point functions that involve the B-mode component of the CMB polarization signal. We focus on bispectra describing the non-Gaussian correlation of the B-mode field and the CMB temperature anisotropies (T) and/or E-mode polarization, i.e. , , and . Such bispectra probe violations of the tensor consistency relation: the model-independent behavior of cosmological correlation functions that involve a large-wavelength tensor mode (gravitational wave). An observed violation of the tensor consistency relation would exclude a large number of inflation models. We describe a generalization of the Komatsu-Spergel-Wandelt (KSW) bispectrum estimator that allows statistical inference on this type of primordial non-Gaussianity with data of the CMB temperature and polarization anisotropies. The generalized estimator shares its statistical properties with the existing KSW estimator and retains the favorable numerical scaling with angular resolution. In this paper we derive the estimator and present a set of Fisher forecasts. We show how the forecasts scale with various experimental parameters such as lower and upper angular band-limit, relevant for e.g. the upcoming ground-based Simons Observatory experiment and proposed LiteBIRD satellite experiment. We comment on possible contaminants due to secondary cosmological and astrophysical sources., Comment: 40 pages, 7 figures

Paleo-detectors are a proposed experimental technique in which one would search for traces of recoiling nuclei in ancient minerals. Natural minerals on Earth are as old as $\mathcal{O}(1)\,$Gyr and, in many minerals, the damage tracks left by recoiling nuclei are also preserved for timescales long compared to 1 Gyr once created. Thus, even reading out relatively small target samples of order 100 g, paleo-detectors would allow one to search for very rare events thanks to the large exposure, $\varepsilon \sim 100\,{\rm g}\,{\rm Gyr} = 10^5\,{\rm t}\,{\rm yr}$. Here, we explore the potential of paleo-detectors to measure nuclear recoils induced by neutrinos from galactic core collapse supernovae. We find that they would not only allow for a direct measurement of the average core collapse supernova rate in the Milky Way, but would also contain information about the time-dependence of the local supernova rate over the past $\sim$1 Gyr. Since the supernova rate is thought to be directly proportional to the star formation rate, such a measurement would provide a determination of the local star formation history. We investigate the sensitivity of paleo-detectors to both a smooth time evolution and an enhancement of the core collapse supernova rate on relatively short timescales, as would be expected for a starburst period in the local group., Comment: 15 pages, 5 figures. code available at https://github.com/tedwards2412/SN-paleology (archived at https://doi.org/10.5281/zenodo.3066206 ) v2: matches the published version

A large experimental program is underway to extend the sensitivity of direct detection experiments, searching for interaction of Dark Matter with nuclei, down to the neutrino floor. However, such experiments are becoming increasingly difficult and costly due to the large target masses and exquisite background rejection needed for the necessary improvements in sensitivity. We investigate an alternative approach to the detection of Dark Matter-nucleon interactions: Searching for the persistent traces left by Dark Matter scattering in ancient minerals obtained from much deeper than current underground laboratories. We estimate the sensitivity of paleo-detectors, which extends far beyond current upper limits for a wide range of Dark Matter masses. The sensitivity of our proposal also far exceeds the upper limits set by Snowden-Ifft et al. more than three decades ago using ancient Mica in an approach similar to paleo-detectors., Comment: 6 pages, 2 figures. v2: significant revisions. v3: matches the published version

The first bright objects to form in the Universe might not have been "ordinary" fusion-powered stars, but "Dark Stars" (DSs) powered by the annihilation of dark matter (DM) in the form of Weakly Interacting Massive Particles (WIMPs). If discovered, DSs can provide a unique laboratory to test DM models. DSs are born with a mass of order $M_\odot$ and may grow to a few million solar masses; in this work we investigate the properties of early DSs with masses up to $\sim \! 1000 \, M_\odot$, fueled by WIMPS weighing $100$ GeV. We improve the previous implementation of the DM energy source into the stellar evolution code MESA. We show that the growth of DSs is not limited by astrophysical effects: DSs up to $\sim \! 1000 \, M_\odot$ exhibit no dynamical instabilities; DSs are not subject to mass-loss driven by super-Eddington winds. We test the assumption of previous work that the injected energy per WIMP annihilation is constant throughout the star; relaxing this assumption does not change the properties of the DSs. Furthermore, we study DS pulsations, for the first time investigating non-adiabatic pulsation modes, using the linear pulsation code GYRE. We find that acoustic modes in DSs of masses smaller than $\sim \! 200 \, M_\odot$ are excited by the $\kappa-\gamma$ and $\gamma$ mechanism in layers where hydrogen or helium is (partially) ionized. Moreover, we show that the mass-loss rates potentially induced by pulsations are negligible compared to the accretion rates., Comment: minor revisions; accepted by MNRAS; 15 pages, 6 Figures

The Event Horizon Telescope (EHT) collaboration has recently released the first image of a black hole (BH), opening a new window onto tests of general relativity in the strong field regime. In this paper, we derive constraints on the nature of M87* (the supermassive object at the centre of the galaxy M87), exploiting the fact that its shadow appears to be highly circular, and using measurements of its angular size. We first consider the simple case where M87* is assumed to be a Kerr BH. We find that the inferred circularity of M87* excludes Kerr BHs with observation angle $\theta_{\rm obs} \gtrsim 45^{\circ}$ for dimensionless rotational parameter $0.95 \lesssim a_* \leq 1$ whereas the observation angle is unbounded for $a_* \lesssim 0.9$. We then consider the possibility that M87* might be a superspinar, i.e. an object described by the Kerr solution and spinning so fast that it violates the Kerr bound by having $|a_*| > 1$. We find that, within certain regions of parameter space, the inferred circularity and size of the shadow of M87* do not exclude the possibility that this object might be a superspinar., Comment: 9 pages, 4 figures, Accepted in PRD

The rotation curve of the Milky Way is commonly used to estimate the local dark matter density $\rho_{{\rm DM},\odot}$. However, the estimates are subject to the choice of the distribution of baryons needed in this type of studies. In this work we explore several Galactic mass models that differ in the distribution of baryons and dark matter, in order to determine $\rho_{{\rm DM},\odot}$. For this purpose we analyze the precise circular velocity curve measurement of the Milky Way up to $\sim 25$ kpc from the Galactic centre obtained from Gaia DR2 [1]. We find that the estimated value of $\rho_{{\rm DM},\odot}$ stays robust to reasonable changes in the spherical dark matter halo. However, we show that $\rho_{{\rm DM},\odot}$ is affected by the choice of the model for the underlying baryonic components. In particular, we find that $\rho_{{\rm DM},\odot}$ is mostly sensitive to uncertainties in the disk components of the Galaxy. We also show that, when choosing one particular baryonic model, the estimate of $\rho_{{\rm DM},\odot}$ has an uncertainty of only about $10\%$ of its best-fit value, but this uncertainty gets much bigger when we also consider the variation of the baryonic model. In particular, the rotation curve method does not allow to exclude the presence of an additional very thin component, that can increase $\rho_{{\rm DM},\odot}$ by more than a factor of 8 (the thin disk could even be made of dark matter). Therefore, we conclude that exclusively using the rotation curve of the Galaxy is not enough to provide a robust estimate of $\rho_{{\rm DM},\odot}$. For all the models that we study without the presence of an additional thin component, our resulting estimates of the local dark matter density take values in the range $\rho_{{\rm DM},\odot} \simeq \text{0.3--0.4}\,\mathrm{GeV/cm^3}$, consistent with many of the estimates in the literature., Comment: 19 pages, 8 figures, 4 tables

The effect of large quantum fluctuations on primordial black-hole formation during an inflationary phase along a plateau is investigated. By using techniques of stochastic inflation in combination with replica field theory and the Feynman-Jensen variational method, it is non-perturbatively demonstrated that the abundance of primordial black holes is amplified by several orders of magnitude as compared to the classical computation., Comment: 5 pages, 3 figures; v2: revision to match published version

Direct detection experiments have set increasingly stringent limits on the cross section for spin-independent dark matter-nucleon interactions. In obtaining such limits, experiments primarily assume the standard halo model (SHM) as the distribution of dark matter in our Milky Way. Three astrophysical parameters are required to define the SHM: the local dark matter escape velocity, the local dark matter density and the circular velocity of the sun around the center of the galaxy. This paper studies the effect of the uncertainties in these three astrophysical parameters on the XENON1T exclusion limits using the publicly available DDCalc code. We compare limits obtained using the widely assumed escape velocity from the RAVE survey and the newly calculated escape velocity by Monari $et$ $al.$ using Gaia data. Our study finds that the astrophysical uncertainties are dominated by the uncertainty in the escape velocity (independent of the best fit value) at dark matter masses below 6 GeV and can lead to a variation of nearly 6 orders of magnitude in the exclusion limits at 4 GeV. Above a WIMP mass of 6 GeV, the uncertainty becomes dominated by the local dark matter density, leading to uncertainties of factors of $\sim$10 (3) at 6 (15) GeV WIMP mass in the exclusion limits. Additionally, this work finds that the updated best fit value for the escape velocity based on Gaia data leads to only very minor changes to the effects of the astrophysical uncertainties on the XENON1T exclusion limits., Comment: 12 pages, 3 figures; author added, references added, minor revisions; title changed, as published in JCAP with minor revisions and references added

One of the most powerful cosmological datasets when it comes to constraining neutrino masses is represented by galaxy power spectrum measurements, $P_{gg}(k)$. The constraining power of $P_{gg}(k)$ is however severely limited by uncertainties in the modeling of the scale-dependent galaxy bias $b(k)$. In this Letter we present a new method to constrain $b(k)$ by using the cross-correlation between the Cosmic Microwave Background (CMB) lensing signal and galaxy maps ($C_\ell^{\rm \kappa g}$) using a simple but theoretically well-motivated parametrization for $b(k)$. We apply the method using $C_\ell^{\rm \kappa g}$ measured by cross-correlating Planck lensing maps and the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 11 (DR11) CMASS galaxy sample, and $P_{gg}(k)$ measured from the BOSS DR12 CMASS sample. We detect a non-zero scale-dependence at moderate significance, which suggests that a proper modeling of $b(k)$ is necessary in order to reduce the impact of non-linearities and minimize the corresponding systematics. The accomplished increase in constraining power of $P_{gg}(k)$ is demonstrated by determining a 95% C.L. upper bound on the sum of the three active neutrino masses $M_{\nu}$ of $M_{\nu}, Comment: 11 pages, 3 figures

International audience; We describe the construction and characterization of the 280 GHz bolometric focal plane units (FPUs) to be deployed on the second flight of the balloon-borne Spider instrument. These FPUs are vital to Spider’s primary science goal of detecting or placing an upper limit on the amplitude of the primordial gravitational wave signature in the cosmic microwave background (CMB) by constraining the B-mode contamination in the CMB from Galactic dust emission. Each 280 GHz focal plane contains a $16 \times 16$ grid of corrugated silicon feedhorns coupled to an array of aluminum–manganese transition-edge sensor (TES) bolometers fabricated on 150 mm diameter substrates. In total, the three 280 GHz FPUs contain 1530 polarization-sensitive bolometers (765 spatial pixels) optimized for the low loading environment in flight and read out by time-division SQUID multiplexing. In this paper, we describe the mechanical, thermal, and magnetic shielding architecture of the focal planes and present cryogenic measurements which characterize yield and the uniformity of several bolometer parameters. The assembled FPUs have high yields, with one array as high as 95% including defects from wiring and readout. We demonstrate high uniformity in device parameters, finding the median saturation power for each TES array to be $\sim $ 3 pW at 300 mK with a less than 6% variation across each array at $1\sigma $ . These focal planes will be deployed alongside the 95 and 150 GHz telescopes in the Spider-2 instrument, slated to fly from McMurdo Station in Antarctica in December 2018.

The pseudoscalar resonance or "A-funnel" in the Minimal Supersymmetric Standard Model~(MSSM) is a widely studied framework for explaining dark matter that can yield interesting indirect detection and collider signals. The well-known Galactic Center excess (GCE) at GeV energies in the gamma ray spectrum, consistent with annihilation of a $\lesssim 40$ GeV dark matter particle, has more recently been shown to be compatible with significantly heavier masses following reanalysis of the background. In this paper, we explore the LHC and direct detection implications of interpreting the GCE in this extended mass window within the MSSM A-funnel framework. We find that compatibility with relic density, signal strength, collider constraints, and Higgs data can be simultaneously achieved with appropriate parameter choices. The compatible regions give very sharp predictions of 200-600 GeV CP-odd/even Higgs bosons at low tan$\beta$ at the LHC and spin-independent cross sections $\approx 10^{-11}$ pb at direct detection experiments. Regardless of consistency with the GCE, this study serves as a useful template of the strong correlations between indirect, direct, and LHC signatures of the MSSM A-funnel region., Comment: 32 pages and 9 figures

Paleo-detectors are a recently proposed method for the direct detection of Dark Matter (DM). In such detectors, one would search for the persistent damage features left by DM--nucleus interactions in ancient minerals. Initial sensitivity projections have shown that paleo-detectors could probe much of the remaining Weakly Interacting Massive Particle (WIMP) parameter space. In this paper, we improve upon the cut-and-count approach previously used to estimate the sensitivity by performing a full spectral analysis of the background- and DM-induced signal spectra. We consider two scenarios for the systematic errors on the background spectra: i) systematic errors on the normalization only, and ii) systematic errors on the shape of the backgrounds. We find that the projected sensitivity is rather robust to imperfect knowledge of the backgrounds. Finally, we study how well the parameters of the true WIMP model could be reconstructed in the hypothetical case of a WIMP discovery., 14 pages, 5 figures, code available at https://github.com/tedwards2412/paleo_detectors/ . v2: Added additional analysis theory details, matches version published in PRD

We explore cosmological constraints on the sum of the three active neutrino masses M-v in the context of dynamical dark energy (DDE) models with equation of state (EoS) parametrized as a function o ...

A significant fraction of cosmological dark matter can be formed by very dense macroscopic objects, for example primordial black holes. Gravitational waves offer a promising way to probe these kinds of dark-matter candidates, in a parameter space region that is relatively untested by electromagnetic observations. In this work we consider an ensemble of macroscopic dark matter with masses in the range $10^{-13}$ - $1\,M_{\odot}$ orbiting a super-massive black hole. While the strain produced by an individual dark-matter particle will be very small, gravitational waves emitted by a large number of such objects will add incoherently and produce a stochastic gravitational-wave background. We show that LISA can be a formidable machine for detecting the stochastic background of such objects orbiting the black hole in the centre of the Milky Way, Sgr\.${\rm A}^{\!*}$, if a dark-matter spike of the type originally predicted by Gondolo and Silk forms near the central black hole., Comment: v2: Update to match published version. arXiv admin note: text overlap with arXiv:1705.10361

The Standard Model Higgs boson, which has previously been shown to develop an effective vacuum expectation value during inflation, can give rise to large particle masses during inflation and reheating, leading to temporary blocking of the reheating process and a lower reheat temperature after inflation. We study the effects on the multiple stages of reheating: resonant particle production (preheating) as well as perturbative decays from coherent oscillations of the inflaton field. Specifically, we study both the cases of the inflaton coupling to Standard Model fermions through Yukawa interactions as well as to Abelian gauge fields through a Chern-Simons term. We find that, in the case of perturbative inflaton decay to SM fermions, reheating can be delayed due to Higgs blocking and the reheat temperature can decrease by up to an order of magnitude. In the case of gauge-reheating, Higgs-generated masses of the gauge fields can suppress preheating even for large inflaton-gauge couplings. In extreme cases, preheating can be shut down completely and must be substituted by perturbative decay as the dominant reheating channel. Finally, we discuss the distribution of reheat temperatures in different Hubble patches, arising from the stochastic nature of the Higgs VEV during inflation and its implications for the generation of both adiabatic and isocurvature fluctuations., Comment: 23 pages, 9 figures. Matches the version published on JCAP

We explore paleo-detectors as an approach to the direct detection of Weakly Interacting Massive Particle (WIMP) dark matter radically different from conventional detectors. Instead of instrumenting a (large) target mass in a laboratory in order to observe WIMP-induced nuclear recoils in real time, the approach is to examine ancient minerals for traces of WIMP-nucleus interactions recorded over timescales as large as 1 Gyr. Here, we discuss the paleo-detector proposal in detail, including background sources and possible target materials. In order to suppress backgrounds induced by radioactive contaminants such as uranium, we propose to use minerals found in marine evaporites or in ultra-basic rocks. We estimate the sensitivity of paleo-detectors to spin-independent and spin-dependent WIMP-nucleus interactions. The sensitivity to low-mass WIMPs with masses $m_\chi \lesssim 10$ GeV extends to WIMP-nucleon cross sections many orders of magnitude smaller than current upper limits. For heavier WIMPs with masses $m_\chi \gtrsim 30$ GeV cross sections a factor of a few to $\sim 100$ smaller than current upper limits can be probed by paleo-detectors., Comment: 22 pages, 7 figures, 1 table. v2: matches the published version

Axion stars are hypothetical objects formed of axions, obtained as localized and coherently oscillating solutions to their classical equation of motion. Depending on the value of the field amplitude at the core $|\theta_0| \equiv |\theta(r=0)|$, the equilibrium of the system arises from the balance of the kinetic pressure and either self-gravity or axion self-interactions. Starting from a general relativistic framework, we obtain the set of equations describing the configuration of the axion star, which we solve as a function of $|\theta_0|$. For small $|\theta_0| \lesssim 1$, we reproduce results previously obtained in the literature, and we provide arguments for the stability of such configurations in terms of first principles. We compare qualitative analytical results with a numerical calculation. For large amplitudes $|\theta_0| \gtrsim 1$, the axion field probes the full non-harmonic QCD chiral potential and the axion star enters the {\it dense} branch. Our numerical solutions show that in this latter regime the axions are relativistic, and that one should not use a single frequency approximation, as previously applied in the literature. We employ a multi-harmonic expansion to solve the relativistic equation for the axion field in the star, and demonstrate that higher modes cannot be neglected in the dense regime. We interpret the solutions in the dense regime as pseudo-breathers, and show that the life-time of such configurations is much smaller than any cosmological time scale., Comment: 11 pages, 4 figures. v2: added references, matches published version

The holy grail of physics has been to merge each of its fundamental branches into a unified "theory of everything" that would explain the functioning and existence of the universe. The last step toward this goal is to reconcile general relativity with the principles of quantum mechanics, a quest that has thus far eluded physicists. Will physics ever be able to develop an all-encompassing theory, or should we simply acknowledge that science will always have inherent limitations as to what can be known? Should new theories be validated solely on the basis of calculations that can never be empirically tested? Can we ever truly grasp the implications of modern physics when the basic laws of nature do not always operate according to our standard paradigms? These and other questions are discussed in this paper.

If dark matter is composed of weakly interacting particles, Earth's orbital motion may induce a small annual variation in the rate at which these particles interact in a terrestrial detector. The DAMA collaboration has identified at a 9.3$\sigma$ confidence level such an annual modulation in their event rate over two detector iterations, DAMA/NaI and DAMA/LIBRA, each with $\sim7$ years of observations. We statistically examine the time dependence of the modulation amplitudes, which "by eye" appear to be decreasing with time in certain energy ranges. We perform a chi-squared goodness of fit test of the average modulation amplitudes measured\ by the two detector iterations which rejects the hypothesis of a consistent modulation amplitude at greater than 80\%, 96\%, and 99.6\% for the 2--4~keVee, 2--5~keVee and 2--6~keVee energy ranges, respectively. We also find that among the 14 annual cycles there are three $\gtrsim 3\sigma$ departures from the average in the 5-6~keVee energy range. In addition, we examined several phenomenological models for the time dependence of the modulation amplitude. Using a maximum likelihood test, we find that descriptions of the modulation amplitude as decreasing with time are preferred over a constant modulation amplitude at anywhere between 1$\sigma$ and 3$\sigma$, depending on the phenomenological model for the time dependence and the signal energy range considered. A time dependent modulation amplitude is not expected for a dark matter signal, at least for dark matter halo morphologies consistent with the DAMA signal. New data from DAMA/LIBRA--phase2 will certainly aid in determining whether any apparent time dependence is a real effect or a statistical fluctuation., Comment: 13 pages, 1 figure

We present a new upper limit on CMB circular polarization from the 2015 flight of SPIDER, a balloon-borne telescope designed to search for $B$-mode linear polarization from cosmic inflation. Although the level of circular polarization in the CMB is predicted to be very small, experimental limits provide a valuable test of the underlying models. By exploiting the non-zero circular-to-linear polarization coupling of the HWP polarization modulators, data from SPIDER's 2015 Antarctic flight provide a constraint on Stokes $V$ at 95 and 150 GHz from $33, Comment: 9 pages, 4 figures, 5 tables - Updated to version accepted for publication in ApJ, modified acknowledgements

We study the collider phenomenology of the extended Higgs sector of the next-to-minimal supersymmetric Standard Model (NMSSM). The region of NMSSM parameter space favored by a 125 GeV SM-like Higgs ...

Constraints on primordial black holes in the range $10^{-18} M_{\odot}$ to $10^{3} M_{\odot}$ are reevaluated for a general class of extended mass functions. Whereas previous work has assumed that PBHs are produced with one single mass, instead there is expected to be a range of masses even in the case of production from a single mechanism; constraints therefore change from previous literature. Although tightly constrained in the majority of cases, it is shown that, even under conservative assumptions, primordial black holes in the mass range $10^{-10} M_{\odot}$ to $10^{-8} M_{\odot}$ could still constitute the entirety of the dark matter. This stresses both the importance for a comprehensive reevaluation of all respective constraints that have previously been evaluated only for a monochromatic mass function, and the need to obtain more constraints in the allowed mass range., Comment: 5 pages, 2 figures; v2: references added, minor changes

Weakly interacting massive particles (WIMPs), which are among the best motivated dark matter (DM) candidates, could make up all or only a fraction of the total DM budget. We consider a scenario in which WIMPs are a subdominant DM component; such a scenario would affect both current direct and indirect bounds on the WIMP-nucleon scattering cross section. In this paper we focus on indirect searches for the neutrino flux produced by annihilation of subdominant WIMPs captured by the Sun or the Earth via either spin-dependent or spin-independent scattering. We derive the annihilation rate and the expected neutrino flux at neutrino observatories. In our computation, we include an updated chemical composition of the Earth with respect to the previous literature, leading to an increase of the Earth's capture rate for spin-dependent scattering by a factor of 3. Results are compared with current bounds from Super-Kamiokande and IceCube. We discuss the scaling of bounds from both direct and indirect detection methods with the WIMP abundance.

Over the past few decades, a consensus picture has emerged in which roughly a quarter of the universe consists of dark matter. I begin with a review of the observational evidence for the existence of dark matter: rotation curves of galaxies, gravitational lensing measurements, hot gas in clusters, galaxy formation, primordial nucleosynthesis and cosmic microwave background observations. Then I discuss a number of anomalous signals in a variety of data sets that may point to discovery, though all of them are controversial. The annual modulation in the DAMA detector and/or the gamma-ray excess seen in the Fermi Gamma Ray Space Telescope from the Galactic Center could be due to WIMPs; a 3.5 keV X-ray line from multiple sources could be due to sterile neutrinos; or the 511 keV line in INTEGRAL data could be due to MeV dark matter. All of these would require further confirmation in other experiments or data sets to be proven correct. In addition, a new line of research on dark stars is presented, which suggests that the first stars to exist in the universe were powered by dark matter heating rather than by fusion: the observational possibility of discovering dark matter in this way is discussed., Proceedings of 14th Marcel Grossman Meeting, MG14, University of Rome "La Sapienza", Rome, July 2015

We develop a method to forecast the outcome of the LHC Run 3 based on the hypothetical detection of $\mathcal{O}(100)$ signal events at XENONnT. Our method relies on a systematic classification of renormalisable single-mediator models for dark matter-quark interactions, and is valid for dark matter candidates of spin less than or equal to one. Applying our method to simulated data, we find that at the end of the LHC Run 3 only two mutually exclusive scenarios would be compatible with the detection of $\mathcal{O}(100)$ signal events at XENONnT. In a first scenario, the energy distribution of the signal events is featureless, as for canonical spin-independent interactions. In this case, if a mono-jet signal is detected at the LHC, dark matter must have spin 1/2 and interact with nucleons through a unique velocity-dependent operator. If a mono-jet signal is not detected, dark matter interacts with nucleons through canonical spin-independent interactions. In a second scenario, the spectral distribution of the signal events exhibits a bump at non zero recoil energies. In this second case, a mono-jet signal can be detected at the LHC Run 3, dark matter must have spin 1/2 and interact with nucleons through a unique momentum-dependent operator. We therefore conclude that the observation of $\mathcal{O}(100)$ signal events at XENONnT combined with the detection, or the lack of detection, of a mono-jet signal at the LHC Run 3 would significantly narrow the range of possible dark matter-nucleon interactions. As we argued above, it can also provide key information on the dark matter particle spin., Comment: 17 pages, 8 figures, updated operator coefficients and figures, version accepted by PRD

Tidally disrupted globular cluster streams are usually observed, and therefore perceived, as narrow, linear and one-dimensional structures in the 6D phase-space. Here we show that the GD-1 stellar stream ($\approx$ 30 pc wide), which is the tidal debris of a disrupted globular cluster, possesses a secondary diffuse and extended stellar component ($\approx$ 100 pc wide) around it, detected at >5$\sigma$ confidence level. Similar morphological properties are seen in synthetic streams that are produced from star clusters that are formed within dark matter sub-halos and then accrete onto a massive host galaxy. This lends credence to the idea that the progenitor of the highly retrograde GD-1 stream was originally formed outside of the Milky Way in a now defunct dark satellite galaxy. We deem that in future studies, this newly found $cocoon$ component may serve as a structural hallmark to distinguish between the in-situ and ex-situ (accreted) formed globular cluster streams., Comment: 11 pages, 10 figures, accepted for publication in The Astrophysical Journal

We present a novel approach to derive constraints on neutrino masses from cosmological data, while taking into account our ignorance of the neutrino mass ordering. We derive constraints from a combination of current and future cosmological datasets on the total neutrino mass $M_\nu$ and on the mass fractions carried by each of the mass eigenstates, after marginalizing over the (unknown) neutrino mass ordering, either normal (NH) or inverted (IH). The bounds take therefore into account the uncertainty related to our ignorance of the mass hierarchy. This novel approach is carried out in the framework of Bayesian analysis of a typical hierarchical problem. In this context, the choice of the neutrino mass ordering is modeled via the discrete hyperparameter $h_{type}$. The preference for either the NH or the IH scenarios is then encoded in the posterior distribution of $h_{type}$ itself. Current CMB measurements assign equal odds to the two hierarchies, and are thus unable to distinguish between them. However, after the addition of BAO measurements, a weak preference for NH appears, with odds of 4:3 from Planck temperature and large-scale polarization in combination with BAO (3:2 if small-scale polarization is also included). Forecasts suggest that the combination of upcoming CMB (COrE) and BAO surveys (DESI) may determine the neutrino mass hierarchy at a high statistical significance if the mass is very close to the minimal value allowed by oscillations, as for NH and $M_\nu=0.06$ eV there is a 9:1 preference of NH vs IH. On the contrary, if $M_\nu$ is of the order of 0.1 eV or larger, even future cosmological observations will be inconclusive. The unbiased limit on $M_\nu$ we obtain with this innovative statistical strategy is crucial for ongoing and planned neutrinoless double beta decay searches., Comment: 19 pages, 7 figures, 3 tables. Abstract abridged. Comments welcome. Added discussion of the results obtained with a log prior on the lightest mass. Matching published version in PLB

We study the impact of new metallicity measurements, from solar wind data, on the solar model. The "solar modelling problem" refers to the persisting discrepancy between helioseismological observations and predictions of solar models computed implementing state-of-the-art photospheric abundances. We critically reassess the problem, in particular considering the new set of abundances of von Steiger \& Zurbuchen 2016, determined through the \textit{in situ} collection of solar wind samples from polar coronal holes. This new set of abundances indicates a solar metallicity $Z_{\odot} \geq 0.0196 \pm 0.0014$, significantly higher than the currently established value. The new values hint at an abundance of volatile elements (i.e. C, N, O, Ne) close to previous results of Grevesse \& Sauval 1998, whereas the abundance of refractory elements (i.e. Mg, Si, S, Fe) is considerably increased. Using the Linear Solar Model formalism, we determine the variation of helioseismological observables in response to the changes in elemental abundances, in order to explore the consistency of these new measurements with constraints from helioseismology. We find that, for observables that are particularly sensitive to the abundance of volatile elements, in particular the radius of the convective zone boundary (CZB) and the sound speed around the radius of CZB, improved agreement over previous models is obtained. Conversely, the high abundance of refractories correlates with a higher core temperature, resulting in an overproduction of neutrinos and a huge increase in the surface Helium abundance. We conclude that the "solar modelling problem" remains unsolved., 10 pages, 7 figures, accepted for publication in ApJ

The spatial and velocity distributions of dark matter particles in the Milky Way Halo affect the signals expected to be observed in searches for dark matter. Results from direct detection experiments are often analyzed assuming a simple isothermal distribution of dark matter, the Standard Halo Model (SHM). Yet there has been skepticism regarding the validity of this simple model due to the complicated gravitational collapse and merger history of actual galaxies. In this paper we compare the SHM to the results of cosmological hydrodynamical simulations of galaxy formation to investigate whether or not the SHM is a good representation of the true WIMP distribution in the analysis of direct detection data. We examine two Milky Way-like galaxies from the MaGICC cosmological simulations (a) with dark matter only and (b) with baryonic physics included. The inclusion of baryons drives the shape of the DM halo to become more spherical and makes the velocity distribution of dark matter particles less anisotropic especially at large heliocentric velocities, thereby making the SHM a better fit. We also note that we do not find a significant disk-like rotating dark matter component in either of the two galaxy halos with baryons that we examine, suggesting that dark disks are not a generic prediction of cosmological hydrodynamical simulations. We conclude that in the Solar neighborhood, the SHM is in fact a good approximation to the true dark matter distribution in these cosmological simulations (with baryons) which are reasonable representations of the Milky Way, and hence can also be used for the purpose of dark matter direct detection calculations., Minor changes to match JCAP version. 21 pages, 9 figures

We study the impact of assumptions about neutrino properties on the estimation of inflationary parameters from cosmological data, with a specific focus on the allowed contours in the $n_s/r$ plane. We study the following neutrino properties: (i) the total neutrino mass $ M_\nu =\sum_i m_i$; (ii) the number of relativistic degrees of freedom $N_{eff}$; and (iii) the neutrino hierarchy: whereas previous literature assumed 3 degenerate neutrino masses or two massless neutrino species (that do not match neutrino oscillation data), we study the cases of normal and inverted hierarchy. Our basic result is that these three neutrino properties induce $< 1 \sigma$ shift of the probability contours in the $n_s/r$ plane with both current or upcoming data. We find that the choice of neutrino hierarchy has a negligible impact. However, the minimal cutoff on the total neutrino mass $M_{\nu,{min}}=0 $ that accompanies previous works using the degenerate hierarchy does introduce biases in the $n_s/r$ plane and should be replaced by $M_{\nu,min}= 0.059$ eV as required by oscillation data. Using current CMB data from Planck and Bicep/Keck, marginalizing over $ M_\nu$ and over $r$ can lead to a shift in the mean value of $n_s$ of $\sim0.3\sigma$ towards lower values. However, once BAO measurements are included, the standard contours in the $n_s/r$ plane are basically reproduced. Larger shifts of the contours in the $n_s/r$ plane (up to 0.8$\sigma$) arise for nonstandard values of $N_{eff}$. We also provide forecasts for the future CMB experiments COrE and Stage-IV and show that the incomplete knowledge of neutrino properties, taken into account by a marginalization over $M_\nu$, could induce a shift of $\sim0.4\sigma$ towards lower values in the determination of $n_s$ (or a $\sim 0.8\sigma$ shift if one marginalizes over $N_{eff}$). Comparison to specific inflationary models is shown., Comment: 24 pages, 12 figures; abstract abridged

We study the capability of theJames Webb Space Telescope (JWST) to detect supermassive dark stars (SMDSs). If the first stars are powered by dark matter (DM) heating in triaxial DM haloes, they may grow to be very large (>10 6 M � ) and very bright (>10 9 L � ). These SMDSs would be visible in deep imaging with JWST and even Hubble Space Telescope (HST). We use sensitivity limits from previous HST surveys to place bounds on the numbers of SMDSs that may be detected in future JWST imaging surveys. We showed that SMDS in the mass range 10 6 -10 7 Mare bright enough to be detected in all the wavelength bands of the NIRCam on JWST (but not in the less sensitive MIRI camera at higher wavelengths). If SMDSs exist at z ∼ 10, 12 and 14, they will be detectable as J-, H -o rK-band dropouts, respectively. With a total survey area of 150 arcmin 2 (assuming a multiyear deep parallel survey with JWST), we find that typically the number of 10 6 MSMDSs found as H -o rK-band dropouts is ∼10 5 f SMDS, where the fraction of early DM haloes hosting DS is likely to be small, f SMDS � 1. If the SMDS survive down to z = 10 where HST bounds apply, then the observable number of SMDSs as H -o rK-band dropouts with JWST is ∼1-30. While individual SMDS are bright enough to be detected by JWST, standard Population III stars (without DM annihilation) are not, and would only be detected in first galaxies with total stellar masses of 10 6 -10 8 M� . Differentiating first galaxies at z > 10 from SMDSs would be possible with spectroscopy: the SMDS (which are too cool produce significant nebular emission) will have only absorption lines, while the galaxies are likely to produce emission lines as well. Of particular interest would be the He II emission lines at λ ∼ 1.6 μ ma s well as Hα lines which would be signatures of early galaxies rather than SMDSs. The detection of SMDSs with JWST would not only provide alternative evidence for weakly interacting massive particles, but also provide a possible pathway for the formation of massive (10 4 -10 6 M� ) seeds for the formation of supermassive black holes that

10 to 10^5 solar mass black holes with dark matter spikes that formed in early minihalos and still exist in our Milky Way Galaxy today are examined in light of recent data from the Fermi Gamma-Ray Space Telescope (FGST). The dark matter spikes surrounding black holes in our Galaxy are sites of significant dark matter annihilation. We examine the signatures of annihilations into gamma-rays, electrons and positrons, and neutrinos. We find that some significant fraction of the point sources detected by FGST might be due to dark matter annihilation near black holes in our Galaxy. We obtain limits on the properties of dark matter annihilations in the spikes using the information in the FGST First Source Catalog as well as the diffuse gamma-ray flux measured by FGST. We determine the maximum fraction of high redshift minihalos that could have hosted the formation of the first generation of stars and, subsequently, their black hole remnants. The strength of the limits depends on the choice of annihilation channel and black hole mass; limits are strongest for the heaviest black holes and annhilation to $b \bar{b}$ and $W^+W^-$ final states. The larger black holes considered in this paper may arise as the remnants of Dark Stars after the dark matter fuel is exhausted and thermonuclear burning runs its course; thus FGST observations may be used to constrain the properties of Dark Stars. Additionally, we comment on the excess positron flux found by PAMELA and its possible interpretation in terms of dark matter annihilation around these black hole spikes., 34 pages, 11 figures. v2: typos corrected, references added. v3: updated to match published version

The first stars in the universe form inside $\sim 10^6 M_\odot$ dark matter (DM) haloes whose initial density profiles are laid down by gravitational collapse in hierarchical structure formation scenarios. During the formation of the first stars in the universe, the baryonic infall compresses the dark matter further. The resultant dark matter density is presented here, using an algorithm originally developed by Young to calculate changes to the profile as the result of adiabatic infall in a spherical halo model; the Young prescription takes into account the non-circular motions of halo particles. The density profiles obtained in this way are found to be within a factor of two of those obtained using the simple adiabatic contraction prescription of Blumenthal et al. Our results hold regardless of the nature of the dark matter or its interactions and rely merely on gravity. If the dark matter consists of weakly interacting massive particles, which are their own antiparticles, their densities are high enough that their annihilation in the first protostars can indeed provide an important heat source and prevent the collapse all the way to fusion. In short, a ``Dark Star'' phase of stellar evolution, powered by DM annihilation, may indeed describe the first stars in the universe., 14 pages, 3 figures, and 1 table

Over the past decade, a consensus picture has emerged in which roughly a quarter of the universe consists of dark matter. The observational evidence for the existence of dark matter is reviewed: rotation curves of galaxies, weak lensing measurements, hot gas in clusters, primordial nucleosynthesis and microwave background experiments. In addition, a new line of research on Dark Stars is presented, which suggests that the first stars to exist in the universe were powered by dark matter heating rather than by fusion: the observational possibilities of discovering dark matter in this way are discussed., Comment: 14 pages, 7 figures, Conference Proceeding for "Dark Matter and Dark Energy" in Lyon, France, July 2008

Dark Stars are the very first phase of stellar evolution in the history of the universe: the first stars to form (typically at redshifts $z \sim 10-50$) are powered by heating from dark matter (DM) annihilation instead of fusion (if the DM is made of particles which are their own antiparticles). We find equilibrium polytropic configurations for these stars; we start from the time DM heating becomes important ($M \sim 1-10 M_\odot$) and build up the star via accretion up to 1000 M$_\odot$. The dark stars, with an assumed particle mass of 100 GeV, are found to have luminosities of a few times $10^6$ L$_\odot$, surface temperatures of 4000--10,000 K, radii $\sim 10^{14}$ cm, lifetimes of at least $ 0.5$ Myr, and are predicted to show lines of atomic and molecular hydrogen. Dark stars look quite different from standard metal-free stars without DM heating: they are far more massive (e.g. $\sim 800 M_\odot$ for 100 GeV WIMPs), cooler, and larger, and can be distinguished in future observations, possibly even by JWST or TMT.

The first phase of stellar evolution in the history of the universe may be Dark Stars, powered by dark matter heating rather than by fusion. Weakly interacting massive particles, which are their own antiparticles, can annihilate and provide an important heat source for the first stars in the the universe. This talk presents the story of these Dark Stars. We make predictions that the first stars are very massive ($\sim 800 M_\odot$), cool (6000 K), bright ($\sim 10^6 L_\odot$), long-lived ($\sim 10^6$ years), and probable precursors to (otherwise unexplained) supermassive black holes. Later, once the initial DM fuel runs out and fusion sets in, DM annihilation can predominate again if the scattering cross section is strong enough, so that a Dark Star is born again., 5 pages, Conference Proceeding for IAU Symposium 255: Low-Metallicity Star Formaion: From the First Stars to Dwarf Galaxies