Уравнения Максвелла

Уравне́ния Ма́ксвелла — система уравнений в дифференциальной или интегральной форме, описывающих электромагнитное поле и его связь с электрическими зарядами и токами в вакууме и сплошных средах.

Вместе с выражением для силы Лоренца, задающим меру воздействия электромагнитного поля на заряженные частицы, эти уравнения образуют полную систему уравнений классической электродинамики, называемую иногда уравнениями Максвелла — Лоренца.

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

История

Джеймс Клерк Максвелл

Уравнения, сформулированные Джеймсом Клерком Максвеллом, возникли на основе ряда важных экспериментальных открытий, которые были сделаны в начале XIX века. В 1820 году Ханс Кристиан Эрстед обнаружил, что пропускаемый через провод гальванический ток заставляет отклоняться магнитную стрелку компаса. Это открытие привлекло широкое внимание учёных того времени. В том же 1820 году Био и Савар экспериментально нашли выражение для порождаемой током магнитной индукции (закон Био — Савара), а Андре Мари Ампер обнаружил также, что взаимодействие на расстоянии возникает между двумя проводниками, по которым пропускается ток. Ампер ввёл термин «электродинамический» и выдвинул гипотезу, что природный магнетизм связан с существованием в магните круговых токов.

Влияние тока на магнит, обнаруженное Эрстедом, привело Майкла Фарадея к идее о том, что должно существовать обратное влияние магнита на токи. После длительных экспериментов, в 1831 году, Фарадей открыл, что перемещающийся возле проводника магнит порождает в проводнике электрический ток. Это явление было названо электромагнитной индукцией. Фарадей ввёл понятие «поля сил» — некоторой среды, находящейся между зарядами и токами. Его рассуждения носили качественный характер, однако они оказали огромное влияние на исследования Максвелла.

После открытий Фарадея стало ясно, что старые модели электромагнетизма (Ампер, Пуассон и др.) неполны. Вскоре появилась теория Вебера, основанная на дальнодействии. Однако к этому моменту вся физика, кроме теории тяготения, имела дело только с близкодействием (оптика, термодинамика, механика сплошных сред и др.). Гаусс, Риман и ряд других учёных высказывали догадки, что свет имеет электромагнитную природу, так что теория электромагнитных явлений тоже должна быть основана на близкодействии. Этот принцип стал существенной особенностью теории Максвелла.

В своём знаменитом «Трактате об электричестве и магнетизме» (1873) Максвелл писал:

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

Заменяя фарадеевский термин «поле сил» на понятие «напряжённость поля», Максвелл сделал его ключевым объектом своей теории:

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

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

Электрический ток создаёт магнитную индукцию (закон Ампера)

Анализируя известные эксперименты, Максвелл получил систему уравнений для электрического и магнитного полей. В 1855 году в своей самой первой статье «О фарадеевых силовых линиях» («On Faraday’s Lines of Force») он впервые записал в дифференциальной форме систему уравнений электродинамики, но не вводя ещё ток смещения. Такая система уравнений описывала все известные к тому времени экспериментальные данные, но не позволяла связать между собой заряды и токи и предсказать электромагнитные волны. Впервые ток смещения был введён Максвеллом в работе «О физических силовых линиях» («On Physical Lines of Force»), состоящей из четырёх частей и опубликованной в 1861—1862 годах. Обобщая закон Ампера, Максвелл вводит ток смещения, вероятно, чтобы связать токи и заряды уравнением непрерывности, которое уже было известно для других физических величин. Следовательно, в этой статье фактически была завершена формулировка полной системы уравнений электродинамики. В статье 1864 года «Динамическая теория электромагнитного поля» («A dynamical theory of the electromagnetic field») рассмотрена сформулированная ранее система уравнений из 20 уравнений для 20 неизвестных. В этой статье Максвелл впервые сформулировал понятие электромагнитного поля как физической реальности, имеющей собственную энергию и конечное время распространения, определяющее запаздывающий характер электромагнитного взаимодействия.

Переменный поток магнитного поля создаёт электрическое поле (закон Фарадея)

Оказалось, что не только ток, но и изменяющееся со временем электрическое поле (ток смещения) порождает магнитное поле. В свою очередь, в силу закона Фарадея, изменяющееся магнитное поле снова порождает электрическое. В результате в пустом пространстве может распространяться электромагнитная волна. Из уравнений Максвелла следовало, что её скорость равна скорости света, поэтому Максвелл сделал вывод об электромагнитной природе света.

Часть физиков выступила против теории Максвелла (особенно много возражений вызвала концепция тока смещения). Гельмгольц предложил свою теорию, компромиссную по отношению к моделям Вебера и Максвелла, и поручил своему ученику Генриху Герцу провести её экспериментальную проверку. Однако опыты Герца однозначно подтвердили правоту Максвелла.

Максвелл не использовал векторных обозначений и записывал свои уравнения в достаточно громоздком компонентном виде. В своём трактате он, кроме того, частично использовал кватернионную формулировку. Современная форма уравнений Максвелла появилась около 1884 года после работ Хевисайда, Герца и Гиббса. Они не только переписали систему Максвелла в векторном виде, но и симметризовали её, переформулировав в терминах поля, избавившись от электрического и магнитного потенциалов, игравших в теории Максвелла существенную роль, поскольку полагали, что эти функции являются лишь ненужными вспомогательными математическими абстракциями. Интересно, что современная физика поддерживает Максвелла, но не разделяет негативное отношение его ранних последователей к потенциалам. Электромагнитный потенциал играет важную роль в квантовой физике и проявляется как физически измеряемая величина в некоторых экспериментах, например, в эффекте Ааронова — Бома.

Система уравнений в формулировке Герца и Хевисайда некоторое время называлась уравнениями Герца — Хевисайда. Эйнштейн в классической статье «К электродинамике движущихся тел» назвал их уравнениями Максвелла — Герца. Иногда в литературе встречается также название уравнения Максвелла — Хевисайда.

Уравнения Максвелла сыграли важную роль при возникновении специальной теории относительности (СТО). Джозеф Лармор (1900 год) и независимо от него Хенрик Лоренц (1904 год) нашли преобразования координат, времени и электромагнитных полей, которые оставляют уравнения Максвелла инвариантными при переходе от одной инерциальной системы отсчёта к другой. Эти преобразования отличались от преобразований Галилея классической механики и, с подачи Анри Пуанкаре, стали называться преобразованиями Лоренца. Они стали математическим фундаментом специальной теории относительности.

Распространение электромагнитных волн со скоростью света первоначально интерпретировалось как возмущения некоторой среды, так называемого эфира. Были предприняты многочисленные попытки (см.исторический обзор) обнаружить движение Земли относительно эфира, однако они неизменно давали отрицательный результат. Поэтому Анри Пуанкаре высказал гипотезу о принципиальной невозможности обнаружить подобное движение (принцип относительности). Ему же принадлежит постулат о независимости скорости света от скорости его источника и вывод (вместе с Лоренцем), исходя из сформулированного так принципа относительности, точного вида преобразований Лоренца (при этом были показаны и групповые свойства этих преобразований). Эти две гипотезы (постулата) легли и в основу статьи Альберта Эйнштейна (1905 год). С их помощью он также вывел преобразования Лоренца и утвердил их общефизический смысл, особо подчеркнув возможность их применения для перехода из любой инерциальной системы отсчёта в любую другую инерциальную. Эта работа фактически ознаменовала собой построение специальной теории относительности. В СТО преобразования Лоренца отражают общие свойства пространства и времени, а модель эфира оказывается ненужной. Электромагнитные поля являются самостоятельными объектами, существующими наравне с материальными частицами.

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

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

На современном уровне понимания все фундаментальные частицы являются квантовыми возбуждениями («квантами») различных полей. Например, фотон — это квант электромагнитного поля, а электрон — квант спинорного поля. Поэтому полевой подход, предложенный Фарадеем и существенно развитый Максвеллом, является основой современной физики фундаментальных частиц, в том числе её стандартной модели.

Исторически несколько раньше он сыграл важную роль в появлении квантовой механики в формулировке Шрёдингера и вообще открытии квантовых уравнений, описывающих движение частиц, в том числе и релятивистских (уравнение Клейна — Гордона, уравнение Дирака), хотя первоначально аналогия с уравнениями Максвелла здесь виделась скорее лишь в общей идее, тогда как впоследствии оказалось, что она может быть понята как более конкретная и детальная (как это описано выше).

Также полевой подход, в целом восходящий к Фарадею и Максвеллу, стал центральным в теории гравитации (включая ОТО).

Запись уравнений Максвелла и системы единиц

Запись большинства уравнений в физике не зависит от выбора системы единиц. Однако в электродинамике это не так. В зависимости от выбора системы единиц в уравнениях Максвелла возникают различные коэффициенты (константы). Международная система единиц (СИ) является стандартом в технике и преподавании, однако споры среди физиков о её достоинствах и недостатках по сравнению с конкурирующей системой единиц СГС не утихают; здесь и всюду далее под СГС подразумевается исключительно симметричная гауссова система СГС. Преимущество системы СГС в электродинамике состоит в том, что все поля в ней имеют одну размерность, а уравнения, по мнению многих учёных, записываются проще и естественней. Поэтому СГС продолжает применяться в научных публикациях по электродинамике и в преподавании теоретической физики, например, в курсе теоретической физики Ландау и Лифшица. Однако для практических применений вводимые в СГС единицы измерений, многие из которых неименованы и неоднозначны, часто неудобны. Система СИ стандартизована и лучше самосогласованна, на этой системе построена вся современная метрология. Кроме того, система СИ обычно используется в курсах общей физики. В связи с этим все соотношения, если они по-разному записываются в системах СИ и СГС, далее приводятся в двух вариантах.

Иногда (например, в некоторых разделах «Фейнмановских лекций по физике», а также в современной квантовой теории поля) применяется система единиц, в которой скорость света, электрическая и магнитная постоянные принимаются за единицу: [math]\displaystyle{ c = \varepsilon_0 = \mu_0 = 1 }[/math]. В такой системе уравнения Максвелла записываются вообще без коэффициентов, все поля имеют единую размерность, а все потенциалы — свою единую. Такая система особенно удобна в ковариантной четырёхмерной формулировке законов электродинамики через 4-потенциал и 4-тензор электромагнитного поля.

Уравнения Максвелла в дифференциальной форме

Уравнения Максвелла представляют собой в векторной записи систему из четырёх уравнений, сводящуюся в компонентном представлении к восьми (два векторных уравнения содержат по три компоненты каждое плюс два скалярных) линейным дифференциальным уравнениям в частных производных первого порядка для 12 компонент четырёх векторных и псевдовекторных функций ([math]\displaystyle{ \mathbf{D},\;\mathbf{E},\;\mathbf{H},\;\mathbf{B} }[/math]):

Название	СГС	СИ	Примерное словесное выражение
Закон Гаусса	[math]\displaystyle{ \nabla\cdot\mathbf{D}=4\pi \rho }[/math]	[math]\displaystyle{ \nabla\cdot\mathbf{D}= \rho }[/math]	Электрический заряд является источником электрической индукции.
Закон Гаусса для магнитного поля	[math]\displaystyle{ \nabla\cdot\mathbf{B}=0 }[/math]	[math]\displaystyle{ \nabla\cdot\mathbf{B}=0 }[/math]	Магнитные заряды не обнаружены.
Закон индукции Фарадея	[math]\displaystyle{ \nabla\times\mathbf{E}=-\frac{1}{c}\,\frac{\partial\mathbf{B}}{\partial t} }[/math]	[math]\displaystyle{ \nabla\times\mathbf{E}=-\frac{\partial\mathbf{B}}{\partial t} }[/math]	Изменение магнитной индукции порождает вихревое электрическое поле.
Теорема о циркуляции магнитного поля	[math]\displaystyle{ \nabla\times\mathbf{H}=\frac{4\pi}{c} \mathbf{j}+\frac{1}{c}\frac{\partial\mathbf{D}}{\partial t} }[/math]	[math]\displaystyle{ \nabla\times\mathbf{H}= \mathbf{j}+\frac{\partial\mathbf{D}}{\partial t} }[/math]	Электрический ток и изменение электрической индукции порождают вихревое магнитное поле

Жирным шрифтом в дальнейшем обозначаются векторные и псевдовекторные величины, курсивом — скалярные.

Введённые обозначения:

• [math]\displaystyle{ \rho\ }[/math] — объёмная плотность стороннего электрического заряда (в единицах СИ — Кл/м³);
• [math]\displaystyle{ \mathbf{j} }[/math] — плотность электрического тока (плотность тока проводимости) (в единицах СИ — А/м²); в простейшем случае — случае тока, порождаемого одним типом носителей заряда, она выражается просто как [math]\displaystyle{ \mathbf{j}=\mathbf{u}\rho_1 }[/math], где [math]\displaystyle{ \mathbf{u} }[/math] — (средняя) скорость движения этих носителей в окрестности данной точки, [math]\displaystyle{ \rho_1 }[/math] — плотность заряда этого типа носителей (она в общем случае не совпадает с [math]\displaystyle{ \rho }[/math]); в общем случае это выражение надо усреднить по разным типам носителей;
• [math]\displaystyle{ c }[/math] — скорость света в вакууме (299 792 458 м/с);
• [math]\displaystyle{ \mathbf E }[/math] — напряжённость электрического поля (в единицах СИ — В/м);
• [math]\displaystyle{ \mathbf H }[/math] — напряжённость магнитного поля (в единицах СИ — А/м);
• [math]\displaystyle{ \mathbf D }[/math] — электрическая индукция (в единицах СИ — Кл/м²);
• [math]\displaystyle{ \mathbf B }[/math] — магнитная индукция (в единицах СИ — Тл = Вб/м² = кг•с⁻²•А⁻¹);
• [math]\displaystyle{ \nabla }[/math] — дифференциальный оператор набла, при этом:

  [math]\displaystyle{ \nabla\times \mathbf{E} \equiv \mathrm{rot}\,\mathbf{E} }[/math] означает ротор вектора [math]\displaystyle{ \mathbf E }[/math],

  [math]\displaystyle{ \nabla\cdot\mathbf{E}\equiv\mathrm{div}\,\mathbf{E} }[/math] означает дивергенцию вектора [math]\displaystyle{ \mathbf E }[/math].

Приведённые выше уравнения Максвелла не составляют ещё полной системы уравнений электромагнитного поля, поскольку они не содержат свойств среды, в которой возбуждено электромагнитное поле. Соотношения, связывающие величины [math]\displaystyle{ \mathbf{E} }[/math], [math]\displaystyle{ \mathbf{B} }[/math], [math]\displaystyle{ \mathbf{D} }[/math], [math]\displaystyle{ \mathbf{H} }[/math] и [math]\displaystyle{ \mathbf{j} }[/math] и учитывающие индивидуальные свойства среды, называются материальными уравнениями.

Уравнения Максвелла в интегральной форме

При помощи формулы Остроградского — Гаусса и теоремы Стокса дифференциальным уравнениям Максвелла можно придать форму интегральных уравнений:

Название	СГС	СИ	Примерное словесное выражение
Закон Гаусса	[math]\displaystyle{ \oint\limits_{\mathbf{s}}\mathbf{D}\cdot d\mathbf{s}=4\pi Q }[/math]	[math]\displaystyle{ \oint\limits_\mathbf{s}\mathbf{D}\cdot d\mathbf{s}= Q }[/math]	Поток электрической индукции через замкнутую поверхность пропорционален величине свободного заряда, находящегося в объёме, ограниченном этой поверхностью.
Закон Гаусса для магнитного поля	[math]\displaystyle{ \oint\limits_\mathbf{s}\mathbf{B}\cdot d\mathbf{s}=0 }[/math]	[math]\displaystyle{ \oint\limits_\mathbf{s}\mathbf{B}\cdot d\mathbf{s}=0 }[/math]	Поток магнитной индукции через замкнутую поверхность равен нулю (магнитные заряды не обнаружены).
Закон индукции Фарадея	[math]\displaystyle{ \oint\limits_\mathbf{l}\mathbf{E}\cdot d\mathbf{l}= }[/math] [math]\displaystyle{ -\frac{1}{c}\frac{d}{d t}\int\limits_\mathbf{s} \mathbf{B}\cdot d\mathbf{s} }[/math]	[math]\displaystyle{ \oint\limits_\mathbf{l}\mathbf{E}\cdot d\mathbf{l}= }[/math] [math]\displaystyle{ -\frac{d}{d t}\int\limits_\mathbf{s}\mathbf{B}\cdot d\mathbf{s} }[/math]	Изменение потока магнитной индукции, проходящего через незамкнутую поверхность [math]\displaystyle{ s }[/math], взятое с обратным знаком, пропорционально циркуляции электрического поля на замкнутом контуре [math]\displaystyle{ l }[/math], который является границей поверхности [math]\displaystyle{ s }[/math].
Теорема о циркуляции магнитного поля	[math]\displaystyle{ \oint\limits_\mathbf{l}\mathbf{H}\cdot d\mathbf{l}= }[/math] [math]\displaystyle{ \frac{4\pi}{c} I+\frac{1}{c}\frac{d}{d t}\int\limits_\mathbf{s}\mathbf{D}\cdot d\mathbf{s} }[/math]	[math]\displaystyle{ \oint\limits_\mathbf{l}\mathbf{H}\cdot d\mathbf{l}= }[/math] [math]\displaystyle{ I+\frac{d}{d t}\int\limits_\mathbf{s}\mathbf{D}\cdot d\mathbf{s} }[/math]	Полный электрический ток свободных зарядов и изменение потока электрической индукции через незамкнутую поверхность [math]\displaystyle{ s }[/math] пропорциональны циркуляции магнитного поля на замкнутом контуре [math]\displaystyle{ l }[/math], который является границей поверхности [math]\displaystyle{ s }[/math].

Поток электрического поля через замкнутую поверхность

Введённые обозначения:

