9 July 2012 Optodynamic energy-conversion efficiency during an Er:YAG-laser-pulse delivery into a liquid through different fiber-tip geometries
Author Affiliations +
J. of Biomedical Optics, 17(7), 075006 (2012). doi:10.1117/1.JBO.17.7.075006
When an erbium-laser pulse is directed into water through a small-diameter fiber tip (FT), the absorption of the laser energy superheats the water and its boiling induces a vapor bubble. We present the influence of different FT geometries and pulse parameters on the vapor-bubble dynamics. In our investigation, we use a free-running erbium: yttrium aluminum garnet (Er:YAG) (λ = 2.94 μm) laser that was designed for laser dentistry. Its pulse is directed into the water through FTs with a flat and conical geometry. Our results show that in the case of the conical FT, a spherical bubble is induced, while a channel-like bubble develops for the flat FT. The ratio between the mechanical energy of the liquid medium and the pulse energy, which we call the optodynamic energy-conversion efficiency, is examined using shadow photography. The results indicate that this efficiency is significantly larger when a conical FT is used and it increases with increasing pulse energy and decreasing pulse duration. The spherical bubbles are compared with the Rayleigh model in order to present the influence of the pulse duration on the dynamics of the bubble's expansion.
Gregorčič, Jezeršek, and Možina: Optodynamic energy-conversion efficiency during an Er:YAG-laser-pulse delivery into a liquid through different fiber-tip geometries



Lasers with wavelengths in the interval 2 to 3 μm are widely used in various medical applications where clinical treatments are usually performed in a liquid environment. Of all the infrared lasers, the erbium: yttrium aluminum garnet (Er:YAG) laser, emitting pulsed light at 2.94 μm, has the highest absorption in water.6 Thus, its radiation delivered through a small-diameter fiber tip (FT) can be used for a wide range of laser-assisted medical applications, such as lithotripsy,7 ophthalmology,8,9 and dentistry.1011.

Due to the strong absorption of erbium radiation in water, a vapor bubble is formed during the erbium-pulse delivery through an optical fiber that is not in contact with any tissue.18,19 The cavitation-bubble dynamics then affects the tissue-ablation efficiency in two ways.20 First, the bubble creates a vapor-water channel for the laser radiation from the FT to the tissue.21,22 Second, the bubble dynamics2324.25.26.27 may contribute to the material’s removal and collateral effects.28,29 The described phenomena represent a typical optodynamic process, where the energy of the laser pulse is coupled with the mechanical energy of the liquid medium.30 The characterization of these phenomena is therefore important for the optimization of biomedical processes based on laser ablation in liquids.

The main goal of our investigation is to study the influence of the FT geometry and the erbium-laser-beam parameters on the optodynamic (OD) energy-conversion efficiency, which we define as the mechanical energy of optically induced motion divided by the pulse energy. Here, the radiation of the Er:YAG laser is delivered into an infinite liquid through FTs with a flat and conical geometry. In many biomedical applications,2,7,10,18,22,31 the flat FTs are applied. In this case, the laser radiation exits the optical fiber in forward direction with a relatively small divergence. On the other hand, in the past few years, the conical FTs are used in dentistry for cleaning, shaping, and sterilization of root canals.1415.16.17 In this case, the conical FT is placed into the coronal portion of the pulpa chamber, filled with a liquid, such as water or sodium hypochlorite.12,17 For these dental treatments, there is a need to maximize the mechanical effects with minimal thermal side effects.16

Using shadow photography, we show that the shape of the bubble and its dynamics depend on the FT geometry. In the case of the flat FT, a channel-like bubble is induced,3,9,21,22,32 while our measurements show that a spherical bubble is developed when erbium radiation is delivered through the conical FT. Thus, we compare the measurements obtained using the conical FT with the Rayleigh model33 in dimensionless form to analyze the influence of the pulse duration on the bubble’s dynamics.

Our measurements indicate that the stress waves are not induced during the bubble’s formation. Therefore, in the case of the free-running mode Er:YAG pulses, the total mechanical energy of the liquid equals the energy of the first bubble’s expansion. With these measurements, we can explain how the FT’s geometry, the laser-pulse energy, and the pulse duration influence the cavitation bubble size, which further determines the acoustic transient amplitudes.32,3435.36.37 The experiments were performed with a commercially available laser for clinical applications in laser dentistry. Thus, the experimental parameters are very similar to the parameters in clinical practice.


Materials and Methods


Experimental Setup

Our experimental setup is schematically shown in Fig. 1. A free-running Er:YAG laser, designed for laser dentistry (λ=2.940μm, Fotona d.d., Slovenia, LightWalker), was used as the excitation source. We used pulses with a duration in the range tL=40160μs and pulse energies (at the output of a FT) between EL=5100mJ. The repetition was set to 5 Hz. The excitation-pulse power as a function of time is detected by a 60-MHz InAs-photodiode (PD-InAs) placed behind the rear mirror of the laser resonator (LR). The delivery of the laser pulse into the vessel with distilled water is performed through an articulated arm. Within the handpiece, the pulse is optically coupled into an interchangeable FT. The FT’s end is placed h=10mm below the water level, as is schematically shown in the inset of Fig. 1. In our experiments, we used FTs with flat and conical terminations. Their geometries are schematically presented in Fig. 2. The FTs are made from fused silica with a low OH-content core, a fluorine silica cladding and a polyimide jacket. Their numerical aperture in air is NA=0.22. Thus, the light exits the flat FT conically with a maximum angle of ψout=12.7deg. Using the core refractive index nc=1.421, the maximum inner angle is calculated to be ψin=8.9deg. The conical FT’s termination subtends an angle at γ=27deg.

