The Epoch of Reionization (EoR) depends on the complex astrophysics governing the birth and evolution of the first galaxies and structures in the intergalactic medium. EoR models rely on cosmic microwave background (CMB) observations, and in particular the large-scale E-mode polarization power spectra (EE PS), to help constrain their highly uncertain parameters. However, rather than directly forward-modelling the EE PS, most EoR models are constrained using a summary statistic -- the Thompson scattering optical depth, $\tau_e$. Compressing CMB observations to $\tau_e$ requires adopting a basis set for the EoR history. The common choice is the unphysical, redshift-symmetric hyperbolic tangent (Tanh) function, which differs in shape from physical EoR models based on hierarchical structure formation. Combining public EoR and CMB codes, 21cmFAST and CLASS, here we quantify how inference using the $\tau_e$ summary statistic impacts the resulting constraints on galaxy properties and EoR histories. Using the last Planck 2018 data release, we show that the marginalized constraints on the EoR history are more sensitive to the choice of the basis set (Tanh vs physical model) than to the CMB likelihood statistic ($\tau_e$ vs PS). For example, EoR histories implied by the growth of structure show a small tail of partial reionization extending to higher redshifts: the lower limit on the volume-averaged neutral hydrogen fraction at $z=10$ changes from $\bar{x}_{\rm HI} \ge 0.925 (1 \sigma)$ using Tanh to $\bar{x}_{\rm HI} \ge 0.849 (1 \sigma)$ using the physical model. However, biases in inference using $\tau_e$ are negligible for the Planck 2018 data.

The very first galaxies that started the cosmic dawn likely resided in so-called "minihaloes", with masses of $\sim10^5$-$10^8\mathrm{M}_\odot$, accreting their gas from the intergalactic medium through H$_2$ cooling. Such molecularly cooled galaxies (MCGs) mostly formed in pristine environments, hosted massive, metal-free stars, and were eventually sterilized by the build-up of a disassociating (Lyman-Werner; LW) background. Therefore, their properties might be very different from the galaxies we see in the later Universe. Although MCGs are probably too faint to be observed directly, we could nevertheless infer their properties from the imprint they leave in the cosmic 21-cm signal. Here we quantify this imprint by extending the public simulation code 21cmFAST to allow for a distinct population of MCGs. We allow MCGs to have different properties from other galaxies, including unique scaling relations for their stellar-to-halo mass ratios, ionizing escape fractions, and spectral energy distributions. We track inhomogeneous recombinations, disassociative LW feedback, and photoheating from reionization. After demonstrating how MCGs can shape the 21-cm signal, we explore to what extent current observations can already place constraints on their properties. The cosmic microwave background optical depth from Planck sets an upper limit on the product of the ionizing escape fraction and the stellar mass in MCGs. When including also the timing of the putative EDGES absorption signal, we find an additional strong degeneracy between the stellar mass and the X-ray luminosity of MCGs. If proven to be of cosmic origin, the timing of the EDGES signal would have been set by MCGs.

Correlations between black holes and their host galaxies provide insight into what drives black hole-host co-evolution. We use the Meraxes semi-analytic model to investigate the growth of black holes and their host galaxies from high redshift to the present day. Our modelling finds no significant evolution in the black hole-bulge and black hole-total stellar mass relations out to a redshift of 8. The black hole-total stellar mass relation has similar but slightly larger scatter than the black hole-bulge relation, with the scatter in both decreasing with increasing redshift. In our modelling the growth of galaxies, bulges and black holes are all tightly related, even at the highest redshifts. We find that black hole growth is dominated by instability-driven or secular quasar-mode growth and not by merger-driven growth at all redshifts. Our model also predicts that disc-dominated galaxies lie on the black hole-total stellar mass relation, but lie offset from the black hole-bulge mass relation, in agreement with recent observations and hydrodynamical simulations.

We study the formation of coalescing binary black holes via the evolution of isolated field binaries that go through the common envelope phase in order to obtain the combined distributions of the main observables of Advanced LIGO. We used a hybrid technique that combines the parametric binary population synthesis code COMPAS with detailed binary evolution simulations performed with the MESA code. We then convolved our binary evolution calculations with the redshift- and metallicity-dependent star-formation rate and the selection effects of gravitational-wave detectors to obtain predictions of observable properties. By assuming efficient angular momentum transport, we are able to present a model that is capable of simultaneously predicting the three main gravitational-wave observables: the effective inspiral spin parameter $\chi_{eff}$, the chirp mass $M_{chirp}$, and the cosmological redshift of merger $z_{merger}$. We find an excellent agreement between our model and the ten events from the first two advanced detector observing runs. We make predictions for the third observing run O3 and for Advanced LIGO design sensitivity. We expect around 80% of events with $\chi_{eff} < 0.1$, while the remaining 20% of events with $\chi_{eff} \ge 0.1$ are split into ~10% with $M_{chirp} < 15$ M$_\odot$ and ~10% with $M_{chirp} \ge 15$ M$_\odot$. In conclusion, the favorable comparison of the existing LIGO/Virgo observations with our model predictions gives support to the idea that the majority, if not all of the observed mergers, originate from the evolution of isolated binaries. The first-born black hole has negligible spin because it lost its envelope after it expanded to become a giant star, while the spin of the second-born black hole is determined by the tidal spin up of its naked helium star progenitor by the first-born black hole companion after the binary finished the common-envelope phase.

Intergalactic medium temperature is a powerful probe of the epoch of reionisation, as information is retained long after reionisation itself. However, mean temperatures are highly degenerate with the timing of reionisation, with the amount heat injected during the epoch, and with the subsequent cooling rates. We post-process a suite of semi-analytic galaxy formation models to characterise how different thermal statistics of the intergalactic medium can be used to constrain reionisation. Temperature is highly correlated with redshift of reionisation for a period of time after the gas is heated. However as the gas cools, thermal memory of reionisation is lost, and a power-law temperature-density relation is formed, $T = T_0(1+\delta)^{1-\gamma}$ with $\gamma \approx 1.5$. Constraining our model against observations of electron optical depth and temperature at mean density, we find that reionisation likely finished at $z_{\rm{reion}} = 6.8 ^{+ 0.5} _{-0.8}$ with a soft spectral slope of $\alpha = 2.8 ^{+ 1.2} _{-1.0}$. By restricting spectral slope to the range $[0.5,2.5]$ motivated by population II synthesis models, reionisation timing is further constrained to $z_{\rm{reion}} = 6.9 ^{+ 0.4} _{-0.5}$. We find that, in the future, the degeneracies between reionisation timing and background spectrum can be broken using the scatter in temperatures and integrated thermal history.

We study the sizes, angular momenta and morphologies of high-redshift galaxies using an update of the Meraxes semi-analytic galaxy evolution model. Our model successfully reproduces a range of observations from redshifts $z=0$-$10$. We find that the effective radius of a galaxy disc scales with UV luminosity as $R_e\propto L_{\textrm{UV}}^{0.33}$ at $z=5$-$10$, and with stellar mass as $R_e\propto M_\ast^{0.24}$ at $z=5$ but with a slope that increases at higher redshifts. Our model predicts that the median galaxy size scales with redshift as $R_e \propto (1+z)^{-m}$, where $m=1.98\pm0.07$ for galaxies with $(0.3$-$1)L^\ast_{z=3}$ and $m=2.15\pm0.05$ for galaxies with $(0.12$-$0.3)L^\ast_{z=3}$. We find that the ratio between stellar and halo specific angular momentum is typically less than one and decreases with halo and stellar mass. This relation shows no redshift dependence, while the relation between specific angular momentum and stellar mass decreases by $\sim0.5$ dex from $z=7$ to $z=2$. Our model reproduces the distribution of local galaxy morphologies, with bulges formed predominantly through galaxy mergers for low-mass galaxies, disc-instabilities for galaxies with $M_\ast\simeq10^{10}$-$10^{11.5}M_\odot$, and major mergers for the most massive galaxies. At high redshifts, we find galaxy morphologies that are predominantly bulge-dominated.

Correlation between the variations of quasar absorption lines and the ionizing continuum have been recently confirmed in systematic studies. However, no convincing individual case is reported. We present a statistical analysis of the variable C iv broad absorption line (BAL) in the quasar SDSS J141007.74+541203.3, which have been observed with 44 epochs by the Sloan Digital Sky Survey Data Release 14. Grier et al. (2015) has recently concluded that the most likely cause of the variability of the BAL in SDSS J141007.74+541203.3 is a rapid response to changes in the incident ionizing continuum. In this paper, we confirm the anticorrelation between the equivalent width of BALs and the flux of the continuum based on the spectra of this quasar that show significant variations, which serve as another independent evidence for the view of Grier et al. (2015).

We present an exploration of the expected detection of the earliest Active Galactic Nuclei (AGN) in the Universe from state-of-art galaxy formation and evolution semi-analytic models and hydro-dynamical simulations. We estimate the number and radiative characteristics of Super Massive Black Holes (SMBHs) at $z\geq 6$, a redshift range that will be intensively explored by the next generation of telescopes, in particular in the radio through the Square Kilometre Array (SKA) and at high energies with ESA's Athena X-ray Observatory. We find that Athena will be able to observe over 5000 AGN/$\rm deg^2$ at the Epoch of Re-ionization (EoR), $6\leq z \leq 10$. Similarly, for the same redshift range the models/simulations suggest that SKA will detect at least 400 AGN/$\rm deg^2$. Additionally, we stress the importance of the volume of the simulation box as well as the initial physical conditions of the models/simulations on their effect on the luminosity functions (LFs) and the creation of the most massive SMBHs that we currently observe at the EoR. Furthermore, following the evolution of the accretion mode of the SMBHs in each model/simulation, we show that, while the quasar dominates over the radio mode at the EoR, detection at radio wavelengths still reaches significant numbers even at the highest redshifts. Finally, we present the effect that the radiative efficiency has on the LFs by comparing results produced with a constant value for the radiative efficiency and more complex calculations based on the spin of each SMBH.

