Low-luminosity active galactic nuclei (LL AGNs) have radiatively inefficient accretion flows (RIAFs), where thermal electrons naturally emit not only synchrotron photons but also soft gamma rays via the Comptonization of their synchrotron photons. We find that without any nonthermal electron population, these upscattered photons from LL AGNs naturally account for the MeV gamma-ray background, whose origin has been a mystery. The model also allows proton acceleration via turbulence or reconnections, producing high-energy neutrinos via hadronuclear interactions. We demonstrate that our RIAF model can simultaneously reproduce the soft (MeV) gamma-ray and high-energy (PeV) neutrino backgrounds. The proposed model is consistent with the latest x-ray observations of nearby LL AGNs and testable by future MeV gamma-ray telescopes.

Some short gamma-ray bursts (SGRBs) show a longer lasting emission phase, called extended emission (EE) lasting $\sim 10^{2-3}\,\rm s$, as well as a plateau emission (PE) lasting $\sim10^{4-5}\,\rm s$. While a long-lasting activity of the central engines is a promising explanation for powering both emissions, their physical origin and their emission mechanisms are still uncertain. In this work, we study the properties of the EEs and their connection with the PEs. First, we constrain the minimal Lorentz factor $\Gamma$ of the outflows powering EEs, using compactness arguments and find that the outflows should be relativistic, $\Gamma\gtrsim10$. We propose a consistent scenario for the PEs, where the outflow eventually catches up with the jet responsible for the prompt emission, injecting energy into the forward shock formed by the prior jet, which naturally results in a PE. We also derive the radiation efficiency of EEs and the Lorentz factor of the outflow within our scenario for 10 well-observed SGRBs accompanied by both EE and PE. The efficiency has an average value of $\sim3\,\%$ but shows a broad distribution ranging from $\sim0.01$ to $\sim100\%$. The Lorentz factor is $\sim20-30$, consistent with the compactness arguments. These results suggest that EEs are produced by a slower outflow via more inefficient emission than the faster outflow which causes the prompt emission with a high radiation efficiency.

We calculate spectra of escaping cosmic rays (CRs) accelerated at shocks produced by expanding Galactic superbubbles powered by multiple supernovae producing a continuous energy outflow in star-forming galaxies. We solve the generalized Kompaneets equations adapted to expansion in various external density profiles, including exponential and power-law shapes, and take into account that escaping CRs are dominated by those around their maximum energies. We find that the escaping CR spectrum largely depends on the specific density profiles and power source properties, and the results are compared to and constrained by the observed CR spectrum. As a generic demonstration, we apply the scheme to a superbubble occurring in the centre of the Milky Way, and find that under specific parameter sets the CRs produced in our model can explain the observed CR flux and spectrum around the second knee at $10^{17}$ eV.

We investigate the blazar contribution to the cumulative neutrino intensity assuming a generic relationship between neutrino and gamma-ray luminosities, $L_{\nu}\propto (L_{\rm ph})^{\gamma_{\rm lw}}$. Using the gamma-ray luminosity functions for blazars including flat spectrum radio quasars (FSRQs) and BL Lac objects, as well as the $Fermi$-LAT detection efficiency, we estimate contributions from {blazars resolved by $Fermi$-LAT as well as the unresolved counterpart.} Combining the existing upper limits from stacking analyses, the cumulative neutrino flux from all blazars ({including $Fermi$-LAT resolved and unresolved ones}) are constrained { in the range $0\lesssim\gamma_{\rm lw}\lesssim2.5$}. We also evaluate the effects of the redshift evolution and the effective local number densities for each class of FSRQs, BL Lacs, and all blazars, by which we place another type of constraints on the blazar contribution using the non-detection of high-energy neutrino multiplets. We demonstrate that these two upper limits are complementary, and that the joint consideration of the stacking and multiplet analyses { not only supports the argument that blazars are disfavored as the dominant sources of the 100-TeV neutrino background, but it extends this argument by including also $Fermi$-LAT-unresolved blazars as well, for a more generic luminosity correlation $L_{\nu}\propto (L_{\rm ph})^{\gamma_{\rm lw}}$}.

We investigate prolonged engine activities of short gamma-ray bursts (SGRBs), such as extended and/or plateau emissions, as high-energy gamma-ray counterparts to gravitational waves (GWs). Binary neutron-star mergers lead to relativistic jets and merger ejecta with $r$-process nucleosynthesis, which are observed as SGRBs and kilonovae/macronovae, respectively. Long-term relativistic jets may be launched by the merger remnant as hinted in X-ray light curves of some SGRBs. The prolonged jets may dissipate their kinetic energy within the radius of the cocoon formed by the jet-ejecta interaction. Then the cocoon supplies seed photons to non-thermal electrons accelerated at the dissipation region, causing high-energy gamma-ray production through the inverse Compton scattering process. We numerically calculate high-energy gamma-ray spectra in such a system using a one-zone and steady-state approximation, and show that GeV--TeV gamma-rays are produced with a duration of $10^2-10^5$ seconds. They can be detected by {\it Fermi}/LAT or CTA as gamma-ray counterparts to GWs.

The cores of active galactic nuclei (AGNs) have been suggested as the sources of IceCube neutrinos, and recent numerical simulations have indicated that hot AGN coronae of Seyfert galaxies and radiatively inefficient accretion flows (RIAFs) of low-luminosity AGNs (LLAGNs) may be promising sites of ion acceleration. We present detailed studies on detection prospects of high-energy multi-messenger emissions from RIAFs in nearby LLAGNs. We construct a model of RIAFs that can reproduce the observational features of the current X-ray observations of nearby LLAGNs. We then calculate the high-energy particle emissions from nearby individual LLAGNs, including MeV gamma rays from thermal electrons, TeV--PeV neutrinos produced by non-thermal protons, and sub-GeV to sub-TeV gamma rays from proton-induced electromagnetic cascades. We find that, although these are beyond the reach of current facilities, proposed future experiments such as e-ASTROGAM and IceCube-Gen2 should be able to detect the MeV gamma rays and the neutrinos, respectively, or else they can place meaningful constraints on the parameter space of the model. On the other hand, the detection of high-energy gamma rays due to the electromagnetic cascades will be challenging with the current and near-future experiments, such as Fermi and Cherenkov Telescope Array. In an accompanying paper, we demonstrate that LLAGNs can be a source of the diffuse soft gamma-ray and TeV--PeV neutrino backgrounds, whereas in the present paper, we focus on the prospects for multi-messenger tests which can be applied to reveal the nature of the high-energy neutrinos and photons from LLAGNs.

Multi-messenger astrophysics, a long-anticipated extension to traditional and multiwavelength astronomy, has recently emerged as a distinct discipline providing unique and valuable insights into the properties and processes of the physical universe. These insights arise from the inherently complementary information carried by photons, gravitational waves, neutrinos, and cosmic rays about individual cosmic sources and source populations. Realizing the observation of astrophysical sources via non-photonic messengers has presented enormous challenges, as evidenced by the fiscal and physical scales of the multi-messenger observatories. However, the scientific payoff has already been substantial, with even greater rewards promised in the years ahead. In this review we survey the current status of multi-messenger astrophysics, highlighting some exciting recent results, and addressing the major follow-on questions they have raised. Key recent achievements include the measurement of the spectrum of ultra-high energy cosmic rays out to the highest observable energies; discovery of the diffuse high energy neutrino background; the first direct detections of gravitational waves and the use of gravitational waves to characterize merging black holes and neutron stars in strong-field gravity; and the identification of the first joint electromagnetic + gravitational wave and electromagnetic + high-energy neutrino multi-messenger sources. We then review the rationales for the next generation of multi-messenger observatories, and outline a vision of the most likely future directions for this exciting and rapidly advancing field.

Fast radio bursts (FRBs) are short duration ($\sim$millisecond) radio transients with cosmological origin. The simple sharp features of the FRB signal have been utilized to probe two fundamental laws of physics, namey, testing Einstein's weak equivalence principle and constraining the rest mass of the photon. Recently, \cite{hessels18} found that after correcting for dispersive delay, some of the bursts in FRB 121102 have complex time-frequency structures that include sub-pulses with a time-frequency downward drifting property. Using the delay time between sub-pulses in FRB 121102, here we show that the parameterized post-Newtonian parameter $\gamma$ is the same for photons with different energies to the level of $\left|\gamma_1-\gamma_2\right|<2.5\times 10^{-16}$, which is 1000 times better than previous constraints from FRBs using similar methods. We also obtain a stringent constraint on the photon mass, $m_{\gamma} < 5.1\times10^{-48}$ g, which is 10 times smaller than previous best limits on the photon mass derived through the velocity dispersion method.

