[论文解读] Generalized Drude-Lorentz Model Complying with the Singularity Expansion Method

Isam Ben Soltane, Félice Dierick+2 more|arXiv (Cornell University)|Jan 11, 2024

Photonic and Optical Devices被引用 4

一句话总结

本论文提出了一种基于奇点展开法（Singularity Expansion Method, SEM）推导的广义德鲁德-洛伦兹（Generalized Drude-Lorentz, GDL）模型，将介电常数视为复频率平面上的亚纯函数。通过利用PyTorch的自动微分技术，该方法能够从实验数据中高精度地反演金属、电介质和二维材料的模型参数，其极点与留数揭示了材料在亚电子伏特（sub-eV）精度下的微观行为特征。

ABSTRACT

Deriving analytical expressions of dielectric permittivities is required for numerical and physical modeling of optical systems and the soar of non-hermitian photonics motivates their prolongation in the complex plane. Analytical models are based on the association of microscopic models to describe macroscopic effects. However, the question is to know whether the resulting Debye Drude Lorentz models are not too restrictive. Here we show that the permittivity must be treated as a meromorphic transfer function that complies with the requirements of complex analysis. This function can be naturally expanded on a set of complex singularities. This singularity expansion of the dielectric permittivity allows us to derive a generalized expression of the Debye Drude Lorentz model that complies with the requirements of complex analysis and the constraints of physical systems. We show that the complex singularities and other parameters of this generalized expression can be retrieved from experimental data acquired along the real frequency axis. The accuracy of this expression is assessed for a wide range of materials including metals, 2D materials and dielectrics, and we show how the distribution of the retrieved poles helps in characterizing the materials.

研究动机与目标

• 为克服经典德鲁德-洛伦兹模型在描述复频率平面上介电常数时的局限性，特别是在非厄米光子学与时域仿真中的不足。
• 建立一个严格的理论框架，将介电常数视为符合复分析与物理因果律的亚纯传递函数。
• 开发一种基于现代优化工具的参数反演方法，能够从实轴上的实验数据中高精度重构介电常数。
• 证明GDL模型中复极点与留数的分布编码了材料在不同频率区域响应的有意义物理信息。

提出的方法

• 在复频率平面上将介电常数形式化为亚纯传递函数，从而通过奇点展开法（SEM）实现精确的奇点展开。
• 将SEM展开式重新表述为包含频率相关虚部留数的广义德鲁德-洛伦兹（GDL）模型，使其与经典洛伦兹模型相区别。
• 引入具有复留数的广义洛伦兹项，以在时域中体现电场一阶导数的影响，从而实现精确的时域建模。
• 采用基于优化的参数反演方法，利用PyTorch实现自动微分，最小化实验数据与GDL模型预测之间的误差。
• 将该方法应用于9种材料（金属、氧化物、二维材料），在宽光谱范围（紫外至近红外）内使用最少数量的极点拟合GDL模型。
• 通过分析极点与留数的分布，重建物理机理，将其与金属、电介质或近似恒定介电常数行为等材料区域关联起来。

实验结果

研究问题

• RQ1 奇点展开是否能为介电常数提供比经典德鲁德-洛伦兹模型更通用且物理上一致的框架？
• RQ2 如何从实频率的实验数据中高精度地反演介电常数函数的复极点与留数？
• RQ3 广义洛伦兹项中留数的虚部具有何种物理意义，其对时域行为有何影响？
• RQ4 复平面上极点的分布在多大程度上能够表征金属、电介质与透明区域之间的转变？
• RQ5 是否可用一组极少数的广义极点准确描述宽光谱范围内的介电常数，包括二硫化钼等复杂材料与氧化物？

主要发现

• GDL模型在1–10 eV光谱范围内对金属（如金）、电介质（如SiO₂）和二维材料（如hBN）的实验介电常数数据拟合精度达到亚电子伏特（sub-eV）级别。
• 对于金（Au），反演得到的德鲁德极点为纯虚数且靠近原点，与自由电子行为一致；而位于2.7–13 eV的洛伦兹极点与已知带隙位置相符。
• 广义洛伦兹项包含复留数，引入了频率相关的虚部贡献，从而通过电场的一阶导数实现精确的时域建模。
• 极点的分布揭示了不同的物理区域：德鲁德极点对应金属行为，洛伦兹极点表示共振跃迁，而实轴附近无极点则对应近似恒定介电常数行为。
• 基于自动微分的反演方法仅使用5–7个极点即实现低误差，展现出高效率与强鲁棒性。
• 该方法成功表明，GDL模型中的德拜型项被纳入广义洛伦兹框架，其留数携带与弛豫动力学相关的物理意义。