Structural transitions of calcium carbonate by molecular dynamics simulation

Elizaveta Sidler,Raffaela Cabriolu

2024-08-08

Abstract:Calcium carbonate plays a crucial role in the global carbon cycle, and its phase diagram has always been of significant scientific interest. In this study, we used molecular dynamics (MD) to investigate several structural phase transitions of calcium carbonate. Using the Raiteri potential model, we explored the structural transitions occurring at a constant pressure of 1 bar with temperatures ranging from 300 K to 2500 K, and at a constant temperature of 1600 K with pressures ranging from 0 to 13 GPa. At increasing temperatures, the transitions calcite, CaCO$_3$-IV, and CaCO$_3$-V are observed and characterized. Within the calcite structure, CO$_3^{2-}$ ions are ordered between layers. As temperature increases, the calcite to CaCO$_3$-IV transition occurs, determining the partial disordering of CO$_3^{2-}$ ions. At a higher temperature, CaCO$_3$-IV transforms into CaCO$_3$-V. By applying free energy analysis, we have classified the last transition as a continuous order-disorder phase transition. At a temperature of 2000 K, it appears a `disordered CaCO$_3$' structure, characterized by low order within the calcium and carbonate sublattices and the free rotation of CO$_3^{2-}$ ions. At increasing pressures, two calcium carbonate transformations were observed. At $P=$ 2 GPa, the CaCO$_3$-IV phase undergoes a phase transition into CaCO$_3$-V, demonstrating that the model can describe the transition between these two phases as pressure and temperature-driven. Another phase transition was detected at $P=$ 4.25 GPa -- CaCO$_3$-V transits into the recently discovered CaCO$_3$-Vb phase. This transition is classified as a first-order phase transition by structural analysis and free energy-based arguments.

Statistical Mechanics,Computational Physics

What problem does this paper attempt to address?

The paper primarily focuses on studying the structural phase transitions of calcium carbonate (CaCO₃) under different temperature and pressure conditions, exploring these phase transition processes and their mechanisms through molecular dynamics simulations. Below is a summary of the key issues the paper attempts to address:

### Issues the Paper Attempts to Address

1. **Investigate the structural phase transitions of calcium carbonate under different temperatures and pressures**:
- Use molecular dynamics (MD) simulation methods to study several structural phase transitions of calcium carbonate.
- Explore phase transitions within the temperature range from 300 K to 2500 K at constant pressure (1 bar) and within the pressure range from 0 to 13 GPa at constant temperature (1600 K).

2. **Describe specific phase transition processes**:
- Observe and characterize the phase transitions from Calcite to CaCO₃-IV to CaCO₃-V as the temperature increases.
- At 2000 K, a "disordered CaCO₃" structure appears, where the sublattice order of calcium and carbonate ions is low, and the carbonate ions rotate freely.

3. **Analyze phase transitions under different pressures**:
- At 2 GPa, the CaCO₃-IV phase transitions to the CaCO₃-V phase.
- At 4.25 GPa, the CaCO₃-V phase transitions to the recently discovered CaCO₃-Vb phase.

4. **Classify phase transition types using free energy analysis**:
- Classify certain phase transitions as second-order order-disorder transitions through free energy analysis.
- Classify other phase transitions as first-order transitions.

5. **Compare simulation results with experimental data**:
- Compare thermodynamic properties such as heat capacity, enthalpy, volume changes, and linear thermal expansion coefficients obtained from simulations with experimental data.

6. **Structural analysis and phase identification**:
- Use tools such as radial distribution function (RDF), Steinhardt's bond order parameters, and synthetic X-ray diffraction (XRD) patterns for structural analysis.
- Determine the unit cells of different phases by analyzing XRD data.

7. **Defect analysis**:
- Use the centrosymmetry parameter (CSP) to analyze crystal symmetry and locate potential defects.

8. **Study of free energy landscapes**:
- Use metadynamics methods to calculate free energy to gain deeper insights into the phase transition processes driven by different temperatures and pressures.

Through the above research, the paper aims to deepen our understanding of the complex phase diagram of calcium carbonate, especially the phase transition behavior under high temperature and high pressure conditions, which is significant for understanding geochemical cycles, geological formation processes, and related industrial applications.