Research Article
Variational inference as an alternative to MCMC for parameter estimation and model selection
- Geetakrishnasai Gunapati, Anirudh Jain, P. K. Srijith, Shantanu Desai
-
- Published online by Cambridge University Press:
- 25 January 2022, e001
-
Most applications of Bayesian Inference for parameter estimation and model selection in astrophysics involve the use of Monte Carlo techniques such as Markov Chain Monte Carlo (MCMC) and nested sampling. However, these techniques are time-consuming and their convergence to the posterior could be difficult to determine. In this study, we advocate variational inference as an alternative to solve the above problems, and demonstrate its usefulness for parameter estimation and model selection in astrophysics. Variational inference converts the inference problem into an optimisation problem by approximating the posterior from a known family of distributions and using Kullback–Leibler divergence to characterise the difference. It takes advantage of fast optimisation techniques, which make it ideal to deal with large datasets and makes it trivial to parallelise on a multicore platform. We also derive a new approximate evidence estimation based on variational posterior, and importance sampling technique called posterior-weighted importance sampling for the calculation of evidence, which is useful to perform Bayesian model selection. As a proof of principle, we apply variational inference to five different problems in astrophysics, where Monte Carlo techniques were previously used. These include assessment of significance of annual modulation in the COSINE-100 dark matter experiment, measuring exoplanet orbital parameters from radial velocity data, tests of periodicities in measurements of Newton’s constant G, assessing the significance of a turnover in the spectral lag data of GRB 160625B, and estimating the mass of a galaxy cluster using weak gravitational lensing. We find that variational inference is much faster than MCMC and nested sampling techniques for most of these problems while providing competitive results. All our analysis codes have been made publicly available.
Massive high-redshift quiescent galaxies with JWST
- Themiya Nanayakkara, James Esdaile, Karl Glazebrook, Juan M. Espejo Salcedo, Mark Durre, Colin Jacobs
-
- Published online by Cambridge University Press:
- 26 January 2022, e002
-
- Article
-
- You have access Access
- Open access
- HTML
- Export citation
- NASA ADS Abstract Service
-
Recent ground-based deep observations of the Universe have discovered large populations of massive quiescent galaxies at $z\sim3\!-\!5$ . With the launch of the James Webb Space Telescope (JWST), the on-board Near-Infrared Spectrograph (NIRSpec) instrument will provide continuous $0.6\!-\!5.3\,\unicode{x03BC}\,\mathrm{m}$ spectroscopic coverage of these galaxies. Here we show that NIRSpec/CLEAR spectroscopy is ideal to probe the completeness of photometrically selected massive quiescent galaxies such as the ones presented by Schreiber et al. (2018b, A&A, 618, A85). Using a subset of the Schreiber et al. (2018b, A&A, 618, A85) sample with deep Keck/MOSFIRE spectroscopy presented by Esdaile J., et al. (2021b, ApJ, 908, L35), we perform a suite of mock JWST/NIRSpec observations to determine optimal observing strategies to efficiently recover the star formation histories (SFHs), element abundances, and kinematics of these massive quiescent galaxies. We find that at $z\sim3$ , medium resolution G235M/FL170LP NIRSpec observations could recover element abundances at an accuracy of ${\sim}15\%$ , which is comparable to local globular clusters. Mimicking ZFOURGE COSMOS photometry, we perform mock spectrophotometric fitting with Prospector to show that the overall shape of the SFHs of our mock galaxies can be recovered well, albeit with a dependency on the number of non-parametric SFH bins. We show that deep high-resolution G235H/FL170LP integral field spectroscopy with a $S/N\sim7$ per spaxel is required to constrain the rotational properties of our sample at $>\!2\sigma$ confidence. Thus, through optimal grism/filter choices, JWST/NIRSpec slit and integral field spectroscopy observations would provide tight constraints to galaxy evolution in the early Universe.
Early-time searches for coherent radio emission from short GRBs with the Murchison Widefield Array
- J. Tian, G. E. Anderson, P. J. Hancock, J. C. A. Miller-Jones, M. Sokolowski, A. Rowlinson, A. Williams, J. Morgan, N. Hurley-Walker, D. L. Kaplan, Tara Murphy, S. J. Tingay, M. Johnston-Hollitt, K. W. Bannister, M. E. Bell, B. W. Meyers
-
- Published online by Cambridge University Press:
- 03 February 2022, e003
-
Many short gamma-ray bursts (GRBs) originate from binary neutron star mergers, and there are several theories that predict the production of coherent, prompt radio signals either prior, during, or shortly following the merger, as well as persistent pulsar-like emission from the spin-down of a magnetar remnant. Here we present a low frequency (170–200 MHz) search for coherent radio emission associated with nine short GRBs detected by the Swift and/or Fermi satellites using the Murchison Widefield Array (MWA) rapid-response observing mode. The MWA began observing these events within 30–60 s of their high-energy detection, enabling us to capture any dispersion delayed signals emitted by short GRBs for a typical range of redshifts. We conducted transient searches at the GRB positions on timescales of 5 s, 30 s, and 2 min, resulting in the most constraining flux density limits on any associated transient of 0.42, 0.29, and 0.084 Jy, respectively. We also searched for dispersed signals at a temporal and spectral resolution of 0.5 s and 1.28 MHz, but none were detected. However, the fluence limit of 80–100 Jy ms derived for GRB 190627A is the most stringent to date for a short GRB. Assuming the formation of a stable magnetar for this GRB, we compared the fluence and persistent emission limits to short GRB coherent emission models, placing constraints on key parameters including the radio emission efficiency of the nearly merged neutron stars ( $\epsilon_r\lesssim10^{-4}$ ), the fraction of magnetic energy in the GRB jet ( $\epsilon_B\lesssim2\times10^{-4}$ ), and the radio emission efficiency of the magnetar remnant ( $\epsilon_r\lesssim10^{-3}$ ). Comparing the limits derived for our full GRB sample (along with those in the literature) to the same emission models, we demonstrate that our fluence limits only place weak constraints on the prompt emission predicted from the interaction between the relativistic GRB jet and the interstellar medium for a subset of magnetar parameters. However, the 30-min flux density limits were sensitive enough to theoretically detect the persistent radio emission from magnetar remnants up to a redshift of $z\sim0.6$ . Our non-detection of this emission could imply that some GRBs in the sample were not genuinely short or did not result from a binary neutron star merger, the GRBs were at high redshifts, these mergers formed atypical magnetars, the radiation beams of the magnetar remnants were pointing away from Earth, or the majority did not form magnetars but rather collapse directly into black holes.
The non-thermal emission from the colliding-wind binary Apep
- S. del Palacio, P. Benaglia, M. De Becker, V. Bosch-Ramon, G. E. Romero
-
- Published online by Cambridge University Press:
- 31 January 2022, e004
-
Therecently discovered massive binary system Apep is the most powerful synchrotron emitter among the known Galactic colliding-wind binaries. This makes this particular system of great interest to investigate stellar winds and the non-thermal processes associated with their shocks. This source was detected at various radio bands, and in addition the wind-collision region was resolved by means of very-long baseline interferometric observations. We use a non-thermal emission model for colliding-wind binaries to derive physical properties of this system. The observed morphology in the resolved maps allows us to estimate the system projection angle on the sky to be $\psi \approx 85^\circ$ . The observed radio flux densities also allow us to characterise both the intrinsic synchrotron spectrum of the source and its modifications due to free–free absorption in the stellar winds at low frequencies; from this, we derive mass–loss rates of the stars of $\dot{M}_\mathrm{WN} \approx 4\times10^{-5}\;\mathrm{M}_\odot\,\mathrm{yr}^{-1}$ and $\dot{M}_\mathrm{WC} \approx 2.9\times10^{-5}\;\mathrm{M}_\odot\,\mathrm{yr}^{-1}$ . Finally, the broadband spectral energy distribution is calculated for different combinations of the remaining free parameters, namely the intensity of the magnetic field and the injected power in non-thermal particles. We show that the degeneracy of these two parameters can be solved with observations in the high-energy domain, most likely in the hard X-rays but also possibly in $\gamma$ -rays under favourable conditions.
GASKAP-HI pilot survey science I: ASKAP zoom observations of Hi emission in the Small Magellanic Cloud
- N. M. Pingel, J. Dempsey, N. M. McClure-Griffiths, J. M. Dickey, K. E. Jameson, H. Arce, G. Anglada, J. Bland-Hawthorn, S. L. Breen, F. Buckland-Willis, S. E. Clark, J. R. Dawson, H. Dénes, E. M. Di Teodoro, B.-Q. For, Tyler J. Foster, J. F. Gómez, H. Imai, G. Joncas, C.-G. Kim, M.-Y. Lee, C. Lynn, D. Leahy, Y. K. Ma, A. Marchal, D. McConnell, M.-A. Miville-Deschènes, V. A. Moss, C. E. Murray, D. Nidever, J. Peek, S. Stanimirović, L. Staveley-Smith, T. Tepper-Garcia, C. D. Tremblay, L. Uscanga, J. Th. van Loon, E. Vázquez-Semadeni, J. R. Allison, C. S. Anderson, Lewis Ball, M. Bell, D. C.-J. Bock, J. Bunton, F. R. Cooray, T. Cornwell, B. S. Koribalski, N. Gupta, D. B. Hayman, L. Harvey-Smith, K. Lee-Waddell, A. Ng, C. J. Phillips, M. Voronkov, T. Westmeier, M. T. Whiting
-
- Published online by Cambridge University Press:
- 07 February 2022, e005
-
- Article
-
- You have access Access
- Open access
- HTML
- Export citation
- NASA ADS Abstract Service
-
We present the most sensitive and detailed view of the neutral hydrogen ( ${\rm H\small I}$ ) emission associated with the Small Magellanic Cloud (SMC), through the combination of data from the Australian Square Kilometre Array Pathfinder (ASKAP) and Parkes (Murriyang), as part of the Galactic Australian Square Kilometre Array Pathfinder (GASKAP) pilot survey. These GASKAP-HI pilot observations, for the first time, reveal ${\rm H\small I}$ in the SMC on similar physical scales as other important tracers of the interstellar medium, such as molecular gas and dust. The resultant image cube possesses an rms noise level of 1.1 K ( $1.6\,\mathrm{mJy\ beam}^{-1}$ ) $\mathrm{per}\ 0.98\,\mathrm{km\ s}^{-1}$ spectral channel with an angular resolution of $30^{\prime\prime}$ ( ${\sim}10\,\mathrm{pc}$ ). We discuss the calibration scheme and the custom imaging pipeline that utilises a joint deconvolution approach, efficiently distributed across a computing cluster, to accurately recover the emission extending across the entire ${\sim}25\,\mathrm{deg}^2$ field-of-view. We provide an overview of the data products and characterise several aspects including the noise properties as a function of angular resolution and the represented spatial scales by deriving the global transfer function over the full spectral range. A preliminary spatial power spectrum analysis on individual spectral channels reveals that the power law nature of the density distribution extends down to scales of 10 pc. We highlight the scientific potential of these data by comparing the properties of an outflowing high-velocity cloud with previous ASKAP+Parkes ${\rm H\small I}$ test observations.
Gravitational clearing of natural satellite orbits
- Roderick J. Hill
-
- Published online by Cambridge University Press:
- 08 February 2022, e006
-
The distribution of diameters and orbital distances from the parent body of 156 named moons of the planets in the Solar System is not random. All 11 moons with diameters larger than $1\,000\,\mathrm{km}$ are positioned between $400\,000\,\mathrm{km}$ and 4 million km from the parent, whereas the far more numerous small moons are distributed on both sides of this central region and are largely absent from the region in between. This small-satellite ‘exclusion region’ is particularly evident for the gas giants since they have multiple satellites spanning a wide range of distances from the parent. Application of mathematical criteria analogous to those that have been used to help define the ‘gravitational clearing’ of planetary orbits around the Sun suggests that the absence of small satellites in this region around the planets may be a result (atleast in part) of gravitational clearing by the large moons present at these distances from the parent. The most significant exception to the observed diameter-distance distribution—Hyperion, on Saturn—is attributed to its 3:4 orbital resonance with Titan, while other obvious exceptions are the Trojan satellites of Saturn’s moons Tethys and Dione. The smallest satellite diameter that seems necessary for clearing of its ‘sphere of influence’ is around $400\,\mathrm{km}$ .
Multimessenger astronomy with a kHz-band gravitational-wave observatory
- Nikhil Sarin, Paul D. Lasky
-
- Published online by Cambridge University Press:
- 09 February 2022, e007
-
- Article
-
- You have access Access
- Open access
- HTML
- Export citation
- NASA ADS Abstract Service
-
Proposed next-generation networks of gravitational-wave observatories include dedicated kilohertz instruments that target neutron star science, such as the proposed Neutron Star Extreme Matter Observatory, NEMO. The original proposal for NEMO highlighted the need for it to exist in a network of gravitational-wave observatories to ensure detection confidence and sky localisation of sources. We show that NEMO-like observatories have significant utility on their own as coincident electromagnetic observations can provide the detection significance and sky localisation. We show that, with a single NEMO-like detector and expected electromagnetic observatories in the late 2020 s and early 2030 s such as the Vera C. Rubin observatory and SVOM, approximately 40% of all binary neutron star mergers detected with gravitational waves could be confidently identified as coincident multimessenger detections. We show that we expect $2^{+10}_{-1}{yr^{-1}}{}$ coincident observations of gravitational-wave mergers with gamma-ray burst prompt emission, $13^{+23}_{-10}{yr^{-1}}{}$ detections with kilonova observations, and $4^{+18}_{-3}{yr^{-1}}{}$ with broadband afterglows and kilonovae, where the uncertainties are 90% confidence intervals arising from uncertainty in current merger-rate estimates. Combined, this implies a coincident detection rate of $14^{+25}_{-11}{yr^{-1}}{}$ out to $300\,\mathrm{Mpc}$ . These numbers indicate significant science potential for a single kilohertz gravitational-wave detector operating without a global network of other gravitational-wave observatories.
A search for technosignatures toward the Galactic Centre at 150 MHz
- Chenoa D. Tremblay, Danny C. Price, Steven J. Tingay
-
- Published online by Cambridge University Press:
- 15 March 2022, e008
-
- Article
-
- You have access Access
- Open access
- HTML
- Export citation
- NASA ADS Abstract Service
-
This paper is the fourth in a series of low-frequency searches for technosignatures. Using the Murchison Widefield Array over two nights, we integrate 7 h of data toward the Galactic Centre (centred on the position of Sagittarius $\mathrm{A}^{*}$ ) with a total field-of-view of $200\,\mathrm{deg}^{2}$ . We present a targeted search toward 144 exoplanetary systems, at our best yet angular resolution (75 arcsec). This is the first technosignature search at a central frequency of 155 MHz toward the Galactic Centre (our previous central frequencies have been lower). A blind search toward in excess of 3 million stars toward the Galactic Centre and Galactic bulge is also completed, placing an equivalent isotropic power limit $<\!1.1\times10^{19}\,\mathrm{W}$ at the distance to the Galactic Centre. No plausible technosignatures are detected.
Understanding the impact of diffusion of CO in the astrochemical models
- Kinsuk Acharyya
-
- Published online by Cambridge University Press:
- 18 March 2022, e009
-
The mobility of lighter species on the surface of interstellar dust grains plays a crucial role in forming simple through complex molecules. Carbon monoxide is one of the most abundant molecules, its surface diffusion on the grain surface is essential to forming many molecules. Recent laboratory experiments found a diverse range of diffusion barriers for CO on the grain surface, their use can significantly impact the abundance of several molecules. The impact of different diffusion barriers of CO, in the astrochemical models, is studied to understand its effect on the abundance of solid CO and the species for which it is a reactant partner. A gas-grain network is used for three different physical conditions; cold core and warm-up models with slow and fast heating rates. Two different ratios (0.3 and 0.5) between diffusion and desorption barrier are utilised for all the species. For each physical condition and ratio, six different models are run by varying diffusion barriers of CO. Solid CO abundance for the models with the lowest diffusion barrier yields less than 0.1% of water ice for cold clouds and a maximum of 0.4% for slow and fast warm-up models. Also, solid $\textrm{CO}_2$ in dense clouds is significantly overproduced ( ${\sim}140\%$ of water). The abundance of H2CO and $\textrm{CH}_3\textrm{OH}$ showed an opposite trend, and HCOOH, $\textrm{CH}_3\textrm{CHO}$ , $\textrm{NH}_2\textrm{CO}$ , and $\textrm{CH}_3\textrm{COCH}_3$ are produced in lower quantities for models with low diffusion barriers for CO. Considerable variation in abundance is observed between models with the high and low diffusion barrier. Models with higher diffusion barriers provide a relatively better agreement with the observed abundances when compared with the models having lower diffusion barriers.
The First Large Absorption Survey in H i (FLASH): I. Science goals and survey design
- James R. Allison, E. M. Sadler, A. D. Amaral, T. An, S. J. Curran, J. Darling, A. C. Edge, S. L. Ellison, K. L. Emig, B. M. Gaensler, L. Garratt-Smithson, M. Glowacki, K. Grasha, B. S. Koribalski, C. del P. Lagos, P. Lah, E. K. Mahony, S. A. Mao, R. Morganti, V. A. Moss, M. Pettini, K. A. Pimbblet, C. Power, P. Salas, L. Staveley-Smith, M. T. Whiting, O. I. Wong, H. Yoon, Z. Zheng, M. A. Zwaan
-
- Published online by Cambridge University Press:
- 21 March 2022, e010
-
We describe the scientific goals and survey design of the First Large Absorption Survey in H i (FLASH), a wide field survey for 21-cm line absorption in neutral atomic hydrogen (H i) at intermediate cosmological redshifts. FLASH will be carried out with the Australian Square Kilometre Array Pathfinder (ASKAP) radio telescope and is planned to cover the sky south of $\delta \approx +40\,\deg$ at frequencies between 711.5 and 999.5 MHz. At redshifts between $z = 0.4$ and $1.0$ (look-back times of 4 – 8 Gyr), the H i content of the Universe has been poorly explored due to the difficulty of carrying out radio surveys for faint 21-cm line emission and, at ultra-violet wavelengths, space-borne searches for Damped Lyman- $\alpha$ absorption in quasar spectra. The ASKAP wide field of view and large spectral bandwidth, in combination with a radio-quiet site, will enable a search for absorption lines in the radio spectra of bright continuum sources over 80% of the sky. This survey is expected to detect at least several hundred intervening 21-cm absorbers and will produce an H i-absorption-selected catalogue of galaxies rich in cool, star-forming gas, some of which may be concealed from optical surveys. Likewise, at least several hundred associated 21-cm absorbers are expected to be detected within the host galaxies of radio sources at $0.4 < z < 1.0$ , providing valuable kinematical information for models of gas accretion and jet-driven feedback in radio-loud active galactic nuclei. FLASH will also detect OH 18-cm absorbers in diffuse molecular gas, megamaser OH emission, radio recombination lines, and stacked H i emission.
Intensity correlation functions of microwave maser sources
- A. V. Shepelev
-
- Published online by Cambridge University Press:
- 04 April 2022, e011
-
Astrophysical sources of microwave radiation with extremely high spectral brightness are interpreted as masers. But by itself, the information about high brightness of radiation does not make it possible to establish whether the radiation is thermal or maser. This can be determined only on the basis of the analysis of high-order correlation functions. A possible measurement procedure for the second-order autocorrelation function (the bunching parameter) for these sources is proposed.
Galaxy evolution through infrared and submillimetre spectroscopy: Measuring star formation and black hole accretion with JWST and ALMA
- Sabrina Mordini, Luigi Spinoglio, Juan Antonio Fernández-Ontiveros
-
- Published online by Cambridge University Press:
- 01 April 2022, e012
-
- Article
-
- You have access Access
- Open access
- HTML
- Export citation
- NASA ADS Abstract Service
-
Rest-frame mid- to far-infrared (IR) spectroscopy is a powerful tool to study how galaxies formed and evolved, because a major part of their evolution occurs in heavily dust enshrouded environments, especially at the so-called Cosmic Noon ( $1< z < 3$ ). Using the calibrations of IR lines and features, recently updated with Herschel and Spitzer spectroscopy, we predict their expected fluxes with the aim to measure the Star Formation (SF) and the Black Hole Accretion (BHA) rates in intermediate to high redshift galaxies. On the one hand, the recent launch of the James Webb Space Telescope (JWST) offers new mid-IR spectroscopic capabilities that will enable for the first time a detailed investigation of both the SF and the BHA obscured processes as a function of cosmic time. We make an assessment of the spectral lines and features that can be detected by JWST-MIRI in galaxies and active galactic nuclei up to redshift $z \sim 3$ . The fine structure lines of [MgIV]4.49 $\unicode{x03BC}\textrm{m}$ and [ArVI]4.53 $\unicode{x03BC}\textrm{m}$ can be used as BHA rate tracers for the $1 \lesssim z \lesssim 3$ range, and we propose the [NeVI]7.65 $\unicode{x03BC}\textrm{m}$ line as the best tracer for $z \lesssim 1.5$ . The [ArII]6.98 $\unicode{x03BC}\textrm{m}$ and [ArIII]8.99 $\unicode{x03BC}\textrm{m}$ lines can be used to measure the SF rate at $z \lesssim 3$ and $z \lesssim 2$ , respectively, while the stronger [NeII]12.8 $\unicode{x03BC}\textrm{m}$ line exits the JWST-MIRI spectral range above $z \gtrsim 1.2$ . At higher redshifts, the PAH features at 6.2 and 7.7 $\unicode{x03BC}\textrm{m}$ can be observed at $z \lesssim 3$ and $z \lesssim 2.7$ , respectively. On the other hand, rest-frame far-IR spectroscopic observations of high redshift galaxies ( $z \gtrsim 3$ ) have been collected with the Atacama Large Millimeter Array (ALMA) in the last few years. The observability of far-IR lines from ALMA depends on the observed frequency, due to the significant decrease of the atmospheric transmission at the highest frequencies ( $\gtrsim420\,\rm{Hz}$ ). The [CII]158 $\unicode{x03BC}\textrm{m}$ line is a reliable tracer of the SF rate and can in most cases ( $0.9 \lesssim z \lesssim 2$ and $2 \lesssim z \lesssim 9$ ) be observed. Additionally, we propose the use of the combination of [OIII]88 $\,\unicode{x03BC}$ m and [OI]145 $\,\unicode{x03BC}$ m lines as an alternative SF rate tracer, that can be detected above $z \gtrsim 3$ . Overall, we emphasize the importance of using multi-feature analysis to measure both BHA and SFR, since individual tracers can be strongly dependent on the local ISM conditions and vary from source to source. However, we conclude that the peak of the obscured SF and BHA activities at Cosmic Noon falls outside the wavelength coverage of facilities currently operating or under development. A new IR space telescope covering the full IR spectral range from about 10 to $300\,\unicode{x03BC}\textrm{m}$ and actively cooled to achieve high sensitivity, will be needed.
Searching for the synchrotron cosmic web again: A replication attempt
- Torrance Hodgson, Melanie Johnston-Hollitt, Benjamin McKinley, Natasha Hurley-Walker
-
- Published online by Cambridge University Press:
- 04 April 2022, e013
-
- Article
-
- You have access Access
- Open access
- HTML
- Export citation
- NASA ADS Abstract Service
-
We follow up on the surprising recent announcement by Vernstrom et al. (2021, MNRAS) of the detection of the synchrotron cosmic web. We attempt to reproduce their detection with new observations with the Phase II, extended configuration of the Murchison Widefield Array at 118.5 MHz. We reproduce their detection methodology by stacking pairs of nearby luminous red galaxies (LRGs)—used as tracers for clusters and galaxy groups—contained in our low-frequency radio observations. We show that our observations are significantly more sensitive than those used in Vernstrom et al. and that our angular sensitivity is sufficient. And yet, we make no statistically significant detection of excess radio emission along the bridge spanning the LRG pairs. This non-detection is true both for the original LRG pair catalogue as used in Vernstrom et al., as well as for other larger catalogues with modified selection criteria. Finally, we return to the original data sets used in Vernstrom et al., and find that whilst we clearly reproduce the excess X-ray emission from ROSAT, we are not able to reproduce any kind of broad and extended excess intercluster filamentary emission using the original 118.5 MHz MWA survey data. In the interests of understanding this result, as part of this paper we release images of the 14 fields used in this study, the final stacked images, as well as key components of our stacking and modelling code.
A geometric view of closure phases in interferometry
- Nithyanandan Thyagarajan, Chris L. Carilli
-
- Published online by Cambridge University Press:
- 07 April 2022, e014
-
- Article
-
- You have access Access
- Open access
- HTML
- Export citation
- NASA ADS Abstract Service
-
Closure phase is the phase of a closed-loop product of spatial coherences formed by a ${\ge}3$ -element interferometer array. Its invariance to phase corruption attributable to individual array elements acquired during the propagation and the measurement processes, subsequent calibration, and errors therein, makes it a valuable tool in interferometry applications that otherwise require high-accuracy phase calibration. However, its understanding has remained mainly mathematical and limited to the aperture plane (Fourier dual of the image plane). Here, we present a geometrical, image domain view of closure phase, which until now has been lacking. Using the principal triangle in a 3-element interference image formed by a triad of interferometer elements, we show that the properties of closure phase, particularly its invariance to multiplicative element-based corruption factors (even of a large magnitude) and to translation, are intricately related to the conserved properties of the triangle, namely, its shape, orientation, and size, which is referred herein as the ‘shape-orientation-size (SOS) conservation principle’. In the absence of a need for element-based amplitude calibration of the interferometer array (as is typical in optical interferometry), the principal triangle in any 3-element interference image formed from phase-uncalibrated spatial coherences is still a true and uncorrupted representation of the source object’s morphology, except for a possible shift. Based on this knowledge of the triangle SOS conservation principle, we present two geometric methods to measure the closure phase directly from a simple 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 the geometric understanding of closure phase in the image plane using observations with the Karl G. Jansky Very Large Array, and the Event Horizon Telescope. These results verify the SOS conservation principle across a wide range of radio interferometric conditions. This geometric insight can be potentially valuable to other interferometric applications, such as optical interferometry. We also generalise these geometric relationships to an N-element interferometer.
What is the SKA-Low sensitivity for your favourite radio source?
- M. Sokolowski, S. J. Tingay, D. B. Davidson, R. B. Wayth, D. Ung, J. Broderick, B. Juswardy, M. Kovaleva, G. Macario, G. Pupillo, A. Sutinjo
-
- Published online by Cambridge University Press:
- 15 April 2022, e015
-
- Article
-
- You have access Access
- Open access
- HTML
- Export citation
- NASA ADS Abstract Service
-
The Square Kilometre Array (SKA) will be the largest radio astronomy observatory ever built, providing unprecedented sensitivity over a very broad frequency band from 50 MHz to 15.3 GHz. The SKA’s low frequency component (SKA-Low), which will observe in the 50–350 MHz band, will be built at the Murchison Radio-astronomy Observatory (MRO) in Western Australia. It will consist of 512 stations each composed of 256 dual-polarised antennas, and the sensitivity of an individual station is pivotal to the performance of the entire SKA-Low telescope. The answer to the question in the title is, it depends. The sensitivity of a low frequency array, such as an SKA-Low station, depends strongly on the pointing direction of the digitally formed station beam and the local sidereal time (LST), and is different for the two orthogonal polarisations of the antennas. The accurate prediction of the SKA-Low sensitivity in an arbitrary direction in the sky is crucial for future observation planning. Here, we present a sensitivity calculator for the SKA-Low radio telescope, using a database of pre-computed sensitivity values for two realisations of an SKA-Low station architecture. One realisation uses the log-periodic antennas selected for SKA-Low. The second uses a known benchmark, in the form of the bowtie dipoles of the Murchison Widefield Array. Prototype stations of both types were deployed at the MRO in 2019, and since then have been collecting commissioning and verification data. These data were used to measure the sensitivity of the stations at several frequencies and over at least 24 h intervals, and were compared to the predictions described in this paper. The sensitivity values stored in the SQLite database were pre-computed for the X, Y, and Stokes I polarisations in 10 MHz frequency steps, $\scriptsize{1/2}$ hour LST intervals, and $5^\circ$ resolution in pointing directions. The database allows users to quickly and easily estimate the sensitivity of SKA-Low for arbitrary observing parameters (your favourite object) using interactive web-based or command line interfaces. The sensitivity can be calculated using publicly available web interface (http://sensitivity.skalow.link) or a command line python package (https://github.com/marcinsokolowski/station_beam), which can also be used to calculate the sensitivity for arbitrary pointing directions, frequencies, and times without interpolations.
HST WFC3/Grism observations of the candidate ultra-high-redshift radio galaxy GLEAM J0917–0012
- N. Seymour, G. Drouart, G. Noirot, J. W. Broderick, R. J. Turner, S. S. Shabala, D. K. Stern, S. Bellstedt, S. Driver, L. Davies, C. A. De Breuck, J. A. Afonso, J. D. R. Vernet, T. J. Galvin
-
- Published online by Cambridge University Press:
- 12 April 2022, e016
-
We present Hubble Space Telescope Wide Field Camera 3 photometric and grism observations of the candidate ultra-high-redshift ( $z>7$ ) radio galaxy, GLEAM J0917–0012. This radio source was selected due to the curvature in its 70–230 MHz, low-frequency Murchison Widefield Array radio spectrum and its faintness in K-band. Follow-up spectroscopic observations of this source with the Jansky Very Large Array and Atacama Large Millimetre Array were inconclusive as to its redshift. Our F105W and F0986M imaging observations detect the host of GLEAM J0917–0012 and a companion galaxy, $\sim$ one arcsec away. The G102 grism observations reveal a single weak line in each of the spectra of the host and the companion. To help identify these lines we utilised several photometric redshift techniques including template fitting to the grism spectra, fitting the ultraviolet (UV)-to-radio photometry with galaxy templates plus a synchrotron model, fitting of the UV-to-near-infrared photometry with EAZY, and fitting the radio data alone with RAiSERed. For the host of GLEAM J0917–0012 we find a line at $1.12\,\mu$ m and the UV-to-radio spectral energy distribution (SED) fitting favours solutions at $z\sim 2$ or $z\sim 8$ . While this fitting shows a weak preference for the lower redshift solution, the models from the higher redshift solution are more consistent with the strength of the spectral line. The redshift constraint by RAiSERed of $>6.5$ also supports the interpretation that this line could be Lyman $-\alpha$ at $z=8.21$ ; however EAZY favours the $z\sim 2$ solution. We discuss the implications of both solutions. For the companion galaxy we find a line at $0.98\,\mu$ m and the SED fitting favours solutions at $z<3$ implying that the line could be the [OII] $\lambda3727$ doublet at $z=1.63$ (although the EAZY solution is $z\sim 2.6\pm 0.5$ ). Further observations are still required to unambiguously determine the redshift of this intriguing candidate ultra-high-redshift radio galaxy.
Imaging the southern sky at 159 MHz using spherical harmonics with the engineering development array 2
- Michael A. Kriele, Randall B. Wayth, Mark J. Bentum, Budi Juswardy, Cathryn M. Trott
-
- Published online by Cambridge University Press:
- 22 April 2022, e017
-
- Article
-
- You have access Access
- Open access
- HTML
- Export citation
- NASA ADS Abstract Service
-
One of the major priorities of international radio astronomy is to study the early universe through the detection of the 21 cm HI line from the epoch of reionisation (EoR). Due to the weak nature of the 21 cm signal, an important part in the detection of the EoR is removing contaminating foregrounds from our observations as they are multiple orders of magnitude brighter. In order to achieve this, sky maps spanning a wide range of frequencies and angular scales are required for calibration and foreground subtraction. Complementing the existing low-frequency sky maps, we have constructed a Southern Sky map through spherical harmonic transit interferometry utilising the Engineering Development Array 2 (EDA2), a Square Kilometre Array (SKA) low-frequency array prototype system. We use the m-mode formalism to create an all-sky map at 159 MHz with an angular resolution of 3 degrees, with data from the EDA2 providing information over +60 degrees to –90 degrees in declination. We also introduce a new method for visualising and quantifying how the baseline distribution of an interferometer maps to the spherical harmonics and discuss how prior information can be used to constrain spherical harmonic components that the interferometer is not sensitive to.
System design and calibration of SITARA—a global 21 cm short spacing interferometer prototype
- Jishnu N. Thekkeppattu, Benjamin McKinley, Cathryn M. Trott, Jake Jones, Daniel C. X. Ung
-
- Published online by Cambridge University Press:
- 21 April 2022, e018
-
Global 21-cm experiments require exquisitely precise calibration of the measurement systems in order to separate the weak 21-cm signal from Galactic and extragalactic foregrounds as well as instrumental systematics. Hitherto, experiments aiming to make this measurement have concentrated on measuring this signal using the single element approach. However, an alternative approach based on interferometers with short baselines is expected to alleviate some of the difficulties associated with a single element approach such as precision modelling of the receiver noise spectrum. Short spacing Interferometer Telescope probing cosmic dAwn and epoch of ReionisAtion (SITARA) is a short spacing interferometer deployed at the Murchison Radio-astronomy Observatory (MRO). It is intended to be a prototype or a test-bed to gain a better understanding of interferometry at short baselines, and develop tools to perform observations and data calibration. In this paper, we provide a description of the SITARA system and its deployment at the MRO, and discuss strategies developed to calibrate SITARA. We touch upon certain systematics seen in SITARA data and their modelling. We find that SITARA has sensitivity to all sky signals as well as non-negligible noise coupling between the antennas. It is seen that the coupled receiver noise has a spectral shape that broadly matches the theoretical calculations reported in prior works. We also find that when appropriately modified antenna radiation patterns taking into account the effects of mutual coupling are used, the measured data are well modelled by the standard visibility equation.
The FAST Ultra-Deep Survey (FUDS): Observational strategy, calibration and data reduction
- Hongwei Xi, Bo Peng, Lister Staveley-Smith, Bi-Qing For, Bin Liu
-
- Published online by Cambridge University Press:
- 26 April 2022, e019
-
- Article
-
- You have access Access
- HTML
- Export citation
- NASA ADS Abstract Service
-
The FAST Ultra-Deep Survey (FUDS) is a blind survey that aims for the direct detection of H i in galaxies at redshifts $z<0.42$ . The survey uses the multibeam receiver on the Five-hundred-metre Aperture Spherical Telescope (FAST) to map six regions, each of size $0.72\ \textrm{deg}^2$ at high sensitivity ( ${\sim}50\,\mu \textrm{Jy}$ ) and high-frequency resolution (23 kHz). The survey will enable studies of the evolution of galaxies and their H i content with an eventual sample size of ${\sim}1\,000$ . We present the science goals, observing strategy, the effects of radio frequency interference at the FAST site, our mitigation strategies and the methods for calibration, data reduction and imaging as applied to initial data. The observations and reductions for the first field, FUDS0, are completed, with around 128 H i galaxies detected in a preliminary analysis. Example spectra are given in this paper, including a comparison with data from the overlapping GAL2577 field of Arecibo Ultra-Deep Survey.
MWA tied-array processing IV: A multi-pixel beamformer for pulsar surveys and ionospheric corrected localisation
- N. A. Swainston, N. D. R. Bhat, I. S. Morrison, S. J. McSweeney, S. M. Ord, S. E. Tremblay, M. Sokolowski
-
- Published online by Cambridge University Press:
- 10 May 2022, e020
-
- Article
-
- You have access Access
- Open access
- HTML
- Export citation
- NASA ADS Abstract Service
-
The Murchison Widefield Array (MWA) is a low-frequency aperture array capable of high-time and frequency resolution astronomy applications such as pulsar studies. The large field-of-view of the MWA (hundreds of square degrees) can also be exploited to attain fast survey speeds for all-sky pulsar search applications, but to maximise sensitivity requires forming thousands of tied-array beams from each voltage-capture observation. The necessity of using calibration solutions that are separated from the target observation both temporally and spatially makes pulsar observations vulnerable to uncorrected, frequency-dependent positional offsets due to the ionosphere. These offsets may be large enough to move the source away from the centre of the tied-array beam, incurring sensitivity drops of ${\sim}30{-}50\%$ in Phase II extended array configuration. We analyse these offsets in pulsar observations and develop a method for mitigating them, improving both the source position accuracy and the sensitivity. This analysis prompted the development of a multi-pixel beamforming functionality that can generate dozens of tied-array beams simultaneously, which runs a factor of ten times faster compared to the original single-pixel version. This enhancement makes it feasible to observe multiple pulsars within the vast field of view of the MWA and supports the ongoing large-scale pulsar survey efforts with the MWA. We explore the extent to which ionospheric offset correction will be necessary for the MWA Phase III and the low-frequency square kilometre array (SKA-low).