On the structure of the turbulent interstellar atomic hydrogen. I- Physical characteristics

{We study in some details the statistical properties of the turbulent 2-phase interstellar atomic gas.{We present high resolution bidimensional numerical simulations of the interstellar atomic hydrogen which describe it over 3 to 4 orders of magnitude in spatial scales.}{The simulations produce naturally small scale structures having either large or small column density. It is tempting to propose that the former are connected to the tiny small scale structures observed in the ISM. We compute the mass spectrum of CNM structures and find that ${\cal N}(M) dM \propto M ^{-1.7} dM$, which is remarkably similar to the mass spectrum inferred for the CO clumps. We propose a theoretical explanation based on a formalism inspired from the Press & Schecter (1974) approach and used the fact that the turbulence within WNM is subsonic. This theory predicts ${\cal N}(M) \propto M ^{-5/3}$ in 2D and ${\cal N}(M) \propto M ^{-16/9}$ in 3D. We compute the velocity and the density power-spectra and conclude that, although the latter is rather flat, as observed in supersonic isothermal simulations, the former follows the Kolmogorov prediction and is dominated by its solenoidal component. This is due to the bistable nature of the flow which produces large density fluctuations even when the rms Mach number (of WNM) is not large. We also find that, whereas the energy at large scales is mainly in the WNM, at smaller scales, it is dominated by the kinetic energy of the CNM fragments.}Comment: Accepted for publication in A&

Synthetic Observations of Carbon Lines of Turbulent Flows in Diffuse Multiphase Interstellar Medium

We examine observational characteristics of multi-phase turbulent flows in the diffuse interstellar medium (ISM) using a synthetic radiation field of atomic and molecular lines. We consider the multi-phase ISM which is formed by thermal instability under the irradiation of UV photons with moderate visual extinction $A_V\sim 1$. Radiation field maps of C$^{+}$, C$^0$, and CO line emissions were generated by calculating the non-local thermodynamic equilibrium (nonLTE) level populations from the results of high resolution hydrodynamic simulations of diffuse ISM models. By analyzing synthetic radiation field of carbon lines of [\ion{C}{2}] 158 $\mu$m, [\ion{C}{1}] $^3P_2-^3P_1$ (809 GHz), $^3P_1-^3P_0$ (492 GHz), and CO rotational transitions, we found a high ratio between the lines of high- and low-excitation energies in the diffuse multi-phase interstellar medium. This shows that simultaneous observations of the lines of warm- and cold-gas tracers will be useful in examining the thermal structure, and hence the origin of diffuse interstellar clouds.Comment: 16 pages, 10 figures : accepted for publication in ApJ. PDF version with high resolution figures is available (http://yso.mtk.nao.ac.jp/~ymasako/paper/ms_hires.pdf

Semi-analytical homologous solutions of the gravo-magnetic contraction

We propose an extension of the semi-analytical solutions derived by Lin et al. (1965) describing the two-dimensional homologous collapse of a self-gravitating rotating cloud having uniform density and spheroidal shape, which includes magnetic field (with important restrictions) and thermal pressure. The evolution of the cloud is reduced to three time dependent ordinary equations allowing to conduct a quick and preliminary investigation of the cloud dynamics during the precollapse phase, for a wide range of parameters. We apply our model to the collapse of a rotating and magnetized oblate and prolate isothermal core. Hydrodynamical numerical simulations are performed and comparison with the semi-analytical solutions is discussed. Under the assumption that all cores are similar, an apparent cloud axis ratio distribution is calculated from the sequence of successive evolutionary states for each of a large set of initial conditions. The comparison with the observational distribution of the starless dense cores belonging to the catalog of Jijina et al. (1999) shows a good agreement for the rotating and initially prolate cores (aspect ratio $\simeq 0.5$) permeated by an helical magnetic field ($\simeq 17-20 \mu$G for a density of $\simeq 10^4$ cm$^{-3}$).Comment: accepted for publication in A&

Two-Fluid MHD Simulations of Converging HI Flows in the Interstellar Medium. I: Methodology and Basic Results

We develop an unconditionally stable numerical method for solving the coupling between two fluids (frictional forces/heatings, ionization, and recombination), and investigate the dynamical condensation process of thermally unstable gas that is provided by the shock waves in a weakly ionized and magnetized interstellar medium by using two-dimensional two-fluid magnetohydrodynamical simulations. If we neglect the effect of magnetic field, it is known that condensation driven by thermal instability can generate high density clouds whose physical condition corresponds to molecular clouds (precursor of molecular clouds). In this paper, we study the effect of magnetic field on the evolution of supersonic converging HI flows and focus on the case in which the orientation of magnetic field to converging flows is orthogonal. We show that the magnetic pressure gradient parallel to the flows prevents the formation of high density and high column density clouds, but instead generates fragmented, filamentary HI clouds. With this restricted geometry, magnetic field drastically diminishes the opportunity of fast molecular cloud formation directly from the warm neutral medium, in contrast to the case without magnetic field.Comment: ApJ accepte

A dynamical model for the dusty ring in the Coalsack

Lada et al. recently presented a detailed near-infrared extinction map of Globule G2 in the Coalsack molecular cloud complex, showing that this starless core has a well-defined central extinction minimum. We propose a model for G2 in which a rapid increase in external pressure is driving an approximately symmetric compression wave into the core. The rapid increase in external pressure could arise because the core has recently been assimilated by the Coalsack cloud complex, or because the Coalsack has recently been created by two large-scale converging flows. The resulting compression wave has not yet converged on the centre of the core, so there is a central rarefaction. The compression wave has increased the density in the swept-up gas by about a factor of ten, and accelerated it inwards to speeds of order $0.4 {\rm km} {\rm s}^{-1}$. It is shown that even small levels of initial turbulence destroy the ring seen in projection almost completely. In the scenario of strong external compression that we are proposing this implies that the initial turbulent energy in this globule is such that $E_{{\rm turb}} / E_{{\rm grav}} \le 2$. Protostar formation should occur in about $40,000 {\rm years}$.Comment: Accepted for publication in A&

H2 distribution during 2-phase Molecular Cloud Formation

We performed high-resolution, 3D MHD simulations and we compared to observations of translucent molecular clouds. We show that the observed populations of rotational levels of H2 can arise as a consequence of the multi-phase structure of the ISM.Comment: 2 pages, 1 figure. Due to appear in the proceedings of the 6th Zermatt ISM Symposium: "Conditions and Impact of Star Formation: From Lab to Space

Disk formation during collapse of magnetized protostellar cores

In the context of star and planet formation, understanding the formation of disks is of fundamental importance. Previous studies found that the magnetic field has a very strong impact on the collapse of a prestellar cloud, particularly in possibly suppressing the formation of a disk even for relatively modest values of the magnetic intensity. Since observations infer that cores have a substantial level of magnetization, this raises the question of how disks form. However, most studies have been restricted to the case in which the initial angle, $\alpha$, between the magnetic field and the rotation axis equals 0$^\circ$. We explore and analyse the influence of non aligned configurations on disk formation. We perform 3D ideal MHD, AMR numerical simulations for various values of $\mu$, the ratio of the mass-to-flux to the critical mass-to-flux, and various values of $\alpha$. We find that disks form more easily as $\alpha$ increases from 0 to 90$^\circ$. We propose that as the magnetized pseudo-disks become thicker with increasing $\alpha$, the magnetic braking efficiency is lowered. We also find that even small values of $\alpha$ ($\simeq$ 10-20$^\circ$) show significant differences with the alligned case. Within the framework of ideal MHD and for our choice of initial conditions, centrifugally supported disks cannot form for values of $\mu$ smaller than $\simeq$3, if the magnetic field and the rotation axis are perpendicular, and smaller than about $\simeq$5-10 when they are perfectly aligned.Comment: accepted for publication in A&

From the warm magnetized atomic medium to molecular clouds

{It has recently been proposed that giant molecular complexes form at the sites where streams of diffuse warm atomic gas collide at transonic velocities.} {We study the global statistics of molecular clouds formed by large scale colliding flows of warm neutral atomic interstellar gas under ideal MHD conditions. The flows deliver material as well as kinetic energy and trigger thermal instability leading eventually to gravitational collapse.} {We perform adaptive mesh refinement MHD simulations which, for the first time in this context, treat self-consistently cooling and self-gravity.} {The clouds formed in the simulations develop a highly inhomogeneous density and temperature structure, with cold dense filaments and clumps condensing from converging flows of warm atomic gas. In the clouds, the column density probability density distribution (PDF) peaks at \sim 2 \times 10^{21} \psc and decays rapidly at higher values; the magnetic intensity correlates weakly with density from $n \sim 0.1$ to 10^4 \pcc, and then varies roughly as $n^{1/2}$ for higher densities.} {The global statistical properties of such molecular clouds are reasonably consistent with observational determinations. Our numerical simulations suggest that molecular clouds formed by the moderately supersonic collision of warm atomic gas streams.}Comment: submitted to A&