Literature DB >> 28179785

Tables of the Inverse Laplace Transform of the Function [Formula: see text].

Menachem Dishon1, John T Bendler2, George H Weiss3.   

Abstract

The inverse transform, [Formula: see text], 0 < β < 1, is a stable law that arises in a number of different applications in chemical physics, polymer physics, solid-state physics, and applied mathematics. Because of its important applications, a number of investigators have suggested approximations to g(t). However, there have so far been no accurately calculated values available for checking or other purposes. We present here tables, accurate to six figures, of g(t) for a number of values of β between 0.25 and 0.999. In addition, since g(t), regarded as a function of β, is uni-modal with a peak occurring at t = tmax we both tabulate and graph tmax and 1/g(tmax) as a function of β, as well as giving polynomial approximations to 1/g(tmax).

Entities:  

Keywords:  numerical inversion of Laplace transforms; relaxation processes; stable laws; stretched exponentials

Year:  1990        PMID: 28179785      PMCID: PMC4922407          DOI: 10.6028/jres.095.036

Source DB:  PubMed          Journal:  J Res Natl Inst Stand Technol        ISSN: 1044-677X


1. Introduction

It has been known for at least 150 years that mechanical relaxation in solids is non-exponential, the decay often being characterized by a fractional power-law or logarithmic function [1,2]. It is also now generally recognized that all glassy materials exhibit non-exponential relaxation behavior both above and below the glass transition temperature, Tg. This is especially clear from measurements obtained from mechanical [3-6], dielectric [7-9], and photon correlation spectroscopy [10,11]. It is also seen in measurements of volumetric [12], and thermal response [13,14]. In recent years theorists have become interested in the possibility that complex disordered systems exhibit universal features in their relaxation and transport properties, possibly arising from self-similar arrangements of obstacles to motion. This has been particularly encouraged by the observation that nearly all glassy relaxation phenomena can be described by the Kohlrausch-Williams-Watts (KWW) function In many physical applications it is convenient to represent ϕ(t; τ) in the form of a Laplace transform, which we write as where Thus, the function h(v) can be found as an inverse Laplace transform of the function ϕ(t;τ). The function h(v) has found application in the context of the theory of trap-controlled hopping in solid state physics [15,16], chromatography [17], and in the study of models for transport in disordered media [18], as well as in the deconvolution of noisy data [19]. A number of approximate algorithms have been proposed in the literature of chemical physics for the numerical evaluation of h(v) [20-25], in addition to a representation of h(v) in terms of a convergent series given by Pollard [26]. Without loss of generality we can set τ = 1 since h(v) can be represented in terms of the inverse transform where Γ is a line to the right of the origin and parallel to the imaginary axis. The convergent series given by Pollard is In an earlier paper we have presented an accurate tabulation of the sine and cosine transforms of the function exp(− t), needed for the analysis of measurements of dielectric properties taken as a function of frequency [27]. In the present paper we tabulate the inverse Laplace transform h(v). These tables may be used directly for the analysis of experimental data, but are also intended for use as a check on more easily programmed approximations, such as those suggested by earlier investigators [28-30].

2. Numerical Analysis

Two techniques were used to generate the tables that follow which provide an internal check on the accuracy of the computation. The first is that of numerical inversion of the Laplace transform, using a method first suggested by Dubner and Abate [31], and later given in an improved version by Crump [32]. The second is that of direct evaluation of the series given in eq (5). The approximate inverse of a Laplace transform ĝ(s) = ℒ{g(t)} can be expressed in the form of a Fourier series: with an error, E(t) = g(t) − g(t), given by The function can be written in terms of the parameters as In eqs (6) and (7) the constants a and T are arbitrary and can be chosen to maximize accuracy in any particular application. In the present instance, in which , the choice of these parameters is quite straightforward as will be shown below. Equation (6) was used to evaluate the inverse transform of ĝ(s) for values of β in the range 0.20≤β≤0.999 and values of t ranging from 10−8 (for selected values of β) to 5, to an accuracy of at least nine significant digits. In these ranges of β and t the choice of parametric ranges aϵ(2.5,5) and Tϵ(4,8) sufficed to produce the stated accuracy. The accuracy of the numerical inversion can be checked in detail for three cases in which the inverse transforms are known exactly, where Ai(x) is an Airy function and U(x, y, z) is a confluent hypergeometric function [29]. Typical results for the relative error are given in table 1.
Table 1

