wolensing: A Python package for computing the amplification factor for gravitational waves with wave-optics effects

Simon M.C. Yeung,Mark H.Y. Cheung,Miguel Zumalacarregui,Otto A. Hannuksela
2024-10-17
Abstract:The wolensing Python package offers a solution for gravitational wave lensing computations within the full wave-optics regime. This tool is primarily designed to calculate the gravitational lensing amplification factor including diffractive effects, an essential component for generating accurate lensed gravitational wave waveforms. These waveforms are integral to astrophysical and cosmological studies related to gravitational-wave lensing. Integrating with lensingGW (Pagano, Hannuksela, and Li 2020), wolensing provides solutions for image positions in the high-frequency regime where wave and geometrical optics converge. This functionality allows the amplification factor to be applicable across a wider frequency range. Another key feature of wolensing is its ability to plot time delay contours on the lens plane, offering researchers a visual tool to better understand the relationship between the lens system and the amplification factor. wolensing is compatible with various lens models in lenstronomy (Birrer et al. 2021). There are also built-in lens models including point mass, singular isothermal sphere (SIS),and nonsingular isothermal ellipsoid (NIE) with jax (Bradbury et al. 2018) supporting GPU computation. Users can accommodate different lens models in the code with jax. wolensing is available as an open-source package on PyPI and can be installed via pip.
Instrumentation and Methods for Astrophysics
What problem does this paper attempt to address?
This paper introduces a Python package named **wolensing**, which aims to solve the problem of calculating the wave - optical effects in the gravitational - wave lensing effect. Specifically, it attempts to solve the following problems: 1. **Limitations of the traditional geometric - optical approximation method**: - Traditional gravitational - wave lensing research mainly focuses on strong lensing effects and uses the geometric - optical approximation method. This method assumes that light travels in a straight line and can predict the time delay and magnification of different images, but ignores the change in frequency evolution. - For lensing celestial bodies with a mass of around 100 solar masses or less (such as stars and the remnants of low - mass compact objects), their gravitational radius (Schwarzschild radius) is comparable to the gravitational - wave wavelength in the sensitive frequency band (10 - 10^3 Hz) of detectors such as LIGO/Virgo/Kagra. In this case, the wave - optical effects (especially the diffraction effect) become significant, and it is necessary to shift from geometric optics to wave optics for accurate modeling. 2. **Complexity introduced by wave - optical effects**: - In the wave - optical regime, the frequency evolution of gravitational waves is affected by the amplification factor, which is calculated through a complex diffraction integral. This leads to a significant increase in computational complexity and cost. - The wolensing package addresses this challenge by providing an efficient method for calculating the amplification factor, especially when dealing with general lensing systems. 3. **Cross - frequency - band applicability**: - wolensing is not only applicable to the high - frequency regime (where geometric optics and wave optics converge), but can also be applied in a wider frequency range, including the sensitive frequency band (mHz - Hz) of LISA (Laser Interferometer Space Antenna), covering lensing celestial bodies with a mass from 10^5 to 10^8 solar masses (such as dark - matter sub - halos and larger compact structures). 4. **Handling of multi - lens systems**: - wolensing is also able to handle complex systems containing multiple lenses, not limited to the case of a single lens. In summary, the wolensing package aims to provide an efficient and accurate method for calculating the amplification factor in the gravitational - wave lensing effect, especially in the wave - optical regime, to support more extensive and accurate astrophysical and cosmological research.