[论文解读] Form factor (Fourier shape transform) of polygon and polyhedron

Joachim Wuttke|arXiv (Cornell University)|Mar 1, 2017

Optical measurement and interference techniques参考文献 16被引用 9

一句话总结

本论文提出了一种数值稳定的闭式解析方法，利用散度定理和斯托克斯定理计算任意多边形和多面体的傅里叶变换（形式因子），并通过级数展开显式处理可去奇点。该方法可在所有复杂波矢下实现高精度计算，包括靠近奇点的区域，并已在BornAgain GISAXS模拟软件中得到验证，相对精度优于10⁻¹⁰。

ABSTRACT

The Fourier transform of the indicator function of arbitrary polygons and polyhedra is computed for complex wavevectors. Using the divergence theorem and Stokes' theorem, closed expressions are obtained. Apparent singularities, all removable, are discussed in detail. Loss of precision due to cancellation near the singularities can be avoided by using series expansions.

研究动机与目标

• 推导三维空间中任意多边形和多面体形式因子的通用闭式解析表达式。
• 解决标准计算中导致数值不稳定和精度损失的形式因子表达式中的表观奇点。
• 开发一种适用于复波矢（包括具有虚部的波矢，与吸收和放大相关）的数值鲁棒算法。
• 确保在不同计算模式（解析 vs. 级数展开）之间保持高数值精度和连续性，以适用于小角散射模拟。
• 为先进GISAXS模拟软件（如BornAgain和IsGISAXS）提供精确且高效的形式因子计算基础。

提出的方法

• 应用散度定理和斯托克斯定理，将多边形和多面体指示函数的体积分和面积分转化为边界上的线积分和面积分。
• 通过参数化波矢分量和边几何的边上线积分，推导多边形形式因子的闭式表达式。
• 在奇点附近（如 q → 0 或 q 平行于法向）使用幂级数展开，避免抵消误差并保持双精度精度。
• 采用混合计算策略，根据波矢大小和方向在解析表达式与级数展开之间切换，使用经优化的阈值参数 c 以确保连续性。
• 通过对称性、特例和连续性测试，在多种粒子形状和波矢配置下验证该方法。
• 将算法集成至BornAgain GISAXS模拟软件中，确保与现有参考代码（如IsGISAXS）的兼容性。

实验结果

研究问题

• RQ1 如何利用向量微积分定理，以闭式表达计算任意多边形的形式因子，避免在奇点处的数值不稳定性？
• RQ2 多面体形式因子中的可去奇点具有何种结构？如何通过解析延拓实现数值稳定性？
• RQ3 如何利用级数展开在标准计算因抵消误差而失效的奇点附近保持双精度精度？
• RQ4 在何种策略下切换解析与基于级数的形式因子计算，可确保所有波矢下的连续性和精度？
• RQ5 所推导方法在多大程度上可与已知解析结果及参考模拟代码（如IsGISAXS和BornAgain）进行验证？

主要发现

• 利用散度定理和斯托克斯定理，任何多边形或多面体的形式因子均可通过闭式表达计算，问题被简化为边界积分。
• 形式因子表达式中的所有表观奇点均为可去奇点，可通过精心构造的幂级数展开处理。
• 在对称性和特例测试中，最坏情况下的相对偏差小于5×10⁻¹⁰，表明具有极佳的数值一致性。
• 解析与基于级数的计算之间的连续性保持良好，相对偏差 δ_cont < 8×10⁻¹⁶，确保在所有波矢区域的鲁棒性。
• 该算法已成功实现于BornAgain GISAXS模拟软件中，并通过与解析积分及参考代码的对比得到验证。
• 该方法可对复杂粒子形状（如平截棱锥、十二面体、截角立方体）在所有波矢方向和大小下实现高精度形式因子计算。