We present the discovery of three sub-Neptune-sized planets transiting the nearby and bright Sun-like star HD 191939 (TIC 269701147, TOI 1339), a $K_{s}=7.18$ magnitude G8 V dwarf at a distance of only 54 parsecs. We validate the planetary nature of the transit signals by combining five months of data from the Transiting Exoplanet Survey Satellite with follow-up ground-based photometry, archival optical images, radial velocities, and high angular resolution observations. The three sub-Neptunes have similar radii ($R_{b} = 3.42^{+0.11}_{-0.11}\,R_{\oplus}$, $R_{c}=3.23_{-0.11}^{+0.11}\,R_{\oplus}$, and $R_{d}=3.16_{-0.11}^{+0.11}\,R_{\oplus}$) and their orbits are consistent with a stable, circular, and co-planar architecture near mean motion resonances of 1:3 and 3:4 ($P_{b}=8.88$ days, $P_{c}=28.58$ days, and $P_{d}=38.35$ days). The HD~191939 system is an excellent candidate for precise mass determinations of the planets with high-resolution spectroscopy due to the host star's brightness and low chromospheric activity. Moreover, the system's compact and near-resonant nature can provide an independent way to measure planetary masses via transit timing variations while also enabling dynamical and evolutionary studies. Finally, as a promising target for multi-wavelength transmission spectroscopy of all three planets' atmospheres, HD 191939 can offer valuable insight into multiple sub-Neptunes born from a proto-planetary disk that may have resembled that of the early Sun.

Although the $TESS$ primary mission observed the northern and southern ecliptic hemispheres, generally avoiding the ecliptic, and the $Kepler$ space telescope during the $K2$ mission could only observe near the ecliptic, many of the $K2$ fields extend far enough from the ecliptic plane that sections overlap with $TESS$ fields. Using photometric observations from both $K2$ and $TESS$, combined with archival spectroscopic observations, we globally modeled four known planetary systems discovered by $K2$ that were observed in the first year of the primary $TESS$ mission. Specifically, we provide updated ephemerides and system parameters for K2-114 b, K2-167 b, K2-237 b, and K2-261 b. These were some of the first $K2$ planets to be observed by $TESS$ in the first year and include three Jovian sized planets and a sub-Neptune with orbital periods less than 12 days. In each case, the updated ephemeris significantly reduces the uncertainty in prediction of future times of transit, which is valuable for planning observations with the James Webb Space Telescope and other future facilities. The $TESS$ extended mission is expected to observe about half of the $K2$ fields, providing the opportunity to perform this type of analysis on a larger number of systems.

We pioneer the black hole mass gap as a powerful new tool for constraining new particles. A new particle that couples to the Standard Model---such as an axion---acts as an additional source of loss in the cores of population-III stars, suppressing mass lost due to winds and quenching the pair-instability. This results in heavier astrophysical black holes. As an example, using stellar simulations we show that the solar axion explanation of the recent XENON1T excess implies astrophysical black holes of ~ 56 MS, squarely within the black hole mass gap predicted by the Standard Model.

New light particles produced in supernovae can lead to additional energy loss and a consequent deficit in neutrino production in conflict with the neutrinos observed from Supernova 1987A (SN1987A). Contrary to the majority of previous SN1987A studies, we examine the impact of $Z'$ bosons, axions, and axion-like particles (ALPs) interacting with the muons produced in SN1987A. For the first time, we find constraints on generic $Z'$ bosons coupled to muons, and apply our results to particle models including gauged $L_\mu-L_\tau$ number, $U(1)_{L_\mu-L_\tau}$, and gauged $B-L$ number, $U(1)_{B-L}$. We constrain $Z'$ bosons with masses up to about 250-500 MeV, and down to about $10^{-9}$ in $Z'$-muon coupling. We also extend previous work on axion-muon couplings by examining the importance of loop-level interactions, as well as performing calculations over a wider range of axion masses. We constrain muon-coupled axions from arbitrarily low masses up to about 200-500 MeV, with bounds extending down to axion-muon couplings of approximately $10^{-8}$ GeV$^{-1}$. We conclude that supernovae broadly provide a sensitive probe of new lightly-coupled particles interacting with muons.

