Researcher profile

Frederick B. Davies

Frederick B. Davies contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
20works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2024arXiv

Recognizing Blazars Using Radio Morphology from the VLA Sky Survey

Blazars are radio-loud Active Galactic Nuclei (AGN) whose jets have a very small angle to our line of sight. Observationally, the radio emission are mostly compact or a compact-core with a 1-sided jet. With 2.5$^{\prime\prime}$ resolution at 3 GHz, the Very Large Array Sky Survey (VLASS) enables us to resolve the structure of some blazar candidates in the sky north of Decl. $-40$ deg. We introduce an algorithm to classify radio sources as either blazar-like or non-blazar-like based on their morphology in the VLASS images. We apply our algorithm to three existing catalogs, including one of known blazars (Roma-BzCAT) and two of blazar candidates identified by WISE colors and radio emission (WIBRaLS, KDEBLLACS). We show that in all three catalogs, there are objects with morphology inconsistent with being blazars. Considering all the catalogs, more than 12% of the candidates are unlikely to be blazars, based on this analysis. Notably, we show that 3% of the Roma-BzCAT "confirmed" blazars could be a misclassification based on their VLASS morphology. The resulting table with all sources and their radio morphological classification is available online.

preprint2022arXiv

(Nearly) Model-Independent Constraints on the Neutral Hydrogen Fraction in the Intergalactic Medium at $z\sim 5-7$ Using Dark Pixel Fractions in Ly$α$ and Ly$β$ Forests

Cosmic reionization was the last major phase transition of hydrogen from neutral to highly ionized in the intergalactic medium (IGM). Current observations show that the IGM is significantly neutral at $z>7$, and largely ionized by $z\sim5.5$. However, most methods to measure the IGM neutral fraction are highly model-dependent, and are limited to when the volume-averaged neutral fraction of the IGM is either relatively low ($\bar{x}_{\rm HI} \lesssim 10^{-3}$) or close to unity ($\bar{x}_{\rm HI}\sim 1$). In particular, the neutral fraction evolution of the IGM at the critical redshift range of $z=6-7$ is poorly constrained. We present new constraints on $\bar{x}_{\rm HI}$ at $z\sim5.1-6.8$, by analyzing deep optical spectra of $53$ quasars at $5.73<z<7.09$. We derive model-independent upper limits on the neutral hydrogen fraction based on the fraction of &#34;dark&#34; pixels identified in the Lyman $α$ (Ly$α$) and Lyman $β$ (Ly$β$) forests, without any assumptions on the IGM model or the intrinsic shape of the quasar continuum. They are the first model-independent constraints on the IGM neutral hydrogen fraction at $z\sim6.2-6.8$ using quasar absorption measurements. Our results give upper limits of $\bar{x}_{\rm HI}(z=6.3) < 0.79\pm0.04$ (1$σ$), $\bar{x}_{\rm HI} (z=6.5) < 0.87\pm0.03$ (1$σ$), and $\bar{x}_{\rm HI} (z=6.7) < 0.94^{+0.06}_{-0.09}$ (1$σ$). The dark pixel fractions at $z>6.1$ are consistent with the redshift evolution of the neutral fraction of the IGM derived from the Planck 2018.

preprint2022arXiv

HeII Ly$α$ Transmission Spikes and Absorption Troughs in Eight High-resolution Spectra Probing the End of HeII Reionization

We present statistics of HeII Lya transmission spikes and large-scale absorption troughs using archival high-resolution ($R=λ/Δλ\simeq 12,500$-$18,000$) far-UV spectra of eight HeII-transparent quasars obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope. The sample covers the redshift range 2.5<z<3.8, thereby probing the rapidly evolving HeII absorption at the end of the HeII reionization epoch. The measured lengths of the troughs decrease dramatically from L>100cMpc at z>3 to L~30cMpc at z~2.7, signaling a significant progression of HeII reionization at these redshifts. Furthermore, unexpectedly long L~65cMpc troughs detected at z~2.9 suggest that the UV background fluctuates at larger scales than predicted by current models. By comparing the measured incidence of transmission spikes to predictions from forward-modeled mock spectra created from the outputs of a (146cMpc)^3 optically thin Nyx hydrodynamical simulation employing different UV background models, we infer the redshift evolution of the HeII photoionization rate $Γ_\mathrm{He\,II}(z)$. The photoionization rate decreases with increasing redshift from $\simeq 4.6\times 10^{-15}\mathrm{\,s^{-1}}$ at z~2.6 to $\simeq 1.2 \times 10^{-15}\mathrm{\,s^{-1}}$ at z~3.2, in agreement with previous inferences from the HeII effective optical depth, and following expected trends of current models of a fluctuating HeII-ionizing background.

preprint2022arXiv

Hydrogen reionisation ends by $z=5.3$: Lyman-$α$ optical depth measured by the XQR-30 sample

The presence of excess scatter in the Ly-$α$ forest at $z\sim 5.5$, together with the existence of sporadic extended opaque Gunn-Peterson troughs, has started to provide robust evidence for a late end of hydrogen reionisation. However, low data quality and systematic uncertainties complicate the use of Ly-$α$ transmission as a precision probe of reionisation&#39;s end stages. In this paper, we assemble a sample of 67 quasar sightlines at $z>5.5$ with high signal-to-noise ratios of $>10$ per $\leq 15$ km s$^{-1}$ spectral pixel, relying largely on the new XQR-30 quasar sample. XQR-30 is a large program on VLT/X-Shooter which obtained deep (SNR $>20$ per pixel) spectra of 30 quasars at $z>5.7$. We carefully account for systematics in continuum reconstruction, instrumentation, and contamination by damped Ly-$α$ systems. We present improved measurements of the mean Ly-$α$ transmission over $4.9<z<6.1$. Using all known systematics in a forward modelling analysis, we find excellent agreement between the observed Ly-$α$ transmission distributions and the homogeneous-UVB simulations Sherwood and Nyx up to $z\leq5.2$ ($<1 σ$), and mild tension ($\sim 2.5 σ$) at $z=5.3$. Homogeneous UVB models are ruled out by excess Ly-$α$ transmission scatter at $z\geq5.4$ with high confidence ($>3.5 σ$). Our results indicate that reionisation-related fluctuations, whether in the UVB, residual neutral hydrogen fraction, and/or IGM temperature, persist in the intergalactic medium until at least $z=5.3$ ($t=1.1$ Gyr after the Big Bang). This is further evidence for a late end to reionisation.

preprint2022arXiv

Improved treatments of the ionizing photon mean free path in semi-numerical simulations of reionization

Efficient and accurate simulations of the reionization epoch are crucial to exploring the vast uncharted parameter space that will soon be constrained by measurements of the 21 cm power spectrum. One of these parameters, $R_{\rm max}$, is meant to characterize the absorption of photons by residual neutral gas inside of ionized regions, but has historically been implemented in a very simplistic fashion acting only as a maximum filtering scale. We leverage the correspondence between excursion set methods and the integrated flux from ionizing sources to define two physically-motivated prescriptions of the mean free path of ionizing photons that smoothly attenuate the contribution from distant sources. Implementation of these methods in semi-numerical reionization codes requires only modest additional computational effort due to the fact that spatial filtering is still performed on scales larger than the characteristic absorption distance. We find that our smoothly-defined mean free path prescriptions more effectively suppress large-scale structures in the ionization field in semi-numerical reionization simulations compared to the standard $R_{\rm max}$ approach, and the magnitude of the mean free path modulates the power spectrum in a much smoother manner. We show that this suppression of large-scale power is significant enough to be relevant for upcoming 21 cm power spectrum observations. Finally, we show that in our model the mean free path plays a larger role in regulating the reionization history than in models using $R_{\rm max}$.

preprint2022arXiv

Long Dark Gaps in the Ly$β$ Forest at $z<6$: Evidence of Ultra Late Reionization from XQR-30 Spectra

