Compact steep-spectrum sources (CSSs) likely represent a population of young radio-loud active galactic nuclei (AGNs) and have been identified as gamma-ray emitting sources. We present a comprehensive analysis of their gamma-ray emission observed with Fermi/LAT and establish their broadband spectral energy distributions (SEDs). We derive their jet properties by the SED fits with a two-zone leptonic model for radiations from the compact core and large-scale extended region, and explore the possible signature of a unification picture of jet radiation among subclasses of AGNs. We show that the observed gamma-rays of CSSs with significant variability are contributed by the radiation of their compact cores via the inverse Compton process of the torus photons. The derived power-law distribution index of the radiating electrons is p_1~1.5-1.8, magnetic field strength is B~0.15-0.6 G, and Doppler boosting factor is \delta~2.8-8.9. Assuming that the jet is composed of electron-positron pairs, the compact cores of CSSs are magnetized and have a high radiation efficiency, similar to that of flat spectrum radio quasars. The six CSSs on average have higher Eddington ratio and black hole mass than those non-GeV-detected CSSs, and they follow the correlation between the jet power in units of Eddington luminosity (P_jet/L_Edd) and Eddington ratio (R_Edd) with other sub-classes of AGNs, P_jet/L_Edd~R_Edd^0.52, indicating that R_Edd would be a key physical driver for the unification scheme of AGN jet radiation.

We propose a compact binary model with an eccentric orbit to explain periodically active fast radio burst (FRB) sources, where the system consists of a neutron star (NS) with strong dipolar magnetic fields and a magnetic white dwarf (WD). In our model, the WD fills its Roche lobe at periastron, and mass transfer occurs from the WD to the NS around this point. The accreted material may be fragmented into a number of parts, which arrive at the NS at different times. The fragmented magnetized material may trigger magnetic reconnection near the NS surface. The electrons can be accelerated to an ultra-relativistic speed, and therefore the curvature radiation of the electrons can account for the burst activity. In this scenario, the duty cycle of burst activity is related to the orbital period of the binary. We show that such a model may work for duty cycles roughly from ten minutes to two days. For the recently reported 16.35-day periodicity of FRB 180916.J0158+65, our model does not naturally explain such a long duty cycle, since an extremely high eccentricity ($e>0.95$) is required.

We have assembled the historical light curves of the BL Lac Object 3FGL J0449.4-4350 at optical and gamma-ray bands, the time spanning about 10 years, analyzed the periodic variability of the light curves by using four different methods (Lomb-Scargle periodogram, REDFIT38, Jurkevich and DACF). We detected a marginally possible quasi-periodic oscillation (QPO) of ~450 days. Assuming it originates from the helical motion jet in a supermassive binary black hole (SMBBH) system undergoing major merger, we estimate the primary black hole mass M~7.7*10^{9} M_sun. To explore the origin of the gamma-ray, we investigated the optical-gamma-ray correlations using discrete correlation function (DCF) method, and found that the correlation between the two bands is very significant. This strong correlation tends to imply lepton self-synchro-Compton (LSSC) model to produce the gamma-ray.

Violent multi-wavelength variabilities are observed in gamma-ray-selected blazars. We present an analysis of long-term light curves for eight bright blazars to explore the co-variation pattern in the gamma-ray and radio bands. We extract their gamma-ray light curves and spectra with data observed by the Fermi/LAT since 2008. We find diverse co-variation patterns between the gamma-ray and radio (at 43 GHz) fluxes in these sources. The gamma-ray and radio fluxes of 3C 454.3 and PKS 1633+382 are correlated without any time-lag, suggesting that they are from the same radiation region. Similar correlation is also observed in 3C 273 and PKS 1222+216, but the radio flux is lag behind the gamma-ray flux approximately ~160 days and ~290 days, respectively. This likely suggests that their gamma-ray emission regions are located at the upstream of their radio cores at 43 GHz. The gamma-ray and radio fluxes of the other four blazars are not correlated, implying that the gamma-ray and radio emission may be from different regions in their jets. The gamma-ray light curves of the eight blazars can be decomposed into some long timescale variability components and fast spike flares. We propose that they may be attributed to the central engine activity and the magnetic reconnection process or turbulence in the local emission region, respectively.

