来源论文: https://arxiv.org/abs/2607.00599v1 生成时间: Jul 02, 2026 15:56

熵驱动的团簇莫特绝缘体 Nb3Cl8 结构相变：结合第一性原理、多体模型与声子热力学的全景深度解析

0. 执行摘要

过渡金属硫族/卤族化合物中的强关联电子器件和新型量子物态（如量子自旋液体 Quantum Spin Liquid, QSL）是当前凝聚态物理与量子化学领域的前沿热点。作为过渡金属卤化物中独具特色的代表，二维过渡金属过渡团簇材料 $\text{Nb}_3\text{Cl}_8$ 展现出了极其典型的“团簇莫特绝缘体”（Cluster Mott Insulator）物理特性。其高温 $\alpha$ 相具有二维三角晶格结构，每个分子团簇拥有一个未配对的 $d$ 电子，是探索量子自旋液体物态的极佳候选平台。然而，实验表明，$\text{Nb}_3\text{Cl}_8$ 在约 90 K 温度下会经历一个显著的一阶结构相变，转变为低温非磁性的 $\beta$ 相，这极大地抑制了 $\alpha$ 相在低温极限下展现量子自旋液体物态的可能性。对该相变的微观物理起源及其热力学驱动力的解析，一直是学界争论的焦点。

本研究提出了一种创新性的、无缝耦合第一性原理密度泛函理论（DFT）与多体模型哈密顿量（Extended Hubbard & Heisenberg Model）的自由能计算框架，对 $\text{Nb}_3\text{Cl}_8$ 系统内的晶格、声子、自旋和关联电子等多自由度耦合进行了精细化剥离与量化分析。计算结果表明，该相变是一个典型的“熵驱动相变”。高温 $\alpha$ 相通过更加柔顺的声子模式（低频声子软化）以及保留的顺磁自旋自由度提供了极大的声子熵与自旋熵贡献，从而在高温下稳定存在；而低温 $\beta$ 相则通过层间二聚化（Interlayer Dimerization）产生强烈的层间化学键合，使体系获得了极大的电子轨道结合能（内能项），但在同时淬灭了磁矩并硬化了声子谱。此外，通过构建零温下的广义焓变（Generalized Enthalpy）曲线，研究进一步成功解释了 $c$ 轴单轴压力如何通过选择性压制 $\beta$ 相的压缩空间来稳定 $\alpha$ 相，为在实验上通过应变调控手段稳定 $\alpha$ 相、进而观测低温量子自旋液体态提供了清晰、定量的理论指引。

1. 核心科学问题，理论基础，技术难点，方法细节

1.1 核心科学问题：$\text{Nb}_3\text{Cl}_8$ 相变的争论与本质

$\text{Nb}_3\text{Cl}_8$ 是典型的范德华力层状材料，其核心结构单元是由三个 $\text{Nb}$ 原子和十三个 $\text{Cl}$ 原子形成的 $\text{Nb}_3\text{Cl}_{13}$ 团簇。在该团簇内部，由于强烈的 $\text{Nb}-\text{Nb}$ 金属键合，$\text{Nb}$ 原子发生三聚化（Trimerization）。这导致每个团簇在分子轨道图像下仅剩下一个未配对电子占据在 $2a_1$ 轨道上，展现出有效自旋 $S = 1/2$。在高温 $\alpha$ 相中，这些团簇在二维平面内排列成高度阻挫的三角晶格，层间耦合较弱，形成一个理想的阻挫强关联莫特物理环境。

当温度降低到约 90 K 时，材料发生一阶相变，磁化率骤降。传统观点对此有以下几种竞争模型：

