Showing 1–50 of 88 results for author: Heitsch, F

Initial Conditions for Star Formation: A Physical Description of the Filamentary ISM

Authors: A. Hacar, S Clark, F. Heitsch, J. Kainulainen, G. Panopoulou, D. Seifried, R. Smith

Abstract: The interstellar medium contains filamentary structure over a wide range of scales. Understanding the role of this structure, both as a conduit of gas across the scales and a diagnostic tool of local physics, is a major focus of star formation studies. We review recent progress in studying filamentary structure in the ISM, interpreting its properties in terms of physical processes, and exploring f… ▽ More The interstellar medium contains filamentary structure over a wide range of scales. Understanding the role of this structure, both as a conduit of gas across the scales and a diagnostic tool of local physics, is a major focus of star formation studies. We review recent progress in studying filamentary structure in the ISM, interpreting its properties in terms of physical processes, and exploring formation and evolution scenarios. We include structures from galactic-scale filaments to tenth-of-a-parsec scale filaments, comprising both molecular and atomic structures, from both observational and theoretical perspectives. In addition to the literature overview, we assemble a large amount of catalogue data from different surveys and provide the most comprehensive census of filamentary structures to date. Our census consists of 22 803 filamentary structures, facilitating a holistic perspective and new insights. We use our census to conduct a meta-analysis, leading to a description of filament properties over four orders of magnitudes in length and eight in mass. Our analysis emphasises the hierarchical and dynamical nature of filamentary structures. Filaments do not live in isolation, nor they generally resemble static structures close to equilibrium. We propose that accretion during filament formation and evolution sets some of the key scaling properties of filaments. This highlights the role of accretion during filament formation and evolution and also in setting the initial conditions for star formation. Overall, the study of filamentary structures during the past decade has been observationally driven. While great progress has been made on measuring the basic properties of filaments, our understanding of their formation and evolution is clearly lacking. In this context, we identify a number of directions and questions we consider most pressing for the field. △ Less

Submitted 17 March, 2022; originally announced March 2022.

Comments: 36 pages, 10 figures. To appear in Protostars and Planets VII, Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura

The width of Herschel filaments varies with distance

Authors: G. V. Panopoulou, S. E. Clark, A. Hacar, F. Heitsch, J. Kainulainen, E. Ntormousi, D. Seifried, R. J. Smith

Abstract: Context: Filamentary structures in nearby molecular clouds have been found to exhibit a characteristic width of 0.1 pc, as observed in dust emission. Understanding the origin of this universal width has become a topic of central importance in the study of molecular cloud structure and the early stages of star formation. Aims: We investigate how the recovered widths of filaments depend on the dista… ▽ More Context: Filamentary structures in nearby molecular clouds have been found to exhibit a characteristic width of 0.1 pc, as observed in dust emission. Understanding the origin of this universal width has become a topic of central importance in the study of molecular cloud structure and the early stages of star formation. Aims: We investigate how the recovered widths of filaments depend on the distance from the observer, by using previously published results from the Herschel Gould Belt Survey. Methods: We obtained updated estimates on the distances to nearby molecular clouds observed with Herschel by using recent results based on 3D dust extinction mapping and Gaia. We examined the widths of filaments from individual clouds separately, as opposed to treating them as a single population. We used these per-cloud filament widths to search for signs of variation amongst the clouds of the previously published study. Results: We find a significant dependence of the mean per-cloud filament width with distance. The distribution of mean filament widths for nearby clouds is incompatible with that of farther away clouds. The mean per-cloud widths scale with distance approximately as 4--5 times the beam size. We examine the effects of resolution by performing a convergence study of a filament profile in the Herschel image of the Taurus Molecular Cloud. We find that resolution can severely affect the shapes of radial profiles over the observed range of distances. Conclusions: We conclude that the data are inconsistent with 0.1 pc being the universal characteristic width of filaments. △ Less

Submitted 8 August, 2022; v1 submitted 15 November, 2021; originally announced November 2021.

Comments: An erratum to this manuscript has been published here: https://www.aanda.org/articles/aa/full_html/2022/07/aa42281e-21/aa42281e-21.html

Journal ref: A&A 657, L13 (2022)

Mass, Morphing, Metallicities: The Evolution of Infalling High Velocity Clouds

Authors: F. Heitsch, A. Marchal, M. -A. Miville-Deschênes, J. M. Shull, A. J. Fox

Abstract: We revisit the reliability of metallicity estimates of high velocity clouds with the help of hydrodynamical simulations. We quantify the effect of accretion and viewing angle on metallicity estimates derived from absorption lines. Model parameters are chosen to provide strong lower limits on cloud contamination by ambient gas. Consistent with previous results, a cloud traveling through a stratifie… ▽ More We revisit the reliability of metallicity estimates of high velocity clouds with the help of hydrodynamical simulations. We quantify the effect of accretion and viewing angle on metallicity estimates derived from absorption lines. Model parameters are chosen to provide strong lower limits on cloud contamination by ambient gas. Consistent with previous results, a cloud traveling through a stratified halo is contaminated by ambient material to the point that <10% of its mass in neutral hydrogen consists of original cloud material. Contamination progresses nearly linearly with time, and it increases from head to tail. Therefore, metallicity estimates will depend on the evolutionary state of the cloud, and on position. While metallicities change with time by more than a factor of 10, well beyond observational uncertainties, most lines-of-sight range only within those uncertainties at any given time over all positions. Metallicity estimates vary with the cloud's inclination angle within observational uncertainties. The cloud survives the infall through the halo because ambient gas continuously condenses and cools in the cloud's wake and thus appears in the neutral phase. Therefore, the cloud observed at any fixed time is not a well-defined structure across time, since material gets constantly replaced. The thermal phases of the cloud are largely determined by the ambient pressure. Internal cloud dynamics evolve from drag gradients caused by shear instabilities, to complex patterns due to ram-pressure shielding, leading to a peloton effect, in which initially lagging gas can catch up to and even overtake the head of the cloud. △ Less

Submitted 31 October, 2021; originally announced November 2021.

Comments: 18 pages, 10 figures. Revised submission to MNRAS

Exploring the Dust Content of Galactic Halos with Herschel. IV. NGC 3079

Authors: S. Veilleux, M. Melendez, M. Stone, G. Cecil, E. Hodges-Kluck, J. Bland-Hawthorn, J. Bregman, F. Heitsch, C. L. Martin, T. Mueller, D. S. N. Rupke, E. Sturm, R. Tanner, C. Engelbracht

Abstract: We present the results from an analysis of deep Herschel far-infrared observations of the edge-on disk galaxy NGC 3079. The PSF-cleaned PACS images at 100 and 160 um display a 25 kpc x 25 kpc X-shape structure centered on the nucleus that is similar in extent and orientation to that seen in Halpha, X-rays, and the far-ultraviolet. One of the dusty filaments making up this structure is detected in… ▽ More We present the results from an analysis of deep Herschel far-infrared observations of the edge-on disk galaxy NGC 3079. The PSF-cleaned PACS images at 100 and 160 um display a 25 kpc x 25 kpc X-shape structure centered on the nucleus that is similar in extent and orientation to that seen in Halpha, X-rays, and the far-ultraviolet. One of the dusty filaments making up this structure is detected in the SPIRE 250 um map out to ~25 kpc from the nucleus. The match between the far-infrared filaments and those detected at other wavelengths suggests that the dusty material has been lifted out of the disk by the same large-scale galactic wind that has produced the other structures in this object. A closer look at the central 10 kpc x 10 kpc region provides additional support for this scenario. The dust temperatures traced by the 100-to-160 um flux ratios in this region are enhanced within a biconical region centered on the active galactic nucleus, aligned along the minor axis of the galaxy, and coincident with the well-known double-lobed cm-wave radio structure and Halpha-X-ray nuclear superbubbles. PACS imaging spectroscopy of the inner 6-kpc region reveals broad [C II] 158 um emission line profiles and OH 79 um absorption features along the minor axis of the galaxy with widths well in excess of those expected from beam smearing of the disk rotational motion. This provides compelling evidence that the cool material traced by the [C II] and OH features directly interacts with the nuclear ionized and relativistic outflows traced by the Halpha, X-ray, and radio emission. △ Less

Submitted 4 October, 2021; originally announced October 2021.

Comments: 17 pages, 12 figures, accepted for publication in MNRAS

A Co-Scaling Grid for Athena++

Authors: Roark Habegger, Fabian Heitsch

Abstract: We present a co-scaling grid formalism and its implementation in the magnetohydrodynamics code Athena++. The formalism relies on flow symmetries in astrophysical problems involving expansion, contraction, and center-of-mass motion. The grid is evolved at the same time order as the fluid variables. The user specifies grid evolution laws, which can be independent of the fluid motion. Applying our im… ▽ More We present a co-scaling grid formalism and its implementation in the magnetohydrodynamics code Athena++. The formalism relies on flow symmetries in astrophysical problems involving expansion, contraction, and center-of-mass motion. The grid is evolved at the same time order as the fluid variables. The user specifies grid evolution laws, which can be independent of the fluid motion. Applying our implementation to standard hydrodynamic test cases leads to improved results and higher efficiency, compared to the fixed-grid solutions. △ Less

Submitted 23 September, 2021; v1 submitted 8 September, 2021; originally announced September 2021.

Comments: 11 pages, 8 figures

Angular Momenta, Magnetization, and Accretion of Protostellar Cores

Authors: Aleksandra Kuznetsova, Lee Hartmann, Fabian Heitsch

Abstract: Building on our previous hydrodynamic study of the angular momenta of cloud cores formed during gravitational collapse of star-forming molecular gas in our previous work, we now examine core properties assuming ideal magnetohydrodynamics (MHD). Using the same sink-patch implementation for the \emph{Athena} MHD code, we characterize the statistical properties of cores, including the mass accretion… ▽ More Building on our previous hydrodynamic study of the angular momenta of cloud cores formed during gravitational collapse of star-forming molecular gas in our previous work, we now examine core properties assuming ideal magnetohydrodynamics (MHD). Using the same sink-patch implementation for the \emph{Athena} MHD code, we characterize the statistical properties of cores, including the mass accretion rates, specific angular momenta, and alignments between the magnetic field and the spin axis of the core on the $0.1 \ \mathrm{pc}$ scale. Our simulations, which reproduce the observed relation between magnetic field strength and gas density, show that magnetic fields can help collimate low density flows and help seed the locations of filamentary structures. Consistent with our previous purely hydrodynamic simulations, stars (sinks) form within the heterogeneous environments of filaments, such that accretion onto cores is highly episodic leading to short-term variability but no long-term monotonic growth of the specific angular momenta. With statistical characterization of protostellar cores properties and behaviors, we aim to provide a starting point for building more realistic and self-consistent disk formation models, helping to address whether magnetic fields can prevent the development of (large) circumstellar disks in the ideal MHD limit. △ Less

Submitted 30 March, 2020; originally announced March 2020.

Comments: 13 pages, 9 figures, accepted by ApJ

Gas Inflow and Star Formation near Supermassive Black Holes: The Role of Nuclear Activity

Authors: Christopher C. Frazer, Fabian Heitsch

