The discovery of high-energy astrophysical neutrinos by IceCube has opened a new window to the Universe. However, the origin of these neutrinos is still a mystery, and some of them could be a result of dark matter interactions such as decay. Next generation gigaton water-Cherenkov neutrino telescope, KM3NeT, is expected to offer significantly improved energy resolution in the cascade channel, and advantageous viewing condition to the Galactic Center; both important for searches of dark matter decay signals. We study the sensitivity of KM3NeT on dark matter decays by performing a mock likelihood analysis for both cascade and track type events, taking into account both angular and energy information. We find that, combining both channels, KM3NeT is expected to produce world leading limits on dark matter decay lifetime in the PeV mass range, and could test some of the dark matter hints in the current IceCube data.

Cosmic rays in galaxy clusters are unique probes of energetic processes operating with large-scale structures in the Universe. Precise measurements of cosmic rays in galaxy clusters are essential for improving our understanding of non-thermal components in the intracluster-medium (ICM) as well as the accuracy of cluster mass estimates in cosmological analyses. In this paper, we perform a cross-correlation analysis with the extragalactic gamma-ray background and the thermal Sunyaev-Zel'dovich (tSZ) effect in the cosmic microwave background. The expected cross-correlation signal would contain rich information about the cosmic-ray-induced gamma-ray emission in the most massive galaxy clusters at $z\sim0.1-0.2$. We analyze the gamma-ray background map with 8 years of data taken by the Large Area Telescope onboard Fermi satellite and the publicly available tSZ map by Planck. We confirm that the measured cross-correlation is consistent with a null detection, and thus it enables us to put the tightest constraint on the acceleration efficiency of cosmic ray protons at shocks in and around galaxy clusters. We find the acceleration efficiency must be below 5\% with a $2\sigma$ confidence level when the hydrostatic mass bias of clusters is assumed to be 30\%, and our result is not significantly affected by the assumed value of the hydrostatic mass bias. Our constraint implies that the non-thermal cosmic-ray pressure in the ICM can introduce only a $\le 3\%$ level of the hydrostatic mass bias, highlighting that cosmic rays alone do not account for the mass bias inferred by the Planck analyses. Finally, we discuss future detectability prospects of cosmic-ray-induced gamma rays from the Perseus cluster for the Cherenkov Telescope Array.

The hypothesis of two different components in the high-energy neutrino flux observed with IceCube has been proposed to solve the tension among different data-sets and to account for an excess of neutrino events at 100 TeV. In addition to a standard astrophysical power-law component, the second component might be explained by a different class of astrophysical sources, or more intriguingly, might originate from decaying or annihilating dark matter. These two scenarios can be distinguished thanks to the different expected angular distributions of neutrino events. Neutrino signals from dark matter are indeed expected to have some correlation with the extended galactic dark matter halo. In this paper, we perform angular power spectrum analyses of simulated neutrino sky maps to investigate the two-component hypothesis with a contribution from dark matter. We provide current constraints and expected sensitivity to dark matter parameters for future neutrino telescopes such as IceCube-Gen2 and KM3NeT. The latter is found to be more sensitive than IceCube-Gen2 to look for a dark matter signal at low energies towards the galactic center. Finally, we show that after 10 years of data-taking, they will firmly probe the current best-fit scenario for decaying dark matter by exploiting the angular information only.

Dwarf spheroidal galaxies that form in halo substructures provide stringent constraints on dark matter annihilation. Many ultrafaint dwarfs discovered with modern surveys contribute significantly to these constraints. At present, because of the lack of abundant stellar kinematic data for the ultrafaints, non-informative prior assumptions are usually made for the parameters of the density profiles. Based on semi-analytic models of dark matter subhalos and their connection to satellite galaxies, we present more informative and realistic satellite priors. We show that our satellite priors lead to constraints on the annihilation rate that are between a factor of 2 and a factor of 7 weaker than under non-informative priors. As a result, the thermal relic cross section can at best only be excluded (with 95% probability) for dark matter masses of $\lesssim 40$ GeV from dwarf spheroidal data, assuming annihilation into $b\bar{b}$.

In recent years, many gamma-ray sources have been identified, yet the unresolved component hosts valuable information on the faintest emission. In order to extract it, a cross-correlation with gravitational tracers of matter in the Universe has been shown to be a promising tool. We report here the first identification of a cross-correlation signal between gamma rays and the distribution of mass in the Universe probed by weak gravitational lensing. We use the Dark Energy Survey Y1 weak lensing catalogue and the Fermi Large Area Telescope 9-year gamma-ray data, obtaining a signal-to-noise ratio of 5.3. The signal is mostly localised at small angular scales and high gamma-ray energies, with a hint of correlation at extended separation. Blazar emission is likely the origin of the small-scale effect. We investigate implications of the large-scale component in terms of astrophysical sources and particle dark matter emission.

The free streaming motion of dark matter particles imprints a cutoff in the matter power spectrum and set the scale of the smallest dark matter halo. Recent cosmological $N$-body simulations have shown that the central density cusp is much steeper in haloes near the free streaming scale than in more massive haloes. Here, we study the abundance and structure of subhaloes near the free streaming scale at very high redshift using a suite of unprecedentedly large cosmological $N$-body simulations, over a wide range of the host halo mass. The subhalo abundance is suppressed strongly below the free streaming scale, but the ratio between the subhalo mass function in the cutoff and no cutoff simulations is well fitted by a single correction function regardless of the host halo mass and the redshift. In subhaloes, the central slopes are considerably shallower than in field haloes, however, are still steeper than that of the NFW profile. Contrary, the concentrations are significantly larger in subhaloes than haloes and depend on the subhalo mass. We compare two methods to extrapolate the mass-concentration relation of haloes and subhaloes to z=0 and provide a new simple fitting function for subhaloes, based on a suite of large cosmological $N$-body simulations. Finally, we estimate the annihilation boost factor of a Milky-Way sized halo to be between 1.8 and 6.2.

Searches for dark matter annihilation signals have been carried out in a number of target regions such as the Galactic Center and Milky Way dwarf spheroidal galaxies (dSphs), among a few others. Here we propose low surface brightness galaxies (LSBGs) asnovel targets for the indirect detection of dark matter emission. In particular, LSBGs are known to have very large dark matter contents and be less contaminated by extragalactic gamma-ray sources (e.g., blazars) compared to star forming galaxies. We report on an analysis that uses eight LSBGs (detected by Subaru Hyper Suprime-Cam survey data) with known redshifts to conduct a search for gamma-ray emission at the positions of these new objects in Fermi Large Area Telescope data. We found no excesses of gamma-ray emission and set constraints on the dark matter annihilation cross-section. We exclude (at the 95% C.L.) dark matter scenarios predicting a cross-section higher than 10^-23[cm^3/s] for dark matter particles of mass 10 GeV self-annihilating in the b_b channel. Although this constraint is weaker than the ones reported in recent studies using other targets, we note that in the near future, the number of detections of new LSBGs will increase by a few orders of magnitude. We forecast that with the use of the full catalog of soon-to-be-detected LSBGs the constraint will reach cross-section sensitivities of ~ 3*10^-25 [cm^3/s] for dark matter particles with masses less than 10 GeV.

We set conservative constraints on decaying dark matter particles with masses spanning a very wide range ($10^4-10^{16}$ GeV). For this we use multimessenger observations of cosmic-ray (CR) protons/antiprotons, electrons/positrons, neutrinos/antineutrinos and gamma rays. Focusing on decays into the $\bar{b}b$ channel, we simulate the spectra of dark matter yields by using the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi equations and the $\texttt{Pythia}$ package. We then propagate the CRs of dark matter origin till Earth by using the state-of-the-art numerical frameworks $\texttt{CRPropa}$, $\texttt{GALPROP}$ and $\texttt{HelMod}$ for the solution of the CR transport equation in the extragalactic, Galactic region and the heliosphere, respectively. Conservative limits are obtained by requiring that the predicted dark matter spectra at Earth be less than the observed CR spectra. Overall, we exclude dark matter lifetimes of $10^{28}$ s or shorter for all the masses investigated in this work. The most stringent constraints reach $10^{30}$ s for very heavy dark matter particles with masses in the range $10^{11}-10^{14}$ GeV.

(Mini) split supersymmetry explains the observed Higgs mass and evades stringent constraints, while keeping good features of TeV-scale supersymmetry other than the little hierarchy problem. Such scenarios naturally predict thermal wino dark matter whose mass is around $3 \, {\rm TeV}$. Its non-perturbatively enhanced annihilation is a promising target of indirect detection experiments. It is known that identifying the smallest halos is essential for reducing an uncertainty in interpreting indirect detection experiments. Despite its importance, the smallest halos of thermal wino dark matter have not been well understood and thus are investigated in this work. In particular, we remark on two aspects: 1) the neutral wino is in kinetic equilibrium with primordial plasma predominantly through inelastic processes involving the slightly heavier charged wino; and 2) the resultant density contrast shows larger powers at dark acoustic oscillation peaks than in cold dark matter, which is known as an overshooting phenomenon. By taking them into account, we provide a rigorous estimate of the boost factor. Our result facilitates accurately pinning down thermal wino dark matter through vigorous efforts in indirect detection experiments.

We study the prospects for the Large Synoptic Survey Telescope (LSST) to find new dwarf spheroidal galaxies in the Milky Way. Adopting models of Milky-Way halo substructure and phenomenological prescriptions connecting subhalos and satellite galaxies, we obtain surface brightness distributions of $V$-band magnitude that lead us to predict that LSST will discover tens to hundreds of dwarf spheroidal galaxies above its sensitivity. The soon-to-be-discovered dwarfs will be interesting targets for indirect searches of dark matter annihilation yields. We forecast the distribution function of gamma-ray emission from dark matter annihilation in these objects, and discuss the detectability of these signals at both Fermi Large Area Telescope (LAT) and Cherenkov Telescope Array (CTA). By combining information from the predicted dwarf galaxies, we obtain an expected sensitivity to the annihilation cross section $\langle \sigma v \rangle$ of $10^{-26}$ cm$^3$ s$^{-1}$ (for dark matter particles of mass 10 GeV with Fermi-LAT) and $5\times 10^{-24}$ cm$^3$ s$^{-1}$ (for dark matter particles of mass 500 GeV with CTA). We find that the current uncertainties in the mass measurement of the Milky-Way halo are relatively minor compared with the Poisson errors associated to drawing the most promising dwarfs from the underlying flux distribution.

The presence of dark matter substructure will boost the signatures of dark matter annihilation. We review recent progress on estimates of this subhalo boost factor---a ratio of the luminosity from annihilation in the subhalos to that originating the smooth component---based on both numerical $N$-body simulations and semi-analytic modelings. Since subhalos of all the scales, ranging from the Earth mass (as expected, e.g., the supersymmetric neutralino, a prime candidate for cold dark matter) to galaxies or larger, give substantial contribution to the annihilation rate, it is essential to understand subhalo properties over a large dynamic range of more than twenty orders of magnitude in masses. Even though numerical simulations give the most accurate assessment in resolved regimes, extrapolating the subhalo properties down in sub-grid scales comes with great uncertainties---a straightforward extrapolation yields a very large amount of the subhalo boost factor of $\gtrsim$100 for galaxy-size halos. Physically motivated theoretical models based on analytic prescriptions such as the extended Press-Schechter formalism and tidal stripping modeling, which are well tested against the simulation results, predict a more modest boost of order unity for the galaxy-size halos. Giving an accurate assessment of the boost factor is essential for indirect dark matter searches and thus, having models calibrated at large ranges of host masses and redshifts, is strongly urged upon.

