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

книги из ГПНТБ / Григоришин, И. Л. Моделирование электроннооптических систем на сетках сопротивлений

.pdf
Скачиваний:
13
Добавлен:
19.10.2023
Размер:
4.89 Mб
Скачать

Дифференциальные уравнения движения электронов в электростатическом поле при скоростях, близких к ско­ рости света, учитывают изменение массы

 

d

am,.

dx

= e,

<5cp

 

d1

 

dt

 

dx

 

d_

am.

dy_

 

<5cp

 

dt

 

dt

 

 

где a =

 

c — скорость

света.

1

c~

 

 

 

 

Как показано в [10, 68], траектории релятивистских электронов с массой m = ain0 в поле потенциала ср совпа­ дают с траекториями иерелятивистских электронов с мас­ сой /«о в поле с эффективным потенциалом

Ф* = ф ( 1 ■1-

т0с-

ф ) =

2

1 Ф-

(2.25)

V 2

)

 

 

Действительно, если в узловых точках создать распреде­ ление эффективного потенциала (2.25), то траектория релятивистского электрона может быть рассчитана так же, как и перелятпвистского. При этом следует иметь в виду, что эффективный потенциал ср* не является гармо­ нической функцией и, следовательно, невозможно полу­ чить его распределение внутри области, задавая его зна­ чение лишь па границе. Помимо задания потенциала и* на границе для получения с помощью сетки сопротивле­ ний распределения эффективного потенциала, необходи­ мо в каждую узловую точку задать дополнительные токи для получения равенства (2.25). Но эта операция доволь­ но трудоемкая и в ней пет необходимости при расчете траекторий релятивистских электронов. Достаточно на сетке сопротивлений получить распределение истинного потенциала ср и перевести его в эффективные значения по

(2.25) [70].

Введем понятие эффективной напряженности поля

Е* = — grad ср*,

60

компоненты которой по аналогии с (2.1) будут:

ГГ*

<Эф

дф*

--•

ф*+1,т— Фй-1,ш

jC Y -■

СС

дх

дх

 

2h

 

 

 

 

 

(2.26)

 

дф

 

дф*

 

ф*

,. — ф*

 

 

 

 

11 ~

 

ду ~

ду ~

 

2h

 

 

1

/

 

 

 

движения релятиви­

Заменой т* = — у

 

t уравнения

стского электрона преобразуются к виду

сРх а дф

dx*“ 2 дх

(2.27)

dhy а дф

dx*2 2 ду

Система уравнений (2.27) по своему виду аналогична (2.2). Решение ее в элементарном квадрате ABCD (см. рис. 2.1) при условии, что компоненты эффективной на­ пряженности в этом квадрате постоянны и определяются по (2.26), записывается в виде, аналогичном (2.3), (2.5) с заменой__истинного потенциала ф на эффективный ф*,

причем 'j/ф* соответствует начальной скорости в поле эф­ фективного потенциала и связан с ней соотношением

Траектория релятивистского электрона может быть вычислена также путем решения уравнения траектории в общем виде (2.13) с использованием распределения по­ тенциала (2.6) в прямоугольнике ABCD (см. рис. 2.2) при замене в (2.6) и (2.13) ф на ф* и выражения (2.14). Одна­ ко, чтобы снизить погрешность в определении эффектив­ ного потенциала внутри прямоугольника, возникающую из-за допущения о линейном изменении ф* на его сторо­ нах, целесообразно это допущение оставить для истинно­ го потенциала, а значение эффективного потенциала в прямоугольнике вычислять по (2.25), учитывая (2.6), т. е.

61

