Photonics Research, Volume. 13, Issue 3, 761(2025)

Compact, folded multi-pass cells for energy scaling of post-compression

Arthur Schönberg1、*, Supriya Rajhans1,2, Esmerando Escoto1, Nikita Khodakovskiy1, Victor Hariton1, Bonaventura Farace1, Kristjan Põder1, Ann-Kathrin Raab3, Saga Westerberg3, Mekan Merdanov3, Anne-Lise Viotti3, Cord L. Arnold3, Wim P. Leemans1, Ingmar Hartl1, and Christoph M. Heyl1,4,5
Author Affiliations
  • 1Deutsches Elektronen-Synchrotron DESY, 22607 Hamburg, Germany
  • 2Friedrich-Schiller-Universität Jena, 07743 Jena, Germany
  • 3Department of Physics, Lund University, SE-221 00 Lund, Sweden
  • 4GSI Helmholtzzentrum für Schwerionenforschung GmbH, 64291 Darmstadt, Germany
  • 5Helmholtz-Institut Jena, 07743 Jena, Germany
  • show less

    Combining high peak power and high average power has long been a key challenge of ultrafast laser technology, crucial for applications such as laser-plasma acceleration and strong-field physics. A promising solution lies in post-compressed ytterbium lasers, but scaling these to high pulse energies presents a major bottleneck. Post-compression techniques, particularly Herriott-type multi-pass cells (MPCs), have enabled large peak power boosts at high average powers but their pulse energy acceptance reaches practical limits defined by setup size and coating damage threshold. In this work, we address this challenge and demonstrate, to our knowledge, a novel type of compact, energy-scalable MPC (CMPC). By employing a novel MPC configuration and folding the beam path, the CMPC introduces a new degree of freedom for downsizing the setup length, enabling compact setups even for large pulse energies. We experimentally and numerically verify the CMPC approach, demonstrating post-compression of 8 mJ pulses from 1 ps down to 51 fs in atmospheric air using a cell roughly 45 cm in length at low fluence values. Additionally, we discuss the potential for energy scaling up to 200 mJ with a setup size reaching 2.5 m. Our work presents a new approach to high-energy post-compression, with up-scaling potential far beyond the demonstrated parameters. This opens new routes for achieving the high peak and average powers necessary for demanding applications of ultrafast lasers.

    1. INTRODUCTION

    Ultrafast laser technology has experienced immense progress within recent years. Ultrashort, high-peak-power lasers are used in a vast range of applications, including attosecond science and high-harmonic generation [14], laser-plasma acceleration [5,6], or high-field science including laser-based nuclear fusion [7]. However, developing a laser source that is simultaneously average and peak power scalable remains a major challenge.

    The invention of mode-locked solid-state laser technology, in particular titanium-doped sapphire (Ti:Sa) lasers in combination with chirped-pulse amplification (CPA), enabled ultrashort, few-cycle pulses with unprecedented pulse energy [810]. Nowadays, peak powers exceeding the terawatt regime are routinely employed [9]. While excelling in peak power performance, Ti:Sa amplifiers are commonly constrained in average power to a few tens of watts, which can be attributed to their large quantum defect [11]. As an alternative to laser amplification in active gain media, optical parametric processes can be employed. In particular optical parametric chirped-pulse amplifiers (OPCPAs) offer broad bandwidths supporting few-cycle pulses and simultaneously high average powers [10,12]. However, OPCPA systems suffer from low pump-to-signal efficiencies typically around 10%–20% for pulses in the range of 10s of femtoseconds (fs) [13]. Ultrafast ytterbium (Yb)-based laser architectures on the other hand provide excellent average power scalability exceeding 10 kW [14], but pulse durations limited to 100s of femtoseconds up to about 1 picosecond (ps). Combining Yb lasers with efficient post-compression methods supporting large (>10) compression factors and high pulse energies can offer an excellent solution to the power scaling challenge.

    In recent years, a number of post-compression techniques have been developed, mostly relying on self-phase modulation (SPM) as the nonlinear process for spectral broadening [15]. In particular, gas-based technologies provide excellent tools for post-compression of high-power lasers. Example systems rely on gas-filled hollow-core fibers (HCF) [1618], cascaded focus and compression (CASCADE) [19], white-light filaments [20], as well as Herriott-type multi-pass cells (MPCs) [15,2123]. In HCFs, post-compression of 70 mJ 220 fs pulses down to 30 fs has been demonstrated in a 3 m long fiber [24]. Post-compression of very high pulse energies in the multiple Joule range has been achieved via thin-film spectral broadening techniques. However, typical compression factors lie in the range of only two to five [25,26]. Similar to HCFs, MPCs enable large compression factors reaching 10–20 or more while supporting a wide range of pulse energies. In addition, MPCs support high average powers [27] and outperform HCFs in system footprint especially for large compression factors [15]. The maximum attainable energy acceptance in a standard, two-mirror MPC is directly proportional to its size [15,28]. A record of 200 mJ has been achieved in a 10 m long MPC setup [29]. Further energy up-scaling leads to MPC sizes that are impractical for standard laboratory settings. The development of a highly efficient post-compression method supporting large compression factors and high pulse energies thus remains a key challenge.

    We here introduce a new MPC type, the compact MPC (CMPC), which possesses a weakly focused fundamental mode as well as a linear beam pattern on the focusing mirrors. This geometry allows us to fold the beams inside the MPC using two additional planar mirrors, thus introducing a new energy scaling parameter, the folding ratio Γ. The CMPC in principle allows for an arbitrary amount of folding and thus, very compact setup sizes while sharing key properties of standard MPCs such as high average power support, excellent beam quality, and efficiency. We experimentally demonstrate spectral broadening of 1030 nm, 8 mJ, 1 ps pulses in a CMPC in atmospheric air using a setup with an effective length of around 45 cm. We keep the maximum mirror fluence at a moderate level of around 170  mJ/cm2 and demonstrate compressibility of 1 ps input pulses down to 51 fs with an MPC throughput reaching 89% while maintaining excellent spatio-temporal pulse characteristics.

    2. CONCEPT

    Most MPCs demonstrated for post-compression to date rely on two identical concave mirrors, resembling the most basic optical cavity arrangement. However, more complex designs employing a convex mirror and a concave mirror [30,31] or even multiple additional mirrors can provide advantageous mode-forming capabilities. MPCs with more than two mirrors have been proposed in previous works focusing on energy scaling of MPCs [31,32]. Other energy scaling methods that do not rely on modifications on the MPC geometry include divided pulse nonlinear post-compression [33,34], or using higher-order Laguerre beams to decrease the fluence on the mirrors [35]. These techniques can be combined with geometrical scaling methods. With the CMPC we demonstrate a new geometrical approach. The concept of the CMPC is based on a weakly focused beam and folding of the beam path via multiple reflections on two additional, planar mirrors in each pass through the cell, as shown in Fig. 1. This provides an additional tuning parameter, namely, the folding ratio Γ, which reduces the length of the CMPC by LeffL/Γ. Figures 1(A) and 1(B) depict the principle of beam folding and the effective size reduction of the cell together with the configuration regimes for standard MPCs and the CMPC. In Fig. 1(B) the original beam path is divided into three parts and then folded such that the setup length reduces by approximately a factor of three and thus Γ=3.

    MPC and CMPC schematic and configuration regimes. (A) Standard MPC mode. (B) CMPC mode without and with folding. In the unfolded geometry the CMPC is typically longer than a comparable MPC, reducing to much shorter length Leff after folding. (C) Typical configuration regimes and corresponding fluence characteristics of MPC and CMPC within the stability range.

    Figure 1.MPC and CMPC schematic and configuration regimes. (A) Standard MPC mode. (B) CMPC mode without and with folding. In the unfolded geometry the CMPC is typically longer than a comparable MPC, reducing to much shorter length Leff after folding. (C) Typical configuration regimes and corresponding fluence characteristics of MPC and CMPC within the stability range.

    In general, the geometry of a symmetric Herriott-type MPC—including a CMPC—is fully determined by three of the following four parameters: radius-of-curvature (ROC) of the mirrors R, the propagation length between the two focusing mirrors L, the number of round-trips N, and the configuration parameter k. These parameters are related by [15] LR=1cos(πkN).Typically, k is chosen to be integer and prime to the number of round-trips N, in order to ensure the re-entrant condition of the MPC. Moreover, the parameter k determines the angular advance of the beam on the MPC mirror, which is equivalent to the Gouy phase per pass ϕG(1)=πk/N [15]. The amount of Gouy phase per pass determines the caustic of the beam in the MPC. For a standard MPC, ϕG(1) approaches π, which leads to a large mode size wm on the mirrors and a small waist w0 in the focus. In the CMPC, the beam propagates within the Rayleigh-range, where the accumulated Gouy phase per pass is ϕG(1)<π/2. This is the case when L/R<1, where the beam radius at the mirror wm is comparable to the beam radius at focus w0. The peak fluence of the beam follows a similar behavior, with F0/Fm<2 for L/R<1 [Fig. 1(C)]. This weak focusing geometry with small fluence variations enables the placement of additional optical components within the beam path of the CMPC and thus folding of the beam. In addition, a weakly focused mode eliminates pulse energy limitations arising due to ionization in standard gas-filled MPCs. Figure 2 illustrates the CMPC scheme. In Fig. 2(A) the setup schematic of a standard MPC with linear pattern alignment is shown. Here, the strongly focused beams propagate directly between the two focusing mirrors FM1 and FM2, without any optical components in-between. In the CMPC [Fig. 2(B)], two additional planar mirrors PM1 and PM2 are placed behind FM1 and FM2. As seen in Fig. 2(B-ii), the beam propagates from FM1 to PM2 and PM1, where it is folded multiple times before reaching FM2. Here, the depicted folding ratio is Γ=9. The total path length between FM1 and FM2 is the length determined by Eq. (1). However, the effective length is now reduced by LeffL/Γ, which under the correct choice of R and Γ can become significantly shorter than that of a standard MPC with similar pulse energy acceptance.

    Schematics of standard MPC and CMPC. (A) Standard MPC geometry and beam pattern considering a linear beam pattern alignment (FM1/2: focusing mirrors). View of the beam path from the front (i) and from the top (ii). (iii) Resulting patterns for circular/linear alignment on FM1/2. (B) CMPC geometry and beam patterns. (i) Front view with additional planar mirrors (PM1/2) behind FM1/2. (ii) Top view onto the CMPC with beam folding ratio Γ=9. The numbers (1–5) are indicating the order of reflections of the first pass. (iii) Beam patterns on the mirrors FM1/2 and PM1/2. The black line follows the beam reflections 1–5 on the mirrors. (C) Photographs of experimental beam patterns with N=11, k=3, and Γ=25. POM: pick-off mirror.

    Figure 2.Schematics of standard MPC and CMPC. (A) Standard MPC geometry and beam pattern considering a linear beam pattern alignment (FM1/2: focusing mirrors). View of the beam path from the front (i) and from the top (ii). (iii) Resulting patterns for circular/linear alignment on FM1/2. (B) CMPC geometry and beam patterns. (i) Front view with additional planar mirrors (PM1/2) behind FM1/2. (ii) Top view onto the CMPC with beam folding ratio Γ=9. The numbers (1–5) are indicating the order of reflections of the first pass. (iii) Beam patterns on the mirrors FM1/2 and PM1/2. The black line follows the beam reflections 1–5 on the mirrors. (C) Photographs of experimental beam patterns with N=11, k=3, and Γ=25. POM: pick-off mirror.

    The pulse energy acceptance of a gas-filled MPC can be estimated considering the mirror fluence and the focus intensity. The beam waist radius in the focus and on the mirrors in an MPC can be calculated using the equations w02=(Rλ/2π)sin(π/kN) and wm2=(Rλ/π)tan(πk/2N), respectively, where λ is the wavelength of the laser [32]. The parameters L/R can be directly mapped to k/N according to Eq. (1). Assuming a laser with Gaussian pulses and pulse energy E, the maximum peak fluence in the focus F0=2E/πw02, as well as the peak fluence on the mirror Fm=2E/πwm2 can be calculated [32]. For a standard MPC consisting of two identical focusing mirrors, the fluence on these mirrors should not exceed the threshold fluence Fth, yielding a maximum pulse energy of EMPCRλFth2tan(πk/2N).For simplicity, the second energy limitation that can arise at high focus intensity due to ionization is omitted in this discussion as it has little or no relevance for the CMPC. In a CMPC, where beam folding is achieved by additional mirrors (Fig. 2) the beam can be reflected at almost any point between the focusing mirrors including the focus. Therefore, the condition for avoiding damage needs to be F0Fth and thus ECMPCRλFth4sin(πk/N).In addition, replacing R with L using Eq. (1) and writing L in terms of LeffL/Γ, we obtain ECMPCλFthΓLeff4tan(πk/2N).

    Equations (2)–(4) illustrate the energy scaling possibilities of MPCs. For standard MPCs [Eq. (2)], the energy scaling options are restricted to maximizing the ratio k/N1, increasing the ROC of the mirrors R (which increases L proportionally at constant k/N), or increasing the threshold fluence Fth. In the case of CMPCs, energy scaling can be achieved differently. Since k/N is generally smaller, the energy acceptance for the same set of parameters R and Fth is typically reduced according to Eq. (3). Thus, the first step of scaling is achieved by adjusting R in Eq. (3) such that ECMPC=EMPC. This typically increases L, i.e., LCMPC>LMPC. The second step is to fold the beam by the factor Γ and decrease the length of the system by exploiting Leff=L/Γ.

    To illustrate the discussion with an example, we use a typical set of parameters for a standard MPC with R=1  m and a λ=1030  nm, 1 ps level laser. The threshold fluence for quarter-wave-stack high-reflectance dielectric coatings can, e.g., be set to Fth=0.17  J/cm2 (corresponding to the linear fluence in our experiment), leaving headroom to damage. With N=15 round-trips and k=14, we arrive at EMPC8.3  mJ and a length of close to 2 m. To match the pulse energy acceptance of the CMPC to the MPC, we first set k=4, and increase the ROC to R=25  m and the length to L=8.27  m. Now folding the beam Γ=25 times, we arrive at ECMPC8.1  mJ with an effective length of Leff=33  cm. A second example is illustrated in Fig. 3. Here we refer to the works reported by Pfaff et al., where 200 mJ pulses have been compressed in a 10 m long MPC, with a fluence on the mirrors of approximately 0.5  J/cm2 [29]. Figure 3 shows that the energy acceptance increases with larger folding ratio and effective system length. At 10 m length, a folding ratio of Γ=25 supports an energy of 700 mJ and at Γ=49, the accepted energy exceeds 1 J. At these parameters, the CMPC with 1 J energy acceptance and Γ=49 would require an ROC of R>1  km with a focusing mirror diameter of 60 cm and a planar mirror size of 60  cm×68  cm [see Eqs. (A3)–(A5) in Appendix A]. In particular the large ROCs can become a limiting factor, as discussed in Section 4. Conversely, for a pulse energy of 200 mJ, the CMPC size can be scaled down to about 3 m or 1.5 m considering Γ=25 and Γ=49, respectively. As visible in Fig. 3, for Γ=1, the fully “unfolded” CMPC is generally longer than the standard MPC considering an identical threshold fluence for both cases. This is due to the fact that the CMPC operates in the L/R<1 regime while the MPC operates at L/R2.

    Energy scaling characteristics: calculated pulse energy acceptance Emax as a function of effective length Leff. The solid black line indicates a standard MPC as reported in Ref. [29] [Eq. (2)]. Each colored solid line represents a CMPC with a folding ratio Γ [Eq. (4)]. The gray line indicates an effective length of 10 m. The dashed lines show the corresponding pulse energies for selected values of Γ. We set the threshold fluence to Fth=0.5 J/cm2.

    Figure 3.Energy scaling characteristics: calculated pulse energy acceptance Emax as a function of effective length Leff. The solid black line indicates a standard MPC as reported in Ref. [29] [Eq. (2)]. Each colored solid line represents a CMPC with a folding ratio Γ [Eq. (4)]. The gray line indicates an effective length of 10 m. The dashed lines show the corresponding pulse energies for selected values of Γ. We set the threshold fluence to Fth=0.5  J/cm2.

    A second limitation in MPCs, which we omitted in the above discussion, typically arises due to the peak intensity of the laser exceeding the ionization threshold of the gas in the MPC [15]. Due to the loose focusing geometry of the CMPC however, this problem is fully avoided here, as the limitation due to the LIDT of the coating typically sets in much earlier. For example, in the above discussed system with k/N=4/15 and R=25  m, ionization of a gas with an ionization threshold of Ith=1×1013  W/cm2 would set in at 478 mJ for 1 ps pulses [32].

    3. EXPERIMENT

    Our experimental setup displayed in Fig. 4 uses a commercial innoslab Yb laser system (AMPHOS GmbH), delivering 1030 nm 1 ps pulses at a repetition rate of 1 kHz with a pulse energy of >8  mJ. After suitable mode-matching, the beam is sent onto the pick-off mirror POM1 and into the CMPC [Fig. 4(B)]. The CMPC consists of two concave 4 inch mirrors with R=25  m and two planar folding mirrors with size 10  cm×10  cm. The configuration is set to N=11 round-trips and k=3, which corresponds, according to Eq. (1), to an unfolded MPC length of L=8.63  m. Following in-coupling at the pick-off mirror POM1, the beam is sent onto the planar folding mirror PM2 and is subsequently reflected by PM1, with both mirrors being aligned in a V-shaped configuration and separated by about 34 cm. The beam continues its path as described in Fig. 2, forming the typical CMPC pattern. A folding ratio of Γ=25 results in an effective CMPC length of Leff=34.5  cm. However, due to spatial constraints, the two focusing mirrors FM1 and FM2 are placed a few centimeters behind PM1 and PM2, making the total CMPC length slightly larger, resulting in around 45 cm [see Fig. 4(B)]. After propagating N=11 round-trips (22 passes), the beam returns on PM1 and is coupled out with a different angle in the vertical direction. The total number of mirror reflections amounts to 2NΓ=550 and the total propagation length is roughly 190 m. We use ambient air at atmospheric pressure (1006 mbar) as the nonlinear medium for spectral broadening. Nevertheless, a chamber or housing is necessary to avoid beam fluctuations due to turbulences in air.

    Experimental setup. The CMPC input beam is indicated in red, and the output beam in orange. FM: focusing mirror. PM: planar mirror. POM: pick-off mirror. Min: input/output mirror. MM: mode-matching. DCM: double-chirped mirror. CCD: camera. z0: reference point for stability measurements. TFP: thin-film polarizer. wdg: wedge. (A) Beam input, telescope, and diagnostics used in the experiment. (B) CMPC setup with vacuum chamber.

    Figure 4.Experimental setup. The CMPC input beam is indicated in red, and the output beam in orange. FM: focusing mirror. PM: planar mirror. POM: pick-off mirror. Min: input/output mirror. MM: mode-matching. DCM: double-chirped mirror. CCD: camera. z0: reference point for stability measurements. TFP: thin-film polarizer. wdg: wedge. (A) Beam input, telescope, and diagnostics used in the experiment. (B) CMPC setup with vacuum chamber.

    After out-coupling, the spectrally broadened beam is separated from the input beam with another pick-off mirror POM2 [Fig. 4(A)]. Here, we measure an output power of around 7.1 W and thus a total transmission of roughly 89%. This corresponds to a reflectance of the mirrors of at least 99.98% per reflection, disregarding clipping losses. Subsequently, a wedge is used for beam sampling, reflecting around 8%, corresponding to roughly 570 mW. The transmission of the wedge is sent onto a power meter. The wedge-reflected pulse is compressed using a chirped-mirror compressor with 26 reflections corresponding to about 5200  fs2 of compensated second-order dispersion.

    We analyze the compressed pulses using frequency-resolved optical gating (FROG) and measure the spectrum. In addition, the beam quality parameter M2 is analyzed using the same wedge reflection. We furthermore record input and output near-field pointing using the input/output mirror (Min) transmission of both input and output beams. For direct comparison of input and output, both cameras are placed in the exact same distance d=42  cm from a common reference point z0, located at the vacuum chamber input window. The spectrum and FROG results are shown in Fig. 5. The input pulse has a pulse duration of about 1.1 ps with a Fourier-transform-limit (FTL) of about 1 ps. After spectral broadening in the CMPC, the output pulse FTL reaches 50 fs. Following compression, we reach 51 fs [Fig. 5(A)]. We simulate the spectral broadening process using the measured input pulse [Fig. 5(A)]. The simulations are conducted using our in-house developed (2 + 1)D radially symmetric simulation code based on Hankel-transforms in the spatial dimension, which is described in Appendix A [Eqs. (A1) and (A2)]. The simulated output pulse agrees very well with the measured output pulse, exhibiting similar temporal characteristics after compression. In both measurement and simulation, we observe a temporal pedestal likely stemming from uncompressed spectral components in the longer wavelength range. On the trailing edge, a post-pulse appears at around 800 fs, which we attribute to a pulse-breakup caused by the delayed response of the Raman-Kerr contribution in air. The measured and simulated spectra [Fig. 5(B)] also show similar characteristics as well as similar broadening, except for the strength of the side lobes compared to the background close to 1030 nm.

    Measured and simulated pulse characteristics after CMPC spectral broadening and subsequent compression in 1 bar of air. (A) Reconstructed input and output temporal pulse. (B) Measured input spectrum and broadened spectrum. The inset shows the measured broadened spectrum along with simulated spectrum. (C) Measured FROG trace (FT). (D) Retrieved FROG trace (FROG error=0.73% on a 1024×1024 grid).

    Figure 5.Measured and simulated pulse characteristics after CMPC spectral broadening and subsequent compression in 1 bar of air. (A) Reconstructed input and output temporal pulse. (B) Measured input spectrum and broadened spectrum. The inset shows the measured broadened spectrum along with simulated spectrum. (C) Measured FROG trace (FT). (D) Retrieved FROG trace (FROG error=0.73% on a 1024×1024 grid).

    The molecular nature of air results in a broadened spectrum that qualitatively deviates from a typical “SPM-like” spectrum, obtained from atomic gases. We therefore conduct spectral broadening measurements in Krypton as well, in order to demonstrate typical SPM broadening and compare it with the simulations. Due to technical constraints originating from mechanical instabilities of the vacuum chamber, the chamber could only be operated at 1 bar pressure. We measure the spectral broadening with 6.5 mJ pulse energy and the results are shown in Fig. 6, where we reach an FTL of approximately 120 fs.

    Spectral broadening of 6.5 mJ pulses in 1 bar of Krypton. CMPC input (gray line) and measured output (dark blue line) spectra are shown. The inset shows the simulated spectrum along with the measured output spectrum.

    Figure 6.Spectral broadening of 6.5 mJ pulses in 1 bar of Krypton. CMPC input (gray line) and measured output (dark blue line) spectra are shown. The inset shows the simulated spectrum along with the measured output spectrum.

    The beam path inside the CMPC involves 550 reflections and a long propagation path of about 190 m. The beam is particularly sensitive to displacements of the folding mirrors PM1 and PM2 (Fig. 4), where most of the reflections take place. Thus, it is important to characterize the stability of the beam. Figures 7(A) and 7(B) show the measured beam position stability for both input and output beams. We conduct the measurements separately within a few minutes. The measurements show only a very slight increase in the RMS beam position deviation of approximately 5% in the output, compared to the input. This demonstrates the spatial stability owing to the input-to-output imaging property of the CMPC, also known from standard MPCs. In addition, we measure the spectral stability of the system over 5 min [Fig. 7(C)]. We measure RMS fluctuations of the FTL of 0.33 fs, with a mean value of 50.9 fs, resulting in a relative RMS deviation of 0.65%, indicating that the spectral broadening is stable.

    Beam position and spectral stability, measured in 1 bar of air. (A), (B) Scatter-plots of the beam position of input (A) and output (B) beams. Measurements are taken separately; the data span for each measurement is 7.5 min and 4000 data points, with a few minutes between measurements. In both measurements, the camera is placed at exactly 42 cm distance from the common reference point z0 (Fig. 4). Circles indicate deviations of 50, 100, and 150 μm from the mean value. The measurements are in accordance with the medium-term beam position stability measurements defined in EN ISO 11670:2003. (C) Measured broadened spectrum recorded over 5 min with 10 points per minute at full broadening (logarithmic scale). The RMS deviation of the FTL amounts to 0.65%.

    Figure 7.Beam position and spectral stability, measured in 1 bar of air. (A), (B) Scatter-plots of the beam position of input (A) and output (B) beams. Measurements are taken separately; the data span for each measurement is 7.5 min and 4000 data points, with a few minutes between measurements. In both measurements, the camera is placed at exactly 42 cm distance from the common reference point z0 (Fig. 4). Circles indicate deviations of 50, 100, and 150 μm from the mean value. The measurements are in accordance with the medium-term beam position stability measurements defined in EN ISO 11670:2003. (C) Measured broadened spectrum recorded over 5 min with 10 points per minute at full broadening (logarithmic scale). The RMS deviation of the FTL amounts to 0.65%.

    We further measure the beam quality parameter M2 of the input beam after the mode-matching telescope and the output beam. We observe a beam quality degradation from M2=1.36 at the input to 1.76 at the output at full power. As this observation may indicate spatio-temporal coupling (STC) effects, we decided to further investigate the spatio-spectral pulse characteristics using spatially resolved Fourier-transform spectrometry [36]. The STC measurements provide information about the spatially dependent spectral homogeneity of the pulse V(x,y), which we calculate via the spectral overlap integral as defined in Eqs. (A6) and (A7) in Appendix A. An average spectral overlap Vavg is then computed by weighting the V(x,y) with the fluence F(x,y) obtained from the same dataset and averaging over an area defined by the measured 1/e2 diameter of the beam. In air (Fig. 12, Appendix A), we measure Vavg=90.6% at full power, confirming a slight degradation of the spatio-spectral homogeneity compared to the input beam, which exhibits an almost perfect homogeneity of Vavg99%. We further investigate more generally how the spatio-spectral homogeneity behaves in atomic gases as a function of the input peak power P of the pulse, compared to the critical power Pc of the gas. The measurement is performed using Krypton as the nonlinear medium at 1 bar pressure and the same CMPC configuration is used as described in the beginning of this section. We carry out STC measurements at four different pulse energies and thus four different values for P/Pc and measure an STC trace for each point. Figure 8 summarizes the STC measurement results. From the first measurement point [Fig. 8(A)] at P/Pc0.4 to the third at P/Pc0.6, we observe a linear decline from 98.3% to 90.5%. The fourth point at P/Pc0.65 exhibits a stronger deterioration. We compare our results with simulations considering again the measured input pulse as input for the simulation. Here we observe an onset of homogeneity reduction [Fig. 8(A)] at approximately P/Pc=0.55 and an overall weaker deterioration compared to the measurements. However, the simulations reproduce the measured behavior qualitatively including the onset of homogeneity degradation at around P/Pc0.6 well. At P/Pc0.8, the simulations predict a beam collapse causing a sharp decline of Vavg. The faster degradation of the spatio-spectral homogeneity in the experimental data might be related to imperfect input beam characteristics in the experiment and/or possible spatial phase distortions arising due to many beam reflections on the CMPC mirrors.

    Measured spatio-spectral beam characteristics. (A) Measured and simulated average spectral overlap Vavg and simulated accumulated nonlinear phase ϕNL(1) as a function of peak power P normalized to the critical power Pc. (B) Example data set of a single measurement point, indicated by a circle in (A). (i), (ii) Spectral distribution of the pulse over the spatial x-direction, integrated over the y-direction and vice versa. (C) Corresponding STC traces for a standard MPC with similar broadening. wx and wy represent the 1/e2 beam radii.

    Figure 8.Measured spatio-spectral beam characteristics. (A) Measured and simulated average spectral overlap Vavg and simulated accumulated nonlinear phase ϕNL(1) as a function of peak power P normalized to the critical power Pc. (B) Example data set of a single measurement point, indicated by a circle in (A). (i), (ii) Spectral distribution of the pulse over the spatial x-direction, integrated over the y-direction and vice versa. (C) Corresponding STC traces for a standard MPC with similar broadening. wx and wy represent the 1/e2 beam radii.

    For an example point at P/Pc0.6, the spatio-spectral distribution obtained via an STC scan is shown in Fig. 8(B). We compare the results with a spectrally broadened pulse from a standard gas-filled MPC with N=17, k=16, and R=1  m at similar spectral broadening characteristics. The corresponding spatio-spectral distribution is shown in Fig. 8(C). Both measurements indicate rather homogeneous spatio-spectral characteristics and thus a similar spectral broadening performance for both MPC types, with Vavg=96.4% for the standard MPC and a slightly reduced Vavg=90.5% for the CMPC. We can further compare our measured Vavg values to another non-standard MPC, where a convex-concave MPC arrangement was used for spectral broadening. In Ref. [31], a value of Vavg=98% was achieved, albeit with a lower spectral broadening factor corresponding to compression from 260 fs to 50 fs (spectral broadening factor of approximately five) and in a solid medium. Another experiment using a convex-concave MPC and a solid medium yielded Vavg=90% [30], with a spectrum that supports 133 fs, starting from 670 fs. In the latter example, the bandwidth of the broadened spectrum is roughly similar to our measurements in air, where we measure Vavg=90.6%.

    4. DISCUSSION

    The CMPC enables post-compression for high pulse energies at compact setup length, tunable via the folding ratio Γ [Eq. (3)] while exhibiting excellent characteristics known from conventional MPCs. These include support for large compression factors, excellent beam quality and pointing properties, as well as high transmission efficiency. In principle, Γ can be increased to arbitrarily large values, enabling very compact systems. There are, however, practical limitations. The mirror size, in particular the length (corresponding to the x-dimension in Fig. 2), sets a limit on Γ [a mirror dimension estimate is provided in Eqs. (A3)–(A5) in Appendix A]. Increasing Γ also gives rise to a larger amount of total reflections. This can lead to a decreased throughput and lower efficiency of the system. The large number of reflections in a CMPC further restricts the mirror coatings to low GDD and at the same time high reflectivity. Typical coatings fulfilling these properties are quarter-wave-stack coatings with bandwidths supporting pulse durations of around 30 fs at 1030 nm. Consequently, post-compression to few-cycles cannot be achieved using the CMPC scheme. In our experiment we measure a throughput of 89% at 550 reflections corresponding to an average mirror reflectivity of >99.98%. Moreover, a large number of reflections can cause wavefront distortions, which we minimize by using mirrors with a high surface flatness of λ/20. Another factor that can restrict energy scalability is the radius-of-curvature of the mirror substrate. While mirrors with an ROC of 2 km can be manufactured [37], large ROCs can become more strongly affected by mirror distortions, such as mechanical stress or thermal effects, which could potentially limit the average power acceptance of the CMPC.

    As discussed in Section 2 the CMPC scheme requires a regime where the beam inside the cell is loosely focused in order to enable beam folding on additional mirrors. This in turn fundamentally limits the amount of nonlinear phase ϕNL(1) that can be accumulated per pass to smaller values compared to a strongly focused geometry. In gas-filled MPCs, the accumulated Gouy phase per pass ϕNL(1) sets a theoretical limit on the nonlinear phase ϕNL(1)ϕG(1) [15], which is directly related to the MPC geometry via ϕG(1)=πk/N. In standard MPCs with k/N1 this limit approaches ϕNL(1)=π, whereas in the CMPC typically k/N is lower. In our experiment, the configuration is set to N=11 and k=3 and thus ϕNL(1)0.27π. Getting close to this limit, the peak power approaches the critical power of the medium Pc and spatio-temporal couplings can become more pronounced leading to a sudden degradation of the spatio-spectral homogeneity. However, this effect is not necessarily a limiting factor for spectral broadening in a CMPC. The generally smaller ϕNL(1) in a CMPC can be compensated by simply employing more passes through the cell, while keeping the ratio k/N approximately constant. For both MPC types, operation at P<Pc ensures that excellent spatio-spectral and thus spatio-temporal pulse characteristics can be reached.

    5. OUTLOOK AND CONCLUSION

    In our proof-of-principle experiment we demonstrate post-compression of 8 mJ pulses in a compact setup. However, we can further scale up the energy by increasing the mirror sizes as well as the radius-of-curvature R of the mirrors. As an outlook for high-energy CMPC operation, we simulate the post-compression of 200 mJ, 1 ps pulses at 1030 nm central wavelength in a CMPC using 50 mbar of argon as the nonlinear medium. We show that post-compression down to 80 fs can be achieved, with a transmission of almost 80%. Figure 9 summarizes the simulation results for this scenario, indicating excellent post-compression performance. Here, we assume a mirror reflectivity of 99.98% and consider the influence of the mirror on the spectral phase of the pulse. We include the effects of all reflections on the mirrors in the simulations, which amount to 2NΓ=1170. Using focusing mirrors with R=300  m and folding mirrors of 35  cm×35  cm, a folding ratio of Γ=39 is achievable, assuming a clear aperture of five times the beam radius wm (see Appendix A Section A.2). In this case an effective length of the setup of 2.5 m can be reached. The calculated fluence on the mirrors is kept at roughly of 0.5  J/cm2, including nonlinear self-focusing effects and nonlinear mode-matching taking into account Kerr lensing in the cell [38]. The peak power is kept at P/Pc<0.6, where we still measure Vavg>90% [Fig. 8(A)], in order to avoid strong spatio-spectral inhomogeneities and prospectively still achieve a good beam quality.

    CMPC simulation for 200 mJ pulses. (A) Input pulse and simulated compressed output pulse. (B) Input spectrum and simulated broadened spectrum.

    Figure 9.CMPC simulation for 200 mJ pulses. (A) Input pulse and simulated compressed output pulse. (B) Input spectrum and simulated broadened spectrum.

    In conclusion, we introduce a novel multi-pass cell scheme enabling post-compression of high-energy laser pulses in a compact setup. The CMPC enables tuning and down-scaling of the setup size via beam folding using additional planar mirrors, using weakly focused cell modes. Instead of increasing the length of the setup, the folding ratio Γ acts as the energy scaling parameter. We demonstrate post-compression in air from 1.1 ps down to 51 fs in a CMPC with an effective length of 45 cm and a folding ratio Γ=25, while keeping the fluence comparable to a standard MPC supporting the same pulse energy but requiring around 2 m cell length. Further up-scaling options promise post-compression of pulses with an energy of 100 mJ and beyond in a table-top setup.

    Acknowledgment

    Acknowledgment. We acknowledge Deutsches Elektronen-Synchrotron DESY (Hamburg, Germany), the Helmholtz Institute Jena (Jena, Germany), and members of the Helmholtz Association HGF for support and/or the provision of experimental facilities.

    APPENDIX A

    Nonlinear Pulse Propagation Model

    The pulse propagation model used for simulations in this work is presented here. We use single-atom gases as well as ambient air, which mainly consists of molecular gases (N2 and O2), as the nonlinear media. In order to conduct simulations, we thus need to take into account time-dependent nonlinear third-order effects, which stem from coupling of the electric field to the rotational states of the gas molecules [39,40]. The full equation for the propagation model used in this work can be written in the frequency and spatial frequency domains as dE(kx,ky,ω)dz=ikzE(kx,ky,ω)+PNL(kx,ky,ω),where kz=k2(ω)kx2ky2, k(ω)=n(ω)ω/c0 with ω denoting the radial frequency, kx and ky the spatial wave-numbers, c0 the speed of light in vacuum, and n(ω) the refractive index. The second term PNL in Eq. (A1) contains all the nonlinear effects used in the model. Here we include the Kerr effect via the nonlinear refractive index n2 up to its first-order derivative, as well as the single damped-oscillator model for the molecular response as described in Refs. [40,41]. The gas-specific single damped-oscillator model is described by the damping time Γ, the frequency Λ, the Raman-Kerr nonlinear refractive index n2R, and the Raman-Kerr fraction fR. In space and time domains, the nonlinear polarization PNL can be written as PNL(x,y,t)=πε0c0n(λ)λ((n2|E|2i(λ2πn2c0+dn2dω)(2E*dEdt+EdE*dt))(1fR)+fR  n2R  Γ2/4+Λ2ΛIm(e(Γ/2iΛ)t  e(Γ/2iΛ)t|E(t)|2dt))E,where E=E(x,y,t) is the electric field of the pulse, E* the complex conjugate of E, ε0 the dielectric constant, n(λ) the refractive index, and n2 the nonlinear refractive index. The first derivative dn2/dω is determined using the scaling formula described in Ref. [42] [Eq. (A7)]. Equation (A2) describes the delayed molecular response of the medium. For the case of 1 bar of air, we use Λ=12  THz, Γ=10  THz, fR=0.6, which we extract from Ref. [40], as well as the Raman-Kerr nonlinear refractive index n2R=42×1024  m2/W [43]. We further use n2=8×1024  m2/W [44] for air and n2=24×1024  m2/W for Krypton [45]. We solve Eq. (A1) using a radially symmetric (2 + 1)D split-step approach with E=E(r,t) and the spatial coordinate r=(x2+y2), where in the spatial domain, the Fourier-transforms are replaced by Hankel-transforms. (In the simulations we use the “pyhank” package by GitHub user etfrogers for Hankel-transforms.)

    Calculation of CMPC Mirror Dimensions

    We here provide some useful equations describing the required minimum mirror dimensions for the CMPC. We consider a pulse energy E, a threshold fluence Fth, as well as the CMPC configuration with k and N. We determine the focusing mirror radius-of-curvature R, the diameter of the focusing mirrors D (which is the same as the height of the planar mirrors), and the width of the planar folding mirrors W. The corresponding dimensions are shown in Fig. 10. In order to calculate R, we use the equation for the fluence in the focus F0=4E/(λRsin(πk/N)) [32], set F0Fth, and re-arrange such that R4EλFth1sin(πk/N).To find out the width of the planar folding mirrors W, we need to ensure that the beam at any reflection has a sufficient free aperture. For this, we define a factor β, where βwm is the distance between the spot on the focusing mirror and the first reflection on the folding mirror (Fig. 10), and wm=((Rλ/π)·tan(πk/2N))1/2 [32] is the 1/e2 beam radius on the focusing mirror. With some basic geometric considerations, we arrive at WΓ4β4EπFthtan(πk/2N)sin(πk/N),defining the minimum width of the planar folding mirrors. Here, Γ is the folding ratio. We typically choose a value β=5 for our size estimations. The height of the mirrors, or equivalently, the minimal diameter of the focusing mirrors D, can be calculated via D(2N+1)βπ3/2Rλtan(πk/2N).

    CMPC geometry.

    Figure 10.CMPC geometry.

    Spectral homogeneity V(x,y) [Eq. (A6)] for each Vavg data point shown in Fig. 8 in the main text. The circles indicate the area of integration for the calculation of Vavg [Eq. (A7)].

    Figure 11.Spectral homogeneity V(x,y) [Eq. (A6)] for each Vavg data point shown in Fig. 8 in the main text. The circles indicate the area of integration for the calculation of Vavg [Eq. (A7)].

    Spectral homogeneity V(x,y) [Eq. (A6)] for spectral broadening in air at 8 mJ. The circle indicates the area of integration for the calculation of Vavg [Eq. (A7)].

    Figure 12.Spectral homogeneity V(x,y) [Eq. (A6)] for spectral broadening in air at 8 mJ. The circle indicates the area of integration for the calculation of Vavg [Eq. (A7)].

    Spectral homogeneity V(x,y) [Eq. (A6)] for the MPC comparison measurement shown in Fig. 8(C) in the main text. The circle indicates the area of integration for the calculation of Vavg [Eq. (A7)].

    Figure 13.Spectral homogeneity V(x,y) [Eq. (A6)] for the MPC comparison measurement shown in Fig. 8(C) in the main text. The circle indicates the area of integration for the calculation of Vavg [Eq. (A7)].

    Tools

    Get Citation

    Copy Citation Text

    Arthur Schönberg, Supriya Rajhans, Esmerando Escoto, Nikita Khodakovskiy, Victor Hariton, Bonaventura Farace, Kristjan Põder, Ann-Kathrin Raab, Saga Westerberg, Mekan Merdanov, Anne-Lise Viotti, Cord L. Arnold, Wim P. Leemans, Ingmar Hartl, Christoph M. Heyl, "Compact, folded multi-pass cells for energy scaling of post-compression," Photonics Res. 13, 761 (2025)

    Download Citation

    EndNote(RIS)BibTexPlain Text
    Save article for my favorites
    Paper Information

    Category: Ultrafast Optics

    Received: Sep. 3, 2024

    Accepted: Jan. 8, 2025

    Published Online: Mar. 3, 2025

    The Author Email: Arthur Schönberg (arthur.schoenberg@desy.de)

    DOI:10.1364/PRJ.540977

    Topics