1. 电荷不成比例分配（Charge Disproportionation）模型：认为相变导致了 $2[\text{Nb}_3]^{8+} \rightarrow [\text{Nb}_3]^{7+} + [\text{Nb}_3]^{9+}$ 的电荷有序转变，磁矩淬灭来源于无配对电子的 $[\text{Nb}_3]^{7+}$ 与 $[\text{Nb}_3]^{9+}$ 的交替排列。然而，不同对称性破缺模型中关于化学键长变化的计算与实验偏差较大，此机制一直缺乏足够的能量学支持。
2. 层间二聚化单态（Interlayer Dimer Singlet）模型：认为相变并未导致电荷分配不均，而是伴随着相邻层间团簇的滑移。在低温 $\beta$ 相中，层间滑移使得某些层间对转变为“重叠（Eclipsed）”排布，从而诱发了极强的层间电子杂化（层间二聚化），形成非磁性价键单态（Valence Bond Singlet），从而淬灭了局域自旋。

本工作的核心科学问题在于：层间二聚化模型是否能够提供足够的热力学驱动力来解释该一阶结构相变？在相变点两端，电子内能、声子振动能、声子熵、自旋熵各自扮演了何种角色？

1.2 理论基础与多自由度剥离

对于典型的强关联体系，标准的单粒子近似密度泛函理论（DFT）往往失效。由于强库仑排斥作用，$\text{Nb}_3\text{Cl}_8$ 的 $2a_1$ 轨道在费米能级附近发生了莫特劈裂。因此，不能直接使用非自旋极化的 DFT 总能量作为判定相变能量基准的唯一标准，更无法利用其获取真实的自旋热力学贡献。本研究提出了如下的 Helmholtz 自由能分解框架：

其中，各项的物理意义及构建方式如下：

• $E_{\text{DFT}}$：由非自旋极化第一性原理计算得到的晶格静态基态能量。它提供了背景晶格格架和全电子（排除了 $2a_1$ 活性关联空间之后）的背景键合内能。
• $F_{\text{phonon}}(T)$：基于谐振近似下的声子自由能。通过计算全声子谱，积分得到晶格零点能（ZPE）及有限温度下的声子热激发和熵效应贡献： $$F_{\text{phonon}}(T) = \sum_{\mathbf{q}\nu} \left[ \frac{1}{2}\hbar\omega_{\mathbf{q}\nu} + \frac{1}{\beta} \ln\left( 1 - e^{-\beta\hbar\omega_{\mathbf{q}\nu}} \right) \right]$$
• $E_{\text{DC}}$（双重计数纠正，Double-Counting Correction）：由于我们在后续的多体模型哈密顿量中显式地计算了 $2a_1$ 活性空间中的电子动能与库仑相互作用，必须系统性地将标准 DFT 中已经近似包含的这部分贡献扣除，以避免双重计数。其形式包含单体动能项纠正与 Hartree 平均场库仑能纠正： $$E_{\text{DC}} = E_{\text{hop}}^{\text{DC}} + E_{U}^{\text{DC}}$$ $$E_{\text{hop}}^{\text{DC}} = \frac{2}{N_k} \sum_{\mathbf{k},\nu}^{\text{occ}} \langle \psi_{\mathbf{k}\nu} | \hat{h}_W(\mathbf{k}) - \hat{h}_{\text{on}} | \psi_{\mathbf{k}\nu} \rangle$$ $$E_{U}^{\text{DC}} = \sum_{i} U_i \langle \hat{n}_{i\uparrow} \rangle \langle \hat{n}_{i\downarrow} \rangle = \sum_i \frac{U_i}{4}$$
• $F_{\text{Hub}}(T)$ 或 $F_{\text{Hei}}(T)$：由活性关联空间（$2a_1$ 轨道）构建的有效 Hubbard 模型或进一步简化的 Heisenberg 自旋模型所提供的有限温度自由能项，用以高精度捕捉多体自旋关联效应和顺磁自旋熵。

1.3 技术难点与解决方法

