Showing 1–50 of 94 results for author: Thyagarajan, N

21cm Epoch of Reionisation Power Spectrum with Closure Phase using the Murchison Widefield Array

Authors: Himanshu Tiwari, Nithyanandan Thyagarajan, Cathryn M. Trott, Benjamin McKinley

Abstract: The radio interferometric closure phases can be a valuable tool for studying cosmological {H\scriptsize{I}}~from the early Universe. Closure phases have the advantage of being immune to element-based gains and associated calibration errors. Thus, calibration and errors therein, which are often sources of systematics limiting standard visibility-based approaches, can be avoided altogether in closur… ▽ More The radio interferometric closure phases can be a valuable tool for studying cosmological {H\scriptsize{I}}~from the early Universe. Closure phases have the advantage of being immune to element-based gains and associated calibration errors. Thus, calibration and errors therein, which are often sources of systematics limiting standard visibility-based approaches, can be avoided altogether in closure phase analysis. In this work, we present the first results of the closure phase power spectrum of {H\scriptsize{I}}~21-cm fluctuations using the Murchison Widefield Array (MWA), with $\sim 12$ hours of MWA-phase II observations centered around redshift, $z\approx 6.79$, during the Epoch of Reionisation. On analysing three redundant classes of baselines -- 14~m, 24~m, and 28~m equilateral triads, our estimates of the $2σ$ ($95\%$ confidence interval) 21-cm power spectra are $\lesssim (184)^2 pseudo \rm ~mK^2$ at ${k}_{||} = 0.36 $ $pseudo~h {\rm Mpc^{-1}}$ in the EoR1 field for the 14~m baseline triads, and $\lesssim (188)^2 pseudo \rm ~mK^2$ at $k_{||} = 0.18 $ $pseudo~h {\rm Mpc^{-1}}$ in the EoR0 field for the 24~m baseline triads. The ``$pseudo$'' units denote that the length scale and brightness temperature should be interpreted as close approximations. Our best estimates are still 3-4 orders high compared to the fiducial 21-cm power spectrum; however, our approach provides promising estimates of the power spectra even with a small amount of data. These data-limited estimates can be further improved if more datasets are included into the analysis. The evidence for excess noise has a possible origin in baseline-dependent systematics in the MWA data that will require careful baseline-based strategies to mitigate, even in standard visibility-based approaches. △ Less

Submitted 4 September, 2024; originally announced September 2024.

Comments: 24 pages, 21 figures, 4 tables (accepted for publication in PASA)

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

Two-dimensional Synchrotron Beam Characterisation from a Single Interferogram

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

Abstract: Double-aperture Young interferometry is widely used in accelerators to provide a one-dimensional beam measurement. We improve this technique by combining and further developing techniques of non-redundant aperture masking and self-calibration from astronomy. Using visible synchrotron radiation, tests at the ALBA synchrotron show that this method provides an accurate two-dimensional beam transverse… ▽ More Double-aperture Young interferometry is widely used in accelerators to provide a one-dimensional beam measurement. We improve this technique by combining and further developing techniques of non-redundant aperture masking and self-calibration from astronomy. Using visible synchrotron radiation, tests at the ALBA synchrotron show that this method provides an accurate two-dimensional beam transverse characterisation, even from a single 1ms interferogram. The technique is resistant to phase fluctuations that might be introduced by vibration of optical components, or in the laboratory atmosphere. △ Less

Submitted 23 May, 2024; v1 submitted 20 May, 2024; originally announced May 2024.

Comments: 5 pages, 5 figures, submitted to PRL

Radio-loud fraction of z>6 quasars

Authors: Pascal M. Keller, Nithyanandan Thyagarajan, Ajay Kumar, Nissim Kanekar, Gianni Bernardi

Abstract: Quasars at redshifts $z>6$ are an excellent probe of the formation and evolution of supermassive black holes in the early Universe. The population of radio-luminous quasars is of particular interest, as such quasars could potentially be used to study the neutral intergalactic medium during cosmic reionisation via H$\,$I 21$\,$cm absorption studies. However, the lack of deep radio observations of… ▽ More Quasars at redshifts $z>6$ are an excellent probe of the formation and evolution of supermassive black holes in the early Universe. The population of radio-luminous quasars is of particular interest, as such quasars could potentially be used to study the neutral intergalactic medium during cosmic reionisation via H$\,$I 21$\,$cm absorption studies. However, the lack of deep radio observations of $z>6$ quasars leaves the population poorly constrained, and suitable candidates for an H$\,$I 21$\,$cm absorption study have yet to be found. In this work, we present Jansky Very Large Array (VLA) 1$-$2 GHz radio continuum observations of 138 quasars at redshifts $6.0 \leq z<7.6$. We detect the radio continuum emission of the $z=6.1$ quasar J1034-1425, with a 1.6 GHz flux density of $170\pm 36\,μ$Jy. This quasar is radio-quiet with radio-loudness, $R \equiv f_{5\text{~GHz}}/f_{ν,\text{4400 A}} = 2.4\pm0.5$. In addition, we detect 7 other quasars at z>6, which have previously been characterised in the literature at these frequencies. Using the full sample, we estimate the radio-loud fraction to be $3.8^{+6.2}_{-2.4}\%$, where the uncertainties are 95% confidence intervals. This is lower than recent estimates of the radio-loud fraction in the literature, but is still marginally consistent with no redshift evolution of the radio-loud fraction. We explore the undetected quasar population by stacking their continuum images at their optical positions and obtain a median stacked flux density of $13.8\pm 3.9~μ$Jy and luminosity of $\log{L_{5~\mathrm{GHz}}/(\mathrm{W~Hz}^{-1})}=24.2\pm0.1$. △ Less

Submitted 13 February, 2024; originally announced February 2024.

