We present numerical experiments aimed to study the correlation between the magnetic field strength, $B$, and the density, $n$, in the cold atomic interstellar medium (CNM). We analyze 24 magneto-hydrodynamic models with different initial magnetic field intensities ($B_0=$0.4, 2.1, 4.2, and 8.3 $\mu$G) and/or mean densities (2, 3, and 4 cm$^{-3}$), in the presence of driven and decaying turbulence, with and without self-gravity, in a cubic computational domain with 100 pc by side. Our main findings are: i) For forced simulations, which reproduce the main observed physical conditions of the CNM in the Solar neighborhood, a positive correlation between $B$ and $n$ develops for all the $B_0$ values. ii) The density at which this correlation becomes significant ($\lesssim 30$ cm$^{-3}$) depends on $B_0$ but is not sensitive to the presence of self-gravity. iii) The effect of self-gravity, when noticeable, consists of producing a shallower correlation at high densities, suggesting that, in the studied regime, self-gravity induces motions along the field lines. iv) Self-gravitating decaying models where the CNM is subsonic and sub-Alfv\'enic with $\beta\lesssim 1$ develop a high density positive correlation whose slopes are consistent with a constant $\beta(n)$. v) Decaying models where the low density CNM is subsonic and sub-Alfv\'enic with $\beta>1$ show a negative correlation at intermediate densities, followed by a high density positive correlation.

To contribute to the understanding of the magnetic field's influence on the segregation of CNM in the solar neighbourhood we analyse MHD simulations which include the main physical characteristics of the local neutral atomic ISM. The simulations have a continuous solenoidal Fourier forcing in a periodic box of 100 pc per side and an initial uniform magnetic field ($\vec{B_0}$) with intensities ranging between $\sim 0.4$ $\mu$G and $\sim 8$ $\mu$G. Our main results are: i) the CNM mass fraction diminishes with the increase in magnetic field intensity. ii) There is a preferred alignment between CNM structures and $\vec{B}$ in all our $B_0$ range but the preference weakens as $B_{0}$ increases. It is worth noticing that this preference is also present in two-dimensional projections making an extreme angle ($0$ or $\pi / 2$) with respect to $\vec{B_0}$ and it is only lost for the strongest magnetic field when the angle of projection is perpendicular to $\vec{B_0}$. iii) The aforementioned results are prevalent despite the inclusion of self-gravity in our continuously forced simulations with a mean density similar to the average value of the solar neighbourhood. iv) Given a fixed $B_0$ and slightly higher mean densities, up to double, the effects of self-gravity are still not qualitatively significant.

We investigate the scale dependence of fluctuations inside a realistic model of an evolving turbulent HII region and to what extent these may be studied observationally. We find that the multiple scales of energy injection from champagne flows and the photoionization of clumps and filaments leads to a flatter spectrum of fluctuations than would be expected from top-down turbulence driven at the largest scales. The traditional structure function approach to the observational study of velocity fluctuations is shown to be incapable of reliably determining the velocity power spectrum of our simulation. We find that a more promising approach is the Velocity Channel Analysis technique of Lazarian & Pogosyan (2000), which, despite being intrinsically limited by thermal broadening, can successfully recover the logarithmic slope of the velocity power spectrum to a precision of +-0.1 from high resolution optical emission line spectroscopy.

We numerically study the volume density probability distribution function (n-PDF) and the column density probability distribution function (Sigma-PDF) resulting from thermally bistable turbulent flows. We analyze three-dimensional hydrodynamic models in periodic boxes of 100pc by side, where turbulence is driven in the Fourier space at a wavenumber corresponding to 50pc. At low densities (n <= 0.6cm^-3) the n-PDF, is well described by a lognormal distribution for average local Mach number ranging from ~0.2 to ~5.5. As a consequence of the non linear development of thermal instability (TI), the logarithmic variance of the distribution for the diffuse gas increases with M faster than in the well known isothermal case. The average local Mach number for the dense gas (n >= 7.1cm^-3) goes from ~1.1 to ~16.9 and the shape of the high density zone of the n-PDF changes from a power-law at low Mach numbers to a lognormal at high M values. In the latter case the width of the distribution is smaller than in the isothermal case and grows slower with M. At high column densities the Sigma-PDF is well described by a lognormal for all the Mach numbers we consider and, due to the presence of TI, the width of the distribution is systematically larger than in the isothermal case but follows a qualitatively similar behavior as M increases. Although a relationship between the width of the distribution and M can be found for each one of the cases mentioned above, these relations are different form those of the isothermal case.

