Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Шпоры (отличные) по МЭСА 7 семестр

.pdf
Скачиваний:
32
Добавлен:
10.05.2014
Размер:
946.18 Кб
Скачать

4.2

 

 

Модель Линхарда LSS

 

 

 

 

В основе – взаимодействие двух атомов, у которых распределение электронов по

 

 

расстоянию от ядра описывается статистической моделью атома Томаса-Ферми.

 

 

 

d 2Φ

 

 

 

 

(x)

=

 

Φ3/ 2

(x)

 

 

 

ТФ

 

 

 

 

 

ТФ

 

 

 

 

 

 

dx

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

его можно решить только численно

 

 

 

 

2

1/3

 

2

 

 

1

 

 

 

аТФ =

 

 

 

 

 

 

 

 

 

 

 

0,8853a0Z

1/3

 

 

 

 

 

 

 

2

 

 

 

1/3

 

 

 

128

 

 

 

 

 

Z

 

где x = r/аТ-Ф

 

 

 

 

mee

 

 

 

 

 

аппроксимация Мольер (баба)

ΦМ(x) = 0,35е-0,3х + 0,55е-1,2х + 0,1е-6х

так вот LSS, принимается

Φ(r / a)= ΦТФ (r / aL ),

aL = 0,8853a0 (Z12/3 + Z22/3 )1/ 2

В случае малых углов, угол рассеяния в СЦМ можно представить в виде

χ = −

dUdr

 

 

1

 

 

dr

μv2

 

 

 

 

 

r

2

2

 

ρ

ρ

 

 

 

dU

= −Z1Z2e2

 

1

 

r

1

Φ'

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Φ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dr

 

 

2

 

 

ar

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

a

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

замена переменной r = ρ/cosα,

 

 

 

dr = (ρtgα/cosα)dα

 

 

 

 

 

 

 

 

 

нижний предел переходит в α = 0, а верхний в α = π/2, тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z

Z

2

e2ρ π/ 2

cos2

α

 

 

 

 

 

ρ

 

 

cosα

 

 

 

ρ

 

 

dα

 

 

χ =

 

1

 

 

 

 

 

 

 

 

 

 

 

 

Φ

 

 

 

 

 

 

 

 

 

 

 

 

Φ'

 

 

 

 

 

 

 

=

μv2

/ 2

 

 

 

ρ2

 

 

 

 

 

 

 

 

 

 

 

aρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a cosα

 

 

 

 

 

 

a cosα cosα

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z1Z2e

2

 

 

 

π/ 2

 

 

 

 

ρ

 

 

 

 

 

 

 

 

 

 

ρ

 

 

 

 

ρ

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Φ

 

 

 

 

 

cosα

 

 

Φ'

 

dα.

 

 

 

 

 

=

 

 

 

2

/ 2

 

ρ

 

 

 

 

 

 

a

 

 

 

 

 

 

aμv

 

 

0

 

 

a cosα

 

 

 

 

 

 

 

 

 

a cosα

 

 

 

 

 

 

