Showing 1–50 of 458 results for author: Carilli, C

A blazar in the epoch of reionization

Authors: Eduardo Banados, Emmanuel Momjian, Thomas Connor, Silvia Belladitta, Roberto Decarli, Chiara Mazzucchelli, Bram P. Venemans, Fabian Walter, Feige Wang, Zhang-Liang Xie, Aaron J. Barth, Anna-Christina Eilers, Xiaohui Fan, Yana Khusanova, Jan-Torge Schindler, Daniel Stern, Jinyi Yang, Irham Taufik Andika, Chris Carilli, Emanuele P. Farina, Andrew Fabian, Joseph F. Hennawi, Antonio Pensabene, Sofia Rojas-Ruiz

Abstract: Relativistic jets are thought to play a crucial role in the formation of massive galaxies and supermassive black holes. Here we report multi-wavelength and multi-epoch observations of the quasar VLASSJ0410-0139 at redshift z=7, powered by a 7e8 solar-mass black hole. Its radio variability, X-ray properties, and compact radio emission on parsec scales reveal that J0410-0139 is a blazar with a relat… ▽ More Relativistic jets are thought to play a crucial role in the formation of massive galaxies and supermassive black holes. Here we report multi-wavelength and multi-epoch observations of the quasar VLASSJ0410-0139 at redshift z=7, powered by a 7e8 solar-mass black hole. Its radio variability, X-ray properties, and compact radio emission on parsec scales reveal that J0410-0139 is a blazar with a relativistic jet aligned with our line of sight. This blazar's existence implies that many more similar (unaligned) jetted sources must exist at z=7. One scenario is that we observe an intrinsically low-power radio jet, but we see it at high luminosity due to relativistic beaming effects. In this case, a large fraction (>80%) of the UV bright quasars must have a similar jet to match the number density expected from the UV quasar luminosity function. These jets can enhance the growth of supermassive black holes and substantially affect their host galaxies. However, the implications would be even more severe if the quasar belongs to the top 10% radio luminous quasars, as measured if the beaming enhancement is less than a factor of 10-15. In this scenario, there should be hundreds to thousands of radio-quiet quasars at z=7 with intrinsic properties similar to J0410-0139 -- in strong tension with the number density of bright quasars derived from their UV luminosity function. To reconcile these results, most black hole growth at z=7 must happen in an obscured phase, as some models predict. The existence of supermassive black holes in the epoch of reionization is facilitated by significant jet-enhanced or obscured super-Eddington accretion. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: Submitted

A quasar-galaxy merger at $z\sim 6.2$: rapid host growth via accretion of two massive satellite galaxies

Authors: Roberto Decarli, Federica Loiacono, Emanuele Paolo Farina, Massimo Dotti, Alessandro Lupi, Romain A. Meyer, Marco Mignoli, Antonio Pensabene, Michael A. Strauss, Bram Venemans, Jinyi Yang, Fabian Walter, Julien Wolf, Eduardo Bañados, Laura Blecha, Sarah Bosman, Chris L. Carilli, Andrea Comastri, Thomas Connor, Tiago Costa, Anna-Christina Eilers, Xiaohui Fan, Roberto Gilli, Hyunsung D. Jun, Weizhe Liu , et al. (16 additional authors not shown)

Abstract: We present JWST/NIRSpec Integral Field Spectroscopy in the rest-frame optical bands of the system PJ308-21, a quasar at $z=6.2342$ caught as its host galaxy interacts with companion galaxies. We detect spatially extended emission of several emission lines (H$α$, H$β$, [OIII], [NII], [SII], HeII), which we use to study the properties of the ionized phase of the interstellar medium: the source and h… ▽ More We present JWST/NIRSpec Integral Field Spectroscopy in the rest-frame optical bands of the system PJ308-21, a quasar at $z=6.2342$ caught as its host galaxy interacts with companion galaxies. We detect spatially extended emission of several emission lines (H$α$, H$β$, [OIII], [NII], [SII], HeII), which we use to study the properties of the ionized phase of the interstellar medium: the source and hardness of the photoionizing radiation field, metallicity, dust reddening, electron density and temperature, and star formation. We also marginally detect continuum starlight emission associated with the companion sources. We find that at least two independent satellite galaxies are part of the system. While the quasar host appears highly enriched and obscured, with AGN-like photoionization conditions, the western companion shows minimal dust extinction, low metallicity ($Z\sim0.4$ Z$_\odot$), and star-formation driven photoionization. The eastern companion shows higher extinction and metallicity ($Z\sim0.8$ Z$_\odot$) compared to the western companion, and it is at least partially photoionized by the nearby quasar. We do not find any indication of AGN in the companion sources. Our study shows that while the quasar host galaxy is already very massive ($M_{\rm dyn}>10^{11}$ M$_\odot$), it is still rapidly building up by accreting two relatively massive ($M_{\rm star}\sim 10^{10}$ M$_\odot$) companion sources. This dataset showcases the power of JWST in exposing the build-up of massive galaxies in the first Gyr of the Universe. △ Less

Submitted 10 June, 2024; originally announced June 2024.

Comments: 15 pages, 16 figures. Accepted for publication in A&A

Deriving the size and shape of the ALBA electron beam with optical synchrotron radiation interferometry using aperture masks: technical choices

Authors: C. L. Carilli, L. Torino, U. Iriso, B. Nikolic, N. Thyagarajan

Abstract: We explore non-redundant aperture masking to derive the size and shape of the ALBA synchrotron light source at optical wavelengths using synchrotron radiation interferometry. We show that non-redundant masks are required due to phase fluctuations arising within the experimental set-up. We also show, using closure phase, that the phase fluctuations are factorizable into element-based errors. We emp… ▽ More We explore non-redundant aperture masking to derive the size and shape of the ALBA synchrotron light source at optical wavelengths using synchrotron radiation interferometry. We show that non-redundant masks are required due to phase fluctuations arising within the experimental set-up. We also show, using closure phase, that the phase fluctuations are factorizable into element-based errors. We employ multiple masks, including 2, 3, 5, and 6 hole configurations. We develop a process for self-calibration of the element-based amplitudes (square root of flux through the aperture), which corrects for non-uniform illumination over the mask, in order to derive visibility coherences and phases, from which the source size and shape can be derived. We explore the optimal procedures to obtain the most reliable results with the 5-hole mask, based on the temporal scatter in measured coherences and closure phases. We find that the closure phases are very stable, and close to zero (within $2^o$). Through uv-modeling, we consider the noise properties of the experiment and conclude that our visibility measurements per frame are likely accurate to an rms scatter of $\sim 1\%$. △ Less

Submitted 4 June, 2024; originally announced June 2024.

Comments: 27 pages, 32 figures, 4 tables. For ALBA Synchrotron Light Source technical memo series

Laboratory Demonstration of Image-Plane Self-Calibration in Interferometry

Authors: Christopher L. Carilli, Bojan Nikolic, Laura Torino, Ubaldo Iriso, Nithyanandan Thyagarajan

Abstract: We demonstrate the Shape-Orientation-Size conservation principle for a 3-element interferometer using aperture plane masking at the ALBA visible synchrotron radiation light source. We then use these data to demonstrate Image Plane Self-Calibration. We demonstrate the Shape-Orientation-Size conservation principle for a 3-element interferometer using aperture plane masking at the ALBA visible synchrotron radiation light source. We then use these data to demonstrate Image Plane Self-Calibration. △ Less

Submitted 20 May, 2024; originally announced May 2024.

Comments: Submitted to the Journal of Optical Society of America A (JOSA A). 8 pages, 9 figures in total

A quasar-galaxy merger at $z\sim6.2$: black hole mass and quasar properties from the NIRSpec spectrum

Authors: Federica Loiacono, Roberto Decarli, Marco Mignoli, Emanuele Paolo Farina, Eduardo Bañados, Sarah Bosman, Anna-Christina Eilers, Jan-Torge Schindler, Michael A. Strauss, Marianne Vestergaard, Feige Wang, Laura Blecha, Chris L. Carilli, Andrea Comastri, Thomas Connor, Tiago Costa, Massimo Dotti, Xiaohui Fan, Roberto Gilli, Hyunsung D. Jun, Weizhe Liu, Alessandro Lupi, Madeline A. Marshall, Chiara Mazzucchelli, Romain A. Meyer , et al. (9 additional authors not shown)

