First-principles study on the high-pressure physical properties of orthocarbonate Ca 2 CO 4

doi:10.21203/rs.3.rs-2822097/v1

Download PDF

Article

First-principles study on the high-pressure physical properties of orthocarbonate Ca 2 CO 4

https://doi.org/10.21203/rs.3.rs-2822097/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 13 Jul, 2023

Read the published version in Scientific Reports →

You are reading this latest preprint version

Orthorhombic Ca₂CO₄(space group Pnma) is a recently discovered orthocarbonate whose high-pressure physical properties are critical for understanding the deep carbon cycle. Here, we study the structure, elastic and seismic properties of Ca₂CO₄-Pnma at 20–140 GPa using first-principles calculations, and compare them with the results of CaCO₃ polymorphs. The results show that the structural parameters of Ca₂CO₄-Pnma are in good agreement with the experimental results. It is the potential host of carbon in the Earth's mantle subduction slab, and its low wave velocity and small anisotropy may be the reason why it cannot be detected in seismic observation. The thermodynamic properties of Ca₂CO₄-Pnma at high temperature and high pressure are obtained using the quasi-harmonic approximation method. This study provides new insights into the behavior of carbonate and carbon cycles in the mantle.

Earth and environmental sciences/Solid earth sciences/Geophysics

Earth and environmental sciences/Solid earth sciences/Mineralogy

As the most important carbonate, CaCO₃ is transported to the deep mantle by subduction slab and plays a crucial role in the global long-term carbon cycle¹. It is also a mineral that plays a key role in biomineralization². However, CaCO₃ undergoes a series of phase transitions under high temperature and high pressure, forming various structures and polymorphs. So far, the predicted structures are mainly calcite, aragonite, post-aragonite, and pyroxene-like^3-9, and these structures and polymorphs have been experimentally verified^{3, 4, 6, 9-15}. Some studies also considered the reaction of calcium carbonate with MgO, SiO₂, and MgSiO₃^{7, 10, 16, 17}, while ignoring the reaction with CaO. Previously, Al-Shemali and Boldyrev¹⁸ mentioned the possible formation of calcium orthocarbonate Ca₂CO₄ in the CaCO₃+CaO system under high pressure. Recently, using AIRSS¹⁹ and USPEX²⁰ crystal structure prediction methods, Sagatova et al.²¹ discovered a new structure of calcium orthocarbonate Ca₂CO₄(space group Pnma) stable at 13-50 GPa and 2000 K, the carbon atoms in this phase are 4-fold coordinated, and the structure is similar to high temperature and high pressure α'_H-Ca₂SiO₄ phase²². Afterward, they found that Ca₂CO₄-Pnma was stable in the temperature range of 20–100 GPa and 1000–2000 K using the density functional theory within quasi-harmonic approximation²³. Subsequently, Binck et al.²⁴ verified the results of Sagatova et al.²³ with single-crystal diffraction experiments. In addition, other alkaline earth orthocarbonates, Mg₂CO₄-Pnma²⁵, Mg₂CO₄-P2₁/c²⁵, Sr₂CO₄-Pnma²⁶, and Ba₂CO₄-Pnma²⁶ have also been predicted, of which Sr₂CO₄-Pnma²⁷ and Mg₂CO₄-P2₁/c²⁸ have been experimentally verified.

The elastic, seismic, and thermodynamic properties of Ca₂CO₄-Pnma under high pressure have not been investigated so far. Even the elastic constants of CaCO₃ polymorphs were only the experimental results of calcite^29-32 and aragonite^{33, 34} at ambient conditions. Using the first-principles method, Belkofsi et al. calculated the elastic constants of three calcite polymorphs(calcite-III, calcite-IIIb, calcite-VI)³⁵, and Huang et al. studied the elastic properties of aragonite, post-aragonite and P2₁/c³⁶. The thermal expansion coefficient^37-42 and heat capacity^{37, 43-45} of calcite and aragonite were measured at ambient conditions, where there was a large difference between the fitted thermal expansion coefficient.

In this work, the structural properties, elastic properties, and seismic properties of Ca₂CO₄-Pnma at 20-140 GPa are studied using the first-principles calculations based on density functional theory and are compared with the results of CaCO₃ polymorphs. The thermodynamic properties of Ca₂CO₄-Pnma are obtained by quasi-harmonic approximation method.

First-principles calculations are done with using the VASP package^{46, 47} with projector-augmented wave⁴⁸. The exchange-correlation interactions adopt the Perdew-Burke-Ernzerhof functional within the generalized gradient approximation⁴⁹. The electronic configurations of the atoms are Ca: 3s²3p⁶4s², C: 2s²2p², O: 2s²2p⁴, respectively. The cutoff energy of the plane-wave basis is set to 900 eV. The k-point mesh generation and data processing are obtained by vaspkit program⁵⁰. The k-points mesh of Ca₂CO₄-Pnma, calcite, aragonite, P2₁/c-l, post-aragonite, P2₁/c-h and C222₁ are set to 5×7×4, 9×9×2, 7×4×6, 7×10×3, 8×7×8, 8×10×4, and 6×5×10 using the Monkhorst–Pack scheme⁵¹, respectively. The convergence criteria for energy and force are 1.0⨯10^–8 eV and 0.02 eV/Å, respectively. Based on the optimized lattice structure, the stress-strain method is used to obtain the elastic stiffness tensor. In order to ensure the accuracy of the elastic constants of Ca₂CO₄-Pnma, the elastic constants of calcite and aragonite are calculated and compared with the available experimental results^{32, 33}. As shown in Table S1(see Supplementary Material), the calculated results are in good agreement with the experimental results^{32, 33}. The thermodynamic properties are calculated using the quasi-harmonic approximation method⁵² of the PHONOPY program^{53, 54}, and the force constants are calculated using the density functional perturbation theory⁵⁵. The supercells of aragonite and Ca₂CO₄-Pnma adopt 2×2×2 and 2×2×1 unit cells, respectively.

Structural properties. Figure 1 shows the crystal structures of Ca₂CO₄-Pnma in the unit cell. The lattice parameters and equations of state for Ca₂CO₄-Pnma are presented in Fig. 2. It is found that the calculated results are in good agreement with the available experimental²⁴ and previous theoretical results^{21, 24}, indicating the validity of the structure. As the pressure increases, the lattice parameter c decreases the fastest, followed by a and b, indicating that the c-axis is the easiest to compress, while the a-axis is the most difficult to compress. The sensitivity of the axis to compression is c > b > a. The unit-cell volume at 0 GPa is 303.38 Å³ and the bulk modulus and its first pressure derivative are K₀ = 113.40 GPa and $K_{0}^{\prime }$ =4.00 by fitting the third-order Birch–Murnaghan equation, respectively, which are consistent with the results(V₀ = 302.0(3) Å³, K₀ = 108(1) GPa, and $K_{0}^{\prime }$= 4.43(3)) of Binck et al.²⁴.

In order to better understand the elastic and seismic properties of Ca₂CO₄-Pnma, the candidate CaCO₃ structures(aragonite, P2₁/c-l, post-aragonite, P2₁/c-h, C222₁, ’-l = low pressure’, ’-h = high pressure’) in the Earth's mantle are considered. The relative stabilities of the CaCO₃ polymorphs considered in this work are evaluated from their enthalpies. According to Fig. S1(see Supplementary Material), P2₁/c-l stabilizes above 30 GPa and retains its stability up to 46 GPa, while P2₁/c-h stabilizes above 75 GPa and retains its stability up to at least 140 GPa, which are consistent with the experimental and previous theoretical results^{3, 5}. CaCO₃-C222₁ above 137 GPa is stable relative to post-aragonite, but this does not make any sense^{5, 57}. Because in this interval, the modification P2₁/c-h is more favorable. For comparison with calcium orthocarbonate, four modifications of CaCO₃ must be considered, namely aragonite (20–35 GPa), P2₁/c-l (35–45 GPa), post-aragonite (45–75 GPa) and P2₁/c-h (75–140 GPa).

