On the crystal structure of the vaterite polymorph of CaCO3: a calcium-43 solid-state NMR and computational assessment

Stacking Structure of Vaterite Revealed by Atomic Imaging and Diffraction Analysis

Anhydrous calcium carbonate crystals exist as three polymorphs: calcite, aragonite, and vaterite. Although vaterite is a metastable phase rarely found in the geological environment, it is intriguing that various biominerals are composed of vaterite. The processes of stable vaterite formation in biological systems cannot be understood without elucidating the nature of vaterite. The crystal structure of vaterite has been discussed for nearly a century but is still an open question. Here we propose the actual structure of vaterite by combining atomic imaging and diffraction analysis with simulations of disordered stacking sequences. Vaterite basically appears as layers of hexagonal calcium planes and carbonate (CO3 2-)-containing sheets stacked with +60°, -60°, or 180° rotations from the underlying layer. However, equivalent carbonate positions in alternating layers are forbidden, and four-layer stacking in which the fourth layer rotates 180° relative to the first layer are predominant, forming an orthogonal reciprocal lattice in diffraction patterns. These stacking characteristics replicate the intensity distribution in the electron and X-ray diffraction patterns. This study has almost completely elucidated the crystal structure and stacking sequence of vaterite. Our findings provide insights into the thermodynamic stability of vaterite, which facilitates comprehension of the biomineralization processes and growth dynamics of calcium carbonate.

Double-rotation (DOR) NMR spectroscopy: Progress and perspectives

Double-rotation (DOR) solid-state NMR spectroscopy is a high-resolution technique developed in the late 1980s. Although multiple-quantum magic-angle spinning (MQMAS) became the most widely used high-resolution method for half-integer spin quadrupoles after 1995, development and application of DOR NMR to a variety of chemical and materials science problems has endured. This Trend article recapitulates the development of DOR NMR, discusses various applications, and describes possible future directions. The main technical limitations specific to DOR NMR are simply related to the size of the double rotor system. The relatively large outer rotor (and thus coil) used for most applications over the past 35 years translates into relatively low rotor spinning frequencies, a low filling factor, and weak radiofrequency powers available for excitation and for proton decoupling. Ongoing developments in NMR instrumentation, including ever-shrinking MAS rotors and spherical NMR rotors, could solve many of these problems and may augur a renaissance for DOR NMR.

A facile approach for estimating radio-frequency field strength of low-receptivity nuclei

Radio-frequency (RF) field calibration is essential in NMR spectroscopy. A common practice is to collect a nutation curve by varying the pulse length in a direct single-pulse excitation experiment or in a cross-polarization magic-angle spinning with a flip-back pulse experiment. From the null points on this curve, one can calculate the RF field strength. Nevertheless, the practical implementation is not always straightforward or can even be unrealizable, especially for low-receptivity nuclei owing to their associated low sensitivity. Several researchers used an approach that involves utilizing other nuclei with more sensitivity but nearly identical Larmor frequencies to that of the nucleus of interest. However, such an approach has not been a common practice so far. In this work, we have systematically revisited this approach using 3.2 mm rotors on different sets of nuclei covering a Larmor frequency range up to 80 MHz. The effect of solid- and solution-states on RF field strength measurements has been investigated. The detection of each set of nuclei is then carried out with a resonant circuit in the NMR probe consisting of identical coils and capacitors. Our methodology is illustrated by recording 135/137Ba NMR spectra of BaTiO3 without prior 135/137Ba RF field calibration.

A 43 Ca nuclear magnetic resonance perspective on octacalcium phosphate and its hybrid derivatives

⁴³ Ca nuclear magnetic resonance (NMR) spectroscopy has been extensively applied to the detailed study of octacalcium phosphate (OCP), Ca₈ (HPO₄ )₂ (PO₄ )₄ .5H₂ O, and hybrid derivatives involving intercalated metabolic acids (viz., citrate, succinate, formate, and adipate). Such phases are of importance in the development of a better understanding of bone structure. High-resolution ⁴³ Ca magic angle spinning (MAS) experiments, including double-rotation (DOR) ⁴³ Ca NMR, as well as ⁴³ Ca{¹ H} rotational echo DOR (REDOR) and ³¹ P{⁴³ Ca} REAPDOR NMR spectra, were recorded on a ⁴³ Ca-labeled OCP phase at very high magnetic field (20 T), and complemented by ab initio calculations of NMR parameters using the Gauge-Including Projector Augmented Wave-density functional theory (GIPAW-DFT) method. This enabled a partial assignment of the eight inequivalent Ca²⁺ sites of OCP. Natural-abundance ⁴³ Ca MAS NMR spectra were then recorded for the hybrid organic-inorganic derivatives, revealing changes in the ⁴³ Ca lineshape. In the case of the citrate derivative, these could be interpreted on the basis of computational models of the structure. Overall, this study highlights the advantages of combining high-resolution ⁴³ Ca NMR experiments and computational modeling for studying complex hybrid biomaterials.

Recent progress in solid-state nuclear magnetic resonance of half-integer spin low-γ quadrupolar nuclei applied to inorganic materials

An overview is presented of recent progress in the solid-state nuclear magnetic resonance (NMR) observation of low-γ nuclei, with a focus on applications to inorganic materials. The technological and methodological advances in the last 20 years, which have underpinned the increased accessibility of low-γ nuclei for study by solid-state NMR techniques, are summarised, including improvements in hardware, pulse sequences and associated computational methods (e.g., first principles calculations and spectral simulation). Some of the key initial observations from inorganic materials of these nuclei are highlighted along with some recent (most within the last 10 years) illustrations of their application to such materials. A summary of other recent reviews of the study of low-γ nuclei by solid-state NMR is provided so that a comprehensive understanding of what has been achieved to date is available.

Predicting hydration layers on surfaces using deep learning

Characterisation of the nanoscale interface formed between minerals and water is essential to the understanding of natural processes, such as biomineralization, and to develop new technologies where function is dominated by the mineral-water interface. Atomic force microscopy offers the potential to characterize solid-liquid interfaces in high-resolution, with several experimental and theoretical studies offering molecular scale resolution by linking measurements directly to water density on the surface. However, the theoretical techniques used to interpret such results are computationally intensive and development of the approach has been limited by interpretation challenges. In this work, we develop a deep learning architecture to learn the solid-liquid interface of polymorphs of calcium carbonate, allowing for the rapid predictions of density profiles with reasonable accuracy.

The atomic structure and dynamics at the CaCO3 vaterite-water interface: A classical molecular dynamics study

Classical molecular and lattice dynamics were applied to explore the structure and dynamics of water on different surfaces of vaterite, the least abundant calcium carbonate polymorph. Surfaces were generated starting from the three possible structural models for vaterite (monoclinic, hexagonal/trigonal, and triclinic) and pre-screened using their surface energies in an implicit solvent. Surfaces with energies lower than 0.55 J/m² were then run in explicit water. The majority of these surfaces dissolve in less than 100 ns, highlighting the low stability of this phase in abiotic environments. Three stable surfaces were identified; they exhibited only minor structural changes when in contact with explicit water and did not show any tendency to dissolve during 1 µs molecular dynamics simulations. The computed water density profiles show that all these surfaces have two distinct hydration layers. The water residence time at the various calcium sites was computed to be within 0.7 and 20.5 ns, which suggests that specific Ca ions will be more readily available to bind with organic molecules present in solution. This analysis is a step forward in understanding the structure of this complex mineral and its role in biomineralization, as it provides a solid theoretical background to explore its surface chemistry. In particular, this study provides realistic surface models and predicts the effect of water exchange at the surface active sites on the adsorption of other molecules.

Stacking sequence variations in vaterite resolved by precession electron diffraction tomography using a unified superspace model

