the formation of high mass stars richard i. klein uc berkeley, department of astronomy and lawrence...

The Formation of High Mass StarsThe Formation of High Mass Stars

Richard I. Klein

UC Berkeley, Department of Astronomy

and

Lawrence Livermore National Laboratory

Collaborators

Mark Krumholz (Princeton University) and Chris McKee (UC Berkeley)

This work was performed under the auspices of the U.S. Department of Energy by the University of California Lawrence LivermoreNational Laboratory under contract No. W-7405-Eng-48. UCRL-PRES-229278

Zurich September 17, 2007Next Generation of Computational Models

UCRL-PRES-229278- 2

Outstanding Challenges of Massive Star Formation

Outstanding Challenges of Massive Star Formation

• What is the formation Mechanism: Gravitational collapse of an unstable turbulent cloud; Competitive Bondi-Hoyle accretion; Collisional Coalescence?

• How can gravitationally collapsing clouds overcome the Eddington limit due to radiation pressure?

• What determines the upper limit for High Mass Stars? (120Msun 150Msun)

• How do feedback mechanisms such as protostellar outflows and radiation affect protostellar evolution? These mechanisms can also have a dramatic effect on cluster formation

• How do the systems in which massive stars are present form?

UCRL-PRES-229278- 3

Theoretical Challenges of High Mass Star FormationTheoretical Challenges of High Mass Star Formation

1. Effects of Strong Radiation Pressure

— Massive stars M 20 M have tK < tform (Shu et al. 1987) and begin nuclear burning during accretion phase

Radiates enormous energy

For M 100 M

however dust >> T

But, observations show M ~ 100 M (Massey 1998, 2003)

Fundamental Problem: How is it possible to sustain a sufficiently high-mass accretion rate onto protostellar core despite “Eddington” barrier?

Does radiation pressure provide a natural limit to the formation of high mass stars?

ÏMMforff gravrad 10>>

ÏL~cMGm

L~LT

pedd*

61034

×

π=

UCRL-PRES-229278- 4

Theoretical Challenges of High Mass Star Formation (cont.)

Theoretical Challenges of High Mass Star Formation (cont.)

2. Effects of Protostellar outflows

— Massive stars produce strong radiation driven stellar winds with momentum fluxes

— Massive YSO have observed (CO) protostellar outflows where (Richer et al. 2000; Cesaroni 2004)

If outflows where spherically symmetric this would create a greater obstacle to massive star formation than radiation pressure

but, flows are found to be collimated with collimation factors 2-10 (Beuther 2002, 2003, 2004)

Fundamental Problem: How do outflows effect the formation of Massive stars? Do outflows limit the mass of a star?

€

˙ M v ≤ L /c

€

˙ M v ~ 100L /c

UCRL-PRES-229278- 5

Physical Effects in High-Mass Star FormationPhysical Effects in High-Mass Star Formation

• Photoionization:

— Effects quenched for moderate accretion rates ~10-4 M /yr

for spherically symmetric infall

— In disk accretion, material above and below disk confine ionized region close to stellar surface

— Outflows are sufficient to quench ionization (Tan and McKee 2003)

Omit Photoionization as a first approximation

• Magnetic fields

— Gravity dominates magnetic fields when M > MB B3/n2 so magnetic fields are dynamically unimportant for high mass cores (Shu et al. 1987)

— At high densities B n1/2 so MB n-1/2 and since n 106 in massive star forming regions, MB is substantially reduced

— Observations of magnetic fields in high mass cores inconlcusive

Neglect of magnetic fields is a reasonable first approximation

€

˙ M crit =4πGM*Sm2

H

α (2)≈ 2 ×10−5 M*

10Msun

⎛

⎝ ⎜

⎞

⎠ ⎟

1/ 2S

1049 s−1

⎛

⎝ ⎜

⎞

⎠ ⎟

1/ 2

Msun yr−1

UCRL-PRES-229278- 6

Physical Effects in High-Mass Star Formation (cont.)Physical Effects in High-Mass Star Formation (cont.)

• Dust

— Critical role in massive star formation couples gas to radiation flux from central star

need radiation transport and multi-species models with good microphysics

• Photostellar outflows

— Molecular outflows in neighborhood of massive stars ~ 10-4 - 10-2 M /yr. Force required to drive such outflows Fco > 10 – 100 LBOL/c

