In this short letter, we find that a magnetic transition dipole moment between tau and sterile neutrinos can account for the XENON1T excess events. Unlike the ordinary neutrino dipole moment, the introduction of the new sterile mass scale allows for astrophysical bounds to be suppressed. Interestingly, the best-fit regions that are compatible with the SN1987A imply either boron-8 or CNO neutrinos as the source flux. We find that sterile neutrinos of either $\sim$ 260 keV or in the $\sim$(500 - 800) keV mass range are capable of evading astrophysical constraints while being able to successfully explain the XENON1T event rate. The sterile neutrino in the best fit parameter space may have significant effects on big bang nucleosynthesis (BBN). We show the region in which a low reheating temperature of the Universe may allow the BBN constraints to be alleviated.

We study a Dark Matter (DM) model in which the dominant coupling to the standard model occurs through a neutrino-DM-scalar coupling. The new singlet scalar will generically have couplings to nuclei/electrons arising from renormalizable Higgs portal interactions. As a result the DM particle $X$ can convert into a neutrino via scattering on a target nucleus $\mathcal{N}$: $ X + \mathcal{N} \rightarrow \nu + \mathcal{N}$, leading to striking signatures at direct detection experiments. Similarly, DM can be produced in neutrino scattering events at neutrino experiments: $ \nu + \mathcal{N} \rightarrow X + \mathcal{N}$, predicting spectral distortions at experiments such as COHERENT. Furthermore, the model allows for late kinetic decoupling of dark matter with implications for small-scale structure. At low masses, we find that COHERENT and late kinetic decoupling produce the strongest constraints on the model, while at high masses the leading constraints come from DM down-scattering at XENON1T and Borexino. Future improvement will come from CE$\nu$NS data, ultra-low threshold direct detection, and rare kaon decays.

Quark nuggets are theoretical objects composed of approximately equal numbers of up, down, and strange quarks. They are also called strangelets, nuclearites, AQNs, slets, Macros, and MQNs. Quark nuggets are a candidate for dark matter, which has been a mystery for decades despite constituting ${\sim}85\%$ of the mass in the universe. Most previous models of quark nuggets have assumed no intrinsic magnetic field. However, Tatsumi found that quark nuggets may exist in magnetars as a ferromagnetic liquid with core magnetic field B between $10^{ 11}$ T and $10^{ 13}$ T. We apply that result to quark-nugget dark-matter and report results on aggregation of magnetized quark nuggets (MQNs) after formation from the quark-gluon plasma until expansion of the universe freezes out the mass distribution of $10^{ -24}$ kg to $10^{ 14}$ kg. Aggregation overcomes decay by weak interaction. Their aggregated mass distributions make MQNs consistent with all requirements for dark matter and fully compliant with Standard Model of Particle Physics. Observations narrow the range of B to between $10^{ 11}$ T and 3 $10^{ 13}$ T and indicate that geologic (craters in peat bogs) and planetary (seismic and radio-frequency emissions during passage through atmosphere) detectors can detect MQN dark matter.

The detection of the high-energy neutrino event, IceCube-170922A, demonstrated that multimessenger particle astrophysics triggered by neutrino alerts is feasible. We consider time delay signatures caused by secret neutrino interactions with the cosmic neutrino background and dark matter and suggest that these can be used as a novel probe of neutrino interactions beyond the Standard Model (BSM). The tests with BSM-induced neutrino echoes are distinct from existing constraints from the spectral modification and will be enabled by multimessenger observations of bright neutrino transients with future experiments such as IceCube-Gen2, KM3Net, and Hyper-Kamiokande. The constraints are complementary to those from accelerator and laboratory experiments and powerful for testing various particle models that explain tensions prevailing in the cosmological data.

This report summarizes the present status of neutrino non-standard interactions (NSI). After a brief overview, several aspects of NSIs are discussed, including connection to neutrino mass models, model-building and phenomenology of large NSI with both light and heavy mediators, NSI phenomenology in both short- and long-baseline neutrino oscillation experiments, neutrino cross-sections, complementarity of NSI with other low- and high-energy experiments, fits with neutrino oscillation and scattering data, DUNE sensitivity to NSI, effective field theory of NSI, as well as the relevance of NSI to dark matter and cosmology. We also discuss the open questions and interesting future directions that can be pursued by the community at large. This report is based on talks and discussions during the Neutrino Theory Network NSI workshop held at Washington University in St. Louis from May 29-31, 2019 (https://indico.cern.ch/event/812851/)

We study direct detection bounds on cosmic ray-upscattered dark matter in simplified models including light mediators. We find that the energy dependence in the scattering cross section is significant, and produces stronger bounds than previously found (which assumed constant cross sections) by many orders of magnitude at low DM mass. Finally, we compute the "neutrino-floor" that will limit future direct detection searches for cosmic ray-upscattered dark matter. While we focus on vector interactions for illustration, we emphasize that the energy dependence is critical in determining accurate bounds on any particle physics model of Dark Matter-CR interactions from experimental data on this scenario.

The ANITA balloon experiment was designed to detect radio signals initiated by neutrinos and cosmic ray air showers. These signals are typically discriminated by the polarization and phase inversions of the radio signal. The reflected signal from cosmic rays suffer phase inversion compared to a direct tau neutrino event. In this paper we study sub-surface reflection, which can occur without phase inversion, in the context of the two anomalous up-going events reported by ANITA. We find that subsurface layers and firn density inversions may plausibly account for the events, while ice fabric layers and wind ablation crusts could also play a role. This hypothesis can be tested with radar surveying of the Antarctic region in the vicinity of the anomalous ANITA events. Future experiments should not use phase inversion as a sole criterion to discriminate between downgoing and upgoing events, unless the subsurface reflection properties are well understood.

The ANITA balloon experiment has recently observed several ~EeV cascade events at an angle below the horizon that renders any Standard Model (SM) interpretation unlikely as the Earth is significantly opaque to all SM particles at such energies. In this paper, we study a sterile neutrino interpretation of these events, calculating the angular acceptance of cascades and the relative sensitivities of several experiments to a cascade initiated by an EeV sterile neutrino. We find that ANITA is uniquely sensitive to this type of upward directed cascade signal over a wide portion of the sky and from the direction of the two observed events has a transient acceptance roughly equivalent to that of the IceCube experiment.

The Laser Interferometer Gravitational-Wave Observatory (LIGO) has recently discovered gravitational waves (GWs) from its first neutron star-neutron star merger at a distance of $\sim 40$~Mpc from the Earth. The associated electromagnetic (EM) detection of the event, including the short gamma-ray burst within $\Delta t \sim 2$~s after the GW arrival, can be used to test various aspects of sources physics and GW propagation. Using GW170817 as the first GW-EM example, we show that this event provides a stringent direct test that GWs travel at the speed of light. The gravitational potential of the Milky Way provides a potential source of Shapiro time delay difference between the arrival of photons and GWs, and we demonstrate that the nearly coincident detection of the GW and EM signals can yield strong limits on anomalous gravitational time delay, through updating the previous limits taking into account details of Milky Way's gravitational potential. Finally, we also obtain an intriguing limit on the size of the prompt emission region of GRB 170817A, and discuss implications for the emission mechanism of short gamma-ray bursts.

A variety of new physics models allows for neutrinos to up-scatter into heavier states. If the incident neutrino is energetic enough, the heavy neutrino may travel some distance before decaying. In this work, we consider the atmospheric neutrino flux as a source of such events. At IceCube, this would lead to a "double-bang" (DB) event topology, similar to what is predicted to occur for tau neutrinos at ultra-high energies. The DB event topology has an extremely low background rate from coincident atmospheric cascades, making this a distinctive signature of new physics. Our results indicate that IceCube should already be able to derive new competitive constraints on models with GeV-scale sterile neutrinos using existing data.

A well-studied possibility is that dark matter may reside in a sector secluded from the Standard Model, except for the so-called photon portal: kinetic mixing between the ordinary and dark photons. Such interactions can be probed at dark matter direct detection experiments, and new experimental techniques involving detection of dark matter-electron scattering offer new sensitivity to sub-GeV dark matter. Typically however it is implicitly assumed that the dark matter is not altered as it traverses the Earth to arrive at the detector. In this paper we study in detail the effects of terrestrial stopping on dark photon models of dark matter, and find that they significantly reduce the sensitivity of XENON10 and DAMIC. In particular we find that XENON10 only excludes masses in the range (5-3000) MeV while DAMIC only probes (20-50) MeV. Their corresponding cross section sensitivity is reduced to a window of cross sections between $(5\times 10^{-38}-10^{-30})~{\rm cm}^{2}$ for XENON10 and a small window around $\sim 10^{-31}~{\rm cm}^{2}$ for DAMIC. We also examine implications for a future DAMIC run.

Sterile neutrinos at the eV scale have long been studied in the context of anomalies in short baseline neutrino experiments. Their cosmology can be made compatible with our understanding of the early Universe provided the sterile neutrino sector enjoys a nontrivial dynamics with exotic interactions, possibly providing a link to the Dark Matter (DM) puzzle. Interactions between DM and neutrinos have also been proposed to address the long-standing "missing satellites" problem in the field of large scale structure formation. Motivated by these considerations, in this paper we discuss realistic scenarios with light steriles coupled to DM. We point out that within this framework active neutrinos acquire an effective coupling to DM that manifests itself as a new matter potential in the propagation within a medium of asymmetric DM. Assuming that at least a small fraction of asymmetric DM has been captured by the Sun, we show that a sizable region of the parameter space of these scenarios can be probed by solar neutrino experiments, especially in the regime of small couplings and light mediators where all other probes become inefficient. In the latter regime these scenarios behave as familiar $3+1$ models in all channels except for solar data, where a Solar Dark MSW effect takes place. Solar Dark MSW is characterized by modifications of the most energetic $^8$B and CNO neutrinos, whereas the other fluxes remain largely unaffected.

We present a targeted search for blazar flux-correlated high-energy ($\varepsilon_\nu > 1$ TeV) neutrinos from six bright northern blazars, using the public database of northern-hemisphere neutrinos detected during "IC40" 40-string operations of the IceCube neutrino observatory (April 2008 to May 2009). Our six targeted blazars are subjects of long-term monitoring campaigns by the VERITAS TeV gamma-ray observatory. We use the publicly-available VERITAS lightcurves to identify periods of excess and flaring emission. These predefined intervals serve as our "active temporal windows" in a search for an excess of neutrinos, relative to Poisson fluctuations of the near-isotropic atmospheric neutrino background which dominates at these energies. After defining the parameters of an optimized search, we confirm the expected Poisson behavior with Monte Carlo simulations prior to testing for excess neutrinos in the actual data. We make two searches: One for excess neutrinos associated with the bright flares of Mrk 421 that occurred during the IC40 run, and one for excess neutrinos associated with the brightest emission periods of five other blazars (Mrk 501, 1ES 0805+524, 1ES 1218+304, 3C66A, and W Comae), all significantly fainter than the Mrk 421 flares. We find no significant excess of neutrinos from the preselected blazar directions during the selected temporal windows. We derive 90%-confidence upper limits on the number of expected flux-associated neutrinos from each search. These limits are consistent with previous point-source searches and Fermi GeV flux-correlated searches. Our upper limits are sufficiently close to the physically-interesting regime that we anticipate future analyses using already-collected data will either constrain models or yield discovery of the first blazar-associated high-energy neutrinos.

Dark Matter (DM) may have a relic density that is in part determined by a particle/antiparticle asymmetry, much like baryons. If this is the case, it can accumulate in stars like the Sun to sizable number densities and annihilate to Standard Model (SM) particles including neutrinos. We show that the combination of neutrino telescope and direct detection data can be used in conjunction to determine or constrain the DM asymmetry from data. Depending on the DM mass, the current neutrino data from Super-K and IceCube give powerful constraints on asymmetric DM unless its fractional asymmetry is $\lesssim 10^{-2}$. Future neutrino telescopes and detectors like Hyper-K and KM3NeT can search for the resulting signal of high-energy neutrinos from the center of the Sun. The observation of such a flux yields information on both the DM-nucleus cross section but also on the relative abundances of DM and anti-DM.

We examine a framework with light new physics, which couples to the Standard Model only via neutrino mixing. Taking the hints from the short-baseline anomalies seriously and combining them with modern cosmological data and recent IceCube measurements, we obtain surprisingly effective constraints on the hidden force: keV $\lesssim M \lesssim0.3$ GeV for the mediator mass and $g_{h}>10^{-6}-10^{-3}$ for the coupling constant. Flavor equilibration between the hidden and active neutrinos can be delayed until temperatures of $\sim 1$ MeV, but not below $\sim 100$ keV. This scenario can be tested with next-generation Cosmic Microwave Background, IceCube, and oscillation experiments.

Recently, the direct detection of gravitational waves from black hole (BH) mergers was announced by the Advanced LIGO Collaboration. Multi-messenger counterparts of stellar-mass BH mergers are of interest, and it had been suggested that a small disk or celestial body may be involved in the binary of two BHs. To test such possibilities, we consider the fate of a wind powered by an active mini-disk in a relatively short, super-Eddington accretion episode onto a BH with ~10-100 solar masses. We show that its thermal emission could be seen as a fast optical transient with the duration from hours to days. We also find that the coasting outflow forms external shocks due to interaction with the interstellar medium, whose synchrotron emission might be expected in the radio band on a time scale of years. Finally, we also discuss a possible jet component and the associated high-energy neutrino emission as well as ultra-high-energy cosmic-ray acceleration.

The flavor of cosmic neutrinos may help unveil their sources and could reveal the presence of new physics in the neutrino sector. We consider the impacts of next-generation neutrino detectors, including the planned upgrade to neutrino detector--IceCube-Gen2, which is well-positioned to make dramatic improvements in both flavor and spectral measurements. We show that various models in neutrino physics beyond the Standard Model, such as neutrino decay, pseudo-Dirac states, and neutrino self-scattering, may be found or strongly constrained at IceCube-Gen2 and KM3NeT. We find that the additional flavor discriminants given by Glashow resonance events and so-called "double-bang" topologies improve the ability to access the flavor of the cosmic high-energy neutrinos and probe the BSM physics. In addition, although the details depend on source properties, Glashow resonance events have the additional feature of being able to inform us of the relative strengths of neutrino and antineutrino emission, which may help us discriminate astrophysical scenarios.

We investigate the direct detection phenomenology of a class of dark matter (DM) models in which DM does not directly interact with nuclei, {but rather} the products of its annihilation do. When these annihilation products are very light compared to the DM mass, the scattering in direct detection experiments is controlled by relativistic kinematics. This results in a distinctive recoil spectrum, a non-standard and or even absent annual modulation, and the ability to probe DM masses as low as a $\sim$10 MeV. We use current LUX data to show that experimental sensitivity to thermal relic annihilation cross sections has already been reached in a class of models. Moreover, the compatibility of dark matter direct detection experiments can be compared directly in $E_{{\rm min}}$ space without making assumptions about DM astrophysics, mass, or scattering form factors. Lastly, when DM has direct couplings to nuclei, the limit from annihilation to relativistic particles in the Sun can be stronger than that of conventional non-relativistic direct detection by more than three orders of magnitude for masses in a 2-7 GeV window.

We study self-interacting dark matter coupled to the Standard Model via the Higgs portal. We consider a scenario where dark matter is a thermal relic with strong enough self interactions that can alleviate the problems of collisionless cold dark matter. We study constraints from direct detection searches, the LHC, and Big Bang nucleosynthesis. We show that the tension between these constraints and the need for sufficiently strong self-interactions with light mediators can be alleviated by coupling the mediator to either active or sterile neutrinos. Future direct detection data offers great potential and can be used to find evidence of a light mediator and verify that dark matter scatters via long-range self-interactions.

It has been suggested that the baseline scenario of collisionless cold dark matter over-predicts the numbers of satellite galaxies, as well as the dark matter (DM) densities in galactic centers. This apparent lack of structure at small scales can be accounted for if one postulates neutrino-DM and DM-DM interactions mediated by light O(MeV) force carriers. In this letter, we consider a simple, consistent model of neutrinophilic DM with these features where DM and a "secluded" SM-singlet neutrino species are charged under a new $U(1)$ gauge symmetry. An important ingredient of this model is that the secluded sector couples to the Standard Model fields only through neutrino mixing. We observe that the secluded and active neutrinos recouple, leading to a large relic secluded neutrino population. This relic population can prevent small-scale halos from collapsing, while at same time significantly modifying the optical depth of ultra-high-energy neutrinos recently observed at Icecube. We find that the bulk of the parameter space accommodating an (a)symmetric thermal relic has potentially observable consequences for the IceCube high energy signal, with some of the parameter ranges already ruled out by the existing data. Future data may confirm this mechanism if either spectral absorption features or correlations with nearby sources are observed.

The relic abundance of particle and antiparticle dark matter (DM) need not be vastly different in thermal asymmetric dark matter (ADM) models. By considering the effect of a primordial asymmetry on the thermal Boltzmann evolution of coupled DM and anti-DM, we derive the requisite annihilation cross section. This is used in conjunction with CMB and Fermi-LAT gamma-ray data to impose a limit on the number density of anti-DM particles surviving thermal freeze-out. When the extended gamma-ray emission from the Galactic Center is reanalyzed in a thermal ADM framework, we find that annihilation into $\tau$ leptons prefer anti-DM number densities 1-4$\%$ that of DM while the $b$-quark channel prefers 50-100$\%$.

We study dark matter that inelastically scatters and de-excites in direct detection experiments, as an interpretation of the CDMS-Si events in light of the recent LUX data. The constraints from LUX and XENON10 require the mass-splitting between the DM excited and de-excited states to be $|\delta| \gtrsim 50$ keV. At the same time, the CDMS-Si data itself do not allow for a consistent DM interpretation for mass splittings larger than $|\delta| \sim $200 keV. We find that a low threshold analysis will be needed to rule out this interpretation of the CDMS-Si events. In a simple model with a kinetically mixed dark photon, we show that the CDMS-Si rate and the thermal relic abundance can both be accommodated.

Isospin-violating dark matter (IVDM) has been proposed as a viable scenario to reconcile conflicting positive and null results from direct detection dark matter experiments. We show that the lowest-order dark matter-nucleus scattering rate can receive large and nucleus-dependent corrections at next-to-leading order (NLO) in the chiral expansion. The size of these corrections depends on the specific couplings of dark matter to quark flavors and gluons. In general the full NLO dark-matter-nucleus cross-section is not adequately described by just the zero-energy proton and neutron couplings. These statements are concretely illustrated in a scenario where the dark matter couples to quarks through scalar operators. We find the canonical IVDM scenario can reconcile the null XENON and LUX results and the recent CDMS-Si findings provided its couplings to second and third generation quarks either lie on a special line or are suppressed. Equally good fits with new values of the neutron-to-proton coupling ratio are found in the presence of nonzero heavy quark couplings. CDMS-Si remains in tension with LUX and XENON10/100 but is not excluded.

We study the capabilities of the MAJORANA DEMONSTRATOR, a neutrinoless double-beta decay experiment currently under construction at the Sanford Underground Laboratory, as a light WIMP detector. For a cross section near the current experimental bound, the MAJORANA DEMONSTRATOR should collect hundreds or even thousands of recoil events. This opens up the possibility of simultaneously determining the physical properties of the dark matter and its local velocity distribution, directly from the data. We analyze this possibility and find that allowing the dark matter velocity distribution to float considerably worsens the WIMP mass determination. This result is traced to a previously unexplored degeneracy between the WIMP mass and the velocity dispersion. We simulate spectra using both isothermal and Via Lactea II velocity distributions and comment on the possible impact of streams. We conclude that knowledge of the dark matter velocity distribution will greatly facilitate the mass and cross section determination for a light WIMP.

The mass of primordial dark matter (DM) protohalos remains unknown. However, the missing satellites problem may be an indication that they are quite large. In this paper, we use effective field theory to map constraints on dark matter-SM interactions into limits on the mass of DM protohalos. Given that leptons remain in the thermal bath until late times, we focus on their interactions with DM. To illustrate the method, we use the null results of LEP missing energy searches along with Fermi-LAT searches for DM annihilation in nearby dwarf galaxies, to derive limits on the protohalo mass, $\lesssim (10^{-6}-10^{-1}) M_{\odot}$, with the range depending on the DM mass and the operator. Thus, if DM is to remain thermally coupled until late times and account for the missing satellites, charged lepton interactions are insufficient. This motivates neutrinophilic DM, which can have protohalo masses orders of magnitude larger, with constraints arising from Planck, IceCube and unpublished Super-K data. We show that effective neutrinophilic models offer a solution to the missing satellites problem for sub-GeV DM masses with larger than WIMP-sized annihilation cross sections.

If dark matter interacts with quarks or gluons, the mediator of these interactions is either directly accessible at the LHC or is so heavy that its effects are encoded in contact operators. We find that the self-consistency of a contact operator description at the LHC implies bounds on the mediator scale stronger than those found from missing energy searches. This translates into spin-independent elastic scattering cross-sections at a level < 10^-41 cm^2, with direct implications for the DAMA, CoGeNT, and CRESST-II anomalies. We then carefully explore the potential of monojet searches in the light mediator limit, focusing on a Z' model with arbitrary couplings to quarks and dark matter. We find that the Tevatron data currently provides the most stringent bounds for dark matter and Z' masses below 100 GeV, and that these searches can constrain models for the DAMA, CoGeNT, and CRESST-II anomalies only if the mediator can decay to a pair of dark matter particles.

The similarity of the visible and dark matter abundances indicates that they may originate via the same mechanism. If both the dark and the visible matter are charged under a generalized baryon number which remains always conserved, then the asymmetry of the visible sector may be compensated by an asymmetry in the dark sector. We show how the separation of baryonic and antibaryonic charge can originate in the vacuum, via the Affleck-Dine mechanism, due to the breaking of a symmetry orthogonal to the baryon number. Symmetry restoration in the current epoch guarantees the individual stability of the two sectors.

In existing dark matter models with global symmetries the relic abundance of dark matter is either equal to that of anti-dark matter (thermal WIMP), or vastly larger, with essentially no remaining anti-dark matter (asymmetric dark matter). By exploring the consequences of a primordial asymmetry on the coupled dark matter and anti-dark matter Boltzmann equations we find large regions of parameter space that interpolate between these two extremes. Interestingly, this new asymmetric WIMP framework can accommodate a wide range of dark matter masses and annihilation cross sections. The present-day dark matter population is typically asymmetric, but only weakly so, such that indirect signals of dark matter annihilation are not completely suppressed. We apply our results to existing models, noting that upcoming direct detection experiments will constrain a large region of the relevant parameter space.

We suggest the existence of a fundamental connection between baryonic and dark matter. This is motivated by both the stability of these two types of matter as well as the observed similarity of their present-day densities. A unified genesis of baryonic and dark matter is natural in models in which the baryon number is promoted to a spontaneously broken local gauge symmetry. This is illustrated in a specific class of SUSY models using the Affleck-Dine mechanism. The dark matter candidate in these scenarios is charged under the baryon gauge symmetry and must have a mass around the GeV scale to give the correct present-day abundance. We discuss constraints from B-factories, LEP, mono-jet searches at the Tevatron, and dark matter direct detection experiments. A baryonic dark force is shown to be consistent with all data for mediators as light as the GeV scale.

Although inflationary models generically predict a flat spectrum of gravitational waves, we point out a general process that produces a sharply peaked spectrum of gravitational radiation. This process is generic for inflationary models with a complex inflaton field which couples to fermions. In particular, for chaotic models these may be the most extreme gravitational waves in the Universe with a very large energy density fraction 10^-9 and ultra-high frequency, 10^10 Hz. Although not amenable to space based interferometers, the signal from this model may be detectable by future table top experiments.

The gauge-mediated model of supersymmetry breaking implies that stable non-topological solitons, Q-balls, could form in the early universe and comprise the dark matter. It is shown that the inclusion of the effects from gravity-mediation set an upper limit on the size of Q-balls. When in a dense baryonic environment Q-balls grow until reaching this limiting size at which point they fragment into two equal-sized Q-balls. This Q-splitting process will rapidly destroy a neutron star that absorbs even one Q-ball. The new limits on Q-ball dark matter require an ultralight gravitino m_3/2 < keV, naturally avoiding the gravitino overclosure problem, and providing the MSSM with a dark matter candidate where gravitino dark matter is not viable.

Supersymmetry implies that stable non-topological solitons, Q-balls, could form in the early universe and could make up all or part of dark matter. We show that the relic Q-balls passing through Earth can produce a detectable neutrino flux. The peculiar zenith angle dependence and a small annual modulation of this flux can be used as signatures of dark-matter Q-balls.