We discuss a minimal realization of the strongly interacting massive particle (SIMP) framework. The model includes a dark copy of QCD with three colors and three light flavors. A massive dark photon, kinetically mixed with the Standard Model hypercharge, maintains kinetic equilibrium between the dark and visible sectors. One of the dark mesons is necessarily unstable but long-lived, with potential impact on CMB observables. We show that an approximate "isospin" symmetry acting on the down-type quarks is an essential ingredient of the model. This symmetry stabilizes the dark matter and allows to split sufficiently the masses of the other states to suppress strongly their relic abundances. We discuss for the first time the SIMP cosmology with sizable mass splittings between all meson multiplets. We demonstrate that the SIMP mechanism remains efficient in setting the dark matter relic density, while CMB constraints on unstable relics can be robustly avoided. We also consider the phenomenological consequences of isospin breaking, including dark matter decay. Cosmological, astrophysical, and terrestrial probes are combined into a global picture of the parameter space. In addition, we outline an ultraviolet completion in the context of neutral naturalness, where confinement at the GeV scale is generic. We emphasize the general applicability of several novel features of the SIMP mechanism that we discuss here.

A hidden sector that kinetically mixes with the Minimal Supersymmetric Standard Model provides simple and well-motivated dark matter candidates that possess many of the properties of a traditional weakly interacting massive particle (WIMP). These supersymmetric constructions can also provide a natural explanation for why the dark matter is at the weak scale - even if it resides in a hidden sector. In the hidden sector, a natural pattern of symmetry breaking generally makes particles and their superpartners lie around the same mass scale, opening novel possibilities for a variety of cosmological histories and complex indirect detection signatures.

We study scenarios where loop processes give the dominant contributions to dark matter decay or annihilation despite the presence of tree level channels. We illustrate this possibility in a specific model where dark matter is part of a hidden sector that communicates with the Standard Model sector via a heavy neutrino portal. We explain the underpinning rationale for how loop processes mediated by the portal neutrinos can parametrically dominate over tree level decay channels, and demonstrate that this qualitatively changes the indirect detection signals in positrons, neutrinos, and gamma rays.

We study the phenomenology of a light (GeV scale) sterile neutrino sector and the pseudo-Goldstone boson (not the majoron) associated with a global symmetry in this sector that is broken at a high scale. Such scenarios can be motivated from considerations of singlet fermions from a hidden sector coupling to active neutrinos via heavy right-handed seesaw neutrinos, effectively giving rise to a secondary, low-energy seesaw framework. Such scenarios allow for rich phenomenology with observable implications for cosmology, dark matter, and direct searches, involving novel sterile neutrino dark matter production mechanisms from the pseudo-Goldstone-mediated scattering or decay, modifications of BBN bounds on sterile neutrinos, suppression of canonical sterile neutrino decay channels at direct search experiments, late injection of an additional population of neutrinos in the Universe after neutrino decoupling, and measurable dark radiation at BBN or CMB decoupling.

Imaging Atmospheric Cherenkov Telescope (IACT) searches for dark matter often perform observations in "wobble mode", i.e. collecting data from the signal region and from a corresponding background control region at the same time, enabling efficient simultaneous determination and subtraction of background. This observation strategy is possibly compromised in scenarios where dark matter annihilates to long-lived mediators that can traverse astrophysical distances before decaying to produce the gamma rays observed by the IACTs. In this paper, we show that this challenge comes with several interesting features and opportunities: in addition to signal contamination in the background control region, the gamma-ray spectrum changes with the observing direction angle and typically exhibits a hard excess at high energies. This affects signal reconstruction via subtraction of the background control region measurements in non-trivial ways. Such features represent a significant departure from the canonical picture, and offer novel handles to identify a dark matter signal and to extract underlying dark matter parameters.

We investigate observable cosmological aspects of sterile neutrino dark matter produced via the freeze-in mechanism. The study is performed in a framework that admits many cosmologically interesting variations: high temperature production via annihilation processes from higher dimensional operators or low temperature production from decays of a scalar, with the decaying scalar in or out of equilibrium with the thermal bath, in supersymmetric or non-supersymmetric setups, thus allowing us to both extract generic properties and highlight features unique to particular variations. We find that while such sterile neutrinos are generally compatible with all cosmological constraints, interesting scenarios can arise where dark matter is cold, warm, or hot, has nontrivial momentum distributions, or provides contributions to the effective number of relativistic degrees of freedom Neff during Big Bang nucleosynthesis large enough to be probed by future measurements.

