BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Christopher F. McKee

Christopher F. McKee appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAastro-ph.SRastro-ph.COastro-phastro-ph.EPastro-ph.IMphysics.comp-phphysics.flu-dyn

Catalog footprint

What is connected

47works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

47 published item(s)

preprint2022arXiv

Magnetic fields in the formation of the first stars.--II Results

Beginning with cosmological initial conditions at z=100, we simulate the effects of magnetic fields on the formation of Population III stars and compare our results with the predictions of Paper I. We use Gadget-2 to follow the evolution of the system while the field is weak. We introduce a new method for treating kinematic fields by tracking the evolution of the deformation tensor. The growth rate in this stage of the simulation is lower than expected for diffuse astrophysical plasmas, which have a very low resistivity (high magnetic Prandtl number); we attribute this to the large numerical resistivity in simulations, corresponding to a magnetic Prandtl number of order unity. When the magnetic field begins to be dynamically significant in the core of the minihalo at z=27, we map it onto a uniform grid and follow the evolution in an adaptive mesh refinement, MHD simulation in Orion2. The nonlinear evolution of the field in the Orion2 simulation violates flux-freezing and is consistent with the theory proposed by Xu & Lazarian. The fields approach equipartition with kinetic energy at densities ~ 10^10 - 10^12 cm^-3. When the same calculation is carried out in Orion2 with no magnetic fields, several protostars form, ranging in mass from ~ 1 to 30 M_sol with magnetic fields, only a single ~ 30 M_sol protostar forms by the end of the simulation. Magnetic fields thus suppress the formation of low-mass Pop III stars, yielding a top-heavy Pop III IMF and contributing to the absence of observed Pop III stars.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

How do bound star clusters form?

Gravitationally-bound clusters that survive gas removal represent an unusual mode of star formation in the Milky Way and similar spiral galaxies. While forming, they can be distinguished observationally from unbound star formation by their high densities, virialised velocity structures, and star formation histories that accelerate toward the present, but extend multiple free-fall times into the past. In this paper we examine several proposed scenarios for how such structures might form and evolve, and carry out a Bayesian analysis to test these models against observed distributions of protostellar age, counts of young stellar objects relative to gas, and the overall star formation rate of the Milky Way. We show that models in which the acceleration of star formation is due either to a large-scale collapse or a time-dependent increase in star formation efficiency are unable to satisfy the combined set of observational constraints. In contrast, models in which clusters form in a "conveyor belt" mode where gas accretion and star formation occur simultaneously, but the star formation rate per free-fall time is low, can match the observations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Magnetic Fields in the Formation of the First Stars. I. Theory vs. Simulation

While magnetic fields are important in contemporary star formation, their role in primordial star formation is unknown. Magnetic fields of order 10^-16 G are produced by the Biermann battery due to the curved shocks and turbulence associated with the infall of gas into the dark matter minihalos that are the sites of formation of the first stars. These fields are rapidly amplified by a small-scale dynamo until they saturate at or near equipartition with the turbulence in the central region of the gas. Analytic results are given for the outcome of the dynamo, including the effect of compression in the collapsing gas. The mass-to-flux ratio in this gas is 2-3 times the critical value, comparable to that in contemporary star formation. Predictions of the outcomes of simulations using smooth particle hydrodynamics (SPH) and grid-based adaptive mesh refinement (AMR) are given. Because the numerical viscosity and resistivity for the standard resolution of 64 cells per Jeans length are several orders of magnitude greater than the physical values, dynamically significant magnetic fields affect a much smaller fraction of the mass in simulations than in reality. An appendix gives an analytic treatment of free-fall collapse, including that in a constant density background. Another appendix presents a new method of estimating the numerical viscosity; results are given for both SPH and grid-based codes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Massive Warm/Hot Galaxy Coronae: II. Isentropic Model

We construct a new analytic phenomenological model for the extended circumgalactic material (CGM) of $L^*$ galaxies. Our model reproduces the OVII/OVIII absorption observations of the Milky Way (MW) and the OVI measurements reported by the COS-Halos and eCGM surveys. The warm/hot gas is in hydrostatic equilibrium in a MW gravitational potential, and we adopt a barotropic equation of state, resulting in a temperature variation as a function of radius. A pressure component with an adiabatic index of $γ=4/3$ is included to approximate the effects of a magnetic field and cosmic rays. We introduce a metallicity gradient motivated by the enrichment of the inner CGM by the Galaxy. We then present our fiducial model for the corona, tuned to reproduce the observed OVI-OVIII column densities, and with a total mass of $M_{\rm gas} \approx 5.5 \times 10^{10}~{\rm M_{\odot}}$ inside $r_{\rm cgm} \approx 280$ kpc. The gas densities in the CGM are low ($n_{\rm H} = 10^{-5} - 3 \times 10^{-4}~{\rm cm^{-3}}$) and its collisional ionization state is modified by the metagalactic radiation field (MGRF). We show that for OVI-bearing warm/hot gas with typical observed column densities $N_{\rm OVI} \sim 3 \times 10^{14}~{\rm cm^{-2}}$ at large ($\gtrsim 100$ kpc) impact parameters from the central galaxies, the ratio of the cooling to dynamical times, $t_{\rm cool}/t_{\rm dyn}$, has a model-independent upper limit of $\lesssim 4$. In our model, $t_{\rm cool}/t_{\rm dyn}$ at large radii is $\sim 2-3$. We present predictions for a wide range of future observations of the warm/hot CGM, from UV/X-ray absorption and emission spectroscopy, to dispersion measure (DM) and Sunyaev-Zeldovich CMB measurements. We provide the model outputs in machine-readable data files, for easy comparison and analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

An Unstable Truth: How Massive Stars get their Mass

The pressure exerted by massive stars' radiation fields is an important mechanism regulating their formation. Detailed simulation of massive star formation therefore requires an accurate treatment of radiation. However, all published simulations have either used a diffusion approximation of limited validity; have only been able to simulate a single star fixed in space, thereby suppressing potentially-important instabilities; or did not provide adequate resolution at locations where instabilities may develop. To remedy this we have developed a new, highly accurate radiation algorithm that properly treats the absorption of the direct radiation field from stars and the re-emission and processing by interstellar dust. We use our new tool to perform three-dimensional radiation-hydrodynamic simulations of the collapse of massive pre-stellar cores with laminar and turbulent initial conditions and properly resolve regions where we expect instabilities to grow. We find that mass is channeled to the stellar system via gravitational and Rayleigh-Taylor (RT) instabilities, in agreement with previous results using stars capable of moving, but in disagreement with methods where the star is held fixed or with simulations that do not adequately resolve the development of RT instabilities. For laminar initial conditions, proper treatment of the direct radiation field produces later onset of instability, but does not suppress it entirely provided the edges of radiation-dominated bubbles are adequately resolved. Instabilities arise immediately for turbulent pre-stellar cores because the initial turbulence seeds the instabilities. Our results suggest that RT features are significant and should be present around accreting massive stars throughout their formation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

What Physics Determines the Peak of the IMF? Insights from the Structure of Cores in Radiation-Magnetohydrodynamic Simulations

