Open Access
Issue
A&A
Volume 678, October 2023
Article Number A72
Number of page(s) 13
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/202347320
Published online 06 October 2023

© The Authors 2023

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model.

Open Access funding provided by Max Planck Society.

1 Introduction

The Galactic mid-plane of the Milky Way is populated with thermal and nonthermal radio sources originating during the evolution of short-lived, high-mass stars (M > 8 M), such as H II regions and supernova remnants (SNRs). The formation of such high-mass stars is controlled by the physical and chemical conditions in the interstellar medium (ISM) that surrounds them, which is affected by these objects themselves, creating a complex feedback loop that is not yet well understood. Observations of the ISM continue to shed light on its structure and composition and help us understand our Galaxy better by constraining theoretical models on the evolution of its contents (e.g., Soler et al. 2020; Katz et al. 2021; Yan et al. 2023). In the past, a plethora of surveys, from sub-millimetre to metre wavelengths, have studied objects arising from various stages of stellar evolution, probing angular scales ranging from a few arcseconds to several degrees. Surveys in the radio regime are particularly helpful as dust is optically thin at low radio frequencies. The Multi-Array Galactic Plane Imaging Survey (MAGPIS; Becker et al. 1994), the APEX Telescope Large Area Survey of the Galaxy (ATLAS-GAL; Schuller et al. 2009), the Coordinated Radio and Infrared Survey for High-Mass Star Formation (CORNISH; Hoare et al. 2012), the Very Large Array Galactic Plane Survey (VGPS; Stil et al. 2006), the 1–2 GHz HI/OH/Recombination line survey (THOR; Beuther et al. 2016; Wang et al. 2020a), and the 4–8 GHz GLObal view on STAR formation survey (GLOSTAR; Medina et al. 2019; Brunthaler et al. 2021) are a few such surveys conducting observations of the northern Galactic plane.

The radio emission at frequencies below 1 GHz is dominated by nonthermal synchrotron radiation, which is observed not only from SNRs (e.g., Wilson et al. 2013) but also from some star forming regions (e.g., Wilner et al. 1999; Veena et al. 2019). The Galactic plane at these frequencies, however, is relatively unexplored at near-arcsecond angular resolution. The GaLactic and Extragalactic All-sky Murchison Widefield Array survey (GLEAM; Hurley-Walker et al. 2019) at 70–230 MHz covers a large portion of the Galactic plane at 2′–4′ resolution and a sensitivity of about 10–20 mJy beam−1. The 150 MHz TIFR-GMRT Sky Survey (TGSS) by Intema et al. (2017) covers the Galactic plane at a better resolution of 25″ and a typical sensitivity of about 5 mJy beam−1, but the extended emission is poorly recovered due to the imaging choices that were made with regards to their survey goals. The 843 MHz Molongolo Galactic plane survey (MGPS; Green et al. 1999) has a sensitivity of about 1–2 mJy beam−1, a resolution of ~43″, and is sensitive to extended structures up to a scale of about 25′, but it covers only the southern Galactic plane, which does not include the first quadrant of the Milky Way. There are multiple ongoing surveys with the Low Frequency Array below 200 MHz, but they cover only the northern sky (declination δ > 0°), which does not include the inner Galactic plane (e.g., Heald et al. 2015; Shimwell et al. 2017; de Gasperin et al. 2021).

The low-frequency regime (< 1 GHz) is particularly helpful to study nonthermal emission from SNRs, since they typically have a spectral index of α~ −0.51 and are brighter at lower frequencies. The number of SNRs identified in the Milky Way so far is lower than 400 (Ferrand & Safi-Harb 2012; Green 2019), which is substantially fewer than the predicted value of ~ 1000–2000 (Li et al. 1991; Ranasinghe & Leahy 2022). Recently, over 150 objects in the first quadrant of the Galactic plane have been identified as SNR candidates (Anderson et al. 2017; Dokara et al. 2021), using data from the THOR and the GLOSTAR surveys. Given that most of these candidates are of small angular size (Dokara et al. 2023) and are expected to have negative spectral indices, typically around −0.5, the next natural step is to search for these objects in the images of high resolution Galactic plane surveys at sub-GHz frequencies. In addition, due to synchrotron-ageing in old SNRs, the radio spectrum becomes steeper as the higher energy photons are lost much more easily than the lower energy photons. At frequencies below 100 MHz, synchrotron self-absorption effects become prominent and make it harder to study nonthermal emission. Some SNRs are known to have turnover frequencies even close to 1 GHz or higher (Sun et al. 2011; Kothes et al. 2020; Dokara et al. 2023). This makes the low-frequency radio window the most suitable to observe and study SNRs. The THOR and the GLOSTAR surveys have also been used to study the warm-ionized medium (Shanahan et al. 2019), cloud formation processes (Wang et al. 2020b), stellar feedback (Rugel et al. 2019), OH and methanol masers (Beuther et al. 2019; Ortiz-León et al. 2021; Nguyen et al. 2022), and the Galactic center star-formation rate (Nguyen et al. 2021), to name a few. However, these centimetre wavelength surveys of the Milky Way lack an equivalent survey below 1 GHz (at wavelengths longer than 30 cm) that is sensitive to large-scale emission at arcsecond-level resolution covering the first quadrant of the Galactic plane.

The Giant Metrewave Radio Telescope (GMRT) is an excellent instrument for this purpose. With its recent upgrade to wide-band receivers and correlator backend (called uGMRT from now on; Gupta et al. 2017), which can provide an instantaneous bandwidth up to 400 MHz, a sensitive survey is now feasible within a reasonable time-frame. The expected sizes of the synthesized beams of uGMRT observations conducted in band-3 (300–500 MHz) and band-4 (550–750 MHz) are about 5″ and 8″, respectively. Due to its unique hybrid configuration (see Fig. 1), extended structures up to scales of ~20′ (at 400 MHz) can also be recovered. Inferring the dominant emission mechanism of an object will be possible by comparing the data from the two frequency bands (or even at different frequencies within one band) and it will also be possible to draw broadband spectral energy distributions for objects with counterparts in the THOR and the GLOSTAR surveys.

To understand the feasibility of such a large-scale survey, we conducted uGMRT observations of the W43/W44 region (29° < l < 35° and |b| < 1°) using the band-3 and band-4 receivers as a part of a “pilot” study. The low-frequency radio regime (below 1 GHz) is particularly susceptible to man-made radio frequency interference (RFI) from consumer electronics and artificial satellites. Moreover, imaging structures larger than a few arcminutes with interferometers, at arcsecond resolution and poor uv-coverage (due to snapshot observations), is quite challenging. This problem gets worse at frequencies below 1 GHz, where the ionospheric conditions greatly affect the phase stability, making it difficult to achieve a good positional accuracy. The chosen region of the Galactic plane for this pilot study covers a region where the Scutum–Centaurus spiral arm meets the bar of the Milky Way. It is home to several bright extended Galactic sources, such as the W44 SNR (G34.7−0.4) and the W43 star-forming “mini-starburst” region (at l ~ 30°, b ~ 0°), and it contains other compact and ultracompact H II regions, in addition to numerous unresolved extragalactic sources (Chakraborty et al. 2020). This region has not yet been explored with the uGMRT. All this together makes it an excellent test-bed for our purposes. The rest of the paper is organized as follows: In Sect. 2, we describe the observations and the data reduction steps. In Sect. 3, we present the continuum and spectral index mosaics, along with the results of the analysis of positions and flux densities. Finally, the conclusions from this study and its implications for future low frequency Galactic plane surveys are discussed in Sect. 4.

thumbnail Fig. 1

Antenna configuration of the uGMRT array. The inset shows the central square (1.6 km × 1.6 km) antennas. North is upwards and east is rightwards.