The optodynamic phenomena in a liquid24,30,38,39 induced by erbium radiation delivered through different FTs were examined using shadowgraphy. Here, a frequency-doubled Nd:YAG laser (Ekspla, Lithuania, PL2250-SH-TH) emitting green (λ=532nm) pulses with a duration of 30 ps was used as the illumination source. The probe is partially reflected by a dielectric plate to the 1-GHz Si-photodiode (PD-Si). A beam expander (BE) is placed in front of the illumination area where the laser-pulse-water interaction takes place. The image is captured through a microscope equipped with a charge-coupled device (CCD) camera (Basler AG, Germany, scA1400-17fm, 1.4 Mpx). The spatial resolution of our optical system was in the range of κ=3.65.7μm/pixel. For a particular series of measurements, the κ value was determined using the images of a calibrating pattern.

Fig. 1

Schematic view of the experimental setup.


Fig. 2

The geometry of the fiber tips.


The lasers and the camera were synchronized with a signal generator (SG; Tektronix, US, AFG 3102) connected to a personal computer (PC). The experimental setup is automatically controlled with custom-developed software that also enables the setting of the excitation-laser parameters, the data acquisition from a digital oscilloscope (OSC; LeCroy, US, 600 MHz Wave Runner 64MXi-A), image acquisition as well as data and image processing.


Data Analysis

Using shadowgraphy, we measured the vapor-bubble dynamics for various parameters EL and tL. A sequence of images from multiple events was recorded by varying the time delay between the excitation pulse and the probe pulse. This delay is set by the SG and measured from the signals of the photodiodes PD-InAs and PD-Si.

The bubble locally changes the refractive index of the liquid. Consequently, a refractive-index gradient appears on its wall, representing the boundary between the steam and the surrounding water.38 This results in the deflection of the probe pulse casting a shadow,40 where the FT and cavitation bubble are visible as a black area on a bright background (e.g., see Fig. 3). The volume of the bubble during its expansion and collapse is estimated using the image processing of the shadowgraphs captured by the CCD. Here, we assume that the cavitation bubble is cylindrically symmetrical around the FT axis.

Fig. 3

Software detection of vertical edges. The yellow points show the very left and the very right edges of an individual line of the image. (a) The FT before the bubble formation. (b) A fiber tip with a cavitation bubble reaching its maximum volume.


First, the grayscale (8-bit) image is converted into a black-and-white (1-bit) image. In each row, l of the 1-bit image the edges are detected using horizontal gradients. Then, the radius rl is calculated as:


where, xRl,Ll, are the x positions of the very left and the very right edge (the yellow points in Fig. 3), separately. From the radii, rl, the volume from the selected area is obtained by summing the one-pixel-wide slices:


where, N is the height of the image in pixels. Finally, the bubble’s volume VB is calculated as:



Here, the volume of the FT, VFT, is determined from the images before the bubble’s formation [e.g., see Fig. 3(a)] using Eqs. (1) and (2).

The pulse energy EL and its duration tL are simultaneously measured during the experiments using the signal U(t) from the photodiode PD-InAs. Assuming a linear response of the photodiode, the pulse power P(t) is proportional to the measured voltage U(t):


and the pulse energy as a function of time is obtained as:



The proportionality constant C is previously determined by employing a calibrated energy meter (Ophir Optronics Inc., Israel, Smart head to USB interface with a pyroelectric sensor PE50BF-DIF-SH-U2) at the output of the FT in air.

Figure 4 shows two typical graphs of power as a function of time for excitation pulses of the same energy EL=20mJ, but different pulse durations: tL=60μs (the solid curve) and tL=120μs (the broken curve). The curves show an average of 50 measurements and the shaded areas are within one standard deviation from the mean.

Fig. 4

Pulse power (an average of 50 measurements) as a function of time for laser energy EL=20mJ, and pulse durations tL=60μs (the solid curve) and tL=120μs (the broken curve), respectively. The shaded region shows a standard deviation.



Results and Discussion


Bubble Formation

The erbium radiation (λ=2.94μm) delivered into liquid water through a submerged FT is completely absorbed right beside the FT due to the high absorption coefficient μa=1.247×106m1.6 Thus, the water is instantly heated over the boiling point and a vapor bubble starts to form at the FT [e.g., see the second images in Fig. 5(a) and 5(b)] within the time interval tvf after the excitation-pulse initiation.

Fig. 5

Typical images of vapor-bubble formation in water in the first 6 μs after the beginning of the excitation-pulse radiation (EL=20mJ, tL=100μs) delivered through the (a) flat and (b) conically terminated FT. The scale is shown by a bar on the first image.


Figure 5 shows typical images of the vapor-bubble formation during the first 6 μs after the beginning of the excitation pulse (EL=20mJ, tL=100μs). The time of the vapor formation (tvf) depends on the pulse energy and the pulse duration. However, in the case of the laser parameters used in our experiments, we found from the shadowgraphs that the vapor formation starts within the time interval tvf<1.5μs, i.e., during the first spike of the laser. This observation is consistent with the experimental results presented by Brinkmann et al.1


