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.

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 identify a strong Ly-alpha damping wing profile in the spectrum of the quasar P183+05 at z=6.4386. Given the detection of several narrow metal absorption lines at z=6.40392, the most likely explanation for the absorption profile is that it is due to a damped Ly-alpha system. However, in order to match the data a contribution of an intergalactic medium 5-38% neutral or additional weaker absorbers near the quasar is also required. The absorption system presented here is the most distant damped Ly-alpha system currently known. We estimate an HI column density ($10^{20.68\pm0.25}\,$cm$^{-2}$), metallicity ([O/H]$=-2.92\pm 0.32$), and relative chemical abundances of a system consistent with a low-mass galaxy during the first Gyr of the universe. This object is among the most metal-poor damped Ly-alpha systems known and, even though it is observed only ~850 Myr after the big bang, its relative abundances do not show signatures of chemical enrichment by Population III stars.

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

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

We present deep MUSE observations of five quasars within the first Gyr of the Universe ($z\gtrsim6$), four of which display extended Ly$\alpha$ halos. After PSF-subtraction, we reveal halos surrounding two quasars for the first time, as well as confirming the presence of two more halos for which tentative detections exist in long-slit spectroscopic observations and narrow-band imaging. The four Ly$\alpha$ halos presented here are diverse in morphology and size, they each display spatial asymmetry, and none are centred on the position of the quasar. Spectra of the diffuse halos demonstrate that none are dramatically offset in velocity from the systemic redshift of the quasars ($\Delta$ v $< 200$ kms$^{-1}$), however each halo shows a broad Ly$\alpha$ line, with a velocity width of order $\sim1000$ kms$^{-1}$. Total Ly$\alpha$ luminosities range between $\sim$ $2 \times 10^{43}$ erg s$^{-1}$ and $\sim$ $2 \times 10^{44}$ erg s$^{-1}$, reaching maximum radial extents of $13 - 30$ pkpc from the quasar positions. We find larger sizes and higher Ly$\alpha$ luminosities than previous literature results at this redshift, but find no correlation between the quasar properties and the Ly$\alpha$ halo, suggesting that the detected emission is most closely related to the physical properties of the circum-galactic medium

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

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

Luminous high-redshift quasars can be used to probe of the intergalactic medium (IGM) in the early universe because their UV light is absorbed by the neutral hydrogen along the line of sight. They help us to measure the neutral hydrogen fraction of the high-z universe, shedding light on the end of reionization epoch. In this paper, we present a discovery of a new quasar (PSO J006.1240+39.2219) at redshift $z=6.61\pm0.02$ from Panoramic Survey Telescope & Rapid Response System 1. Including this quasar, there are nine quasars above $z>6.5$ up to date. The estimated continuum brightness is $M_\text{1450}$=$-25.96\pm0.08$. PSO J006.1240+39.2219 has a strong Ly~$\alpha$ emission compared with typical low-redshift quasars, but the measured near-zone region size is $R_\text{NZ}=3.2\pm1.1$ proper megaparsecs, which is consistent with other quasars at z$\sim$6.

We present ALMA band 8 observations of the [OIII] 88um line and the underlying thermal infrared continuum emission in the z=6.08 quasar CFHQS J2100-1715 and its dust-obscured starburst companion galaxy (projected distance: ~60 kpc). Each galaxy hosts dust-obscured star formation at rates > 100 M_sun/yr, but only the quasar shows evidence for an accreting 10^9 M_sun black hole. Therefore we can compare the properties of the interstellar medium in distinct galactic environments in two physically associated objects, ~1 Gyr after the Big Bang. Bright [OIII] 88um emission from ionized gas is detected in both systems; the positions and line-widths are consistent with earlier [CII] measurements, indicating that both lines trace the same gravitational potential on galactic scales. The [OIII] 88um/FIR luminosity ratios in both sources fall in the upper range observed in local luminous infrared galaxies of similar dust temperature, although the ratio of the quasar is smaller than in the companion. This suggests that gas ionization by the quasar (expected to lead to strong optical [OIII] 5008A emission) does not dominantly determine the quasar's FIR [OIII] 88um luminosity. Both the inferred number of photons needed for the creation of O++ and the typical line ratios can be accounted for without invoking extreme (top-heavy) stellar initial mass functions in the starbursts of both sources.

We present high angular resolution imaging ($23.9 \times 11.3$ mas, $138.6 \times 65.5$ pc) of the radio-loud quasar PSO~J352.4034$-$15.3373 at $z=5.84$ with the Very Long Baseline Array (VLBA) at 1.54 GHz. This quasar has the highest radio-to-optical flux density ratio at such a redshift, making it the radio-loudest source known to date at $z \sim 6$. The VLBA observations presented here resolve this quasar into multiple components with an overall linear extent of 1.62 kpc ($0\rlap{.}{''}28$) and with a total flux density of $6.57 \pm 0.38$ mJy, which is about half of the emission measured at a much lower angular resolution. The morphology of the source is comparable with either a radio core with a one-sided jet, or a compact or a medium-size Symmetric Object (CSO/MSO). If the source is a CSO/MSO, and assuming an advance speed of $0.2c$, then the estimated kinematic age is $\sim 10^4$ yr.

We present the discovery of the radio-loud quasar PSO J352.4034-15.3373 at z=5.84 pm 0.02. This quasar is the radio brightest source known, by an order of magnitude, at z~6 with a flux density in the range of 8-100 mJy from 3GHz to 230MHz and a radio loudness parameter R>~1000. This source provides an unprecedented opportunity to study powerful jets and radio-mode feedback at the highest redshifts, and presents the first real chance to probe deep into the neutral intergalactic medium by detecting 21 cm absorption at the end of cosmic reionization.

We present a survey of the [CII] 158 $\mu$m line and underlying far-infrared (FIR) dust continuum emission in a sample of 27 z>6 quasars using the Atacama Large Millimeter Array (ALMA) at ~1" resolution. The [CII] line was significantly detected (at >5-sigma) in 23 sources (85%). We find typical line luminosities of $L_{\rm [CII]}=10^{9-10}$ L$_\odot$, and an average line width of ~385 km/s. The [CII]-to-far-infrared luminosity ratio ([CII]/FIR) in our sources span one order of magnitude, highlighting a variety of conditions in the star-forming medium. Four quasar host galaxies are clearly resolved in their [CII] emission on a few kpc scales. Basic estimates of the dynamical masses of the host galaxies give masses between $2\times10^{10}$ and $2\times10^{11}$ M$_\odot$, i.e., more than an order of magnitude below what is expected from local scaling relations, given the available limits on the masses of the central black holes ($>3\times10^8$ M$_\odot$, assuming Eddington-limited accretion). In stacked ALMA [CII] spectra of individual sources in our sample, we find no evidence of a deviation from a single Gaussian profile. The quasar luminosity does not strongly correlate with either the [CII] luminosity or equivalent width. This survey (with typical on-source integration times of 8 min) showcases the unparalleled sensitivity of ALMA at millimeter wavelengths, and offers a unique reference sample for the study of the first massive galaxies in the universe.

On 2017 August 17, Swope Supernova Survey 2017a (SSS17a) was discovered as the optical counterpart of the binary neutron star gravitational wave event GW170817. We report time-series spectroscopy of SSS17a from 11.75 hours until 8.5 days after merger. Over the first hour of observations the ejecta rapidly expanded and cooled. Applying blackbody fits to the spectra, we measure the photosphere cooling from $11,000^{+3400}_{-900}$ K to $9300^{+300}_{-300}$ K, and determine a photospheric velocity of roughly 30% of the speed of light. The spectra of SSS17a begin displaying broad features after 1.46 days, and evolve qualitatively over each subsequent day, with distinct blue (early-time) and red (late-time) components. The late-time component is consistent with theoretical models of r-process-enriched neutron star ejecta, whereas the blue component requires high velocity, lanthanide-free material.

