衍射光强分布的数值计算与实验对比
在波动光学的研究与工程应用中，衍射光强分布的解析与表征是核心课题之一。然而，除少数规则孔径（如单缝、圆孔）外，复杂光阑的衍射问题往往难以求得严格的解析解。此时，数值计算成为了不可或缺的手段。将数值计算结果与实际实验测量进行对比，不仅是验证数学模型有效性的关键步骤，也是优化光学系统设计的基础。
波动光学的衍射理论建立在惠更斯-菲涅耳原理之上，其核心思想是将波前上的每一点视为次级子波源，空间中某点的光场是所有子波干涉叠加的结果。
根据观测距离与衍射孔径尺寸的关系，衍射模型通常分为近场（菲涅耳衍射）与远场（夫琅禾费衍射）。在数学表达上，衍射积分可以统一表示为孔径函数与传递函数的卷积形式：
- 菲涅耳衍射积分：保留了二次相位因子，适用于任意距离，但积分求解复杂。
- 夫琅禾费衍射积分：远场近似下，积分简化为孔径函数的傅里叶变换，物理意义明确，计算相对简便。
对于任意形状的孔径，直接计算上述积分极为困难，必须借助离散化与数值算法进行逼近。
数值计算方法与实现策略
将连续的衍射积分转化为离散的数值计算，核心在于离散傅里叶变换（DFT）及其快速算法（FFT）的应用。
- 采样与离散化：将连续的孔径函数和观测面进行网格化采样。采样间隔必须满足奈奎斯特采样定理，以避免高频相位信息的混叠。特别是在菲涅耳衍射中，二次相位因子变化剧烈，对采样密度提出了极高要求。
- 算法选择：
- 直接积分法：通过求和直接逼近积分。适用于孔径极其不规则且采样点较少的情况，但计算复杂度高达 $O(N^2)$。
- FFT 卷积法：利用卷积定理，将空间域的卷积转化为频率域的乘积，再通过逆变换得到光场。计算复杂度降为 $O(N \log N)$，是当前最主流的方法。
- 角谱传播法：将光场分解为不同方向的平面波，在频域乘以传递函数后逆变换回空间域。该方法严格满足标量波动方程，无近场远场之分，且在处理短距离衍射时精度极高。
实验测量与数据提取
数值计算给出了理论预期，而实验测量则提供了物理真实。构建高精度的光强分布采集系统是进行有效对比的前提。
- 光源要求：通常采用单色性良好的激光源，以满足相干条件。同时需配合空间滤波器与准直透镜，获取均匀的平面波前。
- 环境控制：必须隔离环境振动与空气扰动，防止相干光程差的随机漂移。
- 光强采集：传统方法采用照相法结合显微密度计，效率低下。现代实验多采用科学级 CCD 或 CMOS 相机直接采集二维光强分布。需注意相机的线性响应范围，避免像素饱和引入非线性误差，并进行暗电流与背景噪声扣除。
数值计算与实验结果的对比分析
将数值仿真数据与实验图像进行对标，是一个从定性到定量的迭代过程。
- 定性比对（形态与对称性）：首先观察光强分布的整体轮廓。例如，计算结果中主极大与次极大的位置、衍射条纹的走向与对称性，是否与实验图像吻合。若出现整体偏移或不对称，通常提示实验光路存在倾斜或对准误差。
- 定量比对（极值与半高宽）：提取中心剖面的光强曲线，对比主瓣宽度（半高宽 FWHM）与旁瓣电平。数值计算往往基于理想标量波假设，而实验中由于偏振效应、有限尺寸光源的非相干叠加以及探测器的点扩散函数（PSF）卷积，会导致实验测得的峰值变低、极小值不为零（即条纹对比度下降）。
- 误差溯源与模型修正：当定量偏差超出容许范围时，需反向修正数学模型。常见的修正包括：在孔径函数中引入边缘透射率的渐变（模拟实际加工的倒角或蚀刻残留）；在计算中叠加系统的点扩散函数；或引入部分相干度模型以更真实地反映实际光源特性。
应用全景与展望
衍射光强分布的数值计算与实验对比方法，在现代光学工程的各个领域发挥着基石作用：
- 超精密光刻：在极紫外（EUV）光刻机中，掩模图形的衍射效应极为显著。通过高精度数值计算与实验标定，可优化离轴照明与光学邻近效应校正（OPC）策略。
- 微纳光学器件设计：超表面、衍射光学元件（DOE）的相位调控设计，高度依赖严格的耦合波分析（RCWA）或时域有限差分法（FDTD）计算，其实际加工的微观形貌偏差需通过实验对比来评估与补偿。
- 天文光学：大型望远镜的支撑结构会产生衍射尖峰，影响系外行星的观测。通过数值仿真与实测对比，可指导支撑结构的优化以抑制衍射杂散光。
综上所述，数值计算与实验测量并非孤立存在，而是波动光学研究中“理论预测-实验验证-模型优化”闭环的双翼。随着计算算力的提升与高动态范围探测技术的发展，两者的深度融合将持续推动现代光学系统向更高效、更精密的方向演进。