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

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

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

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 study the physical mechanisms that cause the offset between low-redshift and high-redshift galaxies on the [OIII]/H$\beta$ versus [NII]/H$\alpha$ ``Baldwin, Phillips & Terlevich'' (BPT) diagram using a sample of local analogues of high-redshift galaxies. These high-redshift analogue galaxies are selected from the Sloan Digital Sky Survey. Located in the same region on the BPT diagram as the ultra-violet selected galaxies at $z\sim2$, these high-redshift analogue galaxies provide an ideal local benchmark to study the offset between the local and high-redshift galaxies on the BPT diagram. We compare the nitrogen-to-oxygen ratio (N/O), the shape of the ionising radiation field, and ionisation parameters between the high-redshift analogues and a sample of local reference galaxies. The higher ionisation parameter in the high-redshift analogues is the dominant physical mechanism driving the BPT offset from low- to high-redshift, particularly at high {\nii/\ha}. Furthermore, the N/O ratio enhancement also plays a minor role to cause the BPT offset. However, the shape of the ionising radiation field is unlikely to cause the BPT offset because the high-redshift analogues have a similar hard ionising radiation field as local reference galaxies. This hard radiation field cannot be produced by the current standard stellar synthesis models. The stellar rotation and binarity may help solve the discrepancy.

We measure the LyC escape fraction in 54 faint Lyman Alpha Emitters (LAEs) at $z\simeq3.1$ in the GOODS-South field. With the average magnitude of $R=26.7$ AB ($M_{UV}=-18.8$, $L\simeq0.1L^*$), these galaxies represent a population of compact young dwarf galaxies. Their properties are likely to resemble those in the galaxies responsible for reionising the Universe at $z>6$. We do not detect LyC emission in any individual LAEs in the deep {\em HST} F336W images, which covers the rest-frame 820\r{A}. We do not detect the LyC emission of these LAEs in the stacked F336W images, either. The $3\sigma$ upper limit of LyC escape fractions is $f_{\rm esc}<14-32\%$. However, the high Ly$\alpha$ rest-frame equivalent width, low stellar mass and UV luminosity of these LAEs suggest that they should have $f_{\rm esc}>50\%$. The low LyC escape fraction from this work and other stacking analysis suggest that the LyC leaking galaxies with $f_{\rm esc} > 50\%$ at $z=2-3$ do not follow the relation between the $f_{\rm esc}$ and UV luminosity and Ly$\alpha$ equivalent width (EW) derived from typical galaxies at similar redshift. Therefore, the UV luminosity and Ly$\alpha$ equivalent width (EW) are not the best indicators for the LyC escape fraction.

We present results from a low-resolution spectroscopic survey for 21 galaxy clusters at $0.4 < z < 0.8$ selected from the ESO Distant Cluster Survey. We measured spectra using the low-dispersion prism in IMACS on the Magellan Baade telescope and calculate redshifts with an accuracy of $\sigma_z = 0.007$. We find 1763 galaxies that are brighter than $R = 22.9$ in the large-scale cluster environs. We identify the galaxies expected to be accreted by the clusters as they evolve to $z = 0$ using spherical infall models and find that $\sim30\%$ to $\sim70\%$ of the $z = 0$ cluster population lies outside the virial radius at $z \sim 0.6$. For analogous clusters at $z = 0$, we calculate that the ratio of galaxies that have fallen into the clusters since $z \sim 0.6$ to those that were already in the core at that redshift is typically between $\sim0.3$ and $1.5$. This wide range of ratios is due to intrinsic scatter and is not a function of velocity dispersion, so a variety of infall histories is to be expected for clusters with current velocity dispersions of $300 \lesssim\sigma\lesssim 1200$ km s$^{-1}$. Within the infall regions of $z \sim 0.6$ clusters, we find a larger red fraction of galaxies than in the field and greater clustering among red galaxies than blue. We interpret these findings as evidence of "preprocessing", where galaxies in denser local environments have their star formation rates affected prior to their aggregation into massive clusters, although the possibility of backsplash galaxies complicates the interpretation.

The most luminous quasars at high redshift harbour the fastest-growing and most massive black holes in the early Universe. They are exceedingly rare and hard to find. Here, we present our search for the most luminous quasars in the redshift range from $z=4.5$ to $5$ using data from SkyMapper, Gaia and WISE. We use colours to select likely high-redshift quasars and reduce the stellar contamination of the candidate set with parallax and proper motion data. In $\sim$12,500~deg$^2$ of Southern sky, we find 92 candidates brighter than $R_p=18.2$. Spectroscopic follow-up has revealed 21 quasars at $z\ge 4$ (16 of which are within $z=[4.5,5]$), as well as several red quasars, BAL quasars and objects with unusual spectra, which we tentatively label OFeLoBALQSOs at redshifts of $z\approx 1$ to $2$. This work lifts the number of known bright $z\ge 4.5$ quasars in the Southern hemisphere from 10 to 26 and brings the total number of quasars known at $R_p<18.2$ and $z\ge 4.5$ to 42.

