Source author record

Cyril Tasse

Cyril Tasse appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2022arXiv

Faraday tomography of LoTSS-DR2 data: I. Faraday moments in the high-latitude outer Galaxy and revealing Loop III in polarisation

Observations of synchrotron emission at low radio frequencies reveal a labyrinth of polarised Galactic structures. However, the explanation for the wealth of structures remains uncertain due to the complex interactions between the interstellar medium and the magnetic field. A multi-tracer approach to the analysis of large sky areas is needed. This paper aims to use polarimetric images from the LOFAR Two metre Sky Survey (LoTSS) to produce the biggest mosaic of polarised emission in the northern sky at low radio frequencies (150 MHz) to date. The large area this mosaic covers allows for detailed morphological and statistical studies of polarised structures in the high-latitude outer Galaxy, including the well-known Loop III region. We produced a 3100 square degree Faraday tomographic cube using a rotation measure synthesis tool. We calculated the statistical moments of Faraday spectra and compared them with data sets at higher frequencies (1.4 GHz) and with a map of a rotation measure derived from extragalactic sources. The mosaic is dominated by polarised emission connected to Loop III. Additionally, the mosaic reveals an abundance of other morphological structures, mainly {narrow and extended} depolarisation canals, which are found to be ubiquitous. We find a correlation between the map of an extragalactic rotation measure and the LoTSS first Faraday moment image. The ratio of the two deviates from a simple model of a Burn slab (Burn 1966) along the line of sight, which highlights the high level of complexity in the magnetoionic medium that can be studied at these frequencies.

preprint2022arXiv

The discovery of a radio galaxy of at least 5 Mpc

We discover what is in projection the largest known structure of galactic origin: a giant radio galaxy with a projected proper length of $4.99 \pm 0.04\ \mathrm{Mpc}$. The source, named Alcyoneus, was first identified in low-resolution LOFAR Two-metre Sky Survey images from which angularly compact sources had been removed. Being an extreme example in its class, Alcyoneus could shed light on the main mechanisms that drive radio galaxy growth. We find that - beyond geometry - Alcyoneus and its host galaxy appear suspiciously ordinary: the total low-frequency luminosity density, stellar mass and supermassive black hole mass are all lower than, though similar to, those of the medial giant radio galaxy (percentiles $45 \pm 3\%$, $25 \pm 9 \%$ and $23 \pm 11 \%$, respectively). The source resides in a filament of the Cosmic Web, with which it might have significant thermodynamic interaction. At $5 \cdot 10^{-16}\ \mathrm{Pa}$, the pressures in the lobes are the lowest hitherto found, and Alcyoneus therefore represents one of the most promising radio galaxies yet to probe the warm-hot intergalactic medium.

preprint2021arXiv

Spectral analysis of spatially-resolved 3C295 (sub-arcsecond resolution) with the International LOFAR Telescope

3C295 is a bright, compact steep spectrum source with a well-studied integrated radio spectral energy distribution (SED) from 132 MHz to 15 GHz. However, spatially resolved spectral studies have been limited due to a lack of high resolution images at low radio frequencies. These frequencies are crucial for measuring absorption processes, and anchoring the overall spectral modelling of the radio SED. In this paper, we use International LOFAR (LOw-Frequency ARray) Telescope (ILT) observations of 3C295 to study its spatially resolved spectral properties with sub-arcsecond resolution at 132 MHz. Combining our new 132 MHz observation with archival data at 1.6 GHz, 4.8 GHz, and 15 GHz, we are able to carry out a resolved radio spectral analysis. The spectral properties of the hotspots provides evidence for low frequency flattening. In contrast, the spectral shape across the lobes is consistent with a JP spectral ageing model. Using the integrated spectral information for each component, we then fit low-frequency absorption models to the hotspots, finding that both free-free absorption and synchrotron self-absorption models provide a better fit to the data than a standard power law. Although we can say there is low-frequency absorption present in the hotspots of 3C295, future observations with the Low Band Antenna of the ILT at 55 MHz may allow us to distinguish the type of absorption.

preprint2021arXiv

The Galactic Faraday rotation sky 2020