• [math]\displaystyle{ s\ }[/math] — двумерная замкнутая в случае теоремы Гаусса поверхность, ограничивающая объём [math]\displaystyle{ v\ }[/math], и открытая поверхность в случае законов Фарадея и Ампера — Максвелла (её границей является замкнутый контур [math]\displaystyle{ l\ }[/math]).
• [math]\displaystyle{ Q=\int\limits_v \rho\, dv\ }[/math] — электрический заряд, заключённый в объёме [math]\displaystyle{ v\ }[/math], ограниченном поверхностью [math]\displaystyle{ s\ }[/math] (в единицах СИ — Кл);
• [math]\displaystyle{ I=\int\limits_\mathbf{s} \mathbf{j}\cdot d\mathbf{s}\ }[/math] — электрический ток, проходящий через поверхность [math]\displaystyle{ s\ }[/math] (в единицах СИ — А).

При интегрировании по замкнутой поверхности вектор элемента площади [math]\displaystyle{ d\mathbf{s} }[/math] направлен из объёма наружу. Ориентация [math]\displaystyle{ d\mathbf{s} }[/math] при интегрировании по незамкнутой поверхности определяется направлением правого винта, «вкручивающегося» при повороте в направлении обхода контурного интеграла по [math]\displaystyle{ d\mathbf{l} }[/math].

Словесное описание законов Максвелла, например, закона Фарадея, несёт отпечаток традиции, поскольку вначале при контролируемом изменении магнитного потока регистрировалось возникновение электрического поля (точнее электродвижущей силы). В общем случае в уравнениях Максвелла (как в дифференциальной, так и в интегральной форме) векторные функции [math]\displaystyle{ \mathbf{E}, \mathbf{B}, \mathbf{D}, \mathbf{H} }[/math] являются равноправными неизвестными величинами, определяемыми в результате решения уравнений.

Сила Лоренца

Основная статья: Сила Лоренца

При решении уравнений Максвелла распределения зарядов [math]\displaystyle{ \rho }[/math] и токов [math]\displaystyle{ \mathbf{j} }[/math] часто считаются заданными. С учётом граничных условий и материальных уравнений это позволяет определить напряжённость электрического поля [math]\displaystyle{ \mathbf{E} }[/math] и магнитную индукцию [math]\displaystyle{ \mathbf{B} }[/math], которые, в свою очередь, определяют силу, действующую на пробный заряд [math]\displaystyle{ q }[/math], движущийся со скоростью [math]\displaystyle{ \mathbf{u} }[/math]. Эта сила называется силой Лоренца:

СГС	СИ
[math]\displaystyle{ \mathbf{F}=q\,\mathbf{E}+\frac{q}{c}\,[\mathbf{u}\times\mathbf{B}] }[/math]	[math]\displaystyle{ \mathbf{F}=q\,\mathbf{E}+q\,[\mathbf{u}\times\mathbf{B}] }[/math]

Электрическая составляющая силы направлена параллельно электрическому полю, а магнитная — перпендикулярна скорости заряда и магнитной индукции. Впервые выражение для силы, действующей на заряд в магнитном поле (электрическая компонента была известна), получил в 1889 году Хевисайд за три года до Хендрика Лоренца, который вывел выражение для этой силы в 1892 году.

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

Размерные константы в уравнениях Максвелла

В гауссовой системе единиц СГС все поля имеют одинаковую размерность, и в уравнениях Максвелла фигурирует единственная фундаментальная константа [math]\displaystyle{ c }[/math], имеющая размерность скорости, которая сейчас называется скоростью света (именно равенство этой константы скорости распространения света дало Максвеллу основания для гипотезы об электромагнитной природе света).

В системе единиц СИ, чтобы связать электрическую индукцию и напряжённость электрического поля в вакууме, вводится электрическая постоянная [math]\displaystyle{ \varepsilon_0 }[/math] ([math]\displaystyle{ \mathbf{D}=\varepsilon_0\mathbf{E} }[/math]). Магнитная постоянная [math]\displaystyle{ \mu_0 }[/math] является таким же коэффициентом пропорциональности для магнитного поля в вакууме ([math]\displaystyle{ \mathbf{B}=\mu_0\mathbf{H} }[/math]). Названия электрическая постоянная и магнитная постоянная сейчас стандартизованы. Ранее для этих величин также использовались, соответственно, названия электрическая (диэлектрическая) и магнитная проницаемости вакуума.

Скорость электромагнитного излучения в вакууме (скорость света) в СИ появляется при выводе волнового уравнения:

[math]\displaystyle{ c=\frac{1}{\sqrt{\varepsilon_0\mu_0}}. }[/math]

В системе единиц СИ в качестве точной размерной константы определена скорость света в вакууме [math]\displaystyle{ c }[/math], а магнитная постоянная [math]\displaystyle{ \mu_0 }[/math] после изменения 2018—2019 годов является экспериментально определяемой величиной. Через них выражается электрическая постоянная [math]\displaystyle{ \varepsilon_0 }[/math].

Значения скорости света, электрической и магнитной постоянных приведены в таблице:

Символ	Наименование	Численное значение	Единицы измерения СИ
[math]\displaystyle{ c }[/math]	Постоянная скорости света	[math]\displaystyle{ 299\;792\;458 }[/math] (точно)	м/с
[math]\displaystyle{ \mu_0 }[/math]	Магнитная постоянная	[math]\displaystyle{ 1{,}256\;637\times 10^{-6} }[/math]	Гн/м
[math]\displaystyle{ \varepsilon_0=1/(\mu_0 c^2) }[/math]	Электрическая постоянная	[math]\displaystyle{ 8{,}854\;188\times 10^{-12} }[/math]	Ф/м

Иногда вводится величина, называемая «волновым сопротивлением вакуума», или «импедансом» вакуума:

[math]\displaystyle{ Z_0=\sqrt{\frac{\mu_0}{\varepsilon_0}}=\mu_0 c \approx 120\pi }[/math] Ом.

В системе СГС [math]\displaystyle{ Z_0=1 }[/math]. Эта величина имеет смысл отношения амплитуд напряжённостей электрического и магнитного полей плоской электромагнитной волны в вакууме. Однако приписать этой величине физический смысл волнового сопротивления нельзя, поскольку в той же системе СГС её размерность не совпадает с размерностью сопротивления.

Уравнения Максвелла в среде

Чтобы получить полную систему уравнений электродинамики, к системе уравнений Максвелла необходимо добавить материальные уравнения, связывающие величины [math]\displaystyle{ \mathbf{j} }[/math], [math]\displaystyle{ \mathbf{H} }[/math], [math]\displaystyle{ \mathbf{D} }[/math], [math]\displaystyle{ \mathbf{E} }[/math], [math]\displaystyle{ \mathbf{B} }[/math], в которых учтены индивидуальные свойства среды. Способ получения материальных уравнений дают молекулярные теории поляризации, намагниченности и электропроводности среды, использующие идеализированные модели среды. Применяя к ним уравнения классической или квантовой механики, а также методы статистической физики, можно установить связь между векторами [math]\displaystyle{ \mathbf{j} }[/math], [math]\displaystyle{ \mathbf{H} }[/math], [math]\displaystyle{ \mathbf{D} }[/math] с одной стороны и [math]\displaystyle{ \mathbf{E} }[/math], [math]\displaystyle{ \mathbf{B} }[/math] с другой стороны.

Связанные заряды и токи

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

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

Аналогичным образом происходит магнитная поляризация (намагничивание) в материалах, в которых составляющие их атомы и молекулы имеют магнитные моменты, связанные со спином и орбитальным моментом ядер и электронов. Угловые моменты атомов можно представить в виде циркулярных токов. На границе материала совокупность таких микроскопических токов эквивалентна макроскопическим токам, циркулирующим вдоль поверхности, несмотря на то, что движение зарядов в отдельных магнитных диполях происходит лишь в микромасштабе (связанные токи).

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

В среде сторонние электрические и магнитные поля вызывают поляризацию и намагничивание вещества, которые макроскопически описываются соответственно вектором поляризации [math]\displaystyle{ \mathbf P }[/math] и вектором намагниченности [math]\displaystyle{ \mathbf M }[/math] вещества, и вызваны появлением связанных зарядов [math]\displaystyle{ \rho_b\ }[/math] и токов [math]\displaystyle{ \mathbf{j}_b }[/math]. В результате поле в среде оказывается суммой внешних полей и полей, вызванных связанными зарядами и токами.

СГС	СИ
[math]\displaystyle{ \rho_b=-\nabla\cdot\mathbf{P} }[/math] [math]\displaystyle{ \mathbf{j}_b = c \nabla\times\mathbf{M} + \frac{\partial\mathbf{P}}{\partial t} }[/math]	[math]\displaystyle{ \rho_b = -\nabla\cdot\mathbf{P} }[/math] [math]\displaystyle{ \mathbf{j}_b = \nabla\times\mathbf{M} + \frac{\partial\mathbf{P}}{\partial t} }[/math]

Поляризация [math]\displaystyle{ \mathbf P }[/math] и намагниченность вещества [math]\displaystyle{ \mathbf M }[/math] связаны с векторами напряжённости и индукции электрического и магнитного поля следующими соотношениями:

СГС	СИ
[math]\displaystyle{ \mathbf{D}=\mathbf{E}+4\pi \mathbf{P} }[/math] [math]\displaystyle{ \mathbf{B}=\mathbf{H}+4\pi \mathbf{M} }[/math]	[math]\displaystyle{ \mathbf{D}=\varepsilon_0\mathbf{E}+\mathbf{P} }[/math] [math]\displaystyle{ \mathbf{B}=\mu_0(\mathbf{H}+\mathbf{M}) }[/math]

Поэтому, выражая векторы [math]\displaystyle{ \mathbf{D} }[/math] и [math]\displaystyle{ \mathbf{H} }[/math] через [math]\displaystyle{ \mathbf{E} }[/math], [math]\displaystyle{ \mathbf{B} }[/math], [math]\displaystyle{ \rho_b\ }[/math] и [math]\displaystyle{ \mathbf{j}_b }[/math], можно получить математически эквивалентную систему уравнений Максвелла:

СГС	СИ
[math]\displaystyle{ \nabla\cdot\mathbf{E}=4\pi[\rho_f+\rho_b] }[/math]	[math]\displaystyle{ \nabla\cdot\mathbf{E}= \frac{1}{\varepsilon_0}[\rho_f+\rho_b] }[/math]
[math]\displaystyle{ \nabla\cdot\mathbf{B}=0 }[/math]	[math]\displaystyle{ \nabla\cdot\mathbf{B}=0 }[/math]
[math]\displaystyle{ \nabla\times\mathbf{E}=-\frac{1}{c}\,\frac{\partial\mathbf{B}}{\partial t} }[/math]	[math]\displaystyle{ \nabla\times\mathbf{E}=-\frac{\partial\mathbf{B}}{\partial t} }[/math]
[math]\displaystyle{ \nabla\times\mathbf{B}=\frac{4\pi}{c} [\mathbf{j}_b+\mathbf{j}_f]+\frac{1}{c}\frac{\partial\mathbf{E}}{\partial t} }[/math]	[math]\displaystyle{ \nabla\times\mathbf{B}= \mu_0[\mathbf{j}_b+\mathbf{j}_f]+\frac{1}{c^2}\frac{\partial\mathbf{E}}{\partial t} }[/math]

Индексом [math]\displaystyle{ f\ }[/math] здесь обозначены свободные заряды и токи. Уравнения Максвелла в такой форме являются фундаментальными в том смысле, что они не зависят от модели электромагнитного устройства вещества. Разделение зарядов и токов на свободные и связанные позволяет «спрятать» в [math]\displaystyle{ \rho_b\ }[/math], [math]\displaystyle{ \mathbf{j}_b }[/math], а затем в [math]\displaystyle{ \mathbf{P}, \mathbf{M} }[/math] и, следовательно, в [math]\displaystyle{ \mathbf{D},\mathbf{H} }[/math] сложный микроскопический характер электромагнитного поля в среде.

Материальные уравнения

Материальные уравнения устанавливают связь между [math]\displaystyle{ \mathbf{D},\mathbf{H} }[/math] и [math]\displaystyle{ \mathbf{E},\mathbf{B} }[/math]. При этом учитываются индивидуальные свойства среды. На практике в материальных уравнениях обычно используются экспериментально определяемые коэффициенты (зависящие в общем случае от частоты электромагнитного поля), которые собраны в различных справочниках физических величин.

• В слабых электромагнитных полях, сравнительно медленно меняющихся в пространстве и во времени, в случае изотропных, неферромагнитных и несегнетоэлектрических сред справедливо приближение, в котором поляризуемость и намагниченность линейно зависят от приложенных полей:

СГС	СИ
[math]\displaystyle{ \mathbf{P}=\chi_e\mathbf{E} }[/math] [math]\displaystyle{ \mathbf{M}=\chi_m\mathbf{H} }[/math]	[math]\displaystyle{ \mathbf{P}=\varepsilon_0\chi_e\mathbf{E} }[/math] [math]\displaystyle{ \mathbf{M}=\chi_m\mathbf{H}, }[/math]

где введены безразмерные константы: [math]\displaystyle{ \chi_e\ }[/math] — диэлектрическая восприимчивость и [math]\displaystyle{ \chi_m\ }[/math] — магнитная восприимчивость вещества (в системе единиц СИ эти константы в [math]\displaystyle{ 4\pi\ }[/math] раз больше, чем в гауссовой системе СГС). Соответственно, материальные уравнения для электрической и магнитной индукций записываются в следующем виде:

СГС	СИ
[math]\displaystyle{ \mathbf{D}=\varepsilon\mathbf{E}=(1+4\pi\chi_e)\mathbf{E} }[/math] [math]\displaystyle{ \mathbf{B}=\mu\mathbf{H}=(1+4\pi\chi_m)\mathbf{H} }[/math]	[math]\displaystyle{ \mathbf{D}=\varepsilon_0\varepsilon\mathbf{E}=\varepsilon_0(1+\chi_e)\mathbf{E} }[/math] [math]\displaystyle{ \mathbf{B}=\mu_0\mu\mathbf{H}=\mu_0(1+\chi_m)\mathbf{H}, }[/math]

где [math]\displaystyle{ \varepsilon\ }[/math] — относительная диэлектрическая проницаемость, [math]\displaystyle{ \mu\ }[/math] — относительная магнитная проницаемость. Размерные величины [math]\displaystyle{ \varepsilon_0\varepsilon }[/math] (в единицах СИ — Ф/м) и [math]\displaystyle{ \mu_0\mu }[/math] (в единицах СИ — Гн/м), возникающие в системе СИ, называются абсолютная диэлектрическая проницаемость и абсолютная магнитная проницаемость соответственно.

• В проводниках существует связь между плотностью тока и напряжённостью электрического поля, в хорошем приближении выражаемая законом Ома:

[math]\displaystyle{ \mathbf{j}=\sigma\mathbf{E}, }[/math]

где [math]\displaystyle{ \sigma }[/math] — удельная проводимость среды (в единицах СИ — Ом⁻¹•м⁻¹).

• В анизотропной среде [math]\displaystyle{ \varepsilon }[/math], [math]\displaystyle{ \mu }[/math] и [math]\displaystyle{ \sigma }[/math] являются тензорами [math]\displaystyle{ \hat\varepsilon }[/math], [math]\displaystyle{ \hat\mu }[/math] и [math]\displaystyle{ \hat\sigma }[/math]. В системе координат главных осей они могут быть описаны диагональными матрицами. В этом случае связь между напряжённостями полей и индукциями имеют различные коэффициенты по каждой координате. Например, в системе СИ:

[math]\displaystyle{ \begin{array}{lll} D_x = \varepsilon_{0}\varepsilon_{xx}\,E_x,~~~~~ & D_y = \varepsilon_{0}\varepsilon_{yy}\,E_y,~~~~~ & D_z = \varepsilon_{0}\varepsilon_{zz}\,E_z, \\ [3mm] B_x = \mu_0\mu_{xx}\,H_x,~~~~~ & B_y = \mu_0\mu_{yy}\,H_y,~~~~~ & B_z = \mu_0\mu_{zz}\,H_z. \end{array} }[/math]

• Хотя для широкого класса веществ линейное приближение для слабых полей выполняется с хорошей точностью, в общем случае зависимость между [math]\displaystyle{ \mathbf{D},\mathbf{H} }[/math] и [math]\displaystyle{ \mathbf{E},\mathbf{B} }[/math] может быть нелинейной. В этом случае проницаемости среды не являются константами, а зависят от величины поля в данной точке. Кроме того, более сложная связь между [math]\displaystyle{ \mathbf{D},\mathbf{H} }[/math] и [math]\displaystyle{ \mathbf{E},\mathbf{B} }[/math] наблюдается в средах с пространственной или временной дисперсиями. В случае пространственной дисперсии токи и заряды в данной точке пространства зависят от величины поля не только в той же точке, но и в соседних точках. В случае временной дисперсии поляризация и намагниченность среды не определяются только величиной поля в данный момент времени, а зависят также от величины полей в предшествующие моменты времени. В самом общем случае нелинейных и неоднородных сред с дисперсией материальные уравнения в системе СИ принимают интегральный вид:

[math]\displaystyle{ \mathbf{P}(\mathbf{r}, t) = \varepsilon_0 \int\limits_v\!\int\limits_{-\infty}^t \hat{\chi}_e (\mathbf{r},\mathbf{r}' ,t-t', \mathbf{E}, \mathbf{H})\, \mathbf{E}(\mathbf{r}', t')\, d t' d^3\mathbf{r}'; }[/math] [math]\displaystyle{ \mathbf{M}(\mathbf{r}, t) = \int\limits_v\!\int\limits_{-\infty}^t \hat{\chi}_m (\mathbf{r},\mathbf{r}', t-t', \mathbf{E}, \mathbf{H})\, \mathbf{H}(\mathbf{r}', t')\, d t' d^3\mathbf{r}'. }[/math]

Аналогичные уравнения получаются в гауссовой системе СГС (если формально положить [math]\displaystyle{ \varepsilon_0=1 }[/math]).

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

В изотропных и однородных средах без дисперсии уравнения Максвелла принимают следующий вид:

СГС	СИ
[math]\displaystyle{ \nabla\cdot\mathbf{E}=4\pi\,\frac{\rho}{\varepsilon} }[/math] [math]\displaystyle{ \nabla\cdot\mathbf{B}=0 }[/math] [math]\displaystyle{ \nabla\times\mathbf{E}=-\frac{1}{c}\frac{\partial\mathbf{B}}{\partial t} }[/math] [math]\displaystyle{ \nabla\times\mathbf{B}=\frac{4\pi}{c}\,\mu\,\mathbf{j}+\frac{\varepsilon\mu}{c}\frac{\partial\mathbf{E}}{\partial t} }[/math]	[math]\displaystyle{ \nabla\cdot\mathbf{E}=\frac{\rho}{\varepsilon\varepsilon_0} }[/math] [math]\displaystyle{ \nabla\cdot\mathbf{B}=0 }[/math] [math]\displaystyle{ \nabla\times\mathbf{E}=-\frac{\partial\mathbf{B}}{\partial t} }[/math] [math]\displaystyle{ \nabla\times\mathbf{B}={\mu\mu_0}\mathbf{j}+\frac{\varepsilon\mu}{c^2}\frac{\partial\mathbf{E}}{\partial t} }[/math]

В оптическом диапазоне частот вместо диэлектрической проницаемости [math]\displaystyle{ \varepsilon }[/math] используется показатель преломления [math]\displaystyle{ n=\sqrt{\varepsilon\mu} }[/math], показывающий отличие скорости распространения монохроматической световой волны в среде от скорости света в вакууме. При этом в оптическом диапазоне диэлектрическая проницаемость обычно заметно меньше, чем на низких частотах, а магнитная проницаемость большинства оптических сред практически равна единице. Показатель преломления большинства прозрачных материалов составляет от 1 до 2, достигая 5 у некоторых полупроводников. В вакууме и диэлектрическая, и магнитная проницаемости равны единице: [math]\displaystyle{ \varepsilon=\mu=1 }[/math].

Поскольку уравнения Максвелла в линейной среде являются линейными относительно полей [math]\displaystyle{ (\mathbf{E},\mathbf{B}) }[/math] и свободных зарядов и токов [math]\displaystyle{ (\rho,\mathbf{j}) }[/math], справедлив принцип суперпозиции:

Если распределения зарядов и токов [math]\displaystyle{ (\rho_1,\mathbf{j}_1) }[/math] создают электромагнитное поле с компонентами [math]\displaystyle{ (\mathbf{E}_1,\mathbf{B}_1) }[/math], а другие распределения [math]\displaystyle{ (\rho_2,\mathbf{j}_2) }[/math] создают, соответственно, поле [math]\displaystyle{ (\mathbf{E}_2,\mathbf{B}_2) }[/math], то суммарное поле, создаваемое источниками [math]\displaystyle{ (\rho_1+\rho_2,\mathbf{j}_1+\mathbf{j}_2) }[/math], будет равно [math]\displaystyle{ (\mathbf{E}_1+\mathbf{E}_2,\mathbf{B}_1+\mathbf{B}_2) }[/math].

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

Граничные условия

Во многих случаях неоднородную среду можно представить в виде совокупности кусочно-непрерывных однородных областей, разделённых бесконечно тонкими границами. При этом можно решать уравнения Максвелла в каждой области, «сшивая» на границах получающиеся решения. В частности, при рассмотрении решения в конечном объёме необходимо учитывать условия на границах объёма с окружающим бесконечным пространством. Граничные условия получаются из уравнений Максвелла предельным переходом. Для этого проще всего воспользоваться уравнениями Максвелла в интегральной форме.

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

СГС	СИ
[math]\displaystyle{ (\mathbf{E}_1-\mathbf{E}_2)\times \mathbf{n}_{1-2} = \mathbf 0 }[/math], [math]\displaystyle{ (\mathbf{H}_1-\mathbf{H}_2)\times \mathbf{n}_{1-2} = \frac{4\pi}{c}\mathbf{j}_s }[/math],	[math]\displaystyle{ (\mathbf{E}_1-\mathbf{E}_2)\times \mathbf{n}_{1-2} = \mathbf 0 }[/math], [math]\displaystyle{ (\mathbf{H}_1-\mathbf{H}_2)\times \mathbf{n}_{1-2} = \mathbf{j}_s }[/math],

где [math]\displaystyle{ \mathbf{n}_{1-2} }[/math] — единичный вектор нормали к поверхности, направленный из среды 1 в среду 2 и имеющий размерность, обратную длине, [math]\displaystyle{ \mathbf{j}_s }[/math] — плотность поверхностных свободных токов вдоль границы (то есть не включая связанных токов намагничивания, складывающихся на границе среды из микроскопических молекулярных и т.п. токов). Первое граничное условие можно интерпретировать как непрерывность на границе областей тангенциальных компонент напряжённостей электрического поля (из второго следует, что тангенциальные компоненты напряжённости магнитного поля непрерывны только при отсутствии поверхностных токов на границе).

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

СГС	СИ
[math]\displaystyle{ (\mathbf{D}_1-\mathbf{D}_2)\cdot \mathbf{n}_{1-2} = -4\pi\rho_s }[/math], [math]\displaystyle{ (\mathbf{B}_1-\mathbf{B}_2)\cdot \mathbf{n}_{1-2} = 0 }[/math],	[math]\displaystyle{ (\mathbf{D}_1-\mathbf{D}_2)\cdot \mathbf{n}_{1-2} = -\rho_s }[/math], [math]\displaystyle{ (\mathbf{B}_1-\mathbf{B}_2)\cdot \mathbf{n}_{1-2} = 0 }[/math],

где [math]\displaystyle{ \rho_s\ }[/math] — поверхностная плотность свободных зарядов (то есть не включающая в себя связанных зарядов, возникающих на границе среды вследствие диэлектрической поляризации самой среды).

Эти граничные условия показывают непрерывность нормальной компоненты вектора магнитной индукции (нормальная компонента электрической индукции непрерывна только при отсутствии на границе поверхностных зарядов).

Из уравнения непрерывности можно получить граничное условие для токов:

[math]\displaystyle{ (\mathbf{j}_1-\mathbf{j}_2)\cdot \mathbf{n}_{1-2} = \frac{\partial}{\partial t}\rho_s }[/math],

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

СГС	СИ
[math]\displaystyle{ \mathbf{E}_1\times \mathbf{n}_{1-2} = \mathbf 0 }[/math], [math]\displaystyle{ \mathbf{H}_1\times \mathbf{n}_{1-2} = \frac{4\pi}{c}\mathbf{j}_s }[/math], [math]\displaystyle{ \mathbf{D}_1 \cdot \mathbf{n}_{1-2} = -4\pi\rho_s }[/math], [math]\displaystyle{ \mathbf{B}_1\cdot \mathbf{n}_{1-2} = 0 }[/math],	[math]\displaystyle{ \mathbf{E}_1\times \mathbf{n}_{1-2} = \mathbf 0 }[/math], [math]\displaystyle{ \mathbf{H}_1\times \mathbf{n}_{1-2} = \mathbf{j}_s }[/math], [math]\displaystyle{ \mathbf{D}_1\cdot \mathbf{n}_{1-2} = -\rho_s }[/math], [math]\displaystyle{ \mathbf{B}_1\cdot \mathbf{n}_{1-2} = 0 }[/math],

Законы сохранения

Уравнения Максвелла содержат в себе законы сохранения заряда и энергии электромагнитного поля.

Уравнение непрерывности

Источники полей ([math]\displaystyle{ \rho,~\mathbf{j} }[/math]) не могут быть заданы произвольным образом. Применяя операцию дивергенции к четвёртому уравнению (закон Ампера—Максвелла) и используя первое уравнение (закон Гаусса), можно получить уравнение непрерывности для зарядов и токов:

[math]\displaystyle{ \nabla\cdot\mathbf{j}+\frac{\partial\rho}{\partial t}=0. }[/math]

Вывод уравнения непрерывности

Дивергенция от ротора равна нулю, поэтому для четвёртого уравнения Максвелла (Закон Ампера—Максвелла) в системе СИ имеем:

[math]\displaystyle{ 0=\nabla\cdot[\nabla\times \mathbf{H}] = \nabla\cdot\mathbf{j}+\frac{\partial \nabla\cdot \mathbf{D}}{\partial t}= \nabla\cdot\mathbf{j} + \frac{\partial \rho}{\partial t}, }[/math]

где в последнем равенстве подставлено первое уравнение (Закон Гаусса).

Это уравнение при помощи интегральной теоремы Остроградского—Гаусса можно записать в следующем виде:

[math]\displaystyle{ \oint\limits_\mathbf{s} \mathbf{j}\cdot d\mathbf{s} = - \frac{d}{d t}\int\limits_v \rho\, dv. }[/math]

В левой части уравнения находится полный ток, протекающий через замкнутую поверхность [math]\displaystyle{ s\ }[/math]. В правой части — изменение со временем заряда, находящегося внутри объёма [math]\displaystyle{ v\ }[/math]. Таким образом, изменение заряда внутри объёма возможно только при его притоке или оттоке через поверхность [math]\displaystyle{ s\ }[/math], ограничивающую объём.

Уравнение непрерывности, эквивалентное закону сохранения заряда, далеко выходит за пределы классической электродинамики, оставаясь справедливым и в квантовой теории. Поэтому это уравнение само по себе может быть положено в основу электромагнитной теории. Тогда, например, ток смещения (производная по времени электрического поля) должен обязательно присутствовать в законе Ампера.

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

Закон сохранения энергии

Если умножить третье уравнение Максвелла в дифференциальной форме (закон Фарадея) скалярно на [math]\displaystyle{ \mathbf H }[/math], а четвёртое (закон Ампера — Максвелла) на [math]\displaystyle{ -\mathbf E }[/math] и сложить результаты, можно получить теорему Пойнтинга:

[math]\displaystyle{ \nabla\cdot\mathbf{S}+\frac{\partial}{\partial t}\left(w_E+w_H\right)=-\mathbf{E}\mathbf{j}, }[/math]

где

СГС	СИ
[math]\displaystyle{ \mathbf{S}=\frac{c}{4\pi}\mathbf E \times\mathbf H }[/math] [math]\displaystyle{ w_E = \frac{1}{8\pi}\,\mathbf{E}\cdot\mathbf{D} }[/math] [math]\displaystyle{ w_H = \frac{1}{8\pi}\,\mathbf{H}\cdot\mathbf{B} }[/math]	[math]\displaystyle{ \mathbf{S}=\mathbf E \times\mathbf H }[/math] [math]\displaystyle{ w_E = \frac{1}{2}\,\mathbf{E}\cdot \mathbf{D} }[/math] [math]\displaystyle{ w_H = \frac{1}{2}\,\mathbf{H}\cdot\mathbf{B} }[/math]

Вывод теоремы Пойнтинга

При помощи третьего и четвёртого уравнения Максвелла в дифференциальной форме в системе СИ можно получить:

[math]\displaystyle{ \mathbf{E}\cdot[\nabla\times\mathbf{H}] = \mathbf{E}\mathbf{j} + \mathbf{E}\cdot \frac{\partial \mathbf{D}}{\partial t},~~~~~~~~~~~~~ \mathbf{H}\cdot[\nabla\times\mathbf{E}] = -\mathbf{H}\cdot \frac{\partial \mathbf{B}}{\partial t}. }[/math]

Разница левых частей уравнений сворачивается по следующей формуле векторного анализа (производная произведения):

[math]\displaystyle{ \nabla\cdot[\mathbf{E}\times\mathbf{H}]=\mathbf{H}\cdot[\nabla\times\mathbf{E}]-\mathbf{E}\cdot[\nabla\times\mathbf{H}]. }[/math]

В линейных, но, возможно, неизотропных средах между напряжённостями и индукциями существует линейная связь. Например, для электрического поля [math]\displaystyle{ \mathbf{D}=\varepsilon_0\,\hat{\varepsilon}\cdot\mathbf{E} }[/math]. Если [math]\displaystyle{ \hat{\varepsilon} }[/math] — симметричная матрица, не зависящая от времени, то:

[math]\displaystyle{ \mathbf{E}\cdot \frac{\partial \mathbf{D}}{\partial t} = \varepsilon_0\, \mathbf{E}\cdot\hat{\varepsilon}\cdot\frac{\partial \mathbf{E}}{\partial t} =\frac{\varepsilon_0}{2}\,\frac{\partial (\mathbf{E}\cdot\hat{\varepsilon}\cdot\mathbf{E})}{\partial t} =\frac{1}{2}\frac{\partial (\mathbf{E}\cdot\mathbf{D})}{\partial t}. }[/math]

Аналогично для магнитного поля.

Вектор [math]\displaystyle{ \mathbf{S} }[/math] называется вектором Пойнтинга (вектором плотности потока электромагнитной энергии) и определяет количество электромагнитной энергии, переносимой через единицу площади в единицу времени. Интеграл вектора Пойнтинга по сечению распространяющейся волны определяет её мощность. Важно отметить, что, как впервые указал Хевисайд, физический смысл потока энергии имеет только безвихревая часть вектора Пойнтинга. Вихревая часть, дивергенция которой равна нулю, не связана с переносом энергии. Заметим, что Хевисайд получил выражение для закона сохранения независимо от Пойнтинга. В русскоязычной литературе вектор Пойнтинга часто называется также «вектором Умова — Пойнтинга».

Величины [math]\displaystyle{ w_E\ }[/math] и [math]\displaystyle{ w_H\ }[/math] определяют объёмные плотности энергии, соответственно, электрического и магнитного полей. При отсутствии токов и связанных с ними потерь теорема Пойнтинга является уравнением непрерывности для энергии электромагнитного поля. Проинтегрировав его в этом случае по некоторому замкнутому объёму и воспользовавшись теоремой Остроградского — Гаусса, можно получить закон сохранения энергии для электромагнитного поля:

[math]\displaystyle{ \oint\limits_s\mathbf{S}\cdot d\mathbf{s}+\frac{d}{dt}\int\limits_v(w_E+w_H)\,dv=0. }[/math]

Это уравнение показывает, что при отсутствии внутренних потерь изменение энергии электромагнитного поля в объёме происходит только за счёт мощности электромагнитного излучения, переносимого через границу этого объёма.

Вектор Пойнтинга связан с импульсом электромагнитного поля:

[math]\displaystyle{ \mathbf{p} = \frac{1}{c^2}\int \mathbf{S}\,dv, }[/math]

где интегрирование производится по всему пространству. Электромагнитная волна, поглощаясь или отражаясь от некоторой поверхности, передаёт ей часть своего импульса, что проявляется в форме светового давления. Экспериментально этот эффект впервые наблюдался П. Н. Лебедевым в 1899 году.

Потенциалы

Скалярный и векторный потенциалы

Закон Фарадея и закон Гаусса для магнитной индукции выполняются тождественно, если электрическое и магнитное поля выразить через скалярный [math]\displaystyle{ \varphi }[/math] и векторный [math]\displaystyle{ \mathbf{A} }[/math] потенциалы:

СГС	СИ
[math]\displaystyle{ \mathbf{E}=-\nabla \varphi-\frac{1}{c}\,\frac{\partial\mathbf{A}}{\partial t} }[/math] [math]\displaystyle{ \mathbf{B}=\nabla\times\mathbf{A} }[/math]	[math]\displaystyle{ \mathbf{E}=-\nabla \varphi - \frac{\partial \mathbf{A}}{\partial t} }[/math] [math]\displaystyle{ \mathbf{B}=\nabla\times\mathbf{A} }[/math]

Доказательство

