Next Article in Journal
Modeling for Generating Femtosecond Pulses in an Er-Doped Fiber Using Externally Controlled Spectral Broadening and Compression Mechanisms
Next Article in Special Issue
On-Demand Phase Control of a 7-Fiber Amplifiers Array with Neural Network and Quasi-Reinforcement Learning
Previous Article in Journal
Single-Pixel Moving Object Classification with Differential Measuring in Transform Domain and Deep Learning
Previous Article in Special Issue
Design, Fabrication and Characterization of an Adaptive Retroreflector (AR)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Detection of Wave Aberrations Based on the Multichannel Filter

by
Pavel A. Khorin
1,
Alexey P. Porfirev
2 and
Svetlana N. Khonina
1,2,*
1
Samara National Research University, 443086 Samara, Russia
2
Image Processing Systems Institute of RAS—Branch of the FSRC “Crystallography and Photonics” RAS, 443001 Samara, Russia
*
Author to whom correspondence should be addressed.
Photonics 2022, 9(3), 204; https://doi.org/10.3390/photonics9030204
Submission received: 30 January 2022 / Revised: 8 March 2022 / Accepted: 19 March 2022 / Published: 21 March 2022
(This article belongs to the Special Issue Various Applications of Methods and Elements of Adaptive Optics)

Abstract

:
An adaptive method for determining the type and magnitude of aberration in a wide range is proposed on the basis of an optical processing of the analyzed wavefront using a multichannel filter matched to the adjustable Zernike phase functions. The approach is based on an adaptive (or step-by-step) compensation of wavefront aberrations based on a dynamically tunable multichannel filter implemented on a spatial light modulator. For adaptive filter adjustment, a set of criteria is proposed that takes into account not only the magnitude of the correlation peak, but also the maximum intensity, compactness, and orientation of the distribution in each diffraction order. The experimental results have shown the efficiency of the proposed approach for detecting wavefront aberrations in a wide range (from 0.1λ to λ).

1. Introduction

The problem of measuring and correcting wavefront aberrations is often encountered in optics, for example, in the design of ground-based telescopes, in optical communication systems, in industrial laser technology, and in medicine [1,2,3,4,5,6,7,8,9,10,11,12]. Usually, the measurement of wavefront distortions is performed in order to compensate them, in particular, with adaptive or active optics [13,14,15,16,17,18]. The major causes of wavefront aberrations are turbulence of the atmosphere, imperfect shapes of the optical elements of the system, errors in the alignment of the system, etc.
It is known that weak wavefront aberrations (level ≤ 0.4λ) are well detected using spatial filters matched to the basis of Zernike functions [19,20,21,22,23,24,25,26,27] including multichannel diffractive optical elements (DOEs) [21,25,27]. However, with an increase in aberration level, the linear approximation of the wavefront by Zernike functions becomes unacceptable [27]. This is explained by the fact that the contribution of the second and subsequent nonlinear terms of the wavefront expansion to the Taylor series becomes more significant, which leads to the detection of false aberrations.
With high aberrations (level > 0.4λ), when a significant blurring of the focal spot occurs, it makes sense to use methods focused on analyzing the intensity distribution pattern formed by an aberrated optical system in one or several planes. To determine the wavefront in this case, iterative [28,29,30,31,32] and optimization algorithms [10,33] are used, including those with the use of neural networks [34,35,36,37,38,39]. In turn, these approaches demonstrate significant errors for small aberrations, when the point spread function (PSF) is close to the Airy picture of an ideal system [27].
Thus, different methods work at different levels of aberrations, and in order to apply them, it is desirable to determine this level (or magnitude). One of the solutions is the use of additional optical and digital processing, for example, based on a dynamically tunable spatial light modulator (SLM). Previously, we studied the stability of the wavefront expansion coefficients on the basis of Zernike polynomials during the field propagation in free space [40], the application limits of spatial filters matched with the basis of Zernike functions [41,42,43], and the possibility of scaling aberrations levels for testing optical systems [44].
In this paper, to determine the type and magnitude (or level) of aberrations in the investigated wavefront (WF), the application of an adaptive method based on the use of a multichannel filter matched with adjustable Zernike phase functions (Zernike polynomials correspond to the phase of the considered functions) is proposed. In this case, instead of the optical expansion of the field on the basis of the Zernike functions which was realized earlier [21,25,27], we actually perform multichannel aberration compensation based on the Zernike polynomials. The novelty of our study lies in the combination of an adaptive approach and matched filtering based on a multichannel diffractive optical element. Note that phase compensation can occur in each channel in accordance with different types of aberrations, as well as for the same aberration type but with different magnitude (or level). The method is based on a step-by-step compensation of wavefront aberrations based on a dynamically tunable multichannel filter implemented on a spatial light modulator. A set of criteria for adaptive filter tuning is proposed, taking into account not only the correlation peak presence, but also the maximum intensity, compactness, and orientation of the PSF distribution in each diffraction order. The experimental results have shown the efficiency of the proposed approach for detecting wavefront aberrations in a wide range (from 0.1λ to λ).

2. Materials and Methods

2.1. Theoretical Foundations

Wavefront aberrations are usually described in terms of Zernike functions Z n m r , φ [45] in the following way:
W r , φ   = exp i κ α n = 0 N max m = 0 n C n m Z n m r , φ ,
where κ = 2π/λ is the wave number, λ is the radiation wavelength, α is the parameter corresponding to the level of aberration in wavelengths (in this paper it takes values from zero to wavelength λ), and coefficients Cnm are normalized.
The PSF for an aberrated wavefront can be calculated using the Fourier transform:
F ( u , v ) = W ( x , y )   = W ( x , y ) exp i 2 π λ f u x + v y   d x d y ,
where f is the focal length of the lens.
The optical wavefront analyzer considered in this work is based on the addition of a lens with a diffractive multichannel filter [21,25,27] of the following form:
τ ( x , y )   = p = 0 P q = 0 Q Ψ p q * x , y exp i a p q x + b p q y ,
where Ψpq (x, y) are the functions matched to diffraction orders (DOs) with indices (p, q), the position of which in the focal plane is determined by the spatial frequencies apq and bpq.
As a rule, the Zernike polynomials are used as Ψpq (x, y) [7,21,25,27] to analyze a wavefront. However, correct detection in this case is possible only at low levels of aberrations. In [27], it was shown that at an aberration level of more than 0.4λ, the use of such a filter leads to false detection.
In this paper, the use of a filter matched with the Zernike phase functions (Zernike polynomials correspond to the phase of the considered functions) is proposed:
Ψ p q x , y   = exp i κ α 0 Z p q x , y ,
where Z p q x , y are the standard Zernike functions of double indices (p, q) [45]. The correspondence of Zernike functions to various types of aberrations is given in Table 1.

2.2. Scheme for Adaptive Aberration Compensation

