Open access peer-reviewed chapter

Identification, Characterization and Properties of Apatites

Written By

Petr Ptáček

Reviewed: January 8th, 2016 Published: April 13th, 2016

DOI: 10.5772/62211

Chapter metrics overview

2,492 Chapter Downloads

View Full Metrics

Abstract

The parameters of unit cell, structure, refractive index, solubility data, PO4/CO3 ratio, surface area, etc., are important parameters for characterization of phosphate rocks. Third chapter of this book introduces methods for identification, characterization and properties of apatites in four main sections. The first part describes techniques used for identification and investigation of properties of phosphate minerals, including X-ray diffraction analysis, powder neutron diffraction, X-ray fluorescence as well as spectroscopic and microscopic methods. Some of these techniques are then demonstrated on the fluorapatite specimen in the second part. The third part of this chapter deals with thermodynamic properties of apatite-type compounds and introduces some of thermodynamic predictive methods. The fourth part is dedicated to dissolution of apatite, where the reaction between solids and liquids according to different dissolution models is described. Chapter ends with methods for the evaluation of reactivity of phosphate rocks.

Keywords

  • Apatite
  • Carbonate to Phosphate Ratio
  • Dissolution of Apatite
  • Reactivity of Apatite
  • Citrate Solubility

3.1 Techniques used for identification of phosphate minerals

There are several parameters which are used in characterization of phosphate rocks, including the following [1],[2]:

  1. Unit-cell parameters

  2. Refractive index

  3. Solubility data based on chemical extraction methods, e.g. neutral ammonium citrate solubility (NAC, Section 3.4.3)

  4. PO4/CO3 ratio as a measure of carbonate substitution in phosphate minerals

  5. Surface area and pore size distribution indicating the potential reactivity

The most important techniques used for the identification and characterization of phosphate minerals include methods for identification of phase composition, chemical composition, structure, surface properties, etc. A few often applied methods are introduced in this chapter [1],[2].

3.1.1. X-ray diffraction analysis

The diffraction of a beam of X-rays by a crystalline material is the process of beam scattering

Scattering is the process where the beam of radiation or particles is deviated from its initial trajectory by the inhomogeneity in the medium which it transverses [7].

by electrons associated with atoms in the crystal and of the interference of these scattered X-rays because of the periodic arrangement of atoms in the crystal and its symmetry.

Other kinds of radiation commonly used for diffraction analysis are neutrons (Section 3.1.2) and electrons (Chapter 3.1.10).

X-ray diffraction analysis (XRD) is used for the determination of mineralogical composition and quantitative X-ray diffraction analysis (Rietveld method) for the refinements of structure of apatite from measured data using specialized software [3],[4],[5],[6].

Following the discovery of X-rays by Röentgen

Wilhelm Conrad Röentgen (1845–1923) was the rector of the University of Würzburg [8]. The first X-ray photography was published in 1896 [9].

[8],[10],[11] in 1895 and the proof that X-rays have the wave properties and diffract from a periodic atomic array by von Laue [12] and his students in 1912–1913, the analytical application of X-ray diffraction has developed slowly over the next 20 years. Most of the earliness efforts were aimed at the solution of crystal structures of common phases. Debye and Sheerer (1916) and Hull (1917) suggested that powder diffraction patterns could be used for the identification of quantification of crystalline compounds. However, because most of the early developments were directed toward solving single-crystal structures, it was really the middle 1930s when the powder diffraction method began to attract the follower with the publication of the procedure of Hanawalt and Rinn and the database of patterns by Hanawalt, Rinn and Frevel (1938). With the conversion of data sets into the first set of the Powder Diffraction File in 1941, the phase identification applications expanded, and the modern counter diffractometer was developed by Parrish, Hamacher and Lowitzsch [5].

The phase identification was one of the first applications to grow to useful level. Other major applications of diffraction analysis include following phase changes under nonambient conditions and atmospheres. The first diffraction experiments were actually done on single crystals. The method is primarily directed toward determining the crystal periodicity and symmetry and solving the arrangement of atoms in the material because this information is difficult to obtain from powder experiment. Powder diffraction is one of the most important material characterization techniques in the material research and industry [5].

Single-crystal diffraction studies are not limited to crystal structure analysis. The diffraction topography is a large field that has provided much information on the perfection of crystals used in industry as integral parts of devices. The examples include crystals used in sensing and control devices, substrates for electronic components, tools and dies, turbine blades and many other applications [5].

3.1.2. Powder neutron diffraction

Neutron powder diffraction (PND) or elastic neutron scattering enables to determine nuclear and magnetic structure of solids. Most of the information on the nature of ordered magnetic phases or magnetic structures comes from neutron diffraction experiments. Neutrons have no electric charge and interact with the nuclei rather than with the charge distribution of atoms in matter. They have the wavelength in the range of interatomic distances. They have magnetic moment and interact with the magnetic moment of atoms in matter. The mass of neutrons is similar to that of atomic nuclei; hence, they have energy and momentum similar to those of atoms in solid and fluid materials [13],[14]. The first neutron diffraction experiments were performed in 1945 by EO Wollan in the graphite reactor at Oak Ridge National Laboratory, USA [15].

Neutron scattering (NS) results from the interactions with atomic nuclei, i.e. overscattering lengths (distances) of the order of 10−15 m (1 fm).

The interaction with electrons during X-ray analysis takes place over the distances of 10−10 m (1 Å) [14].

Although scattering amplitude decreases greatly with the scattering vector (it is inverse to the scattering length), there are insignificant variations of scattering amplitude in the same range of scattering vector for neutrons. Consequently, powder diffraction with neutrons can resolve very fine structural and textural details of complex atomic structures. Moreover, the weak interaction of neutrons with matter results in very low attenuation offering a unique advantage for nondestructive, in situ work and bulk analysis (for polycrystalline materials, no crushing is required to obtain the patterns) [14].

Neutron diffraction was used to determine atomic arrangement in material. Inelastic neutron scattering measures the vibrations of atoms and small-angle neutron scattering

Small angle scattering (SAS) of X-rays is abbreviated as SAXS [7].

(SANS) is used to study larger structures such as polymers and colloids. The technique of surface reflection (reflectometry) was used to study layered materials. The technique of SANS provides the information about the size, shape and domain orientations; conformational changes and/or flexibility; and molecular associations is solution. For the structural studies, the elastic scattering effects, where there is no energy exchange between the radiation and atoms, are exploited [7],[14],[16].

Neutron powder diffraction is a method often used for the structure refinement of apatite or apatite type compounds from measured data using specialized software [17],[18],[19],[20], [21],[22],[23],[24] and the effect of substitutions in the apatite structure [25]. Neutron diffraction data enable to explain the oxygen over-stoichiometry in the structure of La9.67(SiO4)6O2.5 apatite [17]. Since neutrons make possible the accurate determination of the thermal factors and provide the visualization of the diffusion paths in ionic conductors, powder neutron diffraction is also used for the characterization of solid oxide fuel cell materials [26]. This method is also used to investigate apatite in hard tissues where it provides the evidence about the deficiency of hydroxyl ion in bone apatites [27] and reconfirms that the inorganic portion is basically a hydroxylapatite-like material [27],[28].

3.1.3. X-ray fluorescence analysis and total-reflection X-ray fluorescence analysis

X-ray fluorescent analysis (XRF) is a method for the determination of sample composition [29],[30]. The origin of characteristic X-ray spectra can be described as follows. When sufficient energy is introduced into the atom, the electrons may be knocked out of one the inner shells. The atom is then in an excited (ionized) state and returns to the ground state within 10−8 s. The place of the missing electron is filled by an electron from a neighboring other shell, the place of which, in turn, is filled by an electron from more outer shell. The atom then returns to the ground state in steps. In every step, i.e. in every electron jump, the electron from a higher energy level goes into a lower energy level emitting excess energy in the form of an X-ray quantum. The energy of emitted radiation is characteristic for the atomic number of emitting element as well as for particular electron transition taking place within the electron shell of the atom. By measuring the energy or the wavelength of emitted radiation, the particle element can be identified unambiguously [31].

The energy that is necessary for the atom to get to excited state can be introduced either by the collision with a high-energy electron (sample is bombarded by electrons which are accelerated by high-voltage) or by the absorption of an energy-rich photon, i.e. the X-ray quantum (sample is irradiated by X-ray or gamma rays). In modern X-ray fluorescence analysis, the sample is irradiated by polychromatic radiation from an X-ray tube. In analogy to the optical case, this technique is referred to as fluorescence, which is responsible for the name X-ray fluorescence analysis as the technique of spectrochemical analysis with X-rays [31].

There are two types of instruments (Fig. 1) used for X-ray fluorescence spectrometry [32],[33]:

  1. Wavelength-dispersive XRF (WDXRF) or total reflection XRF (TRXRF): the method is also often abbreviated as XRF. X-rays impinge on the sample (Fig. 1(a)) and generate fluorescent X-rays. These are then diffracted on a crystal. A goniometer selects the geometry between the crystal and detector that controls the detection of X-ray from the element of interest. Different crystals have different sensitivities. Many of commercial WDXRF instruments have two detectors and up to six crystals to optimize the conditions for each element.

  2. Energy-dispersive XRF (EDXRF): the method is also abbreviated as EDX. The EDXRF instruments use much less energetic X-ray tube. Emitted X-ray radiation from the sample impinges directly on a detector, typically Si(Li), which generates pulses on an incident beam. These pulses are sorted and counted by a multichannel analyzer (Fig. 1(b)). Simultaneous determination for all elements, the atomic number of which is greater than Mg is possible.

The resolution and sensitivity of EDXRF is typically an order of magnitude worse than that for WDXRF.

Fig. 1.

Schematic representation of X-ray fluorescence analyzer: (a) wavelength-dispersive (XRF) and (b) energy dispersive (XRF).

Synchrotron radiation X-ray fluorescence (SRXRF) microprobe, a promising technique, is a nondestructive and qualitative to semiquantitative analysis of minerals and single fluid inclusions [34]. Synchrotron radiation (SR) is a powerful advanced light source (synchrotron radiation source, SRS) compared to conventional X-ray tube radiation and has many unique properties, such as high intensity, natural collimations, well-defined polarization, wide spectral range and energy tenability [35]. SRXRF is a widely applied technique for microscopic analysis of chemical elements. High-resolution requirements can be achieved using microbeam synchrotron radiation X-ray fluorescence (μ-SRXRF). Synchrotron radiation X-ray fluorescence can also provide the information about the oxidation state and coordination environment of metals using techniques known as X-ray absorption of near-edge structure (XANEX) or by micro-XANEX spectroscopy [37]. The unique tool for studying, the local structure around selected elements is X-ray absorption fine structure (XAFS) [38].

X-ray fluorescence is usually used to investigate the composition of apatite rocks and minerals for the purpose of their characterization [36],[39],[40], estimation of naturally occurring radionuclides in fertilizers [41] and analysis of phosphate ore at various stage of processing [42], e.g. flotation [43],[44].

3.1.4. Inductively coupled plasma spectrometry

Prior to inductively coupled plasma

Plasma is defined as an electrically neutral gas which consists of positive ions and free electrons. Plasma have sufficiently high energy to atomize, ionize, and excite virtually all elements in the periodic table, which are intentionally introduced into it for the purpose of elemental chemical analysis [45].

(ICP), the flame, direct current-arc and controlled-waveform spark were used for the atomization (i.e. decomposition of sample to individual atoms) and excitation of sample in elemental analysis. ICP denotes the technique that uses atmospheric pressure argon inductively coupled plasma

Although there are many types of plasma, such as direct current, microwave induced, etc., the ICP is considered the most useful technique for analytical spectroscopy [45].

(ICP) for the atomization and excitation of sample. This plasma is a highly energetic media consisting of inert ionized gas with equivalent temperatures from 7000 to 10,000 K. Inductively coupled plasma are formed by coupling energy produced by RF generator (typically 700–1500 W) to the plasma support gas with an electromagnetic fields [45]. The cross-section of typical ICP torch is shown in Fig. 2.

Fig. 2.

Plasma with torch assembly and load coil [45].

The treatment of sample before the quantification includes vaporization, atomization, excitation and ionization.

Each element has characteristics first and second ionization potential, which depends on specific electronic structure of given element. Higher ionization potential means that more externally applied energy is required for ionization (thermal radiation, collision with other ion or electron, or exposure to high-energy photons) [45].

The introduction of analyzed sample into inductively coupled plasma was applied in analytical techniques including [29],[45],[46],[47],[48]:
  1. Inductively coupled plasma atomic emission spectroscopy (ICP-AES): the method connects atomic emission spectrometry with ICP. AES is based on spontaneous emission of free atoms or ions when the excitation is performed by thermal or electric energy. The method can identify and determine the concentration of up to 40 elements simultaneously with the detection limit of parts per billion (ppb). Serious limitations of this technique are the spectral interferences. Despite rapid growth of ICP-MS, ICP-AES still plays a dominant role in elemental analysis of geological, environmental, biological and other materials.

  2. Inductively coupled plasma atomic fluorescence spectroscopy (ICP-AFS): ASF is an analytical technique for the determination of elements in small quantities. It is based on the emission of free atoms when the excitation is performed by radiation energy.

  3. Inductively coupled plasma-mass spectroscopy (ICP-MS): is analytical technique for the determination of elemental composition of virtually any material. A sample, usually in the form of an aqueous solution, is converted to an aerosol by a neutralization process and transported to the plasma by an argon gas stream. In the plasma, the elements of analyte are atomized, followed immediately by ionization. The composition of ion population in the plasma is proportional to the concentration of analyte species in original sample solution. Ions produced by ICP are representatively sampled and extracted from the plasma; next they are separated and measured by a quadrupole or time-of-flight mass spectrometer

The method known as laser ablation—inductively coupled plasma—mass spectrometry (LA-ICP-MS) is a coupling technique of laser ablation with ICP-MS technique [49],[50],[51],[52], [53]. Multicollector inductively coupled plasma mass spectrometry (MC-ICP-MS) was applied as the benchmark method for isotopic analysis [54] and for the determination of heavy rare earth elements in apatites [55],[56],[57].

3.1.5. Thermal ionization mass spectroscopy

Thermal ionization mass spectroscopy (TIMS) is highly specialized technique of mass spectroscopy used for very precise determination of isotope ratios and, as such, is widely used for the determination of stable isotope ratios in isotope geology and for the analysis of nuclear materials [58]. In this method, the solid sample is thermally ionized in solid-source mass spectrometer. Ions are accelerated into the mass analyzer and then transported to the detector [59]. Multicollector thermal ionization mass spectrometry (MC-TIMS) uses laminated magnetic sector filet for high speed peak jumping and low hysteresis for the mass of ion beams [60],[61].

The basic variants of the method of thermal ionization mass spectroscopy are as follows:

  1. Isotope dilution–thermal ionization mass spectrometry method (ID-TIMS) [62],[63] is used for accurate determination of element concentration and is generally considered to be the definitive method to other techniques [58]. The method was applied for direct measurement of uranium, thorium, lead, etc., concentrations in the determination of single grain fission-track ages (Section 7.3.3).

  2. Chemical abrasion–thermal ionization mass spectrometry method (CA-TIMS)

The ability for precise determination of isotope ratio in apatite predetermines this method for geochronological investigations, i.e. studying of chronologic records in accessory minerals of igneous rocks [61],[63],[64],[65].

3.1.6. Secondary ion mass spectrometry

Secondary ion mass spectroscopy (SIMS) is a method

The method was originally developed in the 1950s and 1960s by Herzog et al and Honig et al to analyze metals and oxides [69], but the basis of SIMS can be traced back to the beginning of the twentieth century with the first experimental evidence of secondary ions given by JJ Thomson in 1910 [68].

for the characterization of solid surface elemental composition and isotope distribution. The technique can be applied to all elements and allows quantitative analysis of solid surfaces, including monolayers. Energetic ion bombardment of a solid surface (primary ions, e.g. Ar+, Cs+, O2+, … ) causes that atoms of the sample are shifted from their original states as positive and negative ions which are termed as secondary ions. These ions are then analyzed by mass spectrometer (e.g. quadrupole mass spectrometer) to determine the composition of the surface of sample [66],[67]. The applications of SIMS can be broadly subdivided into static and dynamic SIMS. Static SIMS (SSIMS) is used to investigate the composition of the outermost monolayer on any solid. Dynamic SIMS (DSIMS) examines the concentration profile as the function of depth [68],[69],[70].

Fig. 3.

Main components of SIMS [68] and schematic representation of ion microprobe imaging time-of-flight secondary ion mass spectrometry with reflectron-based mass analyser [69].

The mass filter of mass spectrometer of SIMS instruments defines the type of instruments [68], [70]:

  1. Magnetic sector SIMS (M-SIMS) instruments (Fig. 3(b), similar to those used in original mass spectrometer)

  2. Quadrupole SIMS (Q-SIMS) instruments (Fig. 4(a), first appeared in the 1970)

  3. Time-of-flight SIMS (ToF-SIMS) instruments (Fig. 4(b), first appeared in the 1980)

Fig. 4.

Schematic representation of magnetic sector and quadrupole SIMS instruments [70].

Secondary ion emission begins when the primary ions energy exceeds some threshold level. This energy ranges from 30 to 80 eV and is much higher than the energies for sputtering of neutral particles. The most important features of the method are very low sensitivity limit for majority of elements, the possibility to determine the concentration profile, the possibility of isotope analysis and the identification of all the elements and isotopes starting from hydrogen [66].

The SIMS method is often used for the measurement of deuterium/hydrogen (D/H) ratio and hydroxyl (OH) content in anhydrous minerals and melt inclusions in Martian meteorites [71] and lunar materials (apatite grains, glass beds, melt inclusions and agglutinates in soils) [72]. The method is widely used to investigate biological apatites and collagen apatite composites [73] and the analysis of micrometer-sized samples like, e.g. interplanetary dust, presolar grains and small inclusions in meteorites, has become more and more important in cosmochemistry [74]. The method is also utilized for U-Th-Pb dating of apatites as common accessory in igneous rocks (Chapter 8), based on the radioactive decay of U and Th [75],[76],[77].

3.1.7. Laser secondary neutral mass spectrometry

Laser secondary neutral mass spectroscopy (laser-SNMS) can be further divided to nonresonant laser-SNMS (NR-laser-SNMS) and resonant laser-SNMS (R-laser-SNMS). In NR-laser-SNMS (Fig. 5(b)), an intense laser beam is used to nonselective ionization of all atoms and molecules within the volume intersected by the laser beam. Sufficient laser power density, which is necessary to saturate the ionization process, is typically achieved in a small volume. It limits the sensitivity of the method and leads to the problems with quantification due to the differences between effective ionization volumes of different elements. Laser-SNMS method has significantly improved ionization efficiency over SIMS (a) [78],[79].

Fig. 5.

Comparison of ToF-SIMS and laser-SNMS: (a) direct analysis of secondary ions (ToF-SIMS), (b) nonresonant laser-post-ionization of secondary neutrals (NR-laser-SNMS) and (c) resonant post-ionization of secondary neutrals (R-laser-SNMS) [78].

Resonant laser-SNMS uses a resonance laser ionization process, which selectively and efficiently ionizes atoms and molecules over a relatively large volume (Fig. 5(c)). The method has unit ionization efficiency for over 80% of elements in periodic table, i.e. the overall efficiency is greater than that for NR-laser-SNMS, the quantification is also simpler and extremely high selectivity prevents almost all isobaric and molecular interferences [78],[79].

The method of NR-laser-SNMS was used by Dambach et al [80] to investigate different states of biomineralization in vitro. The results indicate that in the vicinity of single osteoblasts, extracellular enrichment of potassium typically occurs during initial stages of mineralization. Potassium may interact with matrix macromolecules and prevent an uncontrolled apatite deposition. However, apatite biomineral formation is correlated with a potassium release. In conclusion, potassium seems to be involved in the process of extracellular matrix biomineralization.

3.1.8. Electron paramagnetic resonance

The concept of electron paramagnetic resonance (EPR) spectroscopy is very similar to more familiar nuclear magnetic resonance (NMR). Both methods deal with the interaction between electromagnetic radiation and magnetic moments. In the case of EPR, the magnetic moments arise from electron rather than nuclei. The term EPR was introduced as a designation taking into account contributions from electron orbital as well as spin angular momentum. The term electron spin resonance (ESR) was also widely used because in most cases the absorption is linked primarily to the electron-spin angular momentum [82],[83]. EPR spectrum is a diagram in which the absorption of microwave frequency radiation is plotted against the magnetic field intensity [83].

The technique of electron paramagnetic resonance spectroscopy may be regarded as the consequence of the SternGerlach experiment. They showed (in 1920) that an electron magnetic moment in an atom can take on only discrete orientation in a magnetic field, despite the sphericity of the atom. Subsequently, Uhlenbeck and Goudsmit liked the electron magnetic moment with the concept of electron spin angular momentum. In hydrogen atom, there is additional angular momentum arising from the proton nucleus. Breit and Rabi described the resultant energy levels of hydrogen atom in a magnetic field. Rabi et al [81] studied the transition between levels induced by an oscillating magnetic field, and this experiment was the first observation of magnetic resonance. The first observation of electron paramagnetic resonance peak was made in 1945 by Zavoisky, who detected the radiofrequency absorption line from CuCl2·2H2O sample using the radiofrequency (RF) source at 133 MHz [82].

The major components of EPR spectrometer are shown in Fig. 6. The microwave bridge supplies the microwaves at controlled frequency and power, which are transmitted to the sample cavity via the waveguide. The sample cavity is placed perpendicular to applied magnetic field, which can be varied in controlled way. In addition to this main magnetic field, a controlled but smaller oscillating magnetic field is superimposed on the cavity via the Zeeman modulation frequency. The ideal way to perform the experiment would be to apply a fixed magnetic field and vary the microwave frequency. However, microwave generators are only tunable over very limited ranges. Thus, the microwave frequency is fixed and applied magnetic field is varied. The magnetic field is applied until it reaches the value at which the sample will absorb some of the microwave energy, i.e. and EPR transition occurs [84],[85].

Fig. 6.

The block diagram for typical continuous wave EPR spectrometer [84].

Electron paramagnetic resonance (EPR) spectrum of X-irradiated sodium and carbonate containing synthetic apatites has been studied by Moens et al [86]. Observed spectra were decomposed in terms of five theoretical curves representing O radical, two CO3– radicals (surface and bulk) and two CO2− radicals (surface and bulk). These species were also described in A-type and B-type carbonate-apatites [87],[88], tooth enamel [89],[90],[91],[92], [93] and bone [94],[95], apatites, renal stones [96], etc.

3.1.9. Nuclear magnetic resonance

Solid-state nuclear magnetic resonance (NMR) is a technique for accurate measurement of nuclear magnetic moments where the resonance frequency depends on its chemical environment [97],[98],[99]. The method can provide useful information on the number of molecules in the asymmetric unit and on the site symmetry of the molecule in the lattice to assist in the refinement of powder X-ray diffraction (Section 3.1.1) data. The method can distinguish between different polymorphs. Alternatively, solid-state NMR can be used for direct and accurate measurement of internuclear distances. For amorphous and disordered solids, such as inorganic glasses and organic polymers, solid-state NMR provides structural information that cannot be obtained by any other technique [100],[101]. NMR is also the diagnostic method used in veterinary science and medicine particularly in clinical research of human brain by magnetic resonance imaging (MRI) [102].

The solution-state NMR method was developed for the investigation of structure of soluble proteins [103]. Solution and solid-state NMR are both excellent methods for the determination of chemical composition [100].

The structural information of apatites is usually investigated from 1H, 19F and 31P NMR spectra of apatites [104]. The 31P solid-state NMR spectroscopy is a useful tool to investigate structural information about apatites on bone organic and inorganic mineral components, as well as to investigate the crystallinity and compositional changes in carbonated apatites [105]. Intact bone is a demanding tissue for structural studies. Serious experimental problems arise from the morphological diversity of bone and from the co-existence, interrelationship and great complexity of its organic and inorganic components. Furthermore, one has to perform noninvasive analysis because bone samples are very sensitive to physical effects and chemical treatment. Solid-state 31P NMR gives us a unique opportunity to look specifically at the minerals of whole bone without any chemical pretreatment, thus avoiding the intervention into the bone structure [106].

Fig. 7.

Nuclear separation along parallel chains (the crystallographic c-axis) in various apatites (a): FFF group (I), FFH group (II) and HFH group (III). Correlation between observed 19F line width and fluorine content of fluorinated hydroxylapatite (b) [104].

The 19F NMR spectrum of fluorinated calcium hydroxylapatite (Ca10(PO4)6F2x(OH)2–2x, where x is the fraction of OH replaced by F) indicates the correlation between 19F chemical shift tensor parameters and the content of fluorine in apatite. The presence of OH groups induces perturbations of fluorine environments, involving the displacements of both fluorine and hydroxyl groups from their normal positions. This leads to a distortion of the electronic environment with regard to the investigated fluorine nucleus and gives reasons for observed change in the 19F chemical shift tensor of fluoridated hydroxylapatite with different fluorine content. Furthermore, the presence of OH-group destroys the fluoride long-range structure and that results in an isotropic chemical shift distribution. This leads to observed increase in the 19F line width in the case of low fluorine content [104],[107].

3.1.10. Scanning electron microscopy, structure and elemental analysis

The scanning microscope (SEM) permits the observation and characterization of heterogeneous organic and inorganic materials on a nanometer (nm) to micrometer (μm) scale. In SEM, the area to be examined or the volume to be analyzed is irradiated with finely focused electron beam, which may be swept in a raster across the surface of the specimen to form an imager or may be static to obtain the analysis at the position. The type of signals produced from the interaction of the electron beam (primary electron, PE) with the sample (Fig. 8(a)) includes secondary electrons (SE, with energy ≤50 eV), backscattered electrons (BSE, E > 50 eV), Auger electrons (AE), X-ray characteristics (X) and other photons of various energies such as continuum X-rays and heat. Low-loss electrons (LLE) show the energy losses of a few hundreds of eV. These signals are obtained from specific emission volumes within the sample and can be used to examine many characteristics of the sample such as surface topography, crystallography, composition, etc. [108],[109],[110].

Fig. 8.

Electron–specimen interaction (a) and schematic energy spectrum (a) [109].

Secondary and Auger electrons are highly susceptible to elastic and inelastic scattering and can leave the specimen only from a very thin surface layer of the thickness of a few nanometers. The most probable energy of BSE falls into the broad part of the spectrum in Fig. 8(b), but they also show more or less pronounced elastic peak followed by plasnom losses, which depend on the primary energy, the take-off angle and the tilt of the specimen. Continuously slowing-down approximation assumes that the mean electron energy decreases smoothly with decreasing path length of the electron trajectories inside the specimen. The maximum information depth of BSE is of the order of half the electron range. Characteristic X-rays will only be excited in the volume in which the electron energy exceeds the ionization energy of the inner shell involved. Inelastic scattering in semiconductors results in the generation of electron-hole pairs. The recombination can take place without radiation but may result in the emission of light quanta (cathodoluminescence, CL) [111].

The method known as electron backscattering diffraction (EBDS) enables to determine the crystal structure and grain orientation of crystals on the surface of specimen. To collect maximum intensity in the diffraction pattern, the surface of specimen is stipple tilted at an angle of typically 70° from the horizontal (Fig. 22(a)). The intensity of backscatter Kikuchi patterns (please see the pattern of fluorapatite in Fig. 23) is rather low, as is the contrast of the signal, so extremely sensitive cameras and contrast enhancement facilities are required. This pattern allows to identify the phases and shows the misorientation across the grain boundaries [108].

Scanning electron microscope can be also used to determine compositional information using characteristic X-ray. The development of instruments for obtaining localized chemical analysis of solid samples, i.e. electron probe microanalyzer (EMPA), occurred at the same time as the development of SEM.

Scanning electron microscopy (SEM) is used for grain interactions and spot analysis [98],[112], electron microprobe microanalysis (EPMA) for the distribution of elements in the matrix, investigation of the effects of impurities on the properties of apatites and investigation of reaction interface [113],[114],[115],[116].

3.1.11. Fourier transform infrared and Raman spectroscopy

Infrared (IR) spectroscopy is one of the most important analytical techniques that can be used for the investigation of any sample in any state. Liquids, solutions, pastes, powders, films, fibers, gases and surfaces can be examined with judicious choice of sampling technique. Infrared spectrometers have been commercially available since the 1940s [117].

Fourier transform infrared (FT-IR or FTIR) spectroscopy is divided into three regions according to the increasing wavelength [118]:

  1. Near-IR (NIR) spectroscopy, abbreviated as FT-NIR

  2. Mid-IR (MIR) spectroscopy, abbreviated as FT-MIR

  3. Far-IR (FAR) spectroscopy, abbreviated as FT-FAR

The spectral ranges of near-, mid- and far-infrared spectroscopy are shown in Fig. 9(a).

Fig. 9.

Schematic illustration of relationships between the ranges of (a) vibrational spectroscopy and electromagnetic spectrum [118] and (b) spectroscopic transitions underlying several types of vibrational spectroscopy. v0 indicates the laser frequency, while v is the vibrational quantum number. The virtual state is a short-lived distortion of the electron distribution by the electric field of the incident light [119].

The background for Raman spectroscopy was given by the discovery of Raman scattering by Krishna and Raman in 1928. Until approximately 1986 when Fourier transform (FT)–Raman was introduced, physical and structural investigations dominated in literature over relatively few reports of Raman spectroscopy applied in chemical analysis [119],[120].

When monochromatic light with the energy 0 encounters the matter (gas, solid or liquid), there is a small probability that it will be scattered at the same frequency (Fig. 9(b)). If the object, e.g. molecule is much smaller than the wavelength of the light, the scattering is Rayleigh scattering. The “virtual state” is not necessarily a true quantum state of the molecule but can be considered a very short-lived distortion of the electron cloud caused by oscillating electric field of the light. Since blue light is more efficiently scattered than red one, Rayleigh scattering is responsible for the blue color of sky. The electron cloud of the molecule is also perturbed by molecular vibrations, and it is possible for the optical and vibration oscillations to interact, leading to Raman scattering. Raman scattering is shown in (Fig. 9(b)) in which the scattered photon is lower in energy by an amount equal to the vibration transition. Raman spectrum consist of scattered intensity plotted versus energy and each peak corresponds to given Raman shift from the incident light energy 0 [119].

Just like Rayleigh scattering, Raman scattering depends on the polarizability of scattering molecules. IR band, on the other hand, arises from a change in the dipole moment. In many cases, the transitions that are allowed in Raman are forbidden in IR, so these techniques are often complementary (please compare Fig. 18(a) and (b)). In polarizable molecules, incident light can excite the vibrational modes, leading to scattered light diminished in energy by the amount of vibrational transition energies (same as in fluorescence). Scattered light under these conditions reveals the satellite lines below the Rayleigh scattering peak at the incident frequency–Stokes lines (Stokes part of spectrum). If there is enough energy, it is also possible to see anti-Stokes lines. Since anti-Stokes lines are usually weaker than Stokes lines, only the Stokes part of spectrum is usually measured [121].

The method combining Raman spectrometer with microscopic tools, typically an optical microscope, is known as micro-Raman spectroscopy (μRS) or also Raman microscopy. The μRS is nondestructive and noncontact method for the characterization of organic and inorganic materials [122].

Infrared [97],[98],[123],[112],[124],[125],[126], Raman [97],[98] and micro-Raman spectroscopies [125] were often used to identify and investigate the structure and extent of substitution and to optimize the synthesis conditions of minerals from the supergroup of apatite. Since carbonate ions exhibit clear vibrational signature in infrared spectrum, infrared spectroscopy is widely used to investigate the structure and to evaluate the carbonate/phosphate ratio (rc/p) and the amount of carbonate ions in carbonate-apatites [127]:

C O 3 2 wt . % = 28.62 r c / p + 0.0843 E1

The example of infrared and Raman spectrum of fluorapatite is described in Section 3.2.3.

Infrared spectra of phosphate minerals in the pyromorphite series are described by Adler [128]. In the pyromorphite series, the equilibrium internuclear X-O distance in XO43− ion (PO43−, AsO43− and VO43−) is primarily a function of the ionic radium of X atom. Since Pb, in this case, is always the dominant externally coordinated cation, for various members, there is no significant change in the interaction between the molecular vibration and the external environment. Bradger’s equation [128],[129],

k 0 = 1.86 10 5 / R d i j 3 E2

although specifically applicable to internuclear distances in diatomic molecules, reflects generally the inverse relationship between the force constant k0 and the internuclear distance R. Symbol dij denotes the constant the values of which depend on the nature of bonded atoms. The molecular vibration frequency v is dependent on the restoring forces, measured in terms of k0, between participating atoms as well as on the masses of these atoms. The relationship may be expressed approximately by the equation:

v = 1 2 π c k u E3

where the vibration frequency v is a function of the force constant k and the reduced mass u of vibrating atoms, all other terms being invariant.

Fig. 10.

Theoretical effect of change in mass and ionic radius on infrared vibration frequency of tetrahedral XO43− ions, where X = P5+, As5+ of V5+ [128].

The spectral frequency differences between pyromorphite, mimetite and vanadinite are explicable and to a considerable degree predictable in terms of these parameters. On complete substitution of As or V for P the effect of reduced force constants is reinforced by increases in mass, thereby shifting ν3 and ν1 to lower frequencies. Because of opposing mass and force-constant effects and perhaps also because of dissimilarities in orbital configurations, the relative positions of absorption bands are less predictable for mimetite and vanadinite than for pyromorphite and mimetite. The theoretical frequency trends are depicted in Fig. 10 [128].

3.1.12. Thermal analysis

Thermal analysis (TA) refers to a group of techniques

The definition of terms in thermal analysis was developed by ICTAC (Confederation for Thermal Analysis and Calorimetry).

in which the property of a sample is monitored against time or temperature while the temperature of the sample, in a specified atmosphere

Gaseous environment of the sample, which may be controlled by the instrumentation or generated by the sample [130].

is programmed. These methods study the relationship between sample property and its temperature as the sample is heated or cooled in a controlled manner. The individual techniques are divided according to the measured property

Resulting dependence, i.e. any graph of any combination of property vs. time or temperature derived from a thermal analysis technique, should be termed as thermal curve, which is a simplified form of more correct term thermoanalytical curve. The first mathematical derivation of any curve with respect to temperature or time leads to the derivative thermoanalytical curve [130]. Since the name thermogram has medical usage, the thermal analysis curve should not be termed as thermogram [131].

[131],[132], [130] as is introduced in Table 1.

Property Technique and abbreviation Notes
Heat Scanning calorimetry
Temperature Thermometry Also can be described as heating or cooling curve
Temperature difference Differential thermal DMA analysis DTA A technique where the temperature difference between the sample and reference material is measured
Heat flow rate difference Differential scanning calorimetry DSC A technique where the difference between heat flow rates into the sample and reference material is measured
Mass Thermogravimetry or thermogravimetric analysis TG or TGA The abbreviation TG has been used, but should be avoided, so that it was not confused with Tg (temperature of glass transition)
Dimensional and mechanical properties Dynamic mechanical analysis DMA Moduli (storage/loss) are determined
Thermomechanical analysis TMA Deformations are measured
Thermodilatometry TD Dimensions are measured
Electrical properties Dielectric thermal analysis DEA Dielectric constant/dielectric loss are measured
Thermally stimulated current TSC Current is measured
Magnetic properties Thermomagnetometry TM Often combined with TGA
Gas flow Evolved gas analysis EGA The composition and/or the amount of gas/vapor is determined
Emanation thermal analysis ETA Trapped radioactive gas within the sample is released and measured
Pressure Thermomanometry Evolution of gas is detected by pressure change
Thermobarometry Pressure exerted by dense sample on the walls of a constant volume cell is studied
Optical properties Thermoptometry or thermos optical analysis TOA A family of techniques in which optical characteristics or property of the sample is studied
Thermoluminescence TL Light emitted by the sample is measured
Acoustic properties Thermosonimetry or thermoacoustimetry TS Emitted (sonimetry) or absorbed (acoustimetry) sound is measured.
Structure Thermodiffractometry Techniques where the compositional or chemical nature of the sample is studied.
Thermospectrometry

Table 1

Methods of thermal analysis according to measured property or physical quantity.

The measurement should be performed as follows:

  1. Combined: the application of two or more techniques to different samples at the same time. This can include thermal and nonthermal analytical techniques.

  2. Simultaneous: indicates the measurement of two or more properties of a single sample at the same time.

The sample-controlled method where the feedback used to control the heating is the rate of transformation is termed as controlled-rate thermal analysis (CRTA) [130].

Simultaneous thermogravimetry (thermogravimetric analysis) and differential thermal analysis (TG-DTA) are mostly used to investigate the course of synthesis and the characterization of prepared apatites or to investigate the process of thermal decomposition of apatites, i.e. the processes such as dehydroxylation (e.g. Section 1.5.2), defluorination (Section 3.2.4 and 8.6), decarbonation (thermal decomposition of carbonate-apatites, Section 4.6.1), etc.

3.1.13. Optical properties

The analysis of optical properties is essential for each mineral examined, and through the use of microscopy, the optical properties of individual minerals may be interpreted in great detail. Optical mineralogy investigates the interaction of light (usually is limited to visible light) with minerals and rocks. Optical mineralogy concerns mainly the use of polarizing (petrographic) microscope which has two Nicol prisms, polarizer and analyzer (polarizer below the stage and the analyzer above the objective) [133]. Human eye is the most sensitive for viewing a solid in the wavelength symmetrically spread in intensity around 550 nm (Fig. 11(a)). Optical microscopy in visible light (from 700 nm (red) to 420 nm (violet), Fig. 11(b,c)) helps study the objects of smaller sizes up to lower limit of ~ 1 μm (c). Fig. 11(d) shows the scope of spectroscopic techniques associated with electromagnetic spectrum [134].

Fig. 11.

Electromagnetic spectrum of solar radiation and black body (a), the resolution of various techniques (b) and the regions associated with spectroscopic techniques (c) [134].

