Subjects: Astronomy >> Astrophysical processes submitted time 2023-02-21
Yali Shao
Ran Wang
Axel Weiss
Jeff Wagg
Chris L. Carilli
Michael A. Strauss
Fabian Walter
Pierre Cox
Xiaohui Fan
Karl M. Menten
Desika Narayanan
Dominik Riechers
Frank Bertoldi
Alain Omont
Linhua Jiang
Abstract： We present Atacama Large Millimeter/submillimeter Array (ALMA) sub-kiloparsec- to kiloparsec-scale resolution observations of the [C ii], CO (9–8), and OH+ (11–01) lines along with their dust continuum emission toward the far-infrared (FIR) luminous quasar SDSS J231038.88+185519.7 at z = 6.0031, to study the interstellar medium distribution, the gas kinematics, and the quasar-host system dynamics. We decompose the intensity maps of the [C ii] and CO (9–8) lines and the dust continuum with two-dimensional elliptical Sérsic models. The [C ii] brightness follows a flat distribution with a Sérsic index of 0.59. The CO (9–8) line and the dust continuum can be fit with an unresolved nuclear component and an extended Sérsic component with a Sérsic index of ∼1, which may correspond to the emission from an active galactic nucleus dusty molecular torus and a quasar host galaxy, respectively. The different [C ii] spatial distribution may be due to the effect of the high dust opacity, which increases the FIR background radiation on the [C ii] line, especially in the galaxy center, significantly suppressing the [C ii] emission profile. The dust temperature drops with distance from the center. The effective radius of the dust continuum is smaller than that of the line emission and the dust mass surface density, but is consistent with that of the star formation rate surface density. This may indicate that the dust emission is a less robust tracer of the dust and gas distribution but is a decent tracer of the obscured star formation activity. The OH+ (11–01) line shows a P-Cygni profile with an absorption at ∼–400 km/s, which may indicate an outflow with a neutral gas mass of (6.2 ± 1.2) × 108 M along the line of sight. We employed a three-dimensional tilted ring model to fit the [C ii] and CO (9–8) data cubes. The two lines are both rotation dominated and trace identical disk geometries and gas motions. This suggest that the [C ii] and CO (9–8) gas are coplanar and corotating in this quasar host galaxy. The consistent circular velocities measured with [C ii] and CO (9–8) lines indicate that these two lines trace a similar gravitational potential. We decompose the circular rotation curve measured from the kinematic model fit to the [C ii] line into four matter components (black hole, stars, gas, and dark matter). The quasar-starburst system is dominated by baryonic matter inside the central few kiloparsecs. We constrain the black hole mass to be 2.97+0.51 −0.77 × 109 M ; this is the first time that the dynamical mass of a black hole has been measured at z ∼ 6. This mass is consistent with that determined using the scaling relations from quasar emission lines. A massive stellar component (on the order of 109 M ) may have already existed when the Universe was only ∼0.93 Gyr old. The relations between the black hole mass and the baryonic mass of this quasar indicate that the central supermassive black hole may have formed before its host galaxy.
Peer Review Status:
Awaiting Review
Subjects: Astronomy >> Astrophysical processes submitted time 2023-02-19
Zhilei Xu
Jacqueline N. Hewitt
Kai-Feng Chen
Honggeun Kim
Joshua S. Dillon
Nicholas S. Kern
Miguel F. Morales
Bryna J. Hazelton
Ruby Byrne
Nicolas Fagnoni
Eloy de Lera Acedo
Zara Abdurashidova
Tyrone Adams
James E. Aguirre
Paul Alexander
Zaki S. Ali
Rushelle Baartman
Yanga Balfour
Adam P. Beardsley
Gianni Bernardi
Abstract： Motivated by the desire for wide-field images with well-defined statistical properties for 21cm cosmology, we implement an optimal mapping pipeline that computes a maximum likelihood estimator for the sky using the interferometric measurement equation. We demonstrate this direct optimal mapping with data from the Hydrogen Epoch of Reionization (HERA) Phase I observations. After validating the pipeline with simulated data, we develop a maximum likelihood figure-of-merit for comparing four sky models at 166MHz with a bandwidth of 100kHz. The HERA data agree with the GLEAM catalogs to <10%. After subtracting the GLEAM point sources, the HERA data discriminate between the different continuum sky models, providing most support for the model of Byrne et al. 2021. We report the computation cost for mapping the HERA Phase I data and project the computation for the HERA 320-antenna data; both are feasible with a modern server. The algorithm is broadly applicable to other interferometers and is valid for wide-field and non-coplanar arrays.
Peer Review Status:Awaiting Review
Subjects: Astronomy >> Astrophysical processes submitted time 2023-02-19
Adélie Gorce
Samskruthi Ganjam
Adrian Liu
Steven G. Murray
Zara Abdurashidova
Tyrone Adams
James E. Aguirre
Paul Alexander
Zaki S. Ali
Rushelle Baartman
Yanga Balfour
Adam P. Beardsley
Gianni Bernardi
Tashalee S. Billings
Judd D. Bowman
Richard F. Bradley
Philip Bull
Jacob Burba
Steven Carey
Chris L. Carilli
Abstract： Combining the visibilities measured by an interferometer to form a cosmological power spectrum is a complicated process. In a delay-based analysis, the mapping between instrumental and cosmological space is not a one-to-one relation. Instead, neighbouring modes contribute to the power measured at one point, with their respective contributions encoded in the window functions. To better understand the power measured by an interferometer, we assess the impact of instrument characteristics and analysis choices on these window functions. Focusing on the Hydrogen Epoch of Reionization Array (HERA) as a case study, we find that long-baseline observations correspond to enhanced low-k tails of the window functions, which facilitate foreground leakage, whilst an informed choice of bandwidth and frequency taper can reduce said tails. With simple test cases and realistic simulations, we show that, apart from tracing mode mixing, the window functions help accurately reconstruct the power spectrum estimator of simulated visibilities. The window functions depend strongly on the beam chromaticity, and less on its spatial structure - a Gaussian approximation, ignoring side lobes, is sufficient. Finally, we investigate the potential of asymmetric window functions, down-weighting the contribution of low-k power to avoid foreground leakage. The window functions presented here correspond to the latest HERA upper limits for the full Phase I data. They allow an accurate reconstruction of the power spectrum measured by the instrument and will be used in future analyses to confront theoretical models and data directly in cylindrical space.
Peer Review Status:Awaiting Review
Subjects: Astronomy >> Astrophysical processes submitted time 2023-02-19
Pascal M. Keller
Bojan Nikolic
Nithyanandan Thyagarajan
Chris L. Carilli
Gianni Bernardi
Ntsikelelo Charles
Landman Bester
Oleg M. Smirnov
Nicholas S. Kern
Joshua S. Dillon
Bryna J. Hazelton
Miguel F. Morales
Daniel C. Jacobs
Aaron R. Parsons
Zara Abdurashidova
Tyrone Adams
James E. Aguirre
Paul Alexander
Zaki S. Ali
Rushelle Baartman
Abstract： Radio interferometers aiming to measure the power spectrum of the redshifted 21 cm line during the Epoch of Reionisation (EoR) need to achieve an unprecedented dynamic range to separate the weak signal from overwhelming foreground emissions. Calibration inaccuracies can compromise the sensitivity of these measurements to the effect that a detection of the EoR is precluded. An alternative to standard analysis techniques makes use of the closure phase, which allows one to bypass antenna-based direction-independent calibration. Similarly to standard approaches, we use a delay spectrum technique to search for the EoR signal. Using 94 nights of data observed with Phase I of the Hydrogen Epoch of Reionization Array (HERA), we place approximate constraints on the 21 cm power spectrum at $z=7.7$. We find at 95% confidence that the 21 cm EoR brightness temperature is $\le$(372)$^2$ "pseudo" mK$^2$ at 1.14 "pseudo" $h$ Mpc$^{-1}$, where the "pseudo" emphasises that these limits are to be interpreted as approximations to the actual distance scales and brightness temperatures. Using a fiducial EoR model, we demonstrate the feasibility of detecting the EoR with the full array. Compared to standard methods, the closure phase processing is relatively simple, thereby providing an important independent check on results derived using visibility intensities, or related.
Peer Review Status:Awaiting Review
Subjects: Astronomy >> Astrophysical processes submitted time 2023-02-19
The HERA Collaboration
Zara Abdurashidova
Tyrone Adams
James E. Aguirre
Paul Alexander
Zaki S. Ali
Rushelle Baartman
Yanga Balfour
Rennan Barkana
Adam P. Beardsley
Gianni Bernardi
Tashalee S. Billings
Judd D. Bowman
Richard F. Bradley
Daniela Breitman
Philip Bull
Jacob Burba
Steve Carey
Chris L. Carilli
Carina Cheng
Abstract： We report the most sensitive upper limits to date on the 21 cm epoch of reionization power spectrum using 94 nights of observing with Phase I of the Hydrogen Epoch of Reionization Array (HERA). Using similar analysis techniques as in previously reported limits (HERA Collaboration 2022a), we find at 95% confidence that $\Delta^2(k = 0.34$ $h$ Mpc$^{-1}$) $\leq 457$ mK$^2$ at $z = 7.9$ and that $\Delta^2 (k = 0.36$ $h$ Mpc$^{-1}) \leq 3,496$ mK$^2$ at $z = 10.4$, an improvement by a factor of 2.1 and 2.6 respectively. These limits are mostly consistent with thermal noise over a wide range of $k$ after our data quality cuts, despite performing a relatively conservative analysis designed to minimize signal loss. Our results are validated with both statistical tests on the data and end-to-end pipeline simulations. We also report updated constraints on the astrophysics of reionization and the cosmic dawn. Using multiple independent modeling and inference techniques previously employed by HERA Collaboration (2022b), we find that the intergalactic medium must have been heated above the adiabatic cooling limit at least as early as $z = 10.4$, ruling out a broad set of so-called "cold reionization" scenarios. If this heating is due to high-mass X-ray binaries during the cosmic dawn, as is generally believed, our result's 99% credible interval excludes the local relationship between soft X-ray luminosity and star formation and thus requires heating driven by evolved low-metallicity stars.
Peer Review Status:Awaiting Review
Subjects: Astronomy >> Astrophysical processes submitted time 2023-02-19
Seiji Fujimoto
Steven L. Finkelstein
Denis Burgarella
Chris L. Carilli
Véronique Buat
Caitlin M. Casey
Laure Ciesla
Sandro Tacchella
Jorge A. Zavala
Gabriel Brammer
Yoshinobu Fudamoto
Masami Ouchi
Francesco Valentino
M. C. Cooper
Mark Dickinson
Maximilien Franco
Mauro Giavalisco
Taylor A. Hutchison
Jeyhan S. Kartaltepe
Anton M. Koekemoer
Abstract： We present ALMA Band 7 observations of a remarkably bright galaxy candidate at $z_{\rm phot}$=$16.7^{+1.9}_{-0.3}$ ($M_{\rm UV}$=$-21.6$), S5-z17-1, identified in JWST Early Release Observation data of Stephen's Quintet. We do not detect the dust continuum at 866 $\mu$m, ruling out the possibility that S5-z17-1 is a low-$z$ dusty starburst with a star-formation rate (SFR) of $\gtrsim 30\,M_{\odot}$ yr$^{-1}$. We detect a 5.1$\sigma$ line feature at $338.726\pm0.007$ GHz exactly coinciding with the JWST source position, with a 2% likelihood of the signal being spurious. The most likely line identification would be [OIII]52$\mu$m at $z=16.01$ or [CII]158$\mu$m at $z=4.61$, whose line luminosities do not violate the non-detection of the dust continuum in both cases. Together with three other $z\gtrsim$ 11-13 candidate galaxies recently observed with ALMA, we conduct a joint ALMA and JWST spectral energy distribution (SED) analysis and find that the high-$z$ solution at $z\sim$11-17 is favored in every candidate as a very blue (UV continuum slope of $\approx-2.3$) and luminous ($M_{\rm UV}$ $\approx$ [$-$24:$-$21]) system. Still, we find in some candidates that reasonable SED fits ($\Delta$ $\chi^{2}\lesssim4$) are reproduced by type-II quasar and/or quiescent galaxy templates with strong emission lines at $z\sim3$-5, where such populations predicted from their luminosity functions and EW([OIII]+H$\beta$) distributions are abundant in survey volumes used for the $z\sim$11-17 candidates. While these recent ALMA observation results have strengthened the likelihood of the high-$z$ solutions, lower-$z$ possibilities are not completely ruled out in some of the $z\sim$11-17 candidates.
Peer Review Status:Awaiting Review
Subjects: Astronomy >> Astrophysical processes submitted time 2023-02-19
Yuanqi Liu
Ran Wang
Emmanuel Momjian
Jeff Wagg
Xiaolong Yang
Tao An
Yali Shao
Chris L. Carilli
Xuebing Wu
Xiaohui Fan
Fabian Walter
Linhua Jiang
Qiong Li
Jianan Li
Qinyue Fei
Fuxiang Xu
Abstract： 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 $\alpha = -0.52\pm0.18$ ($S \propto \nu^{\alpha}$). 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.
Peer Review Status:Awaiting Review
Subjects: Astronomy >> Astrophysical processes submitted time 2023-02-19
Michael Pagano
Jing Liu
Adrian Liu
Nicholas S. Kern
Aaron Ewall-Wice
Philip Bull
Robert Pascua
Siamak Ravanbakhsh
Zara Abdurashidova
Tyrone Adams
James E. Aguirre
Paul Alexander
Zaki S. Ali
Rushelle Baartman
Yanga Balfour
Adam P. Beardsley
Gianni Bernardi
Tashalee S. Billings
Judd D. Bowman
Richard F. Bradley
Abstract： Radio Frequency Interference (RFI) is one of the systematic challenges preventing 21cm interferometric instruments from detecting the Epoch of Reionization. To mitigate the effects of RFI on data analysis pipelines, numerous inpaint techniques have been developed to restore RFI corrupted data. We examine the qualitative and quantitative errors introduced into the visibilities and power spectrum due to inpainting. We perform our analysis on simulated data as well as real data from the Hydrogen Epoch of Reionization Array (HERA) Phase 1 upper limits. We also introduce a convolutional neural network that capable of inpainting RFI corrupted data in interferometric instruments. We train our network on simulated data and show that our network is capable at inpainting real data without requiring to be retrained. We find that techniques that incorporate high wavenumbers in delay space in their modeling are best suited for inpainting over narrowband RFI. We also show that with our fiducial parameters Discrete Prolate Spheroidal Sequences (DPSS) and CLEAN provide the best performance for intermittent ``narrowband'' RFI while Gaussian Progress Regression (GPR) and Least Squares Spectral Analysis (LSSA) provide the best performance for larger RFI gaps. However we caution that these qualitative conclusions are sensitive to the chosen hyperparameters of each inpainting technique. We find these results to be consistent in both simulated and real visibilities. We show that all inpainting techniques reliably reproduce foreground dominated modes in the power spectrum. Since the inpainting techniques should not be capable of reproducing noise realizations, we find that the largest errors occur in the noise dominated delay modes. We show that in the future, as the noise level of the data comes down, CLEAN and DPSS are most capable of reproducing the fine frequency structure in the visibilities of HERA data.
Peer Review Status:Awaiting Review
Subjects: Astronomy >> Astrophysical processes submitted time 2023-02-19
Yali Shao
Ran Wang
Axel Weiss
Jeff Wagg
Chris L. Carilli
Michael A. Strauss
Fabian Walter
Pierre Cox
Xiaohui Fan
Karl M. Menten
Desika Narayanan
Dominik Riechers
Frank Bertoldi
Alain Omont
Linhua Jiang
Abstract： We present ALMA sub-kpc- to kpc-scale resolution observations of the [CII], CO(9-8), and OH$^{+}$\,($1_{1}$--$0_{1}$) lines along with their dust continuum emission toward the FIR luminous quasar SDSS J231038.88+185519.7 at $z = 6.0031$. The [CII] brightness follows a flat distribution with a Sersic index of 0.59. The CO(9-8) line and the dust continuum can be fit with an unresolved nuclear component and an extended Sersic component with a Sersic index of ~1. The dust temperature drops with distance from the center. The effective radius of the dust continuum is smaller than that of the line emission and the dust mass surface density, but is consistent with that of the star formation rate surface density. The OH$^{+}$\,($1_{1}$--$0_{1}$) line shows a P-Cygni profile with an absorption, which may indicate an outflow with a neutral gas mass of $(6.2\pm1.2)\times10^{8} M_{\odot}$ along the line of sight. We employed a 3D tilted ring model to fit the [CII] and CO(9-8) data cubes. The two lines are both rotation dominated and trace identical disk geometries and gas motions. We decompose the circular rotation curve measured from the kinematic model fit to the [CII] line into four matter components (black hole, stars, gas, and dark matter). The quasar-starburst system is dominated by baryonic matter inside the central few kiloparsecs. We constrain the black hole mass to be $2.97^{+0.51}_{-0.77}\times 10^{9}\,M_{\odot}$; this is the first time that the dynamical mass of a black hole has been measured at $z\sim6$. A massive stellar component (on the order of $10^{9}\,M_{\odot}$) may have already existed when the Universe was only ~0.93 Gyr old. The relations between the black hole mass and the baryonic mass of this quasar indicate that the central supermassive black hole may have formed before its host galaxy. [Abridged version. Please see the full abstract in the manuscript.]
Peer Review Status:Awaiting Review
Hits
Hits
Downloads
Comment