introduction to plasma spectroscopyindico.ictp.it › event › a11188 › session › 5 ›...

68
2370-1 School and Training Course on Dense Magnetized Plasma as a Source of Ionizing Radiations, their Diagnostics and Applications H.-J. Kunze 8 - 12 October 2012 Institute for Experimental Physics V Ruhr-University, 44780 Bochum Germany Introduction to Plasma Spectroscopy

Upload: others

Post on 31-May-2020

1 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

2370-1

School and Training Course on Dense Magnetized Plasma as a Source of Ionizing Radiations, their Diagnostics and Applications

H.-J. Kunze

8 - 12 October 2012

Institute for Experimental Physics V Ruhr-University, 44780

Bochum Germany

Introduction to Plasma Spectroscopy

Page 2: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Introduction to Plasma Spectroscopy

H.-J. Kunze

Institute for Experimental Physics V, Ruhr-University,44780 Bochum, Germany

School and Training Course on Dense Magnetized Plasma as a Source of Ionizing Radiations, their Diagnostics and Applications

ICTP, Trieste, October 8 -12, 2012

Page 3: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

λ (nm)

L λ (W

/cm

nm s

r)2

1012

107

Blackbody

k T B e = 170 eV

Typical Spectrum

Continuum radiationbremsstrahlung free-free transitionsrecombination radiation free-bound transitions

Line radiation bound-bound transitions

Plasmas emit electromagnetic radiation,which contains information on the plasmas

Page 4: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

bremsstrahlung free-free transitionsrecombination radiation free-bound transitions

line radiation bound-bound transitions

Continuum states of electrons

First ionization limit

E(p)

E(q)

E(g)Ground state

E( )∞

Second ionization limit

Doubly excited states

Ekin

Energy levels of an atom or ion

Page 5: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Maximum of bremsstrahlung emissionλ

×max 620nm eV

ekT

kTe = 100 eV λmax = 6.2 nm soft x-ray region

Line radiation from ions present in even in pure hot hydrogen or deuterium plasmas(from walls or deliberately added for diagnostics or for radiation cooling, sometimes done in fusion devices)

50% of light or medium heavy elements with nuclear charge6 ≤ Zn ≤ 26 are completely ionized for

3.40.27eV

en

kT Z≥

Examples: Ar (Zn = 18 ) kTe 5 keVC (Zn = 6) kTe 120 eV

Page 6: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Local quantity which characterizes emission

Spectral emission coefficient , ,λ ν ωε ε ε

defined by3

3d W( , ) , i.e. unit

d d d m sr nmr

Vλε λλ

ΦΩ

is the flux out of the volume element dVΦ

Observation:Emission from the surfaceis a flux

Spectral radiance Lλ

( Intensity ! )

Plasma

dAD

DdΩ

θNormal

Page 7: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

2dd cos d

LA

λλ θ

ΦΩ unit 2

Wm sr nm

How does the flux, respectively the radiance, change through the plasma ?

d ( , ) (x, ) d ( , ) ( , ) dL x x x L x xλ λ λλ ε λ κ λ λ−

L x( ) L x x( +d )

emission absorption

Optical depth d ( , ) dx xτ κ λ−

d ( , ) ( , )( , ) ( , ) ( , )d ( , )

L x xL x L x S xx

λ λλ λ λ

λ ε λλ λ λτ κ λ

− −

