Metaplectic Geometrical Optics

Ray optics is an intuitive and computationally efficient model for wave propagation through nonuniform media. However, the underlying geometrical-optics (GO) approximation of ray optics breaks down at caustics, erroneously predicting the wave intensity to be infinite and thereby limiting the predictive capabilities of GO-based codes. Full-wave modeling can be used instead, but the added computational cost brings its own set of tradeoffs. Developing more efficient caustic remedies has therefore been an active area of research for the past few decades. In this thesis, I present a new ray-based approach called 'metaplectic geometrical optics' (MGO) that can be applied to any linear wave equation. Instead of evolving waves in the usual x (coordinate) or k (spectral) representation, MGO uses a mixed X=Ax+B representation. By continuously adjusting the coefficients A and B along the rays via sequenced metaplectic transforms (MTs) of the wavefield, corresponding to symplectic transformations of the ray phase space, one can ensure that GO remains valid in the X coordinates without caustic singularities. The caustic-free result is then mapped back onto the original x space using metaplectic transforms, as demonstrated on a number of examples. Besides outlining the basic theory, this thesis also presents specialized fast algorithms for MGO. These algorithms focus on the MT, which is a unitary integral mapping that can be considered a generalization of the Fourier transform. First, a discrete representation of the MT is developed that can be computed in linear time when evaluated in the near-identity limit; finite MTs can then be implemented by iterating $K\gg 1$ near-identity MTs. Second, an algorithm based on Gauss--Freud quadrature is developed for efficiently computing finite MTs along their steepest-descent curves, which may be useful in catastrophe-optics applications beyond MGO.