# Long distance laser filamentation using Yb:YAG kHz laser

Pierre Walch, Benoît Mahieu, Victor Moreno, Thomas Produit, Ugo Andral, Yves-Bernard André, Laurent Bizet, Magali Lozano, Clemens Herkommer, Michel Moret, et al.

► To cite this version:

Pierre Walch, Benoît Mahieu, Victor Moreno, Thomas Produit, Ugo Andral, et al. Long distance laser filamentation using Yb:YAG kHz laser. Scientific Reports, 2023, 13 (1), pp.18542. 10.1038/s41598-023-45660-9. hal-04264176

HAL Id: hal-04264176

https://hal.science/hal-04264176v1

Submitted on 30 Oct 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

# Long distance laser filamentation using Yb:YAG kHz laser

Pierre Walch<fnref n="1" />, Benoît Mahieu<fnref n="1" />, Victor Moreno², Thomas Produit², Ugo Andral², Yves-Bernard André<fnref n="1" />, Laurent Bizet<fnref n="1" />, Magali Lozano<fnref n="1" />, Clemens Herkommer³, Michel Moret², Robert Jung³, Robert Bessing³, Sandro Klingebiel³, Yann Bertho⁴, Thomas Metzger³, André Mysyrowicz<fnref n="1" />,⁵, Jean-Pierre Wolf², Jérôme Kasparian²,⁶ & Aurélien Houard<fnref n="1" />

In the framework of the Laser Lightning Rod project, whose aim is to show that laser-induced filaments can guide lightning discharges over considerable distances, we study over a distance of 140 m the filaments created by a laser system with J-range pulses of 1 ps duration at 1 kHz repetition rate. We investigate the spatial evolution of the multiple filamentation regime using the fundamental beam at 1030 nm or using combination with the second and third harmonics. The measurements were made using both a collimated beam and a loosely focused beam.

Laser filamentation has been observed in condensed matter only few years after the invention of the laser<fnref n="1" />. The advent of the chirped pulse amplification strongly boosted the peak power delivered by lasers from the megawatt to the gigawatt², and soon to the TW level, allowing the observation of laser filaments in air³, where the nonlinear refractive index is three orders of magnitude lower. While in condensed matter filamentation is deemed a detrimental effect causing damage to optical components, in air it opens the way to many atmospheric applications, including the control of electric discharges⁴⁻⁸ and lightning⁹⁻¹¹, THz to microwave generation¹²⁻¹⁴, air lasing¹⁵, remote sensing¹⁶⁻¹⁸, optical waveguide¹⁹,²⁰, fog clearing²¹ or laser-induced condensation²²⁻²⁵.

These applications are fueled by the ability of laser filaments to form and propagate at distances up to the kilometer-range²⁶,²⁷, forming long, continuous, slightly ionized plasma channels sustained by a dynamic competition between focusing and defocusing effects of different nonlinear orders²⁸⁻³⁰. The resulting intensity clamping, associated to the short pulse duration that drastically limits avalanche ionization, allows the electron density to stay below the critical density and the filaments remain essentially transparent.

In the context of high-voltage discharge and lightning control, the effect of the laser does not only rely on the free charge carriers (electrons) released by ionization. The energy deposited by the ionization creates a shockwave that results in a transient post-pulse depletion of the air density by up to 70%³¹⁻³⁴, increasing the mean free path of the electrons and reducing their attachment rate to positive ions. This effect, known as Paschen's law³⁵, reduces the breakdown threshold and contributes to the discharges guiding and triggering³⁶.

One main limitation of laser-created filaments to trigger lightning over long distance is the limited lifetime of the plasma (in the nanosecond range)³⁷, which at the typical propagation velocity of a guided leader (10⁶ m s⁻<fnref n="1" />)⁶,³⁸, limits the effect of the free electrons to the meter-scale. The air-depleted channel lasts longer, with a decay over several hundred microseconds, and a remaining of a few percent only after one millisecond³⁴. However, the nanosecond to microsecond lifetime of the ionized and air-depleted channels produced by laser filaments is considered as the main limitation for atmospheric-scale control of lightning⁹. Attempts to circumvent this limitation by using pulse trains have shown limited efficiency³⁹,⁴⁰ due to limited total pulse energy (several mJ) and repetition rate (10–20 Hz), hence an average power limited to tens of mW to the Watt range. The use of laser pulses in the UV⁴¹,⁴², far IR⁴³, or of dual laser pulses⁴⁴ or 100 TW laser pulses⁴⁵ has also been studied in laboratory, but these laser sources do not allow generating long-lived extended filaments with the present technology.

However, the recent advent of ultrashort thin-disk technology⁴⁶ allows to simultaneously reach sub-ps pulse duration, Joule-range pulse energy, and a kHz repetition rate⁴⁷. Therefore, a laser with TW-range peak power together with a kW-level average power raises new perspective for long-distance discharge control¹¹,⁴⁷. Beyond a higher ionization due to longer pulses, a kHz repetition rate allows each pulse to propagate in the remaining of the depleted channel left by the previous ones, allowing the build-up of a cumulative effect that is spectacularly

<footnote n="1">¹Laboratoire d'Optique Appliquée, ENSTA Paris, École Polytechnique, CNRS, IP Paris, 828 boulevard des Maréchaux, 91762 Palaiseau, France. ²Groupe de Physique Appliquée, Université de Genève, Ch. de Pinchat 22, 1211 Geneva 4, Switzerland. ³TRUMPF Scientific Lasers GmbH + Co. KG, Feringastr. 10a, 85774 Unterföhring, Germany. ⁴Université Paris-Saclay, CNRS, Laboratoire FAST, 91405 Orsay, France. ⁵André Mysyrowicz Consultants, 6 Rue Gabriel, 78000 Versailles, France. ⁶Institute for Environmental Sciences, Université de Genève, Bd Carl Vogt 66, 1211 Geneva 4, Switzerland. ⁷email: aurelien.houard@polytechnique.edu</footnote>
illustrated in the case of cloud clearing²¹ and results in the establishment of a quasi-steady-state depleted channel prone to a permanent highly-localized Paschen effect⁴⁸,⁴⁹.