Flat fiber tip

When a FT with the flat end [Fig. 5(a)] is used, the light exits the fiber within the cone defined by the angle ψout (see Fig. 2). Owing to the strong absorption of the erbium radiation in water, the heated layer of a liquid has a thickness of 1/μa, and the initially heated volume can be approximated by Vvf=πRc2μa, where, Rc stands for the radius of the FT core. The characteristic thermal diffusion time, td, is defined as:41


where D=1.4×107m2s1 stands for the thermal diffusivity of water. In our case tvf<td, so the vapor is induced under adiabatic conditions. Therefore, the average local temperature within the heated volume Vvf before the vapor formation can be estimated as:



In the above equation, T0=298K stands for the initial temperature of the water, ρ is the density of the water, cv is the specific heat capacity at a constant volume, EL(tvf) is the irradiated laser energy before the vapor formation, and the factor 0.63 is the portion of energy EL(tvf) absorbed in the layer with the thickness 1/μa.

Figure 6(a) shows a typical temporal power-evolution of the excitation pulse (after the FT) for the case in Fig. 5. The energy as a function of time for the first 14 μs [the shaded area in Fig. 6(a)] is calculated using Eq. (5) and shown in Fig. 6(b). The red point at tvf=1.2μs in Fig. 6(b) corresponds to the time of the boiling point determined from the second image in Fig. 5(a), where the first signs of the bubble formation are visible. Here, a series of images captured in the first 2 μs and separated by 0.2 μs were analyzed. However, an on-line determination of the boiling point can be performed using an optical monitoring proposed by Brinkmann et al.1

Fig. 6

(a) Pulse power as a function of time for the laser energy EL=20mJ and pulse duration tL=100μs. (b) Pulse energy as a function of time for the first 14 μs [the shaded area in (a)]. The red point shows the time when a layer of vapor in the second image of Fig. 5(a) is induced.


Using Eq. (7) and the irradiated laser energy EL(tvf)=0.14mJ [the red point in Fig. 6(b)], we estimate for the flat FT that the average water temperature within the volume Vvf before vaporization is around 230°C. Thus, the bubble is formed by the explosive boiling20 of a superheated water layer.


Conical fiber tip

The situation is more complex when a conical FT is used. Because of the conical geometry, the light experiences several reflections and is gathered at the bottom of the FT’s conical end, as explained below.

Here, the light enters the conical termination within the cone defined by the maximum inner angle ψin (see Fig. 2). This angle is the same as in the case of the flat geometry. Then it reaches the conical-end boundary making an angle α0=90°γ2δ with the boundary’s normal (see Fig. 7). The angle between the entering light and the FT axis, δ, is in the range δ[ψin,+ψin].

Fig. 7

A schematic presentation of internal reflections in the conical FT. For the sake of clarity, the transmitted beams, defined by the Snell and Fresnel laws, are not drawn.


Depending on the refractive indices of the fiber core (nc) and the surrounding medium (nm), and the angles γ and δ, the light is totally or partially reflected. The critical angle for total internal reflection αtot is defined as arcsin(nm/nc). All the rays that are incident at angles greater than αtot are totally reflected and bounded within the FT. However, after the first incidence at an angle α0 with respect to the boundary’s normal, the reflected light is incident upon the conical-end boundary again, but with the smaller angle α1. This repeats several times and the angle of incidence αm after the m-th reflection decreases as:



After the mmax-th reflection, the reflected light returns into the FT (see Fig. 7). During each reflection, the portion of reflected light is defined by the Fresnel laws.9

In our case, the first incidence angle α0 exceeds 67.6 deg, while the angle of the total internal reflection for the silica-water interface equals αtotSW=63.7deg, using the refractive index of water nw=1.274 at λ=2.94μm.6 Thus, the angle of the first incidence α0 is larger than αtotSW, so the whole light is first totally reflected from the conical-end boundary. For this reason, the light gathers at the bottom part of the conical FT and exits the FT at different angles defined by Snell’s law.9,42 Therefore, the power density at the exit of the conical FT is increased in comparison with the flat FT.

Although the condition for total internal reflection is fulfilled for the first incidence of laser radiation upon the silica-water boundary, a vapor layer is observed at t>3μs along the whole conical-termination [see Fig. 5(b)]. We believe that this happens because the conical termination of the FT is not perfectly polished, and due to the evanescent field. Therefore, a portion of the light exits the FT also during the first incidence and turns a thin layer of water into vapor. When the conical-end of the FT is surrounded by a vapor-layer [see Fig. 5(b)] induced by the absorption of Er:YAG light, the angle of total internal reflection is reduced to αtotSV=44.7deg. Moreover, according to the Fresnel laws, the portion of reflected light increases due to the greater difference in the refractive indices nc and nm.


Bubble’s Dynamics

A typical sequence of images of a bubble captured at different times after the beginning of a free-running erbium pulse (EL=40mJ, tL=145μs) is presented in Fig. 8. Our experimental results show that a channel-like bubble is formed in the case of the flat FT [Fig. 8(a)], while a spherical bubble develops when a conical FT is used [Fig. 8(b)]. We assume that the shape of the bubble is correlated with the emission profile at the end of the FT. The conical termination of the optical fiber causes a more divergent output profile,14 which leads to a spherical bubble.

