Recent detection of the neutrino events IceCube-170922A, 13 muon-neutrino events observed in 2014-2015 and IceCube-141209A by IceCube observatory from the Blazars, namely TXS 0506+056, PKS 0502+049/TXS 0506+056 and GB6 J1040+0617 respectively in the state of enhanced gamma-ray emission, indicates the acceleration of cosmic rays in the blazar jets. The photo-meson ($p\gamma$) interaction cannot explain the IceCube observations of 13 neutrino events. The non-detection of broadline emission in the optical spectra of the IceCube blazars, however, question the hadronuclear (pp) interaction interpretation through relativistic jet meets with high density cloud. In this work, we investigate the proton blazar model in which the non-relativistic protons that come into existence under the charge neutrality condition of the blazar jet can offer sufficient target matter for $pp$ interaction with shock-accelerated protons, to describe the observed high-energy gamma-rays and neutrino signal from the said blazars. Our findings suggest that the model can explain consistently the observed electromagnetic spectrum in combination with appropriate number of neutrino events from the corresponding blazars.

Recent detection of the neutrino event, IceCube-170922A by IceCube observatory from the Blazar TXS 0506+056 in the state of enhanced gamma ray emission indicates for acceleration of cosmic rays in the blazar jet. The non-detection of the broadline emission in the optical spectrum of TXS 0506+056 and other BL Lac objects suggests that external photons emissions are weak and hence photo-meson (p-gamma) interaction may not be a favored mechanism for high energy neutrino production. The lack of broadline signatures also creates doubt about the presence of a high density cloud in the vicinity of the super-massive black hole (SMBH) of TXS 0506+056 and consequently raised question on hadronuclear (pp) interaction interpretation like relativistic jet meets with high density cloud. Here we demonstrate that non-relativistic protons in the proton blazar model, those come into existence under charge neutrality condition of the blazar jet, offer sufficient target matter for pp-interaction with shock accelerated protons and consequently the model can describe consistently the observed high energy gamma rays and neutrino signal from the blazar TXS 0506+056.

The galactic cosmic rays are generally believed to be originated in supernova remnants (SNRs), produced in diffusive shock acceleration (DSA) process in supernova blast waves driven by expanding SNRs. One of the key unsettled issue in SNR origin of cosmic ray model is the maximum attainable energy by a cosmic ray particle in the supernova shock. Recently it has been suggested that an amplification of effective magnetic field strength at the shock may take place in young SNRs due to growth of magnetic waves induced by accelerated cosmic rays and as a result the maximum energy achieved by cosmic rays in SNR may reach the knee energy instead of $\sim 200$ TeV as predicted earlier under normal magnetic field situation. In the present work we investigate the implication of such maximum energy scenarios on TeV gamma rays and neutrino fluxes from the molecular clouds interacting with the SNR W28. The authors compute the gamma-ray and neutrino flux assuming two different values for the maximum energy reached by cosmic rays in the SNR, from CR interaction in nearby molecular clouds. Both protons and nuclei are considered as accelerated particles and as target material. Our findings suggest that the issue of the maximum energy of cosmic rays in SNRs can be observationally settled by the upcoming gamma-ray experiment the Large High Altitude Air Shower Observatory (LHAASO). The estimated neutrino fluxes from the molecular clouds are , however, out of reach of the present/near future generation of neutrino telescopes.

The standard perception is that the detection of high energy (TeV energies and above) neutrinos from an astrophysical object is a conclusive evidence for the presence of hadronic cosmic rays at the source. In the present work we demonstrate that TeV neutrinos can also be originated from energetic electrons via electromagnetic interactions in different potential cosmic ray sources with flux levels comparable to that of the hadronic originated neutrinos at high energies. Our findings thus imply that at least a part of the neutrinos observed by Icecube observatory may be originated from energetic electrons. The present analysis further suggests that only a combine study of TeV gamma rays and neutrinos over a wide energy range from an astrophysical object can unambiguously identify the nature of their parents, hadrons or leptons.

It is widely believe that galactic cosmic rays are originated in supernova remnants (SNRs) where they are accelerated by diffusive shock acceleration process at supernova blast waves driven by expanding SNRs. In recent theoretical developments of the diffusive shock acceleration theory in SNRs, protons are expected to accelerate in SNRs at least up to the knee energy. If SNRs are true generator of cosmic rays, they should accelerate not only protons but also heavier nuclei with right proportion and the maximum energy of heavier nuclei should be atomic mass (Z) times that of protons. In this work we investigate the implications of acceleration of heavier nuclei in SNRs on energetic gamma rays those are produced in hadronic interaction of cosmic rays with ambient matter. Our findings suggest that the energy conversion efficiency has to be nearly double for the mixed cosmic ray composition instead of pure protons to explain the observation and secondly the gamma ray flux above few tens of TeV would be significantly higher if cosmic rays particles can attain energies Z times of the knee energy in lieu of 200 TeV, as suggested earlier for non-amplified magnetic fields. The two stated maximum energy paradigm will be discriminated in future by the upcoming gamma ray experiments like Cherenkov Telescope array (CTA).

The possibility of estimating the mass composition of primary cosmic rays above the knee of its energy spectrum through the study of high energy gamma rays, muons and neutrinos produced in the interactions of cosmic rays with the solar ambient matter and radiation has been explored. It is found that the theoretical fluxes of TeV gamma rays, muons and neutrinos from a region around $15^{o}$ of the Sun are sensitive to mass composition of cosmic rays in the PeV energy range. The experimental prospects for detection of such TeV gamma rays/neutrinos by future experiments are discussed.

In this work we examine with the help of Monte Carlo simulation whether a consistent primary energy spectrum of cosmic rays emerges from both the experimentally observed total charged particles and muon size spectra of cosmic ray extensive air showers considering primary composition may or may not change beyond the knee of the energy spectrum. It is found that EAS-TOP observations consistently infer a knee in the primary energy spectrum provided the primary is pure unchanging iron whereas no consistent primary spectrum emerges from simultaneous use of the KASCADE observed total charged particle and muon spectra. However, it is also found that when primary composition changes across the knee the estimation of spectral index of total charged particle spectrum is quite tricky, depends on the choice of selection of points near the knee in the size spectrum.