thumbnail Fig. 2

Pointing configuration map of our pilot survey. Each large red circle marks the field of view of the pointing used for making a mosaic of the band-3 data, while the red cross shows the pointing center. Similarly, the band-4 fields of view and pointing centers are marked in blue.

2 Data

2.1 Observations

The data presented in this paper were obtained from observations with the uGMRT that were carried out under the proposal code 36_061. The total time allotted to the project was 26 h. We observed in multiple observing sessions in 2019 from May to August (see Table 1). We were unable to recover any useful data from the June (band-3) observation, during which we used a 100 MHz bandwidth, due to strong broad-band RFI from artificial satellites near the target sky positions. This prompted us to re-observe in band-3 with a larger bandwidth of 200 MHz to account for the sensitivity loss due to channels corrupted by RFI. These observations were done in August 2019.

The observed region covers the Galactic longitude range l = 29° – 35°, from the Galactic latitude b = −1° to b = +1°. We followed a hexagonal mapping pattern to observe the target region in several individual pointings. The pointing centers are shown in Fig. 2, along with the fields of view used for making a mosaic of the images. The fields of view that were chosen for each pointing in band-3 and band-4 are about 35′ and 18′, respectively, which correspond to primary beam attenuation factors of 0.36 and 0.64, respectively. Each target field was observed for ~14 min in total, done in two separate scans of 7 min each that were spaced across the observation time in order to achieve a slightly better uv-coverage. To establish the flux density scale, the standard primary calibrator 3C 286 was observed at least once during an observation, and the gain calibrator J1851+005 was observed about once every hour. The best continuum point source sensitivity achievable is expected to be about 150 μJy beam−1 and 60 μJy beam−1 in band-3 and band-4, respectively, with our integration time and usable bandwidth, based on the GMRT exposure time calculator2. The expected resolutions of the data are about 8″ and 5″ in band-3 and band-4, respectively. All uGMRT data are recorded in spectral line mode, typically in 8192 channels, which helps in isolating the RFI to a fraction of the total number of channels across the observed bandwidth. The steps to reduce the continuum data and produce mosaics are explained below.

2.2 Data reduction

The raw visibility data measured by an interferometer are corrupted due to variable response from the receiving instruments and the atmosphere (so-called “gains”). The process of calibration estimates the actual visibilities from the observed visibilities, which are related to each other by a measurement equation as described by Hamaker et al. (1996). Imaging the calibrated data constitutes the gridding of the visibilities, performing the inverse Fourier transform of the gridded visibilities, and removing the imprint of the point spread function (PSF, i.e., the “dirty beam”) of the interferometer. These two steps (calibration and imaging) are detailed in the following sections.

Table 1

MeGaPluG pilot survey observations.

thumbnail Fig. 3

650 MHz image of the gain calibrator field (J1851+005) of our observations. Circled in red at the center is the phase calibrator (J1851+005; Galactic coordinates l = 33.5°, b = 0.2°). The large and extended structure to its southeast is the supernova remnant (SNR) G33.6+0.1. The beam size is about 10 arcsec. The intensity scale is saturated to clearly bring out the weak extended structures.

2.2.1 Calibration

The data were calibrated in the CASA software suite (McMullin et al. 2007)3. The compact quasar 3C 286 served as the primary flux density and bandpass calibrator with its flux density scale given by Perley & Butler (2017), and J1851+005 was used as the complex gain calibrator. A routine direction-independent calibration strategy was followed, except for the gain calibrator field. As J1851+005 lies in the Galactic plane with bright and extended sources in its vicinity (see Fig. 3), we first generated a model of this region by self-calibrating this field and then obtained new gain calibration solutions from the model, which were later applied to the target fields. Data corrupted due to RFI were “flagged” (i.e., removed from further use) using automated techniques such as the tfcrop and the rflag modes of the CASA task flagdata. We used some strategies and modules of the VLA scripted pipeline4 to make our calibration scripts fast and reliable. Due to the bright extended structures in the Galactic plane that we intend to detect and study, a direction-dependent calibration is complicated in this case (e.g., Albert et al. 2020). Even with full synthesis observations, imaging extended structures with the GMRT using direction-dependent gains is challenging (e.g., Wykes et al. 2014) and we have not found any previous studies that have successfully imaged large and bright extended sources with direction-dependent calibration. Hence, we decided to proceed with only the standard direction-independent calibration. A flowchart of the calibration process is shown in Fig. 4 and a full description is given in Appendix A.

Preliminary images of the calibrators were produced with the task tclean in order to confirm that the calibration procedure is successful. The task statwt was applied to the target data to weigh down outliers based on local variances. This was helpful to reduce the imaging artefacts caused by residual RFI which may not have been flagged during the previous steps. The corrected data of the target area were then considered ready for further processing.

At the end of this calibration procedure, about 50% of the data in band-3 were flagged, which is typical of uGMRT observations taken at these frequencies due to strong RFI. In band-4, ~15–45%o of the data were flagged, depending on the pointing and the date of the observation (see Table 1). RFI at these frequencies is relatively less of an issue, and a large fraction of the flagging was caused by non-working antennas. The band-4 data from May 27 were the most affected; five antennas were not in working condition during this observation in addition to the intermittent RFI typical at low frequencies, leading to a total of ~43% data loss (see Fig. 5). We also found that several fields had a number of central square antennas flagged for a portion of the scans during the calibration due to RFI, especially in band-4. This meant that the large-scale structure was not well recovered for these fields and the largest angular scales to which the data are sensitive are typically about ~6′ and ~4′ in band-3 and band-4, respectively, significantly smaller than the expected values of ~25′ and ~15′.

thumbnail Fig. 4

Calibration scheme that we followed to reduce the data.

thumbnail Fig. 5

uv-coverage plots of typical fields observed on 2019 August 22 (left, 3SRC105), May 27 (middle, 4SRC001), and May 26 (right, 4SRC104). The asymmetry caused by non-functional antennas (marked in yellow) is clearly visible in the fields observed on August 22 and May 27, which leads to an elongation of the synthesized beam.

2.2.2 Polarization

The band-4 observations (550–750 MHz) were done in full-polarization mode, which allows for images of the Stokes Q, U, and V data to be obtained in addition to Stokes I. Linearly polarized emission is expected from synchrotron emitters such as SNRs, whereas emission from H II regions due to thermal Bremmstrahlung is unpolarized. Hence, the polarization information plays a key role in separating thermal and nonthermal emission and, thus, in identifying new SNRs (e.g., Dokara et al. 2021). The polarization calibration involves obtaining delays of cross-hand polarizations (RL, LR), characterizing instrumental polarization by determining the frequency-dependent leakage terms (colloquially known as the “D-terms”) and calibrating the absolute polarization position angle using a well-known polarized calibrator. The D-terms can be estimated either by observing an unpolarized source for a short interval or by tracking a polarized calibrator with a good coverage of parallactic angle. In order to minimize the observation time, we chose to observe 3C 48, which is unpolarized at these frequencies (Perley & Butler 2013; Farnes et al. 2014), allowing us to determine the leakage terms. For calibrating the polarization position angle, we used 3C 286 following Mohan et al. (2019). However, we obtained a polarization fraction of ~3.5% for 3C 286 at 650 MHz, which agrees with the value of 2.7% at 610 MHz by Farnes (2012) but is in tension with the value of 7.6% at 607 MHz given by Mohan et al. (2019). We note that the polarization properties of 3C 286 are not well known below 1 GHz and it is observed to strongly depolarize (to lower than 1 %) at frequencies below 500 MHz (Perley & Greisen 2019). We imaged the Stokes Q and U data for a field with the W44 supernova remnant, where we expected a polarization signal of at least a few percent, but we were unable to recover any.

