[論文レビュー] Ab initio Simulations of Fe-based Ferric Wheels

A. V. Postnikov, Jens Kortus+1 more|arXiv (Cornell University)|Jul 11, 2003

Magnetic Properties and Applications参考文献 17被引用数 4

ひとこと要約

本研究では、ab initio密度汎関数理論（DFT）を用いて、六核性Feを含むフェリックホイル（MFe6[NC2H4O]6Cl, M = Li, Na）の電子的および磁気的構造を調査した。その結果、磁気モーメントがFeにのみ局在しているのではなく、Fe、O、N原子にわたって非局在化していることが明らかになった。主な発見は、Feにおける局在的モーメントが約4 μBであり、これはFe(II)の性質を示しており、周囲の原子からの寄与を含めることで、1FeサイトあたりS = 5/2のスピンが正確に再現されることである。特に、上面の酸素原子が最大のスピン偏極を示している。

ABSTRACT

Based on first-principles density-functional theory calculations we investigate the electronic structure of hexanuclear "ferric wheels" M Fe_6[N(CH_2 CH_2 O)_3]_6 Cl (M = Li, Na) in their antiferromagnetic ground state. The electronic structure is presented in form of spin- and site-resolved local densities of states. The latter clearly indicate that the magnetic moment is distributed over several sites. The local moment at the iron site is still the largest one with about 4 mu_B, thus indicating the valence state of iron to be closer to Fe(II) than to commonly accepted Fe(III). The local spin of S=5/2 per iron site, following from magnetization measurements, is perfectly reproduced if one takes the moments on the neighbor atoms into account. The largest magnetic polarization is found on the apical oxygen atom, followed by nitrogen bridging oxygens. These findings are confirmed by a map of spatial spin density. A further goal of the present study has been a comparative test of two different DFT implementations, Siesta and NRLMOL. They yield a very good agreement down to small details in the electronic structure.

研究の動機と目的

• 第一原理的手法を用いて、六核性Feを含むフェリックホイル（MFe6[NC2H4O]6Cl, M = Li, Na）の電子的および磁気的構造を理解すること。
• 観測された1FeサイトあたりS = 5/2と予想されるFe(III)状態との乖離を解消するため、系全体におけるモーメント分布を分析すること。
• 開殻遷移金属を含む複雑な分子磁性体に対して、2つのDFT実装（SiestaとNRLMOL）の信頼性と一貫性を検証すること。
• 電子構造と実験的磁化データを結びつけることで、これらの単分子磁性体の磁気的挙動の微視的基盤を提供すること。

提案手法

• 電子相関効果を考慮するため、LDA+U法を用いた一般化勾配近似（GGA）に基づく第一原理DFT計算。
• 電子構造および磁気モーメント分布の分析のため、スピンおよびサイトに分解された局所状態密度（LDOS）の計算。
• NRLMOLコードを用いて空間的スピン密度マップの構築により、磁気偏極の非局在化を可視化。
• 結果の妥当性を確認するため、独立した2つのDFTコード（SiestaとNRLMOL）を用いた相互検証により、電子構造予測の堅牢性を確保。
• 交換結合定数（J）の推定のため、強磁性（FM）および反強磁性（AFM）状態の全エネルギーの計算。
• 構造的緩和およびスピン密度マップの分析により、配位子原子（O、N）が磁気モーメント分布に果たす役割を評価。

実験結果

リサーチクエスチョン

• RQ1 フェリックホイル内のFe、O、N原子にわたる磁気モーメントの分布はどのようになっており、局在的Fe(III)状態からどれほど逸脱しているか？
• RQ2 周囲の原子（O、N）が総スピンにどれほど寄与しているか。また、これらの寄与を含めることで、観測された1FeサイトあたりS = 5/2が再現可能か？
• RQ3 本分類の分子磁性体に対して、2つの異なるDFT実装（SiestaとNRLMOL）の結果の正確性と一貫性はどの程度か？
• RQ4 配位子原子は、磁気的相互作用を媒介し、スピン偏極にどのように寄与しているか？

主な発見

• Feサイトにおける局在的磁気モーメントは約4 μBであり、一般的に想定されるFe(III)状態とは異なり、Fe(II)に近い酸化状態を示している。
• 1FeサイトあたりS = 5/2のスピンが、モーメント解析に周囲の原子の寄与を含めることで定量的に再現可能である。
• スピン密度マップにより、上面の酸素原子が最大の磁気偏極を示しており、次に窒素原子、ブリッジ酸素原子が続く。
• スピン密度マップにより、酸素および窒素原子に顕著なスピン偏極が確認されたが、炭素や水素サイトには検出可能な磁化は認められなかった。
• SiestaとNRLMOLの2つのDFTコードは、局所状態密度および空間的スピン密度分布において優れた定量的一致を示し、窒素原子付近を除いてわずかな差異しか認めなかった。
• 全エネルギー差（FM–AFM）から推定された交換結合定数Jは、実験データと整合的であり、J ≈ -38 Kの値を示し、以前の推定値よりも実験値に近い結果となった。