Kompaneets
Kompaneets
Kompaneets
A.S.KOMPANEETS
Institute of Chemical Physics, Academy of Sciences, USSR
(Submitted to JETP editor October 15, 1955)
J. Exptl. Theoret. Phys. (U .S.S.H.) 31, 876-885 (November, 1956)
The role of the Compton effect in the establishment of equilibrium between quanta and
electrons is considered in the nonrelativistic approximation.
730
THERMAL EQUILIBRIUM 731
(~~)c =-~d-e ~[n(l +n')N(E) (l) Here p 2 /2 m = E, p and p' are the momenta of the
electron before and after the collision, n and n'
~~ t?e d.irecti,ons of propagation of the quantum.
- n' (1 + n) N (E + hw- hw')] dW. Ehrrnnatmg p from these equations, we obtain an
equation which determines w' as a function of
w, p and the angles of scattering. In this equation,
Here N is the distribution function for the free we set w = w '+ 1'1. and limit ourselves to terms
electrons, d 'T is the element of phase space of that are linear in /'1., at least while we determine
the electrons, dQ is the differential probability of the integral which contains /'1.2. After simple
transition from a given state into another, com-
transformations, we get
patible with the laws of conservation of energy
and momentum. The index C on dn/dt denotes h/'1. =_ hcoo (p, n- n') + (hoo)2 (1- nn')
that only Compton processes are considered in the mc 2 [1+ +
(hoo I mc 2 ) (1 nn')- pn I me} • (4 )
equation. Statistical equilibrium among the elec-
trons in the plasma is established very quickly, It is valid to replace the brackets in the denomina-
ind~p~mdently of the radiation; therefore, we can tor of Eq. (4) with unity for the case of interest to
consider the distribution function N (E) to be Max-
us. The first term in the numerator has the order
wellian. Then, if we replace n in Eq. (l) by the
Planck distribution n = ( ehw/ k T - l)- 1 , then the of magnitude (kT/mc 2 ) 3 1 2 and the second term the
order (kT/mc 2 ) 2 • But we shall determine the
right and the left parts vanish, as they must.
We shall consider that the energy of the elec- integral of 1'1. 2 ; therefore, the contribution of the
trons is nonrelativistic; i.e., we assume that the first term in the averaging over the angles does
inequality kT << mc2 holds. Then the energy not vanish. Consequently, upon substitution in
Eq. (2), we must compute the following integral:
transferred in each separate act is small in com-
parison with the energy quantum h w: w' - w 1'1.= I= (hw I mc) 2 ~d-e ~ dW (p, n- n') 2 N (s). (5)
<<w.
Making use of this inequality, we expand the
integrand of Eq. (l) in a power series in 1'1. up to After averaging over all directions of p we get
terms of the second order, inclusively. We set
hw/kT =x:
iJn)
( 7ft [ iJn
c= ax-+n(l
Jh \ \
+n) kT·jd-c ~dWN(s)~
(2)
The first integral is equal to 2m (3kT /2)
= 3 mkT. In the second integral, we replace the
Compton cross section in the nonrelativsitic ap-
proximation by the Thomson cross section, which
is symmetric relative to scattering over the angles
() amd (77- ()),so that J nn'dW = 0. The Thom-
son cross section does not depend on the energy
The second integral on the right side of Eq. (2) is of the quantum; consequently,
much easier to compute than the first ; but it suf-
fices to compute only one of the two integrals. The I= (hw) 2 (kT I mc 2 ) c I l, (7)
other is determined from the condition that the
equation ought to guarantee conservation of the where l is the Compton range, which is determined
total number of quanta in the scattering. by the total cross section ( 811/3 ){e 2/mc 2 ).
732 A. S. KOMPANEETS
co co co
In order that Eq. (2) guarantee the conservation
of the total number of quanta, we must set £~ ~ nx dx
3 = 4 ~ nx dx -3 ~ nx4 dx. (14)
0 0 0
(8)
~ d'" ~ dW' N (s) il While h w is still small in comparison with kT, we
can neglect the second integral on the right in
= (kT I mc 2 ) (c ll) hw (4- hw I kT). . . te 4 Y. This means that the time
(14). Here x = x lnl
for the energy to increase by a factor e is equal to
Equation (8) can he obtained in the following way. '"c = (mc 2 14kT) ll c. (15)
The kinetic equation for the quanta must look like
a sort of conservation law 3, BREMSSTRAHLUNG SPECTRUM
The bremsstrahlung spectrum in the general case
appears rather complicated. We assume that such
where j is the "flow" of quanta in frequency conditions are satisfied in which the Born approxima-
space. Since Eq. (2) is of second order relative tion is applicable for the determination of the
to x and contains the second derivative d2n/dx 2 spectrum. These conditions are satisfied in the
linearly, without a coefficient depending on n, the light elements for sufficiently high temperature
current must have the form of the sum of the first (from several tens of kilovolts and higher). At
derivative an; ax and a certain function depend- the same time, it is taken into account that the
ing on n. But in the state of total equilibrium, inequality kT << mc 2 is satisfied. We shall not
when n = ( e" - 1 )-l, the flow vanishes. In this take into account the scattering of electrons on
case, an/ax= -n( l + n ). It therefore follows electrons, with emission of quanta, since such a
that process makes a significant contribution only in
very light materials and in each case it does not
j=g(x)[anlax+n(I+n)], (10) exceed the bremsstrahlung emission on nuclei,
always remaining smaller than it.
wher~ the function g ( x) must be determined. Sub- In the Born approximation, and in the nonrela-
stituting (10) in Eq. (9), and comparing with (2), tivistic case, the bremsstrahlung cross section is
we find that g = -x2, and the unknown function (Ref. 1, p. 183, Eq. (18)]:
f.. = x ( 4 - x ), which is equivalent to (8). Equa-
dw (16)
tion (8) is in qualitative agreement with Wien'~
law: the energy of the quantum increases so long
as hw < 4kT.
We now introduce the dimensionless time Y by In order to go over to the total number of quanta of
the expression a given frequency that are emitted per second, we
need to multiply the cross section by LN ( i) vd T
t = (mc 2 1kT) (ll c) = y. (ll) ( L is the number of nuclei per unit volume, v is the
In these units, the kinetic equation is written in velocity of the electron), integrate over all states
the following form: of the electron and refer to the number of states of
a quantum of given frequency per unit volume. Thus
on )c ------x
(-ny 1 a ,4 (an
-- x iJx -ax +n+n 2)
. (12) we obtain the equation
2
rate of change of the quanta as a result of brems- not with respect to number ) is in the fashion of
strahlung and free-free absorption: Brownian motion. The limit of the frequency, for
which the Compton process takes up the quanta
(dn I dt)B = (1 I '=V [(I + n)- neh"' I kT]. ( 18) more rapidly than it succeeds in being absorbed,
is determined by the following estimate: Tc/TB rv l.
The derivative vanishes if we put the Planck dis- Substituting TC and TB in it, we get
tribution in place of n
Digressing from the Compton process, and assum- (23)
ign that the total change dn/ dt appears on the left
side of Eq. (18), we can integrate (18). If n (0)
= 0, then we get
The entire process considered (i.e., the estab-
(19)
lishment of equilibrium by the Compton effect) pre-
n = 1
ehw 1 hT -1 (1- exp {- -:...!~-
.. (ehwJkT- 1)}) . sents interest only in the case in which h w << kT
0 '
since only then is there a frequency interval,
w >> w 0 , in which the Compton effect plays an
It then follows that the relaxation time for essential role. For small x 0 , the function K 0 can
thermal equilibrium of the quanta, for a pure brems- be replaced by an expansion in x 0 (Ref. 2, p. 96 ):
strahlung mechanism, is equal to K 0 (xoJ2) = ln(4/yx 0 ), where lny= 0.577. Thus
"B = 't~ 1(ehw 1kT _ 1). (20) x 0 is determined by the equation
The total kinetic equation, which takes into ac- Ax;;-2 In (2.35 I X 0 ) ~ 1. (24)
count the Compton and bremsstrahlung processes,
is now written This equation has meaning only if A is a small
number.
(21) We now compare the transfer of energy from the
electrons to the quanta by pure bremsstrahlung and
= kT
mc2
_c_ _1_ _i_
l x 2 iJx
x (anox 1, n + n2)
4
by the Compton mechanism*.
The total energy of the quanta emitted by the
electrons per unit volume is
quanta for sufficiently low frequencies will always (dE I dt)c--:-- 3kT ~ "il w2 dw 1-r; 2 c3. (26)
undergo transition to statistical equilihirurn in w,
processes of emission and absorption. The quanta
of high frequencies will he taken up by the Comp- The principal contribution to the integral is pro-
ton process and thus increase their frequency, vided by the low frequencies. Therefore, we can re-
approaching the Wien distribution. It should not place the function K 0 under the integral by its
he thought that frequency of the quantum increases
monotonically; Eq. (12) is of second order, of the
diffusion type. The approximation of quanta to * This comparison was pointed out to us by L. D.
incomplete equilibrium ( with respect to energy hut Landau,
734 A. S. KOMPANEETS
If we do not make use of the approximate expres- For large x/.x 0 , the number of quanta is much less
sion for K 0 , then there will still stand a numerical than the equilibrium number. Adding the factor
component in front of the logarithm in Eq. (27). ( 1 - nx) under the integral (27), a factor which
That is, in place of ln 2 ( 4/ yx 0 ), the integral gives takes into account induced emission and absorption
ln 2 ( 4/ y x)- 0.27. But we can hardly consider of quanta, we get a converging expression which
this to be accurate, since the frequency w 0 is can be integrated from w = 0 :
found by an estimation method. 00
(.~£
The ratio of the energy transferred by the Comp- '3kT B \ 1 (33)
ton mechanism to the energy transferred to the
) -
dt c - ' 1t2 c 3h .) -;r- w2 dw
0
quanta upon emission is equal to
for tjJ '. cp satisfies the equation (we are considering the
This can be integrated approximately (by the homogeneous equation)
WKB method) if we consider the coefficient for tjJ'
(36)
to be a large number. With this accuracy,
a a2 ( x 2) _
Ill' - C ex p (l--x-
I -
2 VA vr I n--:;;
2.35} . (31)
~=x---xm+ --:-+2-- Xcp=A(x)cp.
(}y OX2 T 2 X
A
A
Here A ( x) is an operator on the right-hand side.
Some simplifications of the calculations arise be-
z = In x, x= e-z12f ~'- (z) e-(ll."-'l.lY. (40)
cause of the fact that 1 (x) is a Hermitian oper-
ator*. The function f ( z) satisfies an ordinary differ-
Let us assume that the initial photon distribution ential equatio/: of the same form as the Schrodinger
is given in the form of some function n 0 ( x ); then equation for the vibration levels of the diatomic
the distribution at any later moment is written molecules, if the potential function of the nuclei is
symbolically as taken in the form of a Morse potential:
(37)
(44)
~ f(-i!l--3/ 2Jf{ifL-a/2) 0
-ico
(45)
X. 1 12 (x) = "}12 (1 - xj2) e-x/2, x_ 312 = 2-I/2xe-x/2,
We then have x,
ico
The first integral [we call it F ( u)] is computed by
= i~ ~ tg TCS {r (_!s ~2~ s) H (x 0 s)
1 expansion of M in a power series, The numerical
-ioo
values of F ( u) are the following: