Raffaella
Demichelis
*a,
Paolo
Raiteri
a,
Julian D.
Gale
a and
Roberto
Dovesi
b
aNanochemistry Research Institute, Department of Chemistry, Curtin University, GPO Box U1987, Perth, WA 6845, Australia. E-mail: raffaella.demichelis@curtin.edu.au
bDipartimento di Chimica IFM, Università degli Studi di Torino and NIS -Nanostructured Interfaces and Surfaces - Centre of Excellence, Via Giuria 7, 10125, Torino, Italy
First published on 27th October 2011
Both of the previously proposed Pbnm and P6522 ordered structures for vaterite are found to be unstable transition states using first-principles methods. Five stable structures are located, the lowest energy one being of P3221 symmetry. Since interconversion between these structures requires only thermal energy, this provides an additional source of disorder within the vaterite structure.
Whereas the crystal structures of calcite and aragonite are ordered and well defined, that of vaterite remains a matter of debate, in spite of many experimental and theoretical studies carried out for more than half a century (see Table 1 for a concise summary). Generally, different interpretations of X-ray diffraction data collected by various authors have suggested two candidate crystal systems: orthorhombic (Pbnm) and hexagonal (P63/mmc, P63/mmc with partial occupancy of CO32− sites, or P6322). On the basis of Raman spectroscopic studies, Behrens et al.5 concluded that the space group of vaterite can be neither P63/mmc, nor P6322, whereas Gabrielli et al.6 agreed with Meyer,7 proposing a P63/mmc structure with partial occupancy of CO32− sites, though with some residual ambiguities. Recent solid-state NMR data also supports the hexagonal structural model in preference to the orthorhombic one.8 On the contrary, Le Bail et al.9 propose that the structure is indeed ordered and orthorhombic, but with three misaligned domains separated by 120°, thereby leading to an apparent hexagonal character. Based on this interpretation of their X-ray data, they proposed the Ama2 space group, while noting the presence of some unexplained superstructure reflections. Analysis of the vaterite Raman spectrum by Wehrmeister et al.10 shows that none of the current proposed structures are consistent with their data, which indicates the presence of three inequivalent carbonate groups.
Reference | SG | a | b | c |
---|---|---|---|---|
a Partial occupancy. b First-principles calculation. | ||||
Meyer11 | Pbnm | 4.13 | 7.15 | 8.48 |
McConnell12 | P6 3 22 | 7.135 | 7.135 | 8.524 |
Kahmi,13 Sato and Matsuda14 | P6 3 /mmc | 4.13 | 4.13 | 8.49 |
Bradley et al.15 | P6 3 22 | 7.135 | 7.135 | 8.524 |
Meyer,7 Gabrielli et al.6 | P6 3 /mmca | 7.15 | 7.15 | 16.96 |
Dupont et al.16 | P6 3 /mmca | 7.169 | 7.169 | 16.98 |
Le Bail et al.9 | Ama2 | 8.7422 | 7.1576 | 4.1265 |
Medeiros et al.17 | Pbnm b | 4.531 | 6.640 | 8.480 |
Wang and Becker18 | P6 5 22 b | 7.290 | 7.290 | 25.302 |
In the last few years, Medeiros et al.17 and Wang and Becker18 have examined this problem by using first-principles calculations (pseudopotential plane-wave methodology, employing Density Functional Theory, DFT, at the level of the Generalized Gradient Approximation, GGA). In particular, the former study analysed the structural, electronic and optical properties of Pbnm-structured vaterite, whereas the latter investigators proposed the P6522 space group as a result of an energetic analysis of different hexagonal configurations. However, neither of these studies computed the vibrational spectra for the candidate structures. Not only can this be valuable for comparison to experimental data, but also turns out to be very revealing in terms of the stability of the proposed space groups, as will be demonstrated later.
Here we re-examine all of the possible ordered structures of vaterite and thereby arrive at a new proposed model. The present calculations were performed with the quantum-mechanical ab initioCRYSTAL09 program.19 All-electron Gaussian-type basis sets optimised for calcite20 were adopted, i.e. 86–511d(21)G for Ca, 6–311d(11)G for C and 8–411d(11)G for O. Pbnm and P6522-structured vaterite were computed first at the same level of theory as the previous first-principles calculations,17,18,21 obtaining very similar structural results. The structural and energetic analyses were then performed with the PBEsol22 functional, which is a revised version of the widely used PBE functional specifically tailored for solids. This functional has been shown to provide accurate results for H-free systems and, in general, correct relative stabilities.23,24 The corresponding hybrid functionals, PBE025 and PBEsol0, containing 25% of exact Hartree–Fock (HF) exchange, and B3LYP,26 were also used to validate the present results.
The DFT exchange–correlation contribution was evaluated by numerical integration over the unit cell volume, using a pruned grid with 75 radial (Gauss–Legendre radial quadrature) and 974 angular (Lebedev two-dimensional generation) points. Thresholds for the Coulomb and HF exchange series accuracy were set to 10−8, 10−8, 10−8, 10−8 and 10−16 for the exact, and to 10−18 for the bipolar-approximated, integral calculation.19 Reciprocal space was sampled according to a Monkhorst–Pack mesh with shrinking factors of 8 and 4 for the orthorhombic and hexagonal (including subgroups) unit cells, respectively.
Structure optimisations were performed using analytical gradients with respect to atomic coordinates and unit-cell parameters, within a quasi-Newtonian scheme combined with Broyden-Fletcher-Goldfarb-Shanno27 Hessian updating. The default convergence criteria were used for both gradient components and nuclear displacements. Tightening these thresholds by an order of magnitude led to differences of the order of only 10−3 Å in lattice parameters and 10−2 kJ mol−1 in the electronic energy.
The harmonic vibrational frequencies were computed by diagonalising the dynamical matrix, obtained by central finite differencing of the analytic gradients with respect to atomic Cartesian coordinates. Scanning of geometries along imaginary normal modes was performed automatically by use of the SCANMODE option19 and the FINDSYM software was applied for symmetry analysis.28 Tolerances on the Self Consistent Field were set to 10−8 and 10−10 a.u. for optimisation and frequency calculation, respectively.
The calculated structures (PBEsol) are shown in Table 2 and their relative energies in Table 3. As shown by Wang and Becker,18 the P6522-structure for vaterite is more stable than the Pbnm-form. However, while this previous work finds a stability difference of 71.7 kJ mol−1, our results give values of the order of 2 kJ mol−1 with all the functionals examined. Given that the experimental enthalpy differences between the extremes of amorphous calcium carbonate and calcite only span ∼15 kJ mol−1,29 we believe that the present energetics are more reasonable.
P6 5 22 (TS) | P6 5 | P3 2 21 | Pbnm (TS) | C2cm b | P212121 | P112 1 c | |
---|---|---|---|---|---|---|---|
a Lattice parameters [Å], volume per formula unit [Å3], density [g cm−3], minimum, maximum and average bond distances [Å], Ca coordination number (CN). Structural parameters obtained with the other functionals differ by less than 2%. The experimental lattice parameters are given in Table 1. TS implies that a structure is a “transition state”. b Ama2 is permuted into C2cm for ease of comparison to the other orthorhombic structures. c γ = 60.34°. | |||||||
a | 7.1206 | 7.1121 | 7.1239 | 4.3907 | 4.5026 | 4.3668 | 7.1115 |
b | 7.1206 | 7.1121 | 7.1239 | 6.6091 | 6.3905 | 6.5831 | 7.1013 |
c | 25.4390 | 25.4089 | 25.3203 | 8.5029 | 8.4905 | 8.4282 | 25.3601 |
Volume | 62.06 | 61.83 | 61.82 | 61.69 | 61.08 | 60.57 | 61.83 |
Density | 2.674 | 2.684 | 2.684 | 2.691 | 2.71 | 2.740 | 2.684 |
Ca–Omin | 2.304 | 2.284 | 2.280 | 2.346 | 2.334 | 2.344 | 2.269 |
Ca–Omax | 2.423 | 2.687 | 2.622 | 2.367 | 2.422 | 2.603 | 2.625 |
Ca–Oav | 2.359 | 2.384 | 2.392 | 2.358 | 2.375 | 2.407 | 2.385 |
CN(Ca) | 6 | 6 | 6,7 | 6 | 6 | 7 | 6,7,8 |
C–Omin | 1.286 | 1.281 | 1.280 | 1.283 | 1.288 | 1.282 | 1.281 |
C–Omax | 1.297 | 1.306 | 1.306 | 1.296 | 1.295 | 1.309 | 1.307 |
C–Oav | 1.293 | 1.292 | 1.293 | 1.292 | 1.292 | 1.293 | 1.292 |
PBE | PBEsol | PBE0 | PBEsol0 | B3LYP | |
---|---|---|---|---|---|
a TS implies that the structure is a “transition state”. | |||||
P6 5 22 (TS) | −2.05 | −2.04 | −2.36 | −2.31 | — |
P6 5 | −2.07 | −2.58 | — | −2.79 | — |
P3 2 21 | −2.94 | −3.50 | −3.34 | −3.63 | −2.68 |
Pbnm (TS) | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
Ama2 | +16.13 | +15.09 | — | +16.42 | — |
s4 (TS) | −2.35 | −2.80 | — | −3.01 | — |
P212121 | −0.97 | −2.27 | −1.22 | −2.11 | −0.11 |
P112 1 | −2.43 | −3.01 | — | −3.22 | — |
More significantly, both of the previously theoretically examined ordered structures for vaterite turn out to possess imaginary phonons at the Γ point and therefore are unstable transition states, rather than genuine minima. To be precise, four imaginary vibrational frequencies are found for the P6522 structure, in addition to the three translational modes, while the Pbnm configuration possesses one unstable mode. The hexagonal form exhibits imaginary modes at 63i, 47i, 47i and 36i cm−1 (PBEsol, with similar results for PBE), while the unstable mode of the latter falls between 28i and 50i cm−1 (depending on the adopted functional). A small, real frequency (17 cm−1) is obtained as the lowest vibrational mode for Pbnm-vaterite with the B3LYP functional, probably as a result of the larger predicted volume. Despite the marginal stability for this form, this structure remains the highest energy configuration with the B3LYP functional.
The above raises the question of why the instabilities were not previously observed in first-principles studies of vaterite? Aside from the fact that the phonon spectrum has not been previously computed, it might have been expected that the structures could naturally distort to lower symmetries due to the inherent numerical noise in the forces present in first-principles plane-wave calculations. However, both prior studies employed codes that allow the space group symmetry to be constrained by default. Furthermore, the results of such calculations with regard to symmetry breaking have been shown to be sensitive to the details of the pseudopotential employed, whereas the present work utilizes all-electron methods.30 To confirm our findings, we have computed the vibrational frequencies for the P6522 structure of Wang and Becker18 using the same computational approach as given in their paper. This also yields several imaginary modes demonstrating that the structure is only stable due to the space group being automatically constrained.
Having determined that both of the widely considered ordered structures for vaterite are unstable, the next step is to determine the energy minima that are connected via these transition states. First, we performed a scan along the imaginary mode of the Pbnm form, which is symmetry forbidden by the inversion centre. By reducing the point symmetry from D2h to D2, the minimum energy structure was located and determined to belong to the space groupP212121, and is more stable than Pbnm-vaterite by between 0.1 and 2.3 kJ mol−1, depending on the adopted functional. As shown in Table 2, the two orthorhombic structures are very similar, with the P212121 one being denser and allowing Ca and C to assume a more distorted configuration, as shown by the bond length distribution and coordination number. The alternative orthorhombic model (Ama2), recently proposed by Le Bail et al.,9 is indeed a distinct stable energy minimum. However, it is found to be ∼16 kJ mol−1 less stable than Pbnm-vaterite, and so can be discounted as a viable structure for this material.
For the case of the hexagonal P6522-structure the situation is more complicated as there are four different imaginary modes leading to distinct minima. As a result of the scans along these four modes (all corresponding to symmetry-forbidden rotations of CO32−groups), three minimum energy structures, with space groups of P1121, P65 and P3221, were obtained. A fourth structure (labelled in the following discussion as s4) turned out to correspond to a 2nd-order saddle point and so was not further analysed in detail.
P3 2 21-vaterite, shown in Fig. 1, is found to be the lowest energy structure (the unit cells and coordinates of all of the minima are reported in the ESI†). Differences between P6522-vaterite and the structures belonging to its subgroups are quite small, mainly corresponding to more distorted positions of O atoms, whereas Ca and C atoms largely remain on the equivalent higher symmetry sites. This is shown by the varying coordination of Ca atoms and the bond length distribution.
![]() | ||
Fig. 1 Graphical representation of the hexagonal vaterite structures: view of the 18 Ca atoms in the unit cell, with the corresponding coordination polyhedra, projected on to the ac plane. 6-, 7-, and 8-fold coordinated Ca atoms and polyhedra are shown in different shades of green; C and O atoms are grey and red, respectively. |
According to the Ca–O distances in the other CaCO3 polymorphs, we can consider O being part of the first-neighbour shell for Ca when this distance is less than 2.65 Å. As shown in Table 2, the largest Ca–O bond distance within the first-neighbour shell is <2.5 Å when vaterite is optimised in the P6522, Pbnm and Ama2 space groups, the second-neighbour shell being at Ca–O > 3 Å. In our new stable structures there are Ca–O bonds that fall over a broader range, which is qualitatively more consistent with experimental EXAFS data for vaterite.1
The computed IR and Raman spectra of the lowest energy structure are in good agreement with the experimental findings with regard to frequencies,5,6,10,14 though a quantitative comparison of the peak intensities was not possible. Lattice modes and CO3 rigid rotations fall below 430 cm−1, in-plane and out-of-plane CO3 bending modes are found in the 640–740 and 830–840 cm−1 ranges, respectively. Symmetric and asymmetric C–O stretching modes lie in the 1060–1080 and 1380–1550 cm−1 ranges, respectively. However, no major difference was noticed with respect to the vibrational spectra computed for the other vaterite models. Hence this only confirms the validity of the new structure, rather than providing a means of discriminating between candidates.
The vibrational thermodynamic functions, i.e. zero point energy, vibrational entropy and heat capacity contributions, were computed at 298 K for the minimum energy structures with the PBEsol functional. An 80-atom supercell was built for P212121-vaterite in order to allow for phonon dispersion, the other cells being large enough already for the purposes of this comparison. These results predict that P3221-vaterite is more stable than P65-, P212121- and P1121-vaterite by 1.0, 2.6 and 0.5 kJ mol−1, respectively, in terms of free energy, whereas the corresponding enthalpy differences are 0.4, 1.1 and 0.6 kJ mol−1 in favour of P3221-vaterite.
The variation in the relative stability of the various structures with the different functionals considered requires further comment. As already noted in previous works,23 the relative stabilities of phases having different densities are sensitive to the choice of functional due to systematic errors. In particular, the tendency of PBE, B3LYP and PBE0 to over-stabilise less dense structures has been documented. Given this, and the small magnitude of the energy differences involved, it is important to calibrate the quality of the present results. To do this we have computed the thermodynamics of vaterite relative to calcite, as given in Table 4 based both on the lowest energy structure and for the Boltzmann-weighted average of the three distorted structures (i.e. the three minima with P3221, P1121 and P65 symmetries are averaged with weights of 0.40, 0.33 and 0.27, respectively, computed based on their relative free energies). This demonstrates that this free energy difference is in good agreement with experiment for both PBEsol and PBEsol0. To put this in perspective, the free energy difference for the previous ordered literature structures would have been in error by a factor of a third or greater. For further verification, the thermodynamics of aragonite relative to calcite computed at the PBEsol0 level gives ΔH298.15, ΔS298.15 and ΔG298.15 of 0.8 (0.2) kJ mol−1, 3.5 (3.7) J K−1 mol−1 and −0.2 (−0.9) kJ mol−1 (experimental values in parenthesis).31 Regardless of any systematic errors, P3221-vaterite turns out to be the most stable ordered structure with all the functionals considered. Furthermore, although the energy differences are small, the likely accuracy of the quantum mechanical methods is higher than for many chemical reactions because of the close similarity of the states being compared.
Footnote |
† Electronic supplementary information (ESI) available: Fractional coordinates of the minimum energy vaterite structures as obtained with the PBEsol functional. See DOI: 10.1039/c1ce05976a |
This journal is © The Royal Society of Chemistry 2012 |