It has been recently shown that molecular clouds do not exhibit a unique shape for the column density probability distribution function (Npdf). Instead, clouds without star formation seem to possess a lognormal distribution, while clouds with active star formation develope a power-law tail at high column densities. The lognormal behavior of the Npdf has been interpreted in terms of turbulent motions dominating the dynamics of the clouds, while the power-law behavior occurs when the cloud is dominated by gravity. In the present contribution we use thermally bi-stable numerical simulations of cloud formation and evolution to show that, indeed, these two regimes can be understood in terms of the formation and evolution of molecular clouds: a very narrow lognormal regime appears when the cloud is being assembled. However, as the global gravitational contraction occurs, the initial density fluctuations are enhanced, resulting, first, in a wider lognormal Npdf, and later, in a power-law Npdf. We thus suggest that the observed Npdf of molecular clouds are a manifestation of their global gravitationally contracting state. We also show that, contrary to recent suggestions, the exact value of the power-law slope is not unique, as it depends on the projection in which the cloud is being observed.

In this paper we numerically study the behavior of the density power spectrum in turbulent thermally bistable flows. We analyze a set of five three-dimensional simulations where turbulence is randomly driven in Fourier space at a fixed wave-number and with different Mach numbers M (with respect to the warm medium) ranging from 0.2 to 4.5. The density power spectrum becomes shallower as M increases and the same is true for the column density power spectrum. This trend is interpreted as a consequence of the simultaneous turbulent compressions, thermal instability generated density fluctuations, and the weakening of thermal pressure force in diffuse gas. This behavior is consistent with the fact that observationally determined spectra exhibit different slopes in different regions. The values of the spectral indexes resulting from our simulations are consistent with observational values. We do also explore the behavior of the velocity power spectrum, which becomes steeper as M increases. The spectral index goes from a value much shallower than the Kolmogorov one for M=0.2 to a value steeper than the Kolmogorov one for M=4.5.

We present three numerical simulations of randomly driven, isothermal, non-magnetic, self-gravitating turbulence with different rms Mach numbers Ms and physical sizes L, but approximately the same value of the virial parameter, alpha approx 1.2. We obtain the following results: a) We test the hypothesis that the collapsing centers originate from locally Jeans-unstable ("super-Jeans"), subsonic fragments; we find no such structures. b) We find that the fraction of small-scale super-Jeans structures is larger in the presence of self-gravity. c) The velocity divergence of subregions of the simulations exhibits a negative correlation with their mean density. d) The density probability density function (PDF) deviates from a lognormal in the presence of self-gravity. e) Turbulence alone in the large-scale simulation does not produce regions with the same size and mean density as those of the small-scale simulation. Items (b)-(e) suggest that self-gravity is not only involved in causing the collapse of Jeans-unstable density fluctuations produced by the turbulence, but also in their {it formation}. We also measure the star formation rate per free-fall time, as a function of Ms for the three runs, and compare with the predictions of recent semi-analytical models. We find marginal agreement to within the uncertainties of the measurements. However, the hypotheses of those models neglect the net negative divergence of dense regions we find in our simulations. We conclude that a) part of the observed velocity dispersion in clumps must arise from clump-scale inwards motions, and b) analytical models of clump and star formation need to take into account this dynamical connection with the external flow and the fact that, in the presence of self-gravity, the density PDF may deviate from a lognormal.

We discuss molecular cloud formation by large-scale supersonic compressions in the diffuse warm neutral medium (WNM). Initially, a shocked layer forms, and within it, a thin cold layer. An analytical model and high-resolution 1D simulations predict the thermodynamic conditions in the cold layer. After $\sim 1$ Myr of evolution, the layer has column density $\sim 2.5 \times 10^{19} \psc$, thickness $\sim 0.03$ pc, temperature $\sim 25$ K and pressure $\sim 6650$ K $\pcc$. These conditions are strongly reminiscent of those recently reported by Heiles and coworkers for cold neutral medium sheets. In the 1D simulations, the inflows into the sheets produce line profiles with a central line of width $\sim 0.5 \kms$ and broad wings of width $\sim 1 \kms$. 3D numerical simulations show that the cold layer develops turbulent motions and increases its thickness, until it becomes a fully three-dimensional turbulent cloud. Fully developed turbulence arises on times ranging from $\sim 7.5$ Myr for inflow Mach number $\Mr = 2.4$ to $> 80$ Myr for $\Mr = 1.03$. These numbers should be considered upper limits. The highest-density turbulent gas (HDG, $n > 100 \pcc$) is always overpressured with respect to the mean WNM pressure by factors 1.5--4, even though we do not include self-gravity. The intermediate-density gas (IDG, $10 < n [{\rm cm}^ {-3}] < 100$) has a significant pressure scatter that increases with $\Mr$, so that at $\Mr = 2.4$, a significant fraction of the IDG is at a higher pressure than the HDG. Our results suggest that the turbulence and at least part of the excess pressure in molecular clouds can be generated by the compressive process that forms the clouds themselves, and that thin CNM sheets may be formed transiently by this mechanism, when the compressions are only weakly supersonic.