We consider in detail several wave aberrations of different magnitudes and show the scheme for the adaptive method to compensate these aberrations.
We divide the aberrations into even type (n is even) and odd type (n is odd) and consider them separately. Among the even aberrations, one can distinguish axis-symmetric ones (m is equal to 0). Even aberrations such as (2, 0), (4, 0), …, (n, 0) have the wavefront phase and corresponding PSF axis-symmetric distribution regardless of the level of aberration α (see Table 2).
The adaptive method is described as the aberration compensation by complex conjugate functions (4) with different levels of α0. In other words, we calculate the compensated wavefront of the form
Δ W ( x , y )   = W x , y Ψ p q * x , y   = exp i κ α Z n m x , y exp i κ α 0 Z p q x , y
and the corresponding PSF of the form
Δ F ( u , v )   = Δ W ( x , y )   = Δ W ( x , y ) exp i 2 π λ f u x + v y d x d y ,
where W(x, y) is the analyzed wavefront with aberration (n, m) and level α, Ψ p q * x , y is the complex conjugate function (4) with a variable level α0, ΔW(x, y) is the resulting wavefront, and ΔF(x, y) is the corresponding PSF.
With the successful compensation of the specified wave aberration, the phase of the resulting wavefront should be equal to a constant, and the PSF should correspond to the Airy picture (PSF for the plane WF).
The results collected in Figure 1 show the compensating process for the spherical aberration (n, m) = (2, 0) with a level α = 0.5λ by the set of Ψ p q * x , y = Ψ 20 * x , y with different levels of α0 = {0.1λ; 0.25λ; 0.5λ; 0.75λ}. Thus, we construct the resulting wavefront defined by Equation (5) and the corresponding PSF defined by Equation (6).
Figure 1 shows in detail the effect of the compensating wavefront W(x, y) with aberration type (n, m) = (2, 0) and level α = 0.5λ by Ψ p q * x , y = Ψ 20 * x , y with different levels of α0. The compensated phase (equal to 0) is observed at the intersection of column 4 and row 3, and the corresponding PSF (the Airy picture D ( u , v ) ) is observed at the intersection of column 4 and row 4:
Δ W ( x , y )   = W x , y Ψ p q * x , y = exp i 2 π 0.5 Z 20 x , y exp i 2 π 0.5 Z 20 x , y = 1 ,
Δ F ( u , v )   = Δ W ( x , y ) = exp i 2 π λ f u x + v y d x d y = D ( u , v ) .
In addition, it can be seen from Figure 1 that the PSFs observed in row 4, columns 3 and 5, have the same intensity. This is because the phase of the resulting wavefront for these two cases is the same up to sign:
Δ W ( x , y )   = exp i 2 π 0.5 Z 20 x , y exp i 2 π 0.25 Z 20 x , y = exp i 2 π 0.25 Z 20 ,
Δ W ( x , y )   = exp i 2 π 0.5 Z 20 x , y exp i 2 π 0.75 Z 20 x , y = exp i 2 π 0.25 Z 20 .
A change in the sign of the aberration coefficient to the opposite causes a symmetric transformation of the PSF image about the horizontal axis for negative m, since −sin(mφ) = sin(mφ + π), and about the vertical axis for positive m, since −cos(mφ) = cos(mφ + π). For even aberrations, this property is not so noticeable due to their symmetry. Therefore, PSFs for even wave aberrations that differ only in the sign of the weighting coefficient have the same intensity distribution in the focal plane.
Even aberrations with m ≠ 0 have phase images with m symmetry axes and PSFs with 2m symmetry axes. It should be noted that the OX and OY axes are always the PSF symmetry axes (see Table 3).
Analogical results of the compensating process for the astigmatic aberration (n, m) = (2, 2) with a level α = 0.5λ by the set of Ψ p q * x , y   = Ψ 22 * x , y with different levels of α0 = {0.1λ; 0.25λ; 0.5λ; 0.75λ} are shown in Figure 2. The situation is similar to that of Figure 1, and cases at α0 = {0.25λ; 0.75λ} correspond to identical PSFs.
Odd-type aberrations such as (n, m) = (3, 1) and (n, m) = (3, 3) have images of the wavefront phase with m symmetry axes and the corresponding PSFs with 2m symmetry axes (see Table 4).
The results of the compensating process for the coma-type aberration (n, m) = (3, 1) with a level α = 0.5λ by the set of Ψ p q * x , y = Ψ 31 * x , y with different levels of α0 = {0.1λ; 0.25λ; 0.5λ; 0.75λ} are shown in Figure 3. The situation is different from even aberrations.
It can be seen from Figure 3 that in row 4, columns 3 and 5, PSFs are rotated by 180 degrees. This is because the phase of the resulting wavefront for these two cases is the same up to sign:
Δ W ( x , y )   = exp i 2 π 0.5 Z 31 x , y exp i 2 π 0.25 Z 31 x , y   = exp i 2 π 0.25 Z 31 ,
Δ W ( x , y )   = exp i 2 π 0.5 Z 31 x , y exp i 2 π 0.75 Z 31 x , y   = exp i 2 π 0.25 Z 31 .
Changing the sign of the aberration coefficient to the opposite one causes a symmetric transformation of the PSF image about the horizontal axis. As described earlier for positive m, such a symmetric transformation must be performed with respect to the vertical axis, since −cos(mφ) = cos(mφ + π). For odd-type aberrations, this property is evidentially observed in the PSF intensity distribution.
Thus, the introduced empirical criterion for odd-type aberrations is based on the described fundamental properties of the PSF of wave aberrations.
The following combined criteria for correct detection can be used:
  • Formation of a pronounced/evident correlation peak in the focal plane of the filter in the DO with indices (p, q).
  • Change in the orientation of the PSF intensity distribution in the adjacent DOs of the “level and type” (LT) filter within one row (for odd-symmetry aberrations).
  • Presence of the same PSF intensity distributions in the adjacent DOs of the LT-filter within one row (for even-symmetry aberrations).
Based on the analysis of the intensity distribution in the focal plane, taking into account the proposed criteria, it is possible to determine the range of the analyzed aberration level.
Next, we propose using a multichannel filter to implement the processing of the analyzed wavefront simultaneously by several types of aberrations taken into account.

2.3. Multichannel Filter Employment

To determine the type and magnitude of aberration in the analyzed wavefront, we propose combining the adaptive method with the use of a multichannel filter matched to Zernike’s adjustable phase functions. Thus, the principle of operation of our method lies in the parallel application of the adaptive approach by a multichannel diffractive element employment. In this case, a wide range of detectable aberration values and resistance to vibrations are provided.
To realize this approach, the filter should be matched with a set of functions defined by Equation (4) with different levels of αk. We can design a filter tuned to one specific wave aberration with a value ranging from α1 to αKmax.
In order to detect another type of wave aberrations, it is necessary to divide each channel into PmaxQmax additional channels. So, different wave aberrations with different αk should be encoded in different diffraction orders. The transmission function of such an LT-filter will be as follows:
τ L T ( x , y )   = p = 0 P max q = 0 Q max k = 1 K max exp i κ α k Z p q x , y exp i a k p q x + b k p q y .
It should be noted that the accuracy of the detected (and compensated) aberration value directly depends on the number of channels PmaxQmaxKmax; however, the number of filter channels is limited by technological possibility. Since the filter defined by Equation (13) can be implemented using an SLM, the number of channels can be reduced taking into account a certain time delay for rebuilding the filter to different sets of channels. Moreover, in the presence of an additional adapting device, it is possible to organize an iterative process of step-by-step compensation of aberrations [30,31,32].
The formation of a correlation peak in one of the DOs can be considered as a criterion for successful compensation. Additional criteria for adaptive filter tuning are proposed in Section 3.
A special case of the LT-filter defined by Equation (13) is the “type” T-filter tuned for different types of aberrations with the same level α0:
τ T ( x , y )   = p = 0 P max q = 0 Q max exp i κ α 0 Z p q x , y exp i a p q x + b p q y .
If the aberration type is determined and we need to more accurately determine its level, then it makes sense to use another particular version of the filter defined by Equation (13), namely the “level” L-filter tuned to only one aberration type (p0, q0) with different levels of αk:
τ L ( x , y )   = k = 1 K max exp i κ α k Z p 0 q 0 x , y exp i a k x + b k y .
To analyze other aberration types, the filter can be quickly rebuilt using an SLM.
We have designed several variants of multichannel filters defined by Equations (13)–(15). Figure 4 demonstrates three examples of these filters: T-filter matched with different types of aberrations of the same level (Equation (14), Figure 4a), LT-filter matched with different types and levels of aberrations (Equation (13), Figure 4b), and L-filter matched with one aberration type of different levels (Equation (15), Figure 4c). Each of the options is suitable for a specific task. It makes sense to use a hybrid LT-filter for preliminary analysis and a filter consistent with different levels for refinement.
With the help of dynamically tunable filters implemented on an SLM, it becomes possible to specify the level of aberrations of the analyzed WF or to take the average value from the already detected range. Note that the accuracy of this range directly depends on the noise level of the recording device.

2.4. Optical Setup

The optical scheme used in the experiment is shown in Figure 5.
The experimental setup for detecting wavefront aberrations using a diffraction multichannel filter consists of the following elements: a solid-state laser (λ = 532 nm), a pinhole (hole size of 40 μm), spherical lenses (f1 = 350 mm, f2 = 300 mm, f3 = 200 mm, and f4 = 250 mm), the transmissive SLM (HOLOEYE LC 2012), the reflective SLM (HOLOEYE PLUTO VIS), circular apertures, mirrors, and the video camera ToupCam UCMOS08000KPB.
The output laser radiation from a solid-state laser (λ = 532 nm) was collimated using a system consisting of a pinhole (PH) with a hole diameter of 40 μm and a spherical lens (L1) (f1 = 350 mm). A circular diaphragm D1 was used to separate a central light spot from surrounding light and dark rings occurring during diffraction at the pinhole. Then, the laser beam expanded and reflected by mirror M1 passed through a HOLOEYE LC 2012 transmissive spatial light modulator SLM1 with a 1024 × 768 pixel resolution and a pixel size of 36 μm, which was used to form a wavefront with a required combination of aberrations. Lenses L2 (f2 = 300 mm) and L3 (f3 = 200 mm) and diaphragm D2 were used to spatially split the aberration-distorted laser beam formed by the first modulator and the unmodulated zero-order transmitted laser beam.
A mirror M2 was used to direct the formed laser beam to the display of the second modulator. A HOLOEYE PLUTO VIS reflective spatial light modulator SLM2 with a 1920 × 1080 pixel resolution and a pixel size of 8 μm was used to implement a phase mask of a diffractive multichannel filter, which served to decompose the studied light field in terms of the Zernike polynomial basis. The laser beam reflected from this SLM was directed to lens L4 (f4 = 250 mm), which focused it on the matrix of a ToupCam UCMOS08000KPB camera CAM with a 3264 × 2448 resolution and a pixel size of 1.67 μm.

3. Results

3.1. Detection of One Type of Wave Aberration

As a rule, when using matched filters, the appearance of a correlation peak is a criterion for detecting a matched signal [46,47,48,49,50]. Note that the presence of a correlation peak is determined by a nonzero intensity at the center of the DO. However, when analyzing phase distributions, this approach is correct only for small phase variations. In particular, it was shown in [27] that when the aberration level is more than 0.4λ, the use of filters matched with the basis of Zernike functions leads to the appearance of many wrong correlation peaks and false detection.
The adaptive approach considered in Section 2.2 is free from this limitation, but for its application, it is necessary to select suitable compensating functions in accordance with the developed criteria. To speed up the selection process, we use multichannel filters to analyze a wavefront simultaneously by several types of compensating functions.
In this section, we demonstrate the use of the multichannel filters described in Section 2.3 for test examples when aberrations of certain types with a given level are introduced into the wavefront.
Figure 6 and Table 5 show the results of a numerical and experimental analysis of the WF distorted by a coma-type aberration (n, m) = (3, 1) with a weight (level) α = 0.5λ, using various diffractive multichannel filters (the intensity patterns in the focal plane are shown). DOs with correlation peaks are framed by a solid line, and DOs with changed PSF orientation are framed by a dashed line. The selected (framed) orders correspond to the criteria formulated in Section 2.2. A good agreement between the numerical simulation and the optical experiment is seen.
The results for LT-filter show (see Figure 6) that the correlation peak located in the DO corresponding to (p, q) = (3, 1) with αk = 0.5λ and adjacent DOs contains a PSF with changed orientation (see Figure 6, positions are marked with frames). As follows from the given example, for aberrations with odd symmetry, the criterion of the correct detection can be not only the presence of the correlation peak, but also the orientation of PSF distribution in neighboring DOs. A similar result was obtained in [51,52] based on numerical simulation. In contrast to these papers, in this work, we realize experimental investigation of the proposed approach and apply it for a step-by-step compensation of a complex type of aberrations.
T-filters with different levels of α0 also demonstrate (see Table 5) a change of PSF orientation when passing through the level α0 = 0.5λ. A correlation peak is also observed in the DO corresponding to (p, q) = (3, 1) with a level equal to 0.5λ (see Table 5, third row).
The results of simulation and experiment with another test example, when the WF is distorted by an even type of aberration (p, q) = (2, 2) with α = 0.5λ, are presented in Figure 7. Results show that the DOs adjacent to the correlation peak have the same PSF intensity distribution. Thus, when detecting even-symmetry aberrations, in addition to the presence of a correlation peak, the criterion for correct detection can be the similarity of diffraction patterns in neighboring orders.
The results of a test example with a WF distorted by aberration (n, m) = (3, 3) with an increased weight α = 0.7λ are shown in Figure 8 and Table 6.
The increase in weight was considered to show the good performance of filters over a wide range of aberration values (if necessary, the filter can be easily reconfigured for a larger range). In the focal plane of the LT-filter (see Figure 8), the PSF changes orientation in the DOs marked with frames (which corresponds to an aberration of the (p, q) = (3, 3) aberration type); however, no evident correlation peak is observed between these DOs. The absence of such a correlation peak is explained by the inconsistency of the level of aberrations recorded in the filter with the level of the analyzed wave aberration. To refine the weight of the detected aberration, it is required to use an L-filter or T-filter with different levels of aberration (p, q) = (3, 3) (see Table 6).
It can be seen from Table 6 that the use of a T-filter with different levels of matched functions made it possible to refine the weight of the detected aberration (p, q) = (3, 3); thus, the correlation peak is observed at α0 = 0.7λ (see Table 6, second rows). In this case, for the levels α0 = 0.5λ and α0 = 0.9λ (see Table 6, first and third rows), the intensity distribution changes orientation, which is an additional sign of correct detection.

3.2. Detection of Superposition of Wave Aberrations (Complex Aberration)

In this section, we illustrate the action of the adaptive (step-by-step) method for detecting a superposition of wave aberrations, where at each l-th stage the wavefront W l ( x , y ) (at the first iteration it is equal to the investigated wavefront) is analyzed by the dynamically tunable SLM. In the focal plane of the filter, there are different PSF pictures. We recall that the aberration type (p, q) is associated with each string and aberration level with k-column ( α k are fitted in accordance with the combined criteria introduced above) of the diffractive multichannel filter. So, the investigated wavefront is corrected taking into account the detected wave aberration as follows:
W l + 1 ( x , y )   = W l ( x , y ) / W Δ ( x , y )   = W l ( x , y ) exp i κ α k Z p q x , y ,
where W Δ ( x , y )   = exp i κ α k Z p q x , y .
The criterion for correctly detected aberrations is the formation of a pronounced/evident correlation peak in the focal plane of the filter in the DO at one of the l-th steps. To detect the correlation peak formation, we use digital image processing and calculate mean square error (MSE) for W l ( x , y ) in the region of DO (the considered region size is S × S pixels):
M S E W l ( u , v ) , D ( u , v )   = 1 S 2 i = 1 S j = 1 S W l ( u i , v j ) D ( u i , v j ) 2 .
where D ( u , v ) is the Airy picture (PSF for plane/ideal WF).
To illustrate this approach, we consider the WF distorted by the superposition of two aberrations 0.5λ*(2, 2) + 0.5λ*(3, 3). At the first stage, we used an LT-filter (see Figure 9) and revealed that the DO changes orientation in the fourth row, which corresponds to type (p, q) = (3, 3) with α = 0.5λ, and a residual aberration of the type (n, m) = (2, 2) is observed in the corresponding DO. So, at the first iteration, the wavefront aberration (n, m) = (3, 3) with 0.5λ level was compensated by LT-filter. Thus, the problem was reduced by one dimension.
After compensating for one detected aberration, the wavefront has residual aberration (n, m) = (2, 2) with 0.5λ level. The result of the action of the same LT-filter on the obtained wavefront is presented in Figure 10. The LT-filter shows that the DOs are symmetric in the second row with respect to the correlation peak at (p, q) = (2, 2) with 0.5λ level.
In another test example (Figure 11), we considered WF distorted by the superposition of three aberrations 0.5λ*(2, 0) + 0.5λ*(2, 2) + 0.5λ*(3, 3). The LT-filter shows (see Figure 11) that the DO changes orientation in the fourth row, which corresponds to aberration (p, q) = (3, 3) with 0.5λ level. So, at the first iteration, we compensated for one type of aberration and reduced the problem by one dimension. At the second iteration, the residual aberrations are similar to the first example (see Figure 9), which was compensated in a similar way.
To refine the weight of the desired aberration, it is possible to iteratively subtract the detected aberrations with low weights until a pronounced correlation peak is observed in one of the DOs.
Thus, using the adaptive/iterative process and the introduced empirical criteria, the wavefront aberrations are compensated. The criterion for the iterative process completion is the formation of a pronounced correlation peak in one of the DOs of the focal plane of the filter, i.e., when the MSE defined by Equation (17) equals a minimal value.

4. Discussion

In order to show the possibility of using the proposed method for detecting wavefronts with complex aberrations, we investigate how much one wave aberration (n1, m1) with a level α1 affects the wave aberration (n2, m2) with a level α2. For an objective assessment, the analysis of the PSF corresponding to the superposition of aberrations is proposed.
Figure 12 shows the simulation results for the superposition of even- and odd-type aberrations. It is seen that if the ratio of the level of aberrations is more than 10 times, then the effect of an aberration with a lower weight coefficient on the resulting PSF is extremely weak. In addition, the ratio of the level of aberrations being more than 5 times also insignificantly affects the resulting wavefront (with an accuracy to the scale and overall shape of the PSF).
Thus, if the superposition contains aberrations with coefficients different in magnitude, then it makes sense to first compensate for aberration with a large level. This fact allows us to say that the analyzed wavefront with a pronounced single aberration can be successfully detected by the proposed filters defined by Equations (13)–(15) without using an additional iterative algorithm.

5. Conclusions

In this work, we experimentally tested a diffractive multichannel wavefront sensor matched to phase distributions in the form of Zernike functions.
Several variants of filters were designed: filters matched with different types of aberrations of the same level, filters matched with one type of aberration of different levels, and a hybrid filter. Each of the options is convenient for a specific task. It makes sense to apply a hybrid filter for preliminary analysis and then apply a filter matched to different levels for adjusting. When using a dynamic transparency, it is possible to quickly adjust the filters and compensate for the aberration step by step.
Criteria of the proposed adaptive algorithm for determining the aberrations of the analyzed wavefront were formulated: scattering region, PSF compactness, maximum intensity, and orientation of each diffraction order in the focal plane of the multichannel filter. An adaptive compensation of wavefront aberrations was performed using the proposed method and introduced empirical criteria. We show the possibility to reduce a problem dimension at each stage, so the wavefront aberration compensation is provided. We realized an experimental investigation of the proposed approach and applied it for a step-by-step compensation of a complex type of aberrations.
Finally, based on numerical modeling and the results of an optical experiment using a dynamic transparency, we show the efficiency of the proposed filters for detecting wavefront aberrations in a wide range (from 0.1λ to λ).
Considering the wide range of correctly detected aberration values using the proposed method, the field of application can include the measurement and correction of wavefront aberrations in devices such as ground-based telescopes and optical microscopes, in optical communication systems, in industrial laser technology, and in ophthalmology.

Author Contributions

