Violations of both the weak equivalence principle (WEP) and Lorentz invariance can produce vacuum birefringence, which leads to an energy-dependent rotation of the polarization vector of linearly polarized emission from a given astrophysical source. However, the search for the birefringent effect has been hindered by our ignorance concerning the intrinsic polarization angle in different energy bands. Considering the contributions to the observed linear polarization angle from both the intrinsic polarization angle and the rotation angles induced by violations of the WEP and Lorentz invariance, and assuming the intrinsic polarization angle is an unknown constant, we simultaneously obtain robust bounds on possible deviations from the WEP and Lorentz invariance, by directly fitting the multiwavelength polarimetric data of the optical afterglows of gamma-ray burst (GRB) 020813 and GRB 021004. Here we show that at the $3\sigma$ confidence level, the difference of the parameterized post-Newtonian parameter $\gamma$ values characterizing the departure from the WEP is constrained to be $\Delta\gamma=\left(-4.5^{+10.0}_{-16.0}\right)\times10^{-24}$ and the birefringent parameter $\eta$ quantifying the broken degree of Lorentz invariance is limited to be $\eta=\left(6.5^{+15.0}_{-14.0}\right)\times10^{-7}$. These are the first simultaneous verifications of the WEP and Lorentz invariance in the photon sector. More stringent limits can be expected as the analysis presented here is applied to future multiwavelength polarization observations in the prompt gamma-ray emission of GRB

A nonzero-mass hypothesis for the photon can produces a frequency-dependent dispersion of light, which results in arrival-time differences of photons with different frequencies originating from a given transient source. Extragalactic fast radio bursts (FRBs), with their low frequency emissions, short time durations, and long propagation distances, are excellent astrophysical probes to constrain the rest mass of the photon $m_{\gamma}$. However, the derivation of a limit on $m_{\gamma}$ is complicated by the similar frequency dependences of dispersion expected from the plasma and nonzero photon mass effects. If a handful measurements of redshift for FRBs are available, the different redshift dependences of the plasma and photon mass contributions to the dispersion measure (DM) might be able to break dispersion degeneracy in testing the photon mass. For now, nine FRBs with redshift measurements have been reported, which can turn this idea into reality. Taking into account the DM contributions from both the plasma and a possible photon mass, we use the data on the nine FRBs to derive a combined limit of $m_{\gamma}\leq7.1\times10^{-51}\;{\rm kg}$, or equivalently $m_{\gamma}\leq4.0\times10^{-15}\; {\rm eV}/c^{2}$ at 68\% confidence level, which is essentially as good as or represents a factor of 7 improvement over previous limits obtained by the single FRBs. Additionally, a reasonable estimation for the DM contribution from the host galaxy, $\rm DM_{host}$, can be simultaneously achieved in our analysis. The rapid progress in localizing FRBs will further tighten the constraints on both $m_{\gamma}$ and $\rm DM_{host}$.

GWAC will have been built an integrated FOV of 5,000 $degree^2$ and have already built 1,800 square $degree^2$. The limit magnitude of a 10-second exposure image in the moonless night is 16R. In each observation night, GWAC produces about 0.7TB of raw data, and the data processing pipeline generates millions of single frame alerts. We describe the GWAC Data Processing and Management System (GPMS), including hardware architecture, database, detection-filtering-validation of transient candidates, data archiving, and user interfaces for the check of transient and the monitor of the system. GPMS combines general technology and software in astronomy and computer field, and use some advanced technologies such as deep learning. Practical results show that GPMS can fully meet the scientific data processing requirement of GWAC. It can online accomplish the detection, filtering and validation of millions of transient candidates, and feedback the final results to the astronomer in real-time. During the observation from October of 2018 to December of 2019, we have already found 102 transients.

With the distance sum rule in the Friedmann-Lema\^{\i}tre-Robertson-Walker metric, model-independent constraints on both the Hubble constant $H_0$ and spatial curvature $\Omega_{K}$ can be obtained using strong lensing time-delay data and Type Ia supernova (SN Ia) luminosity distances. This method is limited by the relative low redshifts of SNe Ia, however. Here, we propose using quasars as distance indicators, extending the coverage to encompass the redshift range of strong lensing systems. We provide a novel and improved method of determining $H_0$ and $\Omega_{K}$ simultaneously. By applying this technique to the time-delay measurements of seven strong lensing systems and the known ultraviolet versus X-ray luminosity correlation of quasars, we constrain the possible values of both $H_0$ and $\Omega_{K}$, and find that $H_0=75.3^{+3.0}_{-2.9}$ km $\rm s^{-1}$ $\rm Mpc^{-1}$ and $\Omega_{K}=-0.01^{+0.18}_{-0.17}$. The measured $\Omega_{K}$ is consistent with zero spatial curvature, indicating that there is no significant deviation from a flat universe. If we use flatness as a prior, we infer that $H_0=75.3^{+1.9}_{-1.9}$ km $\rm s^{-1}$ $\rm Mpc^{-1}$, representing a precision of 2.5\%. If we further combine these data with the 1048 current Pantheon SNe Ia, our model-independent constraints can be further improved to $H_0=75.3^{+3.0}_{-2.9}$ km $\rm s^{-1}$ $\rm Mpc^{-1}$ and $\Omega_{K}=0.05^{+0.16}_{-0.14}$. In every case, we find that the Hubble constant measured with this technique is strongly consistent with the value ($\sim 74$ km $\rm s^{-1}$ $\rm Mpc^{-1}$) measured using the local distance ladder, as opposed to the value optimized by {\it Planck}.

Five fast radio bursts (FRBs), including three apparently non-repeating ones FRB 180924, FRB 181112, and FRB 190523, and two repeaters, FRB 121102 and FRB 180916.J0158+65, have already been localized so far. We apply a method developed recently by us (Li et al. 2019) to these five localized FRBs to give a cosmology-insensitive estimate of the fraction of baryon mass in the IGM, $f_{\rm IGM}$. Using the measured dispersion measure (DM) and luminosity distance $d_{\rm L}$ data (inferred from the FRB redshifts and $d_{\rm L}$ of type Ia supernovae at the same redshifts) of the five FRBs, we constrain the local $f_{\rm IGM} = 0.84^{+0.16}_{-0.22}$ with no evidence of redshift dependence. This cosmology-insensitive estimate of $f_{\rm IGM}$ from FRB observations is in excellent agreement with previous constraints using other probes. Moreover, using the three apparently non-repeating FRBs only we get a little looser but consistent result $f_{\rm IGM} = 0.74^{+0.24}_{-0.18}$. In these two cases, reasonable estimations for the host galaxy DM contribution (${\rm DM_{host}}$) can be achieved by modelling it as a function of star formation rate. The constraints on both $f_{\rm IGM}$ and ${\rm DM_{host}}$ are expected to be significantly improved with the rapid progress in localizing FRBs.

The discovery of spectral type transition of active galactic nuclei (AGNs), the so-called changing-look(CL) phenomenon, challenges the widely accepted AGN paradigm, not only in the orientation based Unified Model, but also in the standard disk model. In past decades, only a couple of nearby repeat "changing-look" active galactic nuclei (CL-AGNs) have been identified. Here we report spectroscopic observations of UGC 3223 over the course of 18 years, from 2001 onwards. Combining the spectrum taken in 1987 by Stirpe, we have witnessed its type transitions from $1.5\rightarrow2.0\rightarrow1.8$ over 32 years, and captured a long-lived (at least 10 years) thorough "turn-off" state with a spectrum typical of a Seyfert 2 galaxy. The long-term thorough "turn-off" state probably suggests a once-dormant and an awakening central engine in UGC3223. We argue the (dis)appearance of the broad Balmer emission lines can be explained by the disk-wind BLR model given the evolution of the calculated Eddington ratio of accretion of the supermassive black hole.

Rapid spectroscopy of GRB afterglows is an important and hard task. Based on the archival XS-GRB spectral database, we here perform a pilot study on the afterglow's spectroscopy in the forthcoming \it SVOM \rm era in two aspects. At first, our simulation indicates that the color acquired from the \it SVOM\rm/VT blue and red channels is effective in discriminating between low-$z$ ($z\lesssim3$) and distant GRB candidates until $z\sim6$. Secondly, by doubling the sample size, we find that the previously proposed global photoionization response of the circumburst gas to the prompt emission (i.e., the CIV/CII-$L_{\mathrm{iso}}/E^2_{\mathrm{peak}}$ relationship) is roughly confirmed, although the confirmation is dissatisfactory at the large $L_{\mathrm{iso}}/E^2_{\mathrm{peak}}$ end. We believe that this issue can be further addressed in the \it SVOM \rm era by larger spectroscopy sample, thanks to its capability of rapid identification of optical candidates of afterglow and its anti-solar pointing strategy.

We model the cooling of hybrid neutron stars combining a microscopic nuclear equation of state in the Brueckner-Hartree-Fock approach with different quark models. We then analyze the neutron star cooling curves predicted by the different models and single out the preferred ones. We find that the possibility of neutron p-wave pairing can be excluded in our scenario.

Adaptive optics (AO) is critical in modern astronomy, as well as in optical communications and remote sensing, to deal with the rapid blurring caused by the Earth's turbulent atmosphere. But even the best AO systems are limited by their wavefront sensors, which need to be in an optical plane non-common to the science image, unavoidably leading to errors that limit the reach of current astronomy. They are also insensitive to certain wavefront-error modes, and are poorly suited to injecting light into single-mode optical fibres, important for applications such as high-resolution spectroscopy of extra-solar planets. Here we present a new type of wavefront sensor based on a photonic lantern fibre-mode-converter and deep learning. This new wavefront sensor can be placed at the same focal plane as the science image, and is also optimal for single-mode fibre injection. By only measuring the intensities of an array of single-mode outputs, both phase and amplitude information on the incident wavefront can be reconstructed. We demonstrate the concept with both simulations and an experimental realisation of this novel wavefront sensor, wherein Zernike wavefront errors are recovered from focal-plane measurements to a precision of $2.6\times10^{-5}$ radians mean-squared-error.

With its large field of view, GWAC can record hundreds of meteors every day. These meteors are valuable treasures for some meteor research groups. It is therefore very important to accurately find all of these meteors. To address the challenge of precisely distinguishing meteors from other elongated objects in a GWAC-like sky survey system, we design and implement a meteor candidate recognition algorithm, including the recognizing and morphology analysis of the light curves of the meteor candidates. Although the algorithm may filter out some real meteors, it can provide a sample of meteor with high confidence. After processing the images of Mini-GWAC taken in two months, we detect 109,000 elongated objects in which more than 90 percent of objects are not meteor. Among the elongated objects, about 5.9% objects are identified as meteors with high confidence, after the filters based upon an existence in a single frame, a single peak in the light curves, and a slow variation of the light curves.

The ground-based wide-angle camera array (GWAC) generates millions of single frame alerts per night. After the complicated and elaborate filters by multiple methods, a couple of dozens of candidates are still needed to be confirmed by follow-up observations in real-time. In order to free scientists from the complex and high-intensity follow-up tasks, we developed a Real-time Automatic transient Validation System (RAVS), and introduce here its system architecture, data processing flow, database schema, automatic follow-up control flow, and mobile message notification solution. This system is capable of automatically carrying out all operations in real-time without human intervention, including the validation of transient candidates, the adaptive light-curve sampling for identified targets in multi-band, and the pushing of observation results to the mobile client. The running of RAVS shows that an M-type stellar flare event can be well sampled by RAVS without a significant loss of the details, while the observing time is only less than one-third of the time coverage. Because the control logic of RAVS is designed to be independent of the telescope hardware, RAVS can be conveniently transplanted to other telescopes, especially the follow-up system of SVOM. Some future improvements are presented for the adaptive light-curve sampling, after taking into account both the brightness of sources and the evolution trends of the corresponding light-curves.

The observation of the transient sky through a multitude of astrophysical messengers hasled to several scientific breakthroughs these last two decades thanks to the fast evolution ofthe observational techniques and strategies employed by the astronomers. Now, it requiresto be able to coordinate multi-wavelength and multi-messenger follow-up campaign withinstruments both in space and on ground jointly capable of scanning a large fraction of thesky with a high imaging cadency and duty cycle. In the optical domain, the key challengeof the wide field of view telescopes covering tens to hundreds of square degrees is to dealwith the detection, the identification and the classification of hundreds to thousands of opticaltransient (OT) candidates every night in a reasonable amount of time. In the last decade, newautomated tools based on machine learning approaches have been developed to perform thosetasks with a low computing time and a high classification efficiency. In this paper, we presentan efficient classification method using Convolutional Neural Networks (CNN) to discard anybogus falsely detected in astrophysical images in the optical domain. We designed this toolto improve the performances of the OT detection pipeline of the Ground Wide field AngleCameras (GWAC) telescopes, a network of robotic telescopes aiming at monitoring the opticaltransient sky down to R=16 with a 15 seconds imaging cadency. We applied our trainedCNN classifier on a sample of 1472 GWAC OT candidates detected by the real-time detectionpipeline. It yields a good classification performance with 94% of well classified event and afalse positive rate of 4%.

We present a new model-independent method to determine the spatial curvature and to mitigate the circularity problem affecting the use of quasars as distance indicators. The cosmic-chronometer measurements are used to construct the curvature-dependent luminosity distance $D^{\rm cal}_{L}(\Omega_{K},z)$ using a polynomial fit. Based on the reconstructed $D^{\rm cal}_{L}(\Omega_{K},z)$ and the known ultraviolet versus X-ray luminosity correlation of quasars, we simultaneously place limits on the curvature parameter $\Omega_{K}$ and the parameters characterizing the luminosity correlation function. This model-independent analysis suggests that a mildly closed Universe ($\Omega_{K}=-0.918\pm0.429$) is preferred at the $2.1\sigma$ level. With the calibrated luminosity correlation, we build a new data set consisting of 1598 quasar distance moduli, and use these calibrated measurements to test and compare the standard $\Lambda$CDM model and the $R_{\rm h}=ct$ universe. Both models account for the data very well, though the optimized flat $\Lambda$CDM model has one more free parameter than $R_{\rm h}=ct$, and is penalized more heavily by the Bayes Information Criterion. We find that $R_{\rm h}=ct$ is slightly favoured over $\Lambda$CDM with a likelihood of $\sim57.7\%$ versus 42.3\%.

The Type~Ia supernova (SN~Ia) 2017cfd in IC~0511 (redshift z = 0.01209+- 0.00016$) was discovered by the Lick Observatory Supernova Search 1.6+-0.7 d after the fitted first-light time (FFLT; 15.2 d before B-band maximum brightness). Photometric and spectroscopic follow-up observations show that SN~2017cfd is a typical, normal SN~Ia with a peak luminosity MB ~ -19.2+-0.2 mag, Delta m15(B) = 1.16 mag, and reached a B-band maximum ~16.8 d after the FFLT. We estimate there to be moderately strong host-galaxy extinction (A_V = 0.39 +- 0.03 mag) based on MLCS2k2 fitting. The spectrum reveals a Si~II lambda 6355 velocity of ~11,200 kms at peak brightness. The analysis shows that SN~2017cfd is a very typical, normal SN Ia in nearly every aspect. SN~2017cfd was discovered very young, with multiband data taken starting 2 d after the FFLT, making it a valuable complement to the currently small sample (fewer than a dozen) of SNe~Ia with color data at such early times. We find that its intrinsic early-time (B - V)0 color evolution belongs to the "blue" population rather than to the distinct "red" population. Using the photometry, we constrain the companion star radius to be < 2.5 R_sun, thus ruling out a red-giant companion.

We report our photometric and spectroscopic follow-up observations of the superoutbursts of three dwarf novae (GWAC\,180415A, GWAC\,181017A and GWAC\,181211A) identified independently by the Ground Wide-angle Cameras system, one of the ground-based instruments of the China-France SVOM mission. Based on a combination of our photometry and that taken from the AAVSO, our period analysis of the superhumps enables us to determine the mass ratios to be 0.0967-0.1163, 0.1879-0.1883 and 0.0981-0.1173 for GWAC\,180415A, GWAC\,181017A and GWAC\,181211A, respectively. GWAC\,180415A can be firmly identified as a WZ sge-type dwarf novae due to its long duration ($\sim2$ weeks) multiple rebrightenings with amplitudes of 3-4 magnitudes, the early superhump associated with a double-wave modulation and the low mass ratio. The inferred low mass ratio and location in the $\varepsilon-P_{\mathrm{orb}}$ diagram suggest that GWAC\,181211A is a WZ sge-type dwarf novae candidate. The measured Balmer decrements suggest the Balmer line emission is produced from an optical thick region in GWAC\,180415A and GWAC\,181017A, and from an optical thin region in GWAC\,181211A.

We present the Global Rapid Advanced Network Devoted to the Multi-messenger Addicts (GRANDMA). The network consists of 21 telescopes with both photometric and spectroscopic facilities. They are connected together thanks to a dedicated infrastructure. The network aims at coordinating the observations of large sky position estimates of transient events to enhance their follow-up and reduce the delay between the initial detection and the optical confirmation. The GRANDMA program mainly focuses on follow-up of gravitational-wave alerts to find and characterise the electromagnetic counterpart during the third observational campaign of the Advanced LIGO and Advanced Virgo detectors. But it allows for any follow-up of transient alerts involving neutrinos or gamma-ray bursts, even with poor spatial localisation. We present the different facilities, tools, and methods we developed for this network, and show its efficiency using observations of LIGO/Virgo S190425z, a binary neutron star merger candidate. We furthermore report on all GRANDMA follow-up observations performed during the first six months of the LIGO-Virgo observational campaign, and we derive constraints on the kilonova properties assuming that the events' locations were imaged by our telescopes.

The detection of GW170817 and its electromagnetic counterparts allows to constrain the equation of state of dense matter in new and complementary ways. Very stiff equations of state are ruled out by the upper limit on the average tidal deformability, $\tilde\Lambda \lesssim 800$, imposed by the detected gravitational wave signal. A lower limit, $\tilde\Lambda \gtrsim 400$, can also be extracted by considering the large amount of ejected matter which powers the kilonova AT2017gfo. By using several microscopic nucleonic equations of state, we first confirm the existence of a monotonic relation between $R_{1.5}$ (the radius of the $1.5M_{\odot}$ configuration) and $\tilde\Lambda$. This translates the limits on $\tilde\Lambda$ into limits on the radius: $11.8\,\text{km} \lesssim R_{1.5} \lesssim 13.1\,$km. We then show that the monotonic relation is violated, if a second branch of compact stars composed of quark matter exists, as in the two-families or the twin-stars scenarios. In particular, it is possible to fulfill the limits on $\tilde\Lambda$ while having $R_{1.5}$ significantly smaller than $12\,$km. In both those scenarios the event GW170817/AT2017gfo originates from the merger of a hadronic star and a star containing quark matter.

The second observational campaign of gravitational waves organized by the LIGO/Virgo Collaborations has led to several breakthroughs such as the detection of gravitational wave signals from merger systems involving black holes or neutrons stars. During O2,14 gravitational wave alerts were sent to the astronomical community with sky regions covering mostly over hundreds of square degrees. Among them, 6 have been finally confirmed as real astrophysical events. Since 2013, a new set of ground-based robotic telescopes called GWAC and its pathfinder mini-GWAC have been developed to contribute to the various challenges of themulti-messenger and time domain astronomy. The GWAC system is built up in the framework of the ground-segment system of the SVOM mission that will be devoted to the study of the multi-wavelength transient sky in the next decade. During O2, only the mini-GWAC telescopenetwork was fully operational. Due to the wide field of view and fast automatic follow-up capabilities of the mini-GWAC telescopes, they were well adapted to efficiently cover the sky localization areas of the gravitational wave event candidates. In this paper, we present the mini-GWAC pipeline we have set up to respond to the GW alerts and we report our optical follow-up observations of 8 GW alerts detected during the O2 run. Our observations provided the largest coverage of the GW localization areas in a short latency made by any optical facility. We found tens of optical transient candidates in our images, but none of those could be securely associated with any confirmed black hole-black hole merger event. Based on this first experience and the near future technical improvements of our network system, we will be more competitive to detect the optical counterparts from some gravitational wave events that will be detected during the upcoming O3 run, especially those emerging from binary neutron star mergers.

We here report a spectroscopic identification of two new changing-look AGNs (CL-AGNs): SDSS\,J104705.16+544405.8 and SDSS\,J120447.91+170256.8 both with a "turn-off" type transition from type 1 to type 1.8/1.9. The identification is arrived by a follow-up spectroscopic observation of the five changing-look AGN (CL-AGN) candidates that are extracted from the sample recently released in Macleod et al. The candidates are extract by the authors from the Sloan Digit Sky Survey Data Release 7 spectroscopically confirmed quasars with large amplitude variability. By compiling a sample of 26 previously identified CL-AGNs, we confirm the claim in Macleod et al. that CL-AGNs tend to be biased against low Eddington ratio, and identify an overlap between the CL-AGNs at their dim state and the so-called intermediate-type AGNs. The overlap implies that there two populations of the intermediate-type AGNs with different origins. One is due to the torus orientation effect, and the another the intrinsic change of the accretion rate of the central supermassive blackholes.

The precise measurement of the spectrum of protons, the most abundant component of the cosmic radiation, is necessary to understand the source and acceleration of cosmic rays in the Milky Way. This work reports the measurement of the cosmic ray proton fluxes with kinetic energies from 40 GeV to 100 TeV, with two and a half years of data recorded by the DArk Matter Particle Explorer (DAMPE). This is the first time an experiment directly measures the cosmic ray protons up to ~100 TeV with a high statistics. The measured spectrum confirms the spectral hardening found by previous experiments and reveals a softening at ~13.6 TeV, with the spectral index changing from ~2.60 to ~2.85. Our result suggests the existence of a new spectral feature of cosmic rays at energies lower than the so-called knee, and sheds new light on the origin of Galactic cosmic rays.

It is well-known that ordinary supernovae (SNe) are powered by 56Ni cascade decay. Broad-lined type Ic SNe (SNe Ic-BL) are a subclass of SNe that are not all exclusively powered by 56Ni decay. It was suggested that some SNe Ic-BL are powered by magnetar spin-down. iPTF16asu is a peculiar broad-lined type Ic supernova discovered by the intermediate Palomar Transient Factory. With a rest-frame rise time of only 4 days, iPTF16asu challenges the existing popular models, for example, the radioactive heating (56Ni-only) and the magnetar+56Ni models. Here we show that this rapid rise could be attributed to interaction between the SN ejecta and a pre-existing circumstellar medium ejected by the progenitor during its final stages of evolution, while the late-time light curve can be better explained by energy input from a rapidly spinning magnetar. This model is a natural extension to the previous magnetar model. The mass-loss rate of the progenitor and ejecta mass are consistent with a progenitor that experienced a common envelope evolution in a binary. An alternative model for the early rapid rise of the light curve is the cooling of a shock propagating into an extended envelope of the progenitor. It is difficult at this stage to tell which model (interaction+magnetar+56Ni or cooling+magnetar+56Ni) is better for iPTF16asu. However, it is worth noting that the inferred envelope mass in the cooling+magnetar+56Ni is very high.

The quasi-thermal components found in many Fermi gamma-ray bursts (GRBs) imply that the photosphere emission indeed contributes to the prompt emission of many GRBs. But whether the observed spectra empirically fitted by the Band function or cutoff power law, especially the spectral and peak energy ($E_{p}$) evolutions can be explained by the photosphere emission model alone needs further discussion. In this work, we investigate in detail the time-resolved spectra and $E_{p}$ evolutions of photospheric emission from a structured jet, with an inner-constant and outer-decreasing angular Lorentz factor profile. Also, a continuous wind with a time-dependent wind luminosity has been considered. We show that the photosphere spectrum near the peak luminosity is similar to the cutoff power-law spectrum. The spectrum can have the observed average low-energy spectral index $\alpha $ $ \sim -1$, and the distribution of the low-energy spectral index in our photosphere model is similar to that observed ($-2\lesssim $ $\alpha \lesssim 0$). Furthermore, the two kinds of spectral evolutions during the decay phase, separated by the width of the core ($\theta _{c}$), are consistent with the time-resolved spectral analysis results of several Fermi multi-pulse GRBs and single-pulse GRBs, respectively. Also, for this photosphere model we can reproduce the two kinds of observed $E_{p}$ evolution patterns rather well. Thus, by considering the photospheric emission from a structured jet, we reproduce the observations well for the GRBs best fitted by the cutoff power-law model for the peak-flux spectrum or the time-integrated spectrum.

Extragalactic dispersion measures (DMs) obtained from observations of fast radio bursts (FRBs) are an excellent tool for probing intergalactic medium (IGM) and for conducting cosmography. However, the DM contribution from the IGM (${\rm DM_{IGM}}$) depends on the fraction of baryon mass in the IGM, $f_{\rm IGM}$, which is not properly constrained. As $f_{\rm IGM}(z)$ is geometrically related to the Hubble parameter $H(z)$ and ${\rm DM_{IGM}}(z)$, here we propose that combining two independent measurements of FRBs and $H(z)$ in similar redshift ranges provides a novel and cosmology-free method to constrain the evolution of $f_{\rm IGM}(z)$. Under the assumption that $f_{\rm IGM}$ is evolving with redshift in a functional form, we forecast that the evolution of $f_{\rm IGM}(z)$ can be well inferred in a combined analysis of $\sim3000$ ${\rm DM_{IGM}}(z)$ derived from FRBs and $\sim50$ $H(z)$ derived from the Hubble parameter data. Though the efficiency of our method is not as good as that of the other model-independent method involving the joint measurements of DM and luminosity distance of FRBs, our method offers a new model-independent way to constrain $f_{\rm IGM}(z)$.

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

The detection of the GW170817 neutron star merger event has incited an intense research activity towards the understanding of the nuclear matter equation of state. In this paper we compare in particular the pressure-density relation obtained from heavy-ion collisions with the analysis of the NS merger event. Moreover, we present recent calculations of neutron star's moment of inertia and tidal deformability using various microscopic equations of state for nuclear and hybrid star configurations, and confirm several universal relations. We also discuss the recent constraints for the NS radii determined by GW170817, and find compatible radii between 12 and 13 kilometers, thus identifying the suitable equations of state.

The DArk Matter Particle Explorer (DAMPE), a satellite-based cosmic ray and gamma-ray detector, was launched on December 17, 2015, and began its on-orbit operation on December 24, 2015. In this work we document the on-orbit calibration procedures used by DAMPE and report the calibration results of the Plastic Scintillator strip Detector (PSD), the Silicon-Tungsten tracKer-converter (STK), the BGO imaging calorimeter (BGO), and the Neutron Detector (NUD). The results are obtained using Galactic cosmic rays, bright known GeV gamma-ray sources, and charge injection into the front-end electronics of each sub-detector. The determination of the boundary of the South Atlantic Anomaly (SAA), the measurement of the live time, and the alignments of the detectors are also introduced. The calibration results demonstrate the stability of the detectors in almost two years of the on-orbit operation.

Fast radio bursts (FRBs) are short duration ($\sim$millisecond) radio transients with cosmological origin. The simple sharp features of the FRB signal have been utilized to probe two fundamental laws of physics, namey, testing Einstein's weak equivalence principle and constraining the rest mass of the photon. Recently, \cite{hessels18} found that after correcting for dispersive delay, some of the bursts in FRB 121102 have complex time-frequency structures that include sub-pulses with a time-frequency downward drifting property. Using the delay time between sub-pulses in FRB 121102, here we show that the parameterized post-Newtonian parameter $\gamma$ is the same for photons with different energies to the level of $\left|\gamma_1-\gamma_2\right|<2.5\times 10^{-16}$, which is 1000 times better than previous constraints from FRBs using similar methods. We also obtain a stringent constraint on the photon mass, $m_{\gamma} < 5.1\times10^{-48}$ g, which is 10 times smaller than previous best limits on the photon mass derived through the velocity dispersion method.

Fast radio bursts (FRBs) are intense radio transients whose physical origin remains unknown. Therefore, it is of crucial importance to use a model-independent method to obtain the energy function and cosmic formation rate directly from the observational data. Based on current samples from the Parkes and ASKAP telecsopes, we determine, for the first time, the energy function and formation rate of FRBs by using the Lynden-Bell $\rm C^{-}$ method. The energy function derived from the Parkes sample is a broken power law, however it is a simple power law for the ASKAP sample. For Parkes sample, we derive the formation rate which is roughly consistent with the star formation rate up to $z\sim1.7$, with a local formation rate of $\dot \rho (0) \simeq (3.2\pm 0.3)\; \times {10^4}\;{\rm{Gp}}{{\rm{c}}^{ - 3}}{\rm{y}}{{\rm{r}}^{ - 1}}$ above a detection threshold of $2\,{\rm{Jyms}}$. For ASKAP sample, we find that the formation rate evolves much faster than the star formation rate up to $z\sim0.7$, namely $\dot \rho (z) \propto {(1 + z)^{6.9 \pm 1.9}}$, with a local formation rate of $\dot \rho (0) \simeq (4.6\pm 0.8)\; \times {10^3}\;{\rm{Gp}}{{\rm{c}}^{ - 3}}{\rm{y}}{{\rm{r}}^{ - 1}}$ above a detection threshold of $51\,{\rm{Jyms}}$. This might be a important clue for the physical origin of FRBs.

We perform a detailed analysis and simulations on the transient detection capability in the near-ultraviolet (NUV) band by focusing on some major local transient events. These events include the tidal disruption event due to a supermassive blackhole, the shock breakout of a core-collapse supernova and the flare of a late-type star. Our simulations show that a set of small wide-angle NUV cameras can allow us to detect and study numerous galactic and extra-galactic transient events. Based on the analysis and simulations, here we propose a space-based NUV sky patrol mission by updating the proposal that was originally submitted to the Chinese Space Station mission in 2011. The mission proposed here is composed of a set of eight small wide-field NUV cameras each with a diameter of 20cm. The total sky area simultaneously covered by the NUV cameras is as large as 3000$\mathrm{deg^2}$. The survey cadence ranges from 30 to 300s. The transient events are required to be detected by a dedicated on-board software in real time.

If the weak equivalence principle (WEP) is broken, the measured values of the parametrized post-Newtonian parameter $\gamma$ from photons with left- and right-handed circular polarizations should differ slightly, leading to the arrival-time difference of these two circular components. Thus, the polarization vector of a linearly polarized light may rotate during the propagation. The rotation angle of the polarization vector depends on both the photon energy and the distance of the source. It is believed that if the rotation angle differs by more than $\pi/2$ over an energy range, then the net polarization of the signal would be significantly suppressed and could not be as high as the observed level. Thus, the detection of highly polarized photons implies that the relative rotation angle ($\Delta\Theta$) should not be too large. In this paper, we give a detailed calculation on the evolution of gamma-ray burst (GRB) polarization arising from a possible violation of the WEP, and we find that more than $60\%$ of the initial polarization degree can be conserved even if $\Delta\Theta$ is larger than $\pi/2$. In addition, to tightly constrain the WEP violation, GRBs with harder spectra and polarization observations in a wider energy range seem to be favored. Applying our formulas to the measurements of linear polarization from GRB 110721A and GRB 061122, we obtain strict limits on the differences of the $\gamma$ values as low as $\Delta\gamma<1.3\times10^{-33}$ and $\Delta\gamma<0.8\times10^{-33}$. These provide the most stringent limits to date on a deviation from the WEP, improving at least 6 orders of magnitude over previous bounds.

How Supermassive Blackholes (SMBHs) are spun-up is a key issue of modern astrophysics. As an extension of the study in Wang et al. (2016), we here address the issue by comparing the host galaxy properties of nearby ($z<0.05$) radio-selected Seyfert 2 galaxies. With the two-dimensional bulge+disk decompositions for the SDSS $r$-band images, we identify a dichotomy on various host galaxy properties for the radio-powerful SMBHs. By assuming the radio emission from the jet reflects a high SMBH spin, which stems from the well-known BZ mechanism of jet production, high-mass SMBHs (i.e., $M_{\mathrm{BH}}>10^{7.9}M_\odot$) have a preference for being spun-up in classical bulges, and low-mass SMBHs (i.e., $M_{\mathrm{BH}}=10^{6-7}M_\odot$) in pseudo-bulges. This dichotomy suggests and confirms that high-mass and low-mass SMBHs are spun-up in different ways, i.e., a major "dry" merger and a secular evolution.

Lorentz invariance violations (LIV) can yield vacuum birefringence, which results in an energy-dependent rotation of the polarization vector of linearly polarized emission from astrophysical sources. It is believed that if the relative rotation angle ($\Delta\Theta$) of the polarization vector of high energy photons with respect to that of low energy photons is larger than $\pi/2$, then the net polarization of the signal would be significantly depleted and could not be as high as the observed level. Hence, the measurement of high polarization implies that $\Delta\Theta$ should not be too large. In this work, we assemble recent detections of prompt emission polarization in gamma-ray bursts (GRBs), all of whom have high detection significance. Following the method shown in Lin et al. (2016), we give a detailed calculation on the polarization evolution arising from the LIV effect for each GRB, and confirm that, even if $\Delta\Theta$ is approaching to $\pi/2$, the net polarization is not severely suppressed, and more than 60\% of the initial polarization can be conserved. Applying this method to our GRB polarimetric data, we improve existing sensitivities to LIV involving photons by factors ranging from two to ten. In addition, we prove that our constraints are not greatly affected by uncertainties in the spectral parameters of GRBs.

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

The origin of fast radio bursts (FRBs) is still unknown. Multi-wavelength observations during or shortly after the FRB phase would be essential to identify the counterpart of an FRB and to constrain its progenitor and environment. In this work, we investigate the brightness of the "fast optical bursts" (FOBs) associated with FRBs and the prospects of detecting them. We investigate several inverse Compton (IC) scattering processes that might produce an FOB, including both the one-zone and two-zone models. We also investigate the extension of the same mechanism of FRB emission to the optical band. We find that a detectable FOB with the current and forthcoming telescopes is possible under the IC scenarios with very special conditions. In particular, the FRB environment would need to invoke a neutron star with an extremely strong magnetic field and an extremely fast spin, or an extremely young supernova remnant (SNR) surrounding the FRB source. Furthermore, most electrons in the source are also required to have a fine-tuned energy distribution such that most of the IC energy is channeled in the optical band. We conclude that the prospect of detecting FOBs associated with FRBs is low. On the other hand, if FOBs are detected from a small fraction of FRBs, these FOBs would reveal extreme physical conditions in the FRB environments.

The excessive dispersion measure (DM) of fast radio bursts (FRBs) has been proposed to be a powerful tool to study intergalactic medium (IGM) and to perform cosmography. One issue is that the fraction of baryons in the IGM, $f_{\rm IGM}$, is not properly constrained. Here we propose a method of estimating $f_{\rm IGM}$ using a putative sample of FRBs with the measurements of both DM and luminosity distance $d_{\rm L}$. The latter can be obtained if the FRB is associated with a distance indicator (e.g. a gamma-ray burst or a gravitational wave event), or the redshift $z$ of the FRB is measured and $d_{\rm L}$ at the corresponding $z$ is available from other distance indicators (e.g. type Ia supernovae) at the same redshift. Since $d_{\rm L}/{\rm DM}$ essentially does not depend on cosmological parameters, our method can determine $f_{\rm IGM}$ independent of cosmological parameters. We parameterize $f_{\rm IGM}$ as a function of redshift and model the DM contribution from a host galaxy as a function of star formation rate. Assuming $f_{\rm IGM}$ has a mild evolution with redshift with a functional form and by means of Monte Carlo simulations, we show that an unbiased and cosmology-independent estimate of the present value of $f_{\rm IGM}$ with a $\sim 12\%$ uncertainty can be obtained with 50 joint measurements of $d_{\rm L}$ and DM. In addition, such a method can also lead to a measurement of the mean value of DM contributed from the local host galaxy.

The observation of Type Ia supernovae (SNe Ia) plays an essential role in probing the expansion history of the universe. But the possible presence of cosmic opacity can degrade the quality of SNe Ia. The gravitational-wave (GW) standard sirens, produced by the coalescence of double neutron stars and black hole--neutron star binaries, provide an independent way to measure the distances of GW sources, which are not affected by cosmic opacity. In this paper, we first propose that combining the GW observations of third-generation GW detectors with SN Ia data in similar redshift ranges offers a novel and model-independent method to constrain cosmic opacity. Through Monte Carlo simulations, we find that one can constrain the cosmic opacity parameter $\kappa$ with an accuracy of $\sigma_{\kappa}\sim0.046$ by comparing the distances from 100 simulated GW events and 1048 current Pantheon SNe Ia. The uncertainty of $\kappa$ can be further reduced to $\sim0.026$ if 800 GW events are considered. We also demonstrate that combining 2000 simulated SNe Ia and 1000 simulated GW events could result in much severer constraints on the transparent universe, for which $\kappa=0.0000\pm0.0044$. Compared to previous opacity constraints involving distances from other cosmic probes, our method using GW standard sirens and SN Ia standard candles at least achieves competitive results.

The detection of the high-energy ($\sim290$ TeV) neutrino coincident with the flaring blazar TXS 0506+056, the first and only $3\sigma$ neutrino-source association to date, provides new, multimessenger tests of the weak equivalence principle (WEP) and Lorentz invariance. Assuming that the flight time difference between the TeV neutrino and gamma-ray photons from the blazar flare is mainly caused by the gravitational potential of the Laniakea supercluster of galaxies, we show that the deviation from the WEP for neutrinos and photons is conservatively constrained to have an accuracy of $10^{-6}-10^{-7}$, which is 3--4 orders of magnitude better than previous results placed by MeV neutrinos from supernova 1987A. In addition, we demonstrate that the association of the TeV neutrino with the blazar flare sets limits on the energy scales of quantum gravity for both linear and quadratic violations of Lorentz invariance (LIV) to $E_{\rm QG, 1}>3.2\times10^{15}-3.7\times10^{16}$ GeV and $E_{\rm QG, 2}>4.0\times10^{10}-1.4\times10^{11}$ GeV. These improve previous limits on both linear and quadratic LIV energy scales in neutrino propagation by 5--7 orders of magnitude.

We model neutron star cooling with several microscopic nuclear equations of state based on different nucleon-nucleon interactions and three-body forces, and compatible with the recent GW170817 neutron star merger event. They all feature strong direct Urca processes. We find that all models are able to describe well the current set of cooling data for isolated neutron stars, provided that large and extended proton 1S0 gaps and no neutron 3PF2 gaps are active in the stellar matter. We then analyze the neutron star mass distributions predicted by the different models and single out the preferred ones.

We calculate neutron star's moment of inertia and deformabilities using various microscopic equations of state for nuclear and hybrid star configurations. Correlations between the various observables are examined and we confirm several universal relations. We focus in particular on the constraints for the neutron star radii imposed by a determination of the average tidal deformability of the binary neutron star system GW170817. We find compatible radii between 12 and 13 kilometers and identify the suitable equations of state.

We here report a spectroscopic monitor for the supernova SN\,2017iuk associated with the long-duration low-luminosity gamma-ray burst GRB\,171205A at a redshift of 0.037, which is up to now the third GRB-SN event away from us. Our spectroscopic observations and spectral analysis allow us to identify SN\,2017iuk as a typical broad-line type Ic SN. A comparison study suggests that the type-IcBL SN\,2017iuk resembles to SN\,2006aj in following aspects: 1) similar spectra at the nearby epochs, 2) comparable evolution of the photospheric velocity obtained from the measurements based on both \ion{Si}{2}$\lambda$6355 line and spectral modeling, and 3) comparable explosion parameters. This analogy could imply a formation of a neutron star in the core-collapse of GRB\,171205A/SN\,2017iuk as previously suggested in GRB\,060218/SN\,2006aj. The properties of the host galaxy is discussed, which suggests that GRB\,171205A/SN\,2017iuk occurred in an early type (S0), high-mass, starforming galaxy with low specific SFR and solar metallicity.

The detection of gravitational waves (GWs) provides a direct way to measure the luminosity distance, which enables us to probe cosmology. In this paper, we continue to expand the application of GW standard sirens in cosmology, and propose that the spatial curvature can be estimated in a model-independent way by comparing the distances from future GW sources and current cosmic-chronometer observations. We expect an electromagnetic counterpart of the GW event to give the source redshift, and simulate hundreds of GW data from the coalescence of double neutron stars and black hole--neutron star binaries using the Einstein Telescope as reference. Our simulations show that, from 100 simulated GW events and 31 current cosmic-chronometer measurements, the error of the curvature parameter $\Omega_{K}$ is expected to be constrained at the level of $\sim0.125$. If 1000 GW events are observed, the uncertainty of $\Omega_{K}$ would be further reduced to $\sim0.040$. We also find that adding 50 mock $H(z)$ data (consisting of 81 cosmic-chronometer data and 1000 simulated GW events) could result in much tighter constraint on the zero cosmic curvature, for which, $\Omega_{K}=-0.002\pm0.028$. Compared to some actual model-independent curvature tests involving the distances from other cosmic probes, this method with GW data achieves constraints with much higher precision.

We examine whether a comparison between wCDM and R_h=ct using merged Type Ia SN catalogs produces results consistent with those based on a single homogeneous sample. Using the Betoule et al. (2014) joint lightcurve analysis (JLA) of a combined sample of 613 events from SNLS and SDSS-II, we estimate the parameters of the two models and compare them. We find that the improved statistics can alter the model selection in some cases, but not others. In addition, based on the model fits, we find that there appears to be a lingering systematic offset of ~0.04-0.08 mag between the SNLS and SDSS-II sources, in spite of the cross-calibration in the JLA. Treating wCDM, LCDM and R_h=ct as separate models, we find in an unbiased pairwise statistical comparison that the Bayes Information Criterion (BIC) favors the R_h=ct universe with a likelihood of 82.8% versus 17.2% for wCDM, but the ratio of likelihoods is reversed (16.2% versus 83.8%) when w_de=-1 (i.e., LCDM) and strongly reversed (1.0% versus 99.0%) if in addition k=0 (i.e., flat LCDM). We point out, however, that the value of k is a measure of the net energy (kinetic plus gravitational) in the Universe and is not constrained theoretically, though some models of inflation would drive k to 0 due to an expansion-enforced dilution. Since we here consider only the basic LCDM model, the value of k needs to be measured and, therefore, the pre-assumption of flatness introduces a significant bias into the BIC.

The photon zero-mass hypothesis has been investigated for a long time using the frequency-dependent time delays of radio emissions from astrophysical sources. However, the search for a rest mass of the photon has been hindered by the similarity between the frequency-dependent dispersions due to the plasma and nonzero photon mass effects. Considering the contributions to the observed dispersion measure from both the plasma and nonzero photon mass effects, and assuming the dispersion induced by the plasma effect is an unknown constant, we obtain a robust limit on the photon mass by directly fitting a combination of the dispersion measures of radio sources. Using the observed dispersion measures from two statistical samples of extragalactic pulsars, here we show that at the 68\% confidence level, the constraints on the photon mass can be as low as $m_{\gamma}\leq1.51\times10^{-48}~\rm kg\simeq8.47\times10^{-13}~{\rm eV}/c^{2}$ for the sample of 22 radio pulsars in the Large Magellanic Cloud and $m_{\gamma}\leq1.58\times10^{-48}~\rm kg\simeq8.86\times10^{-13} {\rm eV}/c^{2}$ for the other sample of 5 radio pulsars in the Small Magellanic Cloud, which are comparable with that obtained by a single extragalactic pulsar. Furthermore, the statistical approach presented here can also be used when more fast radio bursts with known redshifts are detected in the future.

The first gravitational-wave event from the merger of a binary neutron star system (GW170817) was detected recently. The associated short gamma-ray burst (GRB 170817A) has a low isotropic luminosity ($\sim 10^{47}$ erg s$ ^{-1}$) and a peak energy $E_{p}\sim 145$ keV during the initial main emission between -0.3 and 0.4 s. The origin of this short GRB is still under debate, but a plausible interpretation is that it is due to the off-axis emission from a structured jet. We consider two possibilities. First, since the best-fit spectral model for the main pulse of GRB 170817A is a cutoff power law with a hard low-energy photon index ($\alpha =-0.62_{-0.54}^{+0.49} $), we consider an off-axis photosphere model. We develop a theory of photosphere emission in a structured jet and find that such a model can reproduce a low-energy photon index that is softer than a blackbody through enhancing high-latitude emission. The model can naturally account for the observed spectrum. The best-fit Lorentz factor along the line of sight is $\sim 20$, which demands that there is a significant delay between the merger and jet launching. Alternatively, we consider that the emission is produced via synchrotron radiation in an optically thin region in an expanding jet with decreasing magnetic fields. This model does not require a delay of jet launching but demands a larger bulk Lorentz factor along the line of sight. We perform Markov Chain Monte Carlo fitting to the data within the framework of both models and obtain good fitting results in both cases.

