Paraxial diffraction modeling based on the Fourier transform has seen widespread implementation for simulating the response of a diffraction-limited optical system. For systems where the paraxial ...assumption is not sufficient, a class of algorithms has been developed that employs hybrid propagation physics to compute the propagation of an elementary beamlet along geometric ray paths. These “beamlet decomposition” algorithms include the well-known Gaussian beamlet decomposition (GBD) algorithm, of which several variants have been created. To increase the computational efficiency of the GBD algorithm, we derive an alternative expression of the technique that utilizes the analytical propagation of beamlets to tilted planes. We then use this accelerated algorithm to conduct a parameter-space search to find the optimal combination of free parameters in GBD to construct the analytical Airy function. The experiment is conducted on a consumer-grade CPU, and a high-performance GPU, where the new algorithm is 34 times faster than the previously published algorithm on CPUs, and 67,513 times faster on GPUs.
Context.
Next-generation large segmented mirror telescopes are expected to perform direct imaging and characterization of Earth-like rocky planets, which requires contrast limits of 10
−7
to 10
−8
at ...wavelengths from
I
to
J
band. One critical aspect affecting the raw on-sky contrast are polarization aberrations (i.e., polarization-dependent phase and amplitude patterns in the pupil) arising from the reflection from the telescope’s mirror surfaces and instrument optics. These polarization aberrations induce false signals for polarimetry that can be calibrated to a certain degree, but they can also fundamentally limit the achievable contrast of coronagraphic systems.
Aims.
We simulate the polarization aberrations and estimate their effect on the achievable contrast for three next-generation ground-based large segmented mirror telescopes.
Methods.
We performed ray-tracing in Zemax
®
and computed the polarization aberrations and Jones pupil maps using the polarization ray-tracing algorithm. The impact of these aberrations on the contrast is estimated by propagating the Jones pupil maps through a set of idealized coronagraphs using hcipy, a physical optics-based simulation framework.
Results.
The optical modeling of the giant segmented mirror telescopes (GSMTs) shows that polarization aberrations create significant leakage through a coronagraphic system. The dominant aberration is retardance defocus, which originates from the steep angles on the primary and secondary mirrors. The retardance defocus limits the contrast to 10
−5
to 10
−4
at 1
λ
/
D
at visible wavelengths, and 10
−5
to 10
−6
at infrared wavelengths. The simulations also show that the coating plays a major role in determining the strength of the aberrations.
Conclusions.
Polarization aberrations will need to be considered during the design of high-contrast imaging instruments for the next generation of extremely large telescopes. This can be achieved either through compensation optics, robust coronagraphs, specialized coatings, calibration, and data analysis approaches, or by incorporating polarimetry with high-contrast imaging to measure these effects.
In the pursuit of directly imaging exoplanets, the high-contrast imaging community has developed a multitude of tools to simulate the performance of coronagraphs on segmented-aperture telescopes. As ...the scale of the telescope increases and science cases move toward shorter wavelengths, the required physical optics propagation to optimize high-contrast imaging instruments becomes computationally prohibitive. Gaussian Beamlet Decomposition (GBD) is an alternative method of physical optics propagation that decomposes an arbitrary wavefront into paraxial rays. These rays can be propagated expeditiously using ABCD matrices, and converted into their corresponding Gaussian beamlets to accurately model physical optics phenomena without the need of diffraction integrals. The GBD technique has seen recent development and implementation in commercial software (e.g. FRED, CODE V, ASAP) but appears to lack an open-source platform. We present a new GBD tool developed in Python to model physical optics phenomena, with the goal of alleviating the computational burden for modeling complex apertures, many-element systems, and introducing the capacity to model misalignment errors. This study demonstrates the synergy of the geometrical and physical regimes of optics utilized by the GBD technique, and is motivated by the need for advancing open-source physical optics propagators for segmented-aperture telescope coronagraph design and analysis. This work illustrates GBD with Poisson's spot calculations and show significant runtime advantage of GBD over Fresnel propagators for many-element systems.
For diffraction-limited optical systems an accurate physical optics model is necessary to properly evaluate instrument performance. Astronomical observatories outfitted with coronagraphs for direct ...exoplanet imaging require physical optics models to simulate the effects of misalignment and diffraction. Accurate knowledge of the observatory's PSF is integral for the design of high-contrast imaging instruments and simulation of astrophysical observations. The state of the art is to model the misalignment, ray aberration, and diffraction across multiple software packages, which complicates the design process. Gaussian Beamlet Decomposition (GBD) is a ray-based method of diffraction calculation that has been widely implemented in commercial optical design software. By performing the coherent calculation with data from the ray model of the observatory, the ray aberration errors can be fed directly into the physical optics model of the coronagraph, enabling a more integrated model of the observatory. We develop a formal algorithm for the transfer-matrix method of GBD, and evaluate it against analytical results and a traditional physical optics model to assess the suitability of GBD for high-contrast imaging simulations. Our GBD simulations of the observatory PSF, when compared to the analytical Airy function, have a sum-normalized RMS difference of ~10^-6. These fields are then propagated through a Fraunhofer model of a exoplanet imaging coronagraph where the mean residual numerical contrast is 4x10^-11, with a maximum near the inner working angle at 5x10^-9. These results show considerable promise for the future development of GBD as a viable propagation technique in high-contrast imaging. We developed this algorithm in an open-source software package and outlined a path for its continued development to increase the fidelity and flexibility of diffraction simulations using GBD.
High-contrast imaging in the next decade aims to image exoplanets at smaller angular separations and deeper contrasts than ever before. A problem that has recently garnered attention for telescopes ...equipped with high-contrast coronagraphs is polarization aberration arising from the optics. These aberrations manifest as low-order aberrations of different magnitudes for orthogonal polarization states and spread light into the dark hole of the coronagraph that cannot be fully corrected. The origin of polarization aberrations has been modeled at the telescope level. However, we don't fully understand how polarization aberrations arise at the instrument level. To directly measure this effect, we construct a dual-rotating-retarder polarimeter around the SCoOB high-contrast imaging testbed to measure its Mueller matrix. With this matrix, we directly characterize the diattenuation, retardance, and depolarization of the instrument as a function of position in the exit pupil. We measure the polarization aberrations in the Lyot plane to understand how polarization couples into high-contrast imaging residuals.