Ultrafast laser pulses spatially shaped as Bessel beams in dielectrics create high aspect ratio plasma channels whose relaxation can lead to the formation of nanochannels. We report a strong enhancement of the nanochannel drilling efficiency with illumination by double pulses separated by a delay between 10 and 500 ps. This enables the formation of nanochannels with diameters down to 100 nm. Experimental absorption measurements demonstrate that the increase of drilling efficiency is due to an increase of the confinement of the energy deposition. Nanochannel formation corresponds to a drastic change in absorption of the second pulse, demonstrating the occurrence of a phase change produced by the first pulse. This creates a highly absorbing, long-living state. Our measurements show that it is compatible with the semi-metallization of warm dense glass which takes place within a timescale of <10 ps after the first laser pulse illumination.
Transparent dielectrics are ubiquitous in modern technology, but their structuring at the nanometric scale is a difficult task. The main advantage of ultrafast laser-based processing techniques is the capability of structuring the material in three dimensions as well as processing elongated straight  or curved  voids in single shot. This is possible because of the nonlinearity of the energy deposition process. It allows for a very precise control and versatility of the manufacturing process.
The formation of high aspect ratio voids with diameters in the hundreds of nanometres range using only a single ultrafast laser pulse is particularly useful for the creation of nanophotonic crystals and Bragg gratings , [, 3] or for glass cutting, in the framework of so called “stealth-dicing” technique , , , , , ,  which enables cutting glass using laser illumination speeds exceeding tens of centimetres per second.
Drilling of high aspect ratio nanochannels is made particularly controllable when using zeroth-order Bessel beams , , , , , . For sufficiently high focussing angles, Kerr effect is negligible and energy absorption occurs mainly in the central lobe. Although the mechanism leading to void channel formation in transparent dielectrics after the energy deposition stage is still an open question , , , our problematic here is that depending on the material and illumination geometry, limitations arise on the maximal and minimal nanochannel diameters that can be processed with Bessel beams. This is exemplified in Figure 1(a). It shows scanning electron microscopy (SEM) images of nanochannels processed in Schott D263 glass with single femtosecond pulses. The pulse energy range over which a nanochannel can be opened is very limited. Increasing the pulse energy did not lead to an increase of the channel diameter above ∼300 nm, in contrast with a relatively similar glass, Corning 0211 where nanochannels could be processed with diameters from 200 to 800 nm  with a twice shorter Bessel beam length.
Here, we investigate how double pulses can solve this issue. Surface ablation by double Gaussian pulses conventionally sees a decrease of the ablation efficiency. This has mainly been attributed to a reduction of nonlinear ionization due to the intensity decrease when splitting the pulse in two, as well as to a screening effect from the plasma , , , , . In contrast, the case of in-volume material excitation, at fluences lower than the ablation threshold, shows stronger modifications of fused silica (index of refraction, nanogratings formation) with double pulses in comparison with single pulses of identical total energy, for inter-pulse delay typically on the order of 10–100 ps , , , , , . The stronger modifications were attributed to an increase of the total deposited energy because of a higher absorption of the second pulse by the material softened after the first pulse. The first pulse generates a plasma of free-electrons and holes or ions that decays into self-trapped excitons and colour centres. Those defects increase the absorption efficiency of the second pulse.
In this article, we report that splitting the pulse energy in two equal energy pulses allows for drastically increasing the nanochannel formation efficiency. Nanochannels with ∼100 nm diameters can be reliably processed, and we can increase by a factor two the maximal channel diameter compared to single shot case. We have studied the evolution of the channel morphology with inter-pulse delay and have characterized the evolution of absorption as a function of delay and energy. The overall absorption is slightly less for the double pulse case than for the single pulse one. But our results demonstrate for the first time that pulse splitting with a delay in the 100–500 ps allows for enhancing the confinement of energy deposition.
An analysis of the evolution of the second pulse absorption in time reveals two different dynamics respectively for sub- and above threshold for nanochannel drilling. While at the lowest energies, the absorption dynamics follows the decrease conventionally expected by the formation of self-trapped excitons, our results at higher energies, those which allow nanochannel formation, show a striking new behaviour for bulk excitation. These are compatible with the generation of warm dense silica within ∼10 ps after the first pulse, which has a very localized and high absorptivity because of semi-metallization. The Bessel beam geometry associated to a pulse burst in the GHz repetition rate regime  is therefore expected to lead to highly efficient formation of Warm Dense Matter with tabletop experiments, which is important for the fundamental understanding of astrophysical objects .
Our experimental setup is based on a 5 kHz ultrafast Ti-sapphire laser source with 800 nm central wavelength. Before beam shaping, the ultrafast pulses were split into equal parts and recombined using a Mach–Zehnder interferometer arrangement. The pulses were independently circularly polarized with the same rotation orientation. The inter-pulse delay is controlled using a micrometric motorized translation stage and can be varied up to a maximum of 500 ps. The zero-delay was characterized using the interference between the pulses. Beam shaping was performed using an axicon with base angle 1° and a telescopic arrangement with magnification ∼110 to generate a Bessel beam with a length of ∼100 μm and a cone angle of 26° in air.
A computer-controlled Pockels cell was used for single shot illumination. To avoid any potential cumulative effects due to the irradiation by the small polarization leakage during the dead time between the sample positioning to a fresh area and the effective illumination by the high intensity pulses, two additional fast mechanical shutters were installed with effective opening time of ∼2 ms. In practice, we have observed that the presence of these fast shutters had no effect in the irradiations, ensuring that the leaks from the Pockels cell are absolutely negligible.
The pulse energy and duration were characterized at the sample site. The pulse duration is 100 ± 15 fs. We used recently cleaved Schott D236 borosilicate glass slides as sample. The sample flatness was ensured to be below 0.1 mrad. Beam characterization was performed using a microscope objective with numerical aperture 0.8, which was combined to a lens to form a telescope with a magnification factor of 56 . Figure 2 shows the intensity distribution of the generated beams. We highlight the fact that the high magnification (×110) used for beam shaping also magnifies any angular deviation from a misalignment in the Mach–Zehnder interferometer. Even the small flatness deviation of the high-precision translation stage generates a deviation by more than 1 µrad when translating the mirrors to change the inter-pulse delay from 0 to 100 ps. This deviation is enough to prevent the accurate superposition of the central lobes of the two Bessel beams. The bottom panel in Figure 2 shows the perfect beam superposition, which was corrected for each inter-pulse delay.
The samples were characterized first using conventional optical microscopy with through-illumination using a microscope objective equipped with a variable correction collar. This was followed by Focused Ion Beam (FIB) milling where a careful procedure was developed to avoid modifying the laser processed structures as in reference . This allows side-imaging the nanochannels by SEM.
In a second step, we conducted absorption measurements of the single pump pulse or of the double pump one. Pulse absorption was characterized using the ratio of the signals between two photodiodes positioned before and after the sample. The photodiodes were large photodiodes (100 mm2) to ensure the whole beam was collected. The reference photodiode was placed after a 10% splitter before the first telescope. After the sample, the beam was collimated by a microscope objective with numerical aperture 0.8 and the far-field was imaged onto the second photodiode. The acquisition was performed in single shot using an oscilloscope. The linearity of the acquisition was specifically ensured over the whole energy range investigated. For the absorption measurements, the Bessel beam was fully enclosed within the sample thickness, and the sample was also translated so that each shot illuminates an undamaged area.
Because of the Bessel beam conical structure, reflection and transmission on the plasma are indistinguishable. Both transmitted and reflected beams are directed in the forward direction with the same angle with respect to the optical axis. The overall absorption A is therefore complementary to the measured transmission: A = 1 − E/E0, where E and E0 are the energies collected after the sample with the Bessel beam, respectively, placed inside and outside of the glass sample. Therefore, the Fresnel losses are eliminated from the absorption curves. The experiment was repeated 50 times with the same parameters and the results were averaged.
Figure 1(a) shows SEM images of Schott D263 glass samples after single pulse irradiation. We observe that channels of diameters ∼200 and 300 nm are drilled for energies of, respectively, 3 and 5.2 μJ. However, further increase of the pulse energy does not lead to channel opening but instead to a density modification. This modification extends over a diameter which increases with pulse energy. We infer that at the highest energies shown here, a plasma is formed over a too wide diameter, which limits the thermodynamical gradients needed to generate microexplosions and void formation.
We repeated the same experiment for double pulses with inter-pulse delay of 500 ps and the results are shown in Figure 1(b). We readily observe that the morphology of the channels is relatively similar to the single pulse case, except that, now, nanochannels are opened in all cases. The produced channels’ width extends from ∼100 to more than 700 nm. This is a much wider processing window than in the single pulse case. The small inhomogeneities along the channels are attributed to the intensity inhomogeneities of the beam as can be seen from Figure 2 .
We will now investigate the influence of the inter-pulse delay. For this, we show in Figure 3 images of nanochannels under optical microscopy with five inter-pulse delays. For each case, the drilling experiment has been repeated several times, so as to evaluate the repeatability of the laser-induced material modification. In the figure, we show two images of channels obtained in identical conditions. We note that the distance between two channels is 20 μm to avoid any influence of a material modification by a previous pulse on the next modifications. Under our optical characterization system, voids with diameter ⪆400 nm appear deep black.
Delays below the pulse duration produce results that vary very significantly from shot to shot (not shown). A possible explanation is that the phase between the two pulses cannot be controlled with sufficient accuracy to maintain the same interference state from shot to shot. Strong shot to shot variations were similarly observed in the absorption experiments as described below.
For inter-pulse delays from 1 ps and up to 500 ps, irradiations produced channels over the full tested range of energies (2–16.5 μJ total energy), in high contrast with the single pulse case. This is confirmed from SEM measurements. For inter-pulse delays between 0.2 and ∼10 ps, a transition occurs where the channels appear wider, more uniform along their length and more reproducible on a shot to shot basis. For delays of 10 ps and above, the structures show a higher degree of uniformity and repeatability. We noted that the optimum depends on the total pulse energy. In the case where the inter-pulse delay was of several tens of milliseconds, i.e. when the excitation by the first pulse has fully relaxed before the second pulse arrives (noted as ∞ in the figure), the structures are much more faint and less regular. This indicates the state reached at 500 ps is not identical to state reached after the complete relaxation.
Noticeably, the marks of thermal effects are highly reduced on the optical microscopy images and on the SEM images of Figure 1(b) in comparison with what is conventionally observed with picosecond illumination in glass  or in sapphire , or in comparison with illumination via a second pulse with nanosecond duration or more , [, 35]. We infer that the thermodynamical pathways (ionization dynamics, pressure, temperature and so on) are different between the case of double femtosecond pulses excitation and the case of illumination by picosecond pulses.
Figure 4 quantitatively compares the evolution of the channel diameter with the total input pulse energy for single pulses and three delays. The measurements have been performed using SEM images of the channel cross-sections, after FIB milling. The black markers show the evolution of the diameter of the modification produced in the single shot regime. The error bar on the measurement is due to the accuracy of FIB milling between two SEM images and is ∼40 nm. Voids are shown with black disks, density modifications with the circles. The three other lines show the results for inter-pulse delays of 1, 10 and 500 ps. We first notice that for the three cases in double pulse regime, the evolution is approximately the same, with a quasi-linear increase with pulse energy. This follows the same trend as shown originally in reference . At the lowest energies, the longest delays produce the smallest diameters. This variation can be attributed to a decay of the material excitation by the first pulse. For a delay of 1 ps, at energies above 8 μJ, we observe a fluctuation of the diameter. This variation can be attributed to the strong phase transition that we describe in the next section, which occurs over a scale of 1–10 ps depending on the deposited energy.
The differences between single and double pulses are noticeable. First, double pulses allow reaching extremely small diameters down the ∼100 nm, i.e. nearly eight times lower than the laser wavelength. For the glass used here (Schott D263), the maximal diameter reached with double pulses is almost three times higher than in the single pulse case. These results suggest that wider channels could be produced using higher pulse energies. We expect that the maximal channel diameter achievable will depend on the inter-pulse delay.
A priori, two origins can explain a higher drilling efficiency with double pulses than with single ones. This can be due to a higher overall absorption of energy or to the absorption of a possibly smaller amount of energy but on a significantly smaller diameter (the length of the modifications remains constant between single and double pulse illumination). We demonstrate here that it is the second case.
Figure 5(a) shows the total absorption for single and double pulses as a function of the total input pulse energy. It can be seen that the absorption is overall similar in all cases, thus, splitting the pulse does not lead to a higher absorption. The absorption is in fact higher for the single pulse case (black curve), particularly for energies in the range 2–4 μJ, where the slope is the highest and which also correspond to the threshold energy for channel drilling. Therefore, the highest drilling efficiency observed at 500 ps inter-pulse delay does not correspond to a higher absorption. The higher absorption measured for the single pulse case can be attributed to a higher efficiency for the multiphoton and tunnel ionization processes. Those might take place either out of the central lobe of the Bessel beam or within this central lobe, but on a wider area than in the double pulse case. We conclude that the enhanced efficiency for void channel opening originates from a higher density of deposited energy.
The absorption of the second pulse can be straightforwardly derived from the overall absorption by subtracting the single pulse contribution. Figure 5(b) shows the absorption of the first and second pulses, shown as a function of the individual pulse energy (we recall that first and second pulses have the same energy). Above an individual pulse energy of ∼1 μJ, the second pulse absorption drastically increases with a step-like curve, reaching absorption above 0.6, and well above the first pulse absorption at this energy. This is followed by a smoother increase, for energies greater than 2 μJ. We note that the step in absorption occurs for an energy which nearly corresponds to the threshold for nanochannel drilling (2 μJ total energy). The abrupt change in absorption shows that an important transition is crossed when the energy of the first pulse is sufficiently high. We will see below that this transition corresponds to the formation of warm dense matter.
To obtain a better view of the evolution of absorption in time, we plot in Figure 6 the evolution of the second pulse absorption as a function of the inter-pulse delay for a set of representative energies. Square markers show the first pulse absorption as a reference, for the same pulse energy.
For the lowest individual pulse energies, typically below 1 μJ (i.e. below threshold for nanochannel drilling in double pulse regime) the absorption effectively decreases with the inter-pulse delay (red and blue curves). This is consistent with the decay times of self-trapped excitons (STEs), which are on the order of 30 and 300 ps . The higher absorption of the second pulse in comparison with the first can be attributed in this case to the formation of STEs which have a much higher ionization cross-section than glass because of their intermediate position in the bandgap. The evolution of the second pulse absorption with time at low pulse energies explain the dynamics observed by Wang et al. using Bessel beams in fused silica in double pulse configuration: after chemical etching of the modifications inscribed, channel length highly increased for delays in the pulse separation window 1–50 ps and decreased for longer delays .
In contrast, at energies above the threshold for channel drilling (e.g. green and purple curves for 1.5 and 3 μJ), the absorption curves show only a small decrease with delay up to 10 ps and then stay nearly constant with time. We also note that the curves with different delays are nearly identical within the fluctuations in Figure 5(b) above 1.5 μJ. This indicates that, in the case of nanochannel drilling, a transformation triggered by the first pulse occurs during 10 ps and reaches a state which has a typical lifetime exceeding the nanosecond scale.
Our observations are in high agreement with the results by Garcia-Lechuga et al. of the measurements of the transient reflectivity and transmission of the surface of fused silica after irradiation by high fluence (7 J/cm2) femtosecond pulses . A state with absorptivity close to 1 and negligible reflectivity was characterized to form within the same scale of 1–10 ps. This state of matter was inferred by the authors to be a hot blackbody. In our case, the fluence reached for only 1 μJ single pulse energy exceeds 10 J/cm2 in the central lobe of the Bessel beam because of the extreme focussing angle (without accounting for potential non-linear losses on the beam path). The fluence range in which we operate is therefore very similar to the conditions of the work mentioned above. We note that in our case, the pulse absorption cannot reach 1 because the plasma, with typical size of a few hundreds of nanometres, is smaller than the diameter of the Bessel beam central spot (660 nm FWHM, 1.2 μm in linear regime), so that the interaction cross section between the second pulse and the hot cylinder of matter is necessarily below 1.
For an individual pulse energy of 1.5 μJ, at the onset of the plateau of absorption in Figure 5(a), the energy density absorbed within the central lobe is already exceeding 8 MJ/kg, for a diameter of the excited material estimated to ≲ 600 nm (order of magnitude of the central lobe FWHM of the Bessel beam). This corresponds to the typical energy density of Warm Dense Matter in fused silica . Denoeud et al. have experimentally determined with X-ray absorption near-edge spectroscopy that warm dense glass shows a semi-metallic behaviour . Engelhorn et al. demonstrated that with similar energy density (10 MJ/kg) as in our experiment, warm dense silicon dioxide reaches an ionic temperature of 5000 K and an electronic temperature on the order of 30,000 K within a duration of 20 ps . Our results are fully in line with those, including the typical temporal scale of the transformation.
Therefore, our results can be explained by the metallization of glass occurring within a couple tens of picosecond after the first pulse. Because of the drastic change in drilling efficiency observed between single and double pulses, absorption of the second pulse must occur within a significantly smaller volume than in the case of a single pulse. We interpret the increase of reproducibility as originating from a progressive increase of the confinement of absorption between 10 and 500 ps. The warm dense matter diameter shrinks while its absorption is expected to increase. The diameter shrinkage is associated to some energy loss, which makes, at the lowest energies, the single pulse slightly more efficient than the double pulse with 500 ps, but less efficient than with 1 and 10 ps delays (see Figure 4). We note that modelling the absorption properties of warm dense matter is for now computationally extremely challenging . The large variability observed at delays smaller than 10 ps is attributed to the phase change occurring during this time.
We have demonstrated that double pulse femtosecond illumination can drastically enhance the efficiency of nanochannel formation in comparison with the single pulse case. The morphology of the channels shows a negligible heat affected zone in contrast with conventional illumination with picosecond durations. Absorption measurements show that the increased void channel formation efficiency can be understood as a higher confinement of the energy deposition. At energies below the threshold for nanochannel opening, the absorption of the second pulse decreases with time as can be expected from the formation of self-trapped excitons. At higher energies, where nanochannel formation occurs, our experimental results indicate that the first pulse triggers a phase transition within 1–10 ps, and this state has a lifetime exceeding 1 ns.
The higher efficiency for void formation is explained by the swift generation of the semi-metallic state of warm dense glass generated over a diameter significantly smaller than in the single pulse case, with a diameter typically lower than 600 nm. For all delays above 1 ps, channel formation was observed in double pulse illumination, and a smoother transition operates between 10 and several hundreds of picoseconds. We note that our measurements cannot distinguish between the generation by the first pulse of a cylinder with a uniform density of warm dense matter or if a rarefied region at the centre, surrounded by warm dense matter progressively forms at a scale of several hundreds of picoseconds. Our results are therefore compatible with both expected mechanisms of nanochannel formation , [, 16].
Our results have implications both on applied and fundamental aspects. In the field of nano-machining with ultrafast laser pulses, the best configuration in the investigated delay range is for 10–500 ps depending on the energy. This time interval falls in the GHz repetition range. Our results can also explain the enhancement of localized energy deposition in glass using GHz repetition rate sources . The multiple pulse configuration, at high total energies, can be used to increase the maximal void diameters, while at the lowest energies, fine tuning of the void diameter below ∼100 nm becomes possible. Further research with GHz sources is needed to evaluate an optimal number of pulses. Our results also demonstrate the capability of generating Warm Dense Matter within the bulk of materials by using burst of pulses. We believe a similar state can be created within other transparent materials such as sapphire. The confinement in the material’s bulk using a non-diffracting Bessel beam increases the confinement in comparison with laser impacts on surface. This opens new routes to create and control high temperature and high pressure states within relatively high volumes, with length up to centimetre scale , for tabletop laboratory scale experiments to reproduce astrophysical pressure and temperature conditions.
We acknowledge the support and funding from European Research Council (ERC) 682032-PULSAR, Region Franche-Comte council, the EIPHI Graduate School (ANR-17-EURE-0002), I-SITE BFC (ANR-15-IDEX-0003), French RENATECH network.
Technical assistance by C. Billet, E. Dordor and fruitful discussions with R. Giust are gratefully acknowledged.
 M. K. Bhuyan, F. Courvoisier, P. A. Lacourt, et al.., “High aspect ratio nanochannel machining using single shot femtosecond Bessel beams,” Appl. Phys. Lett., vol. 97, p. 081102, 2010, https://doi.org/10.1063/1.3479419. Search in Google Scholar
 L. Froehly, F. Courvoisier, A. Mathis, et al.., “Arbitrary accelerating micron-scale caustic beams in two and three dimensions,” Opt. Express, vol. 19, p. 16455, 2011, https://doi.org/10.1364/oe.19.016455. Search in Google Scholar
 M. Mikutis, T. Kudrius, G. Šlekys, D. Paipulas, and S. Juodkazis, “High 90% efficiency Bragg gratings formed in fused silica by femtosecond Gauss-Bessel laser beams,” Opt. Mater. Express, vol. 3, p. 862, 2013. Search in Google Scholar
 F. Ahmed, M. Lee, H. Sekita, T. Sumiyoshi, and M. Kamata, “Display glass cutting by femtosecond laser induced single shot periodic void array,” Appl. Phys. A, vol. 93, pp. 189–192, 2008, https://doi.org/10.1007/s00339-008-4672-2. Search in Google Scholar
 W.-J. Tsai, C.-J. Gu, C.-W. Cheng, and J.-B. Horng, “Internal modification for cutting transparent glass using femtosecond Bessel beams,” Opt. Eng., vol. 53, pp. 1–7, 2013, https://doi.org/10.1117/1.oe.53.5.051503. Search in Google Scholar
 M. K. Bhuyan, O. Jedrkiewicz, V. Sabonis, et al.., “High-speed laser-assisted cutting of strong transparent materials using picosecond Bessel beams,” Appl. Phys. A, vol. 120, pp. 443–446, 2015, https://doi.org/10.1007/s00339-015-9289-7. Search in Google Scholar
 K. Mishchik, C. J. Leger, O. D. Caulier, et al.., “Ultrashort pulse laser cutting of glass by controlled fracture propagation,” J. Laser Micro Nanoen., vol. 11, pp. 66–70, 2016, https://doi.org/10.2961/jlmn.2016.01.0012. Search in Google Scholar
 R. Meyer, R. Giust, M. Jacquot, J. M. Dudley, and F. Courvoisier, “Submicron-quality cleaving of glass with elliptical ultrafast Bessel beams,” Appl. Phys. Lett., vol. 111, p. 231108, 2017, https://doi.org/10.1063/1.5008921. Search in Google Scholar
 M. Jenne, D. Flamm, T. Ouaj, et al.., “High-quality tailored-edge cleaving using aberration-corrected Bessel-like beams,” Opt. Lett., vol. 43, p. 3164, 2018, https://doi.org/10.1364/ol.43.003164. Search in Google Scholar
 R. Meyer, L. Froehly, R. Giust, et al.., “Extremely high-aspect-ratio ultrafast Bessel beam generation and stealth dicing of multi-millimeter thick glass,” Appl. Phys. Lett., vol. 114, p. 201105, 2019, https://doi.org/10.1063/1.5096868. Search in Google Scholar
 I. Alexeev, K. H. Leitz, A. Otto, and M. Schmidt, “Application of Bessel beams for ultrafast laser volume structuring of non transparent media,” Phys. Proc., vol. 5, pp. 533–540, 2010, https://doi.org/10.1016/j.phpro.2010.08.177. Search in Google Scholar
 P. Boucher, J. Del Hoyo, C. Billet, O. Pinel, G. Labroille, and F. Courvoisier, “Generation of high conical angle Bessel–Gauss beams with reflective axicons,” Appl. Opt., vol. 57, p. 6725, 2018, https://doi.org/10.1364/ao.57.006725. Search in Google Scholar
 E. N. Glezer and E. Mazur, “Ultrafast-laser driven micro-explosions in transparent materials,” Appl. Phys. Lett., vol. 71, p. 882, 1997, https://doi.org/10.1063/1.119677. Search in Google Scholar
 A. Vailionis, E. G. Gamaly, V. Mizeikis, W. Yang, A. V. Rode, and S. Juodkazis, “Evidence of superdense aluminium synthesized by ultrafast microexplosion,” Nat. Commun., vol. 2, p. 445, 2011, https://doi.org/10.1038/ncomms1449. Search in Google Scholar
 M. K. Bhuyan, M. Somayaji, A. Mermillod-Blondin, et al.., “Ultrafast laser nanostructuring in bulk silica, a “slow” microexplosion,” Optica, vol. 4, p. 951, 2017, https://doi.org/10.1364/optica.4.000951. Search in Google Scholar
 R. Stoian, M. Boyle, A. Thoss, et al.., “Laser ablation of dielectrics with temporally shaped femtosecond pulses,” Appl. Phys. Lett., vol. 80, pp. 353–355, 2002, https://doi.org/10.1063/1.1432747. Search in Google Scholar
 A. Semerok and C. Dutouquet, “Ultrashort double pulse laser ablation of metals,” Thin Solid Films, vols 453–454, pp. 501–505, 2004, https://doi.org/10.1016/j.tsf.2003.11.115. Search in Google Scholar
 I. H. Chowdhury, X. Xu, and A. M. Weiner, “Ultrafast double-pulse ablation of fused silica,” Appl. Phys. Lett., vol. 86, p. 151110, 2005, https://doi.org/10.1063/1.1901806. Search in Google Scholar
 Z. Cao, L. Jiang, S. Wang, et al.., “Influence of electron dynamics on the enhancement of double-pulse femtosecond laser-induced breakdown spectroscopy of fused silica,” Spectrochim. Acta B, vol. 141, pp. 63–69, 2018, https://doi.org/10.1016/j.sab.2018.01.010. Search in Google Scholar
 K. Gaudfrin, J. Lopez, K. Mishchik, L. Gemini, R. Kling, and G. Duchateau, “Fused silica ablation by double femtosecond laser pulses: influence of polarization state,” Opt. Express, vol. 28, p. 15189, 2020, https://doi.org/10.1364/oe.387803. Search in Google Scholar
 T. Nagata, M. Kamata, and M. Obara, “Optical waveguide fabrication with double pulse femtosecond lasers,” Appl. Phys. Lett., vol. 86, p. 251103, 2005, https://doi.org/10.1063/1.1954879. Search in Google Scholar
 D. Wortmann, M. Ramme, and J. Gottmann, “Refractive index modification using fs-laser double pulses,” Opt. Express, vol. 15, p. 10149, 2007, https://doi.org/10.1364/oe.15.010149. Search in Google Scholar
 D. Chu, X. Sun, X. Dong, et al.., “Effect of double-pulse-laser polarization and time delay on laser-assisted etching of fused silica,” J. Phys. D Appl. Phys., vol. 50, p. 465306, 2017, https://doi.org/10.1088/1361-6463/aa8e75. Search in Google Scholar
 Z. Wang, L. Jiang, X. Li, et al.., “High-throughput microchannel fabrication in fused silica by temporally shaped femtosecond laser Bessel-beam-assisted chemical etching,” Opt. Lett., vol. 43, p. 98, 2017, https://doi.org/10.1364/ol.43.000098. Search in Google Scholar
 H. Wang, J. Song, Q. Li, X. Zeng, and Y. Dai, “Formation of nanograting in fused silica by temporally delayed femtosecond double-pulse irradiation,” J. Phys. D Appl. Phys., vol. 51, p. 155101, 2018, https://doi.org/10.1088/1361-6463/aab376. Search in Google Scholar
 V. Stankevič, J. Karosas, G. Račiukaitis, and P. Gečys, “Improvement of etching anisotropy in fused silica by double-pulse fabrication,” Micromachines, vol. 11, p. 483, 2020. Search in Google Scholar
 C. Kerse, H. Kalaycıoğlu, P. Elahi, et al.., “Ablation-cooled material removal with ultrafast bursts of pulses,” Nature, vol. 537, pp. 84–88, 2016, https://doi.org/10.1038/nature18619. Search in Google Scholar
 M. Bonitz, T. Dornheim, Z. A. Moldabekov, et al.., “Ab initio simulation of warm dense matter,” Phys. Plasmas, vol. 27, p. 042710, 2020, https://doi.org/10.1063/1.5143225. Search in Google Scholar
 L. Froehly, M. Jacquot, P. A. Lacourt, J. M. Dudley, and F. Courvoisier, “Spatiotemporal structure of femtosecond Bessel beams from spatial light modulators,” J. Opt. Soc. Am. A, vol. 31, p. 790, 2014, https://doi.org/10.1364/josaa.31.000790. Search in Google Scholar
 L. Rapp, R. Meyer, R. Giust, et al.., “High aspect ratio micro-explosions in the bulk of sapphire generated by femtosecond Bessel beams,” Sci. Rep., vol. 6, p. 34286, 2016, https://doi.org/10.1038/srep34286. Search in Google Scholar
 J. Dudutis, R. Stonys, G. Račiukaitis, and P. Gečys, “Aberration-controlled Bessel beam processing of glass,” Opt. Express, vol. 26, p. 3627, 2018, https://doi.org/10.1364/oe.26.003627. Search in Google Scholar
 N. Götte, T. Winkler, T. Meinl, et al.., “Temporal airy pulses for controlled high aspect ratio nanomachining of dielectrics,” Optica, vol. 3, p. 389, 2016, https://doi.org/10.1364/optica.3.000389. Search in Google Scholar
 C.-H. Lin, Z. H. Rao, L. Jiang, et al.., “Investigations of femtosecond–nanosecond dual-beam laser ablation of dielectrics,” Opt. Lett., vol. 35, p. 2490, 2010, https://doi.org/10.1364/ol.35.002490. Search in Google Scholar
 Y. Ito, R. Yoshizaki, N. Miyamoto, and N. Sugita, “Ultrafast and precision drilling of glass by selective absorption of fiber-laser pulse into femtosecond-laser-induced filament,” Appl. Phys. Lett., vol. 113, p. 061101, 2018, https://doi.org/10.1063/1.5027421. Search in Google Scholar
 D. Grojo, M. Gertsvolf, S. Lei, T. Barillot, D. M. Rayner, and P. B. Corkum, “Exciton-seeded multiphoton ionization in bulk SiO2,” Phys. Rev. B, vol. 81, 2010, https://doi.org/10.1103/physrevb.81.212301. Search in Google Scholar
 M. Garcia-Lechuga, L. Haahr-Lillevang, J. Siegel, P. Balling, S. Guizard, and J. Solis, “Simultaneous time-space resolved reflectivity and interferometric measurements of dielectrics excited with femtosecond laser pulses,” Phys. Rev. B, vol. 95, 2017, https://doi.org/10.1103/physrevb.95.214114. Search in Google Scholar
 K. Engelhorn, V. Recoules, B. I. Cho, et al.., “Electronic structure of warm dense silicon dioxide,” Phys. Rev. B, vol. 91, 2015, https://doi.org/10.1103/physrevb.91.214305. Search in Google Scholar
 A. Denoeud, A. Benuzzi-Mounaix, A. Ravasio, et al.., “Metallization of warm dense SiO2 studied by XANES spectroscopy,” Phys. Rev. Lett., vol. 113, 2014, https://doi.org/10.1103/physrevlett.113.116404. Search in Google Scholar