УДК 534.1 ББК В313
Владимир Евгеньевич Холодовский,
кандидат физико-математических наук, доцент, Брянский государственный университет им. И. Г. Петровского (Брянск, Россия), e-mail: [email protected] Александр Алексеевич Сидоров, кандидат физико-математических наук, доцент, Брянский государственный университет им. И. Г. Петровского (Брянск, Россия), e-mail: [email protected]
Поток энергии и сила реакции на излучение внутриатомного диполя
На основе уравнений Максвелла получена формула, выражающая энергию излуче-ния подвижного диполя, и вычислена сила реакции на излучение диполя, размеры плеча которого сравнимы со среднеквадратичными смещениями атомов кристаллов.
Ключевые слова: уравнение Максвелла, поток вектора Умова-Пойтинга, диполь, энергия излучения, сила реакции на излучение.
Vladimir Yevgenyevich Kholodovskii
Candidate of Physics and Mathematics, Associate Professor, I.G. Petrovsky Bryansk State University, (Bryansk, Russia), e-mail: [email protected] Alexander Alekseevich Sidorov Candidate of Physics and Mathematics, Associate Professor,
I.G. Petrovsky Bryansk State University, (Bryansk, Russia), e-mail: [email protected]
Energy Flow and Force Reaction to the Emission of the Intra-atomic Dipole
Based on Maxwell’s equations, the article presents a formula expressing the radiation energy of the moving dipole, and the calculation of the force reaction to the dipole radiation whose shoulder size is comparable to the mean square displacement of atoms in the crystal.
Keywords: Maxwell’s equations, flow of Umov-Poynting vector, dipole, radiation energy, strength of the reaction to radiation.
В настоящей работе получена формула, выражающая полную энергию излучения диполя за единицу времени через поверхность сферы, радиус которой существенно больше плеча диполя. На основе этой формулы найдено выражение для силы реакции внутриатомного диполя на собственное излучение. При этом под внутриатомным диполем понимается система из двух зарядов, один из которых несёт остов атома, а другой - центр заряда его внешней электронной оболочки. Сле-дует отметить, что в современной литературе, например [1; 2], формула, выражающая энергию излучения диполя применима только для достаточно больших расстояний. В случае малых рассто-яний, например, при исследовании излучения внутриатомных диполей вещества такая формула не годится. В работе [3] показано, что при правильном учёте силы реакции на излучение внутриатомных диполей кубических кристаллов возможно построение динамической модели в адиабатическом приближении, дающей хорошее согласование с экспериментом. В работе [4] с учётом силы реакции внутриатомных диполей произведён расчёт тепловой энергии движения электронов и получена поправка на электронный вклад в теплоёмкость металлов, хорошо согласующаяся с экспериментальными данными.
Рассмотрим диполь, совершающий колебания вблизи своего положения равновесия. Обозначим через р(t), p'(t), р"(t) соответственно плечо и его первую и вторую производные по времени, а через q - заряд диполя; тогда его момент будет равен qp(t). Пусть г - радиус-вектор точки наблюдения относительно положения равновесия диполя, г - его длина, а е = г/г - единичный направляю-щий вектор. Будем считать, что расстояние до точки наблюдения существенно больше возможных
размеров плеча диполя.
© В. Е. Холодовский, А. А. Сидоров, 2012
141
Как следует из уравнений Максвелла, напряжённости электрического и магнитного полей диполя в точке наблюдения в момент времени Ь выражаются соответственно формулами:
где с - скорость света в вакууме, а р, р' и р" вычисляются в момент времени і' = і — г/с. Как известно [1], поток энергии, излучаемой диполем в единицу времени через произвольную ориенти-рованную поверхность, есть поток вектора Умова-Пойнтинга
Следовательно, количество энергии, излучаемой диполем за единицу времени через сферу а радиуса г в направлении внешней нормали равно:
Вычисление вектора 8 сводится к вычислению двойных векторных произведений, в которых первый множитель есть числитель одного из слагаемых формулы (1), а второй множитель - числитель одного из слагаемых формулы (2). Применяя к каждому из шести возможных произведений формулу двойного векторного произведения, приходим к равенствам:
В научной литературе, посвященной подобным задачам, выражение для ТУ содержит лишь первое слагаемое формулы (4). Очевидно, если радиус сферы, через которую вычисляется поток ИЛ достаточно велик, то остальными тремя слагаемыми можно пренебречь. Если же речь идёт о внутриатомном диполе, одним полюсом которого является ионный остов атома, а другим - центр заряда его внешней электронной оболочки и при этом радиус сферы г сравним с радиусом атома, то
Е = д {е(ер") - р" + Зе(ер') - р' + Зе(ер) - р
(1)
(2)
(3)
а
(Зе(ер) - р) х (р х е) = 2р'(ер) - Зе(ер)(ер') + е(рр'),
(Зе(ер') - р') х (р' х е) + (Зе(ер) - р) х (р" х е) = ^(Зе(ер) - р) х (р' х е),
а
Ц>[2р"(еР') “ 4е(ер')(ер") + 2е(р'р")]сго- = ^Цр'2)',
и
и
Таким образом, вычисление интеграла (3) приводит к формуле
(4)
основной вклад в формуле (4) будет уже давать последнее слагаемое. Отметим, что в этом случае можно считать, что Ь' = t.
Действительно, оценим порядки слагаемых в (4), считая, что внутриатомный диполь колеблется по гармоническому закону
р = V зт(о;4 + ф).
где ш - частота, ср - начальная фаза, а V - постоянный множитель, определяющий амплитуду и направление поляризации плеча диполя. Если считать, что внутриатомный диполь совершает теп-ловые колебания, то частота будет измеряться в терагерцах так, что ш « 1012сек-1. Пусть радиус г сферы, через которую вычисляется поток, имеет порядок радиуса атома. Тогда г « 10_1Ом, при этом плечо внутриатомного диполя, сравнимое со среднеквадратичным смещением атома, имеет по-рядок 10~14м, что существенно меньше г. Вычисляя порядки слагаемых, стоящих в скобке формулы (4), получаем
й ^ й 10^, _|_(р«у » * и**
^(р2)'* ш12-2.
Из приведенных оценок видно, что при рассматриваемых условиях подавляющий вклад в выраже-нии (4) дает последнее слагаемое. Поэтому, при расчете излучения внутриатомного диполя можно считать, что
цг=£г^р2у‘_________________________________________®
Пользуясь формулой (5), вычислим силу Ег реакции диполя на его излучение через поверхность сферы радиуса г. Найдем работу, совершаемую этой силой при движении плеча диполя по некоторому пути Ь за промежуток времени от до t. Эта работа равна
_______________Г__________________
Ег(р)сгр = J ¥г(р(і))р'(і)йі.
С другой стороны, энергия, теряемая диполем за счёт излучения через сферу радиуса г за тот же промежуток времени, равна
і
2
іо
Следовательно, справедливо равенство
£ І ~2
Поскольку промежуток интегрирования выбирается произвольно, отсюда следует, что
2 (Р1
^
Для кристаллов величина г, вообще говоря, зависит от температуры и может быть вычислена из уравнения термодинамического равновесия [4].
Список литературы
1. Измайлов С. В. Курс электродинамики. М.: Учпедгиз, 1962.
2. Тамм И. Е. Основы теории электричества. М.: Наука, 1976.
3. Холодовский В. Е., Мачихина И. О. Принцип длинных волн и дисперсионные соотношения для кубических решёток в модели диполь-дипольных взаимодействий // Известия Самарского научного центра РАН. 2009. Т. 11. № 5. С. 49-55.
4. Холодовский В. Е., Мачихина И. О., Кульченков Е. А. Поправка на электрон-ный вклад в теплоёмкость металлов в модели Ван-дер-Ваальсовских взаимодействий // Вестник БГТУ. 2010. № 4. С. 115-123.
Статья поступила в редакцию 28.01.2012 г.