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

Reconstruction techniques for intrinsic quasar continua are crucial for the precision study of Lyman-$\alpha$ (Ly-$\alpha$) and Lyman-$\beta$ (Ly-$\beta$) transmission at $z>5.0$, where the $\lambda<1215 A$ emission of quasars is nearly completely absorbed. While the number and quality of spectroscopic observations has become theoretically sufficient to quantify Ly-$\alpha$ transmission at $5.0<z<6.0$ to better than $1\%$, the biases and uncertainties arising from predicting the unabsorbed continuum are not known to the same level. In this paper, we systematically evaluate eight reconstruction techniques on a unified testing sample of $2.7<z<3.5$ quasars drawn from eBOSS. The methods include power-law extrapolation, stacking of neighbours, and six variants of Principal Component Analysis (PCA) using direct projection, fitting of components, or neural networks to perform weight mapping. We find that power-law reconstructions and the PCA with fewest components and smallest training sample display the largest biases in the Ly-$\alpha$ forest ($-9.58\%/+8.22\%$ respectively). Power-law extrapolations have larger scatters than previously assumed of $+13.1\%/-13.2\%$ over Ly-$\alpha$ and $+19.9\%/-20.1\%$ over Ly-$\beta$. We present two new PCAs which achieve the best current accuracies of $9\%$ for Ly-$\alpha$ and $17\%$ for Ly-$\beta$. We apply the eight techniques after accounting for wavelength-dependent biases and scatter to a sample $19$ quasars at $z>5.7$ with IR X-Shooter spectroscopy, obtaining well-characterised measurements for the mean flux transmission at $4.7<z<6.3$. Our results demonstrate the importance of testing and, when relevant, training, continuum reconstruction techniques in a systematic way.

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.

The mean Galactocentric radial velocities $\langle v_{R}\rangle(R,\phi)$ of luminous red giant stars within the mid-plane of the Milky Way reveal a spiral signature, which could plausibly reflect the response to a non-axisymmetric perturbation of the gravitational potential in the Galactic disk. We apply a simple steady-state toy model of a logarithmic spiral to interpret these observations, and find a good qualitative and quantitative match. Presuming that the amplitude of the gravitational potential perturbation is proportionate to that in the disk's surface mass density, we estimate the surface mass density amplitude to be $\Sigma_{\rm max} (R_{\odot})= 14\pm5\,\rm M_{\odot}\,pc^{-2}$ at the solar radius. Combined with the local disk density, this implies a surface mass density contrast between the arm and inter-arm regions of approximately $\pm 20\%$ at the solar radius, with an increase towards larger Galactocentric radii. Our model constrains the pitch angle of the {\it dynamical} spiral arms to be $12.7^{\circ}\pm1.2^{\circ}$.

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

The supermassive black holes (SMBHs) observed at the centers of all massive galaxies are believed to have grown via luminous accretion during quasar phases in the distant past. The fraction of inflowing rest mass energy emitted as light, the radiative efficiency, has been inferred to be 10\%, in agreement with expectations from thin disk accretion models. But the existence of billion solar-mass SMBHs powering quasars at $z > 7$ challenges this picture: provided they respect the Eddington limit, there is not enough time to grow $z>7$ SMBHs from stellar remnant seeds unless the radiative efficiency is below 10\%. Here we show that one can constrain the radiative efficiencies of the most distant quasars known using foreground neutral intergalactic gas as a cosmological-scale ionizing photon counter. From the Ly$\alpha$ absorption profiles of ULAS J1120+0641 ($z=7.09$) and ULAS J1342+0928 ($z=7.54$), we determine posterior median radiative efficiencies of 0.08\% and 0.1\%, respectively, and the combination of the two measurements rule out the canonical 10\% efficiency at 99.8\% credibility after marginalizing over the unknown obscured fraction. This low radiative efficiency implies rapid mass accretion for the earliest SMBHs, greatly easing the tension between the age of the Universe and the SMBH masses. However, our measured efficiency may instead reflect nearly complete obscuration by dusty gas in the quasar host galaxies over the vast majority of their SMBH growth. Assuming 10\% efficiency during unobscured phases, we find that the obscured fraction would be $>82\%$ at 95\% credibility, and imply a $25.7^{+49.6}_{-16.5}$ times larger obscured than unobscured luminous quasar population at $z>7$.

