We continuously monitored the blazar S5 0716+714 in the optical $g$, $r$ and $i$ bands from Nov. 10, 2017 to Jun. 06, 2019. The total number of observations is 201 nights including 26973 data points. This is a very large quasi-simultaneous multicolor sample for the blazar. The average time spans and time resolutions are 3.4 hours and 2.9 minutes per night, respectively. During the period of observations, the target source in the $r$ band brightens from $14^{\rm m}.16$ to $12^{\rm m}.29$ together with five prominent sub-flares, and then first becomes fainter to $14^{\rm m}.76$ and again brightens to $12^{\rm m}.94$ with seven prominent sub-flares. For the long-term variations, we find a strong flatter when brighter (FWB) trend at a low flux state and then a weak FWB trend at a higher flux state. A weak FWB trend at a low flux state and then a strong FWB trend at a higher flux state are also reported. Most of sub-flares show the strong FWB trends, except for two flares with a weak FWB trend. The particle acceleration and cooling mechanisms together with the superposition of different FWB-slopes from sub-flares are likely to explain the optical color behaviours. A scenario of bent jet is discussed.

We perform a systematic search of quasi-periodic variabilities in $\gamma$-ray active galactic nuclei (AGNs) in the third \emph{Fermi} Large Area Telescope source catalog (3FGL). We employe two techniques, Lomb-Scargle Periodogram (LSP) and Weighted Wavelet Z-transform (WWZ), to obtain power spectra of $\gamma$-ray light curves covering from 2008 August to 2016 December. The results show that besides several objects that have been reported in previous works, an additional source, the FSRQ PKS 0601-70 has a possible quasi-periodic variability of 450 days in its $\gamma$-ray light curves with the significance of $>3\sigma$. The physical implications of our findings are discussed.

The acceleration of electrons in 3C 279 is investigated through analyzing the injected electron energy distribution (EED) in a time-dependent synchrotron self-Compton + external Compton emission model. In this model, it is assumed that relativistic electrons are continuously injected into the emission region, and the injected EED [$Q_e^\prime(\gamma^\prime)$] follows a single power-law form with low- and high-energy cutoffs $\rm \gamma_{min}'$ and $\rm \gamma_{max}'$, respectively, and the spectral index $n$, i.e, $Q_e^\prime(\gamma^\prime)\propto\gamma^{\prime-n}$. This model is applied to 14 quasi-simultaneous spectral energy distributions (SEDs) of 3C 279. The Markov Chain Monte Carlo fitting technique is performed to obtain the best-fitting parameters and the uncertainties on the parameters. The results show that the injected EED is well constrained in each state. The value of $n$ is in the range of 2.5 to 3.8, which is larger than that expected by the classic non-relativistic shock acceleration. However, the large value of $n$ can be explained by the relativistic oblique shock acceleration. The flaring activity seems to be related to an increased acceleration efficiency, reflected in an increased $\gamma'_{\rm min}$ and electron injection power.

Using a significantly enlarged \textit{Fermi}-LAT BL Lac objects (BL Lacs) sample, we construct the gamma-ray luminosity function (GLF) of BL Lacs, by the joint use of the space density distribution and source counts distribution. We use three well-studied forms of the GLF, i.e., the forms of pure density evolution (PDE), pure luminosity evolution (PLE) and luminosity-dependent density evolution (LDDE). The Markov Chain Monte Carlo (MCMC) technique is used to constrain model parameters. Our results suggest that LDDE model can give the best description for the BL Lac GLF. And the model shows that the BL Lacs with a harder GeV spectrum and a less luminosity evolve as strongly as FSRQs, and the evolution decreases as increasing luminosity. We also model the average photon spectra of BL Lacs with a double power laws model. Using this modeled spectra, BL Lacs contribute $ \sim 20\%$ of the total extragalactic gamma-ray background (EGB) at $E>100$ MeV, $\sim 100\%$ of the EGB at $E>50$ GeV, and the unresolved BL Lacs contribute $\sim 20 \%$ of the isotropic diffuse gamma-ray background (IGRB) at $E>100$ MeV. A prediction of the TeV EGB spectra are given, which may be tested by the future detectors.

Considering that the existence of relativistic particles in the protostellar jet has been confirmed by the detection of linearly polarized radio emission from the HH 80-81 jet, we search for $\gamma$-rays (which could be produced by the relativistic particles) from the HH 80-81 system using ten-year {\it Fermi}-LAT observations. A significant $\gamma$-ray excess is found in the direction of the HH 80-81 system with Test-Statistic (TS) value$>$100 , which is likely produced in the HH 80-81 jet. The $\gamma$-ray spectrum extends only to 1 GeV with a photon index of 3.5. No significant variability is found in the $\gamma$-ray emission. Combining the $\gamma$-ray and previous radio observations, we may exclude the relativistic Bremsstrahlung process for the production of the $\gamma$-rays under the assumption of a magnetic field strength of 0.1 mG. $\gamma$-rays can be produced through inverse-Compton (IC) scattering of relativistic electrons or inelastic proton-proton ($pp$) interaction. The maximum energies of electrons and protons are constrained to be $<0.1\ $TeV and $\sim\ $20 GeV, respectively.

