针对引力波探测星间激光干涉仪中望远镜光学参数的综合优化方法：以最小化倾角-长度耦合系数并满足所需波前像差为目标

《Results in Engineering》：Comprehensive Optimization Approach for Telescope Optical Parameters Targeting Minimal Tilt-to-Length Coupling Coefficient with Required Wavefront Aberration in Intersatellite Laser Interferometer for Gravitational Wave Detection

　　

望远镜是确保空间引力波探测中星间激光干涉仪（ILI）皮米级测量灵敏度的关键光学单元。它需要满足倾角-长度耦合系数（TTL coupling coefficient）小于0.6 pm/nrad且波前像差（wavefront aberration）小于λ/40的关键要求。本研究采用高斯光束分解法（Gaussian beam decomposition method）计算望远镜的组合几何和非几何倾角-长度耦合噪声（TTL coupling noise），该方法具有高计算效率和精度，优于现有的衍射积分法（diffraction integral methods）或光线追迹法（ray-tracing methods）。然后，研究人员提出了一种综合优化方法，以同时最小化组合TTL耦合噪声，同时将波前像差维持在λ/40以下，弥补了现有方法仅关注单一性能指标的不足。该方法创新性地将贝叶斯算法（Bayesian algorithm）与Zemax软件相结合，用于望远镜的高效多参数全局优化，前者旨在最小化组合TTL耦合噪声，后者进一步用于波前控制。结果表明，TTL耦合系数为0.134 pm/nrad，显著低于0.6 pm/nrad的要求，且波前像差低于λ/40，优于先前的设计。该方法有效平衡了TTL噪声抑制与波前质量，满足了LISA/Taiji/TianQin等任务中ILI的严苛要求。

**引力波探测星间激光干涉仪中望远镜光学参数的综合优化方法：论文解读**

**研究背景与问题**

引力波是爱因斯坦广义相对论的基本预言，由大质量天体螺旋或并合产生的时空涟漪。2015年，地面激光干涉引力波天文台（LIGO）首次直接探测到双黑洞合并产生的引力波，标志着天文学领域的里程碑。然而，地面探测器受低频地震噪声和地球重力场梯度噪声影响，灵敏度局限于1 Hz以上频段。空间引力波探测则能拓展至更低频段，目前主要国际任务包括欧洲空间局主导的激光干涉空间天线（LISA）、中国科学院主导的太极（Taiji）计划以及中山大学主导的天琴（TianQin）计划。这些任务通过多颗卫星组成星座，形成巨型的星间激光干涉仪（ILI），要求位移测量灵敏度达到10^-12～10^-11 m，距离为10⁸～10⁹ m。望远镜是ILI的关键单元，用于减小出射激光发散角并汇聚入射光束，必须满足倾角-长度耦合系数（TTL coupling coefficient）小于0.6 pm/nrad以及波前像差（wavefront aberration）小于λ/40的严格要求，以保证测量灵敏度。

现有研究对望远镜光学参数优化存在不足：计算TTL耦合噪声的方法分为几何光学法（计算几何TTL噪声）和波动光学法（计算非几何TTL噪声），但多数研究仅针对单一性能指标（如仅优化几何TTL噪声或仅分析非几何TTL噪声），缺乏同时考虑两者并约束波前像差的综合优化框架。此外，衍射积分法计算精度高但计算量巨大，不适用于大孔径望远镜系统；而高斯光束分解法能在保证精度的同时显著提升计算效率。因此，本研究旨在填补这一空白，提出一种综合优化方法，同时最小化组合几何与非几何TTL耦合噪声，并保持波前像差在要求范围内。

**研究概述与结论**

研究人员提出了一种结合贝叶斯算法（Bayesian algorithm）与Zemax软件的综合优化方法，用于望远镜光学参数的多参数全局优化。通过高斯光束分解法计算组合TTL耦合噪声，并以该噪声为优化目标，同时利用Zemax进行波前控制。优化后的望远镜TTL耦合系数降至0.134 pm/nrad，显著低于0.6 pm/nrad要求，波前像差均方根（RMS）最大为0.014λ，满足λ/40要求，性能优于先前设计。该方法有效平衡了TTL噪声抑制与波前质量，满足LISA/Taiji/TianQin等任务中ILI的严苛要求。论文发表在《Results in Engineering》。

**关键技术方法**

本研究主要采用以下关键技术方法：1）高斯光束分解法（Gaussian beam decomposition method）：将初始光场分解为一系列子高斯光束，分别传播通过光学系统后相干叠加，用于精确计算组合几何与非几何TTL耦合噪声，兼顾计算效率与精度。2）贝叶斯优化算法（Bayesian optimization algorithm）：构建目标函数的概率代理模型（高斯过程回归，Gaussian process regression），通过采集函数（采集函数，acquisition function）平衡探索与开发，高效求解多参数、黑箱、计算密集的优化问题。3）Zemax软件：用于优化望远镜的波前像差、出射光束角度和出瞳位置等基本光学性能，作为贝叶斯优化迭代中的辅助步骤。所有计算基于中山大学天琴中心（TianQin Research Center）的光学参数，无样本队列来源。

**研究结果**

**5.1 优化后望远镜的性能**

通过200次迭代优化（前100次采用上置信界（UCB）采集函数探索参数空间，后100次采用期望改进（EI）采集函数精细搜索），最优TTL耦合系数从0.882 pm/nrad降至0.122 pm/nrad，最终稳定在0.134 pm/nrad（考虑指向偏差后的最大TTL耦合系数）。波前像差RMS最大值在无指向偏差时为0.0139λ，有±5 nrad偏差时为0.0140λ，满足λ/40要求。此外，几何TTL耦合系数为0.243 pm/nrad，而总TTL耦合系数为0.122 pm/nrad，表明非几何TTL部分补偿了几何TTL效应。

**5.2 线性验证与工程可行性分析**

通过线性回归分析，光学路径长度信号（LPS）与光束倾斜角在俯仰方向呈良好线性关系（R2=0.992），TTL耦合系数拟合值为0.117 pm/nrad，验证了基于峰值（PV）值的评估一致性。基于蒙特卡洛方法（Monte Carlo analysis）的1000组公差样本分析表明：所有样本的TTL耦合系数均满足0.6 pm/nrad要求，99%分位值为0.275 pm/nrad，显示TTL性能对制造和装配误差具有高鲁棒性；但波前像差RMS仅93%样本满足λ/40要求，表明波前质量对误差更敏感，需在优化中进一步降低灵敏度。

**总结与讨论**

原文未设独立讨论部分，结论（Conclusion）部分对本研究进行了总结：本文提出了一种针对引力波探测星间激光干涉仪中望远镜光学参数的综合优化方法，结合贝叶斯算法与Zemax软件，以最小化TTL耦合系数并满足波前像差要求。优化后，望远镜在光束指向抖动10 nrad/Hz^1/2范围内，TTL耦合系数降至0.134 pm/nrad，波前像差低于λ/40，满足引力波探测中ILI的要求。结果表明所提出的望远镜光学参数优化方法具有显著性能。后续研究将致力于进一步降低由制造误差和装配偏差引起的TTL耦合噪声。