Star Formation1
Kohji Tomisaka
National Astronomical Observatory JapanMay 8, 20071

http://yso.mtk.nao.ac.jp/˜tomisaka/Lecture Notes/StarFormation/3.pdf

cCopyright °2001-2004
by Tomisaka, K.
version 0.5: August 5, 2001
version 1.0: August 30, 2001
version 2.0: November 8, 2003
version 3.0: June 1, 2004

Contents
1 Introduction
1.1 Interstellar Matter . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.2 Case Study — Taurus Molecular Clouds . . . . . . . . . . . . . . . . . . . . . . . . . .
1.3 T Tauri Stars . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.4 Spectral Energy Distribution (SED) . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.5 Protostars . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.5.1 B335 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.5.2 L1551 IRS 5 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.6 L 1544: Pre-protostellar Cores . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.7 Magnetic Fields . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.7.1 Prestellar Core . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.7.2 Cores with Protostars . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.8 Density Distribution . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.9 Mass Spectrum . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.10 Line Width - Size Relation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

1
1
2
4
6
8
8
12
15
15
15
17
22
25
27

2 Physical Background29
2.1 Basic Equations of Hydrodynamics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
2.2 The Poisson Equation of the Self-Gravity . . . . . . . . . . . . . . . . . . . . . . . . . 29
2.3 Free-fall Time . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
2.3.1 Accretion Rate . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
2.4 Gravitational Instability . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
2.4.1 Linear Analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
2.4.2 Sound Wave . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
2.5 Jeans Instability . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
2.6 Gravitational Instability of Thin Disk . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
2.6.1 Rotating Disk . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
2.7 Convective Instability . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
2.8 Super- and Subsonic Flow . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
2.8.1 Flow in the Laval Nozzle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
2.8.2 Steady State Flow under an Influence of External Fields . . . . . . . . . . . . . 41
2.8.3 Stellar Wind — Parker Wind Theory . . . . . . . . . . . . . . . . . . . . . . . 43
2.9 Virial Analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
2.10 Radiative Transfer . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
2.10.1 Radiative Transfer Equation . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
2.10.2 Einstein’s Coefficients . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50i

ii

CONTENTS
2.10.3 Relation of Einstein’s Coefficients to Absorption and Emissivity . . . . . . . . .

51

3 Galactic Scale Star Formation
53
3.1 Schmidt Law . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
3.1.1 Global Star Formation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
3.1.2 Local Star Formation Rate . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
3.2 Gravitational Instability of Rotating Thin Disk . . . . . . . . . . . . . . . . . . . . . . 55
3.2.1 Tightly Wound Spirals . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
3.2.2 Toomre’s Q Value . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
3.3 Spiral Structure . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
3.4 Density Wave Theory . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
3.4.1 Group Velocity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
3.5 Galactic Shock . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
4 Local Star Formation Process69
4.1 Hydrostatic Balance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
4.1.1 Bonnor-Ebert Mass . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
4.1.2 Equilibria of Cylindrical Cloud . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
4.2 Virial Analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
4.2.1 Magnatohydrostatic Clouds . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
4.3 Subcritical Cloud vs Supercritical Cloud . . . . . . . . . . . . . . . . . . . . . . . . . . 75
4.4 Ambipolar Diffusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
4.4.1 Ionization Rate . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
4.4.2 Ambipolar Diffusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
4.5 Dynamical Collapse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
4.5.1 Inside-out Collapse Solution . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
4.5.2 Protostellar Evolution of Supercritical Clouds . . . . . . . . . . . . . . . . . . . 84
4.6 Accretion Rate . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 85
4.7 Outflow . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
4.7.1 Magneto-driven Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
4.7.2 Entrainment Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
4.8 Evolution to Star . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
4.9 Example of Numerical Simulation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 94
4.10 Evolution in the H-R diagram . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
4.10.1 Main Accretion Phase . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
4.10.2 Premain-sequence Evolution . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 102
Bibliography

105

A Basic Equation of Fluid Dynamics
111
A.1 What is fluid? . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 111
A.2 Equation of Motion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 111
A.3 Lagrangian and Euler Equations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
A.4 Continuity Equation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 113
A.4.1 Expression for Momentum Density . . . . . . . . . . . . . . . . . . . . . . . . . 113
A.5 Energy Equation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 114
A.5.1 Polytropic Relation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 114
A.5.2 Energy Equation from the First Law of Theromodynamics . . . . . . . . . . . . 115

CONTENTS

iii

A.6 Shock Wave . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 115
A.6.1 Rankine-Hugoniot Relation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 115
B Basic Equations of Magnetohydrodynamics
117
B.1 Magnetohydrodynamics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 117
B.1.1 Flux Freezing . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 117
B.1.2 Basic Equations of Ideal MHD . . . . . . . . . . . . . . . . . . . . . . . . . . . 118
B.1.3 Axisymmetric Case . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 118
C Hydrostatic Equilibrium
121
C.1 Polytrope . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 121
C.2 Magnetohydrostatic Configuration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 123
D Basic Equations for Radiative Hydrodynamics
125
D.1 Radiative Hydrodynamics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 125
E Random Velocity

127

Chapter 1

Introduction
1.1

Interstellar Matter
8.0
7.0"Coronal" Gas

Stellar Wind

6.0

Supernovalog T (K)

5.0

2.0

Dissociation of H2

HII Regions

4.0
3.0

Star

Interaction with Supernova
Second Collapse

Mass Loss from Red GiantIntercloud Gas

Photo-dissociation

Diffuse Cloud

1.0
0.0
-4.0

Warm or Hot Core

Quasi-static Contraction

Isothermal Collapse

Molecular Cloud
-2.0

0.0

2.0

First Stellar Core

Globule
4.0

6.0

8.0

10.0 12.0
log n (cm-3)

14.0

16.0

Second Stellar Core
18.0

20.0

22.0

Figure 1.1: Multiphases of the interstellar medium. The temperature and number density of gaseous
objects of the interstellar medium in our Galaxy are summarized. Originally made by Myers (1978),reconstructed by Saigo (2000).
Figure 1.1 shows the temperature and number density of gaseous objects in our Galaxy. Coldinterstellar medium forms molecular clouds (T ∼ 10K) and diffuse clouds (T ∼ 100K). Warm inter4
stellar medium 103 K<∼T <∼10 K are thought to be pervasive (wide-spread). HII regions are ionizedby the Ly continuum photons from the nearby early-type stars. There are coronal (hot but tenuous)gases with T ∼ 106 K in the Galaxy, which are heated by the shock fronts of supernova remnants.
4−3
Pressures of these gases are in the range of 102 K cm−3 <∼p<∼10 K cm , except for the HII regions.
This may suggest that gases are in the pressure equilibrium.
1

2

CHAPTER 1. INTRODUCTION

In this figure, a theoretical path from the molecular cloud core to the star is also shown. We will
see the evolution more closely in Chapter 4.
Globally, the molecular form of Hydrogen H2 is much abundant inside the Solar circle, while theatomic hydrogen HI is more abundant than molecular H2 in the outer Galaxy. In Figure 1.2 (left),the radial distributions of molecular and atomic gases are shown. The right panel shows similar
distributions for four typical external galaxies (M51, M101, NGC6946, and IC342). This indicatesthese distributions are similar with each other. HI is distributed uniformly, while H2 density increases
greatly reaching the galaxy center. In other words, only in the region where the total (HI+H2 ) density
exceeds some critical value, H2 molecules are distributed.

Figure 1.2: Radial distribution of H2 (solid line) and HI (dashed line) gas density. (Left:) our Galaxy.
Converting from CO antenna temperature to H2 column deity, n(CO)/n(H2 ) = 6 × 10−5 is assumed.
Taken from Gordon & Burton (1976). (Right:) Radial distribution of H2 and HI gas for external
galaxies. The conversion factor is assumed constant X(H2 /CO) = 3 × 1020 H2 /K km s−1 . Taken from
Honma et al (1995).

1.2

Case Study — Taurus Molecular Clouds

Figure 1.3 (left) shows the 13 CO total column density map of the Taurus molecular cloud (Mizuno etal 1995) whose distance is 140 pc far from the Sun. Since 13 CO contains 13 C, a rare isotope of C, theabundance of 13 CO is much smaller than that of 12 CO. Owing to the low abundance, the emission
lines of 13 CO are relatively optically thiner than that of 12 CO. Using 13 CO line, we can see deep
inside of the molecular cloud. The distributions of T Tauri stars and 13 CO column density coincidewith each other. Since T Tauri stars are young pre-main-sequence stars with M ∼ 1M¯ , which are in
the Kelvin-Helmholtz contraction stage and do not reach the main-sequence Hydrogen burning stage,
it is shown that stars are newly formed in molecular clouds.

1.2. CASE STUDY — TAURUS MOLECULAR CLOUDS

3

Figure 1.3: (Left) 13 CO total column density map of the Taurus molecular cloud (Mizuno et al 1995).
Taken from their home page with URL of http://www.a.phys.nagoya-u.ac.jp/nanten/taurus.html (inJapanese). T Tauri stars, which are thought to be pre-main-sequence stars in the Kelvin-Helmholtz
contraction stage, are indicated by bright spots. (Right) C18 O map of Heiles cloud 2 region in the
Taurus molecular cloud (Onishi et al. 1996). This shows clearly that the cloud is composed of anumber of high-density regions.
Since 18 O is much more rare isotope (18 O/16 O ¿ 13 C/12 C), the distribution of much higherdensity gases is explored using C18 O lines. Figure 1.3 (right) shows C18 O map of Heiles cloud 2 in
the Taurus molecular cloud by Onishi et al (1996). This shows us that there are many molecular
cloud cores which have much higher density than the average. Many of these molecular cloud cores
are associated with IRAS sources and T Tauri stars. It is shown that star formation occurs in themolecular cloud cores in the molecular cloud. They found 40 such cores in the Taurus molecular cloud.
Typical size of the core is ∼ 0.1 pc and the average density of the core is as large as ∼ 104 cm−3 . The
mass of the C18 O cores is estimated as ∼ 1 − 80M¯ .
H13 CO+ ions are excited only after the density is much higher than the density at which CO
molecules are excited. H13 CO+ ions are used to explore the region with higher density than thatobserved by C18 O. Figure 1.4 shows the map of cores observed by H13 CO+ ions. The cores shown in
the lower panels are accompanied with infrared sources. The energy source of the stellar IR radiationis thought to be maintained by the accretion energy. That is, since the gravitational potential energy
at the surface of a protostar with a radius R∗ and a mass M∗ is equal to Φ ' −GM∗ /R∗ , the kinetic
energy of the gas accreting on the stellar surface is approximately equal to ∼ GM∗ /R∗ . The energy
inflow rate owing to the accretion is (∼ GM∗ /R∗ ) × Ṁ ∼ (GM∗ /R∗ ) × A(c3s /G), where Ṁ = A(c3s /G)is the mass accretion rate. In the upper panels, the cores without IR sources are shown. This coredoes not show accretion but collapse. That is, before a protostar is formed, the core itself contract
owing to the gravity, which is explained closely in chapter 4.
In Figure 1.4, H13 CO+ total column density maps of the C18 O cores are shown. Cores in thelower panels have associated IRAS sources, while the cores in the upper panels have no IRAS sources.
Since the IRAS sources are thought to be protostars or objects in later stage, the core seems to evolvefrom that without an IRAS sources to that with an IRAS source. From this, the core with an IRAS
source is called protostellar core, which means that the cores contain protostars. On the other
hand, the core without IRAS source is called pre-protostellar core or, in short, pre-stellar core.
Figure 1.4 shows that the prestellar cores are less dense and more extended than the protostellar
core. This seems to suggest the density distribution around the density peak changes between before

4

CHAPTER 1. INTRODUCTION

Figure 1.4: Pre-protostellar vs protostellar cores (H13 CO+ map). Upper panel shows the C18 O coreswithout associated IRAS sources. Lower panel shows the cores with IRAS sources. Taken from Fig.1
of Mizuno et al. (1994)
and after the protostar formation.

1.3

T Tauri Stars

T Tauri stars are observationally late-type stars with strong emission lines and irregular light variations associated with dark or bright nebulosities. T Tauri stars are thought to be low-mass pre-mainsequence stars, which are younger than the main-sequence stars. Since these stars are connecting
between protostars and main-sequence stars, they attract attention today. More massive counterparts are called Herbig Ae-Be stars. They are doing the Kelvin-Helmholtz contraction in which theown gravitational energies released as it contracts gradually and this is the energy source of the luminosity. Many emission lines are found in the spectra of T Tauri stars. WTTS (Weak emission T Tauri
Stars or Weak line T Tauri Stars) and CTTS (Classical T Tauri Stars) are classified by their equivalent widths (EWs) of emission lines. That is, the objects with an EW of Hα emission < 1nm = 10Å is
usually termed a WTTS. Figure 1.5 is the HR diagram (Teff − Lbol ) of T Tauri stars in Taurus-Auriga
region (Kenyon & Hartman 1995). WTTSs distribute near the main-sequence and CTTSs are foundeven far from the main-sequence. A number of theoretical evolutional tracks for pre-main-sequencestars with M ∼ 0.1M¯ − 2.5M¯ are shown in a solid line, while the isochorones for ages of 105 yr,
106 yr, and 107 yr are plotted in a dashed line. Vertical evolutional paths are the Hayashi convective
track. Since D=2 H has a much lower critical temperature (and density) for a fusion nuclear reaction
to make He than 1 H, Deuterium begins to burn before reaching the zero-age-main sequence. Thisoccurs near the isochrone for the age of 105 yr and some activities related to the ignition of Deuterium
seem to make the central star visible (Stahler 1983).
Disk Frequency
Infrared studies of T Tauri stars in star-forming regions have suggected that initial disk frequency
if rather high and that the disk lifetimes are relatively short 3 − 15Myr. From JHKL photometry,
Haish, Lada, & Lada (2001) obtained the fraction of disk-bearing stars for 6 star firmation regions.
L band excess emission indicates an accompanied disk. The fraction is a decreasing function of age

1.3. T TAURI STARS

5

Figure 1.5: HR diagram of T Tauri stars. Many emission lines are found in the spectra of T Tauristars. WTTS (Weak Emission T Tauri Stars) and CTTS (Classical T Tauri Stars) are classified by
the equivalent widths of emission lines. That is, the objects with an EW of Hα emission < 1nm is
usually termed a WTTS. Taken from Fig.1.2 of Hartmann (1998).

6

CHAPTER 1. INTRODUCTION

Figure 1.6: Fractions of IR excess sources in respective clusters are plotted againt the age of theclusters.
of the cluster as in Figure 1.6. This figure shows clearly that the disk fraction is initially very high(>∼80%) and rapidly decreases with increasing cluster age. In 3 Myr a half of the disk stars lose theirdisks. Overall disk lifetime is estimated as ∼ 6Myr.

1.4

Spectral Energy Distribution (SED)

A tool to know the process of star formation is provided by the spectral energy distribution (SED)
mainly in the near- and mid-infrared light. T Tauri stars and protostars have typical respective SEDs.
IR SEDs of T Tauri stars were classified into three as Class I, Class II, and Class III, from a standpoint of relative importance of the radiation from a dust disk to the stellar black-body radiation.
Today, the classification is extended to the protostars, which is precedence of the T Tauri stars, andthey are called Class 0 objects. (Unfortunately, there is no zero in Roman numerals.) In Figure 1.7,
typical SEDs and models of emission regions are shown.
1. Class III is well fitted by a black-body spectrum, which shows the energy mainly comes froma central star. This SED is observed in the weak-line T Tauri stars. Although T Tauri stars
show emission lines of such as the Hydrogen Balmar sequence, the weak-line T Tauri stars do
not show prominent emission lines, which indicates the amount of gas just outside the star
(this seems to be supplied by the accretion process) is small. In this stage, a disk has been
disappeared or an extremely less-massive disk is still alive.
2. Class II SED is fitted by a single-temperature black-body plus excess IR emission. This showsthat there is a dust disk around a pre-main sequence star and it is heated by the radiationfrom a central star. The width of the spectrum of the disk component is much wider than that
expected from a single-temperature black-body radiation. Thus, the disk has a temperature
gradient which decreases with increasing the distance from the central star. In this stage, thedust disk is more massive than that of Class III sources. Classical T Tauri stars have such
SEDs.

1.4. SPECTRAL ENERGY DISTRIBUTION (SED)

7

Figure 1.7: Spectral Energy Distribution (SED) of young stellar objects (YSOs) and their models.
(Left:) ν − νFν plot taken from Lada (1999). (Right:) λ − λFλ plot taken from André (1994)
3. In Class I SED, the mid infrared radiation which seems to come from the dust envelope ispredominant over the stellar black-body radiation. Since the stellar black-body radiation seems
to escape at least partially from the dust envelope, a relatively large solid angle is expected fora region where the dust envelope does not intervene.
4. Class 0 SED seems to be emitted by an isothermal dust with ∼ 30K. The protostar seems to
be completely covered by gas and dust and is obscured with a large optical depth by the dustenvelope. No contribution can be reached from the stellar-black body radiation.
The reason why the emission from the disk becomes wide in the spectral range is understood(Fig.1.8) as follows: Temperature of the disk is determined by a balance of heating and cooling.
Assuming the disk is geometrically thin but optically thick, the cooling per unit area is given by
the equation of the black-body Planck radiation. Therefore, the temperature is determined by theheating predominantly by viscous heating and extra heating by the radiation from the central star.

