Multi-messenger astrophysics is a fast-growing, interdisciplinary field that combines data, which vary in volume and speed of data processing, from many different instruments that probe the Universe using different cosmic messengers: electromagnetic waves, cosmic rays, gravitational waves and neutrinos. In this Expert Recommendation, we review the key challenges of real-time observations of gravitational wave sources and their electromagnetic and astroparticle counterparts, and make a number of recommendations to maximize their potential for scientific discovery. These recommendations refer to the design of scalable and computationally efficient machine learning algorithms; the cyber-infrastructure to numerically simulate astrophysical sources, and to process and interpret multi-messenger astrophysics data; the management of gravitational wave detections to trigger real-time alerts for electromagnetic and astroparticle follow-ups; a vision to harness future developments of machine learning and cyber-infrastructure resources to cope with the big-data requirements; and the need to build a community of experts to realize the goals of multi-messenger astrophysics.

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.

This report provides an overview of recent work that harnesses the Big Data Revolution and Large Scale Computing to address grand computational challenges in Multi-Messenger Astrophysics, with a particular emphasis on real-time discovery campaigns. Acknowledging the transdisciplinary nature of Multi-Messenger Astrophysics, this document has been prepared by members of the physics, astronomy, computer science, data science, software and cyberinfrastructure communities who attended the NSF-, DOE- and NVIDIA-funded "Deep Learning for Multi-Messenger Astrophysics: Real-time Discovery at Scale" workshop, hosted at the National Center for Supercomputing Applications, October 17-19, 2018. Highlights of this report include unanimous agreement that it is critical to accelerate the development and deployment of novel, signal-processing algorithms that use the synergy between artificial intelligence (AI) and high performance computing to maximize the potential for scientific discovery with Multi-Messenger Astrophysics. We discuss key aspects to realize this endeavor, namely (i) the design and exploitation of scalable and computationally efficient AI algorithms for Multi-Messenger Astrophysics; (ii) cyberinfrastructure requirements to numerically simulate astrophysical sources, and to process and interpret Multi-Messenger Astrophysics data; (iii) management of gravitational wave detections and triggers to enable electromagnetic and astro-particle follow-ups; (iv) a vision to harness future developments of machine and deep learning and cyberinfrastructure resources to cope with the scale of discovery in the Big Data Era; (v) and the need to build a community that brings domain experts together with data scientists on equal footing to maximize and accelerate discovery in the nascent field of Multi-Messenger Astrophysics.

Many motivated extensions of the Standard Model predict the existence of cosmic strings. Gravitational waves originating from the dynamics of the resulting cosmic string network have the ability to probe many otherwise inaccessible properties of the early universe. In this study we show how the spectrum of gravitational waves from a cosmic string network can be used to test the equation of state of the early universe prior to Big Bang Nucleosynthesis (BBN). We also demonstrate that current and planned gravitational wave detectors such as LIGO, LISA, DECIGO/BBO, and ET/CE have the potential to detect signals of a non-standard pre-BBN equation of state and evolution of the early universe (e.g., early non-standard matter domination or kination domination) or new degrees of freedom active in the early universe beyond the sensitivity of terrestrial collider experiments and cosmic microwave background measurements.

Cosmic strings are generic cosmological predictions of many extensions of the Standard Model of particle physics, such as a $U(1)^\prime$ symmetry breaking phase transition in the early universe or remnants of superstring theory. Unlike other topological defects, cosmic strings can reach a scaling regime that maintains a small fixed fraction of the total energy density of the universe from a very early epoch until today. If present, they will oscillate and generate gravitational waves with a frequency spectrum that imprints the dominant sources of total cosmic energy density throughout the history of the universe. We demonstrate that current and future gravitational wave detectors, such as LIGO and LISA, could be capable of measuring the frequency spectrum of gravitational waves from cosmic strings and discerning the energy composition of the universe at times well before primordial nucleosynthesis and the cosmic microwave background where standard cosmology has yet to be tested. This work establishes a benchmark case that gravitational waves may provide an unprecedented, powerful tool for probing the evolutionary history of the very early universe.

We analyse a simple extension of the SM with just an additional scalar singlet coupled to the Higgs boson. We discuss the possible probes for electroweak baryogenesis in this model including collider searches, gravitational wave and direct dark matter detection signals. We show that a large portion of the model parameter space exists where the observation of gravitational waves would allow detection while the indirect collider searches would not.

