Application of Posts formula to optical pulse propagation in dispersive media
In this paper we consider the propagation of optical pulses in dielectric media with nontrivial dispersion relations. In particular, we implement Posts Laplace transform formula to invert in time the FourierLaplace space coefficients which arise from the joint space solution of the optical dispersive wave equation. Due to the inefficiency of a direct application of this formula, we have considered and present here two more efficient implementations. In the first, the GaverPost method, we utilize the well known Gaver functionals but store intermediate calculations to improve efficiency. The second, the BellPost method, involves an analytic reformulation of Posts formula such that knowledge of the dispersion relation and its derivatives are sufficient to invert the coefficients from Laplace space to time. Unlike other approaches to the dispersive wave equation which utilize a DebyeLorentzian assumption for the dispersion relation, our algorithm is applicable to general Maxwell-Hopkinson dielectrics. Moreover, we formulate the approach in terms of the FourierLaplace coefficients which are characteristic of a dispersive medium but are independent of the initial pulse profile. They thus can be precomputed and utilized when necessary in a real-time system. Finally, we present an illustration of the method applied to optical pulse propagation in a free space and in two materials with Cole-type dispersion relations, room temperature ionic liquid (RTIL) hexafluorophosphate and brain white matter. From an analysis of these examples, we find that both methods perform better than a standard Post-formula implementation. The BellPost method is the more robust of the two, while the GaverPost is more efficient at high precision and Post formula approximation orders.
posts formula optical dispersion numerical laplace transform inversion