Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Holographic optical elements with a large adjustable focal length and an aberration correction

Open Access Open Access

Abstract

Holographic optical element (HOE) has the advantages of light weight, small volume and multiple functions, but the fixed focal length limits its application. The current methods of adjusting focal length suffer from the disadvantages of small adjustable range and the introduction of aberration. This paper proposes a design method for HOE, and it can compensate the aberrations of the reconstructed image for each distance in a large adjustable range. The focal length of HOE is adjusted by modulating the incident light through the wavefront modulator, and the aberration of the reconstructed image is corrected by optimizing the phase distribution of HOE and superimposing different compensation phases for different reconstruction distances. The reconstructed image experimentally moves without aberration in a large range, which demonstrates the feasibility of the proposed method. It is expected to be widely used in various optical fields where the focal length of HOE needs to be dynamically tuned.

© 2022 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

Diffractive optical element (DOE) is a kind of element that regulate the light field through microstructure. It is widely used in beam splittin [13], beam shaping [4,5] and other aspects [6]. Focal length is an important parameter of DOE, which affects the distance from the reconstructed image to DOE. Adjusting the focal length of DOE is of great significance: it can expand the function of DOE and improve the performance of optical systems. As a kind of easily integrated focal length adjustable element, DOE is widely required in eye-tracking [79], micro-machining, particle manipulation [10,11] and many other tasks, especially in the systems for holographic projection [12,13], near eye display [1417], imaging [18], microscopy [1921], optical fiber sensing [22], integrated optoelectronic [23] and so on.

The traditional methods of adjusting the focal length can be classified into two categories: manipulating the element or modulating the wavefront. Rotating two moiré elements by a certain angle can adjust the focal length around the designed focus [2428], but this method needs mechanical movement and the adjustable range of focal length is limited. Stretching the stretchable element is another way to manipulating the element [29], besides the above problem, it further suffers from the low stability. Modulating the wavefront can also realize the adjustable focal length [3034]. The usual implementation method is to add a wavefront modulator in the optical path, such as liquid crystal lens or liquid lens, to change the wavefront curvature. However, the wavefront modulator will introduce aberration, and the larger the adjustable range, the greater the aberration.

Holographic optical element (HOE) is a kind of DOE based on holography, which not only has the characteristics of traditional DOE, but also has the advantages of low manufacturing cost and the potential of multiplexing. Correct the aberration while adjusting the focal length in a large range is the hardest part in the design and fabrication of HOE. Researchers have studied the aberrations of HOE in head mounted display and waveguide. The aberrations are mainly analyzed by formula derivation [35,36], ray tracing [3739] and wavefront sensor detection [40], and the aberrations are corrected mainly by introducing cylindrical mirror [36] in the recording optical path or making pre-compensation [3739] in the image source. However, global optimization for multiple aberrations and depths of HOE is still a problem to be solved.

This paper proposes a design method for HOE, which can adjust the focal length of HOE in a wide range and correct the aberrations at the same time. In this paper, formula derivation, numerical simulation and optical experiments are carried out.

This paper is arranged as follows: the second chapter describes the basic ideas and principles, deduces the relationship between the focal length of HOE and the distance of reconstructed image, designs a general algorithm to optimize the phase distribution of HOE, and analyzes the principles and compensation methods of aberration caused by oblique incidence; the third chapter shows and analyzes the experimental results of correcting the aberration at the design focal length and in the process of adjusting the focal length.

2. Principle of adjusting the focal length and correcting aberrations

The function of HOE is to reconstruct the target wavefront at the designed distance. As shown in Fig. 1(a), the HOE is an off-axis reflective DOE made by holographic principle and the fabrication process can be found in [41]. According to the basic principle of holography [42], if the illumination wave is the same as the reference wave, the object will be reconstructed at the designed distance. But when the wavefront of the illumination wave is modulated, the reconstructed wave will also be changed. The movement of the reconstructed image means the change of the focal length. To adjust the focal length of HOE in a wide range, the wavefront of the illumination wave needs to be modulated. As shown in Fig. 1(b), the spatial light modulator (SLM) is selected as the wavefront modulator due to its excellent modulation ability. When the modulation factor is loaded on the SLM, an additional phase distribution will be superimposed on the reconstructed wave. The parameters of modulation factor determine the curvature radius of illumination wave, which is related to the focal length of HOE. When the parameters of the modulation factor change, the position of the reconstructed image will be adjusted to different distances. It is worth mentioning that, for simplicity and clarity, only the positions of three reconstruction planes are drawn in Fig. 1(b). However, due to the high refresh rate of SLM and the continuity of compensation parameters, the number of reconstruction plane is unlimited, which can be continuously adjusted within a range, and the aberration at each position can be compensated.

 figure: Fig. 1.

Fig. 1. (a) shows the reconstruction process of HOE and (b) illustrates the process of focal length adjustment and aberration correction by SLM. $IW$: the illumination wave, $RI$: the reconstructed image, ${l_d}$: the designed focal length, $MF$: the phase distribution of modulation factor, $SLM$: spatial light modulator, ${I_1}{I_3}$: reconstructed images at different distances.

Download Full Size | PDF

When the focal length of HOE is adjusted, the reconstructed image will have aberrations such as distortion, field curvature and astigmatism. Part of it comes from the manufacturing process and can be compensated by optimizing the phase distribution of HOE. The other part comes from the parameters of the modulation factor during reconstruction and the oblique incidence of the illumination wave, which can be corrected by superimposing the compensation phase on the SLM. The following sub-sections will describe the principles of adjusting the focal length of HOE and correcting aberration respectively.

2.1 Adjust the focal length of HOE

In order to describe the relationship between the parameters of illumination wave and the focal length of HOE accurately, the expression is derived based on Fresnel diffraction integral formula. Suppose the coordinates of the target light field are $({{x_1},{y_1}} )$, the coordinates of HOE are $({x,y} )$, and the coordinates of the reconstructed light field are $({x^{\prime},y^{\prime}} )$. Given that the designed focal length is ${l_d}$, the wavefront distribution of the target object is $o({{x_1},{y_1}} )$, the complex amplitude of the reference wave and the illumination wave are $R({x,y} )= |R |\cdot exp({i{\varphi_r}} )$ and $C({x,y} )= |C |\cdot exp({i{\varphi_c}} )$ respectively, where |R| and |C| represent the amplitude while ${\varphi _r}$ and ${\varphi _c}$ represent the phase of wavefront. The complex amplitude of the reconstruction image can be written as

$$\scalebox{0.88}{$\displaystyle E({x^{\prime},y^{\prime}} )= {c_1}{c_2}|R ||C |\cdot \mathrm{{\cal F}}\left\{ {exp\left[ {ik\left( { - \frac{1}{{{l_d}}} - \frac{1}{{{l_c}}} + \frac{1}{{{l_i}}}} \right) \cdot \frac{{{x^2} + {y^2}}}{2}} \right]} \right\} \otimes \mathrm{{\cal F}}\left\{ {{\mathrm{{\cal F}}^{ - 1}}\left\{ {o({{x_1},{y_1}} )exp\left( { - ik\frac{{{x_1}^2 + {y_1}^2}}{{2{l_d}}}} \right)} \right\}} \right\},$}$$
where ${c_1} = {{exp({ - ik{l_d}} )} / {({i\lambda {l_d}} )}}$, ${c_2} = {{exp\{{ik({{l_i} + {{[{{{({x^{\prime}} )}^2} + {{({y^{\prime}} )}^2}} ]} / {({2{l_i}} )}}} )} \}} / {({i\lambda {l_i}} )}}$, $\mathrm{{\cal F}}\{{\cdot} \}$ and ${\mathrm{{\cal F}}^{ - 1}}\{{\cdot} \}$ represent Fourier transform and inverse Fourier transform respectively, ${\otimes} $ represents convolution operation. In order to reconstruct the object wave, the coefficient of the exponential term before the convolution symbol should be zero, which can be deduced
$$\frac{1}{{{l_i}}} = \frac{1}{{{l_c}}} + \frac{1}{{{l_d}}}.$$

Equation (2) shows the relationship between the distance of the reconstructed image and the focal length of HOE, where the object distance is infinite, and $1/({1/{l_c} + 1/{l_d}} )$ is equivalent to the focal length of HOE. The designed focal length ${l_d}$ is determined when designing the phase distribution, so the distance ${l_i}$ of the reconstructed image can be adjusted by changing the parameter ${l_c}$ of the illumination wave.

2.2 Correct the aberrations at the designed focal length

The aberration correction method in this paper mainly has two steps. One is to correct the aberrations at the designed focal length by optimizing the DOE phase distribution before manufacturing, and the other is to compensate the aberrations at each focal length by wavefront modulation element when adjusting the focal length. This section will describe the causes and correction methods of aberrations at the designed focal length.

SLM is used to load the phase distribution, so there exists zero-order light in the reconstructed image due to the dead zone of SLM. It will block the reconstructed image and reduce its diffraction efficiency. However, if the zero-order is directly filtered by high pass filter in the spatial spectrum plane of $4f$ system, the middle region of the reconstructed image will be missing. In this paper, a quadratic phase factor is superimposed on the original phase distribution, resulting in the separation of signal light and zero-order light in depth. Therefore, the zero-order light can be filtered out and has little effect on the signal light.

In addition, distortion is another aberration that affects the quality of reconstructed image. It is mainly introduced by the optical elements and can be corrected by pre-distorting the target image in the process of calculating the phase distribution of DOE [43]. The mathematical models of radial and tangential distortion can be written as

