The identity of dark matter is being sought with increasingly sensitive and voluminous underground detectors. Recently the XENON1T collaboration reported excess electronic recoil events, with most of these having recoil energies around $1-30$ keV. We show that a straightforward model of inelastic dark matter produced via early universe thermal freeze-out annihilation can account for the XENON1T excess. Remarkably, this dark matter model consists of a few simple elements: sub-GeV mass Dirac fermion dark matter coupled to a lighter dark photon kinetically mixed with the Standard Model photon. A scalar field charged under the dark U(1) gauge symmetry can provide a mass for the dark photon and splits the Dirac fermion component state masses by a few keV, which survive in equal abundance and interact inelastically with electrons and nuclei.

We propose a matter effect for the gravitational ultraviolet/infrared (UV/IR) mixing solution to the cosmological constant problem. Previously, the gravitational UV/IR mixing model implied a non-standard equation of state for dark energy, contradicting observation. In contrast, matter effect gravitational UV/IR mixing accommodates a standard $\Lambda$CDM cosmology with constant dark energy. Notably, there are new density-dependent predictions for futuristically precise measurements of fundamental parameters, like the magnetic moments of the muon and electron.

Neutron stars serve as excellent next-generation thermal detectors of dark matter, heated by the scattering and annihilation of dark matter accelerated to relativistic speeds in their deep gravitational wells. However, the dynamics of neutron star cores are uncertain, making it difficult at present to unequivocally compute dark matter scattering in this region. On the other hand, the physics of an outer layer of the neutron star, the crust, is more robustly understood. We show that dark matter scattering solely with the low-density crust still kinetically heats neutron stars to infrared temperatures detectable by forthcoming telescopes. We find that for both spin-independent and spin-dependent scattering on nucleons, the crust-only cross section sensitivity is $10^{-43} - 10^{-41}$~cm$^2$ for dark matter masses of 100 MeV $-$ 1 PeV, with the best sensitivity arising from dark matter scattering with a crust constituent called nuclear pasta (including gnocchi, spaghetti, and lasagna phases). For dark matter masses from 10 eV to 1 MeV, the sensitivity is $10^{-39} - 10^{-34}$~cm$^2$, arising from exciting collective phonon modes in a neutron superfluid in the inner crust. Furthermore, for any $s$-wave or $p$-wave annihilating dark matter, we show that dark matter will efficiently annihilate by thermalizing just with the neutron star crust, regardless of whether the dark matter ever scatters with the neutron star core. This implies efficient annihilation in neutron stars for any electroweakly interacting dark matter with inelastic mass splittings of up to 200 MeV, including Higgsinos. We conclude that neutron star crusts play a key role in dark matter scattering and annihilation in neutron stars.

If dark matter is efficiently captured by a planet, energy released in its annihilation can exceed that planet's total heat output. Building on prior work, we treat Earth's composition and dark matter capture in detail and present improved limits on dark matter-nucleon scattering cross sections for dark matter masses ranging from 0.1 to $10^{10}$ GeV. We also extend Earth limits by applying the same treatment to Mars. The scope of dark matter models considered is expanded to include spin-dependent nuclear interactions including isospin-independent, proton only, and neutron only interactions. We find that Earth and Mars heating bounds are alleviated for dark matter s-wave self-annihilation cross sections $\lesssim 10^{-44}~{\rm cm^2}$.

We show that current and imminent underground detectors are capable of precision astrometry of dark matter. First we show that galactic dark matter velocity distributions can be obtained from reconstructed tracks of dark matter scattering on multiple nuclei during transit; using the liquid scintillator neutrino detector SNO+ as an example, we find that the dark matter velocity vector can be reconstructed event-by-event with such a small uncertainty, that the precision of dark matter astrometry will be limited mainly by statistics. We then determine the number of dark matter events required to determine the dispersion speed, escape speed, and velocity anisotropies of the local dark matter halo, and also find that with as few as $\mathcal{O}(10)$ events, dark matter signals may be discriminated from potential backgrounds arising as power-law distributions. Finally, we discuss the prospects of dark matter astrometry at other liquid scintillator detectors, dark matter experiments, and the recently proposed MATHUSLA detector.

