Electronic structure, bonding, charge distribution, and x-ray absorption spectra of the (001) surfaces of fluorapatite and hydroxyapatite from first principles

Fluorapatite (FAP) and hydroxyapatite (HAP) are two very important bioceramic crystals. The (001) surfaces of FAP and HAP crystals are studied by ab initio density functional calculations using a supercell slab geometry. It is shown that in both crystals, the O-terminated (001) surface is more stable with calculated surface energies of 0.865 and 0.871 J/m(2) for FAP and HAP, respectively. In FAP, the two surfaces are symmetric. In HAP, the orientation of the OH group along the c axis reduces the symmetry such that the top and bottom surfaces are no longer symmetric. It is revealed that the atoms near the surface and subsurface are significantly relaxed especially in the case of HAP. The largest relaxations occurred via the lateral movements of the O ions at the subsurface level. The electronic structures of the surface models in the form of layer-by-layer resolved partial density of states for all the atoms show systematic variation from the surface region toward the bulk region. The calculated Mulliken effective charge on each type of atom and the bond order values between cations (Ca, P) and anions (O, F) show different charge transfers and bond strength variations from the bulk crystal values. Electron charge density calculations show that the surfaces of both FAP and HAP crystals are mostly positively charged due to the presence of Ca ions at the surface. The positively charged surfaces have implications for the absorption on apatite surfaces of water and other organic molecules in an aqueous environment which are an important part of its bioactivity. The x-ray absorption near-edge structure (XANES) spectra (Ca-K, O-K, F-K, P-K, and P-L-3 edges) of both the surface models and the bulk crystals are calculated and compared. The calculations use a supercell approach which takes into account the electron-core-hole interaction. It is shown that the site-specific XANES spectra show significant differences between atoms near the surface and in the bulk and are very sensitive to the local atomic environment of each atom. This information will be very valuable for characterizing the apatite materials and in the interpretation of experimental data. Comparisons of several sets of experimental data with the weighted sums of the calculated spectra at different sites for the same element show very good agreement.