Numerical inversion of the Laplace transform
Abstract
A generalization of Doetsch's formula [Math. Z. 42, 263 (1937)] is derived to develop a stable numerical inversion of the one-sided Laplace transform Ĉ(β). The necessary input is only the values of Ĉ(β) on the positive real axis. The method is applicable provided that the functions Ĉ(β) belong to the function space Lα2 defined by the condition that G(x)=exαĈ(ex), α>0 has to be square integrable. The inversion algorithm consists of two sequential Fourier transforms where the second Fourier integration requires a cutoff, whose magnitude depends on the accuracy of the data. For high accuracy data, the cutoff tends to infinity and the inversion is very accurate. The presence of noise in the signal causes a lowering of the cutoff and a lowering of the accuracy of the inverted data. The optimal cutoff value is shown to be one which leads to an inversion which remains consistent with the original data and its noise level. The method is demonstrated for some model problems: a harmonic partition function, resonant transmission through a barrier, noisy correlation functions, and noisy Monte Carlo generated data for tunneling coefficients obtained via the recently introduced quantum transition state theory (QTST).
- Publication:
-
Journal of Chemical Physics
- Pub Date:
- June 1999
- DOI:
- 10.1063/1.479059
- arXiv:
- arXiv:physics/9807051
- Bibcode:
- 1999JChPh.11011176H
- Keywords:
-
- 02.30.Qy;
- 02.30.Cj;
- 02.60.Jh;
- 02.30.Nw;
- 02.50.Ng;
- 02.70.Lq;
- Measure and integration;
- Numerical differentiation and integration;
- Fourier analysis;
- Distribution theory and Monte Carlo studies;
- Physics - Data Analysis;
- Statistics and Probability;
- Physics - Computational Physics
- E-Print:
- 13 figures