0( , ) ( ', ) d '

xx x xτ λ κ λ−

Equation of radiative transfer source function( , )S xλ λ

Page 8: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Low density plasmasHot plasmas absorption negligible

Spectral radiance at the surface

00

( , )( , ) ( , ) dL S

τ λλ λλ τ λ τ

−0x

λ (nm)

L λ (W

/cm

nm s

r)2

1012

107

Blackbody

k T B e = 170 eV

0( , ) dx xλε λ

Plasmas become optically thick first at long wavelength !Lines ( ) ( )

lower state density oscillator strengthn q fτ λ

×resonance lines are most prone to absorptionbecause of high n(q) !

Page 9: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

00( , ) ( , ) dL x xλ λλ ε λ

Drawback of all spectroscopic measurementsinformation along the line of sight

Computer tomography

measurements of the radiance in a large number of directions

(practically not possible in plasma devices)

Radial symmetry

y

x

y

L y( )r

R

dx

Measurement of Lalong many chords

Abel inversion

Local values of ( , )rλε λ

Page 10: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Applicable also to devices with one hole for the measurement

Plasma Pinhole Detector plane

y’simple transformation

and to tokamaks with a limiter, where density contour linesare very close to a set of nested eccentric circles, the shift beingknown as Shafranov shift

and to plasmas with elliptic cross-section (Rev. Sci. Instrum. 63, 4757 (1992))

Page 11: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Distortion by a thick wall ?

Plasma

(a) (b) yy

L y1( ) L y2( )

Reflection from the wall ?(Wavelength dependent )

Can cause errors !!

Page 12: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Special case

Beam emission spectroscopy

Plasma

Particle beam

Direction ofobservation

Local information

Li-beamsHe- beams boundary diagnostics

Li-beams (energetic) magnetic field measurements

H-, D- beams charge exchange spectroscopy

Page 13: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Let’s now assume, we have the local spectral emission coefficient in the plasma. What does it tell us?λε

Continuum radiation

Recombination radiation at long wavelengths becomes negligible

Bremsstrahlung at long wavelengths for a pure hydrogen plasma 2 1 2

21/( )

( )f

B e

ef nk Tλε λ

λ

Weak temperature dependence, strong electron densitydependence electron density diagnostics

With impurities:2

2 1 21 ,

/( )( )

ie i zff i z

B e

n z nk Tλε λ

λ

All impurity ions of charge z contribute

Page 14: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Deviation of the emission coefficient from that of a pure H-plasma is a measure of the impurity contamination characterized in fusion plasmas by

( ) ( )( )ff ffeffz Hλ λε λ ε λ

At high densities plasmas become optically thick at long wavelengths emission corresponds to thatof a blackbody temperature

At short wavelengthsboth bremsstrahlung and recombination radiationexponential dependence in frequency

110

lg ( )ln B e

h constk Tν

νε ν

Slope of straight line gives eT

effz is some kind of a mean charge

Page 15: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

λ (nm)

L λ (W

/cm

nm s

r)2

1012

107

Blackbody

k T B e = 170 eV

Bremsstrahlung and Recombination radiation

B e

hk Te

ν−

M. Galanti and N. J. PeacockJ. Phys. B 8, 2427 (1975)

J. E. Rice, K. Molvig, H. I. Helava,Phys. Rev. A 25, 1645 (1982)

Laser-producedplasma

Tokamak

Examples

Page 16: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Line radiation

Emission coefficient of a linebetween upper level p and lower level q

4( ) ( )) (z

h A p n pqνε νπ

( ) is the transition probabilityA p q→

The quantity btained thus is the popul (ation density of the upper level)zn po

Theoretical considerations are needed to establish how ( )zn p

is determined by the parameters of the plasma

General approach:

One identifies the relevant processes and sets up the appropriate system of rate equations

Page 17: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

1

21

2( )

( ) ( ) ( )

( ) ( ) ( )

( ) ( ) ( )

(

z zz

z ze z

z ze z

cre z

radiative transitions A p A q h A p q

electron collisional transitions A q e A p e n X q p

electron collisional ionization A q e A g e n S p g

collisional recombination n g

ν

α

− −

− −

→ →

→ →

→1

11

1

( )

( )

)

( ) ( ) ( )

( ) ( *)( ) ( )

z z rre z

z z

z dre z

p

radiative recombination A g e A q h n g p

dielectronic recombination A g e A pA p h n g p

ν α

ν α

→ →

→ →

Proton collisions, photo-excitation and photo-ionization can be neglected in many cases although there are plasmas where they are important

Rate

Page 18: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

For the derivative of the ground state density d ( )

dzn gt

ionization from and recombination into levels of the lower ionization stage (z-1) should be added

, , and are rate coefficients, .

d( )

.

e

X S i e

v vv vf

ασ σ

This means, that in some cases high energy tails of the distribution function may markedly influence population densities and may thus be detected, otherwise they depend on Te

Page 19: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

So, if one were to know all transition probabilities and all rate coefficients as function of Teone should measure a spectrum as function of time, make a fit of the time dependent collisional radiative model to the observations, and one had the electron density ne and the electron temperature Te as function of time.

In general, this is still wishful thinking

However, in practice, simplifications are not only possible but even reasonable.

Page 20: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Looking at time scales allows first simplification in most cases:

Time scales of ionization and recombination are long compared torelaxation of excited level populations, which are at least given by1/A.

One solves separately the rate equation of ionization/recombination

1

e z z e e e z+1 z+1 e e

e z-1 e e e z z e e

d ( ) ( ) ( , ) ( ) ( , )d

( ) ( , ) ( ) ( , )

eff effz

eff effz

n g n n g S T n n n g T nt

n n g S T n n n g T n

α

α−

and0( )zdn p for p g

dt

Metastable levels ?

Excited levels are in quasi-equilibrium with their ground states,they follow instantaneously their ground state density

means

Page 21: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

In stationary plasma

Distribution of ionization stagescalculated for 2 densities:Increasing density shiftsIonization to lower stages

Simplification at low and high densities !

Maximum abundance of an ion gives rough estimate of temperature at low densities

Page 22: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

In transient plasmas: atoms go successively through theirionization stages till they reach equilibrium

1 1

1 1 1 1

e e

rr dr cr rr dr cre e z z

dd

( ) ( )

zz z z z

z z z z z z

n n S n n S nt

n n n nα α α α α α

− − −

In a rapidly ionizing plasma recombination is negligible

11

1e

dd

z zz z

z z

n nn S Sn t n

−− −

decaying part of the abundance 1 1/z zn n−

hence 1e

dd

zz

z

n n Sn t

i.e. it reflects directly the ionization rate to the next ionization stage

Page 23: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Steady state

Page 24: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Dielectronic recombination

High densities: electron collisions dominate and establish a population equilibrium between the levels, which corresponds to that of thermodynamic equilibrium. It is called local thermodynamic equilibrium (LTE)(no equilibrium with the radiation field !!). Explanation!

Dielectronic capture followed by a stabilizing transition (satellites !)

1 10×-15

1 10×-16

1 10×-17

Electron temperature (eV)k TB e

10 100 1000 10000

1 10×-18

1 10×-19

1 10×-20

Rat

ec o

effic

ient

(ms

)3

-1

αrr

(OVII OVI)•

αr(OVII OVI)•

S(OVII OVIII)•E (p )z *

Doubly excited states

δEkin

E ( ) = E (g)Z ∞ z+1

E (g)z

E ( )z p

Page 25: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Advantage: Equilibrium relations can be applied.

Boltzmann distribution:

( ) ( ) ( ) ( )exp( ) ( )

z z z z

z z B e

n p g p E p E qn q g q k T

−−

Saha-Eggert equation:1 2

1 13

1 2/

,( ) ( ) exp mit( ) ( )

z qz zB

z z B

e

e BB

En g g gn q g q k T m k T

n πλλ

∞−

1 13

12 ,( ) exp( )

z gz z

z z BB

e En U Tn U T k T

∞−

gz(p) statistical weight of level (p) of ion (z)Uz(T) partition function of ion (z)• B thermal deBroglie wavelength

Page 26: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

With decreasing electron density collisional coupling ceases to exist first between levels with the largest energy separation,which in many cases is between ground state and excited state

E(g)

E(q)

E( )•

Fast radiative decay

nth thermal level

above nth equilibrium relations still prevail

Partial local thermodynamicequilibrium (PLTE)

Page 27: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Low densities

Excitation of levels is typically from the ground state(maybe from a metastable state) by electron collisions,depopulation by radiative decay

Collisional processes are weak, radiative ones dominate

( ) ( ; ) ( ) ( )e ez zn Tn g X g p n p A p→ →

( ; )( ) , )( )( )