We present Hubble Space Telescope WFC3-IR imaging in the fields of six apparently bright dusty star-forming galaxies (DSFGs) at $z$ = 2-4 identified by their rest-frame far-infrared colors using the Planck and Herschel space facilities. We detect near-infrared counterparts for all six submillimeter sources, allowing us to undertake strong-lensing analyses. One field in particular stands out for its prominent giant arcs, PLCK G165.7+67.0 (G165). After combining the color and morphological information, we identify 11 sets of image multiplicities in this one field. We construct a strong-lensing model constrained by this lensing evidence, which uncovers a bimodal spatial mass distribution, and from which we measure a mass of $(2.6 \pm 0.11)$ $\times$ $10^{14}$ $M_{\odot}$ within $\sim$250 kpc. The bright ($S_{350}$ $\approx$ 750 mJy) DSFG appears as two images: a giant arc with a spatial extent of 4.5" that is merging with the critical curve, and a lower-magnification counterimage that is detected in our new longer-wavelength ground- and space-based imaging data. Using our ground-based spectroscopy, we calculate a dynamical mass of $1.3^{+0.04}_{-0.70} \times 10^{15}$ $M_{\odot}$ to the same fixed radius, although this value may be inflated relative to the true value if the velocity distribution is enhanced in the line-of-sight direction. We suggest that the bimodal mass taken in combination with the weak X-ray flux and low SZ decrement may be explained as a pre-merger for which the intracluster gas is diluted along the line of sight, while the integrated surface mass density is supercritical to strong-lensing effects.

Black hole (BH) spins in low-mass X-ray binaries (LMXBs) cover a range of values that can be explained by accretion after BH birth. In contrast, the three BH spin measurements in high-mass X-ray binaries (HMXBs) show only values near the maximum and likely have a different origin connected to the BH stellar progenitor. We explore here two possible scenarios to explain the high spins of BHs in HMXBs: formation in binaries that undergo mass transfer (MT) during the main sequence (MS; Case-A MT), and very close binaries undergoing chemically homogeneous evolution (CHE). We find that both scenarios are able to produce high-spin BHs if internal angular momentum (AM) transport in the progenitor star after its MS evolution is not too strong (i.e., weak coupling between the stellar core and its envelope). If instead efficient AM transport is assumed, we find that the resulting BH spins are always too low with respect to observations. The Case-A MT model provides a good fit for the BH spins, the masses of the two components, and the final orbital periods for two of the three BHs in HMXBs with measured spins. For one of them, the mass predicted for the BH companion is significantly lower than observed, but this depends strongly on the assumed efficiency of mass transfer. The CHE models predict orbital periods that are too large for all three cases considered here. We expect the Case-A MT to be much more frequent at the metallicities involved, so we conclude that the Case-A MT scenario is preferred. Finally, we predict that the stellar companions of HMXBs formed through the Case-A MT have enhanced nitrogen surface abundances, which can be tested by future observations.

Having completed its commissioning phase, the Advanced Rayleigh guided Ground-layer adaptive Optics System (ARGOS) facility is coming online for scientific observations at the Large Binocular Telescope (LBT). With six Rayleigh laser guide stars in two constellations and the corresponding wavefront sensing, ARGOS corrects the ground-layer distortions for both LBT 8.4m eyes with their adaptive secondary mirrors. Under regular observing conditions, this set-up delivers a point spread function (PSF) size reduction by a factor of ~2--3 compared to a seeing-limited operation. With the two LUCI infrared imaging and multi-object spectroscopy instruments receiving the corrected images, observations in the near-infrared can be performed at high spatial and spectral resolution. We discuss the final ARGOS technical set-up and the adaptive optics performance. We show that imaging cases with ground-layer adaptive optics (GLAO) are enhancing several scientific programmes, from cluster colour magnitude diagrams and Milky Way embedded star formation, to nuclei of nearby galaxies or extragalactic lensing fields. In the unique combination of ARGOS with the multi-object near-infrared spectroscopy available in LUCI over a 4x4 arcmin field of view, the first scientific observations have been performed on local and high-z objects. Those high spatial and spectral resolution observations demonstrate the capabilities now at hand with ARGOS at the LBT.

We investigate the dependence of galaxy clustering at $z \sim 4 - 7$ on UV-luminosity and stellar mass. Our sample consists of $\sim$ 10,000 Lyman-break galaxies (LBGs) in the XDF and CANDELS fields. As part of our analysis, the $M_\star - M_{\rm UV}$ relation is estimated for the sample, which is found to have a nearly linear slope of $d\log_{10} M_\star / d M_{\rm UV} \sim 0.44$. We subsequently measure the angular correlation function and bias in different stellar mass and luminosity bins. We focus on comparing the clustering dependence on these two properties. While UV-luminosity is only related to recent starbursts of a galaxy, stellar mass reflects the integrated build-up of the whole star formation history, which should make it more tightly correlated with halo mass. Hence, the clustering segregation with stellar mass is expected to be larger than with luminosity. However, our measurements suggest that the segregation with luminosity is larger with $\simeq 90\%$ confidence (neglecting contributions from systematic errors). We compare this unexpected result with predictions from the \textsc{Meraxes} semi-analytic galaxy formation model. Interestingly, the model reproduces the observed angular correlation functions, and also suggests stronger clustering segregation with luminosity. The comparison between our observations and the model provides evidence of multiple halo occupation in the small scale clustering.

We directly compare predictions of dwarf galaxy properties in a semi-analytic model (SAM) with those extracted from a high-resolution hydrodynamic simulation. We focus on galaxies with halo masses of 1e9<Mvir/Msol<1e11 at high redshift ($z\ge5$). We find that, with the modifications previously proposed in Qin et al. (2018), including to suppress the halo mass and baryon fraction, as well as to modulate gas cooling and star formation efficiencies, the SAM can reproduce the cosmic evolution of galaxy properties predicted by the hydrodynamic simulation. These include the galaxy stellar mass function, total baryonic mass, star-forming gas mass and star formation rate at $z\sim5-11$. However, this agreement is only possible by reducing the star formation threshold relative to that suggested by local observations. Otherwise, too much star-forming gas is trapped in quenched dwarf galaxies. We further find that dwarf galaxies rapidly build up their star-forming reservoirs in the early universe ($z>10$), with the relevant time-scale becoming significantly longer towards lower redshifts. This indicates efficient accretion in cold mode in these low-mass objects at high redshift. Note that the improved SAM, which has been calibrated against hydrodynamic simulations, can provide more accurate predictions of high-redshift dwarf galaxy properties that are essential for reionization study.

We study dwarf galaxy formation at high redshift ($z\ge5$) using a suite of high- resolution, cosmological hydrodynamic simulations and a semi-analytic model (SAM). We focus on gas accretion, cooling and star formation in this work by isolating the relevant process from reionization and supernova feedback, which will be further discussed in a companion paper. We apply the SAM to halo merger trees constructed from a collisionless N-body simulation sharing identical initial conditions to the hydrodynamic suite, and calibrate the free parameters against the stellar mass function predicted by the hydrodynamic simulations at z = 5. By making comparisons of the star formation history and gas components calculated by the two modelling techniques, we find that semi-analytic prescriptions that are commonly adopted in the literature of low-redshift galaxy formation do not accurately represent dwarf galaxy properties in the hydrodynamic simulation at earlier times. We propose 3 modifications to SAMs that will provide more accurate high-redshift simulations. These include 1) the halo mass and baryon fraction which are overestimated by collisionless N-body simulations; 2) the star formation efficiency which follows a different cosmic evolutionary path from the hydrodynamic simulation; and 3) the cooling rate which is not well defined for dwarf galaxies at high redshift. Accurate semi-analytic modelling of dwarf galaxy formation informed by detailed hydrodynamical modelling will facilitate reliable semi-analytic predictions over the large volumes needed for the study of reionization.

Various binary black hole formation channels have been proposed since the first gravitational event GW150914 was discovered by the Advanced Laser Interferometer Gravitational-Wave Observatory (AdLIGO). For all evolutionary channels based on the evolution of isolated binaries, the immediate progenitor of the binary black hole is a close binary system composed of a black hole and a helium star. We perform detailed binary evolution and systematically explore the parameter space of initial binary properties, including initial black hole and helium star masses, initial rotation of the helium star as well as metallicity. We argue that the spin of the first-born black hole at its birth is negligible ($\lesssim 0.1$), hence the second-born black hole's spin dominates the measured effective spin, $\chi_{\rm eff}$, from gravitational wave events of double black hole mergers. We find that tides can be important only when orbital periods are shorter than 2 days. Upon core collapse, the helium star produces a black hole (the second-born black hole in the system) with a spin that can span the entire range from zero to maximally spinning. We show that the bimodal distribution of the spin of the second-born black hole obtained in recent papers is mainly due to oversimplifying assumptions. We find an anti-correlation between the merging timescale of the two black holes, T$_{\rm merger}$, and the effective spin $\chi_{\rm eff}$. Finally, we provide new prescriptions for the tidal coefficient E$_2$ for both H-rich and the helium-rich stars. We predict that, with future improvements to AdLIGO's sensitivity, the sample of merging binary black hole systems will show an overdensity of sources with positive but small $\chi_{\rm eff}$. originating from lower mass black hole mergers born at low redshift.