Astrophysical neutrino events have been measured in the last couple of years, which show an isotropic distribution, and the current discussion is their astrophysical origin. We use both isotropic and anisotropic components of the diffuse neutrino data to constrain the contribution of a broad number of extra-galactic source populations to the observed neutrino sky. We simulate up-going muon neutrino events by applying statistical distributions for the flux of extragalactic sources, and by Monte Carlo method we exploit the simulation for current and future IceCube, IceCube-Gen2 and KM3NeT exposures. We aim at constraining source populations by studying their angular patterns, for which we assess the angular power spectrum. We leave the characteristic number of sources ($N_{\star}$) as a free parameter, which is roughly the number of neutrino sources over which the measured intensity is divided. With existing two-year IceCube data, we can already constrain very rare, bright sources with $N_{\star}\lesssim$100. This can be improved to $N_{\star}\lesssim 10^4$-$10^5$ with IceCube-Gen2 and KM3NeT with ten-year exposure, constraining the contribution of BL Lacs ($N_{\star}=6\times10^{2}$). On the other hand, we can constrain weak sources with large number densities, like starburst galaxies ($N_{\star} = 10^{7}$), if we measure an anisotropic neutrino sky with future observations.

We present a first joint analysis of the power spectra of the thermal Sunyaev-Zeldovich (tSZ) effect measured by the Planck and the number density fluctuations of galaxies in the 2MASS redshift survey (2MRS) catalog, including their cross-correlation. Combining these measurements with the cosmic microwave background (CMB) data and CMB lensing of Planck assuming a flat $\Lambda$CDM model, we constrain the mass bias parameter as $B = 1.54 \pm 0.098\;(1\sigma)$ [$(1-b) = 0.649 \pm 0.041$, where $(1-b) \equiv B^{-1}$], i.e., the Planck cluster mass should be $35\%$ lower than the true mass. The mass bias determined by the 2MRS-tSZ cross-power spectrum alone is consistent with that determined by the tSZ auto-power spectrum alone, suggesting that this large mass bias is not due to obvious systematics in the tSZ data. We find that the 2MRS-tSZ cross-power spectrum is more sensitive to less massive halos than the tSZ auto-power spectrum and it significantly improves a constraint on the mass dependence of the mass bias. The redshift dependence is not strongly constrained since the multipole range in which high redshift clusters mainly contribute to the tSZ auto is dominated by the contaminating sources. We conclude that no strong mass or redshift evolution of the mass bias is needed to explain the data.

The recent results of the EDGES collaboration indicate that during the era of reionization, the primordial gas was much colder than expected. The cooling of the gas could be explained by interactions between dark matter (DM) and particles in the primordial gas. Constraints from cosmology and particle experiments indicate that this DM should be light ($\sim$10-80 MeV), carry a small charge ($\epsilon\sim 10^{-6}$-$10^{-4}$), and only make up a small fraction of the total amount of DM. Several constraints on the DM parameter space have already been made. We explore the yet unconstrained region in the case that the milli-charged DM makes up for $\sim$2\% of the total dark matter, through the scenario in which this DM annihilates only into mu and tau neutrinos. We set upper limits on the annihilation cross section using the Super-Kamiokande data, and predict the limits that could be obtained through Hyper-Kamiokande, JUNO and DUNE. We find that data from Super-Kamiokande is not yet able to constrain this model, but future experiments might be. We furthermore explore DM annihilation into solely neutrinos in general, giving an update of the current limits, and predict the limits that could be placed with future experiments.

We make a first attempt to find dwarf galaxies in eight \Fermi-LAT extended, unassociated, source fields using \Gaia\ DR2. We probe previously unexplored heliocentric distances of $d<20$~kpc with an extreme-deconvolution (XD) technique. We find no signature of a dwarf galaxy in any of these fields despite \Gaia's excellent astrometric accuracy. We estimate our detection limits by applying the XD method to mock data, obtaining a conservative limit on the stellar mass of $M_* < 10^4$~M$_{\sun}$ for $d < 20$\, kpc. Such a low stellar mass implies either a low-mass subhalo, or a massive stripped-down subhalo. We use an analytic model for stripped subhalos to argue that, given the sizes and fluxes of the \Fermi-LAT sources, we can reject the hypothesis that they owe to dark matter annihilation.

We investigate the radiative decay of the cosmic neutrino background, and its impact on the spectrum of the cosmic microwave background (CMB) that is known to be a nearly perfect black body. We derive exact formulae for the decay of a heavier neutrino into a lighter neutrino and a photon, $\nu_j \to \nu_i + \gamma$, and of absorption as its inverse, $\nu_i + \gamma \to \nu_j$, by accounting for the precise form of the neutrino momentum distribution. Our calculations show that if the neutrinos are heavier than $\mathcal O(0.1)$ eV, the exact formulae give results that differ by $\sim$50%, compared with approximate ones where neutrinos are assumed to be at rest. We also find that spectral distortion due to absorption is more important for heavy neutrino masses (by a factor of $\sim$10 going from a neutrino mass of 0.01 eV to 0.1 eV). By analyzing the CMB spectral data measured with COBE-FIRAS, we obtain lower limits on the neutrino lifetime of $\tau_{12} \gtrsim 4 \times 10^{21}$ s (95% C.L.) for the smaller mass splitting and $\tau_{13} \sim \tau_{23} \gtrsim 10^{19}$ s for the larger mass splitting. These represent up to one order of magnitude improvement over previous CMB constraints. With future CMB experiments such as PIXIE, these limits will improve by roughly 4 orders of magnitude. This translates to a projected upper limit on the neutrino magnetic moment (for certain neutrino masses and decay modes) of $\mu_\nu < 3 \times 10^{-11}\, \mu_B$, where $\mu_B$ is the Bohr magneton. Such constraints would make future precision CMB measurements competitive with lab-based constraints on neutrino magnetic moments.

The origins of the extraterrestrial neutrinos observed in IceCube have yet to be determined. In this study we perform a one-point fluctuation analysis of the six-year high-energy starting event (HESE) shower data, with fixed non-Poissonian contributions from atmospheric, Galactic and some extragalactic components, as well as an isotropic (and weakly non-Poissonian) template. In addition to the star-forming galaxies and blazars, our analysis suggests the presence of an additional isotropic component, not associated with any known class of sources, with best-fit intensity of $(2.8\pm0.2)\times 10^{-18}\,(E/100~{\rm TeV})^{-2.7\pm 0.5}$ ${\mathrm{cm}^{-2}\, \mathrm{s}^{-1}\,\mathrm{sr}^{-1}\,\mathrm{GeV}^{-1}}$. For the first time, we also consider high-energy extrapolations of several phenomenological models for the diffuse Galactic emission (tuned to both local cosmic-ray data and diffuse gamma-ray emission in the GeV-TeV domain). We demonstrate the potential of our framework in discriminating between different scenarios, with possible implications on the physics of cosmic ray transport in the TeV-PeV range.

We study evolution of dark matter substructures, especially how they lose the mass and change density profile after they fall in gravitational potential of larger host halos. We develop an analytical prescription that models the subhalo mass evolution and calibrate it to results of N-body numerical simulations of various scales from very small (Earth size) to large (galaxies to clusters) halos. We then combine the results with halo accretion histories, and calculate the subhalo mass function that is physically motivated down to Earth-mass scales. Our results --- valid for arbitrary host masses and redshifts --- show reasonable agreement with those of numerical simulations at resolved scales. Our analytical model also enables self-consistent calculations of the boost factor of dark matter annhilation, which we find to increase from tens of percent at the smallest (Earth) and intermediate (dwarfs) masses to a factor of several at galaxy size, and to become as large as a factor of $\sim$10 for the largest halos (clusters) at small redshifts. Our analytical approach can accommodate substructures in the subhalos (sub-subhalos) in a consistent framework, which we find to give up to a factor of a few enhancement to the annihilation boost. Presence of the subhalos enhances the intensity of the isotropic gamma-ray background by a factor of a few, and as the result, the measurement by Fermi Large Area Telescope excludes the annihilation cross section greater than $\sim$$4\times 10^{-26}$ cm$^3$ s$^{-1}$ for dark matter masses up to $\sim$200 GeV.

We present the first individual and stacking systematic search for $\gamma$-ray emission in the GeV band with the {\it Fermi} Large Area Telescope (LAT), in the directions of 45 super-luminous supernovae (SLSNe). No excess of \gam rays from the SLSN positions was found. We report $\gamma$-ray luminosity upper limits and discuss the implication of these results on the origin of SLSNe in particular on the scenario of central compact object-aided SNe. From the stacking search, we derive an upper limit at 95\% confidence level (CL) to the $\gamma$-ray luminosity (above 600 MeV) $L_{\gamma}<9.1\times10^{41}$\,erg\,s$^{-1}$ for an assumed $E^{-2}$ photon spectrum, for our full SLSN sample. We conclude that the rate of the neutron stars born with millisecond rotation periods {$P\lesssim 2\,$ms and $B\sim10^{12-13}\,$G} must be lower than the rate of the observed SLSNe. The luminosity limits obtained on individual sources are also constraining: in particular, SN2013fc, CSS140222, SN2010kd, and PTF12dam can be born with millisecond periods only if $B\lesssim 10^{13}\,$G.

If dark energy (DE) is a dynamical field rather than a cosmological constant, an interaction between DE and the neutrino sector could exist, modifying the neutrino oscillation phenomenology and causing CP and apparent Lorentz violating effects. The terms in the Hamiltonian for flavor propagation induced by the DE-neutrino coupling do not depend on the neutrino energy, while the ordinary components decrease as $\Delta m^2/E_{\nu}$. Therefore, the DE-induced effects are absent at lower neutrino energies, but become significant at higher energies, allowing to be searched for by neutrino observatories. We explore the impact of the DE-neutrino coupling on the oscillation probability and the flavor transition in the three-flavor framework, and investigate the CP-violating and apparent Lorentz violating effects. We find that DE-induced effects become observable for $E_{\nu}m_{\text{eff}} \sim 10^{-20}~ \text{GeV}^2$, where $m_{\rm eff}$ is the effective mass parameter in the DE-induced oscillation probability, and CP is violated over a wide energy range. We also show that current and future experiments have the sensitivity to detect anomalous effects induced by a DE-neutrino coupling and probe the new mixing parameters. The DE-induced effects on neutrino oscillation can be distinguished from other new physics possibilities with similar effects, through the detection of the directional dependence of the interaction, which is specific to this interaction with DE. However, current experiments will not yet be able to measure the small changes of $\sim 0.03\%$ in the flavor composition due to this directional effect.

We present the measurement and interpretation of the angular power spectrum of nearby galaxies in the 2MASS Redshift Survey catalog with spectroscopic redshifts up to $z\approx 0.1$. We detect the angular power spectrum up to a multipole of $\ell\approx 1000$. We find that the measured power spectrum is dominated by galaxies living inside nearby galaxy clusters and groups. We use the halo occupation distribution (HOD) formalism to model the power spectrum, obtaining a fit with reasonable parameters. These HOD parameters are in agreement with the 2MASS galaxy distribution we measure toward the known nearby galaxy clusters, confirming validity of our analysis.

Dwarf spheroidals are low-luminosity satellite galaxies of the Milky Way highly dominated by dark matter (DM). Therefore, they are prime targets to search for signals from dark matter annihilation using gamma-ray observations. While the typical assumption is that the dark matter density profile of these satellite galaxies can be described by a spherical symmetric Navarro-Frenk-White (NFW) profile, recent observational data of stellar kinematics suggest that the DM halos around these galaxies are better described by axisymmetric profiles. Motivated by such evidence, we analyse about seven years of PASS8 Fermi data for seven classical dwarf galaxies, including Draco, adopting both the widely used NFW profile and observationally-motivated axisymmetric density profiles. For four of the selected dwarfs (Sextans, Carina, Sculptor and Fornax) axisymmetric mass models suggest a cored density profile rather than the commonly adopted cusped profile. We found that upper limits on the annihilation cross section for some of these dwarfs are significantly higher than the ones achieved using an NFW profile. Therefore, upper limits in the literature obtained using spherical symmetric cusped profiles, such as the NFW, might be overestimated. Our results show that it is extremely important to use observationally motivated density profiles going beyond the usually adopted NFW in order to obtain accurate constraints on the dark matter annihilation cross section.

