Известия РАН. Серия физическая, 2019, T. 83, № 3, стр. 416-419

Анализ сигнала сверхбыстрого оптического эффекта керра с учетом корреляции вращательных откликов молекул в жидкости

В. Г. Никифоров 1*, Д. К. Жарков 1, А. Г. Шмелев 1, А. В. Леонтьев 1, В. С. Лобков 1

1 Казанский физико-технический институт им. Е.К. Завойского – обособленное структурное подразделение Федерального государственного бюджетного учреждения науки “Федеральный исследовательский центр “Казанский научный центр Российской академии наук”
Казань, Россия

* E-mail: vgnik@mail.ru

Поступила в редакцию 03.09.2018
После доработки 10.09.2018
Принята к публикации 22.10.2018

Полный текст (PDF)

Аннотация

Представлен анализ сигнала сверхбыстрого оптического эффекта Керра (ОЭК) в жидкостях: толуол C7H8, нитробензол C6H5NO2 и ацетонитрил C2H3N. Показано, что учет корреляции между ориентационным и либрационным откликами позволяет определить время нарастания ориентационного отклика, связанного с инерционным характером молекулярных вращений на временном интервале до 300–500 фс после действия лазерного импульса, и увеличивает точность определения вида функции рамановской спектральной плотности.

ВВЕДЕНИЕ

Техника фемтосекундной поляризационной спектроскопии с регистраций сверхбыстрого оптического эффекта Керра (ОЭК) позволяет исследовать оптические отклики в изотропных средах, связанные с колебательной и вращательной динамикой молекул, в субпикосекундном диапазоне, что соответствует спектральной области 0–600 см–1 [15]. На сегодняшний день с технической точки зрения она является наиболее доступным методом изучения низкочастотных молекулярных движений в конденсированных средах. Сигнал ОЭК $S(\tau )$ отражает степень наведенной лазерным импульсом анизотропии среды и описывается нелинейным оптическим откликом среды третьего порядка $r(t)\,:$

(1)
$S(\tau ) = \int\limits_{ - \infty }^\infty {r(t)} G(t - \tau )dt,$
где $G(t)$ – автокорреляционная функция огибающей интенсивности лазерного импульса. В литературе описан ряд способов интерпретации сигнала ОЭК [35]. Наиболее распространенный способ состоит в извлечении из сигнала ОЭК редуцированной спектральной плотности ${{T}_{{RSD}}}(\omega ),$ которая отражает частотную зависимость комбинационной активности среды. При этом используется известная процедура деконволюции, которая позволяет разделить электронный и ядерный отклики с учетом длительности лазерных импульсов [68]. Но сначала из сигнала ОЭК требуется вычесть ориентационный вклад и получить редуцированный сигнал ОЭК ${{S}_{R}}(\tau )\,:$

(2)
${{S}_{R}}(\tau ) = S(\tau ) - \int\limits_{ - \infty }^\infty {r(t)} G(t - \tau )dt.$

Традиционно считается, что вид ориентационного отклика $r(t)$ описывается функцией [3, 4]

(3)
$\begin{gathered} r(t) = \quad\,\left( {A\exp \left( {\frac{{ - t}}{{{{\tau }^{{(1)}}}}}} \right) + \left( {1 - A} \right)\exp \left( {\frac{{ - t}}{{{{\tau }^{{(2)}}}}}} \right)} \right) \times \\ \times \,\,\left( {1 - \exp \left( {\frac{{ - t}}{\gamma }} \right)} \right), \\ \end{gathered} $
где ${{\tau }^{{(i)}}}$ – времена релаксации, которая происходит по механизму вращательной диффузии, константа A определяет относительный вклад слагаемых в выражении (2), а параметр γ симулирует нарастание ориентационного отклика, связанное с инерциальным характером молекулярных колебаний. Значения ${{\tau }^{{(i)}}}$ легко получить из анализа спада сигнала ОЭК в пикосекундном диапазоне, где сигнал ОЭК формируется только ориентационной анизотропией молекул. Для точного определения значения параметра γ требуется разделить либрационный и ориентационный вклады в диапазоне сотен фемтосекунд после возбуждающего импульса лазера, которые наряду с колебательным откликом формируют сигнал ОЭК [9]. Интуитивно время нарастания ориентационного отклика считается равным $\gamma = 200$ фс.

В работах [914] мы предложили использовать двухимпульсное возбуждение, когда создаются условия для подавления ориентационного отклика с целью реализации селективной спектроскопии комбинационно-активных межмолекулярных (либрационных) движений. Для корректного извлечения функции ${{T}_{{RSD}}}(\omega )$ из экспериментальных данных в настоящей работе мы предлагаем использовать предположение о корреляции между либрационным и ориентационным откликами. Это позволяет вычислить время нарастания ориентационного отклика ${{\gamma }^{{(i)}}}$ из соотношения:

(4)
${{\tau }_{{lib}}} = {{\gamma }^{{(i)}}}\ln \left( {{{\left( {{{\gamma }^{{(i)}}} + {{\tau }^{{(i)}}}} \right)} \mathord{\left/ {\vphantom {{\left( {{{\gamma }^{{(i)}}} + {{\tau }^{{(i)}}}} \right)} {{{\gamma }^{{(i)}}}}}} \right. \kern-0em} {{{\gamma }^{{(i)}}}}}} \right),$
где ${{\tau }_{{lib}}}$ – длительность либрационного отклика, которую легко оценить из вида сигнала ОЭК. В (4) предполагается, что каждому процессу ориентационной релаксации с временем ${{\tau }^{{(i)}}}$ соответствует свое время нарастания ${{\gamma }^{{(i)}}},$ и выражение для ориентационного отклика приобретает вид:

(5)

Ориентационный отклик в виде (5) мы использовали для вычисления редуцированного сигнала ОЭК ${{S}_{R}}(\tau ).$ Искомая информация о редуцированной спектральной плотности ${{T}_{{RSD}}}(\omega )$ извлекается из ${{S}_{R}}(\tau )$ с помощью процедуры деконволюции:

(6)
${{T}_{{RSD}}}(\omega ) = \operatorname{Im} \left( {\frac{{\int\limits_{ - \infty }^\infty {{{e}^{{i\omega t}}}} {{S}_{R}}(t)dt}}{{\int\limits_{ - \infty }^\infty {{{e}^{{i\omega t}}}} G(t)dt}}} \right).$

Таким образом, на основе описанного анализа (выражения (3)–(5)) мы впервые вычислили длительности нарастания ориентационного отклика ${{\gamma }^{{(i)}}}$ и получили обоснованный вид функции ${{T}_{{RSD}}}(\omega )$ для исследованных жидкостей толуол C7H8, нитробензол C6H5NO2 и ацетонитрил C2H3N.

ЭКСПЕРИМЕНТАЛЬНАЯ ЧАСТЬ

В качестве источника фемтосекундных импульсов использовали титан-сапфировый лазер (производство Авеста-Проект, Троицк, Россия). Средняя энергия лазерных импульсов длительностью 45 фс на полувысоте с длиной несущей волны 790 нм составила 4 × 10–9 Дж. В основе эксперимента лежала традиционная оптическая схема по регистрации сверхбыстрого ОЭК [18]. Лазерный импульс делили на пробный и накачивающий импульсы в соотношении 1 : 9. Угол между линейными поляризациями импульсов составил 45°. Далее лучи фокусировали линзой в кювету с исследуемой жидкостью, после чего с помощью диафрагмы выделяли пробный пучок. Состояние среды после воздействия импульса накачки регистрировали по изменению интенсивности пробного импульса, который проходил через жидкость, расположенную между скрещенными поляризаторами. Таким образом, величина сигнала ОЭК была пропорциональна наведенной оптической анизотропии среды. При регистрации сверхбыстрого оптического эффекта Керра использовали известные методы оптического гетеродинирования и синхронного детектирования. Все эксперименты проводили при комнатной температуре.

Регистрацию с применением синхронного детектирования осуществляли без учета фазы. В этом случае наблюдаемый сигнал ${{S}_{{det}}}(\tau )$ связан с сигналом сверхбыстрого ОЭК $S(\tau )$ соотношением:

(7)
${{S}_{{det}}}(\tau ) = \left| {S(\tau ) + {{C}_{{det}}}} \right|,$
где аппаратная постоянная ${{C}_{{det}}}$ диктуется условиями регистрации.

В качестве среды мы выбрали прозрачные в области 790 нм чистые профильтрованные (размер пор фильтра 0.2 мкм) жидкости толуол C7H8, нитробензол C6H5NO2 и ацетонитрил C2H3N, которые помещали в кварцевую кювету толщиной 2 мм.

РЕЗУЛЬТАТЫ

На рис. 1 и 2 показаны сигналы ОЭК в субпикосекундном и пикосекундном диапазонах длительностей. Из рис. 1 видно, что сигнал формируется электронным вкладом, который имеет вид автокорреляционной функции $G(\tau ),$ либрационным вкладом длительностью ${{\tau }_{{lib}}},$ вкладом внутримолекулярных колебаний и, наконец, ориентационным вкладом. Рис. 2 позволяет определить времена релаксации ориентационного отклика ${{\tau }^{{(1)}}},$ ${{\tau }^{{(2)}}}$ и величину параметра А для исследуемых жидкостей. С помощью выражения (4) можно оценить времена нарастания ориентационного отклика ${{\gamma }^{{(1)}}}$ и ${{\gamma }^{{(2)}}}.$ Таким образом, оказываются известными все необходимые параметры для симуляции ориентационного отклика на основании выражения (5), которые приведены в таблице. Это дает возможность вычислить редуцированный сигнал ОЭК ${{S}_{R}}(\tau )$ (2) и применить к нему процедуру деконволюции (6). Полученные таким образом редуцированные спектральные плотности показаны на рис. 3. Они не содержат вкладов электронного и ориентационного откликов и отражают комбинационно-активные колебания, относящиеся к межмолекулярным вращательным движениям (либрациям). Несмотря на ясное понимание того, что вид низкочастотной области спектра несет информацию о локальной организации жидкости в субпикосекундном диапазоне, его интерпретация на сегодняшний день остается нерешенной задачей. Тем не менее особенности редуцированной рамановской спектральной плотности имеют решающее значение во многих работах по симуляции межмолекулярной динамики в жидкости [1516]. Следует также добавить, что полосы в области 220 и 180 см–1 относятся к комбинационно-активным колебаниям метильной группы и нитрогруппы относительно плоскости молекул толуола и нитробензола [17]. Из рис. 1 видно, что лазерные импульсы длительностью 45 фс оказываются не способными эффективно возбуждать низкочастотные колебания молекул ацетонитрила в области 376 см–1 [13].

Рис. 1.

Субпикосекундный диапазон сигнала сверхбыстрого ОЭК в жидкостях толуол C7H8, нитробензол C6H5NO2 и ацетонитрил C2H3N при возбуждении импульсами длительностью 45 фс.

Рис. 2.

Пикосекундный диапазон сигнала сверхбыстрого ОЭК в жидкостях толуол C7H8, нитробензол C6H5NO2 и ацетонитрил C2H3N при возбуждении импульсами длительностью 45 фс. Сплошной линией показаны экспериментальные данные, пунктиром – ориентационный отклик (5).

Рис. 3.

Функции редуцированной рамановской спектральной плотности, полученные из анализа сигнала ОЭК жидкостей толуол C7H8, нитробензол C6H5NO2 и ацетонитрил C2H3N.

Таблица 1.  

Параметры молекулярной динамики, использованные для симуляции ориентационного отклика r (5)

Среды А ${{\tau }^{{(1)}}},$ пс ${{\tau }^{{(2)}}},$ пс ${{\gamma }^{{(1)}}},$ фс ${{\gamma }^{{(2)}}},$ фс ${{\tau }_{{lib}}}$, фс
Толуол, C7H8 0.25 4.5 0.6 72 240 300
Нитробензол, C6H5NO2 0.25 20 0.9 54 180 320
Ацетонитрил, C2H3N 0.6 1.4 0.5 150 520 350

ЗАКЛЮЧЕНИЕ

Представлен анализ сверхбыстрой молекулярной динамики в жидкостях толуол C7H8, нитробензол C6H5NO2 и ацетонитрил C2H3N путем регистрации сверхбыстрого ОЭК, возбуждаемого нерезонансно лазерными импульсами длительностью 45 фс. Разделение вкладов в сигнал ОЭК, соответствующих разным механизмам наведенной анизотропии среды (электронный молекулярный отклик, ядерные колебательные и вращательные движения молекул), осуществляли на основе процедуры деконволюции с учетом корреляции между либрационным и ориентационным откликами. Это позволило вычислить время нарастания ориентационного отклика и, таким образом, корректно определить редуцированный сигнал ОЭК и редуцированную рамановскую спектральную плотность ${{T}_{{RSD}}}(\omega ),$ что и являлось основной целью настоящего исследования. Следует отметить, что низкочастотная компонента функции ${{T}_{{RSD}}}(\omega )$ отражает информацию о локальной организации жидкости, в связи с чем точное определение данной спектральной компоненты носит фундаментальный характер и является важным шагом в понимании особенностей субпикосекундной межмолекулярной динамики в жидкости.

Экспериментальная часть исследования выполнена при финансовой поддержке РФФИ в рамках научного проекта № 18-32-00076, работа по моделированию откликов и теоретические расчеты выполнены при финансовой поддержке РФФИ в рамках научных проектов № 18-52-00026, № 17-02-00701 и № 19-02-00569.

Список литературы

  1. Kinoshita S., Kai Y., Ariyoshi T., Shimada Y. // Intern. J. Modern Phys. B. 1996. V. 10. P. 1229.

  2. Shirota H., Kato T. // J. Phys. Chem. A. 2011. V. 115. P. 8797.

  3. Zhong Q., Fourkas J.T. // J. Phys. Chem. B. 2008. V. 112. P. 15529.

  4. Hunt N.T., Jaye A.A., Meech S.R. // Phys. Chem. Chem. Phys. 2007. V. 9. P. 2167.

  5. Самарцев В.В., Никифоров В.Г. Фемтосекундная лазерная спектроскопия, М.: Тровант, 2017. 402 с.

  6. Mc Morrow D. // Opt. Commun. 1991. V. 86. P. 236.

  7. Mc Morrow D., Lotshaw W.T. // Chem. Phys. Lett. 1990. V. 174. P. 85.

  8. Mc Morrow D., Lotshaw W.T. // J. Phys. Chem. 1991. V. 95. P. 10395.

  9. Nikiforov V.G., Zharkov D.K., Shmelev A.G. et al. // Appl. Phys. B. 2017. V. 123. P. 209.

  10. Никифоров В.Г., Сафиуллин Г.М., Шмелев А.Г и др. // Письма в ЖЭТФ. 2007. Т. 86. С. 757.

  11. Nikiforov V.G., Shmelev A.G., Safiullin G.M. et al. // Chem. Phys. Lett. 2014. V. 592. P. 196.

  12. Nikiforov V.G., Zharkov D.K., Shmelev A.G. et al. // Appl. Phys. B. 2017. V. 123. № 7. P. 209

  13. Никифоров В.Г., Лобков В.С., Самарцев В.В. // Изв. РАН. Сер. физ. 2018. Т. 82. № 8. С. 1108. Nikiforov V.G., Lobkov V.S., Samartsev V.V. // Bull. RAS: Phys. 2018. V. 82. No 8. P. 1004.

  14. Nikiforov V.G. // Quantum Electron. 2013. V. 43. P. 177.

  15. Bender J.S., Fourkas J.T., Coasne B. // J. Phys. Chem. B. 2017. V. 121. P. 11376.

  16. Lynden-Bell R.M., Xue L., Tamas G., Quitevis E. L. // J. Chem. Phys. 2014. V. 141. Art. no. 044506.

  17. Jaye A.A., Hunt N.T., Meech S.R. // Langmuir. 2005. V. 21. P. 1238.

Дополнительные материалы отсутствуют.