Here we confirm the short periodic variability of a subparsec supermassive binary black hole (SMBBH) candidate Mrk 231 in the extended optical photometric data set collected by the Catalina Real-Time Transient Survey (CRTS) and All-Sky Automated Survey for Supernovae (ASAS-SN). Using the Lomb-Scargle periodogram and 2DHybrid method, we detected the significant periodicity of ~ 1.1 yr beyond a damped random walk model in the CRTS+ASAS-SN optical data set. Mrk 231 has been previously proposed as a SMBBH candidate with a highly unequal mass ratio (q~ 0.03), very tight mutual separation of ~590 AU, and an orbital period of ~1.2 yr. Hence, our result further supports, even though not prove, the intriguing hypothesis that SMBBHs with low mass ratios may be more common than close-equal mass SMBBHs. This result, however, was obtained from the contribution of CRTS data with limited sampling cadence and photometric accuracy, and further monitoring of Mrk 231 is crucial to confirm the periodicity.

We study the prospects of searching for black hole (BH) binary systems with a stellar-mass BH and a non-compact visible companion, by utilizing the spectroscopic data of Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST). We simulate the Galactic BH binary population and determine its optical visibility by considering the stellar synthetic population model and the distributions of binary orbital parameters. By convolving the visibility of BH binaries with the LAMOST detection sensitivity, we predict that $\gtrsim$ 400 candidate BH binaries can be found by the low-resolution, non-time-domain survey, and $\sim$ 50-350 candidates by the LAMOST ongoing medium-resolution, time-domain spectroscopic survey. Most of the candidates are short-period (0.2-2 days) binaries with M-, K-, G-, or F-type companions, in which $\sim$ 47% have a mass function (the lower limit of the BH mass) larger than 3 $M_{\odot}$. By complementing the LAMOST spectroscopic data with other photometric/spectroscopic surveys or follow-up observations, these candidates could be confirmed. Therefore, by exploring the LAMOST data, we can enlarge the sample of dynamically confirmed BH binaries significantly, which can improve our understanding of the mass distribution of BHs and the stellar evolution model.

Most dynamically confirmed stellar-mass black holes and the candidates were originally selected from X-ray outbursts. In the present work, we search for black hole candidates in the LAMOST survey by using the spectra along with photometry from the ASAS-SN survey, where the orbital period of the binary may be revealed by the periodic light curve, such as the ellipsoidal modulation type. Our sample consists of 9 binaries, where each source contains a giant star with large radial velocity variation ($\Delta V_{\rm R} > 70~{\rm km~s^{-1}}$) and periods known from light curves. We focus on the 9 sources with long periods ($T_{\rm ph} > 5$ days) and evaluate the mass $M_2$ of the optically invisible companion. Since the observed $\Delta V_{\rm R}$ from only a few repeating spectroscopic observations is a lower limit of the real amplitude, the real mass $M_2$ can be significantly higher than the current evaluation. It is likely an efficient method to place constraints on $M_2$ by combining $\Delta V_{\rm R}$ from LAMOST and $T_{\rm ph}$ from ASAS-SN, particularly by the ongoing LAMOST Medium Resolution Survey.

We propose the thick-disc model of Gu et al. 2016 to interpret the transition between soft ultraluminous state (SUL) and supersoft ultraluminous (SSUL) state in NGC 247. As accretion rate increases, the inner disc will puff up and act as shield to block the innermost X-ray emission regions and absorb both soft and hard X-ray photons. The absorbed X-ray emission will be re-radiated as a much softer blackbody X-ray spectrum. Hence NGC 247 shows flux dips in the hard X-ray band and transits from the SUL state to the SSUL state. The $\sim 200$s transition timescale can be explained by the viscous timescale. According to our model, the inner disc in the super-soft state is thicker and has smaller viscous timescale than in the soft state. X-ray flux variability, which is assumed to be driven by accretion rate fluctuations, might be viscous time-scale invariant. Therefore, in the SSUL state, NGC 247 is more variable. The bolometric luminosity is saturated in the thick disc; the observed radius-temperature relation can therefore be naturally explained.

iPTF14hls is an unusually bright, long-lived II-P supernova (SN), whose light curve has at least five peaks. We propose that the outflows from the black hole hyperaccretion systems in the center of the collapsars should continuously inject into the envelope. For a jet-driven core-collapsar model, the outflow feedback results in prolonging the accretion timescale and fluctuating accretion rates in our analytic solutions. Thus, the long period of luminous, varying SN iPTF14hls might originate from the choked jets, which are regulated by the feedback of the strong disk outflows in a massive core-collapsar. One can expect that jet-driven iPTF14hls may last no more than approximately 3,000 days, and the luminosity may quickly decrease in the later stages. Moreover, the double-peak light curves in some SNe might be explained by the outflow feedback mechanism.