Fig. 8

Typical images of vapor bubbles induced by erbium radiation (EL=40mJ, tL=145μs) delivered through different fiber tips. (a) Flat FT: the bubble has a channel-like shape. (b) Conical FT: a spherical bubble is developed. The scale is shown in the upper-right corner of images A1. The horizontal arrows indicate the FT’s end.


As we have already explained in the previous subsection, the laser-energy deposition due to the strong absorption of the erbium radiation causes superheated water and its boiling induces a vapor bubble. The expansion, driven by rapid heating, can generate thermo-elastic stress waves.35 This happens when the laser energy is deposited into a medium within a time that is shorter than the characteristic time ts for the propagation of stress waves across the irradiated volume:41,43



In the above equation, cs=1,480ms1 stands for the sound velocity in the water. In our case, the vapor is induced in the time interval tvf<1.5μs, which is three orders of magnitude larger than ts. Therefore, the stress-confinement conditions are not achieved, so the stress waves are not induced. This is proved by the shadowgraphic observation, where shock waves are not observed during the first bubble expansion [e.g., see the images A1-C1 in Fig. 3(a) and 3(b)]. This finding is also consistent with the pressure measurements performed with a hydrophone presented in the paper by Frenz et al.37 However, their laser had a wavelength of λ=2.12μm, so the characteristic time ts=0.2μs in their case is comparable to the time of the vapor formation. For this reason, they observed a slow increase in the pressure during the first bubble expansion.

After the laser-induced boiling, the high pressure of the vapor leads to the bubble’s rapid expansion, observed in images B1-C1 in Fig. 8(a) and B1-D1 in Fig. 8(b). When a vapor bubble is formed the laser radiation is transmitted through the expanding bubble and is absorbed at its opposite side.3,36 Therefore, the bubble is not only vapor-driven; it is also pulse-driven, since its energy increases due to the absorption of laser energy at the bubble’s wall. Here, each individual laser spike contributes to the bubble growth through the evaporation, which happens immediately after the irradiated water reaches the boiling point.22

During the expansion, the bubble passes over the equilibrium state. Thus, at its maximum volume [e.g., see images D1 in Fig. 8(a) and A2 in Fig. 8(b)], the internal pressure is lower than the pressure in the surrounding liquid. This difference in pressures forces the bubble to collapse [see images B2-D2 in Fig. 3(a) and B2-C2 in Fig. 8(b)]. During the collapse, a portion of the energy stored in the bubble is converted into acoustic energy. This results in the emission of acoustic transients, shock waves (SW), labeled in images D2 in Fig. 8(a) and C2 in Fig. 8(b). The bubble’s oscillations are repeated several times and a new shock wave is emitted after each bubble’s rebound [e.g., see images B3 in Fig. 8(a) and 8(b)]. The pressures of the shock waves induced during the bubble’s expansion and rebounds were measured by different authors22,32,35,37 using a hydrophone. Their results show that in the case of free-running-mode pulses, the high-pressure transients exist only during rebounds. This is consistent with the shadowgraphic observations in our experiments.

The horizontal arrows in Fig. 8 indicate the FT’s end. In case of the flat FT [Fig. 8(a)], the laser radiation ablates a cylindrical channel into the water. As is visible in images C1-A2 in Fig. 8(a), this bubble not only grows in the downward direction due to the pulse vaporization, but also expands up along the fiber. Moreover, the top part of this channel-like bubble is wider compared to the rest of the bubble. The reason for this lies in the higher pressure at the beginning of the bubble’s formation.

During the collapse, the bubble is divided into two parts [see image B2 in Fig. 8(a)]. For the laser parameters used in Fig. 8, only the bottom part experiences several rebounds [images B2-C3 in Fig. 8(a)]. In contrast, the top part, which collapses near the FT, already decays during the first collapse.

On the other hand, when the pulse is delivered through the conical FT, we observed a spherical bubble [Fig. 8(b)]. Here, the center of the bubble within the first expansion and the first collapse coincides with the point of the conical FT [images B1-B2 in Fig. 8(b)]. However, the shock wave generated after the first bubble’s collapse [see image C2 in Fig. 8(b)] has its center below the FT’s point. This indicates that the bubble’s center is moved before the second expansion of the bubble, since the acoustic transient is emitted at the end of the bubble’s collapse.44


Volume as a function of time

Applying the image processing described in the Sec. 2.2, the bubble’s volume VB can be obtained from shadowgraphs using Eqs. (1)–(3). A typical bubble’s volume as a function of time for two different pulse parameters and both FT geometries is presented in Fig. 9. Here, we show the results for the same pulse energy EL=20mJ, but different pulse durations: tL=120μs (the squares) and tL=60μs (the circles). The vertical lines mark the end of the excitation pulse. The pulse power as a function of time for the data in Fig. 9 is shown in Fig. 4. The effect of pulse shortening is clearly visible. If the pulse duration is reduced while keeping constant pulse energy, the bubble’s maximum volume increases. It can be also concluded from Fig. 9 that the bubbles induced through the conical FT (the full data marks) reach a larger maximum volume in comparison to those formed through the flat FT (the open data marks).

Fig. 9

Bubble’s volume as a function of time. Measurements are shown for both FT geometries and two pulse durations (tL=60μs and tL=120μs) at an unchanged pulse energy of 20 mJ on the exit of the FT.



