首页/文章/ 详情

基于网格相位模型的环形光设计与仿真

2月前浏览794
本文将介绍一种基于衍射光学元件的环形光源设计方法,旨在解决光学干涉检测中由高空间相干性引起的“牛顿环”噪声问题。研究通过理论建模、相位计算与仿真验证,实现了半径可控、强度均匀的环形光输出,并开发了配套的MATLAB GUI工具,大幅提升设计效率。

一、研究背景:为什么需要环形光源?

在光学干涉测量中,高相干性的激光光源(如氦氖激光器)容易在光学元件表面缺陷(如灰尘、划痕)处产生散射,形成“牛顿环”等相干噪声,严重影响测量精度与图像质量。
传统方法通过光阑滤除杂散光,效果有限。环形光源通过扩大光源有效面积、降低空间相干性,可显著抑制此类噪声,已在高精度干涉仪(如ZYGO Verifire系列)中得到成功应用。

二、设计思路与方法

1. 光学系统建模

基于几何光线追迹,建立平面波与球面波入射的光学模型,推导衍射面相位分布函数。核心思想是:通过衍射面对波前进行相位调制,使出射光在像面上汇聚为指定半径的环形分布。

2. 相位计算与拟合挑战

衍射面相位分布计算,计算公式推导过程(1.1)-(1.10)。
a.平面波入射
设各点坐标分别为
                          (1.1)
                                      (1.2)
                         (1.3)
b.球面波入射
设各点坐标为        
在 中     
                                           (1.4)
 中           
                                                            (1.5)
由折射定律可得
                                     (1.6)
                        (1.7)
                      (1.8)
                                                          (1.9)
                           (1.10) 
衍射面的相位函数在中心区域存在奇点,传统多项式面型(如ZEMAX中的Binary1-4、Zernike面型)难以准确拟合。
解决方案:采用ZEMAX中的网格相位面型(Grid Phase),通过200×200的离散栅格直接映射相位数据,实现高精度拟合。

3. 微纳结构参数设计

二元光学器件的掩膜设计主要有三个过程,首先是求出相位分布函数,接下来对相位函数进行相位压缩和量化,最后将得到的数据进行转化输出。
根据标量衍射理论分析可知,无吸收且易于加工的纯相位型二元光学器件具有较高的衍射效率。在本章已经计算出了相位分布函数,接下来针对二台阶数衍射元件,进行相位压缩和量化。
二元光学器件引入的光程差分布,会如图3所示,被压缩到内,对于经典的二元光学器件,m=1。
为了能够实际进行加工,需要计算表1中的参数。
表1 微纳结构参数

参数名称

含义

刻蚀深度

相位延迟所对应的实际材料高度变化。

衍射级次

用于计算效率的阶次(通常使用第一级衍射)。

衍射效率

各阶衍射光功率分布,决定设计利用率。

占宽比

在二级结构中,高低台阶的比例

最小线宽

加工时能实现的最小 feature size,限制相位分布精度。

3.3.1刻蚀深度的计算

对于相位型光学元件,某一波长下相位延迟φ和材料刻蚀深度d的关系为:
                                                              (1.11)
当入射波长为555 nm,材料选择BK7,折射率=1.5168时,对于二台阶结构,只需设置φ=π,即半波相位延迟即可达到相位反转,有以下计算结果。
                                             (1.12)
故其刻蚀深度为一个波长左右时,对应的衍射效率最高

3.3.2二台阶结构提取(相位二值化)

若要将连续的相位函数转化为二值化的结构,可以编写如下的MATLAB代码,phase_binary = phase_matrix > pi, 超过π为1,低于π为0。从而可以得到一个只有0/1的位图(掩膜),代表低/高台阶区域

3.3.3占宽比计算

通过编写MATLAB程序,从二值化后的结构中直接计算,语句如下。
duty_ratio = sum(phase_binary(:)) / numel(phase_binary),高台阶所占比例占宽比将会控制mask版图的“黑白”面积比例,影响曝光剂量与结构形成[28]。通过对二值化的阈值进行扫描,找到最大的占宽比为50.03%。

3.3.4衍射效率计算

对于二台阶相位型光学元件,傅里叶级数展开中,二台阶相位的衍射效率主峰在m=±1,尝试用1级。可通过傅里叶变换来估计其衍射效率。也可在ZEMAX中分析工具直接查看衍射阶次能量分布。当占宽比为50.03%时,对应的衍射效率为40.53%,如图4所示。

3.3.5最小线宽计算

最小线宽由相位函数的局部梯度决定,通过相位梯度确定相邻环间距Δr,其最小值即为理论最小线宽。具体步骤分四步:1.确定相位函数2.计算相位梯度3.找到最大梯度位置4.计算相邻环间距5.确定最小线宽,见表2。
表2 最小线宽计算

步骤

公式

相位梯度

最大梯度位置

r = R(元件边缘)

最小线宽


三、仿真平台与流程

1. MATLAB GUI 程序

开发了集成化设计工具,支持参数输入、相位计算、数据导出与加工图纸生成,极大提升设计效率。
 

2. ZEMAX 仿真验证

将生成的.dat相位文件导入ZEMAX,在Grid Phase面型中加载,完成光路仿真与结果分析。

四、仿真结果与参数

1. 平面波入射 · 环形光半径0.2 mm

在MATLAB中生成.DAT文件,并将其导入栅格相位面型的附加相位数据中,设置好光学系统各参数,如图6所示。
经过ZEMAX模拟后,在像面生成了半径为0.2mm的环形光,光线通过率为100%,如图7所示。
用于加工的衍射面掩膜图纸以及环带数据如图8所示。
如表3所示,当入射波长为555nm,衍射级次为1级时,经过阈值扫描,以1.2π为界,占宽比为49.99%,衍射效率为40.53%,刻蚀深度为537nm,最大相位梯度90.57 rad/mm,最小线宽0.069mm。
表3 设计报告

项目

参量

项目

参量

波长

555 nm

元件尺寸(直径)

Φ=2.00 mm

最佳阈值

3.2798 rad1.04π

刻蚀深度

537.0 nm

占宽比

49.99%

材料

BK7

理论 ±1  阶效率

40.53%

最小线宽

0.069 mm

3. 光强均匀性分析

通过ZEMAX几何图像分析功能,可观察环形光在X/Y截面上的强度分布,整体均匀性良好,符合干涉检测需求。
当观察X/Y截面的强度分布图时,发现,呈下图11,在环形上强度分布为“针尖状”,中间区域最大,由于输出的环形光,光线宽度约为0.02 mm,所以对应着这种强度分布,整体环上强度分布均匀,符合设计要求。
免责声明:

本文旨在对基于仿真软件的技术研究、分析方法及结果进行非商业性的原理介绍、案例分享与行业科普。文章内容涉及仿真软件的功能描述、应用示例及一般性操作流程,综合参考了公开的学术文献、官方技术资料及行业公认的实践认知,仅供学习、交流与参考之用。

来源:旋算仿真工作室
静力学瞬态动力学振动疲劳碰撞电力电子MATLAB芯片UMSimulink焊接理论电机材料科普数字孪生Zemax
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2026-06-16
最近编辑:2月前
旋算仿真工作室
专业一对一仿真技术服务
获赞 7粉丝 14文章 40课程 0
点赞
收藏
作者推荐

夏克-哈特曼波前传感器:现代光学的精密“探针”

在先进光学系统的研发和应用中,有一类仪器如同一位敏锐的诊断医生,能够精准捕捉并量化光线传播过程中产生的微小“畸变”或“缺陷”,它就是夏克-哈特曼波前传感器。不同于直接成像的相机,它专注于解析光波的“波前”这一物理本质,其测量结果直接决定了自适应光学系统的校正精度、高端光学产品的质量,以及前沿科学实验的可靠性。图1 夏克-哈特曼波前传感器 一、从原理到部件:如何测量光的“平坦度”? 理解夏克-哈特曼传感器的关键在于理解它如何将复杂的整体波前分析转化为大量简单的局部测量。 它的核心工作原理基于几何光学中的光斑偏移追踪法。当一束携带波前畸变信息的光入射到传感器时,首先会经过一个关键部件——微透镜阵列。这个阵列由数万甚至数十万个口径微小(通常直径在几十到几百微米)的透镜规整排列而成。阵列中的每一个微透镜,都会将入射光束中对应孔径范围内的那部分“子波前”独立地会聚到后端的高灵敏度图像传感器(如CCD或CMOS)上,形成一个对应的焦点光斑。 图2 夏克哈特曼测量原理 如果入射光是一个完美的平面波(即无畸变),那么所有子孔径形成的光斑将整齐地排列在一个标准的理想网格位置上。然而,当波前存在局部倾斜或起伏时,每个子波前的传播方向便会发生微小的偏折。这种偏折直接导致其聚焦光斑在焦平面上的位置,相对于理想参考位置产生横向位移。传感器通过精密算法(通常采用质心算法或匹配滤波算法)实时检测出每个光斑在X和Y方向的位移量(Δx, Δy)。这一位移量,在数学上正比于该子孔径区域内波前相位的局部一阶导数(即波前斜率)。 最终,系统通过整合所有子孔径测量到的斜率数据,利用波前重构算法(如经典的Southwell区域法或基于泽尼克多项式的模式法),即可重建出完整的二维波前相位分布图,直观展示出光波经过被测对象后的畸变全貌。整个过程可以在毫秒甚至微秒量级内完成,实现了波前畸变的实时、定量、全场测量。 二、核心性能指标与技术挑战 评价一台夏克-哈特曼传感器的性能,主要看以下几个关键指标: 1. 空间分辨率:取决于微透镜阵列的子孔径数目和排列方式。子孔径越多,对波前的采样越精细,越能反映高频畸变细节。 2. 动态范围:指传感器能够测量的最大 波前斜率或相位变化范围。这与微透镜的焦距(f数)和探测器的像素规模有关。高动态范围与高空间分辨率往往存在设计矛盾。 3. 测量精度与灵敏度:即能够分辨的最小波前变化,可达到纳米量级。这受限于探测器的噪声、光学装调精度以及算法的稳定性。 4. 时间带宽:由探测器的帧率和处理算法的速度共同决定,决定了系统跟踪动态波前变化(如大气湍流)的能力。 当前的技术挑战主要在于扩展动态范围的同时保持高分辨率,以及处理弱光信号(如天文导星)时的信噪比问题。业内解决方案包括采用可编程液晶空间光调制器动态调整子孔径有效焦距,或开发更先进的亚像素光斑定位算法。 三、实际应用 凭借其非干涉、光路简单、实时性强的优点,夏克-哈特曼传感器已成为多个高技术领域的标准诊断工具。 1. 自适应光学系统 这是其最成功的应用领域。在大型地基天文望远镜(如中国的LAMOST、美国的Keck望远镜)中,传感器实时探测由大气湍流造成的星象波前畸变,并将数据反馈给变形镜,驱动其表面进行每秒数百至上千次的形状校正,从而有效克服大气抖动,使望远镜的分辨率逼近理论衍射极限。在惯性约束核聚变装置(如美国国家点火装置NIF)中,它被用于校正高能激光束路中的像差,确保能量精准聚焦于靶丸。 图3 共轭式自适应光学系统结构示意图 2. 精密光学计量 在光学加工与系统集成中,该传感器是检测光学元件面形、光学系统像差的利器。相较于干涉仪,它对环境振动和气流扰动不敏感,更适合在线检测与装调。例如,在极紫外光刻机物镜的集成调试中,就需要使用专门设计的夏克-哈特曼传感器来精确标定和校准复杂多反射镜系统的波前误差。 图4 夏克-哈特曼传感法检测自由曲面光路图 3. 眼科医学诊断 在眼科,它以波前像差仪的形式出现。通过分析从人眼视网膜反射出来的光波前,可以绘制出包含近视、散光等低阶像差以及球差、彗差等高阶像差的个性化眼睛“光学地图”。这张地图是实施个性化准分子激光手术(波前引导LASIK) 和定制高端人工晶体(如多焦点、EDOF晶体) 的根本依据,旨在追求超越常规矫正的“超视力”视觉质量。 图5 屈光测量系统光路原理图 4. 激光光束诊断与质量控制 对于高功率激光器、光纤激光器等,传感器用于直接测量输出激光的波前质量,评估其光束质量因子M²、聚焦特性等,是优化激光器设计、确保加工或医疗应用效果的关键环节。 图6 哈特曼波前传感器测试系统 四、前沿发展与展望 夏克-哈特曼技术本身也在不断演进。一方面,硬件上正与超表面光学结合,利用纳米结构实现更轻薄、功能集成度更高的新型波前传感方案。另一方面,与人工智能的结合正在兴起,利用深度学习算法直接从光斑阵列图像中快速、鲁棒地反演波前或预测光学系统状态,为复杂动态环境下的波前感知开辟了新途径。 图7 基于超表面的Shack-Hartmann波前传感器 此外,其原理也被拓展至更前沿的科学探索中。例如,中国科学技术大学的研究团队曾创新性地将其应用于量子光学领域,实现了对光子纠缠态空间波函数的测量,展现了这一经典技术在量子信息科学中的新生命力。 五、结语 从实验室的精密测量到探索宇宙的巨型望远镜,从改善人类视力的手术台到制造尖端芯片的光刻机,夏克-哈特曼波前传感器作为连接光学现象与定量数据的桥梁,始终在幕后发挥着至关重要的作用。它不仅是工程师解决实际光学问题的“眼睛”,更是科学家探索光与物质相互作用的精密“探针”。随着光学技术不断向更精密、更复杂、更动态的方向发展,对波前传感技术的要求也将水涨船高,持续推动着这一领域向着更高速度、更高精度、更强适应性的未来迈进。 免责声明:本文旨在对夏克-哈特曼波前传感器的原理、应用及发展进行非商业性的技术介绍与行业科普。文章内容综合了公开的技术资料、学术文献及行业实践信息,仅供交流与参考之用。文中所述的具体技术参数、性能指标及实施效果,可能因不同厂商的产品设计、应用场景和测试条件而异,读者在实际科研、工程或商业决策中,应以相关设备官方技术文档和实测数据为准。本文尽力确保内容的准确性,但作者及发布方不对其完整性、时效性或是否适用于特定用途作出任何明示或默示的保证。任何基于本文内容所采取的行动及由此产生的后果,均需由读者自行承担。文中如涉及任何第三方商标、技术或数据,其所有权均归各自权利人所有,本文引用仅为说明之目的,不构成任何授权或关联。来源:旋算仿真工作室

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