Научная статья на тему 'Ряды Каптейна и излучение фотона'

Ряды Каптейна и излучение фотона Текст научной статьи по специальности «Физика»

CC BY
77
22
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
РЯДЫ КАПТЕЙНА / ИЗЛУЧЕНИЕ ФОТОНА

Аннотация научной статьи по физике, автор научной работы — Никишов А.И.

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

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

Похожие темы научных работ по физике , автор научной работы — Никишов А.И.

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

Текст научной работы на тему «Ряды Каптейна и излучение фотона»

УДК 51-73:535.14

РЯДЫ КАПТЕЙНА И ИЗЛУЧЕНИЕ ФОТОНА

А. И. Никиттюв

Получены суммы нескольких видов рядов Каптейна и найдены соотношения, между ними. Особое внимание уделено рядам,, суммы которых выражаются, через трансцендентные функции. Найдены, их асимптотики. В ряде случаев получены интегральные представления, и приведены к виду, удобному для, численны,х расчетов. Найдено время, жизни электрона, в постоянном, электромагнитном поле до излучения фотона. Отмечено, что это время, увеличивается, с энергией ■ультрарелятивистского электрона, в сильном постоянном, поле.

Ключевые слова: ряды Каптейна. излучение фотона.

Ряды Каптейна. ЛИН6ИНЫ6 по функциям Бесселя, мы называем линейными. Ряды, суммы которых даются алгебраическими функциями, называем алгебраическими. В противном случае ряды называем трансцендентными. Линейные трансцендентные ряды. Рассмотрим ряд

П

Ж ^^ г, <Х

V — Jm(mx) = — deY 1[ae-i(e-x sin e)]m

m=l m=l

—П

П

= - — í de ln[l - ae—i(e—xsin0)] = (1)

2n J

— П

П

= - — í de{ln[1 - ae—i(e—xsin0)] + ln[1 - aet{e—xsine)]} = 2n J 0

П

=--deln[1 + a2 - 2a cos(e - x sin e)]. (2)

2n J

0

Левую часть равенства в (1) и правую часть последнего равенства можно представить

x

ФИАН, 119991 Россия, Москва, Ленинский пр-т, 53; e-mail: nikishov@lpi.ru

одинаковыми степенями х даёт выражения для интегралов по в, зависящих от параметра а Подробнее см. [1]. Заметим, что в [1] в разделе 2 в суммах по нечётным п (где суммирование идёт по 2п +1) суммирование должно начи наться от п = 0 а не о тп = 1. Особый интерес представляет случай |а| = 1. Тогда из (1) и (2) имеем

п П

V— Jm(mx) = —- Id- ln[2(l - cos(9 + y))] = - - í d- ln 2sin( - - - sin -' m 2ж / ж / V 2 2

m=l o o

• (3)

Здесь y = — | sin-.Используя разложение Jm (mx) по степен ям x, получим

e l x x2 x3 x4 23x5 llx6 841x7 151x8

2^mJm(mx)= 2 + 22+ 23 + + 27^3 + 24 • 3 • 5 + 29 • 32 • 5 + 24 • 32 • 5 • 7 + '" • (4)

m=1

Дифференцируя это выражение, найдём

23

± J_ OJ_ J_ —jO

2 + 2+1г + У + 3 ' 23 -5 ' 29 -32 -5 ' 2 • 32 -5-7

^ )= 1 x 3x2 x3 23 • 5x4 llx5 841 • 7x6 151x7

Jm(mx) = 2 + 2 + ^ + 27 . 3 + + 29 . 32 . 5 + 2 . 32 . 5.7 + " ' • (5)

m=1

Дифференцируя (3), найдём

e l í (e x \

У J'm(mx) = — d- cot I - — — sin - ) sin-• (6)

, 2ж J \2 2 )

m=1 0

Для 1 — x2 << 1, когда meff >> 1, имеем (см (9.3.43) и (10.4.16) в [2])

j;(mx) = - (m)2/3 A,' ((mУ d - x2)) = ^K2/3 (£(1 - x2)3/2) • (7)

Используя это в левой части (5) и заменяя суммирование интегрированием, получим

J2Jm (mx) = 2^ J2n(2nx) = ц _ x2)1/2 , 1 - x2 << 1 (8)

т — 1 п—1 ^ '

Здесь использовано соотношение (6.561.16) в [3]

dxxa-1Kv(ax) = 2a-2a-ar(г() • (9)

Дифференцируя (8), найдём

<x

Va

J2mJm(mx) = ц _ x2)3/2 ' 1 - x2 << 1 (l0)

m— 1 ^ /

Заметим теперь, что формулу (6) можно записать в виде

п

^ 1 If

Е Jm(mx) = - x + 2nd9 m=1 0

• л , Л 2 ( в 6

sin в I cot f--+—

<<) x \в — x sin в x(c + в2)

+

6 n ex i - x

+--arctan , ю =---sin e, c = 6-. (11)

nx2\Jc \Jc 2 2 x

Интеграл здесь порядка единицы и может быть табулирован. Для 1 - x2 << 1 главный

член в правой части таков

6 п 6 п 3 3 . . arctan ^ & —= ^ & W --2, (12)

п^/с \[с пуС 2 ч/с VI _

как это и должно быть согласно с (8). Аналогично (11) найдём

J]j2„(2nx)= /

n=1 £

п

nx) = I de sin в cot ф =---1--de

2x 2n J 00

sin в ( cot ф — ^ ) +

1 / в 6

+ -

x \e - xsine x(c + e2)

Дифференцируя (6), получим

3 П

+--arctan ^, ф = в — x sin в. (13)

nx2\J c \Jc

5>J¿(mx) =-П /W^Y, < = 2 — fsin в. (14)

4n J0 V sin <fj 2 2

m=l

m

d +— ) J'(mx) = ( — — 1 ) mJ(mx), (15)

dx x x2

получим

d 1\ ^^ , /1

dx + x)EJ'(mx)=( x2 — mJ (mx). (16)

m=l m=l

Аналогично (1) и (2) получим

п

11

V— J2n(2nx) =--Ыв ln(2 sin ф), ф = в — x sin в. (17)

' 2n 2n ./

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

Дифференцируя, имеем

J2n(2nx) = 2- I dв cot ф sin в, ф = в — x sin в. (18)

n=1 о

п

1

n=1 0

Дифференцируя ещё раз. получим

1 Г чт2 в

V 2п3"(2пх) = — ¿в. ^ 2 2ж} чт2 Ь

п=1 о

(19)

Далее легко найдём

V ^ ^ , ч о 7 X 239X

2п32п(2пх) = х2 + — + 22 _ 3 +

7х4 239х6 1481х8 292223х10

+

п=1

3 22 -3-5 22 -32 -7 26 -34 -7

+

(20)

х

1

24

хх

х а 1 , , х х

/-^2п ( х = 22 + 2^3 +

11х6

п=1

+

151х8

+

15619х10

22 22 3 24 3 5 24 32 5 7 28 34 5 7

+

(21)

Дифференцируя это выражение, найдём

. , . х х3 11х5 151х7 15619х9 3-2п(2пх) = - + — + + . „2 . „ + .7 „ + • •

п=1

2 3 23 -5 2- 32 • 5- 7 27 • 34 -7

(22)

Учитывая (8), можно получить формулу, которая при разложении в ряд по степен-х формулу (22):

^ 32п(2пх)

х

п=1

2^1-х2

х2 11х4 59х6 14971х8 1 + + 23 • 3 • 5 + 24 • 32 -7 + 27 • 34 -5-7 +

. (23)

Далее нам потребуется выражение

в

¿х У^ п32п(2пх) = —=--

п=1 2 ^л/3(1 - в2)3/2

1 - в2 « 1.

(24)

Дифференцируя его по в, получим выражение

Е

п=1

п32п(2пв)

л/3

4(1 - в2)5/2;

1 - в2 « 1,

(25)

которое можно проверить, используя формулу (40) ниже.

х

У^ 2п32п(2пх)

х

п=1

(1 - х2)5/2

24

хх

1 - — + +

5х6

103х8

+ ^ „ +

2 3 23 5 24 32 7 27 34 7

(26)

Сходным образом найдём выражение для интеграла суммы трансцендентного ряда

в

ß 3

dxJ2nJ2n (2nx) = ß2)3/2