As star-forming clouds collapse, the gas within them fragments to ever-smaller masses. Naively one might expect this process to continue down to the smallest mass that is able to radiate away its binding energy on a dynamical timescale, the opacity limit for fragmentation, at $\sim 0.01$ $M_\odot$. However, the observed peak of the initial mass function (IMF) lies a factor of $20-30$ higher in mass, suggesting that some other mechanism halts fragmentation before the opacity limit is reached. In this paper we analyse radiation-magnetohydrodynamic simulations of star cluster formation in typical Milky Way environments in order to determine what physical process limits fragmentation in them. We examine the regions in the vicinity of stars that form in the simulations to determine the amounts of mass that are prevented from fragmenting by thermal and magnetic pressure. We show that, on small scales, thermal pressure enhanced by stellar radiation heating is the dominant mechanism limiting the ability of the gas to further fragment. In the brown dwarf mass regime, $\sim 0.01$ $M_\odot$, the typical object that forms in the simulations is surrounded by gas whose mass is several times its own that is unable to escape or fragment, and instead is likely to accrete. This mechanism explains why $\sim 0.01$ $M_\odot$ objects are rare: unless an outside agent intervenes (e.g., a shock strips away the gas around them), they will grow by accreting the warmed gas around them. In contrast, by the time stars grow to masses of $\sim 0.2$ $M_\odot$, the mass of heated gas is only tens of percent of the central star mass, too small to alter its final mass by a large factor. This naturally explains why the IMF peak is at $\sim 0.2$ $M_\odot$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Magnetized Interstellar Molecular Clouds. I. Comparison Between Simulations and Zeeman Observations

The most accurate measurements of magnetic fields in star-forming gas are based on the Zeeman observations analyzed by Crutcher et al. (2010). We show that their finding that the 3D magnetic field scales approximately as density$^{0.65}$ can also be obtained from analysis of the observed line-of-sight fields. We present two large-scale AMR MHD simulations of several thousand $M_\odot$ of turbulent, isothermal, self-gravitating gas, one with a strong initial magnetic field (Alfven Mach number $M_{A,0}= 1$) and one with a weak initial field ($M_{A,0}=10$). We construct samples of the 100 most massive clumps in each simulation and show that they exhibit a power-law relation between field strength and density in excellent agreement with the observed one. Our results imply that the average field in molecular clumps in the interstellar medium is $<B_{tot}> \sim 42 n_{H,4}^{0.65} μ$G. Furthermore, the median value of the ratio of the line-of-sight field to density$^{0.65}$ in the simulations is within a factor of about (1.3, 1.7) of the observed value for the strong and weak field cases, respectively. The median value of the mass-to-flux ratio, normalized to the critical value, is 70% of the line-of-sight value. This is larger than the 50% usually cited for spherical clouds because the actual mass-to-flux ratio depends on the volume-weighted field, whereas the observed one depends on the mass-weighted field. Our results indicate that the typical molecular clump in the ISM is significantly supercritical (~ factor of 3). The results of our strong-field model are in very good quantitative agreement with the observations of Li et al. (2009), which show a strong correlation in field orientation between small and large scales. Because there is a negligible correlation in the weak-field model, we conclude that molecular clouds form from strongly magnetized (although magnetically supercritical) gas.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stars, Gas, and Dark Matter in the Solar Neighborhood

The surface density and vertical distribution of stars, stellar remnants, and gas in the solar vicinity form important ingredients for understanding the star formation history of the Galaxy as well as for inferring the local density of dark matter by using stellar kinematics to probe the gravitational potential. In this paper we review the literature for these baryonic components, reanalyze data, and provide tables of the surface densities and exponential scale heights of main sequence stars, giants, brown dwarfs, and stellar remnants. We also review three components of gas (H2, HI, and HII), give their surface densities at the solar circle, and discuss their vertical distribution. We find a local total surface density of M dwarfs of 17.3 pm 2.3 Mo/pc^2. Our result for the total local surface density of visible stars, 27.0 pm 2.7 Mo/pc^2, is close to previous estimates due to a cancellation of opposing effects: more mass in M dwarfs, less mass in the others. The total local surface density in white dwarfs is 4.9 pm 0.6 Mo/pc^2; in brown dwarfs, it is ~1.2 Mo/pc^2. We find that the total local surface density of stars and stellar remnants is 33.4 pm 3 Mo/pc^2, somewhat less than previous estimates. We analyze data on 21 cm emission and absorption and obtain good agreement with recent results on the local amount of neutral atomic hydrogen obtained with the Planck satellite. The local surface density of gas is 13.7 pm 1.6 Mo/pc^2. The total baryonic mass surface density that we derive for the solar neighborhood is 47.1 pm 3.4 Mo/pc^2. Combining these results with others' measurements of the total surface density of matter within 1-1.1 kpc of the plane, we find that the local density of dark matter is 0.013 pm 0.003Mo/pc^3.The local density of all matter is 0.097 pm 0.013 Mo/pc^3. We discuss limitations on the properties of a possible thin disk of dark matter.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Turbulent Origin of Spin-Orbit Misalignment in Planetary Systems

The turbulent environment from which stars form may lead to misalignment between the stellar spin and the remnant protoplanetary disk. By using hydrodynamic and magnetohydrodynamic simulations, we demonstrate that a wide range of stellar obliquities may be produced as a by-product of forming a star within a turbulent environment. We present a simple semi-analytic model that reveals this connection between the turbulent motions and the orientation of a star and its disk. Our results are consistent with the observed obliquity distribution of hot Jupiters. Migration of misaligned hot Jupiters may, therefore, be due to tidal dissipation in the disk, rather than tidal dissipation of the star-planet interaction.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Bondi-Hoyle Accretion in an Isothermal Magnetized Plasma

In regions of star formation, protostars and newborn stars accrete mass from their natal clouds. These clouds are threaded by magnetic fields with a strength characterized by the plasma beta---the ratio of thermal and magnetic pressures. Observations show molecular clouds have beta <= 1, so magnetic fields can play a significant role in the accretion process. We have carried out a numerical study of the effect of large-scale magnetic fields on the rate of accretion onto a uniformly moving point particle from a uniform, non-self-gravitating, isothermal gas. We consider gas moving with sonic Mach numbers of up M ~ 45, magnetic fields that are either parallel, perpendicular, or oriented 45 degrees to the flow, and beta as low as 0.01. Our simulations utilize AMR to obtain high spatial resolution where needed; this also allows the simulation boundaries to be far from the accreting object. Additionally, we show our results are independent of our exact prescription for accreting mass in the sink particle. We give simple expressions for the steady-state accretion rate as a function of beta, M, and field orientation. Using typical molecular clouds values of M ~ 5 and beta ~ 0.04 from the literature, our fits suggest a 0.4 M_Sun star accretes ~ 4*10^{-9} M_Sun/year, almost a factor of two less than accretion rates predicted by hydrodynamic models. This disparity grows to orders of magnitude for stronger fields and lower Mach numbers. We discuss the applicability of these accretion rates versus accretion rates expected from gravitational collapse, and when a steady state is possible. This reduction in Mdot increases the time required to form stars in competitive accretion models, making such models less efficient. In numerical codes, our results should enable accurate subgrid models of sink particles accreting from magnetized media.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Massive Star Formation