Abstract: Numerical models of gas inflow towards a supermassive black hole (SMBH) show that star formation may occur in such an environment through the growth of a gravitationally unstable gas disc. We consider the effect of nuclear activity on such a scenario. We present the first three-dimensional grid-based radiative hydrodynamic simulations of direct collisions between infalling gas streams and a… ▽ More Numerical models of gas inflow towards a supermassive black hole (SMBH) show that star formation may occur in such an environment through the growth of a gravitationally unstable gas disc. We consider the effect of nuclear activity on such a scenario. We present the first three-dimensional grid-based radiative hydrodynamic simulations of direct collisions between infalling gas streams and a $4 \times 10^6~\text{M}_\odot$ SMBH, using ray-tracing to incorporate radiation consistent with an active galactic nucleus (AGN). We assume inflow masses of $ \approx 10^5~\text{M}_\odot$ and explore radiation fields of 10% and 100% of the Eddington luminosity ($L_\text{edd}$). We follow our models to the point of central gas disc formation preceding star formation and use the Toomre Q parameter ($Q_T$) to test for gravitational instability. We find that radiation pressure from UV photons inhibits inflow. Yet, for weak radiation fields, a central disc forms on timescales similar to that of models without feedback. Average densities of $> 10^{8}~\text{cm}^{-3}$ limit photo-heating to the disc surface allowing for $Q_T\approx1$. For strong radiation fields, the disc forms more gradually resulting in lower surface densities and larger $Q_T$ values. Mass accretion rates in our models are consistent with 1%--60% of the Eddington limit, thus we conclude that it is unlikely that radiative feedback from AGN activity would inhibit circumnuclear star formation arising from a massive inflow event. △ Less

Submitted 26 July, 2019; originally announced July 2019.

Comments: 29 pages, 27 figures, accepted by MNRAS

The Origins of Protostellar Core Angular Momenta

Authors: Aleksandra Kuznetsova, Lee Hartmann, Fabian Heitsch

Abstract: We present the results of a suite of numerical simulations designed to explore the origin of the angular momenta of protostellar cores. Using the hydrodynamic grid code \emph{Athena} with a sink implementation, we follow the formation of protostellar cores and protostars (sinks) from the subvirial collapse of molecular clouds on larger scales to investigate the range and relative distribution of c… ▽ More We present the results of a suite of numerical simulations designed to explore the origin of the angular momenta of protostellar cores. Using the hydrodynamic grid code \emph{Athena} with a sink implementation, we follow the formation of protostellar cores and protostars (sinks) from the subvirial collapse of molecular clouds on larger scales to investigate the range and relative distribution of core properties. We find that the core angular momenta are relatively unaffected by large-scale rotation of the parent cloud; instead, we infer that angular momenta are mainly imparted by torques between neighboring mass concentrations and exhibit a log-normal distribution. Our current simulation results are limited to size scales $\sim 0.05$~pc ($\sim 10^4 \rm AU$), but serve as first steps toward the ultimate goal of providing initial conditions for higher-resolution studies of core collapse to form protoplanetary disks. △ Less

Submitted 22 March, 2019; originally announced March 2019.

Comments: 12 pages, 9 figures, accepted to ApJ

The Role of Gravity in Producing Power-Law Mass Functions

Authors: Aleksandra Kuznetsova, Lee Hartmann, Fabian Heitsch, Javier Ballesteros-Paredes

Abstract: Numerical simulations of star formation have found that a power-law mass function can develop at high masses. In a previous paper, we employed isothermal simulations which created large numbers of sinks over a large range in masses to show that the power law exponent of the mass function, $dN/d\log M \propto M^Γ$, asymptotically and accurately approaches $Γ= -1.$ Simple analytic models show that s… ▽ More Numerical simulations of star formation have found that a power-law mass function can develop at high masses. In a previous paper, we employed isothermal simulations which created large numbers of sinks over a large range in masses to show that the power law exponent of the mass function, $dN/d\log M \propto M^Γ$, asymptotically and accurately approaches $Γ= -1.$ Simple analytic models show that such a power law can develop if the mass accretion rate $\dot{M} \propto M^2$, as in Bondi-Hoyle accretion; however, the sink mass accretion rates in the simulations show significant departures from this relation. In this paper we show that the expected accretion rate dependence is more closely realized provided the gravitating mass is taken to be the sum of the sink mass and the mass in the near environment. This reconciles the observed mass functions with the accretion rate dependencies, and demonstrates that power-law upper mass functions are essentially the result of gravitational focusing, a mechanism present in, for example, the competitive accretion model. △ Less

Submitted 8 October, 2018; originally announced October 2018.

Comments: 11 pages, 10 figures, accepted by ApJ

Exploring the Dust Content of Galactic Winds With Herschel. II. Nearby Dwarf Galaxies

Authors: Alexander McCormick, Sylvain Veilleux, Marcio Meléndez, Crystal L. Martin, Joss Bland-Hawthorn, Gerald Cecil, Fabian Heitsch, Thomas Müller, David S. N. Rupke, Chad Engelbracht

Abstract: We present the results from an analysis of deep Herschel Space Observatory observations of six nearby dwarf galaxies known to host galactic-scale winds. The superior far-infrared sensitivity and angular resolution of Herschel have allowed detection of cold circumgalactic dust features beyond the stellar components of the host galaxies traced by Spitzer 4.5 $μ$m images. Comparisons of these cold du… ▽ More We present the results from an analysis of deep Herschel Space Observatory observations of six nearby dwarf galaxies known to host galactic-scale winds. The superior far-infrared sensitivity and angular resolution of Herschel have allowed detection of cold circumgalactic dust features beyond the stellar components of the host galaxies traced by Spitzer 4.5 $μ$m images. Comparisons of these cold dust features with ancillary data reveal an imperfect spatial correlation with the ionized gas and warm dust wind components. We find that typically $\sim$10-20\% of the total dust mass in these galaxies resides outside of their stellar disks, but this fraction reaches $\sim$60\% in the case of NGC 1569. This galaxy also has the largest metallicity (O/H) deficit in our sample for its stellar mass. Overall, the small number of objects in our sample precludes drawing strong conclusions on the origin of the circumgalactic dust. We detect no statistically significant trends with star formation properties of the host galaxies, as might be expected if the dust were lifted above the disk by energy inputs from on-going star formation activity. Although a case for dust entrained in a galactic wind is seen in NGC 1569, in all cases, we cannot rule out the possibility that some of the circumgalactic dust might be associated instead with gas accreted or removed from the disk by recent galaxy interaction events, or that it is part of the outer gas-rich portion of the disk that lies below the sensitivity limit of the Spitzer 4.5 $μ$m data. △ Less

Submitted 4 April, 2018; v1 submitted 8 March, 2018; originally announced March 2018.

Comments: 39 pages, 10 figures, Accepted for publication in MNRAS

A systematic comparison of two-equation RANS turbulence models applied to shock-cloud interactions

Authors: Matthew D. Goodson, Fabian Heitsch, Karl Eklund, Virginia A. Williams

Abstract: Turbulence models attempt to account for unresolved dynamics and diffusion in hydrodynamical simulations. We develop a common framework for two-equation Reynolds-Averaged Navier-Stokes (RANS) turbulence models, and we implement six models in the Athena code. We verify each implementation with the standard subsonic mixing layer, although the level of agreement depends on the definition of the mixin… ▽ More Turbulence models attempt to account for unresolved dynamics and diffusion in hydrodynamical simulations. We develop a common framework for two-equation Reynolds-Averaged Navier-Stokes (RANS) turbulence models, and we implement six models in the Athena code. We verify each implementation with the standard subsonic mixing layer, although the level of agreement depends on the definition of the mixing layer width. We then test the validity of each model into the supersonic regime, showing that compressibility corrections can improve agreement with experiment. For models with buoyancy effects, we also verify our implementation via the growth of the Rayleigh-Taylor instability in a stratified medium. The models are then applied to the ubiquitous astrophysical shock-cloud interaction in three dimensions. We focus on the mixing of shock and cloud material, comparing results from turbulence models to high-resolution simulations (up to 200 cells per cloud radius) and ensemble-averaged simulations. We find that the turbulence models lead to increased spreading and mixing of the cloud, although no two models predict the same result. Increased mixing is also observed in inviscid simulations at resolutions greater than 100 cells per radius, which suggests that the turbulent mixing begins to be resolved. △ Less

Submitted 25 March, 2017; originally announced March 2017.

Comments: 20 pages, 17 figures, accepted to MNRAS

The impact of magnetic fields on the chemical evolution of the supernova-driven ISM

Authors: Anabele-Linda Pardi, Philipp Girichidis, Thorsten Naab, Stefanie Walch, Thomas Peters, Fabian Heitsch, Simon C. O. Glover, Ralf S. Klessen, Richard Wünsch, Andrea Gatto