Dark photon dark matter will resonantly convert into visible photons when the dark photon mass is equal to the plasma frequency of the ambient medium. In cosmological contexts, this transition leads to an extremely efficient, albeit short-lived, heating of the surrounding gas. Existing work in this field has been predominantly focused on understanding the implications of these resonant transitions in the limit that the plasma frequency of the Universe can be treated as being perfectly homogeneous, \ie neglecting inhomogeneities in the electron number density. In this work we focus on the implications of heating from dark photon dark matter in the presence of inhomogeneous structure (which is particularly relevant for dark photons with masses in the range $10^{-15} \, {\rm eV} \, \lesssim m_{A^\prime} \lesssim 10^{-12}$ eV), emphasizing both the importance of inhomogeneous energy injection, as well as the sensitivity of cosmological observations to the inhomogeneities themselves. More specifically, we derive modified constraints on dark photon dark matter from the Ly-$\alpha$ forest, and show that the presence of inhomogeneities allows one to extend constraints to masses outside of the range that would be obtainable in the homogeneous limit, while only slightly relaxing their strength. We then project sensitivity for near-future cosmological surveys that are hoping to measure the 21cm transition in neutral hydrogen prior to reionization, and demonstrate that these experiments will be extremely useful in improving sensitivity to masses near $\sim 10^{-14}$ eV, potentially by several orders of magnitude. Finally, we discuss implications for reionization, early star formation, and late-time $y$-type spectral distortions, and show that probes which are inherently sensitive to the inhomogeneous state of the Universe could resolve signatures unique to the light dark photon...

As space expands, the energy density in black holes increases relative to that of radiation, providing us with motivation to consider scenarios in which the early universe contained a significant abundance of such objects. In this study, we revisit the constraints on primordial black holes derived from measurements of the light element abundances. Black holes and their Hawking evaporation products can impact the era of Big Bang Nucleosynthesis (BBN) by altering the rate of expansion at the time of neutron-proton freeze-out, as well as by radiating mesons which can convert protons into neutrons and vice versa. Such black holes can thus enhance the primordial neutron-to-proton ratio, and increase the amount of helium that is ultimately produced. Additionally, the products of Hawking evaporation can break up helium nuclei, which both reduces the helium abundance and increases the abundance of primordial deuterium. Building upon previous work, we make use of modern deuterium and helium measurements to derive stringent constraints on black holes which evaporate in $t_{\rm evap} \sim 10^{-1}$ s to $\sim 10^{13}$ s (corresponding to $M \sim 6\times 10^8$ g to $\sim 2 \times 10^{13}$ g, assuming Standard Model particle content). We also consider how physics beyond the Standard Model could impact these constraints. Due to the gravitational nature of Hawking evaporation, the rate at which a black hole evaporates, and the types of particles that are produced through this process, depend on the complete particle spectrum. Within this context, we discuss scenarios which feature a large number of decoupled degrees-of-freedom (\ie~large hidden sectors), as well as models of TeV-scale supersymmetry.

Any abundance of black holes that was present in the early universe will evolve as matter, making up an increasingly large fraction of the total energy density as space expands. This motivates us to consider scenarios in which the early universe included an era that was dominated by low-mass ($M < 5\times 10^8$ g) black holes which evaporate prior to primordial nucleosynthesis. In significant regions of parameter space, these black holes will become gravitationally bound within binary systems, and undergo mergers before evaporating. Such mergers result in three potentially observable signatures. First, any black holes that have undergone one or more mergers will possess substantial angular momentum, causing their Hawking evaporation to produce significant quantities of high-energy gravitons. These products of Hawking evaporation are predicted to constitute a background of hot ($\sim$eV-keV) gravitons today, with an energy density corresponding to $\Delta N_{\rm eff} \sim 0.01-0.03$. Second, these mergers will produce a stochastic background of high-frequency gravitational waves. And third, the energy density of these gravitational waves can be as large as $\Delta N_{\rm eff} \sim 0.3$, depending on the length of time between the mergers and evaporation. These signals are each potentially within the reach of future measurements.

Light dark photons are subject to various plasma effects, such as Debye screening and resonant oscillations, which can lead to a more complex cosmological evolution than is experienced by conventional cold dark matter candidates. Maintaining a consistent history of dark photon dark matter requires ensuring that the super-thermal abundance present in the early Universe $\textit{(i)}$ does not deviate significantly after the formation of the CMB, and $\textit{(ii)}$ does not excessively leak into the Standard Model plasma after BBN. We point out that the role of non-resonant absorption, which has previously been neglected in cosmological studies of this dark matter candidate, produces strong constraints on dark photon dark matter with mass as low as $10^{-22}$ eV. Furthermore, we show that resonant conversion of dark photons after recombination can produce excessive heating of the IGM which is capable of prematurely reionizing hydrogen and helium, leaving a distinct imprint on both the Ly$-\alpha$ forest and the integrated optical depth of the CMB. Our constraints surpass existing cosmological bounds by more than five orders of magnitude across a wide range of dark photon masses.