As a metastable phase, vaterite is involved in the first step of crystallization of several carbonate-forming systems including the two stable polymorphs calcite and aragonite. Its complete structural determination would consequently shed important light to understand scaling formation and biomineralization processes. While vaterite’s hexagonal substructure (a₀ ~ 4.1 Å and c₀ ~ 8.5 Å) and the organization of the carbonate groups within a single layer is known, conflicting interpretations regarding the stacking sequence remain and preclude the complete understanding of the structure. To resolve the ambiguities, we performed precession electron diffraction tomography (PEDT) to collect single crystal data from 100 K to the ambient temperature. The structure was solved ab initio and described over all the temperature range using a unified modulated structure model in the superspace group C12/c1(α0γ)00 with a = a₀ = 4.086(3) Å, b = \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{{\bf{3}}}$$\end{document}3a₀ = 7.089(9) Å, c = c₀ = 8.439(9) Å, α = β = γ = 90° and q = \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tfrac{{\bf{2}}}{{\bf{3}}}$$\end{document}23a* + γc*. At 100 K the model presents a pure 4-layer stacking sequence with γ = \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tfrac{{\bf{1}}}{{\bf{2}}}$$\end{document}12 whereas at the ambient temperature, ordered stacking faults are introduced leading to γ < \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tfrac{{\bf{1}}}{{\bf{2}}}$$\end{document}12. The model was refined against PEDT data using the dynamical refinement procedure including modulation and twinning as well as against x-ray powder data by the Rietveld refinement.

Synthesis Methods and Favorable Conditions for Spherical Vaterite Precipitation: A Review

Vaterite is the least thermodynamically stable anhydrous calcium carbonate polymorph. Its existence is very rare in nature, e.g., in some rock formations or as a component of biominerals produced by some fishes, crustaceans, or birds. Synthetic vaterite particles are proposed as carriers of active substances in medicines, additives in cosmetic preparations as well as adsorbents. Also, their utilization as a pump for microfluidic flow is also tested. In particular, vaterite particles produced as polycrystalline spheres have large potential for application. Various methods are proposed to precipitate vaterite particles, including the conventional solution-solution synthesis, gas-liquid method as well as special routes. Precipitation conditions should be carefully selected to obtain a high concentration of vaterite in all these methods. In this review, classical and new methods used for vaterite precipitation are presented. Furthermore, the key parameters affecting the formation of spherical vaterite are discussed.

A new NMR crystallographic approach to reveal the calcium local structure of atorvastatin calcium

We combine experimental and computational determination of ⁴³Ca solid-state NMR parameters (chemical shift tensors, quadrupolar coupling tensors, and Euler angles) to constrain the structure of the local calcium–ligand coordination environment.

Recent advances in solid-state nuclear magnetic resonance spectroscopy of exotic nuclei

We present a review of recent advances in solid-state nuclear magnetic resonance (SSNMR) studies of exotic nuclei. Exotic nuclei may be spin-1/2 or quadrupolar, and typically have low gyromagnetic ratios, low natural abundances, large quadrupole moments (when I > 1/2), or some combination of these properties, generally resulting in low receptivities and/or prohibitively broad line widths. Some nuclides are little studied for other reasons, also rendering them somewhat exotic. We first discuss some of the recent progress in pulse sequences and hardware development which continues to enable researchers to study new kinds of materials as well as previously unfeasible nuclei. This is followed by a survey of applications to a wide range of exotic nuclei (including e.g., ⁹Be, ²⁵Mg, ³³S, ³⁹K, ⁴³Ca, ^47/49Ti, ⁵³Cr, ⁵⁹Co, ⁶¹Ni, ⁶⁷Zn, ⁷³Ge, ⁷⁵As, ⁸⁷Sr, ¹¹⁵In, ¹¹⁹Sn, ^121/123Sb, ^135/137Ba, ^185/187Re, ²⁰⁹Bi), most of them quadrupolar. The scope of the review is the past ten years, i.e., 2007-2017.

Empirically testing vaterite structural models using neutron diffraction and thermal analysis

Otoliths, calcium carbonate (CaCO₃) ear bones, are among the most commonly used age and growth structures of fishes. Most fish otoliths are comprised of the most dense CaCO₃ polymorph, aragonite. Sturgeon otoliths, in contrast, have been characterized as the rare and structurally enigmatic polymorph, vaterite-a metastable polymorph of CaCO₃. Vaterite is an important material ranging from biomedical to personal care applications although its crystal structure is highly debated. We characterized the structure of Lake Sturgeon otoliths using thermal analysis and neutron powder diffraction, which is used non-destructively. We confirmed that while Lake Sturgeon otoliths are primarily composed of vaterite, they also contain the denser CaCO₃ polymorph, calcite. For the vaterite fraction, neutron diffraction data provide enhanced discrimination of the carbonate group compared to x-ray diffraction data, owing to the different relative neutron scattering lengths, and thus offer the opportunity to uniquely test the more than one dozen crystal structural models that have been proposed for vaterite. Of those, space group P6₅22 model, a = 7.1443(4)Å, c = 25.350(4)Å, V = 1121.5(2)ų provides the best fit to the neutron powder diffraction data, and allows for a structure refinement using rigid carbonate groups.

DFT Calculations with van der Waals Interactions of Hydrated Calcium Carbonate Crystals CaCO3·(H2O, 6H2O): Structural, Electronic, Optical, and Vibrational Properties

The role of hydration on the structural, electronic, optical, and vibrational properties of monohydrated (CaCO3·H2O, hexagonal, P31, Z = 9) and hexahydrated (CaCO3·6H2O, monoclinic, C2/c, Z = 4) calcite crystals is assessed with the help of published experimental and theoretical data applying density functional theory within the generalized gradient approximation and a dispersion correction scheme. We show that the presence of water increases the main band gap of monohydrocalcite by 0.4 eV relative to the anhydrous structure, although practically not changing the hexahydrocalcite band gap. The gap type, however, is modified from indirect to direct as one switches from the monohydrated to the hexahydrated crystal. A good agreement was obtained between the simulated vibrational infrared and Raman spectra and the experimental data, with an infrared signature of hexahydrocalcite relative to monohydrocalcite being observed at 837 cm(-1). Other important vibrational signatures of the lattice, water molecules, and CO3(2-) were identified as well. Analysis of the phonon dispersion curves shows that, as the hydration level of calcite increases, the longitudinal optical-transverse optical phonon splitting becomes smaller. The thermodynamics properties of hexahydrocalcite as a function of temperature resemble closely those of calcite, while monohydrocalcite exhibits a very distinct behavior.

Topological features in crystal structures: a quotient graph assisted analysis of underlying nets and their embeddings

Topological properties of crystal structures may be analysed at different levels, depending on the representation and the topology that has been assigned to the crystal. Considered here is thecombinatorialorbond topologyof the structure, which is independent of its realization in space. Periodic nets representing one-dimensional complexes, or the associated graphs, characterize the skeleton of chemical bonds within the crystal. Since periodic nets can be faithfully represented by their labelled quotient graphs, it may be inferred that their topological features can be recovered by a direct analysis of the labelled quotient graph. Evidence is given for ring analysis and structure decomposition into building units and building networks. An algebraic treatment is developed for ring analysis and thoroughly applied to a description of coesite. Building units can be finite or infinite, corresponding to 1-, 2- or even 3-periodic subnets. The list of infinite units includes linear chains or sheets of corner- or edge-sharing polyhedra. Decomposing periodic nets into their building units relies on graph-theoretical methods classified assurgery techniques. The most relevant operations are edge subdivision, vertex identification, edge contraction and decoration. Instead, these operations can be performed on labelled quotient graphs, evidencing in almost a mechanical way the nature and connection mode of building units in the derived net. Various examples are discussed, ranging from finite building blocks to 3-periodic subnets. Among others, the structures of strontium oxychloride, spinel, lithiophilite and garnet are addressed.

The role of solid-state nuclear magnetic resonance in crystal engineering

This Highlight article discusses the role of solid-state NMR spectroscopy in crystal engineering with the aid of several examples from the literature.

Solid-State NMR at the University of Ottawa

This article describes some highlights of the research which has been carried out in my laboratory at the University of Ottawa over the period covering 2005 to 2014. My research is in the general areas of solid-state NMR, applications of quantum chemistry, and biomolecular NMR. The format will follow that of my 2014 Canadian Society for Chemistry Keith Laidler Award presentation given in Vancouver in June 2014 at the 97th Canadian Chemistry Conference and Exhibition. Following a brief introduction, I will present some of our most interesting and exciting recent advances according to the following six themes: 1. Fundamental solid-state NMR. 2. Materials characterization and NMR crystallography. 3. Pharmaceuticals and polymorphism. 4. Non-covalent interactions: Halogen bonds. 5. Biomolecular NMR. 6. Software development.