Research Article  Open Access
Ricardo D. S. Santos, Marcos V. dos S. Rezende, "Atomistic Simulation of Intrinsic Defects and Trivalent and Tetravalent Ion Doping in Hydroxyapatite", Advances in Condensed Matter Physics, vol. 2014, Article ID 609024, 8 pages, 2014. https://doi.org/10.1155/2014/609024
Atomistic Simulation of Intrinsic Defects and Trivalent and Tetravalent Ion Doping in Hydroxyapatite
Abstract
Atomistic simulation techniques have been employed in order to investigate key issues related to intrinsic defects and a variety of dopants from trivalent and tetravalent ions. The most favorable intrinsic defect is determined to be a scheme involving calcium and hydroxyl vacancies. It is found that trivalent ions have an energetic preference for the Ca site, while tetravalent ions can enter P sites. Charge compensation is predicted to occur basically via three schemes. In general, the charge compensation via the formation of calcium vacancies is more favorable. Trivalent dopant ions are more stable than tetravalent dopants.
1. Introduction
Hydroxyapatite () has been intensively investigated due to the possibility of applying it in a range of biomedical applications. This is attributed to characteristics such as its mechanical properties, biocompatibility, osteoconductivity, nontoxic properties, noninflammatory properties, and the similarity of the mineral constituents to human bone and teeth [1–6]. Due to these properties, it can be applied in the manufacture of artificial bone material, as a coating on surgical implants, and in tissue engineering, drug and gene delivery, and other biomedical areas [7–11]. Hydroxyapatite has a hexagonal crystal structure, with space group P63/m, where the experimental dimensions of the unit cell are , , , and [12]. There are two crystallographically different sites for , which have different coordination numbers. Ca(I) is coordinated with nine oxygen ions and Ca(II) is coordinated with six oxygen ions and one from the hydroxyl group. The average Ca–O distances are lesser than 3 Å.
Many works have investigated the substitution of extrinsic dopant in hydroxyapatite (HAP) in order to improve the material properties [13–17]. The incorporation of dopant in the HAP structure contributes to increasing a range of potential applications in fields such as water purification, bone pathologies, bioceramics, catalysis, and luminescence [18]. For example, it is reported that doped HAP has useful applications in fluorescent lamps [19]. was incorporated as a dopant in the material with the aim of developing a biosensor [20]. HAP can also be doped by trivalent lanthanide ions such as [21, 22], [23, 24], [25], , and [26, 27]. The incorporation of Eu and Gd in a HAP structure can affect the bone remodeling cycle, and they have potential for the treatment of bone density disorders such as osteoporosis [26, 27]. Eu ions exhibit favorable spectroscopic properties and can also be used for biological imaging [28, 29] and in lasers. On the other hand, gadolinium can be used as a contrast agent to provide a brighter magnetic resonance (MR) imaging signal [30].
A systematic investigation of cation substitution in HAP may provide a better understanding of their properties, which is still lacking. Some theoretical studies based on first principle calculations for HAP materials have been realized in order to elucidate some basic physics [31–33]. However, many questions regarding the study of defects remain open. Dopant ions with 3+ and 4+ valence states such as rare earth ions and some transition metals will be incorporated into the Ca or P sites. The compensation mechanism for this doping remains inconclusive. For this reason, this paper makes use of computer modeling methods in order to predict the preferred dopant sites and charge compensation mechanisms in the material.
For Eudoped HAP, for example, an investigation of the structural response is essential due to the intimate relationship between its properties and lattice distortions caused by the Eu substitution and the charge compensation involved in the Eu emission, which are associated with the crystal field [34]. It is well known that the emission spectrum of strongly depends on the symmetry of the site which occupies. If it occupies a site with inversion symmetry, only the magneticdipole transition can be observed; but if there is no inversion symmetry at the site of the ion, the electricdipole transitions can be observed. Also, the transition is sensitive to the coordination environment [35].
In the present work, atomistic simulation based on the lattice energy minimization is used to provide useful information on the energetically favoured intrinsic defects and a variety of dopants from trivalent and tetravalent ions. Specifically, trivalent rare earth ions and some trivalent and tetravalent transition metals were considered due to their influence on various properties [31–33]. A detailed investigation on the defects behavior of HAP is very important in order to understand their importance on the improvement of the material properties. In these methods, detailed estimation of lattice relaxation and the Coulomb energies of a large number of lattice ions around defect species are provided without high computational demanding, on the contrary of the DFT based methods.
2. Methodology
The simulation techniques used in this work are based on energy minimization, with interactions represented by interatomic potentials. Interactions between ions in the solid are represented in terms of a longrange Coulomb term plus a shortrange term, as described by Buckingham potential, that accounts for electron cloud overlap Pauli repulsion and dispersion (Van der Waals) interactions. Polarizability of the oxygen ions is incorporated by means of the DickOverhauser shell model [36]. Most of the interatomic potential models that were obtained empirically for the simulation of hydroxyapatite materials [37–41] showed excellent agreement with the experimental structures and properties. However, in order to improve the fit with respect to structural and mechanical properties, the interatomic potential initially derived by Mostafa and Brown [42] was refitted empirically using the GULP code [43]. In our refit, the and parameters for P–O interaction from Mostafa and Brown [42] were modified. In our work, the tree body potential for P–O interaction used by Mostafa and Brown [42] was not considered. Defects are modeled using the MottLittleton approximation [44], in which a spherical region of lattice surrounding the defect is treated explicitly, with all interactions being considered, and more distant parts of the lattice are treated using a continuum approach.
3. Results
Suitable potential parameters are highly necessary in order to give a good description of the material by atomistic simulation. The refitted potential parameters used in this work are listed in Table 1. The obtained structural parameters are compared with the experimental values and other theoretical calculations in Table 2. It can be seen that the lattice parameters obtained in this work are more accurate than those obtained previously using atomistic simulation with the refit procedure [42] or using other force fields with the partial charge model [45–47] or those obtained using density function (DFT) calculation [48]. Our set of potentials reproduces lattice parameter to within 0.1%, that is, lesser than the other theoretical works. Available elastic properties of HAP materials are used to validate and refine the potential parameters. Table 3 shows the agreement between the calculated and experimental values of the elastic constants for the HAP structure. It can be seen that our obtained values are closer to those of experimental works and other previous theoretical works. The variability in the experimental elastic constants reported in the literature is shown in Table 3.



