Научная статья на тему 'Математические модели, основанные на релятивистском уравнении Гамильтона-Якоби в лазерной физике'

Математические модели, основанные на релятивистском уравнении Гамильтона-Якоби в лазерной физике Текст научной статьи по специальности «Физика»

CC BY
423
65
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ВОЛНОВОЕ УРАВНЕНИЕ ДИРАКА / РЕЛЯТИВИСТСКОЕ УРАВНЕНИЕ ГАМИЛЬТОНА-ЯКОБИ / ВЕКТОР-ПОТЕНЦИАЛ ГАУССОВА ЛАЗЕРНОГО ПУЧКА / DIRAC WAVE EQUATION / RELATIVISTIC HAMILTON-JACOBI EQUATION / VECTOR POTENTIAL OF GAUSSIAN LASER BEAM

Аннотация научной статьи по физике, автор научной работы — Боровский Андрей Викторович, Галкин Андрей Леонидович

В статье исследуется иерархия теоретических моделей, описывающих ускорение электронов в фокусе лазерного импульса с релятивистской интенсивностью излучения. Последовательность моделей включает волновое уравнение Дирака и его классический предел релятивистское уравнение Гамильтона-Якоби (РГЯ). Рассмотрен пи-комплекс для обезразмеренного линейного уравнения Дирака. Показано, что уравнение Дирака имеет смысл применять только на расстояниях порядка комптоновской длины волны. Для описания явлений в фокусе лазерного импульса уравнение Дирака избыточно, достаточно использовать уравнение РГЯ. Последнее представлено в эйлеровой и лагранжевой постановках. При этом переход к лагранжевой постановке осуществляется взятием градиента от уравнения РГЯ, записанного в эйлеровой постановке. Прослежено как уравнение РГЯ в лагранжевой постановке сводится к релятивистскому уравнению движения отдельного электрона под воздействием силы Лоренца. Получено уравнение РГЯ в эйлеровой постановке для гауссова пучка и линейной поляризации лазерного излучения. Разработанную теоретическую модель предлагается использовать в вычислительных задачах лазерной физики.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по физике , автор научной работы — Боровский Андрей Викторович, Галкин Андрей Леонидович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Mathematical models based on the relativistic Hamilton-Jacobi equation in laser physics

The article considers the hierarchy of theoretical models describing the acceleration of electrons in focus of a relativistic intensity laser pulse. The sequence of models includes the Dirac wave equation and its classic limit the relativistic Hamilton-Jacobi equation (RHJ). The ?-complex intended for the dimensionless linear Dirac equation is considered. It is shown that the Dirac equation can be applied only at the Compton wavelength distances. The Dirac equation is surplus for description of the phenomena in a laser pulse focus, and the RHJ equation should be used instead. The last is examined in Euler's and Lagrange's formulations. The transit to Lagrange's formulations realized by taking the gradient from the RHJ equation written in the Euler's formulation. It is shown that the RHJ equation in Lagrange's formulation reduces to a relativistic equation of a separate electron motion under the Lorentz force law. The RHJ equations received in Euler's formulation for a Gaussian beam and linear polarization of laser radiation. It is suggested that the devised theoretical model be used in laser physics calculable tasks.

Текст научной работы на тему «Математические модели, основанные на релятивистском уравнении Гамильтона-Якоби в лазерной физике»

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ, СИСТЕМНЫЙ АНАЛИЗ

MATHEMATICAL MODELING, SYSTEMS ANALYSIS

УДК 535.3;531.1

DOI 10.17150/1993-3541.2015.25(1).141-148

А. В. БОРОВСКИЙ

Байкальский государственный университет экономики и права, г. Иркутск, Российская Федерация А. Л. ГАЛКИН

Институт общей физики РАН им. А. М. Прохорова, г. Москва, Российская Федерация

МАТЕМАТИЧЕСКИЕ МОДЕЛИ, ОСНОВАННЫЕ НА РЕЛЯТИВИСТСКОМ УРАВНЕНИИ ГАМИЛЬТОНА-ЯКОБИ

В ЛАЗЕРНОЙ ФИЗИКЕ