We examine the diagnostic power of rest-frame ultraviolet (UV) nebular emission lines, and compare them to more commonly used rest-frame optical emission lines, using the test case of a single star-forming knot of the bright lensed galaxy RCSGA 032727-132609 at redshift z~1.7. This galaxy has complete coverage of all the major rest-frame UV and optical emission lines from Magellan/MagE and Keck/NIRSPEC. Using the full suite of diagnostic lines, we infer the physical properties: nebular electron temperature (T_e), electron density (n_e), oxygen abundance (log(O/H)), ionisation parameter (log(q)) and interstellar medium (ISM) pressure (log(P/k)). We examine the effectiveness of the different UV, optical and joint UV-optical spectra in constraining the physical conditions. Using UV lines alone we can reliably estimate log(q), but the same is difficult for log(O/H). UV lines yield a higher (~1.5 dex) log(P/k) than the optical lines, as the former probes a further inner nebular region than the latter. For this comparison, we extend the existing Bayesian inference code IZI, adding to it the capability to infer ISM pressure simultaneously with metallicity and ionisation parameter. This work anticipates future rest-frame UV spectral datasets from the James Webb Space Telescope (JWST) at high redshift and from the Extremely Large Telescope (ELT) at moderate redshift.

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

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.

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

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

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

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 study the relation between the metallicity of quasar broad line region (BLR) and black hole (BH) mass ($10^{7.5}M_{\odot} \sim 10^{10}M_{\odot}$) and quasar bolometric luminosity ($10^{44.6}erg/s \sim 10^{48} erg/s$) using a sample of $\sim$130,000 quasars at $2.25\leq z\leq 5.25$ from Sloan Digital Sky Survey (SDSS) Data Release 12 (DR12). We generate composite spectra by stacking individual spectra in the same BH mass (bolometric luminosity) and redshift bins and then estimate the metallicity of quasar BLR using metallicity-sensitive broad emission-line flux ratios based on the photoionization models. We find a significant correlation between quasar BLR metallicity and BH mass (bolometric luminosity) but no correlation between quasar BLR metallicity and redshift. We also compare the metallicity of quasar BLR and that of host galaxies inferred from the mass-metallicity relation of star-forming galaxy at $z\sim2.3$ and $z\sim3.5$. We find quasar BLR metallicity is 0.3 $\sim$ 1.0 dex higher than their host galaxies. This discrepancy cannot be interpreted by the uncertainty due to different metallicity diagnostic methods, mass-metallicity relation of galaxy, metallicity gradient in quasar host galaxies, BH mass estimation, the effect of different spectral energy distribution (SED) models, and a few other potential sources of uncertainties. We propose a possibility that the high metallicity in quasar BLR might be caused by metal enrichment from massive star formation in the nucleus region of quasars or even the accretion disk.

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 a comparative study of the relation between the aperture-based gas-phase metallicity and three structural parameters of star-forming galaxies: mass ($\mathrm{M \equiv M_*}$), average potential ($\Phi \equiv \mathrm{M_*/R_e}$) and average surface mass density ($\Sigma \equiv \mathrm{M_*/R_e^2}$; where $\mathrm{R_e}$ is the effective radius). We use a volume-limited sample drawn from the publicly available SDSS DR7, and base our analysis on aperture-matched sampling by selecting sets of galaxies where the SDSS fibre probes a fixed fraction of $\mathrm{R_e}$. We find that between 0.5 and 1.5 $\mathrm{R_e}$, the gas-phase metallicity correlates more tightly with $\Phi$ than with either $\mathrm{M}$ or $\Sigma$, in that for all aperture-matched samples, the potential-metallicity relation has (i) less scatter, (ii) higher Spearman rank correlation coefficient and (iii) less residual trend with $\mathrm{R_e}$ than either the mass-metallicity relation and the average surface density-metallicity relation. Our result is broadly consistent with the current models of gas enrichment and metal loss. However, a more natural explanation for our findings is a local relation between the gas-phase metallicity and escape velocity.

We study the direct gas-phase oxygen abundance using the well-detected auroral [OIII]4363 line in the stacked spectra of a sample of local analogs of high-redshift galaxies. These local analogs share the same location as $z\sim2$ star-forming galaxies on the Baldwin-Phillips-Terlevich diagram. This type of analog has the same ionized interstellar medium (ISM) properties as high-redshift galaxies. We establish empirical metallicity calibrations between the direct gas-phase oxygen abundances ($7.8<12+\log(\rm{O/H})<8.4$) and the N2/O3N2 indices in our local analogs. We find significant systematic offsets between the metallicity calibrations for our local analogs of high-redshift galaxies and those derived from the local HII regions and a sample of local reference galaxies selected from the Sloan Digital Sky Survey (SDSS). The N2 and O3N2 metallicities will be underestimated by 0.05-0.1~dex relative to our calibration if one simply applies the local metallicity calibration in previous studies to high-redshift galaxies. Local metallicity calibrations also cause discrepancies of metallicity measurements in high-redshift galaxies using the N2 and O3N2 indicators. In contrast, our new calibrations produce consistent metallicities between these two indicators. We also derive metallicity calibrations for R23, O32, log([OIII]/H\beta), and log([NeIII]/[OII]) indices in our local analogs, which show significant offset compared to those in the SDSS reference galaxies. By comparing with MAPPINGS photoionization models, the different empirical metallicity calibration relations in the local analogs and the SDSS reference galaxies can be shown to be primarily due to the change of ionized ISM conditions.

We report the discovery of the ultra-luminous QSO SMSS~J215728.21-360215.1 with magnitude $z=16.9$ and W4$=7.42$ at redshift 4.75. Given absolute magnitudes of $M_{145,\rm AB}=-29.3$, $M_{300,\rm AB}=-30.12$ and $\log L_{\rm bol}/L_{\rm bol,\odot} = 14.84$, it is the QSO with the highest unlensed UV-optical luminosity currently known in the Universe. It was found by combining proper-motion data from Gaia DR2 with photometry from SkyMapper DR1 and the Wide-field Infrared Survey Explorer (WISE). In the Gaia database it is an isolated single source and thus unlikely to be strongly gravitationally lensed. It is also unlikely to be a beamed source as it is not discovered in the radio domain by either NVSS or SUMSS. It is classed as a weak-emission-line QSO and possesses broad absorption line features. A lightcurve from ATLAS spanning the time from October 2015 to December 2017 shows little sign of variability.

We present a search for bright $z\sim5$ quasars using imaging data from SkyMapper Southern Survey, Pan-STARRS1 and the Wide-field Infrared Survey Explorer (WISE). We select two sets of candidates using WISE with optical bands from SkyMapper and alternatively from Pan-STARRS1, limited to a magnitude of $i<18.2$. We follow up several candidates with spectroscopy and find that the four candidates common to both lists are quasars, while others turned out to be cool stars. Two of the four quasars, SMSS J013539.27-212628.4 at $z=4.86$ and SMSS J093032.58-221207.7 at $z=4.94$, are new discoveries and ranked among the dozen brightest known $z>4.5$ QSOs in the $i$-band.

We investigate the evolution of the ionisation parameter of star-forming galaxies using a high-redshift ($z\sim 1.5$) sample from the FMOS-COSMOS survey and matched low-redshift samples from the Sloan Digital Sky Survey. By constructing samples of low-redshift galaxies for which the stellar mass ($\mathrm{M}_{*}$), star formation rate (SFR) and specific star formation rate (sSFR) are matched to the high-redshift sample we remove the effects of an evolution in these properties. We also account for the effect of metallicity by jointly constraining the metallicity and ionisation parameter of each sample. We find an evolution in the ionisation parameter for main-sequence, star-forming galaxies and show that this evolution is driven by the evolution of sSFR. By analysing the matched samples as well as a larger sample of $z<0.3$, star-forming galaxies we show that high ionisation parameters are directly linked to high sSFRs and are not simply the byproduct of an evolution in metallicity. Our results are physically consistent with the definition of the ionisation parameter, a measure of the hydrogen ionising photon flux relative to the number density of hydrogen atoms.

