The VHE emission from the HBLs H 2356-309 and 1ES 1101-232 were observed by HESS telescopes during 2004--2007. Particularly the observation in 2004 from H 2356-309 and during 2004--2005 from 1ES 1101-232 were analyzed to derive strong upper limits on the EBL which was found to be consistent with the lower limits from the integrated light of resolved galaxies. Here we have used the photohadronic model corroborated by two template EBL models to fit the observed VHE gamma-ray data from these two HBLs and to predict their intrinsic spectra. We obtain very good fit to the VHE spectra of these two HBLs. However, the predicted intrinsic spectra are different for each EBL model. For the HBL H 2356-309, we obtain a flat intrinsic spectrum and for 1ES 1101-232 the spectrum is mildly harder than 2 but much softer than 1.5.

Markarian 421 is a high-peaked BL Lac object and it has undergone many strong outbursts since its discovery as a TeV source in 1992. Markarian 421 has been stud- ied intensively and was observed by various Cherenkov tele- scope arrays ever since. The outbursts of April 2004 observed by the Whipple telescope and of February 2010 by the HESS telescopes are explained well in this work by using the photohadronic model. To account for the attenuation of these high- energy gamma-rays by the extragalactic background light (EBL), we use template EBL models. The intrinsic spectrum of each epoch is different even though the high-energy protons have almost the same spectral index. We observe that this difference in intrinsic spectra is due to the change in the spectral index of the low-energy tail of the synchrotron self Compton (SSC) photons during different epochs of flaring. Our results show that the contemporaneous multiwavelength observations, particularly in the low-energy tail region of the SSC emission of the source, are important in explaining the flaring phenomenon.

Markarian 501 is a high-peaked BL Lacertae object and has undergone many major outburst since its discovery in 1996. As a part of the multiwavelength campaign, in the year 2009 this blazar was observed for 4.5 months from March 9 to August 1 and during the period April 17 to May 5 it was observed by both space and ground based observatories covering the entire electromagnetic spectrum. A very strong high energy $\gamma$-ray flare was observed on May 1 by Whipple telescope in the energy range 317 GeV to 5 TeV and the flux was about 10 times higher than the baseline flux. We use the photohadronic model complimented by the extragalactic background radiation (EBL) correction to this very high state flare and have shown that the EBL plays an important role in attenuating the very high energy flux even though Markarian 501 is in the local Universe.

Recently ANTARES collaboration presented a time dependent analysis to a selected number of flaring blazars to look for upward going muon events produced from the charge current interaction of the muon neutrinos. We use the same list of flaring blazars to look for possible positional correlation with the IceCube neutrino events. In the context of photohadronic model we propose that the neutrinos are produced within the nuclear region of the blazar where Fermi accelerated high energy protons interact with the background synchrotron/SSC photons. Although we found that some objects from the ANTARES list are within the error circles of few IceCube events, the statistical analysis shows that none of these sources have a significant correlation.

The extraordinary multi-TeV flare from 1ES 1011 +496 during February-March 2014 was observed by MAGIC telescopes for 17 nights and the average spectrum of the whole period has a non-trivial shape. We have used the photohadronic model and a template EBL model to explain the average spectrum which fits well to the flare data. The spectral index $\alpha$ is the only free parameter in our model. We have also shown that the non-trivial nature of the spectrum is due to the change in the behavior of the optical depth above $\sim 600$ GeV $\gamma$-ray energy accompanied with the high SSC flux.