2024ApJ...961..112A

Extreme scattering events (ESEs) are observed as dramatic (>50%) drops in flux density that occur over an extended period of weeks to months. Discrete plasma lensing structures are theorized to scatter the radio waves produced by distant sources such as pulsars, causing the signature decrease in flux density and characteristic caustic spikes in ESE light curves. While plasma lens models in the extant literature have reproduced key features of ESE light curves, they have all faced the problem of being highly overdense and overpressured relative to the surrounding interstellar medium by orders of magnitude. We model ESEs by numerically ray tracing through analytic, volumetric plasma lens models by solving the eikonal equation. Delaunay triangulation connecting the rays approximates the wave front, generating a mapping from the observer plane to the source plane to account for multiple imaging. This eikonal method of ray tracing is tested against known analytic solutions and is then applied to a three-dimensional Gaussian-distributed electron volume density lens and a filament model inspired by Grafton et al. We find convergence of our numerical results with established analytic solutions, validating our numerical method, and reproduce ESE-like light curves. Our numerical ray-tracing method lends itself well to exploring the lensing effects of volumetric turbulence as well as sheet-like lenses, which is currently in progress.