Structural Chirality and Natural Optical Activity across the $α$-to-$β$ Phase Transition in SiO$_2$ and AlPO$_4$ from first-principles

技术摘要：SiO$_2$ 与 AlPO$_4$ 在 $\alpha$-$\beta$ 相变过程中的结构手性和天然光学活性

问题陈述
天然光学活性（NOA），即平面偏振光的旋转，传统上与结构手性相关。然而，文献中存在一种持久的概念性简化：即假设旋光符号直接由空间群符号的手性（例如 $P6_422$ 与 $P6_222$）决定。这种关联在特定情况下失效，例如石英（SiO$_2$）和白熊矿（AlPO$_4$）的 $\alpha$-$\beta$ 相变。在这些转变中，对称性从高对称性的六方相（$P6_422$ 或 $P6_222$）变为低对称性的三方相（$P3_121$ 或 $P3_221$），导致主螺旋轴的标称手性发生反转。矛盾的是，旋光符号在这一转变过程中保持不变。这种行为的精确微观起源以及决定 SiO$_2$ 和 AlPO$_4$ 系统中 NOA 符号的因素，需要超越简单的空间群分配进行澄清。

方法论
作者利用在 Abinit 代码中实现的密度泛函理论（DFT）和密度泛函微扰理论（DFPT）进行了第一性原理研究。

• 结构建模： 研究采用了带有优化规范赝势（PseudoDojo v.4）的 PBEsol 交换相关泛函，并使用 40 Hartree 的能量截断。结构弛豫采用 $6\times6\times6$ k-网格进行。
• 声子分析： 利用 DFPT 识别了高对称相在布里渊区中心（$\Gamma$ 点）的不稳定声子模。通过凝聚不稳定模的特征向量来模拟相变路径。
• 光学活性计算： 研究重点关注低频极限下的旋光率 $\bar{\rho}(\omega) = \rho(\omega)/(\hbar\omega)^2$，该值是基于 Zabalo 和 Stengel 开发的基于 DFPT 的线性响应方法计算得出的。计算假设光沿光轴（z 方向）传播。
• 对比分析： 为了隔离特定原子种类的贡献，研究进行了一项计算实验：将 SiO$_2$ 中的 Si 原子替换到 AlPO$_4$ 结构的 Al 和 P 位点中，以观察光学响应的变化。

关键结果

1. 相变机制： SiO$_2$ 和 AlPO$_4$ 的 $\alpha$-$\beta$ 转变是由中心 $\Gamma_3$ 声子模的凝聚驱动的。该模式是沉默模（既不是拉曼活性也不是红外活性），并诱导了一种位移相变。这种畸变破坏了高对称相的 $6_4$ 螺旋轴对称性，从而在低对称相中产生 $3_1$ 螺旋轴。
2. 空间群与旋光符号之间的差异：
   • 对于 SiO$_2$，计算预测高对称相（$P6_422$）和低对称相（$P3_121$）均为正旋光率。
   • 对于 AlPO$_4$，计算预测两个相位均为负旋光率。
   • 至关重要的是，尽管在转变过程中螺旋轴的手性发生了反转（例如从 $P6_422$ 到 $P3_121$），但每种材料的旋光符号保持不变。
3. 光学活性符号的起源： 研究表明，NOA 的符号并非由空间群螺旋轴的标称手性决定，而是由结构中极化率最高的原子的原子尺度手性决定的。
   • 在 SiO$_2$ 中，氧原子形成螺旋链，决定了正旋光。
   • 在 AlPO$_4$ 中，Al 和 P 原子的有序化使单胞沿 c 轴加倍，改变了氧原子的相对位置。这导致了氧位点的左手螺旋排列（类似于 $P6_222$ SiO$_2$ 结构），从而决定了负旋光。
   • 替换实验证实，Si、Al 和 P 原子对光学活性的贡献极小；信号主要由氧子晶格主导。
4. 旋光率的演化： 随着 $\Gamma_3$ 畸变幅度从高对称相增加到低对称相，旋光率的量值线性下降。作者指出，虽然进一步的畸变在理论上可能使符号反转，但由于朗道能显著增加，此类构型在物理上是不可达的。

意义与主张
本文声称解决了 SiO$_2$ 和 AlPO$_4$ 中结构对称性变化与保持旋光符号之间看似矛盾的问题。其主要贡献在于澄清了天然光学活性的符号受控于最具极化性的原子种类（这些氧化物中的氧）的手性，而非全局空间群符号。

作者强调，虽然结构手性是 NOA 的充分条件，但不能仅凭空间群螺旋轴的手性来推断具体的旋光符号。这项工作是对仅根据对称性进行光学活性简单解释的一种警示，凸显了分析极化原子微观排列对于理解手性晶体光学性质的必要性。研究结果与 Glazer 和 Stadnicka 关于尽管螺旋轴反转但旋光符号保持不变的先前实验观察结果相一致，并为之提供了第一性原理的微观基础。