Elastic properties. The calculated elastic constants of Ca₂CO₄-Pnma are shown in Fig. 3. Within the studied pressure range, ${c_{11}}>{c_{22}}>{c_{33}}$, indicating that compression is easier along the c-axis than along the a- and b-axes. These results are consistent with those of Fig. 2, where the lattice parameter c decreases faster than the lattice parameters a and b with increasing pressure. The calculated elastic constants of CaCO₃ polymorphs are shown in Fig. 4, Fig. 5, Fig. 6, and Fig. 7, respectively. Therefore, we believe that the calculated elastic constants are correct, but experimental verification is required.

The bulk modulus (B) and shear modulus (G) of Ca₂CO₄-Pnma can be obtained by the Voigt⁵⁸-Reuss⁵⁹-Hill⁶⁰ scheme. As can be seen from Fig. 8, B is greater than G, indicating that with the change of volume, Ca₂CO₄-Pnma is more and more difficult to be compressed, and G is the main factor for the deformation of Ca₂CO₄-Pnma. The B and G of Ca₂CO₄-Pnma at < 75 GPa are larger than those of CaCO₃ polymorphs. The B of Ca₂CO₄-Pnma at 75–140 GPa is equal to that of P2₁/c-h, and the G is slightly larger and almost parallel.

In order to evaluate the elastic anisotropy of Ca₂CO₄-Pnma, we adopt the scheme of Ravindran et al.⁶¹. The shear anisotropic factors of A₁₀₀ in (100) plane, A₀₁₀ in (010) plane, and A₀₀₁ in (001) plane can be obtained from the following expression:

$${A_{100}}=\frac{{4{c_{44}}}}{{{c_{11}}+{c_{33}} - 2{c_{13}}}}$$

$${A_{010}}=\frac{{4{c_{55}}}}{{{c_{22}}+{c_{33}} - 2{c_{23}}}}$$

$${A_{001}}=\frac{{4{c_{66}}}}{{{c_{11}}+{c_{22}} - 2{c_{12}}}}$$

The variation of shear anisotropic factors A₁₀₀, A₀₁₀ and A₀₀₁ of Ca₂CO₄-Pnma with pressure is displayed in Fig. 9. A₀₁₀ and A₀₀₁ gradually decrease with increasing pressure, A₁₀₀ first increases with the increase of pressure, and then gradually decreases at > 40 GPa. It can also be found that the elastic anisotropy of Ca₂CO₄-Pnma in the lower mantle is very small, and the anisotropy of the (010) plane between [101] and [001] directions is the smallest.

Seismic properties. The compressional and shear wave velocities of minerals can be calculated from the elastic constants and densities. The compressional (V_P) and shear (V_S) wave velocities of Ca₂CO₄-Pnma and CaCO₃ polymorphs can be obtained from the Navier's equations⁶²:

$${V_P}=\sqrt {\frac{{3B+4G}}{{3\rho }}} ,\begin{array}{*{20}{c}} {}&{} \end{array}{V_S}=\sqrt {\frac{G}{\rho }}$$

The densities and wave velocities of Ca₂CO₄-Pnma, CaCO₃ polymorphs and the Preliminary Reference Earth Model (PREM)⁶³ are displayed in Fig. 10. From Fig. 10 (a), it is found that the densities of Ca₂CO₄-Pnma in the lower mantle is less than those of PREM, and greater than those of CaCO₃ polymorphs. As shown in Fig. 10(b), the V_P and V_S of CaCO₄-Pnma and CaCO₃ polymorphs are lower than those of PREM, and the V_P and V_S of Ca₂CO₄-Pnma are greater than those of P2₁/c-l and post-aragonite, which are almost the same as those of P2₁/c-h. The wave velocities in various crystallographic directions can be obtained by solving the Christoffel equation $\left| {{C_{ijkl}}{n_j}{n_l} - \rho {V^2}{\delta _{ik}}} \right|=0$⁶⁴. Figure 11 shows the wave velocities of Ca₂CO₄-Pnma along different crystallization directions at various pressures. The V_P of Ca₂CO₄-Pnma propagates the fastest in the [100] direction. The shear fast-wave velocity propagates the slowest in the [001] direction. With the increase of pressure, the propagation in the [100] and [010] directions become slower. The shear slow-wave velocity in [100] direction propagates more and more slowly as pressure increases.

The anisotropy A_P of the compressional waves and the polarization anisotropy A_S of the shear waves are defined as⁶⁶:

$${A_P}=\frac{{{V_{P,\hbox{max} }} - {V_{P,\hbox{min} }}}}{{{V_{P,aggregate}}}} \times 100\%$$

$${A_S}=\frac{{{{\left| {{V_{S1}} - {V_{S2}}} \right|}_{\hbox{max} }}}}{{{V_{S,aggregate}}}} \times 100\%$$

Figure 12 shows the A_P and A_S of Ca₂CO₄-Pnma and CaCO₃ polymorphs. It can be seen that the seismic anisotropy A_P and A_S of Ca₂CO₄-Pnma are less than those of CaCO₃ polymorphs, and decrease with the increase of pressure, and gradually increase at > 45 GPa.

The seismic properties of Ca₂CO₄-Pnma indicate that it is a potential host of carbon in the subduction slab and coexists with CaCO₃ polymorphs, as suggested by Sagatova et al.^{21, 23}. It was also verified by Binck et al.²⁴. The low wave velocity and small anisotropy of Ca₂CO₄-Pnma may be one of the reasons why it is impossible to detect the presence of carbonate in the lower mantle during the seismic observation of the subduction slab.

Thermodynamic properties. The thermodynamic parameters of minerals are a prerequisite for deriving the thermal state of the Earth's interior. In order to obtain the variation of thermodynamic parameters of Ca₂CO₄-Pnma with temperature and pressure, we first verify the constant pressure heat capacity ${C_P}$ of aragonite at 0 GPa, and find that the calculated results are in good agreement with the experimental results⁴⁴(Fig. 13). On this basis, the predicted heat capacity and thermal expansion coefficient $\alpha$ of Ca₂CO₄-Pnma are shown in Figs. 14 and 15, respectively.

On the basis of the determination of the stability for CaCO₃ polymorphs in the lower mantle and the verification of the structural parameters of Ca₂CO₄-Pnma, we study the elastic, seismic and thermodynamic properties of Ca₂CO₄-Pnma, and compared the results with those of CaCO₃ polymorphs. The research shows that the densities of Ca₂CO₄-Pnma in the lower mantle are greater than those of CaCO₃ polymorphs, and the seismic anisotropies are less than those of CaCO₃ polymorphs. The wave velocities of Ca₂CO₄-Pnma and CaCO₃ polymorphs are relatively low, and the wave velocities of Ca₂CO₄-Pnma and CaCO₃-P2₁/c-h are almost the same. This means that the presence of carbonate in the lower mantle is unlikely to be detected by seismic observations of subducted slab. By verifying the constant pressure heat capacity of aragonite at 0 GPa, the thermodynamic properties of Ca₂CO₄-Pnma at high temperature and high pressure are calculated using the quasi-harmonic approximation method. The results of this study are helpful to better understand the behavior of calcium carbonate in the lower mantle.

Data availability

The datasets used and/or analyzed during the current study available from the corresponding author on reasonable request.

Acknowledgements

This work is supported by the Industrial Support and Guidance Project of Colleges and Universities of Gansu Province (No. 2022CYZC-37), the Key Natural Science Foundation of Gansu Province (No. 20JR5RA211) and the Talent Innovation and Entrepreneurship Project of Lanzhou City (No. 2020-RC-18).

Author contributions

Z.-J.L. designed the calculations and wrote the manuscript. X.-W.S., C.-R.Z. and Z.-L.L analyzed the results. T.L., T.S. and X.-D.W. performed partial calculations. All authors reviewed the manuscript.