µv2/2 = m1v2/2(m2/(m1 + m2) = E0M2/(M1 + M2) и вводится

 

 

 

 

 

 

 

 

приведенная энергия Линдхарда:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε = E

 

 

 

 

M2

 

 

 

 

 

 

a

 

 

 

= 0,8853a

 

 

 

M2

 

 

 

 

 

1

 

 

 

 

 

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 M1 + M2

 

 

Z1Z2e2

 

 

 

 

 

 

 

0 M1 + M2

 

Z1Z2e2 Z12/3 + Z22/3

 

0

4.3

 

 

или, после подстановки значений

 

 

 

 

 

 

 

 

ε = 32,52

M2

 

 

1

 

 

 

 

E0 ,

где E0 в кэВ

 

 

 

 

 

 

 

 

 

 

 

 

 

M1 + M2

 

 

 

 

 

 

 

 

 

 

Z1Z2

Z12/3

 

 

 

 

 

 

 

 

 

 

+ Z22/3

 

 

Введя ξ = ρ/а , имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

π/ 2

 

ξ

 

 

 

ξ

1

g(ξ)

χ =

 

 

 

Φ

 

 

cosαξΦ'

 

 

dα =

 

εξ

 

 

 

 

εξ

 

0

cosα

 

 

cosα

 

εχ = g(ξ)/ξ

выражение получено в предположении малых углов χ, при которых χ/2 sin(χ/2), и

εχ = 2εχ/2 2εsin(χ/2) = 2t1/2 = g(ξ)/ξ,

в параметре t1/2 = εsin(χ/2) учитывается и массы с атомными номерами частиц, начальная энергия частицы m1, и угол рассеяния с прицельным параметром.

 

g(ξ)=

π/ 2

 

 

ξ

 

 

ξ

 

 

0

Φ

 

cos αξΦ'

 

dα

 

 

 

функцию

 

 

cos α

cos α

вычислить численно и для каждого ξ из равенства t1/2 = g(ξ)/2ξ определить t1/2 те получить функцию G: ξ = G1/2(t1/2) (Степень ½ у G выбрана для удобства).

далее т.к ρ = аG1/2(t1/2), получим

dσ = 2πρ

 

dρ

 

dχ =

2πρ

 

dρ

 

 

dt1/ 2

 

dχ

 

 

 

 

 

 

 

 

 

 

 

 

 

dχ

 

 

dt1/ 2

 

 

 

dχ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т.к. dρ/dt1/2 = aG'(t1/2)/2G1/2(t1/2),

 

где G'(t1/2) = dG(t1/2)/dt1/2

dt1/2/dχ = εcos(χ/2)/2, а dχ = dt/ε2sin(χ/2)cos(χ/2),

получим Линдхардовское сечение рассеяния

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dσL = πa

2

G'(t1/ 2 )t

dt = πa

2

 

f

L

(t1/ 2 )

dt

 

 

 

2t3/ 2

 

 

 

 

2t3/ 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где fL(t1/2) = tG'(t1/2) – функция Линдхарда.

fL находят по аппроксимации, Винтербона например

4.4

аппроксимация Винтербона

fL(t1/2) = 2,54(t1/2)1/2[1 + 2,164(t1/2)0.71]-2,1

Переход к дифференциальному сечению в СЦМ.

t = ε2sin2(χ/2), dt = ε2sin(χ/2)cos(χ/2)dχ.

 

 

 

 

 

f

 

εsin

 

χ

 

 

 

 

a2

 

 

 

 

 

 

 

 

2

dσ

 

=

 

 

L

 

 

 

 

2πsin χdχ.

L

 

sin

3 χ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в ЛСК

 

 

 

 

 

 

 

 

 

 

 

 

dσ(θ)

 

dσ(χ)

1+ γ2 cos 2θ

d

=

dω

2γcosθ±

 

 

 

 

 

 

 

 

 

1γ2 sin2

 

 

 

θ подставить получить

 

 

 

 

 

 

 

 

 

В диапазоне r/a =0,5÷5 функция Φ(r/a) может быть аппроксимирована прямой. если использовать дважды логарифмический масштаб, тоуравнение прямой

lgΦ(r/a) = a +blg(r/a) = =A + lg(r/a)b = lgA(r/a)b,

т.е. Φ(r/a) = A(r/a)b. А = 0,416; b = -1,и функция экранирования (функция Нильсен)

ΦН = 0,416а/r.

при 0,5а r 5а можно считать, что потенциал U(r) = (Z1Z2e2/r) (0,416а/r) = α2/r2 т.е. является обратноквадратичным потенциалом.

тогда подставляя α2в формулу сечения для обр.кв. потенциала из третьей шпоры получим

dσ =

aZ

Z

 

e2 (M

 

+ M

 

)

π2 0,416

 

πχ

dχ.

1

 

2

 

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

χ2

(2πχ)2

 

 

 

 

 

EM2

 

 

 

 

 

Первая дробь – обратная к энергии Линдхарда. Подставив dχ и учитывая, что

χ 2sin(χ/2) получим

 

 

 

 

 

 

 

 

 

dσ = πa2

 

 

 

1

 

 

 

 

0,416π dt.

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

χ

4

 

 

и сравнив это с последней рамочкой получим

 

 

2 ε sin

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

fL = 0,416π/4 = 0,327 для любых пар ионов

Ядерная и электронная тормозная способность, удельные потери, расчет для 5.1 кулоновского и экранированного кулоновского потенциала. расчет для

соединений.

Ион M1, Z1, движется в твердом теле(ТТ), состоящем из атомов M2, Z2. При взаимодействии иона с атомами происходит передача его энергии атомам, и ион тормозится в ТТ. В каждый рассматриваемый момент энергия иона равна Е, она уменьшается.

σ* – полное (проинтегрированное по всем возможным переданным энергиям) сечение процесса, сопровождающегося передачей энергии движущегося иона атомам ТТ. Если dσ*(E2) – дифф. сеч. передачи энергии в диапазоне E2 ÷ dE2, то

E2 max

σ* = dσ*(E2 ) где E2min – минимально возможная переданная ионом энергия,

E2 min

E2max – максимально возможная.

Выделим в ТТ объем с площадью основания σ* и высотой dl, направленной вдоль траектории движения иона. Концентрация атомов ТТ n*. В объеме находится n*σ* dl атомов (рассеивающих центров). Если E2 – средняя энергия, передаваемая рассеивающему центру в однократном процессе взаимодействия, то на участке траектории dl у иона изменится энергия dE2 = – n*σ*dl. Удельные потери энергии ионом вдоль траектории

dEdl = −n*σ* E2 .

 

 

 

 

E2 max

 

 

E2 max

 

 

 

 

 

 

E2dσ*(E2 )

 

E2dσ*(E2 )

 

 

 

 

 

 

E2 min

 

 

E2 min

 

 

E2 =

 

=

 

 

 

E2 max

 

σ*

 

 

 

 

 

 

dσ*(E2 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E2 min

 

 

 

S = E2dσ *(E2 )

 

тогда

dE = −n* E2dσ *(E2 ) и

-тормозная способность

 

 

 

 

 

E2 max

 

 

E2 max

 

 

 

 

 

dl

E2 min

 

 

E2 min

 

 

 

 

 

 

 

 

 

вещества

Упрощающие предположения:

ТТ представляет собой набор атомов с концентрацией n0, не обр. кристаллическую решетку.

Следуя Нильсу Бору, рассматривать потери энергии ионом как сумму двух независимых процессов, идущих одновременно:

взаимодействие с ядрами атомов, ("ядерное", подстрочный индекс "n"), в ходе этого процесса происходит изменение направления движения иона;

взаимодействие с электронами атомов ТТ, ("электронное", построчный индекс "е"), в ходе этого процесса направления движения иона не меняется, т.к. me << m1.

5.2

 

В этих предположениях

 

 

 

 

 

 

dE

dE

dE

= −(n*)n Sn (n*)e Se

 

 

 

 

dl

=

 

+

 

 

 

 

dl n

 

dl e

 

Для конкретного вида dσ*(E2) необходимо знать потенциал взаимодействия. Кулоновский потенциал.(Применим для быстрых легких ионов ~ МэВ). Дифф. сеч. по переданной энергии в случае кулоновского пот-ла имеет вид:

dσn *(E2 )=

πγ(Z1Z2e2 )2

 

dE2

 

 

 

 

 

 

E22

 

 

 

 

 

 

E

 

 

 

 

 

 

 

Ядерная тормозная способность при энергии иона равной Е

 

 

E2 max

 

2 2

 

 

 

 

 

2 2

 

 

 

Sn (E) =

πγ(Z1Z2e )

dE2

 

 

πγ(Z1Z2e )

E2max

=

ln

 

 

 

 

E2 min

E

 

E2

E

 

E2min

E2max = 4γЕ/(1 + γ)2 (вспомнить кинематический фактор Λ). Точное значение E2min не

существенно, ед. требование E2min 0. примем E2min = Ed = 10 эВ средняя энергия смещения атома из положения равновесия.тогда

Sn (E) =

πγ(Z1Z2e2 )2

 

 

 

E

 

 

 

 

 

 

 

ln

 

 

 

 

 

 

 

 

 

 

E

(1

+ γ)

2

 

 

и

 

 

 

 

 

Ed

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

 

 

 

E

 

 

 

 

 

 

 

 

 

 

 

 

dEn (E) = − πγn0 (Z1Z

2e ) ln

 

 

 

 

 

 

2

 

dl

 

E

 

 

 

 

 

(1+ γ)

Ed

Торможение на электронах можно рассматривать как рассеяние иона на электронах.

При таком рассмотрении потенциал взаимодействия U(r) = Z1e2/r. И, γ = m1/me >> 1.

E2max = 4γЕ/(1 + γ)2 4Е/γ = 4Еme/m1

Для вычисления E2min примем, что в результате взаимодействия иона с электроном переданная ему энергия больше энергии связи, т.е. рассеяние иона на атомном электроне сопровождается ионизацией атома.

Введем среднюю энергию ионизации атома (усредненную по всем оболочкам)

. При таком подходе E2min =

Se (E) =

π(m

/ m )Z 2e4

 

4(m

/ m )E

 

1

e

1

ln

e

1

 

 

 

E

 

 

Iˆ

, n* = n0Z2

 

 

 

 

 

 

5.3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dEe

(E) = −

2π(m1 / me )Z12Z2e4n0

ln

4(me / m1)E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dl

 

 

 

 

 

 

E

 

 

 

 

Iˆ

 

 

 

 

 

 

E = m1v2/2, где v – скорость иона и поэтому

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

4

 

 

2

 

 

 

2

 

 

4

 

 

 

 

 

2

 

 

Se =

Z1 e

 

 

2mev

 

 

dEe

(E) = −

Z1

Z2e

n0

 

2mev

 

 

 

 

mev

2

 

ln

Iˆ

 

 

,

dl

 

mev

2

 

ln

 

Iˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Электронная тормозная способность должна быть положительной величиной.

при каких энергиях иона выполняется это требование? Для ионов гелия (Z1 = 2, М1 = 4),

движущихся в олове (Z2 = 50) получаем условие E > /4(me/m1) 1 МэВ, т.е. условие применимости кулоновского потенциала.

вклад каких потерь является преобладающим?

dEdl n

dEdl e

 

 

 

 

 

4M1

 

 

 

E

 

 

Z2

 

ln

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

me M1

M2 (1

+ M1 / M2 ) Ed

 

 

 

 

4(m / m )E

Z2/М2 ½. , m1 = М1mа.е.м., и

 

M2 m1

 

 

 

 

 

ln

 

e

1

 

 

 

 

 

 

 

 

11,6Z2

 

 

mе/mа.е.м. 5 10-4. Возьмем М1 = 4, Е = 2 МэВ, Еd = 10 эВ. Тогда

dEdl n

dEdl e

 

 

 

 

2M1

 

 

E

 

1

 

ln

 

 

 

 

 

 

 

 

 

 

 

 

f (Z2 ) =

5,5 104

Z2 (1+ M1 / 2Z

2 ) Ed

 

 

 

4 5,5 104

E

 

2M1

 

 

 

 

 

ln

M1 11,6Z2

 

 

 

 

 

 

 

Практически для всех элементов, за исключением трансурановых, отношение

(dE/dl)n/(dE/dl)e < 10-2. т.е ядерными потерями можно попренебречь.

Экранированный кулоновский потенциал

Дифференциальное сечение по переданной энергии.( тот самый кинематический фактор Λ)

E2 = 4M1M2Ecos2Φ/(M1 + M2)2. Угол отдачи Φ = (π χ)/2.

Энергия Е через приведенную энергию Линдхарда

Е = εZ1Z2e2(M1 + M2)/aM2.

Т.к. Линдхарадовская переменная t = ε2sin2(χ/2), то E2 = 4M1Z1Z2e2 t

(M1 + M2 )a ε

прямо пропорциональна t, остальные члены - параметры,. Поэтому Линдхардовское сечение является также сечением по переданной энергии.

dE2 = dt4M1Z1Z2e2/(M1 + M2)aε

 

при t = ε2 переданная энергия E2 = 4M1M2E/(M1 + M2)2 = E2max,

5.4

при t1/2 = 0

E2 = 0 = E2min (для экранированного кулоновского потенциала

 

ограничения для E2min нет).

 

Ядерная тормозная способность в экранированном кулоновском потенциале

 

 

4Z

Z

e2M

1

 

 

ε

 

 

 

f

L

(t1/ 2 )tdt

 

aZ

Z

e2M

1

ε

f

L

(t1/ 2 )d(t1/ 2 )

Sn =

 

1

2

 

 

 

 

πa2

 

 

 

 

 

 

 

 

=

 

1

2

 

 

 

 

 

a(M1 + M2 )

 

 

 

2t

3/ 2

ε

 

M1 + M2

 

 

 

 

 

ε

 

 

0

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

sn (ε) = ε

fL (x)d(x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

ε

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Sn (E) = 8,47 10

15

 

 

 

Z1Z2

 

 

 

 

 

 

M1

 

sn (ε)

 

эВ см2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

атом

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z12/3

+ Z22/3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M1 + M 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dE

 

 

 

 

 

 

 

 

 

 

2

 

 

 

M1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= −an0Z1Z2e

 

 

 

 

 

 

 

 

 

 

 

sn (ε)

 

 

 

 

 

 

 

 

 

 

 

 

 

M1 + M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dl

n

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аппроксимация Юдина sn (ε) = 0,45 ε

0,3 +ε

такие же штуки для электронной тормозной способности.

Se (E) = 8,47

10

15

Z1Z2

 

 

 

 

 

M1

 

se (ε),

 

эВ см2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

атом

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z12/3 + Z22/3 M1 + M 2

 

 

 

 

 

 

 

 

 

0,0793

 

(M1 + M2 )3/ 2

 

 

 

 

 

 

 

 

1/ 6

Z1Z2

 

 

 

 

 

 

 

 

 

 

 

se = ke ε ,

ke = Z1

 

 

 

 

 

 

 

 

 

2/3

 

2/3

3/ 4

3/ 2

1/ 2

 

 

 

 

 

 

 

 

 

 

 

(Z1

+ Z2

)

M1

M2

 

 

 

образец – однородное соединение из N элементов и относительная концентрация каждого элемента Ci = ni /n0 то используются два подхода для вычисления S и dE/dl.

Средний атомный номер

 

 

N

 

N

 

 

2 = Ci Z2i ,

 

2 = Ci M 2i

Z

M

 

 

i

 

i

Правило Брэгга Если Si – тормозная способность i-го элемента, вычисленная по общим правилам, то для многоэлементного образца тормозная способность

N

S = Ci Si

i

6.1

Расчет траекторных пробегов ионов в ТТ и распределение внедренных ионов

по глубине. .

 

пробег -путь, который проходит ион в твердом теле до полной остановки.

До входа ионы имеют энергию Е0 (моноэнерг. пучок). Т.к энергия, теряемая ионом в каждом соударении с атомами ТТ (обозначим ее Т) зависит от прицельного параметра ρ, то эта

величина имеет случайный характер из-за случайности ρ. И угол рассеяния в каждом соударении случайный. →Траектория каждого иона индивидуальна и пробег R различен.

введем функцию распределения ионов по длинам пробега P(R, E0, Z1, Z2, M1, M2, n0) характеризует плотность вероятности остановится после прохождения пути R иону M1, Z1 с начальной энергией Е0 в образце Z2, M2 с концентрацией n0 (пренебрегаем наличием кристаллической решетки и рассматриваем моноэлементный образец).

предположение что распределение ионов по длинам пробега является Гауссовым

 

 

1

 

 

 

(R

 

)2

 

P(R, E0 ) =

 

 

 

R

 

 

 

 

 

 

exp

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2R2

,

R2

 

 

 

 

 

 

 

средний траекторный пробег R = R(E0 , Z1, M1, Z2 , M2 , n0 ) R(E0 )

R2 -среднеквадратичное отклонение

При взаимодействии иона с атомом ТТ есть вероятность, что атому будет передана энергия Tn за счет рассеяния на ядре, а его электронам энергия Te.

Она определяется значением дифф. сеч. dσne(Tn, Te) = dσn(Tn) + dσe(Te) (ядерные и электронные потери независимы)

Рассмотрим на входе иона в ТТ участок его траектории δR такой малый, что на нем происходит только одно взаимодействие с атомом ТТ.

Вероятность что на δR ион потеряет энергию Tn + Te, равна n0 dσne(Tn, Te)δR. энергия иона станет Е0 Tn Te. Чтобы ион имел пробег R, ему надо пройти R δR.

Плотность вероятности прохождения такого пути равна P(R δR, Е0 Tn Te). Произв. P(R δR, Е0 Tn Te) n0 dσne(Tn, Te)δR – вклад рассм. взаим. в полную вероятность.

Для учета всех возможностей передачи энергии проинтегр. по dσne. → вероятность что в

слое δR произойдет взаимодействие

P+ = n0δRP(R δR, E0 Tn Te )dσne

что не произойдет

P= (1n0δRdσne )P(R δR, E0 )

Ну и вероятность иону с нач. энергией Е0 пройти R: Р = Р+ + Р, преобразовывая

P(R, E0 ) = P(R δR, E0 ) +

+n0δR[P(R δR, E0 Tn Te )dσne P(R δR, E0 )dσne ]

P(R, E0 ) P(R δR, E0 ) = n0 [P(R δR, E0 Tn Te )dσne P(R δR, E0 )dσne ]

δR

 

6.2

 

 

 

 

P(R, E0 ) = n0 [[P(R, E0 Tn

Te )dσne

P(R, E0 )]dσne ]

 

 

 

 

 

 

 

 

 

При δR 0

получим уравнение для функции распределения Р

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P(R, E0) определяют с помощью ее моментов распределения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= P(R, E0 )RndR =

 

 

 

 

, a1 =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

an

Rn (E0 )

 

начальный момент n-го порядка

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= P(R, E0 )(R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

μn

 

 

)n dR = (R

 

)n

 

 

 

 

Центральный момент n-го порядка

 

 

 

 

 

 

R

R

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P(R, E0 )dR =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

μ2 = ∆R2 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Умножив обе части уравнения на Rm и проинтег. по R от 0 до

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Rm P(R, E0 )dR = n0 RmdR[P(R, E0 Tn Te ) P(R, E0 )]dσne

 

 

 

 

 

 

 

 

 

0

 

 

 

R

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Rm P(R, E0 )dR = Rm P(R, E0 )

 

0 mRm1P(R, E0 )dR = −mRm1

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n0 dσne Rm P(R, E0 Tn Te )dR Rm P(R, E0 )dR

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

рекуррентное соотношение

 

= −n0 dσne [Rm (E0 ) Rm (E0 Tn Te )]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= n0 dσne [Rm (E0 )

 

 

 

 

 

 

 

]

 

 

 

 

 

 

 

 

 

mRm1(E0 )

Rm (E0 Tn Te )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 = n0 dσne [R(E0 )

 

 

 

 

 

 

 

]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m = 1

 

 

 

 

 

R(E0 Tn Te )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Разность в ряд относительно Е0

по порядку Tn + Te и в первом приближении получаем

 

 

 

 

 

 

 

d

 

 

 

 

(T +T )dσ = n

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 = n

 

 

R(E0 )

R(E0 )

 

(T +T )dσ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dE

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

dE

 

n

e

 

ne

 

 

 

0

 

 

 

 

 

 

n

 

 

 

e

ne

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

=

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

т.к.

(Tn +Te )dσne

= S то

R(E0 )

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dE

 

 

 

 

 

 

n0S(E0 )

 

n0[Sn (E0 ) + Se (E0 )]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

E

 

 

dE

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R(E0 ) =

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

окончательный вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n0 0 Sn (E0 ) + Se (E0 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6.3

 

 

 

 

=

 

 

 

 

2 - средний квадрат траекторного пробега

 

 

R2

R2

 

 

R

 

Из рек. соотношения для m = 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= n0 dσne[

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

R2 (E0 )

 

R2 (E0 Tn Te )

]

 

 

 

 

R(E0 )

вычтем для m = 1, умноженное на 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R(E0 )

 

 

 

 

 

 

 

 

2

 

 

= n0 dσne[2

 

 

2 2

 

 

 

 

 

 

 

 

 

] получим

R(E0 )

R(E0 )

R(E0 )

R(E0 Tn Te )

0 = dσne[

 

 

2

 

2 +2

 

 

 

]

R2 (E0 )

R2 (E0 Tn Te )

R(E0 )

R(E0 )

R(E0 Tn Te )

Прибавим к обеим частям dσne

 

 

 

2

 

 

 

 

 

 

R(E0 Tn Te )

, перегруппируем

{R(E0 )2 R(E0 )2 [R(E0 Tn Te )2 R(E0 Tn Te )2 ]}dσne = = [R(E0 )2 2R(E0 ) R(E0 Tn Te ) + R(E0 Tn Te )2 ]dσne .

[R2 (E0 ) −∆R2 (E0 Tn Te )]dσne = [R(E0 ) R(E0 Tn Te )]2 dσne.

разложим выражения в скобках в правой и левой части в ряд по Tn + Te относительно Е0

 

d(

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

2

 

 

 

 

 

 

R2 (E

))

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R(E

)

 

 

 

 

 

 

 

dE

0

 

 

 

 

(Tn +Te ) dσne =

 

0

 

 

 

(Tn +Te )2 dσne

 

 

 

 

 

 

 

 

 

 

 

 

 

dE

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

введя (E) = (Tn +Te )2 dσne и помня (Tn +Te )dσne

= S(E) получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

d(R2 (E

))

 

 

 

dR(E

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

S(E) =

0

 

 

 

 

(E)

 

 

 

d

 

 

/ dE =1/ n S(E)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т.к

R(E

)

, то

 

 

 

 

dR

 

 

 

 

 

 

 

 

 

dE

 

 

 

 

 

 

 

 

 

 

 

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d(R2 (E

0

)) / dE = Ω(E) / n2S3 (E)

и в итоге

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

E0

 

 

(E)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R2 (E0 ) =

 

 

 

 

 

 

 

dE

 

 

 

 

 

 

n2

[S

n

(E) + S

e

(E)]3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Экранированный кулоновский потенциал

 

 

 

 

 

 

 

 

 

 

вспомним dE/dl = –4πan0Z1Z2e2M1s(ε)/(M1 + M2),

 

 

где s(ε) = sn(ε) + se(ε),

 

dE = dεZ1Z2e2(M1 + M2)/aM2

тогда средний пробег:

 

 

 

 

 

 

= 0

dE

ε

 

= 0

R

(dE / dl)

 

E0

0

 

 

Z1Z2e2 (M1 + M2 ) / aM2

 

 

(M1 + M2 )2

 

ε0

dε

 

 

 

 

 

 

dε =

 

 

 

 

an Z

Z

e2M

s(ε) /(M

+ M

)

a2n M

M

2

s(ε)

0

1

2

1

1

2

 

 

0

1

 

0