The Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) program monitors 849 active galactic nuclei (AGN) both spectroscopically and photometrically. The photometric observations used in this work span over four years and provide an excellent baseline for variability studies of these objects. We present the photometric light curves from 2014 to 2017 obtained by the Steward Observatory's Bok telescope and the CFHT telescope with MegaCam. We provide details on the data acquisition and processing of the data from each telescope, the difference imaging photometry used to produce the light curves, and the calculation of a variability index to quantify each AGN's variability. We find that the Welch-Stetson J-index provides a useful characterization of AGN variability and can be used to select AGNs for further study.

Accurate background estimation is essential for spectral and temporal analysis in astrophysics. In this work, we construct the in-orbit background model for the High-Energy Telescope (HE) of the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT). Based on the two-year blank sky observations of Insight-HXMT/HE, we first investigate the basic properties of the background and find that both the background spectral shape and intensity have long-term evolution at different geographical sites. The entire earth globe is then divided into small grids, each with a typical area of 5x5 square degrees in geographical coordinate system. For each grid, an empirical function is used to describe the long-term evolution of each channel of the background spectrum; the intensity of the background can be variable and a modification factor is introduced to account for this variability by measuring the contemporary flux of the blind detector. For a given pointing observation, the background model is accomplished by integrating over the grids that are passed by the track of the satellite in each orbit. Such a background model is tested with both the blank sky observations and campaigns for observations of a series of celestial sources. The results show an average systematic error of 1.5% for the background energy spectrum (26-100 keV) under a typical exposure of 8 ks, and <3% for background light curve estimation (30-150 keV). Therefore, the background model introduced in this paper is included in the Insight-HXMT software as a standard part specialized for both spectral and temporal analyses.

Fast radio bursts (FRBs) are short pulses observed in radio band from cosmological distances, some of which emit repeating bursts. The physical origins of these mysterious events have been subject to wide speculations and heated debates. One class of models invoke soft gamma-ray repeaters (SGRs), or magnetars, as the sources of FRBs. Magnetars are rotating neutron stars with extremely strong magnetic field and can sporadically emit bursts from X-ray ($\sim$keV) to soft gamma-ray ($\sim$sub-MeV) with duration from $10^{-2}$~s to $10^2$~s. However, even though some bright radio bursts have been observed from some magnetars, no FRB-like events had been detected to be associated with any magnetar burst, including one giant flare, and no radio burst has been associated with any X-ray event from any magnetar. Therefore, there is still no observational evidence for magnetar-FRB association up to today. Recently, a pair of FRB-like bursts (FRB~200428 hereafter) separated by 30 milliseconds (ms) were detected from the general direction of the Galactic magnetar SGR~J1935+2154. Here we report the detection of a non-thermal X-ray burst in the 1-250 keV energy band with the Insight-HXMT satellite, which we identify as emitted from SGR~J1935+2154. The burst showed two hard peaks with a separation of $\sim30$ ms, consistent with the separation between the two bursts in FRB~200428. The delay time between the double radio and X-ray peaks is $\sim8.57$ s, fully consistent with the dispersion delay of FRB~200428. We thus identify the non-thermal X-ray burst is associated with FRB~200428 whose high energy counterpart is the two hard peaks in X-ray. Our results suggest that the non-thermal X-ray burst and FRB~200428 share the same physical origin in an explosive event from SGR~J1935+2154.

We use the In data collected during the 2019 outburst from X-ray pulsar 4U 1901+03 to complement the orbital parameters reported by Fermi/GBM. Using the Insight-HXMT, we examine the correlation between the derivative of the intrinsic spin frequency and bolometric flux based on accretion torque models. It was found that the pulse profiles significantly evolve during the outburst. The existence of two types of the profile's pattern discovered in the Insight-HXMT data indicates that this source experienced transition between a super-critical and a sub-critical accretion regime during its 2019 outburst. Based on the evolution of the pulse profiles and the torque model, we derive the distance to 4U 1901+03 as 12.4+-0.2 kpc.

We study the evolution of the temporal properties of MAXI 1820+070 during the 2018 outburst in its hard state from MJD 58190 to 58289 with Insight-HXMT in a broad energy band 1-150 keV. We find different behaviors of the hardness ratio, the fractional rms and time lag before and after MJD 58257, suggesting a transition occurred at around this point. The observed time lags between the soft photons in the 1-5 keV band and the hard photons in higher energy bands, up to 150 keV, are frequency-dependent: the time lags in the low-frequency range, 2-10 mHz, are both soft and hard lags with a timescale of dozens of seconds but without a clear trend along the outburst; the time lags in the high-frequency range, 1-10 Hz, are only hard lags with a timescale of tens of milliseconds; first increase until around MJD 58257 and decrease after this date. The high-frequency time lags are significantly correlated to the photon index derived from the fit to the quasi-simultaneous NICER spectrum in the 1-10 keV band. This result is qualitatively consistent with a model in which the high-frequency time lags are produced by Comptonization in a jet.

With more than 150 blank sky observations at high Galactic latitude, we make a systematic study to the background of the Low Energy Telescope (LE) of the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT). Both the on-ground simulation and the in-orbit observation indicate that the background spectrum mainly has two components. One is the particle background that dominates above 7 keV and its spectral shape is consistent in every geographical locations. Another is the diffuse X-ray background that dominates below 7 keV and has a stable spectrum less dependent of the sky region. The particle background spectral shape can be obtained from the blind detector data of all the blank sky observations, and the particle background intensity can be measured by the blind detector at 10-12.5 keV. The diffuse X-ray background in the high Galactic latitude can also be obtained from the blank sky spectra after subtracting the particle background. Based on these characteristics, we develop the background model for both the spectrum and the light curve. The systematic error for the background spectrum is investigated with different exposures (T_exp). For the spectrum with T_exp=1 ks, the average systematic errors in 1-7 keV and 1-10 keV are 4.2% and 3.7%, respectively. We also perform the systematic error analyses of the background light curves with different energy bands and time bins. The results show that the systematic errors for the light curves with different time bins are <8% in 1-10 keV.

Glitches correspond to sudden jumps of rotation frequency ($\nu$) and its derivative ($\dot{\nu}$) of pulsars, the origin of which remains not well understood yet, partly because the jump processes of most glitches are not well time-resolved. There are three large glitches of the Crab pulsar, detected in 1989, 1996 and 2017, which were found to have delayed spin-up processes before the normal recovery processes. Here we report two additional glitches of the Crab pulsar occurred in 2004 and 2011 for which we discovered delayed spin up processes, and present refined parameters of the largest glitch occurred in 2017. The initial rising time of the glitch is determined as $<0.48$ hour. We also carried out a statistical study of these five glitches with observed spin-up processes. The two glitches occurred in 2004 and 2011 have delayed spin-up time scales ($\tau_{1}$) of $1.7\pm0.8$\,days and $1.6\pm0.4$\,days, respectively. We find that the $\Delta{\nu}$ vs. $|\Delta{\dot\nu}|$ relation of these five glitches is similar to those with no detected delayed spin-up process, indicating that they are similar to the others in nature except that they have larger amplitudes. For these five glitches, the amplitudes of the delayed spin-up process ($|\Delta{\nu}_{\rm d1}|$) and recovery process ($\Delta{\nu}_{\rm d2}$), their time scales ($\tau_{1}$, $\tau_{2}$), and permanent changes in spin frequency ($\Delta{\nu}_{\rm p}$) and total frequency step ($\Delta{\nu}_{\rm g}$) have positive correlations. From these correlations, we suggest that the delayed spin-up processes are common for all glitches, but are too short and thus difficult to be detected for most glitches.

