# Measurement of ultrasonic phase and group velocities in human dental hard tissue

- Sleiman R Ghorayeb
^{1, 2}Email author, - Panagiotis Petrakis
^{1}, - Michael McGrath
^{1}and - Ben A Scheven
^{3}

**1**:5

https://doi.org/10.1186/2050-5736-1-5

© Ghorayeb et al.; licensee BioMed Central Ltd. 2013

**Received: **7 November 2012

**Accepted: **18 February 2013

**Published: **1 May 2013

## Abstract

### Background

The development of ultrasound for use in dental tissues is hampered by the complex, multilayered nature of the teeth. The purpose of this preliminary study was to obtain the phase and group velocities associated with several directions of ultrasonic wave propagation in relation to the tooth structure, which would then lead to the determination of the elastic constants in dental hard tissue. Knowledge of these elastic constants can be used to feed back into numerical models (such as finite element) in order to simulate/predict ultrasonic wave propagation and behavior in the teeth. This will help to optimize ultrasonic protocols as potential noninvasive therapeutic tools for novel dental regenerative therapies.

### Methods

An extracted human second molar was used to determine time-of-flight information from A-scan signatures obtained at various angles of inclination and rotation using a scanning acoustic microscope at 10 MHz. Phase and group velocities and associated slowness curves were calculated in order to determine the independent elastic constants in the human teeth.

### Results

Results show that as the tooth was inclined at three azimuthal angles (Θ_{in} = 0°, 15°, and 30°) and rotated from Φ_{in} = 0° to 360° in order to cover the whole perimeter of the tooth, slowness curves constructed from the computed phase and group velocities versus angle of rotation confirm the inhomogeneous and anisotropic nature of the tooth as indicated by the nonuniform appearance of uneven circular shape patterns of the measurements when compared to those produced in a control isotropic fused quartz sample.

### Conclusions

This study demonstrates that phase and group velocities of ultrasound as determined by acoustic microscopy change and are dependent on the direction of the tooth structure. Thus, these results confirm that the tooth is indeed a multilayered anisotropic structure underscoring that there is no single elastic constant sufficient to represent the complex structure of the tooth. The findings underline the importance to take into account these crucial characteristics in order to develop and optimize therapeutic as well as diagnostic applications of ultrasound in dental tissue repair, and further studies are warranted to analyze ultrasound transmission at various frequencies and intensities in different teeth to develop appropriate models for ultrasound biophysical behavior in dental tissues.

## Keywords

## Background

Ultrasonic waves such as those used routinely during medical examinations and prognoses propagate in a very complex manner. This level of complexity holds true even when these sound waves propagate in ideal, homogeneous, and uniform media. Despite the fact that diagnostic ultrasound techniques are based on well-established phenomena, there remains the area of dentistry where ultrasound has not played yet a very important role. This is due to the fact that ultrasonic propagation in inhomogeneous/anisotropic media results in further degradations in acoustic field. Scattering effects at the microscale level as well as at the multilayer interfaces are responsible for these degradations.

In earlier diagnostic ultrasound experiments performed on the human teeth, an ultrasonic pulse-echo system was able to detect the enamel-dentin junction and the dentin-pulp interface [1]. Ultrasonic imaging techniques have also been applied to other areas of dentistry, specifically prosthodontics. A B-mode ultrasonic system for determining the thickness of the masticatory mucosa for denture construction has also been developed [2], and the use of the scanning acoustic microscope to image the elastic properties of carious human dental enamel has been described [3]. Acoustic microscopy has been successful for imaging lesions in tooth enamel and detecting the enamel-dentin interface [4]. Finite element models [5, 6] as well as PSpice (OrCAD™, Cadence Design Systems, Inc., San Jose, CA, USA) transmission line (*T*-line) models [7] used to understand ultrasonic wave propagation in the teeth have been established. In order to confirm predictions made by these latter finite element and *T*-line models, experimental data was gathered [8] for the purpose of assessing the viability of ultrasound to detect cavities and fractures in the extracted human teeth. The method showed great promise as a potential improvement upon current dental imaging systems.

