Entropy-Driven Structural Phase Transition in Nb$_3$Cl$_8$ via Density Functional Theory and an Effective Model

技术摘要：通过 DFT 与有效模型研究 Nb$_3$Cl$_8$ 的熵驱动结构相变

问题陈述
Nb$_3$Cl$_8$ 是一个典型的平带簇莫特绝缘体，存在于有效的三角晶格上，被认为是候选量子自旋液体（QSL）态的潜在材料。然而，其低温物理特性被约 90 K 时发生的首次序结构相变所掩盖。该转变将高温顺磁 $\alpha$ 相转变为低温非磁性 $\beta$ 相，从而抑制了候选 QSL 区间的出现。该转变的微观起源仍存在争议，竞争性解释包括电荷离聚（charge disproportionation）或层间二聚单态形成。此外，热力学驱动力以及外部刺激（如单轴压力）如何抑制该转变的机制尚未得到充分理解。

方法论
作者采用了一个统一的自由能框架，结合了第一性原理密度泛函理论（DFT）与扩展哈伯德模型（extended Hubbard model），以解构晶格、电子和自旋的贡献。

1. DFT 计算： 使用非自旋极化 DFT 计算（使用带有 optB88-vdW 和 PBE 泛函的 VASP）来获取 $\alpha$（空间群 $P\bar{3}m1$）和 $\beta$（以 $R\bar{3}m$ 表示）两相的弛豫结构、基态能量和声子谱。构建了两个结构之间的线性插值路径以估算能垒。
2. 声子热力学： 在谐振近似内，利用密度泛函微扰理论（DFPT）计算声子自由能（$F_{\text{phonon}}$）和熵。
3. 有效模型： 为了处理半充满 2$a_1$ 流形中的莫特物理，在基于簇中心轨道的 Wannier 基组下构建了一个扩展哈伯д模型。参数（跳跃项 $t$、在位相互作用 $U$ 以及格点间库仑相互作用 $V$）是通过约束随机相位近似（cRPA）推导得出的。
4. 双计数修正： 实施了哈特里级（Hartree-level）双计数修正（$E_{\text{DC}}$），以减去已被 DFT 参考系捕获的低能 2$a_1$ 贡献，从而确保总自由能（$F_{\text{tot}}$）的一致性。
5. 自旋热力学： 将扩展哈伯德模型简化为有效的海森堡模型。对于 $\alpha$ 相，使用连通簇高温级数展开（HTSE）来计算自旋熵。对于 $\beta$ 相，由于存在强层间耦合，采用了简化的孤立二聚体模型。
6. 压力效应： 在 $T=0$ 时构建了广义焓（$H_{\text{tot}}$）框架，通过引入应力-功项来评估 $c$ 轴单轴压力的影响。

主要结果

• 结构能垒： DFT 计算显示，$\alpha$ 相与 $\beta$ 相之间存在约 300 meV/每两个位点单元格的显著能垒，这与实验观察到的热滞后现象一致。
• 热力学驱动力： 该转变是由内能与熵之间的竞争驱动的：
  • $\beta$ 相稳定性： 低温 $\beta$ 相通过层间二聚化带来的内能增益而稳定。这涉及强层间跳跃（$t^{(\text{inter})}_0 \approx -136.7$ meV），导致巨大的反铁磁交换作用（$J \approx 70$ meV）并形成自旋单态，从而猝灭了自旋熵。此外，$\beta$ 相中更强的层间耦合使声子模式硬化，降低了声子熵。
  • $\alpha$ 相稳定性： 高温 $\alpha$ 相由熵驱动稳定。它拥有较软的声子模式（零点能较低但熵较高），并由于较弱的层间耦合保留了顺磁自旋熵。
• 转变温度： 计算出的总亥姆霍兹自由能交点位于 $T_c \approx 113$ K，这与实验观测到的 $\sim$90 K 转变温度在半定量上是一致的。计算出的跨转变熵变与宏观比热测量结果相符。
• 压力响应： 在 $c$ 轴单轴压力下的广义焓计算表明，由于结构响应的差异，$\alpha$ 相在焓降低方面的压缩性不如 $\beta$ 相。预测该转变在约 2.6 GPa 的临界压力下被完全抑制，从而使 $\alpha$ 相在低温下保持稳定。

意义与主张
本文声称提供了一个关于 Nb$_3$Cl$_8$ 结构转变的统一热力学图景，将微观起源解析为声子熵和自旋熵的共同驱动，而非单一机制。通过显式分离晶格、自旋和关联贡献，研究表明 $c$ 轴单轴压力对转变的抑制是由于在广义焓景观中稳定了 $\alpha$ 相。这为实验观察到的压力抑制转变现象提供了热力学解释，并可能为在低温下获取 $\alpha$ 相的候选 QSL 态提供一条实验途径。作者指出，其结果是半定量的，承认了诸如声子的谐振近似以及压力下有限温度重整化等局限性。