Development of Accurate Force Fields for the Simulation of Biomineralization. Research Methods in Biomineralization Science. Elsevier; 2013. pp. 3-23. [DOI: 10.1016/b978-0-12-416617-2.00001-1]

Simulation of Calcium Phosphate Prenucleation Clusters in Aqueous Solution: Association beyond Ion Pairing

Classical molecular dynamics simulations and free energy methods have been used to obtain a better understanding of the molecular processes occurring prior to the first nucleation event for calcium phosphate biominerals. The association constants for the formation of negatively charged complexes containing calcium and phosphate ions in aqueous solution have been computed, and these results suggest that the previously proposed calcium phosphate building unit, [Ca(HPO4)3]4-, should only be present in small amounts under normal experimental conditions. However, the presence of an activation barrier for the removal of an HPO4 2- ion from this complex indicates that this species could be kinetically trapped. Aggregation pathways involving CaHPO4, [Ca(HPO4)2]2-, and [Ca(HPO4)3]4- complexes have been explored with the finding that dimerization is favorable up to a Ca/HPO4 ratio of 1:2.

Structures and infrared spectra of calcium phosphate clusters by ab initio methods with implicit solvation models

Since the first detection of pre-nucleation clusters during the formation of calcium phosphate minerals, determining such clusters' compositions and structures has become crucial for understanding the early-stage nucleation of these minerals in solutions. In previous experimental studies, the composition and sizes of pre-nucleation clusters have been calculated, but their structural information has been difficult to determine because they are very small (<1 nm). In this study, we examined the structures and infrared spectra of small- and medium-sized calcium phosphate clusters using ab initio calculations combined with implicit solvation models. Adding solvent effects increased the possibility of the existence of alternative configurations of calcium phosphate clusters other than their compact configurations. The calcium atoms had a tendency to be located outside of the clusters to coordinate with water molecules in the aqueous environment. The computed infrared spectra of extended small calcium phosphate clusters captured some of the features measured in the in situ infrared spectra, which supports the network structures proposed by large-scale molecular dynamics studies and X-ray adsorption near-edge spectra. The relative stabilities of medium-sized Ca₉(PO₄)₆ clusters with respect to the stability of Posner's cluster in water were also reviewed. We found that in water, alternative structures with low symmetry or large dipole moments had lower energies than Posner's cluster.