Synchrotron x-ray diffraction and DFT study of non-centrosymmetric EuRhGe3 under high pressure

本研究は、放射光X線回折と密度汎関数理論（DFT）計算を組み合わせることで、非中心対称性を持つEuRhGe3の高圧下における構造挙動を調査し、相転移のない35 GPaまでの滑らかな体積収縮、異方的な格子圧縮、および非整数価のEuに起因する高圧下での実験値と理論値の体積の乖離を明らかにしている。

要約

ユロピウム、ロジウム、ゲルマニウムからなる、極めて精巧で微細な結晶を想像してみてください。この結晶を、微視的な三次元の足場、あるいはレゴブロックのような構造物として考えてみてください。この論文の科学者たちは、この構造を、まるで巨大でハイテクな万力に入れた時のように、猛烈に強く押しつぶしたらどうなるのかを知りたかったのです。

以下に、彼らの実験の物語を分かりやすく解説します。

セットアップ：高額な賭けの「絞り込み」

研究者たちは、EuRhGe3 という結晶を取り上げました。これは単なる結晶ではありません。特殊な「非対称な（中心対称性のない）」形状を持っており、それが興味深い磁気特性を生み出しています。

これをテストするために、彼らは普通の万力は使いませんでした。代わりにダイヤモンドアンビルセルを使用しました。二つのダイヤモンド（地球上で最も硬い物質）が互いに押し付け合っている様子を想像してください。結晶はその間で押しつぶされ、周囲のヘリウムガスによって、高圧の潜水艦のように圧力が均一に保たれます。彼らは、海面での気圧の35,000倍という圧力に達するまで、この結晶を押しつぶしました。

主な発見：スムーズな圧縮、破断なし

通常、物を強く押しつぶしすぎると、パキッと折れたり、壊れたり、あるいは突然形が変わったり（相転移）します。例えば、スポンジが突然岩に変わってしまうようなものです。

しかし、この結晶は驚くほど弾力性がありました。

• 壊れない： あの凄まじい圧力の下でも、結晶は壊れることも、根本的な形を変えることもありませんでした。限界値に達するまで、元の「レゴのパターン」を維持し続けたのです。
• 小さくなる： 壊れる代わりに、ストレスボールを握りつぶす時のように、ただどんどん小さくなっていきました。構造全体がスムーズに収縮したのです。

捻り：片側がより早く縮む

ここからが面白いところです。この結晶は完璧な立方体ではなく、少し背の高い、細長い箱のような形をしています。

• 押しつぶされる際、幅（a軸）は高さ（c軸）よりもずっと速く縮みました。
• 背の高い、細長いソーダ缶を想像してください。それを押しつぶすと、側面はすぐに内側に凹むかもしれませんが、上下の部分はしばらくの間、比較的硬さを保ちます。まさにそれがこの結晶に起きたことです。圧力が上がるにつれて、結晶は「 squat（低くどっしりとした）」形へと変化していきました。

「原子価」の謎（目に見えない重み）

この物語には、隠れた登場人物がいます。それはユロピウム原子です。

• 通常の圧力下では、ユロピウムは「電荷」が約+2（これを**Eu2+**と呼びましょう）であるかのように振る舞います。
• 圧力を上げていくと、科学者たちはユロピウム原子が、まるで電荷が+3（Eu3+）であるかのように振る舞い始めていることに気づきました。
• なぜこれが重要なのか？ 電荷が+3の原子は、+2の原子よりも物理的に小さい（約10%小さい）からです。

科学者たちは、スーパーコンピュータ（DFT計算）を使用して、結晶がどのように縮むはずかを予測しました。

• 13 GPa以下： コンピュータの予測は実際の実験と完璧に一致しました。結晶は、数学が示した通りに縮みました。
• 13 GPa以上： 実際の結晶は、コンピュータの予測よりも速く縮み始めました。
• その説明： コンピュータは、ユロピウム原子が同じサイズ（Eu2+の状態）のままであると想定していました。しかし実際には、原子自体が小さくなっていた（Eu3+へと変化していた）のです。原子そのものが小さくなったため、結晶全体がコンピュータの予想よりも小さくなりました。これは、スーツケースをよりきつく詰めようと予測しているのに、中に入っている服までもが縮んでいることを忘れているようなものです！

「ゴルディロックス（適度な）」比較

論文では、この結晶の親戚にあたる EuCoGe3 および EuNiGe3 と比較を行っています。

• これらの親戚も非常によく似た挙動を示します。彼らもまた、形を崩したり壊れたりすることなく押しつぶされ、ユロピウム原子が完全に小さなバージョンへと変化することなく、ゆっくりと電荷を変化させていきます。
• これは、他の類似した結晶（Eu122系と呼ばれるもの）とは異なります。それらの結晶は、より低い圧力で全く新しい形へと突如変化したり、電荷が劇的に変化したりすることがよくあります。私たちの結晶は、このグループの中で「ゴルディロックス（絶妙なバランス）」なのです。つまり、突然の変化を起こすことなく、ゆっくりと滑らかに変化するのです。

結論

科学者たちは、磁気結晶を極限まで押しつぶし、以下のことを発見しました。

1. この結晶は非常にタフであり、35 GPaの圧力下でも形を変えたり壊れたりすることはありません。
2. 不均一に（幅の方が高さよりも速く）縮みます。
3. 高圧下でコンピュータモデルの予測よりも小さくなる理由は、内部のユロピウム原子がゆっくりとサイズを変えているためであり、これはコンピュータモデルでは完全には考慮されていなかった微妙な変化です。

要するに、この結晶は適応の達人であり、自身のアイデンティティを失うことなく、圧力の下で優雅に縮小していくのです。

技術概要

技術要約：非中心対称性を有するEuRhGe3における高圧下でのシンクロトロンX線回折およびDFT研究

問題と背景
非中心対称なBaNiSn3型構造（EuTX3）を持つ希土類間金属化合物は、磁気特性や圧力誘起現象（超伝導を含む）において豊かな性質を示す。圧力誘起の価数転移と構造崩壊は、同構造のEu122系（ThCr2Si2型）において詳細に報告されているが、そこではEu2+からEu3+への転移に伴う顕著な体積崩壊が伴う。一方、EuTGe3系列における挙動については、まだ十分に理解されていない。具体的には、EuRhGe3で見られる圧力誘起によるEu価数の緩やかな増加（約2.1から約2.4へ）が、特定の構造変化と相関しているのか、あるいは結晶格子が安定したままなのかが不明である。先行研究では、同構造化合物であるEuCoGe3およびEuNiGe3において、安定した反強磁性秩序と連続的な体積収縮が示されており、対称性の変化は見られなかったが、格子進化とEu価数揺らぎの関係を解明するための、EuRhGe3に関する系統的な構造研究は欠けていた。

手法
本研究では、実験と理論を組み合わせたアプローチを採用した：

• 実験: SOLEILのPSICHEビームラインを用い、室温において35 GPaまでのEuRhGe3に対する高圧シンクロトロン粉末X線回折（XRD）を実施した。試料は単結晶を粉末状に研磨して作成し、ヘリウム圧力伝達媒体および圧力基準としての金と共に、膜式ダイヤモンドアンビルセル（DAC）内に封入した。回折データは、格子定数および単位格子体積を抽出するためにRietveld解析（Profexソフトウェア）を用いて解析された。
• 理論: 密度汎関数理論（DFT）計算をQuantum ESPRESSOパッケージを用いて実施した。Perdew–Burke–Ernzerhof (PBE) 交換相関汎関数を用い、Euイオン上の反強磁性秩序を考慮するためにHubbard U補正（DFT+U, U = 3.8 eV）を含めた。計算は、実験データと比較するために、全圧力範囲にわたって静的なEu2+価数を仮定して行われた。
• 解析: 体積の挙動を決定するために、実験および理論の単位格子体積を3次Birch-Murnaghan状態方程式（EOS）に適合させた。これにより、体積弾性率（$B_0$）およびその圧力微分（$B'_0$）を求めた。

主な結果

1. 構造的安定性: EuRhGe3は、35 GPaまでの全圧力範囲にわたって、非中心対称な正方晶対称性（$I4mm$）を維持している。構造相転移は観察されなかった。
2. 異方性圧縮: 単位格子は異方的な圧縮を示す。格子定数 $a$ は、圧力下で $c$ よりも大幅に収縮する。この挙動は、同構造化合物であるEuCoGe3およびEuNiGe3と一致している。
3. 軸比の進化: 軸比（$c/a$）は13 GPa付近で傾きの変化を示し、圧力進化を低圧領域（LP, <13 GPa）と高圧領域（HP, >13 GPa）に分けている。$c/a$ の増加率は、LP領域の0.003 GPa$^{-1}$からHP領域の0.001 GPa$^{-1}$へと減少する。著者らは、Eu122系とは異なり、この傾きの変化は不連続な構造転移や一次の価数転移とは一致しないことを指摘している。
4. 状態方程式: 体積弾性率（$B_0$）およびその圧力微分（$B'_0$）は、それぞれ73 (1) GPaおよび5.5 (2) と決定された。これらの値は、EuCoGe3 (75.6 GPa) および EuNiGe3 (79 GPa) の値に匹敵する。
5. DFTと実験の比較:
   • 低圧域 (<13 GPa): 実験による格子定数および体積は、DFT計算と良好な一致を示した。
   • 高圧域 (>13 GPa): 実験による単位格子体積がDFT予測値よりも小さくなるという偏差が生じた。著者らは、この乖離を、平均Eu価数の圧力誘起による増加（約2.4へ接近）に起因すると考えており、これがイオン半径を減少させている。DFT計算は固定されたEu2+価性を仮定していたため、価数揺らぎに伴うこの体積減少を捉えることができなかった。

意義および主張
本論文は、EuRhGe3が、価数状態の連続的な進化にもかかわらず、構造相転移を起こすことなく、高圧下で単位格子体積の滑らかで連続的な収縮を起こすことを確立した。本研究は、$c$ 軸の圧縮率の変化（$c/a$ の傾きの変化に反映される）が、Eu122系で見られる価数転移のダイナミクスとは独立して起こることを強調している。

主な意義は、BaNiSn3型のEuTGe3系列に対する圧力への構造的応答を明らかにすることにある。著者らは、圧力下でEu価数が上昇するものの、ThCr2Si2型の化合物で見られるような破滅的な体積崩壊や対称性の破れを引き起こさないことを示している。さらに、高圧下における実験とDFTの結果の乖離は、顕著な価数揺らぎを伴う系に適用される際の、静的な価数モデルを用いたDFTの限界を示す指標となっており、これらの間金属の高温圧挙動をモデリングする際には価数変化を考慮する必要があることを裏付けている。