Wave-optical study of the Einstein cross formed by a quadrupole gravitational lens

TL;DR

This paper develops an analytical approach using a quartic equation to model the electromagnetic field in gravitational lensing, specifically for quadrupole lenses, improving computational efficiency and accuracy in simulating Einstein crosses.

Contribution

It introduces a new diffraction integral and an analytical quartic solution for modeling the EM field in quadrupole gravitational lensing, enhancing simulation speed and accuracy.

Findings

Analytical quartic solution approximates EM field well across the image plane.

Solution captures the astroid caustic structure of quadrupole lenses.

Significant computational speed-up in simulating gravitational lensing phenomena.

Abstract

We study imaging of point sources with a quadrupole gravitational lens while focusing on the formation and evolution of the Einstein cross formed on the image sensor of an imaging telescope. We use a new type of a diffraction integral that we developed to study generic, opaque, weakly aspherical gravitational lenses. To evaluate this integral, we use the method of stationary phase that yields a quartic equation with respect to a Cartesian projection of the observer's position vector with respect to the vector of the impact parameter. The resulting quartic equation can be solved analytically using the method first published by Cardano in 1545. We find that the resulting solution provides a good approximation of the electromagnetic (EM) field almost everywhere in the image plane, yielding the well-known astroid caustic of the quadrupole lens. The sole exception is the immediate vicinity of the caustic boundary, where a numerical treatment of the diffraction integral yields better results. We also convolve the quartic solution for the EM field on the image plane with the point-spread function of a thin lens imaging telescope. By doing so, we are able to explore the direct relationship between the algebraic properties of the quartic solution for the EM field, the geometry of the astroid caustic, and the geometry and shape of the resulting Einstein-cross that appear on the image plane of the thin lens telescope. The new quartic solution leads to significant improvements in numerical modeling as evaluation of this solution is computationally far less expensive than a direct numerical treatment of the new diffraction integral. In the case of the solar gravitational lens (SGL), the new results drastically improve the speed of numerical simulations related to sensitivity analysis performed in the context of high-resolution imaging of exoplanets.