We investigate the core mass distribution (CMD) resulting from numerical models of turbulent fragmentation of molecular clouds. In particular we study its dependence on the sonic root-mean-square Mach number $\Ms$. We analyze simulations with $\Ms$ ranging from 1 to 15 to show that, as $\Ms$ increases, the number of cores increases as well while their average mass decreases. This stems from the fact that high-Mach number flows produce many and strong shocks on intermediate to small spatial scales, leading to a highly-fragmented density structure. We also show that the CMD from purely turbulent fragmentation does not follow a single power-law, but it may be described by a function that changes continuously its shape, probably more similar to a log-normal function. The CMD in supersonic turbulent flows does not have a universal slope, and as consequence, cast some doubt on attempts to directly relate the CMD to a universal Initial Mass Function.

We present a systematic numerical study of the effect of turbulent velocity fluctuations on the thermal pressure distribution in thermally bistable flows. The simulations employ a random turbulent driving generated in Fourier space rather than star-like heating. The turbulent fluctuations are characterized by their rms Mach number M and the energy injection wavenumber, k_for. Our results are consistent with the picture that as either of these parameters is increased, the local ratio of turbulent crossing time to cooling time decreases, causing transient structures in which the effective behavior is intermediate between the thermal-equilibrium and adiabatic regimes. As a result, the effective polytropic exponent gamma_ef ranges between ~0.2 to ~1.1. The fraction of high-density zones with P>10^4 Kcm^-3 increases from roughly 0.1% at k_for=2 and M=0.5 to roughly 70% for k_for=16 and M=1.25. A preliminary comparison with the pressure measurements of Jenkins (2004) favors our case with M=0.5 and k_for=2. In all cases, the dynamic range of the pressure summed over the entire density range, typically spans 3-4 orders of magnitude. The total pressure histogram widens as the Mach number is increased, and develops near-power-law tails at high (resp.low) pressures when gamma_ef<~ 0.5 (resp. gamma_ef>~ 1), which occurs at k_for=2 (resp.k_for=16) in our simulations. The opposite side of the pressure histogram decays rapidly, in an approx. lognormal form. Our results show that turbulent advection alone can generate large pressure scatters, with power-law high-P tails for large-scale driving, and provide validation for approaches attempting to derive the shape of the pressure histogram through a change of variable from the known form of the density histogram, such as that performed by MacLow et al.(2004).

This paper is devoted to the study of the gravitational instability of a medium permeated by a uniform magnetic field along which a circularly polarized Alfven wave propagates. We concentrate on the case of perturbations purely transverse to the ambient field by means of direct numerical simulations of the MHD equations and of a linear stability analysis performed on a moderate amplitude asymptotic model. The Alfven wave provides an extra stabilizing pressure when the scale of perturbations is sufficiently large or small compared with the Jeans length $L_J$. However, there is a band of scales around $L_J$ for which the Alfven wave is found to have a destabilizing effect. In particular, when the medium is stable in absence of waves, the gravitational instability can develop when the wave amplitude lies in an appropriate range. This effect appears to be a consequence of the coupling between Alfven and magnetosonic waves. The prediction based on a WKB approach that the Alfven wave pressure tensor is isotropic and thus opposes gravity in all directions is only recovered for large amplitude waves for which the coupling between the different MHD modes is negligible.

We review recent results on the nonlinear development of thermal instability in the context of the turbulent atomic interstellar medium (ISM). First, we pre- sent a brief summary of the linear theory, remarking that, in the atomic ISM, the wave mode is stable at small scales. Next, we revisit the growth of isolated entropy perturbations in initially unstable gas, as a function of the ratio $\eta$ of the cooling to the dynamical crossing times. Third, we consider the evolution of {\it velocity} perturbations. These correspond to the wave mode, and are stable at moderate amplitudes and small scales, as confirmed numerically. Fourth, we consider the behavior of magnetic pressure in turbulent regimes. We propose that recent findings of a poor B-rho correlation at low rho are due to the different B-rho scalings for the slow and fast modes of non- linear MHD waves. This implies that, in fully turbulent regimes, the magnetic field may not be a very efficient source of pressure, and that polytropic de- scriptions of magnetic pressure are probably not adequate. Finally, we discuss simulations of the ISM (and resolution issues) concerned with the possibility of significant amounts of gas being in the ``lukewarm'' temperature range be- tween the warm and cold stable phases. The mass fraction in this range in- creases, and the phase segregation decreases, as smaller scales are considered. We attribute this to the enhanced stability of moderate, adiabatic-like veloc- ity fluctuations with $\eta \gg 1$, to the recycling of gas from the dense to the diffuse phase by stellar energy injection, and to the magnetic field not being strongly turbulent there, possibly providing additional stability. Final- ly, we suggest that the lukewarm gas can be observationally distinguished through simultaneous determination of two of its thermodynamic variables.

