CASA, the Common Astronomy Software Applications package, is the primary data processing software for the Atacama Large Millimeter/submillimeter Array (ALMA) and NSF's Karl G. Jansky Very Large Array (VLA), and is frequently used also for other radio telescopes. The CASA software can process data from both single-dish and aperture-synthesis telescopes, and one of its core functionalities is to support the data reduction and imaging pipelines for ALMA, VLA and the VLA Sky Survey (VLASS). CASA has recently undergone several exciting new developments, including an increased flexibility in Python (CASA 6), support of Very Long Baseline Interferometry (VLBI), performance gains through parallel imaging, data visualization with the new Cube Analysis Rendering Tool for Astronomy (CARTA), enhanced reliability and testing, and modernized documentation. These proceedings of the 2019 Astronomical Data Analysis Software & Systems (ADASS) conference give an update of the CASA project, and detail how these new developments will enhance user experience of CASA.

Imaging in radio astronomy is usually carried out with a single-dish radio telescope doing a raster scan of a region of the sky or with an interferometer that samples the visibility function of the sky brightness. Mosaic observations are the current standard for imaging large fields of view with an interferometer and multi-frequency observations are now routinely carried out with both types of telescopes to increase the continuum imaging sensitivity and to probe spectral structure. This paper describes an algorithm to combine wideband data from these two types of telescopes in a joint iterative reconstruction scheme that can be applied to spectral cube or wideband multi-term imaging both for narrow fields of view as well as mosaics. Our results demonstrate the ability to prevent instabilities and error that typically arise when wide-band or joint mosaicing algorithms are presented with spatial and spectral structure that is inadequetely sampled by the interferometer alone. For comparable noise levels in the single dish and interferometer data, the numerical behaviour of this algorithm is expected to be similar to the idea of generating artificial visibilities from single dish data. However, our discussed implementation is simpler and more flexible in terms of applying relative data weighting schemes to match noise levels while preserving flux accuracy, fits within standard iterative image reconstruction frameworks, is fully compatible with wide-field and joint mosaicing gridding algorithms that apply corrections specific to the interferometer data and may be configured to enable spectral cube and wideband multi-term deconvolution for single-dish data alone.

We have carried out VLA imaging and a Fermi timing analysis of the 115 ms gamma-ray and radio pulsar PSR J0002+6216. We found that the pulsar lies at the apex of a narrowly collimated cometary-like 7 arcmin tail of non-thermal radio emission which we identify as a bow-shock pulsar wind nebula. The tail of the nebula points back toward the geometric center of the supernova remnant CTB 1 (G116.9+0.2) 28 arcmin away, at a position angle $\theta_\mu=113^\circ$. We measure a proper motion with 2.9$\sigma$ significance from a Fermi timing analysis giving $\mu$=115$\pm$33 mas yr$^{-1}$ and $\theta_\mu=121^\circ\pm{13}^\circ$, corresponding to a large transverse pulsar velocity of 1100 km s$^{-1}$ at a distance of 2 kpc. This proper motion is of the right magnitude and direction to support the claim that PSR J0002+6216 was born from the same supernova that produced CTB 1. We explore the implications for pulsar birth periods, asymmetric supernova explosions, and mechanisms for pulsar natal kick velocities.

In this paper we examine to what extent the radio continuum can be used as an extinction free probe of star formation in dwarf galaxies. To that aim we observe $40$ nearby dwarf galaxies with the Very Large Array at 6 cm ($4$-$8$ GHz) in C-configuration. We obtained images with $3$-$8^{\prime\prime}$ resolution and noise levels of $3$-$15{\rm \,\mu Jy\,beam^{-1}}$. We detected emission associated with $22$ of the $40$ dwarf galaxies, $8$ of which are new detections. The general picture is that of an interstellar medium largely devoid of radio continuum emission, interspersed by isolated pockets of emission associated with star formation. We find an average thermal fraction of $\sim 50$-$70$% and an average magnetic field strength of $\sim 5$-$8\,{\rm \mu G}$, only slightly lower than that found in larger, spiral galaxies. At 100 pc scales, we find surprisingly high values for the average magnetic field strength of up to 50$\,{\rm \mu G}$. We find that dwarf galaxies follow the theoretical predictions of the radio continuum-star formation rate relation within regions of significant radio continuum emission but that the non-thermal radio continuum is suppressed relative to the star formation rate when considering the entire optical disk. We examine the far-infrared-star formation rate relation for our sample and find that the far-infrared is suppressed compared to the expected star formation rate. We discuss explanations for these observed relations and the impact of our findings on the radio continuum-far-infrared relation. We conclude that radio continuum emission at centimetre wavelengths has the promise of being a largely extinction-free star formation rate indicator. We find that star formation rates of gas rich, low mass galaxies can be estimated with an uncertainty of $\pm 0.2$ dex between the values of $2 \times 10^{-4}$ and $0.1 {\rm M_\odot\,yr^{-1}}$.

