Профиль Voigt

редактировать
(Центрированный) Voigt
Функция плотности вероятности График центрированного профиля Фойгта для четырех случаев . Построение центрированного профиля Voigt для четырех случаев. Каждый корпус имеет полную ширину на полувысоте почти 3,6. Черный и красный профили - это предельные случаи гауссова (γ = 0) и лоренцевского (σ = 0) профилей соответственно.
Кумулятивная функция распределения По центру Фойгт CDF.
Параметрыγ, σ>0 {\ displaystyle \ gamma, \ sigma>0}\gamma,\sigma>0
Поддержка x ∈ (- ∞, ∞) {\ displaystyle x \ in (- \ infty, \ infty)}x \ in (- \ infty, \ infty)
PDF ℜ [w (z)] σ 2 π, z = x + i γ σ 2 {\ displaystyle {\ frac {\ Re [w (z)]} { \ sigma {\ sqrt {2 \ pi}}}}, ~~~ z = {\ frac {x + i \ gamma} {\ sigma {\ sqrt {2}}}}}\ frac {\ Re [w (z)]} {\ sigma \ sqrt {2 \ pi}}, ~~~ z = \ frac {x + i \ gamma} {\ sigma \ sqrt {2}}
CDF ( сложный - см. текст)
Среднее (не определено)
Медиана 0 {\ displaystyle 0}{\ displaystyle 0}
Режим 0 {\ displaystyle 0}{\ displaystyle 0}
Дисперсия ( не определено)
Асимметрия (не определено)
Пример. эксцесс (не определено)
MGF (не определено)
CF e - γ | t | - σ 2 t 2/2 {\ displaystyle e ^ {- \ gamma | t | - \ sigma ^ {2} t ^ {2} / 2}}e ^ {- \ gamma | t | - \ sigma ^ 2 t ^ 2/2}

Профиль Voigt (назван в честь Woldemar Voigt ) - это распределение вероятностей, полученное посредством свертки распределения Коши-Лоренца. ion и распределение Гаусса. Он часто используется при анализе данных спектроскопии или дифракции.

Содержание
  • 1 Определение
  • 2 История и приложения
  • 3 Свойства
    • 3.1 Кумулятивная функция распределения
    • 3.2 Нецентрированный профиль Фойгта
    • 3.3 Производный профиль
  • 4 Функции Фойгта
    • 4.1 Связь с профилем Фойгта
  • 5 Числовые приближения
    • 5.1 Аппроксимация псевдо-Фойгта
    • 5.2 Ширина профиля Фойгта
  • 6 Ссылки
  • 7 Внешние ссылки
Определение

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

V (x; σ, γ) ≡ ∫ - ∞ ∞ G (x ′; σ) L (x - x ′; γ) dx ′, {\ displaystyle V (x; \ sigma, \ gamma) \ Equiv \ int _ {- \ infty} ^ {\ infty} G (x '; \ sigma) L (x-x'; \ gamma) \, dx ',}{\displaystyle V(x;\sigma,\gamma)\equiv \int _{-\infty }^{\infty }G(x';\sigma)L(x-x';\gamma)\,dx',}

где x - сдвиг от центра линии G (x; σ) {\ displaystyle G (x; \ sigma)}G (x; \ sigma) - это центрированный гауссовский профиль:

G (x; σ) ≡ e - x 2 / (2 σ 2) σ 2 π, {\ Displaystyle G (x; \ sigma) \ Equiv {\ frac {e ^ {- x ^ {2} / (2 \ sigma ^ {2})}} {\ sigma {\ sqrt {2 \ pi}}}},}{\ displaystyle G (x; \ sigma) \ Equiv {\ frac {e ^ { -x ^ {2} / (2 \ sigma ^ {2})}} {\ sigma {\ sqrt {2 \ pi}}}},}

и L (x; γ) {\ displaystyle L (x; \ gamma)}L (x; \ gamma) - это центрированный лоренцевский профиль:

L (x; γ) ≡ γ π (x 2 + γ 2). {\ displaystyle L (x; \ gamma) \ Equiv {\ frac {\ gamma} {\ pi (x ^ {2} + \ gamma ^ {2})}}.}{\ displaystyle L (x; \ gamma) \ Equiv {\ frac {\ gamma} {\ pi ( x ^ {2} + \ gamma ^ {2})}}.}

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

