Formation Of Molecular Clouds Research Articles

Chemical nitrogen fractionation in dense molecular clouds

Nitrogen-bearing molecules display variable isotopic fractionation levels in different astronomical environments such as in the interstellar medium or in the Solar system. Models of interstellar chemistry are unable to induce nitrogen fraction in cold molecular clouds as exchange reactions for 15N are mostly inefficient. Here, we developed a new gas–grain model for nitrogen fractionation including a thorough search for new nitrogen fractionation reactions and a realistic description of atom depletion on to interstellar dust particles. We show that, while dense molecular cloud gas-phase chemistry alone leads to very low fractionation, 14N atoms are preferentially depleted from the gas phase due to a mass-dependent grain surface sticking rate for atomic nitrogen. However, assuming an elementary 14N/15N ratio of 441 (equal to the solar wind value), our model leads to only low 15N enrichment for all N-containing species synthesized in the gas phase with predicted 14N/15N ratios in the range 360–400. Higher enrichment levels can neither be explained by this mechanism, nor through chemistry, with two possible explanations: (i) The elementary 14N/15N ratio in the local ISM is smaller, as suggested by the recent work of Romano et al., with a hypothetic 15NNH+ and 15NNH+ depletion due to variation of the electronic recombination rate constant variation with the isotopes and (ii) N2 photodissociation leads to variable nitrogen fractionation in diffuse molecular clouds where photons play an important role, which is conserved during dense molecular cloud formation as suggested by the work of Furuya & Aikawa.

On the shape and completeness of the column density probability distribution function of molecular clouds

Both observational and theoretical research over the past decade has demonstrated that the probability distribution function (PDF) of the gas density in turbulent molecular clouds is a key ingredient for understanding star formation. It has recently been argued that the PDF of molecular clouds is a pure power-law distribution. It has been claimed that the log-normal part is ruled out when using only the part of the PDF up/down to which it is complete, that is where the column density contours are still closed. By using the results from high-resolution magnetohydrodynamical simulations of molecular cloud formation and evolution, we find that the column density PDF is indeed composed of a log-normal and, if including self-gravity, a power-law part. We show that insufficient sampling of a molecular cloud results in closed contours that cut off the log-normal part. In contrast, systematically increasing the field of view and sampling the entire cloud yields a completeness limit at the lower column densities, which also recovers the log-normal part. This demonstrates that the field of view must be sufficiently large for the PDF to be complete down to its log-normal part, which has important implications for predictions of star-formation activity based on the PDF.

On estimating angular momenta of infalling protostellar cores from observations

We use numerical simulations of molecular cloud formation in the colliding flow scenario to investigate the reliability of observational estimates of the angular momenta of early-state, low-mass protostellar cores. We show that, with suitable corrections for projection factors, molecular line observations of velocity gradients in NH$_3$ can be used to provide reasonable estimates of core angular momenta within a factor of two to three, with a few large underestimates due to unfavorable viewing angles. Our results differ from previous investigations which suggested that observations might overestimate true angular momenta by as much as an order of magnitude; the difference is probably due to the much smoother velocity field on small scales in our simulations, which result from allowing turbulence to decay and gravitational infall to dominate. The results emphasize the importance of understanding the nature of 'turbulent' velocities, with implications for the formation of protostellar disks during core collapse.

The SILCC project – V. The impact of magnetic fields on the chemistry and the formation of molecular clouds

Magnetic fields are ubiquitously observed in the interstellar medium (ISM) of present-day star-forming galaxies with dynamically relevant energy densities. Using three-dimensional magneto-hydrodynamic (MHD) simulations of the supernova (SN) driven ISM in the flux-freezing approximation (ideal MHD) we investigate the impact of the magnetic field on the chemical and dynamical evolution of the gas, fragmentation and the formation of molecular clouds. We follow the chemistry with a network of six species (H$^{+}$, H, H$_2$, C$^+$, CO, free electrons) including local shielding effects. We find that magnetic fields thicken the disc by a factor of a few to a scale height of $\sim100\,\mathrm{pc}$, delay the formation of dense (and molecular) gas by $\sim25\,\mathrm{Myr}$ and result in differently shaped gas structures. The magnetised gas fragments into fewer clumps, which are initially at subcritical mass-to-flux ratios, $M/\Phi\approx0.3(M/\Phi)_\mathrm{crit}$, and accrete gas preferentially parallel to the magnetic field lines until supercritial mass-to-flux ratios of up to order 10 are reached. The accretion rates onto molecular clouds scale with $\dot{M}\propto M^{1.5}$. The median of the inter-cloud velocity dispersion is $\sim2-5\,\mathrm{km\,s}^{-1}$ and lower than the internal velocity dispersion in the clouds ($\sim3-7\,\mathrm{km\,s}^{-1}$). However, individual cloud-cloud collisions occur at speeds of a few $10\,\mathrm{km\,s}^{-1}$.