New results are reported from the operation of the PICO-60 dark matter detector, a bubble chamber filled with 52 kg of C$_3$F$_8$ located in the SNOLAB underground laboratory. As in previous PICO bubble chambers, PICO-60 C$_3$F$_8$ exhibits excellent electron recoil and alpha decay rejection, and the observed multiple-scattering neutron rate indicates a single-scatter neutron background of less than 1 event per month. A blind analysis of an efficiency-corrected 1167-kg-day exposure at a 3.3-keV thermodynamic threshold reveals no single-scattering nuclear recoil candidates, consistent with the predicted background. These results set the most stringent direct-detection constraint to date on the WIMP-proton spin-dependent cross section at 3.4 $\times$ 10$^{-41}$ cm$^2$ for a 30-GeV$\thinspace$c$^{-2}$ WIMP, more than one order of magnitude improvement from previous PICO results.

The gravitational waves measured at LIGO are presumed here to come from merging primordial black holes. We ask how these primordial black holes could arise through inflationary models while not conflicting with current experiments. Among the approaches that work, we investigate the opportunity for corroboration through experimental probes of gravitational waves at pulsar timing arrays. We provide examples of theories that are already ruled out, theories that will soon be probed, and theories that will not be tested in the foreseeable future. The models that are most strongly constrained are those with a relatively broad primordial power spectrum.

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.

We consider various models realizing baryogenesis during the electroweak phase transition (EWBG). Our focus is their possible detection in future collider experiments and possible observation of gravitational waves emitted during the phase transition. We also discuss the possibility of a non-standard cosmological history which can facilitate EWBG. We show how acceptable parameter space can be extended due to such a modification and conclude that next generation precision experiments such as the ILC will be able to confirm or falsify many models realizing EWBG. We also show that, in general, collider searches are a more powerful probe than gravitational wave searches. However, observation of a deviation from the SM without any hints of gravitational waves can point to models with modified cosmological history that generically enable EWBG with weaker phase transition and thus, smaller GW signals.

We study the impact on electroweak baryogenesis from a swifter cosmological expansion induced by dark matter. We detail the experimental bounds that one can place on models that realize it, and we investigate the modifications of these bounds that result from a non-standard cosmological history. The modifications can be sizeable if the expansion rate of the Universe increases by several orders of magnitude. We illustrate the impact through the example of scalar field dark matter, which can alter the cosmological history enough to enable a strong-enough first-order phase transition in the Standard Model when it is supplemented by a dimension six operator directly modifying the Higgs boson potential. We show that due to the modified cosmological history, electroweak baryogenesis can be realized, while keeping deviations of the triple Higgs coupling below HL-LHC sensitivies. The required scale of new physics to effectuate a strong-enough first order phase transition can change by as much as twenty percent as the expansion rate increases by six orders of magnitude.

This draft report summarizes and details the findings, results, and recommendations derived from the ASCR/HEP Exascale Requirements Review meeting held in June, 2015. The main conclusions are as follows. 1) Larger, more capable computing and data facilities are needed to support HEP science goals in all three frontiers: Energy, Intensity, and Cosmic. The expected scale of the demand at the 2025 timescale is at least two orders of magnitude -- and in some cases greater -- than that available currently. 2) The growth rate of data produced by simulations is overwhelming the current ability, of both facilities and researchers, to store and analyze it. Additional resources and new techniques for data analysis are urgently needed. 3) Data rates and volumes from HEP experimental facilities are also straining the ability to store and analyze large and complex data volumes. Appropriately configured leadership-class facilities can play a transformational role in enabling scientific discovery from these datasets. 4) A close integration of HPC simulation and data analysis will aid greatly in interpreting results from HEP experiments. Such an integration will minimize data movement and facilitate interdependent workflows. 5) Long-range planning between HEP and ASCR will be required to meet HEP's research needs. To best use ASCR HPC resources the experimental HEP program needs a) an established long-term plan for access to ASCR computational and data resources, b) an ability to map workflows onto HPC resources, c) the ability for ASCR facilities to accommodate workflows run by collaborations that can have thousands of individual members, d) to transition codes to the next-generation HPC platforms that will be available at ASCR facilities, e) to build up and train a workforce capable of developing and using simulations and analysis to support HEP scientific research on next-generation systems.

New data are reported from a second run of the 2-liter PICO-2L C$_3$F$_8$ bubble chamber with a total exposure of 129$\,$kg-days at a thermodynamic threshold energy of 3.3$\,$keV. These data show that measures taken to control particulate contamination in the superheated fluid resulted in the absence of the anomalous background events observed in the first run of this bubble chamber. One single nuclear-recoil event was observed in the data, consistent both with the predicted background rate from neutrons and with the observed rate of unambiguous multiple-bubble neutron scattering events. The chamber exhibits the same excellent electron-recoil and alpha decay rejection as was previously reported. These data provide the most stringent direct detection constraints on weakly interacting massive particle (WIMP)-proton spin-dependent scattering to date for WIMP masses $<$ 50$\,$GeV/c$^2$.

New data are reported from the operation of the PICO-60 dark matter detector, a bubble chamber filled with 36.8 kg of CF$_3$I and located in the SNOLAB underground laboratory. PICO-60 is the largest bubble chamber to search for dark matter to date. With an analyzed exposure of 92.8 livedays, PICO-60 exhibits the same excellent background rejection observed in smaller bubble chambers. Alpha decays in PICO-60 exhibit frequency-dependent acoustic calorimetry, similar but not identical to that reported recently in a C$_3$F$_8$ bubble chamber. PICO-60 also observes a large population of unknown background events, exhibiting acoustic, spatial, and timing behaviors inconsistent with those expected from a dark matter signal. These behaviors allow for analysis cuts to remove all background events while retaining $48.2\%$ of the exposure. Stringent limits on weakly interacting massive particles interacting via spin-dependent proton and spin-independent processes are set, and most interpretations of the DAMA/LIBRA modulation signal as dark matter interacting with iodine nuclei are ruled out.

The cosmic evolution of a dark matter model which behaves relativistically in the early Universe is explored. Dark matter is described as a complex scalar field, whose earliest evolution is characterized by a stiff equation of state ($p \simeq \rho$). In this phase, it is the dominant component in the Universe. We present constraints from Big Bang nucleosynthesis and primordial gravity waves from inflation. Also, we study how the associated enhanced expansion rate due to the stiff phase might facilitate a first-order electroweak symmetry breaking phase transition, in light of the recently measured value of the Higgs boson mass.

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.

New data are reported from the operation of a 2-liter C$_3$F$_8$ bubble chamber in the 2100 meter deep SNOLAB underground laboratory, with a total exposure of 211.5 kg-days at four different recoil energy thresholds ranging from 3.2 keV to 8.1 keV. These data show that C3F8 provides excellent electron recoil and alpha rejection capabilities at very low thresholds, including the first observation of a dependence of acoustic signal on alpha energy. Twelve single nuclear recoil event candidates were observed during the run. The candidate events exhibit timing characteristics that are not consistent with the hypothesis of a uniform time distribution, and no evidence for a dark matter signal is claimed. These data provide the most sensitive direct detection constraints on WIMP-proton spin-dependent scattering to date, with significant sensitivity at low WIMP masses for spin-independent WIMP-nucleon scattering.

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.

Flat directions are a generic feature of the scalar potential in supersymmetric gauge field theories. They can arise, for example, from D-terms associated with an extra abelian gauge symmetry. Even when supersymmetry is broken softly, there often remain directions in the scalar field space along which the potential is almost flat. Upon breaking a gauge symmetry along one of these almost flat directions, cosmic strings may form. Relative to the standard cosmic string picture based on the abelian Higgs model, these flat-direction cosmic strings have the extreme Type-I properties of a thin gauge core surrounded by a much wider scalar field profile. We perform a comprehensive study of the microscopic, macroscopic, and observational characteristics of this class of strings. We find many differences from the standard string scenario, including stable higher winding mode strings, the dynamical formation of higher mode strings from lower ones, and a resultant multi-tension scaling string network in the early universe. These strings are only moderately constrained by current observations, and their gravitational wave signatures may be detectable at future gravity wave detectors. Furthermore, there is the interesting but speculative prospect that the decays of cosmic string loops in the early universe could be a source of ultra-high energy cosmic rays or non-thermal dark matter. We also compare the observational signatures of flat-direction cosmic strings with those of ordinary cosmic strings as well as (p,q) cosmic strings motivated by superstring theory.

We study the possibility of a first-order electroweak phase transition (EWPT) due to a dimension-six operator in the effective Higgs potential. In contrast with previous attempts to make the EWPT strongly first-order as required by electroweak baryogenesis, we do not rely on large one-loop thermally generated cubic Higgs interactions. Instead, we augment the Standard Model (SM) effective theory with a dimension-six Higgs operator. This addition enables a strong first-order phase transition to develop even with a Higgs boson mass well above the current direct limit of 114 GeV. The phi^6 term can be generated for instance by strong dynamics at the TeV scale or by integrating out heavy particles like an additional singlet scalar field. We discuss conditions to comply with electroweak precision constraints, and point out how future experimental measurements of the Higgs self couplings could test the idea.

Recently the HEAT balloon experiment has confirmed an excess of high-energy positrons in cosmic rays. They could come from annihilation of dark matter in the galactic halo. We discuss expectations for the positron signal in cosmic rays from the lightest superpartner. The simplest interpretations are incompatible with the size and shape of the excess if the relic LSPs evolved from thermal equilbrium. Non-thermal histories can describe a sufficient positron rate. Reproducing the energy spectrum is more challenging, but perhaps possible. The resulting light superpartner spectrum is compatible with collider physics, the muon anomalous magnetic moment, Z-pole electroweak data, and other dark matter searches.

Cosmic ray antiprotons can originate from dark matter annihilating into quarks that subsequently decay into antiprotons. Evaporation of primordial black holes also can produce a significant antiproton flux. Since the spectrum of secondary antiprotons from cosmic ray interactions peaks at ~ 2 GeV and goes down sharply at lower energy, there is a window at energies < 1 GeV in which to look for excess antiprotons as a signature of these exotic antiproton sources. However, in the vicinity of the Earth low energy particles are strongly modulated by the solar wind, which makes any analysis ambiguous. The adverse effects of the solar wind can be avoided by placing a low energy antiproton spectrometer aboard an interstellar probe. The theoretical predictions are reviewed and the preliminary design of a light-weight, low-power instrument to make the measurements and a summary of the anticipated results are given in this paper.

Interstellar antiproton fluxes can arise from dark matter annihilating or decaying into quarks or gluons that subsequently fragment into antiprotons. Evaporation of primordial black holes also can produce a significant antiproton cosmic-ray flux. Since the background of secondary antiprotons from spallation has an interstellar energy spectrum that peaks at $\sim 2\gev$ and falls rapidly for energies below this, low-energy measurements of cosmic antiprotons are useful in the search for exotic antiproton sources. However, measurement of the flux near the earth is challenged by significant uncertainties from the effects of the solar wind. We suggest evading this problem and more effectively probing dark-matter signals by placing an antiproton spectrometer aboard an interstellar probe currently under discussion. We address the experimental challenges of a light, low-power-consuming detector, and present an initial design of such an instrument. This experimental effort could significantly increase our ability to detect, and have confidence in, a signal of exotic, nonstandard antiproton sources. Furthermore, solar modulation effects in the heliosphere would be better quantified and understood by comparing results to inverse modulated data derived from existing balloon and space-based detectors near the earth.

A weakly interacting dark-matter particle may be difficult to discover at an accelerator because it either (1) is too massive, (2) has no standard-model gauge interactions, or (3) is almost degenerate with other states. In each of these cases, searches for annihilation products in the Galactic halo are useful probes of dark-matter properties. Using the example of supersymmetric dark matter, I discuss how astrophysical searches for dark matter may provide discovery and mass information inaccessible to collider physics programs such as the Tevatron and LHC.

Active Galactic Nuclei (AGN) may emit highly collimated and intense jets of relativistic electrons which upscatter ambient photons. These electrons can also scatter off the cold dark matter halo of the galaxy to produce high energy photons which have a more isotropic signature than the upscattered photons from QED processes. We propose to look for these high energy photons coming from AGN as a method to detect dark matter. As a primary example we work out the expected signal from electrons scattering off the lightest supersymmetric partner into a photon plus selectron. Using the optimistic side of astrophysical uncertainties, we still find the signal from M87 or Centaurus A, two close-by AGN, smaller than the sensitivities expected of the currently proposed photon detectors. However, long running photon detectors and future detectors of higher sensitivity might be able to distinguish a signal from AGN sources. In order to have confidence that new physics sources are discernible, we also emphasize the importance of multi-wavelength studies of AGN with varying jet axis orientation to the earth.

Motivated by the supersymmetric interpretation of the CDF ee\gamma\gamma plus missing energy event and the reported Z->bb excess at LEP, we analyze the Higgsino as a cold dark matter candidate. We examine the constraints as implied by the collider experiments, and then calculate its relic density. We find that this Higgsino-like lightest supersymmetric particle is a viable cold dark matter candidate (0.05 < Omega h^2 < 1), and we discuss its favorable prospects for laboratory detection.

In constrained minimal supersymmetry the lightest neutralino of bino-type is the only neutral candidate for dark matter. As a result, one is typically able to restrict all the supersymmetric mass spectra below roughly 1-2\tev\ {\em without} imposing an ill-defined fine-tuning constraint.

In the constrained MSSM one is typically able to restrict the supersymmetric mass spectra below roughly 1-2\tev\ {\em without} resorting to the ambiguous fine-tuning constraint.