In a multi-wavelength survey of $13$ quasars at $5.8\lesssim z\lesssim6.5$, that were pre-selected to be potentially young, we find five objects with extremely small proximity zone sizes that may imply UV-luminous quasar lifetimes of $\lesssim 100,000$ years. Proximity zones are regions of enhanced transmitted flux in the vicinity of the quasars that are sensitive to the quasars' lifetimes because the intergalactic gas has a finite response time to their radiation. We combine sub-mm observations from the Atacama Large Millimetre Array (ALMA) and the NOrthern Extended Millimeter Array (NOEMA), as well as deep optical and near-infrared spectra from medium-resolution spectrograph on the Very Large Telescope (VLT) and on the Keck telescopes, in order to identify and characterize these new young quasars, which provide valuable clues about the accretion behavior of supermassive black holes (SMBHs) in the early universe, and pose challenges on current black hole formation models to explain the rapid formation of billion solar mass black holes. We measure the quasars' systemic redshifts, black hole masses, Eddington ratios, emission line luminosities, and star formation rates of their host galaxies. Combined with previous results we estimate the fraction of young objects within the high-redshift quasar population at large to be $5\%\lesssim f_{\rm young}\lesssim 10\%$. One of the young objects, PSO J158-14, shows a very bright dust continuum flux ($F_{\rm cont}=3.46\pm 0.02\,\rm mJy$), indicating a highly star-bursting host galaxy with a star formation rate of approximately $1420\,M_{\odot}\,\rm yr^{-1}$.

Standard siren cosmology of gravitational wave (GW) merger events relies on the identification of host galaxies and their redshifts. But this can be highly challenging due to numerous candidates of galaxies in the GW localization area. We point out that the number of candidates can be reduced by orders of magnitude for strongly lensed GW events, due to extra observational constraints. For the next-generation GW detectors like Einstein Telescope (ET), we estimate that this number is usually significantly less than one, as long as the GW localization uncertainty is better than $\sim 10\, \rm deg^2$. This implies that the unique identification of the host galaxy of lensed GW event detected by ET and Cosmic Explorer (CE) is possible. This provides us a promising opportunity to measure the redshift of the GW event and facilitate the standard siren cosmology. We also discuss its potential applications in understanding the evolution process and environment of the GW event.

We present deep (9 hours) Gemini-N/GNIRS near-infrared spectroscopic observations of ULAS J1342+0928, a luminous quasar at z=7.54. Various broad emission lines were detected, as well as the underlying continuum and iron forests over the rest-frame wavelength 970-2930A. There is a clear trend that higher-ionization emission lines show larger blueshifts with CIV1549 exhibiting 5510^{+240}_{-110} km s-1 blueshift with respect to the systematic redshift from the far-infrared [CII] 158um emission line. Those high ionization lines have wide profiles with FWHM more than 10000 km s-1. A modest blueshift of 340^{+110}_{-80} km s-1 is also seen in MgII, the lowest ionization line identified in the spectrum. The updated MgII-based black hole mass of M_BH=9.1_{-1.3}^{+1.4} x 10^8 M_sun and the Eddington ratio of L_bol/L_Edd=1.1_{-0.2}^{+0.2} confirm that ULAS J1342+0928 is powered by a massive and actively accreting black hole. There is no significant difference in the emission line ratios such as SiIV/CIV and AlIII/CIV when compared to lower-redshift quasars in a similar luminosity range, which suggests early metal pollution of the broad-line-region clouds. This trend also holds for the FeII/MgII line ratio, known as a cosmic clock that traces the iron enrichment in the early universe. Different iron templates and continuum fitting ranges were used to explore how the FeII/MgII measurement changes as a function of spectral modeling. Quasars at even higher redshift or at fainter luminosity range (L_bol<10^46 erg s-1) are needed to probe the sites of early metal enrichment and a corresponding change in the FeII/MgII ratio.

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.

From near-infrared spectroscopic measurements of the MgII emission line doublet, we estimate the black hole (BH) mass of the quasar, SMSS J215728.21-360215.1, as being (3.4 +/- 0.6) x 10^10 M_sun and refine the redshift of the quasar to be z=4.692. SMSS J2157 is the most luminous known quasar, with a 3000A luminosity of (4.7 +/- 0.5) x 10^47 erg/s and an estimated bolometric luminosity of 1.6 x 10^48 erg/s, yet its Eddington ratio is only ~0.4. Thus, the high luminosity of this quasar is a consequence of its extremely large BH -- one of the most massive BHs at z > 4.

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.

PypeIt is a Python package for semi-automated reduction of astronomical, spectroscopic data. Its algorithms build on decades-long development of previous data reduction pipelines by the developers (Bernstein, Burles, & Prochaska, 2015; Bochanski et al., 2009). The reduction procedure -- including a complete list of the input parameters and available functionality -- is provided as online documentation hosted by Read the Docs, which is regularly updated. (https://pypeit.readthedocs.io/en/latest/). Release v1.0.3 serves the following spectrographs: Gemini/GNIRS, Gemini/GMOS, Gemini/FLAMINGOS 2, Lick/Kast, Magellan/MagE, Magellan/Fire, MDM/OSMOS, Keck/DEIMOS (600ZD, 830G, 1200G), Keck/LRIS, Keck/MOSFIRE (J and Y gratings tested), Keck/NIRES, Keck/NIRSPEC (low-dispersion), LBT/Luci-I, Luci-II, LBT/MODS (beta), NOT/ALFOSC (grism4), VLT/X-Shooter (VIS, NIR), VLT/FORS2 (300I, 300V), WHT/ISIS.

We reveal that the energy spectra of electrons-positrons in primary cosmic rays measured at atmosphere top have double structures: an excess component $\Phi^s_{e^+}(E)=\Phi^s_{e^-}(E)$ around $400 GeV$, which origins from a strong $e^+e^-$-source and the distorted background $\Phi^0_{e^-}(E)$. We supposed that the difference between AMS-CALET and Fermi-LAT-DAMPE data origins from the energy loss of the fluxes due to the anomalous bremsstrahlung effect at a special window. The evolution of spectra under anomalous bremsstrahlung effect satisfies an improved electromagnetic cascade equation. The above spectra are parameterized and they can be regarded as the subjects exploring new physics. We suggest to check the previous applications of the Bethe-Heitler formula in the study of the propagation of high energy electrons and photons.

The primary task of the 1.26-m telescope jointly operated by the National Astronomical Observatory and Guangzhou University is photometric observations of the g, r, and i bands. A data processing pipeline system was set up with mature software packages, such as IRAF, SExtractor, and SCAMP, to process approximately 5 GB of observational data automatically every day. However, the success ratio was significantly reduced when processing blurred images owing to telescope tracking error; this, in turn, significantly constrained the output of the telescope. We propose a robust automated photometric pipeline (RAPP) software that can correctly process blurred images. Two key techniques are presented in detail: blurred star enhancement and robust image matching. A series of tests proved that RAPP not only achieves a photometric success ratio and precision comparable to those of IRAF but also significantly reduces the data processing load and improves the efficiency.

The Hubble parameter $H(z)$ is directly related to the expansion of our Universe. It can be used to study dark energy and constrain cosmology models. In this paper, we propose that $H(z)$ can be measured using fast radio bursts (FRBs) with redshift measurements. We use dispersion measures contributed by the intergalactic medium, which is related to $H(z)$, to measure Hubble parameter. We find that 500 mocked FRBs with dispersion measures and redshift information can accurately measure Hubble parameters using Monte Carlo simulation. The maximum deviation of $H(z)$ from standard $\Lambda$CDM model is about 6\% at redshift $z= 2.4$. We also test our method using Monte Carlo simulation. Kolmogorov-Smirnov (K-S) test is used to check the simulation. The $p$-value of K-S test is 0.23, which confirms internal consistency of the simulation. In future, more localizations of FRBs make it as an attractive cosmological probe.

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.

Significant progress in the classification of Fermi unassociated sources , has led to an increasing number of blazars are being found. The optical spectrum is effectively used to classify the blazars into two groups such as BL Lacs and flat spectrum radio quasars (FSRQs). However, the accurate classification of the blazars without optical spectrum information, i.e., blazars of uncertain type (BCUs), remains a significant challenge. In this paper, we present a principal component analysis (PCA) and machine learning hybrid blazars classification method. The method, based on the data from Fermi LAT 3FGL Catalog, first used the PCA to extract the primary features of the BCUs and then used a machine learning algorithm to further classify the BCUs. Experimental results indicate that the that the use of PCA algorithms significantly improved the classification. More importantly, comparison with the Fermi LAT 4FGL Catalog, which contains the spectral classification of those BCUs in the Fermi-LAT 3FGL Catalog, reveals that the proposed classification method in the study exhibits higher accuracy than currently established methods; specifically, 151 out of 171 BL Lacs and 19 out of 24 FSRQs are correctly classified.

Quasi-biennial oscillations (QBOs) are considered as a fundamental mode of solar magnetic activity at low latitudes ($\leq50^\circ$). However, the evolutionary aspect and the hemispheric distribution of solar QBOs at high latitudes ($\geq60^\circ$) are rarely studied. Here, a relatively novel time-frequency analysis technique, named the synchrosqueezed wavelet transform, is applied to extract the main components of the polar faculae in the northern and southern hemispheres for the time interval from August 1951 to December 1998. It is found as the following: (1) Apart from the 22-year Hale cycle, the 17-year extended activity cycle, and the 11-year Schwabe cycle, the QBOs have been estimated as a prominent timescale of solar magnetic activity at high latitudes; (2) the QBOs of the polar faculae are coherent in the two hemispheres, but the temporal (phase) and the spatial (amplitude) variations of solar QBOs occur unevenly on both hemispheres; and (3) for the 11-year period mode, the northern hemisphere begins three months earlier than that in the southern one. Moreover, the spatial and temporal distributions of the hemispheric QBOs differ from those of the 11-year Schwabe cycle mode in the two hemispheres. Our findings could be helpful to improve our knowledge on the physical origin of the spatial distribution of solar QBOs at high latitudes, and could also provide more constraints on solar dynamo models introduced to characterize the different components of the solar magnetic activity cycle.

The nuclide $^{14}$C can be produced in the atmosphere by high energy particles and $\gamma$-rays from high-energy phenomena. Through the carbon cycle, some of $^{14}$CO$_2$ produced in the atmosphere can be retained in annual tree rings. Four events of rapid increase of the $^{14}$C content occurred in AD 775, AD 994, BC 660 and BC 3371 were found. Recently, the data of Jull et al. (2020) was inconsistent with our records around BC 3371. We measured our sample again and found the $^{14}$C records are consistent with the value in Wang et al. (2017). Therefore, our $^{14}$C records are robust. The inconsistency may be caused by the difference of calendar ages for the wood samples, or the physical origin of the event. First, crossdating on ring width can be performed only between trees whose growth has the same environmental conditions. Because the master tree-ring for dendrochronology is lack for Chinese trees. The master tree-ring from California has to be used. Therefore, the calendar ages derived from dendrochronology may be not precise. Second, the $^{14}$C even may be not global. One evidence is the variation of $^{14}$C content around AD 1006. The $^{14}$C contents of Californian trees increase 12\textperthousand~ in two years, while Japanese trees show no $^{14}$C increase.

In the third catalog of active galactic nuclei detected by the $Fermi$ Large Area Telescope Clean (3LAC) sample, there are 402 blazars candidates of uncertain type (BCU). The proposed analysis will help to evaluate the potential optical classification flat spectrum radio quasars (FSRQs) versus BL Lacertae (BL Lacs) objects of BCUs, which can help to understand which is the most elusive class of blazar hidden in the Fermi sample. By studying the 3LAC sample, we found some critical values of $\gamma$-ray photon spectral index ($\Gamma_{\rm ph}$), variability index (VI) and radio flux (${\rm F_R}$) of the sources separate known FSRQs and BL Lac objects. We further utilize those values to defined an empirical "high-confidence" candidate zone that can be used to classify the BCUs. Within such a zone ($\Gamma_{\rm ph}<2.187$, log${\rm F_R}<2.258$ and ${ \rm logVI <1.702}$), we found that 120 BCUs can be classified BL Lac candidates with a higher degree of confidence (with a misjudged rate $<1\%$). Our results suggest that an empirical "high confidence" diagnosis is possible to distinguish the BL Lacs from the Fermi observations based on only on the direct observational data of $\Gamma_{\rm ph}$, VI and ${\rm F_R}$.

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.

Fast radio bursts (FRBs) are extragalactic, bright pulses of emission at radio frequency with milliseconds duration. Observationally, FRBs can be divided into two classes, repeating FRBs and non-repeating FRBs. At present, twenty repeating FRBs have been discovered with unknown physical origins. Localization of the first repeating FRB 121102 and discovery of an associated persistent radio source support that FRBs are powered by young millisecond magnetars, which could be formed by core-collapses of massive stars or binary neutron stars mergers. These two formation channels can be distinguished by gravitational waves generated by binary neutron stars mergers. We first calculate the lower limit of the local formation rate of repeating FRBs observed by the Canadian Hydrogen Intensity Mapping Experiment (CHIME) . Then we show that only a small fraction ($6\%$) of repeating FRBs is produced by young magnetars from binary neutron star mergers, basing on the gravitational wave detections by the third observing run (O3) of Advanced LIGO/Virgo gravitational-wave detectors. Therefore, we believe that repeating FRBs are more likely produced by the magnetars newborn from the core-collapses of massive stars rather than the magnetars from the binary neutron stars mergers.

Hardening of cosmic proton- and nuclei-spectra is explained by using the gluon condensation (GC) model, which states that a large amount of gluons in proton may condense near the high energy threshold. The results present the universality of the GC-effect, which may be a common origin of a series of anomalous astrophysical phenomena including the broken power-law in gamma-ray spectra, the excess in positron- and electron-spectra and hardening of proton- and nuclei-spectra.

Young neutron stars (NSs) born in core-collapse explosions are promising candidates for the central engines of fast radio bursts (FRBs), since the first localized repeating burst FRB 121102 happens in a star forming dwarf galaxy, which is similar to the host galaxies of superluminous supernovae (SLSNe) and long gamma-ray bursts (LGRBs). However, FRB 180924 and FRB 190523 are localized to massive galaxies with low rates of star formation, compared with the host of FRB 121102. Meanwhile, the offsets between the bursts and host centers are about 4 kpc and 29 kpc for FRB 180924 and FRB 190523, respectively. These properties of hosts are similar to short gamma-ray bursts \textbf{(SGRBs)}, which are produced by mergers of binary neutron star (BNS) or neutron star-black hole (NS-BH). Therefore, the NSs powering FRBs may be formed in BNS mergers. In this paper, we study the BNS merger rates, merger times, and predict their most likely merger locations for different types of host galaxies using population synthesis method. We find that the BNS merger channel is consistent with the recently reported offsets of FRB 180924 and FRB 190523. The offset distribution of short GRBs is well reproduced by population synthesis using galaxy model which is similar to GRB hosts. The event rate of FRBs (including non-repeating and repeating), is larger than those of BNS merger and short GRBs, which requires a large fraction of observed FRBs emitting several bursts. Using curvature radiation by bunches in NS magnetospheres, we also predict the observational properties of FRBs from BNS mergers, including the dispersion measure, and rotation measure. At late times ($t\geq1$yr), the contribution to dispersion measure and rotation measure from BNS merger ejecta could be neglected.

We present a catalog of 5,290 RR Lyrae stars (RRLs) with metallicities estimated from spectra of the LAMOST Experiment for Galactic Understanding and Exploration (LEGUE) and the Sloan Extension for Galactic Understanding and Exploration (SEGUE) surveys. Nearly 70 per cent of them (3,642 objects) also have systemic radial velocities measured. Given the pulsating nature of RRLs, metallicity estimates are based on spectra of individual exposures, by matching them with the synthetic templates. The systemic radial velocities are measured by fitting the observed velocity as a function of phase assuming an empirical pulsating velocity template curve. Various tests show that our analyses yield metallicities with a typical precision of 0.20\,dex and systemic radial velocities with uncertainties ranging from 5 to 21\,km\,s$^{-1}$ (depending on the number of radial velocity measurements available for a given star). Based on the well calibrated near-infrared $PM_{W1}Z$ or $PM_{K_{\rm s}}Z$, and $M_{V}$-[Fe/H] relations, precise distances are derived for these RRLs. Finally, we include Gaia DR2 proper motions in our catalog. The catalog should be very useful for various Galactic studies, especially of the Galactic halo.

Superflares, as strong explosions on stars, have been well studied with the progress of space time-domain astronomy. In this work, we present the study of superflares on solar-type stars using Transiting Exoplanet Survey Satellite ({\em{TESS}}) data. 13 sectors of observations during the first year of the {\em TESS} mission have covered the southern hemisphere of the sky, containing 25,734 solar-type stars. We verified 1,216 superflares on 400 solar-type stars through automatic search and visual inspection with 2-minute cadence data. Our result suggests a higher superflare frequency distribution than the result from {\em Kepler}. The reason may be that the majority of {\em TESS} solar-type stars in our dataset are rapidly rotating stars. The power-law index $\gamma$ of the superflare frequency distribution ($dN/dE\propto E^{-\gamma}$) is constrained to be $\gamma = 2.16\pm 0.10$, which is a little larger than that of solar flares but consistent with the results from {\em Kepler}. Because only 7 superflares of Sun-like stars are detected, we may not give a robust superflare occurrence frequency. And four stars are accompanied by unconfirmed hot planet candidates. Therefore, superflares are possibly caused by stellar magnetic activities instead of planet-star interactions. We also find an extraordinary star TIC43472154, which exhibits about 200 superflares per year. In addition, the correlation between energy and duration of superflares ($T_{\text {duration }} \propto E^{\beta}$) is analyzed. We derive the power-law index to be $\beta=0.42\pm0.01$, which is a little larger than $\beta=1/3$ from the prediction according to magnetic reconnection theory.

The gluon condensation effect in cosmic proton-proton collisions at high energy is used to explain an excess in the positron spectrum observed by AMS. We find that this excess may originate from the GC-effect in Tycho's supernova remnant.

Temporal variation of the solar coronal rotation appears to be very complex and its relevances to the eleven-year solar activity cycle are still unclear. Using the modified coronal index for the time interval from 1939 January 1 to 2019 May 31, the systematic regularities of the solar coronal rotation are investigated. Our main findings are as follows: (1) from a global point of view, the synodic coronal rotation period with a value of 27.5 days is the only significant period at the periodic scales shorter than 64 days; (2) the coronal rotation period exhibit an obviously decreasing trend during the considered time interval, implying the solar corona accelerates its global rotation rate in the long run; (3) there exist significant periods of 3.25, 6.13, 9.53, and 11.13 years in the period length of the coronal rotation, providing an evidence that the coronal rotation should be connected with the quasi-biennial oscillation, the eleven-year solar cycle, and the 22-year Hale cycle (or the magnetic activity reversal); and (4) the phase relationship between the coronal rotation period and the solar magnetic activity is not only time-dependent but also frequency-dependent. For a small range around the 11- year cycle band, there is a systematic trend in the phase, and the small mismatch in this band brings out the phase to drift. The possible mechanism for the above analysis results is discussed.

The discovery of quasars few hundred megayears after the Big Bang represents a major challenge to our understanding of black holes and galaxy formation and evolution. Their luminosity is produced by extreme gas accretion onto black holes, which already reached masses of 10$^9$ M$_\odot$ by z ~ 6. Simultaneously, their host galaxies form hundreds of stars per year, using up gas in the process. To understand which environments are able to sustain the rapid formation of these extreme sources we started a VLT/MUSE effort aimed at characterizing the surroundings of a sample of 5.7 < z < 6.6 quasars dubbed: the Reionization Epoch QUasar InvEstigation with MUSE (REQUIEM) survey. We here present results of our searches for extended Ly-Alpha halos around the first 31 targets observed as part of this program. Reaching 5-sigma surface brightness limits of 0.1-1.1 x 10$^{-17}$ erg/s/cm$^2$/arcsec$^2$ over a 1 arcsec$^2$ aperture, we were able to unveil the presence of 12 Ly-Alpha nebulae, 8 of which are newly discovered. The detected nebulae show a variety of emission properties and morphologies with luminosities ranging from 8 x 10$^{42}$ to 2 x 10$^{44}$ erg/s, FWHMs between 300 and 1700 km/s, sizes < 30 pkpc, and redshifts consistent with those of the quasar host galaxies. As the first statistical and homogeneous investigation of the circum-galactic medium of massive galaxies at the end of the reionization epoch, the REQUIEM survey enables the study of the evolution of the cool gas surrounding quasars in the first 3 Gyr of the Universe. A comparison with the extended Ly-Alpha emission observed around bright (M$_{1450}$ < -25 mag) quasars at intermediate redshift indicates little variations on the properties of the cool gas from z ~ 6 to z ~ 3 followed by a decline in the average surface brightness down to z ~ 2.

The extreme mass ratio inspiral (EMRI), defined as a stellar-mass compact object inspiraling into a supermassive black hole (SMBH), has been widely argued to be a low-frequency gravitational wave (GW) source. EMRIs providing accurate measurements of black hole mass and spin, are one of the primary interests for Laser Interferometer Space Antenna (LISA). However, it is usually believed that there are no electromagnetic (EM) counterparts to EMRIs. Here we show a new formation channel of EMRIs with tidal disruption flares as EM counterparts. In this scenario, flares can be produced from the tidal stripping of the helium (He) envelope of a massive star by an SMBH. The left compact core of the massive star will evolve into an EMRI. We find that, under certain initial eccentricity and semimajor axis, the GW frequency of the inspiral can enter LISA band within 10 $\sim$ 20 years, which makes the tidal disruption flare an EM precursor to EMRI. Although the event rate is just $2\times 10^{-4}~\rm Gpc^{-3}yr^{-1}$, this association can not only improve the localization accuracy of LISA and help to find the host galaxy of EMRI, but also serve as a new GW standard siren for cosmology.

Four significant events of rapid $^{14}$C increase have taken place within the past several thousand years. The physical origin of these rapid increases is still a mystery but must be associated with extremely energetic cosmic processes. Pulsars are highly magnetized neutron stars that emit a beam of electromagnetic radiations. Any sudden release of the energy stored in the magnetic multipole field will trigger outbursts similar to the giant flares of magnetars. Here we show that the relativistic outflow from the outbursts of a nearby pulsar interacting with the interstellar medium generates a shock, which accelerates electrons to trillions of electron volts. The high-energy photons from synchrotron emission of the shock interact with Earth's atmosphere, producing the cosmogenic nuclide $^{14}$C, which can cause the rapid $^{14}$C increases discovered in tree rings. These same relativistic electrons can account for a significant fraction of the cosmic electron spectrum in the trillion electron volts energy range, as observed by space-borne satellites. Since these outburst events can significantly affect our environment, monitoring nearby pulsars for such outbursts may be important in the future.

In this paper, we present statistics of soft gamma repeater (SGR) bursts from SGR J1550-5418, SGR 1806-20 and SGR 1900+14 by adding new bursts from K{\i}rm{\i}z{\i}bayrak et al. (2017) detected with the Rossi X-ray Timing Explorer (RXTE). We find that the fluence distributions of magnetar bursts are well described by power-law functions with indices 1.84, 1.68, and 1.65 for SGR J1550-5418, SGR 1806-20 and SGR 1900+14, respectively. The duration distributions of magnetar bursts also show power-law forms. Meanwhile, the waiting time distribution can be described by a non-stationary Poisson process with an exponentially growing occurrence rate. These distributive features indicate that magnetar bursts can be regarded as a self-organizing critical process. We also compare these distributions with the repeating FRB 121102. The statistical properties of repeating FRB 121102 are similar with magentar bursts, combing with the large required magnetic filed ($B\geq 10^{14}$G) of neutron star for FRB 121102, which indicates that the central engine of FRB 121102 may be a magnetar.

Quasi-biennial oscillation (QBO) of solar magnetic activities is intrinsic to dynamo mechanism, but still far from fully understood. In this work, the phase and amplitude asymmetry of solar QBO of Halpha flare activity in the northern and southern hemispheres is studied by the ensemble empirical mode decomposition, the cross-correlation analysis, and the wavelet transform technique. The following results are found: (1) solar QBO of Halpha flare index in the two hemispheres has a complicated phase relationship, but does not show any systematic regularity; (2) the solar cycle mode of solar Halpha flare index in the northern hemisphere generally leads that in the southern one by 9 months for the time interval from 1966 to 2014. The possible origin of these results is discussed.

Studies on the periodic variation and the phase relationship between different solar activity indicators are useful for understanding the long-term evolution of solar activity cycle. Here we report the statistical analysis of grouped solar flare (GSF) and sunspot number (SN) during the time interval from January 1965 to March 2009. We find that, 1) the significant periodicities of both GSF and SN are related to the differential rotation periodicity, the quasi-biennial oscillation (QBO), and the eleven-year Schwabe cycle (ESC), but the specific values are not absolutely identical; 2) the ESC signal of GSF lags behind that of SN with an average of 7.8 months during the considered time interval, implying that the systematic phase delays between GSF and SN originate from the inter-solar-cycle signal. Our results may provide evidence about the storage of magnetic energy in the corona.