Influence of the Pulse Duration on the Bubble’s Expansion

The expansion of a spherical bubble in an infinitely large, non-viscous, and incompressible liquid can be roughly described by the Rayleigh–Plesset equation:45


where the dot stands for a time derivative, r is the bubble’s radius as a function of time, ρ is the liquid’s density, and pr(t)p(t)=Δp(t) the pressure difference between that at the bubble’s wall pr(t), and the pressure of the surrounding liquid p(t). If the pressure difference Δp is constant, Eq. (10), which is a general equation of motion for a spherical bubble, simplifies into the well-known Rayleigh equation:33


Here, rmax stands for the maximum radius of the bubble.

The spherical bubbles are induced through the conical FT. For this reason, we used these measurements to analyze the influence of pulse duration on the expansion of vapor bubbles. In this case, the bubble’s radius r as a function of time can be calculated from the measured volume VB(t) as:



For a better interpretation of the experimental results, Eq. (11) should be rewritten in a dimensionless form using the following normalization: ξ=r/rmax and τ=t/tmax.39 Here, tmax stands for the time when the bubble reaches its maximum radius rmax.

The integration of Eq. (11) yields the following expression in a dimensionless form:


where B and Iξ3 denote the Beta and the incomplete Beta function,46 respectively.

When cavitation bubbles in an infinite liquid are induced by short (i.e., nanosecond) laser pulses, their dynamics conforms to the Rayleigh model and can be unified in the dimensionless graph ξ(τ).39,47 This confirms that the assumption of a constant pressure difference Δp during the bubble’s expansion is reasonable in the case of short pulses. In contrast, when the pulse duration is comparable to the bubble-expansion time, the evaporation takes place in a liquid shell due to the pulse absorption during the bubble’s growth. We analyzed the influence of the pulse duration on the bubble’s dynamics by a comparison of the measured results with the predictions of the Rayleigh model [Eq. (13)].

Figure 10 shows the dimensionless radius as a function of the dimensionless time for the spherical bubbles induced by the pulses with durations of 60 μs (the crosses), 145 μs (the dots), and 285 μs (the pluses). Here, the measured data are normalized with the maximum measured radii rmax and the measured times tmax. The pulse durations τL in a dimensionless form are also presented and the solid curve shows the Rayleigh model [Eq. (13)]. It can be concluded that the measured bubble’s expansion approaches the Rayleigh model when the pulse duration decreases (i.e., τL0). In contrast, for longer pulses, the “pulse-driven effect” is clearly visible. In this case, the experimental results do not conform to the Rayleigh model, due to the evaporation induced by the light absorption. Therefore, the assumption of a constant pressure difference Δp during the bubble’s expansion is no longer valid.

Fig. 10

Dimensionless radii as a function of the dimensionless time for spherical bubbles induced through the conical FT. τL labels normalized pulse durations.



Optodynamic Energy-Conversion Efficiency

In our case, the shock wave during the bubble’s formation is not present. This finding is important due to two reasons: (i) it shows that at the wavelength of 2.94 μm the shock waves cannot be the main mechanism for cleaning the root canal, as suggested in Refs. 15, 16, 17; and (ii) it means that the total mechanical energy of the liquid medium equals the mechanical energy of the first bubble expansion.

When the bubble reaches its maximum volume, its radial velocity vanishes. The surface tension for millimeter-sized bubbles can be neglected, so a bubble at its maximum radius possesses only a potential energy.48 Thus, the total energy of the bubble EB can be obtained as the product of its maximum volume Vmax and the difference between the hydrostatic and internal (vapor) pressure.49,50 Since the vapor pressure at Vmax is much smaller than the hydrostatic pressure,51 it can be neglected. Thus, the bubble’s energy can be estimated using the following equation:39


where, p=105Pa equals the hydrostatic pressure of the water. The OD energy-conversion efficiency η is therefore defined as:



We determined Vmax from measurements of the volume as a function of time (e.g., see Fig. 9) as follows. First, we fit the quartic function to the measured data corresponding to volumes larger than Vmax/2. Then we determined the time tmax as the time related to the maximum of the fitted function. Finally, Vmax is calculated as an average of the volumes measured within the time interval [0.9,1.1]×tmax.

Using Eq. (15), we have analyzed how (i) the pulse energy, (ii) the pulse duration, and (iii) the FT geometry influence the OD energy-conversion efficiency. Figure 11 shows the measured efficiency as a function of the pulse duration, where the pulse energy (EL=20mJ) was kept constant. On the other hand, Fig. 12 presents the measured efficiency as a function of the pulse energy, where the pulse duration (tL=145μs) was constant during the measurements. The results are presented for the conical (the full circles) and the flat (the open circles) FTs. The error bars correspond to the standard deviation calculated from the volume averaging.

Fig. 11

OD energy-conversion efficiency as a function of pulse duration at a constant pulse energy (EL=20mJ).


Fig. 12

OD energy-conversion efficiency as a function of pulse energy at a constant pulse duration (tL=145μs).


It can be concluded from the results presented in Figs. 11 and 12 that the OD energy-conversion efficiency (i) decreases with an increasing pulse duration, (ii) increases by increasing the pulse energy, and (iii) is significantly larger when a conical FT is used rather than a flat FT. These findings have different meanings for various biomedical applications.