Various particle physics models suggest that, besides the (nearly) cold dark matter that accounts for current observations, additional but sub-dominant dark relics might exist. These could be warm, hot, or even contribute as dark radiation. We present here a comprehensive study of two-component dark matter scenarios, where the first component is assumed to be cold, and the second is a non-cold thermal relic. Considering the cases where the non-cold dark matter species could be either a fermion or a boson, we derive consistent upper limits on the non-cold dark relic energy density for a very large range of velocity dispersions, covering the entire range from dark radiation to cold dark matter. To this end, we employ the latest Planck Cosmic Microwave Background data, the recent BOSS DR11 and other Baryon Acoustic Oscillation measurements, and also constraints on the number of Milky Way satellites, the latter of which provides a measure of the suppression of the matter power spectrum at the smallest scales due to the free-streaming of the non-cold dark matter component. We present the results on the fraction $f_{\rm ncdm}$ of non-cold dark matter with respect to the total dark matter for different ranges of the non-cold dark matter masses. We find that the 2$\sigma$ limits for non-cold dark matter particles with masses in the range 1--10 keV are $f_{\rm ncdm}\leq0.29$ (0.23) for fermions (bosons), and for masses in the 10--100 keV range they are $f_{\rm ncdm}\leq0.43$ (0.45), respectively.

After the discovery of extraterrestrial high-energy neutrinos, the next major goal of neutrino telescopes will be identifying astrophysical objects that produce them. The flux of the brightest source $F_{\rm max}$, however, cannot be probed by studying the diffuse neutrino intensity. We aim at constraining $F_{\rm max}$ by adopting a broken power-law flux distribution, a hypothesis supported by observed properties of any generic astrophysical sources. The first estimate of $F_{\rm max}$ comes from the fact that we can only observe one universe, and hence, the expected number of sources above $F_{\rm max}$ cannot be too small compared with one. For abundant source classes such as starburst galaxies, this one-source constraint yields a value of $F_{\rm max}$ that is an order of magnitude lower than the current upper limits from point-source searches. Then we derive upper limits on $F_{\rm max}$ assuming that the angular power spectrum is consistent with neutrino shot noise yet. We find that the limits obtained with upgoing muon neutrinos in IceCube can already be quite competitive, especially for rare but bright source populations such as blazars. The limits will improve nearly quadratically with exposure, and therefore be even more powerful for the next generation of neutrino telescopes.

We perform the first one-point fluctuation analysis of the high-energy neutrino sky. This method reveals itself to be especially suited to contemporary neutrino data, as it allows to study the properties of the astrophysical components of the high-energy flux detected by the IceCube telescope, even with low statistics and in the absence of point source detection. Besides the veto-passing atmospheric foregrounds, we adopt a simple model of the high-energy neutrino background by assuming two main extra-galactic components: star-forming galaxies and blazars. By leveraging multi-wavelength data from Herschel and Fermi, we predict the spectral and anisotropic probability distributions for their expected neutrino counts in IceCube. We find that star-forming galaxies are likely to remain a diffuse background due to the poor angular resolution of IceCube, and we determine an upper limit on the number of shower events that can reasonably be associated to blazars. We also find that upper limits on the contribution of blazars to the measured flux are unfavourably affected by the skewness of the blazar flux distribution. One-point event clustering and likelihood analyses of the IceCube HESE data suggest that this method has the potential to dramatically improve over more conventional model-based analyses, especially for the next generation of neutrino telescopes.

We measure the cross-correlation between Fermi-LAT gamma-ray photons and over 1000 deg$^2$ of weak lensing data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), the Red Cluster Sequence Lensing Survey (RCSLenS), and the Kilo Degree Survey (KiDS). We present the first measurement of tomographic weak lensing cross-correlations and the first application of spectral binning to cross-correlations between gamma rays and weak lensing. The measurements are performed using an angular power spectrum estimator while the covariance is estimated using an analytical prescription. We verify the accuracy of our covariance estimate by comparing it to two internal covariance estimators. Based on the non-detection of a cross-correlation signal, we derive constraints on weakly interacting massive particle (WIMP) dark matter. We compute exclusion limits on the dark matter annihilation cross-section $\langle\sigma_\rm{ann} v \rangle$, decay rate $\Gamma_\rm{dec}$, and particle mass $m_\rm{DM}$. We find that in the absence of a cross-correlation signal, tomography does not significantly improve the constraining power of the analysis. Assuming a strong contribution to the gamma-ray flux due to small-scale clustering of dark matter and accounting for known astrophysical sources of gamma rays, we exclude the thermal relic cross-section for particle masses of $m_\rm{DM}\lesssim 20$ GeV.

Recently, a new measurement of the auto- and cross-correlation angular power spectrum (APS) of the isotropic gamma-ray background was performed, based on 81 months of data of the Fermi Large-Area Telescope (LAT). Here, we fit, for the first time, the new APS data with a model describing the emission of unresolved blazars. These sources are expected to dominate the anisotropy signal. The model we employ in our analysis reproduces well the blazars resolved by Fermi LAT. When considering the APS obtained by masking the sources in the 3FGL catalogue, we find that unresolved blazars under-produce the measured APS below $\sim$1 GeV. Contrary to past results, this suggests the presence of a new contribution to the low-energy APS, with a significance of, at least, 5$\sigma$. The excess can be ascribed to a new class of faint gamma-ray emitters. If we consider the APS obtained by masking the sources in the 2FGL catalogue, there is no under-production of the APS below 1 GeV, but the new source class is still preferred over the blazars-only scenario (with a significance larger than 10$\sigma$). The properties of the new source class and the level of anisotropies induced in the isotropic gamma-ray background are the same, independent of the APS data used. In particular, the new gamma-ray emitters must have a soft energy spectrum, with a spectral index ranging, approximately, from 2.7 to 3.2. This complicates their interpretation in terms of known sources, since, normally, star-forming and radio galaxies are observed with a harder spectrum. The new source class identified here is also expected to contribute significantly to the intensity of the isotropic gamma-ray background.

It has been long discussed that cosmic rays may contain signals of dark matter. In the last couple of years an anomaly of cosmic-ray positrons has drawn a lot of attentions, and recently an excess in cosmic-ray anti-proton has been reported by AMS-02 collaboration. Both excesses may indicate towards decaying or annihilating dark matter with a mass of around 1-10 TeV. In this article we study the gamma rays from dark matter and constraints from cross correlations with distribution of galaxies, particularly in a local volume. We find that gamma rays due to inverse-Compton process have large intensity, and hence they give stringent constraints on dark matter scenarios in the TeV scale mass regime. Taking the recent developments in modeling astrophysical gamma-ray sources as well as comprehensive possibilities of the final state products of dark matter decay or annihilation into account, we show that the parameter regions of decaying dark matter that are suggested to explain the excesses are excluded. We also discuss the constrains on annihilating scenarios.

Long-duration gamma-ray bursts (GRBs) have been often considered as the natural evolution of some core-collapse supernovae (SNe). While GRBs with relativistic jets emit an electromagnetic signal, GRBs with mildly relativistic jets are opaque to photons and, therefore, could be detectable through neutrinos only. We discuss the possibility that successful GRBs and mildly relativistic jets belong to the same class of astrophysical transients with different Lorentz factor Gamma_b and study the production of high-energy neutrinos as a function of Gamma_b, by including both proton-photon and proton-proton interactions. By assuming a SN-GRB connection, we find that the diffuse neutrino emission from optically thick jets with Lorentz factors lower than the ones of successful GRBs can be one of the main components of the observed IceCube high-energy neutrino flux. Moreover, under the assumption that all these jets belong to the same class of astrophysical transients, we show that the IceCube high-energy neutrino data provide indirect constraints on the rate of non-successful jets, favoring a local rate lower than tens of percent of the local SN rate. These limits are currently comparable to the ones obtained in dedicated searches on choked sources and are expected to become tighter with accumulation of more high-energy neutrino data.

Mounting evidence suggests that the TeV-PeV neutrino flux detected by the IceCube telescope has mainly an extragalactic origin. If such neutrinos are primarily produced by a single class of astrophysical sources via hadronuclear ($pp$) interactions, a similar flux of gamma-ray photons is expected. For the first time, we employ tomographic constraints to pinpoint the origin of the IceCube neutrino events by analyzing recent measurements of the cross correlation between the distribution of GeV gamma rays, detected by the Fermi satellite, and several galaxy catalogs in different redshift ranges. We find that the corresponding bounds on the neutrino luminosity density are up to one order of magnitude tighter than those obtained by using only the spectrum of the gamma-ray background, especially for sources with mild redshift evolution. In particular, our method excludes any hadronuclear source with a spectrum softer than $E^{-2.1}$ as a main component of the neutrino background, if its evolution is slower than $(1+z)^3$. Starburst galaxies, if able to accelerate and confine cosmic rays efficiently, satisfy both spectral and tomographic constraints.

In the cold dark matter paradigm, structures form hierarchically, implying that large structures contain smaller substructures. These subhalos will enhance signatures of dark matter annihilation such as gamma rays. In the literature, typical estimates of this boost factor assume a concentration-mass relation for field halos, to calculate the luminosity of subhalos. However, since subhalos accreted in the gravitational potential of their host loose mass through tidal stripping and dynamical friction, they have a quite characteristic density profile, different from that of the field halos of the same mass. In this work, we quantify the effect of tidal stripping on the boost factor, by developing a semi-analytic model that combines mass-accretion history of both the host and subhalos as well as subhalo accretion rates. We find that when subhalo luminosities are treated consistently, the boost factor increases by a factor 2-5, compared to the typical calculation assuming a field-halo concentration. This holds for host halos ranging from sub-galaxy to cluster masses and is independent of the subhalo mass function or specific concentration-mass relation. The results are particularly relevant for indirect dark matter searches in the extragalactic gamma-ray sky.

We investigate the prospects of indirect and direct dark matter searches within the minimal supersymmetric standard model with nine parameters (MSSM-9). These nine parameters include three gaugino masses, Higgs, slepton and squark masses, all treated independently. We perform a Bayesian Monte Carlo scan of the parameter space taking into consideration all available particle physics constraints such as the Higgs mass of 126 GeV, upper limits on the scattering cross-section from direct-detection experiments, and assuming that the MSSM-9 provides all the dark matter abundance through thermal freeze-out mechanism. Within this framework we find two most probable regions for dark matter: 1-TeV higgsino-like and 3-TeV wino-like neutralinos. We discuss prospects for future indirect (in particular the Cherenkov Telescope Array, CTA) and direct detection experiments. We find that for slightly contracted dark matter profiles in our Galaxy, which can be caused by the effects of baryonic infall in the Galactic center, CTA will be able to probe a large fraction of the remaining allowed region in synergy with future direct detection experiments like XENON-1T.

We study how the kinetic decoupling of dark matter (DM) within a minimal supersymmetric extension of the standard model, by adopting nine independent parameters (MSSM-9), could improve our knowledge of the properties of the DM protohalos. We show that the most probable neutralino mass regions, which satisfy the relic density and the Higgs mass contraints, are those with the lightest supersymmetric neutralino mass around 1 TeV and 3 TeV, corresponding to Higgsino-like and Wino-like neutralino, respectively. The kinetic decoupling temperature in the MSSM-9 scenario leads to a most probable protohalo mass in a range of $M_{\mathrm{ph}}\sim 10^{-12}-10^{-7}\,M_\odot$. The part of the region closer to 2 TeV gives also important contributions from the neutralino-stau co-annihilation, reducing the effective annihilation rate in the early Universe. We also study how the size of the smallest DM substructures correlates to experimental signatures, such as the spin-dependent and spin-independent scattering cross sections, relevant for direct detection of DM. Improvements on the spin-independent sensitivity might reduce the most probable range of the protohalo mass between $\sim$10$^{-9}\,M_\odot$ and $\sim$10$^{-7}\,M_\odot$, while the expected spin-dependent sensitivity provides weaker constraints. We show how the boost of the luminosity due to DM annihilation increases, depending on the protohalo mass. In the Higgsino case, the protohalo mass is lower than the canonical value often used in the literature ($\sim$10$^{-6}\,M_\odot$), while $\langle\sigma v\rangle$ does not deviate from $\langle\sigma v\rangle\sim 10^{-26}$ cm$^3$ s$^{-1}$; there is no significant enhancement of the luminosity. On the contrary, in the Wino case, the protohalo mass is even lighter, and $\langle\sigma v\rangle$ is two orders of magnitude larger; as its consequence, we see a substantial enhancement of the luminosity.