We present the ultraviolet magnitudes for over three million stars in the LAMOST survey, in which 2,202,116 stars are detected by $GALEX$. For 889,235 undetected stars, we develop a method to estimate their upper limit magnitudes. The distribution of (FUV $-$ NUV) shows that the color declines with increasing effective temperature for stars hotter than 7000 K in our sample, while the trend disappears for the cooler stars due to upper atmosphere emission from the regions higher than their photospheres. For stars with valid stellar parameters, we calculate the UV excesses with synthetic model spectra, and find that the (FUV $-$ NUV) vs. $R'_{\mathrm{FUV}}$ can be fitted with a linear relation and late-type dwarfs tend to have high UV excesses. There are 87,178 and 1,498,103 stars detected more than once in the visit exposures of $GALEX$ in the FUV and NUV, respectively. We make use of the quantified photometric errors to determine statistical properties of the UV variation, including intrinsic variability and the structure function on the timescale of days. The overall occurrence of possible false positives is below 1.3\% in our sample. UV absolute magnitudes are calculated for stars with valid parallaxes, which could serve as a possible reference frame in the NUV. We conclude that the colors related to UV provide good criteria to distinguish between M giants and M dwarfs, and the variability of RR Lyrae stars in our sample is stronger than that of other A and F stars.

In this Letter, we present an analysis of the relation between the variability of broad absorption lines (BALs) and that of the continuum. Our sample is multi-epoch observations of 483 quasars by the Sloan Digital Sky Survey-I/II/III (SDSS-I/II/III). We derive the fractional flux variations of the continuum and fractional equivalent width (EW) variations for C IV and Si IV BALs, and explore the correlations between the three. Our results reveal moderate anticorrelations with high significance level between the fractional flux variations of the continuum and fractional EW variations for both C IV and Si IV BALs. We also prove a significant positive correlation between the fractional EW variations for C IV and Si IV BALs, which is in agreement with several previous studies. Our discoveries can serve as evidence for the idea: Change of an ionizing continuum is the primary driver of BAL variability.

The observations and research on the neutrinos provide a kind of indirect way of revealing the properties of dark matter particles. For the detection of muon neutrinos, the main issue is the large atmospheric background, which is caused by the interactions between the cosmic rays and atoms within the atmosphere. Compared with muon neutrinos, tau neutrinos have a smaller atmospheric background especially for the downward-going direction. Except for the classical neutrino sources, dark matter particles can also annihilate into the neutrinos and are the potential high energy astrophysical sources. The annihilation rate of dark matter particles is proportional to the square of number density; therefore, the annihilation rate is large near the center of dark matter halos especially for the new kind of dark matter structures named ultracompact dark matter minihalos (UCMHs). In previous works, we have investigated the potential muon neutrino flux from UCMHs due to dark matter annihilation. Moreover, since the formation of UCMHs is related to the primordial density perturbations of small scales, we get the constraints on the amplitude of the primordial curvature perturbations of small scales, $1 \lesssim k \lesssim 10^{7} ~\rm Mpc^{-1}$. In this work, we focus on the downward-going tau neutrinos from UCMHs due to dark matter annihilation. Compared with the background of tau neutrino flux we get the constraints on the mass fraction of UCMHs. Then using the limits on the mass fraction of UCMHs we got the constraints on the amplitude of the primordial curvature perturbations which are extended to the scale $k \sim 10^{8} ~ \rm Mpc^{-1}$ compared with previous results.

Motivated by recent measurements of the number density of faint AGN at high redshift, we investigate the contribution of quasars to reionization by tracking the growth of central supermassive black holes in an update of the Meraxes semi-analytic model. The model is calibrated against the observed stellar mass function at $z\sim0.6-7$, the black hole mass function at $z\lesssim0.5$, the global ionizing emissivity at $z\sim2-5$ and the Thomson scattering optical depth. The model reproduces a Magorrian relation in agreement with observations at $z<0.5$ and predicts a decreasing black hole mass towards higher redshifts at fixed total stellar mass. With the implementation of an opening angle of 80 deg for quasar radiation, corresponding to an observable fraction of ${\sim}23.4$ per cent due to obscuration by dust, the model is able to reproduce the observed quasar luminosity function at $z\sim0.6-6$. The stellar light from galaxies hosting faint AGN contributes a significant or dominant fraction of the UV flux. At high redshift, the model is consistent with the bright end quasar luminosity function and suggests that the recent faint $z\sim4$ AGN sample compiled by Giallongo et al. (2015) includes a significant fraction of stellar light. Direct application of this luminosity function to the calculation of AGN ionizing emissivity consequently overestimates the number of ionizing photons produced by quasars by a factor of 3 at $z\sim6$. We conclude that quasars are unlikely to make a significant contribution to reionization.

Massive stars with solar metallicity lose important amounts of rotational angular momentum through their winds. When a magnetic field is present at the surface of a star, efficient angular momentum losses can still be achieved even when the mass-loss rate is very modest, at lower metallicities, or for lower-initial-mass stars. In a close binary system, the effect of wind magnetic braking also interacts with the influence of tides, resulting in a complex evolution of rotation. We study the interactions between the process of wind magnetic braking and tides in close binary systems. We discuss the evolution of a 10 M$_\odot$ star in a close binary system with a 7 M$_\odot$ companion using the Geneva stellar evolution code. The initial orbital period is 1.2 days. The 10 M$_\odot$ star has a surface magnetic field of 1 kG. Various initial rotations are considered. We use two different approaches for the internal angular momentum transport. In one of them, angular momentum is transported by shear and meridional currents. In the other, a strong internal magnetic field imposes nearly perfect solid-body rotation. The evolution of the primary is computed until the first mass-transfer episode occurs. The cases of different values for the magnetic fields and for various orbital periods and mass ratios are briefly discussed. We show that, independently of the initial rotation rate of the primary and the efficiency of the internal angular momentum transport, the surface rotation of the primary will converge, in a time that is short with respect to the main-sequence lifetime, towards a slowly evolving velocity that is different from the synchronization velocity. (abridged).

Massive quiescent galaxies (MQGs) are thought to have formed stars rapidly at early times followed by a long period of quiescence. The recent discovery of a MQG, ZF-COSMOS-20115 at $z\sim4$, only 1.5 Gyr after the big bang, places new constraints on galaxy growth and the role of feedback in early star formation. Spectroscopic follow-up confirmed ZF-COSMOS-20115 as a MQG at $z=3.717$ with an estimated stellar mass of ${\sim}10^{11}\mathrm{M}_\odot$, showing no evidence of recent star formation. We use the Meraxes semi-analytic model to investigate how ZF-COSMOS-20115 analogues build stellar mass, and why they become quiescent. We identify three analogue galaxies with similar properties to ZF-COSMOS-20115. We find that ZF-COSMOS-20115 is likely hosted by a massive halo with virial mass of ${\sim}10^{13}\mathrm{M}_\odot$, having been through significant mergers at early times. These merger events drove intense growth of the nucleus, which later prevented cooling and quenched star formation. ZF-COSMOS-20115 likely remained quiescent at $z<3.7$. We find that the analogues host the most massive black holes in our simulation and were luminous quasars at $z\sim5$, indicating that ZF-COSMOS-20115 and other MQGs may be the descendants of high-redshift quasars. In addition, the model suggests that ZF-COSMOS-20115 formed in a region of intergalactic medium that was reionized early.

In this Letter, we report the discovery of a strong correlation between the variability of narrow absorption lines (NALs) and the ionizing continuum from a two-epoch spectra sample of 40 quasars containing 52 variable C iv {\lambda}{\lambda}1548; 1551 absorption doublets. According to the concordance index, this sample is classified into two subsamples. Subsample I shows an anti-correlation between the variations of absorption lines and the continuum, while Subsample II exhibits a positive correlation. These results imply that these variable C iv {\lambda}{\lambda}1548; 1551 absorption doublets are intrinsic to the corresponding quasars and that their variations are caused primarily by the fluctuations of the ionizing continuum. Based on our analysis, we propose that there might be two kinds of absorption gas: one that is very sensitive to the continuum variations, the another that is not. In addition, we suggest that in many cases the emergence or disappearance of NALs is caused by fluctuations of the ionizing continuum.

We investigate how the hydrostatic suppression of baryonic accretion affects the growth rate of dark matter halos during the Epoch of Reionization. By comparing halo properties in a simplistic hydrodynamic simulation in which gas only cools adiabatically, with its collisionless equivalent, we find that halo growth is slowed as hydrostatic forces prevent gas from collapsing. In our simulations, at the high redshifts relevant for reionization (between ${\sim}6$ and ${\sim}11$), halos that host dwarf galaxies ($\lesssim 10^{9} \mathrm{M_\odot}$) can be reduced by up to a factor of 2 in mass due to the hydrostatic pressure of baryons. Consequently, the inclusion of baryonic effects reduces the amplitude of the low mass tail of the halo mass function by factors of 2 to 4. In addition, we find that the fraction of baryons in dark matter halos hosting dwarf galaxies at high redshift never exceeds ${\sim}90\%$ of the cosmic baryon fraction. When implementing baryonic processes, including cooling, star formation, supernova feedback and reionization, the suppression effects become more significant with further reductions of ${\sim}30\%$ to 60\%. Although convergence tests suggest that the suppression may become weaker in higher resolution simulations, this suppressed growth will be important for semi-analytic models of galaxy formation, in which the halo mass inherited from an underlying N-body simulation directly determines galaxy properties. Based on the adiabatic simulation, we provide tables to account for these effects in N-body simulations, and present a modification of the halo mass function along with explanatory analytic calculations.