No robust detection of prompt electromagnetic counterparts to fast radio bursts (FRBs) has yet been obtained, in spite of several multi-wavelength searches carried out so far. Specifically, X/gamma-ray counterparts are predicted by some models. We planned on searching for prompt gamma-ray counterparts in the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) data, taking advantage of the unique combination of large effective area in the keV-MeV energy range and of sub-ms time resolution. We selected 39 FRBs that were promptly visible from the High-Energy (HE) instrument aboard Insight-HXMT. After calculating the expected arrival times at the location of the spacecraft, we searched for a significant excess in both individual and cumulative time profiles over a wide range of time resolutions, from several seconds down to sub-ms scales. Using the dispersion measures in excess of the Galactic terms, we estimated the upper limits on the redshifts. No convincing signal was found and for each FRB we constrained the gamma-ray isotropic-equivalent luminosity and the released energy as a function of emission timescale. For the nearest FRB source, the periodic repeater FRB180916.J0158+65, we find $L_{\gamma,iso}<5.5\times 10^{47}$ erg/s over 1 s, whereas $L_{\gamma,iso}<10^{49}-10^{51}$ erg/s for the bulk of FRBs. The same values scale up by a factor of ~100 for a ms-long emission. Even on a timescale comparable with that of the radio pulse itself no keV-MeV emission is observed. A systematic association with either long or short GRBs is ruled out with high confidence, except for subluminous events, as is the case for core-collapse of massive stars (long) or binary neutron star mergers (short) viewed off axis. Only giant flares from extra-galactic magnetars at least ten times more energetic than Galactic siblings are ruled out for the nearest FRB.

We present the calibration of the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) X-ray satellite, which can be used to perform timing and spectral studies of bright X-ray sources. Insight-HXMT carries three main payloads onboard: the High Energy X-ray telescope (HE), the Medium Energy X-ray telescope (ME) and the Low Energy X-ray telescope (LE). In orbit, the radioactive sources, activated lines, the fluorescence lines and celestial sources are used to calibrate the energy scale and energy resolution of the payloads. The Crab nebular is adopted as the primary effective area calibrator and empirical functions are constructed to modify the simulated effective areas of the three payloads respectively. The systematic errors of HE, compared to the model of the Crab nebular, are less than 2% in 28--120 keV and 2%--10% above 120 keV. The systematic errors of ME are less than 1.5% in 10--35 keV. The systematic errors of LE are less than 1% in 1--7 keV except the Si K--edge (1.839 keV, up to 1.5%) and less than 2% in 7--10 keV.

We report on our analysis of the 2019 outburst of the X-ray accreting pulsar 4U 1901+03 observed with Insight-HXMT and NICER. Both spectra and pulse profiles evolve significantly in the decaying phase of the outburst. Dozens of flares are observed throughout the outburst. They are more frequent and brighter at the outburst peak. We find that the flares, which have a duration from tens to hundreds of seconds, are generally brighter than the persistent emission by a factor of $\sim$ 1.5. The pulse profile shape during the flares can be significantly different than that of the persistent emission. In particular, a phase shift is clearly observed in many cases. We interpret these findings as direct evidence of changes of the pulsed beam pattern, due to transitions between the sub- and super-critical accretion regimes on a short time scale. We also observe that at comparable luminosities the flares' pulse profiles are rather similar to those of the persistent emission. This indicates that the accretion on the polar cap of the neutron star is mainly determined by the luminosity, i.e., the mass accretion rate.

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

A new black hole X-ray binary (BHXRB) MAXI J1535-571 was discovered by MAXI during its outburst in 2017. Using observations taken by the first Chinese X-ray satellite, the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT), we perform a joint spectral analysis (2-150 keV) in both energy and time domains. The energy spectra provide the essential input for probing the intrinsic Quasi-Periodic Oscillation (QPO) fractional rms spectra (FRS). Our results show that during the intermediate state, the energy spectra are in general consistent with those reported by Swift/XRT and NuSTAR. However, the QPO FRS become harder and the FRS residuals may suggest the presence of either an additional power-law component in the energy spectrum or a turn-over in the intrinsic QPO FRS at high energies.

We study the 2018 outburst of Aql X-1 via the monitor of all sky X-ray image (MAXI) data. We show that the outburst starting in February 2018 is a member of short-low class in the frame of outburst duration and the peak count rate although the outburst morphology is slightly different from the other fast-rise-exponential-decay (FRED) type outbursts with a milder rising stage. We study the partial accretion in the weak propeller stage of Aql X-1 via the MAXI data of the 2018 outburst. We report on the spectral analysis of 3 observations of Aquila X-1 obtained by Insight - hard X-ray modulation telescope (Insight-HXMT) during the late decay stage of the 2018 outburst. We discuss that the data taken by Insight-HXMT is just after the transition to the weak propeller stage. Our analysis shows the necessity of a comptonization component to take into account the existence of an electron cloud resulting photons partly up-scattered.

As China's first X-ray astronomical satellite, the Hard X-ray Modulation Telescope (HXMT), which was dubbed as Insight-HXMT after the launch on June 15, 2017, is a wide-band (1-250 keV) slat-collimator-based X-ray astronomy satellite with the capability of all-sky monitoring in 0.2-3 MeV. It was designed to perform pointing, scanning and gamma-ray burst (GRB) observations and, based on the Direct Demodulation Method (DDM), the image of the scanned sky region can be reconstructed. Here we give an overview of the mission and its progresses, including payload, core sciences, ground calibration/facility, ground segment, data archive, software, in-orbit performance, calibration, background model, observations and some preliminary results.

We present a detailed timing study of the brightest persistent X-ray source Sco X-1 using the data collected by the Hard X-ray Modulation Telescope ($Insight$-HXMT) from July 2017 to August 2018. A complete $Z$-track hardness-intensity diagram (HID) is obtained. The normal branch oscillations (NBOs) at $\sim$ 6 Hz in the lower part of the normal branch (NB) and the flare branch oscillations (FBOs) at $\sim$ 16 Hz in the beginning part of the flaring branch (FB) are found in observations with the Low Energy X-ray Telescope (LE) and the Medium Energy X-ray Telescope (ME) of $Insight$-HXMT, while the horizontal branch oscillations (HBOs) at $\sim$ 40 Hz and the kilohertz quasi-periodic oscillations (kHz QPOs) at $\sim$ 800 Hz are found simultaneously up to 60 keV for the first time on the horizontal branch (HB) by the High Energy X-ray Telescope (HE) and ME. We find that for all types of the observed QPOs, the centroid frequencies are independent of energy, while the root mean square (rms) increases with energy; the centroid frequencies of both the HBOs and kHz QPOs increase along the $Z$-track from the top to the bottom of the HB; and the NBOs show soft phase lags increasing with energy. A continuous QPO transition from the FB to NB in $\sim$ 200 s are also detected. Our results indicate that the non-thermal emission is the origin of all types of QPOs, the innermost region of the accretion disk is non-thermal in nature, and the corona is nonhomogeneous geometrically.

The evidences for the influence of thermonuclear (type-I) X-ray bursts upon the surrounding environments in neutron star low-mass X-ray binaries (LMXB) were detected previously via spectral and timing analyses. Benefitting from a broad energy coverage of Insight-HXMT, we analyze one photospheric radius expansion (PRE) burst, and find an emission excess at soft X-rays. Our spectral analysis shows that, such an excess is not likely relevant to the disk reflection induced by the burst emission and can be attributed to an enhanced pre-burst/persistent emission. We find that the burst and enhanced persistent emissions sum up to exceed Eddington luminosity by $\sim$ 40 percentages. We speculate that the enhanced emission is from a region beyond the PRE radius, or through the Comptonization of the corona.