В (Икс; σ, γ) знак равно Re ⁡ [вес (z)] σ 2 π, {\ Displaystyle V (x; \ sigma, \ gamma) = {\ frac {\ operatorname {Re} [w (z)]} {\ sigma {\ sqrt {2 \ pi}}}},}{\ displaystyle V (x; \ sigma, \ gamma) = {\ frac {\ operatorname {Re} [w (z)]} {\ sigma {\ sqrt {2 \ pi}}}},}

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

г = х + я γ σ 2. {\ displaystyle z = {\ frac {x + i \ gamma} {\ sigma {\ sqrt {2}}}}.}z = \ frac {x + i \ gamma} {\ sigma \ sqrt {2}}.

В предельных случаях σ = 0 {\ displaystyle \ sigma = 0 }\ sigma = 0 и γ = 0 {\ displaystyle \ gamma = 0}{\ displaystyle \ gamma = 0} , затем V (x; σ, γ) {\ displaystyle V (x; \ sigma, \ gamma)}{\ displaystyle V (x; \ sigma, \ gamma)} упрощается до L (x; γ) {\ displaystyle L (x; \ gamma)}L (x; \ gamma) и G (x; σ) {\ displaystyle G (x; \ sigma)}G (x; \ sigma) соответственно.

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

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

Свойства

Профиль Фойгта нормализован:

∫ - ∞ ∞ V (x; σ, γ) dx = 1, {\ displaystyle \ int _ {- \ infty} ^ {\ infty} V (x; \ sigma, \ gamma) \, dx = 1,}{\ displaystyle \ int _ {- \ infty} ^ {\ infty} V (x; \ sigma, \ gamma) \, dx = 1,}

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

φ f (t; σ, γ) = E (eixt) = e - σ 2 t 2 / 2 - γ | т |. {\ displaystyle \ varphi _ {f} (t; \ sigma, \ gamma) = E (e ^ {ixt}) = e ^ {- \ sigma ^ {2} t ^ {2} / 2- \ gamma | t |}.}\ varphi_f (t; \ sigma, \ gamma) = E (e ^ {ixt}) = e ^ {- \ sigma ^ 2t ^ 2/2 - \ gamma | t |}.

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

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

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

F (x 0; μ, σ) = ∫ - ∞ x 0 Re ⁡ (w (z)) σ 2 π dx = Re ⁡ (1 π ∫ z (- ∞) z (x 0) w (z) dz). {\ Displaystyle F (x_ {0}; \ mu, \ sigma) = \ int _ {- \ infty} ^ {x_ {0}} {\ frac {\ operatorname {Re} (w (z))} {\ sigma {\ sqrt {2 \ pi}}}} \, dx = \ operatorname {Re} \ left ({\ frac {1} {\ sqrt {\ pi}}} \ int _ {z (- \ infty)} ^ {z (x_ {0})} w (z) \, dz \ right).}{\ displaystyle F (x_ {0}; \ mu, \ sigma) = \ int _ {- \ infty} ^ {x_ { 0}} {\ frac {\ operatorname {Re} (w (z))} {\ sigma {\ sqrt {2 \ pi}}}} \, dx = \ operatorname {Re} \ left ({\ frac {1 } {\ sqrt {\ pi}}} \ int _ {z (- \ infty)} ^ {z (x_ {0})} w (z) \, dz \ right).}

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

1 π ∫ w (z) dz = 1 π ∫ e - z 2 [1 - erf ⁡ (- iz)] dz, {\ displaystyle {\ frac {1} {\ sqrt { \ pi}}} \ int w (z) \, dz = {\ frac {1} {\ sqrt {\ pi}}} \ int e ^ {- z ^ {2}} \ left [1- \ operatorname { erf} (-iz) \ right] \, dz,}{\ displaystyle {\ frac {1} {\ sqrt {\ pi}}} \ int w ( z) \, dz = {\ frac {1} {\ sqrt {\ pi}}} \ int e ^ {- z ^ {2}} \ left [1- \ operatorname {erf} (-iz) \ right] \, dz,}

которое может быть решено, чтобы дать

1 π ∫ w (z) dz = erf ⁡ (z) 2 + iz 2 π 2 F 2 (1, 1; 3 2, 2; - z 2), {\ displaystyle {\ frac {1} {\ sqrt {\ pi}}} \ int w (z) \, dz = {\ frac {\ operatorname {erf} (z)} {2}} + {\ frac {iz ^ {2}} {\ pi}} \, _ {2} F_ {2} \ left (1,1; {\ frac {3} {2} }, 2; -z ^ {2} \ right),}{\ displaystyle {\ frac {1} {\ sqrt {\ pi}}} \ int w (z) \, dz = {\ frac {\ operatorname {erf} (z)} {2}} + {\ frac {iz ^ {2}} {\ pi}} \, _ {2} F_ {2} \ left (1,1; {\ frac {3} {2}}, 2; -z ^ {2} \ right),}

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