We present a new investigation of the intergalactic medium (IGM) near reionization using dark gaps in the Lyman-$β$ (Ly$β$) forest. With its lower optical depth, Ly$β$ offers a potentially more sensitive probe to any remaining neutral gas compared to commonly used Ly$α$ line. We identify dark gaps in the Ly$β$ forest using spectra of 42 QSOs at $z_{\rm em}>5.5$, including new data from the XQR-30 VLT Large Programme. Approximately $40\%$ of these QSO spectra exhibit dark gaps longer than $10h^{-1}{\rm Mpc}$ at $z\simeq5.8$. By comparing the results to predictions from simulations, we find that the data are broadly consistent both with models where fluctuations in the Ly$α$ forest are caused solely by ionizing ultraviolet background (UVB) fluctuations and with models that include large neutral hydrogen patches at $z<6$ due to a late end to reionization. Of particular interest is a very long ($L=28h^{-1}{\rm Mpc}$) and dark ($τ_{\rm eff} \gtrsim 6$) gap persisting down to $z\simeq 5.5$ in the Ly$β$ forest of the $z_{\rm}=5.85$ QSO PSO J025$-$11. This gap may support late reionization models with a volume-weighted average neutral hydrogen fraction of $ \langle x_{\rm HI}\rangle \gtrsim 5\%$ by $z=5.6$. Finally, we infer constraints on $\langle x_{\rm HI}\rangle$ over $5.5 \lesssim z \lesssim 6.0$ based on the observed Ly$β$ dark gap length distribution and a conservative relationship between gap length and neutral fraction derived from simulations. We find $\langle x_{\rm HI}\rangle \leq 0.05$, 0.17, and 0.29 at $z\simeq 5.55$, 5.75, and 5.95, respectively. These constraints are consistent with models where reionization ends significantly later than $z = 6$.

preprint2022arXiv

Staring at the Shadows of Archaic Galaxies: Damped Ly$α$ and Metal Absorbers toward a Young $z \sim 6$ Weak-line Quasar

We characterize the Ly$α$ halo and absorption systems toward PSO J083+11, a unique $z=6.3401$ weak-line quasar, using Gemini/GNIRS, Magellan/FIRE, and VLT/MUSE data. Strong absorptions by hydrogen and several metal lines (e.g., CII, MgII, and OI) are discovered in the spectrum, which indicates the presence of: (i) a proximate sub-damped Ly$α$ (sub-DLA) system at $z=6.314$ and (ii) a MgII absorber at $z=2.2305$. To describe the observed damping wing signal, we model the Ly$α$ absorption with a combination of a sub-DLA with the neutral hydrogen column density of $\log N_\mathrm{HI} = 20.03 \pm 0.30$ cm$^{-2}$ and absorption from the intergalactic medium with a neutral fraction of around 10 percent. The sub-DLA toward PSO J083+11 has an abundance ratio of [C/O] $=-0.04 \pm 0.33$ and metallicity of [O/H] $=-2.19 \pm 0.44$, similar to those of low-redshift metal-poor DLAs. These measurements suggest that the sub-DLA might truncate PSO J083+11&#39;s proximity zone size and complicate the quasar lifetime measurement. However, this quasar shows no sign of a Ly$α$ halo in the MUSE data cube, where the estimated $1σ$ limit of surface brightness is $2.76 \times 10^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ at aperture size of 1 arcsecond, or equivalent to a Ly$α$ luminosity of $\leq 43.46$ erg s$^{-1}$. This non-detection, while being only weak independent evidence on its own, is at least consistent with a young quasar scenario, as expected for a quasar with a short accretion timescale.

preprint2021arXiv

A Luminous Quasar at Redshift 7.642

Distant quasars are unique tracers to study the formation of the earliest supermassive black holes (SMBHs) and the history of cosmic reionization. Despite extensive efforts, only two quasars have been found at $z\ge7.5$, due to a combination of their low spatial density and the high contamination rate in quasar selection. We report the discovery of a luminous quasar at $z=7.642$, J0313$-$1806, the most distant quasar yet known. This quasar has a bolometric luminosity of $3.6\times10^{13} L_\odot$. Deep spectroscopic observations reveal a SMBH with a mass of $(1.6\pm0.4) \times10^9M_\odot$ in this quasar. The existence of such a massive SMBH just $\sim$670 million years after the Big Bang challenges significantly theoretical models of SMBH growth. In addition, the quasar spectrum exhibits strong broad absorption line (BAL) features in CIV and SiIV, with a maximum velocity close to 20% of the speed of light. The relativistic BAL features, combined with a strongly blueshifted CIV emission line, indicate that there is a strong active galactic nucleus (AGN) driven outflow in this system. ALMA observations detect the dust continuum and [CII] emission from the quasar host galaxy, yielding an accurate redshift of $7.6423 \pm 0.0013$ and suggesting that the quasar is hosted by an intensely star-forming galaxy, with a star formation rate of $\rm\sim 200 ~M_\odot ~yr^{-1}$ and a dust mass of $\sim7\times10^7~M_\odot$. Followup observations of this reionization-era BAL quasar will provide a powerful probe of the effects of AGN feedback on the growth of the earliest massive galaxies.