Аннотация. В статье исследуется иерархия теоретических моделей, описывающих ускорение электронов в фокусе лазерного импульса с релятивистской интенсивностью излучения. Последовательность моделей включает волновое уравнение Дирака и его классический предел — релятивистское уравнение Гамильтона-Якоби (РГЯ). Рассмотрен п-комплекс для обезразмеренного линейного уравнения Дирака. Показано, что уравнение Дирака имеет смысл применять только на расстояниях порядка комптоновской длины волны. Для описания явлений в фокусе лазерного импульса уравнение Дирака избыточно, достаточно использовать уравнение РГЯ. Последнее представлено в эйлеровой и лагранжевой постановках. При этом переход к лагранжевой постановке осуществляется взятием градиента от уравнения РГЯ, записанного в эйлеровой постановке. Прослежено как уравнение РГЯ в лагранжевой постановке сводится к релятивистскому уравнению движения отдельного электрона под воздействием силы Лоренца. Получено уравнение РГЯ в эйлеровой постановке для гауссова пучка и линейной поляризации лазерного излучения. Разработанную теоретическую модель предлагается использовать в вычислительных задачах лазерной физики.

Ключевые слова. Волновое уравнение Дирака; релятивистское уравнение Гамильтона-Якоби; вектор-потенциал гауссова лазерного пучка.

Информация о статье. Дата поступления 30 ноября 2014 г.; дата принятия к печати 26 декабря 2014 г.; дата онлайн-размещения 27 февраля 2015 г.

A. V. BOROVSKY

Baikal State University of Economics and Law, Irkutsk, Russian Federation A. L. GALKIN

General Physics Institute of RAS n.a. A.M. Prokhorov, Moscow, Russian Federation

MATHEMATICAL MODELS BASED ON THE RELATIVISTIC HAMILTON-JACOBI EQUATION IN LASER PHYSICS

Abstract. The article considers the hierarchy of theoretical models describing the acceleration of electrons in focus of a relativistic intensity laser pulse. The sequence of models includes the Dirac wave equation and its classic limit — the relativistic Hamilton-Jacobi equation (RHJ). The n-complex intended for the dimen-sionless linear Dirac equation is considered. It is shown that the Dirac equation can be applied only at the Compton wavelength distances. The Dirac equation is surplus for description of the phenomena in a laser pulse focus, and the RHJ equation should be used instead. The last is examined in Euler's and Lagrange's formulations. The transit to Lagrange's formulations realized by taking the gradient from the RHJ equation written in the Euler's formulation. It is shown that the RHJ equation in Lagrange's formulation reduces to a relativistic equation of a separate electron motion under the Lorentz force law. The RHJ equations received in Euler's formulation for a Gaussian beam and linear polarization of laser radiation. It is suggested that the devised theoretical model be used in laser physics calculable tasks.

Keywords. Dirac wave equation; relativistic Hamilton-Jacobi equation; vector potential of Gaussian laser beam.

Article info. Received November 30, 2014; accepted December 26, 2014; available online February 27, 2015.

© А. В. Боровский, А. Л. Галкин, 2015

141

Izvestiya of Irkutsk State Economics Academy, Известия Иркутской государственной экономической академии.

2015, vol. 25, no. 1, pp. 141-148. ISSN 1993-3541 2015. Т. 25, № 1. С. 141-148. ISSN 1993-3541

MATHEMATICAL MODELING, SYSTEMS ANALYSIS

Одной из актуальных проблем лазерной физики является ускорение электронов в фокусе ультракороткого лазерного импульса. Часто лазерный импульс фокусируют в разреженный газ. Электроны, образующиеся в таком газе, взаимодействуют только с электромагнитной волной и не взаимодействуют друг с другом. Поэтому можно применять модели, описывающие движение отдельных электронов в заданном электромагнитном поле. Перспективной математической моделью для описания движения электронов в заданном поле произвольной интенсивности является релятивистское уравнение Гамильтона-Якоби (РГЯ) [5].

Решения уравнения РГЯ для линейно и циркулярно поляризованной плоской волны без ограничения на величину ее амплитуды представлены в исследовании Л. Д. Ландау и Е. М. Лифшица [Там же]. Параметрические уравнения для траектории электрона, движущегося в поле линейно-поляризованной плоской монохроматической волны, получены и исследованы численно А. В. Боровским [3]. Актуальные исследования движения электронов в полях релятивистской интенсивности, опирающиеся на уравнения релятивистской электродинамики бесстолкновительного электронного газа, описаны А. В. Боровским и

A. Л. Галкиным [2]. Там же можно найти библиографию по данному вопросу до 1995 г.

Альтернативное изучение движения электронов в фокусе гауссова импульса предпринято B. Quesnel и P. Mora. Авторы исследуют возможность введения пондеромоторной силы при релятивистских интенсивностях излучения и условия применимости этой концепции, а также делают вывод об усложненном характере движения электронов вне пределов применимости понятия пондеромоторной силы [11].