Competing interests

The authors declare no competing interests.

Additional information

Correspondence and requests for materials should be addressed to Z.-J.L.

Gao, J., Wu, X., Yuan, X. & Su, W. Fate of carbonates in the Earth’s mantle (10–136 GPa). Front. Earth Sci. 10, 837775 (2022).
Lowenstam, H. A. & Weiner, S., On Biomineralization. (Oxford University Press, New York, 1989).
Gavryushkin, P. N., Martirosyan, N. S., Inerbaev, T. M., Popov, Z. I., Rashchenko, S. V., Likhacheva, A. Y., Lobanov, S. S., Goncharov, A. F., Prakapenka, V. B. & Litasov, K. D. Aragonite-II and CaCO₃-VII: New high-pressure, high-temperature polymorphs of CaCO₃. Cryst. Growth Des. 17, 6291–6296 (2017).
Oganov, A. R., Glass, C. W. & Ono, S. High-pressure phases of CaCO₃: Crystal structure prediction and experiment. Earth Planet. Sci. Lett. 241, 95–103 (2006).
Pickard, C. J. & Needs, R. J. Structures and stability of calcium and magnesium carbonates at mantle pressures. Phys. Rev. B 91, 104101 (2015).
Smith, D., Lawler, K. V., Martinez-Canales, M., Daykin, A. W., Fussell, Z., Smith, G. A., Childs, C., Smith, J. S., Pickard, C. J. & Salamat, A. Postaragonite phases of CaCO₃ at lower mantle pressures. Phys. Rev. Mater. 2, 013605 (2018).
Zhang, Z. G., Mao, Z., Liu, X., Zhang, Y. G. & Brodholt, J. Stability and Reactions of CaCO₃ Polymorphs in the Earth's Deep Mantle. J. Geophys. Res. 123, 6491–6500 (2018).
Gavryushkin, P. N., Sagatov, N., Belonoshko, A. B., Banaev, M. V. & Litasov, K. D. Disordered aragonite: The new high-pressure, high-temperature phase of CaCO₃. J. Phys. Chem. C 124, 26467–26473 (2020).
Bayarjargal, L., Fruhner, C. J., Schrodt, N. & Winkler, B. CaCO₃ phase diagram studied with Raman spectroscopy at pressures up to 50 GPa and high temperatures and DFT modeling. Phys. Earth Planet. Inter. 281, 31–45 (2018).
Li, X., Zhang, Z., Lin, J. F., Ni, H., Prakapenka, V. B. & Mao, Z. New high-pressure phase of CaCO₃ at the topmost lower mantle: Implication for the deep‐mantle carbon transportation. Geophys. Res. Lett. 45, 1355–1360 (2018).
Ono, S., Kikegawa, T., Ohishi, Y. & Tsuchiya, J. Post-aragonite phase transformation in CaCO₃ at 40 GPa. Am. Mineral. 90, 667–671 (2005).
Lobanov, S. S., Dong, X., Martirosyan, N. S., Samtsevich, A. I., Stevanovic, V., Gavryushkin, P. N., Litasov, K. D., Greenberg, E., Prakapenka, V. B., Oganov, A. R. & Goncharov, A. F. Raman spectroscopy and x-ray diffraction of sp³ CaCO₃ at lower mantle pressures. Phys. Rev. B 96, 104101 (2017).
Ono, S., Kikegawa, T. & Ohishi, Y. High-pressure transition of CaCO₃. Am. Mineral. 92, 1246–1249 (2007).
Lv, M., Liu, J., Greenberg, E., Prakapenka, V. & Dorfman, S. Thermal equation of state of post-aragonite CaCO₃-Pmmn. Am. Mineral. 105, 1365–1374 (2020).
Druzhbin, D., Rashchenko, S. V., Shatskiy, A. & Crichton, W. A. New high-pressure and high-temperature CaCO₃ polymorph. ACS Earth Space Chem. 6, 1506–1513 (2022).
Seto, Y., Hamane, D., Nagai, T. & Fujino, F. Fate of carbonates within oceanic plates subducted to the lower mantle, and a possible mechanism of diamond formation. Phys. Chem. Miner. 35, 223–229 (2008).
Drewitt, J. W. E., Walter, M. J., Zhang, H., McMahon, S. C., Edwards, D., Heinen, B. J., Lord, O. T., Anzellini, S. & Kleppe, A. K. The fate of carbonate in oceanic crust subducted into earth's lower mantle. Earth Planet. Sci. Lett. 511, 213–222 (2019).
Al-Shemali, M. & Boldyrev, A. I. Search for ionic orthocarbonates: Ab initio study of Na₄CO₄. J. Phys. Chem. A 106, 8951–8954 (2002).
Pickard, C. J. & Needs, R. J. Ab initio random structure searching. J. Phys.: Condens Matter 23, 053201 (2011).
Oganov, A. R. & Glass, C. W. Crystal structure prediction using ab initio evolutionary techniques: principles and applications. J. Chem. Phys. 124, 244704 (2006).
Sagatova, D., Shatskiy, A., Sagatov, N., Gavryushkin, P. N. & Litasov, K. D. Calcium orthocarbonate, Ca₂CO₄-Pnma: A potential host for subducting carbon in the transition zone and lower mantle. Lithos 370–371, 105637 (2020).
Belmonte, D., Ottonello, G. & Zuccolini, M. V. Ab initio-assisted assessment of the CaO-SiO₂ system under pressure. Calphad 59, 12–30 (2017).
Sagatova, D. N., Shatskiy, A. F., Gavryushkin, P. N., Sagatov, N. E. & Litasov, K. D. Stability of Ca₂CO₄-Pnma against the main mantle minerals from ab initio computations. ACS Earth Space Chem. 5, 1709–1715 (2021).
Binck, J., Laniel, D., Bayarjargal, L., Khandarkhaeva, S., Fedotenko, T., Aslandukov, A., Glazyrin, K., Milman, V., Chariton, S., Prakapenka, V. B., Dubrovinskaia, N., Dubrovinsky, L. & Winkler, B. Synthesis of calcium orthocarbonate, Ca₂CO₄-Pnma at P-T conditions of Earth’s transition zone and lower mantle. Am. Mineral. 107, 336–342 (2022).
Gavryushkin, P. N., Sagatova, D. N., Sagatov, N. & Litasov, K. D. Formation of Mg-orthocarbonate through the reaction MgCO₃ + MgO = Mg₂CO₄ at Earth’s lower mantle P–T Conditions. Cryst. Growth Des. 21, 2986–2992 (2021).
Gavryushkin, P. N., Sagatova, D. N., Sagatov, N. & Litasov, K. D. Orthocarbonates of Ca, Sr, and Ba-the appearance of sp³-hybridized carbon at a low pressure of 5 GPa and dynamic stability at ambient pressure. ACS Earth Space Chem. 5, 1948–1957 (2021).
Spahr, D., Binck, J., Bayarjargal, L., Luchitskaia, R., Morgenroth, W., Comboni, D., Milman, V. & Winkler, B. Tetrahedrally coordinated sp³-hybridized carbon in Sr₂CO₄ orthocarbonate at ambient conditions. Inorg. Chem. 60, 5419–5422 (2021).
Gavryushkin, P. N., Martirosyan, N. S., Rashchenko, S. V., Sagatova, D. N., Sagatov, N. E., Semerikova, A. I., Fedotenko, T. M. & Litasov, K. D. First experimental synthesis of Mg orthocarbonate by the MgCO₃ + MgO = Mg₂CO₄ reaction at pressures of the Earth’s lower mantle. JETP Lett. 116, 477–484 (2022).
Chen, C.-C., Lin, C.-C., Liu, L.-G., Sinogeikin, S. V. & Bass, J. D. Elasticity of single-crystal calcite and rhodochrosite by Brillouin spectroscopy. Am. Mineral. 86, 1525–1529 (2001).
Dandekar, D. P. Variation in the elastic constants of calcite with temperature. J. Appl. Phys. 39, 3694–3699 (1968).
Dandekar, D. P. Pressure dependence of the elastic constants of calcite. Phys. Rev. 172, 873–877 (1968).
Hearmon, R. F. S., The elastic constants of crystals and other anisotropic materials. K. H. Hellwege, A. M. Hellwege, Eds., Landolt-Bornstein Tables,III/11 (Springer-Verlag, Berlin, 1979).
Liu, L.-g., Chen, C.-c., Lin, C.-C. & Yang, Y.-j. Elasticity of single-crystal aragonite by Brillouin spectroscopy. Phys. Chem. Minerals 32, 97–102 (2005).
Hearmon, R. F. S. The elastic constants of anisotropic minerals. Rev. Mod. Phys. 8, 409–440 (1946).
Belkofsi, R., Adjaoud, O. & Belabbas, I. Pressure induced phase transitions and elastic properties of CaCO₃ polymorphs: a density functional theory study. Modelling Simul. Mater. Sci. Eng. 26, 065004 (2018).
Huang, D., Liu, H., Hou, M.-Q., Xie, M.-Y., Lu, Y.-F., Liu, L., Yi, L., Cui, Y.-J., Li, Y., Deng, L.-W. & Du, J.-G. Elastic properties of CaCO₃ high pressure phases from first principles. Chin. Phys. B 26, 089101 (2017).
Litasov, K. D., Shatskiy, A., Gavryushkin, P. N., Bekhtenova, A. E., Dorogokupets, P. I., Danilov, B. S., Higo, Y., Akilbekov, A. T. & Inerbaev, T. M. P-V-T equation of state of CaCO₃ aragonite to 29 GPa and 1673 K: In situ X-ray diffraction study. Phys. Earth Planet. Inter. 265, 82–91 (2017).
Palaich, S. E. M., Heffern, R. A., Hanfland, M., Lausi, A., Kavner, A., Manning, C. E. & Merlini, M. High-pressure compressibility and thermal expansion of aragonite. Am. Mineral. 101, 1651–1658 (2016).
Li, Y., Zou, Y., Chen, T., Wang, X., Qi, X., Chen, H., Du, J. & Li, B. P-V-T equation of state and high-pressure behavior of CaCO₃ aragonite. Am. Mineral. 100, 2323–2329 (2015).
Antao, S. M. & Hassan, I. Temperature dependence of the structural parameters in the transformation of aragonite to calcite, as determined from in situ synchrotron powder X-ray-diffraction data. Can. Mineral. 48, 1225–1236 (2010).
Matas, J., Gillet, P., Ricard, Y. & Martinez, I. Thermodynamic properties of carbonates at high pressures from vibrational modelling. Eur. J. Mineral. 12, 703–720 (2000).
Martinez, I., Zhang, J. & Reeder, R. J. In situ X-ray diffraction of aragonite and dolomite at high pressure and high temperature: Evidence for dolomite breakdown to aragonite and magnesit. Am. Mineral. 81, 611–624 (1996).
Lucas, A., Mouallem-Bahout, M., Carel, C., Gaude, J. & Matecki, M. Thermal expansion of synthetic aragonite condensed review of elastic properties. J. Solid State Chem. 146, 73–78 (1999).
Staveley, L. A. K. & Linford, R. G. The heat capacity and entropy of calcite and aragonite, and their interpretation. J. Chem. Thermodynamics 1, 1–11 (1969).
Robie, R. A., Hemingway, B. S. & Fisher, J. R., Thermodynamic properties of minerals and related substances at 298.15 K and 1 bar (105 pascals) pressure and at higher temperatures. Bulletin (U.S. G.P.O., Washington, 1978).
Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54, 11169–11186 (1996).
Kresse, G. & Furthmüller, J. Efficiency of ab-initio total enery calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 6, 15–50 (1996).
Blöchl, P. E. Projector augmented-wave method. Phys. Rev. B 50, 17953–17979 (1994).
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865–3868 (1996).
Wang, V., Xu, N., Liu, J. C., Tang, G. & Geng, W. T. VASPKIT: A user-friendly interface facilitating high-throughput computing and analysis using VASP code. Comput. Phys. Commun. 267, 108033 (2021).
Monkhorst, H. J. & Pack, J. D. Special points for Brillouin-zone integrations. Phys. Rev. B 13, 5188–5192 (1976).
Togo, A., Chaput, L., Tanaka, I. & Hug, G. First-principles phonon calculations of thermal expansion in Ti₃SiC₂, Ti₃AlC₂, and Ti₃GeC₂. Phys. Rev. B 81, 174301 (2010).
Togo, A. & Tanaka, I. First principles phonon calculations in materials science. Scr. Mater. 108, 1–5 (2015).
Togo, A. First-principles phonon calculations with phonopy and phono3py. J. Phys. Soc. Jpn. 92, 012001 (2023).
Gonze, X. & Lee, C. Dynamical matrices, Born effective charges, dielectric permittivity tensors, and interatomic force constants from density-functional perturbation theory. Phys. Rev. B 55, 10355–10368 (1997).
Momma, K. & Izumi, F. VESTA 3 for three-dimensional visualization of crystal, volumetric and morphology data. J. Appl. Crystallogr. 44, 1272–1276 (2011).
Yao, X., Xie, C., Dong, X., Oganov, A. R. & Zeng, Q. Novel high-pressure calcium carbonates. Phys. Rev. B 98, 014108 (2018).
Voigt, W., Lehrbuch der Kristallphysik-mit Ausschluß der Kristalloptik. (Vieweg + Teubner Verlag, Wiesbaden, 1966).
Reuss, A. Calulation of yielding mixed crystals plasticity condition for single crystals. Z. Angew. Math. Mech. 9, 49–58 (1929).
Hill, R. The elastic behavior of a crystalline aggregate. Proc. Phys. Soc., London 65A, 349–354 (1952).
Ravindran, P., Fast, L., Korzhavyi, P. A., Johansson, B., Wills, J. & Eriksson, O. Density functional theory for calculation of elastic properties of orthorhombic crystals: Application to TiSi₂. J. Appl. Phys. 84, 4891–4904 (1998).
Anderson, O. L. A simplified method for calculating the debye temperature from elastic constants. J. Phys. Chem. Solids 24, 909–917 (1963).
Dziewonski, A. M. & Anderson, D. L. Preliminary reference Earth model. Phys. Earth Planet. Inter. 25, 297–356 (1981).
Every, A. G. General closed-form expressions for acoustic waves in elastically anisotropic solids. Phys. Rev. B 22, 1746–1760 (1980).
Muñoz-Santiburcio, D. & Hernández-Laguna, A. AWESoMe 1.1: A code for the calculation of phase and group velocities of acoustic waves in homogeneous solids. Comput. Phys. Commun. 217, 212–214 (2017).
Kiefer, B., Stixrude, L., Hafner, J. & Kresse, G. Structure and elasticity of wadsleyite at high pressures. Am. Mineral. 86, 1387–1395 (2001).

No competing interests reported.

SupplementaryMaterials.docx

Download PDF

Journal Publication

published 13 Jul, 2023

Read the published version in Scientific Reports →

Editorial decision: Major revision
06 Jun, 2023
Reviews received at journal
05 Jun, 2023
Reviews received at journal
17 May, 2023
Reviewers agreed at journal
11 May, 2023
Reviewers invited by journal
11 May, 2023
Editor assigned by journal
09 May, 2023
Editor invited by journal
21 Apr, 2023
Submission checks completed at journal
21 Apr, 2023
First submitted to journal
15 Apr, 2023

You are reading this latest preprint version

First-principles study on the high-pressure physical properties of orthocarbonate Ca 2 CO 4

Status:

Journal Publication

Version 1

Abstract

Figures

Introduction

Methods

Results And Discussion

Conclusions

Declarations

References

Additional Declarations

Supplementary Files

Status:

Journal Publication

Version 1