arxiv: v1 [math-ph] 15 May 2008

Transcription

1 arxiv:85.74v [math-ph 5 May 8 A note on the Voigt profile function G. PAGNINI and R.K. SAXENA ENEA, Centre Ezio Clementel, via Martiri di Monte Sole 4, I-49 Bologna, Italy gianni.pagnini@bologna.enea.it Department of Mathematics and Statistics, Jan Narain Vyas University, Jodhpur 345, India ram.saxena@yahoo.com May 5, 8 Abstract A Voigt profile function emerges in several physical investigations (e.g. atmospheric radiative transfer, astrophysical spectroscopy, plasma waves and acoustics) and it turns out to be the convolution of the Gaussian and the orentzian densities. Its relation with a number of special functions has been widely derived in literature starting from its Fourier type integral representation. The main aim of the present paper is to introduce the Mellin-Barnes integral representation as a useful tool to obtain new analytical results. Here, starting from the Mellin-Barnes integral representation, the Voigt function is expressed in terms of the Fox H-function which includes representations in terms of the Meijer G-function and previously well-known representations with other special functions. PACS numbers: 3.7.-n,.3.Gp,.3.Uu Introduction During the recent past, representations of the Voigt profile function in terms of special functions have been discussed in this Journal [,, 3, 4, 5, 6. With the present paper we aim to continue those researches. In particular, we derive the representation in terms of Fox H-function, which includes previously well-known results. Submitted to: J. Phys. A: Math. Gen.

2 The Voigt profile function emerges in several physical investigations as atmospheric radiative transfer, astrophysical spectroscopy, plasma waves and acoustics and molecular spectroscopy in general. Mathematically, it turns out to be the convolution of the Gaussian and the orentzian densities. Here we are studying the ordinary Voigt function and not its mathematical generalizations, for example, see the papers [7, 8, 9. The computation of the Voigt profile is an old issue in literature and many efforts are directed to evaluate this function with different techniques. In fact, an analytical explicit representation in terms of elementary functions does not exist and it can be considered a special function itself. Moreover it is strictly related to the plasma dispersion function [ and to a number of special functions as, for example, the confluent hypergeometric function, the complex complementary error function, the Dawson function, the parabolic cylinder function and the Whittaker function, see e.g. [,, 3,,, 3, 4. All previous representations are derived starting from the integral formula due to Reiche in 93 [5 that is actually a Fourier type integral. The Voigt profile function remains nowadays a mathematically and computationally interesting problem because computing profiles with high accuracy is still an expensive task. The actual interest on this topic is proven by several recent papers (7) on mathematical [6, 7, 8, 7, 9, and numerical aspects [,, 3. A huge collection exists of published works on this topic, but instead to report it, we give the significative datum that searching in Google Scholar the strings voigt profile function and voigt function the number of files found is 4, 6 and 7,, respectively. The Mellin-Barnes integrals are a family of integrals in the complex plane whose integrand is given by the ratio of products of Gamma functions. Despite of the name, the Mellin-Barnes integrals were initially studied in 888 by the Italian mathematician S. Pincherle [4, 5 in a couple of papers on the duality principle between linear differential equations and linear difference equations with rational coefficients. The Mellin-Barnes integrals are strongly related with the Mellin transform, in particular with the inverse transformation. As shown by O.I. Marichev [6, the problem to evaluate integrals can be successfully faced with a powerful method mainly based on their reduction to functions whose Mellin transform is the ratio of product of Gamma functions and then, after the inversion, the problem consists in the evaluation of Mellin-Barnes integrals. Moreover, they are also the essential tools for treating higher transcendental functions as Fox H-function and Meijer G-function and a useful representation to compute asymptotic behaviour of functions [7. The main object of the present paper is to introduce the Mellin-Barnes integral representation as a useful tool to obtain new analytical results that in the future can lead to efficient numerical algorithms for the Voigt function. A successful application of such approach has been shown in [8, where the parametric evolution equation of the Voigt function (and its probabilistic

3 generalization) is derived and the scaling laws, with respect to the parameter, in the asymptotic regimes are computed using the Mellin-Barnes integral representation. After all, this work can be seen also as an interesting exercise per se of application of the Mellin-Barnes integral method. The rest of the paper is organized as follows. In section the basic definition of the Voigt profile function is given and some classical and recent representations are reviewed. In section 3 the Mellin-Barnes integral representation of the Voigt function is derived and, in section 4 starting from this result, first the Voigt function is expressed in terms of the Fox H-function and later, in cascade, the representation with the Meijer G-function and other special functions are obtained. Finally in section 5 the summary and conclusions are given. The Voigt profile function. Basic definition The Gaussian G(x) and the orentzian (x) profiles are defined as [ G(x) = ( ) x exp, (x) = ω πωg ω G πω x + ω, () where ω G and ω are the corresponding widths. The variable x is a wavenumber and then its physical dimension is a length raised up to. The Voigt profile V (x) is obtained by the convolution of G(x) and (x) e (ξ/ω G) V (x) = (x ξ)g(ξ)dξ = ω /ω G π 3/ (x ξ) + ω dξ. () The comparison between the Voigtian, Gaussian and orentzian profile is shown in figure for different values of the width ratio ω /ω G. et f(κ) be the Fourier transform of f(x) so that then f(κ) = e +iκx f(x)dx, f(x) = π e iκx f(κ)dκ, (3) V (κ) = Ĝ(κ) (κ) = e ω G κ /4 e ω κ, (4) and V (x) = π = π e iκx e ω G κ /4 ω κ dκ e ω κ ω G κ /4 cos(κx)dκ. (5) 3

4 a=. a=.. Voigt Gauss orentz. Voigt Gauss orentz a= a= Voigt Gauss orentz Voigt Gauss orentz Figure : Comparison between the Voigtian, Gaussian and orentzian profiles with a = ω /ω G =.,.,,. 4

5 Formula (5) is the integral representation due to Reiche [5.. Some classical and recent representations et x be the dimensionless variable x = x/ω G, the Voigt function can be re-arranged in the form V (x) = πωg K(x,y), K(x,y) = y π e ξ (x ξ) + ydξ, (6) where y = ω /ω G and from (3) it follows that K(x,y) = e yξ ξ /4 cos(xξ)dξ. (7) π The Voigt function does not possess an explicit representation in terms of elementary functions and several alternatives to () have been given in literature, mainly with the intention to obtain a more efficient numerical computation. Combining x and y in the complex variable z = x iy, the function K(x,y) (6) is K(x,y) = Re[W(z), W(z) = i π e ξ dξ. (8) z ξ where W(z) is strongly related to the plasma dispersion function [ and some classical representations can be found, see e.g. [,, 3. In fact, form the relation of W(z) with the complex complementary error function Erfc( iz) and the Dawson function F(z) = e z z eξ dξ then it follows that K(x,y) = Re[W(z), W(z) = e z Erfc( iz), y >, (9) K(x,y) = Re[W(z), W(z) = e z + i π F(z). () More recent representations are, for example, those derived in by Di Rocco et al [9 K(x,y) = { ( ( ) n n + Γ(n + ) F, ) ;y n= a Γ ( ) n+ F (n +, 3 ) } ;y x n, () 5

6 where F (α,β;z) is the confluent hypergeometric function, and in 7 by Zaghloul [, which completed a previous formula given by Roston & Obaid [3, K(x,y) = [ erf(y)e ( x +y ) cos(xy)+ x e ( x +ξ ) sin[y(x ξ)dξ. π () Recently (7), He & Zhang [7 claimed to have derived an exact calculation of the Voigt profile that is proportional to the product of an exponential and a cosine function. However this representation assumes negative values in contrast with the non-negativity of the Voigt function. For this reason, this result is not correct. Further representations are given in terms of special functions, see for example the one involving the confluent hypergeometric function F [, 3 K(x,y) = e (y x ) cos(xy) π { (y + ix) F (;3/;(y + ix) )+ (y ix) F (;3/;(y ix) ) }, (3) and others involving the Whittaker function W k,m, the Erfc function and the parabolic cylinder function [4, formulae (7,3,6) K(x,y) = { (y ix) / e (y ix)/ W π /4, /4 ((y ix) )+ } (y + ix) / e (y+ix)/ W /4, /4 ((y + ix) ), (4) K(x,y) = { } e (y ix) Erfc(y ix) + e (y+ix) Erfc(y + ix), (5) K(x,y) = e(y x )/ {e ixy D [ (y ix) + e +ixy D [ } (y + ix). π 3 The Mellin-Barnes integral representation (6) et us consider again dimensional variables and the Gauss, orentz and Voigt functions defined as in (-). From (4) we have V (x) = π e iκx e ω G κ /4 ω κ dκ = π {I +(x) + I (x)}, (7) 6

7 where I ± (x) = e (ω ±ix)κ ω G κ /4 dκ. Following [7, the Mellin-Barnes integral representations of I ± (x) can be obtained from the definition of the Gamma function Γ(z) = ξ z e ξ dξ, and the Mellin-Barnes integral representation of e z (z ) e z = Γ(s)z s ds = πi ( ) n z n, n! where denotes a loop in the complex s plane which encircles the poles of Γ(s) (in the positive sense) with endpoints at infinity in Re(s) < and with no restrictions on arg z [7. The functions I ± (x) have the following Mellin-Barnes integral representations I ± (x) = = πi { e ω G κ /4 πi = ω G πi Γ(s) { Γ(s)Γ n= } Γ(s)[(ω ± ix)κ s ds dκ } (ω ± ix) s dκ e ω G κ /4 κ s ds ( s )[ s (ω ± ix) ds. (8) ω G Hence the Mellin-Barnes integral representation of the Voigt function is given by { V (x) = πω G πi πi ( Γ(s)Γ s )[ s (ω + ix) ds + ω G ( Γ(s)Γ s )[ s (ω ix) ds}. (9) ω G However we note that for a complex number z = z e iθ, θ = arctan(im(z)/re(z)), the following rule holds z n + z n = z n e inθ + z n e inθ = z n (e inθ + e inθ ) = z n cos(nθ) = z n cos(n arctan(im(z)/re(z))), 7

8 and formula (9) becomes V (x) = ( Γ(s)Γ πω G πi s ) [ ( x cos s arctan )(4 ω + ) s/ x ω ωg ds. () Consider again (8), changing s s and taking the corresponding integration path as a loop in the complex plane that encircles the poles of Γ( s), an other Mellin-Barnes integral representation of the Voigt function equivalent to (9) is V (x) = { πω G πi πi and in more compact form ( Γ( s)γ + s )[ s (ω + ix) ds + ω G ( Γ( s)γ + s )[ s (ω ix) ds}, () ω G V (x) = ( Γ( s)γ πω G πi + s ) [ ( x cos s arctan )(4 ω + ) s/ x ω ωg ds. () 4 The Fox H- and Meijer G-function representations 4. The Fox H-function representation The Voigt function can be represented in terms of well known special function, see., but its representation in terms of Fox H-function is still not known. The expression in terms of Fox H-function is important because it is the most modern representation method and, actually, it is the most compact form to represent higher transcendental functions. The definition of the Fox H-function is given in Appendix. For ω G and ω fixed, (ω + ix)/ω G and (ω ix)/ω G, from () and (9) we have, respectively, (/,/) (,) (ω ix) ω G (/,/) (,) V (x) = { [ H (ω + ix) πω G ω G [ H + }, (3) 8

9 V (x) = { [ H ω G πω G (ω + ix) [ H ω G (ω ix) (,) (/, /) (,) (/, /) + }, (4) As a consequence of the fact that this is the most comprehensive representation, in cascade, the others with less general functions can be obtained. 4. The Meijer G-function representation The first H-function in (3) can be rewritten in terms of the Meijer G- function. Setting Z = (ω + ix)/ω G, we have [ H Z (/,/) (,) = πi = πi = π πi ( Γ( s)γ + s ) Z s ds ( ) Γ( s)γ + s Z s ds ( ) ( )( Z ) s Γ( s)γ s Γ + s ds 4 = [ Z G / π 4, (5),/ where the change of variable s s and the duplication rule for the Gamma function, Γ(z) = Γ(z)Γ(/ + z) z π /, are applied. The second H- function in (3) follows from the first with Z = (ω ix)/ω G. Finally, the Voigt function in terms of the Meijer G-function is given by V (x) = π 3/ ω G { [ Z G 4 /,/ [ Z + G 4 /,/ }. (6) The Meijer G-function can be reduced to other special functions and, for example, representations given in (4, 5, 6) are straightforwardly recovered. In fact, the Meijer G-function and the Whittaker function W k,m are related by the formula [3, p. 435 [ G z a b,c = Γ(b a + )Γ(c a + )x (b+c )/ e x/ W k,m (k), where k = a (b + c + )/ and m = (b c)/. Thus the G-functions in (6) are expressed in terms of the Whittaker function as 9

10 [ Z G 4 /,/ = π ( ) Z / e Z /8 W /4, /4(Z /4). Hence, the Voigt function in terms of the Whittaker function is { (Z V (x) = ) / e Z /8 W πω G /4, /4(Z /4)+ ( ) / Z e Z /8 W /4, /4( /4)} Z. (7) Moreover, from the identity [ Z G / 4,/ we have and from we have V (x) = V (x) = = π e Z /4 Erfc(Z/), { e Z /4 Erfc(Z/) + e Z /4 Erfc( πω Z/) }, (8) G [ Z G 4 /,/ = π e Z /8 D (Z/ ), {e Z /8 D (Z/ ) + e Z /8 D ( Z/ } ). (9) πωg Setting ω G = in (7,8,9) formulae (4,5,6) are recovered, respectively. 5 Summary and conclusions In the present paper the Mellin-Barnes integral representation of the Voigt profile function is derived. We think that this integral representation is a useful tool to have new analytical and numerical results in the subject. Starting from this, the Voigt function has been expressed in terms of the Fox H-function, which is the most comprehensive representation. In cascade, the expression in terms of the Meijer G-function is obtained and the previous well-known representations with the Whittaker, the Erf c and the Parabolic cylinder functions are recovered.

11 Appendix. The Fox H- and the Meijer G- functions The Fox H-function is defined as [3, 33 [ Hpq mn z (a,a ),,(a p,a p ) (b,b ),,(b q,b q ) where = h(s)z s ds, πi h(s) = m j= Γ(b j B j s) n j= Γ( a j + A j s) q j=m+ Γ( b j + B j s) p j=n+ Γ(a j A j s), an empty product is interpreted as unity, {m, n, p, q} are nonnegative integers so that m q, n p, {A j, B j } are positive numbers and {a j, b j } complex numbers. The H-function is an analytic function of z and makes sense if [3, 33 where i) z, µ > ii) < z < β, µ =, µ = q B j j= p A j, β = j= p j= A A j j q j= B B j j. The Meijer G-function corresponds to the special case A j = B k = (j =,...,p;k =,...,q) [ Hpq mn z (a [,),,(a p,) = G mn pq z (b,),,(b q,) a,,a p. b,,b q References [ Exton H 98 J. Phys. A: Math. Gen [ Katriel J 98 J. Phys. A: Math. Gen [3 Fettis H E 983 J. Phys. A: Math. Gen [4 Keshavamurthy R S 987 J. Phys. A: Math. Gen. 73 [5 Sajo E 993 J. Phys. A: Math. Gen [6 Gubner J A 994 J. Phys. A: Math. Gen [7 Sampoorna M, Nagendra K N and Frisch H 7 J. Quant. Spectrosc. Radiat. Transfer 4 7

12 [8 Srivastava H M and Miller E A 987 Astrophys. Space Sci. 35 [9 Srivastava H M and Chen M P 99 Astrophys. Space Sci [ Fried B D and Conte S D 96 The Plasma Dispersion Function (New York: Academic) [ Armstrong B H 967 J. Quant. Spectrosc. Radiat. Transfer 7 6 [ Schreier F 99 J. Quant. Spectrosc. Radiat. Transfer [3 Shippony Z and Read W G 993 J. Quant. Spectrosc. Radiat. Transfer [4 Yang S 994 Int. J. Math. Educ. Sci. Technol [5 Reiche F 93 Ber. Deutsch. Phys. Ges. 5 3 [6 Dulov E N and Khripunov D M 7 J. Quant. Spectrosc. Radiat. Transfer 7 4 [7 He J and Zhang Q 7 J. Opt. A [8 Petrov E P 7 J. Quant. Spectrosc. Radiat. Transfer 3 7 [9 Srivastava H M and Pogány T K 7 Russian J. Math. Phys [ Zaghloul M R 7 Mon. Not. R. Astron. Soc [ eiweke R J 7 J. Quant. Spectrosc. Radiat. Transfer [ etchworth K and Benner D C 7 J. Quant. Spectrosc. Radiat. Transfer 7 73 [3 Mendenhall M H 7 J. Quant. Spectrosc. Radiat. Transfer 5 59 [4 Pincherle S 888 Atti R. Accad. incei, Rend. Cl. Sci. Fis. Mat. Natur. 4 (4) 694 and 79 [Reprinted in UMI (Unione Matematica Italiana) 954 Salvatore Pincherle: Opere Scelte vol (Roma: Cremonese) p 3 and 3 [5 Mainardi F and Pagnini G 3 J. Comput. Appl. Math [6 Marichev O I 98 Handbook of Integral Transforms of Higher Transcendental Functions: Theory and Algorithmic Tables (Chichester: Ellis Horwood imited) [7 Paris R B and Kaminski D Asymptotics and Mellin-Barnes Integrals (Cambridge: Cambridge University Press) [8 Pagnini G and Mainardi F 7 arxiv:7.446 Submitted

13 [9 Di Rocco H O, Iriarte D I and Pomarico J Applied Spectroscopy 55 8 [3 Roston G D and Obaid F S 5 J. Quant. Spectrosc. Radiat. Transfer [3 Erdélyi A, Magnus W, Oberhettinger F and Tricomi F 954 Tables of Integral Transforms vol New York: McGraw-Hill) [3 Mathai A M and Saxena R K 978 The H-function with Applications in Statistics and Other Disciplines (New Delhi: Wiley Eastern imited) [33 Srivastava H M, Gupta K C and Goyal S P 98 The H-functions of One and Two Variables with Applications (New Delhi: South Asian) 3