Y. I. Salamin, G. R. Mocken, C. H. Keitel анализируют движение одного электрона в трехмерном пространстве вблизи гауссова фокуса [12], а также выводят выражения для полей вблизи фокуса, повторяя исследования

B. Quesnel и P. Mora [11].

Постановка задачи об ускорении электронов, полученных при многофотонной ионизации многозарядных ионов, описана A. Maltsev и T. Ditmire [10]. При этом были использованы формулы для скорости многофотонной ионизации, расчеты проведены методом Монте-Карло. Представлен двухпичковый угловой энергетический спектр. Причину возникновения пичков авторы не обсуждают. Сделан

вывод о необходимости учета продольной компоненты лазерного поля. Приведено сравнение расчетов с этой компонентой и без нее. Кулоновское поле многозарядного иона не учитывается.

S. Masuda, M. Kando, H. Kotaki, K. Nakajima изучают движение одного электрона в трехмерном пространстве вблизи фокуса гауссова импульса, рассматривают необходимость учета продольной компоненты лазерного поля, которая существенно влияет на траекторию движения электронов. Меняют начальную координату электрона. Начальную скорость электрона полагают равной нулю [13]. Набранной статистики в данной работе явно недостаточно, поэтому авторы делают упор на качественных вещах (таких как учет продольной компоненты поля на движение одного электрона), пишут о возможности получения сгустка ускоренных электронов вдоль оси.

S. X. Hu, A. F. Starace cформулировали модель ускорения электрона в суммарном поле лазерного импульса и кулоновском поле многозарядного иона [9]. Наряду с этим изучают начальные условия для численного счета. Расчеты проводятся методом Монте-Карло. Задаются статистические распределения начальных условий для ионизованного электрона. Интервалы начальных параметров разбиваются на равные промежутки, электроны распределяются по интервалам с учетом распределений. Проводится расчет траекторий для каждого электрона и вычисляется остаточная энергия. Смысл работы в том, что многозарядный ион дополнительно удерживает ионизованный электрон в области сильного поля. В результате электрон успевает разогнаться до большей скорости. Модель лазерного поля [9] совпадает с уравнениями для полей в лазерном пучке [11], а в плоском случае — с уравнениями [6]. Зависимость от азимутального угла в данной работе отсутствует.

K. I. Popov, V. Yu. Bychenkov, W. Rozmus, R. D. Sydora изучают движение одного электрона в трехмерном пространстве вблизи фокуса, образованного отражением плоской волны от параболического зеркала [7]. Отраженную волну рассчитывают с использованием метода функций Грина. Такой расчет носит название дифракции Зоммерфельда-Кирхгофа. Для параболической границы получаются интегралы Стрэттона-Чу для полей. Реализован строгий расчет поля вблизи фокуса. Далее авторы рассчитывают движение одного электрона вблизи

A. V. BOROVSKY, A. L. GALKIN

фокуса. Модель движения трехмерная и релятивистская. Делают вывод, что результаты не совпадают с теми, которые получаются для гауссова фокуса.

В работе А. Л. Галкина, В. В. Коробкина, М. Ю. Романовского и О. Б. Ширяева анализируются движение отдельного электрона в поле линейно и циркулярно поляризованного гауссова лазерного импульса, а также особенности траектории электрона, сделан вывод, что траектории имеют мало общего с известными траекториями электрона в поле плоской электромагнитной волны [8]. В последующей их статье описывается излучение электрона, движущегося в поле гауссова лазерного импульса при релятивистских интенсивностях лазерной волны [14].

Далее рассмотрим теоретические основы ускорения электронов лазерным импульсом.

1. Вывод релятивистского уравнения Гамильтона-Якоби из волнового уравнения Дирака. Наиболее точно в современной теоретической физике движение отдельного электрона в заданном поле описывает релятивистское волновое уравнение Дирака [1, с. 146, формула (32.7а)]:

2 / \2

ih d c dt

-| ihV + - Л] - m2c2 +

eh —h + —I.H-1— aE c c

У = 0, (1)

где ф, А — скалярный и векторный потенциалы электромагнитного поля; £ и а — матричные вектора, компоненты которых составлены из матриц Дирака; Е, Н —напряженности электрического и магнитного полей; у — 4-спинор.

Это уравнение получается операторным перемножением двух сопряженных линейных операторов Дирака [1, с.144-148].

Ищем решение уравнения (1) в виде

у = иеь , (2)

где и — медленно изменяющийся 4-спинор; 5(/, г) — функция, которую называют действием.

