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.

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.

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.

Motivated by the recently reported diboson and dijet excesses in Run 1 data at ATLAS and CMS, we explore models of mixed dark matter in left-right symmetric theories. In this study, we calculate the relic abundance and the elastic scattering cross section with nuclei for a number of dark matter candidates that appear within the fermionic multiplets of left-right symmetric models. In contrast to the case of pure multiplets, WIMP-nucleon scattering proceeds at tree-level, and hence the projected reach of future direct detection experiments such as LUX-ZEPLIN and XENON1T will cover large regions of parameter space for TeV-scale thermal dark matter. Decays of the heavy charged W' boson to particles in the dark sector can potentially shift the right-handed gauge coupling to larger values when fixed to the rate of the Run 1 excesses, moving towards the theoretically attractive scenario, gR = gL. This region of parameter space may be probed by future collider searches for new Higgs bosons or electroweak fermions.

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.

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.

Simple models of weakly interacting massive particles (WIMPs) predict dark matter annihilations into pairs of electroweak gauge bosons, Higgses or tops, which through their subsequent cascade decays produce a spectrum of gamma rays. Intriguingly, an excess in gamma rays coming from near the Galactic center has been consistently observed in Fermi data. A recent analysis by the Fermi collaboration confirms these earlier results. Taking into account the systematic uncertainties in the modelling of the gamma ray backgrounds, we show for the first time that this excess can be well fit by these final states. In particular, for annihilations to $(WW, ZZ, hh, tt)$, dark matter with mass between threshold and approximately (165, 190, 280, 310) GeV gives an acceptable fit. The fit range for $b\bar{b}$ is also enlarged to 35 GeV$\lesssim m_\chi \lesssim$ 165 GeV. These are to be compared to previous fits that concluded only much lighter dark matter annihilating into $b$, $\tau$, and light quark final states could describe the excess. We demonstrate that simple, well-motivated models of WIMP dark matter including a thermal-relic neutralino of the MSSM, Higgs portal models, as well as other simplified models can explain the excess.

We analyze the recently released CoGeNT data with a focus on their time-dependent properties. Using various statistical techniques, we confirm the presence of modulation in the data, and find a significant component at high (E_{ee} > 1.5$ keVee) energies. We find that standard elastic WIMPs in a Maxwellian halo do not provide a good description of the modulation. We consider the possibility of non-standard halos, using halo independent techniques, and find a good agreement with the DAMA modulation for Q_{Na} \approx 0.3, but disfavoring interpretations with Q_{Na} = 0.5. The same techniques indicate that CDMS-Ge should see an O(1) modulation, and XENON100 should have seen 10-30 events (based upon the modulation in the 1.5-3.1 keVee range), unless L_{eff} is smaller than recent measurements. Models such as inelastic dark matter provide a good fit to the modulation, but not the spectrum. We note that tensions with XENON could be alleviated in such models if the peak is dominantly in April, when XENON data are not available due to noise.

Underground searches for dark matter involve a complicated interplay of particle physics, nuclear physics, atomic physics and astrophysics. We attempt to remove the uncertainties associated with astrophysics by developing the means to map the observed signal in one experiment directly into a predicted rate at another. We argue that it is possible to make experimental comparisons that are completely free of astrophysical uncertainties by focusing on {\em integral} quantities, such as $g(v_{min})=\int_{v_{min}} dv\, f(v)/v $ and $\int_{v_{thresh}} dv\, v g(v)$. Direct comparisons are possible when the $v_{min}$ space probed by different experiments overlap. As examples, we consider the possible dark matter signals at CoGeNT, DAMA and CRESST-Oxygen. We find that expected rate from CoGeNT in the XENON10 experiment is higher than observed, unless scintillation light output is low. Moreover, we determine that S2-only analyses are constraining, unless the charge yield $Q_y< 2.4 {\, \rm electrons/keV}$. For DAMA to be consistent with XENON10, we find for $q_{Na}=0.3$ that the modulation rate must be extremely high ($\gsim 70%$ for $m_\chi = 7\, \gev$), while for higher quenching factors, it makes an explicit prediction (0.8 - 0.9 cpd/kg) for the modulation to be observed at CoGeNT. Finally, we find CDMS-Si, even with a 10 keV threshold, as well as XENON10, even with low scintillation, would have seen significant rates if the excess events at CRESST arise from elastic WIMP scattering, making it very unlikely to be the explanation of this anomaly.

