Search

CIZA J2242.8+5301, or the Sausage cluster, is well studied over a range of frequencies. Since its first discovery, a lot of interesting features and unique characteristics have been uncovered. In this work, we report some more new morphological features using the uGMRT band-3 and band-4 data. In the north relic, we observe variation in spectral index profiles across the relic width from the east to west, which may indicate a decrease in downstream cooling rate in that direction. We re-confirm the presence of an additional ~ 930 kpc relic in the north. We classify the filamentary source in the downstream region to be a narrow angle tail (NAT) radio galaxy. The bright arc in the east relic shows sub-structure in the spectral index profile, which may indicate the presence of finer filaments. We further report the presence of a double-strand structure in the east relic similar to the 'Toothbrush' relic. We categorize the bright 'L' shaped structure in the southern relic to be a NAT radio galaxy, as well as trace the actual ~ 1.1 Mpc relic component. We re-confirm the existence of the faint southern extent, measuring the relic length to be ~ 1.8 Mpc. Furthermore, we suggest the southern relic to be a union of individual component relics rather than a single giant filamentary relic. Lastly, based on the morphological symmetry between northern and southern relics, we suggest a schematic shock structure associated with the merger event in an attempt to explain their formation scenario., Comment: 23 pages, 18 figures, 3 tables, submitted to ApJ

Context: Interferometric imaging is algorithmically and computationally challenging as there is no unique inversion from the measurement data back to the sky maps, and the datasets can be very large. Many imaging methods already exist, but most of them focus either on the accuracy or the computational aspect. Aims: This paper aims to reduce the computational complexity of the Bayesian imaging algorithm resolve, enabling the application of Bayesian imaging for larger datasets. Methods: By combining computational shortcuts of the CLEAN algorithm with the Bayesian imaging algorithm resolve we developed an accurate and fast imaging algorithm which we name fast-resolve. Results: We validate the accuracy of the presented fast-resolve algorithm by comparing it with results from resolve on VLA Cygnus A data. Furthermore, we demonstrate the computational advantages of fast-resolve on a large MeerKAT ESO 137-006 dataset which is computationally out of reach for resolve. Conclusions: The presented algorithm is significantly faster than previous Bayesian imaging algorithms, broadening the applicability of Bayesian interferometric imaging. Specifically for the single channel VLA Cygnus A datasets fast-resolve is about $144$ times faster than resolve. For the MeerKAT dataset with multiple channels the computational speedup of fast-resolve is even larger., Comment: 24 pages, 14 figures

We present the results of a polarimetric study from our new high-sensitivity L-band (0.8--1.7 GHz) observation of Pictor A with the MeerKAT radio telescope. We confirm the presence of the radio jet extending from the nucleus to the western hotspot of this source. Additionally, we show the radio emission expected to be coincident with previously observed X-ray emission in the radio lobes, confirming that the emission mechanism is of inverse Compton origin, as suggested by a previous study. Our spectropolarimetric analysis using the RM-Synthesis technique reveals a relatively uniform mean RM distribution across the lobes of Pictor A, with most lines-of-sight exhibiting single-peaked Faraday spectra. However, a number of the lines-of-sight exhibit single peaked spectra with a wide base or multiple peaks, suggesting the presence of multiple Faraday components or a Faraday thick structure along Pictor A's lines-of-sight. We also confirm the asymmetry in RM variability and depolarization between the two lobes of this source which were reported in a previous study., Comment: 15 pages,15 figures. This is a pre-copyedited, author-produced PDF of an article accepted for publication in MNRAS following peer review. The version of record [Andati L.A.L. et al., 2024] is available online at: https://doi.org/10.1093/mnras/stae598