We present ALMA observations of the [CII] fine structure line and the underlying far-infrared (FIR) dust continuum emission in J1120+0641, the most distant quasar currently known (z=7.1). We also present observations targeting the CO(2-1), CO(7-6) and [CI] 369 micron lines in the same source obtained at the VLA and PdBI. We find a [CII] line flux of F_[CII]=1.11+/-0.10 Jy km/s and a continuum flux density of S_227GHz=0.53+/-0.04 mJy/beam, consistent with previous unresolved measurements. No other source is detected in continuum or [CII] emission in the field covered by ALMA (~25"). At the resolution of our ALMA observations (0.23", or 1.2 kpc, a factor ~70 smaller beam area compared to previous measurements), we find that the majority of the emission is very compact: a high fraction (~80%) of the total line and continuum flux is associated with a region 1-1.5 kpc in diameter. The remaining ~20% of the emission is distributed over a larger area with radius <4 kpc. The [CII] emission does not exhibit ordered motion on kpc-scales: applying the virial theorem yields an upper limit on the dynamical mass of the host galaxy of (4.3+/-0.9)x10^10 M_sun, only ~20x higher than the central black hole. The other targeted lines (CO(2-1), CO(7-6) and [CI]) are not detected, but the limits of the line ratios with respect to the [CII] emission imply that the heating in the quasar host is dominated by star formation, and not by the accreting black hole. The star-formation rate implied by the FIR continuum is 105-340 M_sun/yr, with a resulting star-formation rate surface density of ~100-350 M_sun/yr/kpc^2, well below the value for Eddington-accretion-limited star formation.

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

We present the selection algorithm and anticipated results for the Time Domain Spectroscopic Survey (TDSS). TDSS is an SDSS-IV eBOSS subproject that will provide initial identification spectra of approximately 220,000 luminosity-variable objects (variable stars and AGN) across 7,500 square degrees selected from a combination of SDSS and multi-epoch Pan-STARRS1 photometry. TDSS will be the largest spectroscopic survey to explicitly target variable objects, avoiding pre-selection on the basis of colors or detailed modeling of specific variability characteristics. Kernel Density Estimate (KDE) analysis of our target population performed on SDSS Stripe 82 data suggests our target sample will be 95% pure (meaning 95% of objects we select have genuine luminosity variability of a few magnitudes or more). Our final spectroscopic sample will contain roughly 135,000 quasars and 85,000 stellar variables, approximately 4,000 of which will be RR Lyrae stars which may be used as outer Milky Way probes. The variability-selected quasar population has a smoother redshift distribution than a color-selected sample, and variability measurements similar to those we develop here may be used to make more uniform quasar samples in large surveys. The stellar variable targets are distributed fairly uniformly across color space, indicating that TDSS will obtain spectra for a wide variety of stellar variables including pulsating variables, stars with significant chromospheric activity, cataclysmic variables and eclipsing binaries. TDSS will serve as a pathfinder mission to identify and characterize the multitude of variable objects that will be detected photometrically in even larger variability surveys such as LSST.

We present a deep and wide $I$ luminosity function for galaxies in Abell 1689 ($z=0.183$) from a mosaic of HST WFPC2 images covering $10'$ on the side. The main result of this work is the detection of a steep upturn in the dwarf galaxy LF, with $\alpha \sim -2$. The dwarf to giant ratio appears to increase outwards, but this is because giant galaxies are missing in the cluster outskirts, indicating luminosity segregation. The red sequence LF has the same parameters, within errors, as the total LF, showing that the faint end upturn consists of red quiescent galaxies. We speculate that the upturn is connected to the `filling-in' of the red sequence at $z < 0.4$ and may represent the latest installment of `downsizing' as the least massive galaxies are being quenched at the present epoch.

We measure alignments on scales of 1 Mpc $h^{-1}_{71}$ for galaxies in Abell 1689 ($z=0.18$) from an existing Hubble Space Telescope mosaic. We find evidence of galaxy alignment in the inner 500 $h^{-1}_{71}$ kpc. The alignment appears to be stronger towards the centre and is mostly present among the fainter galaxies, while bright galaxies are unaligned. This is consistent with a model where alignments originate from tidal locking.