Relative errors in the numerical inversion of ĝ(s) for β = 1/3, 1/2, and 2/3 for different values of a and T

β1/31/22/3
aT  t
2.540.12(−11)3(−11)1(−6)
2.541.02(−10)1(−10)4(−11)
2.545.06(−10)5(−10)4(−10)
2.580.12(−20)2(−20)8(−16)
5.040.013(−20)1(−11)
5.040.14(−20)7(−20)3(−15)
5.041.03(−19)2(−19)9(−20)
5.045.01(−18)1(−18)8(−19)
An alternative approach to the evaluation of g(t) is through the direct series shown in eq (5). The form of the series renders it useful for finding g(t) for large t, but the utility of the series form has occasionally been dismissed because of numerical problems associated with convergence at smaller t. We encountered no difficulties in finding g(t) from eq (5), provided that we used a double precision routine for the gamma functions for k≤22 as well as a Padé correction to Stirling’s approximation at larger k [33]. Thus, we write where The F are constants with the values

3. Tables, Graphs, and Numerical Approximations

The inverse transform of the function ĝ(s) is tabulated in table 2 for the following values of β: 0.25(0.01)0.30(0.02)0.98, 0.99, 0.995, 0.997, 0.998, and 0.999. The finer intervals in β at low values of β are required because of the considerable changes in the function in that neighborhood. Spacings in t vary with β and t in such a way that the peaks of g(t) are most densely covered. There is little need to tabulate g(t) for t > 5 because for these values, the sum of no more than 10 terms of the series in eq (5) suffice to produce g(t) to six-digit accuracy for values of β in the interval (0.05,0.999). For example, if β=0.6 the sum of seven terms of the series gives g(10) to six places, and the sum of four terms gives g(100) to the same accuracy. Figures 1a–c contain graphs of g(t) as a function of t over the entire range of tabulated values of β. Note that for β = 1 g(t) = d(t − 1), a Dirac delta function which is represented as a vertical line in figure 1c.
Table 2

Inverse Laplace transform g(t) of

