Next Generation Of Radio Telescopes Research Articles

Imaging swiFTly: streaming widefield Fourier Transforms for large-scale interferometry

We describe a scalable distributed imaging algorithm framework for next-generation radio telescopes, managing the Fourier transform from apertures to sky (or vice versa) with a focus on minimising memory load, data transfers, and computation. Our algorithm uses smooth window functions to isolate the influence between specific regions of spatial-frequency and image space. This allows the distribution of image data between nodes and the construction of segments of frequency space exactly when and where needed. The developed prototype distributes terabytes of image data across many nodes, while generating visibilities at throughput and accuracy competitive with existing software. Scaling is demonstrated to be better than cubic in problem complexity (for baseline length and field of view), reducing the risk involved in growing radio astronomy processing to large telescopes like the Square Kilometre Array.

Exploring Neutral Hydrogen in Radio MOlecular Hydrogen Emission Galaxies (MOHEGs) and Prospects with the SKA

Empirical studies of cold gas content are essential for comprehending the star formation activities and evolution in galaxies. However, it is not straightforward to understand these processes because they depend on various physical properties of the interstellar medium. Massive Faranoff–Riley I/II type radio galaxies rich in molecular hydrogen with lower star formation activities are known as radio molecular hydrogen emission galaxies (MOHEGs). We present a study of neutral hydrogen-gas-associated radio MOHEGs at redshifts <0.2 probed via the H i 21 cm absorption line. Neutral hydrogen is detected in 70% of these galaxies, which are located at a distance of 8–120 kpc from the neighboring galaxies. These galaxies show a scarcity of H i gas as compared to merging galaxies at similar redshifts. We found no strong correlation between N(H i), NH, and the galaxy properties, regardless of whether the H i is assumed to be cold or warm, indicating that the atomic gas probably plays no important role in star formation. The relation between the total hydrogen gas surface density and the star formation surface density deviates from the standard Kennicutt–Schmidt law. Our study highlights the importance of H i studies and offers insights into the role of atomic and molecular hydrogen gas in explaining the properties of these galaxies. In the upcoming H i 21 cm absorption surveys with next-generation radio telescopes such as the Square Kilometre Array and pathfinder instruments, it may be possible to provide better constraints for these correlations.

Emulation of the cosmic dawn 21-cm power spectrum and classification of excess radio models using an artificial neural network

ABSTRACT The cosmic 21-cm line of hydrogen is expected to be measured in detail by the next generation of radio telescopes. The enormous data set from future 21-cm surveys will revolutionize our understanding of early cosmic times. We present a machine learning approach based on an artificial neural network that uses emulation in order to uncover the astrophysics in the epoch of reionization and cosmic dawn. Using a seven-parameter astrophysical model that covers a very wide range of possible 21-cm signals, over the redshift range 6 to 30 and wavenumber range 0.05 to $1 \ \rm {Mpc}^{-1}$ we emulate the 21-cm power spectrum with a typical accuracy of $10 - 20~{{\ \rm per\ cent}}$. As a realistic example, we train an emulator using the power spectrum with an optimistic noise model of the square kilometre array (SKA). Fitting to mock SKA data results in a typical measurement accuracy of 2.8 per cent in the optical depth to the cosmic microwave background, 34 per cent in the star-formation efficiency of galactic haloes, and a factor of 9.6 in the X-ray efficiency of galactic haloes. Also, with our modelling we reconstruct the true 21-cm power spectrum from the mock SKA data with a typical accuracy of $15 - 30~{{\ \rm per\ cent}}$. In addition to standard astrophysical models, we consider two exotic possibilities of strong excess radio backgrounds at high redshifts. We use a neural network to identify the type of radio background present in the 21-cm power spectrum, with an accuracy of 87 per cent for mock SKA data.

Probing quantum gravity with elastic interactions of ultrahigh-energy neutrinos

