Next Article in Journal
Ductility Index for Refractory High Entropy Alloys
Previous Article in Journal
Microstructure Evolution and Mechanical Properties of TiB2/Al–Cu–Mn–Cd Composite with the Adoption of Two-Stage Solution and Aging Treatment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pressure Effects on the Thermodynamic Properties of MgSiO3 Akimotoite

by
Chang Su
1,2,3,4,
Wei Song
2,
Guang Yang
4,5,
Yonggang Liu
2,* and
Qingyi Li
1
1
School of Earth Sciences, Institute of Disaster Prevention, Sanhe 065201, China
2
Key Laboratory of High Temperature and High Pressure Study of the Earth’s Interior, Institute of Geochemistry, Chinese Academy of Sciences, Guiyang 550081, China
3
Hebei Key Laboratory of Earthquake Dynamics, Sanhe 065201, China
4
Hebei Key Laboratory of Resource and Environmental Disaster Mechanism and Risk Monitoring, Sanhe 065201, China
5
School of Information Engineering, Institute of Disaster Prevention, Sanhe 065201, China
*
Author to whom correspondence should be addressed.
Crystals 2024, 14(10), 837; https://doi.org/10.3390/cryst14100837
Submission received: 17 August 2024 / Revised: 22 September 2024 / Accepted: 25 September 2024 / Published: 26 September 2024
(This article belongs to the Section Mineralogical Crystallography and Biomineralization)

Abstract

:
The thermodynamic properties of MgSiO3 akimotoite at high temperatures and high pressures are important for investigating the phase equilibria of the Earth’s transition zone and the upper part of the lower mantle. In this paper, we present the self-consistent unit-cell volume, elastic properties, and in particular, thermodynamic properties including thermal expansion, heat capacity, entropy, and Grüneisen parameter of MgSiO3 akimotoite at pressures up to 30 GPa and temperatures to 2000 K using an iterative numerical method and available experimental data, which are consistent with the previous studies. The results show that the determined thermal expansion, heat capacity, entropy, and Grüneisen parameter exhibit a nonlinear and negative relationship with increasing pressure. Additionally, the pressure derivatives of these thermodynamic parameters along with the temperature are also presented.

1. Introduction

Akimotoite (Ak) is one of the high-pressure polymorphs of MgSiO3 with the ilmenite-type structure. It is believed to be stable between 18 and 25 GPa, and from 900 to 2200 K, and occurs mainly in cold stagnant slabs at the bottom of the mantle transition zone or in the uppermost region of the lower mantle [1,2,3,4]. Recent studies suggest that akimotoite is stable in the harzburgite layers of subducting slabs rather than in the mid-ocean ridge basalt or pyrolite [4,5]. Akimotoite is stable at the low pressure of bridgmanite (Brg) with the perovskite-type structure and at the low temperature of majorite (Mj) with the garnet-type structure in the pure MgSiO3 system, which is expected to contribute to the structure of the 660 km discontinuity [6,7]. Therefore, the elastic properties of akimotoite have been studied for decades by numerous experimental measurements [3,4,8,9,10], and theoretical calculations [11,12,13,14,15], providing a wealth of data for understanding the seismic discontinuities and density jumps caused by akimotoite-related phase transitions.
On the other hand, knowledge of the thermodynamic properties of MgSiO3 akimotoite is essential for the study of the geotherms and phase equilibria. However, compared with the sufficient elastic moduli and volume data of MgSiO3 akimotoite, the thermodynamic properties of MgSiO3 akimotoite, especially at high-pressure conditions are still lacking. At ambient pressure conditions, the thermal expansion of MgSiO3 akimotoite has been determined from the unit-cell volume by using X-ray diffraction experiments [10,16,17], while the heat capacity and entropy in a wide range of temperatures have been presented from calorimetric measurements [16,17,18]. However, thermodynamic properties such as thermal expansion, heat capacity, entropy, and Grüneisen parameter at simultaneous high temperatures and high pressures have only been obtained by First-principles calculations [11,19].
Decades ago, Davis and Gordon introduced a numerical iterative procedure by which the density, the elastic moduli, and especially the thermal expansion, and heat capacity of the metal melt, at high temperatures and high pressures, could be obtained from adiabatic elastic wave velocities [20]. Recently, the feasibility of this method has been demonstrated for solid minerals [21]. Therefore, self-consistent thermodynamic properties at high-pressure conditions can be conveniently determined from the numerous published experimental data on minerals.
In this paper, we evaluated the unit-cell volume, the elastic properties, and most importantly, the thermodynamic properties of MgSiO3 akimotoite at high temperature and high-pressure, with the combination of the available experimental data and a numerical iterative procedure. Comparisons with previous work have been made to verify the accuracy of our results. In addition, the pressure effects on the thermal expansion, heat capacity, entropy, and Grüneisen parameter at different temperature conditions were presented.

2. Methods

2.1. Calculation Procedure

The basic calculation procedure, which is described in our recent studies [21,22], is presented in the following. Note that the uncertainties of the derived parameters at different temperature and pressure conditions are derived from the propagation of the experimental measurement errors by means of a Taylor series expansion.
The thermal expansion (α) at certain pressure can be calculated with the volume (V) as a function of temperature through Equation (1).
α T = 1 V V T
And the integration of Equation (1) can be expressed as Equation (2).
V T = V 0 e x p T 0 T α T d T
here V0 refers to the volume at ambient conditions. Also, the isothermal pressure derivative of volume can be expressed as Equation (3),
V P T = V 2 1 v Φ 2 + T α 2 C P
where CP is the heat capacity. And vΦ is the bulk sound velocity, which can be determined with the P-wave velocity (vP) and the S-wave velocity (vS) via Equation (4). Equation (5) can be used to evaluate the isothermal derivative of the heat capacity with respect to pressure.
v Φ = v p 2 4 3 v s 2 1 2
C P P T = V T α 2 + α T P
The calculation process starts by using the experimentally measured volume of MgSiO3 akimotoite at ambient pressure conditions (V(P0,T)) to determine the thermal expansion as a function of temperature (α(P0,T)) by Equation (1). Meanwhile, the heat capacity of MgSiO3 akimotoite as a function of temperature (CP(P0,T)) can be fitted from published experimental data with empirical equations. Thus, with the obtained α(P0,T) and CP(P0,T), the approximate volume at an arbitrary reference pressure (V(P0P,T)) can be estimated using Equation (3). Also, the bulk velocity of MgSiO3 akimotoite as a function of temperature and pressure (vΦ(P,T)) can be obtained using the P-wave velocity and S-wave velocity at high temperatures and high pressures (vP(P,T) and vS(P,T)) via Equation (4). Therefore, the resulting V(P0P,T) can be used to update the values of α(P0P,T) and CP(P0P,T) via Equations (1) and (5), respectively. Finally, iteration of this loop results in converged volume, thermal expansion, and heat capacity of MgSiO3 akimotoite as a function of temperature and pressure.
Also, with the combinations of V(P,T), vP(P,T), vS(P,T), the adiabatic bulk modulus (KS), isothermal bulk modulus (KT), and shear modulus (G) of MgSiO3 akimotoite at simultaneous high temperatures and high pressures can be determined through Equations (6)–(8), respectively.
K S = 1 V v p 2 4 3 v s 2
K T = K S 1 + α γ T
G = 1 V v S 2
Itis worth noting that the parameter γ in Equation (7) stands for the Grüneisen parameter, which can be derived from Equation (9). Furthermore, with the entropy measured at ambient conditions, the entropy at high temperature and high-pressure conditions can also be calculated via Equations (10) and (11), respectively [23].
γ = α K S V C P
S T P = S T 0 P + T 0 T C P T d T
S P T = S P 0 T P 0 P α V d P

2.2. Thermoelastic Data of MgSiO3 Akimotoite

To date, there are only a few experimental measurements on the sound velocity of MgSiO3 akimotoite. The sound velocity of MgSiO3 akimotoite at ambient pressure conditions was determined by Weidner and Ito [8] using Brillouin spectroscopy, and the high temperature and high-pressure data were determined only by Zhou et al. [3] using ultrasonic interferometry techniques in the mantle transition zone environment.
In the paper by Zhou et al. [3], they presented a linear pressure dependence and a nonlinear temperature dependence of the sound velocity for MgSiO3 akimotoite, which is in good agreement with the result by Weidner and Ito [8]. Therefore, in this study, we used the proposed fitting equations for vP and vS by Zhou et al. [3] (Table 1), which are shown as Equations (12) and (13), respectively.
v P P , T = 10211 8 + 55.8 4 × P 0.352 12 × T 300 1.22 11 × 10 4 × T 300 2 m / s
v S P , T = 5862 6 + 21.9 3 × P 0.253 9 × T 300 1.39 8 × 10 4 × T 300 2           m / s
In general, the temperature variation of the experimentally measured isobaric volume is needed to derive the thermal expansion (Equation (1)) to proceed with the calculation. The high-temperature unit-cell volume of MgSiO3 akimotoite at ambient pressure conditions was reported to be 876 K by Ashida et al. [16], and the data were recently remeasured by Kojitani et al. [17]. Thus, the available volume data from these papers are analyzed with the EoSFit7 Program [24], and then fitted with the empirical equation proposed by Fei [25]. The obtained equation for the thermal expansion of MgSiO3 akimotoite as a function of temperature is shown in Equation (14). The calculated results via Equation (2) and Equation (14) are shown in Figure 1, which are highly consistent with the experimentally determined data, with the largest difference less than 0.2%.
α T = 2.61 7 × 10 5 + 0.54 14 × 10 8 T 1.28 14 × T 2           K 1
In the studies by Ashida et al. [16] and Kojitani et al. [17], they also reported the high-temperature heat capacity of MgSiO3 akimotoite with calorimetric measurements to ~800 K, and the results were expressed by the Maier-Kelley equation [26] and a polynomial formula shown as Equation (15), respectively. The low-temperature data were accurately determined by Akaogi et al. [18]. These heat capacity values are fitted with Equation (15), as well as Equations (16)–(18), which are three classical empirical formulas proposed by Hass and Fisher [27], Berman and Brown [28], and Richet and Fiquet [29], respectively. The determined fitting coefficients, reduced χ2, and R2 are listed in Table 2. Furthermore, the comparisons between the calculated heat capacity of MgSiO3 akimotoite using Equations (15)–(18) and the available experimental results [16,17,18] are also shown in Figure 2.
The heat capacity of MgSiO3 akimotoite at 2000 K was estimated to be 129.512 J/mol∙K by Hofmeister and Ito (Model C) [30], and 135.83 J/mol∙K by Kojitani et al. [17]. Therefore, after considering the reduced χ2 and R2 of the obtained fitting equations, as well as the previously estimated values, we chose Equation (17) in this study. The largest difference between the value calculated using Equation (17) and the results by Ashida et al. [16], Akaogi et al. [18], and Kojitani et al. [17] are 3.7%, 2.3%, and 2.5% above room temperature, respectively.
Finally, the standard entropy used in the present study was from Akaogi et al. [18] (Table 1), which was determined to be S°298 = 53.7(4) J/mol∙K.
C P = c 0 + c 1 T + c 2 T 0.5 + c 3 T 1 + c 4 T 2 + c 5 T 3 J / m o l · K
C P = c 0 + c 1 T + c 2 T 2 + c 3 T 0.5 + c 4 T 2           J / m o l · K
C P = c 0 + c 1 T 0.5 + c 2 T 2 + c 3 T 3           J / m o l · K
C P = c 0 + c 1 l n T + c 2 T 1 + c 3 T 2 + c 4 T 3           J / m o l · K

3. Results

3.1. High Temperature and High Pressure Unit-Cell Volume of MgSiO3 Akimotoite

