of Importance the Interstellar Medium
Gas has important diagnostic
properties
Role in the star/gas cycle
facilitates ongoing star
formation
repository for element buildup;
integral for chemical evolution
Gas can cool, so its collapse is
dissipational
Hot gas cold gas stars
Star formation cools spiral
disks, leading to arm formation
Gas migrates inwards in the
gravitational potential
Galactic disks are smaller
than dark matter halos
Galaxies have steep density
gradients
Galactic nuclei have high
densities, including massive
black holes
J.M. Lattimer
Doppler eect reveals dynamics of
Galaxy
Abundances show chemical
evolution
Physical conditions can be found
Some emission lines are seen at
cosmological distances
High-redshift absorption lines
reveal galaxy birth & evolution
Can dominate the integrated
spectral energy distribution
Dust: mid-IR to sub-mm
Hot ISM phase (and X-ray
binaries): soft X-rays
HII regions and relativistic
plasmas: cm-radio
Some emission lines (Ly , [CII])
are major coolants
AST 346, Galaxies, Part 3
Activity in the Interstellar Medium
ISM energized by stars
UV light ionizes atoms, dissociates
molecules
photoelectric eect heats gas
SN shocks heat, ionize and accelerate
gas
ISM is inhomogeneous with phases
Hot/warm/cold phases with similar
pressures (P = nkT 1 eV cm3 )
Cloud and intercloud media with huge
density contrasts (102 105 )
Mass and metallicity exchange between
phases
Cooling: hot warm cold stars
SN accelerate gas and rearrange phases
(bubbles and fountains out of disk)
Tidal encounters and resulting starbursts
create bubbles
cycle gas into halos
convert spirals into ellipticals
J.M. Lattimer
Global distribution of ISM
colder phases conned
closely to plane
hotter and turbulent phases
are thicker
ISM in disk is thin at small
radii and ares at large radii
ISM is locally complex
SN create superbubbles
between bubbles are cold,
dense sheets
Equipartition in the ISM
Energy densities of all three
gas phases, starlight,
magnetic elds and cosmic
rays are each 1 eV cm3 .
AST 346, Galaxies, Part 3
Structures in the Interstellar Medium
HII regions
Reection nebulae
Dark nebulae
Photodissociation regions
Supernova remnants
J.M. Lattimer
AST 346, Galaxies, Part 3
HII Regions
Ionized H regions formed by O
and B0-B1 stars with an
abundance of photons with
< 912.
A
nH 10 104 cm3
T 104 K
NASA
Total mass 5 107 M
R 0.5 10 pc
Optical spectra dominated by H
and He recombination and [OII],
[OIII] and [NII] lines.
Strong sources of free-free radio
emission and thermal emission
from warm dust.
Signposts of massive star
formation.
Richard Crisp
J.M. Lattimer
AST 346, Galaxies, Part 3
Reection Nebulae
Bluish dusty nebulae
reect light of nearby
stars later than B1.
nH 103 cm3
Spectrum similar to
illuminating star.
Often seen with HII
regions; both are
diuse nebulae.
NASA
Some dust thermal
emission.
Gas from star
formation, a chance
encounter, or ejecta of
late type stars.
Gary Stephens
J.M. Lattimer
AST 346, Galaxies, Part 3
Dark Nebulae
Made visible by
absence of stars or
when backlighted.
NASA
R 0.501 100 pc
They become bright in
the far-infrared.
J.M. Lattimer
AST 346, Galaxies, Part 3
Photodissociation Regions
Predominately neutral regions in
which penetrating far-UV (6
13.6 eV) radiation dissociates
and ionizes molecules and heats
the gas through the
photo-electric eect.
Bright in IR dust continuum and
atomic and molecular line
emission.
Speck et al. PASP 115, 170 (2003)
Includes neutral atomic gas and
gas in molecular clouds outside
their dense cores.
Typical examples are the gas at
the boundary of a giant
molecular cloud or within
planetary nebulae.
Dominate the sky in the infrared.
J.M. Lattimer
AST 346, Galaxies, Part 3
Supernova Remnants
Ejected material shocks
surrounding ISM
T 106 K
Spectrum is that of a
high-velocity shock.
NASA
Prominent sources of
synchotron radiation and
X-radiation.
Can be compact or wispy.
NASA
Joseph d. Schulman
J.M. Lattimer
AST 346, Galaxies, Part 3
Components of the Interstellar Medium
Intercloud/cloud mass is 1/1; intercloud/cloud volume is 49/1
J.M. Lattimer
AST 346, Galaxies, Part 3
Components of the Interstellar Medium
Neutral atomic gas Dominated by 21 cm emission. Can be in cold
neutral diuse HI clouds (nH 25 cm3 , T 80 K) and warm
intercloud gas (nH 0.3 cm3 , T 8000 K) mixed with ionized
gas. Completely absorbs starlight with > 912 (Lyman edge).
A
Ionized gas Traceable through dispersed pulsar signals. optical and
UV ionic absorption lines, and H recombination line emission.
Most H emission comes from HII regions, but most mass is in
diuse warm ionized medium. Has a complex structure including
laments of up to 1 kpc in length. Source of ionization is uncertain.
Molecular gas Dominated by dense giant molecular clouds of average
size 40 pc, mass 4 104 M , density nH2 200 cm3 and
temperature 10 K, traceable by J = 1 0 CO emission at 2.6 mm.
Smaller between spiral arms. Often surrounded by neutral gas
forming complexes to 100 pc and 107 M . Have high turbulent
pressures but are self-gravitating. The site of active star formation,
stable for about 30 million years. Many rotational lines from over
200 molecules seen. Show structure on all scales, including dense
(104 cm3 ) cores of 1 pc and 10 100 M . H2 /CO 104 .
J.M. Lattimer
AST 346, Galaxies, Part 3
Components of the Interstellar Medium
Coronal gas Hot (106 K) intercloud medium traceable through UV
absorption lines (CIV, SVI, NV, OVI). Emit continuum and line
radiation in far UV and X-rays. Fills most of halo and some of the
disk. Gas heated by stellar winds and supernovae; forms bubbles (in
which the Sun is found) and super bubbles from OB associations
which pump the gas into the halo; it then cools into clouds and rains
back down into the disk.
Interstellar dust Responsible for most extinction, reddening,
scattering and polarization. Dominates IR continuum emission.
Typical sizes of 0.1m, size distribution n(a) a3.5 . Contains half
the mass of heavy elements and 1% of total gas mass. Larger grains
are in radiative equilibrium at 15 K with the stellar radiation eld,
but up to 75 K near massive stars.
Large interstellar molecules Visible at mid-IR in broad emission.
Dominated by polycyclic aromatic hydrocarbons (PAH) materials
containing some 50 C atoms, with densities of 107 nH and locking
up 10% of C. Diuse interstellar bands, of which more than 200 are
known, are attributed to large unsaturated carbon chains. Seem to
be the extension of grains into molecular domain; extra-solar nano
diamonds and silicates have been extracted from meteorites.
J.M. Lattimer
AST 346, Galaxies, Part 3
Observational Considerations
Emission Measure (EM) = < n2 > dz (pc cm6 ), proportional to
surface brightness
Column Density (N) = < n > dz (cm2 ), proportional to
absorption
In an ionized gas, ne is the relevant density.
The ISM is highly opaque in EUV (13.6 100) eV, partially transparent
in soft X-rays ( 0.6 eV), completely transparent by 2 keV.
From www.astro.virginia.edu/class/whittle/astr553
J.M. Lattimer
AST 346, Galaxies, Part 3
Gas Distribution
Molecular gas peaks at 4.5 kpc
Disk of molecular gas is very thin,
thickness 75 pc.
Atomic gas is more uniform.
Hole at center of Galaxy except for
nuclear ring.
Mass of HI gas is about 5 times that
of H2 gas.
Atomic gas has a thickness about
200 pc inside the Sun, and ares to
about 1kpc in outer Galaxy.
The outer disk is warped.
neutral H
J.M. Lattimer
AST 346, Galaxies, Part 3
Gas and Dust Budget in the ISM
i
Mg 1.75 103 M
kpc2 yr1
g 8 106 M kpc2
g 5 Gyr
dust/gas ejecta 1.5%
J.M. Lattimer
AST 346, Galaxies, Part 3
Dust
Dust absorbs about half the Milky
Way stars optical and UV emission.
T 100 K
The typical grain size is 0.1m.
thermal dust
emission
T 30 K
Composition are magnesium and
iron silicates and soot and graphite.
In dense clouds, dust has ice (H2 O,
CO2 , CH4 , NH3 ) mantles.
About 10% of dust is in small (100
carbon atom) PAH particles; emit
from 3 30m.
About 1% of the gas mass is in dust;
half of heavy element mass in dust.
About 1 grain per 1012 H atoms.
Dust lifetime 0.5 Gyr.
Heavy element depletion
in gas phase of ISM
Savage, ARAA 34, 279 (1996)
Molecules like H2 primarily form on
dust surfaces, 108 times faster
than in gas phase.
J.M. Lattimer
AST 346, Galaxies, Part 3
Dust Composition and Condensation Temperatures
C/O = 0.55
Lattimer & Grossman, Moon & Planets 19, 169 (1978)
C/O = 1.2
J.M. Lattimer
AST 346, Galaxies, Part 3
Energy Sources in the ISM
Radiation Fields
Magnetic Fields
Cosmic Rays
Kinetic Energy
ISM is an open system and needs a
continuous energy supply.
The gas layer will cool and dissipate by
random motions of clouds without energy.
In HII regions, the recombination time is a
few thousand years.
In warm ionized gas, it is 2 Myr.
J.M. Lattimer
AST 346, Galaxies, Part 3
Radiation Fields
Coronal gas in
halo, SNRs and
intergalactic medium
J.M. Lattimer
AST 346, Galaxies, Part 3
Magnetic Fields
Important energy and
pressure source.
5G near Sun, 8G
around R = 4 kpc.
Traced by synchotron
emission, extragalactic
and pulsar rotation
measures.
Circular uniform
component of 1.5 G with
two eld reversals inside
and outside solar circle.
Han
Non-uniform eld
connected to
superbubbles and shocks.
In dense clouds,
B 30G.
J.M. Lattimer
AST 346, Galaxies, Part 3
Cosmic Rays
High energy ( 100
MeV/b) particles.
Mostly p (85%) and He
(15%), but also electrons
and heavy nuclei up to U.
Contribute 2 eV cm3 .
Solar have energies up to
1010 eV, ejected during
solar ares.
Galactic 1010 1015 eV.
Galactic may originate
from SN, 10% of KE of
ejecta.
Extragalactic < 10
18
eV.
F ( p ) p q
(p ) p
Axford & Ip, A&A 149, 7 (1985)
Extragalactic may
originate from
supermassive black holes.
J.M. Lattimer
AST 346, Galaxies, Part 3
Kinetic Energy
Winds from early-type stars and supernova explosions
Mechanical energy is just 0.5% of stellar radiation.
Turbulent energy is 6 1051 erg kpc1 near Sun.
Provides support against gravity for HI gas in Galactic plane.
HI gas has ordered ows of 5 km s1 .
Expanding shells from stars and superbubbles from OB associations
sweep up and compress ISM, which becomes unstable to
Rayleigh-Taylor and Kelvin-Helmholtz instabilities which create
turbulence.
Kinetic energy decays through shock waves when clouds collide,
producing line radiation and plasma waves which heat gas.
Turbulence in molecular clouds supports them against gravitational
collapse.
J.M. Lattimer
AST 346, Galaxies, Part 3
Physics and Chemistry of the ISM
The ISM is far from thermodynamic equilibrium.
The velocity distribution of gas well described by a single
temperature.
The excitation, ionization and molecular composition have dierent
characteristic temperatures.
Collisions cannot compete with fast radiative decay rates of atoms
and molecules.
Cosmic rays and a diluted EUV-FUV stellar radiation eld keep
chemical compositions from equilibrium.
The large scale velocity eld greatly inuenced by turbulence.
The level populations, ionization degree, chemical composition and
temperature are determined by heating and cooling rates.
In some environments shocks are important.
Dust grains and large molecules have to be specially treated.
J.M. Lattimer
AST 346, Galaxies, Part 3
Spectroscopy
H, He have lowest-lying transition energies a large
fraction of ionization energies, FUV.
Multi-electron atoms have much smaller transition energies, visible-UV.
Radicals and ions have unpaired electrons in low-lying states, visible.
Molecular vibrational levels lowered by me /M , mid-IR.
Molecular rotational levels lowered by me /M , mm.
J.M. Lattimer
AST 346, Galaxies, Part 3
Hydrogen Atom
=1
2
3
= 1
4
SELECTION RULES
angular momentum
conservation J = 0, 1
electron spin only changed
by magnetic eld
dipole has odd parity
electric dipole allowed
= 1; n arbitrary
magnetic dipole forbidden
= 0; n = 0
elec. quadrupole forbidden
= 0, 2; n arbitrary
J.M. Lattimer
AST 346, Galaxies, Part 3
Statistical Equilibrium in the ISM
Transition strengths are usually expressed in terms of the oscillator
strength fji or the Einstein coecients Aij and Bij and Bji . These are
related by
2hji
4 e
Bij /Bji = gj /gi , Aji = 2 Bij , Bij =
fji .
c
hji me c
The oscillator strength fji is the eective number of classical oscillators
involved in the transition. For an electric dipole
8 2 eme
ji 2
ji
3h
which is about 1 for the strongest allowed transitions. In this case,
Aji 107 s1 . For forbidden transitions, Aji 0.1 s1 .
fji =
Consider a two-level ( , u ) atom in statistical equilibrium in a gas with a
radiation eld whose intensity at the transition energy is Ju . The rate of
collisional excitations and de-excitations are u and u . In statistical
equilibrium, including absorption and stimulated emission,
n (n
u
+ B u Ju ) = nu (nu + Bu Ju ) + nu Au .
J.M. Lattimer
AST 346, Galaxies, Part 3
Cooling Rates
L
= nu Au hu .
V
Two-level system in statistical
equilibrium in the optically thin limit
r=
n n
u
where ncr = Au /u . Collisions
dominate when n >> ncr .
For a species with abundance
A = (nu + n )/n = n (1 + n/ncr )/n,
= nu nu + nu Au
Detailed balance:
u
= u
4
u =
gu hu
e
g
2kT
/kT
r = nu Au hu =
3/2
u (v )v 3 e v
n u
u /u
nu
=
=
n
nu + Au
1 + ncr /n
2
When n << ncr , n >> ncr ,
/2dT
dv
r An2 u hu ,
0
u (v ) v = u T (1+)/2
+
(M+M,M+e ,M +e ) = T
Anncr u hu
.
1 + ncr /n + u /u
(1/2,0,1/2)
J.M. Lattimer
I=
1
2
r Au nu hu
z
n2 (z )dz nN ,
0
AST 346, Galaxies, Part .
Cooling 3
Nu Rates
In the case matter is not optically thin, absorption and stimulated
emission must be considered.
n (n
u
+ B u Ju ) = nu (nu + Bu Ju ) + nu Au .
Suppose (u ) is an escape probability; photons produced locally are only
absorbed locally. Also suppose the local optical depth is global. Then
(n B
u
nu Bu )Ju = nu (1 (u ))Au
n n
u
= nu nu + nu Au (u )
r = nu Au hu =
Anncr (u ) u hu
.
1 + ncr (u )/n + u /u
In the two limits n << ncr , n >> ncr we recover the same results as in
the optically thin case.
J.M. Lattimer
AST 346, Galaxies, Part 3
Cooling Rates
ions
Ly
Atomic
degree of
ionization
C+ ,Si+ ,Fe+
Molecular
Log10
J.M. Lattimer
AST 346, Galaxies, Part 3
Cooling Rates
L
= n2 (T ),
V
T
n(T )
(T ) T
tcool T /n
V
NH = 1019 cm2
tcool
J.M. Lattimer
AST 346, Galaxies, Part 3
Heating Processes
Processes that couple radiation to gas dominate
Photoionization The e gains kinetic energy from the photon.
HII regions: H ionization.
HI regions: photoionization of large molecules and small dust
particles (photoelectric eect).
Molecular regions: photodissociation of molecules (H2 ). Collisional
de-excitation also heats.
Dust-gas heating If dust is warmer than gas (protostar envelope).
Cosmic-ray heating e from ionized gases gain kinetic energy and
can lead to secondary ionizations.
X-ray heating e from ionization can lead to secondary ionizations.
Both galactic and extragalactic contributions.
Turbulent heating Viscous heating from motions.
Ambipolar diusion heating Ions and e develop small dierential
drift velocities from counterplay between magnetic elds and gravity,
leading to frictional heating.
Gravitational heating through compression. Thermal energy nT ,
n1/2 , heating rate n3/2 T .
J.M. Lattimer
AST 346, Galaxies, Part 3
Heating and Cooling Processes
T = 100 K
Diuse hot gas
Heated by supernova shocks, and cools
within 104 105 yr, condensing into cooler
clouds. Near the Sun, a supernova shock
passes every 15 Myr. Heating is primarily
by photoionization of HI. Cooled primarily
T = 8000 K
by recombination at highest temperatures
Diuse warm gas
Photoelectric heating domnates from stellar
UV on small grains. Cooled primarily by e
collisions which excite low lying electronic
states of trace ionized species. Mostly these
are forbidden transitions. In neutral regions,
T = 10 K
cosmic rays and X-rays dominate.
Cool molecular gas
Cosmic ray ionization heating is important
for gas; grains are heated by infrared
photons. Cooling occurs through rotational
transitions of molecules, such as CO.
J.M. Lattimer
AST 346, Galaxies, Part 3
Jeans Mass
Consider a uniform density, isothermal cloud. Potential energy
3
GM
dM =
r
5
=
4
3
1/3
G M5/3
Kinetic energy
3No kT
3No kT
dM =
M
2
2
Virial Theorem in equilibrium: = 2U
U=
M=
3/2
5No kT
G
3
4
Jeans length
J cs
=
G
No kT
G
2.5
T
10 K
100 cm3
pc
n
Jeans mass
MJ
3
=
6J
2
30
3/2
4
M
3
J.M. Lattimer
39
T
10 K
3/2
AST 346, Galaxies, Part 3
100 cm3
M.
n
Pressure-Bounded Stable Mass
Hydrostatic equilibrium can be written, equivalent to Virial Theorem,
1
1 G M(r )
d M(r ) = d
3
r
3
Now suppose the uniform density isothermal cloud of mass M is
bounded by the ISM at a xed pressure po at the radius R .
Vdp =
M
=
0
G M(r )
dM = 3
r
For a constant density gas, M(r ) r
M
0
p
d M 4 po R 3 .
3
3G M2
3No kT M 3G M2
= 3No kT M 4 po R 3 , po =
5R
4 R 3
20 R 4
For small R , po < 0. For large R , po > 0 but tends to 0 for R .
There must be a maximum of po for the radius Rm , ( po / R )M = 0:
=
Rm =
4G M
,
15No kT
M=
po =
15No kT
4G
3/2
J.M. Lattimer
3No kT M
16
15No kT
4G M
3
135
=
4
4 3
AST 346, Galaxies, Part 3
5
MJ
4
3
.
Gravitational Stability
If the cooling is ecient, a cloud that exceeds the Jeans mass will
collapse quickly, on a free-fall time
t = 1/
G
1/2
108 nH
yr.
If the temperature doesnt increase with collapse, the Jeans mass
becomes smaller and the collapsing cloud will fragment until the
fragments become optically thick. They then heat, becoming protostars.
Isothermal hydrostatic equilibrium: p = No kT /.
@
dp (r )/dr = G M(r )(r )/r 2 ,
@
d M(r )/dr = 4(r )r 2
@
@
Let x = r 4 G o /(No kT )
@
y
and y = /o ,
=@
2x
2
@
2
2
y (y ) /y + 2y /x + y = 0.
@
@
Note:
@
x 0, y 1 x 2 /6,
@
x , y 2x 2
J.M. Lattimer
AST 346, Galaxies, Part 3
Bonner-Ebert Sphere
(No kT /)4
,
G 3 M2
G M
R (Xcr ) = 0.411
.
No kT
po (Xcr ) = 1.40
Mass within a radius X
M(X ) = (4o )1/2
3/2
No kT
G
I (X ),
Compare to
X
yx 2 dx
I (X ) =
3
16
0
15 3
4
3.15,
4/15
p0 (X ) = No kT o y (X )/
4
=
R (X ) =
(No kT /) 2
I (X )y (X )
4 G 3 M2
G M X
No kT
X=
4 G o
No kT I (X )
Xcr
Find ( po / X )Xcr = 0: Xcr = 6.5
J.M. Lattimer
AST 346, Galaxies, Part 3
0.267.
Observing the Density Prole of a Cloud
Extinction is proportional to column
density.
Z
Z
N (R ) = Z n(r )dz = 2 0 n(r )dz
2
R 2 + z 2 = r 2,
R 2 + Z 2 = Rm
'$
Case 1: Isothermal singular truncated
sphere
n(r )/no = 2/r 2 ,
r Rm
N (R ) =
4no
tan1
R
2
Rm R 2
R
Case 2: Modied isothermal sphere
n(r )/no = 2/(2 + r 2 )
N (R ) =
4no
tan1
2 + R2
Zq E
z
R r
t
'
Rm
&%
2
Rm R 2
2 + R2
J.M. Lattimer
AST 346, Galaxies, Part 3
Bonner-Ebert Sphere and B 68
For dense clouds, absorption is sensitive
to density. Extinction is wavelength
dependent. A () NH .
V
IR
J.M. Lattimer
AST 346, Galaxies, Part 3
Additional Pressure Support
MBE =
1.4
po
1/2
(No kT /)2
G 3/2
5.9
Like the Jeans mass for cloud
conditions, additional pressure needed
to stabilize massive clouds.
turbulence, MJ v 3
Consistent with observation that
more massive stars are formed in
warmer, more turbulent clouds.
1012 erg cm3
po
1/2
2
T
10 K
M
magnetism dominates at beginning,
they always dominate. Collapse occurs
B 2R 3
3G M2
>
,
5R
3
3/2
5
B3
M>
G
48 2 2
magnetic elds
B
10 G
R
0.1 pc
2
M.
3No kT 3G M2 (BR 2 )2
+
4 R po =
5R
3R
Fields generally not frozen due to
small ionization fraction ( 107 ).
If eld is frozen in matter, magnetic
Ions frozen but neutrals not, leading to
ux BR 2 is conserved and magnetic
drifts. Ambipolar diusion timescale:
term has the same R dependence as
3
gravity term. Thus, if either gravity or
J.M. Lattimer
tad 2 1013 (nion /nneutral ) yr.
AST 346, Galaxies, Part 3
Dynamical Evolution
Euler Equations of Motion in Spherical Symmetry
u
1 p
G M(r )
u
+u
+
+
= 0,
t
r
r
r2
1 2
+2
r u = 0,
t
r r
M
M(r )
M
+u
= 0,
= 4 r 2 .
t
r
r
2
For p = cs , use X = r /[cs ( t )] as a self-similar variable. We will see
that t < 0 corrsponds to times early in the collapse, t > 0 to late times,
and t = 0 is a catastrophe point. Also
u (r , t ) = cs V (X ),
(r , t ) =
D (X )
,
4 G ( t )2
M(r , t ) =
3
cs ( t )
m(X ).
G
Find m = m (X V ) m A(X ), m = X 2 D m = DX 2 A. Other two
Euler equations are
D
V A = DA ,
D
J.M. Lattimer
D
2A
AV =
D
X
AST 346, Galaxies, Part 3
Self-Similar Solutions
For the static case, V = 0, D = 2/X 2 and m = 2X . This corresponds to
the singular isothermal solution.
For t < 0, a solution is the homologous one
D = 2/3,
V = (2/3)X ,
m = (2/9)X 3 .
For t 0+ , X :
D (D 2)
2
.
,
D=
2
X
X
Asymptotically D X 2 , V 2 , m X with 2; V 0 (infall).
X
V =D
Now consider what happens for t +, or X 0. The asymptotic
solutions of the Euler equations are equivalent to free-fall:
m0
2m0
,
V =
,
m0 = (X 2 DV )X 0 .
2X 3
X
The value of mo = 0.975 can be found by integration. There is a singular
point in the ow, however, at the critical point A(Xc )2 = 1, or
Xc V (Xc ) = 1. The aymptotic (X ) solution V = 0 is satised at
X = 1, so Xc = 1 and V (X ) = 0 for X 1.
D
J.M. Lattimer
AST 346, Galaxies, Part 3
Self-Similar Isothermal Case
3
mo cs 1/2 3/2
t
r
32 2 G 2
3
u (r , t ) = 2mo cs t 1/2 r 1/2
(r , t ) =
=
M (t ) =
2G M (t )
r
3
mo cs
t
G
3
mo cs
M = 4 r 2 u =
G
2 106 M yr1
m(Xc ) = 2
J.M. Lattimer
AST 346, Galaxies, Part 3
Pressure-less Case
An analytic collapse solution exists when pressure is negligible compared
to gravity, not a bad approximation once collapse begins. Euler equation
u
GM
d 2r
u
du
+u
= 2 =
= 2.
t
r
r
dt
dt
In spherical collapse or expansion, the mass M(r , t ) internal to a mass
point initially at r0 at t = 0 does not change.
1d
dr d 2 r
=
2
dt dt
2 dt
2
dr
dt
= G M
dr
dt
1
r2
= GM
d
dt
1
r
.
Integrating once and then twice,
dr
dt
2
= 2G M
3/2
2G M t = r0
1
1
r
r0
r
r0
1
r
+ tan1
r0
r0
1
r
3
The time to collapse to r = 0 is t = 2 r0 /(8G M). If is assumed
7
uniform, t = 3/(32G 0 ) = 3.7 10 (cm3 /n)1/2 yr.
J.M. Lattimer
AST 346, Galaxies, Part 3
'$
'$
Observing Collapsing Clouds
Consequence of v r 1/2 collapse is a
double peaked line prole. Along a
line-of-sight, two parts of cloud
contribute at each velocity, but farther
point is at greater optical depth and
obscured. Blue peak is higher than red
peak due to higher temperatures near
cloud center. Simple model:
z 2 = r 2 R 2,
= v0 zr 1 (r + a)1/2
=
n0 (r + a)3/2
vr
n
=
0
= n0 (Ri + a)3/2
Zo
&%
&%
2
2
Zo = Ro R 2
vr
n
dI (v )
dv
vr =
constant z
R
r
s
d
Ri
d
R
o
d
d
d
r < Ri
R i < r < Ro
n(z ) e (z ) e (vr (z )v )
2
/ 2
dz
Zo
J.M. Lattimer
AST 346, Galaxies, Part 3
Observing Collapsing Clouds
Consequence of v r (homologous)
collapse is a single peaked line prole.
Along a line-of-sight, only one part of a
cloud contributes at each velocity.
Simple model:
2
2
Zo = Ro R 2
z 2 = r 2 R 2,
vr =
v 0 z
r < Ri
n = n0 (r + a)3/2
vr
n
=
0
= n0 (Ri + a)3/2
dI (v )
dv
Zo
Ri < r < Ro
n(z ) e (z ) e (vr (z )v )
2
/ 2
dz
Zo
J.M. Lattimer
AST 346, Galaxies, Part 3
Stromgren Spheres
High-energy stellar photons
(E > h1 = 13.6 eV) ionizes
hydrogen in a spherical region
surrounding the star. Ionization
continues and a balance with
recombinations eventually occurs.
Recombinations go into both excited
states and the ground state. A
recombination directly into the
ground state can re-ionize. Since not
all recombinations go directly to the
ground state, the ones that do cant
aect the ultimate equilibrium since
they are eventually used up. We
can make the on-the-spot
approximation that ionizing photons
produced by direct recombination to
the gound state are absorbed where
they are created and can be ignored.
J.M. Lattimer
H = 3 1019 cm2 is the average
ionization coecient.
A (Te ) 4 1013 cm3 s1 is the
total recombination coecient.
B (Te ) 2.6 1013 cm3 s1 is
the net recombination coecient into
excited states.
The ux of ionizing photons is
F d =
F =
1
2
R
=2
r
1
1
2 1012
L ( )
d
4 r 2 h
2 2
1
d
c 2 e h/kT 1
.5 pc
r
2
cm2 s1
The last equality is for an O4 star.
AST 346, Galaxies, Part 3
Ionization Equilibrium
At a distance r from the star, the optical depth is
r
(, r ) =
nH (r )H ( )dr .
0
The ionization balance condition is
F e H d = ne np B (Te ).
nH 0
1
Let x = ne /n = np /n , 1 x = nH 0 /n , and
F H d .
H F =
1
Assuming is not sensitive to ,
n B e
1x
,
2
x
H F
The production rate of ionizing photons from an O5 star is
NLy =
1
L ( )
d 5 1049 s1
h
J.M. Lattimer
AST 346, Galaxies, Part 3
Stromgren Spheres
For an optically thick nebula, the size Rs of the sphere is determined
when the total recombination rate in the sphere equals the ionizing
photon luminosity. For constant density and x 1, multiply both sides of
balance equation by 4 r 2 dr and integrate:
e d d = NLy = 4 n2 B R3 /3
s
(L ( )/h )
1
0
Rs =1.2 103 cm3 /n
Ms =
1x =
NLy
N o n B
n B
H F
155
2/3
NLy /5 1049 s1
NLy
5 1049 s1
4 104
n
3 cm3
10
3
10 cm
n
1/3
pc
3
M
5 1049 s1
NLy
The structure within the Stromgren sphere: use z = r /Rs .
1/3
s = nH Rs 103 n/103 cm3
NLy /5 1049 s1
1x
3z 2 e
=
,
x2
s
= = ln[1 z 3 ],
J.M. Lattimer
d
= (1 x )s
dz
3z 2 x 2
(1 x ) =
s (1 z 3 )
AST 346, Galaxies, Part 3
2
r
1 pc
1/3
,
Ionization Front
The ionization front terminates the HII region and is very narrow as 1 x
quickly becomes unity. The thickness of the front is about 1 mean free
path for the ionizing photons. Use x = 1/2:
(x = 1/2) = (nH H )1 =
2
n H
2 103
103 cm3
n
(r )
1
=
Rs
(1 x )nH Rs
1
3
Rs
r
2
.
J.M. Lattimer
AST 346, Galaxies, Part 3
pc
The Eect of Helium
Element/Ion
H H+
He He+
He+ He++
Ionization Potential
h1 = 13.6 eV
h2 = 24.6 eV
h3 = 54.4 eV
Although He is 10% as abundant as H, it
must be included for hotter stars since the
cross section for ionization of neutral He is
10 times larger than for neutral H at the
threshold energy (25 eV). But the
ionization energy of He+ is so large that
only the hottest stars can ionize it.
> O6
T < 40, 000 K
< O6
W-R
40, 000 K < T < 105 K
J.M. Lattimer
AST 346, Galaxies, Part 3
T > 105 K
Treatment of Helium Ionization
Competition for stellar ionizing photons:
nH 0 H
y=
.
nH 0 H + nHe 0 He
Recombinations into He0 excited states and ground state make photons
capable of ionizing H. Use on-the-spot approximation for H and He.
nHe 0
4 r 2
L
nH 0
H e d + ynHe + ne (H ,A H ,B ) = nH + ne H ,B
2
4 r 1 h
L
He e d + (1 y )nHe + ne (He ,A He ,B ) = nHe + ne H ,A
2 h
d /dr =nH 0 H ,
ne
nH +
d /dr =nH 0 H + nHe 0 He ,
2
Rs
Rs < r < R2
nH ,
ne = nH + + nHe +
L d
4 3
= N2 =
R nHe + ne He ,B ,
h
32
3 NLy
2
4 nH H ,B
1/3
,
R3
s
R3
2
J.M. Lattimer
1
nH + nHe ,
r > R2
L d
4 3
= NLy =
R nH + ne H ,B
h
3s
NLy He ,B nHe nH + nHe
N2 H ,B nH
nH
AST 346, Galaxies, Part 3
R2 < Rs
Structure of H/He Stromgren Spheres
J.M. Lattimer
AST 346, Galaxies, Part 3
Planetary Nebulae
H0
J.M. Lattimer
AST 346, Galaxies, Part 3