This work gives 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, among other regions, further illuminated the previously under-constrained southern celestial sky, as well as parts of the inner disc of the Milky Way. This has culminated in an all-sky data set of 55,190 data points, which is a significant expansion on the 41,330 used in previous works, hence making an updated separation of the Galactic component a promising venture. The increased source density allows us to present our results in a resolution of about $1.3\cdot 10^{-2}\, \mathrm{deg}^2$ ($46.8\,\mathrm{arcmin}^2$), which is a twofold increase compared to previous works. As for previous Faraday rotation sky reconstructions, this work is based on information field theory, a Bayesian inference scheme for field-like quantities which handles noisy and incomplete data. In contrast to previous reconstructions, we find a significantly thinner and pronounced Galactic disc with small-scale structures exceeding values of several thousand $\mathrm{rad}\,\mathrm{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 furthermore give a detailed discussion on statistical properties of the Faraday rotation sky and investigate correlations to other data sets.

preprint2020arXiv

ASKAP reveals giant radio halos in two merging SPT galaxy clusters -- Making the case for a direction-dependent pipeline --

Early science observations from the Australian Square Kilometre Array Pathfinder (ASKAP) have revealed clear signals of diffuse radio emission associated with two clusters detected by the South Pole Telescope via their Sunyaev Zel'dovich signal. SPT CLJ0553-3342 (MACSJ0553.4-3342) and SPT CLJ0638-5358 (Abell S0592) are both high-mass lensing clusters that have undergone major mergers. To improve the data products of these ASKAP early science observations and create science-fidelity images of the galaxy clusters, we performed direction-dependent (DD) calibration and imaging using state-of-the-art software {\sc killMS} and {\sc DDFacet}. We find that artefacts in the ASKAP images are greatly reduced after directional calibration. Here we present our DD calibrated ASKAP radio images of both clusters showing unambiguous giant radio halos with largest linear scales of $\sim1$~Mpc. The halo in MACSJ0553.4-3342 was previously detected with GMRT observations at 323 MHz, but appears more extended in our ASKAP image. Although there is a shock detected in the thermal X-ray emission of this cluster, we find that the particle number density in the shocked region is too low to allow for the generation of a radio shock. The radio halo in Abell S0592 is a new discovery, and the Southwest border of the halo coincides with a shock detected in X-rays. We discuss the origins of these halos considering both the hadronic and turbulent re-acceleration models as well as sources of \textit{seed} electrons. This work gives a positive indication of the potential of ASKAP's Evolutionary Map of the Universe (EMU) survey in detecting intracluster medium radio sources, and showcases the improvement in data products after utilising third-generation calibration techniques.

preprint2015arXiv

KAT-7 detection of radio halo emission in the Triangulum Australis galaxy cluster

We report the presence of high significance diffuse radio emission from the Triangulum Australis cluster using observations made with the KAT-7 telescope and propose that this emission is a giant radio halo. We compare the radio power from this proposed halo with X-ray and SZ measurements and demonstrate that it is consistent with the established scaling relations for cluster haloes. By combining the X-ray and SZ data we calculate the ratio of non-thermal to thermal electron pressure within Triangulum Australis to be $X=0.658\pm0.054$. We use this ratio to constrain the maximum magnetic field strength within the halo region to be $B_{\rm max, halo} = 33.08\,μ$G and compare this with the minimum field strength from equipartition of $B_{\rm min, halo} = 0.77(1+k)^{2/7}\,μ$G to place limits on the range of allowed magnetic field strength within this cluster. We compare these values to those for more well-studied systems and discuss these results in the context of equipartition of non-thermal energy densities within clusters of galaxies.

preprint2015arXiv

Montblanc: GPU accelerated Radio Interferometer Measurement Equations in support of Bayesian Inference for Radio Observations

We present Montblanc, a GPU implementation of the Radio interferometer measurement equation (RIME) in support of the Bayesian inference for radio observations (BIRO) technique. BIRO uses Bayesian inference to select sky models that best match the visibilities observed by a radio interferometer. To accomplish this, BIRO evaluates the RIME multiple times, varying sky model parameters to produce multiple model visibilities. Chi-squared values computed from the model and observed visibilities are used as likelihood values to drive the Bayesian sampling process and select the best sky model. As most of the elements of the RIME and chi-squared calculation are independent of one another, they are highly amenable to parallel computation. Additionally, Montblanc caters for iterative RIME evaluation to produce multiple chi-squared values. Modified model parameters are transferred to the GPU between each iteration. We implemented Montblanc as a Python package based upon NVIDIA's CUDA architecture. As such, it is easy to extend and implement different pipelines. At present, Montblanc supports point and Gaussian morphologies, but is designed for easy addition of new source profiles. Montblanc's RIME implementation is performant: On an NVIDIA K40, it is approximately 250 times faster than MeqTrees on a dual hexacore Intel E5-2620v2 CPU. Compared to the OSKAR simulator's GPU-implemented RIME components it is 7.7 and 12 times faster on the same K40 for single and double-precision floating point respectively. However, OSKAR's RIME implementation is more general than Montblanc's BIRO-tailored RIME. Theoretical analysis of Montblanc's dominant CUDA kernel suggests that it is memory bound. In practice, profiling shows that is balanced between compute and memory, as much of the data required by the problem is retained in L1 and L2 cache.

preprint2015arXiv

Radio interferometric gain calibration as a complex optimization problem

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.

preprint2015arXiv

Sparse representations and convex optimization as tools for LOFAR radio interferometric imaging

Compressed sensing theory is slowly making its way to solve more and more astronomical inverse problems. We address here the application of sparse representations, convex optimization and proximal theory to radio interferometric imaging. First, we expose the theory behind interferometric imaging, sparse representations and convex optimization, and second, we illustrate their application with numerical tests with SASIR, an implementation of the FISTA, a Forward-Backward splitting algorithm hosted in a LOFAR imager. Various tests have been conducted in Garsden et al., 2015. The main results are: i) an improved angular resolution (super resolution of a factor ~2) with point sources as compared to CLEAN on the same data, ii) correct photometry measurements on a field of point sources at high dynamic range and iii) the imaging of extended sources with improved fidelity. SASIR provides better reconstructions (five time less residuals) of the extended emission as compared to CLEAN. With the advent of large radiotelescopes, there is scope for improving classical imaging methods with convex optimization methods combined with sparse representations.

preprint2014arXiv

Applying Wirtinger derivatives to the radio interferometry calibration problem

This paper presents a fast algorithm for full-polarisation, direction dependent calibration in radio interferometry. It is based on Wirtinger's approach to complex differentiation. Compared to the classical case, and under reasonable approximations, the Jacobian appearing in the Levenberg-Maquardt iterative scheme presents a sparse structure, allowing for high gain in terms of algorithmic cost.

preprint2014arXiv

Non-linear Kalman filters for calibration in radio interferometry

We present a new calibration scheme based on a non-linear version of Kalman filter that aims at estimating the physical terms appearing in the Radio Interferometry Measurement Equation (RIME). We enrich the filter's structure with a tunable data representation model, together with an augmented measurement model for regularization. We show using simulations 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 cheap algorithm that we believe to be robust, especially in low signal-to-noise regime. 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.