tβ
0.250.260.270.280.29
0.000010.370151D−040.515902D−060.248661D−08
0.000020.151806D−020.675577D−040.144824D−050.123255D−07
0.000030.853094D−020.649138D−030.275324D−040.560637D−060.451901D−08
0.000040.244726D−010.258459D−020.165894D−030.571718D−050.902078D−07
0.000050.507423D−010.672926D−020.575158D−030.284658D−040.711259D−06
0.000060.873213D−010.137312D−010.145291D−020.940782D−040.330246D−05
0.000070.133437D+000.239884D−010.300022D−020.239622D−030.109589D−04
0.000080.187991D+000.376838D−010.539788D−020.510882D−030.289237D−04
0.000090.249798D+000.548385D−010.879472D−020.958456D−030.647631D−04
0.000100.317705D+000.753611D−010.133042D−010.163414D−020.128251D−03
0.000110.390651D+000.990866D−010.190057D−010.258809D−020.231044D−03
0.000120.467687D+000.125806D+000.259484D−010.386679D−020.386250D−03
0.000130.547978D+000.155287D+000.341556D−010.551183D−020.607956D−03
0.000140.630804D+000.187288D+000.436295D−010.755922D−020.910770D−03
0.000150.715540D+000.221568D+000.543553D−010.100391D−010.130942D−02
0.000160.801653D+000.257891D+000.663051D−010.129761D−010.181839D−02
0.000170.888689D+000.296036D+000.794409D−010.163891D−010.245168D−02
0.000180.976259D+000.335790D+000.937176D−010.202922D−010.322255D−02
0.000190.106404D+010.376957D+000.109085D+000.246947D−010.414334D−02
0.000200.115175D+010.419354D+000.125488D+000.296019D−010.522547D−02
0.000210.123915D+010.462814D+000.142873D+000.350156D−010.647923D−02
0.000220.132605D+010.507183D+000.161181D+000.409343D−010.791385D−02
0.000230.141229D+010.552318D+000.180357D+000.473537D−010.953745D−02
0.000240.149772D+010.598093D+000.200343D+000.542673D−010.113571D−01
0.000250.158223D+010.644389D+000.221084D+000.616668D−010.133787D−01
0.000300.198831D+010.880304D+000.334310D+000.105556D+000.266587D−01
0.000350.236313D+010.111698D+010.459197D+000.159675D+000.453716D−01
0.000400.270528D+010.134838D+010.591057D+000.222165D+000.694073D−01
0.000450.301586D+010.157109D+010.726376D+000.291223D+000.984011D−01
0.000500.329700D+010.178326D+010.862582D+000.365234D+000.131854D+00
0.000600.378101D+010.217324D+010.113082D+010.522815D+000.209913D+00
0.000700.417683D+010.251762D+010.138668D+010.686433D+000.299254D+00
0.000800.450132D+010.281968D+010.162607D+010.850458D+000.396103D+00
0.000900.476810D+010.308392D+010.184753D+010.101129D+010.497430D+00
0.001000.498794D+010.331492D+010.205105D+010.116672D+010.600902D+00
0.001100.516936D+010.351690D+010.223731D+010.131544D+010.704769D+00
0.001200.531912D+010.369363D+010.240740D+010.145678D+010.807740D+00
0.001300.544260D+010.384838D+010.256249D+010.159047D+010.908884D+00
0.001400.554414D+010.398397D+010.270380D+010.171648D+010.100755D+01
0.001500.562724D+010.410281D+010.283253D+010.183497D+010.110327D+01
0.001600.569476D+010.420700D+010.294978D+010.194618D+010.119577D+01
0.001700.574907D+010.429831D+010.305657D+010.205044D+010.128485D+01
0.001800.579212D+010.437828D+010.315385D+010.214809D+010.137044D+01
0.001900.582551D+010.444824D+010.324247D+010.223950D+010.145249D+01
0.002000.585062D+010.450934D+010.332320D+010.232501D+010.153104D+01
0.002100.586859D+010.456258D+010.339676D+010.240499D+010.160613D+01
0.002200.588040D+010.460881D+010.346376D+010.247978D+010.167784D+01
0.002300.588685D+010.464881D+010.352477D+010.254971D+010.174627D+01
0.002400.588868D+010.468323D+010.358032D+010.261507D+010.181152D+01
0.002500.588646D+010.471266D+010.363086D+010.267617D+010.187372D+01
0.003000.583090D+010.480102D+010.382187D+010.292665D+010.214289D+01
0.003500.572936D+010.482136D+010.393647D+010.310461D+010.235284D+01
0.004000.560506D+010.480024D+010.399971D+010.322986D+010.251601D+01
0.004500.547051D+010.475348D+010.402773D+010.331631D+010.264234D+01
0.005000.533269D+010.469083D+010.403120D+010.337384D+010.273959D+01
0.006000.506134D+010.454054D+010.399133D+010.342877D+010.286945D+01
0.007000.480602D+010.437755D+010.391576D+010.343198D+010.293927D+01
0.008000.457076D+010.421433D+010.382270D+010.340442D+010.296975D+01
0.009000.435561D+010.405650D+010.372195D+010.335844D+010.297401D+01
0.010000.415919D+010.390643D+010.361892D+010.330154D+010.296059D+01
0.011000.397972D+010.376497D+010.351666D+010.323842D+010.293520D+01
0.012000.381537D+010.363219D+010.341686D+010.317209D+010.290172D+01
0.013000.366447D+010.350777d+010.332044D+010.310448D+010.286285D+01
0.014000.352550D+010.339122D+010.322787D+010.303685D+010.282046D+01
0.015000.339714D+010.328200D+010.313931D+010.297004D+010.277590D+01
0.016000.327824D+010.317955D+010.305479D+010.290456D+010.273013D+01
0.017000.316778D+010.308333D+010.297421D+010.284075D+010.268383D+01
0.018000.306491D+010.299284D+010.289744D+010.277880D+010.263752D+01
0.019000.296885D+010.290761D+010.282430D+010.271882D+010.259154D+01
0.020000.287894D+010.282722D+010.275460D+010.266086D+010.254614D+01
0.021000.279460D+010.275127D+010.268816D+010.260491D+010.250152D+01
0.022000.271531D+010.267942D+010.262478D+010.255095D+010.245779D+01
0.023000.264063D+010.261135D+010.256428D+010.249893D+010.241504D+01
0.024000.257015D+010.254676D+010.250649D+010.244880D+010.237332D+01
0.025000.250353D+010.248540D+010.245125D+010.240048D+010.233267D+01
0.030000.221824D+010.221946D+010.220816D+010.218367D+010.214538D+01
0.035000.199355D+010.200649D+010.200950D+010.200190D+010.198306D+01
0.040000.181164D+010.183191D+010.184418D+010.184784D+010.184226D+01
0.045000.166111D+010.168604D+010.170445D+010.171578D+010.171946D+01
0.050000.153434D+010.156223D+010.158475D+010.160139D+010.161165D+01
0.060000.133230D+010.136312D+010.139020D+010.141316D+010.143158D+01
0.070000.117811D+010.120973D+010.123868D+010.126467D+010.128738D+01
0.080000.105636D+010.108774D+010.111719D+010.114448D+010.116936D+01
0.090000.957662D+000.988292D+000.101751D+010.104515D+010.107100D+01
0.100000.875953D+000.905594D+000.934205D+000.961646D+000.987760D+00
0.110000.807149D+000.835699D+000.863499D+000.890438D+000.916393D+00
0.120000.748385D+000.775818D+000.802711D+000.828978D+000.854519D+00
0.130000.697593D+000.723924D+000.749875D+000.775379D+000.800357D+00
0.140000.653241D+000.678505D+000.703514D+000.728214D+000.752545D+00
0.150000.614166D+000.638412D+000.662499D+000.686386D+000.710024D+00
0.160000.579473D+000.602753D+000.625950D+000.649031D+000.671960D+00
0.170000.548460D+000.570828D+000.593171D+000.615466D+000.637686D+00
0.180000.520567D+000.542074D+000.563604D+000.585140D+000.606661D+00
0.190000.495343D+000.516040D+000.536798D+000.557605D+000.578444D+00
0.200000.472421D+000.492356D+000.512382D+000.532490D+000.552669D+00
0.210000.451498D+000.470716D+000.490048D+000.509490D+000.529034D+00
0.220000.432322D+000.450866D+000.469541D+000.488347D+000.507281D+00
0.230000.414684D+000.432591D+000.450645D+000.468847D+000.487195D+00
0.240000.398403D+000.415711D+000.433176D+000.450803D+000.468592D+00
0.250000.383330D+000.400072D+000.416980D+000.434059D+000.451313D+00
0.300000.322103D+000.336439D+000.350960D+000.365674D+000.380590D+00
0.350000.277427D+000.289905D+000.302564D+000.315413D+000.328462D+00
0.400000.243395D+000.254402D+000.265577D+000.276930D+000.288471D+00
0.450000.216613D+000.226431D+000.236401D+000.246533D+000.256837D+00
0.500000.194994D+000.203832D+000.212807D+000.221929D+000.231205D+00
0.550000.177181D+000.185199D+000.193341D+000.201614D+000.210025D+00
0.600000.162254D+000.169579D+000.177013D+000.184565D+000.192240D+00
0.650000.149568D+000.156298D+000.163126D+000.170058D+000.177100D+00
0.700000.138657D+000.144872D+000.151174D+000.157569D+000.164062D+00
0.750000.129174D+000.134940D+000.140783D+000.146708D+000.152721D+00
0.800000.120859D+000.126229D+000.131668D+000.137180D+000.142769D+00
0.850000.113510D+000.118529D+000.123610D+000.128756D+000.133970D+00
0.900000.106968D+000.111675D+000.116437D+000.121256D+000.126136D+00
0.950000.101110D+000.105537D+000.110012D+000.114538D+000.119118D+00
1.000000.958339D−010.100008D+000.104226D+000.108488D+000.112798D+00
1.100000.867153D−010.904539D−010.942257D−010.980323D−010.101875D+00
1.200000.791153D−010.824909D−010.858919D−010.893196D−010.927749D−01
1.300000.726869D−010.757563D−010.788447D−010.819528D−010.850815D−01
1.400000.671813D−010.699893D−010.728108D−010.756465D−010.784968D−01
1.500000.624152D−010.649976D−010.675892D−010.701902D−010.728008D−01
1.600000.582505D−010.606366D−010.630282D−010.654252D−010.678277D−01
1.700000.545814D−010.567953D−010.590116D−010.612299D−010.634503D−01
1.800000.513255D−010.533874D−010.554483D−010.575096D−010.595694D−01
1.900000.484176D−010.503443D−010.522683D−010.541892D−010.561066D−01
2.000000.458055D−010.476113D−010.494124D−010.512084D−010.529989D−01
2.100000.434467D−010.451439D−010.468347D−010.485187D−010.501953D−01
2.200000.413066D−010.429058D−010.444971D−010.460801D−010.476542D−01
2.300000.393567D−010.408669D−010.423681D−010.438597D−010.453410D−01
2.400000.375729D−010.390022D−010.404214D−010.418299D−010.432270D−01
2.500000.359352D−010.372906D−010.386350D−010.399677D−010.412880D−01
2.600000.344267D−010.357144D−010.369902D−010.382535D−010.395036D−01
2.700000.330329D−010.342582D−010.354710D−010.366707D−010.378563D−01
2.800000.317409D−010.329092D−010.340640D−010.352050D−010.363314D−01
2.900000.305414D−010.316562D−010.327573D−010.338442D−010.349159D−01
3.000000.294237D−010.304893D−010.315408D−010.325775D−010.335988D−01
3.200000.274043D−010.283816D−010.293441D−010.302911D−010.312219D−01
3.400000.256298D−010.265302D−010.274154D−010.282845D−010.291370D−01
3.600000.240589D−010.248919D−010.257093D−010.265103D−010.272943D−01
3.800000.226589D−010.234324D−010.241901D−010.249311D−010.256550D−01
4.000000.214039D−010.221246D−010.228292D−010.235171D−010.241877D−01
4.200000.202728D−010.209463D−010.216036D−010.222442D−010.228674D−01
4.400000.192485D−010.198796D−010.204945D−010.210927D−010.216735D−01
4.600000.183167D−010.189096D−010.194863D−010.200463D−010.205890D−01
4.800000.174657D−010.180240D−010.185661D−010.190916D−010.196000D−01
5.000000.166856D−010.172124D−010.177231D−010.182174D−010.186946D−01
Figure 1a

Curves of g(t) as a function of t in the neighborhood of the peak values for β=0.25, 0.275, and 0.30.

Figure 1b

Curves of g(t) as a function of t in the neighborhood of the peak values for β=0.35(0.1)0.95.

Figure 1c

Curves of g(t) as a function of fin the neighborhood of the peak values for β=0.99, 0.995, 0.997, 0.998, and 0.999. The delta function at β=1 is indicated by the vertical line.

It is evident, from the curves shown in figure 1, that the g(t) are unimodal. The position of the peak will be denoted by tmax. Table 3 contains some values of tmax and g(tmax) for the values of β for which we performed our tabulations. It is interesting to observe that among the values of g(tmax) there is a minimum value within the interval (0,1). Figure 2a shows graphs of tmax and 1/g(tmax) as functions of β for values of β between 0.15 and 1. The minimum of g(tmax) occurs at tmax=0.252+ and is equal to 0.888+. These values correspond to β=0.567+. Figure 2b contains a plot of 1/g(tmax) as a function of tmax. Finally, we have derived polynomial least-square approximations to 1/g(tmax) as a function of β. The coefficients of the approximating polynomials as well as a graphical indication of the degree of agreement with our more accurately calculated values of this function are shown in figure 3. A good approximation to 1/g(tmax) probably requires fitting some function other than a polynomial.
Table 3

tmax and g(tmax) as a function of β, where g(t) is the inverse Laplace transform of

βtmaxg(tmax)
0.150.00000190503450.9129
0.160.00000630739209.22449
0.170.0000177337108.862395
0.180.0000436172  62.162073
0.190.0000960414  38.311543
0.200.000192795  25.155645
0.210.000357987  17.416529
0.220.000622088  12.610041
0.230.00102137 9.484013
0.240.00159686 7.369062
0.250.00239285 5.888686
0.260.00345532 4.821523
0.270.00483018 4.032270
0.280.00656175 3.435337
0.290.00869134 2.974907
0.300.0112561 2.613581
0.310.0142884 2.325679
0.320.0178147 2.093179
0.330.0218560 1.903181
0.340.0264274 1.746284
0.350.0315381 1.615524
0.360.0371920 1.505665
0.370.0433880 1.412716
0.380.0501201 1.333598
0.390.0573788 1.265907
0.400.0651506 1.207748
0.410.0734193 1.157611
0.420.0821662 1.114283
0.430.0913707 1.076782
0.440.101011 1.044308
0.450.111063 1.016202
0.460.121503 0.991922
0.470.132308 0.971016
0.480.143453 0.953108
0.490.154914 0.937885
0.500.166667 0.925082
0.510.178689 0.914477
0.520.190959 0.905885
0.530.203454 0.899150
0.540.216155 0.894139
0.550.229042 0.890746
0.560.242097 0.888881
0.570.255304 0.888472
0.580.268646 0.889461
0.590.282109 0.891806
0.600.295680 0.895475
0.610.309346 0.900451
0.620.323098 0.906726
0.630.336924 0.914305
0.640.350818 0.923202
0.650.364772 0.933445
0.660.378780 0.945070
0.670.392837 0.958129
0.680.406941 0.972686
0.690.421088 0.988819
0.700.435279 1.006623
0.710.449512 1.026212
0.720.463791 1.047719
0.730.478118 1.071303
0.740.492496 1.097151
0.750.506933 1.125483
0.760.521435 1.156560
0.770.536011 1.190689
0.780.550671 1.228239
0.790.565430 1.269648
0.800.580300 1.315445
0.810.595300 1.366273
0.820.610450 1.422918
0.830.625773 1.486349
0.840.641295 1.557781
0.850.657048 1.638750
0.860.673068 1.731228
0.870.689397 1.837792
0.880.706087 1.961869
0.890.723197 2.108126
0.900.740799 2.283081
0.910.758982 2.496129
0.920.777856 2.761327
0.930.797561 3.100708
0.940.818279 3.550899
0.950.840258 4.177648
0.960.863851 5.112078
0.970.889596 6.659241
0.980.918411 9.731437
0.9900.9522311  18.874343
0.9950.9725934  37.055102
0.9970.9819976  61.229747
0.9980.9871752  91.409296
0.9990.9928860181.858807
1.0001.0000000  INFINITY
  1/30.0233204 1.847590
  2/30.388146 0.953613
0.567874750.25248543 0.888440006391162 =lowest value of g(tmax)
Figure 2a

Curves of tmax and 1/g(tmax) plotted as a function of β.

Figure 2b

Curve of 1/g(tmax) plotted as a function of tmax.

Figure 3

Second through fourth degree polynomial approximations to 1/g(tmax) as a function of β compared to the more accurately calculated value of this quantity.

  2 in total

1.  Phenomenological theory of structural relaxation based on a thermorheologically complex relaxation time distribution.

Authors:  L Andreozzi; M Faetti; M Salmerón Sanchez; J L Gómez Ribelles
Journal:  Eur Phys J E Soft Matter       Date:  2008-09       Impact factor: 1.890

2.  Generalized theory of the kinetics of tracers in biological systems.

Authors:  M Gitterman; G H Weiss
Journal:  Bull Math Biol       Date:  1994-03       Impact factor: 1.758

  2 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.