The past decade has seen significant advances in wide-field cm-wave very long baseline interferometry (VLBI), which is timely given the wide-area, synoptic survey-driven strategy of major facilities across the electromagnetic spectrum. While wide-field VLBI poses significant post-processing challenges that can severely curtail its potential scientific yield, many developments in the km-scale connected-element interferometer sphere are directly applicable to addressing these. Here we present the design, processing, data products, and source counts from a deep (11 $\mu$Jy beam$^{-1}$), quasi-uniform sensitivity, contiguous wide-field (160 arcmin$^2$) 1.6 GHz VLBI survey of the CANDELS GOODS-North field. This is one of the best-studied extragalactic fields at milli-arcsecond resolution and, therefore, is well-suited as a comparative study for our Tera-pixel VLBI image. The derived VLBI source counts show consistency with those measured in the COSMOS field, which broadly traces the AGN population detected in arcsecond-scale radio surveys. However, there is a distinctive flattening in the $ S_{\rm 1.4GHz}\sim$100-500 $\mu$Jy flux density range, which suggests a transition in the population of compact faint radio sources, qualitatively consistent with the excess source counts at 15 GHz that is argued to be an unmodelled population of radio cores. This survey approach will assist in deriving robust VLBI source counts and broadening the discovery space for future wide-field VLBI surveys, including VLBI with the Square Kilometre Array, which will include new large field-of-view antennas on the African continent at $\gtrsim$1000~km baselines. In addition, it may be useful in the design of both monitoring and/or rapidly triggered VLBI transient programmes., Comment: 15 pages, 8 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

Radio phoenices are complex and filamentary diffuse radio sources found in both merging and relaxed clusters. The formation of these sources was proposed to be due to adiabatic compression of old Active Galactic Nucleus (AGN) plasma in shock waves. Most of the previous spectral studies of these sources were limited to integrated spectral indices, which were found to be very steep as well as showing a curved spectrum. Here, we have performed a multi-frequency investigation of the radio phoenix in the Abell 85 cluster. Owing to the sensitive high-resolution observations, we found some of the finer filamentary structures that were previously undetected. We produced resolved spectral index maps of the radio phoenix between 323, 700, and 1280 MHz. The orientation of the filaments, as well as the gradient across the spectral index maps, suggest the possible direction of the shock motion from northeast to southwest. The integrated spectrum of the radio phoenix was found to be very steep and curved toward high-frequencies. Furthermore, the spectral index of the filaments was found to be less steep compared to the non-filamentary regions, implying greater energy injection in the filaments. The observed features in the radio phoenix in the Abell 85 cluster seem to be in support of the adiabatic shock compression mechanism., Comment: 20 pages, 12 figures, 4 tables, accepted for ApJ publication

Galaxy clusters are located at the nodes of cosmic filaments and therefore host a lot of hydrodynamical activity. However, cool core clusters are considered to be relatively relaxed systems without much merging activity. The Abell 85 cluster is a unique example where the cluster hosts both a cool core and multiple ongoing merging processes. In this work, we used 700 MHz uGMRT as well as MeerKAT L-band observations, carried out as part of the MGCLS, of the Abell 85. We reconfirm the presence of a minihalo in the cluster centre at 700MHz that was recently discovered in MGCLS. Furthermore, we discovered a radio bridge connecting the central minihalo and the peripheral radio phoenix. The mean surface brightness, size and flux density of the bridge at 700 MHz is found to be $\sim 0.14\ \mu$Jy/arcsec$^2$, $\sim 220$ kpc and $\sim 4.88$ mJy, respectively, with a spectral index of $\alpha_{700}^{1.28} = -0.92$. Although the origin of the seed relativistic electrons is still unknown, turbulent re-acceleration caused by both the spiralling sloshing gas in the intracluster medium (ICM) and the post-shock turbulence from the outgoing merging shock associated with the phoenix formation may be responsible for the bridge., Comment: 7 pages, 6 Figures, 2 Tables. Accepted for publication in MNRAS-L

