A Comprehensive Catalog of Radio Sources and Rotation Measures in the Perseus Molecular Cloud from Very Large Array Observations

Abstract

We present a comprehensive radio polarization study of the Perseus molecular cloud using wideband L-band observations from the Karl G. Jansky Very Large Array. Our survey covers $\sim13.8$~deg$^2$ with a mean Stokes~$I$ sensitivity of $\sim80~μ$Jy~beam$^{-1}$, enabling the detection of 1410 compact radio sources. From this population, we construct a catalog of source properties, including positions, integrated flux densities, and spectral indices measured across nine spectral windows. The majority of sources exhibit negative spectral indices, consistent with non-thermal synchrotron emission. Using RM Synthesis and RM CLEAN techniques, we detect 205 polarized background sources above an $8σ$ threshold. This corresponds to a sampling density of $\sim14.8$~deg$^{-2}$, representing more than an order-of-magnitude increase compared to previous NVSS-based measurements. The resulting rotation measures exhibit coherent large-scale variations across the surveyed region, with additional small-scale structure superimposed. The enhanced sensitivity, frequency coverage, and sampling density of our observations enable a substantially improved mapping of the line-of-sight magnetic field component toward the Perseus molecular cloud compared to previous surveys.

Figures (15)

• Figure 1: A$_V$ map of the Perseus molecular cloud. Each circle represents a mosaic composed of 55 fields arranged in a hexagonal pattern, with each mosaic observed over five epochs. The four numbers of yellow dashed circles represent the mosaics observed with the VLA in 2019, forming the primary dataset used in this study. The eight white circles indicate the additional regions observed in 2024, aimed at enhancing the completeness of the $B_\parallel$ measurements across the cloud. The color scale in this figure represents the visual extinction ($A_V$, in magnitudes) derived from near-infrared dust extinction maps based on data from the Two Micron All Sky Survey (2MASS) and processed using the NICEST algorithm kainulainen2009probinglombardi2009nicest.
• Figure 2: Stokes I image of the Perseus molecular cloud obtained from the 19B VLA observations in D-configuration. This mosaic represents the total intensity map produced using the multi-frequency synthesis imaging technique with wide-field corrections and deconvolution applied via the tclean algorithm in CASA. The image covers the central region of the molecular cloud with a synthesized beam size of approximately 46$^{\prime\prime}$ and a pixel scale of 2.5$^{\prime\prime}$. The enhanced sensitivity achieved in this dataset ($\sim80~\mu$Jy beam$^{-1}$) allows for the detection of numerous faint background radio sources, which are critical for probing the magnetic field structure via Faraday rotation.
• Figure 3: Comparison between a $1^{\circ} \times 1^{\circ}$ region from this work (left) and the corresponding area in NVSS (right), shown on the same intensity scale. Owing to the higher sensitivity and angular resolution of our observations, the left panel reveals a substantially larger population of compact radio sources. In this region, we detect 120 sources using a $5\sigma$ source-extraction threshold, compared to 34 sources identified in NVSS.
• Figure 4: Example output from PyBDSF showing source detection and Gaussian modeling. Top-left: Original Stokes $I$ image displayed with a logarithmic stretch to emphasize faint emission. Top-right: Detected emission islands outlined in cyan. Individual Gaussian components fitted within each island are shown as magenta ellipses. Blue ellipses represent the final fitted sources, which may consist of one or more Gaussians. Bottom-left: Residual image obtained by subtracting the model from the original image, ideally containing only noise. The remaining bright pixels correspond to emission below the $5\sigma$ threshold. Bottom-right: Gaussian model image synthesized from all fitted components, with white regions indicating modeled emission.
• Figure 5: Positional offsets in RA and DEC between matched sources in our VLA catalog and a reference catalog. Each point represents the angular offset $(\Delta \mathrm{RA}, \Delta \mathrm{DEC})$ for a matched source. The distribution is tightly clustered around the origin, indicating overall good astrometric agreement, with most offsets well within one synthesized beam. The mean offsets are $\langle \Delta \mathrm{RA} \rangle = -1.\!\!^{\prime\prime}160 \pm 0.\!\!^{\prime\prime}305$ and $\langle \Delta \mathrm{DEC} \rangle = 0.\!\!^{\prime\prime}512 \pm 0.\!\!^{\prime\prime}273$. A few sources show larger deviations (up to $\sim 30^{\prime\prime}$).