giant molecular clouds , both galactic and extragalactic , are observed to have velocity gradients that many authors have interpreted as being caused by rotation ( e.g. , kutner et al .
1977 , phillips 1999 ; rosolowsky et al . 2003 ) .
if we start with the premise that these clouds are rotating because they have inherited the angular momentum of the rotating galactic disk out of which they formed , conservation of angular momentum should provide clues that give us insight to the origin of their formation
. however , simple formation theories that assume giant molecular clouds ( gmcs ) form by condensing out of the galactic disk are at odds with some of the observations .
for instance , they do not adequately explain why the directions of gmc velocity gradients are not typically aligned with the direction of galactic rotation , counter to the expectation of conservation of angular momentum .
furthermore , simple formation scenarios tend to overpredict the observed specific angular momentum o f gmcs ( e.g. , blitz 1990 ; rosolowsky et al .
this is the so - called `` angular momentum problem . ''
provided there is no transfer of angular momentum , the angular momentum of a gmc should be equal to that of the gas out of which it formed .
but blitz ( 1990 ) , working under the assumption that the velocity gradients in molecular clouds are due to rotation , showed that the angular momentum due to galactic differential rotation in the solar neighborhood interstellar medium is consistently greater than that contained within molecular clouds .
even molecular clouds with the largest observed velocity gradients  such as the rosette and orion a molecular clouds  have less specific angular momentum compared to the ism from which they presumably formed . blitz ( 1990 ) also pointed out that because the molecular clouds in his sample are rotating in a sense _ opposite _ to that of galactic rotation , they could not have conserved angular momentum from the initial states he calculated , unless the local galactic rotation curve is falling , or unless the clouds always collapsed azimuthall y.    the angular momentum problem also extends to extragalactic molecular clouds .
roso - lowsky et al .
( 2003 ) showed that simple gmc formation theories consistently overestimate the magnitude of the observed angular momentum of molecular clouds in the galaxy m33 . on average
, they found that simple theory overpredicts the observed magnitudes of specific angular momentum by more than a factor of 5 .
furthermore , they found that 40% of the gmcs in m33 are counterrotating with respect to the galactic plane .    in this paper
, we shed light on the angular momentum problem by doing a detailed analysis of the kinematics in gmcs _ and _ the surrounding ism .
whereas previous studies estimated the initial angular momentum imparted to gmcs from the galactic rotation curve ( e.g. , blitz 1990 ; blitz 1993 ) , we compare the velocity fields of gmcs to those of the ism with which they are associated directly from observation . in light of the observation that gmcs have a spatial and kinematic correlation with high - surface density atomic gas (  [ sec : analysis ] ) , we pose the question : does the rotation of the large - scale @xmath0  associated with gmcs mirror that of the gmcs themselves ?
our primary goal is to determine whether or not rotation is the cause of the velocity gradients in gmcs . to that end , we create first - moment maps of the molecular clouds in our sample and of the atomic gas surrounding them are for comparison . in the following section , we describe the and 21-cm data used to conduct this study . in  3 , we create intensity - weighted first moment and @xmath0  maps of five galactic clouds : perseus , orion a , ngc 2264 , monoceros r2 ( monr2 ) , and the rosette .
the results from these measurements are given in  4 , and a summary is provided in  5 .
to measure velocity gradients and other properties across the molecular clouds , we use high - resolution , high - sensitivity published observations . because emission is nearly always optically thin in galactic gmcs , we have the advantage of getting a detailed view of the kinematic structure of the molecular clouds in our sample . and because has narrower line widths than the optically thick , the former permits finer separation of velocity components than the latter .
nevertheless , the large scale velocity gradient of a gmc measured using is generally consistent with that measured using , since the gradient is being measured across several parsecs and small scale variations in the velocity field tend to get averaged out .
the data for perseus , orion a , ngc 2264 , and the monr2 molecular clouds were generously provided by j. bally ( see bally et al .
observations of these clouds were taken at the at&t bell laboratories 7 m telescope and have a beam size of 100 .
perseus , orion a , and ngc
2264 data were resampled onto @xmath2 grids .
the monr2 data were resampled onto a @xmath3 grid .
the spectral resolution of 128 channels at 100 khz corresponds to a velocity resolution of 0.27 . on the @xmath4 scale
, the cubes have rms noise levels of 0.17 k ( perseus ) , 0.32 k ( orion a ) , 0.42 ( ngc 2264 ) and 0.29 k ( monr2 ) .
j. williams and m. heyer graciously provided the fcrao data of the rosette molecular cloud ( see heyer , williams , & brunt 2006 ) . the beam size is @xmath5 and the data were interpolated onto a @xmath6 grid .
the spectral resolution is 59 khz per channel , and the velocity resolution is 0.133 .
the data have an rms noise of 0.21 k in @xmath4 , similar to that of the bell labs data .
the @xmath0  data are obtained from the leiden / argentine / bonn ( lab ) galactic @xmath0  survey ( kalberla et al .
2005 ) , which spans velocities from @xmath7  to @xmath8 .
( the lsr velocities of the molecular clouds in our sample range from 6 to 14  with velocity dispersions of a few . )
the survey has a half power beam width of @xmath9 , velocity resolution of 1.3 , and an rms noise level of 0.07 k. the high sensitivity and resolution of the lab data set enables a detailed study of the atomic gas from which the gmcs formed .
observations in the milky way indicate the molecular clouds are typically associated with high - density atomic gas with column densities around @xmath10 ( e.g. , mckee & ostriker 2007 ) . the rosette molecular cloud is a prototypical example : williams et al .
( 1995 ) measured the mean column of the @xmath0  associated with the gmc to be @xmath11 .
furthermore , in external galaxies , gmcs are often observed to be located on or near bright @xmath0  peaks . in the large magellanic cloud , mizuno et al .
( 2001 ) observed that most gmcs are associated with @xmath0  having column densities greater than @xmath12 .
rosolowsky et al .
( 2003 ) found that every gmc they identified in m33 lies on an over - dense @xmath0  filament , though every over - dense @xmath0  region does not contain a gmc .
this seems to imply that high column density atomic gas is necessary but not sufficient for gmc formation .
it is our goal to do a detailed comparison of the velocity fields in the gmcs and local atomic gas . to that end , we first describe how the physical properties  including column densities , masses , and velocity gradients  in both the atomic and molecular gas are determined .
we then provide our criteria for choosing spatial and kinematic regions of atomic gas associated with the molecular clouds .
lastly , in this section , we discuss how we estimate the specific angular momentum .
these parameters are summarized in tables 1 and 2 .
the @xmath0  column density , @xmath13 , is calculated along each line of sight by integrating the atomic hydrogen emission above a certain background value and over the selected velocity range , ( see  3.1 and table 1 ) , @xmath14 where @xmath15 is the brightness temperature of the @xmath0  observations and @xmath16 is the channel velocity width .
assuming that the @xmath0  is optically thin , equation [ eq : nhi ] provides a lower limit to measured column density , which we convert into units of surface density in figures [ fig : oriona0 ]  [ fig : rosette0 ] .
the total @xmath0  mass is then determined by summing over all pixels in the map where emission is detected , that is , where the emission is at least three times the root - mean - square ( rms ) noise level : @xmath17 where @xmath18 is the correction for helium , @xmath19 is the distance to a given molecular cloud , and @xmath20 is the area of one pixel , which corresponds to one resolution element .    the column density , @xmath21 , is derived assuming that the  emission is optically thin and in local thermodynamic equilibrium . following frerking et al .
( 1982 ) , @xmath22^{-1 }   \nonumber
\\          \times \int_{v_{\mathrm{min}}}^{v_{\mathrm{max}}}{\frac { t_{{\mathrm{b , co}}}~{\ensuremath{\mbox{d}{v}}}}{{\mbox{k km s$^{-1}$ } } } } ~{\mbox{cm$^{-2}$ } } ,      \end{aligned}\ ] ] where @xmath23 is the integrated intensity and @xmath24 is the excitation temperature .
normally , @xmath24 is determined by measuring the radiation temperature toward peaks .
since we lack observations at the same resolution as the data , we use a constant excitation temperature in the calculation of @xmath21 for each of the five gmcs .
based on the following arguments , we adopt a value of 20 k for each of the gmcs . if the actual excitation temperature in a given region is between 10 and 30 k , the derived column density will be in error by less than a factor of 2 .
castets et al .
( 1990 ) showed that the emission in orion a mainly arises from regions where @xmath25 k , and dense cloud cores have temperatures of @xmath26 k. nagahama et al .
( 1998 ) showed that , with the exception of two peaks at @xmath27 and @xmath28 associated with embedded young stellar groups , @xmath24 rises slowly and monotonically in galactic longitude , with an average ranging from @xmath29 k. in his study of molecular clouds , including perseus , orion a , and monr2 , carpenter ( 2000 ) adopts a constant value of @xmath30 k , though the coefficient in his formula for @xmath21 yields slightly higher values than ours in equation [ eq : nco13 ] .
thus , for orion a , perseus , ngc 2264 , and monoceros r2 , we use @xmath31 k in our calculation of @xmath21 .
williams et al .
( 1995 ) showed that @xmath24 in the rosette decreases slowly with increasing distance from the rosette nebula ( centered at @xmath32 ) from @xmath33 to @xmath34 k. in our analysis , we adopt a uniform value of @xmath31 k for the rosette .
our estimate of the cloud s mass using this value is slightly lower than that estimated by williams et al .
( 1995 ) , who measured the mass over a larger surface area ( see below ) .
next , the @xmath35  column density is evaluated assuming a ratio of @xmath36 ( frerking et al .
pixels having values at least three times the rms noise level are counted as detected .
finally , the molecular mass @xmath37 is calculated over the areas where emission is detected ( that is , higher than the 3-@xmath38 level ) using @xmath39 where @xmath40 is the mass of an @xmath35  molecule .
we note that the areas over which emission is detected and , subsequently , the masses we calculate will be smaller than cited in previous studies in which these quantities were measured using emission .
this is because the stronger line is observed over larger areas in gmcs than the line . in the rosette ,
for instance , we measure a projected area of 1500 @xmath41 and a mass of @xmath42 @xmath43 , while williams et al .
( 1995 ) measure 2200 @xmath41 and @xmath44 @xmath43 . in table 1 , both the masses and previously measured masses are listed for the gmcs .
velocity gradients are measured from first moment maps of the atomic and molecular gas .
first , the intensity - weighted average velocity along each line of sight is determined using @xmath45 where @xmath46 and @xmath47 are the velocity and temperature at location @xmath48 . following goodman et al .
( 1993 ) , the uncertainty of a given measurement is @xmath49 where @xmath50 is the spectrum noise , @xmath51 is the maximum temperature along the line of sight , and @xmath52 is the fwhm linewidth of the spectrum along the line of sight .
a plane is then fitted to the first moment map of velocity centroids , as in goodman et al .
( 1993 ) , assuming a linear velocity gradient : @xmath53 where @xmath54 is the mean cloud velocity , and @xmath55 is an arbitrary reference position , which we take to be the center of our maps , and the coefficients are @xmath56 the gradient magnitude and direction , @xmath57 and @xmath58 , are derived from the coefficients to the fit , @xmath59 @xmath60 where @xmath58 , measured in degrees east from north , points in the direction of increasing velocity .
note that since we have no information regarding the inclination of a given cloud , @xmath48 , to our line of sight , gradient measurements are underestimates of the actual values , @xmath61 .
the uncertainties in these values are calculated by propagating the errors in the coefficients . to check whether planes are good fits to the velocity centroid maps , we make plots of the central velocity at a given location in the cloud versus the perpendicular offset from the cloud s rotation axis .
this is done by taking the average velocity along lines parallel to the rotation axis at various distances . in most cases , as will be discussed in
[ sec : results ] , these plots show that planes are good fits to the velocity centroid maps .
we select @xmath0  regions in the position - velocity lab data cube that are centered , spatially and kinematically , on the five gmcs in our sample .
we do not know _ a priori _ the kinematic or spatial extent of the @xmath0  regions that are associated with each gmc nor the extent of the @xmath0  velocity gradients .
molecular clouds have well - defined boundaries at which the molecules are dissociated by uv radiation and where there is a distinct transition from primarily molecular to primarily atomic gas ( e.g. savage et al . 1977 ; blitz & thaddeus 1980 ) .
the atomic gas associated with gmcs does not have such distinct boundaries , however , making it difficult to distinguish @xmath0  that may be related to gmcs from background emission .
thus , we begin by examining regions in position - velocity space that are far from the center of the gmc . while we want to capture the full extent of any linear velocity gradient we may measure in a given region of atomic gas , we do not want to make our aperture so large that we end up including in our measurements too much atomic gas that is unrelated to the molecular clouds .
we start by varying the size of the region ( that is , the subcube extracted from the lab data ) and examine how the velocity gradient magnitude and direction change .
we fix the velocity range ( see below ) and vary the spatial size of the region centered on the gmc from about 10 to several tens of parsecs , in increments of 10 pc .
we find that the gradient direction remains roughly constant until the radius of the region over which it is measured reaches @xmath62 pc , independent of the size of the gmc . beyond this ,
the measurements start to fluctuate , as the gradient in the vicinity of the molecular cloud becomes washed out by unrelated features .
thus for each cloud , we end by selecting a spatial boundary of atomic gas which extends roughly 30  50 pc from the center of the co emission .
because the peaks of the @xmath0  regions are included in the maps we generate , varying the size of the region within @xmath62 pc does not change the gradient direction by more than a few degrees .
the 1-@xmath63 uncertainty level of @xmath58 ranges from 3 to 6 degrees for the clouds in our sample .
our criterion is supported by previous studies such as that of andersson et al .
( 1991 ) who found that the spatial extent of high - intensity @xmath0  halos , measured from the edge of the molecular cloud , ranges 5  10 pc ( see also  [ sec : discussion ] ) .
figures [ fig : oriona0 ] through [ fig : rosette0 ] show the @xmath0  surface density maps derived from the zeroth - moment intensity maps , with the outline of the molecular clouds overlaid .
also overplotted in each figure is a dashed circle indicating the @xmath0  region selected for the analysis .
keep in mind , we are showing the emission of the molecular clouds , and so the maps in figures [ fig : oriona0 ]  [ fig : rosette0 ] do not show the full extent of the molecular emission in the gmcs .    to choose relevant velocities of the atomic gas associated each gmc , we begin by examining the @xmath0  emission in the velocity range @xmath64 centered about the mean lsr velocity of the emission .
again , this is because we want to be sure that our measurements include as much as possible of the associated atomic gas . in the direction of the rosette , for instance
, the @xmath0  has been observed to extend several  beyond the co emission ( williams et al .
studies of both milky way molecular clouds ( e.g. , wannier et al .
1983 , williams et al . 1995 ) and
extragalactic clouds ( e.g. engargiola et al .
2003 ) have shown that the @xmath0  emission line profile tends to peak in the direction of gmcs . in effect
, we are using the @xmath0  velocity as a proxy for distance in order to associate the atomic gas with the gmcs . as figures [ fig : oriona0 ]  [ fig : rosette0 ] show , although the @xmath0  emission line is broader than the co line , the velocity difference between the peaks in the respective lines never exceeds @xmath65 , where @xmath65 is the velocity dispersion of the @xmath0  profile .
the bottom panel in figures [ fig : oriona0 ] through [ fig : rosette0 ] show the average spectrum through each gmc with the @xmath0  spectrum in the same direction ( within the dashed circle ) overplotted .
because the @xmath0  peaks are nearly coincident in each case , this indicates that most of the @xmath0  in the direction of a given cloud is associated with that cloud within the selected velocity range . in the cases of perseus ,
orion a , and monr2 , the @xmath0  emission drops abruptly beyond @xmath66 of the  emission .
ngc 2264 and the rosette have more complicated @xmath0  spectra , each showing double peaks that may be indicative expansion or of a blended , possibly unrelated component .
the latter explanation would not be surprising since , of the five clouds in the sample , ngc 2246 and the rosette are located closest to the galactic plane where line - of - sight blending is more of a problem .    based on figures [ fig : oriona0 ]  [ fig : rosette0 ] , we determine the boundaries of the @xmath0  emission we will use for the subsequent analysis . for perseus , orion a , and monr2 , the @xmath0  line profiles are approximated as gaussians and we assume that @xmath0  emission having velocities within @xmath67 is associated with a given molecular cloud . since both ngc 2264 and the rosette each have a second peak in their @xmath0  temperature profiles ( @xmath68 ) at higher velocities and because we want to be careful to exclude as much unrelated emission as possible , we set a slightly more stringent criterion on the velocities we select . for each cloud ,
a maximum velocity is identified at the local minimum in @xmath15 where the gaussians overlap . in the ngc 2264 spectrum , for instance
, @xmath15 drops to 29 k at 15  and then peaks again at around 20  ( see figure [ fig : ngc22640 ] ) .
thus , we eliminate all emission having velocities above 15 , the locat ion of the local minimum . all of these selections are listed in the third to last column of table 1 .
once the magnitudes of the velocity gradients in the molecular clouds and the surrounding @xmath0  are measured , we may calculate and compare their specific angular momenta , under the assumption that the linear gradients are due to solid body rotation .
the specific angular momentum , @xmath69 , is simply the total angular momentum of a body divided by its mass , @xmath70 where @xmath71 is the radius of the region , and the constant @xmath72 takes into account the moment of inertia of a rotating body . for roughly spherical gmcs having constant surface mass density distributions , @xmath73 . unless otherwise stated , we take the size of a given molecular cloud to be its effective radius , as defined by its projected area : @xmath74
.    we would like to estimate the expected specific angular momentum initially imparted to a gmc by the ism from which it forms .
this depends on the process of gmc formation , and we assume here that gmcs form via a `` top - down '' formation mechanism .
for instance , it has been suggested that molecular cloud formation occurs when an instability triggers collapse or condensation from the galactic disk ( e.g. , mouschovias et al .
1974 ; blitz & shu 1980 ; elmegreen 1982 ; kim et al .
blitz & shu ( 1980 ) show that `` bottom - up '' formation of gmcs via the random agglomeration of pre - existing low - mass clouds is unlikely because of the long timescales for this process .
gravitational instabilities and magneto - gravitational instabilities , however , tend to proceed more quickly .
if the @xmath0  surrounding the gmcs is reflective of the ism out of which the gmcs initially formed , the quantity @xmath75 provides an estimate of the initial specific angular momentum of the gmcs , assuming they are rotating . henceforth , we will often refer to this as the _ expected _ specific angular momentum .
note that this disregards the possible effects of magnetic fields ; that is , we are assuming that the magnetic field strength of the forming cloud is 0 .
the initial angular velocity imparted to a given gmc is the local value of @xmath57 , which we calculate from the first - moment maps of the @xmath0  using equation [ eq : omega ] . in principle , the size of the region from which a forming cloud gathers material , the `` accumulation radius '' @xmath76 , could have a range of values because it depends on details of the formation mechanism , as well as on assumptions regarding the initial surface density and geometry of the gas from which a gmc formed .
nevertheless , we can set an effective lower limit on @xmath76 by calculating the size of the region from which a gmc must contract to obtain its present mass , @xmath77
.    following blitz ( 1993 ) , let us assume that the initial geometry of the collapsing region is a cylinder with a diameter equal to its height .
the size of the cylinder is determined by requiring that the mass contained within it is equal to the present mass of the gmc : @xmath78 where @xmath79 is the mean volume density of the atomic gas from which the gmc formed .
blitz ( 1993 ) estimated @xmath76 using the mean value of @xmath80 in the galactic plane near the sun , 5  ( henderson et al . 1982 ) .
this corresponds to a mass density of @xmath81 ( or a number density of 0.5 @xmath82 ) , using an effective scale height of atomic gas of 200 pc in the solar vicinity ( falgarone & lequeux 1973 ) .
however , using larger values of @xmath80 ( leading to lower estimates of @xmath76 ) might be more appropriate , given the observation that gmcs tend to form in regions of @xmath0  with densities higher than global galactic values ( e.g. , engargiola et al .
2003 ; imara et al .
we estimate the accumulation radii using @xmath83 , a value more in keeping observations of the atomic gas associated with gmcs .
for instance , around the rosette , williams et al .
( 1995 ) measure an @xmath0  column density of @xmath84 , which corresponds to a @xmath85 , twice the mean galactic value . from the sancisi ( 1974 )
study of atomic gas near perseus , the inferred surface density is 11 .
our estimates of @xmath76 are listed in table 2 , as well as the effective radii of the gmcs .
finally , we note that the predicted specific angular momentum in the atomic gas has a marked dependence on the gmc mass . substituting @xmath86 from equation 2.12 for the radius in equation 2.11 yields @xmath87 .
the masses we calculate using the observations underestimate the total molecular mass of the gmcs . for this reason
, we use the larger gmc masses , as previously measured using observations , to predict the initial specific angular momentum imparted to forming molecular clouds .
typically , @xmath88 .
figures [ fig : oriona ] through [ fig : rosette ] show the @xmath0  velocity maps ( grayscale ) overlaid with velocity maps of ( color ) for each gmc .
overplotted on these maps are axes of rotation : the lines perpendicular to the gradient directions , @xmath89 and @xmath90 , as calculated from equation [ eq : theta ] .
these lines are the position angles of the rotation axes , @xmath91 , ( where @xmath92 ) , of the molecular and atomic gas _ if _ the gradients are in fact due to rotation .    in order to check whether planes are
good fits to the first - moment velocity maps , we make position - velocity cuts parallel to the maximum gradient directions and plot the results ( bottom panels of figures [ fig : oriona ]  [ fig : rosette ] ) . we are essentially plotting the central velocity at a given location in the cloud versus displacement along the gradient on a pixel by pixel basis .
figures [ fig : oriona ] , [ fig : monr2 ] , and [ fig : rosette ] show that planes are good fits to the first - moment maps of @xmath0  surrounding orion a , monr2 , and the rosette because there is a clear linear trend in the gradient in the atomic hydrogen .
the velocity fields of the @xmath0  associated with perseus and ngc 2264 appear to have more complex structure . in the case of the ngc 2264 and monr2 molecular clouds ( figures [ fig : ngc2264 ] and [ fig : monr2 ] ) , however , the position - velocity plots do not have monotonically increasing or decreasing slopes , indicating a more complex velocity structure in the _ molecular _ gas ( red points ) . in these two cases , we nevertheless overplot the lines perpendicular to the gradient direction calculated from equation [ eq : theta ] .
the specific angular momenta in the atomic and molecular gas , listed in table 3 , are compared in figure [ fig : jplot ] .
there appears to be a reasonable correlation@xmath93 and @xmath94 increasing together as @xmath95although small number statistics prevent us from making a firm conclusion . in each case , the initially expected specific angular momentum , @xmath93 , is always greater than @xmath94 .
these measurements alone are consistent with a picture whereby gmcs form via some top - down mechanism , such as a gravitational instability , and somehow shed angular momentum in the process .
but since angular momentum is a vector quantity , this scenario is difficult to reconcile with the observation that the gradient position angles in the molecular and atomic gas differ and appear to be uncorrelated ( figure [ fig : theta ] ) .
the key finding of our analysis is that the regions of atomic gas associated with molecular clouds have linear velocity gradients , yet the directions of these gradients are  with one exception  unaligned with the direction of the gradients in the associated gmcs . under the hypothesis that the gradients are caused by solid body rotation
, this would imply that gmcs are _ not _ corotating with the surrounding ism .
the second key result is that the magnitudes of the velocity gradients in the gmcs are larger than the gradient magnitudes in the atomic gas .
thirdly , if the gradients in the molecular and atomic gas are from rotation , the angular momenta in the molecular clouds is less than predicted from calculations of the angular momenta in the associated atomic gas .
below , we describe the results in detail for each gmc .
the gradient directions , magnitudes , and specific angular momenta of both the molecular and atomic gas are given in table 3 .
* *    figure [ fig : oriona ] shows that the orion a molecular cloud has a large - scale gradient whose velocity decreases from about 12  to 3  with increasing galactic longitude . of the gmcs in this study , observations of the kinematics and morphology of the orion a molecular cloud seem to make the best case for a top - down picture of gmc formation .
the direction of the gradient and the long axis of orion a are parallel to each other and to the galactic plane . if the gradient were due to rotation , this would imply that orion a is rotating in a sense nearly opposite to that of the galactic disk .
this result is in agreement with previous studies ( e.g. , kutner et al .
1977 , blitz 1993 ) .
the velocity gradient of the @xmath0  in the immediate surroundings of orion a points in nearly the same direction as that in the molecular cloud .
figure [ fig : oriona ] shows that the gradient in the atomic hydrogen , integrated over the velocity range from about @xmath96 to 22 , differs from th e gradient in the molecular cloud by only @xmath97 .
the origin of the velocity gradient in orion a has been debated ; it has previously been explained by cloud rotation ( kutner et al .
1977 ) , expansion driven by the orion ob association ( bally et al . 1987 ) , and expansion driven by stellar winds from newborn stars ( heyer et al .
if due to rotation , the 0.22  gradient in the molecular cloud implies a specific angular momentum of 42.2 .
however , considering that orion a is filamentary and much more closely resembles a cylinder rotating about its minor axis than a sphere , this estimate of @xmath69 is likely to be low limit , since the configuration of the former has a higher moment of inertia than the latter . for a cylinder rotating about its minor axis , @xmath98 , where @xmath71 is the radius of the cylinder and @xmath99 is the length . the cylindrical radius and length of orion a are approximately @xmath100 and @xmath97 , corresponding to 11 pc and 65 pc .
consequently , @xmath101 , much c loser to , but still a bit less than the estimated angular momentum of the associated atomic gas ( 107 ) .
these two results  that the gradients of the molecular and local atomic gas are in near alignment , and the magnitudes of the specific angular momenta are within range of each other  suggest that orion a may be a case in which the molecular cloud and surrounding @xmath0  are corotating and in which conservation of angular momentum is demonstrated . as we will see below
, it appears to be an exceptional case .    * *    figure [ fig : perseus ] shows that a strong linear gradient exists in the perseus molecular cloud .
the velocities in the gmc range from about 1 to 11  with increasing galactic latitude .
the gradient direction is tilted about @xmath102 to the north - west . of the molecular clouds in the sample
, perseus has the largest velocity gradient at 0.23 .
if the gradient is caused by rotation , from equation [ eq : j ] , the specific angular momentum in perseus is 24.2 , as indicated in table 3 .
examination of the black points in the position velocity - plot shown in figure [ fig : perseus ] indicates that the velocity gradient of the @xmath0  region centered on perseus is not linear . as with each cloud in this study
, we take the center of the @xmath0  field to be located at the position centroid of the gmc in galactic coordinates , ( @xmath103 ) .
this affects the appearance of the position - velocity plot , which shows the velocity of points in the field as a function of distance from the gradient judged from @xmath103 .
the @xmath0  surface density map in figure [ fig : perseus0 ] shows a high - surface density ( @xmath104 @xmath43 ) @xmath0  filament extending vertically near a galactic longitude of @xmath105 .
whereas the other gmcs in this sample have some portion of their molecular material laying directly on top of an @xmath0  peak , this is not strictly the case with perseus , which is tilted toward the east from the bright @xmath0  filament .
thus , we generate another first - moment map of the @xmath0  ce ntered on the filament at @xmath106 , as well as another position - velocity plot centered here .
figure [ fig : perseus2 ] shows that a linear fit is a much more suitable fit to the gradient in the @xmath0  when we shift the reference point . when measured this way , the magnitude of the gradient measured using equation [ eq : omega ] shifts from 0.038 to 0.067 , the value we cite in table 3 .
the gradient direction is nearly @xmath107 , very close to the sense of galactic rotation , in which velocities increase from west to east .
but the directions of the gradients in the molecular cloud and the atomic gas surrounding the gmc differ by nearly @xmath108 .    of the clouds in this analysis , perseus has the largest velocity gradients in the molecular gas .
the magnitude of the gradient in the molecular gas ( 0.23 ) is 3.4 times that in the @xmath0  ( @xmath109 ) .
the estimated specific angular momentum in the atomic gas is @xmath110 , 2.5 times @xmath69 measured in the molecular cloud .
although this is consistent with a scenario in which the perseus molecular cloud formed by collapsing out of the atomic gas , its angular momentum being redistributed somehow in the process , based on the observation that there is such a large difference between the gradient directions in the molecular and atomic gas , it is difficult to see how the gmc could have formed in such a simple way .    * *    the first - moment maps and position - velocity plot in figure [ fig : ngc2264 ] show that the ngc 2264 molecular cloud and surrounding @xmath0  both have complex kinematic features . though there is no significant linear gradient across the entire gmc , there is a weak gradient in the atomic gas .
the gradient in the atomic gas is stronger to the south of the rotation axis , where the velocity decreases from about 8 to 5  as the galactic latitude increases from @xmath111 to @xmath112 .
as seen in figure [ fig : ngc22640 ] , the ngc 2264 molecular cloud is composed of two main structures .
the structure in the north has a larger range of velocities , with @xmath113  10 , than does the structure in the south , which averages around 6 , close to the mean velocity of the molecular cloud as a whole .
this is also the average velocity of the @xmath0  peak located at near @xmath114 , @xmath115 ( figure [ fig : ngc22640 ] ) , which is associated with the southern segment of t he gmc .
the magnitude of the gradient in the atomic gas near ngc 2264 is 0.019 , which is less than @xmath57(@xmath116 ) in the solar vicinity assuming a flat rotation curve .
the magnitude of the gradient fitted across the entire gmc , using equation [ eq : omega ] , is 0.046 .
but the northern segment of the cloud in isolation has a gradient closer to 0.08 .
it seems unlikely that rotation is the origin of the velocity field in the ngc 2264 molecular cloud .
table 3 records the velocity magnitude and direction calculated by fitting a plane to the field , even though the gradient is not linear . using these values , the specific angular momentum of the gmc , @xmath117 , is smaller than the expected value by at least a factor of 6 . the gradient direction in the gmc , pointing nearly @xmath118 east from north , differs from @xmath90 by @xmath119 .
an alternative explanation of the velocity field in ngc 2264 may be , at least in part , the internal stellar activity .
it is noteworthy that the @xmath120  region associated with ngc 2264 ( the cone nebula ) is located in the north , near the part of the gmc that has a larger velocity dispersion than the southern part .
the nebula , at @xmath121 , @xmath122 ( kharchenko et al .
2005 ) is located between two high velocity regions of the gmc that are receding at speeds near 10 .
this morpholog y is suggestive of an expanding ring seen from an edge - on perspective . taken all together , it appears that these features could be causally connected : an @xmath120  region causes the high - speed expansion of the surrounding molecular gas , which subsequently sweeps up atomic gas into a high - density ridge .
it is perhaps significant that we see a similar pattern in the rosette .
yet another explanation of the kinematics of ngc 2264 is put forth by fursz et al .
( 2006 ) , based on their finding that the stars and emission in ngc 2264 are well - correlated in position - velocity space .
they suggest that the velocity field of the gmc is explained by the models of burkert & hartmann ( 2004 ) , in which molecular clouds form from supersonic collisions of gas .
this argument seems to be corroborated by the overall pattern of the cloud s velocity structure : higher velocities at the outskirts of the cloud and lower velocities proceeding toward the center of the cloud .
this pattern could be explained by gravitationally driven infall at the interface of the colliding flows . in this scenario ,
star formation occurs preferentially in the condensations that develop during the collapse of the molecular cloud , thus explaining the position - velocity correlation observed between stars and high - density molecular gas .
the colliding flows hypothesis may also explain the velocity structure of the atomic gas .
figure [ fig : ngc2264 ] shows that the @xmath0  has slightly higher velocities at positions in the field far from the gradient axis and lower velocities close to the axis  again , suggestive of infall accelerated by gravity . if other molecular clouds form by this mechanism , however , it is unclear why we do not see this pattern in the atomic gas surrounding the rest of the gmcs in our sample
. possibly , the ngc 2264 system still retains signs of its early formation history due to the relatively young age of the gmc ( @xmath123  3 myr ; flaccomio et al .
2000 ; ramrez et al .
it is difficult to make a definitive conclusion since we are observing @xmath0  so close to the galactic plane where the problem of line - of - sight blending is exacerbated .    * *    figures [ fig : monr20 ] and [ fig : monr2 ] show that the monr2 molecular cloud sits right on top of a high - density cloud of @xmath0  that has a strong linear velocity gradient .
the gradient , which extends for over 80 pc , has a magnitude of 0.035  and increases from about 7  to 12  from south to north - west .
if the gradient is due to rotation , @xmath124  180 , depending on the size of the region considered . from conservation of angular momentum , one might expect monr2 to have a gradient somewhere between @xmath125  0.84 , pointed about @xmath126 west of north like the gradient in the atomic gas .
this is not what we see in figure [ fig : monr2 ] , however .
the velocity field of the monr2 molecular cloud has a much more complex structure than does the atomic gas .
the outer northern , western , and southern edges of the cloud reach velocities up to 12  13 .
the inner portion of the cloud is moving at slower speeds around 9  11 , which is also the velocity range at which the @xmath0  peak located at @xmath127 , @xmath128 dominates . considering the apparently random nature of the velocity field , it is unlikely that the gmc is undergoing large - scale , coherent rotation .
the position - velocity plot for the molecular cloud in figure [ fig : monr2 ] ( red points ) has a zero slope , as does ngc 2264 , indicating that it has no significant linear gradient .
it is difficult to say whether or not the present velocity field originated during the cloud s formation or during some later stage in its evolution .
if the gmc originally had a more organized velocity field and did _ not _ inherent its present field during formation , perhaps a series of interactions with external forces or internal events  such as turbulence and star formation activity  have washed out any signature of a systematic velocity gradient which was previously present in the molecular cloud .    * *    the separation between the gradient directions in the rosette molecular cloud and the surrounding atomic gas is approximately @xmath1 .
figure [ fig : rosette ] shows that in the velocities of the molecular cloud tend to increase from roughly 7
14  in the north - east to as high as @xmath129  in the south and south - west .
the gradient in the @xmath0  is directed perpendicular to the galactic plane .
we note that , of the clouds in our sample , the @xmath0  maps near the rosette probably suffer the greatest degree of line - of - sight blending .
the cloud is very close to the galactic plane , and at the distance of the rosette , 1600 pc , we have the lowest spatial resolution ( @xmath129 pc ) in the atomic gas .
therefore the surface density and velocity maps displayed in figures [ fig : rosette0 ] and [ fig : rosette ] almost certainly fail to capture many of the local , small - scale variations in the structure of the atomic gas .
we measure a gradient in the molecular cloud of 0.09 , consistent with the 0.08  measured by williams et al .
if the gradient is due to rotation , the specific angular momentum of the cloud is approximately 26 .
based on the estimate of the minimum value of the accumulation radius , the specific angular momentum in the surrounding ism is nearly 3 times larger .
figure [ fig : rosette0 ] shows the outline of the molecular cloud overlaid on a surface density map of the atomic gas , which is integrated over the range @xmath130  27 .
a `` shell '' of @xmath0  containing two high - density peaks is associated with the molecular cloud .
the brightest peak of the shell sits at the southern edge of the molecular cloud near a latitude of @xmath131 , and another peak in the surface density occurs near @xmath132 .
the southern portion of the @xmath0  shell appears to mimic the ring - like structure in the molecular cloud sitting just to the north of it .
when the location of the structures in both the molecular and atomic gas are compared with the first - moment map in figure [ fig : rosette ] , we see that the southern portion of the @xmath0  shell is moving near the same velocity , @xmath133 - 14 , as the south - east segment of the ring in the gmc .
the western half of the ring in the gmc is moving at higher velocities .
also note how the @xmath0  peak locate d near @xmath132 is moving slightly faster ( @xmath134 ) than the molecular gas in that region ( @xmath133 ) .
all of this suggests that both the @xmath0  shell and the ring - like structure in the gmc are expanding .
kuchar & bania ( 1993 ) demonstrated that the @xmath120  region ngc 2244 ( the rosette nebula ) could have given rise to the expansion in the atomic gas . and
williams et al .
( 1995 ) found that certain properties of clumps within the gmc vary with distance from the nebula , centered at @xmath135 , @xmath136 ( celnik 1983 ) .
taken all together , the evidence leads us to suggest that the gradient in the rosette molecular cloud is _ not _ caused by rotation , but by the high - luminosity @xmath120  region , ngc 2244 .
high - energy winds of stars in ngc 2244 may have excavated a hole in the gmc , causing @xmath0  to be swept up into a high - density ridge by the expanding molecular gas ( e.g. , kuchar et al .
1991 , kuchar & bania 1993 ) . in this picture , because the nebula has had less of an impact on distant regions of the cloud , these distant regions are moving at lower velocities compared to molecular gas near the @xmath120  region .
our key findings are that the angular momentum in the gmcs is less than that in the surrounding atomic gas , and the velocity gradient position angles in the molecular and atomic gas are widely divergent  with orion a being the one exception .
this leads us to suggest that rotation may _ not _ be the best explanation of the velocity fields observed in the gmcs .    traditionally , at least three possible solutions have been invoked to resolve the angular momentum problem :    \1 .
one or more of the assumptions in the theory are inappropriate .
for instance , if the average surface density of the precursor atomic gas is underestimated , this will lead to an overestimate in the accumulation radius and , consequently , an overestimate of the angular momentum initially imparted to a gmc . in our analysis , we used 10  for the mean surface density of the precursor gas , twice the average value in the solar vicinity . assuming that the initial gas had an even higher density ,
say 20 , this would not change the main result , namely that @xmath93 is consistently greater than @xmath94 .
this is because the accumulation radius depends weakly on the initial surface density of the gas ( @xmath137 ) . on the other hand ,
if collapsing molecular clouds do not gather material far from the galactic plane , using larger values of @xmath76 might be appropriate . yet
this would only exacerbate the discrepancy between the predicted and obser ved angular momenta .
it might also be argued that using the effective radius of a gmc leads to underestimates in @xmath138 if the gmc is filamentary . in the previous section we took this into consideration with orion a and recalculated @xmath138 assuming a cylindrical morphology .
this raised @xmath138 , which , nevertheless , remained less than @xmath93 .
moreover , varying the size of the gmc or the region from which it gathers material does not solve the problem of the gradient directions in the molecular and atomic gas being unaligned .
an alternative explanation of the angular momentum problem is that though the assumptions regarding the initial conditions may be valid , there may be some kind of external braking force which rapidly reduces the angular momentum of a gmc during its initial condensation .
magnetic braking is often evoked as an angular momentum shedding mechanism ( e.g. , mouschovias 1977 ; fleck & clark 1981 ; mestel & paris 1984 ; rosolowsky et al
zeeman splitting of the oh 18-cm line and the 21-cm line of neutral hydrogen have been used to measure the magnetic field strengths of gmcs .
since gmcs are magnetized and mhd effects are expected to play a significant role in their evolution , magnetic braking  in which magnetic field lines anchoring a gmc to the ambient ism provide the tension necessary to slow down rotation
is a possible solution to the angular momentum problem .
heiles & troland ( 2005 ) measured the mean magnetic field strength in the cold neutral medium of the milky way to be @xmath139 .
the braking time is set by the time it takes for alfvn waves to travel across a region of gas having a moment of inertia comparable to the gmc ( mestel & paris 1984 ) . in order for magnetic braking to be efficient at slowing down cloud rotation
, it would have to occur on timescales no greater than the timescale for cloud collapse . for gas having an initial density of 1 @xmath82 ( = 0.025 @xmath43  @xmath140 ) , the alfvn speed is @xmath141 . in a region
having an accumulation radius of @xmath142 pc , this corresponds to a braking time of roughly @xmath143 myr . by comparison , the timescale for self - gravitational cloud formation in a region having the same density , assuming that this process occurs on a timescale close to the dynamical free - fall time ( mouschovias & paleologou 1979 ; mestel & paris 1984 ; elmegreen 2007 ) , is @xmath144 myr .
however , the effectiveness of magnetic braking largely depends on how the mass of the forming cloud compares to its magnetic critical mass .
mestel & paris ( 1984 ) argue that braking will efficiently slow down rotation only if the mass , @xmath145 , of the forming cloud is much less than its magnetic critical mass , @xmath146 . using crutcher s ( 1999 ) magnetic field strength measurements of molecular clouds ,
mckee & ostriker ( 2007 ) infer that gmcs have approximately @xmath147 , that is , they are magnetically supercritical . elmegreen ( 2007 ) argues that the ism out of which gmcs form has a magnetic field that is `` near - critical , '' i.e. , @xmath148 .
furthermore , magnetic braking slows rotation most efficiently when the cloud s angular momentum vector is perpendicular to the magnetic field ( mouschovias & paleologou 1979 ) . due to the complexity of making precise measurements of the magnetic field in the ism
, it is not established whether or not this is the case .
clearly , more observatio ns are needed in order to determine the importance of magnetic fields in gmc evolution .
finally , perhaps gmcs are not rotating after all , and the `` problem '' is unfounded . as some authors have pointed out , the interpretation that velocity gradients indicate rotation is not unique .
expansion and shear , for instance , also produce velocity gradients .
the numerical simulations of hennebelle et al .
( 2008 ) suggest that the converging flows of atomic gas could produce gmcs .
if shearing occurs at the interface of colliding flows , might this result in an excess of shear , that is , values of @xmath57 that are higher than the shear arising from the local galactic rotation curve ? assuming a flat rotation curve , the local value of @xmath57 in the solar neighborhood is @xmath149 .
this is slightly less than , but basically comparable to the gradient magnitudes we measure in perseus , orion a , and monr2 ( table 3 ) .
as discussed previously , we are underestimating @xmath57 in every case , since we have no information regarding the inclination of clouds .
also , since bl ending of emission along the line - of - sight is most extreme in the cases of ngc 2264 and the rosette , which are situated close to the galactic plane , it is likely that @xmath150 is underestimated to an even greater degree in these cases .
nevertheless , we can not conclude whether or not shear is the cause of the velocity gradients based on these observations .
it is an issue we would like to further investigate with a larger sample of clouds .
burhert & bodenheimer ( 2000 ) show that turbulence may also cause linear velocity gradients .
they demonstrated that the gradient magnitude of turbulent cores scales with size as @xmath151 .
figure [ fig : r_vgrad ] shows @xmath57 as a function of @xmath71 for both the molecular clouds and @xmath0 .
neither the gmcs nor the @xmath0  appear to follow the @xmath151 relationship , though we do not have enough data points to make a definite conslusion .
studying the properties of gmc kinematics in other galaxies will also help to develop our picture of gmc formation . in their analysis of 45 gmcs in m33 , rosolowsky et al .
( 2003 ) did the first systematic , extragalactic study of gmc angular momentum properties .
they showed that simple gmc formation theories consistently overestimate the magnitude of the observed angular momentum .
they measured the velocity gradients from high - resolution ( @xmath152 ) data , finding that the gradients of m33 clouds are similar in magnitude to galactic clouds .
they then tested several formation models by calculating the accumulation radii of the cataloged clouds which are predicted by the various models , including the toomre and parker instabilities . on average , the theories , which do not include the effects of magnetic fields on rotation , overpredict observed magnitudes of velocity gradients by more than a factor of 5 .
and 40% of the gmcs in m33 are counterrotating with respe ct to the sense of galactic rotation .
we extend this study by performing the analysis established here on a much larger sample of gmcs in m33 cataloged by rosolowsky et al .
since m33 has a relatively low inclination of @xmath153 ( corbelli & salucci 2000 ) , such a study has the advantage of bypassing the problem of source confusion along lines of sight that arises when doing comparable surveys of molecular clouds in the milky way .
we have presented a detailed comparison between the kinematics of five galactic molecular clouds and the atomic gas that surrounds them .
we developed a method for selecting regions of @xmath0  that are associated with the gmcs and found that each gmc was nearby high - density @xmath0  which peaked above the mean galactic value .
first - moment maps were created using the emission of the gmcs and the @xmath0  21-cm emission , and then a plane was fit to each map of velocity centroids .
we determined the magnitudes and directions of the velocity gradients from the coefficients to the fits . from these observations and measurements , we arrived at the following conclusions :    \1 .
orion a , perseus , and the rosette each have a significant linear velocity gradient across the face of the molecular clouds , while ngc 2264 and monr2 have complex , non - linear velocity fields .
the perseus molecular cloud has the strongest linear gradient in the sample , with a magnitude of @xmath154 .
the atomic gas associated with orion a , monr2 and the rosette has significant linear velocity gradients , regardless of whether the molecular cloud has one .
the @xmath0  gradients range from 0.019 to 0.067 , or 0.76 to 2.7 times the shear in the solar vicinity as measured by the galactic rotation curve .
if the gradients in the molecular and atomic gas were due to rotation , the specific angular momentum in the gmcs is less than that predicted by the formation scenario in which a gmc preserves angular momentum while undergoing top - down collapse from the surrounding ism by a factor of 1 to 6 . the discrepancy can be narrowed if different assumptions are made regarding the initial density and geometry of the gas . but the direction of the trend  that the observed angular momentum is less than the predicted  remains the same .
we observe large differences between the velocity gradient directions in the gmcs and the atomic gas , with orion a being the one exception . at more than @xmath1 ,
two of the most extreme angle separations can be seen in the monr2 and the rosette systems ( figure [ fig : monr2 ] and [ fig : rosette ] ) .
furthermore , the gradient directions in neither the molecular nor the atomic gas are in alignment with the overall direction of galactic rotation .
if the gradients were due to rotation , this indicates that some gmcs are counterrotating with respect to the galaxy .    \5 .
that the velocity gradient position angles in the atomic and molecular gas are divergent indicates that the gmcs in our sample probably did not inherit their present velocity fields from the atomic gas from which they formed .
finally , in at least two cases , ngc 2264 and the rosette , a good explanation of the morphology and kinematics observed in the gas is that they are caused by stellar winds from o stars in the @xmath120  regions located in the gmcs , not rotation .
lcccccccccc +   + & & & @xmath37 & @xmath77 & & @xmath155 & @xmath156 & + cloud & @xmath157 & @xmath158 & [ pc ] & [ @xmath159 & [ @xmath159 & [ pc@xmath160 & [ ] & [ ] & & + perseus ....... & 159 & -20 & 320 & 1.4 & 1.2 & 825 &  7.5 & @xmath161 & 1 & 6 + orion a ....... & 210 & -19 & 414 & 2.4 & 6.9 & 1510 &  9.6 & @xmath162 & 2 & 7 + ngc 2264 ... & 202 & 1 & 800 & 2.5 & 2.2 & 674 &  6.6 & @xmath163 & 3 & 8 + monr2 ........ & 213 & -13 & 830 & 3.0 & 9.0 & 2940 & 10.8 & @xmath164 & 4 & 9 + rosette ........ & 207 & -2 & 1600 & 4.0 & 7.7 & 2200 & 13.1 & @xmath165 & 5 & 10 +      lcccccc +   + cloud & @xmath167 & @xmath150 & @xmath89 & @xmath90 & @xmath168 & @xmath93 + & [ @xmath169 & [ @xmath169 & [ deg ] & [ deg ] & [ ] & [ ] + perseus ....... & 23.1 & 6.70 & -20.2 & 99.4 & 24.2 & 60.3 + orion a ....... & 22.0 & 3.71 & -111 & -102 & 42.2 [ 84.1 ] & 107 + ngc 2264 ... & 4.60 & 1.87 & 47.8 & -173 &  3.92 & 25.2 + monr2 ........ & 6.68 & 3.48 & -72.4 & -11.2 & 25.5 & 120 + rosette ........ & 9.18 & 2.36 & -132 &  1.12 & 25.8 & 73.4 +