The data produced by the new generation of interferometers are affected by a wide variety of partially unknown complex effects such as pointing errors, phased array beams, ionosphere, troposphere, ...Faraday rotation, or clock drifts. Most algorithms addressing direction-dependent calibration solve for the effective Jones matrices, and cannot constrain the underlying physical quantities of the radio interferometry measurement equation (RIME). A related difficulty is that they lack robustness in the presence of low signal-to-noise ratios, and when solving for moderate to large numbers of parameters they can be subject to ill-conditioning. These effects can have dramatic consequences in the image plane such as source or even thermal noise suppression. The advantage of solvers directly estimating the physical terms appearing in the RIME is that they can potentially reduce the number of free parameters by orders of magnitudes while dramatically increasing the size of usable data, thereby improving conditioning. We present here a new calibration scheme based on a nonlinear version of the Kalman filter that aims at estimating the physical terms appearing in the RIME. We enrich the filter’s structure with a tunable data representation model, together with an augmented measurement model for regularization. Using simulations we show that it can properly estimate the physical effects appearing in the RIME. We found that this approach is particularly useful in the most extreme cases such as when ionospheric and clock effects are simultaneously present. Combined with the ability to provide prior knowledge on the expected structure of the physical instrumental effects (expected physical state and dynamics), we obtain a fairly computationally cheap algorithm that we believe to be robust, especially in low signal-to-noise regimes. Potentially, the use of filters and other similar methods can represent an improvement for calibration in radio interferometry, under the condition that the effects corrupting visibilities are understood and analytically stable. Recursive algorithms are particularly well adapted for pre-calibration and sky model estimate in a streaming way. This may be useful for the SKA-type instruments that produce huge amounts of data that have to be calibrated before being averaged.
Recent developments in optimization theory have extended some traditional algorithms for least-squares optimization of real-valued functions (Gauss–Newton, Levenberg–Marquardt, etc.) into the domain ...of complex functions of a complex variable. This employs a formalism called the Wirtinger derivative, and derives a full-complex Jacobian counterpart to the conventional real Jacobian. We apply these developments to the problem of radio interferometric gain calibration, and show how the general complex Jacobian formalism, when combined with conventional optimization approaches, yields a whole new family of calibration algorithms, including those for the polarized and direction-dependent gain regime. We further extend the Wirtinger calculus to an operator-based matrix calculus for describing the polarized calibration regime. Using approximate matrix inversion results in computationally efficient implementations; we show that some recently proposed calibration algorithms such as StefCal and peeling can be understood as special cases of this, and place them in the context of the general formalism. Finally, we present an implementation and some applied results of CohJones, another specialized direction-dependent calibration algorithm derived from the formalism.
The new generation of radio interferometers is characterized by high sensitivity, wide fields of view and large fractional bandwidth. To synthesize the deepest images enabled by the high dynamic ...range of these instruments requires us to take into account the direction-dependent Jones matrices, while estimating the spectral properties of the sky in the imaging and deconvolution algorithms. In this paper we discuss and implement a wideband wide-field spectral deconvolution framework (ddfacet) based on image plane faceting, that takes into account generic direction-dependent effects. Specifically, we present a wide-field co-planar faceting scheme, and discuss the various effects that need to be taken into account to solve for the deconvolution problem (image plane normalization, position-dependent Point Spread Function, etc). We discuss two wideband spectral deconvolution algorithms based on hybrid matching pursuit and sub-space optimisation respectively. A few interesting technical features incorporated in our imager are discussed, including baseline dependent averaging, which has the effect of improving computing efficiency. The version of ddfacet presented here can account for any externally defined Jones matrices and/or beam patterns.
Abstract
Radio emission is a key indicator of star formation activity in galaxies, but the radio luminosity–star formation relation has to date been studied almost exclusively at frequencies of ...1.4 GHz or above. At lower radio frequencies, the effects of thermal radio emission are greatly reduced, and so we would expect the radio emission observed to be completely dominated by synchrotron radiation from supernova-generated cosmic rays. As part of the LOFAR Surveys Key Science project, the Herschel-ATLAS NGP field has been surveyed with LOFAR at an effective frequency of 150 MHz. We select a sample from the MPA-JHU catalogue of Sloan Digital Sky Survey galaxies in this area: the combination of Herschel, optical and mid-infrared data enable us to derive star formation rates (SFRs) for our sources using spectral energy distribution fitting, allowing a detailed study of the low-frequency radio luminosity–star formation relation in the nearby Universe. For those objects selected as star-forming galaxies (SFGs) using optical emission line diagnostics, we find a tight relationship between the 150 MHz radio luminosity (L150) and SFR. Interestingly, we find that a single power-law relationship between L150 and SFR is not a good description of all SFGs: a broken power-law model provides a better fit. This may indicate an additional mechanism for the generation of radio-emitting cosmic rays. Also, at given SFR, the radio luminosity depends on the stellar mass of the galaxy. Objects that were not classified as SFGs have higher 150-MHz radio luminosity than would be expected given their SFR, implying an important role for low-level active galactic nucleus activity.
ABSTRACT
Measuring the properties of extragalactic magnetic fields through the effect of Faraday rotation provides a means to understand the origin and evolution of cosmic magnetism. Here, we use ...data from the LOFAR Two-Metre Sky Survey (LoTSS) to calculate the Faraday rotation measure (RM) of close pairs of extragalactic radio sources. By considering the RM difference (ΔRM) between physical pairs (e.g. double-lobed radio galaxies) and non-physical pairs (i.e. close projected sources on the sky), we statistically isolate the contribution of extragalactic magnetic fields to ΔRM along the line of sight between non-physical pairs. From our analysis, we find no significant difference between the ΔRM distributions of the physical and non-physical pairs, limiting the excess Faraday rotation contribution to <1.9 rad m−2 (${\sim}95{{\ \rm per\ cent}}$ confidence). We use this limit with a simple model of an inhomogeneous universe to place an upper limit of 4 nG on the cosmological co-moving magnetic field strength on Mpc scales. We also compare the RM data with a more realistic suite of cosmological magnetohydrodynamical simulations that explore different magnetogenesis scenarios. Both magnetization of the large-scale structure by astrophysical processes such as galactic and AGN outflows, and simple primordial scenarios with seed magnetic field strengths <0.5 nG cannot be rejected by the current data; while stronger primordial fields or models with dynamo amplification in filaments are disfavoured.
ABSTRACT
The effects of feedback from high luminosity radio-loud active galactic nuclei (AGN) have been extensively discussed in the literature, but feedback from low-luminosity radio-loud AGN is ...less well understood. The advent of high sensitivity, high angular resolution, large field-of-view telescopes such as LOFAR is now allowing wide-area studies of such faint sources for the first time. Using the first data release of the LOFAR Two Metre Sky Survey (LoTSS) we report on our discovery of a population of 195 radio galaxies with 150-MHz luminosities between 3 × 1022 and 1.5 × 1025 WHz−1 and total radio emission no larger than 80 kpc. These objects, which we term galaxy-scale jets (GSJ), are small enough to be directly influencing the evolution of the host on galaxy scales. We report upon the typical host properties of our sample, finding that 9 per cent are hosted by spirals with the remainder being hosted by elliptical galaxies. Two of the spiral-hosted GSJ are highly unusual with low radio luminosities and FRII-like morphology. The host properties of our GSJ show that they are ordinary AGN observed at a stage in their life shortly after the radio emission has expanded beyond the central regions of the host. Based on our estimates, we find that about half of our GSJ have internal radio lobe energy within an order of magnitude of the ISM energy so that, even ignoring any possible shocks, GSJ are energetically capable of affecting the evolution of the host. The current sample of GSJ will grow in size with future releases of LoTSS and can also form the basis for further studies of feedback from low-luminosity radio sources.
The Galactic Faraday rotation sky 2020 Hutschenreuter, S.; Anderson, C. S.; Betti, S. ...
Astronomy & astrophysics,
01/2022, Volume:
657
Journal Article
Peer reviewed
Open access
Aims.
This work provides an update to existing reconstructions of the Galactic Faraday rotation sky by processing almost all Faraday rotation data sets available at the end of the year 2020. ...Observations of extra-Galactic sources in recent years have further illuminated the previously underconstrained southern celestial sky, as well as parts of the inner disc of the Milky Way, along with other regions. This has culminated in an all-sky data set of 55 190 data points, thereby comprising a significant expansion on the 41 330 used in previous works. At the same time, this novelty makes an updated separation of the Galactic component a promising enterprise. The increased source density allows us to present our results in a resolution of about 1.3 × 10
−2
deg
2
(46.8 arcmin
2
), which is a twofold increase compared to previous works.
Methods.
As for previous Faraday rotation sky reconstructions, this work is based on information field theory, namely, a Bayesian inference scheme for field-like quantities that handles noisy and incomplete data.
Results.
In contrast to previous reconstructions, we find a significantly thinner and pronounced Galactic disc with small-scale structures exceeding values of several thousand rad m
−2
. The improvements can mainly be attributed to the new catalog of Faraday data, but are also supported by advances in correlation structure modeling within numerical information field theory. We also provide a detailed discussion on the statistical properties of the Faraday rotation sky and we investigate correlations with other data sets.
The LOFAR Two-metre Sky Survey Shimwell, T. W.; Röttgering, H. J. A.; Best, P. N. ...
Astronomy and astrophysics (Berlin),
02/2017, Volume:
598
Journal Article
Peer reviewed
Open access
The LOFAR Two-metre Sky Survey (LoTSS) is a deep 120–168 MHz imaging survey that will eventually cover the entire northern sky. Each of the 3170 pointings will be observed for 8 h, which, at most ...declinations, is sufficient to produce ∼5 resolution images with a sensitivity of ∼100 µJy/beam and accomplish the main scientific aims of the survey, which are to explore the formation and evolution of massive black holes, galaxies, clusters of galaxies and large-scale structure. Owing to the compact core and long baselines of LOFAR, the images provide excellent sensitivity to both highly extended and compact emission. For legacy value, the data are archived at high spectral and time resolution to facilitate subarcsecond imaging and spectral line studies. In this paper we provide an overview of the LoTSS. We outline the survey strategy, the observational status, the current calibration techniques, a preliminary data release, and the anticipated scientific impact. The preliminary images that we have released were created using a fully automated but direction-independent calibration strategy and are significantly more sensitive than those produced by any existing large-area low-frequency survey. In excess of 44 000 sources are detected in the images that have a resolution of 25 , typical noise levels of less than 0.5 mJy/beam, and cover an area of over 350 square degrees in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45 • 00 00 to 57 • 00 00).
ABSTRACT
We report the discovery of two new giant radio galaxies (GRGs) using the MeerKAT International GHz Tiered Extragalactic Exploration (MIGHTEE) survey. Both GRGs were found within a ...${\sim}1\,$ deg2 region inside the COSMOS field. They have redshifts of z = 0.1656 and z = 0.3363 and physical sizes of 2.4 and 2.0 Mpc, respectively. Only the cores of these GRGs were clearly visible in previous high-resolution Very Large Array observations, since the diffuse emission of the lobes was resolved out. However, the excellent sensitivity and uv coverage of the new MeerKAT telescope allowed this diffuse emission to be detected. The GRGs occupy an unpopulated region of radio power – size parameter space. Based on a recent estimate of the GRG number density, the probability of finding two or more GRGs with such large sizes at z < 0.4 in a ${\sim}1\,$ deg2 field is only 2.7 × 10−6, assuming Poisson statistics. This supports the hypothesis that the prevalence of GRGs has been significantly underestimated in the past due to limited sensitivity to low surface brightness emission. The two GRGs presented here may be the first of a new population to be revealed through surveys like MIGHTEE that provide exquisite sensitivity to diffuse, extended emission.
Abstract
This paper presents a study of the redshift evolution of radio-loud active galactic nuclei (AGN) as a function of the properties of their galaxy hosts in the Boötes field. To achieve this we ...match low-frequency radio sources from deep 150-MHz LOFAR (LOw Frequency ARray) observations to an I-band-selected catalogue of galaxies, for which we have derived photometric redshifts, stellar masses, and rest-frame colours. We present spectral energy distribution (SED) fitting to determine the mid-infrared AGN contribution for the radio sources and use this information to classify them as high- versus low-excitation radio galaxies (HERGs and LERGs) or star-forming galaxies. Based on these classifications, we construct luminosity functions for the separate redshift ranges going out to z = 2. From the matched radio-optical catalogues, we select a sub-sample of 624 high power (P150 MHz > 1025 W Hz−1) radio sources between 0.5 ≤ z < 2. For this sample, we study the fraction of galaxies hosting HERGs and LERGs as a function of stellar mass and host galaxy colour. The fraction of HERGs increases with redshift, as does the fraction of sources in galaxies with lower stellar masses. We find that the fraction of galaxies that host LERGs is a strong function of stellar mass as it is in the local Universe. This, combined with the strong negative evolution of the LERG luminosity functions over this redshift range, is consistent with LERGs being fuelled by hot gas in quiescent galaxies.