Использование геометризации уравнений Максвелла при расчёте оптических приборов
- Авторы: Кулябов Д.С.1
-
Учреждения:
- Российский университет дружбы народов
- Выпуск: Том 25, № 1 (2017)
- Страницы: 81-90
- Раздел: Математическое моделирование
- URL: https://journals.rudn.ru/miph/article/view/15165
- DOI: https://doi.org/10.22363/2312-9735-2017-25-1-81-90
Цитировать
Полный текст
Аннотация
Развитие физики в XX-м веке было тесно связано с развитием математического аппарата. Общая теория относительности продемонстрировала силу геометрического подхода. К сожалению проникновение этого аппарата в другие области физики происходит достаточно медленно. Например, было несколько попыток внедрения геометрических методов в электродинамику, однако до последнего времени они оставались лишь теоретическими упражнениями. Интерес к геометрическим методам в электродинамике вызван практической необходимостью. Представляется заманчивым следующий алгоритм конструирования электромагнитного прибора. Строятся предполагаемые траектории распространения электромагнитных волн. Затем по этим траекториям вычисляются параметры среды. Также представляет интерес и обратная задача. В работе рассматривается методика расчёта оптических приборов на основе метода геометризации уравнений Максвелла. В основе метода лежит представление материальных уравнений Максвелла в виде эффективной геометрии пространства-времени. Таким образом мы получаем задачу, сходную с некой биметрической теорией гравитации, что позволяет применять хорошо разработанный аппарат дифференциальной геометрии. На основании этого мы можем как исследовать распространение электромагнитного поля по заданным параметрам среды, так и находить параметры среды по заданному закону распространения электромагнитного поля.
Полный текст
1. Введение Аппарат дифференциальной геометрии являлся основным языком физики XX-го века. Его базовые элементы развивались в рамках общей теории относительности. Возникает желание применить этот развитый и к другим областям физики, в частности к оптике. Первые попытки применения методов дифференциальной геометрии в электродинамике следует отнести к публикациям И. Е. Тамма [1-3]. В 1960 году Е. Плебаньский предложил метод геометризации материальных уравнений электромагнитного поля [4-7], ставший классическим. Все последующие работы либо использовали его, либо пытались немного подправить, не меняя идеологии [8]. К сожалению, в статье Плебаньского [4] нет никакого вывода формул, а идеология вывода также не выражена явно. Кроме того, методика Плебаньского выглядит скорее как хитрый трюк. Автор постарался выполнить геометризацию уравнений Максвелла более формально. 2. Обозначения и соглашения • Будем использовать нотацию абстрактных индексов [9]. В данной нотации тензор как целостный объект обозначается просто индексом (например, xi), компоненты обозначаются подчёркнутым индексом (например, xi). • Будем придерживаться следующих соглашений. Греческие индексы (a, ) будут относиться к четырёхмерному пространству и в компонентном виде будут иметь следующие значения: a = 0, 3. Латинские индексы из середины алфавита (i, , k) будут относиться к трёхмерному пространству и в компонентном виде будут иметь следующие значения: i = 1, 3. • Запятой в индексе обозначается частная производная по соответствующей координате (,i := ∂i ); точкой с запятой - ковариантная производная (;i := i). • Для записи уравнений электродинамики в работе используется система СГС симметричная [10]. 3. Представления уравнений Максвелла Будем использовать запись уравнений Максвелла в криволинейных координатах. Более подробное описание дано в статьях [11-15]. Уравнения Максвелла в 3-х мерной форме имеют вид: где eik - альтернирующий тензор. Запишем уравнение Максвелла через тензоры электромагнитного поля Fa и Ga [16-18]: где тензоры Fa и Ga имеют следующие компоненты Здесь Ei, i - компоненты векторов напряжённости электрического и магнитного полей соответственно; Di, i - компоненты векторов электрической и магнитной индукции соответственно. Запишем уравнения (2) и (3) через дифференциальные формы в формализме расслоенных пространств. Будем рассматривать расслоение (5) где = M 4 - четырёхмерное пространство. При этом мы не делаем предположение о метрике данного пространства. Зададим F (2-форма), G (бивектор) и (вектор): Тогда уравнения (2) и (3) примут вид: Здесь d = ♯-1 d ♯ - дивергенция, ♯ : Λk Λn-k задаёт двойственность Пуанкаре. При этом Fa и Ga имеют смысл кривизны в кокасательном ( *) и касательном ( ) расслоениях. Связь между этими величинами задаётся следующим образом: (9) Тогда уравнение (8) примет вид: В линейном случае соотношение (9) можно задать как В этом случае уравнение (8) примет вид: Будем считать, что l задаёт на некоторую эффективную метрику. Тензор проницаемостей Будем считать, что отображение l : Λ2 Λ2 линейное и локальное. Тогда его можно представить в следующем виде: где lagd - тензор проницаемостей, содержащий информацию как об диэлектрической и магнитной проницаемостях, так и об электромагнитной связи [1, 3]. Из (13) видно, что lagd имеет следующую симметрию: Для уточнения симметрии, тензор lagd можно представить в следующем виде [19-22]: Очевидно, что lagd имеет 36 независимых компонент, (1)lagd имеет 20 независимых компонент, (2)lagd имеет 15 независимых компонент, (3)lagd имеет 1 независимую компоненту. Далее будем рассматривать только часть (1)lagd . Запишем материальные уравнения: где i и i - тензоры диэлектрической и магнитной проницаемостей, (1)gi и (2)gi - перекрёстные члены. Учитывая структуру тензоров Fa и Ga (4), а также уравнения связи (16), запишем 5. Линейная локальная геометризация уравнений Максвелла Плебаньским была предложена простейшая геометризация уравнений Максвелла [4, 5, 23, 24]. Несмотря на некоторые недостатки, этот метод нашёл применение, например, в трансформационной оптике, где соотношение i = i является желательным [7, 25, 26]. Основная идея геометризации по Плебаньскому заключается в следующем: • Записать уравнения Максвелла в среде в пространстве Минковского. • Записать вакуумные уравнения Максвелла в эффективном римановом пространстве. • Приравнять соответствующие члены уравнений. • В результате мы получим выражение диэлектрической и магнитной проницаемостей через геометрические объекты. Однако данный подход к геометризации выглядит скорее как трюк. Автор постарался выполнить данные вычисления более формально. 5.1. Вспомогательные соотношения для метрического тензора Нам понадобятся простые соотношения для метрического тензора. Выражение (19) приводит к следующим частным соотношениям: Соотношение (20) перепишем в виде Подставляя (22) в (21), получаем: Это соотношение будет использовано позднее для упрощения записи итоговых уравнений. • Уравнения связи для движущихся сред Минковским были выведены уравнения связи для изотропных движущихся сред [16, 27] (уравнения Минковского для движущихся сред). Пусть ua - 4-скорость среды. Считая диэлектрическую и магнитную проницаемости и скалярами, можно записать В трёхмерном виде уравнения (24) принимают следующий вид: Тамм расширил уравнения (25) для анизотропного случая [1, 3], а именно, считая, что диэлектрическая и магнитная проницаемости имеют вид и вектор скорости ui системы отсчёта - параллельным одной из главных осей анизотропии. Тогда уравнения Минковского для движущихся сред приобретут следующий вид: • Общая геометризация Проведём геометризацию уравнений Максвелла. Введём эффективную метрику в ga . Тогда запишем лагранжиан электромагнитного поля в виде лагранжиана Янга-Миллса: Построим тензор lagd следующим образом: Тогда уравнение (13) примет следующий вид: Для наглядности распишем по компонентам: На основании (17) получим: На основании (17) с учётом соотношения (23) получим: Из (28) можно формально выписать выражение для диэлектрической проницаемости: При этом геометрический смысл второго члена в (28) нуждается в дальнейшем уточнении. Из (29) можно формально выписать выражение для магнитной проницаемости: Таким образом, геометризованные уравнения связи координатах имеют следующий вид: Леонгард предложил интерпретировать перекрёстный член в уравнениях (32) как скорость движения геометризованной системы отсчёта [6]. Действительно, на основании (25) уравнения (32) можно переписать в виде: где ui - трёхмерная скорость движения системы отсчёта. 6. Заключение Автор предложил формальный подход к проблеме геометризации уравнений электромагнитного поля. В качестве иллюстрации представлен метод локальной линейной геометризация уравнений Максвелла. Следует отметить, что данный метод нельзя считать полностью удовлетворительным. Действительно, полный тензор проницаемостей lagd имеет 36 независимых компонент. И даже его основная часть, тензор (1)lagd , имеет 20 независимых компонент. В то время как риманов метрический тензор имеет только 10 независимых компонент. Также следует сказать, что предложенный метод отличается по идеологии и результатам от метода Плебаньского.×
Об авторах
Дмитрий Сергеевич Кулябов
Российский университет дружбы народов
Email: ds@sci.pfu.edu.ru
Кафедра прикладной информатики и теории вероятностей; Лаборатория информационных технологий Объединённый институт ядерных исследований ул. Жолио-Кюри, д. 6, Дубна, Московская область, Россия, 141980 ул. Миклухо-Маклая, д. 6, Москва, Россия, 117198
Список литературы
- Тамм И. Е. Электродинамика анизотропной среды в специальной теории относительности // Журнал Русского физико-химического общества. Часть физическая. - 1924. - Т. 56, № 2-3. - С. 248-262.
- Тамм И. Е. Кристаллооптика теории относительности в связи с геометрией биквадратичной формы // Журнал Русского физико-химического общества. Часть физическая. - 1925. - Т. 57, № 3-4. - С. 209-240.
- Tamm I. E., Mandelstam L. I. Elektrodynamik der anisotropen Medien in der speziellen Relativitatstheorie // Mathematische Annalen. - 1925. - Bd. 95, No. 1. - Ss. 154-160.
- Plebanski J. Electromagnetic Waves in Gravitational Fields // Physical Review. - 1960. - Vol. 118, No 5. - Pp. 1396-1408.
- Felice F. On the Gravitational Field Acting as an Optical Medium // General Relativity and Gravitation. - 1971. - Vol. 2, No 4. - Pp. 347-357.
- Leonhardt U., Philbin T. G., Haugh N. General Relativity in Electrical Engineering. - 2008. - Pp. 1-19.
- Leonhardt U., Philbin T. G. Transformation optics and the geometry of light // Progress in Optics. - 2009. - Vol. 53. - Pp. 69-152.
- Thompson R. T., Cummer S. A., Frauendiener J. A Completely Covariant Approach to Transformation Optics // Journal of Optics. - 2011. - Vol. 13, No 2. - P. 024008.
- Пенроуз Р., Риндлер В. Спиноры и пространство-время. Два-спинорное исчисление и релятивистские поля. - М.: Мир, 1987. - Т. 1.
- Сивухин Д. В. О Международной системе физических величин // Успехи физических наук. - 1979. - Т. 129, № 10. - С. 335-338.
- Kulyabov D. S., Korolkova A. V., Korolkov V. I. Maxwell’s Equations in Arbitrary Coordinate System // Bulletin of Peoples’ Friendship University of Russia. Series “Mathematics. Information Sciences. Physics”. - 2012. - No 1. - Pp. 96-106.
- Korol’kova A. V., Kulyabov D. S., Sevast’yanov L. A. Tensor Computations in Computer Algebra Systems // Programming and Computer Software. - 2013. - Vol. 39, No 3. - Pp. 135-142.
- Kulyabov D. S. Geometrization of Electromagnetic Waves // Mathematical Modeling and Computational Physics. - Dubna: JINR, 2013. - P. 120.
- Кулябов Д. С., Королькова А. В. Уравнения Максвелла в произвольной системе координат // Вестник Тверского государственного университета. Серия: Прикладная математика. - 2013. - № 1 (28). - С. 29-44.
- Кулябов Д. С., Немчанинова Н. А. Уравнения Максвелла в криволинейных координатах // Вестник РУДН. Серия «Математика. Информатика. Физика». - 2011. - № 2. - С. 172-179.
- Minkowski H. Die Grundlagen fu¨r die electromagnetischen Vorg¨onge in bewegten K¨orpern // Nachrichten von der Gesellschaft der Wissenschaften zu G¨ottingen, Mathematisch-Physikalische Klasse. - 1908. - Ss. 53-111.
- Стрэттон Д. А. Теория электромагнетизма. - М.-Л.: ГИТТЛ, 1948.
- Терлецкий Я. П., Рыбаков Ю. П. Электродинамика. - Москва: Высшая школа, 1990. - С. 352.
- Post E. The Constitutive Map and Some of its Ramifications // Annals of Physics. - 1972. - Vol. 71, No 2. - Pp. 497-518.
- Gilkey P. B. Algebraic Curvature Tensors // Geometric Properties of Natural Operators Defined by the Riemann Curvature Tensor. - World Scientific Publishing Company, 2001. - Pp. 1-91.
- Obukhov Y. N., Hehl F. W. Possible Skewon Effects on Light Propagation // Physical Review D - Particles, Fields, Gravitation and Cosmology. - 2004. - Vol. 70, No 12. - Pp. 1-14.
- Hehl F. W., Obukhov Y. N. Linear Media in Classical Electrodynamics and the Post Constraint // Physics Letters, Section A: General, Atomic and Solid State Physics. - 2005. - Vol. 334, No 4. - Pp. 249-259.
- Leonhardt U. Optical Conformal Mapping // Science. - 2006. - Vol. 312. - Pp. 1777-1780.
- Pendry J. B., Schurig D., Smith D. R. Controlling Electromagnetic Fields // Science. - 2006. - Vol. 312, No 5781. - Pp. 1780-1782.
- Nicolet A., Zolla F., Geuzaine C. Transformation optics, generalized cloaking and superlenses // IEEE Transactions on Magnetics. - 2010. - Vol. 46, No 8. - Pp. 2975- 2981.
- Schurig D., Pendry J. B., Smith D. R. Calculation of Material Properties and Ray Tracing in Transformation Media // Optics express. - 2006. - Vol. 14, No 21. - Pp. 9794-9804.
- Зоммерфельд А. Электродинамика. - Москва: Издательство иностранной литературы, 1958.