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.

It is generally believed that the variability of photospheric emission in gamma-ray bursts (GRBs) traces that of the jet power. This work further investigates the variability of photospheric emission in a variable jet. By setting a constant $\eta$ (dimensionless entropy of the jet), we find that the light curve of the photospheric emission shows a ``tracking'' pattern on the time profile of jet power. However, the relative variability is significantly low in the photospheric emission compared with that in the jet power. If the $\eta$ is genetic variable, the variability of the photospheric emission is not only limited by the jet power but also affected by $\eta$ strongly. It becomes complex and is generally different from that of the jet power. Moreover, the opposite phase may stand in the variabilities of the photospheric emission at different photon energies. We also find that the relative variability does not remain constant over the photon energies with an obvious reduction at a certain energy. This is consistent with the analysis of GRB 090902B in which an appreciable thermal component has been detected in a wide energy range. For several other GRBs coupling with the thermal component, we conservatively evaluate the variability of the thermal and non-thermal emission, respectively. Our results show that the relative variability of the thermal emission is likely comparable to that of the non-thermal emission for these bursts. In addition, the analysis of GRB~120323A reveals that the variability of the photospheric emission may be of the opposite phase from that of the non-thermal emission.

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.

Deep observations have detected extended Ly$\alpha$ emission nebulae surrounding tens of quasars at redshift 2 to 6. However, the metallicity of such extended haloes is still poorly understood. We perform a detailed analysis on a large sample of 80 quasars at $z\sim3$ based on MUSE-VLT data. We find clear evidence of extended emission of the UV nebular lines such as CIV or HeII for about 20$\%$ of the sample, while CIII] is only marginally detected in a few objects. By stacking the cubes we detect emission of CIV, HeII and CIII] out to a radius of about 45 kpc. CIV and HeII show a radial decline much steeper than Ly$\alpha$, while CIII] shows a shallower profile similar to Ly$\alpha$ in the inner 45 kpc. We infer that the average metallicity of the circumgalactic gas within the central 30-50~kpc is $\sim$0.5 solar, or even higher. However, we also find evidence of a component of the Ly$\alpha$ haloes, which has much weaker metal emission lines relative to Ly$\alpha$. We suggest that the high metallicity of the circumgalactic medium within the central 30-50 kpc is associated with chemical pre-enrichment by past quasar-driven outflows and that there is a more extended component of the CGM that has much lower metallicity and likely associated with near-pristine gas accreted from the intergalactic medium. We show that our observational results are in good agreement with the expectations of the FABLE zoom-in cosmological simulations.

Quasar broad emission lines are largely powered by photoionization from the accretion continuum. Increased central luminosity will enhance line emissivity in more distant clouds, leading to increased average distance of the broad-line-emitting clouds and decreased averaged line width, known as the broad-line region (BLR) "breathing". However, different lines breathe differently, and some high-ionization lines, such as C IV, can even show "anti-breathing" where the line broadens when luminosity increases. Using multi-year photometric and spectroscopic monitoring data from the Sloan Digital Sky Survey Reverberation Mapping project, we quantify the breathing effect ($\Delta$log W=$\alpha\Delta$log L) of broad H$\alpha$, H$\beta$, Mg II, C IV, and C III] for statistical quasar samples over $z\approx 0.1-2.5$. We found that H$\beta$ displays the most consistent normal breathing expected from the virial relation ($\alpha\sim-0.25$), Mg II and H$\alpha$ on average show no breathing ($\alpha\sim 0$), and C IV (and similarly C III] and Si IV mostly shows anti-breathing ($\alpha>0$). The anti-breathing of C IV can be well understood by the presence of a non-varying core component in addition to a reverberating broad-base component, consistent with earlier findings. The deviation from canonical breathing introduces extra scatter (a luminosity-dependent bias) in single-epoch virial BH mass estimates due to intrinsic quasar variability, which underlies the long argued caveats of C IV single-epoch masses. Using the line dispersion instead of FWHM leads to less, albeit still substantial, deviations from canonical breathing in most cases. Our results strengthen the need for reverberation mapping to provide reliable quasar BH masses, and quantify the level of variability-induced bias in single-epoch BH masses based on various lines.

The first charge-coupled device photometric light curves in B, V, R${_c}$, and I${_c}$ bands of the short-period contact binary star OQ Cam are presented. Through analyzing the light curves with the Wilson-Devinney code, it is found that OQ Cam is a W-type shallow contact binary star with a contact degree of f =20.2%. Based on the O-C analyzing, the orbit period displays a long term increasing with a rate of $dP/dt=4.40\times10{^{-7}} day \cdot year{^{-1}}$. The increasing in orbit period can be interpreted by mass transfer from the less massive component to the more massive one. As the orbit period increasing, OQ Cam would evolve from the present shallow contact configuration to a none contact stage. So it may be a potential candidate to confirm the thermal relaxation oscillation theory.

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.

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.

We present our analysis of the LyC emission and escape fraction of 111 spectroscopically verified galaxies with and without AGN from $2.26<z<4.3$. We extended our ERS sample from Smith et al. (2018; arXiv:1602.01555) with 64 galaxies in the GOODS North and South fields using WFC3/UVIS F225W, F275W, and F336W mosaics we independently drizzled using the HDUV, CANDELS, and UVUDF data. Among the 17 AGN from the 111 galaxies, one provided a LyC detection in F275W at $m_{AB}=23.19$ mag (S/N $\simeq$ 133) and $GALEX$ NUV at $m_{AB}=23.77$ mag (S/N $\simeq$ 13). We simultaneously fit $SDSS$ and $Chandra$ spectra of this AGN to an accretion disk and Comptonization model and find $f_{esc}$ values of $f_{esc}^{F275W}\simeq 28^{+20}_{-4}$% and $f_{esc}^{NUV}\simeq 30^{+22}_{-5}$%. For the remaining 110 galaxies, we stack image cutouts that capture their LyC emission using the F225W, F275W, and F336W data of the GOODS and ERS samples, and both combined, as well as subsamples of galaxies with and without AGN, and $all$ galaxies. We find the stack of 17 AGN dominate the LyC production from $\langle z\rangle\simeq 2.3-4.3$ by a factor of $\sim$10 compared to all 94 galaxies without AGN. While the IGM of the early universe may have been reionized mostly by massive stars, there is evidence that a significant portion of the ionizing energy came from AGN.

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.

We present a cosmic ray muon tomographic experiment for measuring the muon flux and imaging the tunnel overburden structures in Changshu, China. The device used in this study is a tracking detector based on the plastic scintillator with SiPM technology, which can be conveniently operated in field works. The compact system with sensitive area of $6400 cm^2$ can measure the angular distribution of cosmic muons. It's able to image the overburden density length from the surface of overburden to the detector along the muon tracks. The open sky muon flux measurement outside the tunnel has a good agreement with the modified Gassier Formula model. The distributions of muon flux at three positions inside the tunnel are very similar to that of open sky. Assuming the average density of overburden compact sandstone is $2.65 g/cm^3$, the overburden thickness can be obtained from the density length derived from the difference of muon flux inside and outside the tunnel. Moreover, for known penetrated lengths (i.e., topography of overburden), the density anomalies of the overburden can also been obtained. This study suggests a potential application for imaging and detecting subsurface structures in civil engineering, tunnels or caverns with the cosmic ray muon telescope.

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

We report the detection of diffuse Lyman-alpha halos (LAHs) around star-forming galaxies at $z \sim 5.7$ by stacking 310 spectroscopically confirmed Lyman-alpha (LAEs). The majority of the LAEs are identified from our spectroscopic survey of galaxies at $z>5.5$. They are all located in well-studied fields with deep narrowband and broadband imaging data. We combine the LAE sample and its subsamples in the narrowband NB816 (i.e., the Lyman-alpha band) and $z$ band (i.e., the continuum band). By comparing the stacked objects with their corresponding point spread functions, we clearly detect extended LAHs around these LAEs. We perform sophisticated simulations and analyses on statistical and systematic errors, and confirm that the detected halos are not caused by errors. The scale lengths of the LAHs, when described by a double-component model, range from 1.2 to 5.3 kpc with a typical value of $\sim$2 kpc. The LAH sizes from our sample are in agreement with those of individual LAEs at the similar redshift measured by VLT/MUSE, but are relatively smaller than those of photometrically selected LAEs in previous studies. We also find that LAEs with higher Lyman-alpha luminosities, higher UV continuum luminosities or smaller Lyman-alpha equivalent widths tend to have larger LAH sizes. Our results are consistent with a scenario that LAHs originate from the scattered light of the central galaxies by H I gas in the circumgalactic medium.

We study six luminous Lyman-alpha emitters (LAEs) with very blue rest-frame UV continua at $5.7\le z \le 6.6$. These LAEs have previous HST and Spitzer IRAC observations. Combining our newly acquired HST images, we find that their UV-continuum slopes $\beta$ are in a range of $-3.4\le \beta \le -2.6$. Unlike previous, tentative detections of $\beta \simeq -3$ in photometrically selected, low-luminosity galaxies, our LAEs are spectroscopically confirmed and luminous ($M_{\rm UV}<-20$ mag). We model their broadband spectral energy distributions (SEDs), and find that two $\beta\simeq-2.6\pm0.2$ galaxies can be well fitted with young and dust-free stellar populations. However, it becomes increasingly difficult to fit bluer galaxies. We explore further interpretations by including non-zero LyC escape fraction $f_{\rm esc}$, very low metallicities, and/or AGN contributions. Assuming $f_{\rm esc}\simeq0.2$, we achieve the bluest slopes $\beta\simeq-2.7$ when nebular emission is considered. This can nearly explain the SEDs of two galaxies with $\beta\simeq-2.8$ and --2.9 ($\sigma_{\beta}=0.15$). Larger $f_{\rm esc}$ values and very low metallicities are not favored by the strong nebular line emission (evidenced by the IRAC flux) or the observed (IRAC 1 - IRAC 2) color. Finally, we find that the $\beta\simeq-2.9$ galaxy can potentially be well explained by the combination of a very young population with a high $f_{\rm esc}$ ($\ge0.5$) and an old, dusty population. We are not able to produce two $\beta \simeq -3.4 \pm0.4$ galaxies. Future deep spectroscopic observations are needed to fully understand these galaxies.

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.

In the standard cosmological model, dark matter drives the structure formation and constructs potential wells within which galaxies may form. The baryon fraction in dark halos can reach the universal value (15.7%) in massive clusters and decreases rapidly as the mass of the system decreases. The formation of dwarf galaxies is sensitive both to baryonic processes and the properties of dark matter owing to the shallow potential wells in which they form. In dwarf galaxies in the Local Group, dark matter dominates the mass content even within their optical-light half-radii (r_e ~ 1 kpc). However, recently it has been argued that not all dwarf galaxies are dominated by dark matter. Here we report 19 dwarf galaxies that could consist mainly of baryons up to radii well beyond r_e, at which point they are expected to be dominated by dark matter. Of these, 14 are isolated dwarf galaxies, free from the influence of nearby bright galaxies and high dense environments. This result provides observational evidence that could challenge the formation theory of low-mass galaxies within the framework of standard cosmology. Further observations, in particular deep imaging and spatially-resolved kinematics, are needed to constrain the baryon fraction better in such galaxies.

We present the largest-ever sample of 79 Ly$\alpha$ emitters (LAEs) at $z\sim$ 7.0 selected in the COSMOS and CDFS fields of the LAGER project (the Lyman Alpha Galaxies in the Epoch of Reionization). Our newly amassed ultradeep narrowband exposure and deeper/wider broadband images have more than doubled the number of LAEs in COSMOS, and we have selected 30 LAEs in the second field CDFS. We detect two large-scale LAE-overdense regions in the COSMOS that are likely protoclusters at the highest redshift to date. We perform injection and recovery simulations to derive the sample incompleteness. We show significant incompleteness comes from blending with foreground sources, which however has not been corrected in LAE luminosity functions in {the} literature. The bright end bump in the Ly$\alpha$ luminosity function in COSMOS is confirmed with 6 (2 newly selected) luminous LAEs (L$_{Ly\alpha}$ $>$ 10$^{43.3}$ erg s$^{-1}$). Interestingly, the bump is absent in CDFS, in which only one luminous LAE is detected. Meanwhile, the faint end luminosity functions from the two fields well agree with each other. The 6 luminous LAEs in COSMOS coincide with 2 LAE-overdense regions, while such regions are not seen in CDFS. The bright-end luminosity function bump could be attributed to ionized bubbles in a patchy reionization. It appears associated with cosmic overdensities, thus supports an inside-out reionization topology at $z$ $\sim$ 7.0, i.e., the high density peaks were ionized earlier compared to the voids. An average neutral hydrogen fraction of $x_{HI}$ $\sim$ 0.2 -- 0.4 is derived at $z\sim$ 7.0 based on the cosmic evolution of the Ly$\alpha$ luminosity function.

According to the recycling model, neutron stars in low-mass X-ray binaries were spun up to millisecond pulsars (MSPs), which indicates that all MSPs in the Galactic plane ought to be harbored in binaries. However, about $20\%$ Galactic field MSPs are found to be solitary. To interpret this problem, we assume that the accreting neutron star in binaries may collapse and become a strange star when it reaches some critical mass limit. Mass loss and a weak kick induced by asymmetric collapse during the phase transition (PT) from neutron star to strange star can result in isolated MSPs. In this work, we use a population-synthesis code to examine the PT model. The simulated results show that a kick velocity of $\sim60~{\rm km~s}^{-1}$ can produce $\sim6\times10^3$ isolated MSPs and birth rate of $\sim6.6\times10^{-7} {\rm ~yr}^{-1}$ in the Galaxy, which is approximately in agreement with predictions from observations. For the purpose of comparisons with future observation, we also give the mass distributions of radio and X-ray binary MSPs, along with the delay time distribution.

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.