Taipan is a multi-object spectroscopic galaxy survey starting in 2017 that will cover 2pi steradians over the southern sky, and obtain optical spectra for about two million galaxies out to z<0.4. Taipan will use the newly-refurbished 1.2m UK Schmidt Telescope at Siding Spring Observatory with the new TAIPAN instrument, which includes an innovative 'Starbugs' positioning system capable of rapidly and simultaneously deploying up to 150 spectroscopic fibres (and up to 300 with a proposed upgrade) over the 6-deg diameter focal plane, and a purpose-built spectrograph operating from 370 to 870nm with resolving power R>2000. The main scientific goals of Taipan are: (i) to measure the distance scale of the Universe (primarily governed by the local expansion rate, H_0) to 1% precision, and the structure growth rate of structure to 5%; (ii) to make the most extensive map yet constructed of the mass distribution and motions in the local Universe, using peculiar velocities based on improved Fundamental Plane distances, which will enable sensitive tests of gravitational physics; and (iii) to deliver a legacy sample of low-redshift galaxies as a unique laboratory for studying galaxy evolution as a function of mass and environment. The final survey, which will be completed within 5 years, will consist of a complete magnitude-limited sample (i<17) of about 1.2x10^6 galaxies, supplemented by an extension to higher redshifts and fainter magnitudes (i<18.1) of a luminous red galaxy sample of about 0.8x10^6 galaxies. Observations and data processing will be carried out remotely and in a fully-automated way, using a purpose-built automated 'virtual observer' software and an automated data reduction pipeline. The Taipan survey is deliberately designed to maximise its legacy value, by complementing and enhancing current and planned surveys of the southern sky at wavelengths from the optical to the radio.

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

We present a spectroscopic 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.

We report the detection of extended Lyman-$\alpha$ emission from the host galaxy of SDSS~J2222+2745, a strongly lensed quasar at $z = 2.8$. Spectroscopic follow-up clearly reveals extended Lyman-$\alpha$ in emission between two images of the central active galactic nucleus (AGN). We reconstruct the lensed quasar host galaxy in the source plane by applying a strong lens model to HST imaging, and resolve spatial scales as small as $\sim$200 parsecs. In the source plane we recover the host galaxy morphology to within a few hundred parsecs of the central AGN, and map the extended Lyman-$\alpha$ emission to its physical origin on one side of the host galaxy at radii $\sim$0.5-2 kpc from the central AGN. There are clear morphological differences between the Lyman-$\alpha$ and rest-frame ultraviolet stellar continuum emission from the quasar host galaxy. Furthermore, the relative velocity profiles of quasar Lyman-$\alpha$, host galaxy Lyman-$\alpha$, and metal lines in outflowing gas reveal differences in the absorbing material affecting the AGN and host galaxy. These data indicate the presence of patchy local intervening gas in front of the central quasar and its host galaxy. This interpretation is consistent with the central luminous quasar being obscured across a substantial fraction of its surrounding solid angle, resulting in strong anisotropy in the exposure of the host galaxy to ionizing radiation from the AGN. This work demonstrates the power of strong lensing-assisted studies to probe spatial scales that are currently inaccessible by other means.

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 a sample of $\sim 1000$ emission line galaxies at $z=0.4-4.7$ from the $\sim0.7$deg$^2$ High-$z$ Emission Line Survey (HiZELS) in the Bo\"otes field identified with a suite of six narrow-band filters at $\approx 0.4-2.1$ $\mu$m. These galaxies have been selected on their Ly$\alpha$ (73), {\sc [Oii]} (285), H$\beta$/{\sc [Oiii]} (387) or H$\alpha$ (362) emission-line, and have been classified with optical to near-infrared colours. A subsample of 98 sources have reliable redshifts from multiple narrow-band (e.g. [O{\sc ii}]-H$\alpha$) detections and/or spectroscopy. In this survey paper, we present the observations, selection and catalogs of emitters. We measure number densities of Ly$\alpha$, [O{\sc ii}], H$\beta$/{\sc [Oiii]} and H$\alpha$ and confirm strong luminosity evolution in star-forming galaxies from $z\sim0.4$ to $\sim 5$, in agreement with previous results. To demonstrate the usefulness of dual-line emitters, we use the sample of dual [O{\sc ii}]-H$\alpha$ emitters to measure the observed [O{\sc ii}]/H$\alpha$ ratio at $z=1.47$. The observed [O{\sc ii}]/H$\alpha$ ratio increases significantly from 0.40$\pm0.01$ at $z=0.1$ to 0.52$\pm0.05$ at $z=1.47$, which we attribute to either decreasing dust attenuation with redshift, or due to a bias in the (typically) fiber-measurements in the local Universe which only measure the central kpc regions. At the bright end, we find that both the H$\alpha$ and Ly$\alpha$ number densities at $z\approx2.2$ deviate significantly from a Schechter form, following a power-law. We show that this is driven entirely by an increasing X-ray/AGN fraction with line-luminosity, which reaches $\approx 100$ \% at line-luminosities $L\gtrsim3\times10^{44}$ erg s$^{-1}$.

We investigate the energy sources of random turbulent motions of ionised gas from H$\alpha$ emission in eight local star-forming galaxies from the Sydney-AAO Multi-object Integral field spectrograph (SAMI) Galaxy Survey. These galaxies satisfy strict pure star-forming selection criteria to avoid contamination from active galactic nuclei (AGN) or strong shocks/outflows. Using the relatively high spatial and spectral resolution of SAMI, we find that -- on sub-kpc scales our galaxies display a flat distribution of ionised gas velocity dispersion as a function of star formation rate (SFR) surface density. A major fraction of our SAMI galaxies shows higher velocity dispersion than predictions by feedback-driven models, especially at the low SFR surface density end. Our results suggest that additional sources beyond star formation feedback contribute to driving random motions of the interstellar medium (ISM) in star-forming galaxies. We speculate that gravity, galactic shear, and/or magnetorotational instability (MRI) may be additional driving sources of turbulence in these galaxies.

Stars form in cold molecular clouds. However, molecular gas is difficult to observe because the most abundant molecule (H2) lacks a permanent dipole moment. Rotational transitions of CO are often used as a tracer of H2, but CO is much less abundant and the conversion from CO intensity to H2 mass is often highly uncertain. Here we present a new method for estimating the column density of cold molecular gas (Sigma_gas) using optical spectroscopy. We utilise the spatially resolved H-alpha maps of flux and velocity dispersion from the Sydney-AAO Multi-object Integral-field spectrograph (SAMI) Galaxy Survey. We derive maps of Sigma_gas by inverting the multi-freefall star formation relation, which connects the star formation rate surface density (Sigma_SFR) with Sigma_gas and the turbulent Mach number (Mach). Based on the measured range of Sigma_SFR = 0.005-1.5 M_sol/yr/kpc^2 and Mach = 18-130, we predict Sigma_gas = 7-200 M_sol/pc^2 in the star-forming regions of our sample of 260 SAMI galaxies. These values are close to previously measured Sigma_gas obtained directly with unresolved CO observations of similar galaxies at low redshift. We classify each galaxy in our sample as 'Star-forming' (219) or 'Composite/AGN/Shock' (41), and find that in Composite/AGN/Shock galaxies the average Sigma_SFR, Mach, and Sigma_gas are enhanced by factors of 2.0, 1.6, and 1.3, respectively, compared to Star-forming galaxies. We compare our predictions of Sigma_gas with those obtained by inverting the Kennicutt-Schmidt relation and find that our new method is a factor of two more accurate in predicting Sigma_gas, with an average deviation of 32% from the actual Sigma_gas.

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.

We revisit the evolution of the mass-metallicity relation of low- and high-redshift galaxies by using a sample of local analogs of high-redshift galaxies. These analogs share the same location of the UV-selected star-forming galaxies at $z\sim2$ on the [OIII]/H$\beta$ versus [NII]/H$\alpha$ nebular emission-line diagnostic (or BPT) diagram. Their physical properties closely resemble those in $z\sim2$ UV-selected star-forming galaxies being characterized in particular by high ionization parameters ($\log q\approx7.9$) and high electron densities ($n_e\approx100~\rm{cm}^{-3}$). With the full set of well-detected rest-frame optical diagnostic lines, we measure the gas-phase oxygen abundance in the SDSS galaxies and these local analogs using the empirical relations and the photoionization models. We find that the metallicity difference between the SDSS galaxies and our local analogs in the $8.5<log(M_*/M_{\odot})<9.0$ stellar mass bin varies from -0.09 to 0.39 dex, depending on strong-line metallicity measurement methods. Due to this discrepancy the evolution of mass-metallicity should be used to compare with the cosmological simulations with caution. We use the [SII]/H$\alpha$ and [OI]/H$\alpha$ BPT diagram to reduce the potential AGN and shock contamination in our local analogs. We find that the AGN/shock influences are negligible on the metallicity estimation.

Star-forming galaxies at $z > 1$ exhibit significantly different properties to local galaxies of equivalent stellar mass. Not only are high-redshift star-forming galaxies characterized by higher star formation rates and gas fractions than their local counterparts, they also appear to host star-forming regions with significantly different physical conditions, including greater electron densities. To understand what physical mechanisms are responsible for the observed evolution of star-forming conditions we have assembled the largest sample of star-forming galaxies at $z\sim 1.5$ with emission-line measurements of the $\mathrm{[OII]} \lambda \lambda 3726,3729$ doublet. By comparing our $z\sim 1.5$ sample to local galaxy samples with equivalent distributions of stellar mass, star formation rate and specific star formation rate we investigate the proposed evolution in electron density and its dependence on global properties. We measure an average electron density of $114_{-27}^{+28} \, \mathrm {cm}^{-3} $ for our $z\sim 1.5$ sample, a factor of five greater than the typical electron density of local star-forming galaxies. However, we find no offset between the typical electron densities of local and high-redshift galaxies with equivalent star-formation rates. Our work indicates that the average electron density of a sample is highly sensitive to the star formation rates, implying that the previously observed evolution is mainly the result of selection effects.

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.

This paper presents the first major data release and survey description for the ANU WiFeS SuperNovA Program (AWSNAP). AWSNAP is an ongoing supernova spectroscopy campaign utilising the Wide Field Spectrograph (WiFeS) on the Australian National University (ANU) 2.3m telescope. The first and primary data release of this program (AWSNAP-DR1) releases 357 spectra of 175 unique objects collected over 82 equivalent full nights of observing from July 2012 to August 2015. These spectra have been made publicly available via the WISeREP supernova spectroscopy repository. We analyse the AWSNAP sample of Type Ia supernova spectra, including measurements of narrow sodium absorption features afforded by the high spectral resolution of the WiFeS instrument. In some cases we were able to use the integral-field nature of the WiFeS instrument to measure the rotation velocity of the SN host galaxy near the SN location in order to obtain precision sodium absorption velocities. We also present an extensive time series of SN 2012dn, including a near-nebular spectrum which both confirms its "super-Chandrasekhar" status and enables measurement of the sub-solar host metallicity at the SN site.

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

We present a sample of local analogs for high-redshift galaxies selected in the Sloan Digital Sky Survey (SDSS). The physical conditions of the interstellar medium (ISM) in these local analogs resemble those in high-redshift galaxies. These galaxies are selected based on their positions in the [OIII]/H$\beta$ versus [NII]/H$\alpha$ nebular emission-line diagnostic diagram. We show that these local analogs share similar physical properties with high-redshift galaxies, including high specific star formation rates (sSFRs), flat UV continuums and compact galaxy sizes. In particular, the ionization parameters and electron densities in these analogs are comparable to those in $z\simeq2-3$ galaxies, but higher than those in normal SDSS galaxies by $\simeq$0.6~dex and $\simeq$0.9~dex, respectively. The mass-metallicity relation (MZR) in these local analogs shows $-0.2$~dex offset from that in SDSS star-forming galaxies at the low mass end, which is consistent with the MZR of the $z\sim2-3$ galaxies. We compare the local analogs in this study with those in other studies, including Lyman break analogs (LBA) and green pea (GP) galaxies. The analogs in this study share a similar star formation surface density with LBAs, but the ionization parameters and electron densities in our analogs are higher than those in LBAs by factors of 1.5 and 3, respectively. The analogs in this study have comparable ionization parameter and electron density to the GP galaxies, but our method can select galaxies in a wider redshift range. We find the high sSFR and SFR surface density can increase the electron density and ionization parameters, but still cannot fully explain the difference in ISM condition between nearby galaxies and the local analogs/high-redshift galaxies.

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.

Modern cosmology predicts that a galaxy overdensity is associated to a large reservoir of the intergalactic gas, which can be traced by the Ly$\alpha$ forest absorption. We have undertaken a systematic study of the relation between Coherently Strong intergalactic Ly$\alpha$ Absorption systems (CoSLAs), which have highest optical depth ($\tau$) in $\tau$ distribution, and mass overdensities on the scales of $\sim$ 10 - 20 $h^{-1}$ comoving Mpc. On such large scales, our cosmological simulations show a strong correlation between the effective optical depth ($\tau_{\rm{eff}}$) of the CoSLAs and the 3-D mass overdensities. In moderate signal-to-noise spectra, however, the profiles of CoSLAs can be confused with high column density absorbers. For $z>2.6$, where the corresponding Ly$\beta$ is redshifted to the optical, we have developed the technique to differentiate between these two alternatives. We have applied this technique to SDSS-III quasar survey at $z = 2.6$ - 3.3, and we present a sample of five CoSLA candidates with $\tau_{\rm{eff}}$ on 15 $h^{-1}$ Mpc greater than $4.5\times$ the mean optical depth. At lower redshifts of $z < 2.6$, where the background quasar density is higher, the overdensity can be traced by intergalactic absorption groups using multiple sight lines. Our overdensity searches fully utilize the current and next generation of Ly$\alpha$ forest surveys which cover a survey volume of $> (1\ h^{-1}$ Gpc)$^3$. In addition, systems traced by CoSLAs will build a uniform sample of the most massive overdensities at $z > 2$ to constrain the models of structure formation, and offer a unique laboratory to study the interactions between galaxy overdensities and the intergalactic medium.

