# Time-domain THz spectroscopy of the characteristics of hydroxyapatite provides a signature of heating in bone tissue

Marie Plazanet, Jordanka Tasseva, Paolo Bartolini, Andrea Taschin, Renato

Torre, Christèle Combes, Christian Rey, Alessandro Di Michele, Mariana

Verezhak, Aurélien Gourrier

► To cite this version:

Marie Plazanet, Jordanka Tasseva, Paolo Bartolini, Andrea Taschin, Renato Torre, et al.. Time-domain THz spectroscopy of the characteristics of hydroxyapatite provides a signature of heating in bone tissue. PLoS ONE, 2018, 13 (8), pp.1-16. 10.1371/journal.pone.0201745. hal-01899210

HAL Id: hal-01899210

https://hal.science/hal-01899210

Submitted on 11 Jul 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
**To cite this version:**

Plazanet, Marie and Tasseva, Jordanka and Bartolini, Paolo and Taschin, Andrea and Torre, Renato and Combes, Christèle and Rey, Christian and Di Michele, Alessandro and Verezhak, Mariana and Gourrier, Aurélien *Time-domain THz spectroscopy of the characteristics of hydroxyapatite provides a signature of heating in bone tissue.* (2018) PLoS ONE, 13 (8). 1-16. ISSN 1932-6203
RESEARCH ARTICLE

# Time-domain THz spectroscopy of the characteristics of hydroxyapatite provides a signature of heating in bone tissue

Marie Plazanet¹,²*, Jordanka Tasseva³ª, Paolo Bartolini³, Andrea Taschin³, Renato Torre³, Christèle Combes⁴, Christian Rey⁴, Alessandro Di Michele², Mariana Verezhak¹ᵇ, Aurelien Gourrier¹*

¹ Univ. Grenoble Alpes, CNRS, LIPhy, Grenoble, France, ² Dipartimento di Fisica e Geologia, Università degli Studi di Perugia, Perugia, Italy, ³ European Laboratory for Non-Linear Spectroscopy (LENS) and Dip. di Fisica ed Astronomia, Università di Firenze, Sesto Fiorentino, Italy, ⁴ CIRIMAT, Université de Toulouse, CNRS, INPT-ENSIACET, Toulouse, France

¤a Current address: Istituto Nazionale di Fisica Nucleare, Napoli, and Istituto Nazionale di Ottica- CNR, Pozzuoli, Italy

¤b Current address: Paul Scherrer Institut, Villigen, Switzerland

* marie.plazanet@univ-grenoble-alpes.fr (MP); aurelien.gourrier@univ-grenoble-alpes.fr (AG)

<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Plazanet M, Tasseva J, Bartolini P, Taschin A, Torre R, Combes C, et al. (2018) Time-domain THz spectroscopy of the characteristics of hydroxyapatite provides a signature of heating in bone tissue. PLoS ONE 13(8): e0201745. https://doi.org/10.1371/journal.pone.0201745

**Editor:** Yogendra Kumar Mishra, Institute of Materials Science, GERMANY

**Received:** January 9, 2018

**Accepted:** July 20, 2018

**Published:** August 23, 2018

**Copyright:** © 2018 Plazanet et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** We thank the ILL for the use of computational cluster and VASP 5.4 licence. This work was partially funded by Regione Toscana, prog. POR-CROFSE-UNIFI-26 and by Ente Cassa di Risparmio Firenze, prog. 2016-0866, and by the ECs Seventh Framework Programme, Grant Agreement No. 284464 (RT). TD-THz measurements were supported by the program

## Abstract

Because of the importance of bone in the biomedical, forensic and archaeological contexts, new investigation techniques are constantly required to better characterize bone ultrastructure. In the present paper, we provide an extended investigation of the vibrational features of bone tissue in the 0.1-3 THz frequency range by time-domain THz spectroscopy. Their assignment is supported by a combination of X-ray diffraction and DFT-normal modes calculations. We investigate the effect of heating on bone tissue and synthetic calcium-phosphates compounds with close structure and composition to bone mineral, including stoichiometric and non-stoichiometric hydroxyapatite (HA), tricalcium phosphate, calcium pyrophosphate and tetracalcium phosphate. We thus demonstrate that the narrow vibrational mode at 2.1 THz in bone samples exposed to thermal treatment above 750 °C arises from a lattice mode of stoichiometric HA. This feature is also observed in the other synthetic compounds, although weaker or broader, but is completely smeared out in the non-stoichiometric HA, close to natural bone mineral composition, or in synthetic poorly crystalline HA powder. The THz spectral range therefore provides a clear signature of the crystalline state of the investigated bone tissue and could, therefore be used to monitor or identify structural transitions occurring in bone upon heating.

## Introduction

Bone grafts are essential in dental and bone clinical practice to bridge pathological fractures, fill large defects or strengthen a diseased tissue. However, due to numerous problems associated with donor site morbidity and limited availability (autologous grafts) or with immunological response and infectious risk (allografts, xenografts), synthetic bone substitutes have
Laserlab (i.e. Laserlab-Europe, H2020 EC-GA (654148)). M. Verezhak acknowledges support from the Ph.D. program of the Nanosciences Foundation of Grenoble (France). There was no additional external funding received for this study.

**Competing interests:** The authors have declared that no competing interests exist.

become recognized as a valuable alternative [1]. A wide variety of calcium-phosphate based ceramic scaffolds have thus been developed and combined with growth factors and stem cells to improve the osteoconductive and osteoinductive properties [2]. Those mostly consist of hydroxyapatite (HA) and/or β-tricalcium phosphate (β-TCP) manufactured in the form of pastes, granules, coatings or bulk materials [3]. Animal bone, considered as a waste product from the food industry, has extensively been tested as a natural resource for HA production by high temperature heating ( > 1100 °C). Bovine bone, for example, is commercially used to fabricate PepGen P-15®, Endobon® and Cerabone® biomaterials [3] and there is a high potential for the use of other types of bones (fish, ovine, porcine...). Independently of the nature of the starting materials, the precise chemistry and morphometry of those end-products have been shown to influence the tissue formation through the activity of bone cells [4, 5] and therefore require a precise control during (or at least following) the manufacturing process.