We present ALMA observations of the [CII] 158 $\mu$m fine structure line and dust continuum emission from two quasars, SDSS J104433.04-012502.2 and SDSS J012958.51-003539.7, at z=5.78. The ALMA observations at 0.2'' resolution map the dust and gas on kpc scales. The spatially resolved emission shows a similar trend of decreasing [CII]-FIR ratios with increasing FIR surface brightnesses as was found in the infrared luminous galaxies with intense star formation. We confirm the velocity gradients of [CII] emission found previously in SDSS J0129-0035. No clear evidence of order motion is detected in SDSS J1044-0125. The velocity maps and PV-diagrams also suggest turbulent gas clumps in both objects. We tentatively detect a [CII] peak offset 4.9 kpc to the East of SDSS J1044-0125. This may be associated with an infalling companion, or node of gas outflow. All these results suggest significant dynamical evolution of the ISM in the nuclear region of these young quasar-starburst systems. We fit the velocity map of the [CII] emission from SDSS J0129-0035 with a rotating disk model. The result suggests a face-on system with an inclination angle of 16 +/- 20 degree and constrains the lower limit of the host galaxy dynamical mass to be 2.6x10^10 Msun within the [CII] emitting region. It is likely that SDSS J0129-0035, as well as other young quasars with super massive black hole masses on order of 10^7 to 10^8 Msun, falls close to the black hole and host galaxy mass relation defined by local galaxies.

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

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

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

The width of the broad emission lines in quasars is commonly characterized either by the full-width-at-half-maximum (FWHM) or the square root of the second moment of the line profile ($\sigma_{\rm line}$), and used as an indicator of the virial velocity of the broad-line region (BLR) in the estimation of black hole (BH) mass. We measure FWHM and $\sigma_{\rm line}$ for H$\alpha$, H$\beta$ and Mg II broad lines in both the mean and root-mean-square (rms) spectra of a large sample of quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. We introduce a new quantitative recipe to measure $\sigma_{\rm line}$ that is reproducible, less susceptible to noise and blending in the wings, and scales with the intrinsic width of the line. We compare the four definitions of line width (FWHM and $\sigma_{\rm line}$ in mean and rms spectra, respectively) for each of the three broad lines and among different lines. There are strong correlations among different width definitions for each line, providing justification for using the line width measured in single-epoch spectroscopy as a virial velocity indicator. There are also strong correlations among different lines, suggesting alternative lines to H$\beta$ can be used to estimate virial BH masses. We further investigate the correlations between virial BH masses using different line width definitions and the stellar velocity dispersion of the host galaxies, and the dependence of line shape (characterized by the ratio FWHM/$\sigma_{\rm line}$) on physical properties of the quasar. Our results provide further evidence that FWHM is more sensitive to the orientation of a flattened BLR geometry than $\sigma_{\rm line}$, but the overall comparison between the virial BH mass and host stellar velocity dispersion does not provide conclusive evidence that one particular width definition is significantly better than the others.

We report on the highest energy photons from the Crab Nebula observed by the Tibet air shower array with the underground water-Cherenkov-type muon detector array. Based on the criterion of muon number measured in an air shower, we successfully suppress 99.92% of the cosmic-ray background events with energies $E>100$ TeV. As a result, we observed 24 photon-like events with $E>100$ TeV against 5.5 background events, which corresponds to 5.6$\sigma$ statistical significance. This is the first detection of photons with $E>100$ TeV from an astrophysical source.

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.

We present a detailed characterization of the 849 broad-line quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. Our quasar sample covers a redshift range of 0.1<z<4.5 and is flux-limited to i_PSF<21.7 without any other cuts on quasar properties. The main sample characterization includes: 1) spectral measurements of the continuum and broad emission lines for individual objects from the coadded first-season spectroscopy in 2014; 2) identification of broad and narrow absorption lines in the spectra; 3) optical variability properties for continuum and broad lines from multi-epoch spectroscopy. We provide improved systemic redshift estimates for all quasars, and demonstrate the effects of signal-to-noise ratio on the spectral measurements. We compile measured properties for all 849 quasars along with supplemental multi-wavelength data for subsets of our sample from other surveys. The SDSS-RM sample probes a diverse range in quasar properties, and shows well detected continuum and broad-line variability for many objects from first-season monitoring data. The compiled properties serve as the benchmark for follow-up work based on SDSS-RM data. The spectral fitting tools are made public along with this work.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

We present the design for the first narrowband filter NB964 for the Dark Energy Camera (DECam), which is operated on the 4m Blanco Telescope at the Cerro Tololo Inter-American Observatory. The NB964 filter profile is essentially defined by maximizing the power of searching for Lyman alpha emitting galaxies (LAEs) in the epoch of reionization, with the consideration of the night sky background in the near-infrared and the DECam quantum efficiency. The NB964 filter was manufactured by Materion in 2015. It has a central wavelength of 964.2 nm and a full width at half maximum (FWHM) of 9.2 nm. An NB964 survey named LAGER (Lyman Alpha Galaxies in the Epoch of Reionization) has been ongoing since December 2015. Here we report results of lab tests, on-site tests and observations with the NB964 filter. The excellent performances of this filter ensure that the LAGER project is able to detect LAEs at z~7 with a high efficiency.

The ubiquitous variability of quasars across a wide range of wavelengths and timescales encodes critical information about the structure and dynamics of the circumnuclear emitting regions that are too small to be directly resolved, as well as detailed underlying physics of accretion and feedback processes in these active supermassive black holes. We emphasize the importance of studying quasar variability with time-domain spectroscopy, focusing on two science cases: (1) reverberation mapping (RM) to measure the broad-line region sizes and black hole masses in distant quasars; (2) spectroscopic follow-up of extreme variability quasars that dramatically change their continuum and broad-line flux within several years. We highlight the need for dedicated optical-infrared spectroscopic survey facilities in the coming decades to accompany wide-area time-domain imaging surveys, including: (1) the next phase of the Sloan Digital Sky Survey (SDSS-V; ~2020-2025), an all-sky, time-domain multi-object spectroscopic survey with 2.5m-class telescopes; (2) the planned Maunakea Spectroscopic Explorer, a dedicated 10m-class spectroscopic survey telescope with a 1.5 sq. deg field-of-view and multiplex of thousands of fibers in both optical and near-IR (J+H) to begin operations in 2029; (3) the Time-domain Spectroscopic Observatory (TSO), a proposed Probe-class ~1.3m telescope at L2, with imaging and spectroscopy (R=200, 1800) in 4 bands (0.3 - 5 micron) and rapid slew capability to 90% of sky, which will extend the coverage of Hbeta to z=8.

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 present Karl G. Jansky Very Large Array (VLA) observations of the CO (2$-$1) line emission towards three far-infrared luminous quasars at $z\sim6$: SDSS J231038.88$+$185519.7 and SDSS J012958.51$-$003539.7 with $\sim0\farcs6$ resolution and SDSS J205406.42$-$000514.8 with $\sim2\farcs1$ resolution. All three sources are detected in the CO (2$-$1) line emission -- one source is marginally resolved, and the other two appear as point sources. Measurements of the CO (2$-$1) line emission allow us to calculate the molecular gas mass even without a CO excitation model. The inferred molecular gas masses are (0.8$-$4.3) $\times$ 10$^{10}$ $M_{\odot}$. The widths and redshifts derived from the CO (2$-$1) line are consistent with previous CO (6$-$5) and [\ion{C}{2}] measurements. We also report continuum measurements using the Herschel for SDSS J231038.88$+$185519.7 and SDSS J012958.51$-$003539.7, and for SDSS J231038.88+185519.7, data obtained at $\sim140$ and $\sim300$ GHz using the Atacama Large Millimeter/submillimeter Array (ALMA). In the case of SDSS J231038.88+185519.7, we present a detailed analysis of the spectral energy distribution and derive the dust temperature ($\sim40$ K), the dust mass ($\sim10^{9}$ $M_{\odot}$), the far-infrared luminosity (8$-$1000 $\mu$m; $\sim10^{13}$ $ L_{\odot}$) and the star formation rate (2400$-$2700 $M_{\odot}$ yr$^{-1}$). Finally, an analysis of the photo-dissociation regions associated with the three high redshift quasars indicates that the interstellar medium in these sources has similar properties to local starburst galaxies.

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

Enormous Ly$\alpha$ Nebulae (ELANe), unique tracers of galaxy density peaks, are predicted to lie at the nodes and intersections of cosmic filamentary structures. Previous successful searches for ELANe have focused on wide-field narrowband surveys, or have targeted known sources such as ultraluminous quasi-stellar-objects (QSOs) or radio galaxies. Utilizing groups of coherently strong Ly$\alpha$ absorptions (CoSLAs), we have developed a new method to identify high-redshift galaxy overdensities and have identified an extremely massive overdensity, BOSS1441, at $z=2-3$ (Cai et al. 2016a). In its density peak, we discover an ELAN that is associated with a relatively faint continuum. To date, this object has the highest diffuse Ly$\alpha$ nebular luminosity of $L_{\rm{nebula}}=5.1\pm0.1\times10^{44}$ erg s$^{-1}$. Above the 2$\sigma$ surface brightness limit of SB$_{\rm{Ly\alpha}}= 4.8\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$, this nebula has an end-to-end spatial extent of 442 kpc. This radio-quiet source also has extended \civ\ $\lambda1549$ and \heii\ $\lambda1640$ emission on $\gtrsim30$ kpc scales. Note that the Ly$\alpha$, \heii\ and \civ\ emission all have double-peaked line profiles. Each velocity component has a full-width-half-maximum (FWHM) of $\approx700 - 1000$ km s$^{-1}$. We argue that this Ly$\alpha$ nebula could be powered by shocks due to an AGN-driven outflow or/and photoionization by a strongly obscured source.

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

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.

We present reverberation mapping results from the first year of combined spectroscopic and photometric observations of the Sloan Digital Sky Survey Reverberation Mapping Project. We successfully recover reverberation time delays between the $g+i$-band emission and the broad H$\beta$ emission line for a total of 44 quasars, and for the broad H$\alpha$ emission line in 18 quasars. Time delays are computed using the JAVELIN and CREAM software and the traditional interpolated cross-correlation function (ICCF): Using well defined criteria, we report measurements of 32 H$\beta$ and 13 H$\alpha$ lags with JAVELIN, 42 H$\beta$ and 17 H$\alpha$ lags with CREAM, and 16 H$\beta$ and 8 H$\alpha$ lags with the ICCF. Lag values are generally consistent among the three methods, though we typically measure smaller uncertainties with JAVELIN and CREAM than with the ICCF, given the more physically motivated light curve interpolation and more robust statistical modeling of the former two methods. The median redshift of our H$\beta$-detected sample of quasars is 0.53, significantly higher than that of the previous reverberation-mapping sample. We find that in most objects, the time delay of the H$\alpha$ emission is consistent with or slightly longer than that of H$\beta$. We measure black hole masses using our measured time delays and line widths for these quasars. These black hole mass measurements are mostly consistent with expectations based on the local M-sigma relationship, and are also consistent with single-epoch black hole mass measurements. This work increases the current sample size of reverberation-mapped active galaxies by about two-thirds and represents the first large sample of reverberation mapping observations beyond the local universe (z < 0.3).

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

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

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.

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.

We present the rest-frame UV and optical photometry and morphology of low-redshift broad-line quasar host galaxies from the Sloan Digital Sky Survey Reverberation Mapping project. Our sample consists of 103 quasars at $z<0.8$, spanning a luminosity range of $-25\le M_g\le -17$ mag. We stack the multi-epoch images in the $g$ and $i$ bands taken by the Canada-France-Hawaii Telescope. The combined $g$-band ($i$-band) images reach a $5\sigma$ depth of 26.2 (25.2) mag, with a typical PSF size of $0.7"$ ($0.6"$). Each quasar is decomposed into a PSF and a S\'ersic profile, representing the central AGN and the host galaxy components, respectively. The systematic errors of the measured host galaxy flux in the two bands are 0.23 and 0.18 mag. The relative errors of the measured galaxy half-light radii ($R_e$) are about 13%. We estimate the rest-frame $u$- and $g$-band flux of the host galaxies, and find that the AGN-to-galaxy flux ratios in the $g$ band are between 0.9 to 4.4 (68.3% confidence). These galaxies have high stellar masses $M_\ast = 10^{10}\sim10^{11}\, M_\odot$. They have similar color with star-forming galaxies at similar redshifts, in consistent with AGN positive feedback in these quasars. We find that the $M_*-M_\text{BH}$ relation in our sample is shallower than the local $M_\text{Bulge}-M_\text{BH}$ relation. The S\'ersic indices and the $M_*-R_e$ relation indicate that the majority of the host galaxies are disk-like.

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 examine the possible influence of Earth-directed coronal mass ejections (ECMEs) on the Sun's shadow in the 3~TeV cosmic-ray intensity observed by the Tibet-III air shower (AS) array. We confirm a clear solar-cycle variation of the intensity deficit in the Sun's shadow during ten years between 2000 and 2009. This solar-cycle variation is overall reproduced by our Monte Carlo (MC) simulations of the Sun's shadow based on the potential field model of the solar magnetic field averaged over each solar rotation period. We find, however, that the magnitude of the observed intensity deficit in the Sun's shadow is significantly less than that predicted by MC simulations, particularly during the period around solar maximum when a significant number of ECMEs is recorded. The $\chi^2$ tests of the agreement between the observations and the MC simulations show that the difference is larger during the periods when the ECMEs occur, and the difference is reduced if the periods of ECMEs are excluded from the analysis. This suggests the first experimental evidence of the ECMEs affecting the Sun's shadow observed in the 3~TeV cosmic-ray intensity.

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.

