# Spatio-temporal structure of a petawatt femtosecond laser beam

Antoine Jeandet, Antonin Borot, Kei Nakamura, Spencer W. Jolly, Anthony

J. Gonsalves, Csaba Toth, Hann-Shin Mao, Wim P. Leemans, Fabien Quéré

► To cite this version:

Antoine Jeandet, Antonin Borot, Kei Nakamura, Spencer W. Jolly, Anthony J. Gonsalves, et al. Spatio-temporal structure of a petawatt femtosecond laser beam. Journal of Physics: Photonics, 2019, 1 (3), pp.035001. 10.1088/2515-7647/ab250d. cea-02277420

HAL Id: cea-02277420

https://cea.hal.science/cea-02277420

Submitted on 3 Sep 2019

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 Attribution 4.0 International License
PAPER • OPEN ACCESS

Spatio-temporal structure of a petawatt femtosecond laser beam

To cite this article: Antoine Jeandet et al 2019 J. Phys. Photonics 1 035001

View the [article online](#) for updates and enhancements.
<figure><img src="image_2.png" /></figure>

CrossMark

OPEN ACCESS

RECEIVED
27 February 2019

REVISED
19 April 2019

ACCEPTED FOR PUBLICATION
28 May 2019

PUBLISHED
13 June 2019

Original content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence.

Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.

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

PAPER

Spatio-temporal structure of a petawatt femtosecond laser beam

Antoine Jeandet¹,², Antonin Borot¹, Kei Nakamura³, Spencer W Jolly¹, Anthony J Gonsalves³, Csaba Tóth³, Hann-Shin Mao³, Wim P Leemans³,⁴ and Fabien Quéré¹

¹ LIDYL, CEA, CNRS, Université Paris-Saclay, CEA Saclay, F-91191 Gif-sur-Yvette, France

² Amplitude Laser Group, Business Unit Science, 2/4 rue du Bois Chaland, F-91090 Lisses, France

³ Lawrence Berkeley National Laboratory, Berkeley, CA 94720, United States of America

⁴ Currently with Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany.

E-mail: knakamura@lbl.gov and fabien.quere@cea.fr

Keywords: ultrafast lasers, ultrafast metrology, high power lasers, ultrashort pulses

Supplementary material for this article is available online

Abstract

The development of optical metrology suited to ultrafast lasers has played a key role in the progress of these light sources in the last few decades. Measurement techniques providing the complete *E*-field of ultrashort laser beams in both time and space are now being developed. Yet, they had so far not been applied to the most powerful ultrashort lasers, which reach the PetaWatt range by pushing the chirped pulse amplification (CPA) scheme to its present technical limits. This situation left doubts on their actual performance, and in particular on the peak intensity they can reach at focus. In this article we present the first complete spatio-temporal characterization of a PetaWatt femtosecond laser operating at full intensity, the BELLA laser, using two recently-developed independent measurement techniques. Our results demonstrate that, with adequate optimization, the CPA technique is still suitable at these extreme scales, i.e. it is not inherently limited by spatio-temporal couplings. We also show how these measurements provide unprecedented insight into the physics and operation regime of such laser systems.

Introduction

The technology of ultrafast laser sources nowadays makes it possible to generate laser pulses of femtosecond duration with peak powers of up to several PetaWatts (PW) [1–3]. The development of such laser sources, now available or under construction in several laboratories worldwide, is motivated by two main prospects. The first is demonstrating compact particle accelerators for scientific or societal applications, with particle energies up to several GeV [4–8]. The second is exploring the physics of ultra-relativistic laser-matter interactions, and more particularly accessing regimes where highly nonlinear quantum electro-dynamical effects come into play, in order to perform new tests of this fundamental theory [9].

Due to their ultrashort durations, such lasers necessarily have broad spectra, typically covering several tens of nanometers in wavelength. As in any optical system using broadband light, chromaticity—i.e. the spectral dependence of the spatial properties of a beam, or of the spatial response of an optical system—can become a major impediment to optimize their performance. PetaWatt ultrashort lasers are particularly exposed to this issue; the key enabling technology for these systems, chirped pulse amplification (CPA), which is pushed to its present technical limits in these systems, relies on highly chromatic optical elements such as gratings or prisms, used to tailor the pulses temporal properties for the amplification process [10]. The large chromatic effects individually introduced by each of these optical elements should ideally perfectly compensate each other at the system output. But any imperfection or misalignment of these optics results in residual chromaticity, as so does the use of simple chromatic elements such as singlet lenses [11].

In standard optical assemblies, such as photographic lenses, chromaticity degrades the system performance by affecting its spatial response. The consequences are actually much more severe for ultrashort lasers, where all frequencies not only need to be focused on the same focal spot, but also with the appropriate relative phases, i.e.
those leading to the minimal pulse duration. In this case, undesired chromatic effects not only degrade the concentration of energy in space by inducing imperfect overlap of the different frequencies at focus, they can also increase the pulse duration by reducing the spectral width and degrading the spectral phase locally. The combination of these spatial and temporal effects, known in ultrafast optics as spatio-temporal or spatio-spectral couplings, can lead to very significant reductions in peak intensity even for relatively modest chromatic aberrations—especially as the beam size increases and the pulse duration decreases.

The first step in addressing this issue is the development of adequate metrology tools, which can measure the complete *E*-field of the laser beam in space-time, *E(x, y, t)* or space-frequency, *Ẽ(x, y, ω)*. Several suitable techniques have been developed in the last 15 years [12–17], and are now becoming more and more broadly used, especially on ultrashort lasers of moderate power. Only recently, two different and complementary techniques, TERMITES [18] and INSIGHT [19], have been successfully demonstrated for the characterization of laser systems of up to 100 TW peak power. Yet, none of the existing PW-class lasers have ever been characterized in space-time, thus leaving serious doubts on the actual performance of these cutting-edge sources.

In this article, we address this key issue by presenting the first complete spatio-temporal characterization of a 1.1 PW–40 fs system, using both TERMITES and INSIGHT, carried out on the BELLA PW laser at LBNL, operated at full power and intensity on target, thus revealing for the first time the type of spatio-temporal couplings that can affect such cutting-edge laser systems. This advanced characterization provides unprecedented insight into the physics and operation regime of such laser systems, and demonstrates that the CPA technique is still suitable at such extreme pulse powers. We first detail the implementation of these two techniques and summarize the validation tests that have been performed, and then present the main chromatic effects observed on the beam. We quantify the impact of chromaticity on the laser performance, and explain how the measurements provide insight into the laser system operating conditions and the physics of the beam amplification.

## Experimental set-up and measurement procedure

The BELLA laser, acquired from Thales and operated by LBNL for the last 5 years, is based on a double-CPA architecture: broadband pulses originating from an oscillator are amplified in a succession of six Ti:Sapphire amplifiers, and a cross-polarized wave contrast enhancement system is installed in between these two CPA systems. The exact system design is detailed in [1].

After temporal compression in the grating compressor, the beam of 20 cm diameter propagates under vacuum and is focused into the main experimental chamber by an off-axis parabola with a focal length $f = 13.5$ m, producing a focal spot with $w_0 \approx 53 \ \mu m$ (half width at 1/e² of the field intensity). This chamber is generally used for laser-wakefield particle acceleration experiments, with a capillary discharge plasma as a target (T in figure 1). Here this target was removed from the beam path, letting the laser diverge away from the focus, until it gets large enough to avoid damaging subsequent optics.

An all-reflective chromatic-aberration-free telescope is then used to image the best focus, located in the target chamber, to an output diagnostic area, with a magnification factor of 1. This telescope has been extensively used for the diagnostics of all preceding experiments on BELLA (see [1]) and was proven to introduce negligibly small aberrations. The set of relay optics in this telescope includes three uncoated wedges (two in vacuum, one in air) to attenuate laser pulse energy by a factor of $2.7 \times 10^4$. When BELLA is operated above 5 J pulse energy, a thin (5 $\mu$m) self-standing gold-coated pellicle with 6 inch inner diameter (G in figure 1) is inserted into the beam for a further attenuation by a factor 40. This unique attenuation system is what allows the beam to be safely operated at full power inside the vacuum chamber, and therefore also at full intensity at focus, while avoiding optical damage or nonlinear effects in the subsequent diagnostic line (in particular in the exit vacuum window W in figure 1).

The diagnostic area located after this attenuation system incorporates the spatial-only and temporal-only measurement devices typically used for such lasers, as indicated in figure 1. In the present work, we have implemented two additional measurement devices, based on the TERMITES and INSIGHT techniques, in order to determine the spatio-temporal field of the BELLA laser beam.

TERMITES is typically implemented on collimated beams. We therefore recollimated the converging BELLA image beam to a diameter of about 1.8 cm, by inserting a convex spherical mirror prior to the focus image plane (CM in figure 1). TERMITES then uses a Michelson interferometer with a flat mirror in one arm, and a convex spherical mirror in the other. In this second arm, the central part of the beam is converted into a diverging beam, which eventually fully overlaps with the beam to be characterized and is used as a reference, as in point diffraction interferometry. The spatio-spectral comparison of these two beams is achieved by measuring, on each pixel of the detector, the signal resulting from their superposition as a function of delay. Their cross-spectral density function is then obtained by a Fourier transformation of this first order cross-correlation
<figure><img src="image_4.png" /></figure>

function. The design details of TERMITES are described in [18] and its associated supplementary material. In order to fit the device to the experimental conditions of this particular measurement, we chose the reference size to be half the diameter of the incident beam.

INSIGHT also relies on a Michelson interferometer, now with a flat mirror in each arm. This is used to measure the spatially-resolved first order autocorrelation function of the beam at each point of the detector (standard Fourier-transform spectroscopy), from which one can deduce the frequency-resolved spatial intensity profile of the beam. As detailed in [19], any pointing fluctuations must be compensated for in order to calculate the correct spectrum at the positions around focus. This was implemented in these measurements where the fluctuations on the BELLA focus were approximately 16 and 24 µm RMS for the horizontal and vertical axis respectively. This measurement is performed for multiple planes at and around best focus (here, around the image plane of the best focus). These data are then injected into an iterative phase-retrieval algorithm of the Gerchberg-Saxton type, to determine the spatial phase profile of the beam *at each frequency*.

The temporal scans of both techniques were designed so that the whole output spectrum of BELLA (43 nm FWHM centered around 800 nm [1]) could be retrieved with a satisfying resolution. This implied choosing a step size of 150 nm and a scan range of at least 50 µm.

Once the spatial amplitude and phase profiles at each frequency have been determined in a given plane (different for TERMITES and INSIGHT), the beam can be numerically propagated to any arbitrary longitudinal position, using plane wave decomposition [20] at each wavelength. All results presented in this article correspond to the front focal plane of the focusing parabola, i.e. to the plane located at a distance $f = 13.5$ m *before* this optic.

Like most of their relatives, these two techniques should actually be considered as spatio-spectral techniques, which determine $\tilde{E}(x, y, \omega) = \tilde{A}(x, y, \omega)\exp[i\phi(x, y, \omega)]$ up to an unknown spatially-homogeneous spectral phase. An additional measurement of the spectral phase at a single point of the beam is then required to lift this unknown, after which the spatio-temporal field $E(x, y, t)$ can be calculated by a Fourier-transform with respect to frequency. In the present case, this complementary information was provided by the GRENOUILLE device located in the output diagnostic area.

Both techniques require from a few hundreds to a thousand laser shots for each measurement, and all measurements were performed with the laser operated at 1 Hz. Three operation modes of the laser were used: 'low' power mode (225 mJ per pulse, 5.8 TW peak power) with the last amplifier turned off and wedge attenuation in the chain, moderate power mode (8.4 J per pulse, 215 TW peak power) in the same conditions but without attenuation, and full power (43 J per pulse, 1.1 PW peak power) with all amplifiers on and no attenuation before target.

We performed a set of tests with both techniques to verify that the attenuation and transport system does not lead to spatio-temporal distortions of the laser beam. First, to test the potential influence of nonlinear optical effects occurring along this system, we checked that the results of measurements performed at full power with Fourier-transform limited (FTL) pulses were identical to those obtained with pulses of the same energy, but
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Validation test of TERMITES and INSIGHT on the BELLA laser. Results derived from these two techniques (second and third columns, obtained after frequency-averaging of the spectrally-resolved measurement outputs) are compared to the ones provided by a standard Shack–Hartmann wavefront sensor (HASO by Imagine Optics, left column). The top row shows the frequency-averaged spatial intensity profile of the collimated BELLA beam, and the bottom row its frequency-averaged spatial phase profile. Note the small clip on the right edge of the beam: this is due to a small hole in the first wedge of the attenuation line, normally used to collect the accelerated electron beam in laser-wakefield experiments. This hole is usually placed at the center of the laser beam, but was moved to its edge for this measurement campaign. The resulting clip is well retrieved by TERMITES and INSIGHT.</figcaption></figure>

chirped in time to 1.2 ps by moving one of the compressor gratings (see supplementary document). For the consistency of this test, we also verified that chirping the pulse does not modify the spatio-temporal structure of the beam (due e.g. to undesired angular motion of the compressor grating), by comparing measurement results obtained at low power (thus safely avoiding any nonlinear effect in the attenuation system) for FTL and chirped pulses. Finally, at moderate power, we compared the results of measurements carried out with and without the gold foil, which showed that this foil does not introduce any additional coupling on the transmitted beam, but only a smooth spatially-homogeneous modification of the spectral envelope, corresponding to its theoretical transmission curve. With the support of all these tests, all results presented in this article correspond to the optimally-compressed BELLA beam at full power.

Validation tests of the measurement techniques themselves have already been presented, in [18] and [19] for TERMITES and INSIGHT respectively. The measurement configuration used in the present work offered an opportunity for a new type of test, which had not been possible before for practical reasons: with the spatio-temporal measurement devices separated from the Shack–Hartmann wavefront sensor (HASO by Imagine Optics) by only a few plane mirrors of small diameter, it is a safe assumption that the different beams arriving on these devices all have a similar wavefront. Therefore, we can safely compare the spatial intensity and phase profiles provided by this well-established instrument to the information provided by TERMITES and INSIGHT. Since this sensor has no spectral resolution, this is achieved by frequency-averaging the spectrally-resolved spatial profiles determined by these new techniques. This comparison, presented in figure 2, reveals a very good agreement, especially for the phase profiles, which is quite remarkable given the very different principles underlying these three techniques.

## Results and discussion

Spectrally-resolved spatial intensity and phase profiles of the BELLA laser at four different frequencies, measured at full power with INSIGHT, are displayed in figure 3. The complete evolution of these profiles as a function of frequency is displayed in supplementary movie 1, which is available online at stacks.iop.org/JPPHOTON/1/035001/mmmedia.

The first important conclusion is that BELLA does not suffer from major low-order chromatic phase aberration. Such aberrations can be particularly detrimental for ultrashort lasers, and are expected to be of two main types for collimated laser beams: (i) angular dispersion, typically resulting from stretcher or compressor misalignment, associated in the time domain to pulse front tilt (PFT), (ii) frequency-dependent wavefront curvature, typically induced by chromatic lenses, associated in the time domain to pulse front curvature (PFC).
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Frequency-resolved spatial properties of the collimated BELLA beam. Four frequencies $\omega_i$ are selected in the spatially-integrated spectrum (upper panel) and the corresponding spatial intensities $\tilde{A}^2(x, y, \omega_i)$ and phases $\phi(x, y, \omega_i)$ (lower row) are shown. These results have been obtained with INSIGHT.</figcaption></figure>

The magnitude of these couplings is observed to be very weak on BELLA, as clearly appears by considering the coefficients of the frequency-resolved Zernike decomposition of the spatial phase, displayed in figure 4. At focus, the variation of the lateral position of the focal spot, due to frequency-dependent horizontal and vertical wavefront tilts, is less than 10% of the beam waist across the laser spectrum. The wavefront curvature does exhibit a variation with frequency, slightly increasing on the red side of the spectrum (see also figure 3). Yet, at focus, this curvature results in a longitudinal focus shift that never exceeds 6% of the BELLA Rayleigh range. This variation is not linear with frequency, indicating a higher-order spatio-spectral aberration, as opposed to a simple PFC.

Figure 4 also shows that the results obtained with TERMITES for the low-order chromatic phase aberrations are qualitatively consistent with those provided by INSIGHT, especially when considering the fact that these measurements were carried out on different weeks. But TERMITES actually provides finer information on higher-order chromatic aberrations affecting the beam. This clearly appears in supplementary movie 2 (stacks.iop.org/JPPHOTON/1/035001/mmedia), which shows the complete evolution with frequency of the beam spatial phase, provided by TERMITES: on these phase plots, spatial patterns of higher-frequency are observed, which drift along the horizontal direction as frequency varies. This effect is illustrated more quantitatively in figure 5, which displays $(x, \omega)$ and $(y, \omega)$ slices of this spatio-spectral phase, in the horizontal and vertical planes, at the center of the beam ($x = y = 0$). The horizontal plane (x axis) corresponds to the dispersion plane of the compressor. At each wavelength, ‘high-frequency’ (5–6 periods across the full beam) modulations of the spatial phase are observed along x. These modulations linearly shift in space as frequency changes, creating an oblique spatio-spectral pattern in the $(x, \omega)$ space (figure 5(a)). By contrast, such an oblique spatio-spectral pattern is not observed at all in the $(y, \omega)$ space, i.e. in the plane normal to the compressor dispersion direction (figure 5(b)).

A simple interpretation of this feature has been recently analyzed in detail in [21]. In between the first and the last gratings of the compressor the laser beam is laterally chirped, i.e. the central position of the ‘beamlet’ associated to each frequency shifts approximately linearly with frequency. When this occurs, any modulation of the spatial phase induced on the beam, e.g. due to defects of a mirror or a grating, is imprinted at a different position relative to the center of each beamlet. After the compressor, this lateral chirp is removed, and all beamlets are again centered on the same axis. The spatial modulations induced within the compressor are then expected to shift approximately linearly with frequency, precisely as observed in figure 5. The expected slope of the pattern induced in this way can be deduced from the geometry of the compressor and the characteristics of its gratings, and is consistent with the observed effect (red dashed line in figure 5), thus strongly supporting this interpretation.

A simple analytical calculation (see supplementary document) shows that, if the corresponding spatial phase modulation is assumed to be perfectly sinusoidal at each frequency and of weak amplitude, then the scattered energy at focus forms two satellite spots located at a distance of $\lambda_0 f / \Lambda$ ($\Lambda$ is the spatial period of the modulations)
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Low-order coefficients (tilt and curvature) of the Zernike polynomials decomposition of the spectrally-resolved spatial phase retrieved with TERMITES (red) and INSIGHT (blue). Tilt values are converted to lateral displacement of the focal spot, relative to the beam waist. Equivalently, curvature values are expressed in longitudinal displacement of the focal spot, relative to the Rayleigh range of the focused beam. The displacements are calculated relative to the ideal transverse and longitudinal focal positions, i.e. the intersection of the optical axis with the focal plane of the parabola.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Spatio-spectral phase of the collimated BELLA beam, measured with TERMITES. Plots (a) and (b) respectively display (x, ω) and (y, ω) slices of this spatio-spectral phase, across the horizontal and vertical planes, at the center of the beam (see lines in image (c) showing a 2D phase profile at one frequency, with the spectrally-averaged wavefront subtracted). The red dashed line in (a) shows the estimated slope of the oblique spatio-spectral pattern that would be induced by spatial defects of optical elements located within the compressor, where the beam is laterally chirped. The origin of this line is arbitrary.</figcaption></figure>
on each side of the main focal spot, with one satellite pulse advanced in time and the other delayed with respect to the main pulse by $\pm 2\pi s/\Lambda$ ($s$ is the slope of the modulation). For the parameters seen in the measurements presented here ($\lambda_0 = 810$ nm, $f = 13.5$ m, $\Lambda \approx 4$ cm, $s = 0.445$ m fs rad$^{-1}$ or 1.31 mm nm$^{-1}$ as in the figure) this results in an offset of 273 $\mu$m and 71 fs in space and time respectively. With an amplitude of 0.4 rad peak-to-valley over most of the beam, the satellite foci only contain 2% of the beam energy, and are at least 100× less intense than the main pulse.

Because of the relatively low amplitude of these modulations, most of the associated signal falls within the background of images acquired around focus. As a result, they are not retrieved in the near-field when using the INSIGHT measurements, although this might in principle be possible by using a camera with an increased dynamic range. More generally, TERMITES tends to have a much higher spatial resolution than other techniques for reconstruction of collimated beams: this clearly appears when comparing the different profiles in figure 2. This is because resolving the spatial interference fringes generated by the TERMITES interferometer requires a very high spatial sampling of the collimated beam (above 2000 × 2000 pixels within the beam profile here), which is preserved in the final reconstruction of the beam. Use of a larger camera chip for INSIGHT measurements would also increase the resolution in the calculated near-field data, but this is generally impractical when small pixels are also necessary to resolve the focal spot.

To quantify the impact of these various phase aberrations on the laser performance at focus, we define and calculate two different spatio-temporal Strehl ratios of the beam (see the supplementary document for details). The first one, $SR_{Full}$, is the full spatio-temporal Strehl ratio: it is defined as the ratio of the actual peak intensity at focus, derived from the measured *E*-field, to the peak intensity that would be obtained for a beam with the same spatio-spectral amplitude, but with a perfectly constant spatio-spectral phase. This accounts for all phase aberrations present on the beam, both chromatic (i.e. frequency-dependent) and achromatic. We find $SR_{Full} = 0.62$ from TERMITES data, and $SR_{Full} = 0.88$ from the INSIGHT data. We note that the contribution of the oblique spatio-spectral pattern induced by the compressor (detected by TERMITES and not by INSIGHT) to this intensity reduction is actually very weak, so this cannot account for the different values of $SR_{Full}$ predicted by the two techniques.

The second Strehl ratio of interest, $SR_{STC}$, is defined as the ratio of the peak intensity at focus now calculated with the frequency-averaged wavefront distortions removed, to the peak intensity that would be obtained with a perfectly constant spatio-spectral phase. This corresponds to the Strehl ratio that would be achieved if the laser beam wavefront was corrected by a *perfect* adaptive optical system, removing all achromatic phase aberrations—while the deformable mirror actually used in the experiment still leaves some residual frequency-independent wavefront distortions (displayed in figure 2). This optimally-corrected beam would then still suffer from chromatic aberrations, which we detect in the present work. We find $SR_{STC} = 0.85$ for TERMITES, and $SR_{STC} = 0.92$ for INSIGHT. This shows that in the case of BELLA, the intensity reduction induced by purely chromatic effects is limited.

These calculations indicate that the significant difference in values of $SR_{Full}$ derived from the two techniques, mostly originates from differences in the frequency-averaged beam properties, rather than from chromatic effects. The TERMITES and INSIGHT measurements were conducted one week apart, and indeed measurements with a standard wavefront sensor (HASO) indicate a lower frequency-averaged Strehl ratio in the week of the TERMITES measurements (see supplementary document). However, this difference is too weak to fully account for the one found on $SR_{Full}$ between the TERMITES and INSIGHT results. This could be due to differences in the additional aberrations introduced in the transport line or measurement devices themselves, or to a more accurate calculation of the Strehl ratios from the TERMITES data due to the very high spatial resolution of this technique. More systematic and thorough measurement campaigns on smaller laser systems will be required to elucidate this point. In any case, along with the introduction of these calculations, the relatively good agreement and low severity of the $SR_{STC}$ values are key results of this work.

We have so far concentrated on the phase properties of the BELLA beam, but our measurements actually also reveal interesting effects on the spatio-spectral intensity: as observed in figure 3, the beam mode size and shape significantly vary with frequency. As the frequency increases, the beam diameter increases, and the beam profile turns to a ring shape. The fact that the clip imprinted on the right edge of the beam (see figure 2 and explanations in its caption) is observed on the blue side of the spectrum but not on the red one clearly indicates that this is not a measurement artifact. We also emphasize that the reconstruction of such fine features in the near-field beam profile from INSIGHT data is quite remarkable, since the measurement was actually performed around focus, and the beam was then numerically propagated to a plane far from focus. This provides a strong validation of the spatial phase retrieved at focus with INSIGHT.

We can use simple numerical modeling to show that this frequency-dependent intensity profile is an intrinsic feature of high-power CPA lasers, which had not been clearly identified so far, to the best of our knowledge. To this end, we have modeled the amplification of the BELLA beam in the successive amplifiers, using the Frantz–Nodvik equations and considering a positively chirped laser beam (i.e. lower frequencies arrive
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Spatio-spectral intensity of the BELLA beam. The left image displays the (x, ω) slice of this spatio-spectral intensity, across the horizontal plane, at the center of the beam (x = y = 0), deduced from simulations of the CPA process based on the Frantz-Nodvik equations. The right image shows the same quantity, now derived from INSIGHT measurements at full power. To clearly highlight the spectral evolution of the beam spatial profile in these plots, the total energy in the beam at each frequency has been normalized to 1 (i.e. $\int dxdy\tilde{A}^2(x, y, \omega) = 1$ for all ω within the laser spectrum).</figcaption></figure>

first). Details of these simulations are provided in the supplementary document, and essentially start with a measured spectrum after the XPW stage (where the spectral phase is flat), and the spectral phase imparted by the stretcher immediately after. This beam is then amplified in the successive five amplifier stages, with the pump and seed modeled in space, and the full spectral gain cross-section taken into account. In figure 6, the (x, ω) spatio-spectral intensity profile of the beam predicted by these simulations for full power operating conditions (panel a) is compared to the one deduced from our INSIGHT measurements (panel b, same data as in figure 3), revealing a good qualitative agreement.

These simulations can then be exploited to get a simple interpretation of this effect: in the CPA scheme, the lower frequencies arrive first in the amplification medium, and are more efficiently amplified in the center of the beam. These frequencies deplete the amplification medium before the whole chirped pulse has been amplified. As a result, higher frequencies, which arrive later, are more efficiently amplified on the beam edges than at its center. This qualitatively explains the spectral evolution of the beam size, and the ring shape observed at high frequencies.

## Conclusion

In conclusion, we have reported the first spatio-temporal characterization of a PW femtosecond laser, carried out with two independent measurement techniques, with the laser focused at full power in the experimental chamber. These measurements show that the BELLA laser does not suffer from major chromatic aberrations that would severely degrade its performance. We estimate the reduction in peak intensity resulting merely from spatio-temporal couplings to about 10%, thus demonstrating that the CPA technique is still suitable and manageable at these extreme scales. We note however that significant chromatic effects start appearing on the red edge of the spectrum, suggesting that the beam optimization would possibly get more challenging for lasers with broader spectra, such as PW lasers now under construction in different laboratories worldwide.

These measurements provide very fine information on the beam properties and the laser system operation. In particular, we have observed small spatio-spectral modulations that have recently been predicted [21] to result from defects of optics within the grating compressor, as well as effects resulting from the temporal dynamics of the gain during the amplification of the chirped laser pulses. As ultrahigh intensity laser-plasma interaction experiments put more and more stringent requirements on the control of the laser parameters, such information could be used in the future for extremely fine optimization of high-power laser systems.

Laser-matter interaction experiments should ultimately greatly benefit from this new degree of knowledge of ultrashort laser beams. For instance, in the case of BELLA, the full spatio-temporal profile of the laser beam can be calculated in any arbitrary plane (see supplementary movie 3), and this can be used as an input for a realistic description of the laser field in 3D Particle-In-Cell simulations. This opens the way to a deeper understanding of
laser-wakefield electron acceleration experiments on this system, which require guiding of the beam in plasma channels of many centimeters in length. This could ultimately contribute to an unprecedented optimization of laser-plasma accelerators.

## Acknowledgments

The authors gratefully acknowledge Amplitude Laser Group for their support. We also thank Chris Pieronek for helping with the experiment, Donald Syversrud, Zachary Eisentraut, David Evans, Nathan Ybarrolaza and Greg Mannino for technical support.

## Funding information

The research leading to this work has been funded by the ERC (grant ExCoMet number 694596), by Investissements d'Avenir LabEx PALM (projects EXYT and IMAPS), by the CEA's DFR Impulsion program, and has received financial support from Amplitude Laser Group. The work at the BELLA facility was supported by the Director, Office of Science, Office of High Energy Physics of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231.

## ORCID iDs

Antoine Jeandet https://orcid.org/0000-0003-4549-3067

Kei Nakamura https://orcid.org/0000-0001-9842-7114

Spencer W Jolly https://orcid.org/0000-0002-5783-2081

Fabien Quéré https://orcid.org/0000-0002-4343-9780

## References

[1] Nakamura K, Mao H S, Gonsalves A J, Vincenti H, Mittelberger D E, Daniels J, Magana A, Toth C and Leemans W P 2017 *IEEE J. Quantum Electron.* **53** 1200121

[2] Sung J H, Lee H W, Yoo J Y, Yoon J W, Lee C W, Yang J M, Son Y J, Jang Y H, Lee S K and Nam C H 2017 *Opt. Lett.* **42** 2058–61

[3] Zeng X *et al* 2017 *Opt. Lett.* **42** 2014–7

[4] Clayton C E *et al* 2010 *Phys. Rev. Lett.* **105** 105003

[5] Wang X *et al* 2013 *Nat. Commun.* **4** 1988

[6] Leemans W P *et al* 2014 *Phys. Rev. Lett.* **113** 245002

[7] Kim H T *et al* 2017 *Sci. Rep.* **7** 10203

[8] Gonsalves A J *et al* 2019 *Phys. Rev. Lett.* **122** 084801

[9] Marklund M and Shukla P K 2006 *Rev. Mod. Phys.* **78** 591–640

[10] Strickland D and Mourou G 1985 *Opt. Commun.* **56** 219–21

[11] Bor Z 1989 *Opt. Lett.* **14** 119

[12] Bowlan P, Gabolde P, Shreenath A, McGresham K, Trebino R and Akturk S 2006 *Opt. Express* **14** 11892–900

[13] Cousin S L, Bueno J M, Forget N, Austin D R and Biegert J 2012 *Opt. Lett.* **37** 3291

[14] Wyatt A S, Walmsley I A, Stibenz G and Steinmeyer G 2006 *Opt. Lett.* **31** 1914

[15] Gabolde P and Trebino R 2006 *Opt. Express* **14** 11460–7

[16] Dorrer C and Bahk S W 2018 *Opt. Express* **26** 33387

[17] Miranda M, Kotur M, Rudawski P, Guo C, Harth A, L'Huillier A and Arnold C L 2014 *Opt. Lett.* **39** 5142–4

[18] Pariente G, Gallet V, Borot A, Gobert O and Quéré F 2016 *Nat. Photon.* **10** 547–53

[19] Borot A and Quéré F 2018 *Opt. Express* **26** 26444–61

[20] Goodman J W 2005 *Introduction To Fourier Optics* 3rd edn (Englewood, CO: Roberts & Company)

[21] Li Z, Tsubakimoto K, Yoshida H, Nakata Y and Miyanaga N 2017 *Appl. Phys. Express* **10** 102702