Structural formulas of apatites: M(1)4M(2)6(XO4)6Z2 .
The structure of apatites allows large variations of composition given by the generic formula (M10(XO4)6Z2) including partial or complete substitution of both the cationic as well as the anionic sites, formation of nonstoichiometric forms and solid solutions. More than half of naturally occurring elements can be accommodated by apatite structure in a significant extend. The sixth chapter of this book is divided into five sections. The first, second and the third part deals with many examples of substitution including cationic substitution of M sites, anionic substitution of X-site and anionic substitution of Z-site, respectively. The remaining two sections continue with solid solution of apatites and ends with description of trace elements and their isotopes in the structure of apatite, respectively.
- Solid solution
- Trace elements
- Complexation of Metal Cations
Cation substitution, where M = Ca2+, Pb2+, Sr2+, Mg2+, Fe2+, Mn2+, Cd2+, Ba2+, Co2+, Ni2+, Cu2+, Zn2+, Sn2+, Eu2+, Na+, K+, Li+, Rb+, NH4+, La3+, Ce3+, Sm3+, Eu3+, Y3+, Cr3+, Th4+, U4+, U6+ and □.
Substitution for phosphorus by one or two cations, where X = PO43−, AsO43−, SiO43−, VO43−, CrO43−, CrO42−, MnO43−, SO42−, SeO42−, BeF42−, GeO44−, ReO53−, SbO3F4−, SiO3N5−, BO45−, BO33− and CO32−.
Z-site substation, where Z = F−, OH−, Cl−, O2−, O3−, NCO−, BO2−, Br−, I−, NO2−, CO32−, O22−, O2−, S2−, NCN2−, NO22− and □.
where □ represents the vacancy cluster .
Besides the monoionic substitution, the co-substitution and mutual combinations of substitutions in anionic and cationic sites (multi-ionic substitution) were also often reported ,,,. Mutual substitutions of trace elements into apatite structure brought new physicochemical, mechanical and biological properties in comparison with pure apatite or monoionic substituted apatite materials, e.g. hydroxylapatite .
Some substitutions can proceed only at the synthesis stage, while a limited ion exchange between solid apatite and surrounding solution can also occur. Due to their high chemical diversity and ion-exchange capabilities, apatites are considered as materials for toxic waste storage and for wastewater purification. The ion exchange in apatitic structures in human organism also presents an interest for medicine .
Recent studies have shown that a number of alkaline-earth-rare-earth silicates and germanates also have the apatite structure, and these have the cell sizes which span the division between the “apatites” and the “pyromorphites”. Some, particularly barium and lanthanum apatites, have the lattice parameters comparable with the members of the pyromorphite group. Thus, Ba2La8(SiO4)6O2 has the cell parameters a = 9.76 Å and c = 7.30 Å and Pb10(PO4)6F2 shows a = 9.76 and c = 7.29 Å, while Ba3La7(GeO4)6O1.5 has a = 9.99 Å and c = 7.39 Å and Pb10(AsO4)6F2 has a = 10.07 Å and c = 7.42 Å. During synthetic studies, however, it became apparent that the prediction of the composition of compounds with apatite-type structures could not be made solely on the basis of satisfying the valence considerations, since the occurrence of the apatite-type structure also appears to be determined by the ratio of the mean size of “A” ions (i.e. Ca ions in fluorapatite) to the mean size of “X” ions in XO4 ,.
The structure of hydroxyapatite allows large variations from its theoretical composition as well as the formation of nonstoichiometric forms and ionic substitutions. More than half of naturally occurring elements are known to be accommodated in the apatite lattice to significant extent. Ca2+ cation can be substituted by Na+, K+, Mg2+, Sr2+, Pb2+, Mn2+ (Fig. 1(a)) or rare-earth elements - (REE, Fig. 1(b)) and PO43− anions by AsO43−, SO42− or CO32− without destroying the apatite structure. The changes in lattice parameters must be indicative of the type of substitution occurring. For example, Cl− interchange for OH− ions causes a change to lattice parameters from a = 9.4214 Å and c = 6.8814 Å to 6.628 Å and 6.764 Å, respectively. Another example is the Sr2+ substitution for Ca2+, which causes lengthening of a- and c-axes from 9.418 Å and 6.884 Å to 9.76 Å and 7.27 Å, respectively ,,,,.
The substitutions at Z site play a very important role in the crystallography of specific species. The Z site lies in the channel formed by the X sites in fluorapatite and is of just the right size to fit between X atoms, and it lies on (001) mirror planes to yield the space group P63/
Some of the various families of substitutions that were experimentally established in apatites are summarized in Table 1. In general, the ions that substitute for Ca in the A position have the valences from 1 to 3 and the coordination numbers of VII at Ca(2) (6h) site and IX at Ca(1) (4f) site. Table 2 introduces the cation radii of possible apatite substituents at M-site .
|□2La2Ca4La2||□2||M-, Z-site vacancy|
|□2La2La2||Ge6||(OH)2||M-site vacancies (?)|
|Arens||Shannon and Prewitt|
|M+||K||1.33||1.38||1.46 (?)||1.51 (?)||1.55 (?, E)|
|Na||0.97||1.02||1.13 (?)||1.16||1.32 (?, E)|
An example of charge compensating substitution for phosphorus by two cations is the substitution during the synthesis of apatite species of the composition of Ca10(SiO4)3(SO4)3F2 (CSSF, fluorellestadite ) ,,:
These synthetic phases have mineral equivalents in the minerals from the ellestadite group, which are listed in Table 3. Since the mineral with ideal end-member formula Ca5(SiO4)1.5(SO4)1.5Cl is assumed not to exist, the name ellestadite-(Cl) is discredited .
|a [Å]||Lattice constant of hexagonal unit cell|
|c:a||Variable axial ratio|
|rMI [Å]||Shannon’s ionic radii of M(I)-site ion (nine-coordination)|
|rX [Å]||Shannon’s ionic radii of X-site ion|
|rMII [Å]||Shannon’s ionic radii of M(II)-site ion (seven-coordination for X = F−, eight-coordination for X = Cl− and Br−)|
|RZ [Å]||Shannon’s ionic radii of Z-site ion|
|Av CR [Å]||Average crystal radius = [(rM(I)x4) + (rM(II)x6) + (rXx6) + (rOx24) + (rZx2))]/42|
|MEN - OEN||Electronegativity difference between M and O atom|
|XEN - OEN||Electronegativity difference between X and O atom|
|MEN - ZEN||Electronegativity difference between M atom at M(II) site and Z atom|
|MEN - XEN||Electronegativity difference between M atom at M(I) site and X atom|
|M(I) – O(1) [Å]||Distance between M(I) atom and O(1) atom|
|M(I) – O(1)M(I)z = 0 [Å]||Distance between M(I) atom and O(1) atom with the constraint z = 0 at M(I)|
|ΔM(I)−O [Å]||Difference in the lengths M(I) – O(1) and M(I) – O(2).|
|Difference in the lengths M(I) – O(1) and M(I) – O(2) with the constraint z = 0 at M(I)|
|ΨM(I)−O [°]||Angle that M(I)-O(1) bond makes with respect to c|
|[°]||Angle that M(I)-O(1) bond makes with respect to c with the constraint z = 0 at M(I)|
|δM(I) [°]||Counter-rotation angle of M(I)O6 structural unit|
|ϕ M(I) [°]||Metaprism twist angle (π/3 – 2δM(I))|
|αM(I) [°]||Orientation of M(I)6 unit with respect to a|
|<X-O> [Å]||Average X-O bond length|
|<τO−X−O> [°]||Average O-X-O bond-bending angle|
|ρM(II) [Å]||M(II)-M(II) triangular side length|
|M(II) – Z [Å]||Orientation of M(II)-M(II)-M(II) triangles with respect to a|
|M(II) – O(3) [Å]||Distance between M(II) and O3 atoms|
|φO(3)−M(II)−O(3) [°]||O(3) – M(II) – O(3) angle|
|Etotal [eV]||Total energy calculated from ab initio calculations|
The fluorellestadite apatite and its solid solutions are minor components of many fluorine-mineralized clinkers. It is stable to liquidus temperature of 1240°C at which it incongruently melts to dicalcium silicate (2CaO·SiO2) and liquid . The solid-state synthesis and the luminescence properties of europium-doped fluorellestadite (CSSF:Eu2+) cyan-emitting phosphor were described by
The general composition of silico-sulfate apatites, i.e. ellestadites, is Ca10(SiO4)3−x/2(PO4)x(SO4)3−x/2(F,Cl,OH)2 and their structures conform to P63/
The syntheses of Sr and Pb analogues of CSSF are also reported . Strontium silico-sulfate apatite is not stable and decomposes to the mixture of strontium silicate and sulfate when heated to 1130°C for 30 min. Since high temperatures must be avoided, several attempts to prepare cadmium and barium silico-sulfate apatites were unsuccessful and the silicocarnotite-like phase was obtained from a mixture of the composition of Ca10(GeO4)3(SO4)8F2 rather than apatite .
Since there is a huge potential for the substitution in apatite structure (M(1)4M(2)6(XO4)6Z2 and for the formation of solid solution as well, the classification method enables to identify the key crystallographic parameters which can serve as strong classifiers of crystal chemistries. The structure maps for apatite compounds via data mining were reported by
Basically, the structure map approach involves the visualization of the data of known compounds with known crystal structures in a two-dimensional space using two scalar descriptors (normally heuristically chosen), which are associated with physical/chemical properties, crystal chemistry or electronic structure. The objective is to map out the relative geometric position of each structure type from which one tries to discern qualitatively if there are strong associations of certain structure types to certain bivariate combinations of parameters .
1.The αMII (rotation angle of MII-MII-MII triangular units);
The (angle that M(I) – O1 bond makes with the c axis when z = 0 for the M(I) site).
That enables to classify the apatite crystal chemistries based on the site occupancy at M, X and Z sites and this classification is accomplished using the K-means clustering analysis (Fig. 3).
For example, clusters 1 and 2 (k = 1 and k = 2) correspond to F-apatites (fluorapatites). They are well localized in the structure map and are characterized by relatively low αMII and . Two F-apatites which do not belong to the clusters k = 1 and k = 2 are Hg5(PO4)3F (in k = 4) and Zn5(PO4)3F (in k = 7). While the existence of fully stoichiometric Zn5(PO4)3F apatite compound is uncertain due to relatively small ionic size of Zn2+ cations, the relative position of Hg5(PO4)3F suggests some peculiar characteristics .
|k = 1 and k = 2||Site||M||Ba, Pb, Sr, Ca|
|X||P, V, Mn|
|k = 3||Site||M||Ba|
|k = 4||Site||M||Sr, Hg|
|k = 5||Site||M||Ca, Cd, Pb|
|X||V, Cr, As|
|k = 6||Site||M||Ca, Pb|
|k = 7||Site||M||Zn|
Even though Ca2+ and Hg2+ cations have roughly the same ionic size (1.18 and 1.23 Å at M(I) site), their electronegativity data indicates that Hg atoms (electronegativity value of 2 in Pauling scale) are relatively highly covalent compared to Ca atoms (electronegativity value of 1 in Pauling scale). In the structure map, this covalent character is predicted to be manifested in the bond distortion angle .
6.1. Cationic substitution at M(1) and M(2) sites
6.1.1. Strontium-substituted apatites
Sr2+ ion, which is larger than Ca2+, is ordered almost completely into the smaller Ca(2) site in the apatite structure (Fig. 4). The bond valence sums of Sr ions at two sites demonstrate that Sr is severely overbonded at apatite Ca site but less at Ca(2) site. Complete ordering of Sr into Ca(2) sites has important implications for the diffusion of that element in the apatite structure. It is the subject of several recent studies. The diffusion of Sr in (001) was shown to be as rapid or even more rapid than the diffusion parallel to . As there are neither sites available for Sr, which are linked in (001), nor any interstitial sites, which can contain Sr2+ ion, the diffusion mechanism involving the vacancies or defects or both is indicated .
A series of Sr-substituted hydroxyapatites, (SrxCa1−x)5(PO4)3OH, where x = 0.00, 0.25, 0.50, 0.75 and 1.00, was investigated by
Strontium is often substituted for calcium in order to confer the radio-opacity in glasses used for dental cements, biocomposites and bioglass-ceramics. It can be concluded that strontium substitutes for calcium with little change in the glass structure as a result of their similar charge to size ratio. Glasses with low content of strontium nucleate in the bulk to form calcium apatite phase. Glasses with medium strontium content nucleate to mixed calcium-strontium apatite at the surface and glass fully substituted by strontium to strontium fluorapatite .
6.1.2. Magnesium-substituted apatite
Magnesium-substituted hydroxyapatite (MgHAP) powders with different crystallinity levels, prepared at room temperature via a heterogeneous reaction between Mg(OH)2/Ca(OH)2 powders and (NH4)2HPO4 solution using the mechanochemical- hydrothermal route, were reported by
Two effects of different magnesium sources (magnesium nitrate and magnesium stearate) on the synthesis of Mg-substituted hydroxyapatite (Mg-n-HAP) nanoparticles by the co-precipitation method were investigated by
6.1.3. Copper-substituted apatites
Copper-substituted hydroxyapatite (Ca10−xCux(PO4)6(OH)2 (where x = 0.05 – 2.0) and fluorapatite Ca10−xCux(PO4)6F2 (x = 0.05 – 2.0) were synthesized by
6.1.4. Nickel-substituted apatites
In this process, Ni2+ ions are first adsorbed onto the surface of hydroxyapatite (surface complexation, Section 6.5.2) and then the substitution of Ca2+ for Ni2+ ions takes place.
6.1.5. Zinc-substituted apatites
Zinc is a common bioelement. The zinc content in human bones ranges from 0.0126% to 0.0217% by weight . Zinc as a cationic substituent in hydroxyapatite provides the option to counteract the effects of osteoporosis . The incorporation of zinc into the HAP structure (Zn-HAP) was abundantly studied, owing to the key effect of Zn2+ cations in several metabolic processes that makes zinc eligible for use in many biomedical applications and to its possible antimicrobial activity .
The results of structure analysis indicated that Zn ions substituted partially for Ca ions in the apatite structure and the upper limit of Zn substitution for Ca in HA was about 20 mol.%. In general, the HAP lattice parameters, a and c, decreased with Zn addition .
Zn-substituted apatite was synthesized by the precipitation method as follows :
where 0 ≤ x ≤ 1. The pH of the solution was adjusted to 8 by aqueous solution of NH3, and the reaction mixture was kept at 90°C for 5 h with stirring. The resulting suspension was then subjected to suction filtration, and the powdery product was dried at 100°C for 10 h. It is known that the usage of chloride or nitrate of calcium as a starting reagent may cause the incorporation of Cl− or NO3− into the structure of apatite. This can be avoided by the utilization of acetate salts, because acetate ions are not incorporated into the apatite, i.e. they would not affect the apatite structure.
6.1.6. Iron-substituted apatites
The synthesis and the characterization of iron-substituted hydroxyapatite via a simple ion-exchange procedure were described by
6.1.7. Cobalt-substituted apatites
The synthesis and the characterization of cobalt-substituted hydroxyapatite (Co-HAP) powders via the precipitation method were described by
Single crystals of chlorapatite [Ca5(PO4)3Cl] with the substitution of approximately 20% of Ca2+ by Co2+ (space group P63/
6.1.8. Manganese-substituted apatites
Naturally occurring manganese-substituted apatite is known as manganese-bearing apatite (Mn,Ca5(PO4)3F, Section 1.1). According to the findings of
It is interesting to note that apatite acts effectively as a geochemical sieve that traps Mn2+ and excludes Fe2+ elements, which are virtually inseparable in most geochemical systems. The bond valence sums for Fe2+ at apatite Ca sites yield 1.26 and 1.19 valence units for Ca(1) and Ca(2) sites, respectively; large discrepancy from the formal valence prohibits extensive substitution of Fe2+ in the apatite structure (Fig. 5) .
The crystal structure of pale blue transparent Mn-rich fluorapatite (9.79 wt.% of MnO) with optimized formula (Ca8.56M2+1.41Fe2+0.01)P6O24F2 was resolved by
6.1.9. Substitution of REEs in apatite
Crystals of La-, Gd- and Dy-bearing fluorapatite [La-FAP, Gd-FAP, Dy-FAP, Ca10−x−2y NayREEx+y(P1−ySixO4)Z2, where x = 0.24 – 0.29 and y = 0.32 – 0.36] were synthesized by hydrothermal route by
The partitioning of REE between two Ca positions in apatite contradicts usual first-order dependence on spatial accommodation, with LREE - ,,,,, in particular, favoring the smaller Ca(2) position. This behavior was variously ascribed to the control via :
The preference of individual REE among multiple Ca positions in minerals (the site occupancy of individual REE) was not extensively studied because of the inability of conventional diffraction methods to distinguish among individual elements at multiply-occupied sites. The site preference for individual LREE from theoretical bond-valence sums was estimated by
The isomorphic substitutions of neodymium for strontium in the structure of synthetic Sr5(VO4)3OH apatite structure type (P63/
where x = 0, 0.02, 0.08, 0.10, 0.12, 0.14, 0.16, 0.18 and 0.20. The substitution scheme can be expressed as:
The procedure includes three stages:
Preparation of solution;
Treatment of the dry residue.
The solutions for the thermolysis were prepared by dissolving Sr(NO3)2 in water; Nd2O3 was dissolved in water with nitric acid added; NH4VO3 was dissolved in water with hydrogen peroxide added. Dry residues after concentrating the solutions were pestled in an agate mortar and calcined with the temperature steadily raised from 600 to 800°C and intermittent grindings .
6.1.10. Actinides-bearing apatites
Thorium and uranium (actinides , -)-bearing apatites were synthesized by
The structure refinements of U-doped chlorapatites show that U is essentially distributed equally between the two Ca sites with UCa(2)/UCa(1) values, which range from 0.89 to 1.17. The results of Th-doped chlorapatites show that Th substitutes into both Ca(1) and Ca(2) sites with ThCa(2)/ThCa(1) values, which range from 0.61 to 0.67. In Th-doped strontium apatites with F and Cl end-members, Th is incorporated into both Ca(1) and Ca(2) sites. The range of ThCa(2)/ThCa(1) values is 0.56 to 1.00 for the F end-member and 0.39 to 0.94 for the Cl end-member. U-doped samples indicate that U in fluorapatite is tetravalent, whereas, in chlorapatite, it is heterovalent but dominantly hexavalent .
Based on the chemical analyses of U-, Th-doped fluor-, chlor- and strontiumapatite specimens in this study, local charge compensation may be maintained by the following coupled substitutions (M represents U or Th and  represents the vacancy) :
The incorporation of U and Th into fluorapatite results in a decrease of the size of both Ca polyhedra, but the incorporation of U and Th into chlorapatite gives rise to an increase in the volume of both Ca polyhedra. The decrease of both Ca polyhedral volumes in fluorapatite caused by the substitution of U and Th can be explained by the decrease of ionic radius from Ca to U and Th. However, the increase in the volume of both Ca polyhedra in chlorapatite is hard to understand. Because of the effect on Ca(2) polyhedron caused by the replacement of F− by Cl−, it can be explained by the structural distortion of Ca(2) polyhedron .
Uranium-doped oxy-silicophosphates (britholites) of the composition of CaxLay(SiO4)6−u (PO4)uOt:U4+ were synthesized by
where An4+ substitutes for tetravalent U4+ and Th4+.
The incorporation of thorium in the structure is probably possible due to small differences of ionic radius between calcium (1.06 Å), neodymium (1.05 Å) and thorium (1.00 Å). In order to ensure the quantitative incorporation of thorium, it appeared necessary to consider the coupled substitution .
instead of the substitution scheme:
Indeed, in the first way, homogeneous and single-phase solid solutions were prepared from Ca9Nd(PO4)5(SiO4)F2 to Ca9Th(PO4)4(SiO4)2F2 leading to full neodymium substitution. Associated small increase of the unit cell parameters results from the simultaneous replacement of phosphate groups by bigger silicate. It was accompanied by a significant change in the grain morphology. These results contrast with those obtained when the coupled substitution according to Eq. 19 was performed, which confirmed the limitation of about 10 wt.% in the Th substitution . Good resistance of the materials to influence of aqueous solutions enables their utilization for the immobilization of tetravalent actinides in phosphate ceramics .
6.1.11. Silver-substituted apatite
The favorable biocompatibility of hydroxyapatite (HA) makes it a popular bone graft material as well as a coating layer on metallic implant. One common and accepted strategy to prevent the implant-related infections is to create antibacterial properties for the implant. Silver ions can be either incorporated into the apatite during the co-precipitation process (AgHAP-CP) or subjected to the ion exchange with calcium ions in apatite (AgHAP-IE). The incorporation of silver ions into apatite is based on the equation :
where y is the molar amount of silver to be incorporated. However, the distribution of silver ions in AgHAP-CP and AgHAP-IE was different, thus affecting the antibacterial action.
6.1.12. Cadmium-substituted apatites
The absorption of cadmium cations in apatites is relevant both from the medical standpoint of cadmium uptake by human bones, as well as since cadmium migration in nature involves the absorption and desorption equilibria with natural minerals, including apatites. Cadmium has a slight preference for Ca(I) site in fluorapatite and for Ca(II) site in hydroxyapatite ,. The interactions between these two ions (Cd and Ca) during absorption and ionic change processes in apatites present therefore considerable practical and theoretical interest. Cadmium is also a frequent heavy toxic pollutant element in water .
Calculated energy differences (E) between these sites are of 12 and 8 kJ·mol−1 for fluorapatite and hydroxylapatite, respectively. The preference is not strong, and however, a part of the sites of the other type is also occupied by cadmium ions. The relative site occupation can be expressed by the equation :
where E = E(Cd2+ or Zn2+ on Ca(1)) – E(Cd2+ or Zn2+ on Ca(2)). At T = 298 K, P = 85 and 17 Cd2+ in fluorapatite and hydroxylapatite, respectively. From the value of P and from the fact that the sum of the two probabilities is 1, one can calculate that the probability of the lower-energy site occupancy is of 99% and 94%, respectively.
6.2. Anionic substitution at X-site
Pentavalent arsenic, vanadium and chromium substitution can completely replace phosphorus in calcium, strontium and barium fluor- and chlorapatites. Calcium fluor-vanadate, -arsenate and -chromate structures were distorted compared to normal hexagonal apatite. Manganese completely replaced phosphorus only in barium apatites, while chromium and manganese could not be incorporated into lead apatites. Excluding these exceptions, continuous solid solutions were formed between the phosphate and/or vanadate and the chromate or manganese analogues for given divalent and halide ions . The substitution of CO32− ions at X- (carbonate-apatite of A-type) and Z-site (carbonate-apatite of B-type) was already described in Section 4.6.
6.2.1. Arsenate substitution in hydroxylapatite
The arsenate (As5+) substitution in the hydroxylapatite structure was examined by
Complete PO43− ↔ AsO43− substitution was also recognized in experimental studies of apatite analogues, such as in the system Sr5(PO4)3OH-Sr5(AsO4)3OH . The Rietveld refinement of Sr5(AsO4)3Cl (pentastrontium tris[arsenate(V)] chloride, 890.31 g·mol−1) from high-resolution synchrotron data was performed by
6.2.2. Vanadate substitution in hydroxylapatite
The synthesis of synthetic alkaline-earth vanadate hydroxylapatites from hydroxide fluxes was performed by
The compounds (solid solution) of the composition of Pb5(PxV1−xO4)3Cl (0 ≤ x ≤ 1), which are synthetic analogues of minerals pyromorphite, vanadinite and endlichite, were synthesized for the first time by
The variations of unit cell parameters as a function of composition respect Vegard’s law. These compounds are structurally built of discrete phosphate or vanadate tetrahedra linked to one another by lead polyhedra, which form joint layers (Fig. 10). Apatite-type structures offer typically two crystallographic positions for cations differing in the coordination number and local symmetry. Lead atoms occupying the first positions form polyhedra shaped as three-capped trigonal prisms of PbO9 having the symmetry C3, the columns of which run along the threefold axis. Distorted two-capped trigonal prisms of PbO6Cl2 residing in the second positions have local symmetry С1.
6.2.3. Chromium analogues of apatite
Ca5(CrO4)3OH, which is isomorphous to hydroxyapatite (Section 1.5.2): space group P63/
m, a = 9.683 Å and c = 7.010 Å, a:c = 1:0.7239, V = 569.20 Å3, Z = 2.
Sr5(CrO4)3OH with predicted lattice constants a = 9.9561 Å and c = 7.488 Å.
Ca5(CrO4)3F with predicted lattice constants a = 9.733 Å and c = 7.0065 Å
Ca5(CrO4)3Cl with predicted lattice constants a = 10.1288 Å and c = 6.7797 Å.
Sr5(CrO4)3F with predicted lattice constants a = 9.9349 Å and c = 7.5037Å.
Sr5(CrO4)3Cl with lattice constant a = 10.125 Å and c = 7.328Å.
Sr5(CrO4)3Br with predicted lattice constants a = 10.2895 Å and c = 7.2712 Å.
These compounds are in general prepared by the ignition of mixture of alkaline-earth carbonates, hydroxides or oxides with Cr2O3 in the presence of water vapor. Ca3(CrO4)2 compound (orthochromate), which is isomorphous with Ca3(PO4)2, is formed as an intermediate by carrying out the synthesis in dry atmosphere; this compound is often identified as 9CaO·4CrO3·Cr2O3 ,,,.
Theoretical compositions and formula weights of chromium apatite analogues are given in Table 5.
Sr10(CrO4)6F2 possesses typical hexagonal structure of apatite with the space group P63/
The crystal structure (Fig. 12) and the magnetic properties of strontium chromate phase (Sr5(CrO4)3(Cu0.586O)) with apatite-like structure were determined by
This phase crystallizes in the space group P63/
6.2.4. Selenium analogues of apatite
Selenium oxyanion-substituted hydroxyapatite (SeHAP) was synthesized as a promising material for the treatment of bone cancer to reduce the probability of recurrence, because selenium plays an important role in protein functions and it has significant effect on the induction of cancer cell apoptosis . Another study indicated that selenite (SeO42−) or selenite (SeO32−) oxyanions exert their cancer chemopreventive effects by direct oxidation of critical thiol-containing cellular substrates and that they are more efficacious anticarcinogenic agents than selenomethionine or selenomethylselenocysteine with a lack of oxidation capability ,.
Selenium was incorporated into the hydroxyapatite lattice by replacing some of the phosphate groups with selenite groups. SeO42− (selenate) ion has tetrahedral structure like PO43− ion (Fig. 5 and Table 4 in Section 1.2), but it is slightly larger (2.49 Å) in diameter than phosphate ion, which is 2.38 Å in diameter. By contrast, SeO32− (selenite) ion has very similar diameter (2.39 Å), but it has a quite different flat trigonal pyramid geometry. The substitution of bivalent selenium oxyanions forms positively charged vacancy compensated by simultaneous decalcification and dehydroxylation according to the reaction ,:
The SeHAP crystal lattice parameters increased slightly as the Se concentration increased when the Se/P ratios were less than 0.5 . All samples prepared via the precipitation method from aqueous solution by
6.3. Substitution at Z-site
6.3.1. Nitrogen-containing apatites
Nitrogen was incorporated into hydroxyapatite by dry ammonia treatments at temperatures between 900 and 1200°C in the presence of graphite. The process of synthesis of cyanamidapatite (Ca10(PO4)6CN2, Ca10(PO4)6NCN) can be described by the following chemical equations :
Ammonia reacts with graphite during the thermal treatment forming [CN2]2− ions (Eq. 24). These cyanamide ions interchange with moveable OH− ions situated on the sixfold screw axis of apatite to form cyanamidapatite (Eq. 25). A similar reaction is known for the synthesis of calcium cyamide from calcium oxide:
The treatments at temperatures above 1200°C or long-term treatments destroy the apatite lattice completely through the phosphate reduction. Cyanimide ions lose their sites in the apatite lattice and the nitrogen content decreases . The synthesis of Ca10(PO4)6CN2 apatite provides the evidence that the hydroxylapatite structure is able to incorporate larger organic molecules .
Direct transformation of TCP (Ca3(PO4)2) into cyanamidapatite according to the reaction:
was also proposed by
6.3.2. Peroxide-doped apatites
Although “oxygenated” apatites were not much investigated compared to other substituted apatites, some past studies have, however, reported the possibility of apatitic channels to incorporate oxygenated species such as H2O2 or O2 or molecular ions including O22− (the peroxide ion) and superoxide O2−. They are single-phase nanocrystalline apatites, where part of apatitic OH− ions are replaced by oxygenated species. Typically by peroxide ions (quantified) and at least the traces of superoxide ions can be prepared by the precipitation from aqueous calcium and phosphate solutions in the presence of H2O2 under medium room temperature ,.
The local structure of hydroxyl-peroxy apatite was described by
According to the concentration of peroxide ions in hydroxyl-peroxy apatite and the theoretical value, the corresponding formula for the hydroxyl-peroxy apatite is proposed as follows : Ca10(PO4)6(OH)1.34−2x(O2)0.33(O)x□0.33+x.
A scheme of possible configurations of hydroxyl ions, peroxide or oxide ions and vacancies in the channel along the crystallographic c-axis in hydroxyl-peroxy apatite is illustrated in Fig. 13. Peroxide ions incorporated into HAP are located in the channel of apatite structure through the substitution of a portion of OH− radicals, and the material is a solid solution of hydroxyl- and peroxide apatite.
where  was the vacancy. O2− ion was active and could react with O2 to produce O22−.
Peroxide ions associated with the vacancies were situated placed in the channel of HA lattice along the c-axis through the substitution of a portion of OH radicals. The molecular ions constituted a symmetric vibrator with a stretching vibration active in Raman spectrometry. This vibration was recorded at 750 cm−1 in the Raman spectra of O22−-containing HA samples. The final product was a solid solution of hydroxyl- and peroxide-apatite. However, the existence of peroxide ions in the HA lattice caused the contraction of the unit-cell dimensions of HA materials. In addition, a new hydrogen bond was formed between peroxide ions and adjacent OH radicals, which was determined by using molecular spectroscopy analysis. During annealing treatment in air, peroxide ions decomposed and the substituted OH radicals re-entered the HA lattice, resulting in the elimination of the structural aberrations caused by the incorporation of peroxide ions. The concentration of peroxide ions present in HA samples was measured by chemical analysis .
6.3.3. Chalcogenide phosphate apatites
Ca10(PO4)6S: a = 9.4619 Å, c = 9.8342 Å, c:a = 0.7223 and V = 529.88 Å3 (Fig. 14(a));
Sr10(PO4)6S: a = 9.8077 Å, c = 9.2089 Å, c:a = 0.7350 and V = 600.53 Å3 (Fig. 14 (b));
Ba10(PO4)6S: a = 10.2520 Å, c = 7.6590 Å, c:a = 0.7471 and V = 697.14 Å3 (Fig. 14(c));
Ca10(PO4)6Se: a = 9.5007 Å, c = 9.8406 Å, c:a = 0.7200 and V = 534.73 Å3 (Fig. 14 (d)).
were reported by
These four apatites are isostructural and crystallize in the trigonal space group P3 over bar with the chalcogenide ion positioned at (001/2). Sulfoapatites show no ability to absorb H2S in the way that oxyapatite absorbs H2O at elevated temperatures. This can be attributed to the position of sulfide ion and the way it influences the crystal structure around vacant chalcogenide position at (000) .
6.3.4. Metaborate ion-containing apatite phase
Strontium borate-phosphate Sr10(PO4)5.5(BO4)0.5(BO2) - was prepared from SrCO3, NH4H2PO4 and H3BO3 at high temperature (from 1150 to 1550°C) and was found to be free of alkali metal compounds. Sr10(PO4)5.5(BO4)0.5(BO2) phase is a derivative of the apatite crystal structure with metaborate ion at Z-site: space group P3, a = 9.7973 Å, c = 7.3056 Å, V = 607.29 Å3, Z = 1 ,,.
The strontium sites are found to be fully occupied, while [PO4]3ˉ tetrahedra are partially replaced by [BO4]5ˉ groups. The crystal structure contains Sr cations occupying the 6g (Sr(1)) and 2d (Sr(2), Sr(3)) sites, isolated tetrahedral [PO4]3ˉ/[BO4]5ˉ groups and linear [BO2]ˉ groups located in hexagonally shaped (trigonal antiprismatic) channels formed by Sr(1) atoms and running along  (Fig. 15). The space group of the present compound is reduced to P3, because the orientation of the [PO4]/[BO4] tetrahedra destroys the mirror plane characteristic for the apatite crystal structure (P63/
6.4. Solid solutions of apatites
Crystalline solid solutions -  of apatites are frequently encountered where the possibility and type depend on the condition of formation or preparation, thermal history after formation and the end-members of the series . The structure of ternary solid solution of hexagonal (P63/
The solid solution in hexagonal ternary apatite is achieved by a 0.4 Å shift along the c-axis of Cl atom relating to its position in end-member chlorapatite. This adjustment affects the Markovian sequence -  of anions in the (0,0,z) anion columns by providing a structural environment that includes column OH species at the distance of 2.96 Å from Cl. The shift of Cl atom is accompanied by splitting of Ca(2) atoms into two distinct positions as a function of the kind of anion neighbor (F or OH vs. Cl). Additional nonequivalent Cl site, similar to that in end-member chlorapatite, is also present. Those Cl atoms with adjacent OH occupy a site different from Cl atoms adjacent to vacancies in the anion column .
Reduction of symmetry in monoclinic ternary apatite results from the ordering of Cl and OH within the anion columns. The atomic positions of Cl and OH in the anion column are equivalent to those in hexagonal ternary apatite, but each is ordered into only one of the two hexagonal symmetry-equivalent sites .
The apatite supergroup minerals of the solid solution :
where x = 0 – 3 and y = 0 – 1.5 were found in altered calcareous xenoliths within the ignimbrite of the Upper Chegem caldera, Northern Caucasus, Russia. These minerals belonging to the apatite supergroup occur in all zones of skarn from the core to the contact with ignimbrite as follows: brucite-marble, spurrite (Ca5(SiO4)2CO3 ), humite (Mg7(SiO4)3(F,OH)2 ) and larnite - (Ca2SiO4 ,,,) zones. They are associated with both high-temperature minerals: reinhardbraunsite (Ca5(SiO4)2(OH)2 ), chegemite (Ca7(SiO4)3(OH)2 ), wadalite (Ca6Al5Si2O16Cl3 ), rondorfite (Ca8Mg(SiO4)4Cl2 ), cuspidine (Ca4(Si2O7)F2 ), lakargiite (CaZrO3 ) and srebrodolskite (Ca2Fe3+ 2O5 ), corresponding to the sanidinite metamorphic facies, - and secondary low-temperature minerals: calcium hydrosilicates (hillebrandite , awfillite , bultfonteinite ), hydrogarnets  and minerals of the ettringite group .
The minerals of the apatite supergroup often form elongated cracked hexagonal or pseudo-hexagonal crystals up to 250 μm in size as well as grain aggregates. A new solid-solution series was found between ellestadite and svabite-johnbaumite (±apatite) with the ellestadite type isomorphic substitution according the following scheme :
where R = As5+ and P5+. The As content in investigated minerals decreases from the contact skarn zone with the ignimbrite towards the core of altered xenoliths (from 2.11 As pfu - to 0), for example :
Si, S-bearing hydroxylapatite:
The crystals of As-bearing phases belonging to the investigated solid solution are heterogeneous and small in size. Therefore, X-ray single-crystal data were obtained for only Si, S, As-bearing hydroxylapatite (see the formula above): P63/
The hydrothermal synthesis of vanadate/phosphate hydroxyapatite solid solutions of the composition of Ca10(VO4)x(PO4)6−x(OH)2, where x = 0, 1, 2, 3, 4, 5 and 6, was firstly reported by
The crystal structure of 11 samples of synthetic Na-Ca-sulfate apatite systems of the composition of Na6.45Ca3.55(SO4)6(FxCl1−x)1.55, where x = 0 – 1, was refined by
Lead apatites form a family of isomorphous compounds, and well-known members of the group are mimetite (Pb5(AsO4)3Cl, Section 1.6.7) and pyromorphite (Pb5(PO4)3Cl, Section 1.6.4). Isostructural with vanadinite Pb5(VO4)3Cl, these three constituents form a ternary system within the apatite group of P63/
A number of compounds of the mimetite Pb5(AsO4)3Cl-pyromorphite Pb5(PO4)3Cl solid-solution series were synthesized at room temperature by
Solid solutions of Pb8M2(XO4)6 lead alkali apatites were studied by
Some other examples of apatite solid solutions are listed below :
Ca2Y8(SiO4)6O2 – Ca8Y2(PO4)6O2;
Ca2La8(SiO4)6O2 – Ca8La2(PO4)6O2;
Mg2Y8(SiO4)6O2 – Y10(SiO4)4(BO4)2O2;
Pb4+ 3Pb2+ 5Y2(SiO4)6O2 – Pb2+ 2Y8(SiO4)6O2;
Ca10(PO4)6(OH)2 – Ca4Y6(SiO4)6(OH)2;
M10(PO4)6F2 – M10(PO4)6F2 (M = Sr, Ba, Pb);
M10(PO4)6F2 – M10(MnO4)6F2 (M = Sr, Ba, Pb).
6.5. Trace elements and their isotopes
Since apatite is an important accessory mineral in most common rock types, it is often used in trace element and isotope investigations of igneous and metamorphic rocks . Stable isotope compositions of biologically precipitated apatite in bone, teeth and scales are widely used to obtain the information on the diet, behavior and physiology of extinct organisms and to reconstruct past climate in terrestrial and marine conditions .
Broad spectrum of substitutions in the apatite lattice allows the incorporation of various isotopes, which offer a number of instruments for the interpretation of paleoenvironment and diagenesis. The relative stability of francolite compared to other sedimentary minerals led to an enormous number of studies and applications. Various isotopes occupy the Ca2+ and PO43− sites in the lattice of apatite (Fig. 16).
Since the earliest application in deep-time study of Late Cretaceous paleotemperatures in 1950, the oxygen isotope paleothermometry is based on the temperature dependence of oxygen isotope fractionation between authigenic minerals -  and ambient waters. Under the equilibrium conditions, the 18O/16O ratio of sedimentary carbonate and phosphate minerals depends only on the temperature of precipitation and on the 18O/16O ratio of ambient water. Thermodynamic relationships and bond vibrational frequencies can be used to determine the mineral-water isotopic fractionation relations but not with the precision and accuracy necessary for the paleothermometry. Such an application requires the calibrations based on mineral-water oxygen exchange experiments at high temperatures, mineral precipitation experiments at low temperatures and/or natural experiments using minerals grown under known conditions ,.
Carbon, oxygen and sulfur isotopes are used to reconstruct the oxygenation stages of the sediments during organic matter degradation and precipitation of apatite. The application of this method gives good results for modern and Neogene deposits. In older occurrences, the signature of carbon and oxygen composition is commonly overprinted by diagenetic and burial diagenesis ,. The carbon isotope ratios of apatite can be used to interpret the source of carbon in magmas and metamorphic fluids using the assumption that the carbon isotope fractionations between phases are small in igneous and metamorphic systems .
The carbon isotope analysis of bioapatites was first applied to terrestrial mammals in early 1980s ,,. While it is now known that some bones do undergo the C-isotope exchange extremely readily, collagen, bone and enamel record different periods of time during the life of a single individual, and the diet may change. That means, there is a fundamental ambiguity (preservation vs. normal intraindividual differences) in interpreting the isotopic differences among different tissues. Unfortunately, early results were taken to imply that all bioapatites are unreliable, and it was not until the 1990s that it became accepted that the tooth enamel, at least, is a robust recorder of diet. Thus, the early work of
Fossil biogenic apatites display the trace element compositions that can record environmental and biological signals, give insights into past water compositions or be used for dating paleontological and archeological bones and teeth. The mechanisms of the process of trace element and their isotopes incorporation into apatites of skeletal phosphatic tissues are described by
Partitioning between aqueous fluids and crystals;
Surface adsorption, complexation and chelation;
6.5.1. Partitioning of divalent cations
Partitioning of divalent cations is defined by the chemical equilibrium expressing the divalent cation (Me2+) exchange between apatite and aqueous solutions :
where (aq) and (ss) refer to the aqueous solution and to the solid solution, respectively. The equilibrium constant of Eq. 31 can be written as
where x is the molar fraction in apatite solid solution, m is the molality in water, λ is the activity coefficient of the component in the solid solution, K(T) is the solubility product of the end-member at temperature T and γ is the ion activity in the aqueous solution, the ratio of which in water is assumed to be equal to one. The activity coefficients in regular solid-solution model are described by Margules parameters -  and can be approximated by the elastic energy due to the deformation of the host crystal lattice around the substituted cation ,:
where NA is the Avogadro number, E is the Young’s modulus of the crystal, ri is the ionic radius of cation normally occupying the site in the i-compound (Ca in apatite) and rj is the ionic radius of the substituted cation in compound j. The elasticity of hydroxyapatite gives E = 114 ± 2 GPa.
where the term exp (-ΔGideal/RT) is the free enthalpy change of the reaction 31, equivalent to the ratio of end-member solubility products. Unlike carbonates, the data of solubility products and thermodynamics for end-member apatites are scarce. When no data are available for the solubility and enthalpy of formation of the end-members, it is assumed that the elastic energy term dominates over the partitioning, i.e. ΔGideal ≪ WGMeCa. Promising ways for obtaining the enthalpies of formation and the substitution energies are the first-principle calculations  and the atomistic modeling ,.
For heterovalent substitutions, the equilibrium reaction becomes complex since complementary substitutions are necessary to maintain the charge balance in the crystal. Typically, the substitution of trivalent elements of important rare-earth series requires the compensation by Na+ for Ca2+ at an adjacent site or even more complex substitution scheme involving carbonate groups or fluorine. In that case, most thermodynamic data required for the calculation of the equilibrium constant are not available. Among the series of elements with the same charge and substitution scheme, the pattern of equilibrium constants, or of distribution coefficients, can be approximated by combining Eqs. 33 and 34 :
where Eeff is the effective Young’s modulus and r0 is the optimum radius for maximum equilibrium constant KD0, all of which will depend on the charge of the considered series of elements. These parameters can be adjusted to experimental data such as partition coefficients between minerals and liquids and lead to parabola-like curves, the position and curvature of which depend on the charge of the element. This approach was applied so far only to rare-earth elements in apatite, where the relative partition coefficients were extrapolated from magmatic temperatures around 800°C to low temperatures appropriate for fossil diagenesis .
6.5.2. Complexation of metal cations
The complexation of metal cations in aqueous fluids involves binding with a broad range of molecules from simple inorganic ones (e.g. carbonates, phosphates and sulfates) to complex organic ones (humic acids, amino acids, proteins, enzymes, etc.). For molecules with several bonding sites and structural flexibility (e.g. multidentate or chelator), the complexation is thermodynamically favored with respect to the complexation with several monodentates having one bonding site; the process is named chelation. Chelators can be adsorbed on mineral surfaces while remaining complexed to metallic cations. The pattern of the partition coefficients associated with this process was measured for rare-earth elements complexed with humic acids and manganese oxides. It shows null fractionation along the whole series; the effect of chelation is therefore to screen the trace element in the crystal or ligand field and to suppress the fractionation associated with ionic radius variations and tetrad effects, and most of the anomalies associated with redox of Ce . Similar effects might occur for the adsorption of chelated metals on other mineral surfaces and in particular phosphates. In addition to chelators, the transition metals also form complexes with proteins and enzymes that interact with bones and teeth in living organisms and may influence their incorporation in bioapatite .
6.5.3. Diffusion processes
Solid-state diffusion in crystals is a thermally activated process governed by the enthalpy of formation and of migration of defects and usually well described by the Arrhenius relation :
where D0 is the pre-exponential factor corresponding to the diffusion coefficient at infinite temperature and ΔHa is the activation enthalpy (or energy) of the diffusion process. The extrapolation of high-temperature diffusion data of trace elements in apatite shows that these processes are inefficient at temperatures below 300°C, which cover the conditions of diagenetic alteration up to low-grade metamorphism ,.
The differing initial and boundary conditions imposed in three sets of diffusion experiments:
In-diffusion with powder sources experiment;
consequently resulting in different solutions to the diffusion equation. However, in all cases, the process can be described as one-dimensional, concentration-independent diffusion .
A plot of diffusivities of various cations and anions in apatite is shown in Fig. 17(a). The diffusivities of Mn are similar to those of Sr and about an order of magnitude slower than those of Pb. On the other hand, the diffusion of Mn2+ in apatite is about two orders of magnitude faster than the diffusion of (trivalent) REE when coupled substitutions according to Eqs. 4 and 5 are involved ,. The diffusion coefficients of Nd and Sm in various minerals and related oxides are plotted in Fig. 17(b). The diffusion of REE in apatite is relatively fast; when only simple REE exchange is involved, it is among the fastest in rock-forming minerals for which the data exist. Even when the chemical diffusion involving coupled exchange is considered, REE transport in apatite is considerably faster than the REE diffusion in other accessory minerals .
- The minerals with essential rare-earth elements (REE or rare-earth metals, REM) or chemically related elements Y or Sc are termed as rare-earth minerals. They must be named with suffix (Levinson modifier ,), indicating the dominant rare-earth element (some examples can be found in Chapter 1 (Table 3)). Please see also note 2.
- Rare-earth elements or metals (REE or REM) are Sc, Y and lanthanoids . Light rare-earth elements (LREE) are Sc, La, Ce, Pr, Nd, Pm, Sm, Eu and Gd (7 elements from La to Eu are known as the cerium group or cerium-group lanthanides). Heavy rare-earth elements (HREE) are Y, Tb, Dy, Ho, Er, Tm, Yb and Lu. The definition of LREE and HREE is based on the electron configuration. LREEs possess unpaired 4f electron from 0 to 7 (half-filled 4f electron shell). HREEs have paired electron (the clockwise and counterclockwise spinning election) ,. The element with half-filled f-electron shell (Eu) shows enhanced stability of its particular electron configuration . In some cases, REEs are divided into three groups including (1) LREE (La – Pm), (2) MREE (middle rare-earth element, Sm – Dy) and (3) HREE (Ho – Lu) .
- The actinides occupy the second row of the f-block in the periodic table. The actinide group or actinoids (An) include 14 elements with atomic numbers from 90 (Th) to 103 (Lr) . The elements with atomic numbers greater than 92 (U) are termed as transuranes. The elements with atomic numbers greater than 100 are named as the super-heavy elements (SHE). There is also the concept of the periodic table developed by G.T. Seaborg predicting a new inner transition series of 32 elements (from 122 to 153 element), called the superactinite series ,. Only actinium, thorium, uranium and (in trace quantities) protactinium and plutonium are primordial, while the elements from neptunium onwards are present on Earth solely through artificial generation .
- Vegard’s law, first pronounced in 1921, states that the lattice parameter of a solid solution of two phases with similar structures is a linear function of lattice parameters of the two end-members .
- The elements from the chalcogenide group (or oxygen group family) belonging to Group 16 (VI A) of the periodic table: O, S, Se, Te and Po. Elements sulfur, selenium and tellurium are also termed as the elements from the sulfur subgroup .
- See also Section 5.3.
- For the purposes of nomenclature, a complete solid-solution series without structural order of ions defining the end-members is arbitrarily divided at 50 mol.% (“50% rule”) to two portions with different names. Analogously, the 50% rule applied to members of ternary solid-solution series implies that the mineral names should be given only to the three end-members. Each name should be applied to the compositional range from the end-member to the nearest right bisector of the sides of the composition triangle. For example, the apatite series Ca5(PO4)3(F,OH,Cl) is represented by three compositional fields of fluorapatite (A), hydroxylapatite (B) and chlorapatite (C) .
- Statistical model where a random sequence k, the probability distribution p(k) of which satisfies the equation :
- The name is also often used for synthetic phase of the same composition (dicalcium silicate, Ca2SiO4, C2S), which is the main component of belite in Portland cement ,,.
- Since sanidinite facies are formed under conditions of intensive contact metamorphosis at high temperatures and low pressure, volatiles such as carbon dioxide and water are removed from the rock. The components characteristic for the sanidinite facies are sanidine (feldspar, KAlSi3O8 ), corundum (oxide, Al2O3 ), cordierite (cyclosilicate, Mg2Al4Si5O18 ), sillimanite (nesosilicate, Al2OSiO4 ) and glass formed as the product of partial fusion.
- The abbreviation for per formula unit (pfu), see also Footnote 36 in Chapter 1.
- The utilization of apatites as catalysts is described in Section 10.7.
- The minerals of sedimentary rocks are subdivided into two main groups : authigenic (formed on their present position) and allogenic (transported to its current position from elsewhere). Both groups will be further described in Chapter 7.
- The Margules equation expresses the Gibbs free energy (GE) of binary liquid mixture (x1 + x2 = 1) in the symmetric form :