The mechanical effects in water should be minimized when a tissue is ablated in a non-contact mode under water (e.g., see Ref. 22). Here, the water serves as a cooler. Therefore, the longer pulses and a flat FT with a channel-like bubble and lower efficiency are more appropriate for this case.

On the other hand, in the medical treatments such as cleaning root canals, the maximization of the mechanical effects is especially important.16 Causing the same mechanical effects, a higher efficiency allows less collateral heating, which increases the temperature of the surrounding medium and results in collateral tissue damage and increasing pain. In this case, it is thus more convenient to use conical FTs instead of the flat ones and to keep the duration of the laser pulses as short as possible. However, when higher pulse powers are necessary (i.e., the pulse energy exceeds 100 mJ within a pulse duration of less than 100 μs), the flat FT should be used, because at the point of the conical FT the intensity is increased. Thus, a conical termination experiences damage at lower pulse powers than the flat one.



We have studied the influence of the FT geometry and the parameters of the free-running erbium pulse on the vapor-bubble formation and dynamics. The experiments were performed using shadow photography. Two different FT terminations, a flat and a conical, were used. The shadowgraphs showed that the bubble’s shape depends on FT geometry. A spherical bubble develops when the conical FT is used, while a channel-like bubble is obtained with the flat one. In both cases, bubbles experience several oscillations and a shock wave is emitted during each collapse.

During the first bubble expansion, the shock wave is not present. Therefore, the total mechanical energy of the liquid equals the energy of the first bubble’s expansion. This enables us to define the OD energy-conversion efficiency as the ratio between the bubbles’ energy and the energy of the laser pulse. Our results show that the efficiency increases (i) by increasing the pulse energy and (ii) decreasing the pulse duration. Moreover, it is larger when a conical FT is used rather than a flat one. The presented results are important for the optimization of biomedical processes, where the collateral heating of the surrounding medium should be reduced.


We would like to thank Fotona, d.d. for supplying their commercially available laser system for laser cleaning root canals in clinical dentistry applications. Their system, including the Er:YAG laser (LightWalker), handpiece, and interchangeable fiber tips was used as an excitation source in our experiments.


1. R. Brinkmannet al., “Analysis of cavitation dynamics during pulsed laser tissue ablation by optical on-line monitoring,” IEEE J. Sel. Top. Quant. 2(4), 826–835 (1996).IJSQEN1077-260X http://dx.doi.org/10.1109/2944.577305 Google Scholar

2. H. Pratistoet al., “Combination of fiber-guided pulsed erbium and holmium laser radiation for tissue ablation under water,” Appl. Opt. 35(19), 3328–3337 (1996).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.35.003328 Google Scholar

3. T. Lüet al., “Cavitation effect of holmium laser pulse applied to ablation of hard tissue underwater,” J. Biomed. Opt. 15(4), 048002 (2010).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.3470092 Google Scholar

4. T. PerhavecJ. Diaci, “A novel double-exposure shadowgraph method for observation of optodynamic shock waves using fiber-optic illumination,” Strojniški Vestnik—J. Mech. Eng. 56(7–8), 477–482 (2010). Google Scholar

5. A. BizjakK. NemešJ. Možina, “Rotating-mirror Q-switched Er:YAG laser for optodynamic studies,” Strojniški Vestnik—J. Mech. Eng. 57(1), 3–10 (2011). http://dx.doi.org/10.5545/sv-jme.2010.120 Google Scholar

6. G. M. HaleM. R. Querry, “Optical-constants of water in 200-nm to 200-μm wavelength region,” Appl. Opt. 12(3), 555–563 (1973).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.12.000555 Google Scholar

7. K. F. Chanet al., “Erbium:YAG laser lithotripsy mechanism,” J. Urol. 168(2), 436–441 (2002).JOURAA0022-5347 http://dx.doi.org/10.1016/S0022-5347(05)64653-8 Google Scholar

8. H. PetersenM. MrochenT. Seiler, “Comparison of erbium:yttrium-aluminum-garnet-laser vitrectomy and mechanical vitrectomy,” Ophthalmology 107(7), 1389–1392 (2000).OPANEW0743-751X http://dx.doi.org/10.1016/S0161-6420(00)00123-8 Google Scholar

9. M. Mrochenet al., “Erbium:yttrium-aluminum-garnet laser induced vapor bubbles as a function of the quartz fiber tip geometry,” J. Biomed. Opt. 6(3), 344–350 (2001).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.1381052 Google Scholar

10. A. Mehlet al., “Bactericidal effects of 2.94 μm Er:YAG-laser radiation in dental root canals,” J. Endodont. 25(7), 490–493 (1999).0099-2399 http://dx.doi.org/10.1016/S0099-2399(99)80288-X Google Scholar

11. H. TachinamiI. Katsuumi, “Removal of root canal filling materials using Er:YAG laser irradiation,” Dent. Mater. J. 29(3), 246–252 (2010).DMJOD5 http://dx.doi.org/10.4012/dmj.2008-079 Google Scholar

12. H. MatsumotoY. YoshimineA. Akamine, “Visualization of irrigant flow and cavitation induced by Er:YAG laser within a root canal model,” J. Endodont. 37(6), 839–843 (2011).0099-2399 http://dx.doi.org/10.1016/j.joen.2011.02.035 Google Scholar