F (x; μ, σ) = Re ⁡ [1 2 + erf ⁡ (z) 2 + iz 2 π 2 F 2 (1, 1; 3 2, 2; - z 2)]. {\ displaystyle F (x; \ mu, \ sigma) = \ operatorname {Re} \ left [{\ frac {1} {2}} + {\ frac {\ operatorname {erf} (z)} {2}} + {\ frac {iz ^ {2}} {\ pi}} \, _ {2} F_ {2} \ left (1,1; {\ frac {3} {2}}, 2; -z ^ { 2} \ right) \ right].}{\ displaystyle F (x; \ mu, \ sigma) = \ operatorname {Re} \ left [{\ frac {1} {2}} + {\ frac {\ operatorname {erf} (z)} {2}} + {\ frac {iz ^ {2}} {\ pi}} \, _ {2} F_ {2} \ left (1,1; {\ frac {3} {2}}, 2; -z ^ {2} \ right) \ right].}

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

Если гауссовский профиль центрирован на μ G {\ displaystyle \ mu _ {G}}\ mu_G и лоренцевский профиль центрируется на μ L {\ displaystyle \ mu _ {L}}\ mu_L , центр свертки - на μ G + μ L {\ displaystyle \ mu _ { G} + \ mu _ {L}}\mu_G+\mu_Lи характеристическая функция равна

φ f (t; σ, γ, μ G, μ L) = ei (μ G + μ L) t - σ 2 t 2/2 - γ | т |. {\ Displaystyle \ varphi _ {е} (т; \ сигма, \ гамма, \ му _ {\ mathrm {G}}, \ му _ {\ mathrm {L}}) = е ^ {я (\ му _ { \ mathrm {G}} + \ mu _ {\ mathrm {L}}) t- \ sigma ^ {2} t ^ {2} / 2- \ gamma | t |}.}\ varphi_f (t; \ sigma, \ gamma, \ mu_ \ mathrm {G}, \ mu_ \ mathrm {L}) = e ^ {i (\ mu_ \ mathrm {G} + \ mu_ \ mathrm {L}) t- \ sigma ^ 2t ^ 2/2 - \ gamma | t |}.

Мода и медиана оба расположены в μ G + μ L {\ displaystyle \ mu _ {G} + \ mu _ {L}}\mu_G+\mu_L.

Профиль производной

Профили первой и второй производной могут быть выражены в терминах функции Фаддеева следующим образом:

∂ ∂ x V (x; σ, γ) = - Re ⁡ [zw (z)] σ 2 π = - x σ 2 Re ⁡ [w ( z)] σ 2 π + γ σ 2 Im ⁡ [w (z)] σ 2 π; {\ displaystyle {\ frac {\ partial} {\ partial x}} V (x; \ sigma, \ gamma) = - {\ frac {\ operatorname {Re} [z ~ w (z)]} {\ sigma ^ {2} {\ sqrt {\ pi}}}} = - {\ frac {x} {\ sigma ^ {2}}} {\ frac {\ operatorname {Re} [w (z)]} {\ sigma { \ sqrt {2 \ pi}}}} + {\ frac {\ gamma} {\ sigma ^ {2}}} {\ frac {\ operatorname {Im} [w (z)]} {\ sigma {\ sqrt { 2 \ pi}}}};}{\ displaystyle {\ frac {\ partial} {\ partial x}} V (x; \ sigma, \ gamma) = - {\ frac {\ operatorname { Re} [z ~ w (z)]} {\ sigma ^ {2} {\ sqrt {\ pi}}}} = - {\ frac {x} {\ sigma ^ {2}}} {\ frac {\ OperatorName {Re} [w (z)]} {\ sigma {\ sqrt {2 \ pi}}}} + {\ frac {\ gamma} {\ sigma ^ {2}}} {\ frac {\ operatorname {Im } [w (z)]} {\ sigma {\ sqrt {2 \ pi}}}};}
∂ 2 ∂ x 2 V (x; σ, γ) = x 2 - γ 2 - σ 2 σ 4 Re ⁡ [w (z)] σ 2 π - 2 x γ σ 4 Им ⁡ [вес (Z)] σ 2 π + γ σ 4 1 π, {\ Displaystyle {\ frac {\ partial ^ {2}} {\ partial x ^ {2}}} V (x; \ sigma, \ gamma) = {\ frac {x ^ {2} - \ gamma ^ {2} - \ sigma ^ {2}} {\ sigma ^ {4}}} {\ frac {\ operatorname {Re} [w (z)]} {\ sigma {\ sqrt {2 \ pi}}}} - {\ frac {2x \ gamma} {\ sigma ^ {4}}} {\ frac {\ operatorname {Im} [w (z)]} {\ sigma {\ sqrt {2 \ pi}}}} + {\ frac {\ gamma} {\ sigma ^ {4}}} {\ frac {1} {\ pi}},}{\ displaystyle {\ frac {\ partial ^ {2}} {\ partial x ^ {2}}} V (x; \ sigma, \ gamma) = {\ frac {x ^ {2} - \ gamma ^ {2} - \ sigma ^ {2}} {\ sigma ^ {4}} } {\ frac {\ operatorname {Re} [w (z)]} {\ sigma {\ sqrt {2 \ pi}}}} - {\ frac {2x \ gamma} {\ sigma ^ {4}}} { \ frac {\ operatorname {Im} [w (z)]} {\ sigma {\ sqrt {2 \ pi}}}} + {\ frac {\ gamma} {\ sigma ^ {4}}} {\ frac { 1} {\ pi}},}

используя приведенное выше определение для z.

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

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

U (Икс, T) + я В (Икс, T) знак равно π 4 tez 2 erfc ⁡ (z) = π 4 tw (iz), {\ Displaystyle U (x, t) + iV (x, t) = {\ sqrt {\ frac {\ pi} {4t}}} e ^ {z ^ {2}} \ operatorname {erfc} (z) = {\ sqrt {\ frac {\ pi} {4t}}} w (iz),}{\ displaystyle U (x, t) + iV (x, t) = {\ sqrt {\ frac {\ pi} {4t}}} e ^ {z ^ {2}} \ operatorname {erfc} (z) = {\ sqrt {\ frac {\ pi} {4t}}} w (iz),}
ЧАС (a, u) = U (u / a, 1/4 a 2) a π, {\ displaystyle H (a, u) = {\ frac {U (u / a, 1 / 4a ^ {2})} {a {\ sqrt {\ pi}}}},}{\ displaystyle H (a, u) = {\ frac {U (u / a, 1 / 4a ^ {2})} {a {\ sqrt {\ pi}}}},}

где

z = (1 - ix) / 2 t, {\ displaystyle z = (1-ix) / 2 {\ sqrt {t}},}z = (1-ix) / 2 \ sqrt t,

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

Связь с профилем Voigt

V (x ; σ, γ) знак равно ЧАС (a, u) / (2 π σ), {\ Displaystyle V (x; \ sigma, \ gamma) = H (a, u) / ({\ sqrt {2}} {\ sqrt {\ pi}} \ sigma),}{\ displaystyle V (x; \ sigma, \ gamma) = H (а, и) / ({\ sqrt {2}} {\ sqrt {\ pi}} \ sigma),}

с

a = γ / (2 σ) {\ displaystyle a = \ gamma / ({\ sqrt {2}} \ sigma)}{\ displaystyle a = \ gamma / ({\ sqrt {2} } \ sigma)}

и

и = х / (2 σ). {\ displaystyle u = x / ({\ sqrt {2}} \ sigma).}{\ displaystyle u = х / ({\ sqrt {2}} \ sigma).}
Числовые приближения

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

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

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

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

V p (x, f) = η ⋅ L (x, f) + (1 - η) ⋅ G (x, f) {\ displaystyle V_ {p} (x, е) = \ эта \ cdot L (x, f) + (1- \ eta) \ cdot G (x, f)}{\ displaystyle V_ {p} (x, f) = \ eta \ cdot L (x, f) + (1- \ eta) \ cdot G (x, f)} с 0 < η < 1 {\displaystyle 0<\eta <1}0 <\ eta <1 .

η {\ displaystyle \ eta}\ eta - это функция полной ширины на половине максимального значения параметра (FWHM).

Существует несколько возможных вариантов выбора параметра η {\ displaystyle \ eta}\ eta . Простая формула с точностью до 1%:

