[论文解读] User Guide for the Discrete Dipole Approximation Code DDSCAT 6.1

B. T. Draine, Piotr J. Flatau|ArXiv.org|Sep 10, 2004

ZnO doping and properties参考文献 1被引用 31

一句话总结

本论文介绍了DDSCAT 7.0，这是一款免费的开源Fortran-90代码，采用离散偶极近似（DDA）计算任意形状粒子及周期性纳米结构的电磁散射与吸收。该软件可对尺寸参数高达25且折射率不远离单位值的目标实现高精度模拟，支持MPI、OpenMP、MKL以及通过DDfield后处理程序计算近场E/B场等高级功能。

ABSTRACT

DDSCAT 6.1 is a software package which applies the discrete dipole approximation (DDA) to calculate scattering and absorption of electromagnetic waves by targets with arbitrary geometries and complex refractive index. DDSCAT 6.1 allows accurate calculations of electromagnetic scattering from targets with size parameters 2 pi a_eff/lambda < 15 provided the refractive index m is not large compared to unity (|m-1| < 2). DDSCAT 6.1 includes support for MPI and FFTW. We also make available a "plain" distribution of DDSCAT 6.1 that does not include support for MPI, FFTW, or netCDF, but is much simpler to install than the full distribution. The DDSCAT package is written in Fortran and is highly portable. The program supports calculations for a variety of target geometries (e.g., ellipsoids, regular tetrahedra, rectangular solids, finite cylinders, hexagonal prisms, etc.). Target materials may be both inhomogeneous and anisotropic. It is straightforward for the user to import arbitrary target geometries into the code, and relatively straightforward to add new target generation capability to the package. DDSCAT automatically calculates total cross sections for absorption and scattering and selected elements of the Mueller scattering intensity matrix for specified orientation of the target relative to the incident wave, and for specified scattering directions.

研究动机与目标

• 为DDSCAT 7.0提供全面的用户指南，该开源软件包实现了离散偶极近似（DDA）以进行电磁散射计算。
• 使研究人员能够对尺寸参数≤25的任意形状粒子（包括非均匀和各向异性材料）实现高精度的散射与吸收建模。
• 通过支持MPI、OpenMP和Intel MKL等高级计算特性，实现对一维和二维周期性阵列及近场电磁场的模拟。
• 通过取向平均和场后处理，实现散射截面、穆勒矩阵元素及辐射力的精确计算。
• 为用户提供标准化且可扩展的框架，以导入自定义几何结构，并通过新增目标生成例程扩展代码功能。

提出的方法

• 离散偶极近似（DDA）将目标离散化为极化偶极子晶格，利用共轭梯度法等迭代求解器求解耦合偶极子方程。
• 该代码通过快速傅里叶变换（FFT）高效实现迭代求解器中的矩阵-向量乘法，支持多种FFT库（如FFTW和Intel MKL）。
• 对于周期性结构，代码采用Ewald求和技术施加周期性边界条件（PBC），以处理长程偶极-偶极相互作用。
• 用户通过参数文件（ddscat.par）定义目标几何形状与组成，指定偶极子位置、材料属性、取向及入射波参数。
• 通过远场与近场计算相结合，代码计算散射与吸收截面、穆勒矩阵元素及辐射力。
• 通过DDfield代码实现后处理，可在目标附近用户定义的位置计算电场与磁场，支持详细的场分布分析。

实验结果

研究问题

• RQ1 如何在并行、可移植且可扩展的软件包中高效实现适用于任意粒子几何形状的离散偶极近似？
• RQ2 在现代高性能计算（HPC）架构上，求解大规模偶极子相互作用系统时，最优的计算策略是什么，以实现高精度与高性能？
• RQ3 如何将DDA扩展以模拟纳米结构周期性阵列，并计算近场电磁增强效应？
• RQ4 对于具有复杂折射率和各向异性材料的目标，DDA的精度与适用性边界是什么？
• RQ5 用户如何可靠地生成并导入复杂目标几何形状（包括多层与复合结构）至该模拟框架？

主要发现

• DDSCAT 7.0可对尺寸参数高达25且折射率满足|m−1| ≤ 2的目标实现精确的电磁散射与吸收计算。
• 该代码支持通过MPI和OpenMP实现并行执行，并利用Intel MKL优化线性代数运算，在多核与集群系统上显著提升性能。
• DDfield后处理程序的引入使用户能够计算目标附近任意用户定义点的电场与磁场，支持详细的近场分析。
• 该软件通过内置生成器和用户自定义输入文件，支持广泛的几何形状，包括椭球体、圆柱体、棱柱体、多层薄膜及复杂复合结构。
• 参数文件ddscat.par支持灵活配置入射波偏振、散射方向及输出选项，包括对随机取向粒子的取向平均。
• 该代码已通过与解析解及已发表结果的对比完成验证，已知问题已在官方DDSCAT网站上记录并更新。