УДК 535.4
DOI: 10.14529/mmph230209
ПОЛЯРИЗАЦИОННЫЕ ЭФФЕКТЫ ПРИ ОТРАЖЕНИИ СВЕТА ОТ ГРАНИЦЫ ОДНООСНОГО КРИСТАЛЛА
Ю.В. Мухин1'2, Ю.Р. Гатауллина1, Н.Д. Кундикова1'2
Южно-Уральский государственный университет, г. Челябинск, Российская Федерация 2Институт электрофизики УрО РАН, г. Екатеринбург, Российская Федерация e-mail: mukhinyv@susu.ru
Аннотация. Получены и исследованы граничные условия для плоских электромагнитных волн на границе раздела изотропной среды с одноосным кристаллом. Получены элементы матрицы отражения, связывающей компоненты электрического поля в отраженной волне с компонентами электрического поля в падающей волне. Показана возможность сохранения поляризации при отражении от границы с анизотропной средой. Получены условия сохранения векторов Джонса при отражении от поверхности одноосного кристалла. Исследована возможность существования эффекта Брю-стера на границе с такой анизотропной средой. Брюстеровские поляризации оказываются близкими к ^-поляризациям с небольшой примесью s-поляризации.
Ключевые слова: граница раздела с анизотропной средой; вектор Джонса; угол Брюстера.
Введение
Отражение электромагнитных волн от границы с анизотропной средой рассматривалось неоднократно (см. например [1]). Тем не менее детальные исследования поляризационных эффектов при отражении света от таких границ не проводились. Данная работа восполняет частично этот пробел. При отражении от границы двух изотропных сред поляризация волны никогда не сохраняется (за исключением двух тривиальных случаев s- или р-поляризованной падающей волны). Плоскость поляризации после отражения всегда прижимается ближе к границе раздела сред [2]. В случае анизотропных сред в отраженной волне всегда присутствуют обе s- и р-поляризации, что допускает возможность сохранения векторов Джонса при отражении. Чтобы определить, при каких условиях это выполняется, были получены и исследованы граничные условия для плоских электромагнитных волн в случае одноосного кристалла с осью, параллельной границе раздела кристалла и изотропной среды.
Граничные условия
Пусть электромагнитная волна падает под углом падения 9 на границу раздела изотропной среды и одноосного кристалла, как показано на рис. 1. На рис. 2 показаны направления волновых векторов в плоскости падения OXZ для волн в изотропной среде (обозначены как E) и для обыкновенной и необыкновенной волн в анизотропной среде (обозначены как D). Для волн, распространяющихся «сверху вниз», используется индекс «+», а для бегущих «снизу вверх» используется индекс «-».
Запишем проекции векторов напряженности электрического поля Е и напряженности магнитного Н на границу раздела в изотропной среде. Направления векторов для разных поляризаций выбираются так, как показано на рис. 3 и 4, таким образом, чтобы при нормальном падении направления напряженностей электрических полей в падающей и отраженной волнах совпадали.
Суммарные проекции полей записываются отдельно для каждой составляющей: для падающей р-поляризованной волны (см. рис. 3) и s-поляризованной волны (см. рис. 4). Epx = cosd (Ep + E~p) Hpx = 0 Esx = О Rsx = ncosd (Ep -E-)
Epy = 0 Hpy =-n (Ep - Ep) Esy = (К- Щ) Hsy = 0
Ez = sin д (Ep - Ep) Hpz = 0 Esz = 0 Hz = n sin0 (Ep+E-)
Рис. 1. Граница раздела OXY содержит оптическую ось Рис. 2. Расположение направлений волновых векторов
кристалла, расположенную под углом ß к оси OX
в плоскости падения OXZ
Рис. 3. Падающая р-поляризация
Рис. 4. Падающая s-поляризация
Для записи проекций полей в одноосном кристалле удобно пользоваться индукцией электрического поля вместо напряженности. Это обусловлено поперечностью поля D по отношению к волновому вектору. Введём также обозначения углов, связанных с ориентацией полей и волновых векторов, упрощающие вид формул:
sin Qe = n sin Q / ne cosye = cos ß sin Qe cos ae = sin ß sin Qe
sin Q0 = n sinQ/Je^ cosya = cos ß sin Q0 cosa0 = sin ß sin Q0
Перпендикулярные и параллельные оптической оси составляющие электрического E и магнитного H полей необыкновенной и обыкновенной волн выражаются тогда через индукцию электрического поля D следующим образом.
Для необыкновенной волны:
Ее ,=^d;
— i
Н4\ =°
H
D+
в±
НеП =0
H,
_ D-
Для обыкновенной волны:
Ео\\=°
D+
К„=о
F =D-
Fo! =-
D+
Н<А\=*тГ, °
4е!
D+
Ho ! = cosrt 0
4е!
Но\\ =-sin rc
D-
yßl
Ho != cosTo
D-
4е!
Индексы «е» и «о» обозначают компоненты необыкновенной и обыкновенной волн, которые
4ё[ . Здесь п - по-
имеют показатели преломления ne = lei
1 -n¿
S\ i S I
v £II£J- У
sin2 в cos2 ß
и n
казатель преломления изотропной среды.
Граничные условия заключаются в непрерывности тангенциальных составляющих электрического и магнитного полей на границе раздела. Остается только «сшить» все полученные выражения в четыре граничных условия:
n
n
e
e
е
!
!
Для электрического поля, параллельного оптической оси:
cos Pcos0(E+ + E- ) + sin P(E+ + E- ) = ^^ (D+ + D- )
p p
Для магнитного поля, параллельного оптической оси:
-n sin P(E+ - E- ) + n cos P cos 0(E+ - EJ ) =
sin
Фа.
Для электрического поля, перпендикулярного оптической оси: - sin P cos в(E+ + E- ) + cos P(E+ + E- ) =
(d;+ d- )
cos Ге sinae(D + +£>-) + cos aos
(D++ D- )
Для магнитного поля, перпендикулярного оптической оси:
n cos P(Ep+ - E- ) + n sin P cos 0(E+ - E- ) =
cosa
(Dp D- )--
4^2
(1) (2) (3)
e,n+ п-л cos^osinao (d+ _D-) (4)
Матрица отражения
Полученная выше алгебраическая система уравнений позволяет вычислить различные амплитудные коэффициенты отражения для разных поляризаций. В уравнениях (1)-(4) можно избавиться от компонент индукции электрического поля Б и получить элементы матрицы отражения М, связывающей компоненты электрического поля в отраженной волне с компонентами электрического поля в падающей:
С E Л f
r r
pp sp
E„ ^rps rss J{ E,
\ s
=M
\ s J
fE Ep
E+
\ Es J
Здесь элементы матрицы отражения М даются выражениями:
при Ер = 0 ;
pp
ps
E- A1 pC2 ' A2 pC1
E+ B2 pC1 - B1 pC2
E- A1 pB2 p - A2pB1 p
E+ Ep B1pC2 - B2pC1
rsp
, пРи E+= 0; rss =-±=
El= B2sCl - BuC2 , при E+= 0 ; J7+ AC - AC p
E- B1sA2s - B2sA1s
EP A1sC2 A2sC1
при E+ = 0.
Где введены следующие обозначения:
A1p,s = a3p,sd1e2 - a1p,sd3e2 - a2p,sd1e3 A2p,s = a4p,sd1e2 - a1p,sd4e2 - a2p,sd1e4
CC1 — сз^2 с1$зе2 ^2^1ез
B1 p,s = b3 p,sd1e2 b1 p,sd3e2 b2 p,sd1e3 B1 p,s = b3p,sd1e2 - b1 p,sd3e2 - b2p,sd1e3 C2 — С4^е с^4е2 С2^1е4
a1 = cosacosP;a2 = -nsinP;a2s = nsinP;a3 = -a1;a4p = -ncosP; a4s = -a4 ;
b1 p = a; b1s =—;b2p = a2s; К = -na; b3 p =-a; К = a; b4 p = a4s; n n
и o o J sin sin ro
b4s = -na1; c1 = sin P; c2 = -na1; c3 = cos P; c4 = na1; d1 =-e; e2 =-o ■
4^2'
sin a cos^ cosa , - cosa sin a cosr„
И —_e_ ." —_— И —_— ^ —- o
d3 = ; e3 = ; d4 = ; e4 =
пе
В отсутствие поглощения и полного внутреннего отражения матрица М должна обладать определёнными свойствами. Во-первых, она должна быть вещественной матрицей. Помимо этого, по аналогии с линейными эрмитовыми операторами квантовой механики, она обязана быть симметричной. Симметричность М подтверждалась прямыми выкладками и вычислениями.
s
n
e
s
Эффект сохранения поляризации
В изотропных средах поляризация света никогда не сохраняется при отражении от их границ раздела. На границе с анизотропной средой ситуация совсем иная. Представим отраженную в изотропной среде волну в виде
V Е—У
(у
рр
г
V ря
яр
Е"
V
= г
Е+
Е'
V Ея У
(5)
Алгебраическая система уравнений (5) может иметь нетривиальное решение только в случае равенства нулю детерминанта:
<1^
( г — г
рр
г
V ря
\
яр
г — г
= 0
Решение системы уравнений (5) и нахождение коэффициентов у аналогично решению задачи о расщеплении двукратно вырожденного уровня, разобранной в [3].
:0,5(
г + г +
рр яя —
3 =
( Т7+ \
Е'
VЕя У
) , ™ = ^(грр Р гяя )2 + 4| гря I2 вол
Выражения для векторов Джонса У падающей волны имеют при этом следующий вид:
( I-7-Г А
5 1 + -
I г^р IV ' I
Наглядно представить поляризацию электрического поля удобно с помощью годографов нормированных векторов Джонса: точек или линий на поверхности сферы единичного радиуса в используемой системе координат. Волновые векторы при этом будут лежать в плоскости падения перпендикулярно векторам Джонса. В зависимости от угла падения в годограф вектора Джонса для одной из поляризаций прочерчивает на поверхности сферы линию от экватора (нормальное падение) и до полюса, где угол падения света в = 90° . Для различных положений оптической оси, разных углов в, получаются разные линии на сфере. На рис. 5 и 6 приведены годографы векторов Джонса У для сохраняющихся поляризаций: рис. 5 - для углов [ < 45°, рис. 6 - для углов 45° <[< 90° .
Рис. 5. Годографы вектора Джонса для углов в от 0° Рис. 6. Годографы вектора Джонса для углов в от 45° до 45° с шагом 9°. ^ = 1,7;^е2 = 1,55 до 90° с шагом 9°. ^ = 1,7;^/^ = 1,55
На рис. 6 видно, что для 45°<[<90° плоскости сохраняющихся поляризаций могут принимать лишь ограниченные положения: ни одна из линий годографов на рис. 6 не пересекает вертикальные окружности на сфере, параллельные плоскости падения, в то время как для углов [ < 45° каждая из линий пересекает эти окружности. Вертикальные окружности на рис. 5 и 6 соответствуют плоскостям поляризаций, составляющим угол в 45° с плоскостью падения. Данное
свойство не зависит от оптических параметров кристалла, но может быть пояснено с помощью рис. 7, на котором показаны графики зависимости от угла падения в двух коэффициентов матрицы отражения при ¡5 = 10°.
Кривые на графике имеют точку пересечения, которая присутствует для всех углов [ < 45°. Это пересечение соответствует равенству модулей компонент векторов Джонса для сохраняющихся поляризаций, и значит, расположению плоскости поляризации под углом 45° к плоскости падения. Для углов 45° <[<90° такое пересечение кривых отсутствует, именно поэтому положение плоскостей сохраняющихся поляризаций в этом случае ограничено областями, показанными на рис. 6.
Угол Брюстера
Для границы двух изотропных сред угол Брюстера соответствует нулевому отражению p-поляризованного света. На границе с анизотропной средой p-поляризованный свет всегда будет давать s-поляризованную составляющую в отраженной волне. Поэтому сам термин и условия «Брюстера» следует относить здесь уже к другой поляризации, которая, очевидно, будет определяться углом в оптической оси к плоскости падения. Если условия Брюстера выполнены для одной поляризации, то ортогональная ей поляризация должна сохраняться.
Таким образом, условия Брюстера являются частным случаем сохраняющихся поляризаций с нулевым коэффициентом r в системе уравнений (5). Чтобы это было возможным, необходимо равенство нулю детерминанта матрицы отражения в (5): rpprss - rpsrsp = 0.
На рис. 8 показаны зависимости rpprss и rpsrsp от угла падения в и угла в , взятых в полярных для этих переменных координатах: в sin [ и в cos [ . Пересечение двух поверхностей как раз и соответствует выполнению условий Брюстера. При выполнении условия вектор Джонса падающей волны приобретает вид
Г 1 ^
J~ .
— г г
V ps/ SS J
Условия Брюстера являются частным случаем сохранения поляризации при отражении. Другими словами, на линиях, изображенных на рис. 5 и 6, должны быть точки, соответствующие Брюстеровским поляризациям.
На рис. 9 показана макушка сферы единичного радиуса со следами годографов сохраняющихся векторов Джонса, на которых отмечены точки, соответствующие условиям Брюстера. Оказывается, что эти условия реализуются между углами Брюстера для
коэффициентов преломления ne (Р = 0) =
1 -п1
V £II£J- У
sin2 в
и no=*Js^. Брюстеровская
поляризация оказывается при этом близкой к ^-поляризации с небольшой примесью 5-поляризации на уровне нескольких процентов, но которая зависит уже от параметров кристалла.
Рис. 9. Годограф векторов Джонса для углов Брюстера. ^^ = = 1,55
Заключение
В настоящей работе были получены и исследованы граничные условия для плоских электромагнитных волн на границе раздела изотропной среды с одноосным кристаллом. Получены элементы матрицы отражения связывающей компоненты электрического поля в отраженной волне с компонентами электрического поля в падающей волне. Показана возможность сохранения поляризации при отражении от границы с анизотропной средой и получены условия для реализации этого эффекта. Исследована возможность существования эффекта Брюстера на границе с такой анизотропной средой. Брюстеровские поляризации оказываются близкими к ^-поляризациям с небольшой примесью s-поляризаций.
Литература
1. Игнатович, Ф.В. Оптика анизотропных сред / Ф.В. Игнатович, В.К. Игнатович // Успехи физических наук. - 2012. - Т. 182, № 7. - С. 759-771.
2. Ландау, Л.Д. Теоретическая физика: учебное пособие для студентов физических специальностей университетов. В 10 т. Том VIII. Электродинамика сплошных сред / Л.Д. Ландау, Е.М. Лифшиц. - Москва: Физматлит: МАИК Наука, 2005. - 651 с.
3. Ландау, Л.Д. Теоретическая физика: учебное пособие для студентов физических специальностей университетов. В 10 т. Том 3: Квантовая механика. Нерелятивистская теория / Л.Д. Ландау, Е.М. Лифшиц. - Москва: Физматлит: МАИК Наука, 2001. - 803 с.
Поступила в редакцию 25 апреля 2023 г.
Сведения об авторах
Мухин Юрий Викторович - кандидат физико-математических наук, Ph.D., доцент, кафедра оптоинформатики, Южно-Уральский государственный университет, г. Челябинск, Российская Федерация; Институт электрофизики УрО РАН, г. Екатеринбург, Российская Федерация, e-mail: mukhinyv@susu.ru
Гатауллина Юлия Рауфовна - студент, Южно-Уральский государственный университет, г. Челябинск, Российская Федерация, e-mail: gataullinaju@gmail.com
Кундикова Наталия Дмитриевна - доктор физико-математических наук, профессор, кафедра оптоинформатики, Южно-Уральский государственный университет, г. Челябинск, Российская Федерация; Институт электрофизики УрО РАН, г. Екатеринбург, Российская Федерация, е-mail: kundikovand@susu. ru
Bulletin of the South Ural State University Series "Mathematics. Mechanics. Physics" _2023, vol. 15, no. 2, pp. 66-72
DOI: 10.14529/mmph230209
THE POLARISATION EFFECTS OF LIGHT REFLECTION FROM THE INTERFACE WITH A UNIAXIAL CRYSTAL
Y.V. Mukhin12, Y.R. Gataullina1, N.D. Kundikova1,2
South Ural State University, Chelyabinsk, Russian Federation 2 Institute of Electrophysics, Ural Division of the RAS, Ekaterinburg, Russian Federation e-mail: mukhinyv@susu.ru
Abstract. In this paper, the boundary conditions for plane electromagnetic waves at the interface between an isotropic medium and a uniaxial crystal have been obtained and analyzed. Elements of the reflection matrix that relates components of the electric field of the reflected wave back to those of the incident electromagnetic wave have been determined. Using the elements of the reflection matrix, we show that the possibility exists for polarization of light to be preserved upon reflection at the interface with anisotropic media. The conditions for the Jones vectors conservation upon reflection from the surface of a uniaxial crystal are obtained. The existence of the Brewster effect at the interface with such an anisotropic medium has been studied and verified. The Brewster's polarizations have been shown to be almost ^-polarizations with just little admixture of s-polarization.
Keywords: interface with anisotropic media, Jones vector, Brewster angle.
References
1. Ignatovich F.V, Ignatovich V.K. Optics of anisotropic media. Physics-Uspekhi, 2012, Vol. 55, Iss. 7, pp. 709-720. DOI: 10.3367/UFNe.0182.201207f.0759
2. Landau L.D., Lifshits E.M. Teoreticheskaya fizika: uchebnoe posobie dlya studentov fizicheskikh spetsial'nostey universitetov. V10 t. Tom VIII. Elektrodinamika sploshnykh sred (Theoretical physics: a textbook for students of physical specialties of universities. In 10 vols. Vol. VIII. Electrodynamics of continuous media), 2005, Moscow, Fizmatlit Publ., MAIK Nauka Publ., 651 p. (in Russ.).
3. Landau L.D., Lifshits E.M. Teoreticheskaya fizika: uchebnoe posobie dlya studentov fizicheskikh spetsial'nostey universitetov. V10 t. Tom 3: Kvantovaya mekhanika. Nerelyativistskaya teoriya (Theoretical physics: a textbook for students of physical specialties of universities. In 10 vols. Vol. III: Quantum Mechanics. Non-relativistic theory, 2001, Moscow, Fizmatlit Publ., MAIK Nauka Publ., 803 p. (in Russ).
Received April 25, 2023
Information about the authors
Mukhin Yurii Viktorovich is Cand. Sc. (Physics and Mathematics), Ph.D., Associate Professor, Optoinformatics Department, South Ural State University, Chelyabinsk, Russian Federation; Institute of Electrophysics of the Ural Branch of the Russian Academy of Sciences, Yekaterinburg, Russian Federation, e-mail: mukhinyv@susu.ru
Gataullina Julia Raufovna is student, Optoinformatics Department, South Ural State University, Chelyabinsk, Russian Federation, e-mail: gataullinaju@gmail.com
Kundikova Nataliya Dmitrievna is Dr. Sc. (Physics and Mathematics), Professor, Optoinformatics Department, South Ural State University, Chelyabinsk, Russian Federation; Institute of Electrophysics of the Ural Branch of the Russian Academy of Sciences, Yekaterinburg, Russian Federation, e-mail: kundikovand@susu.ru