Our determined unit-cell volume of MgSiO3 akimotoite variation with pressure at different temperatures is shown in Figure 3, together with the results carried out from X-ray diffraction experiments [3,9,10], and First-principle calculations [11]. At room temperature, our result shows good consistency with the X-ray diffraction measurements by Reynard et al. [9], and the difference is less than 0.55%. The excellent agreement is also maintained at high temperature conditions, and the largest differences are 0.43% and 0.53% compared with those of Wang et al. [10] and Zhou et al. [3], respectively. Meanwhile, our obtained unit-cell volume is generally lower than the theoretical result by Karki and Wentzcovitch [11], and the largest differences are 0.87%, 0.98%, and 1.0% at 300 K, 1000 K, and 2000 K, respectively.

3.2. High Temperature and High-Pressure Elastic Properties of MgSiO3 Akimotoite

Our derived elastic moduli of MgSiO3 akimotoite as a function of pressure and temperature are shown in Figure 4. In Figure 4b,c, the calculated adiabatic bulk modulus and shear modulus are in good agreement with the experimental results under the relevant conditions [3], and the largest differences are 2.1%, 1.1%, and 0.4% at 300 K, 700 K, and 1500 K, respectively. Compared with the results obtained from using first-principles calculations [11,14], our adiabatic bulk modulus and shear modulus are generally larger than those of Hao et al. [14], and Karki and Wentzcovitch [11] (shown in Figure 4b,d, respectively). While our shear modulus at 300 K is larger than that by Hao et al. [14], then with increasing temperature, the values become smaller than those by Hao et al. [14] at 1500 K and 2000 K.
As shown in Figure 4, all the elastic moduli generally display a linear dependence on pressure and a nonlinear temperature dependence within the uncertainties. Therefore, we adopt the data to the equation of M = M0 + (∂M/∂P) × P + (∂M/∂T) × (T − 300) + (∂2M/∂T2) × (T − 300)2, where M refers to the elastic moduli, M0 refers to the elastic moduli at ambient conditions, ∂M/∂P is the first pressure derivative at fixed temperature, and ∂M/∂T, ∂2M/∂T2 are the first and second temperature derivatives of M at fixed pressure, respectively. The fitting parameters are listed in Table 3, along with the published results for comparison.
The isothermal bulk modulus, adiabatic bulk modulus, and shear modulus at ambient conditions are determined as KT0 = 221.5 (11) GPa, KS0 = 222.4 (13) GPa, and G0 = 130.8 (5) GPa, respectively. At ambient conditions, the pressure derivatives of the elastic moduli are determined to be ∂KT/∂P = 4.49 (1), ∂KS/∂P = 4.39 (1), and ∂G/∂P = 1.54 (1), and the first temperature derivatives of the elastic moduli are determined as ∂KT/∂T = −2.943 (1) × 10−2 GPa/K, ∂KS/∂T = −1.937 (1) × 10−2 GPa/K, and ∂G/∂T = −1.645 (1) × 10−2 GPa/K, the second temperature derivatives of the elastic moduli are determined as ∂2KT/∂T2 = −0.847 (6) × 10−6 GPa/K2, ∂2KS/∂T2 = −2.060 (8) × 10−6 GPa/K2, and ∂2G/∂T2 = −5.351 (4) × 10−6 GPa/K2. These values are generally within the range of the available values listed in Table 3 [3,4,8,9,10,11,12,13,14].

3.3. High Temperature and High-Pressure Thermodynamic Properties of MgSiO3 Akimotoite

The calculated thermal expansion of MgSiO3 akimotoite as a function of temperature is presented in Figure 5. Using the unit-cell volume data presented by Ashida et al. [16], Saxena et al. [31] obtained the thermal expansion of MgSiO3 akimotoite at ambient pressure conditions. Also, the thermal expansion data has been obtained from experimentally measured entropy by Chopelas [32], as well as from First-principles calculations by Karki and Wentzcovitch [11]. Recently, Kojitani et al. [17] also presented the thermal expansion of MgSiO3 akimotoite as a function of temperature by fitting the X-ray diffraction determined unit-cell volumes. At 0 GPa, our calculated thermal expansion is determined to be α0 = 1.35 (4) × 10−5/K at 300 K, which is slightly lower than the previous results [11,17,31,32]. However, in the temperature range of ~700–2000 K, our data is in great agreement with the results by Karki and Wentzcovitch [11] and Kojitani et al. [17], within 2.9% and 3.4%, respectively, and ~5.7% and ~6.8% larger than those of Saxena et al. [31] and Chopelas [32], respectively. The good consistency also holds at high-pressure conditions and above 700 K. The average differences are ~0.7%, ~2.9%, and ~7.3% at 10 GPa, 20 GPa, and 30 GPa, respectively, compared to the theoretical results [11].
Figure 6 shows the temperature dependence of the calculated heat capacity at 0 GPa, 10 GPa, 20 GPa, and 30 GPa, along with previous results at ambient pressure conditions [16,17,31], and at high-pressure conditions [11]. The heat capacity of MgSiO3 akimotoite was measured by differential scanning calorimetry at ambient pressure conditions [16,17], and by First-principles calculations at high-pressure conditions [11]. Also, Saxena et al. [31] presented the fitting equation for the heat capacity as a function of temperature using the data from Watanabe [33]. As mentioned above, since the heat capacity data at high temperature and ambient pressure conditions used in this study were from the investigations by Ashida et al. [16], Akaogi et al. [18], and Kojitani et al. [17], our calculated heat capacity of MgSiO3 akimotoite at 0 GPa shows good agreement with these published results [16,17]. Our calculated heat capacity is also in good agreement with the results obtained by Saxena et al. [31] and Karki and Wentzcovitch [11], with the largest differences being less than 2.5% and 1.5%, respectively. Compared to the first calculation result at high-pressure conditions [11], our result is ~0.9% lower than that at 10 GPa, while ~1.3%, and ~1.4% higher than that at 20 GPa, and 30 GPa, respectively.
The entropy of minerals at high temperatures and high-pressures is a valuable thermodynamic parameter in determining reliable phase diagrams [34], and its value can be obtained from heat capacity measurements using adiabatic calorimetry. The entropy of MgSiO3 akimotoite at ambient pressure conditions has been reported by Ashida et al. [16], and Saxena et al. [31], and currently, high-pressure values have only been determined by First-principles calculations [11]. Our calculated entropy of MgSiO3 akimotoite at high temperatures and high pressures is shown in Figure 7, which is generally, in good agreement with the previous investigations.
The Grüneisen parameter is of considerable importance in geosciences, as it often appears in equations describing the thermoelastic behavior of minerals at high-pressures and temperatures [35]. With the combination of the obtained thermodynamic and elastic data, our calculated Grüneisen parameter of MgSiO3 akimotoite at high temperatures and high-pressures is illustrated in Figure 8, along with the available results at ambient pressure conditions by Saxena et al. [31], and at high-pressure conditions by Karki and Wentzcovitch [11].
At ambient conditions, the Grüneisen parameter of MgSiO3 akimotoite is determined to be γ0 = 1.03 (25), which is lower than those by Saxena et al. [31] (γ0 = 1.72), and Karki and Wentzcovitch [11] (γ0 = 1.29). As the temperature increases to ~900 K, our Grüneisen parameter increases rapidly, and then the value increases slightly with a very small temperature derivative as 0.00007/K, and the values are generally higher than those presented by Saxena et al. [31], as well as Karki and Wentzcovitch [11]. Also, the differences in the temperature dependencies below ~500 K are probably caused by the different temperature derivatives of the thermal expansion and heat capacity under the low-temperature conditions (Figure 5 and Figure 6). Under high-pressure conditions, the calculated Grüneisen parameter shows a temperature derivative as 0.00002/K at 10 GPa, −0.00002/K at 20 GPa, and −0.00005/K at 30 GPa, which shows a similar pattern to the result of Karki and Wentzcovitch [11].