The enormous radiative and mechanical luminosities of massive stars impact a vast range of scales and processes, from the reionization of the universe, to the evolution of galaxies, to the regulation of the interstellar medium, to the formation of star clusters, and even to the formation of planets around stars in such clusters. Two main classes of massive star formation theory are under active study, Core Accretion and Competitive Accretion. In Core Accretion, the initial conditions are self-gravitating, centrally concentrated cores that condense with a range of masses from the surrounding, fragmenting clump environment. They then undergo relatively ordered collapse via a central disk to form a single star or a small-N multiple. In this case, the pre-stellar core mass function has a similar form to the stellar initial mass function. In Competitive Accretion, the material that forms a massive star is drawn more chaotically from a wider region of the clump without passing through a phase of being in a massive, coherent core. In this case, massive star formation must proceed hand in hand with star cluster formation. If stellar densities become very high near the cluster center, then collisions between stars may also help to form the most massive stars. We review recent theoretical and observational progress towards understanding massive star formation, considering physical and chemical processes, comparisons with low and intermediate-mass stars, and connections to star cluster formation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Massive Protostar Forming by Ordered Collapse of a Dense, Massive Core

We present 30 and 40 micron imaging of the massive protostar G35.20-0.74 with SOFIA-FORCAST. The high surface density of the natal core around the protostar leads to high extinction, even at these relatively long wavelengths, causing the observed flux to be dominated by that emerging from the near-facing outflow cavity. However, emission from the far-facing cavity is still clearly detected. We combine these results with fluxes from the near-infrared to mm to construct a spectral energy distribution (SED). For isotropic emission the bolometric luminosity would be 3.3x10^4 Lsun. We perform radiative transfer modeling of a protostar forming by ordered, symmetric collapse from a massive core bounded by a clump with high mass surface density, Sigma_cl. To fit the SED requires protostellar masses ~20-34 Msun depending on the outflow cavity opening angle (35 - 50 degrees), and Sigma_cl ~ 0.4-1 g cm-2. After accounting for the foreground extinction and the flashlight effect, the true bolometric luminosity is ~ (0.7-2.2)x10^5 Lsun. One of these models also has excellent agreement with the observed intensity profiles along the outflow axis at 10, 18, 31 and 37 microns. Overall our results support a model of massive star formation involving the relatively ordered, symmetric collapse of a massive, dense core and the launching bipolar outflows that clear low density cavities. Thus a unified model may apply for the formation of both low and high mass stars.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Interstellar H$_2$O Masers from J Shocks

We present a model in which the 22 GHz H$_2$O masers observed in star-forming regions occur behind shocks propagating in dense regions (preshock density $n_0 \sim 10^6 - 10^8$ cm$^{-3}$). We focus on high-velocity ($v_s > 30$ km s$^{-1}$) dissociative J shocks in which the heat of H$_2$ re-formation maintains a large column of $\sim 300-400$ K gas; at these temperatures the chemistry drives a considerable fraction of the oxygen not in CO to form H$_2$O. The H$_2$O column densities, the hydrogen densities, and the warm temperatures produced by these shocks are sufficiently high to enable powerful maser action. The observed brightness temperatures (generally $\sim 10^{11} - 10^{14}$ K) are the result of coherent velocity regions that have dimensions in the shock plane that are 10 to 100 times the shock thickness of $\sim 10^{13}$ cm. The masers are therefore beamed towards the observer, who typically views the shock "edge-on", or perpendicular to the shock velocity; the brightest masers are then observed with the lowest line of sight velocities with respect to the ambient gas. We present numerical and analytic studies of the dependence of the maser inversion, the resultant brightness temperature, the maser spot size and shape, the isotropic luminosity, and the maser region magnetic field on the shock parameters and the coherence path length; the overall result is that in galactic H$_2$O 22 GHz masers these observed parameters can be produced in J shocks with $n_0\sim 10^6 - 10^8$ cm$^{-3}$ and $v_s \sim 30 -200$ km s$^{-1}$. A number of key observables such as maser shape, brightness temperature, and global isotropic luminosity depend only on the particle flux into the shock, $j=n_0v_s$, rather than on $n_0$ and $v_s$ separately.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Radiation Transfer of Models of Massive Star Formation. II. Effects of the Outflow

(Abridged) We present radiation transfer simulations of a massive (8 Msun) protostar forming from a massive (Mc=60 Msun) protostellar core, extending the model developed by Zhang & Tan (2011). The two principal improvements are (1) developing a model for the density and velocity structure of a disk wind that fills the bipolar outflow cavities; and (2) solving for the radially varying accretion rate in the disk due to a supply of mass and angular momentum from the infall envelope and their loss to the disk wind. One consequence of the launching of the disk wind is a reduction in the amount of accretion power that is radiated by the disk. For the transition from dusty to dust-free conditions where gas opacities dominate, we now implement a gradual change as a more realistic approximation of dust destruction. We study how the above effects, especially the outflow, influence the SEDs and the images of the protostar. Dust in the outflow cavity significantly affects the SEDs at most viewing angles. It further attenuates the short-wavelength flux from the protostar, controlling how the accretion disk may be viewed, and contributes a significant part of the near- and mid-IR fluxes. These fluxes warm the disk, boosting the mid- and far-IR emission. We find that for near face-on views, the SED from the near-IR to about 60 micron is very flat, which may be used to identify such systems. We show that the near-facing outflow cavity and its walls are still the most significant features in images up to 70 micron, dominating the mid-IR emission and determining its morphology. The thermal emission from the dusty outflow itself dominates the flux at ~20 micron. The detailed distribution of the dust in the outflow affects the morphology, for example, even though the outflow cavity is wide, at 10 to 20 micron, the dust in the disk wind can make the outflow appear narrower than in the near-IR bands.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Fragmentation of Magnetized, Massive Star-Forming Cores with Radiative Feedback

We present a set of 3-dimensional, radiation-magnetohydrodynamic calculations of the gravitational collapse of massive (300 Msun), star-forming molecular cloud cores. We show that the combined effects of magnetic fields and radiative feedback strongly suppress core fragmentation, leading to the production of single star systems rather than small clusters. We find that the two processes are efficient at suppressing fragmentation in different regimes, with the feedback most effective in the dense, central region and the magnetic field most effective in more diffuse, outer regions. Thus, the combination of the two is much more effective at suppressing fragmentation than either one considered in isolation. Our work suggests that typical massive cores, which have mass-to-flux ratios of about 2 relative to critical, likely form a single star system, but that cores with weaker fields may form a small star cluster. This result helps us understand why the observed relationship between the core mass function and the stellar initial mass function holds even for ~100 Msun cores with many thermal Jeans masses of material. We also demonstrate that a ~40 AU Keplerian disk is able to form in our simulations, despite the braking effect caused by the strong magnetic field.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Star Formation Rate of Molecular Clouds

We review recent advances in the analytical and numerical modeling of the star formation rate in molecular clouds and discuss the available observational constraints. We focus on molecular clouds as the fundamental star formation sites, rather than on the larger-scale processes that form the clouds and set their properties. Molecular clouds are shaped into a complex filamentary structure by supersonic turbulence, with only a small fraction of the cloud mass channeled into collapsing protostars over a free-fall time of the system. In recent years, the physics of supersonic turbulence has been widely explored with computer simulations, leading to statistical models of this fragmentation process, and to the prediction of the star formation rate as a function of fundamental physical parameters of molecular clouds, such as the virial parameter, the rms Mach number, the compressive fraction of the turbulence driver, and the ratio of gas to magnetic pressure. Infrared space telescopes, as well as ground-based observatories have provided unprecedented probes of the filamentary structure of molecular clouds and the location of forming stars within them.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ambipolar Diffusion Heating in Turbulent Systems

