Symmetry-Based Design Rules for Second-Harmonic Generation in Stacked and Twisted MoS2 Bilayers

技术摘要：基于对称性的堆叠与扭转 MoS$_2$ 双层二阶谐波产生设计规则

问题陈述
二维过渡金属硫族化合物（TMDs），特别是 MoS$_2$ 的非线性光学响应对堆叠构型和层间旋转高度敏感。虽然第一性原理计算可以模拟理想化系统，但在处理具有实际低对称性特征的大规模莫尔超晶格时，计算成本会变得极高。因此，需要一个严谨的、基于对称性的框架，能够在不依赖于针对每个特定扭转角或堆叠变体进行暴力计算的情况下，将结构构型映射到其二阶极化率张量（$\chi^{(2)}$）上。

方法论
作者建立了一个结合群论分析与针对性第一性原理验证的综合框架：

1. 对称性分析： 利用诺曼对称性原理（Neumann's principle），作者推导了适用于 MoS$_2$ 双层的四种不同点群中允许和禁止的 $\chi^{(2)}$ 张量分量，分别为：$D_{3h}$（单层及 AA 堆叠）、$D_{3d}$（AB(2H) 堆叠）、$C_{3v}$（AB(3R) 堆叠）以及 $C_3$（扭转双层）。他们分析了特定的对称元素（反演中心 $I$、水平镜像 $\sigma_h$、垂直镜像 $\sigma_v$）如何约束张量结构。
2. 第一性原理计算： 通过使用带有 DFT-D3 范德华力修正的 VASP 码中的 PBE 泛函进行的密度泛函理论（DFT）计算来验证理论预测。在独立粒子近似下，利用极大局域化 Wannier 函数（通过 Wannier90 和 postw90）计算非线性光学响应。计算过程包含了自旋轨道耦合，并将响应分解为带间（interband）和带内（intraband）贡献。
3. 角响应建模： 作者推导了各对称类下共偏振二阶谐波产生（SHG）强度 $I_{\parallel}(\theta)$ 随入射偏振角 $\theta$ 变化的解析表达式，特别研究了在扭转构型中相位移动的出现。

核心贡献

• 对称性-响应映射： 本文提供了关于堆叠构型（AA、AB(2H)、AB(3R）及扭转）如何决定 $D_{3h}$、$D_{3d}$、$C_{3v}$ 和 $C_3$ 点群下 $\chi^{(2)}$ 张量结构的确定性映射。
• SHG 抑制与激活的设计规则：
  • 反演对称性（$D_{3d}$）： AB(2H) 堆叠中反演中心的出现严格强制 $\chi^{(2)} = 0$，导致通过破坏性层间干涉实现完全的 SHG 抑制。
  • 水平镜像破缺（$C_{3v}$ 与 $C_3$）： 在 AB(3R) 和扭转双层中，通过打破 $\sigma_h$ 对称性，激活了在 $D_{3h}$ 系统中被禁止的面外张量分量（$\chi_{zzz}^{(2)}, \chi_{zxx}^{(2)}, \chi_{xxz}^{(2)}$）。
  • 垂直镜像破缺（$C_3$）： 在扭转双层中，移除垂直镜像面会激活面内分量 $\chi_{xxx}^{(2)}$，该分量在更高对称性的构型中严格为零。
• SHG 瓣状图的几何锁定： 作者证明了在扭转双层中，$\chi_{xxx}^{(2)}$ 的激活会导致 SHG 极化瓣状图发生刚性方位角旋转。至关重要的是，他们发现这种旋转角 $\phi_0$ 在整个光学光谱范围内是频率无关的，并且通过关系式 $\theta_{tw} = 2\phi_0$ 与结构扭转角 $\theta_{tw}$ 几何锁定。

结果

• 张量验证： 第一性原理计算证实了群论预测。AB(2H) 双层表现出消失的 SHG 信号。AB(3R) 双层恢复了类似于单层的面内响应，同时获得了显著的面外分量。扭转双层（$C_3$）同时表现出非零的 $\chi_{yyy}^{(2)}$ 和 $\chi_{xxx}^{(2)}$。
• 光谱特征： SHG 光谱由对应于 A 和 B 激子跃迁（自旋轨道分裂价带）以及与能带嵌套跃迁相关的更高能量 C 峰的双光子共振主导。
• 极化模式旋转： 对于扭转角 $\theta_{tw} = 21.8^\circ$ 的扭转双层，SHG 极化模式旋转了 $\phi_0 \approx 10.9^\circ$。这种偏移在整个计算频率范围内保持恒定，证实了 $\chi_{xxx}^{(2)}$ 与 $\chi_{yyy}^{(2)}$ 之间的相位关系是由结构几何而非动态光谱变化决定的。
• 手性分量： 虽然手性面外分量 $\chi_{xyz}^{(2)}$ 在 $C_3$ 对称性下是允许的，但由于层间耦合较弱且单个单层内局部保留了镜像对称性，其量级被发现极小。

意义
本文建立了一种直接的、与波长无关的二维材料扭转角测定方法。通过将结构扭转角直接与 SHG 极化瓣状图的旋转（$\theta_{tw} = 2\phi_0$）联系起来，作者提供了一种非破坏性、高通量的工具，用于表征莫尔超晶格和堆叠畴。此外，这项工作为仅通过堆叠控制来工程化二维材料非线性光学响应提供了清晰的路线图，从而能够设计具有定制功能的异质结构，例如完全的 SHG 抑制或特定的张量分量激活。该框架被认为是通用的，适用于 2H-TMD 双层及类似的范德华晶体，为未来高通量筛选非线性光学架构提供了基础。