stringtranslate.com

Профиль Фойгта

Профиль Фойгта (названный в честь Вольдемара Фойгта ) — это распределение вероятностей , заданное сверткой распределения Коши -Лоренца и распределения Гаусса . Его часто используют при анализе данных спектроскопии или дифракции .

Определение

Без ограничения общности можно рассматривать только центрированные профили с нулевым максимумом. Тогда профиль Фойгта

где x — смещение от центра линии, — центрированный гауссов профиль:

и представляет собой центрированный лоренцев профиль:

Определяющий интеграл можно оценить как:

где Re[ w ( z )] — действительная часть функции Фаддеевой , рассчитанная для

В предельных случаях и затем упрощается до и соответственно.

История и приложения

В спектроскопии профиль Фойгта является результатом свертки двух механизмов уширения, один из которых сам по себе создает гауссов профиль (обычно в результате доплеровского уширения ), а другой - лоренцев. Профили Фойгта распространены во многих разделах спектроскопии и дифракции . Из-за затрат на вычисление функции Фаддеевой профиль Фойгта иногда аппроксимируют с помощью профиля псевдо-Фойгта.

Характеристики

Профиль Фойгта нормализуется:

поскольку это свертка нормализованных профилей. Лоренцев профиль не имеет моментов (кроме нулевого), поэтому производящая функция для распределения Коши не определена. Отсюда следует, что профиль Фойгта также не будет иметь производящей момент функции, но характеристическая функция для распределения Коши четко определена, как и характеристическая функция для нормального распределения . Тогда характеристическая функция для (центрированного) профиля Фойгта будет произведением двух:

Поскольку нормальное распределение и распределение Коши являются стабильными распределениями , каждое из них замкнуто при свертке (с точностью до изменения масштаба), а из этого следует, что распределения Фойгта также замкнуты при свертке.

Кумулятивная функция распределения

Используя приведенное выше определение z , кумулятивную функцию распределения (CDF) можно найти следующим образом:

Подстановка определения функции Фаддеевой (масштабированная комплексная функция ошибок ) дает неопределенный интеграл:

которое можно решить, чтобы получить

где – гипергеометрическая функция . Чтобы функция приближалась к нулю по мере того, как x приближается к отрицательной бесконечности (как это должно делать CDF), необходимо добавить константу интегрирования 1/2. Это дает для CDF Фойгта:

Нецентрированный профиль Фойгта

Если гауссов профиль центрирован в, а лоренцев профиль центрирован в , свертка центрируется в и характеристическая функция:

Функция плотности вероятности просто смещается от центрированного профиля на :

где:

И мода, и медиана расположены в точке .

Производные

Профиль Фойгта (здесь при условии , и ) и его первые две частные производные по (первому столбцу) и трем параметрам , , и (второй, третий и четвертый столбцы соответственно), полученные аналитически и численно.

Используя приведенное выше определение для и , первая и вторая производные могут быть выражены через функцию Фаддеевой как

и

соответственно.

Часто один или несколько профилей Фойгта и/или их соответствующие производные необходимо согласовать с измеряемым сигналом с помощью нелинейного метода наименьших квадратов , например, в спектроскопии . Затем для ускорения вычислений можно использовать дополнительные частные производные. Вместо аппроксимации матрицы Якоби по параметрам , , и с помощью конечных разностей можно применить соответствующие аналитические выражения. С и они определяются как:

для исходного профиля Фойгта ;

для частной производной первого порядка ; и

для частной производной второго порядка . Поскольку и играют относительно схожую роль в вычислении , их соответствующие частные производные также выглядят весьма схожими с точки зрения их структуры, хотя они приводят к совершенно разным профилям производных. Действительно, частные производные по и демонстрируют большее сходство, поскольку оба являются параметрами ширины. Все эти производные включают в себя только простые операции (умножения и сложения), поскольку они требуют больших вычислительных затрат и их легко получить при вычислении . Такое повторное использование предыдущих расчетов позволяет получить вывод с минимальными затратами. Это не относится к аппроксимации градиента конечной разностью, поскольку оно требует оценки для каждого градиента соответственно.

Функции Фойгта

Функции Фойгта [1] U , V и H (иногда называемые функцией уширения линии ) определяются формулами

где

erfc — дополнительная функция ошибок , а w ( z ) — функция Фаддеевой .

Связь с профилем Фойгта

с относительными переменными гауссовой сигмы и

Числовые аппроксимации

Функция Теппера-Гарсиа

Функция Теппера-Гарсиа , названная в честь немецко-мексиканского астрофизика Тора Теппера-Гарсии, представляет собой комбинацию экспоненциальной функции и рациональных функций, которая аппроксимирует функцию уширения линии в широком диапазоне ее параметров. [2] Он получается путем разложения в усеченный степенной ряд точной функции уширения линии.

В наиболее эффективной в вычислительном отношении форме функция Теппера-Гарсиа может быть выражена как

где , , и .

Таким образом, функцию уширения линии можно рассматривать в первом порядке как чистую функцию Гаусса плюс поправочный коэффициент, который линейно зависит от микроскопических свойств поглощающей среды (закодирован в ); однако из-за раннего усечения разложения в ряд ошибка аппроксимации все еще имеет порядок , т.е. Это приближение имеет относительную точность

во всем диапазоне длин волн при условии, что . Помимо высокой точности, эта функция проста в реализации и требует быстрых вычислений. Он широко используется в области анализа линий поглощения квазаров. [3]

Псевдо-фойгтовское приближение