Recently, some theoretical models predicted that a small fraction of fast radio bursts (FRBs) could be associated with gravitational waves (GWs). In this work, we discuss the possibility of using GW/FRB association systems, if commonly detected in the future, as a complementary cosmic probe. We propose that upgraded standard sirens can be constructed from the joint measurements of luminosity distances $D_{L}$ derived from GWs and dispersion measures ${\rm DM_{IGM}}$ derived from FRBs (i.e., the combination $D_{L}\cdot{\rm DM_{IGM}}$). Moreover, unlike the traditional standard-siren approach (i.e., the $D_{L}$ method) and the $\rm DM_{IGM}$ method that rely on the optimization of the Hubble constant $H_0$, this $D_{L}\cdot{\rm DM_{IGM}}$ method has the advantage of being independent of $H_0$. Through Monte Carlo simulations, we prove that the $D_{L}\cdot{\rm DM_{IGM}}$ method is more effective to constrain cosmological parameters than $D_{L}$ or ${\rm DM_{IGM}}$ separately, and enables us to achieve accurate multimessenger cosmology from around 100 GW/FRB systems. Additionally, even if GW/FRB associations are not exist, the methodology developed here can still be applied to those GWs and FRBs that occur at the same redshifts.

Exoplanets, or planets outside our own solar system, have long been of interest to astronomers; however, only in the past two decades have scientists had the technology to characterize and study planets so far away from us. With advanced telescopes and spectrometers, astronomers are discovering more exoplanets every year. The two most prolific methods for detecting exoplanets, radial velocity, which measures the wobble of an exoplanet's parent star, and transit, which detects the passage of a planet in front of its star, are responsible for the discovery of the vast majority of exoplanets we know about so far. Other methods such as direct imaging, timing, and gravitational microlensing are less applicable but can sometimes yield accurate data to confirm planets found by radial velocity or transit photometry. Astronomers will continue to find large numbers of new exoplanets in the near future, but it will become increasingly harder to do so after all the planets close to our solar system are discovered.

We present our optical photometric and spectroscopical observations of GRB 140629A. A redshift of $z=2.275\pm0.043$ is measured through the metal absorption lines in our spectroscopic data. Using our photometric data and multiple observational data observed with other telescopes, we show that its optical light curve is well interpreted with the standard forward shock models in the thin shell case. Its optical-X-ray afterglow spectrum is jointly fitted with a single power-law function, yielding a photon index of $-1.90\pm 0.05$. The optical extinction and neutral hydrogen absorption of the GRB host galaxy are negligible. The fit to the light curve with the standard models shows that the ambient density is $60\pm 9$ cm$^{-3}$ and the GRB radiating efficiency is as low as $\sim 0.24\%$, likely indicating a baryonic-dominated ejecta of this GRB. This burst agrees well with the $L_{\rm p, iso}-E_p^{'}-\Gamma_0$ relation, but confidently violates those empirical relations involving geometric corrections (or jet break time). This gives rise to an issue of possible selection effect on these relations since the jet opening angle of this GRB is extremely narrow (0.04 rad).

We here report an identification of SDSSJ141324+530527.0 (SBS1411+533) at $z=0.456344$ as a new "changing-look" quasar with a "turn-on" spectral type transition from Type-1.9/2 to Type-1 within a rest frame time scale of 1-10 yr by a comparison of our new spectroscopic observation and the Sloan Digital Sky Survey (SDSS) archive data base. The SDSS DR7 spectrum taken in 2003 is dominated by a starlight emission from host galaxies redward of the Balmer limit, and has non-detectable broad H$\beta$ line. The new spectrum taken by us on June 1st, 2017 and SDSS DR14 spectrum taken on May 29, 2017 indicate that the object is of a typical quasar spectrum with a blue continuum and strong Balmer broad emission lines. In addition, an intermediate spectral type can be identified in the SDSS DR13 spectrum taken in 2015. The invariability of the line wing of MgII$\lambda2800$ emission and time scale argument (The invariability of [OIII]$\lambda$5007 line blue asymmetry) suggests that a variation of obscuration (an accelerating outflow) is not a favorable scenario. The time scale argument allows us to believe the type transition is possibly caused by either a viscous radial inflow or a disk instability around a $\sim5-9\times10^{7}M_\odot$ black hole.

By using the line ratio \ion{C}{4}$\lambda1549$/\ion{C}{2}$\lambda1335$ as a tracer of ionization ratio of the interstellar medium (ISM) illuminated by a long gamma-ray burst (LGRB), we identify a global photoionization response of the ionization ratio to the photon luminosity of the prompt emission assessed by either $L_{\mathrm{iso}}/E_{\mathrm{peak}}$ or $L_{\mathrm{iso}}/E^2_{\mathrm{peak}}$. The ionization ratio increases with both $L_{\mathrm{iso}}/E_{\mathrm{peak}}$ and $L_{\mathrm{iso}}/E^2_{\mathrm{peak}}$ for a majority of the LGRBs in our sample, although there are a few outliers. The identified dependence of \ion{C}{4}/\ion{C}{2} on $L_{\mathrm{iso}}/E^2_{\mathrm{peak}}$ suggests that the scatter of the widely accepted Amati relation is related with the ionization ratio in ISM. The outliers tend to have relatively high \ion{C}{4}/\ion{C}{2} values as well as relatively high \ion{C}{4}$\lambda1549$/\ion{Si}{4}$\lambda1403$ ratios, which suggests an existence of Wolf-Rayet stars in the environment of these LGRBs. We finally argue that the outliers and the LGRBs following the identified \ion{C}{4}/\ion{C}{2}$-L_{\mathrm{iso}}/E_{\mathrm{peak}}$ ($L_{\mathrm{iso}}/E^2_{\mathrm{peak}}$) correlation might come from different progenitors with different local environments.

Broad-lined type Ic supernovae (SNe Ic-BL) are a subclass of rare core collapse SNe whose energy source is debated in the literature. Recently a series of investigations on SNe Ic-BL with the magnetar (plus 56Ni) model were carried out. Evidence for magnetar formation was found for the well-observed SNe Ic-BL 1998bw and 2002ap. In this paper we systematically study a large sample of SNe Ic-BL not associated with gamma-ray bursts. We use photospheric velocity data determined in a homogeneous way. We find that the magnetar+56Ni model provides a good description of the light curves and velocity evolution of our sample of SNe Ic-BL, although some SNe (not all) can also be described by the pure-magnetar model or by the two-component pure-56Ni model (3 out of 12 are unlikely explained by two-component model). In the magnetar+56Ni model, the amount of 56Ni required to explain their luminosity is significantly reduced, and the derived initial explosion energy is, in general, in accordance with neutrino heating. Some correlations between different physical parameters are evaluated and their implications regarding magnetic field amplification and the total energy reservoir are discussed.

High energy cosmic ray electrons plus positrons (CREs), which lose energy quickly during their propagation, provide an ideal probe of Galactic high-energy processes and may enable the observation of phenomena such as dark-matter particle annihilation or decay. The CRE spectrum has been directly measured up to $\sim 2$ TeV in previous balloon- or space-borne experiments, and indirectly up to $\sim 5$ TeV by ground-based Cherenkov $\gamma$-ray telescope arrays. Evidence for a spectral break in the TeV energy range has been provided by indirect measurements of H.E.S.S., although the results were qualified by sizeable systematic uncertainties. Here we report a direct measurement of CREs in the energy range $25~{\rm GeV}-4.6~{\rm TeV}$ by the DArk Matter Particle Explorer (DAMPE) with unprecedentedly high energy resolution and low background. The majority of the spectrum can be properly fitted by a smoothly broken power-law model rather than a single power-law model. The direct detection of a spectral break at $E \sim0.9$ TeV confirms the evidence found by H.E.S.S., clarifies the behavior of the CRE spectrum at energies above 1 TeV and sheds light on the physical origin of the sub-TeV CREs.

Lorentz invariance violation (LIV) can manifest itself by an energy-dependent vacuum dispersion of light, which leads to arrival-time differences of photons with different energies originating from the same astronomical source. The spectral lags of gamma-ray bursts (GRBs) have been widely used to investigate the possible LIV effect. However, all current investigations used lags extracted in the observer frame only. In this work, we present, for the first time, an analysis of the LIV effect and its redshift dependence in the cosmological rest frame. Using a sample of 56 GRBs with known redshifts, we obtain a robust limit on LIV by fitting their rest-frame spectral lag data both using a maximization of the likelihood function and a minimum $\chi^{2}$ statistic. Our analysis indicates that there is no evidence of LIV. Additionally, we test the LIV in different redshift ranges by dividing the full sample into four redshift bins. We also find no evidence for the redshift variation of the LIV effect.

We perform a systematic study of outflow in the narrow-line region (NLR) of active galactic nuclei (AGNs) at $z\sim0.4-0.8$ basing upon a large sample of $\sim900$ quasars at $z\sim 0.4-0.8$. The sample is extracted from the Sloan Digital Sky Survey by mainly requiring 1) the g-band magnitude is brighter than 19 magnitude; and 2) the [OIII]$\lambda5007$ emission line has a signal-to-noise ration larger than 30. Profiles of multiple emission lines are modeled by a sum of several Gaussian functions. The spectral analysis allows us to identify 1) a prevalence of both [OIII]$\lambda5007$ line blue asymmetry and bulk velocity blueshift of both [NeIII]$\lambda3869$ and [NeV]$\lambda3426$ lines, when the [\ion{O}{2}]$\lambda3727$ line is used as a reference. The velocity offset of [\ion{O}{3}]$\lambda5007$ line is, however, distributed around zero value, except for a few outliers. 2) not only the significant [OIII]$\lambda5007$ line asymmetry, but also the large bulk velocity offsets of [NeIII]$\lambda3869$ and [NeV]$\lambda3426$ emission lines tend to occur in the objects with high $L/L_{\mathrm{Edd}}$, which is considerably consistent with the conclusions based on local AGNs. With three $M_{\mathrm{BH}}$ estimation methods, the significance level of the trend is found to be better than $2.9\sigma$, $3.2\sigma$ and $1.8\sigma$ for [OIII], [NeIII] and [NeV], respectively. \rm After excluding the role of radio jets, the revealed dependence of NLR gas outflow on $L/L_{\mathrm{Edd}}$ allows us to argue that the pressure caused by the wind/radiation launched/emitted from central supermassive black hole is the most likely origin of the outflow in these distant quasars, which implies that the outflow in luminous AGNs up to $z\sim1$ have the same origin.

The coincident detection of a gravitational-wave (GW) event GW170817 with electromagnetic (EM) signals (e.g., a short gamma-ray burst SGRB 170817A or a macronova) from a binary neutron star merger within the nearby galaxy NGC 4933 provides a new, multimessenger test of the weak equivalence principle (WEP), extending the WEP test with GWs and photons. Assuming that the arrival time delay between the GW signals from GW170817 and the photons from SGRB 170817A or the macronova is mainly attributed to the gravitational potential of the Milky Way, we demonstrate that the strict upper limits on the deviation from the WEP are $\Delta \gamma<1.4\times10^{-3}$ for GW170817/macronova and $\Delta \gamma <5.9\times10^{-8}$ for GW170817/SGRB 170817A. A much more severe constraint on the WEP accuracy can be achieved ($\sim0.9\times10^{-10}$) for GW170817/SGRB 170817A when we consider the gravitational potential of the Virgo Cluster, rather than the Milky Way's gravity. This provides the tightest limit to date on the WEP through the relative differential variations of the $\gamma$ parameter for two different species of particles. Compared with other multimessenger (photons and neutrinos) results, our limit is 7 orders of magnitude tighter than that placed by the neutrinos and photons from supernova 1987A, and is almost as good as or is an improvement of 6 orders of magnitude over the limits obtained by the low-significance neutrinos correlated with GRBs and a blazar flare.

Assuming that the early optical emission is dominated by the external reverse shock (RS) in the standard model of gamma-ray bursts (GRBs), we intend to constrain RS models with the initial Lorentz factor $\Gamma_0$ of the outflows based on the ROTSE-III observations. We consider two cases of the RS behavior: the relativistic shock and the non-relativistic shock. For homogeneous interstellar medium (ISM) and wind circum-burst environment, the constraints can be achieved by the fact that the peak flux $F_{\rm \nu}$ at the RS crossing time should be lower than the observed upper limit $F_{\rm \nu,limit}$. We consider the different spectral regimes that the observed optical frequency $\nu_{\rm opt}$ may locate in, which are divided by the orders for the minimum synchrotron frequency $\nu_{\rm m}$ and the cooling frequency $\nu_{\rm c}$. {\bf Considering the homogeneous and wind environment around GRBs, we find that the relativistic RS case can be constrained by the (upper and lower) limits of $\Gamma_0$ in a large range from about hundreds to thousands for 36 GRBs reported by ROTSE-III. The constraints on the non-relativistic RS case are achieved with limits of $\Gamma_0$ for 26 bursts ranging from $\sim 30$ to $\sim 350$.} The lower limits of $\Gamma_0$ achieved for the relativistic RS model is disfavored based on the previously discovered correlation between the initial Lorentz factor $\Gamma_0$ and the isotropic gamma-ray energy $E_{\rm \gamma, iso}$ released in prompt phase.

We here report an identification of SDSS J090152.04+624342.6 as a new "overlapping-trough" iron low-ionization broad absorption line quasar at redshift of $z\sim2.1$. No strong variation of the broad absorption lines can be revealed through the two spectra taken by the Sloan Digital Sky Survey with a time interval of $\sim6$yr. Further optical and infrared spectroscopic study on this object is suggested.

In this paper, we discuss the possibility of using strongly lensed gravitational waves (GWs) and their electromagnetic (EM) counterparts as powerful cosmic rulers. In the EM domain, it has been suggested that joint observations of the time delay ($\Delta\tau$) between lensed quasar images and the velocity dispersion ($\sigma$) of the lensing galaxy (i.e., the combination $\Delta\tau/\sigma^{2}$) are able to constrain the cosmological parameters more strongly than $\Delta\tau$ or $\sigma^{2}$ separately. Here, for the first time, we propose that this $\Delta\tau/\sigma^{2}$ method can be applied to the strongly lensed systems observed in both GW and EM windows. Combining the redshifts, images and $\sigma$ observed in the EM domain with the very precise $\Delta\tau$ derived from lensed GW signals, we expect that accurate multimessenger cosmology can be achieved in the era of third-generation GW detectors. Comparing with the constraints from the $\Delta\tau$ method, we prove that using $\Delta\tau/\sigma^{2}$ can improve the discrimination between cosmological models. Furthermore, we demonstrate that with $\sim50$ strongly lensed GW-EM systems, we can reach a constraint on the dark energy equation of state $w$ comparable to the 580 Union2.1 Type Ia supernovae data. Much more stringent constraints on $w$ can be obtained when combining the $\Delta\tau$ and $\Delta\tau/\sigma^{2}$ methods.

In this paper, we propose to estimate the spatial curvature of the universe and the cosmic opacity in a model-independent way with expansion rate measurements, $H(z)$, and type Ia supernova (SNe Ia). On the one hand, using a nonparametric smoothing method Gaussian process, we reconstruct a function $H(z)$ from opacity-free expansion rate measurements. Then, we integrate the $H(z)$ to obtain distance modulus $\mu_{\rm H}$, which is dependent on the cosmic curvature. On the other hand, distances of SNe Ia can be determined by their photometric observations and thus are opacity-dependent. In our analysis, by confronting distance moduli $\mu_{\rm H}$ with those obtained from SNe Ia, we achieve estimations for both the spatial curvature and the cosmic opacity without any assumptions for the cosmological model. Here, it should be noted that light curve fitting parameters, accounting for the distance estimation of SNe Ia, are determined in a global fit together with the cosmic opacity and spatial curvature to get rid of the dependence of these parameters on cosmology. In addition, we also investigate whether the inclusion of different priors for the present expansion rate ($H_0$: global estimation, $67.74\pm 0.46~\rm km~ s^{-1} ~Mpc^{-1}$, and local measurement, $73.24\pm 1.74~\rm km~ s^{-1} ~Mpc^{-1}$) exert influence on the reconstructed $H(z)$ and the following estimations of the spatial curvature and cosmic opacity. Results show that, in general, a spatially flat and transparent universe is preferred by the observations. Moreover, it is suggested that priors for $H_0$ matter a lot. Finally, we find that there is a strong degeneracy between the curvature and the opacity.

The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives include the study of galactic cosmic rays up to $\sim 10$ TeV and hundreds of TeV for electrons/gammas and nuclei respectively, and the search for dark matter signatures in their spectra. In this paper we illustrate the layout of the DAMPE instrument, and discuss the results of beam tests and calibrations performed on ground. Finally we present the expected performance in space and give an overview of the mission key scientific goals.

Violations of Lorentz invariance can lead to an energy-dependent vacuum dispersion of light, which results in arrival-time differences of photons arising with different energies from a given transient source. In this work, direction-dependent dispersion constraints are obtained on nonbirefringent Lorentz-violating effects, using the observed spectral lags of the gamma-ray burst GRB 160625B. This burst has unusually large high-energy photon statistics, so we can obtain constraints from the true spectral time lags of bunches of high-energy photons rather than from the rough time lag of a single highest-energy photon. Also, GRB 160625B is the only burst to date having a well-defined transition from positive lags to negative lags, which provides a unique opportunity to distinguish Lorentz-violating effects from any source-intrinsic time lag in the emission of photons of different energy bands. Our results place comparatively robust two-sided constraints on a variety of isotropic and anisotropic coefficients for Lorentz violation, including first bounds on Lorentz-violating effects from operators of mass dimension ten in the photon sector.

Einstein's weak equivalence principle (WEP) states that any freely falling, uncharged test particle follows the same identical trajectory independent of its internal structure and composition. Since the polarization of a photon is considered to be part of its internal structure, we propose that polarized photons from astrophysical transients, such as gamma-ray bursts (GRBs) and fast radio bursts (FRBs), can be used to constrain the accuracy of the WEP through the Shapiro time delay effect. Assuming that the arrival time delays of photons with different polarizations are mainly attributed to the gravitational potential of the Laniakea supercluster of galaxies, we show that a strict upper limit on the differences of the parametrized post-Newtonian parameter $\gamma$ value for the polarized optical emission of GRB 120308A is $\Delta\gamma<1.2\times10^{-10}$, for the polarized gamma-ray emission of GRB 100826A is $\Delta\gamma<1.2\times10^{-10}$, and for the polarized radio emission of FRB 150807 is $\Delta\gamma<2.2\times10^{-16}$. These are the first direct verifications of the WEP for multiband photons with different polarizations. In particular, the result from FRB 150807 provides the most stringent limit to date on a deviation from the WEP, improving by one order of magnitude the previous best result based on Crab pulsar photons with different energies.

We calculate the equation of state of strange quark matter and the interior structure of strange quark stars in a Dyson-Schwinger quark model within rainbow or Ball-Chiu vertex approximation. We emphasize constraints on the parameter space of the model due to stability conditions of ordinary nuclear matter. Respecting these constraints, we find that the maximum mass of strange quark stars is about 1.9 solar masses, and typical radii are 9 --11 km. We obtain an energy release as large as $3.6 \times 10^{53}\,\text{erg}$ from conversion of neutron stars into strange quark stars.

We study rapidly rotating hybrid stars with the Dyson-Schwinger model for quark matter and the Brueckner-Hartree-Fock many-body theory with realistic two-body and three-body forces for nuclear matter. We determine the maximum gravitational mass, equatorial radius, and rotation frequency of stable stellar configurations by considering the constraints of the Keplerian limit and the secular axisymmetric instability, and compare with observational data. We also discuss the rotational evolution for constant baryonic mass, and find a spinup phenomenon for supramassive stars before they collapse to black holes.

In this paper, we propose an improved model-independent method to constrain the cosmic curvature by combining the most recent Hubble parameter $H(z)$ and supernovae Ia (SNe Ia) data. Based on the $H(z)$ data, we first use the model-independent smoothing technique, Gaussian processes, to construct distance modulus $\mu_{H}(z)$, which is susceptible to the cosmic curvature parameter $\Omega_{k}$. In contrary to previous studies, the light-curve fitting parameters, which account for distance estimation of SN ($\mu_{SN}(z)$), are set free to investigate whether $\Omega_{k}$ has a dependence on them. By comparing $\mu_{H}(z)$ to $\mu_{SN}(z)$, we put limits on $\Omega_{k}$. Our results confirm that $\Omega_{k}$ is independent of the SN light-curve parameters. Moreover, we show that the measured $\Omega_{k}$ is in good agreement with zero cosmic curvature, implying that there is no significant deviation from a flat Universe at the current observational data level. We also test the influence of different $H(z)$ samples and different Hubble constant $H_{0}$ values, finding that different $H(z)$ samples do not present significant impact on the constraints. However, different $H_{0}$ priors can affect the constraints of $\Omega_{k}$ in some degree. The prior of $H_{0}=73.24\pm1.74$ km $\rm s^{-1}$ $\rm Mpc^{-1}$ gives a value of $\Omega_{k}$ a little bit above $1\sigma$ confidence level away from 0, but $H_{0}=69.6\pm0.7$ km $\rm s^{-1}$ $\rm Mpc^{-1}$ gives it below $1\sigma$.

Broad-lined type Ic supernovae (SNe Ic-BL) are peculiar stellar explosions that distinguish themselves from ordinary SNe. Some SNe Ic-BL are associated with long-duration (\gtrsim 2 s) gamma-ray bursts (GRBs). Black holes and magnetars are two types of compact objects that are hypothesized to be central engines of GRBs. In spite of decades of investigations, no direct evidence for the formation of black holes or magnetars has been found for GRBs so far. Here we report the finding that the early peak (t \lesssim 50 days) and late-time (t \gtrsim 300 days) slow decay displayed in the light curves of both SNe 1998bw (associated with GRB 980425) and 2002ap (not GRB-associated) can be attributed to magnetar spin-down with initial rotation period P0 \sim 20 ms, while the intermediate-time (50 \lesssim t \lesssim 300 days) linear decline is caused by radioactive decay of 56Ni. The connection between the early peak and late-time slow decline in the light curves is unexpected in alternative models. We thus suggest that GRB 980425 and SN 2002ap were powered by magnetars.

Possible violations of Lorentz invariance (LIV) have been investigated for a long time using the observed spectral lags of gamma-ray bursts (GRBs). However, these generally have relied on using a single photon in the highest energy range. Furthermore, the search for LIV lags has been hindered by our ignorance concerning the intrinsic time lag in different energy bands. GRB 160625B, the only burst so far with a well-defined transition from $positive$ lags to $negative$ lags provides a unique opportunity to put new constraints on LIV. Using multi-photon energy bands we consider the contributions to the observed spectral lag from both the intrinsic time lag and the lag by LIV effects, and assuming the intrinsic time lag to have a positive dependence on the photon energy, we obtain robust limits on LIV by directly fitting the spectral lag data of GRB 160625B. Here we show that these robust limits on the quantum gravity energy scales are $E_{\rm QG,1}\geq0.5\times10^{16}$ GeV for the linear, and $E_{\rm QG,2}\geq1.4\times10^{7}$ GeV for the quadratic LIV effects, respectively. In addition, we give for the first time a reasonable formulation of the intrinsic energy-dependent time lag.

Whereas many measurements in cosmology depend on the use of integrated distances or time, galaxies evolving passively on a time scale much longer than their age difference allow us to determine the expansion rate H(z) solely as a function of the redshift-time derivative dz/dt. These model-independent `cosmic chronometers' can therefore be powerful discriminators for testing different cosmologies. In previous applications, the available sources strongly disfavoured models (such as LambdaCDM) predicting a variable acceleration, preferring instead a steady expansion rate over the redshift range 0 < z < 2. A more recent catalog of 30 objects appears to suggest non-steady expansion. In this paper, we show that such a result is entirely due to the inclusion of a high, locally-inferred value of the Hubble constant H_0 as an additional datum in a set of otherwise pure cosmic-chronometer measurements. This H_0, however, is not the same as the background Hubble constant if the local expansion rate is influenced by a Hubble Bubble. Used on their own, the cosmic chronometers completely reverse this conclusion, favouring instead a constant expansion rate out to z ~ 2.

Broad-lined type Ic supernovae (SNe Ic-BL) are of great importance because their association with long-duration gamma-ray bursts (LGRBs) holds the key to deciphering the central engine of LGRBs, which refrains from being unveiled despite decades of investigation. Among the two popularly hypothesized types of central engine, i.e., black holes and strongly magnetized neutron stars (magnetars), there is mounting evidence that the central engine of GRB-associated SNe (GRB-SNe) is rapidly rotating magnetars. Theoretical analysis also suggests that magnetars could be the central engine of SNe Ic-BL. What puzzled the researchers is the fact that light curve modeling indicates that as much as 0.2-0.5 solar mass of 56Ni was synthesized during the explosion of the SNe Ic-BL, which is unfortunately in direct conflict with current state-of-the-art understanding of magnetar-powered 56Ni synthesis. Here we propose a dynamic model of magnetar-powered SNe to take into account the acceleration of the ejecta by the magnetar, as well as the thermalization of the injected energy. Assuming that the SN kinetic energy comes exclusively from the magnetar acceleration, we find that although a major fraction of the rotational energy of the magnetar is to accelerate the SNe ejecta, a tiny fraction of this energy deposited as thermal energy of the ejecta is enough to reduce the needed 56Ni to 0.06 solar mass for both SN 1997ef and SN 2007ru. We therefore suggest that magnetars could power SNe Ic-BL both in aspects of energetics and of 56Ni synthesis.

The origin of spin of low-mass supermassive black hole (SMBH) is still a puzzle at present. We here report a study on the host galaxies of a sample of radio-selected nearby ($z<0.05$) Seyfert 2 galaxies with a BH mass of $10^{6-7} M_\odot$. By modeling the SDSS $r$-band images of these galaxies through a 2-dimensional bulge+disk decomposition, we identify a new dependence of SMBH's radio power on host bulge surface brightness profile, in which more powerful radio emission comes from a SMBH associated with a more disk-like bulge. This result means low-mass and high-mass SMBHs are spun up by two entirely different modes that correspond to two different evolutionary paths. A low-mass SMBH is spun up by a gas accretion with significant disk-like rotational dynamics of the host galaxy in the secular evolution, while a high-mass one by a BH-BH merger in the merger evolution.

To take advantage of the astrophysical potential of Gamma-Ray Bursts (GRBs), Chinese and French astrophysicists have engaged the SVOM mission (Space-based multi-band astronomical Variable Objects Monitor). Major advances in GRB studies resulting from the synergy between space and ground observations, the SVOM mission implements space and ground instrumentation. The scientific objectives of the mission put a special emphasis on two categories of GRBs: very distant GRBs at z$>$5 which constitute exceptional cosmological probes, and faint/soft nearby GRBs which allow probing the nature of the progenitors and the physics at work in the explosion. These goals have a major impact on the design of the mission: the on-board hard X-ray imager is sensitive down to 4 keV and computes on line image and rate triggers, and the follow-up telescopes on the ground are sensitive in the NIR. At the beginning of the next decade, SVOM will be the main provider of GRB positions and spectral parameters on very short time scale. The SVOM instruments will operate simultaneously with a wide range of powerful astronomical devices. This rare instrumental conjunction, combined with the relevance of the scientific topics connected with GRB studies, warrants a remarkable scientific return for SVOM. In addition, the SVOM instrumentation, primarily designed for GRB studies, composes a unique multi-wavelength observatory with rapid slew capability that will find multiple applications for the whole astronomy community beyond the specific objectives linked to GRBs. This report lists the scientific themes that will benefit from observations made with SVOM, whether they are specific GRB topics, or more generally all the issues that can take advantage of the multi-wavelength capabilities of SVOM.

To take advantage of the astrophysical potential of Gamma-Ray Bursts (GRBs), Chinese and French astrophysicists have engaged the SVOM mission (Space-based multi-band astronomical Variable Objects Monitor). Major advances in GRB studies resulting from the synergy between space and ground observations, the SVOM mission implements space and ground instrumentation. The scientific objectives of the mission put a special emphasis on two categories of GRBs: very distant GRBs at z$>$5 which constitute exceptional cosmological probes, and faint/soft nearby GRBs which allow probing the nature of the progenitors and the physics at work in the explosion. These goals have a major impact on the design of the mission: the on-board hard X-ray imager is sensitive down to 4 keV and computes on line image and rate triggers, and the follow-up telescopes on the ground are sensitive in the NIR. At the beginning of the next decade, SVOM will be the main provider of GRB positions and spectral parameters on very short time scale. The SVOM instruments will operate simultaneously with a wide range of powerful astronomical devices. This rare instrumental conjunction, combined with the relevance of the scientific topics connected with GRB studies, warrants a remarkable scientific return for SVOM. In addition, the SVOM instrumentation, primarily designed for GRB studies, composes a unique multi-wavelength observatory with rapid slew capability that will find multiple applications for the whole astronomy community beyond the specific objectives linked to GRBs. This report lists the scientific themes that will benefit from observations made with SVOM, whether they are specific GRB topics, or more generally all the issues that can take advantage of the multi-wavelength capabilities of SVOM.

A conservative constraint on the rest mass of the photon can be estimated under the assumption that the frequency dependence of dispersion from astronomical sources is mainly contributed by the nonzero photon mass effect. Photon mass limits have been earlier set through the optical emissions of the Crab Nebula pulsar, but we prove that these limits can be significantly improved with the dispersion measure (DM) measurements of radio pulsars in the Large and Small Magellanic Clouds. The combination of DM measurements of pulsars and distances of the Magellanic Clouds provide a strict upper limit on the photon mass as low as $m_{\gamma} \leq2.0\times10^{-45}~\rm{g}$, which is at least four orders of magnitude smaller than the constraint from the Crab Nebula pulsar. Although our limit is not as tight as the current best result ($\sim10^{-47}~\rm{g}$) from a fast radio burst (FRB 150418) at a cosmological distance, the cosmological origin of FRB 150418 remains under debate; and our limit can reach the same high precision of FRB 150418 when it has an extragalactic origin ($\sim10^{-45}~\rm{g}$).

We report our observations of very bright prompt optical and reverse shock (RS) optical emission of GRB 140512A and analyze its multi-wavelength data observed with the {\em Swift} and {\em Fermi} missions. It is found that the joint optical-X-ray-gamma-ray spectrum with our first optical detection (R=13.09 mag) at $T_0+136$ seconds during the second episode of the prompt gamma-rays can be fit by a single power-law with index $-1.32\pm 0.01$. Our empirical fit to the afterglow lightcurves indicates that the observed bright optical afterglow with R=13.00 mag at the peak time is consistent with predictions of the RS and forward shock (FS) emission of external shock models. Joint optical-X-ray afterglow spectrum is well fit with an absorbed single power-law, with an index evolving with time from $-1.86\pm 0.01$ at the peak time to $-1.57\pm 0.01$ at late epoch, which could be due to the evolution of the ratio of the RS to FS emission fluxes. We fit the lightcurves with standard external models, and derive the physical properties of the outflow. It is found that the ratio $R_B\equiv\epsilon_{\rm B, r}/\epsilon_{\rm B, f}$ is 8187, indicating a high magnetization degree in the RS region. Measuring the relative radiation efficiency with $R_e\equiv\epsilon_{\rm e, r}/\epsilon_{\rm e, f}$, we have $R_e= 0.02$, implying the radiation efficiency of the RS is much lower than that in FS. We also show that the $R_B$ of GRBs 990123, 090102, and 130427A are similar to that of GRB 140512A and their apparent difference may be mainly attributed to the difference of the jet kinetic energy, initial Lorentz factor, and medium density among them.

Five TeV neutrino events weakly correlated with five gamma-ray bursts (GRBs) were detected recently by IceCube. This work is an attempt to show that if the GRB identifications are verified, the observed time delays between the TeV neutrinos and gamma-ray photons from GRBs provide attractive candidates for testing fundamental physics with high accuracy. Based on the assumed associations between the TeV neutrinos and GRBs, we find that the limiting velocity of the neutrinos is equal to that of photons to an accuracy of $\sim1.9\times10^{-15}-2.5\times10^{-18}$, which is about $10^{4}-10^{7}$ times better than the constraint obtained with the neutrino possibly from a blazar flare. In addition, we set the most stringent limits up to date on the energy scale of quantum gravity for both the linear and quadratic violations of Lorentz invariance, namely $E_{\rm QG, 1}>6.3\times10^{18}-1.5\times10^{21}$ GeV and $E_{\rm QG, 2}>2.0\times10^{11}-4.2\times10^{12}$ GeV, which are essentially as good as or are an improvement of one order of magnitude over the results previously obtained by the GeV photons of GRB 090510 and the PeV neutrino from a blazar flare. Assuming that the Shapiro time delay is caused by the gravitational potential of the Laniakea supercluster of galaxies, we also place the tightest limits to date on Einstein's weak equivalence principle through the relative differential variations of the parameterized post-Newtonian parameter $\gamma$ values for two different species of particles (i.e., neutrinos and photons), yielding $\Delta\gamma \sim 10^{-11}-10^{-13}$. However, it should be emphasized again that these limits here obtained are at best forecast of what could be achieved if the GRB/neutrino correlations would be finally confirmed.

We continue to build support for the proposal to use HII galaxies (HIIGx) and giant extragalactic HII regions (GEHR) as standard candles to construct the Hubble diagram at redshifts beyond the current reach of Type Ia supernovae. Using a sample of 25 high-redshift HIIGx, 107 local HIIGx, and 24 GEHR, we confirm that the correlation between the emission-line luminosity and ionized-gas velocity dispersion is a viable luminosity indicator, and use it to test and compare the standard model $\Lambda$CDM and the $R_{\rm h}=ct$ Universe by optimizing the parameters in each cosmology using a maximization of the likelihood function. For the flat $\Lambda$CDM model, the best fit is obtained with $\Omega_{\rm m}= 0.40_{-0.09}^{+0.09}$. However, statistical tools, such as the Akaike (AIC), Kullback (KIC) and Bayes (BIC) Information Criteria favor $R_{\rm h}=ct$ over the standard model with a likelihood of $\approx 94.8\%-98.8\%$ versus only $\approx 1.2\%-5.2\%$. For $w$CDM (the version of $\Lambda$CDM with a dark-energy equation of state $w_{\rm de}\equiv p_{\rm de}/\rho_{\rm de}$ rather than $w_{\rm de}=w_{\Lambda}=-1$), a statistically acceptable fit is realized with $\Omega_{\rm m}=0.22_{-0.14}^{+0.16}$ and $w_{\rm de}= -0.51_{-0.25}^{+0.15}$ which, however, are not fully consistent with their concordance values. In this case, $w$CDM has two more free parameters than $R_{\rm h}=ct$, and is penalized more heavily by these criteria. We find that $R_{\rm h}=ct$ is strongly favored over $w$CDM with a likelihood of $\approx 92.9\%-99.6\%$ versus only $0.4\%-7.1\%$. The current HIIGx sample is already large enough for the BIC to rule out $\Lambda$CDM/$w$CDM in favor of $R_{\rm h}=ct$ at a confidence level approaching $3\sigma$.

A conservative constraint on the Einstein Weak Equivalence Principle (WEP) can be obtained under the assumption that the observed time delay between correlated particles from astronomical sources is dominated by the gravitational fields through which they move. Current limits on the WEP are mainly based on the observed time delays of photons with different energies. It is highly desirable to develop more accurate tests that include the gravitational wave (GW) sector. The detection by the advanced LIGO/VIRGO systems of gravitational waves will provide attractive candidates for constraining the WEP, extending the tests to gravitational interactions, with potentially higher accuracy. Considering the capabilities of the advanced LIGO/VIRGO network and the source direction uncertainty, we show that the joint detection of GWs and electromagnetic signals could probe the WEP to an accuracy down to $10^{-10}$, which is one order of magnitude tighter than previous limits, and seven orders of magnitude tighter than the multi-messenger (photons and neutrinos) results by supernova 1987A.

We report our very early optical observations of GRB 110530A and investigate its jet properties together with its X-ray afterglow data. A peculiar broad onset bump followed by a plateau is observed in its early R band afterglow lightcurve. The optical data in the other bands and the X-ray data are well consistent with the temporal feature of the R band lightcurve. Our joint spectral fits of the optical and X-ray data show that they are in the same regime, with a photon index of $\sim 1.70$. The optical and X-ray afterglow lightcurves are well fitted with the standard external shock model by considering a delayed energy injection component. Based on our modeling results, we find that the radiative efficiency of the GRB jet is $\sim 1\%$ and the magnetization parameter of the afterglow jet is $<0.04$ with the derived extremely low $\epsilon_B$ (the fraction of shock energy to magnetic field) of $(1.64\pm 0.25)\times 10^{-6}$. These results indicate that the jet may be matter dominated. Discussion on delayed energy injection from accretion of late fall-back material of its pre-supernova star is also presented.

We briefly introduce the disadvantages for Type Ia supernovae (SNe Ia) as standard candles to measure the Universe, and suggest Gamma-ray bursts (GRBs) can serve as a powerful tool for probing the properties of high redshift Universe. We use GRBs as distance indicators in constructing the Hubble diagram at redshifts beyond the current reach of SNe Ia observations. Since the progenitors of long GRBs are confirmed to be massive stars, they are deemed as an effective approach to study the cosmic star formation rate (SFR). A detailed representation of how to measure high-$z$ SFR using GRBs is presented. Moreover, first stars can form only in structures that are suitably dense, which can be parameterized by defining the minimum dark matter halo mass $M_{\rm min}$. $M_{\rm min}$ must play a crucial role in star formation. The association of long GRBs with the collapses of massive stars also indicates that the GRB data can be applied to constrain the minimum halo mass $M_{\rm min}$ and to investigate star formation in dark matter halos.

We present a star catalogue extracted from the Lunar-based Ultraviolet Telescope (LUT) survey program. LUT's observable sky area is a circular belt around the Moon's north pole, and the survey program covers a preferred area for about 2400 deg$^2$ which includes a region of the Galactic plane. The data is processed with an automatic pipeline which copes with stray light contamination, artificial sources, cosmic rays, flat field calibration, photometry and so on. In the first release version, the catalogue provides high confidence sources which have been cross-identified with Tycho-2 catalogue. All the sources have signal-to-noise ratio larger than 5, and the corresponding magnitude limit is typically 14.4 mag, which can be deeper as ~16 mag if the stray light contamination is in the lowest level. A total number of 86,467 stars are recorded in the catalogue. The full catalogue in electronic form is available on line.

High-redshift Gamma-Ray Bursts (GRBs) beyond redshift $\sim6$ are potentially powerful tools to probe the distant early Universe. Their detections in large numbers and at truly high redshifts call for the next generation of high-energy wide-field instruments with unprecedented sensitivity at least one order of magnitude higher than the ones currently in orbit. On the other hand, follow-up observations of the afterglows of high-redshift GRBs and identification of their host galaxies, which would be difficult for the currently operating telescopes, require new, extremely large facilities of at multi-wavelengths. This chapter describes future experiments that are expected to advance this exciting field, both being currently built and being proposed. The legacy of Swift will be continued by SVOM, which is equipped with a set of space-based multi-wavelength instruments as well as and a ground segment including a wide angle camera and two follow-up telescopes. The established Lobster-eye X-ray focusing optics provides a promising technology for the detection of faint GRBs at very large distances, based on which the {THESEUS}, {Einstein Probe} and other mission concepts have been proposed. Follow-up observations and exploration of the reionization era will be enabled by large facilities such as {SKA} in the radio, the 30m class telescopes in the optical/near-IR, and the space-borne {WFIRST} and {JWST} in the optical/near-IR/mid-IR. In addition, the X-ray and $\gamma$-ray polarization experiment POLAR is also introduced.

Based on the star formation histories (SFH) of galaxies in halos of different masses, we develop an empirical model to grow galaxies in dark mattet halos. This model has very few ingredients, any of which can be associated to observational data and thus be efficiently assessed. By applying this model to a very high resolution cosmological $N$-body simulation, we predict a number of galaxy properties that are a very good match to relevant observational data. Namely, for both centrals and satellites, the galaxy stellar mass function (SMF) up to redshift $z\simeq4$ and the conditional stellar mass functions (CSMF) in the local universe are in good agreement with observations. In addition, the 2-point correlation is well predicted in the different stellar mass ranges explored by our model. Furthermore, after applying stellar population synthesis models to our stellar composition as a function of redshift, we find that the luminosity functions in $^{0.1}u$, $^{0.1}g$, $^{0.1}r$, $^{0.1}i$ and $^{0.1}z$ bands agree quite well with the SDSS observational results down to an absolute magnitude at about -17.0. The SDSS conditional luminosity functions (CLF) itself is predicted well. Finally, the cold gas is derived from the star formation rate (SFR) to predict the HI gas mass within each mock galaxy. We find a remarkably good match to observed HI-to-stellar mass ratios. These features ensure that such galaxy/gas catalogs can be used to generate reliable mock redshift surveys.

Fast radio bursts (FRBs) are radio bursts characterized by millisecond durations, high Galactic latitude positions, and high dispersion measures. Very recently, the cosmological origin of FRB 150418 has been confirmed by \cite{kea16}, and FRBs are now strong competitors as cosmological probes. The simple sharp feature of the FRB signal is ideal for them to probe some of the fundamental laws of physics. Here we show that by analyzing the delay time between different frequencies, the FRB data can place stringent upper limits on the rest mass of the photon. For FRB 150418 at $z=0.492$, one can potentially reach $m_{\gamma} \leq 5.2 \times 10^{-47}$ g, which is $10^{20}$ times smaller than the rest mass of electron, and is about $10^3$ times smaller than that obtained using other astrophysical sources with the same method.

Comparing the parameterized post-Newtonian parameter $\gamma$ values for different types of particles, or the same type of particles with different energies is an important method to test the Einstein Equivalence Principle (EEP). Assuming that the observed time delays are dominated by the gravitational potential of the Laniakea supercluster of galaxies, better results of EEP constraints can be obtained. In this paper, we apply photons from three kinds of cosmic transients, including TeV blazars, gamma-ray bursts as well as fast radio bursts to constrain EEP. With a gravitational field far more stronger than a single galaxy, we obtain 4--5 orders of magnitude more stringent than the pervious results.

Millisecond magnetars can be formed via several channels: core-collapse of massive stars, accretion-induced collapse of white dwarfs (WDs), double WD mergers, double neutron star (NS) mergers, and WD-NS mergers. Because the mass of ejecta from these channels could be quite different, their light curves are also expected to be diverse. We evaluate the dynamic evolution of optical transients powered by millisecond magnetars. We find that the magnetar with short spin-down timescale converts its rotational energy mostly into the kinetic energy of the transient, while the energy of a magnetar with long spin-down timescale goes into radiation of the transient. This leads us to speculate that hypernovae could be powered by magnetars with short spin-down timescales. At late times the optical transients will gradually evolve into a nebular phase because of the photospheric recession. We treat the photosphere and nebula separately because their radiation mechanisms are different. In some cases the ejecta could be light enough that the magnetar can accelerate it to a relativistic speed. It is well known that the peak luminosity of a supernova (SN) occurs when the luminosity is equal to the instantaneous energy input rate, as shown by Arnett (1979). We show that photospheric recession and relativistic motion can modify this law. The photospheric recession always leads to a delay of the peak time $t_{\mathrm{pk}}$ relative to the time $t_{\times }$ at which the SN luminosity equals the instantaneous energy input rate. Relativistic motion, however, may change this result significantly.

The cosmic evolution of gamma-ray burst (GRB) luminosity is essential for revealing the GRB physics and for using GRBs as cosmological probes. We investigate the luminosity evolution of long GRBs with a large sample of 258 {\em Swift}/BAT GRBs. Parameterized the peak luminosity of individual GRBs evolves as $L_{\rm p}\propto{\rm }(1+z)^{k}$, we get $k=1.49\pm0.19$ using the non-parametric $\tau$ statistics method without considering observational biases of GRB trigger and redshift measurement. By modeling these biases with the observed peak flux and characterizing the peak luminosity function of long GRBs as a smoothly broken power-law with a break that evolves as $L_{\rm b}\propto (1+z)^{k_{\rm b}}$, we obtain $k_{\rm b}=1.14^{+0.99}_{-0.47}$ through simulations based on assumption that the long GRB rate follows the star formation rate (SFR) incorporating with cosmic metallicity history. The derived $k$ and $k_b$ values are systematically smaller than that reported in previous papers. By removing the observational biases of the GRB trigger and redshift measurement based on our simulation analysis, we generate mock {\em complete} samples of 258 and 1000 GRBs to examine how these biases affects on the $\tau$ statistics method. We get $k=0.94\pm 0.14$ and $k=0.80\pm 0.09$ for the two samples, indicating that these observational biases may lead to overestimate the $k$ value. With the large uncertain of $k_b$ derived from our simulation analysis, one even cannot convincingly argue a robust evolution feature of the GRB luminosity.

Gamma-ray bursts (GRBs) are short and intense flashes at the cosmological distances, which are the most luminous explosions in the Universe. The high luminosities of GRBs make them detectable out to the edge of the visible universe. So, they are unique tools to probe the properties of high-redshift universe: including the cosmic expansion and dark energy, star formation rate, the reionization epoch and the metal evolution of the Universe. First, they can be used to constrain the history of cosmic acceleration and the evolution of dark energy in a redshift range hardly achievable by other cosmological probes. Second, long GRBs are believed to be formed by collapse of massive stars. So they can be used to derive the high-redshift star formation rate, which can not be probed by current observations. Moreover, the use of GRBs as cosmological tools could unveil the reionization history and metal evolution of the Universe, the intergalactic medium (IGM) properties and the nature of first stars in the early universe. But beyond that, the GRB high-energy photons can be applied to constrain Lorentz invariance violation (LIV) and to test Einstein's Equivalence Principle (EEP). In this paper, we review the progress on the GRB cosmology and fundamental physics probed by GRBs.

The observed time delays between different energy bands from TeV blazars provide a new interesting way of testing the Einstein Equivalence Principle (EEP). If the whole time delay is assumed to be dominated by the gravitational field of the Milky Way, the conservative upper limit on the EEP can be estimated. Here we show that the strict limits on the differences of the parameterized post-Newtonian parameter $\gamma$ values are $\gamma_{\rm TeV}-\gamma_{\rm keV}<3.86\times10^{-3}$ for Mrk 421 and $\gamma_{\rm TeV}-\gamma_{\rm keV}<4.43\times10^{-3}$ for Mrk 501, while expanding the scope of the tested EEP energy range out to the TeV--keV range for the first time. With the small time lag from the 0.2--0.8 TeV and $>0.8$ TeV light curves of PKS 2155-304, a much more severe constraint on $\gamma$ differences of $\sim10^{-6}$ can be achieved, although the energy difference is of order of $\sim$ TeV. Furthermore, we can combine these limits on the energy dependence of $\gamma$ with the bound on the absolute $\gamma$ value $\gamma-1\sim0.3\%$ from light deflection measurements at the optical (eV) bands, and conclude that this absolute bound on $\gamma$ can be extended from optical to TeV energies.

Observations of very early multi-wavelength afterglows are critical to reveal the properties of the radiating fireball and its environment as well as the central engine of gamma-ray bursts (GRBs). We report our optical observations of GRB 111228A from 95 sec to about 50 hours after the burst trigger and investigate its properties of the prompt gamma-rays and the ambient medium using our data and the data observed with {\em Swift} and {\em Fermi} missions. Our joint optical and X-ray spectral fits to the afterglow data show that the ambient medium features as low dust-to-gas ratio. Incorporating the energy injection effect, our best fit to the afterglow lightcurves with the standard afterglow model via the Markov Chain Monte Carlo (MCMC) technique shows that $\epsilon_e=(6.9\pm 0.3)\times 10^{-2}$, $\epsilon_B=(7.73\pm 0.62)\times 10^{-6}$, $E_{\rm K}=(6.32\pm 0.86)\times 10^{53}\rm erg$, $n=0.100\pm 0.014$ cm$^{-3}$. The low medium density likely implies that the afterglow jet may be in a halo or in a hot ISM. Achromatic shallow decay segment observed in the optical and X-ray bands is well explained with the long-lasting energy injection from the central engine, which would be a magnetar with a period of about 1.92 ms inferred from the data. The $E_p$ of its time-integrated prompt gamma-ray spectrum is $\sim 26$ KeV. Using the initial Lorentz factor ($\Gamma_0=476^{+225}_{-237}$) derived from our afterglow model fit, it is found that GRB 111228A satisfies the $L_{\rm iso}-E_{\rm p,z}-\Gamma_0$ relation and bridges the typical GRBs and low luminosity GRBs in this relation.

The accuracy of Einstein's Equivalence Principle (EEP) can be tested with the observed time delays between correlated particles or photons that are emitted from astronomical sources. Assuming as a lower limit that the time delays are caused mainly by the gravitational potential of the Milky Way, we prove that fast radio bursts (FRBs) of cosmological origin can be used to constrain the EEP with high accuracy. Taking FRB 110220 and two possible FRB/gamma-ray burst (GRB) association systems (FRB/GRB 101011A and FRB/GRB 100704A) as examples, we obtain a strict upper limit on the differences of the parametrized post-Newtonian parameter $\gamma$ values as low as $\left[\gamma(1.23\; \rm GHz)-\gamma(1.45\; \rm GHz)\right]<4.36\times10^{-9}$. This provides the most stringent limit up to date on the EEP through the relative differential variations of the $\gamma$ parameter at radio energies, improving by 1 to 2 orders of magnitude the previous results at other energies based on supernova 1987A and GRBs.

We briefly present the science capabilities, the instruments, the operations, and the expected performance of the SVOM mission. SVOM (Space-based multiband astronomical Variable Objects Monitor) is a Chinese-French space mission dedicated to the study of Gamma-Ray Bursts (GRBs) in the next decade. The SVOM mission encompasses a satellite carrying four instruments to detect and localize the prompt GRB emission and measure the evolution of the afterglow in the visible band and in X-rays, a VHF communication system enabling the fast transmission of SVOM alerts to the ground, and a ground segment including a wide angle camera and two follow-up telescopes. The pointing strategy of the satellite has been optimized to favor the detection of GRBs located in the night hemisphere. This strategy enables the study of the optical emission in the first minutes after the GRB with robotic observatories and the early spectroscopy of the optical afterglow with large telescopes to measure the redshifts. The study of GRBs in the next decade will benefit from a number of large facilities in all wavelengths that will contribute to increase the scientific return of the mission. Finally, SVOM will operate in the era of the next generation of gravitational wave detectors, greatly contributing to searches for the electromagnetic counterparts of gravitational wave triggers at Xray and gamma-ray energies.

The origin of outflow in narrow-line region (NLR) of active galactic nucleus (AGN) is studied in this paper by focusing on the relationship between the [\ion{O}{3}]$\lambda$5007 line profile and the hard X-ray (in a bandpass of 2-10 keV) emission from the central SMBH in type-I AGNs. A sample of 47 local X-ray selected type-I AGNs at $z<0.2$ is extracted from the 2XMMi/SDSS DR7 catalog that is originally crossmatched by Pineau et al. The X-ray luminosities in an energy band from 2 to 10keV of these luminous AGNs range from $10^{42}$ to $10^{44}\ \mathrm{erg\ s^{-1}}$. A joint spectral analysis is performed on their optical and X-ray spectra, in which the [\ion{O}{3}] line profile is modeled by a sum of several Gaussian functions to quantify its deviation from a pure Gaussian function. The statistics allows us to identify a moderate correlation with a significance level of 2.78$\sigma$: luminous AGNs with stronger [\ion{O}{3}] blue asymmetry tend to have steeper hard X-ray spectra. By identifying a role of $L/L_{\mathrm{Edd}}$ on the correlation at a $2-3\sigma$ significance level in both direct and indirect ways, we argue that the photon index versus asymmetry correlation provides evidence that the AGN's outflow commonly observed in its NLR is related with the accretion process occurring around the central SMBH, which favors the wind/radiation model for the origin of the outflow in luminous AGNs.

First stars can only form in structures that are suitably dense, which can be parametrized by the minimum dark matter halo mass $M_{\rm min}$. $M_{\rm min}$ must plays an important role in star formation. The connection of long gamma-ray bursts (LGRBs) with the collapse of massive stars has provided a good opportunity for probing star formation in dark matter halos. We place some constraints on $M_{\rm min}$ using the latest $Swift$ LGRB data. We conservatively consider that LGRB rate is proportional to the cosmic star formation rate (CSFR) and an additional evolution parametrized as $(1+z)^{\alpha}$, where the CSFR model as a function of $M_{\rm min}$. Using the $\chi^{2}$ statistic, the contour constraints on the $M_{\rm min}$--$\alpha$ plane show that at the $1\sigma$ confidence level, we have $M_{\rm min}<10^{10.5}$ $\rm M_{\odot}$ from 118 LGRBs with redshift $z<4$ and luminosity $L_{\rm iso}>1.8\times10^{51}$ erg $\rm s^{-1}$. We also find that adding 12 high-\emph{z} $(4<z<5)$ LGRBs (consisting of 104 LGRBs with $z<5$ and $L_{\rm iso}>3.1\times10^{51}$ erg $\rm s^{-1}$) could result in much tighter constraints on $M_{\rm min}$, for which, $10^{7.7}\rm M_{\odot}<M_{\rm min}<10^{11.6}\rm M_{\odot}$ ($1\sigma$). Through Monte Carlo simulations, we estimate that future five years of Sino-French spacebased multiband astronomical variable objects monitor (\emph{SVOM}) observations would tighten these constraints to $10^{9.7}\rm M_{\odot}<M_{\rm min}<10^{11.3}\rm M_{\odot}$. The strong constraints on $M_{\rm min}$ indicate that LGRBs are a new promising tool for investigating star formation in dark matter halos.

We here report the photometric performance of Lunar-based Ultraviolet telescope (LUT), the first robotic telescope working on the Moon, for its 18-months operation. In total, 17 IUE standards have been observed in 51 runs until June 2015, which returns a highly stable photometric performance during the past 18 months (i.e., no evolution of photometric performance with time). The magnitude zero point is determined to be $17.53\pm0.05$ mag, which is not only highly consistent with the results based on its first 6-months operation, but also independent on the spectral type of the standard from which the magnitude zero point is determined. The implications of this stable performance is discussed, and is useful for next generation lunar-based astronomical observations.

We describe the data processing pipeline developed to reduce the pointing observation data of Lunar-based Ultraviolet Telescope (LUT), which belongs to the Chang'e-3 mission of the Chinese Lunar Exploration Program. The pointing observation program of LUT is dedicated to monitor variable objects in a near-ultraviolet (245-345 nm) band. LUT works in lunar daytime for sufficient power supply, so some special data processing strategies have been developed for the pipeline. The procedures of the pipeline include stray light removing, astrometry, flat fielding employing superflat technique, source extraction and cosmic rays rejection, aperture and PSF photometry, aperture correction, and catalogues archiving, etc. It has been intensively tested and works smoothly with observation data. The photometric accuracy is typically ~0.02 mag for LUT 10 mag stars (30 s exposure), with errors come from background noises, residuals of stray light removing, and flat fielding related errors. The accuracy degrades to be ~0.2 mag for stars of 13.5 mag which is the 5{\sigma} detection limit of LUT.

We reported the optical observations of GRB 121011A by 0.8-m TNT telescope at Xinglong observatory, China. The light curve of optical afterglow shows a smooth and featureless bump during the epoch of $\sim$130 sec and $\sim$5000 sec with a rising index of $1.57\pm0.28$ before the break time of $539\pm44$ sec, and a decaying index of about $1.29\pm0.07$ up to the end of our observations. Meanwhile, the X-ray light curve decays in a single power-law with a slop of about $1.51\pm0.03$ observed by $XRT$ onboard ${\rm} Swift$ from 100 sec to about 10000 sec after the burst trigger. The featureless optical light curve could be understood as an onset process under the external-shock model. The typical frequency has been below or near the optical one before the deceleration time, and the cooling frequency is located between the optical and X-ray wavelengths. The external medium density has a transition from a mixed stage of ISM and wind-type medium before the peak time to the ISM at the later phase. The joint-analysis of X-ray and optical light curves shows that the emission from both frequencies are consistent with the prediction of the standard afterglow model without any energy injections, indicating that the central engine has stopped its activity and does not restart anymore after the prompt phase.

We use 32 age measurements of passively evolving galaxies as a function of redshift to test and compare the standard model ($\Lambda$CDM) with the $R_{\rm h}=ct$ Universe. We show that the latter fits the data with a reduced $\chi^2_{\rm dof}=0.435$ for a Hubble constant $H_{0}= 67.2_{-4.0}^{+4.5}$ km $\rm s^{-1}$ $\rm Mpc^{-1}$. By comparison, the optimal flat $\Lambda$CDM model, with two free parameters (including $\Omega_{\rm m}=0.12_{-0.11}^{+0.54}$ and $H_{0}=94.3_{-35.8}^{+32.7}$ km $\rm s^{-1}$ $\rm Mpc^{-1}$), fits the age-\emph{z} data with a reduced $\chi^2_{\rm dof}=0.428$. Based solely on their $\chi^2_{\rm dof}$ values, both models appear to account for the data very well, though the optimized $\Lambda$CDM parameters are only marginally consistent with those of the concordance model ($\Omega_{\rm m}=0.27$ and $H_{0}= 70$ km $\rm s^{-1}$ $\rm Mpc^{-1}$). Fitting the age-$z$ data with the latter results in a reduced $\chi^2_{\rm dof}=0.523$. However, because of the different number of free parameters in these models, selection tools, such as the Akaike, Kullback and Bayes Information Criteria, favour $R_{\rm h}=ct$ over $\Lambda$CDM with a likelihood of $\sim 66.5\%-80.5\%$ versus $\sim 19.5\%-33.5\%$. These results are suggestive, though not yet compelling, given the current limited galaxy age-$z$ sample. We carry out Monte Carlo simulations based on these current age measurements to estimate how large the sample would have to be in order to rule out either model at a $\sim 99.7\%$ confidence level. We find that if the real cosmology is $\Lambda$CDM, a sample of $\sim 45$ galaxy ages would be sufficient to rule out $R_{\rm h}=ct$ at this level of accuracy, while $\sim 350$ galaxy ages would be required to rule out $\Lambda$CDM if the real Universe were instead $R_{\rm h}=ct$.

We update gamma-ray burst (GRB) luminosity relations among certain spectral and light-curve features with 139 GRBs. The distance modulus of 82 GRBs at $z>1.4$ can be calibrated with the sample at $z\leq1.4$ by using the cubic spline interpolation method from the Union2.1 Type Ia supernovae (SNe Ia) set. We investigate the joint constraints on the Cardassian expansion model and dark energy with 580 Union2.1 SNe Ia sample ($z<1.4$) and 82 calibrated GRBs data ($1.4<z\leq8.2$). In $\Lambda$CDM, we find that adding 82 high-\emph{z} GRBs to 580 SNe Ia significantly improves the constrain on $\Omega_{m}-\Omega_{\Lambda}$ plane. In the Cardassian expansion model, the best fit is $\Omega_{m}= 0.24_{-0.15}^{+0.15}$ and $n=0.16_{-0.52}^{+0.30}$ $(1\sigma)$, which is consistent with the $\Lambda$CDM cosmology $(n=0)$ in the $1\sigma$ confidence region. We also discuss two dark energy models in which the equation of state $w(z)$ is parametrized as $w(z)=w_{0}$ and $w(z)=w_{0}+w_{1}z/(1+z)$, respectively. Based on our analysis, we see that our Universe at higher redshift up to $z=8.2$ is consistent with the concordance model within $1\sigma$ confidence level.

The use of type Ic Super Luminous Supernovae (SLSN Ic) to examine the cosmological expansion introduces a new standard ruler with which to test theoretical models. The sample suitable for this kind of work now includes 11 SLSNe Ic, which have thus far been used solely in tests involving $\Lambda$CDM. In this paper, we broaden the base of support for this new, important cosmic probe by using these observations to carry out a one-on-one comparison between the $R_{\rm h}=ct$ and $\Lambda$CDM cosmologies. We individually optimize the parameters in each cosmological model by minimizing the $\chi^{2}$ statistic. We also carry out Monte Carlo simulations based on these current SLSN Ic measurements to estimate how large the sample would have to be in order to rule out either model at a $\sim 99.7\%$ confidence level. The currently available sample indicates a likelihood of $\sim$$70-80\%$ that the $R_{\rm h}=ct$ Universe is the correct cosmology versus $\sim$$20-30\%$ for the standard model. These results are suggestive, though not yet compelling, given the current limited number of SLSNe Ic. We find that if the real cosmology is $\Lambda$CDM, a sample of $\sim$$240$ SLSNe Ic would be sufficient to rule out $R_{\rm h}=ct$ at this level of confidence, while $\sim$$480$ SLSNe Ic would be required to rule out $\Lambda$CDM if the real Universe is instead $R_{\rm h}=ct$. This difference in required sample size reflects the greater number of free parameters available to fit the data with $\Lambda$CDM. If such SLSNe Ic are commonly detected in the future, they could be a powerful tool for constraining the dark-energy equation of state in $\Lambda$CDM, and differentiating between this model and the $R_{\rm h}=ct$ Universe.

We study cold dense quark matter and hybrid neutron stars with a Dyson-Schwinger quark model and various choices of the quark-gluon vertex. We obtain the equation of state of quark matter in beta equilibrium and investigate the hadron-quark phase transition in combination with a hadronic equation of state derived within the Brueckner-Hartree-Fock many-body theory. Comparing with the results for quark matter with the rainbow approximation, the Ball-Chiu ansatz and the 1BC ansatz for the quark-gluon vertex lead to a reduction of the effective interaction at finite chemical potential, qualitatively similar to the effect of our gluon propagator. We find that the phase transition and the equation of state of the quark or mixed phase and consequently the resulting hybrid star mass and radius depend mainly on a global reduction of the effective interaction due to effects of both the quark-gluon vertex and gluon propagator, but are not sensitive to the vertex ansatz.

We study an extensive sample of 87 GRBs for which there are well sampled and simultaneous optical and X-ray light-curves. We extract the cleanest possible signal of the afterglow component, and compare the temporal behaviors of the X-ray light-curve, observed by Swift XRT, and optical data, observed by UVOT and ground-based telescopes for each individual burst. Overall we find 62\% GRBs that are consistent with the standard afterglow model. When more advanced modeling is invoked, up to 91\% of the bursts in our sample may be consistent with the external shock model. A large fraction of these bursts are consistent with occurring in a constant interstellar density medium (ISM) (61\%) while only 39\% of them occur in a wind-like medium. Only 9 cases have afterglow light-curves that exactly match the standard fireball model prediction, having a single power law decay in both energy bands which are observed during their entire duration. In particular, for the bursts with chromatic behavior additional model assumptions must be made over limited segments of the light-curves in order for these bursts to fully agree with the external shock model. Interestingly, for 54\% of the X-ray and 40\% of the optical band observations the end of the shallow decay ($t^{\sim-0.5}$) period coincides with the jet break ($t^{\sim-p}$) time, causing an abrupt change in decay slope. The fraction of the burst that consistent with the external shock model is independent of the observational epochs in the rest frame of GRBs. Moreover, no cases can be explained by the cooling frequency crossing the X-ray or optical band.

The content of $\mathrm{OH/H_2O}$ molecules in the tenuous exosphere of the Moon is still an open issue at present. We here report an unprecedented upper limit of the content of the OH radicals, which is obtained from the in-situ measurements carried out \rm by the Lunar-based Ultraviolet Telescope, a payload of Chinese Chang'e-3 mission. By analyzing the diffuse background in the images taken by the telescope, the column density and surface concentration of the OH radicals are inferred to be $<10^{11}\ \mathrm{cm^{-2}}$ and $<10^{4}\ \mathrm{cm^{-3}}$ (by assuming a hydrostatic equilibrium with a scale height of 100km), respectively, by assuming that the recorded background is fully contributed by their resonance fluorescence emission. The resulted concentration is lower than the previously reported value by about two orders of magnitude, and is close to the prediction of the sputtering model. In addition, the same measurements and method allow us to derive a surface concentration of $<10^{2}\ \mathrm{cm^{-3}}$ for the neutral magnesium, which is lower than the previously reported upper limit by about two orders of magnitude. These results are the best known of the OH (MgI) content in the lunar exosphere to date.

The use of Type~Ia SNe has thus far produced the most reliable measurement of the expansion history of the Universe, suggesting that $\Lambda$CDM offers the best explanation for the redshift--luminosity distribution observed in these events. But the analysis of other kinds of source, such as cosmic chronometers, gamma ray bursts, and high-$z$ quasars, conflicts with this conclusion, indicating instead that the constant expansion rate implied by the $R_{\rm h}=ct$ Universe is a better fit to the data. The central difficulty with the use of Type~Ia SNe as standard candles is that one must optimize three or four nuisance parameters characterizing supernova luminosities simultaneously with the parameters of an expansion model. Hence in comparing competing models, one must reduce the data independently for each. We carry~out such a comparison of $\Lambda$CDM and the $R_{\rm h}=ct$ Universe, using the Supernova Legacy Survey (SNLS) sample of 252 SN~events, and show that each model fits its individually reduced data very well. But since $R_{\rm h}=ct$ has only one free parameter (the Hubble constant), it follows from a standard model selection technique that it is to be preferred over $\Lambda$CDM, the minimalist version of which has three (the Hubble constant, the scaled matter density and either the spatial curvature constant or the dark-energy equation-of-state parameter). We estimate by the Bayes Information Criterion that in a pairwise comparison, the likelihood of $R_{\rm h}=ct$ is $\sim 90\%$, compared with only $\sim 10\%$ for a minimalist form of $\Lambda$CDM, in which dark energy is simply a cosmological constant. Compared to $R_{\rm h}=ct$, versions of the standard model with more elaborate parametrizations of dark energy are judged to be even less likely.

We reported the photometric calibration of Lunar-based Ultraviolet telescope (LUT), the first robotic astronomical telescope working on the lunar surface, for its first six months of operation on the lunar surface. Two spectral datasets (set A and B) from near-ultraviolet (NUV) to optical band were constructed for 44 International Ultraviolet Explorer (IUE) standards, because of the LUT's relatively wide wavelength coverage. Set A were obtained by extrapolating the IUE NUV spectra ($\lambda<3200\AA$) to optical band basing upon the theoretical spectra of stellar atmosphere models. Set B were exactly the theoretical spectra from 2000\AA\ to 8000\AA\ extracted from the same model grid. In total, seven standards have been observed in 15 observational runs until May 2014. The calibration results show that the photometric performance of LUT is highly stable in its first six months of operation. The magnitude zero points obtained from the two spectral datasets are also consistent with each other, i.e., $\mathrm{zp=17.54\pm0.09}$mag (set A) and $\mathrm{zp=17.52\pm0.07}$mag (set B).

We use measurements of the galaxy-cluster angular size versus redshift to test and compare the standard model (LCDM) and the R_h=ct Universe. We show that the latter fits the data with a reduced chi^2_dof=0.786 for a Hubble constant H_0= 72.6 (-3.4+3.8) km/s/Mpc, and H_0 is the sole parameter in this model. By comparison, the optimal flat LCDM model, with two free parameters (including Omega_m=0.50 and H_0=73.9 (-9.5+10.6) km/s/Mpc), fits the angular-size data with a reduced chi^2_dof=0.806. On the basis of their chi^2_dof values alone, both models appear to account for the data very well in spite of the fact that the R_h=ct Universe expands at a constant rate, while LCDM does not. However, because of the different number of free parameters in these models, selection tools, such as the Bayes Information Criterion, favour R_h=ct over LCDM with a likelihood of ~86% versus ~14%. These results impact the question of galaxy growth at large redshifts. Previous work suggested an inconsistency with the underlying cosmological model unless elliptical and disk galaxies grew in size by a surprisingly large factor ~6 from z~3 to 0. The fact that both LCDM and R_h=ct fit the cluster-size measurements quite well casts some doubt on the suggestion that the unexpected result with individual galaxies may be due to the use of an incorrect expansion scenario, rather than astrophysical causes, such as mergers and/or selection effects.

Strongly gravitationally lensed quasar-galaxy systems allow us to compare competing cosmologies as long as one can be reasonably sure of the mass distribution within the intervening lens. In this paper, we assemble a catalog of 69 such systems, and carry out a one-on-one comparison between the standard model, LCDM, and the R_h=ct Universe. We find that both models account for the lens observations quite well, though the precision of these measurements does not appear to be good enough to favor one model over the other. Part of the reason is the so-called bulge-halo conspiracy that, on average, results in a baryonic velocity dispersion within a fraction of the optical effective radius virtually identical to that expected for the whole luminous-dark matter distribution. Given the limitations of doing precision cosmological testing using the current sample, we also carry out Monte Carlo simulations based on the current lens measurements to estimate how large the source catalog would have to be in order to rule out either model at a ~99.7% confidence level. We find that if the real cosmology is LCDM, a sample of ~200 strong gravitational lenses would be sufficient to rule out R_h=ct at this level of accuracy, while ~300 strong gravitational lenses would be required to rule out LCDM if the real Universe were instead R_h=ct. The difference in required sample size reflects the greater number of free parameters available to fit the data with LCDM. We point out that, should the R_h=ct Universe eventually emerge as the correct cosmology, its lack of any free parameters for this kind of work will provide a remarkably powerful probe of the mass structure in lensing galaxies, and a means of better understanding the origin of the bulge-halo conspiracy.

We present the optical luminosity function (LF) of gamma-ray bursts (GRBs) estimated from a uniform sample of 58 GRBs from observations with the Robotic Optical Transient Search Experiment III (ROTSE-III). Our GRB sample is divided into two sub-samples: detected afterglows (18 GRBs), and those with upper limits (40 GRBs). The $R$ band fluxes 100s after the onset of the burst for these two sub-samples are derived. The optical LFs at 100s are fitted by assuming that the co-moving GRB rate traces the star-formation rate. The detection function of ROTSE-III is taken into account during the fitting of the optical LFs by using Monte Carlo simulations. We find that the cumulative distribution of optical emission at 100s is well-described with an exponential rise and power-law decay (ERPLD), broken power-law (BPL), and Schechter LFs. A single power-law (SPL) LF, on the other hand, is ruled out with high confidence.

The use of time-delay gravitational lenses to examine the cosmological expansion introduces a new standard ruler with which to test theoretical models. The sample suitable for this kind of work now includes 12 lens systems, which have thus far been used solely for optimizing the parameters of $\Lambda$CDM. In this paper, we broaden the base of support for this new, important cosmic probe by using these observations to carry out a one-on-one comparison between {\it competing} models. The currently available sample indicates a likelihood of $\sim 70-80%$ that the $R_{\rm h}=ct$ Universe is the correct cosmology versus $\sim 20-30%$ for the standard model. This possibly interesting result reinforces the need to greatly expand the sample of time-delay lenses, e.g., with the successful implementation of the Dark Energy Survey, the VST ATLAS survey, and the Large Synoptic Survey Telescope. In anticipation of a greatly expanded catalog of time-delay lenses identified with these surveys, we have produced synthetic samples to estimate how large they would have to be in order to rule out either model at a $\sim 99.7%$ confidence level. We find that if the real cosmology is $\Lambda$CDM, a sample of $\sim 150$ time-delay lenses would be sufficient to rule out $R_{\rm h}=ct$ at this level of accuracy, while $\sim 1,000$ time-delay lenses would be required to rule out $\Lambda$CDM if the real Universe is instead $R_{\rm h}=ct$. This difference in required sample size reflects the greater number of free parameters available to fit the data with $\Lambda$CDM.

The principal goal of this paper is to use attempts at reconciling the Swift long gamma-ray bursts (LGRBs) with the star formation history (SFH) to compare the predictions of $\Lambda$CDM with those in the $R_{\rm h}=ct$ Universe. In the context of the former, we confirm that the latest Swift sample of GRBs reveals an increasing evolution in the GRB rate relative to the star formation rate (SFR) at high redshifts. The observed discrepancy between the GRB rate and the SFR may be eliminated by assuming a modest evolution parameterized as $(1+z)^{0.8}$---perhaps indicating a cosmic evolution in metallicity. However, we find a higher metallicity cut of $Z=0.52Z_{\odot}$ than was seen in previous studies, which suggested that LGRBs occur preferentially in metal poor environments, i.e., $Z\sim0.1-0.3Z_{\odot}$. We use a simple power-law approximation to the high-\emph{z} ($\ga 3.8$) SFH, i.e., $R_{\rm SF}\propto[(1+z)/4.8]^{\alpha}$, to examine how the high-\emph{z} SFR may be impacted by a possible abundance evolution in the Swift GRB sample. For an expansion history consistent with $\Lambda$CDM, we find that the Swift redshift and luminosity distributions can be reproduced with reasonable accuracy if $\alpha=-2.41_{-2.09}^{+1.87}$. For the $R_{\rm h}=ct$ Universe, the GRB rate is slightly different from that in $\Lambda$CDM, but also requires an extra evolutionary effect, with a metallicity cut of $Z=0.44Z_{\odot}$. Assuming that the SFR and GRB rate are related via an evolving metallicity, we find that the GRB data constrain the slope of the high-\emph{z} SFR in $R_{\rm h}=ct$ to be $\alpha=-3.60_{-2.45}^{+2.45}$. Both cosmologies fit the GRB/SFR data rather well.

The Palomar Transient Factory (PTF) team recently reported the discovery of a rapidly fading optical transient source, PTF11agg. A long-lived scintillating radio counterpart was identified, but the search for a high energy counterpart showed negative results. The PTF team speculated that PTF11agg may represent a new class of relativistic outbursts. Here we suggest that a neutron star (NS)-NS merger system with a supra-massive magnetar central engine could be a possible source to power such a transient, if our line of sight is not on the jet axis direction of the system. These systems are also top candidates for gravitational wave sources to be detected in the advanced LIGO/Virgo era. We find that the PTF11agg data could be explained well with such a model, suggesting that at least some gravitational wave bursts due to NS-NS mergers may be associated with such a bright electromagnetic counterpart without a \gamma-ray trigger.

We continue our systematic statistical study on optical afterglow data of gamma-ray bursts (GRBs). We present the apparent magnitude distributions of early optical afterglows at different epochs (t= 10^2 s, t = 10^3 s, and 1 hour) for the optical lightcurves of a sample of 93 GRBs (the global sample), and for sub-samples with an afterglow onset bump or a shallow decay segment. For the onset sample and shallow decay sample we also present the brightness distribution at the peak time t_{p} and break time t_{b}, respectively. All the distributions can be fit with Gaussian functions. We further perform Monte Carlo simulations to infer the luminosity function of GRB optical emission at the rest-frame time 10^3 seconds, t_{p}, and t_{b}, respectively. Our results show that a single power-law luminosity function is adequate to model the data, with indices -1.40+/-0.10, -1.06+/- 0.16, and -1.54\+/- 0.22, respectively. Based on the derived rest-frame 10^3 s luminosity function, we generate the intrinsic distribution of the R-band apparent magnitude M_{R} at the observed time 10^{3} seconds post trigger, which peaks at M_{R}=22.5 mag. The fraction of GRBs whose R-band magnitude is fainter than 22 mag, and 25 mag and at the observer time 10^3 seconds are ~63% and ~25%, respectively. The detection probabilities of the optical afterglows with ground-based robotic telescopes and UVOT onboard {Swift} are roughly consistent with that inferred from this intrinsic M_{R} distribution, indicating that the variations of the dark GRB fraction among the samples with different telescopes may be due to the observational selection effect, although the existence of an intrinsically dark GRB population cannot be ruled out.

We continue our systematic statistical study of various components in gamma-ray burst (GRB) optical lightcurves. We decompose the early onset bump and the late re-brightening bump with empirical fits and analyze their statistical properties. Among the 146 GRBs that have well-sampled optical lightcurves, the onset and re-brightening bumps are observed in 38 and 26 GRBs, respectively. It is found that the typical rising and decaying slopes for both the onset and re-brightening bumps are ~1.5 and -1.15, respectively. No early onset bumps in the X-ray band are detected to be associated with the optical onset bumps, while an X-ray re-brightening bump is detected for half of the re-brightening optical bumps. The peak luminosity is anti-correlated with the peak time, L_p\propto t_{p}^{-1.81+/-0.32} and L_p\propto t_{p}^{-0.83+/-0.17} for the onset and re-brightening bumps, respectively. Both L_p and the isotropic energy release of the onset bumps are correlated with E_{\gamma, iso}, whereas no similar correlation is found for the re-brightening bumps. Taking the onset bumps as probes for the properties of the fireballs and their ambient medium, we find that the typical power-law index of the relativistic electrons is 2.5 and the medium density profile behaves as n\propto r^{-1} within the framework of the synchrotron external shock models. With the medium density profile obtained from our analysis, we also confirm the correlation between initial Lorentz factor (\Gamma_0) and E_{\gamma, iso} in our previous work. The jet component that produces the re-brightening bump seems to be on-axis and independent of the prompt emission jet component. Its typical kinetic energy budget would be about one order of magnitude larger than the prompt emission component, but with a lower \Gamma_0, typically several tens.

In this paper, we continue to build support for the proposal to use gamma-ray bursts (GRBs) as standard candles in constructing the Hubble Diagram at redshifts beyond the current reach of Type Ia supernova observations. We confirm that correlations among certain spectral and lightcurve features can indeed be used as luminosity indicators, and demonstrate from the most up-to-date GRB sample appropriate for this work that the $\Lambda$CDM model optimized with these data is characterized by parameter values consistent with those in the concordance model. Specifically, we find that $(\Omega_m,\Omega_\Lambda)\approx (0.25_{-0.06}^{+0.05}, 0.75_{-0.05}^{+0.06})$, which are consistent, to within $1\sigma$, with $(0.29,0.71)$ obtained from the 9-yr WMAP data. We also carry out a comparative analysis between $\Lambda$CDM and the $R_{\rm h}=ct$ Universe and find that the optimal $\Lambda$CDM model fits the GRB Hubble Diagram with a reduced $\chi^2_{\rm dof}\approx 2.26$, whereas the fit using $R_{\rm h}=ct$ results in a $\chi^2_{\rm dof}\approx 2.14$. In both cases, about 20% of the events lie at least $2\sigma$ away from the best-fit curves, suggesting that either some contamination by non-standard GRB luminosities is unavoidable, or that the errors and intrinsic scatter associated with the data are being underestimated. With these optimized fits, we use three statistical tools---the Akaike Information Criterion (AIC), the Kullback Information Criterion (KIC), and the Bayes Information Criterion (BIC)---to show that, based on the GRB Hubble Diagram, the likelihood of $R_{\rm h}=ct$ being closer to the correct model is $\sim 85-96%$, compared to $\sim 4-15%$ for $\Lambda$CDM.

We study the AGN-host co-evolution issue here by focusing on the correlation between the hard X-ray emission from central AGNs and the stellar populations of the host galaxies. By focusing on the galaxies with strong H$\alpha$ line emission (EW(H$\alpha)>5$\AA), both X-ray and optical spectral analysis are performed on 67 (partially) obscured AGNs that are selected from the \it XMM-Newton\rm\ 2XMMi/SDSS-DR7 catalog originally cross-matched by Pineau et al. The sample allows us to study central AGN activity and host galaxy directly and simultaneously in individual objects. Combining the spectral analysis in both bands reveals that the older the stellar population of the host galaxy, the harder the X-ray emission will be, which was missed in our previous study where the \it ROSAT\rm\ hardness ratios are used. By excluding the contamination from the host galaxies and from the jet beaming emission, the correlation indicates that the Compton cooling in the accretion disk corona decreases with the mean age of the stellar population. We argue that the correlation is related to the correlation of $L/L_{\mathrm{Edd}}$ with the host stellar population. In addition, the [\ion{O}{1}]/H$\alpha$ and [\ion{S}{2}]/H$\alpha$ narrow-line ratios are identified to correlate with the spectral slope in hard X-ray, which can be inferred from the currently proposed evolution of the X-ray emission because of the confirmed tight correlations between the two line ratios and stellar population age.

Well-sampled optical lightcurves of 146 gamma-ray bursts (GRBs) are compiled from the literature. By empirical fitting we identify eight possible emission components and summarize the results in a "synthetic" lightcurve. Both optical flare and early shallow-decay components are likely related to long-term central engine activities. We focus on their statistical properties in this paper. Twenty-four optical flares are obtained from 19 GRBs. The isotropic R-band energy is smaller than 1% of $E_{\gamma, \rm iso}$. The relation between isotropic luminosities of the flares and gamma-rays follows $L^{\rm F}_{\rm R, iso}\propto L_{{\gamma}, \rm iso}^{1.11\pm 0.27}$. Later flares tend to be wider and dimmer, i.e., $w^{\rm F}\sim t^{\rm F}_{\rm p}/2$ and $L^{\rm F}_{\rm R, iso}\propto [t^{\rm F}_{\rm p}/(1+z)]^{-1.15\pm0.15}$. The detection probability of the optical flares is much smaller than that of X-ray flares. An optical shallow decay segment is observed in 39 GRBs. The relation between the break time and break luminosity is a power-law, with an index of $-0.78\pm 0.08$, similar to that derived from X-ray flares. The X-ray and optical breaks are usually chromatic, but a tentative correlation is found. We suggest that similar to the prompt optical emission that tracks $\gamma$-rays, the optical flares are also related to the erratic behavior of the central engine. The shallow decay component is likely related to a long-lasting spinning-down central engine or piling up of flare materials onto the blastwave. Mixing of different emission components may be the reason of the diverse chromatic afterglow behaviors.

We present a time-resolved spectral analysis of 51 long and 11 short bright GRBs observed with the {\em Femri}/GBM, paying special attention to $E_{\rm p}$ evolution within a same burst. Among 8 single-pulse long GRBs, 5 show hard-to-soft evolution, while 3 show intensity-tracking. The multi-pulse long GRBs have more complicated patterns. Among the GRBs whose time-resolved spectrum is available for the first pulse, almost half (15/32 GRBs) show clear hard-to-soft evolution, and the other half (17/32 GRBs) show clear intensity-tracking. Later pulses typically show the tracking behavior, although a hard-to-soft evolution pattern was identified in the 2nd pulse of 2 GRBs whose pulses are well separated. Statistically, the hard-to-soft evolution pulses tend to be more asymmetric than the intensity-tracking ones, with a steeper rising wing than the falling wing. Short GRBs have $E_{\rm p}$ tracking intensity exclusively with the 16ms time resolution analysis. We performed a simulation analysis, and suggest that at least for some bursts, the late intensity-tracking pulses could be a consequence of overlapping hard-to-soft pulses. However, the fact that the intensity-tracking pattern exists in the first pulse of multi-pulse long GRBs and some single-pulse GRBs suggest that intensity tracking is an independent component, which may operate in some late pulses as well. For the GRBs with measured redshifts, we present a time-resolved $E_{\rm p}-L_{\gamma, \rm iso}$ correlation analysis and show that the scatter of the correlation is comparable to that of the global Amati/Yonetoku relation. We discuss the predictions of various radiation models regarding $E_{\rm p}$ evolution, as well as the possibility of a precession jet in GRBs. It seems that the data pose great challenge to all these models, and hold the key to unveil the physics of GRB prompt emission.

Context. We present very early optical lightcurves beginning 10 days before maximum of the Type IIn supernova 1998S, covering the first four months after discovery. Aims. We examine the light evolution and try to compare the lightcurves to two analytical models(Nakar & Sari(2010) and Rabinak & Waxman(2011)) for a red supergiant star. Methods. The photometry was carried at the 60-cm telescope of the Xinglong Station of China. Broadband filters Johnsons B, V and Cousins R were used. Results. The magnitude rose for the first few days and then dropped slowly afterwards. The two different models we use can fit the early lightcurves very well. The explosion date derived from the models is within the range of 1998 March 1.34 - 2.64(JD 2450873.84 - JD 2450875.14.) The radius of the progenitor is found to be ~ 300 Rsun and ~ 2000 Rsun for the model of Nakar & Sari(2010) and Rabinak & Waxman(2011) respectively. The constraint on mass and energy is not strong. The ranges of these two parameters are within that of a red supergiant.

We present the study on GRB 071112C X-ray and optical light curves. In these two wavelength ranges, we have found different temporal properties. The R-band light curve showed an initial rise followed by a single power-law decay, while the X-ray light curve was described by a single power-law decay plus a flare-like feature. Our analysis shows that the observed temporal evolution cannot be described by the external shock model in which the X-ray and optical emission are produced by the same emission mechanism. No significant color changes in multi-band light curves and a reasonable value of the initial Lorentz factor ({\Gamma}0 = 275 \pm 20) in a uniform ISM support the afterglow onset scenario as the correct interpretation for the early R-band rise. The result suggests the optical flux is dominated by afterglow. Our further investigations show that the X-ray flux could be created by an additional feature related to energy injection and X-ray afterglow. Different theoretical interpretations indicate the additional feature in X-ray can be explained by either late internal dissipation or local inverse-Compton scattering in the external shock.

The energy injection model is usually proposed to interpret the shallow-decay phase in Swift GRB X-ray afterglows. However, very few GRBs have simultaneous signatures of energy injection in their optical and X-ray afterglows. Here, we report optical observations of GRB 090529A from 2000 sec to $\sim10^6$ sec after the burst, in which an achromatic decay is seen at both wavelengths. The optical light curve shows a decay from 0.37 to 0.99 with a break at $\sim10^5$ sec. In the same time interval, the decay indices of the X-ray light curve changed from 0.04 to 1.2. Comparing these values with the closure relations, the segment after 3$\times10^{4}$ sec is consistent with the prediction of the forward shock in an ISM medium without any energy injection. The shallow-decay phase between 2000 to 3$\times10^{4}$ sec could be due to the external shock in a wind-type-like medium with an energy injection under the condition of $\nu_o < \nu_c < \nu_x$. However, the constraint of the spectral region is not well consistent with the multi-band observations. For this shallow-decay phase, other models are also possible, such as energy injection with evolving microphysical parameters, or a jet viewed off-axis,etc.

Apparent redshift dependence of the jet opening angles ($\theta_{\rm j}$) of gamma-ray bursts (GRBs) is observed from current GRB sample. We investigate whether this dependence can be explained with instrumental selection effects and observational biases by a bootstrapping method. Assuming that (1) the GRB rate follows the star formation history and the cosmic metallicity history and (2) the intrinsic distributions of the jet-corrected luminosity ($L_{\rm \gamma}$) and $\theta_{\rm j}$ are a Gaussian or a power-law function, we generate a mock {\em Swift}/BAT sample by considering various instrumental selection effects, including the flux threshold and the trigger probability of BAT, the probabilities of a GRB jet pointing to the instrument solid angle and the probability of redshift measurement. Our results well reproduce the observed $\theta_{\rm j}-z$ dependence. We find that in case of $L_{\gamma}\propto \theta_{\rm j}^2$ good consistency between the mock and observed samples can be obtained, indicating that both $L_{\rm \gamma}$ and $\theta_{\rm j}$ are degenerate for a flux-limited sample. The parameter set $(L_{\rm \gamma}, \theta_{\rm j})=(4.9\times 10^{49} \rm {erg\ s}^{-1},\ 0.054 {rad})$ gives the best consistency for the current {\em Swift} GRB sample. Considering the beaming effect, the derived intrinsic local GRB rate accordingly is $2.85\times 10^2$ Gpc$^{-3}$ yr$^{-1}$, inferring that $\sim 0.59%$ of Type Ib/c SNe may be accompanied by a GRB.

The AGN-host co-evolution issue is investigated here by focusing on the evolution of the [\ion{O}{3}]$\lambda5007$ emission-line profile. In order to simultaneously measure both [\ion{O}{3}] line profile and circumnuclear stellar population in individual spectrum, a large sample of narrow emission-line galaxies is selected from the MPA/JHU SDSS DR7 catalog. By requiring that 1) the [\ion{O}{3}] line signal-to-noise ratio is larger than 30, 2) the [\ion{O}{3}] line width is larger than the instrumental resolution by a factor of 2, our sample finally contains 2,333 Seyfert galaxies/LINERs (AGNs), 793 transition galaxies, and 190 starforming galaxies. In additional to the commonly used profile parameters (i.e., line centroid, relative velocity shift and velocity dispersion), two dimensionless shape parameters, skewness and kurtosis, are used to quantify the line shape deviation from a pure Gaussian function. We show that the transition galaxies are systematically associated with narrower line widths and weaker [\ion{O}{3}] broad wings than the AGNs, which implies that the kinematics of the emission-line gas is different in the two kinds of objects. By combining the measured host properties and line shape parameters, we find that the AGNs with stronger blue asymmetries tend to be associated with younger stellar populations. However, the similar trend is not identified in the transition galaxies. The failure is likely resulted from a selection effect in which the transition galaxies are systematically associated with younger stellar populations than the AGNs. The evolutionary significance revealed here suggests that both NLR kinematics and outflow feedback in AGNs co-evolve with their host galaxies.

The potential effect of the lunar exosphere on the near-ultraviolet sky background emission is predicted for Lunar-based Ultraviolet Telescope (LUT: a funded Chinese scientific payload for the Chang'e-III mission). Using the upper limit on the OH concentration inferred from the recent MIP CHACE results, our calculations show that the sky brightness due to the illuminated exosphere is $<8.7\ \mathrm{photon\ s^{-1}\ cm^{-2}\ arcsec^{-2}}$ within the wavelength range 245-340 nm. By evaluating the signal-to-noise ratios of observations of an AB=13 mag point source at a series of sky background levels, our analysis indicates that the detection performance of LUT can be moderately degraded by the lunar exosphere emission in most cases. An AB=13 mag point source can still be detected by the telescope at a signal-to-noise ratio more than 8 when the OH concentration is less than $2\times10^8\ \mathrm{molecules\ cm^{-3}}$. However, the effect on the performance is considerable when the exosphere is as dense as suggested by CHACE.

Fractals are a basic tool to phenomenologically describe natural objects having a high degree of temporal or spatial variability. From a physical point of view the fractal properties of natural systems can also be interpreted by using the standard formalism of thermodynamical fluctuations. In the present paper we introduce and analyze the fractal dimension of the Intra - Day- Variability (IDV) light-curves of the BL Lac objects, in the optical, radio and X-ray bands, respectively. A general description of the fluctuation spectrum of these systems based on general thermodynamical principles is also proposed. Based on the general fractal properties of a given physical system, we also introduce the predictability index for the IDV light curves. We have explicitly determined the fractal dimension for the R-band observations of five blazars, as well as for the radio band observations of the compact extragalactic radio source J 1128+5925, and of several X-ray sources. Our results show that the fractal dimension of the optical and X-ray observations indicates an almost pure "Brownian noise" (random walk) spectrum, with a very low predictability index, while in the radio band the predictability index is much higher. We have also studied the spectral properties of the IDV light curves, and we have shown that their spectral index is very closely correlated with the corresponding fractal dimension.

We present the SVOM (Space-based multi-band astronomical Variable Objects Monitor) mission that the Chinese National Space Agency and the French Space Agency have decided to jointly implement. SVOM has been designed to detect all known types of gamma-ray bursts (GRBs), to provide fast and reliable GRB positions, to measure the broadband spectral shape and temporal properties of the GRB prompt emission, and to quickly identify the optical/near-infrared afterglows of detected GRBs, including high-redshift ones. Scheduled to be in orbit in the second half of the present decade, the SVOM satellite will carry a very innovative scientific payload combining for the first time a wide field X- and gamma-ray coded mask imager for GRB real-time localizations to few arcmin, a non-imaging gamma-ray monitor, and two narrow-field instruments for the study of the GRB early afterglow emission in the X-ray and visible bands. The SVOM payload is complemented by ground-based instruments including a wide-field camera to catch the GRB prompt emission in the visible band and two robotic telescopes to measure the photometric properties of the early afterglow. A particular attention is paid to the GRB follow-up in facilitating the observation of the SVOM detected GRB by the largest ground based telescopes.

An NUV-optical diagram made for sources from the secend Galaxy Evolution Explorer (GALEX) Ultraviolet Variability (GUVV-2) Catalog provide us a method to tentatively classify the unknown GUVV2 sources by their NUV-optical magnitudes. On the purpose of testing the correctness and generality of the method, we carry out a program on the spectroscopic observations of the unidentified GUVV2 sources. The spectroscopic identification for these 37 sources are 19 type -A to -F stars, 10 type -G to -K stars and 7 M dwarf stars together with an AGN. We also present the light curves in R-band for two RR Lyrae star candidates selected from the NUV-optical diagram, both of which perform cyclic variations. Combining there light curves and colors, we classify them as RR Lyrae stars. To confirm the results, we shows a color-color diagram for the 37 newly spectroscopically identified objects compared with the previously identified ones, which manifests good consistence with our former results, indicating that the ultroviolet variable sources can be initially classified by their NUV/optical color-color diagram.

In this paper, we introduce a simplified model for explaining the observations of the optical intraday variability (IDV) of the BL Lac Objects. We assume that the source of the IDV are the stochastic oscillations of an accretion disk around a supermassive black hole. The Stochastic Fluctuations on the vertical direction of the accretion disk are described by using a Langevin type equation with a damping term and a random, white noise type force. Furthermore, the preliminary numerical simulation results are presented, which are based on the numerical analysis of the Langevin stochastic differential equation.

The spectral properties from radio to optical bands are compared between the 18 optically bright Gamma-ray burst afterglows and well established power-spectrum sequence in Blazars. The comparison shows that the afterglows are well agreement with the well known Blazar sequence (i.e., the $\nu L_{\nu}(\mathrm{5GHz})$-$\alpha_{\mathrm{RO}}$ correlation, where $\alpha_{\mathrm{RO}}$ is the broad-band spectral slope from radio to optical bands). The afterglows are, however, clustered at the low luminosity end of the sequence, which is typically occupied by high frequency-peaked BL Lac objects. The correlation suggests that Gamma-ray burst afterglows share the similar emission process with high frequency-peaked BL Lac objects. We further identify a deviation at a significance level larger than 2$\sigma$ from the sequence for three typical optically "dark" bursts. The deviation favors a heavy extinction in optical bands for the "dark" bursts. The extinction $A_V$ is estimated to be larger than 0.5-0.6 magnitude from the $\nu L_{\nu}(\mathrm{5GHz})$-$\alpha_{\mathrm{RO}}$ sequence.

GRB 090426 is a short duration burst detected by Swift ($T_{90}\sim 1.28$ s in the observer frame, and $T_{90}\sim 0.33$ s in the burst frame at $z=2.609$). Its host galaxy properties and some $\gamma$-ray related correlations are analogous to those seen in long duration GRBs, which are believed to be of a massive-star origin (so-called Type II GRBs). We present the results of its early optical observations with the 0.8-m TNT telescope at Xinglong observatory, and the 1-m LOAO telescope at Mt. Lemmon Optical Astronomy Observatory in Arizona. Our well-sampled optical afterglow lightcurve covers from $\sim 90$ seconds to $\sim 10^4$ seconds post the GRB trigger. It shows two shallow decay episodes that are likely due to energy injection, which end at $\sim 230$ seconds and $\sim 7100$ seconds, respectively. The decay slopes post the injection phases are consistent with each other ($\alpha\simeq 1.22$). The X-ray afterglow lightcurve appears to trace the optical, although the second energy injection phase was missed due to visibility constraints introduced by the {\em Swift} orbit. The X-ray spectral index is $\beta_X\sim 1.0$ without temporal evolution. Its decay slope is consistent with the prediction of the forward shock model. Both X-ray and optical emission is consistent with being in the same spectral regime above the cooling frequency ($\nu_c$). The fact that $\nu_c$ is below the optical band from the very early epoch of the observation provides a constraint on the burst environment, which is similar to that seen in classical long duration GRBs. We therefore suggest that death of a massive star is the possible progenitor of this short burst.

As the third paper of our serial studies that are aim at examining the AGN-host coevolution by using partially obscured AGNs, we extend the broad-line composite galaxies (composite AGNs) into ROSAT-selected Seyfert 1.8/1.9 galaxies basing upon the RASS/SDSS-DR5 catalog given by Anderson et al.. The SDSS spectra of in total 92 objects are analyzed by the same method used in our previous studies, after requiring the signal-to-noise ratio in the SDSS r' band is larger than 20. Combing the ROSAT-selected Seyfert galaxies with the composite AGNs reinforces the tight correlation between the line ratio [OI]/H\alpha vs. D_n(4000), and establishes a new tight correlation between [SII]/H\alpha vs. D_n(4000). Both correlations suggest the two line ratios are plausible age indicators of the circumnuclear stellar population for typical type I AGNs in which the stellar populations are difficult to be derived from their optical spectra. The ROSAT-selected Seyfert galaxies show that the two correlations depend on the soft X-ray spectral slope \alpha_X that is roughly estimated from the hardness ratios by requiring the X-ray count rates within 0.1-2.4 keV are larger than 0.02 counts s^-1. However, we fail to establish a relationship between \alpha_X and D_n(4000), which is likely caused by the relatively large uncertainties of both parameters (especially for \alpha_X because of the AGN intrinsic obscuration). The previously established L/L_Edd-D_n(4000) evolutionary sequence is reinforced again by the extension to the ROSAT-selected Seyfert galaxies. These X-ray-selected Seyfert galaxies are, however, biased against the two ends of the sequence, which implies that the X-ray Seyfert galaxies present a population at middle evolutionary stage.

Aims. We have gathered optical spectra of 8 long-duration GRB host galaxies selected from the archival data of VLT/FORS2. We investigated whether or not Wolf-Rayet (WR) stars can be detected in these GRB host galaxies. We also tried to estimate the physical properties of GRB host galaxies, such as metallicity. Methods. We identified the WR features in these spectra by fitting the WR bumps and WR emission lines in blue and red bumps. We also identified the subtypes of the WR stars, and estimated the numbers of stars in each subtype, then calculated the WR/O star ratios. The (O/H) abundances of GRB hosts were estimated from both the electron temperature (Te) and the metallicity-sensitive strong-line ratio (R23), for which we have broken the R23 degeneracy. We compared the environments of long-duration GRB host galaxies with those of other galaxies in terms of their luminosity (stellar mass)-metallicity relations (LZ, MZ). Results. We detected the presence of WR stars in 5 GRB host galaxies having spectra with relatively high signal-to-noise ratios (S/N). In the comparison of LZ, MZ relations, it shows that GRB hosts have lower metallicities than other samples with comparable luminosity and stellar mass. The presence of WR stars and the observed high WR/O star ratio, together with low metallicity, support the core-collapsar model and implie the first stage of star formation in the hosted regions of GRBs.

We present the Space-based multi-band astronomical Variable Object Monitor (SVOM), a future satellite mission for Gamma-Ray Burst (GRB) studies, developed in cooperation between the Chinese National Space Agency (CNSA), the Chinese Academy of Science (CAS), the French Space Agency (CNES) and French research institutes. The scientific objectives of the SVOM GRB studies cover their classification (GRB diversity and unity of the model), their physics (particle acceleration and radiation mechanisms), their progenitors, cosmological studies (host galaxies, star formation history, re-ionization, cosmological parameters), and fundamental physics (origin of cosmic rays, Lorentz invariance, gravitational wave sources). From 2015 on, SVOM will provide fast and accurate localizations of all known types of GRB, and determine the temporal and spectral properties of the GRB emission, thanks to a set of four onboard instruments. The trigger system of the coded-mask telescope ECLAIRs onboard SVOM images the sky in the 4-120 keV energy range, in order to detect and localize GRB in its 2 sr-wide field of view. The low-energy threshold of ECLAIRs is well suited for the detection of highly red-shifted GRB. The high-energy coverage is extended up to 5 MeV thanks to a non-imaging gamma-ray spectrometer. GRB alerts are sent in real-time to the ground observers community, and a spacecraft slew is performed in order to place the GRB within the narrow fields of view of a soft X-ray telescope and a visible-band telescope, to refine the GRB position and study its early afterglow. Ground-based robotic telescopes and wide-angle cameras complement the onboard instruments. A large fraction of GRB will have redshift determinations, thanks to an observing strategy optimized to facilitate follow-up observations by large ground-based spectroscopic telescopes.

SDSS J094857.3+002225 is a very radio-loud narrow-line Seyfert 1 (NLS1) galaxy. Here, we report our discovery of the intranight optical variability (INOV) of this galaxy through the optical monitoring in the B and R bands that covered seven nights in 2009. Violent rapid variability in the optical bands was identified in this RL-NLS1 for the first time, and the amplitudes of the INOV reaches 0.5 mag in both the B and R bands on the timescale of several hours. The detection of the INOV provides a piece of strong evidence supporting the fact that the object carries a relativistic jet with a small viewing angle, which confirms the conclusion drawn from the previous multi-wavelength studies.

The correlation between distant Gamma-Ray Bursts (GRBs) and foreground galaxy clusters is re-examined by using the well localized (with an accuracy down to a few arcseconds) Swift/XRT GRBs. The galaxy clusters are compiled from both X-ray selected ROSAT brightest cluster sample (BCS) and BCS extension by requiring $\delta \geq0\degr$ and $b\geq20\degr$. The Swift/XRT GRBs fulfilling the above selection criteria are cross-correlated with the clusters. Both Nearest-Neighbor Analysis and angular two-point cross-correlation function show that there is no enough evidence supporting the correlation between the GRBs and foreground clusters. We suggest that the non-correlation is probably related to the GRB number-flux relation slope.

Since the launch of Swift satellite, the detections of high-z (z>4) long gamma-ray bursts (LGRBs) have been rapidly growing, even approaching the very early Universe (the record holder currently is z=8.3). The observed high-z LGRB rate shows significant excess over that estimated from the star formation history. We investigate what may be responsible for this high productivity of GRBs at high-z through Monte Carlo simulations, with effective Swif/BAT trigger and redshift detection probabilities based on current Swift/BAT sample and CGRO/BATSE LGRB sample. We compare our simulations to the Swift observations via log N-log P, peak luminosity (L) and redshift distributions. In the case that LGRB rate is purely proportional to the star formation rate (SFR), our simulations poorly reproduce the LGRB rate at z>4, although the simulated log N-log P distribution is in good agreement with the observed one. Assuming that the excess of high-z GRB rate is due to the cosmic metallicity evolution or unknown LGRB rate increase parameterized as (1+z)^delta, we find that although the two scenarios alone can improve the consistency between our simulations and observations, incorporation of them gives much better consistency. We get 0.2<epsilon<0.6 and delta<0.6, where epsilon is the metallicity threshold for the production of LGRBs. The best consistency is obtained with a parameter set (epsilon, delta)=(~0.4, ~0.4), and BAT might trigger a few LGRBs at z~14. With increasing detections of GRBs at z>4 (~15% of GRBs in current Swift LGRB sample based on our simulations), a window for very early Universe is opening by Swift and up-coming SVOM missions.

The normalized radial distribution of young stellar populations (and cold gas) in nearby galactic disks is compared between AGN host galaxies and starforming galaxies (both with Hubble types between S0/a and Scd) by using type II supernovae (SNe) as tracers. A subset of 140 SNe\,II with available supernova position measurements are selected from the SAI-SDSS image catalog by requiring available SDSS spectroscopy data of their host galaxies. Our sample is finally composed of 46 AGNs and 94 starforming galaxies. Both directly measured number distributions and inferred surface density distributions indicate that a) the SNe detected in starforming galaxies follow an exponential law well; b) by contrast, the SNe detected in AGN host galaxies significantly deviate from an exponential law, which is independent of both morphological type and redshift. Specifically, we find a detection deficit around $R_{\mathrm{SN}}/R_{25,\mathrm{cor}}\sim0.5$ and an over-detection at outer region $R_{\mathrm{SN}}/R_{25,\mathrm{cor}}\sim0.6-0.8$. This finding provides a piece of evidence supporting that there is a link between ongoing star formation (and cold gas reservoir) taking place in the extended disk and central AGN activity.

1RXS J201607.0251645 is identified to be an eclipsing binary for the first time. We present the preliminary observations in V band with the 0.6-m telescope for three years and the extensive observations in V and R band with the 0.8-m telescope for six nights respectively. The light curve of the system is EB type. Five light minimum times were obtained and the orbital period of 0.388058d(0.00044d) is determined. The photometric solution given by the 2003-version Wilson-Devinney program suggests the binary is a semi-detached system with the photometric mass ratio 0.895(0.006), which probably comprise a G5 primary and an oversize K5 secondary. The less massive component has completely filled up its Roche lobe, while the other almost fills its Roche lobe with the filling factor of 93.4%. The system shows a varying O'Connell effect in its phase folded diagrams from 2005 to 2007, and is X-ray luminous with log(LXLbol)3.27. Possible mechanisms to account for these two phenomena are discussed. Finally, we infer the binary may be in thermal oscillation predicted by thermal relaxation oscillation (TRO) theory or may evolve into a contact binary.

We present our optical observations of {\em Swift} GRB 070518 afterglow obtained at the 0.8-m Tsinghua University-National Astronomical Observatory of China telescope (TNT) at Xinglong Observatory. Our follow-up observations were performed from 512 sec after the burst trigger. With the upper limit of redshift $\sim$0.7, GRB 070518 is found to be an optically dim burst. The spectra indices $\beta_{ox}$ of optical to X-ray are slightly larger than 0.5, which implies the burst might be a dark burst. The extinction $A_{V}$ of the host galaxy is 3.2 mag inferred from the X-ray hydrogen column density with Galactic extinction law, and 0.3 mag with SMC extinction law. Also, it is similar to three other low-redshift optically dim bursts, which belong to XRR or XRF, and mid-term duration($T_{90}<10$, except for GRB 070419A, $T_{90}$=116s). Moreover, its $R$ band afterglow flux is well fitted by a single power-law with an index of 0.87. The optical afterglow and the X-ray afterglow in the normal segment might have the same mechanism, as they are consistent with the prediction of the classical external shock model. Besides, GRB 070518 agrees with Amati relation under reasonable assumptions. The Ghirlanda relation is also tested with the burst.

We present multi-band results for GRB071010B based on \Swift, \Suzaku, and ground-based optical observations. This burst is an ideal target to evaluate the robustness of the ${\rm E^{src}_{peak}-E_{iso}}$ and ${\rm E^{src}_{peak}-E_{\gamma}}$ relations, whose studies have been in stagnation due to the lack of the combined estimation of $\rm E^{src}_{peak}$ and long term optical monitoring. The joint prompt spectral fitting using \Swift/Burst Alert Telescope and \Suzaku/Wide-band All sky Monitor data yielded the spectral peak energy as E$^{src}_{peak}$ of $86.5^{+6.4}_{-6.3}$ keV and E$_{iso}$ of $2.25^{+0.19}_{-0.16}\times10^{52}$ erg with $z=0.947$. The optical afterglow light curve is well fitted by a simple power law with temporal index $\alpha=-0.60\pm0.02$. The lower limit of temporal break in the optical light curve is 9.8 days. Our multi-wavelength analysis reveals that GRB071010B follows ${\rm E^{src}_{peak}-E_{iso}}$ but violates the ${\rm E^{src}_{peak}-E_{\gamma}}$ and ${\rm E_{iso}-E^{src}_{peak}-t^{src}_{jet}}$ at more than the 3$\sigma$ level.

We present the SVOM (Space-based multi-band astronomical Variable Object Monitor) mission, that is being developed in cooperation between the Chinese National Space Agency (CNSA), the Chinese Academy of Science (CAS) and the French Space Agency (CNES). Its scientific objectives include the study of the GRB phenomenon, GRB physics and progenitors, cosmology, and fundamental physics. SVOM is designed to detect all known types of Gamma-Ray Bursts (GRBs), to provide fast and reliable GRB positions, to measure the broadband spectral characteristics and temporal properties of the GRB prompt emission. This will be obtained in first place thanks to a set of four space flown instruments. A wide field (~2 sr) coded mask telescope (ECLAIRs), operating in the 4-250 keV energy range, will provide the triggers and localizations, while a gamma-ray non-imaging spectrometer (GRM), sensitive in the 50 keV-5 MeV domain, will extend the prompt emission energy coverage. After a satellite slew, in order to place the GRB direction within field of view of the two narrow field instruments - a soft X-ray (XIAO), and a visible telescope (VT) - the GRB position will be refined and the study of the early phases of the GRB afterglow will be possible. A set of three ground based dedicated instruments, two robotic telescopes (GFTs) and a wide angle optical monitor (GWAC), will complement the space borne instruments. Thanks to the low energy trigger threshold (~4 keV) of the ECLAIRs, SVOM is ideally suited for the detection of soft, hence potentially most distant, GRBs. Its observing strategy is optimized to facilitate follow-up observations from the largest ground based facilities.

A gamma-ray burst (GRB) optical photometric follow-up system at the Xinglong Observatory of National Astronomical Observatories of China (NAOC) has been constructed. It uses the 0.8-m Tsinghua-NAOC Telescope (TNT) and the 1-m EST telescope, and can automatically respond to GRB Coordinates Network (GCN)alerts. Both telescopes slew relatively fast, being able to point to a new target field within about 1 min upon a request. Whenever available, the 2.16-m NAOC telescope is also used. In 2006, the system responded to 15 GRBs and detected seven early afterglows. In 2007, six GRBs have been detected among 18 follow-up observations. TNT observations of the second most distant GRB 060927 (z=5.5) are shown, which started as early as 91s after the GRB trigger. The afterglow was detected in the combined image of first 19x20s unfiltered exposures. This GRB follow-up system has joined the East-Asia GRB Follow-up Observation Network (EAFON).

We present our analysis of the three HII-buried-AGN: SDSS J091053+333008, SDSS J121837+091324, and SDSS J153002-020415, by studying their optical spectra extracted from SDSS. The location in the BPT diagnostic diagrams of the three galaxies indicates that the narrow emission lines are mainly exited from HII regions. However, after the removal of the host galaxy's stellar emission, the emission lines display the typical feature of Narrow-line Seyfert 1-like. All of the three objects have large Eddington ratio, small black hole mass, and low star formation rate. We propose that the three galaxies are at the transit stage from the starburst-dominated phase to AGN-dominated phase.

We introduce a new physical parameter, the optical variability amplitude, to the well-established Eigenvector 1 space of quasars and test a sample of long-term B-band light curves of 42 PG quasars monitored by Giveon et al. (1999). We find that the optical variability amplitude strongly correlates with the intensity ratio of Fe II to H$\beta$, H$\beta$ width and peak luminosity at 5007$\rm{\AA}$. We briefly discuss the physical meaning of our findings and suggest that the Eddington ratio may be a key factor in determining a quasar's variability.

We study the issue of AGN-host connection in intermediate-z (1.2>z>0.4) galaxies with hybrid spectra (hybrid QSOs for short). The observed spectra redward of the Balmer limit are dominated by starlight, and the spectra at blue end by both an AGN continuum and a MgII broad emission line. Such unique property allows us to examine both AGN and its host galaxy in individual galaxy simultaneously. At first, 15 hybrid QSOs are selected from the Sloan Digital Sky Survey Data Release 6. The spectra are then analyzed in detail in three objects: SDSS J162446.49+461946.7, SDSS J102633.32+103443.8 and SDSS J090036.44+381353.0. Our spectral analyzing shows that the current star formation activities are strongly suppressed, and that the latest burst ages range from ~400Myr to 1Gyr. Basing upon the MgII-based black hole masses, the three hybrid QSOs are consistent with the D_n(4000)-L/L_{Edd} sequence that was previously established in local AGNs. The three hybrid QSOs are located in the middle range of the sequence, which implies that the hybrid QSOs are at the transition stage not only from young to old AGN, but also from host-dominated phase to AGN-dominated phase.

The nature of Low-ionization Nuclear Emission-line Regions (LINERs) has been an open question for a long time. We study the properties of LINERs from several different aspects. The LINERs are found to consist of two different categories that can be clearly separated in the traditional BPT diagrams, especially in the [OI]/Ha vs. [OIII]/Hb diagram. LINERs with high [O]/Ha ratios (population I) differ from ones with low ratios (population II) in several properties. Broad emission lines are only identified in the spectra of population I LINERs. While only the population II LINERs show luminous infrared emission and occurrence of core-collapse supernovae in the host. Combining these results with the known distribution of stellar populations not only suggests that the two populations have different line excitation mechanisms, but also implies that they are at different evolutionary stages.

We examine the accretion properties in a sample of 42 hard (3-60keV) X-ray selected nearby broad-line AGNs. The energy range in the sample is harder than that usually used in the similar previous studies. These AGNs are mainly complied from the RXTE All Sky Survey (XSS), and complemented by the released INTEGRAL AGN catalog. The black hole masses, bolometric luminosities of AGN, and Eddington ratios are derived from their optical spectra in terms of the broad H$\beta$ emission line. The tight correlation between the hard X-ray (3-20keV) and bolometric/line luminosity is well identified in our sample. Also identified is a strong inverse Baldwin relationship of the H$\beta$ emission line. In addition, all these hard X-ray AGNs are biased toward luminous objects with high Eddington ratio (mostly between 0.01 to 0.1) and low column density ($<10^{22} \mathrm{cm^{-2}}$), which is most likely due to the selection effect of the surveys. The hard X-ray luminosity is consequently found to be strongly correlated with the black hole mass. We believe the sample completeness will be improved in the next few years by the ongoing Swift and INTEGRAL missions, and by the next advanced missions, such as NuSTAR, Simbol-X, and NeXT. Finally, the correlation between RFe (=optical FeII/H$\beta$) and disk temperature as assessed by $T\propto (L/L_{\mathrm{Edd}})M_{\mathrm{BH}}^{-1}$ leads us to suggest that the strength of the FeII emission is mainly determined by the shape of the ionizing spectrum.

We present the Space-based multi-band astronomical Variable Objects Monitor mission (SVOM) decided by the Chinese National Space Agency (CNSA) and the French Space Agency (CNES). The mission which is designed to detect about 80 Gamma-Ray Bursts (GRBs) of all known types per year, will carry a very innovative scientific payload combining a gamma-ray coded mask imagers sensitive in the range 4 keV to 250 keV, a soft X-ray telescope operating between 0.5 to 2 keV, a gamma-ray spectro-photometer sensitive in the range 50 keV to 5 MeV, and an optical telescope able to measure the GRB afterglow emission down to a magnitude limit M$_R=23$ with a 300 s exposure. A particular attention will be also paid to the follow-up in making easy the observation of the SVOM detected GRB by the largest ground based telescopes. Scheduled for a launch in 2013, it will provide fast and reliable GRB positions, will measure the broadband spectral energy distribution and temporal properties of the prompt emission, and will quickly identify the optical afterglows of detected GRBs, including those at very high redshift.

We report on multi-wavelength observations, ranging from the X-ray to radio wave bands, of the IPN-localized gamma-ray burst GRB 070125. Spectroscopic observations reveal the presence of absorption lines due to O I, Si II, and C IV, implying a likely redshift of z = 1.547. The well-sampled light curves, in particular from 0.5 to 4 days after the burst, suggest a jet break at 3.7 days, corresponding to a jet opening angle of ~7.0 degrees, and implying an intrinsic GRB energy in the 1 - 10,000 keV band of around E = (6.3 - 6.9)x 10^(51) erg (based on the fluences measured by the gamma-ray detectors of the IPN network). GRB 070125 is among the brightest afterglows observed to date. The spectral energy distribution implies a host extinction of Av < 0.9 mag. Two rebrightening episodes are observed, one with excellent time coverage, showing an increase in flux of 56% in ~8000 seconds. The evolution of the afterglow light curve is achromatic at all times. Late-time observations of the afterglow do not show evidence for emission from an underlying host galaxy or supernova. Any host galaxy would be subluminous, consistent with current GRB host-galaxy samples. Evidence for strong Mg II absorption features is not found, which is perhaps surprising in view of the relatively high redshift of this burst and the high likelihood for such features along GRB-selected lines of sight.

The goal of our serial papers is to examine the evolutionary connection between AGN and star formation in its host galaxy in the partially obscured AGNs (i.e., Seyfert 1.8 and 1.9 galaxies). Taking advantage of these galaxies, the properties of both components can be studied together by direct measurements. In this paper, we focus on the broad-line composite galaxies (composite AGNs) which are located between the theoretical and empirical separation lines in the [NII]/Ha vs. [OIII]/Hb diagram. These galaxies are searched for from the composite galaxies provided by the SDSS DR4 MPA/JHU catalogs. After re-analyze the spectra, we perform a fine classification for the 85 composite AGNs in terms of the BPT diagrams. All the objects located below the three theoretical separation lines are associated with a young stellar population (<1Gyrs), while either a young or old stellar population is identified in the individual multiply-classified object. The multiply-classified objects with a very old stellar population are located in the LINER region in the [OI]/Ha vs. [OIII]/Hb diagram. We then consider the connection between AGN and star formation to derive the key results. The Eddington ratio inferred from the broad Ha emission, the age of the stellar population of AGN's host as assessed by D_n(4000), and the line ratio [OI]/Ha are found to be related with each other. These relations strongly support the evolutionary scenario in which AGNs evolve from high L/L_Edd state with soft spectrum to low L/L_Edd state with hard spectrum as young stellar population ages and fades. The significant correlation between the line ratio [OI]/Ha and D_n(4000) leads us to suggest that the line ratio could be used to trace the age of stellar population in type I AGNs.

The background-limited spectral imaging of the early Universe requires spaceborne terahertz (THz) detectors with the sensitivity 2-3 orders of magnitude better than that of the state-of-the-art bolometers. To realize this sensitivity without sacrificing operating speed, novel detector designs should combine an ultrasmall heat capacity of a sensor with its unique thermal isolation. Quantum effects in thermal transport at nanoscale put strong limitations on the further improvement of traditional membrane-supported bolometers. Here we demonstrate an innovative approach by developing superconducting hot-electron nanobolometers in which the electrons are cooled only due to a weak electron-phonon interaction. At T<0.1K, the electron-phonon thermal conductance in these nanodevices becomes less than one percent of the quantum of thermal conductance. The hot-electron nanobolometers, sufficiently sensitive for registering single THz photons, are very promising for submillimeter astronomy and other applications based on quantum calorimetry and photon counting.

We study the ``normal'' decay phase of the X-ray afterglows of gamma-ray bursts (GRBs), which follows the shallow decay phase, using the events simultaneously observed in the R-band. The classical external shock model -- in which neither the delayed energy injection nor time-dependency of shock micro-physics is considered -- shows that the decay indices of the X-ray and R-band light curves, $\alpha_{\rm X}$ and $\alpha_{\rm O}$, obey a certain relation, and that in particular, $\alpha_{\rm O}-\alpha_{\rm X}$ should be larger than -1/4 unless the ambient density increases with the distance from the central engine. For our selected 14 samples, we have found that 4 events violate the limit at more than the 3$\sigma$ level, so that a fraction of events are outliers of the classical external shock model at the ``normal'' decay phase.

The nebular spectra of the broad-lined, SN 1998bw-like Type Ic SN 2002ap are studied by means of synthetic spectra. Two different modelling techniques are employed. In one technique, the SN ejecta are treated as a single zone, while in the other a density and abundance distribution in velocity is used from an explosion model. In both cases, heating caused by gamma-ray and positron deposition is computed (in the latter case using a Monte Carlo technique to describe the propagation of gamma-rays and positrons), as is cooling via forbidden-line emission. The results are compared, and although general agreement is found, the stratified models are shown to reproduce the observed line profiles much more accurately than the single-zone model. The explosion produced ~ 0.1 Msun of 56Ni. The distribution in velocity of the various elements is in agreement with that obtained from the early-time models, which indicated an ejected mass of ~ 2.5 Msun with a kinetic energy of 4 x 10^{51} erg. Nebular spectroscopy confirms that most of the ejected mass (~ 1.2 Msun) was oxygen. The presence of an oxygen-rich inner core, combined with that of 56Ni at high velocities as deduced from early-time models, suggests that the explosion was asymmetric, especially in the inner part.

The properties of narrow-line Seyfert 1 (NLS1) galaxies, the links and correlations between them, and the physics behind them, are still not well understood. Apart from accretion rates and black hole masses, density and outflows were speculated to be among the main drivers of the NLS1 phenomenon. Here, we utilize the diagnostic power of the [SII]6716,6731 intensity ratio to measure the density of the NLR systematically and homogeneously for a large sample of NLS1 galaxies, and we perform a comparison with a sample of broad-line type 1 AGN. We report the discovery of a zone of avoidance in density in the sense that AGN with broad lines (FWHM_Hbeta > 2000 km/s) avoid low densities, while NLS1 galaxies show a wider distribution in the NLR density, including a significant number of objects with low densities. A correlation analysis further shows that the Eddington ratio L/L_Edd anti-correlates with density. We investigate a number of different models for the zone of avoidance in density. Supersolar metallicities and temperature effects, a strong starburst contribution in NLS1 galaxies, different NLR extents and selective obscuration are considered unlikely. Possible differences in the fraction of matter-bounded clouds and differences in the interstellar media of the host galaxies of NLS1 galaxies and broad-line Seyfert 1 (BLS1) galaxies can only be tested further with future observations. We tentatively favor the effects of winds/outflows, stronger in NLS1 galaxies than in BLS1 galaxies, to explain the observations.

Observations were made of the optical afterglow of GRB 051028 with the Lulin observatory's 1.0 m telescope and the WIDGET robotic telescope system. R band photometric data points were obtained on 2005 October 28 (UT), or 0.095-0.180 days after the burst. There is a possible plateau in the optical light curve around 0.1 days after the burst; the light curve resembles optically bright afterglows (e.g. GRB 041006, GRB 050319, GRB060605) in shape of the light curve but not in brightness. The brightness of the GRB 051028 afterglow is 3 magnitudes fainter than that of one of the dark events, GRB 020124. Optically dark GRBs have been attributed to dust extinction within the host galaxy or high redshift. However, the spectrum analysis of the X-rays implies that there is no significant absorption by the host galaxy. Furthermore, according to the model theoretical calculation of the Ly$\alpha$ absorption to find the limit of GRB 051028's redshift, the expected $R$ band absorption is not high enough to explain the darkness of the afterglow. While the present results disfavor either the high-redshift hypothesis or the high extinction scenario for optically dark bursts, they are consistent with the possibility that the brightness of the optical afterglow, intrinsically dark.

Surface lithium abundance and rotation velocity can serve as powerful and mutually complementary diagnostics of interior structure of stars. So far, the processes responsible for the lithium depletion during pre-main sequence evolution are still poorly understood. We investigate whether a correlation exists between equivalent widths of Li (EW(Li)) and rotation period (P$_{rot}$) for Weak-line T Tauri stars (WTTSs). We find that rapidly rotating stars have lower EW(Li) and the fast burning of Li begins at the phase when star's P$_{rot}$ evolves towards 3 days among 0.9M$_\odot$ to 1.4M$_\odot$ WTTSs in Taurus-Auriga. Our results support the conclusion by Piau & Turch-Chi\'eze about a model for lithium depletion with age of the star and by Bouvier et al. in relation to rotation evolution. The turn over of the curve for the correlation between EW(Li) and P$_{rot}$ is at the phase of Zero-Age Main Sequence (ZAMS). The EW(Li) decreases with decreasing P$_{rot}$ before the star reaches the ZAMS, while it decreases with increasing P$_{rot}$ (decreasing rotation velocity) for young low-mass main sequence stars. This result could be explained as an age effect of Li depletion and the rapid rotation does not inhibit Li destruction among low mass PMS stars.

We analyze the physical properties of post-starburst AGN SDSS J085338.27+033246.1 according to its optical spectrum and discuss its implications on AGN's evolution. The spectra PCA method is developed to extract emission lines and absorption features from the total light spectrum. The emission-line analysis indicates that the object can be classified as a NLS1 with FeII/H$\beta_{\rm{B}}=2.4\pm0.2$, large Eddington ratio ($\sim0.34$), small black hole mass ($\sim1.1\times10^{7}M_{\odot}$) and intermediately strong radio emission. A simple SSP model indicates that the absorption features are reproduced by a $\sim$100 Myr old starburst with a mass of $\sim7\times10^{9}M_{\odot}$ rather well. The current SFR $\sim3.0M_{\odot}\rm{yr^{-1}}$ inferred from the [OII] emission is much smaller than the past average SFR$\sim70M_{\odot}\rm{yr^{-1}}$, however. The line ratio diagnosis using the BPT diagrams indicates that the narrow emission lines are almost entirely emitted from HII regions. We further discuss a possible evolutionary path that links AGN and starburst phenomena.

We monitored the light variations of 22 weak-line T Tauri stars (WTTS) discovered among the X-ray sources in the field of the Taurus-Auriga cloud. The photometric variability of 12 out of these WTTS samples is confirmed. They are all proved to be in periodic variations. By using the methods of the Phase Dispersion Minimization (PDM) and Fourier analysis, the rotational periods of these stars were determined based on this observation. Most of them are found to be shorter than one day. This gives further evidences for the spin up of solar-type stars as predicted by models of angular momentum evolution of pre-main sequence stars.

The well-documented E1 relationships are first extended to infrared color $\alpha(60,25)$ and flux ratio [OIII]/H$\beta_{\rm{n}}$ by comparing emission line properties to continuum properties in infrared wavelength. Both direct correlations and a principle component analysis are used in a sample of 50 IRAS IR-selected Seyfert 1.5 galaxies. In addition to confirm the correlations of E1 in Boroson & Green (1992), our Eigenvector 1 turns out to be dominated by mid-infrared color $\alpha$(60,25), and most strongly effected by RFe, [OIII]/H$\beta_{\rm{n}}$, and EW(H$\beta_{\rm{b}}$). Our analysis indicate that the objects with large E1 tend to co-existent with relatively young nuclear stellar populations, which implies that the E1 is related with nuclear starformation history. The IR-dominated Eigenvector 1 can be, therefore, inferred to be interpreted as ``age'' of AGN. In confirmation of Xu et al. (2003), it is clear that the extreme Seyfert galaxies with both large RFe and large [OIII]/H$\beta_{\rm{n}}$ are rare in our Universe.

The variability of optical FeII blends in NGC 4051 is examined from spectra extracted from the AGN Watch program.In our analysis, the FeII complex are subtracted and measured with the following results. Firstly, the FeII variations were detected in NGC 4051 during a 3-year period. The identified FeII variations followed the variations in continuum closely. Secondly, the EW of FeII is reported to increase with the rising continuum flux, which is consistent with previous claims that there is no convincing Baldwin Effect in optical FeII. Thirdly, by comparing the variations of H$\beta$ and FeII, we find that RFe scales with continuum flux as $\rm{R_{Fe}\propto (5.0\pm0.8) \log(L/M)}$, which is significantly different from the theoretical expectations. Finally, in six selected Seyferts, four out of five Narrow-Line Seyfert 1 galaxies present positive correlation between RFe and the continuum flux. The negative correlations are identified in the remaindng two objects that have relatively broad profiles of H$\beta$ ($\rm{FWHM>1500 km\ s^{-1}}$).We argue that the difference of electron density of broad line clouds and/or variability behavior of incident high-energy radiation can explain the dichotomy in variability behavior of RFe.

Although the existence of asymmetrical profile of [OIII]$\lambda$5007 has been discovered for ages, its filiation and physics are poorly understood. Two new spectra of I ZWI taken on Nov 16, 2001 and on Dec 3, 2002 were compared with the spectra taken by BG92. Following results are obtained. 1)The certain variations of broad [OIII] during about 10 years separating the observations are identified. The inferred length scale of broad [OIII] emitting region ranges from 0.3pc to 3pc. By assuming a Keplerian motion in emitting region, the material emitting broad [OIII] is likely to be located at transient emission line region, between BLR and NLR. 2)We find a positive relation between the FeII emission and flux of H$\beta$(or continuum). On the other hand, the parameter RFe decreases with ionizing continuum marginally. 3)We detect a low ionized NLR in I ZWI, because of the low flux ratios $\rm{[OIII]_{n}/H\beta_{n}}$($\sim1.7$).

This work studies the optical emission line properties of a sample of 155 low-redshift bright X-ray selected ROSAT Seyfert~1 type AGN for which adequate signal-to-noise ratio spectroscopic observations are available. We measured emission line properties by performing multi-component fits to the emission line profiles, covering the effect of blended iron emission. We also obtained continuum parameters, including 250eV X-ray luminosities derived from the ROSAT database. In addition, the measured properties are gathered for a correlation analysis, which confirms the well-known relations between the strengths of Fe~II, [O III] emission and the X-ray slope. We also detect striking correlations between H$\beta$ redshift (or blueshift), flux ratios of Fe `II to H$\beta$ broad component and [O III] to H$\beta$ narrow component. These trends are most likely driven by the Eddington ratio.