Abstract: We present JWST/NIRSpec integral field data of the quasar PJ308-21 at $z=6.2342$. As shown by previous ALMA and HST imaging, the quasar has two companion sources, interacting with the quasar host galaxy. The high-resolution G395H/290LP NIRSpec spectrum covers the $2.87-5.27\ \rm μm$ wavelength range and shows the rest-frame optical emission of the quasar with exquisite quality ($S/N\sim 100-400$ p… ▽ More We present JWST/NIRSpec integral field data of the quasar PJ308-21 at $z=6.2342$. As shown by previous ALMA and HST imaging, the quasar has two companion sources, interacting with the quasar host galaxy. The high-resolution G395H/290LP NIRSpec spectrum covers the $2.87-5.27\ \rm μm$ wavelength range and shows the rest-frame optical emission of the quasar with exquisite quality ($S/N\sim 100-400$ per spectral element). Based on the H$β$ line from the broad line region, we obtain an estimate of the black hole mass $M_{\rm BH,Hβ}\sim 2.7\times 10^{9}\ \rm M_{\odot}$. This value is within a factor $\lesssim 1.5$ of the H$α$-based black hole mass from the same spectrum ($M_{\rm BH, Hα}\sim 1.93\times 10^{9}\ \rm M_{\odot}$) and is consistent with a previous estimate relying on the MgII $λ2799$ ($M_{\rm BH, MgII}\sim 2.65\times 10^{9}\ \rm M_{\odot}$). All these $M_{\rm BH}$ are within the $\sim 0.5$ dex intrinsic scatter of the adopted mass calibrations. The high Eddington ratio of PJ308-21 $λ_{\rm Edd,Hβ}\sim 0.67$ ($λ_{\rm Edd,Hα}\sim 0.96$) is in line with the overall quasar population at $z \gtrsim 6$. The relative strengths of the [OIII], FeII and H$β$ lines are consistent with the empirical "Eigenvector 1" correlations as observed for low redshift quasars. We find evidence for blueshifted [OIII] $λ5007$ emission with a velocity offset $Δv_{\rm [OIII]}=-1922\pm 39$ km s$^{-1}$ from the systemic velocity and a $\rm FWHM([OIII])=2776^{+75}_{-74}$ km s$^{-1}$. This may be the signature of an outflow from the nuclear region, despite the true values of $Δv_{\rm [OIII]}$ and $\rm FWHM([OIII])$ are likely more uncertain due to the blending with H$β$ and FeII lines. Our study demonstrates the unique capabilities of NIRSpec in capturing quasar spectra at cosmic dawn and studying their properties in unprecedented detail. △ Less

Submitted 20 February, 2024; originally announced February 2024.

Comments: 17 pages, 9 figures, accepted for publication in A&A

Journal ref: A&A 685, A121 (2024)

A demonstration of the effect of fringe-rate filtering in the Hydrogen Epoch of Reionization Array delay power spectrum pipeline

Authors: Hugh Garsden, Philip Bull, Mike Wilensky, Zuhra Abdurashidova, Tyrone Adams, James E. Aguirre, Paul Alexander, Zaki S. Ali, Rushelle Baartman, Yanga Balfour, Adam P. Beardsley, Lindsay M. Berkhout, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Jacob Burba, Steven Carey, Chris L. Carilli, Kai-Feng Chen, Carina Cheng, Samir Choudhuri, David R. DeBoer, Eloy de Lera Acedo, Matt Dexter , et al. (72 additional authors not shown)

Abstract: Radio interferometers targeting the 21cm brightness temperature fluctuations at high redshift are subject to systematic effects that operate over a range of different timescales. These can be isolated by designing appropriate Fourier filters that operate in fringe-rate (FR) space, the Fourier pair of local sidereal time (LST). Applications of FR filtering include separating effects that are correl… ▽ More Radio interferometers targeting the 21cm brightness temperature fluctuations at high redshift are subject to systematic effects that operate over a range of different timescales. These can be isolated by designing appropriate Fourier filters that operate in fringe-rate (FR) space, the Fourier pair of local sidereal time (LST). Applications of FR filtering include separating effects that are correlated with the rotating sky vs. those relative to the ground, down-weighting emission in the primary beam sidelobes, and suppressing noise. FR filtering causes the noise contributions to the visibility data to become correlated in time however, making interpretation of subsequent averaging and error estimation steps more subtle. In this paper, we describe fringe rate filters that are implemented using discrete prolate spheroidal sequences, and designed for two different purposes -- beam sidelobe/horizon suppression (the `mainlobe' filter), and ground-locked systematics removal (the `notch' filter). We apply these to simulated data, and study how their properties affect visibilities and power spectra generated from the simulations. Included is an introduction to fringe-rate filtering and a demonstration of fringe-rate filters applied to simple situations to aid understanding. △ Less

Submitted 13 February, 2024; originally announced February 2024.

Comments: 21 pages, 18 figures, submitted to Monthly Notices of the Royal Astronomical Society

Hydrogen Epoch of Reionization Array (HERA) Phase II Deployment and Commissioning

Authors: Lindsay M. Berkhout, Daniel C. Jacobs, Zuhra 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, Kai-Feng Chen, Carina Cheng, Samir Choudhuri, David R. DeBoer, Eloy de Lera Acedo, Matt Dexter, Joshua S. Dillon , et al. (71 additional authors not shown)

Abstract: This paper presents the design and deployment of the Hydrogen Epoch of Reionization Array (HERA) phase II system. HERA is designed as a staged experiment targeting 21 cm emission measurements of the Epoch of Reionization. First results from the phase I array are published as of early 2022, and deployment of the phase II system is nearing completion. We describe the design of the phase II system an… ▽ More This paper presents the design and deployment of the Hydrogen Epoch of Reionization Array (HERA) phase II system. HERA is designed as a staged experiment targeting 21 cm emission measurements of the Epoch of Reionization. First results from the phase I array are published as of early 2022, and deployment of the phase II system is nearing completion. We describe the design of the phase II system and discuss progress on commissioning and future upgrades. As HERA is a designated Square Kilometer Array (SKA) pathfinder instrument, we also show a number of "case studies" that investigate systematics seen while commissioning the phase II system, which may be of use in the design and operation of future arrays. Common pathologies are likely to manifest in similar ways across instruments, and many of these sources of contamination can be mitigated once the source is identified. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Journal ref: PASP 2024 136 045002

Absorption of Millimeter-band CO and CN in the Early Universe: Molecular Clouds in Radio Galaxy B2 0902+34 at Redshift 3.4

Authors: Bjorn Emonts, Steve Curran, George Miley, Matthew Lehnert, Chris Carilli, Ilsang Yoon, Raffaella Morganti, Reinout van Weeren, Montserrat Villar-Martin, Pierre Guillard, Cristina Cordun, Tom Oosterloo

Abstract: Using the Karl G. Jansky Very Large Array (VLA), we have detected absorption lines due to carbon-monoxide, CO(J=0-1), and the cyano radical, CN(N=0-1), associated with radio galaxy B2 0902+34 at redshift z=3.4. The detection of millimeter-band absorption observed 1.5 Gyr after the Big Bang facilitates studying molecular clouds down to gas masses inaccessible to emission-line observations. The CO a… ▽ More Using the Karl G. Jansky Very Large Array (VLA), we have detected absorption lines due to carbon-monoxide, CO(J=0-1), and the cyano radical, CN(N=0-1), associated with radio galaxy B2 0902+34 at redshift z=3.4. The detection of millimeter-band absorption observed 1.5 Gyr after the Big Bang facilitates studying molecular clouds down to gas masses inaccessible to emission-line observations. The CO absorption in B2 0902+34 has a peak optical depth of $τ$ $\ge$ 8.6% and consists of two components, one of which has the same redshift as previously detected 21-cm absorption of neutral hydrogen (HI) gas. Each CO component traces an integrated H$_2$ column density of N(H2) $\ge$ 3x10$^{20}$ cm$^{-2}$. CN absorption is detected for both CO components, as well as for a blueshifted component not detected in CO, with CO/CN line ratios ranging from $\le$0.4 to 2.4. We discuss the scenario that the absorption components originate from collections of small and dense molecular clouds that are embedded in a region with more diffuse gas and high turbulence, possibly within the influence of the central Active Galactic Nucleus or starburst region. The degree of reddening in B2 0902+34, with a rest-frame color B-K ~ 4.2, is lower than the very red colors (B-K > 6) found among other known redshifted CO absorption systems at z<1. Nevertheless, when including also the many non-detections from the literature, a potential correlation between the absorption-line strength and B-K color is evident, giving weight to the argument that the red colors of CO absorbers are due to a high dust content. △ Less

Submitted 7 January, 2024; originally announced January 2024.

Comments: Accepted by ApJ (submitted 25-Sept-2023, accepted 27-Dec-2023)

matvis: A matrix-based visibility simulator for fast forward modelling of many-element 21 cm arrays

Authors: Piyanat Kittiwisit, Steven G. Murray, Hugh Garsden, Philip Bull, Christopher Cain, Aaron R. Parsons, Jackson Sipple, Zara Abdurashidova, Tyrone Adams, James E. Aguirre, Paul Alexander, Zaki S. Ali, Rushelle Baartman, Yanga Balfour, Adam P. Beardsley, Lindsay M. Berkhout, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Jacob Burba, Steven Carey, Chris L. Carilli, Kai-Feng Chen, Carina Cheng , et al. (73 additional authors not shown)

Abstract: Detection of the faint 21 cm line emission from the Cosmic Dawn and Epoch of Reionisation will require not only exquisite control over instrumental calibration and systematics to achieve the necessary dynamic range of observations but also validation of analysis techniques to demonstrate their statistical properties and signal loss characteristics. A key ingredient in achieving this is the ability… ▽ More Detection of the faint 21 cm line emission from the Cosmic Dawn and Epoch of Reionisation will require not only exquisite control over instrumental calibration and systematics to achieve the necessary dynamic range of observations but also validation of analysis techniques to demonstrate their statistical properties and signal loss characteristics. A key ingredient in achieving this is the ability to perform high-fidelity simulations of the kinds of data that are produced by the large, many-element, radio interferometric arrays that have been purpose-built for these studies. The large scale of these arrays presents a computational challenge, as one must simulate a detailed sky and instrumental model across many hundreds of frequency channels, thousands of time samples, and tens of thousands of baselines for arrays with hundreds of antennas. In this paper, we present a fast matrix-based method for simulating radio interferometric measurements (visibilities) at the necessary scale. We achieve this through judicious use of primary beam interpolation, fast approximations for coordinate transforms, and a vectorised outer product to expand per-antenna quantities to per-baseline visibilities, coupled with standard parallelisation techniques. We validate the results of this method, implemented in the publicly-available matvis code, against a high-precision reference simulator, and explore its computational scaling on a variety of problems. △ Less

Submitted 15 December, 2023; originally announced December 2023.

Comments: 25 pages, 20 figures, submitted to RAS Techniques and Instruments, matvis is publicly available at https://github.com/HERA-Team/matvis

Feeding and feedback processes in the Spiderweb proto-intracluster medium

Authors: M. Lepore, L. Di Mascolo, P. Tozzi, E. Churazov, T. Mroczkowski, S. Borgani, C. Carilli, M. Gaspari, M. Ginolfi, A. Liu, L. Pentericci, E. Rasia, P. Rosati, H. J. A. Röttgering, C. S. Anderson, H. Dannerbauer, G. Miley, C. Norman

Abstract: We present the detailed analysis of the thermal, diffuse emission of the proto-intracluster medium (ICM) detected in the halo of the Spiderweb Galaxy at z=2.16, within a radius of $\sim$ 150 kpc. We combined deep X-ray data from Chandra and millimeter observations of the Sunyaev-Zeldovich (SZ) effect obtained by ALMA. Thanks to independent measurements of the pressure profile from ALMA SZ observat… ▽ More We present the detailed analysis of the thermal, diffuse emission of the proto-intracluster medium (ICM) detected in the halo of the Spiderweb Galaxy at z=2.16, within a radius of $\sim$ 150 kpc. We combined deep X-ray data from Chandra and millimeter observations of the Sunyaev-Zeldovich (SZ) effect obtained by ALMA. Thanks to independent measurements of the pressure profile from ALMA SZ observation and the electron density profile from the available X-ray data, we derived, for the first time, the temperature profile in the ICM of a z>2 protocluster. It reveals the presence of a strong cool core (comparable to the local ones) that may host a significant mass deposition flow, consistent with measured local star formation values. We also find mild evidence of an asymmetry in the X-ray surface brightness distribution, which may be tentatively associated with a cavity carved into the proto-ICM by the radio jets or, alternatively, may be due to the young dynamical status of the halo. The cooling time of baryons in the core of the Spiderweb Protocluster is estimated to be $\sim$ 0.1 Gyr, implying that the baryon cycle in the first stages of the protocluster formation is characterised by a high-duty cycle and a very active environment. In the case of the Spiderweb protocluster, we are witnessing the presence of a strongly peaked core that is possibily hosting a cooling flow with a mass deposition rate up to 250-1000 $M_{\odot}$/yr, responsible for feeding both the central supermassive black hole and the high star formation rate observed in the Spiderweb Galaxy. This phase is expected to be rapidly followed by active galactic nucleus feedback events, whose onset may have already left an imprint in the radio and X-ray appearance of the Spiderweb protocluster, eventually driving the ICM into a self-regulated, long-term evolution in less than one Gyr. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Bayesian estimation of cross-coupling and reflection systematics in 21cm array visibility data

Authors: Geoff G. Murphy, Philip Bull, Mario G. Santos, Zara Abdurashidova, Tyrone Adams, James E. Aguirre, Paul Alexander, Zaki S. Ali, Rushelle Baartman, Yanga Balfour, Adam P. Beardsley, Gianni Bernardi, Tashalee Billings, Judd D. Bowman, Richard F. Bradley, Jacob Burba, Christopher Cain, Steven Carey, Chris L. Carilli, Carina Cheng, David R. DeBoer, Eloy de Lera Acedo, Matt Dexter, Joshua S. Dillon, Nico Eksteen , et al. (54 additional authors not shown)

Abstract: Observations with radio arrays that target the 21-cm signal originating from the early Universe suffer from a variety of systematic effects. An important class of these are reflections and spurious couplings between antennas. We apply a Hamiltonian Monte Carlo sampler to the modelling and mitigation of these systematics in simulated Hydrogen Epoch of Reionisation Array (HERA) data. This method all… ▽ More Observations with radio arrays that target the 21-cm signal originating from the early Universe suffer from a variety of systematic effects. An important class of these are reflections and spurious couplings between antennas. We apply a Hamiltonian Monte Carlo sampler to the modelling and mitigation of these systematics in simulated Hydrogen Epoch of Reionisation Array (HERA) data. This method allows us to form statistical uncertainty estimates for both our models and the recovered visibilities, which is an important ingredient in establishing robust upper limits on the Epoch of Reionisation (EoR) power spectrum. In cases where the noise is large compared to the EoR signal, this approach can constrain the systematics well enough to mitigate them down to the noise level for both systematics studied. Where the noise is smaller than the EoR, our modelling can mitigate the majority of the reflections with there being only a minor level of residual systematics, while cross-coupling sees essentially complete mitigation. Our approach performs similarly to existing filtering/fitting techniques used in the HERA pipeline, but with the added benefit of rigorously propagating uncertainties. In all cases it does not significantly attenuate the underlying signal. △ Less

Submitted 6 December, 2023; originally announced December 2023.

Comments: 19 pages, 14 figures, submitted to MNRAS

Direct Optimal Mapping Image Power Spectrum and its Window Functions

Authors: Zhilei Xu, Honggeun Kim, Jacqueline N. Hewitt, Kai-Feng Chen, Nicholas S. Kern, Eleanor Rath, Ruby Byrne, Adélie Gorce, Robert Pascua, Zachary E. Martinot, Joshua S. Dillon, Bryna J. Hazelton, Adrian Liu, Miguel F. Morales, 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 , et al. (57 additional authors not shown)

Abstract: The key to detecting neutral hydrogen during the epoch of reionization (EoR) is to separate the cosmological signal from the dominating foreground radiation. We developed direct optimal mapping (DOM) to map interferometric visibilities; it contains only linear operations, with full knowledge of point spread functions from visibilities to images. Here, we demonstrate a fast Fourier transform-based… ▽ More The key to detecting neutral hydrogen during the epoch of reionization (EoR) is to separate the cosmological signal from the dominating foreground radiation. We developed direct optimal mapping (DOM) to map interferometric visibilities; it contains only linear operations, with full knowledge of point spread functions from visibilities to images. Here, we demonstrate a fast Fourier transform-based image power spectrum and its window functions computed from the DOM images. We use noiseless simulation, based on the Hydrogen Epoch of Reionization Array Phase I configuration, to study the image power spectrum properties. The window functions show $<10^{-11}$ of the integrated power leaks from the foreground-dominated region into the EoR window; the 2D and 1D power spectra also verify the separation between the foregrounds and the EoR. △ Less

Submitted 5 July, 2024; v1 submitted 17 November, 2023; originally announced November 2023.

Comments: Published in ApJ

Imaging the ring and jet in M87 at 85 GHz with the ngEHT and ngEHT+ngVL

Authors: C. L. Carilli, R. C. Walker, E. Murphy, B. Mason

Abstract: The Global mm-VLBI Array (GMVA) has demonstrated the ability to resolved what may be the general relativistic shadow of the supermassive black hole in M87 at 86 GHz, as well as delineate the inner jet to $\sim 1$~mas distance. We investigate the ability of the planned ngEHT, and the ngEHT + ngVLA at 85 GHz, to image such a nuclear 'ring', and the associated jet, using a constructed model based on… ▽ More The Global mm-VLBI Array (GMVA) has demonstrated the ability to resolved what may be the general relativistic shadow of the supermassive black hole in M87 at 86 GHz, as well as delineate the inner jet to $\sim 1$~mas distance. We investigate the ability of the planned ngEHT, and the ngEHT + ngVLA at 85 GHz, to image such a nuclear 'ring', and the associated jet, using a constructed model based on the current estimate of the ring size, and a scaled version of the best VLBA image of the M87 jet at 43 GHz. While the resolution does not improve due to the limit set by the diameter of the Earth, the ngEHT alone should provide both a higher fidelity image of the ring on scales $\le 0.1$~mas, and a good image of a more extended jet to $\sim 1$~mas. Adding the ngVLA improves substantially the dynamic range (factor 3.5), as well as adds the ability to image structures on larger scales, in this case out to at least 5~mas, and potentially to much larger scales given the $\sim 10^5$ range in spatial scales covered by the ngVLA itself. Both arrays provide good image fidelity ($\le 0.1$), in the inner $\sim 1$~mas, but the ngEHT-only image does not reproduce the outer jet well, or at all, with fidelity values greater than unity. The combined array reproduces much of the outer jet with good fidelity ($\le 0.3$). Adding the ngVLA also decreases the susceptibility to antenna-based phase errors by a similar factor, and improves the ability for fringe fitting and subsequent phase and amplitude self-calibration. As for scales $< 100~μ$as, ie. the ring itself, adding the ngVLA makes little change for very bright sources, where uniform weighting can be employed. But for faint sources, adding the ngVLA adds potentially an order-of magnitude sensitivity improvement (Issaoun et al. 2023). △ Less

Submitted 22 September, 2023; originally announced September 2023.

Comments: 14 pages, 5 figures, ngVLA memo 118

NOEMA observations of GN-z11: Constraining Neutral Interstellar Medium and Dust Formation in the Heart of Cosmic Reionization at $z=10.6$

Authors: Y. Fudamoto, P. A. Oesch, F. Walter, R. Decarli, C. L. Carilli, A. Ferrara, L. Barrufet, R. Bouwens, M. Dessauges-Zavadsky, E. J. Nelson, H. Dannerbauer, G. Illingworth, A. K. Inoue, R. Marques-Chaves, I. Pérez-Fournon, D. A. Riechers, D. Schaerer, R. Smit, Y. Sugahara, P. van der Werf

Abstract: We present results of dust continuum and [CII]$\,158\,{\rm μm}$ emission line observations of a remarkably UV-luminous ($M_{\rm UV}=-21.6$) galaxy at $z=10.603$: GN-z11. Using the Northern Extended Millimeter Array (NOEMA), observations have been carried out over multiple observing cycles. We achieved a high sensitivity resulting in a $λ_{\rm rest}=160\,{\rm μm}$ continuum $1\,σ$ depth of… ▽ More We present results of dust continuum and [CII]$\,158\,{\rm μm}$ emission line observations of a remarkably UV-luminous ($M_{\rm UV}=-21.6$) galaxy at $z=10.603$: GN-z11. Using the Northern Extended Millimeter Array (NOEMA), observations have been carried out over multiple observing cycles. We achieved a high sensitivity resulting in a $λ_{\rm rest}=160\,{\rm μm}$ continuum $1\,σ$ depth of $13.0\,\rm{μJy/beam}$ and a [CII] emission line $1\,σ$ sensitivity of $31\,\rm{mJy/beam\,km/s}$ using $50\,\rm{km/s}$ binning with a $\sim 2\,{\rm arcsec}$ synthesized beam. Neither dust continuum nor [CII]$\,158\,{\rm μm}$ line emission are detected at the expected frequency of $ν_{\rm [CII]} = 163.791\,\rm{GHz}$ and the sky location of GN-z11. The upper limits show that GN-z11 is neither luminous in $L_{\rm IR}$ nor $L_{\rm [CII]}$, with a dust mass $3\,σ$ limit of ${\rm log}(M_{\rm dust}/{\rm M_{\odot}}) < 6.5-6.9$ and with a [CII] based molecular gas mass $3\,σ$ limit of ${\rm log}(M_{\rm mol,[CII]}/{\rm M_{\odot}}) < 9.3$. Together with radiative transfer calculations, we also investigated the possible cause of the dust poor nature of the GN-z11 showed by the blue color in the UV continuum of GN-z11 ($β_{\rm UV}=-2.4$), and found that $\gtrsim3\times$ deeper observations are crucial to study dust production at very high-redshift. Nevertheless, our observations show the crucial role of deep mm/submm observations of very high redshift galaxies to constrain multiple phases in the interstellar medium. △ Less

Submitted 5 September, 2023; originally announced September 2023.

Comments: submitted to MNRAS, 7 pages, 4 figures, 1 table

CO survey of high-z radio galaxies, revisited with ALMA: Jet-cloud Alignments and Synchrotron Brightening by Molecular Gas in the Circumgalactic Environment

Authors: Bjorn Emonts, Matthew Lehnert, Sophie Lebowitz, George K. Miley, Montserrat Villar-Martin, Ray Norris, Carlos De Breuck, Chris Carilli, Ilana Feain

Abstract: Powerful radio sources associated with super-massive black holes are among the most luminous objects in the Universe, and are frequently recognized both as cosmological probes and active constituents in the evolution of galaxies. We present alignments between radio jets and cold molecular gas in the environment of distant radio galaxies, and show that the brightness of the radio synchrotron source… ▽ More Powerful radio sources associated with super-massive black holes are among the most luminous objects in the Universe, and are frequently recognized both as cosmological probes and active constituents in the evolution of galaxies. We present alignments between radio jets and cold molecular gas in the environment of distant radio galaxies, and show that the brightness of the radio synchrotron source can be enhanced by its interplay with the molecular gas. Our work is based on CO J>1 observations with the Atacama Large Millimeter/submillimeter Array (ALMA) of three radio galaxies with redshifts in the range 1.4 < z < 2.1, namely MRC 0114-211 (z = 1.41), MRC 0156-252 (z = 2.02), and MRC 2048-272 (z = 2.05). These ALMA observations support previous work that found molecular gas out to 50 kpc in the circumgalactic environment, based on a CO(1-0) survey performed with the Australia Telescope Compact Array (ATCA). The CO emission is found along the radio axes but beyond the main radio lobes. When compared to a large sample of high-z radio galaxies from the literature, we find that the presence of this cold molecular medium correlates with an increased flux-density ratio of the main vs. counter lobe. This suggest that the radio lobe brightens when encountering cold molecular gas in the environment. While part of the molecular gas is likely related to the interstellar medium (ISM) from either the host or a companion galaxy, a significant fraction of the molecular gas in these systems shows very low excitation, with r$_{2-1/1-0}$ and r$_{3-2/1-0}$ values $\lesssim$0.2. This could be part of the circumgalactic medium (CGM). △ Less

Submitted 21 June, 2023; originally announced June 2023.

Comments: Accepted for publication in ApJ (19 pages, 6 figures)

Forming intracluster gas in a galaxy protocluster at a redshift of 2.16

Authors: Luca Di Mascolo, Alexandro Saro, Tony Mroczkowski, Stefano Borgani, Eugene Churazov, Elena Rasia, Paolo Tozzi, Helmut Dannerbauer, Kaustuv Basu, Christopher L. Carilli, Michele Ginolfi, George Miley, Mario Nonino, Maurilio Pannella Laura Pentericci, Francesca Rizzo

Abstract: Galaxy clusters are the most massive gravitationally bound structures in the Universe, comprising thousands of galaxies and pervaded by a diffuse, hot ``intracluster medium'' (ICM) that dominates the baryonic content of these systems. The formation and evolution of the ICM across cosmic time is thought to be driven by the continuous accretion of matter from the large-scale filamentary surroundings… ▽ More Galaxy clusters are the most massive gravitationally bound structures in the Universe, comprising thousands of galaxies and pervaded by a diffuse, hot ``intracluster medium'' (ICM) that dominates the baryonic content of these systems. The formation and evolution of the ICM across cosmic time is thought to be driven by the continuous accretion of matter from the large-scale filamentary surroundings and dramatic merger events with other clusters or groups. Until now, however, direct observations of the intracluster gas have been limited only to mature clusters in the latter three-quarters of the history of the Universe, and we have been lacking a direct view of the hot, thermalized cluster atmosphere at the epoch when the first massive clusters formed. Here we report the detection (about $6σ$) of the thermal Sunyaev-Zeldovich (SZ) effect in the direction of a protocluster. In fact, the SZ signal reveals the ICM thermal energy in a way that is insensitive to cosmological dimming, making it ideal for tracing the thermal history of cosmic structures. This result indicates the presence of a nascent ICM within the Spiderweb protocluster at redshift $z=2.156$, around 10 billion years ago. The amplitude and morphology of the detected signal show that the SZ effect from the protocluster is lower than expected from dynamical considerations and comparable with that of lower-redshift group-scale systems, consistent with expectations for a dynamically active progenitor of a local galaxy cluster. △ Less

Submitted 28 March, 2023; originally announced March 2023.

Comments: 18 pages (main text + methods), 10 figures, 2 tables; published online in Nature on March 29th, 2023

Search for the Epoch of Reionisation with HERA: Upper Limits on the Closure Phase Delay Power Spectrum

Authors: 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, Yanga Balfour, Adam P. Beardsley, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley , et al. (58 additional authors not shown)

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 standa… ▽ More 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. △ Less

Submitted 15 February, 2023; originally announced February 2023.

Comments: 16 pages, 14 figures, accepted for publication by MNRAS

A NOEMA molecular line scan of the Hubble Deep Field North: Improved constraints on the CO luminosity functions and cosmic density of molecular gas

Authors: Leindert A. Boogaard, Roberto Decarli, Fabian Walter, Axel Weiss, Gergö Popping, Roberto Neri, Manuel Aravena, Dominik Riechers, Richard S. Ellis, Chris Carilli, Pierre Cox, Jérôme Pety

Abstract: We present measurements of the CO luminosity functions (LFs) and the evolution of the cosmic molecular gas density out to z~6 based on an 8.5 arcmin^2 spectral scan survey at 3mm of the iconic Hubble Deep Field North (HDF-N) observed with the NOrthern Extended Millimeter Array (NOEMA). We use matched filtering to search for line emission from galaxies and determine their redshift probability distr… ▽ More We present measurements of the CO luminosity functions (LFs) and the evolution of the cosmic molecular gas density out to z~6 based on an 8.5 arcmin^2 spectral scan survey at 3mm of the iconic Hubble Deep Field North (HDF-N) observed with the NOrthern Extended Millimeter Array (NOEMA). We use matched filtering to search for line emission from galaxies and determine their redshift probability distributions exploiting the extensive multi-wavelength data for the HDF-N. We identify the 7 highest-fidelity sources as CO emitters at 1<z<6, including the well-known submillimeter galaxy HDF850.1 at z=5.18. Four high-fidelity 3mm continuum sources are all found to be radio galaxies at z<=1, plus HDF850.1. We constrain the CO LFs in the HDF-N out to z~6, including a first measurement of the CO(5-4) LF at <z>=5.0. The relatively large area and depth of the NOEMA HDF-N survey extends the existing luminosity functions at 1<z<4 above the knee, yielding a somewhat lower density by 0.15-0.4 dex at the overlap region for the CO(2-1) and CO(3-2) transitions, attributed to cosmic variance. We perform a joint analysis of the CO LFs in the HDF-N and Hubble Ultra Deep Field from ASPECS, finding that they can be well described by a single Schechter function. The evolution of the cosmic molecular gas density from a joint analysis is in good agreement with earlier determinations. This implies that the impact of cosmic field-to-field variance on the measurements is consistent with previous estimates, adding to the challenges for simulations that model galaxies from first principles. △ Less

Submitted 13 January, 2023; originally announced January 2023.

Comments: 19 pages, 13 figures, 7 tables. Accepted for publication in the Astrophysical Journal (ApJ)

ALMA FIR View of Ultra High-redshift Galaxy Candidates at $z\sim$ 11-17: Blue Monsters or Low-$z$ Red Interlopers?

Authors: 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, Takashi Kojima, Rebecca L. Larson, Eric J. Murphy, Casey Papovich, Pablo G. Pérez-González , et al. (28 additional authors not shown)

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~$μ$m, ruling out the possibility that \targb\ is a low-$z$ dusty starburst with a star-formation rate of $\gtrsim 30$~$M_{\odot}$~yr… ▽ More 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~$μ$m, ruling out the possibility that \targb\ is a low-$z$ dusty starburst with a star-formation rate of $\gtrsim 30$~$M_{\odot}$~yr$^{-1}$. We detect a 5.1$σ$ 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$μ$m at $z=16.01$ or [CII]158$μ$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 $\simeq-2.3$) and luminous ($M_{\rm UV}\simeq[-$24:$-21]$) system. Still, we find in several candidates that reasonable SED fits ($Δ$ $χ^{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$β$) distributions are abundant in survey volumes used for the identification of 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 several of the $z\sim$11--17 candidates, indicating the need to consider the relative surface densities of the lower-$z$ contaminants in the ultra high-$z$ galaxy search. △ Less

Submitted 26 July, 2023; v1 submitted 7 November, 2022; originally announced November 2022.

Comments: 25 pages, 8 figures, 7 tables. ApJ in press

Image-Plane Self-Calibration in Interferometry

Authors: C. L. Carilli, B. Nikolic, N. Thyagarajan

Abstract: We develop a new process of image plane self-calibration for interferometric imaging data. The process is based on Shape-Orientation-Size (SOS) conservation for the principal triangle in an image generated from the three fringes made from a triad of receiving elements, in situations where interferometric phase errors can be factorized into element-based terms. The basis of the SOS conservation pri… ▽ More We develop a new process of image plane self-calibration for interferometric imaging data. The process is based on Shape-Orientation-Size (SOS) conservation for the principal triangle in an image generated from the three fringes made from a triad of receiving elements, in situations where interferometric phase errors can be factorized into element-based terms. The basis of the SOS conservation principle is that, for a 3-element array, the only possible image corruption due to an element-based phase screen is a tilt of the aperture plane, leading to a shift in the image plane. Thus, an image made from any 3-element interferometer represents a true image of the source brightness, modulo an unknown translation. Image plane self-calibration entails deriving the unknown translations for each triad image via cross-correlation of the observed triad image with a model image of the source brightness. After correcting for these independent shifts, and summing the aligned triad images, a good image of the source brightness is generated from the full array, recovering source structure at diffraction-limited resolution. The process is iterative, using improved source models based on previous iterations. We demonstrate the technique in the high signal-to-noise context, and include a configuration based on radio astronomical facilities, and simple models of double sources. We show that the process converges for the simple models considered, although convergence is slower than for aperture-plane self-calibration for large-$N$ arrays. As currently implemented, the process is most relevant for arrays with a small number of elements. More generally, the technique provides geometric insight into closure phase and the self-calibration process. The technique is generalizable to non-astronomical interferometric imaging applications across the electromagnetic spectrum. △ Less

Submitted 27 October, 2022; originally announced October 2022.

Comments: 17 pages 4 figures. To appear in Journal Optical Society America - A

Characterization Of Inpaint Residuals In Interferometric Measurements of the Epoch Of Reionization

Authors: 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, Jacob Burba, Steven Carey, Chris L. Carilli, Carina Cheng, David R. DeBoer , et al. (53 additional authors not shown)

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 du… ▽ More 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. △ Less

Submitted 20 February, 2023; v1 submitted 26 October, 2022; originally announced October 2022.

Comments: 21 pages, 13 figures

Spectroscopy of High Redshift Galaxies with the ngVLA

Authors: C. L. Carilli, M. Neeleman

Abstract: We present simulations of the capabilities of the ngVLA to image at $\sim 0.75$ kpc resolution ($0.085"$), molecular line emission from star forming disk galaxies at high redshift. The results are compared to the current capabilities of ALMA. ALMA can detect the integrated emission, and determine the velocity gradient and size across the brighter emission regions of the galaxy. The ngVLA is a fact… ▽ More We present simulations of the capabilities of the ngVLA to image at $\sim 0.75$ kpc resolution ($0.085"$), molecular line emission from star forming disk galaxies at high redshift. The results are compared to the current capabilities of ALMA. ALMA can detect the integrated emission, and determine the velocity gradient and size across the brighter emission regions of the galaxy. The ngVLA is a factor $\sim 6$ more sensitive at the adopted spatial and velocity resolution. This sensitivity is needed to recover the detailed column density distribution, velocity field, and velocity dispersion at full resolution. The ngVLA will enable detailed analysis of spectral line profiles at $0.75$~kpc resolution, even in relatively faint regions. The ngVLA will trace the rotation curves to large radii, and recover sub-structure in the disks, such as clumps, spiral arms, bars, and rings. Detection of these features is crucial in order to assess how cold gas precipitates the formation of stars at high redshift. △ Less

Submitted 24 October, 2022; originally announced October 2022.

Comments: ngVLA memo series 103 https://ngvla.nrao.edu/page/memos

The interstellar medium distribution, gas kinematics, and system dynamics of the far-infrared luminous quasar SDSS J2310+1855 at $z=6.0$

Authors: 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 compon… ▽ More 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.] △ Less

Submitted 21 October, 2022; originally announced October 2022.

Comments: Accepted for publication in A&A

Journal ref: A&A 668, A121 (2022)

ALMA Observation of a $z\gtrsim10$ Galaxy Candidate Discovered with JWST

Authors: Ilsang Yoon, Christopher L. Carilli, Seiji Fujimoto, Marco Castellano, Emiliano Merlin, Paola Santini, Min S. Yun, Eric J. Murphy, Intae Jung, Caitlin M. Casey, Steven L. Finkelstein, Casey Papovich, Adriano Fontana, Tommaso Treu, Jonathan Letai

Abstract: We report the ALMA observation of a $z\gtrsim10$ galaxy candidate (GHZ1) discovered from the GLASS-JWST Early Release Science Program. Our ALMA program aims to detect the [OIII] emission line at the rest-frame 3393.0062 GHz ($88.36μ$m) and far-IR continuum emission with the spectral window setup seamlessly covering a 26.125 GHz frequency range ($10.10<z<11.14$). A total of 7 hours of on-source int… ▽ More We report the ALMA observation of a $z\gtrsim10$ galaxy candidate (GHZ1) discovered from the GLASS-JWST Early Release Science Program. Our ALMA program aims to detect the [OIII] emission line at the rest-frame 3393.0062 GHz ($88.36μ$m) and far-IR continuum emission with the spectral window setup seamlessly covering a 26.125 GHz frequency range ($10.10<z<11.14$). A total of 7 hours of on-source integration was employed, using four frequency settings to cover the full range (1.7 hours per setting), with $0''.7$ angular resolution. No line or continuum is clearly detected, with a 5$σ$ upper limit of the line emission of 0.93 mJy beam$^{-1}$ at 25 km s$^{-1}$ channel$^{-1}$ and of the continuum emission of 30$μ$Jy beam$^{-1}$. We report marginal spectral (at 225 km s$^{-1}$ resolution) and continuum features ($4.1σ$ and $2.6σ$ peak signal-to-noise ratio, respectively), within $0''.17$ from the JWST position of GHZ1. This spectral feature implies $z=10.38$ and needs to be verified with further observations. Assuming that the best photometric redshift estimate ($z=10.60^{+0.52}_{-0.60}$) is correct, the broadband galaxy spectral energy distribution model for the $3σ$ upper limit of the continuum flux from GHZ1 suggests that GHZ1 has a small amount of dust ($M_d\lesssim10^4 M_{\odot}$) with high temperature ($T_d\gtrsim90$K). The $5σ$ upper limit of the [OIII]$_{88μm}$ line luminosity and the inferred star formation rate of GHZ1 is consistent with the properties of the low metallicity dwarf galaxies. We also report serendipitous clear detections of six continuum sources at the locations of the JWST galaxy counterparts in the field. △ Less

Submitted 29 March, 2023; v1 submitted 15 October, 2022; originally announced October 2022.

Comments: Accepted to ApJ after revising the figures and the analysis

Improved Constraints on the 21 cm EoR Power Spectrum and the X-Ray Heating of the IGM with HERA Phase I Observations

Authors: 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, Samir Choudhuri, David R. DeBoer, Eloy de Lera Acedo, Matt Dexter, Joshua S. Dillon , et al. (70 additional authors not shown)

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 $Δ^2(k = 0.34$ $h$ Mpc$^{-1}$) $\leq 457$ mK$^2$ at $z = 7.9$ and that… ▽ More 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 $Δ^2(k = 0.34$ $h$ Mpc$^{-1}$) $\leq 457$ mK$^2$ at $z = 7.9$ and that $Δ^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. △ Less

Submitted 19 January, 2023; v1 submitted 10 October, 2022; originally announced October 2022.

Comments: 57 pages, 37 figures. Updated to match the accepted ApJ version. Corresponding author: Joshua S. Dillon

Journal ref: 2023 ApJ 945 124

Impact of instrument and data characteristics in the interferometric reconstruction of the 21 cm power spectrum

Authors: 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, Carina Cheng, David R. DeBoer, Eloy de Lera Acedo, Matt Dexter, Joshua S. Dillon , et al. (53 additional authors not shown)

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 t… ▽ More 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. △ Less

Submitted 11 January, 2023; v1 submitted 7 October, 2022; originally announced October 2022.

Comments: 18 pages, 19 figures, accepted for publication in MNRAS

The 700 ks Chandra Spiderweb Field II: Evidence for inverse-Compton and thermal diffuse emission in the Spiderweb galaxy

Authors: P. Tozzi, R. Gilli, A. Liu, S. Borgani, M. Lepore, L. Di Mascolo, A. Saro, L. Pentericci, C. Carilli, G. Miley, T. Mroczkowski, M. Pannella, E. Rasia, P. Rosati, C. S. Anderson, A. Calabro', E. Churazov, H. Dannerbauer, C. Feruglio, F. Fiore, R. Gobat, S. Jin, M. Nonino, C. Norman, H. J. A. Rottgering

Abstract: We present the X-ray imaging and spectral analysis of the diffuse emission around the Spiderweb galaxy at z=2.16 and of its nuclear emission, based on a deep (700 ks) Chandra observation. We characterize the nuclear emission and computed the contamination in the surrounding regions due to the wings of the instrument PSF. Then, we quantified the extended emission within 12". We find that the Spider… ▽ More We present the X-ray imaging and spectral analysis of the diffuse emission around the Spiderweb galaxy at z=2.16 and of its nuclear emission, based on a deep (700 ks) Chandra observation. We characterize the nuclear emission and computed the contamination in the surrounding regions due to the wings of the instrument PSF. Then, we quantified the extended emission within 12". We find that the Spiderweb galaxy hosts a mildly absorbed quasar, with modest yet significant variability on a timescale of ~1 year. We find that the emission in the jet regions is well described by a power law with Gamma~2-2.5, and it is consistent with IC upscattering of the CMB photons by the relativistic electrons. We also find a roughly symmetric, diffuse emission within a radius of ~100 kpc. This emission is consistent with thermal bremsstrahlung from a hot ICM with a temperature of kT=2.0_{-0.4}^{+0.7} keV, and a metallicity of Z<1.6Z_sun. The average electron density within 100 kpc is n_e=(1.51+-0.24+-0.14)E-2 cm^{-3}, corresponding to an upper limit for the total ICM mass of <=(1.76+-0.30+-0.17)E+12 M_sun (where error bars are 1 sigma statistical and systematic, respectively). If we apply hydrostatic equilibrium to the ICM, we measure a total gravitational mass M(<100 kpc)=(1.5^{+0.5}_{-0.3})E+13 M_sun and, extrapolating at larger radii, we estimate a total mass M_{500}=(3.2^{+1.1}_{-0.6})E+13 M_sun within a radius of r_{500}=(220+-30) kpc. We conclude that the Spiderweb protocluster shows significant diffuse emission within a radius of 12 arcsec, whose major contribution is provided by IC scattering associated with the radio jets. Outside the jet regions, we also identified thermal emission within a radius of ~100 kpc, revealing the presence of hot, diffuse baryons that may represent the embryonic virialized halo of the forming cluster. △ Less

Submitted 30 September, 2022; originally announced September 2022.

Comments: 19 pages, 16 figures, accepted in Astronomy & Astrophysics

A 4~Kpc Molecular Gas Lane in Cygnus A

Authors: C. L. Carilli, R. A. Perley, D. A. Perley, V. Dhawan, R. Decarli, A. Evans, K. Nyland

Abstract: We present the discovery of a 4 kpc molecular gas lane in the Cygnus A host galaxy, using ALMA CO 2-1 observations. The gas lane is oriented roughly perpendicular to the projected radio jet axis. The CO emission generally follows the clumpy dust lanes seen in HST I-band images. The total molecular gas mass is $30\times 10^8$ M$_\odot$ for Milky Way type clouds, and $3.6 \times 10^8$ M$_\odot$ for… ▽ More We present the discovery of a 4 kpc molecular gas lane in the Cygnus A host galaxy, using ALMA CO 2-1 observations. The gas lane is oriented roughly perpendicular to the projected radio jet axis. The CO emission generally follows the clumpy dust lanes seen in HST I-band images. The total molecular gas mass is $30\times 10^8$ M$_\odot$ for Milky Way type clouds, and $3.6 \times 10^8$ M$_\odot$ for starburst conditions. There is a velocity change from the northern to southern CO peaks of about $\pm 175$~km~s$^{-1}$, and an apparently smooth velocity gradient between the peaks, although the emission in the central region is weak. In the inner $\sim 0.5"$ projected distance from the radio core, comparison of the CO velocities to those observed for H$_2$ 2.1218 $μ$m emission shows higher velocities for the vibrationally excited warm molecular gas than the cooler CO 2-1 line emitting gas at similar projected radii. A possible explanation for these different projected velocities at a given radius is that the cooler CO gas is distributed in a clumpy ring at radius $\sim 1.5"$ to $2"$, while the warm H$_2$ 2.12$μ$m emitting gas is interior to this ring. Of course, the current data cannot rule-out a clumpy, amorphous molecular gas distribution linearly distributed perpendicular to the radio jet axis. We consider surface brightness properties on scales down to $\sim 265$~pc, and discuss the Cygnus A results in the context of other radio galaxies with CO emission. △ Less

Submitted 2 September, 2022; originally announced September 2022.

Comments: 5 Figures, 23 pages, 2 tables, accepted for publication in the ApJ

The Spiderweb proto-cluster is being magnetized by its central radio jet

Authors: Craig S. Anderson, Christopher L. Carilli, Paolo Tozzi, G. K. Miley, S. Borgani, Tracy Clarke, L. Di Mascolo, Ang Liu, Tony Mroczkowski, Maurilio Pannella, L Pentericci, H. J. A. Rottgering, A. Saro

Abstract: We present deep broadband radio polarization observations of the Spiderweb radio galaxy (J1140-2629) in a galaxy proto-cluster at $z=2.16$. These yield the most detailed polarimetric maps yet made of a high redshift radio galaxy. The intrinsic polarization angles and Faraday Rotation Measures (RMs) reveal coherent magnetic fields spanning the $\sim60$ kpc length of the jets, while $\sim50$% fracti… ▽ More We present deep broadband radio polarization observations of the Spiderweb radio galaxy (J1140-2629) in a galaxy proto-cluster at $z=2.16$. These yield the most detailed polarimetric maps yet made of a high redshift radio galaxy. The intrinsic polarization angles and Faraday Rotation Measures (RMs) reveal coherent magnetic fields spanning the $\sim60$ kpc length of the jets, while $\sim50$% fractional polarizations indicate these fields are well-ordered. Source-frame absolute RM values of $\sim1,000$ rad/m/m are typical, and values up to $\sim11,100$ rad/m/m are observed. The Faraday-rotating gas cannot be well-mixed with the synchrotron-emitting gas, or stronger-than-observed depolarization would occur. Nevertheless, an observed spatial coincidence between a localized absolute RM enhancement of $\sim1,100$ rad/m/m, a bright knot of Ly$α$ emission, and a deviation of the radio jet provide direct evidence for vigorous jet-gas interaction. We detect a large-scale RM gradient totaling $\sim1,000$s rad/m/m across the width of the jet, suggesting a net clockwise (as viewed from the AGN) toroidal magnetic field component exists at 10s-of-kpc-scales, which we speculate may be associated with the operation of a Poynting-Robertson cosmic battery. We conclude the RMs are mainly generated in a sheath of hot gas around the radio jet, rather than the ambient foreground proto-cluster gas. The estimated magnetic field strength decreases by successive orders-of-magnitude going from the jet hotspots ($\sim90$ $μ$G) to the jet sheath ($\sim10$ $μ$G) to the ambient intracluster medium ($\sim1$ $μ$G). Synthesizing our results, we propose that the Spiderweb radio galaxy is actively magnetizing its surrounding proto-cluster environment, with possible implications for theories of the origin and evolution of cosmic magnetic fields. △ Less

Submitted 7 July, 2022; originally announced July 2022.

Comments: 12 pages, 7 figures, 2 tables. Accepted for publication in ApJ

Direct Optimal Mapping for 21cm Cosmology: A Demonstration with the Hydrogen Epoch of Reionization Array

Authors: 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, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Philip Bull, Jacob Burba , et al. (56 additional authors not shown)

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 pipe… ▽ More 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. △ Less

Submitted 26 October, 2022; v1 submitted 12 April, 2022; originally announced April 2022.

Comments: 16 pages, 10 figures, 2 tables, published on ApJ

Molecular gas in z~6 quasar host galaxies

Authors: Roberto Decarli, Antonio Pensabene, Bram Venemans, Fabian Walter, Eduardo Banados, Frank Bertoldi, Chris L. Carilli, Pierre Cox, Xiaohui Fan, Emanuele Paolo Farina, Carl Ferkinhoff, Brent A. Groves, Jianan Li, Chiara Mazzucchelli, Roberto Neri, Dominik A. Riechers, Bade Uzgil, Feige Wang, Ran Wang, Axel Weiss, Jan Martin Winters, Jinyi Yang

Abstract: 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 targete… ▽ More 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. △ Less

Submitted 7 March, 2022; originally announced March 2022.

Comments: 13 pages, 7 figures. Accepted for publication in A&A

Journal ref: A&A 662, A60 (2022)

X-ray Emission and Radio Emission from the Jets and Lobes of the Spiderweb Radio Galaxy

Authors: Christopher L. Carilli, Craig S. Anderson, Paolo Tozzi, Maurilio Pannella, Tracy Clarke, L. Pentericci, Ang Liu, Tony Mroczkowski, G. K. Miley, H. J. Rottgering, S. Borgani, Colin Norman, A. Saro, M. Nonino, L. Di Mascolo

Abstract: Deep Chandra and VLA imaging reveals a clear correlation between X-ray and radio emission on scales $\sim 100$~kpc in the Spiderweb radio galaxy at z=2.16. The X-ray emission associated with the extended radio source is likely dominated by inverse Compton up-scattering of cosmic microwave background photons by the radio emitting relativistic electrons. For regions dominated by high surface brightn… ▽ More Deep Chandra and VLA imaging reveals a clear correlation between X-ray and radio emission on scales $\sim 100$~kpc in the Spiderweb radio galaxy at z=2.16. The X-ray emission associated with the extended radio source is likely dominated by inverse Compton up-scattering of cosmic microwave background photons by the radio emitting relativistic electrons. For regions dominated by high surface brightness emission, such as hot spots and jet knots, the implied magnetic fields are $\sim 50~μ$G to $70~μ$G. The non-thermal pressure is these brighter regions is then $\sim 9\times 10^{-10}$ dyne cm$^{-2}$, or three times larger than the non-thermal pressure derived assuming minimum energy conditions, and an order of magnitude larger than the thermal pressure in the ambient cluster medium. Assuming ram pressure confinement implies an average advance speed for the radio source of $\sim 2400$ km s$^{-1}$, and a source age of $\sim 3\times 10^7$ years. Considering the lower surface brightness, diffuse radio emitting regions, we identify an evacuated cavity in the Ly$α$ emission coincident with the tail of the eastern radio lobe. Making reasonable assumptions for the radio spectrum, we find that the relativistic electrons and fields in the lobe are plausibly in pressure equilibrium with the thermal gas, and close to a minimum energy configuration. The radio morphology suggests that the Spiderweb is a high-$z$ example of the rare class of hybrid morphology radio sources (or HyMoRS), which we attribute to interaction with the asymmetric gaseous environment indicated by the Ly$α$ emission. △ Less

Submitted 7 March, 2022; originally announced March 2022.

Comments: 8 figures, 25 pages, to appear in the ApJ

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

Authors: 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 $α= -0.52\pm0.18$ ($S \propto ν^α$). The radio source is unresolved in all VLA bands wi… ▽ More 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. △ Less

Submitted 6 March, 2022; originally announced March 2022.

Comments: 9 pages, 3 figures. Accepted for publication in The Astrophysical Journal

The 700 ks Chandra Spiderweb Field I: evidence for widespread nuclear activity in the Protocluster

Authors: P. Tozzi, L. Pentericci, R. Gilli, M. Pannella, F. Fiore, G. Miley, M. Nonino, H. J. A. Rottgering, V. Strazzullo, C. S. Anderson, S. Borgani, A. Calabro', C. Carilli, H. Dannerbauer, L. Di Mascolo, C. Feruglio, R. Gobat, S. Jin, A. Liu, T. Mroczkowski, C. Norman, E. Rasia, P. Rosati, A. Saro

Abstract: (Abridged) We present an analysis of the 700 ks Chandra ACIS-S observation of the field around the Spiderweb Galaxy at z=2.156, focusing on the nuclear activity in the associated large-scale environment. We identify unresolved X-ray sources down to flux limits of 1.3X10^{-16} and 3.9X10^{-16} erg/s/cm^2 in the soft and hard band, respectively. We search for counterparts in the optical, NIR and sub… ▽ More (Abridged) We present an analysis of the 700 ks Chandra ACIS-S observation of the field around the Spiderweb Galaxy at z=2.156, focusing on the nuclear activity in the associated large-scale environment. We identify unresolved X-ray sources down to flux limits of 1.3X10^{-16} and 3.9X10^{-16} erg/s/cm^2 in the soft and hard band, respectively. We search for counterparts in the optical, NIR and submm bands to identify X-ray sources belonging to the protocluster. We detect 107 X-ray unresolved sources within 5 arcmin (corresponding to 2.5 Mpc) of J1140-2629, among which 13 have optical counterparts with spectroscopic redshift 2.11<z<2.20, and 1 source with photometric redshift consistent with this range. Our X-ray spectral analysis shows that their intrinsic spectral slope is consistent with an average <Γ>~1.84+-0.04. The best-fit intrinsic absorption for 5 protocluster X-ray members is N_H>10^{23} cm^{-2}, while other 6 have upper limits of the order of fewX10^{22} cm^{-2}. Two sources can only be fitted with very flat Γ<=1, and are therefore considered Compton-thick candidates. Their 0.5-10 keV rest frame luminosities are larger than 2X10^{43} erg/s, significantly greater than X-ray luminosities expected from star formation activity. The X-ray luminosity function of AGN in the volume associated to the Spiderweb protocluster in the range 10^{43}<L_X<10^{44.5} erg/s, is at least 10 times higher than that in the field at the same redshift and significantly flatter. The X-ray AGN fraction is measured to be (25.5+-4.5)% in the stellar mass range log(M*/M_sun)>10.5, corresponding to an enhancement of 6.0^{+9.0}_{-3.0} with respect to the COSMOS field at comparable redshifts and stellar mass range. We conclude that the galaxy population in the Spiderweb Protocluster is characterized by enhanced X-ray nuclear activity triggered by environmental effects on Mpc scales. △ Less

Submitted 4 April, 2022; v1 submitted 4 March, 2022; originally announced March 2022.

Comments: Astronomy & Astrophysics. Minor changes in this version. Version accepted in A&A

Journal ref: A&A 662, A54 (2022)

Dense Gas History of the Universe: from ASPECS to the ngVLA

Authors: C. L. Carilli, F. Walter, R. Decarli, M. Aravena, Dominik A. Riechers, J. Gonzalez-Lopez, Yali Shao, L. Boogaard, R. Bouwens, M. Neeleman

Abstract: We review the evolution of the cosmic average molecular gas density to large look-back times, using observations of rotational transitions of CO. Molecular gas is the fuel for star formation in galaxies. Deep searches for CO emission from distant galaxies have delineated the density of molecular gas back to $z \sim 5$, or within 1~Gyr of the Big Bang. The results show a rise and fall in the gas de… ▽ More We review the evolution of the cosmic average molecular gas density to large look-back times, using observations of rotational transitions of CO. Molecular gas is the fuel for star formation in galaxies. Deep searches for CO emission from distant galaxies have delineated the density of molecular gas back to $z \sim 5$, or within 1~Gyr of the Big Bang. The results show a rise and fall in the gas density that parallels, and likely drives, the rise and fall of the cosmic star formation rate density. We present the potential for the next generation Very Large Array to image the distribution and dynamics of the molecular gas in early galaxies, and to make a precise measurement of the dense gas history of the Universe. △ Less

Submitted 28 February, 2022; originally announced March 2022.

Comments: 5 pages. 4 figures. to appear in Radio Science

Microwave Background Temperature at a Redshift of 6.34 from H2O Absorption

Authors: Dominik A. Riechers, Axel Weiss, Fabian Walter, Christopher L. Carilli, Pierre Cox, Roberto Decarli, Roberto Neri

Abstract: Distortions of the observed cosmic microwave background imprinted by the Sunyaev-Zel'dovich effect toward massive galaxy clusters due to inverse Compton scattering of microwave photons by high-energy electrons provide a direct measurement of the microwave background temperature at redshifts from 0 to 1. Some additional background temperature estimates exist at redshifts from 1.8 to 3.3 based on mo… ▽ More Distortions of the observed cosmic microwave background imprinted by the Sunyaev-Zel'dovich effect toward massive galaxy clusters due to inverse Compton scattering of microwave photons by high-energy electrons provide a direct measurement of the microwave background temperature at redshifts from 0 to 1. Some additional background temperature estimates exist at redshifts from 1.8 to 3.3 based on molecular and atomic line excitation temperatures in quasar absorption line systems, but are model dependent. To date, no deviations from the expected (1+z) scaling behavior of the microwave background temperature have been seen, but the measurements have not extended deeply into the matter-dominated era of the universe at redshifts z>3.3. Here we report the detection of sub-millimeter line absorption from the water molecule against the cosmic microwave background at z=6.34 in a massive starburst galaxy, corresponding to a lookback time of 12.8 Gyr. Radiative pumping of the upper level of the ground-state ortho-H2O(110-101) line due to starburst activity in the dusty galaxy HFLS3 results in a cooling to below the redshifted microwave background temperature, after the transition is initially excited by the microwave background. The strength of this effect implies a microwave background temperature of 16.4-30.2 K (1-sigma range) at z=6.34, which is consistent with a background temperature increase with redshift as expected from the standard CDM cosmology. △ Less

Submitted 1 February, 2022; originally announced February 2022.

Comments: 16 pages, 6 figures, 1 table (including Methods and Extended Data). The version of record of this article, first published (open access) in Nature, is available online at: http://dx.doi.org/10.1038/s41586-021-04294-5

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

Authors: Fabian Walter, Marcel Neeleman, Roberto Decarli, Bram Venemans, Romain Meyer, Axel Weiss, Eduardo Banados, Sarah E. I. Bosman, Chris Carilli, Xiaohui Fan, Dominik Riechers, Hans-Walter Rix, Todd A. Thompson

Abstract: We present 0".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 ~4x1… ▽ More We present 0".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'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. △ Less

Submitted 17 January, 2022; originally announced January 2022.

Comments: Accepted for publication in The Astrophysical Journal

Simulations of primary beam effects on the cosmic bispectrum phase observed with the Hydrogen Epoch of Reionization Array

Authors: N. Charles, G. Bernardi, H. L. Bester, O. M. Smirnov, C. Carilli, P. M. Keller, N. Kern, B. Nikolic, N. Thyagarajan, E. de Lera Acedo, N. Fagnoni, M. G. Santos

Abstract: The 21~cm transition from neutral Hydrogen promises to be the best observational probe of the Epoch of Reionisation. The main difficulty in measuring the 21 cm signal is the presence of bright foregrounds that require very accurate interferometric calibration. Closure quantities may circumvent the calibration requirements but may be, however, affected by direction dependent effects, particularly a… ▽ More The 21~cm transition from neutral Hydrogen promises to be the best observational probe of the Epoch of Reionisation. The main difficulty in measuring the 21 cm signal is the presence of bright foregrounds that require very accurate interferometric calibration. Closure quantities may circumvent the calibration requirements but may be, however, affected by direction dependent effects, particularly antenna primary beam responses. This work investigates the impact of antenna primary beams affected by mutual coupling on the closure phase and its power spectrum. Our simulations show that primary beams affected by mutual coupling lead to a leakage of foreground power into the EoR window, which can be up to $\sim4$ orders magnitude higher than the case where no mutual coupling is considered. This leakage is, however, essentially confined at $k < 0.3$~$h$~Mpc$^{-1}$ for triads that include 29~m baselines. The leakage magnitude is more pronounced when bright foregrounds appear in the antenna sidelobes, as expected. Finally, we find that triads that include mutual coupling beams different from each other have power spectra similar to triads that include the same type of mutual coupling beam, indicating that beam-to-beam variation within triads (or visibility pairs) is not the major source of foreground leakage in the EoR window. △ Less

Submitted 19 May, 2022; v1 submitted 23 December, 2021; originally announced December 2021.

The radio spectral turnover of radio-loud quasars at $z>5$

Authors: Yali Shao, Jeff Wagg, Ran Wang, Emmanuel Momjian, Chris L. Carilli, Fabian Walter, Dominik A. Riechers, Huib T. Intema, Axel Weiss, Andreas Brunthaler, Karl M. Menten

Abstract: We present Karl G. Jansky Very Large Array (VLA) S- (2--4 GHz), C- (4--8 GHz), and X-band (8--12 GHz) continuum observations toward seven radio-loud quasars at $z>5$. This sample has previously been found to exhibit spectral peaks at observed-frame frequencies above $\sim$1 GHz. We also present upgraded Giant Metrewave Radio Telescope (uGMRT) band-2 (200 MHz), band-3 (400 MHz), and band-4 (650 MHz… ▽ More We present Karl G. Jansky Very Large Array (VLA) S- (2--4 GHz), C- (4--8 GHz), and X-band (8--12 GHz) continuum observations toward seven radio-loud quasars at $z>5$. This sample has previously been found to exhibit spectral peaks at observed-frame frequencies above $\sim$1 GHz. We also present upgraded Giant Metrewave Radio Telescope (uGMRT) band-2 (200 MHz), band-3 (400 MHz), and band-4 (650 MHz) radio continuum observations toward eight radio-loud quasars at $z>5$, selected from our previous GMRT survey, in order to sample their low-frequency synchrotron emission. Combined with archival radio continuum observations, all ten targets show evidence for spectral turnover. The turnover frequencies are $\sim$1--50 GHz in the rest frame, making these targets gigahertz-peaked-spectrum (GPS) or high-frequency-peaker (HFP) candidates. For the nine well-constrained targets with observations on both sides of the spectral turnover, we fit the entire radio spectrum with absorption models associated with synchrotron self-absorption and free-free absorption (FFA). Our results show that FFA in an external inhomogeneous medium can accurately describe the observed spectra for all nine targets, which may indicate an FFA origin for the radio spectral turnover in our sample. As for the complex spectrum of J114657.79+403708.6 at $z=5.00$ with two spectral peaks, it may be caused by multiple components (i.e., core-jet) and FFA by the high-density medium in the nuclear region. However, we cannot rule out the spectral turnover origin of variability. Based on our radio spectral modeling, we calculate the radio loudness $R_{2500\rm\, Å}$ for our sample, which ranges from 12$^{+1}_{-1}$ to 674$^{+61}_{-51}$. △ Less

Submitted 6 December, 2021; originally announced December 2021.

Comments: Accepted for publication in A&A

Journal ref: A&A 659, A159 (2022)

Hybrid mapping of the Black Hole Shadow in M87

Authors: Christopher L. Carilli, Nithyanandan Thyagarajan

Abstract: We present a reanalysis of the EHT 228 GHz observations of M87. We apply traditional hybrid mapping techniques to the publicly available `network-calibrated' data. We explore the impact on the final image of different starting models, including: a point source, a disk, an annulus, a Gaussian, and an asymmetric double Gaussian. The images converge to an extended source with a size $\sim 44~μ$as. St… ▽ More We present a reanalysis of the EHT 228 GHz observations of M87. We apply traditional hybrid mapping techniques to the publicly available `network-calibrated' data. We explore the impact on the final image of different starting models, including: a point source, a disk, an annulus, a Gaussian, and an asymmetric double Gaussian. The images converge to an extended source with a size $\sim 44~μ$as. Starting with the annulus and disk models leads to images with the lowest noise, smallest off-source artifacts, and better closure residuals. The source appears as a ring, or edge-brightened disk, with higher surface brightness in the southern half, consistent with previous results. Starting with the other models leads to a surface brightness distribution with a similar size, and an internal depression, but not as clearly ring-like. A consideration of visibility amplitudes vs. UV-distance argues for a roughly circularly symmetric structure of $\sim 50~μ$as scale, with a sharp-edge, based on a prominent minimum in the UV-distribution, and the amplitude of the secondary peak in the UV-plot is more consistent with an annular model than a flat disk model. With further processing, we find a possible modest extension from the ring toward the southwest, in a direction consistent with the southern limb of the jet seen on 3mm VLBI images on a factor of few larger scales. However, this extension appears along the direction of one of the principle sidelobes of the synthesized beam, and hence requires testing with better UV-coverage. △ Less

Submitted 30 November, 2021; v1 submitted 22 November, 2021; originally announced November 2021.

Comments: 16 pages, 10 figures, accepted (in press) in ApJ. Added reference. Revised title and added link to digital figure

Automated Detection of Antenna Malfunctions in Large-N Interferometers: A Case Study with the Hydrogen Epoch of Reionization Array

Authors: Dara Storer, Joshua S. Dillon, Daniel C. Jacobs, Miguel F. Morales, Bryna J. Hazelton, Aaron Ewall-Wice, Zara Abdurashidova, James E. Aguirre, Paul Alexander, Zaki S. Ali, 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, Carina Cheng, David R. DeBoer, Eloy de Lera Acedo, Matt Dexter, Scott Dynes , et al. (53 additional authors not shown)

Abstract: We present a framework for identifying and flagging malfunctioning antennas in large radio interferometers. We outline two distinct categories of metrics designed to detect outliers along known failure modes of large arrays: cross-correlation metrics, based on all antenna pairs, and auto-correlation metrics, based solely on individual antennas. We define and motivate the statistical framework for… ▽ More We present a framework for identifying and flagging malfunctioning antennas in large radio interferometers. We outline two distinct categories of metrics designed to detect outliers along known failure modes of large arrays: cross-correlation metrics, based on all antenna pairs, and auto-correlation metrics, based solely on individual antennas. We define and motivate the statistical framework for all metrics used, and present tailored visualizations that aid us in clearly identifying new and existing systematics. We implement these techniques using data from 105 antennas in the Hydrogen Epoch of Reionization Array (HERA) as a case study. Finally, we provide a detailed algorithm for implementing these metrics as flagging tools on real data sets. △ Less

Submitted 4 May, 2022; v1 submitted 26 September, 2021; originally announced September 2021.

Comments: 31 pages, 17 figures

Journal ref: Radio Science, vol. 57, no. 1, 2022

HERA Phase I Limits on the Cosmic 21-cm Signal: Constraints on Astrophysics and Cosmology During the Epoch of Reionization

Authors: The HERA Collaboration, Zara Abdurashidova, James E. Aguirre, Paul Alexander, Zaki Ali, Yanga Balfour, Rennan Barkana, Adam Beardsley, Gianni Bernardi, Tashalee Billings, Judd Bowman, Richard Bradley, Phillip Bull, Jacob Burba, Steven Carey, Christopher Carilli, Carina Cheng, David DeBoer, Matthew Dexter, Eloy de Lera Acedo, Joshua Dillon, John Ely, Aaron Ewall-Wice, Nicolas Fagnoni, Anastasia Fialkov , et al. (59 additional authors not shown)

Abstract: Recently, the Hydrogen Epoch of Reionization Array (HERA) collaboration has produced the experiment's first upper limits on the power spectrum of 21-cm fluctuations at z~8 and 10. Here, we use several independent theoretical models to infer constraints on the intergalactic medium (IGM) and galaxies during the epoch of reionization (EoR) from these limits. We find that the IGM must have been heated… ▽ More Recently, the Hydrogen Epoch of Reionization Array (HERA) collaboration has produced the experiment's first upper limits on the power spectrum of 21-cm fluctuations at z~8 and 10. Here, we use several independent theoretical models to infer constraints on the intergalactic medium (IGM) and galaxies during the epoch of reionization (EoR) from these limits. We find that the IGM must have been heated above the adiabatic cooling threshold by z~8, independent of uncertainties about the IGM ionization state and the nature of the radio background. Combining HERA limits with galaxy and EoR observations constrains the spin temperature of the z~8 neutral IGM to 27 K < T_S < 630 K (2.3 K < T_S < 640 K) at 68% (95%) confidence. They therefore also place a lower bound on X-ray heating, a previously unconstrained aspects of early galaxies. For example, if the CMB dominates the z~8 radio background, the new HERA limits imply that the first galaxies produced X-rays more efficiently than local ones (with soft band X-ray luminosities per star formation rate constrained to L_X/SFR = { 10^40.2, 10^41.9 } erg/s/(M_sun/yr) at 68% confidence), consistent with expectations of X-ray binaries in low-metallicity environments. The z~10 limits require even earlier heating if dark-matter interactions (e.g., through millicharges) cool down the hydrogen gas. Using a model in which an extra radio background is produced by galaxies, we rule out (at 95% confidence) the combination of high radio and low X-ray luminosities of L_{r,ν}/SFR > 3.9 x 10^24 W/Hz/(M_sun/yr) and L_X/SFR<10^40 erg/s/(M_sun/yr). The new HERA upper limits neither support nor disfavor a cosmological interpretation of the recent EDGES detection. The analysis framework described here provides a foundation for the interpretation of future HERA results. △ Less

Submitted 20 December, 2022; v1 submitted 16 August, 2021; originally announced August 2021.

Comments: 40 pages, 19 figures, accepted to ApJ

The Impact of Powerful Jets on the Far-infrared Emission of an Extreme Radio Quasar at z~6

Authors: Sofía Rojas-Ruiz, Eduardo Bañados, Marcel Neeleman, Thomas Connor, Anna-Christina Eilers, Bram P. Venemans, Yana Khusanova, Christopher L. Carilli, Chiara Mazzucchelli, Roberto Decarli, Emmanuel Momjian, Mladen Novak

Abstract: The interactions between radio jets and the interstellar medium play a defining role for the co-evolution of central supermassive black holes and their host galaxies, but observational constraints on these feedback processes are still very limited at redshifts $z > 2$. We investigate the radio-loud quasar PSO J352.4034-15.3373 at $z \sim 6$ at the edge of the Epoch of Reionization. This quasar is… ▽ More The interactions between radio jets and the interstellar medium play a defining role for the co-evolution of central supermassive black holes and their host galaxies, but observational constraints on these feedback processes are still very limited at redshifts $z > 2$. We investigate the radio-loud quasar PSO J352.4034-15.3373 at $z \sim 6$ at the edge of the Epoch of Reionization. This quasar is among the most powerful radio emitters and the first one with direct evidence of extended radio jets ($\sim$1.6 kpc) at these high redshifts. We analyze NOEMA and ALMA millimeter data targeting the CO (6-5) and [CII] far-infrared emission lines, respectively, and the underlying continuum. The broad $440\pm 80$ km s$^{-1}$ and marginally resolved [CII] emission line yields a systemic redshift of $z\!=\!5.832 \pm 0.001$. Additionally, we report a strong 215 MHz radio continuum detection, $88\pm 7$ mJy, using the GMRT. This measurement significantly improves the constraints at the low-frequency end of the spectral energy distribution of this quasar. In contrast to what is typically observed in high-redshift radio-quiet quasars, we show that cold dust emission alone cannot reproduce the millimeter continuum measurements. This is evidence that the strong synchrotron emission from the quasar contributes substantially to the emission even at millimeter (far-infrared in the rest-frame) wavelengths. This quasar is an ideal system to probe the effects of radio jets during the formation of a massive galaxy within the first Gyr of the Universe. △ Less

Submitted 9 August, 2021; originally announced August 2021.

Comments: 14 pages, 5 figures

First Results from HERA Phase I: Upper Limits on the Epoch of Reionization 21 cm Power Spectrum

Authors: The HERA Collaboration, Zara Abdurashidova, James E. Aguirre, Paul Alexander, Zaki S. Ali, Yanga Balfour, Adam P. Beardsley, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Philip Bull, Jacob Burba, Steve Carey, Chris L. Carilli, Carina Cheng, David R. DeBoer, Matt Dexter, Eloy de Lera Acedo, Taylor Dibblee-Barkman, Joshua S. Dillon, John Ely, Aaron Ewall-Wice, Nicolas Fagnoni, Randall Fritz , et al. (52 additional authors not shown)

Abstract: We report upper-limits on the Epoch of Reionization (EoR) 21 cm power spectrum at redshifts 7.9 and 10.4 with 18 nights of data ($\sim36$ hours of integration) from Phase I of the Hydrogen Epoch of Reionization Array (HERA). The Phase I data show evidence for systematics that can be largely suppressed with systematic models down to a dynamic range of $\sim10^9$ with respect to the peak foreground… ▽ More We report upper-limits on the Epoch of Reionization (EoR) 21 cm power spectrum at redshifts 7.9 and 10.4 with 18 nights of data ($\sim36$ hours of integration) from Phase I of the Hydrogen Epoch of Reionization Array (HERA). The Phase I data show evidence for systematics that can be largely suppressed with systematic models down to a dynamic range of $\sim10^9$ with respect to the peak foreground power. This yields a 95% confidence upper limit on the 21 cm power spectrum of $Δ^2_{21} \le (30.76)^2\ {\rm mK}^2$ at $k=0.192\ h\ {\rm Mpc}^{-1}$ at $z=7.9$, and also $Δ^2_{21} \le (95.74)^2\ {\rm mK}^2$ at $k=0.256\ h\ {\rm Mpc}^{-1}$ at $z=10.4$. At $z=7.9$, these limits are the most sensitive to-date by over an order of magnitude. While we find evidence for residual systematics at low line-of-sight Fourier $k_\parallel$ modes, at high $k_\parallel$ modes we find our data to be largely consistent with thermal noise, an indicator that the system could benefit from deeper integrations. The observed systematics could be due to radio frequency interference, cable sub-reflections, or residual instrumental cross-coupling, and warrant further study. This analysis emphasizes algorithms that have minimal inherent signal loss, although we do perform a careful accounting in a companion paper of the small forms of loss or bias associated with the pipeline. Overall, these results are a promising first step in the development of a tuned, instrument-specific analysis pipeline for HERA, particularly as Phase II construction is completed en route to reaching the full sensitivity of the experiment. △ Less

Submitted 4 August, 2021; originally announced August 2021.

Comments: Accepted to ApJ. https://reionization.org/science/public-data-release-1/

Measuring the average molecular gas content of star-forming galaxies at $z=3-4$

Authors: Leindert A. Boogaard, Rychard J. Bouwens, Dominik Riechers, Paul van der Werf, Roland Bacon, Jorryt Matthee, Mauro Stefanon, Anna Feltre, Michael Maseda, Hanae Inami, Manuel Aravena, Jarle Brinchmann, Chris Carilli, Thierry Contini, Roberto Decarli, Jorge González-López, Themiya Nanayakkara, Fabian Walter

Abstract: We study the molecular gas content of 24 star-forming galaxies at $z=3-4$, with a median stellar mass of $10^{9.1}$ M$_{\odot}$, from the MUSE Hubble Ultra Deep Field (HUDF) Survey. Selected by their Lyman-alpha-emission and H-band magnitude, the galaxies show an average EW $\approx 20$ angstrom, below the typical selection threshold for Lyman Alpha Emitters (EW $> 25$ angstrom), and a rest-frame… ▽ More We study the molecular gas content of 24 star-forming galaxies at $z=3-4$, with a median stellar mass of $10^{9.1}$ M$_{\odot}$, from the MUSE Hubble Ultra Deep Field (HUDF) Survey. Selected by their Lyman-alpha-emission and H-band magnitude, the galaxies show an average EW $\approx 20$ angstrom, below the typical selection threshold for Lyman Alpha Emitters (EW $> 25$ angstrom), and a rest-frame UV spectrum similar to Lyman Break Galaxies. We use rest-frame optical spectroscopy from KMOS and MOSFIRE, and the UV features observed with MUSE, to determine the systemic redshifts, which are offset from Lyman alpha by 346 km s$^{-1}$, with a 100 to 600 km s$^{-1}$ range. Stacking CO(4-3) and [CI](1-0) (and higher-$J$ CO lines) from the ALMA Spectroscopic Survey of the HUDF (ASPECS), we determine $3σ$ upper limits on the line luminosities of $4.0\times10^{8}$ K km s$^{-1}$pc$^{2}$ and $5.6\times10^{8}$ K km s$^{-1}$pc$^{2}$, respectively (for a 300 km s$^{-1}$ linewidth). Stacking the 1.2 mm and 3 mm dust continuum flux densities, we find a $3σ$ upper limits of 9 $μ$Jy and $1.2$ $μ$Jy, respectively. The inferred gas fractions, under the assumption of a 'Galactic' CO-to-H$_{2}$ conversion factor and gas-to-dust ratio, are in tension with previously determined scaling relations. This implies a substantially higher $α_{\rm CO} \ge 10$ and $δ_{\rm GDR} \ge 1200$, consistent with the sub-solar metallicity estimated for these galaxies ($12 + \log(O/H) \approx 7.8 \pm 0.2$). The low metallicity of $z \ge 3$ star-forming galaxies may thus make it very challenging to unveil their cold gas through CO or dust emission, warranting further exploration of alternative tracers, such as [CII]. △ Less

Submitted 26 May, 2021; originally announced May 2021.

Comments: 26 pages, 12 figures, 3 tables. Accepted for publication in ApJ

Effects of model incompleteness on the drift-scan calibration of radio telescopes

Authors: Bharat K. Gehlot, Daniel C. Jacobs, Judd D. Bowman, Nivedita Mahesh, Steven G. Murray, Matthew Kolopanis, Adam P. Beardsley, Zara Abdurashidova, James E. Aguirre, Paul Alexander, Zaki S. Ali, Yanga Balfour, Gianni Bernardi, Tashalee S. Billings, Richard F. Bradley, Phil Bull, Jacob Burba, Steve Carey, Chris L. Carilli, Carina Cheng, David R. DeBoer, Matt Dexter, Eloy de Lera Acedo, Joshua S. Dillon, John Ely , et al. (54 additional authors not shown)

Abstract: Precision calibration poses challenges to experiments probing the redshifted 21-cm signal of neutral hydrogen from the Cosmic Dawn and Epoch of Reionization (z~30-6). In both interferometric and global signal experiments, systematic calibration is the leading source of error. Though many aspects of calibration have been studied, the overlap between the two types of instruments has received less at… ▽ More Precision calibration poses challenges to experiments probing the redshifted 21-cm signal of neutral hydrogen from the Cosmic Dawn and Epoch of Reionization (z~30-6). In both interferometric and global signal experiments, systematic calibration is the leading source of error. Though many aspects of calibration have been studied, the overlap between the two types of instruments has received less attention. We investigate the sky based calibration of total power measurements with a HERA dish and an EDGES style antenna to understand the role of auto-correlations in the calibration of an interferometer and the role of sky in calibrating a total power instrument. Using simulations we study various scenarios such as time variable gain, incomplete sky calibration model, and primary beam model. We find that temporal gain drifts, sky model incompleteness, and beam inaccuracies cause biases in the receiver gain amplitude and the receiver temperature estimates. In some cases, these biases mix spectral structure between beam and sky resulting in spectrally variable gain errors. Applying the calibration method to the HERA and EDGES data, we find good agreement with calibration via the more standard methods. Although instrumental gains are consistent with beam and sky errors similar in scale to those simulated, the receiver temperatures show significant deviations from expected values. While we show that it is possible to partially mitigate biases due to model inaccuracies by incorporating a time-dependent gain model in calibration, the resulting errors on calibration products are larger and more correlated. Completely addressing these biases will require more accurate sky and primary beam models. △ Less

Submitted 15 July, 2021; v1 submitted 25 April, 2021; originally announced April 2021.

Comments: 16 pages, 13 figures, 1 table; accepted for publication in MNRAS main journal

Validation of the HERA Phase I Epoch of Reionization 21 cm Power Spectrum Software Pipeline

Authors: James E. Aguirre, Steven G. Murray, Robert Pascua, Zachary E. Martinot, Jacob Burba, Joshua S. Dillon, Daniel C. Jacobs, Nicholas S. Kern, Piyanat Kittiwisit, Matthew Kolopanis, Adam Lanman, Adrian Liu, Lily Whitler, Zara Abdurashidova, Paul Alexander, Zaki S. Ali, Yanga Balfour, Adam P. Beardsley, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Philip Bull, Steve Carey, Chris L. Carilli , et al. (51 additional authors not shown)

Abstract: We describe the validation of the HERA Phase I software pipeline by a series of modular tests, building up to an end-to-end simulation. The philosophy of this approach is to validate the software and algorithms used in the Phase I upper limit analysis on wholly synthetic data satisfying the assumptions of that analysis, not addressing whether the actual data meet these assumptions. We discuss the… ▽ More We describe the validation of the HERA Phase I software pipeline by a series of modular tests, building up to an end-to-end simulation. The philosophy of this approach is to validate the software and algorithms used in the Phase I upper limit analysis on wholly synthetic data satisfying the assumptions of that analysis, not addressing whether the actual data meet these assumptions. We discuss the organization of this validation approach, the specific modular tests performed, and the construction of the end-to-end simulations. We explicitly discuss the limitations in scope of the current simulation effort. With mock visibility data generated from a known analytic power spectrum and a wide range of realistic instrumental effects and foregrounds, we demonstrate that the current pipeline produces power spectrum estimates that are consistent with known analytic inputs to within thermal noise levels (at the 2 sigma level) for k > 0.2 h/Mpc for both bands and fields considered. Our input spectrum is intentionally amplified to enable a strong `detection' at k ~0.2 h/Mpc -- at the level of ~25 sigma -- with foregrounds dominating on larger scales, and thermal noise dominating at smaller scales. Our pipeline is able to detect this amplified input signal after suppressing foregrounds with a dynamic range (foreground to noise ratio) of > 10^7. Our validation test suite uncovered several sources of scale-independent signal loss throughout the pipeline, whose amplitude is well-characterized and accounted for in the final estimates. We conclude with a discussion of the steps required for the next round of data analysis. △ Less

Submitted 19 April, 2021; originally announced April 2021.

Comments: 32 pages, 20 figures. Submitted to the Astrophysical Journal

A Real Time Processing System for Big Data in Astronomy: Applications to HERA

Authors: Paul La Plante, Peter K. G. Williams, Matthew Kolopanis, Joshua S. Dillon, Adam P. Beardsley, Nicholas S. Kern, Michael Wilensky, Zaki S. Ali, Zara Abdurashidova, James E. Aguirre, Paul Alexander, Yanga Balfour, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Phil Bull, Jacob Burba, Steve Carey, Chris L. Carilli, Carina Cheng, David R. DeBoer, Matt Dexter, Eloy de Lera Acedo, John Ely , et al. (50 additional authors not shown)

Abstract: As current- and next-generation astronomical instruments come online, they will generate an unprecedented deluge of data. Analyzing these data in real time presents unique conceptual and computational challenges, and their long-term storage and archiving is scientifically essential for generating reliable, reproducible results. We present here the real-time processing (RTP) system for the Hydrogen… ▽ More As current- and next-generation astronomical instruments come online, they will generate an unprecedented deluge of data. Analyzing these data in real time presents unique conceptual and computational challenges, and their long-term storage and archiving is scientifically essential for generating reliable, reproducible results. We present here the real-time processing (RTP) system for the Hydrogen Epoch of Reionization Array (HERA), a radio interferometer endeavoring to provide the first detection of the highly redshifted 21 cm signal from Cosmic Dawn and the Epoch of Reionization by an interferometer. The RTP system consists of analysis routines run on raw data shortly after they are acquired, such as calibration and detection of radio-frequency interference (RFI) events. RTP works closely with the Librarian, the HERA data storage and transfer manager which automatically ingests data and transfers copies to other clusters for post-processing analysis. Both the RTP system and the Librarian are public and open source software, which allows for them to be modified for use in other scientific collaborations. When fully constructed, HERA is projected to generate over 50 terabytes (TB) of data each night, and the RTP system enables the successful scientific analysis of these data. △ Less

Submitted 30 September, 2021; v1 submitted 8 April, 2021; originally announced April 2021.

Comments: 15 pages, 3 figures, published in Astronomy and Computing

Journal ref: Astronomy and Computing 36 (2021) 100489

Methods of Error Estimation for Delay Power Spectra in $21\,\textrm{cm}$ Cosmology

Authors: Jianrong Tan, Adrian Liu, Nicholas S. Kern, Zara Abdurashidova, James E. Aguirre, Paul Alexander, Zaki S. Ali, Yanga Balfour, Adam P. Beardsley, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Philip Bull, Jacob Burba, Steven Carey, Christopher L. Carilli, Carina Cheng, David R. DeBoer, Matt Dexter, Eloy de Lera Acedo, Joshua S. Dillon, John Ely, Aaron Ewall-Wice, Nicolas Fagnoni , et al. (49 additional authors not shown)

Abstract: Precise measurements of the 21 cm power spectrum are crucial for understanding the physical processes of hydrogen reionization. Currently, this probe is being pursued by low-frequency radio interferometer arrays. As these experiments come closer to making a first detection of the signal, error estimation will play an increasingly important role in setting robust measurements. Using the delay power… ▽ More Precise measurements of the 21 cm power spectrum are crucial for understanding the physical processes of hydrogen reionization. Currently, this probe is being pursued by low-frequency radio interferometer arrays. As these experiments come closer to making a first detection of the signal, error estimation will play an increasingly important role in setting robust measurements. Using the delay power spectrum approach, we have produced a critical examination of different ways that one can estimate error bars on the power spectrum. We do this through a synthesis of analytic work, simulations of toy models, and tests on small amounts of real data. We find that, although computed independently, the different error bar methodologies are in good agreement with each other in the noise-dominated regime of the power spectrum. For our preferred methodology, the predicted probability distribution function is consistent with the empirical noise power distributions from both simulated and real data. This diagnosis is mainly in support of the forthcoming HERA upper limit, and also is expected to be more generally applicable. △ Less

Submitted 25 May, 2021; v1 submitted 17 March, 2021; originally announced March 2021.

Comments: 35 Pages, 9 Figures, 4 Tables. Replaced with accepted ApJ version; some clarifying text added in response to referee comments with no changes to results

Enhanced X-ray Emission from the Most Radio-Powerful Quasar in the Universe's First Billion Years

Authors: Thomas Connor, Eduardo Bañados, Daniel Stern, Chris Carilli, Andrew Fabian, Emmanuel Momjian, Sofía Rojas-Ruiz, Roberto Decarli, Emanuele Paolo Farina, Chiara Mazzucchelli, Hannah P. Earnshaw

Abstract: We present deep (265 ks) Chandra X-ray observations of PSO J352.4034$-$15.3373, a quasar at z=5.831 that, with a radio-to-optical flux ratio of R>1000, is one of the radio-loudest quasars in the early universe and is the only quasar with observed extended radio jets of kpc-scale at $z \gtrsim 6$. Modeling the X-ray spectrum of the quasar with a power law, we find a best fit of… ▽ More We present deep (265 ks) Chandra X-ray observations of PSO J352.4034$-$15.3373, a quasar at z=5.831 that, with a radio-to-optical flux ratio of R>1000, is one of the radio-loudest quasars in the early universe and is the only quasar with observed extended radio jets of kpc-scale at $z \gtrsim 6$. Modeling the X-ray spectrum of the quasar with a power law, we find a best fit of $Γ= 1.99^{+0.29}_{-0.28}$, leading to an X-ray luminosity of $L_{2-10} = 1.26^{+0.45}_{-0.33} \times 10^{45}\ {\rm erg}\ {\rm s}^{-1}$ and an X-ray to UV brightness ratio of $α_{\rm OX} = -1.36 \pm 0.11$. We identify a diffuse structure 50 kpc (${\sim}8^{\prime\prime}$) to the NW of the quasar along the jet axis that corresponds to a $3σ$ enhancement in the angular density of emission and can be ruled out as a background fluctuation with a probability of P=0.9985. While with few detected photons the spectral fit of the structure is uncertain, we find that it has a luminosity of $L_{2-10}\sim10^{44}\ {\rm erg}\ {\rm s}^{-1}$. These observations therefore potentially represent the most distant quasar jet yet seen in X-rays. We find no evidence for excess X-ray emission where the previously-reported radio jets are seen (which have an overall linear extent of $0.^{\prime\prime}28$), and a bright X-ray point source located along the jet axis to the SE is revealed by optical and NIR imaging to not be associated with the quasar. △ Less

Submitted 5 March, 2021; originally announced March 2021.

Comments: 16 pages, 7 Figures. Accepted for publication the Astrophysical Journal