白內障手術是當今最常見的外科手術之一，在該手術中，患者的晶狀體由于光散射增加而變得渾濁，從而被人工晶狀體（IOL）取代。隨著白內障人群越來趨于越年輕化，對優(yōu)質鏡片的需求不斷增長，以提高可實現的圖像質量并解決無需眼鏡聚焦的問題。衍射IOL通過同時創(chuàng)建多個焦點來提供近距離和遠距離的清晰視覺，從而提供了一種可行的解決方案，在本文中我們演示了如何通過使用用戶自定義表面（UDS）DLL來擴展Zemax OpticStudio的功能，以提供衍射式人工晶狀體透鏡的真實模型。最后，我們還將討論應用了區(qū)域分解模型相對于使用內置衍射表面類型的序列分解方法的優(yōu)勢。

簡介

隨著白內障患者的需求不斷增加，人工晶狀體制造商投入更多的時間和資源到高級鏡片的研究和精密設計中。為了在很寬的物體距離范圍內提供良好的圖像質量，最終的目標是再現原始晶狀體的調節(jié)能力。用人造元件直接模仿人眼的自然過程面臨著幾個挑戰(zhàn)，因此這仍然是一個尚未解決的問題。然而，衍射人工晶狀體可以同時為多個觀察距離提供一個易于使用的解決方案。本文展示了如何使用光線追蹤和衍射分析來實現基于真實表面形狀的浮雕型衍射透鏡的真實模型，并展示了該模型在全面評估系統(tǒng)性方面的優(yōu)勢。

基本設計概念

1.階次分解

OpticStudio序列模式下的內置衍射表面模型依賴于階次分解，在此方法中，需要選擇單個衍射順序，然后衍射光焦度( Diffractive Power )由額外的相位貢獻代表，與折射率和表面矢高無關。使用這種方法，階次傳播可以通過從物體到圖像的光線或通過出射瞳孔的標量衍射來建模。這種方法提供了分析單個階次的簡單解決方案，對于使用單個目標衍射階的應用特別有益。使用此方法設計衍射人工晶狀體的工作原理和應用示例可以聯系工作人員了解。

然而，上述分解模型中存在一些不足。首先，由于相位函數只是對經過折射或反射表面的光線施加額外的相位改變，因此該模型不考慮通過衍射元件的真實光線路徑，因此忽略了波長色散以及某些其它像差。此外，這種表面模型沒有考慮衍射效率。最后，必須創(chuàng)建一個多重結構系統(tǒng)來逐個模擬不同的衍射級次。

2.區(qū)域分解

相反，使用區(qū)域分解，可以一次準確地考慮多個衍射階次，并且該方法通過模擬衍射元件的實際形狀來固有地考慮波長色散和衍射效率。這使得創(chuàng)建先進的IOL模型成為可能，其中不同的順序旨在為多個觀看距離提供清晰的視野，從而取代人眼。

區(qū)域分解模型假設分散元件的區(qū)域寬度遠大于波長，并且光學特性在區(qū)域內表現平滑。在這種情況下，幾何光學近似和光線追蹤可用于描述從衍射表面的一側到另一側的傳播。這也意味著這些區(qū)域在近場中可以被視為傳統(tǒng)的折射/反射光學元件，而在遠場光分布中只能通過標量衍射分析來計算。

在OpticStudio 中，PSF計算實現了這一精確過程，幾何光線追跡的結果疊加從出瞳到像面的衍射效應。由于通過衍射元件的相位變化是基于幾何光學計算的，因此當衍射表面位于出瞳或其共軛位置之一（入瞳或孔徑光闌）時，區(qū)域分解最適合的。

人工晶狀體設計和仿真是一個合適的案例，它符合上述標準，因為植入的人工晶狀體通常放置在瞳孔之后（瞳孔作為孔徑光闌）。根據通常的做法，當孔徑光闌和出瞳之間的菲涅爾衍射可忽略時，可以使用區(qū)域分解來有效地模擬衍射IOL。

使用 UDS DLL 的衍射表面模型

為了利用上述區(qū)域分解方法，我們實現了一個新的用戶定義表面DLL，其中可以通過分析描述浮雕型衍射表面的矢高輪廓。除了精確分析衍射光學元件（DOEs）的性能外，使用UDS DLL的參數化形狀表示還可以對這些衍射表面進行優(yōu)化和公差分析。有關如何使用自定義 DLL 擴展 OpticStudio 的功能以及如何編譯新解決方案的更多詳細信息，請聯系工作人員了解。

在使用序列表面 DLL 時，OpticStudio 有 10 種不同的方式與 DLL 交互和交換數據。這些方案表示常規(guī)信息、參數名稱和安全數據傳輸，以及布局圖、近軸和實際光線追蹤計算。不同的功能是在DLL的不同情況下定義的。

在這個模型中，我們應用了一個簡單的旋轉對稱衍射結構，具有統(tǒng)一的浮雕臺階高度，添加在代表基底面的標準表面之上。為了能夠與內置的OpticStudio解決方案進行模擬比較，我們用偶數非球面多項式描述了浮雕形狀。因此，表面矢高由以下公式給出：

在上式中，mod表示取模運算，c是曲率，即半徑的倒數，k是圓錐常數，r是徑向坐標，h是統(tǒng)一的浮雕臺階高度。

ai為偶次非球面系數、h為步高。首先在DLL的Case 1中定義傳播算法，參數列標題名稱。然后，Case 3描述基于上述公式的表面矢高，以便在布局圖中繪制。Case 4 考慮近軸光線追蹤結果，但由于區(qū)域分解方法需要在光線追跡之上進行衍射分析，該方法僅適用于實際光線追蹤，因此我們忽略了這一步。這意味著在近軸近似中，我們的模型表現為標準曲面。最后，Case 5，計算實際光線追跡結果。為此，我們實現了兩個解決方案，一種近似解析算法和一種迭代算法，這將在下面將討論。

光線傳播算法

在復雜表面形狀的情況下，無法通過分析，確定光線-表面相交坐標，因此對于標準表面以外的內置表面類型，OpticStudio 應用迭代算法來查找數值解。這也可以是用戶自定義DLL 的一種方法。然而，由于迭代方法的計算效率低于直接計算，除了通常應用的迭代解決方案外，我們還實現了基于局部線性化的近似閉合解解決方案[1，5]。

>

在后一種替代算法中，我們分別處理基板的矢高和額外的浮雕高度。首先，我們確定與基板準確的光線相交坐標（x0, y0,z0），這可以通過分析來完成，因為基板具有標準表面形狀。然后，作為下一步，我們根據局部浮雕高度（Δz=zDOE(x0,y0)）和給定位置的斜率（x0, y0, z0+Δz)來估計光線浮雕交點。與切平面的估計交點（x，y，z）可以通過求解線性方程再次解析計算。這種直接近似計算可以比默認迭代方法快 30%，而且不會在結果中引起任何重大錯誤。該過程如下圖所示。

在眼睛模型中模擬人工晶狀體性能

雙焦點人工晶狀體設計

為了證明新的衍射表面DLL的適用性和優(yōu)勢，我們基于文獻[1]實現了一個理想的衍射透鏡模型，該模型在設計波長處將光均勻地分布到零階和第一衍射階。（雖然這兩個階的功率相等，但施加的衍射面也會將一小部分功率發(fā)送到更高階。）根據人工晶狀體的ISO標準，我們將中心波長設置為e線，即λ0=546.07nm。人工晶狀體設計屈光度為 P0=22.5 D，衍射附加光焦度為3.5 D。我們使用折射率為n=1.4625的模型材料求解對透鏡材料Benz25進行建模，而周圍的介質天然鹽水則通過折射率為n0=1.3343的模型材料求解來描述。

為了獲得理想的透鏡，其中第一衍射階聚焦在EFL的距離上，相對于球面高斯參考波前，在區(qū)域j和j+1邊界處的光程差必須為jλ，其幾何上可以表示為下面的方程：

這意味著區(qū)域邊界位于以下位置：

由于 EFL>>λ，我們可以忽略后半部分并應用以下近似值：

因此，增加衍射面的光路以實現上述光路rj區(qū)域邊緣和產生100%衍射效率的一階如下：

為了實現零階和一階相等的衍射效率，還必須考慮恒定α=0.5倍增因子。因此，理想的雙焦點衍射人工晶狀體的表面矢高可以描述為：

基于此，在鏡頭數據編輯器中，我們將 2nd階多項式系數a1=αn/2EFL(n-n0)=6.82E-3，步高h=αλ0/（n-n0）=2.13E-3mm。衍射輪廓的矢高增加，即去除了基本半徑的表面矢高圖，如下圖所示：

最后，假設對稱雙凸透鏡具有1.0 mm的厚度和P0=22.5 D的基底光焦度，通過使用造鏡者公式計算透鏡基底的基底半徑，半徑為 11.353 mm。透鏡前表面的圓錐常數設置為k=0，后表面的圓錐常數經過優(yōu)化，以實現零衍射階的衍射極限性能，得到k=-5.8的值。

用于生產線測試的 ISO 標準眼模型

為了驗證人工晶狀體模型，我們將衍射UDS DLL納入ISO 11979-2標準眼模型，該模型專為生產線測試眼科植入物的光學特性而設計[2]。眼模型包含一個幾乎無像差的角膜，然后將人工晶狀體放置在兩個平面窗口之間的液體介質中：

使用快速對焦工具優(yōu)化系統(tǒng)的后焦距，以獲得最小的RMS波前誤差，從而在零階和一階焦點之間產生圖像位置。這也證明了在光線追蹤之上，還需要標量衍射分析來正確描述透鏡的遠場行為。此外，在這個中間位置，波前圖清楚地顯示了相鄰區(qū)域之間的半波差，這與理論預期一致。

當考慮衍射效應時，例如通過FFT PSF或MTF計算，則所有衍射級都在單個配置中被精確建模，并且衍射效率本質上是由模型的性質考慮的。FFT Through Focus MTF分析的結果在50 lp/mm頻率下如下所示。零級和一階焦平面的峰值衍射效率0.34接近理論值。

與內置衍射表面模型的比

為了將新的實際UDS DLL模型與OpticStudio中內置的衍射表面模型進行比較，我們創(chuàng)建了一個多重結構系統(tǒng)，其中我們使用Binary2衍射表面類型來描述針對不同階次的相位疊加。對應于前面討論的相同理想雙焦點透鏡的相位曲線可以用以下公式描述：

了解衍射元件的矢高和相位輪廓之間的轉換請聯系工作人員。

該模型包含 5 種配置，對應于使用Binary2相位表示的0th和 1st的衍射階次，并到 0th和 1st階次，基于浮雕的 UDS DLL 模型的中間幾何焦點。0th和 1st通過基于Binary2模型的最小RMS波前誤差的優(yōu)化來確定焦平面，并在新的UDS DLL模型中拾取相同的位置進行分析。

焦平面上的FFT PSF結果清楚地表明了階次分解模型和區(qū)域分解模型之間的差異。雖然內置的Binary 2 模型僅考慮所選階次，但 UDS dll考慮了多階次的共同作用結果。下面在一階焦點處的對數尺度圖像中可視化了此現象。左側的假彩色FFT PSF圖對應于階次分解模型，而右側的圖顯示了區(qū)域分解結果。相同PSF結果在中間同一行，即Y=0的位置，如下所示。

此外，FFT MTF結果表明，在區(qū)域分解模型中準確考慮了衍射效率，而在階次分解模型中則沒有。

多色結果

最后，為了證明在區(qū)域分解模型中追蹤真實光線路徑也會考慮透鏡中的波長色散，我們分析了可見光譜范圍內的系統(tǒng)性能。為了在擴展范圍時保持設計波長作為參考，我們通過選擇 F’, e, C’ (可見)預設來使用 F’, e, 和 C’線。FFT Through Focus MTF 圖描述了一階焦平面（左側駝峰）如何隨著波長的增加而遠離零階焦平面（右側駝峰）。同時，在較長的波長下，更多的能量被衍射到零階，而更少的能量被衍射到一階，正如理論所預期的那樣。

結果

在本文中，我們演示了如何使用用戶自定義表面DLL來擴展OpticStudio用于模擬浮雕型衍射人工晶狀體的功能。由于人工晶狀體被植入到瞳孔（孔徑光闌）附近，幾何光線追跡結果疊加從出瞳到像面的衍射效應，可以準確地重建多衍射階次的光場分布。我們通過在ISO標準眼方案中模擬理想的雙焦點人工晶狀體設計來測試和驗證我們的模型。最后，我們介紹了新的區(qū)域分解模型相對于內置的階次分解模型的優(yōu)勢。

引用

A . Nemes-Czopf, D. Bercsényi, G. Erdei. Simulation of relief type diffractive lenses in ZEMAX using parametric modelling and scalar diffraction. Applied Optics, 58(32):8931-8942 (2019).

Ophthalmic implants—Intraocular lenses—Part 2: Optical properties and test methods, ISO 11979-2:1999.

A. S. Gutman, I. V. Shchesyuk, V. P. Korolkov. Optical testing of bifocal diffractive-refractive intraocular lenses using Shack-Hartmann wavefront sensor. Proceedings of SPIE - The International Society for Optical Engineering, 7718 (2010).

T. Eppig, K. Scholz, A. Langenbucher. Assessing the optical performance of multifocal (diffractive) intraocular lenses. Ophthalmic and Physiological Optics, 28:467–474 (2008).

H. Sauer, P. Chavel, G. Erdei. Diffractive optical elements in hybrid lenses: modeling and design by zone decomposition. Applied Optics, 38:6482–6486 (1999).

D. A. Atchison, G. Smith. Optics of the Human Eye. Butterworth-Heinemann, UK (2000).