A pair of nonlinear diffusion equations in Fourier space} is used to study the dynamics of strong Alfv\'en-wave turbulence, from MHD to electron scales. Special attention is paid to the regime of imbalance between the energies of counter-propagating waves commonly observed in the solar wind (SW), especially in regions relatively close to the Sun. In the collisionless regime where dispersive effects arise at scales comparable to or larger than those where dissipation becomes effective, the imbalance produced by a given injection rate of generalized cross-helicity (GCH), which is an invariant, is much larger than in the corresponding collisional regime described by the usual (or reduced) magnetohydrodynamics. The combined effect of high imbalance and ion Landau damping induces a steep energy spectrum for the transverse magnetic field at sub-ion scales. This spectrum is consistent with observations in highly Alfvenic regions of the SW, such as trailing edges, but does not take the form of a transition range continued at smaller scales by a shallower spectrum. This suggests that the observed spectra displaying such a transition result from the superposition of contributions originating from various streams with different degrees of imbalance. Furthermore, when imbalanced energy injection is supplemented at small scales in an already fully developed turbulence, for example under the effect of magnetic reconnection, a significant enhancement of the imbalance at all scales is observed.

A two-field Hamiltonian gyrofluid model for kinetic Alfv\'en waves retaining ion finite Larmor radius corrections, parallel magnetic field fluctuations and electron inertia, is used to study turbulent cascades from the MHD to the sub-ion scales. Special attention is paid to the case of imbalance between waves propagating along or opposite to the ambient magnetic field. For weak turbulence in the absence of electron inertia, kinetic equations for the spectral density of the conserved quantities (total energy and generalized cross-helicity) are obtained. They provide a global description, matching between the regimes of reduced MHD at large scales and electron reduced MHD at small scales, previously considered in the literature. In the limit of ultra-local interactions, Leith-type nonlinear diffusion equations in the Fourier space are derived and heuristically extended to the strong turbulence regime by modifying the transfer time appropriately. Relations with existing phenomenological models for imbalanced MHD and balanced sub-ion turbulence are discussed. It turns out that in the presence of dispersive effects, the dynamics is sensitive on the way turbulence is maintained in a steady state. Furthermore, the total energy spectrum at sub-ion scales becomes steeper as the generalized cross-helicity flux is increased.

Reduced fluid models including electron inertia and ion finite Larmor radius corrections are derived asymptotically, both from fluid basic equations and from a gyrofluid model. They apply to collisionless plasmas with small ion-to-electron equilibrium temperature ratio and low $\beta_e$, where $\beta_e$ indicates the ratio between the equilibrium electron pressure and the magnetic pressure exerted by a strong, constant and uniform magnetic guide field. The consistency between the fluid and gyrofluid approaches is ensured when choosing ion closure relations prescribed by the underlying ordering. A two-field reduction of the gyrofluid model valid for arbitrary equilibrium temperature ratio is also introduced, and is shown to have a noncanonical Hamiltonian structure. This model provides a convenient framework for studying kinetic Alfv\'en wave turbulence, from MHD to sub-$d_e$ scales (where $d_e$ holds for the electron skin depth). Magnetic energy spectra are phenomenologically determined within energy and generalized helicity cascades in the perpendicular spectral plane. Arguments based on absolute statistical equilibria are used to predict the direction of the transfers, pointing out that, within the sub-ion range associated with a $k_\perp^{-7/3}$ transverse magnetic spectrum, the generalized helicity could display an inverse cascade if injected at small scales, for example by reconnection processes.

The cascade of kinetic Alfv\'en waves (KAWs) at the sub-ion scales in the solar wind is numerically simulated using a fluid approach that retains ion and electron Landau damping, together with ion finite Larmor radius corrections. Assuming initially equal and isotropic ion and electron temperatures, and an ion beta equal to unity, different simulations are performed by varying the propagation direction and the amplitude of KAWs that are randomly driven at a transverse scale of about one fifth of the proton gyroradius in order to maintain a prescribed level of turbulent fluctuations. The resulting turbulent regimes are characterized by the nonlinearity parameter, defined as the ratio of the characteristic times of Alfv\'en wave propagation and of the transverse nonlinear dynamics. The corresponding transverse magnetic energy spectra display power laws with exponents spanning a range of values consistent with spacecraft observations. The meandering of the magnetic field lines together with the ion temperature homogenization along these lines are shown to be related to the strength of the turbulence, measured by the nonlinearity parameter. The results are interpreted in terms of a recently proposed phenomenological model where the homogenization process along field lines induced by Landau damping plays a central role.

A phenomenological turbulence model for kinetic Alfv\'en waves in a magnetizedcollisionless plasma, able to reproduce the non-universalpower-law spectra observed at the sub-ion scales in the solar wind and the terrestrial magnetosphere, is presented.The process of temperature homogenization along distortedmagnetic field lines, induced by Landau damping,affects the turbulencetransfer time and results in a steepening of the sub-ion power-law spectrumof critically-balanced turbulence, whose exponent is sensitive to the ratio between the Alfv\'en wave period and the nonlinear timescale. Transition from large-scaleweak turbulence to smaller scale strong turbulence is capturedand non local interactions, relevant in the case of steep spectra, are accounted for.

Numerical simulations of dispersive turbulence in magnetized plasmas based on the Hall-MHD description are presented, assuming spatial variations along a unique direction making a prescribed angle with the ambient magnetic field. Main observations concern the energy transfers among the different scales and the various types of MHD waves, together with the conditions for the establishment of pressure-balanced structures. For parallel propagation, Alfv\'en-wave transfer to small scales is strongly inhibited and rather feeds magnetosonic modes, unless the effect of dispersion is strong enough at the energy injection scale. In oblique directions, the dominantly compressible character of the turbulence is pointed out with, for quasi-transverse propagation, the presence of conspicuous kinetic Alfv\'en waves. Preliminary simulations of a Landau fluid model incorporating relevant linear kinetic effects reveal the development of a significant plasma temperature anisotropy leading to recurrent instabilities.

The neutral atomic interstellar medium (HI) is well known to be strongly magnetized. Since HI is a 2-phase medium the questions arise of what is the effect of the magnetic field on a 2-phase medium, how magnetic field affects thermal instability ? Here we summarize analytical and numerical results which have been obtained previously to investigate this question.

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.

The effect of Alfv\'en waves on the thermal instability of the Interstellar Medium (ISM) is investigated both analytically and numerically. A stability analysis of a finite amplitude circularly polarized Alfv\'en wave propagating parallel to an ambient magnetic field in a thermally unstable gas at thermal equilibrium is performed, leading to a dispersion relation which depends on 3 parameters, namely the square ratio of the sonic and Alfv\'en velocities ($\beta$), the wave amplitude and the ratio between the wave temporal period and the cooling time. Depending on the values of these 3 parameters, the Alfv\'en waves can stabilize the large-scale perturbations, destabilize those whose wavelength is a few times the Alfv\'en wavelength $\lambda_{AW}$ or leave the growth rate of the short scales unchanged. In order to investigate the non-linear regime, two different numerical experiments are performed in a slab geometry. (TBC)

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 discuss magnetic pressure and density fluctuations in strongly turbulent isothermal MHD flows in "1+2/3" dimensions. We first consider "simple" nonlinear MHD waves, which allow us show that the slow and fast modes have different asymptotic dependences of the magnetic pressure B^2 vs. rho. For the slow mode, B^2 ~= c_1-c_2 rho, while for the fast mode, B^2 ~= rho^2. We also perform a perturbative analysis to investigate Alfven wave pressure, recovering previous results that B^2 ~= rho^gamma_e, with gamma_e ~= 2, 3/2 and 1/2 at respectively large, moderate and low M_a. This variety of scalings implies that a single polytropic description of magnetic pressure is not possible in general, since the relation between B^2 and rho depends on which mode dominates the density fluctuation production, which in turn depends on the angle between the magnetic field and the direction of wave propagation, and on the Alfvenic Mach number M_a. Typically, at small M_a, the slow mode dominates, and B is ANTIcorrelated with rho. At large M_a, both modes contribute to density fluctuation production, and the magnetic pressure decorrelates from density, exhibiting a large scatter, which however decreases towards higher densities. In this case, the unsystematic behavior of the magnetic pressure causes the density PDF to generally maintain the lognormal shape corresponding to non-magnetic isothermal turbulence, except when the slow mode dominates, in which case the PDF develops an excess at low densities. Our results are consistent with the low values and apparent lack of correlation between the magnetic field strength and density in surveys of the lower-density molecular gas, and also with the recorrelation apparently seen at higher densities, if M_a is relatively large there.

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.

We investigate the behavior of the magnetic pressure, $b^2$, in fully turbulent MHD flows in ``1+2/3'' dimensions by means of its effect on the probability density function (PDF) of the density field. We start by reviewing our previous results for general polytropic flows, according to which the value of the polytropic exponent $\gamma$ determines the functional shape of the PDF. A lognormal density PDF appears in the isothermal ($\gamma=1$) case, but a power-law tail at either large or small densities appears for large Mach numbers when $\gamma >1$ and $\gamma < 1$, respectively. In the isothermal magnetic case, the relevant parameter is the field fluctuation amplitude, $\dbb$. A lognormal PDF still appears for small field fluctuations (generally the case for {\it large mean fields}), but a significant low-density excess appears at large fluctuation amplitudes ({\it weak mean fields}), similar to the behavior at $\gamma > 1$ of polytropic flows. We interpret these results in terms of simple nonlinear MHD waves, for which the magnetic pressure behaves linearly with the density in the case of the slow mode, and quadratically in the case of the fast wave. Finally, we discuss some implications of these results, in particular the fact that the effect of the magnetic field in modifying the PDF is strongest when the mean field is weak.

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.

We review recent results from numerical simulations and related models of MHD turbulence in the interstellar medium (ISM) and in molecular clouds. We discuss the implications of turbulence for the processes of cloud formation and evolution, and the determination of clouds' physical properties. Numerical simulations of the turbulent ISM to date have included magnetic fields, self-gravity, parameterized heating and cooling, modeled star formation and other turbulent inputs. The structures which form reproduce well observed velocity-size scaling properties, while predicting the non-existence of a general density-size scaling law. Criteria for the formation of gravitationally-bound structures by turbulent compression are summarized. For flows with equations of state $P\propto \rho^\gamma$, the statistics of the density field depend on the exponent $\gamma$. Numerical simulations of both forced and decaying MHD compressible turbulence have shown that the decay rate is comparable to the non-magnetic case. For virialized clouds, the turbulent decay time is shorter than the gravitational free-fall time, so wholesale cloud collapse is only prevented by ongoing turbulent inputs and/or a strong mean magnetic field. Finally, perspectives for future work in this field are briefly discussed.

We discuss HD and MHD compressible turbulence as a cloud-forming and cloud-structuring mechanism in the ISM. Results from a numerical model of the turbulent ISM at large scales suggest that the phase-like appearance of the medium, the typical values of the densities and magnetic field strengths in the intercloud medium, as well as the velocity dispersion-size scaling relation in clouds may be understood as consequences of the interstellar turbulence. However, the density-size relation appears to only hold for the densest clouds, suggesting that low-column density clouds, which are hardest to observe, are turbulent transients. We then explore some properties of highly compressible polytropic turbulence, in one and several dimensions, applicable to molecular cloud scales. At low values of the polytropic index $\gamma$, turbulence may induce the gravitational collapse of otherwise linearly stable clouds, except if they are magnetically subcritical. The nature of the density fluctuations in the high Mach-number limit depends on $\gamma$, and in no case resembles that resulting from Burgers turbulence. In the isothermal ($\gamma=1$) case, the dispersion of $\ln (\rho)$ scales like the turbulent Mach number. The latter case is singular with a lognormal density pdf, while power-law tails develop at high (resp. low) densities for $\gamma <1$ (resp. $\gamma >1$).

We discuss the generation and statistics of the density fluctuations in highly compressible polytropic turbulence, based on a simple model and one-dimensional numerical simulations. Observing that density structures tend to form in a hierarchical manner, we assume that density fluctuations follow a random multiplicative process. When the polytropic exponent $\gamma$ is equal to unity, the local Mach number is independent of the density, and our assumption leads us to expect that the probability density function (PDF) of the density field is a lognormal. This isothermal case is found to be singular, with a dispersion $\sigma_s^2$ which scales like the square turbulent Mach number $\tilde M^2$, where $s\equiv \ln \rho$ and $\rho$ is the fluid density. This leads to much higher fluctuations than those due to shock jump relations. Extrapolating the model to the case $\gamma \not =1$, we find that, as the Mach number becomes large, the density PDF is expected to asymptotically approach a power-law regime, at high densities when $\gamma<1$, and at low densities when $\gamma>1$. This effect can be traced back to the fact that the pressure term in the momentum equation varies exponentially with $s$, thus opposing the growth of fluctuations on one side of the PDF, while being negligible on the other side. This also causes the dispersion $\sigma_s^2$ to grow more slowly than $\tilde M^2$ when $\gamma\not=1$. In view of these results, we suggest that Burgers flow is a singular case not approached by the high-$\tilde M$ limit, with a PDF that develops power laws on both sides.

We investigate the form of the one-point probability distribution function (pdf) for the density field of the interstellar medium using numerical simulations that successively reduce the number of physical processes included. Two-dimensional simulations of self-gravitating supersonic MHD and hydrodynamic turbulence, and of decaying Burgers turbulence, produce in all cases filamentary density structures and a power-law density pdf with logarithmic slope around -1.7. This suggests that the functional form of the pdf and the general filamentary morphology are the signature of the nonlinear advection operator. These results do not support previous claims that the pdf is lognormal. A series of 1D simulations of forced supersonic polytropic turbulence is used to resolve the discrepancy. They suggest that the pdf is lognormal only for effective polytropic indices $\gamma=1$ (or nearly lognormal for $\gamma\not=1$ if the Mach number is sufficiently small), while power laws develop at high densities if $\gamma<1$. We evaluate the polytropic index for conditions relevant to the cool interstellar medium using published cooling functions and different heating sources, finding that a lognormal pdf may occur at densities between 10$^3$ and at least 10$^4$ cm$^{-3}$. Several applications are examined. First, we question a recent derivation of the IMF from the density pdf by Padoan, Nordlund & Jones because a) the pdf does not contain spatial information, and b) their derivation produces the most massive stars in the voids of the density distribution. Second, we illustrate how a distribution of ambient densities can alter the predicted form of the size distribution of expanding shells. Finally, a brief comparison is made with the density pdfs found in cosmological simulations.

We investigate the properties of highly compressible turbulence and its ability to produce self-gravitating structures. The compressibility is parameterized by an effective polytropic exponent gama-eff. In the limit of small gama-eff, the density jump at shocks is shown to be of the order of e^{M^2}, and the production of vorticity by the nonlinear terms appears to be negligible. In the presence of self-gravity, we suggest that turbulence can produce bound structures for gama-eff < 2(1-1/n), where 'n' is the typical dimensionality of the turbulent compressions. We show, by means of numerical simulations, that, for sufficiently small gama-eff, small-scale turbulent density fluctuations eventually collapse even though the medium is globally stable. This result is preserved in the presence of a magnetic field for supercritical mass-to-flux ratios.

We investigate the properties of highly compressible turbulence, the compressibility arising from a small effective polytropic exponent $\gamma_e$ due to cooling. In the limit of small $\gamma_e$, the density jump at shocks is shown to be of the order of $e^{M^2}$. Without self-gravity, the density structures arising in the moderately compressible case consist mostly of patches separated by shocks and behaving like waves, while in the highly compressible case clearly defined long-lived object-like clouds emerge. When the forcing in the momentum equation is purely compressible, the rotational energy decays monotonically in time, indicating that the vortex-stretching term is not efficient in transferring energy to rotational modes. This property may be at the origin of the low amount of rotation found in interstellar clouds. Vorticity production is found to rely heavily on the presence of additional terms in the equations. In the presence of self-gravity, we suggest that turbulence can produce bound structures for $\gamma_e < 2(1-1/n)$, where $n$ is the typical dimensionality of the turbulent compressions. We support this result by means of numerical simulations in which, for sufficiently small $\gamma_e$, small-scale turbulent density fluctuations eventually collapse even though the medium is globally stable. This result is preserved in the presence of a magnetic field for supercritical mass-to-flux ratios. At larger polytropic exponents, turbulence alone is not capable of producing bound structures, and collapse can only occur when the medium is globally unstable. This mechanism is a plausible candidate for the differentiation between primordial and present-day stellar-cluster formation and for the low efficiency of star formation.

We discuss the role of turbulence in cloud and star formation, as observed in numerical simulations of the interstellar medium. Turbulent compression at the interfaces of colliding gas streams is responsible for the formation of intermediate ($\simlt 100$ pc) and small clouds (a few tens of pc), although the smallest clouds can also form from fragmentation of expanding shells around stellar heating centers. The largest cloud complexes (several hundred pc) seem to form by slow, gravitational instability-driven merging of individual clouds, which can actually be described as a large-scale tendency towards homogenization of the flow due to gravity rather than cloud collisions. These mechanisms operate as well in the presence of a magnetic field and rotation, although slight variations on the compressibility and cloud morphology are present which depend on the strength and topology of the field. In summary, the role of turbulence in the life-cycle of clouds appears to be twofold: small-scale modes contribute to cloud support, while large-scale modes can both form or disrupt clouds.

We present results from two-dimensional numerical simulations of a supersonic turbulent flow in the plane of the galactic disk, incorporating shear, thresholded and discrete star formation (SF), self-gravity, rotation and magnetic fields. A test of the model in the linear regime supports the results of the linear theory of Elmegreen (1991). In the fully nonlinear turbulent regime, while some results of the linear theory persist, new effects also emerge. Some exclusively nonlinear effects are: a) Even though there is no dynamo in 2D, the simulations are able to maintain or increase their net magnetic energy in the presence of a seed uniform azimuthal component. b) A well-defined power-law magnetic spectrum and an inverse magnetic cascade are observed in the simulations, indicating full MHD turbulence. Thus, magnetic field energy is generated in regions of SF and cascades up to the largest scales. c) The field has a slight but noticeable tendency to be aligned with density features. d) The magnetic field prevents HII regions from expanding freely, as in the recent results of Slavin \& Cox (1993). e) A tendency to exhibit {\it less} filamentary structures at stronger values of the uniform component of the magnetic field is present in several magnetic runs. f) For fiducial values of the parameters, the flow in general appears to be in rough equipartition between magnetic and kinetic energy. There is no clear domination of either the magnetic or the inertial forces. g) A median value of the magnetic field strength within clouds is $\sim 12\mu$G, while for the intercloud medium a value of $\sim 3\mu$G is found. Maximum contrasts of up to a factor of $\sim 10$ are observed.