13. R. GeorgeL. J. Walsh, “Performance assessment of novel side firing safe tips for endodontic applications,” J. Biomed. Opt. 16(4), 048004 (2011).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.3563637 Google Scholar

14. R. GeorgeL. J. Walsh, “Performance assessment of novel side firing flexible optical fibers for dental applications,” Lasers Surg. Med. 41(3), 214–221 (2009).LSMEDI0196-8092 http://dx.doi.org/10.1002/lsm.v41:3 Google Scholar

15. D. L. Gloveret al., “Laser based enhanced generation of photoacoustic pressure waves in dental and medical treatments and procedures,” U.S. Patent No. 7959441 B2 (2011). Google Scholar

16. E. DiVitoO. A. PetersG. Olivi, “Effectiveness of the Erbium:Yag laser and new design radial and stripped tips in removing the smear layer after root canal instrumentation,” Laser. Med. Sci. 27(2), 273–280 (2012).LMSCEZ1435-604X http://dx.doi.org/10.1007/s10103-010-0858-x Google Scholar

17. G. OliviE. DiVito, “Photoacoustic endodontics using PIPS™: experimental background and clinical protocol,” J. Laser Health Acad. 2012(1), 22–25 (2012). Google Scholar

18. C. P. LinD. SternC. A. Puliafito, “High-speed photography of Er-Yag laser ablation in fluid - implication for laser vitreous surgery,” Invest. Ophth. Vis. Sci. 31(12), 2546–2550 (1990).IOVSDA0146-0404 Google Scholar

19. E. D. Jansenet al., “Partial vaporization model for pulsed mid-infrared laser ablation of water,” J. Appl. Phys. 78(1), 564–571 (1995).JAPIAU0021-8979 http://dx.doi.org/10.1063/1.360642 Google Scholar

20. A. VogelV. Venugopalan, “Mechanisms of pulsed laser ablation of biological tissues,” Chem. Rev. 103(2), 577–644 (2003).CHREAY0009-2665 http://dx.doi.org/10.1021/cr010379n Google Scholar

21. H. LoertscherW. Q. ShiW. S. Grundfest, “Tissue ablation through water with erbium:YAG lasers,” IEEE T. Biomed. Eng. 39(1), 86–88 (1992).IEBEAX0018-9294 http://dx.doi.org/10.1109/10.108131 Google Scholar

22. M. Ithet al., “Dynamics of laser-induced channel formation in water and influence of pulse duration on the ablation of biotissue under water with pulsed erbium-laser radiation,” Appl. Phys. B 59(6), 621–629 (1994).APBOEM0946-2171 http://dx.doi.org/10.1007/BF01081183 Google Scholar

23. E. A. Brujanet al., “Dynamics of laser-induced cavitation bubbles near an elastic boundary,” J. Fluid Mech. 433, 251–281 (2001).JFLSA70022-1120 Google Scholar

24. P. Gregorčičet al., “Measurements of cavitation bubble dynamics based on a beam-deflection probe,” Appl. Phys. A 93(4), 901–905 (2008).APSFDB0721-7250 http://dx.doi.org/10.1007/s00339-008-4751-4 Google Scholar

25. R. Evanset al., “Pump-probe imaging of nanosecond laser-induced bubbles in agar gel,” Opt. Express 16(10), 7481–7492 (2008).OPEXFF1094-4087 http://dx.doi.org/10.1364/OE.16.007481 Google Scholar

26. A. B. Karpiouket al., “Quantitative ultrasound method to detect and monitor laser-induced cavitation bubbles,” J. Biomed. Opt. 13(3), 034011 (2008).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.2937478 Google Scholar

27. W. LauterbornT. Kurz, “Physics of bubble oscillations,” Rep. Prog. Phys. 73(10), 106501 (2010).RPPHAG0034-4885 http://dx.doi.org/10.1088/0034-4885/73/10/106501 Google Scholar

28. A. Vogelet al., “Intraocular Nd-Yag laser-surgery—light tissue interaction, damage range, and reduction of collateral effects,” IEEE Quant. Elect. 26(12), 2240–2260 (1990).IEJQA70018-9197 http://dx.doi.org/10.1109/3.64361 Google Scholar

29. F. G. Perez-GutierrezS. Camacho-LopezG. Aguilar, “Time-resolved study of the mechanical response of tissue phantoms to nanosecond laser pulses,” J. Biomed. Opt. 16(11), 115001 (2011).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.3644380 Google Scholar

30. P. GregorčičR. PetkovšekJ. Možina, “Investigation of a cavitation bubble between a rigid boundary and a free surface,” J. Appl. Phys. 102(9), 094904 (2007).JAPIAU0021-8979 http://dx.doi.org/10.1063/1.2805645 Google Scholar

31. T. G. van Leeuwenet al., “Noncontact tissue ablation by holmium:YSGG laser-pulses in blood,” Laser. Surg. Med. 11(1), 26–34 (1991).LSMEDI0196-8092 http://dx.doi.org/10.1002/(ISSN)1096-9101 Google Scholar

32. T. AsshauerK. RinkG. Delacretaz, “Acoustic transient generation by holmium-laser-induced cavitation bubbles,” J. Appl. Phys. 76(9), 5007–5013 (1994).JAPIAU0021-8979 http://dx.doi.org/10.1063/1.357212 Google Scholar

