Search

Hub-filament systems are suggested to be the birth cradles of high-mass stars and clusters. We apply the FILFINDER algorithm to the integrated intensity maps of the 13CO (3-2) line to identify filaments in the G333 complex, and extract the velocity and intensity along the filament skeleton from moment maps. Clear velocity and density fluctuations are seen along the filaments, allowing us to fit velocity gradients around the intensity peaks. The velocity gradients fitted to the LAsMA data and ALMA data agree with each other over the scales covered by ALMA observations in the ATOMS survey. Changes of velocity gradient with scale indicate a ''funnel'' structure of the velocity field in PPV space, indicative of a smooth, continuously increasing velocity gradient from large to small scales, and thus consistent with gravitational acceleration. The typical velocity gradient corresponding to a 1 pc scale is ~1.6km/s/pc. Assuming free-fall, we estimate a kinematic mass within 1 pc of ~1190 M$_\odot$, which is consistent with typical masses of clumps in the ATLASGAL survey. We find direct evidence for gravitational acceleration from comparison of the observed accelerations to those predicted by free-fall onto dense hubs. On large scales, we find that the inflow may be driven by the larger scale structure, consistent with hierarchical structure in the molecular cloud and gas inflow from large to small scales. The hub-filament structures at different scales may be organized into a hierarchical system extending up to the largest scales probed, through the coupling of gravitational centers at different scales. We argue that the ''funnel'' structure in PPV space can be an effective probe for the gravitational collapse motions in molecular clouds. The large scale gas inflow is driven by gravity, implying that the molecular clouds in G333 complex may be in the state of global gravitational collapse.

Whether ionization feedback triggers the formation of massive stars is highly debated. Using ALMA 3 mm observations with a spatial resolution of $\sim 0.05$ pc and a mass sensitivity of 1.1 $\rm M_\odot$ beam$^{-1}$ at 20 K, we investigate the star formation and gas flow structures within the ionizing feedback-driven structure, a clump-scale massive ($\gtrsim 1500$ $\rm M_\odot$) bright-rimmed cloud (BRC) associated with IRAS 18290-0924. This BRC is bound only if external compression from ionized gas is considered. A small-scale ($\lesssim1$ pc) age sequence along the direction of ionizing radiation is revealed for the embedded cores and protostars, which suggests triggered star formation via radiation-driven implosion (RDI). Furthermore, filamentary gas structures converge towards the cores located in the BRC's center, indicating that these filaments are fueling mass towards cores. The local core-scale mass infall rate derived from H$^{13}$CO$^+$ $J=1-0$ blue profile is of the same order of magnitude as the filamentary mass inflow rate, approximately 1 $\rm M_\odot$ kyr$^{-1}$. A photodissociation region (PDR) covering the irradiated clump surface is detected in several molecules, such as CCH, HCO$^+$, and CS whereas the spatial distribution stratification of these molecules is indistinct. CCH spectra of the PDR possibly indicate a photoevaporation flow leaving the clump surface with a projected velocity of $\sim2$ km s$^{-1}$. Our new observations show that RDI accompanied by a clump-fed process is operating in this massive BRC. Whether this combined process works in other massive BRCs is worth exploring with dedicated surveys., Comment: 34 pages with 20 figures, Accepted by MNRAS on 2022 December 28

We present the most sensitive and detailed view of the neutral hydrogen (HI) emission associated with the Small Magellanic Cloud (SMC), through the combination of data from the Australian Square Kilometre Array Pathfinder (ASKAP) and Parkes (Murriyang), as part of the Galactic Australian Square Kilometre Array Pathfinder (GASKAP) pilot survey. These GASKAP-HI pilot observations, for the first time, reveal HI in the SMC on similar physical scales as other important tracers of the interstellar medium, such as molecular gas and dust. The resultant image cube possesses an rms noise level of 1.1 K (1.6 mJy/beam) per 0.98 km s$^{-1}$ spectral channel with an angular resolution of 30$''$ ($\sim$10 pc). We discuss the calibration scheme and the custom imaging pipeline that utilizes a joint deconvolution approach, efficiently distributed across a computing cluster, to accurately recover the emission extending across the entire $\sim$25 deg$^2$ field-of-view. We provide an overview of the data products and characterize several aspects including the noise properties as a function of angular resolution and the represented spatial scales by deriving the global transfer function over the full spectral range. A preliminary spatial power spectrum analysis on individual spectral channels reveals that the power-law nature of the density distribution extends down to scales of 10 pc. We highlight the scientific potential of these data by comparing the properties of an outflowing high velocity cloud with previous ASKAP+Parkes HI test observations., Comment: Accepted for publication in PASA, 34 pages, 18 figures, 5 tables