Using a sample of ~410 000 galaxies to depth I_AB = 24 over 8.26 deg^2 in the Bootes field (~10 times larger than z~1 luminosity function studies in the prior literature), we have accurately measured the evolving B-band luminosity function of red galaxies at z<1.2 and blue galaxies at z<1.0. In addition to the large sample size, we utilise photometry that accounts for the varying angular sizes of galaxies, photometric redshifts verified with spectroscopy, and absolute magnitudes that should have very small random and systematic errors. Our results are consistent with the migration of galaxies from the blue cloud to the red sequence as they cease to form stars, and with downsizing in which more massive and luminous blue galaxies cease star formation earlier than fainter less massive ones. Comparing the observed fading of red galaxies with that to be expected from passive evolution alone, we find that the stellar mass contained within the red galaxy population has increased by a factor of ~3.6 from z~1.1 to z~0.1. The bright end of the red galaxy luminosity function fades with decreasing redshift, the rate of fading increasing from ~0.2 mag per unit redshift at z = 1.0 to ~0.8 at z = 0.2. The overall decrease in luminosity implies that the stellar mass in individual highly luminous red galaxies increased by a factor of ~2.2 from z = 1.1 to z = 0.1.

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.

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.

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

We present new near-infrared (IR) observations of the H$\beta\ \lambda4861$ and MgII $\lambda2798$ lines for 32 luminous quasars with $3.2<z<3.9$ using the Palomar Hale 200 inch telescope and the Large Binocular Telescope. We find that the MgII Full Width at Half Maximum (FWHM) is well correlated with the H$\beta$ FWHM, confirming itself as a good substitute for the H$\beta$ FWHM in the black hole mass estimates. The continuum luminosity at 5100 \AA\ well correlates with the continuum luminosity at 3000 \AA\ and the broad emission line luminosities (H$\beta$ and MgII). With simultaneous near-IR spectroscopy of the H$\beta$ and MgII lines to exclude the influences of flux variability, we are able to evaluate the reliability of estimating black hole masses based on the MgII line for high redshift quasars. With the reliable H$\beta$ line based black hole mass and Eddington ratio estimates, we find that the $z\sim3.5$ quasars in our sample have black hole masses $1.90\times10^{9} M_{\odot} \lesssim M_{\rm BH} \lesssim 1.37\times10^{10} M_{\odot}$, with a median of $\sim 5.14\times10^{9} M_{\odot}$ and are accreting at Eddington ratios between 0.30 and 3.05, with a median of $\sim1.12$. Assuming a duty cycle of 1 and a seed black hole mass of $10^{4} M_{\odot}$, we show that the $z\sim3.5$ quasars in this sample can grow to their estimated black hole masses within the age of the Universe at their redshifts.

Fast radio bursts (FRBs) are one of the most tantalizing mysteries of the radio sky; their progenitors and origins remain unknown and until now no rapid multiwavelength follow-up of an FRB has been possible. New instrumentation has decreased the time between observation and discovery from years to seconds, and enables polarimetry to be performed on FRBs for the first time. We have discovered an FRB (FRB 140514) in real-time on 14 May, 2014 at 17:14:11.06 UTC at the Parkes radio telescope and triggered follow-up at other wavelengths within hours of the event. FRB 140514 was found with a dispersion measure (DM) of 562.7(6) cm$^{-3}$ pc, giving an upper limit on source redshift of $z \lesssim 0.5$. FRB 140514 was found to be 21$\pm$7% (3-$\sigma$) circularly polarized on the leading edge with a 1-$\sigma$ upper limit on linear polarization $<10%$. We conclude that this polarization is intrinsic to the FRB. If there was any intrinsic linear polarization, as might be expected from coherent emission, then it may have been depolarized by Faraday rotation caused by passing through strong magnetic fields and/or high density environments. FRB 140514 was discovered during a campaign to re-observe known FRB fields, and lies close to a previous discovery, FRB 110220; based on the difference in DMs of these bursts and time-on-sky arguments, we attribute the proximity to sampling bias and conclude that they are distinct objects. Follow-up conducted by 12 telescopes observing from X-ray to radio wavelengths was unable to identify a variable multiwavelength counterpart, allowing us to rule out models in which FRBs originate from nearby ($z < 0.3$) supernovae and long duration gamma-ray bursts.

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.

In merger-driven models of massive galaxy evolution, the luminous quasar phase is expected to be accompanied by vigorous star formation in quasar host galaxies. In this paper, we use high column density Damped Lyman Alpha (DLA) systems along quasar sight lines as natural coronagraphs to directly study the far-UV (FUV) radiation from the host galaxies of luminous background quasars. We have stacked the spectra of $\sim$2,000 DLA systems (N_HI>10^{20.6} cm^{-2}) with a median absorption redshift < z > = 2.6 selected from quasars observed in the SDSS-III Baryon Oscillation Spectroscopic Survey. We detect residual flux in the dark troughs of the composite DLA spectra. The level of this residual flux significantly exceeds systematic errors in the SDSS fiber sky subtraction; furthermore, the residual flux is strongly correlated with the continuum luminosity of the background quasar, while uncorrelated with DLA column density or metallicity. We conclude that the flux could be associated with the average FUV radiation from the background quasar host galaxies (with medium redshift < z > = 3.1) that is not blocked by the intervening DLA. Assuming all of the detected flux originates from quasar hosts, for the highest quasar luminosity bin (< L >= 2.5x 10^{13} L_sun), the host galaxy has a FUV intensity of 1.5 +/- 0.2 x 10^{40} erg s^{-1} A^{-1}; this corresponds to an unobscured UV star formation rate of 9 M_sun/yr.

A small fraction of quasars have long been known to show bulk velocity offsets in the broad Balmer lines with respect to the systemic redshift of the host galaxy. Models to explain these offsets usually invoke broad-line region gas kinematics/asymmetry around single black holes (BHs), orbital motion of massive (~sub-pc) binary black holes (BBHs), or recoil BHs, but single-epoch spectra are unable to distinguish between these scenarios. The line-of-sight (LOS) radial velocity (RV) shifts from long-term spectroscopic monitoring can be used to test the BBH hypothesis. We have selected a sample of 399 quasars with offset broad H-beta lines from the SDSS DR7 quasar catalog, and have conducted second-epoch optical spectroscopy for 50 of them. Combined with the existing SDSS spectra, the new observations enable us to constrain the LOS RV shifts of broad H-beta lines with a rest-frame baseline of a few years to nearly a decade. Using cross-correlation analysis, we detect significant (99% confidence) radial accelerations in the broad H-beta lines in 24 of the 50 objects. We suggest that 9 of the 24 detections are sub-pc BBH candidates, which show consistent velocity shifts independently measured from a second broad line (either H-alpha or Mg II) without significant changes in the broad-line profiles. Combining the results on the general quasar population studied in Paper I, we find a tentative anti-correlation between the velocity offset in the first-epoch spectrum and the average acceleration between two epochs, which could be explained by orbital phase modulation when the time separation between two epochs is a non-negligible fraction of the orbital period of the motion causing the line displacement. We discuss the implications of our results for the identification of sub-pc BBH candidates in offset-line quasars and for the constraints on their frequency and orbital parameters. [Abridged]

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

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

This document summarizes the results of a community-based discussion of the potential science impact of the Mayall+BigBOSS highly multiplexed multi-object spectroscopic capability. The KPNO Mayall 4m telescope equipped with the DOE- and internationally-funded BigBOSS spectrograph offers one of the most cost-efficient ways of accomplishing many of the pressing scientific goals identified for this decade by the "New Worlds, New Horizons" report. The BigBOSS Key Project will place unprecedented constraints on cosmological parameters related to the expansion history of the universe. With the addition of an open (publicly funded) community access component, the scientific impact of BigBOSS can be extended to many important astrophysical questions related to the origin and evolution of galaxies, stars, and the IGM. Massive spectroscopy is the critical missing ingredient in numerous ongoing and planned ground- and space-based surveys, and BigBOSS is unique in its ability to provide this to the US community. BigBOSS data from community-led projects will play a vital role in the education and training of students and in maintaining US leadership in these fields of astrophysics. We urge the NSF-AST division to support community science with the BigBOSS multi-object spectrograph through the period of the BigBOSS survey in order to ensure public access to the extraordinary spectroscopic capability.

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.

We present spectroscopic observations in the rest-frame optical and near- to mid-infrared wavelengths of four gravitationally lensed infrared (IR) luminous star-forming galaxies at redshift 1 < z < 3 from the LUCIFER instrument on the Large Binocular Telescope and the Infrared Spectrograph on Spitzer. The sample was selected to represent pure, actively star-forming systems, absent of active galactic nuclei. The large lensing magnifications result in high signal-to-noise spectra that can probe faint IR recombination lines, including Pa-alpha and Br-alpha at high redshifts. The sample was augmented by three lensed galaxies with similar suites of unpublished data and observations from the literature, resulting in the final sample of seven galaxies. We use the IR recombination lines in conjunction with H-alpha observations to probe the extinction, Av, of these systems, as well as testing star formation rate (SFR) indicators against the SFR measured by fitting spectral energy distributions to far-IR photometry. Our galaxies occupy a range of Av from ~0 to 5.9 mag, larger than previously known for a similar range of IR luminosities at these redshifts. Thus, estimates of SFR even at z ~ 2 must take careful count of extinction in the most IR luminous galaxies. We also measure extinction by comparing SFR estimates from optical emission lines with those from far-IR measurements. The comparison of results from these two independent methods indicates a large variety of dust distribution scenarios at 1 < z < 3. Without correcting for dust extinction, the H-alpha SFR indicator underestimates the SFR; the size of the necessary correction depends on the IR luminosity and dust distribution scenario. Individual SFR estimates based on the 6.2 micron PAH emission line luminosity do not show a systematic discrepancy with extinction, although a considerable, ~0.2 dex scatter is observed.

We present first results from the 3D-HST program, a near-IR spectroscopic survey performed with the Wide Field Camera 3 on the Hubble Space Telescope. We have used 3D-HST spectra to measure redshifts and Halpha equivalent widths for a stellar mass-limited sample of 34 galaxies at 1<z<1.5 with M(stellar)>10^11 M(sun) in the COSMOS, GOODS, and AEGIS fields. We find that a substantial fraction of massive galaxies at this epoch are forming stars at a high rate: the fraction of galaxies with Halpha equivalent widths >10 A is 59%, compared to 10% among SDSS galaxies of similar masses at z=0.1. Galaxies with weak Halpha emission show absorption lines typical of 2-4 Gyr old stellar populations. The structural parameters of the galaxies, derived from the associated WFC3 F140W imaging data, correlate with the presence of Halpha: quiescent galaxies are compact with high Sersic index and high inferred velocity dispersion, whereas star-forming galaxies are typically large two-armed spiral galaxies, with low Sersic index. Some of these star forming galaxies might be progenitors of the most massive S0 and Sa galaxies. Our results challenge the idea that galaxies at fixed mass form a homogeneous population with small scatter in their properties. Instead we find that massive galaxies form a highly diverse population at z>1, in marked contrast to the local Universe.

Type IIn SNe show spectral evidence for strong interaction between their blast wave and dense circumstellar material (CSM) around the progenitor star. SN2010jl was the brightest core-collapse SN in 2010, and it was a Type IIn explosion with strong CSM interaction. Andrews et al. recently reported evidence for an IR excess in SN2010jl, indicating either new dust formation or the heating of CSM dust in an IR echo. Here we report multi-epoch spectra of SN2010jl that reveal the tell-tale signature of new dust formation: emission-line profiles becoming systematically more blueshifted as the red side of the line is blocked by increasing extinction. The effect is seen clearly in the intermediate-width (400--4000 km/s) component of H$\alpha$ beginning roughly 30d after explosion. Moreover, we present near-IR spectra demonstrating that the asymmetry in the hydrogen-line profiles is wavelength dependent, appearing more pronounced at shorter wavelengths. This evidence suggests that new dust grains had formed quickly in the post-shock shell of SN 2010jl arising from CSM interaction. Since the observed dust temperature has been attributed to an IR echo and not to new dust, either (1) IR excess emission at $\lambda < 5 \mu$m is not a particularly sensitive tracer of new dust formation in SNe, or (2) some assumptions about expected dust temperatures might require further study. Lastly, we discuss one possible mechanism other than dust that might lead to increasingly blueshifted line profiles in SNeIIn, although the wavelength dependence of the asymmetry argues against this hypothesis in the case of SN2010jl.

The He II \lambda 1640 emission line has been suggested as a direct probe of Population III (Pop III) stars at high-redshift, since it can arise from highly energetic ionizing photons associated with hot, metal free stars. We use the HST WFC3/F130N IR narrowband filter to probe He II \lambda 1640 emission in galaxy IOK-1 at z=6.96. The sensitivity of this measurement is >5x deeper than for previous measurements. From this deep narrowband imaging, combined with broadband observations in the F125W and F160W filters, we find the He II flux to be 1.2+/- 1.0x 10^-18 ergs/s/cm^2, corresponding to a 1\sigma upper limit on the Pop III star formation rate (SFR) of ~ 0.5 M_sun/yr for the case of a Salpeter IMF with 50-500M_sun and mass loss. Given that the broadband measurements can be fit with a UV continuum spectral flux density of ~ 4.85x 10^-10x \lambda ^-2.46 ergs/s/cm^2/A, which corresponds to an overall SFR of ~16+/-2.6 M_sun/yr, massive Pop III stars represent < 6% of the total star formation. This measurement places the strongest limit yet on metal-free star formation at high redshift, although the exact conversion from He II luminosity to Pop III SFR is highly uncertain due to the unknown IMF, stellar evolution, and photoionization effects. Although we have not detected He II \lambda 1640 at more than the 1.2\sigma level, our work suggests that a > 3\sigma level detection is possible with JWST.

The identifications of quasars in the redshift range 2.2<z<3 are known to be very inefficient as their optical colors are indistinguishable from those of stars. Recent studies have proposed to use optical variability or near-IR colors to improve the identifications of the missing quasars in this redshift range. Here we present a case study by combining both factors. We select a sample of 70 quasar candidates from variables in SDSS Stripe 82, which are non-UV excess sources and have UKIDSS near-IR public data. They are clearly separated into two parts on the Y-K/g-z color-color diagram, and 59 of them meet or lie close to a newly proposed Y-K/g-z selection criterion for z<4 quasars. 44 of these 59 sources have been previously identified as quasars in SDSS DR7, and 35 among them are quasars at 2.2<z<3. We present spectroscopic observations of 14 of 15 remaining quasar candidates using the Bok 2.3m telescope and the MMT 6.5m telescope, and successfully identify all of them as new quasars at z=2.36 to 2.88. We also apply this method to a sample of 643 variable quasar candidates with SDSS-UKIDSS nine-band photometric data selected from 1875 new quasar candidates in SDSS Stripe 82 given by Butler & Bloom based on the time-series selections, and find that 188 of them are probably new quasars with photometric redshifts at 2.2<z<3. Our results indicate that the combination of optical variability and optical/near-IR colors is probably the most efficient way in finding 2.2<z<3 quasars and very helpful for constructing a complete quasar sample. We discuss its implications to the ongoing and upcoming large optical and near-IR sky surveys.

[Abridged] We investigate the average physical properties and star formation histories of the most UV-luminous star-forming galaxies at z~3.7. Our results are derived from analyses of the average spectral energy distributions (SEDs), constructed from stacked optical to infrared photometry, of a sample of the 1,902 most UV-luminous star-forming galaxies found in 5.3 square degrees of the NOAO Deep Wide-Field Survey. We bin the sample according to UV luminosity, and find that the shape of the average SED in the rest-frame optical and infrared is fairly constant with UV luminosity: i.e., more UV luminous galaxies are, on average, also more luminous at longer wavelengths. In the rest-UV, however, the spectral slope (measured at 0.13-0.28 um) rises steeply with the median UV luminosity from -1.8 at L L* to -1.2 in the brightest bin (L~4-5L*). We use population synthesis analyses to derive the average physical properties of these galaxies and find that: (1) L_UV, and thus star formation rates (SFRs), scale closely with stellar mass such that more UV-luminous galaxies are also more massive; (2) The median ages indicate that the stellar populations are relatively young (200-400 Myr) and show little correlation with UV luminosity; and (3) More UV-luminous galaxies are dustier than their less-luminous counterparts, such that L~4-5L* galaxies are extincted up to A(1600)=2 mag while L L* galaxies have A(1600)=0.7-1.5 mag. Based on these observations, we argue that the average star formation histories of UV-luminous galaxies are better described by models in which SFR increases with time in order to simultaneously reproduce the tight correlation between the observed SFR and stellar mass, and the universally young ages of these galaxies. We demonstrate the potential of measurements of the SFR-M* relation at multiple redshifts to discriminate between simple models of star formation histories.

We present the first spectroscopic measurements of the [OIII] 5007 A line in two z ~ 3.1 Lyman-alpha emitting galaxies (LAEs) using the new near-infrared instrument LUCIFER1 on the 8.4m Large Binocular Telescope (LBT). We also describe the optical imaging and spectroscopic observations used to identify these Lya emitting galaxies. Using the [OIII] line we have measured accurate systemic redshifts for these two galaxies, and discovered a velocity offset between the [OIII] and Ly-alpha lines in both, with the Lya line peaking 342 and 125 km/s redward of the systemic velocity. These velocity offsets imply that there are powerful outflows in high-redshift LAEs. They also ease the transmission of Lya photons through the interstellar medium and intergalactic medium around the galaxies. By measuring these offsets directly, we can refine both Lya-based tests for reionization, and Lya luminosity function measurements where the Lya forest affects the blue wing of the line. Our work also provides the first direct constraints on the strength of the [OIII] line in high-redshift LAEs. We find [OIII] fluxes of 7 and 36 x 10^-17 erg s^-1 cm^-2 in two z ~ 3.1 LAEs. These lines are strong enough to dominate broad-band flux measurements that include the line (in thiscase, K_s band photometry). Spectral energy distribution fits that do not account for the lines would therefore overestimate the 4000 A (and/or Balmer) break strength in such galaxies, and hence also the ages and stellar masses of such high-z galaxies.

We present rest-frame optical images and spectra of the gravitationally lensed, star-forming galaxy J0900+2234 (z=2.03). The observations were performed with the newly commissioned LUCIFER1 near-infrared instrument mounted on the Large Binocular Telescope (LBT). We fit lens models to the rest-frame optical images and find the galaxy has an intrinsic effective radius of 7.4 kpc with a lens magnification factor of about 5 for the A and B components. We also discovered a new arc belonging to another lensed high-z source galaxy, which makes this lens system a potential double Einstein ring system. Using the high S/N rest-frame optical spectra covering H+K band, we detected Hbeta, OIII, Halpha, NII and SII emission lines. Detailed physical properties of this high-z galaxy were derived. The extinction towards the ionized HII regions (E_g(B-V)) is computed from the flux ratio of Halpha and Hbeta and appears to be much higher than that towards stellar continuum (E_s(B-V)), derived from the optical and NIR broad band photometry fitting. The metallicity was estimated using N2 and O3N2 indices. It is in the range of 1/5-1/3 solar abundance, which is much lower than the typical z~2 star-forming galaxies. From the flux ratio of SII 6717 and 6732, we found that the electron number density of the HII regions in the high-z galaxy were >1000 cm^-3, consistent with other z~2 galaxies but much higher than that in local HII regions. The star-formation rate was estimated via the Halpha luminosity, after correction for the lens magnification, to be about 365\pm69 Msun/yr. Combining the FWHM of Halpha emission lines and the half-light radius, we found the dynamical mass of the lensed galaxy is 5.8\pm0.9x10^10 Msun. The gas mass is 5.1\pm1.1x10^10~Msun from the H\alpha flux surface density by using global Kennicutt-Schmidt Law, indicating a very high gas fraction of 0.79\pm0.19 in J0900+2234.

We report the discovery of a gravitationally lensed quasar identified serendipitously in the Sloan Digital Sky Survey (SDSS). The object, SDSS J094604.90+183541.8, was initially targeted for spectroscopy as a luminous red galaxy, but the SDSS spectrum has the features of both a z=0.388 galaxy and a z=4.8 quasar. We have obtained additional imaging that resolves the system into two quasar images separated by 3.06 arcsec and a bright galaxy that is strongly blended with one of the quasar images. We confirm spectroscopically that the two quasar images represent a single lensed source at z=4.8 with a total magnification of 3.2, and we derive a model for the lensing galaxy. This is the highest redshift lensed quasar currently known. We examine the issues surrounding the selection of such an unusual object from existing data and briefly discuss implications for lensed quasar surveys.

We present the discovery of six new quasars at z~6 selected from the Sloan Digital Sky Survey (SDSS) southern survey, a deep imaging survey obtained by repeatedly scanning a stripe along the celestial equator. The six quasars are about two magnitudes fainter than the luminous z~6 quasars found in the SDSS main survey and one magnitude fainter than the quasars reported in Paper I (Jiang et al. 2008). Four of them comprise a complete flux-limited sample at 21<z_AB<21.8 over an effective area of 195 deg^2. The other two quasars are fainter than z_AB=22 and are not part of the complete sample. The quasar luminosity function at z~6 is well described as a single power law \Phi(L_{1450}) \propto L_{1450}^{\beta} over the luminosity range -28<M_{1450}<-25. The best-fitting slope \beta varies from -2.6 to -3.1, depending on the quasar samples used, with a statistical error of 0.3-0.4. About 40% of the quasars discovered in the SDSS southern survey have very narrow Lya emission lines, which may indicate small black hole masses and high Eddington luminosity ratios, and therefore short black hole growth time scales for these faint quasars at early epochs.

We explore MHD solutions for envelope expansions with core collapse (EECC) with isothermal MHD shocks in a quasi-spherical symmetry and outline potential astrophysical applications of such magnetized shock flows. MHD shock solutions are classified into three classes according to the downstream characteristics near the core. Class I solutions are those characterized by free-fall collapses towards the core downstream of an MHD shock, while Class II solutions are those characterized by Larson-Penston (LP) type near the core downstream of an MHD shock. Class III solutions are novel, sharing both features of Class I and II solutions with the presence of a sufficiently strong magnetic field as a prerequisite. Various MHD processes may occur within the regime of these isothermal MHD shock similarity solutions, such as sub-magnetosonic oscillations, free-fall core collapses, radial contractions and expansions. We can also construct families of twin MHD shock solutions as well as an `isothermal MHD shock' separating two magnetofluid regions of two different yet constant temperatures. The versatile behaviours of such MHD shock solutions may be utilized to model a wide range of astrophysical problems, including star formation in magnetized molecular clouds, MHD link between the asymptotic giant branch phase to the proto-planetary nebula phase with a hot central magnetized white dwarf, relativistic MHD pulsar winds in supernova remnants, radio afterglows of soft gamma-ray repeaters and so forth.

We explore self-similar dynamical processes in a spherical isothermal self-gravitational fluid with an emphasis on shocks and outline astrophysical applications of such shock solutions. The previous similarity shock solutions of Tsai & Hsu and of Shu et al. may be classified into two types: Class I solutions with downstream being free-fall collapses and Class II solutions with downstream being Larson-Penston (LP) type solutions. By the analyses of Lou & Shen and Shen & Lou, we further construct similarity shock solutions in the `semi-complete space'. These general shock solutions can accommodate and model dynamical processes of radial outflows (wind), inflows (accretion or contraction), subsonic oscillations, and free-fall core collapses all with shocks in various settings such as star-forming molecular clouds, `champagne flows' in H{\sevenrm II} regions around luminous massive OB stars or surrounding quasars, dynamical connection between the asymptotic giant branch phase to the proto-planetary nebula phase with a central hot white dwarf as well as accretion shocks around compact objects such as white dwarfs, neutron stars, and supermassive black holes. By a systematic exploration, we are able to construct families of infinitely many discrete Class I and Class II solutions matching asymptotically with a static outer envelope of singular isothermal sphere; the shock solutions of Tsai & Hsu form special subsets. These similarity shocks travel at either subsonic or supersonic constant speeds. We also construct twin shocks as well as an `isothermal shock' separating two fluid regions of two different yet constant temperatures.

It has been known for two decades that a tight correlation exists between far-infrared (FIR) and radio (1.4 and 4.8 GHz) global fluxes/luminosities from galaxies, which may be explained in terms of massive star formation activities in these galaxies. For this very reason, a correlation might also exist between X-ray and FIR/radio global luminosities of galaxies. We analyze data from {\it ROSAT} All-Sky Survey and {\it IRAS} to show that such correlation does indeed exist between FIR (42.5$-122.5\mu$m) and soft X-ray (0.1$-$2.4keV) luminosities in active star-forming galaxies (ASFGs).In order to establish a physical connection between the $L_{\rmX}-L_{\rm FIR}$ correlation and global star formation rate (SFR) in galaxies, we empirically derive both $L_{\rm X}-L_{\rm B}$ and $L_{\rm FIR}-L_{\rm B}$ relations. Futhermore, we propose a relation between soft X-ray luminosity and SFR in star-forming galaxies. To further understand the $L_{\rm X}-L_{\rm FIR}$ relation, we construct an empirical model in which both FIR and X-ray emissions consist of two components with one being closely associated with star formation and the other not.