Acquiring the Lefschetz thimbles: efficient evaluation of the diffraction integral for lensing in wave optics

Xun Shi

DOI: https://doi.org/10.1093/mnras/stae2127

2024-09-19

Abstract:Evaluating the Kirchhoff-Fresnel diffraction integral is essential in studying wave effects in astrophysical lensing, but is often intractable because of the highly oscillatory integrated. A recent breakthrough was made by exploiting the Picard-Lefschetz theory: the integral can be performed along the `Lefschetz thimbles' in the complex domain where the integrand is not oscillatory but rapidly converging. The application of this method, however, has been limited by both the unfamiliar concepts involved and the low numerical efficiency of the method used to find the Lefschetz thimbles. In this paper, we give simple examples of the Lefschetz thimbles and define the `flow lines' that facilitate the understanding of the concepts. Based on this, we propose new ways to obtain the Lefschetz thimbles with high numerical efficiency, which provide an effective tool for studying wave effects in astrophysical lensing.

Instrumentation and Methods for Astrophysics,Cosmology and Nongalactic Astrophysics

What problem does this paper attempt to address?

The paper attempts to address the problem of efficiently evaluating the Kirchhoff-Fresnel diffraction integral in astrophysical lensing systems. The Kirchhoff-Fresnel diffraction integral is crucial for studying wave effects in astrophysical lenses, but it is difficult to compute due to its highly oscillatory integrand. Although Picard-Lefschetz theory offers a method to integrate in the complex domain along the "Lefschetz thimbles" (which are non-oscillatory and rapidly convergent), the application of this method is limited by conceptual complexity and numerical inefficiency.

Specifically, the paper addresses the problem through the following three aspects:
1. **Conceptual Clarification**: Illustrates the concept of Lefschetz thimbles with simple examples and defines "flow lines" to aid in understanding these concepts.
2. **Proposing New Methods**: Introduces new numerical methods to efficiently obtain Lefschetz thimbles.
3. **Numerical Optimization**: Identifies the shortcomings of traditional "flow methods" and proposes improved flow methods and a method based on constant phase contours, the latter being particularly suitable for 1-dimensional integrals and extendable to 2-dimensional integrals.

Through these methods, the paper aims to overcome the limitations of existing methods, thereby providing effective tools for studying wave effects in a wide range of lens configurations.