Ф* (*. У) = (Xхх) (у — y j К + (X — Xj) М + (у—ух) F +ср4+

-г - 6° . [(^ — JCi) (У — Уг) К (х — хг) М + {у—Уг) /Ч ср4]2.

2т0с2

На рис. 2.4 приведем пример численного расчета тра­ ектории релятивистского электрона (кривая 2) в уско­ ряющем однородном поле; точки на ней соответствуют значениям, полученным по теоретической формуле [68].

Рис. 2.4. Траектории электронов в однородном ускоряющем поле: 1 — нерелятнвнстского, 2 — релятивистского при потенциале в начальной точке 100 кВ, £ = 4 0 кВ/А, угол входа у=75°

Для сравнения показана траектория нерелятнвнстского электрона (кривая /).

Используя метод эффективного потенциала, можно рассчитать траектории релятивистских электронов в скрещенном электрическом и магнитном полях. Вводя

обозначение х'*~~^~^> приведем систему уравнений дви-

атп

жения (2.15) для случая релятивистских электронов к виду

d2x

тп Е *

 

dy

j

I♦*

е0В2

Л'

 

dx'*

а т

 

 

 

d2y

тп

р

* __

dx

dx'*'

е0В*

j /♦

 

 

ат

Эта система аналогична (2.15), и решение ее записывает­ ся в виде, аналогичном (2.17) — (2.20).

62

Таким образом, процессы расчета траектории нереля­ тивистских и релятивистских электронов принципиально не отличаются, если известно распределение эффектив­ ного потенциала. Последний может быть рассчитан по значениям истинного потенциала, причем только в тех точках, которые необходимы для расчета траекторий.

§4. РАСЧЕТ ТРАЕКТОРИЙ ЭЛЕКТРОНОВ

БКОМБИНИРОВАННОМ ЭЛЕКТРИЧЕСКОМ И МАГНИТНОМ ПОЛЕ

СОСЕВОЙ СИММЕТРИЕЙ

Электроннооптнческие системы, в которых поля и потоки заряженных частиц симметричны относительно некоторой оси, получили наибольшее практическое при­ менение. С их помощью оказалось наиболее просто со­ здать электрические и магнитные линзы для фокусировки электронных пучков в электронных микроскопах, элек­ троннооптических преобразователях, электроннолуче­ вых трубках и т. д. Электромагнитные поля с осевой сим­ метрией широко применяются для формирования интен­ сивных высокоэиергетических электронных пучков в СВЧ-приборах, в электроннолучевых технологических установках для сварки, фрезерования, плавления и испа­ рения металлов и т. п. Для компенсации расталкивающих сил кулоновского взаимодействия между электронами потока с большой плотностью заряда часто применяется осесимметричное магнитное поле. Возникающее при этом азимутальное движение электронов, обусловленное как радиальной В г, так и осевой Bz компонентами магнитной индукции, при их определенном соотношении компенси­ рует расталкивающее действие пространственного заряда и центробежную силу.

Поведение электрона в комбинированном поле, сим­ метричном относительно оси г, описывается системой уравнений [21]

&_ At

( 2.28)

63

Трудность решения такой системы в элементарных функ­ циях очевидна. Обычно для этой цели используются мето­ ды экспериментального определения траекторий на элек­ тролитической ванне [9] или численные методы.

При дальнейшем рассмотрении этого вопроса будем исходить из того, что распределение потенциала электри­ ческого осесимметричного поля найдено с помощью осе­ симметричной сетки сопротивлений в дискретных точках на плоскости г, z, а созданное системой внешних источни­ ков магнитное поле по крайней мере в этих точках задано [29]. Рассмотрим общий релятивистский случай движения электрона в комбинированном поле. Вводя новую пере­ менную

и используя известные соотношения между компонентами вектора магнитной индукции и магнитным потоком ф [2 1 ]

Вг

 

1

 

дф

в г =

1

 

дф

 

2 яг

 

dz

-

 

dr ’

 

 

 

 

 

2 яг

 

реобразуем систему уравнений (2.28)

 

 

 

 

d2r

 

a

dip

 

 

dQ

V2

 

 

dx*2

~

2

dr

 

 

dx*

)

 

 

 

1

 

/ 2en

 

(Эф

dQ

 

 

 

 

}/

/п»

 

dr

dx*

 

 

d2z

a

dcp

1

1 /

2 е0

<Эф

dQ

dx*"

2

 

dz

к

 

 

dz

dx*

 

2

dQ

 

1

/

ф л . с

 

 

dx*

 

/

 

 

Г

 

4я |/

 

т0

 

 

 

Постоянная Сг для данной траектории определяется из на­

чальных условий т* = 0, г — гу, Ф =ФГ,

=

0! =

____

'

\ dx* / т»=о

 

= ац ^ “5 -

и рав,,а с * = '■’ * ? - ж

V %

гда

0* — приведенная угловая скорость электрона. Рассматривая

64

движение электрона в меридианной плоскости г, z, враща-

d9

ющеися относительно оси z с угловой скоростью ■— , и dx* ’

вводя эффективный потенциал (2.25), представим уравнения движения электрона в виде

dV

a

ф

 

V % - * + c >

dx*2

dr

 

' - M - b

drz

_a_

 

1

(2.29)

,.2

l / ¥ ^ c‘)

dx*2

dz

 

 

 

 

Данная система уравнений аналогична системе (2.2). Следовательно, она описывает движение электрона в не­ котором потенциальном поле, характеризуемом эквива­ лентным потенциалом

ф* = Ф*

(2.30)