Context. Most popular algorithms in use to remove the effects of a telescope's point spread function (PSF) in radio astronomy are variants of the CLEAN algorithm. Most of these algorithms model the sky brightness using the delta-function basis, which results in undesired artefacts when used on image extended emission. The adaptive scale pixel decomposition (Asp-Clean) algorithm models the sky brightness on a scale-sensitive basis and thus gives a significantly better imaging performance when imaging fields that contain both resolved and unresolved emission. Aims. However, the runtime cost of Asp-Clean is higher than that of scale-insensitive algorithms. In this paper, we identify the most expensive step in the original Asp-Clean algorithm and present an efficient implementation of it, which significantly reduces the computational cost while keeping the imaging performance comparable to the original algorithm. The PSF sidelobe levels of modern wide-band telescopes are significantly reduced, allowing us to make approximations to reduce the computing cost, which in turn allows for the deconvolution of larger images on reasonable timescales. Methods. As in the original algorithm, scales in the image are estimated through function fitting. Here we introduce an analytical method to model extended emission, and a modified method for estimating the initial values used for the fitting procedure, which ultimately leads to a lower computational cost. Results.The new implementation was tested with simulated EVLA data and the imaging performance compared well with the original Asp-Clean algorithm. Tests show that the current algorithm can recover features at different scales with lower computational cost.

Young stellar objects are known to exhibit strong radio variability on timescales of weeks to months, and a few reports have documented extreme radio flares, with at least an order of magnitude change in flux density on timescales of hours to days. However, there have been few constraints on the occurrence rate of such radio flares or on the correlation with pre-main sequence X-ray flares, although such correlations are known for the Sun and nearby active stars. Here we report simultaneous deep VLA radio and Chandra X-ray observations of the Orion Nebula Cluster, targeting hundreds of sources to look for the occurrence rate of extreme radio variability and potential correlation with the most extreme X-ray variability. We identify 13 radio sources with extreme radio variability, with some showing an order of magnitude change in flux density in less than 30 minutes. All of these sources show X-ray emission and variability, but only on timescales <1h do we find clear correlations with extreme radio flaring. Strong X-ray variability does not predict the extreme radio sources and vice versa. Radio flares thus provide us with a new perspective on high-energy processes in YSOs and the irradiation of their protoplanetary disks. Finally, our results highlight implications for interferometric imaging of sources violating the constant-sky assumption.

Next generation radio telescope arrays are being designed and commissioned to accurately measure polarized intensity and rotation measures across the entire sky through deep, wide-field radio interferometric surveys. Radio interferometer dish antenna arrays are affected by direction-dependent (DD) gains due to both instrumental and atmospheric effects. In this paper we demonstrate the effect of DD errors for parabolic dish antenna array on the measured polarized intensities of radio sources in interferometric images. We characterize the extent of polarimetric image degradation due to the DD gains through wide-band VLA simulations of representative point source simulations of the radio sky at L-Band(1-2GHz). We show that at the 0.5 gain level of the primary beam (PB) there is significant flux leakage from Stokes $I$ to $Q$, $U$ amounting to 10\% of the total intensity. We further demonstrate that while the instrumental response averages down for observations over large parallactic angle intervals, full-polarization DD correction is required to remove the effects of DD leakage. We also explore the effect of the DD beam on the Rotation Measure(RM) signals and show that while the instrumental effect is primarily centered around 0 rad-m$^{-2}$, the effect is significant over a broad range of RM requiring full polarization DD correction to accurately reconstruct RM synthesis signal.

Many aspects of the Atacama Large Millimeter Array (ALMA) instrument are still unknown due to its young age. One such aspect is the true nature of the primary beam of each baseline, and how changes to the individual primary beams affect astronomical observations when said changes are ignored during imaging. This paper aims to create a more thorough understanding of the strengths and weaknesses of ALMA through realistic modeling of the primary beams and simulated observations, which in turn can inform the user of the necessity of implementing more computationally costly algorithms, such as A-Projection, and when simpler, quicker algorithms will suffice. We quantify our results by examining the dynamic range of each observation, along with the ability to reconstruct the Stokes I amplitude of the test sources. These tests conclude that for dynamic ranges of less than 1000, for point sources and sources much smaller than the main lobe of the primary beam, the accuracy of the primary beam model beyond the physical size of the aperture simply doesn't matter. In observations of large extended sources, deconvolution errors dominate the reconstructed images and the individual primary beam errors were indistinguishable from each other.

