Introduction

Laser-plasma hard x-ray sources are a suitable alternative to conventional accelerator-based x-ray sources of high brightness1,2 and for applications requiring a high x-ray flux and a micrometric size such as x-ray phase contrast imaging3,4. Their characteristics and operating capabilities have been improved for decades5,6,7,8 especially thanks to the continuous progress of femtosecond lasers in terms of intensity, spatial and temporal beam quality, repetition rate, reliability and compactness. In particular, such x-ray sources of ultrashort pulse duration driven by femtosecond laser systems offer an intrinsic synchronization for pump and probe experiments like time-resolved x-ray diffraction9 and time-resolved x-ray absorption spectroscopy10.

Laser plasma hard Kα x-ray sources are generated by the interaction of a high intensity laser pulse (I ≥ 1016 W/cm2) with most generally a high Z solid target. A spectrum composed of a large x-ray Bremsstrahlung emission dominated by a spectrally narrow Kα line is obtained coming from the interaction of hot electrons accelerated in the laser-induced plasma with inner shell electrons.

In first approximation, the laser plasma Kα x-ray source size S at full width at half maximum (FWHM) can be directly related to the size of the FWHM laser focal spot diameter d. Their ratio, S/d, approaches a value reported at best close to one11. This minimum ratio is generally obtained for a rather low laser intensity on target and a high value of the laser temporal contrast ratio (CR). The latter is here classically defined as the main pulse peak intensity (fs) over the background Amplified Spontaneous Emission (ASE, ns) intensity, such as CR = Ipeak/IASE. Many studies have shown that the source size ratio S/d increases with the laser intensity4,12,13,14. The smallest Kα x-ray source size produced by laser plasma reported today is ~ 6 µm with a ratio × 4.6, using a massive silicon target and a p-polarized laser beam of 30 fs duration focused with a nominal incidence angle of 45° and with an intensity up to 5 × 1017 W/cm2 at λ = 800 nm15. By applying cleaning pedestal techniques or second harmonic conversion of the laser, studies have shown that an improvement of the contrast ratio leads to a pronounced decrease of the x-ray source size13,16,17,18. However, at shorter laser wavelength (λ ~ 400 nm), the reduced kinetic energy acquired by the hot electrons is less favorable for x-ray generation since the laser field ponderomotive potential scaled as Iλ2 19. Among these works only very rare studies have been published concerning the impact of different values of the temporal contrast ratio, at the fundamental wavelength of 800 nm, on the evolution of the x-ray source size for a large laser intensity range13.

The mechanisms of x-ray source size enlargement compared to the focal spot size in relation with the driving laser parameters are still not well understood. These mechanisms rely on the transport and energy of the hot electrons which are very complex to study both experimentally and numerically. The principal explanations provided in the literature are the divergence angle of propagation of hot electrons depending mainly on the laser intensity and the laser incidence angle with the target surface, and their larger stopping path inside the target when the laser intensity increases20,21. Other formulated hypotheses argue the wings of the laser focal spot21,22, and/or the surface electromagnetic fields effects23,24 which induce a lateral propagation of hot electrons along the target surface25,26. Besides, space-charge inhibition of hot-electron transport can reduce dramatically their penetration in the target27,28. These collective effects become increasingly important for laser intensities I > 1018 W/cm2 29. It is therefore instructive to study the penetration depth of hot electrons and their trajectories (divergence angle) inside the solid. Hot electron energy and their angular distribution can be measured by imaging electron energy-sensitive plate stacks or an electron spectrometer. They can also be measured by indirect techniques such as coherent transition radiation30 at the rear target surface when using thin foils, or Kα fluorescence measurement from multi-layer targets composed of materials of different atomic numbers Z31,32,33. However, the effect of varying the thickness of a single Z foil target on the x-ray source size and Kα photon number has hardly ever been reported.

The main objective of the present work is to bring original experimental measurements of the size of a Kα x-ray source as a function of the two main laser driving parameters, e.g. the laser intensity and its temporal contrast ratio. Compared to the literature, we intend to go one step further by expanding such parametric study to a wide number of temporal contrast values (4 CR, from 6.7 × 107 < CR < 3.3 × 1010) and a larger range of laser intensity (I ~ 1017–2.8 × 1019 W/cm2) with the same experimental setup, and also by providing dedicated experiments to progress in knowledge and understanding of the characteristics of a Kα x-ray source driven by high intensity ultrashort lasers. So, the correlation between the x-ray source size and the hot electron path inside the solid target is investigated by varying the thickness of a molybdenum foil as target down to 4 µm and the laser focal spot size is varied for two fixed laser energies and one fixed laser intensity to study its impact on the x-ray source size for the best CR. Crossed analysis of obtained results with the help of literature confirm the prevalence of surface electromagnetic fields among others mechanisms for explaining the large size broadening of the x-ray source generated at high laser driving intensity even for temporal contrast as high of 3.3 × 1010. Finally, this large and detailed experimental study of the size broadening of the x-ray source as a function of laser driving characteristics can be used for theoretical future studies and simulations, especially to gain insight into electron transport and for the perspective of optimizing the brilliance of such hard x-ray source induced by laser-driven plasmas at very high laser intensity.

