Dielectric function in WSe2

技术摘要：WSe2 单层维格纳晶体中的介电函数

问题陈述
本文旨在对二维（2D）过渡金属二硫化物（TMD）单层（特别关注单层 WSe2）中的维格纳晶体态进行理论表征。虽然在低密度电子气中存在维格纳结晶已被理论证实，且近期在 TMD 异质结构中得到了实验观察，但仍需要一个微观框架来计算这些晶态的电子能带结构和静态介电响应。其挑战在于，如何在连续平移对称性破缺的状态下，准确模拟动能、库仑相互作用（特别是与 2D 材料相关的 Keldysh 相互作用）以及交换效应之间的相互作用。

方法论
作者开发了一种自洽哈特里-福克（Hartree-Fock, HF）数值方法，用于计算全自旋极化 2D 维格纳晶体的能带结构和静态介电函数。其核心方法组成部分包括：

• 哈特里-福克公式化： 计算是在限制性 HF 框架内进行的，其中假设基态是完全自旋极化的。单粒子轨道使用平面波基组 $|k + G\rangle$ 以布洛赫形式进行展开，其中 $k$ 属于新生成的三角形维格纳晶体第一布里渊区，$G$ 代表倒格矢。
• 相互作用势： 电荷载流子之间的有效静电相互作用采用 Keldhes 相互作用进行建模，该相互作用考虑了周围介质的介电屏蔽以及材料的 2D 特性。
• 自洽循环： 算法从周期性高斯型电荷密度开始初始化。通过迭代更新哈特里势（由密度谐波导出）和福克交换矩阵（由占据本征态导出），直到收敛。密度谐波通过更新后的本征态进行重构，以闭合循环。
• 数值实现：
  • 布里渊区采样利用目标动量的不可约 $\Gamma-K-M$ 楔形区域，而交换项通过由对称性图像重构的全区求和进行评估。
  • 倒空间中的圆形截断法用于选择平面波基组，从而保持三角形晶格的 $D_6$ 点群对称性。
  • 通过在倒空间的小圆形单元内进行解析平均，处理了库仑/Keldysh 势在 $q=0$ 处的奇异性。
• 介电响应计算： 使用适用于绝缘、布里渊区折叠性质的维格纳晶体的 Lindhard 型极化率公式来计算静态介电函数 $\epsilon(q, 0)$。这涉及在占据的 HF 能带与空能带之间的带间跃迁中进行求和，并利用收敛后的本征态和特征值。

主要贡献

1. 开发专用代码： 作者展示了一种能够求解具有 Keldhes 相互作用的 2D 维格纳晶体自洽 HF 方程的数值代码，该代码专门针对 TMD 单层进行了定制。
2. 微观介电响应： 本研究提供了一种直接从 HF 能带结构和特征函数推导静态介电函数的方法，而非依赖于经典的声子模型或均匀电子气近似。
3. 晶体序表征： 研究引入并使用了特定的指标来量化结晶程度：
   • 第一星序参数 ($M_1$)，测量基本密度调制的振幅。
   • 密度对比度比率 ($C_n$)，量化实空间中密度极大值与极小值之间的分离。
   • 高对称点（$\Gamma, K, M$）处能带结构的演化，作为六角形布里渊区对称性的光谱特征。

结果
作者将其方法应用于不同密度参数（$r_s$）范围内的单层 WSe2，代表了从弱结晶到强局域态的转变过程：

• 能带结构演化： 随着 $r_s$ 增加（密度降低），系统从弱调制态转变为稳健的维格纳晶体。在高 $r_s$ 下，能带结构表现出与六角形布里渊区对称性一致的清晰简并模式，且能带沿 $\Gamma-K-M$ 路径平滑演化。
• 密度调制： 序参数 $M_1$ 随 $r_s$ 增加（例如从 $r_s=4$ 时的 0.100 到 $r_s=50$ 时的 0.586），表明电荷局域化增强。密度对比度比率 $C_n$ 也显著增加，证实了晶格位点处清晰电荷密度峰的形成。
• 介电函数： 静态介电函数 $\epsilon_{HF}(q, 0)$ 的峰值随 $r_s$ 增加而呈现系统性下降。这反映了电子极化率的降低；随着电子更加局域化，它们对外部扰动的屏蔽效率降低，导致屏蔽相互作用趋向于裸 Keldhes 相互作用。
• 收敛性： 自洽解表现出稳定性，最终的密度分布从初始的高斯猜测值弛豫到由动能、哈特里能和交换能平衡决定的状态。

意义与主张
本文将其定位为研究低密度 2D 系统中屏蔽、能带结构重构和相互作用效应的理论工具。作者声称，其方法提供了一个关于从弱结晶解到局域化维格纳晶体态演化的连贯微观图景。

这项工作的意义是在近期原子级薄层半导体取得进展的背景下提出的。作者建议，其计算的自旋守恒、冻结自旋/谷扇区可以通过对可调控栅极的 WSe2 单层进行时间分辨、螺旋度分辨的光谱学实验进行测试，该单层可通过光泵浦驱动进入极化态。然而，作者明确指出，其结果是作为晶体特征的“相位敏感指标”，而非热力学意义上的液-晶相变证明，并指出若要进行明确的相变鉴定，则需要能量比较或独立的相稳定性分析。该工作将传统的低密度 2D 维格纳晶体理论与近期关于 TMD 中电荷序的光学研究联系起来，提供了一个适配于单层 TMD 特有的 Keldhes 相互作用和光学响应的框架。