The next generation of radio telescopes will be sensitive to low-scale quantum gravity by measuring ultrahigh-energy neutrinos. In this work, we demonstrate for the first time that neutrino-nucleon soft interactions induced by TeV-scale gravity would significantly increase the number of events detected by the IceCube-Gen2 radio array in the EeV regime. However, we show that these experiments cannot measure the total cross section using only the angular and energy information of the neutrino flux, unless assumptions on the underlying inelasticity distribution of neutral interactions are made.

Introduction to Faraday tomography and its future prospects

Abstract Faraday tomography is a new method of the study of cosmic magnetic fields enabled by broad-band low-frequency radio observations. Using Faraday tomography it is possible to obtain the Faraday dispersion function, which contains information on the line-of-sight distributions of magnetic fields, thermal electron density, and cosmic ray electron density by measuring the polarization spectrum from a source of synchrotron radiation over a wide band. Furthermore, by combining it with two-dimensional imaging, Faraday tomography allows us to explore the three-dimensional structure of polarization sources. Faraday tomography has been active over the last 20 years, since the broad-band observation has become technically feasible, and polarization sources such as interstellar space, supernova remnants, and galaxies have been investigated. However, the Faraday dispersion function is mathematically the Fourier transform of the polarization spectrum. And since the observable band is finite, it is impossible to obtain a complete Faraday dispersion function by performing a Fourier transform. For this purpose, various methods have been developed to accurately estimate the Faraday dispersion function from the observed polarization spectrum. In addition, the Faraday dispersion function does not directly reflect the distribution of magnetic field, thermal electron density, and cosmic ray electron density in the physical space, and its physical interpretation is not straightforward. Despite these two difficult problems, Faraday tomography is attracting much attention because it has great potential as a new method for studying cosmic magnetic fields and magnetized plasmas. In particular, the next-generation radio telescope SKA (Square Kilometre Array) is capable of polarization observation with unprecedented sensitivity and broad bands, and the application of Faraday tomography is expected to make dramatic progress in the field of cosmic magnetic fields. In this review, we explain the basics of Faraday tomography with simple and instructive examples. Representative algorithms to realize Faraday tomography are introduced, and some applications are shown.

面向SKA1时代的科学数据流及阵列模拟分析

After years of planning for the next generation of radio telescopes, the Square Kilometer Array (SKA), construction of SKA phase one (SKA1) started in July 2021. After the formal operation of SKA1, it is expected that 710 petabytes of scientifically processed data will be generated every year. The data will be stored at SKA regional centers around the world for further analysis by scholars. This paper quantitatively analyzes the models of the textcolorblack

Iterative time-domain method for resolving multiple gravitational wave sources in pulsar timing array data

The sensitivity of ongoing searches for gravitational wave (GW) sources in the ultra-low frequency regime ($10^{-9}$ Hz to $10^{-7}$ Hz) using Pulsar Timing Arrays (PTAs) will continue to increase in the future as more well-timed pulsars are added to the arrays. It is expected that next-generation radio telescopes, namely, the Five-hundred-meter Aperture Spherical radio Telescope (FAST) and the Square Kilometer Array (SKA), will grow the number of well-timed pulsars to $O(10^3)$. The higher sensitivity will result in greater distance reach for GW sources, uncovering multiple resolvable GW sources in addition to an unresolved population. Data analysis techniques are, therefore, required that can search for and resolve multiple signals present simultaneously in PTA data. The multisource resolution problem in PTA data analysis poses a unique set of challenges such as non-uniformly sampled data, a large number of so-called pulsar phase parameters that arise from the inaccurately measured distances to the pulsars, and poor separation of signals in the Fourier domain due to a small number of cycles in the observed waveforms.

Radio detections of IR-selected runaway stellar bow shocks

