Invited p aper Mechanisms of femtosecond laser nanosurgery of ...
Invited p aper Mechanisms of femtosecond laser nanosurgery of ...
Invited p aper Mechanisms of femtosecond laser nanosurgery of ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
VOGEL et al. <strong>Mechanisms</strong> <strong>of</strong> <strong>femtosecond</strong> <strong>laser</strong> <strong>nanosurgery</strong> <strong>of</strong> cells and tissues 1035lated using the measured spot size. However, in various studies,the quality <strong>of</strong> the <strong>laser</strong> beam and focusing optics were farfrom optimal and, nevertheless, the diffraction-limited spotsize has been used for the calculation <strong>of</strong> I th . In the presence <strong>of</strong>aberrations, this leads to erroneously large threshold values.Even when the spot size is measured, aberrations will lead t<strong>of</strong>alse results because breakdown occurs in hot spots within thefocal region that are smaller than the measured spot size. Thethreshold values determined in this way are, therefore, smallerthan the actual values [184]. Many breakdown experimentswere performed by focusing a <strong>laser</strong> beam into a glass cellby means <strong>of</strong> external focusing optics. Even when this opticsis <strong>of</strong> high quality, the glass cell itself will introduce sphericalaberrations. This source <strong>of</strong> error can be avoided by usingwater-immersion microscope objectives built into the wall <strong>of</strong>the cell [3].For nanosecond pulses, longitudinal mode beating in the<strong>laser</strong> resonator will usually produce intensity spikes <strong>of</strong> picosecondduration with amplitudes up to 4–5 times as large asthe averaged intensity. Such statistical fluctuations <strong>of</strong> the peakirradiance are a major cause <strong>of</strong> the wide statistical distribution<strong>of</strong> optical breakdown thresholds and substantially reducetheir absolute values [185, 186]. Therefore, single-frequency<strong>laser</strong>s providing a temporally Gaussian pulse shape shouldbe used for threshold measurements. However, to date allmeasurements <strong>of</strong> nanosecond optical breakdown thresholdsin water have been performed with longitudinally multimode<strong>laser</strong> beams.To avoid a change <strong>of</strong> the focal spot size by nonlinearbeam propagation, sufficiently large focusing angles must beused [5, 15, 187]. While this is easy to achieve for nanosecondpulses, a numerical <strong>aper</strong>ture NA ≥ 0.9 is required fora pulse duration <strong>of</strong> 100 fs to be sure to exclude a diminution<strong>of</strong> the spot size by self-focusing and the corresponding apparentreduction <strong>of</strong> the breakdown threshold [5]. To date, nothreshold data obtained with such large NAs are available forultra-short <strong>laser</strong>-induced breakdown in water. We must concludethat the present experimental data basis is not reliableenough to perform a thorough comparison with theoreticalpredictions. Future measurements with aberration-free temporallyGaussian <strong>laser</strong> pulses focused at large NA will have t<strong>of</strong>ill this gap.In addition, a better adjustment <strong>of</strong> the numerical breakdowncriterion to the experimental criterion <strong>of</strong> bubble formationis needed to enable a meaningful comparison <strong>of</strong> experimentaldata with model predictions. While bubble formationrequires an approximately constant energy density within thefocal volume for all <strong>laser</strong> pulse durations and wavelengths,the energy density associated with a fixed value <strong>of</strong> the freeelectrondensity, such as ϱ cr = 10 21 cm −3 , varies considerablywith pulse duration. Thus, the assumption <strong>of</strong> a constant freeelectrondensity as breakdown criterion is quite arbitrary, especiallyfor cases where the threshold is smooth, i.e. whereϱ max increases continuously with irradiance. In these casesit seems more reasonable to relate the critical free-electrondensity to the energy density within the medium that leadsto bubble formation. Equation (17) provides the required linkbetween electron and energy densities, and an analysis <strong>of</strong>ϱ max (I/I rate ) curves (Figs. 3 and 5) then yields the correspondingthreshold irradiance.The new threshold definition for the numerical model willyield a better match between the experimentally observedand numerically predicted wavelength dependence <strong>of</strong> breakdownthresholds than the values in Table 1. Experimentally,the thresholds for ns and ps pulses decrease with decreasingwavelength [3, 15], but values calculated with a constant ϱ crshow the opposite trend. Previous calculations revealed a relativelysmooth threshold for ns and ps breakdown with visibleand UV wavelengths, in contrast to the sharp threshold forIR wavelengths shown in Fig. 3c [91]. Calculated thresholdvalues based on the new definition will for IR wavelengthsthus be the same as those obtained with ϱ cr = 10 21 cm −3 ,butthey will be lower for visible and UV wavelengths, in accordwith experimental observations.6.4 Cavitation bubble dynamicsFigure 18 shows a two-dimensional representation<strong>of</strong> the evolution <strong>of</strong> the thermoelastic stress wave and <strong>of</strong> theregion in which the kinetic spinodal is surpassed (bubble nucleus)for a peak temperature <strong>of</strong> 200 ◦ C, slightly above thethreshold for bubble formation. The subsequent bubble dynamicscalculated for cases 1 and 2 described in Sect. 6.1 ispresented in Fig. 19. It turned out that for a vapor pressurethat changes adiabatically with bubble size (case 2 in Fig. 19b)the bubble motion is very similar to the case where the vaporpressure is not taken into account at all and only the negativepressurepulse drives the bubble expansion (not shown). Figure20 shows the maximum bubble radii for case 1 and case2 as a function <strong>of</strong> the maximum temperature achieved in thecenter <strong>of</strong> the focal volume, together with the radius <strong>of</strong> the nucleus,R 0 . Close to threshold, the actual bubble radius willprobably resemble case 1 in which the bubble nucleus is assumedto be much smaller than the focal region. For energieswell above threshold, the predictions obtained for case 2 becomemore realistic.The most prominent feature <strong>of</strong> the transient bubbles producedclose to the threshold <strong>of</strong> <strong>femtosecond</strong> optical breakdownis their small size and short lifetime. The bubble radiusamounts to only about 200 nm in water, and will be evensmaller in a visco-elastic medium such as the cytoplasm. Thismakes a dissection mechanism associated with bubble formationcompatible with intracellular <strong>nanosurgery</strong>, in contrastto nanosecond optical breakdown (τ L = 6ns,λ= 1064 nm)where the smallest bubble radius in water observed for NA =0.9 was R max = 45 µm [3]. The small bubble size correspondsto a small energy,E B = 4π 3 (p ∞ − p v )R 3 max , (21)<strong>of</strong> the expanded bubble. For the case presented in Fig. 19a, itamounts to 5.3 × 10 −15 J (5.3 femtojoules). The smallness <strong>of</strong>the bubble energy is largely due to the small energy content <strong>of</strong>the stress transient creating the bubble. This energy is equal tothe energy stored in thermoelastic compression <strong>of</strong> the heatedfluid volume [156],E TE = 12ϱ 0 c 2 0∫p 2 0 dV , (22)