To determine Frenkel and Schottky type defect formation energies, isolated point defect (vacancy and interstitial) energies and relevant lattice energies were first calculated. Several possible positions were tested to confirm the optimal position of the interstitial site for defect occupancy, and the positions which have the lowest energy were taken for interstitial cations and oxygen, respectively. We also calculated antisite pair defects, which involve the exchange of a cation with another cation of different species. In addition, two other intrinsic defect schemes are considered. Scheme (i) involves one calcium vacancy and two hydroxyl vacancies, and scheme (ii) involves five calcium vacancies and two phosphorus interstitials. These defects are expressed by KrögerVink notation and are shown in Table 4. The calculated defect formation energies, also listed in Table 4, were obtained by combining the energies of these point defects. The results indicate that the formation of defects involving phosphorus is generally unfavorable due to the high energy required to create this type of defect. The defects involving P ions are highly unfavorable due to strong P–O Coulombic interaction; that is, they would hardly occur compared to the other types of defects. The result also reveals a low energy (−0.71 eV) for scheme (i). This is consistent with the theoretical observation reported by de Leeuw et al. [31] using DFT calculation. The work of de Leeuw et al. [31] reveals that calcium vacancies compensated by hydroxyl are preferred to those compensated by phosphate vacancies. They also show that the calcium vacancies compensated by substitutional defects, such as carbonate groups, are more favorable for compensation by hydroxyl and that if the carbonate defect is accompanied by the substitution of a monovalent or ion for the calcium ion, the defect formation energies are lower. In the more favorable schemes, calcium vacancies are involved, that is, one defect related to the calcium deficiency in the hydroxyapatite materials. Experimental works reveal that apatites are often found to be deficient in calcium compared to the stoichiometric material [51–53]. Other mechanisms of defects, not considered by de Leeuw et al. [31], could be involved in the calcium deficiency in the hydroxyapatite. From our results, it can be seen that calcium pseudoSchottky and calcium Frenkel defects are the next most favorable defects. In both schemes, calcium vacancies are involved. Thus, these defects also contribute to the higher concentration of calcium vacancy in HAP. This suggests that a higher concentration of calcium and hydroxyl vacancies could be present, which may contribute to calcium deficiency in the hydroxyapatite materials.

In the doping process, a description of the favorable substitution site and the charge compensation mechanism is very important information. From atomistic simulation it is possible to obtain quantitative estimates of the relative energies of different modes of dopant substitution. For this, trivalent () and tetravalent () dopants ions substituted at both Ca and P sites were investigated in order to help explain experimental results. The charge compensating mechanisms involved in the incorporation of and ions have not been clearly established, so all possible schemes are considered in this work.
and dopant ions can, in principle, be incorporated into the lattice at either or sites. In both cases, there is more than one possible mode of charge compensation. All possible charge compensation mechanisms were considered and are shown as follows.
Reaction Schemes for Solution of Divalent, Trivalent, Tetravalent, and Pentavalent Dopants. Consider the following mechanism.
Trivalent dopants
Tetravalent dopants In the crystal structure of HAP [12], there are two types of host site. Therefore, the incorporation in each site is considered. Seven reaction schemes have been considered for the incorporation of the trivalent ions and seven for the tetravalent ions.
In (1) to (4) the trivalent cation metals are considered at the site, and in (5) to (7) the trivalent cation metals are considered at the site. On the other hand, in (8) to (11) the tetravalent cation metals are considered at the site, and in (12) to (14) the tetravalent cation metals are considered at the site. In the all case, charge compensation is needed.
The solution energies for the trivalent defect were calculated by combining the appropriate defect and lattice energy terms and are listed in Tables 5–7. It is seen that the solution energy is different for the incorporation of the earth rare/metal transition with different ionic radius. The latter observation is explained by the lattice expansion/contraction, resulting in more/less space to accommodate earth rare/metal transition ions together with any charge compensating defects. From the tables it can be seen that substitution at the Ca site is more favorable for all trivalent dopants. The Ca site preference has been confirmed experimentally by Martin et al. [54] and Wakamura et al. [55]. This behavior can be explained in terms of the difference in ionic radii between the trivalent dopant ions and the host site ions. The ionic radii of transition metal and rare earth ions vary around 0.55–0.75 and 0.86–1.03 Å, respectively, and the ionic radii of and are 1.00 and 0.38 Å, respectively [56]. The small difference between the trivalent dopant ions and ions contributes to a small deformation in the lattice and consequently a small solution energy. Thus, the results show a degree of correlation between the dopant size and the solution energy. Dopants with a large difference in ion size from the host tend to be more energetically unfavorable. Between the two different Ca sites, it seems that, except for Lu, Er, Gd, Eu, Sm, and Nd show a preference for the Ca1 site substitution.



In the Eudoped HAP, the transition of Eu emission is particularly informative. This is because the state is unsplit, and the ground state is also unsplit, so that transitions to it give straightforward information about the excited state. If more than one component is seen for this transition, it shows that there is more than one europium site [34]. The two peaks attributed to transition were observed in the Eu emission spectra by Jagannathan and Kottaisamy [57] and Martin et al. [54], indicating that the Eu ion can be substituted at two sites. From our results, it can be noted that the solution energy for more favorable schemes for Eu substitution at the Ca1 site (1.79 eV) is close to that for the substitution at the Ca2 site (1.50 eV). This small difference justifies the possibility of Eu substitution at both host sites. The results of Martin et al. [54] support our hypothesis that Eu dopant ions prefer to be substituted at the Ca2 site due to the lower solution energy. They also report that the site occupancy ratio between the Ca2 and Ca1 sites is about 80%, with the majority of the ions being in the Ca2 sites.
The charge compensation tends to occur basically via three schemes. In general, charge compensation via the formation of calcium vacancies is more favorable for charge compensation in the trivalent dopants, except for Er, Gd, and Eu, where the charge compensation is more favorable through interstitial hydroxyl than by calcium vacancies, and for Sm and Nd, where charge compensation by interstitial oxygen is more favorable. The charge compensation for Eu ions is in agreement with that proposed by Martin et al. [54] and Ternane et al. [14].
The production of stoichiometric HAP is obtained with a Ca/P ratio of 1.67. Experimental work realized by Mayer et al. [58] reports that this value is below 1.67 in most cases and shows a decrease with increases in the La content in the samples; that is, it indicates the formation of nonstoichiometric HAP. This reduction of the Ca/P ratio is attributed to the formation of the calcium vacancy generated as a charge compensation defect needed in the substitution of the La ion in the Ca site. The occurrence of this scheme (see (2)) is the most favorable, as can be seen in Table 7. Similar behavior can be expected for all the trivalent dopants that are compensated by calcium vacancies.
In Table 8, the solution energies for the tetravalent defect calculated by combining the appropriate defect and lattice energy terms are listed. It can be seen that substitution at the site by a compensated oxygen vacancy is more favorable for and dopants. This behavior also can be explained in terms of similar ionic radii between the tetravalent dopant ions and the host site ions. Mn and Cr dopants can exist as either or . It is necessary to know their precise oxidation state. According to our results, the favorable doping process for these polyvalent dopants should be an form via the formation of calcium vacancies.
4. Conclusion
The basic defect chemistry including intrinsic defects and trivalent and tetravalent extrinsic dopants has been investigated using atomistic simulation. The most favorable intrinsic defect schemes are formed by one calcium vacancy and two hydroxyl vacancies. Trivalent and tetravalent dopant ion substitutions are found to take place preferentially on the Ca and P sites, respectively. In both cases, calcium vacancy defect formation is the more likely charge compensation mechanism for most of the trivalent and tetravalent dopants. In some trivalent ions, charge compensation by interstitial hydroxyl and interstitial oxygen is more favorable. The favorable doping process for these polyvalent dopants should be an form via the formation of calcium vacancies.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors would like to acknowledge the financial support by FINEP, CAPES, and CNPq.
References
 A. Costescu, I. Pasuk, F. Ungureanu et al., “Physicochemical properties of nanosized hexagonal hydroxyapatite powder synthesized by solgel,” Digest Journal of Nanomaterials and Biostructures, vol. 5, no. 4, pp. 989–1000, 2010. View at: Google Scholar
 C. S. Ciobanu, E. Andronescu, A. Stoicu et al., “Influence of annealing treatment of nanohydroxyapatite bioceramics on the vibrational properties,” Digest Journal of Nanomaterials and Biostructures, vol. 6, no. 2, pp. 609–624, 2011. View at: Google Scholar
 C. S. Ciobanu, E. Andronescu, B. S. Vasile, C. M. Valsangiacom, R. V. Ghita, and D. Predoi, “Looking for new synthesis of hydroxyapatite doped with europium,” Optoelectronics and Advanced Materials, Rapid Communications, vol. 4, no. 10, pp. 1515–1519, 2010. View at: Google Scholar
 D. Predoi, S. Derible, and H. Duflo, “Synthesis and ultrasonic characterization of hydroxyapatite ceramic powders,” Journal of Optoelectronics and Advanced Materials, vol. 11, no. 6, pp. 852–856, 2009. View at: Google Scholar
 D. Predoi, R. A. VatasescuBalcan, I. Pasuk, R. Trusca, and M. Costache, “Calcium phosphate ceramics for biomedical applications,” Journal of Optoelectronics and Advanced Materials, vol. 10, no. 8, pp. 2151–2155, 2008. View at: Google Scholar
 T. S. B. Narasaraju and D. E. Phebe, “Some physicochemical aspects of hydroxylapatite,” Journal of Materials Science, vol. 31, no. 1, pp. 1–21, 1996. View at: Google Scholar
 S. I. Stupp and P. V. Braun, “Molecular manipulation of microstructures: biomaterials, ceramics, and semiconductors,” Science, vol. 277, no. 5330, pp. 1242–1248, 1997. View at: Publisher Site  Google Scholar
 W. Suchanek and M. Yoshimura, “Processing and properties of hydroxyapatitebased biomaterials for use as hard tissue replacement implants,” Journal of Materials Research, vol. 13, pp. 94–117, 1998. View at: Publisher Site  Google Scholar
 M. Y. Ma, Y. J. Zhu, L. Li, and S. W. Cao, “Nanostructured porous hollow ellipsoidal capsules of hydroxyapatite and calcium silicate: preparation and application in drug delivery,” Journal of Materials Chemistry, vol. 18, pp. 2722–2727, 2008. View at: Publisher Site  Google Scholar
 A. Maitra, “Calcium phosphate nanoparticles: secondgeneration nonviral vectors in gene therapy,” Expert Review of Molecular Diagnostics, vol. 5, no. 6, pp. 893–905, 2005. View at: Publisher Site  Google Scholar
 Y. R. Cai and R. K. Tang, “Calcium phosphate nanoparticles in biomineralization and biomaterials,” Journal of Materials Chemistry, vol. 18, no. 32, pp. 3775–3787, 2008. View at: Publisher Site  Google Scholar
 Y. Tanaka, T. Iwasaki, K. Katayama, J. Hojo, and K. Yamashita, “Effect of ionic polarization on crystal structure of hydroxyapatite ceramic with hydroxide nonstoichiometry,” Journal of the Japan Society of Powder and Powder Metallurgy, vol. 57, no. 7, pp. 520–528, 2010. View at: Publisher Site  Google Scholar
 D. Shi, J. Lian, W. Wang et al., “Luminescent carbon nanotubes by surface functionalization,” Advanced Materials, vol. 18, no. 2, pp. 189–193, 2006. View at: Publisher Site  Google Scholar
 R. Ternane, M. TrabelsiAyedi, N. KbirAriguib, and B. Piriou, “Luminescent properties of Eu^{3+} in calcium hydroxyapatite,” Journal of Luminescence, vol. 81, no. 3, pp. 165–170, 1999. View at: Publisher Site  Google Scholar
 M. Gaft, R. Reisfeld, G. Panczer, S. Shoval, B. Champagnon, and G. Boulon, “Eu3+ luminescence in highsymmetry sites of natural apatite,” Journal of Luminescence, vol. 72–74, pp. 572–574, 1997. View at: Publisher Site  Google Scholar
 D. Predoi, M. Barsan, E. Andronescu, R. A. VatasescuBalcan, and M. Costache, “Hydroxyapatiteiron oxide bioceramic prepared using nanosize powders,” Journal of Optoelectronics and Advanced Materials, vol. 9, no. 11, pp. 3609–3613, 2007. View at: Google Scholar
 D. Predoi, R. V. Ghita, F. Ungureanu, C. C. Negrila, R. A. VatasescuBalcan, and M. Costache, “Characteristics of hydroxyapatite thin films,” Journal of Optoelectronics and Advanced Materials, vol. 9, no. 12, pp. 3827–3831, 2007. View at: Google Scholar
 A. Zounani, D. Zambon, and J. C. Cousseins, “Optical properties of Eu^{3+} activated Sr_{10}F_{2} (PO_{4}) _{6} elaborated by coprecipitation,” Journal of Alloys and Compounds, vol. 188, pp. 82–86, 1992. View at: Publisher Site  Google Scholar
 R. G. Pappalardo, J. Walsh, and R. B. Hunt Jr., “Cerium‐activated halophosphate phosphors,” Journal of the Electrochemical Society, vol. 130, no. 10, pp. 2087–2096, 1983. View at: Publisher Site  Google Scholar
 T. S. de Araujo, Z. S. MacEdo, P. A. S. C. de Oliveira, and M. E. G. Valerio, “Production and characterization of pure and Cr^{3+}doped hydroxyapatite for biomedical applications as fluorescent probes,” Journal of Materials Science, vol. 42, no. 7, pp. 2236–2243, 2007. View at: Publisher Site  Google Scholar
 M. E. Fleet and Y. Pan, “Site preference of Nd in fluorapatite [Ca_{10}(PO_{4})_{6}F_{2}],” Journal of Solid State Chemistry, vol. 112, no. 1, pp. 78–81, 1994. View at: Publisher Site  Google Scholar
 F. M. Ryan, R. W. Warren, R. H. Hopkins, and J. Murphy, “Selective site laser excitation and ESR studies of Nd^{3+} Ions in Ca_{5} (PO_{4})_{3}F ,” Journal of the Electrochemical Society, vol. 125, no. 9, pp. 1493–1498, 1978. View at: Publisher Site  Google Scholar
 L. D. DeLoach, S. A. Payne, L. L. Chase, L. K. Smith, W. L. Kway, and W. F. Krupke, “Evaluation of absorption and emission properties of Yb^{3+} doped crystals for laser applications,” IEEE Journal of Quantum Electronics, vol. 29, no. 4, pp. 1179–1191, 1993. View at: Publisher Site  Google Scholar
 L. D. DeLoach, S. A. Payne, W. L. Kway, J. B. Tassano, S. N. Dixit, and W. F. Krupke, “Vibrational structure in the emission spectra of Yb3+doped apatite crystals,” Journal of Luminescence, vol. 62, no. 2, pp. 85–94, 1994. View at: Publisher Site  Google Scholar
 K. Spariosu, R. D. Stultz, M. Birnbaum, T. H. Allik, and J. A. Hutchinson, “Er:Ca_{5}(PO_{4})_{3}F saturableabsorber Q switch for the Er:glass laser at 1.53 μm,” Applied Physics Letters, vol. 62, no. 22, pp. 2763–2765, 1993. View at: Publisher Site  Google Scholar
 C. A. Barta, K. SachsBarrable, J. Jia, K. H. Thompson, K. M. Wasan, and C. Orvig, “Lanthanide containing compounds for therapeutic care in bone resorption disorders,” Dalton Transactions, no. 43, pp. 5019–5030, 2007. View at: Publisher Site  Google Scholar
 T. J. Webster, E. A. MassaSchlueter, J. L. Smith, and E. B. Slamovich, “Osteoblast response to hydroxyapatite doped with divalent and trivalent cations,” Biomaterials, vol. 25, no. 11, pp. 2111–2121, 2004. View at: Publisher Site  Google Scholar
 P. P. Yang, Z. W. Quan, C. X. Li, X. J. Kang, H. Z. Lian, and J. Lin, “Bioactive, luminescent and mesoporous europiumdoped hydroxyapatite as a drug carrier,” Biomaterials, vol. 29, no. 32, pp. 4341–4347, 2008. View at: Publisher Site  Google Scholar
 Y. C. Han, X. Y. Wang, and S. P. Li, “Biocompatible europium doped hydroxyapatite nanoparticles as a biological fluorescent probe,” Current Nanoscience, vol. 6, no. 2, pp. 178–183, 2010. View at: Google Scholar
 M. Bottrill, L. Kwok, and N. J. Long, “Lanthanides in magnetic resonance imaging,” Chemical Society Reviews, vol. 35, no. 6, pp. 557–571, 2006. View at: Publisher Site  Google Scholar
 N. H. de Leeuw, J. R. Bowe, and J. A. L. Rabone, “A computational investigation of stoichiometric and calciumdeficient oxy and hydroxyapatites,” Faraday Discussions, vol. 134, pp. 195–214, 2007. View at: Publisher Site  Google Scholar
 R. Snyders, D. Music, D. Sigumonrong, B. Schelnberger, J. Jensen, and J. M. Schneider, “Experimental and ab initio study of the mechanical properties of hydroxyapatite,” Applied Physics Letters, vol. 90, no. 19, Article ID 193902, 2007. View at: Publisher Site  Google Scholar
 A. Pedone, M. Corno, B. Civalleri et al., “An ab initio parameterized interatomic force field for hydroxyapatite,” Journal of Materials Chemistry, vol. 17, no. 20, pp. 2061–2068, 2007. View at: Publisher Site  Google Scholar
 S. Cotton, Lanthanide and Actinide Chemistry, John Wiley & Sons, Chichester, UK, 2006.
 M. V. S. dos Rezende, P. J. Montes, M. E. G. Valerio, and R. A. Jackson, “The optical properties of Eu^{3+} doped BaAl_{2}O_{4}: a computational and spectroscopic study,” Optical Materials, vol. 34, no. 8, pp. 1434–1439, 2012. View at: Publisher Site  Google Scholar
 B. G. Dick and A. W. Overhauser, “Theory of the Dielectric Constants of Alkali Halide Crystals,” Physical Review, vol. 112, no. 1, pp. 90–103, 1958. View at: Publisher Site  Google Scholar
 D. Mkhonto and N. H. de Leeuw, “A computer modelling study of the effect of water on the surface structure and morphology of fluorapatite: Introducing a Ca_{10}(PO_{4})_{6}F_{2} potential model,” Journal of Materials Chemistry, vol. 12, no. 9, pp. 2633–2642, 2002. View at: Publisher Site  Google Scholar
 Z. Du and N. H. de Leeuw, “A combined density functional theory and interatomic potentialbased simulation study of the hydration of nanoparticulate silicate surfaces,” Surface Science, vol. 554, no. 23, pp. 193–210, 2004. View at: Publisher Site  Google Scholar
 N. H. de Leeuw and T. G. Cooper, “The layering effect of water on the structure of scheelite,” Physical Chemistry Chemical Physics, vol. 5, no. 3, pp. 433–436, 2003. View at: Publisher Site  Google Scholar
 T. S. Bush, J. D. Gale, C. R. A. Catlow, and P. D. Battle, “Selfconsistent interatomic potentials for the simulation of binary and ternary oxides,” Journal of Materials Chemistry, vol. 4, pp. 831–837, 1994. View at: Publisher Site  Google Scholar
 A. Pedone, G. Malavasi, M. C. Menziani, A. N. Cormack, and U. Segre, “A new selfconsistent empirical interatomic potential model for oxides, silicates, and silicasbased glasses,” The Journal of Physical Chemistry B, vol. 110, no. 24, pp. 11780–11795, 2006. View at: Publisher Site  Google Scholar
 N. Y. Mostafa and P. W. Brown, “Computer simulation of stoichiometric hydroxyapatite. Structure and substitutions,” Journal of Physics and Chemistry of Solids, vol. 68, no. 3, pp. 431–437, 2007. View at: Google Scholar
 J. D. Gale, “GULP: a computer program for the symmetryadapted simulation of solids,” Journal of the Chemical Society, Faraday Transactions, vol. 93, no. 4, pp. 629–637, 1997. View at: Publisher Site  Google Scholar
 N. F. Mott and M. J. Littleton, “Conduction in polar crystals. I. Electrolytic conduction in solid salts,” Transactions of the Faraday Society, vol. 34, pp. 485–499, 1938. View at: Publisher Site  Google Scholar
 W. T. Lee, M. T. Dove, and E. K. H. Salje, “Surface relaxations in hydroxyapatite,” Journal of Physics Condensed Matter, vol. 12, no. 48, pp. 9829–9841, 2000. View at: Publisher Site  Google Scholar
 J. A. L. Rabone and N. H. de Leeuw, “Interatomic potential models for natural apatite crystals: incorporating strontium and the lanthanides,” Journal of Computational Chemistry, vol. 27, no. 2, pp. 253–266, 2006. View at: Publisher Site  Google Scholar
 S. Hauptmann, H. Dufner, J. Brickmann, S. M. Kast, and R. S. Berry, “Potential energy function for apatites,” Physical Chemistry Chemical Physics, vol. 5, no. 3, pp. 635–639, 2003. View at: Publisher Site  Google Scholar
 N. H. De Leeuw, “Local ordering of hydroxy groups in hydroxyapatite,” Chemical Communications, no. 17, pp. 1646–1647, 2001. View at: Google Scholar
 J. M. Hughes, M. Calmeron, and K. D. Crowley, “Structural variations in natural F, OH and Cl apatites,” American Mineralogist, vol. 74, pp. 870–876, 1989. View at: Google Scholar
 J. L. Katz and K. Ukraincik, “On the anisotropic elastic properties of hydroxyapatite,” Journal of Biomechanics, vol. 4, no. 3, pp. 221–227, 1971. View at: Publisher Site  Google Scholar
 J. C. Elliott, Structure and Chemistry of the Apatite and Other Calcium Orthophosphates, Elsevier, Amsterdam, The Netherlands, 1994.
 R. M. Wilson, J. C. Elliott, S. E. P. Dowker, and L. M. RodriguezLorenzo, “Rietveld refinements and spectroscopic studies of the structure of Cadeficient apatite,” Biomaterials, vol. 26, no. 11, pp. 1317–1327, 2005. View at: Publisher Site  Google Scholar
 L. M. RodriguezLorenzo, “Studies on calcium deficient apatites structure by means of MASNMR spectroscopy,” Journal of Materials Science: Materials in Medicine, vol. 16, no. 5, pp. 393–398, 2005. View at: Google Scholar
 P. Martin, G. Carlot, A. Chevarier, C. DenAuwer, and G. Panczer, “Mechanisms involved in thermal diffusion of rare earth elements in apatite,” Journal of Nuclear Materials, vol. 275, no. 3, pp. 268–276, 1999. View at: Publisher Site  Google Scholar
 M. Wakamura, K. Kandori, and T. Ishikawa, “Surface composition of calcium hydroxyapatite modified with metal ions,” Colloids and Surfaces A: Physicochemical and Engineering Aspects, vol. 142, no. 1, pp. 107–116, 1998. View at: Publisher Site  Google Scholar
 R. D. Shannon, “Revised effective ionic radii and systematic studies of interatomic distances in halides and chalcogenides,” Acta Crystallographica Section A, vol. 32, part 5, pp. 751–767, 1976. View at: Publisher Site  Google Scholar
 R. Jagannathan and M. Kottaisamy, “Eu^{3+} luminescence: a spectral probe in M_{5}(PO_{4})_{3}X apatites (M=Ca or Sr; X=F^{−}, Cl^{−}, Br^{−} or OH^{−}),” Journal of Physics: Condensed Matter, vol. 7, no. 44, pp. 8453–8466, 1995. View at: Publisher Site  Google Scholar
 I. Mayer, J. D. Layani, A. Givan, M. Gaft, and P. Blanc, “La ions in precipitated hydroxyapatites,” Journal of Inorganic Biochemistry, vol. 73, no. 4, pp. 221–226, 1999. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Ricardo D. S. Santos and Marcos V. dos S. Rezende. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.