We propose a method to search for stellar-mass black hole (BH) candidates with giant companions from spectroscopic observations. Based on the stellar spectra of LAMOST Data Release 6, we obtain a sample of seven giants in binaries with large radial velocity variation $\Delta V_R > 80~{\rm km~s^{-1}}$. With the effective temperature, surface gravity, and metallicity provided by LAMOST, and the parallax given by {\it Gaia}, we can estimate the mass and radius of the giant, and therefore evaluate the possible mass of the optically invisible star in the binary. We show that the sources in our sample are potential BH candidates, and are worthy of dynamical measurement by further spectroscopic observations. Our method may be particularly valid for the selection of BH candidates in binaries with unknown orbital periods.

Outflows launched from a remnant disc of compact binary merger may have essential contribution to the kilonova emission. Numerical calculations are conducted in this work to study the structure of accretion flows and outflows. By the incorporation of limited-energy advection in the hyper-accretion discs, outflows occur naturally from accretion flows due to imbalance between the viscous heating and the sum of the advective and radiative cooling. Following this spirit, we revisit the properties of the merger outflow ejecta. Our results show that around $10^{-3} \sim 10^{-1} M_\odot$ of the disc mass can be launched as powerful outflows. The amount of unbound mass varies with the disc mass and the viscosity. The outflow-contributed peak luminosity is around $10^{40} \sim 10^{41}$ erg/s. Such a scenario can account for the observed kilonovae associated with short gamma-ray bursts, including the recent event AT2017gfo (GW170817).

We have studied the variability of \object{S5 0716+714} at radio 15 GHz and $\gamma$-ray band using three different methods. A possible periodicity of $P_{15 GHz}=266.0\pm11.5$ and $P_{\gamma}=344.0\pm16.4$ days are obtained for radio 15 GHz and $\gamma$-ray light curves, respectively. The variability may be related to the intrinsically emission mechanism. The difference between the variability timescales of radio 15 GHz and $\gamma$-ray may be due to that the emission of radio 15 GHz is produced via the synchrotron process, while the $\gamma$-ray is produced by both the SSC and EC processes.

We analyze the long-term lightcurve of 3C 454.3 observed with Fermi/LAT and investigate its relation to the flux in the radio, optical, and X-ray bands. By fitting the 1-day binned GeV lightcurve with multiple Gaussian function (MGF), we propose that the typical variability timescale in the GeV band is 1--10 days. The GeV flux variation is accompanied by the spectral variation characterized as flux-tracking, i.e., "harder when brighter". The GeV flux is correlated with the optical and X-ray fluxes, and a weak correlation between gamma-ray flux and radio flux is also observed. The gamma-ray flux is not correlated with the optical linear polarization degree for the global lightcurves, but they show a correlation for the lightcurves before MJD 56000. The power density spectrum of the global lightcurve shows an obvious turnover at ~7.7 days, which may indicate a typical variability timescale of 3C 454.3 in the gamma-ray band. This is also consistent with the derived timescales by fitting the global lightcurve with MGF. The spectral evolution and an increase of the optical linear polarization degree along with the increase of the gamma-ray flux may indicate that the radiation particles are accelerated and the magnetic field is ordered by the shock processes during the outbursts. In addition, the nature of 3C 454.3 may be consistent with the self-organized criticality system, similar to Sagittarius A, and thus the outbursts are from plasmoid ejections driven by magnetic reconnection. This may further support the idea that the jet radiation regions are magnetized.

