Calcium carbonate is an abundant substance that can be created in several mineral forms by the reaction of dissolved carbon dioxide in water with calcium ions. Through biomineralization, organisms can harness and control this process to form various functional materials that can act as anything from shells through to lenses. The early stages of calcium carbonate formation have recently attracted attention as stable prenucleation clusters have been observed, contrary to classical models. Here we show, using computer simulations combined with the analysis of experimental data, that these mineral clusters are made of an ionic polymer, composed of alternating calcium and carbonate ions, with a dynamic topology consisting of chains, branches and rings. The existence of a disordered, flexible and strongly hydrated precursor provides a basis for explaining the formation of other liquid-like amorphous states of calcium carbonate, in addition to the non-classical behaviour during growth of amorphous calcium carbonate.

CRYSTAL is a periodic ab initio code that uses a Gaussian-type basis set to express crystalline orbitals (i.e., Bloch functions). The use of atom-centered basis functions allows treating 3D (crystals), 2D (slabs), 1D (polymers), and 0D (molecules) systems on the same grounds. In turn, all-electron calculations are inherently permitted along with pseudopotential strategies. A variety of density functionals are implemented, including global and range-separated hybrids of various natures and, as an extreme case, Hartree–Fock (HF). The cost for HF or hybrids is only about 3–5 times higher than when using the local density approximation or the generalized gradient approximation. Symmetry is fully exploited at all steps of the calculation. Many tools are available to modify the structure as given in input and simplify the construction of complicated objects, such as slabs, nanotubes, molecules, and clusters. Many tensorial properties can be evaluated by using a single input keyword: elastic, piezoelectric, photoelastic, dielectric, first and second hyperpolarizabilities, etc. The calculation of infrared and Raman spectra is available, and the intensities are computed analytically. Automated tools are available for the generation of the relevant configurations of solid solutions and/or disordered systems. Three versions of the code exist: serial, parallel, and massive-parallel. In the second one, the most relevant matrices are duplicated on each core, whereas in the third one, the Fock matrix is distributed for diagonalization. All the relevant vectors are dynamically allocated and deallocated after use, making the code very agile. CRYSTAL can be used efficiently on high performance computing machines up to thousands of cores.

Both of the previously proposed Pbnm and P6₅22 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 P3₂21 symmetry. Since interconversion between these structures requires only thermal energy, this provides an additional source of disorder within the vaterite structure.