The transformation and evolution of a galaxy is strongly influenced by interactions with its environment. Neutral hydrogen (\HI) is an excellent way to trace these interactions. Here, we present \HI\ observations of the spiral galaxy NGC~895, which was previously thought to be isolated. High-sensitivity \HI\ observations from the MeerKAT large survey project MIGHTEE reveal possible interaction features, such as extended spiral arms, and the two newly discovered \HI\ companions, that drive us to change the narrative that it is an isolated galaxy. We combine these observations with deep optical images from the Hyper Suprime Camera to show an absence of tidal debris between NGC 895 and its companions. We do find an excess of light in the outer parts of the companion galaxy MGTH$\_$J022138.1-052631 which could be an indication of external perturbation and thus possible sign of interactions. Our analysis shows that NGC~895 is an actively star-forming galaxy with a SFR of $\mathrm{1.75 \pm 0.09 [M_{\odot}/yr]}$, a value typical for high stellar mass galaxies on the star forming main sequence. It is reasonable to state that different mechanisms may have contributed to the observed features in NGC~895 and this emphasizes the need to revisit the target with more detailed observations. Our work shows the high potential and synergy of using state-of-the-art data in both \HI\ and optical to reveal a more complete picture of galaxy environments., Comment: 14 pages, 10 figures. Accepted for publication in MNRAS

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., Comment: 16 pages, 14 figures, accepted for publication by MNRAS

We introduce the first AI-based framework for deep, super-resolution, wide-field radio-interferometric imaging, and demonstrate it on observations of the ESO~137-006 radio galaxy. The algorithmic framework to solve the inverse problem for image reconstruction builds on a recent ``plug-and-play'' scheme whereby a denoising operator is injected as an image regulariser in an optimisation algorithm, which alternates until convergence between denoising steps and gradient-descent data-fidelity steps. We investigate handcrafted and learned variants of high-resolution high-dynamic range denoisers. We propose a parallel algorithm implementation relying on automated decompositions of the image into facets and the measurement operator into sparse low-dimensional blocks, enabling scalability to large data and image dimensions. We validate our framework for image formation at a wide field of view containing ESO~137-006, from 19 gigabytes of MeerKAT data at 1053 and 1399 MHz. The recovered maps exhibit significantly more resolution and dynamic range than CLEAN, revealing collimated synchrotron threads close to the galactic core., Comment: accepted for publication in ApJL

We present Tricolour, a package for Radio Frequency Interference mitigation of wideband finely channelized MeerKAT correlation data. The MeerKAT passband is heavily affected by interference from satellite, mobile, aircraft and terrestrial transponders. Coupled with typical data rates in excess of 100 GiB/hr at 208kHz channelization resolution, mitigation poses a significant processing challenge. Our flagger is highly configurable, parallel and optimized, employing Dask and Numba technologies to implement the widely used SumThreshold and MAD interference detection algorithms. We find that typical 208kHz channelized datasets can be processed at rates in excess of 400 GiB/hr for a typical L-band flagging strategy on a modern dual-socket Intel Xeon server., Comment: This submission to the arXiv is a preprint of a submission to the conference proceedings of the Astronomical Data Analysis Software and Systems XXX to be published in the Astronomical Society of the Pacific conference series (J.E. Ruiz, F. Pierfederici, and P. Teuben, eds.) ). 4 pages, 3 figures (9 eps images total)

Observations of the neutral atomic hydrogen (HI) in the nuclear starburst galaxy NGC 4945 with MeerKAT are presented. We find a large amount of halo gas, previously missed by HI observations, accounting for 6.8% of the total HI mass. This is most likely gas blown into the halo by star formation. Our maps go down to a $3\sigma$ column density level of $5\times10^{18} cm^{-2}$ . We model the HI distribution using tilted-ring fitting techniques and find a warp on the galaxy's approaching and receding sides. The HI in the northern side of the galaxy appears to be suppressed. This may be the result of ionisation by the starburst activity in the galaxy, as suggested by a previous study. The origin of the warp is unclear but could be due to past interactions or ram pressure stripping. Broad, asymmetric HI absorption lines extending beyond the HI emission velocity channels are present towards the nuclear region of NGC 4945. Such broad lines suggest the existence of a nuclear ring moving at a high circular velocity. This is supported by the clear rotation patterns in the HI absorption velocity field. The asymmetry of the absorption spectra can be caused by outflows or inflows of gas in the nuclear region of NGC 4945. The continuum map shows small extensions on both sides of the galaxy's major axis that might be signs of outflows resulting from the starburst activity., Comment: 20 pages, 15 figures, 3 tables, accepted for publication in MNRAS

We report on the detection of a large, extended HI cloud complex in the GAMA G23 field, located at a redshift of $z\,\sim\,0.03$, observed as part of the MeerHOGS campaign (a pilot survey to explore the mosaicing capabilities of MeerKAT). The cloud complex, with a total mass of $10^{10.0}\,M_\odot$, lies in proximity to a large galaxy group with $M_\mathrm{dyn}\sim10^{13.5}\,M_\odot$. We identify seven HI peak concentrations, interconnected as a tenuous 'chain' structure, extending $\sim 400\,\mathrm{kpc}$ from east-to-west, with the largest (central) concentration containing $10{^{9.7}}\,M_\odot$ in HI gas distributed across $50\,\mathrm{kpc}$. The main source is not detected in ultra-violet, optical or infrared imaging. The implied gas mass-to-light ($M_\mathrm{HI}$/$L_\mathrm{r}$) is extreme ($>$1000) even in comparison to other 'dark clouds'. The complex has very little kinematic structure ($110\,\mathrm{km}\,\mathrm{s}^{-1}$), making it difficult to identify cloud rotation. Assuming pressure support, the total mass of the central concentration is $>10^{10.2}\,M_\odot$, while a lower limit to the dynamical mass in the case of full rotational support is $10^{10.4}\,M_\odot$. If the central concentration is a stable structure, it has to contain some amount of unseen matter, but potentially less than is observed for a typical galaxy. It is, however, not clear whether the structure has any gravitationally stable concentrations. We report a faint UV--optical--infrared source in proximity to one of the smaller concentrations in the gas complex, leading to a possible stellar association. The system nature and origins is enigmatic, potentially being the result of an interaction with or within the galaxy group it appears to be associated with., Comment: 25 pages, 10 figures, accepted for publication in ApJ, minor corrections

The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a radio interferometer array currently in development, with an initial 256-element array to be deployed at the South African Radio Astronomy Observatory (SARAO) Square Kilometer Array (SKA) site in South Africa. Each of the 6m, $f/0.23$ dishes will be instrumented with dual-polarisation feeds operating over a frequency range of 400-800 MHz. Through intensity mapping of the 21 cm emission line of neutral hydrogen, HIRAX will provide a cosmological survey of the distribution of large-scale structure over the redshift range of $0.775 < z < 2.55$ over $\sim$15,000 square degrees of the southern sky. The statistical power of such a survey is sufficient to produce $\sim$7 percent constraints on the dark energy equation of state parameter when combined with measurements from the Planck satellite. Additionally, HIRAX will provide a highly competitive platform for radio transient and HI absorber science while enabling a multitude of cross-correlation studies. In this paper, we describe the science goals of the experiment, overview of the design and status of the sub-components of the telescope system, and describe the expected performance of the initial 256-element array as well as the planned future expansion to the final, 1024-element array., Comment: 24 pages, 6 figures. Updated to match published version

Solution intervals are often used to improve the signal-to-noise ratio during radio interferometric gain calibration. This work investigates how factors such as the noise level, intrinsic gain variability, degree of model incompleteness, and the presence of radio frequency interference impact the selection of solution intervals for calibration. We perform different interferometric simulations to demonstrate how these factors, in combination with the choice of solution intervals, affect calibration and imaging outputs and discuss practical guidelines for choosing optimal solution intervals. Furthermore, we present an algorithm capable of automatically selecting suitable solution intervals during calibration. By applying the algorithm to both simulated and real data, we show that it can successfully choose solution intervals that strike a good balance between capturing intrinsic gain variability and not fitting noise as long as the data are not too inhomogeneously flagged. Furthermore, we elaborate on several practical aspects that emphasize the need to develop regularised calibration algorithms that do not require solution intervals., Comment: 18 pages, 14 figures