We present a comprehensive analysis of multiple wavelength observational data of the first GeV-selected narrow-line Seyfert 1 galaxy PMN J0948+0022. We derive its lightcurves in the \gamma-ray and X-ray bands from the data observed with Fermi/LAT and Swift/XRT, and make the optical and radio lightcurves by collecting the data from the literature. These lightcurves show significant flux variations. With the LAT data we show that this source is analogue to typical flat spectrum radio quasars in the $L_\gamma-\Gamma_\gamma$ plane, where $L_\gamma$ and $\Gamma_\gamma$ are the luminosity and spectral index in the LAT energy band. The gamma-ray flux is correlated with the V-band flux with a lag of ~44 days, and a moderate quasi-periodic oscillation (QPO) with a periodicity of ~490 days observed in the LAT lightcurve. A similar QPO signature is also found in the V-band lightcurve. The gamma-ray flux is not correlated with the radio flux in 15 GHz, and no similar QPO signature is found in a confidence level of 95%. Possible mechanisms of the QPO are discussed. We propose that gravitational wave observations in the future may clarify the current plausible models for the QPO.

We collect a sample of 78 Fermi detected blazars with broad line region (BLR) data and the quasi-simultaneous multi-wave band data. By analyzing the sample, we find: (1) For whole blazar sample, there exist significant correlations between radio, $\gamma$-ray luminosity and BLR luminosity. The slope of the best-fit linear regression equation is close to the theoretical value 1. Our results provided another support of jet-disk symbiosis. (2) For FSRQ sample, we find a significant correlation between radio luminosity and optical luminosity, and weaker correlations between X-ray luminosity and radio, optical luminosity. It indicates that the X-ray band radiation mechanism is slightly different from radio and optical bands. The X-ray emission is diluted by inverse Compton scattering (IC). There is a significant correlation between BLR luminosity and $\gamma$-ray luminosity, but no correlations between X-ray luminosity and BLR, $\gamma$-ray luminosity, which suggests that the seed soft photons are different for IC process. From these results we suggest that the X-ray band emission is diluted by synchrotron self-Compton (SSC) process, $\gamma$-ray band emission is produced by external Compton (EC) process.

Bolometric luminosity is a basic physical parameter that is widely used in the study of blazars. Due to the lack of simultaneous full wavelength data, several estimators of the bolometric luminosity are being used in practice. In this paper, we study and evaluate the reliability and significance of six estimators, the $5GHz$ luminosity, the 1keV luminosity, the $\gamma$-ray luminosity, the $5GHz$ luminosity+the 1keV luminosity, the $5GHz$ luminosity+the $\gamma$-ray luminosity and the 1keV luminosity+the $\gamma$-ray luminosity, by analyzing the linear correlations between the integrated bolometric luminosity and them. Our main results are as follows. (i) All the six estimators are reliable in the sense that they are all significant correlated with the bolometric luminosity. (ii) Ranking from the higher significance of the reliability to lower one the six estimators are the the $5GHz$ luminosity+the $\gamma$-ray luminosity, the 1keV luminosity+the $\gamma$-ray luminosity, $\gamma$-ray luminosity, the $5GHz$ luminosity+the 1keV luminosity, the $5GHz$ luminosity and the 1keV luminosity. (iii) We suggest that the bolometric luminosity can be well estimated by the $\gamma$-ray luminosity using the best linear equation that given in this paper for Fermi FSRQs. (iv) According to the linear regressions obtained in the analysis, we provide calibration for each estimator.

We have monitored the quasar 3C 273 in optical $V$, $R$ and $I$ bands from 2005 to 2016. Intraday variability (IDV) is detected on seven nights. The variability amplitudes for most of nights are less than 10\% and four nights more than 20\%. When considering the nights with time spans $>4$ hours, the value of duty cycle (DC) is 14.17 per cent. Over the twelve years, the overall magnitude and color index variabilities are $\bigtriangleup I=0^{\rm m}.67$, $\bigtriangleup R=0^{\rm m}.72$, $\bigtriangleup V=0^{\rm m}.68$, and $\bigtriangleup (V-R)=0^{\rm m}.25$ respectively. The largest clear IDV has an amplitude of 42% over just 5.8 minutes and the weakest detected IDV is 5.4% over 175 minutes. The BWB (bluer when brighter) chromatic trend is dominant for 3C 273 and appears at different flux levels on intraday timescales. The BWB trend exists for short-term timescales and intermediate-term timescales but different timescales have different correlations. There is no BWB trend for our whole time-series data sets. A significant anti-correlation between BWB trend and length of timescales is found. Combining with $V$-band data from previous works, we find a possible quasi-periodicity of $P=3918\pm1112$ days. The possible explanations for the observed variability, BWB chromatic trend and periodicity are discussed.

