[论文解读] Propagation of Electromagnetic Waves on a Rectangular Lattice of Polarizable Points

D. Gutkowicz-Krusin, B. T. Draine|arXiv (Cornell University)|Mar 3, 2004

Electromagnetic Scattering and Analysis被引用 39

一句话总结

本文推导了极化点偶极子矩形晶格上电磁波传播的色散关系，将离散偶极近似（DDA）从立方晶格扩展至矩形晶格。提出了一种新的极化率设定方法，确保晶格能准确匹配连续介质的折射率，纠正了Draine & Goodman（1993）先前的工作，并通过散射计算验证了该方法的精度与原始方法相当。

ABSTRACT

We discuss the propagation of electromagnetic waves on a rectangular lattice of polarizable point dipoles. For wavelengths long compared to the lattice spacing, we obtain the dispersion relation in terms of the lattice spacing and the dipole polarizabilities. We also obtain the polarizabilities required for the lattice to have the same dispersion relation as a continuum medium of given refractive index m.; our result differs slightly from previous work by Draine & Goodman (1993). Our new prescription can be used to assign dipole polarizabilities when the discrete dipole approximation is used to study scattering by finite targets. Results are shown for selected cases.

研究动机与目标

• 将离散偶极近似（DDA）从立方晶格扩展至矩形晶格，以允许各向异性的晶格间距，从而更精确地模拟复杂目标几何结构。
• 在长波长极限（kd ≪ 1）下，推导极化偶极子矩形晶格上电磁波的色散关系。
• 提出一种新的偶极子极化率赋值方法，使得晶格能够再现给定折射率m的连续介质的色散关系。
• 纠正并改进Draine & Goodman（1993）针对立方晶格提出的极化率设定方法，现将其推广至矩形晶格。
• 通过DDA散射计算验证新方法，结果表明其精度与原始DG93方法相当。

提出的方法

• 在Lorentz规范下，建立包含具有张量极化率的晶格位置偶极子源的矢量势波动方程。
• 采用布洛赫波解描述无限晶格中的电磁场，假设谐波时间依赖性和周期性。
• 通过求解原点处偶极子极化率与矢量势的关系，推导出模态方程，同时考虑所有其他偶极子的贡献。
• 将色散关系按(kd)²的幂级数展开，获得O[(kd)²]阶的校正项。
• 通过反演色散关系，确定为使晶格再现指定折射率m的连续介质行为所需的关键偶极子极化率。
• 数值计算涉及倒易晶格矢量和指数截断的晶格和，以处理R₀(i)、R₁、R₂(i)和R₃(i,j)求和中的发散问题。

实验结果

研究问题

• RQ1 在长波长极限下，极化偶极子矩形晶格上电磁波传播的色散关系是什么？
• RQ2 如何在矩形晶格上分配偶极子极化率，使得晶格能模拟给定折射率的连续介质的波传播行为？
• RQ3 在散射计算中，与原始DG93方法相比，矩形晶格的新极化率设定方法在精度上表现如何？
• RQ4 对于不同各向异性的晶格轴比，关键系数R₀(i)、R₁、R₂(i)和R₃(i,j)的数值是多少？
• RQ5 新设定方法在有限目标（如球体）的DDA散射模拟中是否仍保持精度？

主要发现

• 新极化率设定方法确保矩形偶极子晶格能准确再现折射率为m的连续介质的正确色散关系，其结果与DG93方法存在差异。
• 对于立方晶格（d₁=d₂=d₃=1），系数R₀(1)=R₀(2)=R₀(3)=0且R₁=R₂(1)=R₂(2)=R₂(3)=0，与各向同性一致。
• 对于各向异性的晶格，如d₁:d₂:d₃ = 1:1.5:3，R₀(1)=1.38481，R₀(2)=-0.14304，R₀(3)=-1.24176，表明存在显著的方向依赖性。
• 当d₁:d₂:d₃=1:3:3时，R₃(1,1)达到2.04073，显示出张量响应中的显著各向异性。
• R₁的取值范围从-0.53869（1:1:1.5）到-5.47612（1:1:3），表明随着轴比增大，偏离各向同性的程度加剧。
• 使用新设定方法进行的DDA散射计算得到的截面精度与DG93方法相当，验证了其在实际模拟中的适用性。