We measure the intergalactic medium (IGM) opacity in the Ly$\alpha$ as well as in the Ly$\beta$ forest along $19$ quasar sightlines between $5.5\lesssim z_{\rm abs}\lesssim 6.1$, probing the end stages of the reionization epoch. Owing to its lower oscillator strength the Ly$\beta$ transition is sensitive to different gas temperatures and densities than Ly$\alpha$, providing additional constraints on the ionization and thermal state of the IGM. A comparison of our measurements to different inhomogeneous reionization models, derived from post-processing the Nyx cosmological hydrodynamical simulation to include spatial fluctuations in the ultraviolet background (UVB) or the gas temperature field, as well as to a uniform reionization model with varying thermal states of the IGM, leads to two primary conclusions: First, we find that including the effects of spectral noise is key for a proper data to model comparison. Noise effectively reduces the sensitivity to high opacity regions, and thus even stronger spatial inhomogeneities are required to match the observed scatter in the observations than previously inferred. Second, we find that models which come close to reproducing the distribution of Ly$\alpha$ effective optical depths nevertheless underpredict the Ly$\beta$ opacity at the same spatial locations. The origin of this disagreement is not entirely clear but models with an inversion in the temperature-density relation of the IGM just after reionization is completed match our measurements best, although they still do not fully capture the observations at $z\gtrsim 5.8$.

Since the discovery of $z\sim 6$ quasars two decades ago, studies of their Ly$\alpha$-transparent proximity zones have largely focused on their utility as a probe of cosmic reionization. But even when in a highly ionized intergalactic medium, these zones provide a rich laboratory for determining the timescales that govern quasar activity and the concomitant growth of their supermassive black holes. In this work, we use a suite of 1D radiative transfer simulations of quasar proximity zones to explore their time-dependent behaviour for activity timescales from $\sim10^3$ to $10^8$ years. The sizes of the simulated proximity zones, as quantified by the distance at which the smoothed Ly$\alpha$ transmission drops below 10% (denoted $R_p$), are in excellent agreement with observations, with the exception of a handful of particularly small zones that have been attributed to extremely short $\lesssim 10^4$ year lifetimes. We develop a physically motivated semi-analytic model of proximity zones which captures the bulk of their equilibrium and non-equilibrium behaviour, and use this model to investigate how quasar variability on $\lesssim10^5$ year timescales is imprinted on the distribution of observed proximity zone sizes. We show that large variations in the ionizing luminosity of quasars on timescales of $\lesssim10^4$ years are disfavored based on the good agreement between the observed distribution of $R_p$ and our model prediction based on "lightbulb" (i.e. steady constant emission) light curves.

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 measure the circular velocity curve $v_{\rm c}(R)$ of the Milky Way with the highest precision to date across Galactocentric distances of $5\leq R \leq 25$ kpc. Our analysis draws on the $6$-dimensional phase-space coordinates of $\gtrsim 23,000$ luminous red-giant stars, for which we previously determined precise parallaxes using a data-driven model that combines spectral data from APOGEE with photometric information from WISE, 2MASS, and Gaia. We derive the circular velocity curve with the Jeans equation assuming an axisymmetric gravitational potential. At the location of the Sun we determine the circular velocity with its formal uncertainty to be $v_{\rm c}(R_{\odot}) = (229.0\pm0.2)\rm\,km\,s^{-1}$ with systematic uncertainties at the $\sim 2-5\%$ level. We find that the velocity curve is gently but significantly declining at $(-1.7\pm0.1)\rm\,km\,s^{-1}\,kpc^{-1}$, with a systematic uncertainty of $0.46\rm\,km\,s^{-1}\,kpc^{-1}$, beyond the inner $5$ kpc. We exclude the inner $5$ kpc from our analysis due to the presence of the Galactic bar, which strongly influences the kinematic structure and requires modeling in a non-axisymmetric potential. Combining our results with external measurements of the mass distribution for the baryonic components of the Milky Way from other studies, we estimate the Galaxy's dark halo mass within the virial radius to be $M_{\rm vir} = (7.25\pm0.26)\cdot 10^{11}M_{\odot}$ and a local dark matter density of $\rho_{\rm dm}(R_{\odot}) = 0.30\pm0.03\,\rm GeV\,cm^{-3}$.

