  Mineralization of bone and teeth involves interactions between biomolecules
and hydroxyapatite. Associated complex interfaces and processes remain
difficult to analyze at the 1 to 100 nm scale using current laboratory
techniques, and prior models for atomistic simulations are limited in the
representation of chemical bonding, surface chemistry, and interfacial
interactions. This work introduces an accurate force field along with
pH-resolved surface models for hydroxyapatite to represent chemical bonding,
structural, surface, interfacial, and mechanical properties in quantitative
agreement with experiment. The accuracy is orders of magnitude higher in
comparison to earlier models to facilitate quantitative monitoring of
inorganic-biological assembly. The force field is integrated into the CHARMM,
AMBER, OPLS-AA, PCFF, and INTERFACE force fields to enable realistic
simulations of apatite-biological systems of any composition and ionic
strength. Specifically, the parameters reproduce lattice constants (<0.5%
deviation), IR spectrum, cleavage energies, immersion energies in water (<5%
deviation), and elastic constants (<10% deviation) of hydroxyapatite in
comparison to experiment. Interactions between mineral, water, and organic
compounds are represented by standard combination rules in the force field
without additional adjustable parameters and shown to achieve quantitative
accuracy. Surface models for common (001), (010), (020), (101) facets and
nanocrystals are introduced as a function of pH on the basis of extensive
experimental data. New insight into surface and immersion energies, the
structure of aqueous interfaces, density profiles, and superficial dissolution
is described. Mechanisms of specific binding of peptides, drugs, and
mineralization can be analyzed and the force field is extensible to substituted
and defective apatites as well as to other calcium phosphate phases.