Поскольку согласно закону Гаусса дивергенция индукции магнитного поля [math]\displaystyle{ \mathbf{B} }[/math] равна нулю, то по теореме Гельмгольца существует такое векторное поле [math]\displaystyle{ \mathbf{A} }[/math], что [math]\displaystyle{ \mathbf{B}= \nabla\times\mathbf{A}. }[/math] Тогда ротор вектора [math]\displaystyle{ \mathbf{E}'=\mathbf{E}+\frac{1}{c}\frac{\partial\mathbf{A}}{\partial t} }[/math] (в системе СГС) или вектора [math]\displaystyle{ \mathbf{E}'=\mathbf{E}+\frac{\partial\mathbf{A}}{\partial t} }[/math] (в системе СИ) удовлетворяет условию [math]\displaystyle{ \nabla\times\mathbf{E}'=0. }[/math] Например, в системе СИ получаем:

[math]\displaystyle{ \nabla\times\mathbf{E'}=\nabla\times\left(\mathbf{E}+\frac{\partial\mathbf{A}}{\partial t}\right) =\nabla\times\mathbf{E}+\nabla\times\frac{\partial\mathbf{A}}{\partial t} =\nabla\times\mathbf{E}+\frac{\partial}{\partial t}\left(\nabla\times\mathbf{A}\right) = -\frac{\partial\mathbf{B}}{\partial t}+\frac{\partial\mathbf{B}}{\partial t}=0. }[/math]

Из условия же равенства нулю ротора согласно теореме Гельмгольца следует, что существует такая скалярная функция [math]\displaystyle{ \varphi, }[/math] что [math]\displaystyle{ \mathbf{E}'=-\nabla\varphi. }[/math]

Обратная подстановка работает аналогично. Если электрическое и магнитное поля выражены через скалярный [math]\displaystyle{ \varphi }[/math] и векторный [math]\displaystyle{ \mathbf{A} }[/math] потенциалы согласно приведённым выше формулам, то дивергенция индукции магнитного поля автоматически равна нулю:

[math]\displaystyle{ \nabla\cdot\mathbf{B}= \nabla\cdot[\nabla\times\mathbf{A}]=[\nabla\times\nabla]\cdot\mathbf{A} = 0. }[/math]

Для напряжённости же электрического поля [math]\displaystyle{ \mathbf{E} }[/math] будет автоматически выполнен закон Фарадея. Например, в системе СИ получаем:

[math]\displaystyle{ \nabla\times\mathbf{E} =\nabla\times\left(-\nabla \varphi - \frac{\partial \mathbf{A}}{\partial t}\right) = -\nabla\times \frac{\partial \mathbf{A}}{\partial t} = -\frac{\partial\mathbf{B}}{\partial t}. }[/math]

При данных электрическом [math]\displaystyle{ \mathbf{E} }[/math] и магнитном [math]\displaystyle{ \mathbf{B} }[/math] полях скалярный и векторный потенциалы определены неоднозначно. Если [math]\displaystyle{ \psi\ }[/math] — произвольная функция координат и времени, то следующее преобразование не изменит значение полей:

СГС	СИ
[math]\displaystyle{ \mathbf{A} \rightarrow \mathbf{A} +\nabla \psi }[/math] [math]\displaystyle{ \varphi \rightarrow \varphi - \frac{1}{c}\,\frac{\partial \psi}{\partial t} }[/math]	[math]\displaystyle{ \mathbf{A} \rightarrow \mathbf{A} +\nabla \psi }[/math] [math]\displaystyle{ \varphi \rightarrow \varphi - \frac{\partial \psi}{\partial t}. }[/math]

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

Неоднозначность определения потенциалов оказывается удобной для наложения на них дополнительных условий, называемых калибровкой. Благодаря этому уравнения электродинамики принимают более простой вид. Рассмотрим, например, уравнения Максвелла в однородных и изотропных средах с диэлектрической ([math]\displaystyle{ \varepsilon }[/math]) и магнитной ([math]\displaystyle{ \mu\ }[/math]) проницаемостями. Для данных [math]\displaystyle{ \varphi }[/math] и [math]\displaystyle{ \mathbf{A} }[/math] всегда можно подобрать такую функцию [math]\displaystyle{ f }[/math], чтобы выполнялось калибровочное условие Лоренца:

СГС	СИ
[math]\displaystyle{ \nabla\cdot\mathbf{A} + \frac{\varepsilon\mu}{c}\,\frac{\partial \varphi}{\partial t} = 0 }[/math]	[math]\displaystyle{ \nabla\cdot\mathbf{A} + \frac{\varepsilon\mu}{c^2}\,\frac{\partial \varphi}{\partial t} = 0. }[/math]

В этом случае оставшиеся уравнения Максвелла в однородных и изотропных средах могут быть записаны в следующем виде:

СГС	СИ
[math]\displaystyle{ \square\varphi = -4\pi\,\frac{\rho}{\varepsilon} }[/math] [math]\displaystyle{ \square\mathbf{A} = -\frac{4\pi}{c}\,\mu\,\mathbf{j} }[/math]	[math]\displaystyle{ \square\varphi = -\frac{\rho}{\varepsilon\varepsilon_0} }[/math] [math]\displaystyle{ \square\mathbf{A} = -\mu\mu_0\,\mathbf{j}, }[/math]

где [math]\displaystyle{ \square }[/math] — оператор Д’Аламбера, который и в системе СГС, и в системе СИ имеет вид:

[math]\displaystyle{ \square = \Delta -\frac{\varepsilon\mu}{c^2}\frac{\partial^2}{\partial t^2}. }[/math]

Таким образом, 8 уравнений Максвелла (уравнения первого порядка) для компонент электромагнитного поля (2 векторных и 2 скалярных) при помощи потенциалов могут быть сведены к 4 уравнениям, но уже второго порядка (скалярному для [math]\displaystyle{ \varphi }[/math] и векторному для [math]\displaystyle{ \mathbf{A} }[/math]). Решения этих уравнений для произвольно движущегося точечного заряда называются потенциалами Лиенара — Вихерта.

Возможно введение других калибровок. Так, для решения ряда задач удобной оказывается кулоновская калибровка:

[math]\displaystyle{ \nabla\cdot\mathbf{A} = 0 }[/math]

В этом случае:

СГС	СИ
[math]\displaystyle{ \Delta\varphi = -4\pi\,\frac{\rho}{\varepsilon} }[/math] [math]\displaystyle{ \square\mathbf{A} = -\frac{4\pi}{c}\,\mu\,\mathbf{j}_\bot }[/math]	[math]\displaystyle{ \Delta\varphi = -\frac{\rho}{\varepsilon\varepsilon_0} }[/math] [math]\displaystyle{ \square\mathbf{A} = -\mu\mu_0\,\mathbf{j}_\bot, }[/math]

,

где [math]\displaystyle{ \mathbf{j}_\bot }[/math] — соленоидальная часть тока ([math]\displaystyle{ \nabla\mathbf\cdot{j}_\bot=0 }[/math]).

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

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

Векторы Герца

• В 1887 году Генрих Герц предложил вместо непосредственного решения уравнений Максвелла для двух векторных функций электрического и магнитного полей или скалярного и векторного потенциалов перейти к новой единственной векторной функции, которая носит теперь имя электрического вектора Герца [math]\displaystyle{ {\mathbf\Pi}^e }[/math] и позволяет в некоторых случаях упростить решение электродинамических задач, сводя их к решению скалярного волнового уравнения.

СГС	СИ
[math]\displaystyle{ \varphi=-\nabla\cdot\mathbf{\Pi}^e, }[/math] [math]\displaystyle{ \mathbf{A}=\frac{\varepsilon\mu}{c}\frac{\partial\mathbf{\Pi}^e}{\partial t}, }[/math]	[math]\displaystyle{ \varphi=-\nabla\cdot\mathbf{\Pi}^e, }[/math] [math]\displaystyle{ \mathbf{A}=\frac{\varepsilon\mu}{c^2}\frac{\partial \mathbf{\Pi}^e}{\partial t}, }[/math]
[math]\displaystyle{ \mathbf {E}^e = \nabla(\nabla\cdot \mathbf{\Pi}^e) - \frac{\mu\varepsilon}{c^2} \frac{\partial^2 \mathbf{\Pi}^e}{\partial t^2}, }[/math] [math]\displaystyle{ \mathbf {B}^e = \frac{\mu\varepsilon}{c} \nabla\times \frac{\partial \mathbf{\Pi}^e}{\partial t}. }[/math]	[math]\displaystyle{ \mathbf {E}^e = \nabla(\nabla\cdot \mathbf{\Pi}^e) - \frac{\mu\varepsilon}{c^2} \frac{\partial^2 \mathbf{\Pi}^e}{\partial t^2}, }[/math] [math]\displaystyle{ \mathbf {B}^e = \frac{\mu\varepsilon}{c^2} \nabla\times \frac{\partial \mathbf{\Pi}^e}{\partial t}. }[/math]

Заметим, что скалярный [math]\displaystyle{ \varphi }[/math] и векторный [math]\displaystyle{ \mathbf{A} }[/math] потенциалы, выраженные через вектор Герца, автоматически удовлетворяют калибровочному условию Лоренца. Вектор Герца учитывает все поля, связанные со свободными зарядами и их токами.

Подставляя выражения для полей через электрический вектор в два последних уравнения Максвелла, можно получить:

СГС	СИ
[math]\displaystyle{ \Delta \mathbf {\Pi}^e - \frac{\varepsilon\mu}{c^2}\frac{\partial^2{\mathbf {\Pi}^e}}{\partial t^2} =- \frac{4\pi}{\varepsilon}\mathbf{P}_f , }[/math] [math]\displaystyle{ \Delta \mathbf {\Pi}^e - \frac{\mu}{c^2}\frac{\partial^2{\mathbf {\Pi}^e}}{\partial t^2} =- \frac{4\pi}{\varepsilon} \left[\mathbf {P}+\mathbf{P}_f\right]. }[/math]	[math]\displaystyle{ \Delta \mathbf {\Pi}^e - \frac{\varepsilon\mu}{c^2}\frac{\partial^2{\mathbf {\Pi}^e}}{\partial t^2} =- \frac{1}{\varepsilon\varepsilon_0}\mathbf {P}_f, }[/math] [math]\displaystyle{ \Delta \mathbf {\Pi}^e - \frac{\mu}{c^2}\frac{\partial^2{\mathbf {\Pi}^e}}{\partial t^2} =- \frac{1}{\varepsilon\varepsilon_0} \left[\mathbf {P}+\mathbf {P}_f\right]. }[/math]

Здесь введён вектор поляризации свободных зарядов и токов:

[math]\displaystyle{ \mathbf j_f = \frac{\partial }{\partial t} {\mathbf P}_f,~~~~~~~~~~~~~~~~ \rho_f = - \nabla\cdot {\mathbf P}_f, }[/math]

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

Таким образом, электрический вектор Герца определяется волновыми уравнениями, в правой части которых стоит поляризуемость, обусловленная свободными либо свободными и связанными зарядами, то есть электрическими дипольными моментами.

• В 1901 году парный электрическому вектору Герца магнитный вектор, который также традиционно называют именем Герца, ввёл итальянский физик Аугусто Риги.

СГС	СИ
[math]\displaystyle{ \varphi=0, }[/math] [math]\displaystyle{ \mathbf{A}=\nabla\times\mathbf{\Pi}^m, }[/math]	[math]\displaystyle{ \varphi=0, }[/math] [math]\displaystyle{ \mathbf{A}=\frac{1}{c}\nabla\times\mathbf{\Pi}^m, }[/math]
[math]\displaystyle{ \mathbf {D}^m = -\frac{\varepsilon\mu}{c} \nabla\times \frac{\partial \mathbf{\Pi}^m}{\partial t}, }[/math] [math]\displaystyle{ \mathbf {H}^m = \nabla\times[\nabla\times \mathbf{\Pi}^m] . }[/math]	[math]\displaystyle{ \mathbf {D}^m = -\frac{\varepsilon\mu}{c^2} \nabla\times \frac{\partial \mathbf{\Pi}^m}{\partial t}, }[/math] [math]\displaystyle{ \mathbf {H}^m = \nabla\times[\nabla\times \mathbf{\Pi}^m] . }[/math]

Поскольку поля, описываемые магнитным вектором Герца, не зависят от свободных зарядов и токов, а магнитные монополи не обнаружены, потенциалы удовлетворяют калибровке Лоренца в вырожденном виде — так называемой кулоновской калибровке ([math]\displaystyle{ \varphi=0 }[/math], [math]\displaystyle{ \nabla\cdot\mathbf{A}=0 }[/math]).

Аналогичным образом можно получить уравнения для магнитного потенциала Герца, подставляя выраженные через него поля в третье и четвёртое уравнения Максвелла без тока:

СГС	СИ
[math]\displaystyle{ \Delta \mathbf {\Pi}^m - \frac{\varepsilon\mu}{c^2}\frac{\partial^2{\mathbf {\Pi}^m}}{\partial t^2} = 0, }[/math] [math]\displaystyle{ \Delta \mathbf {\Pi}^m - \frac{\varepsilon}{c^2}\frac{\partial^2{\mathbf {\Pi}^m}}{\partial t^2} =- \frac{4\pi}{\mu} \mathbf {M} . }[/math]	[math]\displaystyle{ \Delta \mathbf {\Pi}^m - \frac{\varepsilon\mu}{c^2}\frac{\partial^2{\mathbf {\Pi}^m}}{\partial t^2} =0, }[/math] [math]\displaystyle{ \Delta \mathbf {\Pi}^m - \frac{\varepsilon}{c^2}\frac{\partial^2{\mathbf {\Pi}^m}}{\partial t^2} =- \frac{1}{\mu} \mathbf {M}. }[/math]

Действие сторонних магнитных полей, связанных с внешними источниками, может быть учтено по аналогии с электрическим вектором Герца введением в правые части дополнительной магнитной поляризации [math]\displaystyle{ \mathbf {M}_f }[/math].

Таким образом, выделяется два типа электромагнитных полей, выражающихся через электрический и магнитный потенциалы Герца, а произвольное поле можно представить в виде суммы таких полей. Поля, выражающиеся через электрический вектор Герца, носят название полей электрического типа, или поперечно-магнитных (TM) полей, поскольку магнитное поле для них ортогонально направлению вектора Герца. Соответственно, поля, выражающиеся через магнитный вектор Герца, называют полями магнитного типа, или поперечно-электрическими полями (TE), электрическое поле в которых ортогонально порождающему вектору Герца. Поля TM можно представить как поля, порождаемые распределёнными в пространстве электрическими диполями, а поля TE, соответственно, магнитными. Векторные потенциалы Герца, в свою очередь, могут быть во многих случаях выражены через скалярные потенциалы.

Потенциалы Дебая

В электродинамике широко используются скалярные потенциалы, предложенные Дебаем.

Волновое уравнение представляет собой систему трёх связанных скалярных уравнений, которые распадаются на три скалярных уравнения Гельмгольца только в декартовой системе координат. Для удобства поиска решений, удовлетворяющих граничным условиям, желательно выбирать координатные системы, координатные поверхности которых близки или совпадают с поверхностями границ. Один из подходов к решению векторного уравнения Гельмгольца состоит во введении скалярных функций [math]\displaystyle{ \psi }[/math], удовлетворяющих скалярному волновому уравнению Гельмгольца, через которые затем могут быть выражены векторные поля:

[math]\displaystyle{ \nabla^2 {\psi}+ k^2 \psi=0. }[/math]

[math]\displaystyle{ {\mathbf M}_\psi=\nabla\times ({\mathbf f}\psi), }[/math]

[math]\displaystyle{ {\mathbf N}_\psi = \frac{1}{k}\nabla\times\nabla \times({\mathbf f}\psi), }[/math]

[math]\displaystyle{ {\mathbf L}_\psi=\nabla\psi. }[/math]

Здесь [math]\displaystyle{ \mathbf f }[/math] — некоторая векторная функция координат. Вектор [math]\displaystyle{ {\mathbf L}_\psi }[/math] описывает потенциальную часть поля, и его можно положить равным нулю при отсутствии свободных зарядов.

Если для некоторой ортогональной координатной системы существует функция [math]\displaystyle{ {\mathbf f}({\mathbf r}) }[/math], пропорциональная координатному вектору, то произвольное векторное поле, удовлетворяющее векторному уравнению Гельмгольца в этой системе, можно представить в виде суммы векторных функций, пропорциональных векторам [math]\displaystyle{ {\mathbf M}_\psi }[/math] и [math]\displaystyle{ {\mathbf N}_\psi }[/math]. Как следует из уравнений Максвелла, электрическому полю, пропорциональному [math]\displaystyle{ {\mathbf M}_\psi }[/math], соответствует магнитное поле типа [math]\displaystyle{ {\mathbf N}_\psi }[/math], и наоборот. При этом векторные потенциалы [math]\displaystyle{ {\mathbf f}\psi }[/math] соответствуют векторам Герца. Поскольку в этом случае поле, пропорциональное [math]\displaystyle{ {\mathbf M}_\psi }[/math], нормально вектору [math]\displaystyle{ {\mathbf f} }[/math], его компоненты являются тангенциальными к соответствующей [math]\displaystyle{ {\mathbf f} }[/math] координатной поверхности. Если границы в решаемой задаче совпадают с одной из таких координатных поверхностей, то удовлетворение граничным условиям существенно упрощается.

Такое представление возможно только в ограниченном числе ортогональных координатных систем. В декартовой системе координат в качестве вектора [math]\displaystyle{ {\mathbf f} }[/math] может выступать любой координатный вектор. Соответствующие решения представляют собой плоские волны. Для цилиндрической системы координат [math]\displaystyle{ {\mathbf f}={\mathbf i}_z }[/math], для сферической [math]\displaystyle{ {\mathbf f}={\mathbf r} }[/math]. Кроме того, такое представление возможно в конической, а также относительно оси [math]\displaystyle{ z }[/math] в параболической и эллиптической цилиндрических системах координат.

Векторы Римана — Зильберштейна

Если ввести комплексный вектор Римана — Зильберштейна [math]\displaystyle{ \mathbf F }[/math] и комплексно сопряжённый ему вектор [math]\displaystyle{ \mathbf F^* }[/math]:

СГС	СИ
[math]\displaystyle{ \mathbf F=\frac{1}{\sqrt{8\pi}}\left[\sqrt{\varepsilon}\mathbf E+i\sqrt{\mu}\mathbf H\right], }[/math]	[math]\displaystyle{ \mathbf F=\frac{1}{\sqrt{2}}\left[\sqrt{\varepsilon\varepsilon_0}\mathbf E+i\sqrt{\mu\mu_0}\mathbf H\right], }[/math]

то уравнения Максвелла сводятся к двум:

СГС	СИ
[math]\displaystyle{ \nabla\cdot \mathbf{F}=\sqrt{\frac{2\pi}{\varepsilon}}\rho, }[/math] [math]\displaystyle{ \nabla\times \mathbf{F}=i\sqrt{2\pi\mu}\mathbf{j}+i\frac{\sqrt{\varepsilon\mu}}{c}\frac{\partial \mathbf F}{\partial t}. }[/math]	[math]\displaystyle{ \nabla\cdot \mathbf{F}=\sqrt{\frac{1}{2\varepsilon\varepsilon_0}}\rho, }[/math] [math]\displaystyle{ \nabla\times \mathbf{F}=i\sqrt{\frac{\mu\mu_0}{2}}\mathbf{j}+i\frac{\sqrt{\varepsilon\mu}}{c}\frac{\partial \mathbf F}{\partial t}. }[/math]

При отсутствии сторонних зарядов и токов остаётся только второе уравнение (первое из-за равенства дивергенции ротора нулю в этом случае удовлетворяется автоматически с точностью до не зависящей от времени компоненты):

[math]\displaystyle{ \nabla\times \mathbf F= i\frac{\sqrt{\varepsilon\mu}}{c}\frac{\partial \mathbf F}{\partial t}. }[/math]

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

Для гармонического поля с зависимостью [math]\displaystyle{ \mathbf F= \mathbf F^\pm e^{\pm i\omega t} }[/math] вектор [math]\displaystyle{ \mathbf F }[/math] является собственным вектором оператора ротора:

[math]\displaystyle{ \nabla\times\mathbf F^\pm=\mp k\mathbf F^\pm. }[/math]

При выбранной нормировке [math]\displaystyle{ \mathbf F }[/math] имеет смысл комплексной амплитуды электромагнитного поля, а его квадрат модуля

[math]\displaystyle{ w=|\mathbf F|^2\equiv\mathbf F^*\mathbf F=w_E+w_H }[/math]

имеет смысл плотности энергии поля.

Вектор Пойнтинга:

[math]\displaystyle{ \mathbf S=-\frac{ic}{\sqrt{\varepsilon\mu}}{\mathbf F}^*\times{\mathbf F}. }[/math]

Векторы [math]\displaystyle{ \mathbf F }[/math] и [math]\displaystyle{ \mathbf F^* }[/math] можно интерпретировать как волновые функции циркулярно поляризованных фотонов.

Ковариантная формулировка

С современной точки зрения четырёхмерная ковариантная формулировка электродинамики (и в частности запись уравнений Максвелла в таком виде) является физически наиболее фундаментальной.

Практически она приводит, кроме явной ковариантности, к значительно большей компактности уравнений, а значит — к определённой красоте и в ряде случаев к удобству, и более органично и прямо включает в себя единство электромагнитного поля.

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

• Ниже сначала рассмотрен (как более простой) первый вариант — вариант лоренц-ковариантной формулировки в плоском пространстве-времени.

Четырёхмерные векторы

При ковариантной записи уравнений электродинамики производится переход от трёхмерных векторов и скаляров к четырёхмерным векторам (4-векторы). Независимо от системы единиц, четырёхмерные координаты (4-вектор координат, в компоненты которого входят время и трёхмерные пространственные координаты), производная по этим координатам (4-производная) и плотность тока определяются следующим образом:

[math]\displaystyle{ x^\alpha=(ct, ~\mathbf{r})=(ct,x,y,z), }[/math] [math]\displaystyle{ \partial_\alpha=\frac{\partial }{\partial x^\alpha}=\Bigl(\frac{1}{c}\,\frac{\partial}{\partial t}, ~\nabla\Bigr)=\Bigl(\frac{1}{c}\,\frac{\partial}{\partial t},\frac{\partial}{\partial x},\frac{\partial}{\partial y},\frac{\partial}{\partial z}\Bigr), }[/math] [math]\displaystyle{ j^\alpha=\rho \frac{dx^\alpha}{dt}=(c\rho, ~\mathbf{j})=(c\rho,j_x,j_y,j_z). }[/math]

Индекс 4-вектора принимает значения [math]\displaystyle{ \alpha=0,1,2,3\ }[/math]. В компонентной записи вектора сначала идёт нулевая (временная) компонента, затем — пространственные. Например, время равно [math]\displaystyle{ t=x^0/c\ }[/math], а плотность заряда [math]\displaystyle{ \rho=j^0/c\ }[/math]. В силу этих определений закон сохранения заряда в ковариантной форме принимает следующий вид:

[math]\displaystyle{ \partial_\alpha j^\alpha=0. }[/math]

По повторяющемуся индексу предполагается суммирование от 0 до 3 (правило Эйнштейна).

Пример

Приведенное выше уравнение является компактной записью уравнения непрерывности:

[math]\displaystyle{ \partial_0 j^0+\partial_1 j^1+\partial_2 j^2+\partial_3 j^3= \frac{\partial \rho}{\partial t}+\nabla\mathbf{j}=0. }[/math]

Введём 4-вектор потенциала, имеющий в системах СГС и СИ следующие компоненты:

СГС	СИ
[math]\displaystyle{ A^\alpha = (\varphi, ~~\mathbf{A}) }[/math] [math]\displaystyle{ A_\alpha = (\varphi, -\mathbf{A}) }[/math]	[math]\displaystyle{ A^\alpha = (\varphi/c, ~~\mathbf{A}) }[/math] [math]\displaystyle{ A_\alpha = (\varphi/c, -\mathbf{A}) }[/math]

При ковариантной записи играет роль положение индекса у 4-вектора. Если индекс находится внизу, то такой вектор называется ковариантным вектором (или ковектором), и его пространственные компоненты имеют обратный знак по сравнению с компонентами 4-вектора. Поднятие и опускание индексов проводится при помощи метрического тензора [math]\displaystyle{ A_\alpha=g_{\alpha\beta}A^\beta\ }[/math], который в четырёхмерном пространстве Минковского имеет диагональный вид с сигнатурой: [math]\displaystyle{ g_{\alpha\beta}=g^{\alpha\beta}={\rm diagonal}(1,-1,-1,-1)\ }[/math].

При помощи такого определения 4-вектора потенциала калибровочное условие Лоренца в ковариантной форме можно записать следующим образом:

[math]\displaystyle{ \partial_\alpha A^\alpha=0. }[/math]

Если это условие выполняется, то уравнения Максвелла для потенциалов в вакууме при наличии зарядов и токов принимают вид:

СГС	СИ
[math]\displaystyle{ \partial^2 A^\alpha = \frac{4\pi}{c}\,j^\alpha }[/math]	[math]\displaystyle{ \partial^2 A^\alpha = \mu_0\,j^\alpha }[/math],

где [math]\displaystyle{ \partial^2 }[/math] — оператор Даламбера с обратным знаком:

[math]\displaystyle{ \partial^2 = \partial_\alpha\partial^\alpha = -\square = \frac{1}{c^2} \frac{\partial^2}{\partial t^2} - \Delta. }[/math]

Нулевая компонента уравнений Максвелла для 4-вектора потенциала соответствует уравнению для [math]\displaystyle{ \varphi }[/math], а пространственная — для [math]\displaystyle{ \mathbf{A} }[/math].

Тензор электромагнитного поля

Определим ковариантный тензор электромагнитного поля при помощи производной от 4-вектора потенциала:

[math]\displaystyle{ F_{\alpha\beta}=\partial_\alpha A_\beta - \partial_\beta A_\alpha. }[/math]

Явный вид этого антисимметричного тензора ([math]\displaystyle{ F_{\alpha\beta} =-F_{\beta\alpha}\ }[/math]) может быть представлен в следующем виде:

СГС	СИ
[math]\displaystyle{ F_{\alpha\beta}=\left(\begin{matrix} 0 & E_x & E_y & E_z \\ -E_x & 0 & -B_z & B_y \\ -E_y & B_z & 0 & -B_x \\ -E_z & -B_y & B_x & 0 \end{matrix}\right), }[/math]	[math]\displaystyle{ F_{\alpha\beta}=\left(\begin{matrix} 0 & E_x/c & E_y/c & E_z/c \\ -E_x/c & 0 & -B_z & B_y \\ -E_y/c & B_z & 0 & -B_x \\ -E_z/c & -B_y & B_x & 0 \end{matrix}\right). }[/math]

Временные компоненты тензора составлены из компонент напряжённости электрического поля, а пространственные — магнитного, что может быть записано следующим образом: [math]\displaystyle{ F_{\alpha\beta}=(\mathbf{E},\mathbf{B}) }[/math]. В тензоре электромагнитного поля с верхними индексами изменяется знак у нулевых компонент (то есть перед компонентами электрического поля): [math]\displaystyle{ F^{\alpha\beta}=(-\mathbf{E},\mathbf{B}) }[/math].

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

[math]\displaystyle{ \partial_\alpha F_{\beta \gamma}+\partial_\beta F_{\gamma\alpha}+\partial_\gamma F_{\alpha \beta}=0. }[/math]

Его можно переписать в более компактном виде, введя дуальный тензор электромагнитного поля:

[math]\displaystyle{ \partial_\alpha \tilde{F}^{\alpha \beta} = 0,~~~~~~~~~~~~~~~~~~~\tilde{F}^{\alpha \beta} =\frac{1}{2}\epsilon^{\alpha\beta\gamma\delta}\, F_{\gamma\delta}, }[/math]

где [math]\displaystyle{ \epsilon^{\alpha\beta\gamma\delta}\ }[/math] — антисимметричный символ Леви-Чивиты ([math]\displaystyle{ \epsilon^{0123}=1\ }[/math]). Это уравнение является ковариантной записью закона Гаусса для магнитного поля и закона электромагнитной индукции Фарадея. Компоненты дуального тензора [math]\displaystyle{ \tilde{F}_{\alpha\beta} }[/math] получаются из тензора [math]\displaystyle{ F_{\alpha\beta} }[/math] в результате перестановки электрического и магнитного полей: [math]\displaystyle{ \tilde{F}_{\alpha\beta}=(\mathbf{B},-\mathbf{E}) }[/math], [math]\displaystyle{ \tilde{F}^{\alpha\beta}=(-\mathbf{B},-\mathbf{E}) }[/math].

Полная система уравнений Максвелла в ковариантной форме имеет вид:

СГС	СИ
[math]\displaystyle{ \partial_\alpha \tilde{F}^{\alpha \beta} = 0 }[/math] [math]\displaystyle{ \partial_\alpha F^{\alpha \beta}=\frac{4\pi}{c}j^\beta, }[/math]	[math]\displaystyle{ \partial_\alpha \tilde{F}^{\alpha \beta} = 0 }[/math] [math]\displaystyle{ \partial_\alpha F^{\alpha \beta}=\mu_0\,j^\beta. }[/math]

По повторяющемуся индексу [math]\displaystyle{ \alpha\ }[/math] проводится суммирование от 0 до 3, а в правой части второго уравнения находится 4-вектор тока. Нулевая компонента этого уравнения соответствует закону Гаусса, а пространственные — закону Ампера — Максвелла.

При помощи тензора электромагнитного поля можно получить законы преобразований компонент электрического и магнитного полей, измеряемых относительно различных инерциальных систем отсчёта:

СГС	СИ
[math]\displaystyle{ E'_y = \gamma\, (E_y - {u \over c} \,B_z),~~~ E'_z = \gamma\,(E_z + {u \over c} \,B_y), }[/math] [math]\displaystyle{ B'_y = \gamma\,(B_y + {u \over c}\, E_z),~~~ B'_z = \gamma\,(B_z - {u \over c} \,E_y), }[/math]	[math]\displaystyle{ E'_y = \gamma\,(E_y - u \,B_z),~~~ E'_z = \gamma\,(E_z + u \, B_y), }[/math] [math]\displaystyle{ B'_y = \gamma\,(B_y + {u \over c^2}\,E_z),~~ B'_z = \gamma\,(B_z - {u \over c^2}\,E_y), }[/math]

где «штрихованные» величины измеряются относительно системы отсчёта, движущейся вдоль оси [math]\displaystyle{ x }[/math] со скоростью [math]\displaystyle{ u }[/math] относительно системы, в которой измеряются «не штрихованные» компоненты полей, а [math]\displaystyle{ \gamma=1/\sqrt{1-u^2/c^2} }[/math] — фактор Лоренца. Компоненты полей вдоль направления относительного движения инерциальных систем отсчёта остаются неизменными: [math]\displaystyle{ E'_x = E_x,~~B'_x = B_x }[/math].

Уравнения Максвелла в вакууме инвариантны относительно преобразований Лоренца. Это послужило одним из толчков к созданию специальной теории относительности.

Электрическое и магнитное поля различным образом изменяются при инверсии осей пространственной системы координат. Электрическое поле является полярным вектором, а магнитное — аксиальным вектором. Можно построить две инвариантные относительно преобразований Лоренца величины:

[math]\displaystyle{ F_{\alpha\beta}F^{\alpha\beta}=\mathrm{inv}~~~~~~~~~~~~\varepsilon^{\alpha\beta\gamma\delta}F_{\alpha\beta}F_{\gamma\delta}=\mathrm{inv}. }[/math]

Первый инвариант является скаляром, а второй — псевдоскаляром, то есть изменяет свой знак при инверсии координатных осей.

Лагранжиан

Действие [math]\displaystyle{ S\ }[/math] и лагранжиан (функция Лагранжа) [math]\displaystyle{ L\ }[/math] для пробного заряда, движущегося во внешнем электромагнитном поле, в системе СГС и СИ имеют вид :

СГС	СИ
[math]\displaystyle{ S = \int(-mc\,ds ~-~ \frac{q}{c} \,A_\alpha\, dx^\alpha) }[/math] [math]\displaystyle{ L= - mc^2\sqrt{1-\frac{\mathbf{u^2}}{c^2}} + \frac{q}{c}\, \mathbf{A}\cdot\mathbf{u}-q\varphi }[/math]	[math]\displaystyle{ S = \int(-mc\,ds ~-~ q \,A_\alpha\, dx^\alpha) }[/math] [math]\displaystyle{ L= - mc^2\sqrt{1-\frac{\mathbf{u^2}}{c^2}} + q\, \mathbf{A}\cdot\mathbf{u}-q\varphi }[/math]

где:

• [math]\displaystyle{ m\ }[/math] — масса частицы (в единицах СИ — кг);
• [math]\displaystyle{ \mathbf{u} }[/math] — её скорость (в единицах СИ — м/с);
• [math]\displaystyle{ q\ }[/math] — заряд частицы (в единицах СИ — Кл);
• [math]\displaystyle{ ds=\sqrt{dx_\alpha dx^\alpha} }[/math] — 4-х интервал.

Уравнения движения заряда под воздействием силы Лоренца в ковариантной записи имеют вид:

СГС	СИ
[math]\displaystyle{ mc \frac{d^2x^\alpha}{ds^2} = \frac{q}{c}\, F^{\alpha\beta}\, \frac{dx_\beta}{ds} }[/math]	[math]\displaystyle{ mc \frac{d^2x^\alpha}{ds^2} = q\, F^{\alpha\beta}\, \frac{dx_\beta}{ds}. }[/math]

Уравнения Максвелла получаются из принципа наименьшего действия, в котором динамическими переменными являются 4-потенциалы электромагнитного поля [math]\displaystyle{ A^\alpha\ }[/math]. При этом используется следующее ковариантное выражение для действия:

СГС	СИ
[math]\displaystyle{ S = -\frac{1}{16\pi c} \int F_{\alpha\beta}F^{\alpha\beta}\, d^4x - \frac{1}{c^2} \int A_\alpha\,j^\alpha\,d^4x }[/math]	[math]\displaystyle{ S = -\frac{1}{4\mu_0 c} \int F_{\alpha\beta}F^{\alpha\beta}\, d^4x - \frac{1}{c} \int A_\alpha\,j^\alpha\,d^4x, }[/math]

где производится интегрирование по инвариантному 4-объёму [math]\displaystyle{ d^4x=c\,dt\,dv }[/math].

Запись при помощи дифференциальных форм

Основная статья: Дифференциальные формы в электромагнетизме

Уравнения Максвелла в ковариантной форме, аналогично векторному представлению в трёхмерном пространстве, можно записать в «безындексной форме». Для этого вводится операция внешнего произведения [math]\displaystyle{ \wedge }[/math], обладающая свойством антисимметричности: [math]\displaystyle{ \mathrm{d}x^\alpha\wedge \mathrm{d}x^\beta=-\mathrm{d}x^\beta\wedge \mathrm{d}x^\alpha }[/math]. Внешнее произведение позволяет записывать свёрнутые по всем индексам выражения с антисимметричными тензорами, такими как [math]\displaystyle{ F_{\alpha\beta}\ }[/math]. При этом возникают объекты, называемые дифференциальными формами (или просто формами). 1-форма потенциала поля определяется следующим образом (по индексу [math]\displaystyle{ \alpha\ }[/math] — сумма от 0 до 3):

[math]\displaystyle{ \mathrm{A} = A_\alpha \mathrm{d}x^\alpha\ . }[/math]

Из 1-формы при помощи операции внешнего дифференцирования [math]\displaystyle{ \mathrm{d}\ }[/math] получается 2-форма электромагнитного поля (или 2-форма Фарадея):

[math]\displaystyle{ \mathrm{F} = \mathrm{d}\,\mathrm{A} = \mathrm{d}A_\alpha \wedge \mathrm{d}x^\alpha = \partial_\beta A_\alpha ~\mathrm{d}x^\beta \wedge \mathrm{d}x^\alpha = \frac{1}{2} F_{\alpha\beta} ~\mathrm{d}x^\alpha\wedge \mathrm{d}x^\beta. }[/math]

Операция внешнего дифференцирования обладает свойством [math]\displaystyle{ \mathrm{d}^2=0\ }[/math], что приводит к закону Гаусса для магнитного поля и закону Фарадея:

[math]\displaystyle{ \mathrm{d}\,\mathrm{F} = \mathrm{d}^2\,\mathrm{A} = \frac{1}{2}\,\partial_\gamma F_{\alpha\beta}\,\mathrm{d}x^\gamma \wedge \mathrm{d}x^\alpha \wedge \mathrm{d}x^\beta = \frac{1}{6}\,(\partial_\gamma F_{\alpha\beta}+\partial_\alpha F_{\beta\gamma}+\partial_\beta F_{\gamma\alpha})~ \mathrm{d}x^\alpha \wedge \mathrm{d}x^\beta \wedge \mathrm{d}x^\gamma = 0. }[/math]

Для записи оставшихся уравнений Максвелла вводится дуальная к [math]\displaystyle{ \mathrm{F}\ }[/math] 2-форма [math]\displaystyle{ ^*\mathrm{F}\ }[/math], называемая также 2-формой Максвелла:

[math]\displaystyle{ ^*\mathrm{F} = \frac{1}{2}\, \tilde{F}_{\alpha\beta} ~ \mathrm{d}x^\alpha \wedge \mathrm{d}x^\beta= \frac{1}{4}\, F^{\alpha\beta} \,\epsilon_{\alpha\beta\gamma\delta}~ \mathrm{d}x^\gamma \wedge \mathrm{d}x^\delta, }[/math]

и 3-форма тока:

[math]\displaystyle{ ^*\mathrm{J} = -\frac{1}{3!}\,j^{\alpha} \,\epsilon_{\alpha\beta\gamma\delta}~ \mathrm{d}x^\beta \wedge \mathrm{d}x^\gamma \wedge \mathrm{d}x^\delta, }[/math]

где [math]\displaystyle{ \epsilon_{\alpha\beta\gamma\delta}\ }[/math] — абсолютный антисимметричный символ Леви-Чивиты ([math]\displaystyle{ \epsilon_{0123}=-1\ }[/math]). Свёртка с символом Леви-Чивиты внешнего произведения дифференциалов называется оператором звезды Ходжа.

В этих обозначениях уравнения Максвелла в системах СГС и СИ принимают следующий вид:

СГС	СИ
[math]\displaystyle{ \mathrm{d}\,\mathrm{F}=0, }[/math] [math]\displaystyle{ \mathrm{d}^*\mathrm{F}=\frac{4\pi}{c}\,^*\mathrm{J}, }[/math]	[math]\displaystyle{ \mathrm{d}\,\mathrm{F}=0, }[/math] [math]\displaystyle{ \mathrm{d}^*\mathrm{F}=\mu_0\,^*\mathrm{J}. }[/math]

Доказательство

Чтобы показать эквивалентность этих уравнений уравнениям Максвелла, необходимо записать их в трёхмерной векторной форме. В этом случае в системе СГС ток и 2-форма Максвелла имеют вид:

[math]\displaystyle{ ^*\mathrm{J}=c\rho\mathrm{d}v-\mathbf{j}\mathrm{d}\mathbf{s}\wedge \mathrm{d}t\, c,~~~~~~~~~~~~~~~~~ ^*\mathrm{F}= \mathbf{E}\mathrm{d}\mathbf{s}-\mathbf{B}\mathrm{d}\mathbf{r} \wedge \mathrm{d}t\,c, }[/math]

где [math]\displaystyle{ \mathrm{d}v=\mathrm{d}x\wedge\mathrm{d}y\wedge\mathrm{d}z }[/math] - трёхмерный объём, а [math]\displaystyle{ \mathrm{d}\mathbf{s}=(\mathrm{d}y\wedge\mathrm{d}z,~\mathrm{d}z\wedge\mathrm{d}x, \mathrm{d}x\wedge\mathrm{d}y) }[/math] - вектор площади поверхности в трёхмерном пространстве. Так как:

[math]\displaystyle{ \mathrm{d}^*\mathrm{F}=(\nabla\mathbf{E})\mathrm{d}v +\Bigl( \frac{1}{c}\frac{\partial\mathbf{E}}{\partial t}- \nabla\times\mathbf{B}\Bigr) \mathrm{d}\mathbf{s} \wedge \mathrm{d}t\,c. }[/math]

то с учётом уравнений Максвелла в дифференциальной форме получим [math]\displaystyle{ ^*\mathrm{J}\ }[/math].

С учётом тождества [math]\displaystyle{ \mathrm{d}^2=0\ }[/math] последнее уравнение Максвелла, записанное при помощи дифференциальных форм, сразу приводит к уравнению непрерывности (закону сохранения заряда):

[math]\displaystyle{ \mathrm{d}^*\mathrm{J}=0\ . }[/math]

В такой форме уравнения Максвелла остаются справедливыми и на произвольном 4-мерном многообразии, например, в искривлённом пространстве-времени общей теории относительности. В этом случае в соотношениях дополнительно появляется определитель метрического тензора [math]\displaystyle{ g }[/math]. Например, для тока и внешнего дифференцирования:

[math]\displaystyle{ ^*\mathrm{J} = -\frac{1}{3!}\,j^{\alpha} \sqrt{-g}\,\epsilon_{\alpha\beta\gamma\delta}~ \mathrm{d}x^\beta \wedge \mathrm{d}x^\gamma \wedge \mathrm{d}x^\delta, ~~~~~~~~~~~~~~ \mathrm{d}^*\mathrm{F} = \frac{1}{4}\,F^{\alpha\beta}_{~~;\alpha}\sqrt{-g}\,\epsilon_{\beta\gamma\delta\eta}~ \mathrm{d}x^\gamma \wedge \mathrm{d}x^\delta \wedge \mathrm{d}x^\eta. }[/math]

Общековариантная запись в компонентах

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

В основном рецепт перехода от случая нулевой кривизны пространства-времени и лоренцевых систем отсчёта в нём, подробно описанного выше, к общему случаю состоит в замене обычных производных по координатам на ковариантные производные, учёт того, что метрика в этом случае не постоянна и не имеет специального лоренцева вида (то есть практически произвольна), а также при интегрировании — например, при записи действия — учёт того, что метрика входит в элемент объёма (через множитель [math]\displaystyle{ \sqrt{-g} }[/math] — корень из минус детерминанта метрики).

В общековариантном виде уравнения Максвелла имеют вид:

[math]\displaystyle{ F_{\mu \nu : \sigma} + F_{\nu \sigma : \mu} + F_{\sigma \mu : \nu} = 0 }[/math]

[math]\displaystyle{ F^{\mu \nu}_{: \nu}=4 \pi J^{\mu} }[/math]

Здесь знак «:» означает ковариантную производную, подобно тому, как знак «,» означает обычную производную:

[math]\displaystyle{ A_{\mu : \nu} = A_{\mu , \nu} - \Gamma_{\mu \nu}^{\alpha}A_{\alpha} }[/math],

где [math]\displaystyle{ \Gamma_{\mu \nu}^{\alpha} }[/math] — символ Кристоффеля второго типа.

Закон сохранения электрического заряда в общековариантном виде следует из [math]\displaystyle{ F^{\mu \nu}_{: \nu}=4 \pi J^{\mu} }[/math]. Умножая обе части на [math]\displaystyle{ \sqrt{-g} }[/math] и используя тождество [math]\displaystyle{ F^{\mu \nu}_{: \nu} \sqrt{-g} = ( F^{\mu \nu} \sqrt{-g} )_{, \nu} }[/math] находим [math]\displaystyle{ ( F^{\mu \nu} \sqrt{-g} )_{, \nu}=4 \pi J^{\mu} \sqrt{-g} }[/math].

Отсюда получаем закон сохранения электрического заряда в общековариантном виде:

[math]\displaystyle{ ( J^{\mu} \sqrt{-g} )_{, \mu} = \frac{1}{4 \pi}(F^{\mu \nu}\sqrt{-g})_{, \mu \nu}=0 }[/math].

Спинорная формулировка

Уравнения Максвелла можно записать в спинорной форме:

[math]\displaystyle{ \partial_{\lambda \dot \sigma} f_{\dot \mu}^{\dot \sigma} + \partial_{\dot \mu \sigma} f_{\lambda}^{\sigma} = 2s_{\dot \lambda \mu} }[/math],

[math]\displaystyle{ \partial_{\lambda \dot \sigma} f_{\dot \mu}^{\dot \sigma} - \partial_{\dot \mu \sigma} f_{\lambda}^{\sigma} = 0 }[/math],

где спинор второго ранга [math]\displaystyle{ f }[/math] определяется уравнением [math]\displaystyle{ f_{\lambda \mu}=\frac{1}{2}(\partial_{\lambda \dot \sigma} \varphi_{\mu}^{\dot \sigma}+\partial_{\mu \sigma} \varphi_{\sigma}^{\dot \sigma}) }[/math], [math]\displaystyle{ \varphi_{\mu \nu} }[/math] — четырёхмерный потенциал в форме спинора второго ранга, [math]\displaystyle{ \partial_{\mu \nu} }[/math] — оператор четырёхмерного градиента в спинорной форме, [math]\displaystyle{ s_{\mu \nu} }[/math] — плотность тока в спинорной форме.

Спектральное представление

В электродинамике большое значение имеют гармонические колебания. Такие поля можно представить в виде

[math]\displaystyle{ \mathbf{E}(\mathbf{r},t)=\frac{1}{2}\left[\tilde{\mathbf{E}}(\mathbf{r})e^{-i\omega t}+ \mbox{c. c.}\right]=\Re[\tilde{\mathbf{E}}(\mathbf{r})e^{-i\omega t}], }[/math]

[math]\displaystyle{ \mathbf{H}(\mathbf{r},t)=\frac{1}{2}\left[\tilde{\mathbf{H}}(\mathbf{r})e^{-i\omega t}+\mbox{c. c.}\right]=\Re[\tilde{\mathbf{H}}(\mathbf{r})e^{-i\omega t}]. }[/math]

где [math]\displaystyle{ \omega\ }[/math] — частота колебаний поля. Обозначение «c. c.» означает комплексное сопряжение предыдущего слагаемого. В некоторых работах коэффициент 1/2 в соглашении о гармонических амплитудах не используется, что приводит к соответствующей модификации всех связанных с этим соглашением выражений. В литературе также часто встречается выбор обратного знака в комплексной экспоненте. Рассмотренный здесь вариант согласуется с принятым в квантовой теории в представлении Шрёдингера.

Усреднённые за период плотности энергии электрического и магнитного поля равны соответственно

СГС	СИ
[math]\displaystyle{ \mathbf{S}=\frac{c}{16\pi}\left[{\tilde{\mathbf{E}}^*}\times{\tilde{\mathbf{H}}}+{\tilde{\mathbf{E}}}\times{\tilde{\mathbf{H}}^*}\right], }[/math] [math]\displaystyle{ \langle w_E\rangle = \frac{\tilde{\varepsilon}}{16\pi}\,|\tilde{\mathbf{E}}|^2, }[/math] [math]\displaystyle{ \langle w_H\rangle = \frac{\tilde{\mu}}{16\pi}\,|\tilde{\mathbf{H}}|^2. }[/math]	[math]\displaystyle{ \mathbf{S}=\frac{1}{4}\left[{\mathbf{E}^*}\times{\tilde{\mathbf{H}}}+{\tilde{\mathbf{E}}}\times{\tilde{\mathbf{H}}^*}\right], }[/math] [math]\displaystyle{ \langle w_E\rangle = \frac{\varepsilon_0\tilde{\varepsilon} |{\tilde{\mathbf{E}}}|^2}{4}, }[/math] [math]\displaystyle{ \langle w_H\rangle = \frac{\mu_0\tilde{\mu} |\tilde{\mathbf{H}}|^2}{4}. }[/math]

Используя преобразование Фурье, по гармоническим колебаниям можно разложить поля с произвольной временной зависимостью.

[math]\displaystyle{ \mathbf{E}(r,t)=\int\limits_{-\infty}^{\infty}\!\tilde {\mathbf{E}}(\mathbf{r},\omega)e^{-i\omega t} \frac{d\omega}{2\pi}, }[/math]

[math]\displaystyle{ \mathbf{H}(r,t)=\int\limits_{-\infty}^{\infty}\!\tilde {\mathbf{H}}(\mathbf{r},\omega)e^{-i\omega t} \frac{d\omega}{2\pi}. }[/math]

Переход к спектральным компонентам позволяет сосредоточиться на координатной зависимости полей. Уравнения Максвелла для спектральных компонент в однородных средах при этом принимают вид

СГС	СИ
[math]\displaystyle{ \nabla\cdot\tilde{\mathbf{E}}=\frac{4\pi\tilde {\rho}}{\tilde{\varepsilon}}, }[/math]	[math]\displaystyle{ \nabla\cdot\tilde{\mathbf{E}}= \frac{\tilde{\rho}}{\varepsilon_0\tilde{\varepsilon}}, }[/math]
[math]\displaystyle{ \nabla\cdot\tilde{\mathbf{B}}=0, }[/math]	[math]\displaystyle{ \nabla\cdot\tilde{\mathbf{B}}=0, }[/math]
[math]\displaystyle{ \nabla\times\tilde{\mathbf{E}}=i\frac{\omega}{c}\tilde{\mathbf{B}}, }[/math]	[math]\displaystyle{ \nabla\times\tilde{\mathbf{E}}=i\omega\tilde{\mathbf{B}}, }[/math]
[math]\displaystyle{ \nabla\times\tilde{\mathbf{B}}=-i\tilde{\varepsilon}\tilde{\mu}\frac{\omega}{c}\left[1+i\frac{4\pi \tilde{\sigma}}{\omega\tilde{\varepsilon}}\right]\tilde{\mathbf{E}}. }[/math]	[math]\displaystyle{ \nabla\times\tilde{\mathbf{B}}= -i\tilde{\varepsilon}\tilde{\mu}\frac{\omega}{c^2}\left[1+i\frac{\tilde{\sigma}}{\omega\varepsilon_0\tilde{\varepsilon}}\right]\tilde{\mathbf{E}}. }[/math]

Диэлектрическая и магнитная проницаемости среды в спектральном представлении связаны с восприимчивостями материальных уравнений в интегральном представлении Фурье-преобразованием:

СГС	СИ
[math]\displaystyle{ \tilde{\varepsilon}(\omega)=1+4\pi\int\limits_{-\infty}^{\infty}\!\chi_e(\tau)e^{i\omega\tau}d\tau, }[/math] [math]\displaystyle{ \tilde{\mu}(\omega)=1+4\pi\int\limits_{-\infty}^{\infty}\!\chi_m(\tau)e^{i\omega\tau}d\tau. }[/math]	[math]\displaystyle{ \tilde{\varepsilon}(\omega)=1+\int\limits_{-\infty}^{\infty}\!\chi_e(\tau)e^{i\omega\tau}d\tau, }[/math] [math]\displaystyle{ \tilde{\mu}(\omega)=1+\int\limits_{-\infty}^{\infty}\!\chi_m(\tau)e^{i\omega\tau}d\tau. }[/math]

Уравнения без свободных зарядов и токов

В отсутствие свободных зарядов и токов [math]\displaystyle{ \rho=0 }[/math], [math]\displaystyle{ \mathbf{j}=0 }[/math], в изотропных и однородных средах без дисперсии уравнения Максвелла принимают следующий вид:

СГС	СИ
[math]\displaystyle{ \nabla\cdot\mathbf{E}=0,~~~~~~~~~~~~~~~ \nabla\cdot\mathbf{B}=0, }[/math] [math]\displaystyle{ \nabla\times\mathbf{E}=-\frac{1}{c}\frac{\partial\mathbf{B}}{\partial t},~~~~~~ \nabla\times\mathbf{B}=\frac{\varepsilon\mu}{c}\frac{\partial\mathbf{E}}{\partial t}, }[/math]	[math]\displaystyle{ \nabla\cdot\mathbf{E}=0,~~~~~~~~~~~~~ \nabla\cdot\mathbf{B}=0, }[/math] [math]\displaystyle{ \nabla\times\mathbf{E}=-\frac{\partial\mathbf{B}}{\partial t},~~~~~~ \nabla\times\mathbf{B}=\frac{\varepsilon\mu}{c^2}\frac{\partial\mathbf{E}}{\partial t}. }[/math]

Решениями этих уравнений является напряжённость электрического поля [math]\displaystyle{ \mathbf{E} }[/math] и магнитная индукция [math]\displaystyle{ \mathbf{B} }[/math]. Диэлектрическая [math]\displaystyle{ \varepsilon }[/math] и магнитная [math]\displaystyle{ \mu }[/math] проницаемости определяются свойствами среды. Для вакуума [math]\displaystyle{ \varepsilon=1 }[/math], [math]\displaystyle{ \mu=1 }[/math].

Почтовая марка Никарагуа 1971 года и её оборот. Закон Максвелла (электромагнетизм)

Волновое уравнение

Уравнения Максвелла являются дифференциальными уравнениями первого порядка по координатам и времени. Однако во второй паре в каждое уравнение входят обе неизвестные векторные функции [math]\displaystyle{ \mathbf{E} }[/math] и [math]\displaystyle{ \mathbf{B} }[/math]. При отсутствии зарядов и токов можно перейти к уравнениям второго порядка, каждое из которых зависит только от одного (электрического или магнитного) поля:

[math]\displaystyle{ \Delta\mathbf{E}-\frac{\varepsilon\mu}{c^2}\frac{\partial^2\mathbf E}{\partial t^2}=0, }[/math] [math]\displaystyle{ \Delta\mathbf{B}-\frac{\varepsilon\mu}{c^2}\frac{\partial^2\mathbf B}{\partial t^2}=0. }[/math]

Такие уравнения называются волновыми.

Вывод волнового уравнения

Беря ротор от закона Фарадея и используя закон Ампера—Максвелла, получаем (в системе СИ):

[math]\displaystyle{ \nabla\times[\nabla\times\mathbf E]=\nabla\times\left(-\frac{\partial\mathbf B}{\partial t}\right)=-\frac{\partial}{\partial t}[\nabla\times\mathbf B]=-\frac{\varepsilon\mu}{c^2}\frac{\partial^2\mathbf E}{\partial t^2}. }[/math]

С другой стороны, раскрывая двойное векторное произведение, имеем:

[math]\displaystyle{ \nabla\times[\nabla\times\mathbf E]=\nabla\,(\nabla\cdot\mathbf E)-\Delta\mathbf E = -\Delta\mathbf E, }[/math]

так как дивергенция электрического поля в вакууме равна нулю. Приравнивая эти два выражения, получаем волновое уравнение для электрического поля. Аналогично получается волновое уравнение для магнитного поля.

В лоренцевской калибровке в отсутствие зарядов и токов волновому уравнению удовлетворяют также скалярный и векторный потенциалы:

[math]\displaystyle{ \Delta\varphi-\frac{\varepsilon\mu}{c^2}\frac{\partial^2\varphi}{\partial t^2}=0, ~~~~~~~~~~~~~~~~~~~~~~ \Delta\mathbf{A}-\frac{\varepsilon\mu}{c^2}\frac{\partial^2\mathbf A}{\partial t^2}=0. }[/math]

Величина [math]\displaystyle{ c/\sqrt{\varepsilon\mu} }[/math], входящая в волновые уравнения определяет скорость распространения электромагнитных полей в среде. Её максимальное значение [math]\displaystyle{ c }[/math] достигается в вакууме, когда [math]\displaystyle{ \varepsilon=1 }[/math] и [math]\displaystyle{ \mu=1 }[/math].

Уравнение Гельмгольца

Пусть [math]\displaystyle{ \omega }[/math] — круговая частота гармонического сигнала, и зависимость от времени выбрана в виде [math]\displaystyle{ e^{-i\omega t} }[/math]. При отсутствии электрических зарядов в среде уравнение Гельмгольца принимает вид:

[math]\displaystyle{ \Delta\mathbf E+k^2\mathbf E= 0,~~~~~~~~~~~~~~~~~~~~~~~\Delta\mathbf B+k^2\mathbf B= 0, }[/math]

где [math]\displaystyle{ k^2=\mu\varepsilon\frac{\omega^2}{c^2} }[/math].

Уравнения Максвелла в форме Майорана

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

Уравнения Максвелла в форме Майорана имеют вид:

[math]\displaystyle{ i \frac{\partial \xi}{\partial t} = \vec{s} \vec{p} \xi }[/math], [math]\displaystyle{ \vec{p} \vec{\xi} = 0 }[/math],

[math]\displaystyle{ i \frac{\partial \eta}{\partial t} = \vec{s} \vec{p} \eta }[/math], [math]\displaystyle{ \vec{p} \vec{\eta} = 0 }[/math],

Здесь: [math]\displaystyle{ \vec{\xi} = \vec{E} + i \vec{H} }[/math], [math]\displaystyle{ \vec{\eta} = \vec{E} - i \vec{H} }[/math],

[math]\displaystyle{ \vec{E}, \vec{H} }[/math] - векторы электрического и магнитного поля в уравнениях Максвелла для пустоты (в релятивистской системе единиц):

[math]\displaystyle{ \frac{\partial \vec{H}}{\partial t} = - \operatorname{rot} \vec{E} }[/math], [math]\displaystyle{ \operatorname{div} \vec{H} = 0 }[/math],

[math]\displaystyle{ \frac{\partial \vec{E}}{\partial t} = \operatorname{rot} \vec{H} }[/math], [math]\displaystyle{ \operatorname{div} \vec{E} = 0 }[/math],

[math]\displaystyle{ \vec{p} = ( \frac{1}{i} \frac{\partial}{\partial x_1}, \frac{1}{i} \frac{\partial}{\partial x_2}, \frac{1}{i} \frac{\partial}{\partial x_3}) }[/math] - оператор импульса, [math]\displaystyle{ \vec{s} }[/math] - вектор с матричными компонентами:

[math]\displaystyle{ s_1 = \begin{pmatrix} 0&0&0 \\ 0&0&-i \\ 0&i&0 \end{pmatrix}, }[/math]

[math]\displaystyle{ s_2 = \begin{pmatrix} 0&0&i \\ 0&0&0 \\ -i&0&0 \end{pmatrix}, }[/math]

[math]\displaystyle{ s_3 = \begin{pmatrix} 0&-i&0 \\ i&0&0 \\ 0&0&0 \end{pmatrix}, }[/math]

Некоторые точные решения

Поле движущегося точечного заряда

Если заряд движется с постоянной скоростью [math]\displaystyle{ \mathbf{u} }[/math], то вокруг него возникает магнитное поле [math]\displaystyle{ \mathbf{B} }[/math], а напряжённость электрического перестаёт быть сферически симметричной :

СГС	СИ
[math]\displaystyle{ \mathbf{E}=\frac{Q}{\varepsilon}\, \frac{\mathbf{r}}{r^3}\, \frac{1-u^2/c^2}{(1-[\mathbf{n}\times\mathbf{u}]^2/c^2)^{3/2}} }[/math] [math]\displaystyle{ \mathbf{B}=\frac{\varepsilon\mu}{c}\,[\mathbf{u}\times \mathbf{E}] }[/math]	[math]\displaystyle{ \mathbf{E}=\frac{Q}{4\pi\varepsilon\varepsilon_0}\, \frac{\mathbf{r}}{r^3}\, \frac{1-u^2/c^2}{(1-[\mathbf{n}\times\mathbf{u}]^2/c^2)^{3/2}} }[/math] [math]\displaystyle{ \mathbf{B}=\frac{\varepsilon\mu}{c^2}\,[\mathbf{u}\times \mathbf{E}] }[/math]

Единичный вектор [math]\displaystyle{ \mathbf{n}=\mathbf{r}/r }[/math] направлен от заряда к точке измерения напряжённости поля. [math]\displaystyle{ r }[/math] — модуль вектора [math]\displaystyle{ \mathbf r }[/math]. Если ввести угол [math]\displaystyle{ \theta }[/math] между векторами [math]\displaystyle{ \mathbf{u} }[/math] и [math]\displaystyle{ \mathbf{n} }[/math], то [math]\displaystyle{ [\mathbf{n}\times\mathbf{u}]^2=u^2\sin^2\theta }[/math]. При фиксированном расстоянии от заряда напряжённость электрического поля минимальна в точках, находящихся на линии движения заряда. Максимальное значение достигается в плоскости, проходящей через заряд перпендикулярно его скорости. Магнитная индукция, в силу векторного произведения, перпендикулярна скорости и электрическому полю. Так как заряд движется, в фиксированной точке пространства электрическое и магнитное поля изменяются со временем. Они удовлетворяют уравнениям Максвелла с плотностью заряда и тока, пропорциональных дельта-функции Дирака:

[math]\displaystyle{ \rho(\mathbf{r})=Q\delta(\mathbf{r}-\mathbf{r}_0),~~~~~~~~~~~ \mathbf{j}(\mathbf{r})=\mathbf{u}\,\rho(\mathbf{r}) }[/math]

где [math]\displaystyle{ \mathbf{r}_0 }[/math] — текущее положение заряда.

На пробный заряд [math]\displaystyle{ q }[/math], движущийся в той же системе отсчёта, действует сила Лоренца. Она может быть получена при помощи преобразований Лоренца из закона Кулона и принципа инвариантности заряда. В этом смысле магнитное поле по своей природе является релятивистским эффектом.

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

Плоские электромагнитные волны

• Предположим, что напряжённость электрического поля и магнитная индукция являются произвольными функциями следующей комбинации координат и времени:

[math]\displaystyle{ \mathbf{E}(\mathbf r, t)=\mathbf{E}\Bigl(\mathbf{n}\mathbf{r}-\frac{c}{\sqrt{\varepsilon\mu}} t\Bigr), ~~~~~~~~~~~~~~~ \mathbf{B}(\mathbf r, t)=\mathbf{B}\Bigl(\mathbf{n}\mathbf{r}-\frac{c}{\sqrt{\varepsilon\mu}} t\Bigr), }[/math]

где [math]\displaystyle{ \mathbf{n} }[/math] — некоторый постоянный вектор. В этом случае [math]\displaystyle{ \mathbf{E} }[/math] и [math]\displaystyle{ \mathbf{B} }[/math] удовлетворяют уравнениям Максвелла в отсутствие зарядов и токов, если между ними существует следующая связь:

СГС	СИ
[math]\displaystyle{ \mathbf{B}=\sqrt{\varepsilon\mu}\,[\mathbf{n}\times\mathbf{E}] }[/math]	[math]\displaystyle{ \mathbf{B}=\frac{\sqrt{\varepsilon\mu}}{c}\,[\mathbf{n}\times\mathbf{E}] }[/math]

и они перпендикулярны вектору [math]\displaystyle{ \mathbf{n} }[/math], который должен быть единичным:

[math]\displaystyle{ \mathbf{n}\mathbf{E}=0,~~~~~~~~~~~~~~~\mathbf{n}\mathbf{B}=0,~~~~~~~~~~~~~~~\mathbf{n}^2=1. }[/math]

Вывод решения для плоской волны

Если напряжённость электрического поля зависит от координат и времени в виде следующей их комбинации [math]\displaystyle{ u=\mathbf{n}\mathbf{r}-c t/\sqrt{\varepsilon\mu} }[/math], то для производной [math]\displaystyle{ j\ }[/math]-той компоненты вектора [math]\displaystyle{ \mathbf{E} }[/math] по [math]\displaystyle{ i\ }[/math]-той координате и времени можно записать:

[math]\displaystyle{ \nabla_i E_j = \frac{\partial E_j}{\partial r_i} = \frac{d E_j}{du}\,\frac{\partial u}{\partial r_i} = n_i E_j ,~~~~~~~~~~~~~~~~ \frac{\partial E_j}{\partial t} = -\frac{c}{\sqrt{\varepsilon\mu}} \frac{d E_j}{du}, }[/math]

и аналогично для магнитной индукции. Поэтому уравнения Максвелла в отсутствие зарядов и токов принимают вид (система СИ):

[math]\displaystyle{ \frac{d(\mathbf{n}\mathbf{E})}{du} = 0,~~~~~~~~ \frac{d(\mathbf{n}\mathbf{B})}{du} = 0,~~~~~~~~ \frac{d}{du}\bigl([\mathbf{n}\times \mathbf{E}] - \frac{c}{\sqrt{\varepsilon\mu}}\mathbf{B}\bigr)= 0,~~~~~~~~ \frac{d}{du}\bigl([\mathbf{n}\times \mathbf{B}] + \frac{\sqrt{\varepsilon\mu}}{c} \mathbf{E}\bigr)= 0. }[/math]

Интегрируя эти соотношения по [math]\displaystyle{ u }[/math] и опуская константы интегрирования, которые соответствуют постоянным полям, получаем:

[math]\displaystyle{ \mathbf{n}\mathbf{E} = 0,~~~~~~~~~ \mathbf{n}\mathbf{B} = 0,~~~~~~~~~ \mathbf{B} = \frac{\sqrt{\varepsilon\mu}}{c}\,[\mathbf{n}\times \mathbf{E}],~~~~~~~~~ \mathbf{E}= -\frac{c}{\sqrt{\varepsilon\mu}}\, [\mathbf{n}\times \mathbf{B}]. }[/math]

Подставляя четвёртое уравнение в третье, получаем [math]\displaystyle{ \mathbf{n}^2=1 }[/math].

Физический смысл решения в виде плоской волны состоит в следующем. Выберем ось [math]\displaystyle{ z\ }[/math] декартовой системы координат так, чтобы вектор [math]\displaystyle{ \mathbf{n} }[/math] был направлен вдоль неё. Тогда электрические и магнитные поля волны зависят от координаты [math]\displaystyle{ z\ }[/math] и времени [math]\displaystyle{ t\ }[/math] следующим образом:

[math]\displaystyle{ \mathbf{E}(z, t)=\mathbf{E}\Bigl(z-\frac{c}{\sqrt{\varepsilon\mu}} t\Bigr), ~~~~~~~~~~~~~~~ \mathbf{B}(z, t)=\mathbf{B}\Bigl(z-\frac{c}{\sqrt{\varepsilon\mu}} t\Bigr), }[/math]

Предположим, что в начальный момент времени [math]\displaystyle{ t=0\ }[/math] напряжённость поля является произвольной векторной функцией [math]\displaystyle{ z\ }[/math]. С течением времени эта функция будет сдвигаться в пространстве вдоль оси [math]\displaystyle{ z\ }[/math] со скоростью [math]\displaystyle{ c/\sqrt{\varepsilon\mu} }[/math].

В электромагнитной волне в общем случае напряжённость поля может быть произвольной непериодической функцией [math]\displaystyle{ u=\mathbf{n}\mathbf{r}-c t/\sqrt{\varepsilon\mu} }[/math]. Например, решение в виде плоской волны может описывать электромагнитный импульс, локализованный вдоль направления движения. В плоскости, перпендикулярной [math]\displaystyle{ \mathbf{n} }[/math], электромагнитные поля не изменяются, что означает, что в этой плоскости плоская волна не ограничена и имеет плоский фазовый фронт (именно поэтому волна называется плоской). Так как электрическое и магнитное поля при распространении плоской волны всё время остаются перпендикулярными вектору [math]\displaystyle{ \mathbf{n} }[/math], такие волны называют «поперечными», или «трансверсальными». Векторы [math]\displaystyle{ \mathbf{E} }[/math] и [math]\displaystyle{ \mathbf{B} }[/math] в силу свойств векторного произведения также перпендикулярны друг другу.

[math]\displaystyle{ \bullet }[/math] Плотности энергии электрического и магнитного поля в плоской волне равны друг другу:

СГС	СИ
[math]\displaystyle{ w_E=w_H = \frac{\varepsilon}{8\pi}\, \mathbf{E}^2 }[/math]	[math]\displaystyle{ w_E=w_H = \frac{\varepsilon\varepsilon_0}{2} \, \mathbf{E}^2 }[/math]

Вектор Пойнтинга (плотность потока энергии) независимо от системы единиц связан с полной плотностью энергии следующим образом:

[math]\displaystyle{ \mathbf{S}=\frac{c}{\sqrt{\varepsilon\mu}}\, \mathbf{n}\, (w_E+w_H). }[/math]

Это соотношение соответствует уравнению связи импульса и энергии для безмассовой частицы в релятивистской теории. Однако скорость [math]\displaystyle{ c/\sqrt{\varepsilon\mu} }[/math] в среде меньше, чем скорость света в вакууме [math]\displaystyle{ c }[/math].

Плоская электромагнитная волна с круговой (сверху) и линейной (снизу) поляризацией

Плоские и поперечные волны являются математическими абстракциями. Реальные волны конечной апертуры из-за эффекта дифракции можно считать плоскими и поперечными лишь в некотором приближении.

[math]\displaystyle{ \bullet }[/math] Важный частный случай решения в виде плоских волн возникает, когда напряжённости полей являются гармоническими периодическими функциями. Выберем координатную ось [math]\displaystyle{ z }[/math] вдоль волнового вектора [math]\displaystyle{ \mathbf{k} }[/math]. Тогда вектор электрического поля (как, впрочем, и магнитного) будет лежать в плоскости [math]\displaystyle{ (x,y) }[/math], то есть [math]\displaystyle{ \mathbf{E}=(E_x,E_y,0) }[/math]. Если по каждой проекции в этой плоскости электрическое поле совершает периодические колебания, то такую волну называют монохроматической плоской волной:

[math]\displaystyle{ E_x=a\, \cos(\omega t-\mathbf{k}\mathbf{r}+\phi_1),~~~~~~~~~~~~~~~~~~~ E_y=b\, \cos(\omega t-\mathbf{k}\mathbf{r}+\phi_2). }[/math]

Сравнение с общим решением для плоской волны приводит к следующей связи между вектором [math]\displaystyle{ \mathbf{k} }[/math] и константой [math]\displaystyle{ \omega }[/math], которое называется уравнением дисперсии:

[math]\displaystyle{ \mathbf{k}= \mathbf{n}\, \frac{\sqrt{\varepsilon\mu}}{c}\,\omega,~~~~~~~~~~~~~\mathbf{k}^2 = \frac{\varepsilon\mu}{c^2}\,\omega^2. }[/math]

В этом случае, вектор [math]\displaystyle{ \mathbf{k} }[/math] называется волновым вектором, а [math]\displaystyle{ \omega }[/math] — круговой частотой монохроматической электромагнитной волны. Модуль волнового вектора и круговая частота связаны с длиной волны [math]\displaystyle{ \lambda }[/math] и её частотой [math]\displaystyle{ \nu }[/math] следующим образом:

[math]\displaystyle{ k = \frac{2\pi}{\lambda},~~~~~~~~~\omega=2\pi\nu,~~~~~~~~~~~\lambda\nu = \frac{c}{\sqrt{\varepsilon\mu}}. }[/math]

Константы [math]\displaystyle{ \phi_1 }[/math] и [math]\displaystyle{ \phi_2 }[/math] являются сдвигами фазы, а [math]\displaystyle{ a }[/math] и [math]\displaystyle{ b }[/math] — амплитудами колебаний вдоль каждой оси.

В фиксированной точке пространства ([math]\displaystyle{ \mathbf{k}\mathbf{r}=\mathrm{const} }[/math]) вектор электрического поля, в общем случае, описывает в плоскости [math]\displaystyle{ (x,y)\ }[/math] эллипс, поэтому такие волны называются эллиптически поляризованными. Их частным случаем являются волны, поляризованные по кругу. Вырожденный в прямую эллипс соответствует колебаниям напряжённости поля вдоль одной прямой в плоскости [math]\displaystyle{ (x,y)\ }[/math]. Такие волны называются линейно поляризованными. Аналогична ситуация с вектором магнитной индукции, который всё время перпендикулярен напряжённости электрического поля.

Связь с другими теориями

Уравнения Максвелла полностью совместимы с принципами специальной теории относительности. Они также применимы при микроскопическом описании вещества, когда заряженные частицы подчиняются принципам квантовой механики, а электромагнитное поле остаётся классическим (не квантовым). В этом случае квантовые объекты (например, электроны) описываются уравнением Шрёдингера или уравнением Дирака, однако потенциалы электромагнитного взаимодействия в этих уравнениях определяются классическими уравнениями Максвелла.

Тем не менее существуют явления, для описания которых требуется более последовательное объединение полевого подхода Фарадея — Максвелла с принципами квантовой механики. Оно осуществляется при помощи методов квантовой теории поля в квантовой электродинамике. В этом случае форма уравнений Максвелла (лагранжиан) остаётся неизменной, однако поля становятся операторами, а уравнения Максвелла — операторными уравнениями Гейзенберга. Решение подобных уравнений приводит к появлению новых эффектов, отсутствующих в классической теории поля. Эти эффекты существенны, в частности, в следующих физических ситуациях:

• Сверхсильные поля ([math]\displaystyle{ E \sim m_e^2c^3/e\hbar }[/math]≈1.32×10¹⁸ В/м, где [math]\displaystyle{ m_e\ }[/math] — масса электрона, [math]\displaystyle{ e\ }[/math] — его заряд, [math]\displaystyle{ \hbar }[/math] — постоянная Планка) — работа такого поля на комптоновской длине волны электрона равна по порядку величины энергии покоя электрона, что приводит к самопроизвольной генерации электрон-позитронных пар из вакуума (эффект Швингера). В результате возникает эффективное взаимодействие фотонов, которое отсутствует в классической электродинамике, приводящее к эффективному изменению лагранжиана поля (например, в низкоэнергетическом пределе поле описывается лагранжианом Гейзенберга — Эйлера).
• Сверхслабые поля с энергией [math]\displaystyle{ \sim\hbar\omega }[/math], где [math]\displaystyle{ \omega\ }[/math] — частота поля (см. формула Планка). В этом случае становятся заметными отдельные кванты электромагнитного поля — фотоны.
• Для описания эффектов поглощения и испускания света атомами и молекулами.
• Для описания неклассических, например, сжатых состояний поля.
• На малых расстояниях, сравнимых с комптоновской длиной волны электрона [math]\displaystyle{ \lambda_e=\hbar/m_ec }[/math] ≈ 3,86×10⁻¹³ м, когда в результате вакуумных эффектов модифицируется, например, закон Кулона.

Аксиоматический подход

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

• Постулируются:
  • закон Кулона (сила [math]\displaystyle{ \mathbf{F} }[/math], действующая на пробный заряд [math]\displaystyle{ q }[/math] со стороны неподвижного заряда [math]\displaystyle{ Q }[/math]);
  • инвариантность заряда в различных инерциальных системах отсчёта;
  • принцип суперпозиции.
• При помощи преобразований Лоренца получается значение для вектора силы [math]\displaystyle{ \mathbf{F} }[/math], действующей на пробный заряд со стороны равномерно движущегося со скоростью [math]\displaystyle{ \mathbf{u} }[/math] заряда [math]\displaystyle{ Q }[/math], которое совпадает с силой Лоренца.
• Дивергенция и ротор, вычисленные от электрической ([math]\displaystyle{ \mathbf{E} }[/math]) и магнитной ([math]\displaystyle{ \mathbf{B} }[/math]) составляющих силы, дают уравнения Максвелла для точечного заряда. В силу принципа суперпозиции они записываются для произвольного распределения зарядов и токов. В заключение постулируется применимость этих уравнений и к ускоренному движению зарядов.

Второй подход основан на лагранжевом формализме . При этом постулируется, что электромагнитное поле описывается линейным взаимодействием четырёхмерного потенциала [math]\displaystyle{ A^\alpha }[/math] с четырёхвектором электрического тока [math]\displaystyle{ j^\alpha }[/math], а свободный лагранжиан пропорционален инвариантной свёртке квадрата тензора электромагнитного поля [math]\displaystyle{ F^{\alpha\beta} }[/math].

Как в первом, так и во втором подходе предполагаются установленными принципы теории относительности. Хотя исторически она возникла на основе уравнений Максвелла и второго постулата Эйнштейна, известен восходящий к работам Игнатовского , Франка и Роте аксиоматический способ построения СТО, не использующий постулата об инвариантности скорости света и уравнений Максвелла.

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

Единственность решений уравнений Максвелла

Уравнения Максвелла являются дифференциальными уравнениями в частных производных. Поэтому для их решения необходимо задать начальные и граничные условия. При фиксированных функциях плотности заряда и тока для нестационарных полей получаемое решение единственно. Этот факт формулируется в виде теоремы:

Если напряженности электрического и магнитного полей заданы в начальный момент времени [math]\displaystyle{ t=0 }[/math] в каждой точке некоторой области пространства [math]\displaystyle{ v }[/math], и в течение всего времени [math]\displaystyle{ t\ge 0 }[/math] заданы тангенциальные (касательные) составляющие напряженности электрического или магнитного поля на границе этой области [math]\displaystyle{ s }[/math], то существует единственное решение уравнений Максвелла.

Доказательство

Пусть электрическая и магнитная индукции связаны с напряжённостями полей при помощи следующих материальных уравнений:

[math]\displaystyle{ \mathbf{D}(\mathbf{r}, t)=\hat{\varepsilon}(\mathbf{r})\cdot\mathbf{E}(\mathbf{r}, t), ~~~~~~~~~~~~~~~~~ \mathbf{B}(\mathbf{r}, t)=\hat{\mu}(\mathbf{r})\cdot\mathbf{H}(\mathbf{r}, t), }[/math]

где [math]\displaystyle{ \hat{\varepsilon}(\mathbf{r}) }[/math] и [math]\displaystyle{ \hat{\mu}(\mathbf{r}) }[/math] - положительно определённые, симметричные, стационарные матрицы. Если при данных начальных и граничных условиях существуют два различных решения, то следующие величины будут отличны от нуля:

[math]\displaystyle{ \mathbf{e}(\mathbf{r}, t)=\mathbf{E}_2(\mathbf{r}, t)-\mathbf{E}_1(\mathbf{r}, t), ~~~~~~~~~~~~~~~~~ \mathbf{h}(\mathbf{r}, t)=\mathbf{H}_2(\mathbf{r}, t)-\mathbf{H}_1(\mathbf{r}, t), }[/math]

где индекс указывает номер решения. Так как начальные и граничные условия заданы (одинаковые для обоих возможных решений), то:

[math]\displaystyle{ \mathbf{e}(\mathbf{r}, 0)=\mathbf{h}(\mathbf{r}, 0)=0,~~~~~~~~~~~~~~~ [\mathbf{e}_\tau(\mathbf{r}, t) \cup \mathbf{h}_\tau(\mathbf{r}, t)]\Big|_s=0. }[/math]

Первые соотношения соответствуют начальным условиям, а вторые - граничным условиям на поверхности [math]\displaystyle{ s }[/math], где [math]\displaystyle{ \mathbf{e}=\mathbf{e}_n+\mathbf{e}_\tau }[/math]. (Индекс [math]\displaystyle{ n }[/math] - нормальная составляющая к поверхности, а [math]\displaystyle{ \tau }[/math] - касательная. Аналогично для [math]\displaystyle{ \mathbf{h} }[/math]) Подстановка функций [math]\displaystyle{ \mathbf{e}(\mathbf{r}, t) }[/math] и [math]\displaystyle{ \mathbf{h}(\mathbf{r}, t) }[/math] в уравнения Максвелла для роторов приводит к следующим уравнениям:

[math]\displaystyle{ k\,\nabla\times\mathbf{e}=-\frac{\partial (\hat{\mu}\cdot\mathbf{h})}{\partial t}, ~~~~~~~~~~~~~~~~~ k\,\nabla\times\mathbf{h}=\frac{\partial (\hat{\varepsilon}\cdot\mathbf{e})}{\partial t}, }[/math]

где коэффициент [math]\displaystyle{ k }[/math] равен [math]\displaystyle{ c }[/math] в системе СГС и единице в системе СИ. Если одно из разностных полей [math]\displaystyle{ \mathbf{e} }[/math] или [math]\displaystyle{ \mathbf{h} }[/math] равно нулю, то в силу нулевых начальных условий, соответственно, из первого или второго уравнения следует равенство нулю неопределенного разностного поля, соответственно, [math]\displaystyle{ \mathbf{h} }[/math] или [math]\displaystyle{ \mathbf{e} }[/math], и единственность в этих частных случаях доказана.

Предположим, что не равны нулю оба разностных поля. Если первое уравнение умножить на [math]\displaystyle{ \mathbf{h} }[/math], а второе на [math]\displaystyle{ \mathbf{e} }[/math], и вычесть друг из друга, то получится следующее выражение:

[math]\displaystyle{ -2k\,\nabla\cdot(\mathbf{e}\times\mathbf{h})= \frac{\partial (\mathbf{e}\cdot\hat{\varepsilon}\cdot\mathbf{e})}{\partial t} +\frac{\partial (\mathbf{h}\cdot\hat{\mu}\cdot\mathbf{h})}{\partial t}. }[/math]

Это выражение можно проинтегрировать по объёму [math]\displaystyle{ v }[/math], и применить теорему Гаусса:

[math]\displaystyle{ -2k\oint\limits_s [\mathbf{e}\times\mathbf{h}]\,d\mathbf{s} = \frac{d}{d t} \int\limits_v (\mathbf{e}\cdot\hat{\varepsilon}\cdot\mathbf{e}+\mathbf{h}\cdot\hat{\mu}\cdot\mathbf{h})\, dv. }[/math]

Тангенциальные (касательные) к поверхности [math]\displaystyle{ s }[/math] компоненты векторов [math]\displaystyle{ \mathbf{e}_\tau }[/math] или [math]\displaystyle{ \mathbf{h}_\tau }[/math] при любом [math]\displaystyle{ t\ge 0 }[/math] равны нулю (граничные условия), поэтому равен нулю и интеграл по поверхности. Следовательно:

[math]\displaystyle{ \frac{d}{d t} \int\limits_v (\mathbf{e}\cdot\hat{\varepsilon}\cdot\mathbf{e}+\mathbf{h}\cdot\hat{\mu}\cdot\mathbf{h})\, dv = 0. }[/math]

Полученное соотношение интегрируется по времени. Так как в начальный момент времени [math]\displaystyle{ t=0 }[/math] функции [math]\displaystyle{ \mathbf{e}(\mathbf{r}, 0)=\mathbf{h}(\mathbf{r}, 0)=0 }[/math], то константа интегрирования равна нулю, и при любом [math]\displaystyle{ t }[/math]:

[math]\displaystyle{ \int\limits_v (\mathbf{e}\cdot\hat{\varepsilon}\cdot\mathbf{e}+\mathbf{h}\cdot\hat{\mu}\cdot\mathbf{h})\,dv = 0. }[/math]

Подынтегральная функция является положительно определённой (всегда больше или равна нулю). Интеграл от такой функции равен нулю только в том случае, когда подынтегральная функция тождественно равна нулю. Следовательно, в любой момент времени внутри объёма [math]\displaystyle{ \mathbf{e}(\mathbf{r}, t)=0 }[/math] и [math]\displaystyle{ \mathbf{h}(\mathbf{r}, t)=0 }[/math]. Поэтому решения совпадают.

Для единственности решения уравнений Максвелла вместо задания тангенциальных компонент поля можно потребовать выполнение на границе условия импедансного типа

[math]\displaystyle{ \left(\mathbf{E}_\tau - Z [\mathbf{n}\times \mathbf{H}]\right)\Big|_s=0, }[/math]

где импеданс [math]\displaystyle{ Z }[/math] выбран так, чтобы исключить приток энергии извне. Такое условие позволяет сформулировать теорему единственности и в неограниченном случае, при этом импедансное условие обращается в условие излучения Зоммерфельда на бесконечности.

Для гармонических во времени процессов единственность решения задачи без начальных условий обеспечивается сколь угодно малым поглощением энергии внутри объёма [math]\displaystyle{ v }[/math] или её утечкой через поверхность [math]\displaystyle{ s }[/math] (исключающими собственные колебания на действительных резонансных частотах).

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

Численное решение уравнений Максвелла

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

Основными методами являются проекционные, в которых решение проецируется на какой-либо удобный функциональный базис, и дискретизационные — область пространства разбивается на множество малых конечных областей.

• В проекционном методе Бубнова — Галёркина решение граничной задачи рассматривается в виде приближённого конечного разложения по базисным функциям. После подстановки разложения в исходные уравнения с учётом требования ортогональности получающейся невязки выбранным базисным функциям получается система линейных уравнений для коэффициентов разложения.

Для компьютерных расчётов чаще применяются более универсальные дискретизационные методы:

• Метод конечных элементов (FEM), который используется для решения широкого класса задач, сводящихся к уравнениям в частных производных. В теории электромагнетизма чаще используется для расчёта задач электростатики, магнитостатики, распространения волн и квазистационарных явлений. В методе конечных элементов рассматриваемая область пространства, в которой ищется решение, разбивается на большое число простых дискретных элементов, обычно, но не обязательно, треугольной (в двумерном случае) или тетраэдральной формы (в трёхмерном случае). Форма и плотность элементов адаптируются к требованиям задачи. Поведение отдельных элементов рассматривается как результат линейного взаимодействия соседних узлов решётки разбиения под действием внешних сил и описывается матричными уравнениями. Решение задачи сводится, таким образом, к решению разреженных систем большого числа линейных матричных уравнений. Метод реализован во многих коммерческих и свободных программных пакетах (см. статью Метод конечных элементов).

• Метод конечных разностей во временной области (FDTD) для нахождения временны́х и спектральных зависимостей был разработан специально для решения уравнений Максвелла, в которых изменение электрического и магнитного поля во времени зависит от изменения, соответственно, магнитного и электрического поля в пространстве. В рамках этого метода область пространства и временной интервал подвергаются равномерной дискретизации с заданием начальных условий. Полученные из уравнений Максвелла конечно-разностные уравнения решаются в каждый последующий момент временной сетки, пока не будет получено решение поставленной задачи на всем требуемом временном интервале.

См. также

• Электродинамика
• Квантовая электродинамика
• Максвелл, Джеймс Клерк
• Теория Янга — Миллса
• Потенциалы Лиенара — Вихерта
• Уравнения Ефименко
• Закон электромагнитной индукции

Рецензия

• Рецензия на данную статью Руниверсалис в версии от 17 апреля 2023 года, написанная к.ф.-м.н., доцентом кафедры теоретической физики МФТИ Ивановым Михаилом Геннадьевичем.

Литература

Исторические публикации

• Ампер А. М. Электродинамика. — М.: АН СССР, 1954. — 492 с. — 5000 экз.
• Максвелл Дж. К. Избранные сочинения по теории электромагнитного поля. — М.: ГИТТЛ, 1952. — 687 с.
• Максвелл Дж. К. Статьи и речи. — М.: Наука, 1968. — 423 с.
• Фарадей М. Экспериментальные исследования по электричеству. — М.: АН СССР, 1947—1959. — Т. I—III.
• Maxwell J. C. A dynamical theory of the electromagnetic field // Philosophical Transactions of the Royal Society of London. — 1865. — Т. 155. — С. 459—512.
• Maxwell J. C., A Treatise on Electricity And Magnetism — Volume 1 — 1873 — Posner Memorial Collection — Carnegie Mellon University
• Maxwell J. C., A Treatise on Electricity And Magnetism — Volume 2 — 1873 — Posner Memorial Collection — Carnegie Mellon University
• Максвелл Дж. К. Трактат об электричестве и магнетизме. — М.: Наука, 1989. — Т. 1.
• Максвелл Дж. К. Трактат об электричестве и магнетизме. — М.: Наука, 1989. — Т. 2.
• Maxwell J. C. On Physical Lines of Force // Philosophical Magazine : Ser. 4. — 1861,1862. — Т. 11,13. — С. 161—175, 281—291, 338—347; 12—23, 85—95.
• Maxwell J. C. On Faraday's Lines of Force // Transactions of the Cambridge Philosophical Society. — 1856. — Т. 10, № 1. — С. 155—229.

История развития

• Болотовский Б. М. Оливер Хевисайд. — М.: Наука, 1985. — 260 с.
• Карцев В. П. Максвелл (Жизнь замечательных людей). — М.: Молодая гвардия, 1976. — 336 с.
• Карцев В. П. Приключения великих уравнений. — М.: Знание, 1986. — 288 с.
• Шапиро И. С. К истории открытия уравнений Максвелла (рус.) // Успехи физических наук. — Российская академия наук, 1972. — Т. 108, № 2. — С. 319—333.
• Купер Л. Физика для всех. т.1 = Leon N. Cooper An Introduction to the Meaning and Structure of Phisics, 1968,1970. — М.: Мир, 1973.

Общие курсы физики

• Астахов А. В., Широков Ю. М. Курс физики, Т. II, Электромагнитное поле. — М.: Наука, 1980. — 360 с.
• Баскаков С. И. Основы электродинамики. — М.: Советское радио, 1973. — 248 с.
• Калашников С. Г. Электричество (Общий курс физики, т. 2). — 6-е изд. — М.: Физматлит, 2003. — 624 с. — ISBN 5-9221-0312-1.
• Матвеев А. Н. Электричество и магнетизм. — М.: Высшая школа, 1983.
• Никольский В. В., Никольская Т. И. Электродинамика и распространение радиоволн. М.: Наука, 1989.
• Пеннер Д. И., Угаров В. А. Электродинамика и специальная теория относительности. М.: Просвещение, 1980.
• Парселл Э. Электричество и магнетизм. Берклеевский курс физики. Том 2. М.: Наука, 1971.
• Савельев И. В. Курс общей физики, том II. Электричество. — Наука, 1970.
• Сивухин Д. В. Общий курс физики. — Изд. 3-е, испр. и доп.. — М.: Наука, 1996. — Т. III. Электричество. — 320 с. — ISBN 5-02-014054-6..
• Тоннела М. А. Основы электромагнетизма и теории относительности. М.: ИЛ, 1962.
• Фейнман Р., Лейтон Р., Сэндс М. Фейнмановские лекции по физике. Том 5. Электричество и магнетизм. М.: Мир, 1965
• Фейнман Р., Лейтон Р., Сэндс М. Фейнмановские лекции по физике. Том 6. Электродинамика. М.: Мир, 1966

Курсы теоретической физики

• Баранов А. М., Овчинников С. Г., Золотов О. А., Паклин Н. Н., Титов Л. С. Теоретическая физика: Электродинамика. Электродинамика сплошных сред. Учебное пособие по курсу «Электродинамика и основы электродинамики сплошных сред». — Красноярск: СФУ, 2008. — 198 с.
• Джексон Дж. Классическая электродинамика. — М.: Мир, 1965.
• Ландау Л. Д., Лифшиц Е. М. Теория поля (Теоретическая физика, т. II). — М.: Физматлит, 2003. — 536 с. — ISBN 5-9221-0056-4.
• Ландау Л. Д., Лифшиц Е. М. Электродинамика сплошных сред (Теоретическая физика, т. VIII). — М.: Физматлит, 2005. — 656 с. — ISBN 5-9221-0123-4.
• Батыгин В. В., Топтыгин И. Н. Современная электродинамика. Часть 1. Микроскопическая теория. — М.: НИЦ «Регулярная и хаотическая динамика», 2005. — 736 с. — ISBN 5-93972-492-2.
• Топтыгин И. Н. Современная электродинамика. Часть 2. Теория электромагнитных явлений в веществе. — М.: НИЦ «Регулярная и хаотическая динамика», 2005. — 848 с. — ISBN 5-93972-493-0.
• Стрэттон Дж. А. Теория электромагнетизма. — М.—Л.: ГИТТЛ, 1948. — 539 с. — 8000 экз.
• Тамм И. Е. Основы теории электричества. — М.: Наука, 1989.
• Фущич В. И., Никитин А. Г. Симметрия уравнений Максвелла. — Киев: Наукова думка, 1983. — C. 200.
• Болибрух А. А. Уравнения Максвелла и дифференциальные формы. — МЦНМО, 2002.

Решения уравнений Максвелла

• Баландин М. Ю., Шурина Э. П. Векторный метод конечных элементов: Учебное пособие. — Новосибирск: НГТУ, 2001. — 69 с.
• Вайнштейн Л. А. Электромагнитные волны. — М.: Радио и связь, 1988.
• Вонсовский С. В. Магнетизм. Магнитные свойства диа-, пара-, ферро-, антиферро-, и ферримагнетиков. — М.: Наука, 1971.
• Гинзбург В. Л. Распространение электромагнитных волн в плазме. — 2-е изд. — М.: Наука, 1967.
• Сильвестер П. и Феррари Р. Метод конечных элементов для радиоинженеров и инженеров-электриков. — М.: Мир, 1986. — 336 с.

Ссылки

• Уравнения Максвелла (Научный портал «Элементы большой науки»)
• Основные физические понятия электродинамики