Известия РАН. Серия физическая, 2019, T. 83, № 1, стр. 32-36
Стабилизация оптического импульса, распространяющегося в режиме ионизации и вынужденного комбинационного саморассеянияВ. А. Халяпин
В. А. Халяпин *
БФУ им. И. Канта, ФГБOУ ВО “Калининградский государственный технический университет”
Калининград, Россия
* E-mail: slavasxi@gmail.com
Аннотация
Аналитически исследована динамика солитоноподобных импульсов, распространяющихся в режиме туннельной ионизации и вынужденного комбинационного саморассеяния. Показано, что при определенных условиях возможна взаимная компенсация этих эффектов, которая приводит к стабилизации сигнала. Получено явное решение для длительности сигнала, описывающее данный процесс.
ВВЕДЕНИЕ
В последние годы значительное внимание привлекают новые фотонно-кристаллические волокна (ФКВ) типа “кагоме” [1–3], которые имеют широкую полосу пропускания, включающую всю видимую и инфракрасную области спектра и низкую аномальную дисперсию групповой скорости для широкого спектрального диапазона. Дисперсия таких волокон складывается из дисперсии самого волокна и газа, что позволяет изменять ее за счет давления газа.
В работе [4] было показано, что ионизация приводит к тому, что спектр импульса смещается в сторону высоких частот. Это вызвано тем, что возникающие за счет ионизации свободные электроны вносят отрицательный вклад в показатель преломления [4–12]. Этот эффект противоположен вынужденному комбинационному саморассеянию (ВКС), которое вызывает красное смещение спектра сигнала [13–19]. Уравнение, описывающее распространение импульсов с учетом этих эффектов в ФКВ типа “кагоме” было получено в работах [8, 9]. С помощью теории возмущений [20–22] авторы показали, что при определенных условиях центральная частота солитона может оставаться такой же, какой была на входе в волокно. Если компенсация неполная, то центральная частота смещается по линейному закону в красную или фиолетовую область спектра. В настоящей работе с помощью метода моментов [17, 23, 24] обобщены результаты [8, 9]. Показано, что в рамках приближения нулевого порога ионизации и пренебрежении поглощением стабилизация сигнала будет происходить не только при точной компенсации ВКС и ионизации. Также нами проведен учет поглощения, которое приводит к уменьшению роли ионизации и нарушению стабилизации параметров сигнала.
МЕТОД МОМЕНТОВ И СТАБИЛИЗАЦИЯ ДЛИТЕЛЬНОСТИ СИГНАЛА
Динамика световых импульсов, распространяющихся в ФКВ, описывается уравнением [8, 9]
(1)
$\begin{gathered} \frac{{\partial \psi }}{{\partial z}} + \frac{{i{{\beta }_{2}}}}{2}\frac{{{{\partial }^{2}}\psi }}{{\partial {{\tau }^{2}}}} - \frac{{{{\beta }_{3}}}}{6}\frac{{{{\partial }^{3}}\psi }}{{\partial {{\tau }^{3}}}} - i\gamma \psi {{\left| \psi \right|}^{2}} + \frac{\gamma }{\omega }\frac{\partial }{{\partial \tau }}\left( {\psi {{{\left| \psi \right|}}^{2}}} \right) + \\ + \,\,i\gamma {{T}_{R}}\psi \frac{{\partial {{{\left| \psi \right|}}^{2}}}}{{\partial \tau }} + i\eta \psi \int\limits_{ - \infty }^\tau {{{{\left| \psi \right|}}^{2}}d\tau } + \frac{\beta }{2}\psi = 0, \\ \end{gathered} $Медленно меняющаяся огибающая $\psi $ связана с электрическим полем импульса E соотношением
(2)
$E(z,\tau ) = \frac{1}{2}\psi (z,\tau )\exp \left[ { - i\left( {\omega t - kz} \right)} \right] + {\text{c}}{\text{.c}}{\text{.}}$(3)
$\begin{gathered} \psi = B\operatorname{sech} \left( {\frac{{\tau - T}}{{{{\tau }_{p}}}}} \right) \times \\ \times \,\,\exp \left[ {i\left( {\phi + \Omega \left( {\tau - T} \right) - C\frac{{{{{\left( {\tau - T} \right)}}^{2}}}}{{2\tau _{p}^{2}}}} \right)} \right], \\ \end{gathered} $(5)
$\tilde {C} = \frac{i}{{2E}}\int\limits_{ - \infty }^\infty {\left( {\tau - T} \right)\left( {\psi {\text{*}}\frac{{\partial \psi }}{{\partial \tau }} - \psi \frac{{\partial \psi {\text{*}}}}{{\partial \tau }}} \right)d\tau } ,$(6)
${{\sigma }^{2}} = \frac{1}{E}\int\limits_{ - \infty }^\infty {{{{\left( {\tau - T} \right)}}^{2}}{{{\left| \psi \right|}}^{2}}d\tau } ,$(7)
$T = \frac{1}{E}\int\limits_{ - \infty }^\infty {\left( \tau \right){{{\left| \psi \right|}}^{2}}d\tau } ,$(8)
$\Omega = - \frac{i}{{2E}}\int\limits_{ - \infty }^\infty {\left( {\psi {\text{*}}\frac{{\partial \psi }}{{\partial \tau }} - \psi \frac{{\partial \psi {\text{*}}}}{{\partial \tau }}} \right)d\tau } ,$(9)
$\frac{{dE}}{{dz}} = - \beta \int\limits_{ - \infty }^\infty {{{{\left| \psi \right|}}^{2}}d\tau } ,$(10)
$\begin{gathered} \frac{{d\Omega }}{{dz}} = \frac{{{{T}_{R}}\gamma }}{E}\int\limits_{ - \infty }^\infty {{{{\left( {\frac{{\partial {{{\left| \psi \right|}}^{2}}}}{{\partial \tau }}} \right)}}^{2}}d\tau + } \\ + \,\,\frac{{i\gamma }}{{E\omega }}\int\limits_{ - \infty }^\infty {\left( {\psi {\text{*}}\frac{{\partial \psi }}{{\partial \tau }} - \psi \frac{{\partial \psi {\text{*}}}}{{\partial \tau }}} \right)\frac{{\partial {{{\left| \psi \right|}}^{2}}}}{{\partial \tau }}} d\tau - \frac{\eta }{E}\int\limits_{ - \infty }^\infty {{{{\left| \psi \right|}}^{4}}d\tau } , \\ \end{gathered} $(11)
$\frac{{dT}}{{dz}} = - {{\beta }_{2}}\Omega + \frac{{{{\beta }_{3}}}}{{2E}}\int\limits_{ - \infty }^\infty {{{{\left| {\frac{{\partial \psi }}{{\partial \tau }}} \right|}}^{2}}} d\tau + \frac{{3\gamma }}{{2E\omega }}\int\limits_{ - \infty }^\infty {{{{\left| \psi \right|}}^{4}}} d\tau ,$(12)
$\frac{{d{{\sigma }^{2}}}}{{dz}} = 2{{\beta }_{2}}\tilde {C} + \frac{{{{\beta }_{3}}}}{E}\int\limits_{ - \infty }^\infty {\left( {\tau - T} \right){{{\left| {\frac{{\partial \psi }}{{\partial t}}} \right|}}^{2}}} d\tau ,$(13)
$\begin{gathered} \frac{{d\tilde {C}}}{{dz}} = \Omega \frac{{dT}}{{dz}} + \frac{{{{\beta }_{2}}}}{E}\int\limits_{ - \infty }^\infty {{{{\left| {\frac{{\partial \psi }}{{\partial t}}} \right|}}^{2}}} d\tau + \frac{{i{{\beta }_{3}}}}{{4E}}\int\limits_{ - \infty }^\infty {\left( {\frac{{{{\partial }^{2}}\psi }}{{\partial {{\tau }^{2}}}}\frac{{\partial \psi {\text{*}}}}{{\partial \tau }} - \frac{{{{\partial }^{2}}\psi {\text{*}}}}{{\partial {{\tau }^{2}}}}\frac{{\partial \psi }}{{\partial \tau }}} \right)} d\tau + \frac{\gamma }{{2E}}\int\limits_{ - \infty }^\infty {{{{\left| \psi \right|}}^{4}}} d\tau - \\ - \,\,\frac{{i\gamma }}{{E\omega }}\int\limits_{ - \infty }^\infty {\left( {\tau - T} \right)\frac{{\partial {{{\left| \psi \right|}}^{2}}}}{{\partial \tau }}\left( {\psi {\text{*}}\frac{{\partial \psi }}{{\partial \tau }} - \psi \frac{{\partial \psi {\text{*}}}}{{\partial \tau }}} \right)} d\tau - \frac{{i\gamma }}{{2E\omega }}\int\limits_0^\infty {{{{\left| \psi \right|}}^{2}}\left( {\psi \frac{{\partial \psi {\text{*}}}}{{\partial \tau }} - \psi {\text{*}}\frac{{\partial \psi }}{{\partial \tau }}} \right)} d\tau - \frac{{\gamma {{T}_{R}}}}{E}{{\int\limits_{ - \infty }^\infty {\left( {\tau - T} \right)\left( {\frac{{\partial {{{\left| \psi \right|}}^{2}}}}{{\partial \tau }}} \right)} }^{2}}d\tau . \\ \end{gathered} $Уравнение, определяющее динамику фазы ϕ можно записать в неявном виде
(14)
$\begin{gathered} \int\limits_{ - \infty }^\infty {\left( {\psi {\text{*}}\frac{{\partial \psi }}{{\partial z}} - \psi \frac{{\partial \psi {\text{*}}}}{{\partial z}}} \right.} - i{{\beta }_{2}}{{\left| {\frac{{\partial \psi }}{{\partial \tau }}} \right|}^{2}} - \frac{{{{\beta }_{3}}}}{6}\left( {\psi {\text{*}}\frac{{{{\partial }^{3}}\psi }}{{\partial {{\tau }^{3}}}} - \psi \frac{{{{\partial }^{3}}\psi {\text{*}}}}{{\partial {{\tau }^{3}}}}} \right) - 2i{{\left| \psi \right|}^{4}} + \\ + \,\,2i\eta {{\left| \psi \right|}^{2}}\int\limits_{ - \infty }^\tau {{{{\left| \psi \right|}}^{2}}d\tau {\text{'}}} + \left. {\frac{\gamma }{\omega }\left( {\psi {\text{*}}\frac{{{{\partial }^{3}}\psi }}{{\partial {{\tau }^{3}}}} - \psi \frac{{{{\partial }^{3}}\psi {\text{*}}}}{{\partial {{\tau }^{3}}}}} \right)} \right)d\tau . \\ \end{gathered} $Подставляя пробную функцию (3) в (9)–(13) и интегрируя по τ, получаем систему уравнений для параметров импульса
(16)
$\frac{{d\Omega }}{{dz}} = \frac{{4\gamma E}}{{15\tau _{p}^{3}}}\left( {{{T}_{R}} - \frac{{5C}}{{4\omega }}} \right) - \frac{{\eta E}}{{3{{\tau }_{p}}}},$(17)
$\frac{{dT}}{{dz}} = - {{\beta }_{2}}\Omega + \frac{{{{\beta }_{3}}}}{2}\left( {\frac{{\left( {1 + {{\pi }^{2}}{{{{C}^{2}}} \mathord{\left/ {\vphantom {{{{C}^{2}}} 4}} \right. \kern-0em} 4}} \right)}}{{3\tau _{p}^{2}}} + \Omega _{{}}^{2}} \right) + \frac{{\gamma E}}{{2\omega {{\tau }_{p}}}},$(18)
$\frac{{d\tau _{p}^{2}}}{{dz}} = 2C\left( {{{\beta }_{2}} - {{\beta }_{3}}{{\Omega }_{{}}}} \right),$(19)
$\frac{{dС }}{{dz}} = \left( {{{\beta }_{2}} - {{\beta }_{3}}\Omega } \right)\left( {\frac{{{4 \mathord{\left/ {\vphantom {4 {{{\pi }^{2}}}}} \right. \kern-0em} {{{\pi }^{2}}}} + {{С }^{2}}}}{{\tau _{p}^{2}}}} \right) + \frac{{2\gamma E}}{{{{\pi }^{2}}{{\tau }_{p}}}}\left( {1 - \frac{\Omega }{\omega }} \right),$Будем искать эволюцию параметров ${{\tau }_{p}}$ и $\Omega $ с помощью метода последовательных приближений. Используя (18) и полагая производную от C нулю, из (19) получаем
(20)
$\Omega = \frac{{\left| {{{\beta }_{2}}} \right|}}{{{{\beta }_{3}}}}\left( {\frac{{\gamma {{\tau }_{p}}E}}{{2\left| {{{\beta }_{2}}} \right|}} - 1} \right).$Поскольку на входе в среду смещение центральной частоты равно нулю (${{\Omega }_{0}} = 0$), то из (20) следует, что начальная энергия импульса ${{E}_{0}}$ соответствует известному значению для солитона нелинейного уравнения Шредингера ${{E}_{0}} = {{2\left| {{{\beta }_{2}}} \right|} \mathord{\left/ {\vphantom {{2\left| {{{\beta }_{2}}} \right|} {\gamma {{\tau }_{0}}}}} \right. \kern-0em} {\gamma {{\tau }_{0}}}}.$ Здесь ${{\tau }_{0}}$ – начальное значение длительности сигнала. Подставляя (20) в (16), учитывая (15) и пренебрегая производной нелинейности, получаем уравнение для безразмерной длительности сигнала $u = {{{{\tau }_{p}}} \mathord{\left/ {\vphantom {{{{\tau }_{p}}} {{{\tau }_{0}}}}} \right. \kern-0em} {{{\tau }_{0}}}}$
где $a = {{8{{\beta }_{3}}{{T}_{R}}} \mathord{\left/ {\vphantom {{8{{\beta }_{3}}{{T}_{R}}} {15\tau _{0}^{4}}}} \right. \kern-0em} {15\tau _{0}^{4}}}$ и $b = {{2\eta {{\beta }_{3}}} \mathord{\left/ {\vphantom {{2\eta {{\beta }_{3}}} {3\gamma \tau _{0}^{2}}}} \right. \kern-0em} {3\gamma \tau _{0}^{2}}}$ – обратные характерные длины ВКС и фотоионизации соответственно. Если мы пренебрежем поглощением ($\beta \approx 0$), то уравнение (21) имеет решениеВыражение (22) описывает динамику длительности импульса вплоть до точки стабилизации
Авторы работ [8, 9] с помощью теории возмущений показали, что если $a = b,$ то импульс сразу стабилизируется. Однако из анализа (22) следует, что в рамках используемого в [8, 9] приближения формально режим ионизации будет наступать всегда. Рассмотрим теперь вклад поглощения, которое сопровождает ионизацию. Учет последнего слагаемого в (21) приводит к тому, что при условии, когда правая часть (21) не имеет корней ${{b}^{2}} < 4a\beta $ решение (21) имеет вид:
(24)
$\begin{gathered} 2\beta z = \ln \sqrt {\frac{{\beta {{u}^{4}} - b{{u}^{2}} + a}}{{\beta - b + a}}} + \frac{b}{{\sqrt {4a\beta - {{b}^{2}}} }} \times \\ \times \,\,\left( {{\text{arctg}}\frac{{2\beta {{u}^{2}} - b}}{{\sqrt {4a\beta - {{b}^{2}}} }} - {\text{arctg}}\frac{{2\beta - b}}{{\sqrt {4a\beta - {{b}^{2}}} }}} \right). \\ \end{gathered} $Это решение соответствует случаю отсутствия точки стабилизации и описывает увеличение длительности сигнала до тех пор, пока его интенсивность не приблизится к пороговому значению. Рассмотрим случай, когда ${{b}^{2}} > 4a\beta .$ Тогда решение (21) принимает вид:
(25)
$2\beta z = \frac{{u_{1}^{2}}}{{u_{1}^{2} - u_{2}^{2}}}\ln \frac{{{{u}^{2}} - u_{1}^{2}}}{{1 - u_{1}^{2}}} - \frac{{u_{2}^{2}}}{{u_{1}^{2} - u_{2}^{2}}}\ln \frac{{{{u}^{2}} - u_{2}^{2}}}{{1 - u_{2}^{2}}}.$Здесь ${{u}_{1}}$ и ${{u}_{2}}$ – нули правой части (21)
Из рис. 1 следует, что ${{u}_{2}}$ является точкой стабилизации сигнала. Если выполняется условие ${{b}^{2}} \gg 4a\beta ,$ то ${{u}_{2}}$ переходит в выражение для относительной длительности стаблизации (23), которое получено без учета поглощения. Если перейти непосредственно к длительности сигнала, то из (26) получаем
(27)
${{\tau }_{{1,2}}} = \sqrt {{{\left( {\frac{{2\eta {{\beta }_{3}}}}{{3\gamma }} \pm \sqrt {{{{\left( {\frac{{2\eta {{\beta }_{3}}}}{{3\gamma }}} \right)}}^{2}} - \frac{{16\beta {{\beta }_{3}}{{T}_{R}}}}{{15}}} } \right)} \mathord{\left/ {\vphantom {{\left( {\frac{{2\eta {{\beta }_{3}}}}{{3\gamma }} \pm \sqrt {{{{\left( {\frac{{2\eta {{\beta }_{3}}}}{{3\gamma }}} \right)}}^{2}} - \frac{{16\beta {{\beta }_{3}}{{T}_{R}}}}{{15}}} } \right)} {2\beta }}} \right. \kern-0em} {2\beta }}} .$Если ${{\tau }_{0}} > {{\tau }_{1}},$ то длительность импульса сразу начинает увеличивается. При условии ${{\tau }_{0}} < {{\tau }_{1}}$ длительность стремиться к значению ${{\tau }_{2}}.$ Раcсмотрим динамику импульса, начальная длительность которого лежит вблизи точки стабилизации ${{\tau }_{2}}$ и удовлетворяет условию ${{\tau }_{0}} < {{\tau }_{1}}.$ Используя разложение в ряд Тейлора для (21) или непосредственно из (25) можно получить соответствующее приближенное решение
(28)
$u = {{u}_{2}} + \left( {1 - {{u}_{2}}} \right)\exp \left( { - \frac{{4\beta \sqrt {{{b}^{2}} - 4a\beta } }}{{b - \sqrt {{{b}^{2}} - 4a\beta } }}z} \right).$Как видно из рис. 1, данное приближенное решение описывает процесс стабилизации сигнала. Правая часть (21) имеет один корень (если ${{b}^{2}} = 4a\beta $) при относительной длительности
Решение (21) тогда имеет вид:
(30)
$\ln \frac{{{{u}^{2}} - u_{{st}}^{2}}}{{1 - u_{{st}}^{2}}} + u_{s}^{2}\left( {\frac{1}{{1 - u_{{st}}^{2}}} - \frac{1}{{{{u}^{2}} - u_{{st}}^{2}}}} \right) = 2\beta z.$Явное приближенное решение, справедливое вблизи корня, можно получить, воспользовавшись разложением в ряд Тейлора для (21) или непосредственно из (30), пренебрегая вкладом первого слагаемого по сравнению с остальными. В результате получаем
(31)
$u = {{u}_{{st}}}\left( {1 + \frac{{1 - {{u}_{{st}}}}}{{{{u}_{{st}}} - \left( {1 - {{u}_{{st}}}} \right)4\beta z}}} \right).$Из (31) следует, что если ${{u}_{{st}}} > 1,$ то относительная длительность импульса (начальное значение которой равно единице) будет стремиться к величине ${{u}_{{st}}}$ (рис. 2). В противном случае длительность сигнала будет сразу увеличиваться, минуя режим стабилизации.
ЗАКЛЮЧЕНИЕ
В работе рассмотрен случай распространения светового импульса в режиме ВКС и туннельной ионизации с учетом поглощения, сопровождающего ионизацию. Поскольку процесс ионизации ведет к синему сдвигу спектра импульса, а ВКС вызывает красное смещение, то они могут компенсировать друг друга и стабилизировать динамику импульса. Вслед за авторами работ [8, 9] было использовано приближение, при котором можно пренебречь пороговым значением интенсивности ионизации по сравнению с интенсивностью сигнала. На основе метода моментов проведено аналитическое исследование, позволяющее предсказать режим динамики параметров сигнала. Найдены условия, при которых длительность импульса будет стремиться к стабильному значению, и получено соответствующее аналитическое решение. Показано, что учет поглощения может играть принципиальную роль в ограничении стабилизации сигнала.
Работа выполнена при финансовой поддержке РНФ (грант № 17-11-01157).
Список литературы
Couny F., Benabid F., Light P.S. // Opt. Lett. 2006. V. 31. P. 3574.
Couny F., Roberts P.J., Birks T.A., Benabid F. // Opt. Express 2008. V. 16. P. 20626.
Couny F., Benabid F., Roberts P.J. et al. // Science. 2007. V. 318. P. 1118.
Serebryannikov E.E., Zheltikov A.M. // Phys. Rev. A. 2007. V. 76. P. 013820.
Wood W.M., Siders C.W., Downer M.C. // Phys. Rev. Lett. 1991. V. 67. P. 3523.
Wilks S.C., Dawson J.M., Mori W.B. // Phys. Rev. Lett. 1988. V. 61. P. 337.
Yablonovitch E. // Phys. Rev. A. 1974. V. 10. P. 1888.
Hölzer P., Chang W., Travers J.C. et al. // Phys. Rev. Lett. 2011. V. 107. P. 203902.
Hölzer P., Chang W., Travers J.C. et al. // Phys. Rev. Lett. 2011. V. 107. P. 203901.
Saleh M.F., Biancalana F. // Phys. Rev. A. 2011. V. 84. P. 063838.
Facão M., Carvalho M.I., Almeida P. // Phys. Rev. A. 2013. V. 87. P. 063803.
Facão M., Carvalho M.I. // Appl. Phys. B. 2014. V. 116. P. 353.
Dianov E.M., Karasik A.Ya., Mamyshev P.V. et al. // JETP Lett. 1985. V. 41. P. 294.
Mitschke F.M., Mollenauer L.F. // Opt. Lett. 1986. V. 11. P. 659.
Gordon J.P. // Opt. Lett. 1986. V. 11. P. 662.
Agrawal G.P. Nonlinear Fiber Optics. NewYork: Academic, 2007.
Santhanam J., Agraval G. // Opt. Commun. 2003. V. 222. P. 413.
Bugay A.N., Khalyapin V.A. // Phys. Lett. A. 2017. V. 381. № 4. P. 399.
Bugay A.N., Khalyapin V.A. // Optics and Spectroscopy. 2017. V. 123. № 2. P. 181.
Serkin V.N., Vysloukh V.A. // Wave Phenomena, Techn. Digest, Opt. Soc. of Am. 1993. V. 15. P. 236.
Serkin V.N., Vysloukh V.A., Taylor J.R. // Electron. Lett. 1993. V. 29. P. 12.
Serkin V.N., Belyaeva T.L., Corro G.H., Granados M.A. // Quantum Electron. 2003. V. 33. P. 456.
Lefrancois S., Husko C., Blanco-Redondo A., Eggleton B.J. // J. Opt. Soc. Am. B. 2015. V. 32. P. 218.
Tsoy E.N., Sterke C.M. // J. Opt. Soc. Am. B. 2006. V. 23. P. 2425.
Kozlov S.A., Sazonov S.V. // JETP. 1997. V. 84. P. 221.
Sazonov S.V., Khalyapin V.A. // Opt. Spectr. 2003. V. 95. P. 401.
Дополнительные материалы отсутствуют.
Инструменты
Известия РАН. Серия физическая