Sound velocity and elastic properties of Fe–Ni–S–Si liquid: the effects of pressure and multiple light elements

Fe–Ni–S–Si alloy is considered to be one of the plausible candidates of Mercury core material. Elastic properties of Fe–Ni–S–Si liquid are important to reveal the density profile of the Mercury core. In this study, we measured the P-wave velocity (VP) of Fe–Ni–S–Si (Fe73Ni10S10Si7, Fe72Ni10S5Si13, and Fe67Ni10S10Si13) liquids up to 17 GPa and 2000 K to study the effects of pressure, temperature, and multiple light elements (S and Si) on the VP and elastic properties. The VP of Fe–Ni–S–Si liquids are less sensitive to temperature. The effect of pressure on the VP are close to that of liquid Fe and smaller than those of Fe–Ni–S and Fe–Ni–Si liquids. Obtained elastic properties are KS0 = 99.1(9.4) GPa, KS’ = 3.8(0.1) and ρ0 =6.48 g/cm3 for S-rich Fe73Ni10S10Si7 liquid and KS0 = 112.1(1.5) GPa, KS’ = 4.0(0.1) and ρ0=6.64 g/cm3 for Si-rich Fe72Ni10S5Si13 liquid. The VP of Fe–Ni–S–Si liquids locate in between those of Fe–Ni–S and Fe–Ni–Si liquids. This suggests that the effect of multiple light element (S and Si) on the VP is suppressed and cancel out the effects of single light elements (S and Si) on the VP. The effect of composition on the EOS in the Fe–Ni–S–Si system is indispensable to estimate the core composition combined with the geodesy data of upcoming Mercury mission.