Shocks arising in galaxy mergers could accelerate cosmic-ray (CR) ions to TeV-PeV energies. While propagating in the intergalactic medium, these CRs can produce high-energy neutrinos, electron-positron pairs and gamma-rays. In the presence of intergalactic magnetic fields, the secondary pairs will radiate observable emissions through synchrotron radiation and inverse Compton scattering. In this paper, we demonstrate that these emissions can explain the radio and X-ray fluxes of merging galaxies such as NGC 660 and NGC 3256. Using our model in combination with the observations, we can constrain the gas mass, shock velocity, magnetic field and the CR spectral index $s$ of these systems. For NGC 660 a single-zone model with a spectral index $2.1\lesssim s\lesssim2.2$ is able to reproduce simultaneously the radio and X-ray observations, while a simple one-zone scenario with $s\sim2$ can describe the radio and a large fraction of X-ray observations of NGC 3256. Our work provides a useful approach for studying the dynamics and physical parameters of galaxy mergers, which can play an important part in future multi-messenger studies of similar and related extragalactic sources.

Motivated by the recently reported evidence of an association between a high-energy neutrino and a gamma-ray flare from the blazar TXS 0506+056, we calculate the expected high-energy neutrino signal from past, individual flares, from twelve blazars, selected in declinations favourable for detection with IceCube. To keep the number of free parameters to a minimum, we mainly focus on BL Lac objects and assume the synchrotron self-Compton mechanism produces the bulk of the high-energy emission. We consider a broad range of the allowed parameter space for the efficiency of proton acceleration, the proton content of BL Lac jets, and the presence of external photon fields. To model the expected neutrino fluence we use simultaneous multi-wavelength observations. We find that in the absence of external photon fields and with jet proton luminosity normalised to match the observed production rate of ultra-high-energy cosmic rays, individual flaring sources produce a modest neutrino flux in IceCube, $\lesssim10^{-3}$ muon neutrinos with energy exceeding 100 TeV, stacking ten years of flare periods selected in the >800 MeV Fermi energy range from each source. Under optimistic assumptions about the jet proton luminosity and in the presence of external photon fields, we find that the two most powerful sources in our sample, AO 0235+164, and OJ 287, would produce, in total, $\approx 3$ muon neutrinos during ten years of Fermi flaring periods, in future neutrino detectors with total instrumented volume $\sim$ten times larger than IceCube,or otherwise, constrain the proton luminosity of blazar jets.

I discuss some aspects of the evolution of the standard GRB model, emphasizing various theoretical developments in the last decade, and review the impact of some of the most recent observational discoveries and the new challenges they pose in the expanding realm of multi-messenger astrophysics.

The ANITA collaboration have reported observation of two anomalous events that appear to be $\varepsilon_{\rm cr} \approx 0.6$ EeV cosmic ray showers emerging from the Earth with exit angles of $27^\circ$ and $35^\circ$, respectively. While EeV-scale upgoing showers have been anticipated as a result of astrophysical tau neutrinos converting to tau leptons during Earth passage, the observed exit angles are much steeper than expected in Standard Model (SM) scenarios. Indeed, under conservative extrapolations of the SM interactions, there is no particle that can propagate through the Earth with probability $p > 10^{-6}$ at these energies and exit angles. We explore here whether "beyond the Standard Model" (BSM) particles are required to explain the ANITA events, if correctly interpreted, and conclude that they are. Seeking confirmation or refutation of the physical phenomenon of sub-EeV Earth-emergent cosmic rays in data from other facilities, we find support for the reality of the ANITA events, and three candidate analog events, among the Extremely High Energy Northern Track neutrinos of the IceCube Neutrino Observatory. Properties of the implied BSM particle are anticipated, at least in part, by those predicted for the "stau" slepton ($\tilde{\tau}_R$) in some supersymmetric models of the fundamental interactions, wherein the stau manifests as the next-to-lowest mass supersymmetric partner particle.

The detection of gravitational waves and electromagnetic counterparts from a binary neutron star (BNS) merger confirmed that it is accompanied by the launch of fast merger ejecta. Analogous to supernova remnants, forward shocks formed by the interaction of the ejecta with interstellar material will produce high-energy cosmic rays. We investigate the possibility that Galactic neutron star merger remnants (NSMRs) significantly contribute to the observed cosmic rays in the energy range between the knee and the ankle. Using typical parameters obtained by modeling of GW170817, we find that NSMRs can accelerate iron nuclei up to $\sim500$ PeV. We calculate the cosmic-ray spectrum and composition observed on Earth, and show that the Galactic NSMR scenario can account for the experimental cosmic-ray data in the 20 -- 1000 PeV range. Our model can naturally explain the hardening feature around 20 PeV for total cosmic-ray spectrum, which has been observed by the Telescope Array Low Energy extension and the IceTop air shower array.

The observations of a macronova/kilonova accompanied by gravitational waves from a binary neutron star merger (GW170817) confirmed that neutron star coalescences produce copious ejecta. The coincident gamma-ray detection implies the existence of a relativistic jet in this system. During the jet's propagation within the ejecta, high-energy photons and neutrinos can be produced. The photons are absorbed by the ejecta, while the neutrinos escape and can be detected. Here, we estimate such trans-ejecta neutrino emission, and discuss how neutrino observations could be used to differentiate between gamma-ray burst emission scenarios. We find that neutrinos from the internal shocks inside the ejecta may be detectable by IceCube within a few years of operation, and will likely be detected with IceCube-Gen2. The neutrino signals coincident with gravitational waves would enable us to reveal the physical quantities of the choked jets even without electromagnetic signals.

Recent results from the Pierre Auger Collaboration have shown that the composition of ultrahigh-energy cosmic rays (UHECRs) becomes gradually heavier with increasing energy. Although gamma-ray bursts (GRBs) have been promising sources of UHECRs, it is still unclear whether they can account for the Auger results because of their unknown nuclear composition of ejected UHECRs. In this work, we revisit the possibility that low-luminosity GRBs (LL GRBs) act as the sources of UHECR nuclei, and give new predictions based on the intrajet nuclear composition models considering progenitor dependencies. We find that the nuclear component in the jet can be divided into two groups according to the mass fraction of silicon nuclei, Si-free and Si-rich. Motivated by the connection between LL GRBs and transrelativistic supernovae, we also consider the hypernova ejecta composition. Then, we discuss the survivability of UHECR nuclei in the jet base and internal shocks of the jets, and show that it is easier for nuclei to survive for typical LL GRBs. Finally, we numerically propagate UHECR nuclei ejected from LL GRBs with different composition models and compare the resulting spectra and composition to Auger data. Our results show that both the Si-rich progenitor and hypernova ejecta models match the Auger data well, while the Si-free progenitor models have more difficulty in fitting the spectrum. We argue that our model is consistent with the newly reported cross correlation between the UHECRs and starburst galaxies, since both LL GRBs and hypernovae are expected to be tracers of the star-formation activity. LL GRBs have also been suggested as the dominant origin of IceCube neutrinos in the PeV range, and the LL GRB origin of UHECRs can be critically tested by near-future multimessenger observations.

The merger of dark matter halos and the gaseous structures embedded in them, such as proto-galaxies, galaxies, and groups and clusters of galaxies, results in strong shocks that are capable of accelerating cosmic rays (CRs) to $\sim10~\rm PeV$. These shocks will produce high-energy neutrinos and $\gamma$-rays through inelastic $pp$ collisions with ambient gaseous environments. In this work, we study the contributions of these halo mergers to the diffuse neutrino flux measured in IceCube and to the non-blazar portion of the extragalactic $\gamma$-ray background measured by $Fermi$. In order to calculate them, we formulate the redshift dependence of the shock velocity, galactic radius, halo gas content and galactic/intergalactic magnetic fields over the dark matter halo distribution up to a redshift $z=10$. We find that high-redshift mergers contribute a significant amount of the cosmic-ray energy luminosity density, and the resulting neutrino spectra could explain a large part of the observed diffuse neutrino flux above 0.1 PeV up to several PeV. We also show that our model can somewhat alleviate tensions with the extragalactic $\gamma$-ray background. First, since a larger fraction of the CR energy luminosity density comes from high redshifts, the accompanying $\gamma$-rays are more strongly suppressed through $\gamma\gamma$ annihilations with the cosmic microwave background (CMB) and the extragalactic background light (EBL). Second, mildly radiative-cooled shocks may lead to a harder CR spectrum with spectral indices of $1.5\lesssim s\lesssim2.0$. Our study suggests that halo mergers, a fraction of which may also induce starbursts in the merged galaxies, can be promising neutrino emitters without violating the existing $Fermi$ $\gamma$-ray constraints on the non-blazar component of the extragalactic $\gamma$-ray background.