Experimental setup

The 10 Hz—20 TW—800 nm laser beamline of ASUR facility is used for the study34. The p-polarized, \(\sim\)25 fs laser pulse is focused at an incidence angle of 45° with a f/4.5 off-axis parabolic mirror on a molybdenum (Mo, Z = 42, E = 17.48 keV) solid disk or foils of 100 mm diameter and of variable thickness. Note that a fresh target surface is provided for each laser shot. At 45° incidence angle, laser absorption through vacuum heating mechanism is optimized when the pre-plasma density scale length L is small compared to the laser wavelength λ35. The laser focal spot diameter is dx × dy = 4.9 × 4.6 µm2 at FWHM. The energy is varied from ~ 1 mJ up to ~ 254 mJ on target with ~ 40% of the encircled energy at FWHM. By considering the real 2D beam profile, laser peak intensity on target surface vary from ~ 1017 W/cm2 up to ~ 2.8 × 1019 W/cm2. The temporal contrast ratio CR is modified by using up to three saturable absorbers. It varies from 6.7 × 107 to 3.3 × 1010 for a measured ASE level at t = −480 ps before the peak of the main pulse (see “Methods”, Fig. 6).

The FWHM effective size of the x-ray source is determined by the knife edge technique assuming a gaussian shape12. Due to the geometrical configuration of the targetry system which ensures a fresh surface for each laser shot12, all the defects related to the geometric imperfections of the target and of its displacement system induce a more pronounced widening of the x-ray effective source size in the horizontal dimension compared to the vertical one. This is particularly significant when multipulse laser irradiation (~ 20 to ~ 1000 laser shots) is required to acquire an image of the knife edge. Thus, we only measure the vertical diameter of the x-ray source since the objective is to link the enlargement of the x-ray source size with physics phenomena and not to mechanical displacement and/or geometrical surface defects. The target characteristics are reported in “Methods”. Finally, the absolute number of Kα photons per pulse is deduced from reconstructed x-ray spectrum obtained using the photon counting method36 in which a x-ray CCD camera (PIXIS-XB from Princeton Instruments, 1024 × 1024 pixels of 13 × 13 µm2 pixel size) acts as a dispersive spectrometer.

Results and discussion

Impact of the temporal contrast and laser intensity

The vertical size of the x-ray source and the Kα photon number per steradian and per shot, N, are reported in Fig. 1a,b respectively, for four CR values (3.3 × 1010, 4.2 × 109, 1.4 × 109, and 6.7 × 107), and as a function of laser peak intensity. The peak intensity is varied by changing the laser pulse energy while the duration and the laser spot size on target are kept unchanged.

Figure 1
figure 1

(a) FWHM x-ray vertical diameter versus laser intensity on target, for four temporal contrast ratio values. The error bars correspond to the standard deviation of the average FWHM line spread function measured for three different zones of the knife edge image. The dashed lines for the x-ray source diameter data correspond to a Iε power fit with ε varying with the contrast ratio value. The grey horizontal dashed line represents the FWHM vertical laser focal diameter (~ 4.6 µm). (b) Kα photon number versus laser intensity on target, for four temporal contrast ratio values. The error bars are given by the standard deviation calculated from three independent measurements of NKα.

In Fig. 1a, we first observe that the evolution of the x-ray source size, S, as a function of the laser intensity on target can be fitted using a power law with S \(\propto\) Iε. This behavior has already been put forward in our previous works4,12, but only for one temporal contrast ratio value (CR ~ 1010). Here, we show that it can be extended to all the temporal contrast values explored in the present study and that the ε factor decreases from 0.34 to 0.27 when CR is deteriorated. We also observe that the higher is the temporal contrast value the smaller is the source size for the whole laser intensity range. Secondly, at relatively low laser intensity (~ 1017 W/cm2 to 3.0 × 1017 W/cm2), the smallest x-ray source size, S ~ 5 µm, is obtained for the highest temporal contrast ratio (3.3 × 1010) corresponding to only × 1.1 the laser focal spot size. It is the smallest Mo-based Kα laser plasma x-ray source size reported up to now at the fundamental Ti:Sa wavelength (λ = 800 nm). This is due to the high contrast of the laser driver in addition to the large numerical aperture of the focusing parabola. Inversely, for the lowest CR (6.7 × 107), the x-ray source size reaches 18 µm. This spatial enlargement is correlated to the degradation of the contrast ratio since the measurement of the laser spatial profile for the four CR show no noticeable difference (see “Methods”, Fig. 7). When the laser intensity increases, the source size increases up to ~ 76 µm for I = 2.8 × 1019 W/cm2 and the lowest CR. However, at very high laser intensities, the difference between the source sizes tends to reduce for the different temporal contrasts, as shown by the decrease of the parameter ε for decreasing ASE contrast value (Fig. 1a). Indeed, for I < 3.0 × 1017 W/cm2, SCR = 6.7×107/SCR = 3.3×1010 ~ 4, while this ratio is equal to ~ 2.5 for I > 1.0 × 1019 W/cm2.

Conversely, in Fig. 1b, we observe that the number of produced Kα photons, NKα, is independent of the CR value for lower (I ≤ 1.0 × 1017 W/cm2) and higher (I > 3.0 × 1018 W/cm2) laser intensities, and slightly CR-dependent at intermediate laser intensities where we notice a difference with a lower Kα photon number produced when CR deteriorates. In a previous work we already observed a similar behavior and explained it by the competition and transition between different laser absorption mechanisms governing the generation of hot electrons like resonant absorption, vacuum heating or J × B heating34. Their dominance depends on the pre-plasma density scale length L induced by the ASE and pedestal at the front of the solid target and on the laser intensity35,37.

In support of the experiment, the low fluence Lagrangian one-dimensional hydrodynamic code ESTHER38 is used to simulate the plasma induced by both ASE and laser rising edge up to 1.2 ps before the peak intensity. The ESTHER code is currently used for such applications in this intensity regime39,40,41. This allows us to investigate the influence of the pre-plasma length, L, and therefore of the laser absorption mechanisms on the x-ray source size and the Kα photon number. The simulation considers the temporal evolution of the laser intensity as measured with a third order autocorrelator (see “Methods”). L is defined using the simulated electron density profile ne normalized by the critical density nc (~ 1.72 × 1021 cm−3 at λ = 800 nm) as a function of laser intensity, from 3.0 × 1016 W/cm2 to 1.0 × 1019 W/cm2 and for the four temporal contrast ratios experimentally studied. It corresponds to a range of maximum ASE and pedestal intensities of ~ 3.0 × 1011 W/cm2 to 1.0 × 1014 W/cm2 (see “Methods”, Fig. 6). An example of electron density profile is given in Fig. 8 in “Methods”. Figure 2 gives the different values of L/λ at −1.2 ps before the laser intensity peak for the four experimental temporal contrasts and five laser intensities.

Figure 2
figure 2

Simulated pre-plasma density scale length L normalized by λ for the four experimental CR studied and from I = 3.0 \(\times\) 1016 W/cm2 to 1.0 \(\times\) 1019 W/cm2. At relativistic laser intensity regime (I > 2.2 × 1018 W/cm2), the Lorentz factor γ must be considered in the calculation of the critical density nc, giving a measured L at I = 1019 W/cm2 lower than the one at I = 1018 W/cm2 for CR = 6.7 × 107. The grey horizontal dashed line represents L/λ = 0.1.

We observe that the density scale length L/λ varies almost by three decades from ~ 0.02 at the lowest intensities and CR ≥ 1.4 × 109 up to 10 at the highest intensities and CR < 3.3 × 1010. As described by Gibbon et al.35,42, we deduce in a first approximation that at low laser intensities, the dominant absorption mechanism for at least 3 CR values is vacuum heating (L/λ < < 0.1) which progressively changes to resonant absorption (L/λ > 1) for higher intensities35. However, the light pressure PL, which cannot be included in the ESTHER simulations, becomes to be significant at high laser intensity and has to be taken in consideration. Indeed, it progressively crushes the plasma, thus providing conditions (L/λ < < 0.1) for which the absorption of the laser energy through vacuum heating and J × B heating becomes highly dominant35,37. The laser intensity I0, for which PL = Pth, can be analytically estimated as I0 = \(\frac{{\rho {\rm c}}}{2{{\rm cos}\theta }}\)(9.45 × 103)2 \({\mathrm{I}}_{\rm{ASE}}^{1/2}\) (see “Methods”). When the laser intensity I exceeds I0, the radiative pressure pushes the pre-plasma back to the target, resulting in a steepening of the electronic density gradient and consequently the mechanisms of laser absorption are modified43,44. We evaluate this compression at ~ 500 fs before the time corresponding to the maximum peak intensity (see “Methods”, Fig. 6). For I > 3 × 1018 W/cm2 and for all CR, the calculation gives Ld/λ < < 0.1.

