Melting relation of FeS-H system under high pressure: Implications for the core of Ganymede Earth & Planetary Science 102 Ganymede is one of the Galilean satellites of Jupiter, and is the largest moon in our solar system. The observation data of the Galilean mission suggest that Ganymede has an intrinsic magnetic field and the interior is strongly differentiated into an outermost H 2 O layer, a silicate rock mantle, and an iron-rich core [1]. The core is considered to consist of iron (Fe) and light elements. Sulfur (S) is one of the plausible light elements in the core because iron sulfide is found in many meteorites. Hydrogen (H) is also most probably present in the core. This is because, before the differentiation occurred in an early period, the silicate rock and the iron alloy mixed with and then may have reacted with H 2 O in the interior of Ganymede. As a result, it may well be that the silicate rock layer contains a significant amount of water, and the iron alloy may contain hydrogen. It is known that hydrogen significantly lowers the melting temperature of metals. However, it has not been reported whether hydrogen affects the melting temperature of Fe–S alloy. If hydrogen lowers the melting temperature of Fe–S alloy, the previous models of Ganymede’s core (e.g. [2]) may need to be revised. Therefore, in order to constrain the state and degree of melting of Ganymede’s core, we first determined the melting temperatures of the FeS–H system (FeSH x ) at pressures up to 16.5 GPa [3]. Energy-dispersive X-ray diffraction experiments were carried out at high temperature and high pressure using a Kawai-type multianvil apparatus at BL04B1 . The diffraction patterns were collected for a period of 120 or 300 s at each 50–100 K temperature step and at high pressure. In this study, melting of the sample was determined on the basis of observations of the disappearance of the diffraction peaks from the sample and the appearance of a broad halo peak. The FeS sample powder was packed into a NaCl container along with LiAlH 4 , because NaCl is reported to effectively seal hydrogen at high pressures and temperatures. Hydrogen was supplied to the sample from the thermal decomposition of LiAlH 4 . The observed unit cell volumes of FeS and FeSH x with increasing temperature at 16.5 GPa are shown in Fig. 1. On the basis of the present X-ray diffraction measurements, the low-temperature phases of FeS and FeSH x have a hexagonal NiAs-type superstructure (Phase IV) and the high-temperature phase has a simple NiAs-type structure (Phase V). This is consistent with the previously reported phase relations of pure FeS [4]. However, above 800 K, the observed volume was clearly larger than that of pure FeS up to the melting point as shown in Fig. 1. This expansion in volume can be explained by hydrogenation, i.e., synthesis of FeSH x . The drastic volume increase at high temperature was observed throughout the pressure range studied, except at 1.9 GPa. We considered that the structure of FeSH x is the same as that of FeS whereas the volume of FeSH x is larger than that of FeS owing to the dissolution of hydrogen in the interstitial sites of FeS. 22 23 24 25 26 27 400 600 800 1000 1200 1400 1600 1800 16.5 GPa FeSH x FeS FeS IV (Urakawa et al ., 2004) FeS V (Urakawa et al ., 2004) FeSH x IV (this study) FeSH x V (this study) Temperature (K) Volume per FeS (Å 3 ) 400 600 800 1000 1200 1400 1600 1800 2000 0 5 10 15 20 Pressure (GPa) Temperature (K) 0.02 0.2 0.2 0.2 0.2 0.4 15 20 FeSH x IV FeSH x IV FeSH x V melt melt melting curve of FeSH x melting curve of FeS melt (FeS) FeSH x V × × × × × × × × Fig. 1. Thermal expansion of FeS and FeSH x at 16.5 GPa. The volume per FeS is defined as V (unit cell volume) divided by Z (where Z = 8 for phase IV, and Z = 2 for phase V). The solid diamonds and squares correspond to FeSH x IV and FeSH x V, respectively. The open diamonds and squares correspond to FeS IV and FeS V, respectively, from Ref. [4]. Fig. 2. Pressure-temperature phase diagram of FeSH x , together with Boehler’s melting curve of FeS [5]. The solid diamonds and squares correspond to FeSH x IV and FeSH x V, respectively, as determined in this study. The open circles correspond to the FeSH x –melt. The crosses denote the observed melting temperatures of FeS. The numbers on the melting curve show the hydrogen concentration ( x ) just before the melting of FeSH x , where x is defined as the number of hydrogen atoms per FeS in a unit cell. 103 The observed phase relationships of FeSH x are plotted on the pressure-temperature diagram shown in Fig. 2, together with the melting temperature of pure FeS [5]. The melting temperature of FeSH x was clearly decreased by 150 – 250 K compared with that of FeS between 3.0 and 16.5 GPa owing to hydrogenation. The melting temperature of FeSH x exhibited a minimum and the concentration of hydrogen increased at about 3 GPa ( x ≈ 0.2 at 3.0 GPa, whereas x ≈ 0.02 at 1.9 GPa). This indicates that hydrogen starts to dissolve into the interstitial sites of FeS between 1.9 and 3.0 GPa. As a result, the depression of the melting temperature at about 3.0 GPa occurred. Figure 3 shows the predicted phase diagram of the FeH x –FeSH x system at 5 GPa with model temperatures equivalent to those of Ganymede’s core. We assume that the FeH x –FeSH x system has the same eutectic relation as the Fe–FeS system. Since the melting temperatures of Fe and FeS are depressed due to hydrogenation, the liquidus of the FeH x –FeSH x binary system is lower than the core temperatures over a wide range of sulfur content in the core (Region I in Fig. 3). Then, Ganymede’s core could be completely molten in Region I. If the core contains much sulfur (Region II in Fig. 3), it could be partially molten (solid FeSH x coexisting with liquid Fe–S–H) because the liquidus is higher than the core temperatures. Hauck et al . [2] showed that a completely molten core cannot sustain the Ganymede dynamo up to the present day by thermal convection. Thus, the partially molten core (Region II in Fig. 3) is more reasonable for sustaining the dynamo, suggesting that the core contains more than about 33 wt.% S. The present study suggests that Ganymede’s core could have a solid FeSH x inner core and liquid Fe–S–H outer core if the core consists of the Fe–S–H system (Fig. 4). The results of this study strongly indicate that hydrogen is certainly important for understanding the internal structures of planets and satellites outside the snow line, although additional experiments for the FeH x –FeSH x system are needed in order to discuss the core in more detail. 1000 1200 1400 1600 1800 2000 2200 0 5 10 15 20 25 30 35 Bulk sulfur content in core (wt.%) Temperature (K) FeSH x Liquid 5 GPa FeH x + Liq. FeH x Core temperatures Depression of melting temperature due to hydrogenation I II + Liq. FeSH x Silicate mantle Ice Ice Liquid outer core Fe-S-H Solid inner core FeSH x I n t e r n a l o c e a n Fig. 3. Predicted phase diagram of the FeH x –FeSH x system at 5 GPa, assuming it has the same eutectic relation as the Fe–FeS system. The open and solid squares denote the melting temperatures of Fe (Shen et al ., 1998) and FeH x (Fukai et al ., 2003), respectively. The open and solid triangles denote those of FeS (Boehler, 1992) and FeSH x (this study) [3], respectively. The solid and dotted lines denote the predicted liquidus. The broken line denotes Ganymede’s core temperatures with various bulk sulfur contents in the core, taken from Kimura et al . (2009) (1030 km H 2 O layer). The open circle and diamonds denote the core temperatures reported by Hauck et al . (2006) (S = 10 wt.%) and Bland et al . (2008) (S = 1 and 10 wt.%), respectively . Fig. 4. Schematic figure of the internal structure of Ganymede with Fe–S–H core. References [1] J.D. Anderson et al .: Nature 384 (1996) 541. [2] S.A. Hauck et al .: J. Geophys. Res. 111(2006) doi:10.1029/2005JE002557. [3] Y. Shibazaki, E. Ohtani, H. Terasaki, R. Tateyama, T. Sakamaki, T. Tsuchiya, K. Funakoshi: Earth Planet. Sci. Lett. 301 (2011) 153 . [4] S. Urakawa et al .: Phys. Earth Planet. Inter. 143-144 (2004) 469. [5] R. Boehler: Earth Planet. Sci. Lett. 111 (1992) 217. Yuki Shibazaki a,†, *, Eiji Ohtani a and Hidenori Terasaki b a Department of Earth and Planetary Material Sciences, Tohoku University b Department of Earth and Space Science, Osaka University *E-mail: yshibazaki@ciw.edu † Present Address: Geophysical Laboratory, Carnegie Institution of Washington, USA