Using the LAMOST-Gaia common stars, we demonstrate that the in-plane velocity field for the nearby young stars are significantly different from that for the old ones. For the young stars, the probably perturbed velocities similar to the old population are mostly removed from the velocity maps in the $X$--$Y$ plane. The residual velocity field shows that the young stars consistently move along $Y$ with faster $v_\phi$ at the trailing side of the local arm, while at the leading side, they move slower in azimuth direction. At both sides, the young stars averagely move inward with $v_R$ of $-5\sim-3$ km s$^{-1}$. The divergence of the velocity in $Y$ direction implies that the young stars are associated with a density wave nearby the local arm. We therefore suggest that the young stars may reflect the formation of the local spiral arm by correlating themselves with a density wave. The range of the age for the young stars is around 2 Gyr, which is sensible since the transient spiral arm can sustain that long. We also point out that alternative explanations of the peculiar velocity field for the young population cannot be ruled out only from this observed data.

We explore the kinematics (both the radial velocity and the proper motion) of the vertical X-shaped feature in the Milky Way with an N-body bar/bulge model. From the solar perspective, the distance distribution of particles is double-peaked in fields passing through the X-shape. The separation and amplitude ratio between the two peaks qualitatively match the observed trends towards the Galactic bulge. We confirm clear signatures of cylindrical rotation in the pattern of mean radial velocity across the bar/bulge region. We also find possible imprints of coherent orbital motion inside the bar structure in the radial velocity distribution along l=0 degree, where the near and far sides of the bar/bulge show excesses of approaching and receding particles. The coherent orbital motion is also reflected in the slight displacement of the zero-velocity-line in the mean radial velocity, and the displacement of the maximum/minimum in the mean longitudinal proper motion across the bulge region. We find some degree of anisotropy in the stellar velocity within the X-shape, but the underlying orbital family of the X-shape cannot be clearly distinguished. Two potential applications of the X-shape in previous literature are tested, i.e., bulge rotation and Galactic center measurements. We find that the proper motion difference between the two sides of the X-shape can be used to estimate the mean azimuthal streaming motion of the bulge, but not the pattern speed of the bar. We also demonstrate that the Galactic center can be located with the X-shape, but the accuracy depends on the fitting scheme, the number of fields, and their latitudinal coverage.

As the second work in a series of papers aiming to detect absorption systems in the quasar spectra of the Baryon Oscillation Spectroscopic Survey, we continue the analysis of Paper I by expanding the quasar sample to those quasars with zem>2.4. This yields a sample of 21,963 appropriate quasars to search for narrow C IV absorptions with Wr>=0.2A for both lines. There are 9708 quasars with at least one appropriate absorption system imprinted on their spectra. From these spectra, we detect 13,919 narrow Civabsorption systems whose absorption redshifts cover a range of zabs=1.8784 - 4.3704. In this paper and Paper 1, we have selected 37,241 appropriate quasars with median SNR>= 4 and 1.54<= zem <=5.16 to visually analyze narrow C IV absorption doublets one by one. A total of 15,999 quasars are found to have at least one appropriate absorption system imprinted on their spectra. From these 15,999 quasar spectra, we have detected 23,336 appropriate C IV absorption systems with Wr>=0.2A whose absorption redshifts cover a range of zabs=1.4544 - 4.3704. The largest values of Wr are 3.19 A for the \lambda1548 absorption line and 2.93 A for the lambda1551 absorption line, respectively. We find that only a few absorbers show large values of Wr. About 1.1% absorbers of the total absorbers have Wr>=2.0 A.

We have assembled absorption systems by visually identifying $C~IV\lambda\lambda1548,1551$ absorption doublets in the quasar spectra of the Baryon Oscillation Spectroscopic Survey (BOSS) one by one. This paper is the first of the series work. In this paper, we concern quasars with relatively low redshifts and high signal-to-noise ratios for their spectra, and hence we limit our analysis on quasars with $z_{em}\le2.4$ and on the doublets with $W_r\lambda1548\ge0.2$ \AA. Out of the more than 87,000 quasars in the Data Release 9, we limit our search to 10,121 quasars that have the appropriate redshifts and spectra with high enough signal-to-noise ratios to identify narrow C IV absorption lines. Among them, 5,442 quasars are detected to have at least one $C~IV\lambda\lambda1548,1551$ absorption doublet. We obtain a catalog containing 8,368 $C~IV\lambda\lambda1548,1551$ absorption systems, whose redshifts are within $z_{abs}=1.4544$ - $2.2805$. In this catalog, about $33.7\%$ absorbers have $0.2$ \AA$\le W_r\lambda1548<0.5$ \AA, about $45.9\%$ absorbers have $0.5$ \AA$\le W_r\lambda1548<1.0$ \AA, about $19.2\%$ absorbers have $1.0$ \AA$\le W_r\lambda1548<2.0$ \AA, and about $1.2\%$ absorbers have $W_r\lambda1548\ge2.0$ \AA.

The $L_p - E_p$ relation is a well-known relation in gamma-ray bursts. Its implication remains unclear. We propose to investigate the underlying mechanisms of this relation by considering the corresponding kinetic and dynamic mechanisms separately. In this way, one can tell how much the kinetic or dynamic mechanism contributes to the index of the relationship. Our analysis gives rise to several conclusions. (1) The index of the kinetic effect in the $L_p - E_p$ relation can simply be derived from the theory of special relativity, which is generally larger than 2, depending on the situation concerned. (2) The index of the dynamic effect in the relation can be deduced from observation once a model of jets is adopted. According to current GRB data, we find: the dynamic effect alone tends to make an anti-correlation between $L_p$ and $E_p$; in terms of statistics, the dynamic effect is obviously smaller than the kinetic effect; in the situation of jets with moving discrete radio clouds which moves directly towards the observer, the index of the dynamic effect is currently constrained within $(-1.6, -1)$, while in other situations of jets, the constrains are different; both internal and external shocks can account for the current data.

In this letter, we present for the first time, the discovery of the disappearance of a narrow Mg II $\lambda\lambda2796,2803$ absorption system from the spectra of quasar SDSS J165501.31+260517.4 ($z_{\rm e}=1.8671$). This absorber is located at $z_{\rm abs} =1.7877$, and has a velocity offset of $8,423\rm ~km~s^{-1}$ with respect to the quasar. According to the velocity offset and the line variability, this narrow Mg II $\lambda\lambda2796,2803$ absorption system is likely intrinsic to the quasar. Since the corresponding UV continuum emission and the absorption lines of another narrow Mg II $\lambda\lambda2796,2803$ absorption system at $z_{\rm abs}=1.8656$ are very stable, we think that the disappearance of the absorption system is unlikely to be caused by the change in ionization of absorption gas. Instead, it likely arises from the motion of the absorption gas across the line of sight.

Gamma-ray bursts (GRBs) are believed to belong to two classes and they are conventionally divided according to their durations. This classification scheme is not satisfied due to the fact that duration distributions of the two classes are heavily overlapped. We collect a new sample (153 sources) of GRBs and investigate the distribution of the logarithmic deviation of the $E_p$ value from the Amati relation. The distribution possesses an obvious bimodality structure and it can be accounted for by the combination of two Gaussian curves. Based on this analysis, we propose to statistically classify GRBs in the well-known $E_{p}$ vs. $E_{iso}$ plane with the logarithmic deviation of the $E_p$ value. This classification scheme divides GRBs into two groups: Amati type bursts and non-Amati type bursts. While Amati type bursts well follow the Amati relation, non-Amati type bursts do not. It shows that most Amati type bursts are long duration bursts and the majority of non-Amati type bursts are short duration bursts. In addition, it reveals that Amati type bursts are generally more energetic than non-Amati type bursts. An advantage of the new classification is that the two kinds of burst are well distinguishable and therefore their members can be identified in certainty.

Quasar J152645.61+193006.7 ($z_{\rm e}=2.5771$) was observed by the Sloan Digital Sky Survey (SDSS) on 2006 May 31, and again on 2011 April 9. The time interval of the SDSS two observations is 497.4 days at the quasar rest frame. From the spectra of this quasar, we detect a phenomenon of disappearance of a narrow $\rm Mg II \lambda\lambda2796,2803$ absorption system with a velocity up to $\rm 166,129 km s^{-1}$ with respect to the quasar. This disappearance event can be caused by changes in the ionization state of absorbing gas or by the bulk motion of the absorbing gas across the quasar sightline. The coverage fraction analysis shows that this absorber partially covers the background emission sources with an effective coverage fraction of $C_{\rm f}=0.40\pm0.06$. The time variation analysis and the coverage fraction analysis imply that this absorber might be intrinsic to the quasar. However, the scenario of a cosmologically separated foreground object located at $z=0.9170$ accounting for the phenomenon cannot be ruled out according to current available data.

We explore a second order Hamiltonian vertical resonance model for X-shaped or peanut-shaped galactic bulges. The X-shape is caused by the 2:1 vertical Lindblad resonance with the bar, with two vertical oscillation periods per orbital period in the bar frame. We examine N-body simulations and find that due to the bar slowing down and disk thickening during bar buckling, the resonance and associated peanut-shape moves outward. The peanut-shape is consistent with the location of the vertical resonance, independent of whether the bar buckled or not. We estimate the resonance width from the potential m=4 Fourier component and find that the resonance is narrow, affecting orbits over a narrow range in the angular momentum distribution, dL/L ~ 0.05. As the resonance moves outward, stars originally in the mid plane are forced out of the mid plane into orbits just within the resonance separatrix. The height of the separatrix orbits, estimated from the Hamiltonian model, is approximately consistent with the peanut-shape height. The X-shape is comprised of stars in the vicinity of the resonance separatrix. The velocity distributions from the simulations illustrate that low inclination orbits are depleted within resonance. Within resonance, the vertical velocity distribution is broad, consistent with resonant heating caused by the passage of the resonance through the disk. In the Milky Way bulge we relate the azimuthally averaged mid-plane mass density near the vertical resonance to the rotation curve and bar pattern speed. At an estimated vertical resonance galactocentric radius of ~1.3 kpc, we confirm a mid-plane density of ~5x10^8 Msol/kpc^3, consistent with recently estimated mass distributions. We find that the rotation curve, bar pattern speed, 2:1 vertical resonance location, X-shape tips, and mid-plane mass density, are all self-consistent in the Milky Way galaxy bulge.

The durations (T90) of 315 GRBs detected with Fermi/GBM (8-1000 keV) by 2011 September are calculated using the Bayesian Block method. We compare the T90 distributions between this sample and those derived from previous/current GRB missions. We show that the T90 distribution of this GRB sample is bimodal, with a statistical significance level being comparable to those derived from the BeppoSAX/GRBM sample and the Swift/BAT sample, but lower than that derived from the CGRO/BATSE sample. The short-to-long GRB number ratio is also much lower than that derived from the BATSE sample, i.e., 1:6.5 vs 1:3. We measure T90 in several bands, i.e., 8-15, 15-25, 25-50, 50-100, 100-350, and 350-1000 keV, to investigate the energy-dependence effect of the bimodal T90 distribution. It is found that the bimodal feature is well observed in the 50-100 and 100-350 keV bands, but is only marginally acceptable in the 25-50 keV and 350-1000 keV bands. The hypothesis of the bimodality is confidently rejected in the 8-15 and 15-25 keV bands. The T90 distributions in these bands are roughly consistent with those observed by missions with similar energy bands. The parameter T90 as a function of energy follows \bar T90 \propto E^{-0.20\pm 0.02} for long GRBs. Considering the erratic X-ray and optical flares, the duration of a burst would be even much longer for most GRBs. Our results, together with the observed extended emission of some short GRBs, indicate that the central engine activity time scale would be much longer than T90} for both long and short GRBs and the observed bimodal T90 distribution may be due to an instrumental selection effect.

The \gamma-ray loud blazars (flat spectrum radio quasars--FSRQs and BL Lacertae objects-BLs) are very bright in the \gamma-ray bands, which is perhaps associated with a beaming effect. Therefore, one can expect that the \gamma-ray luminosity is correlated with the beaming factor. In this paper, we investigated the relation between the radio Doppler factors and the gamma-ray luminosities. Our analysis suggests that the \gamma-ray luminosity be strongly correlated with the factor of \delta_R for the whole sample, FSRQs, and BLs. When the effect of a common redshift is excluded, the correlation still exists for the FSRQs sub-sample suggesting that the \gamma-rays are strongly beamed. However, the partial correlation analysis does not show a correlation for the small BL Lac sample.

Using the gamma-ray bursts simultaneously detected by Swift/BAT and Fermi/GBM we performed a joint spectral and temporal analysis of the prompt emission data and confirm the rough correlation between the BAT-band photon index Gamma_BAT and the peak spectral energy Epeak. With the redshift known sub-sample, we derived the isotropic gamma-ray energy E_gamma,iso and also confirm the E_gamma,iso - Epeak,rest relation, with a larger scatter than the Amati sample but consistent with GBM team analyses. We also compare the T_90 values derived in the GBM band with those derived in the BAT band and find that for long GRBs the BAT T_90 is usually longer than the GBM T_90, while for short GRBs the trend reverses. This is consistent with the soft/hard nature of long/short GRBs and suggests the importance of an energy-dependent temporal analysis of GRBs.

Gamma-ray bursts (GRBs) are divided into two classes according to their durations. We investigate if the softness of bursts plays a role in the conventional classification of the objects. We employ the BATSE (Burst and Transient Source Experiment) catalog and analyze the duration distributions of different groups of GRBs associated with distinct softness. Our analysis reveals that the conventional classification of GRBs with the duration of bursts is influenced by the softness of the objects. There exits a bimodality in the duration distribution of GRBs for each group of bursts and the time position of the dip in the bimodality histogram shifts with the softness parameter. Our findings suggest that the conventional classification scheme should be modified by separating the two well-known populations in different softness groups, which would be more reasonable than doing so with a single sample. According to the relation between the dip position and the softness parameter, we get an empirical function that can roughly set apart the short-hard and long-soft bursts: $SP = (0.100 \pm 0.028) T_{90}^{-(0.85 \pm 0.18)}$, where $SP$ is the softness parameter adopted in this paper.

Spectral evolution of gamma-ray burst pulses assumed to arise from emission of fireballs is explored. It is found that, due to the curvature effect, the integrated flux and peak energy are well related by a power law in the decaying phase of pulses, where the index is about 3, being free of the intrinsic emission and the Lorentz factor. The spectrum of a pulse in its decaying phase differs slightly for different intrinsic spectral evolution patterns, indicating that it is dominated by the curvature effect. In the rising phase, the integrated flux keeps increasing whilst the peak energy remains unchanged when the intrinsic emission bears an unchanged spectrum. Within this phase, the flux decreases with the increasing of the peak energy for a hard-to-soft intrinsic spectrum, and for a soft-to-hard-to-soft intrinsic spectrum, the flux generally increases with the increasing of the peak energy. An intrinsic soft-to-hard-to-soft spectral evolution within a co-moving pulse would give rise to a pulse-like evolutionary curve for the peak energy.

The softening process observed in the steep decay phase of early X-ray afterglows of Swift bursts has remained a puzzle since its discovery. The softening process can also be observed in the later phase of the bursts and its cause has also been unknown. Recently, it was suggested that, influenced by the curvature effect, emission from high latitudes would shift the Band function spectrum from higher energy band to lower band, and this would give rise to the observed softening process accompanied by a steep decay of the flux density. The curvature effect scenario predicts that the terminating time of the softening process would be correlated with the duration of the process. In this paper, based on the data from the UNLV GRB group web-site, we found an obvious correlation between the two quantities. In addition, we found that the softening process can be divided into two classes: the early type softening ($t_{s,max}\leq "4000"s$) and the late type softening ($t_{s,max} > "4000"s$). The two types of softening show different behaviors in the duration vs. terminating time plot. In the relation between the variation rates of the flux density and spectral index during the softening process, a discrepancy between the two types of softening is also observed. According to their time scales and the discrepancy between them, we propose that the two types are of different origins: the early type is of internal shock origin and the late type is of external shock origin. The early softening is referred to the steep decay just following the prompt emission, whereas the late decay typically conceives the transition from flat decay to late afterglow decay. We suspect that there might be a great difference of the Lorentz factor in two classes which is responsible for the observed discrepancy.

The curvature effect is explored in the case of extremely short intrinsic emission. Assuming a $\delta$ function emission we get formulas that get rid of the impacts from the intrinsic emission duration, which are applicable to any forms of continuum. The formulas predict that the same form of spectrum could be observed at different times, with the peak energy of the spectrum shifting from higher energy bands to lower bands following $E_{peak}\propto t^{-1}$. When the emission is early enough the light curve in the form $f_{\nu }(t)t^{2}$ will possess exactly the intrinsic spectral form, for which the temporal power law index and the spectral power law index will be related by $\alpha =2+\beta $. The analysis shows that there do exist a temporal steep decay phase and a spectral softening which occur simultaneously, and both are caused by the shifting of the Band function spectrum. According to the analysis, we predict that the softening due to the curvature effect will appear at different frequencies; it occurs earlier for higher frequencies and later for lower frequencies; the maximum spectral index time follows the $t_{b,max}\propto \nu^{-1}$ law. We also find: the softening duration would be linearly correlated with the maximum spectral index time; the observed $\beta_{min}$ and $\beta_{max}$ are determined by the low and high energy indexes of the observed Band function spectrum. We propose to check the curvature effect with the $\log f_{\nu}(t)t^{3}$ vs. $log t$ curve which would be in agreement with the $\log \nu f_{\nu}$ vs. $log \nu$ curve. Applying this to GRB 060614 shows that the peak energy of its observed spectrum is expected to pass through the observation band at $\sim $175 s.

GRB 060124 is the first event that both prompt and afterglow emission were observed simultaneously by the three \emph{Swift} instruments. Its main peak also triggered Konus-Wind and HETE-II. Therefore, investigation on both the temporal and spectral properties of the prompt emission can be extended to X-ray bands. We perform a detailed analysis on the two well identified pulses of this burst, and find that the pulses are narrower at higher energies, and both X-rays and gamma-rays follow the same $w - E$ relation for an individual pulse. However, there is no a universal power-law index of the $w - E$ relation among pulses. We find also that the rise-to-decay ratio $r/d$ seems not to evolve with $E$ and the $r/d$ values are well consistent with that observed in typical GRBs. The broadband spectral energy distribution also suggest that the X-rays are consistent with the spectral behavior of the gamma-rays. These results indicates that the X-ray emission tracks the gamma-ray emission and the emissions in the two energy bands are likely to be originated from the same physical mechanism.