Nuclear Coulomb potential at completely ionized and extremely thin atmosphere can leaks to a macroscopic spatial scale. This effect is used to explain the difference between the energy spectra of cosmic-ray electrons around $1 TeV$ measured by different experimental groups. The result inspires us to review the traditional electromagnetic shower theory at the extreme conditions. It is also reminds us that the energy spectrum of cosmic-ray electrons, which are measured by Fermi-LAT and DAMPE at a higher altitude is more closer to a true signature.

Recent studies have indicated that an anisotropic cosmic expansion may exist. In this paper, we use three datasets of type Ia supernovae (SNe Ia) to probe the isotropy of cosmic acceleration. For the Union2.1 dataset, the direction and magnitude of the dipole are $(l=309.3^{\circ} {}^{+ 15.5^{\circ}}_{-15.7^{\circ}} ,\ b = -8.9^{\circ} {}^{ + 11.2^{\circ}}_{-9.8^{\circ}} ),\ A=(1.46 \pm 0.56) \times 10^{-3}$. For the Constitution dataset, the results are $(l=67.0^{\circ}{}^{+ 66.5^{\circ}}_{-66.2^{\circ}},\ b=-0.6^{\circ}{}^{+ 25.2^{\circ}}_{-26.3^{\circ}}),\ A=(4.4 \pm 5.0) \times 10^{-4}$. For the JLA dataset, no significant dipolar deviation is found. We also explore the effects of anisotropic distributions of coordinates and redshifts on the results using Monte-Carlo simulations. We find that the anisotropic distribution of coordinates can cause dipole directions and make dipole magnitude larger. Anisotropic distribution of redshifts is found to have no significant effect on dipole fitting results.

We present ALMA 0.28"x0.20" (1.4kpc x 1.0kpc) resolution observations of the [CII] 158um line and underlying dust continuum in the host galaxy of the most distant quasar currently known, ULAS J1342+0928 at z=7.5413. Both the [CII] and continuum are detected and spatially resolved. The [CII] line is ~1.5 times more extended than the continuum emission, showing an elongated and complex structure of approximately 3.2kpc x 6.4kpc. Two separate peaks are clearly seen (>6sigma each) in three 100 km/s width [CII] channel maps. The velocity field of the [CII] gas does not show evidence of a coherent rotation field but rather chaotic motion, reminiscent of an on-going merger.

We present the photometric observations of blazars S5 0716+714 and 3C 273 with high temporal resolution (30--60s) in the $I$ or $R$ bands. The observations were performed with a 1.02 m optical telescope from 2007 March 07 to 2012 May 16. The $F$-test, one-way analysis of variance (ANOVA) test, and $z$-transformed discrete correlation function (ZDCF) cross-correlation zero lag test are used to search for intra-day variability (IDV). Four and five reliable IDVs survive three tests for S5 0716+714 and 3C 273, respectively. IDVs are found for S5 0716+714 and 3C 273. A flare on 2008 May 08 has $\Delta I \approx$ 0.06$\pm$0.01 mag in a duration of 0.54 hr for S5 0716+714. A flare on 2011 May 10 shows $\Delta R \approx$ 0.05$\pm$0.01 mag in a duration of 0.40 hr for 3C 273. Sharp dips appear on 2011 May 9 for 3C 273, and show $\Delta R \approx$ 0.05$\pm$0.01 mag. Under the assumptions that the IDV is tightly connected to black hole mass, $M_{\bullet}$, and that the flare durations are representative of the minimum characteristic timescales, we can estimate upper bounds to $M_{\bullet}$. In the case of the Kerr black holes, $M_{\bullet} \la 10^{8.91} M_{\odot}$ and $M_{\bullet} \la 10^{9.02} M_{\odot}$ are given for S5 0716+714 and 3C 273, respectively. These mass measurements are consistent with those measurements reported in the literatures. Also, we discuss the origins of optical variations found in this work.

Fast radio bursts (FRBs) are milliseconds radio transients with large dispersion measures (DMs). An outstanding question is the relation between repeating FRBs and those with a single burst. In this paper, we study the energy distribution of the repeating FRB 121102. For a power-law distribution of energy $dN/dE\propto E^{-\alpha_E}$, we show that the value of $\alpha_E$ is in a narrow range $1.6-1.8$ for bursts observed by different telescopes at different frequencies, which indicates a universal energy distribution for FRB 121102. Interestingly, similar power-law index of energy distribution for non-repeating FRBs observed by Parkes and ASKAP is also found. However, if low-energy bursts below completeness threshold of Arecibo are discarded for FRB 121102, the slope could be up to 2.2. Implications of such a universal energy distribution are discussed.

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.

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 use ALMA observations of the host galaxy of the quasar ULAS-J1342+0928 at z=7.54 to study the dust continuum and far infrared lines emitted from its interstellar medium. The Rayleigh-Jeans tail of the dust continuum is well sampled with eight different spectral setups, and from a modified black body fit we obtain an emissivity coefficient of beta=1.85+-0.3. Assuming a standard dust temperature of 47 K we derive a dust mass of Mdust=0.35x10^8 Msol and a star formation rate of 150+-30 Msol/yr. We have >4sigma detections of the [CII]_158, [OIII]_88 and [NII]_122 atomic fine structure lines and limits on the [CI]_369, [OI]_146 and [NII]_205 emission. We also report multiple limits of CO rotational lines with Jup >= 7, as well as a tentative 3.3sigma detection of the stack of four CO lines (Jup=11, 10, 8, 7). We find line deficits that are in agreement with local ultra luminous infrared galaxies. Comparison of the [NII]_122 and [CII]_158 lines indicates that the [CII]_158 emission arises predominantly from the neutral medium, and we estimate that the photo-disassociation regions in J1342+0928 have densities <~5x10^4cm^-3. The data suggest that ~16% of hydrogen is in ionized form and that the HII regions have high electron densities of n_e>180 cm^-3. Our observations favor a low gas-to-dust ratio of <100 and a metallicity of the interstellar medium comparable to the Solar value. All the measurements presented here suggest that the host galaxy of J1342+0928 is highly enriched in metal and dust, despite being observed just 680 Myr after the Big Bang.

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.

In this paper, we investigate the energy function, formation rate and environment of fast radio bursts (FRBs) using Parkes sample and Australian Square Kilometer Array Pathfinder (ASKAP) sample. For the first time, the metallicity effect on the formation rate is considered. If FRBs are produced by the mergers of compact binaries, the formation rate of FRBs should have a time delay relative to cosmic star formation rate (CSFR). We get the time delay is about 3-5 Gyr and the index of differential energy function $\gamma$ ($dN/dE\propto E^{-\gamma}$) is between 1.6 and 2.0 from redshift cumulative distribution. The value of $\gamma$ is similar to that of FRB 121102, which indicates single bursts may share the same physical mechanism with the repeaters. In another case, if the formation rate of FRB is proportional to the SFR without time delay, the index $\gamma$ is about 2.3. In both cases, we find that FRBs may prefer to occur in low-metallicity environment with $ 12 +\log(\rm{O/H}) \simeq 8.40$, which is similar to those of long gamma-ray bursts (GRBs) and hydrogen-poor superluminous supernovae (SLSNe-I).

Fast radio bursts (FRBs) are bright milliseconds radio transients with large dispersion measures (DMs), and can be used as potential astrophysical and cosmological tools. Many models for FRB progenitors have been proposed. However, none of them can explain all the observational properties, meaning that FRBs remain one of the most intriguing mysteries in astronomy. The statistical properties of radio bursts can unveil the underlying physics. Here, we report statistical results of the repeating FRB 121102, and show that FRB 121102 and solar type III radio bursts share four statistical properties: power-law frequency distributions for energies, fluxes, durations and waiting times. All of the distributions can be explained by avalanche models of self-organized criticality (SOC) systems. It is well known that solar type III radio burst arises from the nonlinear conversion of Langmuir waves generated by two-stream instability of electron beams, which are accelerated by bursty magnetic reconnections. The similarities support that repeating FRBs are coherent emissions from magnetic-reconnection-driven beams in the magnetospheres of magnetars.

The discovery of luminous quasars at redshifts up to 7.5 demonstrates the existence of several billion M_sun supermassive black holes (SMBHs) less than a billion years after the Big Bang. They are accompanied by intense star formation in their host galaxies, pinpointing sites of massive galaxy assembly in the early universe, while their absorption spectra reveal an increasing neutral intergalactic medium (IGM) at the epoch of reionization. Extrapolating from the rapid evolution of the quasar density at z=5-7, we expect that there is only one luminous quasar powered by a billion M_sun SMBH in the entire observable universe at z~9. In the next decade, new wide-field, deep near-infrared (NIR) sky surveys will push the redshift frontier to the first luminous quasars at z~9-10; the combination with new deep X-ray surveys will probe fainter quasar populations that trace earlier phases of SMBH growth. The identification of these record-breaking quasars, and the measurements of their BH masses and accretion properties require sensitive spectroscopic observations with next generation of ground-based and space telescopes at NIR wavelengths. High-resolution integral-field spectroscopy at NIR, and observations at millimeter and radio wavelengths, will together provide a panchromatic view of the quasar host galaxies and their galactic environment at cosmic dawn, connecting SMBH growth with the rise of the earliest massive galaxies. Systematic surveys and multiwavelength follow-up observations of the earliest luminous quasars will strongly constrain the seeding and growth of the first SMBHs in the universe, and provide the best lines of sight to study the history of reionization.

We collected the information of 304 gamma-ray bursts (GRBs) from the literature, and analyzed the correlations among the host galaxy offsets (the distance from the site of the GRB to the center of its host galaxy), $T_{\rm 90,i}$ (the duration $T_{\rm 90}$ in rest-frame), $T_{\rm R45,i}$ (the duration $T_{\rm R45}$ in rest-frame), $E_{\rm \gamma,iso}$ (the isotropic equivalent energy), $L_{\rm \gamma,iso}$ ($=E_{\rm \gamma,iso}/T_{\rm 90,i}$, the isotropic equivalent luminosity) and $L_{\rm pk}$ (peak luminosity). We found that $ T_{\rm 90,i}$, $ T_{\rm R45,i}$, $ E_{\rm \gamma,iso}$, $L_{\rm pk}$ have negative correlation with $ {\rm offset}$, which is consistent with origin of short GRBs (SGRBs) and long GRBs (LGRBs). On separate analysis, we found similar results for $\log E_{\rm \gamma,iso}$ - $\log {\rm (offset)}$ and $\log L_{\rm pk}$ - $\log {\rm (offset)}$ relations in case of SGRBs only, while no obvious relation for LGRBs. There is no correlations between offset and $L_{\rm \gamma,iso}$. We also put the special GRB 170817A {and GRB 060218A} on the plots. {The two GRBs both have low luminosity and small offset.} In the $ \log ({\rm offset})- \log T_{\rm 90,i}$ plot, we found GRB 170817A locates in between the two regions of SGRBs and LGRBs and it is the outlier in the $ {\rm offset}-E_{\rm \gamma, iso}$, $ {\rm offset}-L_{\rm \gamma, iso}$ and $ {\rm offset}-L_{\rm pk}$ plots. Together with GRB 060218A being an outlier in all plots, it indicates the speciality of GRBs 170817A and 060218A, and might imply more subgroups of the GRB samples.

Gamma-ray bursts (GRBs) are a potential tool to probe high-redshift universe. However, the circularity problem enforces people to find model-independent methods to study the luminosity correlations of GRBs. Here, we present a new method which uses gravitational waves as standard sirens to calibrate GRB luminosity correlations. For the third-generation ground-based GW detectors (i.e., Einstein Telescope), the redshifts of gravitational wave (GW) events accompanied electromagnetic counterparts can reach out to $\sim 4$, which is more distant than type Ia supernovae ($z\lesssim 2$). The Amati relation and Ghirlanda relation are calibrated using mock GW catalogue from Einstein Telescope. We find that the $1\sigma$ uncertainty of intercepts and slopes of these correlations can be constrained to less than 0.2\% and 8\% respectively. Using calibrated correlations, the evolution of dark energy equation of state can be tightly measured, which is important for discriminating dark energy models.

In order to obtain an overview of the gamma-ray bursts (GRBs), we need a full sample. In this paper, we collected 6289 GRBs (from GRB 910421 to GRB 160509A) from the literature, including prompt emission, afterglow and host galaxy properties. We hope to use this large sample to reveal the intrinsic properties of GRB. We have listed all the data in machine readable tables, including the properties of the GRBs, correlation coefficients and linear regression results of two arbitrary parameters, and linear regression results of any three parameters. These machine readable tables could be used as a data reservoir for further studies on the classifications or correlations. One may find some intrinsic properties from these statistical results. With this comprehensive table, it is possible to find relations between different parameters, and to classify the GRBs into different kinds of sub-groups. With the completion, it may reveal the nature of GRBs and may be used as tools like pseudo-redshift indicators, standard candles, etc. All the machine readable data and statistical results are available on the website of the journal.

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.

Recently, some divergent conclusions about cosmic acceleration were obtained using type Ia supernovae (SNe Ia), with opposite assumptions on the intrinsic luminosity evolution. In this paper, we use strong gravitational lensing systems to probe the cosmic acceleration. Since the theory of strong gravitational lensing is established certainly, and the Einstein radius is determined by stable cosmic geometry. We study two cosmological models, $\Lambda$CDM and power-law models, through 152 strong gravitational lensing systems, incorporating with 30 Hubble parameters $H(z)$ and 11 baryon acoustic oscillation (BAO) measurements. Bayesian evidence are introduced to make a one-on-one comparison between cosmological models. Basing on Bayes factors $\ln B$ of flat $\Lambda$CDM versus power-law and $R_{h}=ct$ models are $\ln B>5$, we find that the flat $\Lambda$CDM is strongly supported by the combination of the datasets. Namely, an accelerating cosmology with non power-law expansion is preferred by our numeration.

We have designed the Extremely Luminous Quasar Survey (ELQS) to provide a highly complete census of unobscured UV-bright quasars during the cosmic noon, $z=2.8-5.0$. Here we report the discovery of 70 new quasars in the ELQS South Galactic Cap (ELQS-S) quasar sample, doubling the number of known extremely luminous quasars in $4,237.3\,\rm{deg}^2$ of the SDSS footprint. These observations conclude the ELQS and we present the properties of the full ELQS quasar catalog, containing 407 quasars over $11,838.5\,\rm{deg}^2$. Our novel ELQS quasar selection strategy resulted in unprecedented completeness at the bright end and allowed us to discover 109 new quasars in total. This marks an increase of $\sim36\%$ (109/298) to the known population at these redshifts and magnitudes, while we further are able to retain a selection efficiency of $\sim80\%$. On the basis of 166 quasars from the full ELQS quasar catalog, who adhere to the uniform criteria of the 2MASS point source catalog, we measure the bright-end quasar luminosity function (QLF) and extend it one magnitude brighter than previous studies. Assuming a single power law with exponential density evolution for the functional form of the QLF, we retrieve the best fit parameters from a maximum likelihood analysis. We find a steep bright-end slope of $\beta\approx-4.1$ and we can constrain the bright-end slope to $\beta\leq-3.4$ with $99\%$ confidence. The density is well modeled by the exponential redshift evolution, resulting in a moderate decrease with redshift ($\gamma\approx-0.4$).

In this paper, we study the correlation between isotropic energy and duration of gamma-ray bursts (GRBs) for the first time. The correlation is found to be $T_d \propto {E_{iso}}^{0.34\pm 0.03}$ from the {\em Swift} GRB sample. After comparing with solar flares from {\em RHESSI} and stellar superflares from {\em Kepler} satellite, we find that the correlation of GRBs shows similar exponent with those of solar flares and stellar superflares. Inspired by the physical mechanism of solar flares and stellar superflares which is magnetic reconnection, we interpret the correlation using magnetic reconnection theory. This similarity hints that magnetic reconnection may dominate energy releasing process of GRBs.

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.

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.

We report the discovery of a large ($\sim$8500 km diameter) infrared-bright storm at Neptune's equator in June 2017. We tracked the storm over a period of 7 months with high-cadence infrared snapshot imaging, carried out on 14 nights at the 10 meter Keck II telescope and 17 nights at the Shane 120 inch reflector at Lick Observatory. The cloud feature was larger and more persistent than any equatorial clouds seen before on Neptune, remaining intermittently active from at least 10 June to 31 December 2017. Our Keck and Lick observations were augmented by very high-cadence images from the amateur community, which permitted the determination of accurate drift rates for the cloud feature. Its zonal drift speed was variable from 10 June to at least 25 July, but remained a constant $237.4 \pm 0.2$ m s$^{-1}$ from 30 September until at least 15 November. The pressure of the cloud top was determined from radiative transfer calculations to be 0.3-0.6 bar; this value remained constant over the course of the observations. Multiple cloud break-up events, in which a bright cloud band wrapped around Neptune's equator, were observed over the course of our observations. No "dark spot" vortices were seen near the equator in HST imaging on 6 and 7 October. The size and pressure of the storm are consistent with moist convection or a planetary-scale wave as the energy source of convective upwelling, but more modeling is required to determine the driver of this equatorial disturbance as well as the triggers for and dynamics of the observed cloud break-up events.

Using the Hilbert-Huang transform technique, we investigate the midrange periodicities in solar radio flux at 2800 MHz (F10.7) and sunspot areas (SAs) from February 1, 1947 to September 30, 2016. The following prominent results are found: (1) The quasi-periodic oscillations of both data sets are not identical, such as the rotational cycle, the midrange periodicities, and the Schwabe cycle. In particular, the midrange periodicities ranging from 37.9 days to 297.3 days are related to the magnetic Rossby-type waves; 2) The 1.3-year and 1.7-year fluctuations in solar activity indicators are surface manifestations (from photosphere to corona) of magnetic flux changes generated deep inside the Sun; 3) At the timescale of the Schwabe cycle, \textbf{the complicated phase relationships} in the three intervals (1947-1958, 1959-1988, and 1989-2016) agree with the produced periodicities of the magnetic Rossby-type waves. \textbf{The findings indicate that the magnetic Rossby-type waves are the possible physical mechanism behind the midrange periodicities of solar activity indicators. Moreover, the significant change in the relationship between photospheric and coronal activity took place after the maximum of solar cycle 22 could be interpreted by the magnetic Rossby-type waves.

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.

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.

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.

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

A new type of gamma ray spectrum is predicted in a general hadronic framework if considering the gluon condensation effects in proton. The result presents the power-law with a sharp break in the gamma ray spectra at the TeV-band. We suggest to probe this GC-signature in the Earth limb gamma-ray spectra using the DArk Matter Particle Explorer and the CALorimetric Electron Telescope on orbit.

Gluons dominate the proton behavior at high energy collisions, they can be condensed at ultra high energy. The collisions of the accelerated high energy protons with interplanetary matter in cosmic rays will produce a huge number of secondary particles at the gluon condensate energy region, which break the primary power-law of cosmic rays. The above predictions seem to be consistent with the recent DAMPE data concerning the electron plus positron spectra. We find that the smoothly broken power-law at $\sim 0.9 TeV$ and $3\sim 4 TeV$ in the DAMPE data can be understood as the gluon condensation effects in proton.

Surveying dark matter deficient galaxies (those with dark matter mass to stellar mass ratio $M_{\rm dm}/M_{\rm star}<1$) in the Illustris simulation of structure formation in the flat-$\Lambda$CDM cosmogony, we find $M_{\rm star} \approx 2 \times 10^8\, M_\sun$ galaxies that have properties similar to those ascribed by \citet{vanDokkumetal2018a} to the ultra-diffuse galaxy NGC1052-DF2. The Illustris simulation also contains more luminous dark matter deficient galaxies. Illustris galaxy subhalo 476171 is a particularly interesting outlier, a massive and very compact galaxy with $M_{\rm star} \approx 9 \times 10^{10}\, M_\sun$ and $M_{\rm dm}/M_{\rm star} \approx 0.1$ and a half-stellar-mass radius of $\approx 2$ kpc. If the Illustris simulation and the $\Lambda$CDM model are accurate, there are a significant number of dark matter deficient galaxies, including massive luminous compact ones. It will be interesting to observationally discover these galaxies, and to also more clearly understand how they formed, as they are likely to provide new insight into and constraints on models of structure formation and the nature of dark matter.

Peculiar velocities are a precious tool to study the large-scale distribution of matter in the local universe and test cosmological models. However, present measurements of peculiar velocities are based on empirical distance indicators, which introduce large error bars. Here we present a new method to measure the peculiar velocities, by directly estimating luminosity distances through waveform signals from inspiralling compact binaries and measuring redshifts from electromagnetic (EM) counterparts. In the future, with the distance uncertainty of GW events reducing to $0.1$ per cent by future GW detectors, the uncertainty of the peculiar velocity can be reduced to $10$ km/s at 100 mega parsecs. We find that dozens of GW events with EM counterparts can provide a Hubble constant $H_0$ uncertainty of $0.5\%$ and the growth rate of structure with a $0.6\%$ precision in the third-generation ground-base GW detectors, which can reconcile the $H_0$ tension and determine the origins for cosmic accelerated expansion.

Current constraints on Einstein's weak equivalence principle (WEP) utilize the observed time delay between correlated particles of astronomical sources. However, the intrinsic time delay due to particle emission time is impossible to measure and is simply assumed to be zero in previous studies. Here we propose a robust method to test WEP for relativistic messengers using strongly lensed cosmic transients, which can naturally overcome the intrinsic time delay problem. This can be achieved in two ways: by comparing the time delays between lensed images seen in different energy bands or in gravitational waves (GWs) and their electromagnetic (EM) counterparts. The power of our method mainly depends on the timing accuracy of cosmic transient and the strong lensing time delay. If the time delay of cosmic transient can be measured with accuracy about 0.1 s (e.g. gamma-ray bursts), we show that the upper limit on the differences of the parameterized post-Newtonian parameter $\gamma$ value is $\Delta \gamma<10^{-7}$ with a one-month strong lensing time delay event. This accuracy of WEP can be improved by a factor of $10^{7}$, if the lens is galaxy cluster and the strongly lensed cosmic transients have much shorter duration, such as fast radio bursts.

We present the North Galactic Cap sample of the Extremely Luminous Quasar Survey (ELQS-N), which targets quasars with $M_{1450}<-27$ at $2.8 \leq z < 5$ in an area of $\sim7600\,\rm{deg}^2$ of the Sloan Digital Sky Survey (SDSS) footprint with $90\text{\textdegree}<\rm{RA}<270\text{\textdegree}$. Based on a near-infrared/infrared \textit{JKW2} color cut, the ELQS selection efficiently uses random forest methods to classify quasars and to estimate photometric redshifts; this scheme overcomes some of the difficulties of pure optical quasar selection at $z\approx3$. As a result, we retain a completeness of $>70\%$ over $z\sim3.0-5.0$ at $m_{i}\lesssim17.5$, limited toward fainter magnitudes by the depth of the Two Micron All Sky Survey (2MASS). The presented quasar catalog consists of a total of 270 objects, of which 39 are newly identified in this work with spectroscopy obtained at the Vatican Advanced Technology Telescope and the MMT $6.5\,\rm{m}$ telescope. In addition to the high completeness, which allowed us to discover new quasars in the already well-surveyed SDSS North Galactic Cap, the efficiency of our selection is relatively high at $\sim79\%$. Using 120 objects of this quasar sample we are able to extend the previously measured optical quasar luminosity function (QLF) by one magnitude toward the bright end at $2.8 \leq z \leq 4.5$. A first analysis of the QLF suggests a relatively steep bright-end slope of $\beta\approx-4$ for this sample. This result contrasts with previous results in the same redshift range, which find a much flatter slope around $\beta\sim-2.5$, but agrees with recent measurements of the bright-end slope at lower and higher redshifts. Our results constrain the bright-end slope at $z=2.8-4.5$ to $\beta<-2.94$ with a 99\% confidence.

As one of the series of papers reporting on a large reverberation mapping campaign, we apply the maximum entropy method (MEM) to 9 narrow-line Seyfert 1 galaxies with super-Eddington accretion rates observed during 2012-2013 for the velocity-delay maps of their H\beta\ and H\gamma\ emission lines. The maps of 6 objects are reliably reconstructed using MEM. The maps of H\beta\ and H\gamma\ emission lines of Mrk 335 indicate that the gas of its broad-line region (BLR) is infalling. For Mrk 142, its H\beta\ and H\gamma\ lines show signatures of outflow. The H\beta\ and H\gamma\ maps of Mrk 1044 demonstrate complex kinematics -- a virialized motion accompanied by an outflow signature, and the H\beta\ map of IRAS F12397+3333 is consistent with a disk or a spherical shell. The H\beta\ maps of Mrk 486 and MCG +06-26-012 suggest the presence of an inflow and outflow, respectively. These super-Eddington accretors show diverse geometry and kinematics. Brief discussions of their BLRs are provided for each individual object.

The spectral flux density of stars can indicate their atmospheric physical properties. A detector can obtain any band flux density at the design stage. However, the band flux density is confirmed and fixed in the process of operation because of the restriction of filters. Other band flux densities cannot be obtained through the same detector. In this study, a computational model of stellar spectral flux density is established based on basic physical parameters which are effective temperature and angular parameter. The stochastic particle swarm optimization algorithm is adopted to address this issue with appropriately chosen values of the algorithm parameters. Four star catalogues are studied and consist of the Large Sky Area Multi-Object Fibre Spectroscopic Telescope (LAMOST), Wide-field Infrared Survey Explorer (WISE), Midcourse Space Experiment (MSX), and Two Micron All Sky Survey (2MASS). The given flux densities from catalogues are input parameters. Stellar effective temperatures and angular parameters are inversed using the given flux densities according to SPSO algorithm. Then the flux density is calculated according to Planck's law on the basis of stellar effective temperatures and angular parameters. The calculated flux density is compared with the given value from catalogues. It is found that the inversion results are in good agreement for all bands of the MSX and 2MASS catalogues, whereas they do not agree well in some bands of the LAMOST and WISE catalogues. Based on the results, data from the MSX and 2MASS catalogues can be used to calculate the spectral flux density at different wavelengths of given wavelength ranges. The stellar flux density is obtained and can provide data support and an effective reference for detection and recognition of stars.

Quasars are the most luminous non-transient objects known and as a result they enable studies of the Universe at the earliest cosmic epochs. Despite extensive efforts, however, the quasar ULAS J1120+0641 at z=7.09 has remained the only one known at z>7 for more than half a decade. Here we report observations of the quasar ULAS J134208.10+092838.61 (hereafter J1342+0928) at redshift z=7.54. This quasar has a bolometric luminosity of 4e13 times the luminosity of the Sun and a black hole mass of 8e8 solar masses. The existence of this supermassive black hole when the Universe was only 690 million years old---just five percent of its current age---reinforces models of early black-hole growth that allow black holes with initial masses of more than about 1e4 solar masses or episodic hyper-Eddington accretion. We see strong evidence of absorption of the spectrum of the quasar redwards of the Lyman alpha emission line (the Gunn-Peterson damping wing), as would be expected if a significant amount (more than 10 per cent) of the hydrogen in the intergalactic medium surrounding J1342+0928 is neutral. We derive a significant fraction of neutral hydrogen, although the exact fraction depends on the modelling. However, even in our most conservative analysis we find a fraction of more than 0.33 (0.11) at 68 per cent (95 per cent) probability, indicating that we are probing well within the reionization epoch of the Universe.

In this paper, we study the anisotropy of cosmic acceleration the using Pantheon sample, which includes 1048 spectroscopically confirmed Type Ia supernovae (SNe Ia) covering the redshift range $0.01 < z < 2.3$. In hemisphere comparison method, we find the dipole direction is $(l = 37 \pm 40^{\circ}, b = 33 \pm 16^{\circ})$ with the maximum anisotropy level of $\delta=0.136 {}^{+0.009}_{-0.005}$. From the dipole fitting method, we find that the magnitude of anisotropy is $A = (3.7 {}^{+2.5}_{-3.7}) \times 10^{-4}$, and the direction of the dipole $(l = 329^{\circ}{}^{+ 101^{\circ}}_{-28^{\circ}}, b = 37^{\circ}{}^{+ 52^{\circ}}_{-21^{\circ}})$ in the galactic coordinate system. The result is weakly dependent on redshift from the redshift tomography analysis. The anisotropy is small and the isotropic cosmological model is an excellent approximation.

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.

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.

In this paper we report the experimental implementation of a theoretically-proposed technique for creating a photoionized plasma in the laboratory using X-ray line radiation. Using a Sn laser-plasma to irradiate an Ar gas target, the photoionization parameter, \xi = 4{\pi}F/Ne, reached values of order 50 erg cm/s, where F is the radiation flux in erg/cm2/s. The significance of this is that this technique allows us to mimic effective spectral radiation temperatures in excess of 1 keV. We show that our plasma starts to be collisionally dominated before the peak of the X-ray drive. However, the technique is extendable to higher energy laser systems to create plasmas with parameters relevant to benchmarking codes used to model astrophysical objects.

The Einstein's weak equivalence principle (WEP) is one of the foundational assumptions of general relativity and some other gravity theories. In the theory of parametrized post-Newtonian (PPN), the difference between the PPN parameters $\gamma$ of different particles or the same type of particle with different energies, $\Delta \gamma$, represents the violation of WEP. Current constraints on $\Delta \gamma$ are derived from the observed time delay between correlated particles of astronomical sources. However, the observed time delay is contaminated by other effects, such as the time delays due to different particle emission times, the potential Lorentz invariance violation and none-zero photon rest mass. Therefore, current constraints are only upper limits. Here we propose a new method to test WEP based on the fact that the gravitational time delay is direction-dependent while others are not. This is the first method which can naturally correct other time delay effects. Using the time delay measurements of BASTE gamma-ray burst (GRB) sample and the gravitational potential of local super galaxy cluster Laniakea, we find that the constraint on $\Delta \gamma$ of different energy photons can be as low as $10^{-14}$. In the future, if more gravitational wave events and fast radio bursts with much more precise time delay measurements are observed, this method can give a reliable and tight constraint on WEP.

In this paper, we test the cosmic distance duality (CDD) relation using the luminosity distances from joint light-curve analysis (JLA) type Ia supernovae (SNe Ia) sample and angular diameter distance sample from galaxy clusters. The $R_h=ct$ and $\Lambda$CDM models are considered. In order to compare the two models, we constrain the CCD relation and the SNe Ia light-curve parameters simultaneously. Considering the effects of Hubble constant, we find that $\eta\equiv D_A(1+z)^2/D_L=1$ is valid at the 2$\sigma$ confidence level in both models with $H_0 = 67.8\pm0.9$ km/s/Mpc. However, the CDD relation is valid at 3$\sigma$ confidence level with $H_0 = 73.45\pm1.66$ km/s/Mpc. Using the Akaike Information Criterion (AIC) and the Bayesian Information Criterion (BIC), we find that the $\Lambda$CDM model is very strongly preferred over the $R_h=ct$ model with these data sets for the CDD relation test.

This paper presents a fast algorithm for obtaining high-accuracy subpixel translation of low PSNR images. Instead of locating the maximum point on the upsampled images or fitting the peak of correlation surface, the proposed algorithm is based on the measurement of centroid on the cross correlation surface by Modified Moment method. Synthetic images, real solar images and standard testing images with white Gaussian noise added were tested, and the results show that the accuracies of our algorithm are comparable with other subpixel registration techniques and the processing speed is higher. The drawback is also discussed at the end of this paper.

We use high spatial and temporal resolution observations, simultaneously obtained with the New Vacuum Solar Telescope and Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory, to investigate the high-frequency oscillations above a sunspot umbra. A novel time--frequency analysis method, namely the synchrosqueezing transform (SST), is employed to represent their power spectra and to reconstruct the high-frequency signals at different solar atmospheric layers. A validation study with synthetic signals demonstrates that SST is capable to resolving weak signals even when their strength is comparable with the high-frequency noise. The power spectra, obtained from both SST and the Fourier transform, of the entire umbral region indicate that there are significant enhancements between 10 and 14 mHz (labeled as 12 mHz) at different atmospheric layers. Analyzing the spectrum of a photospheric region far away from the umbra demonstrates that this 12~mHz component exists only inside the umbra. The animation based on the reconstructed 12 mHz component in AIA 171 \AA\ illustrates that an intermittently propagating wave first emerges near the footpoints of coronal fan structures, and then propagates outward along the structures. A time--distance diagram, coupled with a subsonic wave speed ($\sim$ 49 km s$^{-1}$), highlights the fact that these coronal perturbations are best described as upwardly propagating magnetoacoustic slow waves. Thus, we first reveal the high-frequency oscillations with a period around one minute in imaging observations at different height above an umbra, and these oscillations seem to be related to the umbral perturbations in the photosphere.

A novel time--frequency technique, called the synchrosqueezing transform (SST), is used to investigate the midterm periodic variations of magnetic fields on the solar surface. The Magnetic Plage Strength Index (MPSI) and the Mount Wilson Sunspot Index (MWSI), measured daily by the Mount Wilson Observatory between 1970 January 19 and 2012 January 22, are selected. The short-, mid, and longer-term periodicities are represented and decomposed by the SST with hardly any mode mixing. This demonstrates that the SST is a useful time--frequency analysis technique to characterize the periodic modes of helioseismic data. Apart from the fundamental modes of the annual periodicity, $\sim$27 day rotational cycle and $\sim$11 year solar cycle, the SST reveals several midterm periodicities in the two magnetic activity indices, specifically, $\sim$157 day (i.e., Rieger-type periodicity), and $\sim$1.3 and 1.7 years. The periodic modes, with 116.4 and 276.2 day periodicity in the MPSI, with 108.5 and 251.6 day periodicity in the MWSI, and the 157.7 day periodicity in the two indices, are in better accord with those significant periodicities derived from the Rossby waves theoretical model. This study suggests that the modes are caused by the Rossby waves. For the 1.30 and 1.71 year periodicity of the MPSI, and the 1.33 and 1.67 year periodicity of the MWSI, our analysis infers that they are related to those periodicity with the same timescale in the interior of the Sun and in the high atmospheric layers.

We present a 45 ks Chandra observation of the quasar ULAS J1342+0928 at z=7.54. We detect 14.0^{+4.8}_{-3.7} counts from the quasar in the observed-frame energy range 0.5-7.0 keV (6-sigma detection), representing the most distant non-transient astronomical source identified in X-rays to date. The present data are sufficient only to infer rough constraints on the spectral parameters. We find an X-ray hardness ratio of HR = -0.51^{+0.26}_{-0.28} between the 0.5-2.0 keV and 2.0-7.0 keV ranges and derive a power-law photon index of Gamma = 1.95^{+0.55}_{-0.53}. Assuming a typical value for high-redshift quasars of Gamma = 1.9, ULAS J1342+0928 has a 2-10 keV rest-frame X-ray luminosity of L_{2-10} = 11.6^{+4.3}_{-3.5} x 10^{44} erg/s. Its X-ray-to-optical power-law slope is alpha_{OX}=-1.67^{+0.16}_{-0.10}, consistent with the general trend indicating that the X-ray emission in the most bolometrically powerful quasars is weaker relative to their optical emission.

Left-Right symmetric model (LRSM) has been an attractive extension of the Standard Model (SM) which can address the origin of parity violation in the SM electroweak (EW) interactions, generate tiny neutrino masses, accommodate dark matter (DM) candidates and provide a natural framework for baryogenesis through leptogenesis. In this work we utilize the minimal LRSM to study the recently reported DAMPE results of cosmic $e^+e^-$ spectrum which exhibits a tentative peak around 1.4 TeV, while satisfying the current neutrino data. We propose to explain the DAMPE peak with a complex scalar DM $\chi$ in two scenarios: 1) $\chi\chi^* \to H_1^{++}H_1^{--} \to \ell_i^+\ell_i^+\ell_j^-\ell_j^-$; 2) $\chi\chi^* \to H_{k}^{++}H_{k}^{--} \to \ell_i^+\ell_i^+\ell_j^-\ell_j^-$ accompanied by $\chi\chi^* \to H_1^+ H_1^- \to \ell_i^+ \nu_{\ell_i} \ell_j^- \nu_{\ell_j}$ with $\ell_{i,j}=e,\mu,\tau$ and $k=1,2$. We fit the theoretical prediction on $e^+e^-$ spectrum to relevant experimental data to determine the scalar mass spectrum favored by the DAMPE excess. We also consider various constraints from theoretical principles, collider experiments as well as DM relic density and direct search experiments. We find that there are ample parameter space which can interpret the DAMPE data while passing the constraints. Our explanations, on the other hand, usually imply the existence of other new physics at the energy scale ranging from $10^7 {\rm GeV}$ to $10^{11} {\rm GeV}$. Collider tests of our explanations are also discussed.

We compile a complete collection of reliable Hubble parameter $H(z)$ data to redshift $z \leq 2.36$ and use them with the Gaussian Process method to determine continuous $H(z)$ functions for various data subsets. From these continuous $H(z)$'s, summarizing across the data subsets considered, we find $H_0\sim 67 \pm 4\,\rm km/s/Mpc$, more consistent with the recent lower values determined using a variety of techniques. In most data subsets, we see a cosmological deceleration-acceleration transition at 2$\sigma$ significance, with the data subsets transition redshifts varying over $0.33<z_{\rm da}<1.0$ at 1$\sigma$ significance. We find that the flat-$\Lambda$CDM model is consistent with the $H(z)$ data to a $z$ of 1.5 to 2.0, depending on data subset considered, with 2$\sigma$ deviations from flat-$\Lambda$CDM above this redshift range. Using the continuous $H(z)$ with baryon acoustic oscillation distance-redshift observations, we constrain the current spatial curvature density parameter to be $\Omega_{K0}=-0.03\pm0.21$, consistent with a flat universe, but the large error bar does not rule out small values of spatial curvature that are now under debate.

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.

During reionization, neutral hydrogen in the intergalactic medium (IGM) imprints a damping wing absorption feature on the spectrum of high-redshift quasars. A detection of this signature provides compelling evidence for a significantly neutral Universe, and enables measurements of the hydrogen neutral fraction $x_{\rm HI}(z)$ at that epoch. Obtaining reliable quantitative constraints from this technique, however, is challenging due to stochasticity induced by the patchy inside-out topology of reionization, degeneracies with quasar lifetime, and the unknown unabsorbed quasar spectrum close to rest-frame Ly$\alpha$. We combine a large-volume semi-numerical simulation of reionization topology with 1D radiative transfer through high-resolution hydrodynamical simulations of the high-redshift Universe to construct models of quasar transmission spectra during reionization. Our state-of-the-art approach captures the distribution of damping wing strengths in biased quasar halos that should have reionized earlier, as well as the erosion of neutral gas in the quasar environment caused by its own ionizing radiation. Combining this detailed model with our new technique for predicting the quasar continuum and its associated uncertainty, we introduce a Bayesian statistical method to jointly constrain the neutral fraction of the Universe and the quasar lifetime from individual quasar spectra. We apply this methodology to the spectra of the two highest redshift quasars known, ULAS J1120+0641 and ULAS J1342+0928, and measured volume-averaged neutral fractions $\langle x_{\rm HI} \rangle(z=7.09)=0.48^{+0.26}_{-0.26}$ and $\langle x_{\rm HI} \rangle(z=7.54)=0.60^{+0.20}_{-0.23}$ (posterior medians and 68% credible intervals) when marginalized over quasar lifetimes of $10^3 \leq t_{\rm q} \leq 10^8$ years.

Using a simplified framework, we attempt to explain the recent DAMPE cosmic $e^+ + e^-$ flux excess by leptophilic Dirac fermion dark matter (LDM). The scalar ($\Phi_0$) and vector ($\Phi_1$) mediator fields connecting LDM and Standard Model particles are discussed. Under constraints of DM relic density, gamma-rays, cosmic-rays and Cosmic Microwave Background (CMB), we find that the couplings $P \otimes S$, $P \otimes P$, $V \otimes A$ and $V \otimes V$ can produce the right bump in $e^+ + e^-$ flux for a DM mass around 1.5 TeV with a natural thermal annihilation cross-section $<\sigma v> \sim 3 \times 10^{-26} cm^3/s$ today. Among them, $V \otimes V$ coupling is tightly constrained by PandaX-II data (although LDM-nucleus scattering appears at one-loop level) and the surviving samples appear in the resonant region, $m_{\Phi_1} \simeq 2m_{\chi}$. We also study the related collider signatures, such as dilepton production $pp \to \Phi_1 \to \ell^+\ell^-$, and muon $g-2$ anomaly. Finally, we present a possible $U(1)_X$ realization for such leptophilic dark matter.

The possibility that a fraction of the dark matter is comprised of massive compact halo objects (MACHOs) remains unclear, especially in the 20-100 $M_{\odot}$ window. MACHOs could make up binaries, whose mergers may be detected by LIGO as gravitational wave events. On the other hand, the cosmological origin of fast radio burst (FRBs) has been confirmed. We investigate the possibility of detecting fast radio bursts (FRBs) gravitational lensed by MACHO binaries to constrain their properties. Since lensing events could generate more than one images, lensing by binaries could cause multiple-peak FRBs. The angular separation between these images is roughly $10^{-3}$ mas, which is too small to be resolved. The typical time interval between different images is roughly 1 millisecond (ms). The flux ratio between different images is roughly from 10 to $10^3$. With the expected detection rate of $10^4$ FRBs per year by the upcoming experiments, we could expect five multi-peak FRBs observed per year with time interval larger than 1 ms and flux ratio less than $10^3$ if the fraction of dark matter in MACHOs is $f\sim0.01$. A null search multiple-peak FRBs for time interval larger than 1 ms and flux ratio less than $10^3$ with $10^4$ FRBs would constrain the fraction $f$ of dark matter in MACHOs to $f<0.001$.

We propose to interpret the DAMPE electron excess at 1.5 TeV through scalar or Dirac fermion dark matter (DM) annihilation with doubly charged scalar mediators that have lepton-specific Yukawa couplings. Hierarchy of such lepton-specific Yukawa couplings is generated through the Froggatt-Nielsen mechanism, so that the dark matter annihilation products can be dominantly electrons. Stringent constraints from LEP2 on intermediate vector boson production can be evaded in our scenarios. In the case of scalar DM, we discuss one scenario with DM annihilating directly to leptons and the other scenario with DM annihilating to scalar mediators followed by their decays. We also discuss the Breit-Wigner resonant enhancement and the Sommerfeld enhancement in case that the s-wave annihilation process is small or helicity suppressed. With both types of enhancement, constraints on the parameters can be relaxed and new ways for model building will be open in explaining the DAMPE results.

Measuring the proximity effect and the damping wing of intergalactic neutral hydrogen in quasar spectra during the epoch of reionization requires an estimate of the intrinsic continuum at rest-frame wavelengths $\lambda_{\rm rest}\sim1200$-$1260$ {\AA}. In contrast to previous works which used composite spectra with matched spectral properties or explored correlations between parameters of broad emission lines, we opted for a non-parametric predictive approach based on principal component analysis (PCA) to predict the intrinsic spectrum from the spectral properties at redder (i.e. unabsorbed) wavelengths. We decomposed a sample of $12764$ spectra of $z\sim2$-$2.5$ quasars from SDSS/BOSS into 10 red-side ($1280$ {\AA} $<\lambda_{\rm rest}<2900$ {\AA}) and 6 blue-side ($1180$ {\AA} $<\lambda_{\rm rest}<1280$ {\AA}) PCA basis spectra, and constructed a projection matrix to predict the blue-side coefficients from a fit to the red-side spectrum. We found that our method predicts the blue-side continuum with $\sim6$-$12\%$ precision and $\lesssim1\%$ bias by testing on the full training set sample. We then computed predictions for the blue-side continua of the two quasars currently known at $z>7$: ULAS J1120+0641 ($z=7.09$) and ULAS J1342+0928 ($z=7.54$). Both of these quasars are known to exhibit extreme emission line properties, so we individually calibrated the precision of the continuum predictions from similar quasars in the training set. We find that both $z>7$ quasars, and in particular ULAS J1342+0928, show signs of damping wing-like absorption at wavelengths redward of Ly$\alpha$.

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.

Recently, DAMPE experiment released the new measurement of the total cosmic $e^+e^-$ flux between 25 GeV and 4.6 TeV which indicates a spectral softening at around 0.9 TeV and a tentative peak at around 1.4 TeV. We utilize the scalar dark matter (DM) annihilation scenario to explain the DAMPE peak by extending $G_{SM}\equiv SU(3)_C \times SU(2)_L \times U(1)_Y$ with additional $U(1)$ gauge symmetries while keeping anomaly free to generate $\chi \chi \to Z^\prime Z^\prime \to \ell\bar{\ell}\ell^\prime\overline{\ell^\prime}$, where $\chi, Z^\prime, \ell^{(^\prime)}$ denote the scalar DM, the new gauge boson and $\ell^{(^\prime)}=e,\mu,\tau$, respectively, with $m_\chi \sim m_{Z^\prime} \sim 2 \times 1.5$ (TeV). We first illustrate that the minimal framework $G_{SM} \times U(1)_{Y^\prime}$ with the above mass choices can explain the DAMPE excess but has been excluded by LHC constraints from the $Z^\prime$ searches. Then we study a non-minimal framework $G_{SM} \times U(1)_{Y^\prime} \times U(1)_{Y^{\prime \prime}}$ in which $U(1)_{Y^{\prime \prime}}$ mixes with $U(1)_{Y^\prime}$. We show that such a framework can interpret the DAMPE data while passing other constraints including the DM relic abundance, DM direct detection and collider bounds. We also investigate the predicted $e^+e^-$ spectrum in this framework and find that the mass splitting $\Delta m = m_\chi - m_{Z'}$ should be less than about 17 GeV to produce the peak-like structure.

As a dedicated solar radio interferometer, the MingantU SpEctral RadioHeliograph (MUSER) generates massive observational data in the frequency range of 400 MHz -- 15 GHz. High-performance imaging forms a significantly important aspect of MUSER's massive data processing requirements. In this study, we implement a practical high-performance imaging pipeline for MUSER data processing. At first, the specifications of the MUSER are introduced and its imaging requirements are analyzed. Referring to the most commonly used radio astronomy software such as CASA and MIRIAD, we then implement a high-performance imaging pipeline based on the Graphics Processing Unit (GPU) technology with respect to the current operational status of the MUSER. A series of critical algorithms and their pseudo codes, i.e., detection of the solar disk and sky brightness, automatic centering of the solar disk and estimation of the number of iterations for clean algorithms, are proposed in detail. The preliminary experimental results indicate that the proposed imaging approach significantly increases the processing performance of MUSER and generates images with high-quality, which can meet the requirements of the MUSER data processing.

In this paper, we study the luminosity function and formation rate of short gamma-ray bursts (sGRBs). Firstly, we derive the $E_p-L_p$ correlation using 16 sGRBs with redshift measurements and determine the pseudo redshifts of 284 Fermi sGRBs. Then, we use the Lynden-Bell c$^-$ method to study the luminosity function and formation rate of sGRBs without any assumptions. A strong evolution of luminosity $L(z)\propto (1+z)^{4.47}$ is found. After removing this evolution, the luminosity function is $ \Psi (L) \propto L_0 ^ {- 0.29 \pm 0.01} $ for dim sGRBs and $ \psi (L) \propto L_0 ^ {- 1.07 \pm 0.01} $ for bright sGRBs, with the break point $8.26 \times 10^{50} $ erg s$^{-1}$. We also find that the formation rate decreases rapidly at $z<1.0$, which is different with previous works. The local formation rate of sGRBs is 7.53 events Gpc$^{-3}$ yr$^{-1}$. Considering the beaming effect, the local formation rate of sGRBs including off-axis sGRBs is $ 203.31^{+1152.09}_{-135.54} $ events Gpc$^{-3}$ yr$^{-1}$. We also estimate that the event rate of sGRBs detected by the advanced LIGO and Virgo is $0.85^{+4.82}_{-0.56} $ events yr$^{-1}$ for NS-NS binary.

GRB 170817A is a weak short gamma-ray burst (GRB) accompanied by the gravitational wave (GW) event GW170817. It is believed, that an off beaming relativistic jet, produces this weak GRB. Here we use the $E_{\rm p,i}- E_{\rm iso}$ and $\Gamma - E_{\rm iso}$ relations to determine the Lorentz factor $\Gamma$ and the viewing angle from the edge of the jet $\theta'_{\rm obs}$ of GRB 170817A. Our result shows, $\Gamma = 13.4 ^{+9.8}_{-5.5}$ and $\theta'_{\rm obs} = 4.3^{+1.8}_{-1.5} \,^{\circ}$. This corresponds to the on-axis $E_{\rm p,i} = 415^{+361}_{-167}$ keV and $E_{\rm iso} = (2.4^{+1.6}_{-1.9}) \times 10^{47}$ ergs of the GRB, which is an intrinsically weak short GRB. Interestingly, the Doppler factor and the luminosity follow a universal relation from GRBs and blazars, which indicates they may share similar radiation mechanism.

GRB 170817A is a weak short gamma-ray burst (GRB) accompanied by the gravitational wave (GW) event GW170817. It is believed, that an off beaming relativistic jet, produces this weak GRB. Here we use the $E_{\rm p,i}- E_{\rm iso}$ and $\Gamma - E_{\rm iso}$ relations to determine the Lorentz factor $\Gamma$ and the viewing angle from the edge of the jet $\theta'_{\rm obs}$ of GRB 170817A. Our result shows, $\Gamma = 13.4 ^{+9.8}_{-5.5}$ and $\theta'_{\rm obs} = 4.3^{+1.8}_{-1.5} \,^{\circ}$. This corresponds to the on-axis $E_{\rm p,i} = 415^{+361}_{-167}$ keV and $E_{\rm iso} = (2.4^{+1.6}_{-1.9}) \times 10^{47}$ ergs of the GRB, which is an intrinsically weak short GRB. Interestingly, the Doppler factor and the luminosity follow a universal relation from GRBs and blazars, which indicates they may share similar radiation mechanism.

Mergers of compact binaries, such as binary neutron stars (BNSs), neutron star-black hole binaries (NSBHs), and binary black holes (BBHs), are expected to be the best candidates for the sources of gravitational waves (GWs) and the leading theoretical models for short gamma-ray bursts (SGRBs). Based on the observations of SGRBs, we could derive the merger rates of these compact binaries, and study the stochastic GW backgrounds (SGWBs) or the co-detection rates of GWs associate with SGRBs (GW-SGRBs). But before that, the most important thing is to derive the GW spectrum from a single GW source. Usually, GW spectrum from a circular orbit binary is assumed. However, observations of the large spatial offsets of SGRBs from their host galaxies imply that SGRB progenitors may be formed by the dynamical processes, and will merge with residual eccentricities. The orbital eccentricity has important effect on GW spectra, and therefore on the SGWB and GW-SGRB co-detection rate. Our results show that the power spectra of the SGWBs from eccentric compact binaries are greatly suppressed at low frequencies. Especially, SGWBs from binaries with high residual eccentricities will hard to be detected (above the detection frequency of $\sim100~\rm Hz$). For the co-detection rates of GW-SGRB events, they could be $\sim1.4$ times higher than the circular case within some particular ranges of $e_{\rm r}$ , but greatly reduced for high residual eccentricities (e.g., $e_{\rm r}>0.1$ for BNSs). In general, the BBH progenitors produce 200 and 10 times higher GW-SGRB events than the BNS and NSBH progenitors, respectively. Therefore, binaries with low residual eccentricities and high total masses will easier to be detected by aLIGO.

We simulated the $R$-band contribution of the host galaxy of TeV $\gamma$-ray BL Lac object Mrk 501 in different aperture sizes and seeing conditions . The intensive observations were run with the 1.02 m optical telescope at Yunnan Observatories from 2010 May 15 to 18. Based on the host subtraction data presented in Nilsson et al. (2007), the subtraction of host galaxy contamination results in significant seeing-brightness correlations. These correlations would lead to illusive large amplitude variations at short timescales, which will mask the intrinsic micro variability, thus gives rise to difficulty in detecting the intrinsic micro variability. Both aperture size and seeing condition influence the flux measurements, but aperture size impact the result more significantly. Based on the parameters of elliptical galaxy provided in Nilsson et al. (1999), we simulated the host contributions of Mrk 501 in the different aperture sizes and seeing conditions. Our simulation data of the host galaxy obviously weaken these significant seeing-brightness correlations for the host-subtracted brightness of Mrk 501, and can help us discover the intrinsic short timescale micro variability. The pure nuclear flux is $\sim$ 8.0 $\rm{mJy}$ in \emph{R} band, i.e., AGN has a magnitude of $R\sim13^{m}_{\cdot}96$.

In this paper, we study an anisotropic universe model with Bianchi-I metric using Joint Light-curve Analysis (JLA) sample of type Ia supernovae (SNe Ia). Because light-curve parameters of SNe Ia vary with different cosmological models and SNe Ia samples, we fit the SNe Ia light-curve parameters and cosmological parameters simultaneously employing Markov Chain Monte Carlo method. Therefore, the results on the amount of deviation from isotropy of the dark energy equation of state ($\delta$), and the level of anisotropy of the large-scale geometry ($\Sigma_0$) at present, are totally model-independent. The constraints on the skewness and cosmic shear are $-0.101<\delta<0.071$ and $-0.007<\Sigma_0<0.008$. This result is consistent with a standard isotropic universe ($\delta=\Sigma_0=0$). However, a moderate level of anisotropy in the geometry of the Universe and the equation of state of dark energy, is allowed. Besides, there is no obvious evidence for a preferred direction of anisotropic axis in this model.

Cosmic rays interact with the Earth's atmosphere to produce $^{14}$C, which can be absorbed by trees. Therefore, rapid increases of $^{14}$C in tree rings can be used to probe previous cosmic-ray events. By this method, three $^{14}$C rapidly increasing events have been found. Plausible causes of these events include large solar proton events, supernovae or short gamma-ray bursts. However, due to the lack of measurements of $^{14}$C by year, the occurrence frequency of such $^{14}$C rapidly increasing events is poorly known. In addition, rapid increases may be hidden in the IntCal13 data with five-year resolution. Here we report the result of $^{14}$C measurements using an ancient buried tree during the period between BC 3388 and 3358. We find a rapid increase of about 9\textperthousand~ in the $^{14}$C content from BC 3372 to BC 3371. We suggest that this event could originate from a large solar proton event.

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.

We present our observations of the optical intra-day variability (IDV) in $\gamma$-ray BL Lac object Mrk 501. The observations were run with the 1.02 m and 2.4 m optical telescopes at Yunnan Observatories from 2005 April to 2012 May. The light curve at the $R$ band on 2010 May 15 passes both variability tests (the $F$ test and the ANOVA test). A flare within the light curve on 2010 May 15 has a magnitude change $\Delta m = 0.03 \pm 0.005_{\rm{stat}} \pm 0.007_{\rm{sys}}$ mag, \textbf{a darkening timescale of $\tau_{\rm{d}}=$ 26.7 minutes}, and an amplitude of IDV $Amp=2.9\% \pm0.7\%$. A decline \textbf{described by 11 consecutive flux measurements} within the flare can be fitted linearly with a Pearson's correlation coefficient $r = 0.945$ at the confidence level of $> 99.99\%$. Under the assumptions that the IDV is tightly connected to the mass of the black hole, \textbf{and that the flare duration, being two times $\tau_{\rm{d}}$, is representative of the minimum characteristic timescale, we can derive upper bounds to the mass of the black hole}. In the case of the Kerr black hole, the timescale of $\Delta t_{\rm{min}}^{\rm{ob}}=$ 0.89 hours gives $M_{\bullet}\la 10^{9.20} M_{\odot}$, which is consistent with measurements reported in the literature. This agreement indicates that the hypothesis about $M_{\bullet}$ and $\Delta t_{\rm{min}}^{\rm{ob}}$ is consistent with the measurements/data.

The cosmic proper distance $d_P$ is a fundamental distance in the Universe. Unlike the luminosity and angular diameter distances, which correspond to the angular size, the proper distance is the length of light path from the source to observer. However, the proper distance has not been measured before. The recent redshift measurement of a repeat fast radio burst (FRB) can shed light on the proper distance. We show that the proper distance-redshift relation can indeed be derived from dispersion measures (DMs) of FRBs with measured redshifts. From Monte Carlo simulations, we find that about 500 FRBs with DM and redshift measurements can tightly constrain the proper distance-redshift relation. We also show that the curvature of our Universe can be constrained with a model-independent method using this derived proper distance-redshift relation and the observed angular diameter distances. Owing to the high event rate of FRBs, hundreds of FRBs can be discovered in the future by upcoming instruments. The proper distance will play an important role in investigating the accelerating expansion and the geometry of the Universe.

Type I X-ray bursts in a low-mass X-ray binary (LMXB) are caused by unstable nuclear burning of accreted materials. Semi-analytical and numerical studies of unstable nuclear burning have successfully reproduced partial properties of this kind of burst. However, some other properties (e.g. the waiting time) are not well explained. In this paper, we find that the probability distributions of fluence, peak count, rise time, duration and waiting time can be described as power-law-like distributions. This indicates that type I X-ray bursts may be governed by a self-organized criticality (SOC) process. The power-law index of waiting time distribution (WTD) is around $-1$, which is not predicted by any current waiting time model. We propose a physical burst rate model, in which the mean occurrence rate is inversely proportional to time $\lambda\propto t^{-1}$. In this case, the WTD is well explained by a non-stationary Poisson process within the SOC theory. In this theory, the burst size is also predicted to follow a power-law distribution, which requires that the emission area possesses only part of the neutron star surface. Furthermore, we find that the WTDs of some astrophysical phenomena can also be described by similar occurrence rate models.

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.

We statistically study gamma-ray burst (GRB) optical flares from the Swift/UVOT catalog. We compile 119 optical flares, including 77 flares with redshift measurements. Some tight correlations among the time scales of optical flares are found. For example, the rise time is correlated with the decay time, and the duration time is correlated with the peak time of optical flares. These two tight correlations indicate that longer rise times are associated with longer decay times of optical flares, and also suggest that broader optical flares peak at later times, which are consistent with the corresponding correlations of X-ray flares. We also study the frequency distributions of optical flare parameters, including the duration time, rise time, decay time, peak time and waiting time. Similar power-law distributions for optical and X-ray flares are found. Our statistic results imply that GRB optical flares and X-ray flares may share the similar physical origin and both of them are possibly related to central engine activities.

Previous works show that a correlation among radio luminosity, X-ray luminosity, and black hole (BH) mass from stellar-mass BHs in X-ray binaries to supermassive BHs in active galactic nuclei (AGNs), which leads to the so-called fundamental plane of BH activity. However, there are two competing explanations for this fundamental plane, including the jet-dominated model and the disk-jet model. Thus, the physical origin of this fundamental plane remains unknown. In this paper, we show that the X-ray luminosities, radio luminosities and BH masses of gamma-ray bursts (GRBs) and M82 X-1 also show a similar distribution. The universal scaling law among stellar-mass, intermediate and supermassive BH systems, together with the fact that radio and X-ray emission of GRBs originates from relativistic jets, reveals that the fundamental plane of BH activity is controlled by a jet, i.e., the radio and X-ray emission is mainly from the jet. Our work also suggests that the jets are scale-invariant with respect to the BH mass.

As a dedicated synthetic aperture radio interferometer, the MingantU SpEctral Radioheliograph (MUSER), initially known as the Chinese Spectral RadioHeliograph (CSRH), has entered the stage of routine observation. More than 23 million data records per day need to be effectively managed to provide high performance data query and retrieval for scientific data reduction. In light of these massive amounts of data generated by the MUSER, in this paper, a novel data management technique called the negative database (ND) is proposed and used to implement a data management system for the MUSER. Based on the key-value database, the ND technique makes complete utilization of the complement set of observational data to derive the requisite information. Experimental results showed that the proposed ND can significantly reduce storage volume in comparison with a relational database management system (RDBMS). Even when considering the time needed to derive records that were absent, its overall performance, including querying and deriving the data of the ND, is comparable with that of an RDBMS. The ND technique effectively solves the problem of massive data storage for the MUSER, and is a valuable reference for the massive data management required in next-generation telescopes.

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.

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

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.

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.

Fast radio bursts (FRBs) are millisecond-duration radio signals occurring at cosmological distances. However the physical model of FRBs is mystery, many models have been proposed. Here we study the frequency distributions of peak flux, fluence, duration and waiting time for the repeating FRB 121102. The cumulative distributions of peak flux, fluence and duration show power-law forms. The waiting time distribution also shows power-law distribution, and is consistent with a non-stationary Poisson process. These distributions are similar as those of soft gamma repeaters (SGRs). We also use the statistical results to test the proposed models for FRBs. These distributions are consistent with the predictions from avalanche models of slowly driven nonlinear dissipative systems.

We report the discovery of a 20-kpc-sized H{\alpha} emission in SDSS J083803.68+540642.0, a ringed dwarf galaxy (M$_V$ = -17.89 mag) hosting an accreting intermediate-mass black hole at z=0.02957. Analysis of the HST images indicates that it is an early-type galaxy with a featureless low-surface brightness disk ({\mu}0 = 20.39 mag arcsec$^{-2}$ in the V band) and a prominent, relatively red bulge (V - I = 2.03, Re = 0.28 kpc or 0."48) that accounts for ~81% of the total light in the I band. A circumgalactic ring of a diameter 16 kpc is also detected, with a disperse shape on its south side. The optical emission lines reveal the nucleus to be a broad-line LINER. Our MMT longslit observation indicates that the kinematics of the extended H{\alpha} emission is consistent with a rotational gaseous disk, with a mean blueshifted velocity of 162 km s$^{-1}$ and mean redshifted velocity of 86 km s$^{-1}$ . According to our photoionization calculations, the large-scale H{\alpha} emission is unlikely to be powered by the central nucleus or by hot evolved (post-AGB) stars interspersed in the old stellar populations, but by in situ star formation; this is vindicated by the line-ratio diagnostic of the extended emission. We propose that both the ring and large-scale H{\alpha}-emitting gas are created by the tidal accretion in a collision-and then merger-with a gas-rich galaxy of a comparable mass.

Standard candles can probe the evolution of dark energy in a large redshift range. But the cosmic opacity can degrade the quality of standard candles. In this paper, we use the latest observations, including type Ia supernovae (SNe Ia) from JLA sample and Hubble parameters, to probe the opacity of the universe. In order to avoid the cosmological dependence of SNe Ia luminosity distances, a joint fitting of the SNe Ia light-curve parameters, cosmological parameters and opacity is used. In order to explore the cosmic opacity at high redshifts, the latest gamma-ray bursts (GRBs) are used. At high redshifts, cosmic reionization process is considered. We find that the sample supports an almost transparent universe for flat $\Lambda$CDM and XCDM models. Meanwhile, free electrons deplete photons from standard candles through the (inverse) Compton scattering, known as an important component of opacity. This Compton dimming may paly an important role in future supernova surveys. From analysis, we find that about a few percent cosmic opacity is caused by Compton dimming in the two models, which can be correctable.

All next generation ground-based and space-based solar telescopes require a good quality assessment metric in order to evaluate their imaging performance. In this paper, a new image quality metric, the median filter gradient similarity (MFGS) is proposed for photospheric images. MFGS is a no-reference/blind objective image quality metric (IQM) by a measurement result between 0 and 1 and has been performed on short-exposure photospheric images captured by the New Vacuum Solar Telescope (NVST) of the Fuxian Solar Observatory and by the Solar Optical Telescope (SOT) onboard the Hinode satellite, respectively. The results show that: (1)the measured value of MFGS changes monotonically from 1 to 0 with degradation of image quality; (2)there exists a linear correlation between the measured values of MFGS and root-mean-square-contrast (RMS-contrast) of granulation; (3)MFGS is less affected by the image contents than the granular RMS-contrast. Overall, MFGS is a good alternative for the quality assessment of photospheric images.

The volume of data generated by modern astronomical telescopes is extremely large and rapidly growing. However, current high-performance data processing architectures/frameworks are not well suited for astronomers because of their limitations and programming difficulties. In this paper, we therefore present OpenCluster, an open-source distributed computing framework to support rapidly developing high-performance processing pipelines of astronomical big data. We first detail the OpenCluster design principles and implementations and present the APIs facilitated by the framework. We then demonstrate a case in which OpenCluster is used to resolve complex data processing problems for developing a pipeline for the Mingantu Ultrawide Spectral Radioheliograph. Finally, we present our OpenCluster performance evaluation. Overall, OpenCluster provides not only high fault tolerance and simple programming interfaces, but also a flexible means of scaling up the number of interacting entities. OpenCluster thereby provides an easily integrated distributed computing framework for quickly developing a high-performance data processing system of astronomical telescopes and for significantly reducing software development expenses.

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.

Many technologies have been developed to help improve spatial resolution of observational images for ground-based solar telescopes, such as adaptive optics (AO) systems and post-processing reconstruction. As any AO system correction is only partial, it is indispensable to use post-processing reconstruction techniques. In the New Vacuum Solar Telescope (NVST), speckle masking method is used to achieve the diffraction limited resolution of the telescope. Although the method is very promising, the computation is quite intensive, and the amount of data is tremendous, requiring several months to reconstruct observational data of one day on a high-end computer. To accelerate image reconstruction, we parallelize the program package on a high performance cluster. We describe parallel implementation details for several reconstruction procedures. The code is written in C language using Message Passing Interface (MPI) and optimized for parallel processing in a multi-processor environment. We show the excellent performance of parallel implementation, and the whole data processing speed is about 71 times faster than before. Finally, we analyze the scalability of the code to find possible bottlenecks, and propose several ways to further improve the parallel performance. We conclude that the presented program is capable of executing in real-time reconstruction applications at NVST.

The New Vacuum Solar Telescope (NVST) is a 1-m solar telescope that aims to observe the fine structures in both the photosphere and the chromosphere of the Sun. The observational data acquired simultaneously from one channel for the chromosphere and two channels for the photosphere bring great challenges to the data storage of NVST. The multi-channel instruments of NVST, including scientific cameras and multi-band spectrometers, generate at least 3 terabytes data per day and require high access performance while storing massive short-exposure images. It is worth studying and implementing a storage system for NVST which would balance the data availability, access performance and the cost of development. In this paper, we build a distributed data storage system (DDSS) for NVST and then deeply evaluate the availability of real-time data storage on a distributed computing environment. The experimental results show that two factors, i.e., the number of concurrent read/write and the file size, are critically important for improving the performance of data access on a distributed environment. Referring to these two factors, three strategies for storing FITS files are presented and implemented to ensure the access performance of the DDSS under conditions of multi-host write and read simultaneously. The real applications of the DDSS proves that the system is capable of meeting the requirements of NVST real-time high performance observational data storage. Our study on the DDSS is the first attempt for modern astronomical telescope systems to store real-time observational data on a low-cost distributed system. The research results and corresponding techniques of the DDSS provide a new option for designing real-time massive astronomical data storage system and will be a reference for future astronomical data storage.

The New Vacuum Solar Telescope (NVST) is a 1-meter vacuum solar telescope that aims to observe the fine structures of active regions on the Sun. The main tasks of the NVST are high resolution imaging and spectral observations, including the measurements of the solar magnetic field. The NVST has been collecting more than 20 million FITS files since it began routine observations in 2012 and produces a maximum observational records of 120 thousand files in a day. Given the large amount of files, the effective archiving and retrieval of files becomes a critical and urgent problem. In this study, we implement a new data archiving system for the NVST based on the Fastbit Not Only Structured Query Language (NoSQL) database. Comparing to the relational database (i.e., MySQL; My Structured Query Language), the Fastbit database manifests distinctive advantages on indexing and querying performance. In a large scale database of 40 million records, the multi-field combined query response time of Fastbit database is about 15 times faster and fully meets the requirements of the NVST. Our study brings a new idea for massive astronomical data archiving and would contribute to the design of data management systems for other astronomical telescopes.

The Chinese Spectral RadioHeliograph (CSRH) is a synthetic aperture radio interferometer built in Inner Mongolia, China. As a solar-dedicated interferometric array, CSRH is capable of producing high quality radio images at frequency range from 400 MHz to 15 GHz with high temporal, spatial, and spectral resolution.To implement high cadence imaging at wide-band and obtain more than 2 order higher multiple frequencies, the implementation of the data processing system for CSRH is a great challenge. It is urgent to build a pipeline for processing massive data of CSRH generated every day. In this paper, we develop a high performance distributed data processing pipeline (DDPP) built on the OpenCluster infrastructure for processing CSRH observational data including data storage, archiving, preprocessing, image reconstruction, deconvolution, and real-time monitoring. We comprehensively elaborate the system architecture of the pipeline and the implementation of each subsystem. The DDPP is automatic, robust, scalable and manageable. The processing performance under multi computers parallel and GPU hybrid system meets the requirements of CSRH data processing. The study presents an valuable reference for other radio telescopes especially aperture synthesis telescopes, and also gives an valuable contribution to the current and/or future data intensive astronomical observations.

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.

We propose a new model-independent method to test the cosmic curvature by comparing the proper distance and transverse comoving distance. Using the measurements of Hubble parameter $H(z)$ and angular diameter distance $d_A$, the cosmic curvature parameter $\Omega_K$ is constrained to be $-0.09\pm0.19$, which is consistent with a flat universe. We also use Monte Carlo simulation to test the validity and efficiency, and find that our method can give a reliable and efficient constraint on cosmic curvature. Compared with other model-independent methods testing the cosmic curvature, our method can avoid some drawbacks and give a better constraint.

Luminous quasars at z>5.6 can be studied in detail with the current generation of telescopes and provide us with unique information on the first gigayear of the universe. Thus far these studies have been statistically limited by the number of quasars known at these redshifts. Such quasars are rare and therefore wide-field surveys are required to identify them and multiwavelength data are needed to separate them efficiently from their main contaminants, the far more numerous cool dwarfs. In this paper, we update and extend the selection for z~6 quasars presented in Banados et al. (2014) using the Pan-STARRS1 (PS1) survey. We present the PS1 distant quasar sample, which currently consists of 124 quasars in the redshift range 5.6<z<6.7 that satisfy our selection criteria. Seventy-seven of these quasars have been discovered with PS1, and 63 of them are newly identified in this paper. We present composite spectra of the PS1 distant quasar sample. This sample spans a factor of ~20 in luminosity and shows a variety of emission line properties. The number of quasars at z>5.6 presented in this work almost double the quasars previously known at these redshifts, marking a transition phase from studies of individual sources to statistical studies of the high-redshift quasar population, which was impossible with earlier, smaller samples.

The Yuan-Tseh Lee Array for Microwave Background Anisotropy (AMiBA) is a co-planar interferometer array operating at a wavelength of 3mm to measure the Sunyaev-Zeldovich effect (SZE) of galaxy clusters. In the first phase of operation -- with a compact 7-element array with 0.6m antennas (AMiBA-7) -- we observed six clusters at angular scales from 5\arcmin to 23\arcmin. Here, we describe the expansion of AMiBA to a 13-element array with 1.2m antennas (AMiBA-13), its subsequent commissioning, and our cluster SZE observing program. The most important changes compared to AMiBA-7 are (1) array re-configuration with baselines ranging from 1.4m to 4.8m covering angular scales from 2\arcmin to 11.5\arcmin, (2) thirteen new lightweight carbon-fiber-reinforced plastic (CFRP) 1.2m reflectors, and (3) additional correlators and six new receivers. From the AMiBA-13 SZE observing program, we present here maps of a subset of twelve clusters. In highlights, we combine AMiBA-7 and AMiBA-13 observations of Abell 1689 and perform a joint fitting assuming a generalized NFW pressure profile. Our cylindrically integrated Compton-y values for this cluster are consistent with the BIMA/OVRA, SZA, and Planck results. We report the first targeted SZE detection towards the optically selected galaxy cluster RCS J1447+0828, and we demonstrate the ability of AMiBA SZE data to serve as a proxy for the total cluster mass. Finally, we show that our AMiBA-SZE derived cluster masses are consistent with recent lensing mass measurements in the literature.

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.

We constrain the warm dark matter (WDM) particle mass with the observations of cosmic reionization and CMB optical depth. We suggest that the GWs from stellar mass black holes (BHs) could give a further constraint on WDM particle mass for future observations. The star formation rates (SFRs) of Population I/II (Pop I/II) and Population III (Pop III) stars are also derived. If the metallicity of the universe have been enriched beyond the critical value of $Z_{\rm crit}=10^{-3.5}Z_{\odot}$, the star formation shift from Pop III to Pop I/II stars. Our results show that the SFRs are quite dependent on the WDM particle mass, especially at high redshifts. Combing with the reionization history and CMB optical depth derived from the recent \emph{Planck} mission, we find that the current data requires the WDM particle mass in a narrow range of $1 < m_x < 3$ keV. Furthermore, we suggest that the stochastic gravitational wave background (SGWB) produced by stellar BHs could give a further constraint on the WDM particle mass for future observations. For $m_{x}=3$ keV with Salpeter (Chabrier) initial mass function (IMF), the SGWB from Pop I/II BHs has a peak amplitude of $\Omega_{\rm GW}\approx2.8\times 10^{-9}~(5.0\times 10^{-9})$ at $f= 316 {\rm Hz}$, while the GW radiation at $f<10$Hz is seriously suppressed. For $m_{\rm x}=1$ keV, the SGWB peak amplitude is the same as that of $m_{\rm x}=1$ keV, but a little lower at low frequencies. Therefore, it is hard to constrain the WDM particle mass by the SGWB from Pop I/II BHs. To assess the detectability of GW signal, we also calculate the signal to noise ratio (SNR), which are $\rm SNR=37.7~ (66.5)$ and $27~(47.7)$ for $m_{\rm x}=3$ keV and $m_{\rm x}=1$ keV for Einstein Telescope (ET) with Salpeter (Chabrier) IMF, respectively.

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.

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.

G-band bright points (GBPs) are regarded as good manifestations of magnetic flux concentrations. We aim to investigate the relationship between the dynamic properties of GBPs and their longitudinal magnetic field strengths. High spatial and temporal resolution observations were recorded simultaneously with G-band filtergrams and Narrow-band Filter Imager (NFI) Stokes I and V images with Hinode /Solar Optical Telescope. The GBPs are identified and tracked in the G-band images automatically, and the corresponding longitudinal magnetic field strength of each GBP is extracted from the calibrated NFI magnetograms by a point-to-point method. After categorizing the GBPs into five groups by their longitudinal magnetic field strengths, we analyze the dynamics of GBPs of each group. The results suggest that with increasing longitudinal magnetic field strengths of GBPs correspond to a decrease in their horizontal velocities and motion ranges as well as by showing more complicated motion paths. This suggests that magnetic elements showing weaker magnetic field strengths prefer to move faster and farther along straighter paths, while stronger ones move more slowly in more erratic paths within a smaller region. The dynamic behaviors of GBPs with different longitudinal magnetic field strengths can be explained by that the stronger flux concentrations withstand the convective flows much better than weaker ones.

X-ray flares are generally supposed to be produced by the later central engine activities, and may share the similar physical origin with prompt emission of gamma-ray bursts (GRBs). In this paper, we have analyzed all significant X-ray flares from the GRBs observed by {\em Swift} from April 2005 to March 2015. The catalog contains 468 bright X-ray flares, including 200 flares with redshifts. We obtain the fitting results of X-ray flares, such as start time, peak time, duration, peak flux, fluence, peak luminosity, and mean luminosity. The peak luminosity decreases with peak time, following a power-law behavior $L_p \propto T_{peak,z}^{-1.27}$. The flare duration increases with peak time. The 0.3-10 keV isotropic energy of X-ray flares distribution is a lognormal peaked at $10^{51.2}$ erg. We also study the frequency distributions of flare parameters, including energies, durations, peak fluxes, rise times, decay times and waiting times. Power-law distributions of energies, durations, peak fluxes, and waiting times are found in GRB X-ray flares and solar flares. These distributions could be well explained by a fractal-diffusive, self-organized criticality model. Some theoretical models basing on magnetic reconnection have been proposed to explain X-ray flares. Our result shows that the relativistic jets of GRBs may be Poynting-flux dominated.

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.

In this paper, we test the discrepancy between the stellar mass density and instantaneous star formation rate in redshift range $0<z<8$ using a large observational data sample. We first compile the measurements of the stellar mass densities up to $z\sim 8$. Comparing the observed stellar mass densities with the time-integral of instantaneous star formation history, we find that the observed stellar mass densities are lower than that implied from star formation history at $z<4$. We also use Markov chain monte carlo method to derive the best-fitting star formation history from the observed stellar mass density data. At $0.5<z<6$, the observed star formation rate densities are larger than the best-fitting one, especially at $z\sim2$ where by a factor of about two. However, at lower ($z<0.5$) and higher redshifts ($z>6$), the derived star formation history is consistent with the observations. This is the first time to test the discrepancy between the observed stellar mass density and instantaneous star formation rate up to very high redshift $z\approx8$ using the Markov chain monte carlo method and a varying recycling factor. Several possible reasons for this discrepancy are discussed, such as underestimation of stellar mass density, initial mass function and cosmic metallicity evolution.

In the sixth of the series of papers reporting on a large reverberation mapping (RM) campaign of active galactic nuclei (AGNs) with high accretion rates, we present velocity-resolved time lags of H$\beta$ emission lines for nine objects observed in the campaign during 2012$-$2013. In order to correct the line-broadening caused by seeing and instruments before the analysis of velocity-resolved RM, we adopt Richardson-Lucy deconvolution to reconstruct their H$\beta$ profiles. The validity and effectiveness of the deconvolution are checked out by Monte Carlo simulation. Five among the nine objects show clear dependence of time delay on velocity. Mrk 335 and Mrk 486 show signatures of gas inflow whereas the clouds in the broad-line regions (BLRs) of Mrk 142 and MCG +06-26-012 tend to be radial outflowing. Mrk 1044 is consistent with the case of virialized motions. The lags of the rest four are not velocity-resolvable. The velocity-resolved RM of super-Eddington accreting massive black holes (SEAMBHs) shows that they have diversity of the kinematics in their BLRs. Comparing with the AGNs with sub-Eddington accretion rates, we do not find significant differences in the BLR kinematics of SEAMBHs.

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.

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

Umbral dots (UDs) are small isolated brightenings observed in sunspot umbrae. They are convective phenomena existing inside umbrae. UDs are usually divided into central UDs (CUDs) and peripheral UDs (PUDs) according to their positions inside an umbra. Our purpose is to investigate UD properties and analyze their relationships, and further to find whether or not the properties depend on umbral magnetic field strengths. Thus, we selected high-resolution TiO images of four active regions (ARs) taken under the best seeing conditions with the \textit{New Vacuum Solar Telescope} in the Fuxian Solar Observatory of the Yunnan Astronomical Observatory, China. The four ARs (NOAA 11598, 11801, 12158, and 12178) include six sunspots. A total of 1220 CUDs and 603 PUDs were identified. Meanwhile, the radial component of the vector magnetic field of the sunspots taken with the \textit{Helioseismic and Magnetic Imager} on-board the \textit{Solar Dynamics Observatory} was used to analyze relationships between UD properties and umbral magnetic field strengths. We find that diameters and lifetimes of UDs exhibit an increasing trend with the brightness, but velocities do not. Moreover, diameters, intensities, lifetimes and velocities depend on the surrounding magnetic field. A CUD diameter was found larger, the CUD brighter, its lifetime longer, and its motion slower in a weak umbral magnetic field environment than in a strong one.

Six high-resolution TiO-band image sequences from the New Vacuum Solar Telescope (NVST) are used to investigate the properties of intergranular bright points (igBPs). We detect the igBPs using a Laplacian and morphological dilation algorithm (LMD) and track them using a three-dimensional segmentation algorithm automatically, and then investigate the morphologic, photometric and dynamic properties of igBPs, in terms of equivalent diameter, the intensity contrast, lifetime, horizontal velocity, diffusion index, motion range and motion type. The statistical results confirm the previous studies based on G-band or TiO-band igBPs from the other telescopes. It illustrates that the TiO data from the NVST have a stable and reliable quality, which are suitable for studying the igBPs. In addition, our method is feasible to detect and track the igBPs in the TiO data from the NVST. With the aid of the vector magnetograms obtained from the Solar Dynamics Observatory /Helioseismic and Magnetic Imager, the properties of igBPs are found to be influenced by their embedded magnetic environments strongly. The area coverage, the size and the intensity contrast values of igBPs are generally larger in the regions with higher magnetic flux. However, the dynamics of igBPs, including the horizontal velocity, the diffusion index, the ratio of motion range and the index of motion type are generally larger in the regions with lower magnetic flux. It suggests that the absence of strong magnetic fields in the medium makes it possible for the igBPs to look smaller and weaker, diffuse faster, move faster and further in a straighter path.

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.

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.

In this paper, we use two model-independent methods to standardize long gamma-ray bursts (GRBs) using the $E_{\rm iso}-E_{\rm p}$ correlation, where $E_{\rm iso}$ is the isotropic-equivalent gamma-ray energy and $E_{\rm p}$ is the spectral peak energy. We update 42 long GRBs and try to make constraint on cosmological parameters. The full sample contains 151 long GRBs with redshifts from 0.0331 to 8.2. The first method is the simultaneous fitting method. The extrinsic scatter $\sigma_{\rm ext}$ is taken into account and assigned to the parameter $E_{\rm iso}$. The best-fitting values are $a=49.15\pm0.26$, $b=1.42\pm0.11$, $\sigma_{\rm ext}=0.34\pm0.03$ and $\Omega_m=0.79$ in the flat $\Lambda$CDM model. The constraint on $\Omega_m$ is $0.55<\Omega_m<1$ at the 1$\sigma$ confidence level. If reduced $\chi^2$ method is used, the best-fit results are $a=48.96\pm0.18$, $b=1.52\pm0.08$ and $\Omega_m=0.50\pm0.12$. The second method is using type Ia supernovae (SNe Ia) to calibrate the $E_{\rm iso}-E_{\rm p}$ correlation. We calibrate 90 high-redshift GRBs in the redshift range from 1.44 to 8.1. The cosmological constraints from these 90 GRBs are $\Omega_m=0.23^{+0.06}_{-0.04}$ for flat $\Lambda$CDM, and $\Omega_m=0.18\pm0.11$ and $\Omega_{\Lambda}=0.46\pm0.51$ for non-flat $\Lambda$CDM. For the combination of GRB and SNe Ia sample, we obtain $\Omega_m=0.271\pm0.019$ and $h=0.701\pm0.002$ for the flat $\Lambda$CDM, and for the non-flat $\Lambda$CDM, the results are $\Omega_m=0.225\pm0.044$, $\Omega_{\Lambda}=0.640\pm0.082$ and $h=0.698\pm0.004$. These results from calibrated GRBs are consistent with that of SNe Ia. Meanwhile, the combined data can improve cosmological constraints significantly, comparing to SNe Ia alone. Our results show that the $E_{\rm iso}-E_{\rm p}$ correlation is promising to probe the high-redshift universe.

We study the redshift evolution of the luminosity function (LF) and redshift selection effect of long gamma-ray bursts (LGRBs). The method is to fit the observed peak flux and redshift distributions, simultaneously. To account for the complex triggering algorithm of Swift, we use a flux triggering efficiency function. We find evidence supporting an evolving LF, where the break luminosity scales as $L_b\propto (1+z)^{\tau}$, with $\tau =3.5^{+0.4}_{-0.2}$ and $\tau =0.8^{+0.1}_{-0.08}$ for two kind of LGRB rate models. The corresponding local GRB rates are $\dot{R}(0)=0.86^{+0.11}_{-0.08} \yr^{-1}\Gpc^{-3}$ and $\dot{R}(0)= 0.54^{+0.25}_{-0.07} \yr^{-1}\Gpc^{-3}$, respectively. Furthermore, by comparing the redshift distribution between the observed one and our mocked one, we find that the redshift detection efficiency of the flux triggered GRBs decreases with redshift. Especially, a great number of GRBs miss their redshifts in the redshift range of $1<z<2.5$, where "redshift desert" effect may be dominated. More interestingly, our results show that the "redshift desert" effect is mainly introduced by the dimmer GRBs, e.g., $P<10^{-7}\ergs /\s/\cm^2$, but has no effect on the brighter GRBs.

In this paper, we obtain the NIRB and SBGWs from the early stars, which are constrained by the observation of reionization and star formation rate. We study the transition from Pop III to Pop II stars via the star formation model of different population, which takes into account the reionization and the metal enrichment evolution. We calculate the two main metal pollution channels arising from the supernova-driven protogalactic outflows and "genetic channel". We obtain the SFRs of Pop III and Pop II and their NIRB and SBGWs radiation. We predict that the upper limit of metallicity in metal-enriched IGM (the galaxies whose polluted via "genetic channel") reaches $Z_{\rm crit}=10^{-3.5}Z_{\odot}$ at $z\sim13$ ($z\sim11$), which is consistent with our star formation model. We constrain on the SFR of Pop III stars from the observation of reionization. The peak intensity of NIRB is about $0.03-0.2~nW m^{-2}{sr}^{-1}$ at $\sim 1 \mu m$ for $z>6$. The prediction of NIRB signal is consistent with the metallicity evolution. We also obtain the gravitational wave background from the black holes formed by these early stars. The predicted gravitational wave background has a peak amplitude of $\Omega_{GW}\simeq8\times10^{-9}$ at $\nu=158$ Hz for Pop II star remnants. However, the background generated by Pop III.2 stars is much less than Pop II stars, with a peak amplitude of $\Omega_{GW}\simeq1.2\times10^{-11}$ at $\nu=28~Hz$. The background of Pop III.1 shifts to lower frequencies, and the amplitude of $\Omega _{GW}$ for Pop III.1 stars shows a minimum value at $\nu\simeq 10$ Hz, due to the lack of gravitational wave signals from the stars with $140~M_{\odot}<M_\ast<260~M_{\odot}$.

Gamma-ray bursts (GRBs), which have isotropic energy up to $10^{54}$ erg, would be the ideal tool to study the properties of early universe: including dark energy, star formation rate, and the metal enrichment history of the Universe. We will briefly review the progress on the field of GRB cosmology. Meanwhile, X-ray flares, which may have important clues to the central engine, are common phenomena in the GRB afterglows. We present statistical results of X-ray flares, i.e., energy, duration time and waiting time distributions, and compare the results with solar flares. The similarity between the two kinds of flares are found, which may indicates that the physical mechanism of GRB X-ray flares is magnetic reconnection.

As one of the probes of universe, strong gravitational lensing systems allow us to compare different cosmological models and constrain vital cosmological parameters. This purpose can be reached from the dynamic and geometry properties of strong gravitational lensing systems, for instance, time-delay $\Delta\tau$ of images, the velocity dispersion $\sigma$ of the lensing galaxies and the combination of these two effects, $\Delta\tau/\sigma^2$. In this paper, in order to carry out one-on-one comparisons between $\Lambda$CDM universe and $R_h=ct$ universe, we use a sample containing 36 strong lensing systems with the measurement of velocity dispersion from the SLACS and LSD survey. Concerning the time-delay effect, 12 two-image lensing systems with $\Delta\tau$ are also used. In addition, Monte Carlo (MC) simulations are used to compare the efficiency of the three methods as mentioned above. From simulations, we estimate the number of lenses required to rule out one model at the $99.7\%$ confidence level. Comparing with constraints from $\Delta\tau$ and the velocity dispersion $\sigma$, we find that using $\Delta\tau/\sigma^2$ can improve the discrimination between cosmological models. Despite the independence tests of these methods reveal a correlation between $\Delta\tau/\sigma^2$ and $\sigma$, $\Delta\tau/\sigma^2$ could be considered as an improved method of $\sigma$ if more data samples are available.

We report the discovery of an ultra-luminous quasar J030642.51+185315.8 (hereafter J0306+1853) at redshift 5.363, which hosts a super-massive black hole (SMBH) with $M_{BH} = (1.07 \pm 0.27) \times10^{10}~M_\odot$. With an absolute magnitude $M_{1450}=-28.92$ and bolometric luminosity $L_{bol}\sim3.4\times10^{14} L_{\odot}$, J0306+1853 is one of the most luminous objects in the early Universe. It is not likely to be a beamed source based on its small flux variability, low radio loudness and normal broad emission lines. In addition, a $z=4.986$ Damped Ly$\alpha$ system (DLA) with $\rm [M/H]=-1.3\pm0.1$, among the most metal rich DLAs at $z \gtrsim 5$, is detected in the absorption spectrum of this quasar. This ultra-luminous quasar puts strong constraint on the bright-end of quasar luminosity function and massive-end of black hole mass function. It will provide a unique laboratory to the study of BH growth and the co-evolution between BH and host galaxy with multi-wavelength follow-up observations. The future high resolution spectra will give more insights to the DLA and other absorption systems along the line-of-sight of J0306+1853.

From generalized gravity mediation we build a SUGRA scenario in which the gluino is much heavier than the electroweak gauginos at the GUT scale. We find that such a non-universal gaugino scenario with very heavy gluino at the GUT scale can be naturally obtained with proper high dimensional operators in the framework of SU(5) GUT. Then, due to the effects of heavy gluino, at the weak scale all colored sparticles are heavy while the uncolored sparticles are light, which can explain the Brookhaven muon g-2 measurement while satisfying the collider constraints (both the 125 GeV Higgs mass and the direct search limits of sparticles) and dark matter requirements. We also find that, in order to explain the muon g-2 measurement, the neutralino dark matter is lighter than 200 GeV in our scenario, which can be mostly covered by the future Xenon1T experiment.

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

X-ray flares were discovered in the afterglow phase of gamma-ray bursts (GRBs) by the {\em Swift} satellite a decade ago and known as a canonical component in GRB X-ray afterglows. In this paper, we constrain the Lorentz factors of GRB X-ray flares using two different methods. For the first method, we estimate the lower limit on the bulk Lorentz factor with the flare duration and jet break time. In the second method, the upper limit on the Lorentz factor is derived by assuming that the X-ray flare jet has undergone saturated acceleration. We also re-estimate the initial Lorentz factor with GRB afterglow onsets, and find the coefficient of the theoretical Lorentz factor is 1.67 rather than the commonly used 2 for interstellar medium (ISM) and 1.44 for the wind case. We find that the correlation between the limited Lorentz factor and the isotropic radiation energy of X-ray flares in the ISM case is more consistent with that of prompt emission than the wind case in a statistical sense. For a comparison, the lower limit on Lorentz factor is statistically larger than the extrapolation from prompt bursts in the wind case. Our results indicate that X-ray flares and prompt bursts are produced by the same physical mechanism.

During the "WISE at 5: Legacy and Prospects" conference in Pasadena, CA -- which ran from February 10 - 12, 2015 -- attendees were invited to engage in an interactive session exploring the future uses of the Wide-field Infrared Survey Explorer (WISE) data. The 65 participants -- many of whom are extensive users of the data -- brainstormed the top questions still to be answered by the mission, as well as the complementary current and future datasets and additional processing of WISE/NEOWISE data that would aid in addressing these most important scientific questions. The results were mainly bifurcated between topics related to extragalactic studies (e.g. AGN, QSOs) and substellar mass objects. In summary, participants found that complementing WISE/NEOWISE data with cross-correlated multiwavelength surveys (e.g. SDSS, Pan-STARRS, LSST, Gaia, Euclid, etc.) would be highly beneficial for all future mission goals. Moreover, developing or implementing machine-learning tools to comb through and understand cross-correlated data was often mentioned for future uses. Finally, attendees agreed that additional processing of the data such as co-adding WISE and NEOWISE and extracting a multi-epoch photometric database and parallax and proper motion catalog would greatly improve the scientific results of the most important projects identified. In that respect, a project such as MaxWISE which would execute the most important additional processing and extraction as well as make the data and catalogs easily accessible via a public portal was deemed extremely important.

The Large sky Area Multi-Object Spectroscopic Telescope (LAMOST) General Survey is a spectroscopic survey that will eventually cover approximately half of the celestial sphere and collect 10 million spectra of stars, galaxies and QSOs. Objects both in the pilot survey and the first year general survey are included in the LAMOST First Data Release (DR1). The pilot survey started in October 2011 and ended in June 2012, and the data have been released to the public as the LAMOST Pilot Data Release in August 2012. The general survey started in September 2012, and completed its first year of operation in June 2013. The LAMOST DR1 includes a total of 1202 plates containing 2,955,336 spectra, of which 1,790,879 spectra have observed signal-to-noise S/N >10. All data with S/N>2 are formally released as LAMOST DR1 under the LAMOST data policy. This data release contains a total of 2,204,696 spectra, of which 1,944,329 are stellar spectra, 12,082 are galaxy spectra and 5,017 are quasars. The DR1 includes not only spectra, but also three stellar catalogues with measured parameters: AFGK-type stars with high quality spectra (1,061,918 entries), A-type stars (100,073 entries), and M stars (121,522 entries). This paper introduces the survey design, the observational and instrumental limitations, data reduction and analysis, and some caveats. Description of the FITS structure of spectral files and parameter catalogues is also provided.

In this paper, we test the dark matter-dark energy interacting cosmological model with a dynamic equation of state $w_{DE}(z)=w_{0}+w_{1}z/(1+z)$, using type Ia supernovae (SNe Ia), Hubble parameter data, baryonic acoustic oscillation (BAO) measurements, and the cosmic microwave background (CMB) observation. This interacting cosmological model has not been studied before. The best-fitted parameters with $1 \sigma$ uncertainties are $\delta=-0.022 \pm 0.006$, $\Omega_{DM}^{0}=0.213 \pm 0.008$, $w_0 =-1.210 \pm 0.033$ and $w_1=0.872 \pm 0.072$ with $\chi^2_{min}/dof = 0.990$. At the $1 \sigma$ confidence level, we find $\delta<0$, which means that the energy transfer prefers from dark matter to dark energy. We also find that the SNe Ia are in tension with the combination of CMB, BAO and Hubble parameter data. The evolution of $\rho_{DM}/\rho_{DE}$ indicates that this interacting model is a good approach to solve the coincidence problem, because the $\rho_{DE}$ decrease with scale factor $a$. The transition redshift is $z_{tr}=0.63 \pm 0.07$ in this model.

Gamma-ray bursts (GRBs) are the most violent explosions in the Universe and can be used to explore the properties of high-redshift universe. It is believed that the long GRBs are associated with the deaths of massive stars. So it is possible to use GRBs to investigate the star formation rate (SFR). In this paper, we use Lynden-Bell's $c^-$ method to study the luminosity function and rate of \emph{Swift} long GRBs without any assumptions. We find that the luminosity of GRBs evolves with redshift as $L(z)\propto g(z)=(1+z)^k$ with $k=2.43_{-0.38}^{+0.41}$. After correcting the redshift evolution through $L_0(z)=L(z)/g(z)$, the luminosity function can be expressed as $\psi(L_0)\propto L_0^{-0.14\pm0.02}$ for dim GRBs and $\psi(L_0)\propto L_0^{-0.70\pm0.03}$ for bright GRBs, with the break point $L_{0}^{b}=1.43\times10^{51}~{\rm erg~s^{-1}}$. We also find that the formation rate of GRBs is almost constant at $z<1.0$ for the first time, which is remarkably different from the SFR. At $z>1.0$, the formation rate of GRB is consistent with the SFR. Our results are dramatically different from previous studies. Some possible reasons for this low-redshift excess are discussed. We also test the robustness of our results with Monte Carlo simulations. The distributions of mock data (i.e., luminosity-redshift distribution, luminosity function, cumulative distribution and $\log N-\log S$ distribution) are in good agreement with the observations. Besides, we also find that there are remarkable difference between the mock data and the observations if long GRB are unbiased tracers of SFR at $z<1.0$.

Gamma-ray bursts (GRBs) are the most luminous electromagnetic explosions in the Universe, which emit up to $8.8\times10^{54}$ erg isotropic equivalent energy in the hard X-ray band. The high luminosity makes them detectable out to the largest distances yet explored in the Universe. GRBs, as bright beacons in the deep Universe, would be the ideal tool to probe the properties of high-redshift universe: including the cosmic expansion and dark energy, star formation rate, the reionization epoch and the metal enrichment history of the Universe. In this article, we review the luminosity correlations of GRBs, and implications for constraining the cosmological parameters and dark energy. Observations show that the progenitors of long GRBs are massive stars. So it is expected that long GRBs are tracers of star formation rate. We also review the high-redshift star formation rate derived from GRBs, and implications for the cosmic reionization history. The afterglows of GRBs generally have broken power-law spectra, so it is possible to extract intergalactic medium (IGM) absorption features. We also present the capability of high-redshift GRBs to probe the pre-galactic metal enrichment and the first stars.

This is the third in a series of papers reporting on a large reverberation-mapping campaign aimed to study the properties of active galactic nuclei (AGNs) with high accretion rates. We present new results on the variability of the optical Fe II emission lines in 10 AGNs observed by the Yunnan Observatory 2.4m telescope during 2012--2013. We detect statistically significant time lags, relative to the AGN continuum, in nine of the sources. This accurate measurement is achieved by using a sophisticated spectral fitting scheme that allows for apparent flux variations of the host galaxy, and several narrow lines, due to the changing observing conditions. Six of the newly detected lags are indistinguishable from the Hbeta lags measured in the same sources. Two are significantly longer and one is slightly shorter. Combining with Fe II lags reported in previous studies, we find a Fe II radius--luminosity relationship similar to the one for Hbeta, although our sample by itself shows no clear correlation. The results support the idea that Fe II emission lines originate in photoionized gas which, for the majority of the newly reported objects, is indistinguishable from the Hbeta-emitting gas. We also present a tentative correlation between the lag and intensity of Fe II and Hbeta and comment on its possible origin.

So far, roughly 40 quasars with redshifts greater than z=6 have been discovered. Each quasar contains a black hole with a mass of about one billion solar masses ($10^9 M_\odot$). The existence of such black holes when the Universe was less than 1 billion years old presents substantial challenges to theories of the formation and growth of black holes and the coevolution of black holes and galaxies. Here we report the discovery of an ultra-luminous quasar, SDSS J010013.02+280225.8, at redshift z=6.30. It has an optical and near-infrared luminosity a few times greater than those of previously known z>6 quasars. On the basis of the deep absorption trough on the blue side of the Ly $\alpha$ emission line in the spectrum, we estimate the proper size of the ionized proximity zone associated with the quasar to be 26 million light years, larger than found with other z>6.1 quasars with lower luminosities. We estimate (on the basis of a near-infrared spectrum) that the black hole has a mass of $\sim 1.2 \times 10^{10} M_\odot$, which is consistent with the $1.3 \times 10^{10} M_\odot$ derived by assuming an Eddington-limited accretion rate.

We propose a direct and model-independent method to constrain the Lorentz factor of a relativistically expanding object, like gamma-ray bursts. Only the measurements, such as thermal component of the emission, the distance and the variable time scale of the light curve, are used. If the uncertainties are considered, we will obtain lower limits of the Lorentz factor instead. We apply this method to GRB 090618 and get a lower limit of the Lorentz factor to be 22. The method can be used to any relativistically moving object, such as gamma-ray bursts, blazars, and soft gamma-ray repeaters, providing the thermal component of the emission being observed.

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

X-ray flares have been discovered in black hole systems, such as gamma-ray bursts, the tidal disruption event Swift J1644+57, the supermassive black hole Sagittarius A$^*$ at the center of our Galaxy, and some active galactic nuclei. Their occurrences are always companied by relativistic jets. However, it is still unknown whether there is a physical analogy among such X-ray flares produced in black hole systems spanning nine orders of magnitude in mass. Here we report the observed data of X-ray flares, and show that they have three statistical properties similar to solar flares, including power-law distributions of energies, durations, and waiting times, which both can be explained by a fractal-diffusive self-organized criticality model. These statistical similarities, together with the fact that solar flares are triggered by a magnetic reconnection process, suggest that all of the X-ray flares are consistent with magnetic reconnection events, implying that their concomitant relativistic jets may be magnetically dominated.

We recalculate the two-loop beta functions for three gauge couplings taking into account all low energy threshold corrections in split supersymmetry (split-SUSY) which assumes a very high scalar mass scale M_S. We find that, in split-SUSY with gaugino mass unification assumption and a large M_S, the gauge coupling unification requires a lower bound on gaugino mass. Combined with the constraints from the dark matter relic density and direct detection limits, we find that split-SUSY is very restricted and for dark matter mass below 1 TeV the allowed parameter space can be fully covered by XENON-1T(2017).

Only very few z>5 quasars discovered to date are radio-loud, with a radio-to-optical flux ratio (radio-loudness parameter) higher than 10. Here we report the discovery of an optically luminous radio-loud quasar, SDSS J013127.34-032100.1 (J0131-0321 in short), at z=5.18+-0.01 using the Lijiang 2.4m and Magellan telescopes. J0131-0321 has a spectral energy distribution consistent with that of radio-loud quasars. With an i-band magnitude of 18.47 and radio flux density of 33 mJy, its radio-loudness parameter is ~100. The optical and near-infrared spectra taken by Magellan enable us to estimate its bolometric luminosity to be L_bol ~ 1.1E48 erg/s, approximately 4.5 times greater than that of the most distant quasar known to date. The black hole mass of J0131-0321 is estimated to be 2.7E9 solar masses, with an uncertainty up to 0.4 dex. Detailed physical properties of this high-redshift, radio-loud, potentially super-Eddington quasar can be probed in the future with more dedicated and intensive follow-up observations using multi-wavelength facilities.

Only very few z>5 quasars discovered to date are radio-loud, with a radio-to-optical flux ratio (radio-loudness parameter) higher than 10. Here we report the discovery of an optically luminous radio-loud quasar, SDSS J013127.34-032100.1 (J0131-0321 in short), at z=5.18+-0.01 using the Lijiang 2.4m and Magellan telescopes. J0131-0321 has a spectral energy distribution consistent with that of radio-loud quasars. With an i-band magnitude of 18.47 and radio flux density of 33 mJy, its radio-loudness parameter is ~100. The optical and near-infrared spectra taken by Magellan enable us to estimate its bolometric luminosity to be L_bol ~ 1.1E48 erg/s, approximately 4.5 times greater than that of the most distant quasar known to date. The black hole mass of J0131-0321 is estimated to be 2.7E9 solar masses, with an uncertainty up to 0.4 dex. Detailed physical properties of this high-redshift, radio-loud, potentially super-Eddington quasar can be probed in the future with more dedicated and intensive follow-up observations using multi-wavelength facilities.

Many dark energy models fail to pass the cosmic age test. In this paper, we investigate the cosmic age problem associated with nine extremely old Global Clusters (GCs) and the old quasar APM $08279+5255$ in the $R_h=ct$ Universe. The age data of these oldest GCs in M31 is acquired from the Beijing-Arizona-Taiwan-Connecticut system with up-to-date theoretical synthesis models. They have not been used to test the cosmic age problem in the $R_h=ct$ Universe in previous literature. By evaluating the age of the $R_h=ct$ Universe with the observational constraints from the type Ia supernovae and Hubble parameter, we find that the $R_h=ct$ Universe can accommodate five GCs and the quasar APM 08279+5255 at redshift $z=3.91$. But for other models, such as $\Lambda$CDM, interacting dark energy model, generalized Chaplygin gas model and holographic dark energy model, can not accommodate all GCs and the quasar APM 08279+5255. It is worthwhile to note that the age estimates of some GCs are controversial. So, unlike other cosmological models, the $R_h=ct$ Universe can marginally solve the cosmic age problem, especially at high redshift.

In this work we show that the general singlet extension of the MSSM can naturally provide a self-interacting singlino dark matter to solve the small cosmological scale anomalies (a large Sommerfeld enhancement factor can also be obtained). However, we find that the NMSSM (the singlet extension of the MSSM with Z3 symmetry) cannot achieve this due to the restricted parameter space. In our analysis we introduce the concept of symmetric and antisymmetric viscosity cross sections to deal with the non-relativistic Majorana-fermion dark matter scattering.

This is the second in a series of papers reporting on a large reverberation mapping (RM) campaign to measure black hole (BH) mass in high accretion rate active galactic nuclei (AGNs). The goal is to identify super-Eddington accreting massive black holes (SEAMBHs) and to use their unique properties to construct a new method for measuring cosmological distances. Based on theoretical models, the saturated bolometric luminosity of such sources is proportional to the BH mass which can be used to obtain their distance. Here we report on five new RM measurements and show that in four of the cases we can measure the BH mass and three of these sources are SEAMBHs. Together with the three sources from our earlier work, we now have six new sources of this type. We use a novel method based on a minimal radiation efficiency to identify nine additional SEAMBHs from earlier RM-based mass measurements. We use a Bayesian analysis to determine the parameters of the new distance expression, and the method uncertainties, from the observed properties of the objects in the sample. The ratio of the newly measured distances to the standard cosmological ones has a mean scatter of 0.14 dex, indicating that SEAMBHs can be use as cosmological distance probes. With their high luminosity, long period of activity and large numbers at high redshifts, SEAMBHs have a potential to extend the cosmic distance ladder beyond the range now explored by type Ia supernovae.

We study the motions of G band bright points (GBPs) in the quiet Sun to obtain the characteristics of different motion types. A high resolution image sequence taken with the Hinode/Solar Optical Telescope (SOT) is used, and GBPs are automatically tracked by segmenting 3D evolutional structures in a space time cube. After putting the GBPs that do not move during their lifetimes aside, the non stationary GBPs are categorized into three types based on an index of motion type. Most GBPs that move in straight or nearly straight lines are categorized into a straight motion type, a few moving in rotary paths into a rotary motion, and the others fall into a motion type we called erratic. The mean horizontal velocity is 2.18 km/s, 1.63 km/s and 1.33 km/s for straight, erratic and rotary motion type, respectively. We find that a GBP drifts at a higher and constant velocity during its whole life if it moves in a straight line. However, it has a lower and variational velocity if it moves in a rotary path. The diffusive process is ballistic, super and sub diffusion for straight, erratic and rotary motion type, respectively. The corresponding diffusion index and coefficients are 2.13 and 850 km2/s, 1.82 and 331 km2/s, 0.73 and 13 km2/s. In terms of direction of motion, it is homogeneous and isotropical, and usually persists between neighbouring frames, no matter what motion type a GBP belongs to.

In this paper, we study the anisotropic expansion of the universe using type Ia supernovae Union 2.1 sample and 116 long gamma-ray bursts. The luminosity distance is expanded with model-independent cosmographic parameters as a function of $z/(1+z)$ directly. Thus the results are independent of cosmology model. We find a dipolar anisotropy in the direction ($l=309.2^\circ \pm 15.8^\circ$, $b=-8.6^\circ \pm 10.5^\circ$) in galactic coordinates with a significant evidence $97.29\%$ (more than $2~\sigma$). The magnitude is $(1.37\pm 0.57) \times 10^{-3}$ for the dipole, and $(2.6\pm 2.1)\times 10^{-4}$ for the monopole, respectively. This dipolar anisotropy is more significant at low redshift from the redshift tomography analysis. We also test whether this preferred direction is caused by bulk flow motion or dark energy dipolar scalar perturbation. We find that the direction and the amplitude of the bulk flow in our results are approximately consistent with the bulk flow surveys. Therefore, bulk flow motion may be the main reason for the anisotropic expansion at low redshift, but the effect of dipolar distribution dark energy can not be excluded, especially at high redshift.

We investigate the mass distribution of long gamma-ray burst (GRB) host galaxies and the redshift distribution of long GRBs by considering that long GRBs occur in low-metallicity environments. We calculate the upper limit on the stellar mass of a galaxy which can produce long GRBs by utilizing the mass-metallicity (M-Z) relation of galaxies. After comparing with the observed GRB host galaxies masses, we find that the observed GRB host galaxy masses can fit the predicted masses well if GRBs occur in low-metallicity $12+\log\rm(O/H)_{\rm KK04}<8.7$. GRB host galaxies have low metallicity, low mass, and high star formation rate compared with galaxies of seventh data release of the Sloan Digital Sky Survey. We also study the cumulative redshift distribution of the latest \emph{Swift} long GRBs by adding dark GRBs and 10 new GRBs redshifts from TOUGH survey. The observed discrepancy between the GRB rate and the star formation history can be reconciled by considering that GRBs tend to occur in low-metallicity galaxies with $12+\log\rm(O/H)_{\rm KK04}<8.7$. We conclude that the metallicity cutoff that can produce long GRBs is about $12+\log\rm(O/H)_{\rm KK04}<8.7$ from the host mass distribution and redshift distribution.

Active galactic nuclei (AGNs) and gamma-ray bursts (GRBs) are powerful astrophysical events with relativistic jets. In this Letter the broadband spectral properties are compared between GRBs and the well-observed blazars. The distribution of GRBs are consistent with the well-known blazar sequence including the $\nu L_\nu(5\rm GHz)-\alpha_{\rm RX}$ and $\nu L_\nu(5\rm GHz)-\nu_{\rm peak}$ correlations, where $\alpha_{\rm RX}$ is defined as the broadband spectral slope in radio-to-X-ray bands, and $\nu_{\rm peak}$ is defined as the spectral peak frequency. Moreover, GRBs occupy the low radio luminosity end of these sequences. These two correlations suggest that GRBs could have a similar radiation process with blazars both in the prompt emission and afterglow phases, i.e., synchrotron radiation.

Several laboratory facilities were used to benchmark theoretical spectral models those extensively used by astronomical communities. However there are still many differences between astrophysical environments and laboratory miniatures that can be archived. Here we setup a spectral analysis system for astrophysical and laboratory (SASAL) plasmas to make a bridge between them, and investigate the effects from non-thermal electrons, contribution from metastable level-population on level populations and charge stage distribution for coronal-like, photoionized, and geocoronal plasmas. Test applications to laboratory measurement (i.e. EBIT plasma) and astrophysical observation (i.e. Comet, Cygnus X-3) are presented. Time evolution of charge stage and level population are also explored for collisional and photoionized plasmas.

The collapsar model of long gamma-ray bursts (GRBs) indicates that they may trace the star formation history. So long GRBs may be a useful tool of measuring the high-redshift star formation rate (SFR). The collapsar model explains GRB formation via the collapse of a rapidly rotating massive star with $M> 30M_\odot$ into a black hole, which may imply a decrease of SFR at high redshift. However, we find that the \emph{Swift} GRBs during 2005-2012 are biased tracing the SFR, including a factor about $(1+z)^{0.5}$, which is in agreement with recent results. After taking this factor, the SFR derived from GRBs does not show steep drop up to $z\sim 9.4$. We consider the GRBs produced by rapidly rotating metal-poor stars with low masses to explain the high-redshift GRB rate excess. The chemically homogeneous evolution scenario (CHES) of rapidly rotating stars with mass larger than $12M_\odot$ is recognized as a promising path towards collapsars in connection with long GRBs. Our results indicate that the stars in the mass range $12M_\odot<M<30M_\odot$ for low enough metallicity $Z\leq 0.004$ with the GRB efficiency factor $10^{-5}$ can fit the derived SFR with good accuracy. Combining these two factors, we find that the conversion efficiency from massive stars to GRBs is enhanced by a factor of 10, which may be able to explain the excess of the high-redshift GRB rate. We also investigate the cosmic reionization history using the derived SFR. The GRB-inferred SFR would be sufficient to maintain cosmic reionization over $6<z<10$ and reproduce the observed optical depth of Thomson scattering to the cosmic microwave background.

Many long-duration gamma-ray bursts (GRBs) were observed by {\it Swift}/XRT to have plateaus in their X-ray afterglow light curves. This plateau phase has been argued to be evidence for long-lasting activity of magnetar (ultra-strongly magnetized neutron stars) central engines. However, the emission efficiency of such magnetars in X-rays is still unknown. Here we collect 24 long GRB X-ray afterglows showing plateaus followed by steep decays. We extend the well-known relationship between the X-ray luminosity ${L_{\mathrm{X}}}$ and spin-down luminosity $L_{\mathrm{sd}}$ of pulsars to magnetar central engines, and find that the initial rotation period $P_{0}$ ranges from 1 ms to 10 ms and that the dipole magnetic field $B$ is centered around $10^{15}$ G. These constraints not only favor the suggestion that the central engines of some long GRBs are very likely to be rapidly rotating magnetars but also indicate that the magnetar plateau emission efficiency in X-rays is close to 100%.

The equation of state (EOS), $w(z)$, is the most important parameter of dark energy. We reconstruct the evolution of this EOS in a model-independent way using the latest cosmic microwave background (CMB) data from Planck and other observations, such as type Ia supernovae (SNe Ia), the baryonic acoustic oscillation measurements (SDSS, 6dF, BOSS, and WiggleZ), and the Hubble parameter value $H(z)$. The results show that the EOS is consistent with the cosmological constant at the $2\sigma$ confidence level, not preferring a dynamical dark energy. The uncorrelated EOS of dark energy constraints from Planck CMB data are much tighter than those from the WMAP 9-year CMB data.

We report first results from a large project to measure black hole (BH) mass in high accretion rate active galactic nuclei (AGNs). Such objects may be different from other AGNs in being powered by slim accretion disks and showing saturated accretion luminosities, but both are not yet fully understood. The results are part of a large reverberation mapping (RM) campaign using the 2.4-m Shangri-La telescope at the Yunnan Observatory in China. The goals are to investigate the gas distribution near the BH and the properties of the central accretion disks, to measure BH mass and Eddington ratios, and to test the feasibility of using such objects as a new type of cosmological candles. The paper presents results for three objects, Mrk 335, Mrk 142 and IRAS F12397+3333 with H$\beta$ time lags relative to the 5100\AA\ continuum of $10.6^{+1.7}_{-2.9}$, $6.4^{+0.8}_{-2.2}$ and $11.4^{+2.9}_{-1.9}$ days, respectively. The corresponding BH masses are $(8.3_{-3.2}^{+2.6})\times 10^6M_{\odot}$, $(3.4_{-1.2}^{+0.5})\times 10^6M_{\odot}$ and $(7.5_{-4.1}^{+4.3})\times 10^6M_{\odot}$, and the lower limits on the Eddington ratios 0.6, 2.3, and 4.6 for the minimal radiative efficiency of 0.038. Mrk 142 and IRAS F12397+333 (extinction corrected) clearly deviate from the currently known relation between H$\beta$ lag and continuum luminosity. The three Eddington ratios are beyond the values expected in thin accretion disks and two of them are the largest measured so far among objects with RM-based BH masses. We briefly discuss implications for slim disks, BH growth and cosmology.

A cosmological preferred direction was reported from the type Ia supernovae (SNe Ia) data in recent years. We use the Union2.1 data to give a simple classification of such studies for the first time. Because the maximum anisotropic direction is independent of isotropic dark energy models, we adopt two cosmological models ($\Lambda$CDM, $w$CDM) for the hemisphere comparison analysis and $\Lambda$CDM model for dipole fit approach. In hemisphere comparison method, the matter density and the equation of state of dark energy are adopted as the diagnostic qualities in the $\Lambda$CDM model and $w$CDM model, respectively. In dipole fit approach, we fit the fluctuation of distance modulus. We find that there is a null signal for the hemisphere comparison method, while a preferred direction ($b=-14.3^\circ \pm 10.1^\circ, l=307.1^\circ \pm 16.2^\circ$) for the dipole fit method. This result indicates that the dipole fit is more sensitive than the hemisphere comparison method.

Identifications of quasars at intermediate redshifts (2.2<z<3.5) are inefficient in most previous quasar surveys as their optical colors are similar to those of stars. The near-IR K-band excess technique has been suggested to overcome this difficulty. Our study also proposed to use optical/near-IR colors for selecting z<4 quasars. To this method, we selected 105 unidentified bright targets with i<18.5 from the quasar candidates of SDSS DR6 with both SDSS ugriz optical and UKIDSS YJHK near-IR photometric data, which satisfy our proposed Y-K/g-z criterion and have photometric redshifts between 2.2 and 3.5 estimated from the 9-band SDSS-UKIDSS data. 43 of them were observed with the 2.16m telescope of NAOC in 2012. 36 of them were identified as quasars at 2.1<z<3.4. High success rate of discovering these quasars in the SDSS spectroscopic surveyed area demonstrates the robustness of both the Y-K/g-z selection criterion and the photometric redshift estimation technique. We also used the above criterion to investigate the possible star contamination rate to the quasar candidates of SDSS DR6, and found that it is much higher in selecting 3<z<3.5 quasar candidates than selecting lower redshift ones (z<2.2). The significant improvement in the photometric redshift estimation by using the 9-band SDSS-UKIDSS data than using the 5-band SDSS data is demonstrated and a catalog of 7,727 unidentified quasar candidates with photometric redshifts between 2.2 and 3.5 is provided. We also tested the Y-K/g-z selection criterion with the SDSS-III/DR9 quasar catalog, and found 96.2% of 17,999 DR9 quasars with UKIDSS Y and K-band data satisfy our criterion. With some samples of red and type II quasars, we found that 88% and 96.5% of them can be selected by the Y-K/g-z criterion respectively, which supports that using the Y-K/g-z criterion we can efficiently select both unobscured and obscured quasars. (abridged)

X-ray flares detected in nearly half of gamma-ray burst (GRB) afterglows are one of the most intriguing phenomena in high-energy astrophysics. All the observations indicate that the central engines of bursts, after the gamma-ray emission has ended, still have long periods of activity, during which energetic explosions eject relativistic materials, leading to late-time X-ray emission. It is thus expected that X-ray flares provide important clues to the nature of the central engines of GRBs, and more importantly, unveil the physical mechanism of the flares themselves, which has so far remained mysterious. Here we report statistical results of X-ray flares of GRBs with known redshifts, and show that X-ray flares and solar flares share three statistical properties: power-law frequency distributions for energies, durations, and waiting times. All of the distributions can be well understood within the physical framework of a self-organized criticality (SOC) system. The statistical properties of X-ray flares of GRBs are similar to solar flares, and thus can both be attributed to a SOC process. Both types of flares may be driven by a magnetic reconnection process, but X-ray flares of GRBs are produced in ultra-strongly magnetized millisecond pulsars or long-term hyperaccreting disks around stellar-mass black holes.

The X-ray emission from Swift J1644+57 is not steadily decreasing instead it shows multiple pulses with declining amplitudes. We model the pulses as reverse shocks from collisions between the late ejected shells and the externally shocked material, which is decelerated while sweeping the ambient medium. The peak of each pulse is taken as the maximum emission of each reverse shock. With a proper set of parameters, the envelope of peaks in the light curve as well as the spectrum can be modelled nicely.

We present a spectroscopic catalog of 67082 M dwarfs from the LAMOST pilot survey. For each spectrum of the catalog, spectral subtype, radial velocity, equivalent width of H${\alpha}$, a number of prominent molecular band indices and the metal sensitive parameter $\zeta$ are provided . Spectral subtype have been estimated by a remedied Hammer program (Original Hammer: Covey et al. 2007), in which indices are reselected to obtain more accurate auto-classified spectral subtypes. All spectra in this catalog have been visually inspected to confirm the spectral subtypes. Radial velocities have been well measured by our developed program which uses cross-correlation method and estimates uncertainty of radial velocity as well. We also examine the magnetic activity properties of M dwarfs traced by H${\alpha}$ emission line. The molecular band indices included in this catalog are temperature or metallicity sensitive and can be used for future analysis of the physical properties of M dwarfs. The catalog is available on the website \url{http://sciwiki.lamost.org/MCatalogPilot/}.

The LAMOST spectroscopic survey for Galactic structure and evolution has been in operation since October 2012, following a one-year pilot survey. The pilot survey produced a data release containing over 600,000 stellar spectra. By cross-checking with a large time series photometric database of RR Lyrae stars in high Galactic latitude regions, we found a total number of 157 stars that have been observed with LAMOST. In this sample, we successfully captured three RR Lyrae stars in the fast expansion phase, all of them showing significant hypersonic shock wave features in the spectra. By fitting the H{\alpha} line shape, we determine that the emission features seen within the broader H{\alpha} absorption line suggest hypersonic relative motions in the atmospheres of these three objects. With a further LAMOST survey of millions of stars, we will capture a large sample of RR Lyrae stars in their hypersonic expansion phase, and therefore provide a large database for the study of the internal structure and pulsation mechanism of RR Lyrae stars.

Spectrograph slit is conventionally used to enhance the spectral resolution and manage how much light can be allowed to enter spectrograph. A narrower slit provides a higher resolution but sacrifices efficiency of spectrograph and results in a low signal to noise ratio (S/N) spectra product. We take GuoShouJing telescope as an example and carry out a series of experiments to study how its 2/3 slit mode affects the precision of stellar radial velocity measurement and atmosphere parameters estimate. By transforming the resolution and adding a Gaussian White Noise to the extremely high quality spectra from the Sloan Digital Sky Survey, we generate synthetic stellar spectra of various brightness with different S/Ns. Comparing the measurements on these noise added spectra with the original high quality ones, we summarize the influences of the 2/3 slit mode on the measurement accuracy of stellar radial velocity and atmospheric parameters.

We present a new way to solve the platform deformation problem of co-planar interferometers. The platform of a co-planar interferometer can be deformed due to driving forces and gravity. A deformed platform will induce extra components into the geometric delay of each baseline, and change the phases of observed visibilities. The reconstructed images will also be diluted due to the errors of the phases. The platform deformations of The Yuan-Tseh Lee Array for Microwave Background Anisotropy (AMiBA) were modelled based on photogrammetry data with about 20 mount pointing positions. We then used the differential optical pointing error between two optical telescopes to fit the model parameters in the entire horizontal coordinate space. With the platform deformation model, we can predict the errors of the geometric phase delays due to platform deformation with given azimuth and elevation of the targets and calibrators. After correcting the phases of the radio point sources in the AMiBA interferometric data, we recover 50% - 70% flux loss due to phase errors. This allows us to restore more than 90% of a source flux. The method outlined in this work is not only applicable to the correction of deformation for other co-planar telescopes but also to single dish telescopes with deformation problems. This work also forms the basis of the upcoming science results of AMiBA-13.

The possibility that we live in a special place in the universe, close to the center of a large, radially inhomogeneous void, has attracted attention recently as an alternative to dark energy or modified gravity to explain the accelerating universe. We show that the distribution of orientations of galaxy pairs can be used to test the Copernican principle that we are not in a central or special region of Universe. The popular void models can not fit both the latest type Ia supernova, cosmic microwave background data and the distribution of orientations of galaxy pairs simultaneously. Our results rule out the void models at the $4\sigma$ confidence level as the origin of cosmic acceleration and favor the Copernican principle.

We present new quasars discovered in the vicinity of the Andromeda and Triangulum galaxies with the LAMOST during the 2010 and 2011 observational seasons. Quasar candidates are selected based on the available SDSS, KPNO 4 m telescope, XSTPS optical, and WISE near infrared photometric data. We present 509 new quasars discovered in a stripe of ~135 sq. deg from M31 to M33 along the Giant Stellar Stream in the 2011 pilot survey datasets, and also 17 new quasars discovered in an area of ~100 sq. deg that covers the central region and the southeastern halo of M31 in the 2010 commissioning datasets. These 526 new quasars have i magnitudes ranging from 15.5 to 20.0, redshifts from 0.1 to 3.2. They represent a significant increase of the number of identified quasars in the vicinity of M31 and M33. There are now 26, 62 and 139 known quasars in this region of the sky with i magnitudes brighter than 17.0, 17.5 and 18.0 respectively, of which 5, 20 and 75 are newly-discovered. These bright quasars provide an invaluable collection with which to probe the kinematics and chemistry of the ISM/IGM in the Local Group of galaxies. A total of 93 quasars are now known with locations within 2.5 deg of M31, of which 73 are newly discovered. Tens of quasars are now known to be located behind the Giant Stellar Stream, and hundreds behind the extended halo and its associated substructures of M31. The much enlarged sample of known quasars in the vicinity of M31 and M33 can potentially be utilized to construct a perfect astrometric reference frame to measure the minute PMs of M31 and M33, along with the PMs of substructures associated with the Local Group of galaxies. Those PMs are some of the most fundamental properties of the Local Group.

Bound electron contribution to the Doppler-shift of gamma-ray spectra in the positron-electron annihilation process of molecular methane has been studied in gas phase. Two accurate ab initio quantum mechanical schemes, i.e. the delocalized molecular orbital (MO) and the localized natural bond orbital (NBO) schemes, are applied to study the multi-centred methane molecule. The present ab initio calculations of methane indicate that the C-H bonds are polarized with the partial negative charge of -0.36 a.u. on the carbon atom and the partial positive charge of +0.09 a.u. on each of the hydrogen atoms. The positively charged hydrogen atoms produce repulsive Coulomb potentials to a positron. Both the MO and NBO schemes further reveal that the 2a1 electrons of methane, that is, the 2a1 electron component of the C-H bonds rather than the whole C-H bonds of methane, predominates the positron-electron annihilation gamma-ray spectra of the molecule. Electrons of a molecule which are dominant the positron-electron annihilation processes are called positrophilic electrons in the present study. It is further shown that the negative electrostatic potential (ESP) of methane facilitates with the density of the positrophilic 2a1 electrons of methane. Other valence electrons (e.g. 1t2) in the C-H bonds play a minor spectator role in the annihilation process of methane.

Theoretical gamma-ray spectra of molecule hexane have been calculated and compared with the experimental results in both gas (Surko et al, 1997) and liquid (Kerr et al, 1965) phases. The present study reveals that in gas phase not all valence electrons of hexane exhibit the same probability to annihilate a positron. Only the positrophilic electrons in the valence space dominate the gamma-ray spectra, which are in good agreement with the gas phase measurement. When hexane is confined in liquid phase, however, the intermolecular interactions ultimately eliminate the free molecular orientation and selectivity for the positrophilic electrons in the gas phase. As a result, the gamma-ray spectra of hexane become an averaged contribution from all valence electrons, which is again in agreement with liquid phase measurement. The roles of the positrophilic electrons in annihilation process for gas and liquid phases of hexane have been recognized for the first time in the present study.

We study the gravitational effect of non-self-annihilating dark matter on compact stellar objects. The self-interaction of condensate dark matter can give high accretion rate of dark matter onto stars. Phase transition to condensation state takes place when the dark matter density exceeds the critical value. A compact degenerate dark matter core is developed and alter the structure and stability of the stellar objects. Condensate dark matter admixed neutron stars is studied through the two-fuid TOV equation. The existence of condensate dark matter deforms the mass-radius relation of neutron stars and lower their maximum baryonic masses and radii. The possible effects on the Gamma-ray Burst rate in high redshift are discussed.

We explore high-redshift gamma-ray bursts (GRBs) as promising tools to probe pre-galactic metal enrichment. We utilize the bright afterglow of a Pop III GRB exploding in a primordial dwarf galaxy as a luminous background source, and calculate the strength of metal absorption lines that are imprinted by the first heavy elements in the intergalactic medium (IGM). To derive the GRB absorption line diagnostics, we use an existing highly-resolved simulation of the formation of a first galaxy which is characterized by the onset of atomic hydrogen cooling in a halo with virial temperature >10^4 K. We explore the unusual circumburst environment inside the systems that hosted Pop III stars, modeling the density evolution with the self-similar solution for a champagne flow. For minihalos close to the cooling threshold, the circumburst density is roughly proportional to (1+z) with values of about a few cm^{-3}. In more massive halos, corresponding to the first galaxies, the density may be larger, n>100 cm^{-3}. The resulting afterglow fluxes may be detectable with the JWST and VLA in the near-IR and radio wavebands, respectively, out to redshift z>20. We predict that the maximum of the afterglow emission shifts from near-IR to millimeter bands with peak fluxes from mJy to Jy at different observed times. GRBs are ideal tools for probing the metal enrichment in the early IGM, due to their high luminosities and featureless power-law spectra. The metals in the first galaxies produced by the first supernova (SN) explosions are likely to reside in low-ionization stages. We show that if the afterglow can be observed sufficiently early, analysis of the metal lines can distinguish whether the first heavy elements were produced in a PISN, or a core-collapse (Type II) SN, thus constraining the initial mass function of the first stars.

Gamma-ray bursts (GRBs) are attractive objects for constraining the nature of dark energy in a way complementary to other cosmological probes, especially at high redshifts. However, the apparent magnitude of distant GRBs can be distorted by the gravitational lensing from the density fluctuations along the line of sight. We investigate the gravitational lensing effect on the cosmological parameters and dark energy equation of state from GRBs. We first calibrated the GRB luminosity relations without assuming any cosmological models. The luminosity distances of low-redshift GRBs were calibrated with the cosmography method using a latest type Ia supernova (SNe Ia) sample. We find that the gravitational lensing has non-negligible effects on the determination of cosmological parameters and dark energy. The gravitational lensing shifts the best-fit constraints on cosmological parameters and dark energy. Because high-redshift GRBs are more likely to be reduced, the most probable value of the observed matter density $\Omega_M$ is slightly lower than its actual value. In the $\Lambda$CDM model, we find that the matter density parameter $\Omega_M$ will shift from 0.30 to 0.26 after including the gravitational lensing effect. The gravitational lensing also affects the dark energy equation of state by shifting it to a more negative value. We constrain the dark energy equation of state out to redshift $z\sim 8$ using GRBs for the first time, and find that the equation of state deviates from $\Lambda$CDM at the $1\sigma$ confidence level, but agrees with $w=-1$ at the $2\sigma$ confidence level.

Swift J1644+57 is an unusual transient event, likely powered by the tidal disruption of a star by a massive black hole. There are multiple short timescales X-ray flares were seen over a span of several days. We propose that these flares could be produced by internal shocks. In the internal shock model, the forward and reverse shocks are produced by collisions between relativistic shells ejected from central engine. The synchrotron emission from the forward and reverse shocks could dominate at two quite different energy bands under some conditions, the relativistic reverse shock dominates the X-ray emission and the Newtonian forward shock dominates the infrared and optical emission. We show that the spectral energy distribution of Swift J1644+57 could be explained by internal shock model.

We present results for the electron-impact excitation of highly-charged sulphur ions (S8+ - S11+) obtained using the intermediate-coupling frame transformation R-matrix approach. A detailed comparison of the target structure has been made for the four ions to assess the uncertainty on collision strengths from the target structure. Effective collision strengths (\Upsilon s) are presented at temperatures ranging from 2\times10^2(z+1)^2 K to 2\times10^6(z+1)^2 K (where z is the residual charge of ions. Detailed comparisons for the \Upsilon are made with the results of previous calculations for these ions, which will pose insight on the uncertainty in their usage by astrophysical and fusion modelling codes.

Several interesting luminosity correlations among gamma-ray burst (GRB) variables have been recently discussed extensively. In this paper, we derive the six luminosity correlations ($\tlag-L$, $V-L$, $\epkk-L$, $\epkk-E_\gamma$, $\trt-L$, $\epkk-E_{\gamma, \mathrm{iso}}$) from the light curves and spectra of the latest 116 long GRBs, including the time lag ($\tlag$) between low and high photon energy light curves, the variability ($V$) of the light curve, the peak energy of the spectrum ($\epkk$), and the minimum rise time ($\trt$) of the peaks. We find that the intrinsic scatter of the $V-L$ correlation is too large and there seems no inherent correlation between the two parameters using the latest GRB data. The other five correlations indeed exist when the sample is enlarged. The $\epkk-E_\gamma$ correlation has a significantly lower intrinsic scatter compared to the other correlations. We divide the full data into four redshift bins when testing possible evolution of the correlations with redshift. We find no statistically significant evidence for the redshift evolution of the luminosity correlations. To avoid the circularity problem when constraining the cosmological parameters, we simultaneously minimize $\chi^2$ with respect to both correlation parameters $a$, $b$ and the cosmological parameters using the maximum likelihood method. For the flat $\Lambda$CDM, the best fit is $\omm=0.31^{+0.13}_{-0.10}$. We also constrain the possible evolution of the equation of state (EOS) of the dark energy using the GRBs together with the Union2 compilation of SNe Ia and the $H(z)$ data. The result is consistent with the cosmological constant at $2 \sigma$ confidence level and mainly due to the GRB data, the dark energy EOS shows slight deviation from -1 at $z \geq0.5$ as was persistently presented with many previous data sets.

Recent studies suggest that Swift gamma-ray bursts (GRBs) may not trace an ordinary star formation history. Here we show that the GRB rate turns out to be consistent with the star formation history with an evolving stellar initial mass function (IMF). We first show that the latest Swift sample of GRBs reveals an increasing evolution in the GRB rate relative to the ordinary star formation rate at high redshifts. We then assume only massive stars with masses greater than the critical value to produce GRBs, and use an evolving stellar IMF suggested by Dav\'{e} (2010) to fit the latest GRB redshift distribution. This evolving IMF would increase the relative number of massive stars, which could lead to more GRB explosions at high redshifts. We find that the evolving IMF can well reproduce the observed redshift distribution of Swift GRBs.

Interferometric millimeter observations of the cosmic microwave background and clusters of galaxies with arcmin resolutions require antenna arrays with short spacings. Having all antennas co-mounted on a single steerable platform sets limits to the overall weight. A 25 kg lightweight novel carbon-fiber design for a 1.2 m diameter Cassegrain antenna is presented. The finite element analysis predicts excellent structural behavior under gravity, wind and thermal load. The primary and secondary mirror surfaces are aluminum coated with a thin TiO$_2$ top layer for protection. A low beam sidelobe level is achieved with a Gaussian feed illumination pattern with edge taper, designed based on feedhorn antenna simulations and verified in a far field beam pattern measurement. A shielding baffle reduces inter-antenna coupling to below $\sim$ -135 dB. The overall antenna efficiency, including a series of efficiency factors, is estimated to be around 60%, with major losses coming from the feed spillover and secondary blocking. With this new antenna, a detection rate of about 50 clusters per year is anticipated in a 13-element array operation.

Clusters of galaxies have been used extensively to determine cosmological parameters. A major difficulty in making best use of Sunyaev-Zel'dovich (SZ) and X-ray observations of clusters for cosmology is that using X-ray observations it is difficult to measure the temperature distribution and therefore determine the density distribution in individual clusters of galaxies out to the virial radius. Observations with the new generation of SZ instruments are a promising alternative approach. We use clusters of galaxies drawn from high-resolution adaptive mesh refinement (AMR) cosmological simulations to study how well we should be able to constrain the large-scale distribution of the intra-cluster gas (ICG) in individual massive relaxed clusters using AMiBA in its configuration with 13 1.2-m diameter dishes (AMiBA13) along with X-ray observations. We show that non-isothermal beta models provide a good description of the ICG in our simulated relaxed clusters. We use simulated X-ray observations to estimate the quality of constraints on the distribution of gas density, and simulated SZ visibilities (AMiBA13 observations) for constraints on the large-scale temperature distribution of the ICG. We find that AMiBA13 visibilities should constrain the scale radius of the temperature distribution to about 50% accuracy. We conclude that the upgraded AMiBA, AMiBA13, should be a powerful instrument to constrain the large-scale distribution of the ICG.

We investigate the contamination of the Sunyaev--Zel'dovich (SZ) effect for six galaxy clusters, A1689, A1995, A2142, A2163, A2261, and A2390, observed by the Y. T. Lee Array for Microwave Background Anisotropy during 2007. With the range of baselines used, we find that the largest effect (of order 13%-50% of the central SZ flux density) comes from primary anisotropies in the cosmic microwave background and exceeds the thermal noise in all six cases. Contamination from discrete radio sources is estimated to be at a level of (3%-60%) of the central SZ flux density. We use the statistics of these contaminating sources to estimate and correct the errors in the measured SZ effects of these clusters.

The Sunyaev-Zel'dovich Effect (SZE) has been observed toward six massive galaxy clusters, at redshifts 0.091 \leq z \leq 0.322 in the 86-102 GHz band with the Y. T. Lee Array for Microwave Background Anisotropy (AMiBA). We modify an iterative method, based on the isothermal \beta-models, to derive the electron temperature T_e, total mass M_t, gas mass M_g, and integrated Compton Y within r_2500, from the AMiBA SZE data. Non-isothermal universal temperature profile (UTP) \beta models are also considered in this paper. These results are in good agreement with those deduced from other observations. We also investigate the embedded scaling relations, due to the assumptions that have been made in the method we adopted, between these purely SZE-deduced T_e, M_t, M_g and Y. Our results suggest that cluster properties may be measurable with SZE observations alone. However, the assumptions built into the pure-SZE method bias the results of scaling relation estimations and need further study.

The high-redshift star formation rate (SFR) is difficult to measure directly even by modern approaches. Long-duration gamma-ray bursts (GRBs) can be detected to the edge of the visible universe because of their high luminorsities. The collapsar model of long gamma-ray bursts indicates that they may trace the star formation history. So long gamma-ray bursts may be a useful tool of measuring the high-redshift SFR. Observations show that long gamma-ray bursts prefer to form in a low-metallicity environment. We study the high-redshift SFR up to z~8.3 considering the Swift GRBs tracing the star formation history and the cosmic metallicity evolution in different background cosmological models including $\Lambda$CDM, quintessence, quintessence with a time-varying equation of state, and brane-world model. We use latest Swift GRBs including two highest-$z$ GRBs, GRB 080913 at $z=6.7$ and GRB 090423 at $z=8.3$. We find that the SFR at $z>4$ shows a steep decay with a slope of $\sim -5.0$ in $\Lambda$CDM. In the other three models, the high-redshift SFR is slightly different from $\Lambda$CDM model, and also shows a steep decay.

In this paper we investigate the deceleration, jerk and snap parameters to distinguish between the dark energy and modified gravity models by using high redshift gamma-ray bursts (GRBs) and supernovae (SNe). We first derive the expressions of deceleration, jerk and snap parameters in dark energy and modified gravity models. In order to constrain the cosmographic parameters, we calibrate the GRB luminosity relations without assuming any cosmological models using SNe Ia. Then we constrain the models (including dark energy and modified gravity models) parameters using type Ia supernovae and gamma-ray bursts. Finally we calculate the cosmographic parameters. GRBs can extend the redshift - distance relation up to high redshifts, because they can be detected to high redshifts.We find that the statefinder pair (r,s) could not be used to distinguish between some dark energy and modified gravity models, but these models could be differentiated by the snap parameter. Using the model-independent constraints on cosmographic parameters, we conclude that the \LambdaCDM model is consistent with the current data.

We investigate the scaling relations between the X-ray and the thermal Sunyaev-Zel'dovich Effect (SZE) properties of clusters of galaxies, using data taken during 2007 by the Y.T. Lee Array for Microwave Background Anisotropy (AMiBA) at 94 GHz for the six clusters A1689, A1995, A2142, A2163, A2261, and A2390. The scaling relations relate the integrated Compton-y parameter Y_{2500} to the X-ray derived gas temperature T_{e}, total mass M_{2500}, and bolometric luminosity L_X within r_{2500}. Our results for the power-law index and normalization are both consistent with the self-similar model and other studies in the literature except for the Y_{2500}-L_X relation, for which a physical explanation is given though further investigation may be still needed. Our results not only provide confidence for the AMiBA project but also support our understanding of galaxy clusters.

The Array for Microwave Background Anisotropy (AMiBA) is a radio interferometer for research in cosmology, currently operating 7 0.6m diameter antennas co-mounted on a 6m diameter platform driven by a hexapod mount. AMiBA is currently the largest hexapod telescope. We briefly summarize the hexapod operation with the current pointing error model. We then focus on the upcoming 13-element expansion with its potential difficulties and solutions. Photogrammetry measurements of the platform reveal deformations at a level which can affect the optical pointing and the receiver radio phase. In order to prepare for the 13-element upgrade, two optical telescopes are installed on the platform to correlate optical pointing tests. Being mounted on different locations, the residuals of the two sets of pointing errors show a characteristic phase and amplitude difference as a function of the platform deformation pattern. These results depend on the telescope's azimuth, elevation and polarization position. An analytical model for the deformation is derived in order to separate the local deformation induced error from the real hexapod pointing error. Similarly, we demonstrate that the deformation induced radio phase error can be reliably modeled and calibrated, which allows us to recover the ideal synthesized beam in amplitude and shape of up to 90% or more. The resulting array efficiency and its limits are discussed based on the derived errors.

Laboratory spectroscopy of non-thermal equilibrium plasmas photoionized by intense radiation is a key to understanding compact objects, such as black holes, based on astronomical observations. This paper describes an experiment to study photoionizing plasmas in laboratory under well-defined and genuine conditions. Photoionized plasma is here generated using a 0.5-keV Planckian x-ray source created by means of a laser-driven implosion. The measured x-ray spectrum from the photoionized silicon plasma resembles those observed from the binary stars Cygnus X-3 and Vela X-1 with the Chandra x-ray satellite. This demonstrates that an extreme radiation field was produced in the laboratory, however, the theoretical interpretation of the laboratory spectrum significantly contradicts the generally accepted explanations in x-ray astronomy. This model experiment offers a novel test bed for validation and verification of computational codes used in x-ray astronomy.

Motivated by the fact that both SNe Ia and GRBs seem to prefer a dark energy EOS greater than -1 at redshifts $z \gtrsim 0.5$, we perform a careful investigation on this situation. We find that the deviation of dark energy from the cosmological constant at redshifts $z \gtrsim 0.5$ is large enough that we should pay close attention to it with future observational data. Such a deviation may arise from some biasing systematic errors in the handling of SNe Ia and/or GRBs or more interestingly from the nature of the dark energy itself.

The Y. T. Lee Array for Microwave Background (AMiBA) has reported the first science results on the detection of galaxy clusters via the Sunyaev Zel'dovich effect. The science objectives required small reflectors in order to sample large scale structures (20') while interferometry provided modest resolutions (2'). With these constraints, we designed for the best sensitivity by utilizing the maximum possible continuum bandwidth matched to the atmospheric window at 86-102GHz, with dual polarizations. A novel wide-band analog correlator was designed that is easily expandable for more interferometer elements. MMIC technology was used throughout as much as possible in order to miniaturize the components and to enhance mass production. These designs will find application in other upcoming astronomy projects. AMiBA is now in operations since 2006, and we are in the process to expand the array from 7 to 13 elements.

The Y.T. Lee Array for Microwave Background Anisotropy (AMiBA) started scientific operation in early 2007. This work describes the optimization of the system performance for the measurements of the Sunyaev-Zel'dovich effect for six massive galaxy clusters at redshifts $0.09 - 0.32$. We achieved a point source sensitivity of $63\pm 7$ mJy with the seven 0.6m dishes in 1 hour of on-source integration in 2-patch differencing observations. We measured and compensated for the delays between the antennas of our platform-mounted interferometer. Beam switching was used to cancel instrumental instabilities and ground pick up. Total power and phase stability were good on time scales of hours, and the system was shown to integrate down on equivalent timescales of 300 hours per baseline/correlation, or about 10 hours for the entire array. While the broadband correlator leads to good sensitivity, the small number of lags in the correlator resulted in poorly measured bandpass response. We corrected for this by using external calibrators (Jupiter and Saturn). Using Jupiter as the flux standard, we measured the disk brightness temperature of Saturn to be $149^{+5}_{-12}$ K.

AMiBA is the largest hexapod astronomical telescope in current operation. We present a description of this novel hexapod mount with its main mechanical components -- the support cone, universal joints, jack screws, and platform -- and outline the control system with the pointing model and the operating modes that are supported. The AMiBA hexapod mount performance is verified based on optical pointing tests and platform photogrammetry measurements. The photogrammetry results show that the deformations in the inner part of the platform are less than 120 micron rms. This is negligible for optical pointing corrections, radio alignment and radio phase errors for the currently operational 7-element compact configuration. The optical pointing error in azimuth and elevation is successively reduced by a series of corrections to about 0.4 arcmin rms which meets our goal for the 7-element target specifications.

We present observations, analysis and results for the first-year operation of AMiBA, an interferometric experiment designed to study cosmology via the measurement of Cosmic Microwave Background (CMB). AMiBA is the first CMB interferometer operating at 3 mm to have reported successful results, currently with seven close-packed antennas of 60-cm diameter giving a synthesized resolution of around 6'. During 2007 AMiBA detected the Sunyaev-Zel'dovich effects (SZE) of six galaxy clusters at redshift 0.091 <= z <= 0.322. An observing strategy with on-off-source switching is used to minimize the effects from electronic offset and ground pickup. Planets were used to test the observational capability of AMiBA and to calibrate the conversion from correlator time-lag data to visibilities. The detailed formalism for data analysis is given. We summarize our early tests including observations of planets and quasars, and present images, visibility profiles, the estimated central coordinates, sizes, and SZE amplitudes of the galaxy clusters. Science results are summarized. We also discuss possible systematic effects in the results.

The Yuan-Tseh Lee Array for Microwave Background Anisotropy (AMiBA) is the first interferometer dedicated to studying the cosmic microwave background (CMB) radiation at 3mm wavelength. The choice of 3mm was made to minimize the contributions from foreground synchrotron radiation and Galactic dust emission. The initial configuration of seven 0.6m telescopes mounted on a 6-m hexapod platform was dedicated in October 2006 on Mauna Loa, Hawaii. Scientific operations began with the detection of a number of clusters of galaxies via the thermal Sunyaev-Zel'dovich effect. We compare our data with Subaru weak lensing data in order to study the structure of dark matter. We also compare our data with X-ray data in order to derive the Hubble constant.

We present a multiwavelength analysis of a sample of four hot (T_X>8keV) X-ray galaxy clusters (A1689, A2261, A2142, and A2390) using joint AMiBA Sunyaev-Zel'dovich effect (SZE) and Subaru weak lensing observations, combined with published X-ray temperatures, to examine the distribution of mass and the intracluster medium (ICM) in massive cluster environments. Our observations show that A2261 is very similar to A1689 in terms of lensing properties. Many tangential arcs are visible around A2261, with an effective Einstein radius \sim 40 arcsec (at z \sim 1.5), which when combined with our weak lensing measurements implies a mass profile well fitted by an NFW model with a high concentration c_{vir} \sim 10, similar to A1689 and to other massive clusters. The cluster A2142 shows complex mass substructure, and displays a shallower profile (c_{vir} \sim 5), consistent with detailed X-ray observations which imply recent interaction. The AMiBA map of A2142 exhibits an SZE feature associated with mass substructure lying ahead of the sharp north-west edge of the X-ray core suggesting a pressure increase in the ICM. For A2390 we obtain highly elliptical mass and ICM distributions at all radii, consistent with other X-ray and strong lensing work. Our cluster gas fraction measurements, free from the hydrostatic equilibrium assumption, are overall in good agreement with published X-ray and SZE observations, with the sample-averaged gas fraction of <f_{gas}> = 0.133 \pm 0.027, for our sample <M_{vir}> = (1.2 \pm 0.1) \times 10^{15} M_{sun} h^{-1}. When compared to the cosmic baryon fraction f_b = \Omega_b/\Omega_m constrained by the WMAP 5-year data, this indicates <f_{gas}>/f_b = 0.78 \pm 0.16, i.e., (22 \pm 16)% of the baryons are missing from the hot phase of clusters.

The behavior of the dark energy equation of state (EOS) is crucial in distinguishing different cosmological models. With a model independent approach, we constrain the possible evolution of the dark energy EOS. Gamma-ray bursts (GRBs) of redshifts up to $z>6$ are used, in addition to type Ia supernovae (SNe Ia). We separate the redshifts into 4 bins and assume a constant EOS parameter for dark energy in each bin. The EOS parameters are decorrelated by diagonalizing the covariance matrix. And the evolution of dark energy is estimated out of the uncorrelated EOS parameters. By including GRB luminosity data, we significantly reduce the confidence interval of the uncorrelated EOS parameter whose contribution mostly comes from the redshift bin of $0.5<z<1.8$. At high redshift where we only have GRBs, the constraints on the dark energy EOS are still very weak. However, we can see an obvious cut at about zero in the probability plot of the EOS parameter, from which we can infer that the ratio of dark energy to matter most probably continues to decrease beyond redshift 1.8. We carried out analyses with and without including the latest BAO measurements, which themselves favor a dark energy EOS of $w<-1$. If they are included, the results show some evidence of an evolving dark energy EOS. If not included, however, the results are consistent with the cosmological constant within $1 \sigma$ for redshift $0<z \lesssim 0.5$ and $2 \sigma$ for $0.5 \lesssim z<1.8$.

We describe methods used to validate data from the Y.T. Lee Array for Microwave Background Anisotropy (AMiBA), an interferometric array designed to measure the Sunyaev-Zel'dovich effect and the anisotropy of the Cosmic Microwave Background (CMB). We perform several statistical tests on data from pointed galaxy cluster observations taken in 2007 and noise data from long-term blank sky observations and measurements with the feeds covered by the absorbers. We apply power spectrum analysis, cross power spectrum analysis among different outputs with different time lags in our analog correlator, and sample variance law tests to noise data. We find that (1) there is no time variation of electronic offsets on the time scale of our two-patch observations (~10 minutes); (2) noise is correlated by less than 10% between different lags; and (3) the variance of noise scales with the inverse of time. To test the Gaussianity of the data, we apply Kolmogorov-Smirnov (K-S) tests to cluster data, and find that a 5% significance level efficiently detects data sets with known hardware problems without rejecting an excess of acceptable data. We also calculate third- and fourth-order moments and cumulants for the noise residual visibilities and find that about 95% of our data are within the 99% confidence regions of Gaussianity.

We explore the gamma-ray bursts' (GRBs') contributions in constraining the dark energy equation of state (EOS) at high ($1.8 < z < 7$) and at middle redshifts ($0.5 < z < 1.8$) and estimate how many GRBs are needed to get substantial constraints at high redshifts. We estimate the constraints with mock GRBs and mock type Ia supernovae (SNe Ia) for comparisons. When constraining the dark energy EOS in a certain redshift range, we allow the dark energy EOS parameter to vary only in that redshift bin and fix EOS parameters elsewhere to -1. We find that it is difficult to constrain the dark energy EOS beyond the redshifts of SNe Ia with GRBs unless some new luminosity relations for GRBs with smaller scatters are discovered. However, at middle redshifts, GRBs have comparable contributions with SNe Ia in constraining the dark energy EOS.

The property of dark energy and the physical reason for acceleration of the present universe are two of the most difficult problems in modern cosmology. The dark energy contributes about two-thirds of the critical density of the present universe from the observations of type-Ia supernova (SNe Ia) and anisotropy of cosmic microwave background (CMB).The SN Ia observations also suggest that the universe expanded from a deceleration to an acceleration phase at some redshift, implying the existence of a nearly uniform component of dark energy with negative pressure. We use the ``gold'' sample containing 157 SNe Ia and two recent well-measured additions, SNe Ia 1994ae and 1998aq to explore the properties of dark energy and the transition redshift. For a flat universe with the cosmological constant, we measure $\Omega_{M}=0.28_{-0.05}^{+0.04}$, which is consistent with Riess et al. The transition redshift is $z_{T}=0.60_{-0.08}^{+0.06}$. We also discuss several dark energy models that define the $w(z)$ of the parameterized equation of state of dark energy including one parameter and two parameters ($w(z)$ being the ratio of the pressure to energy density). Our calculations show that the accurately calculated transition redshift varies from $z_{T}=0.29_{-0.06}^{+0.07}$ to $z_{T}=0.60_{-0.08}^{+0.06}$ across these models. We also calculate the minimum redshift $z_{c}$ at which the current observations need the universe to accelerate.

It has been widely shown that the cosmological parameters and dark energy can be constrained by using data from type-Ia supernovae (SNe Ia), the cosmic microwave background (CMB) anisotropy, the baryon acoustic oscillation (BAO) peak from Sloan Digital Sky Survey (SDSS), the X-ray gas mass fraction in clusters, and the linear growth rate of perturbations at z=0.15 as obtained from the 2dF Galaxy Redshift Survey. Recently, gamma-ray bursts (GRBs) have also been argued to be promising standard candles for cosmography. In this paper, we present constraints on the cosmological parameters and dark energy by combining a recent GRB sample including 69 events with the other cosmological probes. First, we find that for the LambdaCDM cosmology this combination makes the constraints stringent and the best fit is close to the flat universe. Second, we fit the flat Cardassian expansion model and find that this model is consistent with the LambdaCDM cosmology. Third, we present constraints on several two-parameter dark energy models and find that these models are also consistent with the LambdaCDM cosmology. Finally, we reconstruct the dark energy equation-of-state parameter w(z) and the deceleration parameter q(z). We see that the acceleration could have started at a redshift from z_T=0.40_{-0.08}^{+0.14} to z_T=0.65_{-0.05}^{+0.10}. This difference in the transition redshift is due to different dark energy models that we adopt. The most stringent constraint on w(z) lies in the redshift range z\sim 0.3-0.6.

We split the two-Higgs-doublet model by assuming very different vevs for the two doublets: the vev is at weak scale (174 GeV) for the doublet \Phi_1 and at neutrino-mass scale (10^{-2} \sim 10^{-3} eV) for the doublet \Phi_2. \Phi_1 is responsible for giving masses to all fermions except neutrinos; while \Phi_2 is responsible for giving neutrino masses through its tiny vev without introducing see-saw mechanism. Among the predicted five physical scalars H, h, A^0 and H^{\pm}, the CP-even scalar h is as light as 10^{-2} \sim 10^{-3}eV while others are at weak scale. We identify h as the cosmic dark energy field and the other CP-even scalar H as the Standard Model Higgs boson; while the CP-odd A^0 and the charged H^{\pm} are the exotic scalars to be discovered at future colliders. Also we demonstrate a possible dynamical origin for the doublet \Phi_2 from neutrino condensation caused by some unknown dynamics.

A new method of measuring cosmology with gamma-ray bursts(GRBs) has been proposed by Liang and Zhang recently. In this method, only observable quantities including the rest frame peak energy of the \nu F_\nu spectrum (E'_p), the isotropic energy of GRB (E_{\gamma,iso}), and the break time of the optical afterglow light curves in the rest frame (t'_b) are used. By considering this method we constrain the cosmological parameters and the redshift at which the universe expanded from the deceleration to acceleration phase. We add five recently-detected GRBs to the sample and derive E_{\gamma, iso}/10^{52} ergs=(0.93\pm 0.25)\times (E'_p/100 keV)^{1.91\pm 0.32}\times (t'_b/day)^{-0.93\pm 0.38} for a flat cosmology with \Omega_M=0.28 and H_0=71.0 km s^{-1} Mpc^{-1}. This relation is independent of the medium density around bursts and the efficiency of conversion of the explosion energy to gamma-ray energy that are difficult to measure. We regard the E_{\gamma,iso}(E'_p, t'_b) relationship as a standard candle and find 0.05<\Omega_M<0.48 and \Omega_\Lambda<1.15 (at the 1\sigma confidence level). In a flat universe with the cosmological constant we obtain 0.25<\Omega_M <0.46 and 0.54<\Omega_\Lambda<0.78 at the 1\sigma confidence level. The transition redshift is z_T=0.69_{-0.12}^{+0.11}. Combining 20 GRBs and 157 type Ia supernovae, we find \Omega_M=0.29 \pm0.03 for a flat universe and the transition redshift is z_T=0.61_{-0.05}^{+0.06}, which is slightly larger than the value found by considering SNe Ia alone. In particular, We also discuss several dark-energy models in which the equation of state w(z) is arameterized and investigate constraints on the cosmological parameters in detail.

In the Affleck-Dine mechanism of baryogenesis, non-topological solitons called Q-balls can be formed. In this work we study the Q-balls which decay during the Big Bang Nucleosythesis (BBN) era. We show that such late-decaying baryonic Q-balls with lifetime of about 10^3 s can provide a new developing mechanism for the BBN through a rolling baryon to photon ratio \eta, which can naturally explain the discrepancy of BBN prediction with the WMAP data on ^7Li abundance. For the late-decaying leptonic Q-balls with lifetime of about 10^6 s, we find that their decay product, gravitino, can serve as a dark matter candidate without overclosing the universe and, at the same time, can give an explanation of the approximate equality of dark and baryon matter densities.

In this work we perform some studies related to dark energy. Firstly, we propose a dynamical approach to explain the dark energy content of the universe. We assume that a massless scalar field couples to the Hubble parameter with some Planck-mass suppressed interactions. This scalar field developes a Hubble parameter-dependent (thus time-dependent) vacum expectation value, which renders a time-independent relative density for dark energy and thus can explain the coincidence of the dark energy density of the universe. Secondly, we assume the dark matter particle is meta-stable and decays very lately into the dark energy scalar field. Such a conversion of matter to dark energy can give an explanation for the starting time of the accelerating expansion of the universe. Thirdly, we introduce multiple Affleck-Dine fields to the landscape scenario of dark energy in order to have the required baryon-asymmetry universe.

The heavy gravitinos in minimal supergravity (mSUGRA) models couple gravitionally and can naturally be the Superweakly Interacting Massive Particle (SuperWIMP). As plausible candidates for the cold dark matter in the universe, such gravitino SuperWIMPs can naturally evade the current dark matter detection experiments due to their superweak couplings. However, this scenario must subject to the constraints from the Big Bang nucleosynthesis (BBN) predictions for light element abundances as well as the Wilkinson Mircrowave Anisotropy Probe (WMAP) data for the relic density. Our study shows that requiring BBN predictions to agree with the WMAP data of light element abundances can impose upper and lower mass bounds for both the gravitino and the next-to-lightest supersymmetric particle (NLSP). From a scan over the mSUGRA parameter space subjected to the BBN constraints and the WMAP data, we find that in the neutralino NLSP case, the low tan-beta (< 40) region as well as the region accessible at CERN Large Hadron Colldier (LHC) will be severely restrained.

A model is proposed for the observed combination of power-law and thermal emission of keV X-rays from rotationally powered pulsars. For gamma-ray pulsars with accelerators very many stellar radii above the neutron star surface, 100 MeV curvature gamma-rays from $e^{-}$ or $e^{+}$ flowing starward out of such accelerators are converted to electron-positron pairs on closed field lines all around the star. These pairs strongly affect X-ray emission from near the star in two ways. (1) The pairs are a source of synchrotron emission immediately following their creation in regions where $B \sim 10^{10}$ G. This emission, in the photon energy range 0.1 keV $\lesssim E_{X} \lesssim$ 5 MeV, has a power-law spectrum with energy index 0.5 and X-ray luminosity that depends on the backflow current, and is typically $\sim 10^{33}$ \lum. (2) The pairs ultimately form a cyclotron resonance ``blanket'' surrounding the star except for two holes along the open field line bundles which pass through it. In such a blanket the gravitational pull on electron-positron pairs toward the star is balanced by the hugely amplified push of outflowing surface emitted X-rays wherever cyclotron resonance occurs. Because of it the neutron star is surrounded by a leaky ``hohlraum'' of hot blackbody radiation with two small holes, which prevents direct X-ray observation of a heated polar cap of a gamma-ray pulsar. Weakly spin-modulated radiation from the blanket together with more strongly spin-modulated radiation from the holes through it would then dominate observed low energy (0.1--10 keV) emission. For non-gamma-ray pulsars, in which no such accelerators with their accompanying extreme relativistic backflow toward the star are ...