The role of magnetic field in the formation and evolution of filamentary molecular clouds

AbstractRecent observations have emphasized the importance of the formation and evolution of magnetized filamentary molecular clouds in the process of star formation. Theoretical and observational investigations have provided convincing evidence for the formation of molecular cloud cores by the gravitational fragmentation of filamentary molecular clouds. In this review we summarize our current understanding of various processes that are required in describing the filamentary molecular clouds. Especially we can explain a robust formation mechanism of filamentary molecular clouds in a shock compressed layer, which is in analogy to the making of “Sushi.” We also discuss the origin of the mass function of cores.

Observations of magnetic fields in star-forming clouds

AbstractThis review examines observations of magnetic fields in molecular clouds, that is, at spatial scales ranging from tens to tenths of parsecs and densities up to hundreds of particles per cubic centimetre. I will briefly summarize the techniques for observing and mapping magnetic fields in molecular clouds. I will review important examples of observational results obtained using each technique and their implications for our understanding of the role of the magnetic field in molecular cloud formation and evolution. Finally, I will briefly discuss the prospects for advances in our observational capabilities with telescopes and instruments now beginning operation or under construction.

Triggering the formation of the supergiant H ii region NGC 604 in M 33

Abstract Formation mechanism of a supergiant H ii region NGC 604 is discussed in terms of collision of H i clouds in M 33. An analysis of the archival H i data obtained with the Very Large Array (VLA) reveals complex velocity distributions around NGC 604. The H i clouds are composed of two velocity components separated by ∼20 km s−1 for an extent of ∼700 pc, beyond the size of the the H ii region. Although the H i clouds are not easily separated in velocity with some mixed component represented by merged line profiles, the atomic gas mass amounts to $6 \times 10^{6}\,M_{\odot }$ and $9 \times 10^{6}\,M_{\odot }$ for each component. These characteristics of H i gas and the distributions of dense molecular gas in the overlapping regions of the two velocity components suggest that the formation of giant molecular clouds and the following massive cluster formation have been induced by the collision of H i clouds with different velocities. Referring to the existence of a gas bridging feature connecting M 33 with M 31 reported by large-scale H i surveys, the disturbed atomic gas possibly represents the result of past tidal interaction between the two galaxies, which is analogous to the formation of the R 136 cluster in the LMC.

Molecular gas filamentary structures in galaxy clusters

AbstractRecent molecular line observations with ALMA and NOEMA in several Brightest Cluster Galaxies (BCG) have revealed the large-scale filamentary structure at the center of cool core clusters. These filaments extend over 20-100kpc, they are tightly correlated with ionized gas (Hα, [NII]) emission, and have characteristic shapes: either radial and straight, or also showing a U-turn, like a horse-shoe structure. The kinematics is quite regular and laminar, and the derived infall time is much longer than the free-fall time. The filaments extend up to the radius where the cooling time becomes larger than the infall time. Filaments can be perturbed by the sloshing of the BCG in its cluster, and spectacular cooling wakes have been observed. Filaments tend to occur at the border of cavities driven in the X-ray gas by the AGN radio jets. Observations of cool core clusters support the thermal instability scenario, which accounts for the multiphase medium in the upper atmospheres of BCG, where the right balance between heating and cooling is reached, and a chaotic cold gas accretion occurs. Molecular filaments are also seen associated to ram-pressure stripped spiral galaxies in rich galaxy clusters, and in jet-induced star formation, suggesting a very efficient molecular cloud formation even in hostile cluster environments.

Multiarm spirals on the periphery of disc galaxies