8

CHAPTER 1. INTRODUCTION

Figure 1.8: Explanation for the spectral index of the emission from a geometrically thin but opticallythick disk. Taken from Fig.16 of Lada (1999).
The flux density emitted by the disk is given by
Z

νFν ∼

νπBν [T (r)]2πrdr ∼ r(T or ν)2 νBν .

(1.1)

Assuming the radial distribution of temperature asµ

T = T0

R
R0

¶−q

,

(1.2)

(q = 3/4 for the standard accretion disk) and taking notice that each temperature in the disk radiates
at a characteristic frequency ν ∝ T (the Wien’s law for black-body radiation)
νFν ∼ r2 νBν ∝ ν 4 T −2/q ∝ ν 4−2/q ,(1.3)

where we used the fact that the peak value of Bν ∝ ν 3 . Therefore, it is shown thatνFν ∝ ν n ;2
n=4− .
q(1.4)

As shown in the previous section, we have no young stellar objects found by IR before a protostar
is formed. These kind of objects (pre-protostellar core) are often called Class −1. The classification
was originally based on the SED and did not exactly mean an evolution sequence. However, today
YSOs are considered to evolve as the sequence of the classes: Class −1 → Class 0 → Class I → Class
II → Class III → main-sequence star.

1.5

Protostars

1.5.1

B335

B335 is a dark cloud (Fig.1.9) with a distance of D ' 250pc. Inside the dark cloud, a Class 0 IR
source is found. The object is famous for the discovery of gas infall motion. In Figure 1.10, the1.5. PROTOSTARS

9

Figure 1.9: Near infrared images of B335, which is Class 0 source.

Figure 1.10: Line profile of CS J = 2 − 1 line radio emission. Model spectra illustrated in a dashedline (Zhou 1995) are overlaid on to the observed spectra in a solid line (Zhoug et al 1993).

10

CHAPTER 1. INTRODUCTION

Figure 1.11: Explanation of blue-red asymmetry when we observe a spherical symmetric inflow motion. An isovelocity curve for the red-shifted gas is plotted in a solid line. That for the blue-shifted
gas is plotted in a dashed line. Taken from Fig.14 of Lada (1999).
line profiles of CS J = 2 − 1 line emissions are shown (Zhou et al. 1993). The relative position of
the profiles correspond to the position of the beam. (9,9) represents the offset of (9”,9”) from thecenter. At the center (0,0), the spectrum shows two peaks and the blue-shifted peak is brighter than
the red-shifted one. This is believed to be a sign of gas infall motion. The blue-red asymmetry is
explained as follows:
1. Considering a spherical symmetric inflow of gas, whose inflow velocity vr increases with reaching
the center (a decreasing function of r)
2. Considering a gas element at r moving at a speed of vr (r) < 0, the velocity projected on aline-of-sight is equal tovline−of−sight = vsystemic + vr cos θ,(1.5)
where vsystemic represents the systemic velocity of the cloud (line-of-sight velocity of the cloudcenter) and θ is the angle between the line-of-sight and the position vector of the gas element.
The isovector lines, the line which connect the positions whose procession/recession velocities
are the same, become like an ellipse shown in Fig.1.11.
3. An isovelocity curve for the red-shifted gas is plotted in a solid line. That for the blue-shifted
gas is plotted in a dashed line. If the gas is optically thin, the blue-shifted and red-shifted gases
contribute equally to the observed spectrum and the blue- and red-shifted peaks of the emission
line should be the same.
4. In the case that the gas has a finite optical depth, for the red-shifted emission line a cold gas in
the fore side absorbs effectively the emission coming from the hot interior. On the other hand,for the blue-shifted emission line, the emission made by the hot interior gas escapes from thecloud without absorbed by the cold gas (there is no cold blue-shifted gas).

1.5. PROTOSTARS

11

Figure 1.12: C18 O total column density map (left) and H13 CO+ channel map (right) of B335 alongwith the position-velocity maps along the major and minor axes. Taken from Fig.3 of Saito et al(1999).
5. As a result, the blue peak of the emission line becomes more prominent than that of the readshifted emission. This is the explanation of the blue-red asymmetry.
In Figure 1.10, model spectra calculated with the Sobolev approximation (Zhou 1995) are shown.
These show the blue-red asymmetry (the blue line > the red line).
Early observation of star forming regions have revieled fast molecular outflow are often ejected−1from the central protostar with >∼10km s . B335 is also a typical outflow source. In Figure 1.12,
distributions of high density gases traced by the C18 O and H13 CO+ lines are shown as well as the
bipolar outflow whose outline is indicated by a shadow (Hirano et al 1988). Comparing left and rightpanels, it is shown that the distribution of C18 O gas is more extended than that of H13 CO+ which
traces higher-density gas. And the distribution of the H13 CO+ is more compact and the projectedsurface density seems to show the the actual distribution is spherical. And the molecular outflow
seems to be ejected in the direction of the minor axis of the high-density gas. It may suggest that (1) amolecular outflow is focused or collimated by the effect of density distribution or that (2) collimation
is made by the magnetic fields which run preferentially perpendicularly to the gas disk. This gas disk
is observed by these high density tracers.
Combining the C18 O and H13 CO+ distributions, the surface density distribution along the major
axis is obtained by Saito et al (1999). From the lower panel of Figure 1.13, the column density
distribution is well fitted in the range from 7,000 to 42,000 AU in radius,
µ
21

−2

Σ(r) = 6.3 × 10 cm

r
4
10 AU

¶−0.95,(1.6)

where they omitted the data of r<∼7000 since the beam size is not negligible. Similar power-law
density distributions are found by the far IR thermal dust emission.

12

CHAPTER 1. INTRODUCTION

Figure 1.13: Column density distribution NH (r) derived from the H13 CO+ and C18 O data taken by
the Nobeyama 45 m telescope. Taken from Fig.9 of Saito et al (1999).

1.5.2

L1551 IRS 5

Figure 1.14: (Left:) 13 CO column density distribution. The contour lines represent the distributionof 13 CO column density. 2.2 µm infra-red reflection nebula is shown in grey scale which was observedby Hoddap (1994). (Right:) Schematic view of L1551 IRS5 region.
L1551 IRS 5 is one of the most well studied protostellar objects. This has an infra-red emission
nebulosity (Fig.1.14). It is believed that there is a hole perpendicular to the high-density disk and theemission from the central star escapes through the hole and irradiate the nebulosity. In this sense thisis a reflection nebula. L1551 IRS 5 has an elongated structure of dense gas similar to that observedin B335. The gas is extending in the direction from north-west to south-east [Fig.1.14 (left)]. Since
the opposite side of the nebulosity is not observed, the opposite side of nebulosity seems to be locatedbeyond the high-density disk and be obscured by the disk. This is possible if we see the south surfaceof the high-density disk as in Figure 1.14 (right).

1.5. PROTOSTARS

13

Figure 1.15: Isovelocity contours measured by the 13 CO J = 1 − 0 line. It should be noticed that theisovelocity lines run parallelly to the major axis. The north-eastern side shows a red-shift and thesouth-western side shows a blue-shift.
Infall Motion
The inflow motion is measured. Figure 1.15 shows the isovelocity contours measured by the 13 CO
J = 1−0 observation (Ohashi et al 1996). It should be noticed that the isovelocity lines run parallelly
to the major axis. The north-eastern side shows a red-shift and the south-western side shows a blueshift. Considering the configuration of the gas disk shown in Fig.1.14 (right), this pattern of isovelocity
contours indicates not outflow but inflow. That is, the north-east side is a near side of the disk and thesouth-west side is a far side. Since a red-shifted motion is observed in the near side and a blue-shiftedmotion is observed in the far side, it should be concluded that the gas disk of the L1551 IRS5 is now
infalling.
Optical Jet
HST found two optical jets emanating from L1551 IRS5. This has been observed by SUBARU
telescope, which found jet emission is dominated by [FeII] lines in the J- and H-bands. The jet
extents to the south-western direction and disappears at ∼ 10” ' 1400AU from the IRS5. The
width-to-length ratio is very small <∼1/10 or less, while the bipolar molecular outflow shows a lesscollimated flow. As for the origin of the two jets, these two jets might be ejected from a single source.
However, since there are at least two radio continuum sources in IRS5 within the mutual separationof ∼ 0.”5 [see Fig.1.16 (right)], these jets seem to be ejected from the two sources independently.

14

CHAPTER 1. INTRODUCTION