With contemporary infrared spectroscopic surveys like APOGEE, red-giant stars can be observed to distances and extinctions at which Gaia parallaxes are not highly informative. Yet the combination of effective temperature, surface gravity, composition, and age - all accessible through spectroscopy - determines a giant's luminosity. Therefore spectroscopy plus photometry should enable precise spectrophotometric distance estimates. Here we use the APOGEE-Gaia-2MASS-WISE overlap to train a data-driven model to predict parallaxes for red-giant branch stars with $0<\log g\leq2.2$ (more luminous than the red clump). We employ (the exponentiation of) a linear function of APOGEE spectral pixel intensities and multi-band photometry to predict parallax spectrophotometrically. The model training involves no logarithms or inverses of the Gaia parallaxes, and needs no cut on the Gaia parallax signal-to-noise ratio. It includes an L1 regularization to zero out the contributions of uninformative pixels. The training is performed with leave-out subsamples such that no star's astrometry is used even indirectly in its spectrophotometric parallax estimate. The model implicitly performs a reddening and extinction correction in its parallax prediction, without any explicit dust model. We assign to each star in the sample a new spectrophotometric parallax estimate; these parallaxes have uncertainties of a few to 15 percent, depending on data quality, which is more precise than the Gaia parallax for the vast majority of targets, and certainly any stars more than a few kpc distance. We obtain 10-percent distance estimates out to heliocentric distances of $20\,$kpc, and make global maps of the Milky Way's disk.

We publicly release a new sample of $34$ medium resolution quasar spectra at $5.77\leq z_{\rm em}\leq6.54$ observed with the Echellette Spectrograph and Imager (ESI) on the Keck telescope. This quasar sample represents an ideal laboratory to study the intergalactic medium (IGM) during the end stages of the epoch of reionization, and constrain the timing and morphology of the phase transition. For a subset of $23$ of our highest signal-to-noise ratio spectra (S/N$>7$, per $10\,{\rm km\,s^{-1}}$ pixel), we present a new measurement of the Lyman-$\alpha$ (Ly$\alpha$) forest opacity spanning the redshift range $4.8\lesssim z\lesssim6.3$. We carefully eliminate spectral regions that could be causing biases in our measurements due to additional transmitted flux in the proximity zone of the quasars, or extra absorption caused by strong intervening absorption systems along the line of sight. We compare the observed evolution of the IGM opacity with redshift to predictions from a hydrodynamical simulation with uniform ultraviolet background (UVB) radiation, as well as two semi-numerical patchy reionization models, one with a fluctuating UVB and another with a fluctuating temperature field. Our measurements show a steep rise in opacity at $z\gtrsim5.0$ and an increased scatter and thus support the picture of a spatially inhomogeneous reionization process, consistent with previous work. However, we measure significantly higher optical depths at $5.3\lesssim z\lesssim5.7$ than previous studies, which reduces the contrast between the highest opacity Gunn-Peterson troughs and the average opacity trend of the IGM, which may relieve some of the previously noted tension between these measurements and reionization models.

The quasar lifetime $t_{\rm\,Q}$ is one of the most fundamental quantities for understanding quasar evolution and the growth of supermassive black holes (SMBHs), but remains uncertain by several orders of magnitude. In a recent study we uncovered a population of very young quasars ($t_{\rm Q}\lesssim10^4-10^5$ yr), based on the sizes of their proximity zones, which are regions of enhanced Ly$\alpha$ forest transmission near the quasar resulting from its own ionizing radiation. The presence of such young objects poses significant challenges to models of SMBH formation, which already struggle to explain the existence of SMBHs at such early cosmic epochs. We conduct the first comprehensive spectroscopic study of the youngest quasar known, $\rm SDSS\,J1335+3533$ at $z=5.9012$, whose lifetime is $t_{\rm Q}<10^4$ yr ($95\%$ confidence). A careful search of our deep optical and near-infrared spectra for HI and metal absorption lines allows us to convincingly exclude that its small proximity zone results from an associated absorption system rather than a short lifetime. We use the MgII emission line to measure its black hole mass $M_{\rm BH}=(4.09\pm0.58)\times 10^9 M_{\odot}$, implying an Eddington ratio of $0.30\pm0.04$ -- comparable to other co-eval quasars. We similarly find that the relationship between its black hole mass and dynamical mass are consistent with other $z\sim6$ quasars. The only possible anomaly associated with youth are its weak emission lines, but larger samples are needed to shed light on a potential causal connection. We discuss the implications of short lifetimes for various SMBH growth scenarios, and argue that future observations of young quasars with JWST could distinguish between them.