Radio interferometric (RI) data are noisy under-sampled spatial Fourier components of the unknown radio sky affected by direction-dependent antenna gains. Failure to model these antenna gains accurately results in a radio sky estimate with limited fidelity and resolution. The RI inverse problem has been recently addressed via a joint calibration and imaging approach which consists in solving a non-convex minimisation task, involving suitable priors for the DDEs, namely temporal and spatial smoothness, and sparsity for the unknown radio map via an $\ell_1$-norm prior, in the context of realistic RI simulations. Building on these developments, we propose to promote sparsity of the radio map via a log-sum prior, enforcing sparsity more strongly than the $\ell_1$-norm. The resulting minimisation task is addressed via a sequence of non-convex minimisation tasks composed of re-weighted $\ell_1$ image priors, which are solved approximately. We demonstrate the efficiency of the approach on RI observations of the celebrated radio galaxy Cygnus~A obtained with the Karl G. Jansky Very Large Array at X, C, and S bands. More precisely, we showcase that the approach enhances data fidelity significantly while achieving high resolution high dynamic range radio maps, confirming the suitability of the priors considered for the unknown DDEs and radio image. As a clear qualitative indication of the high fidelity achieved by the data and the proposed approach, we report the detection of three background sources in the vicinity of Cyg~A, at S band., Comment: 22 pages,13 figures. Submitted to MNRAS

The celebrated CLEAN algorithm has been the cornerstone of deconvolution algorithms in radio interferometry almost since its conception in the 1970s. For all its faults, CLEAN is remarkably fast, robust to calibration artefacts and in its ability to model point sources. We demonstrate how the same assumptions that afford CLEAN its speed can be used to accelerate more sophisticated deconvolution algorithms., Comment: ADASS 2020 proceedings to be published in ASP conference series

ESO 149-G003 is a close-by, isolated dwarf irregular galaxy. Previous observations with the ATCA indicated the presence of anomalous neutral hydrogen (HI) deviating from the kinematics of a regularly rotating disc. We conducted follow-up observations with the MeerKAT radio telescope during the 16-dish Early Science programme as well as with the MeerLICHT optical telescope. Our more sensitive radio observations confirm the presence of anomalous gas in ESO 149-G003, and further confirm the formerly tentative detection of an extraplanar HI component in the galaxy. Employing a simple tilted-ring model, in which the kinematics is determined with only four parameters but including morphological asymmetries, we reproduce the galaxy's morphology, which shows a high degree of asymmetry. By comparing our model with the observed HI, we find that in our model we cannot account for a significant (but not dominant) fraction of the gas. From the differences between our model and the observed data cube we estimate that at least 7%-8% of the HI in the galaxy exhibits anomalous kinematics, while we estimate a minimum mass fraction of less than 1% for the morphologically confirmed extraplanar component. We investigate a number of global scaling relations and find that, besides being gas-dominated with a neutral gas-to-stellar mass ratio of 1.7, the galaxy does not show any obvious global peculiarities. Given its isolation, as confirmed by optical observations, we conclude that the galaxy is likely currently acquiring neutral gas. It is either re-accreting gas expelled from the galaxy or accreting pristine intergalactic material., Comment: 16 pages, 14 figures, MNRAS, Accepted 2020 December 2, in original form 2019 September 18