We present a detailed study of a unusually bright, lensed galaxy at $z=5.424$ discovered within the CFHTLS imaging survey. With an observed flux of $i_{\rm AB}=23.0$, J141446.82+544631.9 is one of the brightest galaxies known at $z>5$. It is characterized by strong Ly$\alpha$ emission, reaching a peak in (observed) flux density of $>10^{-16}~{\rm erg}~{\rm s}^{-1}~{\rm cm}^{-2}~{\rm \AA}^{-1}$. A deep optical spectrum from the LBT places strong constraints on NV and CIV emission, disfavouring an AGN source for the emission. However, a detection of the NIV]~$\lambda$1486 emission line indicates a hard ionizing continuum, possibly from hot, massive stars. Resolved imaging from HST deblends the galaxy from a foreground interloper; these observations include narrowband imaging of the Ly$\alpha$ emission, which is marginally resolved on $\sim$few~kpc scales and has EW$_0~\sim$ 260\AA. The Ly$\alpha$ emission extends over ~2000~${\rm km}~{\rm s}^{-1}$ and is broadly consistent with expanding shell models. SED fitting that includes Spitzer/IRAC photometry suggests a complex star formation history that include both a recent burst and an evolved population. J1414+5446 lies 30" from the centre of a known lensing cluster in the CFHTLS; combined with the foreground contribution this leads to a highly uncertain estimate for the lensing magnification in the range $5 \leq \mu \leq 25$. Because of its unusual brightness J1414+5446 affords unique opportunities for detailed study of an individual galaxy near the epoch of reionization and a preview of what can be expected from upcoming wide-area surveys that will yield hundreds of similar objects.

We present measurements of the mean and scatter of the IGM Lyman-{\alpha} opacity at 4.9 < z < 6.1 along the lines of sight of 62 quasars at z > 5.7, the largest sample assembled at these redshifts to date by a factor of two. The sample size enables us to sample cosmic variance at these redshifts more robustly than ever before. The spectra used here were obtained by the SDSS, DES-VHS and SHELLQs collaborations, drawn from the ESI and X-Shooter archives, reused from previous studies or observed specifically for this work. We measure the effective optical depth of Lyman-{\alpha} in bins of 10, 30, 50 and 70 cMpc h-1, construct cumulative distribution functions under two treatments of upper limits on flux and explore an empirical analytic fit to residual Lyman-{\alpha} transmission. We verify the consistency of our results with those of previous studies via bootstrap re-sampling and confirm the existence of tails towards high values in the opacity distributions, which may persist down to z = 5.2. Comparing our results with predictions from cosmological simulations, we find further strong evidence against models that include a spatially uniform ionizing background and temperature-density relation. We also compare to IGM models that include either a fluctuating UVB dominated by rare quasars or temperature fluctuations due to patchy reionization. Although both models produce better agreement with the observations, neither fully captures the observed scatter in IGM opacity. Our sample of 62 z > 5.7 quasar spectra opens many avenues for future study of the reionisation epoch.

A small fraction of quasars show an unusually high nitrogen-to-carbon ratio (N/C) in their spectra. These "nitrogen-rich" (N-rich) quasars are a long-standing puzzle because their interstellar medium implies stellar populations with abnormally high metallicities. It has recently been proposed that N-rich quasars may result from tidal disruption events (TDEs) of stars by supermassive black holes. The rapid enhancement of nitrogen and the depletion of carbon due to the carbon--nitrogen--oxygen cycle in supersolar mass stars could naturally produce high N/C. However, the TDE hypothesis predicts that the N/C should change with time, which has never hitherto been observed. Here we report the discovery of the first N-rich quasar with rapid N/C variability that could be caused by a TDE. Two spectra separated by 1.7 years (rest-frame) show that the N III]\lambda 1750/C III]\lambda 1909 intensity ratio decayed by ~86%+/-14% (1 \sigma). Optical (rest-frame UV) light-curve and X-ray observations are qualitatively consistent with the TDE hypothesis; though, the time baseline falls short of a definitive proof. Putting the single-object discovery into context, statistical analyses of the ~80 known N-rich quasars with high-quality archival spectra show evidence (at a 5\sigma\ significance level) of a decrease in N/C on timescales of >1 year (rest-frame) and a constant level of ionization (indicated by the C III]\lambda 1909/C IV\lambda 1549 intensity ratio). If confirmed, our results demonstrate the method of identifying TDE candidates in quasars via abundance ratio variability, opening a new window of TDE observations at high redshift (z>2) with upcoming large-scale time-domain spectroscopic surveys.

The Hubble Space Telescope observation of the double white dwarf (WD) binary SDSS J125733.63$+$542850.5 reveal that the massive WD has a surface gravity log$g_1\sim8.7$ (which implies a mass of $M_1\sim1.06~{\rm M_\odot}$) and an effective temperature $T_1\sim13000$ K, while the effective temperature of the low-mass WD ($M_2<0.24$ {\rm M$_\odot$}) is $T_2\sim6400K$. Therefore, the massive and the low-mass WDs have a cooling age $\tau_1\sim1$ {\rm Gyr} and $\tau_2\geq5$ {\rm Gyr}, respectively. This is in contradiction with traditional binary evolution theory. In this Letter, we propose a strange dwarf scenario to explain the formation of this double WD binary. We assume that the massive WD is a strange dwarf originating from a phase transition in a $\sim1.1$ M$_\odot$ WD, which has experienced accretion and spin-down processes. Its high effective temperature could arise from the heating process during the phase transition. Our simulations suggest that the progenitor of SDSS J125733.63$+$542850.5 can be a binary system consisting of a $0.65~\rm M_{\odot}$ WD and a $1.5~\rm M_{\odot}$ main sequence star in a 1.492 day orbit. Especially, the secondary star (i.e., the progenitor of the low mass WD) is likely to have an ultra-low metallicity of $Z=0.0001$.

Quasars (QSOs) hosting supermassive black holes are believed to reside in massive halos harboring galaxy overdensities. However, many observations revealed average or low galaxy densities around $z\gtrsim6$ QSOs. This could be partly because they measured galaxy densities in only tens of arcmin$^2$ around QSOs and might have overlooked potential larger scale galaxy overdensities. Some previous studies also observed only Lyman break galaxies (LBGs, massive older galaxies) and missed low mass young galaxies like Ly$\alpha$ emitters (LAEs) around QSOs. Here we present observations of LAE and LBG candidates in $\sim700$ arcmin$^2$ around a $z=6.61$ luminous QSO using Subaru Telescope Suprime-Cam with narrow/broadbands. We compare their sky distributions, number densities and angular correlation functions with those of LAEs/LBGs detected in the same manner and comparable data quality in our control blank field. In the QSO field, LAEs and LBGs are clustering in 4-20 comoving Mpc angular scales, but LAEs show mostly underdensity over the field while LBGs are forming $30\times60$ comoving Mpc$^2$ large scale structure containing 3-$7\sigma$ high density clumps. The highest density clump includes a bright (23.78 mag in the narrowband) extended ($\gtrsim 16$ kpc) Ly$\alpha$ blob candidate, indicative of a dense environment. The QSO could be part of the structure but is not located exactly at any of the high density peaks. Near the QSO, LAEs show underdensity while LBGs average to $4\sigma$ excess densities compared to the control field. If these environments reflect halo mass, the QSO may not be in the most massive halo, but still in a moderately massive one.

We analyze the Sun's shadow observed with the Tibet-III air shower array and find that the shadow's center deviates northward (southward) from the optical solar disc center in the "Away" ("Toward") IMF sector. By comparing with numerical simulations based on the solar magnetic field model, we find that the average IMF strength in the "Away" ("Toward") sector is $1.54 \pm 0.21_{\rm stat} \pm 0.20_{\rm syst}$ ($1.62 \pm 0.15_{\rm stat} \pm 0.22_{\rm syst}$) times larger than the model prediction. These demonstrate that the observed Sun's shadow is a useful tool for the quantitative evaluation of the average solar magnetic field.

We present observations of escaping Lyman Continuum (LyC) radiation from 34 massive star-forming galaxies and 12 weak AGN with reliably measured spectroscopic redshifts at $z$$\simeq$2.3-4.1. We analyzed Hubble Space Telescope Wide Field Camera 3 (WFC3) mosaics of the Early Release Science field in three UVIS filters to sample the rest-frame LyC over this redshift range. With our best current assessment of the WFC3 systematics, we provide 1$\sigma$ upper limits for the average LyC emission of galaxies at $z$=2.35, 2.75, and 3.60 to $\sim$28.5, 28.1, and 30.7 mag in image stacks of 11-15 galaxies in the WFC3/UVIS F225W, F275W, and F336W, respectively. The LyC flux of weak AGN at $z$=2.62 and 3.32 are detected at 28.3 and 27.4 mag with SNRs of $\sim$2.7 and 2.5 in F275W and F336W for stacks of 7 and 3 AGN, respectively, while AGN at $z$=2.37 are constrained to $\gtrsim$27.9 mag at 1$\sigma$ in a stack of 2 AGN. The stacked AGN LyC light profiles are flatter than their corresponding non-ionizing UV continuum profiles out to radii of r$\lesssim$0."9, which may indicate a radial dependence of porosity in the ISM. With synthetic stellar SEDs fit to UV continuum measurements longwards of Ly$\alpha$ and IGM transmission models, we constrain the absolute LyC escape fractions to $f_{\rm esc}^{\rm abs}$$\simeq$$22^{+44}_{-22}$% at $z$=2.35 and $\lesssim$55% at $z$=2.75 and 3.60, respectively. All available data for galaxies, including published work, suggests a more sudden increase of $f_{\rm esc}$ with redshift at $z$$\simeq$2. Dust accumulating in (massive) galaxies over cosmic time correlates with increased HI column density, which may lead to reducing $f_{\rm esc}$ more suddenly at $z$$\lesssim$2. This may suggest that star-forming galaxies collectively contributed to maintaining cosmic reionization at redshifts $z$$\gtrsim$2-4, while AGN likely dominated reionization at $z$$\lesssim$2.

Studies of the most luminous quasars at high redshift directly probe the evolution of the most massive black holes in the early Universe and their connection to massive galaxy formation. However, extremely luminous quasars at high redshift are very rare objects. Only wide area surveys have a chance to constrain their population. The Sloan Digital Sky Survey (SDSS) has so far provided the most widely adopted measurements of the quasar luminosity function (QLF) at $z>3$. However, a careful re-examination of the SDSS quasar sample revealed that the SDSS quasar selection is in fact missing a significant fraction of $z\gtrsim3$ quasars at the brightest end. We have identified the purely optical color selection of SDSS, where quasars at these redshifts are strongly contaminated by late-type dwarfs, and the spectroscopic incompleteness of the SDSS footprint as the main reasons. Therefore we have designed the Extremely Luminous Quasar Survey (ELQS), based on a novel near-infrared JKW2 color cut using WISE AllWISE and 2MASS all-sky photometry, to yield high completeness for very bright ($m_{\rm{i}} < 18.0$) quasars in the redshift range of $3.0\leq z\leq5.0$. It effectively uses random forest machine-learning algorithms on SDSS and WISE photometry for quasar-star classification and photometric redshift estimation. The ELQS will spectroscopically follow-up $\sim 230$ new quasar candidates in an area of $\sim12000\,\rm{deg}^2$ in the SDSS footprint, to obtain a well-defined and complete quasars sample for an accurate measurement of the bright-end quasar luminosity function at $3.0\leq z\leq5.0$. In this paper we present the quasar selection algorithm and the quasar candidate catalog.

New Multi-color $B$ $V$ $R_c$ $I_c$ photometric observation are presented for W UMa type eclipsing binary V1073 Cyg. The multi-color light curves analysis with the Wilson-Devinney(W-D) procedure acquired the absolute parameters of this system, showing that V1073Cyg is a shallow contact binary system with fill-out factor $f=0.124(\pm0.011)$. We collected all available times of light minima spanning 119 years including CCD data to construct the O-C curve and made detailed O-C analysis. The O-C diagram shows that the period change is complex. There exist a long-term continuous decrease and a cyclic variation. The period is decreasing at a rate of $\dot P=-1.04(\pm0.18)\times 10^{-10}$ $days\cdot{cycle}^{-1}$, and with the period decrease, V1073 Cyg will evolve to deep contact stage. The cyclic variation with a period of $P_3=82.7(\pm3.6) years$ and an amplitude of $A=0.028(\pm0.002) day$ may be explained by the magnetic activity of one or both components or the light travel time effect (LTTE) caused by a distant third companion with $M_3({i'}=90^{\circ})=0.511M_\odot$.

Aims. We aim to find more eclipsing multiple systems and obtain their parameters, thus increasing our understanding of multiple systems. Methods. The extraneous eclipses on the \textit{kepler} binary light curves indicating extraneous bodies were searched. The binary light curves were analyzed using the binary model, and the extraneous eclipses were studied on their periodicity and shape changes. Results. Three binaries with extraneous eclipses on the binary light curves were found and studied based on the \textit{Kepler} observations. The object KIC 5255552 is an eclipsing triple system with a fast changing inner binary and an outer companion uncovered by three groups of extraneous eclipses of $862.1(\pm0.1)$ d period. The KIC 10091110 is suggested to be a double eclipsing binary system with several possible extraordinary coincidences: the two binaries share similar extremely small mass ratios ($0.060(13)$ and $0.0564(18)$), similar mean primary densities ($0.3264(42)\;\rho_\odot$ and $0.3019(28)\;\rho_\odot$), and, most notably, the ratio of the two binaries' periods is very close to integer 2 (8.5303353/4.2185174 = 2.022). The KIC 11495766 is a probable triple system with a $\sim120.73$ d period binary and (at least) one non-eclipse companion. Furthermore, very close to it in the celestial sphere, there is a blended background stellar binary of 8.3404432 d period. A first list of 25 eclipsing multiple candidates is presented, with the hope that it will be beneficial for study of eclipsing multiples.

We presents results from a spectroscopic survey of $z\sim5$ quasars in the CFHT Legacy Survey (CFHTLS). Using both optical color selection and a likelihood method we select 97 candidates over an area of 105 deg$^2$ to a limit of $i_{\rm AB} < 23.2$, and 7 candidates in the range $23.2 < i_{\rm AB} < 23.7$ over an area of 18.5 deg$^2$. Spectroscopic observations for 43 candidates were obtained with Gemini, MMT, and LBT, of which 37 are $z>4$ quasars. This sample extends measurements of the quasar luminosity function $\sim$1.5 mag fainter than our previous work in SDSS Stripe 82. The resulting luminosity function is in good agreement with our previous results, and suggests that the faint end slope is not steep. We perform a detailed examination of our survey completeness, particularly the impact of the Ly$\alpha$ emission assumed in our quasar spectral models, and find hints that the observed Ly$\alpha$ emission from faint $z\sim5$ quasars is weaker than for $z\sim3$ quasars at a similar luminosity. Our results strongly disfavor a significant contribution of faint quasars to the hydrogen-ionizing background at $z=5$.

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

The SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES) is a three-year JCMT Large Program aiming at reaching the 450 $\mu$m confusion limit in the COSMOS-CANDELS region, to study a representative sample of the high-redshift far-infrared galaxy population that gives rise to the bulk of the far-infrared background. We present the first-year data from STUDIES. We have reached a 450 $\mu$m noise level of 0.91~mJy for point sources at the map center, covered an area of 151 arcmin$^2$, and detected 98 and 141 sources at 4.0 and 3.5 $\sigma$, respectively. Our derived counts are best constrained in the 3.5-25 mJy regime using directly detected sources. Below the detection limits, our fluctuation analysis further constrains the slope of the counts down to 1 mJy. The resulting counts at 1-25 mJy are consistent with a power law having a slope of $-2.59$ ($\pm0.10$ for 3.5-25 mJy, and $^{+0.4}_{-0.7}$ for 1-3.5 mJy). There is no evidence of a faint-end termination or turn-over of the counts in this flux density range. Our counts are also consistent with previous SCUBA-2 blank-field and lensing cluster surveys. The integrated surface brightness from our counts down to 1 mJy is $90.0\pm17.2$ Jy deg$^{-2}$, which can account for up to $83^{+15}_{-16}\%$ of the COBE 450 $\mu$m background. We show that Herschel counts at 350 and 500 $\mu$m are significantly higher than our 450 $\mu$m counts, likely caused by its large beam and source clustering. High-angular resolution instruments like SCUBA-2 at 450 $\mu$m are therefore highly beneficial for measuring the luminosity and spatial density of high-redshift dusty galaxies.

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.

Two twin binaries, KIC 4826439 and KIC 6045264, with very similar component stars were found photometrically based on $\textit{Kepler}$ eclipsing binary light curves. The absolute parameters of the massive components are 1.156(0.03)$M_\odot$, 1.881(0.02)$R_\odot$, 6065K for KIC 4826439, and 0.874(0.3)$M_\odot$, 1.206(0.02)$R_\odot$, 6169(30)K for KIC 6045264. The differences between the components are less than two percents for all the parameters. A very low proportion of the twin binaries ($2/1592\approx0.13\%$) was found, which does not support the previous findings of the excesses of twins on binary mass ratio distribution, but support a deficiently low proportion of twins. A new method is practiced to work out the absolute parameters of the two twins without the radial velocities. This method requires the solution of the light curves, the spectra and the evolutionary isochrones of covering the complete stellar parameter space, simultaneously. We also studied their evolution tracks that: KIC 4826439 will experience an unstable mass transfer stage followed by an unclear ending, and KIC 6045264 will become a single star via an over-contact phase. It seems highly unlikely that the two twin binaries will produce twin degenerate binaries, although they have quite similar components.

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

Narrowband imaging is a highly successful approach for finding large numbers of high redshift Lya emitting galaxies (LAEs) up to z~6.6. However, at z>~7 there are as yet only 3 narrowband selected LAEs with spectroscopic confirmations (two at z~6.9-7.0, one at z~7.3), which hinders extensive studies on cosmic reionization and galaxy evolution at this key epoch. We have selected 23 candidate z~6.9 LAEs in COSMOS field with the large area narrowband survey LAGER (Lyman-Alpha Galaxies at the End of Reionization). In this work we present spectroscopic followup observations of 12 candidates using IMACS on Magellan. For 9 of these, the observations are sufficiently deep to detect the expected lines. Lya emission lines are identified in six sources (yielding a success rate of 2/3), including 3 luminous LAEs with Lya luminosities of L(Lya) ~ 10^{43.5} erg/s, the highest among known spectroscopically confirmed galaxies at >~7.0. This triples the sample size of spectroscopically confirmed narrowband selected LAEs at z>~7, and confirms the bright end bump in the Lya luminosity function we previously derived based on the photometric sample, supporting a patchy reionization scenario. Two luminous LAEs appear physically linked with projected distance of 1.1 pMpc and velocity difference of ~ 170 km/s. They likely sit in a common ionized bubble produced by themselves or with close neighbors, which reduces the IGM attenuation of Lya. A tentative narrow NV${\lambda}$1240 line is seen in one source, hinting at activity of a central massive black hole with metal rich line emitting gas.

We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations of the dust continuum and [C II] 158 $\mu$m fine structure line emission towards a far-infrared-luminous quasar, ULAS J131911.29$+$095051.4 at $z=6.13$, and combine the new Cycle 1 data with ALMA Cycle 0 data. The combined data have an angular resolution $\sim$ $0.3$, and resolve both the dust continuum and the [C II] line emission on few kpc scales. The [C II] line emission is more irregular than the dust continuum emission which suggests different distributions between the dust and [C II]-emitting gas. The combined data confirm the [C II] velocity gradient that we previously detected in lower resolution ALMA image from Cycle 0 data alone. We apply a tilted ring model to the [C II] velocity map to obtain a rotation curve, and constrain the circular velocity to be 427 $\pm$ 55 km s$^{-1}$ at a radius of 3.2 kpc with an inclination angle of 34$^\circ$. We measure the dynamical mass within the 3.2 kpc region to be 13.4$_{-5.3}^{+7.8}$ $\times 10^{10}\,M_{\odot}$. This yields a black hole and host galaxy mass ratio of 0.020$_{-0.007}^{+0.013}$, which is about 4$_{-2}^{+3}$ times higher than the present-day $M_{\rm BH}$/$M_{\rm bulge}$ ratio. This suggests that the supermassive black hole grows the bulk of its mass before the formation of the most of stellar mass in this quasar host galaxy in the early universe.

Recent observations have suggested that the CIII]$\lambda1907/1909$ emission lines could be alternative diagnostic lines for galaxies in the reionization epoch. We use the F128N narrowband filter on the Hubble Space Telescope's ($\it{HST}$) Wide Field Camera 3 (WFC3) to search for CIII] emission in a sample of five galaxies at z = 5.7 in the Subaru Deep Field and the Subaru/XMM-Newton Deep Field. Using the F128N narrowband imaging, together with the broadband imaging, we do not detect CIII] emission for the five galaxies with $J_{\rm{AB}}$ ranging from 24.10 -- 27.00 in our sample. For the brightest galaxy J132416.13+274411.6 in our sample (z = 5.70, $J_{\rm{AB}} = 24.10$), which has a significantly higher signal to noise, we report a CIII] flux of $3.34\pm1.81 \times 10^{-18}$ $\mathrm{erg\ s^{-1}\ cm^{-2}}$, which places a stringent 3-$\rm\sigma$ upper limit of $5.43\times 10^{-18}$ $\mathrm{erg\ s^{-1}\ cm^{-2}}$ on CIII] flux and 6.57 \AA\ on the CIII] equivalent width. Using the stacked image, we put a 3-$\rm\sigma$ upper limit on the mean CIII] flux of $\mathrm{2.55\times10^{-18}\ erg\ s^{-1}\ cm^{-2}}$, and a 3-$\rm\sigma$ upper limit on the mean CIII] equivalent width of 4.20 {\AA} for this sample of galaxies at z = 5.70. Combined with strong CIII] detection reported among high-z galaxies in the literature, our observations suggest that the equivalent widths of CIII] from galaxies at z $>$ 5.70 exhibit a wide range of distribution. Our strong limits on CIII] emission could be used as a guide for future observations in the reionization epoch.

We present the first results from the ongoing LAGER project (Lyman Alpha Galaxies in the Epoch of Reionization), which is the largest narrowband survey for $z \sim$ 7 galaxies to date. Using a specially built narrowband filter NB964 for the superb large-area Dark-Energy Camera (DECam) on the NOAO/CTIO 4m Blanco telescope, LAGER has collected 34 hours NB964 narrowband imaging data in the 3 deg$^2$ COSMOS field. We have identified 23 Lyman Alpha Emitter (LAE) candidates at $z$ = 6.9 in the central 2-deg$^2$ region, where DECam and public COSMOS multi-band images exist. The resulting luminosity function can be described as a Schechter function modified by a significant excess at the bright end (4 galaxies with $L_{Ly\alpha} \sim $ 10$^{43.4\pm0.2}$ erg s$^{-1}$). The number density at $L_{Ly\alpha}\sim$ 10$^{43.4\pm0.2}$ erg s$^{-1}$ is little changed from z= 6.6, while at fainter $L_{Ly\alpha}$ it is substantially reduced. Overall, we see a fourfold reduction in Ly$\alpha$ luminosity density from $z$ = 5.7 to 6.9. Combined with a more modest evolution of the continuum UV luminosity density, this suggests a factor of $\sim 3$ suppression of Ly$\alpha$ by radiative transfer through the $z \sim$ 7 intergalactic medium (IGM). It indicates an IGM neutral fraction $x_{HI}$ $\sim$ 0.4--0.6 (assuming Ly$\alpha$ velocity offsets of 100-200 km s$^{-1}$). The changing shape of the Ly$\alpha$ luminosity function between $z\lesssim 6.6$ and $z=6.9$ supports the hypothesis of ionized bubbles in a patchy reionization at $z\sim$ 7.

We use a pair of high resolution N-body simulations implementing two dark matter models, namely the standard cold dark matter (CDM) cosmogony and a warm dark matter (WDM) alternative where the dark matter particle is a 1.5keV thermal relic. We combine these simulations with the GALFORM semi-analytical galaxy formation model in order to explore differences between the resulting galaxy populations. We use GALFORM model variants for CDM and WDM that result in the same z=0 galaxy stellar mass function by construction. We find that most of the studied galaxy properties have the same values in these two models, indicating that both dark matter scenarios match current observational data equally well. Even in under-dense regions, where discrepancies in structure formation between CDM and WDM are expected to be most pronounced, the galaxy properties are only slightly different. The only significant difference in the local universe we find is in the galaxy populations of "Local Volumes", regions of radius 1 to 8Mpc around simulated Milky Way analogues. In such regions our WDM model provides a better match to observed local galaxy number counts and is five times more likely than the CDM model to predict sub-regions within them that are as empty as the observed Local Void. Thus, a highly complete census of the Local Volume and future surveys of void regions could provide constraints on the nature of dark matter.

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.