We introduce a new method for determining the influence of galaxies and active galactic nuclei (AGN) on the physical state of the intergalactic medium (IGM) at high redshift and illustrate its potential via a first application to the field of the $z=6.42$ QSO J1148+5251. By correlating the spatial positions of spectroscopically-confirmed Lyman break galaxies (LBGs) with fluctuations in the Lyman alpha forest seen in the high signal-to-noise spectrum of a background QSO, we provide a statistical measure of the typical escape fraction of Lyman continuum photons close to the end of cosmic reionisation. Here we use Keck DEIMOS spectroscopy to locate 7 colour-selected LBGs in the redshift range $5.3\lesssim z\lesssim 6.4$ and confirm a faint $z=5.701$ AGN. We then examine the spatial correlation between this sample and Ly$\alpha$/Ly$\beta$ transmission fluctuations in a Keck ESI spectrum of the QSO. Interpreting the statistical HI proximity effect as arising from faint galaxies clustered around the detected LBGs, we translate the observed mean Ly$\alpha$ transmitted flux around an average detected LBG into a constraint on the mean escape fraction $\langle f_{\rm esc}\rangle\geq0.08$ at $z\simeq6$. We also report evidence of the individual transverse HI proximity effect of a $z=6.177$ luminous LBG via a Ly$\beta$ transmission spike and two broad Ly$\alpha$ transmission spikes around the $z=5.701$ AGN. We discuss the possible origin of such associations which suggest that while faint galaxies are primarily driving reionisation, luminous galaxies and AGN may provide important contributions to the UV background or thermal fluctuations of the IGM at $z\simeq6$. Although a limited sample, our results demonstrate the potential of making progress using this method in resolving one of the most challenging aspects of the contribution of galaxies and AGN to cosmic reionisation.

The observable spectrum of an unresolved binary star system is a superposition of two single-star spectra. Even without a detectable velocity offset between the two stellar components, the combined spectrum of a binary system is in general different from that of either component, and fitting it with single-star models may yield inaccurate stellar parameters and abundances. We perform simple experiments with synthetic spectra to investigate the effect of unresolved main-sequence binaries on spectral fitting, modeling spectra similar to those collected by the APOGEE, GALAH, and LAMOST surveys. We find that fitting unresolved binaries with single-star models introduces systematic biases in the derived stellar parameters and abundances that are modest but certainly not negligible, with typical systematic errors of $300\,\rm K$ in $T_{\rm eff}$, 0.1 dex in $\log g$, and 0.1 dex in $[\rm Fe/H]$ for APOGEE-like spectra of solar-type stars. These biases are smaller for spectra at optical wavelengths than in the near-infrared. We show that biases can be corrected by fitting spectra with a binary model, which adds only two labels to the fit and includes single-star models as a special case. Our model provides a promising new method to constrain the Galactic binary population, including systems with single-epoch spectra and no detectable velocity offset between the two stars.

Quasars are galaxies hosting accreting supermassive black holes; due to their brightness, they are unique probes of the early universe. To date, only few quasars have been reported at $z > 6.5$ ($<$800 Myr after the Big Bang). In this work, we present six additional $z \gtrsim 6.5$ quasars discovered using the Pan-STARRS1 survey. We use a sample of 15 $z \gtrsim 6.5$ quasars to perform a homogeneous and comprehensive analysis of this highest-redshift quasar population. We report four main results: (1) the majority of $z\gtrsim$6.5 quasars show large blueshifts of the broad CIV 1549\AA$\,$emission line compared to the systemic redshift of the quasars, with a median value $\sim$3$\times$ higher than a quasar sample at $z\sim$1; (2) we estimate the quasars' black hole masses (M$\rm_{BH}\sim$0.3$-$5 $\times$ 10$^{9}$ M$_{\odot}$) via modeling of the MgII 2798\AA$\,$emission line and rest-frame UV continuum; we find that quasars at high redshift accrete their material (with $\langle (L_{\mathrm{bol}}/L_{\mathrm{Edd}}) \rangle = 0.39$) at a rate comparable to a luminosity-matched sample at lower$-$redshift, albeit with significant scatter ($0.4$ dex); (3) we recover no evolution of the FeII/MgII abundance ratio with cosmic time; (4) we derive near zone sizes; together with measurements for $z\sim6$ quasars from recent work, we confirm a shallow evolution of the decreasing quasar near zone sizes with redshift. Finally, we present new millimeter observations of the [CII] 158 $\mu$m emission line and underlying dust continuum from NOEMA for four quasars, and provide new accurate redshifts and [CII]/infrared luminosities estimates. The analysis presented here shows the large range of properties of the most distant quasars.

We present a new Bayesian algorithm making use of Markov Chain Monte Carlo sampling that allows us to simultaneously estimate the unknown continuum level of each quasar in an ensemble of high-resolution spectra, as well as their common probability distribution function (PDF) for the transmitted Ly$\alpha$ forest flux. This fully automated PDF regulated continuum fitting method models the unknown quasar continuum with a linear Principal Component Analysis (PCA) basis, with the PCA coefficients treated as nuisance parameters. The method allows one to estimate parameters governing the thermal state of the intergalactic medium (IGM), such as the slope of the temperature-density relation $\gamma-1$, while marginalizing out continuum uncertainties in a fully Bayesian way. Using realistic mock quasar spectra created from a simplified semi-numerical model of the IGM, we show that this method recovers the underlying quasar continua to a precision of $\simeq7\%$ and $\simeq10\%$ at $z=3$ and $z=5$, respectively. Given the number of principal component spectra, this is comparable to the underlying accuracy of the PCA model itself. Most importantly, we show that we can achieve a nearly unbiased estimate of the slope $\gamma-1$ of the IGM temperature-density relation with a precision of $\pm8.6\%$ at $z=3$, $\pm6.1\%$ at $z=5$, for an ensemble of ten mock high-resolution quasar spectra. Applying this method to real quasar spectra and comparing to a more realistic IGM model from hydrodynamical simulations would enable precise measurements of the thermal and cosmological parameters governing the IGM, albeit with somewhat larger uncertainties given the increased flexibility of the model.

We study quasar proximity zones in the redshift range $5.77 \leq z \leq 6.54$ by homogeneously analyzing $34$ medium resolution spectra, encompassing both archival and newly obtained data, and exploiting recently updated systemic redshift and magnitude measurements. Whereas previous studies found strong evolution of proximity zone sizes with redshift, and argued that this provides evidence for a rapidly evolving intergalactic medium (IGM) neutral fraction during reionization, we measure a much shallower trend $\propto(1+z)^{-1.44}$. We compare our measured proximity zone sizes to predictions from hydrodynamical simulations post-processed with one-dimensional radiative transfer, and find good agreement between observations and theory irrespective of the ionization state of the ambient IGM. This insensitivity to IGM ionization state has been previously noted, and results from the fact that the definition of proximity zone size as the first drop of the smoothed quasar spectrum below the $10\%$ flux transmission level probes locations where the ionizing radiation from the quasar is an order of magnitude larger than the expected ultraviolet ionizing background that sets the neutral fraction of the IGM. Our analysis also uncovered three objects with exceptionally small proximity zones (two have $R_p < 1$proper Mpc), which constitute outliers from the observed distribution and are challenging to explain with our radiative transfer simulations. We consider various explanations for their origin, such as strong absorption line systems associated with the quasar or patchy reionization, but find that the most compelling scenario is that these quasars have been shining for $\lesssim 10^5$yr.

The amplitude of the ionizing background that pervades the intergalactic medium (IGM) at the end of the epoch of reionization provides a valuable constraint on the emissivity of the sources which reionized the Universe. While measurements of the ionizing background at lower redshifts rely on a simulation-calibrated mapping between the photoionization rate and the mean transmission of the Ly$\alpha$ forest, at $z\gtrsim6$ the IGM becomes increasingly opaque, and transmission arises solely in narrow spikes separated by saturated Gunn-Peterson troughs. In this regime, the traditional approach of measuring the average transmission over large $\sim 50$ Mpc$/h$ regions is less sensitive and sub-optimal. Additionally, the five times smaller oscillator strength of the Ly$\beta$ transition implies the Ly$\beta$ forest is considerably more transparent at $z\gtrsim6$, even in the presence of contamination by foreground $z\sim 5$ Ly$\alpha$ forest absorption. In this work we present a novel statistical approach to analyze the joint distribution of transmission spikes in the co-spatial $z\sim 6$ Ly$\alpha$ and Ly$\beta$ forests. Our method relies on Approximate Bayesian Computation (ABC), which circumvents the necessity of computing the intractable likelihood function describing the highly correlated Ly$\alpha$ and Ly$\beta$ transmission. We apply ABC to mock data generated from a large-volume hydrodynamical simulation combined with a state-of-the-art model of ionizing background fluctuations in the post-reionization IGM, and show that it is sensitive to higher IGM neutral hydrogen fractions than previous techniques. As a proof of concept, we apply this methodology to a real spectrum of a $z=6.54$ quasar and measure the ionizing background from $5.4\leq z \leq 6.4$ along this sightline with $\sim0.2$ dex statistical uncertainties.

We present the first observations of foreground Lyman-$\alpha$ forest absorption from high-redshift galaxies, targeting 24 star-forming galaxies (SFGs) with $z\sim 2.3-2.8$ within a $5' \times 15'$ region of the COSMOS field. The transverse sightline separation is $\sim 2\,h^{-1}\mathrm{Mpc}$ comoving, allowing us to create a tomographic reconstruction of the 3D Ly$\alpha$ forest absorption field over the redshift range $2.20\leq z\leq 2.45$. The resulting map covers $6\,h^{-1}\mathrm{Mpc} \times 14\,h^{-1}\mathrm{Mpc}$ in the transverse plane and $230\,h^{-1}\mathrm{Mpc}$ along the line-of-sight with a spatial resolution of $\approx 3.5\,h^{-1}\mathrm{Mpc}$, and is the first high-fidelity map of large-scale structure on $\sim\mathrm{Mpc}$ scales at $z>2$. Our map reveals significant structures with $\gtrsim 10\,h^{-1}\mathrm{Mpc}$ extent, including several spanning the entire transverse breadth, providing qualitative evidence for the filamentary structures predicted to exist in the high-redshift cosmic web. Simulated reconstructions with the same sightline sampling, spectral resolution, and signal-to-noise ratio recover the salient structures present in the underlying 3D absorption fields. Using data from other surveys, we identified 18 galaxies with known redshifts coeval with our map volume enabling a direct comparison to our tomographic map. This shows that galaxies preferentially occupy high-density regions, in qualitative agreement with the same comparison applied to simulations. Our results establishes the feasibility of the CLAMATO survey, which aims to obtain Ly$\alpha$ forest spectra for $\sim 1000$ SFGs over $\sim 1 \,\mathrm{deg}^2$ of the COSMOS field, in order to map out IGM large-scale structure at $\langle z \rangle \sim 2.3$ over a large volume $(100\,h^{-1}\mathrm{Mpc})^3$.

(abridged) The Hundred-Thousand-Proper-Motion (HTPM) project will determine the proper motions of ~113500 stars using a 23-year baseline. The proper motions will use the Hipparcos data, with epoch 1991.25, as first epoch and the first intermediate-release Gaia astrometry, with epoch ~2014.5, as second epoch. The expected HTPM proper-motion standard errors are 30-190 muas/yr, depending on stellar magnitude. Depending on the characteristics of an object, in particular its distance and velocity, its radial velocity can have a significant impact on the determination of its proper motion. The impact of this perspective acceleration is largest for fast-moving, nearby stars. Our goal is to determine, for each star in the Hipparcos catalogue, the radial-velocity standard error that is required to guarantee a negligible contribution of perspective acceleration to the HTPM proper-motion precision. We employ two evaluation criteria, both based on Monte-Carlo simulations, with which we determine which stars need to be spectroscopically (re-)measured. Both criteria take the Hipparcos measurement errors into account. For each star in the Hipparcos catalogue, we determine the confidence level with which the available radial velocity and its standard error, taken from the XHIP compilation catalogue, are acceptable. We find that for 97 stars, the radial velocities available in the literature are insufficiently precise for a 68.27% confidence level. We also identify 109 stars for which radial velocities are currently unknown yet need to be acquired to meet the 68.27% confidence level. To satisfy the radial-velocity requirements coming from our study will be a daunting task consuming a significant amount of spectroscopic telescope time. Fortunately, the follow-up spectroscopy is not time-critical since the HTPM proper motions can be corrected a posteriori once (improved) radial velocities become available.