In brief, for all CR the same absorption mechanisms are dominant at least for the lowest and highest laser intensities. This analysis explains why the number of Kα photons observed on Fig. 1b is independent of CR at low and high intensities since heating electron mechanisms are the same, as demonstrated. However, the x-ray source size which depends on the trajectories of hot electrons strongly varies with CR values (see Fig. 1a) even for the same laser intensity and as shown for the same absorption mechanism. This preliminary conclusion motivates further dedicated experiments presented hereafter in order to investigate such spatial broadening phenomena of the x-ray source.

Impact of the laser beam spatial contrast

In this section, the increase of the x-ray source size is studied in correlation with the growing importance of the laser beam wing size at very high laser intensity. For this purpose, the spatial distribution of the focal spot intensity is characterized on a large dynamic range (> 104) by removing calibrated neutral density filters in front of the beam analyzer detector (see “Methods”), as shown in Fig. 3. Note that the transmission of a single neutral density filter is ~ 16.7% at λ = 800 nm.

Figure 3
figure 3

(a) Snapshots of the laser spatial beam distribution. Differences between each snapshot correspond to the removal of one neutral density filter. (b) Stacked vertical laser focal spot size profile coming from the red line profile on the images presented in (a). Wings are present and their intensities are I < 5.0 × 1016 W/cm2, even when the laser peak intensity is measured at I ~ 1.0 × 1019 W/cm2.

As show in Fig. 3, low intensity wings are spread over 100 µm from the laser focal spot center. In fact, for a very high laser peak intensity I = 1.0 × 1019 W/cm2, the intensity of these wings is inferior to 5.0 × 1016 W/cm2. If we make the simplified assumption of the locality of the x-ray emission in concordance of the local laser intensity value the Kα photon number produced in the zone excited by the wings of the laser is at least ~ 104 lower than in the zone corresponding to the laser peak intensity (I = 1.0 × 1019 W/cm2) as we can extrapolate from Fig. 1b. Therefore, these photons are not even experimentally detected, and we consider that the contribution of these wings in the x-ray emission can be neglected even at very high laser intensity. These observations rule out this hypothesis based on the spatial beam contrast which was discussed in literature21,22 to explain the enlargement of the x-ray source in our operating conditions.

Impact of the foil target thickness

Hot electron scattering in a solid is dependent on their temperature Th and their energy Eh and therefore on the laser intensity35,45,46. Since the generated hot electron beam presents a divergence angle, an enlargement of the x-ray source size is thus expected at high laser intensities due to their larger penetration depth46,47,48. As a result, experiments where the target thickness is varied would serve to investigate the effects of the hot electron trajectories inside the target on the x-ray source size broadening. More precisely, it would be expected that reducing the thickness of the target to micrometric values should lead to a reduction in the x-ray source size and of the number of Kα photons produced in comparison with a massive target, at least for high laser intensities.

In order to choose the thickness of the target foil, the depth for which a produced Kα photon will not escape from the Mo target due to its absorption in the solid is determined. Practically, one estimates the depth for which the Kα photon number is divided by the factor e. Using the mass attenuation coefficient for a Mo Kα photon propagating in a Mo target (µ/ρ  20 cm2/g yielding to µ  203.4 cm−1)49, the absorption length in Mo is  50 µm. The x-ray source size and the Kα photon number are consequently studied for three molybdenum target thicknesses below this length (25 µm, 10 µm and 4 µm) where re-absorption phenomena do not significantly occur, and compared with the previous results obtained with the 6 mm target presented in Fig. 1. Results are reported in Fig. 4 for two contrast ratio values of 3.3 × 1010 and 1.4 × 109.

Figure 4
figure 4

Comparison of the x-ray vertical source size and the Kα photon number versus laser intensity between a 6 mm massive Mo target and Mo foil targets of 25 µm, 10 µm and 4 µm thickness. For CR = 3.3 \(\times\) 1010 in (a) and for CR = 1.4 × 109 in (b).

In Fig. 4a, we notice that the x-ray source size and NKα evolution are independent of the target thickness for both CR. For the x-ray source size, a maximum variation of × 1.2 between the different target thicknesses is observed in most cases. Thus, these experimental results suggest that the hot electron penetration depth is less than 4 µm. Note also that for the foil targets, we further checked if these results could not be explained by the refluxing effect at the target rear surface8. For inhibiting this effect, we set a SiO2 substrate placed behind and in close contact with the Mo foils. Results of the evolution of S and N are again identical with and without the substrate ruling out this hypothesis.