Comments: 12 pages, 8 figures, 3 tables, published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, stae418 (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

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

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

Interferometric Image Reconstruction using Closure Invariants and Machine Learning

Authors: Nithyanandan Thyagarajan, Lucas Hoefs, O. Ivy Wong

Abstract: Interferometric closure invariants encode calibration-independent details of an object's morphology. Excepting simple cases, a direct backward transformation from closure invariants to morphologies is not well established. We demonstrate using simple Machine Learning models that closure invariants can aid in morphological classification and parameter estimation. We consider six phenomenologically… ▽ More Interferometric closure invariants encode calibration-independent details of an object's morphology. Excepting simple cases, a direct backward transformation from closure invariants to morphologies is not well established. We demonstrate using simple Machine Learning models that closure invariants can aid in morphological classification and parameter estimation. We consider six phenomenologically parametrised morphologies: point-like, uniform circular disc, crescent, dual disc, crescent with elliptical accretion disc, and crescent with double jet lobes. Using logistic regression (LR), multi-layer perceptron (MLP), and random forest models on closure invariants obtained from a sparsely covered aperture, we find that all methods except LR can classify morphologies with $\gtrsim$80% accuracy, which improves with greater aperture coverage. Separately from the classification problem, given an independently confirmed class, we estimate parameters of uniform circular disc, crescent, and dual disc morphologies using simple MLP models, and parametrically reconstruct images. The estimated parameters and images correspond well with inputs, but the accuracy worsens when degeneracies between parameters are present. This independent approach to interferometric imaging under challenging observing conditions such as that faced by the Event Horizon Telescope and Very Long Baseline Interferometry in general can complement other methods in robustly constraining an object's morphology. △ Less

Submitted 25 August, 2024; v1 submitted 10 November, 2023; originally announced November 2023.

Comments: Published in RAS Techniques and Instruments (RASTI) special edition "Next-Generation Interferometric Image Reconstruction". 17 pages including appendix, 38 figures grouped into 13 captioned figures

Simulating the Detection of the Global 21 cm Signal with MIST for Different Models of the Soil and Beam Directivity

Authors: Raul A. Monsalve, Christian H. Bye, Jonathan L. Sievers, Vadym Bidula, Ricardo Bustos, H. Cynthia Chiang, Xinze Guo, Ian Hendricksen, Francis McGee, F. Patricio Mena, Garima Prabhakar, Oscar Restrepo, Nithyanandan Thyagarajan

Abstract: The Mapper of the IGM Spin Temperature (MIST) is a new ground-based, single-antenna, radio experiment attempting to detect the global 21 cm signal from the Dark Ages and Cosmic Dawn. A significant challenge in this measurement is the frequency-dependence, or chromaticity, of the antenna beam directivity. MIST observes with the antenna above the soil and without a metal ground plane, and the beam d… ▽ More The Mapper of the IGM Spin Temperature (MIST) is a new ground-based, single-antenna, radio experiment attempting to detect the global 21 cm signal from the Dark Ages and Cosmic Dawn. A significant challenge in this measurement is the frequency-dependence, or chromaticity, of the antenna beam directivity. MIST observes with the antenna above the soil and without a metal ground plane, and the beam directivity is sensitive to the electrical characteristics of the soil. In this paper, we use simulated observations with MIST to study how the detection of the global 21 cm signal from Cosmic Dawn is affected by the soil and the MIST beam directivity. We simulate observations using electromagnetic models of the directivity computed for single- and two-layer models of the soil. We test the recovery of the Cosmic Dawn signal with and without beam chromaticity correction applied to the simulated data. We find that our single-layer soil models enable a straightforward recovery of the signal even without chromaticity correction. Two-layer models increase the beam chromaticity and make the recovery more challenging. However, for the model in which the bottom soil layer has a lower electrical conductivity than the top layer, the signal can be recovered even without chromaticity correction. For the other two-layer models, chromaticity correction is necessary for the recovery of the signal and the accuracy requirements for the soil parameters vary between models. These results will be used as a guideline to select observation sites that are favorable for the detection of the Cosmic Dawn signal. △ Less

Submitted 23 May, 2024; v1 submitted 10 October, 2023; originally announced October 2023.

Comments: Matches version published in The Astrophysical Journal. Journal reference and DOI provided

Journal ref: 2024, ApJ, 961, 56

Mapper of the IGM spin temperature: instrument overview

Authors: R. A. Monsalve, C. Altamirano, V. Bidula, R. Bustos, C. H. Bye, H. C. Chiang, M. Diaz, B. Fernandez, X. Guo, I. Hendricksen, E. Hornecker, F. Lucero, H. Mani, F. McGee, F. P. Mena, M. Pessoa, G. Prabhakar, O. Restrepo, J. L. Sievers, N. Thyagarajan

Abstract: The observation of the global 21 cm signal produced by neutral hydrogen gas in the intergalactic medium (IGM) during the Dark Ages, Cosmic Dawn, and Epoch of Reionization requires measurements with extremely well-calibrated wideband radiometers. We describe the design and characterization of the Mapper of the IGM Spin Temperature (MIST), which is a new ground-based, single-antenna, global 21 cm ex… ▽ More The observation of the global 21 cm signal produced by neutral hydrogen gas in the intergalactic medium (IGM) during the Dark Ages, Cosmic Dawn, and Epoch of Reionization requires measurements with extremely well-calibrated wideband radiometers. We describe the design and characterization of the Mapper of the IGM Spin Temperature (MIST), which is a new ground-based, single-antenna, global 21 cm experiment. The design of MIST was guided by the objectives of avoiding systematics from an antenna ground plane and cables around the antenna, as well as maximizing the instrument's on-sky efficiency and portability for operations at remote sites. We have built two MIST instruments, which observe in the range 25-105 MHz. For the 21 cm signal, this frequency range approximately corresponds to redshifts 55.5 > z > 12.5, encompassing the Dark Ages and Cosmic Dawn. The MIST antenna is a horizontal blade dipole of 2.42 m in length, 60 cm in width, and 52 cm in height above the ground. This antenna operates without a metal ground plane. The instruments run on 12 V batteries and have a maximum power consumption of 17 W. The batteries and electronics are contained in a single receiver box located under the antenna. We present the characterization of the instruments using electromagnetic simulations and lab measurements. We also show sample sky measurements from recent observations at remote sites in California, Nevada, and the Canadian High Arctic. These measurements indicate that the instruments perform as expected. Detailed analyses of the sky measurements are left for future work. △ Less

Submitted 23 May, 2024; v1 submitted 6 September, 2023; originally announced September 2023.

Comments: Matches version published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 530, Issue 4, June 2024, Pages 4125-4147

Comparing recent PTA results on the nanohertz stochastic gravitational wave background

Authors: The International Pulsar Timing Array Collaboration, G. Agazie, J. Antoniadis, A. Anumarlapudi, A. M. Archibald, P. Arumugam, S. Arumugam, Z. Arzoumanian, J. Askew, S. Babak, M. Bagchi, M. Bailes, A. -S. Bak Nielsen, P. T. Baker, C. G. Bassa, A. Bathula, B. Bécsy, A. Berthereau, N. D. R. Bhat, L. Blecha, M. Bonetti, E. Bortolas, A. Brazier, P. R. Brook, M. Burgay , et al. (220 additional authors not shown)

Abstract: The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTA… ▽ More The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTAs that constitute the International Pulsar Timing Array (IPTA). We show that despite making different modeling choices, there is no significant difference in the GWB parameters that are measured by the different PTAs, agreeing within $1σ$. The pulsar noise parameters are also consistent between different PTAs for the majority of the pulsars included in these analyses. We bridge the differences in modeling choices by adopting a standardized noise model for all pulsars and PTAs, finding that under this model there is a reduction in the tension in the pulsar noise parameters. As part of this reanalysis, we "extended" each PTA's data set by adding extra pulsars that were not timed by that PTA. Under these extensions, we find better constraints on the GWB amplitude and a higher signal-to-noise ratio for the Hellings and Downs correlations. These extensions serve as a prelude to the benefits offered by a full combination of data across all pulsars in the IPTA, i.e., the IPTA's Data Release 3, which will involve not just adding in additional pulsars, but also including data from all three PTAs where any given pulsar is timed by more than as single PTA. △ Less

Submitted 1 September, 2023; originally announced September 2023.

Comments: 21 pages, 9 figures, submitted to ApJ

The Parkes Pulsar Timing Array Third Data Release

Authors: Andrew Zic, Daniel J. Reardon, Agastya Kapur, George Hobbs, Rami Mandow, Małgorzata Curyło, Ryan M. Shannon, Jacob Askew, Matthew Bailes, N. D. Ramesh Bhat, Andrew Cameron, Zu-Cheng Chen, Shi Dai, Valentina Di Marco, Yi Feng, Matthew Kerr, Atharva Kulkarni, Marcus E. Lower, Rui Luo, Richard N. Manchester, Matthew T. Miles, Rowina S. Nathan, Stefan Osłowski, Axl F. Rogers, Christopher J. Russell , et al. (9 additional authors not shown)

Abstract: We present the third data release from the Parkes Pulsar Timing Array (PPTA) project. The release contains observations of 32 pulsars obtained using the 64-m Parkes "Murriyang" radio telescope. The data span is up to 18 years with a typical cadence of 3 weeks. This data release is formed by combining an updated version of our second data release with $\sim 3$ years of more recent data primarily ob… ▽ More We present the third data release from the Parkes Pulsar Timing Array (PPTA) project. The release contains observations of 32 pulsars obtained using the 64-m Parkes "Murriyang" radio telescope. The data span is up to 18 years with a typical cadence of 3 weeks. This data release is formed by combining an updated version of our second data release with $\sim 3$ years of more recent data primarily obtained using an ultra-wide-bandwidth receiver system that operates between 704 and 4032 MHz. We provide calibrated pulse profiles, flux-density dynamic spectra, pulse times of arrival, and initial pulsar timing models. We describe methods for processing such wide-bandwidth observations, and compare this data release with our previous release. △ Less

Submitted 17 October, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 15 pages, 6 figures. Accepted for publication in PASA

Measuring the global 21-cm signal with the MWA-II: improved characterisation of lunar-reflected radio frequency interference

Authors: Himanshu Tiwari, Benjamin McKinley, Cathryn M. Trott, Nithyanandan Thyagarajan

Abstract: Radio interferometers can potentially detect the sky-averaged signal from the Cosmic Dawn (CD) and the Epoch of Reionisation (EoR) by studying the Moon as a thermal block to the foreground sky. The first step is to mitigate the Earth-based RFI reflections (Earthshine) from the Moon, which significantly contaminate the FM band $\approx 88-110$ MHz, crucial to CD-EoR science. We analysed MWA phase-I… ▽ More Radio interferometers can potentially detect the sky-averaged signal from the Cosmic Dawn (CD) and the Epoch of Reionisation (EoR) by studying the Moon as a thermal block to the foreground sky. The first step is to mitigate the Earth-based RFI reflections (Earthshine) from the Moon, which significantly contaminate the FM band $\approx 88-110$ MHz, crucial to CD-EoR science. We analysed MWA phase-I data from $72-180$ MHz at $40$ kHz resolution to understand the nature of Earthshine over three observing nights. We took two approaches to correct the Earthshine component from the Moon. In the first method, we mitigated the Earthshine using the flux density of the two components from the data, while in the second method, we used simulated flux density based on an FM catalogue to mitigate the Earthshine. Using these methods, we were able to recover the expected Galactic foreground temperature of the patch of sky obscured by the Moon. We performed a joint analysis of the Galactic foregrounds and the Moon's intrinsic temperature $(T_{\rm Moon})$ while assuming that the Moon has a constant thermal temperature throughout three epochs. We found $T_{\rm Moon}$ to be at $184.4\pm{2.6}\rm ~K$ and $173.8\pm{2.5}\rm ~K$ using the first and the second methods, respectively, and the best-fit values of the Galactic spectral index $(α)$ to be within the $5\%$ uncertainty level when compared with the global sky models. Compared with our previous work, these results improved constraints on the Galactic spectral index and the Moon's intrinsic temperature. We also simulated the Earthshine at MWA between November and December 2023 to find suitable observing times less affected by the Earthshine. Such observing windows act as Earthshine avoidance and can be used to perform future global CD-EoR experiments using the Moon with the MWA. △ Less

Submitted 1 November, 2023; v1 submitted 31 May, 2023; originally announced June 2023.

Comments: 17 pages, 14 figures and 8 tables, accepted for publication in PASA

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

Optimization and Commissioning of the EPIC Commensal Radio Transient Imager for the Long Wavelength Array

Authors: Hariharan Krishnan, Adam P. Beardsley, Judd D. Bowman, Jayce Dowell, Matthew Kolopanis, Greg Taylor, Nithyanandan Thyagarajan

Abstract: Next generation aperture arrays are expected to consist of hundreds to thousands of antenna elements with substantial digital signal processing to handle large operating bandwidths of a few tens to hundreds of MHz. Conventionally, FX~correlators are used as the primary signal processing unit of the interferometer. These correlators have computational costs that scale as $\mathcal{O}(N^2)$ for larg… ▽ More Next generation aperture arrays are expected to consist of hundreds to thousands of antenna elements with substantial digital signal processing to handle large operating bandwidths of a few tens to hundreds of MHz. Conventionally, FX~correlators are used as the primary signal processing unit of the interferometer. These correlators have computational costs that scale as $\mathcal{O}(N^2)$ for large arrays. An alternative imaging approach is implemented in the E-field Parallel Imaging Correlator (EPIC) that was recently deployed on the Long Wavelength Array station at the Sevilleta National Wildlife Refuge (LWA-SV) in New Mexico. EPIC uses a novel architecture that produces electric field or intensity images of the sky at the angular resolution of the array with full or partial polarization and the full spectral resolution of the channelizer. By eliminating the intermediate cross-correlation data products, the computational costs can be significantly lowered in comparison to a conventional FX~or XF~correlator from $\mathcal{O}(N^2)$ to $\mathcal{O}(N \log N)$ for dense (but otherwise arbitrary) array layouts. EPIC can also lower the output data rates by directly yielding polarimetric image products for science analysis. We have optimized EPIC and have now commissioned it at LWA-SV as a commensal all-sky imaging back-end that can potentially detect and localize sources of impulsive radio emission on millisecond timescales. In this article, we review the architecture of EPIC, describe code optimizations that improve performance, and present initial validations from commissioning observations. Comparisons between EPIC measurements and simultaneous beam-formed observations of bright sources show spectral-temporal structures in good agreement. △ Less

Submitted 23 January, 2023; originally announced January 2023.

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

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

Evaluating the prevalence of spurious correlations in pulsar timing array datasets

Authors: Andrew Zic, George Hobbs, R. M. Shannon, Daniel Reardon, Boris Goncharov, N. D. Ramesh Bhat, Andrew Cameron, Shi Dai, J. R. Dawson, Matthew Kerr, R. N. Manchester, Rami Mandow, Tommy Marshman, Christopher J. Russell, Nithyanandan Thyagarajan, X. -J. Zhu

Abstract: Pulsar timing array collaborations have recently reported evidence for a noise process with a common spectrum among the millisecond pulsars in the arrays. The spectral properties of this common-noise process are consistent with expectations for an isotropic gravitational-wave background (GWB) from inspiralling supermassive black-hole binaries. However, recent simulation analyses based on Parkes Pu… ▽ More Pulsar timing array collaborations have recently reported evidence for a noise process with a common spectrum among the millisecond pulsars in the arrays. The spectral properties of this common-noise process are consistent with expectations for an isotropic gravitational-wave background (GWB) from inspiralling supermassive black-hole binaries. However, recent simulation analyses based on Parkes Pulsar Timing Array data indicate that such a detection may arise spuriously. In this paper, we use simulated pulsar timing array datasets to further test the robustness of the inference methods for spectral and spatial correlations from a GWB. Expanding on our previous results, we find strong support (Bayes factors exceeding $10^5$) for the presence of a common-spectrum noise process in datasets where no common process is present, under a wide range of timing noise prescriptions per pulsar. We show that these results are highly sensitive to the choice of Bayesian priors on timing noise parameters, with priors that more closely match the injected distributions of timing noise parameters resulting in diminished support for a common-spectrum noise process. These results emphasize shortcomings in current methods for inferring the presence of a common-spectrum process, and imply that the detection of a common process is not a reliable precursor to detection of the GWB. Future searches for the nanohertz GWB should remain focussed on detecting spatial correlations, and make use of more tailored specifications for a common-spectrum noise process. △ Less

Submitted 25 July, 2022; originally announced July 2022.

Comments: 12 pages, 12 figures. Accepted for publication in MNRAS

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

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.

High-precision search for dark photon dark matter with the Parkes Pulsar Timing Array

Authors: Xiao Xue, Zi-Qing Xia, Xingjiang Zhu, Yue Zhao, Jing Shu, Qiang Yuan, N. D. Ramesh Bhat, Andrew D. Cameron, Shi Dai, Yi Feng, Boris Goncharov, George Hobbs, Eric Howard, Richard N. Manchester, Aditya Parthasarathy, Daniel J. Reardon, Christopher J. Russell, Ryan M. Shannon, Renée Spiewak, Nithyanandan Thyagarajan, Jingbo Wang, Lei Zhang, Songbo Zhang

Abstract: The nature of dark matter remains obscure in spite of decades of experimental efforts. The mass of dark matter candidates can span a wide range, and its coupling with the Standard Model sector remains uncertain. All these unknowns make the etection of dark matter extremely challenging. Ultralight dark matter, with $m \sim10^{-22}$ eV, is proposed to reconcile the disagreements between observations… ▽ More The nature of dark matter remains obscure in spite of decades of experimental efforts. The mass of dark matter candidates can span a wide range, and its coupling with the Standard Model sector remains uncertain. All these unknowns make the etection of dark matter extremely challenging. Ultralight dark matter, with $m \sim10^{-22}$ eV, is proposed to reconcile the disagreements between observations and predictions from simulations of small-scale structures in the cold dark matter paradigm, while remaining consistent with other observations. Because of its large de Broglie wavelength and large local occupation number within galaxies, ultralight dark matter behaves like a coherently oscillating background field with an oscillating frequency dependent on its mass. If the dark matter particle is a spin-1 dark photon, such as the $U(1)_B$ or $U(1)_{B-L}$ gauge boson, it can induce an external oscillating force and lead to displacements of test masses. Such an effect would be observable in the form of periodic variations in the arrival times of radio pulses from highly stable millisecond pulsars. In this study, we search for evidence of ultralight dark photon dark matter (DPDM) using 14-year high-precision observations of 26 pulsars collected with the Parkes Pulsar Timing Array. While no statistically significant signal is found, we place constraints on coupling constants for the $U(1)_B$ and $U(1)_{B-L}$ DPDM. Compared with other experiments, the limits on the dimensionless coupling constant $ε$ achieved in our study are improved by up to two orders of magnitude when the dark photon mass is smaller than $3\times10^{-22}$~eV ($10^{-22}$~eV) for the $U(1)_{B}$ ($U(1)_{B-L}$) scenario. △ Less

Submitted 26 December, 2021; v1 submitted 14 December, 2021; originally announced December 2021.

Comments: 12 pages, 4 figures

Journal ref: Phys. Rev. Research 4, L012022 (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

Constraining cosmological phase transitions with the Parkes Pulsar Timing Array

Authors: Xiao Xue, Ligong Bian, Jing Shu, Qiang Yuan, Xingjiang Zhu, N. D. Ramesh Bhat, Shi Dai, Yi Feng, Boris Goncharov, George Hobbs, Eric Howard, Richard N. Manchester, Christopher J. Russell, Daniel J. Reardon, Ryan M. Shannon, Renée Spiewak, Nithyanandan Thyagarajan, Jingbo Wang

Abstract: A cosmological first-order phase transition is expected to produce a stochastic gravitational wave background. If the phase transition temperature is on the MeV scale, the power spectrum of the induced stochastic gravitational waves peaks around nanohertz frequencies, and can thus be probed with high-precision pulsar timing observations. We search for such a stochastic gravitational wave backgroun… ▽ More A cosmological first-order phase transition is expected to produce a stochastic gravitational wave background. If the phase transition temperature is on the MeV scale, the power spectrum of the induced stochastic gravitational waves peaks around nanohertz frequencies, and can thus be probed with high-precision pulsar timing observations. We search for such a stochastic gravitational wave background with the latest data set of the Parkes Pulsar Timing Array. We find no evidence for a Hellings-Downs spatial correlation as expected for a stochastic gravitational wave background. Therefore, we present constraints on first-order phase transition model parameters. Our analysis shows that pulsar timing is particularly sensitive to the low-temperature ($T \sim 1 - 100$ MeV) phase transition with a duration $(β/H_*)^{-1}\sim 10^{-2}-10^{-1}$ and therefore can be used to constrain the dark and QCD phase transitions. △ Less

Submitted 25 November, 2022; v1 submitted 6 October, 2021; originally announced October 2021.

Comments: 9 pages, 4 figures, 2 tables

Journal ref: Phys. Rev. Lett. 127, 251303, Published 15 December 2021

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

Invariants in Polarimetric Interferometry: a non-Abelian Gauge Theory

Authors: Joseph Samuel, Rajaram Nityananda, Nithyanandan Thyagarajan

Abstract: The discovery of magnetic fields close to the M87 black hole using Very Long Baseline Interferometry (VLBI) by the Event Horizon Telescope collaboration utilized the novel concept of "closure traces", that are immune to element-based aberrations. We take a fundamentally new approach to this promising tool of polarimetric VLBI, using ideas from the geometric phase and gauge theories. The multiplica… ▽ More The discovery of magnetic fields close to the M87 black hole using Very Long Baseline Interferometry (VLBI) by the Event Horizon Telescope collaboration utilized the novel concept of "closure traces", that are immune to element-based aberrations. We take a fundamentally new approach to this promising tool of polarimetric VLBI, using ideas from the geometric phase and gauge theories. The multiplicative distortion of polarized signals at the individual elements are represented as gauge transformations by general $2\times 2$ complex matrices, so the closure traces now appear as gauge-invariant quantities. We apply this formalism to polarimetric interferometry and generalize it to any number of interferometer elements. Our approach goes beyond existing studies in the following respects: (1) we use triangular combinations of correlations as basic building blocks of invariants, (2) we use well-known symmetry properties of the Lorentz group to transparently identify a complete and independent set of invariants, and (3) we do not need auto-correlations, which are susceptible to large systematic biases, and therefore unreliable. This set contains all the information, immune to corruption, available in the interferometer measurements, thus providing important robust constraints for interferometric studies. △ Less

Submitted 1 February, 2022; v1 submitted 25 August, 2021; originally announced August 2021.

Comments: 8 pages (including references), 1 figure, and 5 appendices. Accepted (in press) in Physical Review Letters. Contains appendices and some text not included in the journal version that give more connections to gauge theory and geometric/Pancharatnam phase. See also companion paper arXiv:2108.11399

Journal ref: Physical Review Letters, 128, 091101 (2022)

Invariants in Co-polar Interferometry: an Abelian Gauge Theory

Authors: Nithyanandan Thyagarajan, Rajaram Nityananda, Joseph Samuel

Abstract: An $N$-element interferometer measures correlations among pairs of array elements. Closure invariants associated with closed loops among array elements are immune to multiplicative, element-based ("local") corruptions that occur in these measurements. Till recently, it has been unclear how a complete set of independent invariants can be analytically determined. We view the local, element-based cor… ▽ More An $N$-element interferometer measures correlations among pairs of array elements. Closure invariants associated with closed loops among array elements are immune to multiplicative, element-based ("local") corruptions that occur in these measurements. Till recently, it has been unclear how a complete set of independent invariants can be analytically determined. We view the local, element-based corruptions in co-polar correlations as gauge tranformations belonging to the gauge group $\textrm{GL}(1,\mathbb{C})$. Closure quantities are then naturally gauge invariant. We use this to provide a simple and effective formalism, and identify the complete set of independent closure invariants from co-polar interferometric correlations using only quantities defined on $(N-1)(N-2)/2$ elementary and independent triangular loops. The $(N-1)(N-2)/2$ closure phases and $N(N-3)/2$ closure amplitudes (totaling $N^2-3N+1$ real invariants), familiar in astronomical interferometry, naturally emerge from this formalism, which unifies what has required separate treatments until now. We do not require auto-correlations, but can easily include them if reliably measured. This unified view clarifies issues relating to noise and inference of object model parameters. It also allows us to extend the rule of parallel transport associated with Pancharatnam phase in optics to apply to amplitudes as well. The framework presented here extends to $\textrm{GL}(2,\mathbb{C}$) for full polarimetric interferometry as presented in a companion paper, which generalizes and clarifies earlier work. Our findings are relevant to state of the art co-polar and full polarimetric very long baseline interferometry measurements to determine features very near the event horizons of blackholes at the centers of M87, Centaurus~A, and the Milky Way. △ Less

Submitted 1 February, 2022; v1 submitted 25 August, 2021; originally announced August 2021.

Comments: 13 pages (including references), 0 figures, 5 appendices, accepted (in press) in Physical Review D. See also companion paper arXiv:2108.11400

Journal ref: Physical Review D, 105, 043019 (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

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/

On the evidence for a common-spectrum process in the search for the nanohertz gravitational-wave background with the Parkes Pulsar Timing Array

Authors: Boris Goncharov, R. M. Shannon, D. J. Reardon, G. Hobbs, A. Zic, M. Bailes, M. Curylo, S. Dai, M. Kerr, M. E. Lower, R. N. Manchester, R. Mandow, H. Middleton, M. T. Miles, A. Parthasarathy, E. Thrane, N. Thyagarajan, X. Xue, X. J. Zhu, A. D. Cameron, Y. Feng, R. Luo, C. J. Russell, J. Sarkissian, R. Spiewak , et al. (4 additional authors not shown)

Abstract: A nanohertz-frequency stochastic gravitational-wave background can potentially be detected through the precise timing of an array of millisecond pulsars. This background produces low-frequency noise in the pulse arrival times that would have a characteristic spectrum common to all pulsars and a well-defined spatial correlation. Recently the North American Nanohertz Observatory for Gravitational Wa… ▽ More A nanohertz-frequency stochastic gravitational-wave background can potentially be detected through the precise timing of an array of millisecond pulsars. This background produces low-frequency noise in the pulse arrival times that would have a characteristic spectrum common to all pulsars and a well-defined spatial correlation. Recently the North American Nanohertz Observatory for Gravitational Waves collaboration (NANOGrav) found evidence for the common-spectrum component in their 12.5-year data set. Here we report on a search for the background using the second data release of the Parkes Pulsar Timing Array. If we are forced to choose between the two NANOGrav models $\unicode{x2014}$ one with a common-spectrum process and one without $\unicode{x2014}$ we find strong support for the common-spectrum process. However, in this paper, we consider the possibility that the analysis suffers from model misspecification. In particular, we present simulated data sets that contain noise with distinctive spectra but show strong evidence for a common-spectrum process under the standard assumptions. The Parkes data show no significant evidence for, or against, the spatially correlated Hellings-Downs signature of the gravitational-wave background. Assuming we did observe the process underlying the spatially uncorrelated component of the background, we infer its amplitude to be $A = 2.2^{+0.4}_{-0.3} \times 10^{-15}$ in units of gravitational-wave strain at a frequency of $1\, \text{yr}^{-1}$. Extensions and combinations of existing and new data sets will improve the prospects of identifying spatial correlations that are necessary to claim a detection of the gravitational-wave background. △ Less

Submitted 11 August, 2021; v1 submitted 26 July, 2021; originally announced July 2021.

Comments: 10 pages, 5 figures

Journal ref: The Astrophysical Journal Letters, Volume 917, Number 2 (2021)

The Parkes pulsar timing array second data release: Timing analysis

Authors: D. J. Reardon, R. M. Shannon, A. D. Cameron, B. Goncharov, G. B. Hobbs, H. Middleton, M. Shamohammadi, N. Thyagarajan, M. Bailes, N. D. R. Bhat, S. Dai, M. Kerr, R. N. Manchester, C. J. Russell, R. Spiewak, J. B. Wang, X. J. Zhu

Abstract: The main goal of pulsar timing array experiments is to detect correlated signals such as nanohertz-frequency gravitational waves. Pulsar timing data collected in dense monitoring campaigns can also be used to study the stars themselves, their binary companions, and the intervening ionised interstellar medium. Timing observations are extraordinarily sensitive to changes in path length between the p… ▽ More The main goal of pulsar timing array experiments is to detect correlated signals such as nanohertz-frequency gravitational waves. Pulsar timing data collected in dense monitoring campaigns can also be used to study the stars themselves, their binary companions, and the intervening ionised interstellar medium. Timing observations are extraordinarily sensitive to changes in path length between the pulsar and the Earth, enabling precise measurements of the pulsar positions, distances and velocities, and the shapes of their orbits. Here we present a timing analysis of 25 pulsars observed as part of the Parkes Pulsar Timing Array (PPTA) project over time spans of up to 24 yr. The data are from the second data release of the PPTA, which we have extended by including legacy data. We make the first detection of Shapiro delay in four Southern pulsars (PSRs J1017$-$7156, J1125$-$6014, J1545$-$4550, and J1732$-$5049), and of parallax in six pulsars. The prominent Shapiro delay of PSR J1125$-$6014 implies a neutron star mass of $M_p = 1.5 \pm 0.2 M_\odot$ (68% credibility interval). Measurements of both Shapiro delay and relativistic periastron advance in PSR J1600$-$3053 yield a large but uncertain pulsar mass of $M_p = 2.06^{+0.44}_{-0.41}$ M$_\odot$ (68% credibility interval). We measure the distance to PSR J1909$-$3744 to a precision of 10 lyr, indicating that for gravitational wave periods over a decade, the pulsar provides a coherent baseline for pulsar timing array experiments. △ Less

Submitted 16 July, 2021; v1 submitted 9 July, 2021; originally announced July 2021.

Comments: 44 pages, 32 figures, 7 tables. Accepted for publication in MNRAS; v2 includes corrected typos in parameter tables and reference added

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

A Geometric View of Closure Phases in Interferometry

Authors: Nithyanandan Thyagarajan, Christopher L. Carilli

Abstract: Closure phase is the phase of a closed-loop product of correlations in a $\ge 3$-element interferometer array. Its invariance to element-based phase corruption makes it invaluable for interferometric applications that otherwise require high-accuracy phase calibration. However, its understanding has remained mainly mathematical and limited to the aperture plane (Fourier dual of image plane). Here,… ▽ More Closure phase is the phase of a closed-loop product of correlations in a $\ge 3$-element interferometer array. Its invariance to element-based phase corruption makes it invaluable for interferometric applications that otherwise require high-accuracy phase calibration. However, its understanding has remained mainly mathematical and limited to the aperture plane (Fourier dual of image plane). Here, we lay the foundations for a geometrical insight. we show that closure phase and its invariance to element-based corruption and to translation are intricately related to the conserved properties (shape, orientation, and size, or SOS) of the principal triangle enclosed by the three fringes formed by a closed triad of array elements, which is referred herein as the "SOS conservation principle". When element-based amplitude calibration is not needed, as is typical in optical interferometry, the 3-element interference image formed from phase-uncalibrated correlations is a true and uncorrupted representation of the source object's morphology, except for a possible shift. Based on this SOS conservation principle, we present two geometric methods to measure the closure phase directly from a 3-element interference image (without requiring an aperture-plane view): (i) the closure phase is directly measurable from any one of the triangle's heights, and (ii) the squared closure phase is proportional to the product of the areas enclosed by the triad of array elements and the principal triangle in the aperture and image planes, respectively. We validate this geometric understanding across a wide range range of interferometric conditions using data from the Very Large Array and the Event Horizon Telescope. This geometric insight can be potentially valuable to other interferometric applications such as optical interferometry. These geometric relationships are generalised for an $N$-element interferometer. △ Less

Submitted 25 February, 2022; v1 submitted 9 December, 2020; originally announced December 2020.

Comments: Accepted in PASA. 18 pages including 3 appendices, 11 captioned figures. Abstract abridged for arXiv. Removed some previous appendices. Parallels of closure phases in crystallography, interferometric synthetic aperture radar, seismic imaging, gravitational wave time-delay interferometry, and Bargmann invariants in quantum mechanics, described in appendices of arXiv:2012.05254v3

Journal ref: Publications of the Astronomical Society of Australia, 39, e014 (2022)

Searching for gravitational wave bursts from cosmic string cusps with the Parkes Pulsar Timing Array

Authors: N. Yonemaru, S. Kuroyanagi, G. Hobbs, K. Takahashi, X. -J. Zhu, W. A. Coles, S. Dai, E. Howard, R. Manchester, D. Reardon, C. Russell, R. Shannon, N. Thyagarajan, R. Spiewak, J. -B. Wang

Abstract: Cosmic strings are potential gravitational wave (GW) sources that can be probed by pulsar timing arrays (PTAs). In this work we develop a detection algorithm for a GW burst from a cusp on a cosmic string, and apply it to Parkes PTA data. We find four events with a false alarm probability less than 1%. However further investigation shows that all of these are likely to be spurious. As there are no… ▽ More Cosmic strings are potential gravitational wave (GW) sources that can be probed by pulsar timing arrays (PTAs). In this work we develop a detection algorithm for a GW burst from a cusp on a cosmic string, and apply it to Parkes PTA data. We find four events with a false alarm probability less than 1%. However further investigation shows that all of these are likely to be spurious. As there are no convincing detections we place upper limits on the GW amplitude for different event durations. From these bounds we place limits on the cosmic string tension of G mu ~ 10^{-5}, and highlight that this bound is independent from those obtained using other techniques. We discuss the physical implications of our results and the prospect of probing cosmic strings in the era of Square Kilometre Array (SKA). △ Less

Submitted 26 November, 2020; originally announced November 2020.

Comments: Accepted for publication by MNRAS

Identifying and mitigating noise sources in precision pulsar timing data sets

Authors: Boris Goncharov, D. J. Reardon, R. M. Shannon, Xing-Jiang Zhu, Eric Thrane, M. Bailes, N. D. R. Bhat, S. Dai, G. Hobbs, M. Kerr, R. N. Manchester, S. Osłowski, A. Parthasarathy, C. J. Russell, R. Spiewak, N. Thyagarajan, J. B. Wang

Abstract: Pulsar timing array projects measure the pulse arrival times of millisecond pulsars for the primary purpose of detecting nanohertz-frequency gravitational waves. The measurements include contributions from a number of astrophysical and instrumental processes, which can either be deterministic or stochastic. It is necessary to develop robust statistical and physical models for these noise processes… ▽ More Pulsar timing array projects measure the pulse arrival times of millisecond pulsars for the primary purpose of detecting nanohertz-frequency gravitational waves. The measurements include contributions from a number of astrophysical and instrumental processes, which can either be deterministic or stochastic. It is necessary to develop robust statistical and physical models for these noise processes because incorrect models diminish sensitivity and may cause a spurious gravitational wave detection. Here we characterise noise processes for the 26 pulsars in the second data release of the Parkes Pulsar Timing Array using Bayesian inference. In addition to well-studied noise sources found previously in pulsar timing array data sets such as achromatic timing noise and dispersion measure variations, we identify new noise sources including time-correlated chromatic noise that we attribute to variations in pulse scattering. We also identify "exponential dip" events in four pulsars, which we attribute to magnetospheric effects as evidenced by pulse profile shape changes observed for three of the pulsars. This includes an event in PSR J1713$+$0747, which had previously been attributed to interstellar propagation. We present noise models to be used in searches for gravitational waves. We outline a robust methodology to evaluate the performance of noise models and identify unknown signals in the data. The detection of variations in pulse profiles highlights the need to develop efficient profile domain timing methods. △ Less

Submitted 1 November, 2020; v1 submitted 12 October, 2020; originally announced October 2020.

Comments: 18 pages, 7 figures, 5 tables

Extremely band-limited repetition from a fast radio burst source

Authors: Pravir Kumar, Ryan M. Shannon, Chris Flynn, Stefan Osłowski, Shivani Bhandari, Cherie K. Day, Adam T. Deller, Wael Farah, Jane F. Kaczmarek, Matthew Kerr, Chris Phillips, Danny C. Price, Hao Qiu, Nithyanandan Thyagarajan

Abstract: The fast radio burst (FRB) population is observationally divided into sources that have been observed to repeat and those that have not. There is tentative evidence that the bursts from repeating sources have different properties than the non-repeating ones. In order to determine the occurrence rate of repeating sources and characterize the nature of repeat emission, we have been conducting sensit… ▽ More The fast radio burst (FRB) population is observationally divided into sources that have been observed to repeat and those that have not. There is tentative evidence that the bursts from repeating sources have different properties than the non-repeating ones. In order to determine the occurrence rate of repeating sources and characterize the nature of repeat emission, we have been conducting sensitive searches for repetitions from bursts detected with the Australian Square Kilometre Array Pathfinder (ASKAP) with the 64-m Parkes radio telescope, using the recently commissioned Ultra-wideband Low (UWL) receiver system, over a band spanning 0.7$-$4.0 GHz. We report the detection of a repeat burst from the source of FRB 20190711A. The detected burst is 1 ms wide and has a bandwidth of just 65 MHz. We find no evidence of any emission in the remaining part of the 3.3 GHz UWL band. While the emission bandwidths of the ASKAP and UWL bursts show $ν^{-4}$ scaling consistent with a propagation effect, the spectral occupancy is inconsistent with diffractive scintillation. This detection rules out models predicting broad-band emission from the FRB 20190711A source and puts stringent constraints on the emission mechanism. The low spectral occupancy highlights the importance of sub-banded search methods in detecting FRBs. △ Less

Submitted 7 December, 2020; v1 submitted 2 September, 2020; originally announced September 2020.

Comments: 7 pages, 5 figures; accepted for publication in MNRAS

Statistical Detection of IGM Structures during Cosmic Reionization using Absorption of the Redshifted 21 cm line by HI against Compact Background Radio Sources

Authors: Nithyanandan Thyagarajan

Abstract: Detecting neutral hydrogen structures in the intergalactic medium (IGM) during cosmic reionization via absorption (21 cm forest) against a background radiation is considered independent and complementary to the three-dimensional tomography and power spectrum techniques. The direct detection of this absorption requires very bright ($\gtrsim 10$-100 mJy) background sources at high redshifts (… ▽ More Detecting neutral hydrogen structures in the intergalactic medium (IGM) during cosmic reionization via absorption (21 cm forest) against a background radiation is considered independent and complementary to the three-dimensional tomography and power spectrum techniques. The direct detection of this absorption requires very bright ($\gtrsim 10$-100 mJy) background sources at high redshifts ($z\gtrsim 8$), which are evidently rare; very long times of integration; or instruments of very high sensitivity. This motivates a statistical one-dimensional (1D) power spectrum approach along narrow sightlines but with fainter background objects ($\sim 1$-10 mJy), which are likely to be more abundant and significant contributors at high redshifts. The 1D power spectrum reduces cosmic variance and improves sensitivity especially on small spatial scales. Using standard radiative transfer and fiducial models for the instrument, the background sources, and the evolution of IGM structures during cosmic reionization, the potential of the 1D power spectrum along selected narrow directions is investigated against uncertainties from thermal noise and the chromatic synthesized point spread function (PSF) response. Minimum requirements on the number of high-redshift background sources, the telescope sensitivity, and the PSF quality are estimated for a range of instrumental, background source, and reionization model parameters. The 1D power spectrum is intrinsically stronger at higher redshifts. A $\sim 1000$ hr observing campaign targeting $\sim 100$ narrow sightlines to radio-faint, high-redshift background objects with modern radio telescopes, especially the Square Kilometre Array, can detect the 1D power spectrum on a range of spatial scales and redshifts, and potentially discriminate between models of cosmic reionization. △ Less

Submitted 7 August, 2020; v1 submitted 17 June, 2020; originally announced June 2020.

Comments: 27 pages, 18 figures (including subfigures), published in the Astrophysical Journal

Journal ref: The Astrophysical Journal, 899, 16 (2020)

Measuring HERA's primary beam in-situ: methodology and first results

Authors: Chuneeta D. Nunhokee, Aaron R. Parsons, Nicholas S. Kern, Bojan Nikolic, Jonathan C. Pober, Gianni Bernardi, Chris L. Carilli, Zara Abdurashidova, James E. Aguirre, Paul Alexander, Zaki S. Ali, Yanga Balfour, Adam P. Beardsley, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Jacob Burba, Carina Cheng, David R. DeBoer, Matt Dexter, Eloy de~Lera~Acedo, Joshua S. Dillon, Aaron Ewall-Wice, Nicolas Fagnoni, Randall Fritz , et al. (42 additional authors not shown)

Abstract: The central challenge in 21~cm cosmology is isolating the cosmological signal from bright foregrounds. Many separation techniques rely on the accurate knowledge of the sky and the instrumental response, including the antenna primary beam. For drift-scan telescopes such as the Hydrogen Epoch of Reionization Array \citep[HERA, ][]{DeBoer2017} that do not move, primary beam characterization is partic… ▽ More The central challenge in 21~cm cosmology is isolating the cosmological signal from bright foregrounds. Many separation techniques rely on the accurate knowledge of the sky and the instrumental response, including the antenna primary beam. For drift-scan telescopes such as the Hydrogen Epoch of Reionization Array \citep[HERA, ][]{DeBoer2017} that do not move, primary beam characterization is particularly challenging because standard beam-calibration routines do not apply \citep{Cornwell2005} and current techniques require accurate source catalogs at the telescope resolution. We present an extension of the method from \citet{Pober2012} where they use beam symmetries to create a network of overlapping source tracks that break the degeneracy between source flux density and beam response and allow their simultaneous estimation. We fit the beam response of our instrument using early HERA observations and find that our results agree well with electromagnetic simulations down to a -20~dB level in power relative to peak gain for sources with high signal-to-noise ratio. In addition, we construct a source catalog with 90 sources down to a flux density of 1.4~Jy at 151~MHz. △ Less

Submitted 25 May, 2020; originally announced May 2020.

Comments: 22 pages, 22 figures, accepted for publication in ApJ

Detection of Cosmic Structures using the Bispectrum Phase. II. First Results from Application to Cosmic Reionization Using the Hydrogen Epoch of Reionization Array

Authors: Nithyanandan Thyagarajan, Chris L. Carilli, Bojan Nikolic, James Kent, Andrei Mesinger, Nicholas S. Kern, Gianni Bernardi, Siyanda Matika, Zara Abdurashidova, James E. Aguirre, Paul Alexander, Zaki S. Ali, Yanga Balfour, Adam P. Beardsley, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Jacob Burba, Steve Carey, Carina Cheng, David R. DeBoer, Matt Dexter, Eloy de Lera Acedo, Joshua S. Dillon, John Ely , et al. (47 additional authors not shown)

Abstract: Characterizing the epoch of reionization (EoR) at $z\gtrsim 6$ via the redshifted 21 cm line of neutral Hydrogen (HI) is critical to modern astrophysics and cosmology, and thus a key science goal of many current and planned low-frequency radio telescopes. The primary challenge to detecting this signal is the overwhelmingly bright foreground emission at these frequencies, placing stringent requirem… ▽ More Characterizing the epoch of reionization (EoR) at $z\gtrsim 6$ via the redshifted 21 cm line of neutral Hydrogen (HI) is critical to modern astrophysics and cosmology, and thus a key science goal of many current and planned low-frequency radio telescopes. The primary challenge to detecting this signal is the overwhelmingly bright foreground emission at these frequencies, placing stringent requirements on the knowledge of the instruments and inaccuracies in analyses. Results from these experiments have largely been limited not by thermal sensitivity but by systematics, particularly caused by the inability to calibrate the instrument to high accuracy. The interferometric bispectrum phase is immune to antenna-based calibration and errors therein, and presents an independent alternative to detect the EoR HI fluctuations while largely avoiding calibration systematics. Here, we provide a demonstration of this technique on a subset of data from the Hydrogen Epoch of Reionization Array (HERA) to place approximate constraints on the brightness temperature of the intergalactic medium (IGM). From this limited data, at $z=7.7$ we infer "$1σ$" upper limits on the IGM brightness temperature to be $\le 316$ "pseudo" mK at $κ_\parallel=0.33$ "pseudo" $h$ Mpc$^{-1}$ (data-limited) and $\le 1000$ "pseudo" mK at $κ_\parallel=0.875$ "pseudo" $h$ Mpc$^{-1}$ (noise-limited). The "pseudo" units denote only an approximate and not an exact correspondence to the actual distance scales and brightness temperatures. By propagating models in parallel to the data analysis, we confirm that the dynamic range required to separate the cosmic HI signal from the foregrounds is similar to that in standard approaches, and the power spectrum of the bispectrum phase is still data-limited (at $\gtrsim 10^6$ dynamic range) indicating scope for further improvement in sensitivity as the array build-out continues. △ Less

Submitted 2 July, 2020; v1 submitted 20 May, 2020; originally announced May 2020.

Comments: 22 pages, 12 figures (including sub-figures). Published in PhRvD. Abstract may be slightly abridged compared to the actual manuscript due to length limitations on arXiv

Journal ref: Phys. Rev. D 102, 022002 (2020)

Detection of Cosmic Structures using the Bispectrum Phase. I. Mathematical Foundations

Authors: Nithyanandan Thyagarajan, Chris Carilli

Abstract: Many low-frequency radio interferometers are aiming to detect very faint spectral signatures from structures at cosmological redshifts, particularly of neutral Hydrogen using its characteristic 21 cm spectral line. Due to the very high dynamic range needed to isolate these faint spectral fluctuations from the very bright foregrounds, spectral systematics from the instrument or the analysis, rather… ▽ More Many low-frequency radio interferometers are aiming to detect very faint spectral signatures from structures at cosmological redshifts, particularly of neutral Hydrogen using its characteristic 21 cm spectral line. Due to the very high dynamic range needed to isolate these faint spectral fluctuations from the very bright foregrounds, spectral systematics from the instrument or the analysis, rather than thermal noise, are currently limiting their sensitivity. Failure to achieve a spectral calibration with fractional inaccuracy $\lesssim 10^{-5}$ will make the detection of the critical cosmic signal unlikely. The bispectrum phase from interferometric measurements is largely immune to this calibration issue. We present a basis to explore the nature of bispectrum phase in the limit of small spectral fluctuations. We establish that they measure the intrinsic dissimilarity in the transverse structure of the cosmic signal relative to the foregrounds, expressed as rotations in the underlying phase angle. Their magnitude is related to the strength of the cosmic signal relative to the foregrounds. Using a range of sky models, we detail the behavior of bispectrum phase fluctuations using standard Fourier-domain techniques and find it comparable to existing approaches, with a few key differences. Mode-mixed foreground contamination is more pronounced than in existing approaches because the bispectrum phase is a product of three individual interferometric phases. The multiplicative coupling of foregrounds in the bispectrum phase fluctuations results in the mixing of foreground signatures with that of the cosmic signal. We briefly outline a variation of this approach to avoid extensive mode-mixing. Despite its limitations, the interpretation of results using bispectrum phase is possible with forward-modeling. Importantly, it is an independent and a viable alternative to existing approaches. △ Less

Submitted 2 July, 2020; v1 submitted 20 May, 2020; originally announced May 2020.

Comments: 27 pages, 22 figures (incl. sub-figures), published in PhRvD. Abstract may be slightly abridged compared to the actual manuscript due to length limitations on arXiv

Journal ref: Phys. Rev. D 102, 022001 (2020)

Foreground modelling via Gaussian process regression: an application to HERA data

Authors: Abhik Ghosh, Florent Mertens, Gianni Bernardi, Mário G. Santos, Nicholas S. Kern, Christopher L. Carilli, Trienko L. Grobler, Léon V. E. Koopmans, Daniel C. Jacobs, Adrian Liu, Aaron R. Parsons, Miguel F. Morales, James E. Aguirre, Joshua S. Dillon, Bryna J. Hazelton, Oleg M. Smirnov, Bharat K. Gehlot, Siyanda Matika, Paul Alexander, Zaki S. Ali, Adam P. Beardsley, Roshan K. Benefo, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley , et al. (48 additional authors not shown)

Abstract: The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize… ▽ More The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize the measurement of the 21-cm signal. In this paper, we use Gaussian Process Regression to model both foreground emission and instrumental systematics in $\sim 2$ hours of data from the Hydrogen Epoch of Reionization Array. We find that a simple co-variance model with three components matches the data well, giving a residual power spectrum with white noise properties. These consist of an "intrinsic" and instrumentally corrupted component with a coherence-scale of 20 MHz and 2.4 MHz respectively (dominating the line of sight power spectrum over scales $k_{\parallel} \le 0.2$ h cMpc$^{-1}$) and a baseline dependent periodic signal with a period of $\sim 1$ MHz (dominating over $k_{\parallel} \sim 0.4 - 0.8$h cMpc$^{-1}$) which should be distinguishable from the 21-cm EoR signal whose typical coherence-scales is $\sim 0.8$ MHz. △ Less

Submitted 12 May, 2020; v1 submitted 13 April, 2020; originally announced April 2020.

Comments: 15 pages, 15 figures, 1 table, Accepted to MNRAS

Redundant-Baseline Calibration of the Hydrogen Epoch of Reionization Array

Authors: Joshua S. Dillon, Max Lee, Zaki S. Ali, Aaron R. Parsons, Naomi Orosz, Chuneeta Devi Nunhokee, Paul La Plante, Adam P. Beardsley, Nicholas S. Kern, 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 , et al. (54 additional authors not shown)

Abstract: In 21 cm cosmology, precision calibration is key to the separation of the neutral hydrogen signal from very bright but spectrally-smooth astrophysical foregrounds. The Hydrogen Epoch of Reionization Array (HERA), an interferometer specialized for 21 cm cosmology and now under construction in South Africa, was designed to be largely calibrated using the self-consistency of repeated measurements of… ▽ More In 21 cm cosmology, precision calibration is key to the separation of the neutral hydrogen signal from very bright but spectrally-smooth astrophysical foregrounds. The Hydrogen Epoch of Reionization Array (HERA), an interferometer specialized for 21 cm cosmology and now under construction in South Africa, was designed to be largely calibrated using the self-consistency of repeated measurements of the same interferometric modes. This technique, known as "redundant-baseline calibration" resolves most of the internal degrees of freedom in the calibration problem. It assumes, however, on antenna elements with identical primary beams placed precisely on a redundant grid. In this work, we review the detailed implementation of the algorithms enabling redundant-baseline calibration and report results with HERA data. We quantify the effects of real-world non-redundancy and how they compare to the idealized scenario in which redundant measurements differ only in their noise realizations. Finally, we study how non-redundancy can produce spurious temporal structure in our calibration solutions--both in data and in simulations--and present strategies for mitigating that structure. △ Less

Submitted 3 November, 2020; v1 submitted 18 March, 2020; originally announced March 2020.

Comments: 24 Pages, 19 Figures. Updated to match the accepted MNRAS version

The ALMA Spectroscopic Survey in the HUDF: Constraining cumulative CO emission at $1 \lesssim z \lesssim 4$ with power spectrum analysis of ASPECS LP data from 84 to 115 GHz

Authors: Bade Uzgil, Chris Carilli, Adam Lidz, Fabian Walter, Nithyanandan Thyagarajan, Roberto Decarli, Manuel Aravena, Frank Bertoldi, Paulo C. Cortes, Jorge González-López, Hanae Inami, Gergö Popping, Paul Van der Werf, Jeff Wagg, Axel Weiss

Abstract: We present a power spectrum analysis of the ALMA Spectroscopic Survey Large Program (ASPECS LP) data from 84 to 115 GHz. These data predominantly probe small-scale fluctuations ($k=10$-$100$ h Mpc$^{-1}$) in the aggregate CO emission in galaxies at $1 \lesssim z \lesssim 4$. We place an integral constraint on CO luminosity functions (LFs) in this redshift range via a direct measurement of their se… ▽ More We present a power spectrum analysis of the ALMA Spectroscopic Survey Large Program (ASPECS LP) data from 84 to 115 GHz. These data predominantly probe small-scale fluctuations ($k=10$-$100$ h Mpc$^{-1}$) in the aggregate CO emission in galaxies at $1 \lesssim z \lesssim 4$. We place an integral constraint on CO luminosity functions (LFs) in this redshift range via a direct measurement of their second moments in the three-dimensional (3D) auto-power spectrum, finding a total CO shot noise power $P_{\textrm{CO,CO}}(k_{\textrm{CO(2-1)}}) \leq 1.9\times10^2$ $μ$K$^2$ (Mpc h$^{-1}$)$^3$. This upper limit ($3σ$) is consistent with the observed ASPECS CO LFs in Decarli et al. 2019, but rules out a large space in the range of $P_{\textrm{CO,CO}}(k_{\textrm{CO(2-1)}})$ inferred from these LFs, which we attribute primarily to large uncertainties in the normalization $Φ_*$ and knee $L_*$ of the Schechter-form CO LFs at $z > 2$. Also, through power spectrum analyses of ASPECS LP data with 415 positions from galaxies with available optical spectroscopic redshifts, we find that contributions to the observed mean CO intensity and shot noise power of MUSE galaxies are largely accounted for by ASPECS blind detections, though there are $\sim20$% contributions to the CO(2-1) mean intensity due to sources previously undetected in the blind line search. Finally, we sum the fluxes from individual blind CO detections to yield a lower limit on the mean CO surface brightness at 99 GHz of $\langle T_{\textrm{CO}} \rangle = 0.55\pm0.02$ $μ$K, which we estimate represents $68$-$80$% of the total CO surface brightness at this frequency. △ Less

Submitted 31 October, 2019; originally announced November 2019.

Comments: 31 pages, 8 Figures, accepted for publication in ApJ