We present the HST WFC3/F275W UV imaging observations of A2218-Flanking, a lensed compact dwarf galaxy at redshift $z\approx2.5$. The stellar mass of A2218-Flanking is $\log(M/M_{\odot}=9.14$ and SFR is 12.5$M_\odot$yr$^{-1}$ after correcting the magnification. This galaxy has a young galaxy age of 127 Myr and a compact galaxy size of $r_{1/2}=2.4$kpc. The HST UV imaging observations cover the rest-frame Lyman continuum (LyC) emission ($\sim800${\AA}) from A2218-Flanking. We firmly detect ($14\sigma$) the LyC emission in A2218-Flanking in the F275W image. Together with the HST F606W images, we find that the absolute escape fraction of LyC is $f_{\rm abs,esc}>28-57\%$ based on the flux density ratio between 1700{\AA} and 800{\AA} ($f_{1700}/f_{800}$). The morphology of the LyC emission in the F275W images is extended and well follows the morphology of the UV continuum morphology in the F606W images, suggesting that the $f_{800}$ is not from foreground contaminants. We find that the region with a high star formation rate surface density has a lower $f_{1700}/f_{800}$ (higher $f_{800}/f_{1700}$) ratio than the diffused regions, suggesting that LyC photons are more likely to escape from the region with the intensive star-forming process. We compare the properties of galaxies with and without LyC detections and find that LyC photons are easier to escape in low mass galaxies.

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

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

It is generally thought that an ultracompact X-ray binary is composed with a neutron star and a helium white dwarf donor star. As one of the most compact binaries, 4U 1820-303 in globular cluster NGC 6624 was predicted an orbital-period derivative at a rate of $\dot{P}/P\sim1.1\times10^{-7}$ ${\rm yr^{-1}}$ if the mass transfer is fully driven by gravitational radiation. However, the recent analysis of the 16 ${\rm yr}$ data from \textit{Rossi X-ray Timing Explorer} and other historical records yielded a negative orbital-period derivative in the past $35$ yr. In this work, we propose an evolutionary circumbinary (CB) disk model to account for this anomalous orbital-period derivative. 4U 1820-30 was known to undergo superbust events caused by runaway thermal nuclear burning on the neutron star. We assume that for a small fraction of the superbursts, part of the ejected material may form a CB disk around the binary. If the recurrence time of such superbursts is $\sim10,000$ yr and $\sim10$% of the ejected mass feeds a CB disk, the abrupt angular-momentum loss causes a temporary orbital shrink, and the donor's radius and its Roche-lobe radius do not keep in step. Driven by mass transfer and angular-momentum loss, the binary would adjust its orbital parameters to recover a new stable stage. Based on the theoretical analysis and numerical simulation, we find that the required feed mass at the CB disk is approximately $\sim 10^{-8}$ ${\rm M_{\odot}}$.

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.

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

Application of fitting techniques to obtain physical parameters---such as ages, metallicities, and $\alpha$-element to iron ratios---of stellar populations is an important approach to understand the nature of both galaxies and globular clusters (GCs). In fact, fitting methods based on different underlying models may yield different results, and with varying precision. In this paper, we have selected 22 confirmed M31 GCs for which we do not have access to previously known spectroscopic metallicities. Most are located at approximately one degree (in projection) from the galactic center. We performed spectroscopic observations with the 6.5 m MMT telescope, equipped with its Red Channel Spectrograph. Lick/IDS absorption-line indices, radial velocities, ages, and metallicities were derived based on the $\rm EZ\_Ages$ stellar population parameter calculator. We also applied full spectral fitting with the ULySS code to constrain the parameters of our sample star clusters. In addition, we performed $\chi^2_{\rm min}$ fitting of the clusters' Lick/IDS indices with different models, including the Bruzual & Charlot models (adopting Chabrier or Salpeter stellar initial mass functions and 1994 or 2000 Padova stellar evolutionary tracks), the GALEV, and the Thomas et al. models. For comparison, we collected their $UVBRIJK$ photometry from the Revised Bologna Catalogue (v.5) to obtain and fit the GCs' spectral energy distributions (SEDs). Finally, we performed fits using a combination of Lick/IDS indices and SEDs. The latter results are more reliable and the associated error bars become significantly smaller than those resulting from either our Lick/IDS indices-only or our SED-only fits.

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.

Cosmological simulations suggest a strong correlation between high optical-depth Ly$\alpha$ absorbers, which arise from the intergalactic medium (IGM), and 3-D mass overdensities on scales of $10-30$ $h^{-1}$ comoving Mpc. By examining the absorption spectra of $\sim$ 80,000 QSO sight-lines over a volume of 0.1 Gpc$^3$ in the Sloan Digital Sky Survey III (SDSS-III), we have identified an extreme overdensity, BOSS1441, which contains a rare group of strong Ly$\alpha$ absorbers at $z=2.32\pm 0.02$. This absorber group is associated with six QSOs at the same redshift on a 30 comoving Mpc scale. Using Mayall/MOSAIC narrowband and broadband imaging, we detect Ly$\alpha$ emitters (LAEs) down to $0.7\times L_{\rm{Ly\alpha}}^*$, and reveal a large-scale structure of Ly$\alpha$ emitters (LAEs) in this field. Our follow-up Large Binocular Telescope (LBT) observations have spectroscopically confirmed 19 galaxies in the density peak. We show that BOSS1441 has an LAE overdensity of $10.8\pm 2.6$ on a 15 comoving Mpc scale which could collapse to a massive cluster with $M\gtrsim10^{15}$ M$_\odot$ at $z\sim0$. This overdensity is among the most massive large-scale structures at $z\sim2$ discovered to date.

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 four-color ($B$ $V$ $R_c$ $I_c$) light curves of V776 Cas are presented and analyzed using the Wilson-Devinney (W-D) method. It is discovered that V776 Cas is an early F-type (F2V) overcontact binary with a very high contact degree ($ f=64.6\,\%$) and an extremely low mass ratio ($q=0.130$), which indicate that it is at the final evolutionary stage of cool short-period binaries. The mass of the primary and secondary stars are calculated to be $M_1 = 1.55(\pm0.04)M_\odot$, $M_2 = 0.20(\pm0.01)M_\odot$. V776 Cas is supposed to be formed from an initially detached binary system via the loss of angular momentum due to the magnetic wind. The initial mass of the present primary and secondary components are calculated to be $M_{1i} = 0.86(\pm0.10)M_\odot$ and $M_{2i} = 2.13(\pm0.04)M_\odot$. The observed-calculated ($O$-$C$) curve exhibits a cyclic period variation, which is due to the light-travel time effect (LTTE) caused by the presence of a third component with a period of 23.7 years. The mass of the third component is estimated to be $M_3 = 1.04(\pm0.03)M_\odot$ and the orbital inclination of the third component is calculated to be $i' = 33.1^{\circ}$. The distance of the binary system to the mass center of the triple system is calculated to be $a'_{12} = 3.45AU$. The presence of the close-in tertiary component may play an important role in the formation and evolution of this binary system by drawing angular momentum from the central system.

The first V Rc Ic bands light curves of MQ UMa are presented and analyzed using the Wilson-Devinney (W-D) program. It is discovered that MQ UMa is an A-subtype contact binary with a high fill-out (f=82 %) and a low mass ratio (q = 0:195), which indicates that it is in the late evolutionary stage of latetype tidal-locked binary stars. The mass of the primary and secondary stars are estimated and the evolutionary status of the two components are placed in the H-R diagram. The W-D solutions also indicates that MQ UMa may be a triple system with an F5V type additional component. A sample of sixteen high fillout, low mass ratio overcontact binaries are collected and their possible evolution scenarios are discussed. Using the five times of light minimum recently observed together with those collected from literatures, the authors find that the observedcalculate (O-C) curve exhibits a cyclic period variation. The cyclic period change also reveals the presence of a tertiary component, which may play an important role in the formation and evolution of this binary system by drawing angular momentum from the central system.

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 report a candidate centi-parsec supermassive black hole binary (SMBHB) in the radio-quiet quasar SDSS J0159+0105 at z=0.217. With a modified lomb-scargle code GLSdeDRW and the auto-correlation analysis ACF, we detect two significant (at P>99%) periodic signals at ~741 day and ~1500 day from the 8.1-year Catalina V-band light curve of this quasar. The period ratio, which is close to 1:2, is typical of a black-hole binary system with a mass ratio of 0.05<q<0.8 according to recent numerical simulations. SDSS J0159+0105 has two SDSS spectroscopic observations separated by ~10 years. There is a significant change in the broad H-beta profile between the two epochs, which can be explained by a single broad-line region (BLR) around the binary system illuminated by the aforementioned mini-disks, or a stream of gas flowing from the circumbinary disk to one of the SMBHs. From the single BLR assumption and the orbital period t_orb ~1500 day, we estimate the total virial masses of M_SMBHB ~ 1.3x10^8 M_sun, the average distances of BLR of ~0.04pc (~50 light-day, with +/-0.3 dex uncertainty), and a SMBHB separation of d= (0.01pc)M_{8,tot}^{1/3} (T_rest/3.3yr)^{2/3} ~ 0.013 pc (15 light-day). Based on analytical work, the postulated circumbinary disk has an inner radius of 2d = 0.026 pc (30 light-day). SDSS J0159+0105 also displays unusual spectral energy distribution. The unique properties of SDSS J0159+0105 are consistent with it being a centi-parsec SMBHB. The GLSdeDRW code link: http://butler.lab.asu.edu/qso_period/

The Baryon Oscillation Spectroscopic Survey (BOSS) has collected more than 150,000 $2.1 \leq z \leq 3.5$ quasar spectra since 2009. Using this unprecedented sample, we create a composite spectrum in the rest-frame of 102,150 quasar spectra from 800 \AA\ to 3300 \AA\ at a signal-to-noise ratio close to 1000 per pixel ($\Delta v$ of 69 km~s$^{-1}$). Included in this analysis is a correction to account for flux calibration residuals in the BOSS spectrophotometry. We determine the spectral index as a function of redshift of the full sample, warp the composite spectrum to match the median spectral index, and compare the resulting spectrum to SDSS photometry used in target selection. The quasar composite matches the color of the quasar population to within 0.02 magnitudes in $g-r$, 0.03 magnitudes in $r-i$, and 0.01 magnitudes in $i-z$ over the redshift range $2.2<z<2.6$. The composite spectrum deviates from the imaging photometry by 0.05 magnitudes around $z = 2.7$, likely due to differences in target selection as the quasar colors become similar to the stellar locus at this redshift. Finally, we characterize the line features in the high signal-to-noise composite and identify nine faint lines not found in the previous composite spectrum from SDSS.

We describe the first data release from the Spitzer-IRAC Equatorial Survey (SpIES); a large-area survey of 115 deg^2 in the Equatorial SDSS Stripe 82 field using Spitzer during its 'warm' mission phase. SpIES was designed to probe sufficient volume to perform measurements of quasar clustering and the luminosity function at z > 3 to test various models for "feedback" from active galactic nuclei (AGN). Additionally, the wide range of available multi-wavelength, multi-epoch ancillary data enables SpIES to identify both high-redshift (z > 5) quasars as well as obscured quasars missed by optical surveys. SpIES achieves 5{\sigma} depths of 6.13 {\mu}Jy (21.93 AB magnitude) and 5.75 {\mu}Jy (22.0 AB magnitude) at 3.6 and 4.5 microns, respectively - depths significantly fainter than WISE. We show that the SpIES survey recovers a much larger fraction of spectroscopically-confirmed quasars (98%) in Stripe 82 than are recovered by WISE (55%). This depth is especially powerful at high-redshift (z > 3.5), where SpIES recovers 94% of confirmed quasars, whereas WISE only recovers 25%. Here we define the SpIES survey parameters and describe the image processing, source extraction, and catalog production methods used to analyze the SpIES data. In addition to this survey paper, we release 234 images created by the SpIES team and three detection catalogs: a 3.6 {\mu}m-only detection catalog containing 6.1 million sources, a 4.5 {\mu}m-only detection catalog containing 6.5 million sources, and a dual-band detection catalog containing 5.4 million sources.

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.

We present a reverberation mapping (RM) experiment that combines broad- and intermediate-band photometry; it is the first such attempt targeting a sample of 13 quasars at $0.2<z<0.9$. The quasars were selected to have strong H$\alpha$ or H$\beta$ emission lines that are located in one of three intermediate bands (with FWHM around 200 \AA) centered at 8045, 8505, and 9171 \AA. The imaging observations were carried out in the intermediate bands and the broad $i$ and $z$ bands using the prime-focus imager 90Prime on the 2.3m Bok telescope. Because of the large ($\sim$1 deg$^2$) field-of-view (FoV) of 90Prime, we were able to include the 13 quasars within only five telescope pointings or fields. The five fields were repeatedly observed over 20--30 epochs that were unevenly distributed over a duration of 5--6 months. The combination of the broad- and intermediate-band photometry allows us to derive accurate light curves for both optical continuum (from the accretion disk) and line (from the broad-line region, or BLR) emission. We detect H$\alpha$ time lags between the continuum and line emission in 6 quasars. These quasars are at a relatively low redshift range $0.2<z<0.4$. The measured lags are consistent with the current BLR size-luminosity relation for H$\beta$ at $z<0.3$. While this experiment appears successful in detecting lags of the bright H$\alpha$ line, further investigation is required to see if it can also be applied to the fainter H$\beta$ line for quasars at higher redshifts. Finally we demonstrate that by using a small telescope with a large FoV, intermediate-band photometric RM can be efficiently executed for a large sample of quasars at $z>0.2$.

Reverberation mapping (RM) measurements of broad-line region (BLR) lags in z>0.3 quasars are important for directly measuring black hole masses in these distant objects, but so far there have been limited attempts and success given the practical difficulties of RM in this regime. Here we report preliminary results of 15 BLR lag measurements from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project, a dedicated RM program with multi-object spectroscopy designed for RM over a wide redshift range. The lags are based on the 2014 spectroscopic light curves alone (32 epochs over 6 months) and focus on the Hbeta and MgII broad lines in the 100 lowest-redshift (z<0.8) quasars included in SDSS-RM; they represent a small subset of the lags that SDSS-RM (including 849 quasars to z~4.5) is expected to deliver. The reported preliminary lag measurements are for intermediate-luminosity quasars at 0.3<~z<0.8, including 9 Hbeta lags and 6 MgII lags, for the first time extending RM results to this redshift-luminosity regime and providing direct quasar black hole mass estimates over ~ half of cosmic time. The MgII lags also increase the number of known MgII lags by several-fold, and start to explore the utility of MgII for RM at high redshift. The location of these new lags at higher redshifts on the observed BLR size-luminosity relationship is statistically consistent with previous Hbeta results at z<0.3. However, an independent constraint on the relationship slope at z>0.3 is not yet possible due to the limitations in our current sample. Our results demonstrate the general feasibility and potential of multi-object RM for z>0.3 quasars.

We present a study of stellar populations in a sample of spectroscopically-confirmed Lyman-break galaxies (LBGs) and Ly$\alpha$ emitters (LAEs) at $5.7<z<7$. These galaxies have deep optical and infrared images from Subaru, $HST$, and $Spitzer$/IRAC. We focus on a subset of 27 galaxies with IRAC detections, and characterize their stellar populations utilizing galaxy synthesis models based on the multi-band data and secure redshifts. By incorporating nebular emission estimated from the observed Ly$\alpha$ flux, we are able to break the strong degeneracy of model spectra between young galaxies with prominent nebular emission and older galaxies with strong Balmer breaks. The results show that our galaxies cover a wide range of ages from several to a few hundred million years (Myr), and a wide range of stellar masses from $\sim10^8$ to $\sim10^{11}$ $M_{\odot}$. These galaxies can be roughly divided into an `old' subsample and a `young' subsample. The `old' subsample consists of galaxies older than 100 Myr, with stellar masses higher than $10^9$ $M_{\odot}$. The galaxies in the `young' subsample are younger than $\sim$30 Myr, with masses ranging between $\sim10^8$ and $\sim3\times10^9$ $M_{\odot}$. Both subsamples display a correlation between stellar mass and star-formation rate (SFR), but with very different normalizations. The average specific SFR (sSFR) of the `old' subsample is 3--4 Gyr$^{-1}$, consistent with previous studies of `normal' star-forming galaxies at $z\ge6$. The average sSFR of the `young' subsample is an order of magnitude higher, likely due to starburst activity. Our results also indicate little or no dust extinction in the majority of the galaxies, as already suggested by their steep rest-frame UV slopes. Finally, LAEs and LBGs with strong Ly$\alpha$ emission are indistinguishable in terms of age, stellar mass, and SFR.

We explore the variability of quasars in the MgII and Hbeta broad emission lines and UV/optical continuum emission using the Sloan Digital Sky Survey Reverberation Mapping project (SDSS-RM). This is the largest spectroscopic study of quasar variability to date: our study includes 29 spectroscopic epochs from SDSS-RM over $6$ months, containing 357 quasars with MgII and 41 quasars with Hbeta . On longer timescales, the study is also supplemented with two-epoch data from SDSS-I/II. The SDSS-I/II data include an additional $2854$ quasars with MgII and 572 quasars with Hbeta. The MgII emission line is significantly variable ($\Delta f/f$ 10% on 100-day timescales), a necessary prerequisite for its use for reverberation mapping studies. The data also confirm that continuum variability increases with timescale and decreases with luminosity, and the continuum light curves are consistent with a damped random-walk model on rest-frame timescales of $\gtrsim 5$ days. We compare the emission-line and continuum variability to investigate the structure of the broad-line region. Broad-line variability shows a shallower increase with timescale compared to the continuum emission, demonstrating that the broad-line transfer function is not a $\delta$-function. Hbeta is more variable than MgII (roughly by a factor of $1.5$), suggesting different excitation mechanisms, optical depths and/or geometrical configuration for each emission line. The ensemble spectroscopic variability measurements enabled by the SDSS-RM project have important consequences for future studies of reverberation mapping and black hole mass estimation of $1<z<2$ quasars.

We identify 885,503 type 1 quasar candidates to i<22 using the combination of optical and mid-IR photometry. Optical photometry is taken from the Sloan Digital Sky Survey-III: Baryon Oscillation Spectroscopic Survey (SDSS-III/BOSS), while mid-IR photometry comes from a combination of data from the Wide-Field Infrared Survey Explorer (WISE) "ALLWISE" data release and several large-area Spitzer Space Telescope fields. Selection is based on a Bayesian kernel density algorithm with a training sample of 157,701 spectroscopically-confirmed type-1 quasars with both optical and mid-IR data. Of the quasar candidates, 733,713 lack spectroscopic confirmation (and 305,623 are objects that we have not previously classified as photometric quasar candidates). These candidates include 7874 objects targeted as high probability potential quasars with 3.5<z<5 (of which 6779 are new photometric candidates). Our algorithm is more complete to z>3.5 than the traditional mid-IR selection "wedges" and to 2.2<z<3.5 quasars than the SDSS-III/BOSS project. Number counts and luminosity function analysis suggests that the resulting catalog is relatively complete to known quasars and is identifying new high-z quasars at z>3. This catalog paves the way for luminosity-dependent clustering investigations of large numbers of faint, high-redshift quasars and for further machine learning quasar selection using Spitzer and WISE data combined with other large-area optical imaging surveys.

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.

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

According to the recycling scenario, millisecond pulsars (MSPs) have evolved from low-mass X-ray binaries (LMXBs). Their orbits are expected to be circular due to tidal interactions during the binary evolution, as observed in most of the binary MSPs. There are some peculiar systems that do not fit this picture. Three recent examples are PSRs J2234$+$06, J1946$+$3417 and J1950$+$2414, all of which are MSPs in eccentric orbits but with mass functions compatible with expected He white dwarf companions. It has been suggested these MSPs may have formed from delayed accretion-induced collapse of massive white dwarfs, or the eccentricity may be induced by dynamical interaction between the binary and a circumbinary disk. Assuming that the core density of accreting neutron stars in LMXBs may reach the density of quark deconfinement, which can lead to phase transition from neutron stars to strange quark stars, we show that the resultant MSPs are likely to have an eccentric orbit, due to the sudden loss of the gravitational mass of the neutron star during the transition. The eccentricities can be reproduced with a reasonable estimate of the mass loss. This scenario might also account for the formation of the youngest known X-ray binary Cir X$-$1, which also possesses a low-field compact star in an eccentric orbit.

We report the discovery of rapid variations of a high-velocity CIV broad absorption line trough in the quasar SDSS J141007.74+541203.3. This object was intensively observed in 2014 as a part of the Sloan Digital Sky Survey Reverberation Mapping Project, during which 32 epochs of spectroscopy were obtained with the Baryon Oscillation Spectroscopic Survey spectrograph. We observe significant (>4sigma) variability in the equivalent width of the broad (~4000 km/s wide) CIV trough on rest-frame timescales as short as 1.20 days (~29 hours), the shortest broad absorption line variability timescale yet reported. The equivalent width varied by ~10% on these short timescales, and by about a factor of two over the duration of the campaign. We evaluate several potential causes of the variability, concluding that the most likely cause is a rapid response to changes in the incident ionizing continuum. If the outflow is at a radius where the recombination rate is higher than the ionization rate, the timescale of variability places a lower limit on the density of the absorbing gas of n_e > 3.9 x 10^5 cm^-3. The broad absorption line variability characteristics of this quasar are consistent with those observed in previous studies of quasars, indicating that such short-term variability may in fact be common and thus can be used to learn about outflow characteristics and contributions to quasar/host-galaxy feedback scenarios.

We present the discovery of eight quasars at z~6 identified in the Sloan Digital Sky Survey (SDSS) overlap regions. Individual SDSS imaging runs have some overlap with each other, leading to repeat observations over an area spanning >4000 deg^2 (more than 1/4 of the total footprint). These overlap regions provide a unique dataset that allows us to select high-redshift quasars more than 0.5 mag fainter in the z band than those found with the SDSS single-epoch data. Our quasar candidates were first selected as i-band dropout objects in the SDSS imaging database. We then carried out a series of follow-up observations in the optical and near-IR to improve photometry, remove contaminants, and identify quasars. The eight quasars reported here were discovered in a pilot study utilizing the overlap regions at high galactic latitude (|b|>30 deg). These quasars span a redshift range of 5.86<z<6.06 and a flux range of 19.3<z_AB<20.6 mag. Five of them are fainter than z_AB=20 mag, the typical magnitude limit of z~6 quasars used for the SDSS single-epoch images. In addition, we recover eight previously known quasars at z~6 that are located in the overlap regions. These results validate our procedure for selecting quasar candidates from the overlap regions and confirming them with follow-up observations, and provide guidance to a future systematic survey over all SDSS imaging regions with repeat observations.

We present host stellar velocity dispersion measurements for a sample of 88 broad-line quasars at 0.1<z<1 (46 at z>0.6) from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. High signal-to-noise ratio coadded spectra (average S/N~30 per 69 km/s pixel) from SDSS-RM allowed decomposition of the host and quasar spectra, and measurement of the host stellar velocity dispersions and black hole (BH) masses using the single-epoch (SE) virial method. The large sample size and dynamic range in luminosity (L5100=10^(43.2-44.7) erg/s) lead to the first clear detection of a correlation between SE virial BH mass and host stellar velocity dispersion far beyond the local universe. However, the observed correlation is significantly flatter than the local relation, suggesting that there are selection biases in high-z luminosity-threshold quasar samples for such studies. Our uniform sample and analysis enable an investigation of the redshift evolution of the M-sigma relation free of caveats by comparing different samples/analyses at disjoint redshifts. We do not observe evolution of the M-sigma relation in our sample, up to z~1, but there is an indication that the relation flattens towards higher redshifts. Coupled with the increasing threshold luminosity with redshift in our sample, this again suggests certain selection biases are at work, and simple simulations demonstrate that a constant M-sigma relation is favored to z~1. Our results highlight the scientific potential of deep coadded spectroscopy from quasar monitoring programs, and offer a new path to probe the co-evolution of BHs and galaxies at earlier times.

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.

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.

We present distributions of orbital parameters of infalling satellites of $\Lambda$CDM haloes in the mass range $10^{12}-10^{14}$M$_\odot$, which represent the initial conditions for the subsequent evolution of substructures within the host halo. We use merger trees constructed in a high resolution cosmological N-body simulation to trace satellite haloes, and identify the time of infall. We find signficant trends in the distribution of orbital parameters with both the host halo mass and the ratio of satellite-to-host halo masses. For all host halo masses, satellites whose infall mass is a larger fraction of the host halo mass have more eccentric, radially biased orbits. At fixed satellite-to-host halo mass ratio, high mass haloes are biased towards accreting satellites on slightly more radial orbits. To charactise the orbital distributions fully requires fitting the correlated bivariate distribution of two chosen orbital parameters (e.g. radial and tangential velocity or energy and angular momentum). We provide simple fits to one choice of the bivariate distributions, which when transformed faithfully, captures the behaviour of any of the projected one-dimensional distributions.

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.

Numerous theoretical models have long proposed that a strong He II 1640 emission line is the most prominent and unique feature of massive Population III (Pop III) stars in high redshift galaxies. The He II 1640 line strength can constrain the mass and IMF of Pop III stars. We use F132N narrowband filter on the Hubble Space Telescope's (HST) Wide Field Camera 3 (WFC3) to look for strong He II lambda 1640 emission in the galaxy BDF-521 at z=7.01, one of the most distant spectroscopically-confirmed galaxies to date. Using deep F132N narrowband imaging, together with our broadband imaging with F125W and F160W filters, we do not detect He II emission from this galaxy, but place a 2-sigma upper limit on the flux of 5.3x10^-19 ergs s^-1 cm^-2. This measurement corresponds to a 2-sigma upper limit on the Pop III star formation rate (SFR_PopIII) of ~ 0.2 M_solar yr^-1, assuming a Salpeter IMF with 50< M/M_solar < 1000. From the high signal-to-noise broadband measurements in F125W and F160W, we fit the UV continuum for BDF-521. The spectral flux density is ~ 3.6x 10^-11 lambda^-2.32 ergs s^-1 cm^-2 A^-1, which corresponds to an overall unobscured SFR of ~ 5 M_solar yr^-1. Our upper limit on SFR_PopIII suggests that massive Pop III stars represent < 4% of the total star formation. Further, the HST high resolution imaging suggests that BDF-521 is an extremely compact galaxy, with a half-light radius of 0.6 kpc.

We report the serendipitous discoveries of companion galaxies to two high-redshift quasars. SDSS J025617.7+001904 is a z=4.79 quasar included in our recent survey of faint quasars in the SDSS Stripe 82 region. The initial MMT slit spectroscopy shows excess Lyman alpha emission extending well beyond the quasar's light profile. Further imaging and spectroscopy with LBT/MODS1 confirms the presence of a bright galaxy (i_AB = 23.6) located 2arcsec (12 kpc projected) from the quasar with strong Lyman alpha emission (EW_0 ~ 100Ang) at the redshift of the quasar, as well as faint continuum. The second quasar, CFHQS J005006.6+344522 (z=6.25), is included in our recent HST SNAP survey of z~6 quasars searching for evidence of gravitational lensing. Deep imaging with ACS and WFC3 confirms an optical dropout ~4.5 mag fainter than the quasar (Y_AB=25) at a separation of 0.9 arcsec. The red i_775-Y_105 color of the galaxy and its proximity to the quasar (5 kpc projected if at the quasar redshift) strongly favor an association with the quasar. Although it is much fainter than the quasar it is remarkably bright when compared to field galaxies at this redshift, while showing no evidence for lensing. Both systems may represent late-stage mergers of two massive galaxies, with the observed light for one dominated by powerful ongoing star formation and for the other by rapid black hole growth. Observations of close companions are rare; if major mergers are primarily responsible for high-redshift quasar fueling then the phase when progenitor galaxies can be observed as bright companions is relatively short.

The Sloan Digital Sky Survey Reverberation Mapping project (SDSS-RM) is a dedicated multi-object RM experiment that has spectroscopically monitored a sample of 849 broad-line quasars in a single 7 deg$^2$ field with the SDSS-III BOSS spectrograph. The RM quasar sample is flux-limited to i_psf=21.7 mag, and covers a redshift range of 0.1<z<4.5. Optical spectroscopy was performed during 2014 Jan-Jul dark/grey time, with an average cadence of ~4 days, totaling more than 30 epochs. Supporting photometric monitoring in the g and i bands was conducted at multiple facilities including the CFHT and the Steward Observatory Bok telescopes in 2014, with a cadence of ~2 days and covering all lunar phases. The RM field (RA, DEC=14:14:49.00, +53:05:00.0) lies within the CFHT-LS W3 field, and coincides with the Pan-STARRS 1 (PS1) Medium Deep Field MD07, with three prior years of multi-band PS1 light curves. The SDSS-RM 6-month baseline program aims to detect time lags between the quasar continuum and broad line region (BLR) variability on timescales of up to several months (in the observed frame) for ~10% of the sample, and to anchor the time baseline for continued monitoring in the future to detect lags on longer timescales and at higher redshift. SDSS-RM is the first major program to systematically explore the potential of RM for broad-line quasars at z>0.3, and will investigate the prospects of RM with all major broad lines covered in optical spectroscopy. SDSS-RM will provide guidance on future multi-object RM campaigns on larger scales, and is aiming to deliver more than tens of BLR lag detections for a homogeneous sample of quasars. We describe the motivation, design and implementation of this program, and outline the science impact expected from the resulting data for RM and general quasar science.

We present ALMA observations of the [CII] line and far-infrared (FIR) continuum of a normally star-forming galaxy in the reionization epoch, the z=6.96 Ly-alpha emitter (LAE) IOK-1. Probing to sensitivities of sigma_line = 240 micro-Jy/beam (40 km/s channel) and sigma_cont = 21 micro-Jy/beam, we found the galaxy undetected in both [CII] and continuum. Comparison of UV - FIR spectral energy distribution (SED) of IOK-1, including our ALMA limit, with those of several types of local galaxies (including the effects of the cosmic microwave background, CMB, on the FIR continuum) suggests that IOK-1 is similar to local dwarf/irregular galaxies in SED shape rather than highly dusty/obscured galaxies. Moreover, our 3 sigma FIR continuum limit, corrected for CMB effects, implies intrinsic dust mass M_dust < 6.4 x 10^7 M_sun, FIR luminosity L_FIR < 3.7 x 10^{10} L_sun (42.5 - 122.5 micron), total IR luminosity L_IR < 5.7 x 10^{10} L_sun (8 - 1000 micron) and dust-obscured star formation rate (SFR) < 10 M_sun/yr, if we assume that IOK-1 has a dust temperature and emissivity index typical of local dwarf galaxies. This SFR is 2.4 times lower than one estimated from the UV continuum, suggesting that < 29% of the star formation is obscured by dust. Meanwhile, our 3 sigma [CII] flux limit translates into [CII] luminosity, L_[CII] < 3.4 x 10^7 L_sun. Locations of IOK-1 and previously observed LAEs on the L_[CII] vs. SFR and L_[CII]/L_FIR vs. L_FIR diagrams imply that LAEs in the reionization epoch have significantly lower gas and dust enrichment than AGN-powered systems and starbursts at similar/lower redshifts, as well as local star-forming galaxies.

We present and release co-added images of the Sloan Digital Sky Survey (SDSS) Stripe 82. Stripe 82 covers an area of 300 deg^2 on the Celestial Equator, and has been repeatedly scanned 70-90 times in the ugriz bands by the SDSS imaging survey. By making use of all available data in the SDSS archive, our co-added images are optimized for depth. Input single-epoch frames were properly processed and weighted based on seeing, sky transparency, and background noise before co-addition. The resultant products are co-added science images and their associated weight images that record relative weights at individual pixels. The depths of the co-adds, measured as the 5 sigma detection limits of the aperture (3.2 arcsec diameter) magnitudes for point sources, are roughly 23.9, 25.1, 24.6, 24.1, and 22.8 AB magnitudes in the five bands, respectively. They are 1.9-2.2 mag deeper than the best SDSS single-epoch data. The co-added images have good image quality, with an average point-spread function FWHM of ~1 arcsec in the r, i, and z bands. We also release object catalogs that were made with SExtractor. These co-added products have many potential uses for studies of galaxies, quasars, and Galactic structure. We further present and release near-IR J-band images that cover ~90 deg^2 of Stripe 82. These images were obtained using the NEWFIRM camera on the NOAO 4-m Mayall telescope, and have a depth of about 20.0--20.5 Vega magnitudes (also 5 sigma detection limits for point sources).

High-redshift quasars are currently the only probes of the growth of supermassive black holes and potential tracers of structure evolution at early cosmic time. Here we present our candidate selection criteria from the Panoramic Survey Telescope & Rapid Response System 1 and follow-up strategy to discover quasars in the redshift range 5.7<z<6.2. With this strategy we discovered eight new 5.7<z<6.0 quasars, increasing the number of known quasars at z>5.7 by more than 10%. We additionally recovered 18 previously known quasars. The eight quasars presented here span a large range of luminosities (-27.3 < M_{1450} < -25.4; 19.6 < z_ps1 < 21.2) and are remarkably heterogeneous in their spectral features: half of them show bright emission lines whereas the other half show a weak or no Ly$\alpha$ emission line (25% with rest-frame equivalent width of the Ly$\alpha$ + Nv line lower than 15{\AA}). We find a larger fraction of weak-line emission quasars than in lower redshift studies. This may imply that the weak-line quasar population at the highest redshifts could be more abundant than previously thought. However, larger samples of quasars are needed to increase the statistical significance of this finding.

We present a new algorithm which groups the subhaloes found in cosmological N- body simulations by structure finders such as SUBFIND into dark matter haloes whose formation histories are strictly hierarchical. One advantage of these `Dhaloes' over the commonly used friends-of-friends (FoF) haloes is that they retain their individual identity in cases when FoF haloes are artificially merged by tenuous bridges of particles or by an overlap of their outer diffuse haloes. Dhaloes are thus well-suited for modelling galaxy formation and their merger trees form the basis of the Durham semi-analytic galaxy formation model, GALFORM. Applying the Dhalo construction to the LambdaCDM Millennium-2 simulation we find that approximately 90% of Dhaloes have a one-to-one, bijective match with a corresponding FoF halo. The remaining 10% are typically secondary components of large FoF haloes. Although the mass functions of both types of haloes are similar, the mass of Dhaloes correlates much more tightly with the virial mass, M200, than FoF haloes. Approximately 80% of FoF and bijective and non-bijective Dhaloes are relaxed according to standard criteria. For these relaxed haloes all three types have similar concentration-M200 relations and, at fixed mass, the concentration distributions are described accurately by log-normal distributions.

In a superhigh magnetic field, direct Urca reactions can proceed for an arbitrary proton concentration. Since only the electrons with high energy $E$ ($E > Q$, $Q$ is the threshold energy of inverse $\beta-$decay) at large Landau levels can be captured, we introduce the Landau level effect coefficient $q$ and the effective electron capture rate $\Gamma_{\rm eff}$. By using $\Gamma_{\rm eff}$, the values of $L_{\rm X}$ and $L_{\rm \nu}$ are calculated, where and $L_{\rm \nu}$, $L_{\rm X}$ are the average neutrino luminosity of Anomalous X-ray Pulsars (AXPs) and the average X-ray luminosity of AXPs, respectively. The complete process of electron capture inside a magnetar is simulated numerically.

In this paper, we consider the effect of Landau levels on the decay of superhigh magnetic fields of magnetars. Applying ${}^3P_2$ anisotropic neutron superfluid theory yield a second-order differential equation for a superhigh magnetic field $B$ and its evolutionary timescale $t$. The superhigh magnetic fields may evolve on timescales $\sim (10^{6}-10^{7})$ yrs for common magnetars. According to our model, the activity of a magnetar may originate from instability caused by the high electron Fermi energy.

We present the Data Release 10 Quasar (DR10Q) catalog from the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey III. The catalog includes all BOSS objects that were targeted as quasar candidates during the first 2.5 years of the survey and that are confirmed as quasars via visual inspection of the spectra. The catalog also includes known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. The catalog contains 166,583 quasars (74,454 are new discoveries since SDSS-DR9) detected over 6,373 deg$^{2}$ with robust identification and redshift measured by a combination of principal component eigenspectra. The number of quasars with $z>2.15$ (117,668) is $\sim$5 times greater than the number of $z>2.15$ quasars known prior to BOSS. Redshifts and FWHMs are provided for the strongest emission lines (CIV, CIII, MgII). The catalog identifies 16,461 broad absorption line quasars and gives their characteristics. For each object, the catalog presents five-band (u, g, r, i, z) CCD-based photometry with typical accuracy of 0.03 mag and information on the optical morphology and selection method. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra cover the wavelength region 3,600-10,500\AA\ at a spectral resolution in the range 1,300$<$R$<$2,500; the spectra can be retrieved from the SDSS Catalog Archive Server. We also provide a supplemental list of an additional 2,376 quasars that have been identified among the galaxy targets of the SDSS-III/BOSS.

IGR J17480$-$2446 is an accreting X-ray pulsar in a low-mass X-ray binary harbored in the Galactic globular cluster Terzan 5. Compared with other accreting millisecond pulsars, IGR J17480$-$2446 is peculiar in its low spin frequency (11 Hz), which suggests that it might be a mildly recycled neutron star at the very early phase of mass transfer. However, this model seems to be in contrast with the low field strength deduced from the kiloHertz quasi-periodic oscillations observed in IGR J17480$-$2446. Here we suggest an alternative interpretation, assuming that the current binary system was formed during an exchange encounter either between a binary (which contains a recycled neutron star) and the current donor, or between a binary and an isolated, recycled neutron star. In the resulting binary, the spin axis of the neutron star could be parallel or anti-parallel with the orbital axis. In the later case, the abnormally low frequency of IGR J17480$-$2446 may result from the spin-down to spin-up evolution of the neutron star. We also briefly discuss the possible observational implications of the pulsar in this scenario.

We present a deep LBT/LBC U-band imaging survey (9 deg2) covering the NOAO Bootes field. A total of 14,485 Lyman Break Galaxies (LBGs) at z~3 are selected, which are used to measure the rest-frame UV luminosity function (LF). The large sample size and survey area reduce the LF uncertainties due to Poisson statistics and cosmic variance by >3 compared to previous studies. At the bright end, the LF shows excess power compared to the best-fit Schechter function, which can be attributed to the contribution of $z\sim3$ quasars. We compute the rest-frame near-infrared LF and stellar mass function (SMF) of z~3 LBGs based on the R-band and IRAC [4.5 micro m]-band flux relation. We investigate the evolution of the UV LFs and SMFs between z~7 and z~3, which supports a rising star formation history in the LBGs. We study the spatial correlation function of two bright LBG samples and estimate their average host halo mass. We find a tight relation between the host halo mass and the galaxy star formation rate (SFR),which follows the trend predicted by the baryonic accretion rate onto the halo, suggesting that the star formation in LBGs is fueled by baryonic accretion through the cosmic web. By comparing the SFRs with the total baryonic accretion rates, we find that cosmic star formation efficiency is about 5%-20% and it does not evolve significantly with redshift, halo mass, or galaxy luminosity.

We present ALMA observations of the [C II] 158 micron fine structure line and dust continuum emission from the host galaxies of five redshift 6 quasars. We also report complementary observations of 250 GHz dust continuum and CO (6-5) line emission from the z=6.00 quasar SDSS J231038.88+185519.7. The ALMA observations were carried out in the extended array at 0.7" resolution. We have detected the line and dust continuum in all five objects. The derived [C II] line luminosities are 1.6x10^{9} to 8.8x10^{9} Lsun and the [C II]-to-FIR luminosity ratios are 3.0-5.6x10^{-4}, which is comparable to the values found in other high-redshift quasar-starburst systems and local ultra-luminous infrared galaxies. The sources are marginally resolved and the intrinsic source sizes (major axis FWHM) are constrained to be 0.3" to 0.6" (i.e., 1.7 to 3.5 kpc) for the [C II] line emission and 0.2" to 0.4" (i.e., 1.2 to 2.3 kpc) for the continuum. These measurements indicate that there is vigorous star formation over the central few kpc in the quasar host galaxies. The ALMA observations also constrain the dynamical properties of the atomic gas in the starburst nuclei. The intensity-weighted velocity maps of three sources show clear velocity gradients. Such velocity gradients are consistent with a rotating, gravitationally bound gas component, although they are not uniquely interpreted as such. Under the simplifying assumption of rotation, the implied dynamical masses within the [C II]-emitting regions are of order 10^{10} to 10^{11} Msun. Given these estimates, the mass ratios between the SMBHs and the spheroidal bulge are an order of magnitude higher than the mean value found in local spheroidal galaxies, which is in agreement with results from previous CO observations of high redshift quasars.

We present a detailed structural and morphological study of a large sample of spectroscopically-confirmed galaxies at z >= 6, using deep HST near-IR broad-band images and Subaru optical narrow-band images. The galaxy sample consists of 51 Lyman-alpha emitters (LAEs) at z ~ 5.7, 6.5, and 7.0, and 16 Lyman-break galaxies (LBGs) at 5.9 < z < 6.5. These galaxies exhibit a wide range of rest-frame UV continuum morphology in the HST images, from compact features to multiple component systems. The fraction of merging/interacting galaxies reaches 40% ~ 50% at the brightest end of M_1500 <= -20.5 mag. The intrinsic half-light radii r_{hl,in}, after correction for PSF broadening, are roughly between r_{hl,in} ~ 0.05" (0.3 kpc) and 0.3" (1.7 kpc) at M_1500 <= -19.5 mag. The median r_{hl,in} value is 0.16" (~0.9 kpc). This is consistent with the sizes of bright LAEs and LBGs at z >= 6 in previous studies. In addition, more luminous galaxies tend to have larger sizes, exhibiting a weak size-luminosity relation r_{hl,in} \propto L^{0.14} at M_1500 <= -19.5 mag. The slope of 0.14 is significantly flatter than those in fainter LBG samples. We discuss the morphology of z >= 6 galaxies with nonparametric methods, including the CAS system and the Gini and M_20 parameters, and demonstrate their validity through simulations. We search for extended Lyman-alpha emission halos around LAEs at z ~ 5.7 and 6.5, by stacking a number of narrow-band images. We do not find evidence of extended halos predicted by cosmological simulations. Such Lyman-alpha halos, if they exist, could be weaker than predicted. Finally, we investigate any positional misalignment between UV continuum and Lyman-alpha emission in LAEs. While the two positions are generally consistent, several merging galaxies show significant positional differences. This is likely caused by a disturbed ISM distribution due to merging activity.

We report on a clear solar-cycle variation of the Sun's shadow in the 10 TeV cosmic-ray flux observed by the Tibet air shower array during a full solar cycle from 1996 to 2009. In order to clarify the physical implications of the observed solar cycle variation, we develop numerical simulations of the Sun's shadow, using the Potential Field Source Surface (PFSS) model and the Current Sheet Source Surface (CSSS) model for the coronal magnetic field. We find that the intensity deficit in the simulated Sun's shadow is very sensitive to the coronal magnetic field structure, and the observed variation of the Sun's shadow is better reproduced by the CSSS model. This is the first successful attempt to evaluate the coronal magnetic field models by using the Sun's shadow observed in the TeV cosmic-ray flux.

We present deep HST near-IR and Spitzer mid-IR observations of a large sample of spectroscopically-confirmed galaxies at z >= 6. The sample consists of 51 Lyman-alpha emitters (LAEs) at z ~ 5.7, 6.5, and 7.0, and 16 Lyman-break galaxies (LBGs) at 5.9 < z < 6.5. The near-IR images were mostly obtained with WFC3 in the F125W and F160W bands, and the mid-IR images were obtained with IRAC in the 3.6um and 4.5um bands. Our galaxies also have deep optical imaging data from Subaru Suprime-Cam. We utilize the multi-band data and secure redshifts to derive their rest-frame UV properties. These galaxies have steep UV continuum slopes roughly between beta ~ -1.5 and -3.5, with an average value of beta ~ -2.3, slightly steeper than the slopes of LBGs in previous studies. The slope shows little dependence on UV continuum luminosity except for a few of the brightest galaxies. We find a statistically significant excess of galaxies with slopes around beta ~ -3, suggesting the existence of very young stellar populations with extremely low metallicity and dust content. Our galaxies have moderately strong rest-frame Lyman-alpha equivalent width (EW) in a range of ~10 to ~200 \AA. The star-formation rates are also moderate, from a few to a few tens solar masses per year. The LAEs and LBGs in this sample share many common properties, implying that LAEs represent a subset of LBGs with strong Lyman-alpha emission. Finally, the comparison of the UV luminosity functions between LAEs and LBGs suggests that there exists a substantial population of faint galaxies with weak Lyman-alpha emission (EW < 20 \AA) that could be the dominant contribution to the total ionizing flux at z >= 6.

We present deep spectroscopic observations of a Lyman-alpha emitter (LAE) candidate at z ~ 7.7 using the infrared spectrograph LUCI on the 2 x 8.4m Large Binocular Telescope (LBT). The candidate is the brightest among the four z ~ 7.7 LAE candidates found in a narrow-band imaging survey by Krug et al. 2012. Our spectroscopic data include a total of 7.5 hours of integration with LBT/LUCI and are deep enough to significantly (3.2-4.9 sigma) detect the Lyman-alpha emission line of this candidate, based on its Lyman-alpha flux 1.2 x 10^{-17} erg s^{-1} cm^{-2} estimated from the narrow-band photometry. However, we do not find any convincing signal at the expected position of its Lyman-alpha emission line, suggesting that this source is not an LAE at z ~ 7.7. The non-detection in this work, together with the previous studies of z ~ 7.7 LAEs, puts a strong constraint on the bright-end Lyman-alpha luminosity function (LF) at z ~ 7.7. We find a rapid evolution of the Lyman-alpha LF from z ~ 6.5 to 7.7: the upper limit of the z ~ 7.7 LF is more than 5 times lower than the z ~ 6.5 LF at the bright end (f > 1.0 x 10^{-17} erg s^{-1} cm^{-2}, or L > 6.9 x 10^{42} erg s^{-1}). This is likely caused by an increasing neutral fraction in the IGM that substantially attenuates Lyman-alpha emission at z ~ 7.7.

We present a measurement of the Type I quasar luminosity function at z=5 using a large sample of spectroscopically confirmed quasars selected from optical imaging data. We measure the bright end (M_1450<-26) with Sloan Digital Sky Survey (SDSS) data covering ~6000 deg^2, then extend to lower luminosities (M_1450<-24) with newly discovered, faint z~5 quasars selected from 235 deg^2 of deep, coadded imaging in the SDSS Stripe 82 region (the celestial equator in the Southern Galactic Cap). The faint sample includes 14 quasars with spectra obtained as ancillary science targets in the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS), and 59 quasars observed at the MMT and Magellan telescopes. We construct a well-defined sample of 4.7<z<5.1 quasars that is highly complete, with 73 spectroscopic identifications out of 92 candidates. Our color selection method is also highly efficient: of the 73 spectra obtained, 71 are high redshift quasars. These observations reach below the break in the luminosity function (M_1450* ~ -27). The bright end slope is steep (beta <~ -4), with a constraint of beta < -3.1 at 95% confidence. The break luminosity appears to evolve strongly at high redshift, providing an explanation for the flattening of the bright end slope reported previously. We find a factor of ~2 greater decrease in the number density of luminous quasars (M_1450<-26) from z=5 to z=6 than from z=4 to z=5, suggesting a more rapid decline in quasar activity at high redshift than found in previous surveys. Our model for the quasar luminosity function predicts that quasars generate ~30% of the ionizing photons required to keep the universe ionized at z=5.

