Researcher profile

Xiaohui Fan

Xiaohui Fan contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
56works
0followers
9topics
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

56 published item(s)

preprint2026arXiv

Debunking Grad-ECLIP: A Comprehensive Study on Its Incorrectness and Fundamental Principles for Model Interpretation

Grad-ECLIP is published at ICML 2024 and represents a new Transformer interpretation technical route (intermediate features-based). First, this paper demonstrates that the intermediate features-based technical route is not a novel one. Based on the existing attention-based route, we have developed Attention-ECLIP, which is completely equivalent to Grad-ECLIP but with simpler computation. Both through formal derivation and experimental validation, we prove that the intermediate feature-based route represented by Grad-ECLIP is actually an equivalent variant of the attention-based route. Next, this paper demonstrates that the Grad-ECLIP method is flawed. The model interpretation results obtained by Grad-ECLIP are not those of the original model, and the interpretation results are misaligned with the model's performance. We analyze the causes of Grad-ECLIP's flaws and propose, or rather, explicitly emphasize two fundamental principles that model interpretation should adhere to in order to avoid similar errors.

preprint2026arXiv

Transformer Interpretability from Perspective of Attention and Gradient

Although researchers' attention is more focused on the performance of Transformer models, the interpretation of Transformer can never be ignored. Gradient is widely utilized in Transformer interpretation. From the perspective of attention and gradient, we conduct an in-depth study of Transformer interpretation and propose a method to achieve it by guiding the gradient direction, or more precisely, the attention direction. The method enables more comprehensive interpretation of feature regions, offers detail interpretation, and helps to better understand Transformer mechanism. Leveraging the difference in how Vision Transformer (ViT) and humans perceive images, we alter the class of an image in a way that is almost imperceptible to the human eye. This class rewriting phenomenon may potentially pose security risks in certain scenarios.

preprint2025arXiv

COSMOS2025: The COSMOS-Web galaxy catalog of photometry, morphology, redshifts, and physical parameters from JWST, HST, and ground-based imaging

We present COSMOS2025, the COSMOS-Web catalog of photometry, morphology, photometric redshifts and physical parameters for more than 700,000 galaxies in the Cosmic Evolution Survey (COSMOS) field. This catalog is based on our \textit{James Webb Space Telescope} 255\,h COSMOS-Web program, which provides deep near-infrared imaging in four NIRCam (F115W, F150W, F277W, F444W) and one MIRI (F770W) filter over the central $\sim 0.54 {\, \rm deg}^2$ ($\sim 0.2 {\, \rm deg}^2$ for MIRI) in COSMOS. These data are combined with ground- and space-based data to derive photometric measurements of NIRCam-detected sources using both fixed-aperture photometry (on the space-based bands) and a profile-fitting technique on all 37 bands spanning 0.3-8 micron. We provide morphology for all sources from complementary techniques including profile fitting and machine-learning classification. We derive photometric redshifts, physical parameters and non-parametric star formation histories from spectral energy distribution (SED) fitting. The catalog has been extensively validated against previous COSMOS catalogs and other surveys. Photometric redshift accuracy measured using spectroscopically confirmed galaxies out to $z\sim9$ reaches $σ_{\rm MAD} = 0.012$ at $m_{\rm F444W}<28$ and remains at $σ_{\rm MAD} \lesssim 0.03$ as a function of magnitude, color, and galaxy type. This represents a factor of $\sim 2$ improvement at 26 AB mag compared to COSMOS2020. The catalog is approximately 80\% complete at $\log(M_{\star}/{\rm M}_{\odot}) \sim 9$ at $z \sim 10$ and at $\log(M_{\star}/{\rm M}_{\odot}) \sim 7$ at $z \sim 0.2$, representing a gain of 1\,dex compared to COSMOS2020. COSMOS2025 represents the definitive COSMOS-Web catalog. It is provided with complete documentation, together with redshift probability distributions, and it is ready for scientific exploitation today.

preprint2023arXiv

Flares in the changing look AGN Mrk 590. I: The UV response to X-ray outbursts suggests a more complex reprocessing geometry than a standard disk

Mrk 590 is a known changing-look AGN which almost turned off in 2012, and then in 2017 partially re-ignited into a repeat flaring state, unusual for an AGN. Our \emph{Swift} observations since 2013 allow us to characterise the accretion-generated emission and its reprocessing in the central engine of a changing-look AGN. The X-ray and UV variability amplitudes are higher than those typically observed in `steady-state&#39; AGN at similar moderate accretion rates; instead, the variability is similar to that of highly accreting AGN. The unusually strong X-ray to UV correlation suggests that the UV-emitting region is directly illuminated by X-ray outbursts. We find evidence that the X-rays are reprocessed by two UV components, with the dominant one at $\sim$3 days and a faint additional reprocessor at near-zero lag. However, we exclude a significant contribution from diffuse broad line region continuum, known to contribute for bona-fide AGN. A near-zero lag is expected for a standard `lamp-post&#39; disk reprocessing model with a driving continuum source near the black hole. That the overall UV response is dominated by the $\sim$3-day lagged component suggests a complicated reprocessing geometry, with most of the UV continuum not produced in a compact disk, as also found in recent studies of NGC 5548 and NGC 4151. Nonetheless, the observed flares display characteristic timescales of $\sim$100 rest-frame days, consistent with the expected thermal timescale in an accretion disk.

preprint2023arXiv

MAMMOTH-Subaru III. Ly$α$ Halo Extended to $\sim200$ kpc Identified by Stacking $\sim 3300$ Ly$α$ Emitters at $z=2.2-2.3$

In this paper, we present a Ly$α$ halo extended to $\sim200$ kpc identified by stacking $\sim 3300$ Ly$α$ emitters at $z=2.2-2.3$. We carry out imaging observations and data reduction with Subaru/Hyper Suprime-Cam (HSC). Our total survey area is $\sim12$ deg$^2$ and imaging depths are $25.5-27.0$ mag. Using the imaging data, we select 1,240 and 2,101 LAE candidates at $z=2.2$ and 2.3, respectively. We carry out spectroscopic observations of our LAE candidates and data reduction with Magellan/IMACS to estimate the contamination rate of our LAE candidates. We find that the contamination rate of our sample is low (8%). We stack our LAE candidates with a median stacking method to identify the Ly$α$ halo at $z=2$. We show that the Ly$α$ halo is extended to $\sim200$ kpc at a surface brightness level of $10^{-20}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. Comparing to previous studies, our Ly$α$ halo is more extended at radii of $\sim25-100$ kpc, which is not likely caused by the contamination in our sample but by different redshifts and fields instead. To investigate how central galaxies affect surrounding LAHs, we divide our LAEs into subsamples based on the Ly$α$ luminosity ($L_{\rm Lyα}$), rest-frame Ly$α$ equivalent width (EW$_0$), and UV magnitude (M$_{\rm uv}$). We stack the subsamples and find that higher $L_{\rm Lyα}$, lower EW$_0$, and brighter M$_{\rm uv}$ cause more extended halos. Our results suggest that more massive LAEs generally have more extended Ly$α$ halos.

preprint2023arXiv

MAMMOTH-Subaru V. Effects of Cosmic Variance on Ly$α$ Luminosity Functions at $z=2.2-2.3$

Cosmic variance introduces significant uncertainties into galaxy number density properties when surveying the high-z Universe with a small volume, such uncertainties produce the field-to-field variance of galaxy number $σ_{g}$ in observational astronomy. This uncertainty significantly affects the Luminosity Functions (LF) measurement of Lya Emitters (LAEs). For most previous Lya LF studies, $σ_{g}$ is often adopted from predictions by cosmological simulations, but barely confirmed by observations. Measuring cosmic variance requires a huge sample over a large volume, exceeding the capabilities of most astronomical instruments. In this study, we demonstrate an observational approach for measuring the cosmic variance contribution for $z\approx2.2$ Lya LFs. The LAE candidates are observed using narrowband and broadband of the Subaru/Hyper Suprime-Cam (HSC), with 8 independent fields, making the total survey area $\simeq11.62$deg$^2$ and a comoving volume of $\simeq8.71\times10^6$Mpc$^3$. These eight fields are selected using the project of MAMMOTH. We report a best-fit Schechter function with parameters $α=-1.75$ (fixed), $L_{Lyα}^{*}=5.18_{-0.40}^{+0.43} \times 10^{42}$erg s$^{-1}$ and $ϕ_{Lya}^{*}=4.87_{-0.55}^{+0.54}\times10^{-4}$Mpc$^{-3}$ for the overall Lya LFs. After clipping out the regions that can bias the cosmic variance measurements, we calculate $σ_{g}$, by sampling LAEs within multiple pointings assigned on the field image. We investigate the relation between $σ_{g}$ and survey volume $V$, and fit a simple power law: $σ_g=k\times(\frac{V_{\rm eff}}{10^5 {\rm Mpc}^3})^β$. We find best-fit values of $-1.209_{-0.106}^{+0.106}$ for $β$ and $0.986_{-0.100}^{+0.108}$ for k. We compare our measurements with predictions from simulations and find that the cosmic variance of LAEs might be larger than that of general star-forming galaxies.

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

A Mock Catalog of Gravitationally Lensed Quasars for the LSST Survey

We present a mock catalog of gravitationally lensed quasars at $z_\text{qso}<7.5$ with simulated images for the Rubin Observatory Legacy Survey of Space and Time (LSST). We adopt recent measurements of quasar luminosity functions to model the quasar population, and use the CosmoDC2 mock galaxy catalog to model the deflector galaxies, which successfully reproduces the observed galaxy velocity dispersion functions up to $z_d\sim1.5$. The mock catalog is highly complete for lensed quasars with Einstein radius $θ_E>0\farcs07$ and quasar absolute magnitude $M_{i}<-20$. We estimate that there are $\sim10^3$ lensed quasars discoverable in current imaging surveys, and LSST will increase this number to $\sim2.4\times10^3$. Most of the lensed quasars have image separation $Δθ>0\farcs5$, which will at least be marginally resolved in LSST images with seeing of $\sim0\farcs7$. There will be $\sim200$ quadruply-lensed quasars discoverable in the LSST. The fraction of quad lenses among all discoverable lensed quasars is about $\sim10\%-15\%$, and this fraction decreases with survey depth. This mock catalog shows a large diversity in the observational features of lensed quasars, in terms of lensing separation and quasar-to-deflector flux ratio. We discuss possible strategies for a complete search of lensed quasars in the LSST era.

preprint2022arXiv

A Spectroscopic Road Map for Cosmic Frontier: DESI, DESI-II, Stage-5

In this white paper, we present an experimental road map for spectroscopic experiments beyond DESI. DESI will be a transformative cosmological survey in the 2020s, mapping 40 million galaxies and quasars and capturing a significant fraction of the available linear modes up to z=1.2. DESI-II will pilot observations of galaxies both at much higher densities and extending to higher redshifts. A Stage-5 experiment would build out those high-density and high-redshift observations, mapping hundreds of millions of stars and galaxies in three dimensions, to address the problems of inflation, dark energy, light relativistic species, and dark matter. These spectroscopic data will also complement the next generation of weak lensing, line intensity mapping and CMB experiments and allow them to reach their full potential.

preprint2022arXiv

ALMA 200 pc imaging of a z~7 quasar reveals a compact, disk-like host galaxy

We present 0&#34;.035 resolution (~200 pc) imaging of the 158 um [CII] line and the underlying dust continuum of the z=6.9 quasar J234833.34-305410.0. The 18 h ALMA observations reveal extremely compact emission (diameter ~1 kpc) that is consistent with a simple, almost face-on, rotation-supported disk with a significant velocity dispersion of ~160 km/s. The gas mass in just the central 200 pc is ~4x10^9 M_sun, about a factor two higher than that of the central supermassive black hole. Consequently we do not resolve the black hole&#39;s sphere of influence, and find no kinematic signature of the central supermassive black hole. Kinematic modeling of the [CII] line shows that the dynamical mass at large radii is consistent with the gas mass, leaving little room for a significant mass contribution by stars and/or dark matter. The Toomre-Q parameter is less than unity throughout the disk, and thus is conducive to star formation, consistent with the high infrared luminosity of the system. The dust in the central region is optically thick, at a temperature >132 K. Using standard scaling relations of dust heating by star formation, this implies an unprecedented high star formation rate density of >10^4 M_sun / yr / kpc^2. Such a high number can still be explained with the Eddington limit for star formation under certain assumptions, but could also imply that the central supermassive black hole contributes to the heating of the dust in the central 110 pc.

preprint2022arXiv

Chemical Abundance of z~6 Quasar Broad-Line Regions in the XQR-30 Sample

The elemental abundances in the broad-line regions of high-redshift quasars trace the chemical evolution in the nuclear regions of massive galaxies in the early universe. In this work, we study metallicity-sensitive broad emission-line flux ratios in rest-frame UV spectra of 25 high-redshift (5.8 < z < 7.5) quasars observed with the VLT/X-shooter and Gemini/GNIRS instruments, ranging over $\log(M_{\rm{BH}}/M_{\odot})= 8.4-9.8$ in black hole mass and $\log(L_{\rm{bol}}/\rm{erg\, s}^{-1})= 46.7-47.7$ in bolometric luminosity. We fit individual spectra and composites generated by binning across quasar properties: bolometric luminosity, black hole mass, and blueshift of the \civ\, line, finding no redshift evolution in the emission-line ratios by comparing our high-redshift quasars to lower-redshift (2.0 < z < 5.0) results presented in the literature. Using Cloudy-based locally optimally-emitting cloud photoionisation model relations between metallicity and emission-line flux ratios, we find the observable properties of the broad emission lines to be consistent with emission from gas clouds with metallicity that are at least 2-4 times solar. Our high-redshift measurements also confirm that the blueshift of the CIV emission line is correlated with its equivalent width, which influences line ratios normalised against CIV. When accounting for the CIV blueshift, we find that the rest-frame UV emission-line flux ratios do not correlate appreciably with the black hole mass or bolometric luminosity.

preprint2022arXiv

Connecting Low- and High-Redshift Weak Emission-Line Quasars via HST Spectroscopy of Ly$α$ Emission

We present ultraviolet spectroscopy covering the Ly$α$ + N V complex of six candidate low-redshift ($0.9 < z < 1.5$) weak emission-line quasars (WLQs) based on observations with the Hubble Space Telescope. The original systematic searches for these puzzling Type 1 quasars with intrinsically weak broad emission lines revealed an $N \approx 100$ WLQ population from optical spectroscopy of high-redshift ($z > 3$) quasars, defined by a Ly$α$ + N V rest-frame equivalent width (EW) threshold $< 15.4$ Å. Identification of lower-redshift ($z < 3$) WLQ candidates, however, has relied primarily on optical spectroscopy of weak broad emission lines at longer rest-frame wavelengths. With these new observations expanding existing optical coverage into the ultraviolet, we explore unifying the low- and high-$z$ WLQ populations via EW[Ly$α$+NV]. Two objects in the sample unify with high-$z$ WLQs, three others appear consistent with the intermediate portion of the population connecting WLQs and normal quasars, and the final object is consistent with typical quasars. The expanded wavelength coverage improves the number of available line diagnostics for our individual targets, allowing a better understanding of the shapes of their ionizing continua. The ratio of EW[Ly$α$+NV] to EW[MgII] in our sample is generally small but varied, favoring a soft ionizing continuum scenario for WLQs, and we find a lack of correlation between EW[Ly$α$+NV] and the X-ray properties of our targets, consistent with a &#34;slim-disk&#34; shielding gas model. We also find indications that weak absorption may be a more significant contaminant in low-$z$ WLQ populations than previously thought.

preprint2022arXiv

Definitive upper bound on the negligible contribution of quasars to cosmic reionization

Cosmic (hydrogen) reionization marks one of the major phase transitions of the universe at redshift z >= 6. During this epoch, hydrogen atoms in the intergalactic medium (IGM) were ionized by Lyman continuum (LyC) photons. However, it remains challenging to identify the major sources of the LyC photons responsible for reionization. In particular, individual contributions of quasars (or active galactic nuclei, AGNs) and galaxies are still under debate. Here we construct the far-ultraviolet (far-UV) luminosity function for type 1 quasars at z >= 6 that spans 10 magnitudes (-19 < M_UV < -29), conclusively showing that quasars made a negligible contribution to reionization. We mainly search for quasars in the low-luminosity range of M_UV > -23 mag that is critical to determine quasars&#39; total LyC photon production but has been barely explored previously. We find that the quasar population can only provide less than 7% (95% confidence level) of the total photons needed to keep the universe ionized at z = 6.0 - 6.6. Our result suggests that galaxies, presumably low-luminosity star-forming systems, are the major sources of hydrogen reionization.

preprint2022arXiv

Exploring the radio spectral energy distribution of the ultraluminous radio-quiet quasar SDSS J0100+2802 at redshift 6.3

We report deep Karl G. Jansky Very Large Array (VLA) observations of the optically ultraluminous and radio-quiet quasar SDSS J010013.02 + 280225.8 (hereafter J0100+2802) at redshift $z=$6.3. We detected the radio continuum emission at 1.5 GHz, 6 GHz, and 10 GHz. This leads to a radio power-law spectral index of $α= -0.52\pm0.18$ ($S \propto ν^α$). The radio source is unresolved in all VLA bands with an upper limit to the size of $0.2^{\prime \prime}$ (i.e., $\sim$ 1.1 kpc) at 10 GHz. We find variability in the flux density (increase by $\sim 33\%$) and the spectral index (steepened) between observations in 2016 and 2017. We also find that the VLA 1.5 GHz flux density observed in the same year is 1.5 times that detected with the Very Long Baseline Array (VLBA) in 2016 at the same frequency. This difference suggests that half of the radio emission from J0100+2802 comes from a compact core within 40 pc, and the rest comes from the surrounding few kpc area which is diffuse and resolved out in the VLBA observations. The diffuse emission is four times brighter than that would be expected if driven by star formation. We conclude that the central active galactic nucleus is the dominant power engine of the radio emission in J0100+2802.

preprint2022arXiv

First Census of Gas-phase Metallicity Gradients of Star-forming Galaxies in Overdense Environments at Cosmic Noon

We report the first spatially resolved measurements of gas-phase metallicity radial gradients in star-forming galaxies in overdense environments at $z\gtrsim2$. The spectroscopic data are acquired by the \mg\ survey, a Hubble Space Telescope (HST) cycle-28 medium program. This program is obtaining 45 orbits of WFC3/IR grism spectroscopy in the density peak regions of three massive galaxy protoclusters (BOSS 1244, BOSS 1542 and BOSS 1441) at $z=2-3$. Our sample in the BOSS 1244 field consists of 20 galaxies with stellar-mass ranging from $10^{9.0}$ to $10^{10.3}$ \Msun\ , star formation rate (SFR) from 10 to 240 \Msun\,yr$^{-1}$, and global gas-phase metallicity (\oh) from 8.2 to 8.6. At $1σ$ confidence level, 2/20 galaxies in our sample show positive (inverted) gradients -- the relative abundance of oxygen increasing with galactocentric radius, opposite the usual trend. Furthermore, 1/20 shows negative gradients and 17/20 are consistent with flat gradients. This high fraction of flat/inverted gradients is uncommon in simulations and previous observations conducted in blank fields at similar redshifts. To understand this, we investigate the correlations among various observed properties of our sample galaxies. We find an anticorrelation between metallicity gradient and global metallicity of our galaxies residing in extreme overdensities, and a marked deficiency of metallicity in our massive galaxies as compared to their coeval field counterparts. We conclude that the cold-mode gas accretion plays an active role in shaping the chemical evolution of galaxies in the protocluster environments, diluting their central chemical abundance, and flattening/inverting their metallicity gradients.

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

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

Molecular Contrastive Learning with Chemical Element Knowledge Graph

Molecular representation learning contributes to multiple downstream tasks such as molecular property prediction and drug design. To properly represent molecules, graph contrastive learning is a promising paradigm as it utilizes self-supervision signals and has no requirements for human annotations. However, prior works fail to incorporate fundamental domain knowledge into graph semantics and thus ignore the correlations between atoms that have common attributes but are not directly connected by bonds. To address these issues, we construct a Chemical Element Knowledge Graph (KG) to summarize microscopic associations between elements and propose a novel Knowledge-enhanced Contrastive Learning (KCL) framework for molecular representation learning. KCL framework consists of three modules. The first module, knowledge-guided graph augmentation, augments the original molecular graph based on the Chemical Element KG. The second module, knowledge-aware graph representation, extracts molecular representations with a common graph encoder for the original molecular graph and a Knowledge-aware Message Passing Neural Network (KMPNN) to encode complex information in the augmented molecular graph. The final module is a contrastive objective, where we maximize agreement between these two views of molecular graphs. Extensive experiments demonstrated that KCL obtained superior performances against state-of-the-art baselines on eight molecular datasets. Visualization experiments properly interpret what KCL has learned from atoms and attributes in the augmented molecular graphs. Our codes and data are available at https://github.com/ZJU-Fangyin/KCL.

preprint2022arXiv

Molecular gas in z~6 quasar host galaxies

We investigate the molecular gas content of z~6 quasar host galaxies using the IRAM / Northern Extended Millimeter Array. We target the 3mm dust continuum, and the line emission from CO(6-5), CO(7-6), [CI]2-1 in 10 infra-red-luminous quasars that have been previously studied in their 1mm dust continuum and [CII] line emission. We detect CO(7-6) at various degrees of significance in all the targeted sources, thus doubling the number of such detections in z~6 quasars. The 3mm to 1mm flux density ratios are consistent with a modified black body spectrum with a dust temperature $T_{dust}$~47 K and an optical depth $τ_ν$=0.2 at the [CII] frequency. Our study provides us with four independent ways to estimate the molecular gas mass, $M_{H2}$, in the targeted quasars. This allows us to set constraints on various parameters used in the derivation of molecular gas mass estimates, such as the mass per luminosity ratios $α_{CO}$ and $α_{[CII]}$, the gas-to-dust ratio $δ_{g/d}$, and the carbon abundance [C]/H2. Leveraging either on the dust, CO, [CI], or [CII] emission yields mass estimates of the entire sample in the range $M_{H2}$~$10^{10}$ to $10^{11}$ M$_{\odot}$. We compare the observed luminosities of dust, [CII], [CI], and CO(7-6) with predictions from photo-dissociation and X-ray dominated regions. We find that the former provide better model fits to our data, assuming that the bulk of the emission arises from dense ($n_H>10^4$ cm$^{-3}$) clouds with a column density $N_{H}$~$10^{23}$ cm$^{-2}$, exposed to a radiation field with intensity $G_0$~$10^3$ (in Habing units). Our analysis reiterates the presence of massive reservoirs of molecular gas fueling star formation and nuclear accretion in $z$~6 quasar host galaxies. It also highlights the power of combined 3mm and 1mm observations for quantitative studies of the dense gas content in massive galaxies at cosmic dawn.

preprint2022arXiv

Radio and Far-IR Emission Associated with a Massive Star-forming Galaxy Candidate at z$\simeq$6.8: A Radio-Loud AGN in the Reionization Era?