Refractive index (n) is related to the angle of incidence (i) and the angle of refraction (r) according to the Snell´s law:

n = sin i / sin r E4

The refractive index increases as the wavelength of light decreases [134]. The absorption coefficient is related to the imaginary part of the refractive index. It was found that since the refractive index of a medium depends on the density of electrons in that medium, the index increases with the density of matter [135].

The luster of mineral (R) depends of the way in which the light is reflected from the surface of a mineral. The reflection is again dependent on the refractive index. Normally, the greater is the index of refraction, the brighter is the luster. The luster is classified into two broad classes [134]:

  1. Nonmetallic luster: results from the interaction of light with dielectric semiconducting and poorly conducting substances

  2. Metallic luster: occurs in minerals having metallic bonding (e.g. native metals) or high degree of covalent bonding (sulfides, sulfosalts, etc.)

Optical properties of apatites can be determined using the complex dielectric function [136], [137],[138]

ε ω = ε 1 ω + i ε 2 ω E5

in the range of linear response. By calculating the wave function matrix and using Kramers–Krönig relations, the imaginary and real part of the dielectric function ε1(ω) and ε2(ω) can be derived respectively as follows:

ε 2 ω = e 2 π m ω 2 B Z d k 3 n , l Ψ n k ¯ r ¯ i Ψ 1 k ¯ r ¯ 2 f 1 k ¯ 1 f n k ¯ δ E n k ¯ E 1 k ¯ ω E6

where f() is the Fermi distribution function, l and n mark occupied state and unoccupied state, respectively. Ψn(,) is the Bloch wave function for the nth band with the energy En() at the Brillouin zone point k. The matrix element momentum transition corresponds to the term Ψ n k ¯ r ¯ i Ψ 1 k ¯ r ¯ 2 .

For each apatite, the real part ε1 and imaginary part ε2 of dielectric function have similar features with some subtle differences. The real part ε1 has two main peaks:

  1. First peak lies near the energy of 7 eV. Obviously, the first peak is caused by the transitions from O–2p and P–3p levels.

  2. The second peak lying at ~25 eV results from the transitions from Ca–3p levels.

The ε2 curve in the energy range from 5 to 15 eV characterizes three main peaks:

  1. The first peak lies in 7.10, 7.16 and 7.43 eV for FAP, ClAP and BrAP, respectively.

  2. The second peak is located at 8.51, 8.62 and 8.76 eV for FAP, ClAP and BrAP, respectively.

  3. The third peak is situated at 10.99, 11.08 and 11.45 eV for FAP, ClAP and BrAP, respectively.

The refractive index n can be obtained by n = √ε1(0). The ε1(0) for FAP, ClAP and BrAP are 1.38, 1.41 and 1.46, respectively. Therefore, the refractive indexes for FAP, ClAP and BrAP are 1.17, 1.19 and 1.21. Moreover, analogous dielectric function curves and similar refractive index values show that the optical property of each apatite has some independence from the c-axis ion [136].

3.1.14. Measuring of surface area, porosity and pore size distribution

There is a conventional mathematical idealization that asserts that a cube of edge length a possesses a surface area of 6a2 and that a sphere of radius r exhibits 4πr2. In reality, however, mathematical perfect or ideal geometric forms are unattainable (Fig. 12) since all real surfaces exhibit flaws under microscopic examination. Real surface irregularities (voids, pores, steps, etc.) make the real surface area greater than corresponding theoretical area. When the cube, real or imaginary, of one meter edge length is subdivided into smaller cubes each one micrometer (10−6 m) in length, there will be 1018 particles formed, each exposing an area of 6 × 10−12 m2. Thus, the total area of all particles is 6 × 106 m2. This operation increases by millionfold the exposed area of fine powder compared to undivided material. Whenever the matter is divided into smaller particles, new surfaces must be produced with corresponding increase in surface area [139],[140].

Fig. 12.

Description of particle shapes [140].

The particle size distribution (PSD) was usually determined by sieve analysis, sedimentation methods (gravitational or centrifugal), microscopic techniques, light scattering, multiangle laser light scattering (MALLS), etc. [141].

The range of specific surface area, i.e. area per gram of matter, can vary widely depending on the particles’ size, shape and porosity. The influence of pores can often overwhelm the size and external shape factors. The powder consisting of spherical particles exhibits total surface (St) and volume (V = M/ρ) [139]:

S t = 4 π r 1 2 N 1 + r 2 2 N 2 + + r n 2 N n = 4 π i = 1 r i 2 N i E7
V = 4 3 π r 1 3 N 1 + r 2 3 N 2 + + r n 3 N n = 4 3 π i = 1 r i 3 N i E8
S = S t M = 3 i = 1 N i r i 2 ρ i = 1 N i r i 3 E9

where ri and Ni are the average radii and numbers of particles in the fraction i. For spheres of uniform size, Eq. 9 becomes the law:

S = 3 ρ r E10

Besides the calculation of specific surface from its geometry and PSD curve, the adsorption isotherm (predominantly Langmuir and BET

BrunauerEmmettTeller (BET) theory of multilayered physical adsorption of gas molecules on a solid surface [143].

), air permeability methods,

For example, the Blaine method, where fixed volume of air passes through the bed at steadily decreasing rate, which is controlled and measured by the movement of oil in a manometer, the time required being measured. The method is widely used for the determination of specific surface area of cements [143].

and monolayer sorption methods were used to determine the specific surface area [142],[143].

The porosity is defined as the ration of pore volume to total volume. Porous material is defined as solids containing pores Fig. 13(a), which are classified into two major types: open and closed pores. Penetrating open pores (interconnected pores) are permeable for fluid and therefore are important in applications such as filters. Pores accessible from only one end are referred to as dead-end pores. Noninterconnected (closed) pores are not accessible at all. The classification of pores according to their size is shown in Fig. 13(b). Pores can be also classified as the pores among agglomerates and pores among primary particles (Fig. 13(c)) [144],[145].

Fig. 13.

Schematic illustration of different morphology of pores (a) and classification of porous materials based on pore size (b) and schematic illustration of pores among agglomerates and primary particles (c) [144].

Pore size distribution and permeability

The flow of fluids through porous materials is of great importance in the fields of industrial chemistry, oil technology, and agriculture. In general, it may be stated that the principal interest is in the transport through reactive materials [147].

are a very desirable quantities for the characterization of structure of porous solids, which can be determined by the following [145],[146],[147]:
  1. Stationary fluid (static) method

  2. Capillary pressure (Hassler

    Experimental techniques for the measurement of relative permeability can be divided to steady- (1) and unsteady-state (displacement) methods (2). The steady-state method was developed by Hassler (1944). Semipermeable membranes are provided at each end which keeps the fluids separated, except inside the core where the fluids flow simultaneously. The pressure is measured in each phase through semipermeable barriers, and the pressure differences between the phases are maintained constant throughout the medium so as to eliminate the capillary end effect as well as to ensure a uniform saturation along the core. The saturation can be altered by applying capillary pressures across the nonwetting phase ports and wetting phase semipermeable membranes [146].

    ) method

  3. Quasi-steady-state methods, including the gas drive method and the solution-gas drive method.

The surface area and the porosity of apatites have effect on the floatability (Section 8.7) of apatites of different type and origin [148]. The surface area, the porosity and the pore size distribution are properties of great importance for the preparation of biological apatites in tissue engineering (Section 10.9) and tailoring their mechanical properties, solubility and bioactivity [149],[150],[151],[152]. Ionic surfactants, such as decyltrimethylammonium bromide (CH3(CH2)9N(CH3)3(Br), C(10)TAB), hexadecyltrimethylammonium bromide (CH3(CH2)15N (Br)(CH3)3, C(16)TAB), as well as nonionic surfactant, can be used to control the pore (pore size and volume) and surface characteristics of mesoporous apatite materials under maintained pH [153]. The porosity also affects electrical properties of oxyapatites (Chapter 5) [154],[155],[156].

Advertisement

3.2. Investigation and characterization of apatite specimen

3.2.1. Identification of the specimen

Some techniques mentioned above will be demonstrated on the specimen of apatite sample (Fig. 14(a)) from Sljudjanka, Bajkal. The translucent specimen with glassy luster is greenish blue colored and brittle as can be seen from large amount of smoothly curving conchoidal fractures on the surface (please see also Fig. 20). When scratched by a single crystal of corundum (Fig. 14(b)), the sample shows white colored scratch. Since the sample surface can be also scratched by feldspar (Fig. 14(c)) but not by fluorite, the hardness in the Mohs scale

The Mohs scale of mineral hardness is graduated as follows: talc (1), gypsum (2), calcite (3), fluorite (4), apatite (5), orthoclase (6), quartz (7), topaz (8), corundum (9), and diamond (10). Apatite should be also scratched by steel knife (up to 5.5) and glass (up to 6).

is equal to 5, i.e. corresponds to apatite.

Fig. 14.

Specimen of apatite from Sljudjanka, Bajkal (a) shows white color of scratch (b). The sample can be also scratched by minerals with hardness ≥6 (feldspar and higher) on the Mohs scale (c).

The microhardness of apatite sample was then determined by Vickers microhardness test.

Hardness tester LECO AMH 43. The method is also known as the Vickers pyramid number (HV) or the diamond pyramid hardness (DPH).

Fig. 15 shows the replica of diamond pyramid base

Diamond pyramid with apical angle of 136°.

on the surface after the indentation of sample. The average hardness of sample was determined to be 552 (±30) HV 0.05/10. The formation of radial cracks on the corners indicates brittle material [157]. According to the mineral hardness conversion charts, the measured value is in good agreement with the tabular value of apatite (535 HV [158]).

Moreover, the sample does not show any luminescence when elucidated by long UV light (Fig. 16(a)). The specific gravity of the specimen was assessed by hydrostatic weighting (b) and pycnometric technique (c) to be 3.18 and 3.16 ± 0.20 g·cm−3, respectively. These values are in a good agreement with average density of apatite (3.19 g·cm−3, Chapter 1). All properties of investigated specimen mentioned above identify it as apatite, but the exact kind of apatite mineral and its chemical composition must still be determined yet.

Fig. 15.

The Vickers microhardness test with the load of 0.05 kgf for the time of 10 s and the correlation of results with the Mohs scale.

Fig. 16.

The specimen of apatite under UV light (compared with the fragment of red luminescence of alumina).

The sample was next treated to fine powder via milling in stain-less steal laboratory vibration mill. The apatite mineral was then determined by X-ray diffraction analysis (Fig. 17), mid-infrared spectroscopy (Fig. 18) and EBDS (Fig. 23) as fluorapatite (Ca5(PO4)3F, ref. [159]) with small amount (1%) of accessory mineral calcite

Since it is present as “free carbonate,” the sample cannot be considered as carbonate-fluorapatite (Section 2.6 and 4.6.1).

(CaCO3). Since the crystallographic parameters of identified hexagonal apatite are a = 9.3917, c = 6.8826 Å and Z = 2, it is possible to calculate the axial ratio (Eq. 11), the volume of cell (Eq. 12) and the density (Eq. 13) as follows:
a : c = 1 : 0.7328 E11
V = a 2 c sin 60 = 9.3917 2 6.8826 sin 60 = 525.74 A 3 E12
ρ calculated = M apatite Z V N A 1 × 10 24 = 504.31 × 2 525.74 6.025 × 10 23 1 10 24 = 3.19 g cm 3 E13

where M is the molecular weight of fluorapatite (Table 7 in Chapter 1) and NA is the Avogadro constant. The reconstruction of the cell of investigated apatite specimen is shown in Fig. 24.

3.2.2. X-ray diffraction analysis

Powder X-ray diffraction analysis of the apatite specimen (Fig. 17) identified it as fluorapatite. According to quantitative Rietveld analysis, the sample contains 99% of fluorapatite. There is also small amount (1%) of calcite

As the results of thermal analysis revealed (Section 3.2.4), the Rietveld analysis slightly overestimates the content of free carbonate in the sample.

that occurs on the surface of apatite specimen. Since the fluorapatite specimen (Fig. 14) is single crystal (Section 3.2.6), the crystal faces cannot be recognized, i.e. the crystal habit is anhedral (refer to Footnote 2 in Chapter 2). Nevertheless, XRD pattern shows that the most intensive diffraction possesses the Miller index (211), which corresponds to the Miller-Bravais indices (21–31), i.e. dihexagonal dipyramid (Chapter 1). There is also basal pinacoid {1000}, first-order {10–10} and second-order {11–20} hexagonal prism and dihexagonal prism {21–30}.

Fig. 17.

X-ray diffraction analysis of investigated specimen of apatite.

It is obvious that fluorapatite belongs to the hexagonal-dipyramidal crystal system, but the estimation of crystal habit of corresponding euhedral crystals from these results is highly speculative due to possible combination of pinacoid (c), first-order (m) and second-order hexagonal prisms (a) and first-order (p) and second-order dipyramids (s) with dihexagonal dipyramid (v) faces in the single crystal.

3.2.3. Infrared and Raman spectroscopy

Infrared (mid-FT-IR

Baseline corrected spectrum measured by KBr technique.

) and Raman spectrum of fluorapatite is shown in Fig. 18(a) and (b), respectively. The most expressive infrared bands are attributed to fundamental frequencies of tetrahedral phosphate ion [PO4]3−. The structure of apatite leads to the reduction of ion symmetry from Td (four fundamental frequencies with IR inactive ν1 mode) to Cs, where ν1 mode becomes IR active [128],[160],[161],[162],[163]:
  1. The ν1(PO4) mode is very weak (vw

    Abbreviation used for the expression of intensity and width of peak in the spectrum: very weak (vw), weak (w), middle (m), strong (s) and very strong (vs), shoulder (sh), broad (b), very broad (vb), and sharp (sp). Spectral bands related to impurities are abbreviated as imp [162].

    ) band that is related to symmetric stretching of phosphate ion.

  2. Bending: ν2 mode (vw).

  3. Anti-symmetric stretching: ν3 mode is the strongest (vs, vb) band in the infrared spectrum appearing in the spectral region from 1000 to 1150 cm−1.

  4. Bending: ν4 mode (m, vb) is observed between 540 and 620 cm−1.

Fig. 18.

Infrared (a) and Raman (b) spectrum of investigated specimen of fluorapatite.

The crystallinity of natural and synthetic apatite samples is often determined from the broadening of ν4(PO4) infrared absorption bands [160]. The assignment of bands in infrared and Raman spectrum of fluorapatite is listed in Table 2.

Vibration mode Assignment C6h factor group symmetry Raman shift IR
[cm−1]
ν1 Ag, E2g, E1u 956 965
ν2 E1g 432
E2g
Ag 449
E1u 460
Au 470
ν3 Au, E2g 1034 1032
E1g, E1u 1042 1040
Ag 1053
E2g 1061
Ag 1081
E1u 1090
ν4 Au 560
E1u 575
E2g 581
Ag, E1g 592
E1u 601
Ag 608
E2g 617

Table 2.

The interpretation of infrared and Raman bands in the spectrum of stoichiometric fluorapatite [164].

Factor group analysis of the hexagonal P63/M space group fluorapatite structure (Z = 2) yields an irreducible representation for optically active vibration of [165]:

Γ = 12 A g R + 7 A u I R + 8 E 1 g R + 11 E 1 u I R + 13 E 2 g R E14

where IR and R denote infrared and Raman activity, respectively. The influence of pressure on the infrared and Raman spectra of fluorapatite was investigated by Willams and Knittle [165]. Fluorapatite remains stable under pressures of at least 25 GPa at 300 K. Local environment of phosphate groups in fluorapatite becomes progressively less distorted from tetrahedral symmetry under the compression, as manifested by progressively smaller site-group.

The Davydov (factor group) splitting also decreases under the compression. This decrease is consistent with nondipole effects playing a primary role in the Davydov splitting of apatite; indeed, the magnitude of the Davydov splitting appears to be modulated by increases of the site symmetry of phosphate group under the compression [165].

The spectrum Fig. 18(a) shows weak peak located in the domain of OH stretching modes (from 3500 to 3600 cm−1) at the wave number of 3535 cm−1. This band belongs to the OH stretching mode in the hydrogen bond F…OH…(F) [166]. According to Freund and Knobel [167], the band at ~744 cm−1 belongs to the vibration of OH…F bond. According to Knubovets [168], the bands in the range from 745 to 720 cm−1 in apatite spectra could also be attributed to symmetric valence oscillations of the P-O-P bridge bonds, formed by the condensation of the PO43− tetrahedron. The presence of calcite causes that antisymmetric stretching mode (ν3) of planar CO32− ion appears in the infrared spectrum of investigated sample [169].

3.2.4. Thermal analysis

The results of simultaneous TG-DTA of investigated fluorapatite specimen are shown in Fig. 19. The sample is heated with the rate of 10°C·min−1 up to the temperature of 1425°C. The mass of sample is reduced by 1.25% during TG-DTA when the final temperature is reached. The most important features are the thermal decomposition of CaCO3 and the thermal decomposition of fluorapatite.

Fig. 19.

Thermal analysis of investigated specimen of fluorapatite.

The content of calcite was verified by thermal analysis. The weight of sample was reduced by 0.26% during the thermal decomposition of calcite. Since the theoretical mass loss of calcite is (100 × 44.09) / 100.086 = 43.97%,

It was calculated as 100 × molar mass (molar weight) of CO2/molar mass of CaCO3.

the content of calcite in the investigated sample of fluorapatite is (100 × 0.26) / 43.97 = 0.59%. Although this value is lower than the content of calcite determined by Rietveld analysis, there is still good agreement of both methods. The DTG curve shows that the process starts at the temperature of 565°C and wide of peak is of 135°C. The maximum rate decarbonation is reached at the temperature of 656°C.

At temperatures higher than 900°C, the weight of sample is reduced by the thermal decomposition of fluorapatite. The extrapolated beginning of the defluorination process (Section 8.6) was determined to be 1199°C. At the temperature of 1425°C, the defluorination process is still not complete. The extrapolation of experimental data

Theoretical mass loss (3.77%) of fluorapatite was used to set the fixed value of parameter A2.

shows that the thermal decomposition is most probably not complete before the temperature of melting point is reached (Table 7 in Chapter 1).

3.2.5. Scanning electron microscopy and WDX analysis

The microphotographs from SEM analysis of break plane of investigated fluorapatite specimen are shown in Fig. 20. The series of conchoidal fractures on the surface show brittle fracture as further characteristic properties of apatite.

Fig. 20.

SEM image of the break plane of investigated fluorapatite specimen: 50× (a), 75× (b), 1 000× (c) and 5 000× (d).

Fig. 21.

Typical results of WDX analysis of the fluorapatite specimen and simulation of electron–specimen interaction.

The chemical composition of sample was determined by WDX analysis (Fig. 21). The average results are introduced in Table 3.

Element Ca P O F Na
Weight [%] 38.64 ± 0.96 16.86 ± 0.47 40.55 ± 1.37 3.75 ± 0.08 0.23 ± 0.03
Ideal composition* 39.74 18.43 38.07 3.77
Atomic [%] 22.67 ± 0.85 12.81 ± 0.52 59.63 ± 1.31 4.65 ± 0.05 0.24 ± 0.03
Atomic ratio 4.91 ± 0.31 2.78 ± 0.18 12.77 ± 0.19 1 0.05 ± 0.01
Ideal composition

Resulting from the apatite formula (Section 1.4).

5 3 12 1

Table 3.

Average chemical composition of investigated specimen of fluorapatite.

The composition of fluorapatite (Ca5(PO4)3F → Ca5P3O12F) corresponds to the element ratio: Ca:P:O:F = 5:3:12:1. The results of WDX analysis of investigated specimen (Table 3) are in good agreement with the composition of ideal fluorapatite and fulfill all the criteria mentioned in Section 1.4.

3.2.6. EBDS analysis

The fragments of investigated fluorapatite specimens (Fig. 22) were further investigated by EBDS analysis.

Fig. 22.

Sample stage (a) and microphotograph of fluorapatite fragments used for EBDS (b).

The Kikuchi patterns confirm the sample as single crystal of fluorapatite Fig. 23.

Fig. 23.

EBDS analysis of investigated fluorapatite specimen.

Fig. 24.

The Kikuchi sphere showing the orientation of apatite crystal break plane.

The Kikuchi sphere in Fig. 24 shows the orientation of the sample that is in agreement with indistinct cleavage of fluorapatite to the direction given by Miller-Bravais indices of {0001}.

Advertisement

3.3. Thermodynamic properties of apatite

Drouet [170] gives the comprehensive review on experimental and predicted thermodynamic properties of phosphate apatites and oxyapatites, where O2− ion replaces 2X in general M10(PO4)6X2 formula of apatite phase and publishes the summary of available thermodynamic data including standard formation Gibbs energy (∆Gf°), ∆Hf° and S° at the temperature of 298 K (25°C) and the pressure of 1 bar (105 Pa), which are listed in the periodic table of phosphate apatites in Fig. 25.

Fig. 25.

Comprehensive periodic table of phosphate apatites provided by C. Drouet [170].

The discrepancies between published thermodynamic data probably arise from variable crystallinity states, polymorphs (either hexagonal or monoclinic, those not being systematically identified in literature reports), nonstoichiometry, hydration state and/or the presence of undetected impurities. When experimental-based data are not available (or are questionable), the so-called prediction of thermodynamic properties of solids becomes relevant. For example, it may allow an understanding of some unsuccessful experiment aiming at obtaining a desired hypothetical composition, or it may fill the gap between reported and needed thermodynamic values for the evaluation of equilibria constants or for the establishment of phase diagrams [170]. There are many methods developed for this purpose the summary of which can be found in works [171],[172].

For double oxides, AxByOz in the system AO-BO was established the dependence [171]:

Δ H f , cc = A B x O z y = f Δ H ¯ f E15

where ∆Hfcc°(AxByOz) is the standard enthalpy of the formation of double oxide AxByOz from the component oxides AO and BO and Δ H ¯ f represents the sum of molar fraction enthalpies of component oxides AO and BO according to the following relationship:

Δ H ¯ f = x A O Δ H f A O + x BO Δ H f BO E16

H°(AO) and ∆H°(BO) are the standard enthalpies of the formation of component oxides from the elements, and XAO and XBO are the molar fractions of component oxides in the double oxide AxByOz with a given composition [171].

The entropy of a solid compound is a function of masses of constituent atoms and the forces acting between these atoms: the greater the mass and the lower the force, the larger the entropy. The entropy of ionic solid will also depend upon the magnitude of the ionic charges. For compounds, the specific heat of which has reached the Dulong and Petit [173] value of 6 cal. per gram-atom [174],[175],

One calorie is 4.184 J (joules). Gram-atom [gm] (and gram-molecule) was used to specify the amount of chemical elements or compound. These units had a direct relation with “atomic weights” and “molecular weights,” which are in fact relative masses. “Atomic weights” were originally referred to the atomic weight of oxygen, by general agreement taken as 16. Although physicists separated the isotopes in a mass spectrometer and attributed the value of 16 to one of the isotopes of oxygen, chemists attributed the same value to the (slightly variable) mixture of isotopes 16, 17, and 18, which was for them naturally occurring element oxygen. Finally, an agreement between the International Union of Pure and Applied Physics (IUPAP) and the International Union of Pure and Applied Chemistry (IUPAC) brought this duality to an end in 1959/1960. Physicists and chemists have ever since agreed to assign the value 12, exactly, to the so-called atomic weight of the isotope of carbon with the mass number 12 (carbon 12, 12C), correctly called the relative atomic mass Ar(12C). The unified scale thus obtained gives the relative atomic and molecular masses, also known as the atomic and molecular weights, respectively [174]. The law is also known as Dulong and Petit principle, which can be expressed in modern unit as: atomic weigh × specific heat ≈ ∂(3kTNA)/∂T 3kNAcV ≈ 25 J·K−1mol−1, i.e. the atomic weight of solid element multiplied by its molar specific heat is a constant [175].

the mass is the principal factor, and in 1921, the authors gave an equation for the contribution of each element to the entropy of the compound [176].
S 298 K = 3 2 R ln at . wt . 0.94 E17

For simple salts, such as alkali halides, the entropy may be estimated with fair accuracy as the sum of the entropies of constituent elements as given by this equation. However, the forces in solid salts are largely the ionic attractions, and the effect of the ionic radii upon the force constants and the vibrational frequencies is appreciable; in general, the entropy of a large ion is increased and the entropy of a small ion is decreased compared to the values given by Eq. 17. [176].

3.3.1. Volume-based thermodynamic predictive method

The volume-based thermodynamic approach (VTB), the so-called first-order method, has especially received much attention because the method is rather easy to use and has been shown in some cases to lead to output data well related to experimental results [170].

3.3.2. Additive estimation methods

Additive estimation or contributive methods are probably the simplest approach based on the following [170],[177]:

  1. Atomic and ionic contribution: the technique based on the method proposed by Kellogg [178]:

    C pm A x B y , 298.15 K = x C pm ° A + y C pm ° B e . g . : C pm BaCl 2 , 298.15 K = C pm ° Ba + 2 C pm ° C l E18

    The approach was later revised by Kubaschewski [179],[180]. These authors also proposed the method for the estimation of parameters A, B and C in the temperature dependence of Cpm°(T)

    The full equation for the temperature dependence is Cpm° = A + BT + C/T2 + DT2 + F/T1/2 [J·K−1·mol−1].

    :
    C pm T = A + B T + C T 2 E19
    A = 10 3 T m C pm 298.15 K + 4.7 n 1.25 n 10 5 T m 2 9.05 n 10 3 T m 0.298 E20
    B = 25.6 n + 4.2 n 10 5 T m 2 C pm ° 298.15 K 10 3 T m 0.298 E21
    C = 4.2 n E22

    where n is the number of ions (contributions) in the formula unit. The described approach is worthy for the substances with melting point temperatures (Tm) bellow 2300 K. The data on cationic and anionic contributions to heat capacity at 298 K are published in works [177], [179],[181],[182],[183].

    For ionic compounds, the entropy can be calculated

    Table 4 [172] refers to the values of entropy contribution from work [176] where some data for PO43− ion are given in brackets. In order to verify this value (for the charge of cation, i.e. Ca2+, the contribution of PO43− anion is 17 [calories] × 4.184 = 71.13 J·K−1·mol−1, it is possible to calculate it from recommended thermodynamic data from Fig. 25. If (for example) three apatite end members were used, it is possible to calculate the contribution to PO43− anion in hydroxylapatite, fluorapatite, chlorapatite, and bromapatite as follows:

    H A P : 10 39.1 + 6 x + 2 18.83 = 780 x = 58.56 J K 1 mol 1 E23
    F A P : 10 39.1 + 6 x + 2 17.00 = 728 x = 50.50 J K 1 mol 1 E24
    ClAP : 10 39.1 + 6 x + 2 31.80 = 835 x = 63.40 J K 1 mol 1 E25
    BrAP : 10 39.1 + 6 x + 2 45.70 = 870 x = 64.60 J K 1 mol 1 E26

    It is also possible to calculate it from the contribution data for Ca3(PO4)2 or Mg3(PO4)2, where S°(298.15K) = 235.998 and 189.2 J·K−1mol−1 (HSC software v.7.1), respectively:

    3 39.1 + 2 x = 235.998 x = 59.35 J K 1 mol 1 E27
    3 23.40 + 2 x = 189.2 x = 59.5 J K 1 mol 1 E28

    Therefore, it is suggested to use average value from these calculations, i.e. PO43− (for M2+ cation) ≈ 59 J·K−1mol−1. It is then possible that the application of contribution techniques to apatite leads to positive error in estimated thermodynamic data.

    from additive data given in Table 4, empirically found for cation and anion constituents of the compound (increments method of Latimer [184]) [172].

    Element Cpm° Sm° Element Cpm° Sm° Element Cpm° Sm°
    Contribution of cations
    [J·K−1·mol−1]
    Ag 25.76 57.6 Hf 25.52 53,0 Pr 24.27 61.1
    Al 19.66 23.4 Hg 25.10 59.4 Rb 26.36 59.2
    As 25.10 45.2 Ho 23.01 56.0 Sb 23.85 58.9
    Au 58.5 In 24.27 55.0 Se 21.34 60.5
    B 23.5 Ir (23.85) 50.0 Si 35.2
    Ba 26.36 62.7 K 25.94 46.4 Sm 25.10 60.2
    Be (9.62) 12.6 La (25.52) 62.3 Sn 23.43 58.2
    Bi 26.78 65.0 Li 19.66 14.6 Sr 25.52 48.7
    Ca 24.69 39.1 Lu 51.5 Ta 23.01 53.8
    Cd 23.01 50.7 Mg 19.66 23.4 Te 69.0
    Ce 23.43 61.9 Mn 23.43 43.8 Th 25.52 59.9
    Co 28.03 34.1 Mo 35.9 Ti 21.76 39.3
    Cr 23.01 32.9 Na 25.94 37.2 Tl 27.61 72.1
    Cs 26.36 67.9 Nb 23.01 48.1 U 26.78 64.0
    Cu 25.10 44.0 Nd 24.27 60.7 V 22.18 36.8
    Fe 24.94 35.0 Ni (27.61) 35.1 Y (25.10) 50.4
    Ga (20.92) 40.0 Os 50.0 Zn 21.76 42.8
    Gd 23.43 56.0 P 14.23 39.5 Zr 23.85 37.2
    Ge 20.08 49.8 Pb 26.78 72.2
    Anion Cpm° Sm° for oxidation charge of cation
    I II III IV
    Contribution of anions
    J·K−1·mol−1]
    H 8.79
    F 22.80 20.8 17.0 18.3 20.4
    Cl 24.69 36.3 31.8 30.3 34.4
    Br 25.94 50.3 45.7 44.7 50.8
    I 26.36 58.3 53.5 54.8 53.9
    O2− 18.41 4.5 2.9 2.4 3.2
    S2− 24.48 20.6 18.4 20.1 17.0
    Se2− 26.78 35.5 32.8 34.1 30.9
    Te2− 27.20 38.1 41.9 44.1 40.1
    OH 30.96 (20.9) (18.8) (12.6)
    SO42− 76.57 80.0 69.5 64.2 (41.8)
    SO32− 42.9
    NO3 64.43 86.0 74.0
    NO2 70.6
    CO32− 58.58 62.4 46.6
    CrO42− 90.76
    MoO42− 90.37
    WO42− 97.49
    SiO32− 59.3 60.7 43.9 29.3
    SiO44− 73.5
    PO43− 73.90 (100.4) (71.1) (50.2)
    UO42− 107.11

    Table 4.

    The contribution of cations and anions to Cpm° and Sm° (298.15 K) [172],[180].

  2. Structural and simple oxide contribution: the value of Cpm°(298.15 K) or the parameters of the temperature dependence in Eq. 19 can be calculated from the contribution of constituent oxides (Neumann–Kopp rule, NKR):

    x A s + y B s = A x B y s , e . g . , S r O s + C u O s = SrCuO 2 s E23
    C pm A x B y , 298.15 K = x C pm ° A , 298.15 K + y C pm ° B , 298.15 K e . g . , C pm SrCuO 2 , 298.15 K = C pm ° S r O , 298.15 K + C pm ° C u O , 298.15 K E24

    or from the structural contribution [179].

  3. Prediction method for homological series and groups of chemically related substances (oxides): based on the approach of Aldabergenov et al [185] and Gospodinov and Mihov [186]. The molar heat capacity in homological series as Am(BxOy)n, is a linear function of n, i.e. the coefficient which specifies the number of complex anions (BxOy)z in the formula unit. For example, for the series of alkaline aluminate, it can be written as

    A l O 2 + A l O 2 A l 2 O 4 2 + A l O 2 A l 3 O 6 3 E25

    Since each higher anion is formed by the addition of primary ion (AlO2) unit, higher anion is considered to form n-multiples of primary ion, the value of which is determined from the available experimental data for KAlO2, LiAlO2 and NaAlO2 and from ions contribution for cation K+, Li+ and Na+ obtained from their standard entropies in an infinitely diluted solution [179].

Apatite phase may be treated in the first approximation as the sum of contributions arising from the constitution of binary oxides/compounds. For example, in the case of fluorapatite (Ca10(PO4)6F2), a decomposition into contribution of 9CaO + 3P2O5 + CaF2 could be considered. It can be generalized to any end-member in the form 9CaO + 3P2O5 + XF2 [170].

Advertisement

3.4. Dissolution of apatites

At the fundamental level, the reactions between solids and liquids involve a coupled sequence of mass transport, adsorption/desorption phenomena, heterogeneous reactions, chemical transformations of intermediates, etc., the identification, separation and kinetic quantification of which are all necessary if the mechanism of the process is to be fully understood and described [111],[187]. It was generally accepted that the process during the dissolution of lattice ions includes the following [187],[188]:

  1. Detachment of species (ion) from a kink site

  2. Surface diffusion from the crystal steps

  3. Desorption from the surface

  4. Diffusion into the bulk solution

The dissolution of apatites under steady-state conditions, in pure water or in aqueous acidic media, includes the following simultaneous steps [187],[189]:

  1. Diffusion of chemical reagents (H+) from bulk solution to the solid/liquid interface. In the case of acidic dissolution, the diffusion of acid anions (An−) must be taken into account as well.

  2. Adsorption of H+ onto the surface of apatite.

  3. Chemical transformation of the surface.

  4. Desorption of ions of fluoride, calcium and phosphate from the crystal surface.

  5. Adsorption of chemicals from solution back onto the surface of apatite.

These steps are likely to be more complicated, e.g. the processes I and V include chemical transformation of ionic species during diffusion because the pH of solution is known to depend on the distance from the solid/liquid interface. In other words, the value of pH is higher near to the surface of apatite and decreases with increasing distance from the surface.

When apatite gets in contact with undersaturated solution, the dissolution states from 1 to 5 mentioned above take place. In order to provide detailed description of the process, the following assumption must be introduced [187]:

  1. Stoichiometric apatite is dissolved, and neither nonstoichiometric layer nor other ions except for calcium, phosphate, hydroxide and fluoride are initially present in the crystal lattice, whereas the volume and surface defects (dislocations and dislocations outlets, respectively) might be present and, when present, they are distributed randomly.

  2. Except for the presence of dislocation outlets, the initial surface of apatite is perfect (molecularly sooth). Otherwise, each imperfection might be the dissolution nucleus.

  3. Despite the limitations and drawbacks discussed above, all models are correct and complementary to each other.

  4. In some cases, anions of acid might have an influence on apatite dissolution due to the specific affinity (e.g. citrate) by means of formation of insoluble compounds (e.g. sulfate). That is the reason why anions are not specified either here or below. It is just assumed that the dissolution of apatite happens in acid HnA, where An− is an anion.

  5. The hydration effect on all ions and molecules involved as well as that on crystal surface of apatite is omitted for simplicity.

  6. All crystal faces of solid apatite are equal and have similar ionic arrangement. There is not any considered specific influence of different crystal faces on the dissolution mechanism.

Since the dissolution models have limitations and drawbacks, none of them was able to describe the dissolution of apatite in general. Furthermore, the most of models were elaborated for the apatite dissolution in slightly acidic or nearly neutral solution (4 < pH < 8), for relatively small values of solution undersaturation and for the temperatures in the range from 25°C to 37°C, nothing is known about their validity for the dissolution of apatite in strong inorganic acid such as HCl, HNO3 and H2SO4 and at temperatures above 70°C [187].

The classification of congruent (stoichiometric)/incongruent (nonstoichiometric) dissolution is based on direct measurements of either ionic concentrations in the solution or the surface composition of apatite during the dissolution [187]:

  1. Congruent dissolution: ions in solids are dissolved simultaneously with the dissolution rates proportional to their molar concentrations, e.g. for Ca5(PO4)3Z it should be written [190],[191]:

    C a t P t P 0 R = 0 E26

    where [Ca]t, [P]t and [P]0 denotes actual (at time t) concentration of calcium and phosphorus in the solution and initial concentration of phosphorus (at time t = 0), respectively. The value of R is given by ideal stoichiometric ration of Ca:P = 5/3 in the formula of apatite (Table 7 in Chapter 1).

  2. Incongruent dissolution: the dissolution rate is different for each ion. That leads to the formation of surface layer with chemical composition different from the bulk of solid apatite phase.

The behavior of surface of apatite during the dissolution according to Dorozhkin [112] is shown in Fig. 26. Fluorine from fluorapatite or hydroxyl from hydroxyapatite dissolves most probably as the first. This can be explained by their position in the channels of crystal lattice. The dissolution starts with replacements of fluorine for water. Proton(s), chemisorbed on the nearest phosphate group(s), most probably catalyze this process. Local positive charge on apatite is formed as the result (b). Obtained local positive charge is removed by the detachment of one of the nearest calcium cations: Ca(2) is more likely to be detached first (c) since Ca(1) is located rather far from the channel. Acidic anions present in the solution most probably participate in this process. Later, proton(s) from the bulk solution replace other calcium cation(s) around the nearest phosphate group. Very thin surface layer of acidic calcium phosphates is formed as the result [112].

Fig. 26.

Schematic illustration of the surface dissolution mechanism of apatite at the nanolevel: (a) part of the initial surface of apatite; (b) replacement of one fluorine (or hydroxyl) anion with water molecule resulting in local positive charge formation; (c) removal of one of the nearest calcium cations; (d) sorption of next proton; (e) removal of another calcium cation with simultaneous formation of acidic calcium phosphate; (f) detachment of one phosphate anion together (or simultaneously) with third calcium cation. A jump-wise shift of the dissolution step occurs simultaneously at stage f. (●) Fluorine for fluorapatite or hydroxyl for hydroxylapatite; (○) Ca(II) on the first plane, (○) Ca(II) on the back plane and Ca(I) on the back plane; (*(+)) molecule of water and local positive charge; (∆) PO43− tetrahedra; H+∆H+ and ∆H+ represent the surface tetrahedral anions of H2PO4 and HPO42−, respectively. Chemisorbed protons, water molecules and acidic anions are omitted for simplicity. Note that crystal structure of apatite is shown very schematically: it should be hexagonal, while here it looks more or less like cubic [112].