n=1 ^ '

1 -

ß2 ß4

19ß6

809ß8

2 • 5 23 • 7 24 -33 -7 27 • 33 -7-11

Билинейные алгебраические ряды. Начнём с соотношения, см. [5]

d d (1 — х2) — Jn (пх) = — [x2Jn2(nx dx dx

(27)

(28)

Оно легко проверяется выполнением в нём дифференцирований и использованием уравнения Бесселя. Умножая (28) на пи и суммируя по п, получим

d

d

(1 — х2) dX Jn (пх) = dx [x2Efnv Jn2(nx

dx

Полагая v = — 2, имеем

(1 — x2) d sr -1 J2(nx) = d

dx n2 dx

В левой части используем формулу Нильсона

х Si 2 Jn (nx) n2

(29)

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

(30)

1ТЛГ 1 j2{ \ _ x_

S1 2 Jn (nx) = А ■

n2 4

Выполняя затем дифференцирование в левой части, получим

2 x d (1 — х ) — = —— 2 dx

х Si Jn (nx) 1 n2 n

(31)

(32)

x

оо 1 а 1 х2

Jn (nx) = 7 n2 4 8

(33)

Заметим, что формулу Нильсона (31) легко получить из формулы (17.33(2)) в [4]:

1 x2

Sr n2 J2n(2nx) = у ■

Заменяя здесь х ^ х sin (интегрируя по р и учитывая (39) (см. ниже), получим (31).

Далее дифференцируем Jn(v)'-

d d2

Jv) = 2Jn(v) Jn (V), ^ Jn (v) = 2J2n(v) + 2Jn(v)JKv) =

p

2J'2n(y) + 2Jn(y)

{% - 0 Jn(y) - 1 J'n(y) Vy2 / y

Отсюда (см. (40.4) в [6]

(34)

J' П(у)

1_(1_ i n2

2 V ydy dy2) V y2

•Ш-

(35)

Полагая y = nx, умножая это соотношение на nv и суммируя по n, получим

Я?п"З'Ку) = 2 (X-X + -X) + (1 - Х2) ^^2(пх). (36)

Билинейные трансцендентные ряды. Имея дело с разложениями в степенные ряды, мы используем соотношение (8.442.1) в [3]:

Jn(y) = s0°°(-i)i

(2n + 2s)!y2(n+s)

s!22(n+s)(2n + s)![(n + s)!]2'

(37)

Аналогичное разложение для .]'П(пх) можно получить из (36), но проще использовать (28) или (29), если нужна сумма. Используя (37), получим

2 _ x2 ( 7x2 239x4 7435x6 292223x8 nJn(nx) = T\1 + 1^ + + + 213 • 32 +

(38)

Значительно проще получить его из (20). используя (6.681.9) в [3], имеющую вид

п/2

J^(nx) = — dpJ2n(2nx cos p). n I

(39)

Для 1 — x2 << 1 используем соотношения (9.3.35) и (10.4.14) в [2]

Jn(nx)

VT—x

л/3

K

1/3

П (1 — x2 )3/2

(40)

заменяем суммирование интегрированием и используем (6.576 4) в [3]. Тогда получим

£°°nJn(nx)l,

1"ж2<<1 пл/33 (1 — x2)2"

(41)

Аналогично (23) получим

S0°nJ2(nx) =

x2

4(1 — x2)2

1- x,-

x4

5x6

23x8

22 25 • 3 28 • 32 213 • 32

(42)

1

Аналогичным образом найдём

,2 ( , 1 5х2 127х4 3133х6 101887х8 Е1 пЗ П{пх) _ — + ~2Т + 27 • 3 + 210 • 32 + 215 • 32 +

Используя (7)сга _ 2п и уравнение (6.576 4) в [3], получим

ЕГпЗ 'П(пх)

1

1-ж2<<1 ^^31 — х'

Аналогично соотношению (42) найдём

ЕГпЗ 'П(пх)

1

4(1 - х2)

С помощью (37) получим

1 х2 7х4 85х6 1631х8 + 2^ + 25 • 3 + 28 • 32 + 213 • 32 +

1 х2

ЕГ - З2(пх) _ п 4

' х2 11х4 151х6 15619х8 + 2^ + 2^3 + 28 • 32 + 213 • 32 • 5 +

Намного проще получить его. используя (21) и (39).

С помощью (39) и (3) получим сумму ряда в интегральной форме:

П П

ЕГ- З2(пх) 1 п п

--2 ¿ф ¿61п[2вт(6 — х вт ф вт!

П2

(43)

(44)

(45)

(46)

(47)

Дифференцируя его. получим

1

П П

ЕГ2Зп(пх)З(пх) _ — ¿ф ¿6 со1(6 — х вт ф вт 6) вт ф вт 6.

П2

(48)

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

00

Вероятность излучения, фотона. В классическом случае электрон в постоянном магнитном поле или в поле циркулярно поляризованной плоской волны движется по круговой орбите (в подходящей системе координат). Интенсивность излучения даётся формулой Шотта. см. §74 в [7]:

е2 . . . . 2 „

(49)

¿1,

п--(пшя№п2 6З2(пв сое 6) +в2З'п(пв сое 6)Ш, в _ -.

2пс с

Проинтегрированная по углам интенсивность равна

/

ш2н

в2пЗ2п (2пв) — (1 — в2) ¿хп2З2п (2пх)

(50)

в

В случае магнитного поля шн = —\/1 — в2 •

тс

Вероятность в единицу времени получается делением интенсивности на Нншн, шн есть частота первой гармоники: ¿Рп = ¿1п/Киш. Так (49), делённая на Кишн и просуммированная по и, приводит к рядам (42) и (45). Полная вероятность в единицу времени равна

2е2

Р = = -— ин

hv

в

в2E?J2п(2пв) - (1 - в2) J dxE^nJ2n(2nx)

0

(51)

и

ческим рядам, а суммирование вероятностей к трансцендентным. Вероятность того, что в течение интервала времени Ь не произойдёт излучения фото на, равна ехр(—РЬ). Используя (23) и (27), получим

Р

2^ ин в! hv - в2 3

1 3^ 41в4 52 • 11в6 28121в8 + ^5 + 23 • 5-7 + 24 • З3 • 7 + 27 • З2 • 5-7-11 +

(52)

В случае плоской циркулярно поляризованной волны ин есть частота волны. Для 1 — в << 1 получим

e2 5 eH e2 , ч

Р = ---=—, — = а = 1/137. (53)

hc 2V3 mc hc

Р

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

e2

длине формирования излучения. Тогда [см. (27') и (27") в [8], где 4—— = а = 1/137]

w 5am V(eF^y Pv )2

W = 273Х' Х = тз . (54)

Здесь W - вероятность в единицу собственного времени, т = \¡1 — в2t - тензор

eHv

поля, p - импульс электрона, h = c = 1, Pt = Wt. Для магнитного поля х =

и мы видим, что (53) согласуется с (54): Ш ^1 — в2 = Р.

Интересно заметить, что в квантовом случае, когда х >> 1

т2л/Г—в2

W = ГО)™ (3х)2/з. (55)

Это означает, что лабораторное время до излучения фотона возрастает в этом случае

[1] A. I. Xikishov, arXiv:l.302.0978 [math.pli] 5 Feb 2013.

[2] M. Абрамовиц. И. Стиган. Справочник по математическим функциям, (М.. 1966).

[3] И. С. Градтптейн. И. М. Рыжик. Таблицы интегралов рядов и произведений (М.. 1962).

[4] G. X. Watson. A Treaty on the Theory of Bessel Functions (Cambridge Univ. Press. Cambridge. 1966).

[5] G. A. Schott. Electromagnetic radiation (Cambridge Univ. Press. Cambridge. 1912).

[6] Д. Д. Иваненко. А. А. Соколов. Классическая, теория, поля, (M.. 1951).

[7] Л. Д. Ландау Е. М. Лифтдид, Теория, поля, (М.. Наука. 1973).

[8] А. И. Никитпов, В. И. Ритус, ЖЭТФ 46, 776 (1964).

как

ЛИТЕРАТУРА

Поступила в редакцию 4 марта 2014 г.

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