Turbulent molecular clouds

Stars form within molecular clouds but our understanding of this fundamental process remains hampered by the complexity of the physics that drives their evolution. We review our observational and theoretical knowledge of molecular clouds trying to confront the two approaches wherever possible. After a broad presentation of the cold interstellar medium and molecular clouds, we emphasize the dynamical processes with special focus to turbulence and its impact on cloud evolution. We then review our knowledge of the velocity, density and magnetic fields. We end by openings towards new chemistry models and the links between molecular cloud structure and star-formation rates.     

Large-Scale Gravitational Instability in the Formation of GMCS

Two models of molecular cloud in disk galaxies are proposed to investigate the formation of giant molecular clouds (GMCs) under the gravitational instability and random collision using PP(Particle–Particle) simulation. Having analysed simulation outputs of the two models and comparing them with observation, we are able to draw some general conclusions, the most significant ones of which are: 1) Similar to results obtained previously, the gravitational instability can make small clouds form large clouds faster than random collision. 2) The differential rotation in gravitational instability model plays a positive role in agglomeration of molecular clouds. This revised version was published online in July 2006 with corrections to the Cover Date.     

What Drives Star Formation?

Current theoretical models for what drives star formation (especially low-mass star formation) are: (1) magnetic support of self-gravitating clouds with ambipolar diffusion removing support in cores and triggering collapse and (2) compressible turbulence forming self-gravitating clumps that collapse as soon as the turbulent cascade produces insufficient turbulent support. Observations of magnetic fields can distinguish between these two models because of different predictions in three areas: (1) magnetic field morphology, (2) the scaling of field strength with density and non-thermal velocities, and (3) the mass to magnetic flux ratio, M/Φ. We first discuss the techniques and limitations of methods for observing magnetic fields in star formation regions, then describe results for the L1544 prestellar core as an exemplar of the observational results. Application of the three tests leads to the following conclusions. The observational data show that both magnetic fields and turbulence are important in molecular cloud physics. Field lines are generally regular rather than chaotic, implying strong field strengths. But fields are not aligned with the minor axes of oblate spheroidal clouds, suggesting that turbulence is important. Field strengths appear to scale with non-thermal velocity widths, suggesting a significant turbulent support of clouds. Giant Molecular Clouds (GMCs) require mass accumulation over sufficiently large volumes that they would likely have an approximately critical M/Φ. Yet H I clouds are observed to be highly subcritical. If self-gravitating (molecular) clouds form with the subcritical M/Φ of H I clouds, the molecular clouds will be subcritical. However, the observations of molecular cloud cores suggest that they are approximately critical, with no direct evidence for subcritical molecular clouds or cloud envelopes. Hence, the observations remain inconclusive in deciding between the two extreme-case models of what drives star formation. What is needed to further advance our understanding of the role of magnetic fields in the star formation process are additional high sensitivity surveys of magnetic field strengths and other cloud properties in order to further refine the assessment of the importance of magnetic fields in molecular cores and envelopes.     

Collapse of primordial clouds — II. The role of dark matter

In this article we extend the study performed in our previous article of the collapse of primordial objects. We here analyse the behaviour of the physical parameters for clouds ranging from 10⁷ to 10¹⁵ M_⊙. We study the dynamical evolution of these clouds in two ways: as purely baryonic clouds and as clouds with non-baryonic dark matter included. We start the calculations at the beginning of the recombination era, following the evolution of the structure until the collapse (which we defined as the time when the density contrast of the baryonic matter is greater than 10⁴). We analyse the behaviour of several physical parameters of the clouds (e.g. the density contrast and the velocities of the baryonic matter and the dark matter) as a function of time and radial position in the cloud. In this study all physical processes that are relevant to the dynamical evolution of the primordial clouds, such as for example photon drag (due to the cosmic background radiation) and hydrogen molecular production, besides the expansion of the Universe, are included in the calculations. In particular we find that the clouds with dark matter collapse at higher redshift when we compare the results with the purely baryonic models. As a general result we find that the distribution of the non-baryonic dark matter is more concentrated than the baryonic one. It is important to stress that we do not take into account the putative virialization of the non-baryonic dark matter; we just follow the time and spatial evolution of the cloud, solving its hydrodynamical equations. We also studied the role of cooling–heating processes in the purely baryonic clouds.     

Hypersonic Swizzle Sticks: Protostellar Turbulence, Outflows and Fossil Outflow Cavities

The expected lifetimes for molecular clouds has become a topic of considerable debate as numerical simulations have shown that MHD turbulence, the nominal means of support for clouds against self-gravity, will decay on short timescales. Thus it appears that either molecular clouds are transient features or they are resupplied with turbulent energy through some means. Jets and molecular outflows are recognized as a ubiquitous phenomena associated with star formation. Stars however form not isolation but in clusters of different density and composion. The ubiquity and high density of outflows from young stars in clusters make them an intriguing candidate for the source of turbulence energy in molecular clouds. In this contribution we present new studies, both observational and theoretical, which address the issue of jet/outflow interactions and their abilityto drive turbulent flows in molecular clouds. Our studies focus on scales associated with young star forming clusters. In particular we first show that direct collisions between active outflows are not effective at stirring the ambient medium. We then show that fossil cavities from “extinct” outflows may provide the missing link in terms of transferring momentum and energy to the cloud.     

行星际磁云研究新进展

从飞船的观测结果、磁云形态及演化的理论模型、磁流体动力学(MHD)数值模拟、激波对磁云的作用、多重磁云等5个方面，评述了行星际磁云的研究成果及最新进展。在太阳峰年，大部分的非重现性地磁暴都与磁云有关。最近的研究表明，压缩后的磁云往往能产生更大的地磁效应。深入研究磁云对空间天气研究有着特殊的价值，特别是对提高大磁暴的预报水平有着重要帮助。     

CO J=2-1 and J=3-2 Line Observations of Molecular Clouds toward the Directions of 59 EGOs in the Northern Sky

In order to investigate the differences between the molecular clouds which are associated with the massive star forming regions and those which are not, we have performed the single-dish simultaneous observations of ¹²CO J=2-1 and J=3-2 lines toward a sample of 59 Spitzer Extended Green Objects (EGOs) as the massive star forming regions in the northern sky. Combining our results with the data of the ¹²CO J=1-0 observations toward the same sample EGOs in the literature, we have made the statistical comparisons on the intensities and linewidths of multiple ¹²CO lines between the molecular clouds associated with EGOs (EGO molecular clouds, in brief) and other non-EGO molecular clouds. On this basis, we have discussed the effects of the gas temperature, density, and velocity field distributions on the statistical characteristics of the two kinds of molecular clouds. It is found that both the EGO molecular clouds and non-EGO molecular clouds have similar mass ranges, hence we conclude that for the formation of massive stars, the key-important factor is probably not the total mass of a giant molecular cloud (GMC), but the volume filling factor of the molecular clumps in the GMC (or the compression extent of the molecular gas in the cloud).     

Noctilucent clouds: Simulation studies of their genesis,properties and global influences

The formation, evolution and properties of noctilucent clouds are studied using a timedependent one-dimensional model of ice particles at mesospheric altitudes. The model treats ice crystals, meteoric dust, water vapor and air ionization as fully interactive cloud elements. For ice particles, the microphysical processes of nucleation, condensation, coagulation and sedimentation are included; the crystal habits of ice are also accounted for. Meteoric dust is analyzed in the manner of Hunten et al. (1980). The simulated particle sizes range from 10 Å to 2.6μm. The chemistry of water vapor and the charge balance of the mesosphere are also analyzed in detail.Based on model calculations, including numerous sensitivity tests, several conclusions are reached. Extremely cold mesopause temperatures (<140K) are necessary to form noctilucent clouds; such temperatures only exist at high latitudes in summer. A water vapor concentration of 4–5 ppmv is sufficient to form a visible cloud. However, a subvisible cloud can exist in the presence of only 1 ppmv of H₂O. Ample cloud condensation nuclei are always present in the mesosphere; at very low temperatures, either meteoric dust or hydrated ions can act as cloud nuclei. To be effective, meteoric dust particles must be larger than 10–15 Å in radius. When dust is present, water vapor supersaturations may be held to such low values that ion nucleation is not possible. Ion nucleation can occur, however, in the absence of dust or at extremely low temperatures (<130K). While dust nucleation leads to a small number (<10cm^?3) of large ice particles (>0.05 μm radius) and cloud optical depths (at 550 nm) ～10^?4, ion nucleation generally leads to a large number (～10³cm^?3) of smaller particles and optical depths ～10^?5). However, because calculated nucleation rates in noctilucent clouds are highly uncertain, the predominant nucleus for the clouds (i.e., dust or ions) cannot be unambiguously established. Noctilucent clouds require several hours-up to a day-to materialize. Once formed, they may persist for several days, depending on local meteorological conditions. However, the clouds can disappear suddenly if the air warms by 10–20 K. The environmental conditions which exist at the high-latitude summer mesopause, together with the microphysics of small ice crystals, dictate that particle sizes will be ? 0.1 μm radius. The ice crystals are probably cubic in structure. It is demonstrated that particles of this size and shape can explain the manifestations of noctilucent clouds. Denser clouds are favored by higher water vapor concentrations, more rapid vertical diffusion and persistent upward convection (which can occur at the summer pole). Noctilucent clouds may also condense in the cold “troughs” of gravity wave trains. Such clouds are bright when the particles remain in the troughs for several hours or more; otherwise they are weak or subvisible.Model simulations are compared with a wide variety of noctilucent cloud data. It is shown that the present physical model is consistent with most of the measurements, as well as many previous theoretical results. Ambient noctilucent clouds are found to have a negligible influence on the climate of Earth. Anthropogenic perturbations of the clouds that are forecast for the next few decades are also shown to have insignificant climatological implications.     

Star formation in the LMC: Comparative CCD observations of young stellar populations in two giant molecular clouds

This work deals with a CCD imaging study at optical and near‐infrared wavelength oftwo giant molecular clouds (plus a control field) in the southern region of the Large Magellanic Cloud, one ofwhich shows multiple signs of star formation, whereas the other does not. The observational data from VLT FORS2 (R band) and NTT SOFI (Ks band) have been analyzed to derive luminosity functions and color‐magnitude diagrams. The young stellar content of these two giant molecular clouds is compared and confirmed to be different, in the sense that the apparently “starless” cloud has so far formed only low‐luminosity, low‐mass stars (fainter than m_Ks ∽ 16.5 mag, not seen by 2MASS), while the other cloud has formed both faint low‐mass and luminous high‐mass stars. The surface density excess oflow‐luminosity stars (∽2 per square arcmin) in the “starless” cloud with respect to the control field is about 20% whereas the excess is about a factor of 3 in the known star‐forming cloud. The difference may be explained theoretically by the gravo‐turbulent evolution of giant molecular clouds, one being younger and less centrally concentrated than the other (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Desorption from interstellar ices

The desorption of molecular species from ice mantles back into the gas phase in molecular clouds results from a variety of very poorly understood processes. We have investigated three mechanisms: desorption resulting from H₂ formation on grains, direct cosmic ray heating and cosmic ray-induced photodesorption. Whilst qualitative differences exist between these processes (essentially deriving from the assumptions concerning the species selectivity of the desorption and the assumed threshold adsorption energies, E _t), all the three processes are found to be potentially very significant in dark cloud conditions. It is therefore important that all three mechanisms should be considered in studies of molecular clouds in which freeze-out and desorption are believed to be important.
Employing a chemical model of a typical static molecular core and using likely estimates for the quantum yields of the three processes, we find that desorption by H₂ formation probably dominates over the other two mechanisms. However, the physics of the desorption processes and the nature of the dust grains and ice mantles are very poorly constrained. We therefore conclude that the best approach is to set empirical constraints on the desorption, based on observed molecular depletions – rather than try to establish the desorption efficiencies from purely theoretical considerations. Applying this method to one such object (L16 89B) yields upper limits to the desorption efficiencies that are consistent with our understanding of these mechanisms.     

Nonlinear thin shell instabilities in molecular clouds

Observations of molecular clouds point to the existence of supersonic, turbulent flows. Therefore, any theory which attempts to describe molecular cloud evolution and star formation must include a consideration of the dynamics of colliding flows. Previous studies have considered the collision of supersonic streams or clouds. The resultant instabilities provide a mechanism which may give rise to observable cloud morphologies and enhance the star formation rate. One such instability is the nonlinear thin shell instability (NTSI) of a shock-bounded slab. This process is driven by ram pressure and efficient cooling. In this study, I use numerical simulations to examine the head-on collision of supersonic gas streams in a cold, molecular gas. A dense slab forms in the collision midplane and is prone to a number of instabilities, including the NTSI. The thermodynamic processes involved are found to have a controlling influence upon the instability and fragmentation of the slab. Although some minimal amount of cooling is needed to drive the instability, too rapid a cooling rate gives rise to smaller wavelength instabilities which wipe out the NTSI. The growth rate of the NTSI in a gas undergoing molecular cooling corresponds to a timescale of order 10¹² s, in general agreement with the theoretical value for an isothermal gas. The NTSI may provide a viable mechanism for the instigation of rapid star formation.     

Can molecular clouds live long?

It is generally accepted that the lifetime of molecular clouds does not exceed 3×10⁷ yr due to disruption by stellar feedback. We put together some arguments giving evidence that a substantial fraction of molecular clouds (primarily in the outer regions of a disc) may avoid destruction process for at least 10⁸ yr or even longer. A molecular cloud can live long if massive stars are rare or absent. Massive stars capable to destroy a cloud may not form for a long time if a cloud is low massive, or stellar initial mass function is top-light, or if there is a delay of the beginning of active star formation. A long duration of the inactive phase of clouds may be reconciled with the low amount of the observed starless giant molecular clouds if to propose that they were preceded by slowly contraction phase of the magnetized dark gas, non-detected in CO-lines.     

Three-dimensional simulations of molecular cloud fragmentation regulated by magnetic fields and ambipolar diffusion

We employ the first fully three-dimensional simulation to study the role of magnetic fields and ion–neutral friction in regulating gravitationally driven fragmentation of molecular clouds. The cores in an initially subcritical cloud develop gradually over an ambipolar diffusion time while the cores in an initially supercritical cloud develop in a dynamical time. The infalling speeds on to cores are subsonic in the case of an initially subcritical cloud, while an extended (≳0.1 pc) region of supersonic infall exists in the case of an initially supercritical cloud. These results are consistent with previous two-dimensional simulations. We also found that a snapshot of the relation between density (ρ) and the strength of the magnetic field ( B ) at different spatial points of the cloud coincides with the evolutionary track of an individual core. When the density becomes large, both the relations tend to   B ∝ρ^0.5  .     

巨分子云聚合形成N体模拟中粒子数的效应

通过N体数值模拟以不同粒子数对星系中巨分子云的形成的影响进行了研究。结果表明：在聚合形成机制下，当对不同数目的分子云取相同的平均密度时，基本分子云的有效半径与其相应质量的立方根成正比，巨分子云的碎裂率与模拟基本分子云的数目无关。     

Dynamics of molecular clouds on the galactic scale

We have carried out an extensive investigation into the dynamics of the molecular clouds in the disk of the Galaxy. We have used both computational methods and physical arguments to try to understand how the ensemble of molecular clouds interacts, how the clouds are affected by the gravitational field of the Galaxy and also the circumstances under which they can aggregate into giant molecular clouds (GMC's).The dynamical model is three dimensional and consists of 120,000 spherical clouds, each having a mass of 10⁴ M . It allows for the mutual gravitation between clouds, up to a cut-off distance; when two clouds collide they rebound, with a specified coefficient of restitutione. We have also developed a physically more realistic model for a cloud, supported by a magnetic field, and used it to select a suitable range of values fore. Our first paper deals with the case of an axially symmetrical galaxy. The clouds are distributed initially in a disk extending 100 pc on either side of the Galactic plane. As it evolves the system of clouds loses energy, and the disk grows thinner at a rate which depends on the value ofe. GMC's start to form once the disk is thin enough. We believe this result to be valid more generally, and that it holds also in models with spiral structure.     

Molecular cores of the high-latitude cloud MBM 40

Towards the high-latitude cloud MBM 40, we identify 3 dense molecular cores of M0.2–0.5 M, and sizes of 0.2 pc in diameter embedded in the H I cloud of 8 M which is observed to be extended along the northeast–southwest direction. The molecular cloud is located almost perpendicularly to the H I emission. We confirm the previous result of Magnani et al. that MBM 40 is not a site for new star formations. We found a very poor correlation between the H I and the IRAS 100 μm emissions, but the CO (1–0) and 100 μm emissions show a better correlation of W_CO/I₁₀₀=1±0.2 K km s⁻¹ (MJy sr⁻¹)⁻¹. This ratio is larger by a factor of ≥5 than in dense dark clouds, which may indicate that the CO is less depleted in MBM 40 than in dense dark clouds.     

The dependence of star formation on initial conditions and molecular cloud structure

We investigate the dependence of stellar properties on the initial kinematic structure of the gas in star-forming molecular clouds. We compare the results from two large-scale hydrodynamical simulations of star cluster formation that resolve the fragmentation process down to the opacity limit, the first of which was reported by Bate, Bonnell & Bromm. The initial conditions of the two calculations are identical, but in the new simulation the power spectrum of the velocity field imposed on the cloud initially and allowed to decay is biased in favour of large-scale motions. Whereas the calculation of Bate et al. began with a power spectrum   P ( k ) ∝ k ⁻⁴  to match the Larson scaling relations for the turbulent motions observed in molecular clouds, the new calculation begins with a power spectrum   P ( k ) ∝ k ⁻⁶  .
Despite this change to the initial motions in the cloud and the resulting density structure of the molecular cloud, the stellar properties resulting from the two calculations are indistinguishable. This demonstrates that the results of such hydrodynamical calculations of star cluster formation are relatively insensitive to the initial conditions. It is also consistent with the fact that the statistical properties of stars and brown dwarfs (e.g. the stellar initial mass function) are observed to be relatively invariant within our Galaxy and do not appear to depend on environment.     

Magnetic Fields in Star-Forming Regions: Theoretical Aspects

We review the basic theoretical elements leading to our current understanding of the role of magnetic fields in the process of star formation. In particular, we concentrate on: (i) the relevance of the mass-to-flux ratio for the stability of molecular clouds; (ii) the consequences of magnetic flux leakage for the evolution of cloud cores; (iii) the phase of anisotropic dynamical collapse following the formation of strongly peaked density distributions; (iv) the mechanism of magnetic braking as a possible solution to the angular momentum problem in star formation.     

Modeling Jupiter's cloud bands and decks: 2. Distribution and motion of condensates

A simple jovian cloud scheme has been developed for the Oxford Planetary Unified model System (OPUS). NH₃-ice, NH₄SH-solid, H₂O-ice and H₂O-liquid clouds have been modeled in Southern hemisphere limited area simulations of Jupiter. We found that either three or four of the condensates existed in the model. For a deep atmospheric water abundance close to solar composition, an NH₃-ice deck above 0.7 bar, an NH₄SH-solid deck above 2.5 bar and a H₂O-liquid deck with a base at about 7.5 bar and frozen cloud tops formed. If a depleted deep water abundance is assumed, however, a very compact cloud structure develops, where an H₂O-ice cloud forms by direct sublimation above 3 bar. The condensates constitute good tracers of atmospheric motion, and we have confirmed that zonal velocities determined from manual feature tracking in the modeled cloud layers agree reasonably well with the modeled zonal velocities. Dense and elevated clouds form over latitudes with strong atmospheric upwelling and depleted clouds exist over areas with strong downwelling. In the NH₃-ice deck this leads to elevated cloud bands over the zones in the domain and thin clouds over the belts, which is consistent with the observationally deduced distribution. Due to changes in the vertical velocity pattern in the deeper atmosphere, the NH₄SH-solid and water cloud decks are more uniform. This modeled cloud structure thus includes the possibility of more frequent water cloud observations in belts, as this deeper deck could be more easily detected under areas with thin NH₃-ice clouds. Large scale vortices appeared spontaneously in the model and were characterized by elevated NH₃-ice clouds, as expected from observations. These eddies leave the most discernible imprint on the lighter condensate particles of the uppermost layer.     

Measuring Venus’ winds using the Absolute Astronomical Accelerometer: Solid super-rotation model of Venus’ clouds

We present a new method of measuring the Venus winds by Doppler velocimetry on the full visible spectrum of solar light scattered by the clouds. In January 2003, we carried out observations to measure the winds of Venus, using the EMILIE high-resolution, cross-dispersed spectrograph and its associated calibrating instrument the Absolute Astronomical Accelerometer (AAA), at Observatoire de Haute-Provence, France. The motivation of this type of measurements is that it measures the actual velocity of cloud particles, while the other method (track of cloud features) may be sensitive to the deformation of the clouds. During observations, Venus was near maximum western elongation, at a phase angle near 90°. The EMILIE-AAA system allows us to measure accurately the Doppler shift induced in the reflected solar spectrum by the radial component of the motion of the clouds of Venus. We present the measurements and compare them with a forward simulation of a solid super-rotation of the atmosphere of Venus. Taking into account the Doppler shift relative to the Sun and that relative to the Earth, the theoretical total Doppler shift induced in the solar spectra is easily computed as a function of the velocity of the reflecting target. A first forward simulation is computed, with a wind model considering a purely horizontal and zonal wind. The magnitude of the wind is assumed to depend on cos(latitude), as for a solid-body rotation. The comparison with the measurements at various points on the illuminated semi-disc allowed us to determine an equatorial velocity of 66, 75, 91 and 85 m/s on 4 consecutive mornings, consistent with previous ultraviolet cloud tracking wind measurements, showing that wave propagation is not a major factor in the apparent motion of the cloud marks. Further, we discuss the effect of the finite angular size of the Sun and its rapid equatorial rotation (that we call the Young effect). It mainly affects measurements taken near the terminator, where the largest discrepancies are found. These discrepancies are alleviated when the Young effect is taken into account in the model but then the retrieved Venus equatorial velocity is reduced to only 48±3 m/s. This is well below classical ultraviolet markings velocities, but the altitude at which the visible photons are scattered (66 km) that we use is 5 km below the UV markings, confirming the vertical gradient of the horizontal winds shown by previous in-situ measurements.