Figure 1.16: (Left:) Infrared image (J- and K-band) of the IR reflection nebula around L1551 IRS5
by SUBARU telescope. Taken from Fig.1 of Itoh et al. (2000). (A jpeg file is available from thefollowing url: http://SubaruTelescope.org/Science/press release/9908/L1551.jpg). (Right:)
Central 100 AU region map of L1551 IRS5. This is taken by the λ = 2.7cm radio continuum
observation. Deconvolved map (lower-left) shows clearly that IRS5 consists of two sources. Taken
from Looney et al. (1997).

1.6. L 1544: PRE-PROTOSTELLAR CORES

15

Although the lengths of these jets are restricted to 10”, Herbig-Haro jets, which are much larger
than the jets in L1551 IRS5, have been found. HH30 has a ∼ 500 AU-scale jet whose emission ismainly from the shock-excited emission lines. One of the largest ones is HH111, which is a member of
the Orion star forming region and whose distance is as large as D ∼ 400pc, and a jet with a length of∼ 4pc is observed. Source of HH111 system is thought to consist of at least binary stars or possibly
triple stars [Reipurth et al (1999)]. Star A, which coincides with a λ = 3.6 cm radio continuum
source (VLA 1), shows an elongation in the VLA map whose direction is parallel to the axis of thejet. Therefore, star A is considered to be a source of the jet. Since the VLA map of star A shows
another elongated structure perpendicular to the jet axis, star A may be a binary composed by two
outflow sources.

1.6

L 1544: Pre-protostellar Cores

L1544 is known as a pre-protostellar core (Taffala et al 1998). That shows an infall motion but thiscontains no IR protostars. In Figure 1.18(left), CCS total column density map is shown, which shows
an elongated structure. Ohashi et al (1999) have found both rotation and infall motion in the cloud.
Position-velocity (PV) diagram along the minor axis shows the infall motion. That along the major
axis indicates a rotational motion, which is shown by a velocity gradient. Due to a finite size of thebeam, a contraction motion is also shown in the PV diagram along the major axis.

1.7

Magnetic Fields

Directions of magnetic field are studied by (1) measuring the polarization of light which is sufferedfrom interstellar absorption. In this case the direction of magnetic field is parallel to the polarization
vector. The reason is explained in Figure 1.19. In the magnetic fields, the dusts are aligned in a
way that their major axes are perpendicular to the magnetic field lines. Such aligned dusts absorb
selectively the radiation whose E-vector is parallel to their major axes. As a result, the detected light
has a polarization parallel to the magnetic field lines.
However, the polarization measurement in the near IR wavelength limited to the region withlow gas density, because background stars suffer severe absorption and becomes hard to be observed
if we want to measure the polarization of the high-density region. More direct method is (2) the
measurement of the linear polarization of the thermal emission from dusts in the mm (or sub-mm)
wavelengths; in this case the direction of magnetic field is perpendicular to the polarization vector.
The mechanism is explained in Figure 1.19b. The aligned dusts, whose major axes are perpendicular
to the magnetic field lines, emit the radiation whose E-vector is parallel to the major axes. Since theabsorption does not have a severe effect in this mm wavelengths, this gives information the magnetic
fields deep inside the clouds.

1.7.1

Prestellar Core

Figure 1.20 illustrates the polarization maps of three prestellar cores (L1544, L183, and L43) donein the 850 µm band by JCMT-SCUBA. In L1544 and L183 the mean magnetic fields are at an angleof 30 deg to the minor axes of the cores. L43 is not a simple object; there is a T Tauri star locatedin the second core which extends to south-western side of the core (an edge of this core is seen near
the western SCUBA frame boundary). And a molecular outflow from the source seems to affectthe core. The magnetic field as well as the gas are swept by the molecular outflow. L43 seems an
exception. The fact that the mean magnetic fields are parallel to the minor axis of the high-density

16

CHAPTER 1. INTRODUCTION

Figure 1.17: A mosaic image of HH 111 based on HST NICMOS images (bottom) and WFPC2 images(top). Taken from Fig.1 of Reipurth et al (1999).

1.7. MAGNETIC FIELDS

17

Figure 1.18: CCS image of prestellar core L1544. (Left:) Total intensity map. (Right:) PV diagrams
along the minor axis (left) and along the major axis (right).

gas distribution seems to mean gas contracts preferentially in the direction parallel to the magnetic
fields.

1.7.2

Cores with Protostars

Bok Globules
Bok globules are isolated dark clouds. Wolf, Launhardt & Henning (2003) have studied the relation
between the magnetic field and outflow directions. Figure 1.21 displays the direction of the magnetic
field obtained with JCMT-SCUBA polarization observation in 850µm as well as the outflow found by
13CO observation of Chandler & Sargent (1993) and 12 CO observation by Yun & Clemens (1994).
In B 335, CB 230 and CB 244 outflows are oriented in the direction perpendicular to the major
axis of the globules. The magnetic field is running parallel to the outflow but perpendicular to themajor axis. This means that a disk is formed with a gas flowing along the magnetic field. Further, theoutflow is generated by a twisted magnetic field which is an outcome of a rotating disk. In globule CB
54, the alignment is not perfect, that is, magnetic field direction is slightly aligned with the outflow(Wolf, Launhardt & Henning 2003). This seems to stengthen the magnetic origin of the outflow andmagnetically guided disk formation.
It is clearly shown that the polarization anti-correlates with the intensity of the thermal emission.
This might be due to the observational error in the low intensity region (or low S/N region). Thismay be related to physical processes such as (1) the alingnment owing to the magnetic field becomesinefficient in high density region or (2) the magnetic field is tangled in the dense region and thepolarization due to the aligned dust is canceled.

18

CHAPTER 1. INTRODUCTION

Figure 1.19: Explanation how the polarized radiation forms. Taken from Weintraub et al.(2000).

1.7. MAGNETIC FIELDS

19

Figure 1.20: Directions of B-Field are shown from the linear polarization observation of 850 µm
thermal emission from dusts by JCMT-SCUBA. L 1544 and L183, the magnetic field and the minor
axis of the molecular gas distribution coincide with each other within ∼ 30deg. Taken from WardThompson et al (2000).

20

CHAPTER 1. INTRODUCTION

Figure 1.21: Directions of magnetic field are shown from the linear polarization observation of 850µm thermal emission from dusts by JCMT-SCUBA. These four objects are known as Bok globules.
T Tauri Disks
Magnetic field at the position of protostars and T Tauri stars are measured for IRAS 16293-2422,
L1551 IRS5, NGC1333 IRAS 4A, and HL Tau (Tamura et al. 1995). Although HL Tau is a T Tauri
star, it has a gas disk. Thus this is a Class I source. The others are believed to be in protostellar
phase (Class 0 sources). It is known that IRAS 16293-2422, L1551 IRS5, and HL Tau have diskswith the radii of 1500-4000 AU from radio observations of molecular lines. Further, near infraredobservations have shown that these objects have 300-1000 AU dust disks. Figure 1.22 shows the
E-vector of polarized light. If this is the dust thermal radiation, the direction of magnetic fields isperpendicular to the polarization E-vector. Figure shows the magnetic fields run almost perpendicular
to the elongation of the gas disk. Global directions of magnetic field outside the gas disk and the
direction of CO outflows are also shown in the figure by arrows. It is noteworthy that the directionsof local magnetic fields, global magnetic fields, and outflows coincide with each other within ∼ 30
deg.

1.7. MAGNETIC FIELDS

21

Figure 1.22: Polarization of the radio continuum λ = 1mm, λ = 0.8mm. IRAS 16293-2422 (upperleft), L1551 IRS5 (upper-right), NGC1333 IRAS 4A (lower-left), and HL Tau (lower-right). Taken
from Tamura et al (1995).

22

1.8

CHAPTER 1. INTRODUCTION

Density Distribution

Motte & André (2001) made λ =1.3 mm continuum mapping survey of the embedded young stellar
objects (YSOs) in the Taurus molecular cloud. Their maps include several isolated Bok globules, aswell as protostellar objects in the Perseus cluster. For the protostellar envelopes mapped in Taurus,
the results are roughly consistent with the predictions of the self-similar inside-out collapse model of
Shu and collaborators (section 4.5.1). The envelopes observed in Bok globules are also qualitatively
consistent with these predictions, providing the effects of magnetic pressure are included in the model.
By contrast, the envelopes of Class 0 protostars in Perseus have finite radii <∼10000 AU and are afactor of 3 to 12 denser than is predicted by the standard model.
Another method to measure the density distribution is to use the near IR extinction. From
(H − K) colors of background stars, the local value of AV in a dark cloud can be obtained using astandard reddening law,
AV = 15.87E(H − K)(1.7)
if the intrinsic colors of background stars are known. Here, the color excess is defined as the difference
between the observed color and the intrinsic color: E(H − K) ≡ (H − K)obs − (H − K)intrinsic . We
can convert the extinction to the column density assuming the gas/dust ratio is constant
N (H + H2 ) = 2 × 1021 cm−2 mag−1 AV .

(1.8)

This is a standard method to obtain the local column density of the dark cloud using the near IR
photometry.
See Figure 1.23. If the density distribution is expressed asµρ(r) = ρ0

r
r0

¶−α

,

(1.9)

where r is a physical distance from the center. The column density distribution against the projecteddistance of the line-of-sight from the center of the cloud is given asNρ (p) = 2

Z (R2 −p2 )1/2
0hi

ρ (s2 + p2 )1/2 ds,(1.10)

where R represents the outer radius of the cloud. Using equation (1.8), this yields AV distribution
AV (p) = 10−23 ρ0 r0α

Z (R2 −p2 )1/2
0

(s2 + p2 )−α/2 ds.

(1.11)

If background stars are uniformly distributed, the number of stars with AV |obs is proportional to thearea which satisfies AV |obs = AV (p). That is, if we plot AV (p) against 2πpdp, this gives the number
distribution of background stars with AV . Figure 1.25 shows the result of L977 dark cloud by Alveset al (1998).
Recently, Alves et al (2001) derived directly the radial distribution of NH by comparing the NH (p)model distribution for B68. They obtained a distribution is well fitted by the Bonner-Ebert sphere inwhich a hydrostatic balance between the self-gravity and the pressure force is achieved (lower panelsof Fig.1.25) (see section 4.1.1).
In this fields, we should pay attention to the density distribution in cylindrical clouds. As seenin the Taurus molecular cloud, there are may filamentary structures in a molecular cloud. In §4.1,
we will give the distribution for a hydrostatic spherically symmetric and that of a cylindrical cloud.
The former is proportional to ρ ∝ r−2 and the latter is ρ ∝ r−4 . Therefore, the distribution ρ ∝ r−4

1.8. DENSITY DISTRIBUTION

23

Figure 1.23: Schematic view to explain an AV distribution.

Figure 1.24: (Left:) Radial intensity profiles of the environment of 7 embedded YSOs (a-g) and 1starless core (h). (Right:) Same as left panel but for 4 isolated globules (a-d) and 4 Perseus protostars(e-h). Taken from Motte & Andre (2001).

24

CHAPTER 1. INTRODUCTION

Figure 1.25: Density distribution of L977 (top) and B68 (bottom) dark clouds. (Top-left:) L977
dark cloud dust extinction map derived from the infrared (H-K) observations. (Top-right:) Observed
frequency distribution of extinction measurements for L977 and the predictions from clouds models
with density structures ρ(r) ∝ r−α having α = 1 (dashed line), 2 (solid line), 3 (dotted line), and4 (dash-dotted line). (Bottom-left:) B68 images (false color images made from B, V, and I images

1.9. MASS SPECTRUM

25

Figure 1.26: A structure of magnetic fields in the L1641 region. Polarization of light from embedded
stars (Vrba et al. 1988) is shown by a bar. The direction of magnetic fields in the line-of-sight isobserved using the HI Zeeman splitting, which is shown by a circle and cross (Heiles 1989).
was expected for cylindrical cloud. From near IR extinctions observation (Alves et al 1998), even if
a cloud is rather elongated [Fig.1.25 (top-left)], the power of the density distribution is equal to not-4 but ' −2. Fiege, & Pudritz (2001) proposed an idea that a toroidal component of the magnetic
field, Bφ , plays an important role in the hydrostatic balance of the cylindrical cloud (Fig.1.26).

1.9

Mass Spectrum

We have seen that a molecular cloud consists in many molecular cloud cores. For many years, there
are attempts to determine the mass spectrum of the cores.
From a radio molecular line survey, a mass of each cloud core is determined. Plotting a histogram
number of cores against the mass, we have found that a mass spectrum can be fitted by a power law
as
dN
= Mn
(1.12)
dM
where dN/dM represents the number of cores per unit mass interval. Many observation indicate thatn ∼ −1.5 as Table 1.1.
Figure 1.27 (Motte et al 2001) shows the cumulative mass spectrum (N (> M ) vs. M ) of the 70
starless condensations identified in NGC 2068/2071. The mass spectrum for the 30 condensations of
the NGC 2068 sub-region is very similar in shape. The best-fit power-law is N (> M ) ∝ M n+1 ∝M −1.1 above M >∼0.8M¯ . That is, n = −2.1. This power derived from the dust thermal emission isdifferent from that derived by the radio molecular emission lines. The power n + 1 = −1.1 which
is close to the Salpeter’s IMF for new born stars, N (> M∗ ) ∝ M∗−1.35 might be meaningful. The
reason of the difference is not clear. For example, Tothill et al. (2002) observed the Lagoon nebula

26

CHAPTER 1. INTRODUCTION

Table 1.1: Mass spectrum indicies derived with molecular line surveys.
Papern
RegionMass range
Loren (1989)
−1.1
ρ Oph
10M¯ <∼M <∼300M¯
3<
Stutzki & Guesten (1990) −1.7 ± 0.15 M17 SW
a few M¯ ∼M <∼a few 10 M¯
Lada et al (1991)
−1.6
L1630
M>
∼20M¯
Nozawa et al (1991)−1.7
ρ Oph North
3M¯ <
∼M <
∼160M¯Tatematsu et al. (1993)−1.6 ± 0.3
Orion AM>
50M¯∼
Dobashi et al (1996)−1.6
CygnusM>∼100M¯Onish et al (1996)−0.9 ± 0.2
Taurus
3M¯ <
∼M <
∼80M¯
4
∗
−4
Kramer et al.(1998)−1.6 ∼ −1.8 L1457 etc
10 M¯ <∼M <∼10 M¯Heithausen et al (2000)−1.84
MCLD123.5+24.9,Polaris Flare MJ <
∼M <
∼10M¯
13
+
Onishi et al. (2002)−2.5
Taurus H CO
3.5M¯ <
∼20.1M¯∼M <∗ MCLD126.6+24.5, NGC 1499 SW, Orion B South, S140, M17 SW, and NGC 7538

Figure 1.27: Cumulative mass distribution of the 70 pre-stellar condensations of NGC 2068/2071. The
dotted and dashed lines are power-laws corresponding to the mass spectrum of CO clumps (Kramer
et al. 1996) and to the IMF of Salpeter (1955), respectively. Taken from Fig.3 of Motte et al (2001).

1.10. LINE WIDTH - SIZE RELATION

27

around the edge of the HII region M8. From the continum emission λ = 450µm, 850µm, 1.3mm,
they obtained index of −1.7 ± 0.6, which is consistent with other molecular line studies.

1.10

Line Width - Size Relation

Larson (1981) compiled the observations for molecular cloud complexes, molecular cloud and molecular clumps published in 1974-1979. He obtained an empirical relation that the size of a structure iswell correlated to the random velocity in the structure which is measured by the width of the emission
line (see Appendix E). Figure 1.28(left) shows this correlation and this is well expressed asµσ ' 1.10km s−1

L
1pc

¶0.38

,

(1.13)

where σ and L represent respectively the three-dimensional random speed of gas and the size of thestructure. A similar correlation is found only for giant molecular clouds (Sanders, Scoville, & Solomon1985) asµ¶1/2µ¶
L 0.62
3−1∆V=(0.535±0.16)kms,(1.14)σ=
FWHM
23 ln 2
1pc
for GMCs whose sizes are larger than 10pc (be careful the typos in their abstract: power was 0.62). He also found another correlation between the mass M and the random velocity like Figure1.28(right), which is expressed asµ−1

σ ' 0.43km s

M
1M¯

¶0.20

.

(1.15)In the next chapter (§2.9), we will see the virial relation, that is, for an isolated system to achievea mechanical equlibrium the gravitaional to thermal energy ratio has to be equal to 2:1 for γ = 5/3
gas. The ratio of the gravitational energy ∼ (3/5)GM 2 /(L/2) to the thermal energy M σ 2 /2 is also
fitted asµ¶
12 GM
L 0.14' 1.1,(1.16)
5 σ2L
1pc
which is weakly dependent of the size or the mass. This seems to mean the ratio is nearly constantirrespective of the mass or size of the clouds.
Since there is a mutual relation between mass, size, and the velocity dispersion to achive a mechanical equlibrium (the Virial relation), there is only one independent correlation in the above twocorrelations (eqs.[1.13] and [1.15]). Although several reasons to explain the correlation are proposed,
we have no consensus yet.

28

CHAPTER 1. INTRODUCTION

Figure 1.28: The left shows the relation between cloud size (holizontal axis) and the three-dimensional
internal velocity (vertical axis). The right shows a similar correlation between mass and the random
velocity.

Chapter 2

Physical Background2.1

Basic Equations of Hydrodynamics

The basic equation of hydrodynamics are (1) the continuity equation of the density [equation (A.11)],∂ρ+ div(ρv) = 0,
∂t

(2.1)

(2) the equation of motion [equation (A.7)]
·

ρ

¸

∂v
+ (v · ∇) v = −∇p + ρg,
∂t

(2.2)

and (3) the equation of energy [equation (A.23)]∂²+ div(² + p)v = ρv · g.
∂t

(2.3)

Occasionally barotropic relation p = P (ρ) substitutes the energy equation (2.3). Especially polytropic relation p = KρΓ is often used on behalf of the energy equation. In the case that the gas is
isothermal or isentropic, the polytropic relations of
p = c2is ρ (isothermal)(2.4)

p = c2s ργ (isentropic)

(2.5)

and
are substitution to equation (2.3). [We can replace equation (2.3) with equations (2.4) and (2.5).]

2.2

The Poisson Equation of the Self-Gravity

In this section, we will show the basic equation describing how the gravity works. First, compare
the gravity and the static electric force. Consider the electric field formed by a point charge Q at a
distance r from the point source as
1 Q
E=
,
(2.6)
4π²0 r2
where ²0 is the electric permittivity of the vacuum. On the other hand, the gravitational accelerationby the point mass of M at the distance r from the point mass is written down asM
g = −G 2 ,
r
29

(2.7)

30

CHAPTER 2. PHYSICAL BACKGROUND

where G = 6.67 × 10−8 kg−1 m3 s−2 is the gravitational constant. Comparing these two, replacing Q
with M and at the same time 1/4π²0 to −G these equations (2.6) and (2.7) are identical with eachother.
The Gauss theorem for electrostatic field as
div E =

ρe
,
²0

(2.8)

and another expression using the electrostatic potential φe as
ρe
,
²0

(2.9)

div g = −4πGρ,(2.10)∇2 φ = 4πGρ,(2.11)∇2 φe = −lead to the equations for the gravity as

andwhere ρe and ρ represent the electric charge density and the mass density. Equation (2.11) is called
the Poisson equation for the gravitational potential and describes how the potential φ is determinedfrom the mass density distribution ρ.
Problem
Consider a spherically symmetric density distribution. Using the Poisson equation, obtain the potential (φ) and the gravitational acceleration (g) for a density distribution shown below.
(ρ= ρ0 for r < R= 0 for r ≥ R

Hint: The Poisson equation (2.11) for the spherically symmetric system is
µ

1 ∂
∂φ
r2
2
r ∂r
∂r

2.3

¶

= 4πGρ.

Free-fall Time

If the pressure force can be neglected in the equation of motion (A.1), the gravitational force remains.
Assuming the spherical symmetry, consider the gravity gr (r) at the position of radial distance fromthe center being equal to r. Using the Gauss’ theorem, gr is related to the mass inside of r, which isexpressed by the equation
Z
Mr =

r

0

ρ4πr2 dr,(2.12)and gr is written asGMr.
(2.13)
r2
This leads to the equation motion for a cold gas under a control of the self-gravity is written
gr (r) = −

d2 r
GMr
=− 2 .
2dt
r(2.14)

Chapter 4

Local Star Formation Process
4.1

Hydrostatic Balance

Consider a hydrostatic balance of isothermal cloud. By the gas density, ρ, the isothermal soundspeed, cis , and the gravitational potential, Φ, the force balance is written as
c2 dρ dΦ
− is
−
= 0,
ρ dr
dr

(4.1)

and the gravity is calculated from a density distribution as
Z r

dΦ
GMr
4πG
−
=− 2 =− 2
dr
r
r

0

ρr2 dr,(4.2)

for a spherical symmetric cloud, where Mr represents the mass contained inside the radius r. The
expression for a cylindrical cloud is
−

Gλr
2πG
dΦ
=−
=−
dr
r
r

Z r
0

ρrdr,

(4.3)

where λr represents the mass per unit length within a cylinder of radius being r.
For the spherical symmetric case, the equation becomes the Lane-Emden equation with the polytropic index of ∞ (see Appendix C.1). This has no analytic solutions. However, the numericalintegration gives us a solution shown in Figure 4.1 (left). Only in a limiting case with the infinite
central density, the solution is expressed as
c2is −2r .
2πG

ρ(r) =(4.4)

Increasing the central density, the solution reaches the above Singular Isothermal Sphere (SIS) solution.
On the other hand, a cylindrical cloud has an analytic solution (Ostriker 1964) asÃρ(r) = ρc

r2
1+
8H 2

!−2

,

(4.5)

where
H 2 = c2is /4πGρc .
69

(4.6)

70

CHAPTER 4. LOCAL STAR FORMATION PROCESS
1

1’fort.4’’1.1.d’’1.d’’0.9.d’

0.1

0.1

0.01

0.01

0.001

0.001

0.00010.0001

1e-05
0.01

0.1

1

10

100

1000

1e-05
0.01

0.1

1

10

100

1000

Figure 4.1: Radial density distribution. A spherical cloud (left) and a cylindrical cloud (right). In theright panel, solutions for polytropic gases with Γ = 1.1 (relatively compact) and Γ = 0.9 (relatively
extended) are plotted as well as the isothermal one.
Far from the cloud symmetric axis, the distribution of equation (4.5) gives
ρ(r) ∝ r−4 ,(4.7)

ρ(r) ∝ r−2(4.8)

while the spherical symmetric cloud has
distribution.
Problem 1
Show that the SIS is a solution of equations (4.1) and (4.2).
Problem 2
Show that the density distribution of [equation (4.5)] is a solution of equations (4.1) and (4.3).

4.1.1

Bonnor-Ebert Mass

In the preceding section [Fig.4.1 (left)], we have seen the radial density distribution of a hydrostatic
configuration of an isothermal gas. Consider a circumstance that such kind of cloud is immersed ina low-density medium with a pressure p0 . To establish a pressure equilibrium, the pressure at thesurface c2is ρ(R) must equal to p0 . This means that the density at the surface is constant ρ(R) = p0 /c2is .
Figure 4.2 (left) shows three models of density distribution, ρc = ρ(r = 0) = 10ρs , 102 ρs , and103 ρs . Comparing these three models, it should be noticed that the cloud size (radius) decreases withincreasing the central density ρc . The mass of the cloud is obtained by integrating the distribution,
which is illustrated against the central-to-surface density√ ratio3 ρc /ρs in Figure 4.2 (right). The y-axis
represents a normalized mass as m = M/[4πρs (cis / 4πGρs ) ]. The maximum value of m = 4.026
means
c2is
Mmax ' 1.14
.
(4.9)
1/2
G3/2 p0

4.1. HYDROSTATIC BALANCE

71

This is the maximum mass which is supported against the self-gravity by the thermal pressure with anisothermal sound speed of cis , when the cloud is immersed in the pressure p0 . This is called BonnorEbert mass [Bonnor (1956), Ebert (1955)]. It is to be noticed that the critical state Mcl = Mmax isachieved when the density contrast is rather low ρc ' 16ρs ≡ ρcr .
Another important result from Figure 4.2 (right) is the stability of an isothermal cloud. Even fora cloud with Mcl < Mmax , any clouds on the part of ∂Mcl /∂ρc < 0 are unstable, whose clouds are
distributed on the branch with ρc > ρcr . This is understood as follows: For a hydrostatic cloud themass should be expressed with the external pressure and the central density [Fig.4.2 (right)] as
Mcl = Mcl (pext , ρc ).

(4.10)In this case, a relation between the partial derivatives such as
µ

∂Mcl
∂pext

¶

µ

∂pext
∂ρc

·
ρc

¶

µ

·
Mcl

∂ρc
∂Mcl

¶

= −1,

(4.11)

pext

is satisfied, unless each term is equal to zero. Figure 4.1 (left) shows that the cloud mass Mcl is adecreasing function of the external pressure pext = ρs c2is , if the central density is fixed. Since this
means
µ
¶
∂Mcl
< 0,
(4.12)∂pext ρcequation (4.11) gives us

µ

∂pext
∂ρc

¶

µ

·
Mcl

∂ρc
∂Mcl

¶

> 0.

(4.13)

pext

For a cloud with ρc < ρcr = 16pext /c2is the mass is an increasing function of the central density as
µ

∂Mcl
∂ρc

¶

> 0.

(4.14)> 0,(4.15)

pext

Thus, equation (4.13) leads to the relation
µ

∂ρc
∂pext

¶
Mcl

for ρ < ρcr . This means that gas cloud contracts (the central density and pressure increase), whenthe external pressure increases. This is an ordinary reaction of a stable gas.
2
On the other hand, the cloud on the part of (∂Mcl /∂ρc )pext < 0 (for 16<∼ρc /(pext /cis )<∼2000)
behaves
¶
µ
∂ρc
< 0,
(4.16)∂pext Mcland this represents that an extra external pressure decreases the central density and thus the pressure.
Pressure encourages expansion. This reaction is unstable.

4.1.2

Equilibria of Cylindrical Cloud

In Figure 4.1 (right) we plotted the structure for a polytropic cloud. Inner structure is not dependentof Γ, it is clear the slope of the outer envelope is dependent on Γ.

72

CHAPTER 4. LOCAL STAR FORMATION PROCESS
1000

5’F1000.d’’F100.d’’F10.d’’rhoc-M.d’

4

100

M

rho

3

2
10

1

1
0.001

0
0.01

0.1

1

10

1

x

10

100

1000

10000

rho

Figure 4.2: (Left) radial density distribution. Each√curve has different ρc . The x-axis denotes theradial distance normalized by a scale-length as c2is / 4πGpext . It is shown that the radius increases
with decreasing ρc in this range (ρc = (10, 102 , 103 ) × pext /c2is ). (Right) The relation between mass
(Mcl ) and the central density (ρc ) is plotted, under the condition of constant external pressure. The
x-axis represents the central √density normalized by ρs ≡ pext /c2is . The y-axis represents the cloudmass normalized by 4πρs (cis / 4πGρs )3 .
1. In the case of the spherically symmetric, consider a polytrope (p ∝ ρΓ ) with Γ < 6/5 (at least
the envelope of Γ = 6/5 cloud extends to ∞.), in which gas extends to ∞. if ρ ∝ r−p , themass inside of r is proportional to Mr ∝ r3−p . Thus, the gravity per unit volume at a radiusr, GM ρ/r2 , is proportional to GM ρ/r2 ∝ r1−2p . On the other hand the pressure force is|(∂p/∂r)| = (∂p/∂ρ)|(∂ρ/∂r)| ∝ (r−p )Γ−1 r−p−1 ∝ r−pΓ−1 . These two powers become the same,
only if p = 2/(2 − Γ).
2. In the case of cylindrical cloud with Γ ≤ 1, the mass per unit length λ ∝ r2−p . The gravity atr, Gλρ/r ∝ r1−2p . Note that the power is the same as the spherical case. Since the power ofthe pressure force should be the same as the spherical case, the resultant p should be the samep = 2/(2 − Γ).
The case of Γ = 0.9, an envelope extending to a large radius indicates the power-law distribution
much shallower than that of the isothermal Γ = 1 one.

4.2

Virial Analysis

Hydrodynamic equation of motion using the Lagrangean time derivative [eq.(A.3)] is
µ

dv
ρ
dt

¶

= −∇p − ρ∇Φ.

(4.17)

Multiplying the position vector r and integrate over a volume of a cloud, we obtain the Virial relationas
1 d2 I= 2(T − T0 ) + W,
(4.18)
2 dt2
where
Z
I = ρr2 dV,(4.19)

4.2. VIRIAL ANALYSIS

73

is an inertia of the cloud,

Z µ

T =

¶

3
1
3
pth + ρv 2 dV = p̄Vcl ,2
2
2(4.20)

is a term corresponding to the thermal pressure plus turbulent pressure,
Z

T0 =

3
Pth r · ndS = P0 Vcl
2S(4.21)

comes from a surface pressure, andZW =−ρr · ∇ΦdV = −3 GM 25 R(4.22)

is a gravitational energy. To derive the last expression in each equation, we have assumed the cloudis spherical and uniform. Here we use a standard notation as the radius R, the volume Vcl = 4πR3 /3,
the average pressure P̄ , and the mass M .
To obtain a condition of mechanical equilibrium, we assume d2 I/dt2 = 0. Equation (4.18) becomes
4π p̄R3 − 4πp0 R3 −3 GM 2= 0.
5 R(4.23)

Assuming the gas is isothermal p = c2is ρ, the average pressure is written asp̄ = c2is ρ̄ = c2is

3M
.
4πR3

(4.24)

Using equation (4.24) to eliminate p̄ from equation (4.18), the external pressure is related to the massand the radius as
3c2 M
3GM 2
p0 = is 3 −.
(4.25)
4πR
20πR4
Keeping M constant and increasing R from zero, p0 increases first, but it takes a maximum, p0,max =3.15c8is /(G3 M 2 ), and finally declines. This indicates that the surface pressure must be smaller thanp0,max for the cloud to be in an equilibrium. In other words, keeping p0 and changing R, it is shownthat M has a maximum value to have a solution. The maximum mass is equal toMmax = 1.77

c4is
1/2

G3/2 p0.

(4.26)

The cloud massive than Mmax cannot be supported against the self-gravity. This corresponds to the
Bonnor-Ebert mass [eq.(4.9)], although the numerical factors are slightly different.

4.2.1

Magnatohydrostatic Clouds

Consider here the effect of the magnetic field. In the magnetized medium, the Lorentz force
F =

1
1
1
(∇ × B)×B = − ∇B 2 +(B · ∇)B
4π
8π
4π(4.27)

works in the ionized medium. The first term of equation (4.27), which is called the magnetic
pressure, has an effect to support the cloud against the self-gravity.

74

CHAPTER 4. LOCAL STAR FORMATION PROCESSThe virial analysis is also applicable to the magnetohydrostatic clouds. The terms related to themagnetic fields are
Z

M

Z

Z

B2
B2
dV + (r · B)B · ndS −r · ndS
8π
S
S! 8π
Ã
Z
B 2 − B02
1
Φ2B
Φ2
'
dV ' 2
− B ,
8π
6π
R
R0
=

(4.28)

where ΦB represents a magnetic flux and it is assumed to be conserved if we change the radius, R,
that is ΦB = πB0 R02 = πBR2 . Equation (4.23) becomes
4π p̄R3 − 4πp0 R3 −3 GM 21 Φ2+ 2 B = 0,
5 R
6π R

(4.29)

Φ2where we ignored the term 6π1 2 RB0 . The last two terms are rewritten as´
3G³ 2
M − MΦ2 ,5R(4.30)

where MΦ is defined as 3GMΦ2 /5 = Φ2B /6π 2 .
This shows the effects of the magnetic fields:
1. B-fields effectively reduce the gravitational mass as M 2 −MΦ2 = M 2 −5Φ2B /(18π 2 G). This plays
a part to support a cloud.
2. However, even a cloud contracts (decreasing its radius from R0 to R), the ratio of the gravitational to the magnetic terms keeps constant since these two terms are proportional to ∝ R−1 .
Thus, if the magnetic term does not work initially, the gravitational term continues to predominate over the magnetic term.
If M < MΦ , a sum of last two terms in equation (4.29) is positive. Since the second term of rhs ofequation (4.25) is positive, there is one R which satisfies equation (4.29) irrespective of the externalpressure p0 . While, if M > MΦ , there is a maximum allowable external pressure p0 . Therefore,
M = MΦ gives a criterion whether the magnetic fields work to support the cloud or not. Morerealistic calculation [Mouschovias (1976a,1976b), Tomisaka et al (1988)] gives us a criterion
G1/2

G1/2 σ1dm== 0.17 '
,
dΦB
B
2π

(4.31)

where σ and B means the column density and the magnetic flux density. A cloud with a massM>ΦB
2πG1/2(4.32)

is sometimes called magnetically supercritical, while that withM<

ΦB
2πG1/2(4.33)

is magnetically subcritical.
More precisely speaking, the criterion showed in equations (4.32) and (4.33) should be applied fora cloud which has a much larger mass than the Bonnor-Ebert mass. That is, even without magnetic
fields, the cloud less-massive than the Bonnor-Ebert mass has a hydrostatic configuration shown in4.3. SUBCRITICAL CLOUD VS SUPERCRITICAL CLOUD

75

Figure 4.2 (left). The cloud with central density of ρc = 10 has a stable density distribution. To fit
the numerical results, Tomisaka et al (1988) obtained an expression for the critical mass when thecloud has a mass-to-flux ratio dm/dΦB , the isothermal sound speed cis , and the external pressure p0
as(·¸2 )−3/2
1/2πc4s
Mcr = 1.3 1 −.
(4.34)
1/2
G1/2 dm/dΦB |r=0
p G3/20This shows that the critical mass is a decreasing function of the mass-to-flux ratio or increasing
function of the magnetic flux. And the critical mass becomes much larger than the Bonnor-Ebert
1/2
mass ' c4s /(p0 G3/2 ) only when the mass-to-flux ratio at the center of the cloud is reaching 1/2π
at which the term in the curry bracket goes to zero. Hereafter, we call here the cloud/cloud corewith mass larger than the critical mass Mcr a supercritical cloud/cloud core. The cloud/cloud coreless-massive than the critical mass is subcritical.

4.3

Subcritical Cloud vs Supercritical Cloud

We have seen there is a critical mass above which the cloud has no (magneto)hydrostatic configuration
but below which the cloud has at least an equilibrium state. This gives us an idea that there are
two kind of clouds/cloud cores: that with a mass larger than the critical mass which has to collapse
dynamically and that with a mass smaller than the critical mass which is in an equilibrium state.
10−3
In the density range of 104 cm−3 <∼10 cm , the interstellar gas is essentially isothermal. In this∼n<region a major cooling agent is dusts; that is, the dust is heated by the collision of molecules. The
excess energy liberated at the collision increases the dust temperature. Finally the thermal emission
from the dust cools down the dust again. By this process, the thermal energy of the gas is reduced.
Therefore, we consider the cloud/cloud core is isothermal and study the collapse of the isothermal
cloud.
Since the supercritical cloud has no hydrostatic configuration, it must evolve in a dynamical way.
On the other hand, since the subcritical cloud is in a static state, it evolves in much longer time-scaleof the free-fall time. Such cloud evolves by the effect called the ambipolar diffusion.

4.4

Ambipolar Diffusion

4.4.1

Ionization Rate

In the dense clouds, the ionization fraction is low. Since the uv/optical radiations from stars can notreach the cloud center, potential ionization comes from the cosmic ray particles. In this case the rateof ionization per volume is given as ζnn , where ζ ∼ 10−17 s−1 . In Figure 4.3, the ionization fraction forvarious density is shown (Nakano, Nishi, & Umebayashi 2002). This clearly shows that the fraction of
−1/2
8
−3
ions decrease approximately in proportion to nH
for nH <∼10 cm . This is understood as follows:+
Equilibrium balance between one kind of ion m and its neutral m0
m 0 ↔ m + + e− ,(4.35)

is considered. The recombination (reaction from the right to left) rate per unit volume is expressed
αne nm+ , while the ionization rate (left to right) per unit volume is ζnm0 , where α means the recom8−3
bination rate coefficient. Nakano (1984) obtained in the range of 102 cm−3 <∼n<∼10 cmρi = Cρ1/2
n ,where the numerical factor C = 4.46 × 10−16 g1/2 cm−3/2 .

(4.36)

76

CHAPTER 4. LOCAL STAR FORMATION PROCESS

Figure 4.3: Ionization fractions X/nH are plotted against the number density nH . The ionization rateof an H2 molecule by cosmic rays outside the cloud is taken ζ0 = 1 × 10−17 s−1 . 20% of C and O and
2% of metallic elements are assumed to remain in the gas phase and the rest in grains.

4.4.2

Ambipolar Diffusion

The ionization fraction of the gas in the cloud is low. Neutral gas and the ions are coupled via
mutual collisions and ions are coupled with the magnetic field. Thus, the neutral molecules, a major
component of the gas, are coupled with the magnetic field indirectly via ionized ions. Equation of
motion for the molecular gas isρdun= −c2s ∇ρ − ρ∇ψ + ρrω 2 er + F ,dt(4.37)

where the forces appeared in the rhs represent, respectively, the pressure force, the self-gravity, thecentrifugal force, and the force exerted on the neutral component per unit volume through two-bodycollision with ions. The friction force has a form
F = ρi ρn γ(v i − v n ).

(4.38)We have assumed ρ ' ρn since the mass density of the charged component is much smaller than thatof the neutral one. On the other hand, that for ions is
ρi

du
1
=
(∇ × B)×B − F ,dt
4π(4.39)

where we ignored the pressure and self-gravity forces compared with the Lorentz force. If the inertiaof the ions are ignored [lhs of eq.(4.39)=0], equations (4.37) and (4.39) give the equation of motion
similar to that of the one-fluid as follows:ρ

du1= −c2s ∇ρ − ρ∇ψ + ρrω 2 er +(∇ × B)×B.
dt
4π(4.40)

4.4. AMBIPOLAR DIFFUSION

77

Figure 4.4: Relation central density ρc against mass-to-flux ratio M/Φ. Plane-parallel and axisymmetric models are illustrated in the left and right panels, respectively.
Ignoring the inertia of the ions, equation (4.39) indicates that the friction force should be balancedwith the Lorentz force as
1(∇ × B)×B.
(4.41)F =4π
Using equations (4.36) and (4.38), this gives the drift velocity of ions against the neutrals as
vD ≡ vi − vn,F,=ρi ρn γ1=(∇ × B)×B,
4πρi ρn γ
1=(∇ × B)×B.
3/2
4πCγρn(4.42)

Since the Lorentz force works outwardly, the drift velocity v D also points outwardly. Viewing fromthe ions and magnetic field lines, v n − v i points inwardly. Thus, the mass in a flux tube varies withtime; that of the cloud center increases while that of the cloud surface decreases. The ambipolar
diffusion is a process to change the mass distribution against the magnetic flux tube.
How is the structure affected by the increase of mass/flux ratio at the cloud center? Here, we will
derive the relation between the mass-to-flux ratio and central density and answer the question. Numerical calculation is necessary for the axisymmetric (realistic) cloud. However, this can be obtained
for the case of a plane-parallel disk in which the magnetic field is running parallel to the disk, For
the plane-parallel disk the hydrostatic balance is achieved as∂∂zwhere σ =

Rz

Ã

B2
p+
8π

!

= −ρ

∂ψ
= −2πGρσ,∂z(4.43)−z ρdz. Equation (4.43) is integrated to obtain the relation between the variables at the

78

CHAPTER 4. LOCAL STAR FORMATION PROCESS

center z = 0 as

Ã

B2
B2
pc + c − ps + s
8π
8π

!

=

πG 2
σ ,
2 s

(4.44)

where we used dσ/dz = 2ρ. This is rewritten as
ρc
1
+
ρs 8π

µ

Bc
ρc

¶2

ρs
c2s

µ

ρc
ρs

¶2

1
−
8π

µ

Bs
ρs

¶2

ρs
πGσs2
=
+ 1.
c2s
2ρs c2s(4.45)

1/2

In the second term of lhs, yc ≡ (ρc /Bc )(cs /ρs ) is a nondimensional mass flux ratio at the center andthis is related to the center-to-surface density ratio ρ0c ≡ ρc /ρs . Using these nondimensional variables,
equation (4.45) is expressed as(ρ0 )21
ρ0c + c 2 −= ρ̄0c ,(4.46)
8πyc
8πys2
where ρ̄0c ≡ πGσs2 /(2ρs c2s ) + 1 represents the central density necessary for the disk to be supported
without a magnetic field (y → ∞). Figure 4.4 (left) illustrates this relation. This shows that if thecentral mass-to-flux ratio increases (moving upward), the central density increases. In low density,
ρ0c ¿ ρ̄0c , equation (4.46) indicates that ρ0c ∝ yc ∝ ρc /Bc . This means that on the low density regime,
the density increases but the strength of magnetic field does not increase. Equation (4.46) also
indicates that yc increase much as yc → ρ̄0c ys when ρ0c → ρ̄0c . This relation between mass-to-flux ratio
and the central density is well fitted to the numerical calculation of plane-parallel disk cloud driven
by the ambipolar diffusion by Mouschovias, Paleologou, & Fiedler (1985). Does this indicate that the
mass-to-flux ratio increases much when the central density increases, for example, from 104 cm−3 to106 cm−3 ?
From this point, quasistatic evolution of the magnetized cloud was first studied by Nakano (1979)
adopting a method seeking for magnetohydrostatic configuration (Mouschovias 1976a, 1976b). He
obtained that the magnetic flux to mass ratio near the center of the cloud does not decrease much
below a critical density above which magnetohydrostatic configurations no more exist. Paleologou &Mouschovias (1983), Shu (1983), and Mouschovias, Paleologou, & Fiedler (1985) gave a completely
different result. That is, the magnetic flux to mass ratio near the center decreases much when the
6−31density becomes larger than n>∼10 cm . Critical discussion has been done between them . However,
we have now realized that the difference comes from the geometry assumed.
Relation of the mass-to-flux ratio to the central density for the axisymmetric realistic (not planeparallel disk) cloud is plotted in Figure 4.4 (right) (Tomisaka, Ikeuchi, & Nakamura 1990). The central
mass-to-flux ratio increases by the ambipolar diffusion. This increase the central density. Lines witharrows indicate the evolutionary paths driven by the ambipolar diffusion. It is to be noticed thatthe relation is completely different from that of the plane-parallel disk. Increase of the mass-flux
ratio is small although the central density increases much. This indicates that in the quasistatic
evolution driven by the ambipolar diffusion the mass-flux ratio does not increase much in contrast to
the plane-parallel case.
In fully ionized plasma, the magnetic field is coupled to the matter. However, in the low ionizedgas the magnetic field is coupled only to the ions. Thus the induction equation of magnatic fieldshould be written∂B+ ∇ × (vi × B) = 0.
(4.47)∂t1

One of the reasons why large increase in the mass/flux ratio is favored is understood as follows: There havebeen a long-standing “magnetic flux problem of stars” in which the magnetic flux of, say, 1M¯ main-sequence starΦ ∼ XXG cm2 is much smaller than that of the parent cloud Φ ∼ XXG cm2 . The magnetic flux must be reduced in
the star formation process. If the mass-flux ratio increases much in these density range, this might resolve the magnetic
flux problem of stars.

4.5. DYNAMICAL COLLAPSE

79This leads to the final expression as∂B+ ∇ × (vn × B) + ∇ ×
∂t

½

¾
4C

1

3n=2

[(∇ × B) ×B] × B= 0,(4.48)

where we used the expression of the drift velocity [equation (4.42)]. In the isothermal cloud, there
is a maximum mass supported by the magnetic field. Since in the subcritical cloud M < Mcr themagnetic flux escapes from the center, the critical mass decreases with time. After the critical mass
becomes smaller than the actual mass (Mcl > Mcr ), there is no hydrostatic configuration. The cloud
evolves into the supercritical cloud region and experiences dynamical collapse.

4.5

Dynamical Collapse

Figure 4.5: Evolution of isothermal clouds massive than the Bonnor-Ebert mass. Density (left) andvelocity (right) distributions are illustrated. Solid lines show the cores of the preprotostellar phase(prestellar core) and dashed lines show those of protostellar phase (protostellar core). The evolutionof the protostellar phase is studied by the sink-cell method, where we assume the gas that entered in
the sink-cells is removed from finite-difference grids and add to the point mass sitting at the centerof the sink-cells which corresponds to a protostar.
In 1969, Larson (1969) and Penston (1969) found a self-similar solution which is suited for the
dynamical contraction. Figure 4.5 is a radial density distribution for a spherical collapse of an
isothermal cloud, where the cloud has a four-times larger mass than that of the Bonnor-Ebert mass.
Although the figure is taken from a recent numerical study by Ogino et al (1999), a similar solution
was obtained in Larson (1969). We can see that the solution has several characteristic points asfollows:
1. The cloud evolves in a self-similar way. That is, the spatial distribution of the density (left) at
t = 0.326 is well fitted by that at t = 0.293 after shifting in the −x and the +y directions. As
for the infall velocity spatial distribution, only a shift in the −x direction is needed.

80

CHAPTER 4. LOCAL STAR FORMATION PROCESS
2. The density distribution in the envelope, which is fitted by ∝ r−2 , is almost unchanged. Only
the central part of the cloud (high-density region) contracts.
3. The time before the core formation epoch (the core formation time t0 is defined as the time atwhich the central density increases greatly) is a good indicator to know how high the central
density is. That is, reading from the figure, at t − t0 = 0.003 (t = 0.326) the central density
reaches ρc ' 104 and at t − t0 = 0.036 (t = 0.293) the density is equal to ρc ' 102 . This shows
the maximum (central) density is approximately proportional to (t − t0 )2 , which is reasonablefrom the description of the free-fall time ∝ ρ−1/2 .

The basic equations of spherical symmetric isothermal flow are∂ρ
1 ∂+ 2 (r2 ρv) = 0,
∂t
r ∂r
∂v
∂v c2is ∂ρ GM
+v
+
+ 2
= 0,
∂t
∂r
ρ ∂r
r
M (r, t) = M (0, t) +(4.49)(4.50)
Z r
0

4πr02 ρ(r0 , t)dr0 ,(4.51)

where M (r, t) represents the mass included in the radius r and M (0, t) denotes the mass of the
protostar. A self-similar solution which has a formΩ(ξ),
4πG(t − t0 )2v(r, t) = cis V (ξ),
c3 |t − t0 |M (r, t) = is
m(ξ)
G
r
ξ =
,
cis |t − t0 |ρ(r, t) =(4.52)(4.53)(4.54)(4.55)

should be found, where Ω and V are functions only on ξ. For example, equation (4.52) asks theshape of the density distribution is the same after resizing of equation (4.55) r → r/cs |t − t0 | andre-normalizing ρ → ρ · 4πG(t − t0 )2 as Ω(ξ). Since
∂
d
=
∂t
dξ
and

µ

∂
d
=
∂r
dξ

∂ξ
∂t

µ

¶

∂ξ
∂r

=

r
d
,
cis |t − t0 |2 dξ

(4.56)=

1d,cis |t − t0 | dξ(4.57)

r¶t

the basic equations for the spherical symmetric model yieldm = (ξ − V )ξ 2 Ω,h(ξ − V )2 − 1

i dV

dξ

·

(4.58)¸= Ω(ξ − V ) −2(ξ − V ),
ξ

·

¸

2
(ξ − V )2 − 1 dΩ= Ω − (ξ − V ) (ξ − V ),Ωdξ
ξ

(4.59)(4.60)

Equations (4.59) and (4.60) have a singular point at which (ξ−V )2 −1 = 0 or V = ξ±1. Since the pointof ξ =const moves with cis , the flow velocity relative to this ξ =const is equal to v − cis ξ = (V − ξ)cis .

4.5. DYNAMICAL COLLAPSE

81

Thus the singular point ξ∗ at which V = ξ ± 1 corresponds to a sonic point. Therefore, since the flow
has to pass the sonic point smoothly, the rhs of equations (4.59) and (4.60) have to be equal to zeroat the singular point ξ∗ . This gives at the sonic pont ±ξ∗ > 0,Ω−2(ξ − V )2
= Ω(ξ − V ) − = 0,ξξ(4.61)

which leads toV∗ = ξ∗ ∓ 1,

(4.62)

2.
ξ∗(4.63)Ω∗ = ±

These equations (4.58), (4.59) and (4.60) have an analytic solution
V (ξ) = 0, Ω =2, m = 2ξ, −∞ < ξ < ∞ξ2(4.64)

This is a solution which agrees with the Chandrasekhar’s SIS. Generally, solutions are obtained only
by numerical integration. |ξ → ∞| the solution have to converge to an asymptotic form of
A − 2 V∞ 4V∞ + (A − 2)(A − 6)+ 2 ++ O(ξ −4 ),
ξ
ξ
6ξ 3
A
Ω∞ (A − 2)−
+ O(ξ −6 ),ξ22ξ 4

V (ξ) = V∞ −(4.65)Ω(ξ) =(4.66)

This shows that for sufficiently large radius the gas flows with a constant inflow velocity V∞ cis .
This has a solution in which the density and the infall velocity should be regular with reaching
the center (ξ ¿ 1). Such kind of solution is plotted in Figure 4.6 (left) and Figure 4.6 (right) witht < 0. This time evolution is expected from the self-similar solution. This shows that
(

ρ
(

v

' ρc
(in the central region)∝ r−2 (in the outer envelope)(4.67)∝r(in the central region)' 3.28cis (in the outer envelope)(4.68)

Reaching the outer boundary the numerical solution (Fig.4.5) differs from the self-similar solution(Fig.4.6). For example, v is reduced to zero in the numerical simulations, while it reaches a finite value
3.28 in the self-similar solution. And as for the density distribution, ρ drops near the outer boundary
in the numerical simulations while it decreases proportional to ∝ r−2 . However, in the region except
for the vicinity of the outer boundary the self-similar solution expresses well the dynamical collapseof the spherical isothermal cloud. This solution gives the evolution of a pre-protostellar core formed
in a supercritical cloud/cloud core.

4.5.1

Inside-out Collapse Solution

In 1977 Shu found another self-similar solution which is realized after a central protostar with infinitesimal mass is formed in the singular isothermal sphere solution. The gas begins to accrete to theprotostar. Outside the region where the accretion occurs, the initial SIS is kept unchanged, since the
SIS is a hydrostatic solution. And the front of accretion expands radially outward in time. Since theinflow region expand outwardly, he called it the inside-out collapse solution. In Figure 4.7, the82

CHAPTER 4. LOCAL STAR FORMATION PROCESS

16

3000
t=0y

14

+104
12 -10vr [m s-1]

log n [m-3]

2000
4

+105
10

-105
+106

8

+104 +105 +106
1000
t=0y

-106

-104
6
13

14

15
16
log r [m]

17

0
13

14

-105

-106

15
16
log r [m]17

Figure 4.6: A self-similar solution indicating a dynamical collapse of isothermal spherical cloud(Larson-Penston solution). Spatial distribution of the density and inflow velocity which are expectedfrom the self-similar solution are plotted. Dashed lines show the evolution prestellar core and solidlines show that of protostellar core. Taken from Hanawa (1999).

4.5. DYNAMICAL COLLAPSE

83

evolution is shown for density and inflow velocity. This solution gives the evolution of a protostellar
core formed in a subcritical cloud/cloud core.
Using equations (4.65) and (4.66) and assuming the inflow velocity should reduce at large radius,
we obtain
A − 2 (A − 2)(A − 6)−
+ ···,(4.69)V (ξ) = −ξ6ξ 3Ω(ξ) =

A(A − 2)
A−
+ ···.
ξ22ξ 4(4.70)

Since ξ → ∞ means t → t0 (if r is finite), Ω(ξ) → A/ξ 2 means thatρ(r, t0 ) =

Ac2is
.
4πGr2

(4.71)

Comparing with the SIS, when A = 2 this gives the SIS and when A > 2 this gives a density
distribution in which the pressure is inefficient and the cloud is contracting. The solution with A > 2
is obtained by a procedure as (1) at a sufficiently large radius ξ1 , calculate V (ξ1 ) and Ω(ξ1 ). (2) fromthese values, integrate equations (4.59) and (4.60) inwardly. Figure 4.7 show the solution of this type.
The solution with A > 2 inflow speed is accelerated towards the center. Decreasing A (A → 2), it+
is shown that an outer part ξ >∼1 reaches V → 0. For A = 2 , the solution reaches the singular line2
V = −ξ + 1 at ξ = 1 (V = 0) . Since V = 0 and Ω = 2 at ξ = 1, this solution with A = 2+ converges
to the SIS at ξ = 1. This means that if there is an infinitesimally small amount of excess mass at thecenter of SIS, the accretion begins from the center while outside a radius the cloud is left static. The
1/2 and Ω ∝ ξ 3/2 .
inner part of the solution ξ <∼1, V and Ω are well expressed as V ∝ ξ

Figure 4.7: Self-similar solution which shows the inside-out collapse. (Left:) Infall velocity (−V =−v/cis ) is plotted against the similarity variable ξ ≡ r/cis (t − t0 ). Three curves correspond to models
A = 2.4, 2.1, and 2.01. (Right:) Density (Ω = 4πGρ(t − t0 )2 ) is plotted against the similarity variableξ ≡ r/cis (t − t0 ).
2This means A = 2 + ² and ² > 0 and → 0

84

CHAPTER 4. LOCAL STAR FORMATION PROCESSThe power-law distributions of V ∝ ξ 1/2 and Ω ∝ ξ 3/2 are explained as follows: Conservation of
the total energy of an inflowing gas shell is expressed as
v 2 GMr
GMr0
−
=−
,
2
r
r0

(4.72)

where r0 denotes the initial radius of a gas shell and Mr0 = Mr represents the mass inside the gas
shell. Neglecting GMr0 /r0 compared with the term GMr /r, this gives
µ

v'

2GMr
r

¶1/2

∝ r−1/2 ,(4.73)

whereR we assumed a major part of Mr comes from the mass of a protostarM∗ , thatis, Mr =R
R
M∗ + 0r ρ4πr2 dr ' M∗ Since the average density inside the radius r ρ̄(< r) = 0r ρ4πr2 dr/ 0r 4πr2 dr =3ρ(r) for SIS distribution, the time necessary for a gas shell to reach the center is proportional to∝ ρ̄(< r)−1/2 ∝ ρ(r)−1/2 ∝ r1 ∝ (t − t0 )1 , where we used a fact that the front of accretion expands
with a constant speed cs is. This means that the time necessary for the gas shell to travel from theradius of the accretion wave front to the center is proportional to t − t0 . Since the mass of the shell
which begins accretion in a unit time is equal to ρ(r)4πr2 cis and is constant irrespective of t − t0 .
These two facts indicate that the mass accretion rate is constant in time. That is,Ṁ = 4πr2 ρv = const.

(4.74)

Since the mass swept by the sound wave per unit time is equal to Ṁsw = 4πr2 ρSIS cs = 2c3s /G, theaccretion rate is proportional to c3s /G. Shu (1977) obtained the accretion ratec3
Ṁ = 0.975 s ,G(4.75)

for SIS A = 2+ . Using this equation, equations (4.73) and (4.74) indicate that the spatial density
distribution is expressed by a power-law as
ρ(r) ∝ r−3/2 ,(4.76)

which is valid for the region except for the vicinity of the front of the accretion.

4.5.2

Protostellar Evolution of Supercritical Clouds

What is a protostellar core formed in a supercritical cloud/cloud core? Is this different from theinside-out solution of Shu (1977)? A solution which corresponds to the protostellar core is obtainedby Hunter (1977) and Whitworth and Summers (1985). This is a solution with t > t0 in equation(4.55). The asymptotic behaviors of the density and infall velocity reaching the center are different
from that of the Larson-Penston self-similar solution for a prestellar collapse. That is,
(

Ω
(

V

→ finite(LP)→ infinite (Inside − out)

(4.77)→ finite(LP)→ infinite (Inside − out)

(4.78)

Using the boundary conditions suitable for the inside-out type solution, another self-similar solution
is obtained. In Figure 4.6, such kind of solution is also plotted for t > 0.

4.6. ACCRETION RATE

85

Take notice that the solutions of t < t0 (prestellar) and t > t0 (protostellar) agree with eachother at t = t0 . Even if the boundary conditions at the center for the similarity variables, V andΩ, are completely different, the difference between the two is small in the physical variable v andρ. Therefore, the evolution of a supercritical core is thought to be expressed by the Larson-Penston
self-similar solution extended to the protostellar core phase by Hunter (1977) and Whitworth andSummers (1985).
Assume that we observe a protostellar core and obtain their density and infall velocity spatial
distributions. Can we distinguish which solution is appropriate for the Shu’s inside-out solution or theextended Larson-Penston solution? This seems hard, because the structure of density and velocity
distributions are similar after the protostar is formed: the density and velocity show almost similar
power-law as ρ ∝ r−3/2 and v ∝ r−1/2 irrespective of the inside-out solution or the extended LarsonPenston solution. The region where the infall velocity is accelerated toward the center (accretiondominated region) is expanding after the protostar is formed. Therefore, to distinguish between thetwo solutions becomes harder and harder after the protostar is formed. The difference would be large
and we would have a definite answer which solution is appropriate to describe the cloud collapse, if we
can observe a very young protostellar core or a preprotostellar core which shows dynamical collapse.
However, since the time-scale of such a phase is much shorter than the evolved protostellar phaseor a younger preprotostellar core, the number of such kind of objects would be small (τff ∝ ρ−1/2 ).
Therefore, we are looking for such objects just before or after the protostar formation.

4.6

Accretion Rate

Using equation (2.26), the necessary time for a mass-shell at R to reach the center (free-fall time) isexpressed asÃ

T (R) ≡

R3
2GM (R)!1/2

π2(4.79)(for detail of this section see Ogino, Tomisaka, & Nakamura 1999).
Consider two shells whose initial radii are R and R + ∆R. The time difference for these two shells
to reach the center ∆T (R) can be written down using equation (4.79) as
·

¸

πR1/2
3
R dM (R)∆T (R) = 3/2
−
∆R.
1/2
2 2M (R) dR
2 (GM (R))

(4.80)

Mass in the shell between R and R + ∆R, ∆M ≡ M (R + ∆R) − M (R) = (dM/dR)∆R, accretes
on the central object in ∆T (R). Thus, mass accretion rate for a pressure-free cloud is expressed as∆M/∆T . This leads to the expression asdM (R)dR.
dM (R)
2M (R) dR

R

dM
23/2 G1/2 M (R)3/2 M (R)
(R) =
3
dT
π
R3/2
− R
2

(4.81)

This gives time variation of the accretion rate. Consider two clouds with the same density distribution ∂ log ρ/∂r but different absolute value. Since these two clouds have the same ∂ log M (R)/∂ log R,
the mass accretion rate depends only on M (R)/R, and is expressed asdM(R) ∝ M (R)3/2 .
dT(4.82)

86

CHAPTER 4. LOCAL STAR FORMATION PROCESS

Figure 4.8: Mass accretion rate against the typical density of the cloud.
This indicates that the accretion rate is proportional to ρ3/2 , while the time scale is to ρ−1/2 . This isconfirmed by hydrodynamical simulations of spherical symmetric isothermal clouds (Ogino et al.1999).
When the initial density distribution is the SIS as ρ ∝ r−2 , the mass included inside R is proportional
to radius M (R) ∝ R. In this case, equation (4.81) gives a constant accretion rate in time. In Figure
4.8 we plot the mass accretion rate against the cloud density. α represents the cloud density relative
to that of a hydrostatic Bonnor-Ebert sphere. This shows clearly that the mass accretion rate is
proportional to α3/2 for massive clouds α > 4. This is natural since the assumption of pressure-lessis valid only for a massive cloud in which the gravity force is predominant against the pressure force.
Similar discussion has been done by Henriksen, André, & Bontemps (1997) to explain a decline
in the accretion rate from Class 0 to Class I IR objects. They assumed initial density distribution ofρ = ρ0 = ρ0³(r ≤ rN ),r
rN´−2/D1(r > rN ),(4.83)as shown in Figure 4.9. Since the free-fall-time of the gas contained in the inner core r ≤ rN is thesame, such gas reaches the center once. It makes a very large accretion rate at t = (3π/32Gρ0 )1/2
1
as Ṁ = (4π/3)ρ0 δ(t − (3π/32Gρ0 )1/2 ). If D1 = 1, ρ ∝ r−2 for r>∼rN . Since M ∝ R and tff ∝ R,
equation(4.81) predicts Ṁ ∝ R0 ∝ t0 . A constant accretion rate is expected for this power-law andthe accretion rate is converged to a constant value after M (rN ) ¿ M (r). If D1 = 2, ρ ∝ r−3 for2
1/2
3/2 ∝ t3 . Theyr>∼rN . Since for this power M ∝ R and tff ∝ R , equation(4.81) predicts Ṁ ∝ R
gave ṁ ∝ t3D1 −3 for D1 > 2/3.
ProblemShow that Ṁ ∝ t3D1 −4/3D1 if ρ ∝ r−2/D1 and D1 < 2/3 [Henriksen, André, & Bontemps (1997)].

4.7. OUTFLOW

87

Figure 4.9: A model proposed to explain time variation in accretion rate by Henriksen, André, &Bontemps (1997). The density distribution ρ(r) at t = 0 (left) and expected accretion rate (right).

4.7

Outflow

In chapter 1, we mentioned ‘outflows’ ejected from the protostars and pre-main-sequence stars. L1551
IRS5 is a typical example which shows a number of outflows are ejected in the course of star formation.
One is molecular outflow, whose lobes extends L ∼ 0.5(D/160pc)pc in two opposite directions fromthe IRS-5 (Snell et al 1980). This is traced by CO J = 1 → 0 emission line. They estimated an
expansion velocity of the CO gas of VCO ∼ 15km s−1 , a mass of 0.3 M¯ . The dynamical age is equal toτdyn = L/VCO ∼ 3 × 104 yr. Inside the CO outflow lobe, several Herbig-Haro objects are found, which
are emission-line nebulosities. Proper motion studies have revealed the motion of such H-H objects.
Cudworth & Herbig (1979) reported HH28 and 29 have 0.15 − 0.20arcsec yr−1 , which corresponds
to 150 − 200km s−1 . In a neighbor of the IR source a stellar jet is found (Mundt & Fried 1983) by
optical emission lines. Its size is equal to 17arcsec = 4 × 1016 cm(D/160pc) and the jet indicatesrather well-collimated shape (the opening angle ' 10◦ ). The relationship between these two outflows(massive molecular outflow and less massive optical jet) are not clear yet.

4.7.1

Magneto-driven Model

To accelerate gas till supersonic speed, the magnetic force is considered to be an important player.
Blandford & Payne (1982) have proposed the “magneto-centrifugal wind model” for the accelerationin two-dimensional axisymmetric configuration.
Figure 4.10 explains the magneto-centrifugal wind mechanism (Fig. 1 of Blandford & Payne 1982).
Consider a disk under the gravity of a central star with mass M∗ . Assume that gas in the disk at theradius r0 rotates with a Kepler speed ofµvφ (r0 ) =

GM∗
r0

¶1/2

.

(4.84)

Consider a gas parcel whose original position was (r0 , z = 0). The effective potential for the particle

88

CHAPTER 4. LOCAL STAR FORMATION PROCESS

2

1.5

1

0.5

0
0

0.5

1

1.5

2

Figure 4.10: Magneto-centrifugal wind model. Figure shows the effective potential for a gas parcel
with conserving the angular rotation speed ω0 of the disk (eq.4.85). The gas parcel is assumed tolocated originally at (r0 , 0). The x- and y-axes mean normalized distance from the central star as
r/r0 and z/r0 . The darker color represents the deeper potential.
rotating with an angular rotation speed of ω0 = vφ (r0 )/r0 is written as
"

GM∗
r0
1
1
φ(r, z) = −− r2 ω02 = −
+
1/2
1/2
2
r0
2
(r2 + z 2 )(r2 + z 2 )
GM∗

µ

r
r0

¶2 #

,

(4.85)

3/2

where ω0 = (GM )1/2 /r0 . (Be careful that ordinary effective potential is calculated for constantangular momentum. However, this is for constant angular rotation speed.) The gas will co-rotatewith the same angular rotation speed as the disk, if the magnetic field is sufficiently strong. Equation(4.85) corresponds the effective potential for the gas co-rotating with the disk at r0 . Figure 4.10
plots the isopotential contour lines for the effective potential [eq.(4.85)]. The contour line passing thepoint (r0 , 0) represents Φ = −3GM/2r0 . The angle between the isopotential line of Φ = −3GM/2r0
and the r-axis is equal to ±60◦ (see problem below). Further, it can be obtained that the effective
potential increases upwardly while it decreases with reaching the central star or radially outwardly.
This leads to an important conclusion that
1. If the angle between the magnetic field line and the r-axis θmag is larger than 60◦ (θmag > 60◦ ),gas parcel moving along the field line must climb the effective potential. It needs extra energy
to depart from the disk.
2. If θmag < 60◦ , gas parcel moving along the field line slides down the effective potential. It can
depart from the disk without extra energy.
Be careful that this is valid if the magnetic field is sufficiently strong. If the Kepler disk is threadedwith the magnetic field lines with a shallow angle, and if the angular momentum is sufficiently
transferred along the magnetic field line, in this case, the gas is accelerated by the extra centrifugal
force and finally escapes from the potential well of the central star. This mechanism is called the“magneto-centrifugal wind” mechanism.

4.7. OUTFLOW

89

Side View
Pole-on View

Figure 4.11: Pole-on (left) and side (right) views of the magnetic field lines, which drive the outflow.
ProblemShow that the angle between the isopotential line which passes the point (r0 , 0) and the r-axis isequal to ±60◦ . Calculate the Taylor expansion of the effective potential near the point (r0 , 0) andshow that to satisfy φ(r0 + ∆r, ∆z) = φ(r0 , 0), (∆z/∆r)2 must be equal to 3.
Angular Momentum Transfer
In the axisymmetric case, the poloidal (z, r) and toroidal (φ) components of the magnetic field B andcurrent j are decoupled with each other. That is, the poloidal (Bp ) and toroidal (Bφ ) magnetic fields
are made by the toroidal (jφ ) and poloidal (jp ) electric currents, respectively. As for the Lorentz force
j × B, the poloidal component comes from jφ Bp or jp Bφ , while the toroidal component does fromj p ×B p . Even if there is no toroidal magnetic field (thus no poloidal electric current), there exists thepoloidal component of the Lorentz force, which acts as a pressure to counter-balance the self-gravity(§4.2). On the other hand, the toroidal component of the Lorentz force appears only the case withthe poloidal electric current and thus toroidal component of magnetic field Bφ . This means thatthe angular momentum is transferred by the magnetic field only when Bφ exists. Equation (B.14)
explains how the angular momentum density ρvφ r is transferred. The left-hand side of equation(B.14) represents the advection of the angular momentum density, while the right-hand side
·

¸

∂Bφ
r 1 ∂
(rBφ )Br +Bz =4π r ∂r
∂z

µ

Bp
4π

¶

· ∇(rBφ ),(4.86)

represents the torque exerted on the gas parcel.
The induction equation of the magnetic field [eq.(B.17)] shows that Bφ is generated from poloidal
magnetic field by the effect of rotational motion vφ . This indicates that the angular momentum istransferred as follows:
1. The toroidal component of magnetic field Bφ is amplified by the dynamo process [equation(B.17)]. The magnetic field lines run like spiral viewing from the top like Figure 4.11.
2. Negative torque works near region attached to the fast-rotating disk, while positive one appears
in the region connected to the radial magnetic field. The angular momentum is transferredfrom inside to outside along the same field line.

90

CHAPTER 4. LOCAL STAR FORMATION PROCESS
3. Angular momentum received by the gas near disk surface leads the excess centrifugal force which
accelerate the gas. The disk inside which is affected by the negative torque loses its angularmomentum and insufficient centrifugal force drives further accretion.

This angular momentum transfer on the same magnetic field has a character of angular momentum
redistribution.
Centrifugal Radius
In a diffuse cloud, the centrifugal force does not play an important role. Consider a gas parcel whose
specific angular momentum is equal to j. When the gas element contracts till the distance r fromthe center, the centrifugal force per mass of
vφ2
j2
= 3,
r
r

(4.87)

works. Contraction stops when the gravitational acceleration GM/r2 is balanced by the centrifugal
force. Its radius (centrifugal radius) is expressed as
rc =

j2
.
GM

(4.88)

Assuming that the disk is near the hydrostatic balance in z-direction, total column density is related
to the volume density on the z = 0 plane asÃσ=

2c2s ρc
πG

!1/2

.

(4.89)

We can show that a nondimensional ratio of specific angular momentum to the mass is equal to theratio of the free-fall time to the rotation period as follow:
cs j
cs ωc r2
ωc
=
=
.
2
GM
Gσπr
(2πGρc )1/2

(4.90)

Many numerical simulations confirm that the contracting disk in the runaway isothermal contraction
phase has some ‘universality.’ The ratio of the free-fall time to the rotation period is approximately
equal toωc
' 0.3,(4.91)(2πGρc )1/2
regardless of the initial conditions (Matsumoto, Nakamura, & Hanawa 1997). Thus, cs j/GM ' 0.3.
Finally, we can see the centrifugal radius asGM
rc = 2
cs

µ

cs j
GM

¶2

' 0.3

GM
,
c2s

(4.92)

increases with time in proportional to the mass M , because gas element with large j contracts later.
As the centrifugal radius increases with time, the launching point of the outflow also expands withtime.

4.8. EVOLUTION TO STAR

91

Figure 4.12: (Left) Relation between density and temperature at the center of the cloud. (Right)
Cooling and heating rates, which are calculated for the cloud center are plotted against the density.

4.7.2

Entrainment Model

Models that a protostar is driving a jet and the jet transfers its momentum to the ambient material
and forms a molecular outflow are called as entrainment models. From hydrodynamical calculations,
De Young (1986) pointed out there are two mechanisms of entrainment. The ejected jet forms a bow
shock and matter hit by the bow shock is accelerated promptly (prompt entrainment). His numericalsimulations predicted that gas within the radius from the jet axis of (1 − 3)× jet radius is entrainedin this mechanism. After the bow shock has passed away, slower mechanism of entrainment works,
in which a turbulent layer is developed around the jet and the Reynolds stress in the layer transfers
the momentum (Stahler 1994). In this mechanism, the ratio of entrained ambient gas to the injectedmass by the jet is estimated as (1/3 − 3) (De Young 1986).

4.8

Evolution to Star

In the proceeding sections we have seen the evolution from gas cloud to stars. However we wererestricted to the isothermal gasses. To understand the change in the temperature, we have to consider the radiation which keeps the gas isothermal. Spherical symmetric radiative hydrodynamical
simulations have been studied to understand the evolutionary path from interstellar cloud to star.
Figure 4.12 (left) illustrates the relation between ρ and T obtained by Masunaga & Inutsuka (2000).
Several characteristic power-laws are seen in the figure. Figure 4.13 plots the density and temperature
distribution.
Isothermal Phase−13−30
Below ρ<∼10 g cm , gas is essentially isothermal T ∝ ρ . This corresponds to points number 1-3
of Figure 4.13. Since internal energy of gas is transferred to the thermal energy of dusts by collisions,the main coolant in this regime is the dust thermal radiation. The cooling rate per mass isΛ = 4κp σT 4 ,(4.93)

where κp (10K) ∼ 0.01cm2 g−1 and σ represent Planck mean absorption coefficient and the Stephan−14−3 is cosmic-ray heating: ²
Boltzman constant. Main heating process for ρ<CR ∼ 2 ×∼10 g cm

92

CHAPTER 4. LOCAL STAR FORMATION PROCESS

Figure 4.13: Dynamical evolution of a spherical cloud studied with RHD simulation. Density (left)
and temperature (right) distributions are plotted against the radius.

4.8. EVOLUTION TO STAR93

10−4 erg s−1 g−1 (Goldsmith & Langer 1978). Balance between these two majors asks that temperature
is constant asµ¶1/4 µ
¶1/2
κp
²CR
)
)
.
(4.94)T ∼ 3K2 × 10−4
0.01cm2 g−1
Heating rate due to dynamical compression,Γ = −pd(1/ρ)' c2s (4πGρ),
dt(4.95)

increases according to the contraction and it balances with the above cooling at the density ρA ∼10−14 g cm−3 (Masunaga, Miyama, & Inutsuka 1998).
The cloud in this phase experiences the dynamical contraction as described in section 4.5. Structure of ρ(r) and vr (r) is well represented by the Larson-Penston self-similar solution.
First Core−13
−3 ρ<10−9 g cm−3 ,After ρ>∼ρA , gas is no more isothermal and becomes adiabatic. Between 10 g cm <∼ ∼gas obeys the γ = 5/3 polytropes. Above the density ρ>∼ρA , the optical depth for the thermal radiation
exceeds unity τ >∼1 and radiative cooling can not compensate the compressional heating. As long as
the temperature is low as T <∼100K, neither rotation nor vibration is excited for H2 molecule. Even
H2 gas behaves like single-atom molecule. Thus γ ' 5/3.
−7.5 g cm−3 , the exponent becomes γ = 7/5, which characterizes that
Between 10−9 g cm−3 <∼10∼ρ<the gas consists of two-atom molecule H2 .
In this phase, relatively large gas pressure supports against the gravity and the cloud becomeshydrostatic (points number 4-6 of Figure 4.13). This is called as “first core” made by the molecularhydrogen. The density structure of the first core is well represented by a polytrope sphere with thespecific heat ratio of γ = 7/5 or the polytropic index n = 2 21 . From equation (C.11) in Appendix
C.1, such a polytrope has a mass-density relation as

Mc1 ∝ ρ1/10,
c(4.96)

where Mc1 and ρc represent, respectively, the mass of the first core and the central density. At the
beginning, the core mass is equal to Mc1 ' 0.01M¯ . As long as the mass increases a factor 3, thecentral density increases 5 orders of magnitude.
Second Collapse−2.5 g cm−3 , the exponent becomes γ = 1.1. In this density range of ρ>10−7.5 g cm−3
For 10−7.5 g cm−3 <
∼ρ<
∼10
∼
the temperature exceeds 103 K. Due to this temperature, the hydrogen molecules begin to dissociate
and the dissociation process absorbs the energy. By this, gas becomes much softer (γ → 1) than thefirst core. This introduces another dynamical contraction, the second collapse (points number 6-7 ofFigure 4.13).

Second Core
After the dissociation is completed, the exponent becomes γ = 5/3. This gas forms a second corecomposed of atomic hydrogen, which is called the second core (points number 7-8 of Figure 4.13).

94

CHAPTER 4. LOCAL STAR FORMATION PROCESS

Figure 4.14: Explanation of nested grid method.
Since the accretion rate [eq.(4.75)] is proportional to c3s or T 3/2 , the accretion rates onto the second
core Ṁ2nd is larger than that of the first core Ṁ1st as¡¢3/2

T /103.5 K
Ṁ2nd' 5 × 103.
Ṁ1st
(T /10K)3/2(4.97)

Thus, the first core disappears quickly and after that the gas begins to accrete onto the second
core, which will be a protostar. However, this is the case of non-rotating, spherical symmetric cloud
collapse. As previously seen, the angular momentum plays a crucial role and forms a disk. The
evolution must be different completely for such a case.

4.9

Example of Numerical Simulation

We have described the cloud run-away collapse and succeeding accretion process In the latter phase,
there is a possibility that the outflow is driven by the effect of magnetic Lorentz force. A restricted
number of numerical simulation can simulate such evolution throughout. The size of the molecularcores ∼ 0.1pc is more than 3 × 104 times as large as the typical size of the first core ∼ 1AU.
Therefore numerical scheme to calculate such process must have a large dynamic-range. In the finitedifference scheme, the spatial dynamic range is restricted by the number of the cells. Although at
least 3 × 105 × 3 × 105 grid points (in 2D) seems necessary to resolve a factor of 3 × 104 , to solve the
MHD equation with the Poisson equation of the self-gravity have not been done. This is done by
the Eularean nested-grid simulation. This method uses a finite number of grid systems with different
spatial resolutions are prepared. Coarse grid covers whole the cloud and is to see a global structure.
Fine grid covers only the central part of the cloud and is to see the fine-scale structure appeared near
the center. Figure 4.14 explains grid cells of the nested grid method. The size of the cells of the n-thlevel grid (Ln) is taken equal to a half of those of the Ln − 1.
A cylindrical cloud with coaxial magnetic field in hydrostatic balance is studied. Rotation vector
ω and the magnetic field B are parallel to the symmetric axis of the cylinder (z-axis). For the4.9. EXAMPLE OF NUMERICAL SIMULATION

(a)

95

(b)

(c)

Figure 4.15: Ωc = 10−13 rad s−1 , Bc = 10µG, ρc = 104 H2 cm−3 , and cs = 190m s−1 (T /10K). Snapshots at the time of 0.6066τff = 1.061Myr represented in different levels are shown: L1 (upper-left), L5(upper-right), and L10 (bottom). The actual size of the frames of L5 (upper-right) and L10 (bottom)
are, respectively, 1/16 and 1/512 smaller than that of L1 (upper-left). Magnetic field lines (dottedlines) and isodensity contours (solid lines) are presented as well as the velocity vectors by arrows.

96

CHAPTER 4. LOCAL STAR FORMATION PROCESS

energy equation, a double-polytropic equation of state is adopted, since the radiation hydrodynamical
calculation (for example Masunaga & Inutsuka 2000) predicts ρ − T relations like shown in Figure
4.12, which consists of a number of power-laws for some density ranges.
p= c2s ρ c2s ρA

(for ρ < ρA ),

³

´
ρ 7/5
ρA

(for ρ > ρA ),(4.98)

where we take ρA = 1010 H2 cm−3 . Figure 4.15 is the model of Ωc = 10−13 rad s−1 , Bc = 10µG, ρc =104 H2 cm−3 , cs = 190m s−1 (T /10K), and scale-length H = 0.341pc(cs /190m s−1 )(ρc /104 H2 cm−3 )−1/2 .
This is a snapshot of 0.6066τff = 1.061Myr. The cylindrical cloud breaks into pieces with a lengthin z-direction equal to the most unstable wavelength of the gravitational instability (Matsumoto,
Nakamura, & Hanawa 1994) as
λmax 'h

2π(1 + α/2 + β)1/2
1/3

0.72 (1 + α/2 + β)

i1/2

− 0.6

cs
.
(4πGρc )1/2

(4.99)The magnetic fields prevent gas motions perpendicular to the field lines B. Further, the centrifugal
force works in the perpendicular direction to ω. A disk forms which running in the perpendicular
direction to B and ω and thus to the z-axis. Figure 4.15 illustrates a snapshot at the time of
0.6066τff = 1.061Myr. Level 1 grid (upper-left) captures the global structure of the contracting disk
perpendicular to the z-axis. Close up view is shown in L5, which has 16 times finer spatial resolution(upper right). Since speed, vz , of infalling gas is larger than the sound speed, a number of shock fronts
are formed near z ' 0.02H and z <∼0.01H. Further inner structure is represented in L10 grid (bottom)
which has 32 times finer resolution than L5. Another discontinuity in density is now forming nearz ∼ 0.0002H. Since gas inside of this discontinuity has the density larger than the critical density
ρ > ρA , in this phase the first core (p.93) begins to be formed. In the isothermal collapse phase, acharacteristic power-law for the Larson-Penston self-similar solution is seen in the radial distributionsof density as ρ(r, z = 0) ∝ r−2 , and Bz (r, z = 0) ∝ r−1 .
After t > 0.6067τff the first core is formed. Radial infall motion is accelerated toward the surfaceof the first core and the rotational motion vφ also increases. This amplification promotes the toroidal
magnetic field, Bφ . Since the torque is proportional to (∇×Bφ eφ )×B p , the amplitude of the torqueincreases after the first core formation.
Structure just after ∼500 yr has passed from Figure 4.15 is shown in Figure 4.15 In τ ∼ 500
yr from the epoch of Figure 4.15, flow is completely changed and outflow is launched. Panel (b),which corresponds to a four-times close-up of panel (a), shows clearly that outflow is ejected from a
region near r ' 1.5 × 10−4 H and z ' 0.5 × 10−4 H. And the outflow vectors and the disk has angngle of approximately 45 deg. The magneto-centrifugal wind model favors a small angle between
the magnetic field and the disk, θmag < 60deg. This model predicts that the angle between the flow
and the disk is also smaller than 60 deg. Numerical results is not inconsistent with this prediction ofmagneto-centrifugal wind model.
To see which force is driving such outflow, amplitudes of three forces are compared: the Lorentz
force F m = (∇ × B)×B/4π, the centrifugal force F c = (ρvφ2 /r)er , and the thermal pressure gradientF t = −∇p. The parallel components to the poloidal magnetic field of the three forces are as
Fm ·Bp|B p |∇ × B)×B B p·,
4π|B p |1 Bp= −· ∇(rBφ )2
8πr2 B p |=(4.100)4.10. EVOLUTION IN THE H-R DIAGRAM(Ustyugova et al. 1999)

97

Fc ·

Bp
Br
= ρvφ2 r ·,|B p ||B p |(4.101)

Fp ·

Bp
Bp
= −∇p ·.
|B p ||B p |(4.102)andFigure 4.16 illustrates the largest force at each grid point. Centrifugal force dominated region (regionC), which is shown by asterisk (*), overlaps the outflow region around point P1 (z = 10−4 H, 2.5 ×
10−4 H). Just radially exterior to this region C, there is another region near point P2 (z = 2 × 10−4 H,
3.5 × 10−4 H) filled with puls signs where the magnetic force is dominated (region M). Comparing
with the flow vectors in the left, the strong outflow coincides with this region M and the above region
C. Therefore, the outflow is driven by the centrifugal force and the Lorentz force (pressure gradientby toroidal magnetic fields).
Difference between models with different magnetic field strength but the same rotation speed isillustrated in Figure 4.17. Magnetic to thermal pressure ratio α is taken 1 [panels (a), and (d)],0.1 [panels (b), and (e)], and 0.01 [panels (c), and (f)]. Panels (a), (b), (c) in the left column aresnapshots when the central density reaches ρA while panels (d), (e), (f) in the right column show
the structure after τ = 4.5 × 10−3 τff = 8000yr passed. In the model with α = 1, the outflow gas
flows through a region whose shape resembles a capital letter U. This is similar to that observed inmodel (α = 1 and Ω = 5) shown in Figure 4.16. Decreasing the magnetic field strength α = 1 → 0.1
and → 0.1, the shape of disk changes from a flat disk to a round ellipsoid. The model with neither
magnetic field nor rotation shows a spherical contracting core. Outflow is also affected by decrease inmagnetic field strength. Globally folded magnetic field lines are seen in the outflow region in panel(e). Smaller-scale structure folding the magnetic field lines develops in the outflow region in panel(f). The strucure seen in the polidal magnetic field is well correlated to the distribution of toroidal
magnetic field strength. That is, a packet of strong Bφ pinch the gas and poloidal magnetic field lines
radially inwardly. This means that the toroidal magnetic field is amplified by the rotation motion andbecomes predominant over the poloidal one in the case with meak initial magnetic field. Thus, thewhoop stress, which is expressed by the first term of the Lorentz force of equation(B.13) and comesfrom the tension of toroidal magnetic field, pinches the plasma contained inside of the loop of Bφ . For
this case, the magnetic force dominant region is widely spread compared with the centrifugal force
dominant region, which is far from the model with strong magnetic field. In this case, the magnetic
pressure gradient plays an important role. Gas of outflow is ejected perpendicularly from the disk
and forms flow similar to the capital letter I. In conclusion, U-type and I-type outflows are generated
for strong magnetic model and weak magnetic model, respectively.

4.10

Evolution in the H-R diagram

4.10.1Main Accretion Phase

Before Deuterium Ignition
Gas accreting onto the protostar generates energy as follows
Lacc =GM∗ Ṁ∗R∗(4.103)where M∗ , R∗ and Ṁ∗ represent the mass and radius of protostar and the mass accretion rate.

98

CHAPTER 4. LOCAL STAR FORMATION PROCESS

(a)

(b)

(c)
P3

P2

P1

Figure 4.16: Continuation of the evolution shown in Fig.4.15. Snapshots at the time of 0.6069τff =1.066Myr just ∼ 500 yr after from Fig.4.15. Panel (a) plots the structure of L10 which is similar to
Fig.4.15 (c). Spherical region inside ∼ 1.2 × 10−3 H has been swept by the outflow. Panel (b) is asnapshot of L12 grid which has 4-times finer spatial resolution than L10 (a). In panel (c), we plottedmagnetic force dominated region with ‘+’, centrifugal force dominated with ‘*’, and thermal pressure
force dominated with blank ’ ’.

4.10. EVOLUTION IN THE H-R DIAGRAM

99

(a)

(d)

(b)

(e)

(c)

(f)

Figure 4.17: Difference in models with the same rotation speed Ω = 2 × 10−14 rad s−1 but different
magnetic field strength ((a) and (d): Bc = 10µG, (b) and (e): 3µG, and (c) and (f): 1µG).

100

CHAPTER 4. LOCAL STAR FORMATION PROCESSThe central temperature increases with time. Finally thermonuclear fusion reaction of Deuterium

2 H(p, γ)3 He begins. Before Deuterium burning begins, the protostar is radiative, that is, the energy

is transported radiatively.
When a fresh gas of ∆m accrets, thermal energy of ∆U = GM∗ ∆m/R∗ increases. Virial theorem(eq.[2.120]) requires the potential energy must decrease (increase in the absolute volume) at∆W = −2∆U ' −GM∗ ∆m,R∗ /2(4.104)to achive a mechanical equilibrium. As a result, the fresh gas must contract to about half the radiusat which the gas first joined the core (Stahler, Shu, & Taam 1980).
When the free-falling fresh gas accrets on the static star, an accretion shock forms. Since the gas
temperature is increased with passing the shock front, temperature of the postshock gas, Tg , is muchhigher than that of the radiation, Tr . Thus, the postshock gas cools very effectively. The postshock
region with Tg > Tr is called radiative relaxation region. The outgoing luminosity at the accretion
shock is much larger than that inside of the radiative relaxation region.
In the case of low mass stars since the Kelvin-Helmholtz contraction time tK−H ' GM∗2 /R∗ L(L represent the luminosity at the base of the radiative relaxation region) is much longer than theacctretion time scale tacc ' M∗ /Ṁ∗ . This gives Lacc ' GM∗ Ṁ∗ /R∗ is much larger than L, which isconsistent with the above statement that in the radiative relaxation region a large amount of accretion
luminosity is radiated away. Since tK−H À tacc , the specific entropy inside the relaxation region isessentially frozen to that when the gas obtained passing through the relaxtaion region. Figure 4.18
taken from Stahler, Shu & Taam (1982) shows the distribution of the specific entropy against theaccumulated mass Mr . The bottom curve corresponds to the state before the nuclear burning begins
when no entropy generation occurs. The temperature increases with mass, bacause the star must becompressed to support an extra mass. After the temperature becomes high enough for Deuterium
burning reaction 2 H(p, γ)3 He as T ∼ 106 K, an extra energy is liberated by the nuclear fusion reaction.
This increases the specific entropy mainly in the offcenter region. Figure 4.18 clearly shows that shell
Deuterium burning occurs at M ∼ 0.025M¯ .
After Deuterium Ignition
One gram of interstellar matter generates energy by the reaction of 2 H(p, γ)3 He as
Ã

δ≡

1
mp µH

!µ

¶

nD
²D→He '
nH

(4.105)where we used an energy per reaction of 2 H(p, γ)3 He as ²D→He = 5.5MeV and a composition of
Deuterium of nD /nH = 2.5 × 10−5 . This gives an estimate of the Deuterium-produced luminosity asÃLD = Ṁ δ = 15L¯

Ṁ
−5
10 M¯ yr−1!

,

(4.106)

for a star with the accretion rate of Ṁ ' 10−5 M¯ yr−1 .
Entropy inversion ∂s/∂r < 0 is unstable for the convection (see §2.7 in p.39). Convection is
generated in the region of ∂s/∂M < 0. Convective motion transfers the energy and the entropy. The
specific entropy peak generated by the Deuterium burning is smoothed out by the convection. The
convection induced with the Deuterium burning brings a fresh gas containing Deuterium into the hot
Deuterium burning region. Thus, although the consent of Deuterium is small and Deuterium burning4.10. EVOLUTION IN THE H-R DIAGRAM

101

Figure 4.18: Specific entropy distribution is plotted against the accumulated mass (Stahler, Shu, &Taam 1982). The bottom curve shows the distribution at t = 2.3 × 104 yr (M∗ = 0.23M¯ ). The top
curve is for t = 3.4 × 104 yr. The dashed curve indicates the entropy distribution if no convection is
generated. While, the solid curves represent that of new distribution achieved with convection.

Figure 4.19: Luminosities in protostars accreting at Ṁ∗ = 10−5 M¯ yr−1 (left) and Ṁ∗ = 10−4 M¯ yr−1(right) (Palla & Staler 1993). Accretion luminosity Lacc , Deuterium burning luminosity LD , maximum
radiatively transferable luminosity Lrad and luminosity at the surface Lsurf are plotted.

102

CHAPTER 4. LOCAL STAR FORMATION PROCESS

could not continue long, Deuterium burning continues actually as long as the convection can mix theouter Deuterium rich fresh gas into the burning region due to the convection.
In Figure 4.19, luminosities in protostars accreting at Ṁ∗ = 10−5 M¯ yr−1 (left) and Ṁ∗ =10−4 M¯ yr−1 (right) are plotted (Palla & Stahler 1993). Lrad represents the maximum luminosity
transferable radiatively for a star with mass M∗ and radius R∗ asµ

Lrad ∼ 1L¯M∗
1M¯¶11/2 µ

R∗
R¯

¶−1/2

,

(4.107)

where the Krammers opacity κ ∝ ρT −3.5 is applied for stars with M <∼5M¯ . Actual luminosity at thesurface of the star4
Lsurf = 4πR∗2 Teff,
(4.108)is larger than Lrad for low mass stars (M∗ < M1 = 2.4M¯ for Ṁ∗ = 10−5 M¯ yr−1 and M1 = 3.4M¯
for Ṁ∗ = 10−4 M¯ yr−1 ). Excess energy is transferred by the convection, Lsurf = Lrad + Lcon . This
figure also shows that except for extremely low-mass object the accretion luminosity is the majorenergy source over Deuterium burning Lacc > LD . Since massive objects with M∗ > M3 satisfies
Lrad ' Lsurf > Lacc , interior of such stars are fully radiative (energy is transferred by the radiation).
M3 = 3.9M¯ for Ṁ∗ = 10−5 M¯ yr−1 and M3 = 5.9M¯ for Ṁ∗ = 10−5 M¯ yr−1 .
Figure 4.20 (left) illustrates how the Deuterium burning in protostars evolves. Accumulating
infalling gas and increasing M∗ , the structure of protostars is changed. (a) in low-mass stars, convection feeds the Deuterium of the accreting matter to the burning center. (b) as increasing themass and thus the temperature, the opacity of the stellar interior decreases. Finally a radiative zone
appears and it prevents from mixing Deuterium. (c) as a consequence, central Deuterium is depleted
and energy generation rate is reduced. Thus the central region becomes radiative. (d) For higher
mass stars, Deuterium ignites as a shell source just outside the radiative central region. Due to thisDeuterium shell burning, which is similar to the Hydrogen shell burning in the red giants, radius of
the protostar expands. Figure 4.20 (right) plots the mass-radius relation of a protostar with accretion
rate of Ṁ∗ = 10−5 M¯ yr−1 . The left open circle represents the beginning of central Deuterium burning (state a). On the second circle, radiative barrier appears (state b) and Deuterium shell burning
begins. Beyond the mass, the radius swells dramatically. Further increase in the mass increases thegravity and the radius shrinks. At the third open circle, Hydrogen begins to burn at the center ofthe protostar.
In the HR diagram, the mass-radius relation for the protostars shown in Figure 4.20(right) corresponds to the dotted line (Fig.4.21). Observed premain-sequence stars: T Tauri stars (M∗ <∼2M¯ )<
5M)
distribute
below
the
dotted
line.
Thus,
this
line
is oftenM
and Herbig Ae/Be stars (2M¯ <∼ ∗∼ ¯called the birth line for stars. In Figure 4.21, evolutionary tracks without mass accretion Ṁ = 0 are
plotted from the birth line. Modulation of the birth line is related to the swell and shrink of radius
shown in Figure 4.20. That is, increase of luminosity around M ' 1M¯ corresponds to the first swellafter Deuterium burning begins (between the first and second open circles). Another increase aroundM ' 2.5M¯ is related to the Deuterium shell burning (after the second open circle).

4.10.2Premain-sequence Evolution

Premain-sequence evolutionary path is understood as follows: Consider a situation that a star withM∗ is left when accretion stops. Figure 4.19 indicates that a less-massive star with M∗ < M1 isfully convective. Thus, this evolves along the convective Hayashi track (Hayashi 1961) down to themain-sequence. In contrast, a massive star with M∗ > M3 (Fig.4.19) is fully radiative and evolves
along the radiative Henyey track (Henyey, LeLevier, & Levee 1955) to the main-sequence. These

4.10. EVOLUTION IN THE H-R DIAGRAM

103

Figure 4.20: Schematic view of Deuterium burning in protostars (left). Evolution of the radius of a
protostar is plotted against its mass (right). Both figures are taken from Palla & Stahler (1990).

104

CHAPTER 4. LOCAL STAR FORMATION PROCESS

Figure 4.21: Premain-sequence evolutionary path in the H-R diagram. (Palla & Staler 1993). Evolution of stars with different masses from 0.6M¯ to 6M¯ are calculated from the birth line (dottedline).
differences are clearly indicated in Figure 4.21. We can see veritical Hayashi track is seen essentiallyfor stars with M <∼2.5M¯ evolve to the upper-left direction along the∼1.5M¯ . Massive stars with M >Henyey track.

Bibliography[1] Alves, J., Lada, C.J., & Lada, E.A. 2001, Internal Structure of a Cold Dark Molecular Cloud
Inferred from Extinction of background Starlight, Nature, 409, 159[2] Alves, J., Lada, C.J., Lada, E.A., Kenyon, S.J., & Phelps, R. 1998, Dust Extinction and Molecular
Cloud Structure: L977, ApJ, 506, 292[3] André, P. 1994, Observations of Protostars and Protostellar Stages in “The Cold Universe,” eds.
T. Montmerle, C.J.Lada, I.F.Mirabel, & J.Tran Thanh Van, Editions Frontières, p.179[4] Binney, J., & Tremaine, S. 1988, Galactic Dynamics, Princeton Univ. Pr. Chapter 6.
[5] Bonnor W. B. 1956, Boyle’s Law and Gravitational Instability, MNRAS, 116, 351[6] Blandford, R. D., & Payne, D. G. 1982, Hydromagnetic Flows from Accretion Discs and theProduction of Radio Jets, MNRAS, 199, 883[7] Chandler, C. J., & Sargent, A. I. 1993, The Small-Scale Structure and Kinematics of B335, ApJ,
414, L29[8] Chandrasekhar, S. 1939, An Introduction to the Study of Stellar Structure, (Univ. of Chicago Pr.;
Chicago)[9] Cudworth, K. M., & Herbig, G. 1979, In Two Large-Proper-Motion Herbig-Haro Objects, AJ, 84,
548[10] De Young, D. S. 1986, Mass Entrainment in Astrophysical Jets, ApJ, 307, 62[11] Ebert R. 1955, er die Verdichtung von H I-Gebieten. Mit 5 Textabbildungen, Z. Astrophys., 37,
217[12] Elmegreen, B.G.; Seiden, P.E. & Elmegreen, D.M. 1989, Spiral Arm Amplitude Variations andPattern Speeds in the Grand Design Galaxies M51, M81, and M100, ApJ, 343, 602[13] Fiege, J.D., & Pudritz, R.E. 2001, Helical Fields and Filamentary Molecular Clouds - I, MNRAS,
311, 85[14] Fridlund, C.V.M., & Liseau, R. 1998, Two Jets from the Protostellar System L1551 IRS5, ApJL,
499, 75[15] Goldsmith, P. F., & Langer, W. D. 1978, Molecular Cooling and Thermal Balance of Dense
Interstellar Clouds, ApJ, 222, 881
105

106

BIBLIOGRAPHY[16] Gordon, M.A., & Burton, W.B. 1976, Carbon Monoxide in the Galaxy. I - The Radial Distribution of CO, H2, and Nucleons, ApJ, 208, 346[17] Haisch, K. E., Lada, E. A., & Lada, C. J. 2001, Disk Frequencies and Lifetimes in Young Clusters,
ApJ, 553, L153
[18] Hanawa, T. 1999, Star Formation, in “Active Universe” in Japanese (Shokabo, Tokyo) ed. by K.
Shibata, J. Fukue, R. Matsumoto, & S. Mineshige[19] Hartmann, L. 1998, Accretion Processes in Star Formation, Cambridge U.Pr.
[20] Hayashi, C. 1961, Stellar Evolution in Early Phases of Gravitational Contraction, PASJ, 13, 450[21] Heithausen, A., Bensch, F., Stutzki, J., Falgarone, E., & Panis, J.F. 1998, The IRAM Key
Project: Small-Scale Structure of Pre-Star Forming Regions Combined Mass Spectra and Scaling
Laws, AAp, 331, L65[22] Henriksen, R., Andre, P., & Bontemps, S. 1997, Time-Dependent Accretion and Ejection Implied
by Pre-Stellar Density Profiles, AAp, 323, 549[23] Henyey, L. G., Lelevier, R. Levee, R. D. 1955, The Early Phases of Stellar Evolution, PASP, 67,154[24] Hirano, N., Kameya, O., Nakayama, M., & Takakubo, K. 1988, Bipolar outflow in B335, ApJL,
327, L69[25] Honma, M., Sofue, Y., & Arimoto, N. 1995, The Molecular Front in Galaxies. II. Galactic-scale
Gas Phase Transition of HI and H2 , AAp, 304, 1[26] Hunter, C. 1977, The Collapse of Unstable Isothermal Spheres, ApJ, 218, 834[27] Itoh, Y. et al. 2000, A Pair of Twisted Jets of Ionized Iron from L 1551 IRS 5, PASJ, 52, 81[28] Kennicutt, R.C. 1998, The Global Schmidt Law in Star-forming Galaxies, ApJ, 498, 541[29] Kennicutt, R.C., Tamblyn, P., & Congdon, C.W. 1994, Past and Future Star Formation in Disk
Galaxies, ApJ, 435, 22[30] Kramer, C., Stutzki, J., Rohrig, R., & Corneliussen, U. 1998, Clump Mass Spectra of Molecular
Clouds, AAp, 329, 249[31] Lada, C.J. 1999, The Formation of Low Mass Stars, in “The Origin of Stars and Planetary
Systems” ed. C.J. Lada & N.D. Kylafis, (Kluwer) p.143[32] Larson, R. B. 1981, Turbulence and Star Formation in Molecular Clouds, MNRAS, 194, 809[33] Loren, R. B. 1989, The Cobwebs of Ophiuchus. I - Strands of 13 CO - The Mass Distribution,
ApJ, 338, 902[34] Masunaga, H. & Inutsuka, S. 2000, A Radiation Hydrodynamic Model for Protostellar Collapse.
II. The Second Collapse and the Birth of a Protostar, ApJ, 531, 350[35] Masunaga, H., Miyama, S. M., & Inutsuka, S. 1998, A Radiation Hydrodynamic Model for
Protostellar Collapse. I. The First Collapse, ApJ, 495, 346

BIBLIOGRAPHY

107[36] Matsumoto, T., Hanawa, T., & Nakamura, F. 1997, Gravitational Contraction of Rotating
Clouds: Formation of Self-similarly Collapsing Disks, ApJ, 478, 569[37] Matsumoto, T., Nakamura, F., Hanawa, T. 1994, Fragmentation of a Magnetized Filamentary
Molecular Cloud Rotating around its Axis, PASJ, 46, 243[38] Mizuno, A., Onishi, T., Hayashi, M., Ohashi, N., Sunada, K., Hasegawa, T., & Fukui, Y. 1994,
Molecular Cloud Condensation as a Tracer of Low Mass Star Formation, Nature, 368, 719[39] Mizuno, A. et al. 1995, Overall Distribution of Dense Molecular Gas and Star Formation in the
Taurus Cloud Complex, ApJL, 445, L161
[40] Mouschovias, T. Ch. 1976a, Nonhomologous contraction and equilibria of self-gravitating, magnetic interstellar clouds embedded in an intercloud medium: Star formation. I Formulation of the
problem and method of solution, ApJ, 206, 753[41] Mouschovias, T. C. 1976b, Nonhomologous Contraction and Equilibria of Self-Gravitating, Magnetic Interstellar Clouds Embedded in an Intercloud Medium: Star Formation. II - Results, ApJ,
207, 141[42] Mouschovias, T. Ch., Paleologou, E. V., & Fiedler, R. A. 1985, The Magnetic Flux Problem andAmbipolar Diffusion during Star Formation - One-Dimensional Collapse. II - Results, ApJ, 291,
772
[43] Motte, F., & André, P. 2001, The Circumstellar Environment of Low-Mass Protostars: A Millimeter Continuum Mapping Survey, AAp, 365, 440[44] Motte, F., André, P., Ward-Thompson, D.,& Bontemps, S. 2001, A SCUBA Survey of the NGC
2068/2071 Protoclusters, AAp, 372, L41[45] Myers, P.C. 1978, A Compilation of Interstellar Gas Properties, ApJ, 225, 380[46] Mundt, R., & Fried, J. W. 1983, Jets from Young Stars, ApJ, 274, 83[47] Nakano, T. 1979, Quasistatic Contraction of Magnetic Protostars due to Magnetic Flux Leakage- Part One - Formulation and an Example, PASJ, 31, 713[48] Nakano, T. Nishi, R., Umebayashi, T. 2002, Mechanism of Magnetic Flux Loss in Molecular
Clouds, ApJ, 572, 199[49] Ogino, S., Tomisaka, K., & Nakamura, F. 1999, Gravitational Collapse of Spherical Interstellar
Clouds PASJ, 51, 637.
[50] Ohashi, N., Hayashi, M., Ho, P.T.P., Momose, & M., Hirano, N. 1996, Possible Infall in the Gas
Disk around L1551 IRS 5, ApJ, 466, 957[51] Ohashi, N. Lee, S. W., Wilner, D. J. & Hayashi, M. 1999, CS Imaging of the Starless Core
L1544: an Envelope with Infall and Rotation, ApJL, 518, L410
[52] Onishi, T., Mizuno, A., Kawamura, A., Ogawa, H., Fukui, Y. 1996, A C18 O Survey of Dense
Cloud Cores in Taurus: Core Properties, ApJ, 468, 815

108

BIBLIOGRAPHY[53] Onishi, T., Mizuno, A., Kawamura, A., Tachihara, K., & Fukui, Y. 2002, A Complete Search for
Dense Cloud Cores in Taurus, ApJ, 575, 950[54] Paleologou, E. V., & Mouschovias, T. Ch. 1983, The Magnetic Flux Problem and Ambipolar
Diffusion during Star Formation - One-Dimensional Collapse. I - Formulation of the Problem andMethod of Solution, ApJ, 275, 838[55] Palla, F. & Stahler, S. W. 1990, The Birthline for Intermediate-Mass Stars, ApJ, 360,47[56] Palla, F. & Stahler, S. W. 1993, The Pre-Main-Sequence Evolution of Intermediate-Mass Stars,
ApJ, 418, 414[57] Reipurth, B., Yu, K.C., Rodrı́guez, L.F., Heathcote, S., & Bally, J. 1999, Multiplicity of the HH
111 Jet Source: it Hubble Space Telescope NICMOS Images and VLA Maps, AAp, 352, L83[58] Saito, M., Sunada, K., Kawabe, R., Kitamura, Y., & Hirano, N. 1999, he Initial Conditions for
Formation of Low-Mass Stars: Kinematics and Density Structure of the Protostellar Envelope in
B335, ApJ, 518, 334[59] Sanders, D.B., Scoville, N.A., & Solomon, P.M. 1985, Giant Molecular Clouds in the Galaxy. II.
Characteristics of Discrete Features, ApJ, 289, 373[60] Schmidt, M. 1959, The Rate of Star Formation, ApJ, 129, 243[61] Shu, F. H. 1978, Self-similar Collapse of Isothermal Spheres and Star Formation, ApJ, 214, 488[62] Shu, F. H. 1983, Ambipolar Diffusion in Self-Gravitating Isothermal Layers, ApJ, 273, 202[63] Shu, F.H., Milione, V., & Roberts, W.W.,Jr. 1973, Nonlinear Gaseous Density Waves and Galactic Shocks, ApJ, 183, 819[64] Snell, R. L., Loren, R. B., & Plambeck, R. L. 1980, Observations of CO in L1551 - Evidence forStellar Wind Driven Shocks, ApJ, 239, 17[65] Spitzer, L. Jr. 1978, Physical Processes in the Interstellar Medium, Wiley, Cahper 13.
[66] Stahler, S. W. 1983, The Birthline for Low-mass Stars, ApJ, 274, 822[67] Stahler, S. W. 1994, The Kinematics of Molecular Outflows, ApJ, 422, 616[68] Stahler, S. W., Shu, F. H., & Taam, R. E. 1980, The Evolution of Protostars. II - The HydrostaticCore, ApJ, 242, 226[69] Stutzki, J. & Guesten, R. 1990, High Spatial Resolution Isotopic CO and CS Observations of
M17 SW - The Clumpy Structure of the Molecular Cloud Core, ApJ, 356, 513[70] Tafalla, M., Mardones, D., Myers, P. C., Caselli, P., Bachiller, R. & Benson, P. J. 1998, L1544:
A Starless Dense Core with Extended Inward Motions, ApJ, 504, 900[71] Tamura, M., Hough, J.H., Hayashi, S.S., 1995, Millimeter Polarimetry of Young Stellar Objects:
Low-Mass Protostars and T Tauri Stars, ApJ, 448, 346[72] Tomisaka, K. 2002, Collapse of Rotating Magnetized Molecular Cloud Cores and Mass Outflows,
ApJ, 575, 306

BIBLIOGRAPHY

109[73] Tomisaka, K., Ikeuchi, S., & Nakamura, T. 1988, The Equilibria and Evolutions of Magnetized,
Rotating, Isothermal Clouds. I - Basic Equations and Numerical Methods, ApJ, 326, 208[74] Tomisaka, K., Ikeuchi, S., & Nakamura, T. 1988, Equilibria and Evolutions of Magnetized, Rotating, Isothermal Clouds. II - The Extreme Case: Nonrotating Cloud, ApJ, 335, 293[75] Tomisaka, K., Ikeuchi, S., & Nakamura, T. The Equilibria and Evolutions of Magnetized, Rotating, Isothermal Clouds. IV - Quasi-Static Evolution, ApJ, 362, 202[76] Toomre, A. 1981, What Amplifies the Spirals, in “The Structure and Evolution of Normal Galaxies” Cambridge U.Pr., 1981, pp.111-136[77] Tothill, N. F. H., White, G. J., Matthews, H. E., McCutcheon, W. H., McCaughrean, M. J.,& Kenworthy, M. A. 2002, The Structure and Evolution of the Lagoon Nebula. I. Submillimeter
Continuum and CO Line Mapping, ApJ, 580, 285[78] Ustyugova, G. V., Koldoba, A. V., Romanova, M. M., Chechetkin, V. M., & Lovelace, R. V.
E. 1999 Magnetocentrifugally Driven Winds: Comparison of MHD Simulations with Theory, ApJ,
516, 221[79] Ward-Thompson, D., Kirk, J.M., Crutcher, R.M., Greaves, J.S., Holland, W.S., André, P. 2000,
First Observations of the Magnetic Field Geometry in Prestellar Cores, ApJL, 537, L135
[80] Weintraub, D.A., Goodman, A.A., & Akeson, R.L. 1999, Polarized Light from Star-Forming
Regions, in Protostars and Planets IV, ed. V.Mannings, A.P.Boss, & S.S. Russell, U. Arizona Pr.
pp.247-271[81] Whitworth, A., & Summers, D. 1985, Self-Similar Condensation of Spherically Symmetric SelfGravitating Isothermal Gas Clouds, MNRAS, 214, 1[82] Wolf, S., Launhardt, R. & Henning, T., 2003, Magnetic Field Evolution in Bok Globules, ApJ,
592, 233[83] Yun, J. L., & Clemens, D. P. 1994, Outflows from Young Stellar Objects in BOK Globules: Maps,
ApJS, 92, 145[84] Zhou, S. 1995, Line Formation in Collapsing Cloud Cores with Rotation and Applications to
B335 and IRAS 16293-2422, ApJ, 442, 685[85] Zhou, S., Evans, N.J.,II, Koempe, C., & Walmsley, C.M. 1993, Evidence for Protostellar Collapsein B335, ApJ, 404, 232