Puneet Varma (Editor)

Voigt profile

Updated on
Edit
Like
Comment
Share on FacebookTweet on TwitterShare on LinkedInShare on Reddit
Mean
  
(not defined)

Voigt profile

Parameters
  
γ , σ > 0 {\displaystyle \gamma ,\sigma >0}

Support
  
x ∈ ( − ∞ , ∞ ) {\displaystyle 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}}}}}

CDF
  
(complicated - see text)

Median
  
0 {\displaystyle 0}

In spectroscopy, the Voigt profile (named after Woldemar Voigt) is a line profile resulting from the convolution of two broadening mechanisms, one of which alone would produce a Gaussian profile (usually, as a result of the Doppler broadening), and the other would produce a Lorentzian profile. Voigt profiles are common in many branches of spectroscopy and diffraction. Due to the computational expense of the convolution operation, the Voigt profile is often approximated using a pseudo-Voigt profile.

Contents

All normalized line profiles can be considered to be probability distributions. The Gaussian profile has a Gaussian, or normal, distribution and a Lorentzian profile has a Lorentz, or Cauchy, distribution. Without loss of generality, we can consider only centered profiles, which peak at zero. The Voigt profile is then a convolution of a Lorentz profile and a Gaussian profile:

V ( x ; σ , γ ) = G ( x ; σ ) L ( x x ; γ ) d x ,

where x is the shift from the line center, G ( x ; σ ) is the centered Gaussian profile:

G ( x ; σ ) e x 2 / ( 2 σ 2 ) σ 2 π ,

and L ( x ; γ ) is the centered Lorentzian profile:

L ( x ; γ ) γ π ( x 2 + γ 2 ) .

The defining integral can be evaluated as:

V ( x ; σ , γ ) = Re [ w ( z ) ] σ 2 π ,

where Re[w(z)] is the real part of the Faddeeva function evaluated for

z = x + i γ σ 2 .

Properties

The Voigt profile is normalized:

V ( x ; σ , γ ) d x = 1 ,

since it is a convolution of normalized profiles. The Lorentzian profile has no moments (other than the zeroth), and so the moment-generating function for the Cauchy distribution is not defined. It follows that the Voigt profile will not have a moment-generating function either, but the characteristic function for the Cauchy distribution is well defined, as is the characteristic function for the normal distribution. The characteristic function for the (centered) Voigt profile will then be the product of the two:

φ f ( t ; σ , γ ) = E ( e i x t ) = e σ 2 t 2 / 2 γ | t | .

Since normal distributions and Cauchy distributions are stable distributions, they are each are closed under convolution (up to change of scale), and it follows that the Voigt distributions are also closed under convolution.

Cumulative distribution function

Using the above definition for z , the cumulative distribution function (CDF) can be found as follows:

F ( x 0 ; μ , σ ) = x 0 R e ( w ( z ) ) σ 2 π d x = R e ( 1 π z ( ) z ( x 0 ) w ( z ) d z ) .

Substituting the definition of the Faddeeva function (scaled complex error function) yields for the indefinite integral:

1 π w ( z ) d z = 1 π e z 2 [ 1 e r f ( i z ) ] d z ,

which may be solved to yield

1 π w ( z ) d z = e r f ( z ) 2 + i z 2 π 2 F 2 ( 1 , 1 ; 3 2 , 2 ; z 2 ) ,

where 2 F 2 ( ) is a hypergeometric function. In order for the function to approach zero as x approaches negative infinity (as the CDF must do), an integration constant of 1/2 must be added. This gives for the CDF of Voigt:

F ( x ; μ , σ ) = R e [ 1 2 + e r f ( z ) 2 + i z 2 π 2 F 2 ( 1 , 1 ; 3 2 , 2 ; z 2 ) ] .

The width of the Voigt profile

The full width at half maximum (FWHM) of the Voigt profile can be found from the widths of the associated Gaussian and Lorentzian widths. The FWHM of the Gaussian profile is

f G = 2 σ 2 ln ( 2 ) .

The FWHM of the Lorentzian profile is

f L = 2 γ .

Define Φ = f L / f G . Then the FWHM of the Voigt profile ( f V ) can be estimated as

f V f G ( 1 c 0 c 1 + ϕ 2 + 2 c 1 ϕ + c 0 2 c 1 2 ) ,

where c 0 = 2.0056 and c 1 = 1.0593. This estimate has a standard deviation of error of about 2.4% for values of φ between 0 and 10. Note that the above equation is exactly correct in the limit of φ = 0 and φ = ∞, that is for pure Gaussian and Lorentzian profiles.

A better approximation with an accuracy of 0.02% is given by

f V 0.5346 f L + 0.2166 f L 2 + f G 2 .

This approximation is exactly correct for a pure Gaussian, but has an error of about 0.000305% for a pure Lorentzian profile.

The uncentered Voigt profile

If the Gaussian profile is centered at μ G and the Lorentzian profile is centered at μ L , the convolution is centered at μ G + μ L and the characteristic function is

φ f ( t ; σ , γ , μ G , μ L ) = e i ( μ G + μ L ) t σ 2 t 2 / 2 γ | t | .

The mode and median are both located at μ G + μ L .

Voigt functions

The Voigt functions U, V, and H (sometimes called the line broadening function) are defined by

U ( x , t ) + i V ( x , t ) = π 4 t e z 2 erfc ( z ) = π 4 t w ( i z ) , H ( a , u ) = U ( u / a , 1 / 4 a 2 ) a π ,

where

z = ( 1 i x ) / 2 t ,

erfc is the complementary error function, and w(z) is the Faddeeva function.

Relation to Voigt profile

V ( x ; σ , γ ) = H ( a , u ) / ( 2 π σ ) ,

with

a = γ / ( 2 σ )

and

u = x / ( 2 σ ) .

Pseudo-Voigt approximation

The pseudo-Voigt profile (or pseudo-Voigt function) is an approximation of the Voigt profile V(x) using a linear combination of a Gaussian curve G(x) and a Lorentzian curve L(x) instead of their convolution.

The pseudo-Voigt function is often used for calculations of experimental spectral line shapes.

The mathematical definition of the normalized pseudo-Voigt profile is given by

V p ( x ) = η L ( x , f ) + ( 1 η ) G ( x , f ) with 0 < η < 1 .

There are several possible choices for the η parameter. A simple formula, accurate to 1%, is

η = 1.36603 ( f L / f ) 0.47719 ( f L / f ) 2 + 0.11116 ( f L / f ) 3 ,

where

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 .

References

Voigt profile Wikipedia