A new EAS hybrid experiment has been designed by constructing a YAC (Yangbajing Air shower Core) detector array inside the existing Tibet-III air shower array. The first step of YAC, called "YAC-I", consists of 16 plastic scintillator units (4 rows times 4 columns) each with an area of 40 cm * 50 cm which is used to check hadronic interaction models used in AS simulations. A Monte Carlo study shows that YAC-I can record high energy electromagnetic component in the core region of air showers induced by primary particles of several tens TeV energies where the primary composition is directly measured by space experiments. It may provide a direct check of the hadronic interaction models currently used in the air shower simulations in the corresponding energy region. In present paper, the method of the observation and the sensitivity of the characteristics of the observed events to the different interaction models are discussed.

A new hybrid experiment has been started by AS{\gamma} experiment at Tibet, China, since August 2011, which consists of a low threshold burst-detector-grid (YAC-II, Yangbajing Air shower Core array), the Tibet air-shower array (Tibet-III) and a large underground water Cherenkov muon detector (MD). In this paper, the capability of the measurement of the chemical components (proton, helium and iron) with use of the (Tibet-III+YAC-II) is investigated by means of an extensive Monte Carlo simulation in which the secondary particles are propagated through the (Tibet-III+YAC-II) array and an artificial neural network (ANN) method is applied for the primary mass separation. Our simulation shows that the new installation is powerful to study the chemical compositions, in particular, to obtain the primary energy spectrum of the major component at the knee.

We present the Data Release 9 Quasar (DR9Q) catalog from the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey III. The catalog includes all BOSS objects that were targeted as quasar candidates during the survey, are spectrocopically confirmed as quasars via visual inspection, have luminosities Mi[z=2]<-20.5 (in a $\Lambda$CDM cosmology with H0 = 70 km/s/Mpc, $\Omega_{\rm M}$ = 0.3, and $\Omega_{\Lambda}$ = 0.7) and either display at least one emission line with full width at half maximum (FWHM) larger than 500 km/s or, if not, have interesting/complex absorption features. It includes as well, known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. This catalog contains 87,822 quasars (78,086 are new discoveries) detected over 3,275 deg$^{2}$ with robust identification and redshift measured by a combination of principal component eigenspectra newly derived from a training set of 8,632 spectra from SDSS-DR7. The number of quasars with $z>2.15$ (61,931) is ~2.8 times larger than the number of z>2.15 quasars previously known. Redshifts and FWHMs are provided for the strongest emission lines (CIV, CIII], MgII). The catalog identifies 7,533 broad absorption line quasars and gives their characteristics. For each object the catalog presents five-band (u,g,r,i,z) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys.

We have identified a very interesting Ly-alpha emitter, whose Ly-alpha emission line has an extremely large observed equivalent width of EW_0=436^{+422}_{-149}A, which corresponds to an extraordinarily large intrinsic rest-frame equivalent width of EW_0^{int}=872^{+844}_{-298}A after the average intergalactic absorption correction. The object was spectroscopically confirmed to be a real Ly-alpha emitter by its apparent asymmetric Ly-alpha line profile detected at z=6.538. The continuum emission of the object was definitely detected in our deep z'-band image; thus, its EW_0 was reliably determined. Follow-up deep near-infrared spectroscopy revealed emission lines of neither He II lambda1640 as an apparent signature of Population III, nor C IV lambda1549 as a proof of active nucleus. No detection of short-lived He II lambda1640 line is not necessarily inconsistent with the interpretation that the underlying stellar population of the object is dominated by Population III. We found that the observed extremely large EW_0 of the Ly-alpha emission and the upper limit on the EW_0 of the He II lambda1640 emission can be explained by population synthesis models favoring a very young age less than 2-4Myr and massive metal-poor (Z<10^{-5}) or even metal-free stars. The observed large EW_0 of Ly-alpha is hardly explained by Population I/II synthesis models with Z>10^{-3}. However, we cannot conclusively rule out the possibility that this object is composed of a normal stellar population with a clumpy dust distribution, which could enhance the Ly-alpha EW_0, though its significance is still unclear.

We present one of the most ultraviolet (UV) luminous Lyman Break Galaxies (LBGs) (J1432+3358) at z=2.78, discovered in the NOAO Deep Wide-Field Survey (NDWFS) Bootes field. The R-band magnitude of J1432+3358 is 22.29 AB, more than two magnitudes brighter than typical L* LBGs at this redshift. The deep z-band image reveals two components of J1432+3358 separated by 1.0" with flux ratio of 3:1. The high signal-to-noise ratio (S/N) rest-frame UV spectrum shows Lya emission line and interstellar medium absorption lines. The absence of NV and CIV emission lines, the non-detection in X-ray and radio wavelengths and mid-infrared (MIR) colors indicate no or weak active galactic nuclei (AGN) (<10%) in this galaxy. The galaxy shows broader line profile with the full width half maximum (FWHM) of about 1000 km/s and larger outflow velocity (~500 km/s) than those of typical z~3 LBGs. The physical properties are derived by fitting the spectral energy distribution (SED) with stellar synthesis models. The dust extinction, E(B-V)=0.12, is similar to that in normal LBGs. The star formation rates (SFRs) derived from the SED fitting and the dust-corrected UV flux are consistent with each other, ~300 Msun/yr, and the stellar mass is 1.3e11 Msun. The SFR and stellar mass in J1432+3358 are about an order of magnitude higher than those in normal LBGs. The SED-fitting results support that J1432+3358 has a continuous star formation history with the star formation episode of 630 Myr. The morphology of J1432+3358 and its physical properties suggest that J1432+3358 is in an early phase of 3:1 merger process. The unique properties and the low space number density (~1e-7 Mpc^{-3})are consistent with the interpretation that such galaxies are either found in a short unobscured phase of the star formation or that small fraction of intensive star-forming galaxies are unobscured.

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

We report the discovery of a protocluster at z~6 containing at least eight cluster member galaxies with spectroscopic confirmations in the wide-field image of the Subaru Deep Field (SDF). The overdensity of the protocluster is significant at the 6 sigma level, based on the surface number density of i'-dropout galaxies. The overdense region covers ~36 sq. arcmin, and includes 30 i'-dropout galaxies. Follow-up spectroscopy revealed that 15 of these are real z~6 galaxies (5.7 < z < 6.3). Eight of the 15 are clustering in a narrow redshift range centered at z=6.01, corresponding to a seven-fold increase in number density over the average in redshift space. We found no significant difference in the observed properties, such as Ly-alpha luminosities and UV continuum magnitudes, between the eight protocluster members and the seven non-members. The velocity dispersion of the eight protocluster members is 647 km/s, which is about three times higher than that predicted by the standard cold dark matter model. This discrepancy could be attributed to the distinguishing three-dimensional distribution of the eight protocluster members. We discuss two possible explanations for this discrepancy: either the protocluster is already mature, with old galaxies at the center, or it is still immature and composed of three subgroups merging to become a larger cluster. In either case, this concentration of z=6.01 galaxies in the SDF may be one of the first sites of formation of a galaxy cluster in the universe.

We present details of the construction and characterization of the coaddition of the Sloan Digital Sky Survey Stripe 82 \ugriz\ imaging data. This survey consists of 275 deg$^2$ of repeated scanning by the SDSS camera of $2.5\arcdeg$ of $\delta$ over $-50\arcdeg \le \alpha \le 60\arcdeg$ centered on the Celestial Equator. Each piece of sky has $\sim 20$ runs contributing and thus reaches $\sim2$ magnitudes fainter than the SDSS single pass data, i.e. to $r\sim 23.5$ for galaxies. We discuss the image processing of the coaddition, the modeling of the PSF, the calibration, and the production of standard SDSS catalogs. The data have $r$-band median seeing of 1.1\arcsec, and are calibrated to $\le 1%$. Star color-color, number counts, and psf size vs modelled size plots show the modelling of the PSF is good enough for precision 5-band photometry. Structure in the psf-model vs magnitude plot show minor psf mis-modelling that leads to a region where stars are being mis-classified as galaxies, and this is verified using VVDS spectroscopy. As this is a wide area deep survey there are a variety of uses for the data, including galactic structure, photometric redshift computation, cluster finding and cross wavelength measurements, weak lensing cluster mass calibrations, and cosmic shear measurements.

Supernova (SN) 2009ig was discovered 17 hours after explosion by the Lick Observatory Supernova Search, promptly classified as a normal Type Ia SN (SN Ia), peaked at V = 13.5 mag, and was equatorial, making it one of the foremost supernovae for intensive study in the last decade. Here, we present ultraviolet (UV) and optical observations of SN 2009ig, starting about 1 day after explosion until around maximum brightness. Our data include excellent UV and optical light curves, 25 premaximum optical spectra, and 8 UV spectra, including the earliest UV spectrum ever obtained of a SN Ia. SN 2009ig is a relatively normal SN Ia, but does display high-velocity ejecta - the ejecta velocity measured in our earliest spectra (v ~ -23,000 km/s for Si II 6355) is the highest yet measured in a SN Ia. The spectral evolution is very dramatic at times earlier than 12 days before maximum brightness, but slows after that time. The early-time data provide a precise measurement of 17.13 +/- 0.07 days for the SN rise time. The optical color curves and early-time spectra are significantly different from template light curves and spectra used for light-curve fitting and K-corrections, indicating that the template light curves and spectra do not properly represent all Type Ia supernovae at very early times. In the age of wide-angle sky surveys, SNe like SN 2009ig that are nearby, bright, well positioned, and promptly discovered will still be rare. As shown with SN 2009ig, detailed studies of single events can provide significantly more information for testing systematic uncertainties related to SN Ia distance estimates and constraining progenitor and explosion models than large samples of more distant SNe.

At high redshift, the universe was so young that core-collapse supernovae (SNe) are suspected to be the dominant source of dust production. However, some observations indicate that the dust production by SNe is an inefficient process, casting doubts on the existence of abundant SNe-dust in the early universe. Recently, Perley et al. (2010) reported that the afterglow of GRB 071025 - an unusually red GRB at z ~ 5 - shows evidence for the SNe-produced dust. Since this is perhaps the only high redshift GRB exhibiting compelling evidence for SNe-dust but the result could easily be affected by small systematics in photometry, we re-examined the extinction properties of GRB 071025 using our own optical/near-infrared data at a different epoch. In addition, we tested SNe-dust models with different progenitor masses and dust destruction efficiencies to constrain the dust formation mechanisms. By searching for the best-fit model of the afterglow spectral energy distribution, we confirm the previous claim that the dust in GRB 071025 is most likely to originate from SNe. We also find that the SNe-dust model of 13 or 25M_{\odot} without dust destruction fits the extinction property of GRB 071025 best, while pair-instability SNe (PISNe) models with a 170M_{\odot} progenitor poorly fit the data. Our results indicate that, at least in some systems at high redshift, SNe with intermediate masses within 10 - 30M_{\odot} were the main contributors for the dust enrichment, and the dust destruction effect due to reverse shock was negligible.

We present Keck spectroscopic observations of z>6 Lyman-break galaxy (LBG) candidates in the Subaru Deep Field (SDF). The candidates were selected as i'-dropout objects down to z'=27 AB magnitudes from an ultra-deep SDF z'-band image. With the Keck spectroscopy we identified 19 LBGs with prominent Ly_alpha emission lines at 6< z < 6.4. The median value of the Ly_alpha rest-frame equivalent widths (EWs) is ~50 A, with four EWs >100 A. This well-defined spectroscopic sample spans a UV-continuum luminosity range of -21.8< M_{UV}<-19.5 (0.6~5 L*_{UV}) and a Ly_alpha luminosity range of (0.3~3) x 10^{43} erg s^{-1} (0.3~3 L*_ {Ly_alpha}). We derive the UV and Ly_alpha luminosity functions (LFs) from our sample at <z>~6.2 after we correct for sample incompleteness. We find that our measurement of the UV LF is consistent with the results of previous studies based on photometric LBG samples at 5<z<7. Our Ly_alpha LF is also generally in agreement with the results of Ly_alpha-emitter surveys at z~5.7 and 6.6. This study shows that deep spectroscopic observations of LBGs can provide unique constraints on both the UV and Ly_alpha LFs at z>6.