For years, the main applications of therapeutic ultrasound in dentistry have been mostly limited to oral surface cleaning (i.e., removal of plaque and calculus) or root canal treatment [9]. In recent years, however, therapeutic ultrasound studies performed on dental tissue have taken a different path [10, 11]. The goal of this research is to develop novel ultrasound therapies that can be easily applied to the teeth to promote tissue healing, more specifically tooth regeneration. Preliminary results showed that this idea of repairing the teeth with ultrasound has potential. The tooth is a hard mineralized structure with a living soft tissue core capable to respond to outside stimuli including physical forces. Research has demonstrated that ultrasound in the low frequency range, generally used in dentistry for dental scaling, stimulated gene expression and the production of tissue regulatory growth factors that are considered to be crucial for odontoblast differentiation, activity, and dentin repair [12, 13]. More recently, it was demonstrated that low-frequency ultrasound stimulated differentiation and *in vitro* mineralization in odontoblast-like cultures [11]. These findings are novel and exciting and support the working hypothesis that ultrasound may be an effective noninvasive tool for novel dental regenerative therapies [14]. Interestingly, these findings were recently corroborated in *in vitro* tooth slice organ cultures where low-intensity pulsed ultrasound was shown to increase pre-dentin thickness and sub-odontoblast cell numbers [10, 15].

However, to develop ultrasound as a therapeutic application, a thorough understanding is required regarding ultrasonic wave propagation and intensity distribution in teeth structures. The tooth is a highly anisotropic medium. This quality stems from two different criteria—layer consistency and structural component. The enamel portion of a tooth mainly consists of calcium hydroxyapatite [16]; the dentin represents a collagenous mineralized substance much like the bone and contains tubular canals, and the tooth's central cavity (pulp) contains a loose connective tissue with various cell types, blood vessels, and nerves. Odontoblasts surround the dental pulp with cellular processes extending into the dentinal tubules which may also contain neurite endings. In addition, the root is enclosed by a thin layer of bone-like material called cementum, which is surrounded by a periodontal ligament. The latter tissue is highly cellular and contains bundles of thick collagenous fibers responsible for attaching the tooth to the alveolar jaw bone. So, as it may be appreciated, the tooth is composed of a very complex mesh of intricate layers, each having different chemical and structural characteristics than the one adjacent to it. As a result of all these effects, ultrasonic waves propagating within a tooth would undergo a fair amount of mode conversion. This, in a way, is similar to the assessment of the structural integrity of anisotropic components used in aircrafts and submarines whereby elastic constants of the host material need to be determined.

In this paper, we propose and test a proof-of-concept methodology based on previous studies [17–21] that would eventually lead to the determination of the elastic constants in dental hard tissue. The study used an extracted human tooth to determine phase and group velocities from measured time-of-flight (TOF) of an ultrasonic wave as it propagates along arbitrary directions through the tooth crown of the specimen using a pulse-echo setup at different inclination angles of incidence.

## Methods

All studies were approved by and in compliance with the guidelines set by the institutional review board committee at Hofstra University (Hempstead, NY, USA).

### Anisotropy consideration

*V*

_{L}), shear horizontal (

*V*

_{SH}), and shear vertical (

*V*

_{SV})) are determined along these directions. Note that all diagonal elastic constants can be determined by propagation waves in symmetric directions. The waves traveling along the lower-symmetry directions are not pure modes. The behavior of elastic waves propagating in anisotropic materials [17] is governed by the following:

*C*

_{ ijkl }is the elastic constant,

*u*

_{ i }is the displacement, and ρ is the mass density of the material. It is assumed that a plane wave solution

*u*

_{ i }, with wave vector

*k*, frequency ω, and amplitude

*A*

_{ i }, of the form

*u*

_{ i }=

*A*

_{ i }

*e*

^{i(kx − ωt)}exists and can be substituted back into the equation of motion (Equation 2) in order to obtain Kelvin-Christoffel elastic stiffness equation as follows:

*Γ*

_{ ik }=

*C*

_{ ijkl }

*α*

_{ i }

*α*

_{ j }are Christoffel elastic stiffnesses [17],

*V*is the phase velocity, and α

_{ i }and α

_{ j }are the direction cosines of propagating wave. A general solution to (Equation 3) gives relationships between the phase velocities

*V*and

*C*

_{ ij }for the wave propagation along arbitrary directions. As an example, in the <0 1 1> direction, we should have the following:

*C*

_{ ii }is known and ρ

*V*

^{2}is measured, then

*C*

_{23}can be calculated as follows:

*C*

_{13}and

*C*

_{12}as follows:

*C*

_{11},

*C*

_{22},

*C*

_{33},

*C*

_{44},

*C*

_{55}, and

*C*

_{66}and measure the wave velocities along the <0 1 1>, <1 0 1>, and <1 1 0> directions for the off-diagonal elements

*C*

_{23},

*C*

_{13}, and

*C*

_{12}. One should appreciate that this is not a trivial procedure. In order to measure the velocities along the aforementioned directions, four similar (if not exact) teeth specimens need to be prepared as follows:

- 1.
All surfaces are perpendicular to the symmetry direction in such a way that the length in each direction is at least five longitudinal wavelengths for proper measurement, and all surfaces are flat.

- 2.
Largest surface perpendicular to the <0 1 1> direction and upper and lower surfaces are flat.

- 3.
Largest surface perpendicular to the <1 0 1> direction and upper and lower surfaces are flat.

- 4.
Largest surface perpendicular to the <1 1 0> direction and upper and lower surfaces are flat.

However, as mentioned earlier, this study is only a proof-of-concept and was performed on only one tooth sample. Future studies will consider the aforementioned implementation to measure the elastic constants from obtained phase velocities.

### Experimental setup

*d*= 6.25 mm in length (as opposed to the 15-mm original length) and 10 mm in width (measured across the breadth of the crown area) as shown in Figure 1. However, some of the indentations due to natural occlusions (OC in Figure 1) were left. Also, one should keep in mind that due to the flattening process, the majority of the enamel layer has been lost, exposing only the dentin portion of the tooth. This in fact is advantageous since the major anisotropy/inhomogeneity is present in the dentin due to the presence of collagenous mineralized substance and tubular canals, as mentioned earlier.

*Φ*

_{in}and azimuthal/inclination

*Θ*

_{in}). During data acquisition, it is important to note that the tooth was scanned with the transducer perpendicular to the tooth surface of the crown as shown in Figure 3.

*Φ*

_{in}and

*Θ*

_{in}), the phase velocity [17] is given by the following:

*t*is the one-way TOF,

*d*is the thickness of the tooth (as shown in Figure 1), and

*V*

_{ w }is the phase velocity in the PBS solution. Furthermore, using Snell's law, the angle of refraction

*Θ*

_{ r }can be determined from the following:

From Equations 9 and 10, the magnitude and directions of phase velocities for all modes of propagation inside the tooth can be determined for any arbitrary angle of incidence, *Θ*_{in}. Once this is done, an inversion routine based on the Newton-Rhapson method [17] can be implemented on the above characteristic polynomial to determine the elastic constants. Knowledge of these elastic constants can be used to feed back into models (such as finite element) in order to simulate/predict ultrasonic wave propagation and behavior in the teeth.

*Θ*

_{ r }has been calculated for each phase velocity, one can determine the various group velocities

*V*

_{ g }which correspond to these respective phase velocities at the calculated angles

*Θ*

_{ r }. These group velocities represent the actual family of longitudinal stress-wave velocities in those particular directions within the tooth. Group velocities can be determined by projecting their vectors, on the surface of the slowness curves (as will be shown later), normal to the corresponding phase velocities. Using simple trigonometry, the group velocities can then be computed from the following:

## Results and discussion

*Θ*

_{in}= 0°, 15°, and 30°. For each of these three angles, the platform was rotated from

*Φ*

_{in}= 0° to about 360° in order to cover the whole perimeter of the tooth. The phase velocity in the PBS solution (

*V*

_{ w }) is assumed to be similar to that in water (1,500 m/s). For each set of incident angle (

*Φ*

_{in},

*Θ*

_{in}), TOF and the phase velocity

*V*

_{ p }were determined. Polar plots of TOF and slowness curves from the computed phase velocities are then constructed. In order to confirm the experimental measurement system, TOF measurements as a function of

*Φ*

_{in}and

*Θ*

_{in}were made on an isotropic, homogeneous square piece of fused quartz plate of 4-mm thickness. The polar plot of TOF values Δ