When all (or almost all) the nearest calcium cations have been replaced with protons according to the reactions [112]:

CaHPO 4 + H + C a 2 + + H 2 P O 4 E27
CaHPO 4 + 2 H + C a 2 + + H 3 P O 4 E28

phosphate anions (H2PO4, CaH2PO4+, or H3PO4) also detach (f). As the result, the dissolution step moves forward jump-wise over a distance equal to the dimension of phosphate anion, of approximately 3 Å. The detachment of phosphate anions and calcium cations results in the formation of hole. The dimension of this hole should be close to the lattice parameters of apatite. Most probably, it is a dissolution nucleus on which the polynuclear dissolution mechanism is based [112].

3.4.1. Dissolution of fluorapatite

The adsorption of H+ onto the surface of apatite (Fig. 26(b)–(d)) resulted in the aqueous pH increasing from 5.60 to 8.45 within the first hour of dissolution. Ions of H+ were adsorbed onto oxygen ions of phosphate groups as well as onto ions of fluoride [187],[191].

The reaction of stoichiometric (congruent) dissolution of pure stoichiometric apatite can be expressed by the reaction

Reaction Eq. 7 and Eq. 12 are widely used for the description of dissolution process of apatite [112] using stoichiometry pertinent single or double apatite formula and Z = F or OH.

[112],[187],[192]:
C a 5 P O 4 3 Z 5 C a 2 + + 3 P O 4 3 + Z E29

where Z = OH and F. Assuming the unit activity of the solid phase (aCa5(PO4)3Z = 1), the equilibrium constant of dissolution (K) can be expressed via the solubility product or ion activity product (IAP) of apatite

Out of equilibrium state, IAP is not equal to K (see the discussion to Eq. 20). Double formula of apatite is assumed then with respect to the apatite stoichiometry; the law for ionic activity product has the following form:

Ksp = a(Ca2+)10 a(PO43−)6 a(Z−)2.

For example, in hydroxylapatite, where Z = OH, the activity of OH anion can be expressed by using ionic product of water (25°C): Kw = a(H+) a(OH) = 1·10−14 (mol·dm−3)2 and then a(OH) = Kw/ a(H+). Since pH = -log [H+] and pOH = -log [OH] the [H+] = 10−pH and [OH] = 10−pOH and pKw = pH + pOH = 14:

Ksp = a(Ca2+)10 a(PO43−)6 10−2pOH = a(Ca2+)10 a(PO43−)6 10−2(14−pH) = a(Ca2+)10 a(PO43−)6 Kw2102pH.

The activity of ionic species is the product of ion molar concentration ([196]) and ion activity coefficient, e.g. (a(Ca2+) = ([Ca2+] / [Ca2+]°) γ(Ca2+), where the standard state [Ca2+]° = 1 mol·dm−3 can be chosen), which can be calculated, e.g. the example via Debye–Hückel, extended Debye–Hückel, or modified Davies equation (in dependence on ionic strength):

log γi = -AZi2I (I < 10−3),

log γi = [-AZi2I / (1+B·αiI)] (10−3 < I < 0.1),

or log γi = [-AZi2I / (1+B·αiI)] + 0.3I (I > 0.1), respectively.

A and B are the temperature-dependent constants, Zi is the charge number of ith ion, αi is the radius of hydrated ith ion, and I is the number of ions: I = ½  ([Ci]Zi2).

The solubility product defined as the product of concentration of compound constituents ions, which are released during the dissociation is often used: Ksp = [Ca2+]10 [PO43−]6 [Z]2.

[187],[193],[194],[195],[196]:
K = I A P = a C a 2 + 5 a P O 4 3 3 a Z eq . E30

where ai denotes the thermodynamic activity of aqueous species. The standard Gibbs (free) energy of the reaction related to the standard temperature (298.15 K) and pressure (0.101 MPa) is given by the formula

rG° = -RT ln K (in the equilibrium state) and ln K = ln 10 log K. From that, it should be derived that ∆rG° = -8.314·298.15 ln 10 log K = -5708 log K [J·mol−1], where ln 10 ≈ 2.303. Since log10K = log K = lg K, it can also be written as ∆rG° = -5708 lg K.

:
Δ r G = 5.707 log K kJ mol 1 E31

For reaction 29, the following equation

rG° = ∑νiG°f,i, where νi denotes the stoichiometric coefficient of given species and ∆G°f,i its standard enthalpy of formation.

can be derived [187]:
Δ r G = 5 Δ G f C a 2 + + 3 Δ G f P O 4 3 + Δ G f Z Δ G f C a 5 P O 4 3 Z E32

Eq. 32 can be further treated as follows:

Δ G f C a 5 P O 4 3 Z = 5 Δ G f C a 2 + + 3 Δ G f P O 4 3 + Δ G f Z Δ r G E33

Harouiya et al [197] assumes that the dissolution of apatite in the temperature range from 5°C to 50°C, and the pH from 1 to 6 can be expressed by the following formula:

C a 5 P O 4 3 F + 3 H + 5 C a 2 + + 3 H P O 4 2 + F E34

With regard to assumed standard state, the equilibrium constant of reaction (Eq. 34) can be written as

Since the saturation of solution with respect to Ca5(PO4)3Z means that K = 0 (equilibrium state), it can be derived that a(Ca2+)5 a(PO43−)3 a(Z−)K´a(H+)3 = 0 and then ∆rG° = -RT ln (a(Ca2+)5 a(PO43−)3 a(Z−) / K´a(H+)3). Since A = -∆rG° (please see note 35), A = -RT ln (K´a(H+)3 / a(Ca2+)5 a(PO43−)3 a(Z−)).

K = a C a 2 + 5 a H P O 4 2 3 a F a H + 3 E35

The chemical affinity

The relationship between the reaction Gibbs energy and chemical affinity: A = -∆rG° was introduced by T. de Donder.

(A) of Eq. 34 is given by the law:
A = R T ln K a H + 3 a C a 2 + 5 a H P O 4 2 3 a F E36

In the closed-system experiment, the dissolution rates are generally obtained from the slope of concentration of reactive solution versus the time:

r = c i t M r ν i s E37

where r refers to the dissolution rate of apatite, ci denotes the concentration of ith element, t is the time, Mr designates the mass of fluid in the reactor, νI is the stoichiometric coefficient and s designates the total surface area of sample in the reactor. The slope of the plot may not be constant and may increase or decrease with time from the following reasons [195],[197]:

  1. Changes in the reactive fluid volume, which may occur due to the evaporation of solvent or regular sampling of reactive fluid

  2. Nonzero order reaction kinetics

  3. Approach to equilibrium, where the dissolution rate decreases and reaches zero at equilibrium. This approach is described by the transition state theory as follows:

r = r + 1 exp A σ R T E38

The symbol r+ symbolizes the far from equilibrium dissolution rate, which may depend on the composition of solution, A is the affinity of reaction of dissolution, σ stands for the Temkin’s average stoichiometric number equal to the ration of rate of destruction of the activated or precursor complex relative to overall rate, R designates universal gas constant and T is the temperature on the absolute scale. Overall rate (r) is equal to forward rate (r+) when A >> σ RT. As one of the approaches of equilibrium, overall rates gradually decrease and reach zero at equilibrium where A = 0. The value of r is within 10% of r+ when A/σ RT > 2.3 which is equivalent to A > 1.36 σ kcal·mol−1. It indicates that the parameter σ plays a crucial role in the variation of dissolution rates at near to equilibrium conditions [195],[197].

The value of r+ depends on the pH according to the following equation [197],[198]:

r + = k p H n E39

where k refers to the tare constant and n stands for the reaction order determined as the slope of linear dependence of ln r+ on pH. The dependence of k on the temperature is given by the Arrhenius law:

k = A A exp E A R T E40

where AA is the preexponential (frequency) factor and EA is the activation energy of the process. The combination of Eqs. 38, 39 and 40 leads to the equation for the dissolution rate as follows [197]:

r = A A a H + n exp E A / R T 1 exp A σ R T E41

For acidic dissolution of calcium fluorapatite, ions of F were found to dissolve faster (or prior to) when compared to calcium and phosphate. A similar phenomenon of prior (or faster) dissolution of calcium when compared to that of phosphate was also found [187],[199]. The release of calcium and phosphate ions from the surface of apatite seems to be affected by the presence of salts, such as Na2SO4, CH3COONa, or NaCl, in the solution. The concentration of phosphate in the solution increases in the following order [187]:

N a 2 SO 4 > C H 3 COON a > N a C 1 . E48

On the contrary, the concentration of calcium ions decreases in the same order.

The undersaturation (US) and relative undersaturation (USr) of apatite solvent dissolved upon is defined as follows

Since the system is not in the equilibrium state IAP ≠ K.

[193]:
US r = 1 US = 1 I A P K 1 / 18 E42

where K is the equilibrium constant of reaction 29 and IAP is the ion activity product. The law is written with regard to the stoichiometry of double formula of apatite, where ∑νi in Eq. 30 is 2 × (5 + 3 + 1) = 18. The value of IAP/K ratio is as follows:

  1. IAP/K > 1, the reaction proceeds to the left (precipitation, supersaturated solution)

  2. IAP/K = 1, the reaction is in the equilibrium state (saturated solution)

  3. IAP/K < 1, the reaction proceeds to the right (dissolution, undersaturated solution)

This ration is also used to calculate the saturation index (SI) for the reaction of dissolution [200]:

S I = log I A P K E43

Depending on the saturation index, the following states of solutions are recognized:

  1. SI < 0 undersaturation

  2. SI = 0 saturation, i.e. mineral

    The ranges of Si near zero are generally considered to be within the equilibrium zone for the mineral. The ranges of SI = 0 ± 0.5 and 0 ± (5%) (lg Kmineral) were used in various studies [200].

    or salt is in equilibrium with solution

  3. SI > 0 supersaturation

3.4.2. Classification of the dissolution models

The models, which are usually applied for the description of dissolution of apatites, include the following [187]:

  1. Diffusion and kinetically controlled models: the dissolution of apatite was found to be the diffusion (transport) controlled in some cases [201],[202], kinetically (surface) controlled in other ones [203],[204], and even intermediate [205], i.e. both diffusion and kinetically controlled. Both models usually operate with the so-called driving force which means either the concentration gradient within the Nernst diffusion layer (the diffusion controlled model) or the gradient of ionic chemical potentials between the apatite crystal surface and bulk solution (the kinetically controlled model). Moreover, the results obtained on these models are valid only within the experimental conditions studied; no extrapolation can be made beyond the tested ranges. For example, after, let´s say, a slight agitation decrease or temperature increase, an initially kinetically controlled dissolution might be controlled by the diffusion. Thus, high sensitivity to applied experimental conditions appears to be the main drawback of these models [187].

  2. Polynuclear model: is based on the study of dissolution and kinetics of growth of apatite under constant composition conditions [191],[193],[206],[207],[208]. Polydispersed samples of apatite were put into a stirred undersaturated (for dissolution experiments) or supersaturated (for those on crystal growth) solutions, and the pH of solution and the amount of added chemicals (an acid for the dissolution experiments and a base for those on crystal growth) were permanently recorded as the functions of time. The results obtained were plotted versus either undersaturation or supersaturation values: straight lines were obtained in the specific logarithmic coordinates typical for this model. According to the model, the dissolution nuclei, i.e., the collections of vacant sites for Ca2+, PO43−, and OH ions, are formed on the crystal surface of apatite and spread over the surface with a definite lateral rate.

    The nucleation rate is assumed as a function of mean ion activity. The lateral growth rate of nuclei is assumed proportional to the difference between total concentration of calcium ions in the saturated solution and in a solution, while the rate constant is related to the frequency for calcium ions to make a diffusion jump into a kink and, simultaneously, partly dehydrate. Recent investigation reveals that the rate-determining step was not the diffusion but two-dimensional surface nucleation [187].

  3. Self-inhibition model: assumes the formation of self-inhibition calcium-rich layer on the surface of apatite during the dissolution. According to this model, apatite is dissolved by ionic detachment of calcium and phosphate ions from the surface to a solution. When an initial portion of apatite has been dissolved, some amount of calcium cations (probably, in connection with anionic counter ions) is returned from the solution and adsorbed back onto the surface of apatite. The latter process results in the formation of a semipermeable ionic membrane, which is formed from positively charged layer containing strongly adsorbed calcium ions, i.e. calcium-rich layer [187], [201],[209],[210],[211].

  4. Congruent and incongruent dissolution model: was already described above.

  5. Chemical model: This model was developed from the self-evident supposition that it would be highly unlikely if apatite were dissolved by the detachment of “single molecules” equal to the unit cells and consisting of 18 ions. Moreover, in the crystal lattice, practically all ions are shared with neighboring unit cells and often cannot be attributed to given “single molecule.” Based on the experimental results obtained on one hand, and on analysis of the data found in references on the other hand, a sequence of four successive chemical reactions was proposed to describe the process of apatite dissolution [112],[187],[212],[213]:

C a 5 P O 4 3 F , O H + H 2 O + H + C a 5 P O 4 3 H 2 O + + H F , H 2 O E44
2 C a 5 P O 4 3 H 2 O + 3 C a 3 P O 4 2 + C a 2 + + 2 H 2 O E45
C a 3 P O 4 2 + 2 H + C a 2 + + 2 CaHPO 4 E46
CaHPO 4 + H + C a 2 + + H 2 P O 4 E47

Eqs. 44–47 can be used instead well-known net reactions [187]:

C a 5 P O 4 3 F , O H + H 2 O + 7 H + 5 C a 2 + + 3 H P O 4 + H F , H 2 O E48
C a 5 P O 4 3 F , O H 5 C a 2 + + 3 P O 4 3 + F , O H E49

In principle, the dissolution process could also happen according to reaction 49 followed by chemical interaction in the solution among ions of apatite and acid near the crystal surface [214]:

5 C a 2 + + 3 P O 4 3 + F , O H + 7 H + 5 C a 2 + + 3 H 2 P O 4 + H F , H 2 O E50
  1. Etch pit formation: the process of etch pit formation describes the dissolution of apatite crystals containing structural defects (dislocations and inclusions). The presence of dislocations accelerates the dissolution because they give rise to continuous steps on the surface and the strain energy, which they cause in crystals, favors the etch pit formation. The pits appear at the dislocation outlets; they are 0.1–10 pm in size (the dimensions depend on the dissolution kinetics and dissolution time: they increase when the dissolution progresses) and usually have a hexagonal shape according to the crystal symmetry P63/M of pure apatite [187],[215],[216],[217],[218],[219].

  2. Ion exchange model: based on a supposition about the adsorption of protons and anions of acid (e.g. citrate anions) from a solution onto the surface of apatite and removing of ions of calcium and phosphate into the solution instead [187].

  3. Hydrogen catalytic model: based on a reasonable suggestion about the adsorption of protons onto negatively charged oxygen ions of phosphate groups of apatite. The sorption of protons results in the transformation of surface PO43− groups into HPO42− and catalyzes the dissolution process [187].

3.4.3 Methods for the evaluation of reactivity of phosphate rocks

Chemical methods are used for the evaluation of reactivity of different phosphate rocks from which the fertilizers are manufactured for their possible direct application as fertilizers via empirical solubility test. Citric acid, formic acid, neutral ammonium citrate,

Neutral ammonium citrate is prepared by dissolving required amount of citric acid and neutralizing it with ammonium hydroxide. The pH of the reagent is adjusted to neutral [221].

and alkaline ammonium citrate are used as solvents for the extraction of P2O5. The latter is used mainly for the evaluation of calcined aluminum phosphates. Most of these reagents were not originally intended to evaluate the reactivity of phosphate rocks. For instance, neutral and alkaline ammonium citrate solutions were originally intended to separate chemical reaction products in superphosphate and other fertilizers from unreacted rock on the assumption that unreacted rock was insoluble in these reagents. The citric acid extraction was developed to evaluate basic slag, a popular fertilizer material in European countries. The formic acid extraction was developed specifically for phosphate rocks [220].

Nearly all extraction methods use the ratio of sample weight to extraction volume 1 g:100 ml.

The amount of used solution is also expressed in the name of the method, e.g. 100 ml method or 150 ml method [222].

The extraction time usually ranges from 30 min to 1 hour. The temperature and the agitation during extraction test may be specified. For example, the AOAC method uses neutral ammonium citrate

The neutral ammonium citrate test was used as the official method in the United States, and the test by acidic acid was developed for the comparison [222].

of specified concentration (1 g of sample and 100 ml of solution) with the extraction time of 30 min at 65°C. The Wagner method uses 2% solution of citric acid, the extraction time of 30 min and the temperature of 17.5°C [220],[221],[222],[223]. Unavailable phosphoric acid is usually expressed as the portion of fertilizer, which is insoluble in neutral ammonium citrate [769]. Phosphate removed during the neutral ammonium test is termed as citrate-soluble. The sum of water-soluble and citrate-soluble phosphate is termed as available [224].

Fused magnesium phosphate (FMP) is highly soluble in 2% citric acid

It is believed that weak citric acid solution imitates the condition near the plant roots [222].

but is less soluble in neutral citrate, while calcined defluorinated phosphate (CDP) and Thomas slag (Thomas phosphate) are fairly soluble in both citric acid and citrate. Actually, FMP dissolves fairly rapidly in neutral citrate at the beginning, but the dissolution is hindered by gelatinous silica, which forms on the surface of the FMP particles. This layer can be broken by vigorous stirring [222].

One disadvantage of all these methods is that the percentage of leached P2O5 depends on the grade of the rock, especially when the rock contains inert gangue minerals such as silica. In order to eliminate the adventitious effect of grade, the concept of absolute citrate solubility index (ASC) was developed [225],[226]:

A S C = AOAC citrate solubility P 2 O 5 % Theoretical P 2 O 5 in apatite % E51

The percentage of dissolved P2O5 is expressed as the gangue-free apatite [220]. If rocks contain free calcium of magnesium carbonate, these carbonates should be removed by the extraction with a suitable reagent before carrying out the test in order to obtain correct indication of reactivity [227].

It was also found that the length of the a-axis in the apatite unit cell (a0) is statistically related to the ASC according to the relationship [225]:

A S C = 421.4 9.369 a 0 E52