We investigate the origin of the observed scaling $j\sim R^{3/2}$ between the specific angular momentum $j$ and the radius $R$ of molecular clouds (MCs) and their their substructures, and of the observed near independence of $\beta$, the ratio of rotational to gravitational energy, from $R$. To this end, we measure the angular momentum (AM) of sets of particles in an SPH simulation of the formation, collapse and fragmentation of giant MCs. The sets of SPH particles are defined either as ``clumps'' (connected particle sets), or as lagrangian sets that conform a connected clump only at a certain time $\tdef$. We find that: {\it i)} Clumps evolve along the observed \jR\ relation at all times, {\it ii)} Lagrangian particle sets evolve along the observed relation when the volume containing them also contains a large number of other ``intruder'' particles. Otherwise, they evolve with $j\sim$ cst. {\it iii)} Tracking lagrangian sets to the future, we find that a subset of the SPH particles participates in the collapse, while another disperses away. {\it iv)} Noting that, under AM conservation, $\beta$ increases during contraction, we suggest that its near independence of radius may arise from the competition between this increase and an increase in the AM exchange rate at higher rotational energy density gradient. {\it v)} We suggest that, if MCs are globally dominated by gravity, the observed \jR\ relation arises because the observational selection of its dense structures amounts to selecting the fragments that have lost AM via Reynolds stresses from their neighbors., Comment: 20 pages, 15 figures, submitted to ApJ, comments welcome

Context: Atomic and molecular cloud formation is a dynamical process. However, kinematic signatures of these processes are still observationally poorly constrained. Methods: Targeting the cloud-scale environment of the prototypical infrared dark cloud G28.3, we employ spectral line imaging observations of the two atomic lines HI and [CI] as well as molecular lines observations in 13CO in the 1--0 and 3--2 transitions. The analysis comprises investigations of the kinematic properties of the different tracers, estimates of the mass flow rates, velocity structure functions, a Histogram of Oriented Gradients (HOG) study as well as comparisons to simulations. Results: The central IRDC is embedded in a more diffuse envelope of cold neutral medium (CNM) traced by HI self-absorption (HISA) and molecular gas. The spectral line data as well as the HOG and structure function analysis indicate a possible kinematic decoupling of the HI from the other gas compounds. Spectral analysis and position-velocity diagrams reveal two velocity components that converge at the position of the IRDC. Estimated mass flow rates appear rather constant from the cloud edge toward the center. The velocity structure function analysis is consistent with gas flows being dominated by the formation of hierarchical structures. Conclusions: The observations and analysis are consistent with a picture where the IRDC G28 is formed at the center of two converging gas flows. While the approximately constant mass flow rates are consistent with a self-similar, gravitationally driven collapse of the cloud, external compression by, e.g., spiral arm shocks or supernovae explosions cannot be excluded yet. Future investigations should aim at differentiating the origin of such converging gas flows., Comment: 13 pages, 15 figures, Astronomy & Astrophysics in press, a higher resolution version can be found at http://www.mpia.de/homes/beuther/papers.html

Observations of molecular clouds reveal a complex structure, with gas and dust often arranged in filamentary rather than spherical geometries. The associations of pre- and proto- stellar cores with the filaments suggest a direct link with the process of star formation. Any study of the properties of such filaments requires a representative samples from different enviroments and so an unbiased detection method. We developed such an approach using the Hessian matrix of a surface-brightness distribution to identify filaments and determine their physical and morphological properties. After testing the method on simulated, but realistic filaments, we apply the algorithms to column-density maps computed from Herschel observations of the Galactic Plane obtained by the Hi-GAL project. We identified ~500 filaments, in the longitude range of l=216.5 to l=225.5, with lengths from ~1 pc up to ~30 pc and widths between 0.1 pc and 2.5 pc. Average column densities are between 10^20 cm^-2 and 10^22 cm^-2. Filaments include the majority of dense material with N_H_2 > 6x10^21cm^-2. We find that the pre- and proto-stellar compact sources already identified in the same region are mostly associated with filaments. However, surface densities in excess of the expected critical values for high-mass star formation are only found on the filaments, indicating that these structures are necessary to channel material into the clumps. Furthermore, we analyze the gravitational stability of filaments and discuss their relationship with star formation., Comment: 32 pages, 23 Figures. Accepted for publication in ApJ

In this contribution, we test our previously published one-dimensional PDR model for deriving total hydrogen volume densities from HI column density measurements in extragalactic regions by applying it to the Taurus molecular cloud, where its predictions can be compared to available data. Also, we make the first direct detailed comparison of our model to CO(1-0) and far-infrared emission. Using an incident UV flux G0 of 4.25 ({\chi} = 5) throughout the main body of the cloud, we derive total hydrogen volume densities of \approx 430 cm-3, consistent with the extensive literature available on Taurus. The distribution of the volume densities shows a log-normal shape with a hint of a power-law shape on the high density end. We convert our volume densities to H2 column densities assuming a cloud depth of 5 parsec and compare these column densities to observed CO emission. We find a slope equivalent to a CO conversion factor relation that is on the low end of reported values for this factor in the literature (0.9 x 1020 cm-2 (K km s-1)-1), although this value is directly proportional to our assumed value of G0 as well as the cloud depth. We seem to under-predict the total hydrogen gas as compared to 100 {\mu}m dust emission, which we speculate may be caused by a higher actual G0 incident on the Taurus cloud than is generally assumed., Comment: 14 pages, 14 figures. Accepted for publication in MNRAS

We present the results of radiation-magnetohydrodynamic simulations of the expansion of HII regions and surrounding photodissociation regions in turbulent, magnetised, molecular clouds on scales of up to 4 parsecs, including the effects of ionising and non-ionising ultraviolet radiation and x rays from young star clusters. We find that HII region expansion reduces the disordered component of the B field, imposing a large-scale order on the field around its border. The field in the neutral gas lies along the ionisation front, while the field in the ionised gas tends to be perpendicular to this. The highest pressure compressed neutral/molecular gas is driven towards approximate equipartition between thermal/magnetic/turbulent energy densities, whereas lower pressure neutral/molecular gas divides into quiescent, magnetically dominated regions, and, on the other hand, turbulent, demagnetised regions. The ionised gas shows approximate thermal/turbulent equipartition, but with magnetic energy densities 1 to 3 orders of magnitude lower. A high velocity dispersion (approx 8 km/s) is maintained in the ionised gas throughout our simulations, despite the mean expansion velocity being significantly lower. The B field does not significantly brake the HII region expansion on the length and timescales accessible to our simulations, but it does tend to suppress the small-scale fragmentation and radiation-driven implosion of neutral/molecular gas that forms globules and pillars at the edge of the HII region. However, the relative luminosity of ionising and non-ionising radiation has a much larger influence than the presence or absence of the B field. When the radiation field is relatively soft (as in the case of a lower mass cluster, with earliest spectral type of B0.5), then fragmentation is less vigorous and a thick, relatively smooth PDR forms. Movies available at http://youtube.com/user/divBequals0, Comment: 23 pages, 23 figures (11 in color), accepted by MNRAS. Minor revisions to match accepted version. Animations at http://youtube.com/user/divBequals0

{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&A

(Abridged) The formation of molecular clouds (MCs) from the diffuse interstellar gas appears to be a necessary step for star formation, as young stars invariably occur within them. However, the mechanisms controlling the formation of MCs remain controversial. In this contribution, we focus on their formation in compressive flows driven by interstellar turbulence and large-scale gravitational instability. Turbulent compression driven by supernovae appears insufficient to explain the bulk of cloud and star formation. Rather, gravity must be important at all scales, driving the compressive flows that form both clouds and cores. Cooling and thermal instability allow the formation of dense gas out of moderate, transonic compressions in the warm diffuse gas, and drive turbulence into the dense clouds. MCs may be produced by an overshoot beyond the thermal-pressure equilibrium between the cold and warm phases of atomic gas, caused by some combination of the ram pressure of compression and the self-gravity of the compressed gas. In this case, properties of the clouds such as their mass, mass-to-magnetic flux ratio, and total kinetic and gravitational energies are in general time-variable quantities. MCs may never enter a quasi-equilibrium or virial equilibrium state but rather continuously collapse to stars., Comment: Proceeding of conference "Structure formation in the universe", held in Chamonix 2007. To be published in Structure formation in Astrophysics, Ed. G. Chabrier edited by Cambridge University Press, 2008

We study the formation of giant dense cloud complexes and of stars within them by means of SPH numerical simulations of the mildly supersonic collision of gas streams (``inflows'') in the warm neutral medium (WNM). The resulting compressions cause cooling and turbulence generation in the gas, forming a cloud that then becomes self-gravitating and undergoes global collapse. Simultaneously, the turbulent, nonlinear density fluctuations induce fast, local collapse events. The simulations show that: a) The clouds are not in a state of equilibrium. Instead, they undergo secular evolution. Initially, their mass and gravitational energy |Eg| increase steadily, while the turbulent energy Ek reaches a plateau. b) When |Eg| becomes comparable to Ek, global collapse begins, causing a simultaneous increase in both that maintains a near-equipartition condition |Eg| ~ 2 Ek. c) Longer inflow durations delay the onset of global and local collapse, by maintaining a higher turbulent velocity dispersion in the cloud over longer times. d) The star formation rate is large from the beginning, without any period of slow and accelerating star formation. e) The column densities of the local star-forming clumps are very similar to reported values of the column density required for molecule formation, suggesting that locally molecular gas and star formation occur nearly simultaneously. The MC formation mechanism discussed here naturally explains the apparent ``virialized'' state of MCs and the ubiquitous presence of HI halos around them. Within their assumptions, our simulations support the scenario of rapid star formation after MCs are formed, although long (>~ 15 Myr) accumulation periods do occur during which the clouds build up their gravitational energy, and which are expected to be spent in the atomic phase., Comment: Abstract matches version to be published. Corrected link to movies: http://www.astrosmo.unam.mx/~e.vazquez/turbulence/movies/VSetal06/VSetal06.html

We review the properties of turbulent molecular clouds (MCs), focusing on the physical processes that influence star formation (SF). MC formation appears to occur during large-scale compression of the diffuse ISM driven by supernovae, magnetorotational instability, or gravitational instability in galactic disks of stars and gas. The compressions generate turbulence that can accelerate molecule production and produce the observed morphology. We then review the properties of MC turbulence, including density enhancements observed as clumps and cores, magnetic field structure, driving scales, the relation to observed scaling relations, and the interaction with gas thermodynamics. We argue that MC cores are dynamical, not quasistatic, objects with relatively short lifetimes not exceeding a few megayears. We review their morphology, magnetic fields, density and velocity profiles, and virial budget. Next, we discuss how MC turbulence controls SF. On global scales turbulence prevents monolithic collapse of the clouds; on small scales it promotes local collapse. We discuss its effects on the SF efficiency, and critically examine the possible relation between the clump mass distribution and the initial mass function, and then turn to the redistribution of angular momentum during collapse and how it determines the multiplicity of stellar systems. Finally, we discuss the importance of dynamical interactions between protostars in dense clusters, and the effect of the ionization and winds from those protostars on the surrounding cloud. We conclude that the interaction of self-gravity and turbulence controls MC formation and behavior, as well as the core and star formation processes within them., Comment: Review chapter in "Protostars and Planets V". 18 pages, 4 figures

We discuss the lifetimes and evolution of dense cores formed as turbulent density fluctuations in magnetized, isothermal molecular clouds. We consider numerical simulations in which we measure the cores' magnetic criticality and Jeans stability in relation to the magnetic criticality of their ``parent clouds'' (the numerical boxes). In subcritical boxes, dense cores do not form, and collapse does not occur. In supercritical boxes, some cores collapse, being part of larger clumps that are supercritical from the start, and whose central, densest regions (the cores) are initially subcritical, but rapidly become supercritical, presumably by accretion along field lines. Numerical artifacts are ruled out. The time scales for cores to go from subcritical to supercritical and then collapse are a few times the core free-fall time, $\tfc$. Our results suggest that cores are out-of-equilibrium, transient structures, rather than quasi-magnetostatic configurations., Comment: Contributed paper for the proceedings of "Gravitational Collapse: From Massive Stars to Planets", a meeting to honor Peter Bodenheimer, ed. G. Garcia-Segura. To appear in Rev. Mex. Astron. Ap., Ser. Conf. Uses rmaa.cls. 5 pages, 3 figures. Animations can be downloaded from http://www.astrosmo.unam.mx/~e.vazquez/turbulence_HP/movies/VKSB04.html

The kinetic energy dissipation rate in the turbulent ISM of disk galaxies is a key ingredient in galaxy evolution models since it determines the effectiveness of large-scale star formation (SF) feedback. Using magneto-hydro-dynamic simulations, we find that the ISM dissipates efficiently the turbulent kinetic energy injected by sources of stellar nature. Thus, the SF process may be self-regulated by an energy balance only at the level of the disk ISM. The use of the self-regulation SF mechanism in galaxy evolutionary models, where disks form inside growing Cold Dark Matter halos, allows to predict the SF history of disk galaxies, including the Milky Way and the solar neighborhood, as well as the contribution of the whole population of disk galaxies to the cosmic SF history. The results are encouraging., Comment: 5 pages, 2 figures. To appear in "Galaxy Evolution: Theory and Observations", eds. V. Avila-Reese, C. Firmani, C. Frenk, & C. Allen, RevMexA&A (SC), v. 17 Small changes in Fig. 2 and in subsection 3.2

We compile a holistic scenario for molecular cloud (MC) evolution and control of the star formation efficiency (SFE), and present a first set of numerical tests of it. A {\it lossy} compressible cascade can generate density fluctuations and further turbulence at small scales from large-scale motions, implying that the turbulence in MCs may originate from the compressions that form them. Below a {\it sonic} scale $\ls$, turbulence cannot induce any further subfragmentation, nor be a dominant support agent against gravity. Since progressively smaller density peaks contain progressively smaller fractions of the mass, we expect the SFE to decrease with decreasing $\ls$, at least when the cloud is globally supported by turbulence. Our numerical experiments confirm this prediction. We also find that the collapsed mass fraction in the simulations always saturates below 100% efficiency. This may be due to the decreased mean density of the leftover interclump medium, which in real clouds (not confined to a box) should then be more easily dispersed, marking the ``death'' of the cloud. We identify two different functional dependences (``modes'') of the SFE on $\ls$, which roughly correspond to globally supported and unsupported cases. Globally supported runs with most of the turbulent energy at the largest scales have similar SFEs to those of unsupported runs, providing numerical evidence of the dual role of turbulence, whereby large-scale turbulent modes induce collapse at smaller scales. We tentatively suggest that these modes may correspond to the clustered and isolated modes of star formation, although here they are seen to form part of a continuum rather than being separate modes. Finally, we compare with previous proposals that the relevant parameter is the energy injection scale., Comment: 6 pages, 3 figures. Uses emulateapj. Accepted in ApJ Letters

(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., Comment: 29 pages, 13 figures, accepted by ApJ, minor changes included

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., Comment: video address has been updated

We discuss the dissipation of turbulent kinetic energy Ek in the global ISM by means of 2-D, MHD, non-isothermal simulations in the presence of model radiative heating and cooling. We argue that dissipation in 2D is representative of that in three dimensions as long as it is dominated by shocks rather than by a turbulent cascade. Energy is injected at a few isolated sites in space, over relatively small scales, and over short time periods. This leads to the coexistence of forced and decaying regimes in the same flow. We find that the ISM-like flow dissipates its turbulent energy rapidly. In simulations with forcing, the input parameters are the radius l_f of the forcing region, the total kinetic energy e_k each source deposits into the flow, and the rate of formation of those regions, sfr_OB. The global dissipation time t_d depends mainly on l_f. In terms of measurable properties of the ISM, t_d >= Sigma_g u_rms^2/(e_k sfr_OB), where Sigma_g is the average gas surface density and u_rms is the rms velocity dispersion. For the solar neighborhood, t_d >= 1.5x10^7 yr. The global dissipation time is consistently smaller than the crossing time of the largest energy-containing scales. In decaying simulations, Ek decreases with time as t^-n, where n~0.8-0.9. This suggests a decay with distance d as Ek\propto d^{-2n/(2-n)} in the mixed forced+decaying case. If applicable to the vertical direction, our results support models of galaxy evolution in which stellar energy injection provides significant support for the gas disk thickness, but not models of galaxy formation in which this energy injection is supposed to reheat an intra-halo medium at distances of up to 10-20 times the optical galaxy size, as the dissipation occurs on distances comparable to the disk height., Comment: 23 pages, including figures. To appear in ApJ. Abstract abridged

We study turbulent dissipation in the ISM and explore some implications for galaxy formation and evolution using 2D MHD numerical simulations of compressible fluids. The turbulent kinetic energy E_k is injected by stellar sources formed self-consistently in the simulation. In the ISM-like fluid, regimes of both forced and decaying turbulence coexist. In the active turbulent regions (forced regime), E_k is dissipated locally and efficiently. In the decaying regime (far from input sources), E_k(t) decays ~(1+t)^-0.8. The residual turbulent motions may propagate distances of the order of the observed disk height, suggesting that turbulence may be the responsible of vertical support and star formation self-regulation at the disk level, but not at the level of the whole cosmological halo, as would be required in some models of galaxy formation., Comment: 3 pages, uses rmaa.cls, to appear in "Astrophysical Plasmas: Codes, Models, and Observations", Eds. J. Franco, J. Arthur, N. Brickhouse, Rev.Mex.AA Conf. Series. An error in the calculation of the dissipation time, t_d, was corrected (it referred to the dissipation in the forced regions). Now we report the global t_d, which is the relevant quantity for the whole flow

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., Comment: 35 pages, 1 ps figure, 1 style file. Invited review for 'Protostars and Planets IV', eds. V. Mannings, A. Boss, S. Russell

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., Comment: Submitted to ApJ. 23 pages, 15 postscript figures

We present numerical simulations of an isothermal turbulent gas undergoing gravitational collapse, aimed at testing for ``logatropic'' behavior of the form $P_t \sim \log \rho$, where $P_t$ is the ``turbulent pressure'' and $\rho$ is the density. To this end, we monitor the evolution of the turbulent velocity dispersion $\sigma$ as the density increases during the collapse. A logatropic behavior would require that $\sigma \propto \rho^{-1/2}$, a result which, however, is not verified in the simulations. Instead, the velocity dispersion increases with density, implying a polytropic behavior of $P_t$. This behavior is found both in purely hydrodynamic as well as hydromagnetic runs. For purely hydrodynamic and rapidly-collapsing magnetic cases, the velocity dispersion increases roughly as $\sigma \propto \rho^{1/2}$, implying $P_t\sim \rho^2$, where $P_t$ is the turbulent pressure. For slowly-collapsing magnetic cases the behavior is close to $\sigma \propto \rho^{1/4}$, which implies $P_t \sim \rho^{3/2}$. We thus suggest that the logatropic ``equation of state'' may represent only the statistically most probable state of an ensemble of clouds in equilibrium between self-gravity and kinetic support, but does not adequately represent the behavior of the ``turbulent pressure'' within a cloud undergoing a dynamic compression due to gravitational collapse. Finally, we discuss the importance of the underlying physical model for the clouds (in equilibrium vs. dynamic) on the results obtained., Comment: Accepted in ApJ. 10 pages, 3 postscript figures

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., Comment: 4 pages, 3 postscript figures. Latex, uses aipproc.sty Contribution to the Conference Proc. of the 7th Annual Astrophysics Conference in Maryland, STAR FORMATION, NEAR AND FAR, eds. Stephen S. Holt and Lee G. Mundy

We present results from a statistical study of clouds in two-dimensional numerical simulations of the interstellar medium. The clouds in the simulations exhibit a differential mass spectrum $dN(M)/dM \sim M^{-1.44 \pm 0.1}$ and a velocity dispersion-size relation $\Delta v \sim R^{0.41 \pm 0.08}$. However, the clouds do {\it not} exhibit a clear density-size relation. At a given mean density, clouds span a range of sizes from the smallest resolved scales up to a maximum given by a Larson-type relation $R_{\rm max} \sim \rho^\alpha$, with $\alpha = -0.81 \pm .15$, although numerical effects cannot be ruled out as responsible for the latter correlation. Clouds span a range of column densities $N$ of two orders of magnitude, supporting the suggestion that the observational density-size relation may be an artifact of survey limitations. In this case, the $\Delta v$--$R$ relation can be interpreted as a direct consequence of a $k^{-2}$ turbulent spectrum, characteristic of a field of shocks, verified in the simulations, rather than of virial equilibrium of clouds with a $\rho \propto R^{-1}$ law. However, we also discuss the possibility that the clouds are in balance between self-gravity and turbulence, but with a scatter of at least a factor of 10 in the $Delta v$--$R$ relation, and of 100 in the density-size relation, according to the equilibrium relation $\Delta v \sim (NR)^{1/2}$. Finally, we compare these results with observational data. We propose a simple model suggesting that recent results finding nearly constant column densities for dark IRAS clouds may be an artifact of a temperature gradient within the clouds induced by external radiative heating. As a consequence, we emphasize that IRAS surface brightness maps are not appropriate for measuring column densities., Comment: Uuencoded, gzipped, tar file with Latex file + 21 postscript figures (20 pages). Gzipped Postscript file with figures (500 KB) available at ftp://kepler.astroscu.unam.mx/incoming/enro/papers/lowdens.ps.gz Mpeg video of Run 28 available at http://www.astroscu.unam.mx/hp/Movies/den.mpg. Accepted in ApJ

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., Comment: 20 pages, 12 Postscript figures. Uses aas2pp4.sty. Accepted in ApJ

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., Comment: Gzipped, tarred Latex file (7 pages), and style file. Uses AAS macros. Gzipped Postscript file with 9 figures available at ftp://kepler.astroscu.unam.mx/incoming/enro/papers/mhdcfsf.ps.gz

We present preliminary results on the energy budgets of clouds in two-dimensional numerical simulations of the interstellar medium. Using an automated cloud-identification algorithm, we calculate the gravitational, internal, kinetic and magnetic energies of the clouds. We find that, within a dispersion of roughly one order of magnitude, the gravitational energy in the clouds is balanced by the remaining energies. Furthermore, within the same dispersion, there appears to be equipartition between the kinetic and magnetic energies., Comment: Gzipped, tarred Latex file (4 pages), 6 Postscript figures and one style file. Uses AAS macros. Gzipped Postscript file also available at ftp://kepler.astroscu.unam.mx/incoming/enro/papers/cloudstat.ps.gz

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., Comment: Plain TeX file, 25 pages. Gzipped, tarred set of Tex file plus 17 figures and 3 tables (Postscript) available at ftp://kepler.astroscu.unam.mx/incoming/enro/papers/mhdgturb.tar.gz