Abstract Spiral patterns in some disc galaxies have two arms in the centre, and three or more arms on the periphery. The same result is also obtained in numerical simulations of stellar and gaseous discs.We argue that such patterns may occur due to fast cooling of the gas, resulting in formation of giant molecular clouds. The timescale of this process is 50 Myr, the factor of 10 shorter than of ordinary secular instability. The giant molecular clouds give rise to multiarm spirals through the mechanism of swing amplification.

A new mechanical stellar wind feedback model for the Rosette Nebula

AbstractThe famous Rosette Nebula has an evacuated central cavity formed from the stellar winds ejected from the 2–6 Myr old codistant and comoving central star cluster NGC 2244. However, with upper age estimates of less than 110 000 yr, the central cavity is too young compared to NGC 2244 and existing models do not reproduce its properties. A new proper motion study herein using Gaia data reveals the ejection of the most massive star in the Rosette, HD 46223, from NGC 2244 occurred 1.73 (+0.34, −0.25) Myr (1σ uncertainty) in the past. Assuming this ejection was at the birth of the most massive stars in NGC 2244, including the dominant centrally positioned HD 46150, the age is set for the famous ionized region at more than 10 times that derived for the cavity. Here, we are able to reproduce the structure of the Rosette Nebula, through simulation of mechanical stellar feedback from a 40 M⊙ star in a thin sheet-like molecular cloud. We form the 135 000 M⊙ cloud from thermally unstable diffuse interstellar medium (ISM) under the influence of a realistic background magnetic field with thermal/magnetic pressure equilibrium. Properties derived from a snapshot of the simulation at 1.5 Myr, including cavity size, stellar age, magnetic field, and resulting inclination to the line of sight, match those derived from observations. An elegant explanation is thus provided for the stark contrast in age estimates based on realistic diffuse ISM properties, molecular cloud formation and stellar wind feedback.

Observations in the UV Band and Problems of Star Formation Studies

In this contribution, I consider those aspects of the modern star formation theory, which can be substantiated with observations in the UV band, paying special attention to early stages of molecular cloud formation and initial conditions for the chemical evolution of starless cores. I describe the main results of available diffuse cloud observations in the UV band, as well as prospective directions for future studies with the instrumentation of the World Space Observatory project.

Magnetic suppression of turbulence and the star formation activity of molecular clouds

We present magnetohydrodynamic simulations aimed at studying the effect of the magnetic suppression of turbulence (generated through various instabilities during the formation of molecular clouds by converging) on the subsequent star formation (SF) activity. We study four magnetically supercritical models with magnetic field strengths $B= 0$, 1, 2, and 3 $\mu$G (corresponding to mass--to--flux ratios of $\infty$, 4.76, 2.38, and 1.59 times the critical value), with the magnetic field, initially being aligned with the flows. We find that, for increasing magnetic field strength, the clouds formed tend to be more massive, denser, less turbulent, and with higher SF activity. This causes the onset of star formation activity in the non--magnetic or more weakly magnetized cases to be delayed by a few Myr in comparison to the more strongly magnetized cases. We attribute this behavior to the suppression of the nonlinear thin shell instability (NTSI) by the magnetic field, previously found by Heitsch and coworkers. This result is contrary to the standard notion that the magnetic field provides support to the clouds, thus reducing their star formation rate (SFR). However, our result is a completely nonlinear one, and could not be foreseen from simple linear considerations.

The numerical modelling of MHD astrophysical flows with chemistry

The new code for numerical simulation of magnetic hydrodynamical astrophysical flows with consideration of chemical reactions is given in the paper. At the heart of the code – the new original low-dissipation numerical method based on a combination of operator splitting approach and piecewise-parabolic method on the local stencil. The chemodynamics of the hydrogen while the turbulent formation of molecular clouds is modeled.

The driving of turbulence in simulations of molecular cloud formation and evolution

Molecular clouds are to a great extent influenced by turbulent motions in the gas. Numerical and observational studies indicate that the star formation rate and efficiency crucially depend on the mixture of solenoidal and compressive modes in the turbulent acceleration field, which can be quantified by the turbulent driving parameter b. For purely solenoidal (divergence-free) driving previous studies showed that b=1/3 and for entirely compressive (curl-free) driving b=1. In this study, we determine the evolution of the turbulent driving parameter b in magnetohydrodynamical simulations of molecular cloud formation and evolution. The clouds form due to the convergence of two flows of warm neutral gas. We explore different scenarios by varying the magnitude of the initial turbulent perturbations in the flows. We show that the driving mode of the turbulence within the cloud strongly fluctuates with time and exhibits no clear correlation with typical cloud properties, such as the cloud mass and the (Alfven) Mach number. We specifically find that $b$ strongly varies from b=0.3 to b=0.8 on timescales t<5 Myr, where the timescale and range of variation can change from cloud to cloud. This rapid change of b from solenoidal to compressive driving is primarily associated with global contraction of the cloud and subsequent onset of star formation. We conclude that the effective turbulence driving parameter should be treated as a free parameter that can vary from solenoidal to compressive in both time and space.

Molecular cloud formation via thermal instability of finite resistive viscous radiating plasma with finite Larmor radius corrections

The effect of radiative heat-loss function and finite ion Larmor radius (FLR) corrections on the thermal instability of infinite homogeneous viscous plasma has been investigated incorporating the effects of thermal conductivity and finite electrical resistivity for the formation of a molecular cloud. The general dispersion relation is derived using the normal mode analysis method with the help of relevant linearized perturbation equations of the problem. Furthermore the wave propagation along and perpendicular to the direction of external magnetic field has been discussed. Stability of the medium is discussed by applying Routh Hurwitz’s criterion and it is found that thermal instability criterion determines the stability of the medium. We find that the presence of radiative heat-loss function and thermal conductivity modify the fundamental criterion of thermal instability into radiatively driven thermal instability criterion. In longitudinal direction FLR corrections, viscosity, magnetic field and finite resistivity have no effect on thermal instability criterion. The presence of radiative heat-loss function and thermal conductivity modify the fundamental thermal instability criterion into radiatively driven thermal instability criterion. Also the FLR corrections modify the growth rate of the Alfven mode. For transverse wave propagation FLR corrections, radiative heat-loss function, magnetic field and thermal conductivity modify the thermal instability criterion. From the curves it is clear that heat-loss function, FLR corrections and viscosity have stabilizing effect, while finite resistivity has destabilizing effect on the thermal modes. Our results show that the FLR corrections and radiative heat-loss functions affect the evolution of interstellar molecular clouds and star formation.

Jeans instability of inhomogeneous dusty plasma with polarization force, ionization and recombination

The self-gravitational Jeans instability has been studied in dusty plasma containing significant background of neutral pressure and recombination of ions and electrons on the dust surface. The full dynamics of charged dust grains, ions and neutral species are employed considering the electrons as Maxwellian. We have derived the general dispersion relation for collisional dusty plasma with ionization, recombination and polarization force. The general dispersion relation describes the effects of considered parameters which are solved in different dusty plasma situations. Further, the dispersion relation is solved numerically. The present work is applicable to understand the structure formation of interstellar molecular clouds in astrophysical plasma.

The role of magnetic fields in the structure and interaction of supershells

Context. Large-scale shocks formed by clustered feedback of young OB stars are considered an important source of mechanical energy for the interstellar medium (ISM) and a trigger of molecular cloud formation. Their interaction sites are locations where kinetic energy and magnetic fields are redistributed between ISM phases. Aims. In this work we address two questions, both involving the role of galactic magnetic fields in the dynamics of supershells and their interactions. On the one hand, we study the effect of the magnetic field on the expansion and fragmentation of supershells and, on the other hand, we look for the signatures of supershell collisions on dense structures and on the kinetic and magnetic energy distribution of the ISM. Methods. We performed a series of high-resolution, three-dimensional simulations of colliding supershells. The shocks are created by time-dependent feedback and evolve in a diffuse turbulent environment that is either unmagnetized or has different initial magnetic field configurations. Results. In the hydrodynamical situation, the expansion law of the superbubbles is consistent with the radius-time relation R ∝ t3/5 that is theoretically predicted for wind-blown bubbles. The supershells fragment over their entire surface into small dense clumps that carry more than half of the total kinetic energy in the volume. However, this is not the case when a magnetic field is introduced, either in the direction of the collision or perpendicular to the collision. In both situations, the shell surfaces are more stable to dynamical instabilities. When the magnetic field opposes the collision, the expansion law of the supershells also becomes significantly flatter than in the hydrodynamical case. Although a two-phase medium arises in all cases, in the magnetohydrodynamical (MHD) simulations the cold phase is limited to lower densities and the cold clumps are located further away from the shocks with respect to the hydrodynamical simulations. Conclusions. For the parameters we explored, self-gravity has no effect on either the superbubble expansion or the shock fragmentation. In contrast, a magnetic field, whether mostly parallel or mostly perpendicular to the collision axis, causes a deceleration of the shocks, deforms them significantly, and largely suppresses the formation of the dense gas on their surface. The result is a multi-phase medium in which the cold clumps are not spatially correlated with the supershells.