We aim to test if a blazar candidate of uncertain-type (BCU) in the third Fermi active galactic nuclei catalog (3LAC) can be potentially classified as a BL Lac object or a flat spectrum radio quasar (FSRQ) by performing a statistical analysis of its broadband spectral properties. We find that 34% of the radio-selected BCUs (583 BCUs) are BL Lac-like and 20% of them are FSRQ-like, which maybe within 90% level of confidence. Similarly, 77.3% of the X-ray selected BCUs (176 BCUs) are evaluated as BL Lac-like and 6.8% of them may be FSRQ-like sources. And 88.7% of the BL Lac-like BCUs that have synchrotron peak frequencies available are high synchrotron peaked BL Lacs in the X-ray selected BCUs. The percentages are accordingly 62% and 7.3% in the sample of 124 optical-selected BCUs. The high ratio of source numbers of the BL Lac-like to the FSRQ-like BCUs in the X-ray and optically selected BCU samples is due to the selection effect. Examining the consistency between our evaluation and spectroscopic identification case by case with a sample of 78 radio-selected BCUs, it is found that the statistical analysis and its resulting classifications agree with the results of the optical follow-up spectroscopic observations. Our observation campaign for high-|\rho_{s}| BCUs , i.e., |\rho_{s}|>0.8, selected with our method is ongoing.

Recent simulations on super-Eddington accretion flows have shown that, apart from the diffusion process, the vertical advection based on magnetic buoyancy can be a more efficient process to release the trapped photons in the optically thick disk. As a consequence, the radiative luminosity from the accretion disk can be far beyond the Eddington value. Following this spirit, we revisit the structure and radiation of hyper-accretion disks with mass accretion rates in the range $10^{-3}\sim 10~M_{\sun}~{\rm s}^{-1}$. Our results show that, due to the strong cooling through the vertical advection, the disk temperature becomes lower than that in the classic model without the vertical advection process, and therefore the neutrino luminosity from the disk is lower. On the other hand, the gamma-ray photons released through the vertical advection can be extremely super-Eddington. We argue that the large amount of escaped gamma-ray photons may have more significant contribution to the primordial fireball than the neutrino annihilation, and may hint a link between gamma-ray bursts and kilonovae in the black hole hyper-accretion scenario.

This work focuses on a sample of seven extremely late-time X-ray flares with peak time $t_{\rm p} > 10^4 {\rm s}$, among which two flares can be confirmed as the late-time activity of central engine. The main purpose is to investigate the mechanism of such late-time flares based on the internal origin assumption. In the hyper-accreting black hole (BH) scenario, we study the possibility of two well-known mechanisms as the central engine to power such X-ray flares, i.e., the neutrino-antineutrino annihilation and the Blandford-Znajek (BZ) process. Our results show that the annihilation luminosity is far below the observational data. Thus, the annihilation mechanism cannot account for such late-time flares. For the BZ process, if the role of outflows is taken into consideration, the inflow mass rate near the horizon will be quite low such that the magnetic field will probably be too weak to power the observed X-ray flares. We therefore argue that, for the late-time flares with internal origin, the central engine is unlikely to be associated with BHs. On the contrary, a fast rotating neutron star with strong bipolar magnetic fields may be responsible for such flares.

We fit the spectral energy distributions (SEDs) of members of a large sample of Fermi 2LAC blazars to synchrotron and inverse Compton (IC) models. Our main results are as follows. (i) As suggested by previous works, the correlation between peak frequency and curvature can be explained by statistical or stochastic particle acceleration mechanisms. For BL Lacs, we find a linear correlation between synchrotron peak frequency and its curvature. The slope of the correlation is consistent with the stochastic acceleration mechanisms and confirm previous studies. For FSRQs, we also find a linear correlation, but its slope cannot be explained by previous theoretical models. (ii) We find a significant correlation between IC luminosity and synchrotron luminosity. The slope of the correlation of FSRQs is consistent with the EC process. And the slope of the correlation of BL Lac is consistent with the SSC process. (iii) We find several significant correlations between IC curvature and several basic parameters of blazars (black hole mass, broad line luminosity, the Lorentz factor of jet). We also find significant correlations between bolometric luminosity and these basic parameters of blazars which suggest that the origin of jet is a mixture of the mechanisms proposed by Blandford $\&$ Znajek and by Blandford $\&$ Payne.

We have assembled the historical variability data of Mrk 421 at radio 15 GHz, X-ray and gamma-ray bands, spanning about 6.3, 10.3 and 7.5 yr, respective- ly. We analyzed the variability by using three methods. The results indicated that there is a period of 287.6 +- 4.4 days for 15 GHz, 309.5 +- 5.8 days for X- ray and 283.4+- 4.7 days for gamma-ray, respectively. This period can be reasonably explained by the nonballistic helical motion of the emitting material. The cor- relation analysis suggested that the variabilities of radio 15 GHz, X-ray and 0-ray are remarkable correlated, and the emission of radio 15 GHz lags behind that of X-ray, and the X-ray flux lags behind the gamma-ray. This suggests that the gamma-ray derives from inverse Compton (IC) scattering of the synchrotron photons, supporting the synchrotron self-Compton (SSC) model. Moreover, the time delay between different wavebands could be explained by the shock-in-jet models, in which a moving emission region produces the radio to gamma-ray activity, implying that the emission region of gamma-ray is closer to the center than ones of X-ray and radio emission.

We have monitored the BL Lac object Mrk 501 in optical $V$, $R$ and $I$ bands from 2010 to 2015. For Mrk 501, the presence of strong host galaxy component can affect the results of photometry. After subtracting the host galaxy contributions, the source shows intraday and long-term variabilities for optical flux and color indices. The average variability amplitudes of $V$, $R$ and $I$ bands are $22.05\%, 22.25\%, 23.82\%$ respectively, and the value of duty cycle 14.87 per cent. A minimal variability timescale of 106 minutes is detected. No significant time lag between $V$ and $I$ bands is found on one night. The bluer-when-brighter (BWB) trend is dominant for Mrk 501 on intermediate, short and intraday timescales which supports the shock-in-jet model. For the long timescale, Mrk 501 in different state can have different BWB trend. The corresponding results of non-correcting host galaxy contributions are also presented.

We compiled the radio, optical, and X-ray data of blazars from the Sloan Digital Sky Survey (SDSS) database, and presented the distribution of luminosities and broad band spectral indices. The distribution of luminosities shows that the averaged luminosity of flat-spectral radio quasars (FSRQs) is larger than that of BL Lacs objects. On the other hand, the broad band spectral energy distribution reveals that FSRQs and low energy peaked BL Lac objects (LBLs) objects have similar spectral properties, but high energy peaked BL Lac objects (HBLs) have a distinct spectral property. This may be due to that different subclasses of blazars have different intrinsic environments and are at different cooling levels. Even so, a unified scheme also is revealed from the color-color diagram, which hints that there are similar physical processes operating in all objects under a range of intrinsic physical conditions or beaming parameter.

We study the correlations of the flux of the broad-line emission ($F_{BLR}$) with the X-ray emission flux, optical emission flux at 5500 \AA and radio emission flux at 5 GHz, respectively, for a large sample of 50 Blazars (39 flat-spectrum radio quasars (FSRQs) and 11 BL Lac objects). Our main results are as follows. There are very strong correlations between $F_{BLR}$ and $F_{X}$ and between $L_{BLR}$ and $L_{X}$ in both states for 39 FSRQs and the slopes of the linear regression equations are almost equal to 1. There are weak correlations between $F_{BLR}$ and $F_{X}$ and between $L_{BLR}$ and $L_{X}$ for 11 BL Lac objects in both states, and the slopes of the linear regression equations are close to 1. There are significant correlations between $F_{BLR}$ and $F_{X}$ and between $L_{BLR}$ and $L_{X}$ for 50 blazars in both states, the slopes of both the linear regression equations are also close to 1. These results support a close link between relativistic jets and accretion on to the central Kerr black hole. On the other hand, we find that BL Lac objects have low accretion efficiency $\eta$, whereas FSRQs have high accretion efficiency $\eta$. The unified model of FSRQs and BL Lac objects is also discussed.

