Recent observations of VER J2227+608 reveal a broken power $\gamma$-ray spectrum with the spectral index increasing from $\sim 1.8$ in the GeV energy range to $\sim 2.3$ in the TeV range. Such a spectral break can be attributed to radiative energy loss of energetic electrons in the leptonic scenario for the $\gamma$-ray emission, which, in combination with characteristic age of the nearby pulsar, can be used to constrain magnetic field in the emission region. We show that the radio and X-ray observations can also be explained in such a scenario. In the hadronic scenario, the spectral break can be attributed to diffusion of energetic ions in a turbulent medium and detailed spectral measurement can be used to constrain the diffusion coefficient. These two models, however, predict drastically different spectra above 100 TeV, which will be uncovered with future high-resolution observations, such as LHAASO and/or CTA.

The spectra fits to a sample of 34 supernova remnants (Zeng et al., 2019) are updated. $\gamma$-ray spectra of 20 supernova remnants (SNRs) with a soft TeV spectrum are further analyzed. We found that 17 of them can be fitted in the hadronic scenario with a single power-law ion distribution with an index of $\sim$ 2.6, which is significantly softer than the ion distribution inferred from $\gamma$-ray observations of star-burst galaxies. If Galactic cosmic rays are mostly produced by SNRs, this result suggests that SNRs in star-burst galaxies may never reach the phase with a soft $\gamma$-ray spectra or escape of high-energy particles from SNRs before they reach this phase with a soft $\gamma$-ray spectrum dominates the contribution of SNRs to Galactic cosmic rays.

We report the detection of GeV $\gamma$-ray emission from the very-high-energy (VHE) $\gamma$-ray source VER J2227+608 associated with the "tail" region of SNR G106.3+2.7. The GeV $\gamma$-ray emission is extended and spatially coincident with molecular clouds traced by CO emission. The broadband GeV to TeV emission of VER J2227+608 can be well fitted by a single power-law function with an index of 1.90$\pm$0.04, without obvious indication of spectral cutoff toward high energies. The pure leptonic model for the $\gamma$-ray emission can be marginally ruled out by the X-ray and TeV data. In the hadronic model, the low energy content of CRs and the hard $\gamma$-ray spectrum, in combination with the center-bright source structure, suggest that VER J2227+608 may be powered by the PWN instead of shocks of the SNR. And the cutoff energy of the proton distribution needs to be higher than $\sim$ 400 TeV, which makes it an attractive PeVatron candidate. Future observations by the upcoming Large High Altitude Air Shower Observatory (LHAASO) and the Cherenkov Telescope Array in the north (CTA-North) could distinguish these models and constrain the maximum energy of cosmic rays in supernova remnants.

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.

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$.

We study the luminosity function and formation rate of long gamma-ray bursts (GRBs) by using a maximum likelihood method. This is the first time this method is applied to a well-defined sample of GRBs that is complete in redshift. The sample is composed of 99 bursts detected by the $Swift$ satellite, 81 of them with measured redshift and luminosity for a completeness level of $82\%$. We confirm that a strong redshift evolution in luminosity (with an evolution index of $\delta=2.22^{+0.32}_{-0.31}$) or in density ($\delta=1.92^{+0.20}_{-0.21}$) is needed in order to reproduce the observations well. But since the predicted redshift and luminosity distributions in the two scenarios are very similar, it is difficult to distinguish between these two kinds of evolutions only on the basis of the current sample. Furthermore, we also consider an empirical density case in which the GRB rate density is directly described as a broken power-law function and the luminosity function is taken to be non-evolving. In this case, we find that the GRB formation rate rises like $(1+z)^{3.85^{+0.48}_{-0.45}}$ for $z\leq2$ and is proportional to $(1+z)^{-1.07^{+0.98}_{-1.12}}$ for $z\geq2$. The local GRB rate is $1.49^{+0.63}_{-0.64}$ Gpc$^{-3}$ yr$^{-1}$. The GRB rate may be consistent with the cosmic star formation rate (SFR) at $z\leq2$, but shows an enhancement compared to the SFR at $z\geq2$.

Relations between radio surface brightness ($\Sigma$) and diameter ($D$) of supernova remnants (SNRs) are important in astronomy. In this paper, following the work Duric \& Seaquist (1986) at adiabatic phase, we carefully investigate shell-type supernova remnants at radiative phase, and obtain theoretical $\Sigma$-$D$ relation at radiative phase of shell-type supernova remnants at 1 GHz. By using these theoretical $\Sigma$-$D$ relations at adiabatic phase and radiative phase, we also roughly determine phases of some supernova remnant from observation data.