При применении оператора, стоящего внутри квадратных скобок в формуле (1), к функции (2) возникает выражение, представляющее собой разложение по степеням А. В нашем случае рассмотрим только слагаемое, не зависящее от А. В силу малости величины последующие слагаемые разложения учитывать не будем. Применим первую круглую скобку к выражению (2) один раз

ih д

c dt c

ф| У =

ih vS du ( dS

=_—h

c dt

Второй раз

ih д

+ еф l-y^-dt J c

dS ] 1

— + еф l-y. dt J c

e | (ih d e

cd~~c v y л7dt~cф'х

д5 ) 1 I 1 ( д5 *<Ч- + ^+ еФ|

Аналогично поступим со второй скобкой

/ЙУ + еА^ -еА^

В результате для слагаемого, не зависящего от А, в разложении уравнения (1) получаем

, 2 / ч 2

+ —ф] -(VS- — Л 1 -m2c2

y = 0.

с- / V с

Перед 4-спинором у стоит скаляр, поэтому

сН! + *ф)2 > -СА12 - т2с2 =0 <3)

е А д5

д = ¥5 = р+— А; е =— + еф. с дt

Мы получили классическое уравнение РГЯ. Уравнение (3) приведено в [5, с.68].

Второе слагаемое в разложении уравнения (1) по степеням представляет собой уравнение Паули для спина [1, с. 150, формула (33.7)].

2. Вектор-потенциал гауссова пучка. Теорию гауссовых пучков можно найти в книгах [4; 6]. Гауссовы пучки — это строгие решения параболического уравнения. Целесообразно записывать параболическое уравнение для векторного потенциала. Дело в том, что векторный потенциал электромагнитного поля в кулоновской калибровке напрямую удовлетворяет волновому уравнению, следовательно, для пучков — параболическому уравнению. Электрическое поле находится дифференцированием векторного потенциала по времени и умножением на —1. Это в частности означает, что направления компонент электрического поля и векторного потенциала совпадают с точностью до знака. Магнитное поле находится как ротор от векторного потенциала. Если известно электрическое поле, то векторный потенциал можно найти, интегрируя по времени электрическое поле в выражении для гауссова

Izvestiya of Irkutsk State Economics Academy, Известия Иркутской государственной экономической академии.

2015, vol. 25, no. 1, pp. 141-148. ISSN 1993-3541 2015. Т. 25, № 1. С. 141-148. ISSN 1993-3541

MATHEMATICAL MODELING, SYSTEMS ANALYSIS

пучка и отбрасывая константу интегрирования, зависящую от координат, поскольку среднее от векторного потенциала должно равняться нулю.

Берем выражение для гауссова пучка [11]:

Ex - Eo^еsin(9g);

Р g

(4)

,xPo

Ez - 2Eoв^е p cos^);

P = PoJ1 + -y;

V ZR

< . gf z ï

Фд - юГ — kz + arctgi — i -

zr

-Фо;

IZR ^ ZRp2 (1) . ( z Л

Фд =фд + агс*9^,

где Е0 — амплитуда напряженности электрического поля в начале координат; ю и к — частота и волновой вектор лазерного излучения; zR =лр2Д — рэлеевская длина; р0 — поперечный размер минимальной гауссовой перетяжки в фокальной плоскости; ф0 — значение фазы в момент времени t = 0, в точке z = 0; е — малый параметр.

При записи выражения (4) выбрана декартова система координат с центром в фокальной плоскости г = д/х2 + у2. Оси х и у системы координат лежат в фокальной плоскости. Ось z перпендикулярна фокальной плоскости и совпадает с осью лазерного пучка.

Поскольку

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

1 д_ c dt

A,

то получаем интегрированием по времени

A - A,^е p2cos^g); Р g

(5)

А =-2Д>в ХР0 в ^¡п^).

Р

Частоту, на которую делится выражение, ввели в А0.

Поскольку фаза произвольна, выделяем из нее слагаемое 3п/2. Используя в выражении (5) формулы приведения, получим в результате следующие формулы, которые по своему виду просто совпадают с формулами для электрического поля в гауссовом пучке:

Ax - A,^е p2sin^g); Р g

= 2A,s XP2L е'^ф«).

(6)

Таким образом, в случае линейно-поляризованной волны, вектор потенциал (6) имеет в гауссовом пучке компоненты Ах (большая) и А1 (малая). Компонента Ау = 0.

3. На каких масштабах применимо уравнение Дирака. Интенсивность в центре фокальной перетяжки для гауссова пучка (максимальная интенсивность) составляет

I = —E 2 Io = 8п E°.

(7)

Максимальная амплитуда векторного потенциала связана с максимальной напряженностью электрического поля соотношением