We investigate for the first time the potential of angular auto- and cross-correlation power spectra in identifying sterile neutrino dark matter in the cosmic X-ray background. We take as reference the performance of the soon-to-be-launched eROSITA satellite. The main astrophysical background sources against sterile neutrino decays are active galactic nuclei, galaxies powered by X-ray binaries, and clusters of galaxies. While sterile neutrino decays are always subdominant in the auto-correlation power spectra, they can be efficiently enhanced when cross-correlating with tracers of the dark matter distribution such as galaxies in the 2MASS catalogues. We show that the planned four-years eROSITA all-sky survey will provide a large enough photon statistics to potentially yield very stringent constraints on the decay lifetime, enabling to firmly test the recently claimed 3.56-keV X-ray line found towards several clusters and galaxies and its decaying dark matter interpretation. However, we also show that in order to fully exploit the potential of eROSITA for dark matter searches, it is vital to overcome the shot-noise limitations inherent to galaxy catalogues as tracers for the dark matter distribution.

Cosmic-ray (CR) protons can accumulate for cosmological times in clusters of galaxies. Their hadronic interactions with protons of the intra-cluster medium (ICM) generate secondary electrons, gamma-rays and high-energy neutrinos. In light of the high-energy neutrino events recently discovered by the IceCube observatory, we estimate the contribution from galaxy clusters to the diffuse gamma-ray and neutrino backgrounds. For the first time, we consistently take into account the synchrotron emission generated by secondary electrons and require the clusters radio counts to be respected. For a choice of parameters respecting current constraints from radio to gamma-rays, and assuming a proton spectral index of -2, we find that hadronic interactions in clusters contribute by less than 10% to the IceCube flux, and much less to the total extragalactic gamma-ray background observed by Fermi. They account for less than 1% for spectral indexes <-2. The high-energy neutrino flux observed by IceCube can be reproduced without violating radio constraints only if a very hard (and speculative) spectral index >-2 is adopted. However, this scenario is in tension with the high-energy IceCube data, which seem to suggest a spectral energy distribution of the neutrino flux that decreases with the particle energy. We stress that our results are valid for all kind of sources injecting CR protons into the ICM, and that, while IceCube can test the most optimistic scenarios for spectral indexes >=-2.2 by stacking few nearby massive objects, clusters of galaxies cannot give any relevant contribution to the extragalactic gamma-ray and neutrino backgrounds in any realistic scenario.

Gamma-ray bursts (GRBs) have been suggested as possible sources of the high-energy neutrino flux recently detected by the IceCube telescope. We revisit the fireball emission model and elaborate an analytical prescription to estimate the high-energy neutrino prompt emission from pion and kaon decays, assuming that the leading mechanism for the neutrino production is lepto-hadronic. To this purpose, we include hadronic, radiative and adiabatic cooling effects and discuss their relevance for long- (including high- and low-luminosity) and short-duration GRBs. The expected diffuse neutrino background is derived, by requiring that the GRB high-energy neutrino counterparts follow up-to-date gamma-ray luminosity functions and redshift evolutions of the long and short GRBs. Although dedicated stacking searches have been unsuccessful up to now, we find that GRBs could contribute up to a few % to the observed IceCube high-energy neutrino flux for sub-PeV energies, assuming that the latter has a diffuse origin. Gamma-ray bursts, especially low-luminosity ones, could however be the main sources of the IceCube high-energy neutrino flux in the PeV range. While high-luminosity and low-luminosity GRBs have comparable intensities, the contribution from the short-duration component is significantly smaller. Our findings confirm the most-recent IceCube results on the GRB searches and suggest that larger exposure is mandatory to detect high-energy neutrinos from high-luminosity GRBs in the near future.

Late time decay of very heavy dark matter is considered as one of the possible explanations for diffuse PeV neutrinos observed in IceCube. We consider implications of multimessenger constraints, and show that proposed models are marginally consistent with the diffuse gamma-ray background data. Critical tests are possible by a detailed analysis and identification of the sub-TeV isotropic diffuse gamma-ray data observed by Fermi and future observations of sub-PeV gamma rays by observatories like HAWC or Tibet AS+MD. In addition, with several-year observations by next-generation telescopes such as IceCube-Gen2, muon neutrino searches for nearby dark matter halos such as the Virgo cluster should allow us to rule out or support the dark matter models, independently of gamma-ray and anisotropy tests.

The one-point function (i.e., the isotropic flux distribution) is a complementary method to (anisotropic) two-point correlations in searches for a gamma-ray dark matter annihilation signature. Using analytical models of structure formation and dark matter halo properties, we compute the gamma-ray flux distribution due to annihilations in extragalactic dark matter halos, as it would be observed by the Fermi Large Area Telescope. Combining the central limit theorem and Monte Carlo sampling, we show that the flux distribution takes the form of a narrow Gaussian of `diffuse' light, with an `unresolved point source' power-law tail as a result of bright halos. We argue that this background due to dark matter constitutes an irreducible and significant background component for point-source annihilation searches with galaxy clusters and dwarf spheroidal galaxies, modifying the predicted signal-to-noise ratio. A study of astrophysical backgrounds to this signal reveals that the shape of the total gamma-ray flux distribution is very sensitive to the contribution of a dark matter component, allowing us to forecast promising one-point upper limits on the annihilation cross section. We show that by using the flux distribution at only one energy bin, one can probe the canonical cross section required for explaining the relic density, for dark matter of masses around tens of GeV.

The intra-cluster medium of several galaxy clusters hosts large-scale regions of diffuse synchrotron radio emission, known as radio halos and relics, which demonstrate the presence of magnetic fields and relativistic electrons in clusters. These relativistic electrons should also emit X-rays through inverse-Compton scattering off of cosmic microwave background photons. The detection of such a non-thermal X-ray component, together with the radio measurement, would permit to clearly separate the magnetic field from the relativistic electron distribution as the inverse-Compton emission is independent from the magnetic field in the cluster. However, non-thermal X-rays have not been conclusively detected from any cluster of galaxies so far. In this paper, for the first time, we model the synchrotron and inverse-Compton emission of all clusters hosting radio halos and relics for which the spectral index can be determined. We provide constraints on the volume-average magnetic field by comparing with current X-ray measurements. We then estimate the maximum volume-average magnetic field that will allow the detection of inverse-Compton hard X-rays by the ASTRO-H satellite. We found that several clusters are good targets for ASTRO-H to detect their inverse-Compton emission, in particular that corresponding to radio relics, and propose a list of promising targets for which ASTRO-H can test $\ge1$~$\mu$G magnetic fields. We conclude that future hard X-ray observations by the already-operating NuSTAR and the soon-to-be-launched ASTRO-H definitely have the potential to shed light on the long-sought non-thermal hard-X-ray emission in clusters of galaxies.

Cosmic-ray protons accumulate for cosmological times in clusters of galaxies as their typical radiative and diffusive escape times are longer than the Hubble time. Their hadronic interactions with protons of the intra-cluster medium generate secondary electrons, gamma-rays and neutrinos. We here estimate the contribution from clusters to the diffuse gamma-ray and neutrino backgrounds. We model the cluster population by means of their mass function, using a phenomenological luminosity-mass relation applied to all clusters, as well as a detailed semi-analytical model. Additionally, we consider observationally-motivated values for the cluster magnetic field. This is a crucial parameter since the observed radio counts due to synchrotron emission by secondary electrons need to be respected. For a choice of parameters respecting all current constraints, and assuming a spectral index of -2, we find that hadronic interactions in clusters contribute by less than 10% to the extragalactic gamma-ray background observed by Fermi and to the IceCube flux. They account for less than 1% for spectral indices <=-2. The neutrino flux observed by IceCube can be reproduced without violating constraints only if a very hard (and speculative) spectral index >-2 is adopted. However, this scenario is in tension with the IceCube data, which seem to suggest a spectral energy distribution of the neutrino flux that decreases with energy. In the case of proton-photon interactions, we find that very likely protons do not reach sufficiently high energies to produce neutrinos in clusters. We argue that our results are optimistic due to our assumptions, and that clusters of galaxies cannot give any relevant contribution to the extragalactic gamma-ray and neutrino backgrounds. Finally, we find that the cluster contribution to the angular fluctuations in the gamma-ray background is subdominant, less than 10%. [abridged]

If dark matter is unstable and the mass is within GeV-TeV regime, its decays produce high-energy photons that give contribution to the extragalactic gamma-ray background (EGRB). We constrain dark matter decay by analyzing the 50-month EGRB data measured with Fermi satellite, for different decay channels motivated with several supersymmetric scenarios featuring R-parity violation. We adopt the latest astrophysical models for various source classes such as active galactic nuclei and star-forming galaxies, and take associated uncertainties properly into account. The lower limits for the lifetime are very stringent for a wide range of dark matter mass, excluding the lifetime shorter than 10^28 s for mass between a few hundred GeV and ~1TeV, e.g., for b\bar{b} decay channel. Furthermore, most dark matter models that explain the anomalous positron excess are also excluded. These constraints are robust, being little dependent on astrophysical uncertainties, unlike other probes such as Galactic positrons or anti-protons.

Cross-correlating the gamma-ray background with local galaxy catalogs potentially gives stringent constraints on dark matter annihilation. We provide updated theoretical estimates of sensitivities to the annihilation cross section from gamma-ray data with Fermi telescope and 2MASS galaxy catalogs, by elaborating the galaxy power spectrum and astrophysical backgrounds, and adopting the Markov-Chain Monte Carlo simulations. In particular, we show that taking tomographic approach by dividing the galaxy catalogs into more than one redshift slice will improve the sensitivity by a factor of a few to several. If dark matter halos contain lots of bright substructures, yielding a large annihilation boost (e.g., a factor of $\sim$100 for galaxy-size halos), then one may be able to probe the canonical annihilation cross section for thermal production mechanism up to masses of $\sim$700 GeV. Even with modest substructure boost (e.g., a factor of $\sim$10 for galaxy-size halos), on the other hand, the sensitivities could still reach a factor of three larger than the canonical cross section for dark matter masses of tens to a few hundreds of GeV.

Star-forming galaxies have been predicted to contribute considerably to the diffuse gamma-ray background as they are guaranteed reservoirs of cosmic rays. Assuming that the hadronic interactions responsible for high-energy gamma rays also produce high-energy neutrinos and that O(100) PeV cosmic rays can be produced and confined in starburst galaxies, we here discuss the possibility that star-forming galaxies are also the main sources of the high-energy neutrinos observed by the IceCube experiment. First, we compute the diffuse gamma-ray background from star-forming galaxies, adopting the latest Herschel PEP/HerMES luminosity function and relying on the correlation between the gamma-ray and infrared luminosities reported by Fermi observations. Then we derive the expected intensity of the diffuse high-energy neutrinos from star-forming galaxies including normal and starburst galaxies. Our results indicate that starbursts, including those with active galactic nuclei and galaxy mergers, could be the main sources of the high-energy neutrinos observed by the IceCube experiment. We find that assuming a cosmic-ray spectral index of 2.1-2.2 for all starburst-like galaxies, our predictions can be consistent with both the Fermi and IceCube data, but larger indices readily fail to explain the observed diffuse neutrino flux. Taking the starburst high-energy spectral index as free parameter, and extrapolating from GeV to PeV energies, we find that the spectra harder than E^(-2.15) are likely to be excluded by the IceCube data, which can be more constraining than the Fermi data for this population.