While filamentation is well-known, including on atmospheric ranges, for 10 Hz Ti:Sa lasers at 800 nm, thin-disk Yb lasers simultaneously feature a slightly higher wavelength, a longer pulse duration and an associated narrower spectrum, and, more important, a kHz repetition rate allowing thermal coupling between successive pulses. They therefore deserve a full characterization in view of atmospheric applications. Here, we characterize the propagation of such high-peak-power, high-average-power pulses over atmospheric-range distances (140 m), in the multiple filamentation regime, for the fundamental wavelength (1030 nm) as well as the second and third harmonics. We show that the onset distance of filaments can be controlled by both chirp and geometrical focusing, and that filament can propagate up to 140 m distance. We also show that the beam tends to self-organize as a filament ring unless forced into a superfilament by geometrical focusing. Furthermore, the partial frequency-doubling of the pulses perturbs the filamentation of the fundamental wavelength through coupling effect between the two colors. We finally characterize the air density channel in different focusing conditions.

## Experimental setup

We used the thin-disk Yb:YAG laser source described in⁴⁷, which provides pulses at 1030 nm center wavelength with an energy of 500 mJ after compression and a minimum pulse duration of 1 ps, at a repetition rate of 1 kHz. The laser output beam has a beam diameter of about 30 mm (cross section of 700 mm²), with an average $M^2$ value of 2.1. The use of lithium triborate (LBO) crystals enables us to produce 295 mJ (59% conversion efficiency) and 119 mJ (27% conversion efficiency) at respectively the second (515 nm) and the third harmonic (343 nm)⁵⁰. A specially-designed telescope can expand the beam by a factor of 7.14 and focus it at a distance ranging from 55 m to infinity¹⁰. In our measurement, the propagation distance $z$ along the laser axis is defined as the distance between the output of the telescope and the diagnostics. Our experimental setup was installed in a 150 m long hall at the University Paris-Saclay (Orsay, France).

To characterize the presence of filaments along the laser path, beam impact measurements were realized every meter along the laser propagation using photosensitive paper (Ilford Multigrade IV RC 44M). The paper was swept sufficiently fast across the beam path to distinguish the impact of individual shots, while the laser was continuous working at 1 kHz. In addition, to quantify the energy deposition in air by the filaments we characterized the acoustic emission along the laser propagation path. During ionization, the filamentary laser pulses quickly deposit energy in the medium resulting in the creation of a lateral pressure wave that can be measured by a microphone⁵¹⁻⁵³. The evolution of the acoustic signal along the laser path was measured for the same pulse durations used previously by means of a microphone (model G.R.A.S. 46 HB) amplified by an amplifier (model G.R.A.S. 12AK) with a gain of 40 (detection bandwidth 20 Hz–20 kHz). We considered the integrated magnitude of the absolute value of the signal over 0.2 ms as representative of the pressure wave magnitude.

## Results

### Long distance filamentation using a collimated beam

#### Effect of the pulse duration

We first investigated the filamentation of the free-propagating collimated laser beam. In such a configuration, the control of the filamentation distance can be achieved simply by adjusting the laser pulse duration⁵⁴.

The evolution of the number of filaments $N$ as a function of the propagation distance $z$ was measured in the case of 1030 nm 500 mJ pulses at several pulse durations ranging from 1 to 7 ps. The results, displayed in Fig. 1, show that filaments are present over several tens of meters for all pulse durations. As expected from previous observations at 800 nm²⁷,⁵⁴,⁵⁵, and the Marburger formula, longer pulses generate filaments at a longer distance. Indeed, as shown in Table 1, the position of the maximum number of observed filaments matches well the Marburger formula for self-focusing distance $L_c$ applied to the whole beam⁵⁶:

$$L_c = \frac{0.367 z_r}{\sqrt{\left(\frac{P_{in}}{P_{cr}} - 0.852\right)^2 - 0.0219}}, \quad (1)$$

with $P_{in}$ the initial peak power of the laser pulse, $P_{cr} = 4.9$ GW the critical power for self-focusing in air at 1030 nm⁵⁷ and $z_r$ the Rayleigh length ($z_r = 690$ m at 1030 nm).

Furthermore, the onset of filamentation is close to the local application of the Marburger formula on each filamenting cell within the beam, showing that the nucleation of the observed small scale filamentation is governed by the local splitting of the beam profile into self-focusing cells influenced by inhomogeneities originating from the initial beam profile and from the effect of turbulence in the experimental hall.

Finally, let us note that the filamenting region is also longer for longer pulses (Fig. 1), although the distance available for the experiment was limited to 140 m, preventing to measure the full length of the filament region for the pulse duration of 7 ps.

We compared the measured acoustic signal with a simulated signal based on the sum of the pressure waves emitted by all individual filaments detected with the photosensitive paper, each one being located at a different distance from the microphone. This propagation-corrected simulated signal $S$, can be expressed as

$$S = \frac{CN}{2\pi R^2} \int_{0}^{2\pi} \int_{0}^{R} \frac{r dr d\theta}{(D - r \cos\theta)^2 + (r \sin\theta)^2}, \quad (2)$$