ABSTRACT Massive stars moving at supersonic peculiar velocities through the interstellar medium (ISM) can create bow shocks, arc-like structures at the interface between the stellar wind and the ISM. Many such bow shocks have been detected and catalogued at IR wavelengths, but detections in other wavebands remain rare. Strikingly, while electrons are expected to be accelerated in the bow shock and their non-thermal emission may include synchrotron emission at low frequencies, only two massive runaway stellar bow shocks have to date been detected in the radio band. Here, we examine a sample of fifty IR-detected bow shocks from the E-BOSS catalogues in recently released radio images from the Rapid ASKAP Continuum Survey. We identify three confident and three likely counterparts, as well as three inconclusive candidates requiring confirmation via follow-up observations. These detections significantly increase the number of known radio massive stellar bow shocks and highlight the advantage of dedicated searches with current and next-generation radio telescopes. We investigate the underlying radio emission mechanism for these radio sources, finding a mix of free–free-dominated and synchrotron-dominated systems. We also discuss the non-detected targets by putting constraints on their emission properties and investigating their detectability with future observations. Finally, we propose several future avenues of research to advance the study and understanding of bow shocks at radio frequencies.

Search for Continuous Gravitational-wave Signals in Pulsar Timing Residuals: A New Scalable Approach with Diffusive Nested Sampling

Abstract Detecting continuous nanohertz gravitational waves (GWs) generated by individual close binaries of supermassive black holes (CB-SMBHs) is one of the primary objectives of pulsar timing arrays (PTAs). The detection sensitivity is slated to increase significantly as the number of well-timed millisecond pulsars will increase by more than an order of magnitude with the advent of next-generation radio telescopes. Currently, the Bayesian analysis pipeline using parallel tempering Markov Chain Monte Carlo has been applied in multiple studies for CB-SMBH searches, but it may be challenged by the high dimensionality of the parameter space for future large-scale PTAs. One solution is to reduce the dimensionality by maximizing or marginalizing over uninformative parameters semianalytically, but it is not clear whether this approach can be extended to more complex signal models without making overly simplified assumptions. Recently, the method of diffusive nested (DNest) sampling has shown capability in coping with high dimensionality and multimodality effectively in Bayesian analysis. In this paper, we apply DNest to search for continuous GWs in simulated pulsar timing residuals and find that it performs well in terms of accuracy, robustness, and efficiency for a PTA including  ( 10 2 ) pulsars. DNest also allows a simultaneous search of multiple sources elegantly, which demonstrates its scalability and general applicability. Our results show that it is convenient and also highly beneficial to include DNest in current toolboxes of PTA analysis.

Using the sample variance of 21 cm maps as a tracer of the ionisation topology

Intensity mapping of the 21 cm signal of neutral hydrogen will yield exciting insights into the Epoch of Reionisation and the nature of the first galaxies. However, the large amount of data that will be generated by the next generation of radio telescopes, such as the Square Kilometre Array, as well as the numerous observational obstacles to overcome, require analysis techniques tuned to extract the reionisation history and morphology. In this context, we introduce a one-point statistic, which we refer to as the local variance,σloc, that describes the distribution of the mean differential 21 cm brightness temperatures measured in two-dimensional maps along the frequency direction of a light cone. The local variance takes advantage of what is usually considered an observational bias, the sample variance. We find the redshift-evolution of the local variance to not only probe the reionisation history of the observed patches of the sky, but also trace the ionisation morphology. This estimator provides a promising tool to constrain the midpoint of reionisation as well as gain insight into the ionising properties of early galaxies.

Multi-epoch searches for relativistic binary pulsars and fast transients in the Galactic Centre