л = —Eo = fo = ^ р A - - , - т Eo.

ю k 2п

(8)

Поскольку

Eo =

8п1п

то, используя выражения (7) и (8), найдем

A-А

o 2п

8п1,

o

(9)

Интенсивность в формуле (9) целесообразно сравнивать с релятивистским значением интенсивности

I-

2 3 2

m c ю

4пе

(io)

Напомним, что при I = 1г масса электрона, осциллирующего в поле циркулярно поляризованной волны, увеличивается в раз [3]. Введем параметр

|П (11)

а-

Ir

который показывает, во сколько раз максимальная интенсивность в фокальной перетяжке гауссова пучка превышает релятивистское значение интенсивности. В уравнении (9) значение 10 разделим и умножим на 1г. В результате, используя формулы (10) и (11), получим

Ao-А

o 2п

-V2 m^ ,/а-

- 2,411-io—

эрг

ед. заряда в СГСЭ

r

r

r

r

2

c

c

r

r

2

A. V. BOROVSKY, A. L. GALKIN

Выпишем далее линейное уравнение Дирака

Y + eAoa(x) - mc dx c

W (x) = 0,

где x = x0, xc, x2, x3; x0 = ct. Вводим нормированную переменную

~ x _ „

X — I X — xxn I

W(x) — 0

и получаем п-комплекс задачи:

y| i dx+^ia(x)l -п2

v(x) — 0;

(13)

_ еХН А _ тСХН

П1 _ ^ П2 _ г. .

П

Целесообразно выбрать в (13)

exн л 1

П _—ГA0 _

СП

отсюда с учетом выражения (12)

_ cй _ Хг:

н eA0 -у/2а'

где Хс _ П/mc _ 3,862-10-11 — комптоновская длина волны. Подставив хн в выражение для п2, найдем

mcx

mc

1

1^2 — — — I-■

П eA0 л/2а Таким образом, п-комплекс для линейного уравнения Дирака (13) составляет

1

п _1; П2 ■ (14)

При условии выполнения (14) константа обезразмеривания хн сравнима с компто-новской длиной волны. Отсюда следует, что уравнение Дирака применимо только на расстояниях порядка комптоновской длины волны. Размеры атома порядка 10-8 см. Даже для атомов релятивистские эффекты будут порядка 10-3 (отношение комптоновской длины волны к размерам атома). В фокусе лазерного импульса размеры составляют порядок длины волны излучения X = 10-4 см. Для таких задач релятивистские поправки к движению электронов, обусловленные уравнением Дирака, будут пренебрежимо малы. Это означает, что рассматривая ускорение электронов в фокусе лазерного импульса, можно применять классическую (некванто-

вую) электродинамику, в частности, уравнение РГЯ (3).

4. Уравнение РГЯ в лазерной физике. Лазерное излучение представляет собой электромагнитное вихревое поле, поэтому-потенциал ф = 0. Уравнение РГЯ (3) принимает вид:

^ _0; (15)

^ д5 д_У5 _ р + -А; е_—.

С ОТ

Уравнение (15) можно разрешить относительно производной по времени

с dS \lVS сA

2 2 + m c .

(16)

Два знака означают возможность существования частиц с положительной и отрицательной энергией. Мы не имеем возможности рассматривать теорию электрон-позитронно-го мира, предложенную Дираком. Рассмотрим движение в лазерном поле свободных электронов. В этом случае в выражении (16) следует выбрать знак «+».

Квадрат вектора означает скалярное перемножение вектора на самого себя. Это есть квадрат модуля вектора, который равен сумме квадратов его компонент. Получается дифференциальное уравнение в частных производных первого порядка по времени со сложной зависимостью от поперечной координаты. В уравнении (16) можно раскрыть квадрат под радикалом

1 dS