Sources with an absolute value of the rotation measure (RM)5 higher than ~7 rad m−2 would appear depolarized in our images, since their polarization angle rotates multiple times over the 200 MHz bandwidth. Galactic sources can have RMs much larger than ~7 rad m−2 and RMs may even be in the range of hundreds of rad m−2 (e.g., Gaensler et al. 2000). To be able to detect the linear polarization from sources with RMs as high as ~100 rad m−2, we split the data into 20 frequency bins, imaged each bin, and made the linearly polarized intensity image of each bin separately; finally, we combined them by averaging to form a single output image. However, we were unable to obtain any useful Stokes Q and U images, which is due to the poor signal-to-noise ratio (S/N) in each bin. This is not unexpected, since the total time spent on each pointing is less than 15 min and the linearly polarized intensity is expected to be only a fraction of the total intensity, meaning the image reconstruction in each frequency bin will be quite poor. Hence, we imaged only the Stokes I data.

2.2.3 Imaging

Generally, imaging is the most intensive part of the reduction of radio interferometric data, both in terms of CPU time and manual effort. The task tclean of CASA was employed for imaging. Several deconvolution algorithms based on the Cotton-Schwab CLEAN are available in this task. It was used in parallel mode on four cores in order to speed up the imaging process. The complex structures in the target region, the wide bandwidth of the receivers, and the wide field of view at low frequencies pose significant challenges in imaging. We briefly explain our imaging strategies below:

w-projection. To account for the non-coplanar baseline effect, we used w-projection (Cornwell et al. 2008). The number of planes was set to be automatically calculated using the parameter wprojplanes=−1.

Outliers. We imaged a region that extends at least out to the first null of the primary beam (~2.3° and ~1.4° at 400 MHz and 650 MHz, respectively) and, in some cases, we used outlier fields as well. This was done in order to image very bright sources that are located well outside the primary beam. Such sources, if not CLEANed, can produce noticeable features and significantly increase the root mean square (rms) background noise.

Wide bandwidth. In order to account for the frequency-dependent variations across the 200 MHz bandwidth, we used the Multi-Term Multi-Frequency synthesis (MT-MFS) algorithm (specmode=’mfs’ and deconvolver=’mtmfs’; Rau et al. 2009; Rau & Cornwell 2011) with two Taylor coefficients in the spectral model (nterms=2).

Multiple scales. To model emission structures of various sizes, we used the multi-scale CLEAN (Cornwell 2008) along with the aforementioned MT-MFS algorithm. The scales were chosen such that they roughly correspond to the sizes of structures expected to be detected in our observations of a given region.

Thresholds. The MT-MFS algorithm in the task tclean is prone to divergence if the stopping criteria and the scales are not carefully set, especially if the noise threshold to be reached is set to a value that is lower than the achievable noise levels. Hence, we used a dynamic threshold based on the local image statistics, fixed the maximum number of minor cycle iterations to 2000 (with the parameter cycleniter), and lowered the loop gain from the default of 0.1 to 0.06. This came at the cost of a larger computational footprint and more processing time, due to major cycles being visited more often. However, we found that the algorithm does not diverge with these parameters, and the images were also more accurate due to better residual image construction.

Weighting. The Briggs weighting scheme (Briggs 1995) with a robust parameter of 0.5 was used, which is a compromise between the best achievable resolution and the best sensitivity to large-scale structures.

Pixel size. The pixel size was chosen such that there are at least three pixels across the minor axis of the expected synthesized beam. This allows a good Gaussian fit for the central lobe of the point spread function.

Masking. It is impractical to manually identify real emission for CLEAN masks since we covered a large region. Therefore, we made use of the auto-masking feature of tclean (usemask=’auto-multithresh’; Kepley et al. 2020). This also helps to keep the imaging process reproducible and minimize human errors and bias.

Self-calibration. By modelling the phases from an initial shallow CLEAN, we self-calibrated the target fields and repeated this process twice to generate a better model and a better final image. Although the improvement in dynamic range was quite low (~15%, which was expected due to the poor uv-coverage), we found that it improved the image fidelity by reducing the sidelobe artefacts caused by imperfect deconvolution.

Wide-band primary beam correction. The attenuation caused due to the primary beam response was corrected for using the contributed CASA task wbpbgmrt6.

Uniform resolution. Due to the time-variable RFI and the non-working antennas, the resolution achieved greatly depended on the field and the time of an observation. For the sake of uniformity and to allow comparison of the images of both bands, we restricted the uv-range to 10 to make images with a common circular beam. The “full-resolution” images were made separately with no restriction on the uv-range.

Mosaics. We imaged each field individually first and combined them later to form a mosaic using the tool linearmosaic of CASA, where the pixels in the overlapping regions are weighted by using the primary beam response as the sensitivity. We note that while wide-band mosaic imaging is available in CASA7, the vastly varying beam sizes across the target fields observed during different nights (see Sect. 3) make it unsuitable for imaging our data.

Coordinate transformation. Finally, the mosaic was regridded from equatorial to Galactic coordinates using the software Montage (Berriman et al. 2003).

We found that the synthesized beams obtained after imaging depended strongly on the day of the observation (due to the varying number of working antennas and the time-variable RFI) and also on the software being used. AIPS, and the versions of CASA prior to v5.8, have a PSF-fitting algorithm that works reliably only if the range of over-sampling factor (number of pixels across an axis of the expected synthesized beam) is from three to five. Because of the non-functional outer antennas during the May 27 and August 22 observations (see Fig. 5), employing an over-sampling factor of 3–5 was not simultaneously possible on both the major and minor axes of the synthesized beam when we produced images using visibilities from all functional baselines. This led to an incorrect beam size estimation if AIPS or an older version of CASA was used for imaging (see Appendix B for more details). Hence, we used a version of CASA, v5.8, in which the PSF fitting algorithm was updated8 so that it is able to deliver correct results even in extreme cases such as ours. Unfortunately, the non-working antennas meant that the final common beam size for both bands we were able to obtain is 25″, which is about three times larger than the best possible resolution of about 8″. In Fig. 6, we show an image resulting from a pointing observed in band-4 on a night with all working antennas, with no primary beam correction and at its native resolution. An rms noise of about 0.3 mJy beam−1 and a beam-size of 12″ × 7″ were achieved, which are close to the theoretical values for this observation.

Background noise maps corresponding to each band were created using the software SExtractor (Bertin & Arnouts 1996). A spectral index map was created using the images from the two bands, after reprojecting them onto the same grid. This map was obtained on a pixel-by-pixel basis, and the calculation was done only if the S/N of a pixel in both images was greater than three. We note that although the multi-frequency synthesis (MFS) algorithm in CASA (Rau et al. 2009) is capable of producing spectral index maps directly during imaging, we find it to deliver un-physical spectral indices for our observations. We routinely obtained spectral index values of the order of tens even after accounting for the primary beam attenuation. This was true not just for faint or extended emission, but also for point-like sources with S/N values of about 50. A recent study by Rashid et al. (2023) analyzed this issue using CASA simulations and also came to the same conclusion that the spectral index maps produced by the MFS algorithm are not reliable for uGMRT data with a S/N below 100. Hence, for this work, we use only the spectral index map derived using the band-3 and the band-4 mosaics together.

2.3 Missing flux density

All interferometers spatially filter out a fraction of the flux density of extended sources. The angular scale at which such filtering starts to affect the flux densities depends on the uv-coverage and the structure of the extended source, but an upper limit of the largest angular scale recovered in a full-synthesis observation (θLAS) can be estimated using θLASλbmin,${\theta _{{\rm{LAS}}}} \mathbin{\lower.3ex\hbox{$\buildrel<\over {\smash{\scriptstyle\sim}\vphantom{_x}}$}} {\lambda \over {{b_{\min }}}},$