The Fermi-LAT collaboration has recently released a new point source catalog, referred to as 4FGL. For the first time, we perform a template fit using information from this new catalog and find that the Galactic center excess is still present. On the other hand, we find that a wavelet-based search for point sources is highly sensitive to the use of the 4FGL catalog: no excess of bright regions on small angular scales is apparent when we mask out 4FGL point sources. We postulate that the 4FGL catalog contains the large majority of bright point sources that have previously been suggested to account for the excess in gamma rays detected at the Galactic center in Fermi-LAT data. Furthermore, after identifying which bright sources have no known counterpart, we place constraints on the luminosity function necessary for point sources to explain the smooth emission seen in the template fit.

Large, non-standard neutrino self-interactions have been shown to resolve the $\sim 4\sigma$ tension in Hubble constant measurements and a milder tension in the amplitude of matter fluctuations. We demonstrate that interactions of the necessary size imply the existence of a force-carrier with a large neutrino coupling ($> 10^{-4}$) and mass in the keV -- 100 MeV range. This mediator is subject to stringent cosmological and laboratory bounds, and we find that nearly all realizations of such a particle are excluded by existing data unless it carries spin 0 and couples almost exclusively to $\tau$-flavored neutrinos. Furthermore, we find that the light neutrinos must be Majorana, and that a UV-complete model requires a non-minimal mechanism to simultaneously generate neutrino masses and appreciable self-interactions.

We study the gravothermal evolution of dark matter halos in the presence of dissipative dark matter self-interactions. Dissipative interactions are present in many particle-physics realizations of the dark-sector paradigm and can significantly accelerate the gravothermal collapse of halos compared to purely elastic dark matter self-interactions. This is the case even when the dissipative interaction timescale is longer than the free-fall time of the halo. Using a semianalytical fluid model calibrated with isolated and cosmological $N$-body simulations, we calculate the evolution of the halo properties -- including its density profile and velocity dispersion profile -- as well as the core-collapse time as a function of the particle model parameters that describe the interactions. A key property is that the inner density profile at late times becomes cuspy again. Using 18 dwarf galaxies that exhibit a corelike dark matter density profile, we derive constraints on the strength of the dissipative interactions and the energy loss per collision.

We consider the Big Bang Nucleosynthesis (BBN) bounds on light dark matter whose cross section off nucleons is sufficiently large to enable acceleration by scattering off of cosmic rays in the local galaxy. Such accelerated DM could then deposit energy in terrestrial detectors. Since this signal involves DM of mass ~ keV - 100 MeV and requires large cross sections > 10^-31 cm^2 in a relativistic kinematic regime, we find that the DM population in this scenario is generically equilibrated with Standard Model particles in the early universe. For sufficiently low DM masses < 10 MeV, corresponding to the bulk of the favored region of many cosmic-ray upscattering studies, this equilibrated DM population adds an additional component to the relativistic energy density around T ~ few MeV and thereby spoils the successful predictions of BBN. In the remaining ~ 10-100 MeV mass range, the large couplings required in this scenario are either currently excluded or within reach of current or future accelerator-based searches.

We report the detection of a transiting Earth-size planet around GJ 357, a nearby M2.5V star, using data from the Transiting Exoplanet Survey Satellite (TESS). GJ 357 b (TOI-562.01) is a transiting, hot, Earth-sized planet (Teq=525+-11 K) with a radius of Rb=1.217+-0.084 Re and an orbital period of Pb=3.93 d. Precise stellar radial velocities from CARMENES and PFS, as well as archival data from HIRES, UVES, and HARPS also display a 3.93-day periodicity, confirming the planetary nature and leading to a planetary mass of Mb=1.84+-0.31 Me. In addition to the radial velocity signal for GJ 357 b, more periodicities are present in the data indicating the presence of two further planets in the system: GJ 357 c, with a minimum mass of Mc=3.40+-0.46 Me in a 9.12 d orbit, and GJ 357 d, with a minimum mass of Md=6.1+-1.0 Me in a 55.7 d orbit inside the habitable zone. The host is relatively inactive and exhibits a photometric rotation period of Prot=78+-2 d. GJ 357 b is to date the second closest transiting planet to the Sun, making it a prime target for further investigations such as transmission spectroscopy. Therefore, GJ 357 b represents one of the best terrestrial planets suitable for atmospheric characterization with the upcoming JWST and ground-based ELTs.

