The first ever identification of a cosmic ray accelerator as the consequence of spacial and temporal correlation of IceCube event 170922A with flaring of a blazar TXS 0506+056 motivated us to look for other flaring blazars in Fermi-LAT 3FGL catalog, which could be correlated with IceCube high energy track events. We have studied the Fermi-LAT light curves of blazars correlated with neutrino track events. Among the eight sources identified within 2$\sigma$ angular uncertainty of the IceCube track events selected in our study, we find only one source 3FGL J2255+2409 was in flaring state during the neutrino detection. We have carried out a time dependent modelling of the multi-wavelength data from this blazar, and the neutrino event including leptonic energy losses and proton-proton interactions in its jet to determine whether it could be the origin of the neutrino event. Our lepto-hadronic model estimates a jet luminosity of $L_j = 3.6 \times10^{47}$ erg/sec during the neutrino phase of 3FGL/4FGL J2255+2409.

The multi-wavelength photon spectrum from the BL Lac object AP Librae extends from radio to TeV gamma rays. The X-ray to very high energy gamma ray emission from the extended jet of this source has been modeled earlier with inverse Compton (IC) scattering of relativistic electrons off the CMB photons. The IC/CMB model requires the kpc scale extended jet to be highly collimated with bulk Lorentz factor close to 10. Here we discuss the possibility of proton synchrotron origin of X-rays and gamma-rays from the extended jet with bulk Lorentz factor 3. This scenario requires extreme energy of protons and high magnetic field 1 mG of the extended jet with very high jet power in particles and magnetic field (which is more than 100 times the Eddington luminosity of AP Librae) to explain the very high energy gamma ray emission. Moreover, we have shown that X-ray emission from the extended jets of 3C 273 and PKS 0637-752 could be possible by proton synchrotron emission with jet powers comparable to their Eddington luminosities.

PKS 1510-089 is one of the most variable blazars. Very high energy gamma ray emission from this source was observed by H.E.S.S. during March-April 2009 and by MAGIC from February 3 to April 3, 2012 quasi-simultaneously with multi-wavelength flares. The spectral energy distributions of these flares have been modeled earlier with the external Compton mechanism which depends on our knowledge of the densities of the seed photons in the broad line region, the dusty infrared torus or a hypothetical slow sheath surrounding the jet around the radio core. Here we show that to explain the multi-wavelength data with synchrotron emission of electrons and protons the jet power should be of the order of 10^48 ergs/sec.