where λ is the wavelength and bmin is the length of the shortest available baseline (~100 m for the uGMRT). We found that the values of θLAS in our observations are about 6′ and 4′ in band 3 and band 4, respectively, and in some regions it is even smaller. These are significantly lower than the estimates for two reasons. One, we observed in the snapshot-mode, implying that the uv-coverage is quite poor compared to a full-synthesis observation. Two, the RFI situation at the observatory is such that it affects the central square antennas much more than the outer antennas. This further degrades the quality of image reconstruction of extended structures. In order to mitigate this issue of missing flux density, zero-spacing information can be added from a dataset at the same frequency that recovers the extended emission, such as the data from a single-dish telescope. This can be done using a method called “feathering”, or by joint deconvolution (e.g., Vogel et al. 1984; Koda et al. 2011; Rau et al. 2019; Dokara et al. 2023). In any case, the diameter of the single-dish telescope must be greater than the shortest baseline length for the combination to work well, since, only then an overlap in the uv-coverage is guaranteed. Only two steerable single-dish telescopes exist that are suitable for such a zero-spacing correction of the GMRT data: the 100 m Effelsberg telescope and the 100 m Robert C. Byrd Green Bank Telescope. While there are receivers on both the telescopes that can complement the frequency coverage in our study, the RFI from man-made electronics at both the observatories renders a large fraction of the bandwidth useless. It currently remains unfeasible to perform a large-scale survey using the available bandwidth on these receivers to reach the sensitivity required to be added to our uGMRT data.

thumbnail Fig. 6

Example field observed in band-4 during a night with all antennas functioning, uncorrected for primary beam attenuation, in the native equatorial coordinates. The white circle marks the region used for mosaicking. The beam size is 12″ × 7″, and the rms noise is about 0.3 mJy beam−1. The colormap and the limits are chosen such that the low-level noise features and the negative sidelobes are seen clearly.

3 Results

The mosaics of band-3 and band-4 at a resolution of 25″, along with the spectral index map, are shown in Fig. 7. All the uGMRT images in the figure are dynamic range limited, with rms noise levels of about ~5 mJy beam−1 in band 3 and ~2 mJy beam−1 in band 4. The distributions of the brightness values of the pixels in the mosaics and the cumulative distributions of the values of pixels in their respective background noise maps are shown in Fig. 8. We notice sidelobe artefacts near a few bright point sources in both bands, which are likely due to the lack of direction-dependent calibration and imperfect deconvolution in conjunction with the sparse uv-coverage. Bright sources that are not close to the center of the pointing corrupt the image, since the gains are truly applicable only at the phase center of the field. In addition, the regions close to the SNR W44 (G34.7−0.4) have a larger rms noise because the SNR is quite bright and extended with a complicated structure. In the band-3 mosaic, in the region enclosed between the Galactic coordinates l = 34°–35° and b = 0°–0.5°, there appears to be a very faint, large-scale structure bending away from the Galactic mid-plane. This feature, however, does not have a counterpart in the 70–230 MHz GLEAM images (Hurley-Walker et al. 2019) or the THOR+VGPS data (Beuther et al. 2016) and its surface brightness is close to the local noise level as well. So, this feature is likely to be an artifact of the imaging process.

3.1 Comparison with other surveys

To date, there is no survey at high angular resolution that covers the region surveyed in our pilot study at frequencies of 300–750 MHz. Hence, we used data from surveys conducted at other frequencies to compare the positions and flux densities of the sources detected in our survey.

The 150 MHz TIFR-GMRT sky survey (TGSS; Intema et al. 2017) covers the region we observed at the same resolution (25″) and a comparable sensitivity (~10 mJy beam−1 in the Galactic plane). Since the inner baselines were severely down-weighted during the data reduction, only the brightest edges of extended sources were detected in their images. The catalog of this survey lists 50 sources in the region we observed. We detect all the sources except one, named J185141.6+003739 in the TGSS catalog. Its S/N is about nine, but this source has not been detected in any radio survey other than the TGSS, according to the SIM-BAD database. There are two bright sources very close to this object: J185146.7+003531, which is, in fact, our gain calibrator, and SNR G33.6+0.1, which is an extended source (Fig. 3) and thus not well sampled in the TGSS images. Hence, the sole TGSS source that we do not detect in our data is likely to be a “ghost” (Grobler et al. 2014)9. Of the remaining 49 TGSS sources, 17 are clearly either a part of a supernova remnant or are double-lobed sources that are almost certainly distant radio galaxies, and the other 32 are compact sources.

The NRAO-VLA Sky Survey (NVSS; Condon et al. 1998) is a 1.4 GHz survey of the whole sky north of δ = −40°, at a resolution of 45″. While the NVSS catalog contains over 1300 sources in the area we covered, a visual inspection reveals that most of these sources are obviously false. This was noted by earlier studies as well (Bihr et al. 2016) and this is not unexpected. The techniques used for data reduction and cataloging of the NVSS data are best suited for imaging and detecting compact extragalactic objects and do not perform well in very crowded and complex fields such as ours. However, of the 32 compact sources detected in TGSS, 28 were detected in the NVSS as well, using which we made position and flux density comparisons to validate our mosaics. For these sources, we evaluated their peak intensities and integrated flux densities, along with their positions, by performing two-dimensional Gaussian fits on our images. We compared these values to the values from the TGSS and the NVSS source catalogs, and the results are shown in Fig. 9.

While the position accuracies of both TGSS and NVSS are reported to be better than 2″ (Condon et al. 1998; Intema et al. 2017), we find that there are significant offsets between the peak positions of a few sources in this field (Fig. 9, left panel). The complex background in the region in question moves the peak of a source along the local noise gradient, which makes many positions listed in these source catalogs less reliable compared to the regions outside the Galactic plane. In addition to background contamination, the emission at different frequencies may be dominated by different mechanisms, which may result in shifting peak positions, or they may be entirely different sources altogether. While the ionosphere is expected to contribute to propagation delay at these frequencies, each individual antenna of the GMRT sees a roughly constant ionosphere above ~400 MHz. Hence, antenna-based phase-only self-calibration is usually sufficient (Lonsdale 2005; Intema et al. 2009) to deal with the resulting phase errors. Since we have not performed direction-dependent calibration, calibration issues may also play a role in shifting the peak positions. However, we note that most of the offsets are within 1–2 pixels, and all are smaller than the beam size.

Assuming that a single power-law index holds between 150 MHz and 1.4 GHz and that the sources are not variable, a spectral index catalog was prepared by de Gasperin et al. (2018) using the data from the TGSS and the NVSS surveys. It has already been shown to be useful to discriminate between thermal and nonthermal emission from and H II regions and SNRs (Dokara et al. 2018), in addition to pinpointing extragalactic sources. Here we use this spectral index catalog to derive an expected flux density of the sources at 400 MHz and 650 MHz. We compared these values with the integrated flux densities from the Gaussian fits to band-3 and band-4 sources, and the results are shown in the middle and right panels of Fig. 9.

Although there appears to be a general agreement between the values observed by us and the derived estimations, the number of bright sources that fall below the 1:1 line is larger than the number above the 1:1 line; namely, the number of bright sources with observed flux density lower than the expected flux density is visibly high. While this may be caused by physical reasons such as the assumptions of a single power-law spectral index being invalid and/or variability, it may also point to a systematic error within one of the TGSS, NVSS, and our band-3 and band-4 images. The NVSS flux densities are based on the scale of Baars et al. (1977), while we used the recent Perley & Butler (2017) scale. On the other hand, the TGSS is scaled using the low frequency scale of Scaife & Heald (2012), which is not based on an absolute and independent standard, but it is still consistent to within 2% of the other two scales. Hurley-Walker (2017) finds systematic position-based amplitude offsets on the angular scale of degrees, which can reach as high as 40% in some regions (see their Fig. 3). It had been noted that the higher system temperatures (Tsys) in Galactic plane observations may lead to the receivers being pushed into a non-linear regime10. Since our observations were done with the upgraded receivers, this is unlikely to affect the results, and we find that the flux densities of the primary and secondary calibrators are a good match with the literature values.

thumbnail Fig. 7

Radio images of the W43/W44 region. From top to bottom: band-3 and band-4 mosaics from this project convolved to the common beam size of 25″, their resulting spectral index map, the 1.4 GHz THOR image of the same area combined with the VGPS data (Beuther et al. 2016; Stil et al. 2006), which is also at 25″ resolution, and the 5.8 GHz GLOSTAR data (VLA D-configuration image combined with the Effelsberg image), which is at a resolution of 18″.

thumbnail Fig. 8

Cumulative distribution functions of the noise maps (left), the histograms of values of the pixels in the continuum mosaics of band-3 (middle), and band-4 (right). Gaussian least-squares fits are performed on the histograms, which gave standard deviations as ~4.9 mJy beam−1 for band-3 and ~1.8 mJy beam−1 for band-4.

thumbnail Fig. 9

Position offsets of the TGSS sources seen in our survey and the NVSS, with obvious extended sources excluded (left). Observed and expected flux density comparison is shown for sources in band-3 (middle) and band-4 (right), where the expected flux density is obtained using spectral indices from de Gasperin et al. (2018).

thumbnail Fig. 10

Probability density functions obtained using kernel density estimations of the values of spectral index pixels belonging to SNRs and H II regions.

3.2 Nonthermal emission

H II regions are routinely misidentified as SNRs in the Galactic plane due to their similar radio morphology (see Dokara et al. 2021, for instance); however, the spectral index is an excellent discriminant between the two. H II regions typically have a “flat” or “rising” spectrum at radio wavelengths (−0.1 < α < 2) depending on their optical depth (Wilson et al. 2013). In contrast, shell-type SNRs have α ≲ −0.5 and Crab Nebula-like SNRs (pulsar wind nebulae, PWNe) or composite SNRs have α ~ −0.2 (Bietenholz et al. 1997; Dubner & Giacani 2015). In Fig. 10, we show the probability distribution functions of the pixels in the spectral index map within the angular extents of known SNRs and H II regions. The list of H II regions is taken from the WISE catalog of H II regions (Anderson et al. 2014) and that of the SNRs from the latest version of Green’s catalog (Green 2019). Accurate spectral index determination based on our uGMRT data is hampered by the following problems. One, the central frequencies of the two bands covered in our survey are not very different (400 MHz and 650 MHz). Two, the data in the two bands have a non-trivial overlap of uv-coverage. Three, a significant fraction of the flux density remains undetected due to the missing short-spacing information (see also Sect. 2.3). Nevertheless, we see the spectral indices of SNRs to be significantly lower than those of H II regions, which is what we expected.

The measured flux densities of SNRs detected in our images are presented in Table 2 compared with expected values extrapolated from the 1.4 GHz flux densities from the THOR+VGPS data (Beuther et al. 2016) and the broad-band spectral indices derived by Dokara et al. (2023). Large-scale emission is filtered out by the uGMRT, while the THOR+VGPS data have zero-spacing information added. Consequently, the expected flux densities are anticipated to be larger than the observed values, which is indeed found to be the case for all the remnants except for G31.9+0.0. In band 4, its measured and expected flux densities are equal within uncertainties, but in band 3, they are significantly different. This is perhaps because the spectral indices were calculated by Dokara et al. (2023) using a single power-law spectrum. While G31.9+0.0 is known to show a turnover in its spectrum below 100 MHz (Brogan et al. 2005), there is a wide variation in the reported flux densities below 1 GHz as well (see Kassim 1989, 1992; Brogan et al. 2005; Sun et al. 2011).

We searched for the recently identified SNR candidates in the THOR and GLOSTAR surveys (Anderson et al. 2017; Dokara et al. 2021) in our images, but at a 3σ-level, we cannot positively identify any of the 15 candidates that are located in the covered region. Considering the fact that SNRs generally do not show spectral indices steeper than ~−0.8 and that the largest angular scale detected in our data is less than 10′, this inability to identify their counterparts in our images is not unexpected. The sensitivity to extended emission must be generally better than 1–2 mJy beam−1 in our data in order to detect these faint SNR candidates. For the PWN-like GLOSTAR SNR candidate G031.256-0.041 and THOR SNR candidate G32.22-0.21, we are able to derive a lower limit of the spectral index (using the average surface brightness) of ~−0.1 and ~−0.5, respectively. We also searched for shell-shaped objects that may have been undetected in the THOR and the GLOSTAR surveys because of synchrotron ageing at higher frequencies, but we did not find any such objects. Longer integration time is required to constrain the spectral index for more SNR candidates and the identification of new candidates.

4 Prospects for a future work

Recently, based on the new and the revised distances to several SNRs, Ranasinghe & Leahy (2022) provided evidence that the number of SNRs in the Milky Way must be at least over 2000 and may be as high as 5000 – and not just 1000, as proposed by Li et al. (1991). These numbers are glaringly inconsistent with that of the only 400 SNRs detected so far (Green 2019; Ferrand & Safi-Harb 2012). Since a large fraction of the SNRs are expected to be only a few arcminutes in extent (Dokara et al. 2023) and since SNRs are brighter at lower frequencies, significant progress can be made toward solving the problem of “missing” SNRs via a high-resolution, sub-GHz, sensitive survey of the Milky Way. Such a survey will also be useful to understand the spectral turnover caused by synchrotron self-absorption observed in many SNRs (e.g., Sun et al. 2011; Dokara et al. 2023), as the turnover frequency is generally lower than 1 GHz. In addition, it is possible to also study the nonthermal emission that arises from star-forming regions (Vig et al. 2018; Veena et al. 2019). The wide-band receivers on the uGMRT appear promising for a large-scale survey of the Milky Way. The promise of detecting SNRs in such low-frequency surveys is illustrated by the fact that Brogan et al. (2006) detected 35 new SNRs in a 44 square degree area of the Galactic plane at 75 MHz (90 cm) with the classic VLA.

We performed this pilot study as a precursor to a larger survey of the Milky Way, the Metrewave Galactic Plane with the uGMRT (MeGaPluG) Survey, which we wish to undertake in the coming years. We have developed automated calibration and imaging pipelines to reduce the large amount of data in an efficient manner. The primary goal was to study the feasibility of imaging extended structures with the wide-band uGMRT at a high resolution, using only snapshot observations. Despite losing sensitivity and resolution due to the problem of non-working antennas, as is evident from Fig. 7, we recovered the extended emission reasonably. This makes this a feasibly sensitive study for imaging extended structures at low frequencies (< 1 GHz) and high resolution (25″). Based on the experiences from this study, for the upcoming larger survey, we have designed an optimum strategy for overcoming the shortcomings of the pilot study:

  • Undertake multiple scans of each field instead of just two scans as was done for this work. The uv-coverage can be improved if we use four scans of 4–5 min spread over an observing session. This also helps in minimizing the uv-coverage loss due to non-functional antennas during the observation;

  • The pointing configurations chosen (see Fig. 2) were not optimal given the size of primary beams. A denser grid of pointings, as done in the GLOSTAR survey (Medina et al. 2019; Brunthaler et al. 2021), will offer a better sensitivity, although this will increase the observation time;

  • In targeted observations, we can afford to observe the gain calibrator only once in an hour or so, since self-calibration is almost always guaranteed to enhance the dynamic range. However, in our snapshot survey mode, self-calibration results in only a marginal improvement. Hence, it may be wise to observe the gain calibrator more often (e.g., once every half an hour);

  • While polarization data are incredibly useful to study nonthermal emission, it is not possible to obtain any meaningful Stokes Q and U images of the target fields in our snapshot observations. The polarization angles of sources with |RM| > 7 rad m−2 rotate over the 200 MHz bandwidth, causing depolarization. While this limitation may be overcome by dividing the bandwidth into multiple frequency bins, the resulting deterioration of the S/N makes the images resulting from the deconvolution process unusable. Also, given that a large number of sources are located away from the center of the pointing, it is necessary to know the position dependence of the instrumental polarization effects, which has not yet been investigated for the uGMRT. Hence, it is prudent to image and study only the Stokes I data under these circumstances. We note that since the cross-hand polarizations are helpful in eliminating the RFI, they may still be recorded.