We explore the influence of the full curvature effect on the flux of early X-ray afterglow of gamma-ray bursts (GRBs) in cases when the spectrum of the intrinsic emission is a power-law. We find that the well-known $t^{-(2+\beta)}$ curve is present only when the intrinsic emission is extremely short or the emission arises from an exponential cooling. The time scale of this curve is independent of the Lorentz factor. The resulting light curve would contain two phases when the intrinsic emission has a power-law spectrum and a temporal power-law profile with infinite duration. The first phase is a rapid decay one where the light curve well follows the $t^{-(2+\beta)}$ curve. The second is a shallow decay phase where the power-law index of the light curve is obviously smaller than that in the first phase. The start of the shallow phase is strictly constrained by the fireball radius, which in turn, can put a lower limit to the latter. In the case when the temporal power-law emission lasts a limited interval of time, there will be a third phase after the $t^{-(2+\beta)}$ curve and the shallow decay phase, which is much steeper than the shallow phase. As an example of application, we fit the XRT data of GRB 050219A with our model and show that the curvature effect alone can roughly account for this burst. Although fitting parameters can not be uniquely determined due to various choices of fitting, a lower limit of the fireball radius of this burst can be estimated, which is $\sim 10^{14}cm$.

In this paper a relative number density parameter, called the neighborhood function, is introduced so that the crowded nature of the neighborhood of individual sources can be described. With this parameter one can determine the probability of forming a cluster by chance. A method is proposed to identify large-scale structure in the cosmological comoving frame. The scale used to sort out clustered sources is determined by the mean local number density, and the sample adopted is required to have regular borders. The method is applied to a sample drawn from the Two Degree Field survey. We find from the analysis that the probability of forming the resulting large-scale structures by chance is very small, and the phenomenon of clustering is dominant in the local universe. Within a 3\sigma confidence level, a coherent cluster with a scale as large as 357h^{-1}Mpc is identified from the sample. There exist some galaxies which are not affected by the gravitation of clusters, and hence are suspected to be at rest in the comoving frame of the universe. Voids are likely volumes within which very crowded sources are not present, and they are likely formed in embryo by fluctuations in the very early epochs of the universe. In addition, we find that large-scale structures are coral-like and are likely made up of smaller structures; sources with large values of the neighborhood function are mainly distributed within the structure of prominent clusters, and they form the framework of the large-scale structure.

We employed a sample compiled by Norris et al. (2005, ApJ, 625, 324) to study the dependence of the pulse temporal properties on energy in long-lag, wide-pulse gamma-ray bursts. Our analysis shows that the pulse peak time, rise time scale and decay time scale are power law functions of energy, which is a preliminary report on the relationships between the three quantities and energy. The power law indexes associated with the pulse width, rise time scale and decay time scale are correlated and the correlation between the indexes associated with the pulse width and the decay time scale is more obvious. In addition, we have found that the pulse peak lag is strongly correlated with the CCF lag, but the centroid lag is less correlated with the peak lag and CCF lag. Based on these results and some previous investigations, we tend to believe that all energy-dependent pulse temporal properties may come from the joint contribution of both the hydrodynamic processes of the outflows and the curvature effect, where the energy-dependent spectral lag may be mainly dominated by the dynamic process and the energy-dependent pulse width may be mainly determined by the curvature effect.

We have investigated the gamma-ray bursts pulses assumed to arise from relativistically expending fireballs. We find, due to the curvature effect, the evolutionary curve of the raw hardness ratio (when the background count is included) varies significantly, depending on the hardness of a burst. For a soft burst, an upside-down pulse of its raw hardness ratio curve (HRC) would be observed; for a hard burst, its raw HRC shows a pulse-like profile with a sinkage in its decaying phase; for a very hard burst, the raw HRC possesses a pulse-like profile without a sinkage in its decaying phase. For a pulse-like raw HRC as shown in the case of the hard and very hard bursts, its peak would appear in advance of that of the corresponding light curve. We have studied here the HRC of GRB 920216, GRB 920830 and GRB 990816 in detail. The features of the raw HRC expected in the hard burst are observed in these bursts. A fit to the three bursts shows that the curvature effect alone could indeed account for the predicted characteristics of HRCs. In addition, we find that the observed hardness ratio tends to be harder at the beginning of the pulses than what the curvature effect could predict and be softer at the late time of the pulses. We believe this is an evidence showing the existence of intrinsic hard-to-soft radiation.

Recently, under the assumption that the Doppler effect of the relativistically expanding fireball surface is important, Qin et al. showed that in most cases the power law relationship between the pulse width and energy of gamma-ray bursts (GRBs)would exist in a certain energy range. We check this prediction with two GRB samples which contain well identified pulses. A power law anti-correlation between the full pulse width and energy and a power law correlation between the pulse width ratio and energy are seen in the light curves of the majority (around 65%) of bursts of the two samples within the energy range of BATSE, suggesting that these bursts are likely to arise from the emission associated with the shocks occurred on a relativistically expanding fireball surface. For the rest of the bursts, the relationships between these quantities were not predicted previously. We propose to consider other spectral evolutionary patterns or other radiation mechanisms such as a varying synchrotron or Comptonized spectrum to check if the observed relationships for these rest bursts can also be accounted for by the Doppler model. In addition, we find that the upper limits of the width ratio for the two samples do not exceed 0.9, in agrement with what predicted previously by the Doppler model. The plateau/power law/plateau and the peaked features predicted and detected previously by Qin et al. are generally observed, with the exceptions being noticed only in a few cases. According to the distinct values of two power law indices of FWHM and ratio and energy, we divide the bursts into three subsets which are located in different areas of the two indices plane. We suspect that different locations of the two indices might correspond to different mechanisms.

We present in this paper an approach to estimate the initial Lorentz factor of gamma-ray bursts (GRBs) without referring to the delayed emission of the early afterglow. Under the assumption that the afterglow of the bursts concerned occurs well before the prompt emission dies away, the Lorentz factor measured at the time when the duration of the prompt emission is ended could be estimated by applying the well-known relations of GRB jets. With the concept of the efficiency for converting the explosion energy to radiation, this Lorentz factor can be related to the initial Lorentz factor of the source. The corresponding rest frame peak energy can accordingly be calculated. Applying this method, we estimate the initial Lorentz factor of the bulk motion and the corresponding rest frame spectral peak energy of GRBs for a new sample where the redshift and the break time in the afterglow are known. Our analysis shows that, in the circumstances, the initial Lorentz factor of the sample would peak at 200 and would be distributed mainly within $(100,400)$, and the peak of the distribution of the corresponding rest frame peak energy would be $0.8keV$ and its main region would be $(0.3keV,3keV)$.

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.

In this paper, we contrast an empirical pulse function of gamma-ray bursts with a theoretical model in order to provide a physical explanation to the parameters of the KRL function and directly relate parameters in both models. We find the ratio of the rise index r to the decay index d, derived when fitted by the KRL function, increases quickly first and then remains nearly invariant with the relative width of local pulses when the width exceeds 2. The rise and decay times of pulses are found to be related to the Lorentz factor by a power law, where the power law index associated with the rise time is less than that of the decay time, and both are close to -2. In addition, the mean asymmetry shows a slight trend of decreasing with Lorentz factors. In plots of decay indices versus asymmetry, there is a descending phase and after this phase there is a rising portion. We find that these long bursts of the KRL sample are mainly associated with those light curves arising from co-moving pulses with the relative width being larger than 0.1. We show that the effect of the co-moving pulse shape on the KRL function parameters of the resulting pulses is considerable and can be distinguished by the decay index d when the relative co-moving pulse width is less than 2.

Qin et al. (2004) have derived a formula of count rates based on a model of highly symmetric expanding fireballs, where the Doppler effect is the key factor to be concerned. In this paper, we employ the formula to both the Qin model and the uniform jet model to study how the rising timescale, $\Delta\tau_{\theta,r}$, and the decay timescale, $\Delta\tau_{\theta,d}$, of a local pulse affect the light curve of GRBs. Our analysis shows that they do make contributions to both rising and decay portions of the light curve of GRBs. Associated with a local pulse with both rising and decay portions, the light curve of GRBs in the rising portion is expected to undergo a concave phase and then a convex one, which we call a "concave-to-convex" character, whereas that in the decay portion is expected to evolve an opposite process which is called a "convex-to-concave" character, regardless of being in Qin model or in the uniform jet model. These characteristics are independent of local pulse shapes and their rest frame radiation forms. The diversity of light curves are caused by different forms of local pulses and different ratios of $\Delta\tau_{\theta,r}$ to $\Delta\tau_{\theta,d}$. We study a sample of 86 GRB pulses detected by the BATSE instrument on board the Compton Gamma Ray Observatory and find that the "concave-to-convex" (expected in the rising portion of the light curve) and "convex-to-concave" (expected in the decay portion) characters, which result from the curvature effect of spherical surface (or the Doppler effect), do exist in the light curve of some GRBs.

As revealed previously, under the assumption that some pulses of gamma-ray bursts are produced by shocks in spherical fireballs or uniform jets of large opening angles, there exists a standard decay form of the profile of pulses arising from very narrow or suddenly dimming local (or intrinsic) pulses due to the relativistic curvature effect (the Doppler effect over the spherical shell surface). Profiles of pulses arising from other local pulses were previously found to possess a reverse S-feature deviation from the standard decay form. We show in this paper that, in addition to the standard decay form shown in Qin et al. (2004), there exists a marginal decay curve associated with a local $\delta$ function pulse with a mono-color radiation. We employ the sample of Kocevski et al. (2003) to check this prediction and find that the phenomenon of the reverse S-feature is common, when compared with both the standard decay form and the marginal decay curve. We accordingly propose to take the marginal decay curve (whose function is simple) as a criteria to check if an observed pulse could be taken as a candidate suffered from the curvature effect. We introduce two quantities $A_1$ and $A_2$ to describe the mentioned deviations within and beyond the $FWHM$ position of the decay phase, respectively. The values of $A_1$ and $A_2$ of pulses of the sample are calculated, and the result suggests that for most of these pulses their corresponding local pulses might contain a long decay time relative to the time scale of the curvature effect.

The FWHM of gamma-ray burst (GRB) pulses is known to be related with energy by a power-law. We wonder if the power-law index $\alpha$ is related with the corresponding local pulse width $FWHM_0$. Seven FRED (fast rise and exponential decay) pulse GRBs are employed to study this issue, where six of them were interpreted recently by the relativistic curvature effect (the Doppler effect of fireballs) and the corresponding local pulses were intensely studied. A regression analysis shows an anti-correlation between $log \alpha$ and $log FWHM_0$ with a slope of $-0.37\pm0.13$. This suggests that, for the class of the GRB pulses which are consequences of the curvature effect, the difference of the local pulse width might lead to the variation of the power law index, where the smaller the width the larger the value of $\alpha$. Since the number of sources employed in this analysis is small, our result is only a preliminary one which needs to be confirmed by larger samples.

We study in details how the pulse width of gamma-ray bursts is related with energy under the assumption that the sources concerned are in the stage of fireballs. Due to the Doppler effect of fireballs, there exists a power law relationship between the two quantities within a limited range of frequency. The power law range and the power law index depend strongly on the observed peak energy $E_p$ as well as the rest frame radiation form, and the upper and lower limits of the power law range can be determined by $E_p$. It is found that, within the same power law range, the ratio of the $FWHM$ of the rising portion to that of the decaying phase of the pulses is also related with energy in the form of power laws. A platform-power-law-platform feature could be observed in the two relationships. In the case of an obvious softening of the rest frame spectrum, the two power law relationships also exist, but the feature would evolve to a peaked one. Predictions on the relationships in the energy range covering both the BATSE and Swift bands for a typical hard burst and a typical soft one are made. A sample of FRED (fast rise and exponential decay) pulse bursts shows that 27 out of the 28 sources belong to either the platform-power-law-platform feature class or the peaked feature group, suggesting that the effect concerned is indeed important for most of the sources of the sample. Among these bursts, many might undergo an obvious softening evolution of the rest frame spectrum.

In this paper, we randomly selected a sample (sample 1) of long gamma-ray bursts, with its size being the same as that of the short burst group (sample 2) (N=500), from the Current BATSE Bursts Catalog. We randomly selected a short burst and assigned its T90 to each of these long bursts. We thus constructed a long burst sample with both the sample size and the distribution of T90 being the same as those of the short burst sample obtained from the Current BATSE Bursts Catalog. Then we calculated the hardness ratio ($hr_{T}$) over the assigned T90 for the long bursts and over their own T90 for the shout bursts, and studied the relation between the hardness ratio and the corresponding T90 for these two samples. We also calculated the hardness ratio ($hr_{t}$) over the randomly selected 64 ms time intervals within the T90, and investigated the relation between this hardness ratio and the selected 64 ms time interval. In addition, the $hr_{t}$ within and beyond the first 2 seconds for all the long bursts (sample 3; N=1541) were also investigated. We found that the KS probabilities of the distributions of the $hr_{T}$ (7.15337E-15) and $hr_{t}$ (9.54833E-10) for samples 1 and 2 are very small, and the average value of $hr_{T}$ and $hr_{t}$ of short bursts are obviously larger than that of the long bursts. The correlations between log$hr_{T}$ and logT90, and between log$hr_{t}$ and log t, for samples 1 and 2 are different. These show that short and long bursts in the first 2 seconds have different characters and they probably originate from different progenitors. For sample 3, for the two time intervals, the KS probability is 5.35828E-5, which suggests that the hardness ratios in different time intervals for long bursts are also different.

Starting from XMM-Newton EPIC-PN data, we present the X-ray variability characteristics of PKS 2155-304 using a simple analysis of the excess variance, \xs, and of the fractional rms variability amplitude, fvar. The scatter in \xs\ and \fvar, calculated using 500 s long segments of the light curves, is smaller than the scatter expected for red noise variability. This alone does not imply that the underlying process responsible for the variability of the source is stationary, since the real changes of the individual variance estimates are possibly smaller than the large scatters expected for a red noise process. In fact the averaged \xs and \fvar, reducing the fluctuations of the individual variances, chang e with time, indicating non-stationary variability. Moreover, both the averaged \sqxs (absolute rms variability amplitude) and \fvar show linear correlation with source flux but in an opposite sense: \sqxs correlates with flux, but \fvar anti-correlates with flux. These correlations suggest that the variability process of the source is strongly non-stationary as random scatters of variances should not yield any correlation. \fvar spectra were constructed to compare variability amplitudes in different energy bands. We found that the fractional rms variability amplitude of the source, when significant variability is observed, increases logarithmically with the photon energy, indicating significant spectral variability. The point-to-point variability amplitude may also track this trend, suggesting that the slopes of the power spectral density of the source are energy-independent. Using the normalized excess variance the black hole mass of \pks was estimated to be about $1.45 \times 10^8 M_{\bigodot}$. This is compared and contrasted with the estimates derived from measurements of the host galaxies.

In the present paper, we investigated the distribution of hardness ratio (HR) for short and long gamma-ray bursts (GRBs) in different time scales for the first two seconds. After including and subtracting the background count, we performed a Kolmogorov--Smirnov (K-S) test to the HR distributions of the two classes of GRBs in each time interval. Our analysis shows that the probabilities of the KS test to the distributions are very small, suggesting that the two classes of bursts are unlikely to arise from the same HR distributions, and indicating that they probably originate from the different physical processes and central engine. In addition, we found that the hardness ratio of short bursts within the time interval of 0$-$0.96 seconds changes hard-to-soft, on the other hand long bursts do not. The two kinds of bursts have different characteristics in the first 2 seconds which might be associated with different physical mechanisms.

In this paper, we employ the peak count rate{\bf \}$C_p$ and the total count $C_{total}$ of light curves to study in the corresponding aspects the relationship between different channel light curves. To make a direct comparison between count rates of different channel light curves we introduce a plot of $C(\tau)$ versus $C_H(\tau)$, where $C(\tau)$ is the count rate of a channel and $C_H(\tau)$ is the count rate of a definite cannel, channel H (see the text). According to the plot we define $\Delta C_{\max}$ as the maximum deviation of the two count rate values of $C(\tau)$ associated with a same count rate value of $% C_H(\tau)$ and define $\Delta S$ as the area confined by the close curve of $C(\tau)$ in the plot to measure the difference of the rising and decaying portions of a light curve relative to the count rate of channel H. Under the assumption that some GRBs observed are in the stage of fireballs which expand relativistically, predictions on the relationships between the four quantities{\bf (}$C_p$, $C_{total}$, $\Delta C_{\max}$, and $\Delta S$) and energy within a wide band, calculated with different rest frame radiation forms and two typical Lorentz factors ($\Gamma =20$ and 200), are made and presented, which would make the test of our model with the coming Swift data easier. Interpretations to the relationships within the mechanism of fireballs are also presented.

In this paper, we derive in a much detail the formula of count rates, in terms of the integral of time, of gamma-ray bursts in the framework of fireballs, where the Doppler effect of the expanding fireball surface is the key factor to be concerned. Effects arising from the limit of the time delay due to the limited regions of the emitting areas in the fireball surface and other factors are investigated. Our analysis shows that the formula of the count rate of fireballs can be expressed as a function of $\tau $ which is the observation time scale relative to the dynamical time scale of the fireball. The profile of light curves of fireballs depends only on the relative time scale, entirely independent of the real time scale and the real size of the objects. It displays in detail how a cutoff tail, or a turn over, feature (called a cutoff tail problem) in the decay phase of a light curve can be formed. This feature is a consequence of a hot spot in the fireball surface, moving towards the observer, and was observed in a few cases previously. By performing fits to the count rate light curves of six sample sources, we show how to obtain some physical parameters from the observed profile of the count rate of GRBs. In addition, the analysis reveals that the Doppler effect of fireballs could lead to a power law relationship between the $FWHM$ of pulses and energy, which were observed previously by many authors.

We make an analysis of the distribution of the extensive Doppler redshift defined as $\widetilde{z}_{Dopp}=(z_{abs}-z_{em})/(1+z_{em})$ with a sample of 1317 absorption redshifts available from 401 quasars. The analysis reveals a bi-peak structure in the distribution, with one component located at $\widetilde{z}_{Dopp}\simeq 0.00$ and the other at $\widetilde{z}% _{Dopp}\simeq -0.01$. A study of some possible causes suggests that the structure can be well interpreted: while the absorbers inside the same cluster as the quasar concerned could contribute to the first component and the less populated space between clusters could explain the gap between the two peaks, the typical distance between clusters could account for the second component. If the bi-peak structure is true, which needs to be confirmed by complete samples, one would be able to obtain some useful cosmological information.

