Поверхность. Рентгеновские, синхротронные и нейтронные исследования, 2020, № 7, стр. 79-84
Управляющие свойства границ раздела в нелинейных структурах типа сэндвича с дефокусирующим внутренним слоем
С. Е. Савотченко *
Белгородский государственный технологический университет
им. В.Г. Шухова
308012 Белгород, Россия
* E-mail: savotchenkose@mail.ru
Поступила в редакцию 27.12.2019
После доработки 26.01.2020
Принята к публикации 28.01.2020
Аннотация
Рассмотрена модель трехслойной оптической структуры, в которой плоскопараллельные границы обладают собственными нелинейными свойствами. Внутренний слой конечной толщины представляет собой оптически прозрачную среду с дефокусирующей керровской нелинейностью, снаружи контактирующий с диэлектрическими линейными полупространствами. Математическая формулировка модели сводится к нелинейному уравнению Шредингера с положительным коэффициентом кубической нелинейности и нелинейным самосогласованным потенциалом. Аналитически показано, что в системе существует нелинейная световая волна, распространяющаяся вдоль оптического слоя и локализованная в диэлектрических обкладках. Получены частоты локализации светового поля в данной структуре и определены условия их существования при различных характеристиках сред и границ их раздела. Показано, что локализация светового поля вдоль слоев может происходить при различных знаках нелинейного отклика границ раздела слоев трехслойной структуры, когда одна из них характеризуется фокусирующей нелинейностью, а другая – дефокусирующей.
ВВЕДЕНИЕ
Изучение многообразия свойств нелинейных поверхностных волн представляет интерес в связи с их широким применением в различных технических системах, основанных на волноводных свойствах многослойных гетероструктур [1, 2]. Нелинейные поверхностные волны оптического диапазона, распространяющиеся вдоль границ раздела нелинейных сред в слоистых структурах, в том числе и трехслойных (так называемых “сэндвичах”), аналитически исследовали многие авторы [3–7]. В данных работах искомое поле и его нормальные производные вблизи границ раздела слоистой структуры удовлетворяли условиям их непрерывности, что означало отсутствие взаимодействия волны с границей раздела как с плоским дефектом.
В [8] на основе нелинейного уравнения Шредингера с керровской нелинейностью учитывали взаимодействие возбуждений с двумя плоскими границами раздела трех нелинейных сред, характеризующимися одним параметром. Для теоретического описания локального взаимодействия нелинейных возбуждений с границами раздела слоев их моделировали короткодействующим потенциалом в нелинейном уравнении Шредингера, который в одномерном случае для трехслойной структуры записывается в виде:
где δ(x) – дельта-функция Дирака, U0 – интенсивность взаимодействия возбуждения с границей, 2а – расстояние между симметрично расположенными границами. Как отмечалось в [9], короткодействующий потенциал с одним параметром не всегда в полной мере позволяет проанализировать влияние характеристик границ раздела слоев на особенности локализации возбуждений. Нелинейные свойства внутри ультратонкой границы раздела слоев учитывали в [10–12] с использованием нелинейного потенциала $U \propto \delta (x){{\left| \psi \right|}^{2}}.$Следует отметить, что нелинейное уравнение Шредингера широко используется при теоретическом описании нелинейных волн в кристаллах [13–16]. Хорошо разработаны вопросы локализации возбуждений различной физической природы вблизи дефектов и границ раздела в нелинейных средах [17–22], а также локализации состояний на границе нелинейных и линейных сред в различных моделях [23–25], в том числе с учетом внутренних нелинейных свойств дефектов [10, 12, 20, 26, 27]. Изучение особенностей взаимодействия возбуждений с границей раздела сред при учете ее внутренних характеристик в многослойных структурах представляется важным, поскольку они могут быть выбраны в качестве управляющих параметров, контролирующих локализацию и волноводные свойства. Такие управляющие параметры необходимы для определения требуемых значений пропускных характеристик границ при определенных частотах в оптических устройствах, использующих волноводные свойства многослойных систем.
В настоящей работе предложено аналитическое описание локализации светового поля, возникающей в трехслойной структуре, в которой границы раздела представляют собой плоские дефекты с нелинейными свойствами, разделяющие нелинейную пластину конечной толщины и линейные полупространства [27]. Основной целью работы было нахождение профиля и частот локализации поля в явном аналитическом виде, а также условий их реализации в рассматриваемой системе.
УРАВНЕНИЯ МОДЕЛИ
В [8] было показано, что динамика светового поля в слоистой структуре описывается нелинейным уравнением Шредингера. Основываясь на данных результатах, рассмотрим систему плоскопараллельных чередующихся немагнитных узких и широких слоев. Пусть ось х направлена перпендикулярно плоскости слоев, а плоскость yz параллельна им. В плоско-поляризованной монохроматической волне, распространяющейся вдоль слоев, вектор напряженности электрического поля Е (он выбран параллельным орту ey = {0, 1, 0}) подчиняется уравнению Максвелла с показателем преломления, зависящим от координаты х в поперечном слоям направлении:
где n0 и n1 − линейные показатели преломления широких и узких слоев соответственно. В случае сред с эффектом Керра нелинейный показатель преломления зависит от квадрата амплитуды поля: α(x) и β(x) − коэффициенты керровской нелинейности сред в широких и узких слоях соответственно [20].Монохроматическая волна с волновым вектором k = exkz (где ex = {1, 0, 0} – орт) и частотой ω0 = = ckz/n0 в адиабатическом приближении представима в виде:
Теперь предлагается учесть нелинейность и внутри узких слоев по аналогии с [10, 11, 20, 26, 27]. Поэтому потенциал в (1) будет иметь вид:
Рассмотрим трехслойную структуру, в которой внутренний оптический слой толщиной 2а с нелинейностью керровского типа разделяет два диэлектрических (линейных) кристалла без эффекта Керра. Пусть границы раздела сред плоские и много меньше характерного масштаба локализации возмущений параметров среды, создаваемых ими, а ширина внутреннего слоя существенно больше ширины одной границ раздела слоев. Систему координат выберем так, чтобы средина нелинейного слоя проходила через начало координат. Границы раздела слоев лежали в плоскостях x = ±а перпендикулярно оси x. Линейные среды занимают полупространства $\left| x \right| > a,$ а нелинейный оптический слой расположен в области $\left| x \right| < a.$ Для описания интересующих нас свойств и условий их проявления используем одномерное стационарное нелинейное уравнение Шредингера, которому подчиняется поле ψ:
где ω – частота стационарных колебаний поля. Коэффициенты нелинейного уравнения Шредингера (1) определяются следующим образом:Потенциал, описывающий нелинейные свойства границ раздела, имеет вид:
(3)
$U\left( {x,{{{\left| \psi \right|}}^{2}}} \right) = {{F}_{1}}\left( {x + a,{{{\left| \psi \right|}}^{2}}} \right) + {{F}_{2}}\left( {x - a,{{{\left| \psi \right|}}^{2}}} \right),$Решение нелинейного уравнения Шредингера (2) с потенциалом (3) эквивалентно решению нелинейного уравнения Шредингера без потенциала с граничными условиями:
(5)
$\begin{gathered} \psi {\kern 1pt} '( \pm a + 0) - \psi {\kern 1pt} '( \pm a - 0) = \\ = 2m\{ {{U}_{j}} + {{W}_{j}}{{\left| {\psi ( \pm a)} \right|}^{2}}\} \psi ( \pm a). \\ \end{gathered} $В случае одной границы раздела в плоскости x = 0 при Uj = U0, Wj = 0 и а = 0 из (5) получается одно граничное условие, приведенное в [18] для линейного плоского дефекта, а при Uj = 0, Wj = W0 и а = 0 в случае нелинейного плоского дефекта − условие [10–12]. Для системы двух плоскопараллельных дефектов с линейным взаимодействием при Uj = U0, Wj = 0 и а ≠ 0 из (5) получаются граничные условия [8]. При Uj = 0, Wj ≠ 0, и а ≠ 0 из (5) имеют место граничные условия, использованные в [27] для двух плоскопараллельных дефектов с преобладающим нелинейным откликом. В настоящей работе проанализирован только случай одинаковых значений соответствующих характеристик слоев и их границ раздела.
Если частота волны лежит в диапазоне Ω0 < ω < < min{Ωj}, то нелинейное уравнение Шредингера (2) имеет решение:
(6)
$\psi (x) = \left\{ \begin{gathered} {{\psi }_{{0j}}}{\text{exp(}} \mp {{q}_{j}}(x \mp a)),\,\,\,\,\left| x \right| > a, \hfill \\ {{A}_{s}}{\text{sn(}}{{q}_{s}}(x - {{x}_{s}}),k),\,\,\,\,\,\,\,\,\left| x \right| < a, \hfill \\ \end{gathered} \right.$Подстановка (6) в (4) позволяет получить выражения для амплитуд затухающего в линейных полупространствах поля (т.е. амплитуд колебаний поля на границах раздела):
(7)
${{\psi }_{{0j}}} = {{ \pm {{q}_{s}}{\text{sn(}}{{q}_{s}}(a \mp {{x}_{s}}),k)} \mathord{\left/ {\vphantom {{ \pm {{q}_{s}}{\text{sn(}}{{q}_{s}}(a \mp {{x}_{s}}),k)} {{{{(mg)}}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}}}} \right. \kern-0em} {{{{(mg)}}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}}}.$Подстановка (6) в (5) с учетом (7) приводит к паре дисперсионных соотношений, определяющих зависимость частоты нелинейной волны от параметров слоев и их границ раздела:
(8)
$\begin{gathered} {{Q}_{j}}(\omega ) + {{D}_{j}}(\omega ,{{U}_{j}},{{W}_{j}}) = 0, \\ {{Q}_{j}}(\omega ) = {{q}_{s}}{\text{cn(}}{{q}_{s}}(a \mp {{x}_{s}}),k) \times \\ \times \,\,{{dn{\text{(}}{{q}_{s}}(a \mp {{x}_{s}}),k)} \mathord{\left/ {\vphantom {{dn{\text{(}}{{q}_{s}}(a \mp {{x}_{s}}),k)} {{\text{sn(}}{{q}_{s}}(a \mp {{x}_{s}}),k) - {{q}_{j}}}}} \right. \kern-0em} {{\text{sn(}}{{q}_{s}}(a \mp {{x}_{s}}),k) - {{q}_{j}}}}, \\ {{D}_{j}}(\omega ,{{U}_{j}},{{W}_{j}}) = 2m\{ {{U}_{j}} + {{V}_{j}}q_{s}^{2}{\text{s}}{{{\text{n}}}^{{\text{2}}}}{\text{(}}{{q}_{s}}(a \mp {{x}_{s}}),k)\} , \\ \end{gathered} $РЕЗУЛЬТАТЫ И ИХ ОБСУЖДЕНИЕ
В настоящей работе интересны эффекты несимметричной локализации, возникающие вследствие различия параметров слоев и их границ раздела, когда Ω1 ≠ Ω2, U1 ≠ U2, W1 ≠ W2. Будем рассматривать состояния специального вида, для которых xs = 0. Для таких состояний из (7) следует, что они соответствуют противофазным колебаниям: ψ01 = −ψ02. В длинноволновом приближении qsа $ \ll $ 1, когда частота близка к краю спектра при $\left| {\omega - {{\Omega }_{0}}} \right| \ll {{(1 + {{k}^{2}})} \mathord{\left/ {\vphantom {{(1 + {{k}^{2}})} {2m{{a}^{2}}}}} \right. \kern-0em} {2m{{a}^{2}}}},$ из дисперсионных уравнений (8) можно получить выражение для частоты в явном виде:
(9)
$\begin{gathered} \omega = {{\Omega }_{0}} + (1 + {{k}^{2}}) \times \\ \times \,\,{{\{ - 1 \pm {{{[{{1 - 4{{a}_{j}}{{c}_{j}}} \mathord{\left/ {\vphantom {{1 - 4{{a}_{j}}{{c}_{j}}} {{{{(1 + {{k}^{2}})}}^{2}}}}} \right. \kern-0em} {{{{(1 + {{k}^{2}})}}^{2}}}}]}}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}\} } \mathord{\left/ {\vphantom {{\{ - 1 \pm {{{[{{1 - 4{{a}_{j}}{{c}_{j}}} \mathord{\left/ {\vphantom {{1 - 4{{a}_{j}}{{c}_{j}}} {{{{(1 + {{k}^{2}})}}^{2}}}}} \right. \kern-0em} {{{{(1 + {{k}^{2}})}}^{2}}}}]}}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}\} } {4m{{a}_{j}}}}} \right. \kern-0em} {4m{{a}_{j}}}}, \\ \end{gathered} $Следует отметить, что выражение (9) справедливо при условии, что параметры слоев и их границ раздела связаны соотношением:
Проанализируем далее дисперсионные соотношения (8) в других частных случаях.
1. Слабо нелинейный отклик границ, когда в пределе можно положить W0 = 0. В этом случае частота локализованной около внутреннего слоя нелинейной волны в длинноволновом приближении определяется выражением:
(10)
$\omega = {{\Omega }_{j}} + {{{{{({1 \mathord{\left/ {\vphantom {1 a}} \right. \kern-0em} a} + 2m{{U}_{j}})}}^{2}}} \mathord{\left/ {\vphantom {{{{{({1 \mathord{\left/ {\vphantom {1 a}} \right. \kern-0em} a} + 2m{{U}_{j}})}}^{2}}} {2m}}} \right. \kern-0em} {2m}}.$(11)
${{\psi }_{{0j}}} = \pm \frac{{a\{ 2m({{\Omega }_{j}} - {{\Omega }_{0}}) - {{{({1 \mathord{\left/ {\vphantom {1 a}} \right. \kern-0em} a} + 2m{{U}_{j}})}}^{2}}\} }}{{(1 + {{k}^{2}})\sqrt {mg} }}.$2. Преобладающий нелинейный отклик границ, когда в пределе можно положить U0 = 0. В этом случае частота локализованной около внутреннего слоя нелинейной волны в длинноволновом приближении определяется выражением:
(12)
$\omega = {{\Omega }_{0}} + \frac{{1 + {{k}^{2}}}}{{2a}}\sqrt {\frac{{{{\Omega }_{2}} - {{\Omega }_{1}}}}{{ma({{V}_{2}} - {{V}_{1}})}}} .$(13)
${{\psi }_{{0j}}} = \pm \sqrt {\frac{{{{\Omega }_{2}} - {{\Omega }_{1}}}}{{ga({{V}_{2}} - {{V}_{1}})}}} .$3. Теперь рассмотрим случай, когда параметры слоев одинаковы Ω1 = Ω2, а характеристики границ по-прежнему различные: U1 ≠ U2, W1 ≠ W2. В этом случае частота локализованной около внутреннего слоя нелинейной волны в длинноволновом приближении определяется выражением:
(14)
$\omega = {{\Omega }_{0}} + \frac{{1 + {{k}^{2}}}}{{2a}}\sqrt {\frac{{{{U}_{1}} - {{U}_{2}}}}{{m({{V}_{2}} - {{V}_{1}})}}} .$ЗАКЛЮЧЕНИЕ
Таким образом, установлено, что в трехслойной структуре типа сэндвича, в которой внутренний оптический нелинейный слой с дефокусировкой находится между диэлектрическими полупространствами, разделенными границами раздела с нелинейными свойствами, вдоль слоев может распространяться нелинейная локализованная волна возмущения напряженности эклектического поля. Найдены частоты длинноволновых колебаний в явном аналитическом виде, и проанализированы условия их существования.
Описано новое свойство, заключающееся в том, что локализация поля вдоль слоев может происходить при различных знаках нелинейного отклика обеих границ раздела слоев: одна из них может обладать фокусирующей нелинейностью, а другая – дефокусирующей. Также определена частота, при которой локализация волны будет происходить, только если границы раздела будут обладать выраженным нелинейным откликом (в случае оптической системы – сильным эффектом Керра с нелинейными показателями преломления, значения которых могут быть различны).
Представляется важным то, что локализация светового поля возможна при различных нелинейных показателях преломления внутри тонких граничных слоев, разделяющих широкие слои трехслойной структуры. Управляя их значениями в ходе создания оптической системы типа сэндвича, можно получать заданные частоты локализации светового поля вдоль слов. Полученные результаты способствуют развитию и совершенствованию технологий разработки оптических систем, основанных на многослойных структурах [1, 2], использующих свойства границ раздела слоев.
Список литературы
Zhang D., Li Z., Hu W., Cheng B. // Appl. Phys. Lett. 1995. V. 67. P. 2431. https://doi.org/10.1063/1.114597
Naim B.A. // Chin. J. Phys. 2017. V. 55. P. 2384. https://doi.org/10.1016/j.cjph.2017.10.008
Ахмедиев Н.Н. // ЖЭТФ. 1982. Т. 83. С. 545.
Boardman A.D., Shabat M.M., Wallis R.F. // J. Phys. D. 1991. V. 24. P. 1702. https://doi.org/10.1088/0022-3727/24/10/002
Ashour H.S., Assa’d A.I. // Turk. J. Phys. 2012. V. 36. P. 207. https://doi.org/10.3906/fiz-1106-8
Михалаке Д., Назмитдинов Р.Г., Федянин В.К. // Физика элементарных частиц и атомного ядра. 1989. Т. 20. № 1. С. 198.
Коровай О.В., Хаджи П.И. // ФТТ. 2003. Т. 45. С. 364. https://doi.org/10.1134/1.1553548
Герасимчук И.В., Ковалев А.С. // ФНТ. 2000. Т. 26. № 8. С. 799. https://doi.org/10.1063/1.1289129
Савотченко С.Е. // ЖЭТФ. 2018. Т. 153. № 2. С. 339. https://doi.org/10.1134/S1063776118020061
Gerasimchuk I.V., Gorbach P.K., Dovhopolyi P.P. // Ukr. J. Phys. 2012. V. 57. P. 678.
Герасимчук И.В. // ЖЭТФ. 2015. Т. 121. № 4. С. 596. https://doi.org/10.1134/S1063776115100076
Savotchenko S.E. // Mod. Phys. Lett. B. 2018. V. 32. № 10. P. 1 850 120. https://doi.org/10.1142/S0217984918501208
Косевич А.М., Ковалев А.С. Введение в нелинейную физическую механику. Киев: Наукова думка, 1989. 304 с.
Kivshar Y.S., Agrawal G.P. Optical Solitons: From Fibers to Photonic Crystals. San Diego: Academic Press, 2003. 540 p.
Fibich G., Sivan Y., Weinstein M.I. // Phys. D. 2006. V. 217. P. 31. https://doi.org/10.1016/j.physd.2006.03.009
Murali R., Senthilnathan K., Porsezian K. // J. Phys. B. 2008. V. 41. P. 025401. https://doi.org/10.1088/0953-4075/41/2/025401
Kivshar Y.S., Kosevich A.M., Chubykalo O.A. // Phys. Rev. A. 1990. V. 41. P. 1677. https://doi.org/10.1103/PhysRevA.41.1677
Богдан М.М., Герасимчук И.В., Ковалев А.С. // Физика низких температур. 1997. Т. 23. С. 197. https://doi.org/10.1063/1.593346
Савотченко С.Е. // Изв. вузов. Физика. 2004. Т. 47. № 5. С. 79. https://doi.org/10.1023/B:RUPJ.0000046330.92744.73
Sukhorukov A.A., Kivshar Yu.S. // Phys. Rev. Lett. 2001. V. 87. P. 083901. https://doi.org/10.1103/PhysRevLett.87.083901
Kartashov Y.V., Malomed B.A., Torner L. // Rev. Mod. Phys. 2011. V. 83. P. 247. https://doi.org/10.1103/RevModPhys.83.247
Савотченко С.Е. // Конденсированные среды и межфазные границы. 2017. № 2. С. 291. https://doi.org/10.17308/kcmf.2017.19/205
Савотченко С.Е. // Конденсированные среды и межфазные границы. 2017. Т. 19. № 4. С. 567. https://doi.org/10.17308/kcmf.2017.19/238
Савотченко С.Е. // ЖТФ. 2017. Т. 62. № 12. С. 1776. https://doi.org/10.1134/S1063784217120210
Савотченко С.Е. // Вестн. Воронеж. гос. ун-та. Сер. Физика. Математика. 2018. № 1. С. 44.
Савотченко С.Е. // Письма в ЖЭТФ. 2018. Т. 107. № 8. С.481. https://doi.org/10.7868/S0370274X18080027
Савотченко С.Е. // Нелинейный мир. 2018. № 3. С. 25.
Дополнительные материалы отсутствуют.
Инструменты
Поверхность. Рентгеновские, синхротронные и нейтронные исследования