We present observations and models of the kinematics and the distribution of the neutral hydrogen (HI) in the isolated dwarf irregular galaxy, Wolf-Lundmark-Melotte (WLM). We observed WLM with the Green Bank Telescope (GBT) and as part of the MeerKAT Early Science Programme, where 16 dishes were available. The HI disc of WLM extends out to a major axis diameter of 30 arcmin (8.5 kpc), and a minor axis diameter of 20 arcmin (5.6 kpc) as measured by the GBT. We use the MeerKAT data to model WLM using the TiRiFiC software suite, allowing us to fit different tilted-ring models and select the one that best matches the observation. Our final best-fitting model is a flat disc with a vertical thickness, a constant inclination and dispersion, and a radially-varying surface brightness with harmonic distortions. To simulate bar-like motions, we include second-order harmonic distortions in velocity in the tangential and the vertical directions. We present a model with only circular motions included and a model with non-circular motions. The latter describes the data better. Overall, the models reproduce the global distribution and the kinematics of the gas, except for some faint emission at the 2-sigma level. We model the mass distribution of WLM with a pseudo-isothermal (ISO) and a Navarro-Frenk-White (NFW) dark matter halo models. The NFW and the ISO models fit the derived rotation curves within the formal errors, but with the ISO model giving better reduced chi-square values. The mass distribution in WLM is dominated by dark matter at all radii., Comment: Accepted for publication in MNRAS, 25 pages, 21 figures, 5 tables

MeerKATHI is the current development name for a radio-interferometric data reduction pipeline, assembled by an international collaboration. We create a publicly available end-to-end continuum- and line imaging pipeline for MeerKAT and other radio telescopes. We implement advanced techniques that are suitable for producing high-dynamic-range continuum images and spectroscopic data cubes. Using containerization, our pipeline is platform-independent. Furthermore, we are applying a standardized approach for using a number of different of advanced software suites, partly developed within our group. We aim to use distributed computing approaches throughout our pipeline to enable the user to reduce larger data sets like those provided by radio telescopes such as MeerKAT. The pipeline also delivers a set of imaging quality metrics that give the user the opportunity to efficiently assess the data quality., Comment: 4 pages, 1 figure, in conference proceedings of the ADASS XXIX, to appear in ASPC. CARACal (formerly MeerKATHI) can be accessed at: https://caracal.readthedocs.io

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., Comment: 15 pages, 15 figures, 1 table, Accepted to MNRAS

This paper presents the results from one of the first observations of ionospheric scintillation taken using the Low-Frequency Array (LOFAR). The observation was of the strong natural radio source Cas A, taken overnight on 18-19 August 2013, and exhibited moderately strong scattering effects in dynamic spectra of intensity received across an observing bandwidth of 10-80MHz. Delay-Doppler spectra (the 2-D FFT of the dynamic spectrum) from the first hour of observation showed two discrete parabolic arcs, one with a steep curvature and the other shallow, which can be used to provide estimates of the distance to, and velocity of, the scattering plasma. A cross-correlation analysis of data received by the dense array of stations in the LOFAR "core" reveals two different velocities in the scintillation pattern: a primary velocity of ~30m/s with a north-west to south-east direction, associated with the steep parabolic arc and a scattering altitude in the F-region or higher, and a secondary velocity of ~110m/s with a north-east to south-west direction, associated with the shallow arc and a scattering altitude in the D-region. Geomagnetic activity was low in the mid-latitudes at the time, but a weak sub-storm at high latitudes reached its peak at the start of the observation. An analysis of Global Navigation Satellite Systems (GNSS) and ionosonde data from the time reveals a larger-scale travelling ionospheric disturbance (TID), possibly the result of the high-latitude activity, travelling in the north-west to south-east direction, and, simultaneously, a smaller--scale TID travelling in a north-east to south-west direction, which could be associated with atmospheric gravity wave activity. The LOFAR observation shows scattering from both TIDs, at different altitudes and propagating in different directions. To the best of our knowledge this is the first time that such a phenomenon has been reported., Comment: 24 pages, 16 figures. Accepted for open-access publication in the Journal of Space Weather and Space Climate. For associated movie file, see https://www.swsc-journal.org/10.1051/swsc/2020010/olm