The electron penetration range is generally limited both by collisional and collective effects. Collisional stopping power is related to the electron beam collisions with the propagation medium material while the collective effects come from the existence of self-generated return currents provided by the hot electrons31,50. Based on Bell et al. studies24, Volpe et al. have shown for an aluminum target that the collective effects start to compete with collisional effects for laser intensities ~ 1017 W/cm2 and become dominant, over collisional processes with the ions and electrons of the target medium, for intensities ~ 1018 W/cm2 28. Here we use the analytical description proposed by Volpe et al. to calculate the total electron penetration depth, Rtotal, taking into account collisional and collective effects from the harmonic mean, \({\mathrm{R}}_{\rm{total}}= \frac{{\rm{R}}_{\rm{collisional}}{\mathrm{R}}_{\rm{collective}}}{{\mathrm{R}}_{\rm{collisional}}+{\mathrm{R}}_{\rm{collective}}}\) (see “Methods”). For I = 1.0 × 1017 W/cm2 and 2.0 × 1019 W/cm2, the maximum electron penetration depth Rtotal is 0.13 µm and 2 µm respectively, while Rcollisional would be ~ 0.35 µm and 7.6 µm without considering the collective effects. Rtotal is lower than the thinnest target thickness used in the present experiment (4 µm). This electron penetration depth limitation induced by collective effects may therefore explain why in our experimental condition N is independent of the target thickness. However, this phenomenon cannot explain the observed enlargement of the x-ray source size.

Reich et al. have observed similar experimental results, consisting of a large broadening of the Kα source size compared to the laser focal spot size. Their explanation is the existence of strong self-induced fields which redirect the hot electrons over the target surface21. Later for experimental conditions very similar to ours, Li et al. observed fast electrons emitted along the surface (SFE) with a relatively large ratio between the number of SFE to the total number of hot electrons (17% to 28%)51. Their 2D-PIC numerical simulations reveal that they are induced by surface quasistatic electromagnetic fields which confine them at the target surface. In addition, it has been shown that these SFE are dependent of laser intensity and temporal contrast52. In view of these observations and of the analyses of the evolution of the x-ray source size as a function of the target thickness, we thus suggest that the trajectory of hot electrons propagating through the solid target might be bent by magnetic fields in the interaction region and the overdense plasma.

Impact of the laser focal spot size

The divergence of the hot electron beam formed by the laser-induced plasma is in part due to the variation of the laser intensity along the laser focal spot profile53. Considering also that the surface magnetic field is dependent on the focal spot size54, it is instructive to evaluate the influence of the laser focal spot size on the x-ray source size. In this perspective, the laser spot size is varied by moving the parabola along the focal axis toward the target or away from it, corresponding to a defocus range of ± 100 µm. This is superior to the experimental Rayleigh length zR = 67 µm when the laser beam is approximated as gaussian. However, when the beam is not at the best focus (in intermediate field), its spatial profile can no longer be considered as gaussian. Nevertheless, the beam diameter can still be measured, and the laser spot diameter is varied between ~ 5 and ~ 17 µm in our experimental conditions. This experiment is performed for the highest temporal contrast ratio (CR = 3.3 × 1010) and for three laser conditions: two laser energies fixed at EL = 40 mJ and EL = 180 mJ respectively, and a constant maximum peak intensity set at Ifixed = 4.0 × 1018 W/cm2 by compensating the increase of the focal spot size with the increase of the laser energy, considering the real 2D laser beam profile. Figure 5 shows the evolution of the ratio between the size of the x-ray source S and the size of the laser focal spot d as a function of d.

Figure 5
figure 5

Ratio between the x-ray source size (S) and the laser focal spot size (d) versus laser focal spot size for two constant laser energies (40 mJ—black square and 180 mJ—red point) and for a fixed maximum laser intensity on target (4.0 × 1018 W/cm2—blue triangle).

In Fig. 5, we show that S/d increases by a factor of ~ 2 when the laser energy increases by 4.5 whatever the size of the laser spot. Moreover, the smaller the focal spot size the higher is S/d, and finally when the focal spot size increases, S/d tends towards one. These experimental observations are in accordance with a scenario where the self-induced electromagnetic fields increase with laser intensity and are higher with tighter focal spots54. These fields will have the tendency to confine the hot electrons within the target along the surface as discussed above and thus to increase the size of the x-ray source.

Finally, our result presented in Fig. 5 may explain the low ratio S = 1.7d with larger laser spot size d = 15 µm, obtained by Fourmaux et al. for similar laser intensity (I = 6.0 × 1018 W/cm2)15. In this work, it was only attributed to the improvement of the CR. However, thanks to the Fig. 5, we show that low S/d ratios found in the literature can also be explained by the laser focal spot size itself with reduced effects of self-induced electromagnetic fields on hot electron trajectories for larger laser focal spots compared to smaller ones.

Conclusion

In this paper, the x-ray source size enlargement is experimentally studied by varying the laser intensity on target, the temporal contrast of the laser pulse, the target thickness and the laser focal spot size. We observe that the x-ray source size is increasing with the laser intensity on target and the degradation of the temporal contrast. Since the ASE and pedestal intensity in a time scale of ns to ps is crucial in tailoring the plasma properties and electron heating, a 1D numerical simulation (ESTHER) is used to estimate the electron density scale length L/λ as a function of laser intensity and CR. At non-relativistic laser intensity regime, L/λ is extremely small (L/λ < < 0.1) except for the lowest CR = 6.7 × 107. As the laser intensity increases, the laser pressure radiation limits the expansion of the electron density scale length and even leads to a gradient steepening in such a way that L/λ becomes very small at high laser peak intensity for all CR. As already demonstrated in previous works, the parameter L/λ is of prime importance in high intensity laser-plasma interaction because it allows to explain with confidence and in relation with the different laser absorption mechanisms the evolution of the number of Kα photons produced as a function of laser intensity and temporal contrast ratio. However, this parameter is not relevant to unambiguously describe the evolution of the x-ray source size because the transport of hot electrons also depends for a fixed laser intensity on collisions, collective effects and self-generated electromagnetic fields. In an original study, where the thickness of the Mo target is varied down to 4 µm and for a large range of laser intensity (~ 1017 W/cm2 to 2.0 × 1019 W/cm2), we show that both the evolution of N and the x-ray source size S are outright independent of the target thickness (from 6 mm down to 4 µm), whatever the temporal contrast studied. The limited penetration depth (< 4 µm) of hot electrons experimentally observed may be explained by considering the collective effects of the hot electrons scattered inside the solid target. In a last experiment, the laser focal spot size d is varied to study its influence on the x-ray source size enlargement. The reduction of the laser focus spot size increases the ratio between the size of the x-ray source and the size of the laser focus spot. The cross-interpretation of all our experimental results suggests that the enlargement of the size of the x-ray source is strongly linked to the presence of self-generated electromagnetic fields which tend to bent a significant part of the hot electrons generated by J × B heating or vacuum heating along the surface of the target.

Considering the optimization of the peak brilliance of such Kα x-ray source in relation with source size minimization, we show that trying to limit the x-ray source size enlargement by using target foils as thin as 4 µm does not bring any improvement in such experimental conditions even with a high temporal contrast. Working at high laser intensity allows to increase dramatically the number of produced x-ray photons. Additionally, working with a very high temporal and spatial contrast of the laser beam allows to limit the x-ray source size broadening. The transport of the hot electron to a much larger zone than the laser focal aera induced by the self-electromagnetic fields is intrinsic of such laser interaction with a solid target but this is limited with high temporal contrast ratio. The improvement of temporal contrast ratio to higher value than ~ 3.0 × 1010 by various techniques55,56 as already demonstrated may be relevant to still enhance the peak brilliance of the x-ray source.

Moreover, this wide set of experimental data on the evolution of the x-ray source size as a function of a large range of parameters (temporal contrast ratio, laser intensity, target thickness and laser focus beam size) with the same experimental setup can be useful for the development and validation of numerical studies. This includes a combination of PIC and Monte-Carlo approaches, which are far out of the scope of the present paper, to study in details the energy distribution and trajectories of hot electrons which remains still today not well understood. Indeed, when using such tightly laser focus and a laser pulse including ASE and pre-pulses at different time scales from ns to fs, 2D to 3D simulations are required. As an example, such simulation studies will be useful to try to explain the dependence of the x-ray source size to the temporal contrast ratio, experimentally shown as S  Iε, where ε varies with the temporal contrast ratio.

Methods

Targets

The 6 mm massive Mo target thickness comes from Kurt J. Lesker Company. Its characteristics are: 99.95% purity, rugosity < 1 µm, and a flatness of ± 20 µm overall the 100 mm disk diameter. The 25 µm, 10 µm and 4 µm Mo targets come from GoodFellow with thickness variation of ± 25% below 10 µm, and ± 15% for 10 µm and 25 µm thickness. The molybdenum foils are held only in their center and on a few millimeters in the periphery of the outer diameter of the disk by a rotating support partially hollowed out. This latter also ensures the obtention of vacuum at the rear side of the targets. The deviation to a perfect flatness of the thinnest target (4 µm) is measured to be < 40 µm, which is lower than the Rayleigh range of the laser (67 µm), ensuring approximatively the same laser intensity on target between two consecutive shots. More details about the target positioning can be found elsewhere12.

Measurement and modification of the temporal contrast ratio