-^_\Р5>2 - 2-(А7Я + — А + т2^ . (17) c дТ \ c с

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

В ультрарелятивистском случае, когда в формуле (17) массой покоя частицы можно пренебречь, приходим к уравнению

1 д5

c dt

VS - С A

5. Уравнение РГЯ для гауссова пучка с линейной поляризацией. В гауссовом пучке с линейной поляризацией существуют только две компоненты вектор-потенциала Ах, Аг:

1 д5 _

c дТ "

dS с х 1 (dS1 (dS с х

---Ax | + e | — | + ez |---Az

dx c ) \dy) \dz c

2 2 + m c .

x

H

Izvestiya of Irkutsk State Economics Academy, Известия Иркутской государственной экономической академии.

2015, vol. 25, no. 1, pp. 141-148. ISSN 1993-3541 2015. Т. 25, № 1. С. 141-148. ISSN 1993-3541

MATHEMATICAL MODELING, SYSTEMS ANALYSIS

Раскрывая квадрат, получим

1öS _

с dt ~

( dS еА 12

"läz" - J

с dS s _ — dt

(18)

_2_2 .

Представим также еще одну полезную для дальнейших исследований форму уравнения РГЯ. Для упрощения записи положим далее c = 1, e = 1. Возведем в квадрат круглые скобки в выражении (18)

, 2 2 / 2 , 2

dS_ ду

+ 2ЛХ + 2AZ ^ = т2 + (Лх )2 + (Лг )2.

дх ^

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

6. Вывод из уравнения РГЯ выражения для силы Лоренца, действующей на электрон в электромагнитном поле. Используем систему единиц c = 1, e = 1. Уравнение РГЯ (15) запишем в следующем виде:

§)2-(У5-А)2 -m2 = 0,

где 5 и А — скаляр и вектор, зависящие от координат и времени; 5 = 5(г, I); А = А(г, I).

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

12 = ГЭ51(-ЗУ51 = 2бад кз») КдТ дТ ) дТ Обозначим

е(г, Т) = Щ~Г) и д(г, Т) = У5(г, Т).

Далее

31

У(У5 - А)2 = У ((У 5 - А)(У5 - А)) = = У[(У5)2 - 2(АУ5) + А2];

V (У5)2 =У (У5У5). Для вычисления градиента от скалярного произведения векторов воспользуемся формулой векторного анализа

У(аЬ) = (в-У)Ь+ (ЬУ)а + + ах (Ух Ь) + Ь х (V х а).

В данном случае а = Ь = У5 = д. Ротор от градиента равен нулю УхУ5 = 0. В результате получим У(У5)2 =У(У5-У5) = 2(д У)д. Затем V (А У5) = (А У)д+ (д • У)А + дх (Ух А). Одно из слагаемых выпадает, поскольку ротор от градиента равен нулю УА2 = У(АА) = 2(А У)А + 2А х (Ух А).

Собирая вычисленные производные, найдем

д-д = (д • У)д-(А У)д-(д • У)А-дх дТ

х (Ух А) + (А У )А + 2 А х (Ух А).

В теории электромагнетизма в отсутствие потенциальных полей

Н=Ух А; £=-—.

дт

Вычтем из левой и правой частей уравнения выражение

ЗА

е-

дТ

и введем комбинацию д— А, получим

д(д- А) ЗА

——-- = -е-+

дТ дТ

+ ((д- А)У)(д-А) - (д-А) х Н.

Согласно курсу теории поля [1, с. 63-64] величина представляет собой обобщенный импульс д = У5 = р + А, где р— импульс частицы. Подставив в выражение, получим

др дА

е^7 = -е^7 + (Р^)р-рх Н. дТ дТ

Скорость частицы связана с импульсом соотношением

г = Р.

е

Поделив на е, найдем др дА

-!— =--+ (V• У)р- vхH.

дТ дТ

В уравнении РГЯ можно одновременно поменять знаки у У5 и вектора А или, что то же самое, у векторов р и А, при этом уравнение РГЯ не изменится. Учитывая это, получим

-РР + (V • У)р = E+vх Н; дТ

dp

dt

_ E+vx H.

(19)

Правая часть релятивистского уравнения (19) представляет собой выражение для силы Лоренца, действующей на заряженную частицу.

Таким образом, в данной статье рассмотрена иерархия теоретических моделей,

A. V. BOROVSKY, A. L. GALKIN

описывающих движение электронов в фокусе лазерного импульса. Эта иерархия включает релятивистское волновое уравнение Дирака, описывающее электрон во внешнем поле. Классическим пределом уравнения Дирака является уравнение РГЯ. Показано, что уравнение Дирака для решения задачи является избыточным, так как достаточно анализировать и решать уравнение РГЯ. Предлагается вид уравнения РГЯ для электромагнитного поля лазерного излучения. Такое поле является вихревым и в нем отсутствует электрический потенциал. Выписано уравнение РГЯ для линейно поляри-

зованного лазерного излучения, представляющего собой гауссов пучок. Уравнение РГЯ может изучаться в эйлеровой и лагранжевой постановках. Показан теоретический переход от эйлеровой постановки к лагранжевой, который осуществляется взятием градиента от уравнения РГЯ. Строгое уравнение РГЯ в лагранжевой постановке сводится к релятивистскому уравнению движения для электрона, на который действует сила Лоренца.