The statistical properties of X-ray flares from two separate locations (nucleus and HST-1) in the M87 jet are investigated to reveal the physical origin of the flares. We analyse the archival \textit{Chandra} data for M87, and identify 14 flares in the nucleus and 9 flares in HST-1. The peak intensity ($I_{\rm{P}}$) and the flaring duration time ($T_{\rm{fl}}$) for each flare are obtained. It is found that the distributions of both $I_{\rm{P}}$ and $T_{\rm{fl}}$ for the nucleus obey a power-law form with a similar index. A similar result is also obtained for HST-1, and no significant inconsistency between the nucleus and HST-1 is found for the indices. Similar to solar X-ray flares, the power-law distributions of the flare event parameters can be well explained by a self-organized criticality (SOC) system, which are triggered by magnetic reconnection. Our results suggest that the flares from nucleus and HST-1 are possibly triggered by magnetic reconnection process. The consistent indices for the distributions of $I_{\rm{P}}$ and $T_{\rm{fl}}$ in the CORE and HST-1 indicate that the dimensions of the energy dissipation of the magnetic reconnection is identical in the two regions. A strong correlation between the flares in the two regions also suggests a similar physical origin for the flares.

The attenuation produced by extragalactic background light (EBL) in $\gamma$-ray spectra of blazars has been used to constrain the Hubble constant ($H_0$) and matter density ($\Omega_{\rm m}$) of the Universe. We propose to estimate $H_0$ and $\Omega_{\rm m}$ using the well measured $>$10 GeV extragalactic $\gamma$-ray background (EGB). This suggestion is based on the facts that the $>$10 GeV EGB is totally explained by the emissions from blazars, and an EBL-absorption cutoff occurs at $\sim$50 GeV in the EGB spectrum. We fit the $>$10 GeV EGB data with modeled EGB spectrum. This results in $H_0=64.9^{+4.6}_{-4.3}\rm \ km\ s^{-1}\ Mpc^{-1}$ and $\Omega_{\rm m}=0.31^{+0.13}_{-0.14}$. Note that the uncertainties may be underestimated due to the limit of our realization for EBL model. $H_0$ and $\Omega_{\rm m}$ are degenerate in our method. Independent determination of $\Omega_{\rm m}$ by other methods would improve the constraint on $H_0$.

Since the mass loss rates are the function of the mean density of the planet and the stellar irradiation, we calculated about 450 models covering planets with different densities and stellar irradiation. Our results show that the mass loss rates are dependent on the stellar irradiation and the mean density. However, the mass loss rates predicted by the energy-limited equation are higher than that of hydrodynamic model when the XUV integrated flux is higher than $\sim$2$\times$10$^{4}$ erg/cm$^{2}$/s. The overestimation can be revised if the kinetic and thermal energy of the escaping atmosphere is included in the energy-limited equation. We found that the heating efficiencies are proportional to the product of the gravitational potential of the planet and the stellar irradiation. The mean absorption radii of stellar irradiation are 1.1-1.2 R$_{p}$ for the Jupiter-like planets while they vary in the range of 1.1-1.7 R$_{p}$ for the planets with smaller sizes. We evaluated the absorption of stellar Ly$\alpha$ by planetary atmosphere and found that the deeper Ly$\alpha$ absorptions tend to locate in the high stellar irradiation and low planetary mean density regions, and vice versa. Moreover, planets with mass loss rates higher than 10$^{11} g/s$ are likely to exhibit obvious absorptions. Finally, we suggested that the absorption levels are related to the inherent properties of the exoplanets. The planets with larger sizes (or lower mean density) show strong Ly$\alpha$ absorptions. Neptune-like and Earth-like planets tend to have weak Ly$\alpha$ absorptions because of their small sizes (or high densities).

The Large High Altitude Air Shower Observatory (LHAASO) project is a new generation multi-component instrument, to be built at 4410 meters of altitude in the Sichuan province of China, with the aim to study with unprecedented sensitivity the spec trum, the composition and the anisotropy of cosmic rays in the energy range between 10$^{12}$ and 10$^{18}$ eV, as well as to act simultaneously as a wide aperture (one stereoradiant), continuously-operated gamma ray telescope in the energy range between 10$^{11}$ and $10^{15}$ eV. The experiment will be able of continuously surveying the TeV sky for steady and transient sources from 100 GeV to 1 PeV, t hus opening for the first time the 100-1000 TeV range to the direct observations of the high energy cosmic ray sources. In addition, the different observables (electronic, muonic and Cherenkov/fluorescence components) that will be measured in LHAASO will allow to investigate origin, acceleration and propagation of the radiation through a measurement of energy spec trum, elemental composition and anisotropy with unprecedented resolution. The remarkable sensitivity of LHAASO in cosmic rays physics and gamma astronomy would play a key-role in the comprehensive general program to explore the High Energy Universe. LHAASO will allow important studies of fundamental physics (such as indirect dark matter search, Lorentz invariance violation, quantum gravity) and solar and heliospheric physics. In this document we introduce the concept of LHAASO and the main science goals, providing an overview of the project.