The influence of the Doppler effect in the fireball framework on the spectrum of gamma-ray bursts is investigated. The study shows that the shape of the expected spectrum of an expanding fireball remains almost the same as that of the corresponding rest frame spectrum for constant radiations of the bremsstrahlung, Comptonized, and synchrotron mechanisms as well as for that of the GRB model. The peak flux spectrum and the peak frequency are obviously correlated. When the value of the Lorentz factor becomes 10 times larger, the flux of fireballs would be several orders of magnitude larger. The expansion speed of fireballs is a fundamental factor of the enhancement of the flux of gamma-ray bursts.

We display by numerical calculation how rest frame spectral lines appear in the observed spectrum of gamma-ray bursts due to the Doppler effect in the fireball framework. The analysis shows that: a) in the spectrum of a relativistically expanding fireball, all rest frame lines would shift to higher energy bands and would be significantly smoothed; b) rest frame weak narrow emission lines as well as narrow absorption lines and absorption line forests would be smoothed and would hardly be detectable; c) the features of rest frame broad emission lines as well as both strong and weak broad absorption lines would remain almost unchanged and therefore would be easier to detect; d) deep gaps caused by rest frame broad absorption lines would be significantly filled; e) a rest frame emission line forest would form a single broad line feature; f) the observed relative width of the rest frame very narrow line would approach $ 0.162$; g) when the Lorentz factor $\Gamma $ is large enough, the observed line frequency $\nu_{line}$ and the rest frame line frequency $\nu_{0,line} $ would be related by $\nu_{line}\approx 2\Gamma \nu_{0,line}$. We also investigate the effect of time dependence of the line intensity and the effect of variation of $\Gamma $. We find that the feature of rest frame dimming narrow emission lines would disappear when $\Gamma $ is very large. The form of emission lines would be sharp on both edges when $\Gamma $ varies with time. This phenomenon depends not only on the initial Lorentz factor but also on the observation time.

The statistical shapes of the highest pulse have been studied by aligned method. A wavelet package analysis technique and a developed pulse--finding algorithm have been applied to select the highest pulse from burst profiles observed by BATSE on board CGRO from 1991 April 21 to 1999 January 26. The results of this work show that the statistical shapes of the highest pulses are related to energy: the higher the energy, the narrower the pulse. However, the characteristic structures of the pulses have nothing to do with energy, which strongly supports the previous conclusion that the temporal profiles in different channels are self--similar. The characteristic structures of the pulses can be well described by a model proposed by Norris et al. (1996). The fitting parameters are: $t_r$=0.12, $t_d$=0.16, $\upsilon =1.09$, the ratio of $t_r$ to $t_d$ for the pulse is 0.75. The result leads to our conjecture that the mechanisms of bursts in different gamma-ray bands might be the same. The shock, either an internal or an external one, producing the pulse, might emit photons over the four energy channels in the same way.

A preview study of using observational data of proper motions and Doppler shifted velocities of twin-jets to determine the distance of sources inside and outside our galaxy is made. We investigate the feasibility of this method by studying the uncertainty of the distance caused by the uncertainties of the measured quantities. It shows that, when the motion of components of the jet is relativistic, then the uncertainty of the distance is within the same order of the uncertainties of the measured values of proper motions and Doppler shifted velocities. In particular, when assuming the pattern speed equals the flow speed in the jet, for 10% uncertainties of the measured quantities, the uncertainty of the distance caused by them would be well within 13%. With current technique, this method is realizable. For the convenience of choosing sources to observe, some sources as potential targets are also listed in this paper.

In this paper we show explicit and direct relations between the expected quantities and the observed quantities for the components of Galactic superluminal sources. Basic formulas for calculating the real speed and the angle to the line of sight of the components from the data of proper motions and the distance of the source are presented. We point out that the real speed and the angle to the line of sight of components can be uniquely and directly determined from the observed values of the proper motions and the distance of the source. It is not necessary to calculate an intermediate quantity first and then using the resulted value to calculate the two required quantities. The process of the calculation is simple, and in this way, some extra uncertainties are avoided.

Bidirectional relativistic proper motions of radio components of nearby extragalactic sources give a strong constraint on the determination of the Hubble constant $H_0$. Under the assumption that the real velocity of radio components of extragalactic sources is not less than that of Galactic sources, the value of $H_0$ can be estimated at a high level of accuracy. The assumption is reasonable due to the general belief that the activity in the core of galaxies must be more powerful than that of stars. This method is simple and with only one uncertainty --- the real velocity of components. This uncertainty is related to the value of the real velocity of componenets of Galactic sources and the latter is always well-determined (note that the determination is independent of $H_0$ and the distance of Galactic sources can be directly measured at a rather high level of accuracy). Hopefully the method will play an important role in future research to fix the value of $H_0$. With the data of the three sources available so far and the assumption that the real velocity of componenets of at least one of the sources is not less than a known velocity of componenets of a Galactic source, the constant is estimated to be within $27.08kms^{-1}Mpc^{-1} < H_0 \leq 53.15kms^{-1}Mpc^{-1}$ with this method.

We investigate quantitatively how the peak emission frequency of the overall energy spectrum is at work in distinguishing RBL-like and XBL-like objects. We employ the sample of Giommi et al. (1995) to study the distribution of BL Lacertae objects with various locations of the cutoff of the overall energy spectrum. We find that the sources with the cutoff located at lower frequency are indeed sited in the RBL region of the $\alpha_{ro}-\alpha_{ox}$ plane, while those with the cutoff located at higher frequency are distributed in the XBL region. For a more quantitative study, we employ the BL Lacertae samples presented by Sambruna et al. (1996), where, the peak emission frequency, $\nu _p$, of each source is estimated by fitting the data with a parabolic function. In the plot of $\alpha_{rx}-\log \nu_p$ we find that, in the four different regions divided by the $\alpha_{rx}=0.75$ line and the $\log \nu_p=14.7$ line, all the RBL-like objects are inside the upper left region, while most XBL-like objects are within the lower right region. A few sources are located in the lower left region. No sources are in the upper right region. This result is rather quantitative. It provides an evidence supporting what Giommi et al. (1995) suggested: RBL-like and XBL-like objects can be distinguished by the difference of the peak emission frequency of the overall energy spectrum.

We study correlations between various hardness ratios of gamma-ray bursts (GRBs) and investigate if there are any differences between the two classes of the objects in the distributions of the ratios. The results suggest that, statistically, the slope of the higher part of the spectrum of the long duration bursts has nothing to do with that of the lower part; emissions at higher energy bands from the bursts of both short and long duration classes must be significantly different for different sources, while radiations at lower energy bands of the objects are relatively similar; the spectrum of the short duration bursts must be harder than that of the long duration bursts, confirming what the well-known hardness-duration correlation reveals; the profiles of the spectra between the long duration bursts must be more similar than that between the short duration bursts. The long duration bursts would share more common properties than the short duration bursts. A possible interpretation is proposed with the concept of the Doppler boosting in the relativistic beaming model in AGNs.

The well-known hardness-duration correlation of gamma-ray bursts (GRBs) is investigated with the data of the 4B catalog. We find that, while the hardness ratio and the duration are obviously correlated for the entire set of the 4B catalog, they are not at all correlated for the two subsets divided at the duration of 2 seconds. However, for other subsets with comparable sizes, the two quantities are significantly correlated. The following conclusions are then reached: (1) the existence of two classes of GRBs is confirmed; (2) the hardness ratio and the duration are not at all correlated for any of the two classes; (3) different classes of GRBs have different distributions of the hardness ratio and the duration and it is this difference that causes the correlation between the two quantities for the entire set of the bursts.

In this letter, the old issue of whether redshifts of quasars are of cosmological origin is investigated. We make a plot of absorption redshifts versus emission redshifts for quasars with large amounts of data. Our study shows that, almost all absorption redshifts are smaller than the corresponding emission redshifts. The relation between the absorption and emission redshifts predicted by current cosmological models is well obeyed. The result confirms that redshifts of quasars are indeed distance indicators. It might be the most obvious evidence found so far to be consistent with the cosmological interpretation of quasar redshifts.

In this paper, we will discuss the long-term variations in the optical and the infrared bands of some AGNs. It is interesting to note that the reported periods of AGNs are of the similar values (about 10 years, see Fan et al. 1998a, ApJ, 507 in press, and references therein). DCF method shows the optical and infrared bands are strongly correlated suggesting that the emission mechanisms in the two bands have a common origin.

In this paper, the long-term historical optical (UBVRI) and near-infrared (JHK) data are presented with some new observations in the optical (February 1994-January 1995) and near-infrared (November 1995) bands included for BL Lac object OJ287. The new optical data in V-band are in agreement with the results reported by other authors (Sillanpaa et al. 1996a; Arimoto et al. 1997), a close correlation between the color index of B-V and the magnitude V has been obtained from our new observations; The new infrared observations presented here indicate that the source was at a high level in the infrared band during the observation period; From the available literature, we have got that the largest variations for UBVRIJHK bands are respectively: 4.72mag, 5.93mag, 5.18mag, 4.45mag, 4.07mag, 3.87mag, 3.78mag, and 3.54mag. A strong correlation is found between the optical and near-infrared bands when the DCF method is used, which suggests that these two bands have the same emission mechanism.