*t*((

*Φ*

_{in}, 0° → 2

*π*),

*Θ*

_{in}) determined from around the quartz plate is perfectly circular for propagation of longitudinal wave, with a constant longitudinal velocity of about 5,714 m/s. This observation confirms the homogenous characteristic of the control quartz plate. Figure 5 shows the results for both TOF and longitudinal velocities obtained for the quartz sample at an angle of inclination (

*Θ*

_{in}= 0°). Figure 6 shows the TOF measurements for the tooth for all three angles of inclination (

*Θ*

_{in}= 0°, 15°, and 30°), and Figure 7 illustrates the corresponding phase velocities in the tooth for the same angles.

The observations highlight polar/slowness curves (Figures 6 and 7) which confirm the inhomogeneous and anisotropic nature of the tooth. This is indicated by the nonuniform appearance of uneven circular shape patterns of both TOF and *V*_{
p
} measurements when compared with the perfect circular measurements made in the control isotropic fused quartz sample shown in Figure 5. This outcome has obvious implications for the development of ultrasound as a diagnostic tool, and in this respect, further analyses need to be conducted to determine the phase velocity of ultrasound at different frequencies.

*V*

_{ p }measurements shown in Figures 7 and 8, the group velocities can be calculated. Figure 9 shows the corresponding polar plots.

As mentioned earlier, since the flattening process of the tooth caused the majority of the enamel layer to be lost, exposing only the dentin portion, the average group velocities attained were approximately 3,178, 3,088, and 3,080 m/s for each of the three azimuthal angles (*Θ*_{in} = 0°, 15°, and 30°), respectively. These values confirm what was reported in earlier studies [1, 5–8].

Moreover, regarding the potential use of ultrasound for dental tissue repair, it is important to relate the physical behavior of ultrasound through the different layers of the teeth to those in a viable pulp cavity and relate these characteristics to localized ultrasound intensity in order to assess the biological effects of ultrasound in the pulp region. In other words, knowledge of the elastic constants can be used to feed back into numerical models (such as finite element) in order to simulate/predict ultrasonic wave propagation and behavior in the teeth. This will help optimize ultrasonic protocols as potential noninvasive tools for novel dental regenerative therapies.

## Conclusions

This proof-of-concept study demonstrated that phase velocities (and corresponding group velocities) of ultrasound as determined by acoustic microscopy change and are dependent on the direction of the tooth structure. The data underscore the anisotropic structure and demonstrate that no single elastic constant will be sufficient to represent the complex nature of the tooth. The findings have obvious implications for the development of diagnostic and therapeutic ultrasonic devices for dental tissues, and further correlation studies that will examine other samples are warranted to analyze ultrasound transmission at various frequencies and intensities in different teeth to develop appropriate models for ultrasound biophysical behavior in dental tissues.

## Declarations

### Acknowledgements

This research was performed within the framework of the Capstone Senior Design component in the School of Engineering and Applied Sciences, Ultrasound Research Laboratory, at Hofstra University.

## Authors’ Affiliations

## References

