A Critical Evaluation and Modification of the Padé-Laplace Method for Deconvolution of Viscoelastic Spectra

This paper shows that using the Padé-Laplace (PL) method for deconvolution of multi-exponential functions (stress relaxation of polymers) can produce ill-conditioned systems of equations. Analysis of different sets of generated data points from known multi-exponential functions indicates that by increasing the level of Padé approximants, the condition number of a matrix whose entries are coefficients of a Taylor series in the Laplace space grows rapidly. When higher levels of Padé approximants need to be computed to achieve stable modes for separation of exponentials, the problem of generating matrices with large condition numbers becomes more pronounced. The analysis in this paper discusses the origin of ill-posedness of the PL method and it was shown that ill-posedness may be regularized by reconstructing the system of equations and using singular value decomposition (SVD) for computation of the Padé table. Moreover, it is shown that after regularization, the PL method can deconvolute the exponential decays even when the input parameter of the method is out of its optimal range.