Molecular Orbital Study of Apatite (Ca5(PO4)3OH) Nucleation at Silica Bioceramic Surfaces

Silica bioceramics are used as prosthetic bone and dental implants because they promote apatite formation at their surfaces when immersed in simulated body fluids of composition similar to human blood plasma. Apatite formation occurs in stages, but the reaction pathway remains unresolved. We have used molecular orbital calculations to model the interactions of Ca 2+, H 2 PO 4 -, HPO 4 2-, and H 2 O with bioceramic surface sites represented by Si 7 O 12 H 10, Si 4 O 8 H 8, and Si 3 O 6 H 6 to determine the reaction sequence for apatite nucleation. Predicted reaction energies, vibrational frequencies, and 29 Si, 31 P, and 43 Ca NMR shifts were used as probes of the stable intermediates during reaction progress. Energies and vibrational frequencies were calculated using effective core potentials and valence double-ζ basis sets. Most NMR shifts were calculated using the 3-21G* basis, and using the 6-31G* basis for selected clusters. With increasing pH above the point of zero charge of silica, (pzc ≈ 3) the stability of the calcium surface complexes on the Si 3-ring ([Si 3 O 6 H 6 ]) is predicted to be outer-sphere < inner-sphere < surface-hydroxide. Comparison of calculated 29 Si and 31 P NMR shielding trends of Si 3- and 4-rings (Si 3 O 6 H 6 and Si 4 O 8 H 8, respectively) with experimental trends for bioceramic reacting with SBF suggest that (i) the partially deprotonated 3-ring is the active surface site, because it promotes calcium ion dehydration as calcium adsorbs, (ii) formation of [SiO−Ca−OPO 3 H] bonds is energetically preferred over direct Si−O−P bonds, and (iii) an acidic precursor with bidentate Ca>OPO 3 H bonds nucleates rapidly (minutes to 1 h) at the bioceramic surface. The precursor oligomer is modeled by the complex [Si 3 O 6 H 5 CaHPO 4 (H 2 O) 3 ] - . The actual precursor is proposed to contain two or three calcium hydrogen phosphate units with a stoichiometry such as [Si 3 O 6 H 5 CaHPO 4 (H 2 O) n ] -, where n = 2 or 3. The complex [Si 3 O 6 H 5 CaHPO 4 (H 2 O) 3 ] -, young bone, and bioceramic reacting with simulated body fluid share unique infrared/Raman bands at 631 and 1125−1145 cm -1, which are distinct from the bands observed in crystalline apatite and mature bone. Predicted 43 Ca NMR shifts provide avenues for future experiments.