Если распределение потенциала ф в дискретных точ­ ках области известно, то по (2.30) может быть вычислено распределение эквивалентного потенциала Ф* для каж­ дой траектории. В отличие от рассмотренных ранее слу­ чаев градиент эквивалентного потенциала в элементарной ячейке ABCD (рис. 2.5) может оказаться слишком боль­ шим, так как он растет с увеличением магнитного поля. Поэтому предположение о постоянстве компонент напря­ женности поля эквивалентного потенциала в квадрате ABCD может привести к грубым ошибкам. Чтобы избе­ жать их, поступим следующим образом. Разложим коор­ динату г траектории в ряд по степеням Дz, ограничиваясь членами второго порядка,

где Az — z zp. В

исходной точке

Р (рис.

2.5) значения

dr

d2r

Zp), ty{rp, zp), а

также

истинный по-

-----,

--------- , cp (rp,

dz

dz2

 

 

 

тенциал в узловых точках сетки предполагаются известны­ ми. Если А2 и соответствующее приращение координаты г малы, то можно считать, что на интервале Да величины

5. За к. 596

65

dr

dh

постоянны. Следовательно, координату гРг

dz

dz1

точки Р,

можно определить по (2.31).

Для

расчета коор-

динаты гр найдем величину

(

dr

\

 

 

 

 

\

-----

jpt

 

 

 

 

 

 

 

 

dz

 

 

 

 

dr

 

dr_

 

 

d1r

 

(zp,

zp)

(2.32)

 

dz p,

dz ; p

 

~~df

 

 

 

 

 

/ p

 

 

 

 

 

I

r'\

U

 

!v,

 

I

 

 

 

 

I

 

 

 

 

 

 

 

 

I

 

i

 

 

I

 

 

 

 

I

A

 

 

 

I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

\ь,п

 

jlа ,«‘/

\Ук,п*г

 

 

 

 

 

г,1 —•» " га " t “

"

 

 

 

 

 

I

р'Т>. Ъ рг.

 

I

 

 

- +

 

 

 

 

■— ч

L-------------

 

 

 

 

|

 

I

c

'

 

t,nrl

 

 

 

 

 

 

[fh-l.n

!'A

 

 

 

 

 

l

 

I

 

I

 

 

 

 

 

 

 

I

 

I

 

I

 

 

 

 

 

a ~

Рис. 2.5. К выводу формул расчета траектории

релятивистского электрона в комбинированном

 

 

 

 

осесимметричном поле

и величину

(

^ Г

\

, для чего воспользуемся уравнением

\

-------

 

 

dz-

/я,

траектории в общем виде

ч‘г _

,+ Ы )

 

 

(2МД33) Л

dz2

2Ф*

{

dr

dz dz I

где Ф* определяется по (2.30).

Рассмотрим, каким обра-

'зом можно определить величину эквивалентного потен­ циала Ф* в произвольной точке элементарного квадрата ABCD. Предположим, что компоненты градиента потен-

G6

циала ср в данном квадрате постоянны и определяются по (2.1). Тогда для потенциала в произвольной точке квад­ рата можно записать

ср(Л г)

известен

С р = ф { Г р ,

Подставив

Ф/r+i.n Фа-1,

Фа,п+1

Фа,п—1 z -j-C0, (2.34)

Ш

 

 

 

/’г>’ 2d ~-координаты

точки

 

константа С0 определяется

 

Vр> z p )

потенциал

ср (гр, <

Z p)~

 

1-1,71

Фй,П+1

r p

2h

 

 

(2.35)

в (2.34),

получим

 

,

, Фл+1,П

Фа—1,71 ,

.

Ф (г, z ) = ---------------—{Г—Гр) X

Ф/1,П-1

Zp

(2.35)

,

Фа,7i+i

Фа,тг-1

.

. ,

,

,

/Г1

Ч------

2-----

2 /Г ^ ------

(z — Zp) +

ф {Гр'

Zp)■

(2'36)

С учетом (2.25), (2.30) и (2.36) запишем выражение для эквивалентного потенциала в произвольной точке квадрата

Ф* (г,

г) =

Фа +1,71

Фа-1,77 {г — гр)

 

 

 

Фа,?1+1

Фа, 71-1

(z

Z p ) * ; (р (Гру Z p) X

 

2!г

 

 

 

Х 11+

ео

Ф&+1»71

Ф/1-1,71 {Г — Гр) +

 

 

 

2Ь

Фа,77+1

Фа,77-i

(г — Zp) + ф (Гр,

 

 

 

 

 

1

 

2е„

 

 

/

mn

Ф (/•. z) + Сг

67

Из этого выражения определяются компоненты градиента

эквивалентного

потенциала

 

 

 

 

 

 

 

ЭФ*

Фа+1,п

Ф;1- 1,77

|, ,

ео

Фа. 1,71

 

Ф/|-1,71

v

дг

 

2h

 

 

1^ т0с2

 

 

2h

 

 

 

 

 

 

,

Ф/ ,77+1

Ф/7,71-1

 

 

 

 

 

 

х

(/- -

" 1 Р) ' Г ~

*

2к

 

-(Z -~ Zp) 'I '

 

 

 

 

1)

2

Г

1

/

2еп

 

 

 

 

 

 

ф('>> z p )

Л ■7 г _ 4л |/

 

то- Ф ,

z) -j-Cx X

 

 

 

 

1/

~еп

 

 

д\|>

 

 

 

 

 

 

\ 4л

[ф ('">

z) — r а г

г С ,},

 

 

 

то

 

 

 

,ф*

 

Ф/|,71+1

Ф/1,,1-1

| .

 

 

 

 

 

 

 

dz

 

 

2h

 

\ '

тпс2

 

 

 

X

Ф л + 1,77

Ф/<— 1,71

(г — гр)

Ф/<,7! I-1

‘ Ф/1 ,71 - 1

(z — Zp)

 

 

 

2/г

 

 

 

р> 1

 

2/г

 

 

 

 

4" Ф

{ Г р > Z p )

 

 

 

2ё7 ^

Г 1

 

 

 

 

 

2лг2 |/

т0

dz

|4я

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя в (2.33) значения

dr

Ф*,

дФ*

дФ*

dz

dr

dz

 

 

 

 

 

 

 

 

 

найденные в точке

Р1 (г = rp ,

z — zp )

при

известных

\|),

<Эгр>

dil)

 

 

 

'

 

/

d1/-

\

по (2.31)

 

—— , — в той же точке, получаем ------- и

 

dr

dz

 

 

 

 

 

 

\

dz1

j Pi

 

 

определяем координату траектории в точке Р2.

Расчет

ко­

ординат точек

Р3,

Р4

и

т. д.

выполняется

аналогичным

образом. Для расчета траектории в следующей ячейке

предварительно по (2.36)

определяется

потенциал ср(л, г) в

точке Q, и процесс расчета

координат траектории при но­

вых исходных данных повторяется.

 

В заключение кратко рассмотрим расчет траекторий не­

релятивистских электронов,

движущихся в однородном маг­

нитном поле с нулевыми начальными

скоростями. В этом

случае 0О= О, С1 = —

1

у r w -фк,

я|) = лг2В и система

уравнений движения (2.29)

принимает вид

68.

 

dV

 

 

a

 

gpfi2

( r - r ^ f

 

dx2

 

 

5r

 

 

 

 

Ф — 8m0r2

 

d2z

 

 

_5_

 

e B°-

 

 

dx2

 

 

52

Cp------ (r2 _ r2)2

 

 

 

 

8m/~

где т

/

2ё^

Эквивалентный потенциал в любой точке

1/'

=

/»п

 

 

 

 

 

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

 

Ф (г,

г)

= Ф/.+1.П

Фй-1,п

(г — гр) +

 

 

 

 

 

2h

 

 

 

Фл,Т!+1

Ф*,П-1

(2 — 2p) +

q>(/>, 2 )

 

 

 

2h

 

 

 

 

 

 

 

gpfi2

 

(2.37)

 

 

 

 

8m0r2

 

 

 

 

 

 

 

откуда

находятся производные

 

 

 

 

Фй+1.П— фй-1,п

,

{ r l - r %

 

дг

 

 

2/i

 

 

 

 

 

 

4tu0r3

 

 

фй,/г-{-1

фк,п—'

 

(2.38)

 

 

 

 

 

 

 

2/i

 

 

 

необходимые для расчета траектории.

Подчеркнем еще раз, что при расчете траекторий изло­ женным способом интервал Az необходимо выбирать ма­ лым. Если в процессе расчета окажется, что при некото­ ром отклонении электрона в радиальном направлении эквивалентный потенциал получается отрицательным, то это свидетельствует о большой величине интервала Az.

§5. РАСЧЕТ ТРАЕКТОРИЙ ЭЛЕКТРОНОВ

ВНЕСТАЦИОНАРНОМ ЭЛЕКТРИЧЕСКОМ ПОЛЕ

Задача расчета траекторий электронов в переменных во времени полях связана с рядом затруднений. Извест­ ные способы ее решения основаны на использовании ме­ тода электролитической ванны с автоматическим траектографом, представляющим собой динамическую модель

(59

Соседние файлы в папке книги из ГПНТБ