Voigt profile

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

${\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).}$

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

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

which may be solved to yield

${\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),}$

where ${\displaystyle {}_{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:

${\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].}$

If the Gaussian profile is centered at ${\displaystyle \mu _{G}}$ and the Lorentzian profile is centered at ${\displaystyle \mu _{L}}$, the convolution is centered at ${\displaystyle \mu _{G}+\mu _{L}}$ and the characteristic function is

${\displaystyle \varphi _{f}(t;\sigma ,\gamma ,\mu _{\mathrm {G} },\mu _{\mathrm {L} })=e^{i(\mu _{\mathrm {G} }+\mu _{\mathrm {L} })t-\sigma ^{2}t^{2}/2-\gamma |t|}.}$

The mode and median are both located at ${\displaystyle \mu _{G}+\mu _{L}}$.

The first and second derivative profiles can be expressed in terms of the Faddeeva function as follows:

${\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 ^{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 }},}$

using the above definition for z.

${\displaystyle V(x;\sigma ,\gamma )=H(a,u)/({\sqrt {2}}{\sqrt {\pi }}\sigma ),}$

with

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

and

${\displaystyle u=x/({\sqrt {2}}\sigma ).}$

The Tepper-García function, named after German-Mexican Astrophysicist Thor Tepper-García, is a combination of an exponential function and rational functions that approximates the line broadening function ${\displaystyle H(a,u)}$ over a wide range of its parameters.[2] It is obtained from a truncated power series expansion of the exact line broadening function.

In its most computationally efficient form, the Tepper-García function can be expressed as

${\displaystyle T(a,u)=R-\left(a/{\sqrt {\pi }}P\right)~\left[R^{2}~(4P^{2}+7P+4+Q)-Q-1\right]\,,}$

where ${\displaystyle P\equiv u^{2}}$, ${\displaystyle Q\equiv 3/(2P)}$, and ${\displaystyle R\equiv e^{-P}}$.

Thus the line broadening function can be viewed, to first order, as a pure Gaussian function plus a correction factor that depends linearly on the microscopic properties of the absorbing medium (encoded in ${\displaystyle a}$); however, as a result of the early truncation in the series expansion, the error in the approximation is still of order ${\displaystyle a}$, i.e. ${\displaystyle H(a,u)\approx T(a,u)+{\mathcal {O}}(a)}$. This approximation has a relative accuracy of

${\displaystyle \epsilon \equiv {\frac {\vert H(a,u)-T(a,u)\vert }{H(a,u)}}\lesssim 10^{-4}}$

over the full wavelength range of ${\displaystyle H(a,u)}$, provided that ${\displaystyle a\lesssim 10^{-4}}$. In addition to its high accuracy, the function ${\displaystyle T(a,u)}$ is easy to implement as well as computationally fast. It is widely used in the field of quasar absorption line analysis.[3]

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

${\displaystyle V_{p}(x,f)=\eta \cdot L(x,f)+(1-\eta )\cdot G(x,f)}$ with ${\displaystyle 0<\eta <1}$.

${\displaystyle \eta }$ is a function of full width at half maximum (FWHM) parameter.

There are several possible choices for the ${\displaystyle \eta }$ parameter.[4][5][6][7] A simple formula, accurate to 1%, is[8][9]

${\displaystyle \eta =1.36603(f_{L}/f)-0.47719(f_{L}/f)^{2}+0.11116(f_{L}/f)^{3},}$

where now, ${\displaystyle \eta }$ is a function of Lorentz (${\displaystyle f_{L}}$), Gaussian (${\displaystyle f_{G}}$) and total (${\displaystyle f}$) Full width at half maximum (FWHM) parameters. The total FWHM (${\displaystyle f}$) parameter is described by:

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

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

${\displaystyle f_{\mathrm {G} }=2\sigma {\sqrt {2\ln(2)}}.}$

The FWHM of the Lorentzian profile is

${\displaystyle f_{\mathrm {L} }=2\gamma .}$

A rough approximation for the relation between the widths of the Voigt, Gaussian, and Lorentzian profiles is:

${\displaystyle f_{\mathrm {V} }\approx f_{\mathrm {L} }/2+{\sqrt {f_{\mathrm {L} }^{2}/4+f_{\mathrm {G} }^{2}}}.}$

This approximation is exactly correct for a pure Gaussian.

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

${\displaystyle f_{\mathrm {V} }\approx 0.5346f_{\mathrm {L} }+{\sqrt {0.2166f_{\mathrm {L} }^{2}+f_{\mathrm {G} }^{2}}}.}$

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