4. Discussions

The pressure effects on the thermodynamic properties of minerals are prerequisites for deriving the thermal state of the Earth’s interior. Therefore, the calculated thermal expansion, heat capacity, entropy, and Grüneisen parameter of MgSiO3 akimotoite as a function of pressure at various temperatures are illustrated in Figure 9a–d, respectively.
As shown in Figure 9a–d, it is obvious that all these thermodynamic parameters tend to be negative and nonlinear with increasing pressure. Thus, we fitted our calculated results to a polynomial equation of N = N0 + (∂N/∂P)T × P + (∂2N/∂P2)T × P2 at fixed temperatures, where N refers to the thermodynamic properties, N0 refers to M at ambient pressure conditions, (∂N/∂P)T and (∂2N/∂P2)T refer to the first and second pressure derivatives of N, respectively. The fitting coefficients are listed in Table 4. Additionally, the first and second pressure derivatives of thermal expansion, heat capacity, entropy, and Grüneisen parameter of MgSiO3 akimotoite as a function of temperature are shown in Figure 10a–d, respectively.
In Table 4 and Figure 10a, the first pressure derivative of thermal expansion is ∂α/∂P = −5.84 (2) × 10−7/K∙GPa at 300 K, then the value increases to ~−5.53 × 10−7/K∙GPa at ~450 K, and decreases with the temperature to −9.42 (5) × 10−7/K∙GPa at 2000 K. The value of ∂CP/∂P is calculated to be −0.8066 (6) J/mol∙K∙GPa at 300 K, and the value increases rapidly to ~−0.2077 J/mol∙K∙GPa at ~1000 K, and then decreases slowly to −0.3178 (28) J/mol∙K∙GPa at 2000 K (Table 4, Figure 10b). The pressure effect on the entropy of MgSiO3 akimotoite is shown in Table 4 and Figure 10c, which generally increases with the temperature, and ∂S/∂P decreases from −0.315 (2) J/mol∙K∙GPa to −0.929 (4) J/mol∙K∙GPa in the temperature range of 300–2000 K. In addition, ∂γ/∂P is determined to be −2.200 (7) × 10−2/GPa at 300 K, then, similar to the heat capacity, the value tends to increase to a maximum of ~−0.979 × 10−2/GPa at ~650 K, and finally decreases to −1.643 (6) × 10−2/GPa at 2000 K (Table 4, Figure 10d). Also, with the combination of the temperature variation of ∂α/∂P and ∂CP/∂P (Figure 10a,b, respectively), it seems that the pressure derivative of the Grüneisen parameter is largely dependent on the pressure effects on the thermal expansion and heat capacity.
The phase transition boundaries and the Clapeyron slopes of the MgSiO3 Majorite-akimotoite-bridgmanite system are of great geophysical interest as they are thought to be one of the causes of the boundary between the Earth’s transition zone and lower mantle. Although plenty of work has been done to determine the phase transition diagram of this MgSiO3 system, the results are still controversial [19,36,37,38]. In the meantime, the phase transition boundaries in the MgSiO3 system are usually derived using Equation (19):
G P , T = H ° T T S ° T + 1   a t m P V P , T d P = 0
where ΔG(P,T) and ΔV(P,T) are the Gibbs free energy and the volume changes of two coexisting phases, respectively, ΔH°(T) and ΔS°(T) are the enthalpy and the entropy changes of two coexisting phases at ambient pressure conditions (1 atm) and temperature T, respectively. This equation provides a convenient way to estimate the phase transition boundaries, avoiding the use of thermodynamic properties at high-pressure conditions, which have been employed by a majority of previous investigations [19,36,37].
On the other hand, the Clausius-Clapeyron equation (Equation (20)) is another classical equation to derive the phase equilibrium:
d P d T = S P , T V P , T
here dP/dT is the Clapeyron slope, and ΔS(P,T) and ΔV(P,T) are the entropy and volume changes between the two different phases, respectively. Therefore, the Clapeyron slope of the MgSiO3 system can be determined from the entropy data of related minerals at high temperatures and high-pressures [38]. However, similar to akimotoite, the experimentally measured thermodynamic properties of majorite and bridgmanite at high temperatures and high-pressures have only been presented under ambient pressure conditions. Therefore, more data is needed for further studies.