We present $\simeq$0$.\!\!^{\prime\prime}4$-resolution extinction-independent distributions of star formation and dust in 11 star-forming galaxies (SFGs) at $z = 1.3-3.0$. These galaxies are selected from sensitive, blank-field surveys of the $2' \times 2'$ Hubble Ultra-Deep Field at $\lambda = 5$ cm and 1.3 mm using the Karl G. Jansky Very Large Array (VLA) and Atacama Large Millimeter/submillimeter Array (ALMA). They have star-formation rates (SFRs), stellar masses, and dust properties representative of massive main-sequence SFGs at $z \sim 2$. Morphological classification performed on spatially-resolved stellar mass maps indicates a mixture of disk and morphologically disturbed systems; half of the sample harbor X-ray active galactic nuclei (AGN), thereby representing a diversity of $z \sim 2$ SFGs undergoing vigorous mass assembly. We find that their intense star formation most frequently occurs at the location of stellar-mass concentration and extends over an area comparable to their stellar-mass distribution, with a median diameter of $4.2 \pm 1.8$ kpc. This provides direct evidence for galaxy-wide star formation in distant, blank-field-selected main-sequence SFGs. The typical galactic-average SFR surface density is 2.5 M$_{\odot}$yr$^{-1}$kpc$^{-2}$, sufficiently high to drive outflows. In X-ray-selected AGN where radio emission is enhanced over the level associated with star formation, the radio excess pinpoints the AGN, which are found to be co-spatial with star formation. The median extinction-independent size of main-sequence SFGs is two times larger than those of bright submillimeter galaxies whose SFRs are $3-8$ times larger, providing a constraint on the characteristic SFR ($\sim300$ M$_{\odot}$yr$^{-1}$) above which a significant population of more compact star-forming galaxies appears to emerge.

Many deep wide-band wide-field radio interferometric surveys are being designed to accurately measure intensities, spectral indices and polarization properties of faint source populations. In this paper we compare various wideband imaging methods to evaluate the accuracy to which intensities and spectral indices of sources close to the confusion limit can be reconstructed. We simulated a wideband single-pointing (C-array, L-Band (1-2GHz)) and 46-pointing mosaic(D-array, C-Band (4-8GHz)) JVLA observation using realistic brightness distribution ranging from $1\mu$Jy to $100m$Jy and time-,frequency-, polarization- and direction-dependent instrumental effects. The main results from these comparisons are (a) errors in the reconstructed intensities and spectral indices are larger for weaker sources even in the absence of simulated noise, (b) errors are systematically lower for joint reconstruction methods (such as MT-MFS) along with A-Projection for accurate primary beam correction, and (c) use of MT-MFS for image reconstruction eliminates Clean-bias (which is present otherwise). Auxiliary tests include solutions for deficiencies of data partitioning methods (e.g. the use of masks to remove clean bias and hybrid methods to remove sidelobes from sources left undeconvolved), the effect of sources not at pixel centers and the consequences of various other numerical approximations within software implementations. This paper also demonstrates the level of detail at which such simulations must be done in order to reflect reality, enable one to systematically identify specific reasons for every trend that is observed and to estimate scientifically defensible imaging performance metrics and the associated computational complexity of the algorithms/analysis procedures.

We present a deep centimeter-wavelength catalog of the Orion Nebula Cluster (ONC), based on a 30h single-pointing observation with the Karl G. Jansky Very Large Array in its high-resolution A-configuration using two 1 GHz bands centered at 4.7 GHz and 7.3 GHz. A total of 556 compact sources were detected in a map with a nominal rms noise of 3 muJy/bm, limited by complex source structure and the primary beam response. Compared to previous catalogs, our detections increase the sample of known compact radio sources in the ONC by more than a factor of seven. The new data show complex emission on a wide range of spatial scales. Following a preliminary correction for the wideband primary-beam response, we determine radio spectral indices for 170 sources whose index uncertainties are less than +/-0.5. We compare the radio to the X-ray and near-infrared point-source populations, noting similarities and differences.

Superbubbles are crucial for stellar feedback, with supposedly high (of the order of 10 per cent) thermalization rates. We combined multiband radio continuum observations from the Very Large Array (VLA) with Effelsberg data to study the non-thermal superbubble (NSB) in IC 10, a starburst dwarf irregular galaxy in the Local Group. Thermal emission was subtracted using a combination of Balmer H$\alpha$ and VLA 32 GHz continuum maps. The bubble's non-thermal spectrum between 1.5 and 8.8 GHz displays curvature and can be well fitted with a standard model of an ageing cosmic ray electron population. With a derived equipartition magnetic field strength of $44\pm 8 \rm\mu G$, and measuring the radiation energy density from Spitzer MIPS maps as $5\pm 1\times 10^{-11} \rm erg\, cm^{-3}$, we determine, based on the spectral curvature, a spectral age of the bubble of $1.0\pm 0.3 \rm Myr$. Analysis of the LITTLE THINGS HI data cube shows an expanding HI hole with 100 pc diameter and a dynamical age of $3.8\pm 0.3 \rm Myr$, centred to within 16 pc on IC 10 X-1, a massive stellar mass black hole ($M > 23 M_\odot$). The results are consistent with the expected evolution for a superbubble with a few massive stars, where a very energetic event like a Type Ic supernova/hypernova has taken place about 1 Myr ago. We discuss alternatives to this interpretation.

Synthesis imaging in radio astronomy is affected by instrumental and atmospheric effects which introduce direction-dependent (DD) gains.The antenna power pattern varies both as a function of time and frequency. The broad band time varying nature of the antenna power pattern when not corrected leads to gross errors in full Stokes imaging and flux estimation. In this poster we explore the errors that arise in image deconvolution while not accounting for the time and frequency dependence of the antenna power pattern. Simulations were conducted with the wide-band full Stokes power pattern of the Karl G. Jansky Very Large Array (VLA) antennas to demonstrate the level of errors arising from direction-dependent gains and their non-neglegible impact on upcoming sky surveys such as the VLASS. DD corrections through hybrid projection algorithms are computationally expensive to perform. A highly parallel implementation through high performance computing architectures is the only feasible way of applying these corrections to the large data sizes of these upcoming surveys.

We report new observations of Abell 2256 with the Karl G. Jansky Very Large Array (VLA) at frequencies between 1 and 8 GHz. These observations take advantage of the 2:1 bandwidths available for a single observation to study the spectral index, polarization and Rotation Measure as well as using the associated higher sensitivity to image total intensity features down to ~0.5" resolution. We find the Large Relic, which dominates the cluster, is made up of a complex of filaments which show correlated distributions in intensity, spectral index, and fractional polarization. The Rotation Measure varies across the face of the Large Relic but is not well correlated with the other properties of the source. The shape of individual filaments suggests that the Large Relic is at least 25 kpc thick. We detect a low surface brightness arc connecting the Large Relic to the Halo and other radio structures suggesting a physical connection between these features. The center of the F-complex is dominated by a very steep-spectrum, polarized, ring-like structure, F2, without an obvious optical identification, but the entire F-complex has interesting morphological similarities to the radio structure of NGC1265. Source C, the Long Tail, is unresolved in width near the galaxy core and is </~100pc in diameter there. This morphology suggests either that C is a one-sided jet or that the bending of the tails takes place very near the core, consistent with the parent galaxy having undergone extreme stripping. Overall it seems that many of the unusual phenomena can be understood in the context of Abell 2256 being near the pericenter of a slightly off-axis merger between a cluster and a smaller group. Given the lack of evidence for a strong shock associated with the Large Relic, other models should be considered, such as reconnection between two large-scale magnetic domains.

A large number of deep and wide-field radio interferometric surveys are being designed to measure accurate statistics of faint source populations. Most require mosaic observations, and expect to benefit from the sensitivity provided by broad-band instruments. In this paper, we present preliminary results from a comparison of several wideband imaging methods in the context of how accurately they reconstruct the intensities and spectral indices of micro-Jy level sources.

Variations of the antenna primary beam (PB) pattern as a function of time, frequency and polarization form one of the dominant direction-dependent effects at most radio frequency bands. These gains may also vary from antenna to antenna. The A-Projection algorithm, published earlier, accounts for the effects of the narrow-band antenna PB in full polarization. In this paper we present the Wide-Band A-Projection algorithm (WB A-Projection) to include the effects of wide bandwidth in the A-term itself and show that the resulting algorithm simultaneously corrects for the time, frequency and polarization dependence of the PB. We discuss the combination of the WB A-Projection and the Multi-term Multi Frequency Synthesis (MT-MFS) algorithm for simultaneous mapping of the sky brightness distribution and the spectral index distribution across a wide field of view. We also discuss the use of the narrow-band A-Projection algorithm in hybrid imaging schemes that account for the frequency dependence of the PB in the image domain.