Abstract: We present three-dimensional magneto-hydrodynamical simulations of the self-gravitating interstellar medium (ISM) in a periodic (256 pc)$^3$ box with a mean number density of 0.5 cm$^{-3}$. At a fixed supernova rate we investigate the multi-phase ISM structure, H$_{2}$ molecule formation and density-magnetic field scaling for varying initial magnetic field strengths (0, $6\times 10^{-3}$, 0.3, 3… ▽ More We present three-dimensional magneto-hydrodynamical simulations of the self-gravitating interstellar medium (ISM) in a periodic (256 pc)$^3$ box with a mean number density of 0.5 cm$^{-3}$. At a fixed supernova rate we investigate the multi-phase ISM structure, H$_{2}$ molecule formation and density-magnetic field scaling for varying initial magnetic field strengths (0, $6\times 10^{-3}$, 0.3, 3 $μ$G). All magnetic runs saturate at mass weighted field strengths of $\sim$ 1 $-$ 3 $μ$G but the ISM structure is notably different. With increasing initial field strengths (from $6\times 10^{-3}$ to 3 $μ$G) the simulations develop an ISM with a more homogeneous density and temperature structure, with increasing mass (from 5% to 85%) and volume filling fractions (from 4% to 85%) of warm (300 K $<$ T $<$ 8000 K) gas, with decreasing volume filling fractions (VFF) from $\sim$ 35% to $\sim$ 12% of hot gas (T $> 10^5$ K) and with a decreasing H$_{2}$ mass fraction (from 70% to $<$ 1%). Meanwhile the mass fraction of gas in which the magnetic pressure dominates over the thermal pressure increases by a factor of 10, from 0.07 for an initial field of $6\times 10^{-3}$ $μ$G to 0.7 for a 3 $μ$G initial field. In all but the simulations with the highest initial field strength self-gravity promotes the formation of dense gas and H$_{2}$, but does not change any other trends. We conclude that magnetic fields have a significant impact on the multi-phase, chemical and thermal structure of the ISM and discuss potential implications and limitations of the model. △ Less

Submitted 2 November, 2016; originally announced November 2016.

Comments: submitted

Structure, Dynamics and Deuterium Fractionation of Massive Pre-Stellar Cores

Authors: Matthew D. Goodson, Shuo Kong, Jonathan C. Tan, Fabian Heitsch, Paola Caselli

Abstract: High levels of deuterium fraction in N$_2$H$^+$ are observed in some pre-stellar cores. Single-zone chemical models find that the timescale required to reach observed values ($D_{\rm frac}^{{\rm N}_2{\rm H}^+} \equiv {\rm N}_2{\rm D}^+/{\rm N}_2{\rm H}^+ \gtrsim 0.1$) is longer than the free-fall time, possibly ten times longer. Here, we explore the deuteration of turbulent, magnetized cores with… ▽ More High levels of deuterium fraction in N$_2$H$^+$ are observed in some pre-stellar cores. Single-zone chemical models find that the timescale required to reach observed values ($D_{\rm frac}^{{\rm N}_2{\rm H}^+} \equiv {\rm N}_2{\rm D}^+/{\rm N}_2{\rm H}^+ \gtrsim 0.1$) is longer than the free-fall time, possibly ten times longer. Here, we explore the deuteration of turbulent, magnetized cores with 3D magnetohydrodynamics simulations. We use an approximate chemical model to follow the growth in abundances of N$_2$H$^+$ and N$_2$D$^+$. We then examine the dynamics of the core using each tracer for comparison to observations. We find that the velocity dispersion of the core as traced by N$_2$D$^+$ appears slightly sub-virial compared to predictions of the Turbulent Core Model of McKee & Tan, except at late times just before the onset of protostar formation. By varying the initial mass surface density, the magnetic energy, the chemical age, and the ortho-to-para ratio of H$_2$, we also determine the physical and temporal properties required for high deuteration. We find that low initial ortho-to-para ratios ($\lesssim 0.01$) and/or multiple free-fall times ($\gtrsim 3$) of prior chemical evolution are necessary to reach the observed values of deuterium fraction in pre-stellar cores. △ Less

Submitted 30 October, 2016; v1 submitted 22 September, 2016; originally announced September 2016.

Comments: 20 pages, 18 figures; accepted for publication in ApJ

Reorienting MHD Colliding Flows: A Shock Physics Mechanism for Generating Filaments Normal to Magnetic Fields

Authors: Erica L. Fogerty, Jonathan Carroll-Nellenback, Adam Frank, Fabian Heitsch, Andy Pon

Abstract: We present numerical simulations of reorienting oblique shocks that form in the collision layer between magnetized colliding flows. Reorientation aligns parsec-scale post-shock filaments normal to the background magnetic field. We find that reorientation begins with pressure gradients between the collision region and the ambient medium. This drives a lateral expansion of post-shock gas, which reor… ▽ More We present numerical simulations of reorienting oblique shocks that form in the collision layer between magnetized colliding flows. Reorientation aligns parsec-scale post-shock filaments normal to the background magnetic field. We find that reorientation begins with pressure gradients between the collision region and the ambient medium. This drives a lateral expansion of post-shock gas, which reorients the growing filament from the outside-in (i.e. from the flow-ambient boundary, toward the colliding flows axis). The final structures of our simulations resemble polarization observations of filaments in Taurus and Serpens South, as well as the integral-shaped filament in Orion A. Given the ubiquity of colliding flows in the interstellar medium, shock reorientation may be relevant to the formation of filaments normal to magnetic fields. △ Less

Submitted 7 June, 2017; v1 submitted 9 September, 2016; originally announced September 2016.

Comments: In press, MNRAS, 11 pages, 13 figures

Scaling Relations of Starburst-Driven Galactic Winds

Authors: Ryan Tanner, Gerald Cecil, Fabian Heitsch

Abstract: Using synthetic absorption lines generated from 3D hydro-dynamical simulations we explore how the velocity of a starburst-driven galactic wind correlates with the star formation rate (SFR) and SFR density. We find strong correlations until the scaling relations flatten abruptly at a point set by the mass loading of the starburst. Below this point the scaling relation depends on the temperature reg… ▽ More Using synthetic absorption lines generated from 3D hydro-dynamical simulations we explore how the velocity of a starburst-driven galactic wind correlates with the star formation rate (SFR) and SFR density. We find strong correlations until the scaling relations flatten abruptly at a point set by the mass loading of the starburst. Below this point the scaling relation depends on the temperature regime being probed by the absorption line, not on the mass loading. The exact scaling relation depends on whether the maximum or mean velocity of the absorption line is used. We find that the outflow velocity of neutral gas is four to five times lower than the average velocity of the hottest gas, with the difference in velocity between the neutral and ionized gas increasing with gas ionization. Thus, absorption lines of neutral or low ionized gas will underestimate the outflow velocity of hot gas, severely underestimating outflow energetics. △ Less

Submitted 18 August, 2016; originally announced August 2016.

Comments: 12 pages, 8 figures, Submitted to ApJ

Chemical enrichment of the pre-solar cloud by supernova dust grains

Authors: Matthew D. Goodson, Ian Luebbers, Fabian Heitsch, Christopher C. Frazer

Abstract: The presence of short-lived radioisotopes (SLRs) in solar system meteorites has been interpreted as evidence that the solar system was exposed to a supernova shortly before or during its formation. Yet results from hydrodynamical models of SLR injection into the proto-solar cloud or disc suggest that gas-phase mixing may not be efficient enough to reproduce the observed abundances. As an alternati… ▽ More The presence of short-lived radioisotopes (SLRs) in solar system meteorites has been interpreted as evidence that the solar system was exposed to a supernova shortly before or during its formation. Yet results from hydrodynamical models of SLR injection into the proto-solar cloud or disc suggest that gas-phase mixing may not be efficient enough to reproduce the observed abundances. As an alternative, we explore the injection of SLRs via dust grains as a way to overcome the mixing barrier. We numerically model the interaction of a supernova remnant containing SLR-rich dust grains with a nearby molecular cloud. The dust grains are subject to drag forces and both thermal and non-thermal sputtering. We confirm that the expanding gas shell stalls upon impact with the dense cloud and that gas-phase SLR injection occurs slowly due to hydrodynamical instabilities at the cloud surface. In contrast, dust grains of sufficient size (> 1 micron) decouple from the gas and penetrate into the cloud within 0.1 Myr. Once inside the cloud, the dust grains are destroyed by sputtering, releasing SLRs and rapidly enriching the dense (potentially star-forming) regions. Our results suggest that SLR transport on dust grains is a viable mechanism to explain SLR enrichment. △ Less

Submitted 26 August, 2016; v1 submitted 22 July, 2016; originally announced July 2016.

Comments: 15 pages, 10 figures, Accepted for publication in MNRAS. Movies can be found here: http://user.physics.unc.edu/~mdgood86/research.html

Three-Dimensional Orientation of Compact High Velocity Clouds

Authors: F. Heitsch, B. Bartell, S. E. Clark, J. E. G. Peek, D. Cheng, M. Putman

Abstract: We present a proof-of-concept study of a method to estimate the inclination angle of compact high velocity clouds (CHVCs), i.e. the angle between a CHVC's trajectory and the line-of-sight. The inclination angle is derived from the CHVC's morphology and kinematics. We calibrate the method with numerical simulations, and we apply it to a sample of CHVCs drawn from HIPASS. Implications for CHVC dista… ▽ More We present a proof-of-concept study of a method to estimate the inclination angle of compact high velocity clouds (CHVCs), i.e. the angle between a CHVC's trajectory and the line-of-sight. The inclination angle is derived from the CHVC's morphology and kinematics. We calibrate the method with numerical simulations, and we apply it to a sample of CHVCs drawn from HIPASS. Implications for CHVC distances are discussed. △ Less

Submitted 21 June, 2016; originally announced June 2016.

Comments: 5 pages, 4 figures, accepted by MNRAS Letters

Molecular cloud formation in high-shear, magnetized colliding flows

Authors: Erica Fogerty, Adam Frank, Fabian Heitsch, Jonathan Carroll-Nellenback, Christina Haig, Marissa Adams

Abstract: The colliding flows (CF) model is a well-supported mechanism for generating molecular clouds. However, to-date most CF simulations have focused on the formation of clouds in the normal-shock layer between head-on colliding flows. We performed simulations of magnetized colliding flows that instead meet at an oblique-shock layer. Oblique shocks generate shear in the post-shock environment, and this… ▽ More The colliding flows (CF) model is a well-supported mechanism for generating molecular clouds. However, to-date most CF simulations have focused on the formation of clouds in the normal-shock layer between head-on colliding flows. We performed simulations of magnetized colliding flows that instead meet at an oblique-shock layer. Oblique shocks generate shear in the post-shock environment, and this shear creates inhospitable environments for star formation. As the degree of shear increases (i.e. the obliquity of the shock increases), we find that it takes longer for sink particles to form, they form in lower numbers, and they tend to be less massive. With regard to magnetic fields, we find that even a weak field stalls gravitational collapse within forming clouds. Additionally, an initially oblique collision interface tends to reorient over time in the presence of a magnetic field, so that it becomes normal to the oncoming flows. This was demonstrated by our most oblique shock interface, which became fully normal by the end of the simulation. △ Less

Submitted 30 May, 2016; v1 submitted 1 February, 2016; originally announced February 2016.

Comments: MNRAS, in press, Monthly Notices of the Royal Astronomical Society 2016

On the Metallicity and Origin of the Smith High-Velocity Cloud

Authors: Andrew J. Fox, Nicolas Lehner, Felix J. Lockman, Bart P. Wakker, Alex S. Hill, Fabian Heitsch, David V. Stark, Kathleen A. Barger, Kenneth R. Sembach, Mubdi Rahman

Abstract: The Smith Cloud is a gaseous high-velocity cloud (HVC) in an advanced state of accretion, only 2.9 kpc below the Galactic plane and due to impact the disk in 27 Myr. It is unique among HVCs in having a known distance (12.4+/-1.3 kpc) and a well-constrained 3D velocity (296 km/s), but its origin has long remained a mystery. Here we present the first absorption-line measurements of its metallicity,… ▽ More The Smith Cloud is a gaseous high-velocity cloud (HVC) in an advanced state of accretion, only 2.9 kpc below the Galactic plane and due to impact the disk in 27 Myr. It is unique among HVCs in having a known distance (12.4+/-1.3 kpc) and a well-constrained 3D velocity (296 km/s), but its origin has long remained a mystery. Here we present the first absorption-line measurements of its metallicity, using HST/COS UV spectra of three AGN lying behind the Cloud together with Green Bank Telescope 21 cm spectra of the same directions. Using Voigt-profile fitting of the S II 1250, 1253, 1259 triplet together with ionization corrections derived from photoionization modeling, we derive the sulfur abundance in each direction; a weighted average of the three measurements gives [S/H]=-0.28+/-0.14, or 0.53+0.21-0.15 solar metallicity. The finding that the Smith Cloud is metal-enriched lends support to scenarios where it represents recycled Galactic material rather than the remnant of a dwarf galaxy or accreting intergalactic gas. The metallicity and trajectory of the Cloud are both indicative of an origin in the outer disk. However, its large mass and prograde kinematics remain to be fully explained. If the cloud has accreted cooling gas from the corona during its fountain trajectory, as predicted in recent theoretical work, its current mass would be higher than its launch mass, alleviating the mass concern. △ Less

Submitted 15 December, 2015; originally announced December 2015.

Comments: Accepted for publication in ApJL

Starburst Driven Galactic Superbubbles Radiating to 10 K

Authors: Ryan Tanner, Gerald Cecil, Fabian Heitsch

Abstract: Our three-dimensional hydro-dynamical simulations of starbursts examine the formation of superbubbles over a range of driving luminosities and mass loadings that determine superbubble growth and wind velocity. From this we determine the relationship between the velocity of a galactic wind and the power of the starburst. We find a threshold for the formation of a wind, above which the speed of the… ▽ More Our three-dimensional hydro-dynamical simulations of starbursts examine the formation of superbubbles over a range of driving luminosities and mass loadings that determine superbubble growth and wind velocity. From this we determine the relationship between the velocity of a galactic wind and the power of the starburst. We find a threshold for the formation of a wind, above which the speed of the wind is not affected by grid resolution or the temperature floor of our radiative cooling. We investigate the effect two different temperature floors in our radiative cooling prescription have on wind kinematics and content. We find that cooling to $10$ K instead of to $10^4$ K increases the mass fraction of cold neutral and hot X-ray gas in the galactic wind while halving that in warm H$α$. Our simulations show the mass of cold gas transported into the lower halo does not depend on the starburst strength. Optically bright filaments form at the edge of merging superbubbles, or where a cold dense cloud has been disrupted by the wind. Filaments formed by merging superbubbles will persist and grow to $>400$ pc in length if anchored to a star forming complex. Filaments embedded in the hot galactic wind contain warm and cold gas that moves $300-1200$ km s$^{-1}$ slower than the surrounding wind, with the coldest gas hardly moving with respect to the galaxy. Warm and cold matter in the galactic wind show asymmetric absorption profiles consistent with observations, with a thin tail up to the wind velocity. △ Less

Submitted 18 August, 2016; v1 submitted 16 November, 2015; originally announced November 2015.

Comments: 14 pages, 14 figures. Published in ApJ

Journal ref: ApJ, 821, 7 2016

THOR - The HI, OH, Recombination Line Survey of the Milky Way - The pilot study: HI observations of the giant molecular cloud W43

Authors: S. Bihr, H. Beuther, J. Ott, K. G. Johnston, A. Brunthaler, L. D. Anderson, F. Bigiel, P. Carlhoff, E. Churchwell, S. C. O. Glover, P. F. Goldsmith, F. Heitsch, T. Henning, M. H. Heyer, T. Hill, A. Hughes, R. S. Klessen, H. Linz, S. N. Longmore, N. M. McClure-Griffiths, K. M. Menten, F. Motte, Q. Nguyen-Lu'o'ng, R. Plume, S. E. Ragan , et al. (8 additional authors not shown)

Abstract: To study the atomic, molecular and ionized emission of Giant Molecular Clouds (GMCs), we have initiated a Large Program with the VLA: 'THOR - The HI, OH, Recombination Line survey of the Milky Way'. We map the 21cm HI line, 4 OH lines, 19 H_alpha recombination lines and the continuum from 1 to 2 GHz of a significant fraction of the Milky Way (l=15-67deg, |b|<1deg) at ~20" resolution. In this paper… ▽ More To study the atomic, molecular and ionized emission of Giant Molecular Clouds (GMCs), we have initiated a Large Program with the VLA: 'THOR - The HI, OH, Recombination Line survey of the Milky Way'. We map the 21cm HI line, 4 OH lines, 19 H_alpha recombination lines and the continuum from 1 to 2 GHz of a significant fraction of the Milky Way (l=15-67deg, |b|<1deg) at ~20" resolution. In this paper, we focus on the HI emission from the W43 star-formation complex. Classically, the HI 21cm line is treated as optically thin with column densities calculated under this assumption. This might give reasonable results for regions of low-mass star-formation, however, it is not sufficient to describe GMCs. We analyzed strong continuum sources to measure the optical depth, and thus correct the HI 21cm emission for optical depth effects and weak diffuse continuum emission. Hence, we are able to measure the HI mass of W43 more accurately and our analysis reveals a lower limit of M=6.6x10^6 M_sun, which is a factor of 2.4 larger than the mass estimated with the assumption of optically thin emission. The HI column densities are as high as N(HI)~150 M_sun/pc^2 ~ 1.9x10^22 cm^-2, which is an order of magnitude higher than for low mass star formation regions. This result challenges theoretical models that predict a threshold for the HI column density of ~10 M_sun/pc^2, at which the formation of molecular hydrogen should set in. By assuming an elliptical layered structure for W43, we estimate the particle density profiles. While at the cloud edge atomic and molecular hydrogen are well mixed, the center of the cloud is strongly dominated by molecular hydrogen. We do not identify a sharp transition between hydrogen in atomic and molecular form. Our results are an important characterization of the atomic to molecular hydrogen transition in an extreme environment and challenges current theoretical models. △ Less

Submitted 19 May, 2015; originally announced May 2015.

Journal ref: 2015, A&A, 580, A112

Exploring the Dust Content of Galactic Winds with Herschel. I. NGC 4631

Authors: M. Meléndez, S. Veilleux, C. Martin, C. Engelbracht, J. Bland-Hawthorn, G. Cecil, F. Heitsch, A. McCormick, T. Müller, D. Rupke, S. H. Teng

Abstract: We present a detailed analysis of deep far-infrared observations of the nearby edge-on star-forming galaxy NGC 4631 obtained with the Herschel Space Observatory. Our PACS images at 70 and 160 um show a rich complex of filaments and chimney-like features that extends up to a projected distance of 6 kpc above the plane of the galaxy. The PACS features often match extraplanar Halpha, radio-continuum,… ▽ More We present a detailed analysis of deep far-infrared observations of the nearby edge-on star-forming galaxy NGC 4631 obtained with the Herschel Space Observatory. Our PACS images at 70 and 160 um show a rich complex of filaments and chimney-like features that extends up to a projected distance of 6 kpc above the plane of the galaxy. The PACS features often match extraplanar Halpha, radio-continuum, and soft X-ray features observed in this galaxy, pointing to a tight disk-halo connection regulated by star formation. On the other hand, the morphology of the colder dust component detected on larger scale in the SPIRE 250, 350, and 500 um data matches the extraplanar H~I streams previously reported in NGC 4631 and suggests a tidal origin. The PACS 70/160 ratios are elevated in the central ~3.0 kpc region above the nucleus of this galaxy (the "superbubble"). A pixel-by-pixel analysis shows that dust in this region has a higher temperature and/or an emissivity with a steeper spectral index (beta > 2) than the dust in the disk, possibly the result of the harsher environment in the superbubble. Star formation in the disk seems energetically insufficient to lift the material out of the disk, unless it was more active in the past or the dust-to-gas ratio in the superbubble region is higher than the Galactic value. Some of the dust in the halo may also have been tidally stripped from nearby companions or lifted from the disk by galaxy interactions. △ Less

Submitted 26 February, 2015; originally announced February 2015.

Comments: Accepted for publication in The Astrophysical Journal

Accretion and Diffusion Timescales in Sheets and Filaments

Authors: Fabian Heitsch, Lee Hartmann

Abstract: A comparison of accretion and (turbulent) magnetic diffusion timescales for sheets and filaments demonstrates that dense star-forming clouds generally will -- under realistic conditions -- become supercritical due to mass accretion on timescales at least an order of magnitude shorter than ambipolar and/or turbulent diffusion timescales. Thus, ambipolar or turbulent diffusion -- while present -- is… ▽ More A comparison of accretion and (turbulent) magnetic diffusion timescales for sheets and filaments demonstrates that dense star-forming clouds generally will -- under realistic conditions -- become supercritical due to mass accretion on timescales at least an order of magnitude shorter than ambipolar and/or turbulent diffusion timescales. Thus, ambipolar or turbulent diffusion -- while present -- is unlikely to control the formation of cores and stars. △ Less

Submitted 9 June, 2014; originally announced June 2014.

Comments: 12 pages, 6 figures, accepted by MNRAS

The formation of the W43 complex: constraining its atomic-to-molecular transition and searching for colliding clouds

Authors: F. Motte, Q. Nguyen Luong, N. Schneider, F. Heitsch, S. Glover, P. Carlhoff, T. Hill, S. Bontemps, P. Schilke, F. Louvet, M. Hennemann, P. Didelon, H. Beuther

Abstract: Numerical simulations have explored the possibility to form molecular clouds through either a quasi-static, self-gravitating mechanism or the collision of gas streams or lower-density clouds. They also quantitatively predict the distribution of matter at the transition from atomic to molecular gases. We aim to observationally test these models by studying the environment of W43, a molecular cloud… ▽ More Numerical simulations have explored the possibility to form molecular clouds through either a quasi-static, self-gravitating mechanism or the collision of gas streams or lower-density clouds. They also quantitatively predict the distribution of matter at the transition from atomic to molecular gases. We aim to observationally test these models by studying the environment of W43, a molecular cloud complex near the tip of the Galactic long bar. Using Galaxy-wide HI and 12CO surveys we searched for gas flowing toward the W43 molecular cloud complex. We also estimated the HI and H2 mass surface densities to constrain the transition from atomic to molecular gas around and within W43. We found 3 cloud ensembles within the position-velocity diagrams of 12CO and HI gases. They are separated by 20km/s along the line of sight and extend into the 13CO velocity structure of W43. Since their velocity gradients are consistent with free-fall, they could be nearby clouds attracted by, and streaming toward, the W43 10^7Msun potential well. We show that the HI surface density, Sigma_HI=45-85Msun/pc2, does not reach any threshold level but increases when entering the 130pc-wide molecular complex previously defined. This suggests that an equilibrium between H2 formation and photodissociation has not yet been reached. The H2-to-HI ratio measured over the W43 region and its surroundings, R_H2~3.5, is high, indicating that most of the gas is already in molecular form in W43 and in structures several hundreds of parsecs downstream along the Scutum-Centaurus arm. The W43 molecular cloud complex may have formed, and in fact may still be accreting mass from the agglomeration of clouds. Already in the molecular-dominated regime, most of these clouds are streaming from the Scutum-Centaurus arm. This is in clear disagreement with quasi-static and steady-state models of molecular cloud formation. △ Less

Submitted 16 April, 2014; originally announced April 2014.

Comments: 12 pages, 11 figures, accepted by A&A

Journal ref: A&A 571, A32 (2014)

Large scale IRAM 30m CO-observations in the giant molecular cloud complex W43

Authors: P. Carlhoff, Q. Nguyen Luong, P. Schilke, F. Motte, N. Schneider, H. Beuther, S. Bontemps, F. Heitsch, T. Hill, C. Kramer, V. Ossenkopf, F. Schuller, R. Simon, F. Wyrowski

Abstract: We aim to give a full description of the distribution and location of dense molecular clouds in the giant molecular cloud complex W43. It has previously been identified as one of the most massive star-forming regions in our Galaxy. To trace the moderately dense molecular clouds in the W43 region, we initiated an IRAM 30m large program, named W43-HERO, covering a large dynamic range of scales (from… ▽ More We aim to give a full description of the distribution and location of dense molecular clouds in the giant molecular cloud complex W43. It has previously been identified as one of the most massive star-forming regions in our Galaxy. To trace the moderately dense molecular clouds in the W43 region, we initiated an IRAM 30m large program, named W43-HERO, covering a large dynamic range of scales (from 0.3 to 140 pc). We obtained on-the-fly-maps in 13CO (2-1) and C18O (2-1) with a high spectral resolution of 0.1 km/s and a spatial resolution of 12". These maps cover an area of ~1.5 square degrees and include the two main clouds of W43, as well as the lower density gas surrounding them. A comparison with Galactic models and previous distance calculations confirms the location of W43 near the tangential point of the Scutum arm at a distance from the Sun of approximately 6 kpc. The resulting intensity cubes of the observed region are separated into sub-cubes, centered on single clouds which are then analyzed in detail. The optical depth, excitation temperature, and H2 column density maps are derived out of the 13CO and C18O data. These results are then compared with those derived from Herschel dust maps. The mass of a typical cloud is several 10^4 solar masses while the total mass in the dense molecular gas (>100 cm^-3) in W43 is found to be about 1.9e6 solar masses. Probability distribution functions obtained from column density maps derived from molecular line data and Herschel imaging show a log-normal distribution for low column densities and a power-law tail for high densities. A flatter slope for the molecular line data PDF may imply that those selectively show the gravitationally collapsing gas. △ Less

Submitted 16 September, 2013; v1 submitted 28 August, 2013; originally announced August 2013.

Gravitational Infall onto Molecular Filaments II. Externally Pressurized Cylinders

Authors: Fabian Heitsch

Abstract: In an extension of Fischera & Martin (2012a) and Heitsch (2013), two aspects of the evolution of externally pressurized, hydrostatic filaments are discussed. (a) The free-fall accretion of gas onto such a filament will lead to filament parameters (specifically, FWHM--column density relations) inconsistent with the observations of Arzoumanian et al. (2011), except for two cases: For low-mass, isoth… ▽ More In an extension of Fischera & Martin (2012a) and Heitsch (2013), two aspects of the evolution of externally pressurized, hydrostatic filaments are discussed. (a) The free-fall accretion of gas onto such a filament will lead to filament parameters (specifically, FWHM--column density relations) inconsistent with the observations of Arzoumanian et al. (2011), except for two cases: For low-mass, isothermal filaments, agreement is found as in the analysis by Fischera & Martin (2012b). Magnetized cases, for which the field scales weakly with the density as $B\propto n^{1/2}$, also reproduce observed parameters. (b) Realistically, the filaments will be embedded not only in gas of non-zero pressure, but also of non-zero density. Thus, the appearance of sheet-embedded filaments is explored. Generating a grid of filament models and comparing the resulting column density ratios and profile shapes with observations suggests that the three-dimensional filament profiles are intrinsically flatter than isothermal, beyond projection and evolution effects. △ Less

Submitted 26 August, 2013; originally announced August 2013.

Comments: 15 pages, 7 figures, accepted by ApJ

Low-velocity shocks traced by extended SiO emission along the W43 ridges: witnessing the formation of young massive clusters

Authors: Q. Nguyen Luong, F. Motte, P. Carlhoff, F. Louvet, P. Lesaffre, P. Schilke, T. Hill, M. Hennemann, A. Gusdorf, N. Schneider, S. Bontemps, A. Duarte-Cabral, K. M. Menten, P. G. Martin, F. Wyrowski, G. Bendo, H. Roussel, J-P. Bernard, L. Bronfman, T. Henning, C. Kramer, F. Heitsch

Abstract: The formation of high-mass stars is tightly linked to that of their parental clouds. We here focus on the high-density parts of W43, a molecular cloud undergoing an efficient event of formation. The cloud structure is studied with a column density image derived from Herschel continuum maps obtained at 70, 160, 250, 350, and 500 micron. We identify two high-column density filamentary clouds, quoted… ▽ More The formation of high-mass stars is tightly linked to that of their parental clouds. We here focus on the high-density parts of W43, a molecular cloud undergoing an efficient event of formation. The cloud structure is studied with a column density image derived from Herschel continuum maps obtained at 70, 160, 250, 350, and 500 micron. We identify two high-column density filamentary clouds, quoted as the W43-MM1 and W43-MM2 ridges, which both account for 1.5x10^4 Msun gas mass above 10^23 cm-2 and within areas of 5 and 14pc^2, respectively. We used the N_2H^+ 1--0 line to confirm that the W43-MM1 and W43-MM2 ridges are structures coherent in velocity and gravitationally bound, despite their large velocity dispersion and ~5 kms line widths. The most intriguing result of the W43 large program is the bright wide-spread SiO 2--1 emission: 1--11 K kms$ stretching an area of ~28 pc^2. Concentrated toward the W43-MM1 and W43-MM2 ridges and their immediate surroundings, it leads to a total luminosity of L_SiO 2-1 ~4 10^4 K kms kpc^2pc^2. We measured a steep relation between the luminosity and velocity extent of the SiO~2--1 lines and propose to use it to distinguish the low-velocity shocks observed here from the more classical high-velocity ones associated with outflows of high-mass young stellar objects. We used state-of-the-art shock models to demonstrate that low-velocity (<10 kms^-1) shocks with a small amount (10%) of Si atoms initially in gas phase or in grain mantles can explain the observed SiO column density in W43. The spatial and velocity overlaps between the ridges high-density gas (n_H2>10^4-10^5 cm^-3) and the shocked SiO gas suggests that ridges could be forming via colliding flows driven by gravity and accompanied by low-velocity shocks. This mechanism may be the initial conditions for the formation of young massive clusters in these ridges. △ Less

Submitted 17 July, 2013; v1 submitted 3 June, 2013; originally announced June 2013.

Comments: accepted for publication to ApJ in July 2013

Gravitational Infall onto Molecular Filaments

Authors: Fabian Heitsch

Abstract: Two aspects of filamentary molecular cloud evolution are addressed: (1) Exploring analytically the role of the environment for the evolution of filaments demonstrates that considering them in isolation (i.e. just addressing the fragmentation stability) will result in unphysical conclusions about the filament's properties. Accretion can also explain the observed decorrelation between FWHM and peak… ▽ More Two aspects of filamentary molecular cloud evolution are addressed: (1) Exploring analytically the role of the environment for the evolution of filaments demonstrates that considering them in isolation (i.e. just addressing the fragmentation stability) will result in unphysical conclusions about the filament's properties. Accretion can also explain the observed decorrelation between FWHM and peak column density. (2) Free-fall accretion onto finite filaments can lead to the characteristic "fans" of infrared-dark clouds around star-forming regions. The fans may form due to tidal forces mostly arising at the ends of the filaments, consistent with numerical models and earlier analytical studies. △ Less

Submitted 16 April, 2013; originally announced April 2013.

Comments: 14 pages, 10 figures, accepted by ApJ

The Effects of Inhomogeneities within Colliding Flows on the Formation and Evolution of Molecular Clouds

Authors: Jonathan Carroll-Nellenback, Adam Frank, Fabian Heitsch

Abstract: Observational evidence from local star-forming regions mandates that star formation occurs shortly after, or even during, molecular cloud formation. Models of the formation of molecular clouds in large-scale converging flows have identified the physical mechanisms driving the necessary rapid fragmentation. They also point to global gravitational collapse driving supersonic turbulence in molecular… ▽ More Observational evidence from local star-forming regions mandates that star formation occurs shortly after, or even during, molecular cloud formation. Models of the formation of molecular clouds in large-scale converging flows have identified the physical mechanisms driving the necessary rapid fragmentation. They also point to global gravitational collapse driving supersonic turbulence in molecular clouds. Previous cloud formation models have focused on turbulence generation, gravitational collapse, magnetic fields, and feedback. Here, we explore the effect of structure in the flow on the resulting clouds and the ensuing gravitational collapse. We compare two extreme cases, one with a collision between two smooth streams, and one with streams containing small clumps. We find that structured converging flows lead to a delay of local gravitational collapse ("star formation"). Thus, more gas has time to accumulate, eventually leading to a strong global collapse, and thus to a high star formation rate. Uniform converging flows fragment hydrodynamically early on, leading to the rapid onset of local gravitational collapse and an overall low sink formation rate. △ Less

Submitted 9 April, 2013; v1 submitted 4 April, 2013; originally announced April 2013.

Comments: Submitted to ApJ

Aspect Ratio Dependence of the Free-Fall Time for Non-Spherical Symmetries

Authors: Andy Pon, Jesús A. Toalá, Doug Johnstone, Enrique Vázquez-Semadeni, Fabian Heitsch, Gilberto C. Gómez

Abstract: We investigate the collapse of non-spherical substructures, such as sheets and filaments, which are ubiquitous in molecular clouds. Such non-spherical substructures collapse homologously in their interiors but are influenced by an edge effect that causes their edges to be preferentially accelerated. We analytically compute the homologous collapse timescales of the interiors of uniform-density, sel… ▽ More We investigate the collapse of non-spherical substructures, such as sheets and filaments, which are ubiquitous in molecular clouds. Such non-spherical substructures collapse homologously in their interiors but are influenced by an edge effect that causes their edges to be preferentially accelerated. We analytically compute the homologous collapse timescales of the interiors of uniform-density, self-gravitating filaments and find that the homologous collapse timescale scales linearly with the aspect ratio. The characteristic timescale for an edge driven collapse mode in a filament, however, is shown to have a square root dependence on the aspect ratio. For both filaments and circular sheets, we find that selective edge acceleration becomes more important with increasing aspect ratio. In general, we find that lower dimensional objects and objects with larger aspect ratios have longer collapse timescales. We show that estimates for star formation rates, based upon gas densities, can be overestimated by an order of magnitude if the geometry of a cloud is not taken into account. △ Less

Submitted 27 July, 2012; v1 submitted 12 July, 2012; originally announced July 2012.

Comments: 10 pages, 2 figures, accepted by ApJ, minor grammatical errors fixed

Journal ref: ApJ 756 (2012) 145-153

Complex Structure in Class 0 Protostellar Envelopes III: Velocity Gradients in Non-Axisymmetric Envelopes, Infall or Rotation?

Authors: John J. Tobin, Lee Hartmann, Edwin A. Bergin, Hsin-Fang Chiang, Leslie W. Looney, Claire J. Chandler, Sebastien Maret, Fabian Heitsch

Abstract: We present an interferometric kinematic study of morphologically complex protostellar envelopes based on observations of the dense gas tracers N2H+ and NH3. The strong asymmetric nature of most envelopes in our sample leads us to question the common interpretation of velocity gradients as rotation, given the possibility of projection effects in the observed velocities. Several "idealized" sources… ▽ More We present an interferometric kinematic study of morphologically complex protostellar envelopes based on observations of the dense gas tracers N2H+ and NH3. The strong asymmetric nature of most envelopes in our sample leads us to question the common interpretation of velocity gradients as rotation, given the possibility of projection effects in the observed velocities. Several "idealized" sources with well-ordered velocity fields and envelope structures are now analyzed in more detail. We compare the interferometric data to position-velocity diagrams of kinematic models for spherical rotating collapse and filamentary rotating collapse. For this purpose, we developed a filamentary parametrization of the rotating collapse model to explore the effects of geometric projection on the observed velocity structures. We find that most envelopes in our sample have PV structures that can be reproduced by an infalling filamentary envelope projected at different angles within the plane of the sky. The infalling filament produces velocity shifts across the envelope that can mimic rotation, especially when viewed at single-dish resolutions and the axisymmetric rotating collapse model does not uniquely describe any dataset. Furthermore, if the velocities are assumed to reflect rotation, then the inferred centrifugal radii are quite large in most cases, indicating significant fragmentation potential or more likely another component to the line-center velocity. We conclude that ordered velocity gradients cannot be interpreted as rotation alone when envelopes are non-axisymmetric and that projected infall velocities likely dominate the velocity field on scales larger than 1000 AU. △ Less

Submitted 10 January, 2012; originally announced January 2012.

Comments: 37 pages, 15 Figures, 2 Tables, Accepted to ApJ

Very Large Array Observations of Ammonia in Infrared-Dark Clouds II: Internal Kinematics

Authors: Sarah E. Ragan, Fabian Heitsch, Edwin A. Bergin, David Wilner

Abstract: Infrared-dark clouds (IRDCs) are believed to be the birthplaces of rich clusters and thus contain the earliest phases of high-mass star formation. We use the Green Bank Telescope (GBT) and Very Large Array (VLA) maps of ammonia (NH3) in six IRDCs to measure their column density and temperature structure (Paper 1), and here, we investigate the kinematic structure and energy content. We find that IR… ▽ More Infrared-dark clouds (IRDCs) are believed to be the birthplaces of rich clusters and thus contain the earliest phases of high-mass star formation. We use the Green Bank Telescope (GBT) and Very Large Array (VLA) maps of ammonia (NH3) in six IRDCs to measure their column density and temperature structure (Paper 1), and here, we investigate the kinematic structure and energy content. We find that IRDCs overall display organized velocity fields, with only localized disruptions due to embedded star formation. The local effects seen in NH3 emission are not high velocity outflows but rather moderate (few km/s) increases in the line width that exhibit maxima near or coincident with the mid-infrared emission tracing protostars. These line width enhancements could be the result of infall or (hidden in NH3 emission) outflow. Not only is the kinetic energy content insufficient to support the IRDCs against collapse, but also the spatial energy distribution is inconsistent with a scenario of turbulent cloud support. We conclude that the velocity signatures of the IRDCs in our sample are due to active collapse and fragmentation, in some cases augmented by local feedback from stars. △ Less

Submitted 6 January, 2012; originally announced January 2012.

Comments: 15 pages, 12 figures, accepted for publication in ApJ

Rapid Star Formation and Global Gravitational Collapse

Authors: Lee Hartmann, Javier Ballesteros-Paredes, Fabian Heitsch

Abstract: Most young stars in nearby molecular clouds have estimated ages of 1-2 Myr, suggesting that star formation is rapid. However, small numbers of stars in these regions with inferred ages of >= 5-10 Myr have been cited to argue that star formation is instead a slow, quasi-static process. When considering these alternative pictures it is important to recognize that the age spread in a given star-formi… ▽ More Most young stars in nearby molecular clouds have estimated ages of 1-2 Myr, suggesting that star formation is rapid. However, small numbers of stars in these regions with inferred ages of >= 5-10 Myr have been cited to argue that star formation is instead a slow, quasi-static process. When considering these alternative pictures it is important to recognize that the age spread in a given star-forming cloud is necessarily an upper limit to the timescales of local collapse, as not all spatially-distinct regions will start contracting at precisely the same instant. Moreover, star-forming clouds may dynamically evolve on timescales of a few Myr; in particular, global gravitational contraction will tend to yield increasing star formation rates with time due to generally increasing local gas densities. We show that two different numerical simulations of dynamic, flow-driven molecular cloud formation and evolution 1) predict age spreads for the main stellar population roughly consistent with observations, and 2) raise the possibility of forming small numbers of stars early in cloud evolution, before global contraction concentrates the gas and the bulk of the stellar population is produced. In general, the existence of a small number of older stars among a generally much-younger population is consistent with the picture of dynamic star formation, and may even provide clues to the time evolution of star-forming clouds. △ Less