ABSTRACT The high stellar density in the central parsecs around the Galactic Centre makes it a seemingly favourable environment for finding relativistic binary pulsars. These include pulsars orbiting other neutron stars, stellar-mass black holes, or the central supermassive black hole, Sagittarius A*. Here, we present multi-epoch pulsar searches of the Galactic Centre at four observing frequencies, $4.85,\, 8.35,\, 14.6,\, 18.95\, {\rm GHz}$, using the Effelsberg 100-m radio telescope. Observations were conducted 1 yr prior to the discovery of, and during monitoring observations of, the Galactic Centre magnetar PSR J1745−2900. Our data analysis features acceleration searches on progressively shorter time series to maintain sensitivity to relativistic binary pulsars. The multi-epoch observations increase the likelihood of discovering transient or nulling pulsars, or ensure orbital phases are observed at which acceleration search methods work optimally. In ${\sim}147\, {\rm h}$ of separate observations, no previously undiscovered pulsars have been detected. Through calibration observations, we conclude this might be due to insufficient instantaneous sensitivity, caused by the intense continuum emission from the Galactic Centre, its large distance, and, at higher frequencies, the aggregate effect of steep pulsar spectral indices and atmospheric contributions to the system temperature. Additionally, we find that for millisecond pulsars in wide circular orbits $({\lesssim}800\, {\rm d})$ around Sagittarius A*, linear acceleration effects cannot be fully corrected in deep observations $(9\, {\rm h})$ with existing software tools. Pulsar searches of the Galactic Centre with the next generation of radio telescopes – such as MeerKat, ngVLA, and SKA1-mid – will have improved chances of uncovering this elusive population.

Research on Pulsar Candidate Identification Method Based on Deep Residual Neural Network

As the next-generation radio astronomical telescopes continuously improve and develop, the pulsar survey will produce millions of pulsar candidates, which pose considerable challenges for pulsar identification and classification. The rapidly evolving artificial intelligence (AI) techniques are being used for pulsar identification and discovery of new pulsars. Using the pulsar data set observed with the Parkes telescope, namely the High Time Resolution Universe Survey (HTRUS), a 14-layer deep residual network has been designed (called the Residual Network, ResNet) for pulsar candidate classifications. In the HTRUS sample data, the number of non-pulsar candidates (i.e., negative samples) is much larger than that of pulsar candidates (i.e., positive samples). The imbalance between the positive and negative samples is prone to result in model misjudgement. By using the over-sampling technique to enhance the data of positive samples in the training set and adjusting the ratio of positive and negative samples, we have solved this imbalance problem. In the training process, the hyperparameters are adjusted by means of 5-fold cross validation to build the model. The test results indicate that the model can achieve a high precision (98%) and recall (100%), the F1-score can reach 99%, and that the implementation of each sample test needs only 7 ms, it has provided a feasible approach for the future big-data analysis of pulsars.

Modelling neutron star–black hole binaries: future pulsar surveys and gravitational wave detectors

ABSTRACT Binaries comprised of a neutron star (NS) and a black hole (BH) have so far eluded observations as pulsars and with gravitational waves (GWs). We model the formation and evolution of these NS+BH binaries – including pulsar evolution – using the binary population synthesis code COMPAS. We predict the presence of a total of 50–2000 binaries containing a pulsar and a BH (PSR+BHs) in the Galactic field. We find the population observable by the next generation of radio telescopes, represented by the SKA and MeerKAT, current (LIGO/Virgo) and future (LISA) GW detectors. We conclude that the SKA will observe 1–80 PSR+BHs, with 0–4 binaries containing millisecond pulsars. MeerKAT is expected to observe 0–40 PSR+BH systems. Future radio detections of NS+BHs will constrain uncertain binary evolution processes such as BH natal kicks. We show that systems in which the NS formed first (NSBH) can be distinguished from those where the BH formed first (BHNS) by their pulsar and binary properties. We find 40 per cent of the LIGO/Virgo observed NS+BHs from a Milky Way like field population will have a chirp mass ≥3.0 M⊙. We estimate the spin distributions of NS+BHs with two models for the spins of BHs. The remnants of BHNS mergers will have a spin of ∼0.4, whilst NSBH merger remnants can have a spin of ∼0.6 or ∼0.9 depending on the model for BH spins. We estimate that approximately 25–1400 PSR+BHs will be radio alive whilst emitting GWs in the LISA frequency band, raising the possibility of joint observation by the SKA and LISA.