Since atomic or molecular scale sensitive techniques are required to monitor the structural evolution of bone mineral or synthetic HA, vibrational spectroscopy is widely used to characterize the effect of heating on bone and calcium phosphate minerals, often in conjunction with X-ray diffraction (XRD), scanning electron microscopy (SEM) and thermogravimetric analysis [6–12]. Interestingly, similar studies were also reported in the fields of forensic science and archeology [13–16]. Although undertaken in very different contexts, those studies provide very complementary mechanistic information. The infra-red (IR) and Raman spectroscopic investigations were based on well established assignment of the vibrational bands from 400 up to 3600 cm⁻¹ [17–19], or even down to ~100 cm⁻¹ [20]. Both collagen and mineral phases present characteristic peaks in this spectral range.

Those studies provided a fine characterization of bone transformations upon heating. It was shown indeed that, up to 300 °C, the mineral crystalline phase is only slightly affected by the loss of hydrogen phosphate and carbonates ions, as well as by a moderate crystal growth [21–25]. In this temperature range, the most important changes concern the collagen degradation from denaturation to calcination above 300 °C, leading to a brown and black coloration of the samples. Beyond 600 °C, the color fades and turns white around 800 °C indicating a transformation into a pure mineral phase. The loss of organic matter is accompanied by important crystallographic changes in the mineral phase starting around 750 °C, which were accurately characterized by X-rays and TEM investigations [24, 26, 27]. Beyond this characteristic temperature, the crystallite average dimension then increases from tens to thousands of nm in conjunction with sintering processes [12]. This structural transition is critical for the mechanical, microstructural and surface properties of the material [28, 29] and its precise identification would therefore provide a very important cue of the transformation process induced by heating.

One drawback of these studies is that they are generally performed on complex instruments or involve specific sample preparation, which is relatively impractical for fast and cost-effective screening of a large number of samples. Hence, a method that could provide a vibrational signature of this transition would prove highly beneficial. Up to now, no striking features were revealed in the low frequency (far-IR) range which can be advantageously investigated by THz time-domain spectroscopy (THz-TDS). This technique was mostly employed for imaging purposes [[30, 31]]. So far, THz-TDS measurements mainly showed a continuous increase of the absorption coefficient as a function of frequency [32–34] and did not reveal any underdamped (absorption) feature in bones or HA. However, the characterization of phonon lattice modes which are expected in the THz frequency range, could provide valuable information since they are particularly sensitive to crystalline environment.

Because THz-TDS can easily be implemented at relatively low cost using commercial instruments, we undertook an investigation of heated cortical bone of a bovine femur heated
up to 1000 °C, covering the frequency range from 0.1 to 3 THz. An additional set of synthetic calcium phosphate powders, of direct or related interests for bone-like materials were also studied: stoichiometric HA, α and β-tricalcium phosphate, β-calcium pyrophosphate and tetracalcium phosphate. To complete the work, we also characterized a non-stoichiometric commercial HA as received and heated above 800°C. We observed the appearance of a well defined peak at ~ 2.1 THz in the absorption spectra of bone above 750 °C, simultaneous to the well characterized recrystallisation of the mineral part. We then combined the spectroscopic measurements with X-ray diffraction, enabling structural characterization, and density functional theory (DFT) simulations of the crystal phonon modes of the bone mineral components that are responsible for the absorption spectra. As demonstrated in the following sections, we could identify the signature of the crystalline state of bone tissue from the THz phonon spectrum.

# Materials and methods

## Bone samples

The samples were prepared from a bovine femur obtained from the local slaughterhouse, 38120 Fontanil-Cornillon, France. The periosteum and marrow were mechanically removed and a transverse cross-section of ~ 15 mm in thickness was sawed in the diaphysis, fixed with ethanol 70% v/v, dehydrated in a graded series of ethanol solutions of 80% v/v, 90% v/v and 100% v/v. A cortical block of 10x10x10 mm³ was selected in the medial quadrant. Light microscopy examination revealed a homogeneous fibro-lamellar structure with alternating layers of fibrous and lamellar bone with characteristic dimensions of 100 µm. Longitudinal sections of approximately 100 to 150 µm in thickness were cut from this block using a high precision low-speed diamond saw (Mecatome T210, PRESI) (Fig 1a).

One section was used as a control, while the others were heated to 100, 200, 300, 400, 500, 600, 650 700, 750, 800, 850 and 1000 °C in an unvented oven (Carbolite Gero, ELF) and cooled in air. The samples were treated sequentially: each sample was inserted at the final heating temperature of the previous one, the next temperature was reached, followed by curing for 10 minutes. The temperature precision of the thermocouple was ~ 2-3 °C and the sample reached the given temperature in less than a minute.

<figure><img src="image_4.png" /><figcaption>**Fig 1. Macroscopic and nanoscale evolution of HA upon heating.** a) sample preparation scheme; b) color changes of bone sticks as a function of temperature; scale bar 10 mm; c) SEM images showing the nanoscale heating effect on ns-HA, and compared to s-HA and β-TCP. Scale bar is 200 nm for all the images.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0201745.g001
## Synthetic and commercial compounds

We investigated a series of calcium phosphate model compounds synthesized in CIRIMAT (Toulouse, France). These include: the α and β phases of tricalcium phosphate (α- and β-TCP, Ca₃(PO₄)₂), β phase of calcium pyrophosphate (β-CPP, Ca₂P₂O₇), tetracalcium phosphate (TTCP, Ca₄(PO₄)₂O) and stoichiometric hydroxyapatite (s-HA, Ca₅(PO₄)₃OH).

Briefly, α-TCP was prepared by heating β-TCP at 1350°C for 1 hour in a platinum crucible and quenching in liquid nitrogen. β-TCP was synthesized according to the method described by Heughebaert and Montel [35]. β-CPP was obtained by heating dicalcium phosphate dihydrate at 900°C for 3 hours [36]. TTCP was prepared by heating a thoroughly ground mixture of equimolar quantities of calcium carbonate and dicalcium phosphate dihydrate at 1400°C, in a platinum crucible, for 6 hours, and quenching in liquid nitrogen. The stoichiometric HA was synthesized by double decomposition of aqueous solutions of calcium nitrate and diammonium phosphate at boiling temperature; after filtration the precipitate was heated at 1000°C for 15 hours [37]. All samples were controlled by X-ray diffraction and FTIR spectroscopy. They were not found to contain any detectable foreign crystalline phases.

Eventually, we studied a powder of non-stoichiometric hydroxyapatite (referred as ns-HA) ref. 55497 from Sigma-Aldrich. This last powder was also heat treated in the same conditions applied to the bone slices, at temperatures of 650, 750, 800 and 1000°C.

## THz Time Domain Spectroscopy setup

For measuring the frequency dependence of the refractive index, *n*(*ω*), and the absorption coefficient, *α*(*ω*), of bone samples, we used a THz-TDS set-up in transmission configuration, which enabled us to study the optical parameters in the range of frequency 0.1-4 THz. The THz pulses are produced by a photoconductive antenna excited by a femtosecond optical laser pulse (λ = 780 nm, Δt = 120 fs at 100 MHz) and biased with a sinusoidal voltage for a lock-in detection. The emitted THz radiation is collimated by an off-axis parabolic mirror then focused on the sample by a second parabolic mirror. By other two off-axis parabolic mirrors, the THz beam is again collimated and finally focused on a second photoconductive antenna for its detection. This second antenna works as current generator: free electron-couple carriers are produced into the antenna gap by a second optical laser pulse whilst the THz field now acts as a bias. Thus, the temporal evolution of the photocurrent amplitude, acquired by the changing time delay between excitation and gate pulses, is directly related to the electric field amplitude of the THz radiation. This current signal is amplified by a lock-in amplifier and digitized by an acquisition board. By a homemade software we simultaneously acquire the processed signal and delay line encoder, and retrace the final time-dependent THz field. The whole THz set-up is enclosed in a nitrogen purged chamber for removing the water vapor contribution present at the THz frequencies spanned by the experiment. Exhaustive description about our THz-TDS set-up can be found in previously published papers [38–40].

The ratio between the THz field transmitted after the sample, *E*<sub>t</sub>(*ω*), and the incident field, *E*<sub>i</sub>(*ω*), is named transfer function of the material, *H*(*ω*) [38]. The optical properties of the material, *n*(*ω*) and *α*(*ω*), can be fully extracted by the experimental transfer function, *H*<sub>exp</sub>(*ω*), which can be obtained by the ratio of the complex Fourier transform of the sample and the reference signal. However, the analytical expression for *H*(*ω*) is not given in a closed-form and the optical parameters, *n*(*ω*) and *α*(*ω*), can not be calculated without an iterative algorithm. The complexity of this process strongly depends on the nature of the THz signal and the sample thickness. Recently, we have implemented an innovative experimental procedure and numerical method to analyze the transmission THz-TDS signal of samples composed of multiple thin layers [39, 40]. The numerical method is quite complex and is based on a polynomial
fit of the optical parameters. By this procedure, we can safely extract the real frequency depen-
dence of the transmission parameters, totally removing the multiple reflections due to the
Fabry-Perot effect. This data analysis enables the measurement of the n(ω) and α(ω) in an
absolute scale, obtaining a complete quantitative set of the THz transmission parameters for
bone samples.

X-ray diffraction

The X-ray diffraction data were recorded on a laboratory diffractometer at the Physics and Geology department of the University of Perugia, using the radiation of the Cu K<sub>α</sub> (1.5418 Å). A monochromatic incoming beam was obtained using a pyrolytic graphite monochromator which allows removing the K<sub>β</sub> line. Higher order wavelength contamination was avoided using a scintillation detector having adequate energy resolution. The bone sections were mounted on a holed aluminium plate.

Scanning electron microscopy

The SEM characterization was performed on FE-SEM LEO1525 Zeiss. Powder samples were glued on a graphite layer and metalized by depositing a Chromium layer of 8 nm thickness.

DFT calculations

The phonon density of states was calculated using density functional theory (DFT) imple-
mented in the VASP 5.4 software package. For each compound, the crystalline structure from
experimental diffraction data (HA [41], β-TCP [42], β-CPP [43], TTCP [44]) was used as the
initial configuration. Atomic positions were optimized without any change of cell volume or
shape in order to find the minimum energy structure using PAW-PBE ultrasoft pseudopoten-
tials and an energy cut-off of 700 eV for the plane wave basis. The optimized structures are in
good agreement with the experimental ones, with averaged discrepancies in atomic distances
smaller than 1 %. The normal modes calculation was then performed using the linear response
theory (IBRION = 7), and the IR absorption intensities were calculated according to the meth-
ods proposed by Karhanek et al. [45]. In all the calculations, the three translational modes
appear at negative frequencies, although not lower than -0.1 THz. All the other frequencies
were positive, except for β-TCP where two additional modes have frequencies between 0 and
-1.0 THz. The calculated frequencies are in good agreement with the calculations of reference
[46]. The spectral characteristics were convolved with a FWHM 0.1 THz Gaussian to account
for instrumental resolution function. As a second test, we computed the full IR spectrum of
HA and β-TCP, and got results in good agreement with experimental data in the 400–1400
cm⁻¹ frequency range, as published in references [14] and [47] respectively.

Results

SEM

As previously reported [21], no significant color change was observed until 150-200 °C where the bone blocks progressively became brown (200 °C), black (300-600 °C) before turning gray (600-700 °C) and white (> 800 °C). Note the gray speckles at 800 °C which tend to indicate a heterogenous transformation (See Fig 1b). At the nanoscale, the ns-HA which is close to bone mineral structure and composition is formed of ~ 20 nm diameter nanocrystals which fuse and grow up to hundreds of nanometers above 750 °C (See Fig 1c).
THz-TDS

The absorption spectra measured by THz-TDS on control and samples thermally treated within a temperatures range of 100–800 °C are presented in Fig 2a. The overall decrease of the absorption coefficient with temperature can be assigned to the progressive loss of collagen content in the bone samples. The spectra measured on the samples treated at temperatures

