Erik Rosolowsky edited section_Analysis_subsection_GMC_properties__.tex  about 8 years ago

Commit id: cdc297c532e7396900742993cac87441ad2d446d

deletions | additions      

       

\section{Analysis}  \subsection{GMC properties}  We characterize the properties of the extracted clouds in the CO data to assess whether we are identifying clouds that can be compared to GMCs seen in the Milky Way, or whether the emission structures in the M83 data are better described as Giant Molecular Associations (GMAs), which are larger scale structures of molecular gas (Rand et al.).   We identify molecular clouds in the the CO emission line data using the CPROPS algorithm (Rosolowsky & Leroy, 2006)\footnote{We use the \texttt{cpropstoo} implementation at \url{http://github.com/akleroy/cpropstoo}}. We utilize their recommended algorithm for identifying GMCs in interferometer data described as follows. The algorithm first calculates a spatially and varying estimate of the noise in the map by calculating the rms ($\sigma(\alpha,\delta)$) of signal-free channels. Emission is then identified as those pixels in the (three-dimensional) data cube that are larger than $4\sigma(\alpha,\delta)$ in two adjacent velocity channels. This emission mask is then extended to include all connected pixels which are larger than $2\sigma(\alpha,\delta)$ in two adjacent channels. The masked emission is then divided into individual molecular clouds using a seeded watershed algorithm, with individual clouds being defined by local maxima that are separated by at least XXX pc spatially or XXX km/s in projected velocity. Any pair of local maxima in the same contiguous region of the mask are also required to be at least XXX K above the saddle point of emission connecting those maxima.   We determine the macroscopic properties of the GMCs in the system by calculating moments of the emission line data. To account for emission below the edge of the emission mask, the values of the moments are extrapolated as a function of pixel value to the 0 K brightness threshold (see Rosolowsky \& Leroy 2006) for details. This extrapolation is necessary to avoid bias in low signal-to-noise data, but it can introduce substantial uncertainty (up to 100\% for clouds near the signal-to-noise cut). We calculate the CO luminosity of the molecular clouds ($L_{\mathrm{CO}}$) by integrating the emission associated with each cloud, with a quadratic extrapolation. The luminous mass is calculated by scaling by a single CO-to-H$_2$ conversion factor:  \begin{equation}  M_{\mathrm{lum}} = \alpha_{\mathrm{CO}} L_{\mathrm{CO}},  \end{equation}  where $\alpha_{CO}=4.35 M_{\odot}\mbox{ pc}^{-2}\mbox{(K km s^{-1})^{-1}}$. The velocity dispersion is the linearly-extrapolated, emission-weighted second moment of the velocity axis, corrected for the channel width. Similarly radius of the cloud is the root-mean-square of the linearly-extrapolated, emission-weighted second moments of the major and minor spatial axis of the emission. The radius is also corrected for the instrumental response by assuming an elliptical beam and subtracting its width in quadrature. See Rosolowsky \& Leroy (2006) for details. The algorithm corrects for the case where the beam and cloud position angles are not aligned.   The virial mass of the cloud is calculated from the radius and line width of the molecular cloud: $M_{\mathrm{vir}} = 5 \sigma^2 R/G$. Comparing the virial and luminous masses gives insight to the dynamical nature of the molecular clouds identified in the data. The average surface density is calculated from the luminous mass $\Sigma = M_{\mathrm{lum}}/(\pi R^2)$.  Correlations of these macroscopic properties give clues to the nature of the molecular medium. We compare the properties of the molecular clouds to those seen in the Milky Way study of \citet[][;S87]{Solomon_1987} because that work measured GMC properties using similar techniques as we do here. In Figure XXX, we correlate the GMC properties and compare the result to the trend lines seen in the S87 data. First, we see that there is good agreement between the virial and luminous masses in these clouds, and this is seen throughout the system. We code each datum with the galactocentric assuming each cloud   The figure shows that the objects we identify can be associated with Milky Way molecular clouds. They have comparable sizes and mass scales, and they show the same underlying  \citet{Solomon_1987} examined GMCs in the Milky Way and described a few key relationships that characterize the properties of these clouds. We derived these properties for the M83 clouds to be confident that they are GMCs.  The M83 GMCs exhibit a similar relationship for velocity line width $\mathrm{\sigma}$ to radius $R$ as the Milky Way GMCs, which have a fit of $\sigma=\left(\frac{\pi^{1/2}R}{3.4}\right)^{0.5}$ kms$^{-1}$ (\ref{fig:rdv}). This allows us to conclude the clouds are in virial equilibrium \cite{Solomon_1987}. The virial mass $M_\mathrm{vir}$ can then be calculated from the size and velocity dispersion, while the luminous mass $M_\mathrm{lum}$ can be calculated from the luminosity and X$_\mathrm{CO}$ \cite{Rosolowsky_2006}. These masses for M83 conform to a 1:1 ratio as expected (\ref{fig:mlummvir}). Finally, the luminous mass to radius relationships are also consistent with previous findings, with \citet{Solomon_1987} finding a fit of $M=540R^2 (M_\odot)$ (\ref{fig:rm}). From this we can infer that the observations in M83 are of GMCs.