Dark matter pair production at high energy colliders may leave observable signatures in the energy and momentum spectra of the objects recoiling against the dark matter. We use LEP data on mono-photon events with large missing energy to constrain the coupling of dark matter to electrons. Within a large class of models, our limits are complementary to and competitive with limits on dark matter annihilation and on WIMP-nucleon scattering from indirect and direct searches. Our limits, however, do not suffer from systematic and astrophysical uncertainties associated with direct and indirect limits. For example, we are able to rule out light (< 10 GeV) thermal relic dark matter with universal couplings exclusively to charged leptons. In addition, for dark matter mass below about 80 GeV, LEP limits are stronger than Fermi constraints on annihilation into charged leptons in dwarf spheroidal galaxies. Within its kinematic reach, LEP also provides the strongest constraints on the spin-dependent direct detection cross section in models with universal couplings to both quarks and leptons. In such models the strongest limit is also set on spin independent scattering for dark matter masses below ~4 GeV. Throughout our discussion, we consider both low energy effective theories of dark matter, as well as several motivated renormalizable scenarios involving light mediators.

Direct detection of dark matter (DM) requires an interaction of dark matter particles with nucleons. The same interaction can lead to dark matter pair production at a hadron collider, and with the addition of initial state radiation this may lead to mono-jet signals. Mono-jet searches at the Tevatron can thus place limits on DM direct detection rates. We study these bounds both in the case where there is a contact interaction between DM and the standard model and where there is a mediator kinematically accessible at the Tevatron. We find that in many cases the Tevatron provides the current best limit, particularly for light dark matter, below 5 GeV, and for spin dependent interactions. Non-standard dark matter candidates are also constrained. The introduction of a light mediator significantly weakens the collider bound. A direct detection discovery that is in apparent conflict with mono-jet limits will thus point to a new light state coupling the standard model to the dark sector. Mono-jet searches with more luminosity and including the spectrum shape in the analysis can improve the constraints on DM-nucleon scattering cross section.

String and M-theory compactifications generically have compact moduli which can potentially act as the QCD axion. However, as demonstrated here, such a compact modulus can not play the role of a QCD axion and solve the strong CP problem if gravitational waves interpreted as arising from inflation with Hubble constant $H_inf \gsim 10^{13}$ GeV are observed by the PLANCK polarimetry experiment. In this case axion fluctuations generated during inflation would leave a measurable isocurvature and/or non-Gaussian imprint in the spectrum of primordial temperature fluctuations. This conclusion is independent of any assumptions about the initial axion misalignment angle, how much of the dark matter is relic axions, or possible entropy release by a late decaying particle such as the saxion; it relies only on the mild assumption that the Peccei-Quinn symmetry remains unbroken in the early universe.

The decay constant of the QCD axion is required by observation to be small compared to the Planck scale. In theories of "natural inflation," and certain proposed anthropic solutions of the cosmological constant problem, it would be interesting to obtain a large decay constant for axion-like fields from microscopic physics. String theory is the only context in which one can sensibly address this question. Here we survey a number of periodic fields in string theory in a variety of string vacua. In some examples, the decay constant can be parameterically larger than the Planck scale but the effective action then contains appreciable harmonics of order $f_A/M_p$. As a result, these fields are no better inflaton candidates than Planck scale axions.

We study the possibility of emission of the saxion, a superpartner of the axion, from SN1987A. The fact that the observed neutrino pulse from SN1987A is in excellent agreement with the current theory of supernovae places a strong bound on the energy loss into any non-standard model channel, therefore enabling bounds to be placed on the decay constant, f_a, of a light saxion. The low-energy coupling of the saxion, which couples at high energies to the QCD gauge field strength, is expected to be enhanced from QCD scaling, making it interesting to investigate if the saxion could place stronger bounds on f_a than the axion itself. Moreover, since the properties of the saxion are determined by f_a, a constraint on this parameter can be translated into a constraint on the supersymmetry breaking scale. We find that the bound on f_a from saxion emission is comparable with the one derived from axion emission due to a cancellation of leading-order terms in the soft-radiation expansion.

In an attempt to understand the properties of convective energy transport in the solar convection zone, a numerical model has been constructed for turbulent flows in a compressible, radiation-coupled, non-magnetic, gravitationally stratified medium using a realistic equation of state and realistic opacities. The time-dependent, three-dimensional hydrodynamic equations are solved with minimal simplifications. The statistical information obtained from the present simulation provides an improved understanding of solar photospheric convection. The characteristics of solar convection in shallow regions is parameterized and compared with the results of Chan and Sofia's simulations of deep and efficient convection (Chan and Sofia 1989). We assess the importance of the zones of partial ionization in the simulation, and confirm that the radiative energy transfer is negligible throughout the region except in the uppermost scale heights of the convection zone, a region of very high super-adiabaticity. When the effects of partial ionization are included, the dynamics of flows are altered significantly. However, we confirm the Chan and Sofia result that kinetic energy flux is non-negligible and can have a negative value in the convection zone.