1. 强关联体系中局域轨道的精确抽取：利用极宽能区投影往往由于混杂导致物理图像模糊。本工作使用 Wannier90 软件包，基于 $2a_1$ 分子轨道特征构建了最大局域化 Wannier 函数（MLWFs），从而在能量窗口内获得一个精准单能带的紧束缚模型（Tight-Binding Model）。
2. 库仑排斥参数 $U$ 和 $V$ 的定量计算：不能使用经验值。本工作基于受约束随机相位近似（constrained Random Phase Approximation, cRPA）高精度评估了层内及层间的动态屏蔽库仑参数，消除了模型中的拟合参数依赖。
3. 强关联自旋阻挫体系的自旋自由能计算：在过渡温度 $T_c \sim 100$ K 下，量子自旋效应依然显著。对于二维强阻挫三角晶格自旋系统，蒙特卡洛面临严重的负符号问题。研究团队采用了**高温级数展开（High-Temperature Series Expansion, HTSE）**方法，将自旋哈密顿量的配分函数展开至第 10 阶，并结合巴德逼近（Padé Approximants）进行解析延拓，获得了 $\alpha$ 相极高精度的顺磁自旋熵和自由能曲线。对于低温 $\beta$ 相，由于其极强的层间反铁磁交换作用（二聚化），自旋自由能采用基于隔离二聚体模型（Isolated Dimer Model）的微扰理论进行解析求解。

2. 关键 Benchmark 体系、计算所得数据与性能数据

2.1 结构模型与能量学分析

研究团队首先对多种候选的低温 $\beta$ 相结构（如 $C2/m$、$R3$、$R\bar{3}m$）进行了精细的 DFT 结构弛豫与基态能量 Benchmark 分析：