A remarkable prediction of the Standard Model is that, in the absence of corrections lifting the energy density, the Higgs potential becomes negative at large field values. If the Higgs field samples this part of the potential during inflation, the negative energy density may locally destabilize the spacetime. We use numerical simulations of the Einstein equations to study the evolution of inflation-induced Higgs fluctuations as they grow towards the true (negative-energy) minimum. These simulations show that forming a single patch of true vacuum in our past light cone during inflation is incompatible with the existence of our Universe; the boundary of the true vacuum region grows outward in a causally disconnected manner from the crunching interior, which forms a black hole. We also find that these black hole horizons may be arbitrarily elongated---even forming black strings---in violation of the hoop conjecture. By extending the numerical solution of the Fokker-Planck equation to the exponentially suppressed tails of the field distribution at large field values, we derive a rigorous correlation between a future measurement of the tensor-to-scalar ratio and the scale at which the Higgs potential must receive stabilizing corrections in order for the Universe to have survived inflation until today.

Sterile neutrino dark matter, a popular alternative to the WIMP paradigm, has generally been studied in non-supersymmetric setups. If the underlying theory is supersymmetric, we find that several interesting and novel dark matter features can arise. In particular, in scenarios of freeze-in production of sterile neutrino dark matter, its superpartner, the sterile sneutrino, can play a crucial role in early Universe cosmology as the dominant source of cold, warm, or hot dark matter, or of a subdominant relativistic population of sterile neutrinos that can contribute to the effective number of relativistic degrees of freedom Neff during Big Bang nucleosynthesis.

The pseudoscalar resonance or "A-funnel" in the Minimal Supersymmetric Standard Model~(MSSM) is a widely studied framework for explaining dark matter that can yield interesting indirect detection and collider signals. The well-known Galactic Center excess (GCE) at GeV energies in the gamma ray spectrum, consistent with annihilation of a $\lesssim 40$ GeV dark matter particle, has more recently been shown to be compatible with significantly heavier masses following reanalysis of the background. In this paper, we explore the LHC and direct detection implications of interpreting the GCE in this extended mass window within the MSSM A-funnel framework. We find that compatibility with relic density, signal strength, collider constraints, and Higgs data can be simultaneously achieved with appropriate parameter choices. The compatible regions give very sharp predictions of 200-600 GeV CP-odd/even Higgs bosons at low tan$\beta$ at the LHC and spin-independent cross sections $\approx 10^{-11}$ pb at direct detection experiments. Regardless of consistency with the GCE, this study serves as a useful template of the strong correlations between indirect, direct, and LHC signatures of the MSSM A-funnel region.

A sterile neutrino is a well-motivated and widely studied dark matter candidate. The most straightforward realization of sterile neutrino dark matter, through the Dodelson-Widrow mechanism, is now ruled out by a combination of X-ray and Lyman-\alpha measurements. An alternative production mechanism that is becoming increasingly popular in the literature is the freeze-in mechanism, involving frameworks where a feeble coupling to a particle - usually a scalar beyond the Standard Model - in the thermal bath results in a gradual accumulation of the sterile neutrino dark matter abundance. This article reviews the various motivations for realizing such frameworks in the literature, their common characteristic features, and phenomenological signatures.

Neutrino masses and light (keV-GeV) sterile neutrinos can arise naturally via a modified, low energy seesaw mechanism if the right-handed neutrinos are charged under a new symmetry broken by a PeV scale vacuum expectation value, presumably tied to supersymmetry breaking. The additional field content also allows for freeze-in production of sterile neutrino dark matter. This framework can accommodate the recently observed 3.5 keV X-ray line, while a straightforward extension of the framework, using the new symmetry and the PeV energy scale, can explain the PeV energy neutrino events at IceCube. Together, these can therefore be taken as hints of the existence of a PeV scale supersymmetric neutrino sector.

We show that active neutrino masses and a keV-GeV mass sterile neutrino dark matter candidate can result from a modified, low energy seesaw mechanism if right-handed neutrinos are charged under a new symmetry broken by a scalar field vacuum expectation value at the PeV scale. The dark matter relic abundance can be obtained through active-sterile oscillation, freeze-in through the decay of the heavy scalar, or freeze-in via non-renormalizable interactions at high temperatures. The low energy effective theory maps onto the widely studied nuMSM framework.

Recent measurements of PeV energy neutrinos at IceCube and a 3.5 keV X-ray line in the spectra of several galaxies are both tantalizing signatures of new physics. This paper shows that one or both of these observations can be explained within an extended supersymmetric neutrino sector. Obtaining light active neutrino masses as well as phenomenologically interesting (keV-GeV) sterile neutrino masses without any unnaturally small parameters hints at a new symmetry in the neutrino sector that is broken at the PeV scale, presumably tied to supersymmetry breaking. The same symmetry and structure can sufficiently stabilize an additional PeV particle, produce its abundance through the freeze-in mechanism, and lead to decays that can give the energetic neutrinos observed by IceCube. The lightest sterile neutrino, if at 7 keV, is a non-resonantly produced fraction of dark matter, and can account for the 3.5 keV X-ray line. The two signals could therefore be the first probes of an extended supersymmetric neutrino sector.

Measurements of the Higgs boson and top quark masses indicate that the Standard Model Higgs potential becomes unstable around $\Lambda_I \sim 10^{11}$ GeV. This instability is cosmologically relevant since quantum fluctuations during inflation can easily destabilize the electroweak vacuum if the Hubble parameter during inflation is larger than $\Lambda_I$ (as preferred by the recent BICEP2 measurement). We perform a careful study of the evolution of the Higgs field during inflation, obtaining different results from those currently in the literature. We consider both tunneling via a Coleman-de Luccia or Hawking-Moss instanton, valid when the scale of inflation is below the instability scale, as well as a statistical treatment via the Fokker-Planck equation appropriate in the opposite regime. We show that a better understanding of the post-inflation evolution of the unstable AdS vacuum regions is crucial for determining the eventual fate of the universe. If these AdS regions devour all of space, a universe like ours is indeed extremely unlikely without new physics to stabilize the Higgs potential; however, if these regions crunch, our universe survives, but inflation must last a few e-folds longer to compensate for the lost AdS regions. Lastly, we examine the effects of generic Planck-suppressed corrections to the Higgs potential, which can be sufficient to stabilize the electroweak vacuum during inflation.

In a recent paper arXiv:1107.5048, we discussed the correlation between the elastic neutralino-nucleon scattering cross section, constrained by dark matter direct detection experiments, and fine-tuning at tree-level in the electroweak symmetry breaking sector of the Minimal Supersymmetric Standard Model (MSSM). Here, we show that the correlation persists in the Next-to-Minimal Supersymmetric Standard Model (NMSSM), and its variant, lambda-SUSY. Both models are strongly motivated by the recent discovery of a 125 GeV Higgs-like particle. We also discuss the implications of the recently published bound on the direct detection cross section from 225 live days of XENON100 experiment. In both the MSSM and the NMSSM, most of the parameter space with fine-tuning less than 10% is inconsistent with the XENON100 bound. In lambda-SUSY, on the other hand, large regions of completely natural electroweak symmetry breaking are still allowed, primarily due to a parametric suppression of fine-tuning with large \lambda. The upcoming XENON1T experiment will be able to probe most of the parameter space with less than 1% fine-tuning in all three models.

The viability of neutralino dark matter as an explanation of the 130 GeV gamma ray signal from the Galactic Center recently observed by the Fermi Large Area Telescope is examined. It is found that the signal can be compatible with a sharp feature from internal bremsstrahlung from a mostly bino dark matter particle of mass around 145 GeV, augmented by a contribution from annihilation into gamma+Z via a small wino admixture. This scenario circumvents the problematic overproduction of lower energy continuum photons that plague line interpretations of this signal. Sleptons approximately degenerate in mass with the neutralino are required to enhance the internal bremsstrahlung feature.

We study theoretical implications of direct dark matter searches in the minimal supersymmetric standard model (MSSM). We assume that no accidental cancellations occur in the spin-independent elastic neutralino-quark scattering cross section, but do not impose any relations among the weak-scale MSSM parameters. We show that direct detection cross section below 10^-44 cm^2 requires the lightest supersymmetric particle (LSP) neutralino to be close to either a pure gaugino or pure Higgsino limit, with smaller cross sections correlated with smaller admixture of the subdominant components. The current XENON100 bound rules out essentially all models in which the lightest neutralino has the Higgsino fraction between 0.2 and 0.8. Furthermore, smaller direct detection cross sections correlate with stronger fine-tuning in the electroweak symmetry breaking sector. In the gaugino LSP scenario, the current XENON100 bound already implies some fine-tuning: for example, at least 10% tuning is required if the LSP mass is above 70 GeV. In both gaugino and Higgsino LSP scenarios, the direct dark matter searches currently being conducted and designed should lead to a discovery if no accidental cancellations or fine-tuning at a level below 1% is present.

Energetic antiprotons in cosmic rays can serve as an important indirect signature of dark matter. Conventionally, the antiproton flux from dark matter decays or annihilations is calculated by solving the transport equation with a space-independent diffusion coefficient within the diffusion zone of the galaxy, and assuming free propagation outside this zone. Antiproton sources outside of the diffusion zone are ignored. In reality, it is far more likely that the diffusion coefficient increases smoothly with distance from the disk, and the outlying part of the dark matter halo ignored in the conventional approach can be significant, containing as much as 90% of the galactic dark matter by mass in some models. We extend the conventional approach to address these issues. We obtain analytic approximations and numerical solutions for antiproton flux assuming that the diffusion coefficient increases exponentially with the distance from the disk, and including contributions from dark matter annihilations/decays in essentially the full dark matter halo. We find that the antiproton flux predicted in this model deviates from the conventional calculation for the same dark matter parameters by up to about 25%.

If the positron fraction and combined electron-positron flux excesses recently observed by PAMELA, Fermi and HESS have a dark matter origin, final state radiation (FSR) photons from dark matter annihilation into lepton-rich final states may be detected with observations of satellite dwarf galaxies of the Milky Way by ground-based atmospheric Cherenkov telescopes (ACTs). We find that current and near-future ACTs have excellent potential for such detection, although a discovery cannot be guaranteed due to large uncertainties in the distribution of dark matter within the dwarfs. We find that models predicting dark matter annihilation into two-lepton final states and those favoring four-lepton final states (as in, for example, "axion portal" models) can be reliably distinguished using the FSR photon spectrum once measured, and the dark matter particle mass can also be accurately determined.

If the positron fraction and combined electron-positron flux excesses recently observed by PAMELA, FERMI and HESS are due to dark matter annihilation into lepton-rich final states, the accompanying final state radiation (FSR) photons may be detected by ground-based atmospheric Cherenkov telescopes (ACTs). Satellite dwarf galaxies in the vicinity of the Milky Way are particularly promising targets for this search. We find that current and near-future ACTs have an excellent potential for discovering the FSR photons from dwarfs, although a discovery cannot be guaranteed due to large uncertainties in the fluxes resulting from lack of precise knowledge of dark matter distribution within the dwarfs. We also investigate the possibility of discriminating between different dark matter models based on the measured FSR photon spectrum. For typical parameters, we find that the ACTs can reliably distinguish models predicting dark matter annihilation into two-lepton final states from those favoring four-lepton final states (as in, for example, "axion portal" models). In addition, we find that the dark matter particle mass can also be determined from the FSR spectrum.

Energetic positrons produced in annihilation or decay of dark matter particles in the Milky Way can serve as an important indirect signature of dark matter. Computing the positron flux expected in a given dark matter model involves solving transport equations, which account for interaction of positrons with matter and galactic magnetic fields. Existing calculations solve the equations inside the diffusion zone, where galactic magnetic fields confine positrons, and assume vanishing positron density on the boundaries of this zone. However, in many models, a substantial fraction of the dark matter halo lies outside the diffusion zone. Positrons produced there can then enter the diffusion zone and get trapped, potentially reaching the Earth and increasing the expected flux. We calculate this enhancement for a variety of models. We also evaluate the expected enhancement of the flux of energetic photons produced by the inverse Compton scattering of the extra positrons on starlight and cosmic microwave background. We find maximal flux enhancements of order 20% in both cases.