Corona cooling was detected previously from stacking a series of short type-I bursts occurred during the low/had state of atoll outburst. Type-I bursts are hence regarded as sharp probe to our better understanding on the basic property of the corona. The launch of the first Chinese X-ray satellite Insight-HXMT has large detection area at hard X-rays which provide almost unique chance to move further in this research field. We report the first detection of the corona cooling by Insight-HXMT from single short type-I burst showing up during {\bf flare} of 4U 1636-536. This type-I X-ray burst has a duration of $\sim$13 seconds and hard X-ray shortage is detected with significance 6.2~$\sigma$ in 40-70 keV. A cross-correlation analysis between the lightcurves of soft and hard X-ray band, shows that the corona shortage lag the burst emission by 1.6 $\pm$1.2~s. These results are consistent with those derived previously from stacking a large amount of bursts detected by RXTE/PCA within a series of {\bf flares} of 4U 1636-536. Moreover, the broad bandwidth of Insight-HXMT allows as well for the first time to infer the burst influence upon the continuum spectrum via performing the spectral fitting of the burst, which ends up with the finding that hard X-ray shortage appears at around 40 keV in the continuum spectrum. These results suggest that the evolution of the corona along with the outburst{\bf /flare} of NS XRB may be traced via looking into a series of embedded type-I bursts by using Insight-HXMT.

The Hard X-ray Modulation Telescope HXMT or also dubbed as Insight HXMT is China s first astronomical satellite . It was launched on 15 th June 2017 in JiuQuan, China and is currently in service smoothly. It was designed to perform point ing , scan ning and gamma ray burst (GRB) observations and , based on the Direct Demodulation Method (DDM), the image of the scanned sky region can be reconstructed. Here we introduce the mission and its progresses in aspects of payload, core sciences, ground calibration/facility , ground segment, data archive, software, in orbit performance, calibration, background model, observations and preliminary results .

The Hard X-ray Modulation Telescope (HXMT) named Insight is China's first X-ray astronomical satellite. The Low Energy X-ray Telescope (LE) is one of its main payloads onboard. The detectors of LE adopt swept charge device CCD236 with L-shaped transfer electrodes. Charges in detection area are read out continuously along specific paths, which leads to a time response distribution of photons readout time. We designed a long exposure readout mode to measure the time response distribution. In this mode, CCD236 firstly performs exposure without readout, then all charges generated in preceding exposure phase are read out completely. Through analysis of the photons readout time in this mode, we obtained the probability distribution of photons readout time.

We present a study of timing properties of the accreting pulsar 2S 1417-624 observed during its 2018 outburst, based on Swift/BAT, Fermi/GBM, Insight-HXMT and NICER observations. We report a dramatic change of the pulse profiles with luminosity. The morphology of the profile in the range 0.2-10.0keV switches from double to triple peaks at $\sim2.5$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$ and from triple to quadruple peaks at $\sim7$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$. The profile at high energies (25-100keV) shows significant evolutions as well. We explain this phenomenon according to existing theoretical models. We argue that the first change is related to the transition from the sub to the super-critical accretion regime, while the second to the transition of the accretion disc from the gas-dominated to the radiation pressure-dominated state. Considering the spin-up as well due to the accretion torque, this interpretation allows to estimate the magnetic field self-consistently at $\sim7\times 10^{12}$G.

The long-term evolution of the centroid energy of the CRSF in Her X-1 is still a mystery. We report a new measurement from a campaign between {\sl Insight}-HXMT and {\sl NuSTAR} performed in February 2018. Generally, the two satellites show well consistent results of timing and spectral properties. The joint spectral analysis confirms that the previously observed long decay phase has ended, and that the line energy instead keeps constant around 37.5 keV after flux correction.

We report on analysis of observations of the bright transient X-ray pulsar \src obtained during its 2017-2018 giant outburst with Insight-HXMT, \emph{NuSTAR}, and \textit{Swift} observatories. We focus on the discovery of a sharp state transition of the timing and spectral properties of the source at super-Eddington accretion rates, which we associate with the transition of the accretion disk to a radiation pressure dominated (RPD) state, the first ever directly observed for magnetized neutron star. This transition occurs at slightly higher luminosity compared to already reported transition of the source from sub- to super-critical accretion regime associate with onset of an accretion column. We argue that this scenario can only be realized for comparatively weakly magnetized neutron star, not dissimilar to other ultra-luminous X-ray pulsars (ULPs), which accrete at similar rates. Further evidence for this conclusion is provided by the non-detection of the transition to the propeller state in quiescence which strongly implies compact magnetosphere and thus rules out magnetar-like fields.

In this work, we report the in-orbit demonstration of X-ray pulsar navigation with Insight-Hard X-ray Modulation Telescope (Insight-HXMT), which was launched on Jun. 15th, 2017. The new pulsar navigation method 'Significance Enhancement of Pulse-profile with Orbit-dynamics' (SEPO) is adopted to determine the orbit with observations of only one pulsar. In this test, the Crab pulsar is chosen and observed by Insight-HXMT from Aug. 31th to Sept. 5th in 2017. Using the 5-day-long observation data, the orbit of Insight-HXMT is determined successfully with the three telescopes onboard - High Energy X-ray Telescope (HE), Medium Energy X-ray Telescope (ME) and Low Energy X-ray Telescope (LE) - respectively. Combining all the data, the position and velocity of the Insight-HXMT are pinpointed to within 10 km (3 sigma) and 10 m/s (3 sigma), respectively.

We investigate the effects of extended multi-year light curves (9-year photometry and 5-year spectroscopy) on the detection of time lags between the continuum variability and broad-line response of quasars at z>~1.5, and compare with the results using 4-year photometry+spectroscopy presented in a companion paper. We demonstrate the benefits of the extended light curves in three cases: (1) lags that are too long to be detected by the shorter-duration data but can be detected with the extended data; (2) lags that are recovered by the extended light curves but missed in the shorter-duration data due to insufficient light curve quality; and (3) lags for different broad line species in the same object. These examples demonstrate the importance of long-term monitoring for reverberation mapping to detect lags for luminous quasars at high-redshift, and the expected performance of the final dataset from the Sloan Digital Sky Survey Reverberation Mapping project that will have 11-year photometric and 7-year spectroscopic baselines.

The recently discovered neutron star transient Swift J0243.6+6124 has been monitored by {\it the Hard X-ray Modulation Telescope} ({\it Insight-\rm HXMT). Based on the obtained data, we investigate the broadband spectrum of the source throughout the outburst. We estimate the broadband flux of the source and search for possible cyclotron line in the broadband spectrum. No evidence of line-like features is, however, found up to $\rm 150~keV$. In the absence of any cyclotron line in its energy spectrum, we estimate the magnetic field of the source based on the observed spin evolution of the neutron star by applying two accretion torque models. In both cases, we get consistent results with $B\rm \sim 10^{13}~G$, $D\rm \sim 6~kpc$ and peak luminosity of $\rm >10^{39}~erg~s^{-1}$ which makes the source the first Galactic ultraluminous X-ray source hosting a neutron star.

We present reverberation-mapping lags and black-hole mass measurements using the CIV 1549 broad emission line from a sample of 349 quasars monitored as a part of the Sloan Digital Sky Survey Reverberation Mapping Project. Our data span four years of spectroscopic and photometric monitoring for a total baseline of 1300 days. We report significant time delays between the continuum and the CIV 1549 emission line in 52 quasars, with an estimated false-positive detection rate of 10%. Our analysis of marginal lag measurements indicates that there are on the order of 100 additional lags that should be recoverable by adding more years of data from the program. We use our measurements to calculate black-hole masses and fit an updated CIV radius-luminosity relationship. Our results significantly increase the sample of quasars with CIV RM results, with the quasars spanning two orders of magnitude in luminosity toward the high-luminosity end of the CIV radius-luminosity relation. In addition, these quasars are located at among the highest redshifts (z~1.4-2.8) of quasars with black hole masses measured with reverberation mapping. This work constitutes the first large sample of CIV reverberation-mapping measurements in more than a dozen quasars, demonstrating the utility of multi-object reverberation mapping campaigns.

The DESI Legacy Imaging Surveys are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image approximately 14,000 deg^2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12 and 22 micorons) observed by the Wide-field Infrared Survey Explorer (WISE) satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.

We present long-term, multi-color photometric monitoring and spectroscopic observations of the short-period M-type eclipsing binary BX Tri. Six flare events were recorded in 4 years from 2014 to 2017. Three of them were detected on one night within an orbital cycle. The strongest one was identified on December 23, 2014. With the amplitudes of $\Delta$B=0.48 mag, $\Delta$V=0.28 mag, $\Delta$R=0.10 mag, $\Delta$I=0.02 mag, the total energy due to the flare event was measured to be 4.08 ($\pm$0.24) $\times$$10^{34}$ erg, exceeding the superflare level ($10^{34}$). Based on the observations, the evolutionary status of the binary system as well as the long-term orbital period changes were analyzed. It reveals that BX Tri is probably a semi-detached system with the primary component filling its Roche lobe. The extremely high occurrence of flare events of the binary could be related to the rapid mass transfer between components.

In this paper, we present extensive optical and ultraviolet (UV) observations of the type Ia supernova (SN Ia) 2013gs discovered during the Tsinghua-NAOC Transient Survey. The photometric observations in the optical show that the light curves of SN 2013gs is similar to that of normal SNe Ia, with an absolute peak magnitude of $M_{B}$ = $-$19.25 $\pm$ 0.15 mag and a post-maximum decline rate $\Delta$m$_{15}$(B) = 1.00 $ \pm $ 0.05 mag. \emph{Gehrels Swift} UVOT observations indicate that SN 2013gs shows unusually strong UV emission (especially in the $uvw1$ band) at around the maximum light (M$_{uvw1}$ $\sim$ $-$18.9 mag). The SN is characterized by relatively weak Fe~{\sc ii} {\sc iii} absorptions at $\sim$ 5000{\AA} in the early spectra and a larger expansion velocity ($v_{Si}$ $\sim$ 13,000 km s$^{-1}$ around the maximum light) than the normal-velocity SNe Ia. We discuss the relation between the $uvw1-v$ color and some observables, including Si~{\sc ii} velocity, line strength of Si~{\sc ii} $\lambda$6355, Fe~{\sc ii}/{\sc iii} lines and $\Delta m_{15}$(B). Compared to other fast-expanding SNe Ia, SN 2013gs exhibits Si and Fe absorption lines with similar strength and bluer $uvw1-v$ color. We briefly discussed the origin of the observed UV dispersion of SNe Ia.

We present the X-ray timing results of the new black hole candidate (BHC) MAXI J1535-571 during its 2017 outburst from Hard X-ray Modulation Telescope (\emph{Insight}-HXMT) observations taken from 2017 September 6 to 23. Following the definitions given by \citet{Belloni2010}, we find that the source exhibits state transitions from Low/Hard state (LHS) to Hard Intermediate state (HIMS) and eventually to Soft Intermediate state (SIMS). Quasi-periodic oscillations (QPOs) are found in the intermediate states, which suggest different types of QPOs. With the large effective area of \emph{Insight}-HXMT at high energies, we are able to present the energy dependence of the QPO amplitude and centroid frequency up to 100 keV which is rarely explored by previous satellites. We also find that the phase lag at the type-C QPOs centroid frequency is negative (soft lags) and strongly correlated with the centroid frequency. By assuming a geometrical origin of type-C QPOs, the source is consistent with being a high inclination system.

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

