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領域を抑制してしまう。この転移の微視的な起源については、電荷不均一化や層間ダイマー・シングレット形成といった競合する説明が存在し、議論が続いている。さらに、熱力学的な駆動力や、外部刺激（一軸圧力など）がこの転移を抑制するメカニズムについても、完全には理解されていない。

手法
著者らは、格子、電子、およびスピンの寄与を分離するために、第一原理密度汎関数理論（DFT）と拡張ハバードモデルを組み合わせた統一的な自由エネルギー・フレームワークを用いている。

1. DFT計算: 非スピン分極DFT計算（VASPを用い、optB88-vdWおよびPBE汎関数を使用）を行い、$\alpha$相（空間群 $P\bar{3}m1$）および$\beta$相（$R\bar{3}m$で表される）の両方について、緩和された構造、基底状態エネルギー、およびフォノンスペクトルを得た。両構造間のエネルギー障壁を推定するために、線形補間パスを構築した。
2. フォノン熱力学: 調和近似の範囲内で、密度汎関数摂動論（DFPT）を用いてフォノン自由エネルギー（$F_{\text{phonon}}$）およびエントロピーを計算した。
3. 有効モデル: 半充填の2$a_1$多様体におけるモット物理学に対処するため、クラスター中心軌道のワニエ基底を用いて拡張ハバードモデルを構築した。パラメータ（ホッピング $t$、オンサイト $U$、およびサイト間 $V$ クーロン相互作用）は、制限されたランダム位相近似（cRPA）を用いて導出した。
4. 二重計上補正: DFT参照系ですでに捉えられている低エネルギーの2$a_1$寄与を差し引くために、ハートリーレベルの二重計上補正（$E_{\text{DC}}$）を実装し、一貫した全自由エネルギー（$F_{\text{tot}}$）を確保した。
5. スピン熱力学: 拡張ハバードモデルを有効ハイゼンベルクモデルへと簡約した。$\alpha$相については、連結クラスター高温展開（HTSE）を用いてスピンエントロピーを計算した。$\beta$相については、強い層間結合があるため、簡略化された孤立ダイマーモデルを採用した。
6. 圧力効果: $c$軸一軸圧力下でのエンタルピー減少を評価するため、$T=0$ における一般化エンタルピー（$H_{\text{tot}}$）のフレームワークを構築し、応力-仕事項を組み込んだ。

主な結果

• 構造障壁: DFT計算により、$\alpha$相と$\beta$相の間には、2サイト単位胞あたり約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領域へアクセスするための実験的な経路を提供する可能性がある。著者らは、フォノンの調和近似や、圧力下での非調和効果および有限温度での繰り込みの無視といった限界を認め、結果が半定量的であることを注記している。