For the first time, a short gamma-ray burst (GRB) was unambiguously associated with a gravitational wave (GW) observation from a binary neutron star (NS) merger. This allows us to link the details of the central engine properties to GRB emission models. We find that photospheric models (both dissipative and non-dissipative variants) have difficulties accounting for the observations. Internal shocks give the most natural account of the observed peak energy, viewing angle and total energy. We also show that a simple external shock model can reproduce the observed GRB pulse with parameters consistent with those derived from the afterglow modeling. We find a simple cocoon shock breakout model is in mild tension with the observed spectral evolution, however it cannot be excluded based on gamma-ray data alone. Future joint observations of brighter GRBs will pose even tighter constraints on prompt emission models.

There are indications that $\gamma$-ray dark objects such as supernovae (SNe) with choked jets, and the cores of active galactic nuclei may contribute to the diffuse flux of astrophysical neutrinos measured by the IceCube observatory. In particular, stripped-envelope SNe have received much attention since they are capable of producing relativistic jets and could explain the diversity in observations of collapsar explosions (e.g., gamma-ray bursts (GRBs), low-luminosity GRBs, and Type Ibc SNe). We use an unbinned maximum likelihood method to search for spatial and temporal coincidences between Type Ibc core-collapse SNe, which may harbor a choked jet, and muon neutrinos from a sample of IceCube up-going track-like events measured from May 2011-May 2012. In this stacking analysis, we find no significant deviation from a background-only hypothesis using one year of data, and are able to place upper limits on the total amount of isotropic equivalent energy that choked jet core-collapse SNe deposit in cosmic rays ${\mathcal E}_{\rm cr}$ and the fraction of core-collapse SNe which have a jet pointed towards Earth $f_{\rm jet}$. This analysis can be extended with yet to be made public IceCube data, and the increased amount of optically detected core-collapse SNe discovered by wide field-of-view surveys such as the Palomar Transient Factory and All-Sky Automated Survey for Supernovae. The choked jet SNe/high-energy cosmic neutrino connection can be more tightly constrained in the near future.

On February 17 2016, the IceCube real-time neutrino search identified, for the first time, three muon neutrino candidates arriving within 100 s of one another, consistent with coming from the same point in the sky. Such a triplet is expected once every 13.7 years as a random coincidence of background events. However, considering the lifetime of the follow-up program the probability of detecting at least one triplet from atmospheric background is 32%. Follow-up observatories were notified in order to search for an electromagnetic counterpart. Observations were obtained by Swift's X-ray telescope, by ASAS-SN, LCO and MASTER at optical wavelengths, and by VERITAS in the very-high-energy gamma-ray regime. Moreover, the Swift BAT serendipitously observed the location 100 s after the first neutrino was detected, and data from the Fermi LAT and HAWC observatory were analyzed. We present details of the neutrino triplet and the follow-up observations. No likely electromagnetic counterpart was detected, and we discuss the implications of these constraints on candidate neutrino sources such as gamma-ray bursts, core-collapse supernovae and active galactic nucleus flares. This study illustrates the potential of and challenges for future follow-up campaigns.

The coincident detection of a gravitational-wave (GW) event GW170817 with electromagnetic (EM) signals (e.g., a short gamma-ray burst SGRB 170817A or a macronova) from a binary neutron star merger within the nearby galaxy NGC 4933 provides a new, multimessenger test of the weak equivalence principle (WEP), extending the WEP test with GWs and photons. Assuming that the arrival time delay between the GW signals from GW170817 and the photons from SGRB 170817A or the macronova is mainly attributed to the gravitational potential of the Milky Way, we demonstrate that the strict upper limits on the deviation from the WEP are $\Delta \gamma<1.4\times10^{-3}$ for GW170817/macronova and $\Delta \gamma <5.9\times10^{-8}$ for GW170817/SGRB 170817A. A much more severe constraint on the WEP accuracy can be achieved ($\sim0.9\times10^{-10}$) for GW170817/SGRB 170817A when we consider the gravitational potential of the Virgo Cluster, rather than the Milky Way's gravity. This provides the tightest limit to date on the WEP through the relative differential variations of the $\gamma$ parameter for two different species of particles. Compared with other multimessenger (photons and neutrinos) results, our limit is 7 orders of magnitude tighter than that placed by the neutrinos and photons from supernova 1987A, and is almost as good as or is an improvement of 6 orders of magnitude over the limits obtained by the low-significance neutrinos correlated with GRBs and a blazar flare.

We discuss the electromagnetic radiation from newborn binary black holes (BBHs). As a consequence of the evolution of massive stellar binaries, a binary consisting of a primary black hole (BH) and a secondary Wolf-Rayet star is expected as a BBH progenitor system. We investigate optical transients from the birth of BBHs powered by the Bondi-Hoyle-Lyttleton accretion onto the primary BH, which occur $\sim1-10$ Gyr earlier than gravitational wave signals at the BH-BH merger. When the secondary massive star collapses into a BH, it may eject a fraction of its outer material and may form a disk around the primary BH and induces a powerful disk wind. These primary-induced winds can lead to optical transients with a kinetic energy of $\sim10^{47}-3\times10^{48}$ erg, an ejecta velocity of $10^8-10^9\rm~cm~s^{-1}$, a duration of a few days, and an absolute magnitude ranging from about $-11$ to $-14$. The light curves and late-time spectra of these transients are distinctive from those of ordinary supernovae, and detection of this type of transient is possible by future optical transient surveys if the event rate of this transient is comparable to the merger rate of BBHs. This paper focuses on the emissions from disk-driven transients induced by the primary BH, different from the first paper that focuses on wind-driven transients from the tidally-locked secondary massive star.

We propose a new type of electromagnetic transients associated with the birth of binary black holes (BBHs), which may lead to merger events accompanied by gravitational waves in $\sim0.1-1$ Gyr. We consider the newborn BBHs formed through the evolution of isolated massive stellar binaries. For a close massive binary, consisting of a primary black hole (BH) and a secondary Wolf-Rayet (WR) star that are orbiting around each other, the spin period of the WR star can be tidally synchronized to its orbital period. Then, the angular momentum of the outer material of the WR star is large enough to form an accretion disk around a newborn, secondary BH, following its core-collapse. This disk may produce an energetic outflow with a kinetic energy of $\sim10^{50}-10^{52}$ erg and an outflow velocity of $\sim10^{10}\rm~cm~s^{-1}$, resulting in an optical transient with an absolute magnitude from $\sim -14$ to $\sim-17$ with a duration of around a day. This type of transient also produces detectable radio signals $\sim1-10$ years after the birth of BBHs, via synchrotron emission from non-thermal electrons accelerated at external shocks. The predicted optical transients have a shorter duration than ordinary core-collapse supernovae. Dedicated optical transient surveys can detect them, and distinguish it from ordinary SNe using the different features of its light curve and late-time spectrum. In the paper I, we investigate disk-driven outflows from the secondary BH, whereas possible signals from the primary BH will be discussed in the paper II.

We investigate current and future prospects for coincident detection of high-energy neutrinos and gravitational waves (GWs). Short gamma-ray bursts (SGRBs) are believed to originate from mergers of compact star binaries involving neutron stars. We estimate high-energy neutrino fluences from prompt emission, extended emission, X-ray flares, and plateau emission, and show that neutrino signals associated with the extended emission are the most promising. Assuming that the cosmic-ray loading factor is $\sim10$ and the Lorentz factor distribution is lognormal, we calculate the probability of neutrino detection from extended emission by current and future neutrino detectors, and find that the quasi-simultaneous detection of high-energy neutrinos, gamma rays, and GWs is possible with future instruments or even with current instruments for nearby SGRBs having extended emission. We also discuss stacking analyses that will also be useful with future experiments such as IceCube-Gen2.

Violations of Lorentz invariance can lead to an energy-dependent vacuum dispersion of light, which results in arrival-time differences of photons arising with different energies from a given transient source. In this work, direction-dependent dispersion constraints are obtained on nonbirefringent Lorentz-violating effects, using the observed spectral lags of the gamma-ray burst GRB 160625B. This burst has unusually large high-energy photon statistics, so we can obtain constraints from the true spectral time lags of bunches of high-energy photons rather than from the rough time lag of a single highest-energy photon. Also, GRB 160625B is the only burst to date having a well-defined transition from positive lags to negative lags, which provides a unique opportunity to distinguish Lorentz-violating effects from any source-intrinsic time lag in the emission of photons of different energy bands. Our results place comparatively robust two-sided constraints on a variety of isotropic and anisotropic coefficients for Lorentz violation, including first bounds on Lorentz-violating effects from operators of mass dimension ten in the photon sector.

