Source author record

Oleg M. Smirnov

Oleg M. Smirnov 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

16works
6topics
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

16 published item(s)

preprint2022arXiv

The extended HI halo of NGC 4945 as seen by MeerKAT

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σ$ 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.

preprint2022arXiv

The Hydrogen Intensity and Real-time Analysis eXperiment: 256-Element Array Status and Overview

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.

preprint2022arXiv

Tricolour: an optimized SumThreshold flagger for MeerKAT

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.

preprint2021arXiv

A practical preconditioner for wide-field continuum imaging of radio interferometric data

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.

preprint2021arXiv

Anomalous gas in ESO 149-G003: A MeerKAT-16 View

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.

preprint2021arXiv

The detection of a massive chain of dark HI clouds in the GAMA G23 Field

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.

preprint2020arXiv

A LOFAR Observation of Ionospheric Scintillation from Two Simultaneous Travelling Ionospheric Disturbances

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.

preprint2020arXiv

Foreground modelling via Gaussian process regression: an application to HERA data

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.

preprint2020arXiv

MeerKAT-16 HI observation of the dIrr galaxy WLM

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.

preprint2020arXiv

MeerKATHI -- an end-to-end data reduction pipeline for MeerKAT and other radio telescopes

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.

preprint2015arXiv

Generalized Formalisms of the Radio Interferometer Measurement Equation

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.

preprint2011arXiv

Revisiting the radio interferometer measurement equation. I. A full-sky Jones formalism

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.

preprint2011arXiv

Revisiting the radio interferometer measurement equation. II. Calibration and direction-dependent effects

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.

preprint2011arXiv

Revisiting the radio interferometer measurement equation. III. Addressing direction-dependent effects in 21 cm WSRT observations of 3C 147

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.

preprint2011arXiv

Revisiting the radio interferometer measurement equation. IV. A generalized tensor formalism

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.

preprint2011arXiv

The MeqTrees software system and its use for third-generation calibration of radio interferometers

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.