1RXS J232953.9+062814 was identified as a cataclysmic variable by Wei et al. (1999). Four low-resolution spectra of 1RXS J232953.9+062814 were obtained by using the 2.16-m telescope of the National Astronomical Observatories, in which two of them were at outburst, and the other two were at quiescence. The system is about 16.8 B and 16.5 V at quiescence, and 12.6 B and 12.6 V at outburst. The quiescent spectra were dominated by double-peaked Balmer emissions, which indicates a hydrogen-rich system with a high-inclination accretion disc. MgH and TiO absorption bands appeared in the quiescent spectrum imply a companion with a spectral type of early M dwarf. If we take it as a M0 dwarf, the system is located at a distance of 350 pc with a proper motion velocity 150 km s$^{-1}$. The superhump period of 0.046311 days (Uemura et al. 2001) was confirmed by our V photometry. The short period and the hydrogen-rich nature reveal that this system is another SU Ursae Majoris-type dwarf nova below the period minimum after V485 Centauri. 1RXS J232953.9+062814 is one of the most important systems for studying the evolutionary scenario of cataclysmic variables since it is much brighter than V485 Cen.

In paper 1 (Cheng et al. 2001) we have shown that the soft X-ray spectra of two types of Seyfert 1 galaxies statistically vary differently with increasing intensity. In order to understand how the spectrum of blazars changes, the spectral shape variability of 18 blazars observed by ROSAT/PSPC mode are studied by presenting the correlation of Hardness Ratio 1 versus Count Rates (HR1-CTs). According to our criteria, 10 blazars show a positive HR1-CTs relation, and only 2 blazars display an anti-correlation of HR1 versus CTs. The rest 6 blazars do not indicate any clear correlation. From these we can see that most blazars of our sample statistically show a hardening spectrum during overall flux increase, though some vary randomly. By investigating the photon index of these objects and different radiation theories, we argue that the dominance of the synchrotron or inverse Compton emission in the soft X-ray band may interpret the dichotomy of spectral variability well, and that different spectral variations might represent a sequence of synchrotron peaked frequency.

