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

Richard I. Klein

Richard I. Klein appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRastro-ph.GAastro-ph.COastro-ph.EPastro-ph.IMphysics.comp-phastro-phastro-ph.HEphysics.flu-dynphysics.plasm-ph

Catalog footprint

What is connected

25works
10topics
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

25 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
preprint2022arXiv

The Davis-Chandrasekhar-Fermi Method Revisited

Despite the rich observational results on interstellar magnetic fields in star-forming regions, it is still unclear how dynamically significant the magnetic fields are at varying physical scales, because direct measurement of the field strength is observationally difficult. The Davis-Chandrasekhar-Fermi (DCF) method has been the most commonly used method to estimate the magnetic field strength from polarization data. It is based on the assumption that gas turbulent motion is the driving source of field distortion via linear Alfvén waves. In this work, using MHD simulations of star-forming clouds, we test the validity of the assumption underlying the DCF method by examining its accuracy in the real 3D space. Our results suggest that the DCF relation between turbulent kinetic energy and magnetic energy fluctuation should be treated as a statistical result instead of a local property. We then develop and investigate several modifications to the original DCF method using synthetic observations, and propose new recipes to improve the accuracy of DCF-derived magnetic field strength. We further note that the biggest uncertainty in the DCF analysis may come from the linewidth measurement instead of the polarization observation, especially since the line-of-sight gas velocity can be used to estimate the gas volume density, another critical parameter in the DCF method.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Role of Magnetic Fields in the Stability and Fragmentation of Filamentary Molecular Clouds: Two Case Studies at OMC-3 and OMC-4

We present the stability analysis of two regions, OMC-3 and OMC-4, in the massive and long molecular cloud complex of Orion A. We obtained $214~μ$m HAWC+/SOFIA polarization data, and we make use of archival data for the column density and C$^{18}$O (1-0) emission line. We find clear depolarization in both observed regions and that the polarization fraction is anti-correlated with the column density and the polarization-angle dispersion function. We find that the filamentary cloud and dense clumps in OMC-3 are magnetically supercritical and strongly subvirial. This region should be in the gravitational collapse phase and is consistent with many young stellar objects (YSOs) forming in the region. Our histogram of relative orientations (HROs) analysis shows that the magnetic field is dynamically sub-dominant in the dense gas structures of OMC-3. We present the first polarization map of OMC-4. We find that the observed region is generally magnetically subcritical except for an elongated dense core, which could be a result of projection effect of a filamentary structure aligned close to the line-of-sight. The relative large velocity dispersion and the unusual positive shape parameters at high column densities in the HROs analysis suggest that our viewing angle may be close to axes of filamentary substructures in OMC-4. The dominating strong magnetic field in OMC-4 is unfavorable for star formation and is consistent with much fewer YSOs than in OMC-3.

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

Hybrid Adaptive Ray-Moment Method (HARM$^2$): A Highly Parallel Method for Radiation Hydrodynamics on Adaptive Grids

We present a highly-parallel multi-frequency hybrid radiation hydrodynamics algorithm that combines a spatially-adaptive long characteristics method for the radiation field from point sources with a moment method that handles the diffuse radiation field produced by a volume-filling fluid. Our Hybrid Adaptive Ray-Moment Method (HARM$^2$) operates on patch-based adaptive grids, is compatible with asynchronous time stepping, and works with any moment method. In comparison to previous long characteristics methods, we have greatly improved the parallel performance of the adaptive long-characteristics method by developing a new completely asynchronous and non-blocking communication algorithm. As a result of this improvement, our implementation achieves near-perfect scaling up to $\mathcal{O}(10^3)$ processors on distributed memory machines. We present a series of tests to demonstrate the accuracy and performance of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Three temperature plasma shock solutions with gray radiation diffusion

The effects of radiation on the structure of shocks in a fully-ionized plasma are investigated by solving the steady-state fluid equations for ions, electrons, and radiation. The electrons and ions are assumed to have the same bulk velocity but separate temperatures, and the radiation is modeled with the gray-diffusion approximation. Both electron and ion conduction are included, as well as ion viscosity. When the material is optically thin, three-temperature behavior occurs. When the diffusive flux of radiation is important but radiation pressure is not, two-temperature behavior occurs, with the electrons strongly coupled to the radiation. Since the radiation heats the electrons on length scales that are much longer than the electron-ion Coulomb coupling length scale, these solutions resemble radiative shock solutions rather than plasma shock solutions that neglect radiation. When radiation pressure is important, all three components are strongly coupled. Results with constant values for the transport and coupling coefficients are compared to a full numerical simulation with a good match between the two, demonstrating that steady shock solutions constitute a straightforward and comprehensive verification test methodology for multi-physics numerical algorithms.

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

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

Star Cluster Formation and Feedback

Stars do not generally form in isolation. Instead, they form in clusters, and in these clustered environments newborn stars can have profound effects on one another and on their parent gas clouds. Feedback from clustered stars is almost certainly responsible for a number of otherwise puzzling facts about star formation: that it is an inefficient process that proceeds slowly when averaged over galactic scales; that most stars disperse from their birth sites and dissolve into the galactic field over timescales $\ll 1$ Gyr; and that newborn stars follow an initial mass function (IMF) with a distinct peak in the range $0.1 - 1$ $M_\odot$, rather than an IMF dominated by brown dwarfs. In this review we summarize current observational constraints and theoretical models for the complex interplay between clustered star formation and feedback.

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

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

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
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
preprint2009arXiv

On the Role of Disks in the Formation of Stellar Systems: A Numerical Parameter Study of Rapid Accretion

We study rapidly accreting, gravitationally unstable disks with a series of global, three dimensional, numerical experiments using the code ORION. In this paper we conduct a numerical parameter study focused on protostellar disks, and show that one can predict disk behavior and the multiplicity of the accreting star system as a function of two dimensionless parameters which compare the disk's accretion rate to its sound speed and orbital period. Although gravitational instabilities become strong, we find that fragmentation into binary or multiple systems occurs only when material falls in several times more rapidly than the canonical isothermal limit. The disk-to-star accretion rate is proportional to the infall rate, and governed by gravitational torques generated by low-m spiral modes. We also confirm the existence of a maximum stable disk mass: disks that exceed ~50% of the total system mass are subject to fragmentation and the subsequent formation of binary companions.

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