<figure><img src="image_5.png" /><figcaption>**Fig 2.** THz absorption spectra of bone and synthetic reference calcium phosphate compounds. a) absorption coefficients of an untreated bone sample and a series of temperature-treated samples up to 800 °C; b) absorption coefficients from a series of high-temperature (from 800 to 900 °C) treated bone samples; c) absorption coefficients of commercial ns-HA powder samples, untreated and treated at different temperatures; d) absorption coefficients from a series of reference calcium phosphate compounds: stoichiometric hydroxyapatite s-HA, TTCP, α and β-TCP, and β-CPP.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0201745.g002
≥800°C show a new absorption peak at about 2.1 THz (see Fig 2b), which intensity slightly increases up to 900 °C. After storing the sample for two weeks at room temperature, the spectrum slightly evolved, suggesting a slow relaxation of the bone structure in the samples (sample labeled as “aged” in Fig 2b).

In Fig 2c, we report the THz spectra of commercial HA powders, as received and treated at two different temperatures. Similarly to the bone samples, the ns-HA exhibits the 2.1 THz peak after heating at 800 °C and higher. We can therefore rule out the hypothesis that this peak is hidden in the samples heated at temperature lower than 750°C by the absorption background produced by the presence of organic components. The peak growth can therefore only be assigned to an increase of the size of the crystallites during the recrystallization process at ~750°C and/or to a simultaneous ordering of the crystalline structure.

Since it is known that phase transitions and decomposition can take place in hydroxyapatite at different temperatures [48], we also investigated several compounds that may be present in the heated bone: α- and β-TCP, β-CPP, TTCP and s-HA. The corresponding THz absorption spectra (reported in Fig 2d) clearly show that the 2.1 THz vibrational band is particularly intense and narrow in s-HA, and is also well defined in β-CPP. The spectra of α-, β-TCP and TTCP are not featureless but only exhibit broad bands.

## X-ray investigation of heated HA

X-ray diffraction was used to assess the existence of phase transitions previously reported in heat-treated samples. X-ray diffractograms, measured on bone samples and commercial hydroxyapatite (ns-HA) are presented in Fig 3a and 3b, respectively. Additionally, the simulated patterns from refined structures based on experimental data available from the literature (see Materials and methods section) are shown in Fig 3c. The comparison between experimental and calculated patterns has to be done keeping in mind the finite resolution of the diffractometer of about 0.5 °. According to the Scherrer equation that links the width of the diffraction peaks and the extension of the monocrystalline domains, our resolution corresponds to a characteristic size of a few nanometers, indicating that we are in any case limited by the instrument resolution and cannot evaluate the size of our crystallites.

Structural transitions clearly occur in bone tissue and ns-HA around 750 °C, although they do not result in the same crystalline phase. The ns-HA heated at 1000 °C exhibits Bragg peaks at 27.9, 31.0 and 34.5° that could be assigned to the presence of a significant quantity of β-TCP, overlapped with a HA ‘background’. Moreover, the broad feature around 38° can be assigned to the formation of CaO [49]. This observation is consistent with numerous studies on the thermal stability of non-stoichiometric HA, that are known to undergo various transformations upon heating [48, 50]. Indeed, the FTIR spectrum of the ns-HA before heating exhibits two bands at 870 and 1140 cm⁻¹ that could be assigned to the presence of HPO₄²⁻ ions, in agreement with the calcium deficiency of the apatite. Upon heating, the ns-HA decomposes in various compounds depending on its Ca/P ratio, in this case into TCP and HA. The FTIR spectrum of the heated ns-HA sample also reflects the TCP and HA decomposition (see S2 Fig). While in stoichiometric HA (s-HA), we only expect a phase transition from a monoclinic (space group P2₁/b) to a hexagonal (P6₃/m) lattice but no decomposition.

The phase transformation in heated bone tissue appears to be different from the one occurring in ns-HA. The X-rays diffractograms of the bone samples heated above 750 °C do not indicate the formation of β-TCP nor of CaO. The observed differences in diffraction pattern with the s-HA could simply be caused by the remaining substituting ions and, even after crystall growth, a less crystalline phase than what would be measured on perfect s-HA.
<figure><img src="image_6.png" /><figcaption>**Fig 3.** X-ray diffraction patterns collected from a) untreated and heat-treated bone tissues, and b) untreated and heat-treated commercial hydroxyapatite powder (ns-HA). c) Simulated patterns for synthetic calcium phosphate compounds based on refined structures [41–44].</figcaption></figure>

https://doi.org/10.1371/journal.pone.0201745.g003

DFT calculation of normal modes of vibration

In Fig 4 we compare the DFT-calculated absorption spectra of some of the investigated com-
pounds (solid lines) with the THz experimental data (dashed lines). The experimental data
have been rescaled in frequency to match the peak positions of the simulated spectra (factor
1.33, so that 2.1 THz measured aligns with 2.8 calculated). It is important to note that the num-
ber of atoms per unit cell, N, is very different for each compound: N = 44, 60, 88 and 276 for
<figure><img src="image_7.png" /><figcaption>**Fig 4.** Simulation of THz spectra (continuous lines) by calculation of the normal modes of vibrational with DFT and experimental THz absorption data (dotted lines). In order to facilitate the visual comparison of these spectra, we adapted the experimental amplitude and frequency scales to match the 2.1 THz measured peak with the 2.8 THz calculated peak.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0201745.g004

HA, TTCP, β-CPP and β-TCP, respectively. 3N defines the total number of vibrational modes present in each compound, that are in all compounds approximately spread over the frequency range of 0-30 THz (without taking into account the isolated O-H stretching mode). We therefore observe, as expected, that the number of modes varies greatly between the considered compounds in the investigated frequency range.

We only find a qualitative agreement between the calculated and measured spectra. However, some interesting observations can be made. All compounds present vibrational modes in the 0.1-3.0 THz range, arising from lattice distortions. Among the compounds, the HA presents the most intense optical phonon at 2.8 THz, that we readily assign to the one observed in the experimental spectrum. The presence of an isolated single mode comes from the HA lattice structure: it has a small asymmetric unit cell enabling a low inhomogeneous broadening. This is also the case for the β-CPP compound which has a limited number of modes in this spectral region. The other compounds, TTCP and β-TCP, have more vibrational modes as can be
expected from the number of atoms per unit cell. Moreover, in real samples, the modes are broadened by crystalline defects, so the overlapping of these broad peaks should lead to the observed continuum.