preprint2021arXiv

IGM damping wing constraints on reionisation from covariance reconstruction of two $z\gtrsim7$ QSOs

Bright, high redshift ($z>6$) QSOs are powerful probes of the ionisation state of the intervening intergalactic medium (IGM). The detection of Ly$α$ damping wing absorption imprinted in the spectrum of high-z QSOs can provide strong constraints on the epoch of reionisation (EoR). In this work, we perform an independent Ly$α$ damping wing analysis of two known $z>7$ QSOs; DESJ0252-0503 at $z=7.00$ (Wang et al.) and J1007+2115 at $z=7.51$ (Yang et al.). For this, we utilise our existing Bayesian framework which simultaneously accounts for uncertainties in: (i) the intrinsic Ly$α$ emission profile (reconstructed from a covariance matrix of measured emission lines; extended in this work to include NV) and (ii) the distribution of ionised (H\,{\scriptsize II}) regions within the IGM using a $1.6^3$ Gpc$^3$ reionisation simulation. This approach is complementary to that used in the aforementioned works as it focuses solely redward of Ly$α$ ($1218 < λ< 1230$Å) making it more robust to modelling uncertainties while also using a different methodology for (i) and (ii). We find, for a fiducial EoR morphology, $\bar{x}_{\rm HI} = 0.64\substack{+0.19 \\ -0.23}$ (68 per cent) at $z=7$ and $\bar{x}_{\rm HI} = 0.27\substack{+0.21 \\ -0.17}$ at $z=7.51$ consistent within $1σ$ to the previous works above, though both are slightly lower in amplitude. Following the inclusion of NV into our reconstruction pipeline, we perform a reanalysis of ULASJ1120+0641 at $z=7.09$ (Mortlock et al.) and ULASJ1342+0928 at $z=7.54$ (Bañados et al.) finding $\bar{x}_{\rm HI} = 0.44\substack{+0.23 \\ -0.24}$ at $z=7.09$ and $\bar{x}_{\rm HI} = 0.31\substack{+0.18 \\ -0.19}$ at $z=7.54$. Finally, we combine the QSO damping wing constraints for all four $z\gtrsim7$ QSOs to obtain a single, unified constraint of $\bar{x}_{\rm HI} = 0.49\substack{+0.11 \\ -0.11}$ at $z=7.29$.

preprint2021arXiv

Probing Early Super-massive Black Hole Growth and Quasar Evolution with Near-infrared Spectroscopy of 37 Reionization-era Quasars at 6.3 < z <= 7.64

We report the results of near-infrared spectroscopic observations of 37 quasars in the redshift range $6.3< z\le7.64$, including 32 quasars at $z>6.5$, forming the largest quasar near-infrared spectral sample at this redshift. The spectra, taken with Keck, Gemini, VLT, and Magellan, allow investigations of central black hole mass and quasar rest-frame ultraviolet spectral properties. The black hole masses derived from the MgII emission lines are in the range $(0.3-3.6)\times10^{9}\,M_{\odot}$, which requires massive seed black holes with masses $\gtrsim10^{3-4}\,M_{\odot}$, assuming Eddington accretion since $z=30$. The Eddington ratio distribution peaks at $λ_{\rm Edd}\sim0.8$ and has a mean of 1.08, suggesting high accretion rates for these quasars. The CIV - MgII emission line velocity differences in our sample show an increase of CIV blueshift towards higher redshift, but the evolutionary trend observed from this sample is weaker than the previous results from smaller samples at similar redshift. The FeII/MgII flux ratios derived for these quasars up to $z=7.6$, compared with previous measurements at different redshifts, do not show any evidence of strong redshift evolution, suggesting metal-enriched environments in these quasars. Using this quasar sample, we create a quasar composite spectrum for $z>6.5$ quasars and find no significant redshift evolution of quasar broad emission lines and continuum slope, except for a blueshift of the CIV line. Our sample yields a strong broad absorption line quasar fraction of $\sim$24%, higher than the fractions in lower redshift quasar samples, although this could be affected by small sample statistics and selection effects.

preprint2020arXiv

A metal-poor damped Ly-alpha system at redshift 6.4

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.

preprint2020arXiv

A Significantly Neutral Intergalactic Medium Around the Luminous z=7 Quasar J0252-0503

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$α$ 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 $λ_{\rm Edd}=0.7\pm0.1$, consistent with the values seen in other luminous $z\sim 7$ quasars. Furthermore, the Ly$α$ 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.

preprint2020arXiv

Detecting and Characterizing Young Quasars I: Systemic Redshifts and Proximity Zones Measurements

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&#39; 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&#39; 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}$.

preprint2020arXiv

No Redshift Evolution in the Broad Line Region Metallicity up to z=7.54: Deep NIR Spectroscopy of ULAS J1342+0928

We present deep (9 hours) Gemini-N/GNIRS near-infrared spectroscopic observations of ULAS J1342+0928, a luminous quasar at z=7.54. Various broad emission lines were detected, as well as the underlying continuum and iron forests over the rest-frame wavelength 970-2930A. There is a clear trend that higher-ionization emission lines show larger blueshifts with CIV1549 exhibiting 5510^{+240}_{-110} km s-1 blueshift with respect to the systematic redshift from the far-infrared [CII] 158um emission line. Those high ionization lines have wide profiles with FWHM more than 10000 km s-1. A modest blueshift of 340^{+110}_{-80} km s-1 is also seen in MgII, the lowest ionization line identified in the spectrum. The updated MgII-based black hole mass of M_BH=9.1_{-1.3}^{+1.4} x 10^8 M_sun and the Eddington ratio of L_bol/L_Edd=1.1_{-0.2}^{+0.2} confirm that ULAS J1342+0928 is powered by a massive and actively accreting black hole. There is no significant difference in the emission line ratios such as SiIV/CIV and AlIII/CIV when compared to lower-redshift quasars in a similar luminosity range, which suggests early metal pollution of the broad-line-region clouds. This trend also holds for the FeII/MgII line ratio, known as a cosmic clock that traces the iron enrichment in the early universe. Different iron templates and continuum fitting ranges were used to explore how the FeII/MgII measurement changes as a function of spectral modeling. Quasars at even higher redshift or at fainter luminosity range (L_bol<10^46 erg s-1) are needed to probe the sites of early metal enrichment and a corresponding change in the FeII/MgII ratio.

preprint2020arXiv

Pōniuā&#39;ena: A Luminous $z=7.5$ Quasar Hosting a 1.5 Billion Solar Mass Black Hole

We report the discovery of a luminous quasar, J1007+2115 at $z=7.515$ (&#34;Pōniuā&#39;ena&#34;), 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&#39;s Ly$α$ 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 $μ$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.

preprint2020arXiv

PypeIt: The Python Spectroscopic Data Reduction Pipeline

PypeIt is a Python package for semi-automated reduction of astronomical, spectroscopic data. Its algorithms build on decades-long development of previous data reduction pipelines by the developers (Bernstein, Burles, & Prochaska, 2015; Bochanski et al., 2009). The reduction procedure -- including a complete list of the input parameters and available functionality -- is provided as online documentation hosted by Read the Docs, which is regularly updated. (https://pypeit.readthedocs.io/en/latest/). Release v1.0.3 serves the following spectrographs: Gemini/GNIRS, Gemini/GMOS, Gemini/FLAMINGOS 2, Lick/Kast, Magellan/MagE, Magellan/Fire, MDM/OSMOS, Keck/DEIMOS (600ZD, 830G, 1200G), Keck/LRIS, Keck/MOSFIRE (J and Y gratings tested), Keck/NIRES, Keck/NIRSPEC (low-dispersion), LBT/Luci-I, Luci-II, LBT/MODS (beta), NOT/ALFOSC (grism4), VLT/X-Shooter (VIS, NIR), VLT/FORS2 (300I, 300V), WHT/ISIS.

preprint2020arXiv

Revealing the Accretion Physics of Supermassive Black Holes at Redshift z~7 with Chandra and Infrared Observations

X-ray emission from quasars has been detected up to redshift $z=7.5$, although only limited to a few objects at $z>6.5$. In this work, we present new Chandra observations of five $z>6.5$ quasars. By combining with archival Chandra observations of six additional $z>6.5$ quasars, we perform a systematic analysis on the X-ray properties of these earliest accreting supermassive black holes (SMBHs). We measure the black hole masses, bolometric luminosities ($L_{\rm bol}$), Eddington ratios ($λ_{\rm Edd}$), emission line properties, and infrared luminosities ($L_{\rm IR}$) of these quasars using infrared and sub-millimeter observations. Correlation analysis indicates that the X-ray bolometric correction (the factor that converts from X-ray luminosity to bolometric luminosity) decreases with increasing $L_{\rm bol}$, and that the UV/optical-to-X-ray ratio, $α_{\rm ox}$, strongly correlates with $L_{\rm 2500}$, and moderately correlates with $λ_{\rm Edd}$ and blueshift of CIV emission lines. These correlations are consistent with those found in lower-$z$ quasars, indicating quasar accretion physics does not evolve with redshift. We also find that $L_{\rm IR}$ does not correlate with $L_{\rm 2-10 keV}$ in these luminous distant quasars, suggesting that the ratio of the SMBH growth rate and their host galaxy growth rate in these early luminous quasars are different from those of local galaxies. A joint spectral analysis of the X-ray detected $z>6.5$ quasars yields an average X-ray photon index of $ Γ=2.32^{+0.31}_{-0.30}$, steeper than that of low-$z$ quasars. By comparing it with the $Γ-λ_{\rm Edd}$ relation, we conclude that the steepening of $Γ$ for quasars at $z>6.5$ is mainly driven by their higher Eddington ratios.

preprint2019arXiv

Ionization bias and the ghost proximity effect near $z\gtrsim6$ quasars in the shadow of proximate absorption systems

The larger-than-expected scatter in the opacity of the Ly$α$ forest suggests that the metagalactic ionizing background is strongly fluctuating at $z > 5.5$. Models for ionizing background fluctuations predict a strong positive bias on large scales, so the environments of massive $>10^{12}\,{\rm M}_\odot$ dark matter halos, e.g. $z\sim6$ quasar hosts, would be ideal laboratories to constrain the sources of ionizing photons. While the quasars themselves should overwhelm any plausible ionizing photon contribution from neighboring galaxies, proximate damped Ly$α$ absorbers (DLAs) have recently been discovered in the foreground of $z\sim6$ quasars, and the Ly$α$ forest in the shadow of these DLAs could probe the local ionization environment. Using Gpc$^3$ simulations of $z=6$ ionizing background fluctuations, we show that while the Ly$α$ forest signal from ionization bias around a quasar host halo should be strong, it is likely suppressed by the associated intergalactic matter overdensity. We also show that the quasar itself may still overwhelm the clustering signal via a &#34;ghost&#34; of the proximity effect from the quasar radiation causing a large-scale bias in the ionizing photon mean free path. This ghost proximity effect is sensitive to the lifetime and geometry of quasar emission, potentially unlocking a new avenue for constraining these fundamental quasar properties. Finally, we present observations of a $z\sim6$ quasar with a proximate DLA which shows a strong excess in Ly$α$ forest transmission at the predicted location of the ghost proximity effect.

preprint2019arXiv

The REQUIEM Survey I: A Search for Extended Ly-Alpha Nebular Emission Around 31 z>5.7 Quasars

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.

preprint2019arXiv

Time-dependent behaviour of quasar proximity zones at $z \sim 6$

Since the discovery of $z\sim 6$ quasars two decades ago, studies of their Ly$α$-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$α$ 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 &#34;lightbulb&#34; (i.e. steady constant emission) light curves.