介观光学的理论方法.PDF

第三章介观光学的理论方法 由于介观体系的细微结构可以和波长相比拟，甚至比波 长小得多，在处理此类问题时，边界尤为复杂，可解析 求解的两类特例已不能满足许多实际问题的需要。 所以，目前发展出多种数值方法和理论方法，如： 纯数值计算：有限时域差分FDTD 有限元方法FEM 电偶极近似的算法：离散偶极近似DDA 耦合偶极子近似CDA 周期性结构中常用算法：多重散射MMS和平面波展开PWE 处理不规则纳米结构：格林函数方法GFM ， 转移矩阵方法TMM 1 古英 ygu@pku.edu.cn 2012年4 月 有限时域差分（FDTD ） Finite Difference Time Domain 在电磁领域用的最广泛的纯数值算法 －popular 能算任意电磁波段，任意结构，任意介电常数 －powerful 研究物理问题的手段：实验、解析和数值计算 FDTD是纯数值计算法，如格林函数和DDA都是半解析 半数值的算法，所以FDTD的适用性更强。原理上说， FDTD须用更多的计算资源和计算时间。目前快速发展 的计算机技术和计算方法可解决这一问题。 FDTD的计算结果需要清楚的头脑和物理思想去分析和 提炼。 2 2 FDTD 的实质就是数值解满足边界条件的maxwell 方程。 从maxwell 方程出发： ∇× =−∂ ∂ （1 ） E B / t ∇× =∂ ∂+ （2 ） H D / t J ∇⋅ ρ D ∇⋅B 0 and D εE B µH 加上初始和边条件的考虑后，只有（1 ），(2) 起作用 E 和H 可写成以下形式： ∇×E total =−µ∂H total / ∂−t σ* H total ∇×H total ε=∂E total / ∂+t σE total 3 Finite difference time domain method for elctromagnetics, Karl S. Kunz and Raymond J. Luebbers,1993 CRC press. E total ≡E incident +E scattered H total ≡H incident +H scattered 将E 和H 代入后得： inc inc ∇× =− ∂ ∂