A relativistic electron-positron pair beam can be produced in the interaction of TeV photons from a blazar with the extragalactic background light (EBL). The relativistic $e^{\pm}$ pairs would loss energy through inverse-Compton scattering (ICS) photons of cosmic microwave background (CMB) or plasma instabilities. The dominant energy-loss process is under debate. Based on the assumption that the dominant energy-loss process is ICS, the resulted cascade GeV radiation is usually used to constrain the intergalactic magnetic field (IGMF). Here, we include the energy-loss due to plasma oblique instability in the calculation of cascade gamma-ray flux, and investigate the impact of the plasma instability on the constraint of IGMF. The up-to-date GeV data and archival TeV data of the blazar 1ES 0229+200 are used. The results indicate that even if the oblique instability cooling is dominating over ICS cooling, the cascade flux could be still used to constrain the IGMF. It is found that with the ratio between the cooling rates of the oblique instability and the ICS varying from 0.1, 1 to 10, the lower limit of the IGMF putted by the cascade flux and the gamma-ray data changes from $8\times10^{-18}\ $G, $5\times10^{-18}\ $G to $10^{-18}\ $G. If the ratio between the two cooling rates is 30, the estimate of IGMF based on the cascade flux is invalid.

The Globular Cluster 47 Tucanae was firstly detected in gamma-rays by the Large Area Telescope (LAT) onboard the \emph{Fermi} Gamma-ray Space Telescope, and the gamma-ray emission has been widely attributed to the millisecond pulsars. In this work, we analyze the Fermi-LAT pass 8 data ranging from 2008 August to 2017 May and report the detection of a modulation with a period of $18.416\pm0.008$ hours at a significance level of $\sim 4.8\sigma$. This is the first time to detect a significant modulation with a period much longer than that of millisecond pulsars in gamma-rays from Globular Clusters. The periodic modulation signal displays in the {\it Swift}-BAT data as well. The phase-folded Chandra X-ray light curve of a point source may have provided an additional clue. These facts indicate the detection of the first young pulsar binary in gamma-rays in Globular Clusters.

The mechanism for producing gamma-ray quasi-periodic oscillation (QPO) in blazar is unknown. One possibility is the geometric model, in which without the need for intrinsic quasi-periodic variation, the relativistic Doppler factor changes periodically, resulting in observed gamma-ray QPO. We propose a method to test this geometric model. We analyze the {\it Fermi}-LAT data of PG 1553+113 spanning from 2008 August until 2018 February. According to 29 four-month average spectral energy distributions (SEDs) in the energy range of 0.1-300 GeV, we split the {\it Fermi}-LAT energy range into three bands: 0.1-1 GeV, 1-10 GeV, and 10-300 GeV. The spectrum in each energy range can be successfully fitted by a power-law. The light curves and photon indices in the three energy ranges are obtained. Then, light curves in three narrow energy ranges, i.e., 0.2-0.5 GeV, 2-5 GeV and 20- 40 GeV, are constructed, and the relative variability amplitudes in the three narrow energy ranges are calculated. A discrete-correlation analysis is performed for the light curves. Our results indicate that (i) the light curves in the different energy ranges follow the same pattern showed in the light curve above 0.1 GeV; (ii) the three groups of photon indices in the energy ranges of 0.1-1 GeV, 1-10 GeV, and 10-300 GeV keep nearly constant; (iii) the ratio between relative variability amplitudes in different narrow energy ranges are equal (within their errors) to the prediction by the Doppler effect. Our results support the scenario of the relativistic boost producing the gamma-ray QPO for PG 1553+113.

The energy dissipation mechanism in blazar jet is unknown. Blazar's flares could provide insights on this problem. Here we report statistical results of XMM-Newton X-ray flares of Mrk 421. We analyze all public XMM-Newton X-ray observations for Mrk 421, and construct the light curves. Through fitting light curves, we obtain the parameters of flare-profiles, such as peak flux ($F_{\rm p}$) and flaring time duration ($T_{\rm fl}$). It is found that both the distributions of $F_{\rm p}$ and $T_{\rm fl}$ obey a power-law form, with the same index of $\alpha_{\rm F}=\alpha_{\rm T}\approx1$. The statistical properties are consistent with the predictions by a self-organized criticality (SOC) system with energy dissipation in one-dimensional space. This is similar to solar flare, but with different space dimensions of the energy dissipation domain. This suggests that X-ray flaers of Mrk 421 are possibly driven by a magnetic reconnection mechanism. Moreover, in the analysis, we find that variability on timescale of $\sim1000\ $s frequently appears. Such rapid variability indicates a magnetic field of $\geq 2.1\delta_{\rm D}^{-1/3}$ G ($\delta_{\rm D}$ is the Doppler factor) in emission region.