It was recently demonstrated that asymmetric dark matter can ignite supernovae by collecting and collapsing inside lone sub-Chandrasekhar mass white dwarfs, and that this may be the cause of Type Ia supernovae. A ball of asymmetric dark matter accumulated inside a white dwarf and collapsing under its own weight, sheds enough gravitational potential energy through scattering with nuclei, to spark the fusion reactions that precede a Type Ia supernova explosion. In this article we elaborate on this mechanism and use it to place new bounds on interactions between nucleons and asymmetric dark matter for masses $m_{X} = 10^{6}-10^{16}$ GeV. Interestingly, we find that for dark matter more massive than $10^{11}$ GeV, Type Ia supernova ignition can proceed through the Hawking evaporation of a small black hole formed by the collapsed dark matter. We also identify how a cold white dwarf's Coulomb crystal structure substantially suppresses dark matter-nuclear scattering at low momentum transfers, which is crucial for calculating the time it takes dark matter to form a black hole. Higgs and vector portal dark matter models that ignite Type Ia supernovae are explored.

We show that dark matter with a per-nucleon scattering cross section $\gtrsim 10^{-28}~{\rm cm^2}$ could be discovered by liquid scintillator neutrino detectors like BOREXINO, SNO+, and JUNO. Due to the large dark matter fluxes admitted, these detectors could find dark matter with masses up to $10^{21}$ GeV, surpassing the mass sensitivity of current direct detection experiments (such as XENON1T and PICO) by over two orders of magnitude. We derive the spin-independent and spin-dependent cross section sensitivity of these detectors using existing selection triggers, and propose an improved trigger program that enhances this sensitivity by two orders of magnitude. We interpret these sensitivities in terms of three dark matter scenarios: (1) effective contact operators for scattering, (2) QCD-charged dark matter, and (3) a recently proposed model of Planck-mass baryon-charged dark matter. We calculate the flux attenuation of dark matter at these detectors due to the earth overburden, taking into account the earth's density profile and elemental composition, and nuclear spins.

We demonstrate that dark matter heating of gas clouds hundreds of parsecs from the Milky Way Galactic center provides a powerful new test of dark matter interactions. To illustrate, we set a leading bound on nucleon scattering for 10-100 MeV mass dark matter. We also constrain millicharged dark matter models, including those proposed to match the recent EDGES 21 cm absorption anomaly. For Galactic center gas clouds, galactic fields' magnetic deflection of electromagnetically charged dark matter is mitigated, because the magnetic fields around the Galactic center are poloidal, as opposed to being aligned parallel to the Milky Way disk. We discuss prospects for detecting dark matter using a population of Galactic center gas clouds warmed by dark matter.

We identify new astrophysical signatures of dark matter that implodes neutron stars (NSs), which could decisively test whether NS-imploding dark matter is responsible for missing pulsars in the Milky Way galactic center, the source of some $r$-process elements, and the origin of fast-radio bursts. First, NS-imploding dark matter forms $\sim 10^{-10}$ solar mass or smaller black holes inside neutron stars, which proceed to convert neutron stars into $\sim$1.5 solar mass BHs. This decreases the number of neutron star mergers seen by LIGO/Virgo (LV) and associated merger kilonovae seen by telescopes like DES, BlackGEM, and ZTF, and instead, producing a population of "black mergers" containing $\sim$1.5 solar mass black holes. Second, dark matter-induced neutron star implosions may create a new kind of kilonovae that lacks a detectable, accompanying gravitational signal, which we call "quiet kilonovae." Using DES data and the Milky Way's r-process abundance, we constrain quiet kilonovae. Third, the spatial distribution of neutron star merger kilonovae and quiet kilonovae in galaxies can be used to detect dark matter. NS-imploding dark matter destroys most neutron stars at the centers of disc galaxies, so that neutron star merger kilonovae would appear mostly in a donut at large radii. We find that as few as ten neutron star merger kilonova events, located to $\sim$1 kpc precision could validate or exclude dark matter-induced neutron star implosions at $2 \sigma$ confidence, exploring dark matter-nucleon cross-sections 4-10 orders of magnitude below current direct detection experimental limits. Similarly, NS-imploding dark matter as the source of fast radio bursts can be tested at $2 \sigma$ confidence once 20 bursts are located in host galaxies by radio arrays like CHIME and HIRAX.

