Voigt profile

From Wikipedia, the free encyclopedia

(Centered) Voigt
Probability density function
Plot of the centered Voigt profile for four cases
Plot of the centered Voigt profile for four cases. Each case has a full width at half-maximum of very nearly 3.6. The black and red profiles are the limiting cases of the Gaussian (γ =0) and the Lorentzian (σ =0) profiles respectively.
Cumulative distribution function
Centered Voigt CDF.
Parameters γ,σ > 0
Support x\in(-\infty,\infty)
Probability density function (pdf) \frac{\Re[w(z)]}{\sigma\sqrt{2\pi}},               ~~~z=\frac{x+i\gamma}{\sigma\sqrt{2}}
Cumulative distribution function (cdf) (complicated - see text)
Mean (not defined)
Median 0
Mode 0
Variance (not defined)
Skewness (not defined)
Excess Kurtosis (not defined)
Entropy
mgf (not defined)
Char. func. e^{-\gamma|t|-\sigma^2 t^2/2}

In spectroscopy, the Voigt profile is a spectral line profile named after Woldemar Voigt and found in all branches of spectroscopy in which a spectral line is broadened by two types of 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.

All normalized line profiles can be considered to be probability distributions. The Gaussian profile is equivalent to a Gaussian or normal distribution and a Lorentzian profile is equivalent to 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;\sigma,\gamma)=\int_{-\infty}^\infty G(x';\sigma)L(x-x';\gamma)\, dx'

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

G(x;\sigma)\equiv\frac{e^{-x^2/2\sigma^2}}{\sigma \sqrt{2\pi}}

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

L(x;\gamma)\equiv\frac{\gamma}{\pi(x^2+\gamma^2)}.

The defining integral can be evaluated as:

V(x;\sigma,\gamma)=\frac{\textrm{Re}[w(z)]}{\sigma\sqrt{2 \pi}}

where Re[w(z) ] is the real part of the complex error function of z  and

z=\frac{x+i\gamma}{\sigma\sqrt{2}}.

Contents

[edit] Properties

The Voigt profile is normalized:

\int_{-\infty}^\infty V(x;\sigma,\gamma)\,dx = 1

since it is the 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:

\varphi_f(t;\sigma,\gamma) = E(e^{ixt}) = e^{-\sigma^2t^2/2 - \gamma |t|}.

[edit] Cumulative distribution function

Using the above definition for z , the CDF can be found as follows:

F(x_0;\mu,\sigma) =\int_{-\infty}^{x_0} \frac{\mathrm{Re}(w(z))}{\sigma\sqrt{2\pi}}\,dx =\mathrm{Re}\left(\frac{1}{\sqrt{\pi}}\int_{z(-\infty)}^{z(x_0)} w(z)\,dz\right)

Substituting the definition of the complex error function yields for the indefinite integral:

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

Which may be solved (see e.g. the Mathematica integrator) to yield:

\frac{1}{\sqrt{\pi}}\int w(z)\,dz = \frac{\mathrm{erf}(z)}{2} +\frac{iz^2}{\pi}\,_2F_2\left(1,1;\frac{3}{2},2;-z^2\right)

where \,_2F_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:

F(x;\mu,\sigma)=\mathrm{Re}\left[\frac{1}{2}+ \frac{\mathrm{erf}(z)}{2} +\frac{iz^2}{\pi}\,_2F_2\left(1,1;\frac{3}{2},2;-z^2\right)\right]

[edit] 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_\mathrm{G}=2\sigma\sqrt{2\ln(2)}.\,

The FWHM of the Lorentzian profile is just fL = 2γ. Define φ = fL / fG. Then the FWHM of the Voigt profile (fV ) can be estimated as:

f_\mathrm{V}\approx f_\mathrm{G}\left(1-c_0c_1+\sqrt{\phi^2+2c_1\phi+c_0^2c_1^2}\right)

where c0 = 2.0056 and c1 = 1.0593. This estimate will have a standard deviation of error of about 2.4 percent for values of φ between 0 and 10. Note that the above equation will have the proper behavior in the limit of φ = 0 and φ = ∞.

[edit] The uncentered Voigt profile

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

\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|}.

The mode and median will then both be located at μG + μL.

[edit] See also

Image:Bvn-small.png Probability distributionsview  talk  edit ]
Univariate Multivariate
Discrete: BernoullibinomialBoltzmanncompound PoissondegenerateGauss-Kuzmingeometrichypergeometriclogarithmicnegative binomialparabolic fractalPoissonRademacherSkellamuniformYule-SimonzetaZipfZipf-Mandelbrot Ewensmultinomial
Continuous: BetaBeta primeCauchychi-squareDirac delta functionErlangexponentialexponential powerFfadingFisher's zFisher-TippettGammageneralized extreme valuegeneralized hyperbolicgeneralized inverse GaussianHalf-LogisticHotelling's T-squarehyperbolic secanthyper-exponentialhypoexponentialinverse chi-squareinverse Gaussianinverse gammaKumaraswamyLandauLaplaceLévyLévy skew alpha-stablelogisticlog-normalMaxwell-BoltzmannMaxwell speednormal (Gaussian)ParetoPearsonpolarraised cosineRayleighrelativistic Breit-WignerRiceStudent's ttriangulartype-1 Gumbeltype-2 GumbeluniformVoigtvon MisesWeibullWigner semicircleWilks' lambda DirichletKentmatrix normalmultivariate normalvon Mises-FisherWigner quasiWishart
Miscellaneous: Cantorconditionalexponential familyinfinitely divisiblelocation-scale familymarginalmaximum entropyphase-typeposteriorpriorquasisamplingsingular
In other languages