Access the full text.
Sign up today, get DeepDyve free for 14 days.
E Goulielmakis, ZH Loh, A Wirth, R Santra, N Rohringer (2010)
Real-time observation of valence electron motionNature, 466
N Leng, L Jiang, X Li, CC Xu, PJ Liu (2012)
Femtosecond laser processing of fused silica and aluminum based on electron dynamics control by shaping pulse trainsAppl Phys A, 109
HF Hu, XL Wang, HC Zhai, N Zhang, P Wang (2010)
Generation of multiple stress waves in silica glass in high fluence femtosecond laser ablationAppl Phys Lett, 97
L Jiang, DW Ying, X Li, YF Lu (2012)
Two-step femtosecond laser pulse train fabrication of nanostructured substrates for highly surface-enhanced Raman scatteringOpt Lett, 37
R Fazeli (2017)
Enhanced X-ray emission from laser-produced gold plasma by double pulses irradiation of nano-porous targetsPhys Lett A, 381
D Yu, L Jiang, F Wang, LT Qu, YF Lu (2016)
First-principles calculation of multiphoton absorption cross section of α-quartz under femtosecond laser irradiationAppl Phys A, 122
G Sansone, F Kelkensberg, JF Pérez-Torres, F Morales, MF Kling (2010)
Electron localization following attosecond molecular photoionizationNature, 465
Mengmeng Wang, Sumei Wang, Zhitao Cao, Peng Wang, Cong Wang (2015)
Investigation of double-pulse femtosecond laser induced breakdown spectroscopy of polymethyl methacrylate (PMMA), 9351
W Lee, SA Pruzinsky, PV Braun (2002)
Multi-photon polymerization of waveguide structures within three-dimensional photonic crystalsAdv Mater, 14
L Jiang, HL Tsai (2005)
Improved two-temperature model and its application in ultrashort laser heating of metal filmsJ Heat Transfer, 127
P Russo, R Liang, E Jabari, E Marzbanrad, E Toyserkani (2016)
Single-step synthesis of graphene quantum dots by femtosecond laser ablation of graphene oxide dispersionsNanoscale, 8
MA Bandres, BM Rodríguez-Lara (2013)
Nondiffracting accelerating waves: weber waves and parabolic momentumNew J Phys, 15
J Reif, F Costache, M Henyk, SV Pandelov (2002)
Ripples revisited: non-classical morphology at the bottom of femtosecond laser ablation craters in transparent dielectricsAppl Surface Sci, 197-198
RD Murphy, B Torralva, DP Adams, SM Yalisove (2013)
Pump-probe imaging of laser-induced periodic surface structures after ultrafast irradiation of SiAppl Phys Lett, 103
CC Xu, L Jiang, N Leng, YP Yuan, PJ Liu (2013)
Ultrafast laser ablation size and recast adjustment in dielectrics based on electron dynamicscontrol by pulse train shapingChin Opt Lett, 11
CH Lin, L Jiang, YH Chai, H Xiao, SJ Chen (2009)
One-step fabrication of nanostructures by femtosecond laser for surface-enhanced Raman scatteringOpt Express, 17
AY Vorobyev, CL Guo (2008)
Colorizing metals with femtosecond laser pulsesAppl Phys Lett, 92
I Choi, HY Jeong, H Shin, G Kang, M Byun (2016)
Laser-induced phase separation of silicon carbideNat Commun, 7
CSR Nathala, A Ajami, AA Ionin, SI Kudryashov, SV Makarov (2015)
Experimental study of fs-laser induced sub-100-nm periodic surface structures on titaniumOpt Express, 23
E Axente, S Noël, J Hermann, M Sentis, IN Mihailescu (2008)
Correlation between plasma expansion and damage threshold by femtosecond laser ablation of fused silicaJ Phys D Appl Phys, 41
B Xia, L Jiang, XW Li, XL Yan, YF Lu (2015)
Mechanism and elimination of bending effect in femtosecond laser deep-hole drillingOpt Express, 23
WK Zhang, BN Markiewicz, RS Doerksen, AB Smith, F Gai (2016)
C≡N stretching vibration of 5-cyanotryptophan as an infrared probe of protein local environment: what determines its frequency?Phys Chem Chem Phys, 18
B Dusser, Z Sagan, H Soder, N Faure, JP Colombier (2010)
Controlled nanostructrures formation by ultra fast laser pulses for color markingOpt Express, 18
J. Parkinson, D. Farnell (2010)
Antiferromagnetic Spin Waves
R Le Harzic, D Dörr, D Sauer, M Neumeier, M Epple (2011)
Large-area, uniform, high-spatial-frequency ripples generated on silicon using a nanojoule-femtosecond laser at high repetition rateOpt Lett, 36
K Toyoda, K Miyamoto, N Aoki, R Morita, T Omatsu (2012)
Using optical vortex to control the chirality of twisted metal nanostructuresNano Lett, 12
P Hockett, CZ Bisgaard, OJ Clarkin, A Stolow (2011)
Time-resolved imaging of purely valence-electron dynamics during a chemical reactionNat Phys, 7
J Gottmann, M Hermans, N Repiev, J Ortmann (2017)
Selective laser-induced etching of 3D precision quartz glass components for microfluidic applications—up-scaling of complexity and speedMicromachines, 8
F Courvoisier, PA Lacourt, M Jacquot, MK Bhuyan, L Furfaro (2009)
Surface nanoprocessing with nondiffracting femtosecond Bessel beamsOpt Lett, 34
YP Yuan, L Jiang, X Li, C Wang, YF Lu (2012)
Adjustment of ablation shapes and subwavelength ripples based on electron dynamics control by designing femtosecond laser pulse trainsJ Appl Phys, 112
B Xia, L Jiang, X Li, X Yan, W Zhao (2015)
High aspect ratio, high-quality microholes in PMMA: a comparison between femtosecond laser drilling in air and in vacuumAppl Phys A, 119
A Borowiec, HK Haugen (2003)
Subwavelength ripple formation on the surfaces of compound semiconductors irradiated with femtosecond laser pulsesAppl Phys Lett, 82
R Kammel, R Ackermann, J Thomas, J Götte, S Skupin (2014)
Enhancing precision in fs-laser material processing by simultaneous spatial and temporal focusingLight Sci Appl, 3
Y Hayasaki, K Iwata, S Hasegawa, A Takita, S Juodkazis (2011)
Time-resolved axial-view of the dielectric breakdown under tight focusing in glassOpt Mater Express, 1
L Jiang, HL Tsai (2006)
Energy transport and nanostructuring of dielectrics by femtosecond laser pulse trainsJ Heat Transfer, 128
AH Zewail (1988)
Laser femtochemistryScience, 242
AQ Wu, XF Xu (2007)
2007. First International Conference on Integration and Commercialization of MICRO and Nanosystems; 10–13 January 2007; Sanya, Hainan, China
K Okamuro, M Hashida, Y Miyasaka, Y Ikuta, S Tokita (2010)
Laser fluence dependence of periodic grating structures formed on metal surfaces under femtosecond laser pulse irradiationPhys Rev B, 82
J Bonse, A Rosenfeld, J Krüger (2009)
On the role of surface plasmon polaritons in the formation of laser-induced periodic surface structures upon irradiation of silicon by femtosecond-laser pulsesJ Appl Phys, 106
P Nagpal, NC Lindquist, SH Oh, DJ Norris (2009)
Ultrasmooth patterned metals for plasmonics and metamaterialsScience, 325
EAA Pogna, M Marsili, D De Fazio, S Dal Conte, C Manzoni (2016)
Photo-induced bandgap renormalization governs the ultrafast response of single-layer MoS2ACS Nano, 10
L Jiang, HL Tsai (2005)
Proceedings Volume 5713, Photon Processing in Microelectronics and Photonics IV; 12 April 2005; San Jose, California, United States
DC Papazoglou, S Tzortzakis (2008)
In-line holography for the characterization of ultrafast laser filamentation in transparent mediaAppl Phys Lett, 93
MJ Zhao, J Hu, L Jiang, KH Zhang, PJ Liu (2015)
Controllable high-throughput high-quality femtosecond laser-enhanced chemical etching by temporal pulse shaping based on electron density controlSci Rep, 5
YC Chen, PS Salter, S Knauer, LY Weng, AC Frangeskou (2016)
Laser writing of coherent colour centres in diamondNat Photon, 11
XF Li, CY Zhang, H Li, QF Dai, S Lan (2014)
Formation of 100-nm periodic structures on a titanium surface by exploiting the oxidation and third harmonic generation induced by femtosecond laser pulsesOpt Express, 22
DW Li, YS Zhou, X Huang, L Jiang, JF Silvain (2015)
In situ imaging and control of layer-by-layer femtosecond laser thinning of grapheneNanoscale, 7
H Park, X Wang, S Nie, R Clinite, J Cao (2005)
Mechanism of coherent acoustic phonon generation under nonequilibrium conditionsPhys Rev B, 72
W Han, YF Yang, W Cheng, QW Zhan (2013)
Vectorial optical field generator for the creation of arbitrarily complex fieldsOpt Express, 21
D Wu, J Xu, LG Niu, SZ Wu, K Midorikawa (2015)
In-channel integration of designable microoptical devices using flat scaffold-supported femtosecond-laser microfabrication for coupling-free optofluidic cell countingLight Sci Appl, 4
AD Wang, L Jiang, XW Li, Y Liu, XZ Dong (2015)
Mask-free patterning of high-conductivity metal nanowires in open air by spatially modulated femtosecond laser pulsesAdv Mater, 27
O Smirnova, Y Mairesse, S Patchkovskii, N Dudovich, D Villeneuve (2009)
High harmonic interferometry of multi-electron dynamics in moleculesNature, 460
RV Mikhaylovskiy, E Hendry, A Secchi, JH Mentink, M Eckstein (2015)
Ultrafast optical modification of exchange interactions in iron oxidesNat Commun, 6
S Höhm, A Rosenfeld, J Krüger, J Bonse (2013)
Femtosecond diffraction dynamics of laser-induced periodic surface structures on fused silicaAppl Phys Lett, 102
YP Deng, XH Xie, H Xiong, YX Leng, CF Cheng (2005)
Optical breakdown for silica and silicon with double femtosecond laser pulsesOpt Express, 13
S Höhm, A Rosenfeld, J Krüger, J Bonse (2012)
Femtosecond laser-induced periodic surface structures on silicaJ Appl Phys, 112
DG Papazoglou, S Suntsov, D Abdollahpour, S Tzortzakis (2010)
Tunable intense Airy beams and tailored femtosecond laser filamentsPhys Rev A, 81
M Malinauskas (2016)
e16133Light Sci Appl, 5
JF Li, YF Huang, Y Ding, ZL Yang, SB Li (2010)
Shell-isolated nanoparticle-enhanced Raman spectroscopyNature, 464
I Martín-Fabiani, E Rebollar, S Pérez, DR Rueda, MC García-Gutiérrez (2012)
Laser-induced periodic surface structures nanofabricated on poly(trimethylene terephthalate) spin-coated filmsLangmuir, 28
Dong Wu, Jian Xu, L. Niu, Si-zhu Wu, K. Midorikawa, K. Sugioka (2015)
Correction: Corrigendum: In-channel integration of designable microoptical devices using flat scaffold-supported femtosecond-laser microfabrication for coupling-free optofluidic cell countingLight: Science & Applications, 4
I Paradisanos, E Kymakis, C Fotakis, G Kioseoglou, E Stratakis (2014)
Intense femtosecond photoexcitation of bulk and monolayer MoS2Appl Phys Lett, 105
SK Sundaram, E Mazur (2002)
Inducing and probing non-thermal transitions in semiconductors using femtosecond laser pulsesNat Mater, 1
M Shinoda, RR Gattass, E Mazur (2009)
Femtosecond laser-induced formation of nanometer-width grooves on synthetic single-crystal diamond surfacesJ Appl Phys, 105
J Bonse, J Krüger (2010)
Pulse number dependence of laser-induced periodic surface structures for femtosecond laser irradiation of siliconJ Appl Phys, 108
XS Shi, L Jiang, X Li, SM Wang, YP Yuan (2013)
Femtosecond laser-induced periodic structure adjustments based on electron dynamics control: from subwavelength ripples to double-grating structuresOpt Lett, 38
G Dumitru, V Romano, HP Weber, M Sentis, W Marine (2002)
Femtosecond ablation of ultrahard materialsAppl Phys A, 74
M Renard, E Hertz, B Lavorel, O Faucher (2004)
Controlling ground-state rotational dynamics of molecules by shaped femtosecond laser pulsesPhys Rev A, 69
A. Wu, Xianfan Xu (2007)
Coherent Phonon Excitation by Ultrafast Pulse Trains
XD Chen, W Xin, WS Jiang, ZB Liu, YS Chen (2016)
High-precision twist-controlled bilayer and trilayer grapheneAdv Mater, 28
R Stoian, M Boyle, A Thoss, A Rosenfeld, G Korn (2002)
Laser ablation of dielectrics with temporally shaped femtosecond pulsesAppl Phys Lett, 80
G Seifert, M Kaempfe, F Syrowatka, C Harnagea, D Hesse (2005)
Self-organized structure formation on the bottom of femtosecond laser ablation craters in glassAppl Phys A, 81
T Brixner, G Gerber (2001)
Femtosecond polarization pulse shapingOpt Lett, 26
SS Hou, YY Huo, PX Xiong, Y Zhang, SA Zhang (2011)
Formation of long- and short-periodic nanoripples on stainless steel irradiated by femtosecond laser pulsesJ Phys D Appl Phys, 44
JL Wang, L Guo, CH Liu, X Xu, YF Chen (2015)
Influence of coherent optical phonon on ultrafast energy relaxationAppl Phys Lett, 107
D Esser, S Rezaei, JZ Li, PR Herman, J Gottmann (2011)
Time dynamics of burst-train filamentation assisted femtosecond laser machining in glassesOpt Express, 19
H. Tsai, Lan Jiang (2005)
Fundamentals of energy cascade during ultrashort laser-material interactions (Invited Paper), 5713
C Kerse, H Kalaycıoğlu, P Elahi, B Çetin, DK Kesim (2016)
Ablation-cooled material removal with ultrafast bursts of pulsesNature, 537
TQ Jia, HX Chen, M Huang, FL Zhao, JR Qiu (2005)
Formation of nanogratings on the surface of a ZnSe crystal irradiated by femtosecond laser pulsesPhys Rev B, 72
A Mathis, F Courvoisier, L Froehly, L Furfaro, M Jacquot (2012)
Micromachining along a curve: femtosecond laser micromachining of curved profiles in diamond and silicon using accelerating beamsAppl Phys Lett, 101
J-P Colombier, P Combis, A Rosenfeld, IV Hertel, E Audouard (2006)
Optimized energy coupling at ultrafast laser-irradiated metal surfaces by tailoring intensity envelopes: Consequences for material removal from Al samplesPhys Rev B, 74
V Zorba, E Stratakis, M Barberoglou, E Spanakis, P Tzanetakis (2008)
Biomimetic artificial surfaces quantitatively reproduce the water repellency of a lotus leafAdv Mater, 20
K Moore, H Rabitz (2012)
Laser control: manipulating moleculesNat Chem, 4
K Sokolowski-Tinten, J Bialkowski, A Cavalleri, D von der Linde, A Oparin (1998)
Transient states of matter during short pulse laser ablationPhys Rev Lett, 81
L Rapp, R Meyer, R Giust, L Furfaro, M Jacquot (2016)
High aspect ratio micro-explosions in the bulk of sapphire generated by femtosecond Bessel beamsSci Rep, 6
TJY Derrien, J Krüger, TE Itina, S Höhm, A Rosenfeld (2014)
Rippled area formed by surface plasmon polaritons upon femtosecond laser double-pulse irradiation of silicon: the role of carrier generation and relaxation processesAppl Phys A, 117
CH Lin, ZH Rao, L Jiang, WJ Tsai, PH Wu (2010)
Investigations of femtosecond-nanosecond dual-beam laser ablation of dielectricsOpt Lett, 35
Y Shimotsuma, PG Kazansky, JR Qiu, K Hirao (2003)
Self-Organized nanogratings in glass irradiated by ultrashort light pulsesPhys Rev Lett, 91
WW Zhao, XW Li, B Xia, XL Yan, WN Han (2014)
Proceedings Volume 9296, International Symposium on Optoelectronic Technology and Application 2014: Advanced Display Technology; Nonimaging Optics: Efficient Design for Illumination and Solar Concentration; 21 November 2014; Beijing, China
X Li, L Jiang, H-L Tsai (2009)
Phase change mechanisms during femtosecond laser pulse train ablation of nickel thin filmsJ Appl Phys, 106
P. Alvarez (2018)
First Principles Calculations
K Sokolowski-Tinten, J Solis, J Bialkowski, J Siegel, C Afonso (1998)
Dynamics of ultrafast phase changes in amorphous GeSb filmsPhys Rev Lett, 81
C Coleman, R Erasmus, S Bhattacharyya (2016)
Nanoscale deformations in graphene by laser annealingAppl Phys Lett, 109
RJ Levis, GM Menkir, H Rabitz (2001)
Selective bond dissociation and rearrangement with optimally tailored, strong-field laser pulsesScience, 292
LA Girifalco, VG Weizer (1959)
Application of the Morse potential function to cubic metalsPhys Rev, 114
D Yu, L Jiang, F Wang, X Li, LT Qu (2015)
Electron ionization and spin polarization control of Fe atom adsorbed graphene irradiated by a femtosecond laserPhys Lett A, 379
YW Yu, L Jiang, Q Cao, XS Shi, QS Wang (2016)
Ultrafast imaging the light-speed propagation of a focused femtosecond laser pulse in air and its ionized electron dynamics and plasma-induced pulse reshapingAppl Phys A, 122
M Dantus, RM Bowman, AH Zewail (1990)
Femtosecond laser observations of molecular vibration and rotationNature, 343
TJY Derrien, J Krüger, TE Itina, S Höhm, A Rosenfeld (2013)
Rippled area formed by surface plasmon polaritons upon femtosecond laser double-pulse irradiation of siliconOpt Express, 21
CX Sheng, C Zhang, YX Zhai, K Mielczarek, WW Wang (2015)
Exciton versus free carrier photogeneration in organometal trihalide perovskites probed by broadband ultrafast polarization memory dynamicsPhys Rev Lett, 114
QQ Yang, X Li, L Jiang, N Zhang, GM Zhang (2015)
Nanopillar arrays with nanoparticles fabricated by a femtosecond laser pulse train for highly sensitive SERRSOpt Lett, 40
WN Han, L Jiang, XW Li, QS Wang, H Li (2014)
Anisotropy modulations of femtosecond laser pulse induced periodic surface structures on silicon by adjusting double pulse delayOpt Express, 22
MM Wang, SM Wang, ZT Cao, P Wang, C Wang (2015)
Proceedings Volume 9351, Laser-based Micro- and Nanoprocessing IX; 12 March 2015; San Francisco, California, United States
L Jiang, PJ Liu, XL Yan, N Leng, CC Xu (2012)
High-throughput rear-surface drilling of microchannels in glass based on electron dynamics control using femtosecond pulse trainsOpt Lett, 37
S Baierl, M Hohenleutner, T Kampfrath, AK Zvezdin, AV Kimel (2016)
Nonlinear spin control by terahertz-driven anisotropy fieldsNat Photon, 10
N Zhang, XN Zhu, JJ Yang, XL Wang, MW Wang (2007)
Time-resolved shadowgraphs of material ejection in intense femtosecond laser ablation of aluminumPhys Rev Lett, 99
G Auböck, M Chergui (2015)
Sub-50-fs photoinduced spin crossover in [Fe(bpy)3]2+Nat Chem, 7
M Huang, FL Zhao, Y Cheng, NS Xu, ZZ Xu (2009)
Origin of laser-induced near-subwavelength ripples: interference between surface plasmons and incident laserACS Nano, 3
S Wienholdt, D Hinzke, U Nowak (2012)
THz switching of antiferromagnets and ferrimagnetsPhys Rev Lett, 108
L Jiang, HL Tsai (2005)
Repeatable nanostructures in dielectrics by femtosecond laser pulse trainsAppl Phys Lett, 87
(2016)
Ultrafast laser processing of materials: from science to industry, 5
H Raether (1988)
Surface Plasmons on Smooth and Rough Surfaces and on Gratings
T Kampfrath, A Sell, G Klatt, A Pashkin, S Mahrlein (2011)
Coherent terahertz control of antiferromagnetic spin wavesNat Photon, 5
L Jiang, XS Shi, X Li, YP Yuan, C Wang (2012)
Subwavelength ripples adjustment based on electron dynamics control by using shaped ultrafast laser pulse trainsOpt Express, 20
R Fazeli, MH Mahdieh (2015)
Comparison of line X-ray emission from solid and porous nano-layer coated targets irradiated by double laser pulsesPhys Plasmas, 22
M Barberoglou, D Gray, E Magoulakis, C Fotakis, PA Loukakos (2013)
Controlling ripples’ periodicity using temporally delayed femtosecond laser double pulsesOpt Express, 21
JT Chen, WC Lai, YJ Kao, YY Yang, JK Sheu (2012)
Laser-induced periodic structures for light extraction efficiency enhancement of GaN-based light emitting diodesOpt Express, 20
A Marcinkevičius, S Juodkazis, S Matsuo, V Mizeikis, H Misawa (2001)
Application of Bessel beams for microfabrication of dielectrics by femtosecond laserJpn J Appl Phys, 40
N Zhang, X Li, L Jiang, XS Shi, C Li (2013)
Femtosecond double-pulse fabrication of hierarchical nanostructures based on electron dynamics control for high surface-enhanced Raman scatteringOpt Lett, 38
H Rabitz, R de Vivie-Riedle, M Motzkus, K Kompa (2000)
Whither the future of controlling quantum phenomena?Science, 288
AH Zewail (1996)
Femtochemistry: recent progress in studies of dynamics and control of reactions and their transition statesJ Phys Chem, 100
B Poumellec, M Lancry, R Desmarchelier, E Hervé, B Bourguignon (2016)
Parity violation in chiral structure creation under femtosecond laser irradiation in silica glass?Light Sci Appl, 5
F Liang, R Vallée, SL Chin (2012)
Mechanism of nanograting formation on the surface of fused silicaOpt Express, 20
HN Wang, CJ Zhang, F Rana (2014)
Ultrafast dynamics of defect-assisted electron-hole recombination in monolayer MoS2Nano Lett, 15
E Le Ru, P Etchegoin (2008)
Principles of Surface-Enhanced Raman Spectroscopy: and Related Plasmonic Effects
SW Bahk, P Rousseau, TA Planchon, V Chvykov, G Kalintchenko (2004)
Generation and characterization of the highest laser intensities (1022 W/cm2Opt Lett, 29
T Kampfrath, K Tanaka, KA Nelson (2013)
Resonant and nonresonant control over matter and light by intense terahertz transientsNat Photon, 7
L Jiang, HL Tsai (2008)
A plasma model combined with an improved two-temperature equation for ultrafast laser ablation of dielectricsJ Appl Phys, 104
A Otto, I Mrozek, H Grabhorn, W Akemann (1992)
Surface-enhanced Raman scatteringJ Phys Condens Matter, 4
B Barwick, DJ Flannigan, AH Zewail (2009)
Photon-induced near-field electron microscopyNature, 462
XQ Liu, L Yu, QD Chen, HB Sun (2017)
Mask-free construction of three-dimensional silicon structures by dry etching assisted gray-scale femtosecond laser direct writingAppl Phys Lett, 110
F Carbone, OH Kwon, AH Zewail (2009)
Dynamics of chemical bonding mapped by energy-resolved 4D electron microscopyScience, 325
A Tarun, MRH Daza, N Hayazawa, Y Inouye, S Kawata (2002)
Apertureless optical near-field fabrication using an atomic force microscope on photoresistsAppl Phys Lett, 80
S Richter, M Heinrich, S Döring, A Tünnermann, S Nolte (2012)
Nanogratings in fused silica: formation, control, and applicationsJ Laser Appl, 24
YP Yuan, L Jiang, X Li, C Wang, LT Qu (2013)
Simulation of rippled structure adjustments based on localized transient electron dynamics control by femtosecond laser pulse trainsAppl Phys A, 111
AH Zewail (2000)
Femtochemistry: atomic-scale dynamics of the chemical bondJ Phys Chem A, 104
P Zuo, L Jiang, X Li, B Li, YD Xu (2017)
Shape-controllable gold nanoparticles-MoS2 hybrids prepared by tuning edge-active sites and surface structures of MoS2 via temporally shaped femtosecond pulsesACS Appl Mater Interfaces, 9
PX Fan, BF Bai, JY Long, DF Jiang, GF Jin (2015)
Broadband high-performance infrared antireflection nanowires facilely grown on ultrafast laser structured Cu surfaceNano Lett, 15
S Döring, J Szilagyi, S Richter, F Zimmermann, M Richardson (2012)
Evolution of hole shape and size during short and ultrashort pulse laser deep drillingOpt Express, 20
R Le Harzic, D Dörr, D Sauer, F Stracke, H Zimmermann (2011)
Generation of high spatial frequency ripples on silicon under ultrashort laser pulses irradiationAppl Phys Lett, 98
M Duocastella, CB Arnold (2012)
Bessel and annular beams for materials processingLaser Photon Rev, 6
AA Andreev, J Limpouch, AB Iskakov, H Nakano (2002)
Enhancement of X-ray line emission from plasmas produced by short high-intensity laser double pulsesPhys Rev E, 65
C Wang, L Jiang, F Wang, X Li, YP Yuan (2012)
First-principles electron dynamics control simulation of diamond under femtosecond laser pulse train irradiationJ Phys Condens Matter, 24
M Sakakura, M Terazima (2005)
Initial temporal and spatial changes of the refractive index induced by focused femtosecond pulsed laser irradiation inside a glassPhys Rev B, 71
ME Corrales, J González-Vázquez, G Balerdi, IR Solá, R De Nalda (2014)
Control of ultrafast molecular photodissociation by laser-field-induced potentialsNat Chem, 6
J Guo, TF Wang, JF Shao, T Sun, R Wang (2012)
Emission enhancement ratio of the metal irradiated by femtosecond double-pulse laserOpt Commun, 285
ZH Wang, B Zeng, GH Li, HQ Xie, W Chu (2015)
Time-resolved shadowgraphs of transient plasma induced by spatiotemporally focused femtosecond laser pulses in fused silica glassOpt Lett, 40
Q Xie, XW Li, L Jiang, B Xia, XL Yan (2016)
High-aspect-ratio, high-quality microdrilling by electron density control using a femtosecond laser Bessel beamAppl Phys A, 122
A Žukauskas, M Malinauskas, E Brasselet (2013)
Monolithic generators of pseudo-nondiffracting optical vortex beams at the microscaleAppl Phys Lett, 103
HC Yuan, VE Yost, MR Page, P Stradins, DL Meier (2009)
Efficient black silicon solar cell with a density-graded nanoporous surface: optical properties, performance limitations, and design rulesAppl Phys Lett, 95
RR Fang, A Vorobyev, CL Guo (2017)
Direct visualization of the complete evolution of femtosecond laser-induced surface structural dynamics of metalsLight Sci Appl, 6
EG Gamaly (2011)
Femtosecond Laser-Matter Interaction: Theory, Experiments and Applications
S Odoulov, A Shumelyuk, H Badorreck, S Nolte, KM Voit (2015)
Interference and holography with femtosecond laser pulses of different coloursNat Commun, 6
J Kim, S Na, S Cho, W Chang, K Whang (2008)
Surface ripple changes during Cr film ablation with a double ultrashort laser pulseOpt Lasers Eng, 46
A Assion, T Baumert, M Bergt, T Brixner, B Kiefer (1998)
Control of chemical reactions by feedback-optimized phase-shaped femtosecond laser pulsesScience, 282
P Lorazo, LJ Lewis, M Meunier (2003)
Short-pulse laser ablation of solids: from phase explosion to fragmentationPhys Rev Lett, 91
MK Bhuyan, F Courvoisier, PA Lacourt, M Jacquot, R Salut (2010)
High aspect ratio nanochannel machining using single shot femtosecond Bessel beamsAppl Phys Lett, 97
(2014)
Ultrafast lasers—reliable tools for advanced materials processing, 3
YW Yu, L Jiang, Q Cao, B Xia, QS Wang (2015)
Pump-probe imaging of the fs-ps-ns dynamics during femtosecond laser Bessel beam drilling in PMMAOpt Express, 23
S Rezaei, JZ Li, PR Herman (2015)
Burst train generator of high energy femtosecond laser pulses for driving heat accumulation effect during micromachiningOpt Lett, 40
A Konar, YN Shu, VV Lozovoy, JE Jackson, BG Levine (2014)
Polyatomic molecules under intense femtosecond laser irradiationJ Phys Chem A, 118
B Wetzel, C Xie, PA Lacourt, JM Dudley, F Courvoisier (2013)
Femtosecond laser fabrication of micro and nano-disks in single layer graphene using vortex Bessel beamsAppl Phys Lett, 103
ZW Liu, QH Wei, X Zhang (2005)
Surface plasmon interference nanolithographyNano Lett, 5
R Improta, F Santoro, L Blancafort (2016)
Quantum mechanical studies on the photophysics and the photochemistry of nucleic acids and nucleobasesChem Rev, 116
HJ Kwon, S Chung, J Jang, CP Grigoropoulos (2016)
Laser direct writing and inkjet printing for a sub-2 μm channel length MoS2 transistor with high-resolution electrodesNanotechnology, 27
PM Kraus, OI Tolstikhin, D Baykusheva, A Rupenyan, J Schneider (2015)
Observation of laser-induced electronic structure in oriented polyatomic moleculesNat Commun, 6
JS Preston, HM Van Driel, JE Sipe (1987)
Order-disorder transitions in the melt morphology of laser-irradiated siliconPhys Rev Lett, 58
H He, SY Li, SY Wang, ML Hu, YJ Cao (2012)
Manipulation of cellular light from green fluorescent protein by a femtosecond laserNat Photon, 6
MK Bhuyan, PK Velpula, JP Colombier, T Olivier, N Faure (2014)
Single-shot high aspect ratio bulk nanostructuring of fused silica using chirp-controlled ultrafast laser Bessel beamsAppl Phys Lett, 104
K Sugioka (2014)
e149Light Sci Appl, 3
J Bonse, M Munz, H Sturm (2005)
Structure formation on the surface of indium phosphide irradiated by femtosecond laser pulsesJ Appl Phys, 97
RS Taylor, C Hnatovsky, E Simova, DM Rayner, VR Bhardwaj (2003)
Femtosecond laser fabrication of nanostructures in silica glassOpt Lett, 28
M Domke, S Rapp, M Schmidt, HP Huber (2012)
Ultrafast pump-probe microscopy with high temporal dynamic rangeOpt Express, 20
L Wang, B-B Xu, X-W Cao, Q-K Li, W-J Tian (2017)
Competition between subwavelength and deep-subwavelength structures ablated by ultrashort laser pulsesOptica, 4
A Lindinger, C Lupulescu, M Plewicki, F Vetter, A Merli (2004)
Isotope selective ionization by optimal control using shaped femtosecond laser pulsesPhys Rev Lett, 93
JZ Li, S Ho, M Haque, PR Herman (2012)
Nanograting Bragg responses of femtosecond laser written optical waveguides in fused silica glassOpt Mater Express, 2
D. Meisel, M. Deubel, M. Hermatschweiler, K. Busch, W. Koch, G. Freymann, A. Blanco, C. Enkrich, M. Wegener (2004)
Three-Dimensional Photonic CrystalsSolid State Phenomena, 99-100
RR Gattass, E Mazur (2008)
Femtosecond laser micromachining in transparent materialsNat Photon, 2
D Dufft, A Rosenfeld, SK Das, R Grunwald, J Bonse (2009)
Femtosecond laser-induced periodic surface structures revisited: a comparative study on ZnOJ Appl Phys, 105
X Zhao, YC Shin (2011)
Femtosecond laser drilling of high-aspect ratio microchannels in glassAppl Phys A, 104
L Gao, F Wang, L Jiang, LT Qu, YF Lu (2015)
Controlling the excitation process of free electrons by a femtosecond elliptically polarized laserInt J Modern Phys B, 29
KH Zhang, L Jiang, X Li, XS Shi, D Yu (2014)
Femtosecond laser pulse-train induced breakdown in fused silica: the role of seed electronsJ Phys D Appl Phys, 47
X Zhao, YC Shin (2013)
Coulomb explosion and early plasma generation during femtosecond laser ablation of silicon at high laser fluenceJ Phys D Appl Phys, 46
S Han, S Hong, J Ham, J Yeo, J Lee (2014)
Fast plasmonic laser nanowelding for a Cu-nanowire percolation network for flexible transparent conductors and stretchable electronicsAdv Mater, 26
PDC King, HI Wei, YF Nie, M Uchida, C Adamo (2014)
Atomic-scale control of competing electronic phases in ultrathin LaNiO3Nat Nanotechnol, 9
X Li, L Jiang (2012)
Size distribution control of metal nanoparticles using femtosecond laser pulse train: a molecular dynamics simulationAppl Phys A, 109
F Ladieu, P Martin, S Guizard (2002)
Measuring thermal effects in femtosecond laser-induced breakdown of dielectricsAppl Phys Lett, 81
G Miyaji, K Miyazaki (2008)
Origin of periodicity in nanostructuring on thin film surfaces ablated with femtosecond laser pulsesOpt Express, 16
C Wang, L Jiang, F Wang, X Li, YP Yuan (2012)
Transient localized electron dynamics simulation during femtosecond laser tunnel ionization of diamondPhys Lett A, 376
C Wang, L Jiang, X Li, F Wang, YP Yuan (2013)
Frequency dependence of electron dynamics during femtosecond laser resonant photoionization of Li4 clusterJ Appl Phys, 114
J Jersch, K Dickmann (1996)
Nanostructure fabrication using laser field enhancement in the near field of a scanning tunneling microscope tipAppl Phys Lett, 68
SS Harilal, PK Diwakar, A Hassanein (2013)
Electron-ion relaxation time dependent signal enhancement in ultrafast double-pulse laser-induced breakdown spectroscopyAppl Phys Lett, 103
LL Gao, F Wang, L Jiang, LT Qu, YF Lu (2015)
Optical-induced electrical current in diamond switched by femtosecond-attosecond laser pulses by ab initio simulationsJ Phys D Appl Phys, 49
HY Shi, RS Yan, S Bertolazzi, J Brivio, B Gao (2013)
Exciton dynamics in suspended monolayer and few-layer MoS2 2D crystalsACS Nano, 7
GS Su, L Jiang, F Wang, LT Qu, YF Lu (2016)
First-principles simulations for excitation of currents in linear carbon chains under femtosecond laser pulse irradiationPhys Lett A, 380
C Brif, R Chakrabarti, H Rabitz (2010)
Control of quantum phenomena: past, present, and futureNew J Phys, 12
M Straub, M Afshar, D Feili, H Seidel, K König (2012)
Surface plasmon polariton model of high-spatial frequency laser-induced periodic surface structure generation in siliconJ Appl Phys, 111
C Wang, L Jiang, F Wang, X Li, YP Yuan (2011)
First-principles calculations of the electron dynamics during femtosecond laser pulse train material interactionsPhys Lett A, 375
XP Hong, J Kim, SF Shi, Y Zhang, CH Jin (2014)
Ultrafast charge transfer in atomically thin MoS2/WS2 heterostructuresNat Nanotechnol, 9
C Wang, L Jiang, X Li, F Wang, YP Yuan (2013)
Nonlinear ionization mechanism dependence of energy absorption in diamond under femtosecond laser irradiationJ Appl Phys, 113
F De Angelis, F Gentile, F Mecarini, G Das, M Moretti (2011)
Breaking the diffusion limit with super-hydrophobic delivery of molecules to plasmonic nanofocusing SERS structuresNat Photon, 5
YY Jau, A Hankin, T Keating, IH Deutsch, G Biedermann (2016)
Entangling atomic spins with a Rydberg-dressed spin-flip blockadeNat Phys, 12
SE Clark, DC Emmony (1989)
Ultraviolet-laser-induced periodic surface structuresPhys Rev B, 40
CJR Sheppard, T Wilson (1978)
Gaussian-beam theory of lenses with annular apertureIEE J Microw Opt Acoustics, 2
A Rosenfeld, M Rohloff, S Höhm, J Krüger, J Bonse (2012)
Formation of laser-induced periodic surface structures on fused silica upon multiple parallel polarized double-femtosecond-laser-pulse irradiation sequencesAppl Surface Sci, 258
S Karimelahi, L Abolghasemi, PR Herman (2014)
Rapid micromachining of high aspect ratio holes in fused silica glass by high repetition rate picosecond laserAppl Phys A, 114
VV Temnov, K Sokolowski-Tinten, P Zhou, A El-Khamhawy, D von der Linde (2006)
Multiphoton ionization in dielectrics: comparison of circular and linear polarizationPhys Rev Lett, 97
U Quentin, KH Leitz, L Deichmann, I Alexeev, M Schmidt (2012)
Optical trap assisted laser nanostructuring in the near-field of microparticlesJ Laser Appl, 24
Weiwei Zhao, Xiaowei Li, B. Xia, Xueliang Yan, W. Han, Y. Lu, Lan Jiang (2014)
Single-pulse femtosecond laser Bessel beams drilling of high-aspect-ratio microholes based on electron dynamics control, 9296
A Ben-Yakar, A Harkin, J Ashmore, RL Byer, HA Stone (2007)
Thermal and fluid processes of a thin melt zone during femtosecond laser ablation of glass: the formation of rims by single laser pulsesJ Phys D Appl Phys, 40
YL Sun, WF Dong, LG Niu, T Jiang, DX Liu (2014)
Protein-based soft micro-optics fabricated by femtosecond laser direct writingLight Sci Appl, 3
F He, Y Cheng, ZZ Xu, Y Liao, J Xu (2010)
Direct fabrication of homogeneous microfluidic channels embedded in fused silica using a femtosecond laserOpt Lett, 35
YP Yuan, L Jiang, X Li, C Wang, H Xiao (2012)
Formation mechanisms of sub-wavelength ripples during femtosecond laser pulse train processing of dielectricsJ Phys D Appl Phys, 45
S Richter, F Jia, M Heinrich, S Döring, U Peschel (2012)
The role of self-trapped excitons and defects in the formation of nanogratings in fused silicaOpt Lett, 37
L Jiang, HL Tsai (2005)
Energy transport and material removal in wide bandgap materials by a femtosecond laser pulseInt J Heat Mass Transfer, 48
V Hommes, M Miclea, R Hergenröder (2006)
Silicon surface morphology study after exposure to tailored femtosecond pulsesAppl Surface Sci, 252
S Juodkazis, K Nishimura, S Tanaka, H Misawa, EG Gamaly (2006)
Laser-induced microexplosion confined in the bulk of a sapphire crystal: evidence of multimegabar pressuresPhys Rev Lett, 96
S Haessler, J Caillat, W Boutu, C Giovanetti-Teixeira, T Ruchon (2010)
Attosecond imaging of molecular electronic wavepacketsNat Phys, 6
Light: Science & Applications (2018) 7, 17134; doi:10.1038/lsa.2017.134 OPEN Official journal of the CIOMP 2047-7538/18 www.nature.com/lsa REVIEW ARTICLE Electrons dynamics control by shaping femtosecond laser pulses in micro/nanofabrication: modeling, method, measurement and application 1 1 1 2 3 Lan Jiang , An-Dong Wang ,BoLi , Tian-Hong Cui and Yong-Feng Lu During femtosecond laser fabrication, photons are mainly absorbed by electrons, and the subsequent energy transfer from elec- trons to ions is of picosecond order. Hence, lattice motion is negligible within the femtosecond pulse duration, whereas femtose- cond photon-electron interactions dominate the entire fabrication process. Therefore, femtosecond laser fabrication must be improved by controlling localized transient electron dynamics, which poses a challenge for measuring and controlling at the elec- tron level during fabrication processes. Pump-probe spectroscopy presents a viable solution, which can be used to observe elec- tron dynamics during a chemical reaction. In fact, femtosecond pulse durations are shorter than many physical/chemical characteristic times, which permits manipulating, adjusting, or interfering with electron dynamics. Hence, we proposed to control localized transient electron dynamics by temporally or spatially shaping femtosecond pulses, and further to modify localized tran- sient materials properties, and then to adjust material phase change, and eventually to implement a novel fabrication method. This review covers our progresses over the past decade regarding electrons dynamics control (EDC) by shaping femtosecond laser pulses in micro/nanomanufacturing: (1) Theoretical models were developed to prove EDC feasibility and reveal its mechanisms; (2) on the basis of the theoretical predictions, many experiments are conducted to validate our EDC-based femtosecond laser fabrication method. Seven examples are reported, which proves that the proposed method can significantly improve fabrication precision, quality, throughput and repeatability and effectively control micro/nanoscale structures; (3) a multiscale measurement system was proposed and developed to study the fundamentals of EDC from the femtosecond scale to the nanosecond scale and to the millisecond scale; and (4) As an example of practical applications, our method was employed to fabricate some key struc- tures in one of the 16 Chinese National S&T Major Projects, for which electron dynamics were measured using our multiscale measurement system. Light: Science & Applications (2018) 7, 17134; doi:10.1038/lsa.2017.134; published online 9 February 2018 Keywords: electrons dynamics control; femtosecond laser; micro/nano fabrication; pulse shaping INTRODUCTION the initial defects of the target materials. Hence, femtosecond laser Because of their ultrashort irradiation periods and ultrahigh inten- ablation is deterministic and reproducible , and almost any material 8–10 sities, femtosecond laser pulses in some aspects fundamentally change can be machined using a femtosecond laser, including metals , 11,12 13,14 15 the laser-material interactions mechanisms compared with long laser semiconductors , dielectrics , polymers , two-dimensional 16–20 21 22,23 pulses, which has created wide-range and exciting new possibilities in materials , ultrahard materials and biological tissues .In 1–3 micro/nanoscale fabrication . addition, ionization mechanism can be adjusted by changing femto- The ultrahigh intensity makes femtosecond laser-material interac- second laser energy and its temporal/spatial distribution to control 4,5 24 tions a strongly nonlinear process . Linear ionization generally laser-material interactions . dominates in free electron generations for continuous or long-pulse The ultrashort irradiation period of a femtosecond laser also makes laser processing of wide bandgap materials . The intensities of femtosecond laser-material interactions a strongly nonequilibrium 12 − 2 3,25 femtosecond lasers can easily exceed 10 Wcm , thus nonlinear process . The duration of a femtosecond laser pulse is much shorter 12 − 2 − 10 − 12 ionization mechanism such as avalanche ionization (~10 Wcm ), than the electron-lattice energy relaxation time (10 − 10 s). 13 14 − 2 multiphoton ionization (~10 –10 Wcm ), and tunnel ionization Therefore, laser energy absorption is completed before the lattice 15 − 2 (410 Wcm ; Ref. 3), can occur in the femtosecond laser fabrica- changes, resulting in a significantly nonequilibrium state between tion processes. The nonlinear ionizations are almost independent on electrons and lattices. Hydrodynamic motion and heat conduction 1 2 Laser Micro/Nano-Fabrication Laboratory, School of Mechanical Engineering, Beijing Institute of Technology, Beijing 100081, China; Department of Mechanical Engineering, University of Minnesota, Minneapolis, MN 55455, USA and Department of Electrical Engineering, University of Nebraska-Lincoln, Lincoln, NE 68588-0511, USA Correspondence: L Jiang, Email: [email protected] Received 12 April 2017; revised 27 August 2017; accepted 28 August 2017; accepted article preview online 30 August 2017 The accepted article preview was available with the details: Light: Science & Applications (2018) 7, e17134; doi: 10.1038/lsa.2017.134. Electrons dynamics control LJiang et al Table 1 The estimations of optical and thermal properties in classical and improved two-temperature models Property Classical TTM Improved TTM ∂hi ε Electron heat capacity C = γT C ¼ n e e e e ∂T e pffiffiffiffiffi 3=2 3 m ðk T Þ e B pffiffiffi t ¼ 3mek e 2 2 2pðZ Þ n e lnL Electron relaxation time t ¼ e e 2 p n k T e B e fg 1 þ exp½ mðÞ n ; T =k T F e B 1=2 T 1 2 Electron heat conductivity k ¼ k k ¼ n t C e eq e e e T 3 e 2 2 2 2 o ðÞ n t o ðÞ n t e e p e p e εðÞ t; r ; z ¼ 1 þ i DR DT e 2 2 2 Optical property D 1þo t o 1þo t e ðÞ e ðÞ DR ðÞ DT max e max þ DεðÞ contribution of the interband transition ab c 14 000 14 000 Electron temperature Electron temperature Experiments Lattice temperature Lattice temperature 12 000 12 000 Existing model 0.8 Proposed model 10 000 10 000 0.6 8000 8000 6000 6000 0.4 4000 4000 0.2 2000 2000 0 0 0 2 3 4 5 –5 0 5 1015202530 –5 0 5 10 15 20 25 30 10 10 10 10 Pulse duration (fs) Time (ps) Time (ps) − 2 Figure 1 The calculated electron and phonon temperatures of 200 nm gold film irradiated by a 140 fs, 1053 nm pulse at 0.2 J cm by (a) the classical model and (b) the improved model. (c) The predicted damage threshold fluences of 200 nm gold film processed by a 1053 nm laser at different pulse. Reproduced from Ref. 117 (with the permission of SPIE). through lattices are negligible during the femtosecond pulse duration. from experiments may contribute to theoretical comprehension and Thus, recast, thermal damage (microcracks) and heat-affected-zone prediction. are greatly reduced . Because of the nonequilibrium laser-material Many innovations in measurement technologies could also be interactions, including phase change and material removal, are converted to new methods in fabrication technologies. In the pump- 26,27 essentially determined by laser-electron interactions . Hence, fem- probe experiments, the observation of electron dynamics is achieved by analyzing the probe pulse disturbed by electrons. By increasing the tosecond laser fabrication must be improved by controlling electron energy of the probe laser, the disturbing process can be transformed dynamics during the fabrication process, which poses a challenge for measuring and controlling at the electron level during fabrication into certain control if it is purposely designed. Much more complex processes. pulses other than merely two pulses can be designed to control the Pump-probe spectroscopy, which has been widely used in fields electron dynamics. Recent developments of optical devices substan- 28–33 34–37 38–40 58–63 including physics ,chemistry and materials science , tially enhanced the capability to shape laser pulses. The amplitude 64,65 66,67 68,69 presents a viable solution for detecting electron dynamics. The basic /phase /polarization /frequency can be easily manipulated in mechanism involves splitting one femtosecond pulse into two sub- both temporal and spatial domains. By temporally/spatially shaping pulses. One sub-pulse is the pump pulse and the other is the probe the femtosecond laser, the local transient electron dynamics can be pulse, and the time delay between them can be precisely controlled by precisely controlled. Many studies have demonstrated that properties adjusting the optical path difference. On the basis of this principle, of atoms and molecules can be controlled by temporal/spatially shaped 41 70–72 Zewail and other pioneers developed the discipline of femtochem- ultrafast laser pulses . For example, using a shaped ultrafast laser 73–75 istry, which elucidates ultrafast chemical reaction processes at the pulse, atoms can be selectively ionized ; spin states can be 42–44 76–80 electron level . Their research demonstrates how chemical bonds dynamically controlled ; molecular ground state rotational 45,46 81,82 change and electrons transfer during chemical reactions . Recently, dynamics and vibrational modes can be manipulated ;chemical 83,84 attosecond laser pulses have been used to probe electron dynamics in reactions can be controlled and X-ray line emission from plasma 47–49 85–87 greater detail , and even the electronic charge distribution within under the femtosecond pulse can be significantly enhanced .In the 34,50,51 5,88–90 molecules has been detected . In addition, electron dynamics field of ultrafast laser fabrication, Herman and Ilday et al during laser-material interactions have been captured through the demonstrated the heat accumulation effect can be controlled by pump-probe technique. For instance, researchers have imaged electron harnessing high-repetition-rate burst trains of ultrafast laser pulses, excitation, electron spatiotemporal distribution and electron decay on achieving significant benefits compared with unshaped pulses. Xu 52–54 24,91,92 the femtosecond scale . Furthermore, the electron collision time et al revealed the conceptual significance of controlling coherent and laser-induced plasma lifetime can be determined using pump- phonon dynamics for controlling phase change by specifically designed 55–57 3 probe measurement . The information on electrons measured ultrafast pulse trains with different pulse delays. Gamaly reported Light: Science & Applications doi:10.1038/lsa.2017.134 Temperature (K) Temperature (K) –2 Threshold fluence (J cm ) Electrons dynamics control LJiang et al a c 0.2 0.2 0.1 0.1 0.0 0.0 –0.1 –0.1 –0.2 0.8 –0.2 010 200 10 20 0.6 Time (fs) Time (fs) 0.4 0.2 0.0 0.2 0.1 0.0 –1.1e–2 1.5e–2 –1.1e–2 1.5e–2 –9e–4 9e–4 –0.1 Two pulses 1:1 Two pulses 4:1 Difference –0.2 400nm+800nm 800nm+400nm 800nm+800nm 0 510010 20 30 0 25 50 75 100 125025 50 75 100 1250 25 50 75 100 125 0 20 40 Time (fs) Time (fs) Time (fs) Time (fs) Time (fs) Time (fs) Figure 2 Schemes of electrons dynamics adjusted by shaped femtosecond laser pulses. Electric fields of the applied laser pulse and time-dependent excited 99 99 electrons of diamond with different (a) pulse delays and (b) sub-pulse numbers .(c)Electric fields of the applied laser pulse (top panel) and the electron density change of diamond from that in the ground state after laser termination with different pulse energy ratios .(d) Time-dependent electron temperature of fused silica with different pulse dual-wavelengths . Reproduced from Ref. 99 (with the permission of IOP) and Ref. 111 (with the permission of AIP publishing). controlling over the ablation rate and phase state of laser-produced domains, we are able to control photon–electron interactions, and plume by using temporal/sptial pulse shaping. Sheppard and Wilson then to control the localized transient electron dynamics (including produced the Bessel beam with an annular lens, and Marcinkevičius electron density, temperature and excited-state distribution), and 94 95–97 et al and Courvoisier et al obtained some micro/nanostructures further to modify localized transient materials properties, and then by spatially shaping ultrafast laser pulses. to adjust material phase change, and finally to implement the novel Although various results have been obtained to control the laser- fabrication method. We have devoted our efforts on using EDC to material interactions, electron dynamics have not been extensively improve femtosecond micro/nano fabrications for more than a studied or deliberately controlled, which is vital for laser micro/nano decade. This review summarizes our recent progresses based on fabrication. In femtosecond laser micro/nano fabrication process, EDC by shaping femtosecond pulses in micro/nanofabrication: photons are mainly absorbed by electrons, and the subsequent energy transfer from electrons to ions is of picosecond order 1) Theoretical modeling for EDC feasibility: four models were − 12 − 10 (10 − 10 s). Hence, lattice motion is negligible within the developed, which consists of the ab initio calculations for electron femtosecond pulse duration, and femtosecond photon − electron 98–107 dynamics , a revised molecular dynamics simulation for phase interactions are the only factor to be considered, which dominates 108,109 110–115 change , a plasma model for ionization processes and − 10 − 3 the subsequent fabrication processes (10 –10 s). Therefore, 116–118 an improved two-temperature model for energy transport . femtosecond laser fabrication must be improved by controlling 2) Experiments for EDC validation: By shaping a femtosecond pulse localized transient electron dynamics. The key challenge is to measure in temporal or spatial domains, the photon-electron interactions and control at the electron level during fabrication processes. Pump- can be controlled to adjust the localized transient electron probe spectroscopy presents a viable solution, which can be used to dynamics to implement the novel fabrication methods. Various observe electron dynamics during a chemical reaction. In the pump- experiments were conducted to validate the effectiveness of EDC. probe experiments, the observation of electron dynamics is achieved Basedon EDC, weproposedto: by analyzing the probe pulse disturbed by electrons. By increasing the energy of the probe laser, the disturbing process can be transformed ▪ control the localized transient electron density to induce reso- into certain control if it is purposely designed. Therefore, we propose nance absorption, by which microchannel processing efficiency the core idea of electrons dynamics control (EDC; Supplementary was increased by 56 times and the maximum aspect ratio was Movie 1, Supplementary Information): By shaping the amplitude, 119,120 extended by three times ; phase and polarization of femtosecond pulses in temporal and spatial doi:10.1038/lsa.2017.134 Light: Science & Applications Excitation energy Electric amplitude (a.u.) Electric amplitude (a.u.) Electric amplitude (a.u.) (eV per atom) Electric temperature (10 K) Electric amplitude (a.u.) 50 –25 –50 –25 –50 –25 –50 –25 –50 –25 –50 –25 –50 –25 –50 –25 –50 –25 –50 –25 –50 Depth (nm) Depth (nm) Depth (nm) Time (fs) Depth (nm) Time (fs) Time (fs) Time (fs) Time (fs) Depth (nm) Electrons dynamics control LJiang et al a b c d e Pulse delay: 0 fs Pulse delay: 100 fs Pulse delay: 200 fs Pulse delay: 300 fs Pulse delay: 400 fs 1.0 0.8 0.6 0.4 0.2 0.0 fg h i j Pulse delay: 0 fs Pulse delay: 100 fs Pulse delay: 200 fs Pulse delay: 300 fs Pulse delay: 400 fs 75.0 50.0 25.0 0.0 Figure 3 Simulation of materials properties adjusted by varying pulse delay within a femtosecond laser pulse train. (a–e)Reflectivity of the material surface and (f–j) the corresponding peak laser intensity distribution at different pulse delays. ▪ modify free electron density and the corresponding photon Furthermore, we revealed the forming mechanisms of the high- absorption efficiency, then modulate material properties, by quality and high aspect-ratio microholes using the multiscale time- which laser-assisted chemical etching rate was enhanced by 37 resolved measurement system . times ; ▪ adjust electron generation on fabricated material surfaces, through which the periods, orientations and structures of the MODELING FOR EDC FEASIBILITY 122,123 surface ripples can be effectively adjusted ; Modeling method ▪ control electron density and its distribution, through which we Femtosecond photon–electron interactions dominate the entire none- obtained controllable micro/nano hierarchical structures on quilibrium and nonlinear laser fabrication processes, including the material surfaces and enhancement factors of surface-enhanced absorption of laser energy by electrons, energy transfer from electrons 9 136–138 Raman scattering (SERS) up to 1.1 × 10 (Refs. 124–126); to lattices, plasma generation, phase change and material ▪ control electron density distribution using temporally shaped modification. These interactions range from nanometers to milli- femtosecond laser pulse to modify chemical and physical proper- meters spatially and from femtoseconds to microseconds temporally. ties of materials, by which polymorphous Au-MoS hybrids were Although a growing body of experimental observations exists, a prepared ; comprehensive model remains unavailable. Such a model is essential ▪ adjust the phase change mechanism by changing photon-electron for revealing the fundamental science underlying ultrafast laser- interactions, which reduced the recast layer thickness by 60% material interactions. According to the applicable time and space (Ref. 128); scales of each involved process, four fundamental models were used: ▪ manipulate electron density distribution using spatially modu- lated femtosecond laser pulse, by which deep subwavelength 1) The ab initio model based on the time-dependent density func- (B1/14 of the laser wavelength) and high conductivity (B1/4 of tional theory (TDDFT) was employed to understand initial non- the bulk gold) nanowires were fabricated in the open air . linear laser radiation absorption through photon-electron 98,99 interactions . 3) Measurements for EDC fundamentals: A multiscale measurement In TDDFT, the fundamental variable is not the many-body wave system (from femtosecond scale to millisecond scale) was proposed function in quantum mechanics, but the electronic density. This and developed for understanding the electron dynamics during time-dependent electron density is determined by solving an femtosecond laser ablation. It comprises a pump-probe shadow- auxiliary set of noninteracting Kohn-Sham equations, where the 130,131 graph imaging technique , time-resolved plasma laser fields are treated as time-dependent spatial uniform external 132,133 photography , laser-induced breakdown spectroscopy vector potentials. Describing the motion of the electrons in the (LIBS) and commercialized fast imaging device (CCD). We system, the time-dependent Kohn-Sham equation for single reveal the multiple time scale fundamentals during femtosecond particle orbitals is laser-material interactions, including the femtosecond-scale propa- gation of a laser pulse, picosecond-scale generation/evolution of ! ! ! laser-induced plasma, nanosecond-scale plasma ejection/expansion, i_ c r ; t ¼ H r ; t c r ; t ð1Þ KS ∂t and millisecond-scale hole formation. 4) Applications for the EDC method: Our proposed method was used to fabricate some key structures. Using spatially-shaped femtose- ! ! cond pulses, we optimized electron density distribution in plasma nr ; t ¼ c r ; t ð2Þ at the focus point, and then manipulated plasma expansion and phase change, which was used to drill microholes with a diameter ! ! of 1.6 μm and an aspect ratio of 330:1 (Refs. 132,135). wherenr ; t is the electron density and H r ; t is the Kohn- KS Light: Science & Applications doi:10.1038/lsa.2017.134 1.0 1.0 0.8 1.0 0.8 1.0 0.6 0.8 1.0 0.6 0.8 0.4 0.6 0.8 0.4 0.6 0.2 0.4 0.2 0.6 0 0.4 0.2 0.4 0.2 0.2 80 120 60 100 40 80 Radius (μm) Radius (μm) Radius (μm) Radius (μm) Radius (μm) Radius (μm) Radius (μm) Radius (μm) Radius (μm) Radius (μm) Laser intenstiy Reflectivity Laser intenstiy Reflectivity Laser intenstiy Reflectivity Laser intenstiy Reflectivity Laser intenstiy Reflectivity Electrons dynamics control LJiang et al –2 –2 Laser fluence=0.28 J cm 1pulse 15ps Laser fluence=0.28 J cm 20pulse 15ps ae 80 80 40 40 0 0 –40 –40 –80 –80 700 350 0 –350 –700 –1050 –1400 –1750 –2100 –2450 –2800 700 350 0 –350 –700 –1050 –1400 –1750 –2100 –2450 –2800 Z (Å) Z (Å) –2 –2 Laser fluence=0.28 J cm 1pulse 50ps Laser fluence=0.28 J cm 20pulse 50ps b f 80 80 40 40 0 0 –40 –40 –80 –80 700 350 0 –350 –700 –1050 –1400 –1750 –2100 –2450 –2800 700 350 0 –350 –700 –1050 –1400 –1750 –2100 –2450 –2800 Z (Å) Z (Å) –2 –2 Laser fluence=0.28 J cm 1pulse 92ps Laser fluence=0.28 J cm 20pulse 92ps c g 80 80 40 40 0 0 –40 –40 –80 –80 700 350 0 –350 –700 –1050 –1400 –1750 –2100 –2450 –2800 700 350 0 –350 –700 –1050 –1400 –1750 –2100 –2450 –2800 Z (Å) Z (Å) –2 –2 Laser fluence=0.28 J cm 1pulse 200ps Laser fluence=0.28 J cm 20pulse 200ps d h 80 80 40 40 –40 –40 –80 –80 700 350 0 –350 –700 –1050 –1400 –1750 –2100 –2450 –2800 700 350 0 –350 –700 –1050 –1400 –1750 –2100 –2450 –2800 Z (Å) Z (Å) ik m n 1.0×10 4 4 1.2×10 1.2×10 3 15 15 ps 4 4 8.0×10 1.0×10 1.0×10 20 ps 50 ps 3 8.0×10 3 10 8.0×10 6.0×10 92 ps 15 ps 3 3 200 ps 20 ps 6.0×10 6.0×10 4.0×10 50 ps 92 ps 4.0×10 4.0×10 Binodal Binodal 3 200 ps 2.0×10 Spinodal Spinodal 3 3 2.0×10 2.0×10 Dynamic film surface z=–7 nm 0.0 –5 0.0 0.0 3 3 3 3 3 3 700 500 300 100 –100 –300 –500 –700 700 500 300 100 –100 –300 –500 –700 00 3.0×10 6.0×10 9.0×10 3.0×10 6.0×10 9.0×10 –3 –3 Density (kg m ) Density (kg m ) Z (Å) Z (Å) jl o p 4 20 1.0×10 4 4 1.2×10 1.2×10 3 4 15 4 8.0×10 15 ps 1.0×10 1.0×10 20 ps 3 3 50 ps 8.0×10 8.0×10 3 10 6.0×10 15 ps 92 ps 3 3 20 ps 200 ps 6.0×10 6.0×10 4.0×10 50 ps 3 3 92 ps 4.0×10 4.0×10 Binodal Binodal 200 ps 0 Spinodal 2.0×10 3 Spinodal 3 2.0×10 2.0×10 z=–7nm Dynamic film surface 0.0 –5 0.0 0.0 3 3 3 3 3 3 700 500 300 100 –100 –300 –500 –700 700 500 300 100 –100 –300 –500 –700 0 3.0×10 6.0×10 9.0×10 0 3.0×10 6.0×10 9.0×10 –3 –3 Z (Å) Z (Å) Density (kg m ) Density (kg m ) Figure 4 Schemes of phase change controlled by varying the pulse delay in a fs pulse train . Snapshots of nickel thin films irradiated by femtosecond laser − 2 (a–d) single pulse and (e–h) 20 pulse trains with the total fluence of 0.28 J cm ,where X is in the direction of Ni (100) surface and Z is in the direction of laser irradiance. Lattice temperature and stress distributions at different times for (i and k) the single pulse and (j and l) the 20 pulse trains. Time evolution of the system in the ρ-T plane for different regions for (m and n) the single pulse and (o and p) the 20 pulse trains. Arrows indicate the time evolution. Reproduced from Ref. 108 (with the permission of AIP publishing). Sham Hamiltonian given by 2) The molecular dynamics model was employed to reveal the phase change resulting from electron-ion interactions . In the mole- cular dynamics model, Newtonian equations of motion of a set of ! ! ! 1 e H r ; t ¼ p þ A ðÞ t þ V r ; t KS tot ion 2m c N particles are solved to describe the phase change process. ! ð3Þ n r ;t 2 0 þe dr þ V r ; t XC ! ! r r d r ∂U m ¼ f ; f ¼ ð5Þ i t t dt ∂t ! ! where V r ; t is the electron-ion potential and V r ; t is ion XC where r is the coordinate of the ith atom, and f is the force the exchange-correlation (XC) potential.The time evolution of the i i acting on the ith atom that is usually derived from potential wave function over a short period, Δt, can be approximately energy U(r) (Ref. 139): calculated as follows: hi q q ^ Dt ^ Dt 1 2 gðÞ 1r=r g=21ðÞ r=r 0 0 ! iH ðÞ tþDt iH ðÞ t ! KS KS 2 2 UrðÞ ¼ þ D e 2e ð6Þ c r ; t þ Dt ¼ e e c r ; t ð4Þ 4pεr doi:10.1038/lsa.2017.134 Light: Science & Applications Temperature (K) Temperature (K) X (Å) X (Å) X (Å) X (Å) Stress (GPa) Stress (GPa) Temperature (K) Temperature (K) X (Å) X (Å) X (Å) X (Å) Temperature (K) Temperature (K) Electrons dynamics control LJiang et al ac –2 5.0 J cm –2 4.0 J cm 75 fs 50 fs 25 fs AB CD 0 fs Energy ratio 1:2 –2 5.0 J cm –2 4.0 J cm 40 20 0 20 40 Radius (μm) 40 AB 50 CD Energy ratio 1:1 –2 5.0 J cm –2 4.0 J cm CD AB 200 80 400nm+800nm 800nm+400nm 800nm+800nm Energy ratio 2:1 1.0 0.5 0.0 0.5 1.0 1.0 0.5 0.0 0.5 1.0 Radius (μm) Radius (μm) 110–112 Figure 5 Schemes of ablation crater shape controlled by shaped femtosecond laser pulse .(a) Ablation crater shapes created by femtosecond pulse − 2 trains consisting of double pulses with different pulse delays at a total fluence of 5 J cm and central wavelength of 780 nm. (b) Ablation crater shapes − 2 created by femtosecond laser pulse trains with different wavelength composition at the total fluence of 5.0 J cm and the pulse delay of 50 fs. (c)Ablation crater shapes created by 800 nm femtosecond pulse trains consisting of double pulses with three different energy ratios at the pulse delay of 50 fs. Reproduced from Ref. 110 (with the permission of IOP), 111 (with the permission of AIP publishing) and 112 (with the permission of Springer). 3) We proposed the plasma model with quantum treatment to investigate plasma generation and changes (ionization and recom- 110,113–115 (μs~ms) bination) . 114,115 In the plasma model , the free electron density distribution in dielectrics under a femtosecond laser pulse is obtained by solving the Fokker − Planck Equation: ∂nðÞ t; r; z (ps~fs) e ¼ a ItðÞ ; r; z nðÞ t; r; z þ dðÞ ItðÞ ; r; z t e N ∂t nðÞ t; r; z ð7Þ where t is time, r is the distance to the Gaussian beam axis, z is the (ps~fs) depth from the surface of the bulk material, n (t,r,z) is the free electron density, τ is the free electron recombination time, α is the impact ionization constant, and δ is the cross-section of N- photon absorption. (ps~fs) The optical properties of ionized dielectrics are calculated using the Drude model. The spatial and temporal dependence of the complex dielectric function for the plasma is expressed as: 2 2 2 2 o ðÞ n t o ðÞ n t e e p e p e εðÞ t; r; z ¼ 1 þ i ð8Þ 2 2 2 2 1 þ o t o 1 þ o t e e pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi where o ðÞ n ¼ nðÞ t; r; z e =m ε is the plasma frequency and p e e e 0 τ is the free electron relaxation time, calculated as, follows by applying formulas derived from the Boltzmann transport Equation: Figure 6 Temporally shaping femtosecond pulses: (a) a conventional pffiffiffiffiffiffi 3=2 femtosecond pulse is temporally shaped into a pulse train; (b) the number of 3 mðÞ k TðÞ t; r; z e B e tðÞ t; r; z ¼ pffiffiffi sub-pulses within a train can be controlled; (c) the delay between sub-pulses e 2 2pðÞ Z nðÞ t; r; z e lnL ð9Þ can be controlled; (d) the energy ratio of sub-pulses can be controlled. ðÞ 1 þ expðÞ mðÞ n ; T =k TðÞ t; r; z F e e B e 1=2 where Z* is the ionization state, m is the electron mass, k is the e B Light: Science & Applications doi:10.1038/lsa.2017.134 Depth (nm) Depth (nm) Depth (nm) Depth (nm) Depth (nm) Electrons dynamics control LJiang et al CCD Reflector lens Grating Lens SLM Lens Grating Reflector Reflector Fs laser Polarizer Pulse shaper Shutter Dichroic mirror Half wave Attenuator plate Objective lens Sample Computer control Hexapod stage Halogen light Figure 7 Schematic of the experimental setup for temporal shaping of femtosecond pulses. Reproduced from Ref. 141(with the permission of Springer). a b c d D1: 185.6441 um 1 μm 1 μm 1 μm D1: 12.28527 um Figure 8 Holes machined using (a) five 355-nm nanosecond pulses (duration of 30 ns, energy of 0.16 μJ) ,(b) five 800-nm femtosecond pulses (duration 119 119 of 120 fs, energy of 0.12 μJ) ,(c) five femtosecond-nanosecond pulse pairs ,(d) an 800 nm femtosecond laser double-pulse train (duration of 50 fs, energy of 20 μJ) . Reproduced from Ref. 119 (with the permission of OSA) and Ref. 120 (with the permission of OSA). Boltzmann constant, F denotes the Fermi-Dirac integrals, μ is transport equation for dense plasmas; and (3) the free-electron heating 1/2 the chemical potential and ln Λ is the Coulomb logarithm.In and interband transition were taken into account using a modified addition, the wave-particle duality of photon is also considered to Drude model as Equations (8) and (9) to calculate the reflectivity and predict the formation mechanism of laser-induced periodic surface absorption coefficient. Table 1 compares the estimates of optical and structures . thermal properties between classical and improved two-temperature 4) We proposed the improved two-temperature model with full-run models. Figure 1 demonstrates the differences in (a, b) electron and quantum treatment to predict the electron and lattice temperature phonon temperatures between the classical approach and the 116,117 distributions . In the improved two-temperature model, improved model. The improved model significantly increases the electron and lattice temperatures are given by prediction precisions of the damage thresholds compared with the classical model, as shown in Figure 1c . ∂T C ðÞ T ¼ ∇ ½ kðÞ T ∇T GTðÞ T þ SzðÞ ; t ð10Þ e e e e e e l The four theoretical models are suitable for systems with different ∂t temporal and spatial scales and must be employed in combination to compensate for each other’s limitations. In our improved two- ∂T CðÞ T ¼ GTðÞ T ð11Þ l l e l temperature model, the electron density is assumed to be constant. ∂t Consequently, the model is not suitable when the bound electrons in where T is the electron temperature, T is the lattice temperature, e l non-metals are substantially ionized by a femtosecond laser pulse. The C is the electron heat capacity, C is the lattice heat capacity, k is e l e plasma model is proposed to consider the ionization process in the electron heat conductivity, G is the electron-lattice coupling nonmetals, but material ablation is assumed to commence when the factor, and S is the laser source term. free electron density reaches the critical density without the con- sideration of phase changes, for example, melting, vaporization, The two-temperature model was improved by us as follows :(1) Coulomb explosion and electrostatic repulsion. Moreover, the mole- using the Fermi distribution, the heat capacity of the free electrons was cular dynamics can be used to identify the relative roles of different calculated; (2) free-electron relaxation time and electron conductivity were determined using a quantum model derived from the Boltzmann phase change mechanisms. However, the interatomic potential doi:10.1038/lsa.2017.134 Light: Science & Applications Electrons dynamics control LJiang et al 0 0 (ns~ms) 0 fs 50 fs a b cd e k l –3 –100 –100 n (cm ) 2.000E+21 –200 –200 2 μm 2 μm 5 μm 10 μm Single pulse 20 μm t = 0 t = 20 t = 45 t = 90 t = 150 e e e e e –300 –300 1.500E+21 (fs~ps) f gh ij –400 –400 –1.5 –1.0 –0.5 0.0 0.5 1.0 1.5 –1.5 –1.0 –0.5 0.0 0.5 1.0 1.5 1.000E+21 0 0 100 fs 350 fs Double pulse m n 2 μm 2 μm 5 μm 10 μm 20 μm t = 0 t = 20 t = 45 t = 90 t = 150 –100 –100 =350 fs e e e e e 5.000E+20 –200 –200 0.000 –300 –300 –400 –400 –1.5 –1.0 –0.5 0.0 0.5 1.0 1.5 –1.5 –1.0 –0.5 0.0 0.5 1.0 1.5 Radius (μm) op 20 20 Original beam Original beam Transmitted beam Transmitted beam O 16 16 Si 12 12 8 8 t u 4 4 D D 1 2 3.0 0 0 0 100 200 300 400 500 600 0 100 200 300 400 500 600 Time (fs) Time (fs) 2.6 qr 20 20 Original beam Original beam 2.2 Transmitted beam Double pulse 16 Transmitted beam 16 12 12 1.8 8 8 Single pulse 1.4 200 400 600 800 1000 1200 0 500 1000 1500 2000 0 0 –1 0 100 200 300 400 500 600 0 100 200 300 400 500 600 Wave number (cm ) Pulse delay (fs) Time (fs) Time (fs) Figure 9 Morphology evolution of the sample surface exposed by (a–e) a femtosecond laser single pulse; (f–j) femtosecond laser double pulses at different stages of the etching process, the pulse delay is 350 fs. The SEM images have varying scale bars. (k–n) Simulation of the free electron density distributions and (o–r) center laser intensity distributions in fused silica irradiated by femtosecond double pulses at different pulse delays. (s) Schematic diagram of the manufacturing and Si-O bond structure. (t) Normalized Raman spectra of modified regions irradiated using femtosecond laser single and double pulses (the pulse delay is 350 fs) in fused silica. Dashed lines below the D2 peaks are baselines used in the peak area measurement in u. Inset is the schematic diagram of 4- and 3-membered ring structures. (u) Percent area of the total reduced Raman spectrum under the D2 line versus different pulse delays. The − 2 femtosecond laser with wavelength of 800 nm, duaration of 50 fs and repetition rate up to 1 KHz. The laser fluence is fixed at 9.46 J cm in all experiments and the energy distribution ratio is 1:1. Reproduced from Ref. 121 (with the permission of NPG). remains constant in the molecular dynamics model, even though it in validated the concept of EDC and predicted its potentials in ultrafast 88–107 laser micro/nanofabrication . fact changes dramatically during strong ionizations. The ab initio method is therefore used to describe a system with no parameteriza- Simulation results tion, which facilitates the investigation of processes whose mechanisms In this section, we report several representative simulation results. As are not fully understood. In a typical ab initio method, the potential shown in Figure 2, the time-dependent electrons dynamics in energy for a particular electronic state is defined by the electronic materials irradiated by temporally shaped femtosecond laser pulse Schrödinger equation. However, the computational demand of ab train have been calculated based on TDDFT and the improved two- initio calculation is substantial, and impact ionization is not considered temperature model . The simulation results indicated that increasing because the theory is based on orbitals that interact only through the pulse delay from 10 to 30 fs markedly increased the number of excited mean field. The exact time-dependent xc potential and the functional electrons and absorbed energy (Figure 2a). However, the excitation for physical observables are still being studied. energy decreased as the sub-pulse number per train increased from 1 On the basis of the aforementioned factors, during femtosecond to 4 (Figure 2b). Furthermore, the shaped pulse energy and laser-material interactions, the target material can be conceptually wavelength distribution were also determined to be key parameters divided into three systems: the electron system, atom or molecule for changing the electron density and temperature (Figure 2c and 2d). system and plasma system outside the bulk material. The ab initio The results indicated that localized transient electrons dynamics, model can be applied to sub-nanometers sampling areas to determine including the number of excited electrons, excitation energy, electron the ionization mechanisms and multiphoton absorption cross density distribution, and electron temperature, can be controlled by section . The revised molecular dynamics can be employed to shaping femtosecond pulses (pulse delay, sub-pulse number, pulse determine the phase change mechanisms in areas of tens to hundreds energy ratio and dual-wavelength distribution). of nanometers within an electric field formed by the plasma system We theoretically demonstrated that localized transient materials outside the bulk material, which is calculated using our plasma properties can be effectively modulated through localized transient model . The plasma model can be used to describe the electron EDC by temporally shaping femtosecond laser pulses (Figure 3). The system for photon absorption and plasma generation (ionization) and reflectivity and corresponding peak laser intensity distribution were 110,114 recombination before the lattice actually changes .Our improved completely changed by varying the pulse delays. Moreover, the two-temperature model can be employed to calculate energy transport electron density was adjusted by varying the pulse delay, and 116,117 through electron-phonon interactions .Using themodels,we reflectivity decreased as the pulse delay increased because of the free Light: Science & Applications doi:10.1038/lsa.2017.134 Intensity (arbitrary units) Percent area of the reduced spectrum –1 under D =605 cm line Depth (nm) 13 –2 13 –2 Laser intensity (10 W cm ) Laser intensity (10 W cm ) 13 –2 13 –2 Laser intensity (10 W cm ) Laser intensity (10 W cm ) Electrons dynamics control LJiang et al 1.05 F 1.16 F th th Δt = 0 fs Δt = 50 fs Δt = 0 fs Δt = 100 fs Δt = 300 fs a b c eg 2 μm2 μm 1 μm 4 μm 1 μm 1 μm Δt = 50 fs Δt = 200 fs Δt = 400 fs df h 1 μm 1 μm 1 μm 1 μm 4 μm 1 μm 1.8F 2F 0.6F 1.2F ij k th l th th th 2 μm 2 μm 2 μm 2 μm AFM profile 2F m 2F n 2F o th th th –0.3 –0.4 LSFL(//E ) HSFL(⊥E ) –0.5 –0.6 –0.7 0.0 0.5 1.0 1.5 2.0 2 μm 2 μm 2 μm Distance (μm) Figure 10 Control of period, orientation and topology of the LIPSS on the surface of fused silica via symmetrically shaped femtosecond pulses. (a–h)LSFL changes to HSFL with different orientation under certain pulse fluences and pulse delays on fused silica .(i–o) SEM images of representative HSFL (i–k), LSFL (m–o) and double-grating structure (l) on fused silica for single-pulse trains (Ns= 1), double-pulse trains (Ns= 2), triple-pulse trains (Ns= 3) and quadruple-pulse trains (Ns= 4), respectively . Reproduced from Ref. 122 (with the permission of OSA) and Ref. 123 (with the permission of OSA). electron density was lower. Thus, the original laser intensity distribu- slightly above the volumetric phase change threshold fluence of single tion was considerably reshaped, providing substantial control over the pulse ablation. Therefore, the results demonstrated that the phase energy absorption process. The results demonstrated that by control- change mechanism can be controlled by carefully designing the pulse ling localized transient electron dynamics through shaping femtose- train parameters. cond pulses, the localized transient material properties, such as the The ablation shape modulated by designed femtosecond pulse trains 110–112 on fused silica was validated using our plasma model . The crater reflectivity and absorption coefficient, can be modified. Meanwhile, we studied the phase change mechanisms in materials shapes of the fused silica and the spacing of the subwavelength ripples such as nickel thin films under irradiation by femtosecond laser pulse differed greatly among sub-pulses with different delay times, wave- trains by using molecular dynamics simulations and the two- lengths and energy distributions. The ablation depths were approxi- temperature model . The theoretical simulation revealed that more mately 170, 145, 102 and 91 nm for pulse delays of 0, 25, 50 and 75 fs, and smaller high-quality uniform nanoparticles can be obtained on the respectively (Figure 5a). The ablation depths were approximately 186, nickel thin film by femtosecond laser pulse trains compared with 43 and 139 nm, and the ablation radii were approximately 860, 630 conventional pulses (Figure 4a–4h). The use of pulse train reduced the and 410 nm for pulse trains of 400+800, 800+400 and 800+800 nm, − 2 electron temperature and electron thermal conductivity dramatically respectively (Figure 5b). For 4.0 J cm , the depths of the microholes because of the lower intensity of the sub-pulses or higher transient were 77, 84 and 71 nm (depth from the surface to the position of line surface temperatures, which left the absorbed energies deposited AB), and the corresponding radii of the ablation craters were 250, 160 mainly within the nanoscale layers of the dynamic film surface and 270 nm for the energy ratios of 1:2, 1:1 and 2:1, respectively; For − 2 (Figure 4i and 4j). By designing the pulse train, smaller film 5.0 J cm , the depths of the microholes were 80, 95 and 70 nm compressive stresses and tensile stresses can be obtained, which (depth from the surface to the position of line CD), and the reduced microcracks (Figure 4k and 4l). Furthermore, a transition corresponding radii of the ablation craters were 430, 410 and from phase explosion to critical point phase separation (Figure 4m– 440 nm for the energy ratios of 1:2, 1:1 and 2:1, respectively 4p) enabled small uniform nanoparticle generation. The modulated (Figure 5c). Calculations based on the plasma model revealed that phase change mechanism was considered the main factor for the the time-dependent free electron density differed substantially among morphology control. Both the compressive and tensile stresses can be laser pulses with different parameters. Changing the ionized electron reduced by the pulse trains, leading to the critical point phase density distributions significantly modified the optical and thermal separation within the uppermost films and no liquid-vapor phase properties of the material. This interaction process greatly altered separation within the subsurface films when total laser fluence was ablation shape and subwavelength ripple. The results indicated that by doi:10.1038/lsa.2017.134 Light: Science & Applications Ns=1 Ns=1 Ns=2 Ns=2 Ns=3 Ns=3 Ns=4 Depth (μm) Electrons dynamics control LJiang et al 1.0 μm 1.0 μm a b f Δt Δt ~ ms Y 2nd 1st 2nd 1st hv + + – Ag +H O Ag+H +OH Ripple Fs laser Fs laser pulse train Nanopillar array AgNO Z Water 3 Y Si Si Y STEP TWO STEP ONE Reference 0 fs Δt =0 fs 400 fs Δt =800 fs 600 fs 800 fs 500 nm 500 nm 18 1000 fs ik 60 8 0 fs 600 fs 800 fs 1000 fs Reference 0 40 80 120 160 180 240 600 1000 1400 1800 –1 Partical diameter (nm) Raman shift (cm ) 8 500 nm 500 nm h 6 –1 Raman shift (cm ) 0 fs 400 fs 600 fs 800 fs 1000 fs 610 9.2E+06 6.2E+07 4.5E+07 2.2E+07 1.9E+07 1310 7.1E+06 2.7E+07 1.9E+07 8.4E+07 1.2E+07 1360 1.1E+07 4.8E+07 3.3E+07 1.6E+08 2.0E+07 600 800 1000 1200 1400 1600 1800 1510 1.1E+07 5.2E+07 3.5E+07 1.6E+08 2.3E+07 –1 Raman shift (cm ) 1650 8.0E+06 4.7E+07 3.0E+07 1.5E+08 1.8E+07 Figure 11 (a–e) One-step method fabrication of controllable SERS substrates .(a and b) SEM images of the silicon irradiated at pulse delays of 0 fs, and 800 fs, respectively. (c) Size distribution of silver nanoparticles in a (black), and b (red). (d) SERS signals of R6G molecules on the as-prepared substrates at various pulse delay. (e) Enhancement factors with different pulse delays. (f–k) Two-step method fabrication of controllable SERS substrate . f Schematic diagram of the SERRS substrate fabrication process. (g–j) SEM images of silicon substrates irradiated at pulse delays of g 0 fs, and i 1000 fs in deionized water, (h) 0 fs and (j) 1000 fs in 10-mM silver nitrate solution. (k) SERS spectrum of substrates fabricated at different pulse delays. Reproduced from Ref. 125 (with the permission of OSA), and Ref. 126 (with the permission of OSA). controlling femtosecond pulses trains, the final material modification us to further modify the chemical and physical properties of the can be improved. materials . By spatially shaping a femtosecond laser, a metal The model simulations validate theoretical feasibility of EDC, nanowire with a super sub-diffraction-limit precision (1/14th of the providing a theoretical prediction and guide for practical experiments wavelength) can be achieved . and applications. Furthermore, our models and subsequent multiscale measurements can form a mutually-supporting system. According to Femtosecond laser temporal pulse shaping the theoretical prediction, a novel fabrication method based on During femtosecond laser fabrication, most of the photon energy is temporally and spatially shaped pulse trains was proposed. Localized initially absorbed by electrons. However, conventional femtose- cond laser pulses are separated by a time scale that ranges from transient electron dynamics and the corresponding material properties microseconds to milliseconds (Figure 6), which is much longer can be actively controlled based on the proposed method. We discuss than the time scale of electron-lattice coupling (typically a few several pieces of experimental evidences in detail in the following sections. picoseconds to tens of picoseconds). Temporal pulse shaping enables sub-pulse generation with a pulse delay shorter than the EXPERIMENTS: EDC-BASED NOVEL FABRICATION METHOD characteristic time scale of electron-lattice coupling so that we can To validate the proposed method of EDC-based fabrication method control femtosecond laser photon–electron–phonon interactions. through temporal and spatial shaping of femtosecond laser pulses, As Figure 6a shows, a conventional femtosecond laser pulse can be 119–128 many experiments have been conducted . By designing tempo- split into several sub-pulses with a delay, which is called the pulse rally shaped pulse trains, we can control the localized transient train. The separation between the sub-pulses occurs in the time electron density to induce resonance absorption; thus, high efficiency scale—tens of femtoseconds to several picoseconds—which is 119,120 fabrication can be achieved . We can adjust the phase change generally similar to the characteristic time scale of electron- mechanism from thermal phase change to nonthermal phase change; lattice coupling. Before the photon energy that is absorbed by subsequently, high-quality fabrication can be performed .In addi- the electrons transfers to lattice, the subsequent sub-pulses con- tion, we can modify the free electron density, the corresponding tinuously interact with the materials so that controlling ultrafast photon-absorption efficiency and the material properties, enabling us photon-electron-phonon interactions is possible. Moreover, a to further control the chemical reaction . We can adjust the electron conventional femtosecond laser single pulse can be shaped into density and its distribution so that the periods, orientations and almost any arbitrary pulse shapes; for example, (1) a pulse can be 122,123,140 structures of the surface ripples can be effectively modulated split into a pulse train with different numbers of sub-pulses and a high sensitivity of SERS substrate (controllable micro/nano (Figure 6b); (2) the delay between the sub-pulses can be controlled 124–125 hierarchical structures on materials’ surfaces) can be achieved . (Figure 6c); and (3) the energy ratio of the sub-pulses can be Furthermore, we can control the electron density distribution controlled (Figure 6d). By shaping a femtosecond laser pulse in to induce surface chemical reduction activity in materials, enabling temporal domains to obtain a specific pulse shape, we can control Light: Science & Applications doi:10.1038/lsa.2017.134 Number Intensity (×10 a.u.) Intensity (×10 a.u.) Electrons dynamics control LJiang et al ab c d (ns~ms) 500 nm Single pulse 500 nm 500 nm 500 nm t = 0 t = 20 s t = 1 min t = 30 min r r r r ef g h (fs~ps) Double pulses ( = 5ps) 500 nm 500 nm 500 nm 500 nm t = 0 t = 20 s t = 1 min t = 30 min r r r r ij k MoS modified by temporally MoS modified by conventional 2 S 2p 2 1/2 S 2p 1/2 shaped femtosecond pulses femtosecond pulses Mo/S=1:4.28 Mo/S=1:2.85 Unbound Unbound sulfur sulfur (27.78%) (27.74%) S 2p 3/2 S 2p 3/2 4+ 4+ MoS film 172 171 170 169 168 167 166 165 164 163 162 161 160 159 158 172 171 170 169 168 167 166 165 164 163 162 161 160 159 158 Binding energy (eV) Binding energy (eV) l m n 3+ Au Shape-controlled Dangling bond Au NPs Unsaturated Saturated S atoms 500 nm Fs laser treated MoS S atoms Mo S 2 Figure 12 Morphology evolution of gold NPs reduced on MoS surfaces irradiated by femtosecond (a–d) single and (e–h) double pulses, at different stages of the reduction process, where tr represents the chemical reduction time, the pulse delay is 5 ps. (i) Schematic diagram of the manufacturing and Mo-S bond structure. XPS S 2p spectra of modified regions irradiated by femtosecond laser (j) double and (k) single pulses on MoS , where the percentage value represents the content of unbound sulfur and the atomic ratio represents the relative atomic concentration ratio of Mo and S atoms. (l) AFM image and (m) atomic scale schematic of the laser-broken micro/nano MoS debris. (n) Mechanism of chemical reduction of gold cations on laser-treated MoS 2 2 (Ref. 127). Reproduced from Ref. 127 (with the permission of ACS). the photon-electron interactions, which allow us to control the Temporal shaping of the femtosecond laser pulses can be achieved localized transient electron dynamics, including electron density, by using a commercial 4f-configuration-based pulse shaper (BSI temperature, excited state distribution and further modify loca- MIIPS BOX 640, Biophotonic Solutions, Inc., East Lansing, MI, lized transient material properties, adjust material phase changes USA), which allows us to split each conventional single pulse into a and ultimately implement the novel fabrication methods. pulse train and control the number of sub-pulses, delay between sub- pulses and energy ratio of the sub-pulses within a pulse train. The Experimental setup for temporal pulse shaping.The schematic of irradiation time (that is, number of pulse bursts) is precisely controlled the experimental setup for temporal shaping of femtosecond laser by using an electromechanical shutter. The sample is mounted on a pulses is shown in Figure 7 (Ref. 141). An amplified Ti: sapphire laser computer-controlled, six-axis translation stage (M-840.5DG, PI, Inc., system (Spectra Physics Inc., Santa Clara, CA, USA) is used to generate Karlsruhe, Germany) with a positioning resolution of 1 μm. The entire 35 fs (full width at half maximum) linearly polarized laser pulses on a fabrication process can be observed by using a charge-coupled device central wavelength of 800 nm with a repetition rate of 1 kHz. In (CCD) camera along with a white-light source irradiate on the sample addition, the pulse energy can be up to 3.5 W and can be continuously surface. adjusted by combining a half-wave plate with a polarizer. Pulse energy can also be reduced to the desired values according to specific High-efficiency fabrication by temporal pulse shaping based on experimental conditions by using a neutral density (ND) filters. EDC. By designing femtosecond laser pulse trains on the basis of doi:10.1038/lsa.2017.134 Light: Science & Applications Intensity (a.u.) Intensity (a.u.) Electrons dynamics control LJiang et al 0.40 WS CCD 0.36 BS 0.32 SLM 0.28 DM 0.24 0.20 Phase pattern 0 200 400 600 800 Lens Pulse delay (fs) bc 0 3 μm μm –200 –200 Translation stage Figure 15 Schematic diagram of the experimental setup. WS: white light nm nm source; BS: beam splitter; DM: dichroic mirror; L1, L2: two convex lenses consisting of a 4f relay system; L3: convex lens; Inset: section of the de samples and the focusing laser in the yz plane . Reproduced from Ref. 129 (with the permission of Wiley). Figure 8a, when five 355 nm nanosecond laser pulses irradiated on fused silica, no apparent damage was observed because the photon energy (B3.5 eV) of the 355 nm wavelength is lower than the bandgap (B8.9–9.3 eV) of fused silica, resulting in little absorption of laser energy .For five femtosecond laser pulses, a shallow hole was machined (Figure 8b), because of the low-efficient 800 nm (B1.55 eV) photon absorption through multiphoton ionization in Figure 13 (a) Dependence between the recast ratio (recast area/ablation wide bandgap dielectrics. In Figure 8c, femtosecond-nanosecond dual- area) and pulse delay on fused silica fabrication using femtosecond laser beam laser manufacturing revealed that a much higher fabrication pulse trains consisting of two identical sub-pulses with an identical total fluence. AFM profiles of the structures of the fused silica fabrication using efficiency (that is, a 50.7-fold enhancement in material removal (b) a conventional single pulse and (c) femtosecond laser pulse train with a volume) was obtained. This enhancement was attributed to the high pulse delay of 300 fs. SEM images of the structures on fused silica free electron density generated by the femtosecond laser pulses, which fabrication using a femtosecond laser pulse train with different energy ratios leading to the significantly-increased absorption of the nanosecond between the two sub-pulses: (d)1:1 and(e) 2:1. The femtosecond laser laser pulses energy. However, the femtosecond-nanosecond dual-beam with wavelength of 800 nm, duration of 35 fs and repetition rate up to − 2 system was complex, and the quality of the as-fabricated holes could 1 KHz. The total fluence of the pulse trains in all experiments is 5 J cm ; not be guaranteed. To date, temporally shaped femtosecond laser the scale bar in d and e is 500 nm. Reproduced from Ref. 128 (with the permission of OSA). double-pulse train was used to manufacture high-quality microholes with high-efficiency . Through temporal pulse shaping, free electron density can be adjusted to be around the critical point, at which the laser frequency equal to the plasma frequency, nearly optimizing to the resonance absorption so that the fabrication efficiency is enhanced 56 folds and the aspect-ratio is enhanced 3 fold (Figure 8d). Chemical etching controlled by temporal pulse shaping based on EDC. By designing femtosecond pulse trains on the basis of EDC, Flat-top beam Vortex beam the free electron density and corresponding photon-absorption efficiency can be modified. Thus, the material properties can be modified to improve the etching rate of fused silica (Figure 9). Compared with conventional femtosecond laser pulses, femtosecond Gaussian beam laser double-pulse trains achieved a 37-fold enhancement in the laser- assisted chemical etching rate. Simulations indicated that by optimiz- Bessel beam Airy beam ing the pulse delay between the two sub-pulses, the free-electron density can be modified, leading to the change of localized transient Figure 14 Shaping conventional Gaussian beam into different beam types. material properties, such as the physical properties (that is, reflectiv- ity), so that the laser field was reshaped, then the free electron density EDC, we can control the localized transient electron density to induce distributions (Figure 9k–9n) and absorbed laser intensity distributions resonance absorption between laser and its generated plasma; thus, (Figure 9o–9r) can be controlled, contributing to the enhancement of 119,120 high efficiency fabrication can be achieved .As illustrated in the photon absorption efficiency and result in the modification Light: Science & Applications doi:10.1038/lsa.2017.134 Recast ratio Electrons dynamics control LJiang et al ae Ablation area Metal nanowire 112 nm 190 nm 56 nm 1 μm 40 nm 20 nm 10 nm 0.0 0.5 1.0 1.5 2.0 2.5 Energy (keV) 5 μm cd f –200 –100 0 100 200 X (nm) (nm) Figure 16 Fabrication of nanowire by spatial pulse shaping. (a) Single spot fabricated by the shaped beam. (b) Scanning electron microscope (SEM) images of nanowire. (c and d) AFM images of the nanowire and its cross section. (e) EDXS measurements of the metal nanowire and the ablation area. (f) Five-ring patterns fabricated by the proposed methods . Reproduced from Ref. 129 (with the permission of Wiley). a b Electrode pads –2 220 J m –2 260 J m –2 –4 305 J m Resistor –2 345 J m –2 –8 390 J m –1000 –500 0 500 1000 Current (μA) Figure 17 Testing of the resistivity of the nanowire. (a) The SEM images of the nanowire and the electrode pads. (b) The volt-ampere characteristics of the nanowires. Reproduced from Ref. 129 (with the permission of Wiley). improvement in the irradiated zone. We also conducted micro-Raman which in turn improve the manufacturing efficiency. Overall, femto- spectroscopy to characterize the internal structure of the sample. As second laser temporal pulse shaping fabrication that is based on EDC Figure 9s–9u shows, by optimizing the pulse delays between the two represents a preliminary attempt to control the chemical reaction. sub-pulses, we can adjust the chemical properties (such as the Si-O bond structure) in irradiated material, which results in a higher Modulation of femtosecond laser-induced periodic surface structures based on EDC. On the basis of the aforementioned theory on EDC, number of 3- and 4-membered ring structures in double pulses we demonstrated that femtosecond laser-induced periodic surface modified regions than that in single pulse modified regions. These structures (LIPSS, also referred as ripples) can be deliberately results lead to the increases in the reactivity of the oxygen atoms, modulated by controlling the electron density and its distribution which contributes to a higher etching rate induced by femtosecond via designed femtosecond laser pulse trains. LIPSS have been studied laser double-pulse train. In short, by varying the pulse delay of pulse 142–144 extensively in various materials, including semiconductors , trains, we can control the localized transient free electron dynamics, 145 110,146,147 including bound electron ionization, free-electron density, tempera- metals and dielectrics , because of their promising 124,148–154 ture and excited state distribution, and further modify localized applications . The periodicity, orientation and structure are transient material properties, such as physical and chemical properties, the typical parameters in the study of ripples. According to its doi:10.1038/lsa.2017.134 Light: Science & Applications Height (mm) Voltage/mV Intensity/(Arb. Unt.) Fs laser system Electrons dynamics control LJiang et al Coulomb explosion Electrostatic ablation Plasma expansion Impact ionization Shockwave Radiation Electron heating Structure evolution Photoionization Vaporization/Melting Time scale (s) –15 –14 –13 –12 –11 –10 –9 –8 –7 –6 –5 –4 –3 –2 –1 0 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 Time-resolved plasma ICCD/CCD Femtosecond laser pump-probe photography and LIBS imagery Multipulse induced Plasma expansion Laser propagation Early plasma and plasma evolution shockwave expansion Ejected material composition Laser-induced free Multipulse induced elctrons evolution Energy deposition Plasma density and temperature structure evolution a b Beam splitter Half wave plate Fs laser CCD Mirror Shutter Polarizer Illuminator Shutter Half wave plate Mirror Delay Pump stage Polarizer Lens Lens Lens Objective filter BBO Filter Probe Z 20× X Y PMMA CCD PMMA c d Lens Lens ICCD CCD Illuminator Zoom lens Zoom lens Objective Objective PMMA PMMA Figure 18 Schematic of multiscale measurement of femtosecond laser drilling, including laser propagation and laser-induced material excitation, plasma and shockwave evolution and hole formation and so on. (a) Pump-probe shadowgraph imaging technique. (b) Laser-induced breakdown spectroscopy (LIBS). (c) Time-resolved plasma photography with gated intensified charge-coupled device (ICCD). (d) Industrial continuous imagery. 159,170–172 periodicity, LIPSS can be divided into low spatial frequency LIPSS state into a metallic state . Subsequently, at the interface (LSFL) and high spatial frequency LIPSS (HSFL). It is now widely between the metallic state surface and air, SPs can be excited by the accepted that the excitation and propagation of surface plasmon coupling between the surface electrons of the irradiated sample and 155–157 polaritons (SPPs) plays a crucial role in LSFL formation .The the incident field when the real part of the dielectric function is less formation of LSFL is affected by the initial laser-SPPs interference and than − 1 (Ref. 173). The SPs are characterized by surface electro- 158,159 the subsequent grating-assisted SPPs-laser coupling effect .Up to magnetic waves so that the coupling field is a superposition of the now, the formation mechanism of HSFL is still under investigation. incident field and the SP field. When the free electron density reaches Recently, Wang et al demonstrated that structure evolution of LSFL 21 − 3 the critical density (B1.74 × 10 cm for the wavelength of and HSFL is highly dependent on the localized effective laser fluence, 800 nm), the SPs can be resonantly excited . The SPPs excitation which determines the instantaneous optical permittivity by the laser- and resonance can reshape the laser intensity distribution in the excited electrons creating an active plasma layer. In general, the material and affect the subsequent linear/nonlinear ionization process. 147,161 formation mechanisms include self-organization , second har- Thus, the transient free electron density and its distribution is the key 162–164 165 monic generation (SHG) third harmonic generation (THG) , factor that affects SP excitation and properties, and ultimately, the 166 167 168 excitation of SPPs ,split , Coulomb explosion and cavitation 110,147,173 corresponding ripple formation . Meanwhile, the increased instability and so on. When a femtosecond laser irradiates the localized electron density can be further affected by trapping, diffusion surface of dielectric/semiconductor materials, free electrons can be and recombination with a time scale of several hundred femtose- generated, leading to the formation of electron-hole plasma (surface conds. Therefore, the time delay within the picosecond timescale is plasma) with time scales shorter than the electron-phonon relaxation proposed to control the electron dynamics to modulate the electron time. The localized transient free electron density is rapidly increased density and distribution, thus to modulate the resulted ripple through linear and nonlinear (multiphoton and avalanche) ionization, structures. leading to the material transforming from a dielectric/semiconducting Light: Science & Applications doi:10.1038/lsa.2017.134 Electrons dynamics control LJiang et al Femtosecond-picosecond Picosecond-nanosecond Air probe 600 ps 4 ns 16 ns ab c d e fg h 1.0 SF 0.9 100 fs 300 fs 500 fs 700 fs 900 fs 0.8 Air 0.7 PMMA 60 μm 0.6 60 μm PMMA 0.5 Laser propagation and maetrial excitation Plasma and shockwave expansion Nanosecond-submillisecond Millisecond-second Air Air Air 200 μm 500 μm PMMA PMMA PMMA Air Air Air 0 ms 0.5 s 5 s 20 s 33 s 45 s 60 s PMMA PMMA PMMA Plasma images above the surface Plasma images in the hole Multiple pulse induced structure in second scale 200 μm 0 s 10 s 20 s 30 s 40 s 50 s 60 s 70 s 80 s 90 s 100 s 110 s 120 s Figure 19 Multiscale measurement results of deep-hole drilling process in PMMA. Multiscale measurement results of deep-hole drilling process in PMMA with 100 μJ pulse energy focused by plano-convex lens (f= 100 mm). The dynamics include femtosecond-picosecond electron excitation, picosecond– nanosecond plasma and shockwave evolution and multiple pulse-induced structure in second scale. Optimal EDC using suitably shaped temporal pulse trains thus orientation parallel to the laser polarization, were obtained by gives the possibility to modulate the LIPSS artificially, offering changing the pulse delay (Δt) and pulse fluence (Figure 10 upper extended flexibility in material processing. Studies show that for panel). Thus, three types of LIPSS under specific conditions can be femtosecond (fs) laser pulse train processing of materials, the pulse obtained. (1) LSFL with orientation parallel to the laser polarization delay between sub-pulses strongly impacts the formation of direction; (2) HSFL with orientation parallel to the laser polarization 7,176–179 176–178 nanostructures , especially the morphology of LIPSS . direction with low pulse fluence; and (3) HSFL with orientation Here, we performed relevant experiments as examples on the surface perpendicular to the laser polarization direction at higher pulse 122,123 of dielectrics to illustrate the aforementioned mechanisms . fluence. The experimental results indicate that: (1) at lower pulse Fused silica was used as a dielectric material in a case study on the fluences, a transition from LSFL to HSFL occurred at a pulse delay of control of the LIPSS period, area and orientation. For conventional 50 fs with a decrease in area (Figure 10a and 10b,); (2) whereas, at femtosecond laser irradiation, only LSFL with an orientation parallel higher fluences, LSFL were replaced by another type of HSFL with an to the laser polarization and HSFL with an orientation perpendicular orientation perpendicular to the laser polarization at Δt4100 fs to the laser polarization were obtained on fused silica depending on (Figure 10c–10h). The average periods of LSFL and HSFL were the laser fluences (F) or pulse number (N). Nevertheless, compared 560± 8 nm and 255± 30 nm, respectively. with the conventional situation, both types of ripples with controllable During processing, the second sub-pulse of the doublepulse train periods, areas and orientations, especially the HSFL with an significantly affects the free electron density and distribution generated doi:10.1038/lsa.2017.134 Light: Science & Applications Commercialized CCD imaging Plasma photography Pump and probe Electrons dynamics control LJiang et al 3 3 2 2 C swan (d Π -a Π ) CN (B Σ–X Σ) 2 g u SP DP @ 80 ps Δv =0 Δv =0 Na I Ca I Δv =+1 Ca II 2 2 CH (A Δ–X Π,Δv =0) Δv =–1 Δv =–1 380 400 420 440 460 480 500 520 540 560 580 Wavelength (nm) bc 75 8500 70 8000 65 7500 60 7000 55 6500 50 6000 45 5500 0.1 1 10 20 40 60 80 100 120 0.1 1 10 20 40 60 80 100 120 DP delay (ps) DP delay (ps) − 2 Figure 20 Typical spectra of PMMA plasma irradiated by a single pulse and a double pulse with the same total fluence of 7.8 J cm . Reproduced from Ref. 134 (with the permission of SPIE). by the first sub-pulse, thereby influencing the mechanism of LIPSS recombination, which is determined by the electron density. By formation and the surface morphology. Consequently, by controlling controlling the pulse delay, the electron occupation is adjusted , the pulse delay and pulse fluence, we can control excited electron leading to the manipulation of electron density, thus facilitating production, distribution, motivation and the interaction between SHG and resulting in a 50% cut in LIPSS periods; surface plasmon (SP) and the incident laser, then to control the 3) The strong decrease in the ripple area was due to the electron periodicity, orientation and morphology of the LIPSS. decay. The free electrons in the conduction band excited by the first sub-pulse relaxed and returned to the valence band during the 1) LSFL obtained here at a low pulse fluence were oriented parallel to time interval between the two pulses, in terms of diffusion and 184–186 the laser polarization, which cannot be explained by SPs recombination , leading to the reduced electron density; 147,180 excitation . According to calculations based on the Sipe- therefore, the energy coupling by the second sub-pulse to the 180 53 Drude model and pump-probe results , femtosecond energy excited material decreased, leading to the decay of ablation, which deposition by the first/previous pulses can occur at specific (LSFL) resulted in the reduced rippled area; spatial frequencies, reshaping the electron density distribution 4) The transition of HSFL at a higher pulse fluence was mainly along the polarization direction, and then determines the forma- attributed to the periodic plasma enhancement of the incident laser tion of LSFL with an orientation parallel to the polarization field, which is related to the excited free electron density. When the direction. Meanwhile, thermal effects also play a critical role in pulse delay was low (o100 fs), the effects of induced SPs on the subsequent material removals at the fluences for LSFL incident laser were insufficient due to the low absorbed intensity. 181–183 generation ; With pulse delays of 100–500 fs, however, SPs excitation can easily 2) The periods of LSFL and HSFL (with orientation parallel to the be achieved at the initial stage of the second sub-pulse due to the polarization direction) were close to the fundamental (λ/ accumulation of the first sub-pulse (see modeling section n= 551 nm) and second-harmonic (λ/2n= 275 nm) wavelengths for details). The interaction between SPPs and the incident laser in fused silica with lower pulse fluence. Thus, SHG plays a key role field resulted in the periodic modulated intensity enhancement at in HSFL formation with an orientation parallel to the polarization the surface. The incubation effects with multiple bursts irradiation 155,156 direction at a lower pulse fluence. SHG is a result of electron led to the evolution of the local intensity distribution along the Light: Science & Applications doi:10.1038/lsa.2017.134 Intensity (a.u.) 16 –3 Electron density (×10 cm ) Temperature (K) Electrons dynamics control LJiang et al 0 ps 0.2 ps 0.3 ps 0.4 ps 0.5 ps 0.6 ps 0.7 ps 0.8 ps 2 ps 6 ps 10 ps 40 ps 80 ps 120 ps 1mm –2 3.3 J cm –2 4.4 J cm –2 6.6 J cm –2 8.8 J cm –2 11 J cm bc 2.6 –2 –2 6.6 J cm 6.6 J cm –2 –2 8.8 J cm 8.8 J cm –2 2.4 –2 11 J cm 11 J cm 2.2 2.0 1.8 1.6 1.4 0 0.2 1 10 40 80 120 0 0.2 1 10 40 80 120 Interpulse delay (ps) Interpulse delay (ps) Figure 21 Characterizations of femtosecond laser double pulse induced plasma of fused silica (a) plasma images, (b) electron densities and (c)plasma temperature. electric field direction , resulting in HSFL formation at a specific analysis of chemical and biological species because of its high 188,189 period with an orientation perpendicular to the electric field sensitivity and fingerprint-identification features . Previous stu- direction. dies have demonstrated that surface morphologies (for example, 190,191 ripples, nanoparticles and nanopillars) play key roles in SERS More complicated morphology control of LIPSS on fused silica can enhancement in which the enhanced electromagnetic field on surface be obtained by EDC via shaping the conventional femtosecond laser nanostructures induced by the localized surface plasmons resonance 191–193 pulses into symmetrical triple- and quadruple-pulse sequences, (LSPR) effect dominated . Consequently, tailoring the surface structures into different morphologies and sizes, in this aspect, is especially, double-grating structures and an HSFL period as small as significant for tuning LSPR features to further improve the sensitivity 190nmwereobtained(seeFigure10i–10o) . In addition, the and push SERS devices into practical applications. As shown in geometric morphology modulation on Si can also be obtained by Figure 11d, the SERS intensity gradually reached a maximum when EDC via adjusting double pulse delay . Therefore, the aforemen- the pulse delay was increased from 0 to 800 fs, and then decreased tioned experimental studies demonstrate that by designing a femto- when the pulse delay was further increased to 1000 fs . Although the second laser pulse train, the electron dynamics can be controlled, that small changes (that is, the modulation of pulse delays) occurred in the is, electron density, distribution, thus to control the coupling between incident pulses, the obtained surface structures significantly improved the SPP and incident laser, and finally guide the material response the signal intensity (Figure 11e). Compared with the conventional (LIPSS morphology) towards user-designed directions with various femtosecond laser ablation (Δt= 0 fs), the designed pulse train could morphologies (for example, periods, orientations, distributions and reduce more silver nanoparticles (Figure 11a and 11b) and lead to a geometric morphologies). more uniform distribution of the nanoparticles deposited on the Detection sensitivity improvement of SERS based on EDC.By subwavelength ripples (Figure 11c). Consequently, the SERS sensitivity designing femtosecond laser pulse trains on the basis of EDC, we was improved. can control the localized transient electron density and its distribution To further confirm that we controlled the properties of SPs, a two- to modulate the properties of SPs, then promote the energy transfer to step experiment was conducted using a designed pulse train, as materials and control the surface structures and photochemical schematically shown in Figure 11f . The laser polarization and pulse reduction process; thus, high detection sensitivity SERS substrates delay were synergistically controlled. Compared with the ripples formed can be achieved. SERS has been recognized as the most promising by conventional femtosecond laser pulses (Figure 11g and 11h), regular trace analyte detection method for rapid and accurate label-free nanopillar arrays were generated by a double pulse train (i.e., Δt40fs) doi:10.1038/lsa.2017.134 Light: Science & Applications 16 3 Electron density (×10 cm ) Electron temperature ( °C) Electrons dynamics control LJiang et al M M Bessel beams ac 2 3 CCD2 0 1 0 0.5 –5 –5 Initial Bessel beams –10 0 500 1000 1500 –10 –5 05 10 DM y (μm) z (μm) PBS Axicon Gaussian beams 10 1 Micro-Bessel beams MO S2 HWP z 0 0.5 –5 –5 CCD1/ICCD L x 0 –10 0 0 500 1000 1500 –10 –5 05 10 Fs laser y M y (μm) z (μm) S1 b d Bessel beams Gaussian beams 1 μm 5 μm 5 μm 5 μm 3.5 2.5 5 μm 50 μm μm 10 μm 123456 μm Figure 22 (a) Schematic of the spatial shaping femtosecond laser pulses microdrilling setups; (b) Morphology images of microholes drilled with a single pulse Bessel beam and Gaussian beam, respectively; (c) Intensity distribution simulations of the Bessel beam and Gaussian beam; (d) The hollowness characterization of microholes drilled by single-pulse Bessel beam. Reproduced from Ref. 132 (with the permission of Springer). with the polarization direction rotated by 90°, as shown in Figure 11i on the basis of EDC, bound electrons can be ionized, chemical bonds and 11j. In the case of conventional femtosecond laser ablation between atoms can be interrupted, atoms can be selectively removed and material properties can be modified so that gold cations can be (Δt= 0), SP induced by the previous pulse of the subwavelength ripples would fast damp to its original state before the subsequent pulse spontaneously reduced on MoS surface , as shown in Figure 12. Compared with the conventional femtosecond laser pulses, the arrived because of the long-time interval (millisecond scale), leading to chemical reduction rate of gold cations on MoS surface modified the formation of large and non-uniformly distributed silver nanopar- by femtosecond double-pulse trains was significantly enhanced (Figure ticles on conventional grating-like ripples. By contrast, because the 12a–12h). We conducted X-ray photoelectron spectroscopy (XPS) and pulse delay of a pulse train is shorter than the damping time of SP, the atomic force microscope (AFM) to characterize and analyze the properties of SPs would be significantly controlled to change the energy element valence, atomic ratio and morphology of the sample, as transfer efficiency, which could contribute to the formation of showninFigure 12i–12m. When material was irradiated by the abundant small and uniformly distributed silver nanoparticles on the femtosecond laser pulses, electrons were excited from bonding to anti- nanopillar arrays. In addition, the larger enhancement of the incident bonding states (bound electrons were ionized), hence abundant laser electric field on nanopillar arrays lead to generation of much more chemical bonds between atoms (such as Mo-S bonds) were instanta- silver nanoparticles, which resulted in a lager SERS signals (maximum neously weaken or even broken entirely, resulting in an integral enhancement factor up to 2.2 × 10 ). Overall, by designing a pulse train increase in the binding energy of S atoms and the appearance of on the basis of EDC, the electron density and its distribution induced unbound sulfur (unsaturated S atoms with dangling bonds). Mean- by laser irradiation can be controlled to modulate the properties of SPs, while, abundant Mo atoms were selectively removed, and the surface resulting in a more effective energy transfer and changes in the resulting lattice structure of the material was non-perfect and broken into a structures. These findings provide new insights regarding tuning LSPR large number of micro/nano debris terminating with unsaturated S features for related applications. atoms (edge active sites). These unsaturated-terminal S atoms induced on MoS can reduce gold cations to gold atoms (Figure 12n). Surface chemical reduction activity of MoS controlled by temporal 2 2 pulse shaping based on EDC. By designing femtosecond pulse trains Compared with conventional femtosecond laser pulses, more Mo Light: Science & Applications doi:10.1038/lsa.2017.134 r (μm) r (μm) x (μm) x (μm) Electrons dynamics control LJiang et al a c 0~5 ns 5~10 ns 10~15 ns 15~20 ns 20~25 ns 25~30 ns 30~35 ns 200 fs c1 400 fs 0.5 ns c2 c11 Air 600 fs c3 PMMA 800 fs 1 ns c12 c4 1 ps c5 100 μm 1st pulse 1.2 ps 3.6 ns c6 c13 1.4 ps c7 0~5 ns 5~10 ns 10~15 ns 15~20 ns 20~25 ns 25~30 ns 30~35 ns 1.6 ps 8.6 ns c8 c14 10 ps Air c9 13.6 ns PMMA 100 ps c15 c10 The final hole fs laser 18.6 ns c16 100 μm c17 200th pulse 485 μm Air PMMA Figure 23 Time-resolved images of the plasma expansion during the Gauss beam drilling for the 1st pulse (a) and 200th pulse (b), respectively ,pump- probe shadowgraph study (c) of Bessel beam drilling . c1–c16 Time-resolved images of femtosecond-picosecond-nanosceond dynamics of Bessel beam drilling. c17 The final hole morphology drilled by Bessel beam. Reproduced from Ref. 133 (with the permission of OSA) and Ref. 131 (with the permission of OSA). atoms were selectively removed and more unbound sulfur was including melting (various phase change mechanisms coexisted during formed, which might result from the significantly enhanced ionization the process) dominated the ablation process, which resulted in inferior of bound electrons by femtosecond double-pulse trains (first, bound fabrication quality. While in the case of femtosecond laser pulse train electrons were ionized by the first sub-pulse and the electron-hole with the identical total fluence, a single pulse was split into two sub- pairs separated; second, the photoexcited electron-hole pairs recom- pulses with the identical fluence lower than the ablation threshold of bined, which began within 500 fs and would last for more than a fused silica. Therefore, through optimizing the pulse delay between the hundred picoseconds; last, the second sub-pulse could further ionize two sub-pulses, the laser-induced electron density can be controlled to the recombining electron-hole pairs, which might further facilitate the be slightly higher than the critical density. Subsequently, the non- 38,194,195 photochemical bond breaking) and inevitably resulted in thermal phase-change mechanisms mainly dominate the fabrication stronger chemical reduction activity and higher chemical reduction process, contributing to less recast and high fabrication quality. rate on laser-treated MoS . In conclusion, by designing pulse trains Furthermore, by adjusting of the energy ratio between the two sub- based on EDC, the bound electrons ionized by laser irradiation can be pulses, the ionization processes were altered to change the pulse energy controlled to modulate chemical bond cleavage and increase chemical absorption. Consequently, the free electron distribution can be reduction ability, which can reduce gold cations to obtain metal-MoS adjusted to make the phase change process a nonthermal one. hybrids for relevant applications. Therefore, we achieved much higher fabrication quality and more controllable structures (Figure 13d and 13e). Overall, through High-quality fabrication by temporal pulse shaping based on EDC. femtosecond laser temporal pulse shaping, the localized electron By designing femtosecond laser pulse trains on the basis of EDC, the dynamics of the materials can be controlled, thus the phase-change phase change mechanism can be controlled to achieve high-quality mechanisms can be adjusted to be dominated through nonthermal fabrication . As Figure 13a shows, the recast ratio (recast area/ phase change process so that higher fabrication quality and more ablation area) of the ablation structures on the fused silica decreased as controllable structures can be achieved. the pulse delay increased within a femtosecond laser double-pulse train. Compared with a conventional single pulse, the recast height Femtosecond laser spatial pulse shaping surrounding the ablation spot decreased by 60% when fabrication was Spatial laser shaping is another key aspect of EDC. The spatial states of performed using a femtosecond laser pulse train (see Figure 13b and electron dynamics are closely related to the spatial distribution of laser 13c, for atomic force microscope (AFM) profiles). In this case, the − 2 energy. Conventionally, most commercial laser systems provide a total fluence of the pulse train was 5 J cm (greater than the ablation fundamental Gaussian intensity profile, which has a short Rayleigh threshold of fused silica). The conventional femtosecond laser single length, small beam waist of focusing spot, homogeneous phase/ pulse induced an electron density much higher than the critical density polarization state and limited numbers of laser spots. These char- and a higher Coulomb barrier, therefore leading to significant electron acteristics of Gaussian laser cannot satisfy the increasingly higher screening effects. Consequently, the accumulation of positive charge demands for flexibility, precision and efficiency in high-end laser during ablation was reduced, thus the electric field was weakened so fabrication. Thus, it is essential to spatially shape the laser pulses to that the effectiveness of the nonthermal processes (Coulomb repulsion and/or electrostatic ablation) was reduced. The thermal phase change control the spatial electron density/temperature distribution. By spatial doi:10.1038/lsa.2017.134 Light: Science & Applications Electrons dynamics control LJiang et al a c 50 μm bd 2.23 μm 5 μm 20 μm Figure 24 (a and b) The image of a microhole array throughout a 1 cm × 1 cm large area using the flying punch method for Bessel beam; (c) the local magnified OTM image of a microhole array; (d) the local magnified OTM image of microholes in the central array under 100 × microscope objective. Reproduced from Ref. 132 (with the permission of Springer). pulse shaping, the Gaussian beam can be converted into different By using a spatially modulated femtosecond beam based on spatial beam types, such as flat top beam, vortex beam and Bessel beam, Airy EDC, we achieved high-resolution nanowire patterning that breaks the beam and so on (Figure 14) Some unique characteristics of these beam light diffraction limit . The initial beam had a Gaussian intensity types can lead to special functionality for laser fabrication, which has a profile with an even phase wave front. The phase of the initial beam large advantage over conventional Gaussian beams. was modulated using a liquid crystal on silicon spatial light modulator. After modulation, a relative phase difference was created between two Spatial pulse shaping experimental setup. The experimental setup is equal parts of the incident beam. Subsequently, the modulated beam shown in Figure 15. An amplified Ti: sapphire laser system (Spitfire was focused by an objective lens, forming a dual-peak focusing spot Ace-35F Spectra Physics Inc.) provides a fundamental Gaussian mode with an intensity valley in the center because of a diffraction effect. with a central wavelength of 800 nm and a pulse duration of 35 fs. The spatially modulated beam was used to pattern a gold thin film, Phase patterns were generated using the liquid crystal on a silicon which was deposited on a silica substrate by electron beam evaporator. spatial light modulator (SLM, Holoeye Pluto). The size of the liquid The central thin film was preserved due to the intensity valley thus a crystal screen was 15.36 × 8.64 mm. The modulated laser beam passes nanowire was formed in the beam center. Arbitrary nanowire can be through a 4f relay system, which consisted of two plano-convex lenses generated on the substrate by dynamically adjusting the orientation of (L1, L2). The distance between the SLM and lens L1 was equal to the the intensity valley. A minimum nanowire width of approximately L1 focal length. The distance between the two lenses was the sum of their focal lengths. The laser beam can then be transmitted to the focal 56 nm (B1/14 of the laser wavelength) can be achieved (Figure 16a– plane of the second lens L2 without any distortion. The focal length 16d). The high resolution is achieved by combining the ultrashort was 100 mm (L1) and 150 mm (L2), respectively. The samples were nature of the femtosecond and the low thermal conductivity of the mounted on a nanometer-precision stage (Newport, NPXYZ100) has a thin film. We tested the amount of Au residue using an electron resolution of 0.2 nm in the X, Y and Z planes. The fabrication process diffraction X-ray spectrum (EDXS) experiment (Thermo Scientific, was monitored using a CCD camera and white light source (WS). MA, USA; Figure 16e). By changing the direction of the phase pattern loaded on the SLM, the center nanowire changed its orientation High precision nanowire fabrication by spatial pulse shaping based accordingly. By dynamically adjusting the direction of the phase on EDC. Diffraction is a universal phenomenon in wave optics, pattern to tangent of the scanning route, arbitrary curves can be which greatly limits the resolution of laser fabrication to half of the fabricated. Figure 16f shows Olympic rings fabricated using this wavelength level. One of the greatest challenges in laser micro/nano method, which evidence its effective curve patterning ability. fabrication is ‘overcoming’ the diffraction limit. Several methods have To investigate the electronic characteristics of the nanowire, the volt- been developed to solve this problem using different mechanisms, such 196,197 198 ampere characteristics curve of the nanowire fabricated by different as near-field fabrication , two/multi-photon polymerizations energy pulses was measured. (Figure 17) No inner nanopores and and plasmonics-based fabrications .However,these methods all have 196,197 particle intervals are generated inside the nanowire, thus endowing the different disadvantages, including the low efficiency ,complex 196,197 199 procedures ,weak flexibility . and limited types of applicable nanowire with favorable electronic characteristics: the conductivity of 198 7 − 1 the nanowires wasashighas1.2×10 Sm ,and themaximum materials . A novel simple, repeatable, mask-free, high-throughput, 8 − 2 broad-applicability and high-flexibility method is highly desired. current density was up to 1.66 × 10 Am . This approach offers a Light: Science & Applications doi:10.1038/lsa.2017.134 Electrons dynamics control LJiang et al simple, robust alternative for high-quality nanowire fabrication as a strongly depended on the pulse delay; the signal was stronger than that complementary method to conventional lithography methods. of single-pulse irradiation when the pulse delay exceeded 10 ps. The maximum enhancement value was about 7 at the pulse delay of MULTISCALE MEASUREMENT OF ELECTRON DYNAMICS approximately 80 ps. To determine the electron dynamics difference, DURING FEMTOSECOND LASER-MATERIAL INTERACTIONS the plasma temperature and electron density were calculated using the To comprehensively understand the electron dynamics during femto- Saha–Boltzmann plot and the Stark broadening method, respectively. second laser micro/nano fabrications, a multiscale measurement The plasma temperature variation was correlated with the signal system was developed to monitor the spatiotemporal electron enhancements achieved using the double-pulse delay, implying that dynamics of laser-material interactions. This system integrates the the plasma reheating effect of the second pulse was the main 198,199 widely applied ultrafast pump-probe microscopy , time-resolved mechanism of the enhancement effect. By contrast, the maximum plasma photography with a gated intensified charge-coupled device enhancements were different for molecules, atoms and ions, which 200,201 202–204 205 (ICCD) ,LIBS and industrial continuous imagery were related to the different upper excitation energy of the emission 200,201 (Figure 18). These techniques have different characteristic time transition and the ionization stage of the electrons . resolutions, ranging from femtoseconds to seconds. By virtue of the A comparison of fused silica plasma emission induced by a single multiscale ability, the electron dynamics of femtosecond laser micro/ pulse and that induced by a double pulse is presented in Figure 21a nano fabrications can be revealed at different scales. and 21b. The plasma emissions at the pulse delay below 10 ps were As a case study, the deep-hole drilling process in poly (methyl stronger than they were in the single-pulse case. This phenomenon methacrylate) (PMMA) was investigated using the multiscale measure- was explained by the existence of free electrons and self-trapped ment system, as shown in Figure 19. In the femtosecond to picosecond excitons at this pulse delay. The free electrons and self-trapped time scale, focused laser propagation in the material, as well as free excitons left by the first pulse increased the absorption efficiency of electron generation, diffusion and recombination, was detected the second pulse. Moreover, extraordinarily high enhancement factors through pump-probe shadowgraphy (Supplementary Movie 2, were observed at a pulse delay above 10 ps. The maximum enhance- Supplementary Information). Nonlinear phenomena were observed ment factor of double-pulse irradiation wasB35 times at a pulse delay − 2 at a high pulse energy as a result of the strong self-focusing effect of of 120 ps for a fluence of 11 J cm . The plasma consisted of a fast the intense femtosecond laser. In the picosecond to nanosecond time component (ionized atoms and ions) and a slow plume component scale, the early plasma expansion and following shockwave evolution (partially ionized nanoparticles), with the slow part contributing little in the atmosphere can also be studied through pump-probe shadow- to plasma emission . The ionization of the slow part by the second graphy (Supplementary Movie 3, Supplementary Information). Shock- pulse greatly increased the plasma quantity and was demonstrated to wave expansion properties were investigated on the basis of the Sedov- be the main cause of the high enhancement factor when the pulse Taylor solution, revealing the change of environment and its effect on delay was larger than 10 ps. The plasma temperature was also laser ablation. For larger time scales, time-resolved plasma photo- calculated to ascertain its changes with respect to the pulse delay graphy with the gated ICCD and LIBS were employed to collect the (Figure 21c). Hence, the electron dynamics (free electrons, self- plasma expansion morphology (Supplementary Movie 4, trapped excitons or nanoparticles) left by the first pulse dominated Supplementary Information) and plasma emission spectroscopy, the plasma emission intensity of the femtosecond double pulse LIBS. respectively. These results not only provided the expansion dynamics of plasma but also the plasma intrinsic information, including the APPLICATIONS OF EDC IN HIGH QUALITY AND ASPECT-RATIO material composition, plasma density and temperature. In the milli- MICROHOLES DRILLING second to second time scale, the plasma evolution induced by multiple As discussed before, because of the significant electron-lattice none- pulses was studied using time-resolved plasma photography, revealing quilibrium state in the femtosecond laser fabrication process, the laser- the effect of the prior structure on the plasma intensity and material interactions, including phase change and material removal, 26,27 distribution (Supplementary Movie 5, Supplementary Information). are actually determined by the initial photon–electron interactions . Furthermore, microhole formation was studied using industrial Furthermore, the electron density distribution can be manipulated by continuous imagery (Supplementary Movie 6, Supplementary modifying the laser intensity distribution. By implementing EDC Information), providing information concerning the evolution of the through spatial pulse shaping, forming an intense, long and uniform depth, diameter and quality of the high-aspect-ratio microhole. Bessel beam, to adjust the localized transient electron density By using the measurement system, we determined the nanosecond- distribution, and thus control phase change, we fabricated high quality scale electron temperature and density evolution under double pulse and high aspect-ratio microholes. The technique was applied in key irradiation. Large differences in plasma plume ejection of PMMA and structure fabrication in one of the 16 Chinese National S&T Major fused silica were demonstrated using the LIBS technique, following the Projects. different electron dynamics designed by femtosecond laser double- Microholes fabrication in the key structure faced many challenges, pulsetrain. Bychanging the pulsedelay, wecould modify the including the high aspect-ratio (420:1), small diameter (o10 μm), generation and distribution of electrons and thus effectively improve taper-free, high quality, reduced recast/ejected materials and minimized the electron density and temperature. Finally, we could improve the in-cavity residues. Although various novel methods had been proposed plasma enhancement factor and optimize the ablation accuracy for microholes fabrication, there were many challenges that limited through EDC. their application in the key structure fabrication. By using a particle The plasma spectra of PMMA induced by single pulse and double based near-field nanostructuring to overcome the optical diffraction pulse irradiation at the same laser fluence were compared to explore limit, Quentin et al achieved nanoholes with diameters below the LIBS enhancement mechanism. The plasma spectra of PMMA 200 nm, but the structure was a shallow crater and the aspect-ratio consisted of emission peaks of molecular species (CN, CH and C )as was below 1:1. He et al reported that high aspect-ratio microholes well as atomic and ionic species (Ca I, Na I and Ga II) .AsFigure 20 were fabricated in fused silica with femtosecond laser transverse direct- shows, double-pulse-induced plasma emission signal enhancement writing followed by wet chemical etching, but the cross-sections of the doi:10.1038/lsa.2017.134 Light: Science & Applications Electrons dynamics control LJiang et al microholes fabricated by this method are usually in poor shape Bessel beam and a focused Gaussian beam in the longitudinal and ascribing to the asymmetric shape of the focal spot. Gottmann transverse plane, as shown in Figure 22c. The micro-Bessel beam et al constructed a selective laser-induced etching system of 3D exhibited a large focal depth (Z = 597 μm) in comparison with the max precision quartz glass components. By the combination of a three-axis focused Gaussian beam (Raleigh range R= 18.5 μm). According to the system to move the glass sample and a fast 3D system to move the laser simulation results, the single-pulse Bessel beam exhibited an intense, focus, the selective laser-induced etching process (LightFab 3D Printer) long, consistent and stable intensity distribution without any nonlinear is suitable to produce more complex structures in a shorter time. beam distortion along the propagation direction. These unique However, the over-etching for some complex structures and the properties of the spatial shaping pulse allow for uniform energy formation of cracks during the etching process limit advanced deposition over extended propagation lengths and then adjust free applications. In addition to the aforementioned methods, percussion electron density distribution to be intense, long and uniform through 204,206,207 drilling has attracted the most interests . Nevertheless, the photo–electron interactions. aspect-ratio of microholes with percussion drilling is usually smaller Using the aforementioned multiscale time-resolved measurement than 10:1 in air because of the saturation effect and the occurrence of system developed by us, we further revealed the forming mechanisms the bending effect. We recently reduced the ambient pressure from of the high aspect-ratio and high quality microholes using a approximately 10 Pa (air) to B1 Pa (rough vacuum) to control the femtosecond laser spatial pulse shaping beam (Bessel beam) by 133,208 expansion dynamic of the ablated plasma/material .The aspect comparing it with the Gaussian beam percussion drilling. Figure 23a ratio of the microholes was significantly improved from ~ 40:1 (in air) shows the time-resolved plasma photography of a Gaussian beam to approximately 100:1 (in vacuum), and the bending effect was drilling with a single-pulse, the ablation plume was detected only simultaneously eliminated. Nonetheless, conventional Gaussian beam above the PMMA surface. As shown in Figure 22c, the focused drilling could barely achieve a much higher aspect-ratio of taper-free Gaussian beam was condensed, thus, the laser energy could only be microholes with diameters below 10 μm. localized within several microns below the surface. Such a high-energy To overcome the disadvantages of the traditional Gaussian beam concentration would make free electron density much higher than the drilling, we proposed a new processing method for high aspect-ratio critical density. Therefore, the Coulomb explosion and electrostatic and high-quality microdrilling by optimizing the localized transient ablation are weakened, and the material removal is mainly attributed electron density distribution in plasma in the focal spot through spatial to melting and evaporation, resulting in a large heat-affected zone. pulse shaping . Taper-free microholes with a diameter of approxi- Meanwhile, the shielding-effect induced by the ultrahigh free electron mately 1.6 μm and an aspect ratio of up to 330:1 using a single spatial density suppresses the laser propagation, leading to a short optical shaping pulse (Bessel beam). The aspect ratio of these fabricated penetration depth; consequently, only a shallow crater can be microholes was 52 times larger than those formed by using a Gaussian fabricated with a single-pulse Gaussian beam. Furthermore, during beam in similar focusing conditions. The formation of these high the multipulses percussion drilling process (Figure 23b), the subse- aspect-ratio microholes is attributed to the intense, long and uniform quent pulses energy was not sufficiently absorbed by the material transient localized electron density distribution, which is adjusted on because of the strong reflection of the generated dense ablated plasma. the basis of the unique intensity distribution and propagation stability The generated plasma would disturb the propagation of the subse- of the Bessel beam through spatial pulse shaping. quent laser pulses so that leads to an unstable filament bending slightly In our experiments, we normally used an axicon (Edmund Inc., deep inside the hole, then the laser beam energy was deposited along a Barrington, NJ, USA, base angle α= 2°, refractive index= 1.45) to bent direction. Moreover, for the deeper microholes, the plasma transform a Gaussian beam into a Bessel beam for high aspect-ratio, would cool down gradually, then adhered to the microholes sidewalls, high quality and high efficiency microdrilling. This is a convenient and which negatively affected the quality of the microholes. effective method for spatial pulse shaping. Figure 22a shows the By contrast, spatial shaping pulse (Bessel beam) drilling exhibits a schematic of the experimental setup used for femtosecond laser completely different fabrication mechanism. Figure 23c displays the microhole drilling . Subsequently, we performed the experiments time-resolved images of Bessel beam drilling in PMMA. As the Bessel to prove that high aspect-ratio and high quality microholes can be beam entered into the PMMA, the rising edge induced electron obtained by using a single spatial shaping pulse. As Figure 22b shows, excitation, which formed a dark strip along the light path that a microhole with a mean diameter of 1.5 μmand adepth of indicated the plasma channel dynamics. The maximum length (optical approximately 523 μm can be drilled using a single-pulse Bessel beam penetration depth) of the Bessel plasma channel was reached at 1.6 ps at a pulse energy (E) of 20 μJ in PMMA. Under identical processing and effectively coincided well with the final hole-depth, indicating the conditions, the diameter and depth of a microhole fabricated using the dominant role of initial electron excitation and the free electron Gaussian beam were 7.2 and 41 μm, respectively. Subsequently, to distribution on the final structure formation. In the following tens of further demonstrate the hollowness of the microholes through single picoseconds, there existed electron-ion energy transfer occurred that spatial shaping pulse drilling, the two different methods were mainly induced an extremely high pressure and temperature in the focal area, adopted, namely the liquid infiltration method and the cross-section resulting in an explosive and supersonic expansion of the material in profile test, as shown in Figure 22d. The microholes drilled using this the nanosecond domain, as shown in Figure 23c11–23c16. In contrast method were small in diameter and exhibited a high aspect-ratio, a to the pressure waves in air and some other cases, the pressure wave taper-free sidewall, a highly circular entrance and fewer surface-ejected induced by the Bessel beam in PMMA was a cylindrical shockwave and materials. gradually expanded outward along the radial direction. This phenom- Moreover, we performed a series of theoretical investigations and enon suggests that the formation of microholes is an extrusion effect, simulations on the spatial intensity distribution of Bessel beams and leading to the formation of high quality microholes with taper-free, Gaussian beams to more thoroughly understand the physical mechan- reduced recast/ejected materials and minimized in-cavity residues. isms of pulse microdrilling with spatial shaping. On the basis of the After the understanding of the processing mechanisms of the Bessel theory of spatial pulse shaping using the axicon, we respectively beam, we proposed the flying punch method for machining large-area simulated the spatial intensity distribution of single-pulse, micro- microhole arrays in PMMA by using a Bessel beam, as shown in Light: Science & Applications doi:10.1038/lsa.2017.134 Electrons dynamics control LJiang et al Figure 24. A 1 cm × 1 cm microhole array (with 251 001 ultrahigh- system limited the spatial/temporal resolution and the sensitivity aspect-ratio holes in total, at a processing speed of 100 holes of the measurement, which shall be improved to monitor the per second) was fabricated within 42 min, indicating the high- electron dynamics in more detail. In addition, more time-resolved efficiency and high repeatability of the fabrication process using the measurement techniques should be integrated into the multiscale flying punch method . system, so that it can provide more information of electrons from different aspects. Furthermore, the measurement results should CONCLUSIONS AND OUTLOOKS correspond with the theoretical models, such as providing the In this paper, we comprehensively reviewed our decade-long efforts on characteristic values of electrons for theoretical calculations. four parts of EDC in femtosecond laser micro/nano fabrications: the (4) Broader applications: High-aspect-ratio and high-quality micro- theoretical fundamentals, experiments, multiscale measurements and holes by using spatial shaping pulses have been applied in applications. Theoretically, based on the four models with different fabricating some key structure fabrications. However, the current − 3 − 15 − 3 − 10 time scales (10 –10 s) and space scales (10 –10 m), we method is limited to a few transparent materials. We will extend demonstrated that the localized transient electron dynamics (including the range of materials and then explore novel applications of electron density, temperature and excited state distribution), and microholes, such as for microfluidic devices and three- subsequent phase change can be controlled by temporally/spatially dimensional integrated chip packaging. In addition, we will shaping femtosecond laser pulses. Experimentally, seven experiments substantially expand the applications of the novel method, such were reported as examples to validate the feasibility of EDC by as for the adjustment of chemical reaction pathways by ultrafast temporally/spatially shaping femtosecond pulses in micro/nanofabri- laser EDC, transient or permanent adjustment of material proper- cation. The experiments revealed that the precisions, efficiencies and ties through ultrafast laser EDC. qualities can be significantly improved and that various surface micro/ nano-structures can be effectively modulated by the proposed EDC- based methods. Additionally, multiscale measurements further directly CONFLICT OF INTEREST demonstrated the fundamentals of EDC from femtosecond scale to The authors declare no conflict of interest. nanosecond scale and to millisecond scale. Finally, EDC was applied in high aspect-ratio (330:1) and high-quality microholes drilling at the ACKNOWLEDGEMENTS speed of 100 holes per second (251 001 holes fabricated in 1 cm × 1 cm This research was supported by the National Natural Science Foundation of area within 42 min), in which multiscale measurements were used to China (NSFC) (Grant Nos. 90923039, 91323301, 50705009, 51105037, analyze and optimize the electron dynamics. The high aspect-ratio 51322511 and 51025521), National Basic Research Program of China (973 microholes drilling was applied to key structure fabrication in one of Program) (Grant No. 2011CB013000), the 863 Project of China under Grant the 16 Chinese National S&T Major Projects. No. 2008AA03Z301, the Cultivation Fund of the Key Scientific and Technical Innovation Project, Ministry of Education of China (No. 708018), the 111 We have devoted the past ten years to studying the mechanisms, Project of China (Grant No. B08043), Multidisciplinary University Research methodologies and applications of femtosecond laser micro/nano Initiative (MURI) program of USA under Grant No. N00014-05-1-0432 and fabrications. However, many challenges still remain, especially on National Science Foundation of USA under Grant No. 0423233. We thank Prof the following topics: Xin Li, Prof Jie Hu, Prof Sumei Wang, Prof Xiaowei Li, Prof Jingya Sun, Prof Jianfeng Yan, Dr Xiaoxing Su, Dr Weina Han, Dr Qingsong Wang, Dr Zhitao (1) Comprehensive models: Femtosecond laser-material interactions Cao, Mr Zhi Wang, Dr Changji Pan, Dr Peng Ran and Dr Pei Zuo for dis- are a comprehensive nonlinear, nonequilibrium process ranging cussing, revising and proofreading this manuscript. from a nanometer scale to a millimeter scale and from a femtosecond scale to a millisecond scale. However, our present models, including the plasma and improved two-temperature models, are not applicable to some materials. Furthermore, a 1 Malinauskas M, Žukauskas A, Hasegawa S, Hayasaki Y, Mizeikis V et al. Ultrafast laser comprehensive, integrated multiscale physical-chemical modeling, processing of materials: from science to industry. Light Sci Appl 2016; 5: e16133, doi:10.1038/lsa.2016.133. from a nanometer scale to a millimeter scale and from a 2 Sugioka K, Cheng Y. Ultrafast lasers—reliable tools for advanced materials processing. femtosecond scale to a microsecond scale, shall be developed Light Sci Appl 2014; 3: e149, doi:10.1038/lsa.2014.30. and improved to describe femtosecond laser-material interactions. 3 Gamaly EG. Femtosecond Laser-Matter Interaction: Theory, Experiments and Applica- tions. Boca Raton: CRC Press. 2011. Attosecond laser-material interactions shall also be studied. 4 Levis RJ, Menkir GM, Rabitz H. Selective bond dissociation and rearrangement with (2) Method improvements: In our previous work, by designing and optimally tailored, strong-field laser pulses. Science 2001; 292:709–713. 5 Rezaei S, Li JZ, Herman PR. Burst train generator of high energy femtosecond laser shaping femtosecond laser pulses in temporal/spatial domain(s) to pulses for driving heat accumulation effect during micromachining. Opt Lett 2015; adjust localized transient electron dynamics, the throughput, 40:2064–2067. quality and aspect-ratio limit was greatly improved/extended. 6 Gattass RR, Mazur E. Femtosecond laser micromachining in transparent materials. Nat Photon 2008; 2:219–225. However, we shall: (i) further optimize the experiment parameters 7 Jiang L, Tsai HL. Repeatable nanostructures in dielectrics by femtosecond laser to enhance/improve/extend the fabrication efficiency, quality and pulse trains. Appl Phys Lett 2005; 87: 151104. limit; and (ii) employ these novel fabrication methodologies in 8 Han S, Hong S, Ham J, Yeo J, Lee J et al. Fast plasmonic laser nanowelding for a Cu- nanowire percolation network for flexible transparent conductors and stretchable other fabrication areas, such as three-dimensional laser bio- electronics. Adv Mater 2014; 26: 5808–5814. printing fabrication, green energy fabrication and biomimetic 9 Fang RR, Vorobyev A, Guo CL. Direct visualization of the complete evolution of femtosecond laser-induced surface structural dynamics of metals. Light Sci Appl materials fabrication. 2017; 6: e16256, doi:10.1038/lsa.2016.256. (3) Multiscale measurement system improvements: By using the 10 Fan PX, Bai BF, Long JY, Jiang DF, Jin GF et al. Broadband high-performance infrared current multiscale measurement system, the panoramic dynamics antireflection nanowires facilely grown on ultrafast laser structured Cu surface. Nano Lett 2015; 15:5988–5994. of laser ablation from a femtosecond scale to a second scale have 11 Liu XQ, Yu L, Chen QD, Sun HB. Mask-free construction of three-dimensional silicon been revealed. However, the intrinsic characteristics of a probe structures by dry etching assisted gray-scale femtosecond laser direct writing. Appl beam (wavelength and pulse duration) and the signal collecting Phys Lett 2017; 110: 091602. doi:10.1038/lsa.2017.134 Light: Science & Applications Electrons dynamics control LJiang et al 12 Choi I, Jeong HY, Shin H, Kang G, Byun M et al. Laser-induced phase separation of 48 Hockett P, Bisgaard CZ, Clarkin OJ, Stolow A. Time-resolved imaging of purely silicon carbide. Nat Commun 2016; 7:13562. valence-electron dynamics during a chemical reaction. Nat Phys 2011; 7:612–615. 13 Poumellec B, Lancry M, Desmarchelier R, Hervé E, Bourguignon B. Parity violation in 49 Konar A, Shu YN, Lozovoy VV, Jackson JE, Levine BG et al. Polyatomic molecules chiral structure creation under femtosecond laser irradiation in silica glass? Light Sci under intense femtosecond laser irradiation. J Phys Chem A 2014; 118: Appl 2016; 5: e16178. 11433–11450. 14 Chen YC, Salter PS, Knauer S, Weng LY, Frangeskou AC et al. Laser writing of 50 Sansone G, Kelkensberg F, Pérez-Torres JF, Morales F, Kling MF et al. Electron coherent colour centres in diamond. Nat Photon 2016; 11:77–80. localization following attosecond molecular photoionization. Nature 2010; 465: 15 Wu D, Xu J, Niu LG, Wu SZ, Midorikawa K et al. In-channel integration of designable 763–766. microoptical devices using flat scaffold-supported femtosecond-laser microfabrication 51 Haessler S, Caillat J, Boutu W, Giovanetti-Teixeira C, Ruchon T et al. Attosecond for coupling-free optofluidic cell counting. Light Sci Appl 2015; 4:e228. imaging of molecular electronic wavepackets. Nat Phys 2010; 6:200–206. 16 Li DW, Zhou YS, Huang X, Jiang L, Silvain JF et al. In situ imaging and control of layer- 52 Wang ZH, Zeng B, Li GH, Xie HQ, Chu W et al. Time-resolved shadowgraphs of by-layer femtosecond laser thinning of graphene. Nanoscale 2015; 7:3651–3659. transient plasma induced by spatiotemporally focused femtosecond laser pulses in 17 Chen XD, Xin W, Jiang WS, Liu ZB, Chen YS et al. High-precision twist-controlled fused silica glass. Opt Lett 2015; 40: 5726–5729. bilayer and trilayer graphene. Adv Mater 2016; 28:2563–2570. 53 Höhm S, Rosenfeld A, Krüger J, Bonse J. Femtosecond diffraction dynamics of laser- 18 Coleman C, Erasmus R, Bhattacharyya S. Nanoscale deformations in graphene by induced periodic surface structures on fused silica. Appl Phys Lett 2013; 102: laser annealing. Appl Phys Lett 2016; 109: 253102. 054102. 19 Russo P, Liang R, Jabari E, Marzbanrad E, Toyserkani E et al. Single-step synthesis of 54 Murphy RD, Torralva B, Adams DP, Yalisove SM. Pump-probe imaging of laser- graphene quantum dots by femtosecond laser ablation of graphene oxide dispersions. induced periodic surface structures after ultrafast irradiation of Si. Appl Phys Lett Nanoscale 2016; 8: 8863–8877. 2013; 103: 141104. 20 Kwon HJ, Chung S, Jang J, Grigoropoulos CP. Laser direct writing and inkjet printing 55 Hayasaki Y, Iwata K, Hasegawa S, Takita A, Juodkazis S. Time-resolved axial-view of for a sub-2 μm channel length MoS2 transistor with high-resolution electrodes. the dielectric breakdown under tight focusing in glass. Opt Mater Express 2011; 1: Nanotechnology 2016; 27: 405301. 1399–1408. 21 Dumitru G, Romano V, Weber HP, Sentis M, Marine W. Femtosecond ablation of 56 Papazoglou DC, Tzortzakis S. In-line holography for the characterization of ultrafast ultrahard materials. Appl Phys A 2002; 74:729–739. laser filamentation in transparent media. Appl Phys Lett 2008; 93: 041120. 22 Sun YL, Dong WF, Niu LG, Jiang T, Liu DX et al. Protein-based soft micro-optics fabricated by 57 Temnov VV, Sokolowski-Tinten K, Zhou P, El-Khamhawy A, von der Linde D. femtosecond laser direct writing. Light Sci Appl 2014; 3: e129, doi:10.1038/lsa.2014.10. Multiphoton ionization in dielectrics: comparison of circular and linear polarization. 23 He H, Li SY, Wang SY, Hu ML, Cao YJ et al. Manipulation of cellular light from green Phys Rev Lett 2006; 97: 237403. fluorescent protein by a femtosecond laser. Nat Photon 2012; 6:651–656. 58 Duocastella M, Arnold CB. Bessel and annular beams for materials processing. Laser 24 Park H, Wang X, Nie S, Clinite R, Cao J. Mechanism of coherent acoustic phonon Photon Rev 2012; 6:607–621. generation under nonequilibrium conditions. Phys Rev B 2005; 72: 100301. 59 Papazoglou DG, Suntsov S, Abdollahpour D, Tzortzakis S. Tunable intense Airy beams 25 Juodkazis S, Nishimura K, Tanaka S, Misawa H, Gamaly EG et al. Laser-induced and tailored femtosecond laser filaments. Phys Rev A 2010; 81: 061807. microexplosion confined in the bulk of a sapphire crystal: evidence of multimegabar 60 Mathis A, Courvoisier F, Froehly L, Furfaro L, Jacquot M et al. Micromachining along a pressures. Phys Rev Lett 2006; 96: 166101. curve: femtosecond laser micromachining of curved profiles in diamond and silicon 26 Sundaram SK, Mazur E. Inducing and probing non-thermal transitions in semicon- using accelerating beams. Appl Phys Lett 2012; 101: 071110. ductors using femtosecond laser pulses. Nat Mater 2002; 1:217–224. 61 Žukauskas A, Malinauskas M, Brasselet E. Monolithic generators of pseudo- 27 Sakakura M, Terazima M. Initial temporal and spatial changes of the refractive index nondiffracting optical vortex beams at the microscale. Appl Phys Lett 2013; 103: induced by focused femtosecond pulsed laser irradiation inside a glass. Phys Rev B 181122. 2005; 71: 024113. 62 Bhuyan MK, Velpula PK, Colombier JP, Olivier T, Faure N et al. Single-shot high 28 Zhang N, Zhu XN, Yang JJ, Wang XL, Wang MW. Time-resolved shadowgraphs of aspect ratio bulk nanostructuring of fused silica using chirp-controlled ultrafast laser material ejection in intense femtosecond laser ablation of aluminum. Phys Rev Lett Bessel beams. Appl Phys Lett 2014; 104: 021107. 2007; 99: 167602. 63 Bandres MA, Rodríguez-Lara BM. Nondiffracting accelerating waves: weber waves and 29 Hu HF, Wang XL, Zhai HC, Zhang N, Wang P. Generation of multiple stress waves in silica parabolic momentum. New J Phys 2013; 15: 013054. glass in high fluence femtosecond laser ablation. Appl Phys Lett 2010; 97: 061117. 64 Wetzel B, Xie C, Lacourt PA, Dudley JM, Courvoisier F. Femtosecond laser fabrication 30 Zhao X, Shin YC. Coulomb explosion and early plasma generation during femtosecond of micro and nano-disks in single layer graphene using vortex Bessel beams. Appl Phys laser ablation of silicon at high laser fluence. J Phys D Appl Phys 2013; 46: 335501. Lett 2013; 103: 241111. 31 Domke M, Rapp S, Schmidt M, Huber HP. Ultrafast pump-probe microscopy with high 65 Toyoda K, Miyamoto K, Aoki N, Morita R, Omatsu T. Using optical vortex to control the temporal dynamic range. Opt Express 2012; 20:10330–10338. chirality of twisted metal nanostructures. Nano Lett 2012; 12:3645–3649. 32 Goulielmakis E, Loh ZH, Wirth A, Santra R, Rohringer N et al. Real-time observation of 66 Brixner T, Gerber G. Femtosecond polarization pulse shaping. Opt Lett 2001; 26: valence electron motion. Nature 2010; 466:739–743. 557–559. 33 Smirnova O, Mairesse Y, Patchkovskii S, Dudovich N, Villeneuve D et al. High harmonic 67 Han W, Yang YF, Cheng W, Zhan QW. Vectorial optical field generator for the creation interferometry of multi-electron dynamics in molecules. Nature 2009; 460: 972–977. of arbitrarily complex fields. Opt Express 2013; 21:20692–20706. 34 Improta R, Santoro F, Blancafort L. Quantum mechanical studies on the photophysics 68 Kammel R, Ackermann R, Thomas J, Götte J, Skupin S et al. Enhancing precision in and the photochemistry of nucleic acids and nucleobases. Chem Rev 2016; 116: fs-laser material processing by simultaneous spatial and temporal focusing. Light Sci 3540–3593. Appl 2014; 3: e169, doi:10.1038/lsa.2014.50. 35 Zhang WK, Markiewicz BN, Doerksen RS, Smith AB III, Gai F. C≡N stretching 69 Odoulov S, Shumelyuk A, Badorreck H, Nolte S, Voit KM et al. Interference and vibration of 5-cyanotryptophan as an infrared probe of protein local environment: what holography with femtosecond laser pulses of different colours. Nat Commun 2015; 6: determines its frequency? Phys Chem Chem Phys 2016; 18:7027–7034. 5866. 2+ 36 Auböck G, Chergui M. Sub-50-fs photoinduced spin crossover in [Fe(bpy) ] . Nat 70 Rabitz H, de Vivie-Riedle R, Motzkus M, Kompa K. Whither the future of controlling Chem 2015; 7:629–633. quantum phenomena? Science 2000; 288:824–828. 37 Kraus PM, Tolstikhin OI, Baykusheva D, Rupenyan A, Schneider J et al. Observation of 71 Brif C, Chakrabarti R, Rabitz H. Control of quantum phenomena: past, present, laser-induced electronic structure in oriented polyatomic molecules. Nat Commun and future. New J Phys 2010; 12: 075008. 2015; 6:7039. 72 Colombier J-P, Combis P, Rosenfeld A, Hertel IV, Audouard E et al. Optimized energy 38 Wang HN, Zhang CJ, Rana F. Ultrafast dynamics of defect-assisted electron-hole coupling at ultrafast laser-irradiated metal surfaces by tailoring intensity envelopes: recombination in monolayer MoS . Nano Lett 2014; 15:339–345. Consequences for material removal from Al samples. Phys Rev B 2006; 74: 224106. 39 Pogna EAA, Marsili M, De Fazio D, Dal Conte S, Manzoni C et al. Photo-induced 73 Lindinger A, Lupulescu C, Plewicki M, Vetter F, Merli A et al. Isotope selective bandgap renormalization governs the ultrafast response of single-layer MoS . ACS ionization by optimal control using shaped femtosecond laser pulses. Phys Rev Lett Nano 2016; 10:1182–1188. 2004; 93: 033001. 40 Hong XP, Kim J, Shi SF, Zhang Y, Jin CH et al. Ultrafast charge transfer in atomically 74 Jau YY, Hankin A, Keating T, Deutsch IH, Biedermann G. Entangling atomic spins thin MoS2/WS2 heterostructures. Nat Nanotechnol 2014; 9:682–686. with a Rydberg-dressed spin-flip blockade. Nat Phys 2016; 12:71–74. 41 Zewail AH. Femtochemistry: atomic-scale dynamics of the chemical bond. JPhys 75 King PDC, Wei HI, Nie YF, Uchida M, Adamo C et al. Atomic-scale control of competing Chem A 2000; 104: 5660–5694. electronic phases in ultrathin LaNiO . Nat Nanotechnol 2014; 9: 443–447. 42 Dantus M, Bowman RM, Zewail AH. Femtosecond laser observations of molecular 76 Wienholdt S, Hinzke D, Nowak U. THz switching of antiferromagnets and ferri- vibration and rotation. Nature 1990; 343:737–739. magnets. Phys Rev Lett 2012; 108: 247207. 43 Zewail AH. Laser femtochemistry. Science 1988; 242: 1645–1653. 77 Kampfrath T, Tanaka K, Nelson KA. Resonant and nonresonant control over matter 44 Zewail AH. Femtochemistry: recent progress in studies of dynamics and control of and light by intense terahertz transients. Nat Photon 2013; 7:680–690. reactions and their transition states. J Phys Chem 1996; 100:12701–12724. 78 Mikhaylovskiy RV, Hendry E, Secchi A, Mentink JH, Eckstein M et al. Ultrafast optical 45 Barwick B, Flannigan DJ, Zewail AH. Photon-induced near-field electron microscopy. modification of exchange interactions in iron oxides. Nat Commun 2015; 6: 8190. Nature 2009; 462:902–906. 79 Baierl S, Hohenleutner M, Kampfrath T, Zvezdin AK, Kimel AV et al. Nonlinear spin 46 Carbone F, Kwon OH, Zewail AH. Dynamics of chemical bonding mapped by energy- control by terahertz-driven anisotropy fields. Nat Photon 2016; 10:715–718. resolved 4D electron microscopy. Science 2009; 325:181–184. 80 Kampfrath T, Sell A, Klatt G, Pashkin A, Mahrlein S et al. Coherent terahertz control of 47 Sheng CX, Zhang C, Zhai YX, Mielczarek K, Wang WW et al. Exciton versus free carrier antiferromagnetic spin waves. Nat Photon 2011; 5:31–34. photogeneration in organometal trihalide perovskites probed by broadband ultrafast 81 Renard M, Hertz E, Lavorel B, Faucher O. Controlling ground-state rotational dynamics polarization memory dynamics. Phys Rev Lett 2015; 114: 116601. of molecules by shaped femtosecond laser pulses. Phys Rev A 2004; 69: 043401. Light: Science & Applications doi:10.1038/lsa.2017.134 Electrons dynamics control LJiang et al 82 Moore K, Rabitz H. Laser control: manipulating molecules. Nat Chem 2012; 4: 111 Yuan YP, Jiang L, Li X, Wang C, Lu YF. Adjustment of ablation shapes and 72–73. subwavelength ripples based on electron dynamics control by designing femtosecond 83 Assion A, Baumert T, Bergt M, Brixner T, Kiefer B et al. Control of chemical reactions laser pulse trains. J Appl Phys 2012; 112: 103103. by feedback-optimized phase-shaped femtosecond laser pulses. Science 1998; 282: 112 Yuan YP, Jiang L, Li X, Wang C, Qu LT et al. Simulation of rippled structure 919–922. adjustments based on localized transient electron dynamics control by femtosecond 84 Corrales ME, González-Vázquez J, Balerdi G, Solá IR, De Nalda R et al. Control of laser pulse trains. Appl Phys A 2013; 111:813–819. ultrafast molecular photodissociation by laser-field-induced potentials. Nat Chem 113 Zhang KH, Jiang L, Li X, Shi XS, Yu D et al. Femtosecond laser pulse-train induced 2014; 6:785–790. breakdown in fused silica: the role of seed electrons. J Phys D Appl Phys 2014; 47: 85 Andreev AA, Limpouch J, Iskakov AB, Nakano H. Enhancement of X-ray line emission 435105. from plasmas produced by short high-intensity laser double pulses. Phys Rev E 2002; 114 Jiang L, Tsai HL. Energy transport and material removal in wide bandgap materials by 65: 026403. a femtosecond laser pulse. Int J Heat Mass Transfer 2005; 48:487–499. 86 Fazeli R, Mahdieh MH. Comparison of line X-ray emission from solid and porous nano- 115 Jiang L, Tsai HL. Energy transport and nanostructuring of dielectrics by femtosecond layer coated targets irradiated by double laser pulses. Phys Plasmas 2015; 22: 113303. laser pulse trains. J Heat Transfer 2006; 128:926–933. 87 Fazeli R. Enhanced X-ray emission from laser-produced gold plasma by double pulses 116 Jiang L, Tsai HL. Improved two-temperature model and its application in ultrashort irradiation of nano-porous targets. Phys Lett A 2017; 381:467–471. laser heating of metal films. JHeat Transfer 2005; 127:1167–1173. 88 Esser D, Rezaei S, Li JZ, Herman PR, Gottmann J. Time dynamics of burst-train 117 Jiang L, Tsai HL. Fundamentals of energy cascade during ultrashort laser-material filamentation assisted femtosecond laser machining in glasses. Opt Express 2011; 19: interactions. Proceedings Volume 5713, Photon Processing in Microelectronics and 25632–25642. Photonics IV; 12 April 2005; San Jose, California, United States. San Jose, California, 89 Karimelahi S, Abolghasemi L, Herman PR. Rapid micromachining of high aspect ratio United States: SPTE, 2005. holes in fused silica glass by high repetition rate picosecond laser. Appl Phys A 2014; 118 Jiang L, Tsai HL. A plasma model combined with an improved two-temperature 114:91–111. equation for ultrafast laser ablation of dielectrics. J Appl Phys 2008; 104: 093101. 90 Kerse C, Kalaycıoğlu H, Elahi P, Çetin B, Kesim DK et al. Ablation-cooled material 119 Jiang L, Liu PJ, Yan XL, Leng N, Xu CC et al. High-throughput rear-surface drilling of removal with ultrafast bursts of pulses. Nature 2016; 537:84–88. microchannels in glass based on electron dynamics control using femtosecond 91 Wu AQ, Xu XF. Coherent phonon excitation by ultrafast pulse trains, 2007. First pulse trains. Opt Lett 2012; 37:2781–2783. International Conference on Integration and Commercialization of MICRO and 120 Lin CH, Rao ZH, Jiang L, Tsai WJ, Wu PH et al. Investigations of femtosecond- nanosecond dual-beam laser ablation of dielectrics. Opt Lett 2010; 35:2490–2492. Nanosystems; 10–13 January 2007; Sanya, Hainan, China. Sanya, Hainan, China: 121 Zhao MJ, Hu J, Jiang L, Zhang KH, Liu PJ et al. Controllable high-throughput high- ASME. 2007, pp955–pp957. 92 Wang JL, Guo L, Liu CH, Xu X, Chen YF. Influence of coherent optical phonon on quality femtosecond laser-enhanced chemical etching by temporal pulse shaping ultrafast energy relaxation. Appl Phys Lett 2015; 107: 063107. based on electron density control. Sci Rep 2015; 5:13202. 93 Sheppard CJR, Wilson T. Gaussian-beam theory of lenses with annular aperture. IEE J 122 Jiang L, Shi XS, Li X, Yuan YP, Wang C et al. Subwavelength ripples adjustment based on electron dynamics control by using shaped ultrafast laser pulse trains. Opt Express Microw Opt Acoustics 1978; 2:105–112. 94 Marcinkevičius A, Juodkazis S, Matsuo S, Mizeikis V, Misawa H. Application of Bessel 2012; 20: 21505–21511. beams for microfabrication of dielectrics by femtosecond laser. Jpn J Appl Phys 2001; 123 Shi XS, Jiang L, Li X, Wang SM, Yuan YP et al. Femtosecond laser-induced periodic 40:L1197–L1199. structure adjustments based on electron dynamics control: from subwavelength ripples 95 Courvoisier F, Lacourt PA, Jacquot M, Bhuyan MK, Furfaro L et al. Surface to double-grating structures. Opt Lett 2013; 38:3743–3746. nanoprocessing with nondiffracting femtosecond Bessel beams. Opt Lett 2009; 34: 124 Jiang L, Ying DW, Li X, Lu YF. Two-step femtosecond laser pulse train fabrication of 3163–3165. nanostructured substrates for highly surface-enhanced Raman scattering. Opt Lett 96 Bhuyan MK, Courvoisier F, Lacourt PA, Jacquot M, Salut R et al. High aspect ratio 2012; 37: 3648–3650. nanochannel machining using single shot femtosecond Bessel beams. Appl Phys Lett 125 Zhang N, Li X, Jiang L, Shi XS, Li C et al. Femtosecond double-pulse fabrication of hierarchical nanostructures based on electron dynamics control for high surface- 2010; 97: 081102. 97 Rapp L, Meyer R, Giust R, Furfaro L, Jacquot M et al. High aspect ratio micro- enhanced Raman scattering. Opt Lett 2013; 38:3558–3561. 126 Yang QQ, Li X, Jiang L, Zhang N, Zhang GM et al. Nanopillar arrays with nanoparticles explosions in the bulk of sapphire generated by femtosecond Bessel beams. Sci Rep 2016; 6:34286. fabricated by a femtosecond laser pulse train for highly sensitive SERRS. Opt Lett 98 Wang C, Jiang L, Wang F, Li X, Yuan YP et al. First-principles calculations of the 2015; 40: 2045–2048. 127 Zuo P, Jiang L, Li X, Li B, Xu YD et al. Shape-controllable gold nanoparticles-MoS electron dynamics during femtosecond laser pulse train material interactions. Phys 2 hybrids prepared by tuning edge-active sites and surface structures of MoS via Lett A 2011; 375:3200–3204. 99 Wang C, Jiang L, Wang F, Li X, Yuan YP et al. First-principles electron dynamics temporally shaped femtosecond pulses. ACS Appl Mater Interfaces 2017; 9: control simulation of diamond under femtosecond laser pulse train irradiation. JPhys 7447–7455. 128 Xu CC, Jiang L, Leng N, Yuan YP, Liu PJ et al. Ultrafast laser ablation size and recast Condens Matter 2012; 24: 275801. adjustment in dielectrics based on electron dynamicscontrol by pulse train shaping. 100 Wang C, Jiang L, Wang F, Li X, Yuan YP et al. Transient localized electron dynamics simulation during femtosecond laser tunnel ionization of diamond. Phys Lett A 2012; Chin Opt Lett 2013; 11: 041403. 376: 3327–3331. 129 Wang AD, Jiang L, Li XW, Liu Y, Dong XZ et al. Mask-free patterning of high- 101 Wang C, Jiang L, Li X, Wang F, Yuan YP et al. Frequency dependence of electron conductivity metal nanowires in open air by spatially modulated femtosecond laser pulses. Adv Mater 2015; 27:6238–6243. dynamics during femtosecond laser resonant photoionization of Li cluster. J Appl 130 Yu YW, Jiang L, Cao Q, Shi XS, Wang QS et al. Ultrafast imaging the light-speed Phys 2013; 114: 143105. 102 Wang C, Jiang L, Li X, Wang F, Yuan YP et al. Nonlinear ionization mechanism propagation of a focused femtosecond laser pulse in air and its ionized electron dependence of energy absorption in diamond under femtosecond laser irradiation. dynamics and plasma-induced pulse reshaping. Appl Phys A 2016; 122: 205. 131 Yu YW, Jiang L, Cao Q, Xia B, Wang QS et al. Pump-probe imaging of the fs-ps-ns J Appl Phys 2013; 113: 143106. dynamics during femtosecond laser Bessel beam drilling in PMMA. Opt Express 2015; 103 Yu D, Jiang L, Wang F, Li X, Qu LT et al. Electron ionization and spin polarization control of Fe atom adsorbed graphene irradiated by a femtosecond laser. Phys Lett A 23:32728–32735. 2015; 379:2615–2618. 132 Xie Q, Li XW, Jiang L, Xia B, Yan XL et al. High-aspect-ratio, high-quality microdrilling 104 Yu D, Jiang L, Wang F, Qu LT, Lu YF. First-principles calculation of multiphoton by electron density control using a femtosecond laser Bessel beam. Appl Phys A 2016; 122:136. absorption cross section of α-quartz under femtosecond laser irradiation. Appl Phys A 133 Xia B, Jiang L, Li XW, Yan XL, Lu YF. Mechanism and elimination of bending effect in 2016; 122:494. 105 Su GS, Jiang L, Wang F, Qu LT, Lu YF. First-principles simulations for excitation of femtosecond laser deep-hole drilling. Opt Express 2015; 23:27853–27864. currents in linear carbon chains under femtosecond laser pulse irradiation. Phys Lett A 134 Wang MM, Wang SM, Cao ZT, Wang P, Wang C. Investigation of double-pulse femtosecond laser induced breakdown spectroscopy of polymethyl methacrylate 2016; 380:2453–2457. (PMMA). Proceedings Volume 9351, Laser-based Micro- and Nanoprocessing IX; 12 106 Gao L, Wang F, Jiang L, Qu LT, Lu YF. Controlling the excitation process of free electrons by a femtosecond elliptically polarized laser. Int J Modern Phys B 2015; 29: March 2015; San Francisco, California, United States. San Francisco, California, 1550033. United States: SPIE, 2015. 107 Gao LL, Wang F, Jiang L, Qu LT, Lu YF. Optical-induced electrical current in diamond 135 Zhao WW, Li XW, Xia B, Yan XL, Han WN et al. Single-pulse femtosecond laser Bessel beams drilling of high-aspect-ratio microholes based on electron dynamics control- switched by femtosecond-attosecond laser pulses by ab initio simulations. JPhys D Appl Phys 2015; 49: 025102. Proceedings Volume 9296, International Symposium on Optoelectronic Technology 108 Li X, Jiang L, Tsai H-L. Phase change mechanisms during femtosecond laser pulse and Application 2014: Advanced Display Technology; Nonimaging Optics: Efficient train ablation of nickel thin films. J Appl Phys 2009; 106: 064906. Design for Illumination and Solar Concentration; 21 November 2014; Beijing, China. 109 Li X, Jiang L. Size distribution control of metal nanoparticles using femtosecond Beijing, China: SPIE. 2014. laser pulse train: a molecular dynamics simulation. Appl Phys A 2012; 109: 136 Lorazo P, Lewis LJ, Meunier M. Short-pulse laser ablation of solids: from phase 367–376. explosion to fragmentation. Phys Rev Lett 2003; 91: 225502. 110 YuanYP, Jiang L,Li X,WangC, Xiao H et al. Formation mechanisms of sub- 137 Sokolowski-Tinten K, Solis J, Bialkowski J, Siegel J, Afonso C et al. Dynamics of wavelength ripples during femtosecond laser pulse train processing of dielectrics. ultrafast phase changes in amorphous GeSb films. Phys Rev Lett 1998; 81: J Phys D Appl Phys 2012; 45: 175301. 3679–3682. doi:10.1038/lsa.2017.134 Light: Science & Applications Electrons dynamics control LJiang et al 138 Sokolowski-Tinten K, Bialkowski J, Cavalleri A, von der Linde D, Oparin A et al. 168 Bahk SW, Rousseau P, Planchon TA, Chvykov V, Kalintchenko G et al. Generation and 22 2 Transient states of matter during short pulse laser ablation. Phys Rev Lett 1998; 81: characterization of the highest laser intensities (10 W/cm ). Opt Lett 2004; 29: 224–227. 2837–2839. 139 Girifalco LA, Weizer VG. Application of the Morse potential function to cubic metals. 169 Nathala CSR, Ajami A, Ionin AA, Kudryashov SI, Makarov SV et al. Experimental study Phys Rev 1959; 114:687–690. of fs-laser induced sub-100-nm periodic surface structures on titanium. Opt Express 140 Han WN, Jiang L, Li XW, Wang QS, Li H et al. Anisotropy modulations of femtosecond 2015; 23:5915–5929. laser pulse induced periodic surface structures on silicon by adjusting double 170 Straub M, Afshar M, Feili D, Seidel H, König K. Surface plasmon polariton model of pulse delay. Opt Express 2014; 22: 15820–15828. high-spatial frequency laser-induced periodic surface structure generation in silicon. 141 Leng N, Jiang L, Li X, Xu CC, Liu PJ et al. Femtosecond laser processing of fused silica J Appl Phys 2012; 111: 124315. and aluminum based on electron dynamics control by shaping pulse trains. Appl Phys 171 Derrien TJY, Krüger J, Itina TE, Höhm S, Rosenfeld A et al. Rippled area formed by A 2012; 109:679–684. surface plasmon polaritons upon femtosecond laser double-pulse irradiation of silicon. 142 Preston JS, Van Driel HM, Sipe JE. Order-disorder transitions in the melt morphology Opt Express 2013; 21: 29643–29655. of laser-irradiated silicon. Phys Rev Lett 1987; 58:69–72. 172 Derrien TJY, Krüger J, Itina TE, Höhm S, Rosenfeld A et al. Rippled area formed by 143 Jia TQ, Chen HX, Huang M, Zhao FL, Qiu JR et al. Formation of nanogratings on the surface plasmon polaritons upon femtosecond laser double-pulse irradiation of silicon: surface of a ZnSe crystal irradiated by femtosecond laser pulses. Phys Rev B 2005; the role of carrier generation and relaxation processes. Appl Phys A 2014; 117:77–81. 72: 125429. 173 Raether H. Surface Plasmons on Smooth and Rough Surfaces and on Gratings. 144 Bonse J, Munz M, Sturm H. Structure formation on the surface of indium phosphide Heidelberg: Springer, 1988. irradiated by femtosecond laser pulses. J Appl Phys 2005; 97: 013538. 174 Clark SE, Emmony DC. Ultraviolet-laser-induced periodic surface structures. Phys Rev 145 Okamuro K, Hashida M, Miyasaka Y, Ikuta Y, Tokita S et al. Laser fluence dependence B 1989; 40:2031–2041. of periodic grating structures formed on metal surfaces under femtosecond laser pulse 175 Stoian R, Boyle M, Thoss A, Rosenfeld A, Korn G et al. Laser ablation of dielectrics irradiation. Phys Rev B 2010; 82: 165417. with temporally shaped femtosecond pulses. Appl Phys Lett 2002; 80:353–355. 146 Taylor RS, Hnatovsky C, Simova E, Rayner DM, Bhardwaj VR et al. Femtosecond laser 176 Hommes V, Miclea M, Hergenröder R. Silicon surface morphology study after exposure fabrication of nanostructures in silica glass. Opt Lett 2003; 28: 1043–1045. to tailored femtosecond pulses. Appl Surface Sci 2006; 252: 7449–7460. 147 Shimotsuma Y, Kazansky PG, Qiu JR, Hirao K. Self-Organized nanogratings in glass 177 Kim J, Na S, Cho S, Chang W, Whang K. Surface ripple changes during Cr film irradiated by ultrashort light pulses. Phys Rev Lett 2003; 91: 247405. ablation with a double ultrashort laser pulse. Opt Lasers Eng 2008; 46:306–310. 148 Yuan HC, Yost VE, Page MR, Stradins P, Meier DL et al. Efficient black silicon solar 178 Rosenfeld A, Rohloff M, Höhm S, Krüger J, Bonse J. Formation of laser-induced cell with a density-graded nanoporous surface: optical properties, performance periodic surface structures on fused silica upon multiple parallel polarized double- limitations, and design rules. Appl Phys Lett 2009; 95: 123501. femtosecond-laser-pulse irradiation sequences. Appl Surface Sci 2012; 258: 149 Li JZ, Ho S, Haque M, Herman PR. Nanograting Bragg responses of femtosecond laser 9233–9236. written optical waveguides in fused silica glass. Opt Mater Express 2012; 2: 179 Barberoglou M, Gray D, Magoulakis E, Fotakis C, Loukakos PA et al. Controlling 1562–1570. ripples’ periodicity using temporally delayed femtosecond laser double pulses. Opt 150 Dusser B, Sagan Z, Soder H, Faure N, Colombier JP et al. Controlled nanostructrures Express 2013; 21:18501–18508. formation by ultra fast laser pulses for color marking. Opt Express 2010; 18: 180 Höhm S, Rosenfeld A, Krüger J, Bonse J. Femtosecond laser-induced periodic surface 2913–2924. structures on silica. J Appl Phys 2012; 112: 014901. 151 Vorobyev AY, Guo CL. Colorizing metals with femtosecond laser pulses. Appl Phys Lett 181 Seifert G, Kaempfe M, Syrowatka F, Harnagea C, Hesse D et al. Self-organized 2008; 92: 041914. structure formation on the bottom of femtosecond laser ablation craters in glass. Appl 152 Chen JT, Lai WC, Kao YJ, Yang YY, Sheu JK. Laser-induced periodic structures for Phys A 2005; 81:799–803. light extraction efficiency enhancement of GaN-based light emitting diodes. Opt 182 Ben-Yakar A, Harkin A, Ashmore J, Byer RL, Stone HA. Thermal and fluid processes of Express 2012; 20:5689–5695. a thin melt zone during femtosecond laser ablation of glass: the formation of rims by 153 Zorba V, Stratakis E, Barberoglou M, Spanakis E, Tzanetakis P et al. Biomimetic single laser pulses. J Phys D Appl Phys 2007; 40:1447–1459. artificial surfaces quantitatively reproduce the water repellency of a lotus leaf. Adv 183 Ladieu F, Martin P, Guizard S. Measuring thermal effects in femtosecond laser- Mater 2008; 20:4049–4054. induced breakdown of dielectrics. Appl Phys Lett 2002; 81:957–959. 154 Martín-Fabiani I, Rebollar E, Pérez S, Rueda DR, García-Gutiérrez MC et al. Laser- 184 Richter S, Heinrich M, Döring S, Tünnermann A, Nolte S et al. Nanogratings in fused induced periodic surface structures nanofabricated on poly(trimethylene terephtha- silica: formation, control, and applications. J Laser Appl 2012; 24: 042008. late) spin-coated films. Langmuir 2012; 28:7938–7945. 185 Richter S, Jia F, Heinrich M, Döring S, Peschel U et al. The role of self-trapped 155 Le Harzic R, Dörr D, Sauer D, Stracke F, Zimmermann H. Generation of high spatial excitons and defects in the formation of nanogratings in fused silica. Opt Lett 2012; frequency ripples on silicon under ultrashort laser pulses irradiation. Appl Phys Lett 37:482–484. 2011; 98: 211905. 186 Deng YP, Xie XH, Xiong H, Leng YX, Cheng CF et al. Optical breakdown for silica and 156 Shinoda M, Gattass RR, Mazur E. Femtosecond laser-induced formation of nanometer- silicon with double femtosecond laser pulses. Opt Express 2005; 13:3096–3103. width grooves on synthetic single-crystal diamond surfaces. JAppl Phys 2009; 105: 187 Liang F, Vallée R, Chin SL. Mechanism of nanograting formation on the surface of 053102. fused silica. Opt Express 2012; 20:4389–4396. 157 Bonse J, Krüger J. Pulse number dependence of laser-induced periodic surface structures 188 Li JF, Huang YF, Ding Y, Yang ZL, Li SB et al. Shell-isolated nanoparticle-enhanced for femtosecond laser irradiation of silicon. J Appl Phys 2010; 108: 034903. Raman spectroscopy. Nature 2010; 464:392–395. 158 Bonse J, Rosenfeld A, Krüger J. On the role of surface plasmon polaritons in the 189 De Angelis F, Gentile F, Mecarini F, Das G, Moretti M et al. Breaking the diffusion limit formation of laser-induced periodic surface structures upon irradiation of silicon by with super-hydrophobic delivery of molecules to plasmonic nanofocusing SERS femtosecond-laser pulses. J Appl Phys 2009; 106: 104910. structures. Nat Photon 2011; 5:682–687. 159 Huang M, Zhao FL, Cheng Y, Xu NS, Xu ZZ. Origin of laser-induced near- 190 Lin CH, Jiang L, Chai YH, Xiao H, Chen SJ et al. One-step fabrication of subwavelength ripples: interference between surface plasmons and incident laser. nanostructures by femtosecond laser for surface-enhanced Raman scattering. Opt ACS Nano 2009; 3:4062–4070. Express 2009; 17:21581–21589. 160 Wang L, Xu B-B, Cao X-W, Li Q-K, Tian W-J et al. Competition between subwavelength 191 Nagpal P, Lindquist NC, Oh SH, Norris DJ. Ultrasmooth patterned metals for and deep-subwavelength structures ablated by ultrashort laser pulses. Optica 2017; 4: plasmonics and metamaterials. Science 2009; 325:594–597. 637–642. 192 Le Ru E, Etchegoin P. Principles of Surface-Enhanced Raman Spectroscopy: and 161 Reif J, Costache F, Henyk M, Pandelov SV. Ripples revisited: non-classical morphology Related Plasmonic Effects. Amsterdam: Elsevier, 2008. at the bottom of femtosecond laser ablation craters in transparent dielectrics. Appl 193 Otto A, Mrozek I, Grabhorn H, Akemann W. Surface-enhanced Raman scattering. Surface Sci 2002; 197-198:891–895. J Phys Condens Matter 1992; 4:1143–1212. 162 Le Harzic R, Dörr D, Sauer D, Neumeier M, Epple M et al. Large-area, uniform, high- 194 Paradisanos I, Kymakis E, Fotakis C, Kioseoglou G, Stratakis E. Intense femtosecond spatial-frequency ripples generated on silicon using a nanojoule-femtosecond laser at photoexcitation of bulk and monolayer MoS . Appl Phys Lett 2014; 105: 041108. high repetition rate. Opt Lett 2011; 36:229–231. 195 Shi HY, Yan RS, Bertolazzi S, Brivio J, Gao B et al. Exciton dynamics in suspended 163 Dufft D, Rosenfeld A, Das SK, Grunwald R, Bonse J. Femtosecond laser-induced monolayer and few-layer MoS 2D crystals. ACS Nano 2013; 7:1072–1080. periodic surface structures revisited: a comparative study on ZnO. JAppl Phys 2009; 196 Jersch J, Dickmann K. Nanostructure fabrication using laser field enhancement in the 105: 034908. near field of a scanning tunneling microscope tip. Appl Phys Lett 1996; 68: 164 Borowiec A, Haugen HK. Subwavelength ripple formation on the surfaces of 868–870. 197 Tarun A, Daza MRH, Hayazawa N, Inouye Y, Kawata S. Apertureless optical near-field compound semiconductors irradiated with femtosecond laser pulses. Appl Phys Lett 2003; 82: 4462–4464. fabrication using an atomic force microscope on photoresists. Appl Phys Lett 2002; 165 Li XF, Zhang CY, Li H, Dai QF, Lan S et al. Formation of 100-nm periodic structures 80:3400–3402. on a titanium surface by exploiting the oxidation and third harmonic generation 198 Lee W, Pruzinsky SA, Braun PV. Multi-photon polymerization of waveguide structures induced by femtosecond laser pulses. Opt Express 2014; 22:28086–28099. within three-dimensional photonic crystals. Adv Mater 2002; 14:271–274. 166 Miyaji G, Miyazaki K. Origin of periodicity in nanostructuring on thin film surfaces 199 Liu ZW, Wei QH, Zhang X. Surface plasmon interference nanolithography. Nano Lett ablated with femtosecond laser pulses. Opt Express 2008; 16:16265–16271. 2005; 5:957–961. 167 Hou SS, Huo YY, Xiong PX, Zhang Y, Zhang SA et al. Formation of long- and short- 200 Harilal SS, Diwakar PK, Hassanein A. Electron-ion relaxation time dependent signal periodic nanoripples on stainless steel irradiated by femtosecond laser pulses. JPhys enhancement in ultrafast double-pulse laser-induced breakdown spectroscopy. Appl D Appl Phys 2011; 44: 505401. Phys Lett 2013; 103: 041102. Light: Science & Applications doi:10.1038/lsa.2017.134 Electrons dynamics control LJiang et al 201 Guo J, Wang TF, Shao JF, Sun T, Wang R et al. Emission enhancement ratio of the 207 Döring S, Szilagyi J, Richter S, Zimmermann F, Richardson M et al. Evolution of hole shape and size during short and ultrashort pulse laser deep drilling. Opt Express 2012; metal irradiated by femtosecond double-pulse laser. Opt Commun 2012; 285: 20:27147–27154. 1895–1899. 208 Xia B, Jiang L, Li X, Yan X, Zhao W et al. High aspect ratio, high-quality microholes in 202 Axente E, Noël S, Hermann J, Sentis M, Mihailescu IN. Correlation between plasma PMMA: a comparison between femtosecond laser drilling in air and in vacuum. Appl expansion and damage threshold by femtosecond laser ablation of fused silica. JPhys Phys A 2015; 119:61–68. D Appl Phys 2008; 41: 105216. 203 Quentin U, Leitz KH, Deichmann L, Alexeev I, Schmidt M. Optical trap assisted laser nanostructuring in the near-field of microparticles. J Laser Appl 2012; 24: 042003. This work is licensed under a Creative Commons Attribution 4.0 204 He F, Cheng Y, Xu ZZ, Liao Y, Xu J et al. Direct fabrication of homogeneous International License. The images or other third party material in this microfluidic channels embedded in fused silica using a femtosecond laser. Opt Lett article are included in the article’s Creative Commons license, unless indicated otherwise 2010; 35:282–284. in the credit line; if the material is not included under the Creative Commons license, users 205 Gottmann J, Hermans M, Repiev N, Ortmann J. Selective laser-induced etching of 3D will need to obtain permission from the license holder to reproduce the material. To view a precision quartz glass components for microfluidic applications—up-scaling of copy of this license, visit http://creativecommons.org/licenses/by/4.0/ complexity and speed. Micromachines 2017; 8: 110. 206 Zhao X, Shin YC. Femtosecond laser drilling of high-aspect ratio microchannels in glass. Appl Phys A 2011; 104:713–719. r The Author(s) 2018 Supplementary Information for this article can be found on the Light: Science & Applications’ website (http://www.nature.com/lsa). doi:10.1038/lsa.2017.134 Light: Science & Applications
Light Science & Applications – Springer Journals
Published: Aug 30, 2017
You can share this free article with as many people as you like with the url below! We hope you enjoy this feature!
Read and print from thousands of top scholarly journals.
Already have an account? Log in
Bookmark this article. You can see your Bookmarks on your DeepDyve Library.
To save an article, log in first, or sign up for a DeepDyve account if you don’t already have one.
Copy and paste the desired citation format or use the link below to download a file formatted for EndNote
Access the full text.
Sign up today, get DeepDyve free for 14 days.
All DeepDyve websites use cookies to improve your online experience. They were placed on your computer when you launched this website. You can change your cookie settings through your browser.