We revisit constraints on dark photons with masses below ~ 100 MeV from the observations of Supernova 1987A. If dark photons are produced in sufficient quantity, they reduce the amount of energy emitted in the form of neutrinos, in conflict with observations. For the first time, we include the effects of finite temperature and density on the kinetic-mixing parameter, epsilon, in this environment. This causes the constraints on epsilon to weaken with the dark-photon mass below ~ 15 MeV. For large-enough values of epsilon, it is well known that dark photons can be reabsorbed within the supernova. Since the rates of reabsorption processes decrease as the dark-photon energy increases, we point out that dark photons with energies above the Wien peak can escape without scattering, contributing more to energy loss than is possible assuming a blackbody spectrum. Furthermore, we estimate the systematic uncertainties on the cooling bounds by deriving constraints assuming one analytic and four different simulated temperature and density profiles of the proto-neutron star. Finally, we estimate also the systematic uncertainty on the bound by varying the distance across which dark photons must propagate from their point of production to be able to affect the star. This work clarifies the bounds from SN1987A on the dark-photon parameter space.

If even a relatively small number of black holes were created in the early universe, they will constitute an increasingly large fraction of the total energy density as space expands. It is thus well-motivated to consider scenarios in which the early universe included an era in which primordial black holes dominated the total energy density. Within this context, we consider Hawking radiation as a mechanism to produce both dark radiation and dark matter. If the early universe included a black hole dominated era, we find that Hawking radiation will produce dark radiation at a level $\Delta N_{\rm eff} \sim 0.03-0.2$ for each light and decoupled species of spin 0, 1/2, or 1. This range is well suited to relax the tension between late and early-time Hubble determinations, and is within the reach of upcoming CMB experiments. The dark matter could also originate as Hawking radiation in a black hole dominated early universe, although such dark matter candidates must be very heavy ($m_{\rm DM} >10^{11}$ GeV) if they are to avoid exceeding the measured abundance.

Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ .

The future of exoplanet science is bright, as TESS once again demonstrates with the discovery of its longest-period confirmed planet to date. We hereby present HD 21749b (TOI 186.01), a sub-Neptune in a 36-day orbit around a bright (V = 8.1) nearby (16 pc) K4.5 dwarf. TESS measures HD21749b to be 2.61$^{+0.17}_{-0.16}$ $R_{\oplus}$, and combined archival and follow-up precision radial velocity data put the mass of the planet at $22.7^{+2.2}_{-1.9}$ $M_{\oplus}$. HD 21749b contributes to the TESS Level 1 Science Requirement of providing 50 transiting planets smaller than 4 $R_{\oplus}$ with measured masses. Furthermore, we report the discovery of HD 21749c (TOI 186.02), the first Earth-sized ($R_p = 0.892^{+0.064}_{-0.058} R_{\oplus}$) planet from TESS. The HD21749 system is a prime target for comparative studies of planetary composition and architecture in multi-planet systems.

Cosmological observations offer unique and robust avenues for probing the fundamental nature of dark matter particles-they broadly test a range of compelling theoretical scenarios, often surpassing or complementing the reach of terrestrial and other experiments. We discuss observational and theoretical advancements that will play a pivotal role in realizing a strong program of cosmological searches for the identity of dark matter in the coming decade. Specifically, we focus on measurements of the cosmic-microwave-background anisotropy and spectral distortions, and tracers of structure (such as the Lyman-$\alpha$ forest, galaxies, and the cosmological 21-cm signal).

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

We study the effect of a dibaryon, S, in the mass range 1860 MeV < m_S < 2054 MeV, which is heavy enough not to disturb the stability of nuclei and light enough to possibly be cosmologically metastable. Such a deeply bound state can act as a baryon sink in regions of high baryon density and temperature. We find that the ambient conditions encountered inside a newly born neutron star are likely to sustain a sufficient population of hyperons to ensure that a population of S dibaryons can equilibrate in less than a few seconds. This would be catastrophic for the stability of neutron stars and the observation of neutrino emission from the proto-neutron star of Supernova 1987A over ~ O(10)s. A deeply bound dibaryon is therefore incompatible with the observed supernova explosion, unless the cross section for S production is severely suppressed.

Sasankan et al, have recently claimed that there are significant deviations in the phase-space distributions of the kinetic energies of nuclei from the Maxwell-Boltzmann form usually assumed in BBN, and further, that these deviations lead to big changes in the predicted light-element abundances. Solving the relativistic Boltzmann equation perturbatively, we explicitly show that these deviations are not 20% as claimed, but rather are about 10^-17 in size and hence cannot significantly alter the predicted light-element abundances. We discuss and compute two related effects: O(0.1%) corrections to the kinetic distributions of nuclei that arise from small relativistic corrections to the MB distribution and a much smaller effect, nuclear kinetic drag, which arises from the heat transferred from the EM plasma to nuclei that is needed to maintain kinetic equilibrium.