where $r$ and $\theta$ are the cylindrical coordinates, $R$ the radius of the filamentation zone, $D$ the distance between the center of the beam and the microphone, $N$ the number of filaments and $C$ is the signal emitted by a single
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Evolution of the number of filaments *N* along the laser propagation axis *z* of a 500 mJ 1030 nm laser pulse, for pulse durations of (a) 1 ps, (b) 3 ps, (c) 5 ps and (d) 7 ps. Examples of photosensitive paper measurements performed at the distances (e) *z* = 19 m, (f) *z* = 25 m, (g) *z* = 29 m for a pulse duration of 1 ps. (h) Comparison between the microphone signal (orange dashed line) and the simulated signal *S* (blue solid line) along the laser propagation axis *z* for a laser pulse of 500 mJ at 1030 nm with a duration of 1 ps.</figcaption></figure>

<table><caption>**Table 1.** Comparison between the distance at which a maximum of filaments is detected and the whole-beam self-focusing distance *L*<sub>c</sub>, for pulse durations of 1 ps, 3 ps, 5 ps and 7 ps.</caption><thead><tr><th></th><th>1</th><th>3</th><th>5</th><th>7</th></tr></thead><tbody><tr><td>Laser pulse duration (ps)</td><td></td><td></td><td></td><td></td></tr><tr><td>Peak power (GW)</td><td>500</td><td>167</td><td>100</td><td>71</td></tr><tr><td>Position of the peak of filament number (m)</td><td>33</td><td>56</td><td>74</td><td>98</td></tr><tr><td>Calculated filamentation distance L<sub>c</sub> (m)</td><td>30</td><td>56</td><td>76</td><td>94</td></tr><tr><td>Length of filamentation region (m)</td><td>50</td><td>70</td><td>80</td><td>&gt; 90</td></tr></tbody></table>

filament. The integral corresponds to the averaging over a homogeneous probability distribution of filaments within the filamentation zone, considering 1000 shots generating independent filament patterns. Since distances between filaments are in the cm-range, well shorter than the acoustic wavelength, we integrate the sound intensity rather than the signal amplitudes. The simulated signal *S* is calculated using the values of *R* and *N* measured from the beam characterization with the photosensitive paper.

The comparison between the experimental acoustic signal and the simulated signal calculated from the photosensitive paper measurements for a laser pulse of 500 mJ at 1030 nm with a duration of 1 ps is shown in Fig. 1h. The clear proportionality between the two signals confirms that the power carried and deposited by each filament is constant over the propagation distance. The same observation is made for all tested pulse durations.

### Multi-color filamentation

We used the frequency doubling LBO crystal to study the filamentation of a dual frequency beam consisting of 285 mJ energy at 515 nm and 215 mJ energy at 1030 nm<sup>50</sup>, by repeating the beam impact measurements. The fundamental laser pulse duration was set to 1 ps. A consequence of the frequency doubling is the modification of the filament diameter *d*<sub>fil</sub>. The proportionality of *d*<sub>fil</sub> with the wavelength is given by<sup>58,59</sup>:

$$d_{fil} = \frac{1.22 \lambda}{\sqrt{8n_0n_2I_{max}}}, \quad (3)$$

with *λ* the beam wavelength, *n*<sub>0</sub> the air refractive index, *n*<sub>2</sub> the nonlinear Kerr refractive index and *I*<sub>max</sub> = 2 10<sup>13</sup> W cm<sup>-2</sup> the intensity of the laser in the filament<sup>57</sup>. The filament diameter being directly proportional to the wavelength, the filaments of the frequency-doubled beam have half the diameter of the ones from the fundamental beam. The two populations of filaments can thus be distinguished on images of the beam impacts on photosensitive paper.

The filament size distributions for the fundamental and dual-color beams indeed display a factor of two in their diameters (Fig. 2). Note that the impacts on the photosensitive paper is typically 5 times larger than the filament diameter reported experimentally and predicted by the theory (120 µm and 62 µm at 1030 nm and 515 nm, respectively)<sup>29,30</sup>. This is explained by the low threshold of our photosensitive paper that is darkened not only by the filament core but also by the photon bath around it, resulting in a burned mark larger than the
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Experimental distribution of the filament widths produced by (a) 500 mJ 1 ps laser pulses at 1030 nm and (b) 500 mJ 1 ps laser pulses at 1030+515 nm measured on photosensitive paper. Insert: Example of impacts of single filaments on photosensitive paper.</figcaption></figure>

filament itself. Remarkably, no filament with a diameter above 0.35 mm (which should be the case for a filament at 1030 nm) is formed while using the dual frequency beam, suggesting that the unconverted fraction of the beam at 1030 nm does not produce filaments despite its 214 mJ energy per pulse (43 $P_{cr}$).

To confirm and understand the latter observation, additional measurements were made. Due to technical issues, the laser was repaired between these two sets of measurements and the beam profile was heavily modified, hence comparison should only be done within the same measurement set. To characterize the filamentation of each wavelength, we measured the number of filaments created by the 1030 nm and by the 515 nm beam by identifying them from impact beam measurements. Filaments up to a diameter of 0.4 mm are identified as 515 nm filaments, while any filament above this diameter is deemed a 1030 nm filament. Four configurations were used: the fundamental beam, the dual-frequency beam, and the individual fundamental and second harmonic beams separated from the dual-frequency one by means of a dichroic beam splitter. The results are presented in Fig. 3. The depletion of the fundamental by the frequency doubling results in a much weaker and later filamentation (compare panels a and c). This difference is mainly explained by the fact that the non-linear depletion of the pump beam tends to self-regularize the beam, i.e., smooth the gradients within the beam profile due to higher conversion efficiency (hence higher losses) in the most intense regions. Such self-regularization is expected to hinder filamentation by reducing the seeding of spatial modulational instability. Modeling this effect based on the depleted pump formalism$^{60}$, in our conditions featuring 57% conversion efficiency into the SHG$^{50}$, shows that the gradients are reduced by a factor of $\sqrt{2}$.

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Evolution of the number of filaments *N* at 1030 nm (red solid line) and at 515 nm (green blue dashed line) along the laser propagation, created by (a) the fundamental 1030 nm only, (b) 1030+515 nm pulses, (c) the isolated fundamental part of 1030 + 515 nm pulses and (d) the isolated second harmonic part of the 1030 + 515 nm pulses. The dotted line corresponds to the total number of filaments at 1030 nm and 515 nm.</figcaption></figure>
On the other hand, frequency-doubling the beam (compare panels a and d) results in more than twice as many filaments, with an earlier onset. This increase is in line with the fact that the SHG conversion efficiency is beyond 50%, which, together with the $\lambda^2$ dependence of the critical power<sup>29</sup> results in more than doubling the normalized power $P/P_{cr}$. More specifically, the critical power at 1030 nm amounts to 4.9 GW, so that 500 mJ in 1 ps corresponds to 100$P_{cr}$. At a rate of 1 filament per $5P_{cr}$<sup>61,62</sup>, we can therefore expect up to 20 filaments, about twice the number observed. The 12 filaments observed correspond to $8P_{cr}$/filament. In the case of the frequency-doubled beam at 515 nm, the critical power amounts to 1.24 GW, so that the 285 mJ, corresponding to 285 GW peak power if we assume a 1-ps pulse duration, corresponds to 230$P_{cr}$. The 26 filaments observed correspond to approximately 9$P_{cr}$/ filament, again slightly above the ratio expected from measurements for shorter pulses. The beam diameter of 30 mm has a cross section of 700 mm<sup>2</sup>, corresponding to more than 26 mm<sup>2</sup>/filament. This value lies far from the geometrical saturation, which amounts to only a few mm<sup>2</sup>/filament<sup>63</sup>. Therefore, spatial congestion is not the factor limiting the number of filaments.

Finally comparing the dual-color beam (panel b) with the individual colors (panels c and d) evidences a nonlinear interaction between the two frequencies during the propagation. This interaction (1) reduces the total number of filaments (compare panel b with panels c and d), (2) increases the total filamenting length, (3) restores / anticipates the filamentation of the remaining fundamental.

One other important observation is that the isolated fundamental part and second harmonic part of the dual frequency beam show a different filamentation process when they are isolated from each other, as seen in Fig. 3c,d. When isolated, the fundamental part of the beam only starts filamentation after a propagation of more than 135 m (as compared to 40 m in the dual frequency beam). This can be explained by the regularization effect and the reduction of the available power induced by the SHG crystal on the fundamental beam. On the other hand, the isolated second harmonic part of the beam starts the filamentation process after the same propagation distance of 13 m but the maximum number of filaments is higher and their spread is limited to a smaller distance compared to what is measured when using the dual frequency beam. These two observations suggest that the two parts of the beam interact during the propagation and improve the filamentation of the fundamental part of the beam while hindering the filamentation of the second harmonic part. It is important to observe that, due to a phase velocity mismatch in air of 46.5 fs m<sup>-1</sup> between these two spectral components, they do not temporally overlap anymore once the filament process begins after about 20 m of propagation, so the optical interaction of the two parts of the beam only occurs over this distance. But they can also interact through the ionization or the molecular orientation left by the first IR pulse, or through cumulative thermal effects that persists over very long timescale.

Another aspect of the beam affected by the nonlinear propagation is its spectrum. The high beam intensity preventing direct measurement of its spectrum, we characterized the laser spectrum at different distances by measuring with a spectrometer the spectrum of the light scattered on a ceramic plate placed on the beam path, using the triple frequency-converted pulse at 1 ps<sup>50</sup>.

The results, presented in Fig. 4a-c show that the spectral broadening of the mixed beam is much more significant for the second harmonic than for the fundamental and the third harmonic. In the second harmonic component, the spectral width rises from 1.0 to 7.6 nm, which would allow the Fourier-limited pulse durations to drop from 389 to 51 fs. Note that the broadening we observe occurs during the propagation, well after the SHG crystal. The broadening during the propagation before and in the crystal is already taken as an initial condition.

The redshift of the 515 nm beam can be explained by the retarded Kerr effect in air<sup>64</sup> also referred as Stimulated Raman Scattering<sup>42,65</sup>. Based on the model presented in<sup>27</sup>, the red shift $\Delta\lambda$ is calculated along the laser propagation using the equation:

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Evolution of the spectrum of the triple frequency-converted pulse along its propagation. (a-c) being the spectral components around 343, 515 and 1030 nm, respectively. (d) Comparison between the measured and calculated redshift along the laser propagation, for the pulse spectral component at 515 nm.</figcaption></figure>
$$ \Delta\lambda(z) = \lambda_0 \frac{n_2 I_0}{\tau \sqrt{2c}} \frac{\Delta t_0}{\Delta t} (z - z_0), \quad (4) $$

with $\lambda_0$ the initial wavelength, $n_2$ the instantaneous nonlinear index of air, $I_0$ the laser intensity, $\tau$ the characteristic time of the partial alignment of the molecules (1.5 ps), $c$ the speed of light, $\Delta t_0$ the minimum pulse duration, $\Delta t$ the pulse duration, $z$ the distance of propagation and $z_0$ the onset position of the filament. The calculated values, presented in Fig. 4d, show a good agreement between the theory and the measurement.

The same measurement was repeated for pulse durations of 1 ps, 3 ps and 5 ps. For all tested pulse durations, the spectral broadening of the mixed beam is more important for the second harmonic than for the fundamental and the third harmonic. Furthermore, longer pulse durations delay the onset of the spectral broadening and decrease its importance. The distances where the spectral broadening starts are consistent with the distances of onset of filaments for all pulse durations. Moreover, the broadening only takes place over a limited distance corresponding to the filamentation range. Beyond, no further evolution of the spectrum is observed.

## Long distance filamentation using a focused beam

In a second time, we studied the filamentation induced after our large-aperture focusing telescope (see description in$^{10}$). This configuration is more relevant for applications such as the laser lightning rod$^{11}$ because it allows a control over the collapse distance as well as stronger energy deposition through superfilamentation regime$^{66}$.

### Filamentation for different focusing distances

By changing the focal length of the telescope, we studied the filamentation process when the beam is focused at 55 m (Numerical Aperture = 0.0023) or 115 m (N.A. = 0.0011). As for the collimated beam, beam impact on photosensitive paper were used to characterize the presence of filaments along the laser path. However, in this configuration individual filaments cannot be distinguished because they are squeezed by the overall focusing of the beam. Therefore, the total surface occupied by filaments is used as a proxy for the number of filaments. This measurement was repeated for the two focusing distances using the fundamental and the dual-color beams, for pulse durations varying from 1 to 7 ps.

The results, presented in Fig. 5, show that the use of the frequency doubled beam leads to a higher surface occupied by filaments compared to the fundamental beam, albeit the difference is small. Furthermore, only small differences were observed between measurements at different pulse durations. This is an important observation because it allows us to use long pulse duration to reduce the laser intensity and prevent damage on the optics of our setup while maintaining a high number of filaments. On the contrary, an important difference is observed when comparing the two focusing distances. The longer focusing distance of 115 m results in a lower maximum of the surface occupied by filaments but in a longer filamentation length (defined as the length over which at least one filament is detected), of more than 30 m compared to 12 m for the focusing distance of 55 m. We also note the apparition of a ring like organization of the filaments in the case where the beam is focused at 115 m, while a dense pack of filaments is observed with the focusing distance of 55 m as illustrated in Fig. 5c,d. The latter can be interpreted as a superfilamentation regime$^{66}$.

In addition, we characterized the low-density air channel resulting from the energy deposition by the filamentation process$^{34,67}$. By measuring the phase shift induced by this density depletion on a transverse probe beam (see$^{48}$ for a description of the setup), we characterized the low-density channel present at the focal point for the two focusing distances at $z = 55$ m and $z = 115$ m, using the fundamental and frequency doubled beam. The measurement was performed after the shockwave has developed, 1 µs after the laser pulse, corresponding to the peak of the air density depletion$^{34}$.

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Spatial evolution of the surface occupied by filaments created by the 500 mJ 1 ps fundamental beam (red solid line) or the frequency doubled beam (green dashed line), focused at (a) $z = 55$ m or (b) $z = 115$ m. Examples of photosensitive paper measurements performed at the respective focal point (c) $z = 55$ m and (d) $z = 115$ m.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 6.** (a) Phase shift $\Delta\phi$ induced by the frequency doubled 500 mJ, 1 ps pulse, focused at $z = 55$ m (red solid line) and $z = 115$ m (blue solid line), measured at the geometrical focus of the beam 1 µs after the last laser pulse. (b) Phase shift $\Delta\phi$ induced by the 500 mJ, 1 ps frequency doubled beam (blue solid line) and the 500 mJ, 1 ps fundamental beam (red dashed line), focused at $z = 55$ m, measured at the focal distance 1 µs after the last laser pulse.</figcaption></figure>

As seen in Fig. 6a, we observe an asymmetric phase shift, characteristic of filamentation at high repetition rate⁴⁸. This shape is due to the accumulation of successive low-density channels going progressively upward due to Archimedes thrust. But the transverse size of this phase shift (10–15 mm) is significantly larger than the one produced by a single filament (100 µm FWHM) or by a kHz laser at 800 nm with an energy of 2 mJ (3 mm FWHM). We also note that the measured phase shift is more important when the beam is focused at $z = 55$ m compared to $z = 110$ m, as expected since the transversely integrated density of deposited energy increases with the laser numerical aperture. This observation is consistent with the previously measured higher maximum of surface occupied by the filament, when the beam is focused at 55 m. However, while the previous measurements of the occupied surface have shown only minor differences between the fundamental and frequency doubled beam, the measured phase shifts $\Delta\phi$ show an important increase when using the frequency doubled beam, as seen in Fig. 6b. This suggests that while the number of filaments is not higher when using the frequency doubled beam, each filament deposits more energy that result in a greater depletion of the air density.

## Conclusion

We have characterized the filamentation of a 500 mJ, 1 ps, 1 kHz laser system in different focusing conditions. Using a collimated beam, the presence of filaments over 50 m has been observed and the onset of the filamentation was shifted from 18 to 60 m by changing the pulse duration from 1 to 7 ps. Using a large-aperture focusing telescope, we characterized the filaments created by the beam focused at 55 m and 115 m. Filaments created at 55 m are more locally dense than the ones created at 115 m, but they spread over a distance of 12 m compared to the 30 m covered by the ones created at 115 m.

We characterized the low-density air channel resulting from the filamentation and have shown that the one resulting from a tighter focus is more important. The use of frequency doubling crystal also resulted in a more important reduction of the air density.

These measurements were crucial for the main experiment of the Laser Lightning Rod project¹⁰,¹¹ where this laser system was used to guide ascending lightning leader starting at the top of a tower. Our ability to create filaments at a distance of 115 m was a requirement as this is the tower height, and the measured presence of filaments over 30 m would be an important parameter to scale the guiding effect.

## Data availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Received: 6 July 2023; Accepted: 22 October 2023

Published online: 29 October 2023

## References

1. Maiman, T. H. Stimulated optical radiation in ruby. *Nature* https://doi.org/10.1038/187493a0 (1960).

2. Strickland, D. & Mourou, G. Compression of amplified chirped optical pulses. *Opt. Commun.* **56**, 219. https://doi.org/10.1016/0030-4018(85)90120-8 (1985).

3. Braun, A. et al. Self-channeling of high-peak-power femtosecond laser pulses in air. *Opt. Lett.* **20**, 73. https://doi.org/10.1364/OL.20.000073 (1995).

4. Zhao, X. M., Diels, J.-C., Wang, C. Y. & Elizondo, J. M. Femtosecond ultraviolet laser pulse induced lightning discharges in gases. *IEEE J. Quantum Electron.* **31**, 599. https://doi.org/10.1109/3.364418 (1995).

5. Comtois, D. et al. Triggering and guiding leader discharges using a plasma channel created by an ultrashort laser pulse. *Appl. Phys. Lett.* **76**, 819. https://doi.org/10.1063/1.125595 (2000).

6. Rodriguez, M. et al. Triggering and guiding megavolt discharges by use of laser-induced ionized filaments. *Opt. Lett.* **27**, 772. https://doi.org/10.1364/OL.27.000772 (2002).

7. Kosareva, O. G. et al. Remote triggering of air-gap discharge by a femtosecond laser filament and postfilament at distances up to 80 m. *Appl. Phys. Lett.* **119**, 041103. https://doi.org/10.1063/5.0057544 (2021).

8. Pei, Z. et al. The contribution of femtosecond laser filaments to positive and negative breakdown discharge in a long air gap. *Phys. Plasmas* **30**, 043511. https://doi.org/10.1063/5.0138646 (2023).
9. Kasparian, J. et al. Electric events synchronized with laser filaments in thunderclouds. *Opt. Express* **16**, 5757. https://doi.org/10.1364/OE.16.005757 (2008).

10. Produit, T. et al. The laser lightning rod project. *Eur. Phys. J. Appl. Phys.* **93**, 10504. https://doi.org/10.1051/epjap/2020200243 (2021).

11. Houard, A. et al. Laser-guided lightning. *Nat. Photon.* **17**, 231. https://doi.org/10.1038/s41566-022-01139-z (2023).

12. D'Amico, C. et al. Conical forward THz emission from femtosecond-laser-beam filamentation in air. *Phys. Rev. Lett.* **98**, 235002. https://doi.org/10.1103/PhysRevLett.98.235002 (2007).

13. Englesbe, A. et al. Ultrabroadband microwave radiation from near- and mid-infrared laser-produced plasmas in air. *Phys. Rev. A* **104**, 013107. https://doi.org/10.1103/PhysRevA.104.013107 (2021).

14. Mitrofanov, A. V. et al. Broadband ultrawide-angle laser-plasma microwave antennas. *Phys. Rev. A* **105**, 053503. https://doi.org/10.1103/PhysRevA.105.053503 (2022).

15. Dogariu, A., Michael, J. B., Scully, M. O. & Miles, R. B. High-gain backward lasing in air. *Science* **331**, 442. https://doi.org/10.1126/science.1199492 (2011).

16. Kasparian, J. et al. White-light filaments for atmospheric analysis. *Science* **301**, 61. https://doi.org/10.1126/science.1085020 (2003).

17. Galvez, M. C. et al. Three-wavelength backscatter measurement of clouds and aerosols using a white light lidar system. *Jpn. J. Appl. Phys.* **41**, L284. https://doi.org/10.1143/JJAP.41.L284 (2002).

18. Stelmaszczyk, K. et al. Long-distance remote laser-induced breakdown spectroscopy using filamentation in air. *Appl. Phys. Lett.* **85**, 3977. https://doi.org/10.1063/1.1812843 (2004).

19. Jhajj, N., Rosenthal, E. W., Birnbaum, R., Wahlstrand, J. K. & Milchberg, H. M. Demonstration of long-lived high-power optical waveguides in air. *Phys. Rev. X* **4**, 011027. https://doi.org/10.1103/PhysRevX.4.011027 (2014).

20. Fu, S., Mahieu, B., Mysyrowicz, A. & Houard, A. Femtosecond filamentation of optical vortices for the generation of optical air waveguides. *Opt. Lett.* **47**, 5228. https://doi.org/10.1364/OL.47.2143 (2022).

21. Schimmel, G., Produit, T., Mongin, D., Kasparian, J. & Wolf, J.-P. Free space laser telecommunication through fog. *Optica* **5**, 1338. https://doi.org/10.1364/OPTICA.5.001338 (2018).

22. Henin, S. et al. Field measurements suggest the mechanism of laser-assisted water condensation. *Nat. Commun.* **2**, 456. https://doi.org/10.1038/ncomms1462 (2011).

23. Ju, J. et al. Laser-filamentation-induced condensation and snow formation in a cloud chamber. *Opt. Lett.* **37**, 1214. https://doi.org/10.1364/OL.37.001214 (2012).

24. Yoshihara, K. et al. UV-light-induced water condensation in air and the role of hydrogen peroxide. *Bull. Chem. Soc. Jpn.* **87**, 593. https://doi.org/10.1246/bcsj.20130322 (2014).

