Chapter 2 23
2.2 Theoretial models for the Loal Bubble
2.2.1 Collisional Ionization Equilibrium Plasmas
From an historial perspetive, the rst type of model proposed to explain the
available EUV and soft x-ray data, was a model in whih the ionization ours
by ollisional proesses and in equilibrium. This type of models are ontrolled by
ollisional ionization equilibrium (CIE) physis. They are also known as oronal
plasmas. Adetaileddisussionofthesemattersisfarbeyondthesopesofthiswork.
However, averyshortdesriptionofthephysialonditions,inwhihtheCIEplasma
holds, is given. For a more detailed physial desription of these models, and for
someof their appliationssee the artilesby Böhringer [121℄,Mewe [122℄, Kallman
and Palmeri [123℄ and Kaastra[124℄, respetively.
Collisionalionizationours,whenaneletronollideswithanionwitharelative
kineti energy greater than the ionization potential of that ion. As a result of
this physial proess, an amount of energy equal to the ionization potential of the
ion is removed from the eletron plasma and transfered as kineti energy to the
removed eletron. Therefore resulting in a net ooling of the eletron plasma. In
the oronal approximation it is assumed, that the eletrons and ions have reahed
thethermalisationstate, meaningthat they havereahed the sametemperature. It
istherefore assumed that a Maxwellian distribution desribestheir energies.
The ollisionalionizationproess is then balaned by areverse proess of equal
magnitude, the reombination of an eletron with an ion. This proess leads to a
radiativeproesswithtwomainhannels. Theradiativeproess ours,beausethe
eletrons are mostly aptured inan exited state of the ion. The eletron then
de-exitesvia a asade through permittedradiativetransitions intothe groundstate.
Itishoweverpossible,that theaptured eletron exites aninnershell eletron. As
onsequene, the exited ion will deay via a two-step proess, alled di-eletroni
reombination. Oneofthevaleneeletronsradiativelyde-exitesinduingaasade
ofphotons similartoradiativereombination. Athigh enoughtemperatures the
di-eletronireombinationdominates overthe radiative reombination. Sinethe hot
gas is a tenuous one, the photons produed in these reombination proesses will
havealowprobabilitytointeratwiththeplasma. Therefore, theyleavetheplasma
withoutany futher interation.
Several models, e.g. from Hayakawa [125℄ and Sanders et al. [103℄, using the
assumption of CIE were developed to explain the LB UV and x-ray emission. The
rstseriousattempttoexplaintheavailableEUVandsoftx-raydatawiththemodel
from Cox and Anderson [120℄. In their model a very large supernova blast wave,
propagatingintothe hot low-density omponent of the interstellar (ISM) medium,
was invoked to explain the soft x-ray observations with typial parameters of the
gas: a temperature of
∼ 10 6
K, ap/k B ∼ 10 4
m− 3
K and an extent ranging fromone toseveral hundred p. In this modelthe bubblehad anage of
∼ 10 5
years.Innes and Hartquist [126℄ built a model in whih the H i avity was arved
by very old supernova remnants and superbubbles. On their model a superbubble
reated roughly
4 × 10 6
yr ago was responsible for injeting∼ 10 52
erg into theinterstellarmediumsurroundingtheSolarsystem,whihouldproduetheobserved
softx-ray and extreme ultravioletradiation bakground.
A deade later Edgar and Cox [127℄ investigated the possibility,that the Loal
Bubble had been reated either by a single supernova or by a series of supernovae.
In this way, they have tried tond the physial parameters that ould parametrize
hot bubbleswith atypialradiusof 100 p and with the average ount rates inthe
B and Cbands mathing the observed ones. Withthe right modelthen they ould
proeed inexploringthe featuresof thismodel,suhasthe x-ray and EUVspetra,
the emissivity and O vi olumn densities. In the models developed by Innes and
Hartquist[126℄ andEdgar and Cox[127℄ the age ofthe superbubble ismuh larger.
However, later it was realized by Breitshwerdt [128℄, that these CIE models
produe a high pressurean order of magnitude higherwhih is in onit with
thepresene ofloallouds,whihrequire muhlowerpressure inordertosurvive 5
.
Anotherdiultyofthiskindofmodelresidesonitsveryfoundations,see
Breitshw-erdt [130℄. The assumption, that athree body proess ours in a tenuous and hot
optially thin plasma is not likely to happen in real hot and rareed astrophysial
plasmas. Moreover, the assumption, that the ionization rate by inelasti ollisions
isexatlybalaned by reombinationproesses, requiresathirdbody. Therefore, in
this proess the temperature of the plasma inreases. The other diulty of these
models resides in the assumption, that reombination proesses dominate in CIE
plasmas. Sine radiative reombination is a ooling proess, the plasma will
even-tuallydepartfrom equilibriumand enter in anon-equilibriumionizationphase, see
Kafatos [131℄,Shapiro and Moore [132℄, Shmutzler and Tsharnuter [133℄.
There are two importanttime saleswhihontrolthe evolution of the thermal
plasmatheooling time (
t cool
)and the reombinationtime (t rec
).t cool
is the timethat a plasma needs to onvert all its internal energy into radiative energy. The
typial ooling time for a thermal plasma an be estimated by using the following
equation:
t cool = 3 2
nk B T
L = 3 2
(n e + n p )k B T n e n p Λ = 3
2 n p n e
n p + 1 k B T n e n p Λ =
= 6.57 × 10 4
n e
n p + 1 T 6 n e Λ 22
yr
(2.1)where
L
is the interstellarooling funtion,Λ 22 = Λ/10 −22
denotes the ooling rateinunits of
10 − 22
erg m− 3
s− 1
, and T6
isthe temperaturein units of T/10 6
K.The reombinationtime, dened by
t rec = 1/(n e α total rec )
(2.2)5
There are several louds in the solar neighbourhood. One of the best studied loud is the
loalinterstellarloud(LIC). ForareentreviewontheLoalInterstellarMediumseetheartile
fromRedeld[129℄. ThederivedtemperatureanddensityfortheLICare6300K,and0.33m
− 3
,
respetively. These values imply a
p/k B ∼ 3200
K m− 3
for the LIC. The pressure imbalanebetweentheLB(
p/k B ∼ 10 4
Km− 3
)and theLICwouldmaketheexisteneofthese loudsnotis the time required for a given ion to reombine with a free eletron plasma, in a
plasma with eletron density
n e
and whereα total rec
is total reombination rate for a given ion ata given temperature. Aftert rec
the ion should not be present anymoreinthe plasma, orit shouldhave a redued ionizationstage.
Thesetwotimesalesan beeasilyalulated. UsingthevaluesfortheLB from
Snowdenetal.[116℄,assuminga
Λ = 1.081 × 10 −22
erg m−3
s−1
andn e ≃ n p
,at cool
of 31.6 Myr is obtained. Comparing this value with the time required for the, e.g.,
Ovii toreombine, wenda
t rec
of4.8Myr,withtheα total rec
of1.47 × 10 − 12
m3
s− 1
,a oeient taken from Nahar [134℄. Doing the same exerise for the C vi ion
one nds a
t rec
of 7.2 Myr, using the reombination rates oeients of Nahar and Pradhan[135℄.The ooling time establishes howlong a hot plasmaan exist inCIE atagiven
temperature and eletron density. The time sale required for a plasma to ahieve
aCIE state is given by the followingexpression:
τ
ie= 1
( α coll + α rec ) n e
(2.3)
where
α coll
andα rec
are the ollisional ionization and reombination rates. This equation gives an estimate of the time sale neessary for a given ion speies toahieve a CIE state. It is normally assumed that the CIE ondition for a thermal
plasmais ahieved, when
τ
ie is muh smallerthan thet cool
:τ
ie<< t cool
(2.4)There areurrentlyalarge numberofnumerialodes,that alulatethe
emer-gent spetrum of hot plasmas in CIE and whih are used to t x-ray spetra from
interstellarplasmas. Mostited in the literature are the odes fromRaymond and
Smith[111℄, Mewe etal. [136℄and Kaastraand Mewe [137℄.
The so-alledRaymond-Smithollisionalionizedplasma have been widely used
totlowresolution spetra obtained duringmissionslikeROSAT and others. This
plasmaodehasevolvedtotheAstrophysialPlasmaEmissionCode(APEC),whih
ombines a sophistiated ode along with the APED database, omposed by more
thanone milliontransitionemission lines [138℄.
2.2.2 Non-Equilibrium Ionization Plasmas
A ompletely dierent approah to model the EUV and the soft x-ray bakground
emission was proposed by Breitshwerdt and Shmutzler [92℄. The authors have
inorporated the atomi physis ode developed by Shmutzler and Tsharnuter
[133℄intoanalgorithmwherethe atomiphysisevolvesaording tothedynamis
of the astrophysialsystem in whih the atomi proesses take plae.
The alulationsbyKafatos[131℄andtheShapiroandMoore [132℄havealready
shown, that even bystarting fromanequilibriumondition, non-equilibriumeets
our, leadingtoa dierent ooling funtionfor ahot plasma. In theiralulations
thenon-equilibriumeetshappeninthetemperaturerangefrom3
× 10 5
to5× 10 5
K.i
ii iii
iv
v vi vii
-4.5 -4.0 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0
CarboninCIECarboninCIE
log(IonFration)log(IonFration)
i ii iii iv v vi vii
viii ix
OxygeninCIEOxygeninCIE
i ii iii iv
v vi vii
-4.5 -4.0 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0
4.0 4.5 5.0 5.5 6.0 6.5
log(T)[K℄
log(T)[K℄
log(IonFration)log(IonFration) CarboninNIECarboninNIE
i ii iii iv v vi vii
viii
ix
4.0 4.5 5.0 5.5 6.0 6.5
log(T)[K℄
log(T)[K℄
OxygeninNIEOxygeninNIE
Figure 2.3: Ionization frations of arbon and oxygen. The left top plot shows the
arbon ionization frations assuming a plasma in CIE. The left bottom plot shows the
arbon ionization frations assuming a plasma in a NEI state. Onthe right top plot the
ionization frations ofoxygen inCIEstate are shown. The ionization frations of oxygen
inaNEIonditionareshownonthebottomrightplot. Thesealulationswereperformed
by Shmutzler and Tsharnuter [133 ℄. The NEI ondition allows the presene of an ion
overa large rangeoftemperatures.
In this regime, more highly ionized speies are found, ontributing then with more
strongemission lines.
The rst real attempttosimultaneously alulatethe emissionof an
astrophys-ial optially thin plasma and to take into onsideration the time-dependent
ion-izationstagesdue totheir thermodynamialhistory path, wasdone by Shmutzler
andTsharnuter [133℄with the purpose ofinvestigatingthe emissionlineregions in
Ative GalatiNulei.
isshownthe ionizedarbonand theionized oxygen alulatedin CIEand NEI,
taken fromShmutzlerand Tsharnuter [133℄. Asitan beseen, thesetwodierent
physial onditions result in two distint ionizationstrutures. The top plots show
the alulationsassuming the ondition of CIE, and the bottom ones the ondition
of NEI.
Shmutzler and Tsharnuter [133℄ had alreadyargued that the CIEassumption
for the Galati halo and hot ISM plasmas with temperatures of 10
5
K [139℄ and
10
6
K [140℄, respetively, had a serious aw. No global astrophysial mehanism,
apableof permanently maintainingthis wide temperaturerange is known.
There-fore,itwassuggested,that wewereobservingthesegasesinapointoftheirthermal
history wherethe gas isout of the CIE ondition.
In the Breitshwerdt and Shmutzler [92℄ elegant model, the dynamialhistory
of the physial system is taken into onsideration. This model assumes, that
sev-eral supernovae explosions ourred inthe past, reatinga superbubble, in adense
environment.
The most important onept of this model is the assumption, that during a
phase of fast adiabati expansion, the temperature of the plasma dereases very
fast, leaving all high ionization stages frozen into the plasma. This happens,
be-ause the reombination time is muh longer than the ooling time of the plasma.
Therefore,alargefrationoftheinternalenergyoftheplasmaisstoredinthehighly
ionizedspeies,whihansubsequentlybereleasedbyreombinationradiation. The
dynamis of the plasma also hanges the density of the plasma and the ionization
struture, whih hanges the ooling funtion, and therefore the dynamis of the
plasma itself. With this new lass of model, Breitshwerdt and Shmutzler [92℄
derived a muh lower temperature T
≈ 4.2 × 10 4
K for the LB. From this modelthey alsoinferred an emissionmeasure EM
≈ 5.8 × 10 −2
m−6
p, to whihorre-sponds an eletron density n
e ≈ 2.4 × 10 −2
m−3
along a line of sight of 100 p.From the two physial parameters, temperature and eletron density, a pressure of
p/k
B
= 2nT≈ 2000
K m−3
isthen obtained, whih is four times less than the onederivedfromCIEmodels. ThisthenmakestheveryexisteneofloudsinsidetheLB
mediumpossible. Theestimatedageofthe LBforthismodelisabout
4 × 10 6
years.Theirmodel also predits the loss of some of the hot ISM into agalati wind.
The emerging spetrum of a LB in a NEI state is haraterized by an unusual
largeontributionfromthe reombinationlines, insuhaway that thehigh energy
part of the spetrum is ompletely dominated by reombination photons. Due to
this highenergy reombinationtail,aNEI plasma witha temperature of
4 × 10 4
Kan mimi the spetrum of
∼ 10 6
K CIE plasma in shape and in hardness, seegure2afromBreitshwerdt and Shmutzler [92℄. Also,this delayed reombination
emission lines leave a strong singnature in the plasma spetrum with a saw-tooth
prolein the upper energy part of the spetrum.
In this type of model the existene of ion speies like O vi, N v and C iv
intheinterstellarmediuman beeasilyaommodatedwithinthe lowtemperature
plasmas. Beauseofthehighlyionizedspeiesandbeausetheydidnothaveenough
timetoreombine,they are abletoexistwithinalargerangeof temperatures. This
regimes.