Magnetohydrodynamic Simulations of the Formation of Molecular Clouds toward the Stellar Cluster Westerlund 2: Interaction of a Jet with a Clumpy Interstellar Medium

Abstract The formation mechanism of CO clouds observed with the NANTEN2 and Mopra telescopes toward the stellar cluster Westerlund 2 is studied by 3D magnetohydrodynamic simulations, taking into account the interstellar cooling. These molecular clouds show a peculiar shape composed of an arc-shaped cloud on one side of the TeV γ-ray source HESS J1023-575 and a linear distribution of clouds (jet clouds) on the other side. We propose that these clouds are formed by the interaction of a jet with clumps of interstellar neutral hydrogen (H i). By studying the dependence of the shape of dense cold clouds formed by shock compression and cooling on the filling factor of H i clumps, we found that the density distribution of H i clumps determines the shape of molecular clouds formed by the jet–cloud interaction: arc clouds are formed when the filling factor is large. On the other hand, when the filling factor is small, molecular clouds align with the jet. The jet propagates faster in models with small filling factors.

Galaxy chemical evolution models: the role of molecular gas formation

In our grid of multiphase chemical evolution models (Moll\'a & D\'iaz, 2005), star formation in the disk occurs in two steps: first, molecular gas forms, and then stars are created by cloud-cloud collisions or interactions of massive stars with the surrounding molecular clouds. The formation of both molecular clouds and stars are treated through the use of free parameters we refer to as efficiencies. In this work we modify the formation of molecular clouds based on several new prescriptions existing in the literature, and we compare the results obtained for a chemical evolution model of the Milky Way Galaxy regarding the evolution of the Solar region, the radial structure of the Galactic disk, and the ratio between the diffuse and molecular components, HI/H$_2$. Our results show that the six prescriptions we have tested reproduce fairly consistent most of the observed trends, differing mostly in their predictions for the (poorly-constrained) outskirts of the Milky Way and the evolution in time of its radial structure. Among them, the model proposed by Ascasibar et al. (2017), where the conversion of diffuse gas into molecular clouds depends on the local stellar and gas densities as well as on the gas metallicity, seems to provide the best overall match to the observed data.

Velocity Anisotropy in Self-gravitating Molecular Clouds. I. Simulation

Abstract The complex interplay between turbulence, magnetic fields, and self-gravity leads to the formation of molecular clouds out of the diffuse interstellar medium (ISM). One avenue of studying this interplay is by analyzing statistical features derived from observations, where the interpretation of these features is greatly facilitated by comparisons with numerical simulations. Here we focus on the statistical anisotropy present in synthetic maps of velocity centroid data, which we derive from three-dimensional magnetohydrodynamic simulations of a turbulent, magnetized, self-gravitating patch of ISM. We study how the orientation and magnitude of the velocity anisotropy correlate with the magnetic field and with the structures generated by gravitational collapse. Motivated by recent observational constraints, our simulations focus on the supersonic (sonic Mach number  ≈ 2 – 17 ) but sub- to trans-alfvénic (alfvénic Mach number  A ≈ 0.2 – 1.2 ) turbulence regime, and we consider clouds that are barely to mildly magnetically supercritical (mass-to-flux ratio equal to once or twice the critical value). Additionally we explore the impact of the turbulence driving mechanism (solenoidal or compressive) on the velocity anisotropy. While we confirm previous findings that the velocity anisotropy generally aligns well with the plane-of-sky magnetic field, our inclusion of the effects of self-gravity reveals that in regions of higher column density, the velocity anisotropy may be destroyed or even reoriented to align with the gravitationally formed structures. We provide evidence that this effect is not necessarily due to the increase of  A inside the high-density regions.