Introduction
The cores of terrestrial planets such as Mercury, Earth, and Mars are considered to consist of Fe-Ni alloy with some amount of light elements, such as S, Si, C, O, and H (e.g., Hauck et al. 2013;Hirose et al. 2013;Yoshizaki and McDonough 2020). Mercury is thought to have large core fraction (58-68 mass%), which is much larger than those of other terrestrial planets, from the bulk density and moment of inertia of Mercury (e.g., Hauck et al. 2013) and have a molten core based on Earth-based radar measurements (Margot et al. 2007). X-ray spectrometer in MESSENGER (MErcury Surface ENvironment GEochemistry and Ranging) spacecraft found that Mercury surface has high Mg/Si ratio and contains 1-4 wt% of sulfur (Nittler et al. 2011). High Mg/Si ratio reflects Mercury formed from highly reduced precursor materials, suggesting that Mercury core contains substantial Si as well as S (Smith et al. 2012). Furthermore, if the S abundance of Mercury surface (1-4 wt%) is representative of bulk silicate composition, metal-silicate partition experiments suggest that the core-forming melts in early stage magma ocean contain both S and Si (Chabot et al. 2014). Therefore, Mercury's core most likely contains both S and Si.
In the Fe-S-Si system, there are immiscible region, i.e., miscibility gap, in the liquid phase at ambient pressure. The miscibility gap shrinks with increasing pressure (Morard and Katsura 2010). The estimated possible range of Mercury core composition in the Fe-S-Si system corresponds to be the miscible region at the core conditions (6-40 GPa) as shown in Fig. 1 (Chabot et al. 2014). Therefore, the equation of state (EOS) of the Fe-Ni-S-Si miscible liquid is required for understanding the size and composition of Mercury core combined with observed geodesy data.
The effects of pressure and temperature on the density and sound velocity of Fe-Ni-S-Si liquid are needed to determine its EOS. The effect of a single alloying element on the P-wave velocity (V P ) of Fe and Fe-Ni liquid has been studied at 1 atm and high pressures in various systems (Fe-(Ni)-S, Fe-(Ni)-C, Fe-(Ni)-Si) (Nasch et al. 1997;Nishida et al. 2013Nishida et al. , 2016Nishida et al. , 2020Jing et al. 2014;Kuwabara et al. 2016;Shimoyama et al. 2016;Nakajima et al. 2015;Kawaguchi et al. 2017;Terasaki et al. 2019). Ni reduces the V P of pure Fe liquid (Kuwabara et al.2016). Si increases the V P of Fe-Ni liquid (Terasaki et al. 2019). On the contrary, the V P of pure Fe and Fe-Ni liquid decreases with increasing S content (Nishida et al. 2013(Nishida et al. , 2016(Nishida et al. , 2020Jing et al. 2014;Kawaguchi et al. 2017;Terasaki et al. 2019). However, the effect of simultaneous dissolution of S and Si in liquid Fe and Fe-Ni on the V P or elastic property has not been investigated. In this study, we measured the V P of liquid Fe-Ni-S-Si up to 17 GPa to study the effect of both pressure and multiple light elements (S and Si) on the V P and elastic properties and discuss the EOS of Fe-Ni-S-Si liquid.

Experimental
We used Fe-Ni-S-Si sample with three different compositions (S-rich: Fe 73 Ni 10 S 10 Si 7 , Si-rich: Fe 72 Ni 10 S 5 Si 13 , and S&Si-rich: Fe 67 Ni 10 S 10 Si 13 in at%). The samples were powder mixtures of Fe (99.9% in purity, Wako Pure Chemical Co., Ltd), Ni (99.9%, Rare Metallic Co., Ltd), FeS (99.9%, Rare Metallic Co., Ltd), and FeSi (99.9%, Rare Metallic Co., Ltd). The used compositions are in the range of possible Mercury core composition suggested by Chabot et al. (2014) as shown in Fig. 1. S-rich and Si-rich compositions of this study locate near the liquid miscible-immiscible boundary at ambient pressure and S&Sirich composition locates at the immiscible boundary at 6 GPa. Since the liquid immiscible region shrinks at higher pressure (Morard and Katsura. 2010), the sample used in this study melts into a miscible liquid at the pressure conditions of this study (S-rich: 3-16 GPa, Si-rich: 4-13 GPa, S&Si-rich: 17 GPa).
The P-wave velocity (V P ) measurements at high pressure were performed using 180 ton DIA-type cubic multianvil press (SMAP-180) at BL22XU beamline (Shobu et al. 2007) or using 1500ton Kawai-type multianvil press (SPEED-Mk. II) at BL04B1 (Higo et al. 2009), SPring-8 synchrotron facility. DIA-type press with an anvil top size of 6 mm was used in the measurement below 5 GPa and Kawai-type press with a truncated edge length of 2nd stage anvil of 5 mm was used in the measurement above 5 GPa. Details of used cell assemblies are described in Terasaki et al. (2019). The sample was sandwiched by the single-crystal sapphire buffer rod for V P measurement. High temperature was generated using cylindrical graphite heater below 5 GPa and TiC-Al 2 O 3 composite heater above 5 GPa. The experimental pressures and temperatures were determined from lattice volumes of two pressure markers (BN and MgO below 5 GPa and NaCl and MgO above 5 GPa) with their equations of state (BN: Wakabayashi and Funamori. 2015, MgO: Tange et al. 2009, NaCl: Matsui 2009). The lattice volumes were determined from X-ray diffraction (XRD) measurement. The uncertainties of experimental pressures and temperatures are derived from errors of lattice volumes of the pressure markers. At BL22XU, a monochromatic X-ray of 35 keV was used and XRD was collected using a flat panel detector (C7942CK-12, Hamamatsu Co Ltd.) with an angle-dispersive method. At BL04B1, white X-ray was used and XRD was collected using a Ge solid-state detector with an energy-dispersive method at 2θ angle of 6°. In the present study, all the V P measurements were performed in the molten state of the sample. The sample state was checked from disappearance of sample diffraction peaks and appearance of diffuse scattering signal. The experimental conditions of this study are listed in Table 1. Fe 73 Ni 10 S 10 Si 7 Fe 72 Ni 10 S 5 Si 13 and Fe 72 Ni 10 S 10 Si 13 in atomic ratio correspond to Fe-11.3wt%Ni-6.2wt%S-3.8wt%Si, Fe-11.4wt%Ni-3.1wt%S-7.1wt%Si, Fe-11.7wt%Ni-6.4wt%S-7.3wt%Si, respectively. Dashed curves show the miscible/immiscible boundary in the Fe-S-Si system at 1 atm, 6, and 10 GPa (Morard and Katsura 2010). Shaded area denotes possible range of Mercury core composition with producing 1-4 wt% S in the silicate, as measured on the Mercury surface (Chabot et al. 2014) Page 3 of 10 19 V P was measured using ultrasonic pulse-echo overlapping method (e.g., Higo et al. 2009). The V P is determined from V P = 2L/∆t, where L is a sample length and ∆t is two-way travel time of the ultrasonic wave in the sample. LiNbO 3 transducer was attached on the back side of the tungsten carbide anvil to generate and receive ultrasonic waves. The electric signals of sine waves with frequencies of 37-55 MHz were formed using a waveform generator (AFG3251C, Tektronix Inc.) and were input to the transducer. The ultrasonic echo signals from the sample interfaces were detected by the high resolution digital oscilloscope (DPO5054 and MDO4104C, Tektronix Inc.) with a sampling rate of 5 × 10 9 sample/s. Figure 2(a)(b) is typical example of measured echo signals from the sample interfaces. The difference in arrival times of echo signals from the sample front and back interfaces corresponds two-way travel time in the sample (∆t). Details of waveform analysis are described elsewhere (Kono et al. 2012). The sample length was measured from X-ray radiography image (inset of Fig. 2(a)(b)). X-ray radiography image was observed using CMOS camera (orca-flashV4, Hamamatsu Co Ltd.) with a YAG scintillator. The pixel sizes of the image were 1.5-3.0 μm.
In order to obtain constraint for equation of state of the liquid Fe-Ni-S-Si samples, we also measured sample densities at a certain pressure. The density was measured using X-ray absorption method at BL22XU beamline (e.g., Shimoyama et al. 2016). In this method, intensities of incident X-ray (I 0 ) and transmitted X-ray (I) were measured using ion chambers and X-ray absorption profile along the radial direction of cylindrical sample were obtained as shown in Fig. 2(c). Then, following Beer-Lambert law is used to deduce sample density; where μ, ρ and t denote mass absorption coefficient, density, and thickness in X-ray direction, respectively. The ρ of sample was obtained by fitting the measured X-ray absorption profile of both sample and sample capsule using Eq. (1). The third "Environment" term represents the contribution of the materials surrounding the sample capsule. The second and third terms of Eq. (1) can be obtained from X-ray intensity at the sample/capsule boundary. The μ of the sample was obtained from X-ray absorption profiles of the cold pressed solid samples and their densities determined from XRD. The μ of the capsule was referred from NIST Standard Reference database (https:// doi. org/ 10. 18434/ T4HS32). Details of X-ray absorption method are described elsewhere Terasaki and Nishida 2018).
Recovered samples were mounted in epoxy resin and polished for textural observation and chemical composition analysis. Textural observation and chemical analysis were carried out using a scanning electron microscope (SEM, JSM-6010LA, JEOL Ltd.) equipped with energy-dispersive spectroscopy (EDS) detector (UltraDry, Thermo Fischer Scientific K.K.) installed at Osaka university. (1)

Results and discussion
The effects of temperature and pressure The experimental conditions of this study and measured P-wave velocity (V P ) and density of liquid Fe-Ni-S-Si are summarized in Table 1. The error of the V P is derived from error of the travel time and that of sample length. The error of travel time is 0.3-0.7%, caused by overlapping echo signals between sample front and back. The error of sample length is 1.3-7.9%, caused by sharpness of a boundary between sample and buffer rod and also by variation in sample length. Figure 3 shows a backscattered electron image of the recovered Fe-Ni-S-Si sample from 13 GPa and 1910 K. The sample shows homogeneous quench texture, indicating that the sample was a miscible liquid at the present experimental conditions. A boundary between the sample and single-crystal Al 2 O 3 buffer rod was clear as shown in Fig. 3. Although small bits of Al 2 O 3 were observed near the boundary between the sample and buffer rod, Al 2 O 3 component was not detected in Fe-Ni-S-Si sample from EDS measurement suggesting that there was no chemical reaction between the sample and buffer rod. The effect of temperature on the V P of two sample compositions (S-rich: Fe 73 Ni 10 S 10 Si 7 , Si-rich: Fe 72 Ni 10 S 5 Si 13 ) is plotted in Fig. 4. The temperature dependence on V P is quite small in both S-rich and Si-rich liquids. It has been reported that the V P of Fe-10 wt%Ni is insensitive to temperature ) and that of Fe 57 S 43 and Fe 84 S 16 nearly independent of temperature (Nishida et al. 2013;Jing et al. 2014). Although liquid Fe-Ni-Si (Si = 6-20 wt%) and Fe has negative gentle temperature dependence (Fe-Ni-Si: The V P of liquid Fe-Ni-S-Si are plotted as a function of pressure in Fig. 5. The V P of both S-rich and Si-rich Fe-Ni-S-Si liquids increase with pressure. The pressure dependence, dV P /dP, of Fe-Ni-S-Si liquids (both S-rich and Si-rich) are smaller than those of Fe-Ni-S and Fe-Ni-Si liquids and close to that of Fe (Anderson and Ahrens 1994) and Fe-Ni 28 -S 25 (Kawaguchi et al. 2017) liquids as shown in Fig. 5. The V P of Si-rich liquid is larger than that of S-rich liquid in present pressure range. It has been reported that addition of S decreases the V P , while addition of Si increases it in the present pressure range at least below 10-14 GPa (e.g., Jing et al. 2014;Terasaki et al. 2019). Compared to previous studies, V P of Fe-Ni-S-Si liquids of this study locate in between the V P of Fe 60 Ni 10 S 30 and Fe 52 Ni 10 Si 38 liquids (Terasaki et al. 2019) as shown in a Fig. 5. The V P of S-rich liquid of this study is approximately 2% smaller and that of Si-rich liquid is ~ 3% larger than the V P of liquid Fe-10wt%Ni at 4.3-4.4 GPa . Thus, the effect of S or Si on the V P of liquid Fe-Ni-S-Si  is suppressed compared to those of liquid Fe-Ni-S and Fe-Ni-Si. This suggests that the effect of S and Si on the V P could be cancelled out each other. Detailed effect of light elements will be discussed in later section.

Elastic properties
To consider the pressure dependence of the V P and estimate elastic properties of the liquids, we fitted the V P data using second-order and third-order Birch-Murnaghan equations of state (EOS) as shown below. Since the effect of temperature on the V P of the present liquids are negligibly small, the V P data at all the temperature conditions are fitted together. The relationship between V P and density (ρ) in liquid is expressed as, where K S is adiabatic bulk modulus . The third-order Birch-Murnaghan EOS (3BM) is described as  Nishida et al. (2020) where K S0 and K S ' are adiabatic bulk modulus at ambient pressure and its pressure derivative, respectively. Then, from Eqs. (2) and (4), V P can be expressed as a function of pressure. The second-order Birch-Murnaghan EOS (2BM) can be obtained if we assume K S0 ' = 4. In this study, we directly measured density at low pressures (3-5 GPa, Table 1) and this measured density is used as a reference density (ρ ref ) at reference pressure instead of ρ 0 at ambient pressure in Eqs.
(3),(4) (e.g., Sata et al. 2002). Then, we fitted the V P data using 2BM or 3BM EOS with ρ ref to obtain K S0 and K S '. Obtained elastic properties and those of previous studies are listed in Table 2. In terms of K S0 , it is reported that K S0 of Fe-Ni reduces (approximately 63%) by addition of 30 at%S, while it stays constant (~ 1%) even by addition 38 at% Si (Terasaki et al. 2019). To understand the effect of S and Si on the K S0 correctly, we need to compare the K S0 with using same K S ' since there is tradeoff relation between K S0 and K S '. Figure 6 shows the K S0 of this study plotted as a function of difference between S and Si contents (X S -X Si ), i.e., relative concentration of light elements. The K S0 of previous studies were re-estimated by fitting the data with K S ' = 4 and also plotted in Fig. 6. The K S0 tends to decrease significantly with increasing S content and increases with increasing Si content. It is found that Fe-Ni-S-Si liquid of this study align in the (4) same trend of the Fe-Ni liquids with single light element (Fe-Ni-S and Fe-Ni-Si). This suggests a variation of K S0 of Fe-Ni-S-Si liquid simply follows concentrations of S and Si.

The effect of light elements
To evaluate the effect of two light elements (S and Si) alloying on the density and V P quantitatively, we compare the results of this study with reported mixing model of Fe-FeS binary liquid alloy system. In the ideal mixing model, volume (V) of liquid alloy is described as, where V i and x, respectively, represent a partial molar volume of end-member component i (Fe and FeS) and a mole fraction of end-member component 2. Using the mixing model, the density (ρ), thermal expansivity (α), and isothermal bulk modulus (K T ) of liquid alloy can be expressed using those of the end-member component as follows (e.g., Rivoldini et al. 2011), where M i is molar mass of end-member component i. Then, the V P can be calculated from the K S (= (1 + αγT)K T , where γ denotes Grüneisen parameter) and ρ using Eq. (2). Used parameters of the end-member components in the mixing model are referred from Table 5 in Terasaki et al. (2019). The V P as a function of S content of the Fe-FeS liquids at 1 atm and 20 GPa and at 1811 K, calculated based on the mixing model, are shown in Fig. 7. The V P from the Fe-FeS mixing model at 20 GPa is consistent with that of Fe 57 S 43 liquid measured at 20 GPa (Nishida et al. 2020). The V P of Fe-Ni-S-Si liquids at 1 atm and 20 GPa and at 1811 K, calculated by Eqs. (3) and (4) using the elastic data from this study (Table 2), are also plotted. At 1 atm, the V P gap between the Fe-Ni-S-Si liquid and Fe-FeS binary liquid are quite large. At 20 GPa, the V P of Fe-Ni-S-Si liquid becomes closer to the Fe-FeS binary liquid. Since the V P gap between the Fe-Ni-S-Si and Fe-FeS binary liquids could be an influence of Si dissolution, the results may suggest that the effect of Si on the V P of Fe-Ni-S liquid becomes smaller at higher pressure.
In the local structure of liquid Fe-alloys, S incorporates as an interstitial atom in Fe lattice, while Si substitutes the Fe site up to at least 5 GPa (Shibazaki and Kono 2018). Below 16 at% of S, S incorporates interstitially, increasing neighboring atomic distances and causing softening of the Fe-S liquid. On the contrary, when Si substitutes Fe atoms, The effect of S content on the V P at 1 atm (red) and 20 GPa (blue) at 1811 K. Dashed curves represent Fe-FeS ideal mixing model. Circle and square symbols, respectively, indicate the data of Fe 73 Ni 10 S 10 Si 7 and Fe 72 Ni 10 S 5 Si 13 calculated from the obtained EOS in this study (see Table 2) (Red open and blue closed symbols denote the data at 1 atm and 20 GPa). Parameters used in the mixing model are given in Table 5 in Terasaki et al. (2019) the nearest and 2nd neighbor distances are shortened causing stiffening of the Fe-Si liquid. Existence of liquid immiscibility in the Fe-S-Si system could be attributed to the strong covalency of Fe-S bond, rendering difficult S incorporation to Fe-Si liquid structure (Morard et al. 2008). Even in the Fe-Ni-S-Si liquid in the miscible region, the strong interaction between Fe and S is also likely to exist and the interaction between S and Si to be weak. The K S0 of Fe-Ni-S-Si liquids decrease with increasing the relative concentration of light elements, X S -X Si (Fig. 6). If the Fe-Ni-S-Si liquid is a mixture of Fe-Ni-S and Fe-Ni-Si liquids, this trend can be interpreted by the idea that the elastic moduli of an aggregate of isotropic materials can be expressed as an integral of the modulus of elasticity of each component through volume (e.g., Hill 1952).

Implication to mercury core
Fe-Ni-S-Si alloy is the most plausible candidate that consists of the Mercury core (Malavergne et al. 2010;Smith et al. 2012;Hauck et al. 2013;Chabot et al. 2014). The size and chemical composition of the Mercury core can be constrained by searching the EOS of core material that matches the observed planet mass and moment of inertia (MOI). In previous studies, the size and composition of the Mercury core were estimated based on the EOS of binary Fe-alloys with using observed gravity field and spin state data. Estimated core radius and composition are 2004 ± 39 km for Fe-S core with 4.5 ± 1.8 wt%S (7.6 ± 2.9 at%S) (Rivoldini and Van Hoolst, 2013), and 1965-2050 km for Fe-S core with 4.6 + 2.5/− 2.0 wt%S (7.7 + 4.0/-3.3 at%S) or Fe-Si core with 10.5 + 3.3/− 3.7 wt%Si (18.9 + 5.2/-6.2 at%Si) (Terasaki et al. 2019). Hauck et al. (2013) estimated the size of Mercury Fe-S-Si core to be 2017 ± 29 km with outer core density ranging 3377 ± 201 kg/m 3 from Monte Carlo internal structure model. In this study, we determined the EOSs of Fe-Ni-S-Si (Fe 73 Ni 10 S 10 Si 7 and Fe 72 Ni 10 S 5 Si 13 ) liquids, which are in the suggested composition range of the Mercury core (Chabot et al. 2014). Based on these results, the effect of multiple light elements (S and Si) on the V P and K S0 tends to be cancelled out each other in the range of S and Si contents of this study. This suggests that the V P and K S0 are similar to those of liquid Fe and Fe-Ni. Therefore, if both S and Si are dissolved into the core, it could be difficult to identify the light element from V P and elastic properties. On the other hand, if the amounts of S and Si are considerably different each other, i.e., the amount of S or Si is much more than the other light element, the difference in the V P and elastic properties between S-rich liquid and Si-rich liquid could be identified. According to the possible compositional range of the Mercury core (Chabot et al. 2014), i.e., shaded area in Fig. 1, the Mercury core composition also covers the region where a large difference between S and Si amounts, such as S-rich side (S > 20 wt% [29 at%], Si < 5 wt% [8 at%]), or Si-rich side (S < 2 wt% [3 at%], Si > 10 wt% [18 at%]). If the Mercury core is in these compositional ranges, we can probably distinguish whether S-rich or Si-rich core from V P and elastic properties. The EOS of Fe-Ni-S-Si liquids with these compositions are required to be determined in future. In addition, other pairs of multiple light elements, such as S-C, Si-C, S-H, and Si-H, are still unclear for influence on the V P and elastic properties.
If we can obtain the EOS of Fe-Ni-light elements liquids in arbitrary compositions, it is possible to strictly constrain the size and composition of the Mercury core by combining the measured data of Mercury mass and MOI. In late 2025, two BepiColombo space craft will arrive at Mercury. In this mission, more detailed and precise data of Mercury mass and MOI are also expected to be measured. To constrain the core size and composition by using these latest dataset, we need to establish the compositional dependence of the EOS of liquid Fe-Ni with multiple light elements in near future.

Conclusion
We measure the V P of Fe-Ni-S-Si (Fe 73 Ni 10 S 10 Si 7 , Fe 72 Ni 10 S 5 Si 13 , Fe 67 Ni 10 S 10 Si 13 ) liquids up to 17 GPa and 2000 K using the ultrasonic pulse-echo overlapping method with the multianvil press. The effect of temperature on the V P is quite small in Fe-Ni-S-Si liquids. This trend is similar to Fe-S liquid. The effect of pressure on the V P of Fe-Ni-S-Si liquids are smaller than those of Fe-Ni-S and Fe-Ni-Si liquids and close to that of liquid Fe. Obtained elastic properties using 3rd-order Birch-Murnaghan EOS are K S0 = 99.1(9.4) GPa, K S ' = 3.8(0.1) and ρ 0 =6.48 g/cm 3 for S-rich Fe 73 Ni 10 S 10 Si 7 liquid and K S0 = 112.1(1.5) GPa, K S ' = 4.0(0.1) and ρ 0 =6.64 g/ cm 3 for Si-rich Fe 72 Ni 10 S 5 Si 13 liquid. In the present pressure range, the V P of Fe-Ni-S-Si liquids of this study locate in between those of Fe 60 Ni 10 S 30 and Fe 52 Ni 10 Si 38 liquids and the V P of Si-rich liquid is larger than that of S-rich liquid.
Comparing with the binary Fe-FeS ideal mixing model, the V P of Fe-Ni-S-Si liquid deviate from that of Fe-FeS binary liquid at 1 atm and difference becomes to be small at 20 GPa. The K S0 of Fe-Ni-S-Si liquids decrease with increasing relative concentration of S and Si. This tendency is likely to reflect the mechanical mixing of Fe-S and Fe-Si liquids deriving from strong covalency of Fe-S bond.
If S and Si are dissolved simultaneously into the Mercury core in the compositional range of this study, the effect of S and Si on the elastic properties can be cancelled out each other, suggesting that it could be difficult to identify the light element in the core from V P and elastic properties. However, if the amounts of S and Si are considerably different each other, we can probably distinguish whether S-rich or Si-rich core. Thus, we need to determine the EOS of Fe-Ni-S-Si liquids in the compositions of S-rich side (S > 20 wt%, Si < 5 wt%), and Si-rich side (S < 2 wt%, Si > 10 wt%) to obtain the EOS in the Fe-Ni-S-Si system which covers entire range of possible Mercury core composition. This can be useful dataset to estimate the core composition combined with the geodesy data of upcoming Mercury mission.