[論文レビュー] 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行列計算手法を提示する。微小な物体（例：マイクロローター）に対して、メモリ使用量および計算時間を最大2桁削減する。相互作用行列の圧縮とベクトル球面波関数におけるモードの冗長性を活用することで、任意の形状の高対称性粒子に対する効率的な散乱計算が可能になる。

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 $ を用いてデカルト座標系と球座標系との間の座標変換を実施し、球座標系におけるベクトル場の整合性を確保する。
• 1つの回転単位のみを用いて縮小された連立一次方程式系を構築し、対称性を介して双極子モーメントを計算・伝搬する。
• 鏡対称性を活用して計算負荷をさらに4倍低減し、$ xy $ 平面における反射におけるVSWFsの偶奇性（偶関数・奇関数）を応用する。
• 近場および遠方場の点一致をVSWFsを用いて実装し、T行列要素 $ p_{nm} $ および $ q_{nm} $ を計算する。位相および振幅一致を用いて散乱係数を抽出する。

実験結果

リサーチクエスチョン

• RQ1 大規模で対称的なマイクロオブジェクトに対するT行列計算の計算コストを、正確性を損なわずにどのように低減できるか。
• RQ2 マイクロローターに見られる回転対称性および鏡対称性を、DDA相互作用行列の圧縮にどの程度活用できるか。
• RQ3 VSWF展開におけるモードの冗長性が、対称性を持つ粒子に対するT行列構築の効率性に与える影響は何か。
• RQ4 近場および遠方場の一致戦略は、T行列抽出において正確性と計算効率の点でどのように比較できるか。
• RQ5 DDA、対称性の活用、およびモードの冗長性の組み合わせにより、10–20波長を超えるサイズの粒子に対するスケーラブルなシミュレーションが可能になるか。

主な発見

• 相互作用行列は $ Q^2 $ 倍に圧縮され、ここで $ Q $ は回転対称性の順序を表す。これにより、メモリ使用量および計算時間の大幅な削減が達成された。
• 4回回転対称性（例：立方体型マイクロローター）を有する粒子に対して、本手法は計算時間およびメモリ使用量を最大2桁削減した。
• 鏡対称性により、計算負荷がさらに4倍低減され、対称的マイクロ構造物の効率的シミュレーションが可能になった。
• モードの冗長性は、4重対称性を持つ粒子では $ m = 0, 4, 8, \text{etc.} $ のみが寄与することに着目することで活用され、重複計算が削減された。
• 遠方場および近場一致手法は一貫したT行列要素を生成し、異なる場の表現においても本手法の妥当性が検証された。
• クロス形状ローターのトルク計算において、解像度を $ \frac{\theta}{5} $ にまで低下させても収束が達成され、参照理論結果と高い精度で一致した。