The ASE and pre-pulses are measured up to 480 ps before the main pulse by a 3rd order autocorrelator (Sequoia, from Amplitude Technologies, see Fig. 6). Note that three pre-pulses are also present few ns ahead of the main pulse. However, they do not influence the x-ray generation as their contrast value is higher than 108 34. A pre-pulse around −12 ps is a measurement artefact coming from optics inside the Sequoia. The ASE duration is also measured up to 3 ns ahead of the main pulse thanks to the implementation of a ~ 1 m delay line added to the Sequoia autocorrelator. One SA (always present) is positioned between the booster amplifier (a multipass amplifier located after the oscillator) and the stretcher, to seed the regenerative amplifier with a clean and energetic pulse. Two extra saturable absorbers (SA, RG-850 long pass filter) are inserted in the laser chain between the different preamplifier stages located before the compressor in order to decrease the ASE level and obtain different contrast ratio values for the present study (CR = Ipeak/IASE). Four experimental configurations can be achieved by inserting or removing these SA:

  • Two extra SA in the laser chain corresponding to the black curve of laser intensity time evolution in Fig. 6, CR = 3.3 × 1010.

  • Only one extra SA in the laser chain, which increases the ASE level (red curve in Fig. 6) and so decreases the CR down to 4.2 × 109.

  • Zero SA between the preamplifiers. The ASE is again increased (blue curve in Fig. 6). In this case, CR = 1.4 × 109.

  • The same configuration as for CR = 1.4 × 109 (only the SA located in the Booster amplifier is present), but the pump energy in the Booster amplifier is reduced. The regenerative amplifier is thus seeded by a train of pulses with moderate energy and the pumping level of the regenerative amplifier is increased to correctly saturate the amplifier. The ASE level is therefore increased, giving CR = 6.7 × 107 (green curve in Fig. 6).

Figure 6
figure 6

Temporal characterization of the laser pulses interacting with the Mo target. The main pulse corresponds to the time t = 0 ps and the ASE measurement goes up to −480 ps before the main pulse as shown in (a). The CR is defined by Ipeak/IASE. (b) Zoom up to −50 ps before the main pulse, corresponding to the grey dashed rectangular line in (a).

Characteristics of the laser focal spot

The laser focal spot is imaged by a relay imaging system coupled to a beam analyzer with a magnification of 8.4. The detector is a 14 bits CCD device from Gentec-EO with 1024 × 1024 pixels of 4.65 × 4.65 µm2 surface. Two parallel planar optical uncoated substrates at 45° incident angle are placed after the beam focus to attenuate the reflected beam for spatial characterization, even at full laser energy. A combination of five neutral density (ND) filters NE10A-B from Thorlabs placed in front of the CCD detector allows to precisely characterize the focal spot up to the maximum pulse energy of 254 mJ. More details about the beam imagery system can be found elsewhere57. We also thoroughly characterized the focal spot profile for the different contrast configurations used in this study. These measurements are presented in Fig. 7, here for the vertical profile and for I ~ 1.0 × 1017 W/cm2.

Figure 7
figure 7

Spatial vertical profile of the laser focal spot for the four temporal contrast values studied. Profiles are CR independent.

As shown in Fig. 7, the spatial profile is independent of the temporal contrast configuration. To precisely characterize the laser focal spot wings, some ND filters are removed (up to 3), in order to increase the detection dynamic. The maximal pixel intensity is then saturated revealing lower intensity levels coming from the wings, as we can see in Fig. 3. The profile obtained in this figure is obtained by considering the respective attenuation coefficients of the calibrated neutral density filters that were removed.

Electron density gradient simulations

The ESTHER code is used to simulate the electron density gradients obtained after the interaction of the ASE and pulse pedestal with the molybdenum target38. The experimentally characterized ASE distribution was used in the numerical simulation up to −1.2 ps before the main pulse, which corresponds to the time before the rising edge of the main pulse. Note that the rising time is the same for the four CR. Optical indices for molybdenum are given by Palik in the solid state58 and by Faussurier et al. in the plasma state59. The equations of state (EOS) used come from SESAME, a database developed by the Los Alamos National Laboratory60. We also confronted these EOS with those from work of Bushman, Lomonosov and Fortov (BLF)61. No significant differences were observed between them concerning the generation of pre-plasma on a molybdenum target. We simulated electron density gradients for the four CR studied in the present paper, and for five laser intensities (3.0 × 1016–1.0 × 1017–3.0 × 1017–1.0 × 1018–1.0 × 1019 W/cm2). As an example, Fig. 8 shows the raw electron density extracted from ESTHER, normalized by the critical density nc for CR = 6.7 × 107. The pre-plasma length L was determined when ne/nc = 0.5γ/e, with γ = \(\sqrt{1+{\mathrm{a}}_{0}^{2}}\). The dimensionless intensity parameter a0 is defined as a0 = \(\frac{{\mathrm{eE}}_{\rm{L}}}{\mathrm{m}{\upomega }_{\mathrm{L}}\mathrm{c}}\) with e and m respectively the electric charge and the electron mass, c the speed of light, ωL the laser frequency and EL the laser electric field. Note that a0 > 1 represents the relativistic laser intensity regime.

Figure 8
figure 8

Normalized electron density ne/nc as a function of the position from the target surface for CR = 6.7 \(\times\) 107 and for the five laser intensities simulated.

Estimation of the laser intensity I0 and of the steepening of the normalized electron density gradient L/λ

The laser intensity I0, for which PL = Pth, can be analytically estimated. The plasma pressure is defined as Pth = \(\uprho \frac{{\mathrm{v}}_{\mathrm{s}}^2}{2}\) where vs is the ionic sound speed. The latter is taken as vs = 1.37 × 106Z1/8ϒ-9/16(Iλτ1/2)1/462 with Z the target atomic number, τ the laser pulse duration and ϒ = \(\frac{Z}{2{[{Z}^{*2}\left({Z}^{*}+1\right)]}^{1/3}}\) where Z* is the effective charge that is determined using the code FLYCHK63. For our experimental conditions (Z = 42, λ = 0.8 µm, τ ~ 25 fs and Z* = 6.5), vs = 9.45 × 103 \({\mathrm{I}}_{\rm{ASE}}^{1/4}\) with IASE the intensity of the ASE. Thus, I0 = \(\frac{{\rho {\rm c}}}{2{{\rm cos}\theta }}\)(9.45 × 103)2 \({\mathrm{I}}_{\rm{ASE}}^{1/2}\) with c the speed of light, \(\uptheta =45{ \%}\), and ρ = meγnccos2θ the mass density, with me the electron mass, nc the critical electron density, and γ the Lorentz factor.

When the laser intensity I exceeds I0, the radiative pressure pushes the pre-plasma back to the target. It is important to evaluate this phenomenon at a crucial time of the interaction corresponding just before the arrival of the main pulse and of the deposition of the laser incident energy into the target. Arbitrarily, we take the time t0 equal to ~500 fs before the time corresponding to the maximum peak intensity (see “Methods”, Fig. 6). Indeed, in the temporal interval between t0 and the peak of the main pulse, plasma expansion into vacuum related to the rising edge of the main pulse is very limited. For operating conditions corresponding to ours, an ion sound speed of vs = 37 nm/ps was experimentally determined for a pre-pulse developing in the 0–4 ps temporal range before the peak of the main pulse64, yielding an estimate of ΔL/λ  0.02 during 500 fs. This justifies neglecting the influence of any further pre-plasma developing in that temporal range. Note that other choice of t0 is not critical and would yield to very similar results and conclusions. Then, an intensity averaged between I0 and up to 0.01Imax on the rising edge of the main pulse is taken to estimate the squashing of the pre-plasma induced by the laser pulse at a given laser peak intensity I. This choice results in a compromise between a minimal effect related to the application of the laser optical pressure PL at I0 and an intensity corresponding to a time ~ 500 fs (or slightly shorter) before the maximum peak intensity (see “Methods”, Fig. 6). We can thus estimate how far the laser pressure crushes the density gradient just before the main pulse thanks to the relation d = vL × t, with vL the velocity associated to the laser pressure and t the rising edge time of the main pulse taken from I0 and up to 0.01Imax. Then, this displacement d is subtracted to the simulated values of L presented in Fig. 2, such as Ld = L − d.

Calculation of the total penetration depth Rtotal

\({\mathrm{R}}_{\rm{total}}= \frac{{\mathrm{R}}_{\rm{collisional}}{\mathrm{R}}_{\rm{collective}}}{{\mathrm{R}}_{\rm{collisional}}+{\mathrm{R}}_{\rm{collective}}}\) with Rcollisonal= \(0.4\times {\mathrm{E}}^{1.265-0.0954\mathrm{ln}(\mathrm{E})}/\uprho\)47 and Rcollective = A(Th) × z0. A(Th) is a factor depending to the Maxwell-Jüttner distribution and z0 = 3.0 × 1016 \(\frac{\upsigma {\mathrm{T}}_{\rm{h}}^{2}}{{\upeta }_{\rm{L}\to \mathrm{e}}{\mathrm{I}}_{18}}\), with \({\upeta }_{\rm{L}\to \mathrm{e}}\approx\) 3.8 × 10–2 \({\mathrm{I}}_{18}^{3/4}\) the laser to electron conversion efficiency16 and σ = 9.7 × 103 \(\frac{{\mathrm{T}}_{\rm{e}}^{3/2}}{{\mathrm{Z}}^{*}\mathrm{ln}(\Lambda )}\) the Spitzer conductivity with ln(Λ) the coulomb logarithm, depending on Te. Th is determined according to Th ≈ 7 × (I16[W/cm2]\({\uplambda }_{\upmu {\rm{m}}}^{2}\))1/3 keV since vacuum heating is the dominant absorption mechanism42. The temperature of the plasma layer is calculated using \({\mathrm{T}}_{\rm{e}}= 0.53{\times \mathrm{Z}}^{*1/6}{\mathrm{I}}_{18}^{1/3}{\uplambda }^{-1/6}{\uptau }^{1/6}\mathrm{ keV}\)65.