No Arabic abstract
Column density (N) PDFs serve as a powerful tool to characterize the physical processes that influence the structure of molecular clouds. Star-forming clouds can best be characterized by lognormal PDFs for the lower N range and a power-law tail for higher N, commonly attributed to turbulence and self-gravity and/or pressure, respectively. We report here on PDFs obtained from observations of 12CO, 13CO, C18O, CS, and N2H+ in the Cygnus X North region and compare to a PDF derived from dust observations with the Herschel satellite. The PDF of 12CO is lognormal for Av~1-30, but is cut for higher Av due to optical depth effects. The PDFs of C18O and 13CO are mostly lognormal up for Av~1-15, followed by excess up to Av~40. Above that value, all CO PDFs drop, most likely due to depletion. The high density tracers CS and N2H+ exhibit only a power law distribution between Av~15 and 400, respectively. The PDF from dust is lognormal for Av~3-15 and has a power-law tail up to Av~500. Absolute values for the molecular line column densities are, however, rather uncertain due to abundance and excitation temperature variations. Taken the dust PDF face value, we calibrate the molecular line PDF of CS to the one of the dust and determined an abundance [CS]/[H2] of 10^-9. The slopes of the power-law tails of the CS, N2H+, and dust PDFs are consistent with free-fall collapse of filaments and clumps. A quasi static configuration of filaments and clumps can possibly also account for the observed N-PDFs, as long as they have a sufficiently condensed density structure and external ram pressure by gas accretion is provided. The somehow flatter slopes of N2H+ and CS can reflect an abundance change and/or subthermal excitation at low column densities.
We analyse column density and temperature maps derived from Herschel dust continuum observations of a sample of massive infrared dark clouds (G11.11-0.12, G18.82-0.28, G28.37+0.07, G28.53-0.25). We disentangle the velocity structure of the clouds using 13CO 1-0 and 12CO 3-2 data, showing that these IRDCs are the densest regions in massive giant molecular clouds and not isolated features. The probability distribution function (PDF) of column densities for all clouds have a power-law distribution over all (high) column densities, regardless of the evolutionary stage of the cloud: G11.11-0.12, G18.82-0.28, and G28.37+0.07 contain (proto)-stars, while G28.53-0.25 shows no signs of star formation. This is in contrast to the purely log-normal PDFs reported for near/mid-IR extinction maps. We only find a log-normal distribution for lower column densities, if we perform PDFs of the column density maps of the whole GMC in which the IRDCs are embedded. By comparing the PDF slope and the radial column density profile, we attribute the power law to the effect of large-scale gravitational collapse and to local free-fall collapse of pre- and protostellar cores. Independent from the PDF analysis, we find infall signatures in the spectral profiles of 12CO for G28.37+0.07 and G11.11-0.12, supporting the scenario of gravitational collapse. IRDCs are the densest regions within GMCs, which may be the progenitors of massive stars or clusters. At least some of the IRDCs are probably the same features as ridges (high column density regions with N>1e23 cm-2 over small areas), which were defined for nearby IR-bright GMCs. Because IRDCs are only confined to the densest (gravity dominated) cloud regions, the PDF constructed from this kind of a clipped image does not represent the (turbulence dominated) low column density regime of the cloud.
We test some ideas for star formation relations against data on local molecular clouds. On a cloud by cloud basis, the relation between the surface density of star formation rate and surface density of gas divided by a free-fall time, calculated from the mean cloud density, shows no significant correlation. If a crossing time is substituted for the free-fall time, there is even less correlation. Within a cloud, the star formation rate volume and surface densities increase rapidly with the corresponding gas densities, faster than predicted by models using the free-fall time defined from the local density. A model in which the star formation rate depends linearly on the mass of gas above a visual extinction of 8 mag describes the data on these clouds, with very low dispersion. The data on regions of very massive star formation, with improved star formation rates based on free-free emission from ionized gas, also agree with this linear relation.
We characterize the column density probability distributions functions (PDFs) of the atomic hydrogen gas, HI, associated with seven Galactic molecular clouds (MCs). We use 21 cm observations from the Leiden/Argentine/Bonn Galactic HI Survey to derive column density maps and PDFs. We find that the peaks of the HI PDFs occur at column densities ranging from ~1-2$times 10^{21}$ cm$^2$ (equivalently, ~0.5-1 mag). The PDFs are uniformly narrow, with a mean dispersion of $sigma_{HI}approx 10^{20}$ cm$^2$ (~0.1 mag). We also investigate the HI-to-H$_2$ transition towards the cloud complexes and estimate HI surface densities ranging from 7-16 $M_odot$ pc$^{-2}$ at the transition. We propose that the HI PDF is a fitting tool for identifying the HI-to-H$_2$ transition column in Galactic MCs.
We investigate Schmidts conjecture (i.e., that the star formation rate scales in a power-law fashion with the gas density) for four well-studied local molecular clouds (GMCs). Using the Bayesian methodology we show that a local Schmidt scaling relation of the form Sigma*(A_K) = kappa x (A_K)^{beta} (protostars pc^{-2}) exists within (but not between) GMCs. Further we find that the Schmidt scaling law, by itself, does not provide an adequate description of star formation activity in GMCs. Because the total number of protostars produced by a cloud is given by the product of Sigma*(A_K) and S(> A_K), the differential surface area distribution function, integrated over the entire cloud, the clouds structure plays a fundamental role in setting the level of its star formation activity. For clouds with similar functional forms of Sigma*(A_K), observed differences in their total SFRs are primarily due to the differences in S(> A_K) between the clouds. The coupling of Sigma*(A_K) with the measured S(> A_K) in these clouds also produces a steep jump in the SFR and protostellar production above A_K ~ 0.8 magnitudes. Finally, we show that there is no global Schmidt law that relates the star formation rate and gas mass surface densities between GMCs. Consequently, the observed Kennicutt-Schmidt scaling relation for disk galaxies is likely an artifact of unresolved measurements of GMCs and not a result of any underlying physical law of star formation characterizing the molecular gas.
Supernovae from core-collapse of massive stars drive shocks into the molecular clouds from which the stars formed. Such shocks affect future star formation from the molecular clouds, and the fast-moving, dense gas with compressed magnetic fields is associated with enhanced cosmic rays. This paper presents new theoretical modeling, using the Paris-Durham shock model, and new observations, using the Stratospheric Observatory for Infrared Astronomy (SOFIA), of the H$_2$ S(5) pure rotational line from molecular shocks in the supernova remnant IC443. We generate MHD models for non-steady-state shocks driven by the pressure of the IC443 blast wave into gas of densities $10^3$ to $10^5$ cm$^{-3}$. We present the first detailed derivation of the shape of the velocity profile for emission from H$_2$ lines behind such shocks, taking into account the shock age, preshock density, and magnetic field. For preshock densities $10^3$-$10^5$ cm$^{-3}$, the the predicted shifts of line centers, and the line widths, of the H$_2$ lines range from 20-2, and 30-4 km/s, respectively. The a priori models are compared to the observed line profiles, showing that clumps C and G can be explained by shocks into gas with density 10$^3$ to $2times 10^4$ cm$^{-3}$ and strong magnetic fields. For clump B2 (a fainter region near clump B), the H$_2$ spectrum requires a J-type shock into moderate density (~100 cm$^{-3}$) with the gas accelerated to 100 km/s from its pre-shock location. Clump B1 requires both a magnetic-dominated C-type shock (like for clumps C and G) and a J-type shock (like for clump B1) to explain the highest observed velocities. The J-type shocks that produce high-velocity molecules may be locations where the magnetic field is nearly parallel to the shock velocity, which makes it impossible for a C-type shock (with ions and neutrals separated) to form.