The production site of gamma-rays in blazar jet is an unresolved problem. We present a method to locate gamma-ray emission region in the framework of one-zone emission model. From measurements of core-shift effect, the relation between the magnetic field strengths ($B'$) in the radio cores of jet and the distances ($R$) of these radio cores from central supermassive black hole (SMBH) can be inferred. Therefore once the magnetic field strength in gamma-ray emission region ($B'_{\rm diss}$) is obtained, one can use the relation of $B'$-$R$ to derive the distance ($R_{\rm diss}$) of gamma-ray emission region from SMBH. Here we evaluate the lower limit of $B'_{\rm diss}$ by using the criteria that the optical variability timescale $t_{\rm var}$ should be longer or equal to the synchrotron radiation cooling timescale of the electrons that emit optical photons. We test the method with the observations of PSK 1510-089 and BL Lacertae, and derive $R_{\rm diss}<0.15\delta_{\rm D}^{1/3}(1+A)^{2/3}\ $pc for PSK 1510-089 with $t_{\rm var}\sim$ a few hours, and $R_{\rm diss}<0.003\delta_{\rm D}^{1/3}(1+A)^{2/3}\ $pc for BL Lacertae with $t_{\rm var}\sim$ a few minutes. Here $\delta_{\rm D} $ is the Doppler factor and $A$ is the Compton dominance (i.e., the ratio of the Compton to the synchrotron peak luminosities).

We present our photometric monitoring of a flat spectrum radio quasar (FSRQ) 3C 454.3 at Yunnan observatories from 2006 to 2011. We find that the optical color of 3C 454.3 shows obvious redder-when-brighter trend, which reaches a saturation stage when the source is brighter than 15.15 mag at V band. We perform a simulation with multiple values of disk luminosity and spectral index to reproduce the magnitude-color diagram. The results show that the contamination caused by the disk radiation alone is difficult to produce the observed color variability. The variability properties during the outburst in December 2009 are also compared with $\gamma$-ray data derived from Fermi $\gamma$-ray space telescope. The flux variation of these two bands follow a linear relation with $F_{\gamma} \propto F_R^{1.14\pm0.07}$, which provides an observational evidence for external Compton process in 3C 454.3. Meanwhile, this flux correlation indicates that electron injection is the main mechanism for variability origin. We also explore the variation of the flux ratio $F_{\gamma}/F_R$ and the detailed structures in the lightcurves, and discuss some possible origins for the detailed variability behaviors.

We model simultaneous or quasi-simultaneous multi-band spectral energy distributions (SEDs) for a sample of 25 blazars that have radio core-shift measurements, where a one-zone leptonic model and Markov chain Monte Carlo technique are adopted. In the SED fitting for 23 low-synchrotron-peaked (LSP) blazars, the seed photons from broad line region (BLR) and molecular torus are considered respectively in the external Compton process. We find that the SED fitting with the seed photons from the torus are better than those utilizing BLR photons, which suggests that the $\gamma$-ray emitting region may be located outside the BLR. Assuming the magnetic field strength in the $\gamma$-ray emitting region as constrained from the SED fitting follows the magnetic field distribution as derived from the radio core-shift measurements (i.e., $B(R)\simeq B_{\rm 1 pc}(R/\rm 1pc)^{-1}$, $R$ is the distance from the central engine and $B_{\rm 1 pc}$ is the magnetic field strength at 1 pc), we further calculate the location of the $\gamma$-ray emitting region, $R_{\gamma}$, for these blazars. We find that $R_{\gamma}\sim2\times10^{4}R_{\rm S}\simeq10 R_{\rm BLR}$ ($R_{\rm S}$ is the Schwarzschild radius and $R_{\rm BLR}$ is the BLR size), where $R_{\rm BLR}$ is estimated from the broad line luminosities using the empirical correlations obtained using the reverberation mapping methods.

We present a model-dependent method to estimate the redshifts of three TeV BL Lac objects (BL Lacs) through fltting their (quasi-) simultaneous multi-waveband spectral energy distributions (SEDs) by one-zone leptonic synchrotron self-Compton (SSC) model. Considering the impact of electron energy distributions (EEDs) on the results, we use three types of EEDs, such as the power-law EED with exponential cut-ofi (PLC), the log-parabola (PLLP) and the broken power-law (BPL) EEDs, to flt the SEDs. We also use a parameter fi to describe the uncertainties of the extragalactic background light (EBL) models, as in (Abdo et al. 2010b). We then use Markov Chain Monte Carlo (MCMC) method to explore multi-dimensional parameter space and obtain the uncertainties of the model parameters based on the observational data. We apply our method to obtain the redshifts of three TeV BL Lac objects in the marginalized 68% confldence, and flnd that the PLC EED does not flt the SEDs. For 3C66A, the redshift is 0.14 - 0.31 and 0.16 - 0.32 in the BPL and PLLP EEDs; for PKS1424+240, the redshift is 0.55 - 0.68 and 0.55 - 0.67 in the BPL and PLLP EEDs; for PG1553+113, the redshift is 0.22 - 0.48 and 0.22 - 0.39 in the BPL and PLLP EEDs. We also estimate the redshift of PKS1424+240 in the high stage to be 0.46 - 0.67 in the PLLP EED, roughly consistent with that in the low stage.

We report a nominally high-confidence $\gamma$-ray quasi-periodic modulation in the blazar PKS 0301$-$243. For this target, we analyze its \emph{Fermi}-LAT Pass 8 data covering from 2008 August to 2017 May. Two techniques, i.e., the maximum likelihood optimization and the exposure-weighted aperture photometry, are used to build the $\gamma$-ray light curves. Then both the Lomb-Scargle Periodogram and the Weighted Wavelet Z-transform are applied to the light curves to search for period signals. A quasi-periodicity with a period of $2.1\pm0.3$ yr appears at the significance level of $\sim5\sigma$, although it should be noted that this putative quasi-period variability is seen in a data set barely four times longer. We speculate that this $\gamma$-ray quasi-periodic modulation might be evidence of a binary supermassive black hole.