η = 1,36603 (f L / f) - 0,47719 (f L / f) 2 + 0,11116 (f L / f) 3, {\ displaystyle \ eta = 1,36603 (f_ {L} / f) -0,47719 (f_ {L} / f) ^ {2} +0,11116 (f_ {L} / f) ^ {3},}{\ displaystyle \ eta = 1,36603 (f_ {L } / f) -0,47719 (f_ {L} / f) ^ {2} +0,11116 (f_ {L} / f) ^ {3},}

где теперь, η {\ displaystyle \ eta}\ eta является функцией Лоренца (f L {\ displaystyle f_ {L}}{\ displaystyle f_ {L}} ), Gaussian (f G {\ displaystyle f_ {G}}{\ displaystyle f_ {G}} ) и всего (f {\ displaystyle f}f ) Параметры полной ширины на половине максимума (FWHM). Общая ширина FWHM (f {\ displaystyle f}f ) параметр описывается следующим образом:

f = [f G 5 + 2.69269 f G 4 f L + 2.42843 f G 3 f L 2 + 4.47163 f G 2 f L 3 + 0,07842 f G f L 4 + f L 5] 1/5. {\ displaystyle f = [f_ {G} ^ {5} + 2.69269f_ {G} ^ {4} f_ {L} + 2.42843f_ {G} ^ {3} f_ {L} ^ {2} + 4,47163f_ {G} ^ {2} f_ {L} ^ {3} + 0,07842f_ {G} f_ {L} ^ {4} + f_ {L} ^ {5}] ^ {1 / 5}.}{\ displaystyle f = [f_ {G} ^ {5} + 2.69269f_ {G} ^ {4} f_ {L} + 2.42843f_ {G} ^ {3} f_ {L} ^ {2} + 4,47163f_ {G} ^ {2} f_ {L} ^ {3} + 0,07842f_ {G} f_ {L} ^ {4} + f_ {L} ^ {5}] ^ {1/5}.}

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

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

f G = 2 σ 2 ln ⁡ (2). {\ displaystyle f _ {\ mathrm {G}} = 2 \ sigma {\ sqrt {2 \ ln (2)}}.}{\ displaystyle f _ {\ mathrm {G} } = 2 \ sigma {\ sqrt {2 \ ln (2)}}.}

FWHM лоренцевского профиля составляет

f L = 2 γ. {\ displaystyle f _ {\ mathrm {L}} = 2 \ gamma.}{\ displaystyle f _ {\ mathrm {L}} = 2 \ gamma.}

Грубое приближение для соотношения между ширинами профилей Фойгта, Гаусса и Лоренца:

f V ≈ f L / 2 + f L 2/4 + f G 2. {\ displaystyle f _ {\ mathrm {V}} \ приблизительно f _ {\ mathrm {L}} / 2 + {\ sqrt {f _ {\ mathrm {L}} ^ {2} / 4 + f _ {\ mathrm {G} } ^ {2}}}.}{\ displaystyle f _ {\ mathrm {V}} \ приблизительно f _ {\ mathrm {L}} / 2 + {\ sqrt {f _ {\ mathrm {L}} ^ {2 } / 4 + f _ {\ mathrm {G}} ^ {2}}}.}

Это приближение в точности верно для чистого гауссовского.

Лучшее приближение с точностью 0,02% дает

f V ≈ 0,5346 f L + 0,2166 f L 2 + f G 2. {\ displaystyle f _ {\ mathrm {V}} \ приблизительно 0,5346f _ {\ mathrm {L}} + {\ sqrt {0,2166f _ {\ mathrm {L}} ^ {2} + f _ {\ mathrm {G}} ^ {2}}}.}f_ \ mathrm {V} \ приблизительно 0,5346 f_ \ mathrm {L} + \ sqrt {0,2166f_ \ mathrm {L} ^ 2 + f_ \ mathrm {G} ^ 2}.

Это приближение точно верно для чистого гауссовского профиля, но имеет ошибку около 0,000305% для чистого лоренцевского профиля.

Ссылки
Внешние ссылки
  • http://jugit.fz-juelich.de/mlz/libcerf, числовая библиотека C для сложных функций ошибок, предоставляет функцию voigt ( x, sigma, gamma) с точностью приблизительно 13–14 знаков.
  • Оригинальная статья: Voigt, Woldemar, 1912, «Das Gesetz der Intensitätsverteilung innerhalb der Linien eines Gasspektrums», Sitzungsbericht der Bayerisr der Bayerischen Akademie Wissenschaften, 25, 603 (см. Также: http://publikationen.badw.de/de/003395768)
Последняя правка сделана 2021-06-18 04:57:00
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте