14
.doc
5.1 Уравнения Минковского Формально уравнения Максвелла для движущейся среды могут быть получены усреднением соотношений (, (1.1) , (1.2) , (1.3) , (1.4) . (1.5)) с учетом преобразований Лоренца . (2.5) , (2.6) . (2.8) для скоростей, полей, плотностей заряда и тока. Однако даже при движении с нерелятивистскими скоростями, удобнее, как это было отмечено Г. Минковским, рассматривать релятивистскую задачу, а затем, при необходимости, переходить к нерелятивистскому пределу. |
5.2элементами которого являются компоненты макроскопических полей. Для этого тензора справедливо уравнение . (2.13) , которое эквивалентно первой паре (1.28) и (1.29) уравнений Максвелла в среде. С другой стороны, будучи справедливой как для неподвижных, так и для движущихся сред, вторая пара уравнений (1.30) ,(1.31) должна сохранять свой вид при преобразованиях Лоренца. Для поля в пустоте векторы D и H совпадают с векторами Е и В. Поэтому для обеспечения релятивистской инвариантности уравнений (1.30) и (1.31) необходимо, чтобы компоненты векторов D и H преобразовывались как компоненты 4-тензора Нik, построенного аналогично 4-тензору Fik вида (2.10): . (2.15) Введенный тензор удовлетворяет уравнению, аналогичному уравнению . (2.14) для поля в пустоте и являющемуся |
5.3 Граничные условия Кроме материальных уравнений в случае движущейся среды нужно переформулировать и граничные условия к уравнениям поля. Граничные условия в сопутствующей системе координат имеют вид (1.32) Е1 = Е2. [n(H1 – H2)] = 4i/c. (1.33) B1n = B2n. (1.34) D2n – D1n = 4. (1.35)только в том случае, когда частные производные по времени, входящие в уравнения Максвелла, являются ограниченными величинами, как это имеет место в неподвижной среде. Если же тело движется, то в тот момент, когда его граница проходит через точку наблюдения, поле в этой точке меняется скачком, и частные производные D/t и В/t обращаются в бесконечность. Для того чтобы получить граничные условия для тангенциальных компонент напряженностей электрического и магнитного полей, перейдем к системе отсчета, движущейся вместе с данным элементом поверхности тела. Проекцию скорости этого элемента на нормаль n к поверхности тела обозначим vn. В сопутствующей системе справедливы обычные условия (1.32) и (1.33) для тангенциальных компонент Е и Н. E + [v´B]/c и H + [D´v]/c. Проецируя их на плоскость, перпендикулярную к n и учитывая условия (1.34), (1.35), а также формулы (2.21) и (2.22) в первом порядке малости по отношению v/c, получим: [n(Е1 – Е2)] = vn(B2 – B1)/c vn(2 – 1)H/c (2.23) [n(H1 – H2)] = 4pi/c –vn(D2 – D1)/c 4pi/c –vn(2 – 1)E/c. (2.24) Отметим, что если тело движется так, что его поверхность не смещается в перпендикулярном к самой себе направлении, например при повороте тела вращения вокруг оси, то vn = 0. В этом случае граничные условия (2.23) и (2.24) сводятся к обычным условиям для тангенциальных компонент (1.32) и (1.33). |
|
|
6.1 Комплексная диэлектрическая проницаемость Дисперсия, как правило, связана с внутренними свойствами материальной среды, обычно выделяются частотная (временная) дисперсия, когда поляризация в диспергирующей среде зависит от значений поля в предшествующие моменты времени (память), и пространственная дисперсия, когда поляризация в данной точке зависит от значений поля в некоторой области (нелокальность). В среде с пространственной и временной дисперсией материальные уравнения (Di = D0i + ikEk, ik = ik + 4ik, (1.37)) и ji = ikEk, (1.40) имеют операторный вид: Это – наиболее общая форма линейных материальных уравнений, учитывающая нелокальность, запаздывание и анизотропию. Для однородной и стационарной среды материальные характеристики , и , называемые функциями отклика, должны зависеть только от разностей координат и времени R = r – r1, = t – t1: , (3.1) |
6.2 оказывается нелокальной. Эффекты нелокальности тем больше, чем сильнее поле меняется в пространстве. Обычно с ростом R функции отклика ij(R), ij(R), ij(R) достаточно быстро убывают. Представляет интерес вопрос о совместимости дисперсионных явлений с макроскопическим описанием полей в среде. Наиболее быстрый механизм установления поляризации – электронный. Его время релаксации порядка a/v, где а – характерный размер атома (элементарной ячейки), v – характерная скорость движения электронов в атоме. Дисперсия становится существенной, если >> v/a, где – частота изменения электромагнитного поля. Условие же применения макроскопического описания требует, чтобы длина, на которой существенно меняется напряженность макроскопического поля, значительно превосходила атомные размеры, то есть с/ >> a, или << c/a. Так как c/v ~ 137, существует область частот, при которой оба условия выполняются Рассмотрим однородную непироэлектрическую и неферромагнитную среду без пространственной дисперсии и без сторонних токов и зарядов. Материальные уравнения (3.1) – (3.3) с учетом формулы (1.38) в этом случае принимают вид свертки , (3.5) , 3.6) |
6.3 Взяв от соотношений (3.5) – (3.7) преобразование Фурье, получим: Di() = ij()Ej(), (3.9) Вi() = ij()Hj(), (3.10) ji() = ij()Ej(). (3.11) Здесь , , , , – спектры компонент электрического и магнитного поля, электрической и магнитной индукции и плотности тока соответственно, , (3.12) , (3.13) – (3.14) |
|
|
7.1 Соотношения Крамерса – Кронига Рассмотрим аналитические свойства комплексной проницаемости. Для простоты будем считать ее скалярной величиной, тензорный характер проницаемости не меняет полученных выводов. Интегрируя по частям соотношение (2.11) и учитывая, что функция (t) гладкая и () = 0, получим: . (3.17) Следовательно, , . Аналогично показывается, что , , , , , Если рассматривать частоту как комплексную переменную = + i, то функция () не имеет полюсов в верхней полуплоскости, то есть является там аналитической. Действительно, экспонента в подынтегральном выражении формулы , |
7.2. Устремим R к бесконечности, а r к нулю. В силу соотношения (3.17) подынтегральное выражение стремится к нулю при || быстрее, чем 1/||, следовательно интеграл по дуге С1 при R стремится к нулю. При r 0 интегралы по дугам С2 и С3 соответственно равны –42(0)/ и –. Поэтому . Разделяя в этой формуле действительную и мнимую части, получим соотношения Крамерса – Кронига: , (3.19) . (3.20) Несобственные интегралы в этих формулах понимаются в смысле главного значения. Аналогично легко получить соотношения Крамерса – Кронига для диэлектрической |
|
|
|
8.1 Пространственная дисперсия Связь между воздействием поля и реакцией среды является интегральной не только во времени. Эта связь может быть нелокальной, то есть реакция системы в некоторой точке может определяться не только напряженностью полей в той же самой точке, но и полями в окрестности данной точки. Этому соответствует интеграл по объему в соотношениях , (3.1) , (3.2) . (3.3). Аналогичное соотношение в случае немагнитной среды можно записать и для обобщенной индукции , (3.4) .(3.26) |
8.2ний (3.30) видно, что , , Если связь реакции среды с электрическим полем носит локальный характер, то есть ij(r, t) = (r)ij(t),ij(r, t) = (r)ij(t) то из уравнений (3.26) и (3.30) следует , , , то есть функция не зависит от волнового вектора k и может быть записана только как функция одной частоты , положим в ней k = 0. Таким образом, переход к пределу соответствует пренебрежению пространственной дисперсией, а учет зависимости комплексной проницаемости от волнового вектора – нелокальному характеру связи между полем и реакцией среды. |
8.3(3.32) Здесь обозначено E|| = k(kE)/k2 – проекция вектора Е на направление волнового вектора k, соответственно Е = Е – E|| – составляющая вектора Е, перпендикулярная волновому вектору k. Нетрудно видеть, что, если вектор напряженности электрического поля перпендикулярен волновому вектору, то kE = 0 и . Таким образом, величина является комплексной проницаемостью для поперечного (по отношению к волновому вектору) поля. Если же вектора k и E параллельны, то , то есть величина является комплексной проницаемостью для продольного поля. Рассмотрим изотропную магнитную среду. Из Фурье-образов уравнения (1.30) с материальным уравнением (Di = D0i + ikEk, ik = ik + 4ik, (1.37)) ikE = 4( + e), уравнения непрерывности . (1.21) и закона Ома ji = ikEk, (1.40) i = kЕ получим ikE( + 4i/) = 4e. Сравнивая это уравнение и Фурье-образ уравнения , (1.43) с материальным уравнением (3.32) ikE = 4e, получаем . |
8.4 получившееся двойное векторное произведение и сокращая на Е обе части, получим . (3.34) Соответственно, из уравнений (3.31) и (3.34) получаем: . (3.35) Таким образом, в изотропной среде различие продольной и поперечной проницаемостей обусловлено магнитными свойствами среды. Для переменных полей, пока вектор В можно выразить через вектор Е с помощью уравнения (1.28) как , для полного описания и электрических и магнитных свойств изотропной среды достаточно знать 2 зависящие от k и от величины – диэлектрические проницаемости и . Для постоянных полей при = 0 вектор В уже нельзя выразить через вектор Е, поэтому в этом случае для полного описания электромагнитных свойств однородных сред требуется знание и зависящей только от волнового вектора магнитной проницаемости (k). |
|
|
|
плотностью . С учетом тока связанных зарядов (1.18) полное выражение для средней плотности микроскопического тока может быть записано в виде: . (1.19) При вычислении средней плотности заряда следует учесть, что теперь электронейтральность тела в целом обеспечивается совокупностью как свободных, так и связанных зарядов, поэтому первое слагаемое в формуле (1.14) отлично от нуля, физически малый объем при этом может уже не быть электронейтральным: , (1.20) где – макроскопическая плотность свободных зарядов, непрерывная функция, в отличие от микроскопической плотности, определенной соотношением (1.6). Подставляя выражения (1.19) и (1.20) в усредненное микроскопическое уравнение непрерывности (1.13), получим макроскопическое уравнение непрерывности: . (1.21) Подставим теперь выражения (1.19) и (1.20) в усредненные уравнения (1.10) и (1.11): , (1.22) . (1.23) Вместо |
их ядра . (1.14) Для электронейтральной среды первое слагаемое в правой части уравнения (1.14) равно нулю, а второе может быть выражено через среднюю плотность дипольного момента Р: . (1.15) Поэтому средняя плотность связанных зарядов в этом случае равна: . (1.16) Аналогично можно усреднить и выражение (1.7) для плотности тока. Ограничиваясь первыми членами разложения по малому параметру |ri|/|r| с учетом формулы (1.15) получим: . Здесь a(b) |
Здесь k – номер атома (молекулы или элементарной ячейки) в твердом теле, Rk – радиус-вектор центра инерции k-го атома, ri(t) – радиус-вектор i-го заряда qi,k, отсчитанный от центра инерции соответствующего атома, r – радиус-вектор точки наблюдения. динамика носителей заряда может быть в классическом приближении описана уравнениями Ньютона: . (1.8) В правую часть уравнения (1.8) помимо сил Кулона и Лоренца входит также сила fe неэлектромагнитного происхождения. Такие силы могут быть связаны с гравитационным или ядерным взаимодействием, а также описывать столкновения частиц или квантовые переходы. В принципе, система уравнений (1.1) – (1.8) полностью описывает микроскопическое электромагнитное поле. Однако столь детальное описание поля в средах обычно невозможно из-за невообразимо большого количества уравнений вида (1.8), поскольку такое уравнение нужно написать для каждой частицы вещества. При этом все микроскопические величины e, b, и v, входящие в уравнения (1.1) – (1.8), испытывают существенные изменения на расстояниях порядка атомных размеров.Введем обозначения E = e – напряженность электрического поля в среде, В = b – магнитная индукция в среде. , где V(r) – физически малый объем с центром в точке r. Уравнения для введенных макроскопических (усредненных) полей получаются усреднением линейных микроскопических уравнений (1.1) – (1.5): , (1.9) |
|
|
|
(1.35) Соотношения (1.32) – (1.35) должны выполняться в любой точке границы раздела сред и представляют собой граничные условия, с помощью которых должно производиться сшивание решений системы уравнений Максвелла, получаемых в соприкасающихся между собой различных средах. Отметим, что поверхностная плотность заряда и плотность поверхностного тока i в соотношениях (1.33) и (1.34) соответственно учитывают как наведенные, так и сторонние заряды. Материальные уравнения Система уравнений Максвелла (1.28) – (1.31), . (1.31) даже дополненная граничными условиями (1.32) – (1.35) не является замкнутой. Ее нужно дополнить еще и материальными уравнениями, связывающими между собой величины H, D, j c B и Е. При этом удобно установить связь напряженностей поля с векторами Р и М, так как эти величины имеют более наглядный физический смысл.: Pi = P0i + ikEk, Mi = M0i + ikHk, i, k = x, y, z, (1.36) где ik и ik – тензоры диэлектрической и магнитной восприимчивости. Здесь и далее по повторяющимся тензорным индексам подразумевается суммирование по всему диапазону принимаемых значений (правило Эйнштейна). Физически более последовательно было бы записать в формуле (1.36) выражение для М в виде линейной функции В. Однако, в рассматриваемом приближении, как будет показано ниже, В и Н линейно связаны друг с другом, поэтому обе формы записи математически эквивалентны. |
например, уравнение (1.28) , по поверхности бесконечно малого прямоугольника, расположенного в переходном слое так, что его высота есть толщина переходного слоя, а основание длиной l параллельно возникающей при 0 границе раздела. Поскольку размеры прямоугольника малы в сравнении с расстояниями, на которых изменяются макроскопические поля, при вычислении интегралов можно воспользоваться теоремой о среднем. С учетом теоремы Стокса, устремляя высоту прямоугольника к нулю, то есть, переходя к резкой границе, получим: . Здесь Bn1 – проекция вектора магнитной индукции на границе раздела на нормаль n1 к поверхности прямоугольника, Е1 и Е2 – проекции векторе Е в первой и второй средах соответственно на направление, касательное к границе раздела. В силу произвольной ориентации вектора l, касательного границе раздела, получим: Е1 = Е2. (1.32) Аналогично, интегрируя уравнение (1.31) по той же области, получаем: |
|
|
|
Величина w, определенная соотношением (1.47), является суммой плотностей энергии электрического и магнитного полей, то есть плотностью энергии электромагнитного поля. Выражение (j + je)E в правой части уравнения (1.46) представляет собой работу, совершаемую в единицу времени силами электромагнитного поля над зарядами в единице объема. Соответственно вектор S, определенный соотношением (1.48) и называемый обычно вектором Пойнтинга, является плотностью потока энергии электромагнитного поля, а само уравнение (1.46) называется теоремой Пойнтинга. Отметим, что при выводе теоремы Пойнтинга (1.46) использовалось предположение о постоянстве диэлектрической и магнитной проницаемостей среды (1.37), справедливое лишь для медленно меняющихся полей и в линейном приближении. В высокочастотном случае и для нелинейных сред выражение для энергии электромагнитного поля и потока энергии более сложное. |
(1.30) и (1.31) принимают вид: , (1.43) . (1.44) Совместно с уравнениями (1.28) и (1.29) уравнения (1.43) и (1.44) образуют систему уравнений Максвелла, описывающую электромагнитное поле в веществе с помощью трех векторов Е, В и D, а не четырех, как система (1.28) – (1.31) , (1.29) . (1.31) . Свойства среды при этом описываются одним материальным уравнением (1.42) с помощью одной обобщенной проницаемости . D2n – D1n = 4ps. Граничное же условие (1.33) в случае ограниченности производной по времени напряженности электрического поля Е принимает вид [n´(В1 – В2)] = 4pi/c, (1.45) где i – полная плотность поверхностного тока. Энергия электромагнитного поля Одним из важнейших следствий уравнений макроскопической электродинамики является закон, связывающий плотность энергии и плотность потока энергии электромагнитного поля в макроскопических телах. пусть (Di = D0i + ikEk, ik = ik + 4ik, (1.37)) D0 = 0 и В0 = 0. |
|
|
. Важнейшим примером 4-тензора является тензор электромагнитного поля , (2.10) объединяющий компоненты векторов электрического е и магнитного b полей. То, что величины Fik образуют 4-тензор, следует из закона преобразования электромагнитного поля (2.5) и определения тензора. Антисимметричный тензор поля Fik можно представить в виде: , (2.11) где 4-вектор Ai имеет вид {А, i}, А – вектор-потенциал магнитного поля, – электростатический потенциал. Сравнивая выражение (2.11) с определением (2.10), получим . (2.12) Важность 4-векторов и 4-тензоров для современной физики заключается в том, что все физические величины могут быть объединены в такие совокупности, которые являются либо 4-векторами, либо 4-тензорами, либо 4-спинорами, квадратичные комбинации которых образуют 4-векторы. При этом иных математических объектов, связанных с преобразованием Лоренца, |
. (1.8) в виде: . Можно показать, что закону (2.5) сложения скоростей удовлетворяет правило Лоренца преобразования электрического и магнитного полей: , (2.6) где e|| и b|| – составляющие полей вдоль скорости V. Из формул (2.6) следует, что составляющие полей, параллельные вектору относительной скорости систем при переходе из одной системы в другую не меняются. При V << c из преобразования Лоренца (2.6) получаются известные в электричестве соотношения: . (2.7) Тензор электромагнитного поля Поскольку электрическое и магнитное поля при переходе к движущейся системе отсчета связаны преобразованиями Лоренца (2.6), описывающими поворот системы координат в 4-х мерном пространстве Минковского, их удобно описывать, пользуясь понятиями 4-х мерных векторов и тензоров, сокращенно 4-векторов и 4-тензоров. 4-вектором называется упорядоченная совокупность четырех величин {А1, А2, А3, А4}, которая при повороте системы координат преобразуется так же как координаты мировой точки {х1, х2, х3, х4} в псевдоевклидовом |
r = r – Vt, t = t, относительно которого инвариантны уравнения Ньютона, этому условию не удовлетворяют. Рассмотрим плоскую электромагнитную волну, распространяющуюся в вакууме вдоль оси х, тогда уравнение (2.1), описывающее ее, становится скалярным и одномерным: . (2.2) Пусть сопутствующая система координат (x, y, z) движется относительно лабораторной системы (x, y, z) с постоянной скоростью V, направленной вдоль оси х, причем соответствующие оси остаются параллельными. В сопутствующей системе волна должна остаться плоской и направленной вдоль оси х. Нетрудно показать, что при этом y = y, z = z. Рассмотрим линейное преобразование координаты х вида: x = ax + bt, t= fx + ht, где в силу однородности пространства и времени коэффициенты a, b, f и h не зависят ни от координат, ни от времени, но могут зависеть от скорости V сопутствующей (штрихованной) системы относительной лабораторной. Легко видеть, что V = b/a. Учитывая, что , требование неизменности, или инвариантности, волнового уравнения приводит к условиям: –f 2 + h2/c2 = 1, a2 – b2/c2 = 1, af – bh/c2 = 0, откуда следует преобразование Лоренца, связывающие координаты инерциальных систем |