We search for $\gamma$-ray and optical periodic modulations in a distant flat spectrum radio quasar (FSRQ) PKS 0426-380 (the redshift $z=1.1$). Using two techniques (i.e., the maximum likelihood optimization and the exposure-weighted aperture photometry), we obtain $\gamma$-ray light curves from \emph{Fermi}-LAT Pass 8 data covering from 2008 August to 2016 December. We then analyze the light curves with the Lomb-Scargle Periodogram (LSP) and the Weighted Wavelet Z-transform (WWZ). A $\gamma$-ray quasi-periodicity with a period of 3.35 $\pm$ 0.68 years is found at the significance-level of $\simeq3.6\ \sigma$. The optical-UV flux covering from 2005 August to 2013 April provided by ASI SCIENCE DATA CENTER is also analyzed, but no significant quasi-periodicity is found. It should be pointed out that the result of the optical-UV data could be tentative because of the incomplete of the data. Further long-term multiwavelength monitoring of this FSRQ is needed to confirm its quasi-periodicity.

We examine the gamma-ray quasi-periodic variability of PKS 2155-304 with the latest publicly available \emph{Fermi}-LAT Pass 8 data which covers the years from 2008 August to 2016 October. We produce the light curves in two ways, i.e., the exposure-weighted aperture photometry and the maximum likelihood optimization. The light curves are then analyzed by using Lomb-Scargle Periodogram (LSP) and Weighted Wavelet Z-transform (WWZ); and the results reveal a significant quasi-periodicity with a period of $1.74\pm0.13$ years and a significance of $\sim$ 4.9 $\sigma$. The constraint of multifrequencies quasi-periodic variabilities on blazar emission model is discussed.

We test one-zone synchrotron self-Compton (SSC) models with high-quality multiwavelength spectral energy distribution (SED) data of Mrk 421. We use Markov chain Monte Carlo (MCMC) technique to fit twelve day-scale SEDs of Mrk 421 with one-zone SSC models. Three types of electron energy distribution (EED), a log-parabola (LP) EED, a power-law log-parabola (PLLP) EED and a broken power-law (BPL) EED, are assumed in fits. We find that the one-zone SSC model with the PLLP EED provides successful fits to all the twelve SEDs. However, the one-zone SSC model with the LP and BPL EEDs fail to provide acceptable fits to the highest energy X-ray data or GeV data in several states. We therefore conclude that the one-zone SSC model works well in explaining the SEDs of Mrk 421, and the PLLP EED is preferred over the LP and BPL EEDs for Mrk 421 during the flare in March 2010. We derive magnetic field $B'\sim0.01$\ G, Doppler factor $\delta_{\rm D}\sim$30--50, and the curvature parameter of EED $r\sim1$--$10$ in the model with the PLLP EED. The evolutions of model parameters are explored. The physical implications of our results are discussed.

By modeling the broadband spectral energy distributions (SEDs) of a typical flat spectrum radio quasar (FSRQ, 3C 279) and two GeV narrow-line Seyfert 1 galaxies (NLS1s, PMN J0948+0022 and 1H 0323+342) in different flux stages with the one-zone leptonic models, we find a universal correlation between their Doppler factors ($\delta$) and peak luminosities ($L_{\rm c}$) of external Compton scattering bumps. Compiling a combined sample of FSRQs and GeV NLS1s, it is found that both FSRQs and GeV NLS1s in different stages and in different sources well follow the same $\delta$-$L_{\rm c}$ correlation. This indicates that the variations of observed luminosities may be essentially due to the Doppler boosting effect. And the universal $\delta$-$L_{\rm c}$ relation between FSRQs and GeV NLS1s in different stages may be further evidence that the particle acceleration and radiation mechanisms for the two kinds of sources are similar. In addition, by replacing $L_{\rm c}$ with the observed luminosity in the Fermi/LAT band ($L_{\rm LAT}$), this correlation holds, and it may serve as an empirical indicator of $\delta$. We estimate the $\delta$ values with $L_{\rm LAT}$ for 484 FSRQs in the Fermi/LAT Catalog and they range from 3 to 41, with a median of 16, which are statistically consistent with the values derived by other methods.

In external Compton scenario, we investigate the formation of the very hard electron spectrum in the fast-cooling regime, using a time-dependent emission model. It is shown that a very hard electron distribution $N'_{\rm e}(\gp)\propto\gp^{-p}$ with the spectral index $p\sim1.3$ is formed below the minimum energy of injection electron when inverse Compton scattering takes place in the Klein-Nishina regime, i.e., inverse Compton scattering of relativistic electrons on broad-line region radiation in flat spectrum radio quasars. This produces a very hard gamma-ray spectrum, and can reasonably explain the very hard \emph{Fermi}-LAT spectrum of the flat spectrum radio quasar 3C 279 during the extreme gamma-ray flare in 2013 December.

We present a method using principal component analysis (PCA) to process x-ray pulses with severe shape variation where traditional optimal filter methods fail. We demonstrate that PCA is able to noise-filter and extract energy information from x-ray pulses despite their different shapes. We apply this method to a dataset from an x-ray thermal kinetic inductance detector which has severe pulse shape variation arising from position-dependent absorption.

We study the dynamic changes of electron energy distribution (EED) through systematically analysing the quasi-simultaneous spectral energy distributions (SEDs) of the flat spectrum radio quasar 3C 279 in different states. With Markov chain Monte Carlo (MCMC) technique we model fourteen SEDs of 3C 279 using a leptonic model with a three-parameter log-parabola electron energy distribution (EED). The 14 SEDs can be satisfactorily fitted with the one-zone leptonic model. The observed $\gamma$ rays in 13 states are attributed to Compton scattering of external infrared photons from a surrounding dusty torus. The curved $\gamma$-ray spectrum observed during 2-8 April 2014 is well explained by the external Compton of dust radiation. It is found that there is a clear positive correlation between the curvature parameter $b$ of the EED and the electron peak energy $\gamma'_{\rm pk}$. No clear correlation between $b$ and the synchrotron peak frequency $\nu_{\rm s}$ is found, due to the varied product of Doppler factor and fluid magnetic field from state to state. We interpret the correlation of $b-\gamma'_{\rm pk}$ in a stochastic acceleration scenario. This positive correlation is in agreement with the prediction in the stage when the balance between acceleration and radiative cooling of the electrons is nearly established in the case of the turbulence spectral index $q=2$.

Precise spectra of 3C 279 in the 0.5-70 keV range, obtained during two epochs of \emph{Swift} and \emph{NuSTAR} observations, are analyzed using a near-equipartition model. We apply a one-zone leptonic model with a three-parameter log-parabola electron energy distribution (EED) to fit the \emph{Swift} and \emph{NuSTAR} X-ray data, as well as simultaneous optical and \emph{Fermi}-LAT $\gamma$-ray data. The Markov Chain Monte Carlo (MCMC) technique is used to search the high-dimensional parameter space and evaluate the uncertainties on model parameters. We show that the two spectra can be successfully fit in near-equipartition conditions, defined by the ratio of the energy density of relativistic electrons to magnetic field $\zeta_{\rm e}$ being close to unity. In both spectra, the observed X-rays are dominated by synchrotron-self Compton photons, and the observed $\gamma$ rays are dominated by Compton scattering of external infrared photons from a surrounding dusty torus. Model parameters are well constrained. From the low state to the high state, both the curvature of the log-parabola width parameter and the synchrotron peak frequency significantly increase. The derived magnetic fields in the two states are nearly identical ($\sim1$\ G), but the Doppler factor in the high state is larger than that in the low state ($\sim$28 versus $\sim$18). We derive that the gamma-ray emission site takes place outside the broad-line region, at $\gtrsim$ 0.1 pc from the black hole, but within the dusty torus. Implications for 3C 279 as a source of high-energy cosmic-rays are discussed.

Fermi-LAT analyses show that the gamma-ray photon spectral indices Gamma_gamma of a large sample of blazars correlate with the vFv peak synchrotron frequency v_s according to the relation Gamma_gamma = d - k log v_s. The same function, with different constants d and k, also describes the relationship between Gamma_gamma and peak Compton frequency v_C. This behavior is derived analytically using an equipartition blazar model with a log-parabola description of the electron energy distribution (EED). In the Thomson regime, k = k_EC = 3b/4 for external Compton processes and k = k_SSC = 9b/16 for synchrotron self-Compton (SSC) processes, where b is the log-parabola width parameter of the EED. The BL Lac object Mrk 501 is fit with a synchrotron/SSC model given by the log-parabola EED, and is best fit away from equipartition. Corrections are made to the spectral-index diagrams for a low-energy power-law EED and departures from equipartition, as constrained by absolute jet power. Analytic expressions are compared with numerical values derived from self-Compton and external Compton scattered gamma-ray spectra from Ly alpha broad-line region and IR target photons. The Gamma_gamma vs. v_s behavior in the model depends strongly on b, with progressively and predictably weaker dependences on gamma-ray detection range, variability time, and isotropic gamma-ray luminosity. Implications for blazar unification and blazars as ultra-high energy cosmic-ray sources are presented. Arguments by Ghisellini et al. (2014) that the jet power exceeds the accretion luminosity depend on the doubtful assumption that we are viewing at the Doppler angle.

We propose a scenario for self-consistent interpretation of GeV - TeV spectrum of a distant blazar PKS 1424+240. In this scenario, ultra-high energy (UHE) protons are assumed to exist in the blazar jet and produce gamma rays through synchrotron emission emitted by relativistic protons and pair cascades (resulted from $p\gamma$ interaction); meanwhile, some of these UHE protons may escape from the jet and are injected into intergalactic space. Therefore, we assume that UHE cosmic rays (CR) originate from relativistic protons in the jet and use energy-independent escape timescale to obtain UHECR injection spectrum. Both contributions of gamma rays injected by the source and secondary gamma rays produced in interactions of UHECRs emitted by the blazar with photon background during their propagation through intergalactic space are calculated. Our results show that this scenario is able to reproduce the GeV-TeV spectrum of PKS 1424+240 self-consistently in a broad range of redshifts $0.6<z<1.3$. The required relativistic jet power $L_p\simeq3\times10^{46}\ \rm erg\ s^{-1}$ only moderately depends on the assumed source redshift, while the proton escape timescale (and the injected UHECR luminosity) strongly depends on $z$. We compare the integral TeV fluxes predicted in our scenario with the sensitivity of the planned Cherenkov Telescope Array (CTA), and discuss the implications to future observations.