Table 2

Supernova remnant flux densities in the MeGaPluG pilot.

Acknowledgements

We thank the anonymous referee for their comments and suggestions. R.D. is a member of the International Max Planck Research School (IMPRS) for Astronomy and Astrophysics at the Universities of Bonn and Cologne. N.R. acknowledges support from Max-Planck-Gesellschaft through Max-Planck India partner group grant. This work partially uses images from the GLOSTAR survey supported by the MPIfR, Bonn (https://glostar.mpifr-bonn.mpg.de). The GLOSTAR data products are the results of observations conducted with the Karl G. Jansky Very Large Array operated by the National Radio Astronomy Observatory and the Effelsberg 100 meter radio telescope operated by the MPIfR. This research has made use of NASA’s Astrophysics Data System and the SIMBAD database. We have used the softwares Astropy (Astropy Collaboration 2013), APLpy (Robitaille & Bressert 2012), Montage (Berriman et al. 2003), and DS9 (Joye & Mandel 2003) at various stages of this research.

Appendix A Calibration scheme

  1. FITS to CASA Measurement Set (MS) conversion: The data for each observing session are provided on the NAPS website in a standard FITS format, which is converted to a MS using the CASA task importgmrt.

  2. Hanning smoothing: In order to prevent Gibbs ringing, the task hanningsmooth is used to smooth the data with a Hanning window.

  3. Initial flagging: “Flagging” is the process of masking the data that cannot be recovered (generally due to poor instrument performance or RFI). The task flagdata is used for this purpose. We flag dead antennas, shadowed antennas, and edge channels with low amplitude gains. Then, using the tfcrop mode of the task flagdata, we performed a round of automated RFI excision. The tfcrop mode uses an algorithm that identifies and flags outliers on the 2D time-frequency plane based on local statistics. It can operate on un-calibrated data as it can account for the bandpass shape, and it is especially helpful in flagging short duration RFI and time-persistent narrow-band RFI.

  4. Set models: The task setjy is used to set the model of the primary flux calibrator, 3C 286, based on the flux density scale of Perley & Butler (2017). The secondary calibrator, J1851+005, is initially assumed to be a point source with an amplitude 1 Jy before a model of the field is produced.

  5. Delay and bandpass solutions: Phase-corrected data of the primary flux calibrator are used to derive the corrections due to instrumental delays and the bandpass shape.

  6. Apply, flag, and repeat twice: The delay and bandpass solutions are applied to the calibrator fields and the corrected data are flagged with tfcrop. New solutions are then derived keeping these flags. This was done twice in order to ensure that the calibration tables are obtained only from reliable data.

  7. Prepare for gain calibration: The final calibration tables for delay and bandpass are applied on the primary and secondary gain calibrators. Their corrected data are split to a new MS with appropriate channel-averaging. This step significantly speeds up11 the upcoming temporal gain calibration steps.

  8. Gain calibration: Using the data of primary and secondary calibrator fields, time-dependent solutions are derived for gain amplitude and phase variations caused by instrumental response and ionospheric effects. The gain calibration tables are applied to the secondary calibrator field, and outliers in the corrected data column are flagged using the tfcrop mode of the task flagdata. New gain solutions are obtained keeping the latest flags. The amplitude corrections to the secondary calibrator are then scaled using the gains from the primary calibrator, with the task fluxscale.

  9. A better gain calibration by modeling the secondary calibrator: If the secondary calibrator has been modeled as a point source, then a correct model needs to be made by self-calibration12. The CLEAN components of the final self-calibrated image are used as the model of J1851+005. This step is done in order to account for the nearby extended sources in the model of the gain calibrator field. New gain solutions are obtained using this model of the gain calibrator field.

  10. Apply calibration tables: The latest set of tables are then applied to the target fields. The time-dependent gain solutions are linearly interpolated between adjacent scans and the data that could not be calibrated and those with low S/N values are flagged.

  11. Flag RFI on targets: We perform a step of automated flagging on the target fields with the task flagdata in the modes tfcrop and rflag.

  12. Split target fields: The calibrated data of target fields are split into separate MSs with the task split and averaged in bandwidth in order to make the imaging process faster and convenient. The data are averaged only to some extent such that no bandwidth smearing occurs (channel width must be ≲0.7 MHz in the observed frequency regime).

Appendix B Beam size estimation

The point spread function (PSF, dirty beam, or Bdirty) of an interferometric observation depends on the array configuration and the sky position of the target being observed. The raw image (dirty image, or Iduty) made by an interferometer before the CLEAN process is the sky brightness distribution (Isky) convolved with the PSF. The final reconstructed image Ireconstr is made by adding the residuals (Ires) to the convolution of a Gaussian CLEAN beam Bclean with the model of the sky image obtained by deconvolving the PSF from the dirty image, where the Gaussian CLEAN beam is estimated by fitting a Gaussian to the central lobe of the dirty beam. In effect, the dirty beam is replaced by the CLEAN beam: Idirty=IskyBdirty=ImodelBdirty+IresIreconstr=ImodelBdirty+Ires,$\matrix{{{I_{{\rm{dirty}}}}} \hfill & = \hfill & {{I_{{\rm{sky}}}}\oast {B_{{\rm{dirty}}}}} \hfill \cr {} \hfill & = \hfill & {{I_{{\rm{model}}}}\oast{B_{{\rm{dirty}}}} + {I_{{\rm{res}}}}} \hfill \cr {{I_{{\rm{reconstr}}}}} \hfill & = \hfill & {{I_{{\rm{model}}}}\oast{B_{{\rm{dirty}}}} + {I_{{\rm{res}}}},} \hfill \cr}$(B.1)

In CASA, all the above is taken care of by a single task tclean. In the version 5.8, a major upgrade was made to the PSF fitting algorithm13, which changes the results of the Gaussian CLEAN beam fitting. The images from the observations of May 27, during which one arm of antennas was completely flagged are drastically affected by this change. The older CASA versions before 5.8 (ie ≤5.7) gave the beam size for the May 27 observations as ~ 9″ × 8″, whereas we obtain a size of ~ 20″ × 5″ when using v5.8. The result from v5.7 is incorrect because the southern arm was completely flagged during that observation and the PSF must be highly asymmetric for snapshot observations such as ours.

To make sure that this discrepancy was not a result of improper calibration, we imaged a simulated GMRT dataset made with CASA simulations without any corrupting gains. We flagged the antennas on one arm and imaged the data with the same parameters in AIPS, CASA v5.7 and CASA v5.8. The results were similar to the beam sizes obtained in the real data, and hence we confirm that the PSF fitting algorithm used in AIPS and older versions of CASA delivers incorrect results when applied to our observations. All the imaging took place in CASA v5.8 in order to ensure that the beam sizes were correctly estimated. Incorrect beam size estimation leads to unnatural flux densities of imperfectly deconvolved sources.

References

  1. Albert, J. G., van Weeren, R. J., Intema, H. T., & Röttgering, H. J. A. 2020, A & A, 635, A147 [CrossRef] [EDP Sciences] [Google Scholar]
  2. Anderson, L. D., Bania, T. M., Balser, D. S., et al. 2014, ApJS, 212, 1 [Google Scholar]
  3. Anderson, L. D., Wang, Y., Bihr, S., et al. 2017, A & A, 605, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Astropy Collaboration (Robitaille, T. P., et al.) 2013, A & A, 558, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Baars, J. W. M., Genzel, R., Pauliny-Toth, I. I. K., & Witzel, A. 1977, A & A, 500, 135 [Google Scholar]
  6. Becker, R. H., White, R. L., Helfand, D. J., & Zoonematkermani, S. 1994, ApJS, 91, 347 [Google Scholar]
  7. Berriman, G. B., Good, J. C., Curkendall, D. W., et al. 2003, ASP Conf. Ser., 295, 343 [Google Scholar]
  8. Bertin, E., & Arnouts, S. 1996, A & AS, 117, 393 [NASA ADS] [Google Scholar]
  9. Beuther, H., Bihr, S., Rugel, M., et al. 2016, A & A, 595, A32 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Beuther, H., Walsh, A., Wang, Y., et al. 2019, A & A, 628, A90 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Bietenholz, M. F., Kassim, N., Frail, D. A., et al. 1997, ApJ, 490, 291 [NASA ADS] [CrossRef] [Google Scholar]
  12. Bihr, S., Johnston, K. G., Beuther, H., et al. 2016, A & A, 588, A97 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Briggs, D. S. 1995, AAS Meet. Abstr., 187, 112.02 [NASA ADS] [Google Scholar]
  14. Brogan, C. L., Lazio, T. J., Kassim, N. E., & Dyer, K. K. 2005, AJ, 130, 148 [NASA ADS] [CrossRef] [Google Scholar]
  15. Brogan, C. L., Gelfand, J. D., Gaensler, B. M., Kassim, N. E., & Lazio, T. J. W. 2006, ApJ, 639, L25 [Google Scholar]
  16. Brunthaler, A., Menten, K. M., Dzib, S. A., et al. 2021, A & A, 651, A85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Chakraborty, A., Roy, N., Wang, Y., et al. 2020, MNRAS, 492, 2236 [Google Scholar]
  18. Condon, J. J., Cotton, W. D., Greisen, E. W., et al. 1998, AJ, 115, 1693 [Google Scholar]
  19. Cornwell, T. J. 2008, IEEE J. Sel. Top. Signal Process., 2, 793 [Google Scholar]
  20. Cornwell, T. J., Golap, K., & Bhatnagar, S. 2008, IEEE J. Sel. Top. Signal Process., 2, 647 [Google Scholar]
  21. de Gasperin, F., Intema, H. T., & Frail, D. A. 2018, MNRAS, 474, 5008 [Google Scholar]
  22. de Gasperin, F., Williams, W. L., Best, P., et al. 2021, A & A, 648, A104 [CrossRef] [EDP Sciences] [Google Scholar]
  23. Dokara, R., Roy, N., Beuther, H., et al. 2018, ApJ, 866, 61 [Google Scholar]
  24. Dokara, R., Brunthaler, A., Menten, K. M., et al. 2021, A & A, 651, A86 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Dokara, R., Gong, Y., Reich, W., et al. 2023, A & A, 671, A145 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Dubner, G., & Giacani, E. 2015, A & ARv, 23, 3 [NASA ADS] [Google Scholar]
  27. Farnes, J. S. 2012, Ph.D. Thesis, University of Cambridge, UK [Google Scholar]
  28. Farnes, J. S., Green, D. A., & Kantharia, N. G. 2014, MNRAS, 437, 3236 [CrossRef] [Google Scholar]
  29. Ferrand, G., & Safi-Harb, S. 2012, Adv. Space Res., 49, 1313 [Google Scholar]
  30. Gaensler, B. M., Dickel, J. R., & Green, A. J. 2000, ApJ, 542, 380 [Google Scholar]
  31. Green, D. A. 2019, JApA, 40, 36 [NASA ADS] [Google Scholar]
  32. Green, A. J., Cram, L. E., Large, M. I., & Ye, T. 1999, ApJS, 122, 207 [NASA ADS] [CrossRef] [Google Scholar]
  33. Grobler, T. L., Nunhokee, C. D., Smirnov, O. M., van Zyl, A. J., & de Bruyn, A. G. 2014, MNRAS, 439, 4030 [CrossRef] [Google Scholar]
  34. Gupta, Y., Ajithkumar, B., Kale, H. S., et al. 2017, Curr. Sci., 113, 707 [NASA ADS] [CrossRef] [Google Scholar]
  35. Hamaker, J. P., Bregman, J. D., & Sault, R. J. 1996, A & AS, 117, 137 [NASA ADS] [Google Scholar]
  36. Heald, G. H., Pizzo, R. F., Orrú, E., et al. 2015, A & A, 582, A123 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Hoare, M. G., Purcell, C. R., Churchwell, E. B., et al. 2012, PASP, 124, 939 [Google Scholar]
  38. Hurley-Walker, N. 2017, ArXiv e-prints, [arXiv:1703.06635] [Google Scholar]
  39. Hurley-Walker, N., Hancock, P. J., Franzen, T. M. O., et al. 2019, PASA, 36, e047 [Google Scholar]
  40. Intema, H. T., van der Tol, S., Cotton, W. D., et al. 2009, A & A, 501, 1185 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Intema, H. T., Jagannathan, P., Mooley, K. P., & Frail, D. A. 2017, A & A, 598, A78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Joye, W. A., & Mandel, E. 2003, ASP Conf. Ser., 295, 489 [Google Scholar]
  43. Kassim, N. E. 1989, ApJS, 71, 799 [NASA ADS] [CrossRef] [Google Scholar]
  44. Kassim, N. E. 1992, AJ, 103, 943 [NASA ADS] [CrossRef] [Google Scholar]
  45. Katz, D., Gomez, A., Haywood, M., Snaith, O., & Di Matteo, P. 2021, A & A, 655, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Kepley, A. A., Tsutsumi, T., Brogan, C. L., et al. 2020, PASP, 132, 024505 [Google Scholar]
  47. Koda, J., Sawada, T., Wright, M. C. H., et al. 2011, ApJS, 193, 19 [NASA ADS] [CrossRef] [Google Scholar]
  48. Kothes, R., Reich, W., Safi-Harb, S., et al. 2020, MNRAS, 496, 723 [Google Scholar]
  49. Li, Z., Wheeler, J. C., Bash, F. N., & Jefferys, W. H. 1991, ApJ, 378, 93 [Google Scholar]
  50. Lonsdale, C. J. 2005, ASP Conf. Ser., 345, 399 [NASA ADS] [Google Scholar]
  51. McMullin, J. P., Waters, B., Schiebel, D., Young, W., & Golap, K. 2007, ASP Conf. Ser., 376, 127 [Google Scholar]
  52. Medina, S. N. X., Urquhart, J. S., Dzib, S. A., et al. 2019, A & A, 627, A175 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Mohan, N., Chellappan Pillai, S. R., Swarup, G., & Oberoi, D. 2019, MNRAS, 487, 4819 [NASA ADS] [CrossRef] [Google Scholar]
  54. Nguyen, H., Rugel, M. R., Menten, K. M., et al. 2021, A & A, 651, A88 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Nguyen, H., Rugel, M. R., Murugeshan, C., et al. 2022, A & A, 666, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Ortiz-León, G. N., Menten, K. M., Brunthaler, A., et al. 2021, A & A, 651, A87 [CrossRef] [EDP Sciences] [Google Scholar]
  57. Perley, R. A., & Butler, B. J. 2013, ApJS, 206, 16 [Google Scholar]
  58. Perley, R. A., & Butler, B. J. 2017, ApJS, 230, 7 [NASA ADS] [CrossRef] [Google Scholar]
  59. Perley, R., & Greisen, E. 2019, EVLA Memos, 207 [Google Scholar]
  60. Ranasinghe, S., & Leahy, D. 2022, ApJ, 940, 63 [NASA ADS] [CrossRef] [Google Scholar]
  61. Rashid, M., Roy, N., Pandian, J. D., et al. 2023, MNRAS, submitted [Google Scholar]
  62. Rau, U., & Cornwell, T. J. 2011, A & A, 532, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Rau, U., Bhatnagar, S., Voronkov, M. A., & Cornwell, T. J. 2009, IEEE Proc., 97, 1472 [NASA ADS] [CrossRef] [Google Scholar]
  64. Rau, U., Naik, N., & Braun, T. 2019, AJ, 158, 3 [NASA ADS] [CrossRef] [Google Scholar]
  65. Robitaille, T., & Bressert, E. 2012, Astrophysics Source Code Library [record ascl:1208.017] [Google Scholar]
  66. Rugel, M. R., Rahner, D., Beuther, H., et al. 2019, A & A, 622, A48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  67. Scaife, A. M. M., & Heald, G. H. 2012, MNRAS, 423, L30 [Google Scholar]
  68. Schuller, F., Menten, K. M., Contreras, Y., et al. 2009, A & A, 504, 415 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  69. Shanahan, R., Lemmer, S. J., Stil, J. M., et al. 2019, ApJ, 887, L7 [NASA ADS] [CrossRef] [Google Scholar]
  70. Shimwell, T. W., Röttgering, H. J. A., Best, P. N., et al. 2017, A & A, 598, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  71. Soler, J. D., Beuther, H., Syed, J., et al. 2020, A & A, 642, A163 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  72. Stil, J. M., Taylor, A. R., Dickey, J. M., et al. 2006, AJ, 132, 1158 [NASA ADS] [CrossRef] [Google Scholar]
  73. Sun, X. H., Reich, P., Reich, W., et al. 2011, A & A, 536, A83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  74. Veena, V. S., Vig, S., Sebastian, B., et al. 2019, MNRAS, 482, 4630 [CrossRef] [Google Scholar]
  75. Vig, S., Veena, V. S., Mandal, S., Tej, A., & Ghosh, S. K. 2018, MNRAS, 474, 3808 [CrossRef] [Google Scholar]
  76. Vogel, S. N., Wright, M. C. H., Plambeck, R. L., & Welch, W. J. 1984, ApJ, 283, 655 [NASA ADS] [CrossRef] [Google Scholar]
  77. Wang, Y., Beuther, H., Rugel, M. R., et al. 2020a, A & A, 634, A83 [CrossRef] [EDP Sciences] [Google Scholar]
  78. Wang, Y., Bihr, S., Beuther, H., et al. 2020b, A & A, 634, A139 [CrossRef] [EDP Sciences] [Google Scholar]
  79. Wilner, D. J., Reid, M. J., & Menten, K. M. 1999, ApJ, 513, 775 [NASA ADS] [CrossRef] [Google Scholar]
  80. Wilson, T. L., Rohlfs, K., & Hüttemeister, S. 2013, Tools of Radio Astronomy (Berlin: Springer Berlin Heidelberg) [Google Scholar]
  81. Wykes, S., Intema, H. T., Hardcastle, M. J., et al. 2014, MNRAS, 442, 2867 [NASA ADS] [CrossRef] [Google Scholar]
  82. Yan, Y. T., Henkel, C., Kobayashi, C., et al. 2023, A & A, 670, A98 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

1

α is the spectral index defined using a power-law relation between the flux density, Sν, and the frequency, ν: Sννα.

5

RM is defined using Δχ = RM · λ2, where Δχ is the change in the electric vector position angle and λ is the wavelength.

9

Ghost sources are calibration artefacts but show a PSF similar to the instrumental PSF. They arise due to imperfect sky models and direction-dependent effects.

11

All GMRT data are taken in spectral line mode (typically 8192 or 16384 channels) in order to isolate the narrow-band RFI which is ubiquitous below 1 GHz at the observatory. Because of the large number of channels, applying delay and bandpass tables ‘on the fly’ (see https://casadocs.readthedocs.io/en/stable/notebooks/synthesis_calibration.html?highlight=fly#Solve-for-Calibration) in order to obtain temporal gain solutions takes a much longer time compared to the method we used.

12

The CLEAN parameters used for the self-calibration of the phase calibrator field are the same as those for target fields, and are discussed in the imaging section.

All Tables

Table 1

MeGaPluG pilot survey observations.

Table 2

Supernova remnant flux densities in the MeGaPluG pilot.

All Figures

thumbnail Fig. 1

Antenna configuration of the uGMRT array. The inset shows the central square (1.6 km × 1.6 km) antennas. North is upwards and east is rightwards.

In the text
thumbnail Fig. 2

Pointing configuration map of our pilot survey. Each large red circle marks the field of view of the pointing used for making a mosaic of the band-3 data, while the red cross shows the pointing center. Similarly, the band-4 fields of view and pointing centers are marked in blue.

In the text
thumbnail Fig. 3

650 MHz image of the gain calibrator field (J1851+005) of our observations. Circled in red at the center is the phase calibrator (J1851+005; Galactic coordinates l = 33.5°, b = 0.2°). The large and extended structure to its southeast is the supernova remnant (SNR) G33.6+0.1. The beam size is about 10 arcsec. The intensity scale is saturated to clearly bring out the weak extended structures.

In the text
thumbnail Fig. 4

Calibration scheme that we followed to reduce the data.

In the text
thumbnail Fig. 5

uv-coverage plots of typical fields observed on 2019 August 22 (left, 3SRC105), May 27 (middle, 4SRC001), and May 26 (right, 4SRC104). The asymmetry caused by non-functional antennas (marked in yellow) is clearly visible in the fields observed on August 22 and May 27, which leads to an elongation of the synthesized beam.

In the text
thumbnail Fig. 6

Example field observed in band-4 during a night with all antennas functioning, uncorrected for primary beam attenuation, in the native equatorial coordinates. The white circle marks the region used for mosaicking. The beam size is 12″ × 7″, and the rms noise is about 0.3 mJy beam−1. The colormap and the limits are chosen such that the low-level noise features and the negative sidelobes are seen clearly.

In the text
thumbnail Fig. 7

Radio images of the W43/W44 region. From top to bottom: band-3 and band-4 mosaics from this project convolved to the common beam size of 25″, their resulting spectral index map, the 1.4 GHz THOR image of the same area combined with the VGPS data (Beuther et al. 2016; Stil et al. 2006), which is also at 25″ resolution, and the 5.8 GHz GLOSTAR data (VLA D-configuration image combined with the Effelsberg image), which is at a resolution of 18″.

In the text
thumbnail Fig. 8

Cumulative distribution functions of the noise maps (left), the histograms of values of the pixels in the continuum mosaics of band-3 (middle), and band-4 (right). Gaussian least-squares fits are performed on the histograms, which gave standard deviations as ~4.9 mJy beam−1 for band-3 and ~1.8 mJy beam−1 for band-4.

In the text
thumbnail Fig. 9

Position offsets of the TGSS sources seen in our survey and the NVSS, with obvious extended sources excluded (left). Observed and expected flux density comparison is shown for sources in band-3 (middle) and band-4 (right), where the expected flux density is obtained using spectral indices from de Gasperin et al. (2018).

In the text
thumbnail Fig. 10

Probability density functions obtained using kernel density estimations of the values of spectral index pixels belonging to SNRs and H II regions.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.