Рассеяние света уверенно зарекомендовало себя как инструмент для исследования флуктуаций плотности и их динамики вблизи критической точки фазового перехода газ-жидкость. В последнее время несколько научных коллективов осуществили измерения ширины комбинационно активных переходов разных сред в ее окрестности [1, 2, 19, 17, 27, 22, 6, 7, 8, 18, 24]. В отличие от O2 [6], N2 [7, 8, 18], и HD [24]; двуокись углерода [1, 2, 19] (Q-полоса 1388 см-1) наряду с N2O [17], этаном [27] и CHF3 [22] не проявляет явным образом критического уширения. Этому факту до сих пор не удалось получить объяснения. В настоящей работе будет показано, как можно оценить критическое уширение и учесть его роль в формировании КАРС-спектра Q-полосы фермиевского дублета ν1/2ν2 двуокиси углерода.
Критическое уширение рассматриваемой Q-полосы уже было исследовано в работах Х.Л. Штрауса и Ш. Мукамела [25]. Диапазоном исследования данной работы являются области, прилегающие к критической точке. Также в данной работе предложен и использован альтернативный подход для исследования неупругого канала уширения в широком диапазоне плотностей.
Полный спектральный отклик формируется под влиянием следующих факторов: эффекта Доплера, вращательно-колебательного взаимодействия, ориентационно-колебательного взаимодействия, неупругой передачи энергии, упругой дефазировки. Было обнаружено [4] и определено [3], что вращательно-колебательное взаимодействие Q-полосы 1285 см-1 существенно отличается от такового Q-полосы 1388 см-1, что явилось неожиданным результатом. Тем не менее, при плотностях более 0.2ρc (ρc - критическая плотность), двумя первыми факторами, формирующими спектральный отклик, можно пренебречь в связи с коллапсом. [4] Ориентационно-колебательное взаимодействие растет с плотностью и является заметным, в данном случае, при плотностях выше критической. Упругая дефазировка при плотностях, близких к ρc, обусловлена преимущественно флуктуациями плотности. Особенности этого канала дефазировки вблизи критической точки хорошо известны [16].
Для расчета упругой дефазировки флюид представляется в виде N двухуровневых частиц, взаимодействующих парными аддитивными силами. Основное состояние обозначим как |a>, возбужденное - как |b >. Мы обозначим взаимодействие двух атомов в состоянии |a> с помощью Va(Q), а одного атома в состоянии |a> и другого - в состоянии |b > с помощью Vb(Q); тогда взаимодействие, ответственное за уширение спектра, можно записать в виде
.
Здесь Q - расстояние между атомами. В этом выражении сделано предположение, что состояние частицы не влияет на ее траекторию. Форма спектрального отклика данной системы может быть представлена в следующем виде [16]
T2i = 2T1 ,
где T2i = 2T1 - время дефазировки, обусловленной неупругой передачей энергии.
где Ф - случайный фазовый сдвиг дефазируемого перехода.
В пределе однородного уширения, член
в 1 уменьшается в e раз за время τ, при котором φ(τ) значительно меньше 1. Для много больших значений τ можно записать
(2)
φ(τ) будет рассчитан с помощью уравнения (9) в [16]. Для верхнего предела k использована следующая оценка
, (3)
где n - средняя численная плотность.
Уравнения состояния [23, 9, 26] позволяют рассчитать не только сжимаемость, но и термическую диффузию. Уравнение состояния [28], использованное для определения самодиффузии, дает точность 7 %. U0 (ур-е (4c) в [16]) было оценено из спектроскопических наблюдений,
,
где Ω - средняя частота спектрального профиля.
(∂Ω/∂n)2 изменяется, когда плотность пересекает околокритическую область.
Рис. 1. Длина корреляции в CO2 как функция приведенной плотности для температур 33°C (рассчитанные результаты показаны сплошной линией) и 37,1°C (1.02 T) (рассчитанные и экспериментальные результаты показаны штриховой линией и кружками, соответственно).
Рис. 2. Рассчитанные логарифмические отклики -φ Q-полос 1285 см-1 (а) и 1388 см-1 (б) фермиевского дублета ν1/2ν2 CO2 для температуры 33°C.
Графики надписаны значениями приведенной плотности и сдвинуты по вертикали для удобства. Штриховыми и пунктирными линиями показаны отклики для температур 32°C и 36°C, соответственно.
К сожалению, на данный момент, не получено точных данных для околокритических температур. Таким образом, значения (∂Ω/∂n)2 полагаются независимыми от плотности для каждой Q-полосы и более детально обсуждаются в параграфе III.
Уравнение (9) в [16] содержит длину корреляции. Значения и формулы расчета корреляционной длины для CO2 были опубликованы в нескольких статьях. Одни работы [15, 14, 12] ограничены рассмотрением исключительно случая критической изохоры, в то время как область исследования других [21, 13, 11, 20] распространяется вне изохоры. Уравнение (6.10) из [13] позволяет рассчитать длину корреляции,
,
где ξ0=1.5Å, Г=0.0481, Tref - опорная температура, подобранная таким образом, что критическое увеличение теплопроводности пренебрежимо мало для температур, больших Tref,
симметризованная изотермическая сжимаемость, γ=1.239 и ν=0.630 - критические экспоненты. В настоящей работе использовано Tref = (3/2)Tc. Примеры расчета корреляционной длины показаны на рис. 1.
Полученные результаты количественно согласуются с экспериментальными данными [20] (рис. 1, кружки) в диапазоне, где аппроксимирующая функция статического структурного фактора Фишера-Бурфорда (уравнение (8) в [16], а также [10, 5]) чувствительна к длине корреляции (т. е. где последняя больше среднего межатомного расстояния ≈5 Å).
Рассчитанные логарифмические отклики -φ Q-полос 1285 см-1 и 1388 см-1 фермиевского дублета ν1/2ν2 для температуры 33○C показаны на рис. 2.
Как ожидалось [16], отклики сильнее проявляют неоднородные свойства вблизи критической точки. Отклики при температуре 32°C не демонстрируют весомых различий на временах, при которых уровень отклика больше, чем 1/e от начального. Данный факт свидетельствует о том, что критическое уширение не чувствительно к малым отстройкам температур относительно критической. То же самое справедливо, когда смещение относительно критической точки берется по плотности. Таким образом, температуру 33°C можно считать «эквивалентной критической» с точки зрения спектроскопии. Существует еще одна причина выбора данной температуры. Она состоит в более высокой точности эксперимента на шкале плотности, чем для температур, более близких к критической.
В настоящей работе рассчитаны отклики Q-полос 1285 см-1 и 1388 см-1 фермиевского дублета ν1/2ν2 двуокиси углерода. Новые данные термодинамических свойств и корреляционной длины, доступные на сегодняшний день, позволяют расширить область исследования по плотности и температуре вдаль от критической точки. Показано, что критическое уширение не чувствительно к температуре и плотности вблизи критической точки.
Список литературы