We leverage the Sparsity Averaging Reweighted Analysis (SARA) approach for interferometric imaging, that is based on convex optimisation, for the super-resolution of Cyg A from observations at the frequencies 8.422GHz and 6.678GHz with the Karl G. Jansky Very Large Array (VLA). The associated average sparsity and positivity priors enable image reconstruction beyond instrumental resolution. An adaptive Preconditioned Primal-Dual algorithmic structure is developed for imaging in the presence of unknown noise levels and calibration errors. We demonstrate the superior performance of the algorithm with respect to the conventional CLEAN-based methods, reflected in super-resolved images with high fidelity. The high resolution features of the recovered images are validated by referring to maps of Cyg A at higher frequencies, more precisely 17.324GHz and 14.252GHz. We also confirm the recent discovery of a radio transient in Cyg A, revealed in the recovered images of the investigated data sets. Our matlab code is available online on GitHub., Comment: 14 pages, 7 figures (3/7 animated figures), accepted for publication in MNRAS

The Radio Interferometer Measurement Equation (RIME) is a matrix-based mathematical model that describes the response of a radio interferometer. The Jones calculus it employs is not suitable for describing the analogue components of a telescope. This is because it does not consider the effect of impedance mismatches between components. This paper aims to highlight the limitations of Jones calculus, and suggests some alternative methods that are more applicable. We reformulate the RIME with a different basis that includes magnetic and mixed coherency statistics. We present a microwave network inspired 2N-port version of the RIME, and a tensor formalism based upon the electromagnetic tensor from special relativity. We elucidate the limitations of the Jones-matrix-based RIME for describing analogue components. We show how measured scattering parameters of analogue components can be used in a 2N-port version of the RIME. In addition, we show how motion at relativistic speed affects the observed flux. We present reformulations of the RIME that correctly account for magnetic field coherency. These reformulations extend the standard formulation, highlight its limitations, and may have applications in space-based interferometry and precise absolute calibration experiments., Comment: Accepted for publication in MNRAS

The radio interferometer measurement equation (RIME), especially in its 2x2 form, has provided a comprehensive matrix-based formalism for describing classical radio interferometry and polarimetry, as shown in the previous three papers of this series. However, recent practical and theoretical developments, such as phased array feeds (PAFs), aperture arrays (AAs) and wide-field polarimetry, are exposing limitations of the formalism. This paper aims to develop a more general formalism that can be used to both clearly define the limitations of the matrix RIME, and to describe observational scenarios that lie outside these limitations. Some assumptions underlying the matrix RIME are explicated and analysed in detail. To this purpose, an array correlation matrix (ACM) formalism is explored. This proves of limited use; it is shown that matrix algebra is simply not a sufficiently flexible tool for the job. To overcome these limitations, a more general formalism based on tensors and the Einstein notation is proposed and explored both theoretically, and with a view to practical implementations. The tensor formalism elegantly yields generalized RIMEs describing beamforming, mutual coupling, and wide-field polarimetry in one equation. It is shown that under the explicated assumptions, tensor equations reduce to the 2x2 RIME. From a practical point of view, some methods for implementing tensor equations in an optimal way are proposed and analysed. The tensor RIME is a powerful means of describing observational scenarios not amenable to the matrix RIME. Even in cases where the latter remains applicable, the tensor formalism can be a valuable tool for understanding the limits of such applicability., Comment: 16 pages, no figures, accepted by A&A

Paper I of the series re-derived the radio interferometry measurement equation (RIME) from first principles, and extended the Jones formalism to the full-sky case, incorporating direction-dependent effects (DDEs). This paper aims to describe both classical radio interferometric calibration (selfcal and related methods), and the recent developments in the treatment of DDEs, using the RIME-based mathematical framework developed in Paper I. It also aims to demonstrate the ease with which the various effects can be described and understood. The first section of this paper uses the RIME formalism to describe self-calibration, both with a full RIME, and with the approximate equations of older software packages, and shows how this is affected by DDEs. The second section gives an overview of real-life DDEs and proposed methods of dealing with them. This results in a formal RIME-based description and comparison of existing and proposed approaches to the problem of DDEs., Comment: 11 pages, no figures, accepted by A&A

Since its formulation by Hamaker et al., the radio interferometer measurement equation (RIME) has provided a rigorous mathematical basis for the development of novel calibration methods and techniques, including various approaches to the problem of direction-dependent effects (DDEs). This series of papers aims to place recent developments in the treatment of DDEs into one RIME-based mathematical framework, and to demonstrate the ease with which the various effects can be described and understood. It also aims to show the benefits of a RIME-based approach to calibration. Paper I re-derives the RIME from first principles, extends the formalism to the full-sky case, and incorporates DDEs. Paper II then uses the formalism to describe self-calibration, both with a full RIME, and with the approximate equations of older software packages, and shows how this is affected by DDEs. It also gives an overview of real-life DDEs and proposed methods of dealing with them. Applying this to WSRT data (Paper III) results in a noise-limited image of the field around 3C 147 with a very high dynamic range (1.6 million), and none of the off-axis artifacts that plague regular selfcal. The resulting differential gain solutions contain significant information on DDEs, and can be used for iterative improvements of sky models. Perhaps most importantly, sources as faint as 2 mJy have been shown to yield meaningful differential gain solutions, and thus can be used as potential calibration beacons in other DDE-related schemes., Comment: 12 pages, no figures, published in A&A

The formulation of the radio interferometer measurement equation (RIME) by Hamaker et al. has provided us with an elegant mathematical apparatus for better understanding, simulation and calibration of existing and future instruments. The calibration of the new radio telescopes (LOFAR, SKA) would be unthinkable without the RIME formalism, and new software to exploit it. MeqTrees is designed to implement numerical models such as the RIME, and to solve for arbitrary subsets of their parameters. The technical goal of MeqTrees is to provide a tool for rapid implementation of such models, while offering performance comparable to hand-written code. We are also pursuing the wider goal of increasing the rate of evolution of radio astronomical software, by offering a tool for rapid experimentation and exchange of ideas. MeqTrees is implemented as a Python-based front-end called the meqbrowser, and an efficient (C++-based) computational back-end called the meqserver. Numerical models are defined on the front-end via a Python-based Tree Definition Language (TDL), then rapidly executed on the back-end. The use of TDL facilitates an extremely short turn-around time for experimentation with new ideas. This is also helped by unprecedented visualization capabilities for all final and intermediate results. A flexible data model and a number of important optimizations in the back-end ensures that the numerical performance is comparable to that of hand-written code. MeqTrees is already widely used as the simulation tool for new instruments (LOFAR, SKA) and technologies (focal plane arrays). It has demonstrated that it can achieve a noise-limited dynamic range in excess of a million, on WSRT data. It is the only package that is specifically designed to handle what we propose to call third-generation calibration (3GC), which is needed for the new generation of giant radio telescopes., Comment: 15 pages; 14 figures

Papers I and II of this series have extended the radio interferometry measurement equation (RIME) formalism to the full-sky case, and provided a RIME-based description of calibration and the problem of direction-dependent effects (DDEs). This paper aims to provide a practical demonstration of a RIME-based approach to calibration, via an example of extremely high-dynamic range calibration of WSRT observations of 3C 147 at 21 cm, with full treatment of DDEs. A version of the RIME incorporating differential gains has been implemented in MeqTrees, and applied to the 3C 147 data. This was used to perform regular selfcal, then solve for interferometer-based errors and for differential gains. The resulting image of the field around 3C 147 is thermal noise-limited, has a very high dynamic range (1.6 million), and none of the off-axis artefacts that plague regular selfcal. The differential gain solutions show a high signal-to-noise ratio, and may be used to extract information on DDEs and errors in the sky model. The differential gain approach can eliminate DDE-related artefacts, and provide information for iterative improvements of sky models. Perhaps most importantly, sources as faint as 2 mJy have been shown to yield meaningful differential gain solutions, and thus can be used as potential calibration beacons in other DDE-related schemes., Comment: 11 pages, 17 figures, accepted by A&A