Similarity of jet radiation between flat spectrum radio quasars and GeV narrow-line Seyfert 1 galaxies: a universal δ-L_c correlation

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 1H0323+342) in different flux stages with one-zone leptonic models, we find a universal correlation between their Doppler factors(δ) and peak luminosities(L_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 follow the same δ-L_c correlation well. This indicates that the variations of observed luminosities may be essentially due to the Doppler boosting effect. The universal δ-L_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_c with the observed luminosity in the Fermi/LAT band(L_(LAT)), this correlation holds and it may serve as an empirical indicator of δ. We estimate the δ values with L_(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.     

The updated Bulk Lorentz Factors of Gamma-Ray Burst X-Ray Flares

X-ray flares are the most common phenomena in the afterglow phase of gamma-ray bursts(GRBs) in the Swift era,and are known as a canonical component in X-ray afterglows.In this work,we constrain the Lorentz factor of X-ray flares with an updated sample.We extensively search for X-ray light curves showing flare and jet break simultaneously.A smooth broken power law function is used to fit the jet breaks in 11 GRBs.We also use a smooth broken power law function to fit the profile of X-ray flares,and the total number of the flares is 20.We obtain the lower and upper limits of Lorentz factor(Γ_X) with the timescale,half-opening angle and mean luminosity of the X-ray flares for interstellar medium(ISM) and wind cases.The lower limits on Γ_X range from tens to a few hundred,and the upper limits are mainly about a few hundred.We also apply the limited Lorentz factor to test correlations of Γ_0-E_(γ,iso) and Γ_0-L_(γ,iso) for GRBs,and find X-ray flares in the ISM case are much more consistent with those of prompt emission than the wind case in a statistical sense for both correlations.X-ray flares are almost consistent with the trend in the correlations of Γ_0-E__(γ,iso)(L_(γ,iso)) for prompt GRBs,indicating X-ray flares and prompt bursts may have the same physical origin.     

The isotropic energy function and formation rate of short gamma-ray bursts

Gamma-ray bursts (GRBs) are brief,intense,gamma-ray flashes in the universe,lasting from a few milliseconds to a few thousand seconds.For short gamma-ray bursts (s GRBs) with duration less than 2 seconds,the isotropic energy (E_(iso)) function may be more scientifically meaningful and accurately measured than the luminosity (L_p) function.In this work we construct,for the first time,the isotropic energy function of s GRBs and estimate their formation rate.First,we derive the L_p-E_pcorrelation using22 s GRBs with known redshifts and well-measured spectra and estimate the pseduo redshifts of 334 Fermi s GRBs.Then,we adopt the Lynden-Bell c~-method to study isotropic energy functions and formation rate of s GRBs without any assumption.A strong evolution of isotropic energy E_(iso)∝(1+z)~(5.79)is found,which is comparable to that between L_pand z.After removing effect of the cosmic evolution,the isotropic energy function can be reasonably fitted by a broken power law,which isφ(E_(iso,0))∝E~(-0)_(iso),._0~(45)for dim s GRBs andφ(E_(iso,0))∝E~(-1)_(iso),._0~(11)for bright s GRBs,with the break energy 4.92×10~(49)erg.We obtain the local formation rate of s GRBs is about 17.43 events Gpc~(-3)yr~(-1).If assuming a beaming angle is 6?to 26?,the local formation rate including off-axis s GRBs is estimated asρ_(0,all)=155.79-3202.35 events Gpc~(-3)yr~(-1).     

What determines the observational differences of blazars?

We examine the scenario that the Doppler factor determines the observational differences of blazars. Significantly negative correlations are found between the observational synchrotron peak frequency and the Doppler factor. After correcting the Doppler boosting, the intrinsic peak frequency has a tight linear relation with the Doppler factor. It is interesting that this relation is consistent with the scenario that the black hole mass governs both the bulk Lorentz factor and the synchrotron peak frequency. In addition, the distinction between the kinetic jet powers of BL Lac objects and flat spectrum radio quasars disappears after the boosting factor δ~2 is considered. The negative correlation between the peak frequency and the observational isotropic luminosity, known as the blazar sequence, also disappears after the Doppler boosting is corrected. We also find that the correlation between the Compton dominance and the Doppler factor exists for all types of blazars. Therefore, this correlation is unsuitable for examining the external Compton emission dominance.     

The Relationship between the Rise Width and the Full Width of γ-ray Burst Pulses and Its Implications

We investigate the relationship between the rise width and the full width of gamma-ray burst pulses.Theoretical analysis shows that either width is proportional to Γ~(-2)Δ_(Тθ),FWHM(RRc)/c(Γ the Lorentz factor of the bulk motion,Δ_(Тθ),FWHM a local pulse's width,R_c the radius of fireballs and c the velocity of light).We study the relationship for four samples of observed pulses.We find:(1)merely the curvature effect could reproduce the relationship between the rise and full widths with the same slope as derived from the model of Qin et al.;(2)gamma-ray burst pulses,selected from both the short and long GRBs,follow the same sequence in the rise width vs. full width diagram,with the shorter pulses at one end;(3)all GRBs may intrinsically result from local Gaussian pulses.These features place constraints on the physical mechanism(s)for producing long and short GRBs.     

Comparisons of Jet Properties between GeV Radio Galaxies and Blazars

We compile a sample of spectral energy distributions(SEDs) of 12 GeV radio galaxies(RGs),including eight FR Ⅰ RGs and four FR Ⅱ RGs. These SEDs can be represented with the one-zone leptonic model. No significant unification, as expected in the unification model, is found for the derived jet parameters between FR Ⅰ RGs and BL Lacertae objects(BL Lacs) and between FR Ⅱ RGs and flat spectrum radio quasars(FSRQs). However, on average FR Ⅰ RGs have a larger γb(break Lorentz factor of electrons) and lower B(magnetic field strength) than FR Ⅱ RGs, analogous to the differences between BL Lacs and FSRQs. The derived Doppler factors(δ) of RGs are on average smaller than those of blazars, which is consistent with the unification model such that RGs are the misaligned parent populations of blazars with smaller δ. On the basis of jet parameters from SED fits, we calculate their jet powers and the powers carried by each component, and compare their jet compositions and radiation efficiencies with blazars. Most of the RG jets may be dominated by particles, like BL Lacs, not FSRQs.However, the jets of RGs with higher radiation efficiencies tend to have higher jet magnetization. A strong anticorrelation between synchrotron peak frequency and jet power is observed for GeV RGs and blazars in both the observer and co-moving frames, indicating that the "sequence" behavior among blazars, together with the GeV RGs, may be intrinsically dominated by jet power.     

Inclination of Broad Line Region in Narrow Line and Broad Line Seyfert 1 Galaxies

The sizes of the Broad Line Region (BLR) of some Seyfert 1 galax-ies and nearby quasars can be determined by the reverberation mapping method.Combining with the observed FWHM of Hβ emission line and assuming that themotion of BLR clouds is virialized, the black hole masses of these objects have beenestimated. However, this method strongly depends on the poorly-understood geom-etry and inclination of the BLR. On the other hand, a tight correlation between theblack hole mass and the bulge velocity dispersion was recently found for both activeand nearby inactive galaxies. This may provide another method, independent of theBLR geometry, for estimating the black hole mass. Using this method for estimatingthe black hole mass and combining with the measured BLR size and FWHM of Hβemission line, we derived the BLR inclination angles for 20 Seyfert I galaxies underthe assumption that the BLR is disk-like. The derived inclination angles agree wellwith those derived previously by fitting the UV continuum and Hβ emission lineprofiles. Adopting a relation between the FWHMs of [OⅢ]λ5007 forbidden line andthe stellar velocity dispersion, we also estimated the BLR inclinations for 50 nar-row line Seyfert 1 galaxies (NLSls). We found that the inclinations of broad LineSeyfert 1 galaxies (BLS1s) are systematically greater than those of NLS1s, whichseldom exceed 30. This may be an important factor that leads to the differencesbetween NLS1s and BLS1s if the BLR of NLS1s is really disk-like.     

On the Black Hole - Bulge Mass Ratios in Narrow-Line Seyfert 1 Galaxies

We present estimated ratios of the central black hole mass to the bulge mass (Mbh/Mbulge) for 15 Narrow Line Seyfert 1 galaxies (NLS1s). It is found that NLS1s apparently have lower mass ratios: the average mass ratio is about 1 × 10-4 with a spread of 2, which is one order of magnitude lower than for Broad Line AGNs and quiescent galaxies. This lower value, as compared to that established essentially for all other types of galaxies, can be accounted for by an underestimation of the black hole masses and an overestimation of the bulge masses in the NLS1s.     

GRB 170817A: a short GRB seen off-axis

The angular distribution of gamma-ray burst(GRB)jets is not yet clear.The observed luminosity of GRB 170817A is the lowest among all known short GRBs,which is best explained by the fact that our line of sight is outside of the jet opening angle,θ_(obs)θ_j,whereθ_(obs) is the angle between our line of sight and the jet axis.As inferred by gravitational wave observations,as well as radio and X-ray afterglow modeling of GRB 170817A,it is likely that θ_(obs)~20°–28°.In this work,we quantitatively consider two scenarios of angular energy distribution of GRB ejecta:a top-hat jet and a structured jet with a power law index s.For the top-hat jet model,we get a large θ_j(e.g.,θ_j10°),a rather high local (i.e., z 0.01) short GRB rate ~8–15×10~3 Gpc~(-3)yr~(-1((estimated to be 90~1850 Gpc~(-3)yr~(-1) in Fong et al.)and an extremely high(on-axis,V(~500 ke V for a typical short GRB).For the structured jet model,we use θ_(obs) to give limits on s and θ_j for typical on-axis luminosity of a short GRB(e.g.,10~(49)erg s(-1) 1051erg s(-1)),and a low on-axis luminosity case(e.g.,1049erg s(-1))gives more reasonable values of s.The structured jet model is more feasible for GRB170817A than the top-hat jet model due to the rather high local short GRB rate,and the extremely high on-axis E_(peak,0) almost rules out the top-hat jet model.GRB 170817A is likely a low on-axis luminosity GRB(1049erg s(-1))with a structured jet.     

Sequences in the Hardness Ratio-Peak Energy Plane of Gamma-Ray Bursts

The narrowness of the distribution of the peak energy of the vFv spectrum of gamma-ray bursts (GRBs) and the unification of GRB populations are great puzzles yet to be solved. We investigate the two puzzles based on the global spectral behaviors of different GRB populations, the long GRBs, the short GRBs, and the X-ray flashes (XRFs), in the HR - Ep plane (HR the spectral hardness ratio) with BATSE and HETE-2 observations. It is found that the long GRBs and the XRFs observed by HETE-2 seem to follow the same sequence in the HR-Ep plane, with the XRFs at the low end of this sequence. We fit the sequence by a universal Band function, and find that this sequence is mainly defined by the low energy index α, and is insensitive to the high energy index,β. With fixed β=-5, a best fit is given by α=-1.00 with X^2min/dof=2.2. The long and short GRBs observed by BATSE follow significantly different sequences in the HR-Ep plane, with most of the short GRBs having a larger hardness ratio than the long GRBs at a given Ep. For the long GRBs a best-fit yields α=-0.30 and β=-2.05. For the short GRBs, a best fit gives α=-0.60 with x^2min=1.1 (withβfixed at-2.0 because it is numerically unstable). The α value for the short GRBs is significantly greater than that for the long GRBs. These results indicate that the global spectral behaviors of the long GRB sample and the XRF sample are similar, while that of the short GRBs is different. The short GRBs seem to be a unique subclass of GRBs, and they are not the higher energy extension of the long GRBs.     

Morphology and structure of BzK-selected galaxies at z ～ 2 in the CANDELS-COSMOS field

Utilizing a Bz K-selection technique, we obtain 14 550 star-forming galaxies(sB z Ks) and 1763 passive galaxies(p Bz Ks) at z ～ 2 from the K-selected(KAB 22.5) catalog in the COSMOS/Ultra VISTA field. The differential number counts of sB z Ks and p Bz Ks are consistent with the results from the literature.Compared to the observed results, semi-analytic models of galaxy formation and evolution provide too few(many) galaxies at the high(low) mass end. Moreover, we find that the star formation rate and stellar mass of sB z Ks follow the relation of the main sequence. Based on HST/Wide Field Camera 3 F160 W imaging, we find a wide range of morphological diversities for sB z Ks, from diffuse to early-type spiral structures, with relatively high M20, large size and low G, while p Bz Ks have elliptical-like compact morphologies with lower M20, smaller size and higher G, indicating a more concentrated and symmetric spatial extent of stellar population distribution in p Bz Ks than sB z Ks. Furthermore, the sizes of p Bz Ks(sB z Ks) at z ～ 2 are on average two to three(one to two) times smaller than those of local early-type(late-type) galaxies with similar stellar mass. Our findings imply that the two classes have different evolution models and mass assembly histories.     

Polarization of GRB prompt emission and its application to POLAR’s data

Synchrotron emission polarization is very sensitive to the magnetic field configuration.Recently,polarization of synchrotron emission with a mixed(SM)magnetic field in the gamma-ray burst(GRB)afterglow phase has been developed.Here,we apply these SM models to the GRB prompt phase and compare their polarization properties with that of synchrotron emission in purely ordered(SO)magnetic field.We find that the polarization properties in an SM model are very similar to these in a corresponding SO model(e.g.,synchrotron emission in a mixed magnetic field with an aligned ordered part(SMA)and synchrotron emission with a purely ordered aligned magnetic field(SOA)),only with a lower polarization degree(PD).We also discuss the statistical properties of the models.We find PDs of the simulated bursts are concentrated around 25%for both SOA and synchrotron emission in a purely ordered toroidal magnetic field(SOT),while they can range from 0%to 25%for SMA and synchrotron emission in a mixed magnetic field with a toroidal ordered part(SMT),depending onξ_B value,i.e.,the ratio of magnetic reduction of the ordered magnetic field over that of random magnetic field.From statistics,if PDs of majority GRBs are non-zero,then it favors SO and SM models.Further,if there are some bright GRBs with prominently lower PDs than that of the majority GRBs,it favors SOT(SMT)models;if all the bright GRBs have comparable PDs with the majority ones,it favors SOA(SMA)models.Finally,we apply our results to POLAR’s data and find that~10%time-integrated PDs of the observed bursts favor SMA and SMT models,and theξ_B parameter of these bursts is constrained to be around 1.135.     

The Early X-ray Afterglows of Optically Bright and Dark Gamma-Ray Bursts

A systematic study on the early X-ray afterglows of both optically bright and dark gamma-ray bursts (B-GRBs and D-GRBs) observed by Swift is presented. Our sample includes 25 GRBs of which 13 are B-GRBs and 12 are D-GRBs. Our results show that the distributions of the X-ray afterglow fluxes (Fx), the gamma-ray fluxes (5r), and the ratio (Rr,x.) are similar for the two kinds of GRBs, that any observed differences should be simply statistical fluctuation. These results indicate that the progenitors of the two kinds of GRBs are of the same population with comparable total energies of explosion. The suppression of optical emission in the D-GRBs should result from circumburst but not from their central engine.     

The Role of T_(50) in the Classification of Gamma-ray Bursts

The role of T50 in classifying gamma-ray bursts (GRBs) is investigated. We take T50=0.7s as the line of division and find that some bursts belonging to the class of long bursts defined by T90≥2s now become short bursts (sample 1), while some belonging to the class of short bursts defined by T90 < 2 s now become long bursts (sample 2). We study how these sources are affected by the two methods of classification and find the change of classes of sample 1 is due to some peculiar properties of the light curves. Based on their characters, most of the bursts of sample 1 should be taken as short bursts.     

Spectral energy distributions for TeV blazars

In this work, we collected a sample of 69 TeV blazars from TeVCat, obtained their multiwavelength observations, and fitted their spectral energy distributions by using a second degree polynomial function. The structure parameters of synchrotron bumps for 68 blazars and those of inverseCompton bumps for 56 blazars were derived. Then, we conducted statistical analysis on the parameters(curvature, peak frequency, peak luminosity, bolometric luminosity and X/γ-ray spectral indexes).From our analysis and discussions, we can conclude the following:(1) There is a clear positive correlation between the synchrotron peak frequency, log ν_p~s, and the inverse-Compton peak frequency, log ν_p~(IC),and between the synchrotron peak luminosity, log ν_p~sL_(νp)~s, and the inverse-Compton peak luminosity,log ν_p~(IC)L_(νp)~(IC).(2) The correlation between the peak frequency and the curvature of synchrotron bump is clearly different from that of the inverse-Compton bump, which further indicates that there are different emission mechanisms between them.(3) There is a correlation between log νIC pand γ-ray spectral index, α_γ, for the TeV blazars: log ν_p~(IC)=-(4.59 ± 0.30)α_γ+(32.67 ± 0.59), which is consistent with previous work.(4) An "L-shape" relation is found between log ν_p~s and α_X for both TeV blazars and Fermi blazars. A significant correlation between log νs pand X-ray photon index(α_X) is found for the TeV blazars with high synchrotron peak frequency: log ν_p~s=-(3.20 ± 0.34)α_X+(24.33 ± 0.79),while the correlation is positive for low synchrotron peaked TeV sources.(5) In the αX-αγdiagram,there is also an "L-shape." The anti-correlation is consistent with results available in the literature, and we also find a positive correlation between them.(6) Inverse-Compton dominant sources have luminous bolometric luminosities.     

Highly variable γ-ray emission of CTD 135 and implications for its compact symmetric structure

The γ-ray emission properties of CTD 135,a typical compact symmetric object(CSO),are investigated with ～11-year Fermi/LAT observations.We show that it has bright and significantly variable GeV emission,with the γ-ray luminosity of L_γ～ 10~(47) erg s~(-1) and a variation index of TS_(var)=1002.A quasi-periodic oscillation(QPO) with a periodicity of ～460 days is detected in the global 95% false-alarm level.These γ-ray emission features are similar to that of blazars.Its broadband spectral energy distribution(SED) can be attributed to the radiations of the relativistic electrons accelerated in the core region and the extended region.The SED modeling shows that the γ-rays are from the core region,which has a Doppler boosting factor of δ～10.8 and relativistically moves with a small viewing angle,being similar to blazar jets.On the base of the analysis results,we propose that the episodic activity of the central engine in CTD135 results in a blazar-like jet and the bubble-like lobes as the Fermi bubbles in the Galaxy.The strong 7-ray emission with obvious variability is from the jet radiations and the symmetric radio structure is attributed to the bubbles.The jet radiation power and disk luminosity in units of Eddington luminosity of CTD 135 follow the same relation as other young radio sources,indicating that its jet radiation may also be driven by the Eddington ratio.     

Do Radio-loud Active Galactic Nuclei really follow the same MBH-σ＊ Relation as Normal Galaxies？

In an examination of the relationship between the black hole mass MBH and stellar velocity dispersionσ* in radio-loud active galactic nuclei (AGNs), we studied two effects which may cause uncertainties in the black hole mass estimates of radio-loud AGNs: the relativistic beaming effect on the observed optical continuum radiation and the orientation effect on the broad emission line width. After correcting these two effects, we re-examined the MBH-σ[OIII] relation for a sample of radio-loud and radio-quiet AGNs, and found the relation for radio-loud AGNs still deviated from that for nearby normal galaxies and radio-quiet AGNs. We also found there is no significant correlation between radio jet power and narrow [OIII] line width, indicating absence of strong interaction between radio jet and narrow line region. It may be that the deviation of the MBH-σ* relation of radio-loud AGNs is intrinsic, or that the [OIII] line width is not a good indicator ofσ* for radio-loud AGNs.     

Possible γ-ray emission of radio intermediate AGN Ⅲ Zw 2 and its implication on the evolution of jets in AGNs

AGNs with hard γ-ray emission identified so far are radio-loud.III Zw 2 is a radio intermediate AGN with a relativistic jet.We study its spectral energy distribution (SED) and find that the broad band emissions are dominated by the non-thermal emissions from the jet.We model its SED through a synchrotron + inverse Compton (IC) model.The results show that the IC component of III Zw 2 peaks at a few MeV, and the flux density drops rapidly at higher energy with photon index Γ≈3.3 above 0.1GeV.The predicted flu...     

Luminosity Function of the Cluster of Galaxies Abell 566

We investigate the Luminosity Function(LF)of the cluster of galaxies Abell 566.The photometric data of 15 intermediate-bands are obtained from the Beijing-Arizona- Taiwan-Connecticut(BATC)photometric sky survey.For each of the 15 wavebands,the LF of cluster galaxies is well modelled by the Schechter function,with characteristic luminosi- ties from-18.0 to-21.9 magnitude,from the a- to the p-band.Morphological dependence of the LF is investigated by separating the cluster members into‘red’and‘blue’subsamples. It is clear that late type galaxies have a steeper shape of LF than the early type galaxies.We also divided the sample galaxies by their local environment.It was found that galaxies in the sparser region have steeper shape of LF than galaxies in the denser region.Combining the results of morphological and environmental dependence of LFs,we show that Abell 566 is a well relaxed cluster with positive evidence of galaxy interaction and merger,and excess number of bright early type galaxies located in its denser region.     

The intrinsic γ-ray emissions of Fermi blazars

The beaming effect is important for understanding the observational properties of blazars.In this work, we collect 91 Fermi blazars with available radio Doppler factors. γ-ray Doppler factors are estimated and compared with radio Doppler factors for some sources. The intrinsic(de-beamed)γ-ray flux density(f_γ~(in)), intrinsic γ-ray luminosity(L_γ~(in)) and intrinsic synchrotron peak frequency(V_P~(in))are calculated. Then we study the correlations between finγand redshift and find that they follow the theoretical relation: log f =-2.0 log z + const. When the subclasses are considered, we find that stationary jets are perhaps dominant in low synchrotron peaked blazars. Sixty-three Fermi blazars with both available short variability time scales(?T) and Doppler factors are also collected. We find that the intrinsic relationship between L_γ~(in) and ?Tinobeys the Elliot Shapiro and Abramowicz Nobili relations. Strong positive correlation between f_γ~(in)and V_P~(in)is found, suggesting that synchrotron emissions are highly correlated with γ-ray emissions.