We consider the constraints from Supernova 1987A on particles with small couplings to the Standard Model. We discuss a model with a fermion coupled to a dark photon, with various mass relations in the dark sector; millicharged particles; dark-sector fermions with inelastic transitions; the hadronic QCD axion; and an axion-like particle that couples to Standard Model fermions with couplings proportional to their mass. In the fermion cases, we develop a new diagnostic for assessing when such a particle is trapped at large mixing angles. Our bounds for a fermion coupled to a dark photon constrain small couplings and masses <200 MeV, and do not decouple for low fermion masses. They exclude parameter space that is otherwise unconstrained by existing accelerator-based and direct-detection searches. In addition, our bounds are complementary to proposed laboratory searches for sub-GeV dark matter, and do not constrain several "thermal" benchmark-model targets. For a millicharged particle, we exclude charges between 10^(-9) to a few times 10^(-6) in units of the electron charge; this excludes parameter space to higher millicharges and masses than previous bounds. For the QCD axion and an axion-like particle, we apply several updated nuclear physics calculations and include the energy dependence of the optical depth to accurately account for energy loss at large couplings. We rule out a hadronic axion of mass between 0.1 and a few hundred eV, or equivalently bound the PQ scale between a few times 10^4 and 10^8 GeV, closing the hadronic axion window. For an axion-like particle, our bounds disfavor decay constants between a few times 10^5 GeV up to a few times 10^8 GeV. In all cases, our bounds differ from previous work by more than an order of magnitude across the entire parameter space. We also provide estimated systematic errors due to the uncertainties of the progenitor.

We propose a class of models in which a stable inflaton is produced as a thermal relic in the early universe and constitutes the dark matter. We show that inflaton annihilations can efficiently reheat the universe, and identify several examples of inflationary potentials that can accommodate all cosmic microwave background observables and in which the inflaton dark matter candidate has a weak scale mass. As a simple example, we consider annihilations that take place through a Higgs portal interaction, leading to encouraging prospects for future direct detection experiments.

Due to shielding, direct detection experiments are in some cases insensitive to dark matter candidates with very large scattering cross sections with nucleons. In this paper, we revisit this class of models, and derive a simple analytic criterion for conservative but robust direct detection limits. While large spin-independent cross sections seem to be ruled out, we identify potentially viable parameter space for dark matter with a spin-dependent cross section with nucleons in the range of $10^{-27} {\rm cm}^2 < \sigma_{{\rm DM}-p} < 10^{-24} \, {\rm cm}^{2}$. With these parameters, cosmic-ray scattering with dark matter in the extended halo of the Milky Way could generate a novel and distinctive gamma-ray signal at high galactic latitudes. Such a signal could be observable by Fermi or future space-based gamma-ray telescopes.

The EDGES Collaboration has recently reported the detection of a stronger-than-expected absorption feature in the global 21-cm spectrum, centered at a frequency corresponding to a redshift of z ~ 17. This observation has been interpreted as evidence that the gas was cooled during this era as a result of scattering with dark matter. In this study, we explore this possibility, applying constraints from the cosmic microwave background, light element abundances, Supernova 1987A, and a variety of laboratory experiments. After taking these constraints into account, we find that the vast majority of the parameter space capable of generating the observed 21-cm signal is ruled out. The only range of models that remains viable is that in which a small fraction, ~ 0.3-2%, of the dark matter consists of particles with a mass of ~ 10-80 MeV and which couple to the photon through a small electric charge, epsilon ~ 10^{-6}-10^{-4}. Furthermore, in order to avoid being overproduced in the early universe, such models must be supplemented with an additional depletion mechanism, such as annihilations through a L_{\mu}-L_{\tau} gauge boson or annihilations to a pair of rapidly decaying hidden sector scalars.

We analyze the gamma-ray sky at energies of 0.5 to 50 GeV using the undecimated wavelet transform on the sphere. Focusing on the inner $60^{\circ} \times 60^{\circ}$ of the sky, we identify and characterize four separate residuals beyond the expected Milky Way diffuse emission. We detect the \textit{Fermi} Bubbles, finding compelling evidence that they are diffuse in nature and contain very little small-scale structure. We detect the "cocoon" inside the Southern Bubble, and we also identify its northern counterpart above 2 GeV. The Northern Cocoon lies along the same axis but is $\sim 30 \%$ dimmer than the southern one. We characterize the Galactic center excess, which we find extends up to $20^{\circ}$ in $|b|$. At latitudes $|b| \leq 5^{\circ}$ we find evidence for power in small angular scales that could be the result of point-source contributions, but for $|b| \geq 5^{\circ}$ the Galactic center excess is dominantly diffuse in its nature. Our findings show that either the Galactic center excess and {\it Fermi} Bubbles connect smoothly or that the Bubbles brighten significantly below $15^\circ$ in latitude. We find that the Galactic center excess appears off-center by a few degrees towards negative $\ell$. Additionally, we find and characterize two emissions along the Galactic disk centered at $\ell \simeq +25^{\circ}$ and $-20^{\circ}$. These emissions are significantly more elongated along the Galactic disk than the Galactic center excess.

We suggest that two-to-two dark matter fusion may be the relaxation process that resolves the small-scale structure problems of the cold collisionless dark matter paradigm. In order for the fusion cross section to scale correctly across many decades of astrophysical masses from dwarf galaxies to galaxy clusters, we require the fractional binding energy released to be greater than v^n ~ [10^{-(2-3)}]^n, where n=1,2 depends on local dark sector chemistry. The size of the dark-sector interaction cross sections must be sigma ~ 0.1-1 barn, moderately larger than for Standard Model deuteron fusion, indicating a dark nuclear scale Lambda ~ O(100 MeV). Dark fusion firmly predicts constant sigma v below the characteristic velocities of galaxy clusters. Observations of the inner structure of galaxy groups with velocity dispersion of several hundred kilometer per second, of which a handful have been identified, could differentiate dark fusion from a dark photon model.

This white paper summarizes the workshop "U.S. Cosmic Visions: New Ideas in Dark Matter" held at University of Maryland on March 23-25, 2017.

We calculate the exact width for a dark photon decaying to three photons at one loop order for dark photon masses m' below the e+e- production threshold of 2m_e. We find substantial deviations from previous results derived from the lowest order Euler-Heisenberg effective Lagrangian in the range m_e < m' < 2m_e, where higher order terms in the derivative expansion are nonnegligible. This mass range is precisely where the three photon decay takes place on cosmologically relevant timescales. Our improved analysis opens a window for dark photons in the range 850 keV < m' < 2m_e, 10^-5 < epsilon < 10^-4.

It has recently been demonstrated that, in the event of a putative signal in dark matter direct detection experiments, properly identifying the underlying dark matter-nuclei interaction promises to be a challenging task. Given the most optimistic expectations for the number counts of recoil events in the forthcoming Generation 2 experiments, differentiating between interactions that produce distinct features in the recoil energy spectra will only be possible if a strong signal is observed simultaneously on a variety of complementary targets. However, there is a wide range of viable theories that give rise to virtually identical energy spectra, and may only differ by the dependence of the recoil rate on the dark matter velocity. In this work, we investigate how degeneracy between such competing models may be broken by analyzing the time dependence of nuclear recoils, i.e. the annual modulation of the rate. For this purpose, we simulate dark matter events for a variety of interactions and experiments, and perform a Bayesian model-selection analysis on all simulated data sets, evaluating the chance of correctly identifying the input model for a given experimental setup. We find that including information on the annual modulation of the rate may significantly enhance the ability of a single target to distinguish dark matter models with nearly degenerate recoil spectra, but only with exposures beyond the expectations of Generation 2 experiments.

We demonstrate how the image analysis technique of wavelet decomposition can be applied to the gamma-ray sky to separate emission on different angular scales. New structures on scales that differ from the scales of the conventional astrophysical foreground and background uncertainties can be robustly extracted, allowing a model-independent characterization with no presumption of exact signal morphology. As a test case, we generate mock gamma-ray data to demonstrate our ability to extract extended signals without assuming a fixed spatial template. For some point source luminosity functions, our technique also allows us to differentiate a diffuse signal in gamma-rays from dark matter annihilation and extended gamma-ray point source populations in a data-driven way.

Identifying the true theory of dark matter depends crucially on accurately characterizing interactions of dark matter (DM) with other species. In the context of DM direct detection, we present a study of the prospects for correctly identifying the low-energy effective DM-nucleus scattering operators connected to UV-complete models of DM-quark interactions. We take a census of plausible UV-complete interaction models with different low-energy leading-order DM-nuclear responses. For each model (corresponding to different spin-, momentum-, and velocity-dependent responses), we create a large number of realizations of recoil-energy spectra, and use Bayesian methods to investigate the probability that experiments will be able to select the correct scattering model within a broad set of competing scattering hypotheses. We conclude that agnostic analysis of a strong signal (such as Generation-2 would see if cross sections are just below the current limits) seen on xenon and germanium experiments is likely to correctly identify momentum dependence of the dominant response, ruling out models with either "heavy" or "light" mediators, and enabling downselection of allowed models. However, a unique determination of the correct UV completion will critically depend on the availability of measurements from a wider variety of nuclear targets, including iodine or fluorine. We investigate how model-selection prospects depend on the energy window available for the analysis. In addition, we discuss accuracy of the DM particle mass determination under a wide variety of scattering models, and investigate impact of the specific types of particle-physics uncertainties on prospects for model selection.

We describe a general scenario, dubbed "Inflatable Dark Matter", in which the density of dark matter particles can be reduced through a short period of late-time inflation in the early universe. The overproduction of dark matter that is predicted within many otherwise well-motivated models of new physics can be elegantly remedied within this context, without the need to tune underlying parameters or to appeal to anthropic considerations. Thermal relics that would otherwise be disfavored can easily be accommodated within this class of scenarios, including dark matter candidates that are very heavy or very light. Furthermore, the non-thermal abundance of GUT or Planck scale axions can be brought to acceptable levels, without invoking anthropic tuning of initial conditions. A period of late-time inflation could have occurred over a wide range of scales from ~ MeV to the weak scale or above, and could have been triggered by physics within a hidden sector, with small but not necessarily negligible couplings to the Standard Model.

Dark matter particles annihilating into Standard Model fermions may be able to explain the recent observation of a gamma-ray excess in the direction of the Galactic Center. Recently, a hidden photon model has been proposed to explain this signal. Supplementing this model with a dipole moment operator and a small dark sector mass splitting allows a large cross section to a photon line while avoiding direct detection and other constraints. Comparing the line and continuum cross sections, we find that the line is suppressed only by the relative scales and couplings. Given current constraints on this ratio, a line discovery in the near future could point to a new scale Lambda ~ O(1 TeV), where we would expect to discover new charged particles. Moreover, such a line would also imply that dark matter can be visible in near-future direct detection experiments.

Motivated by the gamma-ray excess observed from the region surrounding the Galactic Center, we explore particle dark matter models that could potentially account for the spectrum and normalization of this signal. Taking a model-independent approach, we consider an exhaustive list of tree-level diagrams for dark matter annihilation, and determine which could account for the observed gamma-ray emission while simultaneously predicting a thermal relic abundance equal to the measured cosmological dark matter density. We identify a wide variety of models that can meet these criteria without conflicting with existing constraints from direct detection experiments or the Large Hadron Collider (LHC). The prospects for detection in near future dark matter experiments and/or the upcoming 14 TeV LHC appear quite promising.

The gamma-ray excess observed from the Galactic Center can be interpreted as dark matter particles annihilating into Standard Model fermions with a cross section near that expected for a thermal relic. Although many particle physics models have been shown to be able to account for this signal, the fact that this particle has not yet been observed in direct detection experiments somewhat restricts the nature of its interactions. One way to suppress the dark matter's elastic scattering cross section with nuclei is to consider models in which the dark matter is part of a hidden sector. In such models, the dark matter can annihilate into other hidden sector particles, which then decay into Standard Model fermions through a small degree of mixing with the photon, Z, or Higgs bosons. After discussing the gamma-ray signal from hidden sector dark matter in general terms, we consider two concrete realizations: a hidden photon model in which the dark matter annihilates into a pair of vector gauge bosons that decay through kinetic mixing with the photon, and a scenario within the generalized NMSSM in which the dark matter is a singlino-like neutralino that annihilates into a pair of singlet Higgs bosons, which decay through their mixing with the Higgs bosons of the MSSM.

We present constraints on decaying and annihilating dark matter (DM) in the 4 keV to 10 GeV mass range, using published results from the satellites HEAO-1, INTEGRAL, COMPTEL, EGRET, and the Fermi Gamma-ray Space Telescope. We derive analytic expressions for the gamma-ray spectra from various DM decay modes, and find lifetime constraints in the range 10^24-10^28 sec, depending on the DM mass and decay mode. We map these constraints onto the parameter space for a variety of models, including a hidden photino that is part of a kinetically mixed hidden sector, a gravitino with R-parity violating decays, a sterile neutrino, DM with a dipole moment, and a dark pion. The indirect constraints on sterile-neutrino and hidden-photino DM are found to be more powerful than other experimental or astrophysical probes in some parts of parameter space. While our focus is on decaying DM, we also present constraints on DM annihilation to electron-positron pairs. We find that if the annihilation is p-wave suppressed, the galactic diffuse constraints are, depending on the DM mass and velocity at recombination, more powerful than the constraints from the Cosmic Microwave Background.

