-
The AGORA high-resolution galaxy simulations comparison project: Public data release
Authors:
Santi Roca-Fàbrega,
Ji-hoon Kim,
Joel R. Primack,
Michael J. Butler,
Daniel Ceverino,
Jun-Hwan Choi,
Robert Feldmann,
Ben W. Keller,
Alessandro Lupi,
Kentaro Nagamine,
Thomas R. Quinn,
Yves Revaz,
Romain Teyssier,
Spencer C. Wallace
Abstract:
As part of the AGORA High-resolution Galaxy Simulations Comparison Project (Kim et al. 2014, 2016) we have generated a suite of isolated Milky Way-mass galaxy simulations using 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical galaxy formation community. In these simulations we adopted identical galactic disk initial conditions, and common physics models (e.g., radiative…
▽ More
As part of the AGORA High-resolution Galaxy Simulations Comparison Project (Kim et al. 2014, 2016) we have generated a suite of isolated Milky Way-mass galaxy simulations using 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical galaxy formation community. In these simulations we adopted identical galactic disk initial conditions, and common physics models (e.g., radiative cooling and ultraviolet background by a standardized package). Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production were carefully constrained. Here we release the simulation data to be freely used by the community. In this release we include the disk snapshots at 0 and 500Myr of evolution per each code as used in Kim et al. (2016), from simulations with and without star formation and feedback. We encourage any member of the numerical galaxy formation community to make use of these resources for their research - for example, compare their own simulations with the AGORA galaxies, with the common analysis yt scripts used to obtain the plots shown in our papers, also available in this release.
△ Less
Submitted 16 January, 2020; v1 submitted 13 January, 2020;
originally announced January 2020.
-
Core Emergence in a Massive Infrared Dark Cloud: A Comparison Between Mid-IR Extinction and 1.3 mm Emission
Authors:
Shuo Kong,
Jonathan C. Tan,
Héctor G. Arce,
Paola Caselli,
Francesco Fontani,
Michael J. Butler
Abstract:
Stars are born from dense cores in molecular clouds. Observationally, it is crucial to capture the formation of cores in order to understand the necessary conditions and rate of the star formation process. The {\it Atacama Large Mm/sub-mm Array} (ALMA) is extremely powerful for identifying dense gas structures, including cores, at mm wavelengths via their dust continuum emission. Here we use ALMA…
▽ More
Stars are born from dense cores in molecular clouds. Observationally, it is crucial to capture the formation of cores in order to understand the necessary conditions and rate of the star formation process. The {\it Atacama Large Mm/sub-mm Array} (ALMA) is extremely powerful for identifying dense gas structures, including cores, at mm wavelengths via their dust continuum emission. Here we use ALMA to carry out a survey of dense gas and cores in the central region of the massive ($\sim10^5\:M_\odot$) Infrared Dark Cloud (IRDC) G28.37+0.07. The observation consists of a mosaic of 86 pointings of the 12m-array and produces an unprecedented view of the densest structures of this IRDC. In this first paper about this data set, we focus on a comparison between the 1.3 mm continuum emission and a mid-infrared (MIR) extinction map of the IRDC. This allows estimation of the "dense gas" detection probability function (DPF), i.e., as a function of the local mass surface density, $Σ$, for various choices of thresholds of mm continuum emission to define "dense gas". We then estimate the dense gas mass fraction, $f_{\rm dg}$, in the central region of the IRDC and, via extrapolation with the DPF and the known $Σ$ probability distribution function, to the larger-scale surrounding regions, finding values of about 5\% to 15\% for the fiducial choice of threshold. We argue that this observed dense gas is a good tracer of the protostellar core population and, in this context, estimate a star formation efficiency per free-fall time in the central IRDC region of $ε_{\rm ff}\sim$10\%, with approximately a factor of two systematic uncertainties.
△ Less
Submitted 20 February, 2018;
originally announced February 2018.
-
Kiloparsec-Scale Simulations of Star Formation in Disk Galaxies. IV. Regulation of Galactic Star Formation Rates by Stellar Feedback
Authors:
Michael J. Butler,
Jonathan C. Tan,
Romain Teyssier,
Joakim Rosdahl,
Sven Van Loo,
Sarah Nickerson
Abstract:
Star formation from the interstellar medium of galactic disks is a basic process controlling the evolution of galaxies. Understanding the star formation rate in a local patch of a disk with a given gas mass is thus an important challenge for theoretical models. Here we simulate a kiloparsec region of a disk, following the evolution of self-gravitating molecular clouds down to subparsec scales, as…
▽ More
Star formation from the interstellar medium of galactic disks is a basic process controlling the evolution of galaxies. Understanding the star formation rate in a local patch of a disk with a given gas mass is thus an important challenge for theoretical models. Here we simulate a kiloparsec region of a disk, following the evolution of self-gravitating molecular clouds down to subparsec scales, as they form stars that then inject feedback energy by dissociating and ionizing UV photons and supernova explosions. We assess the relative importance of each feedback mechanism. We find that $\rm H_2$-dissociating feedback results in the largest absolute reduction in star formation compared to the run with no feedback. Subsequently adding photoionization feedback produces a more modest reduction. Our fiducial models that combine all three feedback mechanisms yield, without fine-tuning, star formation rates that are in excellent agreement with observations, with $\rm H_2$-dissociating photons playing a crucial role. Models that only include supernova feedback---a common method in galaxy evolution simulations---settle to similar star formation rates, but with very different temperature and chemical states of the gas, and with very different spatial distributions of young stars.
△ Less
Submitted 28 April, 2017; v1 submitted 13 March, 2017;
originally announced March 2017.
-
Zooming in to Massive Star Birth
Authors:
Shuo Kong,
Jonathan C. Tan,
Paola Caselli,
Francesco Fontani,
Ke Wang,
Michael J. Butler
Abstract:
We present high resolution (0.2", 1000 AU) 1.3 mm ALMA observations of massive infrared dark cloud clump, G028.37+00.07-C1, thought to harbor the early stages of massive star formation. Using $\rm N_2D^+$(3-2) we resolve the previously identified C1-S core, separating the bulk of its emission from two nearby protostellar sources. C1-S is thus identified as a massive ($\sim50\:M_\odot$), compact (…
▽ More
We present high resolution (0.2", 1000 AU) 1.3 mm ALMA observations of massive infrared dark cloud clump, G028.37+00.07-C1, thought to harbor the early stages of massive star formation. Using $\rm N_2D^+$(3-2) we resolve the previously identified C1-S core, separating the bulk of its emission from two nearby protostellar sources. C1-S is thus identified as a massive ($\sim50\:M_\odot$), compact ($\sim0.1\:$pc diameter) starless core, e.g., with no signs of outflow activity. Being highly deuterated, this is a promising candidate for a pre-stellar core on the verge of collapse. An analysis of its dynamical state indicates a sub-virial velocity dispersion compared to a trans-Alfvénic turbulent core model. However, virial equilibrium could be achieved with sub-Alfvénic conditions involving $\sim2\:$mG magnetic field strengths.
△ Less
Submitted 16 August, 2018; v1 submitted 20 January, 2017;
originally announced January 2017.
-
The AGORA High-Resolution Galaxy Simulations Comparison Project. II: Isolated Disk Test
Authors:
Ji-hoon Kim,
Oscar Agertz,
Romain Teyssier,
Michael J. Butler,
Daniel Ceverino,
Jun-Hwan Choi,
Robert Feldmann,
Ben W. Keller,
Alessandro Lupi,
Thomas Quinn,
Yves Revaz,
Spencer Wallace,
Nickolay Y. Gnedin,
Samuel N. Leitner,
Sijing Shen,
Britton D. Smith,
Robert Thompson,
Matthew J. Turk,
Tom Abel,
Kenza S. Arraki,
Samantha M. Benincasa,
Sukanya Chakrabarti,
Colin DeGraf,
Avishai Dekel,
Nathan J. Goldbaum
, et al. (18 additional authors not shown)
Abstract:
Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the st…
▽ More
Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.
△ Less
Submitted 8 March, 2018; v1 submitted 10 October, 2016;
originally announced October 2016.
-
A Hunt for Massive Starless Cores
Authors:
Shuo Kong,
Jonathan C. Tan,
Paola Caselli,
Francesco Fontani,
Mengyao Liu,
Michael J. Butler
Abstract:
We carry out an ALMA $\rm N_2D^+$(3-2) and 1.3~mm continuum survey towards 32 high mass surface density regions in seven Infrared Dark Clouds with the aim of finding massive starless cores, which may be the initial conditions for the formation of massive stars. Cores showing strong $\rm N_2D^+$(3-2) emission are expected to be highly deuterated and indicative of early, potentially pre-stellar stag…
▽ More
We carry out an ALMA $\rm N_2D^+$(3-2) and 1.3~mm continuum survey towards 32 high mass surface density regions in seven Infrared Dark Clouds with the aim of finding massive starless cores, which may be the initial conditions for the formation of massive stars. Cores showing strong $\rm N_2D^+$(3-2) emission are expected to be highly deuterated and indicative of early, potentially pre-stellar stages of star formation. We also present maps of these regions in ancillary line tracers, including C$^{18}$O(2-1), DCN(3-2) and DCO$^+$(3-2). Over 100 $\rm N_2D^+$ cores are identified with our newly developed core-finding algorithm based on connected structures in position-velocity space. The most massive core has $\sim70\:M_\odot$ (potentially $\sim170\:M_\odot$) and so may be representative of the initial conditions for massive star formation. The existence and dynamical properties of such cores constrain massive star formation theories. We measure the line widths and thus velocity dispersion of six of the cores with strongest $\rm N_2D^+$(3-2) line emission, finding results that are generally consistent with virial equilibrium of pressure confined cores.
△ Less
Submitted 28 November, 2016; v1 submitted 19 September, 2016;
originally announced September 2016.
-
Mid-J CO Shock Tracing Observations of Infrared Dark Clouds III: SLED fitting
Authors:
A. Pon,
M. J. Kaufman,
D. Johnstone,
P. Caselli,
F. Fontani,
M. J. Butler,
I. Jiménez-Serra,
A. Palau,
J. C. Tan
Abstract:
Giant molecular clouds contain supersonic turbulence that can locally heat small fractions of gas to over 100 K. We run shock models for low-velocity, C-type shocks propagating into gas with densities between 10^3 and 10^5 cm^(-3) and find that CO lines are the most important cooling lines. Comparison to photodissociation region (PDR) models indicates that mid-J CO lines (J = 8-7 and higher) shoul…
▽ More
Giant molecular clouds contain supersonic turbulence that can locally heat small fractions of gas to over 100 K. We run shock models for low-velocity, C-type shocks propagating into gas with densities between 10^3 and 10^5 cm^(-3) and find that CO lines are the most important cooling lines. Comparison to photodissociation region (PDR) models indicates that mid-J CO lines (J = 8-7 and higher) should be dominated by emission from shocked gas. In Papers I and II we presented CO J = 3-2, 8-7, and 9-8 observations toward four primarily quiescent clumps within infrared dark clouds. Here we fit PDR models to the combined spectral line energy distributions and show that the PDR models that best fit the low-J CO emission underpredict the mid-J CO emission by orders of magnitude, strongly hinting at a hot gas component within these clumps. The low-J CO data clearly show that the integrated intensities of both the CO J = 8-7 and 9-8 lines are anomalously high, such that the line ratio can be used to characterize the hot gas component. Shock models are reasonably consistent with the observed mid-J CO emission, with models with densities near 10^(4.5) cm^(-3) providing the best agreement. Where this mid-J CO is detected, the mean volume filling factor of the hot gas is 0.1%. Much of the observed mid-J CO emission, however, is also associated with known protostars and may be due to protostellar feedback.
△ Less
Submitted 25 July, 2016; v1 submitted 9 June, 2016;
originally announced June 2016.
-
The Distribution of Mass Surface Densities in a High-Mass Protocluster
Authors:
Wanggi Lim,
Jonathan C. Tan,
Jouni Kainulainen,
Bo Ma,
Michael J. Butler
Abstract:
We study the probability distribution function (PDF) of mass surface densities, $Σ$, of infrared dark cloud (IRDC) G028.37+00.07 and its surrounding giant molecular cloud. This PDF constrains the physical processes, such as turbulence, magnetic fields and self-gravity, that are expected to be controlling cloud structure and star formation activity. The chosen IRDC is of particular interest since i…
▽ More
We study the probability distribution function (PDF) of mass surface densities, $Σ$, of infrared dark cloud (IRDC) G028.37+00.07 and its surrounding giant molecular cloud. This PDF constrains the physical processes, such as turbulence, magnetic fields and self-gravity, that are expected to be controlling cloud structure and star formation activity. The chosen IRDC is of particular interest since it has almost 100,000 solar masses within a radius of 8 parsecs, making it one of the most massive, dense molecular structures known and is thus a potential site for the formation of a "super star cluster." We study $Σ$ in two ways. First, we use a combination of NIR and MIR extinction maps that are able to probe the bulk of the cloud structure up to $Σ\sim1\:{\rm{g\:cm}^{-2}}\:$($A_V\simeq200$ mag). Second, we study the FIR and sub-mm dust continuum emission from the cloud utilizing Herschel PACS and SPIRE images and paying careful attention to the effects of foreground and background contamination. We find that the PDFs from both methods, applied over a $\sim20^\prime$(30 pc)-sized region that contains $\simeq1.5\times10^5\:M_\odot$ and encloses a minimum closed contour with $Σ\simeq0.013\:{\rm{g\:cm}^{-2}}\:$($A_V\simeq3$ mag), shows a log-normal shape with the peak measured at $Σ\simeq0.021\:{\rm{g\:cm}^{-2}}\:$($A_V\simeq4.7$ mag). There is tentative evidence for the presence of a high-$Σ$ power law tail that contains from $\sim3\%$ to 8\% of the mass of the cloud material. We discuss the implications of these results for the physical processes occurring in this cloud.
△ Less
Submitted 7 September, 2016; v1 submitted 30 May, 2016;
originally announced May 2016.
-
An Ordered Bipolar Outflow from a Massive Early-Stage Core
Authors:
Jonathan C. Tan,
Shuo Kong,
Yichen Zhang,
Francesco Fontani,
Paola Caselli,
Michael J. Butler
Abstract:
We present ALMA follow-up observations of two massive, early-stage core candidates, C1-N & C1-S, in Infrared Dark Cloud (IRDC) G028.37+00.07, which were previously identified by their N2D+(3-2) emission and show high levels of deuteration of this species. The cores are also dark at far infrared wavelengths up to ~100 microns. We detect 12CO(2-1) from a narrow, highly-collimated bipolar outflow tha…
▽ More
We present ALMA follow-up observations of two massive, early-stage core candidates, C1-N & C1-S, in Infrared Dark Cloud (IRDC) G028.37+00.07, which were previously identified by their N2D+(3-2) emission and show high levels of deuteration of this species. The cores are also dark at far infrared wavelengths up to ~100 microns. We detect 12CO(2-1) from a narrow, highly-collimated bipolar outflow that is being launched from near the center of the C1-S core, which is also the location of the peak 1.3mm dust continuum emission. This protostar, C1-Sa, has associated dense gas traced by C18O(2-1) and DCN(3-2), from which we estimate it has a radial velocity that is near the center of the range exhibited by the C1-S massive core. A second outflow-driving source is also detected within the projected boundary of C1-S, but appears to be at a different radial velocity. After considering properties of the outflows, we conclude C1-Sa is a promising candidate for an early-stage massive protostar and as such it shows that these early phases of massive star formation can involve highly ordered outflow, and thus accretion, processes, similar to models developed to explain low-mass protostars.
△ Less
Submitted 2 March, 2016; v1 submitted 24 January, 2016;
originally announced January 2016.
-
Mid-J CO shock tracing observations of infrared dark clouds II Low-J CO constraints on excitation, depletion, and kinematics
Authors:
A. Pon,
D. Johnstone,
P. Caselli,
F. Fontani,
A. Palau,
M. J. Butler,
M. Kaufman,
I. Jiménez-Serra,
J. C. Tan
Abstract:
Infrared dark clouds are kinematically complex molecular structures in the interstellar medium that can host sites of massive star formation. We present 4 square arcminute maps of the 12CO, 13CO, and C18O J = 3 to 2 lines from selected locations within the C and F (G028.37+00.07 and G034.43+00.24) infrared dark clouds (IRDCs), as well as single pointing observations of the 13CO and C18O J = 2 to 1…
▽ More
Infrared dark clouds are kinematically complex molecular structures in the interstellar medium that can host sites of massive star formation. We present 4 square arcminute maps of the 12CO, 13CO, and C18O J = 3 to 2 lines from selected locations within the C and F (G028.37+00.07 and G034.43+00.24) infrared dark clouds (IRDCs), as well as single pointing observations of the 13CO and C18O J = 2 to 1 lines towards three cores within these clouds. We derive CO gas temperatures throughout the maps and find that CO is significantly frozen out within these IRDCs. We find that the CO depletion tends to be the highest near column density peaks, with maximum depletion factors between 5 and 9 in IRDC F and between 16 and 31 in IRDC C. We also detect multiple velocity components and complex kinematic structure in both IRDCs. Therefore, the kinematics of IRDCs seem to point to dynamically evolving structures yielding dense cores with considerable depletion factors.
△ Less
Submitted 15 December, 2015;
originally announced December 2015.
-
The Deuterium Fraction in Massive Starless Cores and Dynamical Implications
Authors:
Shuo Kong,
Jonathan C. Tan,
Paola Caselli,
Francesco Fontani,
Thushara Pillai,
Michael J. Butler,
Yoshito Shimajiri,
Fumitaka Nakamura,
Takeshi Sakai
Abstract:
We study deuterium fractionation in two massive starless/early-stage cores C1-N and C1-S in Infrared Dark Cloud (IRDC) G028.37+00.07, first identified by Tan et al. (2013) with ALMA. Line emission from multiple transitions of $\rm N_2H^+$ and $\rm N_2D^+$ were observed with the ALMA, CARMA, SMA, JCMT, NRO 45m and IRAM 30m telescopes. By simultaneously fitting the spectra, we estimate the excitatio…
▽ More
We study deuterium fractionation in two massive starless/early-stage cores C1-N and C1-S in Infrared Dark Cloud (IRDC) G028.37+00.07, first identified by Tan et al. (2013) with ALMA. Line emission from multiple transitions of $\rm N_2H^+$ and $\rm N_2D^+$ were observed with the ALMA, CARMA, SMA, JCMT, NRO 45m and IRAM 30m telescopes. By simultaneously fitting the spectra, we estimate the excitation conditions and deuterium fraction, $D_{\rm frac}^{\rm N_2H^+} \equiv [\rm N_2D^+]/[N_2H^+]$, with values of $D_{\rm frac}^{\rm N_2H^+} \simeq 0.2$--$0.7$, several orders of magnitude above the cosmic [D]/[H] ratio. Additional observations of o-H$_2$D$^+$ are also presented that help constrain the ortho-to-para ratio of $\rm H_2$, which is a key quantity affecting the degree of deuteration. We then present chemodynamical modeling of the two cores, exploring especially the implications for the collapse rate relative to free-fall, $α_{\rm ff}$. In order to reach the high level of observed deuteration of $\rm N_2H^+$, we find that the most likely evolutionary history of the cores involves collapse at a relatively slow rate, $\lesssim1/10$th of free-fall.
△ Less
Submitted 7 May, 2016; v1 submitted 29 September, 2015;
originally announced September 2015.
-
Mid-J CO Shock Tracing Observations of Infrared Dark Clouds I
Authors:
A. Pon,
P. Caselli,
D. Johnstone,
M. Kaufman,
M. J. Butler,
F. Fontani,
I. Jiménez-Serra,
J. C. Tan
Abstract:
Infrared dark clouds (IRDCs) are dense, molecular structures in the interstellar medium that can harbour sites of high-mass star formation. IRDCs contain supersonic turbulence, which is expected to generate shocks that locally heat pockets of gas within the clouds. We present observations of the CO J = 8-7, 9-8, and 10-9 transitions, taken with the Herschel Space Observatory, towards four dense, s…
▽ More
Infrared dark clouds (IRDCs) are dense, molecular structures in the interstellar medium that can harbour sites of high-mass star formation. IRDCs contain supersonic turbulence, which is expected to generate shocks that locally heat pockets of gas within the clouds. We present observations of the CO J = 8-7, 9-8, and 10-9 transitions, taken with the Herschel Space Observatory, towards four dense, starless clumps within IRDCs (C1 in G028.37+00.07, F1 and F2 in G034.43+0007, and G2 in G034.77-0.55). We detect the CO J = 8-7 and 9-8 transitions towards three of the clumps (C1, F1, and F2) at intensity levels greater than expected from photodissociation region (PDR) models. The average ratio of the 8-7 to 9-8 lines is also found to be between 1.6 and 2.6 in the three clumps with detections, significantly smaller than expected from PDR models. These low line ratios and large line intensities strongly suggest that the C1, F1, and F2 clumps contain a hot gas component not accounted for by standard PDR models. Such a hot gas component could be generated by turbulence dissipating in low velocity shocks.
△ Less
Submitted 7 May, 2015; v1 submitted 2 March, 2015;
originally announced March 2015.
-
Kiloparsec-Scale Simulations of Star Formation in Disk Galaxies III. Structure and Dynamics of Filaments and Clumps in Giant Molecular Clouds
Authors:
Michael J. Butler,
Jonathan C. Tan,
Sven Van Loo
Abstract:
We present hydrodynamic simulations of self-gravitating dense gas in a galactic disk, exploring scales ranging from 1 kpc down to $\sim 0.1$~pc. Our primary goal is to understand how dense filaments form in Giant Molecular Clouds (GMCs). These structures, often observed as Infrared Dark Clouds (IRDCs) in the Galactic plane, are thought to be the precursors to massive stars and star clusters, so th…
▽ More
We present hydrodynamic simulations of self-gravitating dense gas in a galactic disk, exploring scales ranging from 1 kpc down to $\sim 0.1$~pc. Our primary goal is to understand how dense filaments form in Giant Molecular Clouds (GMCs). These structures, often observed as Infrared Dark Clouds (IRDCs) in the Galactic plane, are thought to be the precursors to massive stars and star clusters, so their formation may be the rate limiting step controlling global star formation rates in galactic systems as described by the Kennicutt-Schmidt relation. Our study follows on from Van Loo et al. (2013, Paper I), which carried out simulations to 0.5~pc resolution and examined global aspects of the formation of dense gas clumps and the resulting star formation rate. Here, using our higher resolution, we examine the detailed structural, kinematic and dynamical properties of dense filaments and clumps, including mass surface density ($Σ$) probability distribution functions, filament mass per unit length and its dispersion, lateral $Σ$ profiles, filament fragmentation, filament velocity gradients and infall, and degree of filament and clump virialization. Where possible, these properties are compared to observations of IRDCs. By many metrics, especially too large mass fractions of high $Σ>1\:{\rm g\:cm^{-2}}$ material, too high mass per unit length dispersion due to dense clump formation, too high velocity gradients and too high velocity dispersion for a given mass per unit length, the simulated filaments differ from observed IRDCs. We thus conclude that IRDCs do not form from global fast collapse of GMCs. Rather, we expect IRDC formation and collapse is slowed significantly by the influence of dynamically important magnetic fields, which may thus play a crucial role in regulating galactic star formation rates.
△ Less
Submitted 16 April, 2015; v1 submitted 21 October, 2014;
originally announced October 2014.
-
The Darkest Shadows: Deep Mid-Infrared Extinction Mapping of a Massive Protocluster
Authors:
Michael J. Butler,
Jonathan C. Tan,
Jouni Kainulainen
Abstract:
We use deep $8\:μm$ Spitzer-IRAC imaging of a massive Infrared Dark Cloud (IRDC) G028.37+00.07 to construct a Mid-Infrared (MIR) extinction map that probes mass surface densities up to $Σ\:\sim 1\:\rm{g~cm^{-2}}$ ($A_V\sim200\:$mag), amongst the highest values yet probed by extinction mapping. Merging with a NIR extinction map of the region, creates a high dynamic range map that reveals structures…
▽ More
We use deep $8\:μm$ Spitzer-IRAC imaging of a massive Infrared Dark Cloud (IRDC) G028.37+00.07 to construct a Mid-Infrared (MIR) extinction map that probes mass surface densities up to $Σ\:\sim 1\:\rm{g~cm^{-2}}$ ($A_V\sim200\:$mag), amongst the highest values yet probed by extinction mapping. Merging with a NIR extinction map of the region, creates a high dynamic range map that reveals structures down to $A_V\sim1\:$mag. We utilize the map to: (1) Measure a cloud mass $\sim7\times10^4\:M_\odot$ within a radius of $\sim8\:$pc. $^{13}$CO kinematics indicate that the cloud is gravitationally bound. It thus has the potential to form one of the most massive young star clusters known in the Galaxy. (2) Characterize the structures of 16 massive cores within the IRDC, finding they can be fit by singular polytropic spheres with $ρ\propto{r}^{-k_ρ}$ and $k_ρ=1.3\pm0.3$. They have $\overlineΣ\simeq0.1-0.4\:\rm{g~cm^{-2}}$ --- relatively low values that, along with their measured cold temperatures, suggest magnetic fields, rather than accretion-powered radiative heating, are important for controlling fragmentation of these cores. (3) Determine the $Σ$ (equivalently column density or $A_V$) probability distribution function (PDF) for a region that is near complete for $A_V>3\:$mag. The PDF is well fit by a single log-normal with mean $\overline{A}_V\simeq9\:$mag, high compared to other known clouds. It does not exhibit a separate high-end power law tail, which has been claimed to indicate the importance of self-gravity. However, we suggest that the PDF does result from a self-similar, self-gravitating hierarchy of structure being present over a wide range of scales in the cloud.
△ Less
Submitted 27 January, 2014;
originally announced January 2014.
-
The Dynamics of Massive Starless Cores with ALMA
Authors:
Jonathan C. Tan,
Shuo Kong,
Michael J. Butler,
Paola Caselli,
Francesco Fontani
Abstract:
How do stars that are more massive than the Sun form, and thus how is the stellar initial mass function (IMF) established? Such intermediate- and high-mass stars may be born from relatively massive pre-stellar gas cores, which are more massive than the thermal Jeans mass. The Turbulent Core Accretion model invokes such cores as being in approximate virial equilibrium and in approximate pressure eq…
▽ More
How do stars that are more massive than the Sun form, and thus how is the stellar initial mass function (IMF) established? Such intermediate- and high-mass stars may be born from relatively massive pre-stellar gas cores, which are more massive than the thermal Jeans mass. The Turbulent Core Accretion model invokes such cores as being in approximate virial equilibrium and in approximate pressure equilibrium with their surrounding clump medium. Their internal pressure is provided by a combination of turbulence and magnetic fields. Alternatively, the Competitive Accretion model requires strongly sub-virial initial conditions that then lead to extensive fragmentation to the thermal Jeans scale, with intermediate- and high-mass stars later forming by competitive Bondi-Hoyle accretion. To test these models, we have identified four prime examples of massive (~100Msun) clumps from mid-infrared extinction mapping of infrared dark clouds (IRDCs). Fontani et al. found high deuteration fractions of N2H+ in these objects, which are consistent with them being starless. Here we present ALMA observations of these four clumps that probe the N2D+(3-2) line at 2.3" resolution. We find six N2D+ cores and determine their dynamical state. Their observed velocity dispersions and sizes are broadly consistent with the predictions of the Turbulent Core model of self-gravitating, magnetized (with Alfven Mach number m_A~1) and virialized cores that are bounded by the high pressures of their surrounding clumps. However, in the most massive cores, with masses up to ~60Msun, our results suggest that moderately enhanced magnetic fields (so that m_A~0.3) may be needed for the structures to be in virial and pressure equilibrium. Magnetically regulated core formation may thus be important in controlling the formation of massive cores, inhibiting their fragmentation, and thus helping to establish the stellar IMF.
△ Less
Submitted 19 November, 2013; v1 submitted 18 March, 2013;
originally announced March 2013.
-
Kiloparsec-Scale Simulations of Star Formation in Disk Galaxies. I. The unmagnetized and zero-feedback limit
Authors:
Sven Van Loo,
Michael J. Butler,
Jonathan C. Tan
Abstract:
We present hydrodynamic simulations of the evolution of self-gravitating dense gas on scales of 1 kiloparsec down to < parsec in a galactic disk, designed to study dense clump formation from giant molecular clouds (GMCs). These structures are expected to be the precursors to star clusters and this process may be the rate limiting step controling star formation rates in galactic systems as describe…
▽ More
We present hydrodynamic simulations of the evolution of self-gravitating dense gas on scales of 1 kiloparsec down to < parsec in a galactic disk, designed to study dense clump formation from giant molecular clouds (GMCs). These structures are expected to be the precursors to star clusters and this process may be the rate limiting step controling star formation rates in galactic systems as described by the Kennicutt-Schmidt relation. We follow the thermal evolution of the gas down to ~5K using extinction-dependent heating and cooling functions. We do not yet include magnetic fields or localized stellar feedback, so the evolution of the GMCs and clumps is determined solely by self-gravity balanced by thermal and turbulent pressure support and the large scale galactic shear. While cloud structures and densities change significantly during the simulation, GMC virial parameters remain mostly above unity for time scales exceeding the free-fall time of GMCs indicating that energy from galactic shear and large-scale cloud motions continuously cascades down to and within the GMCs. We implement star formation at a slow, inefficient rate of 2% per local free-fall time, but even this yields global star formation rates that are about two orders of magnitude larger than the observed Kennicutt-Schmidt relation due to over-production of dense gas clumps. We expect a combination of magnetic support and localized stellar feedback is required to inhibit dense clump formation to ~1% of the rate that results from the nonmagnetic, zero-feedback limit.
△ Less
Submitted 12 December, 2012; v1 submitted 7 November, 2012;
originally announced November 2012.
-
Kiloparsec-Scale Simulations of Magnetised Molecular Clouds in Disk Galaxies
Authors:
Sven Van Loo,
Michael J. Butler,
Jonathan C. Tan,
Sam A. E. G. Falle
Abstract:
We present simulations of the evolution of self-gravitating dense gas on kiloparsec-size scales in a galactic disk, designed to study dense clump formation from giant molecular clouds (GMCs). These dense clumps are expected to be the precursors to star clusters and this process may be the rate limiting step controling star formation rates in galactic systems as described by the Kennicutt-Schmidt r…
▽ More
We present simulations of the evolution of self-gravitating dense gas on kiloparsec-size scales in a galactic disk, designed to study dense clump formation from giant molecular clouds (GMCs). These dense clumps are expected to be the precursors to star clusters and this process may be the rate limiting step controling star formation rates in galactic systems as described by the Kennicutt-Schmidt relation. The evolution of these simulated GMCs and clumps is determined by self-gravity balanced by turbulent pressure support and the large scale galactic shear. While the cloud structures and densities significantly change during their evolution, they remain roughly in virial equilibrium for time scales exceeding the free-fall time of GMCs, indicating that energy from the galactic shear continuously cascades down. We implement star formation at a slow, inefficient rate of 2% per local free-fall time, but this yields global star formation rates that are more than ~two orders of magnitude larger than the observed Kennicutt-Schmidt relation due to the over-production of dense clump gas. To explain this discrepancy, we anticipate magnetic fields to provide additional support. Low-resolution simulations indeed show that the magnetic field reduces the star formation rate.
△ Less
Submitted 28 November, 2012; v1 submitted 2 November, 2012;
originally announced November 2012.
-
A Virialized Filamentary Infrared Dark Cloud
Authors:
Audra K. Hernandez,
Jonathan C. Tan,
Jouni Kainulainen,
Paola Caselli,
Michael J. Butler,
Izaskun Jimenez-Serra,
Francesco Fontani
Abstract:
The initial conditions of massive star and star cluster formation are expected to be cold, dense and high column density regions of the interstellar medium, which can reveal themselves via near, mid and even far-infrared absorption as Infrared Dark Clouds (IRDCs). Elucidating the dynamical state of IRDCs thus constrains theoretical models of these complex processes. In particular, it is important…
▽ More
The initial conditions of massive star and star cluster formation are expected to be cold, dense and high column density regions of the interstellar medium, which can reveal themselves via near, mid and even far-infrared absorption as Infrared Dark Clouds (IRDCs). Elucidating the dynamical state of IRDCs thus constrains theoretical models of these complex processes. In particular, it is important to assess whether IRDCs have reached virial equilibrium, where the internal pressure balances that due to the self-gravitating weight of the cloud plus the pressure of the external environmental. We study this question for the filamentary IRDC G035.39-00.33 by deriving mass from combined NIR & MIR extinction maps and velocity dispersion from C18O (1-0) & (2-1) line emission. In contrast to our previous moderately super-virial results based on 13CO emission and MIR-only extinction mapping, with improved mass measurements we now find that the filament is consistent with being in virial equilibrium, at least in its central parsec-wide region where ~1000 M_Sun snakes along several parsecs. This equilibrium state does not require large-scale net support or confinement by magnetic fields.
△ Less
Submitted 30 July, 2012; v1 submitted 17 July, 2012;
originally announced July 2012.
-
Mid-Infrared Extinction Mapping of Infrared Dark Clouds II. The Structure of Massive Starless Cores and Clumps
Authors:
Michael J. Butler,
Jonathan C. Tan
Abstract:
(abridged) We develop the mid-infrared extinction (MIREX) mapping technique of Butler & Tan (2009, Paper I), presenting a new method to correct for the Galactic foreground emission based on observed saturation in independent cores. Using Spitzer GLIMPSE 8 micron images, this allows us to accurately probe mass surface densities, Sigma, up to ~0.5g/cm^2 with 2" resolution. We then characterize the s…
▽ More
(abridged) We develop the mid-infrared extinction (MIREX) mapping technique of Butler & Tan (2009, Paper I), presenting a new method to correct for the Galactic foreground emission based on observed saturation in independent cores. Using Spitzer GLIMPSE 8 micron images, this allows us to accurately probe mass surface densities, Sigma, up to ~0.5g/cm^2 with 2" resolution. We then characterize the structure of 42 massive starless and early-stage IRDC cores and their surrounding clumps, measuring Sigma_cl(r) from the core/clump centers. We first assess the properties of the core/clump at a scale where the total enclosed mass as projected on the sky is M_cl=60Msun. We find these objects have a mean radius of R_cl~0.1pc, mean Sigma_cl=0.3g/cm^2 and, if fit by a power law density profile rho_cl ~ r^{-k_{rho,cl}}, a mean value of k_{rho,cl}=1.1. If we assume a core is embedded in each clump and subtract the surrounding clump envelope to derive the core properties, we find a mean core density power law index of k_{rho,c} = 1.6. We repeat this analysis as a function of radius and derive the best-fitting power law plus uniform clump envelope model for each of the 42 core/clumps. The cores have typical masses of M_c~100Msun and mean Sigma_c~0.1g/cm^2, and are embedded in clumps with comparable mass surface densities. We conclude massive starless cores exist and are well-described by singular polytropic spheres. Their relatively low values of Sigma and the fact that they are IR dark may imply that their fragmentation is inhibited by magnetic fields rather than radiative heating. Comparing to massive star-forming cores, there is tentative evidence for an evolution towards higher densities and steeper density profiles as star formation proceeds.
△ Less
Submitted 10 May, 2012;
originally announced May 2012.
-
Mapping Large-Scale CO Depletion in a Filamentary Infrared Dark Cloud
Authors:
Audra K. Hernandez,
Jonathan C. Tan,
Paola Caselli,
Michael J. Butler,
Izaskun Jimenez-Serra,
Francesco Fontani,
Peter Barnes
Abstract:
Infrared Dark Clouds (IRDCs) are cold, high mass surface density and high density structures, likely to be representative of the initial conditions for massive star and star cluster formation. CO emission from IRDCs has the potential to be useful for tracing their dynamics, but may be affected by depleted gas phase abundances due to freeze-out onto dust grains. Here we analyze C18O J=1-0 and J=2-1…
▽ More
Infrared Dark Clouds (IRDCs) are cold, high mass surface density and high density structures, likely to be representative of the initial conditions for massive star and star cluster formation. CO emission from IRDCs has the potential to be useful for tracing their dynamics, but may be affected by depleted gas phase abundances due to freeze-out onto dust grains. Here we analyze C18O J=1-0 and J=2-1 emission line data, taken with the IRAM 30m telescope, of the highly filamentary IRDC G035.39.-0033. We derive the excitation temperature as a function of position and velocity, with typical values of ~7K, and thus derive total mass surface densities, Sigma_C18O, assuming standard gas phase abundances and accounting for optical depth in the line, which can reach values of ~1. The mass surface densities reach values of ~0.07 g/cm^2. We compare these results to the mass surface densities derived from mid-infrared (MIR) extinction mapping, Sigma_SMF, by Butler & Tan, which are expected to be insensitive to the dust temperatures in the cloud. With a significance of >10sigma, we find Sigma_C18O/Sigma_SMF decreases by about a factor of 5 as Sigma increases from ~0.02 to ~0.2 g/cm^2, which we interpret as evidence for CO depletion. Several hundred solar masses are being affected, making this one of the most massive clouds in which CO depletion has been observed directly. We present a map of the depletion factor in the filament and discuss implications for the formation of the IRDC.
△ Less
Submitted 13 June, 2011;
originally announced June 2011.
-
Deuteration as an evolutionary tracer in massive-star formation
Authors:
F. Fontani,
Aina Palau,
P. Caselli,
A. Sanchez-Monge,
M. J. Butler,
J. C. Tan,
I. Jimenez-Serra,
G. Busquet,
S. Leurini,
M. Audard
Abstract:
Theory predicts, and observations confirm, that the column density ratio of a molecule containing D to its counterpart containing H can be used as an evolutionary tracer in the low-mass star formation process. Since it remains unclear if the high-mass star formation process is a scaled-up version of the low-mass one, we investigated whether the relation between deuteration and evolution can be app…
▽ More
Theory predicts, and observations confirm, that the column density ratio of a molecule containing D to its counterpart containing H can be used as an evolutionary tracer in the low-mass star formation process. Since it remains unclear if the high-mass star formation process is a scaled-up version of the low-mass one, we investigated whether the relation between deuteration and evolution can be applied to the high-mass regime. With the IRAM-30m telescope, we observed rotational transitions of N2D+ and N2H+ and derived the deuterated fraction in 27 cores within massive star-forming regions understood to represent different evolutionary stages of the massive-star formation process. Results. Our results clearly indicate that the abundance of N2D+ is higher at the pre-stellar/cluster stage, then drops during the formation of the protostellar object(s) as in the low-mass regime, remaining relatively constant during the ultra-compact HII region phase. The objects with the highest fractional abundance of N2D+ are starless cores with properties very similar to typical pre-stellar cores of lower mass. The abundance of N2D+ is lower in objects with higher gas temperatures as in the low-mass case but does not seem to depend on gas turbulence. Our results indicate that the N2D+-to-N2H+ column density ratio can be used as an evolutionary indicator in both low- and high-mass star formation, and that the physical conditions influencing the abundance of deuterated species likely evolve similarly during the processes that lead to the formation of both low- and high-mass stars.
△ Less
Submitted 29 March, 2011;
originally announced March 2011.
-
Parsec-scale SiO Emission in an Infrared Dark Cloud
Authors:
I. Jimenez-Serra,
P. Caselli,
J. C. Tan,
A. K. Hernandez,
F. Fontani,
M. J. Butler,
S. van Loo
Abstract:
We present high-sensitivity 2'x4' maps of the J=2-1 rotational lines of SiO, CO, 13CO and C18O, observed toward the filamentary Infrared Dark Cloud (IRDC) G035.39-00.33. Single-pointing spectra of the SiO J=2-1 and J=3-2 lines toward several regions in the filament, are also reported. The SiO images reveal that SiO is widespread along the IRDC (size >2 pc), showing two different components: one br…
▽ More
We present high-sensitivity 2'x4' maps of the J=2-1 rotational lines of SiO, CO, 13CO and C18O, observed toward the filamentary Infrared Dark Cloud (IRDC) G035.39-00.33. Single-pointing spectra of the SiO J=2-1 and J=3-2 lines toward several regions in the filament, are also reported. The SiO images reveal that SiO is widespread along the IRDC (size >2 pc), showing two different components: one bright and compact arising from three condensations (N, E and S), and the other weak and extended along the filament. While the first component shows broad lines (linewidths of ~4-7 kms-1) in both SiO J=2-1 and SiO J=3-2, the second one is only detected in SiO J=2-1 and has narrow lines (~0.8 kms-1). The maps of CO and its isotopologues show that low-density filaments are intersecting the IRDC and appear to merge toward the densest portion of the cloud. This resembles the molecular structures predicted by flow-driven, shock-induced and magnetically-regulated cloud formation models. As in outflows associated with low-mass star formation, the excitation temperatures and fractional abundances of SiO toward N, E and S, increase with velocity from ~6 to 40 K, and from ~1E-10 to >1E-8 respectively, over a velocity range of ~7 kms-1. Since 8 micron sources, 24 micron sources and/or extended 4.5 micron emission are detected in N, E and S, broad SiO is likely produced in outflows associated with high-mass protostars. The excitation temperatures and fractional abundances of the narrow SiO lines, however, are very low (~9 K and ~1E-11, respectively), and consistent with the processing of interstellar grains by the passage of a shock with vs~12 kms-1. This emission could be generated i) by a large-scale shock, perhaps remnant of the IRDC formation process; ii) by decelerated or recently processed gas in large-scale outflows driven by 8 micron and 24 micron sources; or iii) by an undetected and widespread population of lower mass protostars. High-angular resolution observations are needed to disentangle between these three scenarios.
△ Less
Submitted 17 March, 2010;
originally announced March 2010.
-
Mid-Infrared Extinction Mapping of Infrared Dark Clouds: Probing the Initial Conditions for Massive Stars and Star Clusters
Authors:
Michael J. Butler,
Jonathan C. Tan
Abstract:
(Abridged) We use 8 micron Spitzer GLIMPSE images to make extinction maps of 10 IRDCs, selected to be relatively nearby and massive. The extinction mapping technique requires modeling the IR background intensity behind the cloud, which is achieved by correcting for foreground emission and then interpolating from the surrounding regions. The correction for foreground emission can be quite large,…
▽ More
(Abridged) We use 8 micron Spitzer GLIMPSE images to make extinction maps of 10 IRDCs, selected to be relatively nearby and massive. The extinction mapping technique requires modeling the IR background intensity behind the cloud, which is achieved by correcting for foreground emission and then interpolating from the surrounding regions. The correction for foreground emission can be quite large, thus restricting the utility of this technique to relatively nearby clouds. We investigate three methods for the interpolation, finding systematic differences at about the 10% level, which, for fiducial dust models, corresponds to a mass surface density Sigma = 0.013 g cm^-2, above which we conclude this extinction mapping technique attains validity. We examine the probability distribution function of Sigma in IRDCs. From a qualitative comparison with numerical simulations of astrophysical turbulence, many clouds appear to have relatively narrow distributions suggesting relatively low (<5) Mach numbers and/or dynamically strong magnetic fields. Given cloud kinematic distances, we derive cloud masses. Rathborne, Jackson & Simon identified cores within the clouds and measured their masses via mm dust emission. For 43 cores, we compare these mass estimates with those derived from our extinction mapping, finding good agreement: typically factors of <~2 difference for individual cores and an average systematic offset of <~10% for the adopted fiducial assumptions of each method. We find tentative evidence for a systematic variation of these mass ratios as a function of core density, which is consistent with models of ice mantle formation on dust grains and subsequent grain growth by coagulation, and/or with a temperature decrease in the densest cores.
△ Less
Submitted 14 April, 2009; v1 submitted 15 December, 2008;
originally announced December 2008.