Профиль псевдо-Фойгта (или функция псевдо-Фойгта ) — это аппроксимация профиля Фойгта V ( x ) с использованием линейной комбинации гауссовой кривой G ( x ) и лоренцевой кривой L ( x ) вместо их свертки .

Функция псевдо-Фойгта часто используется для расчета экспериментальных форм спектральных линий .

Математическое определение нормализованного профиля псевдо-Фойгта дается формулой

с .

является функцией параметра полной ширины на половине максимума (FWHM).

Существует несколько возможных вариантов выбора параметра. [4] [5] [6] [7] Простая формула с точностью до 1%: [8] [9]

где теперь является функцией параметров Лоренца ( ), Гаусса ( ) и полной ( ) Полная ширина на половине высоты (FWHM). Общий параметр FWHM ( ) описывается следующим образом:

Ширина профиля Фойгта

Полную ширину на половине высоты (FWHM) профиля Фойгта можно найти по ширине соответствующих гауссовских и лоренцевых ширин. Полуширина гауссова профиля равна

Полуширина лоренцева профиля равна

Примерное соотношение (с точностью до 1,2%) между шириной профилей Фойгта, Гаусса и Лоренца: [10]

По построению это выражение является точным для чистого гауссиана или лоренциана.

Лучшее приближение с точностью 0,02% дает [11] (первоначально найдено Килкопфом [12] ).

Опять же, это выражение точно для чистого гауссова или лоренциана. В той же публикации [11] можно найти несколько более точное (в пределах 0,012%), но значительно более сложное выражение.

Рекомендации

  1. ^ Темме, Нью-Мексико (2010), «Функция Фойгта», в Олвере, Фрэнке В.Дж .; Лозье, Дэниел М.; Буасверт, Рональд Ф.; Кларк, Чарльз В. (ред.), Справочник NIST по математическим функциям , издательство Кембриджского университета, ISBN 978-0-521-19225-5, МР  2723248.
  2. ^ Теппер-Гарсия, Торстен (2006). «Подбор профиля Фойгта к линиям поглощения квазаров: аналитическое приближение к функции Фойгта-Хьертинга». Ежемесячные уведомления Королевского астрономического общества . 369 (4): 2025–2035. arXiv : astro-ph/0602124 . Бибкод : 2006MNRAS.369.2025T. дои : 10.1111/j.1365-2966.2006.10450.x. S2CID  16981310.
  3. ^ Список цитат, найденных в Системе астрофизических данных SAO/NASA (ADS): https://ui.adsabs.harvard.edu/abs/2006MNRAS.369.2025T/citations.
  4. ^ Вертхайм Г.К., Батлер М.А., Вест К.В., Бьюкенен Д.Н. (1974). «Определение гауссовского и лоренцевого содержания экспериментальных форм линий». Обзор научных инструментов . 45 (11): 1369–1371. Бибкод : 1974RScI...45.1369W. дои : 10.1063/1.1686503.
  5. ^ Санчес-Бахо, Ф.; Флорида Камбрера (август 1997 г.). «Использование функции псевдо-Фойгта в дисперсионном методе анализа уширения рентгеновских линий». Журнал прикладной кристаллографии . 30 (4): 427–430. дои : 10.1107/S0021889896015464.
  6. ^ Лю Ю, Линь Дж, Хуан Г, Го Ю, Дуань С (2001). «Простое эмпирически-аналитическое приближение профиля Фойгта». ЖОСА Б. 18 (5): 666–672. Бибкод : 2001JOSAB..18..666L. дои : 10.1364/josab.18.000666.
  7. ^ Ди Рокко Х.О., Крузадо А (2012). «Профиль Фойгта как сумма функций Гаусса и Лоренца, когда весовой коэффициент зависит только от соотношения ширин». Acta Physica Polonica А. 122 (4): 666–669. Бибкод : 2012AcPPA.122..666D. дои : 10.12693/APhysPolA.122.666 . ISSN  0587-4246.
  8. ^ Ида Т, Андо М, Торая Х (2000). «Расширенная функция псевдо-Фойгта для аппроксимации профиля Фойгта». Журнал прикладной кристаллографии . 33 (6): 1311–1316. дои : 10.1107/s0021889800010219. S2CID  55372305.
  9. ^ П. Томпсон, Д.Э. Кокс и Дж.Б. Гастингс (1987). «Уточнение Ритвельдом данных синхротронного рентгеновского излучения Дебая-Шеррера от Al 2 O 3 ». Журнал прикладной кристаллографии . 20 (2): 79–83. дои : 10.1107/S0021889887087090.
  10. ^ Уайтинг, Э.Э. (июнь 1968 г.). «Эмпирическое приближение к профилю Фойгта». Журнал количественной спектроскопии и переноса излучения . 8 (6): 1379–1384. Бибкод : 1968JQSRT...8.1379W. дои : 10.1016/0022-4073(68)90081-2. ISSN  0022-4073.
  11. ^ Аб Оливеро, Джей Джей; Р. Л. Лонгботум (февраль 1977 г.). «Эмпирическое соответствие ширине линии Фойгта: краткий обзор». Журнал количественной спектроскопии и переноса излучения . 17 (2): 233–236. Бибкод : 1977JQSRT..17..233O. дои : 10.1016/0022-4073(77)90161-3. ISSN  0022-4073.
  12. ^ Джон Ф. Килкопф (1973), «Новое приближение функции Фойгта с применением к анализу профиля спектральных линий», Журнал Оптического общества Америки , 63 (8): 987, Бибкод : 1973JOSA...63.. 987К, номер документа :10.1364/JOSA.63.000987

Внешние ссылки