Supernova remnants (SNRs) have been considered as the dominant contributors to Galactic cosmic rays. However, the relation between high-energy particles trapped in SNRs and cosmic rays observed at the Earth remains obscure. In this paper, we fit the spectral energy distributions of 35 SNRs with a simple one-zone emission model and analyze correlations of model parameters to uncover the evolution of high-energy particle distribution in SNRs. We find that (1) the particle distribution in general can be described by a broken power-law function with a high-energy cutoff for all SNRs; (2) the low-energy spectrum becomes harder and the break energy decreases with aging of SNRs, (3) for most middle-age SNRs, the energy loss timescale of electrons at the high-energy cutoff is approximately equal to the age of the corresponding remnant implying quenching of very high-energy electron acceleration; for young SNRs, this energy loss timescale is shorter than the age of SNRs implying continuous electon acceleration at the cutoff energy; and for a few old age SNRs, the energy loss timescale is longer than the corresponding age which may suggest escaping of higher energy particles from SNRs. Finally, we comment on the implications of these results on the SNR origin of Galactic cosmic rays.

Multi-wavelength observations of mature supernova remnants (SNRs), especially with recent advances in gamma-ray astronomy, make it possible to constrain energy distribution of energetic particles within these remnants. In consideration of the SNR origin of Galactic cosmic rays and physics related to particle acceleration and radiative processes, we use a simple one-zone model to fit the nonthermal emission spectra of three shell-type SNRs located within 2 degrees on the sky: RX J1713.7-3946, CTB 37B, and CTB 37A. Although radio images of these three sources all show a shell (or half-shell) structure, their radio, X-ray, and gamma-ray spectra are quite different, offering an ideal case to explore evolution of energetic particle distribution in SNRs. Our spectral fitting shows that 1) the particle distribution becomes harder with aging of these SNRs, implying a continuous acceleration process, and the particle distributions of CTB 37A and CTB 37B in the GeV range are harder than the hardest distribution that can be produced at a shock via the linear diffusive shock particle acceleration process, so spatial transport may play a role; 2) the energy loss timescale of electrons at the high-energy cutoff due to synchrotron radiation appears to be always a bit (within a factor of a few) shorter than the age of the corresponding remnant, which also requires continuous particle acceleration; 3) double power-law distributions are needed to fit the spectra of CTB 37B and CTB 37A, which may be attributed to shock interaction with molecular clouds.

The extensive catalog of $\gamma$-ray selected flat-spectrum radio quasars (FSRQs) produced by \emph{Fermi} during a four-year survey has generated considerable interest in determining their $\gamma$-ray luminosity function (GLF) and its evolution with cosmic time. In this paper, we introduce the novel idea of using this extensive database to test the differential volume expansion rate predicted by two specific models, the concordance $\Lambda$CDM and $R_{\rm h}=ct$ cosmologies. For this purpose, we use two well-studied formulations of the GLF, one based on pure luminosity evolution (PLE) and the other on a luminosity-dependent density evolution (LDDE). Using a Kolmogorov-Smirnov test on one-parameter cumulative distributions (in luminosity, redshift, photon index and source count), we confirm the results of earlier works showing that these data somewhat favour LDDE over PLE; we show that this is the case for both $\Lambda$CDM and $R_{\rm h}=ct$. Regardless of which GLF one chooses, however, we also show that model selection tools very strongly favour $R_{\rm h}=ct$ over $\Lambda$CDM. We suggest that such population studies, though featuring a strong evolution in redshift, may nonetheless be used as a valuable independent check of other model comparisons based solely on geometric considerations.

With the standard deviation for the logarithm of the re-scaled range $\langle |F(t+\tau)-F(t)|\rangle$ of simulated fractal Brownian motions $F(t)$ given in a previous paper \cite{q14}, the method of least squares is adopted to determine the slope, $S$, and intercept, $I$, of the log$(\langle |F(t+\tau)-F(t)|\rangle)$ vs $\rm{log}(\tau)$ plot to investigate the limitation of this procedure. It is found that the reduced $\chi^2$ of the fitting decreases with the increase of the Hurst index, $H$ (the expectation value of $S$), which may be attributed to the correlation among the re-scaled ranges. Similarly, it is found that the errors of the fitting parameters $S$ and $I$ are usually smaller than their corresponding standard deviations. These results show the limitation of using the simple least square method to determine the dimension of a fractal time series. Nevertheless, they may be used to reinterpret the fitting results of the least square method to determine the dimension of fractal Brownian motions more self-consistently. The currency exchange rate between Euro and Dollar is used as an example to demonstrate this procedure and a fractal dimension of 1.511 is obtained for spans greater than 30 transactions.

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 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 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).