CC BY
86
УДК 537.876.22+535.56 В.В. Фисанов
Электромагнитные волны в изотропной метакиральной среде
Выполнена декомпозиция электромагнитного поля в прозрачной метакиральной биизотроп-ной среде на два поля Бельтрами, которые характеризуются неравными положительными волновыми числами и общим волновым импедансом. Знак волнового импеданса определяется знаками диэлектрической и магнитной проницаемостей среды. Получено условие различения прямых и обратных нормальных волн круговой поляризации.
Ключевые слова: метакиральная среда, диэлектрическая проницаемость, магнитная проницаемость, параметр киральности, электромагнитные поля Бельтрами, волновой импеданс, положительные волновые числа, прямые и обратные плоские волны.
Искусственным композитным электромагнитным материалам с экстремальными (в том числе отрицательными) значениями диэлектрической и магнитной проницаемостей в настоящее время уделяется повышенное внимание вследствие многообразия связанных с ними необычных волновых явлений, таких как отрицательное преломление, плоская линза, маскировка объектов. При наличии дополнительного параметра киральности такие материалы называют метакиральными средами. В отличие от обычных киральных сред в них допускается возможность существования не только прямых, но и обратных нормальных волн круговой поляризации. Обратная волна возникает даже при неотрицательных значениях проницаемостей [1, 2]. Следуя В.Г. Веселаго [3], в электродинамике метаматериалов при описании обратных волн применяют термины «отрицательный коэффициент (показатель) преломления», «отрицательное волновое число». Между тем, согласно каноническому определению, волновое число является длиной (модулем) волнового вектора и по этой причине не может принимать отрицательные значения. Налицо присутствует методический парадокс, который требует разрешения.
С этой целью в данной работе для исследования прямых и обратных нормальных волн в биизо-тропной метакиральной среде применяется новый подход, при котором волновые числа в отсутствие диссипативных потерь полагаются безусловно положительными величинами и применяется уточнённая трактовка волнового импеданса.
Постановка задачи, декомпозиция уравнений поля. Для феноменологической характеристики киральных сред используют различные системы материальных уравнений, но применительно к метакиральной среде наиболее употребительным является ЕН-представление (представление Тел-легена), согласно которому векторы индукций Б и В выражаются через векторы напряжённостей Е и Н электрического и магнитного поля в среде (например, [1, 4, 5]). Это представление происходит от предложенных Кондоном [6, 7] и усовершенствованных Силверманом [8] уравнений, предназначенных для описания явления оптической активности. В случае монохроматического поля круговой частоты га и временного фактора ехр(-/га?) система материальных уравнений имеет вид
Б = еЕ+/кН, В = цН -/кЕ , (1)
где параметр киральности к = ю£ пропорционален коэффициенту гирации g Кондона-Силвермана,
/=7-1 . Диэлектрическая проницаемость е, магнитная проницаемость ц и параметр к принимаются здесь и далее действительными величинами. Энантиоморфные разновидности киральной среды различаются по знаку параметра киральности к . Уравнения Максвелла после подстановки в них уравнений (1) принимают вид
УхЕ = /га(цН-/кЕ), УхН = -/га(еЕ + /кН), V-Е = 0, V-Н = 0. (2)
Далее обычно переходят к дифференциальным уравнениям второго порядка относительно вектора Е или вектора Н [4]. Борен [9] предложил иной подход: перейти к уравнениям первого порядка для волновых полей левой круговой поляризации (комплексная амплитуда 01) и правой круговой поляризации (комплексная амплитуда 02 ) посредством линейного преобразования
Е = 01 -гл02, Н = 02 -/Ш, (3)
где п = С 1 - волновой импеданс. Поля и О2 пропорциональны своим роторам, т.е. подчиняются дифференциальным уравнениям
Ух О! =у101, V-О! = 0, Ух02 = -^2, V-02 = 0 (4)
и поэтому называются электромагнитными полями Бельтрами. В обычной киральной среде положительные величины У1 и у 2 являются волновыми числами соответствующих полей Бельтрами. Действительно, в силу своей соленоидальности поля Бельтрами удовлетворяют уравнению Гельмгольца
2-2 ~ вида V О у + у у О у = 0 , где у = 1,2 . В метакиральной среде величины у у могут изменить знак, поэтому их следует считать просто собственными значениями дифференциального оператора вихря. В терминах гидродинамики величины у у являются постоянными коэффициентами пропорциональности между завихренностью (Ух Оу ) и скоростью (0у) векторных потоков Бельтрами-Тркала. Волновые числа у у определим через них по формулам
У1 = У1 вБпУ1, у2 = У2 вБпу2 . (5)
Для отыскания волновых чисел подставим формулы (3) в вихревые уравнения Максвелла (2), заменим роторы согласно (4) и перегруппируем слагаемые так, чтобы однотипные поля располагались по разные стороны от знака равенства:
(у 1 -ю к-ю ц^) =-(пУ2 -ю ц+ю кп) О2 , (6)
-г(СУ1 - юе-юкС)О1 = (у 2 + юк-юеп)О 2. (7)
Приравнивая к нулю коэффициенты при полях и О2 , получим следующие соотношения:
У1 =ю(к + ц^) = ю(к + еп), у2 =-ю(к-еп) = -ю(к-ц^). (8)
2
Из формул (8) следует, что еп = |С и п = |/е . Далее, при к = 0 (некиральная среда) величины у1 и у 2 становятся равными волновому числу для полей круговой поляризации в обычной изотропной среде, которое обозначим символом к, так что
к = юеп = юцС . (9)
Так как волновое число к является положительной величиной, из формулы (9) непосредственно следует, что знак волнового импеданса п (и волнового адмитанса ^) определяется знаком прони-цаемостей. В обычном диэлектрике (положительные е и ц) волновой импеданс положительный,
равен арифметическому квадратному корню из отношения проницаемостей: п = л/ц/е . Однако в среде Веселаго (т.е. в среде с проницаемостями е<0 и ц<0) импеданс является отрицательной величиной. Обобщённо он вычисляется по формуле
п=Л/Це8§пе = ч/|Ц8 8§пц, (10)
при этом имеется в виду, что в прозрачной среде знаки диэлектрической и магнитной проницаемо-стей совпадают. Формула (10) является справедливой не только для изотропной среды, но и для би-изотропной метакиральной среды с материальными уравнениями вида (1). С учётом (9) формулы (8) принимают вид
У1 = к + юк, у2 = к-юк . (11)
Как следует из формулы (11), случай метакиральности реализуется при выполнении неравенства ец<к2 : при любом знаке параметра киральности к одна из величин у у оказывается отрицательной. Вводя показатель преломления п = к/к0 , где к0 - волновое число для вакуума, указанное неравенство можно преобразовать к виду п <|^/л/|0е0 , что соответствует либо сильной киральности, либо малому значению показателя преломления эквивалентной некиральной среды (80, Ц0 - электрическая и магнитная постоянные).
Электромагнитные поля Бельтрами О у распространяются без взаимодействия в однородной и
безграничной среде. Переносимая ими энергия характеризуется векторами Пойнтинга
S j - Sj s '' -
§ j = (1/2)Яе(Е j х Н* ) , где Sj - длина вектора Пойнтинга ] -го поля (положительная величина), §j - его орт. Определим плоские волны Бельтрами как Рj(г) = QQjехр(/'удj-г), где уj - орт волнового вектора уj в направлении движения фазового фронта волны, Роj - комплексная амплитуда. Вместо системы (4) получается совокупность векторных алгебраических уравнений
71 х р01 =-?'р018Еп7Ь '71 ^01 = 0; 72 х р02 = ^02 8Еп72, "72 ' р02 = 0. Для действительных значений импеданса п вектор Пойнтинга левой волны Бельтрами 81 = SlSl определяется выражением
81 = (]/2)Ке(/^01 х Р01) = (1/2) С р012 71 ®§п 71. (12)
Соответственно для вектора Пойнтинга правой волны Бельтрами 82 = S2S2 имеем
82 = (1/2)1^-/^02 хР02) = (1/2)пр02|272®§п72 . (13)
В метакиральной изотропной среде направления переноса энергии (луча) и движения фронта волны могут либо совпадать (прямая волна), либо быть противоположными (обратная волна), т.е.
у j = (§ j у j)§ j. Подставляя это соотношение в (12) и (13) и учитывая, что Sj > 0, получим универсальное неравенство п(§ j 7 j)пуj > 0 . Из него сразу следует ключевая формула
у j-7 j = «Еп (п7 j ) (14)
для определения типа плоских волн Бельтрами: прямой волне соответствует значение знаковой функции +1, обратной волне - значение -1. Следует учитывать, что согласно (8) оба числа уj с
номерами j = 1 и j = 2 не могут одновременно принимать отрицательные значения. Формула (14) даёт возможность уточнить классификацию нормальных волн круговой поляризации в изотропной киральной среде. В обычной киральной среде с положительными проницаемостями, называемой также средой Пастера [10], положительными являются волновой импеданс п и оба числа у j =уj .
Метакиральную среду с двумя отрицательными проницаемостями было предложено назвать средой Пастера-Веселаго [11]. Она характеризуется отрицательным волновым импедансом п и допускает
режим с двумя обратными волнами, если оба числа уj являются положительными. Смешанные режимы с прямой и обратной нормальными волнами реализуются как при положительном, так и отрицательном волновом импедансе при условии, что числа уj имеют противоположные знаки. Если
в среде изменяется только знак параметра киральности, то смешанный режим сохраняется, но каждая нормальная волна изменяет свой первоначальный тип противоположным образом.
Отказ от концепции отрицательного волнового числа позволяет устранить отмеченное ранее [11, 12] различие в знаках диэлектрической и магнитной проницаемостей при переходе от системы материальных уравнений в представлении Теллегена к уравнениям в форме Друде-Борна-Фёдорова, которое происходит при достижении сильной киральности.
Заключение. Описание обратной волны в метакиральной среде не требует введения для неё нового понятия «отрицательное волновое число», но нуждается в уточнении трактовки волнового импеданса в случае отрицательных значений и диэлектрической, и магнитной проницаемостей. Также следует делать различие между волновым числом и соответствующим ему параметром (собственным числом) уравнения для поля Бельтрами. Плоская волна Бельтрами является обратной при условии, что знаки у волнового импеданса и у собственного числа оказываются противоположными.
Литература
1. Waves and energy in chiral nihility / S. Tretyakov, I. Nefedov, A. Sihvola, S. Maslovski, C. Simovski // J. of Electromagn. Waves and Appl. - 2003. -Vol. 17, № 5. - P. 695-706.
2. Pendry J. A chiral route to negative refraction // Science. - 2004. - Vol. 306, № 5700. - P. 13531355.
3. Веселаго В.Г. Электродинамика веществ с одновременно отрицательными значениями s и ц // УФН. - 1967. - Т. 92, вып. 3. - С. 517-526.
4. Chern R.-L. Wave propagation in chiral media: composite Fresnel equations // J. Opt. - 2013. -Vol. 15, № 7. - P. 075702-1-075702-7.
5. Cao Y. Complete band gaps in one-dimensional photonic crystals with negative refraction arising from strong chirality / Y. Cao, J. Li // Phys. Rev. B. - 2014. - Vol. 89, № 11. - P. 115420-1-115420-5.
6. Condon E.U. Theories of optical rotatory power // Rev. Mod. Phys. - 1937. - Vol. 9, № 4. - P. 432-457.
7. Кондон Е. Теория оптической вращающей способности // УФН. - 1938. - Т. 19, вып. 3. -С.380-430.
8. Silverman M.P. Reflection and refraction at the surface of a chiral medium: comparison of gyrotropic constitutive relations invariant or noninvariant under a duality transformation // J. Opt. Soc. Am. A. - 1986. - Vol. 3, № 6. - P. 830-837.
9. Bohren C.F. Light scattering by an optically active sphere // Chem. Phys. Lett. - 1974. - Vol. 21, № 3. - P. 458-462.
10. Sihvola A.H. Bi-isotropic constitutive relations / A.H. Sihvola, I.V. Lindell // Microwave Opt. Technol. Lett. - 1991. - Vol. 4, № 8. - P. 295-297.
11. Фисанов В.В. Инварианты изотропной киральной среды // Радиотехника и электроника. -2007. - Т. 52, № 9. - С. 1089-1091.
12. Zhang C. Spatial dispersion and energy in a strong chiral medium / C. Zhang, T.J. Cui // Optics Express. - 2007. - Vol. 15, № 8. - P. 5114-5119.
Фисанов Василий Васильевич
Д-р физ.-мат. наук, вед. науч. сотрудник СФТИ при НИТГУ, профессор каф. радиофизики
Национального исследовательского Томского государственного университета
Тел.: 8 (382-2) 41-20-78
Эл. почта: [email protected]
Fisanov V.V.
Electromagnetic waves in an isotropic metachiral medium
Electromagnetic-field decomposition in a transparent metachiral bi-isotropic medium is performed onto two Beltrami fields, which are characterized by unequal positive wave numbers and common wave impedance. Sign of the wave impedance is determined by signs of permittivity and permeability of the medium. The condition of distinction between forward and backward circularly polarized eigenwaves is obtained.
Keywords: metachiral medium, dielectric permittivity, magnetic permeability, chirality parameter, electromagnetic Beltrami fields, wave impedance, positive wave numbers, forward and backward plane waves.