Outflows may be important to protostellar evolution

• Three-Dimensional Effects

— Interaction of radiation with infalling envelope subject to radiation driven instabilities

— Interaction of protostellar outflow with infalling envelope possibly unstable

— Accretion disks develop non-axisymmetric structures in turbulent flows

Three dimensional simulations are crucial

UCRL-PRES-229278- 7

Equations of Gravito-Radiation Hydrodynamics to order v/c (Krumholz, Klein & McKee 2007a)

Equations of Gravito-Radiation Hydrodynamics to order v/c (Krumholz, Klein & McKee 2007a)

( ) 0=ρν⋅∇+∂

ρ∂

t

( ) ( ) Fc

Pt

Rκ+φ∇ρ−−∇=ρνν⋅∇+ρν

∂

∂

( ) ( )[ ] ( ) Fc

ceet

R

R

P ⋅νκ

⎟⎟⎠

⎞⎜⎜⎝

⎛−

κ

κ−Ε−Βπκ−φ∇⋅ρν−=νρ+ρ⋅∇+ρ

∂

∂Ρ 124

( )⎥⎦

⎤⎢⎣

⎡−δ+ρπ=φ∇ ∑

iii xxMG42

( ) ( ) ⎟⎠

⎞⎜⎝

⎛ ν−

⋅∇−∇⋅ν⎟⎟⎠

⎞⎜⎜⎝

⎛−

κκ

λ−−Βπκ+−δ=⎟⎟⎠

⎞⎜⎜⎝

⎛∇

κλ−

⋅∇+∂∂

Ρ∑ ER

EcExxLEc

t

E

R

P

iii

R 2

3124 2

( )E

EcF

R

∇κ

λ−=

(Continuity)

(Gas momentum)

(Gas energy)

(Poisson)

(Radiation energy)

(Flux-limited diffusion approximation)

Equations exact to (v/c) in static diffusion regime.

UCRL-PRES-229278- 8

High Mass Star Formation Simulation PhysicsHigh Mass Star Formation Simulation Physics

• Euler equations of compressible gas dynamics with gravity

• Radiative transfer and radiation pressure in the gray, flux-limited diffusion approximation radiative feedback

• Model of dust opacity based on Pollack et al. (1994) (6 species)

• Outflows: hydromagnetic outflow models

Dynamical Feedback

• Eulerian sink particles:

— Created when the density in a cell exceeds the local Jeans density (Krumholz, McKee, & Klein 2004)

— Free to move through the grid and continue to accrete gas

— Sink particles feed radiation and (for some runs) winds back into the grid based on a protostellar model

— Model includes accretion, KH contraction, deuterium and hydrogen burning (McKee & Tan 2003), x-winds

• Capability to handle the enormous range of scales involved AMR

UCRL-PRES-229278- 9

Physics ImplementationPhysics Implementation

Our AMR code is a combination of C++ and FORTRAN 90 Uses parallel MPI-based Box Lib Library

ORION is our magneto-radiation-hydrodynamics AMR code

We Solve: Parallel, 3-D coupled self-gravitating-Radiation-Hydrodynamics on Adaptive Meshes Multi-Scale Physics

I. Hydrodynamics: is solved with conservative, high order, time explicit Godunov scheme with Approximate Riemann Solver

Multi-fluid Hydrodynamics

II. Self-Gravity: We employ parallel, scalable, multi-grid solution algorithms

We use implicit multi-grid iteration to first solve Poisson Equation on a single level

UCRL-PRES-229278- 10

Physics Implementation (cont.)Physics Implementation (cont.)

Level solutions are then coupled and iterated to convergence to obtain solution for gravitational potential on all levels

III. Radiation Transfer: Non-Equilibrium Flux-Limited diffusion including important O(v/c) terms - Radiation solved implicity with parallel multi-grid, iteration scheme taking into account multi-level solves

Solutions must be obtained which couple all grids at a single refinement level, or even across multiple level

IV. Ideal MHD: fully 2nd order unsplit Godunov MHD

We are now implementing this in our AMR self-gravity rad-hydro code

Much better dissipation properties than split staggered mesh schemes (Crockett, Collella, Fisher, Klein & McKee JCP 2004)


HMSF Initial Conditions: Non-TurbulentHMSF Initial Conditions: Non-Turbulent