We identify a largely model-independent signature of dark matter interactions with nucleons and electrons. Dark matter in the local galactic halo, gravitationally accelerated to over half the speed of light, scatters against and deposits kinetic energy into neutron stars, heating them to infrared blackbody temperatures. The resulting radiation could potentially be detected by the James Webb Space Telescope, the Thirty Meter Telescope, or the European Extremely Large Telescope. This mechanism also produces optical emission from neutron stars in the galactic bulge, and X-ray emission near the galactic center, because dark matter is denser in these regions. For GeV - PeV mass dark matter, dark kinetic heating would initially unmask any spin-independent or spin-dependent dark matter-nucleon cross-sections exceeding $2 \times 10^{-45}$ cm$^2$, with improved sensitivity after more telescope exposure. For lighter-than-GeV dark matter, cross-section sensitivity scales inversely with dark matter mass because of Pauli blocking; for heavier-than-PeV dark matter, it scales linearly with mass as a result of needing multiple scatters for capture. Future observations of dark sector-warmed neutron stars could determine whether dark matter annihilates in or only kinetically heats neutron stars. Because inelastic inter-state transitions of up to a few GeV would occur in relativistic scattering against nucleons, elusive inelastic dark matter like pure Higgsinos can also be discovered.

Observations by HAWC and Milagro have detected bright and spatially extended TeV gamma-ray sources surrounding the Geminga and Monogem pulsars. We argue that these observations, along with a substantial population of other extended TeV sources coincident with pulsar wind nebulae, constitute a new morphological class of spatially extended TeV halos. We show that HAWCs wide field-of-view unlocks an expansive parameter space of TeV halos not observable by atmospheric Cherenkov telescopes. Under the assumption that Geminga and Monogem are typical middle-aged pulsars, we show that ten-year HAWC observations should eventually observe 37$^{+17}_{-13}$ middle-aged TeV halos that correspond to pulsars whose radio emission is not beamed towards Earth. Depending on the extrapolation of the TeV halo efficiency to young pulsars, HAWC could detect more than 100 TeV halos from mis-aligned pulsars. These pulsars have historically been difficult to detect with existing multiwavelength observations. TeV halos will constitute a significant fraction of all HAWC sources, allowing follow-up observations to efficiently find pulsar wind nebulae and thermal pulsar emission. The observation and subsequent multi-wavelength follow-up of TeV halos will have significant implications for our understanding of pulsar beam geometries, the evolution of PWN, the diffusion of cosmic-rays near energetic pulsars, and the contribution of pulsars to the cosmic-ray positron excess.

The early universe could feature multiple reheating events, leading to jumps in the visible sector entropy density that dilute both particle asymmetries and the number density of frozen-out states. In fact, late time entropy jumps are usually required in models of Affleck-Dine baryogenesis, which typically produces an initial particle-antiparticle asymmetry that is much too large. An important consequence of late time dilution, is that a smaller dark matter annihilation cross section is needed to obtain the observed dark matter relic density. For cosmologies with high scale baryogenesis, followed by radiation-dominated dark matter freeze-out, we show that the perturbative unitarity mass bound on thermal relic dark matter is relaxed to $10^{10}$ GeV. We proceed to study superheavy asymmetric dark matter models, made possible by a sizable entropy injection after dark matter freeze-out, and identify how the Affleck-Dine mechanism would generate the baryon and dark asymmetries.

Recent observations of Reticulum II have uncovered an overabundance of r-process elements, compared to similar ultra-faint dwarf spheroidal galaxies (UFDs). Because the metallicity and star formation history of Reticulum II appear consistent with all known UFDs, the high r-process abundance of Reticulum II suggests enrichment through a single, rare event, such as a double neutron star (NS) merger. However, we note that this scenario is extremely unlikely, as binary stellar evolution models require significant supernova natal kicks to produce NS-NS or NS-black hole mergers, and these kicks would efficiently remove compact binary systems from the weak gravitational potentials of UFDs. We examine alternative mechanisms for the production of r-process elements in UFDs, including a novel mechanism wherein NSs in regions of high dark matter density implode after accumulating a black-hole-forming mass of dark matter. We find that r-process proto-material ejection by tidal forces, when a single neutron star implodes into a black hole, can occur at a rate matching the r-process abundance of both Reticulum II and the Milky Way. Remarkably, dark matter models which collapse a single neutron star in observed UFDs also solve the missing pulsar problem in the Milky Way Galactic center. We propose tests specific to dark matter r-process production which may uncover, or rule out, this model.

Inflation occurring at energy densities less than (10$^{14}$ GeV)$^4$ produces tensor perturbations too small to be measured by cosmological surveys. However, we show that it is possible to probe low scale inflation by measuring the mass of the inflaton at low energy experiments. Detection prospects and cosmological constraints are determined for low scale quartic hilltop models of inflation paired with a curvaton field, which imprints the spectrum of scalar perturbations observed in large scale structure and on the cosmic microwave background. With cosmological constraints applied, low scale quartic inflation at energies GeV--PeV, can be mapped to an MeV--TeV mass inflaton resonance, discoverable through a Higgs portal coupling at upcoming collider and meson decay experiments. It is demonstrated that low scale inflatons can have detectably large couplings to Standard Model particles through a Higgs portal, permitting prompt reheating after inflation, without spoiling, through radiative corrections to the inflaton's self-coupling, the necessary flatness of a low scale inflationary potential. A characteristic particle spectrum for a quartic inflaton-curvaton pair is identified: to within an order of magnitude, the mass of the curvaton can be predicted from the mass of the inflaton, and vice-versa. Low scale inflation Higgs portal sensitivity targets are found for experiments like the LHC, SHiP, BEPC, and KEKB.

Recent studies of low redshift type Ia supernovae (SNIa) indicate that half explode from less than Chandrasekhar mass white dwarfs, implying ignition must proceed from something besides the canonical criticality of Chandrasekhar mass SNIa progenitors. We show that $1-100$ PeV mass asymmetric dark matter, with imminently detectable nucleon scattering interactions, can accumulate to the point of self-gravitation in a white dwarf and collapse, shedding gravitational potential energy by scattering off nuclei, thereby heating the white dwarf and igniting the flame front that precedes SNIa. We combine data on SNIa masses with data on the ages of SNIa-adjacent stars. This combination reveals a $2.8 \sigma$ inverse correlation between SNIa masses and ignition ages, which could result from increased capture of dark matter in 1.4 versus 1.1 solar mass white dwarfs. Future studies of SNIa in galactic centers will provide additional tests of dark-matter-induced type Ia ignition. Remarkably, both bosonic and fermionic SNIa-igniting dark matter also resolve the missing pulsar problem by forming black holes in $\gtrsim 10$ Myr old pulsars at the center of the Milky Way.

If forthcoming measurements of cosmic photon polarization restrict the primordial tensor-to-scalar ratio to $r < 0.01$, small field inflation will be a principal candidate for the origin of the universe. Here we show that small multifield inflation, without the hybrid mechanism, typically results in large squeezed nongaussianity. Small multifield potentials contain multiple flat field directions, often identified with the gauge invariant field directions in supersymmetric potentials. We find that unless these field directions have equal slopes, large nongaussianity arises. After identifying relevant differences between large and small two-field potentials, we demonstrate that the latter naturally fulfill the Byrnes-Choi-Hall large nongaussianity conditions. Computations of the primordial power spectrum, spectral index, and squeezed bispectrum, reveal that small two-field models which otherwise match observed primordial perturbations, produce excludably large nongaussianity if the inflatons' field directions have unequal slopes.

Pulsars apparently missing from the galactic center could have been destroyed by asymmetric fermionic dark matter ($m_X = 1-100$ GeV) coupled to a light scalar ($m_{\phi}= 5-20$ MeV), which mixes with the Higgs boson. We point out that this pulsar-collapsing dark sector can resolve the core-cusp problem and will either be excluded or discovered by upcoming direct detection experiments. Another implication is a maximum pulsar age curve that increases with distance from the galactic center, with a normalization that depends on the couplings and masses of dark sector particles. In addition, we use old pulsars outside the galactic center to place bounds on asymmetric Higgs portal models.

We map the parameter space for MSSM neutralino dark matter which freezes out to the observed relic abundance, in the limit that all superpartners except the neutralinos and charginos are decoupled. In this space of relic neutralinos, we show the dominant dark matter annihilation modes, the mass splittings among the electroweakinos, direct detection rates, and collider cross-sections. The mass difference between the dark matter and the next-to-lightest neutral and charged states is typically much less than electroweak gauge boson masses. With these small mass differences, the relic neutralino surface is accessible to a future 100 TeV hadron collider, which can discover inter-neutralino mass splittings down to 1 GeV and thermal relic dark matter neutralino masses up to 1.5 TeV with a few inverse attobarns of luminosity. This coverage is a direct consequence of the increased collider energy: the Standard Model events with missing transverse momentum in the TeV range have mostly hard electroweak radiation, distinct from the soft radiation shed in compressed electroweakino decays. We exploit this kinematic feature in final states including photons and leptons, tailored to the 100 TeV collider environment.

We examine models and prospects for proton annihilation to dileptons, a process which violates baryon and lepton number each by two. We determine that currently Super-Kamiokande would place the most draconian bound on $pp \rightarrow \ell^+ \ell^+$, ruling out new physics below a scale of $\sim 1.6$ TeV. We also find present and future hadron collider sensitivity to these processes. While 8 TeV LHC data excludes new physics at a scale below $\sim 800$ GeV, the reach of a 14 TeV LHC run is $\sim 1.8$ TeV, putting it on par with the sensitivity of Super-Kamiokande. On the other hand, a 100 TeV proton-proton collider would be sensitive to proton annihilation at a scale up to 10 TeV, allowing it to far exceed the reach of both Super-Kamiokande and the projected 2 TeV reach of Hyper-Kamiokande. Constraints from neutron star observation and cosmological evolution are not competitive. Therefore, although high-luminosity water Cherenkov experiments currently place the leading bounds on baryon and lepton number violation, next generation high-energy hadron colliders will begin surpassing them in sensitivity to some $B/L$-violating processes.

The paucity of old millisecond pulsars observed at the galactic center of the Milky Way could be the result of dark matter accumulating in and destroying neutron stars. In regions of high dark matter density, dark matter clumped in a pulsar can exceed the Schwarzschild limit and collapse into a natal black hole which destroys the pulsar. We examine what dark matter models are consistent with this hypothesis and find regions of parameter space where dark matter accumulation can significantly degrade the neutron star population within the galactic center while remaining consistent with observations of old millisecond pulsars in globular clusters and near the solar position. We identify what dark matter couplings and masses might cause a young pulsar at the galactic center to unexpectedly extinguish. Finally, we find that pulsar collapse age scales inversely with the dark matter density and linearly with the dark matter velocity dispersion. This implies that maximum pulsar age is spatially dependent on position within the dark matter halo of the Milky Way. In turn, this pulsar age spatial dependence will be dark matter model dependent.

By incorporating both the tensor-to-scalar ratio and the measured value of the spectral index, we set a bound on solo small field inflation of $\Delta \phi / m_{Pl} \geq 1.00 \sqrt{r/0.1}$. Unlike previous bounds which require monotonic $\epsilon_V,$ $|\eta_V| <1$, and 60 e-folds of inflation, the bound remains valid for non-monotonic $\epsilon_V$, $|\eta_V| \gtrsim 1$, and for inflation which occurs only over the 8 e-folds which have been observed on the cosmic microwave background. The negative value of the spectral index over the observed 8 e-folds is what makes the bound strong; we illustrate this by surveying single field models and finding that for $r \gtrsim 0.1$ and 8 e-folds of inflation, there is no simple potential which reproduces observed CMB perturbations and remains sub-Planckian. Models that are sub-Planckian after 8 e-folds must be patched together with a second epoch of inflation that fills out the remaining $\sim 50$ e-folds. This second, post-CMB epoch is characterized by extremely small $\epsilon_V$ and therefore an increasing scalar power spectrum. Using the fact that large power can overabundantly produce primordial black holes, we bound the maximum energy level of the second phase of inflation.

The existence of old neutron stars deeply constrains self-interacting fermion dark matter, which can form star-killing black holes. We quantify this constraint on dark matter-nucleon scattering, considering collapse scenarios that broaden bounds over intermediate masses. We then find the self- and co-annihilation rates necessary to lift these dark matter-nucleon scattering bounds. For Yukawa-coupled dark matter that fits dwarf galaxy halo profiles with a coupling $\alpha = 10^{-1}-10^{-4}$, a scalar mediator mass $m_\phi = 1-500$ MeV, and DM mass $m_X = 0.1-10^7$ GeV, we show that fermion dark matter is unconstrained if it self-annihilates at a rate greater than $10^{-40} ~ \rm{cm^3/s}$ or co-annihilates with baryons at a rate greater than $10^{-50} ~ \rm{cm^3/s}$.

If cosmological perturbations in our Hubble sized volume are nongaussian, then they will be coupled to any larger perturbation modes outside our Hubble volume. This has important consequences for modeling inflation: the scalar power spectrum and spectral index measured in our Hubble volume would depend on an adjacent background of super Hubble perturbations. In other words, a detection of nongaussianity implies a possible portal to the superhorizon. By the same token, ruling out nongaussianity would rule out the possibility that the power spectrum's size and running are accidents of super cosmic variance. In this note, we provide a compact derivation of super cosmic variance, survey recent results, and show how experimental limits on nongaussianity help to rule it out.

We demonstrate that local, scale-dependent non-Gaussianity can generate cosmic variance uncertainty in the observed spectral index of primordial curvature perturbations. In a universe much larger than our current Hubble volume, locally unobservable long wavelength modes can induce a scale-dependence in the power spectrum of typical subvolumes, so that the observed spectral index varies at a cosmologically significant level (|\Delta ns | ~ O(0.04)). Similarly, we show that the observed bispectrum can have an induced scale dependence that varies about the global shape. If tensor modes are coupled to long wavelength modes of a second field, the locally observed tensor power and spectral index can also vary. All of these effects, which can be introduced in models where the observed non-Gaussianity is consistent with bounds from the Planck satellite, loosen the constraints that observations place on the parameters of theories of inflation with mode coupling. We suggest observational constraints that future measurements could aim for to close this window of cosmic variance uncertainty.

We calculate non-Gaussianities in the bispectrum and trispectrum arising from the cubic term in the local expansion of the scalar curvature perturbation. We compute to three-loop order and for general momenta. A procedure for evaluating the leading behavior of the resulting loop-integrals is developed and discussed. Finally, we survey unique non-linear signals which could arise from the cubic term in the squeezed limit. In particular, it is shown that loop corrections can cause $f_{NL}^{sq.}$ to change sign as the momentum scale is varied. There also exists a momentum limit where $\tau_{NL} <0$ can be realized.