How Accurately Can We Describe Spin Crossover?

技術要約：スピンクロスオーバーをいかに正確に記述できるか？

問題提起
金属錯体における熱的なスピンクロスオーバー（SCO）の転移温度（$T_{1/2}$）を予測することは、日常的な電子状態計算において依然として大きな課題である。核心的な困難は、分光学的な性質である分子のスピンクロスオーバーエネルギー（$\Delta E_{SCO}$）と、凝縮相の熱力学的性質である $T_{1/2}$ との間の複雑な関連性に存在する。$\Delta E_{SCO}$ は、転移を支配する極めて重要な分子量であるが、実験的に抽出することは困難である。その結果、報告されている実験値は稀であり、プロトコルに依存している。標準的な密度汎関数理論（DFT）のアプローチは、スピン対形成エネルギーと配位子場の分裂との間の微妙なバランスにより、計算されたエンタルピー（$\Delta H$）およびエントロピー（$\Delta S$）項に大きな不確実性が生じるため、正確な $T_{1/2}$ 値を予測できないことが多い。さらに、結合クラスター（CC）法のような高レベルの波動関数法は、特に従来のハートリー・フォック参照を用いる場合、これらの系に対して深刻な計算スケーリングの制限に直面する。

手法
著者らは、実験的な $T_{1/2}$ データから導出された $\Delta E_{SCO}$ のための参照データセットを生成するための、系統的な逆データ駆動型プロトコルを提案している。この手法は以下のステップで構成される：

1. データセットの構築： 熱駆動型のSCOを示すことが知られている24種類の金属錯体（Cr, Mn, Fe, Co）を収集した。
2. LC-$\omega$PBEによる逆エンジニアリング： 厳密な範囲分離パラメータ $\omega$ を持つことから、範囲分離ハイブリッド汎関数である LC-$\omega$PBE を選択した。各錯体について、実験的な $T_{1/2}$ を可能な限り正確に再現するように $\omega$ をチューニングした。これには、低スピン（LS）状態および高スピン（HS）状態の両方に対する構造緩和と調和振動解析が含まれる。
3. 参照データの抽出： $T_{1/2}$ に一致する最適な $\omega$ が決定された後、対応する分子スピンクロスオーバーエネルギー（$\Delta E_{SCO} = \Delta E_{HL} + \Delta E_{ZPE}$）を抽出した。このデータセットは、他の手法をベンチマークするための「参照」として機能する。
4. 感度分析： $\omega$ パラメータの微小な摂動（$\pm 1\%$）に対する $T_{1/2}$ および $\Delta E_{SCO}$ の感度を分析し、誤差範囲を確立した。
5. ベンチマーク： 作成した参照データセットを用いて、GGA、meta-GGA、グローバルハイブリッド、範囲分離ハイブリッド、およびダブルハイブリッドを含む様々な他の汎関数の性能を評価した。
6. 高レベル検証： r2SCAN 汎関数で最適化された構造を用い、8つの小さな錯体のサブセットに対して全電子結合クラスター計算（CCSD および CCSD(T)）を実施した。これらの計算では、参照状態の選択および励起の打ち切りによる影響を評価するために、非制限ハートリー・フォック（UHF）参照状態を利用した。

主な結果