- Barber FE, Lees S, Lobene RR: Ultrasonic pulse-echo measurements in teeth. Arch Oral Biol. 1969, 14: 745-60. 10.1016/0003-9969(69)90166-6.View ArticlePubMedGoogle Scholar
- Uchida H, Kobayashi K, Nagao M: Measurement
*in vivo*of masticatory mucosal thickness with 20 MHz B-mode ultrasonic diagnostic equipment. J Dent Res. 1989, 68: 95-100. 10.1177/00220345890680021501.View ArticlePubMedGoogle Scholar - Peck SD, Rooe JM, Briggs GAD: Studies on sound and carious enamel with the quantitative acoustic microscope. J Dent Res. 1989, 68: 107-12. 10.1177/00220345890680020201.View ArticlePubMedGoogle Scholar
- Briggs A: Acoustic Microscopy. 1992, New York: Oxford University Press, 165-202.Google Scholar
- Ghorayeb SR, Xue T, Lord W: A finite element study of ultrasonic wave propagation in a tooth phantom. J Dent Res. 1998, 77 (1): 39-49. 10.1177/00220345980770010301.View ArticlePubMedGoogle Scholar
- Ghorayeb SR, Lord W, Xue T: In
*Proceedings of the Third Annual ASNT Research Symposium*. The study of ultrasonic wave propagation in a tooth phantom. 1994, New Orleans: ASNT, 27-30.Google Scholar - Ghorayeb SR, Maione E, La Magna V: Modeling of ultrasonic wave propagation in teeth using PSpice: a comparison with finite element models. IEEE Trans Ultrason Ferroelec Freq Contr. 2001, 48 (4): 1124-31. 10.1109/58.935731.View ArticleGoogle Scholar
- Ghorayeb SR, Valle T: Experimental evaluation of human teeth using noninvasive ultrasound: echodentography. IEEE Trans Ultrason Ferroelec Freq Contr. 2002, 49 (10): 1437-43.View ArticleGoogle Scholar
- Walmsley AD: Applications of ultrasound in dentistry. Ultrasound Med Biol. 1988, 14: 7-14. 10.1016/0301-5629(88)90158-5.View ArticlePubMedGoogle Scholar
- El-Bialy T, Lam B, Aldaghreer S, Sloan AJ: The effect of low intensity pulsed ultrasound in a 3D
*ex vivo*orthodontic model. J Dent. 2011, 39 (10): 693-9. 10.1016/j.jdent.2011.08.001.View ArticlePubMedGoogle Scholar - Scheven BAA, Shelton RM, Cooper PR, Walmsley AD, Smith AJ: Therapeutic ultrasound for dental tissue repair. Med Hypoth. 2009, 73: 591-3. 10.1016/j.mehy.2009.05.032.View ArticleGoogle Scholar
- Scheven BAA, Millard JL, Cooper PR, Lea SC, Walmsley AD, Smith AJ: Short-term
*in vitro*effects of low frequency ultrasound on odontoblast-like cells. Ultrasound Med Biol. 2007, 33: 1475-82. 10.1016/j.ultrasmedbio.2007.03.010.View ArticlePubMedGoogle Scholar - Scheven BAA, Man J, Millard JL, Cooper PR, Lea SC, Walmsley AD, Smith AJ: VEGF and odontoblast-like cells: stimulation by low frequency ultrasound. Arch Oral Biol. 2009, 54: 185-91. 10.1016/j.archoralbio.2008.09.008.View ArticlePubMedGoogle Scholar
- Man J, Shelton RM, Cooper PR, Scheven BAA: Low-intensity low-frequency ultrasound promotes proliferation and differentiation of odontoblast-like cells. J Endod. 2012, 38: 608-13. 10.1016/j.joen.2012.01.015.View ArticlePubMedGoogle Scholar
- Al-Daghreer S, Doschak M, Sloan AJ, Major PW, Heo G, Scurtescu C, Tsui YY, El-Bialy T: Long term effect of low intensity pulsed ultrasound on a human tooth slice organ culture. Arch Oral Biol. 2012, 57 (6): 760-8. 10.1016/j.archoralbio.2011.11.010.View ArticlePubMedGoogle Scholar
- Hole JW: Human Anatomy and Physiology. 1984, Dubuque: WCB Publishers, 485-7.Google Scholar
- Batra NK, Simmonds KE, Mignogna RB: Determination of elastic constants of anisotropic materials from propagation of ultrasonic waves along arbitrary symmetry directions. 14th World Conf On Non Destruct Testing. 1996, 4: 2275-8.Google Scholar
- Rose JL, Ditri JJ, Huang Y, Dandekar DP, Chou S-C: One-sided ultrasonic inspection technique for the elastic constant determination of advanced anisotropic materials. J Nondestr Eval. 1991, 10 (4): 159-66. 10.1007/BF00567098.View ArticleGoogle Scholar
- Chedid-Helou FA, Hemann JH: Mathematical modeling of wave propagation in anisotropic media. Mat Eval. 1991, 49: 708-15.Google Scholar
- Mignogna RB: General analytic solution of the wave equation for anisotropic materials using symbolic manipulation, in review of Progress in Quantitative Nondestructive. Evaluation. 1989, 8A: 133-40.Google Scholar
- Pearson LH, Murri WJ: Measurement of ultrasonic wavespeeds in off-axis directions of composite materials. Rev Prog Quant Nondestructive Eval. 1987, 6B: 1093-101.View ArticleGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.