r –3/2 density profile, r = 0.1–0.2 pc, M = 100–200 Msun, slow solid-body rotation: = 0.02, dynamic range = 8192


Non-Turbulent IC: Early EvolutionNon-Turbulent IC: Early Evolution

At early stages the star accretes steadily and a Keplerian disk forms. Cylindrical symmetry is maintained.


Non-Turbulent IC: Radiation Bubble FormationNon-Turbulent IC: Radiation Bubble Formation

At higher luminosities, radiation pressure forms bubbles above and below the accretion disk. Bubble growth is up-down and cylindrically asymmetric.


Continued Expansion of Radiation BubbleContinued Expansion of Radiation Bubble


High Mass Disk and Formation of Expanding Radiation Driven BubbleHigh Mass Disk and Formation of Expanding Radiation Driven Bubble


Rayleigh-Taylor Instability in Radiation Driven BubbleRayleigh-Taylor Instability in Radiation Driven Bubble


Collapse of radiation driven bubbleCollapse of radiation driven bubble


HMSF: Turbulent Initial ConditionsHMSF: Turbulent Initial Conditions

r –3/2 density profile, Gaussian random velocity field with power on large scales, kinetic energy ~ potential energy Mach number ~8.5, dynamic range = 16,348

QuickTime™ and aYUV420 codec decompressor

are needed to see this picture.


HMSF Protostellar EvolutionHMSF Protostellar Evolution

Turbulent ICs

Non Turbulent ICs


Radius, Accretion Rate and Luminosity of Primary StarRadius, Accretion Rate and Luminosity of Primary Star

start ofDeuterium

burning

accretionluminosity

Principal source of raising temperature in the core is accretion luminosity which is the dominant source of energy prior to nuclear burning


Temperature Distribution in 100 Solar Mass CoreTemperature Distribution in 100 Solar Mass Core

T>100 K, BAR T>300 K, BAR

T>300 K, RTT>100 K, RTT>50 K, RTALL

T>50 K, BAR

Accretion luminosity transported by radiation heats a radius of 1000 AU of the core to > 50K and substantial parts of the core to > 100K


Evolution of 100 Solar Mass Turbulent Protostellar CoreEvolution of 100 Solar Mass Turbulent Protostellar Core

Radiative heating results in the formation of a primary high mass star and 2 low mass stars in the disk


Evolution of 100 Solar Mass Turbulent Isothermal Protostellar CoreEvolution of 100 Solar Mass Turbulent Isothermal Protostellar Core

Isothermal or barotropic models result in the formation of a multitude of low mass stars only erroneous fragmentation


Observing Massive Disks with ALMAObserving Massive Disks with ALMA

Integrated TB in simulated 1000 s / pointing ALMA observation of disk at 0.5 kpc in CH3CN 220.7472 GHz (KKM 2007c, ApJ,)Integrated TB in simulated 1000 s / pointing ALMA observation of disk at 0.5 kpc in CH3CN 220.7472 GHz (KKM 2007c, ApJ,)


Effects of Protostellar OutflowsEffects of Protostellar Outflows

• High mass protostars have outflows that look like larger versions of low mass protostellar outflows (Beuther et al. 2004)

• Outflows are launched inside star’s dust destruction radius

• Due to high outflow velocities, there is no time for dust grains to regrow inside outflow cavities. Grains reach only ~10–3m by the time they escape the core.

• Because grains are small, outflow cavities are optically thin.

• Thin cavities can be very effective at collimating protostellar radiation, reducing the radiation pressure force in the equatorial plane

• Krumholz, McKee & Klein, (2005) using toy Monte-Carlo radiative transfer calculations find outflows cause a factor of 5 – 10 radiation pressure force reduction

• Outflows may be responsible for driving turbulence in clumps (Li & Nakamura 2006)


Protostellar Outflows in High-Mass Star FormationProtostellar Outflows in High-Mass Star Formation

Temperature Distribution Radiation and Gravitational Forces

• Temperature distribution from Monte Carlo diffusion (Whitney, et. al. 2003); Radiation transfer with ray solution to get radiative forces

• Envelope rotationally flattened density dist.; cavity shape Z=ab; M=50M ZAMS; 50M envelope

• With no wind cavity; frad > fgrav everywhere except inside the accretion disk accretion halted

• With wind cavity, frad < fgrav outside disk radius accretion can continue


HMSF with Outflows: Very Early 3-D EvolutionHMSF with Outflows: Very Early 3-D Evolution