• 参照データセットと不確実性： フィッティングされた $\omega$ パラメータは $0.177 \le \omega \le 0.528$ Bohr$^{-1}$ の範囲にあり、金属イオンの種類との単純な相関は見られなかったが、配位子のサイズが特定の傾向（例：マンガノセン類似体における傾向）に影響を与えていた。感度分析により、$\omega$ の $\pm 1\%$ の変動は、$T_{1/2}$ においては約 $\pm 50$ K、$\Delta E_{SCO}$ においては約 $\pm 2$ kJ mol$^{-1}$ の不確実性をもたらすことが明らかになった。特に、軟らかい調和振動モードを持つ錯体は、$\omega$ の摂動に対して高い感度を示した。
• DFT 汎関数の性能： 逆エンジニアリングされた参照データと比較した際、一般的な交換相関汎関数は顕著な偏差を示した。
  • GGA および Meta-GGA： PBE や TPSS などの標準的な汎関数は、約 57–59 kJ mol$^{-1}$ の中央偏差を示した。Non-hybrid の中では meta-GGA である r2SCAN が最も優れた性能を示し、中央偏差は 6.25 kJ mol$^{-1}$ であったが、それでも一部の種に対して非現実的な負の $\Delta E_{SCO}$ 値を生成した。
  • ハイブリッド： グローバルハイブリッド（PBE0, B3LYP）は、一般に非物理的な負のスピンクロスオーバーエネルギーを与え、偏差はしばく 40 kJ mol$^{-1}$ を超えた。範囲分離ハイブリッドおよびダブルハイブリッド（PWPB95, DSD-PBEB95）も必要な精度を達成できず、ダブルハイブリッドは最大の偏差（約 73 kJ mol$^{-1}$）を示した。
  • DFT に関する結論： テストされたどの汎関数も、$\pm 15$ kJ mol$^{-1}$（ターゲット範囲である $0-15$ kJ mol$^{-1}$ に対して約 100% の相対誤差）以内の精度を達成できなかった。著者らは、化学的に類似した集合内での傾向を研究するには、PBE+U$_{eff}$ および r2SCAN がコストとベネフィットの最良の妥協案であると述べている。
• 結合クラスターの性能： 8つの錯体のサブセットに対する CCSD(T) 計算は、UHF 参照を用いた場合に定性的な失敗を示すことを明らかにした。
  • ハートリー・フォックはスピンギャップを過小評価し（負の値）、
  • CCSD は改善は見られたものの、8つのうち6つの錯体に対して依然として負のスピンギャップを与えた。
  • CCSD(T) は、4つの錯体に対してのみ符号を修正した。残りの3つは、依然として大きな負の誤差（$>100$ kJ mol$^{-1}$）を保持していた。
  • 励起系列（CCSD から CCSD(T)）の非収束的な挙動は、UHF 参照がこれらの系に対して不適切であることを示唆しており、これは多参照性（multi-reference character）に起因すると考えられる。精度を得るには、より高次の励起や代替の参照状態が必要である。

意義と主張
本論文は、提案された逆エンジニアリング・プロトコルが、実験的な $T_{1/2}$ に制約された $\Delta E_{SCO}$ のための、系統的かつ再現可能な参照データセットを正常に生成したと主張している。このデータセットは、現在の電子状態手法の深刻な限界を浮き彫りにしている：

1. DFT の限界： 標準的、あるいは洗練された汎関数であっても、現実的な熱力学的記述に必要な予測精度を欠いている。関連するエネルギー尺度（$0-15$ kJ mol$^{-1}$）は室温における熱エネルギーと同程度であり、小さな誤差が $T_{1/2}$ 予測において壊滅的な結果をもたらす。
2. CC の限界： 高レベルの波動関数法は、理論的には優れているものの、実用的な障壁に直面している。UHF 参照への依存と励起系列の打ち切りは、非系統的な収束と定性的な失敗を引き起こしている。
3. ベンチマークの必要性： 著者らは、遷移金属に関するベンチマークデータセットの不足を強調し、将来の汎関数開発には、これらのバイアスに対処するために現実的な遷移金属テストセットを組み込む必要があると主張している。彼らは、逆エンジニアリングされた $\Delta E_{SCO}$ 値が、経験的なチューニングなしでは正確な予測を提供できない現在の標準的な DFT の能力を鑑みると、高レベルの波動関数法の信頼性を評価するための必要なベンチマークとして機能すると考えている。

著者らは、自身のプロトコルが参照エネルギーを抽出するものである一方で、固有の不確実性（$\pm 2$ kJ mol$^{-1}$）や、DFT および CC の両方がこれらの値を一貫して再現できないという事実は、正確で日常的な予測が依然として計算化学における未解決の課題であることを示している、として控えめな立場をとっている。