We present a detail analysis on the spectral lags of the short gamma-ray bursts (GRBs) and compare them with that of the long GRBs by using the CGRO/BATSE GRB Catalog. Our results are as follows. (1)The spectral lag distribution of the short GRBs is significantly different from that of the long GRBs. Excluding the statistical fluctuation effect, a proportion of ~17% of the short GRBs have a negative spectral lag, i.e., the hard photons being lag behind the soft photons. We do not find any peculiar features from their light curves to distinguish these bursts from those with a positive spectral lag. We argue that a more physical mechanism dominated the hard lag may be hid behind the morphological features of the light curves. This should be a great challenge to the current GRB models. We notice that this proportion is consistent with the proportion of short GRBs correlated with nearby galaxies newly discovered by Tanvir et al., although it is unclear if these short GRBs are indeed associated with the sources originated at low redshift. (2)While the spectral lags of the long GRBs are strongly correlated with the pulse durations, they are not for the short GRBs. However, the ratios of the spectral lag to the pulse duration for the short and long GRBs normally distributed at 0.023 and 0.046, respectively, with the sample width, indicating that the curvature effect alone could not explain the difference of the spectral lags between the two types of GRBs. The hydrodynamic timescales of the outflows and the radiative processes at work in GRBs might also play an important role as suggested by Daigne and Mochkovitch.

Recently, Shen et al. (2005) studied the contributions of curvature effect of fireballs to the spectral lag and showed that the observed lags could be accounted for by the effect. Here we check their results by performing a more precise calculation with both the formulas presented in Shen et al. (2005) and Qin et al. (2004). Several other aspects which were not considered in Shen et al. (2005) are investigated. We find that in the case of ultra-relativistic motions, both formulas are identical as long as the whole fireball surface is concerned. In our analysis, the previous conclusion that the detected spectral lags could be accounted for by the curvature effect is confirmed, while the conclusion that the lag has no dependence on the radius of fireballs is not true. We find that introducing extreme physical parameters is not the only outlet to explain those observed large lags. Even for the larger lags ($\sim5s$) (see Norris et al. 2005), a wider local pulse ($\Delta t_{\theta,FWHM} = 10^{7}s$) could account for it. Some conclusions not presented in Shen et al. (2005) or those modified in our analysis are listed below: a) lag $\propto \Gamma^{-\epsilon}$ with $\epsilon>$ 2; b) lag is proportional to the local pulse width and the $FWHM$ of the observed light curves; c) a large lag requires a large $\alpha_0$ and a small $\beta_0$ as well as a large $E_{0p}$; d) when the rest frame spectrum varies with time, the lag would become larger; e) lag decreases with the increasing of $ R_c$; f) lag $\propto$ E within the certain energy range for a given Lorentz factor; g) lag is proportional to the opening angle of uniform jets when $\theta_j < 0.6\Gamma^{-1}$.

In this paper, the effect of the intrinsic distribution of cosmological candles is investigated. We find that, in the case of a narrow distribution, the deviation of the observed modulus of sources from the expected central value could be estimated within a ceratin range. We thus introduce a lower and upper limits of $\chi ^{2}$, $\chi_{\min}^{2}$ and $ \chi_{\max}^{2}$, to estimate cosmological parameters by applying the conventional minimizing $\chi ^{2}$ method. We apply this method to a gamma-ray burst (GRB) sample as well as to a combined sample including this GRB sample and an SN Ia sample. Our analysis shows that: a) in the case of assuming an intrinsic distribution of candles of the GRB sample, the effect of the distribution is obvious and should not be neglected; b) taking into account this effect would lead to a poorer constraint of the cosmological parameter ranges. The analysis suggests that in the attempt of constraining the cosmological model with current GRB samples, the results tend to be worse than what previously thought if the mentioned intrinsic distribution does exist.

Kocevski et al. (2003) found that there is a linear relation between the rise width and the full width of gamma-ray burst pulses detected by the BATSE instrument based on their empirical functions. Motivated by this, we investigate the relationship based on Qin et al. (2004) model. Theoretical analysis shows that each of the two quantities, the rise width and the full width of observed pulse, is proportional to $\Gamma^{-2}\Delta\tau_{\theta,\rm FWHM}\frac{R_c}{c}$, where $\Gamma$ is the Lorentz factor for the bulk motion, $\Delta \tau_{\theta,\rm FWHM}$ is a local pulse's width, $R_{c}$ is the radius of fireballs and c is the velocity of light. We employ the observed pulses coming from four samples to study the relationship and find that: (1) Merely the curvature effect could produce the relationship with the same slope as those derived from Qin et al. (2004) model in the rise width vs. the full width panel. (2) Gamma-ray burst pulses, long or short ones (selected from the short and long GRBs), follow the same sequence in the rise width vs. the full width panel, with the shorter pulses at the end of this sequence. (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.