$$\begin{array}{l} \left\{ \begin{array}{l} {x_d} = x({1 + {k_1}{r^2} + {k_2}{r^4} + {k_3}{r^6}} )\\ {y_d} = y({1 + {k_1}{r^2} + {k_2}{r^4} + {k_3}{r^6}} )\end{array} \right.\\ \left\{ \begin{array}{l} {x_d} = x + [{2{p_1}y + {p_2}({{r^2} + 2{x^2}} )} ]\\ {y_d} = y + [{2{p_2}x + {p_1}({{r^2} + 2{y^2}} )} ]\end{array} \right., \end{array}$$
where ${r^2} = {x^2} + {y^2}$ and $[{{k_1},{k_2},{k_3},{p_1},{p_2}} ]$ constitutes the distortion vector. Since the tangential distortion introduced in the actual system is not obvious and the radial distortion is relatively weak, only ${k_1}$ needs to be determined. The initial value of the correction parameter is calculated from the coordinates of the distortion image and the original image, and the evaluation function is defined as
$${E_{dist}} = \sum\limits_i {\sqrt {{{({{x_i} - {x_c}} )}^2} + {{({{y_i} - {y_c}} )}^2}} } ,$$
where $({{x_c},{y_c}} )$ is the central coordinate of the reconstructed image and ${\; }({{x_i},{y_i}} )$ is the coordinate of the $\mathop i\nolimits^{th}$ grid intersection. The larger the ${E_{dist}}$, the more serious the distortion. Therefore, the initial parameters are optimized to reduce the value of the evaluation function.

Astigmatism and field curvature are also aberrations that can’t be ignored. Astigmatism is introduced by the eccentricity and inclination of optical elements, which will make the horizontal and vertical lines in the reconstructed image not clear at the same time. Field curvature is introduced by the wavefront distortion caused by optical elements, and the existence of field curvature makes the center and edge of the reconstructed image not clear at the same time. Although they are formed for different reasons, astigmatism and field curvature can be corrected by superimposing compensation phase on the phase distribution of DOE. The compensated phase can be expressed as

$$exp\left[ {\frac{{i\pi }}{\lambda }\left( {\frac{{{x^2}}}{{{l_{ax}}}} + \frac{{{y^2}}}{{{l_{ay}}}}} \right)} \right] \cdot exp\left[ {\frac{{i\pi }}{{\lambda {l_v}}}({{x^2} + {y^2}} )} \right].$$

The two exponential factors in Eq. (5) are used to correct astigmatism and field curvature respectively. Where ${l_{ax}}$ and ${l_{ay}}$ are the correction parameters in the x-z plane and y-z plane respectively, ${l_v}$ is the correction parameter of field curvature.

The evaluation function of astigmatism and field curvature can be defined as Eq. (6) and Eq. (7) respectively.

$${E_{asti}} = 1 - \frac{{|{\overline {{h_x}} - \overline {{h_y}} } |}}{{\overline {{h_x}} + \overline {{h_y}} }}.$$
$\overline {{h_x}} $ and $\overline {{h_y}} $ are the average widths of the horizontal and vertical lines of the reconstructed image, respectively.
$${E_{curv}} = 1 - \frac{{|{\overline {{h_{ex}}} - \overline {{h_{cx}}} } |}}{{\overline {{h_{ex}}} + \overline {{h_{cx}}} }}.$$
$\overline {{h_{ex}}} $ and ${h_{cx}}$ are the average width of the line at the edge and center of the reconstructed image, respectively. When optimizing the parameters, the larger the values of ${E_{asti}}$ and ${E_{curv}}$, the weaker the aberration.

Since the correction of distortion, astigmatism and field curvature involves the optimization of compensation parameters, the general process of determining compensation parameters is introduced.

Figure 2 shows the process of optimizing parameters by dichotomy. The search range as input can be obtained by two methods, one is to carry out several pre-experiments and take the parameters over and under corrected aberration, the other is to introduce the camera and image processing program into the optimization process. Firstly, an initial value k and a search step s are selected to calculate value of the evaluation function when the parameters are k, $k + s$ and $k - s$ respectively, and then the direction in which the evaluation function decreases is selected as the search direction. When the evaluation function increases, $k \pm ns$ and k are taken as the two endpoints of the search range, where n is the number of optimization cycle. In the optimization process, the parameters for correcting aberrations can be found.

 figure: Fig. 2.

Fig. 2. Flow chart for determining correction parameters.

Download Full Size | PDF

2.3 Correct the aberrations during focal length adjustment

Because the phase factors with different curvature radii are superimposed on the x-z and y-z planes during recording, if the parameters required to adjust the focal length are still calculated by Eq. (2), astigmatism will appear, and the astigmatism corresponding to different distances is different. To solve this problem, Eq. (2) is deformed to

$$\frac{1}{{{l_i}}} = \frac{1}{{{l_d}}} + \frac{1}{{{l_{dst(m )}}}} + \frac{1}{{{l_{dcv(m )}}}} + \frac{1}{{{l_{c(m )}}}}.$$
$${f_{(m )}} = {1 / {\left( {\frac{1}{{{l_d}}} + \frac{1}{{{l_{dst(m )}}}} + \frac{1}{{{l_{dcv(m )}}}} + \frac{1}{{{l_{c(m )}}}}} \right)}}.$$

The focal length f of HOE can be calculated by Eq. (9), where m can be x or y, representing the parameters in the x-z or y-z plane respectively. The parameter ${l_c}$ of the modulated illumination wave can be calculated according to the Eq. (8). The parameters in the x-z and y-z planes vary with the distance, and the corresponding parameters in the two planes are different at the same distance.

High diffraction efficiency can be obtained by using reflective HOE, but the illumination wave needs to be incident obliquely. If the illumination wave is an inclined spherical wave, the wavefront curvature on the x-z and y-z planes are different when incident on HOE, which will cause astigmatism. The greater the tilt angle and the larger the size of HOE, the more obvious the astigmatism will be.

As shown in Fig. 3, it is assumed that the transverse dimension of HOE is ${d_{hy}}$, the curvature radius of the illumination wave corresponding to the center is ${l_c}$, the included angle between the illumination wave and the z-axis is $\theta $.

 figure: Fig. 3.

Fig. 3. The oblique incidence of illumination wave into HOE introduces astigmatism. (b) is the top view of (a).

Download Full Size | PDF

As described in 2.1, to adjust the focal length of HOE, a quadratic factor needs to be superimposed on the original phase distribution. The quadratic factor is equivalent to the phase distribution of a spherical wave perpendicular to the surface of HOE. However, the actual illumination wave is the inclined spherical wave ${\varphi _{ts}}$. when it subtracts the inclined phase factor ${\varphi _{tp}}$ of the reference wave, it is not equal to the spherical wave ${\varphi _s}$ which perpendicular to the surface of HOE, and the phase difference between them is ${\varphi _d}$.

$$\left\{ \begin{array}{l} {\varphi_{ts}} ={-} k\sqrt {{{({{y_h} + {l_{c - SLM}}sin\theta } )}^2} + {{({{l_{c - SLM}}cos\theta } )}^2}} \\ {\varphi_{tp}} = k\left( {\frac{{{d_{hy}}}}{2} - {y_h}} \right)sin\theta \\ {\varphi_s} ={-} k\sqrt {{x_h}^2 + {y_h}^2 + {l_{c - SLM}}^2} \\ {\varphi_d} = {\varphi_{ts}} - {\varphi_{tp}} - {\varphi_s} \end{array} \right..$$

As shown in Fig. 4, the numerical simulation results show that the actual and ideal illumination wavefront differ by a cylindrical phase factor.

 figure: Fig. 4.

Fig. 4. Phase difference caused by oblique incidence of illumination light

Download Full Size | PDF

To pre-compensate the phase difference caused by oblique incidence, the phase of illumination wave corresponding to the center of HOE is taken out, converted into height value, and then the illumination light curvature in this direction is calculated through geometric relationship. According to the above steps, the curve of the curvature radius of HOE in the y-z plane with the curvature radius of illumination light in the y-z plane can be drawn by numerical simulation.

First, a series of discrete data are calculated, as shown by the red dot in Fig. 5. Then, these discrete data are fitted to obtain the mathematical expression between $l{}_{c - HOE}$ and $l{}_{c - SLM}$. With this mathematical expression, phase parameter $l{}_{c - SLM}$ that needs to be loaded on SLM for any known $l{}_{c - HOE}$ can be directly calculated, instead of calling the numerical simulation program again when the focal length needs to be adjusted, which can improve the speed of adjustment. Since the distribution of discrete data is approximately linear, we choose the first-order equation for fitting and the coefficient can be obtained.

 figure: Fig. 5.

Fig. 5. The actual radius of curvature on HOE varies with the radius of curvature of the illumination light with the influence of oblique incidence angle.

Download Full Size | PDF

As shown in Fig. 5, the slope is related to the inclination angle of the illumination wave. The greater the inclination angle is, the greater the slope is. The slope of y-z plane is expected to be close to that of x-z plane, which means that the astigmatism introduced by oblique incidence is small and the compensation effect is better. However, the inclination angle of the illumination wave cannot be too small, which will overlap with the area where the reconstructed image is located. Considering the above two points and the convenience of setting the recording light path, 45° is choosed as the inclination angle of the illumination wave.

3. Experimental results

In the second section, we divide the aberration correction of HOE into the correction at the design focal length and the correction in the adjustment process. This is because the recording and reconstruction process of HOE is independent. The recording process is to “freeze” the wavefront of the object wave in the holographic material through interference, and then “activate” the wavefront of the object wave by the illumination wave in the reconstruction process, to make it continue to propagate. The aberration correction at the design focal length is the correction in the recording process. After the correction, the aberration will not be introduced into the subsequent reconstruction process. However, the reconstruction process itself will also introduce new aberrations due to oblique incidence, and this aberration will change with the adjustment of focal length. Therefore, in the reconstruction process, it is necessary to compensate for each reconstruction distance separately. Therefore, the following sub-sections will correct the aberration in the order of first at the design focal length and then in the adjustment process according to the above description.

3.1 Correct the aberrations at the designed focal length

Figure 6 shows the process and results of correcting the aberrations at the designed focal length. Figure 6(b) is a reconstructed image before correction. Through the method described in section 2, zero-order light, distortion, astigmatism and field curvature are eliminated in Fig. 6(c)-(g) respectively. This can be seen more clearly through the partially enlarged picture. The difference between the width of horizontal and vertical line in Fig. 6(j) is obviously less than Fig. 6(h), indicating that astigmatism has been corrected. The comparison between Fig. 6(i) and Fig. 6(j) shows that the difference between the width of horizontal and vertical line at the edge is greater than that at the center, that is, there is still field curve in the reconstructed image, Fig. 6(k) and Fig. 6(l) show that the width of horizontal and vertical line at the edge and center is basically the same, indicating that the field curve is corrected.

 figure: Fig. 6.

Fig. 6. Aberration correction at the designed focal length. (a) Original image. (b) the reconstructed image before correction. (c) the reconstructed image after zero-order light elimination. (d) the original image with pre-distortion. (e) - (g) represents the reconstructed image after correcting distortion, astigmatism and field curvature respectively. (h) - (l) is a partial enlarged view of (e) - (g).

Download Full Size | PDF

The reconstructed image before and after correction is evaluated by image processing program. The value of distortion evaluation function decreases from 4337.9 to 2865.6, the value of astigmatism evaluation function increases from 0.72 to 0.93, and the field curve evaluation function increases from 0.89 to 0.97. It shows that the distortion, astigmatism and field curve of the reconstructed image have been well corrected.

3.2 Correct the aberrations during focal length adjustment

Figure 1 shows the optical path of the experiment. The illumination wave is inclined and its wavefront is modulated by SLM. The reconstructed wave is emitted perpendicular to the surface of HOE, and a clear reconstructed image is formed at the position satisfying Eq. (8).

In the experiment, the designed reconstruction distance ${l_d}$ is 40cm, the transverse dimension ${d_{hy}}$ of HOE is 2.15mm, and the included angle $\theta $ between illumination light and z-axis is 45°. The expression fitted from the above data is

$${l_{c - HOE}} = 1.22 \times {l_{c - SLM}} - 1.35 \times {10^{ - 4}}.$$

The curves of illumination wave parameters of HOE and SLM with distance are drawn respectively, as shown in Fig. 7.

 figure: Fig. 7.

Fig. 7. Parameters for compensating phase distribution

Download Full Size | PDF

As shown in Fig. 8, the reconstructed image is adjusted to 30 cm, 40 cm and 50 cm respectively, and the astigmatism introduced by focal length adjustment and inclined illumination wave is well corrected. To present the compensation results more clearly, the part at the same position in the reconstructed image is enlarged as a comparison.

 figure: Fig. 8.

Fig. 8. Reconstructed images at different distances after aberration compensation. The reconstructed images before and after compensation are shown in (a) - (c) and (d) - (f) respectively, where (a) (d), (b) (e) and (c) (f) correspond to the reconstructed images adjusted to 30 cm, 40 cm and 50 cm respectively.

Download Full Size | PDF

In the experiment, the reconstructed image is moved in the range of 30∼50 cm by changing the focal lengths of x-z and y-z planes at 54.55∼200.00 cm and 23.35∼37.71 cm respectively, but due to the excellent wavefront modulation ability of SLM, theoretically, the adjustment range of the focal length of HOE is unlimited. However, some practical factors need to be considered in the experiment. For example, in the process of gradually increasing the reconstruction distance, the energy of the reconstructed pattern will gradually decrease due to the loss of power in the process of propagation and the expansion of the image size. In addition, the reduction of reconstruction distance will also lead to the reduction of reconstruction image size. Therefore, the adjustable range of focal length needs to be comprehensively considered according to the needs of the optical system. The actual optical system usually requires an adjustable range of no more than a few meters, and the proposed method can meet the needs of most actual optical systems. After the adjustable range of focal length is determined, the designed focal length of HOE can be selected in the middle of the adjustable range for subsequent adjustment and compensation.

Since there is no lens in the reconstruction process, distortion and field curvature are not introduced. Among the aberrations introduced in the reconstruction process, the astigmatism at different reconstruction distances caused by the oblique incidence of illumination light is the aberration that has the greatest impact on the quality of the reconstructed image, which is the main aberration concerned in this paper. The first row and the second row in Fig. 8 are the reconstructed images before and after correction, respectively. To fully illustrate the effect of astigmatism correction, the change of its evaluation function is calculated through the image processing program. The value of astigmatism evaluation function changes from 0.85 to 0.90 at 30cm, from 0.83 to 0.93 at 40cm and from 0.81 to 0.89 at 50cm. When the value of astigmatism evaluation function is greater than 0.9, it can be considered that astigmatism is not obvious. The experiment results show that the method in this paper can well compensate the aberrations at each focal length.

4. Conclusion

This paper proposes a method to correct the aberrations of HOE in the process of large-scale focal length adjustment. The relationship between the focal length of HOE and the parameters of modulation factor is derived. Combined with SLM which has excellent wavefront modulation ability, the theoretically designed adjustable range of focal length can be selected according to the requirements of optical system. The aberration introduced in the fabrication process is corrected by optimizing the phase distribution. The values of the evaluation functions of distortion, astigmatism and field curvature change experimentally by 33.94%, 29.17% and 8.99% respectively in the direction of aberration reduction. The effects of the parameters of modulation factor on aberration in the reconstruction process are analyzed. The value of astigmatism evaluation function of corrected reconstructed image increases by 6.00%, 11.94% and 10.44% at 30 cm, 40 cm and 50 cm in optical experiment, respectively. The increased values are close to or more than 0.9, which shows that the proposed method can realize the HOE with large range adjustment of focal length and aberration correction. In future, the design method for HOE is expected to match various optical fields with multifunctional diffractive optical elements, such as the holographic projection and the multi-depth structured light generation.

Funding

National Natural Science Foundation of China (61975014, 62035003); Beijing Municipal Science and Technology Commission, Administrative Commission of Zhongguancun Science Park (Z211100004821012).

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. L. L. Doskolovich, N. L. Kazanskiy, V. A. Soifer, P. Perlo, and P. Repetto, “Design of DOEs for wavelength division and focusing,” J Mod Opt 52(6), 917–926 (2005). [CrossRef]  

2. A. Hermerschmidt, S. Kruger, and G. Wernicke, “Binary diffractive beam splitters with arbitrary diffraction angles,” Opt. Lett. 32(5), 448–450 (2007). [CrossRef]  

3. S. R. Kitchen, C. Dam-Hansen, and M. L. Jakobsen, “Quasi-achromatic laser Doppler anemometry systems based on a diffractive beam splitter,” Appl. Opt. 42(28), 5642–5648 (2003). [CrossRef]  

4. R. de Saint Denis, N. Passilly, M. Laroche, T. Mohammed-Brahim, and K. Ait-Ameur, “Beam-shaping longitudinal range of a binary diffractive optical element,” Appl. Opt. 45(31), 8136–8141 (2006). [CrossRef]  

5. C. Yu, Y. Peng, T. Guo, Q. Fu, H. Li, and X. Liu, “Beam shaping for multicolour light-emitting diodes with diffractive optical elements,” J Mod Opt 64(4), 388–395 (2017). [CrossRef]  

6. H. Sasaki, K. Kotani, H. Wada, T. Takamori, and T. Ushikubo, “Scalability analysis of diffractive optical element-based free-space photonic circuits for interoptoelectronic chip interconnections,” Appl. Opt. 40(11), 1843–1855 (2001). [CrossRef]  

7. C. Jang, K. Bang, S. Moon, J. Kim, S. Lee, and B. Lee, “Retinal 3D: Augmented Reality Near-Eye Display Via Pupil-Tracked Light Field Projection on Retina,” ACM Trans. Graph. 36(6), 1–13 (2017). [CrossRef]  

8. J. Zhao, B. Chrysler, and R. K. Kostuk, “Design of a high-resolution holographic waveguide eye-tracking system operating in near-infrared with conventional optical elements,” Opt. Express 29(15), 24536–24551 (2021). [CrossRef]  

9. R. K. Kostuk, J. Zhao, and B. Chrysler, “Design of a waveguide eye-tracking system operating in near-infrared with holographic optical elements,” Opt. Eng. 60(08), 085101 (2021). [CrossRef]  

10. L. Gong, Z. Zhu, X. Wang, Y. Li, M. Wang, and S. Nie, “Changeable focused field distribution of double-ring-shaped cylindrical vector beams,” Opt. Commun. 342(1), 204–213 (2015). [CrossRef]  

11. M. A. Rykov and R. V. Skidanov, “Modifying the laser beam intensity distribution for obtaining improved strength characteristics of an optical trap,” Appl. Opt. 53(2), 156–164 (2014). [CrossRef]  

12. Y. Su, Z. Cai, W. Zou, L. Shi, F. Zhou, P. Guo, Y. Lu, and J. Wu, “Viewing angle enlargement in holographic augmented reality using an off-axis holographic lens,” Optik 172(11), 462–469 (2018). [CrossRef]  

13. K. Wakunami, P. Y. Hsieh, R. Oi, T. Senoh, H. Sasaki, Y. Ichihashi, M. Okui, Y. P. Huang, and K. Yamamoto, “Projection-type see-through holographic three-dimensional display,” Nat. Commun. 7(1), 12954 (2016). [CrossRef]  

14. M. H. Choi, K. S. Shin, J. Jang, W. Han, and J. H. Park, “Waveguide-type Maxwellian near-eye display using a pin-mirror holographic optical element array,” Opt. Lett. 47(2), 405–408 (2022). [CrossRef]  

15. L. Hwang, G. Hur, J. Kim, P. Gentet, S. Kwon, and S. Lee, “Uniformity improvement of a reconstructed-holographic image in a near-eye display system using off-axis HOE,” Opt. Express 30(12), 21439–21454 (2022). [CrossRef]  

16. K. S. Shin, M. H. Choi, J. Jang, and J. H. Park, “Waveguide-type see-through dual focus near-eye display with a polarization grating,” Opt. Express 29(24), 40294–40309 (2021). [CrossRef]  

17. K. Wu, H. Zhang, Y. Chen, Q. Luo, and K. Xu, “All-Silicon Microdisplay Using Efficient Hot-Carrier Electroluminescence in Standard 0.18 µ m CMOS Technology,” IEEE Electron Device Lett. 42(4), 541–544 (2021). [CrossRef]  

18. I. Sieber, R. Thelen, and U. Gengenbach, “Assessment of high-resolution 3D printed optics for the use case of rotation optics,” Opt. Express 28(9), 13423–13431 (2020). [CrossRef]  

19. Y. Luo, C. Chu, S. Vyas, H. Kuo, Y. Chia, M. Chen, X. Shi, T. Tanaka, H. Misawa, Y.-Y. Huang, and D. P. Tsai, “Varifocal Metalens for Optical Sectioning Fluorescence Microscopy,” Nano Lett. 21(12), 5133–5142 (2021). [CrossRef]  

20. M. Bawart, M. A. May, T. Oettl, C. Roider, S. Bernet, M. Schmidt, M. Ritsch-Marte, and A. Jesacher, “Diffractive tunable lens for remote focusing in high-NA optical systems,” Opt. Express 28(18), 26336–26347 (2020). [CrossRef]  

21. M. A. May, M. Bawart, M. Langeslag, S. Bernet, M. Kress, M. Ritsch-Marte, and A. Jesacher, “High-NA two-photon single cell imaging with remote focusing using a diffractive tunable lens,” Biomed. Opt. Express 11(12), 7183–7191 (2020). [CrossRef]  

22. A. Leal-Junior, L. Avellar, A. Frizera, and C. Marques, “Smart textiles for multimodal wearable sensing using highly stretchable multiplexed optical fiber system,” Sci. Rep 10(1), 13867 (2020). [CrossRef]  

23. K. Xu, “Silicon electro-optic micro-modulator fabricated in standard CMOS technology as components for all silicon monolithic integrated optoelectronic systems,” J. Micromech. Microeng. 31(5), 054001 (2021). [CrossRef]  

24. S. G. Ghebjagh, A. Behrens, P. Fesser, and S. Sinzinger, “Rotationally tunable multi-focal diffractive moire lenses,” Appl. Opt. 60(17), 5145–5152 (2021). [CrossRef]  

25. N. Bregenzer, T. Oettl, M. Zobernig, M. Bawart, S. Bernet, and M. Ritsch-Marte, “Demonstration of a multi-color diffractive lens with adjustable focal length,” Opt. Express 28(20), 30150–30163 (2020). [CrossRef]  

26. A. Grewe, P. Fesser, and S. Sinzinger, “Diffractive array optics tuned by rotation,” Appl. Opt. 56(1), A89–A96 (2017). [CrossRef]  

27. S. Bernet, W. Harm, and M. Ritsch-Marte, “Demonstration of focus-tunable diffractive Moire-lenses,” Opt. Express 21(6), 6955–6966 (2013). [CrossRef]  

28. S. Bernet and M. Ritsch-Marte, “Adjustable refractive power from diffractive moiré elements,” Appl. Opt. 47(21), 3722–3730 (2008). [CrossRef]  

29. P. Lin, Y. Lin, J. Lin, and B. Yang, “Stretchable metalens with tunable focal length and achromatic characteristics,” Results Phys. 31(1), 105005 (2021). [CrossRef]  

30. D. Wang, C. Liu, N. Li, and Q. Wang, “Holographic Zoom System With Large Focal Depth Based on Adjustable Lens,” IEEE Access 8, 85784–85792 (2020). [CrossRef]  

31. C. Kim, D. Shin, G. Koo, and Y. H. Won, “Fabrication of an electrowetting liquid microlens array for a focus tunable integral imaging system,” Opt. Lett. 45(2), 511–514 (2020). [CrossRef]  

32. C. Liu, D. Wang, and Q. Wang, “Holographic display system with adjustable viewing angle based on multi-focus optofluidic lens,” Opt. Express 27(13), 18210–18221 (2019). [CrossRef]  

33. Y. Wang and Y. Lin, “An optical system for augmented reality with electrically tunable optical zoom function and image registration exploiting liquid crystal lenses,” Opt. Express 27(15), 21163–21172 (2019). [CrossRef]  

34. H. Chen, Y. Wang, P. Chen, and Y. Lin, “Electrically adjustable location of a projected image in augmented reality via a liquid-crystal lens,” Opt. Express 23(22), 28154–28162 (2015). [CrossRef]  

35. J. Atencia and M. Quintanilla, “Ray tracing for holographic optical element recording with non-spherical waves,” J. Opt. A: Pure Appl. Opt. 3(5), 387–397 (2001). [CrossRef]  

36. S. Lee, B. Lee, J. Cho, C. Jang, J. Kim, and B. Lee, “Analysis and Implementation of Hologram Lenses for See-Through Head-Mounted Display,” IEEE Photonics Technol. Lett. 29(1), 82–85 (2017). [CrossRef]  

37. W.-K. Lin, O. Matoba, B.-S. Lin, and W.-C. Su, “Astigmatism and deformation correction for a holographic head-mounted display with a wedge-shaped holographic waveguide,” Appl. Opt. 57(25), 7094–7101 (2018). [CrossRef]  

38. W.-K. Lin, O. Matoba, B.-S. Lin, and W.-C. Su, “Astigmatism correction and quality optimization of computer-generated holograms for holographic waveguide displays,” Opt. Express 28(4), 5519–5527 (2020). [CrossRef]  

39. H.-J. Yeom, H.-J. Kim, S.-B. Kim, H. Zhang, B. Li, Y.-M. Ji, S.-H. Kim, and J.-H. Park, “3D holographic head mounted display using holographic optical elements with astigmatism aberration compensation,” Opt. Express 23(25), 32025–32034 (2015). [CrossRef]  

40. T. Lloret, V. Navarro-Fuster, M. G. Ramírez, M. Morales-Vidal, A. Beléndez, and I. Pascual, “Aberration-Based Quality Metrics in Holographic Lenses,” Polymers 12(4), 993 (2020). [CrossRef]  

41. Z. Lv, J. Liu, and L. Xu, “A Multi-Plane Augmented Reality Head-Up Display System Based on Volume Holographic Optical Elements With Large Area,” IEEE Photonics J. 13(5), 1–8 (2021). [CrossRef]  

42. M. Yu, Optical holography and information processing (National Defense Industry Press, 1984), Chap.2.

43. W. J. Smith, Modern Optical Engineering (Chemical Industry Press, 2011), Chap. 5.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (8)

Fig. 1.
Fig. 1. (a) shows the reconstruction process of HOE and (b) illustrates the process of focal length adjustment and aberration correction by SLM. $IW$: the illumination wave, $RI$: the reconstructed image, ${l_d}$: the designed focal length, $MF$: the phase distribution of modulation factor, $SLM$: spatial light modulator, ${I_1}{I_3}$: reconstructed images at different distances.
Fig. 2.
Fig. 2. Flow chart for determining correction parameters.
Fig. 3.
Fig. 3. The oblique incidence of illumination wave into HOE introduces astigmatism. (b) is the top view of (a).
Fig. 4.
Fig. 4. Phase difference caused by oblique incidence of illumination light
Fig. 5.
Fig. 5. The actual radius of curvature on HOE varies with the radius of curvature of the illumination light with the influence of oblique incidence angle.
Fig. 6.
Fig. 6. Aberration correction at the designed focal length. (a) Original image. (b) the reconstructed image before correction. (c) the reconstructed image after zero-order light elimination. (d) the original image with pre-distortion. (e) - (g) represents the reconstructed image after correcting distortion, astigmatism and field curvature respectively. (h) - (l) is a partial enlarged view of (e) - (g).
Fig. 7.
Fig. 7. Parameters for compensating phase distribution
Fig. 8.
Fig. 8. Reconstructed images at different distances after aberration compensation. The reconstructed images before and after compensation are shown in (a) - (c) and (d) - (f) respectively, where (a) (d), (b) (e) and (c) (f) correspond to the reconstructed images adjusted to 30 cm, 40 cm and 50 cm respectively.

Equations (11)

Equations on this page are rendered with MathJax. Learn more.

$$\scalebox{0.88}{$\displaystyle E({x^{\prime},y^{\prime}} )= {c_1}{c_2}|R ||C |\cdot \mathrm{{\cal F}}\left\{ {exp\left[ {ik\left( { - \frac{1}{{{l_d}}} - \frac{1}{{{l_c}}} + \frac{1}{{{l_i}}}} \right) \cdot \frac{{{x^2} + {y^2}}}{2}} \right]} \right\} \otimes \mathrm{{\cal F}}\left\{ {{\mathrm{{\cal F}}^{ - 1}}\left\{ {o({{x_1},{y_1}} )exp\left( { - ik\frac{{{x_1}^2 + {y_1}^2}}{{2{l_d}}}} \right)} \right\}} \right\},$}$$
$$\frac{1}{{{l_i}}} = \frac{1}{{{l_c}}} + \frac{1}{{{l_d}}}.$$
$$\begin{array}{l} \left\{ \begin{array}{l} {x_d} = x({1 + {k_1}{r^2} + {k_2}{r^4} + {k_3}{r^6}} )\\ {y_d} = y({1 + {k_1}{r^2} + {k_2}{r^4} + {k_3}{r^6}} )\end{array} \right.\\ \left\{ \begin{array}{l} {x_d} = x + [{2{p_1}y + {p_2}({{r^2} + 2{x^2}} )} ]\\ {y_d} = y + [{2{p_2}x + {p_1}({{r^2} + 2{y^2}} )} ]\end{array} \right., \end{array}$$
$${E_{dist}} = \sum\limits_i {\sqrt {{{({{x_i} - {x_c}} )}^2} + {{({{y_i} - {y_c}} )}^2}} } ,$$
$$exp\left[ {\frac{{i\pi }}{\lambda }\left( {\frac{{{x^2}}}{{{l_{ax}}}} + \frac{{{y^2}}}{{{l_{ay}}}}} \right)} \right] \cdot exp\left[ {\frac{{i\pi }}{{\lambda {l_v}}}({{x^2} + {y^2}} )} \right].$$
$${E_{asti}} = 1 - \frac{{|{\overline {{h_x}} - \overline {{h_y}} } |}}{{\overline {{h_x}} + \overline {{h_y}} }}.$$
$${E_{curv}} = 1 - \frac{{|{\overline {{h_{ex}}} - \overline {{h_{cx}}} } |}}{{\overline {{h_{ex}}} + \overline {{h_{cx}}} }}.$$
$$\frac{1}{{{l_i}}} = \frac{1}{{{l_d}}} + \frac{1}{{{l_{dst(m )}}}} + \frac{1}{{{l_{dcv(m )}}}} + \frac{1}{{{l_{c(m )}}}}.$$
$${f_{(m )}} = {1 / {\left( {\frac{1}{{{l_d}}} + \frac{1}{{{l_{dst(m )}}}} + \frac{1}{{{l_{dcv(m )}}}} + \frac{1}{{{l_{c(m )}}}}} \right)}}.$$
$$\left\{ \begin{array}{l} {\varphi_{ts}} ={-} k\sqrt {{{({{y_h} + {l_{c - SLM}}sin\theta } )}^2} + {{({{l_{c - SLM}}cos\theta } )}^2}} \\ {\varphi_{tp}} = k\left( {\frac{{{d_{hy}}}}{2} - {y_h}} \right)sin\theta \\ {\varphi_s} ={-} k\sqrt {{x_h}^2 + {y_h}^2 + {l_{c - SLM}}^2} \\ {\varphi_d} = {\varphi_{ts}} - {\varphi_{tp}} - {\varphi_s} \end{array} \right..$$
$${l_{c - HOE}} = 1.22 \times {l_{c - SLM}} - 1.35 \times {10^{ - 4}}.$$
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.