We report the identification of radio (0.144-3 GHz), mid-IR, far-IR, and sub-mm (24-850$μ$m) emission at the position of one of 41 UV-bright (M$_\mathrm{UV}^{}\lesssim-21.25$) $z\simeq6.6-6.9$ Lyman-break galaxy candidates in the 1.5 deg$^2$ COSMOS field. This source, COS-87259, exhibits a sharp flux discontinuity (factor $>$3) between two narrow/intermediate bands at 9450 and 9700 Angstroms and is undetected in all nine bands blueward of 9600 Angstroms, as expected from a Lyman-alpha break at $z\simeq6.8$. The full multi-wavelength (X-ray through radio) data of COS-87529 can be self-consistently explained by a very massive (M$_{\ast}=10^{10.8}$ M$_{\odot}$) and extremely red (rest-UV slope $β=-0.59$) $z\simeq6.8$ galaxy with hyperluminous infrared emission (L$_{\mathrm{IR}}=10^{13.6}$ L$_{\odot}$) powered by both an intense burst of highly-obscured star formation (SFR$\approx$1800 M$_{\odot}$ yr$^{-1}$) and an obscured ($τ_{\mathrm{9.7μm}}=7.7\pm2.5$) radio-loud (L$_{\mathrm{1.4\ GHz}}\approx10^{25.4}$ W Hz$^{-1}$) AGN. The radio emission is compact (1.04$\pm$0.12 arcsec) and exhibits an ultra-steep spectrum between 1.32-3 GHz ($α=-1.57^{+0.22}_{-0.21}$) that flattens at lower frequencies ($α=-0.86^{+0.22}_{-0.16}$ between 0.144-1.32 GHz), consistent with known $z>4$ radio galaxies. We also demonstrate that COS-87259 may reside in a significant (11$\times$) galaxy overdensity at $z\simeq6.6-6.9$, as common for systems hosting radio-loud AGN. Nonetheless, a spectroscopic redshift will ultimately be required to establish the true nature of COS-87259 as we cannot yet completely rule out low-redshift solutions. If confirmed to lie at $z\simeq6.8$, the properties of COS-87259 would be consistent with a picture wherein AGN and highly-obscured star formation activity are fairly common among very massive (M$_{\ast}>10^{10}$ M$_{\odot}$) reionization-era galaxies.

preprint2022arXiv

Submillimetre galaxies in two massive protoclusters at z = 2.24: witnessing the enrichment of extreme starbursts in the outskirts of HAE density peaks

Submillimetre galaxies represent a rapid growth phase of both star formation and massive galaxies. Mapping SMGs in galaxy protoclusters provides key insights into where and how these extreme starbursts take place in connections with the assembly of the large-scale structure in the early Universe. We search for SMGs at 850$\,μm$ using JCMT/SCUBA-2 in two massive protoclusters at $z=2.24$, BOSS1244 and BOSS1542, and detect 43 and 54 sources with $S_{850}>4\,$mJy at the $4σ$ level within an effective area of 264$\,$arcmin$^2$, respectively. We construct the intrinsic number counts and find that the abundance of SMGs is $2.0\pm0.3$ and $2.1\pm0.2$ times that of the general fields, confirming that BOSS1244 and BOSS1542 contain a higher fraction of dusty galaxies with strongly enhanced star formation. The volume densities of the SMGs are estimated to be $\sim15-$30 times the average, significantly higher than the overdensity factor ($\sim 6$) traced by H$α$ emission-line galaxies (HAEs). More importantly, we discover a prominent offset between the spatial distributions of the two populations in these two protoclusters -- SMGs are mostly located around the high-density regions of HAEs, and few are seen inside these regions. This finding may have revealed for the first time the occurrence of violent star formation enhancement in the outskirts of the HAE density peaks, likely driven by the boosting of gas supplies and/or starburst triggering events. Meanwhile, the lack of SMGs inside the most overdense regions at $z\sim2$ implies a transition to the environment disfavouring extreme starbursts.

preprint2022arXiv

The Identification of a Dusty Multiarm Spiral Galaxy at $z=3.06$ with JWST and ALMA

Spiral arms serve crucial purposes in star formation and galaxy evolution. In this paper, we report the identification of A2744-DSG-$z3$, a dusty, multiarm spiral galaxy at $z=3.059$ using the James Webb Space Telescope (JWST) NIRISS imaging and grism spectroscopy. A2744-DSG-$z3$ was discovered as a gravitationally lensed sub-millimeter galaxy with ALMA. This is the most distant stellar spiral structure seen thus far, consistent with cosmological simulations which suggest $z\approx3$ as the epoch when spirals emerge. Thanks to the gravitational lensing and excellent spatial resolution of JWST, the spiral arms are resolved with a spatial resolution of $\approx290$\,pc. Based on SED fitting, the spiral galaxy has a de-lensed star formation rate of $85\pm30 \ M_{\odot}$ yr$^{-1}$, and a stellar mass of $\approx10^{10.6}\ M_{\odot}$, indicating that A2744-DSG-$z3$ is a main-sequence galaxy. After fitting the spiral arms, we find a stellar effective radius ($R_{e, \rm{star}}$) of $5.0\pm1.5$ kpc. Combing with ALMA measurements, we find that the effective radii ratio between dust and stars is $\approx0.4$, similar to {those} of massive SFGs at $z\sim2$, indicating a compact dusty core in A2744-DSG-$z3$. Moreover, this galaxy appears to be living in a group environment: including A2744-DSG-$z3$, at least three galaxies at $z=3.05 - 3.06$ {are} spectroscopically confirmed by JWST/NIRISS and ALMA, residing within a lensing-corrected projected scale of $\approx 70$ kpc. This, along with the asymmetric brightness profile, further suggests that the spiral arms may be triggered by minor merger events at $z\gtrsim3$.

preprint2022arXiv

The MegaMapper: A Stage-5 Spectroscopic Instrument Concept for the Study of Inflation and Dark Energy

In this white paper, we present the MegaMapper concept. The MegaMapper is a proposed ground-based experiment to measure Inflation parameters and Dark Energy from galaxy redshifts at $2<z<5$. In order to achieve path-breaking results with a mid-scale investment, the MegaMapper combines existing technologies for critical path elements and pushes innovative development in other design areas. To this aim, we envision a 6.5-m Magellan-like telescope, with a newly designed wide field, coupled with DESI spectrographs, and small-pitch robots to achieve multiplexing of at least 26,000. This will match the expected achievable target density in the redshift range of interest and provide a 10x capability over the existing state-of the art, without a 10x increase in project budget.

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

Deep XMM-Newton Observations of an X-ray Weak, Broad Absorption Line Quasar at $z=6.5$

We report X-ray observations of the most distant known gravitationally lensed quasar, J0439+1634 at $z=6.52$, which is also a broad absorption line (BAL) quasar, using the XMM-Newton Observatory. With a 130 ks exposure, the quasar is significantly detected as a point source at the optical position with a total of 358$^{+19}_{-19}$ net counts using the EPIC instrument. By fitting a power-law plus Galactic absorption model to the observed spectra, we obtain a spectral slope of $Γ=1.45^{+0.10}_{-0.09}$. The derived optical-to-X-ray spectral slope $α_{\rm{ox}}$ is $-2.07^{+0.01}_{-0.01}$, suggesting that the X-ray emission of J0439+1634 is weaker by a factor of 18 than the expectation based on its 2500 Angstrom luminosity and the average $α_{\rm{ox}}$ vs. luminosity relationship. This is the first time that an X-ray weak BAL quasar at $z>6$ has been observed spectroscopically. Its X-ray weakness is consistent with the properties of BAL quasars at lower redshift. By fitting a model including an intrinsic absorption component, we obtain intrinsic column densities of $N_{\rm{H}}=2.8^{+0.7}_{-0.6}\times10^{23}\,\rm{cm}^{-2}$ and $N_{\rm{H}}= 4.3^{+1.8}_{-1.5}\times10^{23}\,\rm{cm}^{-2}$, assuming a fixed $Γ$ of 1.9 and a free $Γ$, respectively. The intrinsic rest-frame 2--10 keV luminosity is derived as $(9.4-15.1)\times10^{43}\,\rm{erg\,s}^{-1}$, after correcting for lensing magnification ($μ=51.3$). The absorbed power-law model fitting indicates that J0439+1634 is the highest redshift obscured quasar with a direct measurement of the absorbing column density. The intrinsic high column density absorption can reduce the X-ray luminosity by a factor of $3-7$, which also indicates that this quasar could be a candidate of intrinsically X-ray weak quasar.

preprint2021arXiv

Measuring the Density Fields around Bright Quasars at $z\sim 6$ with XQR-30 Spectra

Measuring the density of the intergalactic medium using quasar sightlines in the epoch of reionization is challenging due to the saturation of Lyman-$α$ absorption. Near a luminous quasar, however, the enhanced radiation creates a proximity zone observable in the quasar spectra where the Lyman-$α$ absorption is not saturated. In this study, we use $10$ high-resolution ($R\gtrsim 10,000$) $z\sim 6$ quasar spectra from the extended XQR-30 sample to measure the density field in the quasar proximity zones. We find a variety of environments within $3$ pMpc distance from the quasars. We compare the observed density cumulative distribution function (CDF) with models from the $\textit{Cosmic Reionization on Computers}$ simulation, and find a good agreement between $1.5$ to $3$pMpc from the quasar. This region is far away from the quasar hosts and hence approaching the mean density of the universe, which allows us to use the CDF to set constraints on the cosmological parameter $σ_8=0.6\pm0.3$. The uncertainty is mainly due to the limited number of high-quality quasar sightlines currently available. Utilizing the more than $>200$ known quasars at $z\gtrsim 6$, this method will allow us in the future to tighten the constraint on $σ_8$ to the percent level. In the region closer to the quasar within $1.5$ pMpc, we find the density is higher than predicted in the simulation by $1.23 \pm 0.17$, suggesting the typical host dark matter halo mass of a bright quasar ($M_{\rm 1450}<-26.5$) at $z\sim 6$ is $\log_{\rm 10} (M/M_\odot)=12.5^{+0.4}_{-0.7}$.

preprint2021arXiv

Metallicity in Quasar Broad Line Regions at Redshift $\sim$ 6

Broad line regions (BLRs) in high-redshift quasars provide crucial information of chemical enrichment in the early universe. Here we present a study of BLR metallicities in 33 quasars at redshift $5.7<z<6.4$. Using the near-IR spectra of the quasars obtained from the Gemini telescope, we measure their rest-frame UV emission line flux and calculate flux ratios. We then estimate BLR metallicities with empirical calibrations based on photoionization models. The inferred median metallicity of our sample is a few times the solar value, indicating that the BLR gas had been highly metal-enriched at $z\sim6$. We compare our sample with a low-redshift quasar sample with similar luminosities and find no evidence of redshift evolution in quasar BLR metallicities. This is consistent with previous studies. The Fe II$/$Mg II flux ratio, a proxy for the Fe$/α$ element abundance ratio, shows no redshift evolution as well, further supporting rapid nuclear star formation at $z\sim6$. We also find that the black hole mass-BLR metallicity relation at $z\sim6$ is consistent with the relation measured at $2<z<5$, suggesting that our results are not biased by a selection effect due to this relation.

preprint2021arXiv

Paving the Way for Euclid and JWST via Optimal Selection of High-z Quasars

We introduce a probabilistic approach to select 6<z<8 quasar candidates for spectroscopic follow-up, which is based on density estimation in the high-dimensional space inhabited by the optical and near-infrared photometry. Density distributions are modeled as Gaussian mixtures with principled accounting of errors using the extreme deconvolution (XD) technique, generalizing an approach successfully used to select lower redshift (z<3) quasars. We train the probability density of contaminants on 733,694 7-d flux measurements from the 1076 square degrees overlapping area from the DECaLS (z), VIKING (YJHK), and unWISE (W1W2) imaging surveys, after requiring they dropout of DECaLS g and r, whereas the distribution of high-z quasars are trained on synthetic model photometry. Extensive simulations based on these density distributions and current estimates of the quasar luminosity function indicate that this method achieves a completeness of >75% and an efficiency of >15% for selecting quasars at 6<z<8 with J<21.5. Among the classified sources are 8 known 6<z<7 quasars, of which 2/8 are selected suggesting a completeness ~25%, whereas classifying the 6 known (J<21.5) quasars at z>7 from the entire sky, we select 5/6 or a completeness of ~80%.The failure to select the majority of 6<z<7 quasars arises because our model of quasar SEDs underestimates the scatter in the distribution of fluxes. This new optimal approach to quasar selection paves the way for efficient spectroscopic follow-up of Euclid quasar candidates with ground based telescopes and JWST.

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.

preprint2021arXiv

Quasar UV Luminosity Function at $3.5<z<5.0$ from SDSS Deep Imaging Data

We present a well-designed sample of more than 1000 type 1 quasars at $3.5<z<5$ and derive UV quasar luminosity functions (QLFs) in this redshift range. These quasars were selected using the Sloan Digital Sky Survey (SDSS) imaging data in SDSS Stripe 82 and overlap regions with repeat imaging observations. They are about one magnitude fainter than those found using the SDSS single-epoch data. The spectroscopic observations were conducted by the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) as one of the BOSS ancillary programs. This quasar sample reaches $i\sim21.5$ mag and bridges previous samples from brighter surveys and deeper surveys. We use a $1/V_\mathrm{a}$ method to derive binned QLFs at $3.6<z<4.0$, $4.0<z<4.5$, and $4.5<z<4.9$, and use a double-power law model to parameterize the QLFs. We also combine our data with those in the literature to better constrain the QLFs in the context of a much wider luminosity baseline. We find that the faint-end and bright-end slopes of the QLFs in this redshift range are around $-1.7$ and $-3.7$, respectively, with uncertainties from 0.2$-$0.3 to $>0.5$. The evolution of the QLFs from $z\sim5$ to $3.5$ can be described by a pure density evolution model ($\propto10^{kz}$) and the parameter $k$ is similar to that at $5<z<7$, suggesting a nearly uniform evolution of the quasar density at $z=3.5-7$.

preprint2021arXiv

Revisiting the Lensed Fraction of High-Redshift Quasars

The observed lensed fraction of high-redshift quasars $(\sim0.2\%)$ is significantly lower than previous theoretical predictions $(\gtrsim4\%)$. We revisit the lensed fraction of high-redshift quasars predicted by theoretical models, where we adopt recent measurements of galaxy velocity dispersion functions (VDFs) and explore a wide range of quasar luminosity function (QLF) parameters. We use both analytical methods and mock catalogs which give consistent results. For ordinary QLF parameters and the depth of current high-redshift quasar surveys $(m_z\lesssim22)$, our model suggests a multiply-imaged fraction of $F_\text{multi}\sim 0.4\%-0.8\%$. The predicted lensed fraction is $\sim1\%-6\%$ for the brightest $z_s\sim6$ quasars $(m_z\lesssim19)$, depending on the QLF. The systematic uncertainties of the predicted lensed fraction in previous models can be as large as $2-4$ times and are dominated by the VDF. Applying VDFs from recent measurements decreases the predicted lensed fraction and relieves the tension between observations and theoretical models. Given the depth of current imaging surveys, there are $\sim15$ lensed quasars at $z_s>5.5$ detectable over the sky. Upcoming sky surveys like the LSST survey and the {\em Euclid} survey will find several tens of lensed quasars at this redshift range.

preprint2021arXiv

The mass-metallicity relation at cosmic noon in overdense environments: first results from the MAMMOTH-Grism HST slitless spectroscopic survey

The MAMMOTH-Grism slitless spectroscopic survey is a Hubble Space Telescope (HST) cycle-28 medium program, which is obtaining 45 orbits of WFC3/IR grism spectroscopy in the density peak regions of three massive galaxy protoclusters at $z=2-3$ discovered using the MAMMOTH technique. We introduce this survey by presenting the first measurement of the mass-metallicity relation (MZR) at high redshift in overdense environments via grism spectroscopy. From the completed MAMMOTH-Grism observations in the field of the BOSS1244 protocluster at $z=2.24\pm0.02$, We secure a sample of 36 protocluster member galaxies at $z\sim2.24$, showing strong nebular emission lines ([O III], H$β$ and [O II]) in their G141 spectra. Using the multi-wavelength broad-band deep imaging from HST and ground-based telescopes, we measure their stellar masses in the range of $[10^{9},10^{10.4}]M_\odot$, instantaneous star formation rates (SFR) from 10 to 240$M_\odot yr^{-1}$, and global gas-phase metallicities [$\frac{1}{3}$,1] of solar. Compared with similarly selected field galaxy sample at the same redshift, our galaxies show on average increased SFRs by $\sim$0.06dex and $\sim$0.18dex at $\sim$10$^{10.1}M_\odot$ and $\sim$10$^{9.8}M_\odot$, respectively. Using the stacked spectra of our sample galaxies, we derive the MZR in the BOSS1244 protocluster core as $12+\log({\rm O/H})=(0.136\pm0.018)\times\log(M_\ast/M_\odot)+(7.082\pm0.175)$, showing significantly shallower slope than that in the field. This shallow MZR slope is likely caused by the combined effects of efficient recycling of feedback-driven winds and cold-mode gas accretion in protocluster environments. The former effect helps low-mass galaxies residing in overdensities retain their metal production, whereas the latter effect dilutes the metal content of high-mass galaxies, making them more metal poor than their coeval field counterparts.

preprint2021arXiv

Ultra-luminous high-redshift quasars from SkyMapper -- II. New quasars and the bright end of the luminosity function

We search for ultra-luminous Quasi-Stellar Objects (QSOs) at high redshift using photometry from the SkyMapper Southern Survey Data Release 3 (DR3), in combination with 2MASS, VHS DR6, VIKING DR5, AllWISE, and CatWISE2020, as well as parallaxes and proper motions from Gaia DR2 and eDR3. We report 142 newly discovered Southern QSOs at $3.8<z<5.5$, of which 126 have $M_{145} <-27$ ABmag and are found in a search area of 14,486 deg$^2$. This Southern sample, utilising the Gaia astrometry to offset wider photometric colour criteria, achieves unprecedented completeness for an ultra-luminous QSO search at high redshift. In combination with already known QSOs, we construct a sample that is $>80$ per cent complete for $M_{145}<-27.33$ ABmag at $z=4.7$ and for $M_{145}<-27.73$ ABmag at $z=5.4$. We derive the bright end of the QSO luminosity function at restframe 145 nm for $z=4.7-5.4$ and measure its slope to be $β= -3.60\pm0.37$ and $β= -3.38\pm0.32$ for two different estimates of the faint-end QSO density adopted from the literature. We also present the first $z\sim 5$ QSO luminosity function at restframe 300 nm.

preprint2020arXiv

A Closer Look at Two of the Most Luminous Quasars in the Universe

Ultra-luminous quasars ($M_{1450} \leq -29$) provide us with a rare view into the nature of the most massive and most rapidly accreting supermassive black holes (SMBHs). Following the discovery of two of these extreme sources, J0341${+}$1720 ($M_{1450}=-29.56$, $z=3.71$) and J2125${-}$1719 ($M_{1450}=-29.39$, $z=3.90$), in the Extremely Luminous Quasar Survey (ELQS) and its extension to the Pan-STARRS\,1 footprint (PS-ELQS), we herein present an analysis of their rest-frame UV to optical spectroscopy. Both quasars harbor very massive SMBHs with $M_{\rm{BH}}=6.73_{-0.83}^{+0.75}\times10^{9}\,M_{\odot}$ and $M_{\rm{BH}}=5.45_{-0.55}^{+0.60}\times10^{9}\,M_{\odot}$, respectively, showing evidence of accretion above the Eddington limit ($L_{\rm{bol}}/L_{\rm{Edd}}=2.74_{-0.27}^{+0.39}$ and $L_{\rm{bol}}/L_{\rm{Edd}}=3.01_{-0.30}^{+0.34}$). NOEMA 3 millimeter observations of J0341${+}$1720 reveal a highly star-forming ($\rm{SFR}\approx1500\,M_{\odot}\,\rm{yr}^{-1}$), ultra-luminous infrared galaxy ($L_{\rm{TIR}}\approx1.0\times10^{13}\,L_{\odot}$) host, which, based on an estimate of its dynamical mass, is only ${\sim}30$ times more massive than the SMBH it harbors at its center. As examples of luminous super-Eddington accretion, these two quasars provide support for theories, which explain the existence of billion solar mass SMBHs ${\sim}700$ million years after the Big Bang by moderate super-Eddington growth from standard SMBH seeds.

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

A Thirty-Four Billion Solar Mass Black Hole in SMSS J2157-3602, the Most Luminous Known Quasar

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

preprint2020arXiv

Chandra Detection of Three X-ray Bright Quasars at z>5

We report Chandra detection of three UV bright radio quiet quasars at $z\gtrsim5$. We have collected a sufficient number of photons to extract an X-ray spectrum of each quasar to measure their basic X-ray properties, such as the X-ray flux, power law photon index ($Γ$), and optical-to-X-ray spectral slope ($α_{\rm OX}$). J074749+115352 at $z=5.26$ is the X-ray brightest radio-quiet quasar at $z>5$. It may have a short timescale variation (on a timescale of $\sim3800\rm~s$ in the observer&#39;s frame, or $\sim600\rm~s$ in the rest frame) which is however largely embedded in the statistical noise. We extract phase folded spectra of this quasar. There are two distinguishable states: a &#34;high soft&#34; state with an average X-ray flux $\sim2.7$ times of the &#34;low hard&#34; state, and a significantly steeper X-ray spectral slope ($Γ=2.40_{-0.32}^{+0.33}$ vs $1.78_{-0.24}^{+0.25}$). We also compare the three quasars detected in this paper to other quasar samples. We find that J074749+115352, with a SMBH mass of $M_{\rm SMBH}\approx1.8\times10^9\rm~M_\odot$ and an Eddington ratio of $λ_{\rm Edd}\approx2.3$, is extraordinarily X-ray bright. It has an average $α_{\rm OX}=-1.46\pm0.02$ and a 2-10 keV bolometric correction factor of $L_{\rm bol}/L_{\rm2-10keV}=42.4\pm5.8$, both significantly depart from some well defined scaling relations. We compare $Γ$ of the three quasars to other samples at different redshifts, and do not find any significant redshift evolution based on the limited sample of $z>5$ quasars with reliable measurements of the X-ray spectral properties.

preprint2020arXiv

CIV emission line properties and uncertainties in black hole mass estimates of z ~ 3.5 quasars

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

preprint2020arXiv

Deep Hubble Space Telescope Imaging on the Extended Ly$α$ Emission of a QSO at $z = 2.19$ with Damped Lyman Alpha System as a Natural Coronagraph

Recent observations suggest that proximate damped Ly$α$ (PDLA) systems can be used to study the host galaxies of Quasi-stellar objects (QSOs), because the PDLAs can block the bright point-spread-function (PSF) from central QSOs. Using six-orbits of narrowband imaging with $\it{HST}$/WFC3, we present the first high resolution narrowband image of the Ly$α$ emission in the PDLA trough of the QSO SDSSJ115432.67-021537.9 . We detect one major component and one minor component in the narrowband imaging. Combining the $\it{HST}$/WFC3 imaging with deep Magellan/MagE spectra, we measure that the Ly$α$ flux F$_{\rm{Lyα}}$ = 1.56$\pm0.10 \times10^{-16}$ $\rm{erg\ s^{-1}\ cm^{-2}}$, which is among the luminous ($\approx$ 2.7 L$^{\star}_{\rm{Lyα}}$) Ly$α$ emitters at $\it{z} =$ 2.19. The $\rm{Lyα}$-based star formation rate (SFR) is $\sim$ 7 $\textrm{M_{\sun} \ yr^{-1}}$. These observational results favor that the star formation from the host galaxy could be the main mechanism to power the Ly$α$ emission. This new method sheds new light on the study of the kinematic structure and the spatial distribution of the extended Ly$α$ emitting regions around the QSO host.

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

Ionized and atomic interstellar medium in the z = 6.003 quasar SDSS J2310+1855

Observing the interstellar medium (ISM) in $z \gtrsim 6$ quasars host galaxies is essential for understanding the co-evolution between the supermassive black holes and their hosts. To probe the gas physical conditions and search for imprints of Active Galactic Nuclei (AGN) on the ISM, we report ALMA observations of the $\rm [N\ II]_{122 μm}$ and $\rm [O\ I]_{146 μm}$ lines and the underlying continuum from the $z=6.003$ quasar SDSS J231038.88+185519.7. Together with previous $\rm [C\ II]_{158 μm}$ and $\rm [O\ III]_{88 μm}$ observations, we use the ratios of these fine-structure lines to probe the ISM properties. Similar to other high-$z$ systems, this object exhibits a $\rm [C\ II]_{158 μm}$/$\rm [O\ I]_{146 μm}$ ratio comparable to the lowest values found in local (Ultra) luminous infrared galaxies, suggesting a &#34;warmer&#34; and &#34;denser&#34; gas component compared to typical local systems. The $\rm [O\ III]_{88 μm}$/$\rm [O\ I]_{146 μm}$ ratio is lower than that of other local and high-$z$ systems, indicating a smaller ionized gas fraction in this quasar. The $\rm [O\ III]_{88 μm}$/$\rm [N\ II]_{122 μm}$ ratio is comparable to that of local systems, and suggests a metallicity of $Z/Z_{\odot}$=1.5$-$2.1. Based on the $\rm [N\ II]_{122 μm}$ detection, we estimate that $17\%$ of the $\rm [C\ II]_{158 μm}$ emission is associated with ionized gas. The $\rm [N\ II]_{122 μm}$ line shows a &#34;flux deficit&#34; comparable to local systems. The $\rm [O\ I]_{146 μm}$ line, with a $\rm [O\ I]_{146 μm}$/FIR ratio $\ge 2\times$ than expected from the local relation, indicates no $\rm [O\ I]_{\rm 146 μm}$ deficit. The low $\rm [C\ II]_{158 μm}$/$\rm [O\ I]_{146 μm}$ ratio, together with the high $\rm [O\ I]_{146 μm}$/FIR ratio in J2310+1855, reveals that the warm and dense gas is likely a result of AGN heating to the ISM.

preprint2020arXiv

Limits to Rest-Frame Ultraviolet Emission From Far-Infrared-Luminous z~6 Quasar Hosts

We report on a Hubble Space Telescope search for rest-frame ultraviolet emission from the host galaxies of five far-infrared-luminous $z\simeq{}6$ quasars and the $z=5.85$ hot-dust free quasar SDSS J0005-0006. We perform 2D surface brightness modeling for each quasar using a Markov-Chain Monte-Carlo estimator, to simultaneously fit and subtract the quasar point source in order to constrain the underlying host galaxy emission. We measure upper limits for the quasar host galaxies of $m_J>22.7$ mag and $m_H>22.4$ mag, corresponding to stellar masses of $M_\ast<2\times10^{11}M_\odot$. These stellar mass limits are consistent with the local $M_{\textrm{BH}}$-$M_\ast$ relation. Our flux limits are consistent with those predicted for the UV stellar populations of $z\simeq6$ host galaxies, but likely in the presence of significant dust ($\langle A_{\mathrm{UV}}\rangle\simeq 2.6$ mag). We also detect a total of up to 9 potential $z\simeq6$ quasar companion galaxies surrounding five of the six quasars, separated from the quasars by 1.4&#39;&#39;-3.2&#39;&#39;, or 8.4-19.4 kpc, which may be interacting with the quasar hosts. These nearby companion galaxies have UV absolute magnitudes of -22.1 to -19.9 mag, and UV spectral slopes $β$ of -2.0 to -0.2, consistent with luminous star-forming galaxies at $z\simeq6$. These results suggest that the quasars are in dense environments typical of luminous $z\simeq6$ galaxies. However, we cannot rule out the possibility that some of these companions are foreground interlopers. Infrared observations with the James Webb Space Telescope will be needed to detect the $z\simeq6$ quasar host galaxies and better constrain their stellar mass and dust content.

preprint2020arXiv

Lyman Continuum Escape Fraction in Ly$α$ Emitters at $z\simeq3.1$

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

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

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.

preprint2020arXiv

SCUBA2 High Redshift Bright Quasar Survey: Far-infrared Properties and Weak-line Features

We present a submillimetre continuum survey (&#39;SCUBA2 High rEdshift bRight quasaR surveY&#39;, hereafter SHERRY) of 54 high redshift quasars at $5.6<z<6.9$ with quasar bolometric luminosities in a range of (0.2$-$$ 5)\times10^{14}\,L_{\odot}$, using the Submillimetre Common-User Bolometer Array-2 (SCUBA2) on the James Clerk Maxwell Telescope. About 30% (16/54) of the sources are detected with a typical 850$μ$m rms sensitivity of 1.2 $\rm mJy\,beam^{-1}$ ($S\rm _{ν,850\,μm} = 4$-5 mJy, at $>3.5σ$). The new SHERRY detections indicate far-infrared (FIR) luminosities of $\rm 3.5\times10^{12}$ to $\rm 1.4\times10^{13}$ $L_{\odot}$, implying extreme star formation rates of 90 to 1060 $M_{\odot}$ yr$^{-1}$ in the quasar host galaxies. Compared with $z =$ 2$-$5 samples, the FIR luminous quasars ($L_{\rm FIR} > 10^{13}\,L_{\odot}$) are more rare at $z \sim 6$. The optical/near-infrared (NIR) spectra of these objects show 11% (6/54) of the sources have weak Ly$α$, emission line features, which may relate to different sub-phases of the central active galactic nuclei (AGNs). Our SCUBA2 survey confirms the trend reported in the literature that quasars with submillimeter detections tend to have weaker ultraviolet (UV) emission lines compared to quasars with nondetections. The connection between weak UV quasar line emission and bright dust continuum emission powered by massive star formation may suggest an early phase of AGN-galaxy evolution, in which the broad line region is starting to develop slowly or is shielded from the central ionization source, and has unusual properties such as weak line features or bright FIR emission.

preprint2020arXiv

Spectroscopy of broad absorption line quasars at $3\lesssim z \lesssim 5$ -- I: evidence for quasar winds shaping broad/narrow emission line regions

We present an observational study of 22 broad absorption line quasars (BAL QSOs) at $3\lesssim z \lesssim5$ based on optical/near-IR spectroscopy, aiming to investigate quasar winds and their effects. The near-IR spectroscopy covers the \hb\ and/or \mgii\ broad emission lines (BELs) for these quasars, allowing us to estimate their central black hole (BH) masses in a robust way. We found that our BAL QSOs on average do not have a higher Eddington ratio than that from non-BAL QSOs matched in redshift and/or luminosity. In a subset consisting of seven strong BAL QSOs possessing sub-relativistic BAL outflows, we see the prevalence of large \civ-BEL blueshift ($\sim$3100 km s$^{-1}$) and weak \oiii\ emission (particularly the narrow \oiii$λ$5007 component), indicative of nuclear outflows affecting the narrow emission-line (NEL) regions. In another subset consisting of thirteen BAL QSOs having simultaneous observations of \mgii\ and \hb, we found a strong correlation between 3000~Å and 5000~Å monochromatic luminosity, consistent with that from non-BAL QSOs matched in redshift and luminosity; however, there is no correlation between \mgii\ and \hb\ in FWHM, likely due to nuclear outflows influencing the BEL regions. Our spectroscopic investigations offer strong evidence that the presence of nuclear outflows plays an important role in shaping the BEL/NEL regions of these quasars and possibly, regulating the growth of central supermassive black holes (SMBHs). We propose that BEL blueshift and BAL could be different manifestations of the same outflow system viewed at different sightlines and/or phases.

preprint2020arXiv

Statistical correlation between the distribution of Ly$α$ emitters and IGM HI at $z\sim2.2$ mapped by Subaru/Hyper Suprime-Cam

The correlation between neutral Hydrogen (HI) in the intergalactic medium (IGM) and galaxies now attracts great interests. We select four fields which include several coherently strong Ly$α$ absorption systems at $z\sim2.2$ detected by using background quasars from the whole SDSS/(e)BOSS database. Deep narrow-band and $g$-band imaging are performed using the Hyper Suprime-Cam on the Subaru Telescope. We select out 2,642 Ly$α$ emitter (LAE) candidates at $z=2.177\pm0.023$ down to the Ly$α$ luminosity of $L_{\text{Ly}α}\approx 2 \times 10^{42} {\rm erg~s}^{-1}$ to construct the galaxy overdensity maps, covering an effective area of 5.39 deg$^2$. Combining the sample with the Ly$α$ absorption estimated from 64 (e)BOSS quasar spectra, we find a moderate to strong correlation between the LAE overdensity $δ_{\rm LAE}$ and the effective optical depth $τ_{\rm LoS}$ in line-of-sights, with $P$-value$=0.09\%$ ($<0.01\%$) when the field that contains a significant quasar overdensity is in(ex)cluded. The cross-correlation analysis also clearly suggests that up to $4\pm1$ pMpc, LAEs tend to cluster in the regions rich in HI gas, indicated by the high $τ_{\rm LoS}$, and avoid the low $τ_{\rm LoS}$ region where the HI gas is deficient. By averaging the $τ_{\rm LoS}$ as a function of the projected distance ($d$) to LAEs, we find a $30\%$ excess signal at $2σ$ level at $d<200$ pkpc, indicating the dense HI in circumgalactic medium, and a tentative excess at $400<d<600$ pkpc in IGM regime, corroborating the cross-correlation signal detected at about $0.5$ pMpc. These statistical analyses indicate that galaxy$-$IGM HI correlations exist on scales ranging from several hundred pkpc to several pMpc at $z\sim2.2$.

preprint2020arXiv

Strong Mg II and Fe II Absorbers at 2.2 < z < 6.0

We present a study of strong intervening absorption systems in the near-IR spectra of 31 luminous quasars at $z>5.7$. The quasar spectra were obtained with {\it Gemini} GNIRS that provide continuous wavelength coverage from $\sim$0.9 to $\sim$2.5 $μ$m. We detect 32 strong Mg II doublet absorbers with rest-frame equivalent width $W_r$ ($\lambda2796$) $>1.0$ Åat $2.2 < z < 6.0$. Each Mg II absorber is confirmed by at least two associated Fe II absorption lines in the rest-frame wavelength range of $\sim 1600-2600$ Å. We find that the comoving line density ($dN/dX$) of the strong Fe II-bearing Mg II absorbers decreases towards higher redshift at $z>3$, consistent with previous studies. Compared with strong Mg II absorbers detected in damped Ly$α$ systems at 2 $<z<$ 4, our absorbers are potentially less saturated and show much larger rest-frame velocity widths. This suggests that the gas traced by our absorbers are potentially affected by galactic superwinds. We analyze the {\it Hubble Space Telescope} near-IR images of the quasars and identify possible associated galaxies for our strong absorbers. There are a maximum of two galaxy candidates found within 5&#34; radius of each absorber. The median F105W-band magnitude of these galaxy candidates is 24.8 mag, which is fainter than the $L^*$ galaxy luminosity at $z\sim$ 4. By using our observed $dN/dX$ of strong Mg II absorbers and galaxy candidates median luminosity, we suggest that at high redshift, strong Mg II absorbers tend to have a more disturbed environment but smaller halo size than that at $z <$ 1.

preprint2020arXiv

The Sixteenth Data Release of the Sloan Digital Sky Surveys: First Release from the APOGEE-2 Southern Survey and Full Release of eBOSS Spectra

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library &#34;MaStar&#34;). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

preprint2020arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Photometric g and i Light Curves

The Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) program monitors 849 active galactic nuclei (AGN) both spectroscopically and photometrically. The photometric observations used in this work span over four years and provide an excellent baseline for variability studies of these objects. We present the photometric light curves from 2014 to 2017 obtained by the Steward Observatory&#39;s Bok telescope and the CFHT telescope with MegaCam. We provide details on the data acquisition and processing of the data from each telescope, the difference imaging photometry used to produce the light curves, and the calculation of a variability index to quantify each AGN&#39;s variability. We find that the Welch-Stetson J-index provides a useful characterization of AGN variability and can be used to select AGNs for further study.

preprint2020arXiv

The Third Data Release of the Beijing-Arizona Sky Survey

The Beijing-Arizona Sky Survey (BASS) is a wide and deep imaging survey to cover a 5400 deg$^2$ area in the Northern Galactic Cap with the 2.3m Bok telescope using two filters ($g$ and $r$ bands). The Mosaic $z$-band Legacy Survey (MzLS) covers the same area in $z$ band with the 4m Mayall telescope. These two surveys will be used for spectroscopic targeting of the Dark Energy Spectroscopic Instrument (DESI). The BASS survey observations were completed in 2019 March. This paper describes the third data release (DR3) of BASS, which contains the photometric data from all BASS and MzLS observations between 2015 January and 2019 March. The median astrometric precision relative to {\it Gaia} positions is about 17 mas and the median photometric offset relative to the PanSTARRS1 photometry is within 5 mmag. The median $5σ$ AB magnitude depths for point sources are 24.2, 23.6, and 23.0 mag for $g$, $r$, and $z$ bands, respectively. The photometric depth within the survey area is highly homogeneous, with the difference between the 20\% and 80\% depth less than 0.3 mag. The DR3 data, including raw data, calibrated single-epoch images, single-epoch photometric catalogs, stacked images, and co-added photometric catalogs, are publicly accessible at \url{http://batc.bao.ac.cn/BASS/doku.php?id=datarelease:home}.

preprint2020arXiv

X-ray Observations of a [C II]-bright, z=6.59 Quasar/Companion System

We present deep Chandra observations of PSO J231.6576$-$20.8335, a quasar at redshift z=6.59 with a nearby (${\sim}8$ proper kpc) companion galaxy. ALMA observed both the quasar and companion to be bright in [C II], and the system has significant extended Ly$α$ emission around the quasar, suggesting that a galaxy merger is ongoing. Unlike previous studies of two similar systems, and despite observing the system with Chandra for 140 ks, we do not detect the companion in X-rays. The quasar itself is detected, but only $13.3^{+4.8}_{-3.7}$ net counts are observed. From a basic spectral analysis, the X-ray spectrum of the quasar is soft (hardness ratio of $\mathcal{HR} = -0.60_{-0.27}^{+0.17}$, power-law index of $Γ=2.6^{+1.0}_{-0.9}$), which results in a rest-frame X-ray luminosity comparable to other bright quasars ($L_{2-10} = 1.09^{+2.20}_{-0.70}\times 10^{45}\ \textrm{erg}\ \textrm{s}^{-1}$) despite the faint observed X-ray flux. We highlight two possible interpretations of this result: the quasar has a steep value of $Γ$ -- potentially related to observed ongoing Eddington accretion -- thereby pushing much of the emission out of our observed band, or the quasar has a more normal spectrum ($Γ{\sim}2$) but is therefore less X-ray luminous ($L_{2-10} \sim 0.6 \times 10^{45}\ \textrm{ erg}\ \textrm{ s}^{-1}$).

preprint2019arXiv

Probing the full CO spectral line energy distribution (SLED) in the nuclear region of a quasar-starburst system at $z=6.003$

We report Atacama Large Millimeter/submillimeter Array (ALMA) observations of CO $(8-7)$, $(9-8)$, $\rm H_{2}O (2_{0,2}-1_{1,1})$ and $\rm OH^{+} (1_{1}-0_{1})$ and NOrthern Extended Millimeter Array (NOEMA) observations of CO $(5-4)$, $(6-5)$, $(12-11)$ and $(13-12)$ towards the $z = 6.003$ quasar SDSS J231038.88+185519.7, aiming to probe the physical conditions of the molecular gas content of this source. We present the best sampled CO spectral line energy distribution (SLED) at $z = 6.003$, and analyzed it with the radiative transfer code MOLPOP-CEP. Fitting the CO SLED to a one-component model indicates a kinetic temperature $T_{\rm kin} = 228 \ \rm K$, molecular gas density $log (n(\rm H_{2})/\rm cm^{-3}$ )=4.75, and CO column density $log(N(\rm CO)/\rm cm^{-2}) =17.5$, although a two-component model better fits the data. In either case, the CO SLED is dominated by a &#34;warm&#34; and &#34;dense&#34; component. Compared to samples of local (Ultra) Luminous Infrared Galaxies ((U)LIRGs), starburst galaxies and high redshift Submillimeter Galaxies (SMGs), J2310+1855 exhibits higher CO excitation at ($J \geq 8$), like other high redshift quasars. The high CO excitation, together with the enhanced $L_{\rm H_{2}O}/ L_{IR} $, $L_{\rm H_{2}O}/ L_{CO} $ and $L_{OH^{+}}/L_{\rm H_{2}O} $ ratios, suggests that besides the UV radiation from young massive stars, other mechanisms such as shocks, cosmic rays and X-rays might also be responsible for the heating and ionization of the molecular gas. In the nuclear region probed by the molecular emissions lines, any of these mechanisms might be present due to the powerful quasar and the starburst activity.

preprint2019arXiv

Resolving the interstellar medium in the nuclear region of two z=5.78 quasar host galaxies with ALMA

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

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.