25. Shumakova, V. et al. Laser induced aerosol formation mediated by resonant excitation of volatile organic compounds. *Optica* **8**, 1256. https://doi.org/10.1364/OPTICA.43.4659 (2021).

26. Rodriguez, M. et al. Kilometer-range nonlinear propagation of femtosecond laser pulses. *Phys. Rev. E* **69**, 036607. https://doi.org/10.1103/PhysRevE.69.036607 (2004).

27. Durand, M. et al. Kilometer range filamentation. *Opt. Express* **21**, 26836. https://doi.org/10.1364/OE.21.026836 (2013).

28. Chin, S. L. et al. The propagation of powerful femtosecond laser pulses in opticalmedia: Physics, applications, and new challenges. *Can. J. Phys.* **83**, 863. https://doi.org/10.1139/p05-048 (2005).

29. Couairon, A. & Mysyrowicz, A. Femtosecond filamentation in transparent media. *Phys. Rep.* **441**, 47. https://doi.org/10.1016/j.physrep.2006.12.005 (2007).

30. Bergé, L., Skupin, S., Nuter, R., Kasparian, J. & Wolf, J.-P. Ultrashort filaments of light in weakly ionized, optically transparent media. *Rep. Prog. Phys.* **70**, 1633. https://doi.org/10.1088/0034-4885/70/10/R03 (2007).

31. Vidal, F. et al. Modeling the triggering of streamers in air by ultrashort laser pulses. *IEEE Trans. Plasma Sci.* **28**, 418. https://doi.org/10.1109/27.848101 (2000).

32. Cheng, Y.-H., Wahlstrand, J. K., Jhajj, N. & Milchberg, H. M. The effect of long timescale gas dynamics on femtosecond filamentation. *Opt. Express* **21**, 4740. https://doi.org/10.1364/OE.21.004740 (2013).

33. Lahav, O. et al. Long-lived waveguides and sound-wave generation by laser filamentation. *Phys. Rev. A* **90**, 021801(R). https://doi.org/10.1103/PhysRevA.90.021801 (2014).

34. Point, G., Milián, C., Couairon, A., Mysyrowicz, A. & Houard, A. Generation of long-lived underdense channels using femtosecond filamentation in air. *J. Phys. B: At. Mol. Opt. Phys.* **48**, 094009. https://doi.org/10.1088/0953-4075/48/9/094009 (2015).

35. Tirumala, R. & Go, D. B. An analytical formulation for the modified Paschen's curve. *Appl. Phys. Lett.* **97**, 151502. https://doi.org/10.1063/1.3497231 (2010).

36. Tzortzakis, S. et al. Femtosecond laser-guided electric discharge in air. *Phys. Rev. E* **64**, 057401. https://doi.org/10.1103/PhysRevE.64.057401 (2001).

37. Tzortzakis, S., Prade, B., Franco, F. & Mysyrowicz, A. Time-evolution of the plasma channel at the trail of a self-guided IR femtosecond laser pulse in air. *Opt. Commun.* **181**, 123. https://doi.org/10.1016/S0030-4018(00)00734-3 (2000).

38. La Fontaine, B. et al. Guiding large-scale spark discharges with ultrashort pulse laser filaments. *J. Appl. Phys.* **88**, 610. https://doi.org/10.1063/1.373710 (2000).

39. Lu, X. et al. Quasi-steady-state air plasma channel produced by a femtosecond laser pulse sequence. *Sci. Rep.* **5**, 15515. https://doi.org/10.1038/srep15515 (2015).

40. Schubert, E. et al. Optimal laser-pulse energy partitioning for air ionization. *Phys. Rev. A* **94**, 033824. https://doi.org/10.1103/PhysRevA.94.033824 (2016).

41. Rambo, P., Schwarz, J. & Diels, J.-C. High-voltage electrical discharges induced by an ultrashort-pulse UV laser system. *J. Opt. A: Pure Appl. Opt.* **3**, 146. https://doi.org/10.1088/1464-4258/3/2/309 (2001).

42. Zvorykin, V. et al. Range of multiple filamentation of a terawatt-power large-aperture krT laser beam in atmospheric air. *J. Opt. Soc. Am. B* **36**, G25-G32. https://doi.org/10.1364/JOSAB.36.000G25 (2019).

43. Tochitsky, S. et al. Megafilament in air formed by self-guided terawatt long-wavelength infrared laser. *Nat. Photon.* **13**, 41. https://doi.org/10.1038/s41566-018-0315-0 (2019).

44. Scheller, M. et al. Externally refuelled optical filaments. *Nat. Photon.* **8**, 297. https://doi.org/10.1038/nphoton.2014.47 (2014).

45. Henin, S. E. A. Saturation of the filament density of ultrashort intense laser pulses in air. *Appl. Phys. B* **100**, 77-84. https://doi.org/10.1007/s00340-010-3941-x (2010).

46. Saraceno, C. J., Sutter, D., Metzger, T. & Abdou Ahmed, M. The amazing progress of high-power ultrafast thin-disk lasers. *J. Euro. Opt. Soc. Rapid Publ.* **15**, 15. https://doi.org/10.1186/s41476-019-0108-1 (2019).

47. Herkommer, C. et al. Ultrafast thin-disk multipass amplifier with 720 mJ operating at kilohertz repetition rate for applications in atmospheric research. *Opt. Express* **28**, 30164. https://doi.org/10.1364/OE.40.4185 (2020).

48. Walch, P. et al. Cumulative air density depletion during high repetition rate filamentation of femtosecond laser pulses: Application to electric discharge triggering. *Appl. Phys. Lett.* **119**, 264101. https://doi.org/10.1063/5.0077635 (2021).