In order to understand how the soft X-ray spectra vary we present the Hardness Ratio 1 versus Count Rates (HR1-CTs) correlation of 8 Narrow-line Seyfert 1 Galaxies (NLS1s) and 14 Broad-line Seyfert1 Galaxies (BLS1s) obtained during the ROSAT PSPC pointing observations. According to our criteria, six of the NLS1s show a positive HR1-CTs correlation, and seven of the BLS1s display an anti-correlation of HR1 versus CTs. The other 2 NLS1s and 7 BLS1s do not show a clear HR1-CTs correlation. From these we can see that the NLS1s statistically show a different spectral shape variability with flux change from the BLS1s: the spectra of NLS1s become harder as total flux increases while those of BLS1s soften. We attribute the different spectral variations to a strong stable 'soft excess' in NLS1s, while it is weak in BLS1s. For two types of objects, the power law component similarly becomes softer with increasing intensity. These imply that the soft excess originates from the Big Blue Bump and power law emission is from Compton upscattering of UV or Soft X-ray photons. Our results are consistent with what is widely accepted that NLS1s have smaller black hole masses and higher accretion rates than BLS1s.

New optical and radio data on the BL Lacertae object AO 0235+16 have been collected in the last four years by a wide international collaboration, which confirm the intense activity of this source. The optical data also include the results of the Whole Earth Blazar Telescope (WEBT) first-light campaign organized in November 1997. The optical spectrum is observed to basically steepen when the source gets fainter. We have investigated the existence of typical variability time scales and of possible correlations between the optical and radio emissions by means of visual inspection, Discrete Correlation Function analysis, and Discrete Fourier Transform technique. The major radio outbursts are found to repeat quasi-regularly with a periodicity of about 5.7 years; this period is also in agreement with the occurrence of some of the major optical outbursts, but not all of them.

We present the first ROSAC results of an AGN clustering analysis. This study comprises a sample of 200 AGNs, 75% of which being at low redshifts z<0.5, in the Ursa Major constellation. The spatial 2-point-correlation function (SCF) as well as the minimal spanning tree (MST) technique were applied. Some evidence for clustering is found in the SCF, although with low significance. Using the MST technique, we could find two AGN groups. This result is preliminary and the exact significance will be tested with careful simulations.

We are searching for large-scale structures in the distribution of AGN discovered by the ROSAT All-Sky Survey. The RASS detected > 60000 X-ray objects, of which approximately 35% are AGN at z < 0.5. The surface density in the extragalactic sky is approximately 0.5 AGN/deg^2, which has not been reached until now by any other survey for this redshift range. We efficiently single out the AGN among all the RASS sources using the Hamburg/RASS database of optical identifications, which contains presently 13867 entries. Follow-up spectroscopy of RASS AGN candidates identified in selected areas of the northern sky is underway to determine the spatial distribution of the AGN in these areas. In every area structures reminiscent of clusters and filaments are found on scales 50 - 100 h^{-1}Mpc. These structures have a similar size as the serendipitously discovered groups of AGN in optical surveys at higher redshifts. Samples of low redshift AGN drawn from the RASS are large enough to be used with samples of higher redshift to study evolutionary effects in the large-scale distribution of AGN.

We report the discovery of two new Wolf-Rayet (WR) galaxies: Mrk~1039, and F08208$+$2816. Two broad WR bumps at 5808\AA~ and 4650\AA~ indicate the presence of WCE and WNL star population in all two sources. We also confirm the presenceof WR features in Mrk~35, previously detected in a different position. The observed equivalent width of the WR bump at 4650\AA~ and the derived number ratios of WR/(WR$+$O) imply that star formation in these sources takes place inshort burst duration. Comparisons with the recent models of WR populations in young starbursts with the observed EW(\HeII)/EW(\CIV)/EW(WRbump) and their relative intensitie provide an indication that the stellar initial mass function in some WR galaxies might not be Salpeter-like. It is interesting to find that the luminous IRAS source, F08208$+$2816, has little dust reddening, probably because of the existence of a powerful superwind. By comparisons with other starbursts observed with the Hopkins Ultraviolet Telescope, F08208$+$2816 as a merging system renders a chance to study the contribution from young starbursts to the UV background radiation in universe.

We present the results of three years of ground-based observations of the Seyfert 1 galaxy NGC 5548, which combined with previously reported data, yield optical continuum and broad-line H-beta light curves for a total of eight years. The light curves consist of over 800 points, with a typical spacing of a few days between observations. During this eight-year period, the nuclear continuum has varied by more than a factor of seven, and the H-beta emission line has varied by a factor of nearly six. The H-beta emission line responds to continuum variations with a time delay or lag of 10-20 days, the precise value varying somewhat from year to year. We find some indications that the lag varies with continuum flux in the sense that the lag is larger when the source is brighter.

We report redshift measurements of 5 X-ray emitting blue stellar objects (BSO's) located less than 12 arc min from the X-ray Seyfert galaxy, NGC 3516. We find these quasars to be distributed along the minor axis of the galaxy and to show a very good correlation between their redshift and their angular distance from NGC 3516. Moreover the redshifts of these 5 quasars are: 0.33, 0.69, 0.93, 1.40 and 2.10 which are very near the peaks of the redshift periodicity distribution (i.e. z=0.3, 0.6, 0.96, 1.41 and 1.96). All these observed properties strikingly confirm, around this single example of a Seyfert, the composite picture derived from previous physical associations of quasars with low redshift, active galaxies.