Einstein's weak equivalence principle (WEP) states that any freely falling, uncharged test particle follows the same identical trajectory independent of its internal structure and composition. Since the polarization of a photon is considered to be part of its internal structure, we propose that polarized photons from astrophysical transients, such as gamma-ray bursts (GRBs) and fast radio bursts (FRBs), can be used to constrain the accuracy of the WEP through the Shapiro time delay effect. Assuming that the arrival time delays of photons with different polarizations are mainly attributed to the gravitational potential of the Laniakea supercluster of galaxies, we show that a strict upper limit on the differences of the parametrized post-Newtonian parameter $\gamma$ value for the polarized optical emission of GRB 120308A is $\Delta\gamma<1.2\times10^{-10}$, for the polarized gamma-ray emission of GRB 100826A is $\Delta\gamma<1.2\times10^{-10}$, and for the polarized radio emission of FRB 150807 is $\Delta\gamma<2.2\times10^{-16}$. These are the first direct verifications of the WEP for multiband photons with different polarizations. In particular, the result from FRB 150807 provides the most stringent limit to date on a deviation from the WEP, improving by one order of magnitude the previous best result based on Crab pulsar photons with different energies.

We discuss likely sources of cosmic rays in the $10^{15}-10^{20}$ eV range and their possible very high energy neutrino and gamma-ray signatures which could serve to identify these sources and constrain their physics. Among these sources we discuss in particular low luminosity gamma-ray bursts, including choked and shock-breakout objects, tidal disruption events and white dwarf mergers. Among efforts aimed at simultaneous secondary multi-messenger detections we discuss the AMON program.

We review the physics of GRB production by relativistic jets that start highly opaque near the central source and then expand to transparency. We discuss dissipative and radiative processes in the jet and how radiative transfer shapes the observed nonthermal spectrum released at the photosphere. A comparison of recent detailed models with observations gives estimates for important parameters of GRB jets, such as the Lorentz factor and magnetization. We also discuss predictions for GRB polarization and neutrino emission.

Possible violations of Lorentz invariance (LIV) have been investigated for a long time using the observed spectral lags of gamma-ray bursts (GRBs). However, these generally have relied on using a single photon in the highest energy range. Furthermore, the search for LIV lags has been hindered by our ignorance concerning the intrinsic time lag in different energy bands. GRB 160625B, the only burst so far with a well-defined transition from $positive$ lags to $negative$ lags provides a unique opportunity to put new constraints on LIV. Using multi-photon energy bands we consider the contributions to the observed spectral lag from both the intrinsic time lag and the lag by LIV effects, and assuming the intrinsic time lag to have a positive dependence on the photon energy, we obtain robust limits on LIV by directly fitting the spectral lag data of GRB 160625B. Here we show that these robust limits on the quantum gravity energy scales are $E_{\rm QG,1}\geq0.5\times10^{16}$ GeV for the linear, and $E_{\rm QG,2}\geq1.4\times10^{7}$ GeV for the quadratic LIV effects, respectively. In addition, we give for the first time a reasonable formulation of the intrinsic energy-dependent time lag.

We report our discovery in Swift satellite data of a transient gamma-ray counterpart (3.2$\sigma$ confidence) to the fast radio burst FRB131104, the first such counterpart to any FRB. The transient has duration $T_{90} \gtrsim 100$s and fluence $S_\gamma\approx 4\times 10^{-6}$ erg cm$^{-2}$, increasing the energy budget for this event by more than a billion times; at the nominal $z\approx 0.55$ redshift implied by its dispersion measure, the burst's gamma-ray energy output is $E_\gamma \approx 5\times 10^{51}$ erg. The observed radio to gamma-ray fluence ratio for FRB131104 is consistent with a lower limit we derive from Swift observations of another FRB, which is not detected in gamma-rays, and with an upper limit previously derived for the brightest gamma-ray flare from SGR 1806-20, which was not detected in the radio. X-ray, ultraviolet, and optical observations beginning two days after the FRB do not reveal any associated afterglow, supernova, or transient; Swift observations exclude association with the brightest 65% of Swift gamma-ray burst X-ray afterglows, while leaving the possibility of an associated supernova at much more than 10% the FRB's nominal distance, $D\gtrsim 320$ Mpc, largely unconstrained. Transient high-luminosity gamma-ray emission arises most naturally in a relativistic outflow or shock breakout, as for example from magnetar flares, gamma-ray bursts, relativistic supernovae, and some types of galactic nuclear activity. Our discovery thus bolsters the case for an extragalactic origin for some FRBs and suggests that future rapid-response observations might identify long-lived counterparts, resolving the nature of these mysterious phenomena and realizing their promise as probes of cosmology and fundamental physics.

The merger of two white dwarfs is expected to result in a central fast rotating core surrounded by a debris disk, in which magnetorotational instabilities give rise to a hot magnetized corona and a magnetized outflow. The dissipation of magnetic energy via reconnection could lead to the acceleration of cosmic-rays in the expanding material, which would result in high energy neutrinos. We discuss the possibility of using these neutrino signals as probes of the outflow dynamics, magnetic energy dissipation rate and cosmic-ray acceleration efficiency. Importantly, the accompanying high-energy gamma-rays are absorbed within these sources because of the large optical depth, so these neutrino sources can be regarded as hidden cosmic-ray accelerators that are consistent with the non-detection of gamma-rays with Fermi-LAT. While the cosmic-ray generation rate is highly uncertain, if it reaches $\sim10^{45}\,\rm erg\,Mpc^{-3}\,yr^{-1}$, the diffuse neutrino flux could contribute a substantial fraction of the IceCube observations. We also evaluate the prospect of observing individual merger events, which provides a means for testing such sources in the future.

Five TeV neutrino events weakly correlated with five gamma-ray bursts (GRBs) were detected recently by IceCube. This work is an attempt to show that if the GRB identifications are verified, the observed time delays between the TeV neutrinos and gamma-ray photons from GRBs provide attractive candidates for testing fundamental physics with high accuracy. Based on the assumed associations between the TeV neutrinos and GRBs, we find that the limiting velocity of the neutrinos is equal to that of photons to an accuracy of $\sim1.9\times10^{-15}-2.5\times10^{-18}$, which is about $10^{4}-10^{7}$ times better than the constraint obtained with the neutrino possibly from a blazar flare. In addition, we set the most stringent limits up to date on the energy scale of quantum gravity for both the linear and quadratic violations of Lorentz invariance, namely $E_{\rm QG, 1}>6.3\times10^{18}-1.5\times10^{21}$ GeV and $E_{\rm QG, 2}>2.0\times10^{11}-4.2\times10^{12}$ GeV, which are essentially as good as or are an improvement of one order of magnitude over the results previously obtained by the GeV photons of GRB 090510 and the PeV neutrino from a blazar flare. Assuming that the Shapiro time delay is caused by the gravitational potential of the Laniakea supercluster of galaxies, we also place the tightest limits to date on Einstein's weak equivalence principle through the relative differential variations of the parameterized post-Newtonian parameter $\gamma$ values for two different species of particles (i.e., neutrinos and photons), yielding $\Delta\gamma \sim 10^{-11}-10^{-13}$. However, it should be emphasized again that these limits here obtained are at best forecast of what could be achieved if the GRB/neutrino correlations would be finally confirmed.

The possible short gamma-ray burst (GRB) observed by {\it Fermi}/GBM in coincidence with the first gravitational wave (GW) detection, offers new ways to test GRB prompt emission models. Gravitational wave observations provide previously unaccessible physical parameters for the black hole central engine such as its horizon radius and rotation parameter. Using a minimum jet launching radius from the Advanced LIGO measurement of GW~150914, we calculate photospheric and internal shock models and find that they are marginally inconsistent with the GBM data, but cannot be definitely ruled out. Dissipative photosphere models, however have no problem explaining the observations. Based on the peak energy and the observed flux, we find that the external shock model gives a natural explanation, suggesting a low interstellar density ($\sim 10^{-3}$ cm$^{-3}$) and a high Lorentz factor ($\sim 2000$). We only speculate on the exact nature of the system producing the gamma-rays, and study the parameter space of a generic Blandford Znajek model. If future joint observations confirm the GW-short GRB association we can provide similar but more detailed tests for prompt emission models.

A conservative constraint on the Einstein Weak Equivalence Principle (WEP) can be obtained under the assumption that the observed time delay between correlated particles from astronomical sources is dominated by the gravitational fields through which they move. Current limits on the WEP are mainly based on the observed time delays of photons with different energies. It is highly desirable to develop more accurate tests that include the gravitational wave (GW) sector. The detection by the advanced LIGO/VIRGO systems of gravitational waves will provide attractive candidates for constraining the WEP, extending the tests to gravitational interactions, with potentially higher accuracy. Considering the capabilities of the advanced LIGO/VIRGO network and the source direction uncertainty, we show that the joint detection of GWs and electromagnetic signals could probe the WEP to an accuracy down to $10^{-10}$, which is one order of magnitude tighter than previous limits, and seven orders of magnitude tighter than the multi-messenger (photons and neutrinos) results by supernova 1987A.