(ee

ez ez

n X g p T f n Tn p n gA p

→→

Excited state populations are low, most ions are in the ground state or in a metastable state

Low density limit is also called corona approximation

neX versus A

Page 28: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Distribution of the ions at low densities:

1 1 1

1

1 1

( ) ( )( ) ( )

( ) ( )

rr drz e z e z e z z

z z err dr

z z e z ee

n n S T n nn S T fn T T

T

α α

α α

known as corona distribution

Page 29: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Examples of spectroscopic methods

employing the fundamental equation4

( ) ( )) (pqzn p

hp q A p q

νε

π→ →

4 ( )( )( )pq z

zp qn

hp

A p qπ

νε

→→

At high densities with LTE satisfied for the ion:

Particle densities nz

4

( ) ( ) ( )( ) exp

( )

( )

( ) exp() )(

z z zz e

z B e

pqz e

zp e

z

z Bq

p q

n p E p E qU Tg p k T

hUh A p q

n

Tg p k T

π νεν

→→

Te and partition function Uz(Te) must be known to obtain nz!

Page 30: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

In case the upper level is in PLTE, the Saha Eggert equation yields3 2

11 2

2 12

4 /,( )( ) e( ) p

()x

)(z pz e B e

zz B

z

p eq z e

Eg g mp qh A p

k Tn gg p k Tnq

π επν

∞→→

Now also the electron density ne must be known, in addition.

At low densities in the coronal regime

4 (

( )( ) ( )( )

()( ) (

); )e z e

zz z

z

zz

z

pq

e

zp

A pn n g n pn X g

qh n

pA pX g pp q T

qA

π εν

→ →

→→→

Forbidden transitions of the M1-type between fine structure levelsof the ground state of highly ionized atoms.

Page 31: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Electron and proton collisions are fast enough to establishequilibrium population distribution given simply by the statistical weigthsThus one line suffices to calculate nz,

neither ne nor Te must be known

Further advantage: lines of many ions are in the visible.

Examples are the green and red iron lines of FeXIV and FeXVobserved first in the solar corona, now used in many tokamaks

Charge exchange with injected fast beams of hydrogen is an indispensable tool to measure the • -particles produced in a burning plasma (local).

1( ) ( ) ( )cxz z bn p A p v v n nσ→

Charge exchange is selectively into levels of high principal quantum number and there into high l, decay via Δnp= -1 transitions at long wavelengths

Page 32: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Helium-like ions

0.394 0.395 0.396 0.397 0.398 0.397 0.400

Wavelength (nm) • •

Rad

ianc

e (a

.u)

L

10000

0

20000

30000

40000

x

1 2 n , n>2s l l’

w

y

qr

k

z

FitExp

Ionization limit of 1 states s nl

Ionization limit of

Ionization limit of

Lithium-like

Helium-like

1 2s s2

1s2

1 3s

1 statess nl2

1 2 statess p nl

f Ee kin( )

Ekin

1 2s p

1s S2 1

0

1s2s S 1

0

1s2p P1

1

1s2p P3

21s2p P

31

1s2p P3

0

1s2s S3

1

1s S2 1

01s S2 1

0

w yx z

E1E1

M2

M1

2E1

2E1

E1+M

1

w: resonance liney: intercombination linex,z: forbidden lines

Dielectronic and innershell excitedsatellites

all close in wavelength

Page 33: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

( ) ( )w e w e Hen X T n gεResonance line:

Dielectronic satellite ( ) ( )dr drsat e He e Hen T n gε α

Innershell excitation ( )in insat e e Lin X T nε

1( )drsat

weF Tε

ε

2( )insat

ew

Li

He

nFn

Tεε

The z-line is essentially populated by recombination of thehydrogenlike ion:

3( ) H

He

w e

z F T nn

εε

Such spectra of heavy helium-like ions are now much used in tokamak plasmas (up to Ni). They are in the x-ray region.Look at the nice fit of the previous example!! They certainly are observed also in high-density hot plasmas

Page 34: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Electron temperature

One example for low densities we just have seen.

Taking line ratios is convenient since the ground state density and the electron density drop

' '( ) ( ) ( ' ) ( ) ( )( ' ') ( ' ') ( ) ( ')

p qz z z z

z p ze

q z z

p q A p q A p X g p fp q A p q A p X g p

Tλε

ε λ→ → → →→ → → →

You need to know the excitation functions !

With increasing density, ratios become also dependent on the electron density and one has to use, in general, collisional-radiative models

Ratios become simple at high densities, when the upper levelpopulations are described by a Boltzmann distributionCheck for condition !

Page 35: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

' '( ) ( ) ( ) ( ) ( ')exp( ' ') ( ' ') ( ')

p qz z z z z

z pq z z B e

p q A p q g p E p E pp Tq A p q g p k

λεε λ

→ → −−→ →

Instead of several line pairs:

calculate the population densities nz(p) of the upper levels, and plot lg[nz(p)] versus the energy Ez(p)

Boltzmann plotstraight line, the slope gives Te

For high sensitivity energy separation of levels should be large !

Lines from successive ionization stages

For example: Saha-Eggert equation connects population ofupper level (p’) of ion (z) with population of ground state of ion (z+1), and Boltzmann relation connects that with the Population density nz+1(p) of level (p) in that ion

Page 36: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

A number of examples are found in the literature. Drawback: valid only in some high density range

Electron densityFrom line ratiosSelect one line with an upper level (p) decaying essentially onlyby radiation to the ground state, and another line where the upper level (p’) is also depopulated by collisions

e

( ) ( ) ( ) ( )( ) ( ') ( ') ( ' ) ( ' )

z e z z

z e z z z

n g n X g p n p A pn g n X g p n p A p X pn

→ →→ → →

1

' '( ) ( ) ( ' ) ( )( ' ') ( ) ( ' ') ( ')

( ' )( ' )

p qz z z z

z pq z z z

ze

z

p q A p q A p X g pp q A p A p q X g p

An X p

p

λεε λ

→ → → →≅→ → → →

→×→

Page 37: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

ExampleRatio of intercombination and resonance line in helium-like ions(often used in laser-produced plasmas)

Spectral lines are broadenednatural broadeningDoppler broadeningStark broadeningvan der Waals broadeningresonance broadening

1

( ) ( ) ( )

( ) d

p qλε λ ε λ

λ λ∞

−∞

Natural broadening is typically negligibly small

( )λ line profile, line shape function

Stark broadening (pressure broadening) is by charged plasma particles (electrons and ions). Their electric field leads to broadening of the levels.

Further advantage: lines are close-by in wavelength!

Page 38: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Two approximations are usually employed to calculate the broadening. Very difficult to calculate in general, so this was to benefit of the experimentalists who had to furnish experimental verifications

Impact approximation and quasistatic approximation

Impact approximation: electric field produced by the electronsvaries so rapidly at the position of the radiator, that the effect canbe treated as collision.

In cases, where the broadening effect on the lower state can be neglected (one-state version), the half-width of a line is

where X(p) stands for the sum of elastic and inelastic collisions.

)p(e2/1 Xn≅Δω

Page 39: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Collisions with neighboring fine structure levels usually being strongest leads to

1 2 0 2 0 5/ . .( )

ex

e

with xTnωΔ

(weak temperature dependence)

Lines are also shifted! In principle, it also can be used for density diagnostics, but it is less accurate!

Broadening data (including the corrections by the ion contributions)are found in the book by Griem (Spectral Broadening by Plasmas, Academic Press, 1974)A comprehensive bibliographic data bank is maintained at NISThttp://physics.nist.gov/cgi-bin/ASBib1/LineBroadBib.cgi

Profiles have a Lorentzian shape

Page 40: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Quasistatic approximation:

The field produced by the ions is considered static during the emission process. This microfield splits and shifts upper and lower levels by the Stark effect. Perturbation theory gives

( ) ,k mpq pqE C EωΔ

where m=1 and m=2 holds for linear and quadratic Stark effect and(k) designates the Stark component. First calculation of the microfieldwas by Holtsmark, new calculations are mostly for high densities. The line profiles mirror the microfield distribution, certainly slightly modified by electron collision contributions

Because the linear Stark effect is much stronger than the quadratic Stark effect, the quasistatic approximation is appropriatefor hydrogen and hydrogenic ions.

Theory is still being improved, and good measurements are needed for comparison

Page 41: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Wavelength (nm)•460 470 480

4000

8000

6000

Spe

ctra

l rad

ianc

e (O

MA

cou

nts) Paschen-α line of HeII

at 468.6 nm measuredfor ne= 3.5 x 1024 m-3

kBTe= 6.8 eV emitted from a well diagnosed pinchdevice

0 831201 2 2 74 10

./

-3.nm m

enλ −Δ ×

The most widely used line for density measurements isthe H -line of hydrogen at 486.13nm with a central dip

Page 42: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

With increasing principal quantum number np of upper levels lines become broader and finally merge Inglis Teller limit

H Balmer series He Lyman series

Allows estimate of the perturber density

29 12 4 5 1 5 7 5 ,max-3lg . . lg . lg . lgm

zpZn nz≈ − −

Page 43: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Lines with forbidden components (e.g. HeI lines)

Electric field increasing with density mixes wave-functions of upper levels forbidden lines increase

ΔΔ• f• aλa λf

Page 44: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

When the self-absorption of a line becomes important (the optical depth τ grows)

Flat maximumtemperature measurement

In the wings,radiation canescape !

In the case of cooler boundary self-reversal

Page 45: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Transition 2p 2sof CIV

120 000

100 2

KkWcm

L

Page 46: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Particle temperature Tz

Motion of an emitter results in Doppler shift of emitted lines

,pq x

pq

vc

ω ωω−

where vx is the velocity component in the direction x of emission, and c is the speed of light.

Line profile and velocity distribution function are related by

D( ) d ( ) dz x xf v vω ω

The line profile thus mirrors fz(vx), and if this is Maxwellian, the profile has a Gaussian shape half width gives Tz

51 2 7 715 10/ z

z

/ eV./u

B

pq

k Tm

λλ

−Δ ×

Page 47: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Magnetic field measurements

Magnetic fields split lines into their Zeeman componentswhich gives, in principle, the magnetic field, if Doppler broadening does not dominate

1 21 2

// /

split

pqB z z

Bk T m

λ λλ

ΔΔ

This suggests to use long wavelength transitions (large λpq)in heavy ions (large mz)

magnetic dipole transitions between fine structure levelof the ground state

22P 22S transition of injected high-energy Li-beams at 670.8 nm

Page 48: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Motional Stark effect

Fast beams (e.g. injected hydrogen heating beams) experiencethe Lorentz electric field

in their rest frame. bE v B×

This gives a Stark-Zeeman pattern, where the Zeeman splitting can be smaller than the Stark splitting.

For diagnostics usually the Balmer-α line of hydrogen (deuterium)is employed. Excitation is primarily by proton collisionsfrom the ground state.

Use is made also of the polarization characteristics

Page 49: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Levinton et al. Phys. Rev. Lett. 63, 2060 (1989)

Early example

Page 50: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Some possible experimental mistakes, some of these are even found in publications

No accurate wavelength calibration of the spectrographic system,people used only data from the manufacturer

No absolute spectral sensitivity calibration (authors simply used reflectivity data of grating and mirror, sensitivity of detector from manufacturer, area of entrance slit, solid angle),

Required is in most cases the calibration of the completespectrographic system since alignment, imaging, solid angle, apertures, etc. influence the total sensitivity

A number of sources emitting well-known emission lines are available

From the infrared into vacuum-ultraviolet spectral region variousprimary and secondary standard are available

Page 51: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

When studying two lines close in wavelength, no sensitivity calibration may be necessary

Difficulties with sources come up in the EUV and X-ray regions

Branching ratio methodemploys two lines from the same upper level

Calibration at long wavelength (e.g. in the visible) yields calibration at short wavelength providedBoth transition probabilities Aare known with sufficient accuracy

herep

q

q’

λ2

λ1

Page 52: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Transmission through air may have to be considered,in some cases also transmission through optical components, windows, lenses, filters

Wavelength (nm)•

Tran

smitt

ance

()%

T• (a)(c) (b) (b)(c)

0.1 1 10 1000.1

1

10

100

Transmission through 100 cm of gas at atmospheric pressureand at a temperature of 20°C: (a) He, (b) air, (c) Ar

In some spectral regions you have to work in vacuum !

Page 53: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Neglect of sensitivity change across a two-dimensionaldetector (pixel to pixel variation)

Irising effect in case of very fast gating a two-dimensional detectorsdue to finite travel time of gating pulse delay in switching from center to edge)