Conceptualization: S.N.K. and P.A.K.; methodology, S.N.K. and P.A.K.; software, P.A.K.; validation, S.N.K., P.A.K. and A.P.P.; formal analysis, S.N.K. and P.A.K.; investigation, S.N.K., P.A.K. and A.P.P.; resources, S.N.K.; data curation, S.N.K., P.A.K. and A.P.P.; writing—original draft preparation, S.N.K. and P.A.K.; writing—review and editing, S.N.K. and P.A.K.; visualization, S.N.K., P.A.K. and A.P.P.; supervision, S.N.K.; project administration, S.N.K.; funding acquisition, S.N.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the Russian Foundation for Basic Research (grant No. 20-37-90129) for numerical calculations and analysis and by the Ministry of Science and Higher Education of the Russian Federation under the FSRC “Crystallography and Photonics” of the Russian Academy of Sciences (the state task No. 007-GZ/Ch3363/26) for experimental research.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Camacho, L.; Mico, V.; Zalevsky, Z.; Garcia, J. Quantitative phase microscopy using defocusing by means of a spatial light modulator. Opt. Express 2010, 18, 6755–6766. [Google Scholar] [CrossRef] [PubMed]
  2. Lombardo, M.; Lombardo, G. Wave aberration of human eyes and new descriptors of image optical quality and visual performance. J. Cataract Refract. Surg. 2010, 36, 313–320. [Google Scholar] [CrossRef] [PubMed]
  3. Zhao, Q.; Fan, H.; Hu, S.; Zhong, M.; Baida, L. Effect of optical aberration of telescopes to the laser radar. Proc. SPIE 2010, 7656, 76565Z. [Google Scholar]
  4. González-Núñez, H.; Prieto-Blanco, X.; De la Fuente, R. Pupil aberrations in Offner spectrometers. J. Opt. Soc. Am. A 2011, 29, 442–449. [Google Scholar] [CrossRef] [PubMed]
  5. Khonina, S.N.; Ustinov, A.V.; Pelevina, E.A. Analysis of wave aberration influence on reducing focal spot size in a high-aperture focusing system. J. Opt. 2011, 13, 095702. [Google Scholar] [CrossRef]
  6. Booth, M.; Andrade, D.; Burke, D.; Patton, B.; Zurauskas, M. Aberrations and adaptive optics in super-resolution microscopy. Microscopy 2015, 64, 251–261. [Google Scholar] [CrossRef] [Green Version]
  7. Khorin, P.A.; Khonina, S.N.; Karsakov, A.V.; Branchevskiy, S.L. Analysis of corneal aberration of the human eye. Comput. Opt. 2016, 40, 810–817. [Google Scholar] [CrossRef]
  8. Wilby, M.J.; Keller, C.U.; Haert, S.; Korkiakoski, V.; Snik, F.; Pietrow, A.G.M. Designing and testing the coronagraphic modal wavefront sensor: A fast non-common path error sensor for high-contrast imaging. Proc. SPIE 2016, 9909, 990921. [Google Scholar]
  9. Klebanov, I.M.; Karsakov, A.V.; Khonina, S.N.; Davydov, A.N.; Polyakov, K.A. Wavefront aberration compensation of space telescopes with telescope temperature field adjustment. Comput. Opt. 2017, 41, 30–36. [Google Scholar] [CrossRef] [Green Version]
  10. Rastorguev, A.A.; Kharitonov, S.I.; Kazanskiy, N.L. Modeling of arrangement tolerances for the optical elements in a spaceborne Offner imaging hyperspectrometer. Comput. Opt. 2018, 42, 424–431. [Google Scholar] [CrossRef]
  11. Martins, A.C.; Vohnsen, B. Measuring ocular aberrations sequentially using a digital micromirror device. Micromachines 2019, 10, 117. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Baum, O.I.; Omel’chenko, A.I.; Kasianenko, E.M.; Skidanov, R.V.; Kazanskiy, N.L.; Sobol’, E.N.; Bolshunov, A.V.; Avetisov, S.E.; Panchenko, V.Y. Control of laser-beam spatial distribution for correcting the shape and refraction of eye cornea. Quantum Electron. 2020, 50, 87–93. [Google Scholar] [CrossRef]
  13. Mu, Q.; Cao, Z.; Hu, L.; Li, D.; Xuan, L. Adaptive optics imaging system based on a high-resolution liquid crystal on silicon device. Opt. Express 2006, 14, 8013–8018. [Google Scholar] [CrossRef] [PubMed]
  14. Ellerbroek, B.L.; Vogel, C.R. Inverse problems in astronomical adaptive optics. Inverse Probl. 2009, 25, 063001. [Google Scholar] [CrossRef]
  15. Esposito, S.; Riccardi, A.; Pinna, E.; Puglisi, A.; Quirós-Pacheco, F.; Arcidiacono, C.; Xompero, M.; Briguglio, R.; Agapito, G.; Busoni, L.; et al. Large binocular telescope adaptive optics system: New achievements and perspectives in adaptive optics. Proc. SPIE 2011, 8149, 814902. [Google Scholar]
  16. Lukin, V.P. Adaptive optics in the formation of optical beams and images. Phys. Uspekhi 2014, 57, 556–592. [Google Scholar] [CrossRef]
  17. Ji, N. Adaptive optical fluorescence microscopy. Nat. Methods 2017, 14, 374–380. [Google Scholar] [CrossRef] [PubMed]
  18. Bond, C.Z.; Wizinowich, P.; Chun, M.; Mawet, D.; Lilley, S.; Cetre, S.; Jovanovic, N.; Delorme, J.; Wetherell, E.; Jacobson, S.M.; et al. Adaptive optics with an infrared pyramid wavefront sensor. Proc. SPIE 2018, 10703, 107031Z. [Google Scholar] [CrossRef]
  19. Mahajan, V.N. Zernike circle polynomials and optical aberration of system with circular pupils. Appl. Opt. 1994, 33, 8121–8124. [Google Scholar] [CrossRef] [PubMed]
  20. Love, G.D. Wavefront correction and production of Zernike modes with a Liquid crystal spatial light modulator. Appl. Opt. 1997, 36, 1517–1525. [Google Scholar] [CrossRef] [PubMed]
  21. Khonina, S.N.; Kotlyar, V.V.; Soifer, V.A.; Wang, Y.; Zhao, D. Decomposition of a coherent light field using a phase Zernike filter. Proc. SPIE 1998, 3573, 550–553. [Google Scholar]
  22. Neil, M.A.A.; Booth, M.J.; Wilson, T. New modal wave-front sensor: A theoretical analysis. J. Opt. Soc. Am. A 2000, 17, 1098–1107. [Google Scholar] [CrossRef]
  23. Booth, M.J. Direct measurement of Zernike aberration modes with a modal wavefront sensor. Proc. SPIE 2003, 5162, 79–90. [Google Scholar]
  24. Sheppard, C.J.R. Zernike expansion of pupil filters: Optimization of the signal concentration factor. J. Opt. Soc. Am. A 2015, 32, 928–933. [Google Scholar] [CrossRef] [PubMed]
  25. Porfirev, A.P.; Khonina, S.N. Experimental investigation of multi-order diffractive optical elements matched with two types of Zernike functions. Proc. SPIE 2016, 9807, 98070E. [Google Scholar]
  26. Wilby, M.J.; Keller, C.U.; Snik, F.; Korkiakoski, V.; Pietrow, A.G.M. The coronagraphic Modal Wavefront Sensor: A hybrid focal-plane sensor for the high-contrast imaging of circumstellar environments. Astron. Astrophys. 2017, 597, A112. [Google Scholar] [CrossRef] [Green Version]
  27. Khonina, S.N.; Karpeev, S.V.; Porfirev, A.P. Wavefront aberration sensor based on a multichannel diffractive optical element. Sensors 2020, 20, 3850. [Google Scholar] [CrossRef]
  28. Gerchberg, R.; Saxton, W. Phase determination for image and diffraction plane pictures in the electron microscope. Optik 1971, 34, 275–284. [Google Scholar]
  29. Fienup, J.R. Reconstruction of an object from the modulus of its Fourier transform. Opt. Lett. 1978, 3, 27–29. [Google Scholar] [CrossRef]
  30. Elser, V. Phase retrieval by iterated projections. J. Opt. Soc. Am. A 2003, 20, 40–55. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Marchesini, S.A. Unified evaluation of iterative projection algorithms for phase retrieval. Rev. Sci. Instrum. 2007, 78, 011301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Cheng, Z.; Meiqin, W.; Qianwen, C.; Dong, W.; Sui, W. Two-step phase retrieval algorithm using single-intensity measurement. Int. J. Opt. 2018, 2018, 8643819. [Google Scholar] [CrossRef] [Green Version]
  33. Tokovinin, A.; Heathcote, S. DONUT: Measuring optical aberrations from a single extrafocal image. Publ. Astron. Soc. Pac. 2006, 118, 1165–1175. [Google Scholar] [CrossRef] [Green Version]
  34. Guo, H.; Korablinova, N.; Ren, Q.; Bille, J. Wavefront reconstruction with artificial neural networks. Opt. Express 2006, 14, 6456–6462. [Google Scholar] [CrossRef] [PubMed]
  35. Paine, S.W.; Fienup, J.R. Machine learning for improved image-based wavefront sensing. Opt. Lett. 2018, 43, 1235–1238. [Google Scholar] [CrossRef] [PubMed]
  36. Rivenson, Y.; Zhang, Y.; Günaydın, H.; Teng, D.; Ozcan, A. Phase recovery and holographic image reconstruction using deep learning in neural networks. Light Sci. Appl. 2018, 7, 17141. [Google Scholar] [CrossRef]
  37. Nishizaki, Y.; Valdivia, M.; Horisaki, R.; Kitaguchi, K.; Saito, M.; Tanida, J.; Vera, E. Deep learning wavefront sensing. Opt. Express 2019, 27, 240–251. [Google Scholar] [CrossRef]
  38. Rodin, I.A.; Khonina, S.N.; Serafimovich, P.G.; Popov, S.B. Recognition of wavefront aberrations types corresponding to single Zernike functions from the pattern of the point spread function in the focal plane using neural networks. Comput. Opt. 2020, 44, 923–930. [Google Scholar] [CrossRef]
  39. Khorin, P.A.; Dzyuba, A.P.; Serafimovich, P.G.; Khonina, S.N. Neural networks application to determine the types and magnitude of aberrations from the pattern of the point spread function out of the focal plane. J. Phys. Conf. Ser. 2021, 2086, 012148. [Google Scholar] [CrossRef]
  40. Khorin, P.A. Analysis wavefront propagating in free space based on the Zernike polynomials and Gauss-Laguerre modes expansion. J. Phys. Conf. Ser. 2019, 1096, 012104. [Google Scholar] [CrossRef]
  41. Kirilenko, M.S.; Khorin, P.A.; Porfirev, A.P. A wavefront analysis based on Zernike polynomials. CEUR Workshop Proc. 2016, 1638, 66–75. [Google Scholar]
  42. Khorin, P.A.; Degtyarev, S.A. Wavefront aberration analysis with a multi-order diffractive optical element. CEUR Workshop Proc. 2017, 1900, 28–33. [Google Scholar]
  43. Khorin, P.A.; Volotovskiy, S.G. Analysis of the threshold sensitivity of a wavefront aberration sensor based on a multi-channel diffraction optical element. Proc. SPIE 2021, 11793, 117930B. [Google Scholar]
  44. Khorin, P.A.; Podlipnov, V.V.; Khonina, S.N. Generation of scalable wavefront for testing optical systems. Proc. SPIE 2020, 11516, 115161K. [Google Scholar]
  45. Born, M.; Wolf, E. Principles of Optics: Electromagnetic Theory of Propagation, Interference and Diffraction of Light, 7th ed.; Cambridge University Press: Cambridge, UK, 1999. [Google Scholar]
  46. Weaver, C.S.; Goodman, J.W. A technique for optically convolving two functions. Appl. Opt. 1966, 5, 1248–1249. [Google Scholar] [CrossRef] [PubMed]
  47. Horner, J.L.; Gianino, P.D. Phase-only matched filtering. Appl. Opt. 1984, 23, 812–816. [Google Scholar] [CrossRef]
  48. Millán, M.S. Advanced optical correlation and digital methods for pattern matching—50th anniversary of Vander Lugt matched filter. J. Opt. 2012, 14, 103001. [Google Scholar] [CrossRef]
  49. Khonina, S.N.; Karpeev, S.V.; Paranin, V.D. A technique for simultaneous detection of individual vortex states of Laguerre-Gaussian beams transmitted through an aqueous suspension of microparticles. Opt. Lasers Eng. 2018, 105, 68–74. [Google Scholar] [CrossRef]
  50. Khonina, S.N.; Karpeev, S.V.; Butt, M.A. Spatial-light-modulator-based multichannel data transmission by vortex beams of various orders. Sensors 2021, 21, 2988. [Google Scholar] [CrossRef] [PubMed]
  51. Khorin, P.A.; Volotovskiy, S.G.; Khonina, S.N. Optical detection of values of separate aberrations using a multi-channel filter matched with phase Zernike functions. Comput. Opt. 2021, 45, 525–533. [Google Scholar] [CrossRef]
  52. Khorin, P.A. Iterative algorithm for wavefront correction based on optical decomposition in wave aberrations. In Proceedings of the 2021 International Conference on Information Technology and Nanotechnology (ITNT), Samara, Russia, 20–24 September 2021; pp. 1–6. [Google Scholar]
Figure 1. Illustration of the compensating process for wave aberration W(x, y) of type (n, m) = (2, 0) with α = 0.5λ level (green highlights) by complex conjugate wave aberration Ψ p q * x , y   = Ψ 20 * x , y with different levels of α0 (orange highlights). Results of compensation (blue highlights) are placed at the intersection of corresponding rows and columns.
Figure 1. Illustration of the compensating process for wave aberration W(x, y) of type (n, m) = (2, 0) with α = 0.5λ level (green highlights) by complex conjugate wave aberration Ψ p q * x , y   = Ψ 20 * x , y with different levels of α0 (orange highlights). Results of compensation (blue highlights) are placed at the intersection of corresponding rows and columns.
Photonics 09 00204 g001
Figure 2. Illustration of the compensating process for wave aberration W(x, y) of type (n, m) = (2, 2) with α = 0.5λ level by complex conjugate wave aberration Ψ p q * x , y   = Ψ 22 * x , y with different levels of α0 (description for pictures are the same as for Figure 1).
Figure 2. Illustration of the compensating process for wave aberration W(x, y) of type (n, m) = (2, 2) with α = 0.5λ level by complex conjugate wave aberration Ψ p q * x , y   = Ψ 22 * x , y with different levels of α0 (description for pictures are the same as for Figure 1).
Photonics 09 00204 g002
Figure 3. Illustration of the compensating process for wave aberration W(x, y) of type (n, m) = (3, 1) with α = 0.5λ level (green highlights) by complex conjugate wave aberration Ψ p q * x , y   = Ψ 31 * x , y with different levels of α0 (orange highlights). Results of compensation (blue highlights) are placed at the intersection of corresponding rows and columns.
Figure 3. Illustration of the compensating process for wave aberration W(x, y) of type (n, m) = (3, 1) with α = 0.5λ level (green highlights) by complex conjugate wave aberration Ψ p q * x , y   = Ψ 31 * x , y with different levels of α0 (orange highlights). Results of compensation (blue highlights) are placed at the intersection of corresponding rows and columns.
Photonics 09 00204 g003
Figure 4. Examples of intensity patterns for different filters under plane wavefront illumination: (a) T-filter defined by Equation (14) at α0 = 0.5; (b) LT-filter defined by Equation (13), and (c) L-filter defined by Equation (15) at (p0, q0) = (3, 1).
Figure 4. Examples of intensity patterns for different filters under plane wavefront illumination: (a) T-filter defined by Equation (14) at α0 = 0.5; (b) LT-filter defined by Equation (13), and (c) L-filter defined by Equation (15) at (p0, q0) = (3, 1).
Photonics 09 00204 g004
Figure 5. Schematic of the experimental setup for detection of wavefront aberrations and their superpositions using a diffractive multichannel filter. Laser is a solid-state laser (λ = 532 nm); PH is a pinhole (hole size of 40 μm); L1, L2, L3, and L4 are spherical lenses (f1 = 350 mm, f2 = 300 mm, f3 = 200 mm, and f4 = 250 mm); SLM1 is a transmissive spatial light modulator (HOLOEYE LC 2012); SLM2 is a reflective spatial light modulator (HOLOEYE PLUTO VIS); D1 and D2 are circular apertures; M1 and M2 are mirrors; CAM is a ToupCam UCMOS08000KPB video camera.
Figure 5. Schematic of the experimental setup for detection of wavefront aberrations and their superpositions using a diffractive multichannel filter. Laser is a solid-state laser (λ = 532 nm); PH is a pinhole (hole size of 40 μm); L1, L2, L3, and L4 are spherical lenses (f1 = 350 mm, f2 = 300 mm, f3 = 200 mm, and f4 = 250 mm); SLM1 is a transmissive spatial light modulator (HOLOEYE LC 2012); SLM2 is a reflective spatial light modulator (HOLOEYE PLUTO VIS); D1 and D2 are circular apertures; M1 and M2 are mirrors; CAM is a ToupCam UCMOS08000KPB video camera.
Photonics 09 00204 g005
Figure 6. Detailed results of the analysis of the WF distorted by the odd-type aberration (n, m) = (3, 1) with α = 0.5λ, using the LT-filter: (a) the phase of the analyzed WF; (b) simulation of the filter action; (c) optical experiment.
Figure 6. Detailed results of the analysis of the WF distorted by the odd-type aberration (n, m) = (3, 1) with α = 0.5λ, using the LT-filter: (a) the phase of the analyzed WF; (b) simulation of the filter action; (c) optical experiment.
Photonics 09 00204 g006
Figure 7. Detailed results of the analysis of the WF distorted by the even type of aberration (n, m) = (2, 2) with α = 0.5λ, using the LT-filter: (a) the phase of the analyzed WF; (b) simulation of the filter action; (c) optical experiment.
Figure 7. Detailed results of the analysis of the WF distorted by the even type of aberration (n, m) = (2, 2) with α = 0.5λ, using the LT-filter: (a) the phase of the analyzed WF; (b) simulation of the filter action; (c) optical experiment.
Photonics 09 00204 g007
Figure 8. Results of the analysis of the WF distorted by the aberration (n, m) = (3, 3) with α = 0.7λ, using the LT-filter: (a) the phase of the analyzed WF; (b) simulation of the filter action; (c) optical experiment.
Figure 8. Results of the analysis of the WF distorted by the aberration (n, m) = (3, 3) with α = 0.7λ, using the LT-filter: (a) the phase of the analyzed WF; (b) simulation of the filter action; (c) optical experiment.
Photonics 09 00204 g008
Figure 9. Results of modeling and optical experiment in detecting WF distorted by the superposition of aberrations 0.5λ*(2, 2) + 0.5λ*(3, 3) using the LT-filter (change in the DO orientation—dashed line): (a) the phase of the analyzed WF; (b) simulation of the filter action; (c) optical experiment.
Figure 9. Results of modeling and optical experiment in detecting WF distorted by the superposition of aberrations 0.5λ*(2, 2) + 0.5λ*(3, 3) using the LT-filter (change in the DO orientation—dashed line): (a) the phase of the analyzed WF; (b) simulation of the filter action; (c) optical experiment.
Photonics 09 00204 g009
Figure 10. Results of modeling and optical experiment in detecting WF with residual aberration (n, m) = (2, 2) with α = 0.5λ level using the LT-filter (change in DO orientation—dashed line; correlation peak—solid line): (a) the phase of the analyzed WF; (b) simulation of the filter action; (c) optical experiment.
Figure 10. Results of modeling and optical experiment in detecting WF with residual aberration (n, m) = (2, 2) with α = 0.5λ level using the LT-filter (change in DO orientation—dashed line; correlation peak—solid line): (a) the phase of the analyzed WF; (b) simulation of the filter action; (c) optical experiment.
Photonics 09 00204 g010
Figure 11. Results of modeling and optical experiment in detecting WF distorted by the superposition of aberrations 0.5λ*(2, 0) + 0.5λ*(2, 2) + 0.5λ*(3, 3) using the LT-filter (change in the DO orientation—dashed line): (a) the phase of the analyzed WF; (b) simulation of the filter action; (c) optical experiment.
Figure 11. Results of modeling and optical experiment in detecting WF distorted by the superposition of aberrations 0.5λ*(2, 0) + 0.5λ*(2, 2) + 0.5λ*(3, 3) using the LT-filter (change in the DO orientation—dashed line): (a) the phase of the analyzed WF; (b) simulation of the filter action; (c) optical experiment.
Photonics 09 00204 g011
Figure 12. The results of modeling the superposition of the odd wave aberration of the type (n1, m1) = (3, 1) with α1 = 1.0λ level (green highlights) and even wave aberrations of the type (n2, m2) = (2, 2) with different levels of α2 (orange highlights). Results of compensation (blue highlights) are placed at the intersection of corresponding rows and columns.
Figure 12. The results of modeling the superposition of the odd wave aberration of the type (n1, m1) = (3, 1) with α1 = 1.0λ level (green highlights) and even wave aberrations of the type (n2, m2) = (2, 2) with different levels of α2 (orange highlights). Results of compensation (blue highlights) are placed at the intersection of corresponding rows and columns.
Photonics 09 00204 g012
Table 1. Correspondence of Zernike functions of double indices (n, m) to the aberration type.
Table 1. Correspondence of Zernike functions of double indices (n, m) to the aberration type.
(n, m)(0, 0)(1, 1)(2, 2)(2, 0)(3, 3)(3, 1)(4, 4)(4, 2)(4, 0)
Z n m x , y 1 2 r cos ( φ ) 6 r 2 cos ( 2 φ ) 3 ( 2 r 2 1 ) 2 2 r 3 cos ( 3 φ ) 2 2 ( 3 r 3 2 r ) cos ( φ ) 10 r 4 cos ( 4 φ ) 10 ( 4 r 4 3 r 2 ) cos ( 2 φ ) 5 ( 6 r 4 6 r 2 + 1 )
AberrationtypeIdealTiltAstigmatismDefocusComa (Trefoil)Pure comaQuadrofoil2nd order astigmatismSpherical
Table 2. Simulation results for radial wave aberration of type (n, m) = (2, 0) and (n, m) = (4, 0).
Table 2. Simulation results for radial wave aberration of type (n, m) = (2, 0) and (n, m) = (4, 0).
Aberration Typeα0.10λ0.25λ0.50λ0.75λ1.00λ
Defocus
(n, m) = (2, 0)
WF phase Photonics 09 00204 i001 Photonics 09 00204 i002 Photonics 09 00204 i003 Photonics 09 00204 i004 Photonics 09 00204 i005
PSF Photonics 09 00204 i006 Photonics 09 00204 i007 Photonics 09 00204 i008 Photonics 09 00204 i009 Photonics 09 00204 i010
Spherical
(n, m) = (4, 0)
WF phase Photonics 09 00204 i011 Photonics 09 00204 i012 Photonics 09 00204 i013 Photonics 09 00204 i014 Photonics 09 00204 i015
PSF Photonics 09 00204 i016 Photonics 09 00204 i017 Photonics 09 00204 i018 Photonics 09 00204 i019 Photonics 09 00204 i020
Table 3. Results of modeling an even wave aberration of the type (n, m) = (2, 2) and (n, m) = (4, 4).
Table 3. Results of modeling an even wave aberration of the type (n, m) = (2, 2) and (n, m) = (4, 4).
Aberration Typeα0.10λ0.25λ0.50λ0.75λ1.00λ
Astigmatism
(n, m) = (2, 2)
WF phase Photonics 09 00204 i021 Photonics 09 00204 i022 Photonics 09 00204 i023 Photonics 09 00204 i024 Photonics 09 00204 i025
PSF Photonics 09 00204 i026 Photonics 09 00204 i027 Photonics 09 00204 i028 Photonics 09 00204 i029 Photonics 09 00204 i030
Quadrofoil
(n, m) = (4, 4)
WF phase Photonics 09 00204 i031 Photonics 09 00204 i032 Photonics 09 00204 i033 Photonics 09 00204 i034 Photonics 09 00204 i035
PSF Photonics 09 00204 i036 Photonics 09 00204 i037 Photonics 09 00204 i038 Photonics 09 00204 i039 Photonics 09 00204 i040
Table 4. Modeling results for odd wave aberration of type (n, m) = (3,1) and (n, m) = (3, 3).
Table 4. Modeling results for odd wave aberration of type (n, m) = (3,1) and (n, m) = (3, 3).
Aberration Typeα0.10λ0.25λ0.50λ0.75λ1.00λ
Pure coma
(n, m) = (3, 1)
WF phase Photonics 09 00204 i041 Photonics 09 00204 i042 Photonics 09 00204 i043 Photonics 09 00204 i044 Photonics 09 00204 i045
PSF Photonics 09 00204 i046 Photonics 09 00204 i047 Photonics 09 00204 i048 Photonics 09 00204 i049 Photonics 09 00204 i050
Coma
(Trefoil)
(n, m) = (3, 3)
WF phase Photonics 09 00204 i051 Photonics 09 00204 i052 Photonics 09 00204 i053 Photonics 09 00204 i054 Photonics 09 00204 i055
PSF Photonics 09 00204 i056 Photonics 09 00204 i057 Photonics 09 00204 i058 Photonics 09 00204 i059 Photonics 09 00204 i060
Table 5. Results of modeling and optical experiment of the operation of multichannel T-filters in detecting WF distorted by coma-type aberration (n, m) = (3, 1) with α = 0.5λ.
Table 5. Results of modeling and optical experiment of the operation of multichannel T-filters in detecting WF distorted by coma-type aberration (n, m) = (3, 1) with α = 0.5λ.
WFType of the FilterIntensity in the Focal Plane of the T-Filter Equation (14)
ModelingExperiment
Photonics 09 00204 i061
(n, m) = (0, 0);
α = 0
T-filter
at α0 = 0.5λ
Photonics 09 00204 i062 Photonics 09 00204 i063
Photonics 09 00204 i064
(n, m) = (3, 1);
α = 0.5λ
T-filter
at α0 = 0.3λ
Photonics 09 00204 i065 Photonics 09 00204 i066
T-filter
at α0 = 0.5λ
Photonics 09 00204 i067 Photonics 09 00204 i068
T-filter
at α0 = 0.7λ
Photonics 09 00204 i069 Photonics 09 00204 i070
Table 6. Results of modeling and optical experiment of the operation of multichannel T-filters in detecting WF distorted by aberration (n, m) = (3, 3) with α = 0.7λ.
Table 6. Results of modeling and optical experiment of the operation of multichannel T-filters in detecting WF distorted by aberration (n, m) = (3, 3) with α = 0.7λ.
Aberrated WFType of the FilterIntensity in the Focal Plane of the T-Filter Equation (14)
ModelingExperiment
Photonics 09 00204 i071
(n, m) = (3, 3);
α = 0.5λ
T-filter
at α0 = 0.5λ
Photonics 09 00204 i072 Photonics 09 00204 i073
T-filter
at α0 = 0.7λ
Photonics 09 00204 i074 Photonics 09 00204 i075
T-filter
at α0 = 0.9λ
Photonics 09 00204 i076 Photonics 09 00204 i077
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Khorin, P.A.; Porfirev, A.P.; Khonina, S.N. Adaptive Detection of Wave Aberrations Based on the Multichannel Filter. Photonics 2022, 9, 204. https://doi.org/10.3390/photonics9030204

AMA Style

Khorin PA, Porfirev AP, Khonina SN. Adaptive Detection of Wave Aberrations Based on the Multichannel Filter. Photonics. 2022; 9(3):204. https://doi.org/10.3390/photonics9030204

Chicago/Turabian Style

Khorin, Pavel A., Alexey P. Porfirev, and Svetlana N. Khonina. 2022. "Adaptive Detection of Wave Aberrations Based on the Multichannel Filter" Photonics 9, no. 3: 204. https://doi.org/10.3390/photonics9030204

APA Style

Khorin, P. A., Porfirev, A. P., & Khonina, S. N. (2022). Adaptive Detection of Wave Aberrations Based on the Multichannel Filter. Photonics, 9(3), 204. https://doi.org/10.3390/photonics9030204

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop