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

Design of adjustable multifocal diffractive optical elements with an improved smooth phase profile by continuous variable curve with multi-subperiods method

Open Access Open Access

Abstract

Multifocal diffractive optical elements (MDOEs), which produce arbitrary light distribution, are widely used in lightweight and compact optical systems. MDOEs that are combined with multiple functions tend to have complex step structures, limiting their applications. We propose a facile method named continuous variable curve with multi-subperiods (CVCMS) to design adjustable multifocal single-layer diffractive optical elements. Through the analysis, the model achieved arbitrary diffraction efficiency distribution with an improved smooth continuous phase profile in each diffractive ring while retaining the periodicity. To display the high design freedom of the method, we utilized this method to design and discuss a broadband multifocal intraocular lens (MIOL) focused on the optimization of far focal point. Finally, the method was compared with other multifocal design methods. The results show that the CVCMS method achieved adjustable multifocal design with better performance and smoother profile than other MDOE design techniques. The proposed model can be applied to multifocal ophthalmic lens designs.

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

1. Introduction

Diffractive optical elements (DOEs) produce the ideal diffraction efficiency distribution by generating phase change [1]. Whereas refractive optical elements are too bulky to be used in lightweight and compact systems, DOEs—with their arbitrary phase distribution [2], thin element property [3], and special dispersion ability [4,5]—are applied in infrared systems [6], microscopes systems [7], and vision systems [810].

The phase distribution of DOEs causes light beams to split into multiple diffraction orders. Based on this performance, Multifocal diffractive optical elements (MDOEs) can produce different focal lengths, thus achieving multifocal imaging. MDOEs are commonly utilized in the design of ophthalmic lenses, especially intraocular lenses (IOLs) [1114]. Conventional MDOE consist of superimposed diffractive structures with different period radii from the inside to the outside of the lens, which yield a varying diopter and thus produce a multifocal image [1517].

However, the structural design of MDOEs is responsible for imaging defects because each diffraction ring can only take part of the imaging effect when an MDOE with a combination of different diffractive structures is applied. The transition region leads to stray light and to interference between adjacent focal points, causing image quality to deteriorate. The imaging in this type of design is limited by the size of the aperture: when the aperture of the system is reduced, this type of design will not make all the diffractive regions effective [17]. For example, in diffractive IOL designs, when the pupil contracts in a bright environment, the edge diffractive structures do not participate in imaging, which may cause incomplete realization of the design function [18]. When the pupil is dilated under low-light conditions, more diffractive regions with different designs direct light to the retina simultaneously; therefore, diffractive multifocal lenses have the potential to produce halos and glare [19]. They also lead to a decrease in contrast sensitivity [20]. Overall, a method to design high-diffraction-efficiency MDOEs without destroying the diffraction periodicity must be discovered.

Furthermore, owing to their superior performance and thin element property, MDOEs can be designed with other functions, such as broadband. Typically, the DOE method increases the design freedom required to realize broadband by designing multilayer diffractive optical elements. We previously studied the imaging characteristics of broadband multilayer diffractive optics elements and realized high-diffraction-efficiency imaging with broadband [21]. However, for applications such as ophthalmology systems, multilayer MDOEs cannot be used because of the limited design space. Hence, it is necessary to create a model to design single-layer DOEs with a high degree of freedom to achieve functions such as broadband imaging. Peng et al. [22] used the modified particle swarm optimization algorithm to identify the diffraction phase structure of the max point spread function (PSF) and achieve high-quality broadband imaging. Kim et al. [23] proposed to divide the diffraction plane into M × N small regions. This approach provided high-quality broadband imaging based on the direct binary search (DBS) method using the Fresnel diffraction principle. Doskolovich et al. [24] proposed a diffraction design that generates multiple fixed focal points at different wavelengths, by dividing the diffraction plane into N small regions. Each region used an algorithm to obtain the optimal solution, thus, the design model of broadband imaging was obtained. Overall, most studies achieved excellent broadband imaging on single-layer diffractive optical elements through direct laser writing technique either with complex step structures changing in one period or with dividing small regions discarded periodicity. However, for typical applications of single-layer MDOEs such as ophthalmic lenses, the most common manufacturing method is still single point diamond machining. Unfortunately, as we illustrated in another study, step DOE processing with single-point diamond has high requirements in terms of tool type and radius; further limitations arise from the additional shadow effects caused by the diffractive step structure [25]. Therefore, processing of the complex diffractive step structures with single point diamond is extremely difficult. A smooth continuous diffraction structure can be the key to improve this situation. Moreover, a smooth curved design is also more advantageous for lenses implanted in the human eye such as IOL. Owing to there are no visible protrusions which may cause postoperative adverse effects, the design is more compatible with the human eyes [26]. In addition, the accumulation of removed lens epithelial cells (LECs) and inflammatory cells is the main cause of recurrence after IOL surgery [27]. A smoother curved phase profile is less likely to accumulate impurities in the grooves.

Considering the above views, a facile method of MDOE design with high freedom and a smooth profile while maintaining the periodicity would be desirable. Accordingly, we proposed the continuous variable curve with multi-subperiods (CVCMS) method. The detailed design process of broadband MDOE systems using CVCMS method is proposed and discussed through model and parametric analyses. The proposed method offers significantly better freedom of design, high diffraction efficiency, excellent broadband performance, and smooth diffraction profile. The method can also be used for realizing other functions that are difficult to achieve with single-layer MDOEs. The remainder of the paper is structured as follows. Section 2 analyzes and discusses the scalar diffraction theory (SDT) to establish the basis of the study. Section 3 details the complete design process of the CVCMS method based on the SDT and analyzes the method. Section 4 derives the evaluation function of the diffraction efficiency for broadband MDOEs. Subsequently, a multifocal intraocular lens with broadband capability is designed, and the MDOE produced using the CVCMS method is compared with conventional multifocal designs. Section 5 summarizes the study.

2. Principle of diffractive optical elements (DOEs)

In this section, we analyzed the SDT as the basis for the rest of the study from the perspective of wavefront phase transformation. Figure 1(a) illustrates the wavefront transformation of DOE. The wavefront determines the focal length, which is expressed as:

$$f = \frac{{{r^2}}}{{2j{\lambda _0}}},$$
where r is the period radius, λ0 is the design wavelength, and j is the diffraction period number.

 figure: Fig. 1.

Fig. 1. Schematic of diffractive optical elements (DOEs): (a) Wavefront of DOEs; (b) actual profile of DOEs; (c) optical path diagram; (d) diffraction efficiency of monofocal DOEs with different harmonic numbers at different wavelengths; (e) diffraction efficiencies of bifocal DOEs with different phase delays at different wavelengths; (f) diffraction efficiency of an example multifocal DOE at different wavelengths.

Download Full Size | PDF

The diffractive optical element achieves an arbitrary light distribution through a period of phase change. The phase profile is expressed as:

$$\varphi = 2\pi \mu \beta \gamma , $$
where $\gamma = j - \frac{{{r^2}}}{{2{\lambda _0}f}}$ is the phase curve, µ is the dispersion constant, β is the phase delay which determines the diffraction step height. Step height (Fig. 1(b)) is expressed as:
$$h = \frac{{{\lambda _0}\beta }}{{n^{\prime} - n}}, $$
where h is the step height, n’ and n are the refractive indices of the lens and medium, respectively. The dispersion constant µ is determined by the dispersion coefficient, which is expressed as:
$$l = \frac{{{\varphi _{{\lambda _{\max }}}} - {\varphi _{{\lambda _{\min }}}}}}{{{\varphi _{{\lambda _0}}}}}, $$
where ${\varphi _{{\lambda _{\max }}}}$ and ${\varphi _{{\lambda _{\min }}}}$ are the phases at maximum and minimum wavelengths, respectively, and ${\varphi _{{\lambda _0}}}$ is the phase at the design wavelength. The inverse of the dispersion coefficient is called the Abbe number. The two quantities are generally used to evaluate the dispersion ability of a DOE structure. The dispersion constant at any point in this interval can be written as:
$$\mu = \frac{{{\varphi _\lambda }}}{{{\varphi _{{\lambda _0}}}}} = \frac{{{\lambda _0}({{{n^{\prime}}_\lambda } - {n_\lambda }} )}}{{\lambda ({{{n^{\prime}}_{{\lambda_0}}} - {n_{{\lambda_0}}}} )}}, $$
where λ is the actual wavelength. ${n^{\prime}_\lambda }$ and ${n_\lambda }$ are the refractive indices of the lens and medium at the actual wavelength, respectively. ${n^{\prime}_{{\lambda _0}}}$ and ${n_{{\lambda _0}}}$ are the refractive indices of the lens and medium at the design wavelength, respectively.

The diffraction efficiency can be derived from the Fourier expansion of the lens transmittance function as follows:

$${c_m} = \frac{1}{{2{\lambda _0}f}}\int_0^{2{\lambda _0}f} {\exp ( - \frac{{{\lambda _0}({{{n^{\prime}}_\lambda } - {n_\lambda }} )}}{{\lambda ({{{n^{\prime}}_{{\lambda_0}}} - {n_{{\lambda_0}}}} )}}\frac{{i\pi \beta {r^2}}}{{{\lambda _0}f}})\exp (\frac{{im\pi {r^2}}}{{{\lambda _0}f}})} d{r^2}. $$

Then, the diffraction efficiency is expressed as:

$${\eta _m} = {|{{c_m}} |^2} = {\left|{\frac{1}{{2{\lambda_0}f}}\int_0^{2{\lambda_0}f} {\exp ( - \frac{{{\lambda_0}({{{n^{\prime}}_\lambda } - {n_\lambda }} )}}{{\lambda ({{{n^{\prime}}_{{\lambda_0}}} - {n_{{\lambda_0}}}} )}}\frac{{i\pi \beta {r^2}}}{{{\lambda_0}f}})\exp (\frac{{im\pi {r^2}}}{{{\lambda_0}f}})} d{r^2}} \right|^2}. $$

The variation of the refractive index of most diffractive optical element materials with the wavelength is slow, i.e., $\frac{{({{{n^{\prime}}_\lambda } - {n_\lambda }} )}}{{({{{n^{\prime}}_{{\lambda_0}}} - {n_{{\lambda_0}}}} )}} \to 1$. Thus, the effect of the material can be neglected. Dispersion constant can be simplified as $\mu = \frac{{{\lambda _0}}}{\lambda }$. Then, Eq. (7) is rewritten as:

$${\eta _m} = {\left|{\frac{1}{{2{\lambda_0}f}}\int_0^{2{\lambda_0}f} {\exp ( - \frac{{i\pi \mu \beta {r^2}}}{{{\lambda_0}f}})\exp (\frac{{im\pi {r^2}}}{{{\lambda_0}f}})} d{r^2}} \right|^2}. $$

As depicted in Fig. 1(c), different orders of light with different diffraction efficiencies are directed to their respective focal points.

For the conventional DOE design, the phase distribution is determined by the phase delay β. Different phase delays correspond to different light distribution (Figs. 1(d)–(f)). The diffraction efficiency at different orders changes with β. In Fig. 1(d), when β is 1, the diffraction efficiency at the design wavelength reaches 1. When the phase delay is set to β+M, at this time the phase $\varphi = 2\pi (\beta + M)\gamma $, where M is called the harmonic number. Accordingly, such DOEs are known as harmonic diffractive optical elements (HDOEs), which have a smaller dispersion range and are typically used in optical systems to correct chromatic aberration. Several studies have used HDOEs to design diffractive lenses that performed better than their conventional counterparts. However, as Fig. 1(d) shows, the use of HDOEs also poses problems for imaging over a wide range of wavelength. When the value β varies between 0 and 1, the distribution of light changes between two diffraction orders. When β is 0.5, the distribution between diffraction orders 0 and 1 is 1:1, which is the common bifocal DOE design with equal diffraction efficiency. As indicated in Fig. 1(e), the broadband nature of the bifocal design is degraded compared to the monofocal design, leading to poor imaging over the broadband spectrum. Figure 1(f) illustrates a trifocal design with design ideas from Ref. [28]. However, as Fig. 1(f) shown, this design is inadequate for imaging over the broadband, and will cause degradation of imaging quality when the wavelength change. Considering Fig. 1(d)-(f), the broadband performance of DOEs declines as more focal points are applied. Thus, we have described and analyzed monofocal, bifocal and multifocal diffractive optical elements under SDT. Based on these principles, we designed the CVCMS method in Section 3 and a specific system using the CVCMS method in Section 4.

3. Construction and analysis of CVCMS method

In this section, we detail the design method, CVCMS, for adding degrees of freedom and smoothening the phase profile of single-layer DOEs. As depicted in Fig. 2, the construction of the model is divided into three steps. First, the diffraction phase-profile model of a single-layer DOE must be improved. An efficient initial structural model is essential for the subsequent optimization of the structure. Therefore, in Step 1, we transformed the diffraction phase profile to smoothen it and make the profile parameters more favorable for subsequent derivations.

 figure: Fig. 2.

Fig. 2. Continuous variable curve with multi-subperiods (CVCMS) design process.

Download Full Size | PDF

Second, we performed a regional division within the same period. As described in Section 1, multifocal design based on the superposition of different DOE structures causes degradation of the image quality. Therefore, in Step 2, we establish a subperiods optimization method for the structure proposed in Step 1 to improve the freedom of design while maintaining periodicity.

Finally, for a diffractive structure, as noted in Section 1, a continuous smooth phase profile in each period is beneficial both in manufacturing and applications. Hence, in Step 3, we propose a continuous curve design method to smoothen the curve obtained in Step 2 and make the connections of each subperiod more continuous.

In summary, in Section 3.1, we establish a method to change the diffraction phase profile. In Section 3.2, we provide a method to create multiple subperiods within the same period. In Section 3.3, we provide a method to make phase curve continuous and smooth, and demonstrate the CVCMS method. The entire design model is constructed and analyzed in three steps.

3.1 Step 1-Change of diffraction phase profile

From Section 2, the Fourier series of transmittance function of DOEs can be expressed as:

$${c_m} = \frac{1}{{2{\lambda _0}f}}\int_0^{2{\lambda _0}f} {\exp (\frac{{ - i\pi \mu \beta {r^2}}}{{{\lambda _0}f}})\exp (\frac{{m\pi {r^2}}}{{{\lambda _0}f}})} d{r^2}. $$

The choice of the initial structure is crucial. At this point, as observed in Eq. (9), the diffraction efficiency is determined by the maximum phase delay and phase profile. The phase profile of a conventional diffractive structure is a fixed step curve, we transformed the profile curve into a complex function using more operands. The first step is to transform the conventional step structure curve into a Sin/Cos-type function, which has a period of 2π and a more curved phase profile. It is consistent with the diffraction periodicity and easily accepts more degrees of freedom. In summary, the Fourier series for the change curve is expressed as:

$${c_m} = \frac{1}{{2\pi }}\int_0^{2\pi } {\exp ( - i\mu \beta \cos (\theta ))\exp (im\theta )} d\theta. $$

However, for this model, the problems stemming from the limitation of the structural degrees of freedom must be solved. the equal energy distribution in the –1 and +1 orders by Eq. (10) leads to a restricted design; meanwhile, when we attempted to increase the diffraction efficiency at ±1 orders, the diffraction efficiencies at ±1 orders could only be increased to approximately 0.338, and lead to a significant decrease of diffraction efficiency at 0 order. For example, if β = 1.7, the diffraction efficiency of the ±1 order reaches the almost maximum designable value. However, the diffraction efficiency at the 0th order is only 0.15 and the broadband performance is poor. Hence, it is also inconvenient to design a DOE with a high diffraction efficiency at an arbitrary order.

3.2 Step 2-Method of multi-subperiods

As explained in Section 1, the diffraction efficiency of traditional multifocal DOEs is limited by the size of the applied diffraction area. Therefore, only structural changes within a period can achieve consistent imaging effects from inside to outside. Accordingly, we improved the form of the phase expression within the period mentioned in Section 3.1 to maintain periodicity. We divided the whole period into k weighted regions. w represents the design proportion of each subperiod. The selection of w depends on the design and solution requirements. cm can be expressed as:

$${c_m} = \frac{1}{{2\pi }}\left[ \begin{array}{l} \int_0^{2{w_1}\pi } {\exp ( - i\mu {\beta_1}\cos (\theta ))\exp (im\theta )} d\theta + \int_{2{w_1}\pi }^{2{w_2}\pi } {\exp ( - i\mu {\beta_2}\cos (\theta ))\exp (im\theta )} d\theta \\ + \int_{2{w_2}\pi }^{2{w_3}\pi } {\exp ( - i\mu {\beta_3}\cos (\theta ))\exp (im\theta )} d\theta + \cdots \cdots \\ + \int_{2{w_{k - 1}}\pi }^{2{w_k}\pi } {\exp ( - i\mu {\beta_k}\cos (\theta ))\exp (im\theta )} d\theta \end{array} \right], $$
where $0 \le {w_1} \le {w_2} \le {w_3} \cdots \le 1$, and ${w_k} = 1$.

Because the diffraction efficiency is expressed as ${|{{c_m}} |^2}$, in Eq. (11), the subperiods are not independent and all affect each other. Moreover, no regions can be arbitrarily adjusted, which means that the sequential order of adjustment has an impact on the diffraction efficiency.

3.3 Step 3-Method of continuous curve

Application of the method derived in Section 3.2 would inevitably cause discontinuities between subperiods. Combining the analysis of the above two steps, we proposed a continuous curve approach. Under Cos transformation, we defined tk as the continuous factor and added the corresponding tk to each subperiod in Eq. (11), obtaining Eq. (12). After change, the maximum and minimum values of Cos remained ±1, showing that tk has no impact on the maximum phase delay of each subperiod in Eq. (11). Therefore, the continuous-curve design can be achieved while still supplementing more operands. cm that with continuous factor tk can be expressed as:

$$\scalebox{0.96}{${c_m} = \frac{1}{{2\pi }}\left[ \begin{array}{@{}l@{}} \int_0^{2{w_1}\pi } {\exp ( - i\mu {\beta_1}\cos (\theta + {t_1}))\exp (im\theta )} d\theta + \int_{2{w_1}\pi }^{2{w_2}\pi } {\exp ( - i\mu {\beta_2}\cos (\theta + {t_2}))\exp (im\theta )} d\theta \\ + \int_{2{w_2}\pi }^{2{w_3}\pi } {\exp ( - i\mu {\beta_3}\cos (\theta + {t_3}))\exp (im\theta )} d\theta + \cdots \cdots \\ + \int_{2{w_{k - 1}}\pi }^{2{w_k}\pi } {\exp ( - i\mu {\beta_k}\cos (\theta + {t_k}))\exp (im\theta )} d\theta \end{array} \right].$}$$

The first period should preferably be set in the center, so that the initial point is well determined and processed, i.e., t1 = 0. We make the diffraction phase profile have equal values at each subperiod edge. Hence,

$$\left\{ {\begin{array}{{c}} {{\beta_1}\cos (2\pi {w_1}) = {\beta_2}\cos (2\pi {w_1} + {t_2})}\\ {{\beta_2}\cos (2\pi {w_2} + {t_2}) = {\beta_3}\cos (2\pi {w_2} + {t_3})}\\ {{\beta_{k - 1}}\cos (2\pi {w_{k - 1}} + {t_{k - 1}}) = {\beta_k}\cos (2\pi {w_{k - 1}} + {t_k})} \end{array}} \right.. $$

Then, t2, t3, …tkcan be expressed as:

$$\left\{ {\begin{array}{{c}} {{t_2} = \arccos \left( {\frac{{{\beta_1}}}{{{\beta_2}}}\cos ({2\pi {w_1}} )} \right) - 2\pi {w_1}}\\ {{t_3} = \arccos \left( {\frac{{{\beta_2}}}{{{\beta_3}}}\cos ({2\pi {w_2} + {t_2}} )} \right) - 2\pi {w_2}}\\ {{t_k} = \arccos \left( {\frac{{{\beta_{k - 1}}}}{{{\beta_k}}}\cos ({2\pi {w_{k - 1}} + {t_{k - 1}}} )} \right) - 2\pi {w_{k - 1}}} \end{array}} \right.. $$

To ensure the validity of the equation, we need to ensure that

$$- 1 \le \left( {\frac{{{\beta_{k - 1}}}}{{{\beta_k}}}\cos ({2\pi {w_{k - 1}} + {t_{k - 1}}} )} \right) \le 1. $$

Values outside this region are discarded. Substituting Eq. (14) into Eq. (12) yields

$$\begin{array}{l} {c_m} = \frac{1}{{2\pi }}\\ \left[ \begin{array}{l} \int_0^{2{w_1}\pi } {\exp ( - i\mu {\beta_1}\cos (\theta ))\exp (im\theta )} d\theta \\ + \int_{2{w_1}\pi }^{2{w_2}\pi } {\exp ( - i\mu {\beta_2}\cos (\theta + \arccos \left( {\frac{{{\beta_1}}}{{{\beta_2}}}\cos ({2\pi {w_1}} )} \right) - 2\pi {w_1}))\exp (im\theta )} d\theta \\ + \int_{2{w_2}\pi }^{2{w_3}\pi } {\exp ( - i\mu {\beta_3}\cos (\theta + \arccos \left( {\frac{{{\beta_2}}}{{{\beta_3}}}\cos ({2\pi {w_2} + {t_2}} )} \right) - 2\pi {w_2}))\exp (im\theta )} d\theta + \cdots \cdots \\ + \int_{2{w_{k - 1}}\pi }^{2{w_k}\pi } {\exp ( - i\mu {\beta_k}\cos (\theta + \arccos \left( {\frac{{{\beta_{k - 1}}}}{{{\beta_k}}}\cos ({2\pi {w_{k - 1}} + {t_{k - 1}}} )} \right) - 2\pi {w_{k - 1}}))\exp (im\theta )} d\theta \end{array} \right] \end{array}, $$
$${\eta _m} = {|{{c_m}} |^2}. $$

Thus, we completed the construction of the CVCMS method. To demonstrate the degrees of freedom of the model, we illustrate an example of two subperiods in Fig. 3. As indicated in Eq. (1), different orders correspond to different distances. We utilized this approach to design a multifocus IOL, which propose a distribution of the diffraction efficiencies at the far, intermediate, and near focal points with the orders of –1, 0, and 1 (The specific system design is explained in Section 4.2). We derived the relationship between the phase delays of two subperiods and the diffraction efficiency at subperiod division proportion w1of 0.3, 0.5, or 0.7. The three various w1 in Fig. 3 are analyzed to illustrate the diverse diffraction efficiency distributions at different subperiod division proportions. The design results, compared to the problems listed with Eq. (10), break the diffraction efficiency value of ±1 orders that must be equal, and lead to the arbitrary diffraction efficiency distribution at different orders. Unlike structures designed with conventional design, the structure designed with CVCMS has continuity, i.e., the curve remains continuous over a period (examples are available in Figs. 2 and 7). The brown region in Fig. 3 represents the area where Eq. (15) is not satisfied, i.e., the region where a continuous curve could not be designed. The diffraction efficiency distribution varies with the size of the diffraction subregion and phase delay distribution, which corroborates the assertion in Section 3.2 that the subperiods affect one another and that the sequence of each subperiod cannot be manipulated. By accepting different input values, the design model can achieve an adjustable phase distribution with the continuous curve. A further increase in the number of subperiods helps increase the degrees of freedom for further design. This will be discussed in the subsequent section.

 figure: Fig. 3.

Fig. 3. Results of CVCMS-designed MDOE between phase delays of two subperiods and diffraction efficiency. (a) Subperiod division proportion w1 = 0.3; (b) subperiod division proportion w1 = 0.5; (c) subperiod division proportion w1 = 0.7.

Download Full Size | PDF

4. Design example and discussion

In this section, we demonstrated a specific application of the proposed method. We first propose the evaluation method for broadband MDOEs. Next, a broadband multifocal IOL is designed based on the CVCMS method, and CVCMS method is discussed with other design methods.

4.1 Evaluation of diffraction efficiency with broadband multifocal design

Evaluation of broadband MDOEs differs from that of conventional single-focus diffractive optical elements in that excessive diffraction efficiency at a single order is considered detrimental because it affects the magnitude of the efficiency at other orders. MDOEs need to achieve a balanced high diffraction efficiency for broadband imaging. With the proposed design model, an arbitrary diffraction-efficiency distribution can be achieved by applying a simple loop solution. Nevertheless, the design results for the DOE should be evaluated. Owing to the large variation in diffraction efficiency, conventional evaluation methods like the polychromatic integral diffraction efficiency (PIDE) are not applicable to broadband MDOEs. Therefore, we propose using two equations for evaluation instead. First, the degree of separation is expressed as:

$$\sigma = \left|{\frac{{{\eta_\lambda } - {\eta_{{\lambda_0}}}}}{{{\eta_{{\lambda_0}}}}}} \right|, $$
where ${\eta _\lambda }$ is the diffraction efficiency at actual wavelength, ${\eta _{{\lambda _0}}}$ is the diffraction efficiency at design wavelength. Equation (18) relates the degree of separation to the diffraction efficiency at the design wavelength and the actual wavelength. By controlling the differences, we selected the appropriate structural parameters. We also provide a formula to evaluate the overall broadband stability of the diffraction efficiency in this wavelength range:
$$\begin{aligned} {{\bar{p}}_m} &= \left|{\frac{1}{{{\lambda_{\max }} - {\lambda_{\min }}}}\int_{{\lambda_{\min }}}^{{\lambda_{\max }}} {\frac{{{\eta_\lambda }}}{{{\eta_{{\lambda_{\max }}}} - {\eta_{{\lambda_{\min }}}}}}d\lambda } } \right|\\ &= \left|{\frac{1}{{{\lambda_{\max }} - {\lambda_{\min }}}}\int_{{\lambda_{\min }}}^{{\lambda_{\max }}} {\frac{{{{|{{c_m}_{_\lambda }} |}^2}}}{{{{|{{c_m}_{_{{\lambda_{\max }}}}} |}^2} - {{|{{c_m}_{_{{\lambda_{\min }}}}} |}^2}}}d\lambda } } \right|\end{aligned}, $$
where λmax, λmin are the maximum and minimum wavelengths, ${\eta _m}_{_{{\lambda _{\max }}}}$, ${\eta _m}_{_{{\lambda _{\min }}}}$, ${c_m}_{_{{\lambda _{\max }}}}$, ${c_m}_{_{{\lambda _{\min }}}}$ are the diffraction efficiency and the Fourier expansion functions at maximum and minimum wavelengths, ${c_m}_{_\lambda }$ is the Fourier expansion function at actual wavelength, we call ${\bar{p}_m}$ as the stability factor. Equation (19) evaluates the stability of the diffraction efficiency over the entire waveband. A larger ${\bar{p}_m}$ indicates no excessive fall rate or rise rate in the entire waveband (An excessive rise rate inevitably decreases the diffraction efficiency at other focal points). Hence, the evaluation function must also include:
$${\bar{p}_m} \to \max, $$
$$\sum\limits_N {\eta {{(\lambda )}_{{m_N}}}} \to \max \textrm{ }, $$
where N is the diffraction order utilized.

4.2 Evaluation of broadband diffractive optical system design

We designed a multifocal intraocular lens system based on the structure of the human eye using the Liou–Brennan model [29]. The Liou–Brennan model consists of a cornea, pupil, crystalline lens, and retina. In this model, the crystalline lens uses the gradient refractive index to better simulate the refractive index of the natural crystalline lens in the real human eye. The specific parameters related to this model are listed in Table 1. The human eye focuses on the far field, while focusing on the other distances depend on the modulation of the crystalline lens.

Tables Icon

Table 1. Liou–Brennan model of the human eyea

When the natural crystalline lens aged, it was removed and placed into the designed diffractive multifocal intraocular lenses. The initial design parameters of the IOL are listed in Table 2. After adding the trifocal MDOE using ±1 and 0 orders, we set the additional diffractive diopters to –1.66D, 0D, 1.66D at the far, intermediate, and near points. Then, the diffractive IOL has powers of 18.7D at the far point (similar to the crystalline lens in the Liou–Brennan model), 20.36D at the intermediate point, and 22.02D at the near point. Different design bands will cause different optimization effects. In this design, we analyzed the wavelength that the human eye normally accepts (500–600 nm). The design process is illustrated in Fig. 4.

 figure: Fig. 4.

Fig. 4. Intraocular lens design (a) Liou–Brennan model of the eye (adapted from Ref. [29], Fig. 4). (b) Designed optical path diagram (red: Far point, green: Intermediate point, blue: Near point). (c) MTF results at design wavelength.

Download Full Size | PDF

Tables Icon

Table 2. Intraocular lens (IOL) design parameters

The final optimization result is summarized in Table 3. The modulation transfer function (MTF) at 80lp/mm of this model reached 0.3295, 0.3705, and 0.3945 at the design wavelength and achieved broadband MTF at the intermediate point as well as the near point; the broadband performance was relatively poor at the far point but still within the acceptable range.

Tables Icon

Table 3. Modulation transfer function (MTF) design

4.3 Results and discussion

During the design of multifocal lenses, the allocation of multifocal points depends on the requirement of the application. For the design in this section, we aimed at optimizing the far point to demonstrate the arbitrary manipulability of the proposed method, which helps people whose work requires far-point vision. The adjustable diffraction efficiency at specific orders is useful for different design objectives. However, arbitrary phase distribution also poses certain design challenges.

In addition, for a diffractive optical system, the PSF is determined by the diffraction efficiency size and MTF size. For MDOEs, the PSF is expressed as:

$$PSF(m,\lambda ) = \eta (m,\lambda ) \cdot PS{F_{doe}}(m,\lambda ), $$
$$MTF = |{FFT(PSF(m,\lambda ))} |. $$

As indicated in Eqs. (22) and (23), the final PSF is related to the diffraction efficiency. Optical software such as Code V and Zemax OpticsStudio normally set the diffraction efficiency to 1, because diffraction efficiencies with different orders are not considered in conventional monofocal lens designs. However, for MDOEs, as indicated in Eq. (22), a reasonable allocation of multifocus diffraction efficiency can help improve the overall imaging quality. In the system design proposed in Section 4.2, the MTF of the far point decreased more than that of the other two focal points when the system deviated from the design wavelength. Therefore, for this system, the diffraction efficiency at the far point should be increased.

To this end, we describe the design of a three-subperiods MDOE in this section. Different waveband requirements lead to different solution results. In this design, we solved for the best broadband design over the 500-600 nm waveband. The broadband evaluation function Eqs. (18), (19), (20) and (21) are given in Section 4.1 as the solution conditions. The parameters were solved by a simple loop solution. We first determined three equal subperiods in the split ratios of 0.33, 0.66, and 1, respectively. We then set the diffraction efficiencies at the intermediate point (0 order) and near point (1 order) to be equal as possible while maintaining the maximum at the optimized far point (–1 order), which is impossible to design through Eq. (10), and maintaining the broadband performance at all three focal points. The result parameters are shown in Table 4. The diffraction efficiency distribution is presented in Fig. 5 and Table 6.

 figure: Fig. 5.

Fig. 5. Diffraction efficiency of MDOE design using the CVCMS method.

Download Full Size | PDF

Tables Icon

Table 4. MDOE parameter results with CVCMS design

From Fig. 5, we achieved the case of almost equal diffraction efficiency at the intermediate and near focal points, with a focus on the optimization of far focal point. Then, to demonstrate the superiority of the model, we compared the proposed design with other multifocal design methods in terms of the diffraction efficiency, broadband performance evaluation, and actual profile in the following analysis.

To validate the comparison, we attempted to design similar diffraction efficiencies at the far point using different design methods. The comparison group design A was based on the superposition of different diffraction structures designed from inside to outside, which are common [14,16,17]. We simulated the ratio of incident light energy on the central and edge areas, which was 1:1. The region for far- and near-point imaging at the center radius spanned 0–2.121 mm. The region for the far- and intermediate-point imaging at the edge radius spanned 2.121–3 mm. The specific parameters are listed in Table 5. The comparison group design B was solved through a step structure through Eq. (24) using the design method of Ref. [28]. The specific parameters related to this are also listed in Table 5.

$$\begin{aligned} {c_m} &= \frac{1}{{2{\lambda _0}f}}\sum\limits_k {\int_{2{\lambda _0}{w_{k - 1}}f}^{2{\lambda _0}{w_k}f} {\exp ( - i\frac{{\pi {\beta _k}{r^2}}}{{{\lambda _0}f}})\exp (i\frac{{m\pi {r^2}}}{{{\lambda _0}f}})d{r^2}} } \\ &\quad {w_0} \le {w_{k - 1}} \le {w_k} \cdots \le {w_T},{w_0} = 0,{w_T} = 1 \end{aligned}. $$

First, we compared design A with the proposed design. We programmed a distribution with a similar diffraction efficiency as the proposed design at the far point. We observed that, at the other two focal points, the diffraction efficiency decreased significantly with the proposed design, and the overall diffraction efficiency utilization rate also differed notably. For multifocal diffraction structures, on the one hand, the overall utilization rate of diffraction efficiency determines the light energy distribution at the utilized focal points. On the other hand, because the unutilized light leads to phenomena such as glare, which also lead to the degradation of imaging quality, the overall diffraction efficiency utilization rate must be improved. For broadband performance, according to Table 6, the ${\bar{p}_m}$ values at the near and far points of design A were significantly lower than those with the proposed design but at the intermediate points were marginally higher. Notably, as the diffraction efficiency of design A was lower than that in the proposed design at design wavelength, the actual effect of intermediate points was smaller. In summary, the CVCMS design and design A exhibited similar broadband performances at the intermediate point. At the other two points, the proposed design was significantly better than design A in terms of both diffraction efficiency and broadband performance. Furthermore, the size of the utilized diffraction area significantly affected the performance of design A: when the radius was smaller than a certain threshold (< 2.121 mm), the diffraction structure at the edge responsible for the intermediate focal point could not be utilized, i.e., the intermediate diffraction efficiency was 0. Moreover, the overall performance resembled that of a bifocal lens. Combined with the above, our design’s performance is superior to that of design A.

Tables Icon

Table 5. Parameters of comparative designs A and B

Second, design B was compared with the proposed design. A diffraction efficiency similar to that achieved by the proposed method was also designed at the far point. As indicated in Table 6, the near- and intermediate-point diffraction efficiencies of design B were smaller than ours, and the overall diffraction efficiency was lower. Moreover, the diffraction efficiency of design B decreased rapidly in the waveband. Especially, at the intermediate and near points, design B produced weak imaging results in the broadband range. The evaluation function (Eq. (19)) reveals that its far, intermediate, and near ${\bar{p}_m}$ values were small compared to those of the proposed design. Design B exemplifies a general problem with this class of design with weak broadband imaging, which is also illustrated in Fig. 1(f). This phenomenon is visually represented in Fig. 6 through an evaluation of the degree of separation. As illustrated in Fig. 6, the separation of the two designs at the far point focus was extremely similar. This was probably because we made design B to closely emulate our design on the far point at the beginning of the design. At the other two distances, the degree of separation of our design was much less than that of design B. Combined with the analysis of ${\bar{p}_m}$, the results demonstrate simultaneous stable broadband imaging at multifocal points achieved by the proposed design compared to design B.

Tables Icon

Table 6. Results with different design approaches

 figure: Fig. 6.

Fig. 6. Degree of separation comparison between design B and CVCMS method.

Download Full Size | PDF

We evaluated the profile optimization between design B and the proposed design. The phase profile using the proposed CVCMS model at design wavelength can be expressed as

$$\begin{aligned} &{\varphi _k} = {\beta _k}\cos (\theta + {t_k}),\\ &{w_0} \le {w_{k - 1}} \le {w_k} \cdots \le {w_T},2\pi {w_{k - 1}} \le \theta \le 2\pi {w_k},{w_0} = 0,{w_T} = 1 \end{aligned}$$
where ${\varphi _k}$ is the phase profile in the kth subperiod. T is the total number of subperiods. The step height using the proposed model can be derived from Eq. (25) as ${h_{cvcms}} = \frac{{{\lambda _0}({\textstyle{{{\varphi _k}} \over {2\pi }}})}}{{n^{\prime} - n}}$. According to the expression of ${h_{cvcms}}$, the feature size of the design diffractive optical element meets the requirement of single-point diamond turning manufacturing. As shown in Fig. 7(a), the traditional design led to abrupt changes in height in a single period and thus to a discontinuous and jagged structure. In another study, we expounded that excessive shading would create more process difficulties when single-point diamond turning was used [25]. As the design requirements increased, when the number of subregions increase, single point diamond machining will become challenging. As shown in Fig. 7(b), within one period, the proposed profile is a smooth continuous curve. It should be noted that our structure, with optimized performance, will have an up and down height variation between each period but we can control this by controlling the phase of the last subperiod if the design so requires. Owing to the smoothness of the curve and the high degree of freedom in the proposed model, our model can be a way to replace the conventional diffraction structure model in order to expect a more optimized design. Moreover, the advantage of the smooth curve will continue to be extended as the number of interval divisions continues to increase over a period. As mentioned in Section 1, this architecture has unique advantages for systems such as IOL, and its application is yet to be extended. In summary, our design surpasses design B in two respects: performance and smoothness.

 figure: Fig. 7.

Fig. 7. DOE phase profile comparison between design B and CVCMS method. (a) Result by design B. (b) Result by CVCMS method.

Download Full Size | PDF

Finally, we elaborate on the adjustable diffraction efficiency improvement on both the PSF at any given order and the systems imaging quality. For the same DOE system, first, the design of three uniformly irradiated focal points was selected using Eq. (10) with β = 1.4376 as the control group. We compared it with the proposed CVCMS design in term of the PSF at far points. As illustrated in Fig. 8, the PSF of the proposed design at each wavelength is improved than that of control group. The growth rates of the entire pixel area with the PSF value at the far focal points were 1.1940, 1.1887, 1.1854, and 1.1777 at 500, 530, 550, and 600 nm, respectively. Figure 8 yields PSF is enhanced by the diffraction efficiency. Combined with Eqs. (22) and (23), the enhanced PSF, i.e., the enhanced MTF alleviate the relatively poor imaging performance at far point which shows in Section 4.2, as well as the lack of consideration of diffraction efficiencies at different orders by optical software. Meanwhile, we reiterate that our model can achieve a high diffraction efficiency at a far point without affecting the performance of the other two focal points (as a counterexample, the poorer broadband performance at the other two focal points of design B causes a reduction in the broadband PSF at the other two focal points). As the proposed design has a relatively high diffraction efficiency at the other two focal points and the IOL system itself has a high MTF at them (Table 3), through design, the quality of imaging at three points is high enough to satisfy the design requirements. In other words, the imaging quality of system is enhanced by the modulation of the diffraction efficiency. This also demonstrates that designing MDOEs that can achieve an arbitrary diffraction efficiency distribution can help systems achieve higher imaging quality.

 figure: Fig. 8.

Fig. 8. Schematic of the PSF at far focal point of an actual diffraction system (normalized by the maximum value), PSF values from left to right are compared at 500, 530, 550, and 600 nm. The PSF is calculated using a Fourier transform with a trifocal uniform design (a) or the proposed CVCMS design (b). The optical system uses the design in Section 4.2.

Download Full Size | PDF

In summary, based on the multifocal diffraction efficiency distribution required by the actual application and the evaluation of imaging quality of a diffractive optical system, we determined the diffraction efficiency distribution needed to design MDOEs. Next, according to the CVCMS method, we built the specific model design of MDOE. Finally, the design results were evaluated on five criteria: diffraction efficiency value of each wavelength, degree of separation, broadband stability, overall diffraction efficiency value, and optimal profile. After the three steps, we established the overall design process of MDOEs, which exhibited better performance than the conventional ones, and obtained a multifocal IOL design with high broadband diffraction efficiency and optimized profile. The method has excellent generality and plays an exemplary role in the design prospects of MDOEs.

5. Conclusion

In this study, we proposed the CVCMS method for designing MDOEs with arbitrary diffraction efficiency distribution and a smooth phase profile in each diffractive ring while retaining the periodicity. The method was implemented on a broadband multifocal intraocular lens design for demonstration and analyzed. The study was conducted as follows:

  • (1) We implemented the CVCMS method in three steps. To identify the initial structure, we first changed the phase profile. Then, to enhance the degrees of freedom, we create a sub-region within the same period. Finally, to smoothen and simplify the diffractive profile, we optimized the phase profile.
  • (2) MDOEs with two subperiods using the CVCMS method were analyzed with different period divisions and phase delays. The model achieved arbitrary diffraction efficiency modulation with a smooth continuous curve while retaining its periodicity.
  • (3) MDOEs with multi-subperiods using the CVCMS method were designed and discussed. First, the evaluation function of the broadband diffraction efficiency for MDOEs was proposed. Then, a broadband multifocal IOL using CVCMS was designed for demonstration. Conventional MDOE designs were compared with the proposed design in terms of diffraction efficiency, broadband performance, and actual profile optimization. The proposed method can achieve adjustable phase distribution and better imaging performance with a smoother profile than other designs. Based on its well performance, large degrees of freedom, and smooth phase profile, the method can be applied to multifocal ophthalmic lens designs.

Funding

Natural Science Foundation of Jilin Province (20220101124JC).

Acknowledgments

Authors thank Prof. Changxi Xue for his academic advice.

Disclosures

The authors declare that there are no conflicts of interest related to the work in this paper.

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. D. C. O’Shea, T. J. Sulski, A. D. Kathman, and D. W. Prather, Diffractive Optics: Design, Fabrication, and Test (SPIE, 2003).

2. H. Liu, Z. W. Lu, Q. Sun, and H. Zhang, “Design of multiplexed phase diffractive optical elements for focal depth extension,” Opt. Express 18(12), 12798–12806 (2010). [CrossRef]  

3. Z. Y. Hu, T. Jiang, Z. N. Tian, L. G. Niu, J. W. Mao, Q. D. Chen, and H. B. Sun, “Broad-bandwidth micro-diffractive optical elements,” Laser Photonics Rev. 16(3), 2270015 (2022). [CrossRef]  

4. M. Singh, J. Tervo, and J. Turunen, “Broadband beam shaping with harmonic diffractive optics,” Opt. Express 22(19), 22680–22688 (2014). [CrossRef]  

5. M. S. Millán, F. Vega, and I. Ríos-López, “Polychromatic image performance of diffractive bifocal intraocular lenses: longitudinal chromatic aberration and energy efficiency,” Invest. Ophthalmol. Vis. Sci. 57(4), 2021–2028 (2016). [CrossRef]  

6. Y. Hu, Q. F. Cui, L. D. Zhao, and M. X. Piao, “PSF model for diffractive optical elements with improved imaging performance in dual-waveband infrared systems,” Opt. Express 26(21), 26845–26857 (2018). [CrossRef]  

7. F. Reda, M. Salvatore, F. Borbone, P. Maddalena, A. Ambrosio, and S. L. Oscurato, “Varifocal diffractive lenses for multi-depth microscope imaging,” Opt. Express 30(8), 12695–12711 (2022). [CrossRef]  

8. F. Castignoles, M. Flury, and T. Lepine, “Comparison of the efficiency, MTF and chromatic properties of four diffractive bifocal intraocular lens designs,” Opt. Express 18(5), 5245–5256 (2010). [CrossRef]  

9. W. D. Furlan, A. Martínez-Espert, D. Montagud-Martínez, V. Ferrando, S. García-Delpech, and J. A. Monsoriu, “Optical performance of a new design of a trifocal intraocular lens based on the Devil’s diffractive lens,” Biomed. Opt. Express 14(5), 2365–2374 (2023). [CrossRef]  

10. N. Lopez-Gil and R. Montes-Mico, “New intraocular lens for achromatizing the human eye,” Journal of Cataract and Refractive Surgery 33(7), 1296–1302 (2007). [CrossRef]  

11. S. Ravikumar, A. Bradley, and L. N. Thibos, “Chromatic aberration and polychromatic image quality with diffractive multifocal intraocular lenses,” Journal of Cataract and Refractive Surgery 40(7), 1192–1204 (2014). [CrossRef]  

12. E. Philippaki, L. Gobin, J. Mandoda, S. Lamy, and F. Castignoles, “Optical evaluation of new-design multifocal IOLs with extended depth of focus,” J. Opt. Soc. Am. A 36(5), 759–767 (2019). [CrossRef]  

13. D. Gatinel, C. Pagnoulle, Y. Houbrechts, and L. Gobin, “Design and qualification of a diffractive trifocal optical profile for intraocular lenses,” Journal of Cataract and Refractive Surgery 37(11), 2060–2067 (2011). [CrossRef]  

14. L. Zeng and F. Z. Fang, “Advances and challenges of intraocular lens design [Invited],” Appl. Opt. 57(25), 7363–7376 (2018). [CrossRef]  

15. J. L. Alio, A. B. Plaza-Puche, R. Fernandez-Buenaga, J. Pikkel, and M. Maldonado, “Multifocal intraocular lenses: An overview,” Surv. Ophthalmol. 62(5), 611–634 (2017). [CrossRef]  

16. J. A. Davison and M. J. Simpson, “History and development of the apodized diffractive intraocular lens,” Journal of Cataract and Refractive Surgery 32(5), 849–858 (2006). [CrossRef]  

17. M. S. Millan and F. Vega, “Extended depth of focus intraocular lens: Chromatic performance,” Biomed. Opt. Express 8(9), 4294–4309 (2017). [CrossRef]  

18. J. M. Artigas, J. L. Menezo, C. P. Martinez, and M. Llopis, “Image quality with multifocal intraocular lenses and the effect of pupil size: Comparison of refractive and hybrid refractive-diffractive designs,” Journal of Cataract and Refractive Surgery 33(12), 2111–2117 (2007). [CrossRef]  

19. R. Montés-Micó, E. Espaa, I. Bueno, W. N. Charman, and J. L. Menezo, “Visual performance with multifocal intraocular lenses: mesopic contrast sensitivity under distance and near conditions,” Ophthalmology 111(1), 85–96 (2004). [CrossRef]  

20. N. E. de Vries and R. Nuijts, “Multifocal intraocular lenses in cataract surgery: Literature review of benefits and side effects,” Journal of Cataract and Refractive Surgery 39(2), 268–278 (2013). [CrossRef]  

21. C. X. Xue and Q. F. Cui, “Design of multilayer diffractive optical elements with polychromatic integral diffraction efficiency,” Opt. Lett. 35(7), 986–988 (2010). [CrossRef]  

22. Y. F. Peng, Q. Fu, F. Heide, and W. Heidrich, “The diffractive achromat full spectrum computational imaging with diffractive optics,” ACM Trans. Graph. 35(4), 1–11 (2016). [CrossRef]  

23. G. Kim, J. A. Dominguez-Caballero, and R. Menon, “Design and analysis of multi-wavelength diffractive optics,” Opt. Express 20(3), 2814–2823 (2012). [CrossRef]  

24. L. L. Doskolovich, E. A. Bezus, A. A. Morozov, V. Osipov, J. S. Wolffsohn, and B. Chichkov, “Multifocal diffractive lens generating several fixed foci at different design wavelengths,” Opt. Express 26(4), 4698–4709 (2018). [CrossRef]  

25. X. Gao, C. X. Xue, C. Yang, and Y. Liu, “Optimization method of manufacturing for diamond turning soft-brittle materials’ harmonic diffractive optical elements,” Appl. Opt. 60(1), 162–171 (2021). [CrossRef]  

26. R. Mukherjee, K. Chaudhury, S. Das, S. Sengupta, and P. Biswas, “Posterior capsular opacification and intraocular lens surface micro-roughness characteristics: an atomic force microscopy study,” Micron 43(9), 937–947 (2012). [CrossRef]  

27. M. Lombardo, M. P. De Santo, G. Lombardo, R. Barberi, and S. Serrao, “Analysis of intraocular lens surface properties with atomic force microscopy,” Journal of Cataract and Refractive Surgery 32(8), 1378–1384 (2006). [CrossRef]  

28. A. Z. Zhang, “Multifocal diffractive lens design in ophthalmology,” Appl. Opt. 59(31), 9807–9823 (2020). [CrossRef]  

29. H. L. Liou and N. A. Brennan, “Anatomically accurate, finite model eye for optical modeling,” J. Opt. Soc. Am. A 14(8), 1684–1695 (1997). [CrossRef]  

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. Schematic of diffractive optical elements (DOEs): (a) Wavefront of DOEs; (b) actual profile of DOEs; (c) optical path diagram; (d) diffraction efficiency of monofocal DOEs with different harmonic numbers at different wavelengths; (e) diffraction efficiencies of bifocal DOEs with different phase delays at different wavelengths; (f) diffraction efficiency of an example multifocal DOE at different wavelengths.
Fig. 2.
Fig. 2. Continuous variable curve with multi-subperiods (CVCMS) design process.
Fig. 3.
Fig. 3. Results of CVCMS-designed MDOE between phase delays of two subperiods and diffraction efficiency. (a) Subperiod division proportion w1 = 0.3; (b) subperiod division proportion w1 = 0.5; (c) subperiod division proportion w1 = 0.7.
Fig. 4.
Fig. 4. Intraocular lens design (a) Liou–Brennan model of the eye (adapted from Ref. [29], Fig. 4). (b) Designed optical path diagram (red: Far point, green: Intermediate point, blue: Near point). (c) MTF results at design wavelength.
Fig. 5.
Fig. 5. Diffraction efficiency of MDOE design using the CVCMS method.
Fig. 6.
Fig. 6. Degree of separation comparison between design B and CVCMS method.
Fig. 7.
Fig. 7. DOE phase profile comparison between design B and CVCMS method. (a) Result by design B. (b) Result by CVCMS method.
Fig. 8.
Fig. 8. Schematic of the PSF at far focal point of an actual diffraction system (normalized by the maximum value), PSF values from left to right are compared at 500, 530, 550, and 600 nm. The PSF is calculated using a Fourier transform with a trifocal uniform design (a) or the proposed CVCMS design (b). The optical system uses the design in Section 4.2.

Tables (6)

Tables Icon

Table 1. Liou–Brennan model of the human eyea

Tables Icon

Table 2. Intraocular lens (IOL) design parameters

Tables Icon

Table 3. Modulation transfer function (MTF) design

Tables Icon

Table 4. MDOE parameter results with CVCMS design

Tables Icon

Table 5. Parameters of comparative designs A and B

Tables Icon

Table 6. Results with different design approaches

Equations (25)

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

$$f = \frac{{{r^2}}}{{2j{\lambda _0}}},$$
$$\varphi = 2\pi \mu \beta \gamma , $$
$$h = \frac{{{\lambda _0}\beta }}{{n^{\prime} - n}}, $$
$$l = \frac{{{\varphi _{{\lambda _{\max }}}} - {\varphi _{{\lambda _{\min }}}}}}{{{\varphi _{{\lambda _0}}}}}, $$
$$\mu = \frac{{{\varphi _\lambda }}}{{{\varphi _{{\lambda _0}}}}} = \frac{{{\lambda _0}({{{n^{\prime}}_\lambda } - {n_\lambda }} )}}{{\lambda ({{{n^{\prime}}_{{\lambda_0}}} - {n_{{\lambda_0}}}} )}}, $$
$${c_m} = \frac{1}{{2{\lambda _0}f}}\int_0^{2{\lambda _0}f} {\exp ( - \frac{{{\lambda _0}({{{n^{\prime}}_\lambda } - {n_\lambda }} )}}{{\lambda ({{{n^{\prime}}_{{\lambda_0}}} - {n_{{\lambda_0}}}} )}}\frac{{i\pi \beta {r^2}}}{{{\lambda _0}f}})\exp (\frac{{im\pi {r^2}}}{{{\lambda _0}f}})} d{r^2}. $$
$${\eta _m} = {|{{c_m}} |^2} = {\left|{\frac{1}{{2{\lambda_0}f}}\int_0^{2{\lambda_0}f} {\exp ( - \frac{{{\lambda_0}({{{n^{\prime}}_\lambda } - {n_\lambda }} )}}{{\lambda ({{{n^{\prime}}_{{\lambda_0}}} - {n_{{\lambda_0}}}} )}}\frac{{i\pi \beta {r^2}}}{{{\lambda_0}f}})\exp (\frac{{im\pi {r^2}}}{{{\lambda_0}f}})} d{r^2}} \right|^2}. $$
$${\eta _m} = {\left|{\frac{1}{{2{\lambda_0}f}}\int_0^{2{\lambda_0}f} {\exp ( - \frac{{i\pi \mu \beta {r^2}}}{{{\lambda_0}f}})\exp (\frac{{im\pi {r^2}}}{{{\lambda_0}f}})} d{r^2}} \right|^2}. $$
$${c_m} = \frac{1}{{2{\lambda _0}f}}\int_0^{2{\lambda _0}f} {\exp (\frac{{ - i\pi \mu \beta {r^2}}}{{{\lambda _0}f}})\exp (\frac{{m\pi {r^2}}}{{{\lambda _0}f}})} d{r^2}. $$
$${c_m} = \frac{1}{{2\pi }}\int_0^{2\pi } {\exp ( - i\mu \beta \cos (\theta ))\exp (im\theta )} d\theta. $$
$${c_m} = \frac{1}{{2\pi }}\left[ \begin{array}{l} \int_0^{2{w_1}\pi } {\exp ( - i\mu {\beta_1}\cos (\theta ))\exp (im\theta )} d\theta + \int_{2{w_1}\pi }^{2{w_2}\pi } {\exp ( - i\mu {\beta_2}\cos (\theta ))\exp (im\theta )} d\theta \\ + \int_{2{w_2}\pi }^{2{w_3}\pi } {\exp ( - i\mu {\beta_3}\cos (\theta ))\exp (im\theta )} d\theta + \cdots \cdots \\ + \int_{2{w_{k - 1}}\pi }^{2{w_k}\pi } {\exp ( - i\mu {\beta_k}\cos (\theta ))\exp (im\theta )} d\theta \end{array} \right], $$
$$\scalebox{0.96}{${c_m} = \frac{1}{{2\pi }}\left[ \begin{array}{@{}l@{}} \int_0^{2{w_1}\pi } {\exp ( - i\mu {\beta_1}\cos (\theta + {t_1}))\exp (im\theta )} d\theta + \int_{2{w_1}\pi }^{2{w_2}\pi } {\exp ( - i\mu {\beta_2}\cos (\theta + {t_2}))\exp (im\theta )} d\theta \\ + \int_{2{w_2}\pi }^{2{w_3}\pi } {\exp ( - i\mu {\beta_3}\cos (\theta + {t_3}))\exp (im\theta )} d\theta + \cdots \cdots \\ + \int_{2{w_{k - 1}}\pi }^{2{w_k}\pi } {\exp ( - i\mu {\beta_k}\cos (\theta + {t_k}))\exp (im\theta )} d\theta \end{array} \right].$}$$
$$\left\{ {\begin{array}{{c}} {{\beta_1}\cos (2\pi {w_1}) = {\beta_2}\cos (2\pi {w_1} + {t_2})}\\ {{\beta_2}\cos (2\pi {w_2} + {t_2}) = {\beta_3}\cos (2\pi {w_2} + {t_3})}\\ {{\beta_{k - 1}}\cos (2\pi {w_{k - 1}} + {t_{k - 1}}) = {\beta_k}\cos (2\pi {w_{k - 1}} + {t_k})} \end{array}} \right.. $$
$$\left\{ {\begin{array}{{c}} {{t_2} = \arccos \left( {\frac{{{\beta_1}}}{{{\beta_2}}}\cos ({2\pi {w_1}} )} \right) - 2\pi {w_1}}\\ {{t_3} = \arccos \left( {\frac{{{\beta_2}}}{{{\beta_3}}}\cos ({2\pi {w_2} + {t_2}} )} \right) - 2\pi {w_2}}\\ {{t_k} = \arccos \left( {\frac{{{\beta_{k - 1}}}}{{{\beta_k}}}\cos ({2\pi {w_{k - 1}} + {t_{k - 1}}} )} \right) - 2\pi {w_{k - 1}}} \end{array}} \right.. $$
$$- 1 \le \left( {\frac{{{\beta_{k - 1}}}}{{{\beta_k}}}\cos ({2\pi {w_{k - 1}} + {t_{k - 1}}} )} \right) \le 1. $$
$$\begin{array}{l} {c_m} = \frac{1}{{2\pi }}\\ \left[ \begin{array}{l} \int_0^{2{w_1}\pi } {\exp ( - i\mu {\beta_1}\cos (\theta ))\exp (im\theta )} d\theta \\ + \int_{2{w_1}\pi }^{2{w_2}\pi } {\exp ( - i\mu {\beta_2}\cos (\theta + \arccos \left( {\frac{{{\beta_1}}}{{{\beta_2}}}\cos ({2\pi {w_1}} )} \right) - 2\pi {w_1}))\exp (im\theta )} d\theta \\ + \int_{2{w_2}\pi }^{2{w_3}\pi } {\exp ( - i\mu {\beta_3}\cos (\theta + \arccos \left( {\frac{{{\beta_2}}}{{{\beta_3}}}\cos ({2\pi {w_2} + {t_2}} )} \right) - 2\pi {w_2}))\exp (im\theta )} d\theta + \cdots \cdots \\ + \int_{2{w_{k - 1}}\pi }^{2{w_k}\pi } {\exp ( - i\mu {\beta_k}\cos (\theta + \arccos \left( {\frac{{{\beta_{k - 1}}}}{{{\beta_k}}}\cos ({2\pi {w_{k - 1}} + {t_{k - 1}}} )} \right) - 2\pi {w_{k - 1}}))\exp (im\theta )} d\theta \end{array} \right] \end{array}, $$
$${\eta _m} = {|{{c_m}} |^2}. $$
$$\sigma = \left|{\frac{{{\eta_\lambda } - {\eta_{{\lambda_0}}}}}{{{\eta_{{\lambda_0}}}}}} \right|, $$
$$\begin{aligned} {{\bar{p}}_m} &= \left|{\frac{1}{{{\lambda_{\max }} - {\lambda_{\min }}}}\int_{{\lambda_{\min }}}^{{\lambda_{\max }}} {\frac{{{\eta_\lambda }}}{{{\eta_{{\lambda_{\max }}}} - {\eta_{{\lambda_{\min }}}}}}d\lambda } } \right|\\ &= \left|{\frac{1}{{{\lambda_{\max }} - {\lambda_{\min }}}}\int_{{\lambda_{\min }}}^{{\lambda_{\max }}} {\frac{{{{|{{c_m}_{_\lambda }} |}^2}}}{{{{|{{c_m}_{_{{\lambda_{\max }}}}} |}^2} - {{|{{c_m}_{_{{\lambda_{\min }}}}} |}^2}}}d\lambda } } \right|\end{aligned}, $$
$${\bar{p}_m} \to \max, $$
$$\sum\limits_N {\eta {{(\lambda )}_{{m_N}}}} \to \max \textrm{ }, $$
$$PSF(m,\lambda ) = \eta (m,\lambda ) \cdot PS{F_{doe}}(m,\lambda ), $$
$$MTF = |{FFT(PSF(m,\lambda ))} |. $$
$$\begin{aligned} {c_m} &= \frac{1}{{2{\lambda _0}f}}\sum\limits_k {\int_{2{\lambda _0}{w_{k - 1}}f}^{2{\lambda _0}{w_k}f} {\exp ( - i\frac{{\pi {\beta _k}{r^2}}}{{{\lambda _0}f}})\exp (i\frac{{m\pi {r^2}}}{{{\lambda _0}f}})d{r^2}} } \\ &\quad {w_0} \le {w_{k - 1}} \le {w_k} \cdots \le {w_T},{w_0} = 0,{w_T} = 1 \end{aligned}. $$
$$\begin{aligned} &{\varphi _k} = {\beta _k}\cos (\theta + {t_k}),\\ &{w_0} \le {w_{k - 1}} \le {w_k} \cdots \le {w_T},2\pi {w_{k - 1}} \le \theta \le 2\pi {w_k},{w_0} = 0,{w_T} = 1 \end{aligned}$$
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.