No check of the optical thickness of spectral lines; authors even claimed as consequence of such wrong interpretation that thetransition probability A changes at high electron densities. Several papers in top journal !

Saturation of detectors

Neglect of damping on fast signals in cables

Vignetting in the optical path

Page 54: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Signal from a photodiode recorded with a fastoscilloscope with a cable RG 58 C/Ua) 1 m long b) 11 m long c) 79 m long

Page 55: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Phys. Rev. A 19, 2260 (1978)

Wrong interpretation of Rayleigh scattering signals from a ruby laser beam with varying pulse duration t

Page 56: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Use of a system outside its sensitivity/design range !

Hydrino !!!

Page 57: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Some comments on spectrographic systems

Prism instruments are relatively easy to construct and to alignand have no overlapping orders

Exit plane

Entrance slit

1L 2L y’ ( ) •

1L

Entrance

Exitslit Littrow mounting

Page 58: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Instruments with gratings.

The most common mounting is the Czerny-Turner mount employing two mirrors;They are off-axis spherical and properly tilted and coma can be canceled. Astigmatism is rather small as well.

Entrance slit

Grating

Mirror

MirrorImage plane

Exit slit

A very modern variant of this mounting employs an echellegrating close to the Littrow mount and allows high resolution in high orders. An additional prism separates orders.

fEAE

Page 59: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

It is obvious that the radiation flux entering the instrument isproportional to AE E,where AE is the area of the entrance slit and E is the solid angle,AE E is called throughput.Since the width of the entrance slit E usually has to be small to obtaining reasonable spectral resolution (depending on the goal to measure total line intensities or line profiles), E is relevant forthe maximum achievable flux.Usual characterization of the radiation-collecting power by theso-called f-number f/D, where D is the dimension of the grating

With decreasing wavelength the reflectivity of the metal-coatedmirrors and gratings decreases,and the use of a single concave mirror becomes necessary, which combines both dispersion and focusing.The standard concave grating is of spherical shape,

has, of course, strong astigmatism

Page 60: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Grating of radius RRowland circle of radius R/2entrance slit and grating positioned on Rowland circlefocusing is on this circle for all wavelength.

Nowadays, the above and other aberrations can be reduced, by employing toroidal gratings, variation of groove spacingand groove curving

Page 61: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

The reflectivity of gratings is increased by going to small angles of incidencegrazing incidence mount

Wavelength range typically between 1.5 and 100 nm

Varying line spacing focal surface is flat !!

Such flat-field instruments are very compact

Page 62: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Spectroscopy at wavelength, say shorter than 1.5 nm, requires crystals as dispersing elements.

2 sinhkl

Bragg conditionm dλ θ

Radiation scattered coherently at the atoms of a crystal interfereconstructively only for one direction given by the famous

Selection of the crystal spacing dhklfixes the relation between • andBragg angle

Reflection from a plane crystalRadiation forms a cone for each wavelengthSpectral lines on a plane detectorare thus sections of a circle,an ellipse, a parabola or a hyperbola

Page 63: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Curved crystals focus like spherical gratings employing also theRowland mount

Here we see a mount used on large tokamaks. The plasma is outside the Rowland circle

Examples of other mounts are:

Page 64: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Detector surface

Source

Crystal

Axis

Johansson Cauchois mount

Van Hamos mountemploys a cylindricallybent crystal

Page 65: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

A very simple method the determination of the electron temperature provides the two foils method

L λ (W

/cm

nm s

r)2

1012 Blackbody

k T B e = 170 eV

One detects the radiation transmitted through two thin metallic foilsof the same metal but of different thickness or of different metals:

long wavelength radiation is absorbedshort wavelength radiation is transmitted and detected

The ratio of the radiationtransmitted throughproperly selected foilsis a strong function of the electron temperature Te.It is usually very convenient for plasma above 100 eV.One must be sure that there is no line radiation !

Transmission curves: http://www.-cxro.lbl.gov/optical_constantsfilter2.html

Page 66: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

How do you get started with the spectroscopic investigation of a plasma?

First you must estimate the temperature range you expect the plasma to have, since this determines the optimalspectral range of line emission and thus the type of instrument needed.

Each element has a specific spectrum for each ionization stagebut the spectra of ions along an isoelectronic sequence show similarities which change slowly with Z.

The spectra emitted by the working gas, which you know, are usually the ones you will use and analyze to obtain the plasma parameters.

In fusion plasma, where the working gas is deuterium, impurity species from the walls are utilized, in other plasmassmall amounts of impurities may be injected simply for diagnostic purposes.

Page 67: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

The selection of the detector in the exit plane of the instrument depends on

Do you want to record a spectrum at several times during the life time of the plasma?

You probably select a gated CCD camera and repeat the measurement at different times or for very short-lived plasmas you use a faststreak camera

or Do you want to record the time evolution of single lines?You employ photomultipliers or photodiodes behind an exit slit or even behind multiple slits, maybe by using optical fibers

Now you have to collect the relevant atomic data A, X, S, α

Although the available data collected in data centers are steadily growing, the shortage is still large.One continuously has to search the new literature

Page 68: Introduction to Plasma Spectroscopyindico.ictp.it › event › a11188 › session › 5 › contribution › 3 › material › 0 … · Introduction to Plasma Spectroscopy H.-J

Introduction to Plasma SpectroscopyH.-J. Kunze

Springer, Heidelberg, 2009

Collisional radiative models of varying complexity have been setup by several groups and are reported in the literature.

An extremely useful tool is provided by the FLYCHK code.It is maintained by National Institute of Science and Technology, NIST.It can be used interactively on the internethttp://physics.nist.gov/PhysRefData/contents.html

Lecture is based on