The Water Cherenkov Detector Array (WCDA) is a major component of the Large High Altitude Air Shower Array Observatory (LHAASO), a new generation cosmic-ray experiment with unprecedented sensitivity, currently under construction. The WCDA is aimed at the study of TeV $\gamma$-rays. In order to evaluate the prospects of searching for TeV $\gamma$-ray sources with the WCDA, we present in this paper a projection for the one-year sensitivity of the WCDA to TeV $\gamma$-ray sources from TeVCat (footnote: http://tevcat.uchicago.edu) using an all-sky approach. Out of 128 TeVCat sources observable to the WCDA up to a zenith angle of $45^\circ$, we estimate that 42 would be detectable for one year of observations at a median energy of 1 TeV. Most of them are Galactic sources, and the extragalactic sources are Active Galactic Nuclei (AGN).

Recent HESS observations of the ~200 pc scale diffuse gamma-ray emission from the central molecular zone (CMZ) suggest the presence of a PeV cosmic-ray accelerator (PeVatron) located in the inner 10 pc region of the Galactic Center. Interestingly, the gamma-ray spectrum of the point-like source (HESS J1745-290) in the Galactic Center shows a cutoff at ~10 TeV, implying a cutoff around 100 TeV in the cosmic-ray proton spectrum. Here we propose that the gamma-ray emission from the inner and the outer regions may be explained self-consistently by run-away protons from a single, yet fading accelerator. In this model, gamma rays from the CMZ region are produced by protons injected in the past, while gamma rays from the inner region are produced by protons injected more recently. We suggest that the blast wave formed in a tidal disruption event (TDE) caused by the supermassive black hole (Sgr A*) could serve as such a fading accelerator. With typical parameters of the TDE blast wave, gamma-ray spectra of both the CMZ region and HESS J1745-290 can be reproduced simultaneously. Meanwhile, we find that the cosmic-ray energy density profile in the CMZ region may also be reproduced in the fading accelerator model when appropriate combinations of the particle injection history and the diffusion coefficient of cosmic rays are adopted.

The nearly isotropic distribution of the TeV-PeV neutrinos recently detected by IceCube suggests that they come from sources at distance beyond our Galaxy, but how far they are is largely unknown due to lack of any associations with known sources. In this paper, we propose that the cumulative TeV gamma-ray emission accompanying the production of neutrinos can be used to constrain the distance of these neutrino sources, since the opacity of TeV gamma rays due to absorption by the extragalactic background light (EBL) depends on the distance that these TeV gamma rays have travelled. As the diffuse extragalactic TeV background measured by \emph{Fermi} is much weaker than the expected cumulative flux associated with IceCube neutrinos, the majority of IceCube neutrinos, if their sources are transparent to TeV gamma rays, must come from distances larger than the horizon of TeV gamma rays. We find that above 80\% of the IceCube neutrinos should come from sources at redshift $z>0.5$. Thus, the chance for finding nearby sources correlated with IceCube neutrinos would be small. We also find that, to explain the flux of neutrinos under the TeV gamma-ray emission constraint, the redshift evolution of neutrino source density must be at least as fast as the the cosmic star-formation rate.

Star-forming galaxies, due to their high star-formation rates and hence large number of supernova remnants therein, are huge reservoirs of cosmic rays (CRs). These CRs collide with gases in the galaxies and produce high-energy neutrinos through $pp$ collisions. In this paper, we calculate the neutrino production efficiency in star-forming galaxies by considering realistic galaxy properties, such as the gas density and galactic wind in star-forming galaxies. To calculate the accumulated neutrino flux, we use the infrared luminosity function of star-forming galaxies obtained by {\em Herschel} PEP/HerMES survey recently. The intensity of CRs producing PeV neutrinos in star-forming galaxies is normalized with the observed CR flux at EeV ({1\,EeV=$10^{18}\,$eV}), assuming that supernova remnants or hypernova remnants in star-forming galaxies can accelerate protons to EeV energies. Our calculations show that the accumulated neutrino emission produced by CRs in star-forming galaxies can account for the flux and spectrum of the sub-PeV/PeV neutrinos under reasonable assumptions on the CR confinement time in these galaxies.

One attractive scenario for the excess of sub-PeV/PeV neutrinos recently reported by IceCube is that they are produced by cosmic rays in starburst galaxies colliding with the dense interstellar medium. These proton-proton ($pp$) collisions also produce high-energy gamma-rays, which finally contribute to the diffuse high-energy gamma-ray background. We calculate the diffuse gamma-ray flux with a semi-analytic approach and consider that the very high energy gamma-rays will be absorbed in the galaxies and converted into electron-position pairs, which then lose almost all their energy through synchrotron radiation in the strong magnetic fields in the starburst region. Since the synchrotron emission goes into energies below GeV, this synchrotron loss reduces the diffuse high-energy gamma-ray flux by a factor of about two, thus leaving more room for other sources to contribute to the gamma-ray background. For a $E_\nu^{-2}$ neutrino spectrum, we find that the diffuse gamma-ray flux contributes about 20% of the observed diffuse gamma-ray background in the 100 GeV range. However, for a steeper neutrino spectrum, this synchrotron loss effect is less important, since the energy fraction in absorbed gamma-rays becomes lower.