We investigate the formation of TeV spectrum of a distant blazar PKS 1424+240 residing at a redshift $z\geq0.6$ in two scenarios in the frame of a lepto-hadronic jet model taking both the uncertainties of the extragalactic background light (EBL) and its redshift into account. In the first scenario, TeV emission is attributed to the synchrotron emission of pair cascades resulting from $p\gamma$ interaction; in the second scenario, TeV emission is attributed to the proton-synchrotron emission, and an internal absorption due to interaction with the photons around the jet is included. The results show that in the first scenario the 68\% upper limit of its redshift within which this scenario can explain the VERITAS TeV spectrum in 2009 well is $\sim0.75$; in the second scenario, this upper limit of the redshift becomes $\sim1.03$. However, the second scenario can be excluded because it requires an unreasonable photon field around the jet with a luminosity of $\sim10^{43}\rm \ erg\ s^{-1}$. In conclusion, the jet model can explain its TeV spectrum with a low EBL density if $0.6<z<0.75$.

We investigate the physical properties of BL Lac objects (BL Lacs) jets by modeling the quasi-simultaneous spectral energy distributions (SEDs) of 22 {\it Fermi} BL Lac objects in the frame of a simple one-zone synchrotron self-Compton (SSC) model. We obtained the best-fit model parameters and their uncertainties for each BL Lac object through the $\chi^2$-minimization procedure and discussed their implications on the physical processes. The modeling results show that the one-zone SSC model can successfully fit the SEDs of high-synchrotron-peaked BL Lacs (HBLs) and intermediate-synchrotron-peaked BL Lacs (IBLs), but fails to explain the SEDs of low-synchrotron-peaked BL Lacs (LBLs). The statistical analysis results for model parameters are summarized as follows. (1) No correlation is found between magnetic field ($B$) and the broken energy of relativistic electrons distribution ($\gamma^{\prime}_{\rm b}$) for HBLs and IBLs, but there are inverse correlations between $\gamma^{\prime}_{\rm b}$ and the radius of emitting blob ($R^{\prime}_{\rm b}$) as well as the electrons number $K^{\prime}_{\rm e}$ for HBLs and IBLs. It's therefore concluded that the variation of $\gamma^{\prime}_{\rm b}$ is mainly caused by that of $R^{\prime}_{\rm b}$ rather than $B$ for HBLs and IBLs. (2) The Poynting flux in jets can not account for the observed radiations since the power in Poynting flux is smaller than the radiative power, and the cold protons could be the primary energy carrier in the jets.

We investigate the electron energy distributions (EEDs) and the corresponding acceleration processes in the jet of PKS 0447$-$439 and estimate its redshift through modeling its observed spectral energy distribution (SED) in the frame of a one-zone synchrotron-self Compton (SSC) model. Three EEDs formed in different acceleration scenarios are assumed: the power-law with exponential cut-off (PLC) EED (shock-acceleration scenario or the case of the EED approaching equilibrium in the stochastic-acceleration scenario), the log-parabolic (LP) EED (stochastic-acceleration scenario and the acceleration dominating) and the broken power law (BPL) EED (no acceleration scenario), and then the corresponding fluxes of both synchrotron and SSC are calculated. The model is applied to PKS 0447-439 and modeling SEDs are compared to the observed SED of this object by using the Markov Chain Monte Carlo (MCMC) method. Calculating results show that PLC model fails to fit the observed SED well, while the LP and BPL models give comparably good fits for the observed SED. The results indicate that it is possible that stochastic acceleration process acts in the emitting region of PKS 0447-439 and the EED is far from equilibrium (acceleration dominating) or no acceleration process works (in the emitting region). The redshift of PKS 0447-439 is also estimated in our fitting, and $z=0.16\pm0.05$ for LP case and $z=0.17\pm0.04$ for BPL case.

We investigate the electron energy distributions (EEDs) and the acceleration processes in the jet of Mrk 421 through fitting the spectral energy distributions (SEDs) in different active states in the frame of a one-zone synchrotron self-Compton (SSC) model. After assuming two possible EEDs formed in different acceleration models: the shock accelerated power-law with exponential cut-off (PLC) EED and the stochastic turbulence accelerated log-parabolic (LP) EED, we fit the observed SEDs of Mrk 421 in both low and giant flare states by using the Markov Chain Monte Carlo (MCMC) method which constrains the model parameters in a more efficient way. Our calculating results indicate that (1) the PLC and LP models give comparably good fits for the SED in low state, but the variations of model parameters from low state to flaring can be reasonably explained only in the case of the PLC in low state; and (2) the LP model gives better fits compared to the PLC model for the SED in flare state, and the intra-day/night variability observed at GeV-TeV bands can be accommodated only in the LP model. The giant flare may be attributed to the stochastic turbulence re-acceleration of the shock accelerated electrons in low state. Therefore, we may conclude that shock acceleration is dominant in low state, while stochastic turbulence acceleration is dominant in flare state. Moreover, our result shows that the extrapolated TeV spectra from the best-fit SEDs from optical through GeV with the two EEDs are different. It should be considered in caution when such extrapolated TeV spectra are used to constrain extragalactic background light (EBL) models.