Early results show that radiation is collimated effectively by outflow cavities: radiation energy density is factor of ~5 higher inside cavity


Advances Necessary in Algorithmic Performance and Scalability for High Mass Star Formation

Advances Necessary in Algorithmic Performance and Scalability for High Mass Star Formation

• State-of-the-art simulations follow collapse from the scale of turbulent cores to stars (KKM 2007)

Dynamic range > 104

• Simulations will require more realistic initial conditions in core derived from the outer scale imposed by turbulent clumps (M ~ several X 103 M )

Simulations are just beginning to follow collapse

from turbulent Clumps Cores Stars with radiative feedback and AMR

Dynamic Range > 105

• Current state-of-the-art (Krumholz, Klein & McKee 2006, 07) require months to evolve high mass stars on parallel machines (~ 256 processors) with Grey Radiation Transfer multi-frequency will be several times more expensive

• Future simulations will evolve GMCs Clumps Cores Stars

Dynamic Range > 106 - 107

• For galaxy simulations to incorporate star formation

Galaxy GMCs Clumps Cores Stars

Dynamic Range > 3x108 - 3x1010


Summary and Future DirectionsSummary and Future Directions

• 3-D high resolution AMR simulations with ORION achieves protostellar masses considerably above previous 2-D axisymmetric gray simulations

• Two new mechanisms have been shown to overcome radiation pressure barrier to achieve high mass star formation

— 3-D Rayleigh-Taylor instabilities in radiation driven bubbles appear to be important in allowing accretion onto protostellar core

— Protostellar outflows resulting in optically thin cavities promote focusing of radiation and reduction of radiation pressure enhances accretion

— Radiation feedback from accreting protostars inhibits fragmentation

• ALMA observations will help distinguish between competing models of high mass star formation gravitational core collapse predicts large scale disks

Future Directions

• Multi-frequency radiation-hydrodynamics and inclusion of ionization

• Improvement in flux limited diffusion (Monte-Carlo; Sn transport; Variable Edd Tensor)

• Improvement in dust physics (e.g. shattering; coagulation; multi-species)

• Evolution of wind outflow models and interaction with infalling envelope

• Self consistent evolution of high mass turbulent cores from large scale turbulent clump

• Inclusion of MHD can launch hydromagnetic wind; possible photon bubble instab.


Back up SlidesBack up Slides


Radiation Transport Results in Suppression of Large Scale Fragmentation in Massive Star Formation

Radiation Transport Results in Suppression of Large Scale Fragmentation in Massive Star Formation

0.26 pc 6700 AU

Most of the available mass in turbulent cloud goes into one massive star.

QuickTime™ and aYUV420 codec decompressor

are needed to see this picture.


High Mass Disk at 27,000 yr. (Krumholz, Klein & McKee 2007b) High Mass Disk at 27,000 yr. (Krumholz, Klein & McKee 2007b)

0.6 pc radius Disk radius 3000 AU


Observing Massive DisksObserving Massive Disks

Integrated TB in simulated 1000 s / pointing ALMA observation of disk at 0.5 kpc in CH3CN 220.7472 GHz (KKM 2007c, ApJ, in press)Integrated TB in simulated 1000 s / pointing ALMA observation of disk at 0.5 kpc in CH3CN 220.7472 GHz (KKM 2007c, ApJ, in press)


Adaptive Mesh OverviewAdaptive Mesh Overview

• A block-structured refinement strategy combines the advantages of adaptive mesh refinement with the efficiencies provided by uniform grids

• For hyperbolic systems, such as the advection component of fluid dynamics, explicit difference schemes can be used which minimize communication

— A serial algorithm can proceed one grid at a time

— A parallel algorithm can process many grids at once. Library support for this approach is provided by CCSE at LBL.

• For parabolic and elliptic systems, such as those associated with radiation diffusion, implicit difference schemes must be used. Solutions must be obtained which couple all grids at a single refinement level, or even across multiple levels.

— Interactive solvers based on multigrid provide efficient solutions.

— We use the hypre parallel multigrid library developed in CASC for this part of the algorithm

the formation of high mass stars richard i. klein uc berkeley, department of astronomy and lawrence...

Documents

formation of massive

highmass accretion rate

high mass cores shu

formation mechanism

cluster formation

high densities b n

approximation slide

disk accretion