Submitted 10 November, 2011; originally announced November 2011.

Comments: 5 pages, 4 figures, accepted by MNRAS

Modes of Star Formation in Finite Molecular Clouds

Authors: A. Pon, D. Johnstone, F. Heitsch

Abstract: We analytically investigate the modes of gravity-induced star formation possible in idealized finite molecular clouds where global collapse competes against both local Jeans instabilities and discontinuous edge instabilities. We examine these timescales for collapse in spheres, discs, and cylinders, with emphasis on the structure, size, and degree of internal perturbations required in order for lo… ▽ More We analytically investigate the modes of gravity-induced star formation possible in idealized finite molecular clouds where global collapse competes against both local Jeans instabilities and discontinuous edge instabilities. We examine these timescales for collapse in spheres, discs, and cylinders, with emphasis on the structure, size, and degree of internal perturbations required in order for local collapse to occur before global collapse. We find that internal, local collapse is more effective for the lower dimensional objects. Spheres and discs, if unsupported against global collapse, must either contain strong perturbations or must be unrealistically large in order for small density perturbations to collapse significantly faster than the entire cloud. We find, on the other hand, that filamentary geometry is the most favorable situation for the smallest perturbations to grow before global collapse overwhelms them and that filaments containing only a few Jeans masses and weak density perturbations can readily fragment. These idealized solutions are compared with simulations of star-forming regions in an attempt to delineate the role of global, local, and edge instabilities in determining the fragmentation properties of molecular clouds. The combined results are also discussed in the context of recent observations of Galactic molecular clouds. △ Less

Submitted 5 August, 2011; originally announced August 2011.

Comments: 42 pages, 6 figures, 1 table, accepted by ApJ

Journal ref: ApJ 740 (2011) 88-97

Complex Structure in Class 0 Protostellar Envelopes II: Kinematic Structure from Single-Dish and Interferometric Molecular Line Mapping

Authors: John J. Tobin, Lee Hartmann, Hsin-Fang Chiang, Leslie W. Looney, Edwin A. Bergin, Claire J. Chandler, Josep M. Masque, Sebastien Maret, Fabian Heitsch

Abstract: We present a study of dense molecular gas kinematics in seventeen nearby protostellar systems using single-dish and interferometric molecular line observations. The non-axisymmetric envelopes around a sample of Class 0/I protostars were mapped in the N2H+ (J=1-0) tracer with the IRAM 30m, CARMA and PdBI as well as NH3 (1,1) with the VLA. The molecular line emission is used to construct line-center… ▽ More We present a study of dense molecular gas kinematics in seventeen nearby protostellar systems using single-dish and interferometric molecular line observations. The non-axisymmetric envelopes around a sample of Class 0/I protostars were mapped in the N2H+ (J=1-0) tracer with the IRAM 30m, CARMA and PdBI as well as NH3 (1,1) with the VLA. The molecular line emission is used to construct line-center velocity and linewidth maps for all sources to examine the kinematic structure in the envelopes on spatial scales from 0.1 pc to ~1000 AU. The direction of the large-scale velocity gradients from single-dish mapping is within 45 degrees of normal to the outflow axis in more than half the sample. Furthermore, the velocity gradients are often quite substantial, the average being ~2.3 km\s\pc. The interferometric data often reveal small-scale velocity structure, departing from the more gradual large-scale velocity gradients. In some cases, this likely indicates accelerating infall and/or rotational spin-up in the inner envelope; the median velocity gradient from the interferometric data is ~10.7 km/s/pc. In two systems, we detect high-velocity HCO+ (J=1-0) emission inside the highest-velocity \nthp\ emission. This enables us to study the infall and rotation close to the disk and estimate the central object masses. The velocity fields observed on large and small-scales are more complex than would be expected from rotation alone, suggesting that complex envelope structure enables other dynamical processes (i.e. infall) to affect the velocity field. △ Less

Submitted 21 July, 2011; originally announced July 2011.

Comments: 85 Pages, 31 Figures, 11 Tables, Accepted to ApJS

The Coalsack near and far

Authors: H. Beuther, J. Kainulainen, Th. Henning, R. Plume, F. Heitsch

Abstract: Context: The large Coalsack dark cloud is one of the rare molecular clouds without clear signs of star formation. Aims: We investigate the dynamical properties of the gas within the Coalsack. Methods: The two highest extinction regions were mapped with the APEX telescope in 13CO(2-1) comprising a region of ~1 square degree. Results: In addition to the well-known, nearby gas component around -4km/s… ▽ More Context: The large Coalsack dark cloud is one of the rare molecular clouds without clear signs of star formation. Aims: We investigate the dynamical properties of the gas within the Coalsack. Methods: The two highest extinction regions were mapped with the APEX telescope in 13CO(2-1) comprising a region of ~1 square degree. Results: In addition to the well-known, nearby gas component around -4km/s, we identified additional molecular gas components - in particular a second extended molecular cloud at a velocity of ~-30km/s and an estimated distance of ~3.1kpc - that dominate the column density and visual extinction distributions in the northeastern part of the Coalsack. Although comprising ~2600Msun, the mass of this distant cloud is distributed over an extent of ~73pc, much larger than typical high-mass infrared dark clouds. Its filamentary structure is consistent with a compressible gaseous self-gravitating cylinder, and its low mass per length indicates that it may be stable against gravitational collapse. We find barely any mid-infrared emission in archival MSX data, which is indicative of almost no star-formation activity in the near and far cloud complexes. The nearby clouds have narrow, almost thermal velocity dispersions with median values between 0.2 and 0.4km/s, which is also consistent with low star-formation activity. Only Tapia's Globule 2 exhibits a velocity dispersion increase toward the extinction peak and peak-velocity gradients over the core, which is indicative of a state of elevated dynamical properties. Conclusions: The Coalsack is not one single coherent structure, but consists of several cloud complexes nearby as well as at several kpc distance. All studied clouds appear as starless low-turbulence regions that may not even collapse in the future. Only one globule exhibits more dynamical signatures and is a good candidate for present/future star formation. △ Less

Submitted 12 July, 2011; originally announced July 2011.

Comments: 9 pages, 11 figures accepted for A&A, available also at http://www.mpia-hd.mpg.de/homes/beuther/papers.html

Gravity or turbulence? II. Evolving column density PDFs in molecular clouds

Authors: Javier Ballesteros-Paredes, Enrique Vazquez-Semadeni, Adriana Gazol, Lee W. Hartmann, Fabian Heitsch, Pedro Colin

Abstract: 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 motio… ▽ More 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. △ Less

Submitted 26 May, 2011; originally announced May 2011.

Comments: 9 pages, 4 figures. MNRAS, accepted

Flow-Driven Cloud Formation and Fragmentation: Results From Eulerian and Lagrangian Simulations

Authors: Fabian Heitsch, Thorsten Naab, Stefanie Walch

Abstract: The fragmentation of shocked flows in a thermally bistable medium provides a natural mechanism to form turbulent cold clouds as precursors to molecular clouds. Yet because of the large density and temperature differences and the range of dynamical scales involved, following this process with numerical simulations is challenging. We compare two-dimensional simulations of flow-driven cloud formation… ▽ More The fragmentation of shocked flows in a thermally bistable medium provides a natural mechanism to form turbulent cold clouds as precursors to molecular clouds. Yet because of the large density and temperature differences and the range of dynamical scales involved, following this process with numerical simulations is challenging. We compare two-dimensional simulations of flow-driven cloud formation without self-gravity, using the Lagrangian Smoothed Particle Hydrodynamics (SPH) code VINE and the Eulerian grid code Proteus. Results are qualitatively similar for both methods, yet the variable spatial resolution of the SPH method leads to smaller fragments and thinner filaments, rendering the overall morphologies different. Thermal and hydro-dynamical instabilities lead to rapid cooling and fragmentation into cold clumps with temperatures below 300K. For clumps more massive than 1 Msun/pc, the clump mass function has an average slope of -0.8. The internal velocity dispersion of the clumps is nearly an order of magnitude smaller than their relative motion, rendering it subsonic with respect to the internal sound speed of the clumps, but supersonic as seen by an external observer. For the SPH simulations most of the cold gas resides at temperatures below 100K, while the grid-based models show an additional, substantial component between 100 and 300K. Independently of the numerical method our models confirm that converging flows of warm neutral gas fragment rapidly and form high-density, low-temperature clumps as possible seeds for star formation. △ Less

Submitted 11 March, 2011; originally announced March 2011.

Comments: 9 pages, 8 figures, MNRAS accepted

W43: the closest molecular complex of the Galactic Bar?

Authors: Q. Nguyen Luong, F. Motte, F. Schuller, N. Schneider, S. Bontemps, P. Schilke, K. M. Menten, F. Heitsch, F. Wyrowski, P. Carlhoff, L. Bronfman, T. Henning

Abstract: We used a large database extracted from Galaxy-wide surveys of H {\scriptsize I}, $^{13}$CO 1-0, $8 \micron$ and $870 \micron$ continuum to trace diffuse atomic gas, low- to medium-density molecular gas, high-density molecular gas, and star formation activity which we complemented by dedicated $^{12}$CO 2--1, 3--2 observations of the W43 region. From the detailed 3D (space-space-velocity) analysis… ▽ More We used a large database extracted from Galaxy-wide surveys of H {\scriptsize I}, $^{13}$CO 1-0, $8 \micron$ and $870 \micron$ continuum to trace diffuse atomic gas, low- to medium-density molecular gas, high-density molecular gas, and star formation activity which we complemented by dedicated $^{12}$CO 2--1, 3--2 observations of the W43 region. From the detailed 3D (space-space-velocity) analysis of the molecular and atomic cloud tracers through the region and despite its wide velocity range (\emph{FWHM}$\sim22.3 \kms$ around $\vlsr\sim 95.9 \kms$), we identified W43 as a large (equivalent diameter $\sim 140$ pc) and coherent complex of molecular clouds which is surrounded by an atomic gas envelope (equivalent diameter $\sim 290$ pc). We measured the total mass of this newly-identified molecular complex ($M_{\tiny total}\sim 7.1 \times 10^6 \msun$), the mass contained in dense $870 \micron$ clumps ($<5$ pc dense cloud structures, $M_{\tiny clumps}\sim 8.4 \times 10^5 \msun$) and conclude that W43 is particularly massive and concentrated. The distance we assume for the W43 complex is 6 kpc from the Sun, which may place it at the meeting point of the Scutum-Centaurus (or Scutum-Crux) Galactic arm and the Bar, a dynamically complex region where high-velocity streams could easily collide. The star formation rate of W43 is suggested not to be steady but it is increasing from $\, \sim 0.01 \msun\, {yr}^{-1}$ (measured from its $8\,μ$m luminosity) to $\sim 0.1 \msun\, {yr}^{-1}$ (measured from its molecular content). From the global properties of W43, we claim it is an extreme molecular complex in the Milky Way and it could even be forming starburst clusters in the near future. △ Less

Submitted 19 February, 2011; v1 submitted 16 February, 2011; originally announced February 2011.

Comments: 16 pages, accepted for A&A (14 Feb 2011)

Formation of Cold Filamentary Structure from Wind Blown Superbubbles

Authors: Evangelia Ntormousi, Andreas Burkert, Katharina Fierlinger, Fabian Heitsch

Abstract: The expansion and collision of two wind-blown superbubbles is investigated numerically. Our models go beyond previous simulations of molecular cloud formation from converging gas flows by exploring this process with realistic flow parameters, sizes and timescales. The superbubbles are blown by time-dependent winds and supernova explosions, calculated from population synthesis models. They expand i… ▽ More The expansion and collision of two wind-blown superbubbles is investigated numerically. Our models go beyond previous simulations of molecular cloud formation from converging gas flows by exploring this process with realistic flow parameters, sizes and timescales. The superbubbles are blown by time-dependent winds and supernova explosions, calculated from population synthesis models. They expand into a uniform or turbulent diffuse medium. We find that dense, cold gas clumps and filaments form naturally in the compressed collision zone of the two superbubbles. Their shapes resemble the elongated, irregular structure of observed cold, molecular gas filaments and clumps. At the end of the simulations, between 65 and 80 percent of the total gas mass in our simulation box is contained in these structures. The clumps are found in a variety of physical states, ranging from pressure equilibrium with the surrounding medium to highly under-pressured clumps with large irregular internal motions and structures which are rotationally supported. △ Less

Submitted 26 November, 2010; originally announced November 2010.

Comments: Submitted to ApJ

Physical Properties of Complex C Halo Clouds

Authors: W. -H. Hsu, M. E. Putman, F. Heitsch, S. Stanimirovic, J. E. G. Peek, S. E. Clark

Abstract: Observations from the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) Survey of the tail of Complex C are presented and the halo clouds associated with this complex cataloged. The properties of the Complex C clouds are compared to clouds cataloged at the tail of the Magellanic Stream to provide insight into the origin and destruction mechanism of Complex C. Magellanic Stream and Complex C clouds… ▽ More Observations from the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) Survey of the tail of Complex C are presented and the halo clouds associated with this complex cataloged. The properties of the Complex C clouds are compared to clouds cataloged at the tail of the Magellanic Stream to provide insight into the origin and destruction mechanism of Complex C. Magellanic Stream and Complex C clouds show similarities in their mass distributions (slope = -0.7 and -0.6, respectively) and have a common linewidth of 20 - 30 km/s (indicative of a warm component), which may indicate a common origin and/or physical process breaking down the clouds. The clouds cataloged at the tail of Complex C extend over a mass range of 10^1.1 to 10^4.8 solar masses, sizes of 10^1.2 to 10^2.6 pc, and have a median volume density of 0.065 cm^(-3) and median pressure of (P/k) = 580 K cm^{-3}. We do not see a prominent two-phase structure in Complex C, possibly due to its low metallicity and inefficient cooling compared to other halo clouds. From assuming the Complex C clouds are in pressure equilibrium with a hot halo medium, we find a median halo density of 5.8 x 10^(-4) cm^(-3), which given a constant distance of 10 kpc, is at a z-height of ~3 kpc. Using the same argument for the Stream results in a median halo density of 8.4 x 10^(-5) x (60kpc/d) cm^(-3). These densities are consistent with previous observational constraints and cosmological simulations. We also assess the derived cloud and halo properties with three dimensional grid simulations of halo HI clouds and find the temperature is generally consistent within a factor of 1.5 and the volume densities, pressures and halo densities are consistent within a factor of 3. △ Less

Submitted 29 October, 2010; originally announced November 2010.

Comments: Accepted for publication in AJ. 54 pages, including 6 tables and 16 figures

Journal ref: The Astronomical Journal, Volume 141, Issue 2, article id. 57 (2011)

Gravity or Turbulence? The velocity dispersion-size relation

Authors: Javier Ballesteros-Paredes, Lee W. Hartmann, Enrique Vázquez-Semadeni, Fabian Heitsch, Manuel A. Zamora-Avilés

Abstract: We discuss the nature of the velocity dispersion vs. size relation for molecular clouds. In particular, we add to previous observational results showing that the velocity dispersions in molecular clouds and cores are not purely functions of spatial scale but involve surface gas densities as well. We emphasize that hydrodynamic turbulence is required to produce the first condensations in the progen… ▽ More We discuss the nature of the velocity dispersion vs. size relation for molecular clouds. In particular, we add to previous observational results showing that the velocity dispersions in molecular clouds and cores are not purely functions of spatial scale but involve surface gas densities as well. We emphasize that hydrodynamic turbulence is required to produce the first condensations in the progenitor medium. However, as the cloud is forming, it also becomes bound, and gravitational accelerations dominate the motions. Energy conservation in this case implies $|E_g| \sim E_k$, in agreement with observational data, and providing an interpretation for two recent observational results: the scatter in the $δv-R$ plane, and the dependence of the velocity dispersion on the surface density ${δv^2/ R} \propto Σ$. We argue that the observational data are consistent with molecular clouds in a state of hierarchical gravitational collapse, i.e., developing local centers of collapse throughout the whole cloud while the cloud itself is collapsing, and making equilibrium unnecessary at all stages prior to the formation of actual stars. Finally, we discuss how this mechanism need not be in conflict with the observed star formation rate. △ Less

Submitted 8 September, 2010; originally announced September 2010.

Comments: Accepted by MNRAS. 7 pages, 3 figures

Morphological Complexity of Protostellar Envelopes

Authors: John J. Tobin, Lee Hartmann, Edwin Bergin, Leslie W. Looney, Hsin-Fang Chiang, Fabian Heitsch

Abstract: Extinction maps at 8 micron from the Spitzer Space Telescope show that many Class 0 protostars exhibit complex, irregular, and on-axisymmetric structure within the densest regions of their dusty envelopes. Many of the systems have highly irregular and on-axisymmetric morphologies on scales $\sim$1000 AU, with a quarter of the sample exhibiting filamentary or flattened dense structures. Complex env… ▽ More Extinction maps at 8 micron from the Spitzer Space Telescope show that many Class 0 protostars exhibit complex, irregular, and on-axisymmetric structure within the densest regions of their dusty envelopes. Many of the systems have highly irregular and on-axisymmetric morphologies on scales $\sim$1000 AU, with a quarter of the sample exhibiting filamentary or flattened dense structures. Complex envelope structure is observed in regions spatially distinct from outflow cavities, and the densest structures often show no systematic alignment perpendicular to the cavities. We suggest that the observed envelope complexity is the result of collapse from protostellar cores with initially non-equilibrium structures. The striking non-axisymmetry in many envelopes could provide favorable conditions for the formation of binary systems. We then show that the kinematics around L1165 as probed with N2H+ are indicative of asymmetric infall; the velocity gradient is not perpendicular to the outflow. △ Less

Submitted 19 August, 2010; originally announced August 2010.

Comments: 4 pages; To appear in the proceedings for IAU Symposium 270: Computational Star Formation

Competitive Accretion in Sheet Geometry and the Stellar IMF

Authors: Wen-hsin Hsu, Lee Hartmann, Fabian Heitsch, Gilberto C. Gómez

Abstract: We report a set of numerical experiments aimed at addressing the applicability of competitive accretion to explain the high-mass end of the stellar initial mass function in a sheet geometry with shallow gravitational potential, in contrast to most previous simulations which have assumed formation in a cluster gravitational potential. Our flat cloud geometry is motivated by models of molecular clou… ▽ More We report a set of numerical experiments aimed at addressing the applicability of competitive accretion to explain the high-mass end of the stellar initial mass function in a sheet geometry with shallow gravitational potential, in contrast to most previous simulations which have assumed formation in a cluster gravitational potential. Our flat cloud geometry is motivated by models of molecular cloud formation due to large-scale flows in the interstellar medium. The experiments consisted of SPH simulations of gas accretion onto sink particles formed rapidly from Jeans-unstable dense clumps placed randomly in the finite sheet. These simplifications allow us to study accretion with a minimum of free parameters, and to develop better statistics on the resulting mass spectra. We considered both clumps of equal mass and gaussian distributions of masses, and either uniform or spatially-varying gas densities. In all cases, the sink mass function develops a power law tail at high masses, with $dN/dlog M \propto M^{-Γ}$. The accretion rates of individual sinks follow $\dot{M} \propto M^2$ at high masses; this results in a continual flattening of the slope of the mass function towards an asymptotic form $Γ\sim 1$ (where the Salpeter slope is $Γ= 1.35$). The asymptotic limit is most rapidly reached when starting from a relatively broad distribution of initial sink masses. In general the resulting upper mass slope is correlated with the maximum sink mass; higher sink masses are found in simulations with flatter upper mass slopes. Although these simulations are of a highly idealized situation, the results suggest that competitive accretion may be relevant in a wider variety of environments than previously considered, and in particular that the upper mass distribution may generally evolve towards a limiting value of $Γ\sim 1$. △ Less

Submitted 18 August, 2010; originally announced August 2010.

Comments: 20 pages, 12 figures

Spitzer View of Young Massive Stars in the LMC HII Complexes. II. N159

Authors: C. -H. Rosie Chen, Remy Indebetouw, You-Hua Chu, Robert A. Gruendl, Gerard Testor, Fabian Heitsch, Jonathan P. Seale, Margaret Meixner, Marta Sewilo

Abstract: The HII complex N159 in the Large Magellanic Cloud (LMC) is used to study massive star formation in different environments, as it contains three giant molecular clouds (GMCs) that have similar sizes and masses but exhibit different intensities of star formation. We identify candidate massive young stellar objects (YSOs) using infrared photometry, and model their SEDs to constrain mass and evolutio… ▽ More The HII complex N159 in the Large Magellanic Cloud (LMC) is used to study massive star formation in different environments, as it contains three giant molecular clouds (GMCs) that have similar sizes and masses but exhibit different intensities of star formation. We identify candidate massive young stellar objects (YSOs) using infrared photometry, and model their SEDs to constrain mass and evolutionary state. Good fits are obtained for less evolved Type I, I/II, and II sources. Our analysis suggests that there are massive embedded YSOs in N159B, a maser source, and several ultracompact HII regions. Massive O-type YSOs are found in GMCs N159-E and N159-W, which are associated with ionized gas, i.e., where massive stars formed a few Myr ago. The third GMC, N159-S, has neither O-type YSOs nor evidence of previous massive star formation. This correlation between current and antecedent formation of massive stars suggests that energy feedback is relevant. We present evidence that N159-W is forming YSOs spontaneously, while collapse in N159-E may be triggered. Finally, we compare star formation rates determined from YSO counts with those from integrated H-alpha and 24 micron luminosities and expected from gas surface densities. Detailed dissection of extragalactic GMCs like the one presented here is key to revealing the physics underlying commonly used star formation scaling laws. △ Less

Submitted 29 July, 2010; originally announced July 2010.

Comments: 60 pages, 11 figures. Accepted for publication in Astrophysical Journal

Journal ref: ApJ (2010) 721, 1206-1232

Modelling Shear Flows with SPH and Grid Based Methods

Authors: Veronika Junk, Stefanie Walch, Fabian Heitsch, Andreas Burkert, Markus Wetzstein, Marc Schartmann, Daniel Price

Abstract: Given the importance of shear flows for astrophysical gas dynamics, we study the evolution of the Kelvin-Helmholtz instability (KHI) analytically and numerically. We derive the dispersion relation for the two-dimensional KHI including viscous dissipation. The resulting expression for the growth rate is then used to estimate the intrinsic viscosity of four numerical schemes depending on code-specif… ▽ More Given the importance of shear flows for astrophysical gas dynamics, we study the evolution of the Kelvin-Helmholtz instability (KHI) analytically and numerically. We derive the dispersion relation for the two-dimensional KHI including viscous dissipation. The resulting expression for the growth rate is then used to estimate the intrinsic viscosity of four numerical schemes depending on code-specific as well as on physical parameters. Our set of numerical schemes includes the Tree-SPH code VINE, an alternative SPH formulation developed by Price (2008), and the finite-volume grid codes FLASH and PLUTO. In the first part, we explicitly demonstrate the effect of dissipation-inhibiting mechanisms such as the Balsara viscosity on the evolution of the KHI. With VINE, increasing density contrasts lead to a continuously increasing suppression of the KHI (with complete suppression from a contrast of 6:1 or higher). The alternative SPH formulation including an artificial thermal conductivity reproduces the analytically expected growth rates up to a density contrast of 10:1. The second part addresses the shear flow evolution with FLASH and PLUTO. Both codes result in a consistent non-viscous evolution (in the equal as well as in the different density case) in agreement with the analytical prediction. The viscous evolution studied with FLASH shows minor deviations from the analytical prediction. △ Less

Submitted 12 April, 2010; originally announced April 2010.

Comments: 16 pages, 17 figures

Spitzer and HHT observations of starless cores: masses and environments

Authors: Amelia M. Stutz, George H. Rieke, John H. Bieging, Zoltan Balog, Fabian Heitsch, Miju Kang, William L. Peters, Yancy L. Shirley, Michael W. Werner

Abstract: We present Spitzer observations of a sample of 12 starless cores selected to have prominent 24 micron shadows. The Spitzer images show 8 and 24 micron shadows and in some cases 70 micron shadows; these spatially resolved absorption features trace the densest regions of the cores. We have carried out a 12CO (2-1) and 13CO (2-1) mapping survey of these cores with the Heinrich Hertz Telescope (HHT)… ▽ More We present Spitzer observations of a sample of 12 starless cores selected to have prominent 24 micron shadows. The Spitzer images show 8 and 24 micron shadows and in some cases 70 micron shadows; these spatially resolved absorption features trace the densest regions of the cores. We have carried out a 12CO (2-1) and 13CO (2-1) mapping survey of these cores with the Heinrich Hertz Telescope (HHT). We use the shadow features to derive optical depth maps. We derive molecular masses for the cores and the surrounding environment; we find that the 24 micron shadow masses are always greater than or equal to the molecular masses derived in the same region, a discrepancy likely caused by CO freeze--out onto dust grains. We combine this sample with two additional cores that we studied previously to bring the total sample to 14 cores. Using a simple Jeans mass criterion we find that ~ 2/3 of the cores selected to have prominent 24 micron shadows are collapsing or near collapse, a result that is supported by millimeter line observations. Of this subset at least half have indications of 70 micron shadows. All cores observed to produce absorption features at 70 micron are close to collapse. We conclude that 24 micron shadows, and even more so the 70 micron ones, are useful markers of cloud cores that are approaching collapse. △ Less

Submitted 3 November, 2009; v1 submitted 29 October, 2009; originally announced October 2009.

Comments: 41 pages, 28 figures, 5 tables; accepted by ApJ

Journal ref: Astrophys.J.707:137-166,2009

Gravitational Collapse and Filament Formation: Comparison with the Pipe Nebula

Authors: Fabian Heitsch, Javier Ballesteros-Paredes, Lee Hartmann

Abstract: Recent models of molecular cloud formation and evolution suggest that such clouds are dynamic and generally exhibit gravitational collapse. We present a simple analytic model of global collapse onto a filament and compare this with our numerical simulations of the flow-driven formation of an isolated molecular cloud to illustrate the supersonic motions and infall ram pressures expected in models… ▽ More Recent models of molecular cloud formation and evolution suggest that such clouds are dynamic and generally exhibit gravitational collapse. We present a simple analytic model of global collapse onto a filament and compare this with our numerical simulations of the flow-driven formation of an isolated molecular cloud to illustrate the supersonic motions and infall ram pressures expected in models of gravity-driven cloud evolution. We apply our results to observations of the Pipe Nebula, an especially suitable object for our purposes as its low star formation activity implies insignifcant perturbations from stellar feedback. We show that our collapsing cloud model can explain the magnitude of the velocity dispersions seen in the $^{13}$CO filamentary structure by Onishi et al. and the ram pressures required by Lada et al. to confine the lower-mass cores in the Pipe nebula. We further conjecture that higher-resolution simulations will show small velocity dispersions in the densest core gas, as observed, but which are infall motions and not supporting turbulence. Our results point out the inevitability of ram pressures as boundary conditions for molecular cloud filaments, and the possibility that especially lower-mass cores still can be accreting mass at significant rates, as suggested by observations. △ Less

Submitted 10 September, 2009; originally announced September 2009.

Comments: 10 pages, 4 figures, accepted by ApJ

Journal ref: Astrophys.J.704:1735-1742,2009

The Accretion of Fuel at the Disk-Halo Interface

Authors: M. E. Putman, J. E. G. Peek, F. Heitsch

Abstract: We discuss the support for the cooling of gas directly at the disk-halo interface as a mechanism to continually fuel a galaxy. This may be an important reservoir as there is not enough cold gas observed in galaxy halos and simulations indicate the existing cold clouds will be rapidly destroyed as they move through the surrounding halo medium. We show possible evidence for a net infall of the WIM… ▽ More We discuss the support for the cooling of gas directly at the disk-halo interface as a mechanism to continually fuel a galaxy. This may be an important reservoir as there is not enough cold gas observed in galaxy halos and simulations indicate the existing cold clouds will be rapidly destroyed as they move through the surrounding halo medium. We show possible evidence for a net infall of the WIM layer in the Milky Way, simulation results showing the recooling of warm clouds at the disk-halo interface, and GALFA HI data of small, cold HI clouds that could represent this recooling. △ Less

Submitted 6 July, 2009; originally announced July 2009.

Comments: 7.5 pages, 4 figures, invited review for The Role of the Disk-Halo Interaction in Galaxy Evolution: Outflow vs. Infall?, ed. M. A. de Avillez, European Astronomical Society Publications, 2009 in press

Numerical Star-Formation Studies -- A Status Report

Authors: Ralf S. Klessen, Mark R. Krumholz, Fabian Heitsch

Abstract: The formation of stars is a key process in astrophysics. Detailed knowledge of the physical mechanisms that govern stellar birth is a prerequisite for understanding the formation and evolution of our galactic home, the Milky Way. A theory of star formation is an essential part of any model for the origin of our solar system and of planets around other stars. Despite this pivotal importance, and… ▽ More The formation of stars is a key process in astrophysics. Detailed knowledge of the physical mechanisms that govern stellar birth is a prerequisite for understanding the formation and evolution of our galactic home, the Milky Way. A theory of star formation is an essential part of any model for the origin of our solar system and of planets around other stars. Despite this pivotal importance, and despite many decades of research, our understanding of the processes that initiate and regulate star formation is still limited. Stars are born in cold interstellar clouds of molecular hydrogen gas. Star formation in these clouds is governed by the complex interplay between the gravitational attraction in the gas and agents such as turbulence, magnetic fields, radiation and thermal pressure that resist compression. The competition between these processes determines both the locations at which young stars form and how much mass they ultimately accrete. It plays out over many orders of magnitude in space and time, ranging from galactic to stellar scales. In addition, star formation is a highly stochastic process in which rare and hard-to-predict events, such as the formation of very massive stars and the resulting feedback, can play a dominant role in determining the evolution of a star-forming cloud. As a consequence of the wide range of scales and processes that control star formation, analytic models are usually restricted to highly idealized cases. These can yield insight, but the complexity of the problem means that they must be used in concert with large-scale numerical simulations. Here we summarize the state of modern star formation theory and review the recent advances in numerical simulation techniques. △ Less

Submitted 24 June, 2009; originally announced June 2009.

Comments: Invited Review to appear on Advanced Science Letters (ASL), Special Issue on Computational Astrophysics, edited by Lucio Mayer (32 pages, 10 figures)

Journal ref: Advanced Science Letters, 4, 258-285 (2011)