33. L. Rayleigh, “On the pressure developed in a liquid during the collapse of a spherical cavity,” Philos. Mag. 34(200), 94–98 (1917).PHMAA40031-8086 http://dx.doi.org/10.1080/14786440808635681 Google Scholar

34. E. A. BrujanA. Vogel, “Stress wave emission and cavitation bubble dynamics by nanosecond optical breakdown in a tissue phantom,” J. Fluid Mech. 558, 281–308 (2006).JFLSA70022-1120 http://dx.doi.org/10.1017/S0022112006000115 Google Scholar

35. E. D. Jansenet al., “Effect of pulse duration on bubble formation and laser-induced pressure waves during holmium laser ablation,” Laser. Surg. Med. 18(3), 278–293 (1996).LSMEDI0196-8092 http://dx.doi.org/10.1002/(ISSN)1096-9101 Google Scholar

36. T. Asshaueret al., “Pulsed holmium laser ablation of tissue phantoms: correlation between bubble formation and acoustic transients,” Appl. Phys. B 65(4–5), 647–657 (1997).APBOEM0946-2171 http://dx.doi.org/10.1007/s003400050327 Google Scholar

37. M. Frenzet al., “Starting mechanisms and dynamics of bubble formation induced by a Ho:Yttrium aluminum garnet laser in water,” J. Appl. Phys. 84(11), 5905–5912 (1998).JAPIAU0021-8979 http://dx.doi.org/10.1063/1.368906 Google Scholar

38. R. PetkovšekP. GregorčičJ. Možina, “A beam-deflection probe as a method for optodynamic measurements of cavitation bubble oscillations,” Meas. Sci. Technol. 18(9), 2972–2978 (2007).MSTCEP0957-0233 http://dx.doi.org/10.1088/0957-0233/18/9/030 Google Scholar

39. R. PetkovšekP. Gregorčič, “A laser probe measurement of cavitation bubble dynamics improved by shock wave detection and compared to shadow photography,” J. Appl. Phys. 102(4), 044909 (2007).JAPIAU0021-8979 http://dx.doi.org/10.1063/1.2774000 Google Scholar

40. G. S. Settles, Schlieren and Shadowgraph Techniques, Springer-Verlag, Berlin (2001). Google Scholar

41. S. L. Jacques, “Role of tissue optics and pulse duration on tissue effects during high-power laser irradiation,” Appl. Opt. 32(13), 2447–2454 (1993).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.32.002447 Google Scholar

42. O. Fohnet al., "Side-firing fiber device for underwater tissue ablation with Ho:YAG and Er:YAG laser radiation," J. Biomed. Opt. 3(1), 112–122 (1998).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.429868 Google Scholar

43. R. O. Esenalievet al., “Studies of acoustical and shock-waves in the pulsed-laser ablation of biotissue,” Laser. Surg. Med. 13(4), 470–484 (1993).LSMEDI0196-8092 http://dx.doi.org/10.1002/(ISSN)1096-9101 Google Scholar

44. I. Akhatovet al., “Collapse and rebound of a laser-induced cavitation bubble,” Phys. Fluids 13(10), 2805–2819 (2001).PHFLE61070-6631 http://dx.doi.org/10.1063/1.1401810 Google Scholar

45. M. S. Plesset, “The dynamics of cavitation bubbles,” J. Appl. Mech. 16, 228–231 (1949).JAMCAV 0021-8936 Google Scholar

46. M. AbramowitzI. A. Stegun, Handbook of Mathematical Functions, Dover Publications, New York (1972). Google Scholar

47. S. W. Gonget al., “Scaling law for bubbles induced by different external sources: Theoretical and experimental study,” Phys. Rev. E 81(5), 056317 (2010).PLEEE81063-651X http://dx.doi.org/10.1103/PhysRevE.81.056317 Google Scholar

48. W. Kreideret al., “Observations of the collapses and rebounds of millimeter-sized lithotripsy bubbles,” J. Acoust. Soc. Am. 130(5), 3531–3540 (2011).JASMAN0001-4966 http://dx.doi.org/10.1121/1.3626157 Google Scholar

49. A. VogelW. LauterbornR. Timm, “Optical and acoustic investigations of the dynamics of laser-produced cavitation bubbles near a solid boundary,” J. Fluid Mech. 206, 299–338 (1989).JFLSA70022-1120 http://dx.doi.org/10.1017/S0022112089002314 Google Scholar

50. S. BuogoG. B. Cannelli, “Implosion of an underwater spark-generated bubble and acoustic energy evaluation using the Rayleigh model,” J. Acoust. Soc. Am. 111(6), 2594–2600 (2002).JASMAN0001-4966 http://dx.doi.org/10.1121/1.1476919 Google Scholar

51. A. Vogelet al., “Energy balance of optical breakdown in water at nanosecond to femtosecond time scales,” Appl. Phys. B 68(2), 271–280 (1999).APBOEM0946-2171 http://dx.doi.org/10.1007/s003400050617 Google Scholar

© 2012 Society of Photo-Optical Instrumentation Engineers (SPIE)
Peter Gregorcic, Matija Jezersek, Janez Mozina, "Optodynamic energy-conversion efficiency during an Er:YAG-laser-pulse delivery into a liquid through different fiber-tip geometries," Journal of Biomedical Optics 17(7), 075006 (9 July 2012). https://doi.org/10.1117/1.JBO.17.7.075006

Back to Top