Следовательно, уравнение РГЯ может выступать в качестве основы численного моделирования ускорения электронов лазерным импульсом.

Список использованной литературы

1. Берестецкий В. Б. Квантовая электродинамика / В. Б. Берестецкий, Е. М. Лифшиц, Л. П. Питаевский. — М. : Наука, 1980. — 704 с.

2. Боровский А. В. Лазерная физика: рентгеновские лазеры, ультракороткие импульсы, мощные лазерные системы / А. В. Боровский, А. Л. Галкин. — М. : ИздАт, 1996. — 496 с.

3. Боровский А. В. Релятивистские спирали / А. В. Боровский // Применение математических методов и информационных технологий в экономике : сб. науч. тр. — Иркутск : Изд-во БГУЭП, 2011. — Вып. 10. — С. 14.

4. Вайнштейн Л. А. Электромагнитные волны, радио и связь / Л. А. Вайнштейн. — М. : Радио и связь, 1990. — 440 с.

5. Ландау Л. Д. Теория поля / Л. Д. Ландау, Е. М. Лифшиц. — М. : Наука, 1973. — 504 с.

6. Ярив А. Квантовая электроника / А. Ярив. — М. : Советское радио, 1980. — 488 с.

7. Electron vacuum acceleration by a tightly focused laser pulse / K. I. Popov, V. Yu. Bychenkov, W. Rozmus, R. D. Sy-dora // Physics of Plasmas. — 2008. — Vol. 15. — Р. 013108. — DOI : 10.1063/1.2830651.

8. Dynamicsofan electron driven by relativistically intense laser radiation / A. L. Galkin, V. V. Korobkin, M. Yu. Roma-novsky, O. B. Shiryaev // Physics of Plasmas. — 2008. — Vol. 15, no. 2. — P. 023104. — DOI :10.1063/1.2839349.

9. Hu S. X. Laser acceleration of electrons to giga-electron-volt energies using highly charged ions / S. XHu,A. F. Sta-race // Physical Review E. — 2006. — Vol. 73, no. 6. — P. 066502. — DOI : 10.1103/PhysRevE.73.066502.

10. Maltsev A. Above Threshold Ionization in Tightly Focused, Strongly Relativistic Laser Fields / A. Maltsev, T. Dit-mire // Physical Review Letters. — 2003. —Vol. 90, no. 5.

11. Quesnel B. Theory and simulation of the interaction of ultraintense laser pulses with electrons in vacuum / B. Quesnel, P. Mora // Physical Review. — 1998. —Vol. 58, no. 3. — P. 3719-3732.

12. Salamin Y. I. Electron scattering and acceleration by a tightly focused laser beam / Y. I. Salamin, G. R. Mocken, C. H. Keitel // Physical Review Special — Accelerators and Beams. — 2002. — Vol. 5, no. 10. — P. 9-22.

13. Suppression of electron scattering by the longitudinal components of tightly focused laser fields / S. Masuda, M. Kan-do, H. Kotaki, K. Nakajima // Physics of Plasmas. — 2005. — Vol. 12. — P. 013102. — DOI : 10.1063/1.1815000.

14. Temporal and Spectral Characteristics of the Radiation Generated by an Electron Driven by a Relativistically Intense Laser Field / A. L. Galkin, M. P. Kalashnikov, V. V. Korobkin, M. Yu. Romanovsky, O. B. Shiryaev, V. A. Trofimov // Contributions to Plasma Physics. — 2011. — Vol. 51, no. 5. — P. 482-490. — DOI : 10.1002/ctpp.201110020.

References

1. Berestetskii V. B., Lifshits E. M., Pitaevskii L. P. Kvantovaya elektrodinamika [Quantum electrodynamics]. Moscow, Nauka Publ., 1980. 704p.

2. Borovskii A. V., Galkin A. L. Lazernaya fizika: rentgenovskie lazery, ul'trakorotkie impul'sy, moshchnye lazernye sistemy [Laser Physics: X-ray lasers, ultrashort pulses, high-power laser systems]. Moscow, IzdAt Publ., 1996. 496p.

3. Borovskii A. V. Relativistic Spirals. Primenenie matematicheskikh metodov i informatsionnykh tekhnologiy v ekono-mike [Mathematical methods and information technologies application in economics]. Irkutsk, Baikal State University of Economics and Law Publ., 2011. Iss. 10, pp. 14. (InRussian).

4. Vainshtein L. A. Elektromagnitnye volny, radio isvyaz' [Electromagnetic waves, radio and communication]. Moscow, Radio I svyaz Publ., 1990. 440p.

5. Landau L. D., Lifshits E. M. Teoriya polya [Field theory]. Moscow, Nauka Publ., 1973. 504 p.

6. Yariv A. Kvantovaya elektronika [Quantum Electronics]. Moscow, Sovetskoe radio Publ., 1980. 488 p.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

7. Popov K. I., Bychenkov V. Yu., Rozmus W., Sydora R. D. Electron vacuum acceleration by a tightly focused laser pulse. Physics of Plasmas, 2008, vol. 15, pp. 013108. DOI: 10.1063/1.2830651.

8. Galkin A. L., Korobkin V. V., Romanovsky M. Yu., Shiryaev O. B. Dynamics of an electron driven by relativistically intense laser radiation. Physics of Plasmas, 2008, vol. 15, no. 2, pp. 023104. DOI : 10.1063/1.2839349.

9. Hu S. X., Starace A. F. Laser acceleration of electrons to giga-electron-volt energies using highly charged ions. Physical Review E, 2006, vol. 73, no. 6, p. 066502. DOI : 10.1103/PhysRevE.73.066502.

Izvestiya of Irkutsk State Economics Academy, Известия Иркутской государственной экономической академии.

2015, vol. 25, no. 1, pp. 141-148. ISSN 1993-3541 2015. Т. 25, № 1. С. 141-148. ISSN 1993-3541

MATHEMATICAL MODELING, SYSTEMS ANALYSIS

10. Maltsev A., Ditmire T. Above Threshold Ionization in Tightly Focused, Strongly Relativistic Laser Fields. Physical Review Letters, 2003, vol. 90, no. 5.

11. Quesnel B., Mora P. Theory and simulation of the interaction of ultraintense laser pulses with electrons in vacuum. Physical Review, 1998, vol. 58, no. 3, pp. 3719-3732.

12. Salamin Y. I., Mocken G. R., KeitelC. H. Electron scattering and acceleration by a tightly focused laser beam. Physical Review Special — Accelerators and Beams,2002, vol. 5, no. 10, pp. 9-22.

13. Masuda S., Kando M., Kotaki H., Nakajima K.Suppression of electron scattering by the longitudinal components of tightly focused laser fields.Physics of Plasmas, 2005, vol. 12, pp. 013102. DOI : 10.1063/1.1815000.

14. Galkin A.L., Kalashnikov M. P., Korobkin V. V., Romanovsky M. Yu., Shiryaev O. B., Trofimov V. A. Temporal and Spectral Characteristics of the Radiation Generated by an Electron Driven by a Relativistically Intense Laser Field. Contributions to Plasma Physics, 2011, vol. 51, no. 5, pp. 482-490. DOI : 10.1002/ctpp.201110020.

Информация об авторах Боровский Андрей Викторович — доктор физико-математических наук, профессор, кафедра информатики и кибернетики, Байкальский государственный университет экономики и права, 664003, г. Иркутск, ул. Ленина, 11, e-mail: [email protected].

Галкин Андрей Леонидович — доктор физико-математических наук, заведующий сектором, Институт общей физики РАН им. А. М. Прохорова, 119991, г. Москва, ул. Вавилова, 38, e-mail: [email protected].

Библиографическое описание статьи Боровский А. В. Математические модели, основанные на релятивистском уравнении Гамильтона-Якоби в лазерной физике / А. В. Боровский, А. Л. Галкин // Известия Иркутской государственной экономической академии. — 2015. — Т. 25, № 1. — С. 141-148. — DOI : 10.17150/1993-3541.2015.25(1).141-148.

Authors

Andrei V. Borovsky — Doctor habil. (Physical and Mathematical Sciences), Professor, Baikal State University of Economics and Law, 11 Lenin St., 664003, Irkutsk, Russian Federation,e-mail: [email protected].

Andrei L. Galkin — Doctor habil. (Physical and Mathematical Sciences), Professor, General Physics Institute of RAS n.a. A. M. Prokhorov, 38 Vavilova St., 119991, Moscow, Russian Federation,e-mail: [email protected].

Reference to article

Borovsky A. V., Galkin A. L. Mathematical models based on the relativistic Hamilton-Jacobi equation in laser physics. Izvestiya Irkutskoy gosudarstvennoy ekonomicheskoy ak-ademii = Izvestiya of Irkutsk State Economics Academy, 2015, vol. 25, no. 1, pp. 141-148. DOI: 10.17150/1993-3541.2015.25П).141-148. (In Russian).

i Надоели баннеры? Вы всегда можете отключить рекламу.