Discussion

Our THz investigation of bone tissue clearly shows the appearance and increase of a well
defined absorption peak at 2.1 THz induced by the thermal treatment, which is a characteristic
lattice mode of HA crystals. This spectral change arises from the modifications in the mineral
phase of bone tissue, induced in our study by the heating process. Indeed, in the low tempera-
ture range (≤ 750 °C) the organic components are gradually destroyed and the mineral phase
present in bone tissue is composed mainly of nanometric crystals that do not contribute to any
specific signature in the THz spectrum. The THz absorption spectrum therefore follows a
monotonic slope that continuously lowers. In the high temperature range (≥ 800 °C) instead,
only the mineral components are left in bone tissue and the THz measurements exhibit a well
defined phonon vibration of 2.1 THz frequency that can be assigned to HA crystals. The cen-
tral frequency and the width of the absorption peak prove that the HA crystals must be charac-
terized by a well defined lattice structure (i.e. extended lattice, not affected by numerous
impurities/defects). The X-ray diffraction investigations confirm that the thermally treated
bone has a lattice structure similar to HA crystals. This scenario is in full agreement with other
experimental results [24, 51].

A deeper analysis was carried out by fitting the spectra with a polynomial background and a Voigt function to extract the characteristics of the peak, i.e. frequency, shape and width. Their evolution as a function of temperature is plotted in Fig 5. From the observation of the spectra in Fig 5a, we note that the absorption coefficient in the 1.5-2.0 THz region, before the peak, is identical (within spectrum reproducibility) for bones samples heated at temperatures greater than 800°C and s-HA. This is in agreement with the view that bone HA becomes more stoichiometric after the heating transition at ~750 °C. The non-stoichiometric heated powders indeed show a higher absorption coefficient, witnessing of a different composition, as for example mixed HA and TCP phases. Then, the fitting parameters of the peak give quantitative informations. The clearest behavior is shown by the peak area (Fig 5c), that increases with the heating temperature, up to the highest value in s-HA. The frequency, also plotted in Fig 5c, is another indicator of the crystalline state, since it shifts toward lowest values upon increasing the heating temperature. It even surprisingly reaches a value lower than s-HA, which may however be assigned to a difference in profile resulting from surface effects in the bone mineral nanocrystals with respect to s-HA powder. Eventually, the heated ns-HA powders exhibit a higher frequency and broader peak at ~ 2.2 THz, confirming the presence of several phases in the samples identified in the diffraction patterns. The widths of the lorentzian and gaussian contributions to the Voigt profile also follow interesting trends, as represented in Fig 5b. The overall (Voigt) width is found to decrease upon increasing temperature, which is associated with crystal growth. Besides, the width of the gaussian contribution decreases with the heating temperature, until a purely lorentzian profile for the s-HA sample. Such a profile evolution clearly indicates a heterogeneous (gaussian) broadening in the bone spectra with respect to synthetic crystals. This indicates that, even at high temperature, a non-negligible degree of disorder remains in bone mineral when compared to the idealized hydroxyapatite structure. Only the 900°C-aged sample shows significant deviations from those trends. At this stage, it is important to report on the aging of the sample and underline that further characterizations are needed to conclude on the nature of the process.
<figure><img src="image_8.png" /><figcaption>**Fig 5.** (a) Magnification of THz-TDS spectra of bones and HA samples around the peak at 2.1 THz (dots) and fit with a Voigt peak profile on top of the background (lines). For ns-HA, only the experimental data are shown as lines; (b) evolution of the gaussian and lorentzian contributions to the Voigt profile for the heated bone samples (left), as well as aged sample and s-HA (right, as indicated along the x-axis.); (c) evolution of the peak frequency (left axis) and Voigt width (right axis).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0201745.g005

All those results are in perfect agreement with X-ray diffraction investigations widely
spread in the literature, therefore proving THz-TDS to be a quantitative technique for phase
and crystallinity investigations in bone tissues.

The investigation of commercial HA and of a series of synthetic calcium phosphate com-
pounds related to apatite, with close chemical composition and structure, enables a deeper
understanding and shows the broad utility of the technique. DFT normal mode calculation
proves that all the investigated apatite-related compounds have normal mode in the THz
range. Furthermore, other experimental investigations by inelastic neutron scattering in cal-
cium-related compounds (calcite [52], CaCO₃, calcium hydroxide [53], Ca(OH)₂, and brushite
[54, 55], CaHPO₄·2H₂O) show the presence of phonons of similar frequency. However, in the
THz-TDS spectrum, only β-CPP and especially s-HA exhibit a narrow vibrational mode at the
characteristic frequency of 2.1 THz. Similarly, a peak at ~2 THz has been observed in cuttle-
fish bone [56] due to the mineral phase made of aragonite calcium carbonate (CaCO₃). The
normal mode displacements analysis based on our DFT calculations indicates an optical
phonon, in good agreement with previous calculations [46] and phonon measurement [52]. These two studies indicate a very low group velocity $v_g = dw/dk$ for this mode in the Brillouin zone center (Q lower than the inverse micron in our measurements), so that our measured phonon width, or inverse of the phonon lifetime, is poorly sensitive to the crystallite size. Both effects, size and disorder, are related to each other and are responsible for the peak broadening in ns-HA and cannot be disentangled in such measurements.

## Conclusion

In conclusion, we must understand that the presence of a 2.1 THz peak in the absorption spectra is not sufficient to prove the existence of HA crystals, since calcium—oxygen—phosphorous motions occur in a broad variety of crystal structures in a lattice mode. However, THz analysis can provide a simple and quantitative signature of a bone tissue that was heated at a temperature higher than ~750°C, or that contains apatite crystallites larger than 100 nm. Thanks to its ease of implementation, THz-TDS could therefore be envisaged to screen large sample batches, either as a quality control during a graft manufacturing stage, or to evidence/ confirm high temperature heating of bone fragments collected on archaeological sites. The technique could also be of great help in quantifying residuals of HA in biomaterials during their resorption without calcination, as for example in BCP® implants.

## Supporting information

**S1 Fig. Refractive index from THz-TDS.** Absorption coefficient and refractive index in bones at 1.5 THz as a function of heating temperature measured by THz-TDS.
(EPS)

**S2 Fig. FTIR spectroscopy.** FTIR spectra of non stoichiometric hydroxyapatites (ns-HA). The samples were made of pellets of ns-HA powder as received and heated at 1000 °C in KBr.
(EPS)

**S3 Fig. SEM characterization.** SEM images of s-HA, ns-HA as received and heated at 1000°C, CPP-β, β-TCP and TTCP. at three different scales: 100 µm (left), 100 nm (center) and 20 nm (right).
(EPS)

**S4 Fig. EDX analysis.** EDX characterization of three samples: s-HAP, ns-HA as received and ns-HA heated at 1000°C. The intensities are normalized on the Ca peak to highlight the calcium deficiency of the ns-HA with respect to s-HA. After heat treatment, the ratio C/P is changed but does not equal the one of s-HA because of the formation of β-TCP.
(EPS)

## Acknowledgments

We thank the ILL for the use of computational cluster and VASP 5.4 licence.

## Author Contributions

**Conceptualization:** Marie Plazanet, Renato Torre, Christian Rey, Aurelien Gourrier.

**Data curation:** Jordanka Tasseva, Paolo Bartolini, Andrea Taschin.

**Formal analysis:** Marie Plazanet, Jordanka Tasseva, Paolo Bartolini, Andrea Taschin, Mariana Verezhak.
**Funding acquisition:** Marie Plazanet, Renato Torre, Aurelien Gourrier.

**Methodology:** Marie Plazanet, Paolo Bartolini, Andrea Taschin, Renato Torre, Christian Rey, Aurelien Gourrier.

**Project administration:** Aurelien Gourrier.

**Resources:** Christèle Combes, Alessandro Di Michele.

**Software:** Jordanka Tasseva, Paolo Bartolini, Andrea Taschin.

**Supervision:** Paolo Bartolini, Renato Torre.

**Validation:** Christian Rey.

**Visualization:** Mariana Verezhak.

**Writing – original draft:** Marie Plazanet.

**Writing – review & editing:** Renato Torre, Christèle Combes, Christian Rey, Mariana Verezhak, Aurelien Gourrier.

## References

1. Zimmermann G, Moghaddam A Allograft bone matrix versus synthetic bone graft substitutes. Injury. 2011; 42(S2):S16-S21. https://doi.org/10.1016/j.injury.2011.06.199 PMID: 21889142

2. Janicki P, Schmidmaier G What should be the characteristics of the ideal bone graft substitute? Combining scaffolds with growth factors and/or stem cells. Injury. 2011; 42(S2):S77-S81. https://doi.org/10.1016/j.injury.2011.06.014 PMID: 21724186

3. Tadic D, Epple M A thorough physicochemical characterisation of 14 calcium phosphate-based bone substitution materials in comparison to natural bone. Biomaterials. 2004; 25(6):987-994. https://doi.org/10.1016/S0142-9612(03)00621-5 PMID: 14615163

4. Li X, Van Blitterswijk CA, Feng Q, Cui F, Watari F The effect of calcium phosphate microstructure on bone-related cells in vitro. Biomaterials. 2008; 29(33):3306-3316. https://doi.org/10.1016/j.biomaterials.2008.04.039 PMID: 18485476

5. Seebach C, Schultheiss J, Wilhelm K, Frank J, Henrich D Comparison of six bone-graft substitutes regarding to cell seeding efficiency, metabolism and growth behaviour of human mesenchymal stem cells (MSC) in vitro. Injury. 2010; 41(7):731-738. https://doi.org/10.1016/j.injury.2010.02.017 PMID: 20233614

6. Mkukuma LD, Skakle JMS, Gibson IR, Imrie CT, Aspden RM, Hukins DWL Effect of the proportion of organic material in bone on thermal decomposition of bone mineral: an investigation of a variety of bones from different species using thermogravimetric analysis coupled to mass spectrometry, high-temperature X-ray diffraction, and Fourier transform infrared spectroscopy. Calcified Tissue International. 2004; 75(4):321-328. https://doi.org/10.1007/s00223-004-0199-5 PMID: 15549647

7. Haberko K, Bućko MM, Brzezińska-Miecznik J, Haberko M, Mozgawa W, Panz T, et al. Natural hydroxyapatite—its behaviour during heat treatment. Journal of the European Ceramic Society. 2006; 26(4):537-542. https://doi.org/10.1016/j.jeurceramsoc.2005.07.033

8. Ooi CY, Hamdi M, Ramesh S Properties of hydroxyapatite produced by annealing of bovine bone. Ceramics international. 2007; 33(7):1171-1177. https://doi.org/10.1016/j.ceramint.2006.04.001

9. Figueiredo M, Fernando A, Martins G, Freitas J, Judas F, Figueiredo H Effect of the calcination temperature on the composition and microstructure of hydroxyapatite derived from human and animal bone. Ceramics international. 2010; 36(8):2383-2393. https://doi.org/10.1016/j.ceramint.2010.07.016

10. Nasiri-Tabrizi B, Fahami A, Ebrahimi-Kahrizsangi R A comparative study of hydroxyapatite nanostructures produced under different milling conditions and thermal treatment of bovine bone. Journal of Industrial and Engineering Chemistry. 2014; 20(1):245-258. https://doi.org/10.1016/j.jiec.2013.03.041

11. Sofronia AM, Baies R, Anghel EM, Marinescu CA, Tanasescu S Thermal and structural characterization of synthetic and natural nanocrystalline hydroxyapatite. Materials Science and Engineering: C. 2014; 43:153-163. https://doi.org/10.1016/j.msec.2014.07.023

12. Niakan A, Ramesh S, Ganesan P, Tan CY, Purbolaksono J, Chandran H, et al. Sintering behaviour of natural porous hydroxyapatite derived from bovine bone. Ceramics International. 2015; 41(2):3024-3029. https://doi.org/10.1016/j.ceramint.2014.10.138
13. Lebon M, Reiche I, Fröhlich F, Bahain JJ, Falguères C Characterization of archaeological burnt bones: contribution of a new analytical protocol based on derivative FTIR spectroscopy and curve fitting of the v1, v3 PO4 domain. Analytical and Bioanalytical Chemistry. 2008; 392(7-8):1479–1488. https://doi.org/10.1007/s00216-008-2469-y PMID: 18972105

14. Thompson TJU, Gauthier M, Islam M The application of a new method of Fourier Transform Infrared Spectroscopy to the analysis of burned bone. Journal of Archaeological Science. 2009; 36(3):910–914. https://doi.org/10.1016/j.jas.2008.11.013

15. Thompson TJU, Islam M, Piduru K, Marcel A An investigation into the internal and external variables acting on crystallinity index using Fourier Transform Infrared Spectroscopy on unaltered and burned bone. Palaeogeography, Palaeoclimatology, Palaeoecology. 2011; 299(1-2):168–174. https://doi.org/10.1016/j.palaeo.2010.10.044

16. Piga G, Gonçalves D, Thompson TJU, Brunetti A, Malgosa A, Enzo S Understanding the crystallinity indices behavior of burned bones and teeth by ATR-IR and XRD in the presence of bioapatite mixed with other phosphate and carbonate phases. International Journal of Spectroscopy. 2016;4810149. https://doi.org/10.1155/2016/4810149

17. Wopenka B, Pasteris JD. A mineralogical perspective on the apatite in bone. In: Materials Science and Engineering C. vol. 25; 2005. p. 131–143. https://doi.org/10.1016/j.msec.2005.01.008

18. Paschalis EP, Betts F, DiCarlo E, Mendelsohn R, Boskey AL. FTIR microspectroscopic analysis of human iliac crest biopsies from untreated osteoporotic bone. Calcified Tissue International. 1997; 61 (6):487–492. https://doi.org/10.1007/s002239900372 PMID: 9383276

19. Freeman JJ, Wopenka B, Silva MJ, Pasteris JD. Raman spectroscopic detection of changes in bioapatite in mouse femora as a function of age and in vitro fluoride treatment. Calcified Tissue International. 2001; 68(3):156–162. https://doi.org/10.1007/s002230001206 PMID: 11351499

20. Schuetz R, Fix D, Schade U, Aziz EF, Timofeeva N, Weinkamer R, et al. Anisotropy in bone demineralization revealed by polarized far-IR spectroscopy. Molecules (Basel, Switzerland). 2015; 20(4):5835–50. https://doi.org/10.3390/molecules20045835

21. Reiche I, Chadefaux C. Archaeological Bone from Macro- to Nanoscale: Heat-Induced Modifications at Low Temperatures. Journal of Nano Research. 2009; 8:157–172. https://doi.org/10.4028/www.scientific.net/JNanoR.8.157

22. Collins MJ, Nielsen-Marsh CM, Hiller J, Smith CI, Roberts JP, Prigodich RV, et al. The survival of organic matter in bone: a review. Archaeometry. 2002; 44(3):383–394. https://doi.org/10.1111/1475-4754.t01-1-00071

23. Ellingham STD, Thompson TJU, Islam M. The Effect of Soft Tissue on Temperature Estimation from Burnt Bone Using Fourier Transform Infrared Spectroscopy. Journal of Forensic Sciences. 2016; 61 (1):153–159. https://doi.org/10.1111/1556-4029.12855 PMID: 26275238

24. Piga G, Malgosa A, Thompson TJU, Enzo S. A new calibration of the XRD technique for the study of archaeological burned human remains. Journal of Archaeological Science. 2008; 35(8):2171–2178. https://doi.org/10.1016/j.jas.2008.02.003

25. Gourrier A, Chadefaux C, Lemaitre E, Bellot-Gurlet L, Reynolds M, Burghammer M, et al. Nanoscale modifications in the early heating stages of bone are heterogeneous at the microstructural scale. PLOS ONE. 2017; 12(4):e0176179. https://doi.org/10.1371/journal.pone.0176179 PMID: 28423023

26. Kalsbeek N, Richter J. Preservation of Burned Bones: An Investigation of the Effects of Temperature and pH on Hardness. Studies in Conservation. 2006; 51(2):123–138. https://doi.org/10.1179/sic.2006.51.2.123

27. Rogers KD, Daniels P. An X-ray diffraction study of the effects of heat treatment on bone mineral microstructure. Biomaterials. 2002; 23(12):2577–2585. http://dx.doi.org/10.1016/S0142-9612(01)00395-7. PMID: 12033606

28. Woodard JR, Hilldore AJ, Lan SK, Park CJ, Morgan AW, Eurell JAC, et al. The mechanical properties and osteoconductivity of hydroxyapatite bone scaffolds with multi-scale porosity. Biomaterials. 2007; 28 (1):45–54. https://doi.org/10.1016/j.biomaterials.2006.08.021 PMID: 16963118

29. Weiner S, Wagner HD. The Material Bone: Structure-Mechanical Function Relations. Annual Review of Materials Science. 1998; 28(1):271–298. https://doi.org/10.1146/annurev.matsci.28.1.271

30. Bessou M, Chassagne B, Caumes JP, Pradère C, Maire P, Tondusson M, et al. Three-dimensional terahertz computed tomography of human bones. Applied Optics. 2012; 51(28):6738. https://doi.org/10.1364/AO.51.006738 PMID: 23033089

31. Kim SM, Kung P. Terahertz Functional Tissue Imaging. PIERS Proceedings. 2013; p. 434–437.

32. Jung E, Choi HJ, Lim M, Kang H, Park H, Han H, et al. Quantitative analysis of water distribution in human articular cartilage using terahertz time-domain spectroscopy. Biomedical Optics Express. 2012; 3(5):1110. https://doi.org/10.1364/BOE.3.001110 PMID: 22567600
33. Stringer MR, Lund DN, Foulds aP, Uddin a, Berry E, Miles RE, et al. The analysis of human cortical bone by terahertz time-domain spectroscopy. Physics in medicine and biology. 2005; 50(14):3211-3219. https://doi.org/10.1088/0031-9155/50/14/001 PMID: 16177504

