离散偶极近似

✍ dations ◷ 2025-07-11 06:56:51 #计算科学,电动力学,散射

离散偶极近似(Discrete Dipole Approximation DDA)是一种用来求解物体散射电磁波的计算方法。它使用大量偶极子组成的阵列来模仿连续的物体,通过求解这些偶极子在入射电磁波照射下的极化度来获得物体吸收、散射电磁波的性质。

在电动力学的框架内,求解物体对电磁波的吸收、散射情况,实质上是计算物体内部和周围空间的电磁场分布。理论上,有介质存在情况下的电磁场分布可以通过求解麦克斯韦方程组获得。然而,由于麦克斯韦方程组的复杂性,它只能在具有独特对称性的体系中求得解析解。那么,对于一般形状的物体,通常采用数值方法近似求解其周围电磁场分布。离散偶极近似就是这样一种方法,它假设物体的电磁波散射特性是由其电子对于入射电磁波的反馈作用形成——电子在电磁波的作用下发生受迫振动,而与其正电荷中心分离形成振荡电偶极,它们在振动时能够辐射电磁波并作用于其它电偶极。进一步,设想物体是由大量的电偶极组成,则由电动力学理论可以建立起描述所有偶极子相互影响的线性方程组,求解该方程组获得偶极电磁场。最后,把所有偶极的电场作用叠加后就获得了整个物体内部以及周围空间的电磁场。

1964年,Howard DeVoe在其论文中建立了DDA方法的基本框架。DeVoe发展出一种经典物理模型,该模型可以由单体(比如分子)的光学性质出发求解出聚集体(比如分子晶体)的光学性质。DeVoe认为,如果求出聚集体的总偶极矩,就可以导出其折射率和消光系数等性质。为此,需要将所有的分子偶极矩叠加从而获得总偶极矩。然而,分子偶极矩是由入射电磁波诱导产生,同时这些分子偶极还会产生电场,进一步影响其它分子偶极。那么,分子偶极矩中就包含了入射电场和其它偶极电场两部分,由此可以建立起包含所有偶极矩的线性方程组加以求解。DeVoe奠定了由偶极间相互作用求解物体光学性质的方法。但是,该论文中描述偶极间电场作用时使用的是电偶极的静电场表达式,而非振荡偶极子的电场。

1973年,Purcell在其论文中确立了DDA方法的基本原理并计算了任意形状颗粒的光散射特性。他是在关于星际尘埃散射星光的研究中发展出这一方法的,因而并未引用DeVoe的研究成果。然而,相比较DeVoe旨在建立一种经典物理的模型来描述单体与聚集体光学性质之间的关系,Purcell则明确了DDA方法的一些基本概念,以及它用于计算颗粒光散射的用途。他首先提出使用3维的偶极子阵列模拟颗粒的散射行为,并且采用振荡偶极子的电磁波公式描述偶极子电场的影响,相比与DeVoe更为准确。其次,Purcell通过Clausius-Mossotti关系式确定了偶极阵列中偶极子极化率(polarizability) α {\displaystyle \alpha } 与颗粒材料的介电常数 ϵ {\displaystyle \epsilon } 的关系(DeVoe使用分子消光系数导出分子偶极极化率)。由此求出极化率就解决了DDA方法中的基本问题:偶极矩与局部电场之间的数值关系。这样就能够构建描述偶极矩之间关系的线性方程组 P i = α E i {\displaystyle {\boldsymbol {P}}_{i}=\alpha {\boldsymbol {E}}_{i}} 。最后,他设计了一个递推关系求解该方程组,并且根据电动力学原理求出颗粒的吸收系数和消光系数。

下面主要以Draine编写的软件DDSCAT为例,结合他的研究工作,介绍DDA方法的发展过程。

1988年,Draine在其论文中对DDA方法做出了几项重要改进(使用Fortran语言编写了计算程序,成为后来的DDSCAT)。首先,Draine认为Clausius-Mossotti关系式描述的是静电场环境下极化度与介电常数的关系,并不完全适用于电磁波条件。因此,他引入radiation reaction对其进行了修正。其次,通过分析计算结果关于偶极子阵列粒度的变化情况,他提出了DDA方法的误差表达式,确定了偶极子数量N的取值标准。最后,Draine采用复数共轭梯度算法迭代求解偶极矩方程组,获得更好的收敛性。

1991年,Goodman等人在论文中指出,当偶极子阵列具有空间周期性时,复数共轭梯度算法中的矩阵乘法实质上是卷积运算,因此可以使用快速傅立叶变换(Fast-Fourier Transform FFT)技术对其加速,节省运算时间。作者在其DDSCAT软件中加入了这一功能并对其时间复杂性进行了研究。

1993年,Draine与Goodman指出使用Lattice Dispersion Relation (LDR)关系描述偶极子极化度能够使DDA方法的解与mie散射的结果更为接近,因而取代了Clausius-Mossotti plus radiative reaction (CMRR)。

2008年,Draine与Flatau对DDA方法进行了改进,使其可以计算二维周期性结构或者一维无限长物体的散射,扩展了DDA方法的应用范围。

此外,Yurkin和Maltsev在其文章中分析了DDA方法的误差与偶极间距d的关系。作者从电磁场散射理论出发,推导出偶极离散化误差的表达式,同时分析了正方体偶极阵列和实际散射体形状差别造成的误差。

2007年,Penttila等人发表文章,从计算速度、内存消耗等方面比较了SIRRI、DDSCAT、ADDA、ZDD这4种DDA软件各自的优劣。

2007年,Yurkin和Hoekstra发表综述,对DDA方法进行了较为全面的总结,涵盖了基本原理,计算方法,求解技术的诸多方面的研究成果和发展。

实质上,DDA方法是电磁散射公式积分形式的离散化表达。

相关

  • 饥饿饥饿是动物的一种感受,产生的原因是肝的糖原水平下降到低于一个阈值。这种不快的感觉是在下丘脑产生,由肝和胃的感受器感应。一般人可以超过30天不进食而不至于死亡,但在缺水的
  • 自然资源自然资源,亦称天然资源,是指在其原始状态下就有价值的物产。一般来说假如获取这个物产的主要工程是收集和纯化,而不是生产的话,那么这个物产是一种自然资源。采矿、采油、渔业和
  • δ-3-蒈烯3-蒈烯(Carene, 或delta-3-carene)是一种天然存在的单萜烯,是松节油的组分之一。视来源不同,松节油中含量可达42%。3-蒈烯有强烈香气,不溶于水,与油脂混溶。3-蒈烯的天然来源包括
  • 钻石号A型运载火箭钻石号系列运载火箭(法语:Diamant 意为钻石),为法国第一个完全自行研制的火箭,是一不可重复使用的卫星发射系统且不由美国或苏联协助建造,钻石号系列运载火箭为其后的欧罗巴运载火
  • 仓鸮仓鸮(学名:Tyto alba)是分布最广的猫头鹰种类,也是所有鸟类中分布最广的之一,属于鸮形目中两科之一的草鸮科。它广布于世界各地,除了极地、荒漠、喜马拉雅山脉以北、印度尼西亚的
  • 武田信廉武田信廉(日语:武田 信廉,1532年-1582年),武田信虎三男,母亲为正室大井夫人。乳名孙六,又名刑部少辅,出家后号称武田逍遥轩、逍遥轩信廉,是武田二十四将中的一人。永禄四年(1561年)的第4
  • 椒盐椒盐为中国各个地区常见的调味料。可以用小中火将花椒粒与盐炒约一两分钟至花椒香气溢出,盛起待凉即可。在中国华北地区,椒盐的做法是将花椒炸酥,碾碎成粉,搀在细盐之中,一般作为
  • 特奥巴尔德·冯·贝特曼-霍尔韦格特奥巴尔德·冯·贝特曼-霍尔韦格(德语:Theobald von Bethmann-Hollweg,1856年11月29日-1921年1月1日),德国政治家,曾于1909年至1917年间任德意志帝国总理。特奥巴尔德·冯·贝特曼
  • 人情鬼蜮《人情鬼蜮》(英语:)是一部1935年的电影,主演是艾琳·麦克马洪(英语:Aline MacMahon)和贝锡·罗斯本。本片改编自爱德华·乔多洛夫(英语:Edward Chodorov)的同名戏剧,以及休·沃波尔的
  • 赵羾赵羾(1364年-1436年),字云翰,中书省晋宁路解州夏县(今山西省夏县)人,明朝政治人物,官至礼部尚书、兵部尚书。洪武年间,其因乡举进入太学,授兵部职方司主事。后上疏图天下要害厄塞,并屯戍