• 高温 $\alpha$ 相空间群：$P\bar{3}m1$（No. 164），其层间滑移展现出高度对称的均匀交错 Stacking ($AA'$) 结构。
• 实验提出的低温 $\beta$ 相三种候选模型在 DFT 弛豫后的能量差极小：
  • 相较于高对称性的 $R\bar{3}m$（No. 166），低对称性的 $R3$ 结构在能量上仅低了约 $4\text{ meV/f.u.}$；
  • $C2/m$ 结构在能量上低了约 $2\text{ meV/f.u.}$。
  • 鉴于这几类结构内部的层间排列（Stacking sequence）完全一致（均为 $AA'BB'CC'$ 的二聚化特征），且局域畸变产生的能量差极微，研究团队选取极具代表性且易于计算的 $R\bar{3}m$ 结构作为 $\beta$ 相的标准代表体系开展后续物理建模。

层间滑移（Slip）是该一阶相变的控制参数。通过在线性插值路径上移动滑移参量，DFT 能量表现出两相之间高达 $\sim 300\text{ meV}$（每两个 Nb 团簇单元）的势垒，这完美解释了实验中观测到的剧烈一阶相变特性和显著的热滞回线行为。

2.2 紧束缚模型与库仑参数的 Benchmark 结果

利用 cRPA 结合 MLWFs 在 DFT 计算基础上构建的有效 Extended Hubbard 模型的关键参数如下表所示（单位均为 $\text{meV}$）：

数据物理透析：

• 层间跳符的骤变：在 $\alpha$ 相中，由于 Stacking 的几何位置，层间跳符 $t_0^{\text{inter}}$ 仅为 $-17.3\text{ meV}$；而转变至 $\beta$ 相后，部分层间键发生“重叠”（Eclipsed），这导致层间跳符暴增至 $-136.7\text{ meV}$（提升了近 8 倍）。这一跳符骤增是层间二聚化单态形成的最直接物理驱动力。
• 有效自旋 Heisenberg 模型交换常数 $J$ 的映射：根据二阶超交换作用（Superexchange）公式： $$J = \frac{4t^2}{U-V}$$ 对上述参数进行映射，得到体系层内和层间自旋反铁磁交换参数（单位为 $\text{meV}$）：

• 在 $\alpha$ 相中，层内和层间交换作用极其微弱（均为 $1-2\text{ meV}$），说明系统自旋在温度高于数个开尔文时表现出极强的顺磁性，拥有大量的自旋自由度。其自旋熵极限高达 $S_{\text{spin}} = k_B \ln 4$（每双自旋单元）。
• 在 $\beta$ 相中，层间强二聚键处的 $J_0^{\text{inter}}$ 跃升至高达 $70.34\text{ meV}$ （折合约 $815\text{ K}$）。在 90 K 转变点附近，该交换作用远超热涨落能级（$k_B T \sim 8\text{ meV}$），从而将相邻层间的自旋强力锁死在单态上，淬灭了几乎所有的顺磁自旋熵。

2.3 声子动力学与有限温度热力学 Benchmark

• 全谱声子谱计算：计算表明，$\alpha$ 相中层间剪切和层间呼吸等低频剪切模式（如 $\Gamma$ 点的层间呼吸模，约 $1.729\text{ THz}$）相较于 $\beta$ 相（约 $1.826\text{ THz}$）发生显著的“红移”（软化）。
• 声子热力学平衡：由于 $\alpha$ 相整体声子频带更宽且在低频区更为柔顺（软化），其振动自由能曲线随温度上升下降得更为迅速。在 $T = 0\text{ K}$ 时，$\beta$ 相由于强层间二聚键形成的更低基态能占优；但随着温度提升，$\alpha$ 相的声子熵和自旋熵双重累积，最终在理论预测的 $T_c \approx 113\text{ K}$ 处与 $\beta$ 相自由能线发生交叉，与实验所得的 $T_c \sim 100\text{ K}$ 高度契合。

3. 代码实现细节、复现指南及开源软件包指南

为方便相关领域的量子化学和计算凝聚态物理科研人员复现该工作，以下整理了完整的计算流程图景与核心数据提取方法：

+-----------------------+      +---------------------------+
| 1. VASP 晶格弛豫 (DFT) | ---> | 2. VASP + Phonopy (DFPT)  |
+-----------------------+      +---------------------------+
            |                                |             
            v                                v             
+-----------------------+      +---------------------------+
| 3. Wannier90 (MLWFs)  |      |  获得声子 Helmholtz 自由能  |
+-----------------------+      +---------------------------+
            |                                              
            v                                              
+-----------------------+                                  
| 4. VASP (cRPA)        |                                  
+-----------------------+                                  
            |                                              
            v                                              
+-----------------------+                                  
| 5. 模型哈密顿量解析建模| ---> 联合构筑 F_tot (T) 自由能曲线
+-----------------------+                                  

3.1 核心复现步骤与设置详解

步骤 1：VASP 静态电荷密度与晶格弛豫

• 使用非自旋极化设置开展计算（ISPIN = 1）。
• 泛函选择：必须考虑范德华修正以获得准确的层状滑移能。推荐采用 optB88-vdW 泛函（设置 GGA = BO 和 LUSE_VDW = .TRUE.）或标准的 PBE 泛函进行对比 Benchmark。
• 参数推荐：

  ENCUT = 450
  EDIFF = 1E-8
  EDIFFG = -1E-4
  ISMEAR = 0; SIGMA = 0.05
  KPOINTS = 4 x 4 x 4 (对于单胞)
  

步骤 2：利用 Wannier90 构建 $2a_1$ 轨道紧束缚模型

• 在 VASP 运行完成后，开启 LWANNIER90 = .TRUE.。
• wannier90.win 中主要投影参数设置：

  begin projections
  Nb:d_z2  # 或根据团簇对称性构筑本地化的 2a1 分子轨道线性组合
  end projections
  dis_win_min = -1.5
  dis_win_max = 1.0
  

  通过拟合获得精确对应 $2a_1$ 单带的紧束缚哈密顿量实空间矩阵元素 $t_{ij}$。

步骤 3：利用 cRPA 抽取关联相互作用 $U$ 和 $V$

• 在 VASP 中使用 GW 近似下的 cRPA 计算。需要包含足够的能带数（论文中设置 NBANDS = 288）以保证极佳的空带收敛性。
• 设置 ALGO = CRPA，并限制投影体系为上面抽取的 $2a_1$ 活性轨道空间（Wannier 轨道）。

步骤 4：声子动力学复现（VASP + Phonopy）

• 使用密度泛函微扰理论（DFPT）方法（在 VASP 中设置 IBRION = 8）。
• 构筑 $2\times 2\times 2$ 超胞以完全克服动力学力常数实空间的截断误差。
• 运行完毕后，利用 Phonopy 提取实空间力常数矩阵，并进行声子态密度的细密网格积分采样（设置密集 $q$ 网格为 $31\times 31\times 31$）：

  phonopy -p --factor=524.497 -c POSCAR-unitcell mesh.conf
  

  通过 mesh.conf 导出热力学函数 thermal_properties.yaml，提取有限温度下的声子 Helmholtz 自由能 $F_{\text{phonon}}(T)$、内能及声子熵。

3.2 自旋哈密顿量的 HTSE 解析求解与 Padé 延拓算法

研究团队利用开源 HTSE-code 软件包求解 $S = 1/2$ Heisenberg 模型。以下展示一段用于实现 HTSE 自由能插值与 Padé 逼近的外推逻辑伪代码，供量子化学工作者构建自由能脚本时参考：

import numpy as np
from scipy.interpolate import interp1d

def pade_approximant_4_4(beta, coeffs):
    """
    计算 P[4, 4] 形式的巴德逼近
    f(\beta) \approx (p0 + p1*\beta + p2*\beta^2 + p3*\beta^3 + p4*\beta^4) / 
                     (1 + q1*\beta + q2*\beta^2 + q3*\beta^3 + q4*\beta^4)
    coeffs 为通过 HTSE 代码算出的各阶泰勒系数 [a0, a1, a2, a3, a4, a5, a6, a7, a8]
    """
    # 构造线性方程组求解 P[4, 4] 的分子分母多项式系数
    # 此处省略具体求解线性系统 Ap = B 的过程
    # p_coeffs, q_coeffs = solve_pade_system(coeffs, L=4, M=4)
    
    # 示例返回
    num = sum(p_coeffs[i] * (beta**i) for i in range(5))
    denom = 1.0 + sum(q_coeffs[j] * (beta**(j+1)) for j in range(4))
    return num / denom

def calculate_spin_entropy_alpha(T, J_intra, J_inter):
    """
    输入交换常数，基于级数项外推得到阿尔法相的高温顺磁自旋熵
    """
    # KB 为玻尔兹曼常数
    KB = 8.617333262145E-5 # eV/K
    beta = 1.0 / (KB * T)
    
    # 这里的泰勒展开系数由 HTSE 计算得出
    htse_coeffs = get_htse_coefficients_for_triangular(J_intra, J_inter)
    
    # 使用 P[4, 6] 或 P[4, 4] 进行自旋自由能外推
    f_spin = pade_approximant_4_4(beta, htse_coeffs)
    
    # 熵为自由能对温度的负导数
    # S = -dF/dT
    return f_spin

3.3 开源工具链及 Repo 推荐：

1. VASP（第一性原理平台，商业收费）：https://www.vasp.at/
2. Wannier90（MLWFs 投影）：https://github.com/wannier-developers/wannier90
3. Phonopy（晶格振动热力学）：https://github.com/phonopy/phonopy
4. HTSE-code（自旋高温展开）：https://github.com/scipost-phys/htse-code

4. 关键引用文献与局限性批判评论

4.1 关键引用文献

本研究所立足的关键基础科学文献包含以下领域：

• 团簇莫特物理奠基与 Nb3Cl8 体系发现：
  • Imada, A. et al., Rev. Mod. Phys. 70, 1039 (1998). (强关联体系综述)
  • Sheckelton, J. P. et al., Inorg. Chem. Front. 4, 481 (2017). (报道 $\text{Nb}_3\text{Cl}_8$ 的一阶相变与自旋单态形成的实验先驱工作)
• HREELS 谱学实验证据：
  • Su, Z. et al., Phys. Rev. Lett. 136, 1065 (2026). (利用高分辨率电子能量损失谱证实 $\alpha-\beta$ 相变伴随着层间激子能带自二维线性色散向三维抛物线色散的重构，强力支持层间二聚化物理图像)
• 模型哈密顿量与计算物理方法：
  • Aryasetiawan, F. et al., Phys. Rev. B 70, 195104 (2004). (cRPA 理论基础)
  • Pierre, L., Bernu, B., and Messio, L., SciPost Phys. 17, 105 (2024). (开发了计算各向异性自旋 1/2 三角晶格的 HTSE 高级算法工具)

4.2 本工作局限性批判

尽管本工作成功预测了相变温度及压力对相变的压制效应，但在严苛的量子化学与关联电子学标准下，本工作仍存在以下局限性：

1. 谐振声子近似（Harmonic Phonon Approximation）的局限： 研究在计算 $F_{\text{phonon}}(T)$ 时采用了静态零温结构下的谐振声子谱。实际上，$\text{Nb}_3\text{Cl}_8$ 作为范德华层状材料，其层间滑移低频振动具有极强的非谐振（Anharmonic）效应。在相变转变点（$\sim 100\text{ K}$）附近，温度导致的声子重整化（Phonon Renormalization）可能非常显著，这将改变低频剪切声子模随温度变化的斜率，对相变温度的定量预测带来修正。
2. 静态库仑屏蔽近似： cRPA 计算出的 $U$ 和 $V$ 为静态极限（$\omega \rightarrow 0$）。然而，对于低能有效模型，高能区的轨道杂化会引入频率依赖的屏蔽相互作用 $U(\omega)$。直接忽略频变库仑作用将导致在计算自旋超交换参数 $J$ 时存在约 $10\% - 15\%$ 的系统偏差。
3. 非非自旋极化背景下的电子-声子耦合遗失： 在声子谱计算中，DFT 参考态为非自旋极化的非磁性状态，未能将 $2a_1$ 轨道上真实的莫特局域自旋和莫特能隙对声子晶格力常数的影响完全自洽地考虑进去（即缺乏自洽的 Spin-Phonon Coupling 修正）。
4. $\beta$ 相低温极限的极简化处理： 计算将低温 $\beta$ 相的低能激发简化为隔离二聚体模型。尽管 $\beta$ 相中 $J_0^{\text{inter}}$ 极强（$\sim 70\text{ meV}$），但层内剩余的微弱反铁磁相互作用（$\sim 1.58\text{ meV}$）在极低温极限下（$T < 10\text{ K}$）依然会引发精细的磁性相行为或复杂的集体自旋轨道涨落，隔离二聚体模型在此能区将不再适用。

5. 其他必要的学术补充与推导

5.1 二位三角晶格上超交换作用 $J$ 的微观解析推导

为了让量子化学工作者更直观地理解如何从单带 Hubbard 模型跃迁至 Heisenberg 磁学模型，以下给出由 Supplementary Material 所阐述的、基于两中心体系的二阶微扰详细数学推导过程：

考虑单带 Extended Hubbard 模型在半填充状态下的双中心体系（点 $1$ 和点 $2$），其低能单占基态空间由 4 个状态构成：

高能双占受激态空间为：

其完整的系统哈密顿量 $H$ 在此基底下的矩阵表达形式为：

在强关联极限下（$U - V \gg t$），利用 Löwdin 简并微扰理论（Degenerate Perturbation Theory），将高能受激态空间投影回低能有效自旋空间，可求得低能有效哈密顿量矩阵元：

对于对角项，如 $\langle 1\uparrow, 2\downarrow| H_{\text{eff}} |1\uparrow, 2\downarrow\rangle$，两个中间虚跃迁过程各自贡献 $\frac{t^2}{-(U-V)}$，累计为 $\frac{2t^2}{-(U-V)}$。对于非对角项，由于电子自旋交换跃迁，累计贡献为 $\frac{2t^2}{U-V}$。

将此有效矩阵与 Heisenberg 自旋耦合算符 $\hat{H}_{\text{Heisenberg}} = J \left( \mathbf{S}_1 \cdot \mathbf{S}_2 - \frac{1}{4} \right)$ 在对应基底下的矩阵表达进行一一映射配对：

通过对角元或非对角元矩阵系数的精准匹配，可完美导出超交换作用常数与模型微观参数的关系：

这一经典物理公式不仅被成功应用于评估层内和层间的磁性强度，更深刻揭示了为什么层间滑移导致的层间跳符 $t_0^{\text{inter}}$ 暴增，会以二次方的速度级联放大层间磁性自旋耦合，直接拉升了 $\beta$ 相中自旋单态锁死的发生温度。

5.2 团簇莫特绝缘体（Cluster Mott）与常规莫特绝缘体（Conventional Mott）的区别

在量子化学的视野中，通常的莫特绝缘体（如典型的过渡金属单氧化物 $\text{NiO}$ 或过渡金属二硫化物 $\text{TaS}_2$ 的某些异质相）其未配对的电荷定位在单一过渡金属原子的局域 $d$ 轨道上。

而对于以 $\text{Nb}_3\text{Cl}_8$ 为代表的团簇莫特绝缘体：

1. 局域化中心是分子/团簇：电荷局域化的物理实体并非单个 $\text{Nb}$ 阳离子，而是由强共价键三聚化形成的整个 $[\text{Nb}_3]^{8+}$ 团簇核。其未配对电子填充在一个完全离域于这三个金属核顶点的分子轨道（大自旋中心，$\text{molecular Wannier orbital}$）上。
2. 对晶格畸变的敏感性显著放大：由于活性 Wannier 轨道是由分子级分子轨道线性组合而来，其层间电子云的形状、对称性与层内和层间的物理滑移和畸变高度敏感。这也是为什么通常的莫特绝缘体相变常伴随轨道有序，而 $\text{Nb}_3\text{Cl}_8$ 却展现出剧烈层间滑移（二聚化）相变的原因。

5.3 单轴压力调控机制的深度剖析

通过计算 $T = 0$ 时的广义焓变（Generalized Enthalpy）：

结果发现，$\alpha$ 相与 $\beta$ 相展现出了截然不同的可压缩性（Compressibility）：

• 在 $\beta$ 相中，层间由于已经形成了强二聚化的共价价键单态，其沿着 $c$ 轴的原子排布高度紧凑，层间距离极短，因此其对抗外加单轴压力的可压缩空间极小，晶格极其坚硬。
• 在 $\alpha$ 相中，层间仅由相对较弱的范德华作用和轻微的静电引力维持，晶格在层间方向上高度柔顺、易被压缩。

当施加 $c$ 轴单轴压力时，晶格在压应力下收缩。由于 $\alpha$ 相更容易被压缩，其力学响应能（Stress-work Term, $-P_z A^{(xy)}_0(z_0 - z)$）下降得远比 $\beta$ 相要迅速。根据计算结果，当单轴应力超过约 $2.60\text{ GPa}$ 的临界阈值后，$\alpha$ 相的广义焓变将始终低于 $\beta$ 相。这就意味着在单轴压力调控下，一阶结构相变将被彻底消除，$\alpha$ 相可以稳定保持到极低温极限。这一结论极大地鼓舞了实验物理学家通过基底应变或高压实验，将优秀的量子自旋液体候选相稳定至低至 $\text{mK}$ 的温区，从而在该体系中寻找分数化激发的决定性物理证据。

本研究所构建的这一套“第一性原理+关联哈密顿量+级数展开自旋热力学”的复合物理分析范式，未来在计算量子化学和固体关联物性理论领域，必将发挥越来越关键的基准作用。