References

  1. 1. Straaten PV. Rocks for Crops, Agro Minerals of Sub-Sahara Africa. (CD), ICRAF, Nairobi, Kenya, 2002. ISBN: 0-88955-512-5.
  2. 2. Abouzeid Abdel-Z.M. Physical and thermal treatment of phosphate ores—an overview. International Journal of Mineral Processing 2008;85(4) 59–84.
  3. 3. Comodi P, Liu Y, Zanazzi PF, Montagnoli M. Structural and vibrational behaviour of fluorapatite with pressure. Part 1: in situ single-crystal X-ray diffraction investigation. Physics and Chemistry of Minerals 2001;28(4) 219–224.
  4. 4. Wilson RM, Elliott JC, Dowker SEP, Smith RI. Rietveld structure refinement of precipitated carbonate apatite using neutron diffraction data. Biomaterials 2004;25(11) 2205–2213.
  5. 5. Smith F. Industrial Applications of X-Ray Diffraction. CRC Press, 1999. ISBN: 978-0824719920
  6. 6. Elliott JC. Structure and Chemistry of the Apatites and Other Calcium Orthophosphates, Elsevier Press, Amsterdam, 1994. ISBN: 978-0-444-81582-8
  7. 7. Svergun DI, Koch MHJ, Timmins PA, May RP. Small Angle X-Ray and Neutron Scattering from Solutions of Biological Macromolecules. IUCr (International Union of Crystallography) texts on crystallography—Volume 19. Oxford University Press, 2013. ISBN: 978-0199639533
  8. 8. Nitske WR. The Life of Wilhelm Conrad Röntgen, Discoverer of the X ray. University of Arizona Press, 1971. ISBN: 978-0816502592
  9. 9. Daubrée A. Experiences sur la production artficielle de l’apatite, de la topaze, at de quelques autres métax fruorifés. Comptes rendus de l’Académie des Sciences 1851;32 625.
  10. 10. Röntgen WC. On a new kind of rays. Nature 1896;53 274–276. doi: 10.1038/053274b0
  11. 11. Weber AS. Nineteenth-Century Science: An Anthology. Broadview Press, 2000. ISBN: 978-1551111650
  12. 12. Rosenthal-Schneider I. Begegnungen mit Einstein, von Laue und Planck: Realität und wissenschaftliche Wahrheit, Facetten der Physik—Volume 12. Springer-Verlag, 2013. ISBN: 978-3663018841
  13. 13. Chatterji T. Neutron Scattering from Magnetic Materials. Gulf Professional Publishing, 2005. ISBN: 978-0080457055
  14. 14. Liang L, Rinaldi R, Schober H. Neutron Applications in Earth, Energy and Environmental Sciences. Neutron Scattering Applications and Techniques. Springer Science & Business Media, 2008. ISBN: 978-0387094168
  15. 15. Kolb U, Shankland K, Meshi L, Avilov A, David WIF. Uniting Electron Crystallography and Powder Diffraction, NATO Science for Peace and Security Series B: Physics and Biophysics. Springer, 2012. ISBN: 978-9400755802
  16. 16. Fitter J, Gutberlet T, Katsaras J. Neutron Scattering in Biology: Techniques and Applications. Biological and Medical Physics, Biomedical Engineering. Springer Science & Business Media, 2006. ISBN: 978-3540291114
  17. 17. Guillot S, Beaudet-Savignat S, Lambert S, Vannier Rose-N, Roussel P, Porcher F. Evidence of local defects in the oxygen excess apatite La9.67(SiO4)6O2.5 from high resolution neutron powder diffraction. Journal of Solid State Chemistry 2009;182(12) 3358–3364.
  18. 18. Lambert S, Vincent A, Bruneton E, Beaudet-Savignat S, Guillet F, Minot B, Bouree F. Structural investigation of La9.33Si6O26- and La9AESi6O26+δ-doped apatites-type lanthanum silicate (AE = Ba, Sr and Ca) by neutron powder diffraction. Journal of Solid State Chemistry 2006;179(8) 2602–2608.
  19. 19. Sansom JHJ, Richings D, Slater PR. A powder neutron diffraction study of the oxide-ion-conducting apatite-type phases, La9.33Si6O26 and La8Sr2Si6O26. Solid State Ionics 2001;139(3–4) 205–210
  20. 20. Matsushita Y, Izumi F, Kobayashi K, Igawa N, Kitazawa H, Oyama Y, Miyoshi S, Yamaguchi S. Powder neutron diffraction of La-apatite under low temperature. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment 2009;600(1) 319–321.
  21. 21. Tolchard JR, Slater PR. A high temperature powder neutron diffraction structural study of the apatite-type oxide ion conductor, La9.67Si6O26.5. Journal of Physics and Chemistry of Solids 2008;69(10) 2433–2439
  22. 22. Mazza D, Tribaudino M, Delmastro A, Lebech B. Synthesis and Neutron Diffraction Study of La5Si2BO13, an Analog of the Apatite Mineral. Journal of Solid State Chemistry 2000;155(2) 389–393.
  23. 23. O’Donnell MD, Hill RG, Fong SK. Neutron diffraction of chlorine substituted fluorapatite. Materials Letters 2009;63(15) 1347–1349.
  24. 24. Masubuchi Y, Higuchi M, Takeda T, Kikkawa S. Oxide ion conduction mechanism in RE9.33(SiO4)6O2 and Sr2RE8(SiO4)6O2 (RE = La, Nd) from neutron powder diffraction. Solid State Ionics 2006;177(3–4) 263–268.
  25. 25. Leventouri Th, Bunaciu CE, Perdikatsis V. Neutron powder diffraction studies of silicon-substituted hydroxyapatite. Biomaterials 2003;24(23) 4205–4211.
  26. 26. Alonso JA, Martínez-Lope MJ, Aguadero A, Daza L. Neutron powder diffraction as a characterization tool of solid oxide fuel cell materials. Progress in Solid State Chemistry 2008;36(1–2) 134–150.
  27. 27. Loong C-K, Rey C, Kuhn LT, Combes C, Wu Y, Chen S-H, Glimcher MJ. Evidence of hydroxyl-ion deficiency in bone apatites: an inelastic neutron-scattering study. Bone 2000;26(6) 599–602.
  28. 28. Young RA, Spooner S. Neutron diffraction studies of human tooth enamel. Archives of Oral Biology 1970;15(1) 47–63.
  29. 29. Ibrahim DM, Mostafa AA, Korowash SI. Chemical characterization of some substituted hydroxyapatites. Chemistry Central Journal 2011;5(1) 74.
  30. 30. Lachance GR, Claisse F. Quantitative X-ray fluorescence analysis: theory and application. Wiley, 1995. ISBN: 978-0471951674
  31. 31. Muller R. Spectrochemical Analysis by X-Ray Fluorescence. Springer Science & Business Media, 2012. ISBN: 978-1468417975
  32. 32. Crompton TR. Polymer Reference Book. Smithers Rapra Publishing, 2006. ISBN: 978-1859574928
  33. 33. Nadkarni RA. Modern Instrumental Methods of Elemental Analysis of Petroleum Products and Lubricants, ASTM special technical publication—Volume 1109. ASTM (American Society for Testing and Materials) International, 1991. ISBN: 978-0803114166
  34. 34. Tang YH, Han CM, Bao ZK, Huang YY, He W, Hua W. Analysis of apatite crystals and their fluid inclusions by synchrotron radiation X-ray flourescence microprobe. Spectrochimica Acta Part B: Atomic Spectroscopy 2005;60(4) 439–446.
  35. 35. Uda M, Demortier G, Nakai I. X-rays for Archaeology. NATO Science Series: Mathematics, Physics, and Chemistry—Volume 224. Springer Science & Business Media, 2005. ISBN: 978-1402035814
  36. 36. Nutman AP, Friend CRL. Petrography and geochemistry of apatites in banded iron formation, Akilia, W. Greenland: consequences for oldest life evidence. Precambrian Research 2006;147(1–2) 100–106.
  37. 37. Chen Ch, Chai Z, Gao Y. Nuclear Analytical Techniques for Metallomics and Metalloproteomics. Royal Society of Chemistry, 2010. ISBN: 978-1847559012
  38. 38. Bunker G. Introduction to XAFS: A Practical Guide to X-ray Absorption Fine Structure Spectroscopy. Cambridge University Press, 2010. ISBN: 978-0521767750
  39. 39. Parat F, Dungan MA, Streck MJ. Anhydrite, pyrrhotite, and sulfur-rich apatite: tracing the sulfur evolution of an Oligocene andesite (Eagle Mountain, CO, USA). Lithos 2002;64(3–4) 63–75.
  40. 40. Klemm W. Anomale Wertigkeiten. Angewandte Chemie 1951;63(17–18) 396–402.
  41. 41. Chauhan P, Chauhan RP, Gupta M. Estimation of naturally occurring radionuclides in fertilizers using gamma spectrometry and elemental analysis by XRF and XRD techniques. Microchemical Journal 2013;106; 73–78.
  42. 42. Hasikova J, Sokolov A, Titov V, Dirba A. On-line XRF analysis of phosphate materials at various stages of processing. Procedia Engineering 2014;83 455–461.
  43. 43. Zhou F, Wang L, Xu Z, Liu Q, Chi R. Reactive oily bubble technology for flotation of apatite, dolomite and quartz. International Journal of Mineral Processing 2015;134 74–81.
  44. 44. Oliveira MS, Santana RC, Ataíde CH, Barrozo MAS. Recovery of apatite from flotation tailings. Separation and Purification Technology 2011;79(1) 79–84.
  45. 45. Taylor HE. Inductively Coupled Plasma–Mass Spectrometry: Practices and Techniques. Academic Press, 2001. ISBN: 978-0126838657
  46. 46. Montaser A. Inductively Coupled Plasma Mass Spectrometry. John Wiley & Sons, 1998. ISBN: 978-0471186205
  47. 47. Levinson R. More Modern Chemical Techniques. Royal Society of Chemistry, 2001. ISBN: 978-0854049295
  48. 48. Lajunen LHJ, Perämäki P. Spectrochemical Analysis by Atomic Absorption and Emission. Royal Society of Chemistry, 2004. ISBN: 978-0854046249
  49. 49. Hasebe N, Barbarand J, Jarvis K, Carter A, Hurford AJ. Apatite fission-track chronometry using laser ablation ICP-MS. Chemical Geology 2004;207(3–4) 135–145.
  50. 50. Johnstone S, Hourigan J, Gallagher Ch. LA-ICP-MS depth profile analysis of apatite: Protocol and implications for (U–Th)/He thermochronometry. Geochimica et Cosmochimica Acta 2013;109 143–161.
  51. 51. Guzmics T, Zajacz Z, Kodolányi J, Halter W, Szabó C. LA-ICP-MS study of apatite- and K feldspar-hosted primary carbonatite melt inclusions in clinopyroxenite xenoliths from lamprophyres, Hungary: Implications for significance of carbonatite melts in the Earth’s mantle. Geochimica et Cosmochimica Acta 2008;72(7) 1864–1886.
  52. 52. Chew DM, Sylvester PJ, Tubrett MN. U–Pb and Th–Pb dating of apatite by LA-ICPMS. Chemical Geology 2011;280(1–2) 200–216.
  53. 53. Donelick RA, O’Sullivan PB, Ketcham RA, Hendriks BWH, Redfield TF. Relative U and Th concentrations from LA-ICP-MS for apatite fission-track grain-age dating. Geochimica et Cosmochimica Acta 2006;70(18) A143.
  54. 54. Fisher ChM, McFarlane ChRM, Hanchar JM, Schmitz MD, Sylvester PJ, Lam R, Longerich HP. Sm-Nd isotope systematics by laser ablation-multicollector-inductively coupled plasma mass spectrometry: methods and potential natural and synthetic reference materials. Chemical Geology 2001;284(1–2) 1–20.
  55. 55. Yang Yue-H Wu Fu-Y, Xie Lie-W, Chu Zhu-Y, Yang Jin-H. Re-evaluation of interferences of doubly charged ions of heavy rare earth elements on Sr isotopic analysis using multi-collector inductively coupled plasma mass spectrometry. Spectrochimica Acta Part B: Atomic Spectroscopy 2014;97 118–123.
  56. 56. Yang Yue-H Wu Fu-Y, Chu Zhu-Y, Xie Lie-W, Yang Jin-H. High-precision simultaneous determination of 147Sm/144Nd and 143Nd/144Nd ratios in Sm-Nd mixtures using multi-collector inductively coupled plasma mass spectrometry and its comparison to isotope dilution analysis. Spectrochimica Acta Part B: Atomic Spectroscopy 2013;79–80 82–87.
  57. 57. Pourmand A, Dauphas N, Ireland TJ. A novel extraction chromatography and MC-ICP-MS technique for rapid analysis of REE, Sc and Y: revising CI-chondrite and Post-Archean Australian Shale (PAAS) abundances. Chemical Geology 2012;291 38–54.
  58. 58. Boutton T. Mass Spectrometry of Soils (Books in Soils, Plants, and the Environment). CRC Press, 1996. ISBN: 978-0824796990
  59. 59. Casas JS, Sordo J. Lead: Chemistry, analytical aspects, environmental impact and health effects. Elsevier, 2011. ISBN: 978-0080463889
  60. 60. Becker S. Inorganic Mass Spectrometry: Principles and Applications. John Wiley & Sons, 2008. ISBN: 978-0470517208
  61. 61. Hellebrand SJE. Textural, geochronological and chemical constraints from polygenetic titanite and monogenetic apatite from a mid-crustal shear zone: an integrated EPMA, SIMS, and TIMS study. Chemical Geology 2007;241(1–2) 88–107.
  62. 62. Schoene B, Latkoczy C, Schaltegger U, Günther D. A new method integrating high-precision U–Pb geochronology with zircon trace element analysis (U–Pb TIMS-TEA). Geochimica et Cosmochimica Acta 2010;74(24) 7144–7159.
  63. 63. Teixeira W, D’Agrella-Filho MS, Hamilton MA, Ernst RE, Girardi VAV, Mazzucchelli M, Bettencourt JS. U–Pb (ID-TIMS) baddeleyite ages and paleomagnetism of 1.79 and 1.59 Ga tholeiitic dyke swarms, and position of the Rio de la Plata Craton within the Columbia supercontinent. Lithos 2013;174 157–174.
  64. 64. Metcalfe I, Crowley JL, Nicoll RS, Schmitz M. High-precision U-Pb CA-TIMS calibration of middle Permian to lower Triassic sequences, mass extinction and extreme climate-change in eastern Australian Gondwana. Gondwana Research, In Press, Corrected Proof, Available online 16 September 2014. doi:10.1016/j.gr.2014.09.002
  65. 65. Yang Yue-H Wu Fu-Y, Yang Jin-H, Chew DM, Xie Lie-W, Chu Zhu-Y, Zhang Yan-B, Huang Ch. Sr and Nd isotopic compositions of apatite reference materials used in U–Th–Pb geochronology. Chemical Geology 2014;385 35–55.
  66. 66. Cherepin VT. Secondary Ion Mass Spectroscopy of Solid Surfaces. VSP, 1987. ISBN: 978-9067640787
  67. 67. O’Connor J, Sexton BA, Smart RSC. Surface Analysis Methods in Materials Science Physics and Astronomy. Springer Series in Surface Sciences—Volume 23. Springer Science & Business Media, 2003. ISBN: 978-3540413301
  68. 68. van der Heide P. Secondary Ion Mass Spectrometry: An Introduction to Principles and Practices. John Wiley & Sons, 2014. ISBN: 978-1118916773
  69. 69. Kuo J. Electron Microscopy: Methods and Protocols—Volume 369. Springer Science & Business Media, 2007. ISBN: 978-1597452946
  70. 70. Microelectronics Failure Analysis: Desk Reference. ASM International, 2004. ISBN: 978-0871708045
  71. 71. Watson LL, Hutcheon ID, Epstein S, Stolper EM, Water on Mars: clues from deuterium/hydrogen andwater contents of hydrous phases in SNCmeteorites. Science 1994;265 86–90.
  72. 72. Barnes JJ, Franchi IA, Anand M, Tartèse R, Starkey NA, Koike M, Sano Y, Russell SS. Accurate and precise measurements of the D/H ratio and hydroxyl content in lunar apatites using NanoSIMS. Chemical Geology 2013;337–338 48–55.
  73. 73. Gotliv Bat-A, Robach JS, Veis A. The composition and structure of bovine peritubular dentin: mapping by time of flight secondary ion mass spectroscopy. Journal of Structural Biology 2006;156(2) 320–333.
  74. 74. Stephan T. TOF-SIMS in cosmochemistry. Planetary and Space Science 200149(1) 859–906.
  75. 75. Li Qiu-L, Li Xian-H, Wu Fu-Y, Yin Qing-Z, Ye Hai-M, Liu Y, Tang Guo-Q, Zhang Chuan-L. In-situ SIMS U-Pb dating of phanerozoic apatite with low U and high common Pb. Gondwana Research 2012;21(4) 745–756.
  76. 76. Schaltegger U, Schmitt AK, Horstwood MSA. U-Th-Pb zircon geochronology by ID-TIMS, SIMS, and laser ablation ICP-MS: recipes, interpretations, and opportunities. Chemical Geology 2015;402 89–110.
  77. 77. Ling Xiao-X, Schmädicke E, Li Qiu-L, Gose J, Wu Rui-H, Wang Shi-Q, Liu Y, Tang Guo-Q, Li Xian-H. Age determination of nephrite by in-situ SIMS U–Pb dating syngenetic titanite: a case study of the nephrite deposit from Luanchuan, Henan, China. Lithos 2015;220–223 289–299.
  78. 78. Friedbacher G, Bubert H. Surface and Thin Film Analysis: A Compendium of Principles, Instrumentation, and Applications. John Wiley & Sons, 2011. ISBN: 978-3527636945
  79. 79. Lee MS. Mass Spectrometry Handbook. Wiley Series on Pharmaceutical Science and Biotechnology: Practices, Applications and Methods—Volume 7. John Wiley & Sons, 2012. ISBN: 978-1118180723
  80. 80. Dambach S, Fartmann M, Kriegeskotte C, Brüning C, Wiesmann HP, Lipinsky D, Arlinghaus HF. Laser-SNMS analysis of apatite formation in vitro. Applied Surface Science 2004;231–232 506–509.
  81. 81. Rabi II, Zacharias S, Millman P. Kusch. A new method of measuring nuclear magnetic moment. Physical Review 1938;53(4) 318.
  82. 82. Weil JA, Bolton JR. Electron Paramagnetic Resonance: Elementary Theory and Practical Applications. 2nd ed., John Wiley & Sons, 2007. ISBN: 978-0470084984
  83. 83. Brustolon M. Electron Paramagnetic Resonance: A Practitioner’s Toolkit. John Wiley & Sons, 2009. ISBN: 978-0470432228
  84. 84. Mabbs FE, Collison D. Electron Paramagnetic Resonance of d Transition Metal Compounds Studies in Inorganic Chemistry. Elsevier, 2013. ISBN: 978-1483291499
  85. 85. Misra SK. Handbook of Multifrequency Electron Paramagnetic Resonance: Data and Techniques. John Wiley & Sons, 2014. ISBN: 978-3527672455
  86. 86. Moens PDW, Callens FJ, Verbeeck RMH, Naessens DE. An EPR spectrum decomposition study of precipitated carbonated apatites (NCAP) dried at 25°C: adsorption of molecules from the atmosphere on the apatite powders. Applied Radiation and Isotopes 1993;44(1–2) 279–285.
  87. 87. de Oliveira LM, Rossi AM, Lopes RT. Dose response of A-type carbonated apatites prepared under different conditions. Radiation Physics and Chemistry 2001;61(3–6) 485–487.
  88. 88. Schramm DU, Rossi AM. Electron spin resonance (ESR) studies of CO2 radicals in irradiated A and B-type carbonate-containing apatites. Applied Radiation and Isotopes 2000;52(5) 1085–1091.
  89. 89. Vugman NV, Rossi AM, Rigby SEJ. EPR dating CO2 sites in tooth enamel apatites by ENDOR and triple resonance. Applied Radiation and Isotopes 1995;46(5) 311–315.
  90. 90. Fattibene P, Callens F. EPR dosimetry with tooth enamel: a review. Applied Radiation and Isotopes 2010;68(11) 2033–2116.
  91. 91. Rudko VV, Ishchenko SS, Vorona IP, Baran NP. Isotropic radical image in biological apatites. Radiation Measurements 2007;42(9) 1580–1582.
  92. 92. Ishchenko SS, Vorona IP, Okulov SM, Baran NP. 13C hyperfine interactions of CO2 in irradiated tooth enamel as studied by EPR. Applied Radiation and Isotopes 2002;56(6) 815–819.
  93. 93. Nosenko VV, Vorona IP, Baran NP, Ishchenko SS, Vysotskyi BV, Krakhmalnaya TV, Semenov YA. Comparative EPR study CO2 radicals in modern and fossil tooth enamel. Radiation Measurements 2015;78 53–57.
  94. 94. Rudko VV, Vorona IP, Baran NP, Ishchenko SS. Thermally stimulated transformation of the EPR spectra in γ-irradiated bone tissue. Radiation Measurements 2009;44(3) 239–242.
  95. 95. Ishchenko SS, Vorona IP, Okulov SM, Baran NP, Rudko VV. ENDOR study of CO2 radicals in hydroxyapatite of γ-irradiated bone. Radiation Measurements 2011;46(1) 37–39
  96. 96. Köseogˇlu R, Köseogˇlu E, Köksal F, Başaran E, Demirci D. EPR of some irradiated renal stones. Radiation Measurements 2005;40(1) 65–68.
  97. 97. Fleet ME. Infrared spectra of carbonate apatites: ν2-Region bands. Biomaterials 2009;30(8) 1473–1481.
  98. 98. Yu H, Zhang H, Wang X, Gu Z, Li X, Deng F. Local structure of hydroxy–peroxy apatite: a combined XRD, FT-IR, Raman, SEM, and solid-state NMR study. Journal of Physics and Chemistry of Solids 2007;68(10) 1863–1871.
  99. 99. Klinowski J. New Techniques in Solid-State NMR, Volume 246. Springer Science & Business Media, 2004. ISBN: 978-3540221685
  100. 100. Duer MJ. Solid State NMR Spectroscopy: Principles and Applications. John Wiley & Sons, 2008. ISBN: 978-0470999387
  101. 101. Ohannesian L, Streeter A. Handbook of Pharmaceutical Analysis Drugs and the Pharmaceutical Sciences. CRC Press, 2001. ISBN: 978-0824741945
  102. 102. Bakhmutov VI. Solid-State NMR in Materials Science: Principles and Applications. CRC Press, 2011. ISBN: 978-1439869635
  103. 103. Burnell E, de Lange CA. NMR of Ordered Liquids. Springer Science & Business Media, 2003. ISBN: 978-1402013430
  104. 104. Braun M, Jana C. 19F NMR spectroscopy of fluoridated apatites. Chemical Physics Letters 1995;245(1) 19–22.
  105. 105. McElderry John-DP, Zhu P, Mroue KH, Xu J, Pavan B, Fang M, Zhao G, McNerny E, Kohn DH, Franceschi RT, Holl MMB, Tecklenburg MMJ, Ramamoorthy A, Morris MD. Crystallinity and compositional changes in carbonated apatites: Evidence from 31P solid-state NMR, Raman, and AFM analysis. Journal of Solid State Chemistry 2013;206 192–198.
  106. 106. Kaflak A, Kolodziejski W. Phosphorus-31 spin-lattice NMR relaxation in bone apatite and its mineral standards. Solid State Nuclear Magnetic Resonance 2007;31(4) 174–183.
  107. 107. Yesinowski JP, Eckert H. Hydrogen environments in calcium phosphates: proton MAS NMR at high spinning speeds. Journal of the American Chemical Society 1987;109(21) 6274–6282.
  108. 108. Goldstein J, Newbury DE, Joy DC, Lyman ChE, Echlin P, Lifshin E, Sawyer L, Michael JR. Scanning Electron Microscopy and X-ray Microanalysis: Third Edition. Springer Science & Business Media, 2012. ISBN: 978-1461502159
  109. 109. Reimer L. Scanning Electron Microscopy: Physics of Image Formation and Microanalysis. Springer Science & Business Media, 1998. ISBN: 978-3540639763
  110. 110. Krinsley DH, Pye K, Boggs S, Tovey NK Jr. Backscattered Scanning Electron Microscopy and Image Analysis of Sediments and Sedimentary Rocks. Cambridge University Press, 2005. ISBN: 978-0521019743
  111. 111. Compton RG, Harding MS, Pluck MR, Atherton JH, Brennan CM. Journal of Physical Chemistry 1993;97(40) 10416–10420.
  112. 112. Dorozhkin SV. Surface Reactions of Apatite Dissolution. Journal of Colloid and Interface Science 1997;191(2) 489–497.
  113. 113. Xu Wen-G, Fan Hong-R, Hu Fang-F, Santosh M, Yang Kui-F, Lan Ting-G. In situ chemical and Sr–Nd–O isotopic compositions of apatite from the Tongshi intrusive complex in the southern part of the North China Craton: implications for petrogenesis and metallogeny. Journal of Asian Earth Sciences 2015;105 208–222.
  114. 114. Xing Chang-M, Wang ChY, Li C. Trace element compositions of apatite from the middle zone of the Panzhihua layered intrusion, SW China: insights into the differentiation of a P- and Si-rich melt. Lithos 2014;204 188–202.
  115. 115. Roman-Lopez J, Correcher V, Garcia-Guinea J, Prado-Herrero P, Rivera T, Lozano IB. Effect of the chemical impurities on the luminescence emission of natural apatites. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy 2014;126 142–147.
  116. 116. Kamiishi E, Utsunomiya S. Nano-scale reaction processes at the interface between apatite and aqueous lead. Chemical Geology 2013;340 121–130.
  117. 117. Stuart BH. Infrared Spectroscopy: Fundamentals and Applications. Vydavatel John Wiley & Sons, 2004. ISBN: 978-0470011133
  118. 118. Chalmers JM, Edwards HGM, Hargreaves MD. Infrared and Raman Spectroscopy in Forensic Science. John Wiley & Sons, 2012. ISBN: 978-0470749067
  119. 119. McCreery RL. Raman Spectroscopy for Chemical Analysis. Chemical Analysis: A Series of Monographs on Analytical Chemistry and Its Applications—Volume 225. John Wiley & Sons, 2005. ISBN: 978-0471231875
  120. 120. Ferraro JR, Nakamoto K. Introductory Raman Spectroscopy. Academic Press, 2012. ISBN: 978-0323137904
  121. 121. Madou MJ. Solid-State Physics, Fluidics, and Analytical Techniques in Micro- and Nanotechnology. Fundamentals of Microfabrication and Nanotechnology—Volume 1. CRC Press, 2011. ISBN: 978-1420055115
  122. 122. Kumar ChSSR. Raman Spectroscopy for Nanomaterials Characterization. Springer Science & Business Media, 2012. ISBN: 978-3642206207
  123. 123. Fleet ME, Liu X. Local structure of channel ions in carbonate apatite. Biomaterials 2005;26(36) 7548–7554.
  124. 124. Errassifi F, Sarda S, Barroug A, Legrouri A, Sfihi H, Rey C. Infrared, Raman and NMR investigations of risedronate adsorption on nanocrystalline apatites. Journal of Colloid and Interface Science 2014;420 101–111.
  125. 125. Antonakos A, Liarokapis E, Leventouri T. Micro-Raman and FTIR studies of synthetic and natural apatites. Biomaterials 2007;28(19) 3043–3054.
  126. 126. Gómez-Morales J, Iafisco Mi, Delgado-López JM, Sarda S, Drouet C. Progress on the preparation of nanocrystalline apatites and surface characterization: overview of fundamental and applied aspects. Progress in Crystal Growth and Characterization of Materials 2013;59(1) 1–46.
  127. 127. Grunenwald A, Keyser C, Sautereau AM, Crubézy E, Ludes B, Drouet C. Revisiting carbonate quantification in apatite (bio)minerals: a validated FTIR methodology. Journal of Archaeological Science 2014;49 134–141.
  128. 128. Adler HH. Infrared spectra of phosphate minerals: symmetry and substitutional effects in the pyromorphite series. American Mineralogist 1964;49 1002–1015.
  129. 129. Bencnn RM. A relation between internuclear distances and bond force constants. Journal of Chemical Physiscs 1934;2 128–131.
  130. 130. Lever T, Haines P, Rouquerol J, Charsley EL, Van Eckeren P, Burlett DJ. ICTAC nomenclature of thermal analysis, 2014 (IUPAC Recommendations 2014). Pure and Applied Chemistry 2014;86(4) 545–553. DOI 10.1515/pac-2012-0609
  131. 131. Brown ME. Introduction to Thermal Analysis: Techniques and Applications. Hot Topics in Thermal Analysis and Calorimetry—Volume 1. Springer Science & Business Media, 2001. ISBN: 978-1402002113
  132. 132. Brown ME, Gallagher PK. Handbook of Thermal Analysis and Calorimetry: Recent Advances, Techniques and Applications. Handbook of Thermal Analysis and Calorimetry—Volume 5. Elsevier, 2011. ISBN: 978-0080556314
  133. 133. Kerr PF. Optical Mineralogy. 4nd ed., McGraw-Hill, 1977. ISBN: 978-0070342187
  134. 134. Mitra S. Fundamentals Of Optical, Spectroscopic And X-Ray Mineralogy. New Age International, 1996. ISBN: 978-8122409826
  135. 135. Jonasz M, Fournier GR. Chapter 6—Refractive indices and morphologies of aquatic particles. Light Scattering by Particles in Water, 2007; 447–558.
  136. 136. Li Cheng-X, Duan Yong-H, Hu Wen-Ch. Electronic structure, elastic anisotropy, thermal conductivity and optical properties of calcium apatite Ca5(PO4)3X (X = F, Cl or Br). Journal of Alloys and Compounds 2015;619 66–77.
  137. 137. Duan YH, Sun Y. First-principles calculations of optical properties of Mg2Pb. Science China Physics, Mechanics and Astronomy 2014;57(2) 233–238.
  138. 138. Rulis P, Ouyang L, Ching WY. Electronic structure and bonding in calcium apatite crystals: Hydroxyapatite, fluorapatite, chlorapatite, and cromapatite. Physical Review B 70, 155104 (2004). DOI: 10.1103/PhysRevB.70.155104.
  139. 139. Lowell S, Shields JE, Thomas MA, Thommes M. Characterization of Porous Solids and Powders: Surface Area, Pore Size and Density. Particle Technology Series—Volume 16. Springer Science & Business Media, 2012. ISBN: 978-1402023033
  140. 140. Merkus HG. Particle Size Measurements: Fundamentals, Practice, Quality Particle Technology Series—Volume 17. Springer Science & Business Media, 2009. ISBN: 978-1402090165
  141. 141. Allen T. Powder Sampling and Particle Size Determination. Elsevier, 2003. ISBN: 978-0080539324
  142. 142. LeRoy Jackson M, Barak P. Soil Chemical Analysis: Advanced Course. UW-Madison Libraries Parallel Press, 2005. ISBN: 978-1893311473
  143. 143. Taylor HFW. Cement Chemistry. Revised edition, Thomas Telford, 1997. ISBN: 978-0727725929
  144. 144. Ishizaki K, Komarneni S, Nanko M. Porous Materials: Process Technology and Applications. Materials Technology Series—Volume 4. ISBN: 978-1461558118
  145. 145. Borch J, Lyne MB, Mark RE, Habeger Ch. Handbook of Physical Testing of Paper, Volume 2. 2nd ed., CRC Press, 2001. ISBN: 978-0203910498
  146. 146. Donaldson EC, Chilingarian GV, Yen TF. Enhanced Oil Recovery, I: Fundamentals and Analyses. Developments in Petroleum Science. Elsevier, 1985. ISBN: 978-0080868721
  147. 147. Millington RJ, Quirk JP. Permeability of porous solids. Transactions of the Faraday Society 1961;57 1200–1207.
  148. 148. Zhong K, Vasudevan TV, Somasundaran P. Floatability of apatites of different type and origin: role of surface area and porosity. International Journal of Mineral Processing 1993;38(3–4) 177–188.
  149. 149. Cardoso L, Fritton SP, Gailani G, Benalla M, Cowin SC. Advances in assessment of bone porosity, permeability and interstitial fluid flow. Journal of Biomechanics 2013;46(2) 253–265.
  150. 150. Tang H, Guo Y, Jia D, Zhou Y. High bone-like apatite-forming ability of mesoporous titania films. Microporous and Mesoporous Materials 2010;131(1–3) 366–372.
  151. 151. Habibovic P, Juhl MV, Clyens S, Martinetti R, Dolcini L, Theilgaard N, van Blitterswijk CA. Comparison of two carbonated apatite ceramics in vivo. Acta Biomaterialia 2010;6(6) 2219–2226.
  152. 152. Liu Dean-M. Influence of porosity and pore size on the compressive strength of porous hydroxyapatite ceramic. Ceramics International 1997;23(2) 135–139.
  153. 153. Mohammad NF, Othman R, Yeoh FY. Controlling the pore characteristics of mesoporous apatite materials: hydroxyapatite and carbonate apatite. Ceramics International, In Press, Accepted Manuscript, Available online 9 May 2015. doi:10.1016/j.ceramint.2015.04.162
  154. 154. Panteix PJ, Julien I, Abélard P, Bernache-Assollant D. Influence of porosity on the electrical properties of La9.33(SiO4)6O2 oxyapatite. Ceramics International 2008;34(7) 1579–1586.
  155. 155. Pons A, Jouin J, Béchade E, Julien I, Masson O, Geffroy PM, Mayet R, Thomas P, Fukuda K, Kagomiya I. Study of the formation of the apatite-type phases La9.33+x(SiO4)6O2+3x/2 synthesized from a lanthanum oxycarbonate La2O2CO3. Solid State Sciences 2014;38 150–155.
  156. 156. Li B, Liu W, Pan W. Synthesis and electrical properties of apatite-type La10Si6O27. Journal of Power Sources 2010;195(8) 2196–2201.
  157. 157. Medeiros EE, Dias AMS. Experimental and numerical analysis of Vickers hardness testing. IJRRAS 2013;17(1) 9–18.
  158. 158. CiDRA Precision Services, LLC. Mineral Hardness Conversion Chart. http://www.cidraprecisionservices.com/mohs-conversion.html
  159. 159. Saenger AT, Kuhs WF. Structural disorder in hydroxyapatite. Zeitschrift für Kristallographie 1992;199 123–148.
  160. 160. Bhatnagar VM. IR-spectra of fluorapatite and fluorchlorapatite. Experientia 1967;24(1) 10–12
  161. 161. Balan E, Delattre S, Roche D, Segalen L, Morin G, Guillaumet M, Blanchard M, Lazzeri M, Brouder Ch, Salje EKH. Line-broadening effects in the powder infrared spectrum of apatite. Physics and Chemistry of Minerals 2011;38(2) 111–112.
  162. 162. Foil AM, Wilkins ChH. infrared spectra and characteristic frequencies of inorganic ions. Their use in qualitative analysis. Analytical Chemistry 1952;24(8) 1253–1294.
  163. 163. Bonner CE Jr, Chess ChC, Meegoda Ch, Stefanos S, Loutts GB. Volume effects on the Raman frequencies of phosphate ions in fluorapatite crystals. Optical Materials 2004;26 17–22.
  164. 164. Tressaud A, Haufe G. Fluorine and Health: Molecular Imaging, Biomedical Materials and Pharmaceuticals. 1st ed., Elsevier, 2008. ISBN: 978–0444530868
  165. 165. Williams Q, Knittle E. Infrared and raman spectra of Ca5(PO4)3F2-fluorapatite at high pressures: compression-induced changes in phosphate site and Davydov splittings. Journal of Physics and Chemistry of Solids 199657(4) 417–422.
  166. 166. Veiderma M, Knubovets R, Tõnsuaadu K. Structural properties of apatites from Finland studied by FTIR spectroscopy. Bulletin of the Geological Society of Finland 1998;1–2 69–75.
  167. 167. Freund F, Knobel RM. Distribution of fluorine in hydroxyapatite studied by infrared spectroscopy. Journal of the Chemical Society, Dalton Transactions 1977;11 1136–1140.
  168. 168. Knubovets R. Structural mineralogy and properties of natural phosphates. Reviews in Chemical Engineering 1993;9(3–4) 161–216.
  169. 169. Ptáček P, Bartoníčková E, Švec J, Opravil T, Šoukal F, Frajkorová F. The kinetics and mechanism of thermal decomposition of SrCO3 polymorphs. Ceramics International 2015;41(1) 115–126.
  170. 170. Drouet Ch. A comprehensive guide to experimental and predicted thermodynamic properties of phosphate apatite minerals in view of applicative purposes. Journal of Chemical Thermodynamics 2015;81 143–159.
  171. 171. Moiseev G, Leitner J, Sestak J, Zhukovsky V. Empirical dependences of the standard enthalpy of formation for related inorganic compounds enhancing glass formers. Thermochimica Acta 1996;280–281 511–521.
  172. 172. Moiseev GK, Sestak J. Some calculations methods for estimation of thermodynamical and thermochemical properties of inorganic compounds. Progress in Crystal Growth and Characterization of Materials 1995;30(1) 23–81.
  173. 173. Petit A.-T, Dulong P-L. Recherches sur quelques points importants de la Théorie de la Chaleur. Annales de Chimie et de Physique (in French) 1819;10 395–413.
  174. 174. Bureau International des Poids et Mesures the International System of Units (SI). 8th ed., Organisation Intergouvernementale de la Convention du Mètre, STEDI Media, Paris, 2006. ISBN: 92-822-2213-6
  175. 175. Rogers D. Einstein’s Other Theory: The Planck-Bose-Einstein Theory of Heat Capacity. Princeton University Press, 2005. ISBN: 978-0691118260
  176. 176. Latimer WM. Methods of estimating the entropies of solid compounds. Journal of the American Chemical Society 1951;73(4) 1480–1482.
  177. 177. Leitner J, Chuchvalec P, Sedmidubský D, Strejc A, Abrman P. Estimation of heat capacities of solid mixed oxides. Thermochemica Acta 2003;395 47–46.
  178. 178. Kellogg HH, in Fitterer GR (Ed.). Problems in the non-ferrous industries and a note on estimation of heat capacity, in: Applications of Fundamental Thermodynamics to Metallurgical Processes. Gorden and Breach, London, 1967, p.357
  179. 179. Kubaschewski O, Ünal H. An empirical estimation of the heat capacities of inorganic compounds. High Temperatures—High Pressures 1977;9(3) 361–365.
  180. 180. Kubaschewski O, Alock CB, Spencer PJ. Materials thermochemistry. Materials Science Monographs. 6th ed., Pergamon Press, 1993. ISBN: 9780080418896
  181. 181. Spencer PJ. Estimation of thermodynamic data for metallurgical applications. Thermochimica Acta 1988;314(1–2) 1–21.
  182. 182. Moiseev GK, Šesták J. Some calculations methods for estimation of thermodynamical and thermochemical properties of inorganic compounds. Progress in Crystal Growth and Characterization of Materials 1995;30 23–81.
  183. 183. Kumok VN. Problem of Correlation Methods for Evaluation of Thermodynamic and Characteristics, Direct and Reverse Task of Chemical Thermodynamics. Nauka, Novosibirsk, 1987, p.108.
  184. 184. Latimer WM. Methods of estimating the entropies of solid compounds. Journal of the American Chemical Society 1951;73 1480–1482.
  185. 185. Aldabergenov MK, Balakaeva GT, Kokibasova GT. Calculation of thermodynamic functions of aluminates of alkaline and alkaline earth elements. Journal of Physical chemistry 1998;72(5) 808–811.
  186. 186. Gospodinov GG, Mihov DI. Molar heat capacities of the tellurites MnTeO3, MnTe2O5, and Mn2Te3O8. Journal of Chemical Thermodynamics 1993;25(10) 1249–1252.
  187. 187. Dorozhkin SV. A review on the dissolution models of calcium apatites. Progress in Crystal Growth and Characterization of Materials 2001;44(1) 45–61.
  188. 188. Zhang J, Nancolias GH. Unexpected pH dependence of dissolution kinetics of dicalcium phosphate dihydrate. Journal of Physical Chemistry 1994;98(6) 1689–1694.
  189. 189. Zhu Y, Zhang X, Chen Y, Xie Q, Lan J, Qian M, He N. A comparative study on the dissolution and solubility of hydroxylapatite and fluorapatite at 25°C and 45°C. Chemical Geology 2009;268(1–2) 89–96.
  190. 190. Patel PR, Gregory TM, Brown WE. Solubility of CaHPO4·2H2O in the quaternary system Ca(OH)2-H3PO4-NaCl-H2O at 25°C. Journal of Research of the National Bureau of Standards—A. Physics and Chemistry 1974;78(6) 675–681.
  191. 191. Christoffersen J, Christoffersen MR. Kinetics of dissolution of calcium hydroxyapatite V. The acidity constant for the hydrogen phosphate surface complex. Journal of Crystal Growth 1982;57 21–26.
  192. 192. McDowell H, Gregory TM, Brown WE. Solubility of Ca5(PO4)3OH in the System Ca(OH)2-H3PO4-H2O at 5, 15, 25 and 37°C. Journal of Research of the National Bureau of Standards—A. Physics and Chemistry 1977;81(2–3) 273–281.
  193. 193. Tang R, Henneman ZJ, Nancollas GH. Constant composition kinetics study of carbonated apatite dissolution. Journal of Crystal Growth 2003;249(3–4) 614–624.
  194. 194. Shellis RP, Lee AR, Wilson RM. Observations on the apparent solubility of carbonate-apatites. Journal of Colloid and Interface Science 1999;218(2) 351–358.
  195. 195. Chaïrat C, Schott J, Oelkers EH, Lartigue Jean-E, Harouiya N. Kinetics and mechanism of natural fluorapatite dissolution at 25°C and pH from 3 to 12. Geochimica et Cosmochimica Acta 2007;71(24) 5901–5912.
  196. 196. Wang LK, Hung Yung-T, Shammas NK. Physicochemical Treatment Processes. Handbook of Environmental Engineering—Volume 3. Springer Science & Business Media, 2007. ISBN: 978-1592598205
  197. 197. Harouiya N, Chaïrat C, Köhler SJ, Gout R, Oelkers EH. The dissolution kinetics and apparent solubility of natural apatite in closed reactors at temperatures from 5 to 50°C and pH from 1 to 6. Chemical Geology 2007;244(3–4) 554–568.
  198. 198. Chaïrat C, Oelkers EH, Schott J. Lartigue J-E Harouiya N. Fluorapatite surface composition in aqueous solution deduced from potentiometric, electrokinetic, and solubility measurements, and spectroscopic observations. Geochimica Cosmochimica Acta 2007;71(24) 5888–5900.
  199. 199. Smith AN, Posner AM, Quirk JP. Incongruent dissolution and surface complexes of hydroxyapatite. Journal of Colloid and Interface Science 1974;48(3) 442–449.
  200. 200. Deutsch WJ, Siegel R. Groundwater Geochemistry: Fundamentals and Applications to Contamination. CRC Press, 1997. ISBN: 978-0873713085
  201. 201. Thomann JM, Voegel JC, Gramain P. Kinetics of dissolution of calcium hydroxyapatite powder. III: pH and sample conditioning effects. Calcified Tissue International 1990;46(2) 121–129.
  202. 202. Chin KOA, Nancollas GH. Dissolution of fluorapatite. A constant-composition kinetics study. Langmuir 1991;7 2175–2179.
  203. 203. Margolis HC, Moreno EC. Kinetics of hydroxyapatite dissolution in acetic, lactic, and phosphoric acid solutions. Calcified Tissue International 1992;50(2) 137–143.
  204. 204. Paschalis EP, Wikiel K, Nancollas GH. Dual constant composition kinetics characterization of apatitic surfaces. Journal of biomedical materials research 1994;28(12) 1411–1418.
  205. 205. Nelson DGA, Featherstone JDB, Duncan JF, Cutress TW. Effect of carbonate and fluoride on the dissolution behaviour of synthetic apatites. Caries Research 1983;17(3) 200–211.
  206. 206. Christoffersen MR, Christoffersen J. Possible mechanisms for the growth of the biomaterial, calcium hydroxyapatite microcrystals. Journal of Crystal Growth 1992;121(4) 617–630.
  207. 207. Christoffersen J, Christoffersen MR, Johansen T. Kinetics of growth and dissolution of fluorapatite. Journal of Crystal Growth 1996;163(3) 295–303.
  208. 208. Christoffersen MR, Dohrup J, Christoffersen J. Kinetics of growth and dissolution of calcium hydroxyapatite in suspensions with variable calcium to phosphate ratio. Journal of Crystal Growth 1998;186 283–290.
  209. 209. Mafe S, Manzanares JA, Reis H, Thomann JM, Gramain P. Model for the dissolution of calcium hydroxyapatite powder. Journal of Physical Chemistry 1992;96(2) 861–866.
  210. 210. Gasser P, Voegel JC, Gramain P. Surface reactions on hydroxyapatite in the presence of fluoride ions 1. Saturated and congruent conditions. Colloids and Surfaces A: Physicochemical and Engineering Aspects 1993;74(2–3) 275–286.
  211. 211. Schaad P, Poumier F, Voegel JC, Gramain P. Analysis of calcium hydroxyapatite dissolution in non-stoichiometric solutions. Colloids and Surfaces A: Physicochemical and Engineering Aspects 1997;121(2–3) 217–228.
  212. 212. Dorozhkin SV. Acidic dissolution mechanism of natural fluorapatite. II. Nanolevel of investigations. Journal of Crystal Growth 1997;182(1–2) 133–140.
  213. 213. Dorozhkin SV. Inorganic chemistry of the dissolution phenomenon: the dissolution mechanism of calcium apatites at the atomic (ionic) level. Comments on Inorganic Chemistry 1999;20(4–6) 285–299.
  214. 214. Pearce EIF. On the dissolution of hydroxyapatite in acid solutions. Journal of Dental Research 1988;67(7) 1056–1058.
  215. 215. Lowell LC. Dislocation etch pits in apatite. Acta Metallurgica 1958;6(12) 775–778.
  216. 216. Jongebloed WL, Molenaar I, Arends J. Orientation-dependent etchpit penetration and dissolution of fluoroapatite. Caries Research 1973;7(2) 154–165.
  217. 217. Jongebloed WL, van den Berg PJ, Arends J. The dissolution of single crystals of hydroxyapatite in citric and lactic acids. Calcified Tissue Research1974;15 1–9.
  218. 218. Arends J, Jongebloed WL. Dislocations and dissolution in apatites: theoretical considerations. Caries Research 1977;11(3) 186–188.
  219. 219. Arends J, Jongebloed WL. Ultrastructural studies of synthetic apatite crystals Journal of Dental Research 1979;58 (Spec. Issue B) 837–843.
  220. 220. Hignett TP. Developments in Plant and Soil Sciences—Volume 15. Springer Science & Business Media, 2013. ISBN: 978-9401715386
  221. 221. Somani L.L. Dictionary of agricultural and allied sciences—Volume 4. Mittal Publications, 1989. ISBN: 978-8170991410
  222. 222. Nielsson FT. Manual of Fertilizer Processing. Fertilizer Science and Technology—Volume 5. CRC Press, 1986. ISBN: 978-0824775223
  223. 223. Manual Fertilizer, UN Industrial Development Organization, International Fertilizer Development Center. Springer Science & Business Media, 1998. ISBN: 978-0792350323
  224. 224. Naidu S. International Encyclopaedia of Agricultural Science and Technology: Soils and composts. International Encyclopaedia of Agricultural Science and Technology—Volume 4. Mittal Publications, 2008. ISBN: 978-8183242233
  225. 225. Espinosa J. Influence of Rock Phosphate on Available Phosphorus as Measured by Plant Uptake and Soil Extractants. Thesis, INIAP Archivo Historico, 1979
  226. 226. Espinosa J. Abstract Influence of Rock Phosphate on Available Phosphorus as Measured by Plant Uptake and Soil Extractants. Thesis, INIAP Archivo Historico.
  227. 227. Shemilt LW. Chemistry and World Food Supplies: The New Frontiers, Chemrawn II: Invited papers presented at the International Conference on Chemistry and World Food Supplies, Manila, Philippines, 1982. ISBN: 0-08-029243-7.

Notes

  • Scattering is the process where the beam of radiation or particles is deviated from its initial trajectory by the inhomogeneity in the medium which it transverses [7].
  • Other kinds of radiation commonly used for diffraction analysis are neutrons (Section 3.1.2) and electrons (Chapter 3.1.10).
  • Wilhelm Conrad Röentgen (1845–1923) was the rector of the University of Würzburg [8]. The first X-ray photography was published in 1896 [9].
  • The interaction with electrons during X-ray analysis takes place over the distances of 10−10 m (1 Å) [14].
  • Small angle scattering (SAS) of X-rays is abbreviated as SAXS [7].
  • Plasma is defined as an electrically neutral gas which consists of positive ions and free electrons. Plasma have sufficiently high energy to atomize, ionize, and excite virtually all elements in the periodic table, which are intentionally introduced into it for the purpose of elemental chemical analysis [45].
  • Although there are many types of plasma, such as direct current, microwave induced, etc., the ICP is considered the most useful technique for analytical spectroscopy [45].
  • Each element has characteristics first and second ionization potential, which depends on specific electronic structure of given element. Higher ionization potential means that more externally applied energy is required for ionization (thermal radiation, collision with other ion or electron, or exposure to high-energy photons) [45].
  • The method was originally developed in the 1950s and 1960s by Herzog et al and Honig et al to analyze metals and oxides [69], but the basis of SIMS can be traced back to the beginning of the twentieth century with the first experimental evidence of secondary ions given by JJ Thomson in 1910 [68].
  • The definition of terms in thermal analysis was developed by ICTAC (Confederation for Thermal Analysis and Calorimetry).
  • Gaseous environment of the sample, which may be controlled by the instrumentation or generated by the sample [130].
  • Resulting dependence, i.e. any graph of any combination of property vs. time or temperature derived from a thermal analysis technique, should be termed as thermal curve, which is a simplified form of more correct term thermoanalytical curve. The first mathematical derivation of any curve with respect to temperature or time leads to the derivative thermoanalytical curve [130]. Since the name thermogram has medical usage, the thermal analysis curve should not be termed as thermogram [131].
  • Brunauer–Emmett–Teller (BET) theory of multilayered physical adsorption of gas molecules on a solid surface [143].
  • For example, the Blaine method, where fixed volume of air passes through the bed at steadily decreasing rate, which is controlled and measured by the movement of oil in a manometer, the time required being measured. The method is widely used for the determination of specific surface area of cements [143].
  • The flow of fluids through porous materials is of great importance in the fields of industrial chemistry, oil technology, and agriculture. In general, it may be stated that the principal interest is in the transport through reactive materials [147].
  • Experimental techniques for the measurement of relative permeability can be divided to steady- (1) and unsteady-state (displacement) methods (2). The steady-state method was developed by Hassler (1944). Semipermeable membranes are provided at each end which keeps the fluids separated, except inside the core where the fluids flow simultaneously. The pressure is measured in each phase through semipermeable barriers, and the pressure differences between the phases are maintained constant throughout the medium so as to eliminate the capillary end effect as well as to ensure a uniform saturation along the core. The saturation can be altered by applying capillary pressures across the nonwetting phase ports and wetting phase semipermeable membranes [146].
  • The Mohs scale of mineral hardness is graduated as follows: talc (1), gypsum (2), calcite (3), fluorite (4), apatite (5), orthoclase (6), quartz (7), topaz (8), corundum (9), and diamond (10). Apatite should be also scratched by steel knife (up to 5.5) and glass (up to 6).
  • Hardness tester LECO AMH 43. The method is also known as the Vickers pyramid number (HV) or the diamond pyramid hardness (DPH).
  • Diamond pyramid with apical angle of 136°.
  • Since it is present as “free carbonate,” the sample cannot be considered as carbonate-fluorapatite (Section 2.6 and 4.6.1).
  • As the results of thermal analysis revealed (Section 3.2.4), the Rietveld analysis slightly overestimates the content of free carbonate in the sample.
  • Baseline corrected spectrum measured by KBr technique.
  • Abbreviation used for the expression of intensity and width of peak in the spectrum: very weak (vw), weak (w), middle (m), strong (s) and very strong (vs), shoulder (sh), broad (b), very broad (vb), and sharp (sp). Spectral bands related to impurities are abbreviated as imp [162].
  • It was calculated as 100 × molar mass (molar weight) of CO2/molar mass of CaCO3.
  • Theoretical mass loss (3.77%) of fluorapatite was used to set the fixed value of parameter A2.
  • One calorie is 4.184 J (joules). Gram-atom [gm] (and gram-molecule) was used to specify the amount of chemical elements or compound. These units had a direct relation with “atomic weights” and “molecular weights,” which are in fact relative masses. “Atomic weights” were originally referred to the atomic weight of oxygen, by general agreement taken as 16. Although physicists separated the isotopes in a mass spectrometer and attributed the value of 16 to one of the isotopes of oxygen, chemists attributed the same value to the (slightly variable) mixture of isotopes 16, 17, and 18, which was for them naturally occurring element oxygen. Finally, an agreement between the International Union of Pure and Applied Physics (IUPAP) and the International Union of Pure and Applied Chemistry (IUPAC) brought this duality to an end in 1959/1960. Physicists and chemists have ever since agreed to assign the value 12, exactly, to the so-called atomic weight of the isotope of carbon with the mass number 12 (carbon 12, 12C), correctly called the relative atomic mass Ar(12C). The unified scale thus obtained gives the relative atomic and molecular masses, also known as the atomic and molecular weights, respectively [174]. The law is also known as Dulong and Petit principle, which can be expressed in modern unit as: atomic weigh × specific heat ≈ ∂(3kTNA)/∂T 3kNA ≈ cV ≈ 25 J·K−1mol−1, i.e. the atomic weight of solid element multiplied by its molar specific heat is a constant [175].
  • The full equation for the temperature dependence is Cpm° = A + BT + C/T2 + DT2 + F/T1/2 [J·K−1·mol−1].
  • Table 4 [172] refers to the values of entropy contribution from work [176] where some data for PO43− ion are given in brackets. In order to verify this value (for the charge of cation, i.e. Ca2+, the contribution of PO43− anion is 17 [calories] × 4.184 = 71.13 J·K−1·mol−1, it is possible to calculate it from recommended thermodynamic data from Fig. 25. If (for example) three apatite end members were used, it is possible to calculate the contribution to PO43− anion in hydroxylapatite, fluorapatite, chlorapatite, and bromapatite as follows: H A P : 10 ⋅ 39.1 + 6 x + 2 ⋅ 18.83 = 780 ⇒ x = 58.56 J ⋅ K − 1 ⋅ mol − 1 F A P : 10 ⋅ 39.1 + 6 x + 2 ⋅ 17.00 = 728 ⇒ x = 50.50 J ⋅ K − 1 ⋅ mol − 1 ClAP : 10 ⋅ 39.1 + 6 x + 2 ⋅ 31.80 = 835 ⇒ x = 63.40 J ⋅ K − 1 ⋅ mol − 1 BrAP : 10 ⋅ 39.1 + 6 x + 2 ⋅ 45.70 = 870 ⇒ x = 64.60 J ⋅ K − 1 ⋅ mol − 1
  • Reaction Eq. 7 and Eq. 12 are widely used for the description of dissolution process of apatite [112] using stoichiometry pertinent single or double apatite formula and Z = F or OH.
  • Out of equilibrium state, IAP is not equal to K (see the discussion to Eq. 20). Double formula of apatite is assumed then with respect to the apatite stoichiometry; the law for ionic activity product has the following form:
  • ∆rG° = -RT ln K (in the equilibrium state) and ln K = ln 10 log K. From that, it should be derived that ∆rG° = -8.314·298.15 ln 10 log K = -5708 log K [J·mol−1], where ln 10 ≈ 2.303. Since log10K = log K = lg K, it can also be written as ∆rG° = -5708 lg K.
  • ∆rG° = ∑νi ∆G°f,i, where νi denotes the stoichiometric coefficient of given species and ∆G°f,i its standard enthalpy of formation.
  • Since the saturation of solution with respect to Ca5(PO4)3Z means that K′ = 0 (equilibrium state), it can be derived that a(Ca2+)5 a(PO43−)3 a(Z−) – K´a(H+)3 = 0 and then ∆rG° = -RT ln (a(Ca2+)5 a(PO43−)3 a(Z−) / K´a(H+)3). Since A = -∆rG° (please see note 35), A = -RT ln (K´a(H+)3 / a(Ca2+)5 a(PO43−)3 a(Z−)).
  • The relationship between the reaction Gibbs energy and chemical affinity: A = -∆rG° was introduced by T. de Donder.
  • Since the system is not in the equilibrium state IAP ≠ K.
  • The ranges of Si near zero are generally considered to be within the equilibrium zone for the mineral. The ranges of SI = 0 ± 0.5 and 0 ± (5%) (lg Kmineral) were used in various studies [200].
  • Neutral ammonium citrate is prepared by dissolving required amount of citric acid and neutralizing it with ammonium hydroxide. The pH of the reagent is adjusted to neutral [221].
  • The amount of used solution is also expressed in the name of the method, e.g. 100 ml method or 150 ml method [222].
  • The neutral ammonium citrate test was used as the official method in the United States, and the test by acidic acid was developed for the comparison [222].
  • It is believed that weak citric acid solution imitates the condition near the plant roots [222].

Written By

Petr Ptáček

Reviewed: January 8th, 2016 Published: April 13th, 2016