This paper presents the second data release (DR2) of the Beijing-Arizona Sky Survey (BASS). BASS is an imaging survey of about 5400 deg$^2$ in $g$ and $r$ bands using the 2.3 m Bok telescope. DR2 includes the observations as of July 2017 obtained by BASS and Mayall $z$-band Legacy Survey (MzLS). This is our first time to include the MzLS data covering the same area as BASS. BASS and MzLS have respectively completed about 72% and 76% of their observations. The two surveys will be served for the spectroscopic targeting of the upcoming Dark Energy Spectroscopic Instrument. Both BASS and MzLS data are reduced by the same pipeline. We have updated the basic data reduction and photometric methods in DR2. In particular, source detections are performed on stacked images, and photometric measurements are co-added from single-epoch images based on these sources. The median 5$\sigma$ depths with corrections of the Galactic extinction are 24.05, 23.61, and 23.10 mag for $g$, $r$, and $z$ bands, respectively. The DR2 data products include stacked images, co-added catalogs, and single-epoch images and catalogs. The BASS website (http://batc.bao.ac.cn/BASS/) provides detailed information and links to download the data.

We present the photometric calibration of the Beijing-Arizona Sky Survey (BASS) and Mayall z-band Legacy Survey (MzLS), which are two of the three wide-field optical legacy imaging surveys to provide the baseline targeting data for the Dark Energy Spectroscopic Instrument (DESI) project. The method of our photometric calibration is subdivided into the external and internal processes. The former utilizes the point-source objects of Pan-STARRS1 survey (PS1) as the reference standards to achieve the zero points of the absolute flux for individual exposures. And then the latter revise the zero points to make them consistent across the survey based on multiple tilings and large offset overlaps. Our process achieves a homogeneous photometric calibration over most of the sky with precision better than 10 mmag for g and r bands, 15 mmag for z band. The accuracy of the calibration is better than 1% at the bright end (16-18 mag) over most of the survey area.

We investigate the hypothesized tidal disruption of the Hercules ultra-faint dwarf galaxy (UFD). Previous tidal disruption studies of the Hercules UFD have been hindered by the high degree of foreground contamination in the direction of the dwarf. We bypass this issue by using RR Lyrae stars, which are standard candles with a very low field-volume density at the distance of Hercules. We use wide-field imaging from the Dark Energy Camera on CTIO to identify candidate RR Lyrae stars, supplemented with observations taken in coordination with the Beijing-Arizona Sky Survey on the Bok Telescope. Combining color, magnitude, and light-curve information, we identify three new RR Lyrae stars associated with Hercules. All three of these new RR Lyrae stars lie outside its published tidal radius. When considered with the nine RR Lyrae stars already known within the tidal radius, these results suggest that a substantial fraction of Hercules' stellar content has been stripped. With this degree of tidal disruption, Hercules is an interesting case between a visibly disrupted dwarf (such as the Sagittarius dwarf spheroidal galaxy) and one in dynamic equilibrium. The degree of disruption also shows that we must be more careful with the ways we determine object membership when estimating dwarf masses in the future. One of the three discovered RR Lyrae stars sits along the minor axis of Hercules, but over two tidal radii away. This type of debris is consistent with recent models that suggest Hercules' orbit is aligned with its minor axis.

Finding the electromagnetic (EM) counterpart of binary compact star merger, especially the binary neutron star (BNS) merger, is critically important for gravitational wave (GW) astronomy, cosmology and fundamental physics. On Aug. 17, 2017, Advanced LIGO and \textit{Fermi}/GBM independently triggered the first BNS merger, GW170817, and its high energy EM counterpart, GRB 170817A, respectively, resulting in a global observation campaign covering gamma-ray, X-ray, UV, optical, IR, radio as well as neutrinos. The High Energy X-ray telescope (HE) onboard \textit{Insight}-HXMT (Hard X-ray Modulation Telescope) is the unique high-energy gamma-ray telescope that monitored the entire GW localization area and especially the optical counterpart (SSS17a/AT2017gfo) with very large collection area ($\sim$1000 cm$^2$) and microsecond time resolution in 0.2-5 MeV. In addition, \textit{Insight}-HXMT quickly implemented a Target of Opportunity (ToO) observation to scan the GW localization area for potential X-ray emission from the GW source. Although it did not detect any significant high energy (0.2-5 MeV) radiation from GW170817, its observation helped to confirm the unexpected weak and soft nature of GRB 170817A. Meanwhile, \textit{Insight}-HXMT/HE provides one of the most stringent constraints (~10$^{-7}$ to 10$^{-6}$ erg/cm$^2$/s) for both GRB170817A and any other possible precursor or extended emissions in 0.2-5 MeV, which help us to better understand the properties of EM radiation from this BNS merger. Therefore the observation of \textit{Insight}-HXMT constitutes an important chapter in the full context of multi-wavelength and multi-messenger observation of this historical GW event.

We present ten $R$-band photometric observations of eight different transits of the hot Jupiter HAT-P-33b, which has been targeted by our Transiting Exoplanet Monitoring Project (TEMP). The data were obtained by two telescopes at the Xinglong Station of National Astronomical Observatories of China (NAOC) from 2013 December through 2016 January, and exhibit photometric scatter of $1.6-3.0\,\rm{mmag}$. After jointly analyzing the previously published photometric data, radial-velocity (RV) measurements, and our new light curves, we revisit the system parameters and orbital ephemeris for the HAT-P-33b system. Our results are consistent with the published values except for the planet-to-star radius ratio ($R_{P}/R_{*}$), the ingress/egress duration ($\tau$) and the total duration ($T_{14}$), which together indicate a slightly shallower and shorter transit shape. Our results are based on more complete light curves, whereas the previously published work had only one complete transit light curve. No significant anomalies in Transit Timing Variations (TTVs) are found, and we place upper mass limits on potential perturbers, largely supplanting the loose constraints provided by the extant RV data. The TTV limits are stronger near mean-motion resonances, especially for the low-order commensurabilities. We can exclude the existence of a perturber with mass larger than 0.6, 0.3, 0.5, 0.5, and $0.3\,{\rm M_\oplus}$ near the 1:3, 1:2, 2:3, 3:2, and 2:1 resonances, respectively.

GC-1 and GC-2 are two globular clusters (GCs) in the remote halo of M81 and M82 in the M81 group discovered by Jang et al. using the {\it Hubble Space Telescope} ({\it HST}) images. These two GCs were observed as part of the Beijing--Arizona--Taiwan--Connecticut (BATC) Multicolor Sky Survey, using 14 intermediate-band filters covering a wavelength range of 4000--10000 \AA. We accurately determine these two clusters' ages and masses by comparing their spectral energy distributions (from 2267 to 20000~{\AA}, comprising photometric data in the near-ultraviolet of the {\it Galaxy Evolution Explorer}, 14 BATC intermediate-band, and Two Micron All Sky Survey near-infrared $JHK_{\rm s}$ filters) with theoretical stellar population-synthesis models, resulting in ages of $15.50\pm3.20$ for GC-1 and $15.10\pm2.70$ Gyr for GC-2. The masses of GC-1 and GC-2 obtained here are $1.77-2.04\times 10^6$ and $5.20-7.11\times 10^6 \rm~M_\odot$, respectively. In addition, the deep observations with the Advanced Camera for Surveys and Wide Field Camera 3 on the {\it HST} are used to provide the surface brightness profiles of GC-1 and GC-2. The structural and dynamical parameters are derived from fitting the profiles to three different models; in particular, the internal velocity dispersions of GC-1 and GC-2 are derived, which can be compared with ones obtained based on spectral observations in the future. For the first time, in this paper, the $r_h$ versus $M_V$ diagram shows that GC-2 is an ultra-compact dwarf in the M81 group.

The Beijing-Arizona Sky Survey (BASS) is a wide-field two-band photometric survey of the Northern Galactic Cap using the 90Prime imager on the 2.3 m Bok telescope at Kitt Peak. It is a four-year collaboration between the National Astronomical Observatory of China and Steward Observatory, the University of Arizona, serving as one of the three imaging surveys to provide photometric input catalogs for target selection of the Dark Energy Spectroscopic Instrument (DESI) project. BASS will take up to 240 dark/grey nights to cover an area of about 5400 deg$^2$ in the $g$ and $r$ bands. The 5$\sigma$ limiting AB magnitudes for point sources in the two bands, corrected for the Galactic extinction, are 24.0 and 23.4 mag, respectively. BASS, together with other DESI imaging surveys, will provide unique science opportunities that cover a wide range of topics in both Galactic and extragalactic astronomy.

In this paper, we collect a sample of 81 ellipsoidal red giant binaries in the Large Magellanic Cloud (LMC), and we study their orbital natures individually and statistically. The sample contains 59 systems with circular orbits and 22 systems with eccentric orbits. We derive orbital solutions using the 2010 version of the Wilson-Devinney code. The sample is selection-bias corrected, and the orbital parameter distributions are compared to model predictions for the LMC and to observations in the solar vicinity. The masses of the red giant primaries are found to range from about 0.6 to 9 Msun with a peak at around 1.5 Msun, in agreement with studies of the star formation history of the LMC, which find a burst of star formation beginning around 4 Gyr ago. The observed distribution of mass ratios q=m2/m1 is more consistent with the flat q distribution derived for the solar vicinity by Raghavan et al. than it is with the solar vicinity q distribution derived by Duquennoy & Mayor. There is no evidence for an excess number of systems with equal mass components. We find that about 20% of the ellipsoidal binaries have eccentric orbits, twice the fraction estimated by Soszynski et al. Our eccentricity evolution test shows that the existence of eccentric ellipsoidal red giant binaries on the upper parts of the red giant branch (RGB) can only be explained if tidal circularization rates are ~1/100 the rates given by the usual theory of tidal dissipation in convective stars.

We study the integral Galactic extinction and reddening based on the galaxy catalog of the South Galactic Cap U-band Sky Survey (SCUSS), where $u$ band galaxy number counts and $u-r$ color distribution are used to derive the Galactic extinction and reddening respectively. We compare these independent statistical measurements with the reddening map of \citet{Schlegel1998}(SFD) and find that both the extinction and reddening from the number counts and color distribution are in good agreement with the SFD results at low extinction regions ($E(B-V)^{SFD}<0.12$ mag). However, for high extinction regions ($E(B-V)^{SFD}>0.12$ mag), the SFD map overestimates the Galactic reddening systematically, which can be approximated by a linear relation $\Delta E(B-V)= 0.43[E(B-V)^{SFD}-0.12$]. By combing the results of galaxy number counts and color distribution together, we find that the shape of the Galactic extinction curve is in good agreement with the standard $R_V=3.1$ extinction law of \cite{ODonnell1994}.

We present and analyze the possibility of using optical ${\it u}$-band luminosities to estimate star-formation rates (SFRs) of galaxies based on the data from the South Galactic Cap ${\it u }$ band Sky Survey (SCUSS), which provides a deep ${\it u}$-band photometric survey covering about 5000 $deg^2$ of the South Galactic Cap. Based on two samples of normal star-forming galaxies selected by the BPT diagram, we explore the correlations between ${\it u}$-band, H$\alpha$, and IR luminosities by combing SCUSS data with the Sloan Digital Sky Survey (SDSS) and ${\it Wide}$-${\it field\ Infrared\ Survey\ Explorer}$ (${\it WISE}$). The attenuation-corrected ${\it u}$-band luminosities are tightly correlated with the Balmer decrement-corrected H$\alpha$ luminosities with an rms scatter of $\sim$ 0.17 dex. The IR-corrected ${\it u }$ luminosities are derived based on the correlations between the attenuation of ${\it u}$-band luminosities and ${\it WISE}$ 12 (or 22) $\mu$m luminosities, and then calibrated with the Balmer-corrected H$\alpha$ luminosities. The systematic residuals of these calibrations are tested against the physical properties over the ranges covered by our sample objects. We find that the best-fitting nonlinear relations are better than the linear ones and recommended to be applied in the measurement of SFRs. The systematic deviations mainly come from the pollution of old stellar population and the effect of dust extinction; therefore, a more detailed analysis is needed in the future work.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

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

Context: We present a newly discovered class of low-luminosity, dusty, evolved objects in the Magellanic Clouds. These objects have dust excesses, stellar parameters, and spectral energy distributions similar to those of dusty post-asymptotic giant branch (post-AGB) stars. However, they have lower luminosities and hence lower masses. We suggest that they have evolved off the red giant branch (RGB) instead of the AGB as a result of binary interaction. Aims: In this study we aim to place these objects in an evolutionary context and establish an evolutionary connection between RGB binaries (such as the sequence E variables) and our new sample of objects. Methods: We compared the theoretically predicted birthrates of the progeny of RGB binaries to the observational birthrates of the new sample of objects. Results: We find that there is order-of-magnitude agreement between the observed and predicted birthrates of post-RGB stars. The sources of uncertainty in the birthrates are discussed; the most important sources are probably the observational incompleteness factor and the post-RGB evolution rates. We also note that mergers are relatively common low on the RGB and that stars low on the RGB with mid-IR excesses may recently have undergone a merger. Conclusions: Our sample of dusty post-RGB stars most likely provides the first observational evidence for a newly discovered phase in binary evolution: post-RGB binaries with circumstellar dust.

The SCUSS is a deep $u$-band imaging survey in the south Galactic cap using the 2.3m Bok telescope. The survey observations were completed in the end of 2013, covering an area of about 5000 square degrees. We release the data in the region with an area of about 4000 deg$^2$ that is mostly covered by the Sloan digital sky survey. The data products contain calibrated single-epoch images, stacked images, photometric catalogs, and a catalog of star proper motions derived by Peng et al, 2015. The median seeing and magnitude limit ($5\sigma$) are about 2".0 and 23.2 mag, respectively. There are about 8 million objects having measurements of absolute proper motions. All the data and related documentations can be accessed through the SCUSS data release website of \url{http://batc.bao.ac.cn/Uband/data.html}.

We present a new selection technique of producing spectroscopic target catalogues for massive spectroscopic surveys for cosmology. This work was conducted in the context of the extended Baryon Oscillation Spectroscopic Survey (eBOSS), which will use ~200 000 emission line galaxies (ELGs) at 0.6<zspec<1.0 to obtain a precise baryon acoustic oscillation measurement. Our proposed selection technique is based on optical and near-infrared broad-band filter photometry. We used a training sample to define a quantity, the Fisher discriminant (linear combination of colours), which correlates best with the desired properties of the target: redshift and [OII] flux. The proposed selections are simply done by applying a cut on magnitudes and this Fisher discriminant. We used public data and dedicated SDSS spectroscopy to quantify the redshift distribution and [OII] flux of our ELG target selections. We demonstrate that two of our selections fulfil the initial eBOSS/ELG redshift requirements: for a target density of 180 deg^2, ~70% of the selected objects have 0.6<zspec<1.0 and only ~1% of those galaxies in the range 0.6<zspec<1.0 are expected to have a catastrophic zspec estimate. Additionally, the stacked spectra and stacked deep images for those two selections show characteristic features of star-forming galaxies. The proposed approach using the Fisher discriminant could, however, be used to efficiently select other galaxy populations, based on multi-band photometry, providing that spectroscopic information is available. This technique could thus be useful for other future massive spectroscopic surveys such as PFS, DESI, and 4MOST.

In this paper, we report the detections of stellar variabilities from the first 2-year observations of sky area of about 1300 square degrees from the Tsinghua University-NAOC Transient Survey (TNTS). A total of 1237 variable stars (including 299 new ones) were detected with brightness < 18.0 mag and magnitude variation >= 0.1 mag on a timescale from a few hours to few hundred days. Among such detections, we tentatively identified 661 RR Lyrae stars, 431 binaries, 72 Semiregular pulsators, 29 Mira stars, 11 slow irregular variables, 11 RS Canum Venaticorum stars, 7 Gamma Doradus stars, 5 long period variables, 3 W Virginis stars, 3 Delta Scuti stars, 2 Anomalous Cepheids, 1 Cepheid, and 1 nove-like star based on their time-series variability index Js and their phased diagrams. Moreover, we found that 14 RR Lyrae stars show the Blazhko effect and 67 contact eclipsing binaries exhibit the O'Connell effect. Since the period and amplitude of light variations of RR Lyrae variables depend on their chemical compositions, their photometric observations can be used to investigate distribution of metallicity along the direction perpendicular to the Galactic disk. We find that the metallicity of RR Lyrae stars shows large scatter at regions closer to the Galactic plane (e.g., -3.0 < [Fe/H] < 0) but tends to converge at [Fe/H]~ -1.7 at larger Galactic latitudes. This variation may be related to that the RRAB Lyrae stars in the Galactic halo come from globular clusters with different metallicity and vertical distances, i.e. OoI and OoII populations, favoring for the dual-halo model.

The South Galactic Cap u-band Sky Survey (SCUSS) is a deep u-band imaging survey in the Southern Galactic Cap, using the 90Prime wide-field imager on the 2.3m Bok telescope at Kitt Peak. The survey observations started in 2010 and ended in 2013. The final survey area is about 5000 deg2 with a median 5-sigma point source limiting magnitude of about 23.2. This paper describes the survey data reduction process, which includes basic imaging processing, astrometric and photometric calibrations, image stacking, and photometric measurements. Survey photometry is performed on objects detected both on SCUSS u-band images and in the SDSS database. Automatic, aperture, point-spread function (PSF), and model magnitudes are measured on stacked images. Co-added aperture, PSF, and model magnitudes are derived from measurements on single-epoch images. We also present comparisons of the SCUSS photometric catalog with those of the SDSS and CFHTLS.

SCUSS is a u-band photometric survey covering about 4000 square degree of the South Galactic Cap, reaching depths of up to 23 mag. By extending around 1.5 mag deeper than SDSS single-epoch u data, SCUSS is able to probe much a larger volume of the outer halo, i.e. with SCUSS data blue horizontal branch (BHB) stars can trace the outer halo of the Milky Way as far as 100-150 kpc. Utilizing this advantage we combine SCUSS u band with SDSS DR9 gri photometric bands to identify BHB stars and explore halo substructures. We confirm the existence of the Pisces overdensity, which is a structure in the outer halo (at around 80 kpc) that was discovered using RR Lyrae stars. For the first time we are able to determine its spatial extent, finding that it appears to be part of a stream with a clear distance gradient. The stream, which is ~5 degrees wide and stretches along ~25 degrees, consists of 20-30 BHBs with a total significance of around 6sigma over the background. Assuming we have detected the entire stream and that the progenitor has fully disrupted, then the number of BHBs suggests the original system was similar to smaller classical or a larger ultra-faint dwarf galaxy. On the other hand, if the progenitor still exists, it can be hunted for by reconstructing its orbit from the distance gradient of the stream. This new picture of the Pisces overdensity sheds new light on the origin of this intriguing system.

Based on SDSS g, r and SCUSS (South Galactic Cap of u-band Sky Survey) $u$ photometry, we develop a photometric calibration for estimating the stellar metallicity from $u-g$ and $g-r$ colors by using the SDSS spectra of 32,542 F- and G-type main sequence stars, which cover almost $3700$ deg$^{2}$ in the south Galactic cap. The rms scatter of the photometric metallicity residuals relative to spectrum-based metallicity is $0.14$ dex when $g-r<0.4$, and $0.16$ dex when $g-r>0.4$. Due to the deeper and more accurate magnitude of SCUSS $u$ band, the estimate can be used up to the faint magnitude of $g=21$. This application range of photometric metallicity calibration is wide enough so that it can be used to study metallicity distribution of distant stars. In this study, we select the Sagittarius (Sgr) stream and its neighboring field halo stars in south Galactic cap to study their metallicity distribution. We find that the Sgr stream at the cylindrical Galactocentric coordinate of $R\sim 19$ kpc, $\left| z\right| \sim 14$ kpc exhibits a relative rich metallicity distribution, and the neighboring field halo stars in our studied fields can be modeled by two-Gaussian model, with peaks respectively at [Fe/H]$=-1.9$ and [Fe/H]$=-1.5$.

Ellipsoidal red giant binaries are close binary systems where an unseen, relatively close companion distorts the red giant, leading to light variations as the red giant moves around its orbit. These binaries are likely to be the immediate evolutionary precursors of close binary planetary nebula and post-asymptotic giant branch and post-red giant branch stars. Due to the MACHO and OGLE photometric monitoring projects, the light variability nature of these ellipsoidal variables has been well studied. However, due to the lack of radial velocity curves, the nature of their masses, separations, and other orbital details has so far remained largely unknown. In order to improve this situation, we have carried out spectral monitoring observations of a large sample of 80 ellipsoidal variables in the Large Magellanic Cloud and we have derived radial velocity curves. At least 12 radial velocity points with good quality were obtained for most of the ellipsoidal variables. The radial velocity data are provided with this paper. Combining the photometric and radial velocity data, we present some statistical results related to the binary properties of these ellipsoidal variables.

Absolute proper motions for $\sim$ 7.7 million objects were derived based on data from the South Galactic Cap u-band Sky Survey (SCUSS) and astrometric data derived from uncompressed Digitized Sky Surveys that the Space Telescope Science Institute (STScI) created from the Palomar and UK Schmidt survey plates. We put a great deal of effort into correcting the position-, magnitude-, and color-dependent systematic errors in the derived absolute proper motions. The spectroscopically confirmed quasars were used to test the internal systematic and random error of the proper motions. The systematic errors of the overall proper motions in the SCUSS catalog are estimated as -0.08 and -0.06 mas/yr for {\mu}{\alpha} cos {\delta} and {\mu}{\delta}, respectively. The random errors of the proper motions in the SCUSS catalog are estimated independently as 4.2 and 4.4 mas/yr for {\mu}{\alpha} cos {\delta} and {\mu}{\delta}. There are no obvious position-, magnitude-, and color-dependent systematic errors of the SCUSS proper motions. The random error of the proper motions goes up with the magnitude from about 3 mas/yr at u < 18.0 mag to about 7 mas/yr at u = 22.0 mag. The proper motions of stars in SCUSS catalog are compared with those in the SDSS catalog, and they are highly consistent.

The South Galactic Cap $u$-band Sky Survey (SCUSS) provides a deep $u$-band imaging of about 5000 deg$^2$ in south Galactic cap. It is about 1.5 mag deeper than the SDSS $u$-band. In this paper we evaluate the capability of quasar selection using both SCUSS and SDSS data, based on considerations of the deep SCUSS $u$-band imaging and two-epoch $u$-band variability. We find that the combination of the SCUSS $u$-band and the SDSS $griz$ band allows us to select more faint quasars and more quasars at redshift around 2.2 than the selection only with the SDSS $ugriz$ data. Quasars have significant $u$-band variabilities. The fraction of quasars with large two-epoch variability is much higher than that of stars. The selection by variability can select both low-redshift quasars with ultraviolet excess and mid-redshift ($2 < z <3.5$) quasars where quasar selection by optical colors is inefficient. The above two selections are complementary and make full use of the SCUSS u-band advantages.

We present the [OII] luminosity function measured in the redshift range 0.1<z<1.65 with unprecedented depth and accuracy. Our measurements are based on medium resolution flux-calibrated spectra of emission line galaxies with the FORS2 instrument at VLT and with the SDSS-III/BOSS spectrograph. The FORS2 spectra and the corresponding catalog containing redshifts and line fluxes are released along with this paper. In this work we use a novel method to combine the aforementioned surveys with GAMA, zCOSMOS and VVDS, which have different target selection, producing a consistent weighting scheme to derive the [OII] luminosity function. The measured luminosity function is in good agreement with previous independent estimates. The comparison with two state-of-the-art semi-analytical models is good, which is encouraging for the production of mock catalogs of [OII] flux limited surveys. We observe the bright end evolution over 8.5 Gyr: we measure the decrease of log L* from 42.4 erg/s at redshift 1.44 to 41.2 at redshift 0.165 and we find that the faint end slope flattens when redshift decreases. This measurement confirms the feasibility of the target selection of future baryonic acoustic oscillation surveys aiming at observing [OII] flux limited samples.

The Tsinghua University-National Astronomical Observatories of China (NAOC) Transient Survey (TNTS) is an automatic survey for a systematic exploration of optical transients (OTs), conducted with a 60/90 cm Schmidt telescope at Xinglong station of NAOC. This survey repeatedly covers ~ 1000 square degrees of the north sky with a cadence of 3-4 days. With an exposure of 60 s, the survey reaches a limited unfiltered magnitude of about 19.5 mag. This enables us to discover supernovae at their relatively young stages. In this paper, we describe the overall performance of our survey during the first year and present some preliminary results.

This paper presents CCD multicolor photometry for 304 old star clusters in the nearby spiral galaxy M31. Of which photometry of 55 star clusters is first obtained. The observations were carried out as a part of the Beijing--Arizona--Taiwan--Connecticut (BATC) Multicolor Sky Survey from 1995 February to 2008 March, using 15 intermediate-band filters covering 3000--10000 \AA. Detailed comparisons show that our photometry is in agreement with previous measurements. Based on the ages and metallicities from Caldwell et al. and the photometric measurements here, we estimated the clusters' masses by comparing their multicolor photometry with stellar population synthesis models. The results show that the sample clusters have masses between $\sim 3\times10^4 M_\odot$ and $\sim 10^7 M_\odot$ with the peak of $\sim 4\times10^5 M_\odot$. The masses here are in good agreement with those in previous studies. Combined with the masses of young star clusters of M31 from Wang et al., we find that the peak of mass of old clusters is ten times that of young clusters.

The Graphics Processing Unit (GPU) has become an integral part of astronomical instrumentation, enabling high-performance online data reduction and accelerated online signal processing. In this paper, we describe a wide-band reconfigurable spectrometer built using an off-the-shelf GPU card. This spectrometer, when configured as a polyphase filter bank (PFB), supports a dual-polarization bandwidth of up to 1.1 GHz (or a single-polarization bandwidth of up to 2.2 GHz) on the latest generation of GPUs. On the other hand, when configured as a direct FFT, the spectrometer supports a dual-polarization bandwidth of up to 1.4 GHz (or a single-polarization bandwidth of up to 2.8 GHz).

Based on SDSS and South Galactic Cap of u-band Sky Survey (SCUSS) early data, we use star counts method to estimate the Galactic structure parameters in an intermediate latitude with 10,180 main-sequence (MS) stars in absolute magnitude interval of $4 \leq M_r \leq 13$. We divide the absolute magnitude into five intervals:$4 \leq M_r < 5$, $5 \leq M_r < 6$, $6 \leq M_r < 8$, $8 \leq M_r < 10$, $10 \leq M_r \leq 13$, and estimate the Galactic structure parameters in each absolute magnitude interval to explore their possible variation with the absolute magnitude. Our study shows the parameters depend on absolute magnitude. For the thin disk, the intrinsic faint MS stars have large local space density and they tend to stay close to the Galactic plane. A plausible explanation is that faint MS stars with long lifetime experience long gravitational interaction time result in a short scaleheight. However, for the thick disk, the parameters show a complex trend with absolute magnitude, which may imply the complicated original of the thick disk. For the halo, the intrinsic faint MS stars have large local density and small axial ratio, which indicate a flattened inner halo and a more spherical outer halo.

Sequence E variables are close binary red giants that show ellipsoidal light variations. They are likely the immediate precursors of planetary nebulae (PNe) with close binary central. We have made a Monte Carlo simulation to determine the fraction of red giant binaries that go through a common envelope (CE) event leading to the production of a close binary system or a merged star. The novel aspect of this simulation is that we use the observed frequency of sequence E binaries in the LMC to normalize our calculations. In our standard model, we find that in the LMC today the fraction of PNe with close binary central stars is 7-9%, the fraction of PNe with intermediate period binary central stars having separations capable of influencing the nebula shape (P<500 yrs) is 23-27%, the fraction of PNe containing wide binaries that are unable to influence the nebula shape (P>500 yrs) is 46-55%, the fraction of PNe derived from single stars is 3-19% and 5-6% of PNe are produced by previously merged stars. We also predict that the birthrate of post-RGB stars is ~4% of the total PN birthrate, equivalent to ~50% of the production rate of PNe with close binary central stars. These post-RGB stars most likely appear initially as luminous low-mass helium white dwarf binaries. We use our model and the observed number of red giant stars in the top one magnitude of the RGB in the LMC to predict the number of PNe in the LMC. We predict 548 PNe in good agreement with the 541+/-89 PNe observed by Reid & Parker (2006). Since most of these PNe come from single or non-interacting binary stars in our model, this means that most such stars produce PNe contrary to the "Binary Hypothesis" which suggests that binary interaction is required to produce a PN.

In this paper, the rms-flux (root mean square-flux) relation along the Z-track of the bright Z-Source Cyg X-2 is analyzed using the observational data of Rossi X-ray Timing Explorer (RXTE). Three types of rms-flux relations, i.e. positive, negative, and 'arch'-like correlations are found in different branches. The rms is positively correlated with flux in normal branch (NB), but anti-correlated in the vertical horizontal branch (VHB). The rms-flux relation shows an 'arch'-like shape in the horizontal branch (HB). We also try to explain this phenomenon using existing models.

Sequence E variables are close binary red giants that show ellipsoidal light variations. They are likely to terminate their red giant evolution by a common envelope (CE) event when the red giant fills its Roche lobe, and produce close binary Planetary Nebulae (PNe). We made a Monte Carlo simulation to predict the fraction of Planetary Nebulae Nuclei (PNNe) that are post-CE binaries, using the observed frequency of sequence E binaries in the LMC to normalize our calculations. We find that 10-16% of PNNe should be short period, post-CE binaries.

The multi-epoch YJKs sub-arcsecond photometry of the VMC survey provides a long anticipated deep near-infrared (NIR) window into further understanding the stellar populations of the Magellanic Clouds. The first year of observations consisted of six tiles covering ~9% of the Large Magellanic Cloud (LMC) survey region and contains 102 objects previously classified as planetary nebulae (PNe). A large proportion of the sample were found to be contaminated by non-PNe. These initial results underline the importance of establishing a clean catalogue of LMC PNe before they are applied in areas such as the planetary nebula luminosity function (PNLF) and searches for binary central stars. As the VMC survey progresses it will play a fundamental role in cleaning extant PN catalogues and a complementary role in the discovery of new PNe.

The Large Magellanic Cloud (LMC) contains the nearest large extragalactic population of planetary nebulae (PNe). A shallow viewing angle and low interstellar reddening towards the LMC potentially means a larger, more complete flux-limited population can be assembled than for any other galaxy. These advantages appear to be reflected by the small gap between the catalogued ($\sim$700 PNe) and estimated ($1000\pm250$ PNe) population size. With more detailed multi-wavelength studies the catalogued number of LMC PNe may fall, potentially widening this gap. We demonstrate here that the gap can be further bridged with improved optical and near-infrared imaging surveys. We present three [O III]-selected PNe discovered from ESO WFI observations of the 30 Doradus region and one serendipitous discovery from near-infrared Vista Magellanic Cloud (VMC) survey observations. The WFI PNe have resolved [O III] and H$\alpha$ nebulae that verify their PN nature and their [O III] fluxes place them 6--7 mag ($m_{5007}=20$--21 mag) fainter than the bright-end of the planetary nebula luminosity function (PNLF). Their faintness, small angular size and surrounding complex emission-line background explains why previous H$\alpha$ surveys of the region did not select them. We estimate there may be as many as 50--75 similar PNe awaiting discovery in the central $5\times5$ degrees of the LMC. The VMC survey routinely detects PNe as red resolved nebulae that may allow some of this expected population to be recovered without traditional narrow-band imaging surveys. We demonstrate this potential with the first new VMC-selected PN which has a rare Wolf-Rayet [WC9]--[WC11] central star.

The origin of the long secondary periods (LSPs) in red variables remains a mystery up to now, although there exist many models. The light curves of some LSPs stars mimic an eclipsing binary with a pulsating red giant component. To test this hypothesis, the observational data of two LSP variable red giants, 77.7795.29 and 77.8031.42, discovered by the MACHO project from the LMC, are collected and analyzed. The probable eclipsing features of the light curves are simulated by the Wilson-Devinney (W-D) method. The simulation yields a contact and a semidetached geometry for the two systems, respectively. In addition, the pulsation constant of the main pulsating component in each binary system is derived. By combining the results of the binary model and the pulsation component, we investigate the feasibility of the pulsating binary model. It is found that the radial velocity curve expected from the binary model has a much larger amplitude than the observed one and a period double the observed one. Furthermore, the masses of the components based on the density derived from the binary orbit solution are too low to be compatible with both the evolutionary stage and the high luminosity. Although the pulsation mode identified by the pulsation constant which is dependent on the density from the binary-model is consistent with the first or second overtone radial pulsation, we conclude that the pulsating binary model is a defective model for the LSP.

A remarkable similarity between the large-scale non-Gaussian pattern of cosmic microwave background (CMB) temperatures obtained by Wilkinson Microwave Anisotropy Probe (WMAP) mission and the distribution feature of observation numbers is noted. Motivated from such a similarity, in this work we check the WMAP data for the correlation between pixel temperature t and observation number N. Systematic effect of imbalance differential observation and significant t-N correlation in magnitude, distribution non-Gaussianity and north-south asymmetry are found. Our results indicate that, for precision cosmology study based on WMAP observations, the observation effect on released WMAP temperature maps has to be further carefully studied.

We apply a Cross-correlation (CC) method developed previously for detecting gamma-ray point sources to the WMAP first year data by using the Point-Spread Function of WMAP and obtain a full sky CC coefficient map. Analyzing this map, we find that the CC method is a powerful tool to examine the WMAP foreground residuals which can be further cleaned accordingly. Evident foreground signals are found in WMAP foreground cleaned maps and Tegmark cleaned map. In this process 101 point-sources are detected, and 26 of them are new sources besides the originally listed WMAP 208 sources. We estimate the flux of these new sources and verify them by another method. As a result, a revised mask file based on the WMAP first year data is produced by including these new sources.