Cross-correlating gamma-ray maps with locations of galaxies in the low-redshift Universe vastly increases sensitivity to signatures of annihilation of dark matter particles. Low-redshift galaxies are ideal targets, as the largest contribution to anisotropy in the gamma-ray sky from annihilation comes from $z\lesssim 0.1$, where we expect minimal contributions from astrophysical sources such as blazars. Cross-correlating the five-year data of Fermi-LAT with the redshift catalog of the 2MASS survey can detect gamma rays from annihilation if dark matter has the canonical annihilation cross section and its mass is smaller than $\sim$100 GeV.

We analyze 5-year (63 months) data of the Large Area Telescope on board Fermi satellite from the Coma galaxy cluster in the energy range between 100 MeV and 100 GeV. The likelihood analyses are performed with several templates motivated by models predicting gamma-ray emission due to structure formation processes. We find no excess emission and derive the most stringent constraints to date on the Coma cluster above 100 MeV, and on the tested scenarios in general. The upper limits on the integral flux range from 10^-10 to 10^-9 cm^-2s^-1, and are stringent enough to challenge different scenarios. We find that the acceleration efficiency of cosmic ray protons and electrons at shocks must be below approximately 15% and 1%, respectively. Additionally, we argue that the proton acceleration efficiency should be lower than 5% in order to be consistent with radio data. This, however, relays on magnetic field estimates in the cluster. In particular, this implies that the contribution to the diffuse extragalactic gamma-ray background due to gamma-rays from structure formation processes in clusters of galaxies is negligible, below 1%. Finally, we discuss future detectability prospects for Astro-H, Fermi after 10-yr of operation, and the Cherenkov Telescope Array.

Many of the astrophysical sources and violent phenomena observed in our Universe are potential emitters of gravitational waves (GW) and high-energy neutrinos (HEN). Both GWs and HENs may escape very dense media and travel unaffected over cosmological distances, carrying information from the innermost regions of the astrophysical engines. Such messengers could also reveal new, hidden sources that have not been observed by conventional photon-based astronomy. Coincident observation of GWs and HENs may thus play a critical role in multimessenger astronomy. This is particularly true at the present time owing to the advent of a new generation of dedicated detectors: IceCube, ANTARES, VIRGO and LIGO. Given the complexity of the instruments, a successful joint analysis of this data set will be possible only if the expertise and knowledge of the data is shared between the two communities. This review aims at providing an overview of both theoretical and experimental state-of-the-art and perspectives for such a GW+HEN multimessenger astronomy.

Annihilation of dark matter particles in cosmological halos (including a halo of the Milky Way) contributes to the diffuse gamma-ray background (DGRB). As this contribution will appear anisotropic in the sky, one can use the angular power spectrum of anisotropies in DGRB to constrain properties of dark matter particles. By comparing the updated analytic model of the angular power spectrum of DGRB from dark matter annihilation with the power spectrum recently measured from the 22-month data of Fermi Large Area Telescope (LAT), we place upper limits on the annihilation cross section of dark matter particles as a function of dark matter masses. We find that the current data exclude <\sigma v> >~ 10^{-25} cm^3 s^{-1} for annihilation into b\bar{b} at the dark matter mass of 10 GeV, which is a factor of three times larger than the canonical cross section. The limits are weaker for larger dark matter masses. The limits can be improved further with more Fermi-LAT data as well as by using the power spectrum at lower multipoles (l <~ 150), which are currently not used due to a potential Galactic foreground contamination.

We analyze 2.8-yr data of 1-100 GeV photons for clusters of galaxies, collected with the Large Area Telescope onboard the Fermi satellite. By analyzing 49 nearby massive clusters located at high Galactic latitudes, we find no excess gamma-ray emission towards directions of the galaxy clusters. Using flux upper limits, we show that the Fornax cluster provides the most stringent constraints on the dark matter annihilation cross section. Stacking a large sample of nearby clusters does not help improve the limit for most dark matter models. This suggests that a detailed modeling of the Fornax cluster is important for setting robust limits on the dark matter annihilation cross section based on clusters. We therefore perform the detailed mass modeling and predict the expected dark matter annihilation signals from the Fornax cluster, by taking into account effects of dark matter contraction and substructures. By modeling the mass distribution of baryons (stars and gas) around a central bright elliptical galaxy, NGC 1399, and using a modified contraction model motivated by numerical simulations, we show that the dark matter contraction boosts the annihilation signatures by a factor of 4. For dark matter masses around 10 GeV, the upper limit obtained on the annihilation cross section times relative velocity is <\sigma v> <~ (2-3)x10^{-25} cm^3 s^{-1}, which is within a factor of 10 from the value required to explain the dark matter relic density. This effect is more robust than the annihilation boost due to substructure, and it is more important unless the mass of the smallest subhalos is much smaller than that of the Sun.

We report a measurement of intergalactic magnetic fields using combined data from Atmospheric Cherenkov Telescopes and Fermi Gamma-Ray Space Telescope, based on the spectral data alone. If blazars are assumed to produce both gamma rays and cosmic rays, the observed spectra are not sensitive to the intrinsic spectrum of the source, because, for a distant blazar, secondary photons produced in line-of-sight cosmic-ray interactions dominate the signal. In this case, we find 0.01 fG < B < 30 fG. If one excludes the cosmic-ray component, the 0.01 fG lower limit remains, but the upper limit depends on the spectral properties of the source. We present the allowed ranges for a variety of model parameters.

The legacy of solar neutrinos suggests that large neutrino detectors should be sited underground. However, to instead go underwater bypasses the need to move mountains, allowing much larger water Cherenkov detectors. We show that reaching a detector mass scale of ~5 Megatons, the size of the proposed Deep-TITAND, would permit observations of neutrino "mini-bursts" from supernovae in nearby galaxies on a roughly yearly basis, and we develop the immediate qualitative and quantitative consequences. Importantly, these mini-bursts would be detected over backgrounds without the need for optical evidence of the supernova, guaranteeing the beginning of time-domain MeV neutrino astronomy. The ability to identify, to the second, every core collapse in the local Universe would allow a continuous "death watch" of all stars within ~5 Mpc, making practical many previously-impossible tasks in probing rare outcomes and refining coordination of multi-wavelength/multi-particle observations and analysis. These include the abilities to promptly detect otherwise-invisible prompt black hole formation, provide advance warning for supernova shock-breakout searches, define tight time windows for gravitational-wave searches, and identify "supernova impostors" by the non-detection of neutrinos. Observations of many supernovae, even with low numbers of detected neutrinos, will help answer questions about supernovae that cannot be resolved with a single high-statistics event in the Milky Way.

Intergalactic magnetic fields (IGMF) can cause the appearance of halos around the gamma-ray images of distant objects because an electromagnetic cascade initiated by a high-energy gamma-ray interaction with the photon background is broadened by magnetic deflections. We report evidence of such gamma-ray halos in the stacked images of the 170 brightest active galactic nuclei (AGN) in the 11-month source catalog of the Fermi Gamma-Ray Space Telescope. Excess over point spread function in the surface brightness profile is statistically significant at 3.5\sigma (99.95% confidence level), for the nearby, hard population of AGN. The halo size and brightness are consistent with IGMF, B_{IGMF} ~ 10^{-15} G. The knowledge of IGMF will facilitate the future gamma-ray and charged-particle astronomy. Furthermore, since IGMF are likely to originate from the primordial seed fields created shortly after the Big Bang, this potentially opens a new window on the origin of cosmological magnetic fields, inflation, and the phase transitions in the early Universe.

If dark energy (DE) couples to neutrinos, then there may be apparent violations of Lorentz/CPT invariance in neutrino oscillations. The DE-induced Lorentz/CPT violation takes a specific form that introduces neutrino oscillations that are energy independent, differ for particles and antiparticles, and can lead to novel effects for neutrinos propagating through matter. We show that ultra-high-energy neutrinos may provide one avenue to seek this type of Lorentz/CPT violation in \nu_\mu-\nu_\tau oscillations, improving the current sensitivity to such effects by seven orders of magnitude. Lorentz/CPT violation in electron-neutrino oscillations may be probed with the zenith-angle dependence for high-energy atmospheric neutrinos. The ``smoking gun,'' for DE-neutrino coupling would, however, be a dependence of neutrino oscillations on the direction of the neutrino momentum relative to our peculiar velocity with respect to the CMB rest frame. While the amplitude of this directional dependence is expected to be small, it may nevertheless be worth seeking in current data and may be a target for future neutrino experiments.

We present analytic approximations to the optically thin synchrotron and synchrotron self-Compton (SSC) spectra when Klein-Nishina (KN) effects are important and pair production and external radiation fields can be neglected. This theory is useful for analytical treatment of radiation from astrophysical sources, such as gamma-ray bursts (GRBs), active galactic nuclei and pulsar wind nebula, where KN effects may be important. We consider a source with a continuous injection of relativistic electrons with a power-law energy distribution above some typical injection energy. We find that the synchrotron-SSC spectra can be described by a broken power-law, and provide analytic estimates for the break frequencies and power-law indices. In general, we show that the dependence of the KN cross-section on the energy of the upscattering electron results in a hardening of the energy distribution of fast cooling electrons and therefore in a hardening of the observed synchrotron spectrum. As a result the synchrotron spectrum of fast cooling electrons, below the typical injection energy, can be as hard as $F_\nu \propto \nu^0$, instead of the classical $\nu^{-1/2}$ when KN effects are neglected. The synchrotron energy output can be dominated by electrons with energy above the typical injection energy. We solve self-consistently for the cooling frequency and find that the transition between synchrotron and SSC cooling can result in a discontinuous variations of the cooling frequency and the synchrotron and SSC spectra. We demonstrate the application of our results to theory by applying them to prompt and afterglow emission models of GRBs.

Star-forming galaxies are predicted to contribute considerably to the cosmic gamma-ray background (CGB) as they are confirmed gamma-ray emitters and are the most numerous population of gamma-ray sources, although individually faint. Even though the Fermi Gamma-ray Space Telescope will be able to resolve few star-forming galaxies individually, their fractional contribution to the CGB should become far more significant than it was for past measurements of the CGB as many of the brighter, formerly unresolved sources will be resolved out. Thus, the clustering feature of galaxies imprinted on the CGB might be detectable by Fermi. In anticipation of such measurements, we calculate the predicted angular auto-power and cross-power spectra of the CGB from normal galaxies. We find that the amplitude of the auto-power spectrum is smaller than that for other sources such as blazars and dark-matter annihilation; the shape is also characteristic. We also show that the cross-power spectrum with galaxy surveys features larger amplitude. Fermi should be able to detect the correlation signature in both the auto-power and cross-power spectra at angular scales of ~1-10 deg after 5-yr of operation. Such a detection would be valuable in confirming the level of the star-forming galaxy contribution to the CGB, and more importantly, in serving as a tool in the effort to discriminate between possible origins of the CGB.

Dark matter annihilation in galactic substructure would imprint characteristic angular signatures on the all-sky map of the diffuse gamma-ray background. We study the gamma-ray background anisotropy due to the subhalos and discuss detectability at Fermi Gamma-ray Space Telescope. We derive analytic formulae that enable to directly compute the angular power spectrum, given parameters of subhalos. As our fiducial subhalo models, we adopt M^{-1.9} mass spectrum, subhalos radial distribution suppressed towards the galactic center, and luminosity profile of each subhalo dominated by its smooth component. We find that, for multipole regime corresponding to \theta <~ 5 deg, the angular power spectrum is dominated by a noise-like term, with suppression due to internal structure of relevant subhalos. If the mass spectrum extends down to Earth-mass scale, then the subhalos would be detected in the anisotropy with Fermi at angular scales of ~10 deg, if their contribution to the gamma-ray background is larger than ~20%. If the minimum mass is around 10^4 M_{sun}, on the other hand, the relevant angular scale for detection is ~1 deg, and the anisotropy detection requires that the subhalo contribution to the gamma-ray background intensity is only ~4%. These can be achieved with a modest boost for particle physics parameters. We also find that the anisotropy analysis could be a more sensitive probe for the subhalos than individual detection. We also study dependence on model parameters, where we reach the similar conclusions for all the models investigated. The analytic approach should be very useful when Fermi data are analyzed and the obtained angular power spectrum is interpreted in terms of subhalo models.

One of the targets of the recently launched Fermi Gamma-ray Space Telescope is a diffuse gamma-ray background from dark-matter annihilation or decay in the Galactic halo. N-body simulations and theoretical arguments suggest that the dark matter in the Galactic halo may be clumped into substructure, rather than smoothly distributed. Here we propose the gamma-ray-flux probability distribution function (PDF) as a probe of substructure in the Galactic halo. We calculate this PDF for a phenomenological model of halo substructure and determine the regions of the substructure parameter space in which the PDF may be distinguished from the PDF for a smooth distribution of dark matter. In principle, the PDF allows a statistical detection of substructure, even if individual halos cannot be detected. It may also allow detection of substructure on the smallest microhalo mass scales, $\sim M_\oplus$, for weakly-interacting massive particles (WIMPs). Furthermore, it may also provide a method to measure the substructure mass function. However, an analysis that assumes a typical halo substructure model and a conservative estimate of the diffuse background suggests that the substructure PDF may not be detectable in the lifespan of Fermi in the specific case that the WIMP is a neutralino. Nevertheless, for a large range of substructure, WIMP annihilation, and diffuse background models, PDF analysis may provide a clear signature of substructure.

Many of the astrophysical sources and violent phenomena observed in our Universe are potential emitters of gravitational waves (GWs) and high-energy neutrinos (HENs). A network of GW detectors such as LIGO and Virgo can determine the direction/time of GW bursts while the IceCube and ANTARES neutrino telescopes can also provide accurate directional information for HEN events. Requiring the consistency between both, totally independent, detection channels shall enable new searches for cosmic events arriving from potential common sources, of which many extra-galactic objects.

Dark matter (DM) annihilation could in principle contribute to the diffuse cosmic gamma-ray back- ground (CGB). While with standard assumptions for cosmological and particle physics parameters this contribution is expected to be rather small, a number of processes could boost it, including a larger-than-expected DM annihilation cross-section, or the occurance of DM substructures such as DM mini-spikes around intermediate-mass black holes. We show that angular correlations of the CGB provide a tool to disentangle the signal induced by DM annihilation in mini-spikes from a conventional astrophysical component. Treating blazars as a known background, we study the prospects for detecting DM annihilations with the Fermi Gamma-Ray Space Telescope for different choices of DM mass and annihilation channels.

Under the fossil field hypothesis of the origin of magnetar magnetic fields, the magnetar inherits its magnetic field from its progenitor. We show that during the supernova of such a progenitor, protons may be accelerated to \sim 10^4 GeV as the supernova shock propagates in the magnetic stellar envelope. Inelastic nuclear collisions of these protons produce a flash of high-energy neutrinos arriving a few hours after thermal (10 MeV) neutrinos. The neutrino flash is characterized by energies up to O(100) GeV and durations seconds to hours, depending on the progenitor: those from smaller Type Ibc progenitors are typically shorter in duration and reach higher energies compared to those from larger Type II progenitors. A Galactic Type Ib supernova leaving behind a magnetar remnant will yield up to \sim 160 neutrino induced muon events in Super-Kamiokande, and up to \sim 7000 in a km^3 class detector such as IceCube, providing a means of probing supernova models and the presence of strong magnetic fields in the stellar envelope.

One of the goals of NASA's Fermi Gamma-ray Space Telescope (formerly GLAST) will be detection of gamma rays from dark-matter annihilation in the Galactic halo. Theoretical arguments suggest that dark matter may be bound into subhalos with masses as small as 10^{-4}--10^2 Earth mass. If so, it may be possible to detect individual subhalos as point sources in the Fermi Telescope. It has further been argued that some of these point sources may exhibit proper motions. Here we show that upper limits to the diffuse gamma-ray background constrain the number of subhalos close enough to exhibit proper motions to be less than one.

A class of very energetic supernovae (hypernovae) is associated with long gamma-ray bursts, in particular with a less energetic but more frequent population of gamma-ray bursts. Hypernovae also appear to be associated with mildly relativistic jets or outflows, even in the absence of gamma-ray bursts. Here we consider radiation from charged particles accelerated in such mildly relativistic outflows with kinetic energies of ~10^{50} erg. The radiation processes of the primarily accelerated electrons considered are synchrotron radiation and inverse-Compton scattering of synchrotron photons (synchrotron self-Compton; SSC) and of supernova photons (external inverse-Compton; EIC). In the soft X-ray regime, both the SSC and EIC flux can be the dominant component, but due to their very different spectral shapes it should be easy to distinguish between them. When the fraction of the kinetic energy going into the electrons (\epsilon_e) is large, the SSC is expected to be important; otherwise the EIC will dominate. The EIC flux is quite high, almost independently of \epsilon_e, providing a good target for X-ray telescopes such as XMM-Newton and Chandra. In the GeV gamma-ray regime, the EIC would be the dominant radiation process and the Gamma-ray Large Area Space Telescope (GLAST) should be able to probe the value of \epsilon_e, the spectrum of the electrons, and their maximum acceleration energy. Accelerated protons also lead to photon radiation through the secondary electrons produced by the photopion and photopair processes. We find that over a significant range of parameters the proton component is generally less prominent than the primary electron component. We discuss the prospects for the detection of the X-ray and GeV signatures of the mildly relativistic outflow of hypernovae.

We investigate the GeV emission from gamma-ray bursts (GRBs), using the results from the Energetic Gamma Ray Experimental Telescope (EGRET), and in view of the Gamma-ray Large Area Space Telescope (GLAST). Assuming that the conventional prompt and afterglow photons originate from synchrotron radiation, we compare an accompanying inverse-Compton component with EGRET measurements and upper limits on GeV fluence, taking Klein-Nishina feedback into account. We find that EGRET constraints are consistent with the theoretical framework of the synchrotron self-Compton model for both prompt and afterglow phases, and discuss constraints on microphysical parameters in both phases. Based on the inverse-Compton model and using EGRET results, we predict that GLAST would detect GRBs with GeV photons at a rate >~20 yr^{-1} from each of the prompt and afterglow phases. This rate applies to the high-energy tail of the prompt synchrotron emission and to the inverse-Compton component of the afterglow. Theory predicts that in a large fraction of the cases where synchrotron GeV prompt emission would be detected by GLAST, inverse-Compton photons should be detected as well at high energies >~10 GeV. Therefore GLAST will enable a more precise test of the high-energy emission mechanism. Finally, we show that the contribution of GRBs to the flux of the extragalactic gamma-ray background measured with EGRET is at least 0.01% and likely around 0.1%.

As suggested by some extensions of the Standard Model of particle physics, dark matter may be a super-weakly interacting lightest stable particle, while the next-to-lightest particle (NLP) is charged and meta-stable. One could test such a possibility with neutrino telescopes, by detecting the charged NLPs produced in high-energy neutrino collisions with Earth matter. We study the production of charged NLPs by both atmospheric and astrophysical neutrinos; only the latter, which is largely uncertain and has not been detected yet, was the focus of previous studies. We compute the resulting fluxes of the charged NLPs, compare those of different origins, and analyze the dependence on the underlying particle physics setup. We point out that even if the astrophysical neutrino flux is very small, atmospheric neutrinos, especially those from the prompt decay of charmed mesons, may provide a detectable flux of NLP pairs at neutrino telescopes such as IceCube. We also comment on the flux of charged NLPs expected from proton-nucleon collisions, and show that, for theoretically motivated and phenomenologically viable models, it is typically sub-dominant and below detectable rates.

Highly relativistic jets are a key element of current gamma-ray burst models, where the jet kinetic energy is converted to radiation energy at optically thin shocks. High-energy neutrinos are also expected, from interactions of protons accelerated in the same shocks. Here we revisit the early evolution of a relativistic jet, while the jet is still inside the star, and investigate its neutrino emission. In particular we study propagation of mildly relativistic and ultrarelativistic jets through a type Ib progenitor, and follow reverse shocks as the jets cross the star. We show that protons can be accelerated to 10^4-10^5 GeV at reverse shocks, and efficiently produce mesons. The mesons experience significant cooling, suppressing subsequent neutrino emission. We show, however, that the neutrino yield from the reverse shock is still reasonably large, especially for low-luminosity and long-duration jets, where meson cooling is less severe. We discuss implications of our results in the context of neutrinos from choked jets, which are completely shock heated and do not break out of the star. From a choked jet with isotropic equivalent energy of 10^{53} erg at 10 Mpc, we expect ~20 neutrino events at IceCube.

The proton-proton fusion reaction, $pp\to de^+\nu$, is studied in pionless effective field theory (EFT) with di-baryon fields up to next-to leading order. With the aid of the di-baryon fields, the effective range corrections are naturally resummed up to the infinite order and thus the calculation is greatly simplified. Furthermore, the low-energy constant which appears in the axial-current-di-baryon-di-baryon contact vertex is fixed through the ratio of two- and one-body matrix elements which reproduces the tritium lifetime very precisely. As a result we can perform a parameter free calculation for the process. We compare our numerical result with those from the accurate potential model and previous pionless EFT calculations, and find a good agreement within the accuracy better than 1%.

Cosmic rays produced in cluster accretion and merger shocks provide pressure to the intracluster medium (ICM) and affect the mass estimates of galaxy clusters. Although direct evidence for cosmic-ray ions in the ICM is still lacking, they produce gamma-ray emission through the decay of neutral pions produced in their collisions with ICM nucleons. We investigate the capability of the Gamma-ray Large Area Space Telescope (GLAST) and imaging atmospheric Cerenkov telescopes (IACTs) for constraining the cosmic-ray pressure contribution to the ICM. We show that GLAST can be used to place stringent upper limits, a few per cent for individual nearby rich clusters, on the ratio of pressures of the cosmic rays and thermal gas. We further show that it is possible to place tight (<~10%) constraints for distant (z <~ 0.25) clusters in the case of hard spectrum, by stacking signals from samples of known clusters. The GLAST limits could be made more precise with the constraint on the cosmic-ray spectrum potentially provided by IACTs. Future gamma-ray observations of clusters can constrain the evolution of cosmic-ray energy density, which would have important implications for cosmological tests with upcoming X-ray and Sunyaev-Zel'dovich effect cluster surveys.

In the indirect detection of dark matter through its annihilation products, the signals depend on the square of the dark matter density, making precise knowledge of the distribution of dark matter in the Universe critical for robust predictions. Many studies have focused on regions where the dark matter density is greatest, e.g., the Galactic Center, as well as on the cosmic signal arising from all halos in the Universe. We focus on the signal arising from the whole Milky Way halo; this is less sensitive to uncertainties in the dark matter distribution, and especially for flatter profiles, this halo signal is larger than the cosmic signal. We illustrate this by considering a dark matter model in which the principal annihilation products are neutrinos. Since neutrinos are the least detectable Standard Model particles, a limit on their flux conservatively bounds the dark matter total self-annihilation cross section from above. By using the Milky Way halo signal, we show that previous constraints using the cosmic signal can be improved on by 1-2 orders of magnitude; dedicated experimental analyses should be able to improve both by an additional 1-2 orders of magnitude.