Several classes of astrophysical sources contribute to the approximately isotropic gamma-ray background measured by the Fermi Gamma-Ray Space Telescope. In this paper, we use Fermi's catalog of gamma-ray sources (along with corresponding source catalogs at infrared and radio wavelengths) to build and constrain a model for the contributions to the extragalactic gamma-ray background from astrophysical sources, including radio galaxies, star-forming galaxies, and blazars. We then combine our model with Fermi's measurement of the gamma-ray background to derive constraints on the dark matter annihilation cross section, including contributions from both extragalactic and galactic halos and subhalos. The resulting constraints are competitive with the strongest current constraints from the Galactic Center and dwarf spheroidal galaxies. As Fermi continues to measure the gamma-ray emission from a greater number of astrophysical sources, it will become possible to more tightly constrain the astrophysical contributions to the extragalactic gamma-ray background. We project that with 10 years of data, Fermi's measurement of this background combined with the improved constraints on the astrophysical source contributions will yield a sensitivity to dark matter annihilations that exceeds the strongest current constraints by a factor of ~ 5 - 10.

Evidence continues to grow in the MiniBooNE (MB) antineutrino mode supporting a low-energy excess compatible with the MB neutrino mode and possibly also confirming the results of the LSND experiment. At least one sterile neutrino is required to explain the anomalies consistent with the observations of other experiments. At the same time, there is a strong tension between the positive signals of LSND and MB and the null results of nu_e and nu_mu disappearance experiments. We explore a scenario, first proposed in \cite{Nelson:2010hz}, where the presence of an additional heavy sterile neutrino (with mass well above an eV) can alleviate tension between LSND, MB and the null results of disappearance experiments. We compare and contrast this 3+1+1 scenario with the more standard 3+1 scenario and carry out global fits to all oscillation data including new 2011 MB anti-nu data. We find that the tension can be somewhat alleviated and that a phenomenologically viable window for the heavy neutrino, consistent with rare decays and BBN constraints, can be found if the fifth neutrino has a mass of order 0.3 - 10 GeV. We also find, however, that the 2011 MB anti-nu data exacerbates the tension with null experiments in both the 3+1 and 3+1+1 models when the lowest energy bins are included, resulting in little improvement in the global fit. We also discuss the implications of an additional neutrino for the reactor and gallium anomalies, and show that an oscillation explanation of the anomalies is disfavored by cosmological considerations, direct searches, and precision electroweak tests.

We consider possibly observable effects of asymmetric dark matter (ADM) in neutron stars. Since dark matter does not self-annihilate in the ADM scenario, dark matter accumulates in neutron stars, eventually reaching the Chandrasekhar limit and forming a black hole. We focus on the case of scalar ADM, where the constraints from Bose-Einstein condensation and subsequent black hole formation are most severe due to the absence of Fermi degeneracy pressure. We also note that in some portions of this constrained parameter space, non-trivial effects from Hawking radiation can modify our limits. We find that for scalar ADM with mass between 100 keV and 10^5 GeV, the constraint from pulsars in globular clusters on the scattering cross-section with neutrons ranges from \sigma_n < 10^{-45} cm^2 to 10^{-52} cm}^2. In particular, for scalar ADM with mass between 1 GeV and 1 TeV (in the case where black hole evaporation due to Hawking radiation is unimportant), the constraint on the scattering cross-section is below what is reachable with ton scale direct detection experiments.

We consider current observational constraints on the electromagnetic charge of dark matter. The velocity dependence of the scattering cross-section through the photon gives rise to qualitatively different constraints than standard dark matter scattering through massive force carriers. In particular, recombination epoch observations of dark matter density perturbations require that $\epsilon$, the ratio of the dark matter to electronic charge, is less than $10^{-6}$ for $m_X = 1 GeV$, rising to $\epsilon < 10^{-4}$ for $m_X = 10 TeV$. Though naively one would expect that dark matter carrying a charge well below this constraint could still give rise to large scattering in current direct detection experiments, we show that charged dark matter particles that could be detected with upcoming experiments are expected to be evacuated from the Galactic disk by the Galactic magnetic fields and supernova shock waves, and hence will not give rise to a signal. Thus dark matter with a small charge is likely not a source of a signal in current or upcoming dark matter direct detection experiments.