49. Goffin, A., Tartaro, A. & Milchberg, H. M. Quasi-steady-state air waveguide. *Optica* **10**, 505. https://doi.org/10.1364/OPTICA.48.292 (2023).

50. Andral, U. et al. Second and third harmonic generation from simultaneous high peak- and high average-power thin disk laser. *Appl. Phys. B* **128**, 177. https://doi.org/10.1007/s00340-022-07887-8 (2022).

51. Yu, J. et al. Sonographic probing of laser filaments in air. *Appl. Opt.* **42**, 7117. https://doi.org/10.1364/AO.42.007117 (2003).
52. Point, G., Thouin, E., Mysyrowicz, A. & Houard, A. Energy deposition from focused terawatt laser pulses in air undergoing multifilamentation. *Opt. Express* **24**, 6271. https://doi.org/10.1364/OE.24.006271 (2016).

53. Larkin, I. et al. Full path single-shot imaging of femtosecond pulse collapse in air turbulence. *Opt. Lett.* **45**, 2518. https://doi.org/10.1364/OL.389495 (2020).

54. Wille, H. et al. Teramobile: A mobile femtosecond-terawatt laser and detection system. *Eur. Phys. J. Appl. Phys.* **20**, 183. https://doi.org/10.1051/epjap:2002090 (2002).

55. Méchain, G. et al. Range of plasma filaments created in air by a multi-terawatt femtosecond laser. *Opt. Commun.* **247**, 171. https://doi.org/10.1016/j.optcom.2004.11.052 (2005).

56. Marburger, J. H. Self-focusing: Theory. *Prog. Quantum Electron.* **4**, 35. https://doi.org/10.1016/0079-6727(75)90003-8 (1975).

57. Houard, A. et al. Study of filamentation with a high power high repetition rate ps laser at 1.03 µm. *Opt. Express* **24**, 7437. https://doi.org/10.1364/OE.24.007437 (2016).

58. Brewer, R. G., Lifsitz, J. R., Garmire, E., Chiao, R. Y. & Townes, C. H. Small-scale trapped filaments in intense laser beams. *Phys. Rev.* **166**, 326. https://doi.org/10.1103/PhysRev.166.326 (1968).

59. Chiao, R. et al. A new class of trapped light filaments. *IEEE J. Quantum Electron.* **2**, 467. https://doi.org/10.1109/JQE.1966.1074117 (1966).

60. Boyd, R. W. *Non-linear optics* 3rd edn. (Elsevier, 2008).

61. Campillo, A. J., Shapiro, S. L. & Suydam, B. R. Periodic breakup of optical beams due to self-focusing. *Appl. Phys. Lett.* **23**, 628. https://doi.org/10.1063/1.1654772 (1973).

62. Méjean, G. et al. Multifilamentation transmission through fog. *Phys. Rev. E* **72**, 026611. https://doi.org/10.1103/PhysRevE.72.026611 (2005).

63. Petrarca, M. et al. White-light femtosecond lidar at 100 TW power level. *Appl. Phys. B* **114**, 319. https://doi.org/10.1007/s00340-013-5741-6 (2014).

64. Ripoche, J.-F. et al. Determination of the time dependence of n2 in air. *Opt. Commun.* **135**, 310–314. https://doi.org/10.1016/S0030-4018(96)00675-X (1997).

65. Smetanin, I., Levchenko, A., Shutov, A., Ustinovskii, N. & Zvorykin, V. Role of coherent resonant nonlinear processes in the ultrashort KrF laser pulse propagation and filamentation in air. *Nucl. Instrum. Methods Phys. Res., Sect. B* **369**, 87–91. https://doi.org/10.1016/j.nimb.2015.10.032 (2016).

66. Point, G. et al. Superfilamentation in air. *Phys. Rev. Lett.* **112**, 223902. https://doi.org/10.1103/PhysRevLett.112.223902 (2014).

67. Elias, P.-Q. et al. Improving supersonic flights with femtosecond laser filamentation. *Sci. Adv.* **4**, eaau5239. https://doi.org/10.1126/sciadv.aau5239 (2018).

## Acknowledgements

We acknowledge the strong support from Xavier Maître, Moana Pittman, IJCLab and University Paris-Saclay. We acknowledge Rolland Pellet, Laurent Reuenberger, Christophe Vassaux, and their colleagues at the mechanical workshops of the University of Geneva for the construction of the telescope.

## Author contributions

A.H., J.P.W., A.M., T.M. obtained the funding and supervised the project. P.W., T.P., B.M., U.A., Y.B.A., M.M., J.K., J.P.W., B.M., P.W., A.H. designed and installed the experiment. Y.B. provided technical assistance in the experimental hall. C.H., R.J., R.B., S.K., T.M. developed the laser system. P.W., B.M., T.P., V.M., U.A., J.K., J.P.W., Y.B.A., M.L., L.B., A.H., J.P.W. conducted the experiments on site. P.W., J.K., B.M., V.M., T.P. analyzed the data. P.W., J.K., A.H. prepared the original draft and all authors critically reviewed and approved the manuscript.

## Funding

European Union Horizon 2020 Research and innovation programme FET-OPEN grant (737033-LLR). U. Andral and J.-P. Wolf acknowledge financial support from the Swiss National Science Foundation (Grant No. 200021-178926). L. Bizet, P. Walch, and Y.-B. André acknowledge financial support from the French Direction Générale de l'Armement.

## Competing interests

The authors declare no competing interests.

## Additional information

Correspondence and requests for materials should be addressed to A.H.

Reprints and permissions information is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_9.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2023