ISSN 00213640, JETP Letters, 2015, Vol. 101, No. 6, pp. 394–397. © Pleiades Publishing, Inc., 2015. Original Russian Text © D.A. Zayarny, A.A. Ionin, S.I. Kudryashov, S.V. Makarov, A.A. Rudenko, S.G. Bezhanov, S.A. Uryupin, A.P. Kanavin, V.I. Emel’yanov, S.V. Alferov, S.N. Khonina, S.V. Karpeev, A.A. Kuchmizhak, O.B. Vitrik, Yu.N. Kulchin, 2015, published in Pis’ma v Zhurnal Eksperimental’noi i Teoreticheskoi Fiziki, 2015, Vol. 101, No. 6, pp. 428–432.
Nanoscale Boiling during SingleShot Femtosecond Laser Ablation of Thin Gold Films¶ D. A. Zayarnya, A. A. Ionina, S. I. Kudryashova, b, *, S. V. Makarova, c, A. A. Rudenkoa, S. G. Bezhanova, b, S. A. Uryupina, b, A. P. Kanavina, b, V. I. Emel’yanovd, S. V. Alferove, S. N. Khoninae, S. V. Karpeeve, A. A. Kuchmizhakf, O. B. Vitrik f, g, and Yu. N. Kulchin f, g a
Lebedev Physical Institute, Russian Academy of Sciences, Leninskii pr. 53, Moscow, 119991 Russia *email:
[email protected] b National Research Nuclear University MEPhI (Moscow Engineering Physics Institute), Kashirskoe sh. 31, Moscow, 115409 Russia c ITMO University, St. Petersburg, 197101 Russia d Moscow State University, Moscow, 119899 Russia e Image Processing Systems Institute, Russian Academy of Sciences, Samara, 443001 Russia f Institute for Automation and Control Processes, Far Eastern Branch, Russian Academy of Science, Vladivostok, 690041 Russia g Far Eastern Federal University, Vladivostok, 690041 Russia Received January 26, 2015
A nanoscale chaotic relief structure appears as a result of subthreshold singleshot femtosecond laser ablation of gold films in the regimes of fabrication of microbumps and nanospikes, but only for a relatively thick film. The observed nanoablation tendency versus film thickness makes it possible to suppose the existence of a sub surface temperature maximum in thicker gold films and its absence within thinner film, which results from competing evaporative cooling and electronic heat conduction, as demonstrated by numerical simulations of the thermal dynamics. DOI: 10.1134/S0021364015060132
1. Despite numerous preceding observations of various surface mesostructures—nanocrowns, microbumps, nanojets, and nanoholes—in thin films and on surfaces of bulk materials during their single and multishot nanoscale ablation by nanosecond [1– 5] and femtosecond [6–8] laser pulses, their formation mechanisms are not well understood yet. A rather good approximation in understanding the mecha nisms of nanoscale ablation is the description of mac roscopic femtosecond (fs) laser ablation provided in a number of recent studies [9–15]. In particular, it was demonstrated that spallative femtosecond laser abla tion of thin films and bulk materials is initiated by delayed subsurface boiling of a molten layer [13, 14], rather than thermoelastic stresses, with picosecond evaporative cooling of the layer suppressing its phase explosion [15]. As a result, nanoscale spallation in thin films proceeds in two stages: (1) irreversible formation of a microbump apparently owing to high vapor pres sure in a vapor cavity between the film and its substrate [2, 10, 12] and (2) complete film liftoff during hydro dynamic expulsion of a molten nanojet [11]. ¶The article was translated by the authors.
In the present work, we report results of a compar ative study of singleshot fs laser ablation of thinner and thicker gold films, indicating similarity of their relief evolution versus fs laser fluence, with exception of a subthreshold nanoroughness on the thicker film. 2. In our studies, laser irradiation of fresh spots of samples was performed by single femtosecond, second harmonic (SH) pulses of a Yb3+doped fiber laser (Satsuma, Amplitude Systemes) [16]: SH wavelength is 515 nm, FWHM duration is 200 fs, maximum pulse energy is 4 μJ in TEM00 mode, repetition rate is 0– 2 MHz. Laser radiation was focused onto the sample surface in air by a microscope objective with numerical aperture NA = 0.25 into a circular spot with the radius R1/e ≈ 0.45 μm (Fig. 1). The sample was arranged on a threedimensional motorized translation stage with minimum step of 150 nm, moved from pulse to pulse. The resulting ablative relief was visualized using a JEOL 7001F scanning electron microscope (SEM). As samples, we used thin films of gold (80%) and pal ladium (20%) alloy with their thicknesses of h ≈ 60 and ≈200 nm, deposited in an argon atmosphere on dielec tric (silica glass) substrates by means of a magnetron sputtering system (SC7620, Quorum Technologies).
394
NANOSCALE BOILING DURING SINGLESHOT FEMTOSECOND LASER
395
Fig. 1. (Color online) Experimental setup for nanoscale laser ablation of material surfaces: (RA) reflection attenuator, (AC) auto correlator, (EM) pyroelectric energy meter, (WL) whitelight halogen illumination lamp, (BS) beam splitter, (CCD) CCD cam era, and (PC) laptop for control of the laser and 3D motorized translation stage.
3. Singleshot fs laser exposure of the thinner film at the peak laser fluence F ≈ 0.5 J/cm2 resulted in its partial melting and resolidification, becoming evident in the form of significantly larger recrystallized metal lic nanograins (Fig. 2a). In contrast, lowfluence fs laser modification of the thicker film at F ≈ 1.5 J/cm2 occurred as surface nanofoaming with the nanorough ness scale increasing toward the spot center (Fig. 2d). The origin of such nanoscale surface foaming becomes clear at higher laser fluences: for the thinner film at F ≈ 0.7 J/cm2 a central microbump appears, besides the film melting (Fig. 2b). Similarly, for the thicker film at F ≈ 1.8 J/cm2, such a central microbump also emerges on the ablated spot, but the nanofoam structure is now smoothened on its remolten top (Fig. 2e). Finally, at a further increase in fluence, a central nanospike (“fro zen nanojet”) appears for both the thinner (F ≈ 0.9 J/cm2) and thicker (F ≈ 2.5 J/cm2) films (Figs. 2c, 2f), with one or a few nanoparticles atop. In the latter regime, the foamlike structure on the thicker film is displaced to the crater periphery, appearing as a regu lar nanocrown with one or two circular rows of tips (Fig. 2f). The observed evolution of the surface topog raphy of the thicker film—nanofoam, microbump, nanospike—demonstrates obvious consistency not only in the peak fluence domain but, apparently, also in terms of the film temperature increasing to near critical magnitudes to initiate explosive growth of a nanojet [11, 17]. Importantly, the characteristic flu JETP LETTERS
Vol. 101
No. 6
2015
ence values for the appearance of the different nanofeatures in the thicker films are well scalable to the thinner film, accounting for the ratio of their thicknesses. In our opinion, these facts indicate that, for the thicker film, subsurface boiling takes place, which is, among other physical and chemical factors, the main driving process not only for microbump and nanojet formation on surfaces of thin films [2, 10, 12], but also for surface spallation in bulk materials [13, 14]. One possible reason for subsurface boiling can be the sub surface temperature maximum [5, 18–23], emerging on a subnanosecond timescale within a thick film owing to a balance of heat conduction and intense evaporative cooling, as experimentally demonstrated in [15]. Otherwise, in the case of thin film, heat con duction is fast and evaporative cooling shifts the corre sponding temperature maximum toward the unirradi ated film edge. In order to calculate temperature profiles in these films, numerical simulations were undertaken for temperatures of electrons and lattice for the condi tions of our experiments. The equation for electron temperature accounted for a heating source term cor responding to optical absorption of the incident inho mogeneous electromagnetic field in the films, elec tronic heat conduction term, and electron–lattice coupling term. The latter term was also included in the equation for lattice temperature. The corresponding
396
ZAYARNY et al.
Fig. 3. Temperature profile in the 200nm thick gold film 13 ps after single 515nm, 200fs, 10TW/cm2 pump pulse. Inset: Temperature profile in the 60nmthick gold film 11 ps after single 515nm, 200fs, 7TW/cm2 pump pulse. The displacement of each curve relative to the beginning of the x axis indicates the thickness of the vapor ized layer at the corresponding time instant.
Fig. 2. Obliqueincidence (≈45°) SEM images of 60nm and 200nm gold films ablated by single fs laser pulses at F ≈ (a) 0.5, (b) 0.7, (c) 0.9, (d) 1.5, (e) 1.8, and (f) 2.5 J/cm2. The remolten central zone is shown by the dashed circle in image (a). The scale bars are different for different images.
temperature dependences of heat capacities, elec tron–lattice coupling, and effective collision frequen cies were taken from [24–26], which were previously tested in simulations of gold film absorbance [27]. Moreover, the fivefold lower thermal conductivity of the gold alloy, as compared to that of pure gold, was also taken into account [28]. The set of temperature equations was solved until equilibration of electron and lattice temperatures at ≈7 ps, corresponding to the electron–lattice energy exchange timescale for the electron temperature of 0.6 eV. As a result, at that time, the transient lattice temperature exceeded the melting temperature of gold, Tmelt ≈ 1337 K [29]. The laserinduced film melting resulted in significant reduction of its thermal conductivity [30], accounted for in an additional equation for melt temperature. Also, on this timescale, the film was intensely cooled by surface evaporation [5, 18–23], as was accounted for by the boundary conditions [21]. The temperature profiles calculated for the 60nm and 200nm films for the 200fs laser pulse are shown in Fig. 3. In the case of the 200nm film, at the laser intensity of 10 TW/cm2 and the time instant of 13 ps, the temperature profile demonstrates a pronounced maximum at a distance of 50–60 nm from the outer surface with its magnitude
considerably exceeding the corresponding surface temperature, as well as the normal boiling temperature of gold Tboil ≈ 3150 K [29], in agreement with surface nanotopographies presented in Figs. 2a and 2d. The results of the numerical calculation for the 60nm film at the laser intensity of 7 TW/cm2 at the time instant of 11 ps are presented in the inset to Fig. 3. As seen from the inset, the temperature maximum within the thin ner film is absent owing to less intense evaporative cooling and efficient heat conduction across the film [25, 31]. 4. Therefore, in this study, pronounced effects of subsurface boiling were observed for the first time dur ing singleshot femtosecond laser ablation of sup ported gold films in the regimes of fabrication of microbumps and nanospikes, supporting the evapora tive origin of the latter surface features. Our theoretical analysis and numerical simulations demonstrate that subsurface boiling becomes possible in relatively thick films owing to subnanosecond appearance of a sub surface temperature maximum as a result of a balance between electronic heat conduction and surface evap orative cooling in the films. This work was supported by the Russian Founda tion for Basic Research (project nos. 130200971a, 130201377a, 140200460a, 140200748a, 14 3250026 mol_nr), by the Presidium of the Russian Academy of Sciences, by the ITMO PostDoctoral Fellowship Program (government subsidy no. 074 U01), and by the Ministry of Education and Science of the Russian Federation (resolution no. P218, con tract no. 02.G25.31.0116 on August 14, 2014, between Open Joint Stock Company Ship Repair Center JETP LETTERS
Vol. 101
No. 6
2015
NANOSCALE BOILING DURING SINGLESHOT FEMTOSECOND LASER
Dalzavod and the Ministry of Education and Science of the Russian Federation). REFERENCES 1. Yu. N. Kulchin, O. B. Vitrik, A. A. Kuchmizhak, V. I. Emel’yanov, A. A. Ionin, S. I. Kudryashov, and S. V. Makarov, Phys. Rev. E 90, 023017 (2014). 2. Yu. N. Kulchin, O. B. Vitrik, A. A. Kuchmizhak, A. G. Savchuk, A. A. Nepomnyashchii, P. A. Danilov, D. A. Zayarniy, A. A. Ionin, S. I. Kudryashov, S. V. Makarov, A. A. Rudenko, V. I. Yurovskikh, and A. A. Samokhin, J. Exp. Theor. Phys. 119, 15 (2014). 3. J. P. Moening, S. S. Thanawala, and D. G. Georgiev, Appl. Phys. A 95, 635 (2009). 4. Yu. N. Kulchin, O. B. Vitrik, A. A. Kuchmizhak, A. V. Nepomnyashchii, A. G. Savchuk, A. A. Ionin, S. I. Kudryashov, and S. V. Makarov, Opt. Lett. 38, 1452 (2013). 5. V. I. Mazhukin, A. A. Samokhin, M. M. Demin, and A. V. Shapranov, Quantum Electron. 44, 283 (2014). 6. Y. Nakata, N. Miyanaga, and T. Okada, Appl. Surf. Sci. 253, 6555 (2007). 7. A. I. Kuznetsov, J. Koch, and B. N. Chichkov, Opt. Express 17, 18820 (2009). 8. C. Unger, J. Koch, L. Overmeyer, and B. N. Chichkov, Opt. Express 20, 24864 (2012). 9. V. I. Emel’yanov, D. A. Zayarniy, A. A. Ionin, S. I. Kud ryashov, S. V. Makarov, A. A. Rudenko, D. I. Shikunov, and V. I. Yurovskikh, JETP Lett. 100, 145 (2014). 10. P. A. Danilov, D. A. Zayarniy, A. A. Ionin, S. I. Kudr yashov, S. V. Makarov, A. A. Rudenko, V. I. Yurovskikh, Yu. N. Kulchin, O. B. Vitrik, A. A. Kuchmizhak, E. A. Drozdova, and S. B. Odinokov, Quantum Elec tron. 44, 540 (2014). 11. V. I. Emel’yanov, D. A. Zayarniy, A. A. Ionin, I. V. Kise leva, S. I. Kudryashov, S. V. Makarov, A. A. Rudenko, and Ch. T. Kh. Nguen, JETP Lett. 99, 518 (2014). 12. P. A. Danilov, E. A. Drozdova, A. A. Ionin, S. I. Kudr yashov, S. B. Odinokov, A. A. Rudenko, V. I. Yurovskih, and D. A. Zayarny, Appl. Phys. A 117, 981 (2014). 13. A. A. Ionin, S. I. Kudryashov, L. V. Seleznev, and D. V. Sinitsyn, JETP Lett. 94, 753 (2011).
JETP LETTERS
Vol. 101
No. 6
2015
397
14. A. A. Ionin, S. I. Kudryashov, L. V. Seleznev, D. V. Sin itsyn, A. F. Bunkin, V. N. Lednev, and S. M. Pershin, J. Exp. Theor. Phys. 116, 347 (2013). 15. I. A. Artyukov, D. A. Zayarniy, A. A. Ionin, S. I. Kudr yashov, S. V. Makarov, and P. N. Saltuganov, JETP Lett. 99, 51 (2014). 16. http://www.amplitudesystemes.com/satsumafiberlaser. html 17. M. A. Gubko, A. A. Ionin, S. I. Kudryashov, S. V. Ma karov, A. A. Rudenko, L. V. Seleznev, and D. V. Sin itsyn, JETP Lett. 97, 599 (2013). 18. S. I. Anisimov and Ya. A. Imas, Action of HighPower Radiation on Metals (Nauka, Moscow, 1970; National Technical Information Service, Springfield, VA, 1971). 19. F. W. Dabby and U. C. Paek, IEEE J. Quantum. Elec tron. 8, 106 (1972). 20. A. A. Samokhin, Kvant. Elektron. 1, 2056 (1974). 21. R. K. Singh, D. Bhattacharya, and J. Narayan, Appl. Phys. Lett. 57, 2022 (1990). 22. Y. V. Afanasiev, B. N. Chichkov, V. A. Isakov, A. P. Ka navin, and S. A. Uryupin, J. Russ. Laser Res. 20, 189 (1999). 23. V. I. Emel’yanov, in Laser Ablation in Liquids, Ed. by G. Yang (Pan Stanford, Singapore, 2012), ch. 1. 24. L. V. Zhigilei and V. Celli, Phys. Rev. B 77, 075113 (2008). 25. S. G. Bezhanov, A. P. Kanavin, and S. A. Uryupin, Quantum Electron. 44, 859 (2014). 26. N. A. Inogamov, V. V. Zhakhovskii, S. I. Ashitkov, V. A. Khokhlov, Yu. V. Petrov, P. S. Komarov, M. B. Ag ranat, S. I. Anisimov, and K. Nishihara, Appl. Surf. Sci. 255, 9712 (2009). 27. S. G. Bezhanov, A. A. Ionin, A. P. Kanavin, S. I. Kudr yashov, S. V. Makarov, L. V. Seleznev, D. V. Sinitsyn, and S. A. Uryupin, Zh. Eksp. Teor. Fiz. 147 (2015, in press). 28. C. Y. Ho, M. W. Ackerman, K. Y. Wu, S. G. Oh, and T. N. Havill, J. Phys. Chem. Ref. Data 7, 959 (1978). 29. I. S. Grigor’ev and E. Z. Meilikhov, Physical Values (Energoatomizdat, Moscow, 1991) [in Russian]. 30. C. Y. Ho, R. W. Powell, and P. E. Liley, J. Phys. Chem. Ref. Data 1, 279 (1972). 31. S. G. Bezhanov, A. P. Kanavin, and S. A. Uryupin, Phys. Lett. A 378, 975 (2014).