The spatial-frequency coverage of a radio interferometer is increased by combining samples acquired at different times and observing frequencies. However, astrophysical sources often contain complicated spatial structure that varies within the time-range of an observation, or the bandwidth of the receiver being used, or both. Image reconstruction algorithms can been designed to model time and frequency variability in addition to the average intensity distribution, and provide an improvement over traditional methods that ignore all variability. This paper describes an algorithm designed for such structures, and evaluates it in the context of reconstructing three-dimensional time-varying structures in the solar corona from radio interferometric measurements between 5 GHz and 15 GHz using existing telescopes such as the EVLA and at angular resolutions better than that allowed by traditional multi-frequency analysis algorithms.

We exploit the vastly increased sensitivity of the Expanded Very Large Array (EVLA) to study the radio continuum and polarization properties of the post-starburst, dwarf irregular galaxy IC10 at 6 cm, at a linear resolution of ~50 pc. We find close agreement between radio continuum and Halpha emission, from the brightest HII regions to the weaker emission in the disk. A quantitative analysis shows a strictly linear correlation, where the thermal component contributes 50% to the total radio emission, the remainder being due to a non-thermal component with a surprisingly steep radio spectral index of between -0.7 and -1.0 suggesting substantial radiation losses of the cosmic-ray electrons. We confirm and clearly resolve polarized emission at the 10-20% level associated with a non-thermal superbubble, where the ordered magnetic field is possibly enhanced due to the compression of the expanding bubble. A fraction of the cosmic-ray electrons has likely escaped because the measured radio emission is a factor of 3 lower than what is suggested by the Halpha inferred SFR.

The radio continuum emission from the Galaxy has a rich mix of thermal and non-thermal emission. This very richness makes their interpretation challenging since the low radio opacity means that a radio image represents the sum of all emission regions along the line-of-sight. These challenges make the existing narrow-band radio surveys of the Galactic plane difficult to interpret: e.g. a small region of emission might be a supernova remnant (SNR) or an HII region, or a complex combination of both. Instantaneous wide bandwidth radio observations in combination with the capability for high resolution spectral index mapping, can be directly used to disentangle these effects. Here we demonstrate simultaneous continuum and spectral index imaging capability at the full continuum sensitivity and resolution using newly developed wide-band wide-field imaging algorithms. Observations were done in the L- and C-Band with a total bandwidth of 1 and 2 GHz respectively. We present preliminary results in the form of a full-field continuum image covering the wide-band sensitivity pattern of the EVLA centered on a large but poorly studied SNR (G55.7+3.4) and relatively narrower field continuum and spectral index maps of three fields containing SNR and diffused thermal emission. We demonstrate that spatially resolved spectral index maps differentiates regions with emission of different physical origin (spectral index variation across composite SNRs and separation of thermal and non-thermal emission), superimposed along the line of sight. The wide-field image centered on the SNR G55.7+3.4 also demonstrates the excellent wide-field wide-band imaging capability of the EVLA.

Aims : We describe MS-MFS, a multi-scale multi-frequency deconvolution algorithm for wide-band synthesis-imaging, and present imaging results that illustrate the capabilities of the algorithm and the conditions under which it is feasible and gives accurate results. Methods : The MS-MFS algorithm models the wide-band sky-brightness distribution as a linear combination of spatial and spectral basis functions, and performs image-reconstruction by combining a linear-least-squares approach with iterative $\chi^2$ minimization. This method extends and combines the ideas used in the MS-CLEAN and MF-CLEAN algorithms for multi-scale and multi-frequency deconvolution respectively, and can be used in conjunction with existing wide-field imaging algorithms. We also discuss a simpler hybrid of spectral-line and continuum imaging methods and point out situations where it may suffice. Results : We show via simulations and application to multi-frequency VLA data and wideband EVLA data, that it is possible to reconstruct both spatial and spectral structure of compact and extended emission at the continuum sensitivity level and at the angular resolution allowed by the highest sampled frequency.

This paper summarizes some of the major calibration and image reconstruction techniques used in radio interferometry and describes them in a common mathematical framework. The use of this framework has a number of benefits, ranging from clarification of the fundamentals, use of standard numerical optimization techniques, and generalization or specialization to new algorithms.