There has been growing interest in the possibility of testing more precisely the assumption of statistical isotropy of primordial density perturbations. If it is to be tested with galaxy surveys at distance scales <~ 10 Mpc, then nonlinear evolution of anisotropic power must be understood. To this end, we calculate the angular dependence of the power spectrum to third order in perturbation theory for a primordial power spectrum with a quadrupole dependence on the wavevector direction. Our results suggest that primordial power anisotropies will be suppressed by <~ 7% in the quasilinear regime. We also show that the skewness in the statistically anisotropic theory differs by no more than 1% from that in the isotropic theory.

The origin of the cosmic gamma-ray background (CGB) is a longstanding mystery in high-energy astrophysics. Possible candidates include ordinary astrophysical objects such as unresolved blazars, as well as more exotic processes such as dark matter annihilation. While it would be difficult to distinguish them from the mean intensity data alone, one can use anisotropy data instead. We investigate the CGB anisotropy both from unresolved blazars and dark matter annihilation (including contributions from dark matter substructures), and find that the angular power spectra from these sources are very different. We then focus on detectability of dark matter annihilation signals using the anisotropy data, by treating the unresolved blazar component as a known background. We find that the dark matter signature should be detectable in the angular power spectrum of the CGB from two-year all-sky observations with the Gamma Ray Large Area Space Telescope (GLAST), as long as the dark matter annihilation contributes to a reasonable fraction, e.g., >~ 0.3, of the CGB at around 10 GeV. We conclude that the anisotropy measurement of the CGB with GLAST should be a powerful tool for revealing the CGB origin, and potentially for the first detection of dark matter annihilation.

Blazars, a beamed population of active galactic nuclei, radiate high-energy gamma-rays, and thus are a good target for the Gamma Ray Large Area Space Telescope (GLAST). As the blazars trace the large-scale structure of the universe, one may observe spatial clustering of blazars. We calculate the angular power spectrum of blazars that would be detected by GLAST. We show that we have the best chance of detecting their clustering at large angular scales, \theta >~ 10 deg, where shot noise is less important, and the dominant contribution to the correlation comes from relatively low redshift, z <~ 0.1. The GLAST can detect the correlation signal, if the blazars detected by GLAST trace the distribution of low-z quasars observed by optical galaxy surveys, which have the bias of unity. If the bias of blazars is greater than 1.5, GLAST will detect the correlation signal unambiguously. We also find that GLAST may detect spatial clustering of clusters of galaxies in gamma-rays. The shape of the angular power spectrum is different for blazars and clusters of galaxies; thus, we can separate these two contributions on the basis of the shape of the power spectrum.

We investigate contributions to the extragalactic gamma-ray background (EGB) due to neutralino dark matter (DM) pair-annihilation into photons, from DM density enhancements (minispikes) surrounding intermediate-mass black holes (IMBHs). We focus on two IMBH formation scenarios; our conservative scenario where IMBHs are remnants of Population-III stars, and our optimistic scenario here IMBHs are formed in protogalactic disks. In both scenarios, their formation in pregalactic halos at high redshift lead to the formation of minispikes that are bright sources of gamma-ray photons. Taking into account minispike depletion processes, we only sum contributions from a cosmological distribution of IMBHs with maintained minispikes. Our conservative scenario (BH mass 10^2 M_sun with a r^{-3/2} minispike) predicts gamma-ray fluxes that are an order larger than the equivalent flux, using the same DM parameters (mass 100 GeV and annihilation cross-section 3 \times 10^{-26} cm^3 s^{-1}, from the host halo without IMBH minispikes. Our optimistic scenario (BH mass 10^5 M_sun with a r^{-7/3} minispike) predicts fluxes that are three orders larger, that can reach current EGB observations taken by EGRET (DM parameters as above). This fact may serve interesting consequences for constraining DM parameters and elucidating the true nature of IMBHs. Additionally, we determine the spectra of DM annihilation into monochromatic gamma-rays, and show that its flux can be within observational range of GLAST, providing a potential `smoking-gun' signature of DM.

The total cross section for radiative neutron capture on a proton, $np \to d \gamma$, is evaluated at big bang nucleosynthesis (BBN) energies. The electromagnetic transition amplitudes are calculated up to next-to leading order within the framework of pionless effective field theory with dibaryon fields. We also calculate the $d\gamma\to np$ cross section and the photon analyzing power for the $d\vec{\gamma}\to np$ process from the amplitudes. The values of low energy constants that appear in the amplitudes are estimated by a Markov Chain Monte Carlo analysis using the relevant low energy experimental data. Our result agrees well with those of other theoretical calculations except for the $np\to d\gamma$ cross section at some energies estimated by an R-matrix analysis. We also study the uncertainties in our estimation of the $np\to d\gamma$ cross section at relevant BBN energies and find that the estimated cross section is reliable to within $\sim$1% error.

The time-integrated luminosity and average energy of the neutrino emission spectrum are essential diagnostics of core-collapse supernovae. The SN 1987A electron antineutrino observations by the Kamiokande-II and IMB detectors are only roughly consistent with each other and theory. Using new measurements of the star formation rate history, we reinterpret the Super-Kamiokande upper bound on the electron antineutrino flux from all past supernovae as an excluded region in neutrino emission parameter space. A gadolinium-enhanced Super-Kamiokande should be able to jointly measure these parameters, and a future megaton-scale detector would enable precision studies.

High-energy photons from dark matter annihilation contribute to the cosmic gamma-ray background (CGB). Since dark matter particles are weakly interacting, annihilation can happen only in high density regions such as dark matter halos. The precise shape of the energy spectrum of CGB depends on the nature of dark matter particles, as well as the cosmological evolution of dark matter halos. In order to discriminate between the signals from dark matter annihilation and other astrophysical sources, however, the information from the energy spectrum may not be sufficient. We show that dark matter annihilation also produces a characteristic anisotropy of the CGB, which provides a powerful tool for testing the origin. We develop the formalism based on a halo model approach to calculate the three-dimensional power spectrum of dark matter clumping, which determines the power spectrum of annihilation signals. We show that the sensitivity of future gamma-ray detectors such as GLAST should allow us to measure the angular power spectrum of CGB anisotropy, if dark matter particles are supersymmetric neutralinos and they account for most of the observed mean intensity of CGB in GeV region. On the other hand, if dark matter has a relatively small mass, and accounts for most of the CGB in MeV region, then the future Advanced Compton Telescope should be able to measure the anisotropy in MeV region. As the intensity of photons from annihilation is proportional to the density squared, we show that the predicted shape of the angular power spectrum of gamma rays from dark matter annihilation is different from that due to other astrophysical sources such as blazars. Therefore, the angular power spectrum of the CGB provides a smoking-gun signature of dark matter annihilation.

We review our recent calculations of electroweak processes involving the deuteron, based on pionless effective field theory with dibaryon fields. These calculations are concerned with neutron-neutron fusion and np -> d gamma at BBN energies.

While existing detectors would see a burst of many neutrinos from a Milky Way supernova, the supernova rate is only a few per century. As an alternative, we propose the detection of ~ 1 neutrino per supernova from galaxies within 10 Mpc, in which there were at least 9 core-collapse supernovae since 2002. With a future 1-Mton scale detector, this could be a faster method for measuring the supernova neutrino spectrum, which is essential for calibrating numerical models and predicting the redshifted diffuse spectrum from distant supernovae. It would also allow a > 10^4 times more precise trigger time than optical data alone for high-energy neutrinos and gravitational waves.

Gamma-ray bursts (GRBs) are rare but powerful explosions displaying highly relativistic jets. It has been suggested that a significant fraction of the much more frequent core-collapse supernovae are accompanied by comparably energetic but mildly relativistic jets, which would indicate an underlying supernova--GRB connection. We calculate the neutrino spectra from the decays of pions and kaons produced in jets in supernovae, and show that the kaon contribution is dominant and provides a sharp break near 20 TeV, which is a sensitive probe of the conditions inside the jet. For a supernova at 10 Mpc, 30 events above 100 GeV are expected in a 10 s burst in the IceCube detector.

Annihilating dark matter (DM) has been discussed as a possible source of gamma-rays from the galactic center (GC) and contributing to the extragalactic gamma-ray background (EGB). Assuming universality of the density profile of DM halos, we show that it is quite unlikely that DM annihilation is a main constituent of EGB, without exceeding the observed gamma-ray flux from the GC. This argument becomes stronger when we include enhancement of the density profiles by supermassive black holes or baryon cooling. The presence of substructure may loosen the constraint, but only if a very large cross section as well as the rather flat profile are realized.

Present and future observations of supernova relic neutrinos (SRNs), i.e., a cosmological neutrino background from past core-collapse supernova explosions, potentially give us useful information concerning various fields of astrophysics, cosmology and particle physics. We review recent progress of theoretical and observational studies of SRNs, particularly focusing on the detectability and also on implications for cosmic star formation history and neutrino physics.

Nonradiative neutrino decay, which is not satisfactorily constrained, possibly and significantly changes the detected neutrino signal from galactic supernova explosions. We focus on the appearance of a sharp peak due to a neutronization burst in the time profile; this phenomenon would occur if the original \nu_e, produced at the neutrinosphere and becoming \nu_{2,3} at the stellar surface, decays into a lighter antineutrino state such as \bar\nu_{1,2} coupled to \bar\nu_e. This is due to the fact that the signature of the neutronization burst is common to all numerical simulations. Therefore, the appearance of a neutronization peak in the \bar\nu_e signal, if it were detected, would clearly indicate the nonstandard properties of neutrinos. Using a newly developed formulation that includes flavor conversions inside the supernova envelope and neutrino decay during propagation in a vacuum, we calculate the expected neutrino signal at the detectors; the lifetimes of three modes \tau_{12}, \tau_{13}, and \tau_{23} are taken to be free parameters. We further introduce simple quantities, which represent a peak sharpness of time profile and spectral hardness, and investigate the parameter dependence of these quantities. As the result, it is found that they are quite dependent on the relevant parameters, but it would be quite difficult to distinguish models using the signal obtained by the Super-Kamiokande detector; the future megaton-class detectors would have sufficient sensitivity. We also compare the neutrino decay model with another mechanism, i.e., resonant spin-flavor conversion, which also may give the appearance of a neutronization peak, and conclude that these two independent mechanisms give a different signal and one can be distinguished from the other.

We study regenerated high-energy emission from the gamma-ray burst (GRB) afterglows, and compare its flux with the direct component from the same afterglow. When the intrinsic emission spectrum extends to TeV region, these very high-energy photons are significantly absorbed by the cosmic infrared background (CIB) radiation field, creating electron/positron pairs; since these pairs are highly energetic, they can scatter the cosmic microwave background radiation up to GeV energies, which may change the intrinsic afterglow light curve in the GeV region. Using the theoretical modeling given in literature and the reasonable choice of relevant parameters, we calculate the expected light curve due to the regeneration mechanism. As the result, we find that the regenerated emission could only slightly change the original light curve, even if we take a rather large value for the CIB density, independently of the density profile of surrounding medium, i.e., constant or wind-like profile. This ensures us the reliable estimation of the intrinsic GRB parameters when the high-energy observation is accessible, regardless of a large amount of uncertainty concerning the CIB density as well as extragalactic magnetic field strength.

We study the properties of short gamma-ray bursts (GRBs), based on the assumption that they are all connected to the binary neutron star mergers, whose formation mechanism contains a large amount of uncertainty. In particular, the merger time distribution is focused on, which is essential for short GRB rate estimate; we simply parameterize it as P_m(t)\propto t^\alpha with lower cutoff time-scale \tau. The short GRB rate, calculated using such a parameterization with a specific model of cosmic star formation rate, is found to be quite insensitive to the value of \tau, but considerably sensitive to \alpha. With such formation rate models, the peak flux distribution observed by BATSE is fitted. In the fitting process, the local short GRB rate \rho_0 is treated as one of free parameters, giving the best fit when \rho_0 = 0.1--2 yr^{-1} Gpc^{-3}, depending on the other free parameters as well as the formation rate model represented by (\tau,\alpha). We also discuss several implications for the intrinsic properties of possible sources of short GRBs, such as jet structure, intrinsic luminosity function, and intrinsic duration distribution.

We investigate the flux and event rate of supernova relic neutrinos (SRNs) and discuss their implications for the cosmic star formation rate. As reference models, we adopt the supernova rate model based on recent observations and the supernova neutrino spectrum numerically calculated by several groups. In the detection energy range E_e>10 MeV, which will possibly be a background-free region in the near future, the SRN event rate is found to be 1-2 yr^{-1} at a water Cerenkov detector with a fiducial volume of 22.5 kton, depending on the adopted neutrino spectrum. We also simulate the expected signal with one set of the reference models by using the Monte Carlo method and then analyze these pseudodata with several free parameters, obtaining the distribution of the best-fit values for them. In particular, we use a parameterization such that R_{SN}(z)=R_{SN}^0 (1+z)^\alpha, where R_{SN}(z) is the comoving supernova rate density at redshift z and R_{SN}^0 and \alpha are free parameters, assuming that the supernova neutrino spectrum and luminosity are well understood by way of a future Galactic supernova neutrino burst or the future development of numerical supernova simulations. The obtained 1\sigma errors for these two parameters are found to be \delta\alpha /\alpha = 30% (7.8%) and \delta R_{SN}^0/ R_{SN}^0 = 28% (7.7%) for a detector with an effective volume of 22.5 kton 5 yr (440 kton 5 yr), where one of the parameters is fixed. On the other hand, if we fix neither of the values for these two parameters, the expected errors become rather large, \delta\alpha /\alpha = 37% and \delta R_{SN}^0/ R_{SN}^0 = 55%, even with an effective volume of 440 kton 5 yr.

Resonant spin-flavour (RSF) conversions of supernova neutrinos, which is induced by the interaction between the nonzero neutrino magnetic moment and supernova magnetic fields, are studied for both normal and inverted mass hierarchy. As the case for the pure matter-induced neutrino oscillation (Mikheyev--Smirnov--Wolfenstein (MSW) effect), we find that the RSF transitions are strongly dependent on the neutrino mass hierarchy as well as the value of $\theta_{13}$. Flavour conversions are solved numerically for various neutrino parameter sets, with presupernova profile calculated by Woosley and Weaver. In particular, it is very interesting that the RSF-induced $\nu_\rme\to\bar\nu_\rme$ transition occurs, if the following conditions are all satisfied: the value of $\mu_\nu B$ ($\mu_\nu$ is the neutrino magnetic moment, and $B$ is the magnetic field strength) is sufficiently strong, the neutrino mass hierarchy is inverted, and the value of $\theta_{13}$ is large enough to induce adiabatic MSW resonance. In this case, the strong peak due to original $\nu_\rme$ emitted from neutronization burst would exist in time profile of the neutrino events detected at the Super-Kamiokande detector. If this peak were observed in reality, it would provide fruitful information on the neutrino properties. On the other hand, characters of the neutrino spectra are also different between the neutrino models, but we find that there remains degeneracy among several models. Dependence on presupernova models is also discussed.

We propose that supernova relic neutrino (SRN) observation can be used to set constraints on the neutrino decay models. Because of the long distance scale from cosmological supernovae to the Earth, SRN have possibility to provide much stronger limit than the present one obtained from solar neutrino observation. Since the currently available data are only the upper limit on the flux integrated over $E_{\bar\nu_e}>19.3$ MeV, the decay models on which we can set constraints is quite restricted; they must satisfy specific conditions such that the daughter neutrinos are active species, the neutrino mass spectrum is quasi-degenerate, and the neutrino mass hierarchy is normal. Our numerical calculation clearly indicates that the neutrino decay model with $(\tau_2/m,\tau_3/m)\alt (10^{10},10^{10})$ [s/eV], where $\tau_i$ represents the lifetime of mass eigenstates $\bar\nu_i$, appears to give the SRN flux that is larger than the current upper limit. However, since the theoretical SRN prediction contains many uncertainties concerning a supernova rate in the universe or simulation of supernova explosions, we cannot conclude that there exists the excluded parameter region of the neutrino lifetime. In the near future, further reduced upper limit is actually expected, and it will provide more severe constraints on the neutrino decay models.

We derive the cross section of neutrino-nucleon scatterings in supernova magnetic fields, including weak-magnetism and recoil corrections. Since the weak interaction violates the parity, the scattering cross section asymmetrically depends on the directions of the neutrino momenta to the magnetic field; the origin of pulsar kicks may be explained by the mechanism. An asymmetric neutrino emission (a drift flux) due to neutrino-nucleon scatterings is absent at the leading level of $\mathcal O(\mu_BB/T)$, where $\mu_B$ is the nucleon magneton, $B$ is the magnetic field strength, and $T$ is the matter temperature at a neutrinosphere. This is because at this level the drift flux of the neutrinos are exactly canceled by that of the antineutrinos. Hence, the relevant asymmetry in the neutrino emission is suppressed by much smaller coefficient of $\mathcal O(\mu_BB/m)$, where $m$ is the nucleon mass; detailed form of the relevant drift flux is also derived from the scattering cross section, using a simple diffusion approximation. It appears that the asymmetric neutrino emission is too small to induce the observed pulsar kicks. However, we note the fact that the drift flux is proportional to the deviation of the neutrino distribution function from the value of thermal equilibrium at neutrinosphere. Since the deviation can be large for non-electron neutrinos, it is expected that there occurs cancellation between the deviation and the small suppression factor of $\mathcal O(\mu_BB/m)$. Using a simple parameterization, we show that the drift flux due to neutrino-nucleon scatterings may be comparable to the leading term due to beta processes with nucleons, which has been estimated to give a relevant kick velocity when the magnetic field is sufficiently strong as $10^{15}$--$10^{16}$ G.

We study the resonant spin-flavor (RSF) conversion of supernova neutrinos, which is induced by the interaction between the nonzero neutrino magnetic moment and the supernova magnetic fields, and its dependence on presupernova models. As the presupernova models, we adopt the latest ones by Woosley, Heger, and Weaver, and, further, models with both solar and zero metallicity are investigated. Since the $(1-2Y_e)$ profile of the new presupernova models, which is responsible for the RSF conversion, suddenly drops at the resonance region, the completely adiabatic RSF conversion is not realized, even if $\mu_\nu B_0= (10^{-12}\mu_B)(10^{10} G)$, where $B_0$ is the strength of the magnetic field at the surface of the iron core. In particular for the model with zero metallicity, the conversion is highly nonadiabatic in the high energy region, reflecting the $(1-2Y_e)$ profile of the model. In calculating the flavor conversion, we find that the shock wave propagation, which changes density profiles drastically, is a much more severe problem than it is for the pure Mikheyev-Smirnov-Wolfenstein (MSW) conversion case. This is because the RSF effect occurs at a far deeper region than the MSW effect. To avoid the uncertainty concerning the shock propagation, we restrict our discussion to 0.5 s after the core bounce (and for more conservative discussion, 0.25 s), during which the shock wave is not expected to affect the RSF region. We also evaluate the energy spectrum at the Super-Kamiokande detector for various models using the calculated conversion probabilities, and find that it is very difficult to obtain useful information on the supernova metallicities and magnetic fields or on the neutrino magnetic moment from the supernova neutrino observation. Future prospects are also discussed.

We investigate the flux of supernova relic neutrinos (SRN) for several neutrino oscillation models with parameters inferred from recent experimental results. In the calculation, we adopt the realistic {\it time-dependent} supernova density profile, which is very different from the static progenitor profile owing to shock propagation. The Earth matter effect is also included appropriately using realistic density profile of the Earth. As a result, these two effects are found to induce the flux difference by a few % in the detection energy range ($E_\nu > 19.3$ MeV). We also set 90% C.L. upper limit on SRN flux for various oscillation models using the recently released observational result by Super-Kamiokande (SK). In the near future, further reduced upper limit is actually expected, when the current SK data are reanalyzed using some technique to reduce background events against SRN signals. It is expected that the reduced upper limit is sufficient to provide useful implications for neutrino oscillation.

In interpreting the SNO experiments, accurate estimates of the \nu d reaction cross sections are of great importance. In our recent work, we have improved our previous calculation by updating some of its inputs and by incorporating the results of a recent effective-field-theoretical calculation. The new cross sections are slightly (\sim 1%) larger than the previously reported values. It is reasonable to assign 1% uncertainty to the \nu d cross sections reported here; this error estimate does not include radiative corrections.

We investigate resonant spin-flavor (RSF) conversions of supernova neutrinos which are induced by the interaction of neutrino magnetic moment and supernova magnetic fields. From the formulation which includes all three-flavor neutrinos and anti-neutrinos, we give a new crossing diagram that includes not only ordinary MSW resonance but also magnetically-induced RSF effect. With the diagram, it is found that four conversions occur in supernova, two are induced by the RSF effect and two by the pure MSW. We also numerically calculate neutrino conversions in supernova matter, using neutrino mixing parameters inferred from recent experimental results and a realistic supernova progenitor model. The results indicate that until 0.5 seconds after core bounce, the RSF-induced $\bar{\nu}_e \leftrightarrow \nu_\tau$ transition occurs efficiently (adiabatic resonance), when $\mu_\nu \agt 10^{-12} \mu_B (B_0 / 5 \times 10^{9} \mathrm G)^{-1}$, where $B_0$ is the strength of the magnetic field at the surface of iron core. We also evaluate the energy spectrum as a function of $\mu_\nu B_0$ at the SuperKamiokande detector and the Sudbury Neutrino Observatory using the calculated conversion probabilities, and find that the spectral deformation might have possibility to provide useful information on neutrino magnetic moment as well as magnetic field strength in supernovae.

In interpreting the SNO experiments, very accurate estimates of the \nu d reaction cross sections are of great importance. We improve the previous estimates of our group by updating some of its inputs and by taking into account the results of a recent effective-field-theoretical calculation. The new cross sections are slightly (\sim 1%) larger than the previously reported values. We present arguments that lead to the conclusion that it is reasonable to assign 1% uncertainty to the \nu d cross sections reported here.

We investigate the flux and the event rate of the supernova relic neutrino background (SRN) at the SuperKamiokande detector for various neutrino oscillation models with parameters inferred from recent experimental results. A realistic model of neutrino emission from supernova explosions and several models of the cosmic star formation history are adopted in the calculation. The number flux over entire energy range is found to be $11-15 \mathrm{cm^{-2}s^{-1}}$. We discuss the detection possibility of SRN at SuperKamiokande, comparing this SRN flux with other background neutrinos in more detail than previous studies. Even though there is no energy window in which SRN is dominant, we might detect it as the distortion of the other background event. We found in the energy range $17-25 \mathrm{MeV}$ the expected event rate at SuperKamiokande $0.4-0.8 ~\mathrm{yr^{-1}}$. In this range, ten-year observation might enable us to detect SRN signal (at one sigma level) in the case of LMA solar neutrino solution. We also investigate event rate at SNO and KamLAND. Although we can find energy window, the expected event rate is rather small (0.03 yr$^{-1}$ for SNO, 0.1 yr$^{-1}$ for KamLAND).

We estimate the flux of the GRB neutrino background and compute the event rate at SK and TITAND in the collapsar model, assuming that GRB formation rate is proportional to the star formation rate. We find that the flux and the event rate depend sensitively on the mass-accretion rate. Although the detection of signals from GRBs seems to be difficult by SK, we find that we can detect them by TITAND (optimistically about 20 events per year in the energy range where the energy of anti-electron neutrino is larger than 50 MeV). Therefore, we show that we will obtain the informations on the mass-accretion rate of collapsar by observing the GRB neutrino background, which, in turn, should give us informations on the total explosion energy of GRBs.