日本地球惑星科学連合2023年大会

講演情報

[E] 口頭発表

セッション記号 S (固体地球科学) » S-IT 地球内部科学・地球惑星テクトニクス

[S-IT18] 惑星中心核:内部構造・形成・進化

2023年5月26日(金) 13:45 〜 15:00 102 (幕張メッセ国際会議場)

コンビーナ:飯塚 理子(東京大学大学院理学系研究科地殻化学実験施設)、寺崎 英紀(岡山大学理学部)、大谷 栄治(東北大学大学院理学研究科地学専攻)、William F McDonough(Department of Earth Science and Research Center for Neutrino Science, Tohoku University, Sendai, Miyagi 980-8578, Japan)、座長:寺崎 英紀(岡山大学理学部)、飯塚 理子(東京大学大学院理学系研究科地殻化学実験施設)


14:45 〜 15:00

[SIT18-05] First-principles calculations and neutron diffraction measurements on hydrogenation of FeS at high-pressure conditions

*高野 将大1鍵 裕之1森 悠一郎1柿澤 翔2佐野 亜沙美3土屋 卓久4 (1.東京大学大学院理学系研究科、2.高輝度光科学研究センター、3.日本原子力研究開発機構 J-PARCセンター 物質・生命科学ディビジョン 中性子利用セクション、4.愛媛大学地球深部ダイナミクス研究センター)


キーワード:FeS、水素化、中性子回折、第一原理計算

The cores of terrestrial planets consist mainly of iron and small amounts of light elements. Among the candidates for light elements, hydrogen is one of the most promising elements because of its significant effect on the volume expansion of iron by entering the interstitial sites to form iron hydride, FeHx.

The hydrogen-induced volume expansion coefficient ΔVH, representing the volume expansion per hydrogen atom in the unit cell, is an important parameter to estimate the amount of hydrogen present in the cores. To date, the hydrogen-induced volume expansion coefficients for hcp and fcc phases of iron have been determined by in-situ neutron diffraction experiments at high-temperature and high-pressure conditions.

However, there have been not many studies of hydrogenation of iron in ternary systems. Because the cores of terrestrial planets can also contain light elements other than hydrogen, it is necessary to elucidate how the presence of other light elements affects the hydrogen-induced volume expansion coefficient.

As an additional element in the system, we focus on sulfur, which is commonly found in primitive meteorites as troilite (FeS). In this study, the hydrogen-induced volume expansion coefficients of FeS, a typical iron sulfide, have been determined by static (T = 0 K) first-principles calculations using the Generalized Gradient Approximation (GGA) method. In the calculations, the hydrogen occupancies in the octahedral sites of phase V of FeS (NiAs-type phase, hereafter FeS-V) were fixed to 1 (FeSH2 in stoichiometry). We fully optimized all structural degree of freedom, the lattice parameters and the atomic coordinates of FeSH2 including stable hydrogen positions, up to 150 GPa.

The obtained ΔVH of FeS-V was notably larger than that of pure hcp iron and showed a clear negative correlation with pressure. In addition, the positions of the H atom notably deviated from the center of the octahedral sites. Because FeS is known to have a strongly correlated electron system at pressures up to ~70 GPa at 0 K (e.g., Terranova and Leeuw, 2017), where the iron magnetization survives, the GGA method may not be the best approach to estimate the positions and the volume expansion effect by the H atoms incorporation into FeS. Furthermore, the hydrogen incorporation might extend the stability field of iron magnetizations to higher pressures. In this case, using the GGA combined with the Hubbard U correction (GGA+U) is expected to improve the electronic and the atomic structures of hydrogenated magnetic FeS.

We plan to perform neutron diffraction measurements on hydrogenation of FeS above 10 GPa and determine experimentally the hydrogen-induced volume expansion coefficient of FeS. The results of the experiments and a comparison between the calculations and the experiments will be discussed in the presentation.