Recrystallization Behavior of a Pure Cu Connection Interface with Ultrasonic Welding

Three-dimensional metal waveguide components are key components in the next generation of radio telescopes. Ultrasonic additive manufacturing technology combining ultrasonic welding and micro electrical discharge machining (micro-EDM) provides a new method for the overall manufacturing of waveguide elements, and the effective welding of Electrolytic Tough Pitch copper (Cu-ETP) sheets is the key process of this method. This study demonstrates that the orthogonal test optimization method is used to conduct ultrasonic welding tests on Cu-ETP. Specifically, electron backscattered diffraction (EBSD) technology is used to analyze the crystal grains, grain boundary types and texture changes during interface recrystallization. In addition, the finite element software ABAQUS 6.13 is employed to calculate the temperature field in order to determine the possibility of recrystallization of the welding interface. The results showed that the average grain size of the welding interface decreased from 20 to 1~2 μm. The Cu-ETP matrix is mainly composed of coarse grains with high-angle grain boundaries (HAGBs), while a large number of low-angle grain boundaries (LAGBs), subcrystals and fine equiaxed grains appear in the welded joint. At the same time, discontinuous dynamic recrystallization (DDRX) occurs in the less strained area, and continuous dynamic recrystallization (CDRX) is predominant in the greater strain area. The temperature field calculation shows that the peak temperature of the welding interface exceeds the recrystallization temperature of Cu-ETP from 379.05 to 433.2 °C.

Lux Ex Tenebris: The Imprint of Annihilating Dark Matter on the Intergalactic Medium during Cosmic Dawn

Abstract Upcoming measurements of the highly redshifted 21 cm line with next-generation radio telescopes such as the Hydrogen Epoch of Reionization Array (HERA) and Square Kilometer Array will provide the intriguing opportunity to probe dark matter (DM) physics during the Epoch of Reionization, Cosmic Dawn, and the Dark Ages. With HERA already under construction, there is a pressing need to thoroughly understand the impact of DM physics on the intergalactic medium (IGM) during these epochs. We present first results of a hydrodynamic simulation suite with 2 × 5123 particles in a (100 h −1 Mpc)3 box with DM annihilation and baryonic cooling physics. We focus on redshift z ∼ 11, just before reionization starts in our simulations, and discuss the imprint of DM annihilation on the IGM and on structure formation. We find that whereas structure formation is not affected by thermal weakly interacting massive particles heavier than m χ ≳ 100 MeV, heating from (GeV) DM particles may leave a significant imprint on the IGM that alters the 21 cm signal. Cold gas in low-density regions is particularly susceptible to the effects of DM heating. We note, however, that delayed energy deposition is not currently accounted for in our simulations.

Phase synchronization system receiver module for the mid-frequency Square Kilometre Array

Next-generation radio telescopes, such as the Square Kilometre Array (SKA) and Next Generation Very Large Array (ngVLA), require precise microwave frequency reference signals to be transmitted over fiber links to each dish to coherently sample astronomical signals. Such telescopes employ phase stabilization systems to suppress the phase noise imparted on the reference signals by environmental perturbations on the links; however, the stabilization systems are bandwidth limited by the round-trip time of light traveling on the fiber links. A phase-locked receiver module (RM) is employed on each dish to suppress residual phase noise outside the round-trip bandwidth. The SKA RM must deliver a 3.96-GHz output signal with 4 MHz of tuning range and <100 fs of timing jitter. We present an RM architecture to meet both requirements. Analytical modeling of the RM predicts 30 fs of output jitter when the reference signal is integrated between 1 Hz and 2.8 GHz. The proposed RM was conceived with best practice electromagnetic compatibility in mind and to meet size, weight, and power requirements for the SKA dish indexer. As the ngVLA reference design also incorporates a round-trip phase stabilization system, this RM may be applicable to future ngVLA design.

HIR4: cosmology from a simulated neutral hydrogen full sky using Horizon Run 4

ABSTRACT The distribution of cosmological neutral hydrogen will provide a new window into the large-scale structure of the Universe with the next generation of radio telescopes and surveys. The observation of this material, through 21 cm line emission, will be confused by foreground emission in the same frequencies. Even after these foregrounds are removed, the reconstructed map may not exactly match the original cosmological signal, which will introduce systematic errors and offset into the measured correlations. In this paper, we simulate future surveys of neutral hydrogen using the Horizon Run 4 (HR4) cosmological N-body simulation. We generate H i intensity maps from the HR4 halo catalogue, and combine with foreground radio emission maps from the Global Sky Model, to create accurate simulations over the entire sky. We simulate the H i sky for the frequency range 700–800 MHz, matching the sensitivity of the Tianlai pathfinder. We test the accuracy of the fastICA, PCA, and log-polynomial fitting foreground removal methods to recover the input cosmological angular power spectrum and measure the parameters. We show the effect of survey noise levels and beam sizes on the recovered the cosmological constraints. We find that while the reconstruction removes power from the cosmological 21 cm distribution on large scales, we can correct for this and recover the input parameters in the noise-free case. However, the effect of noise and beam size of the Tianlai pathfinder prevents accurate recovery of the cosmological parameters when using only intensity mapping information.

Unique Multimessenger Signal of QCD Axion Dark Matter.

We propose a multimessenger probe of QCD axion dark matter based on observations of black hole-neutron star binary inspirals. It is suggested that a dense dark matter spike may grow around intermediate mass black holes (10^{3}-10^{5} M_{⊙}). The presence of such a spike produces two unique effects: a distinct phase shift in the gravitational wave strain during the inspiral and an enhancement of the radio emission due to the resonant axion-photon conversion occurring in the neutron star magnetosphere throughout the inspiral and merger. Remarkably, the observation of the gravitational wave signal can be used to infer the dark matter density and, consequently, to predict the radio emission. We study the projected reach of the LISA interferometer and next-generation radio telescopes such as the Square Kilometre Array. Given a sufficiently nearby system, such observations will potentially allow for the detection of QCD axion dark matter in the mass range 10^{-7} eV to 10^{-5} eV.

Dark Matter searches through cross-correlations of gamma rays with neutral-hydrogen intensity mapping

In this work we derive the first theoretical prediction of the cross-correlation signal between the unresolved gamma-ray background and the 21-cm line originated by the spin-flip transition of neutral hydrogen atoms, by taking as benchmark experiment the space telescope Fermi-LAT for gamma rays and the next-generation radio telescope Square Kilometer Array (SKA) as well as its precursor MeerKAT for the 21-cm emission. The attainable bounds in the dark matter (DM) parameter space are envisioned to be competitive already with the combination Fermi-LAT × MeerKAT, but SKA will allow to go deeper and probe a thermal DM particle up to masses of 130 GeV. A future gamma-ray detector with better angular resolution and larger exposure together with an ungraded radio telescope will have the potentiality to probe a DM candidate with thermal annihilation cross-section and masses up to the TeV scale.

Detecting axion-like particles at radio frequencies

Dark matter in the form of axion-like particles can decay into monochromatic photons. We investigate the prospects for detection of these signals with radio telescopes. We show that the presence of an ambient radiation field induces a stimulated enhancement of the decay rate. Depending on the environment, this effect can boost the photon signal by several orders of magnitude. We analyze different astrophysical targets, namely dwarf spheroidal galaxies, the Galactic Center and galaxy clusters. For axion-photon couplings allowed by astrophysical and laboratory constraints, the signal is within the reach of next-generation radio telescopes such as the Square Kilometer Array.