Star-forming and starburst galaxies are considered as one of the viable candidate sources of the high-energy cosmic neutrino background detected in IceCube. We revisit contributions of supernova remnants (SNRs) and hypernova remnants (HNRs) in such galaxies to the diffuse high-energy neutrino and gamma-ray backgrounds, in light of the latest Fermi data above 50GeV. We also take into account possible time dependent effects of the cosmic-ray (CR) acceleration during the SNR evolution. CRs accelerated by the SNR shocks can produce high-energy neutrinos up to $\sim100$ TeV energies, but CRs from HNRs can extend the spectrum up to PeV energies. We show that, only if HNRs are dominant over SNRs, the diffuse neutrino background above 100 TeV can be explained without contradicting the gamma-ray data. However, the neutrino data around 30 TeV remain unexplained, which might suggest a different population of gamma-ray dark CR sources. Alternatively, we consider possible contributions of Pop-III HNRs up to $z\lesssim10$, and show that they are not constrained by the gamma-ray data, and thus could contribute to the diffuse high-energy backgrounds if their explosion energy reaches ${\mathcal E}_{\rm POP-III}\sim{\rm a~few}\times10^{53}$erg. More conservatively, our results suggest that the explosion energy of POP-III HNRs is ${\mathcal E}_{\rm POP-III}\lesssim7\times{10}^{53}$erg.

Fast radio bursts (FRBs) are radio bursts characterized by millisecond durations, high Galactic latitude positions, and high dispersion measures. Very recently, the cosmological origin of FRB 150418 has been confirmed by \cite{kea16}, and FRBs are now strong competitors as cosmological probes. The simple sharp feature of the FRB signal is ideal for them to probe some of the fundamental laws of physics. Here we show that by analyzing the delay time between different frequencies, the FRB data can place stringent upper limits on the rest mass of the photon. For FRB 150418 at $z=0.492$, one can potentially reach $m_{\gamma} \leq 5.2 \times 10^{-47}$ g, which is $10^{20}$ times smaller than the rest mass of electron, and is about $10^3$ times smaller than that obtained using other astrophysical sources with the same method.

The accuracy of Einstein's Equivalence Principle (EEP) can be tested with the observed time delays between correlated particles or photons that are emitted from astronomical sources. Assuming as a lower limit that the time delays are caused mainly by the gravitational potential of the Milky Way, we prove that fast radio bursts (FRBs) of cosmological origin can be used to constrain the EEP with high accuracy. Taking FRB 110220 and two possible FRB/gamma-ray burst (GRB) association systems (FRB/GRB 101011A and FRB/GRB 100704A) as examples, we obtain a strict upper limit on the differences of the parametrized post-Newtonian parameter $\gamma$ values as low as $\left[\gamma(1.23\; \rm GHz)-\gamma(1.45\; \rm GHz)\right]<4.36\times10^{-9}$. This provides the most stringent limit up to date on the EEP through the relative differential variations of the $\gamma$ parameter at radio energies, improving by 1 to 2 orders of magnitude the previous results at other energies based on supernova 1987A and GRBs.

Gamma-ray burst sources appear to fulfill all the conditions for being efficient cosmic ray accelerators, and being extremely compact, are also expected to produce multi-GeV to PeV neutrinos. I review the basic model predictions for the expected neutrino fluxes in classical GRBs as well as in low luminosity and choked bursts, discussing the recent IceCube observational constraints and implications from the observed diffuse neutrino flux.

We have performed a search for muon neutrinos from dark matter annihilation in the center of the Sun with the 79-string configuration of the IceCube neutrino telescope. For the first time, the DeepCore sub-array is included in the analysis, lowering the energy threshold and extending the search to the austral summer. The 317 days of data collected between June 2010 and May 2011 are consistent with the expected background from atmospheric muons and neutrinos. Upper limits are therefore set on the dark matter annihilation rate, with conversions to limits on spin-dependent and spin-independent WIMP-proton cross-sections for WIMP masses in the range 20 - 5000 GeV. These are the most stringent spin-dependent WIMP-proton cross-sections limits to date above 35 GeV.

The Einstein Equivalence Principle (EEP) can be probed with astrophysical sources emitting simultaneously different types of neutral particles, or particles with varying energies, by testing their time of flight through the same gravitational field. Here we use the time delays between correlated photons from cosmological transients to constrain the accuracy of the EEP. We take data from two gamma-ray bursts as an example, and use, as a lower limit to the theoretical time delays between different energies, delays arising from only the gravitational field of our own galaxy. We then show that the parameterized post-Newtonian parameter $\gamma$ is the same for photons over energy ranges between eV and MeV and between MeV and GeV to a part in $10^{-7}$, which is at least one order of magnitude better than previous limits. Combining this bound on the wavelength dependence of $\gamma$ with the absolute bound $|\gamma-1|<0.3\%$ from light-deflection measurements at optical (eV) wavelengths, we thus extend this absolute bound on $\gamma$ to GeV energies.

A generic synchrotron external shock model is the widely preferred paradigm used to interpret the broad-band afterglow data of gamma-ray bursts (GRBs), including predicted observable signatures from a reverse shock which have been confirmed by observations. Investigations of the nature of the reverse shock emission can provide valuable insights into the intrinsic properties of the GRB ejecta. Here we briefly review the standard and the extended models of the reverse shock emission, discussing the connection between the theory and observations, including the implications of the latest observational advances.

Within the framework of the external shock model of gamma-ray bursts (GRBs) afterglows, we perform a morphological analysis of the early optical lightcurves to directly constrain model parameters. We define four morphological types, i.e. the reverse shock dominated cases with/without the emergence of the forward shock peak (Type I/ Type II), and the forward shock dominated cases without/with $\nu_m$ crossing the band (Type III/IV). We systematically investigate all the Swift GRBs that have optical detection earlier than 500 s and find 3/63 Type I bursts (4.8%), 12/63 Type II bursts (19.0%), 30/63 Type III bursts (47.6%), 8/63 Type IV bursts (12.7%) and 10/63 Type III/IV bursts (15.9%). We perform Monte Carlo simulations to constrain model parameters in order to reproduce the observations. We find that the favored value of the magnetic equipartition parameter in the forward shock ($\epsilon_B^f$) ranges from $10^{-6}$ to $10^{-2}$, and the reverse-to-forward ratio of $\epsilon_B$ ($R_B$) is about 100. The preferred electron equipartition parameter $\epsilon_e^{r,f}$ value is 0.01, which is smaller than the commonly assumed value, e.g., 0.1. This could mitigate the so- called "efficiency problem" for the internal shock model, if $\epsilon_e$ during the prompt emission phase (in the internal shocks) is large (say, $\sim 0.1$). The preferred $R_B$ value is in agreement with the results in previous works that indicates a moderately magnetized baryonic jet for GRBs.

Gamma-Ray Bursts are the most energetic explosions in the Universe, and are among the most promising for detecting multiple non-electromagnetic signals, including cosmic rays, high energy neutrinos and gravitational waves. The multi-GeV to TeV gamma-ray range of GRB could have significant contributions from hadronic interactions, mixed with more conventional leptonic contributions. This energy range is important for probing the source physics, including overall energetics, the shock parameters and the Lorentz factor. We discuss some of the latest observational and theoretical developments in the field.

In light of the latest IceCube data, we discuss the implications of the cosmic ray energy input from hypernovae and supernovae into the Universe, and their propagation in the hosting galaxy and galaxy clusters or groups. The magnetic confinement in these environments may lead to efficient $pp$ collisions, resulting in a diffuse neutrino spectrum extending from PeV down to 10 TeV energies, with a spectrum and flux level compatible with that recently reported by IceCube. If the diffuse 10 TeV neutrino background largely comes from such the CR reservoirs, the corresponding diffuse gamma-ray background should be compatible with the recent \textit{Fermi} data. In this scenario, the CR energy input from hypernovae should be dominant over that of supernovae, implying that the starburst scenario does not work if the supernova energy budget is a factor of two larger than the hypernova energy budget. Thus, this strong case scenario can be supported or ruled out in near future.

The GRB central engine intrinsic activity time $T_{\rm ce}$ is usually described through either the $\gamma$-ray duration $T_{90}$ or through a generalized burst duration $t_{\rm burst}$ which includes both the $\gamma$-ray emission and (when present) an extended flaring X-ray plateau. Here, we define a more specific operational description of $T_{\rm ce}$, and within the framework of the internal-external shock model, we develop a numerical code to study the relationship between $T_{90}$ and $T_{\rm ce}$, as well as between $t_{\rm burst}$ and $T_{\rm ce}$, for different initial conditions. We find that when $T_{\rm ce}\lesssim 10^4$ s, late internal collisions or refreshed external collisions result in values of $T_{\rm 90}$ and $t_{\rm burst}$ larger than $T_{\rm ce}$, usually by factors of $2-3$. For $T_{\rm ce}\gtrsim 10^4$ s, the $t_{\rm burst}$ is always a good estimator for $T_{\rm ce}$, while $T_{90}$ can underpredict $T_{\rm ce}$ when the late central engine activity is moderate. We find a clear bimodal distribution for $T_{\rm ce}$, based on our simulations as well as on the observational data for $T_{90}$ and $t_{\rm burst}$. We suggest that $t_{\rm burst}$ is a reliable measure for defining "ultra-long" GRBs. Bursts with $T_{90}$ of order $10^3$ s need not belong to a special population, while bursts with $t_{\rm burst} > 10^4$ s, where the late central engine activity is more moderate and shows up in X-rays, may represent a new population. These conclusions are insensitive to the initial conditions assumed in the models.

I review gamma-ray burst models (GRBs) and observations, and discuss the possible production of ultra-high energy cosmic rays and neutrinos in both the standard internal shock models and the newer generation of photospheric and hadronic GRB models, in the light of current constraints imposed by IceCube, Auger and TA observations. I then discuss models that have been proposed to explain the recent astrophysical PeV neutrino observations, including star-forming and star-burst galaxies, hypernovae and galaxy accretion and merger shocks.

We investigate the shock acceleration of particles in massive galaxy mergers or collisions, and show that cosmic rays (CRs) can be accelerated up to the second knee energy ~0.1-1 EeV and possibly beyond, with a hard spectral index Gamma ~ 2. Such CRs lose their energy via hadronuclear interactions within a dynamical timescale of the merger shock, producing gamma rays and neutrinos as a by-product. If ~ 10 % of the shock dissipated energy goes into CR acceleration, some local merging galaxies will produce gamma-ray counterparts detectable by CTA. Also, based on the concordance cosmology, where a good fraction of the massive galaxies experience a major merger in a cosmological timescale, the neutrino counterparts can constitute ~ 20-60 % of the isotropic background detected by IceCube.

We discuss the GeV to TeV photon emission of gamma-ray bursts (GRBs) within the refreshed shock and the continuous injection scenarios, motivated by the observation of extended emission in a substantial fraction of short GRBs. In the first model we assume that the central engine emits promptly material with a range of Lorentz factors. When the fastest shell starts to decelerate, it drives a forward shock into the ambient medium and a reverse shock in the ejecta. These shocks are reenergized by the slower and later arriving material. In the second model we assume that there is a continued ejection of material over an extended time, and the continuously arriving new material keeps reenergizing the shocks formed by the preceding shells of ejecta. We calculate the synchrotron and synchrotron self-Compton radiation components for the forward and reverse shocks and find that prospective and current GeV to TeV range instruments such as CTA, HAWC, VERITAS, MAGIC and HESS have a good chance to detect afterglows of short bursts with extended emission, assuming a reasonable response time.

Papers on atmospheric and diffuse UHE neutrino searches of all flavors submitted to the 33nd International Cosmic Ray Conference (Rio de Janeiro 2013) by the IceCube Collaboration.

Recently, Thornton et al. reported the detection of four fast radio bursts (FRBs). The dispersion measures indicate that the sources of these FRBs are at cosmological distance. Given the large full sky event rate ~ 10^4 sky^-1 day^-1, the FRBs are a promising target for multi-messenger astronomy. Here we propose double degenerate, binary white-dwarf (WD) mergers as the source of FRBs, which are produced by coherent emission from the polar region of a rapidly rotating, magnetized massive WD formed after the merger. The basic characteristics of the FRBs, such as the energetics, emission duration and event rate, can be consistently explained in this scenario. As a result, we predict that some FRBs can accompany type Ia supernovae (SNe Ia) or X-ray debris disks. Simultaneous detection could test our scenario and probe the progenitors of SNe Ia, and moreover would provide a novel constraint on the cosmological parameters. We strongly encourage future SN and X-ray surveys that follow up FRBs.

We study a type of particle acceleration that operates via neutron-proton conversion in inelastic nuclear collisions. This mechanism can be expected for relativistic shocks at subphotospheres if relativistic outflows contain neutrons. Using a test-particle approximation, we numerically calculate the energy spectrum and the efficiency of accelerated particles, and show that a good energy fraction of the nucleons can be accelerated. This mechanism may especially be relevant if the shock is radiation-mediated, and it would enhance the detectability of GeV-TeV neutrinos.

Papers on point source searches submitted to the 33nd International Cosmic Ray Conference (Rio de Janeiro 2013) by the IceCube Collaboration.

Papers on neutrino oscillation and supernova searches submitted to the 33nd International Cosmic Ray Conference (Rio de Janeiro 2013) by the IceCube Collaboration.

Papers on searches for Dark Matter and exotic particles submitted to the 33nd International Cosmic Ray Conference (Rio de Janeiro 2013) by the IceCube Collaboration.

Papers on ice properties, reconstruction and future developments submitted to the 33nd International Cosmic Ray Conference (Rio de Janeiro 2013) by the IceCube Collaboration.

Papers on cosmic rays submitted to the 33nd International Cosmic Ray Conference (Rio de Janeiro 2013) by the IceCube Collaboration.

We present a dissipative photospheric model for gamma-ray bursts where the usual Band peak around MeV energies arises as synchrotron emission from around the photosphere. We treat the initial acceleration in a general way and the GeV emission arises as the interaction of photospheric radiation and the shocked electrons at the deceleration radius. We show some applications of this model.

Thermal photons from the photosphere may be the primary source of the observed prompt emission of gamma-ray bursts (GRBs). In order to produce the observed non-thermal spectra, some kind of dissipation mechanism near the photosphere is required. In this paper we numerically simulate the evolution of the photon spectrum in a relativistically expanding shell with a time-dependent numerical code. We consider two basic models. One is a leptonic model, where a dissipation mechanism heats the thermal electrons maintaining their high temperature. The other model involves a cascade process induced by $pp$($pn$)-collisions which produce high-energy electrons, modify the thermal spectrum, and emit neutrinos. The qualitative properties of the photon spectra are mainly determined by the optical depth at which the dissipation mechanism sets in. Too large optical depths lead to a broad and curved spectrum contradicting the observations, while for optical depths smaller than unity the spectral hardness becomes softer than observed. A significant shift of the spectral peak energy to higher energies due to a large energy injection can lead to an overly broad spectral shape. We show ideal parameter ranges for which these models are able to reproduce the observed spectra. For the $pn$-collision model, the neutrino fluence in the 10-100 GeV range is well above the atmospheric neutrino fluence, but its detection is challenging for presently available detectors.

We report on the observation of two neutrino-induced events which have an estimated deposited energy in the IceCube detector of 1.04 $\pm$ 0.16 and 1.14 $\pm$ 0.17 PeV, respectively, the highest neutrino energies observed so far. These events are consistent with fully contained particle showers induced by neutral-current $\nu_{e,\mu,\tau}$ ($\bar\nu_{e,\mu,\tau}$) or charged-current $\nu_{e}$ ($\bar\nu_{e}$) interactions within the IceCube detector. The events were discovered in a search for ultra-high energy neutrinos using data corresponding to 615.9 days effective livetime. The expected number of atmospheric background is $0.082 \pm 0.004 \text{(stat)}^{+0.041}_{-0.057} \text{(syst)}$. The probability to observe two or more candidate events under the atmospheric background-only hypothesis is $2.9\times10^{-3}$ ($2.8\sigma$) taking into account the uncertainty on the expected number of background events. These two events could be a first indication of an astrophysical neutrino flux, the moderate significance, however, does not permit a definitive conclusion at this time.

We report on the observation of a significant deficit of cosmic rays from the direction of the Moon with the IceCube detector. The study of this "Moon shadow" is used to characterize the angular resolution and absolute pointing capabilities of the detector. The detection is based on data taken in two periods before the completion of the detector: between April 2008 and May 2009, when IceCube operated in a partial configuration with 40 detector strings deployed in the South Pole ice, and between May 2009 and May 2010 when the detector operated with 59 strings. Using two independent analysis methods, the Moon shadow has been observed to high significance (> 6 sigma) in both detector configurations. The observed location of the shadow center is within 0.2 degrees of its expected position when geomagnetic deflection effects are taken into account. This measurement validates the directional reconstruction capabilities of IceCube.

Gamma-Ray Bursts (GRBs) have been proposed as a leading candidate for acceleration of ultra high-energy cosmic rays, which would be accompanied by emission of TeV neutrinos produced in proton-photon interactions during acceleration in the GRB fireball. Two analyses using data from two years of the IceCube detector produced no evidence for this neutrino emission, placing strong constraints on models of neutrino and cosmic-ray production in these sources.

Relativistic shocks that accompany supernovae (SNe) produce X-ray burst emissions as they break out in the dense circumstellar medium around the progenitors. This phenomenon is sometimes associated with peculiar low-luminosity gamma-ray bursts (LL GRBs). Here, we investigate the high energy neutrino and gamma-ray counterparts of such a class of SNe. Just beyond the shock breakout radius, particle acceleration in the collisionless shock starts to operate in the presence of breakout photons. We show that protons may be accelerated to sufficiently high energies and produce high energy neutrinos and gamma rays via the photomeson interaction. These neutrinos and gamma rays may be detectable from 10 Mpc away by IceCube/KM3Net as multi-TeV transients almost simultaneously with the X-ray burst emission, and even from 100 Mpc away with follow-up observations by CTA using a wide-field sky monitor like Swift as a trigger. A statistical technique using a stacking approach could also be possible for the detection, with the aid of the SN optical/infrared counterparts. Such multi-messenger observations offer the possibility to probe the transition of relativistic shocks from radiation-mediated to collisionless ones, and would also constrain the mechanisms of particle acceleration and emission in LL GRBs.

A search for muon neutrinos from dark matter annihilations in the Galactic Center region has been performed with the 40-string configuration of the IceCube Neutrino Observatory using data collected in 367 days of live-time starting in April 2008. The observed fluxes were consistent with the atmospheric background expectations. Upper limits on the self-annihilation cross-section are obtained for dark matter particle masses ranging from 100 GeV to 10 TeV. In the case of decaying dark matter, lower limits on the lifetime have been determined for masses between 200 GeV and 20 TeV.

The IceCube Neutrino Observatory, approximately 1 km^3 in size, is now complete with 86 strings deployed in the Antarctic ice. IceCube detects the Cherenkov radiation emitted by charged particles passing through or created in the ice. To realize the full potential of the detector, the properties of light propagation in the ice in and around the detector must be well understood. This report presents a new method of fitting the model of light propagation in the ice to a data set of in-situ light source events collected with IceCube. The resulting set of derived parameters, namely the measured values of scattering and absorption coefficients vs. depth, is presented and a comparison of IceCube data with simulations based on the new model is shown.

We report on searches for neutrino sources at energies above 200 GeV in the Northern sky of the galactic plane, using the data collected by the South Pole neutrino telescopes IceCube and AMANDA. The galactic region considered here includes the Local Arm towards the Cygnus region and our closest approach to the Perseus Arm. The data have been collected between 2007 and 2009 when AMANDA was an integrated part of IceCube, which was still under construction and operated with 22-strings (2007-8) and 40-strings (2008-9) of optical modules deployed in the ice. By combining the larger IceCube detector with the lower energy threshold of the more compact AMANDA detector, we obtain an improved sensitivity at energies below $\sim$10 TeV with respect to previous searches. The analyses presented here are: a scan for point sources within the galactic plane; a search optimized for multiple and extended sources in the Cygnus region, which might be below the sensitivity of the point source scan; and studies of seven pre-selected neutrino source candidates. For one of them, Cygnus X-3, a time-dependent search for neutrinos in coincidence with observed radio and X-ray flares has been performed. No evidence of a signal is found, and upper limits are reported for each of the searches. We investigate neutrino spectra proportional to E$^{-2}$ and E$^{-3}$ to cover the entire range of possible spectra. The soft E$^{-3}$ spectrum results in an energy distribution similar to a source with cut-off below $\sim$50 TeV. For the considered region of the galactic plane, the 90% confidence level muon neutrino flux upper limits are in the range E$^3$dN/dE$\sim 5.4 - 19.5 \times 10^{-11} \rm{TeV^{2} cm^{-2} s^{-1}}$ for point-like neutrino sources in the energy region [180.0 GeV - 20.5 TeV]. These represent the most stringent upper limits for soft-spectra neutrino sources within the Galaxy reported to date.

We consider gamma-ray burst models where the radiation is dominated by a photospheric region providing the MeV Band spectrum, and an external shock region responsible for the GeV radiation via inverse Compton scattering. We parametrize the initial dynamics through an acceleration law Gamma \propto r^\mu, with \mu between 1/3 and 1 to represent the range between an extreme magnetically dominated and a baryonically dominated regime. We compare these models to several bright Fermi LAT bursts, and show that both the time integrated and the time resolved spectra, where available, can be well described by these models. We discuss the parameters which result from these fits, and discuss the relative merits and shortcomings of the two models.

We report on the observation of anisotropy in the arrival direction distribution of cosmic rays at PeV energies. The analysis is based on data taken between 2009 and 2012 with the IceTop air shower array at the South Pole. IceTop, an integral part of the IceCube detector, is sensitive to cosmic rays between 100 TeV and 1 EeV. With the current size of the IceTop data set, searches for anisotropy at the 10^-3 level can, for the first time, be extended to PeV energies. We divide the data set into two parts with median energies of 400 TeV and 2 PeV, respectively. In the low energy band, we observe a strong deficit with an angular size of about 30 degrees and an amplitude of (-1.58 +/- 0.46 (stat) +/- 0.52 (sys)) x 10^(-3) at a location consistent with previous observations of cosmic rays with the IceCube neutrino detector. The study of the high energy band shows that the anisotropy persists to PeV energies and increases in amplitude to (-3.11 +/- 0.38 (stat) +/- 0.96 (sys)) x 10^(-3).

In cosmic ray air showers, the muon lateral separation from the center of the shower is a measure of the transverse momentum that the muon parent acquired in the cosmic ray interaction. IceCube has observed cosmic ray interactions that produce muons laterally separated by up to 400 m from the shower core, a factor of 6 larger distance than previous measurements. These muons originate in high pT (> 2 GeV/c) interactions from the incident cosmic ray, or high-energy secondary interactions. The separation distribution shows a transition to a power law at large values, indicating the presence of a hard pT component that can be described by perturbative quantum chromodynamics. However, the rates and the zenith angle distributions of these events are not well reproduced with the cosmic ray models tested here, even those that include charm interactions. This discrepancy may be explained by a larger fraction of kaons and charmed particles than is currently incorporated in the simulations.

We summarize the science opportunity, design elements, current and projected partner observatories, and anticipated science returns of the Astrophysical Multimessenger Observatory Network (AMON). AMON will link multiple current and future high-energy, multimessenger, and follow-up observatories together into a single network, enabling near real-time coincidence searches for multimessenger astrophysical transients and their electromagnetic counterparts. Candidate and high-confidence multimessenger transient events will be identified, characterized, and distributed as AMON alerts within the network and to interested external observers, leading to follow-up observations across the electromagnetic spectrum. In this way, AMON aims to evoke the discovery of multimessenger transients from within observatory subthreshold data streams and facilitate the exploitation of these transients for purposes of astronomy and fundamental physics. As a central hub of global multimessenger science, AMON will also enable cross-collaboration analyses of archival datasets in search of rare or exotic astrophysical phenomena.

We present the first results in the search for relativistic magnetic monopoles with the IceCube detector, a subsurface neutrino telescope located in the South Polar ice cap containing a volume of 1 km$^{3}$. This analysis searches data taken on the partially completed detector during 2007 when roughly 0.2 km$^{3}$ of ice was instrumented. The lack of candidate events leads to an upper limit on the flux of relativistic magnetic monopoles of $\Phi_{\mathrm{90%C.L.}}\sim 3\e{-18}\fluxunits$ for $\beta\geq0.8$. This is a factor of 4 improvement over the previous best experimental flux limits up to a Lorentz boost $\gamma$ below $10^{7}$. This result is then interpreted for a wide range of mass and kinetic energy values.

The Fermi observations of GRB 110721A (Axelsson et al. 2012) have revealed an unusually high peak energy ~ 15 MeV in the first time bin of the prompt emission. We find that an interpretation is unlikely in terms of internal shock models, and confirm that a standard black-body photospheric model also falls short. We show that dissipative photospheric synchrotron models ranging from extreme magnetically dominated to baryon dominated dynamics, on the other hand, are able to accommodate such high peak values.

The measurement of muon energy is critical for many analyses in large Cherenkov detectors, particularly those that involve separating extraterrestrial neutrinos from the atmospheric neutrino background. Muon energy has traditionally been determined by measuring the specific energy loss (dE/dx) along the muon's path and relating the dE/dx to the muon energy. Because high-energy muons (E_mu > 1 TeV) lose energy randomly, the spread in dE/dx values is quite large, leading to a typical energy resolution of 0.29 in log10(E_mu) for a muon observed over a 1 km path length in the IceCube detector. In this paper, we present an improved method that uses a truncated mean and other techniques to determine the muon energy. The muon track is divided into separate segments with individual dE/dx values. The elimination of segments with the highest dE/dx results in an overall dE/dx that is more closely correlated to the muon energy. This method results in an energy resolution of 0.22 in log10(E_mu), which gives a 26% improvement. This technique is applicable to any large water or ice detector and potentially to large scintillator or liquid argon detectors.

The mass composition of high energy cosmic rays depends on their production, acceleration, and propagation. The study of cosmic ray composition can therefore reveal hints of the origin of these particles. At the South Pole, the IceCube Neutrino Observatory is capable of measuring two components of cosmic ray air showers in coincidence: the electromagnetic component at high altitude (2835 m) using the IceTop surface array, and the muonic component above ~1 TeV using the IceCube array. This unique detector arrangement provides an opportunity for precision measurements of the cosmic ray energy spectrum and composition in the region of the knee and beyond. We present the results of a neural network analysis technique to study the cosmic ray composition and the energy spectrum from 1 PeV to 30 PeV using data recorded using the 40-string/40-station configuration of the IceCube Neutrino Observatory.

Gamma-ray induced air showers are notable for their lack of muons, compared to hadronic showers. Hence, air shower arrays with large underground muon detectors can select a sample greatly enriched in photon showers by rejecting showers containing muons. IceCube is sensitive to muons with energies above ~500 GeV at the surface, which provides an efficient veto system for hadronic air showers with energies above 1 PeV. One year of data from the 40-string IceCube configuration was used to perform a search for point sources and a Galactic diffuse signal. No sources were found, resulting in a 90% C.L. upper limit on the ratio of gamma rays to cosmic rays of 1.2 x 10^(-3)for the flux coming from the Galactic Plane region (-80 deg < l < -30 deg; -10 deg < b < 5 deg) in the energy range 1.2 - 6.0 PeV. In the same energy range, point source fluxes with E^(-2) spectra have been excluded at a level of (E/TeV)^2 d\Phi/dE ~ 10^(-12)-10^(-11) cm^2/s/TeV depending on source declination. The complete IceCube detector will have a better sensitivity, due to the larger detector size, improved reconstruction and vetoing techniques. Preliminary data from the nearly-final IceCube detector configuration has been used to estimate the 5 year sensitivity of the full detector. It is found to be more than an order of magnitude better, allowing the search for PeV extensions of known TeV gamma-ray emitters.

We present a fast likelihood method for including event-level neutrino telescope data in parameter explorations of theories for new physics, and announce its public release as part of DarkSUSY 5.0.6. Our construction includes both angular and spectral information about neutrino events, as well as their total number. We also present a corresponding measure for simple model exclusion, which can be used for single models without reference to the rest of a parameter space. We perform a number of supersymmetric parameter scans with IceCube data to illustrate the utility of the method: example global fits and a signal recovery in the constrained minimal supersymmetric standard model (CMSSM), and a model exclusion exercise in a 7-parameter phenomenological version of the MSSM. The final IceCube detector configuration will probe almost the entire focus-point region of the CMSSM, as well as a number of MSSM-7 models that will not otherwise be accessible to e.g. direct detection. Our method accurately recovers the mock signal, and provides tight constraints on model parameters and derived quantities. We show that the inclusion of spectral information significantly improves the accuracy of the recovery, providing motivation for its use in future IceCube analyses.

Gamma-ray bursts are the most concentrated explosions in the Universe. They have been detected electromagnetically at energies up to tens of GeV, and it is suspected that they could be active at least up to TeV energies. It is also speculated that they could emit cosmic rays and neutrinos at energies reaching up to the $10^{18}-10^{20}$ eV range. Here we review the recent developments in the photon phenomenology in the light of \swift and \fermi satellite observations, as well as recent IceCube upper limits on their neutrino luminosity. We discuss some of the theoretical models developed to explain these observations and their possible contribution to a very high energy cosmic ray and neutrino background.

Gamma-ray bursts are the most luminous explosions in the Universe, whose origin and mechanism is the focus of intense interest. They appear connected to supernova remnants from massive stars or the merger of their remnants, and their brightness makes them temporarily detectable out to the larges distances yet explored in the Universe. After pioneering breakthroughs from space and ground experiments, their study is entering a new phase with observations from the recently launched \fermi satellite, as well as the prospect of detections or limits from large neutrino and gravitational wave detectors. The interplay between such observations and theoretical models of gamma-ray bursts is reviewed, as well as their connections to supernovae and cosmology.

The temporal--spectral evolution of the prompt emission of gamma-ray bursts (GRBs) is simulated numerically for both leptonic and hadronic models. For weak enough magnetic fields, leptonic models can reproduce the few seconds delay of the onset of GeV photon emission observed by Fermi-LAT, due to the slow growth of the target photon field for inverse Compton scattering. However, even for stronger magnetic fields, the GeV delay can be explained with hadronic models, due to the long acceleration timescale of protons and the continuous photopion production after the end of the particle injection. While the FWHMs of the MeV and GeV lightcurves are almost the same in one-zone leptonic models, the FWHM of the 1--30 GeV lightcurves in hadronic models are significantly wider than those of the 0.1--1 MeV lightcurves. The amount of the GeV delay depends on the importance of the Klein--Nishina effect in both the leptonic and hadronic models. In our examples of hadronic models the energies of the escaped neutrons are comparable to the gamma-ray energy, although their contribution to the ultra high-energy cosmic rays is still subdominant. The resulting neutrino spectra are hard enough to avoid the flux limit constraint from IceCube. The delay of the neutrino emission onset is up to several times longer than the corresponding delay of the GeV photon emission onset. The quantitative differences in the lightcurves for various models may be further tested with future atmospheric Cherenkov telescopes whose effective area is larger than that of Fermi-LAT, such as CTA.

IceTop, the surface component of the IceCube Neutrino Observatory at the South Pole, is an air shower array with an area of 1 km2. The detector allows a detailed exploration of the mass composition of primary cosmic rays in the energy range from about 100 TeV to 1 EeV by exploiting the correlation between the shower energy measured in IceTop and the energy deposited by muons in the deep ice. In this paper we report on the technical design, construction and installation, the trigger and data acquisition systems as well as the software framework for calibration, reconstruction and simulation. Finally the first experience from commissioning and operating the detector and the performance as an air shower detector will be discussed.

The first dedicated search for ultra-high energy (UHE) tau neutrinos of astrophysical origin was performed using the IceCube detector in its 22-string configuration with an instrumented volume of roughly 0.25 km^3. The search also had sensitivity to UHE electron and muon neutrinos. After application of all selection criteria to approximately 200 live-days of data, we expect a background of 0.60 +/- 0.19 (stat.) $^{+0.56}_{-0.58}$ (syst.) events and observe three events, which after inspection emerge as being compatible with background but are kept in the final sample. Therefore, we set an upper limit on neutrinos of all-flavors from UHE astrophysical sources at 90% CL of $E^{2} \Phi(\nu_{x}) < 16.3 * 10^-8 GeV cm^-2 sr^-1 s^-1 over an estimated primary neutrino energy range of 340 TeV to 200 PeV.

Most Fermi GRB spectra appear as either a broken power law extending to GeV energies or as a broken power with a separate GeV power law component. Here we show that such spectra can be understood in terms of magnetically dominated relativistic jets where a dissipative photosphere produces the prompt MeV emission, which is extended into the GeV range by inverse Compton scattering in the external shock, with possible contributions from a reverse shock as well. The bulk Lorentz factors required in these models are in the range of 300-600, and the MeV-GeV time delays arise naturally. In some cases an optical flash and a sub-dominant thermal component are also present.

Recently, Fermi-LAT detected GeV emission during the X-ray flaring activity in GRB 100728A. We study various scenarios for its origin. The hard spectrum of the GeV emission favors the external inverse-Compton origin in which X-ray flare photons are up-scattered by relativistic electrons in the external forward shock. This external IC scenario, with anisotropic scattering effect taken into account, can reproduce the temporal and spectral properties of the GeV emission in GRB 100728A.

Gamma-ray bursts have been detected at photon energies up to tens of GeV. We review some recent developments in the X-ray to GeV photon phenomenology in the light of Swift and Fermi observations, and some of the theoretical models developed to explain them, with a view towards implications for C.T.A.