34. Yatongchai C, Wren aW, Sundaram SK. Characterization of Hydroxyapatite-Glass Composites Using Terahertz Time-Domain Spectroscopy. Journal of Infrared, Millimeter, and Terahertz Waves. 2015; 36 (1):81-93. https://doi.org/10.1007/s10762-014-0103-y

35. Heughebaert J, Montel G. Preparation de l'orthophosphate tricalcique pur. Bull Soc chim Fr. 1970; 89:2923-2934.

36. Jacob S, Schmitt ML, Schlesinger ME. Use of high-temperature differential scanning calorimetry to investigate the β-α transition in calcium pyrophosphate. Journal of Thermal Analysis and Calorimetry. 2010; 100(3):897-900. https://doi.org/10.1007/s10973-009-0373-0

37. Arends J, Christoffersen J, Christoffersen M, Eckert H, Fowler B, Heughebaert J, et al. A calcium hydroxyapatite precipitated from an aqueous solution: an international multimethod analysis. Journal of Crystal Growth. 1987; 84(3):515-532. https://doi.org/10.1016/0022-0248(87)90284-3

38. Tasseva J, Taschin A, Bartolini P, Striova J, Fontana R, Torre R. Thin layered drawing media probed by THz time-domain spectroscopy. Analyst. 2017; 142:42-47. https://doi.org/10.1039/C6AN02113A

39. Taschin A, Bartolini P, Tasseva J, Torre R. THz time-domain spectroscopic investigations of thin films. Measurement 2018; 118: 282-288. https://doi.org/10.1016/j.measurement.2017.05.074

40. Taschin A, Bartolini P, Tasseva J, Striova J, Fontana R, Riminesi C, et al. Drawing materials studied by THz spectroscopy. Acta-Imeko. 2017; 6:12-17. https://doi.org/10.21014/acta_imeko.v6i3.447

41. Hughes JM, Cameron M, Crowley KD. Structural variations in natural F, OH, and Cl apatites. American Mineralogist. 1989; 74(7-8):870-876.

42. Yashima M, Sakai A, Kamiyama T, Hoshikawa A. Crystal structure analysis of β-tricalcium phosphate Ca3(PO4)2 by neutron powder diffraction. Journal of Solid State Chemistry. 2003; 175(2):272-277. http://dx.doi.org/10.1016/S0022-4596(03)00279-2.

43. Boudin S, Grandin a, Borel MM, Leclaire a, Raveau B. Redetermination of the β-Ca2P2O7 structure. Acta Crystallographica Section C Crystal Structure Communications. 1993; 49(12):2062-2064. https://doi.org/10.1107/S0108270193005608

44. Dickens B, Brown WE, Kruger GJ, Stewart JM. Ca4(PO4)2O, tetracalcium diphosphate monoxide. Crystal structure and relationships to Ca5(PO4)3OH and K3Na(SO4)2. Acta Crystallographica Section B. 1973; 29(10):2046-2056. https://doi.org/10.1107/S0567740873006102

45. Karhánek D, Bučko T, Hafner J. J Phys: Condens Matter. 2010; 22:265006.

46. Bhat SS, Waghmare UV, Ramamurty U. First-Principles Study of Structure, Vibrational, and Elastic Properties of Stoichiometric and Calcium-Deficient Hydroxyapatite. Crystal Growth & Design. 2014; 14 (6):3131-3141. https://doi.org/10.1021/cg5004269

47. De Aza PN, Santos C, Pazo A, de Aza S, Cuscó R, Artús L. Vibrational properties of calcium phosphate compounds. 1. Raman spectrum of b-tricalcium phosphate. Chemical Materials. 1997; 9 (10):15-912.

48. Rey C, Combes C, Drouet C, Grossin D. 1.111-Bioactive Ceramics: Physical Chemistry. vol. 1. Elsevier Ltd.; 2011. Available from: http://dx.doi.org/10.1016/B978-0-08-055294-1.00178-1.

49. Gross Ka, Berndt CC. Thermal processing of hydroxyapatite for coating production, Master of Philosophy in Physics. Journal of Biomedical Materials Research. 1998; 39(4):580-587. https://doi.org/10.1002/(SICI)1097-4636(19980315)39:4%3C580::AID-JBM12%3E3.0.CO;2-B PMID: 9492219

50. Piga G, Solinas G, Thompson TJU, Brunetti A, Malgosa A, Enzo S. Is X-ray diffraction able to distinguish between animal and human bones? Journal of Archaeological Science. 2013; 40(1):778-785. https://doi.org/10.1016/j.jas.2012.07.004

51. Verezhak M, Rauch EF, Véron M, Lancelon-Pin C, Putaux JL, Plazanet M, et al. Ultrafine heat-induced structural perturbations of bone mineral at the individual nanocrystal level. Acta Biomaterialia. 2018; 73:500-508. https://doi.org/10.1016/j.actbio.2018.04.004. PMID: 29649638

52. Cowley ER, Pant AK. Lattice Dynamics of Calcite. Phys Rev B. 1973; 8:4795-4800. https://doi.org/10.1103/PhysRevB.8.4795

53. Baddour-Hadjean R, Fillaux F, Floquet N, Belushkin S, Natkaniec I, Desgranges L, et al. Inelastic neutron scattering study of proton dynamics in Ca(OH)2 at 20 K. Chemical Physics. 1995; 197(1):81-90. https://doi.org/10.1016/0301-0104(95)00155-H

54. Taylor MG, Simkiss K, Parker F, Mitchell PCH, Rg UK, Appleton L, et al. Inelastic neutron scattering studies of synthetic calcium phosphates. Physical Chemistry Chemical Physics. 1999; p. 3141-3144. https://doi.org/10.1039/a902511a
55. Taylor MG, Parker SF, Simkiss K, Mitchell PCH. Bone mineral: evidence for hydroxy groups by inelastic neutron scattering. Physical Chemistry Chemical Physics. 2001; 3(8):1514–1517. https://doi.org/10.1039/b005666i

56. Mizuno M, Fukunaga K. Analysis of tissue condition based on interaction between inorganic and organic matter in cuttlefish bone. Journal of Biological Physics. 2013; 39(1):123–130. https://doi.org/10.1007/s10867-012-9290-8 PMID: 23860837