The quasi-simultaneous multi-wavelength emission of TeV blazar 3C 66A is studied by using a one-zone multi-component leptonic jet model. It is found that the quasi-simultaneous spectral energy distribution (SED) of 3C 66A can be well reproduced, especially its {\it Fermi}-LAT first 3 months average spectrum can be well reproduced by the synchrotron-self Compton (SSC) component plus external Compton (EC) component of the broad line region (BLR). Clues on its redshift and gamma-ray emission location are obtained. The results indicate the following. (i) On the redshift; The theoretical intrinsic TeV spectra can be predicted by extrapolating the reproduced GeV spectra. Through comparing this extrapolated TeV spectra with the extragalactic background light (EBL) corrected observed TeV spectra, it is suggested that the redshift of 3C 66A could be between 0.1 and 0.3, the most likely value is $\sim$ 0.2. (ii) On the gamma-ray emission location; To well reproduce the GeV emission of 3C 66A under different assumptions on BLR, the gamma-ray emission region is always required to be beyond the inner zone of BLR. The BLR absorption effect on gamma-ray emission confirms this point.

The detections of X-ray emission from the kiloparsec-scale jets of blazars and radio galaxies may imply the existence of high energy electrons in these extended jets, and these electrons could produce high energy emission through inverse Compton (IC) process. In this paper we study the non-variable hard TeV emission from a blazar. The multi-band emission consists of two components: one is the traditional synchrotron self-Compton (SSC) emission from the inner jet, and the other is the emission produced via SSC and IC scattering of cosmic microwave background (CMB) photons (IC/CMB) and extragalactic background light (EBL) photons (IC/EBL) by relativistic electrons in the extended jet under the stochastic acceleration scenario. Such a model is applied to 1ES 1101-232. The results indicate that (1) the non-variable hard TeV emission of 1ES 1101-232 can be reproduced well, which is dominated by IC/CMB emission from the extended jet, using three characteristic values of Doppler factor ($\delta_{\rm D}=5,10,15$) for the TeV emitting region in the extended jet; and (2) in the cases of $\delta_{\rm D}=15$ and 10, the physical parameters can achieve the equipartition (or quasi-equipartition) between the relativistic electrons and the magnetic field; In contrast, the physical parameters largely deviate from the equipartition for the case of $\delta_{\rm D}=5$. It is therefore concluded that the TeV emission region of 1ES 1101-232 in the extended jet should be moderately or highly beamed.

We estimate the contribution to the extragalactic gamma-ray background (EGRB) from both intrinsic and cascade emissions produced by blazars using a simple semi- analysis method for two models of the blazar gamma-ray luminosity function (GLF). For the cascade emission, we consider two possible contributions: one is due to that the flux of the cascade emission is lower than the Fermi LAT sensitivity (case I), which is independent on the extragalactic magnetic field (EGMF), another is due to the fact that the flux of the cascade emission is larger than the Fermi LAT sensitivity but the emission angle is larger than LAT point-spread-function (PSF) angle (case II), which depends on the EGMF. Our results indicate that (1) blazar contribution to the EGRB is dominant although it depends on the GLF model and the EGMF; (2) the EGMF plays an important role in estimating the contribution from the cascade emission produced by blazars, the contribution from the cascade emission will significantly alter the EGRB spectrum when the strength of the EGMF is large enough (say BEGMF > 10-12 G); and (3) since the cascade emission in case II reaches a saturation when the strength of the EGMF is ? 10-11 G, it is very possible that the contribution from the cascade emission produced by blazars can be considered as another method to probe the upper limit of the strength of the EGMF.

We study the locations of emission regions through modelling the quasi-simultaneous multi-frequency spectral energy distributions of 21 {\it Fermi} flat spectrum radio quasars (FSRQs) in the frame of a multi-component one-zone leptonic model. In our calculations, we take the detailed broad line region (BLR) structure into account and discuss the effect of the uncertainty of the BLR structure on constraining the location of the emission regions for each FSQR, meanwhile we also include both the internal and external absorptions. Our results indicate that (1) the contribution of external Compton-BLR component is important to $\gamma$-ray emission, and the energy density of external target photon fields depends on the location of the emission region, which can be derived through reproducing the observed $\gamma$-ray emission, and (2) the emission regions of FSRQs with relative low accretion disk luminosity lie in the region of $7.9\times10^{16}$\ --\ $1.3\times10^{18}$\ cm (300 -- 4300 Schwarzschild radii) from central black hole, and for FSRQs with high accretion disk luminosity, the emission regions locate in the larger region of $2.6\times10^{17}$\ --\ $4.2\times10^{18}$\ cm (300 -- 5600 Schwarzschild radii).