The temperature of the gas in molecular clouds is a key determinant of the characteristic mass of star formation. Ambipolar diffusion (AD) is considered one of the most important heating mechanisms in weakly ionized molecular clouds. In this work, we study the AD heating rate using 2-fluid turbulence simulations and compare it with the overall heating rate due to turbulent dissipation. We find that for observed molecular clouds, which typically have Alfven Mach numbers of ~1 (Crutcher 1999) and AD Reynolds numbers of ~20 (McKee et al. 2010), about 70% of the total turbulent dissipation is in the form of AD heating. AD has an important effect on the length scale where energy is dissipated: when AD heating is strong, most of the energy in the cascade is removed by ion-neutral drift, with a comparatively small amount of energy making it down to small scales. We derive a relation for the AD heating rate that describes the results of our simulations to within a factor of two. Turbulent dissipation, including AD heating, is generally less important that cosmic-ray heating in molecular clouds, although there is substantial scatter in both.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Feedback Effects on Low-Mass Star Formation

Protostellar feedback, both radiation and bipolar outflows, dramatically affects the fragmentation and mass accretion from star-forming cores. We use ORION, an adaptive mesh refinement gravito-radiation-hydrodynamics code, to simulate the formation of a cluster of low-mass stars, including both radiative transfer and protostellar outflows. We ran four simulations to isolate the individual effects of radiation feedback and outflow feedback as well as the combination of the two. Outflows reduce protostellar masses and accretion rates each by a factor of three and therefore reduce protostellar luminosities by an order of magnitude. Thus, while radiation feedback suppresses fragmentation, outflows render protostellar radiation largely irrelevant for low-mass star formation above a mass scale of 0.05 M_sun. We find initial fragmentation of our cloud at half the global Jeans length, ~ 0.1 pc. With insufficient protostellar radiation to stop it, these 0.1 pc cores fragment repeatedly, forming typically 10 stars each. The accretion rate in these stars scales with mass as predicted from core accretion models that include both thermal and turbulent motions. We find that protostellar outflows do not significantly affect the overall cloud dynamics, in the absence of magnetic fields, due to their small opening angles and poor coupling to the dense gas. The outflows reduce the mass from the cores by 2/3, giving a core to star efficiency ~ 1/3. The simulation with radiation and outflows reproduces the observed protostellar luminosity function. All of the simulations can reproduce observed core mass functions, though they are sensitive to telescope resolution. The simulation with both radiation and outflows reproduces the galactic IMF and the two-point correlation function of the cores observed in rho Oph.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Observing Simulated Protostars with Outflows: How Accurate are Protostellar Properties Inferred from SEDs?

The properties of unresolved protostars and their local environment are frequently inferred from spectral energy distributions (SEDs) using radiative transfer modeling. We perform synthetic observations of realistic star formation simulations to evaluate the accuracy of properties inferred from fitting model SEDs to observations. We use ORION, an adaptive mesh refinement (AMR) three-dimensional gravito-radiation-hydrodynamics code, to simulate low-mass star formation in a turbulent molecular cloud including the effects of protostellar outflows. To obtain the dust temperature distribution and SEDs of the forming protostars, we post-process the simulations using HYPERION, a state-of-the-art Monte-Carlo radiative transfer code. We find that the ORION and HYPERION dust temperatures typically agree within a factor of two. We compare synthetic SEDs of embedded protostars for a range of evolutionary times, simulation resolutions, aperture sizes, and viewing angles. We demonstrate that complex, asymmetric gas morphology leads to a variety of classifications for individual objects as a function of viewing angle. We derive best-fit source parameters for each SED through comparison with a pre-computed grid of radiative transfer models. While the SED models correctly identify the evolutionary stage of the synthetic sources as embedded protostars, we show that the disk and stellar parameters can be very discrepant from the simulated values. Parameters such as the stellar accretion rate, stellar mass, and disk mass show better agreement, but can still deviate significantly, and the agreement may in some cases be artificially good due to the limited range of parameters in the set of model SEDs. Lack of correlation between the model and simulation properties in many individual instances cautions against over-interpreting properties inferred from SEDs for unresolved protostellar sources. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Photometric Redshifts of Submillimeter Galaxies

We use the photometric redshift method of Chakrabarti & McKee (2008) to infer photometric redshifts of submillimeter galaxies with far-IR (FIR) $\it{Herschel}$ data obtained as part of the PACS Evolutionary Probe (PEP) program. For the sample with spectroscopic redshifts, we demonstrate the validity of this method over a large range of redshifts ($ 4 \ga z \ga 0.3$) and luminosities, finding an average accuracy in $(1+z_{\rm phot})/(1+z_{\rm spec})$ of 10%. Thus, this method is more accurate than other FIR photometric redshift methods. This method is different from typical FIR photometric methods in deriving redshifts from the light-to-gas mass ($L/M$) ratio of infrared-bright galaxies inferred from the FIR spectral energy distribution (SED), rather than dust temperatures. Once the redshift is derived, we can determine physical properties of infrared bright galaxies, including the temperature variation within the dust envelope, luminosity, mass, and surface density. We use data from the GOODS-S field to calculate the star formation rate density (SFRD) of sub-mm bright sources detected by AzTEC and PACS. The AzTEC-PACS sources, which have a threshold $850 \micron$ flux $\ga 5 \rm mJy$, contribute 15% of the SFRD from all ULIRGs ($L_{\rm IR} \ga 10^{12} L_{\odot}$), and 3% of the total SFRD at $z \sim 2$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Radiation-Hydrodynamic Simulations of the Formation of Orion-Like Star Clusters II. The Initial Mass Function from Winds, Turbulence, and Radiation

[abridged] We report a series of simulations of the formation of a star cluster similar to the Orion Nebula Cluster (ONC), including both radiative transfer and protostellar outflows, and starting from both smooth and self-consistently turbulent initial conditions. Each simulation forms >150 stars and brown dwarfs, yielding a stellar mass distribution from < 0.1 to > 10 Msun. We show that a simulation that begins with self-consistently turbulence embedded in a larger turbulent volume, and that includes protostellar outflows, produces an initial mass function (IMF) consistent both with that of the ONC and the Galactic field. This is the first simulation published to date that reproduces the observed IMF in a cluster large enough to contain massive stars, and where the result is determined by a fully self-consistent calculation of gas thermodynamics. This simulation also produces a star formation rate that, while still somewhat too high, is much closer to observed values than if we omit either the larger turbulent volume or the outflows. Moreover, we show that the combination of outflows, self-consistently turbulent initial conditions, and turbulence continually fed by motions on scales larger than that of the protocluster yields an IMF that is in agreement with observations and invariant with time, resolving the "overheating" problem in which simulations without these features have an IMF peak that shifts to progressively higher masses over time. The simulation that matches the observed IMF also reproduces the observed trend of stellar multiplicity strongly increasing with mass. We show that this simulation produces massive stars from distinct massive cores whose properties are consistent with those of observed massive cores. However, the stars formed in these cores also undergo dynamical interactions that naturally produce Trapezium-like hierarchical multiple systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Radiatively Efficient Magnetized Bondi Accretion

We have carried out a numerical study of the effect of large scale magnetic fields on the rate of accretion from a uniform, isothermal gas onto a resistive, stationary point mass. Only mass, not magnetic flux, accretes onto the point mass. The simulations for this study avoid complications arising from boundary conditions by keeping the boundaries far from the accreting object. Our simulations leverage adaptive refinement methodology to attain high spatial fidelity close to the accreting object. Our results are particularly relevant to the problem of star formation from a magnetized molecular cloud in which thermal energy is radiated away on time scales much shorter than the dynamical time scale. Contrary to the adiabatic case, our simulations show convergence toward a finite accretion rate in the limit in which the radius of the accreting object vanishes, regardless of magnetic field strength. For very weak magnetic fields, the accretion rate first approaches the Bondi value and then drops by a factor ~ 2 as magnetic flux builds up near the point mass. For strong magnetic fields, the steady-state accretion rate is reduced by a factor ~ 0.2 β^{1/2} compared to the Bondi value, where βis the ratio of the gas pressure to the magnetic pressure. We give a simple expression for the accretion rate as a function of the magnetic field strength. Approximate analytic results are given in the Appendixes for both time-dependent accretion in the limit of weak magnetic fields and steady-state accretion for the case of strong magnetic fields.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Stable, Accurate Methodology for High Mach Number, Strong Magnetic Field MHD Turbulence with Adaptive Mesh Refinement: Resolution and Refinement Studies

Performing a stable, long duration simulation of driven MHD turbulence with a high thermal Mach number and a strong initial magnetic field is a challenge to high-order Godunov ideal MHD schemes because of the difficulty in guaranteeing positivity of the density and pressure. We have implemented a robust combination of reconstruction schemes, Riemann solvers, limiters, and Constrained Transport EMF averaging schemes that can meet this challenge, and using this strategy, we have developed a new Adaptive Mesh Refinement (AMR) MHD module of the ORION2 code. We investigate the effects of AMR on several statistical properties of a turbulent ideal MHD system with a thermal Mach number of 10 and a plasma $β_0$ of 0.1 as initial conditions; our code is shown to be stable for simulations with higher Mach numbers ($M_rms = 17.3$) and smaller plasma beta ($β_0 = 0.0067$) as well. Our results show that the quality of the turbulence simulation is generally related to the volume-averaged refinement. Our AMR simulations show that the turbulent dissipation coefficient for supersonic MHD turbulence is about 0.5, in agreement with unigrid simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Universal, Local Star Formation Law in Galactic Clouds, Nearby Galaxies, High-Redshift Disks, and Starbursts

[abridged] While observations of Local Group galaxies show a very simple, local star formation law in which the star formation rate per unit area in each patch of a galaxy scales linearly with the molecular gas surface density, recent observations of both Milky Way molecular clouds and high redshift galaxies apparently show a more complicated relationship, in which regions of equal surface density can form stars at quite different rates. These data have been interpreted as implying either that different star formation laws apply in different circumstances, that the star formation law is sensitive to large-scale galaxy properties rather than local properties, or that there are high density thresholds for star formation. Here we collate resolved observations of Milky Way molecular clouds, kpc-scale observations of Local Group galaxies, and unresolved observations of both disk and starburst galaxies in the local universe and at high redshift. We show that all of these data are in fact consistent with a simple, local, volumetric star formation law. The apparent variations stem from the fact that the observed objects have a wide variety of 3D size scales and degrees of internal clumping, so even at fixed gas column density the regions being observed can have wildly varying volume densities. We provide a simple theoretical framework to remove this projection effect, and we show that all the data, from small Solar neighborhood clouds with masses ~10^3 Msun to sub-mm galaxies with masses ~10^11 Msun, fall on a single star formation law in which the SFR is simply ~1% of the molecular gas mass per local free-fall time. In contrast, proposed star formation laws in which the star formation timescale is set by the galactic rotation period or the SFR is linearly proportional to the gas mass above some density threshold fail to match at least some of the data.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

IRAS 15099-5856: Remarkable Mid-Infrared Source with Prominent Crystalline Silicate Emission Embedded in the Supernova Remnant MSH15-52

We report new mid-infrared observations of the remarkable object IRAS 15099-5856 using the space telescopes AKARI and Spitzer, which demonstrate the presence of prominent crystalline silicate emission in this bright source. IRAS 15099-5856 has a complex morphology with a bright central compact source (IRS1) surrounded by knots, spurs, and several extended (~4') arc-like filaments. The source is seen only at >= 10 um. The Spitzer MIR spectrum of IRS1 shows prominent emission features from Mg-rich crystalline silicates, strong [Ne II] 12.81 um and several other faint ionic lines. We model the MIR spectrum as thermal emission from dust and compare with the Herbig Be star HD 100546 and the luminous blue variable R71, which show very similar MIR spectra. Molecular line observations reveal two molecular clouds around the source, but no associated dense molecular cores. We suggest that IRS1 is heated by UV radiation from the adjacent O star Muzzio 10 and that its crystalline silicates most likely originated in a mass outflow from the progenitor of the supernova remnant (SNR) MSH 15-52. IRS1, which is embedded in the SNR, could have been shielded from the SN blast wave if the progenitor was in a close binary system with Muzzio 10. If MSH15-52 is a remnant of Type Ib/c supernova (SN Ib/c), as has been previously proposed, this would confirm the binary model for SN Ib/c. IRS1 and the associated structures may be the relics of massive star death, as shaped by the supernova explosion, the pulsar wind and the intense ionizing radiation of the embedded O star.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Metallicity and the Universality of the IMF

The stellar initial mass function (IMF), along with the star formation rate, is one of the fundamental properties that any theory of star formation must explain. An interesting feature of the IMF is that it appears to be remarkably universal across a wide range of environments. Particularly, there appears to be little variation in either the characteristic mass of the IMF or its high-mass tail between clusters with different metallicities. Previous attempts to understand this apparent independence of metallicity have not accounted for radiation feedback from high-mass protostars, which can dominate the energy balance of the gas in star-forming regions. We extend this work, showing that the fragmentation of molecular gas should depend only weakly on the amount of dust present, even when the primary heating source is radiation from massive protostars. First, we report a series of core collapse simulations using the ORION AMR code that systematically vary the dust opacity and show explicitly that this has little effect on the temperature or fragmentation of the gas. Then, we provide an analytic argument for why the IMF varies so little in observed star clusters, even as the metallicity varies by a factor of 100.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Radiation-Hydrodynamic Simulations of Massive Star Formation with Protostellar Outflows

We report the results of a series of AMR radiation-hydrodynamic simulations of the collapse of massive star forming clouds using the ORION code. These simulations are the first to include the feedback effects protostellar outflows, as well as protostellar radiative heating and radiation pressure exerted on the infalling, dusty gas. We find that that outflows evacuate polar cavities of reduced optical depth through the ambient core. These enhance the radiative flux in the poleward direction so that it is 1.7 to 15 times larger than that in the midplane. As a result the radiative heating and outward radiation force exerted on the protostellar disk and infalling cloud gas in the equatorial direction are greatly diminished. The simultaneously reduces the Eddington radiation pressure barrier to high-mass star formation and increases the minimum threshold surface density for radiative heating to suppress fragmentation compared to models that do not include outflows. The strength of both these effects depends on the initial core surface density. Lower surface density cores have longer free-fall times and thus massive stars formed within them undergo more Kelvin contraction as the core collapses, leading to more powerful outflows. Furthermore, in lower surface density clouds the ratio of the time required for the outflow to break out of the core to the core free-fall time is smaller, so that these clouds are consequently influenced by outflows at earlier stages of collapse. As a result, outflow effects are strongest in low surface density cores and weakest in high surface density one. We also find that radiation focusing in the direction of outflow cavities is sufficient to prevent the formation of radiation pressure-supported circumstellar gas bubbles, in contrast to models which neglect protostellar outflow feedback.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Radiation-Hydrodynamic Simulations of the Formation of Orion-Like Star Clusters I. Implications for the Origin of the Initial Mass Function

One model for the origin of typical galactic star clusters such as the Orion Nebula Cluster (ONC) is that they form via the rapid, efficient collapse of a bound gas clump within a larger, gravitationally-unbound giant molecular cloud. However, simulations in support of this scenario have thus far have not included the radiation feedback produced by the stars; radiative simulations have been limited to significantly smaller or lower density regions. Here we use the ORION adaptive mesh refinement code to conduct the first ever radiation-hydrodynamic simulations of the global collapse scenario for the formation of an ONC-like cluster. We show that radiative feedback has a dramatic effect on the evolution: once the first ~10-20% of the gas mass is incorporated into stars, their radiative feedback raises the gas temperature high enough to suppress any further fragmentation. However, gas continues to accrete onto existing stars, and, as a result, the stellar mass distribution becomes increasingly top-heavy, eventually rendering it incompatible with the observed IMF. Systematic variation in the location of the IMF peak as star formation proceeds is incompatible with the observed invariance of the IMF between star clusters, unless some unknown mechanism synchronizes the IMFs in different clusters by ensuring that star formation is always truncated when the IMF peak reaches a particular value. We therefore conclude that the global collapse scenario, at least in its simplest form, is not compatible with the observed stellar IMF. We speculate that processes that slow down star formation, and thus reduce the accretion luminosity, may be able to resolve the problem.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Sub-Alfvenic Non-Ideal MHD Turbulence Simulations with Ambipolar Diffusion: III. Implications for Observations and Turbulent Enhancement

Ambipolar diffusion (AD) is believed to be a crucial process for redistributing magnetic flux in the dense molecular gas that occurs in regions of star formation. We carry out numerical simulations of this process in regions of low ionization using the heavy ion approximation. The simulations are for regions of strong field (plasma β=0.1) and mildly supersonic turbulence (M=3, corresponding to an Alfven mach number of 0.67). The velocity power spectrum of the neutral gas changes from an Iroshnikov-Kraichnan spectrum in the case of ideal MHD to a Burgers spectrum in the case of a shock-dominated hydrodynamic system. The magnetic power spectrum shows a similar behavior. We use a 1D radiative transfer code to post-process our simulation results; the simulated emission from the CS J=2-1 and H13CO+ J=1-0 lines shows that the effects of AD are observable in principle. Linewidths of ions are observed to be less than those of neutrals, and we confirm previous suggestions that this is due to AD. We show that AD is unlikely to affect the Chandrasekhar-Fermi method for inferring field strengths unless the AD is stronger than generally observed. Finally, we present the first fully 3D study of the enhancement of AD by turbulence, finding that AD is accelerated by factor 2-4.5 for non self-gravitating systems with the level of turbulence we consider.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Global Evolution of Giant Molecular Clouds II: The Role of Accretion

We present virial models for the global evolution of giant molecular clouds. Focusing on the presence of an accretion flow, and accounting for the amount of mass, momentum, and energy supplied by accretion and star formation feedback, we are able to follow the growth, evolution, and dispersal of individual giant molecular clouds. Our model clouds reproduce the scaling relations observed in both galactic and extragalactic clouds. We find that accretion and star formation contribute contribute roughly equal amounts of turbulent kinetic energy over the lifetime of the cloud. Clouds attain virial equilibrium and grow in such a way as to maintain roughly constant surface densities, with typical surface densities of order 50 - 200 Msun pc^-2, in good agreement with observations of giant molecular clouds in the Milky Way and nearby external galaxies. We find that as clouds grow, their velocity dispersion and radius must also increase, implying that the linewidth-size relation constitutes an age sequence. Lastly, we compare our models to observations of giant molecular clouds and associated young star clusters in the LMC and find good agreement between our model clouds and the observed relationship between H ii regions, young star clusters, and giant molecular clouds.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Protostellar Luminosity Function

The protostellar luminosity function (PLF) is the present-day luminosity function of the protostars in a region of star formation. It is determined using the protostellar mass function (PMF) in combination with a stellar evolutionary model that provides the luminosity as a function of instantaneous and final stellar mass. As in McKee & Offner (2010), we consider three main accretion models: the Isothermal Sphere model, the Turbulent Core model, and an approximation of the Competitive Accretion model. We also consider the effect of an accretion rate that tapers off linearly in time and an accelerating star formation rate. For each model, we characterize the luminosity distribution using the mean, median, maximum, ratio of the median to the mean, standard deviation of the logarithm of the luminosity, and the fraction of very low luminosity objects. We compare the models with bolometric luminosities observed in local star forming regions and find that models with an approximately constant accretion time, such as the Turbulent Core and Competitive Accretion models, appear to agree better with observation than those with a constant accretion rate, such as the Isothermal Sphere model. We show that observations of the mean protostellar luminosity in these nearby regions of low-mass star formation suggest a mean star formation time of 0.3$\pm$0.1 Myr. Such a timescale, together with some accretion that occurs non-radiatively and some that occurs in high-accretion, episodic bursts, resolves the classical "luminosity problem" in low-mass star formation, in which observed protostellar luminosities are significantly less than predicted. An accelerating star formation rate is one possible way of reconciling the observed star formation time and mean luminosity.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

What Phase of the Interstellar Medium Correlates with the Star Formation Rate?

Nearby spiral galaxies show an extremely tight correlation between tracers of molecular hydrogen (H_2) in the interstellar medium (ISM) and tracers of recent star formation, but it is unclear whether this correlation is fundamental or accidental. In the galaxies that have been surveyed to date, H_2 resides predominantly in gravitationally bound clouds cooled by carbon monoxide (CO) molecules, but in galaxies of low metal content the correlations between bound clouds, CO, and H_2 break down, and it is unclear if the star formation rate will then correlate with H_2 or with some other quantity. Here we show that star formation will continue to follow H_2 independent of metallicity. This is not because H_2 is directly important for cooling, but instead because the transition from predominantly atomic hydrogen (HI) to H_2 occurs under the same conditions as a dramatic drop in gas temperature and Bonnor-Ebert mass that destabilizes clouds and initiates collapse. We use this model to compute how star formation rate will correlate with total gas mass, with mass of gas where the hydrogen is H_2, and with mass of gas where the carbon is CO in galaxies of varying metallicity, and show that preliminary observations match the trend we predict.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

An Initial Mass Function for Individual Stars in Galactic Disks: I. Constraining the Shape of the IMF

We derive a semi-empirical galactic initial mass function (IMF) from observational constraints. We assume that the star formation rate in a galaxy can be expressed as the product of the IMF, $ψ(m)$, which is a smooth function of mass $m$ (in units of \msun), and a time- and space-dependent total rate of star formation per unit area of galactic disk. The mass dependence of the proposed IMF is determined by five parameters: the low-mass slope $γ$, the high-mass slope $-Γ$, the characteristic mass $m_{ch}$ (which is close to the mass $m_{\rm peak}$ at which the IMF turns over), and the lower and upper limits on the mass, $m_l$ (taken to be 0.004) and $m_u$ (taken to be 120). The star formation rate in terms of number of stars per unit area of galactic disk per unit logarithmic mass interval, is proportional to $m^{-Γ} \left\{1-\exp\left[{-(m/m_{ch})^{γ+Γ}}\right]\right\}$, where $\cal N_*$ is the number of stars, $m_l<m<m_u$ is the range of stellar masses. The values of $γ$ and $\emch$ are derived from two integral constraints: i) the ratio of the number density of stars in the range $m=0.1-0.6$ to that in the range $m=0.6-0.8$ as inferred from the mass distribution of field stars in the local neighborhood, and ii) the ratio of the number of stars in the range $m=0.08 - 1$ to the number of brown dwarfs in the range $m=0.03-0.08$ in young clusters. The IMF satisfying the above constraints is characterized by the parameters $γ=0.51$ and $\emch=0.35$ (which corresponds to $m_{\rm peak}=0.27$). This IMF agrees quite well with the Chabrier (2005) IMF for the entire mass range over which we have compared with data, but predicts significantly more stars with masses $< 0.03\, M_\odot$; we also compare with other IMFs in current use.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Metal-Ion Absorption in Conductively Evaporating Clouds

We present computations of the ionization structure and metal-absorption properties of thermally conductive interface layers that surround evaporating warm spherical clouds, embedded in a hot medium. We rely on the analytical formalism of Dalton & Balbus to calculate the temperature profile in the evaporating gas, and explicitly solve the time-dependent ionization equations for H, He, C, N, O, Si, and S in the interface. We include photoionization by an external field. We estimate how departures from equilibrium ionization affect the resonance-line cooling efficiencies in the evaporating gas, and determine the conditions for which radiative losses may be neglected in the solution for the evaporation dynamics and temperature profile. Our results indicate that non-equilibrium cooling significantly increases the value of the saturation parameter at which radiative losses begin to affect the flow dynamics. As applications we calculate the ion fractions and projected column densities arising in the evaporating layers surrounding dwarf-galaxy-scale objects that are also photoionized by metagalactic radiation. We compare our results to the UV metal-absorption column densities observed in local highly-ionized metal-absorbers, located in the Galactic corona or intergalactic medium. Conductive interfaces significantly enhance the formation of high-ions such as C^3+, N^4+, and O^5+ relative to purely photoionized clouds, especially for clouds embedded in a high-pressure corona. However, the enhanced columns are still too low to account for the O VI columns (~1e14 cm^-2) observed in the local high-velocity absorbers. We find that O VI columns larger than ~1e13 cm^-2 cannot be produced in evaporating clouds. Our results do support the conclusion of Savage & Lehner, that absorption due to evaporating O VI likely occurs in the local interstellar medium, with characteristic columns of ~1e13 cm^-2.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Radiation Feedback, Fragmentation, and the Environmental Dependence of the Initial Mass Function

The fragmentation of star-forming interstellar clouds, and the resulting stellar initial mass function (IMF), is strongly affected by the temperature structure of the collapsing gas. Since radiation feedback from embedded stars can modify this as collapse proceeds, feedback plays an important role in determining the IMF. However, the effects and importance of radiative heating are likely to depend strongly on the surface density of the collapsing clouds, which determines both their effectiveness at trapping radiation and the accretion luminosities of the stars forming within them. In this paper we report a suite of adaptive mesh refinement radiation-hydrodynamic simulations using the ORION code in which we isolate the effect of column density on fragmentation by following the collapse of clouds of varying column density while holding the mass, initial density and velocity structure, and initial virial ratio fixed. We find that radiation does not significantly modify the overall star formation rate or efficiency, but that it suppresses fragmentation more and more as cloud surface densities increase from those typical of low mass star-forming regions like Taurus, through the typical surface density of massive star-forming clouds in the Galaxy, up to conditions found only in super star clusters. In regions of low surface density, fragmentation during collapse leads to the formation of small clusters rather than individual massive star systems, greatly reducing the fraction of the stellar population with masses >~ 10 Msun. Our simulations have important implications for the formation of massive stars and the universality of the IMF.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Regulation of Star Formation Rates in Multiphase Galactic Disks: a Thermal/Dynamical Equilibrium Model

We develop a model for regulation of galactic star formation rates Sigma_SFR in disk galaxies, in which ISM heating by stellar UV plays a key role. By requiring simultaneous thermal and (vertical) dynamical equilibrium in the diffuse gas, and star formation at a rate proportional to the mass of the self-gravitating component, we obtain a prediction for Sigma_SFR as a function of the total gaseous surface density Sigma and the density of stars + dark matter, rho_sd. The physical basis of this relationship is that thermal pressure in the diffuse ISM, which is proportional to the UV heating rate and therefore to Sigma_SFR, must adjust to match the midplane pressure set by the vertical gravitational field. Our model applies to regions where Sigma < 100 Msun/pc^2. In low-Sigma_SFR (outer-galaxy) regions where diffuse gas dominates, the theory predicts Sigma_SFR \propto Sigma (rho_sd)^1/2. The decrease of thermal equilibrium pressure when Sigma_SFR is low implies, consistent with observations, that star formation can extend (with declining efficiency) to large radii in galaxies, rather than having a sharp cutoff. The main parameters entering our model are the ratio of thermal pressure to total pressure in the diffuse ISM, the fraction of diffuse gas that is in the warm phase, and the star formation timescale in self-gravitating clouds; all of these are (in principle) direct observables. At low surface density, our model depends on the ratio of the mean midplane FUV intensity (or thermal pressure in the diffuse gas) to the star formation rate, which we set based on Solar neighborhood values. We compare our results to recent observations, showing good agreement overall for azimuthally-averaged data in a set of spiral galaxies. For the large flocculent spiral galaxies NGC 7331 and NGC 5055, the correspondence between theory and observation is remarkably close.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Dark Molecular Gas

The mass of molecular gas in an interstellar cloud is often measured using line emission from low rotational levels of CO, which are sensitive to the CO mass, and then scaling to the assumed molecular hydrogen H_2 mass. However, a significant H_2 mass may lie outside the CO region, in the outer regions of the molecular cloud where the gas phase carbon resides in C or C+. Here, H_2 self-shields or is shielded by dust from UV photodissociation, where as CO is photodissociated. This H_2 gas is "dark" in molecular transitions because of the absence of CO and other trace molecules, and because H_2 emits so weakly at temperatures 10 K < T < 100 K typical of this molecular component. This component has been indirectly observed through other tracers of mass such as gamma rays produced in cosmic ray collisions with the gas and far-infrared/submillimeter wavelength dust continuum radiation. In this paper we theoretically model this dark mass and find that the fraction of the molecular mass in this dark component is remarkably constant (~ 0.3 for average visual extinction through the cloud with mean A_V ~ 8) and insensitive to the incident ultraviolet radiation field strength, the internal density distribution, and the mass of the molecular cloud as long as mean A_V, or equivalently, the product of the average hydrogen nucleus column and the metallicity through the cloud, is constant. We also find that the dark mass fraction increases with decreasing mean A_V, since relatively more molecular H_2 material lies outside the CO region in this case.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Protostellar Mass Function

The protostellar mass function (PMF) is the Present-Day Mass Function of the protostars in a region of star formation. It is determined by the initial mass function weighted by the accretion time. The PMF thus depends on the accretion history of protostars and in principle provides a powerful tool for observationally distinguishing different protostellar accretion models. We consider three basic models here: the Isothermal Sphere model (Shu 1977), the Turbulent Core model (McKee & Tan 2003), and an approximate representation of the Competitive Accretion model (Bonnell et al. 1997, 2001a). We also consider modified versions of these accretion models, in which the accretion rate tapers off linearly in time. Finally, we allow for an overall acceleration in the rate of star formation. At present, it is not possible to directly determine the PMF since protostellar masses are not currently measurable. We carry out an approximate comparison of predicted PMFs with observation by using the theory to infer the conditions in the ambient medium in several star-forming regions. Tapered and accelerating models generally agree better with observed star-formation times than models without tapering or acceleration, but uncertainties in the accretion models and in the observations do not allow one to rule out any of the proposed models at present. The PMF is essential for the calculation of the Protostellar Luminosity Function, however, and this enables stronger conclusions to be drawn (Offner & McKee 2010).

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Atomic-to-Molecular Transition in Galaxies. III. A New Method for Determining the Molecular Content of Primordial and Dusty Clouds

Understanding the molecular content of galaxies is a critical problem in star formation and galactic evolution. Here we present a new method, based on a Stromgren-type analysis, to calculate the amount of HI that surrounds a molecular cloud irradiated by an isotropic radiation field. We consider both planar and spherical clouds, and H_2 formation either in the gas phase or catalyzed by dust grains. Under the assumption that the transition from atomic to molecular gas is sharp, our method gives the solution without any reference to the photodissociation cross section. We test our results for the planar case against those of a PDR code, and find typical accuracies of about 10%. Our results are also consistent with the scaling relations found in Paper I of this series, but they apply to a wider range of physical conditions. We present simple, accurate analytic fits to our results that are suitable for comparison to observations and to implementation in numerical and semi-analytic models.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Formation of Massive Star Systems by Accretion

Massive stars produce so much light that the radiation pressure they exert on the gas and dust around them is stronger than their gravitational attraction, a condition that has long been expected to prevent them from growing by accretion. We present three-dimensional radiation-hydrodynamic simulations of the collapse of a massive prestellar core and find that radiation pressure does not halt accretion. Instead, gravitational and Rayleigh-Taylor instabilities channel gas onto the star system through non-axisymmetric disks and filaments that self-shield against radiation, while allowing radiation to escape through optically-thin bubbles. Gravitational instabilities cause the disk to fragment and form a massive companion to the primary star. Radiation pressure does not limit stellar masses, but the instabilities that allow accretion to continue lead to small multiple systems.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The formation of the first stars and galaxies

Observations made using large ground-based and space-borne telescopes have probed cosmic history all the way from the present-day to a time when the Universe was less than a tenth of its present age. Earlier on lies the remaining frontier, where the first stars, galaxies, and massive black holes formed. They fundamentally transformed the early Universe by endowing it with the first sources of light and chemical elements beyond the primordial hydrogen and helium produced in the Big Bang. The interplay of theory and upcoming observations promises to answer the key open questions in this emerging field.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Star Formation Law in Atomic and Molecular Gas

We propose a simple theoretical model for star formation in which the local star formation rate in a galaxy is determined by three factors. First, the interplay between the interstellar radiation field and molecular self-shielding determines what fraction of the gas is in molecular form and thus eligible to form stars. Second, internal feedback determines the properties of the molecular clouds that form, which are nearly independent of galaxy properties until the galactic ISM pressure becomes comparable to the internal GMC pressure. Above this limit, galactic ISM pressure determines molecular gas properties. Third, the turbulence driven by feedback processes in GMCs makes star formation slow, allowing a small fraction of the gas to be converted to stars per free-fall time within the molecular clouds. We combine analytic estimates for each of these steps to formulate a single star formation law, and show that the predicted correlation between star formation rate, metallicity, and surface densities of atomic, molecular, and total gas agree well with observations.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

First Stars III Conference Summary

The understanding of the formation, life, and death of Population III stars, as well as the impact that these objects had on later generations of structure formation, is one of the foremost issues in modern cosmological research and has been an active area of research during the past several years. We summarize the results presented at "First Stars III," a conference sponsored by Los Alamos National Laboratory, the Kavli Institute for Particle Astrophysics and Cosmology, and the Joint Institute for Nuclear Astrophysics. This conference, the third in a series, took place in July 2007 at the La Fonda Hotel in Santa Fe, New Mexico, U.S.A.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

A General Theory of Turbulence-Regulated Star Formation, From Spirals to ULIRGs

We derive an analytic prediction for the star formation rate in environments ranging from normal galactic disks to starbursts and ULIRGs in terms of the observables of those systems. Our calculation is based on three premises: (1) star formation occurs in virialized molecular clouds that are supersonically turbulent; (2) the density distribution within these clouds is lognormal, as expected for supersonic isothermal turbulence; (3) stars form in any sub-region of a cloud that is so overdense that its gravitational potential energy exceeds the energy in turbulent motions. We show that a theory based on this model is consistent with simulations and with the observed star formation rate in the Milky Way. We use our theory to derive the Kennicutt-Schmidt Law from first principles, and make other predictions that can be tested by future observations. We also provide an algorithm for estimating the star formation rate that is suitable for inclusion in numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Massive star formation in 100,000 years from turbulent and pressurized molecular clouds

Massive stars (with mass m_* > 8 solar masses) are fundamental to the evolution of galaxies, because they produce heavy elements, inject energy into the interstellar medium, and possibly regulate the star formation rate. The individual star formation time, t_*f, determines the accretion rate of the star; the value of the former quantity is currently uncertain by many orders of magnitude, leading to other astrophysical questions. For example, the variation of t_*f with stellar mass dictates whether massive stars can form simultaneously with low-mass stars in clusters. Here we show that t_*f is determined by conditions in the star's natal cloud, and is typically ~10^5 yr. The corresponding mass accretion rate depends on the pressure within the cloud - which we relate to the gas surface density - and on both the instantaneous and final stellar masses. Characteristic accretion rates are sufficient to overcome radiation pressure from ~100 solar mass protostars, while simultaneously driving intense bipolar gas outflows. The weak dependence of t_*f on the final mass of the star allows high- and low-mass star formation to occur nearly simultaneously in clusters.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Multi--Pressure Polytropes as Models for the Structure and Stability of Molecular Clouds. I. Theory

Molecular clouds are supported by thermal pressure, magnetic pressure, and turbulent pressure. Each of these can be modeled with a polytropic equation of state, so that overall the total pressure is the sum of the individual components. We model the turbulent pressure as being due to a superposition of Alfven waves. The theory of polytropes is generalized to allow for the flow of entropy in response to a perturbation, as expected for the entropy associated with wave pressure. The equation of state of molecular clouds is "soft", so that the properties of the clouds are generally governed by the conditions at the surface. In general, the polytropes are not isentropic, and this permits large density and pressure drops to occur between the center and the edge of the polytropes, as is observed.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Gravitational Collapse and Fragmentation in Molecular Clouds with Adaptive Mesh Refinement

We describe a powerful methodology for numerical solution of 3-D self-gravitational hydrodynamics problems with extremely high resolution. Our method utilizes the technique of local adaptive mesh refinement (AMR), employing multiple grids at multiple levels of resolution. These grids are automatically and dynamically added and removed as necessary to maintain adequate resolution. This technology allows for the solution of problems in a manner that is both more efficient and more versatile than other fixed and variable resolution methods. The application of AMR to simulate the collapse and fragmentation of a molecular cloud, a key step in star formation, is discussed. Such simulations involve many orders of magnitude of variation in length scale as fragments form. In this paper we briefly describe the methodology and present an illustrative application for nonisothermal cloud collapse. We describe the numerical Jeans condition, a criterion for stability of self-gravitational hydrodynamics problems. We show the first well-resolved nonisothermal evolutionary sequence beginning with a perturbed dense molecular cloud core that leads to the formation of a binary system consisting of protostellar cores surrounded by distinct protostellar disks. The scale of the disks, of order 100 AU, is consistent with observations of gaseous disks surrounding single T-Tauri stars and debris disks surrounding systems such as $β$ Pictoris.

0 citations0 reviewsNo code linkedOpen access