Showing 1–50 of 217 results for author: Freese, K

Minimal Dark Matter Freeze-in with Low Reheating Temperatures and Implications for Direct Detection

Authors: Kimberly K. Boddy, Katherine Freese, Gabriele Montefalcone, Barmak Shams Es Haghi

Abstract: We investigate the influence of the reheating temperature of the visible sector on the freeze-in dark matter (DM) benchmark model for direct detection experiments, where DM production is mediated by an ultralight dark photon. Here we consider a new regime for this benchmark: we take the initial temperature of the thermal Standard Model (SM) bath to be below the DM mass. The production rate from th… ▽ More We investigate the influence of the reheating temperature of the visible sector on the freeze-in dark matter (DM) benchmark model for direct detection experiments, where DM production is mediated by an ultralight dark photon. Here we consider a new regime for this benchmark: we take the initial temperature of the thermal Standard Model (SM) bath to be below the DM mass. The production rate from the SM bath is drastically reduced due to Boltzmann suppression, necessitating a significant increase in the portal coupling between DM and the SM to match the observed relic DM abundance. This enhancement in coupling strength increases the predicted DM-electron scattering cross section, making freeze-in DM more accessible to current direct detection experiments. △ Less

Submitted 9 May, 2024; originally announced May 2024.

Comments: 5 pages, 3 figures

Report number: UTWI-16-2024, NORDITA-2024-016

Mineral Detection of Neutrinos and Dark Matter 2024. Proceedings

Authors: Sebastian Baum, Patrick Huber, Patrick Stengel, Natsue Abe, Daniel G. Ang, Lorenzo Apollonio, Gabriela R. Araujo, Levente Balogh, Pranshu Bhaumik Yilda Boukhtouchen, Joseph Bramante, Lorenzo Caccianiga, Andrew Calabrese-Day, Qing Chang, Juan I. Collar, Reza Ebadi, Alexey Elykov, Katherine Freese, Audrey Fung, Claudio Galelli, Arianna E. Gleason, Mariano Guerrero Perez, Janina Hakenmüller, Takeshi Hanyu, Noriko Hasebe, Shigenobu Hirose , et al. (35 additional authors not shown)

Abstract: The second "Mineral Detection of Neutrinos and Dark Matter" (MDvDM'24) meeting was held January 8-11, 2024 in Arlington, VA, USA, hosted by Virginia Tech's Center for Neutrino Physics. This document collects contributions from this workshop, providing an overview of activities in the field. MDvDM'24 was the second topical workshop dedicated to the emerging field of mineral detection of neutrinos a… ▽ More The second "Mineral Detection of Neutrinos and Dark Matter" (MDvDM'24) meeting was held January 8-11, 2024 in Arlington, VA, USA, hosted by Virginia Tech's Center for Neutrino Physics. This document collects contributions from this workshop, providing an overview of activities in the field. MDvDM'24 was the second topical workshop dedicated to the emerging field of mineral detection of neutrinos and dark matter, following a meeting hosted by IFPU in Trieste, Italy in October 2022. Mineral detectors have been proposed for a wide variety of applications, including searching for dark matter, measuring various fluxes of astrophysical neutrinos over gigayear timescales, monitoring nuclear reactors, and nuclear disarmament protocols; both as paleo-detectors using natural minerals that could have recorded the traces of nuclear recoils for timescales as long as a billion years and as detectors recording nuclear recoil events on laboratory timescales using natural or artificial minerals. Contributions to this proceedings discuss the vast physics potential, the progress in experimental studies, and the numerous challenges lying ahead on the path towards mineral detection. These include a better understanding of the formation and annealing of recoil defects in crystals; identifying the best classes of minerals and, for paleo-detectors, understanding their geology; modeling and control of the relevant backgrounds; developing, combining, and scaling up imaging and data analysis techniques; and many others. During the last years, MDvDM has grown rapidly and gained attention. Small-scale experimental efforts focused on establishing various microscopic readout techniques are underway at institutions in North America, Europe and Asia. We are looking ahead to an exciting future full of challenges to overcome, surprises to be encountered, and discoveries lying ahead of us. △ Less

Submitted 2 May, 2024; originally announced May 2024.

Comments: Summary and proceedings of the MDvDM'24 conference, Jan 8-11 2024

Dark Matter production during Warm Inflation via Freeze-In

Authors: Katherine Freese, Gabriele Montefalcone, Barmak Shams Es Haghi

Abstract: We present a novel perspective on the role of inflation in the production of Dark Matter (DM). Specifically, we explore the DM production during Warm Inflation via ultraviolet Freeze-In (WIFI). We demonstrate that in a Warm Inflation (WI) setting the persistent thermal bath, sustained by the dissipative interactions with the inflaton field, can source a sizable DM abundance via the non-renormaliza… ▽ More We present a novel perspective on the role of inflation in the production of Dark Matter (DM). Specifically, we explore the DM production during Warm Inflation via ultraviolet Freeze-In (WIFI). We demonstrate that in a Warm Inflation (WI) setting the persistent thermal bath, sustained by the dissipative interactions with the inflaton field, can source a sizable DM abundance via the non-renormalizable interactions that connect the DM with the bath. Compared to the (conventional) radiation-dominated (RD) UV freeze-in scenario for the same reheat temperature (after inflation), the resulting DM yield in WIFI is always enhanced, showing a strongly positive dependence on the mass dimension of the non-renormalizable operator. Of particular interest, for a sufficiently large mass dimension of the operator, the entirety of the DM abundance of the Universe can be created during the inflationary phase. For the specific models we study, we find an enhancement in DM yield of up to 30 orders of magnitude relative to RD UV freeze-in for the same reheat temperature. Our findings also suggest a broader applicability for producing other cosmological relics, which may have a substantial impact on the evolution of the early Universe. △ Less

Submitted 30 January, 2024; originally announced January 2024.

Comments: 7 pages, 4 figures

Report number: UTWI-01-2024, NORDITA 2024-005

Origin of the Stochastic Gravitational Wave Background: First-Order Phase Transition vs. Black Hole Mergers

Authors: Martin Wolfgang Winkler, Katherine Freese

Abstract: The NANOGrav, Parkes and European Pulsar Timing Array (PTA) experiments have collected strong evidence for a stochastic gravitational wave background in the nHz-frequency band. In this work we perform a detailed statistical analysis of the signal in order to elucidate its physical origin. Specifically, we test the standard explanation in terms of supermassive black hole mergers against the promine… ▽ More The NANOGrav, Parkes and European Pulsar Timing Array (PTA) experiments have collected strong evidence for a stochastic gravitational wave background in the nHz-frequency band. In this work we perform a detailed statistical analysis of the signal in order to elucidate its physical origin. Specifically, we test the standard explanation in terms of supermassive black hole mergers against the prominent alternative explanation in terms of a first-order phase transition. By means of a frequentist hypothesis test we find that the observed gravitational wave spectrum prefers a first-order phase transition at $2-3σ$ significance compared to black hole mergers (depending on the underlying black hole model). This mild preference is linked to the relatively large amplitude of the observed gravitational wave signal (above the typical expectation of black hole models) and to its spectral shape (which slightly favors the phase-transition spectrum over the predominantly single power-law spectrum predicted in black hole models). The best fit to the combined PTA data set is obtained for a phase transition which dominantly produces the gravitational wave signal by bubble collisions (rather than by sound waves). The best-fit (energy-density) spectrum features, within the frequency band of the PTA experiments, a crossover from a steeply rising power law (causality tail) to a softly rising power law; the peak frequency then falls slightly above the PTA-measured range. Such a spectrum can be obtained for a strong first-order phase transition in the thick-wall regime of vacuum tunneling which reheats the Universe to a temperature of $T_*\sim \text{GeV}$. A dark sector phase transition at the GeV-scale provides a comparably good fit. △ Less

Submitted 24 January, 2024; originally announced January 2024.

Comments: 26 pages, 8 figures

Report number: NORDITA-2024-002, UT-WI-02-2024

UHZ1 and the other three most distant quasars observed: possible evidence for Supermassive Dark Stars

Authors: Cosmin Ilie, Katherine Freese, Andreea Petric, Jillian Paulin

Abstract: The James Webb Space Telescope (JWST) has recently uncovered a new record-breaking quasar, UHZ1, at a redshift of $z\sim10$. This discovery continues JWST's trend of confronting the expectations from the standard $Λ$CDM model of cosmology with challenges. Namely, too many very massive galaxies and quasars have been observed at very high redshifts, when the universe was only a few hundred million y… ▽ More The James Webb Space Telescope (JWST) has recently uncovered a new record-breaking quasar, UHZ1, at a redshift of $z\sim10$. This discovery continues JWST's trend of confronting the expectations from the standard $Λ$CDM model of cosmology with challenges. Namely, too many very massive galaxies and quasars have been observed at very high redshifts, when the universe was only a few hundred million years old. We have previously shown that Supermassive Dark Stars (SMDSs) may offer a solution to this puzzle. These fascinating objects would be the first stars in the universe, growing to be $\sim 10^5-10^7 M_{\odot}$ and shining as bright as $10^9$ suns. Unlike Population III stars (the major alternative proposed model for the first stars in the universe, which would also have zero metallicity and would be powered by nuclear fusion), SMDSs would be powered by dark matter heating (e.g. dark matter annihilation) and would be comparatively cooler. At the ends of their lives (when they run out of dark matter fuel), SMDSs would directly collapse into black holes, thus providing possible seeds for the first quasars. Previous papers have shown that to form at $z\sim10$, UHZ1 would require an incredibly massive seed ($\sim 10^4 -10^5 M_{\odot}$), which was assumed to be a Direct Collapse Black Hole (DCBH). In this paper, we demonstrate that Supermassive Dark Stars (SMDSs) offer an equally valid solution to the mystery of the first quasars, by examining the four most distant known quasars: UHZ1, J0313-1806, J1342+0928, and J1007+2115, with particular emphasis on UHZ1. △ Less

Submitted 21 December, 2023; originally announced December 2023.

Defying eternal inflation in warm inflation with a negative running

Authors: Gabriele Montefalcone, Rudnei O. Ramos, Gustavo S. Vicente, Katherine Freese

Abstract: It was pointed out previously~\cite{Kinney:2014jya} that a sufficiently negative running of the spectral index of curvature perturbations from (ordinary i.e. cold) inflation is able to prevent eternal inflation from ever occurring. Here, we reevaluate those original results, but in the context of warm inflation, in which a substantial radiation component (produced by the inflaton) exists throughou… ▽ More It was pointed out previously~\cite{Kinney:2014jya} that a sufficiently negative running of the spectral index of curvature perturbations from (ordinary i.e. cold) inflation is able to prevent eternal inflation from ever occurring. Here, we reevaluate those original results, but in the context of warm inflation, in which a substantial radiation component (produced by the inflaton) exists throughout the inflationary period. We demonstrate that the same general requirements found in the context of ordinary (cold) inflation also hold true in warm inflation; indeed an even tinier amount of negative running is sufficient to prevent eternal inflation. This is particularly pertinent, as models featuring negative running are more generic in warm inflation scenarios. Finally, the condition for the existence of eternal inflation in cold inflation -- that the curvature perturbation amplitude exceed unity on superhorizon scales -- becomes more restrictive in the case of warm inflation. The curvature perturbations must be even larger, i.e. even farther out on the potential, away from the part of the potential where observables, e.g. in the Cosmic Microwave Background, are produced. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Comments: 14 pages, Preprint Numbers: UTWI-40-2023, NORDITA 2023-067

The Gravitational Wave Spectrum of Chain Inflation

Authors: Katherine Freese, Aliki Litsa, Martin Wolfgang Winkler

Abstract: Chain inflation is an alternative to slow-roll inflation in which the inflaton tunnels along a large number of consecutive minima in its potential. In this work we perform the first comprehensive calculation of the gravitational wave spectrum of chain inflation. In contrast to slow-roll inflation the latter does not stem from quantum fluctuations of the gravitational field during inflation, but ra… ▽ More Chain inflation is an alternative to slow-roll inflation in which the inflaton tunnels along a large number of consecutive minima in its potential. In this work we perform the first comprehensive calculation of the gravitational wave spectrum of chain inflation. In contrast to slow-roll inflation the latter does not stem from quantum fluctuations of the gravitational field during inflation, but rather from the bubble collisions during the first-order phase transitions associated with vacuum tunneling. Our calculation is performed within an effective theory of chain inflation which builds on an expansion of the tunneling rate capturing most of the available model space. The effective theory can be seen as chain inflation's analogue of the slow-roll expansion in rolling models of inflation. We show that chain inflation produces a very characteristic double-peak spectrum: a faint high-frequency peak associated with the gravitational radiation emitted during inflation, and a strong low-frequency peak associated with the graceful exit from chain inflation (marking the transition to the radiation-dominated epoch). There exist very exciting prospects to test the gravitational wave signal from chain inflation at the aLIGO-aVIRGO-KAGRA network, at LISA and /or at pulsar timing array experiments. A particularly intriguing possibility we point out is that chain inflation could be the source of the stochastic gravitational wave background recently detected by NANOGrav, PPTA, EPTA and CPTA. We also show that the gravitational wave signal of chain inflation is often accompanied by running/ higher running of the scalar spectral index to be tested at future Cosmic Microwave Background experiments. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Comments: 33 pages, 7 figures

Report number: UTWI-30-2023, NORDITA-2023-045

WarmSPy: a numerical study of cosmological perturbations in warm inflation

Authors: Gabriele Montefalcone, Vikas Aragam, Luca Visinelli, Katherine Freese

Abstract: We present WarmSPy, a numerical code in Python designed to solve for the perturbations' equations in warm inflation models and compute the corresponding scalar power spectrum at CMB horizon crossing. In models of warm inflation, a radiation bath of temperature $T$ during inflation induces a dissipation (friction) rate of strength $Q \propto T^c/φ^m$ in the equation of motion for the inflaton field… ▽ More We present WarmSPy, a numerical code in Python designed to solve for the perturbations' equations in warm inflation models and compute the corresponding scalar power spectrum at CMB horizon crossing. In models of warm inflation, a radiation bath of temperature $T$ during inflation induces a dissipation (friction) rate of strength $Q \propto T^c/φ^m$ in the equation of motion for the inflaton field $φ$. While for a temperature-independent dissipation rate ($c=0$) an analytic expression for the scalar power spectrum exists, in the case of a non-zero value for $c$ the set of equations can only be solved numerically. For $c>0$ ($c<0$), the coupling between the perturbations in the inflaton field and radiation induces a growing (decaying) mode in the scalar perturbations, generally parameterized by a multiplicative function $G(Q)$ which we refer to as the scalar dissipation function. Using WarmSPy, we provide an analytic fit for $G(Q)$ for the cases of $c=\{3,1,-1\}$, corresponding to three cases that have been realized in physical models. Compared to previous literature results, our fits are more robust and valid over a broader range of dissipation strengths $Q\in[10^{-7},10^{4}]$. Additionally, for the first time, we numerically assess the stability of the scalar dissipation function against various model parameters, inflationary histories as well as the effects of metric perturbations. As a whole, the results do not depend appreciably on most of the parameters in the analysis, except for the dissipation index $c$, providing evidence for the universal behaviour of the scalar dissipation function $G(Q)$. △ Less

Submitted 23 July, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 37 pages, 8 figures, 1 table, Replacement due to typo on table 1 and figure 2

Report number: UTWI-25-2023, NORDITA-2023-031

Journal ref: JCAP 2401, 032 (2024)

Comparing representations of high-dimensional data with persistent homology: a case study in neuroimaging

Authors: Ty Easley, Kevin Freese, Elizabeth Munch, Janine Bijsterbosch

Abstract: Despite much attention, the comparison of reduced-dimension representations of high-dimensional data remains a challenging problem in multiple fields, especially when representations remain high-dimensional compared to sample size. We offer a framework for evaluating the topological similarity of high-dimensional representations of very high-dimensional data, a regime where topological structure i… ▽ More Despite much attention, the comparison of reduced-dimension representations of high-dimensional data remains a challenging problem in multiple fields, especially when representations remain high-dimensional compared to sample size. We offer a framework for evaluating the topological similarity of high-dimensional representations of very high-dimensional data, a regime where topological structure is more likely captured in the distribution of topological "noise" than a few prominent generators. Treating each representational map as a metric embedding, we compute the Vietoris-Rips persistence of its image. We then use the topological bootstrap to analyze the re-sampling stability of each representation, assigning a "prevalence score" for each nontrivial basis element of its persistence module. Finally, we compare the persistent homology of representations using a prevalence-weighted variant of the Wasserstein distance. Notably, our method is able to compare representations derived from different samples of the same distribution and, in particular, is not restricted to comparisons of graphs on the same vertex set. In addition, representations need not lie in the same metric space. We apply this analysis to a cross-sectional sample of representations of functional neuroimaging data in a large cohort and hierarchically cluster under the prevalence-weighted Wasserstein. We find that the ambient dimension of a representation is a stronger predictor of the number and stability of topological features than its decomposition rank. Our findings suggest that important topological information lies in repeatable, low-persistence homology generators, whose distributions capture important and interpretable differences between high-dimensional data representations. △ Less

Submitted 23 November, 2023; v1 submitted 23 June, 2023; originally announced June 2023.

Comments: Revision of preprint of submission for NeurIPS 2023 conference

Detectability of Supermassive Dark Stars with the Roman Space Telescope

Authors: Saiyang Zhang, Cosmin Ilie, Katherine Freese

Abstract: The first bright objects to form in the Universe at redshift $z \sim 10-20$ might have been Dark Stars, made primarily of hydrogen and helium but powered by dark matter. In this study, we investigate the detectability of Supermassive Dark Stars (SMDS) by the Roman Space Telescope. RST will be able to detect SMDSs at redshifts as high as $z\simeq 14$. In cases with gravitational lensing factors of… ▽ More The first bright objects to form in the Universe at redshift $z \sim 10-20$ might have been Dark Stars, made primarily of hydrogen and helium but powered by dark matter. In this study, we investigate the detectability of Supermassive Dark Stars (SMDS) by the Roman Space Telescope. RST will be able to detect SMDSs at redshifts as high as $z\simeq 14$. In cases with gravitational lensing factors of $μ\sim 100$, RST will be able to find SMDS as small as $\sim10^4 M_{\odot}$ at $z\sim 12$ with $\sim 10^6$ s of exposure. To differentiate SMDS from early galaxies containing zero metallicity stars at similar redshifts, we compare their spectra, photometry in RST bands, color indexes and image morphology. With RST alone, the differentiation is possible only for limited cases: SMDS formed via "adiabatic contraction" (DM pulled into the star via gravity alone) with $M\gtrsim 10^5M_{\odot}$ and lensed by $μ\gtrsim 30$ have distinct photometric signatures from those of the first galaxies. For SMDSs formed via "dark matter capture," their spectra are degenerate to those of many galaxies with little to no nebular emission. Thus with RST alone, the only way to tell them apart from first galaxies would be via image morphology: i.e. point object (SMDSs) vs. extended object (sufficiently magnified galaxies). However, if the same objects are further examined by JWST spectroscopy, a "smoking gun" for detection of SMDS is the HeII $λ$1640 absorption line. While RST does not cover the wavelength band required to find this line (for $z_{\rm emi}\gtrsim 10$), JWST does. Hence the two detectors can be used together in identifying SMDS. The confirmed detection of any SMDSs will provide evidence for a new type of star, powered by dark matter. Moreover, such massive stars can also be natural progenitors of the supermassive black holes powering the extremely bright quasars observed at $z\gtrsim 6$. △ Less

Submitted 12 April, 2024; v1 submitted 20 June, 2023; originally announced June 2023.

Comments: 49 pages, 24 figures

Report number: NORDITA-2023-027, UTWI-24-2023

Journal ref: ApJ 965 (2024) 121

Higgs Inflation and the Electroweak Gauge Sector

Authors: Stephon Alexander, Cyril Creque-Sarbinowski, Humberto Gilmer, Katherine Freese

Abstract: We introduce a new method that allows for the Higgs to be the inflaton. That is, we let the Higgs be a pseudo-Nambu-Goldstone (pNG) boson of a global coset symmetry $G/H$ that spontaneously breaks at an energy scale $\sim 4πf$ and give it a suitable $SU(2) \subset G$ Chern-Simons interaction, with $β$ the dimensionless Chern-Simons coupling strength and $f$ an $SU(2)$ decay constant. As a result,… ▽ More We introduce a new method that allows for the Higgs to be the inflaton. That is, we let the Higgs be a pseudo-Nambu-Goldstone (pNG) boson of a global coset symmetry $G/H$ that spontaneously breaks at an energy scale $\sim 4πf$ and give it a suitable $SU(2) \subset G$ Chern-Simons interaction, with $β$ the dimensionless Chern-Simons coupling strength and $f$ an $SU(2)$ decay constant. As a result, slow-roll inflation occurs via $SU(2)$-induced friction down a steep sinusoidal potential. In order to obey electroweak $SU(2)_{\rm L}\times U(1)_Y$ symmetry, the lowest-order Chern-Simons interaction is required to be quadratic in the Higgs with coupling strength $\propto β^2/f^2$. Higher-order interaction terms keep the full Lagrangian nearly invariant under the approximate pNG shift symmetry. Employing the simplest symmetry coset $SU(5)/SO(5)$, $N$ $e$-folds of inflation occur when $N \approx 60 \left(g/0.64\right)^2\left[β/\left(3\times 10^6\right)\right]^{8/3}\left[f/\left(5\times 10^{11}\ {\rm GeV}\right)\right]^{2/3}$, with $g$ the weak isospin gauge coupling constant. Small values of the decay constant, $f \lesssim 5 \times 10^{11} {\rm GeV}$, which are needed to address the Higgs hierarchy problem, are ruled out by electric dipole measurements and so successfully explaining inflation requires large $β$. We discuss possible methods to achieve such large couplings and other alternative Higgs inflation scenarios outside the standard modified-gravity framework. △ Less

Submitted 26 July, 2023; v1 submitted 7 June, 2023; originally announced June 2023.

Comments: 11 pages, including appendices, 4 figures

Supermassive Dark Star candidates seen by JWST?

Authors: Cosmin Ilie, Jillian Paulin, Katherine Freese

Abstract: The first generation of stars in the Universe is yet to be observed. There are two leading theories for those objects that mark the beginning of the cosmic dawn: hydrogen burning Population~III stars and Dark Stars, made of hydrogen and helium but powered by Dark Matter heating. The latter can grow to become supermassive ($M_\star\sim 10^6\Msun$) and extremely bright ($L\sim 10^9L_\odot$). We show… ▽ More The first generation of stars in the Universe is yet to be observed. There are two leading theories for those objects that mark the beginning of the cosmic dawn: hydrogen burning Population~III stars and Dark Stars, made of hydrogen and helium but powered by Dark Matter heating. The latter can grow to become supermassive ($M_\star\sim 10^6\Msun$) and extremely bright ($L\sim 10^9L_\odot$). We show that each of the following three objects: JADES-GS-z13-0, JADES-GS-z12-0, and JADES-GS-z11-0 (at redshifts $z\in[11,14]$) are consistent with a Supermassive Dark Star interpretation, thus identifying, for the first time, Dark Star candidates. △ Less

Submitted 12 April, 2023; v1 submitted 3 April, 2023; originally announced April 2023.

Comments: Added complete affiliations for KF. Added preprint numbers from UT Austin and Nordita. Edited acknowledgments to include former collaborators on Dark Stars

Report number: NORDITA-2023-025,UTWI-13-2023

Dark Matter and Gravity Waves from a Dark Big Bang

Authors: Katherine Freese, Martin Wolfgang Winkler

Abstract: The Hot Big Bang is often considered as the origin of all matter and radiation in the Universe. Primordial nucleosynthesis (BBN) provides strong evidence that the early Universe contained a hot plasma of photons and baryons with a temperature $T>\text{MeV}$. However, the earliest probes of dark matter originate from much later times around the epoch of structure formation. In this work we describe… ▽ More The Hot Big Bang is often considered as the origin of all matter and radiation in the Universe. Primordial nucleosynthesis (BBN) provides strong evidence that the early Universe contained a hot plasma of photons and baryons with a temperature $T>\text{MeV}$. However, the earliest probes of dark matter originate from much later times around the epoch of structure formation. In this work we describe a scenario in which dark matter (and possibly dark radiation) can be formed around or even after BBN in a second Big Bang which we dub the ``Dark Big Bang''. The latter occurs through a phase transition in the dark sector which transforms dark vacuum energy into a hot dark plasma of particles; in this paper we focus on a first-order phase transition for the Dark Big Bang. The correct dark matter abundance can be set by dark matter cannibalism or by pair-annihilation within the dark sector followed by a thermal freeze-out. Alternatively ultra-heavy ``dark-zilla'' dark matter can originate directly from bubble collisions during the Dark Big Bang. We will show that the Dark Big Bang is consistent with constraints from structure formation and the Cosmic Microwave Background (CMB) if it occurred when the Universe was less than one month old, corresponding to a temperature in the visible sector above $\mathcal{O}$(keV). While the dark matter evades direct and indirect detection, the Dark Big Bang gives rise to striking gravity wave signatures to be tested at pulsar timing array experiments. Furthermore, the Dark Big Bang allows for realizations of self-interacting and/or warm dark matter which suggest exciting discovery potential in future small-scale structure observations. △ Less

Submitted 22 February, 2023; originally announced February 2023.

Comments: 46 pages, 9 figures

Report number: UTWI-06-2023, NORDITA-2023-003

Mineral Detection of Neutrinos and Dark Matter. A Whitepaper

Authors: Sebastian Baum, Patrick Stengel, Natsue Abe, Javier F. Acevedo, Gabriela R. Araujo, Yoshihiro Asahara, Frank Avignone, Levente Balogh, Laura Baudis, Yilda Boukhtouchen, Joseph Bramante, Pieter Alexander Breur, Lorenzo Caccianiga, Francesco Capozzi, Juan I. Collar, Reza Ebadi, Thomas Edwards, Klaus Eitel, Alexey Elykov, Rodney C. Ewing, Katherine Freese, Audrey Fung, Claudio Galelli, Ulrich A. Glasmacher, Arianna Gleason , et al. (44 additional authors not shown)

Abstract: Minerals are solid state nuclear track detectors - nuclear recoils in a mineral leave latent damage to the crystal structure. Depending on the mineral and its temperature, the damage features are retained in the material from minutes (in low-melting point materials such as salts at a few hundred degrees C) to timescales much larger than the 4.5 Gyr-age of the Solar System (in refractory materials… ▽ More Minerals are solid state nuclear track detectors - nuclear recoils in a mineral leave latent damage to the crystal structure. Depending on the mineral and its temperature, the damage features are retained in the material from minutes (in low-melting point materials such as salts at a few hundred degrees C) to timescales much larger than the 4.5 Gyr-age of the Solar System (in refractory materials at room temperature). The damage features from the $O(50)$ MeV fission fragments left by spontaneous fission of $^{238}$U and other heavy unstable isotopes have long been used for fission track dating of geological samples. Laboratory studies have demonstrated the readout of defects caused by nuclear recoils with energies as small as $O(1)$ keV. This whitepaper discusses a wide range of possible applications of minerals as detectors for $E_R \gtrsim O(1)$ keV nuclear recoils: Using natural minerals, one could use the damage features accumulated over $O(10)$ Myr$-O(1)$ Gyr to measure astrophysical neutrino fluxes (from the Sun, supernovae, or cosmic rays interacting with the atmosphere) as well as search for Dark Matter. Using signals accumulated over months to few-years timescales in laboratory-manufactured minerals, one could measure reactor neutrinos or use them as Dark Matter detectors, potentially with directional sensitivity. Research groups in Europe, Asia, and America have started developing microscopy techniques to read out the $O(1) - O(100)$ nm damage features in crystals left by $O(0.1) - O(100)$ keV nuclear recoils. We report on the status and plans of these programs. The research program towards the realization of such detectors is highly interdisciplinary, combining geoscience, material science, applied and fundamental physics with techniques from quantum information and Artificial Intelligence. △ Less

Submitted 16 May, 2023; v1 submitted 17 January, 2023; originally announced January 2023.

Comments: 115 pages, many pictures of tracks. Please see the source file for higher resolution versions of some plots. v2: matches the published version

Journal ref: Phys.Dark Univ. 41 (2023) 101245

A Gap No More: Mechanism for Non-Nuclear Energy to Fill in the Black Hole Mass Gap

Authors: Joshua Ziegler, Katherine Freese

Abstract: Standard stellar evolution models predict that black holes in the range of approximately $50 - 140 M_\odot$ should not exist directly from stellar evolution. This gap appears because stars with masses between 100 and 240 $M_\odot$ are expected to undergo a pair instability supernova and leave behind no remnant, or a pulsational pair instability supernova and leave behind a remnant much smaller tha… ▽ More Standard stellar evolution models predict that black holes in the range of approximately $50 - 140 M_\odot$ should not exist directly from stellar evolution. This gap appears because stars with masses between 100 and 240 $M_\odot$ are expected to undergo a pair instability supernova and leave behind no remnant, or a pulsational pair instability supernova and leave behind a remnant much smaller than their initial stellar mass. However, black holes have been discovered by the LIGO/Virgo collaboration within this mass range. In previous work (arXiv:2010.00254), we used the stellar evolution code MESA to show that the addition of non-nuclear energy (such as from annihilation of dark matter) could alter the evolution of a 180 $M_\odot$ star so that the observed black holes could be produced from isolated stars. In this paper, we extend this analysis to stars of other masses, and find that sufficient amounts of non-nuclear energy can allow any star to avoid pair instability, and could produce a black hole of mass comparable to the initial stellar mass. In addition, we produce examples of the type of black hole initial mass function that can be produced from this mechanism. These illustrative examples suggest that adding non-nuclear energy to stars offers a way to fully close the mass gap. △ Less

Submitted 11 January, 2023; v1 submitted 21 December, 2022; originally announced December 2022.

Comments: 34 pages, 7 figures

Report number: UTWI-25-2022, NORDITA-2022-101

Observational Constraints on Warm Natural Inflation

Authors: Gabriele Montefalcone, Vikas Aragam, Luca Visinelli, Katherine Freese

Abstract: Warm natural inflation is studied for the case of the original cosine potential. The radiation bath during inflation induces a dissipation (friction) rate in the equation of motion for the inflaton field, which can potentially reduce the field excursion needed for an observationally viable period of inflation. We examine if the dissipation thus provides a mechanism to avoid the large decay constan… ▽ More Warm natural inflation is studied for the case of the original cosine potential. The radiation bath during inflation induces a dissipation (friction) rate in the equation of motion for the inflaton field, which can potentially reduce the field excursion needed for an observationally viable period of inflation. We examine if the dissipation thus provides a mechanism to avoid the large decay constant $f \gtrsim M_{\mathrm{pl}}$ of cold cosine natural inflation. Whereas temperature independent dissipation has previously been shown to alleviate the need for a trans-Planckian decay constant $f$, we illustrate here the difficulties of accommodating a significantly sub-Planckian decay constant ($f<10^{-1}M_{\mathrm{pl}}$) in the case of the following temperature dependent dissipation rates, $Γ\propto T^c$, with $c=\{1,3\}$. Such dissipation rates represent physically well-motivated constructions in the literature. For each model, we map its location in the $r$-$n_s$ plane and compare with Cosmic Microwave Background data. For $c=1 \, (c=3)$, we find that agreement with CMB data requires that dissipation be in the weak (moderate) regime and that the minimum allowed value of the decay constant in the potential is $f_{\rm min} = 0.3 \, (0.8)\,M_{\mathrm{pl}}$ respectively. △ Less

Submitted 17 January, 2023; v1 submitted 8 December, 2022; originally announced December 2022.

Comments: 23 pages, 3 figures, 1 table, Replacement due to typo on DOE acknowledgement number

Journal ref: JCAP 2303, 002 (2023)

Dimensionally Reduced Waveforms for Spin-Induced Quadrupole Searches

Authors: Horng Sheng Chia, Thomas D. P. Edwards, Richard N. George, Aaron Zimmerman, Adam Coogan, Katherine Freese, Cody Messick, Christian N. Setzer

Abstract: We present highly accurate, dimensionally-reduced gravitational waveforms for binary inspirals whose components have large spin-induced quadrupole moments. The spin-induced quadrupole of a body first appears in the phase of a waveform at the early inspiral stage of the binary coalescence, making it a relatively clean probe of the internal structure of the body. However, for objects with large quad… ▽ More We present highly accurate, dimensionally-reduced gravitational waveforms for binary inspirals whose components have large spin-induced quadrupole moments. The spin-induced quadrupole of a body first appears in the phase of a waveform at the early inspiral stage of the binary coalescence, making it a relatively clean probe of the internal structure of the body. However, for objects with large quadrupolar deviations from Kerr, searches using binary black hole (BBH) models would be ineffective. In order to perform a computationally-feasible search, we present two dimensionally-reduced models which are derived from the original six-dimensional post-Newtonian waveform for such systems. Our dimensional reduction method is guided by power counting in the post-Newtonian expansion, suitable reparameterizations of the source physics, and truncating terms in the phase that are small in most physically well-motivated regions of parameter space. In addition, we note that large quadrupolar deviations cause the frequency at which a binary system reaches its minimum binding energy to be reduced substantially. This minimum signals the end of the inspiral regime and provides a natural cutoff for the PN waveform. We provide accurate analytic estimates for these frequency cutoffs. Finally, we perform injection studies to test the effectualness of the dimensionally reduced waveforms. We find that over $80\%$ of the injections have an effectualness of $\varepsilon > 0.999$, significantly higher than is typically required for standard BBH banks, for systems with component spins of $|χ_i| \lesssim 0.6$ and dimensionless quadrupole of $κ_i \lesssim 10^3$. Importantly, these waveforms represent an essential first step towards enabling an effective search for astrophysical objects with large quadrupoles. △ Less

Submitted 31 October, 2022; originally announced November 2022.

Comments: 18 pages, 10 figures

Report number: UTWI-13-2022, LIGO-P2200325

Constraints on the scalar-field potential in warm inflation

Authors: Gabriele Montefalcone, Vikas Aragam, Luca Visinelli, Katherine Freese

Abstract: We quantify the degree of fine tuning required to achieve an observationally viable period of inflation in the strongly dissipative regime of warm inflation. The ``fine-tuning'' parameter $λ$ is taken to be the ratio of the change in the height of the potential $ΔV$ to the change in the scalar field $(Δφ)^{4}$, i.e. the width of the potential, and therefore measures the requisite degree of flatnes… ▽ More We quantify the degree of fine tuning required to achieve an observationally viable period of inflation in the strongly dissipative regime of warm inflation. The ``fine-tuning'' parameter $λ$ is taken to be the ratio of the change in the height of the potential $ΔV$ to the change in the scalar field $(Δφ)^{4}$, i.e. the width of the potential, and therefore measures the requisite degree of flatness in the potential. The best motivated warm inflationary scenarios involve a dissipation rate of the kind $Γ\propto T^c$ with $c\geq 0$, and for all such cases, the bounds on $λ$ are tighter than those for standard cold inflation by at least 3 orders of magnitude. In other words, these models require an even flatter potential than standard inflation. On the other hand for the case of warm inflation with $c< 0$, we find that in a strongly dissipative regime the bound on $λ$ can significantly weaken with respect to cold inflation. Thus, if a warm inflation model can be constructed in a strongly dissipative, negatively temperature-dependent regime, it accommodates steeper potentials otherwise ruled out in standard inflation. △ Less

Submitted 17 January, 2023; v1 submitted 29 September, 2022; originally announced September 2022.

Comments: 6 pages, Preprint Numbers: UTWI-03-2022, NORDITA 2022-068, Replacement due to typo on DOE acknowledgement number

Journal ref: Phys. Rev. D 107, 063543 (2023)

Have Pulsar Timing Arrays detected the Hot Big Bang? Gravitational Waves from Strong First Order Phase Transitions in the Early Universe

Authors: Katherine Freese, Martin Wolfgang Winkler

Abstract: The origins of matter and radiation in the universe lie in a Hot Big Bang. We present a number of well-motivated cosmologies in which the Big Bang occurs through a strong first order phase transition -- either at the end of inflation, after a period of kination ("Kination-Induced Big Bang"), or after a second period of vacuum-domination in the early universe ("Supercooled Big Bang"); we also propo… ▽ More The origins of matter and radiation in the universe lie in a Hot Big Bang. We present a number of well-motivated cosmologies in which the Big Bang occurs through a strong first order phase transition -- either at the end of inflation, after a period of kination ("Kination-Induced Big Bang"), or after a second period of vacuum-domination in the early universe ("Supercooled Big Bang"); we also propose a "Dark Big Bang" where only the dark matter in the Universe is created in a first-order phase transition much after inflation. In all of these scenarios, the resulting gravitational radiation can explain the tentative signals reported by the NANOGrav, Parkes and European Pulsar Timing Array experiments if the reheating temperature of the Hot Big Bang, and correspondingly the energy scale of the false vacuum, falls in the range $T_* \sim ρ_{\rm vac}^{1/4} $= MeV--100 GeV. All the same models at higher reheating temperatures will be of interest to upcoming ground- and space-based interferometer searches for gravitational waves at larger frequency. △ Less

Submitted 5 August, 2022; originally announced August 2022.

Comments: 41 pages, 11 figures

Report number: UTTG 11-2022, NORDITA 2022-056

Updated neutrino mass constraints from galaxy clustering and CMB lensing-galaxy cross-correlation measurements

Authors: Isabelle Tanseri, Steffen Hagstotz, Sunny Vagnozzi, Elena Giusarma, Katherine Freese

Abstract: We revisit cosmological constraints on the sum of the neutrino masses $Σm_ν$ from a combination of full-shape BOSS galaxy clustering [$P(k)$] data and measurements of the cross-correlation between Planck Cosmic Microwave Background (CMB) lensing convergence and BOSS galaxy overdensity maps [$C^{κ\text{g}}_{\ell}$], using a simple but theoretically motivated model for the scale-dependent galaxy bia… ▽ More We revisit cosmological constraints on the sum of the neutrino masses $Σm_ν$ from a combination of full-shape BOSS galaxy clustering [$P(k)$] data and measurements of the cross-correlation between Planck Cosmic Microwave Background (CMB) lensing convergence and BOSS galaxy overdensity maps [$C^{κ\text{g}}_{\ell}$], using a simple but theoretically motivated model for the scale-dependent galaxy bias in auto- and cross-correlation measurements. We improve upon earlier related work in several respects, particularly through a more accurate treatment of the correlation and covariance between $P(k)$ and $C^{κ\text{g}}_{\ell}$ measurements. When combining these measurements with Planck CMB data, we find a 95% confidence level upper limit of $Σm_ν<0.14\,{\rm eV}$, while slightly weaker limits are obtained when including small-scale ACTPol CMB data, in agreement with our expectations. We confirm earlier findings that (once combined with CMB data) the full-shape information content is comparable to the geometrical information content in the reconstructed BAO peaks given the precision of current galaxy clustering data, discuss the physical significance of our inferred bias and shot noise parameters, and perform a number of robustness tests on our underlying model. While the inclusion of $C^{κ\text{g}}_{\ell}$ measurements does not currently appear to lead to substantial improvements in the resulting $Σm_ν$ constraints, we expect the converse to be true for near-future galaxy clustering measurements, whose shape information content will eventually supersede the geometrical one. △ Less

Submitted 22 July, 2022; v1 submitted 5 July, 2022; originally announced July 2022.

Comments: 52 pages, 7 sub-figures arranged in 6 figures, 2 tables. v2: additional references added, version accepted for publication in JHEAp

Journal ref: JHEAp 36 (2022) 1-26

Dark Stars Powered by Self-Interacting Dark Matter

Authors: Youjia Wu, Sebastian Baum, Katherine Freese, Luca Visinelli, Hai-Bo Yu

Abstract: Dark matter annihilation might power the first luminous stars in the Universe. These types of stars, known as dark stars, could form in $(10^6\mathrm{-}10^8)\,M_\odot$ protohalos at redshifts $z \sim 20$, and they could be much more luminous and larger in size than ordinary stars powered by nuclear fusion. We investigate the formation of dark stars in the self-interacting dark matter (SIDM) scenar… ▽ More Dark matter annihilation might power the first luminous stars in the Universe. These types of stars, known as dark stars, could form in $(10^6\mathrm{-}10^8)\,M_\odot$ protohalos at redshifts $z \sim 20$, and they could be much more luminous and larger in size than ordinary stars powered by nuclear fusion. We investigate the formation of dark stars in the self-interacting dark matter (SIDM) scenario. We present a concrete particle physics model of SIDM that can simultaneously give rise to the observed dark matter density, satisfy constraints from astrophysical and terrestrial searches, and address the various small-scale problems of collisionless dark matter via the self-interactions. In this model, the power from dark matter annihilation is deposited in the baryonic gas in environments where dark stars could form. We further study the evolution of SIDM density profiles in the protohalos at $z \sim 20$. As the baryon cloud collapses due to the various cooling processes, the deepening gravitational potential can speed up gravothermal evolution of the SIDM halo, yielding sufficiently high dark matter densities for dark stars to form. We find that SIDM-powered dark stars can have similar properties, such as their luminosity and size, as dark stars predicted in collisionless dark matter models. △ Less

Submitted 29 August, 2022; v1 submitted 22 May, 2022; originally announced May 2022.

Comments: 11 pages, 4 figures. Matches the published version

Report number: UTTG 12-2021; LCTP-21-18; NORDITA 2021-068

Journal ref: Phys. Rev. D 106 (4) 043028 (2022)

Non-Universal Stellar Initial Mass Functions: Large Uncertainties in Star Formation Rates at $z\approx 2-4$ and Other Astrophysical Probes

Authors: Joshua J. Ziegler, Thomas D. P. Edwards, Anna M. Suliga, Irene Tamborra, Shunsaku Horiuchi, Shin'ichiro Ando, Katherine Freese

Abstract: We explore the assumption, widely used in many astrophysical calculations, that the stellar initial mass function (IMF) is universal across all galaxies. By considering both a canonical Salpeter-like IMF and a non-universal IMF, we are able to compare the effect of different IMFs on multiple observables and derived quantities in astrophysics. Specifically, we consider a non-universal IMF which var… ▽ More We explore the assumption, widely used in many astrophysical calculations, that the stellar initial mass function (IMF) is universal across all galaxies. By considering both a canonical Salpeter-like IMF and a non-universal IMF, we are able to compare the effect of different IMFs on multiple observables and derived quantities in astrophysics. Specifically, we consider a non-universal IMF which varies as a function of the local star formation rate, and explore the effects on the star formation rate density (SFRD), the extragalactic background light, the supernova (both core-collapse and thermonuclear) rates, and the diffuse supernova neutrino background. Our most interesting result is that our adopted varying IMF leads to much greater uncertainty on the SFRD at $z \approx 2-4$ than is usually assumed. Indeed, we find a SFRD (inferred using observed galaxy luminosity distributions) that is a factor of $\gtrsim 3$ lower than canonical results obtained using a universal Salpeter-like IMF. Secondly, the non-universal IMF we explore implies a reduction in the supernova core-collapse rate of a factor of $\sim2$, compared against a universal IMF. The other potential tracers are only slightly affected by changes to the properties of the IMF. We find that currently available data do not provide a clear preference for universal or non-universal IMF. However, improvements to measurements of the star formation rate and core-collapse supernova rate at redshifts $z \gtrsim 2$ may offer the best prospects for discernment. △ Less

Submitted 16 May, 2022; originally announced May 2022.

Comments: 15 pages, 11 figures, 1 appendix

Report number: UTTG-05-2022, N3AS-22-007, NORDITA 2022-026

Journal ref: MNRAS 517, 2471 (2022)

Neutrino mass and mass ordering: No conclusive evidence for normal ordering

Authors: Stefano Gariazzo, Martina Gerbino, Thejs Brinckmann, Massimiliano Lattanzi, Olga Mena, Thomas Schwetz, Shouvik Roy Choudhury, Katherine Freese, Steen Hannestad, Christoph A. Ternes, Mariam Tórtola

Abstract: The extraction of the neutrino mass ordering is one of the major challenges in particle physics and cosmology, not only for its implications for a fundamental theory of mass generation in nature, but also for its decisive role in the scale of future neutrinoless double beta decay experimental searches. It has been recently claimed that current oscillation, beta decay and cosmological limits on the… ▽ More The extraction of the neutrino mass ordering is one of the major challenges in particle physics and cosmology, not only for its implications for a fundamental theory of mass generation in nature, but also for its decisive role in the scale of future neutrinoless double beta decay experimental searches. It has been recently claimed that current oscillation, beta decay and cosmological limits on the different observables describing the neutrino mass parameter space provide robust decisive Bayesian evidence in favor of the normal ordering of the neutrino mass spectrum [arXiv:2203.14247]. We further investigate these strong claims using a rich and wide phenomenology, with different sampling techniques of the neutrino parameter space. Contrary to the findings of Jimenez et al [arXiv:2203.14247], no decisive evidence for the normal mass ordering is found. Neutrino mass ordering analyses must rely on priors and parameterizations that are ordering-agnostic: robust results should be regarded as those in which the preference for the normal neutrino mass ordering is driven exclusively by the data, while we find a difference of up to a factor of 33 in the Bayes factors among the different priors and parameterizations exploited here. An ordering-agnostic prior would be represented by the case of parameterizations sampling over the two mass splittings and a mass scale, or those sampling over the individual neutrino masses via normal prior distributions only. In this regard, we show that the current significance in favor of the normal mass ordering should be taken as $2.7σ$ (i.e. moderate evidence), mostly driven by neutrino oscillation data. △ Less

Submitted 4 May, 2022; originally announced May 2022.

Comments: 10 pages, 4 figures

Journal ref: JCAP 10 (2022) 010

Snowmass2021 Cosmic Frontier: Cosmic Microwave Background Measurements White Paper

Authors: Clarence L. Chang, Kevin M. Huffenberger, Bradford A. Benson, Federico Bianchini, Jens Chluba, Jacques Delabrouille, Raphael Flauger, Shaul Hanany, William C. Jones, Alan J. Kogut, Jeffrey J. McMahon, Joel Meyers, Neelima Sehgal, Sara M. Simon, Caterina Umilta, Kevork N. Abazajian, Zeeshan Ahmed, Yashar Akrami, Adam J. Anderson, Behzad Ansarinejad, Jason Austermann, Carlo Baccigalupi, Denis Barkats, Darcy Barron, Peter S. Barry , et al. (107 additional authors not shown)

Abstract: This is a solicited whitepaper for the Snowmass 2021 community planning exercise. The paper focuses on measurements and science with the Cosmic Microwave Background (CMB). The CMB is foundational to our understanding of modern physics and continues to be a powerful tool driving our understanding of cosmology and particle physics. In this paper, we outline the broad and unique impact of CMB science… ▽ More This is a solicited whitepaper for the Snowmass 2021 community planning exercise. The paper focuses on measurements and science with the Cosmic Microwave Background (CMB). The CMB is foundational to our understanding of modern physics and continues to be a powerful tool driving our understanding of cosmology and particle physics. In this paper, we outline the broad and unique impact of CMB science for the High Energy Cosmic Frontier in the upcoming decade. We also describe the progression of ground-based CMB experiments, which shows that the community is prepared to develop the key capabilities and facilities needed to achieve these transformative CMB measurements. △ Less

Submitted 15 March, 2022; originally announced March 2022.

Comments: contribution to Snowmass 2021

Constraining Axion-Like Particles with HAWC Observations of TeV Blazars

Authors: Sunniva Jacobsen, Tim Linden, Katherine Freese

Abstract: Axion-like particles (ALPs) are a broad class of pseudo-scalar bosons that generically arise from broken symmetries in extensions of the standard model. In many scenarios, ALPs can mix with photons in regions with high magnetic fields. Photons from distant sources can mix with ALPs, which then travel unattenuated through the Universe, before they mix back to photons in the Milky Way galactic magne… ▽ More Axion-like particles (ALPs) are a broad class of pseudo-scalar bosons that generically arise from broken symmetries in extensions of the standard model. In many scenarios, ALPs can mix with photons in regions with high magnetic fields. Photons from distant sources can mix with ALPs, which then travel unattenuated through the Universe, before they mix back to photons in the Milky Way galactic magnetic field. Thus, photons can traverse regions where their signals would normally be blocked or attenuated. In this paper, we study TeV $γ$-ray observations from distant blazars, utilizing the significant $γ$-ray attenuation expected from such signals to look for excess photon fluxes that may be due to ALP-photon mixing. We find no such excesses among a stacked population of seven blazars and constrain the ALP-photon coupling constant to fall below $\sim$3$\times$10$^{-11}$ GeV$^{-1}$ for ALP masses below 300 neV. These results are competitive with, or better than, leading terrestrial and astrophysical constraints in this mass range. △ Less

Submitted 8 March, 2022; originally announced March 2022.

Efficient Gravitational Wave Template Bank Generation with Differentiable Waveforms

Authors: Adam Coogan, Thomas D. P. Edwards, Horng Sheng Chia, Richard N. George, Katherine Freese, Cody Messick, Christian N. Setzer, Christoph Weniger, Aaron Zimmerman

Abstract: The most sensitive search pipelines for gravitational waves from compact binary mergers use matched filters to extract signals from the noisy data stream coming from gravitational wave detectors. Matched-filter searches require banks of template waveforms covering the physical parameter space of the binary system. Unfortunately, template bank construction can be a time-consuming task. Here we pres… ▽ More The most sensitive search pipelines for gravitational waves from compact binary mergers use matched filters to extract signals from the noisy data stream coming from gravitational wave detectors. Matched-filter searches require banks of template waveforms covering the physical parameter space of the binary system. Unfortunately, template bank construction can be a time-consuming task. Here we present a new method for efficiently generating template banks that utilizes automatic differentiation to calculate the parameter space metric. Principally, we demonstrate that automatic differentiation enables accurate computation of the metric for waveforms currently used in search pipelines, whilst being computationally cheap. Additionally, by combining random template placement and a Monte Carlo method for evaluating the fraction of the parameter space that is currently covered, we show that search-ready template banks for frequency-domain waveforms can be rapidly generated. Finally, we argue that differentiable waveforms offer a pathway to accelerating stochastic placement algorithms. We implement all our methods into an easy-to-use Python package based on the jax framework, diffbank, to allow the community to easily take advantage of differentiable waveforms for future searches. △ Less

Submitted 30 November, 2022; v1 submitted 18 February, 2022; originally announced February 2022.

Comments: 15 pages, 6 figures. Comments welcome! Code can be found at https://github.com/adam-coogan/diffbank. Updated to match PRD version

Report number: UTTG 26-2021, NORDITA-2022-004

Neutrino point source searches for dark matter spikes

Authors: Katherine Freese, Irina Galstyan, Pearl Sandick, Patrick Stengel

Abstract: Any dark matter spikes surrounding black holes in our Galaxy are sites of significant dark matter annihilation, leading to a potentially detectable neutrino signal. In this paper we examine $10-10^5 M_\odot$ black holes associated with dark matter spikes that formed in early minihalos and still exist in our Milky Way Galaxy today, in light of neutrino data from the ANTARES and IceCube detectors. I… ▽ More Any dark matter spikes surrounding black holes in our Galaxy are sites of significant dark matter annihilation, leading to a potentially detectable neutrino signal. In this paper we examine $10-10^5 M_\odot$ black holes associated with dark matter spikes that formed in early minihalos and still exist in our Milky Way Galaxy today, in light of neutrino data from the ANTARES and IceCube detectors. In various regions of the sky, we determine the minimum distance away from the solar system that a dark matter spike must be in order to have not been detected as a neutrino point source for a variety of representative dark matter annihilation channels. Given these constraints on the distribution of dark matter spikes in the Galaxy, we place significant limits on the formation of the first generation of stars in early minihalos -- stronger than previous limits from gamma-ray searches in Fermi Gamma-Ray Space Telescope data. The larger black holes considered in this paper may arise as the remnants of Dark Stars after the dark matter fuel is exhausted; thus neutrino observations may be used to constrain the properties of Dark Stars. The limits are particularly strong for heavier WIMPs. For WIMP masses $\sim 5 \,$TeV, we show that $\lesssim 10 \%$ of minihalos can host first stars that collapse into BHs larger than $10^3 M_\odot$. △ Less

Submitted 8 September, 2022; v1 submitted 2 February, 2022; originally announced February 2022.

Comments: 19 pages, 3 figures, 1 table; references added; minor revisions to match JCAP version

Report number: NORDITA-2022-002; UTTG-32-2022

New constraints on the dark matter density profiles of dwarf galaxies from proper motions of globular cluster streams

Authors: Khyati Malhan, Monica Valluri, Katherine Freese, Rodrigo A. Ibata

Abstract: The central density profiles in dwarf galaxy halos depend strongly on the nature of dark matter. Recently, in Malhan et al. (2021), we employed N-body simulations to show that the cuspy cold dark matter (CDM) subhalos predicted by cosmological simulations can be differentiated from cored subhalos using the properties of accreted globular cluster (GC) streams since these GCs experience tidal stripp… ▽ More The central density profiles in dwarf galaxy halos depend strongly on the nature of dark matter. Recently, in Malhan et al. (2021), we employed N-body simulations to show that the cuspy cold dark matter (CDM) subhalos predicted by cosmological simulations can be differentiated from cored subhalos using the properties of accreted globular cluster (GC) streams since these GCs experience tidal stripping within their parent halos prior to accretion onto the Milky Way. We previously found that clusters that are accreted within cuspy subhalos produce streams with larger physical widths and higher dispersions in line-of-sight velocity and angular momentum than streams that are accreted within cored subhalos. Here, we use the same suite of simulations to demonstrate that the dispersion in the tangential velocities of streams ($σ_{v_\mathrm{Tan}}$) is also sensitive to the central DM density profiles of their parent dwarfs and GCs that were accreted from: cuspy subhalos produce streams with larger $σ_{v_\mathrm{Tan}}$ than those accreted inside cored subhalos. Using Gaia EDR3 observations of multiple GC streams we compare their $σ_{v_\mathrm{Tan}}$ values with simulations. The measured $σ_{v_\mathrm{Tan}}$ values are consistent with both an ``in situ'' origin and with accretion inside cored subhalos of $M\sim 10^{8-9}M_{\odot}$ (or very low-mass cuspy subhalos of mass $\sim 10^8M_{\odot}$). Despite the large current uncertainties in $σ_{v_\mathrm{Tan}}$, we find a low probability that any of the progenitor GCs were accreted from cuspy subhalos of $M_{\rm subhalo}\buildrel > \over \sim$ $10^9 M_{\odot}$. The uncertainties on Gaia tangential velocity measurements are expected to decrease in future and will allow for stronger constraints on subhalo DM density profiles. △ Less

Submitted 29 November, 2022; v1 submitted 10 January, 2022; originally announced January 2022.

Comments: 14 pages, 5 figures, accepted for publication in ApJ Letters

Estimating the local dark matter density in a non-axisymmetric wobbling disc

Authors: S. Sivertsson, J. I. Read, H. Silverwood, P. F. de Salas, K. Malhan, A. Widmark, C. F. P. Laporte, S. Garbari, K. Freese

Abstract: The density of dark matter near the Sun is important for experiments hunting for dark matter particles in the laboratory, and for constraining the local shape of the Milky Way's dark matter halo. Estimates to date have typically assumed that the Milky Way's stellar disc is axisymmetric and in a steady-state. Yet the Milky Way disc is neither, exhibiting prominent spiral arms and a bar, and vertica… ▽ More The density of dark matter near the Sun is important for experiments hunting for dark matter particles in the laboratory, and for constraining the local shape of the Milky Way's dark matter halo. Estimates to date have typically assumed that the Milky Way's stellar disc is axisymmetric and in a steady-state. Yet the Milky Way disc is neither, exhibiting prominent spiral arms and a bar, and vertical and radial oscillations. We assess the impact of these assumptions on determinations of the local dark matter density by applying a free-form, steady-state, Jeans method to two different N-body simulations of Milky Way-like galaxies. In one, the galaxy has experienced an ancient major merger, similar to the hypothesized Gaia-Sausage-Enceladus; in the other, the galaxy is perturbed more recently by the repeated passage and slow merger of a Sagittarius-like dwarf galaxy. We assess the impact of each of the terms in the Jeans-Poisson equations on our ability to correctly extract the local dark matter density from the simulated data. We find that common approximations employed in the literature - axisymmetry and a locally flat rotation curve - can lead to significant systematic errors of up to a factor ~1.5 in the recovered surface mass density ~2kpc above the disc plane, implying a fractional error on the local dark matter density of order unity. However, once we add in the tilt term and the rotation curve term in our models, we obtain an unbiased estimate, consistent with the true value within our 95% confidence intervals for realistic 20% uncertainties on the baryonic surface density of the disc. Other terms - the axial tilt, 2:nd Poisson and time dependent terms - contribute less than 10% to the local dark matter density (given current data) and can be safely neglected for now. In the future, as more data become available, these terms will need to be included in the analysis. △ Less

Submitted 5 January, 2022; originally announced January 2022.

Comments: 15 pages, 16 figues. Submitted to MNRAS. Comments welcome!

In-flight gain monitoring of SPIDER's transition-edge sensor arrays

Authors: J. P. Filippini, A. E. Gambrel, A. S. Rahlin, E. Y. Young, P. A. R. Ade, M. Amiri, S. J. Benton, A. S. Bergman, R. Bihary, J. J. Bock, J. R. Bond, J. A. Bonetti, S. A. Bryan, H. C. Chiang, C. R. Contaldi, O. Dore, A. J. Duivenvoorden, H. K. Eriksen, M. Farhang, A. A. Fraisse, K. Freese, M. Galloway, N. N. Gandilo, K. Ganga, R. Gualtieri , et al. (45 additional authors not shown)

Abstract: Experiments deploying large arrays of transition-edge sensors (TESs) often require a robust method to monitor gain variations with minimal loss of observing time. We propose a sensitive and non-intrusive method for monitoring variations in TES responsivity using small square waves applied to the TES bias. We construct an estimator for a TES's small-signal power response from its electrical respons… ▽ More Experiments deploying large arrays of transition-edge sensors (TESs) often require a robust method to monitor gain variations with minimal loss of observing time. We propose a sensitive and non-intrusive method for monitoring variations in TES responsivity using small square waves applied to the TES bias. We construct an estimator for a TES's small-signal power response from its electrical response that is exact in the limit of strong electrothermal feedback. We discuss the application and validation of this method using flight data from SPIDER, a balloon-borne telescope that observes the polarization of the cosmic microwave background with more than 2000 TESs. This method may prove useful for future balloon- and space-based instruments, where observing time and ground control bandwidth are limited. △ Less

Submitted 16 June, 2022; v1 submitted 1 December, 2021; originally announced December 2021.

Comments: 7 pages, 3 figures; Proceedings of the 19th International Workshop on Low Temperature Detectors (LTD19); Minor updates to match published version

Journal ref: Journal of Low Temperature Physics (2022)

Faint objects in motion: the new frontier of high precision astrometry

Authors: Fabien Malbet, Céline Boehm, Alberto Krone-Martins, Antonio Amorim, Guillem Anglada-Escudé, Alexis Brandeker, Frédéric Courbin, Torsten Enßlin, Antonio Falcão, Katherine Freese, Berry Holl, Lucas Labadie, Alain Léger, Gary Mamon, Barbara Mcarthur, Alcione Mora, Mike Shao, Alessandro Sozzetti, Douglas Spolyar, Eva Villaver, Ummi Abbas, Conrado Albertus, João Alves, Rory Barnes, Aldo Stefano Bonomo , et al. (61 additional authors not shown)

Abstract: Sky survey telescopes and powerful targeted telescopes play complementary roles in astronomy. In order to investigate the nature and characteristics of the motions of very faint objects, a flexibly-pointed instrument capable of high astrometric accuracy is an ideal complement to current astrometric surveys and a unique tool for precision astrophysics. Such a space-based mission will push the front… ▽ More Sky survey telescopes and powerful targeted telescopes play complementary roles in astronomy. In order to investigate the nature and characteristics of the motions of very faint objects, a flexibly-pointed instrument capable of high astrometric accuracy is an ideal complement to current astrometric surveys and a unique tool for precision astrophysics. Such a space-based mission will push the frontier of precision astrometry from evidence of Earth-mass habitable worlds around the nearest stars, to distant Milky Way objects, and out to the Local Group of galaxies. As we enter the era of the James Webb Space Telescope and the new ground-based, adaptive-optics-enabled giant telescopes, by obtaining these high precision measurements on key objects that Gaia could not reach, a mission that focuses on high precision astrometry science can consolidate our theoretical understanding of the local Universe, enable extrapolation of physical processes to remote redshifts, and derive a much more consistent picture of cosmological evolution and the likely fate of our cosmos. Already several missions have been proposed to address the science case of faint objects in motion using high precision astrometry missions: NEAT proposed for the ESA M3 opportunity, micro-NEAT for the S1 opportunity, and Theia for the M4 and M5 opportunities. Additional new mission configurations adapted with technological innovations could be envisioned to pursue accurate measurements of these extremely small motions. The goal of this White Paper is to address the fundamental science questions that are at stake when we focus on the motions of faint sky objects and to briefly review instrumentation and mission profiles. △ Less

Submitted 16 November, 2021; originally announced November 2021.

Comments: arXiv admin note: substantial text overlap with arXiv:1910.08028, arXiv:1707.01348

Journal ref: Experimental Astronomy, Springer Link, 2021, 51 (3), pp.845-886

A Simulation-Based Method for Correcting Mode Coupling in CMB Angular Power Spectra

Authors: J. S. -Y. Leung, J. Hartley, J. M. Nagy, C. B. Netterfield, J. A. Shariff, P. A. R. Ade, M. Amiri, S. J. Benton, A. S. Bergman, R. Bihary, J. J. Bock, J. R. Bond, J. A. Bonetti, S. A. Bryan, H. C. Chiang, C. R. Contaldi, O. Doré, A. J. Duivenvoorden, H. K. Eriksen, M. Farhang, J. P. Filippini, A. A. Fraisse, K. Freese, M. Galloway, A. E. Gambrel , et al. (45 additional authors not shown)

Abstract: Modern CMB analysis pipelines regularly employ complex time-domain filters, beam models, masking, and other techniques during the production of sky maps and their corresponding angular power spectra. However, these processes can generate couplings between multipoles from the same spectrum and from different spectra, in addition to the typical power attenuation. Within the context of pseudo-… ▽ More Modern CMB analysis pipelines regularly employ complex time-domain filters, beam models, masking, and other techniques during the production of sky maps and their corresponding angular power spectra. However, these processes can generate couplings between multipoles from the same spectrum and from different spectra, in addition to the typical power attenuation. Within the context of pseudo-$C_\ell$ based, MASTER-style analyses, the net effect of the time-domain filtering is commonly approximated by a multiplicative transfer function, $F_{\ell}$, that can fail to capture mode mixing and is dependent on the spectrum of the signal. To address these shortcomings, we have developed a simulation-based spectral correction approach that constructs a two-dimensional transfer matrix, $J_{\ell\ell'}$, which contains information about mode mixing in addition to mode attenuation. We demonstrate the application of this approach on data from the first flight of the SPIDER balloon-borne CMB experiment. △ Less

Submitted 21 April, 2022; v1 submitted 1 November, 2021; originally announced November 2021.

Comments: 14 pages, 7 figures; updated to match published version

Journal ref: ApJ 928(2):109, 2022

Natural Chain Inflation

Authors: Katherine Freese, Aliki Litsa, Martin Wolfgang Winkler

Abstract: In Chain Inflation the universe tunnels along a series of false vacua of ever-decreasing energy. The main goal of this paper is to embed Chain Inflation in high energy fundamental physics. We begin by illustrating a simple effective formalism for calculating Cosmic Microwave Background (CMB) observables in Chain Inflation. Density perturbations seeding the anisotropies emerge from the probabilisti… ▽ More In Chain Inflation the universe tunnels along a series of false vacua of ever-decreasing energy. The main goal of this paper is to embed Chain Inflation in high energy fundamental physics. We begin by illustrating a simple effective formalism for calculating Cosmic Microwave Background (CMB) observables in Chain Inflation. Density perturbations seeding the anisotropies emerge from the probabilistic nature of tunneling (rather than from quantum fluctuations of the inflation). To obtain the correct normalization of the scalar power spectrum and the scalar spectral index, we find an upper limit on the scale of inflation at horizon crossing of CMB scales, $V_*^{1/4}< 10^{12}$~GeV. We then provide an explicit realization of chain inflation, in which the inflaton is identified with an axion in supergravity. The axion enjoys a perturbative shift symmetry which is broken to a discrete remnant by instantons. The model, which we dub `natural chain inflation' satisfies all cosmological constraints and can be embedded into a standard $Λ$CDM cosmology. Our work provides a major step towards the ultraviolet completion of chain inflation in string theory. △ Less

Submitted 23 September, 2021; originally announced September 2021.

Comments: 7 pages, 4 figures

Report number: UTTG-18-2021, NORDITA-2021-078

New Projections for Dark Matter Searches with Paleo-Detectors

Authors: Sebastian Baum, Thomas D. P. Edwards, Katherine Freese, Patrick Stengel

Abstract: 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 nuc… ▽ More 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. △ Less

Submitted 11 June, 2021; originally announced June 2021.

Comments: 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

Journal ref: Instruments 2021, 5(2), 21

Using Action Space Clustering to Constrain the Accretion History of Milky Way like Galaxies

Authors: Youjia Wu, Monica Valluri, Nondh Panithanpaisal, Robyn E. Sanderson, Katherine Freese, Andrew Wetzel, Sanjib Sharma

Abstract: 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 accr… ▽ More 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. △ Less

Submitted 3 May, 2021; v1 submitted 16 April, 2021; originally announced April 2021.

Comments: 16 pages, 9 figures, fixed some typos and added two references

Report number: NORDITA-2020-094, UTTG-14-2020, LCTP-21-10

The XFaster Power Spectrum and Likelihood Estimator for the Analysis of Cosmic Microwave Background Maps

Authors: A. E. Gambrel, A. S. Rahlin, X. Song, C. R. Contaldi, P. A. R. Ade, M. Amiri, S. J. Benton, A. S. Bergman, R. Bihary, J. J. Bock, J. R. Bond, J. A. Bonetti, S. A. Bryan, H. C. Chiang, A. J. Duivenvoorden, H. K. Eriksen, M. Farhang, J. P. Filippini, A. A. Fraisse, K. Freese, M. Galloway, N. N. Gandilo, R. Gualtieri, J. E. Gudmundsson, M. Halpern , et al. (42 additional authors not shown)

Abstract: We present the XFaster analysis package. XFaster is a fast, iterative angular power spectrum estimator based on a diagonal approximation to the quadratic Fisher matrix estimator. XFaster uses Monte Carlo simulations to compute noise biases and filter transfer functions and is thus a hybrid of both Monte Carlo and quadratic estimator methods. In contrast to conventional pseudo-$C_\ell$ based method… ▽ More We present the XFaster analysis package. XFaster is a fast, iterative angular power spectrum estimator based on a diagonal approximation to the quadratic Fisher matrix estimator. XFaster uses Monte Carlo simulations to compute noise biases and filter transfer functions and is thus a hybrid of both Monte Carlo and quadratic estimator methods. In contrast to conventional pseudo-$C_\ell$ based methods, the algorithm described here requires a minimal number of simulations, and does not require them to be precisely representative of the data to estimate accurate covariance matrices for the bandpowers. The formalism works with polarization-sensitive observations and also data sets with identical, partially overlapping, or independent survey regions. The method was first implemented for the analysis of BOOMERanG data, and also used as part of the Planck analysis. Here, we describe the full, publicly available analysis package, written in Python, as developed for the analysis of data from the 2015 flight of the SPIDER instrument. The package includes extensions for self-consistently estimating null spectra and for estimating fits for Galactic foreground contributions. We show results from the extensive validation of XFaster using simulations, and its application to the SPIDER data set. △ Less

Submitted 24 May, 2021; v1 submitted 2 April, 2021; originally announced April 2021.

Comments: 18 pages, 11 figures

A Constraint on Primordial $B$-Modes from the First Flight of the SPIDER Balloon-Borne Telescope

Authors: SPIDER Collaboration, P. A. R. Ade, M. Amiri, S. J. Benton, A. S. Bergman, R. Bihary, J. J. Bock, J. R. Bond, J. A. Bonetti, S. A. Bryan, H. C. Chiang, C. R. Contaldi, O. Doré, A. J. Duivenvoorden, H. K. Eriksen, M. Farhang, J. P. Filippini, A. A. Fraisse, K. Freese, M. Galloway, A. E. Gambrel, N. N. Gandilo, K. Ganga, R. Gualtieri, J. E. Gudmundsson , et al. (46 additional authors not shown)

Abstract: We present the first linear polarization measurements from the 2015 long-duration balloon flight of SPIDER, an experiment designed to map the polarization of the cosmic microwave background (CMB) on degree angular scales. Results from these measurements include maps and angular power spectra from observations of 4.8% of the sky at 95 and 150 GHz, along with the results of internal consistency test… ▽ More We present the first linear polarization measurements from the 2015 long-duration balloon flight of SPIDER, an experiment designed to map the polarization of the cosmic microwave background (CMB) on degree angular scales. Results from these measurements include maps and angular power spectra from observations of 4.8% of the sky at 95 and 150 GHz, along with the results of internal consistency tests on these data. While the polarized CMB anisotropy from primordial density perturbations is the dominant signal in this region of sky, Galactic dust emission is also detected with high significance; Galactic synchrotron emission is found to be negligible in the SPIDER bands. We employ two independent foreground-removal techniques in order to explore the sensitivity of the cosmological result to the assumptions made by each. The primary method uses a dust template derived from Planck data to subtract the Galactic dust signal. A second approach, employing a joint analysis of SPIDER and Planck data in the harmonic domain, assumes a modified-blackbody model for the spectral energy distribution of the dust with no constraint on its spatial morphology. Using a likelihood that jointly samples the template amplitude and $r$ parameter space, we derive 95% upper limits on the primordial tensor-to-scalar ratio from Feldman-Cousins and Bayesian constructions, finding $r<0.11$ and $r<0.19$, respectively. Roughly half the uncertainty in $r$ derives from noise associated with the template subtraction. New data at 280 GHz from SPIDER's second flight will complement the Planck polarization maps, providing powerful measurements of the polarized Galactic dust emission. △ Less

Submitted 24 March, 2021; originally announced March 2021.

Comments: 29 pages, 13 figures

Chain Early Dark Energy: Solving the Hubble Tension and Explaining Today's Dark Energy

Authors: Katherine Freese, Martin Wolfgang Winkler

Abstract: We propose a new model of Early Dark Energy (EDE) as a solution to the Hubble tension in cosmology, the apparent discrepancy between local measurements of the Hubble constant $H_0\simeq 74$ km s$^{-1}$ Mpc$^{-1}$ and $H_0\simeq 67$ km s$^{-1}$ Mpc$^{-1}$ inferred from the Cosmic Microwave Background (CMB). In Chain EDE, the Universe undergoes a series of first order phase transitions, starting at… ▽ More We propose a new model of Early Dark Energy (EDE) as a solution to the Hubble tension in cosmology, the apparent discrepancy between local measurements of the Hubble constant $H_0\simeq 74$ km s$^{-1}$ Mpc$^{-1}$ and $H_0\simeq 67$ km s$^{-1}$ Mpc$^{-1}$ inferred from the Cosmic Microwave Background (CMB). In Chain EDE, the Universe undergoes a series of first order phase transitions, starting at a high energy vacuum in a potential, and tunneling down through a chain of every lower energy metastable minima. As in all EDE models, the contribution of the vacuum energy to the total energy density of the universe is initially negligible, but reaches $\sim 10\%$ around matter-radiation equality, before cosmological data require it to redshift away quickly -- at least as fast as radiation. We indeed obtain this required behavior with a series of $N$ tunneling events, and show that for $N>600$ the phase transitions are rapid enough to allow fast percolation and thereby avoid large scale anisotropies in the CMB. We construct a specific example of Chain EDE featuring a scalar field in a quasiperiodic potential (a tilted cosine), which is ubiquitous in axion physics and, therefore, carries strong theoretical motivation. Interestingly, the energy difference between vacua can be roughly the size of today's Dark Energy (meV scale). Therefore, the end result of Chain EDE could provide a natural explanation of Dark Energy, if the tunneling becomes extremely slow in the final step before the field reaches zero (or negative) energy. We discuss a simple mechanism which can stop the scalar field in the desired minimum. Thus Chain EDE offers the exciting prospect to explain EDE and Dark Energy by the same scalar field. △ Less

Submitted 19 April, 2021; v1 submitted 26 February, 2021; originally announced February 2021.

Comments: 20 pages, 6 figures

Report number: UTTG-02-2021, NORDITA-2021-020

Inelastic dark matter scattering off Thallium cannot save DAMA

Authors: Sunniva Jacobsen, Katherine Freese, Chris Kelso, Pearl Sandick, Patrick Stengel

Abstract: We study the compatibility of the observed DAMA modulation signal with inelastic scattering of dark matter (DM) off of the $0.1\%$ Thallium (Tl) dopant in DAMA. In this work we test whether there exist regions of parameter space where the Tl interpretation gives a good fit to the most recent data from DAMA, and whether these regions are compatible with the latest constraints from other direct dete… ▽ More We study the compatibility of the observed DAMA modulation signal with inelastic scattering of dark matter (DM) off of the $0.1\%$ Thallium (Tl) dopant in DAMA. In this work we test whether there exist regions of parameter space where the Tl interpretation gives a good fit to the most recent data from DAMA, and whether these regions are compatible with the latest constraints from other direct detection experiments. Previously, Chang et al. in 2010, had proposed the Tl interpretation of the DAMA data, and more recently (in 2019) the DAMA/LIBRA collaboration found regions in parameter space of Tl inelastic scattering that differ by more than $10σ$ from a no modulation hypothesis. We have expanded upon their work by testing whether the regions of parameter space where inelastic DM-Tl scattering gives a good fit to the most recent DAMA data survive the constraints placed by the lack of a DM signal in XENON1T and CRESST-II. In addition, we have tested how these regions change with the main sources of uncertainty: the Tl quenching factor, which has never been measured directly, and the astrophysical uncertainties in the DM distribution. We conclude that inelastic DM scattering off Tl cannot explain the DAMA data in light of null results from other experiments. △ Less

Submitted 17 March, 2021; v1 submitted 16 February, 2021; originally announced February 2021.

Comments: Added references in sections 1 and 2. Added a comment on the recent results from ANAIS-112

Design and pre-flight performance of SPIDER 280 GHz receivers

Authors: E. C. Shaw, P. A. R. Ade, S. Akers, M. Amiri, J. Austermann, J. Beall, D. T. Becker, S. J. Benton, A. S. Bergman, J. J. Bock, J. R. Bond, S. A. Bryan, H. C. Chiang, C. R. Contaldi, R. S. Domagalski, O. Doré, S. M. Duff, A. J. Duivenvoorden, H. K. Eriksen, M. Farhang, J. P. Filippini, L. M. Fissel, A. A. Fraisse, K. Freese, M. Galloway , et al. (57 additional authors not shown)

Abstract: In this work we describe upgrades to the Spider balloon-borne telescope in preparation for its second flight, currently planned for December 2021. The Spider instrument is optimized to search for a primordial B-mode polarization signature in the cosmic microwave background at degree angular scales. During its first flight in 2015, Spider mapped ~10% of the sky at 95 and 150 GHz. The payload for th… ▽ More In this work we describe upgrades to the Spider balloon-borne telescope in preparation for its second flight, currently planned for December 2021. The Spider instrument is optimized to search for a primordial B-mode polarization signature in the cosmic microwave background at degree angular scales. During its first flight in 2015, Spider mapped ~10% of the sky at 95 and 150 GHz. The payload for the second Antarctic flight will incorporate three new 280 GHz receivers alongside three refurbished 95- and 150 GHz receivers from Spider's first flight. In this work we discuss the design and characterization of these new receivers, which employ over 1500 feedhorn-coupled transition-edge sensors. We describe pre-flight laboratory measurements of detector properties, and the optical performance of completed receivers. These receivers will map a wide area of the sky at 280 GHz, providing new information on polarized Galactic dust emission that will help to separate it from the cosmological signal. △ Less

Submitted 22 December, 2020; originally announced December 2020.

Comments: 13 pages, 8 figures; as published in the conference proceedings for SPIE Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy X (2020)

Journal ref: "Design and pre-flight performance of SPIDER 280 GHz receivers," Proc. SPIE 11453, Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy X, 114532F (13 December 2020)

The Power Spectrum of Density Perturbations in Chain Inflation

Authors: Martin Wolfgang Winkler, Katherine Freese

Abstract: 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 $Λ$CDM cosmol… ▽ More 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 $Λ$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. △ Less

Submitted 25 November, 2020; originally announced November 2020.

Comments: 19 pages, 9 figures

Report number: NORDITA-2020-107, UTTG 22-2020

Journal ref: Phys. Rev. D 103, 043511 (2021)

Primordial non-Gaussianity from the Effects of the Standard Model Higgs during Reheating after Inflation

Authors: Aliki Litsa, Katherine Freese, Evangelos I. Sfakianakis, Patrick Stengel, Luca Visinelli

Abstract: We propose a new way of studying the Higgs potential at extremely high energies. The Standard Model (SM) Higgs boson, as a light spectator field during inflation in the early Universe, can acquire large field values from its quantum fluctuations which vary among different causal (Hubble) patches. Such a space dependence of the Higgs after the end of inflation leads to space-dependent SM particle m… ▽ More We propose a new way of studying the Higgs potential at extremely high energies. The Standard Model (SM) Higgs boson, as a light spectator field during inflation in the early Universe, can acquire large field values from its quantum fluctuations which vary among different causal (Hubble) patches. Such a space dependence of the Higgs after the end of inflation leads to space-dependent SM particle masses and hence variable efficiency of reheating, when the inflaton decays to Higgsed SM particles. Inhomogeneous reheating results in (observable) temperature anisotropies. Further, the resulting temperature anisotropy spectrum acquires a significant non-Gaussian component, which is constrained by $\textit{Planck}$ observations of the Cosmic Microwave Background (CMB) and potentially detectable in next-generation experiments. Constraints on this non-Gaussian signal largely exclude the possibility of the observed temperature anisotropies arising primarily from Higgs effects. Hence, in principle, observational searches for non-Gaussianity in the CMB can be used to constrain the dynamics of the Higgs boson at very high (inflationary) energies. △ Less

Submitted 15 March, 2023; v1 submitted 23 November, 2020; originally announced November 2020.

Comments: 7 pages, 3 figures, joint submission with arXiv:2009.14218, figures added, important changes, version published in JCAP

Report number: NORDITA-2020-083; UTTG-17-2020; Nikhef-2020-037

Journal ref: JCAP 2303, 033 (2023)

Stability and pulsation of the first dark stars

Authors: Tanja Rindler-Daller, Katherine Freese, Richard H. D. Townsend, Luca Visinelli

Abstract: 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… ▽ More 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 $κ-γ$ and $γ$ 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. △ Less

Submitted 12 February, 2021; v1 submitted 31 October, 2020; originally announced November 2020.

Comments: minor revisions; accepted by MNRAS; 15 pages, 6 Figures

Report number: NORDITA-2020-078; LCTP-20-21; UTTG-20-2020

Journal ref: Mon. Not. Roy. Astron. Soc. 503 3, 3677-3691 (2021)

Filling the Black Hole Mass Gap: Avoiding Pair Instability in Massive Stars through Addition of Non-Nuclear Energy

Authors: Joshua Ziegler, Katherine Freese

Abstract: 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… ▽ More 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$. △ Less

Submitted 1 October, 2020; originally announced October 2020.

Comments: 20 pages, 9 figures, 2 tables

Report number: UTTG-12-2020, NORDITA 2020-088

Journal ref: Phys. Rev. D 104, 043015 (2021)

Large Density Perturbations from Reheating to Standard Model particles due to the Dynamics of the Higgs Boson during Inflation

Authors: Aliki Litsa, Katherine Freese, Evangelos I. Sfakianakis, Patrick Stengel, Luca Visinelli

Abstract: Cosmic Microwave Background (CMB) observations are used to constrain reheating to Standard Model (SM) particles after a period of inflation. As a light spectator field, the SM Higgs boson acquires large field values from its quantum fluctuations during inflation, gives masses to SM particles that vary from one Hubble patch to another, and thereby produces large density fluctuations. We consider bo… ▽ More Cosmic Microwave Background (CMB) observations are used to constrain reheating to Standard Model (SM) particles after a period of inflation. As a light spectator field, the SM Higgs boson acquires large field values from its quantum fluctuations during inflation, gives masses to SM particles that vary from one Hubble patch to another, and thereby produces large density fluctuations. We consider both perturbative and resonant decay of the inflaton to SM particles. For the case of perturbative decay from coherent oscillations of the inflaton after high scale inflation, we find strong constraints on the reheat temperature for the inflaton decay into heavy SM particles. For the case of resonant particle production (preheating) to (Higgsed) SM gauge bosons, we find temperature fluctuations larger than observed in the CMB for a range of gauge coupling that includes those found in the SM and conclude that such preheating cannot be the main source of reheating the Universe after inflation. △ Less

Submitted 21 December, 2020; v1 submitted 29 September, 2020; originally announced September 2020.

Comments: 24 pages, 10 figures, references added, minor changes

Report number: UT-03-2020; UTTG-11-2020; Nikhef 2020-029; NORDITA-2020-082

Journal ref: Phys. Rev. D 104 12, 123546 (2021)

Cornering (quasi) degenerate neutrinos with cosmology

Authors: Massimiliano Lattanzi, Martina Gerbino, Katherine Freese, Gordon Kane, José W. F. Valle

Abstract: 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 o… ▽ More 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 $σ$'s, while adding Baryon acoustic oscillations (BAO) data brings the rejection to 5.9$σ$'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 $Σm_ν= 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$σ$. % We discuss the robustness of these projections with respect to assumptions on the underlying cosmological model, and also compare them with bounds from $β$ decay endpoint and neutrinoless double beta decay studies. △ Less

Submitted 3 July, 2020; originally announced July 2020.

Comments: 19 pages, 6 figures, 1 table

Report number: IFIC/20-XXX

Journal ref: JHEP 10 (2020) 213

The First Three Seconds: a Review of Possible Expansion Histories of the Early Universe

Authors: Rouzbeh Allahverdi, Mustafa A. Amin, Asher Berlin, Nicolás Bernal, Christian T. Byrnes, M. Sten Delos, Adrienne L. Erickcek, Miguel Escudero, Daniel G. Figueroa, Katherine Freese, Tomohiro Harada, Dan Hooper, David I. Kaiser, Tanvi Karwal, Kazunori Kohri, Gordan Krnjaic, Marek Lewicki, Kaloian D. Lozanov, Vivian Poulin, Kuver Sinha, Tristan L. Smith, Tomo Takahashi, Tommi Tenkanen, James Unwin, Ville Vaskonen , et al. (1 additional authors not shown)

Abstract: 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 theref… ▽ More 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. △ Less

Submitted 27 January, 2021; v1 submitted 29 June, 2020; originally announced June 2020.

Comments: 67 pages, 18 figures. v2: Discussion and references added. Accepted for publication in The Open Journal of Astrophysics

Report number: FERMILAB-PUB-20-242-A, KCL-PH-TH/2020-33, KEK-Cosmo-257, KEK-TH-2231, IPMU20-0070, PI/UAN-2020-674FT, RUP-20-22

Journal ref: Open J. Astrophys. Vol. 4, 2021

Probing the nature of dark matter with accreted globular cluster streams

Authors: Khyati Malhan, Monica Valluri, Katherine Freese

Abstract: 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… ▽ More 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. △ Less

Submitted 16 December, 2020; v1 submitted 26 May, 2020; originally announced May 2020.

Comments: 24 pages, 15 figures, 1 table, accepted for publication in MNRAS

Bounds on light sterile neutrino mass and mixing from cosmology and laboratory searches

Authors: Steffen Hagstotz, Pablo F. de Salas, Stefano Gariazzo, Martina Gerbino, Massimiliano Lattanzi, Sunny Vagnozzi, Katherine Freese, Sergio Pastor

Abstract: We provide a consistent framework to set limits on properties of light sterile neutrinos coupled to all three active neutrinos using a combination of the latest cosmological data and terrestrial measurements from oscillations, $β$-decay and neutrinoless double-$β$ decay ($0νββ$) experiments. We directly constrain the full $3+1$ active-sterile mixing matrix elements $|U_{\alpha4}|^2$, with… ▽ More We provide a consistent framework to set limits on properties of light sterile neutrinos coupled to all three active neutrinos using a combination of the latest cosmological data and terrestrial measurements from oscillations, $β$-decay and neutrinoless double-$β$ decay ($0νββ$) experiments. We directly constrain the full $3+1$ active-sterile mixing matrix elements $|U_{\alpha4}|^2$, with $α\in ( e,μ,τ)$, and the mass-squared splitting $Δm^2_{41} \equiv m_4^2-m_1^2$. We find that results for a $3+1$ case differ from previously studied $1+1$ scenarios where the sterile is only coupled to one of the neutrinos, which is largely explained by parameter space volume effects. Limits on the mass splitting and the mixing matrix elements are currently dominated by the cosmological data sets. The exact results are slightly prior dependent, but we reliably find all matrix elements to be constrained below $|U_{\alpha4}|^2 \lesssim 10^{-3}$. Short-baseline neutrino oscillation hints in favor of eV-scale sterile neutrinos are in serious tension with these bounds, irrespective of prior assumptions. We also translate the bounds from the cosmological analysis into constraints on the parameters probed by laboratory searches, such as $m_β$ or $m_{ββ}$, the effective mass parameters probed by $β$-decay and $0νββ$ searches, respectively. When allowing for mixing with a light sterile neutrino, cosmology leads to upper bounds of $m_β< 0.09$ eV and $m_{ββ} < 0.07$ eV at 95\% C.L, more stringent than the limits from current laboratory experiments. △ Less

Submitted 4 March, 2020; originally announced March 2020.

Comments: 29 pages, 7 figures

Journal ref: Phys. Rev. D 104 (2021) 123524

CMB B-mode non-Gaussianity: optimal bispectrum estimator and Fisher forecasts

Authors: Adriaan J. Duivenvoorden, P. Daniel Meerburg, Katherine Freese

Abstract: 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. <TTB>, <EEB>, and <TEB>. Such bispectra probe violations of the tensor consisten… ▽ More 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. <TTB>, <EEB>, and <TEB>. 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. △ Less

Submitted 26 November, 2019; originally announced November 2019.

Comments: 40 pages, 7 figures

Journal ref: Phys. Rev. D 102, 023521 (2020)