[论文解读] T-matrix calculation via discrete-dipole approximation, point matching and exploiting symmetry

Vincent L. Y. Loke, Timo A. Nieminen+2 more|ArXiv.org|Dec 10, 2008

Orbital Angular Momentum in Optics参考文献 19被引用 7

一句话总结

本论文提出了一种基于离散偶极近似（DDA）并结合旋转对称性和镜像对称性利用的优化T-矩阵计算方法，使微米尺度物体（如微马达）的内存占用和计算时间减少达两个数量级。通过压缩相互作用矩阵并利用矢量球面波函数中的模态冗余，该方法实现了对具有高对称性的任意形状粒子的高效散射计算。

ABSTRACT

We present a method of incorporating the discrete dipole approximation (DDA) method with the point matching method to formulate the T-matrix for modeling arbitrarily shaped micro-sized objects. The \emph{T}-matrix elements are calculated using point matching between fields calculated using vector spherical wave functions and DDA. When applied to microrotors, their discrete rotational and mirror symmetries can be exploited to reduce memory usage and calculation time by orders of magnitude; a number of optimization methods can be employed based on the knowledge of the relationship between the azimuthal mode and phase at each discrete rotational point, and mode redundancy from Floquet's theorem. A 'reduced-mode' T-matrix can also be calculated if the illumination conditions are known.

研究动机与目标

• 开发基于离散偶极近似（DDA）的高效T-矩阵计算技术，用于任意形状的微米级粒子。
• 降低具有旋转对称性的微结构（如微马达）在T-矩阵模拟中的计算成本与内存占用。
• 利用离散旋转对称性和镜像对称性压缩相互作用矩阵，最小化冗余计算。
• 实现对尺寸超过10–20倍波长的粒子的精确且可扩展的散射模拟。
• 结合近场与远场点匹配方法，利用矢量球面波函数（VSWFs）实现稳健的T-矩阵构建。

提出的方法

• 通过利用离散旋转对称性压缩DDA相互作用矩阵，使其尺寸缩小 $ Q^2 $ 倍，其中 $ Q $ 为对称阶数。
• 使用相位因子 $ \text{exp}(im\theta) $ 和旋转矩阵，关联对称单元之间的偶极矩，避免完整矩阵重建。
• 利用正交矩阵 $ S $ 和 $ C_q $ 在笛卡尔坐标系与球坐标系之间进行坐标变换，实现在球坐标系中矢量场的对齐。
• 仅使用一个旋转单元构建简化线性方程组，求解偶极矩并利用对称性传播结果。
• 利用镜像对称性进一步将计算负载降低4倍，利用VSWFs在 $ xy $-平面上反射时的偶/奇对称性。
• 实现近场与远场点匹配，结合VSWFs计算T-矩阵元素 $ p_{nm} $ 和 $ q_{nm} $，通过相位与振幅匹配提取散射系数。

实验结果

研究问题

• RQ1 如何在不损失精度的前提下，降低大型对称微结构的T-矩阵计算成本？
• RQ2 在微马达中，旋转对称性和镜像对称性在多大程度上可用于压缩DDA相互作用矩阵？
• RQ3 VSWF展开中的模态冗余对对称粒子T-矩阵构建效率有何影响？
• RQ4 近场与远场匹配策略在T-矩阵提取中的精度与计算效率方面有何差异？
• RQ5 DDA、对称性利用与模态冗余的结合是否能够实现尺寸超过10–20倍波长粒子的可扩展模拟？

主要发现

• 相互作用矩阵被压缩了 $ Q^2 $ 倍，其中 $ Q $ 为旋转对称阶数，显著降低了内存占用与计算时间。
• 对于具有四阶旋转对称性的粒子（如立方体微马达），该方法使计算时间与内存使用减少达两个数量级。
• 镜像对称性进一步将计算负载降低4倍，实现了对称微结构的高效模拟。
• 通过识别四重对称粒子中仅 $ m = 0, 4, 8, \text{等。} $ 模态有贡献，有效利用了模态冗余，减少了冗余计算。
• 近场与远场匹配方法得到的T-矩阵元素一致，验证了该方法在不同场表示下的可靠性。
• 该方法在仅使用 $ \frac{\theta}{5} $ 分辨率下即实现了对十字形转子的力矩计算收敛，与参考理论结果高度一致。