Source author record

Michael J. Butler

Michael J. Butler appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2020arXiv

The AGORA high-resolution galaxy simulations comparison project: Public data release

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.

preprint2016arXiv

An Ordered Bipolar Outflow from a Massive Early-Stage Core

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.

preprint2016arXiv

The Deuterium Fraction in Massive Starless Cores and Dynamical Implications

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.

preprint2016arXiv

The Distribution of Mass Surface Densities in a High-Mass Protocluster

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.

preprint2015arXiv

Kiloparsec-Scale Simulations of Star Formation in Disk Galaxies III. Structure and Dynamics of Filaments and Clumps in Giant Molecular Clouds

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.

preprint2014arXiv

The Darkest Shadows: Deep Mid-Infrared Extinction Mapping of a Massive Protocluster

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.

preprint2013arXiv

The Dynamics of Massive Starless Cores with ALMA

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.

preprint2012arXiv

A Virialized Filamentary Infrared Dark Cloud

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.

preprint2012arXiv

Kiloparsec-Scale Simulations of Magnetised Molecular Clouds in Disk Galaxies

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.

preprint2012arXiv

Kiloparsec-Scale Simulations of Star Formation in Disk Galaxies. I. The unmagnetized and zero-feedback limit

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.

preprint2012arXiv

Mid-Infrared Extinction Mapping of Infrared Dark Clouds II. The Structure of Massive Starless Cores and Clumps

(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.

preprint2011arXiv

Mapping Large-Scale CO Depletion in a Filamentary Infrared Dark Cloud

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.

preprint2009arXiv

Mid-Infrared Extinction Mapping of Infrared Dark Clouds: Probing the Initial Conditions for Massive Stars and Star Clusters

(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.