5. Conclusions

To summarize, we used a numerical iterative method and published experimental data to determine the self-consistent thermodynamic parameters of MgSiO3 akimotoite in a wide temperature and pressure range. All the obtained thermodynamic parameters, including thermal expansion, heat capacity, entropy, and Grüneisen parameter show a nonlinear and negative behavior with the increasing pressure. The pressure derivatives of the thermal expansion, heat capacity, and Grüneisen parameter generally increase with the temperature to ~450 K, ~1000 K, and ~650 K, respectively, and then decrease with the temperature. The pressure effect on entropy, on the other hand, generally increases with the temperature in the temperature range studied. This work provides a basic understanding of the thermodynamic properties of MgSiO3 akimotoite under high temperature and high-pressure conditions. Further work is needed to be done to obtain more accurate thermodynamic properties when considering the phase equilibria associated with MgSiO3 akimotoite.

Author Contributions

Conceptualization, Y.L.; methodology, C.S. and Y.L.; software, G.Y.; resources, C.S. and Q.L.; data curation, Q.L.; writing—Original draft preparation, C.S.; writing—Review and editing, W.S.; supervision, Y.L.; project administration, C.S.; funding acquisition, C.S. and Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Fundamental Research Funds for the Central Universities (Grant number: ZY20215133), the National Natural Science Foundation of China (Grant number: 42174225, 42404114), and the Natural Science Foundation of Hebei Province of China (Grant number: D2022512007).

Data Availability Statement

Data is contained within the article.

Acknowledgments

We thank Kexin Dong, Jingxian Fang, Xiaozhe Liu, and Yuanyuan Ye at the School of Geosciences, Institute of Disaster Prevention for helping the data analysis.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Irifune, T.; Ringwood, A.E. Phase transformations in a harzburgite composition to 26 GPa: Implications for dynamical behavior of the subducting slab. Earth Planet. Sci. Lett. 1987, 86, 365–376. [Google Scholar] [CrossRef]
  2. Ringwood, A.E.; Irifune, T. Nature of the 650–km seismic discontinuity: Implications for mantle dynamics and differentiation. Nature 1988, 331, 131–136. [Google Scholar] [CrossRef]
  3. Zhou, C.; Gréaux, S.; Nishiyama, N.; Irifune, T.; Higo, Y. Sound velocities measurement on MgSiO3 akimotoite at high pressures and high temperatures with simultaneous in situ X-ray diffraction and ultrasonic study. Phys. Earth Planet. Inter. 2014, 228, 97–105. [Google Scholar] [CrossRef]
  4. Siersch, N.C.; Kurnosov, A.; Criniti, G.; Ishii, T.; Boffa Ballaran, T.; Frost, D.J. The elastic properties and anisotropic behavior of MgSiO3 akimotoite at transition zone pressures. Phys. Earth Planet. Inter. 2021, 320, 106786. [Google Scholar] [CrossRef]
  5. Ishii, T.; Kojitani, H.; Akaogi, M. Phase relations of harzburgite and MORB up to the uppermost lower mantle conditions: Precise comparison with pyrolite by multisample cell high-pressure experiments with implication to dynamics of subducted slabs. J. Geophys. Res. Solid Earth 2019, 124, 3491–3507. [Google Scholar] [CrossRef]
  6. Gasparik, T. Phase relations in the transition zone. J. Geophys. Res. Solid Earth 1990, 95, 15751–15769. [Google Scholar] [CrossRef]
  7. Chanyshev, A.; Ishii, T.; Bondar, D.; Bhat, S.; Kim, E.J.; Farla, R.; Nishida, K.; Liu, Z.; Wang, L.; Nakajima, A.; et al. Depressed 660-km discontinuity caused by akimotoite–bridgmanite transition. Nature 2022, 601, 69–73. [Google Scholar] [CrossRef]
  8. Weidner, D.; Ito, E. Elasticity of MgSiO3 in the ilmenite phase. Phys. Earth Planet. Inter. 1985, 40, 65–70. [Google Scholar] [CrossRef]
  9. Reynard, B.; Fiquet, G.; Itie, J.-P.; Rubie, D.C. High-pressure X-ray diffraction study and equation of state of MgSiO3 ilmenite. Am. Mineral. 1996, 81, 45–50. [Google Scholar] [CrossRef]
  10. Wang, Y.; Uchida, T.; Zhang, J.; Rivers, M.L.; Sutton, S.R. Thermal equation of state of akimotoite MgSiO3 and effects of the akimotoite–garnet transformation on seismic structure near the 660 km discontinuity. Phys. Earth Planet. Inter. 2004, 143–144, 57–80. [Google Scholar] [CrossRef]
  11. Karki, B.; Wentzcovitch, R. First-principles lattice dynamics and thermoelasticity of MgSiO3 ilmenite at high pressure. J. Geophys. Res. 2002, 107, ECV 2-1–ECV 2-6. [Google Scholar] [CrossRef]
  12. Zhang, Y.; Zhao, D.; Matsui, M. Anisotropy of akimotoite: A molecular dynamics study. Phys. Earth Planet. Inter. 2005, 151, 309–319. [Google Scholar] [CrossRef]
  13. Li, L.; Weidner, D.J.; Brodholt, J.; Alfè, D.; Price, G.D. Ab initio molecular dynamics study of elasticity of akimotoite MgSiO3 at mantle conditions. Phys. Earth Planet. Inter. 2009, 173, 115–120. [Google Scholar] [CrossRef]
  14. Hao, S.; Wang, W.; Qian, W.; Wu, Z. Elasticity of akimotoite under the mantle conditions: Implications for multiple discontinuities and seismic anisotropies at the depth of ∼600–750 km in subduction zones. Earth Planet. Sci. Lett. 2019, 528, 115830. [Google Scholar] [CrossRef]
  15. Zhao, Y.; Wu, Z.; Hao, S.; Wang, W.; Deng, X.; Song, J. Elastic properties of Fe-bearing Akimotoite at mantle conditions: Implications for composition and temperature in lower mantle transition zone. Fundam. Res. 2022, 2, 570–577. [Google Scholar] [CrossRef]
  16. Ashida, T.; Kume, S.; Ito, E.; Navrotsky, A. MgSiO3 ilmenite: Heat capacity, thermal expansivity, and enthalpy of transformation. Phys. Chem. Miner. 1988, 16, 239–245. [Google Scholar] [CrossRef]
  17. Kojitani, H.; Yamazaki, M.; Tsunekawa, Y.; Katsuragi, S.; Noda, M.; Inoue, T.; Inaguma, Y.; Akaogi, M. Enthalpy, heat capacity and thermal expansivity measurements of MgSiO3 akimotoite: Reassessment of its self-consistent thermodynamic data set. Phys. Earth Planet. Inter. 2022, 333, 106937. [Google Scholar] [CrossRef]
  18. Akaogi, M.; Kojitani, H.; Morita, T.; Kawaji, H.; Atake, T. Low-temperature heat capacities, entropies and high-pressure phase relations of MgSiO3 ilmenite and perovskite. Phys. Chem. Miner. 2008, 35, 287–297. [Google Scholar] [CrossRef]
  19. Hernández, E.; Brodholt, J.; Alfè, D. Structural, vibrational and thermodynamic properties of Mg2SiO4 and MgSiO3 minerals from first-principles simulations. Phys. Earth Planet. Inter. 2014, 240, 1–24. [Google Scholar] [CrossRef]
  20. Davis, L.A.; Gordon, R.B. Compression of mercury at high pressure. J. Chem. Phys. 1967, 46, 2650–2660. [Google Scholar] [CrossRef]
  21. Su, C.; Liu, Y.; Fan, D.; Song, W.; Yang, G. Self-consistent thermodynamic parameters of pyrope and almandine at high-temperature and high-pressure conditions: Implication on the adiabatic temperature gradient. Phys. Earth Planet. Inter. 2022, 322, 106789. [Google Scholar] [CrossRef]
  22. Qiu, W.; Su, C.; Liu, Y.; Song, W. Thermodynamic properties of MgAl2O4 spinel at high temperatures and high pressures. Crystals 2023, 13, 240. [Google Scholar] [CrossRef]
  23. Johari, G.P. Entropy, enthalpy and volume of perfect crystals at limiting high pressure and the third law of thermodynamics. Thermochim. Acta 2021, 698, 178891. [Google Scholar] [CrossRef]
  24. Gonzalez-Platas, J.; Alvaro, M.; Nestola, F.; Angel, R. EosFit7-GUI: A new graphical user interface for equation of state calculations, analyses and teaching. J. Appl. Crystallogr. 2016, 49, 1377–1382. [Google Scholar] [CrossRef]
  25. Fei, Y. Thermal expansion. In Mineral Physics & Crystallography; AGU Reference Shelf; Wiley: Hoboken, NJ, USA, 1995; pp. 29–44. [Google Scholar]
  26. Maier, C.G.; Kelley, K.K. An equation for the representation of high-temperature heat content data. J. Am. Chem. Soc. 2002, 54, 3243–3246. [Google Scholar] [CrossRef]
  27. Haas, J.L.; Fisher, J.R. Simultaneous evaluation and correlation of thermodynamic data. Am. J. Sci. 1976, 276, 525–545. [Google Scholar] [CrossRef]
  28. Berman, R.G.; Brown, T.H. Heat capacity of minerals in the system Na2O-K2O-CaO-MgO-FeO-Fe2O3-Al2O3-SiO2-TiO2-H2O-CO2: Representation, estimation, and high temperature extrapolation. Contrib. Mineral. Petrol. 1985, 89, 168–183. [Google Scholar] [CrossRef]
  29. Richet, P.; Fiquet, G. High-temperature heat capacity and premelting of minerals in the system MgO-CaO-Al2O3-SiO2. J. Geophys. Res. 1991, 96, 445–456. [Google Scholar] [CrossRef]
  30. Hofmeister, A.M.; Ito, E. Thermodynamic properties of MgSiO3 ilmenite from vibrational spectra. Phys. Chem. Miner. 1992, 18, 423–432. [Google Scholar] [CrossRef]
  31. Saxena, S.K.; Chatterjee, N.; Fei, Y.; Shen, G. Thermodynamic Data on Oxides and Silicates, 1st ed.; Springer: Berlin/Heidelberg, Germany, 1993. [Google Scholar]
  32. Chopelas, A. Thermal expansivity of mantle relevant magnesium silicates derived from vibrational spectroscopy at high pressure. Am. Mineral. 2000, 85, 270–278. [Google Scholar] [CrossRef]
  33. Watanabe, H. Thermochemical Properties of Synthetic High-Pressure Compounds Relevant to the Earth’s Mantle; Springer: Dordrecht, The Netherlands, 1982; pp. 441–464. [Google Scholar]
  34. Benisek, A.; Dachs, E. The accuracy of standard enthalpies and entropies for phases of petrological interest derived from density-functional calculations. Contrib. Mineral. Petrol. 2018, 173, 90. [Google Scholar] [CrossRef] [PubMed]
  35. Vočadlo, L.; Poirer, J.P.; Price, G.D. Grüneisen parameters and isothermal equations of state. Am. Mineral. 2000, 85, 390–395. [Google Scholar] [CrossRef]
  36. Yusa, H.; Akaogi, M.; Ito, E. Calorimetric study of MgSiO3 garnet and pyroxene: Heat capacities, transition enthalpies, and equilibrium phase relations in MgSiO3 at high pressures and temperatures. J. Geophys. Res. Solid Earth 1993, 98, 6453–6460. [Google Scholar] [CrossRef]
  37. Yu, Y.G.; Wentzcovitch, R.M.; Vinograd, V.L.; Angel, R.J. Thermodynamic properties of MgSiO3 majorite and phase transitions near 660 km depth in MgSiO3 and Mg2SiO4: A first principles study. J. Geophys. Res. 2011, 116, B02208. [Google Scholar] [CrossRef]
  38. Wentzcovitch, R.M.; Stixrude, L.; Karki, B.B.; Kiefer, B. Akimotoite to perovskite phase transition in MgSiO3. Geophys. Res. Lett. 2004, 31, L10611. [Google Scholar] [CrossRef]
Figure 1. The unit-cell volume of MgSiO3 akimotoite as a function of temperature at ambient pressure conditions [16,17].
Figure 1. The unit-cell volume of MgSiO3 akimotoite as a function of temperature at ambient pressure conditions [16,17].
Crystals 14 00837 g001
Figure 2. The heat capacity of MgSiO3 akimotoite as a function of temperature at ambient pressure conditions [16,17,18].
Figure 2. The heat capacity of MgSiO3 akimotoite as a function of temperature at ambient pressure conditions [16,17,18].
Crystals 14 00837 g002
Figure 3. The unit-cell volumes of MgSiO3 akimotoite as a function of pressure at different temperatures [3,9,10,11].
Figure 3. The unit-cell volumes of MgSiO3 akimotoite as a function of pressure at different temperatures [3,9,10,11].
Crystals 14 00837 g003
Figure 4. The isothermal bulk modulus (a), adiabatic bulk modulus (b), and shear modulus (c) of MgSiO3 akimotoite as a function of pressure at different temperatures, and the bulk modulus (d), adiabatic bulk modulus (e), and shear modulus (f) of MgSiO3 akimotoite as a function of temperature at different pressures [3,11,14].
Figure 4. The isothermal bulk modulus (a), adiabatic bulk modulus (b), and shear modulus (c) of MgSiO3 akimotoite as a function of pressure at different temperatures, and the bulk modulus (d), adiabatic bulk modulus (e), and shear modulus (f) of MgSiO3 akimotoite as a function of temperature at different pressures [3,11,14].
Crystals 14 00837 g004aCrystals 14 00837 g004b
Figure 5. The thermal expansion of MgSiO3 akimotoite as a function of temperature at different pressures [11,17,31,32].
Figure 5. The thermal expansion of MgSiO3 akimotoite as a function of temperature at different pressures [11,17,31,32].
Crystals 14 00837 g005
Figure 6. The heat capacity of MgSiO3 akimotoite as a function of temperature at different pressures [11,16,17,31].
Figure 6. The heat capacity of MgSiO3 akimotoite as a function of temperature at different pressures [11,16,17,31].
Crystals 14 00837 g006
Figure 7. The entropy of MgSiO3 akimotoite as a function of temperature at different pressures [11,16,31].
Figure 7. The entropy of MgSiO3 akimotoite as a function of temperature at different pressures [11,16,31].
Crystals 14 00837 g007
Figure 8. The Grüneisen parameter of MgSiO3 akimotoite as a function of temperature at different pressures [11,31].
Figure 8. The Grüneisen parameter of MgSiO3 akimotoite as a function of temperature at different pressures [11,31].
Crystals 14 00837 g008
Figure 9. The thermal expansion (a), heat capacity (b), entropy (c), and Grüneisen parameter (d) of MgSiO3 akimotoite as a function of pressure at different temperatures.
Figure 9. The thermal expansion (a), heat capacity (b), entropy (c), and Grüneisen parameter (d) of MgSiO3 akimotoite as a function of pressure at different temperatures.
Crystals 14 00837 g009
Figure 10. The first and second pressure derivatives of thermal expansion (a), heat capacity (b), entropy (c), and Grüneisen parameter (d) of MgSiO3 akimotoite as a function of temperature.
Figure 10. The first and second pressure derivatives of thermal expansion (a), heat capacity (b), entropy (c), and Grüneisen parameter (d) of MgSiO3 akimotoite as a function of temperature.
Crystals 14 00837 g010aCrystals 14 00837 g010b
Table 1. Sound velocity, unit-cell volume, heat capacity and entropy data of MgSiO3 akimotoite used in this study.
Table 1. Sound velocity, unit-cell volume, heat capacity and entropy data of MgSiO3 akimotoite used in this study.
Data typeTemperature RangePressure RangeMethodReferences
KGPa
Sound velocity~1500~25.7Ultrasonic interferometryZhou et al. (2014) [3]
Unit-cell volume298–876AmbientX-ray diffractionAshida et al. (1988) [16]
298–773AmbientX-ray diffractionKojitani et al. (2022) [17]
Heat capacity170–700AmbientCalorimetryAshida et al. (1988) [16]
1.90–302.43AmbientThermal relaxation methodAkaogi et al. (2008) [18]
300–820AmbientCalorimetryKojitani et al. (2022) [17]
Standard entropy298.15AmbientCalorimetryAkaogi et al. (2008) [18]
Table 2. Fitting coefficients for the heat capacity of MgSiO3 akimotoite using different empirical equations.
Table 2. Fitting coefficients for the heat capacity of MgSiO3 akimotoite using different empirical equations.
Fitting EquationsReduced χ2R2References of the Empirical Equations
C P = c 0 + c 1 T + c 2 T 0.5 + c 3 T 1 + c 4 T 2 + c 5 T 3
c 0 = 247 433 × 10 2
c 1 = 8 9 × 10 2
c 2 = 2 2 × 10 4
c 3 = 2 2 × 10 5
c 4 = 1 2 × 10 7
c 5 = 3 9 × 10 8
1.190.996Kojitani et al. (2022) [17]
C P = c 0 + c 1 T + c 3 T 2 + c 4 T 0.5 + c 5 T 2
c 0 = 494 54
c 1 = 2.8 6 × 10 1
c 2 = 1.2 3 × 10 4
c 3 = 6.4 7 × 10 3
c 4 = 2.4 4 × 10 6
1.210.996Hass and Fisher (1976) [27]
C P = c 0 + c 1 T 0.5 + c 2 T 2 + c 3 T 3
c 0 = 146 4
c 1 = 6 1 × 10 2
c 2 = 5.3 4 × 10 6
c 3 = 6.1 5 × 10 8
1.180.996Berman and Brown (1985) [28]
C P = c 0 + c 1 l n T + c 2 T 1 + c 3 T 2 + c 4 T 3
c 0 = 155 133
c 1 = 5 39
c 2 = 1 2 × 10 4
c 3 = 4 3 × 10 6
c 4 = 5 2 × 10 7
1.190.996Richet and Fiquet (1991) [29]
Table 3. Isothermal bulk modulus, adiabatic bulk modulus, and shear modulus of MgSiO3 akimotoite at ambient conditions, and their pressure and temperature derivatives.
Table 3. Isothermal bulk modulus, adiabatic bulk modulus, and shear modulus of MgSiO3 akimotoite at ambient conditions, and their pressure and temperature derivatives.
TemperatureKTKSGReferences
KGPaGPaGPa
M0 (GPa)
300--212132Weidner and Ito (1985) [8]
298212----Reynard et al. (1996) [9]
300201----Karki & Wentzcovitch (2002) [11]
298210----Wang et al. (2004) [10]
300221226136Zhang et al. (2005) [12]
2000158.1 (6)--85.7Li et al. (2009) [13]
300207 (3)218.9 (6)131.8 (3)Zhou et al. (2014) [3]
300202204127Hao et al. (2019) [14]
300205 (1)209 (2)--Siersch et al. (2021) [4]
300221.5 (11)222.4 (13)130.8 (5)This study
M/∂P
3004.64----Karki & Wentzcovitch (2002) [11]
2985.6 (8)----Wang et al. (2004) [10]
3003.943.851.04Zhang et al. (2005) [12]
20003.7 (2)--4.5 (3)Li et al. (2009) [13]
3004.64.62 (3)1.64 (1)Zhou et al. (2014) [3]
3004.404.391.64Hao et al. (2019) [14]
300--4.4--Siersch et al. (2021) [4]
3004.49 (1)4.39 (1)1.54 (1)This study
M/∂T (×10−2 GPa/K)
300−2.5----Karki & Wentzcovitch (2002) [11]
298−4.0 (1)----Wang et al. (2004) [10]
300−3.0−2.16−1.79Zhang et al. (2005) [12]
300--−1.99(9)−1.58(4)Zhou et al. (2014) [3]
300−2.3−1.719−1.242Hao et al. (2019) [14]
300−2.943 (1)−1.937 (1)−1.645 (1)This study
2M/∂T2 (×10−6 GPa/K2)
300--−2.9 (8)−6.7 (4)Zhou et al. (2014) [3]
300--−1.26−1.94Hao et al. (2019) [14]
300−0.847 (6)−2.060 (8)−5.351 (4)This study
Table 4. Thermal expansion, heat capacity, entropy, and Grüneisen parameter of MgSiO3 akimotoite at ambient conditions, and their first and second pressure derivatives at different temperatures.
Table 4. Thermal expansion, heat capacity, entropy, and Grüneisen parameter of MgSiO3 akimotoite at ambient conditions, and their first and second pressure derivatives at different temperatures.
T
K
α
10−5/K
CP
J/mol∙K
S
J/mol∙K
γ
N0
10−5/KJ/mol∙KJ/mol∙K
3001.35 (4)77 (17)54(1)1.03 (25)
7001.44 (3)116 (9)138 (10)1.35 (2)
10003.03 (20)124 (8)181 (14)1.37 (2)
15003.37 (27)130 (7)233 (16)1.41 (7)
20003.66 (35)133 (6)271 (18)1.44 (12)
N/∂P
10−7/K∙GPaJ/mol∙K∙GPaJ/mol∙K∙GPa10−2/GPa
300−5.84 (2)−0.8066 (6)−0.315 (2)−2.200 (7)
700−5.79 (2)−0.2356 (6)−0.681 (2)−0.981 (3)
1000−6.64 (3)−0.2078 (11)−0.758 (2)−1.096 (4)
1500−8.01 (4)−0.2472 (20)−0.848 (3)−1.377 (5)
2000−9.42 (5)−0.3178 (28)−0.929 (4)−1.643 (6)
2N/∂P2
10−9/K∙GPa210−3J/mol∙K∙GPa210−3J/mol∙K∙GPa210−4/GPa2
3006.40 (6)0.53 (2)5.14 (6)−1.210 (2)
7006.29 (6)1.60 (2)5.99 (6)0.762 (10)
10007.13 (10)2.36 (3)6.69 (7)0.945 (13)
15009.12 (14)3.72 (6)7.90 (9)1.243 (18)
200011.19 (18)5.20 (9)9.17 (11)1.482 (21)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Su, C.; Song, W.; Yang, G.; Liu, Y.; Li, Q. Pressure Effects on the Thermodynamic Properties of MgSiO3 Akimotoite. Crystals 2024, 14, 837. https://doi.org/10.3390/cryst14100837

AMA Style

Su C, Song W, Yang G, Liu Y, Li Q. Pressure Effects on the Thermodynamic Properties of MgSiO3 Akimotoite. Crystals. 2024; 14(10):837. https://doi.org/10.3390/cryst14100837

Chicago/Turabian Style

Su, Chang, Wei Song, Guang Yang, Yonggang Liu, and Qingyi Li. 2024. "Pressure Effects on the Thermodynamic Properties of MgSiO3 Akimotoite" Crystals 14, no. 10: 837. https://doi.org/10.3390/cryst14100837

APA Style

Su, C., Song, W., Yang, G., Liu, Y., & Li, Q. (2024). Pressure Effects on the Thermodynamic Properties of MgSiO3 Akimotoite. Crystals, 14(10), 837. https://doi.org/10.3390/cryst14100837

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop