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

Measuring the cosmic ray flux over timescales comparable to the age of the solar system, $\sim 4.5\,$Gyr, could provide a new window on the history of the Earth, the solar system, and even our galaxy. We present a technique to indirectly measure the rate of cosmic rays as a function of time using the imprints of atmospheric neutrinos in paleo-detectors, natural minerals which record damage tracks from nuclear recoils. Minerals commonly found on Earth are $\lesssim 1\,$Gyr old, providing the ability to look back across cosmic ray history on timescales of the same order as the age of the solar system. Given a collection of differently aged samples dated with reasonable accuracy, this technique is particularly well-suited to measuring historical changes in the cosmic ray flux at Earth and is broadly applicable in astrophysics and geophysics.

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

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

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

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

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

Many extensions of the Standard Model predict large numbers of additional unstable particles whose decays in the early universe are tightly constrained by observational data. For example, the decays of such particles can alter the ratios of light-element abundances, give rise to distortions in the cosmic microwave background, alter the ionization history of the universe, and contribute to the diffuse photon flux. Constraints on new physics from such considerations are typically derived for a single unstable particle species with a single well-defined mass and characteristic lifetime. In this paper, by contrast, we investigate the cosmological constraints on theories involving entire ensembles of decaying particles --- ensembles which span potentially broad ranges of masses and lifetimes. In addition to providing a detailed numerical analysis of these constraints, we also formulate a set of simple analytic approximations for these constraints which may be applied to generic ensembles of unstable particles which decay into electromagnetically-interacting final states. We then illustrate how these analytic approximations can be used to constrain a variety of toy scenarios for physics beyond the Standard Model. For ease of reference, we also compile our results in the form of a table which can be consulted independently of the rest of the paper. It is thus our hope that this work might serve as a useful reference for future model-builders concerned with cosmological constraints on decaying particles, regardless of the particular model under study.

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

We investigate constraints on the interactions of light dark matter with Standard Model quarks in a framework with effective contact operators mediating the decay of heavy flavor bound state quarkonium to dark matter and a photon. When considered in combination with decays to purely invisible final states, constraints from heavy quarkonium decays at high intensity electron-positron colliders can complement missing energy searches at high energy colliders and provide sensitivity to dark matter masses difficult to probe at direct and indirect detection experiments. We calculate the approximate limits on the branching fraction for $\Upsilon (1 S)$ decays to dark matter and a photon. Given the approximate limits on the branching fractions for all dimension 6 or lower contact operators, we present the corresponding limits on the interaction strength for each operator and the inferred limits on dark matter-nucleon scattering. Complementary constraints on dark matter annihilation from gamma-ray searches from dwarf spheroidal galaxies are also considered.

We investigate constraints on the properties of light dark matter which can be obtained from analysis of invisible quarkonium decays at high intensity electron-positron colliders in the framework of a low energy effective field theory. A matrix element analysis of all contact operators pertinent for these meson decays allows for a model-independent calculation of associated dark matter-nucleon scattering and dark matter annihilation cross sections. Assuming dark matter couples universally to all quark flavors, we then obtain bounds on nucleon scattering which complement direct dark matter detection searches. In contrast to similar analyses of monojet searches at high energy colliders, B and charm factories are more suitable probes of light dark matter interactions with less massive mediators. Relevant bounds on dark matter annihilation arising from gamma ray searches of dwarf spheroidal galaxies are also presented.

We consider a class of leptogenesis models in which the lepton asymmetry arises from dark matter annihilation processes which violate CP and lepton number. Importantly, a necessary one-loop contribution to the annihilation matrix element arises from absorptive final state interactions. We elucidate the relationship between this one-loop contribution and the CP-violating phase. As we show, the branching fraction for dark matter annihilation to leptons may be small in these models, while still generating the necessary asymmetry.

We consider a class of leptogenesis models in which the lepton asymmetry arises from dark matter annihilation processes which violate CP and lepton number. Importantly, a necessary one-loop contribution to the annihilation matrix element arises from absorptive final state interactions. We elucidate the relationship between this one-loop contribution and the CP-violating phase. As we show, the branching fraction for dark matter annihilation to leptons may be small in these models, while still generating the necessary asymmetry.