We use observational measurements of thermal and kinetic Sunyaev-Zel'dovich effects, as well as soft X-ray emission associated with galaxy groups to constrain the gas density and temperature in the circumgalactic medium (CGM) for dark matter halos with masses above $10^{12.5}M_{\odot}$. A number of generic models are used together with a Bayesian scheme to make model inferences. We find that gas with a single temperature component cannot provide a consistent model to match the observational data. A simple two-phase model assuming a hot component and an ionized warm component can accommodate all the three observations. The total amount of the gas in individual halos is inferred to be comparable to the universal baryon fraction corresponding to the halo mass. The inferred temperature of the hot component is comparable to the halo virial temperature. The fraction of the hot component increases from $(15-40)\%$ for $10^{12.5}{M}_\odot$ halos to $(40-60)\%$ for $10^{14.5}{M}_\odot$ halos, where the ranges reflect uncertainties in the assumed gas density profile. Our results suggest that a significant fraction of the halo gas is in a non-thermalized component with temperature much lower than the virial temperature.

SiTian is an ambitious ground-based all-sky optical monitoring project, developed by the Chinese Academy of Sciences. The concept is an integrated network of dozens of 1-m-class telescopes deployed partly in China and partly at various other sites around the world. The main science goals are the detection, identification and monitoring of optical transients (such as gravitational wave events, fast radio bursts, supernovae) on the largely unknown timescales of less than 1 day; SiTian will also provide a treasure trove of data for studies of AGN, quasars, variable stars, planets, asteroids, and microlensing events. To achieve those goals, SiTian will scan at least 10,000 square deg of sky every 30 min, down to a detection limit of $V \approx 21$ mag. The scans will produce simultaneous light-curves in 3 optical bands. In addition, SiTian will include at least three 4-m telescopes specifically allocated for follow-up spectroscopy of the most interesting targets. We plan to complete the installation of 72 telescopes by 2030 and start full scientific operations in 2032.

We report the discovery of a luminous quasar, J1007+2115 at $z=7.515$ ("P\={o}niu\={a}'ena"), from our wide-field reionization-era quasar survey. J1007+2115 is the second quasar now known at $z>7.5$, deep into the reionization epoch. The quasar is powered by a $(1.5\pm0.2)\times10^9$ $M_{\odot}$ supermassive black hole (SMBH), based on its broad MgII emission-line profile from Gemini and Keck near-IR spectroscopy. The SMBH in J1007+2115 is twice as massive as that in quasar J1342+0928 at $z=7.54$, the current quasar redshift record holder. The existence of such a massive SMBH just 700 million years after the Big Bang significantly challenges models of the earliest SMBH growth. Model assumptions of Eddington-limited accretion and a radiative efficiency of 0.1 require a seed black hole of $\gtrsim 10^{4}$ $M_{\odot}$ at $z=30$. This requirement suggests either a massive black hole seed as a result of direct collapse or earlier periods of rapid black hole growth with hyper-Eddington accretion and/or a low radiative efficiency. We measure the damping wing signature imprinted by neutral hydrogen absorption in the intergalactic medium (IGM) on J1007+2115's Ly$\alpha$ line profile, and find that it is weaker than that of J1342+0928 and two other $z\gtrsim7$ quasars. We estimate an IGM volume-averaged neutral fraction $\langle x\rm_{HI}\rangle=0.39^{+0.22}_{-0.13}$. This range of values suggests a patchy reionization history toward different IGM sightlines. We detect the 158 $\mu$m [C II] emission line in J1007+2115 with ALMA; this line centroid yields a systemic redshift of $z=7.5149\pm0.0004$ and indicates a star formation rate of $\sim210$ $M_{\odot}$ yr$^{-1}$ in its host galaxy.

By assuming the formation of a black hole soon after the merger event of GW170817, Shibata et al. updated the constraints on the maximum mass ($M_\textrm{max}$) of a stable neutron star within $\lesssim$ 2.3 $M_{\odot}$, but there is no solid evidence to rule out $M_\textrm{max}>2.3~M_{\odot}$ from the point of both microphysical and astrophysical views. In order to explain massive pulsars, it is naturally expected that the equation of state (EOS) would become stiffer beyond a specific density. In this paper, we consider the possibility of EOSs with $M_\textrm{max}>2.3~M_{\odot}$, investigating the stiffness and the transition density in a polytropic model. Two kinds of neutron stars are considered, i.e., normal neutron stars (the density vanishes on gravity-bound surface) and strange stars (a sharp density discontinuity on self-bound surface). The polytropic model has only two parameter inputs in both cases: ($\rho_{\rm t}$, $\gamma$) for gravity-bound objects, while ($\rho_{\rm s}$, $\gamma$) for self-bound ones, with $\rho_{\rm t}$ the transition density, $\rho_{\rm s}$ the surface density and $\gamma$ the polytropic exponent. In the matter of $M_\textrm{max}>2.3~M_{\odot}$, it is found that the smallest $\rho_{\rm t}$ and $\gamma$ should be $\sim 0.50~\rho_0$ and $\sim 2.65$ for normal neutron stars, respectively, whereas for strange star, we have $\gamma > 1.40$ if $\rho_{\rm s} > 1.0~\rho_0$ and $\rho_{\rm s} < 1.58~\rho_0$ if $\gamma <2.0$ ($\rho_0$ is the nuclear saturation density). These parametric results could guide further research of the real EOS with any foundation of microphysics if a pulsar mass higher than $2.3~M_{\odot}$ is measured in the future. We also derive rough results of common neutron star radius range, which is $9.8~\rm{km} < R_{1.4} < 13.8~\rm{km}$ for normal neutron stars and $10.5~\rm{km} < R_{1.4} < 12.5~\rm{km}$ for strange stars.

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

A Kepler solver is an analytical method used to solve a two-body problem. In this paper, we propose a new correction method by slightly modifying the Kepler solver. The only change to the analytical solutions is that the obtainment of the eccentric anomaly relies on the true anomaly that is associated to a unit radial vector calculated by an integrator. This scheme rigorously conserves all integrals and orbital elements except the mean longitude. However, the Kepler energy, angular momentum vector and Laplace-Runge-Lenz vector for perturbed Kepler problems are slowly-varying quantities. However, their integral invariant relations give the quantities high-precision values that directly govern five slowly-varying orbital elements. These elements combined with the eccentric anomaly determine the desired numerical solutions. The newly proposed method can considerably reduce various errors for a post-Newtonian two-body problem compared with an uncorrected integrator, making it suitable for a dissipative two-body problem. Spurious secular changes of some elements or quasi-integrals in the outer solar system may be caused by short integration times of the fourth-order Runge-Kutta algorithm. However, they can be eliminated in a long integration time of 10^8 years by the proposed method, similar to Wisdom-Holman second-order symplectic integrator. The proposed method has an advantage over the symplectic algorithm in the accuracy but gives a larger slope to the phase error growth.

When a magnetar magnetosphere is trigged by crustal deformations, an electric field $E_\parallel$ parallel to the magnetic field line would appear via Alvf\'en waves in the charge starvation region. If electron-positron pair bunches pre-exist, e.g., via some possible plasma instabilities, in the magnetosphere, these pairs will undergo charge separation in the $E_\parallel$ and in the meantime emit coherent curvature radiation. Following the approach of Yang & Zhang (2018), we find that the superposed curvature radiation becomes narrower due to charge separation, with the width of spectrum depending on the separation between the electron and positron clumps. This mechanism can interpret the narrow spectra of FRBs, in particular, the Galactic FRB 200428 recently detected in association with a hard X-ray burst from the Galactic magnetar SGR J1935+2154.

Pulsar's position, proper motion and parallax are important parameters in timing equations. It is a really challenging work to fit astrometric parameters accurately through pulsar timing, especially for pulsars that show irregular timing properties. As the fast development of related techniques, it is possible to measure astrometric parameters of more and more pulsars in a model$\textrm{-}$independent manner with the Very Long Baseline Interferometry (VLBI). In this work, we select 16 normal pulsars, whose parallax and proper motion have not been successfully fitted with timing observations or show obvious differences with corresponding latest VLBI solutions, and do further studies on their timing properties. After updating astrometric parameters in pulsar ephemerides with the latest VLBI measurements, we derive the latest rotation solutions of these pulsars with observation data at S and C$\textrm{-}$band obtained from the Shanghai Tian Ma Radio Telescope (TMRT). Compared with spin frequency $\nu$ inferred from previous rotation solutions, the newly$\textrm{-}$fitted $\nu$ show differences larger than 10$^{-9}$ Hz for most pulsars. The contribution of the Shklovsky effect to period derivative $\dot{P}$ can be properly removed taking advantages of accurate proper motion and distance of target pulsars measured by VLBI astrometry. This further leads to a precise estimate of intrinsic characteristic age $\tau_{\rm c}$. Differences between the newly$\textrm{-}$measured $\tau_{\rm c}$ and corresponding previous results are as large as 2% for some pulsars. VLBI astrometric parameter solutions also lead to better measurements of timing irregularities. For PSR B0154$+$61, the glitch epoch (MJD 58279.5) measured with previous ephemeris is about 13 d later than the result (MJD 58266.4) obtained with VLBI astrometric parameter solutions.

Fast radio bursts (FRBs) are mysterious millisecond-duration radio transients observed from extragalactic distances. Magnetars have been long speculated as the possible engine to power repeating bursts from FRB sources, but no convincing evidence has been collected so far. Recently, a Galactic magnetar dubbed Soft Gamma-ray Repeater (SGR) J1935+2154 entered an active phase by emitting intense soft gamma-ray bursts. One fast radio burst with two peaks (hereafter FRB 200428) and a luminosity comparable to those of extragalactic FRBs was detected from the source in association with a soft $\gamma$-ray / hard X-ray flare. Here we report an eight-hour targeted observational campaign with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) in four observing sessions, assisted by multi-wavelength (optical and hard X-rays) observations. During the third session, twenty-nine SGR bursts were detected by the {\it Fermi} Gamma-ray Burst Monitor (GBM). Throughout the observing period, especially during the scrutinized epochs when SGR bursts reached the Earth, no single dispersed pulsed emission was detected by FAST. This places a fluence upper limit eight orders of magnitude deeper than the fluence of FRB 200428. Our results suggest that FRB -- SGR burst associations are rather rare. FRBs may be highly relativistic and geometrically beamed or may carry narrow spectra with the characteristic frequencies of the most bursts far outside the FAST band. More intriguingly, it is also possible that the physical condition to achieve coherent radiation in SGR bursts is contrived. Only under extreme conditions could an FRB be made in association with an SGR burst.

Using a high luminosity ($L_{\rm bol} \sim 10^{47.5} - 10^{48.3}$ erg s$^{-1}$), high redshift ($3.2 < z < 3.8$) quasar sample of 19 quasars with optical and near-infrared spectroscopy, we investigate the reliability of the CIV-based black hole mass estimates ($M_{\rm BH}$). The median logarithm of the CIV- and H$\beta$-based $M_{\rm BH}$ ratios is 0.110 dex with the scatter of 0.647 dex. The CIV-to-H$\beta$ BH mass differences are significantly correlated with the CIV FWHMs, blueshifts and asymmetries. Corrections of the CIV FWHM using the blueshift and asymmetry reduce the scatter of the mass differences by $\sim$ 0.04-0.2 dex. Quasars in our sample accrete at the Eddington ratio $R_{\rm EDD}>0.3$ and cover a considerable range of blueshifts, with 18/19 of the quasars showing CIV blueshifts (with the median value of 1126 km s$^{-1}$) and 14/19 of the quasars showing CIV blueshifts larger than 500 km s$^{-1}$. It suggests that not all quasars with high Eddington ratios show large blueshifts. The Baldwin effect between the CIV rest-frame equivalent width (REW) and the continuum luminosity at 1350 $\overset{\lower.5em\circ}{\mathrm{A}}$ is not seen, likely due to the limited luminosity range of our sample. We find a lack of flux in the red wing of the composite spectrum with larger CIV blueshift, and detect a higher ratio of [OIII] quasars with REW$_{\rm [OIII]}>5$ $\overset{\lower.5em\circ}{\mathrm{A}}$ in the subsample with lower CIV blueshift. It is more likely that they are caused by the combination of the Eddington ratio and the orientation effect.

We propose a geometrical explanation for the periodically and nonperiodicly repeating fast radio bursts (FRBs) under neutron star (NS)-companion system. We suggest a critical binary separation, $r_{\rm c}$, within which the interaction between the NS and its companion can trigger FRBs. On one hand, for an elliptic orbit with the minimum and maximum binary separations, $r_{\rm min}$ and $r_{\rm max}$ (satisfying the condition $r_{\rm min}<r_{\rm c}<r_{\rm max}$), this binary system could reproduce periodical FRBs. A small separation $r$ ($<r_{\rm c}$) of such an elliptic orbit would result in an active period of periodical FRBs. On the other hand, if $r_{\rm max}<r_{\rm c}$, the modulation of orbital motion will not work due to persistent interaction, and this kind of repeating FRBs should be nonperiodic. It is estimated that the critical separation $r_{\rm c}\sim 10^{12}$ cm for stellar binary system.

Luminous $z\ge7$ quasars provide direct probes of the evolution of supermassive black holes (SMBHs) and the intergalactic medium (IGM) during the epoch of reionization (EoR). The Ly$\alpha$ damping wing absorption imprinted by neutral hydrogen in the IGM can be detected in a single EoR quasar spectrum, allowing the measurement of the IGM neutral fraction towards that line of sight. However, damping wing features have only been detected in two $z>7$ quasars in previous studies. In this paper, we present new high quality optical and near-infrared spectroscopy of the $z=7.00$ quasar DES J025216.64--050331.8 obtained with Keck/NIRES and Gemini/GMOS. By using the MgII single-epoch virial method, we find that it hosts a $\rm (1.39\pm0.16) \times10^{9} ~M_\odot$ SMBH accreting at an Eddington ratio of $\lambda_{\rm Edd}=0.7\pm0.1$, consistent with the values seen in other luminous $z\sim 7$ quasars. Furthermore, the Ly$\alpha$ region of the spectrum exhibits a strong damping wing absorption feature. The lack of associated metal absorption in the quasar spectrum indicates that this absorption is imprinted by a neutral IGM. Using a state-of-the-art model developed by Davies et al., we measure a volume-averaged neutral hydrogen fraction at $z=7$ of $\langle x_{\rm HI} \rangle = 0.70^{+0.20}_{-0.23} (^{+0.28}_{-0.48})$ within 68% (95%) confidence intervals when marginalizing over quasar lifetimes of $10^3\le t_{\rm Q}\le10^8$ yr. This is the highest IGM neutral fraction yet measured using reionization-era quasar spectra.

The reverse shock (RS) model is generally introduced to interpret the optical afterglows with the rapid rising and decaying, such as the early optical afterglow of GRB 990123 (which is also called optical flash). In this paper, we collected 11 gamma-ray bursts (GRBs) early optical afterglows, which have such signatures of dominant RS emission. Since the temporal slopes of the optical flashes are determined by both the medium density distribution index $k$ and the electron spectral index $p$, we apply the RS model of the thin shell case to the optical flashes and determine the ambient medium of the progenitors. We find that the $k$ value is in the range of 0 - 1.5. The $k$ value in this paper is consistent with the result in Yi et al. (2013), where the forward shock (FS) model was applied to some onset bumps. However, the method adopted in this paper is only applicable to GRB afterglows with significant sharp rising and decaying RS emission. Our results indicate that the RS model can also be applied to confirm the circumburst medium, further implying that GRBs may have diverse circumburst media.

A large number of observations from the Parkes 64\,m-diameter radio telescope, recorded with high time resolution, are publicly available. We have re-processed all of the observations obtained during the first four years (from 1997 to 2001) of the Parkes Multibeam receiver system in order to identify transient events and have built a database that records the 568,736,756 pulse candidates generated during this search. We have discovered a new fast radio burst (FRB), FRB~010305, with a dispersion measure of 350$\pm$5\,\,cm$^{-3}\,$pc and explored why so few FRBs have been discovered in data prior to 2001. After accounting for the dispersion smearing across the channel bandwidth and the sky regions surveyed, the number of FRBs is found to be consistent with model predictions. We also present five single pulse candidates from unknown sources, but with Galactic dispersion measures. We extract a diverse range of sources from the database, which can be used, for example, as a training set of data for new software being developed to search for FRBs in the presence of radio frequency interference.

We study the photometric properties and sizes of the reionization-epoch galaxies in high-resolution Simba cosmological hydrodynamical simulations with box sizes of $[25,50]~h^{-1}$Mpc. Assuming various attenuation laws, we compute photometry by extincting each star particle's spectrum using the line-of-sight gas metal column density. The predicted ultraviolet luminosity function (UVLF) generally agrees with observations at $z=6$, owing to a partial cancellation between the high metallicities of the simulated galaxies and lower dust-to-metal ratios. The simulated $z=8$ UVLF is low compared to observations, likely owing to excessive dust extinction. Simba predicts UV continuum slopes ($\beta$) in agreement with the $z=6$ observations, with the best agreement obtained using a Calzetti extinction law. Interestingly, the gas-phase mass-metallicity relation in Simba is higher at $z\sim 6$ than at $z\sim 2$, suggesting that rapid early enrichment (and dust growth) might be necessary to match the observed $\beta$. We find that $\beta$ is more sensitive to the dust extinction law than the UVLF. By generating mock James Webb Space Telescope (JWST) images and analysing in a manner similar to observations, we show that Simba's galaxy size-luminosity relation well reproduces the current $z=6$ Hubble observations. Unlike observations at lower redshifts, Simba predicts similar rest-UV and rest-optical sizes of $z=6$ galaxies, owing to weak age gradients and dust extinction in star-forming regions counteract each other to weaken the color gradients within galaxies. These predictions will be testable with JWST.

We present an observational study of 22 broad absorption line quasars (BAL QSOs) at $3\lesssim z \lesssim5$ based on optical/near-IR spectroscopy, aiming to investigate quasar winds and their effects. The near-IR spectroscopy covers the \hb\ and/or \mgii\ broad emission lines (BELs) for these quasars, allowing us to estimate their central black hole (BH) masses in a robust way. We found that our BAL QSOs on average do not have a higher Eddington ratio than that from non-BAL QSOs matched in redshift and/or luminosity. In a subset consisting of seven strong BAL QSOs possessing sub-relativistic BAL outflows, we see the prevalence of large \civ-BEL blueshift ($\sim$3100 km s$^{-1}$) and weak \oiii\ emission (particularly the narrow \oiii$\lambda$5007 component), indicative of nuclear outflows affecting the narrow emission-line (NEL) regions. In another subset consisting of thirteen BAL QSOs having simultaneous observations of \mgii\ and \hb, we found a strong correlation between 3000~\AA\ and 5000~\AA\ monochromatic luminosity, consistent with that from non-BAL QSOs matched in redshift and luminosity; however, there is no correlation between \mgii\ and \hb\ in FWHM, likely due to nuclear outflows influencing the BEL regions. Our spectroscopic investigations offer strong evidence that the presence of nuclear outflows plays an important role in shaping the BEL/NEL regions of these quasars and possibly, regulating the growth of central supermassive black holes (SMBHs). We propose that BEL blueshift and BAL could be different manifestations of the same outflow system viewed at different sightlines and/or phases.

We present the observational results from the 43-GHz Very Long Baseline Array (VLBA) observations of 124 compact radio-loud active galactic nuclei (AGNs) that were conducted between 2014 November and 2016 May. The typical dimensions of the restoring beam in each image are about 0.5 mas $\times$ 0.2 mas. The highest resolution of 0.2 mas corresponds to a physical size of 0.02 pc for the lowest redshift source in the sample. The 43-GHz very long baseline interferometry (VLBI) images of 97 AGNs are presented for the first time. We study the source compactness on milli-arcsec (mas) and sub-mas scales, and suggest that 95 sources in our sample are suitable for future space VLBI observations. By analyzing our data supplemented with other VLBA AGN surveys from literature, we find that the core brightness temperature increases with increasing frequency below a break frequency ~ 7 GHz, and decreases between ~7--240~GHz but increases again above~240 GHz in the rest frame of the sources. This indicates that the synchrotron opacity changes from optically thick to thin. We also find a strong statistical correlation between radio and $\gamma$-ray flux densities. Our correlation is tighter than those in literature derived from lower-frequency VLBI data, suggesting that the $\gamma$-ray emission is produced more co-spatially with the 43-GHz VLBA core emission. This correlation can also be extrapolated to the un-beamed AGN population, implying that a universal $\gamma$-ray production mechanism might be at work for all types of AGNs.

Ground based optical telescopes are seriously affected by atmospheric turbulence induced aberrations. Understanding properties of these aberrations is important both for instruments design and image restoration methods development. Because the point spread function can reflect performance of the whole optic system, it is appropriate to use the point spread function to describe atmospheric turbulence induced aberrations. Assuming point spread functions induced by the atmospheric turbulence with the same profile belong to the same manifold space, we propose a non-parametric point spread function -- PSF-NET. The PSF-NET has a cycle convolutional neural network structure and is a statistical representation of the manifold space of PSFs induced by the atmospheric turbulence with the same profile. Testing the PSF-NET with simulated and real observation data, we find that a well trained PSF--NET can restore any short exposure images blurred by atmospheric turbulence with the same profile. Besides, we further use the impulse response of the PSF-NET, which can be viewed as the statistical mean PSF, to analyze interpretation properties of the PSF-NET. We find that variations of statistical mean PSFs are caused by variations of the atmospheric turbulence profile: as the difference of the atmospheric turbulence profile increases, the difference between statistical mean PSFs also increases. The PSF-NET proposed in this paper provides a new way to analyze atmospheric turbulence induced aberrations, which would be benefit to develop new observation methods for ground based optical telescopes.

The inclination distribution of circumbinary planets (CBPs) is an important scientific issue. It is of great significance in estimating the occurrence rate of CBPs and studying their formation and evolution. Although the CBPs currently discovered by the transit method are nearly coplanar, the true inclination distribution of CBPs is still unknown. Previous researches on CBPs mostly regarded them as an isolated binary-planet system, without considering the birth environment of their host binaries. It is generally believed that almost all stars are born in clusters. Therefore, it is necessary to consider the impact of the close encounters of stars on the CBP systems. This article discusses how the close encounters of fly-by stars affect the inclination of CBP. Based on extensive numerical simulations, we found that CBPs in close binary with a spacing of about 0.2 au are almost unaffected by star fly-bys. Their orbits remain coplanar. However, when the spacing of the binary stars is greater than 1 au, the 2-3 fly-bys of the intruding star can excite a considerable inclination even for the CBP near the unstable boundary of the binary. For the planets in the outer region, a single star fly-by can excite an inclination to more than 5 degrees. Especially, CBPs in near polar or retrograde orbits can naturally form through bianry-star encounters. If close binaries are born in open clusters, our simulations suggest that there may be high-inclination CBPs in binaries with spacing > 1 au.

We analyze the X-ray, optical, and mid-infrared data of a "changing-look" Seyfert galaxy \sdssj15 at $z\simeq0.086$. Over a period of one decade (2009 - 2018), its broad H$\alpha$ line intensity increased by a factor of $\sim$4. Meanwhile, the X-ray emission in 2014 as observed by \chandra\ was about five times brighter than that in 2010 by {\it Suzaku}, and the corresponding emissions in V-band, mid-infrared W1 band brighten by $\sim$ 0.18, 0.32 mag, respectively. Moreover, the absorption in X-rays is moderate and stable, i.e. ${\rm N_{H}}\sim 10^{21} {\rm cm^{-2}}$, but the X-ray spectrum becomes harder in the 2014 \chandra\ bright state (i.e. photon index $\Gamma = 1.52^{+0.06}_{-0.06}$) than that of the 2010 \suzaku\ low state ($\Gamma=2.03^{+0.22}_{-0.21}$). With an Eddington ratio being lower than a few percent, the inner region of the accretion disk in \sdssj15\ is likely a hot accretion flow. We then compile from literature the X-ray data of "changing-look" AGNs, and find that they generally follow the well-established "V"-shaped correlation in AGNs, that is, above a critical turn-over luminosity the X-ray spectra soften with the increasing luminosity, and below that luminosity the trend is reversed in a way of "harder when brighter". This presents a direct evidence that CL-AGNs have distinctive changes in not only the optical spectral type, but also the X-ray spectral shape. The similarity in the X-ray spectral evolution between CL-AGNs and black hole X-ray binaries indicates that the observed CL-AGNs phenomena may relate to the state transition in accretion physics.

The Beijing-Arizona Sky Survey (BASS) is a wide and deep imaging survey to cover a 5400 deg$^2$ area in the Northern Galactic Cap with the 2.3m Bok telescope using two filters ($g$ and $r$ bands). The Mosaic $z$-band Legacy Survey (MzLS) covers the same area in $z$ band with the 4m Mayall telescope. These two surveys will be used for spectroscopic targeting of the Dark Energy Spectroscopic Instrument (DESI). The BASS survey observations were completed in 2019 March. This paper describes the third data release (DR3) of BASS, which contains the photometric data from all BASS and MzLS observations between 2015 January and 2019 March. The median astrometric precision relative to {\it Gaia} positions is about 17 mas and the median photometric offset relative to the PanSTARRS1 photometry is within 5 mmag. The median $5\sigma$ AB magnitude depths for point sources are 24.2, 23.6, and 23.0 mag for $g$, $r$, and $z$ bands, respectively. The photometric depth within the survey area is highly homogeneous, with the difference between the 20\% and 80\% depth less than 0.3 mag. The DR3 data, including raw data, calibrated single-epoch images, single-epoch photometric catalogs, stacked images, and co-added photometric catalogs, are publicly accessible at \url{http://batc.bao.ac.cn/BASS/doku.php?id=datarelease:home}.

The characteristics of the spectral evolution of the prompt emission of gamma-ray bursts (GRBs), which are closely related to the radiation mechanism (synchrotron or photosphere), are still an unsolved subject. Here, by performing the detailed time-resolved spectral fitting of GRB 131231A, which has a very bright and well-defined single pulse, some interesting spectral evolution features have been found. (i) Both the low-energy spectral index $\alpha$ and the peak energy $E_{\rm p}$ exhibit the "flux-tracking" pattern ("double-tracking" characteristics). (ii) The parameter relations, i.e., $F$ (the energy flux)-$\alpha$, $F$-$E_{\rm p}$, and $E_{\rm p}$-$\alpha$, along with the analogous Yonetoku $E_{\rm p}$-$L_{\gamma,\rm iso}$ relation for the different time-resolved spectra, show strong monotonous (positive) correlations, both in the rising and the decaying phases. (iii) The values of $\alpha$ do not exceed the synchrotron limit ($\alpha$= -2/3) in all slices across the pulse, favoring the synchrotron origin. We argue that the one-zone synchrotron emission model with the emitter streaming away at a large distance from the central engine can explain all of these special spectral evolution characteristics.

In order to find clues to the origin of the "winged" or "X-shaped" radio galaxies (XRGs) we investigate here the parent galaxies of a large sample of 106 XRGs for optical-radio axes alignment, interstellar medium, black hole mass, and large-scale environment. For 41 of the XRGs it was possible to determine the optical major axis and the primary radio axis and the strong tendency for the two axes to be fairly close is confirmed. However, several counter-examples were also found and these could challenge the widely discussed backflow diversion model for the origin of the radio wings. Comparison with a well-defined large sample of normal FR II radio galaxies has revealed that: (i) XRGs possess slightly less massive central black holes than the normal radio galaxies (average masses being log$M_{\rm BH} \sim$ 8.81 $M_{\odot}$ and 9.07 $M_{\odot}$, respectively); (ii) a much higher fraction of XRGs ($\sim$ 80%) exhibits red mid-IR colors ($W2 - W3 > 1.5$), indicating a population of young stars and/or an enhanced dust mass, probably due to relatively recent galaxy merger(s). A comparison of the large-scale environment (i.e., within $\sim$ 1 Mpc) shows that both XRGs and FRII radio galaxies inhabit similarly poor galaxy clustering environments (medium richness being 8.94 and 11.87, respectively). Overall, the origin of XRGs seems difficult to reconcile with a single dominant physical mechanism and competing mechanisms seem prevalent.

The POLAR detector is a space based Gamma-Ray Burst (GRB) polarimeter sensitive in the 15-500 keV energy range. Apart from its main scientific goal as a Gamma-Ray Burst polarimeter it is also able to detect photons from pulsars in orbit. By using the six-months in-orbit observation data, significant pulsation from the PSR B0531+21 (Crab pulsar) was obtained. In this work, we present the precise timing analysis of the Crab pulsar, together with a phase-resolved spectroscopic study using a joint-fitting method adapted for wide field of view instruments like POLAR. By using single power law fitting over the pulsed phase, we obtained spectral indices ranging from 1.718 to 2.315, and confirmed the spectral evolution in a reverse S shape which is homogenous with results from other missions over broadband. We will also show, based on the POLAR in-orbit performance and Geant4 Monte-Carlo simulation, the inferred capabilities of POLAR-2, the proposed follow-up mission of POLAR on board the China Space Station (CSS), for pulsars studies.

Reionization leads to large spatial fluctuations in the intergalactic temperature that can persist well after its completion. We study the imprints of such fluctuations on the $z\sim5$ Ly$\alpha$ forest flux power spectrum using a set of radiation-hydrodynamic simulations that model different reionization scenarios. We find that large-scale coherent temperature fluctuations bring $\sim20-60\%$ extra power at $k\sim0.002$ s/km, with the largest enhancements in the models where reionization is extended or ends the latest. On smaller scales ($k\gtrsim0.1$ s/km), we find that temperature fluctuations suppress power by $\lesssim10\%$. We find that the shape of the power spectrum is mostly sensitive to the reionization midpoint rather than temperature fluctuations from reionization's patchiness. However, for all of our models with reionization midpoints of $z\le 8$ ($z\le 12$) the shape differences are $\lesssim20\%$ ($\lesssim40\%$) because of a surprisingly well-matched cancellation between thermal broadening and pressure smoothing that occurs for realistic thermal histories. We also consider fluctuations in the ultraviolet background, finding their impact on the power spectrum to be much smaller than temperature fluctuations at $k\gtrsim0.01$ s/km. Furthermore, we compare our models to power spectrum measurements, finding that none of our models with reionization midpoints of $z<8$ is strongly preferred over another and that all of our models with midpoints of $z\geq8$ are excluded at $2.5\sigma$. Future measurements may be able to distinguish between viable reionization models if they can be performed at lower $k$ or, alternatively, if the error bars on the high-$k$ power can be reduced by a factor of $1.5$.

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.

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

Jupiter has the most energetic and complex radiation belts in our solar system. Their hazardous environment is the reason why so many spacecraft avoid rather than investigate them, and explains how they have kept many of their secrets so well hidden, despite having been studied for decades. In this White Paper we argue why these secrets are worth unveiling. Jupiter's radiation belts and the vast magnetosphere that encloses them constitute an unprecedented physical laboratory, suitable for both interdisciplinary and novel scientific investigations: from studying fundamental high energy plasma physics processes which operate throughout the universe, such as adiabatic charged particle acceleration and nonlinear wave-particle interactions; to exploiting the astrobiological consequences of energetic particle radiation. The in-situ exploration of the uninviting environment of Jupiter's radiation belts present us with many challenges in mission design, science planning, instrumentation and technology development. We address these challenges by reviewing the different options that exist for direct and indirect observation of this unique system. We stress the need for new instruments, the value of synergistic Earth and Jupiter-based remote sensing and in-situ investigations, and the vital importance of multi-spacecraft, in-situ measurements. While simultaneous, multi-point in-situ observations have long become the standard for exploring electromagnetic interactions in the inner solar system, they have never taken place at Jupiter or any strongly magnetized planet besides Earth. We conclude that a dedicated multi-spacecraft mission to Jupiter's radiation belts is an essential and obvious way forward and deserves to be given a high priority in ESA's Voyage 2050 programme.

We report the discovery of an exceptional MIR flare in a Type 2 AGN, SDSS J165726.81+234528.1, at $z=0.059$. This object brightened by 3 mag in the Wide-field Infrared Survey Explorer (WISE) $W1$ and $W2$ bands between 2015 and 2017 (and is fading since 2018), without significant changes ($\lesssim$ 0.2 mag) in the optical over the same period of time. Based on the WISE light curves and near-IR imaging, the flare is more significant at longer wavelengths, suggesting an origin of hot dust emission. The estimated black hole mass ($\sim 10^{6.5}\,M_\odot$) from different methods places its peak bolometric luminosity around the Eddington limit. The high luminosity of the MIR flare and its multi-year timescale suggest that it most likely originated from reprocessed dust radiation in an extended torus surrounding the AGN, instead of from stellar explosions. The MIR color variability is consistent with known changing-look AGN and tidal disruption events (TDEs), but inconsistent with normal supernovae. We suggest that it is a turning-on Type 2 AGN or TDE, where the optical variability is obscured by the dust torus during the transition. This MIR flare event reveals a population of dramatic nuclear transients that are missed in the optical.

The Giant Radio Array for Neutrino Detection (GRAND) is a planned large-scale observatory of ultra-high-energy (UHE) cosmic particles, with energies exceeding 10^8 GeV. Its goal is to solve the long-standing mystery of the origin of UHE cosmic rays. To do this, GRAND will detect an unprecedented number of UHE cosmic rays and search for the undiscovered UHE neutrinos and gamma rays associated to them with unmatched sensitivity. GRAND will use large arrays of antennas to detect the radio emission coming from extensive air showers initiated by UHE particles in the atmosphere. Its design is modular: 20 separate, independent sub-arrays, each of 10 000 radio antennas deployed over 10 000 km^2. A staged construction plan will validate key detection techniques while achieving important science goals early. Here we present the science goals, detection strategy, preliminary design, performance goals, and construction plans for GRAND.

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 Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

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.

We present IRAM/NOEMA, JCMT/SCUBA-2 and VLA observations of the most distant known gravitationally lensed quasar J0439+1634 at z = 6.5. We detect strong dust emission, [CII] 158 $\mu$m, [CI] 369 $\mu$m, [OI] 146 $\mu$m, CO(6-5), CO(7-6), CO(9-8), CO(10-9), H$_{\rm 2}$O $3_{\rm 1,2}-2_{\rm 2,1}$, and H$_{\rm 2}$O $3_{\rm 2,1}-3_{\rm 1,2}$ lines as well as a weak radio continuum. The strong [CII] line yields a systemic redshift of the host galaxy to be z=6.5188$\pm$0.0001. The magnification makes J0439+1634 the far-infrared (FIR) brightest quasar at z > 6 known, with the brightest [CII] line yet detected at this redshift. The FIR luminosity is (3.4$\pm$0.2)$\times$10$^{13}$ $\mu^{-1}$ $L_{\odot}$, where $\mu$ $\sim$ 2.6 - 6.6 is the magnification of the host galaxy, estimated based on the lensing configuration from HST imaging. We estimate the dust mass to be (2.2$\pm$0.1)$\times$10$^{9}$ $\mu^{-1}$ $M_{\odot}$. The CO Spectral Line Energy Distribution using four CO lines is best fit by a two-component model of the molecular gas excitation. The estimates of molecular gas mass derived from CO lines and atomic carbon mass are consistent, in the range of 3.9 - 8.9 $\times$10$^{10} \mu^{-1}$ $M_{\odot}$. The [CII]/[CI], [CII]/CO, and [OI]/[CII] line luminosity ratios suggest a photodissociation region model with more than one component. The ratio of H$_{\rm 2}$O $3_{\rm 2,1}-3_{\rm 1,2}$ line luminosity to $L_{\rm TIR}$ is consistent with values in local and high redshift ultra-/hyper-luminous infrared galaxies. The VLA observations reveal an unresolved radio continuum source, and indicate that J0439+1634 is a radio quiet quasar with R = 0.05 - 0.17.

We present self-consistent radiation hydrodynamic simulations of hydrogen reionization performed with Arepo-RT complemented by a state-of-the-art galaxy formation model. We examine how photoheating feedback, due to reionization, shapes the galaxies properties. Our fiducial model completes reionization by $z\approx6$ and matches observations of the Ly$\alpha$ forest, the CMB electron scattering optical depth, the high-redshift UV luminosity function, and stellar mass function. Contrary to previous works, photoheating suppresses star formation rates by more than $50\%$ only in halos less massive than $\sim10^{8.4}\ M_\odot$ ($\sim10^{8.8}\ M_\odot$) at $z=6$ $(z=5)$, suggesting inefficient photoheating feedback from photons within galaxies. The use of a uniform UV background that heats up the gas at $z\approx10.7$ generates an earlier onset of suppression of star formation compared to our fiducial model. This discrepancy can be mitigated by adopting a UV background model with a more realistic reionization history. In the absence of stellar feedback, photoheating alone is only able to quench halos less massive than $\sim10^9\ M_\odot$ at $z\gtrsim5$, implying that photoheating feedback is sub-dominant in regulating star formation. In addition, stellar feedback, implemented as a non-local galactic wind scheme in the simulations, weakens the strength of photoheating feedback by reducing the amount of stellar sources. Most importantly, photoheating does not leave observable imprints in the UV luminosity function, stellar mass function, or the cosmic star formation rate density. The feasibility of using these observables to detect imprints of reionization therefore requires further investigation.

Among more than 200 quasars known at $z\gtrsim6$, only one object, J0100+2802 (z=6.327), was found hosting a $>10^{10}M_\odot$ super-massive black hole (SMBH). In order to investigate the host galaxy properties of J0100+2802, we performed multi-band ALMA observations, aiming at mapping the dust continuum, [CII] and CO(6-5) emission lines with sub-kiloparsec scale resolution, as well as detecting high-J CO lines in CO(11-10), CO(10-9), and CO(7-6). The galaxy size is measured to be $R_{\rm major}=3.6\pm0.2$ kpc from the high resolution continuum observations. No ordered motion on kilo-parsec scales was found in both [CII] and CO(6-5) emissions. The velocity dispersion is measured to be 161$\pm$7 km s$^{-1}$, $\sim$3 times smaller than that estimated from the local M-$\sigma$ relation. In addition, we found that the CO emission is more concentrate (a factor of 1.8$\pm$0.4) than the [CII] emission. Together with CO(2-1) detected by VLA, we measured the CO Spectral Line Energy Distribution (SLED), which is best fit by a two-components model, including a cool component at $\sim24$ K with a density of $n_{\rm (H_2)}$=10$^{4.5}$ cm$^{-3}$, and a warm component at $\sim224$ K with a density of $n_{\rm (H_2)}$=10$^{3.6}$ cm$^{-3}$. We also fitted the dust continuum with a graybody model, which indicates that it has either a high dust emissivity $\beta\gtrsim2$ or a hot dust temperature $T_{\rm dust}\gtrsim60$ K, or a combination of both factors. The highly excited CO emission and hot dust temperature suggest that the powerful AGN in J0100+2802 could contribute to the gas and dust heating although future observations are needed to confirm this.

The overwhelming foreground contamination is one of the primary impediments to probing the Epoch of Reionization (EoR) through measuring the redshifted 21 cm signal. Among various foreground components, radio halos are less studied and their impacts on the EoR observations are still poorly understood. In this work, we employ the Press-Schechter formalism, merger-induced turbulent re-acceleration model, and the latest SKA1-Low layout configuration to simulate the SKA "observed" images of radio halos. We calculate the one-dimensional power spectra from simulated images and find that radio halos can be about $10^4$, $10^3$ and $10^{2.5}$ times more luminous than the EoR signal on scales of $0.1\,\text{Mpc}^{-1} < k < 2\,\text{Mpc}^{-1}$ in the 120-128, 154-162, and 192-200 MHz bands, respectively. By examining the two-dimensional power spectra inside properly defined EoR windows, we find that the power leaked by radio halos can still be significant, as the power ratios of radio halos to the EoR signal on scales of $0.5\,\text{Mpc}^{-1} \lesssim k \lesssim 1\,\text{Mpc}^{-1}$ can be up to about 230-800%, 18-95%, and 7-40% in the three bands, when the 68% uncertainties caused by the variation of the number density of bright radio halos are considered. Furthermore, we find that radio halos located inside the far side-lobes of the station beam can also impose strong contamination within the EoR window. In conclusion, we argue that radio halos are severe foreground sources and need serious treatments in future EoR experiments.

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.

In this paper we report on the characterization of the Hamamatsu VUV4 (S/N: S13370-6152) Vacuum Ultra-Violet (VUV) sensitive Silicon Photo-Multipliers (SiPMs) as part of the development of a solution for the detection of liquid xenon scintillation light for the nEXO experiment. Various SiPM features, such as: dark noise, gain, correlated avalanches, direct crosstalk and Photon Detection Efficiency (PDE) were measured in a dedicated setup at TRIUMF. SiPMs were characterized in the range $163 \text{ } \text{K} \leq \text{T}\leq 233 \text{ } \text{K}$. At an over voltage of $3.1\pm0.2$ V and at $\text{T}=163 \text{ }\text{K}$ we report a number of Correlated Avalanches (CAs) per pulse in the $1 \upmu\text{s}$ interval following the trigger pulse of $0.161\pm0.005$. At the same settings the Dark-Noise (DN) rate is $0.137\pm0.002 \text{ Hz/mm}^{2}$. Both the number of CAs and the DN rate are within nEXO specifications. The PDE of the Hamamatsu VUV4 was measured for two different devices at $\text{T}=233 \text{ }\text{K}$ for a mean wavelength of $189\pm7\text{ nm}$. At $3.6\pm0.2$ V and $3.5\pm0.2$ V of over voltage we report a PDE of $13.4\pm2.6\text{ }\%$ and $11\pm2\%$, corresponding to a saturation PDE of $14.8\pm2.8\text{ }\%$ and $12.2\pm2.3\%$, respectively. Both values are well below the $24\text{ }\%$ saturation PDE advertised by Hamamatsu. More generally, the second device tested at $3.5\pm0.2$ V of over voltage is below the nEXO PDE requirement. The first one instead yields a PDE that is marginally close to meeting the nEXO specifications. This suggests that with modest improvements the Hamamatsu VUV4 MPPCs could be considered as an alternative to the FBK-LF SiPMs for the final design of the nEXO detector.

The ultra-diffuse dwarf galaxy NGC 1052-DF2 (DF2) has ten (eleven) measured globular clusters (GCs) with a line-of-sight velocity dispersion of $\sigma=7.8^{+5.2}_{-2.2}\,$km/s ($\sigma=10.6^{+3.9}_{-2.3}\,$km/s). Our conventional statistical analysis of the original ten GCs gives $\sigma=8.0^{+4.3}_{-3.0}\,$km/s. The overall distribution of velocities agrees well with a Gaussian of this width. Due to the non-linear Poisson equation in MOND, a dwarf galaxy has weaker self-gravity when in close proximity to a massive host. This external field effect is investigated using a new analytic formulation and fully self-consistent live $N$-body models in MOND. Our formulation agrees well with that of Famaey and McGaugh (2012). These new simulations confirm our analytic results and suggest that DF2 may be in a deep-freeze state unique to MOND. The correctly calculated MOND velocity dispersion agrees with our inferred dispersion and that of van Dokkum et al. (2018b) if DF2 is within 150 kpc of NGC 1052 and both are 20 Mpc away. The GCs of DF2 are however significantly brighter and larger than normal GCs, a problem which disappears if DF2 is significantly closer to us. A distance of 10-13 Mpc makes DF2 a normal dwarf galaxy even more consistent with MOND and the 13 Mpc distance reported by Trujillo et. al. (2019). We discuss the similar dwarf DF4, finding good agreement with MOND. We also discuss possible massive galaxies near DF2 and DF4 along with their distances and peculiar velocities, noting that NGC 1052 may lie at a distance near 10 Mpc.

As hybrids of narrow-line Seyfert 1 (NLS1) galaxies and blazars, {\gamma}-ray emitting NLS1s are important probes of jet physics in the high Eddington-ratio regime. Only very few of them are known to date; the majority of them below redshift z = 0.5. Here we present the identification of the {\gamma}-ray emitting AGN TXS 0943+105 (SDSS J094635.06+101706.1) as a high-redshift NLS1 galaxy. It turns out to be one of the radio-loudest NLS1s known, highly variable at all wavelengths, and shows widely extended radio emission at a (projected) > 100 kpc scale. It is a known strong {\gamma}-ray emitter with a luminous flare reported previously. At redshift z=1.004, this is the most distant {\gamma}-NLS1 known to date.

Space very long baseline interferometry (VLBI) has unique applications in high-resolution imaging of fine structure of astronomical objects and high-precision astrometry due to the key long space-Earth or space-space baselines beyond the Earth's diameter. China has been actively involved in the development of space VLBI in recent years. This review briefly summarizes China's research progress in space VLBI and the future development plan.

We study the kinematics of stars in clusters undergoing gas expulsion in standard Newtonian dynamics and also in Milgromian dynamics (MOND). Gas expulsion can explain the observed line-of-sight (LoS) velocity dispersion profile of NGC 2419 in Newtonian dynamics. For a given star formation efficiency (SFE), the shapes of the velocity dispersion profiles, which are normalised by the velocity dispersion at the projected half-mass radius, are almost indistinguishable for different SFE models in Newtonian dynamics. The velocity dispersion of a star cluster in the outer halo of a galaxy can indeed have a strong radial anisotropy in Newtonian dynamics after gas expulsion. MOND displays several different properties from Newtonian dynamics. In particular, the slope of the central velocity dispersion profile is less steep in MOND for the same SFE. Moreover, for a given SFE, more massive embedded cluster models result in more rapidly declining central velocity dispersion profiles for the final star clusters, while less massive embedded cluster models lead to flatter velocity dispersion profiles for the final products. The onset of the radial-orbit instability in post-gas-expulsion MOND models is discussed. SFEs as low as a few percent, typical of molecular clouds, lead to surviving ultra-diffuse objects. Gas expulsion alone is unlikely the physical mechanism for the observed velocity dispersion profile of NGC 2419 in MOND.

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.

All-Sky-ASTROGAM is a gamma-ray observatory operating in a broad energy range, 100 keV to a few hundred MeV, recently proposed as the "Fast" (F) mission of the European Space Agency for a launch in 2028 to an L2 orbit. The scientific payload is composed of a unique gamma-ray imaging monitor for astrophysical transients, with very large field of view (almost 4$\pi$ sr) and optimal sensitivity to detect bright and intermediate flux sources (gamma-ray bursts, active galactic nuclei, X-ray binaries, supernovae and novae) at different timescales ranging from seconds to months. The mission will operate in a maturing gravitational wave and multi-messenger epoch, opening up new and exciting synergies.

We present a catalog of 290 `winged' or X-shaped radio galaxies (XRGs) extracted from the latest (17 Dec. 2014) data release of the `Very Large Array Faint Images of the Radio Sky at Twenty centimeter' (VLA FIRST survey). We have combined these radio images with their counterparts in the TIFR GMRT sky survey (TGSS) at 150 MHz (TGSS alternative data release, TGSS_ADR1), in an attempt to identify any low surface-brightness radio emission present in these sources. This has enabled us to assemble a sample of 106 `strong' XRG candidates and 184 `probable' XRG candidates whose XRG designation must be verified by further observations. The present sample of 290 XRG candidates is almost twice as large as the number of XRGs presently known. Twenty-five of our 290 XRG candidates (9 `strong' and 16 `probable') are identified as quasars. Double-peaked narrow emission lines are seen in the optical spectra of three of the XRG candidates (2 `strong' and 1 `probable'). Nearly 90% of the sample is located in the FR II domain of the Owen-Ledlow diagram. A few of the strong XRG candidates have a rather flat radio spectrum (spectral index $\alpha$ flatter than -0.3) between 150 MHz and 1.4 GHz, or between 1.4 GHz and 5 GHz. Since this is not expected for lobe-dominated extragalactic radio sources (like nearly all known XRGs), these sources are particularly suited for follow-up radio imaging and near-simultaneous measurement of the radio spectrum.

A central compact object (CCO, e.g. a black hole) with an accretion disk has been suggested as the common central engine of various astrophysical phenomena, such as gamma-ray bursts (GRBs), tidal disruption events (TDEs) and active galactic nuclei (AGNs). A jet powered by such a system might precess due to the misalignment of the angular momenta of the CCO and accretion disk. Some quasi-periodic behaviors observed in the light curves of these phenomena can be well interpreted within the framework of a precessing jet model. In this paper, we study the emission polarization of precessing jets in the three kinds of phenomena. The polarization angle also shows a gradual change for the synchrotron emission in both the random and toroidal magnetic field configurations with the precessing jet, while it can only change abruptly by $90^\circ$ for the non-precessing top-hat jet. Polarization properties are periodic due to the assumptions made in our model. The polarization observations are crucial to confirm the precession nature of jets in GRBs, TDEs and AGNs.

w-Projection is a wide-field imaging technique that is widely used in radio synthesis arrays. Processing the wide-field big data generated by the future Square Kilometre Array (SKA) will require significant updates to current methods to significantly reduce the time consumed on data processing. Data loading and gridding are found to be two major time-consuming tasks in w-projection. In this paper, we investigate two parallel methods of accelerating w-projection processing on multiple nodes: the hybrid Message Passing Interface (MPI) and Open Multi-Processing (OpenMP) method based on multicore Central Processing Units (CPUs) and the hybrid MPI and Compute Unified Device Architecture (CUDA) method based on Graphics Processing Units (GPUs). Both methods are successfully employed and operated in various computational environments, confirming their robustness. The experimental results show that the total runtime of both MPI + OpenMP and MPI + CUDA methods is significantly shorter than that of single-thread processing. MPI + CUDA generally shows faster performance when running on multiple nodes than MPI + OpenMP, especially on large numbers of nodes. The single-precision GPU-based processing yields faster computation than the double-precision processing; while the single- and doubleprecision CPU-based processing shows consistent computational performance. The gridding time remarkably increases when the support size of the convolution kernel is larger than 8 and the image size is larger than 2,048 pixels. The present research offers useful guidance for developing SKA imaging pipelines.

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.

Neutron star-neutron star mergers are known to be associated with short gamma-ray bursts. If the neutron star equation of state is sufficiently stiff, at least some of such mergers will leave behind a supramassive or even a stable neutron star that spins rapidly with a strong magnetic field (i.e., a magnetar). Such a magnetar signature may have been observed as the X-ray plateau following a good fraction (up to 50%) of short gamma-ray bursts, and it has been expected that one may observe short gamma-ray burst-less X-ray transients powered by double neutron star mergers. A fast X-ray transient (CDF-S XT1) was recently found to be associated with a faint host galaxy whose redshift is unknown. Its X-ray and host-galaxy properties allow several possibleexplanations including a short gamma-ray burst seen off axis, a low-luminosity gamma-ray burst at high redshift, or a tidal disruption event involving an intermediate mass black hole and a white dwarf. Here we report a second X-ray transient, CDF-S XT2, that is associated with a galaxy at redshift z = 0.738. The light curve is fully consistent with being powered by a millisecond magnetar. More intriguingly, CDF-S XT2 lies in the outskirts of its star-forming host galaxy with a moderate offset from the galaxy center, as short bursts often do. The estimated event rate density of similar X-ray transients, when corrected to the local value, is consistent with the double neutron star merger rate density inferred from the detection of GW170817.

A great challenge in present-day physics is to understand whether the observed internal dynamics of galaxies is due to dark matter matter or due to a modification of the law of gravity. Recently, van Dokkum et al. reported that the ultra-diffuse dwarf galaxy NGC1052-DF2 lacks dark matter, and they claimed that this would -- paradoxically -- be problematic for modified gravity theories like Milgromian dynamics (MOND). However, NGC1052-DF2 is not isolated, so that a valid prediction of its internal dynamics in MOND cannot be made without properly accounting for the external gravitational fields from neighbouring galaxies. Including this external field effect following Haghi et al. shows that NGC1052-DF2 is consistent with MOND.

TREND50 is a radio detection setup of 50 self-triggered antennas working in the 50-100MHz frequency range and deployed in a radio-quiet valley of the Tianshan mountains (China). TREND50 achieved its goal: the autonomous radiodetection and identification of air showers. Thanks to a dedicated offine selection algorithm, 564 air shower candidates were indeed selected out of $7\cdot10^8$ transient radio signals recorded during the 314 live days of data taken during the first two years of operation of this setup (2011 and 2012). This event rate, as well as the distribution of the candidate directions of arrival, is consistent with what is expected from cosmic-ray-induced air showers according to simulations, assuming an additional $\sim$20% contamination of the final sample by background events. This result is obtained at the cost of a reduced air shower detection efficiency, estimated to be $\sim$3%. This low efficiency is mostly due to the large amount of dead time of the setup. This result paves the way for the GRANDProto35 experiment, the first stage of the GRAND project.

We present optical luminosity functions (LFs) of galaxies in the $^{0.1}g$, $^{0.1}r$, and $^{0.1}i$ bands, calculated using data in $\sim 40$ $deg^{2}$ sky area of LAMOST Complete Spectroscopic Survey of Pointing Area (LaCoSSPAr) in Southern Galactic Cap. Redshifts for galaxies brighter $r = 18.1$ were obtained mainly with LAMOST. In each band, LFs derived using both parametric and non-parametric maximum likelihood methods agree well with each other. In the $^{0.1}r$ band, our fitting parameters of the Schechter function are $\phi_{*}=(1.65\pm0.36)\times10^{-2}h^{3}Mpc^{-3}$, $M_{*}=-20.69\pm0.06$ mag, and $\alpha=-1.12\pm0.08$, in agreements with previous studies. Separate LFs are also derived for emission line galaxies and absorption line galaxies, respectively. The LFs of absorption line galaxies show a dip at $^{0.1}r \sim 18.5$ and can be well fitted by a double-Gaussian function, suggesting a bi-modality in passive galaxies.

We aim to select quasar candidates based on the two large survey databases, Pan-STARRS and AllWISE. Exploring the distribution of quasars and stars in the color spaces, we find that the combination of infrared and optical photometry is more conducive to select quasar candidates. Two new color criterions (yW1W2 and izW1W2) are constructed to distinguish quasars from stars efficiently. With izW1W2, 98.30% of star contamination is eliminated, while 99.50% of quasars are retained, at least to the magnitude limit of our training set of stars. Based on the optical and infrared color features, we put forward an efficient schema to select quasar candidates and high redshift quasar candidates, in which two machine learning algorithms (XGBoost and SVM) are implemented. The XGBoost and SVM classifiers have proven to be very effective with accuracy of 99.46% when 8Color as input pattern and default model parameters. Applying the two optimal classifiers to the unknown Pan-STARRS and AllWISE cross-matched data set, a total of 2,006,632 intersected sources are predicted to be quasar candidates given quasar probability larger than 0.5 (i.e. P_QSO>0.5). Among them, 1,201,211 have high probability (P_QSO>0.95). For these newly predicted quasar candidates, a regressor is constructed to estimate their redshifts. Finally 7,402 z>3.5 quasars are obtained. Given the magnitude limitation and site of the LAMOST telescope, part of these candidates will be used as the input catalogue of the LAMOST telescope for follow-up observation, and the rest may be observed by other telescopes.

When applying the foreground removal methods to uncover the faint cosmological signal from the epoch of reionization (EoR), the foreground spectra are assumed to be smooth. However, this assumption can be seriously violated in practice since the unresolved or mis-subtracted foreground sources, which are further complicated by the frequency-dependent beam effects of interferometers, will generate significant fluctuations along the frequency dimension. To address this issue, we propose a novel deep-learning-based method that uses a 9-layer convolutional denoising autoencoder (CDAE) to separate the EoR signal. After being trained on the SKA images simulated with realistic beam effects, the CDAE achieves excellent performance as the mean correlation coefficient ($\bar{\rho}$) between the reconstructed and input EoR signals reaches $0.929 \pm 0.045$. In comparison, the two representative traditional methods, namely the polynomial fitting method and the continuous wavelet transform method, both have difficulties in modelling and removing the foreground emission complicated with the beam effects, yielding only $\bar{\rho}_{\text{poly}} = 0.296 \pm 0.121$ and $\bar{\rho}_{\text{cwt}} = 0.198 \pm 0.160$, respectively. We conclude that, by hierarchically learning sophisticated features through multiple convolutional layers, the CDAE is a powerful tool that can be used to overcome the complicated beam effects and accurately separate the EoR signal. Our results also exhibit the great potential of deep-learning-based methods in future EoR experiments.

We report the discovery of a new fast radio burst (FRB), FRB~010312, in archival data from a 1.4\,GHz survey of the Magellanic Clouds using the multibeam receiver on the Parkes 64\,m-diameter radio telescope. These data sets include the Lorimer burst (FRB~010724), which it pre-dates and which we also re-detect. The new burst has a much higher dispersion measure of 1187\,cm$^{-3}$pc. The burst is one of the broadest found to date, the second earliest FRB known, and the ninth FRB discovered with a dispersion measure larger than 1000\,cm$^{-3}\,$pc. Our discovery indicates that there are likely to be more burst events still to be found in the existing Parkes data archive.

We present a morphological classification of 14,245 radio active galactic nuclei (AGNs) into six types, i.e., typical Fanaroff--Riley Class I / II (FRI/II), FRI/II-like bent-tailed, X-shaped radio galaxy, and ringlike radio galaxy, by designing a convolutional neural network (CNN) based autoencoder, namely MCRGNet, and applying it to a labeled radio galaxy (LRG) sample containing 1442 AGNs and an unlabeled radio galaxy (unLRG) sample containing 14,245 unlabeled AGNs selected from the Best--Heckman sample. We train MCRGNet and implement the classification task by a three-step strategy, i.e., pre-training, fine-tuning, and classification, which combines both unsupervised and supervised learnings. A four-layer dichotomous tree is designed to classify the radio AGNs, which leads to a significantly better performance than the direct six-type classification. On the LRG sample, our MCRGNet achieves a total precision of $\sim 93\%$ and an averaged sensitivity of $\sim 87\%$, which are better than those obtained in previous works. On the unLRG sample, whose labels have been human-inspected, the neural network achieves a total precision of $\sim 80\%$. Also, using the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) to calculate the $r$-band absolute magnitude ($M_\mathrm{opt}$) and using the flux densities to calculate the radio luminosity ($L_\mathrm{radio}$), we find that the distributions of the unLRG sources on the $L_\mathrm{radio}$--$M_\mathrm{opt}$ plane do not show an apparent redshift evolution and could confirm with a sufficiently large sample that there could not exist an abrupt separation between FRIs and FRIIs as reported in some previous works.

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 present integrated pulse profiles at 5~GHz for 71 pulsars, including eight millisecond pulsars (MSPs), obtained using the Shanghai Tian Ma Radio Telescope (TMRT). Mean flux densities and pulse widths are measured. For 19 normal pulsars and one MSP, these are the first detections at 5~GHz and for a further 19, including five MPSs, the profiles have a better signal-to-noise ratio than previous observations. Mean flux density spectra between 400~MHz and 9~GHz are presented for 27 pulsars and correlations of power-law spectral index are found with characteristic age, radio pseudo-luminosity and spin-down luminosity. Mode changing was detected in five pulsars. The separation between the main pulse and interpulse is shown to be frequency independent for six pulsars but a frequency dependence of the relative intensity of the main pulse and interpulse is found. The frequency dependence of component separations is investigated for 20 pulsars and three groups are found: in seven cases the separation between the outmost leading and trailing components decreases with frequency, roughly in agreement with radius-to-frequency mapping; in eleven cases the separation is nearly constant; in the remain two cases the separation between the outmost components increases with frequency. We obtain the correlations of pulse widths with pulsar period and estimate the core widths of 23 multi-component profiles and conal widths of 17 multi-component profiles at 5.0~GHz using Gaussian fitting and discuss the width-period relationship at 5~GHz compared with the results at at 1.0~GHz and 8.6~GHz.

We report initial results from a large Gemini program to observe z>~5.7 quasars with GNIRS near-IR spectroscopy. Our sample includes 50 quasars with simultaneous ~0.85-2.5 micron spectra covering the rest-frame ultraviolet and major broad emission lines from Ly-alpha to MgII. We present spectral measurements for these quasars and compare to their lower-redshift counterparts at z=1.5-2.3. We find that when quasar luminosity is matched, there are no significant differences between the rest-UV spectra of z>~5.7 quasars and the low-z comparison sample. High-z quasars have similar continuum and emission line properties and occupy the same region in the black hole mass and luminosity space as the comparison sample, accreting at an average Eddington ratio of ~0.3. There is no evidence for super-Eddington accretion or hypermassive (>10^10 Msun) black holes within our sample. We find a mild excess of quasars with weak CIV lines relative to the control sample. Our results, corroborating earlier studies but with better statistics, demonstrate that these high-z quasars are already mature systems of accreting supermassive black holes operating with the same physical mechanisms as those at lower redshifts.

PAN is a scientific instrument suitable for deep space and interplanetary missions. It can precisely measure and monitor the flux, composition, and direction of highly penetrating particles ($> \sim$100 MeV/nucleon) in deep space, over at least one full solar cycle (~11 years). The science program of PAN is multi- and cross-disciplinary, covering cosmic ray physics, solar physics, space weather and space travel. PAN will fill an observation gap of galactic cosmic rays in the GeV region, and provide precise information of the spectrum, composition and emission time of energetic particle originated from the Sun. The precise measurement and monitoring of the energetic particles is also a unique contribution to space weather studies. PAN will map the flux and composition of penetrating particles, which cannot be shielded effectively, precisely and continuously, providing valuable input for the assessment of the related health risk, and for the development of an adequate mitigation strategy. PAN has the potential to become a standard on-board instrument for deep space human travel. PAN is based on the proven detection principle of a magnetic spectrometer, but with novel layout and detection concept. It will adopt advanced particle detection technologies and industrial processes optimized for deep space application. The device will require limited mass (~20 kg) and power (~20 W) budget. Dipole magnet sectors built from high field permanent magnet Halbach arrays, instrumented in a modular fashion with high resolution silicon strip detectors, allow to reach an energy resolution better than 10\% for nuclei from H to Fe at 1 GeV/n.

Gamma-ray bursts are the strongest explosions in the Universe since the Big Bang, believed to be produced either in forming black holes at the end of massive star evolution or merging of compact objects. Spectral and timing properties of gamma-ray bursts suggest that the observed bright gamma-rays are produced in the most relativistic jets in the Universe; however, the physical properties, especially the structure and magnetic topologies in the jets are still not well known, despite several decades of studies. It is widely believed that precise measurements of the polarization properties of gamma-ray bursts should provide crucial information on the highly relativistic jets. As a result there have been many reports of gamma-ray burst polarization measurements with diverse results, see, however many such measurements suffered from substantial uncertainties, mostly systematic. After the first successful measurements by the GAP and COSI instruments, here we report a statistically meaningful sample of precise polarization measurements, obtained with the dedicated gamma-ray burst polarimeter, POLAR onboard China's Tiangong-2 spacelab. Our results suggest that the gamma-ray emission is at most polarized at a level lower than some popular models have predicted; although our results also show intrapulse evolution of the polarization angle. This indicates that the low polarization degrees could be due to an evolving polarization angle during a gamma-ray burst.

Gamma-ray bursts are the strongest explosions in the Universe since the Big Bang, believed to be produced either in forming black holes at the end of massive star evolution or merging of compact objects. Spectral and timing properties of gamma-ray bursts suggest that the observed bright gamma-rays are produced in the most relativistic jets in the Universe; however, the physical properties, especially the structure and magnetic topologies in the jets are still not well known, despite several decades of studies. It is widely believed that precise measurements of the polarization properties of gamma-ray bursts should provide crucial information on the highly relativistic jets. As a result there have been many reports of gamma-ray burst polarization measurements with diverse results, see, however many such measurements suffered from substantial uncertainties, mostly systematic. After the first successful measurements by the GAP and COSI instruments, here we report a statistically meaningful sample of precise polarization measurements, obtained with the dedicated gamma-ray burst polarimeter, POLAR onboard China's Tiangong-2 spacelab. Our results suggest that the gamma-ray emission is at most polarized at a level lower than some popular models have predicted; although our results also show intrapulse evolution of the polarization angle. This indicates that the low polarization degrees could be due to an evolving polarization angle during a gamma-ray burst.

We report first results from the CHinese Exoplanet Searching Program from Antarctica (CHESPA)---a wide-field high-resolution photometric survey for transiting exoplanets carried out using telescopes of the AST3 (Antarctic Survey Telescopes times 3) project. There are now three telescopes (AST3-I, AST3-II, and CSTAR-II) operating at Dome A---the highest point on the Antarctic Plateau---in a fully automatic and remote mode to exploit the superb observing conditions of the site, and its long and uninterrupted polar nights. The search for transiting exoplanets is one of the key projects for AST3. During the Austral winters of 2016 and 2017 we used the AST3-II telescope to survey a set of target fields near the southern ecliptic pole, falling within the continuous viewing zone of the TESS mission \citep{Ricker10}. The first data release of the 2016 data, including images, catalogs and lightcurves of 26578 bright stars ($7.5\le i \le15$) was presented in \citet{Zhang18}. The best precision, as measured by the RMS of the lightcurves at the optimum magnitude of the survey ($i=10$), is around 2\,mmag. We detect 222 objects with plausible transit signals from these data, 116 of which are plausible transiting exoplanet candidates according to their stellar properties as given by the TESS Input Catalog \citep{Stassun17}, Gaia DR2 \citep{Gaia18} and TESS-HERMES spectroscopy \citep{Sharma18}. With the first data release from TESS expected in late 2018, this candidate list will be a timely for improving the rejection of potential false positives.

Located at Dome A, the highest point of the Antarctic plateau, the Chinese Kunlun station is considered to be one of the best ground-based photometric sites because of its extremely cold, dry, and stable atmosphere(Saunders et al. 2009). A target can be monitored from there for over 40 days without diurnal interruption during a polar winter. This makes Kunlun station a perfect site to search for short-period transiting exoplanets. Since 2008, an observatory has been built at Kunlun station and three telescopes are working there. Using these telescopes, the AST3 project has been carried out over the last six years with a search for transiting exoplanets as one of its key programs (CHESPA). In the austral winters of 2016 and 2017, a set of target fields in the Southern CVZ of TESS (Ricker et al. 2009) were monitored by the AST3-II telescope. In this paper, we introduce the CHESPA and present the first data release containing photometry of 26,578 bright stars (m_i < 15). The best photometric precision at the optimum magnitude for the survey is around 2 mmag. To demonstrate the data quality, we also present a catalog of 221 variables with a brightness variation greater than 5 mmag from the 2016 data. Among these variables, 179 are newly identified periodic variables not listed in the AAVSO databasea), and 67 are listed in the Candidate Target List(Stassun et al. 2017). These variables will require careful attention to avoid false-positive signals when searching for transiting exoplanets. Dozens of new transiting exoplanet candidates will be also released in a subsequent paper(Zhang et al. 2018b).

In this paper we propose the model that the coalescence of primordial black holes (PBHs) binaries with equal mass $M \sim 10^{28}$g can emit luminous gigahertz (GHz) radio transient, which may be candidate sources for the observed fast radio bursts (FRBs), if at least one black hole holds appropriate amount of net electric charge $Q$. Using a dimensionless quantity for the charge $q = Q/\sqrt{G}M$, our analyses infer that $q\sim O(10^{-4.5})$ can explain the FRBs with released energy of order $O(10^{40}) {\rm ergs}$. With the current sample of FRBs and assuming a distribution of charge $\phi(q)$ for all PBHs, we can deduce that its form is proportional to $q^{-3.0\pm0.1}$ for $q\geq 7.2\times10^{-5}$ if PBHs are sources of the observed FRBs. Furthermore, with the proposed hypothetical scenario and by estimating the local event rate of FRBs $\sim 2.6 \times 10^3 {\rm Gpc}^{-3} {\rm yr}^{-1}$, one derives a lower bound for the fraction of PBHs (at the mass of $10^{28}$g) against that of matter $f_{\rm PBH}(10^{28}{\rm g})$ $\gtrsim 10^{-5}$ needed to explain the rate. With this inspiring estimate, we expect that future observations of FRBs can help to falsify their physical origins from the PBH binaries coalescences. In the future, the gravitational waves produced by mergers of small black holes can be detected by high frequency gravitational wave detectors. We believe that this work would be a useful addition to the current literature on multimessenger astronomy and cosmology.

The central engine of gamma-ray bursts (GRBs) is poorly constrained. There exist two main candidates: a fast-rotating black hole and a rapidly spinning magnetar. Furthermore, the X-ray plateaus are widely accepted by the energy injection into the external shock. In this paper, we systematically analyze the \emph{Swift}/XRT light curves of 101 GRBs having plateau phases and known redshifts (before May 2017). Since a maximum energy budget ($\sim2\times10^{52}$ erg) exists for magnetars but not for a black hole, this provides a good clue to identify the type of GRB central engine. We calculate the isotropic kinetic energy $E_{\rm K,iso}$ and the isotropic X-ray energy release $E_{\rm X,iso}$ for individual GRB. We identify three categories based on how likely a black hole harbor at central engine: 'Gold' (9 out of 101, both $E_{\rm X,iso}$ and $E_{\rm K,iso}$ exceed the energy budget), 'Silver' (69 out of 101, $E_{\rm X,iso}$ less than the limit but $E_{\rm K,iso}$ is greater than the limit), and 'Bronze' (23 out of 101, the energies are not above the limit). We then derive and test the black hole parameters with the Blandford Znajek mechanism, and find that the observations of the black hole candidate ('Gold'+'Silver') samples are consistent with the expectations of the black hole model. Furthermore, we also test the magnetar candidate ('Bronze') sample with the magnetar model, and find that the magnetar surface magnetic led ($B_{p}$) and initial spin period ($P_{0}$) are consistent with the expectations of the magnetar model. Our analysis indicates that, if the magnetar wind is isotropic, a magnetar central engine is possible for 20\% of the analyzed GRBs. For most GRBs a black hole is most likely operating.

The spectrum of gamma-ray burst (GRB) afterglows can be studied with color indices. Here we present a large comprehensive catalogue of 70 GRBs with multi-wavelength optical transient data on which we perform a systematic study to find the temporal evolution of color indices. We categorize them into two samples based on how well the color indices are evaluated. The Golden sample includes 25 bursts mostly observed by GROND, and the Silver sample includes 45 bursts observed by other telescopes. For the Golden sample, we find that 95\% of the color indices do not vary over time. However, the color indices do vary during short periods in most bursts. The observed variations are consistent with effects of (i) the cooling frequency crossing the studied energy bands in a wind medium (43\%) and in a constant density medium (30\%), (ii) early dust extinction (12\%), (iii) transition from reverse shock to forward shock emission (5\%), or (iv) an emergent supernova emission (10\%). We also study the evolutionary properties of the mean color indices for different emission episodes. We find that 86\% of the color indices in the 70 bursts show constancy between consecutive ones. The color index variations occur mainly during the late GRB-SN bump, the flare and early reversed-shock emission components. We further perform a statistical analysis of various observational properties and model parameters (spectral index $\beta_{o}^{\rm CI}$, electron spectral indices $p^{\rm CI}$, etc.) using color indices. Overall, we conclude that $\sim$ 90\% of colors are constant in time and can be accounted for by the simplest external forward shock model, while the varying color indices call for more detailed modeling.

In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

In this paper we present the enhanced X-ray Timing and Polarimetry mission - eXTP. eXTP is a space science mission designed to study fundamental physics under extreme conditions of density, gravity and magnetism. The mission aims at determining the equation of state of matter at supra-nuclear density, measuring effects of QED, and understanding the dynamics of matter in strong-field gravity. In addition to investigating fundamental physics, eXTP will be a very powerful observatory for astrophysics that will provide observations of unprecedented quality on a variety of galactic and extragalactic objects. In particular, its wide field monitoring capabilities will be highly instrumental to detect the electro-magnetic counterparts of gravitational wave sources. The paper provides a detailed description of: (1) the technological and technical aspects, and the expected performance of the instruments of the scientific payload; (2) the elements and functions of the mission, from the spacecraft to the ground segment.

The Square Kilometre Array (SKA) will be the first low-frequency instrument with the capability to directly image the structures of the Epoch of Reionization (EoR). Indeed, deep imaging of the EoR over 5 targeted fields of 20 square degrees each has been selected as the highest priority science objective for SKA1, probably the most exciting scientific outcome to be expected for the SKA1-Low. Aiming at preparing for this highly challenging observation, we perform an extensive pre-selection of the 'quietest' and 'cleanest' candidate fields in the southern sky to be suited for deep imaging of the EoR using existing catalogs and observations over a broad band. The candidate fields should meet a number of strict criteria to avoid contamination from foreground structures and sources such as the Galactic plane, LMC/SMC, clusters of galaxies and bright radio sources. The candidate fields should also exhibit both the lowest average surface brightness and smallest variance to ensure uniformity and high quality deep imaging over the fields. We use two sets of field of view in our selection: the exactly same size of 20 square degrees as defined in current SKA1 documentation and 40 square degrees - an area of twice big as the EoR imaging fields. The choice of a larger area allows us to maximally reduce the environmental effect from radio structures/sources in the vicinity of the candidate field. Our selection yields a sample of 29 (9) `ideal' fields of 20 (40) square degrees in the southern sky that could be targeted for deep imaging of the EoR. We have launched a campaign to observe these candidate fields with some of the SKA pathfinders to set further constraints on these candidates. These pre-selected candidate fields can provide a valuable guidance for eventual selection of the five targeted fields for imaging the EoR structures with the future SKA1-low.

We report a synchronized kinematic shift of CIV and SiIV broad absorption lines (BAL) in a high-ionization, radio-loud, and X-ray bright quasar SDSS-J092345+512710 (at $z_{em} \sim 2.1627$). This quasar shows two broad absorption components (blue component at $v \sim 14,000\ km\ s^{-1}$, and red component at $v \sim 4,000\ km\ s^{-1}$ with respect to the quasars systemic redshift). The absorption profiles of CIV and SiIV BAL of the blue component show decrease in outflow velocity with an average deceleration rate of $-1.62_{-0.05}^{+0.04}\ cm\ s^{-2}$ and $-1.14^{+0.21}_{-0.22}\ cm\ s^{-2}$ over a rest-frame time-span of 4.15 yr. We do not see any acceleration-like signature in the red component. This is consistent with dramatic variabilities usually seen at high velocities. During our monitoring period the quasar has shown no strong continuum variability. We suggest the observed variability could be related to the time dependent changes in disk wind parameters like launching radius, initial flow velocity or mass outflow rate.

A large-scale ordered magnetic field plays a very important role in the formation and acceleration of a gamma-ray burst (GRB) jet. During the GRB prompt phase, some dissipation processes may happen and disturb the magnetic field, making the field become random to some extent. Therefore, a mixed magnetic field consisting of an ordered component and a random component is plausible for the early afterglow phase. Here we investigate the polarization evolution and light curve of an afterglow under such a 3-dimensional mixed magnetic field. Three kinds of ordered component (i.e., aligned, toroidal and radial) are discussed. We find that the three cases are distinguishable through jet polarization evolution. The polarization angle for a 3D mixed magnetic field with an aligned ordered component can evolve gradually but only changes abruptly by $90^\circ$ in the toroidal and radial cases. Furthermore, during the reverse shock crossing time, the polarization degree (PD) can be non-zero for the toroidal case but roughly zero for the radial case. Since an aligned component in a jet corresponds to a magnetar central engine and a toroidal component corresponds to a black hole, GRB central engines are distinguishable through polarization observations even if the magnetic field is mixed in a jet. In addition, our polarization calculation can be applied to GRB 170817A associated with GW170817. Using the recently-observed PD upper limit $12\%$ of GRB 170817A at $t=244$ days, the magnetic field strength ratio of the ordered to random components in this afterglow is constrained to be $\lesssim0.9$.

This is the fourth paper in a series of publications aiming at discovering quasars at the epoch of reionization. In this paper, we expand our search for $z\sim 7$ quasars to the footprint of the Dark Energy Survey (DES) Data Release One (DR1), covering $\sim 5000$ deg$^2$ of new area. We select $z\sim 7$ quasar candidates using deep optical, near-infrared (near-IR) and mid-IR photometric data from the DES DR1, the VISTA Hemisphere Survey (VHS), the VISTA Kilo-degree Infrared Galaxy (VIKING) survey, the UKIRT InfraRed Deep Sky Surveys -- Large Area Survey (ULAS) and the unblurred coadds from the Wide-field Infrared Survey Explore ($WISE$) images (unWISE). The inclusion of DES and unWISE photometry allows the search to reach $\sim$ 1 magnitude fainter, comparing to our $z \gtrsim 6.5$ quasar survey in the northern sky (Wang et al. 2018). We report the initial discovery and spectroscopic confirmation of six new luminous quasars at $z>6.4$, including an object at $z=7.02$, the fourth quasar yet known at $z>7$, from a small fraction of candidates observed thus far. Based on the recent measurement of $z \sim 6.7 $ quasar luminosity function using the quasar sample from our survey in the northern sky, we estimate that there will be $\gtrsim$ 55 quasars at $z > 6.5$ at $M_{1450} < -24.5$ in the full DES footprint.

A central problem of contemporary physics is whether the law of gravity is non-Newtonian on galaxy scales. Rodrigues et al. argue that Milgromian gravitation, which solves the flat rotation curve problem without the need for dark matter particles, is ruled out at > 10{\sigma} significance. To a large extent, this conclusion relies on galaxies with very uncertain distances and/or nearly edge-on orientations, where dust obscuration often becomes significant. Applying appropriate quality cuts to the data leaves only a handful of outliers to the predictions of Milgromian gravitation according to the analysis of Rodrigues et al., but even these outliers can be explained with Milgromian gravitation.

One large glitch was detected in PSR B1737$-$30 using data spanning from MJD 57999 to 58406 obtained with the newly built Shanghai Tian Ma Radio Telescope (TMRT). The glitch took place at the time around MJD 58232.4 when the pulsar underwent an increase in the rotation frequency of $\Delta \nu$ about 1.38$\times 10^{-6}$ Hz, corresponding to a fractional step change of $\Delta \nu / \nu$ $\thicksim$ 8.39$\times 10^{-7}$. Post$\textrm{-}$glitch $\nu$ gradually decreased to the pre$\textrm{-}$glitch value. The frequency derivative was observed to undergo a step change of about $-$9$\times 10^{-16}$ s$^{-2}$. Since July 1987, there are 36 glitches already reported in PSR B1737$-$30 including this one. According to our analysis, the glitch size distribution is well described by the power law with index of 1.13. The distribution of the interval between two adjacent glitches (waiting time $\Delta T$) follows a Poissonian probability density function. For PSR B1737$-$30, the interval is prone to be long after a large glitch. But no correlation is found between glitch size and the interval since previous glitch.

Despite extensive efforts, only two quasars have been found at $z>7$ to date due to a combination of low spatial density and high contamination from more ubiquitous Galactic cool dwarfs in quasar selection. This limits our current knowledge of the super-massive black hole (SMBH) growth mechanism and reionization history. In this letter, we report the discovery of a luminous quasar at $z=7.021$, DELS J003836.10$-$152723.6 (hereafter J0038$-$1527), selected using photometric data from DESI Legacy imaging Survey (DELS), Pan-STARRS1 (PS1) imaging Survey, as well as Wide-field Infrared Survey Explore ($WISE$) mid-infrared all-sky survey. With an absolute magnitude of $M_{1450}$=$-$27.1 and bolometric luminosity of $L_{\rm Bol}$=5.6$\times$10$^{13}$ $L_\odot$, J0038$-$1527 is the most luminous quasar known at $z>7$. Deep optical to near infrared spectroscopic observations suggest that J0038-1527 hosts a 1.3 billion solar mass BH accreting at the Eddington limit, with an Eddington ratio of 1.25$\pm$0.19. The CIV broad emission line of J0038$-$1527 is blue-shifted by more than 3000 km s$^{-1}$ to the systemic redshift. More detailed investigations of the high quality spectra reveal three extremely high velocity CIV broad absorption lines (BALs) with velocity from 0.08 to 0.14 times the speed of light and total balnicity index of more than 5000 km s$^{-1}$, suggesting the presence of relativistic outflows. J0038$-$1527 is the first quasar found at the epoch of reionization (EoR) with such strong outflows and provides a unique laboratory to investigate AGN feedback on the formation and growth of the most massive galaxies in the early universe.

We choose the bright compact group HCG 62, which was found to exhibit both excess X-ray emission and high Fe abundance to the southwest of its core, as an example to study the impact of mergers on chemical enrichment in the intragroup medium. We first reanalyze the high-quality Chandra and XMM-Newton archive data to search for the evidence for additional SN II yields, which is expected as a direct result of the possible merger-induced starburst. We reveal that, similar to the Fe abundance, the Mg abundance also shows a high value in both the innermost region and the southwest substructure, forming a high-abundance plateau, meanwhile all the SN Ia and SN II yields show rather flat distributions in $>0.1r_{200}$ in favor of an early enrichment. Then we carry out a series of idealized numerical simulations to model the collision of two initially isolated galaxy groups by using the TreePM-SPH GADGET-3 code. We find that the observed X-ray emission and metal distributions, as well as the relative positions of the two bright central galaxies with reference to the X-ray peak, can be well reproduced in a major merger with a mass ratio of 3 when the merger-induced starburst is assumed. The `best-match' snapshot is pinpointed after the third pericentric passage when the southwest substructure is formed due to gas sloshing. By following the evolution of the simulated merging system, we conclude that the effects of such a major merger on chemical enrichment are mostly restricted within the core region when the final relaxed state is reached.

We present the Data Release 4&5 quasar catalog from the quasar survey by Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST), which includes quasars observed between September 2015 and June 2017. There are a total of 19,253 quasars identified by visual inspections of the spectra. Among them, 11,458 are independently discovered by LAMOST, in which 3296 were reported by SDSS DR12 and DR14 quasar catalog after our survey began, while the rest 8162 are new discoveries of LAMOST. We provide the emission line measurements for the Halpha, Hbeta, MgII and/or CIV for 18100 quasars. Since LAMOST does not have absolute flux calibration information, we obtain the monochromatic continuum luminosities by fitting the SDSS photometric data using the quasar spectra, and then estimate the black hole masses. The catalog and spectra for these quasars are available online. This is the third installment in the series of LAMOST quasar survey which has released spectra for totally ~43,000 quasars hitherto. There are 24,772 independently discovered quasars, 17,128 of which are newly discovered. In addition to the great supplement to the new quasar discoveries, LAMOST has also provided a large database (overlapped with SDSS) for investigating the quasar spectral variability and discovering unusual quasars, including changing-look quasars, with ongoing and upcoming large surveys.

We present the final results from our survey of luminous $z \sim $ 5.5 quasars. This is the first systematic quasar survey focusing on quasars at $z \sim$ 5.5, during the post-reionization epoch. It has been challenging to select quasars at $5.3 < z < 5.7$ using conventional color selections, due to their similar optical colors to those of late-type stars, especially M dwarfs. We developed a new selection technique for $z \sim$ 5.5 quasars based on optical, near-IR, and mid-IR photometry, using data from the Sloan Digital Sky Survey (SDSS), PanSTARR1 (PS1), the UKIRT Infrared Deep Sky Surveys - Large Area Survey, the UKIRT Hemisphere Survey, the VISTA Hemisphere Survey, and the Wide Field Infrared Survey Explorer ($WISE$), covering $\sim$ 11000 deg$^2$ of high galactic latitude sky. In this paper, we present the discovery of 15 new quasars at $z\sim 5.5$. Together with results from Yang et al. (2017), our survey provides a complete, flux-limited sample of 31 quasars at $5.3 \le z \le 5.7$. We measure the quasar spatial density at $z \sim 5.5$ and $M_{1450} < -26.2$. Our result is consistent with the rapid decline of quasar spatial density from $z = 5$ to 6, with $k=-0.66 \pm 0.05$ ($\rho(z) \propto 10^{kz}$). In addition, we present a new survey using optical colors only from the full PS1 area for luminous quasars at $z = 5.0 - 5.5$, beyond the SDSS footprint, and report the preliminary results from this survey, including 51 new quasars discovered at $4.61\le z \le5.71$.

Strong gravitational lensing provides a powerful probe of the physical properties of quasars and their host galaxies. A high fraction of the most luminous high-redshift quasars was predicted to be lensed due to magnification bias. However, no multiple imaged quasar was found at z>5 in previous surveys. We report the discovery of J043947.08+163415.7, a strongly lensed quasar at z=6.51, the first such object detected at the epoch of reionization, and the brightest quasar yet known at z>5. High-resolution HST imaging reveals a multiple imaged system with a maximum image separation theta ~ 0.2", best explained by a model of three quasar images lensed by a low luminosity galaxy at z~0.7, with a magnification factor of ~50. The existence of this source suggests that a significant population of strongly lensed, high redshift quasars could have been missed by previous surveys, as standard color selection techniques would fail when the quasar color is contaminated by the lensing galaxy.

This is the third paper in a series aims at finding reionzation-era quasars with the combination of DESI Legacy imaging Surveys (DELS) and near-infrared imaging surveys, such as the UKIRT Hemisphere Survey (UHS), as well as the Wide-field Infrared Survey Explore ($WISE$) mid-infrared survey. In this paper, we describe the updated quasar candidate selection procedure, report the discovery of 16 quasars at $6.4\lesssim z \lesssim6.9$ from area of $\sim$13,020 deg$^2$, and present the quasar luminosity function (QLF) at $z\sim6.7$. The measured QLF follows $\Phi(L_{1450})\propto L_{1450}^{-2.35}$ in the magnitude range $27.6<M_{1450}<-25.5$. We determine the quasar comoving spatial density at $\langle z \rangle$=6.7 and $M_{1450}<-26.0$ to be $\rm 0.39\pm0.11 Gpc^{-3}$ and find that the exponential density evolution parameter to be $k=-0.78\pm0.18$ from $z\sim6$ to $z\sim6.7$, corresponding to a rapid decline by a factor of $\sim 6$ per unit redshift towards earlier epoch, a rate significantly faster than that at $z\sim 3- 5$. The cosmic time between $z\sim6$ and $z\sim6.7$ is only 121 Myrs. The quasar density declined by a factor of more than three within such short time requires that SMBHs must grow rapidly or they are less radiatively efficient at higher redshifts. We measured quasar comoving emissivity at $z\sim6.7$ which indicate that high redshift quasars are highly unlikely to make a significant contribution to hydrogen reionization. The broad absorption line (BAL) quasar fraction at $z\gtrsim6.5$ is measured to be $\gtrsim$22%. In addition, we also report the discovery of additional five quasars at $z\sim6$ in the appendix.

Galaxy clusters trace the largest structures of the Universe and provide ideal laboratories for studying galaxy evolution and cosmology. Clusters with extended X-ray emission have been discovered at redshifts up to z ~ 2.5. Meanwhile, there has been growing interest in hunting for protoclusters, the progenitors of clusters, at higher redshifts. It is, however, very challenging to find the largest protoclusters at early times when they start to assemble. Here we report a giant protocluster of galaxies at redshift z = 5.7, when the Universe was only one billion years old. This protocluster occupies a volume of about 35x35x35 cubic co-moving megaparsecs. It is embedded in an even larger overdense region with at least 41 spectroscopically confirmed, luminous Lyman-alpha emitting galaxies (Lyman-alpha Emitters, or LAEs), including several previously reported LAEs. Its LAE density is 6.6 times the average density at z ~ 5.7. It is the only one of its kind in an LAE survey in four square degrees on the sky. Such a large structure is also rarely seen in current cosmological simulations. This protocluster will collapse into a galaxy cluster with a mass of (3.6+/-0.9) x 10^{15} solar masses, comparable to those of the most massive clusters or protoclusters known to date.

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.

Collisional ring galaxies (CRGs) are formed through off-center collisions between a target galaxy and an intruder dwarf galaxy. We study the mass distribution and kinematics of the CRGs by tuning the bulge-to-disk mass ratio ($B/D$) for the progenitor; i.e., the target galaxy. We find that the lifetime of the ring correlates with the initial impact velocity vertical to the disk plane (i.e., $v_{z0}$). Three orbits for the collisional galaxy pair, on which clear and asymmetric rings form after collisions, are selected to perform the \textit{N}-body simulations at different values of $B/D$ for the progenitor. It is found that the ring structures are the strongest for the CRGs with small values of $B/D$. The S\'{e}rsic index, $n$, of the central remnant in the target galaxy becomes larger after collision. Moreover, the S\'{e}rsic index of a central remnant strongly correlates with the initial value of $B/D$ for the progenitor. A bulge-less progenitor results in a late-type object in the center of the ring galaxy, whereas a bulge-dominated progenitor leads to an early-type central remnant. Progenitors with $B/D\in [0.1,~0.3]$ (i.e., minor bulges) leave central remnants with $n\approx 4$. These results provide a possible explanation for the formation of a recently observed CRG with an early-type central nucleus, SDSS J1634+2049. In addition, we find that the radial and azimuthal velocity profiles for a ring galaxy are more sensitive to the $B/D$ than the initial relative velocity of the progenitor.

We present our deep $Y$-band imaging data of a two square degree field within the F22 region of the VIMOS VLT Deep Survey. The observations were conducted using the WIRCam instrument mounted at the Canada--France--Hawaii Telescope (CFHT). The total on-sky time was 9 hours, distributed uniformly over 18 tiles. The scientific goals of the project are to select faint quasar candidates at redshift $z>2.2$, and constrain the photometric redshifts for quasars and galaxies. In this paper, we present the observation and the image reduction, as well as the photometric redshifts that we derived by combining our $Y$-band data with the CFHTLenS $u^*g'r'i'z'$ optical data and UKIDSS DXS $JHK$ near-infrared data. With $J$-band image as reference total $\sim$80,000 galaxies are detected in the final mosaic down to $Y$-band $5\sigma$ point source limiting depth of 22.86 mag. Compared with the $\sim$3500 spectroscopic redshifts, our photometric redshifts for galaxies with $z<1.5$ and $i'\lesssim24.0$ mag have a small systematic offset of $|\Delta{z}|\lesssim0.2$, 1$\sigma$ scatter $0.03<\sigma_{\Delta z} < 0.06$, and less than 4.0% of catastrophic failures. We also compare to the CFHTLenS photometric redshifts, and find that ours are more reliable at $z\gtrsim0.6$ because of the inclusion of the near-infrared bands. In particular, including the $Y$-band data can improve the accuracy at $z\sim 1.0-2.0$ because the location of the 4000\AA-break is better constrained. The $Y$-band images, the multi-band photometry catalog and the photometric redshifts are released at \url{http://astro.pku.edu.cn/astro/data/DYI.html}.

The radiation mechanism of the prompt emission of gamma-ray bursts (GRBs) remains an open question. Although their spectra are usually well fitted with the empirical Band function, which is widely believed to be fully non-thermal and interpreted as an optically thin synchrotron emission, accumulating evidence shows that a thermal component actually exists. In this paper, a multi-color blackbody (mBB) model is proposed for the time-integrated spectrum of GRB 081221 by assuming a power-law distribution of the thermal luminosities with temperature, which manifests photospheric emissions from a different radius and/or angle. The effects of the minimum temperature $kT_{\min }$, the maximum temperature $kT_{\max }$ and the power law index $m$ of the luminosity distribution of an mBB are discussed. The fitting to the time-integrated spectrum during the bright phase (from 20s to 30s since the trigger) of GRB 081221 by the mBB model yields $kT_{\min }$ = $4.4\pm 0.3$ keV, $kT_{\max }$ = 57.0 $_{-1.4}^{+1.6}$ keV, and $m=-0.46_{-0.06}^{+0.05}$. When the time bin is small enough, the time-resolved spectra of GRB 081221 are well fitted with a series of single-temperature blackbodies. Our results imply the prompt emission of GRB 081221 is dominated by the photosphere emission and its time-integrated spectrum is a superposition of pure blackbody components at different times, indicating that some empirical Band spectra may be interpreted as mBB if the temperature is widely distributed.

The Beijing-Arizona Sky Survey (BASS) is a new wide-field legacy imaging survey in the northern Galactic cap using the 2.3m Bok telescope. The survey will cover about 5400 deg$^2$ in the $g$ and $r$ bands, and the expected 5$\sigma$ depths (corrected for the Galactic extinction) in the two bands are 24.0 and 23.4 mag, respectively. BASS started observations in January 2015, and has completed about 41% of the whole area as of July 2016. The first data release contains both calibrated images and photometric catalogs obtained in 2015 and 2016. The depths of single-epoch images in the two bands are 23.4 and 22.9 mag, and the full depths of three epochs are about 24.1 and 23.5 mag, respectively.

e-ASTROGAM (enhanced ASTROGAM) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV - the lower energy limit can be pushed to energies as low as 150 keV for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and LISA.

The eXTP (enhanced X-ray Timing and Polarimetry) mission is a major project of the Chinese Academy of Sciences (CAS) and China National Space Administration (CNSA) currently performing an extended phase A study and proposed for a launch by 2025 in a low-earth orbit. The eXTP scientific payload envisages a suite of instruments (Spectroscopy Focusing Array, Polarimetry Focusing Array, Large Area Detector and Wide Field Monitor) offering unprecedented simultaneous wide-band X-ray spectral, timing and polarimetry sensitivity. A large European consortium is contributing to the eXTP study and it is expected to provide key hardware elements, including a Large Area Detector (LAD). The LAD instrument for eXTP is based on the design originally proposed for the LOFT mission within the ESA context. The eXTP/LAD envisages a deployed 3.4 m2 effective area in the 2-30 keV energy range, achieved through the technology of the large-area Silicon Drift Detectors - offering a spectral resolution of up to 200 eV FWHM at 6 keV - and of capillary plate collimators - limiting the field of view to about 1 degree. In this paper we provide an overview of the LAD instrument design, including new elements with respect to the earlier LOFT configuration.

The eXTP (enhanced X-ray Timing and Polarimetry) mission is a major project of the Chinese Academy of Sciences (CAS) and China National Space Administration (CNSA) currently performing an extended phase A study and proposed for a launch by 2025 in a low-earth orbit. The eXTP scientific payload envisages a suite of instruments (Spectroscopy Focusing Array, Polarimetry Focusing Array, Large Area Detector and Wide Field Monitor) offering unprecedented simultaneous wide-band X-ray timing and polarimetry sensitivity. A large European consortium is contributing to the eXTP study and it is expected to provide key hardware elements, including a Wide Field Monitor (WFM). The WFM instrument for eXTP is based on the design originally proposed for the LOFT mission within the ESA context. The eXTP/WFM envisages a wide field X-ray monitor system in the 2-50 keV energy range, achieved through the technology of the large-area Silicon Drift Detectors. The WFM will consist of 3 pairs of coded mask cameras with a total combined Field of View (FoV) of 90x180 degrees at zero response and a source localization accuracy of ~1 arcmin. In this paper we provide an overview of the WFM instrument design, including new elements with respect to the earlier LOFT configuration, and anticipated performance.

e-ASTROGAM is a concept for a breakthrough observatory space mission carrying a gamma-ray telescope dedicated to the study of the non-thermal Universe in the photon energy range from 0.15 MeV to 3 GeV. The lower energy limit can be pushed down to energies as low as 30 keV for gamma-ray burst detection with the calorimeter. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with remarkable polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous and current generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will be a major player of the multiwavelength, multimessenger time-domain astronomy of the 2030s, and provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LISA, LIGO, Virgo, KAGRA, the Einstein Telescope and the Cosmic Explorer, IceCube-Gen2 and KM3NeT, SKA, ALMA, JWST, E-ELT, LSST, Athena, and the Cherenkov Telescope Array.

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.

Well-sampled optical light curves of 50 gamma-ray bursts (GRBs) with plateau features are compiled from the literature. By empirical fitting, we obtained the parameters of the optical plateaus, such as the decay slopes ($\alpha_{\rm 1}$ and $\alpha_{\rm 2}$), the break times ($T_{\rm b}$), and the corresponding optical fluxes ($F_{\rm b}$) at the break times. The break time of optical plateaus ranges from tens of seconds to $10^6$ seconds, with a typical value about $10^4$ seconds. We have calculated the break luminosity, and it mainly ranges from $10^{44}$ erg $s^{-1}$ to $10^{47}$ erg $s^{-1}$, which is generally two or three orders of magnitude less than the corresponding break luminosity of the X-ray afterglow plateaus. We reanalyzed the optical plateaus and also found that a significantly tighter correlation exists when we added the isotropic equivalent energy of GRBs $E_{\rm \gamma,iso}$ into the $L_{\rm b,z}-T_{\rm b,z}$ relation. The best fit correlation is obtained to be $L_{\rm b,z}\propto T_{\rm b,z}^{-0.9}E_{\rm \gamma,iso}^{0.4}$. We next explored the possible correlations among $L_{\rm b,z}$, $T_{\rm b,z}$ and $E_{\rm p,i}$, and found there is also a tight correlation between them, which takes the form of $L_{\rm b,z}\propto T_{\rm b,z}^{-0.9}E_{\rm p,i}^{0.5}$. We argue that these two tight $L_{\rm b,z}-T_{\rm b,z}-E_{\rm \gamma,iso}$ and $L_{\rm b,z}-T_{\rm b,z}-E_{\rm p,i}$ correlations are more physical, and it may be directly related to radiation physics of GRBs. The tight correlations are possible to be used as standard candles.

DAMPE (DArk Matter Particle Explorer) is a spaceborne high-energy cosmic ray and gamma-ray detector, successfully launched in December 2015. It is designed to probe astroparticle physics in the broad energy range from few GeV to 100 TeV. The scientific goals of DAMPE include the identification of possible signatures of Dark Matter annihilation or decay, the study of the origin and propagation mechanisms of cosmic-ray particles, and gamma-ray astronomy. DAMPE consists of four sub-detectors: a plastic scintillator strip detector, a Silicon-Tungsten tracKer-converter (STK), a BGO calorimeter and a neutron detector. The STK is composed of six double layers of single-sided silicon micro-strip detectors interleaved with three layers of tungsten for photon conversions into electron-positron pairs. The STK is a crucial component of DAMPE, allowing to determine the direction of incoming photons, to reconstruct tracks of cosmic rays and to estimate their absolute charge (Z). We present the in-flight performance of the STK based on two years of in-flight DAMPE data, which includes the noise behavior, signal response, thermal and mechanical stability, alignment and position resolution.

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.

X-ray flares, lasting for $\sim 100 - 1000$ s in the X-ray band, are often observed following gamma-ray bursts (GRBs). The physical origin of X-ray flares is still unknown merely with the temporal/spectral information. On the other hand, some polarimeters are expected to be launched within several years thanks to the increasing interest on astronomical X-ray polarimetry. Here, by assuming that X-ray flares are synchrotron radiation from relativistic spherical shells, we show that the linear polarization degree during the rising phase of an X-ray flare is much higher for the emitting region with toroidal magnetic fields than that with random magnetic fields. In the decay phase of the flare, the evolution of the polarization degree is determined by the curvature effect of the emitting shell, which is a natural feature of jet scenarios for flares. Therefore, the measurement of the polarization of X-ray flares would provide a useful tool to probe the configuration of magnetic fields in the emission region, and may even help to test the curvature effect. The information on the magnetic configuration can further help us to understand the properties of GRB jets.

We have observed the remnant of supernova SN~1987A (SNR~1987A), located in the Large Magellanic Cloud (LMC), to search for periodic and/or transient radio emission with the Parkes 64\,m-diameter radio telescope. We found no evidence of a radio pulsar in our periodicity search and derived 8$\sigma$ upper bounds on the flux density of any such source of $31\,\mu$Jy at 1.4~GHz and $21\,\mu$Jy at 3~GHz. Four candidate transient events were detected with greater than $7\sigma$ significance, with dispersion measures (DMs) in the range 150 to 840\,cm$^{-3}\,$pc. For two of them, we found a second pulse at slightly lower significance. However, we cannot at present conclude that any of these are associated with a pulsar in SNR~1987A. As a check on the system, we also observed PSR~B0540$-$69, a young pulsar which also lies in the LMC. We found eight giant pulses at the DM of this pulsar. We discuss the implications of these results for models of the supernova remnant, neutron star formation and pulsar evolution.

POLAR is a compact wide-field space-borne detector for precise measurements of the linear polarisation of hard X-rays emitted by transient sources in the energy range from 50 keV to 500 keV. It consists of a 40$\times$40 array of plastic scintillator bars used as a detection material. The bars are grouped in 25 detector modules. The energy range sensitivity of POLAR is optimized to match with the prompt emission photons from the gamma-ray bursts (GRBs). Polarization measurements of the prompt emission would probe source geometries, emission mechanisms and magnetic structures in GRB jets. The instrument can also detect hard X-rays from solar flares and be used for precise measurement of their polarisation. POLAR was launched into a low Earth orbit on-board the Chinese space-lab TG-2 on September 15th, 2016. To achieve high accuracies in polarisation measurements it is essential to assure both before and during the flight a precise energy calibration. Such calibrations are performed with four low activity $^{22}$Na radioactive sources placed inside the instrument. Energy conversion factors are related to Compton edge positions from the collinear annihilation photons from the sources. This paper presents main principles of the in-flight calibration, describes studies of the method based on Monte Carlo simulations and its laboratory verification and finally provides some observation results based on the in-flight data analysis.

The rare case of changing-look (CL) AGNs, with the appearance or disappearance of broad Balmer emission lines within a few years, challenges our understanding of the AGN unified model. We present a sample of 21 new CL AGNs at $0.08<z<0.58$, which doubles the number of such objects known to date. These new CL AGNs were discovered by several ways, from (1) repeat spectra in the SDSS, (2) repeat spectra in the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) and SDSS, and (3) photometric variability and new spectroscopic observations. We use the photometric data from surveys, including the SDSS imaging survey, the Pan-STARRS1, the DESI Legacy imaging survey, the Wide-field Infrared Survey Explorer (WISE), the Catalina Real-time Transient Survey, and the Palomar Transient Factory. The estimated upper limits of transition timescale of the CL AGNs in this sample spans from 0.9 to 13 years in the rest frame. The continuum flux in the optical and mid-infrared becomes brighter when the CL AGNs turn on, or vice versa. Variations of more than 0.2 mag in $W1$ band were detected in 15 CL AGNs during the transition. The optical and mid-infrared variability is not consistent with the scenario of variable obscuration in 10 CL AGNs at more than $3\sigma$ confidence level. We confirm a bluer-when-brighter trend in the optical. However, the mid-infrared WISE colors $W1-W2$ become redder when the objects become brighter in the $W1$ band, possibly due to a stronger hot dust contribution in the $W2$ band when the AGN activity becomes stronger. The physical mechanism of type transition is important for understanding the evolution of AGNs.

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.

POLAR is a compact space-borne detector designed to perform reliable measurements of the polarization for transient sources like Gamma-Ray Bursts in the energy range 50-500keV. The instrument works based on the Compton Scattering principle with the plastic scintillators as the main detection material along with the multi-anode photomultiplier tube. POLAR has been launched successfully onboard the Chinese space laboratory TG-2 on 15th September, 2016. In order to reliably reconstruct the polarization information a highly detailed understanding of the instrument is required for both data analysis and Monte Carlo studies. For this purpose a full study of the in-orbit performance was performed in order to obtain the instrument calibration parameters such as noise, pedestal, gain nonlinearity of the electronics, threshold, crosstalk and gain, as well as the effect of temperature on the above parameters. Furthermore the relationship between gain and high voltage of the multi-anode photomultiplier tube has been studied and the errors on all measurement values are presented. Finally the typical systematic error on polarization measurements of Gamma-Ray Bursts due to the measurement error of the calibration parameters are estimated using Monte Carlo simulations.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

We present the ensemble variability analysis results of quasars using the Dark Energy Camera Legacy Survey (DECaLS) and the Sloan Digital Sky Survey (SDSS) quasar catalogs. Our dataset includes 119,305 quasars with redshifts up to 4.89. Combining the two datasets provides a 15-year baseline and permits analysis of the long timescale variability. Adopting a power-law form for the variability structure function, $V=A(t/1yr)^{\gamma}$, we use the multi-dimensional parametric fitting to explore the relationships between the quasar variability amplitude and a wide variety of quasar properties, including redshift (positive), bolometric luminosity (negative), rest-frame wavelength (negative), and black hole mass (uncertain). We also find that $\gamma$ can be also expressed as a function of redshift (negative), bolometric luminosity (positive), rest-frame wavelength (positive), and black hole mass (positive). Tests of the fitting significance with the bootstrap method show that, even with such a large quasar sample, some correlations are marginally significant. The typical value of $\gamma$ for the entire dataset is $\gtrsim 0.25$, consistent with the results in previous studies on both the quasar ensemble variability and the structure function. A significantly negative correlation between the variability amplitude and the Eddington ratio is found, which may be explained as an effect of accretion disk instability.

Gamma-ray bursts (GRBs) have been widely argued to originate from binary compact object mergers or core collapses of massive stars. Jets from these systems may have two components, an inner, narrow sub-jet and an outer, wider sub-jet. Such a jet subsequently interacts with its ambient gas, leading to a reverse shock (RS) and a forward shock (FS). The magnetic field in the narrow sub-jet is very likely to be mixed by an ordered component and a random component during the afterglow phase. In this paper, we calculate light curves and polarization evolution of optical afterglows with this mixed magnetic field in the RS region of the narrow sub-jet in a two-component jet model. The resultant light curve has two peaks: an early peak arises from the narrow sub-jet and a late-time rebrightening is due to the wider sub-jet. We find the polarization degree (PD) evolution under such a mixed magnetic field confined in the shock plane is very similar to that under the purely ordered magnetic field condition. The two-dimensional `mixed' magnetic fields confined in the shock plane are essentially the ordered magnetic fields only with different configurations. The position angle (PA) of the two-component jet can change gradually or abruptly by $90^\circ$. In particular, an abrupt $90^\circ$ change of the PA occurs when the PD changes from its decline phase to rise phase.

In this paper, we explore the physics of the accretion and jet in narrow-line Seyfert 1 galaxies (NLS1). Specifically, we compile a sample composed of 16 nearby NLS1 with $L_{\rm bol}/L_{\rm Edd} \gtrsim 0.1$. We investigate the mutual correlation between their radio luminosity $L_{\rm R}$, X-ray luminosity $L_{\rm X}$, optical luminosity $L_{\rm 5100}$ and black hole mass $M_{\rm BH}$. By adopting partial correlation analysis: (1) we find a positive correlation between $L_{\rm X}$ and $M_{\rm BH}$, and (2) we find a weak positive correlation between $L_{\rm R}$ and $L_{5100}$. However, we don't find significant correlations between $L_{\rm R}$ and $L_{\rm X}$ or between $L_{\rm X}$ and $L_{5100}$ after considering the effect of the black hole mass, which leads to a finding of the independence of $L_{\rm X}/L_{\rm Edd}$ on $L_{5100}/L_{\rm Edd}$. Interestingly, the findings that $L_{\rm X}$ is correlated with $M_{\rm BH}$ and $L_{\rm X}/L_{\rm Edd}$ is not correlated with $L_{5100}/L_{\rm Edd}$ support that the X-ray emission is saturated with increasing $\dot{M}$ for $L_{\rm bol}/L_{\rm Edd} \gtrsim 0.1$ in NLS1s, which may be understood in the framework of slim disc scenario. Finally, we suggest that a larger NLS1 sample with high quality radio and X-ray data is needed to further confirm this result in the future.

Recent follow-up observations of the binary neutron star (NS) merging event GW170817/SGRB 170817A reveal that its X-ray/optical/radio emissions are brightening continuously up to $\sim 100$ days post-merger. This late-time brightening is unexpected from the kilonova model or the off-axis top-hat jet model for gamma-ray burst afterglows. In this paper, by assuming that the merger remnant is a long-lived NS, we propose that the interaction between an electron-positron-pair ($e^+e^-$) wind from the central NS and the jet could produce a long-lived reverse shock, from which a new emission component would rise and can interpret current observations well. The magnetic-field-induced ellipticity of the NS is taken to be $4 \times 10^{-5}$ in our modeling, so that the braking of the NS is mainly through the gravitational wave (GW) radiation rather than the magnetic dipole radiation, and the emission luminosity at early times would not exceed the observational limits. In our scenario, since the peak time of the brightening is roughly equal to the spin-down time scale of the NS, the accurate peak time may help constrain the ellipticity of the remnant NS. We suggest that radio polarization observations of the brightening would help to distinguish our scenario from other scenarios. Future observations on a large sample of short gamma-ray burst afterglows or detections of GW signals from merger remnants would test our scenario.

This is the second installment for the Large Sky Area Multi-Object Fibre Spectroscopic Telescope (LAMOST) Quasar Survey, which includes quasars observed from September 2013 to June 2015. There are 9024 confirmed quasars in DR2 and 10911 in DR3. After cross-match with the SDSS quasar catalogs and NED, 12126 quasars are discovered independently. Among them 2225 quasars were released by SDSS DR12 QSO catalogue in 2014 after we finalised the survey candidates. 1801 sources were identified by SDSS DR14 as QSOs. The remaining 8100 quasars are considered as newly founded, and among them 6887 quasars can be given reliable emission line measurements and the estimated black hole masses. Quasars found in LAMOST are mostly located at low-to-moderate redshifts, with a mean value of 1.5. The highest redshift observed in DR2 and DR3 is 5. We applied emission line measurements to H$\alpha$, H$\beta$, Mg{\sc ii} and C{\sc iv}. We deduced the monochromatic continuum luminosities using photometry data, and estimated the virial black hole masses for the newly discovered quasars. Results are compiled into a quasar catalog, which will be available online.

We have investigated the mode-changing properties of PSR B0329+54 using 31 epochs of simultaneous 13 cm/3 cm single-pulse observations obtained with Shanghai Tian Ma 65 m telescope. The pulsar was found in the abnormal emission mode 17 times, accounting for ~13% of the 41.6 hours total observation time. Single pulse analyses indicate that mode changes took place simultaneously at 13 cm/3 cm within a few rotational periods. We detected occasional bright and narrow pulses whose peak flux densities were 10 times higher than that of the integrated profile in both bands. At 3 cm, about 0.66% and 0.27% of single pulses were bright in the normal mode and abnormal mode respectively, but at 13 cm the occurrence rate was only about 0.007%. We divided the pulsar radiation window into three components (C1, C2 and C3) corresponding to the main peaks of the integrated profile. The bright pulses preferentially occurred at pulse phases corresponding to the peaks of C2 and C3. Fluctuation spectra showed that C2 had excess red noise in the normal mode, but broad quasi-periodic features with central frequencies around 0.12 cycles/period in the abnormal mode. At 3 cm, C3 had a stronger quasi-periodic modulation centered around 0.06 cycles/period in the abnormal mode. Although there were some asymmetries in the two-dimensional fluctuation spectra, we found no clear evidence for systematic subpulse drifting. Consistent with previous low-frequency observations, we found a very low nulling probability for B0329+54 with upper limits of 0.13% and 1.68% at 13 cm/3 cm respectively.

Limited by the sensitivities of the current gravitational wave (GW) detectors, the central remnant of the binary neutron star (NS) merger associated with GW170817 remains an open question. Considering the relatively large total mass, it is generally proposed that the merger of GW170817 would lead to a shortly lived hypermassive NS or directly produce a black hole (BH). There is no clear evidence to support or rule out a long-lived NS as the merger remnant. Here we utilize the GW and electromagnetic (EM) signals to comprehensively investigate the parameter space that allows a long-lived NS to survive as the merger remnant of GW170817. We find that for some stiff equations of state, the merger of GW170817 could, in principle, lead to a massive NS, which has a millisecond spin period. The post-merger GW signal could hardly constrain the ellipticity of the NS. If the ellipticity reaches $10^{-3}$, in order to be compatible with the multi-band EM observations, the dipole magnetic field of the NS ($B_p$) is constrained to the magnetar level of $\sim10^{14}$ G. If the ellipticity is smaller than $10^{-4}$, $B_p$ is constrained to the level of $\sim10^{10}-10^{12}\,$G. These conclusions weakly depend on the adoption of equations of state.

We present the discovery of strong Balmer line absorption in H$\alpha$ to H$\gamma$ in two luminous low-ionization broad absorption line quasars (LoBAL QSOs) at z~1.5, with black hole masses around $10^{10}$ $M_\odot$ from near-IR spectroscopy. There are only two previously known quasars at z>1 showing Balmer line absorption. SDSS J1019+0225 shows blueshifted absorption by ~1400 km/s with an H$\alpha$ rest-frame equivalent width of 13 \AA{}. In SDSS J0859+4239 we find redshifted absorption by ~500 km/s with an H$\alpha$ rest-frame equivalent width of 7 \AA{}. The redshifted absorption could indicate an inflow of high density gas onto the black hole, though we cannot rule out alternative interpretations. The Balmer line absorption in both objects appears to be saturated, indicating partial coverage of the background source by the absorber. We estimate the covering fractions and optical depth of the absorber and derive neutral hydrogen column densities, $N_{\rm{H\,I}}\sim1.3\times 10^{18}$ cm$^{-2}$ for SDSS J1019+0225 and $N_{\rm{H\,I}}\sim9\times 10^{17}$ cm$^{-2}$ for SDSS J0859+4239, respectively. In addition, the optical spectra reveal also absorption troughs in HeI* $\lambda3889$ and $\lambda3189$ in both objects.

Double neutron star (DNS) merger events are promosing candidates of short Gamma-ray Burst (sGRB) progenitors as well as high-frequecy gravitational wave (GW) emitters. On August 17, 2017, such a coinciding event was detected by both the LIGO-Virgo gravitational wave detector network as GW170817 and Gamma-Ray Monitor on board NASA's {\it Fermi} Space Telescope as GRB 170817A. Here we show that the fluence and spectral peak energy of this sGRB fall into the lower portion of the distributions of known sGRBs. Its peak isotropic luminosity is abnormally low. The estimated event rate density above this luminosity is at least $190^{+440}_{-160} {\rm Gpc^{-3} \ yr^{-1}}$, which is close to but still below the DNS merger event rate density. This event likely originates from a structured jet viewed from a large viewing angle. There are similar faint soft GRBs in the {\it Fermi} archival data, a small fraction of which might belong to this new population of nearby, low-luminosity sGRBs.

We report the result of a faint quasar survey in a one square degree field. The aim is to test the Y-K/g-z and J-K/i-Y color selection criteria for quasars at faint magnitude, to obtain a complete sample of quasars based on deep optical and near-infrared color-color selection, and to measure the faint end of quasar luminosity function (QLF) over a wide redshift range. We carried out a quasar survey based on the Y-K/g-z and J-K/i-Y quasar selection criteria, using the deep Y-band data obtained from our CFHT/WIRCam Y-band images in a two-degree field within the F22 field of the VIMOS VLT deep survey, optical co-added data from Sloan Digital Sky Survey Stripe 82 and deep near-infrared data from the UKIDSS Deep Extragalactic Survey in the same field. We discovered 25 new quasars at 0.5 < z < 4.5 and i < 22.5 mag within one square degree field. The survey significantly increases the number of faint quasars in this field, especially at z ~ 2-3. It confirms that our color selections are highly complete in a wide redshift range (z < 4.5), especially over the quasar number density peak at z ~ 2-3, even for faint quasars. Combining all previous known quasars and new discoveries, we construct a sample with 109 quasars, and measure the binned QLF and parametric QLF. Although the sample is small, our results agree with a pure luminosity evolution at lower redshift and luminosity evolution and density evolution model at redshift z > 2.5.

We have analyzed multi-passband photometric observations, obtained with the {\it Hubble Space Telescope}, of the massive ($1.8 \times 10^5 M_\odot$), intermediate-age (1.8 Gyr-old) Large Magellanic Cloud star cluster NGC 1783. The morphology of the cluster's red giant branch does not exhibit a clear broadening beyond its intrinsic width; the observed width is consistent with that owing to photometric uncertainties alone and independent of our photometric selection boundaries applied to obtain our sample of red-giant stars. The color dispersion of the cluster's red-giant stars around the best-fitting ridgeline is $0.062 \pm 0.009$ mag, which is equivalent to the width of $0.080 \pm 0.001$ mag derived from artificial simple stellar population tests, that is, tests based on single-age, single-metallicity stellar populations. NGC 1783 is comparably massive as other star clusters that show clear evidence of multiple stellar populations. After incorporating mass-loss recipes from its current age of 1.8 Gyr to an age of 6 Gyr, NGC 1783 is expected to remain as massive as some other clusters that host clear multiple populations at these intermediate ages. If we were to assume that mass is an important driver of multiple population formation, then NGC 1783 should have exhibited clear evidence of chemical abundance variations. However, our results support the absence of any chemical abundance variations in NGC 1783.

We modified the broadband photometric reverberation mapping (PRM) code, JAVELIN, and tested the availability to get broad line region (BLR) time delays that are consistent with the spectroscopic reverberation mapping (SRM) project SDSS-RM. The broadband light curves of SDSS-RM quasars produced by convolution with the system trans- mission curves were used in the test. We found that under similar sampling conditions (evenly and frequently sampled), the key factor determining whether the broadband PRM code can yield lags consistent with the SRM project is the flux ratio of the broad emission line to the reference continuum, which is in line with the previous findings. We further found a critical line-to-continuum flux ratio, about 6%, above which the mean of the ratios between the lags from PRM and SRM becomes closer to unity, and the scatter is pronouncedly reduced. We also tested our code on a subset of SDSS Stripe 82 quasars, and found that our program tends to give biased lag estimations due to the observation gaps when the R-L relation prior in Markov Chain Monte Carlo (MCMC) is discarded. The performance of the damped random walk (DRW) model and the power-law (PL) structure function model on broadband PRM were compared. We found that given both SDSS-RM-like or Stripe 82-like light curves, the DRW model performs better in carrying out broadband PRM than the PL model.

Radio-loud active galactic nuclei (AGNs), hosting powerful relativistic jet outflows, provide an excellent laboratory for studying jet physics. Very long baseline interferometry (VLBI) enables high-resolution imaging on milli-arcsecond (mas) and sub-mas scales, making it a powerful tool to explore the inner jet structure, shedding light on the formation, acceleration and collimation of AGN jets. In this paper, we present Very Long Baseline Array (VLBA) observations of ten radio-loud AGNs at 43 and 86~GHz, which were selected from the {\it Planck} catalogue of compact sources and are among the brightest in published VLBI images at and below 15 GHz. The image noise levels in our observations are typically 0.3 mJy beam$^{-1}$ and 1.5 mJy beam$^{-1}$ at 43 and 86 GHz, respectively. Compared with the VLBI data observed at lower frequencies from the literature, our observations with higher resolution (the highest resolution up to 0.07 mas at 86 GHz and 0.18 mas at 43 GHz) and at higher frequencies detected new jet components at sub-parsec scales, offering valuable data for studies of the physical properties of innermost jets. These include compactness factor of the radio structure (the ratio of core flux density to total flux density), and core brightness temperature ($T_{\rm b}$). In all these sources, the compact core accounts for a significant fraction ($> 60\%$) of the total flux density. Their correlated flux density at the longest baselines is higher than 0.16 Jy. The compactness of these sources make them good phase calibrators of mm-wavelength ground-based and space VLBI.

We study the evolution of star clusters located in the outer regions of a galaxy undergoing a sudden mass loss through gas expulsion in the framework of Milgromian dynamics (MOND) by means of N-body simulations. We find that, to leave a bound star cluster, the star formation efficiency (SFE) of an embedded cluster dominated by deep MOND gravity can be reduced down to $2.5\%$. For a given SFE, the star clusters that survive in MOND can bind a larger fraction of mass compared to the Newtonian dynamics. Moreover, the more diffuse the embedded cluster is, the less substantial the size expansion of the final star cluster is. The density profiles of a surviving star cluster are more cuspy in the centre for more massive embedded clusters, and the central density profiles are flatter for less massive embedded clusters or for lower SFE. This work may help to understand the low concentration and extension of the distant low-density globular clusters (GCs) and ultra-faint and diffuse satellite galaxies around the Milky Way.

Very recently, the gravitational-wave (GW) event GW170817 was discovered to be associated with the short gamma-ray burst (GRB) 170817A. Multi-wavelength follow-up observations were carried out, and X-ray, optical, and radio counterparts to GW170817 were detected. The observations undoubtedly indicate that GRB 170817A originates from a binary neutron star merger. However, the GRB falls into the low-luminosity class that could have a higher statistical occurrence rate and detection probability than the normal (high-luminosity) class. This implies the possibility that GRB 170817A is intrinsically powerful, but we are off-axis and only observe its side emission. In this Letter, we provide a timely modeling of the multi-wavelength afterglow emission from this GRB and the associated kilonova signal from the merger ejecta, under the assumption of a structured jet, a two-component jet, and an intrinsically less-energetic quasi-isotropic fireball, respectively. Comparing the afterglow properties with the multi-wavelength follow-up observations, we can distinguish between these three models. Furthermore, a few model parameters (e.g., the ejecta mass and velocity) can be constrained.

The low-energy spectra of gamma-ray bursts' (GRBs) prompt emission are closely related to the energy distribution of electrons, which is further regulated by their cooling processes. We develop a numerical code to calculate the evolution of the electron distribution with given initial parameters, in which three cooling processes (i.e., adiabatic, synchrotron and inverse Compton cooling) and the effect of decaying magnetic field are coherently considered. A sequence of results are presented by exploring the plausible parameter space for both the fireball and the Poynting-flux-dominated regime. Different cooling patterns for the electrons can be identified and they are featured by a specific dominant cooling mechanism. Our results show that the hardening of the low-energy spectra can be attributed to the dominance of synchrotron self-Compton cooling within the internal shock model, or to decaying synchrotron cooling within the Poynting-flux-dominated jet scenario. These two mechanisms can be distinguished by observing the hard low-energy spectra of isolated short pulses in some GRBs. The dominance of adiabatic cooling can also lead to hard low-energy spectra when the ejecta is moving at an extreme relativistic speed. The information from the time-resolved low-energy spectra can help to probe the physical characteristics of the GRB ejecta via our numerical results.

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.

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 present near-IR spectroscopy of 22 luminous low-ionization broad absorption line quasars (LoBAL QSOs) at redshift 1.3<z<2.5, with 12 objects at z~1.5 and 10 at z~2.3. The spectra cover the rest-frame H$\alpha$ and H$\beta$ line regions, allowing us to obtain robust black hole mass estimates based on the broad H$\alpha$ line. We use these data, augmented by a lower-redshift sample from the Sloan Digital Sky Survey, to test the proposed youth scenario for LoBALs, which suggests that LoBALs constitute an early short-lived evolutionary stage of quasar activity, by probing for any difference in their masses, Eddington ratios, or rest-frame optical spectroscopic properties compared to normal quasars. In addition, we construct the UV to mid-IR spectral energy distributions (SEDs) for the LoBAL sample and a matched non-BAL quasar sample. We do not find any statistically significant difference between LoBAL QSOs and non-BAL QSOs in their black hole mass or Eddington ratio distributions. The mean UV to mid-IR SED of the LoBAL QSOs is consistent with non-BAL QSOs, apart from their stronger reddening. At z>1 there is no clear difference in their optical emission line properties. We do not see particularly weak [OIII] or strong FeII emission. The LoBAL QSOs do not show a stronger prevalence of ionized gas outflows as traced by the [OIII] line, compared to normal QSOs of similar luminosity. We conclude that the optical-MIR properties of LoBAL QSOs are consistent with the general quasar population and do not support them to constitute a special phase of active galactic nucleus evolution.

We present a new algorithm to estimate quasar photometric redshifts (photo-$z$s), by considering the asymmetries in the relative flux distributions of quasars. The relative flux models are built with multivariate Skew-t distributions in the multi-dimensional space of relative fluxes as a function of redshift and magnitude. For 151,392 quasars in the SDSS, we achieve a photo-$z$ accuracy, defined as the fraction of quasars with the difference between the photo-$z$ $z_p$ and the spectroscopic redshift $z_s$, $|\Delta z| = |z_s-z_p|/(1+z_s)$ within 0.1, of 74%. Combining the WISE W1 and W2 infrared data with the SDSS data, the photo-$z$ accuracy is enhanced to 87%. Using the Pan-STARRS1 or DECaLS photometry with WISE W1 and W2 data, the photo-$z$ accuracies are 79% and 72%, respectively. The prior probabilities as a function of magnitude for quasars, stars and galaxies are calculated respectively based on (1) the quasar luminosity function; (2) the Milky Way synthetic simulation with the Besan\c{c}on model; (3) the Bayesian Galaxy Photometric Redshift estimation. The relative fluxes of stars are obtained with the Padova isochrones, and the relative fluxes of galaxies are modeled through galaxy templates. We test our classification method to select quasars using the DECaLS $g$, $r$, $z$, and WISE W1 and W2 photometry. The quasar selection completeness is higher than 70% for a wide redshift range $0.5<z<4.5$, and a wide magnitude range $18<r<21.5$ mag. Our photo-$z$ regression and classification method has the potential to extend to future surveys. The photo-$z$ code will be publicly available.

The LIGO detection of gravitational waves (GW) from merging black holes in 2015 marked the beginning of a new era in observational astronomy. The detection of an electromagnetic signal from a GW source is the critical next step to explore in detail the physics involved. The Antarctic Survey Telescopes (AST3), located at Dome A, Antarctica, is uniquely situated for rapid response time-domain astronomy with its continuous night-time coverage during the austral winter. We report optical observations of the GW source (GW~170817) in the nearby galaxy NGC 4993 using AST3. The data show a rapidly fading transient at around 1 day after the GW trigger, with the $i$-band magnitude declining from $17.23\pm0.13$ magnitude to $17.72\pm0.09$ magnitude in $\sim 1.8$ hour. The brightness and time evolution of the optical transient associated with GW~170817 are broadly consistent with the predictions of models involving merging binary neutron stars. We infer from our data that the merging process ejected about $\sim 10^{-2}$ solar mass of radioactive material at a speed of up to $30\%$ the speed of light.

The discovery of the first electromagnetic counterpart to a gravitational wave signal has generated follow-up observations by over 50 facilities world-wide, ushering in the new era of multi-messenger astronomy. In this paper, we present follow-up observations of the gravitational wave event GW170817 and its electromagnetic counterpart SSS17a/DLT17ck (IAU label AT2017gfo) by 14 Australian telescopes and partner observatories as part of Australian-based and Australian-led research programs. We report early- to late-time multi-wavelength observations, including optical imaging and spectroscopy, mid-infrared imaging, radio imaging, and searches for fast radio bursts. Our optical spectra reveal that the transient source afterglow cooled from approximately 6400K to 2100K over a 7-day period and produced no significant optical emission lines. The spectral profiles, cooling rate, and photometric light curves are consistent with the expected outburst and subsequent processes of a binary neutron star merger. Star formation in the host galaxy probably ceased at least a Gyr ago, although there is evidence for a galaxy merger. Binary pulsars with short (100 Myr) decay times are therefore unlikely progenitors, but pulsars like PSR B1534+12 with its 2.7 Gyr coalescence time could produce such a merger. The displacement (about 2.2 kpc) of the binary star system from the centre of the main galaxy is not unusual for stars in the host galaxy or stars originating in the merging galaxy, and therefore any constraints on the kick velocity imparted to the progenitor are poor.

We present a spectroscopic redshift catalog from the LAMOST Complete Spectroscopic Survey of Pointing Area (LaCoSSPAr) in the Southern Galactic Cap (SGC), which is designed to observe all sources (Galactic and extra-galactic) by using repeating observations with a limiting magnitude of $r=18.1~mag$ in two $20~deg^2$ fields. The project is mainly focusing on the completeness of LAMOST ExtraGAlactic Surveys (LEGAS) in the SGC, the deficiencies of source selection methods and the basic performance parameters of LAMOST telescope. In both fields, more than 95% of galaxies have been observed. A post-processing has been applied to LAMOST 1D spectrum to remove the majority of remaining sky background residuals. More than 10,000 spectra have been visually inspected to measure the redshift by using combinations of different emission/absorption features with uncertainty of $\sigma_{z}/(1+z)<0.001$. In total, there are 1528 redshifts (623 absorption and 905 emission line galaxies) in Field A and 1570 redshifts (569 absorption and 1001 emission line galaxies) in Field B have been measured. The results show that it is possible to derive redshift from low SNR galaxies with our post-processing and visual inspection. Our analysis also indicates that up to 1/4 of the input targets for a typical extra-galactic spectroscopic survey might be unreliable. The multi-wavelength data analysis shows that the majority of mid-infrared-detected absorption (91.3%) and emission line galaxies (93.3%) can be well separated by an empirical criterion of $W2-W3=2.4$. Meanwhile, a fainter sequence paralleled to the main population of galaxies has been witnessed both in $M_r$/$W2-W3$ and $M_*$/$W2-W3$ diagrams, which could be the population of luminous dwarf galaxies but contaminated by the edge-on/highly inclined galaxies ($\sim30\%$).

The ejecta composition is an open question in gamma-ray bursts (GRB) physics. Some GRBs possess a quasi-thermal spectral component in the time-resolved spectral analysis, suggesting a hot fireball origin. Others show a featureless non-thermal spectrum known as the "Band" function, consistent with a synchrotron radiation origin and suggesting that the jet is Poynting-flux-dominated at the central engine and likely in the emission region as well. There are also bursts showing a sub-dominant thermal component and a dominant synchrotron component suggesting a likely hybrid jet composition. Here we report an extraordinarily bright GRB 160625B, simultaneously observed in gamma-rays and optical wavelengths, whose prompt emission consists of three isolated episodes separated by long quiescent intervals, with the durations of each "sub-burst" being $\sim$ 0.8 s, 35 s, and 212 s, respectively. Its high brightness (with isotropic peak luminosity L$_{\rm p, iso}\sim 4\times 10^{53}$ erg/s) allows us to conduct detailed time-resolved spectral analysis in each episode, from precursor to main burst and to extended emission. The spectral properties of the first two sub-bursts are distinctly different, allowing us to observe the transition from thermal to non-thermal radiation between well-separated emission episodes within a single GRB. Such a transition is a clear indication of the change of jet composition from a fireball to a Poynting-flux-dominated jet.

A binary neutron star (BNS) merger has been widely argued to be one of the progenitors of a short gamma-ray burst (SGRB). This central engine can be verified if its gravitational-wave (GW) event is detected simultaneously. Once confirmed, this kind of association will be a landmark in multi-messenger astronomy and will greatly enhance our understanding of the BNS merger processes. Due to the limited detection horizon of BNS mergers for the advanced LIGO/Virgo GW observatories, we are inclined to local SGRBs within few hundreds of mega-parsecs. Since normal SGRBs rarely fall into such a close range, to make it more observationally valuable, we have to focus on low-luminosity SGRBs which have a higher statistical occurrence rate and detection probability. However, there is a possibility that an observed low-luminosity SGRB is intrinsically powerful but we are off-axis and only observe its side emission. In this paper, we provide some theoretical predictions of both the off-axis afterglow emission from a nearby SGRB under the assumption of a structured jet and the macronova signal from the ejecta of this GW-detectable BNS merger. From the properties of the afterglow emission, we could distinguish an off-axis normal SGRB from an intrinsically low-energy quasi-isotropic class. Furthermore, with follow-up multi-wavelength observations, a few parameters for BNS mergers (e.g. the medium density and the ejecta mass and velocity) would be constrained.

A hyperaccreting stellar-mass black hole has been proposed as the candidate central engine of gamma-ray bursts (GRBs). The rich observations of GRBs by \textit{Fermi} and \textit{Swift} make it possible to constrain the central engine model by comparing the model predications against data. This paper is dedicated to studying the temporal evolution of central engine parameters for both prompt emission and afterglow phases. We consider two jet launching mechanisms, i.e., $\nu\bar{\nu}$ annihilations and the Blandford-Znajek (BZ) processe, and obtain analytical solutions to these two models. We then investigate the black hole central engine parameters, such as the jet power, the dimensionless entropy $\eta$, and the central engine parameter $\mu_0=\eta (1+\sigma_0)$ (where $\sigma_0$ is the initial magnetization of the engine) at the base of the jet. The black hole may be spun up by accretion, or spun down by the BZ process, leaving imprints in GRB lightcurves. Usually, a BZ jet is more powerful and is likely responsible for the late time central engine activities. However, an initially non-spinning black hole central engine may first launch a thermal "fireball" via neutrino annihilations, and then launch a Poynting-flux-dominated jet via the BZ process. Multiple flares, giant bumps, and plateaus in GRB afterglows can be well produced as the result of late time accretion onto the black hole.

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.

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.

We present a spectroscopic survey of high-redshift, luminous galaxies over four square degrees on the sky, aiming to build a large and homogeneous sample of Ly$\alpha$ emitters (LAEs) at $z\approx5.7$ and 6.5, and Lyman-break galaxies (LBGs) at $5.5<z<6.8$. The fields that we choose to observe are well-studied, such as SXDS and COSMOS. They have deep optical imaging data in a series of broad and narrow bands, allowing efficient selection of galaxy candidates. Spectroscopic observations are being carried out using the multi-object spectrograph M2FS on the Magellan Clay telescope. M2FS is efficient to identify high-redshift galaxies, owing to its 256 optical fibers deployed over a circular field-of-view 30 arcmin in diameter. We have observed $\sim2.5$ square degrees. When the program is completed, we expect to identify more than 400 bright LAEs at $z\approx5.7$ and 6.5, and a substantial number of LBGs at $z\ge6$. This unique sample will be used to study a variety of galaxy properties and to search for large protoclusters. Furthermore, the statistical properties of these galaxies will be used to probe cosmic reionization. We describe the motivation, program design, target selection, and M2FS observations. We also outline our science goals, and present a sample of the brightest LAEs at $z\approx5.7$ and 6.5. This sample contains 32 LAEs with Ly$\alpha$ luminosities higher than 10$^{43}$ erg s$^{-1}$. A few of them reach $\ge3\times10^{43}$ erg s$^{-1}$, comparable to the two most luminous LAEs known at $z\ge6$, `CR7' and `COLA1'. These LAEs provide ideal targets to study extreme galaxies in the distant universe.

The Giant Radio Array for Neutrino Detection (GRAND) aims at detecting ultra-high energy extraterrestrial neutrinos via the extensive air showers induced by the decay of tau leptons created in the interaction of neutrinos under the Earth's surface. Consisting of an array of $\sim10^5$ radio antennas deployed over $\sim 2\times10^5\,\rm {km}^2$, GRAND plans to reach, for the first time, an all-flavor sensitivity of $\sim1.5\times10^{-10} \,\rm GeV\, cm^{-2} \,s^{-1}\, sr^{-1}$ above $5\times10^{17}$ eV and a sub-degree angular resolution, beyond the reach of other planned detectors. We describe here preliminary designs and simulation results, plans for the ongoing, staged approach to the construction of GRAND, and the rich research program made possible by GRAND's design sensitivity and angular resolution.

POLAR is a new satellite-born detector aiming to measure the polarization of an unprecedented number of Gamma-Ray Bursts in the 50-500 keV energy range. The instrument, launched on-board the Tiangong-2 Chinese Space lab on the 15th of September 2016, is designed to measure the polarization of the hard X-ray flux by measuring the distribution of the azimuthal scattering angles of the incoming photons. A detailed understanding of the polarimeter and specifically of the systematic effects induced by the instrument's non-uniformity are required for this purpose. In order to study the instrument's response to polarization, POLAR underwent a beam test at the European Synchrotron Radiation Facility in France. In this paper both the beam test and the instrument performance will be described. This is followed by an overview of the Monte Carlo simulation tools developed for the instrument. Finally a comparison of the measured and simulated instrument performance will be provided and the instrument response to polarization will be presented.

Integrated pulse profiles at 8.6~GHz obtained with the Shanghai Tian Ma Radio Telescope (TMRT) are presented for a sample of 26 pulsars. Mean flux densities and pulse width parameters of these pulsars are estimated. For eleven pulsars these are the first high-frequency observations and for a further four, our observations have a better signal-to-noise ratio than previous observations. For one (PSR J0742-2822) the 8.6~GHz profiles differs from previously observed profiles. A comparison of 19 profiles with those at other frequencies shows that in nine cases the separation between the outmost leading and trailing components decreases with frequency, roughly in agreement with radius-to-frequency mapping, whereas in the other ten the separation is nearly constant. Different spectral indices of profile components lead to the variation of integrated pulse profile shapes with frequency. In seven pulsars with multi-component profiles, the spectral indices of the central components are steeper than those of the outer components. For the 12 pulsars with multi-component profiles in the high-frequency sample, we estimate the core width using gaussian fitting and discuss the width-period relationship.

We carry out a systematical study of the spectral lag properties of 50 single-pulsed Gamma-Ray Bursts (GRBs) detected by Fermi/GBM. By dividing the light curves into multiple consecutive energy channels we provide a new measurement of the spectral lag which is independent on energy channel selections. We perform a detailed statistical study of our new measurements. We find two similar power-law energy dependencies of both the pulse arrival time and pulse width. Our new results on the power-law indices would favor the relativistic geometric effects for the origin of spectral lag. However, a complete theoretical framework that can fully account for the diverse energy dependencies of both arrival time and pulse width revealed in this work is still missing. We also study the spectral evolution behaviors of the GRB pulses. We find that the GRB pulse with negligible spectral lag would usually have a shorter pulse duration and would appear to have a "hardness-intensity tracking" (HIT) behavior and the GRB pulse with a significant spectral lag would usually have a longer pulse duration and would appear to have a "hard-to-soft" (HTS) behavior.

Adopting Schwarzschild's orbit-superposition technique, we construct a series of self-consistent galaxy models, embedded in the external field of galaxy clusters in the framework of Milgrom's MOdified Newtonian Dynamics. These models represent relatively massive ellipticals with a Hernquist radial profile at various distances from the cluster centre. Using $N$-body simulations, we perform a first analysis of these models and their evolution. We find that self-gravitating axisymmetric density models, even under a weak external field, lose their symmetry by instability and generally evolve to triaxial configurations. A kinematic analysis suggests that the instability originates from both box and non-classified orbits with low angular momentum. We also consider a self-consistent isolated system which is then placed in a strong external field and allowed to evolve freely. This model, just as the corresponding equilibrium model in the same external field, eventually settles to a triaxial equilibrium as well, but has a higher velocity radial anisotropy and is rounder. The presence of an external field in MOND universe generically predicts some lopsidedness of galaxy shapes.

A software system has been developed for the DArk Matter Particle Explorer (DAMPE) mission, a satellite-based experiment. The DAMPE software is mainly written in C++ and steered using Python script. This article presents an overview of the DAMPE offline software, including the major architecture design and specific implementation for simulation, calibration and reconstruction. The whole system has been successfully applied to DAMPE data analysis, based on which some results from simulation and beam test experiments are obtained and presented.

e-ASTROGAM (`enhanced ASTROGAM') is a breakthrough Observatory mission dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. In the largely unexplored MeV-GeV domain, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on Galactic ecosystems. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and the promise of eLISA. Keywords: High-energy gamma-ray astronomy, High-energy astrophysics, Nuclear Astrophysics, Compton and Pair creation telescope, Gamma-ray bursts, Active Galactic Nuclei, Jets, Outflows, Multiwavelength observations of the Universe, Counterparts of gravitational waves, Fermi, Dark Matter, Nucleosynthesis, Early Universe, Supernovae, Cosmic Rays, Cosmic antimatter.

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.

The DArk Matter Particle Explorer (DAMPE) is one of the four satellites within Strategic Pioneer Research Program in Space Science of the Chinese Academy of Science (CAS). DAMPE can detect electrons, photons in a wide energy range (5 GeV to 10 TeV) and ions up to iron (100GeV to 100 TeV). Silicon-Tungsten Tracker (STK) is one of the four subdetectors in DAMPE, providing photon-electron conversion, track reconstruction and charge identification for ions. Ion beam test was carried out in CERN with 60GeV/u Lead primary beams. Charge reconstruction and charge resolution of STK detectors were investigated.

A brief Chandra observation of the ultraluminous quasar, SDSS J010013.02+280225.8 at redshift 6.326, showed it to be a relatively bright, soft X-ray source with a count rate of about 1 ct/ks. In this paper we present results for the quasar from a 65ks XMM-Newton observation, which well constrains its spectral shape. The quasar is clearly detected with a total of $\sim$ 460 net counts in the 0.2-10 keV band. The spectrum is characterised by a simple power-law model with photon index of $\Gamma = 2.30^{+0.10}_{-0.10}$, and the intrinsic 2-10 keV luminosity is $3.14\times10^{45}$ erg $\text{s}^{-1}$. The 1 $\sigma$ upper limit to any intrinsic absorption column density is $N_{H} = 6.07\times 10^{22} {\text{cm}}^{-2}$. No significant iron emission lines were detected. We derive the X-ray-to-optical flux ratio $\alpha_{\text{ox}}$ of $-1.74\pm$0.01, consistent with the values found in other quasars of comparable ultraviolet luminosity. We did not detect significant flux variations either in the XMM-Newton exposure or between XMM-Newton and XMM-Newton observations, which are separated by $\sim$ 8 months. The X-ray observation enables the bolometric luminosity to be calculated after modelling the spectral energy distribution: the accretion rate is found to be sub-Eddington.

The Bose-Einstein condensation of $\alpha$ partciles in the multicomponent environment of dilute, warm nuclear matter is studied. We consider the cases of matter composed of light clusters with mass numbers $A\leq 4$ and matter that in addition these clusters contains $\isotope[56]{Fe}$ nuclei. We apply the quasiparticle gas model which treats clusters as bound states with infinite life-time and binding energies independent of temperature and density. We show that the $\alpha$ particles can form a condensate at low temperature $T\le 2$ MeV in such matter in the first case. When the $\isotope[56]{Fe}$ nucleus is added to the composition the cluster abundances are strongly modified at low temperatures, with an important implication that the $\alpha$ condensation at these temperatures is suppressed.

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.

POLAR is space-borne detector designed for a precise measurement of gamma-ray polarization of the prompt emissions of Gamma-Ray Bursts in the energy range 50 keV - 500 keV. POLAR is a compact Compton polarimeter consisting of 40$\times$ 40 plastic scintillator bars read out by 25 multi-anode PMTs. In May 2015, we performed a series of tests of the POLAR flight model with 100\% polarized x-rays beams at the European Synchrotron Radiation Facility beam-line ID11 aming to study thresholds, crosstalk between channels and responses of POLAR flight model to polarized X-ray beams. In this paper we present the data analysis method and some analysis results. According the results, POLAR FM has good polarimetric capabilities.

We present the first discoveries from a survey of $z\gtrsim6$ quasars using imaging data from the DECam Legacy Survey (DECaLS) in the optical, the UKIRT Deep Infrared Sky Survey (UKIDSS) and a preliminary version of the UKIRT Hemisphere Survey (UHS) in the near-IR, and ALLWISE in the mid-IR. DECaLS will image 9000 deg$^2$ of sky down to $z_{\rm AB}\sim23.0$, and UKIDSS and UHS, which will map the northern sky at $0<DEC<+60^{\circ}$, reaching $J_{\rm VEGA}\sim19.6$ (5-$\sigma$). The combination of these datasets allows us to discover quasars at redshift $z\gtrsim7$ and to conduct a complete census of the faint quasar population at $z\gtrsim6$. In this paper, we report on the selection method of our search, and on the initial discoveries of two new, faint $z\gtrsim6$ quasars and one new $z=6.63$ quasar in our pilot spectroscopic observations. The two new $z\sim6$ quasars are at $z=6.07$ and $z=6.17$ with absolute magnitudes at rest-frame wavelength 1450 \AA\ being $M_{1450}=-25.83$ and $M_{1450}=-25.76$, respectively. These discoveries suggest that we can find quasars close to or fainter than the break magnitude of the Quasar Luminosity Function (QLF) at $z\gtrsim6$. The new $z=6.63$ quasar has an absolute magnitude of $M_{1450}=-25.95$. This demonstrates the potential of using the combined DECaLS and UKIDSS/UHS datasets to find $z\gtrsim7$ quasars. Extrapolating from previous QLF measurements, we predict that these combined datasets will yield $\sim200$ $z\sim6$ quasars to $z_{\rm AB} < 21.5$, $\sim1{,}000$ $z\sim6$ quasars to $z_{\rm AB}<23$, and $\sim 30$ quasars at $z>6.5$ to $J_{\rm VEGA}<19.5$.

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

As a space-borne detector POLAR is designed to conduct hard X-ray polarization measurements of gamma-ray bursts on the statistically significant sample of events and with an unprecedented accuracy. During its development phase a number of tests, calibrations runs and verification measurements were carried out in order to validate instrument functionality and optimize operational parameters. In this article we present results on gain optimization togeter with verification data obtained in the course of broad laboratory and environmental tests. In particular we focus on exposures to the $^{137}$Cs radioactive source and determination of the gain dependence on the high voltage for all 1600 detection channels of the polarimeter. Performance of the instrument is described in detail with respect to the dynamic range, energy resolution and temperature dependence. Gain optimization algorithms and response non-uniformity studies are also broadly discussed. Results presented below constitute important parts for development of the POLAR calibration and operation database.

The emission geometries, e.g. the emission region height, the beam shape, and radius-to-frequency mapping, are important predictions of pulsar radiation model. The multi-band radio observations carry such valuable information. In this paper, we study two bright pulsars, (PSRs B0329+54 and B1642-03) and observe them in high frequency (2.5 GHz, 5 GHz, and 8 GHz). The newly acquired data together with historical archive provide an atlas of multi-frequency profiles spanning from 100 MHz to 10 GHz. We study the frequency evolution of pulse profiles and the radiation regions with the these data. We firstly fit the pulse profiles with Gaussian functions to determine the phase of each component, and then calculate the radiation altitudes of different emission components and the radiation regions. We find that the inverse Compton scattering (ICS) model can reproduce the radiation geometry of these two pulsars. But for PSR B0329+54 the radiation can be generated in either annular gap (AG) or core gap (CG), while the radiation of PSR B1642-03 can only be generated in the CG. This difference is caused by the inclination angle and the impact angle of these two pulsars. The relation of beaming angle (the angle between the radiation direction and the magnetic axis) and the radiation altitudes versus frequency is also presented by modelling the beam-frequency evolution in the ICS model. The multi-band pulse profiles of these two pulsars can be described well by the ICS model combined with the CG and AG.

We present initial results from the first systematic survey of luminous $z\sim 5.5$ quasars. Quasars at $z \sim$ 5.5, the post-reionization epoch, are crucial tools to explore the evolution of intergalactic medium, quasar evolution and the early super-massive black hole growth. However, it has been very challenging to select quasars at redshifts 5.3 $\le z \le$ 5.7 using conventional color selections, due to their similar optical colors to late-type stars, especially M dwarfs, resulting in a glaring redshift gap in quasar redshift distributions. We develop a new selection technique for $z \sim$ 5.5 quasars based on optical, near-IR and mid-IR photometric data from Sloan Digital Sky Survey (SDSS), UKIRT InfraRed Deep Sky Surveys - Large Area Survey (ULAS), VISTA Hemisphere Survey (VHS) and Wide field Infrared Survey Explorer (WISE). From our pilot observations in SDSS-ULAS/VHS area, we have discovered 15 new quasars at 5.3 $\le z \le$ 5.7 and 6 new lower redshift quasars, with SDSS z band magnitude brighter than 20.5. Including other two $z \sim$ 5.5 quasars already published in our previous work, we now construct an uniform quasar sample at 5.3 $\le z \le$ 5.7 with 17 quasars in a $\sim$ 4800 square degree survey area. For further application in a larger survey area, we apply our selection pipeline to do a test selection by using the new wide field J band photometric data from a preliminary version of the UKIRT Hemisphere Survey (UHS). We successfully discover the first UHS selected $z \sim$ 5.5 quasar.

Very few low-ionization broad absorption line (LoBAL) QSOs have been found at high redshifts to date. One high-redshift LoBAL QSO, J0122+1216, was recently discovered at the Lijiang 2.4-m Telescope with an initial redshift determination of 4.76. Aiming to investigate its physical properties, we carried out follow-up observations in the optical and near-IR spectroscopy. Near-IR spectra from UKIRT and P200 confirms that it is a LoBAL, with a new redshift determination of $4.82\pm0.01$ based on the \mgii~ emission-line. The new \mgii~ redshift determination reveals strong blueshifts and asymmetry of the high-ionization emission lines. We estimated a black hole mass of $\sim 2.3\times 10^9 M_\odot$ and Eddington ratio of $\sim 1.0$ according to the empirical \mgii-based single-epoch relation and bolometric correction factor. It is possible that strong outflows are the result of an extreme quasar environment driven by the high Eddington ratio. A lower limit on the outflowing kinetic power ($>0.9\% L_{Edd}$) was derived from both emission and absorption lines, indicating these outflows play a significant role in the feedback process to regulate the growth of its black hole as well as host galaxy evolution.

The Giant Radio Array for Neutrino Detection (GRAND) is a planned array of ~200 000 radio antennas deployed over ~200 000 km2 in a mountainous site. It aims primarly at detecting high-energy neutrinos via the observation of extensive air showers induced by the decay in the atmosphere of taus produced by the interaction of cosmic neutrinos under the Earth surface. GRAND aims at reaching a neutrino sensitivity of 5.10$^{11}$ E$^{-2}$ GeV$^{-1}$cm$^{-2}$s$^{-1}$sr$^{-1}$ above 3.10$^{16}$ eV. This ensures the detection of cosmogenic neutrinos in the most pessimistic source models, and ~50 events per year are expected for the standard models. The instrument will also detect UHECRs and possibly FRBs. Here we show how our preliminary design should enable us to reach our sensitivity goals, and discuss the steps to be taken to achieve GRAND.

The e-ASTROGAM is a gamma-ray space mission to be proposed as the M5 Medium-size mission of the European Space Agency. It is dedicated to the observation of the Universe with unprecedented sensitivity in the energy range 0.2 - 100 MeV, extending up to GeV energies, together with a groundbreaking polarization capability. It is designed to substantially improve the COMPTEL and Fermi sensitivities in the MeV-GeV energy range and to open new windows of opportunity for astrophysical and fundamental physics space research. e-ASTROGAM will operate as an open astronomical observatory, with a core science focused on (1) the activity from extreme particle accelerators, including gamma-ray bursts and active galactic nuclei and the link of jet astrophysics to the new astronomy of gravitational waves, neutrinos, ultra-high energy cosmic rays, (2) the high-energy mysteries of the Galactic center and inner Galaxy, including the activity of the supermassive black hole, the Fermi Bubbles, the origin of the Galactic positrons, and the search for dark matter signatures in a new energy window; (3) nucleosynthesis and chemical evolution, including the life cycle of elements produced by supernovae in the Milky Way and the Local Group of galaxies. e-ASTROGAM will be ideal for the study of high-energy sources in general, including pulsars and pulsar wind nebulae, accreting neutron stars and black holes, novae, supernova remnants, and magnetars. And it will also provide important contributions to solar and terrestrial physics. The e-ASTROGAM telescope is optimized for the simultaneous detection of Compton and pair-producing gamma-ray events over a large spectral band. It is based on a very high technology readiness level for all subsystems and includes many innovative features for the detectors and associated electronics.

A vast wealth of literature exists on the topic of rocket trajectory optimisation, particularly in the area of interplanetary trajectories due to its relevance today. Studies on optimising interstellar and intergalactic trajectories are usually performed in flat spacetime using an analytical approach, with very little focus on optimising interstellar trajectories in a general relativistic framework. This paper examines the use of low-acceleration rockets to reach galactic destinations in the least possible time, with a genetic algorithm being employed for the optimisation process. The fuel required for each journey was calculated for various types of propulsion systems to determine the viability of low-acceleration rockets to colonise the Milky Way. The results showed that to limit the amount of fuel carried on board, an antimatter propulsion system would likely be the minimum technological requirement to reach star systems tens of thousands of light years away. However, using a low-acceleration rocket would require several hundreds of thousands of years to reach these star systems, with minimal time dilation effects since maximum velocities only reached about 0.2c. Such transit times are clearly impractical, and thus, any kind of colonisation using low acceleration rockets would be difficult. High accelerations, on the order of 1g, are likely required to complete interstellar journeys within a reasonable time frame, though they may require prohibitively large amounts of fuel. So for now, it appears that humanity's ultimate goal of a galactic empire may only be possible at significantly higher accelerations, though the propulsion technology requirement for a journey that uses realistic amounts of fuel remains to be determined.

We report Very Long Baseline Array (VLBA) observations of the 1.5 GHz radio continuum emission of the {\it z=6.326} quasar SDSS J010013.02+280225.8 (hereafter J0100+2802). J0100+2802 is, by far the most optically luminous, and radio-quiet quasar with the most massive black hole known at z>6. The VLBA observations have a synthesized beam size of 12.10 mas $\times$5.36 mas (FWHM), and detected the radio continuum emission from this object with a peak surface brightness of 64.6+/-9.0 micro-Jy/beam and a total flux density of 88+/-19 micro-Jy. The position of the radio peak is consistent with that from SDSS in the optical and Chandra in the X-ray. The radio source is marginally resolved by the VLBA observations. A 2-D Gaussian fit to the image constrains the source size to $\rm (7.1+/-3.5) mas x (3.1+/-1.7) mas. This corresponds to a physical scale of (40+/-20) pc x (18+/-10) pc. We estimate the intrinsic brightness temperature of the VLBA source to be T_{B}=(1.6 +/- 1.2) x 10^{7} K. This is significantly higher than the maximum value in normal star forming galaxies, indicating an AGN origin for the radio continuum emission. However, it is also significantly lower than the brightness temperatures found in highest redshift radio-loud quasars. J0100+2802 provides a unique example to study the radio activity in optically luminous and radio quiet active galactic nuclei in the early universe. Further observations at multiple radio frequencies will accurately measure the spectral index and address the dominant radiation mechanism of the radio emission.

We report new IRAM/PdBI, JCMT/SCUBA-2, and VLA observations of the ultraluminous quasar SDSSJ010013.02+280225.8 (hereafter, J0100+2802) at z=6.3, which hosts the most massive supermassive black hole (SMBH) of 1.24x10^10 Msun known at z>6. We detect the [C II] 158 $\mu$m fine structure line and molecular CO(6-5) line and continuum emission at 353 GHz, 260 GHz, and 3 GHz from this quasar. The CO(2-1) line and the underlying continuum at 32 GHz are also marginally detected. The [C II] and CO detections suggest active star formation and highly excited molecular gas in the quasar host galaxy. The redshift determined with the [C II] and CO lines shows a velocity offset of ~1000 km/s from that measured with the quasar Mg II line. The CO (2-1) line luminosity provides direct constraint on the molecular gas mass which is about (1.0+/-0.3)x10^10 Msun. We estimate the FIR luminosity to be (3.5+/-0.7)x10^12 Lsun, and the UV-to-FIR spectral energy distribution of J0100+2802 is consistent with the templates of the local optically luminous quasars. The derived [C II]-to-FIR luminosity ratio of J0100+2802 is 0.0010+/-0.0002, which is slightly higher than the values of the most FIR luminous quasars at z~6. We investigate the constraint on the host galaxy dynamical mass of J0100+2802 based on the [C II] line spectrum. It is likely that this ultraluminous quasar lies above the local SMBH-galaxy mass relationship, unless we are viewing the system at a small inclination angle.

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.

We work on a GRB sample whose initial Lorentz factors ($\Gamma_0$) are constrained with the afterglow onset method and the jet opening angles ($\theta_{\rm j}$) are determined by the jet break time. We confirm the $\Gamma_0$ - $E_{\gamma,\rm iso}$ correlation by Liang et al. (2010), and the $\Gamma_0$ - $L_{\gamma,\rm iso}$ correlation by L{\"u} et al. (2012). Furthermore, we find correlations between $\Gamma_0$ and the beaming corrected $\gamma$-ray energy ($E_{\gamma}$) and mean $\gamma$-ray luminosity ($L_{\gamma}$). By also including the kinetic energy of the afterglow, we find rough correlations (with larger scatter) between $\Gamma_0$ and the total ($\gamma$-ray plus kinetic) energy and the total mean luminosity, both for isotropic values and beaming corrected values: these correlations allow us to test the data with GRB central engine models. Limiting our sample to the GRBs that likely have a black hole central engine, we compare the data with theoretical predictions of two types of jet launching mechanisms from BHs, i.e. the non-magnetized $\nu \bar\nu$ - annihilation mechanism, and the strongly magnetized Blandford-Znajek (BZ) mechanism. We find that the data are more consistent with the latter mechanism, and discuss the implications of our findings for GRB jet composition.

Gamma-ray polarimetry is a new powerful tool to study the processes responsible for the emission from astrophysical sources and the environments in which this emission takes place. Few successful polarimetric measurements have however been performed thus far in the gamma-ray energy band due to the difficulties involved. POLAR is a dedicated polarimeter designed to perform high precision measurements of the polarization of the emission from gamma-ray burst in the 50-500 keV energy range. This new polarimeter is expected to detect approximately 50 gamma-ray bursts per year while performing high precision polarization measurements on approximately 10 bursts per year. The instrument was launched into lower earth orbit as part of the second Chinese space lab, the Tiangong-2, on September 15th 2016 and has been taking data successfully since being switched on one week after. The instrument uses a segmented scintillator array consisting of 1600 plastic scintillator bars, read out by 25 flat-panel multi-anode photomultipliers, to measure the Compton scattering angles of incoming photons. The small segmentation and relatively large uniform effective area allow the instrument to measure the polarization of a large number of transient events, such as gamma-ray bursts, with an unprecedented precision during its two year life-time. The final flight model underwent detailed calibration prior to launch as well as intensive space qualification tests, a summary of which will be presented in this paper. The instrument design will be discussed first followed by an overview of the on-ground tests, finally the in-orbit behavior as measured during the first weeks of the mission will be presented.

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.

Recently, DeLaunay et al. (2016) discovered a gamma-ray transient, Swift J0644.5-5111, associated with the fast radio burst (FRB) 131104. They also reported follow-up broadband observations beginning two days after the FRB and provided upper limits on a putative afterglow of this transient. In this paper, we show that if such a transient drives a relativistic shock as in a cosmological gamma-ray burst (GRB), these upper limits are consistent with an environment of which density is much less than that of an interstellar medium but typical for the outskirts' density of a galaxy when the typical values of three microphysical parameters of the shock are taken. This appears to be inconsistent with the catastrophic event models in which the central engine of Swift J0644.5-5111 is surrounded by an interstellar medium, but together with the properties of the gamma-ray transient, favors the binary neutron star merger origin. We further constrain the physical parameters of the postmerger object by assuming that Swift J0644.5-5111 results from internal dissipation of a spinning-down pulsar wind, and we find that the postmerger object is an ultra-strongly magnetized, very rapidly rotating pulsar. This merger event should have given birth to a gravitational wave burst, an FRB, and a short GRB or an extended X-ray/gamma-ray emission if a relativistic jet of the GRB is missed. Such "triplets" would be testable in the near future with the advanced LIGO and Virgo gravitational-wave observatories.

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.

We present the discovery of nine quasars at $z\sim6$ identified in the Sloan Digital Sky Survey (SDSS) imaging data. This completes our survey of $z\sim6$ quasars in the SDSS footprint. Our final sample consists of 52 quasars at $5.7<z\le6.4$, including 29 quasars with $z_{\rm AB}\le20$ mag selected from 11,240 deg$^2$ of the SDSS single-epoch imaging survey (the main survey), 10 quasars with $20\le z_{\rm AB}\le20.5$ selected from 4223 deg$^2$ of the SDSS overlap regions (regions with two or more imaging scans), and 13 quasars down to $z_{\rm AB}\approx22$ mag from the 277 deg$^2$ in Stripe 82. They span a wide luminosity range of $-29.0\le M_{1450}\le-24.5$. This well-defined sample is used to derive the quasar luminosity function (QLF) at $z\sim6$. After combining our SDSS sample with two faint ($M_{1450}\ge-23$ mag) quasars from the literature, we obtain the parameters for a double power-law fit to the QLF. The bright-end slope $\beta$ of the QLF is well constrained to be $\beta=-2.8\pm0.2$. Due to the small number of low-luminosity quasars, the faint-end slope $\alpha$ and the characteristic magnitude $M_{1450}^{\ast}$ are less well constrained, with $\alpha=-1.90_{-0.44}^{+0.58}$ and $M^{\ast}=-25.2_{-3.8}^{+1.2}$ mag. The spatial density of luminous quasars, parametrized as $\rho(M_{1450}<-26,z)=\rho(z=6)\,10^{k(z-6)}$, drops rapidly from $z\sim5$ to 6, with $k=-0.72\pm0.11$. Based on our fitted QLF and assuming an IGM clumping factor of $C=3$, we find that the observed quasar population cannot provide enough photons to ionize the $z\sim6$ IGM at $\sim90$\% confidence. Quasars may still provide a significant fraction of the required photons, although much larger samples of faint quasars are needed for more stringent constraints on the quasar contribution to reionization.

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 present a catalog of 624 radio sources detected around the North Celestial Pole (NCP) with the 21 Centimeter Array (21CMA), a radio interferometer dedicated to the statistical measurement of the epoch of reionization (EoR). The data are taken from a 12 h observation made on 2013 April 13, with a frequency coverage from 75 to 175 MHz and an angular resolution of ~ 4 arcmin. The catalog includes flux densities at eight sub-bands across the 21CMA bandwidth and provides the in-band spectral indices for the detected sources. To reduce the complexity of interferometric imaging from the so-called "w" term and ionospheric effects, the present analysis are restricted to the east-west baselines within 1500 m only. The 624 radio sources are found within 5 degrees around the NCP down to ~ 0.1 Jy. Our source counts are compared, and also exhibit a good agreement, with deep low-frequency observations made recently with the GMRT and MWA. In particular, for fainter radio sources below ~ 1 Jy, we find a flattening trend of source counts towards lower frequencies. While the thermal noise (~0.4 mJy) is well controlled to below the confusion limit, the dynamical range (~10^4) and sensitivity of current 21CMA imaging is largely limited by calibration and deconvolution errors, especially the grating lobes of very bright sources, such as 3C061.1, in the NCP field which result from the regular spacings of the 21CMA. We note that particular attention should be paid to the extended sources, and their modeling and removals may constitute a large technical challenge for current EoR experiments. Our analysis may serve as a useful guide to design of next generation low-frequency interferometers like the Square Kilometre Array.

The X-ray source CXO J133815.6+043255 has counterparts in the UV, optical, and radio bands. Based on the multi-band investigations, it has been recently proposed by Kim et al. (2015) as a rarely-seen off-nucleus ultraluminous X-ray (ULX) source with a black hole mass of >= 10^4 solar mass in the nearby Seyfert galaxy NGC 5252. To explore its radio properties at very high angular resolution, we performed very long-baseline interferometry (VLBI) observations with the European VLBI Network (EVN) at 1.7 GHz. We find that the radio counterpart is remarkably compact among the known ULXs. It does not show a resolved structure with a resolution of a few milliarcsecond (mas), and the total recovered flux density is comparable to that measured in earlier sub-arcsecond-resolution images. The compact radio structure, the relatively flat spectrum, and the high radio luminosity are consistent with a weakly accreting supermassive black hole in a low-luminosity active galactic nucleus. The nucleus of NGC 5252 itself has similar radio properties. We argue that the system represents a relatively rare pair of active galactic nuclei, where both components emit in the radio.

We present the X-ray luminosity (L) versus dynamical mass (M) relation for 63 nearby clusters in the HIFLUGCS. The luminosity measurements are obtained based on ~1.3 Ms of clean XMM data and ROSAT pointed observations. The masses are estimated using optical spectroscopic redshifts of 13647 cluster galaxies in total. Given sufficient numbers of member galaxies in computing the dynamical masses, the L-M relations agree between the disturbed and undisturbed clusters. The cool-core clusters still dominate the scatter in the L-M relation even when a core corrected X-ray luminosity is used, which indicates that the scatter mainly reflects the structure formation history of the clusters. As shown by the clusters with a small number of redshifts, the dynamical masses can be underestimated leading to a biased scaling relation. To investigate the potential of spectroscopic surveys to follow up high-redshift galaxy clusters/groups observed in X-ray surveys for the identifications and mass calibrations, we carried out Monte-Carlo re-sampling of the cluster galaxy redshifts and calibrated the uncertainties of the redshift and dynamical mass estimates when only reduced numbers of galaxy redshifts per cluster are available. The re-sampling considers the SPIDERS and 4MOST configurations, designed for the follow-up of the eROSITA clusters, and was carried out for each cluster at the actual cluster redshift as well as at z=0.2, 0.4, 0.6, and 0.8. For following up very distant cluster/groups, we carried out the mass calibration based on the re-sampling with only 10zs/cluster, and redshift calibration based on the re-sampling with only 5zs/cluster and 10zs/cluster, respectively. Our results demonstrate the power of combining upcoming X-ray and optical spectroscopic surveys for mass calibration. The scatter in the dynamical mass estimates for the clusters with at least ten members is within 50%.

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.

The high Antarctic plateau provides exceptional conditions for conducting infrared observations of the cosmos on account of the cold, dry and stable atmosphere above the ice surface. This paper describes the scientific goals behind the first program to examine the time-varying universe in the infrared from Antarctica - the Kunlun Infrared Sky Survey (KISS). This will employ a small (50 cm aperture) telescope to monitor the southern skies in the 2.4um Kdark window from China's Kunlun station at Dome A, on the summit of the Antarctic plateau, through the uninterrupted 4-month period of winter darkness. An earlier paper discussed optimisation of the Kdark filter for the best sensitivity (Li et al 2016). This paper examines the scientific program for KISS. We calculate the sensitivity of the camera for the extrema of observing conditions that will be encountered. We present the parameters for sample surveys that could then be carried out for a range of cadences and sensitivities. We then discuss several science programs that could be conducted with these capabilities, involving star formation, brown dwarfs and hot Jupiters, exoplanets around M dwarfs, the terminal phases of stellar evolution, discovering fast transients as part of multi-wavelength campaigns, embedded supernova searches, reverberation mapping of active galactic nuclei, gamma ray bursts and the detection of the cosmic infrared background. Accepted for publication in PASA, 04/08/16.

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 selected 82 u-band variable objects based on the u-band photometry data from SCUSS and SDSS, in the field of LAMOST Complete Spectroscopic Survey of Pointing Area at Southern Galactic Cap. The magnitude variation of the targets is restricted to larger than 0.2 mag and limiting magnitude down to 19.0 mag in u-band. According to the spectra from LAMOST, there are 11 quasars with red-shift between 0.4 and 1.8, 60 variable stars and 11 unidentified targets. The variable stars include one active M-dwarf with a series of Balmer emission lines, seven Horizontal Branch stars containing six RR Lyrae stars matching with SIMBAD, and one giant, one AGB star and two RR Lyrae candidates by different colour selections. All these variable stars mainly locate near the main sequence in the g-r verse u-g diagram. The quasars are well distinguished from stars by both u-g colour and variation in u-band.

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of about 0.9 m^2 and 0.6 m^2 at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering <180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of about 3.4 m^2, between 6 and 10 keV, and spectral resolution <250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm^2 at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees FoV. The eXTP international consortium includes mostly major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025.

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of about 0.9 m^2 and 0.6 m^2 at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering <180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of about 3.4 m^2, between 6 and 10 keV, and spectral resolution <250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm^2 at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees FoV. The eXTP international consortium includes mostly major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025.

In this review we present a progress report of the connection between long-duration gamma-ray bursts (GRBs) and their accompanying supernovae (SNe). The analysis is from the point of view of an observer, with much of the emphasis placed on how observations, and the modelling of observations, have constrained what we known about GRB-SNe. We discuss their photometric and spectroscopic properties, their role as cosmological probes, including their measured luminosity$-$decline relationships, and how they can be used to measure the Hubble constant. We present a statistical analysis of their bolometric properties, and use this to determine the properties of the "average" GRB-SNe: which has a kinetic energy of $E_{\rm K} \approx 2.5\times10^{52}$ erg, an ejecta mass of $M_{\rm ej} \approx 6$ M$_{\odot}$, a nickel mass of $M_{\rm Ni} \approx 0.4$ M$_{\odot}$, a peak photospheric velocity of $v_{\rm ph} \approx 21,000$ km s$^{-1}$, a peak bolometric luminosity of $L_{\rm p} \approx 1\times10^{43}$ erg s$^{-1}$, and it reaches peak bolometric light in $t_{\rm p} \approx 13$ days. We discuss their geometry, consider the various physical processes that are thought to power the luminosity of GRB-SNe, and whether differences exist between GRB-SNe and the SNe associated with ultra-long duration GRBs. We discuss how observations of the environments of GRB-SNe further constrain the physical properties of their progenitor stars, and give an overview of the current theoretical paradigms of their suspected central engines. We also present an overview of the radioactively powered transients that have been photometrically associated with short-duration GRBs. We conclude the review by discussing what additional research is needed to further our understanding of GRB-SNe, in particular the role of binary-formation channels and the connection of GRB-SNe with superluminous SNe (abridged).

This is the second paper in a series on a new luminous z ~ 5 quasar survey using optical and near-infrared colors. Here we present a new determination of the bright end of the quasar luminosity function (QLF) at z ~ 5. Combined our 45 new quasars with previously known quasars that satisfy our selections, we construct the largest uniform luminous z ~ 5 quasar sample to date, with 99 quasars in the range 4.7 <= z < 5.4 and -29 < M1450 <= -26.8, within the Sloan Digital Sky Survey (SDSS) footprint. We use a modified 1/Va method including flux limit correction to derive a binned QLF, and we model the parametric QLF using maximum likelihood estimation. With the faint-end slope of the QLF fixed as alpha = -2.03 from previous deeper samples, the best fit of our QLF gives a flatter bright end slope beta = -3.58+/-0.24 and a fainter break magnitude M*1450 = -26.98+/-0.23 than previous studies at similar redshift. Combined with previous work at lower and higher redshifts, our result is consistent with a luminosity evolution and density evolution (LEDE) model. Using the best fit QLF, the contribution of quasars to the ionizing background at z ~ 5 is found to be 18% - 45% with a clumping factor C of 2 - 5. Our sample suggests an evolution of radio loud fraction with optical luminosity but no obvious evolution with redshift.

One of the cornerstones of modern physics is Einstein's special relativity, with its constant speed of light and zero photon mass assumptions. Constraint on the rest mass m_{\gamma} of photons is a fundamental way to test Einstein's theory, as well as other essential electromagnetic and particle theories. Since non-zero photon mass can give rise to frequency-(or energy-) dependent dispersions, measuring the time delay of photons with different frequencies emitted from explosive astrophysical events is an important and model-independent method to put such a constraint. The cosmological gamma-ray bursts (GRBs), with short time scales, high redshifts as well as broadband prompt and afterglow emissions, provide an ideal testbed for m_{\gamma} constraints. In this paper we calculate the upper limits of the photon mass with GRB early time radio afterglow observations as well as multi-band radio peaks, thus improve the results of Schaefer (1999) by nearly half an order of magnitude.

Almost all superluminous supernovae (SLSNe) whose peak magnitudes are $\lesssim -21$ mag can be explained by the $^{56}$Ni-powered model, magnetar-powered (highly magnetized pulsar) model or ejecta-circumstellar medium (CSM) interaction model. Recently, iPTF13ehe challenges these energy-source models, because the spectral analysis shows that $\sim 2.5M_\odot$ of $^{56}$Ni have been synthesized but are inadequate to power the peak bolometric emission of iPTF13ehe, while the rebrightening of the late-time light-curve (LC) and the H$\alpha$ emission lines indicate that the ejecta-CSM interaction must play a key role in powering the late-time LC. Here we propose a triple-energy-source model, in which a magnetar together with some amount ($\lesssim 2.5M_\odot$) of $^{56}$Ni may power the early LC of iPTF13ehe while the late-time rebrightening can be quantitatively explained by an ejecta-CSM interaction. Furthermore, we suggest that iPTF13ehe is a genuine core-collapse supernova rather than a pulsational pair-instability supernova candidate. Further studies on similar SLSNe in the future would eventually shed light on their explosion and energy-source mechanisms.

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.

While the extended main-sequence turn-offs (eMSTOs) found in almost all 1--2 Gyr-old star clusters in the Magellanic Clouds are often explained by postulating extended star-formation histories, the tight subgiant branches (SGBs) seen in some clusters challenge this popular scenario. Puzzlingly, the SGB of the eMSTO cluster NGC 419 is significantly broader at bluer than at redder colors. We carefully assess and confirm the reality of this observational trend. If we would assume that the widths of the features in color--magnitude space were entirely owing to a range in stellar ages, the star-formation histories of the eMSTO stars and the blue SGB region would be significantly more prolonged than that of the red part of the SGB. This cannot be explained by assuming an internal age spread. We show that rotational deceleration of a population of rapidly rotating stars, a currently hotly debated alternative scenario, naturally explains the observed trend along the SGB. Our analysis shows that a `converging' SGB could be produced if the cluster is mostly composed of rapidly rotating stars that slow down over time owing to the conservation of angular momentum during their evolutionary expansion from main-sequence turn-off stars to red giants.

Very recently Spitler et al. (2016) and Scholz et al. (2016) reported their detections of sixteen additional bright bursts from the direction of the fast radio burst (FRB) 121102. This repeating FRB is inconsistent with all the catastrophic event models put forward previously for hypothetically non-repeating FRBs. Here we propose a different model, in which highly magnetized pulsars travel through asteroid belts of other stars. We show that a repeating FRB could originate from such a pulsar encountering lots of asteroids in the belt. During each pulsar-asteroid impact, an electric field induced outside the asteroid has such a large component parallel to the stellar magnetic field that electrons are torn off the asteroidal surface and accelerated to ultra-relativistic energies instantaneously. Subsequent movement of these electrons along magnetic field lines will cause coherent curvature radiation, which can account for all the properties of an FRB. In addition, this model can self-consistently explain the typical duration, luminosity, and repetitive rate of the seventeen bursts of FRB 121102. The predicted occurrence rate of repeating FRB sources may imply that our model would be testable in the next few years.

The Silicon Tracker (STK) is a detector of the DAMPE satellite to measure the incidence direction of high energy cosmic ray. It consists of 6 X-Y double layers of silicon micro-strip detectors with 73,728 readout channels. It's a great challenge to readout the channels and process the huge volume of data in the critical space environment. 1152 Application Specific Integrated Circuits (ASIC) and 384 ADCs are adopted to readout the detector channels. The 192 Tracker Front-end Hybrid (TFH) modules and 8 identical Tracker Readout Board (TRB) modules are designed to control and digitalize the front signals. In this paper, the design of the readout electronics for STK and its performance will be presented in detail.

The inner boundary of a black hole accretion disk is often set to the marginally stable circular orbit (or the innermost stable circular orbit, ISCO) around the black hole. It is important for the theories of black hole accretion disks and their applications to astrophysical black hole systems. Traditionally, the marginally stable circular orbit is obtained by considering the equatorial motion of a test particle around a black hole. However, in reality the accretion flow around black holes consists of fluid, in which the pressure often plays an important role. Here we consider the influence of fluid pressure on the location of marginally stable circular orbit around black holes. It is found that when the temperature of the fluid is so low that the thermal energy of a particle is much smaller than its rest energy, the location of marginally stable circular orbit is almost the same as that in the test particle case. However, we demonstrate that in some special cases the marginally stable circular orbit can be different when the fluid pressure is large and the thermal energy becomes non-negligible comparing with the rest energy. We present our results for both the cases of non-spinning and spinning black holes. The influences of our results on the black hole spin parameter measurement in X-ray binaries and the energy release efficiency of accretion flows around black holes are discussed.

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

The X-ray afterglows of almost one half of gamma-ray bursts have been discovered to have a shallow decay phase by the {\em Swift} satellite, whose origin remains mysterious. Two main models have been proposed to explain this phase, relativistic wind bubbles (RWBs) and structured ejecta, which could originate from millisecond magnetars and rapidly-rotating black holes, respectively. Based on these models, we here investigate polarization evolution in the shallow decay phase of X-ray and optical afterglows. We find that in the RWB model, a significant bump of the polarization degree evolution curve appears during the shallow decay phase of both optical and X-ray afterglows, while the polarization position angle changes its direction by $90^\circ$ abruptly. In the structured ejecta model, however, the polarization degree does not evolve significantly during the shallow decay phase of afterglows, no matter whether the magnetic field configuration in the ejecta is random or globally large-scale. Therefore, we conclude that these two models for the shallow decay phase and relevant central engines would be testable with future polarization observations.

We report exploratory \chandra\ observation of the ultraluminous quasar SDSS J010013.02+280225.8 at redshift 6.30. The quasar is clearly detected by \chandra\ with a possible component of extended emission. The rest-frame 2-10 keV luminosity is 9.0$^{+9.1}_{-4.5}$ $\times$ 10$^{45}$ erg s$^{-1}$ with inferred photon index of $\Gamma$ = 3.03$^{+0.78}_{-0.70}$. This quasar is X-ray bright, with inferred X-ray-to-optical flux ratio \aox\ $=-1.22^{+0.07}_{-0.05}$, higher than the values found in other quasars of comparable ultraviolet luminosity. The properties inferred from this exploratory observation indicate that this ultraluminous quasar might be growing with super-Eddington accretion and probably viewed with small inclination angle. Deep X-ray observation will help to probe the plausible extended emission and better constraint the spectral features for this ultraluminous quasar.

Optical re-brightenings in the afterglows of some gamma-ray bursts (GRBs) are unexpected within the framework of the simple external shock model. While it has been suggested that the central engines of some GRBs are newly born magnetars, we aim to relate the behaviors of magnetars to the optical re-brightenings. A newly born magnetar will lose its rotational energy in the form of Poynting-flux, which may be converted into a wind of electron-positron pairs through some magnetic dissipation processes. As proposed by Dai (2004), this wind will catch up with the GRB outflow and a long-lasting reverse shock would form. By applying this scenario to GRB afterglows, we find that the reverse shock propagating back into the electron-positron wind can lead to an observable optical re-brightening and a simultaneous X-ray plateau (or X-ray shallow decay). In our study, we select four GRBs, i.e., GRB 080413B, GRB 090426, GRB 091029, and GRB 100814A, of which the optical afterglows are well observed and show clear re-brightenings. We find that they can be well interpreted. In our scenario, the spin-down timescale of the magnetar should be slightly smaller than the peak time of the re-brightening, which can provide a clue to the characteristics of the magnetar.

In this work, we present a study of 207 quasars selected from the Sloan Digital Sky Survey quasar catalogs and the Herschel Stripe 82 survey. Quasars within this sample are high luminosity quasars with a mean bolometric luminosity of $10^{46.4}$ erg s$^{-1}$. The redshift range of this sample is within $z<4$, with a mean value of $1.5\pm0.78$. Because we only selected quasars that have been detected in all three Herschel-SPIRE bands, the quasar sample is complete yet highly biased. Based on the multi-wavelength photometric observation data, we conducted a spectral energy distribution (SED) fitting through UV to FIR. Parameters such as active galactic nucleus (AGN) luminosity, FIR luminosity, stellar mass, as well as many other AGN and galaxy properties are deduced from the SED fitting results. The mean star formation rate (SFR) of the sample is 419 $M_{\odot}$ yr$^{-1}$ and the mean gas mass is $\sim 10^{11.3}$ $M_{\odot}$. All these results point to an IR luminous quasar system. Comparing with star formation main sequence (MS) galaxies, at least 80 out of 207 quasars are hosted by starburst galaxies. It supports the statement that luminous AGNs are more likely to be associated with major mergers. The SFR increases with the redshift up to $z=2$. It is correlated with the AGN bolometric luminosity, where $L_{\rm FIR} \propto L_{\rm Bol}^{0.46\pm0.03}$. The AGN bolometric luminosity is also correlated with the host galaxy mass and gas mass. Yet the correlation between $L_{\rm FIR}$ and $L_{\rm Bol}$ has higher significant level, implies that the link between AGN accretion and the SFR is more primal. The $M_{\rm BH}/M_{\ast}$ ratio of our sample is 0.02, higher than the value 0.005 in the local Universe. It might indicate an evolutionary trend of the $M_{\rm BH} - M_{\ast}$ scaling relation.

The shallow decay phase and flares in the afterglows of gamma-ray bursts (GRBs) is widely believed to be associated with the later activation of central engine. Some models of energy injection involve with a continuous energy flow since the GRB trigger time, such as the magnetic dipole radiation from a magnetar. However, in the scenario involving with a black hole accretion system, the energy flow from the fall-back accretion may be delayed for a fall-back time $\sim t_{\rm fb}$. Thus we propose a delayed energy injection model, the delayed energy would cause a notable rise to the Lorentz factor of the external shock, which will "generate" a bump in the multiple band afterglows. If the delayed time is very short, our model degenerates to the previous models. Our model can well explain the significant re-brightening in the optical and infrared light curves of GRB 081029 and GRB 100621A. A considerable fall-back mass is needed to provide the later energy, this indicates GRBs accompanied with fall-back material may be associated with a low energy supernova so that fraction of the envelope can be survived during eruption. The fall-back time can give meaningful information of the properties of GRB progenitor stars.

SDSS J1634+2049 is a local (z = 0.1293) infrared-luminous quasar with LIR= 10^11.91 Lsun. We present a detailed multiwavelength study of both the host galaxy and the nucleus. The host galaxy demonstrates violent, obscured star formation activities with SFR ~ 140 Msun yr^-1, estimated from either the PAH emission or IR luminosity. The optical to NIR spectra exhibit a blueshifted narrow cuspy component in Hb, HeI5876,10830 and other emission lines consistently with an offset velocity of ~900 km/s, as well as additional blueshifting phenomena in high-ionization lines , while there exist blueshifted broad absorption lines (BALs) in NaID and HeI*3889,10830, indicative of the AGN outflows producing BALs and emission lines. Constrained mutually by the several BALs with CLOUDY, the physical properties of the absorption-line outflow are derived as follows: 10^4 < n_H <= 10^5 cm^-3, 10^-1.3 <= U <= 10^-0.7 and 10^22.5<= N_H <= 10^22.9 cm^-2 , similar to those derived for the emission-line outflows. The similarity suggests a common origin. Taking advantages of both the absorption lines and outflowing emission lines, we find that the outflow gas is located at a distance of 48 - 65 pc from the nucleus, and that the kinetic luminosity of the outflow is 10^44-10^46 erg s^-1. J1634+2049 has a off-centered galactic ring on the scale of ~ 30 kpc that is proved to be formed by a recent head-on collision by a nearby galaxy. Thus this quasar is a valuable object in the transitional phase emerging out of dust enshrouding as depicted by the co-evolution scenario.

In this work we propose a new diagnostic to segregate cool core (CC) clusters from non-cool core (NCC) clusters by studying the two-dimensional power spectra of the X-ray images observed with the Chandra X-ray observatory. Our sample contains 41 members ($z=0.01\sim 0.54$), which are selected from the Chandra archive when a high photon count, an adequate angular resolution, a relatively complete detector coverage, and coincident CC-NCC classifications derived with three traditional diagnostics are simultaneously guaranteed. We find that in the log-log space the derived image power spectra can be well represented by a constant model component at large wavenumbers, while at small wavenumbers a power excess beyond the constant component appears in all clusters, with a clear tendency that the excess is stronger in CC clusters. By introducing a new CC diagnostic parameter, i.e., the power excess index (PEI), we classify the clusters in our sample and compare the results with those obtained with three traditional CC diagnostics. We find that the results agree with each other very well. By calculating the PEI values of the simulated clusters, we find that the new diagnostic works well at redshifts up to 0.5 for intermediately sized and massive clusters with a typical Chandra or XMM pointing observation. The new CC diagnostic has several advantages over its counterparts, e.g., it is free of the effects of the commonly seen centroid shift of the X-ray halo caused by merger event, and the corresponding calculation is straightforward, almost irrelevant to the complicated spectral analysis.

In this paper we propose that a fast radio burst (FRB) could originate from the magnetic interaction between double neutron stars (NSs) during their final inspiral within the framework of a unipolar inductor model. In this model, an electromotive force is induced on one NS to accelerate electrons to an ultra-relativistic speed instantaneously. We show that coherent curvature radiation from these electrons moving along magnetic field lines in the magnetosphere of the other NS is responsible for the observed FRB signal, that is, the characteristic emission frequency, luminosity, duration and event rate of FRBs can be well understood. In addition, we discuss several implications of this model, including double-peaked FRBs and possible associations of FRBs with short-duration gamma-ray bursts and gravitational wave events.

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.

There is growing evidence that a stable magnetar could be formed from the coalescence of double neutron stars. In previous papers, we investigated the signature of formation of stable millisecond magnetars in radio and optical/ultraviolet bands by assuming that the central rapidly rotating magnetar deposits its rotational energy in the form of a relativistic leptonized wind. We found that the optical transient PTF11agg could be the first evidence for the formation of post-merger millisecond magnetars. To enhance the probability of finding more evidence for the post-merger magnetar formation, it is better to extend the observational channel to other photon energy bands. In this paper we propose to search the signature of post-merger magnetar formation in X-ray and especially gamma-ray bands. We calculate the SSC emission of the reverse shock powered by post-merger millisecond magnetars. We find that the SSC component peaks at $1\,{\rm GeV}$ in the spectral energy distribution and extends to $\gtrsim 10\,{\rm TeV}$ for typical parameters. These energy bands are quite suitable for Fermi/LAT and CTA, which, with their current observational sensitivities, can detect the SSC emission powered by post-merger magnetars up to $1\,{\rm Gpc}$. NuSTAR, sensible in X-ray bands, can detect the formation of post-merger millisecond magnetars at redshift $z\sim 1$. Future improvement in sensitivity of CTA can also probe the birth of post-merger millisecond magnetars at redshift $z\sim 1$. However, because of the $\gamma$-$\gamma$ collisions, strong high-energy emission is clearly predicted only for ejecta masses lower than $10^{-3}M_\odot$.

High-redshift quasars are important tracers of structure and evolution in the early universe. However, they are very rare and difficult to find when using color selection because of contamination from late-type dwarfs. High-redshift quasar surveys based on only optical colors suffer from incompleteness and low identification efficiency, especially at $z\gtrsim4.5$. We have developed a new method to select $4.7\lesssim z \lesssim 5.4$ quasars with both high efficiency and completeness by combining optical and mid-IR Wide-field Infrared Survey Explorer (WISE) photometric data, and are conducting a luminous $z\sim5$ quasar survey in the whole Sloan Digital Sky Survey (SDSS) footprint. We have spectroscopically observed 99 out of 110 candidates with $z$-band magnitudes brighter than 19.5 and 64 (64.6\%) of them are quasars with redshifts of $4.4\lesssim z \lesssim 5.5$ and absolute magnitudes of $-29\lesssim M_{1450} \lesssim -26.4$. In addition, we also observed 14 fainter candidates selected with the same criteria and identified 8 (57.1\%) of them as quasars with $4.7<z<5.4$ . Among 72 newly identified quasars, 12 of them are at $5.2 < z < 5.7$, which leads to an increase of $\sim$36\% of the number of known quasars at this redshift range. More importantly, our identifications doubled the number of quasars with $M_{1450}<-27.5$ at $z>4.5$, which will set strong constraints on the bright end of the quasar luminosity function. We also expand our method to select quasars at $z\gtrsim5.7$. In this paper we report the discovery of four new luminous $z\gtrsim5.7$ quasars based on SDSS-WISE selection.

Detection and mitigation of radio frequency interference (RFI) is the first and also the key step for data processing in radio observations, especially for ongoing low frequency radio experiments towards the detection of the cosmic dawn and epoch of reionization (EoR). In this paper we demonstrate the technique and efficiency of RFI identification and mitigation for the 21 Centimeter Array (21CMA), a radio interferometer dedicated to the statistical measurement of EoR. For terrestrial, man-made RFI, we concentrate mainly on a statistical approach by identifying and then excising non-Gaussian signatures, in the sense that the extremely weak cosmic signal is actually buried under thermal and therefore Gaussian noise. We also introduce the so-called visibility correlation coefficient instead of conventional visibility, which allows a further suppression of rapidly time-varying RFI. Finally, we briefly discuss removals of the sky RFI, the leakage of sidelobes from off-field strong radio sources with time-invariant power and a featureless spectrum. It turns out that state of the art technique should allow us to detect and mitigate RFI to a satisfactory level in present low frequency interferometer observations such as those acquired with the 21CMA, and the accuracy and efficiency can be greatly improved with the employment of low-cost, high-speed computing facilities for data acquisition and processing.

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.

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.

The central engine and jet composition of gamma-ray bursts (GRBs) remain mysterious. Here we suggest that observations on polarization evolution of early optical afterglows may shed light on these questions. We first study the dynamics of a reverse shock and a forward shock that are generated during the interaction of a relativistic jet and its ambient medium. The jet is likely magnetized with a globally large-scale magnetic field from the central engine. The existence of the reverse shock requires that the magnetization degree of the jet should not be high ($\sigma\leq 1$), so that the jet is mainly composed of baryons and leptons. We then calculate the light curve and polarization evolution of an early optical afterglow, and find that when the polarization position angle changes by $90^\circ$ during the early afterglow, the polarization degree is zero for a toroidal magnetic field but is very likely to be non-zero for an aligned magnetic field. This result would be expected to provide a probe for the central engine of GRBs, because an aligned field configuration could originate from a magnetar central engine and a toroidal field configuration could be produced from a black hole via the Blandford-Znajek mechanism. Finally, for such two kinds of magnetic field configurations, we fit the observed data of the early optical afterglow of GRB 120308A equally well.

A small fraction($<10\%$) of SDSS main sample galaxies(MGs) have not been targeted with spectroscopy due to the the fiber collision effect. These galaxies have been compiled into the input catalog of the LAMOST extra-galactic survey and named as the complementary galaxy sample. In this paper, we introduce the project and the status of the spectroscopies of the complementary galaxies in the first two years of the LAMOST spectral survey(till Sep. of 2014). Moreover, we present a sample of 1,102 galaxy pairs identified from the LAMOST complementary galaxies and SDSS MGs, which are defined as that the two members have a projected distance smaller than 100 kpc and the recessional velocity difference smaller than 500 $\rm kms^{-1}$. Compared with the SDSS only selected galaxy pairs, the LAMOST-SDSS pairs take the advantages of not being biased toward large separations and therefor play as a useful supplement to the statistical studies of galaxy interaction and galaxy merging.

The ultraviolet broad absorption lines have been seen in the spectra of quasars at high redshift, and are generally considered to be caused by outflows with velocities from thousands kilometers per second to one tenth of the speed of light. They provide crucial implications for the cosmological structures and physical evolutions related to the feedback of active galactic nuclei (AGNs). Recently, through a dedicated program of optically spectroscopic identifications of selected quasar candidates at redshift 5 by using the Lijiang 2.4 m telescope, we discovered two luminous broad absorption line quasars (BALQSOs) at redshift about 4.75. One of them may even have the potentially highest absorption Balnicity Index (BI) ever found to date, which is remarkably characterized by its deep, broad absorption lines and sub-relativistic outflows. Further physical properties, including the metal abundances, variabilities, evolutions of the supermassive black holes (SMBH) and accretion disks associated with the feedback process, can be investigated with multi-wavelength follow-up observations in the future.

The particle mass used in cosmology N-body simulations is close to $10^{8}M_{\odot}$, which is about $10^{65}$ times larger than the GeV scale expected in particle physics. However, self-gravity interacting particle systems made up of different particles mass have different statistical and dynamical properties. Here we demonstrate that, due to this particle mass difference, the nowaday cosmology N-body simulations can have introduced an excessive core collapse process, especially for the low mass halos at high redshift. Such dynamical effect introduces an excessive cuspy center for these small halos, and it implies a possible connection to the so called "small scale crisis" for CDM models. Our results show that there exist a physical limit in cosmological simulations, by using about $10^3$ particles to describe smallest halos, and we provide a simple suggestion based on it to relieve those effects from the bias.

In order to investigate the spatial distribution of the ICM temperature in galaxy clusters in a quantitative way and probe the physics behind, we analyze the X-ray spectra of a sample of 50 galaxy clusters, which were observed with the Chandra ACIS instrument in the past 15 years, and measure the radial temperature profiles out to $0.45r_{500}$. We construct a physical model that takes into account the effects of gravitational heating, thermal history (such as radiative cooling, AGN feedback, and thermal conduction) and work done via gas compression, and use it to fit the observed temperature profiles by running Bayesian regressions. The results show that in all cases our model provides an acceptable fit at the 68% confidence level. To further validate this model we select nine clusters that have been observed with both Chandra (out to $\gtrsim 0.3r_{500}$) and Suzaku (out to $\gtrsim 1.5r_{500}$), fit their Chandra spectra with our model, and compare the extrapolation of the best-fits with the Suzaku measurements. We find that the model profiles agree with the Suzaku results very well in seven clusters. In the rest two clusters the difference between the model and observation is possibly caused by local thermal substructures. Our study also implies that for most of the clusters the assumption of hydrostatic equilibrium is safe out to at least $0.5r_{500}$, and the non-gravitational interactions between dark matter and its luminous counterpart is consistent with zero.

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 present preliminary results of the quasar survey in Large Sky Area Multi- Object Fiber Spectroscopic Telescope (LAMOST) first data release (DR1), which includes pilot survey and the first year regular survey. There are 3921 quasars identified with reliability, among which 1180 are new quasars discovered in the survey. These quasars are at low to median redshifts, with highest z of 4.83. We compile emission line measurements around the H{\alpha}, H{\beta}, Mg II, and C IV regions for the new quasars. The continuum luminosities are inferred from SDSS photo- metric data with model fitting as the spectra in DR1 are non-flux-calibrated. We also compile the virial black hole mass estimates, and flags indicating the selec- tion methods, broad absorption line quasars. The catalog and spectra for these quasars are available online. 28% of the 3921 quasars are selected with optical- infrared colours independently, indicating that the method is quite promising in completeness of quasar survey. LAMOST DR1 and the on-going quasar survey will provide valuable data in the studies of quasars.