(Abridged) We discuss the nonlinear development of the isobaric mode of thermal instability (TI) in the context of the atomic interstellar medium (ISM), in order to assess the ability of TI to establish a well-segregated multi-phase structure in the turbulent ISM. The key parameter is the ratio of the cooling time to the dynamical crossing time $\eta$. The isobaric density perturbations of sizes $\gtrsim 15$ pc in media with mean density $\sim 1$ cm$^{-3}$ develop inflow motions with Mach number larger than 0.5 and a shock that propagates outwards from the condensation and brings the surrounding medium out of thermal equilibrium. The time for the dynamical transient state to subside ranges from 4 to 30 Myr for initial density perturbations of 20% and sizes 3 to 75 pc. By the time the condensations have formed, a substantial fraction of the mass is still traversing the unstable range and, therefore, clouds formed by TI should be bounded by accreting gas traversing the unstable range, rather than by sharp transitions to the stable warm phase. In the presence of random velocity forcing the condensation process can be suppressed for arbitrarily long times if the forcing causes a moderate rms Mach number ($\gtrsim 0.3$) and extends to small enough scales that $\eta>1$. We suggest that these mechanisms may be at the origin of the relatively large amounts of gas mass in the unstable regime. We remark that in the (stable) warm diffuse medium, $\eta>1$ for velocity perturbations of scales up to several pc. The flow's response to them is thus nearly adiabatic, and relatively weakly compressible, consistent with observations that suggest a nearly Kolmogorov power spectrum in this medium.

We discuss the temperature distribution in a two-dimensional, thermally unstable numerical simulation of the warm and cold gas in the Galactic disk, including the magnetic field, self-gravity, the Coriolis force, stellar energy injection and a realistic cooling function. We find that ~50% of the turbulent gas mass has temperatures in what would be the thermally unstable range if thermal instability were to be considered alone. This appears to be a consequence of there being many other forces at play than just thermal pressure. We also point out that a bimodal temperature pdf is a simple consequence of the form of the interstellar cooling function and is not necessarily a signature of discontinuous phase transitions.

We investigate numerically the role of thermal instability (TI) as a generator of density structures in the interstellar medium (ISM), both by itself and in the context of a globally turbulent medium. Simulations of the instability alone show that the condenstion process which forms a dense phase (``clouds'') is highly dynamical, and that the boundaries of the clouds are accretion shocks, rather than static density discontinuities. The density histograms (PDFs) of these runs exhibit either bimodal shapes or a single peak at low densities plus a slope change at high densities. Final static situations may be established, but the equilibrium is very fragile: small density fluctuations in the warm phase require large variations in the density of the cold phase, probably inducing shocks into the clouds. This result suggests that such configurations are highly unlikely. Simulations including turbulent forcing show that large- scale forcing is incapable of erasing the signature of the TI in the density PDFs, but small-scale, stellar-like forcing causes erasure of the signature of the instability. However, these simulations do not reach stationary regimes, TI driving an ever-increasing star formation rate. Simulations including magnetic fields, self-gravity and the Coriolis force show no significant difference between the PDFs of stable and unstable cases, and reach stationary regimes, suggesting that the combination of the stellar forcing and the extra effective pressure provided by the magnetic field and the Coriolis force overwhelm TI as a density-structure generator in the ISM. We emphasize that a multi-modal temperature PDF is not necessarily an indication of a multi-phase medium, which must contain clearly distinct thermal equilibrium phases.

We study the effects of projection of three-dimensional (3D) data onto the plane of the sky by means of numerical simulations of turbulence in the interstellar medium including the magnetic field, parameterized cooling and diffuse and stellar heating, self-gravity and rotation. We compare the physical-space density and velocity distributions with their representation in position-position-velocity (PPV) space (``channel maps''), noting that the latter can be interpreted in two ways: either as maps of the column density's spatial distribution (at a given line-of-sight (LOS) velocity), or as maps of the spatial distribution of a given value of the LOS velocity (weighted by density). This ambivalence appears related to the fact that the spatial and PPV representations of the data give significantly different views. First, the morphology in the channel maps more closely resembles that of the spatial distribution of the LOS velocity component than that of the density field, as measured by pixel-to-pixel correlations between images. Second, the channel maps contain more small-scale structure than 3D slices of the density and velocity fields, a fact evident both in subjective appearance and in the power spectra of the images. This effect may be due to a pseudo-random sampling (along the LOS) of the gas contributing to the structure in a channel map: the positions sampled along the LOS (chosen by their LOS velocity) may vary significantly from one position in the channel map to the next.