
==== Front
Nat Commun
Nat Commun
Nature Communications
2041-1723
Nature Publishing Group UK London

39237494
51084
10.1038/s41467-024-51084-4
Article
Spatiotemporal dynamics of fast electron heating in solid-density matter via XFEL
http://orcid.org/0000-0002-7972-9894
Sawada H. hsawada@unr.edu

1
http://orcid.org/0000-0003-0638-7640
Yabuuchi T. 23
Higashi N. 4
Iwasaki T. 4
Kawasaki K. 4
Maeda Y. 4
Izumi T. 4
Nakagawa Y. 4
http://orcid.org/0000-0002-3978-8427
Shigemori K. 4
Sakawa Y. 4
http://orcid.org/0000-0001-8756-181X
Curry C. B. 56
http://orcid.org/0000-0001-6879-0422
Frost M. 5
http://orcid.org/0000-0002-2823-1842
Iwata N. 4
http://orcid.org/0000-0002-4838-976X
Ogitsu T. 7
Sueda K. 3
http://orcid.org/0000-0001-7615-7213
Togashi T. 23
http://orcid.org/0000-0003-2465-3818
Hu S. X. 8
http://orcid.org/0000-0001-9112-0558
Glenzer S. H. 5
Kemp A. J. 7
http://orcid.org/0000-0002-4879-9072
Ping Y. 7
Sentoku Y. 4
1 https://ror.org/01keh0577 grid.266818.3 0000 0004 1936 914X Department of Physics, University of Nevada, Reno, Reno, NV USA
2 https://ror.org/01xjv7358 grid.410592.b 0000 0001 2170 091X Japan Synchrotron Radiation Research Institute, Hyogo, Japan
3 grid.472717.0 RIKEN SPring-8 Center, Hyogo, Japan
4 https://ror.org/035t8zc32 grid.136593.b 0000 0004 0373 3971 Institute of Laser Engineering, Osaka University, Suita, Osaka Japan
5 https://ror.org/05gzmn429 grid.445003.6 0000 0001 0725 7771 SLAC National Accelerator Laboratory, Menlo Park, CA USA
6 https://ror.org/0160cpw27 grid.17089.37 Department of Electrical and Computer Engineering, University of Alberta, Edmonton, AB Canada
7 https://ror.org/041nk4h53 grid.250008.f 0000 0001 2160 9702 Lawrence Livermore National Laboratory, Livermore, CA USA
8 https://ror.org/022kthw22 grid.16416.34 0000 0004 1936 9174 Laboratory for Laser Energetics, University of Rochester, Rochester, NY USA
5 9 2024
5 9 2024
2024
15 752814 3 2022
29 7 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
High-intensity, short-pulse lasers are crucial for generating energetic electrons that produce high-energy-density (HED) states in matter, offering potential applications in igniting dense fusion fuels for fast ignition laser fusion. High-density targets heated by these electrons exhibit spatially non-uniform and highly transient conditions, which have been challenging to characterize due to limitations in diagnostics that provide simultaneous high spatial and temporal resolution. Here, we employ an X-ray Free Electron Laser (XFEL) to achieve spatiotemporally resolved measurements at sub-micron and femtosecond scales on a solid-density copper foil heated by laser-driven fast electrons. Our X-ray transmission imaging reveals the formation of a solid-density hot plasma localized to the laser spot size, surrounded by Fermi degenerate, warm dense matter within a picosecond, and the energy relaxation occurring within the hot plasma over tens of picoseconds. These results validate 2D particle-in-cell simulations incorporating atomic processes and provide insights into the energy transfer mechanisms beyond current simulation capabilities. This work significantly advances our understanding of rapid fast electron heating and energy relaxation in solid-density matter, serving as a key stepping stone towards efficient high-density plasma heating and furthering the fields of HED science and inertial fusion energy research using intense, short-pulse lasers.

Intense, short-pulse laser irradiation generates energetic electrons that heat targets to extreme conditions relevant to laser fusion. Here, authors used an X-ray Free Electron Laser to perform spatiotemporal measurements in solid-density copper foil with sub-micron and femtosecond resolutions.

Subject terms

Laser-produced plasmas
High-field lasers
Free-electron lasers
https://doi.org/10.13039/100000001 National Science Foundation (NSF) PHY-1707357 PHY-2010502 Sawada H. https://doi.org/10.13039/501100001691 MEXT | Japan Society for the Promotion of Science (JSPS) 19K03788 22K03571 17H02996 17H02996 17H02996 17H06202 JPJSCCB20190003 17H06202 JPJSCCB20190003 17H02996 17H06202 JPJSCCB20190003 JP20H00140 20K14439 JP19KK0072 JP20H00140 20K14439 JP19KK0072 Yabuuchi T. Iwasaki T. Kawasaki K. Maeda Y. Izumi T. Nakagawa Y. Shigemori K. Sakawa Y. Iwata N. Sentoku Y. https://doi.org/10.13039/100006207 DOE | SC | Fusion Energy Sciences (FES) FWP 100182 FWP 100182 FWP 100182 Curry C. B. Frost M. Glenzer S. H. https://doi.org/10.13039/100000015 U.S. Department of Energy (DOE) DE-AC52-07NA27344 DE-NA0004144 DE-AC52-07NA27344 DE-AC52-07NA27344 Ogitsu T. Hu S. X. Kemp A. J. Ping Y. https://doi.org/10.13039/100006132 DOE | Office of Science (SC) SCW 1265 SCW 1420 Kemp A. J. Ping Y. issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

High-intensity, short-pulse lasers, with intensities exceeding 1018 W/cm2, efficiently generate high-charge, Mega electron-volt (MeV) electron currents during laser-matter interactions1,2. These electron beams not only advance the development of laser-driven electron accelerators3–5 and the generation of secondary X-ray and γ-ray sources6,7 but also create matter under extreme conditions. The properties of such materials, resulting from interactions with non-thermal electrons ranging from keV to MeV energies, are of great interest across various disciplines, including astrophysics8,9, aerospace engineering10, and both magnetic and inertial confinement fusion (ICF) research11–13. Rapidly heated matter by intense short-pulse lasers, with pulse durations less than a picosecond, undergoes an instantaneous increase in electron temperatures while the ions remain cold, leading to a non-equilibrium state14–16. This rapid heating process, known as isochoric heating, is pivotal for the fast ignition laser fusion approach17,18 in ICF19,20, employing an additional high-intensity laser to generate a beam of energetic charged particles (electrons, protons21,22, or other ions23). These particles deposit their energy directly into a separately compressed fuel core, initiating fusion reactions for ignition before fuel disassembly. Moreover, matter heated by laser-driven charged particles provides a crucial platform for benchmarking transport and radiative property models, including resistivity24, opacity25,26, and ion-stopping powers27 within the warm dense matter (WDM)28 regime. This regime lies at the boundaries among classical plasmas, Fermi degenerate matter, and strongly coupled matter. Studying the properties of degenerate matter across different temperatures and densities is crucial for achieving high fusion gains in inertial fusion energy29, as compressed cold fuels inherently exist in such states30,31.

Understanding fast electron isochoric heating in solid and high-density plasmas faces long-standing challenges due to the presence of preformed plasma and multiple, competing heating and ionization processes that vary across both space and time. Unlike interactions at non-relativistic intensities32, a low-intensity pedestal or the rising edge of the main pulse produces a thin plasma layer, preventing the high-intensity pulse from directly interacting with the solid target. Instead, laser-plasma interactions accelerate free electrons in the preplasma predominantly via the J×B force. As illustrated in Fig. 1a, these fast electrons, upon penetrating and moving transversely through a micron-thick solid-density foil, initiate heating and ionization through various mechanisms33,34: Joule heating occurs as the fast electrons induce a return current of cooler electrons, leading to resistive heating; drag heating results from direct collisions among non-thermal and background electrons, especially at high densities above 1025 cm−3; and diffusive heating takes place at the interface of the hot preplasma and the cold solid target, primarily near the laser interaction area. Additionally, the mean ionization state (denoted as Z¯) is influenced not only by thermal ionization but also by electron impact ionization, meaning that Z¯ varies independently from the electron temperature (Te). Consequently, a laser-irradiated thin solid foil exhibits at least three different plasma conditions: a solid-density region near the laser spot primarily heated by Joule and diffusive heating and ionized by thermal and electron impact ionization processes (Region 1), the outside region in the foil affected by Joule heating and impact ionization (Region 2), and the low-density preplasma. Thus, inferring electron temperatures from space-integrated X-ray emission spectroscopy could be misleading unless the origin of the X-ray sources is identified35. To date, indirect measurements and simulations suggest the spatial non-uniformity of dense plasmas within laser-irradiated solid-density foils36–38. However, a comprehensive understanding of the dynamics within a solid-density foil under fast electron heating—through the rapid heating and energy relaxation phases between electrons and ions—necessitates direct simultaneous measurements of Te and Z¯ with sufficient spatial and temporal resolutions.Fig. 1 Fast electron isochoric heating mechanisms in a solid foil and experimental layout.

a Schematic illustration of a high-intensity, short-pulse laser interacting with a thin, solid-density foil. Energetic electrons generated in a preplasma traverse the foil due to sheath potentials, leading to the creation of two distinct regions attributed to different heating and ionization mechanisms: a high-temperature, solid-density area in Region 1 by Joule and diffusive heating and thermal and impact ionization, and Fermi degenerate warm dense matter in Region 2 by Joule heating and electron impact ionization. b Layout of the experimental setup featuring a high-intensity femtosecond (fs) laser and an XFEL pulse. The imaging detector positioned outside the vacuum chamber, ~4.1 m away from the target. Examples of X-ray transmission images recorded (c) before and (d) after laser irradiation with a delay of 0.35 ps. e A ratio of these images provides a raw transmission ratio. An averaged horizontal lineout of the enhanced transmission area is superposed in (e).

Here, we present the investigation of the spatiotemporal dynamics of laser-driven fast electron heating in a solid-density copper foil, achieving sub-micron and femtosecond resolutions using an X-ray Free Electron Laser (XFEL). The hard X-ray pulses from the XFEL penetrate the preplasma and directly probe the solid-density region of the foil. Leveraging the tunability of XFEL photon energy and the spatial coherence of its beam, we developed a novel X-ray transmission imaging technique to diagnose the solid copper foil heated by an expanding fast electron heat front and the onset of electron-to-ion energy relaxation. Specifically, we inferred the electron temperature of the heated copper foil to be between 7 and 18 eV, based on the smearing observed in an X-ray transmission spectrum near the Cu K-edge. The ionization state, estimated to be between 2.0 and 4.0, was deduced by comparing the time evolution of the heat front with two-dimensional particle-in-cell (PIC) simulations incorporating an electron impact ionization model. Our spatiotemporal measurements and simulations reveal the formation of a solid-density hot plasma the size of the laser spot, surrounded by Fermi degenerate, WDM within a picosecond, and the energy relaxation occurring within the hot plasma over tens of picoseconds. These results not only validate our 2D PIC simulations but also provide important insights into the energy transfer processes in solid-density matter beyond current simulation capabilities. Experimental validation of PIC simulations in fast-electron-heated solid-density targets is a key stepping stone for efficient heating of high-density fusion fuels for fast ignition. This work also underscores the XFEL-based spatiotemporally resolved capability in high-energy-density science using intense short-pulse lasers, paving the way for generating well-characterized WDM and fusion-relevant hot dense matter at temperatures above 1 keV by laser-driven fast electrons.

Results

High-intensity femtosecond laser experiment at SACLA

We conducted a pump-probe experiment at Experimental Hutch 6 of the SACLA beamline39,40 using a high-intensity femtosecond laser and an XFEL beam. As shown in Fig. 1b, the experimental setup consisted of a 2 μm thick, solid-density copper foil irradiated by the laser at ~2 × 1018 W/cm2 for fast electron isochoric heating. A collimated X-ray pulse subsequently probed the laser-irradiated foil within a diameter spot of ~1 mm. The resultant 2D X-ray transmission image was captured by a scintillator-based CMOS detector located outside the vacuum chamber. Fig. 1c and d display the X-ray images taken at a probe energy of 9.05 keV before and after laser irradiation, respectively. The raw image (Fig. 1d) exhibits only marginal changes in transmission, which becomes more pronounced in the image obtained by taking the ratio of the two images (Fig. 1e). The non-uniform background was attributed to spatial and intensity fluctuations originating from the self-amplified spontaneous emission source. After correcting for uneven signals, we extracted a horizontal lineout of the transmission-affected area (indicated by the white line in Fig. 1e) to compare with 2D simulations that incorporate the 45° laser incident angle. By varying X-ray probe energies (9.05, 8.92, 9.12, and 8.05 keV) and delay timings from −1 ps to 200 ps, we measured the temporal evolution of the size of the electron-heated area. These transmission ratio images were used to deduce the electron temperature and ionization state within the heated foil. It is noted that the laser irradiation did not alter the solid copper foil, i.e., by punching a hole or evaporating it, as the X-ray images clearly show the edges of the foil even at a 500 ps timing delay and no diffraction fringes at earlier delay timings. Furthermore, the generation of fast electrons was monitored using a time-integrated X-ray spectrometer, monochromatic X-ray imager, and an electron spectrometer. The mean energy of the fast electron spectrum generated in this experiment was estimated to be ~100 keV, as determined by monochromatic X-ray and XFEL imaging, as reported in our previous publication41.

X-ray transmission imaging results

Figure 2 presents a series of transmission ratio images captured at three different X-ray energies (9.05, 8.92, and 8.05 keV), corresponding to timing delays of ~0.4 ps, 1.0 ps, and ~100 ps post-laser irradiation. Notably, the transmission ratio in the area affected by electron heating shows a clear dependence on probe energy: it increases at 9.05 keV (above the K-edge) and decreases at 8.92 keV (below the K-edge), as illustrated by the lineouts. The experimental feature observed at both 9.05 and 8.92 keV expands until ~1 ps, beyond which no significant growth is detected. The analysis of these lineouts and the estimation of width are detailed in the “Method” section of Experimental and Data Analysis Details. The width of the transmission-affected area was found to be consistent at 9.05 and 8.92 keV for the same delay timings. In contrast, at 8.05 keV and 9.12 keV (see below), no change in transmission was observed. The monochromatic Kα imager and X-ray spectrometer verified that the quantity and energy of the fast electrons generated remained consistent across all probe energies. These experimental images were used to infer the electron temperature and ionization state, a discussion of which follows in subsequent sections. Concentric diffraction fringes, visible in Fig. 2c, f, and i, were observed around 50 ps and at later delays regardless of the probe energies. These fringes likely result from the spatially coherent X-ray beam illuminating a locally ablated target surface. Qualitative comparisons between the measured and calculated fringe patterns suggest that the size of the ablated region is on the order of the laser spot (~26 µm in diameter). This finding, much smaller than the measured transmission-affected area of ~120 µm presented later, implies that only the region near the laser spot is significantly heated, radiatively cooled, and thermalized, causing the ablation over tens of picoseconds, while the ions in the outside region remain within the initial foil (see the “Diffraction fringe analysis” of the “Methods” section). The ratio of ~5 between the size of the hot plasma region and the transmission-affected area is consistent with past literature36. It is important to note that the 50 ps threshold represents an upper bound for the electron-ion energy relaxation and may be specific to this particular laser-target interaction. Further insights into the size of the deformed area, along with the ion temperature and collision frequency in the intensely heated region, are detailed in the “Methods” section (collision frequency in WDM).Fig. 2 Measured X-tray transmission images and estimating the electron temperature.

a–i Experimental transmission ratio images at 9.05, 8.92, and 8.05 keV for various timing delays. Averaged horizontal lineouts of the image around Y = 0 µm are superposed in each image. j A transmission profile of a 2 μm thick cold solid Cu foil and the Fermi-Dirac distribution at electron temperatures of 10 and 30 eV. The shaded vertical bars indicate the XFEL’s photon energies used with a spectral bandwidth of 50 eV. Transmission changes increase at 9.05 and 9.12 keV, while it decreases at 8.92 keV. k Relative changes in transmission are calculated using the Fermi–Dirac distribution as a function of Te. The shading indicates the minimum and maximum temperature range set by the spectral bandwidth for each photon energy.

Electron temperature estimation

The transmission ratio images in Fig. 2, varying with probe energy, enable the determination of the electron temperature in the heated area. This process leverages the sensitivity of the K-edge transmission (absorption) spectrum to electron temperature in Fermi degenerate materials like metals—a phenomenon explored in previous research42–45. Theoretically, the transmission slope near the K-edge is predicted to follow a Fermi–Dirac distribution, expressed as fE=1+e(E−μ)/kTe−1, where μ, the chemical potential, is approximately equal to the Fermi energy (EF). Figure 2j shows the X-ray transmission profile for a 2 µm thick, cold solid-density Cu foil alongside the Fermi-Dirac distribution for 10 and 30 eV electron temperatures. Our measurements of the changes in transmission at 8.92 and 9.05 keV reflect the changes in electron temperature of the distribution, underscoring that the target remains in a degenerate state. The absence of transmission change at 9.12 keV suggests an upper-temperature limit. To quantitatively determine a range of electron temperatures, we calculated the relative changes in transmission at each X-ray probe energy by varying electron temperature, as illustrated in Fig. 2k. Given a detection threshold at 5% of the peak signal, we determined that our imaging technique can detect relative signal changes above 0.15%—details are provided in the “Experimental and data analysis details” of the “Method” section. This analysis allowed us to deduce a lower temperature limit of 7.3 eV, necessary to observe transmission changes at both 8.92 and 9.05 keV, and to identify an upper temperature bound of 17.9 eV, the minimum temperature at which transmission changes at 9.12 keV become discernible. This marks the first determination of electron temperature bounds from the slope of the smeared K-edge profile in fast electron heated targets.

2D particle-in-cell simulations

To understand the dynamics of fast electron heating mechanisms and determine the ionization state of the electron-heated area, we performed two-dimensional collisional PIC simulations using the PICLS code46,47. Figures 3a and b illustrate these simulations, providing time-dependent spatial profiles of electron temperature (Te) and ionization state (Z¯), averaged over a 2 μm thick Cu target, with further details available in the “Methods” section (2D PIC Simulation Details). The evolution of the Te and Z¯ profiles proceeds as follows: immediately after the laser irradiation onto the foil, which generates fast electrons, Joule heating occurs near the focal spot (around Y =~0 μm) as these electrons penetrate the depth of the target. Upon reaching the rear of the target, most electrons are confined by the sheath potential and traverse laterally, raising the temperature to a range of 6–10 eV. This simulated peripheral temperature is consistent with the experimentally deduced temperature presented in Fig. 2k. Despite multiple collisions between fast and bulk electrons due to electron recirculation in the thin foil, the temperature outside the intensely heated region remains unchanged, indicating a negligible impact from drag heating in solid-density matter. After 0.56 ps, an increase in electron temperature near the laser interaction region suggests that diffusive heating is occurring, caused by the thermal gradient between the colder solid and the warmer preplasma. The diffusion process begins after the laser pulse has ended because laser-induced magnetic fields trap electrons within the preplasma. Additionally, Fig. 3a and b reveal a noticeable discrepancy between the spatial distributions of Te and Z¯, attributed to electron impact ionization.Fig. 3 Results of 2-D particle-in-cell simulations.

a Simulated electron temperature and (b) average ionization state profiles averaged over a 2 μm thick solid Cu foil from 2D PIC simulations (see Fig. 9). The evolution of these profiles is presented from 0.16 ps to 1.44 ps. Deviations of the ionization state profile from the temperature profile indicate the contribution of electron impact ionization. c Comparison of the widths of the measured transmission-affected area with simulated ionization states. The circle and square symbols represent the data at 9.05 keV and 8.92 keV, respectively. A dashed line indicates the speed of light (v = c) for reference. The error bars in time are ± 0.1 ps based on the uncertainty of the XFEL and the optical laser timing jitter, while the errors in width are ± 7.9 μm due to averaging of the experimental line profiles. The solid lines are results from the PIC simulation, and a constant extension of the simulation is assumed as extended dashed lines after ~1.4 ps.

Ionization state estimation

As the simulations suggest, the observed changes in transmission can be more accurately explained by the evolution of the ionization profile, which broadens over time, rather than by the temperature profiles. Figure 3c illustrates the progression of the transmission-affected area within the copper foil, as recorded at 9.05 keV and 8.92 keV, alongside simulated ionization states. During heating, the expansion of the heated area approaches the speed of light, indicating that the heat front is driven by near-relativistic electrons41. This rapid expansion stagnates after ~0.5 ps, and the experimental feature remains unchanged until the appearance of diffraction fringes around 50 ps (with further details up to 100 ps provided in the “Experimental and data analysis details” in “Methods”). In Fig. 3c, the experimental measurements are bounded by the simulations with an ionization state (Z¯) ranging from 2.0 to 4.0. Note that the inferred ionization state occurs behind the propagating heat front, with the overall spatial ionization profile peaking in the high-temperature region, as shown in Fig.3b. The relatively large measurement uncertainties in width, particularly noticeable during the heating phase (e.g., at ~0.3 ps), may be due to laser arrival time drift over several hours40.

Discussion

Our study, utilizing spatiotemporally resolved measurements alongside 2D PIC simulations, has unveiled the rapid transitions of solid-density cold copper foils to highly charged hot plasma and WDM. Figure 4 illustrates the temporal evolution of the foil’s conditions on an electron temperature-density diagram, categorized by the electron coupling parameter, Γe, and the degree of degeneracy, Θ, which are defined as Γe=(4πne/3)1/3e2/kBTe and Θ=kBTe/kBTF, where kB is the Boltzmann constant and the Fermi temperature, TF, is (ℏ2/2mekB)(3π2ne)2/3. Our measurements show the formation of a highly ionized hot region near the laser focal spot, surrounded by Fermi degenerate, WDM. The transition to these conditions is described by the trajectories derived from the simulated spatial profiles at Y = 0 µm and Y = 60 µm in Fig. 3a and b, as representatives for Region 1 and Region 2. Near the laser focal area (Y = 0 µm), the foil undergoes rapid heating to high temperatures through Joule heating (indicated by the purple arrow), evolving into diffusive heating above 100 eV. This region experiences changes in ionization state due to both thermal and electron impact ionization, consequently increasing the electron density through ne=ni⋅Z¯ even though the ion density remains constant. The plasma condition in Region 1 is non-degenerate (Te > TF), so the smearing of the X-ray spectrum near the K-edge was not applicable. In contrast, the peripheral region at Y = 60 µm is subject to rapid Joule heating to several eV in ~0.1 ps, followed by electron-impact ionization (shown by the red arrow), reaching up to Z¯ ~ 3, with a slight subsequent temperature increase. Our measurements of this peripheral region agree with the simulation results, as indicated by a white box in the figure. The simulations indicate that measuring the initial rise in the electron temperature due to Joule heating requires a faster resolution than our current 0.1 ps. The target conditions of the other parts of the foil fall between these two extremes. This comparison validates the Joule heating and electron impact ionization physics models incorporated in our numerical code. However, further experiments are needed to verify the high-temperature plasma and to examine the effects of spatial gradients in the depth direction of the foil.Fig. 4 Electron temperature and density contour for isochorically heated matter using an intense short-pulse laser.

The temperature and density map is divided using the coupling constant, Γe, and the degree of degeneracy, Θ, both of which are defined in the text. The shading indicates that matter at high temperatures and low densities with Γe<< 1 and Θ > 1 is treated as classical plasmas, while cold, dense matter (Γe> 1 and Θ < 1) is considered to be strongly coupled and Fermi degenerate matter. The warm dense matter regime is loosely defined in the vicinity of Γe~1 and Θ~1. The experimentally deduced conditions are depicted with a white box. Trajectories represented by red and purple lines correspond to PIC simulations at different distances from the laser interaction region, aligning with the axes in Fig. 3a and b. The evolving conditions at Y = 0 µm and Y = 60 µm are presented as representatives of conditions near the laser interaction in Region 1 and further from the center in Region 2. Filled circles on the solid lines indicate the simulation time of 0.15 ps.

The temperature-density diagram visualizes the two distinctly different high-energy-density states accessible with an intense short-pulse laser for unique applications. Diffusive heating plays a pivotal role in generating high-temperature, highly charged plasmas. In this experiment, using a 0.64 J, 40 fs optical laser, we achieved an efficiency of ~0.6% in heating the solid-density copper foil within a 4 µm spot diameter above 100 eV. Further research with state-of-the-art petawatt-class, short-pulse lasers could facilitate the creation of hot dense matter exceeding temperatures of 1 keV48, providing crucial insights into the underlying physics for fast ignition laser fusion. Additionally, Fermi-degenerate, WDM, produced by significant electron impact ionization, exhibits continuously evolving ionization states in the transverse direction of the foil, while the electron temperature remains nearly constant as shown in Figs. 3a and 4. These regions, where ionization states are the primary variable, could serve as an ideal testbed for benchmarking quantum molecular dynamics and atomic physics calculations, including models of ionization potential depression49 for dense matter.

This work has demonstrated the XFEL’s ability to perform detailed, spatiotemporally resolved measurements at the sub-micron scale and femtosecond resolutions in high-intensity short-pulse laser experiments. These diagnostic capabilities, vastly more precise than those typically employed in high-energy-density and ICF experiments29, play a crucial role in closing existing knowledge gaps. Importantly, the versatility of this measurement technique across various FEL beam wavelengths paves the way for future investigations with next-generation XFEL facilities, such as MEC-U50.

Methods

Experimental and data analysis details

In experimental hutch 6 of the SACLA XFEL facility, the high-intensity femtosecond laser system comprises two beams, each with a potential of 0.5 PW. One beam, delivering 12.5 J in 25 fs, has successfully commissioned and operates in conjunction with the XFEL beam, producing 8 J in 40 fs. However, due to damage in the laser optics, the energy delivered to the target in this experiment was limited to 0.6 ± 0.1 J. The pulse duration and laser spot size were 40 fs (full width at half maximum, FWHM) and 15 µm × 20 µm (FWHM), respectively. This resulted in a peak intensity on the target of I = 2 × 1018 W/cm2 and a relativistic factor (γ) of 1.2, where γ is defined as 1+a02/2 and a0 is 0.85Iλμm2/1018.

For the pump-probe experiment, the incident angles for the high-intensity optical pump laser and the XFEL probe laser were set at 45° and 27°, respectively, relative to the normal of a 2 μm-thick, 0.5 × 3.0 mm2 copper strip. The X-ray pulse had a nominal duration of ~10 fs and an energy of 0.6 mJ. The timing delay of the XFEL beam with respect to the fs laser varied from −0.1 ps to 200 ps. Additionally, X-ray transmission images of the laser-irradiated copper strips were recorded at the photon energies of 8.92, 9.05, 8.05, and 9.12 keV with a spectral bandwidth of ~50 eV. Figure 5 shows the transmission ratio images with the horizontal lineouts at 9.12 keV for delays of 0.4, 1.0, and 100 ps. At earlier delays, no clear variations in transmission were observed, whereas concentric fringe patterns appeared after ~50 ps, similar to the case at 8.05 keV. The shot rate for the experimental campaign was approximately one shot every 3 min.Fig. 5 Measured X-tray transmission ratio images at 9.12 keV.

The transmission ratio images are recorded using a photon energy of 9.12 keV at timing delays of (a) 0.4 ps, (b) 1.0 ps, and (c) 100 ps.

Time-integrated 2D monochromatic Kα images and X-ray spectral measurements were used as follows: The mean energy of the fast electrons (~100 keV) was estimated based on the radius of the Kα emission spot compared to the electron’s collisional stopping range, calculated using the continuous-slowing-down approximation51. With the X-ray spectrometer, we recorded characteristic emission lines at 8.05 keV (Cu Kα) and 8.9 keV (Cu Kβ). However, no line emissions from highly ionized copper ions, such as Heα, Lyα, or ionized Kα lines, were observed. This information indicates an upper bound of the electron temperature within the copper foil.

The lineouts of the transmission ratio image, as presented in Fig. 2a and d, are further detailed in Fig. 6 to illustrate how the width of the transmission-affected area was estimated. After correcting for the non-uniform background and setting the base signal to be zero, we first calculated a threshold level at 5% of the averaged peak signals41. This detection threshold corresponds to thresholds in transmission changes after the base adjustment to be between 0.15% and 0.23% in this experiment. With these thresholds, we determined not only the width of the experimental signature as shown in Fig. 6a and b but also the minimum thickness variation detectable by our imaging technique. Transmission through a 2 µm thick solid-density copper foil at 8.05 keV is 91.3%. Using a tabulated absorption coefficient52, we calculated relative changes in X-ray transmission, as illustrated in Fig. 6c, indicating that our measurement technique is sensitive to foil thickness changes greater than 32.4 nm. This minimum thickness threshold was later used to discuss the collision frequency and ion temperature.Fig. 6 Analysis of the experimental lineouts and detection threshold.

The lineouts of the experimental features shown in Fig. 2a and d for 9.05 keV and 8.92 keV are presented in (a) and (b). The width of the transmission-affected region is estimated using a threshold set at 5% of the peak signal. Due to fluctuations in the measured signals, the detection threshold varied between 0.15% and 0.23% in transmission changes over the data analyzed. c An estimate of the minimum changes in foil thickness is considered by calculating changes in X-ray transmission at 8.05 keV through a solid copper foil with a various thickness of (2 µm – ∆L). The shading indicates the detection threshold range (0.15–0.23%). With these thresholds, this figure suggests that our imaging technique detects changes in transmission when the foil thickness varies from the initial 2 µm thickness by more than 32.4 nm.

Figure 7 presents an extended comparison of the width of the experimental features with the simulated ionization states over a period of up to 100 ps from Fig. 3c. The experimental widths at 1 ps and 10–20 ps are similar, suggesting that no significant target disassembly occurs before 10 ps. This observation is supported by the appearance of concentric diffraction fringes around 50 ps. Data that are strongly interfered with diffraction fringes for delays between 50 and 200 ps are excluded. The notable shot-to-shot variations observed, particularly at 10 and 100 ps, are likely due to fluctuations in laser conditions and target positions. The laser arrival time is known to drift over several hours, even within a single shot cycle40, which also likely contributes to the data scatter observed at 0.2–0.4 ps. Additionally, these variations could occur when measurements are taken in different shot cycles, such as one measurement during a single shot cycle and another following the venting of the vacuum chamber and adjustment of the laser focus spot. The experimental data, particularly after 1.5 ps, are bounded by ionization states between Z¯ = 3.0 and Z¯ = 4.0. However, extrapolating the simulated ionization state as constant over 100 ps oversimplifies the ionization process. Further details such as radiative cooling, recombination, and cooling by plasma expansion need to be considered.Fig. 7 Comparison of the extended experimental data upto 100 ps.

The time history of the widths of the measured transmission-affected (ionized) regions is compared with simulated ionization states ranging between Z¯ = 2.0 and 4.0. The circle and square symbols correspond to the data at 9.05 keV and 8.92 keV, respectively. The error bars in time and in width are ± 0.1 ps and ± 7.9 μm, respectively. Solid lines represent the PIC simulations up to 1.4 ps. Beyond this point, the simulation results are assumed to remain constant, as indicated by dashed lines.

Fresnel diffraction calculations

Concentric fringe patterns consistently appear at delays exceeding 50 ps, as shown in Fig. 2c, f, and i. These fringes are likely caused by the 2D coherent X-ray pulse illuminating a planar foil, with the deformation of the target surface resulting from ion motion (ablation), rather than from surface plasmons or laser-induced modulations. This is supported by Fig. 2f, where at least two sets of circular fringe patterns are overlapped at 150 ps. Here, we present Fresnel diffraction calculation results to estimate the approximate size of the ablated region within a 2 µm thick solid-density copper foil by qualitatively comparing the measured and calculated diffraction patterns. The experimental setup is essentially equivalent to a classical diffraction experiment, wherein a monochromatic, coherent plane wave strikes a semi-transparent flat target featuring a locally ablated region, differing from the scenario with an opaque aperture. Given the detector distance (L) of 4.12 ± 0.01 m from the target and the light wavelength (λ) of 0.154 nm (corresponding to 8.05 keV), the Fresnel number (NF), defined as πa2/(L⋅λ), was calculated to range from 0.12 and 12.38 for a characteristic aperture size (radius, a) ranging from 5 to 50 µm. Therefore, Fresnel diffraction calculations are applicable in this near-field diffraction regime, where NF is greater than approximately 153.

We modeled two-dimensional Fresnel diffraction of a circular transparent aperture using the transfer function method54. This approach was specifically utilized to compare with the image captured at 8.05 keV with a 50 ps timing delay, as shown in Fig. 8a. In our model, a plane wave with a wavelength of 0.154 nm illuminates a two-dimensional plane containing a circular spot. This spot has a transmission slightly different from the cold transmission of a 2 μm solid Cu foil at 8.05 keV. The model calculations for aperture radii of a = 5, 13, and 25 μm, alongside comparisons with the experimental measurements, are displayed in Fig. 8b–e. These comparisons reveal that the position of the first bright ring, as indicated by arrows in the figures, shifts outward with increasing aperture size, consequently affecting the diffraction intensity patterns. The calculations with an aperture radius of 13 ± 2 μm showed a qualitative agreement with the experimental measurement, as marked by the vertical dashed lines in Fig. 8e. The inferred size of the ablated region ( ~ 26 µm in diameter) corresponds closely to the laser spot size and a high-temperature region (near Y = 0 μm or Region 1), as illustrated in Figs. 3a and 1a. This observation provides direct evidence of which part of the heated target begins to disassemble first following equilibration.Fig. 8 Results of Fresnel diffraction calculations.

a A measured transmission ratio image at 8.05 keV at a timing delay of 50 ps. Fresnel diffraction intensity patterns are calculated with various aperture sizes, specifically aperture radii of (b) 5 μm, (c) 13 μm, and (d) 25 μm. e Comparison of the azimuthally averaged measured fringe profile with the calculations. The arrows indicate the spatial positions of the first bright fringe for the different aperture radii.

The timing of the electron-ion equilibration from our measurements is consistent with findings from other studies that have established an upper limit for the surface ablation timing55,56. In particular, our measurements are close to the duration of time-resolved Kα emission measurements (10–20 ps)57 rather than that of K-shell emissions (several ps)36,37,58, indicating that the hot plasma that produces K-shell emissions (Heα, Heβ, Lyα, etc) during a high-intensity, laser-solid interaction is radiatively cooled, which diminishes the X-ray emission, before target expansion occurs. In the next section, we calculate the evolution of electron temperature and ionization state due to radiative cooling and consider it in estimating an electron-ion collision frequency in warm dense copper.

2D PIC simulation details

We performed two-dimensional PIC simulations using the fully relativistic, collisional PIC code PICLS, as depicted in Fig. 3. The simulation setup and contours of electron temperature and mean ionization state at 0.16 ps is shown in Fig. 9. PICLS self-consistently solves the relativistic equations of motion for numerous particles and Maxwell’s equations, addressing laser-target interactions, high-energy electron generation, and electron transport. It includes modules for Coulomb collisions and dynamic ionizations, the primary ionization mechanisms. The simulations used the experimental laser conditions: an 800 nm wavelength, a 40 fs Gaussian pulse (FWHM), P-polarized light at a 45° incidence angle, and a peak intensity of 2 × 1018 W/cm2. The simulation includes a preplasma with a scale length of 0.43 μm within a 2 μm layer. The simulation encompassed a 40 × 160 µm2 system box with 2000 × 8000 cells, running up to ~1.5 ps. The simulation predicts that the high-temperature solid-density plasma in Region 1 reaches Te = 120 eV and Z¯ = 20. The collision frequency and electron-to-ion energy relaxation time are calculated to be 4.9 × 1016 s−1 and 1–2 ps, respectively. The latter value is an order of magnitude smaller than the observation 50 ps in this experiment. This discrepancy provides further evidence that radiative cooling must occur prior to energy relaxation.Fig. 9 PIC simulation setup.

a A schematic of the 2D PIC simulation setup. Contours of (b) the simulated electron temperature on a log scale and (c) the ionization state at 0.16 ps.

Radiative cooling calculations

The evolution of the target conditions due to radiative cooling is calculated over tens of picoseconds, beyond the capabilities of PIC simulations. We computed the radiative cooling of non-degenerate plasmas over time by numerically solving the equation netdTetdt=ηTet,ni), where ηTet,ni) represents the emissivity. This emissivity, accounting for free-free, bound-free, and bound-bound emissions across various electron temperatures at the solid copper ion density, was precomputed using the collisional-radiative atomic physics code FLYCHK59 for temperatures between 100 eV and 1 keV. The electron density, ne(t)=Z¯(Te)⋅ni, was updated based on FLYCHK results. Figure 10 illustrates the time evolution of Te and Z¯ for plasmas with initial electron temperatures of 100, 300, and 1000 eV. Notably, Te significantly decreases within several picoseconds after the Gaussian laser pulse ends around 0.1 ps, as indicated by a dashed vertical line.Fig. 10 Radiative cooling calculations.

Radiative energy loss of a solid-density hot Cu plasma is calculated for initial electron temperatures of 100, 300, and 1000 eV up to 50 ps. The temporal evolution of (a) the electron temperatures and (b) the ionization states shows significant decreases in both quantities between 1 and 10 ps. c Comparison of the temporal decay in the electron temperature and Fermi temperature with constant density. The rapid drop in the electron temperature leads to the transition of the plasma to a degenerate state (TF > Te). The vertical dashed line indicates the transition time of 14.3 ps for the 100 eV case.

In high-temperature, thermally ionized plasma (e.g., around Y = 0 μm in Fig. 3a), the ionization state decreases over time with Te. For example, in the 100 eV case, Te falls below the Fermi temperature (TF) after 14.3 ps, as shown in Fig. 10c, because the decrease in Te is faster than that in the ionization state. The corresponding transition times from non-degenerate to degenerate plasmas for 300 eV and 1000 eV are 18.4 ps and 20.1 ps, respectively. These results suggest that even a highly charged, thermally ionized hot plasma, typically considered classical, can transition into Fermi degenerate matter due to radiative cooling over tens of picoseconds, irrespective of the initial temperature. Considering multi-dimensional effects, this transition could occur even more rapidly. Thus, it is plausible that the target near the laser interaction region (∣Y∣<∼13 μm in Fig. 3a) was in a degenerate state when the surface ablation occured. The radiative cooling calculations show that a plasma with an initial condition of Te of 100 eV and Z¯ of 12 decreases to Te of 6.7 eV and Z¯ of 4.3 at 50 ps. These values at 50 ps are used in the following section to infer the ion temperature and an empirical parameter for the electron-ion collision frequency in Fermi degenerate copper plasma.

Collision frequency in warm dense matter

The electron-ion collision frequency for WDM (νwdm)60,61 significantly differs from that for ideal plasmas and is given by1 νwdm=Kwdme2ℏvFkBTiℏ

where Ti is the ion temperature, kB is the Boltzmann constant, vF is the Fermi velocity, and Kwdm is an empirical constant. In this section, we determined the ion temperature from the velocity of ablated ions and an empirical constant, Kwdm, in the above equation.

The characteristic time, τ, of ~50 ps when the significant ion motion is observed allows us to estimate the collision frequency at 1.2 × 1015 s−1 through τ=Mi/2meνei, where Mi and me are the ion and electron masses, respectively. The Fermi velocity, defined as vF=3π2ne1/3ℏ/me, changes in proportion to the ionization state at a constant density. In the previous section, the radiatively cooled plasma condition shows a Te of 6.7 eV and Z¯ of 4.3 at 50 ps. Using this ionization state and the initial ion density, we calculate the Fermi velocity to be 2.6 × 108 cm/s.

Our imaging technique, as shown in Fig. 6c, detects changes in thickness once the 2 µm thick foil thickness changes more than 32.4 nm. We hypothesize that all ions within a 32.4 nm depth (and a 13 µm radius) are ablated around 50 ps, leading to an ion velocity of 6.5 × 104 cm/s and a corresponding ion temperature of 0.28 eV. Applying these values to Eq. (1), we found the Kwdm constant to be ~3.2. Several values of the empirical constant have been reported: 4.6 for solid-density, room-temperature aluminum, and 18.8 for warm dense aluminum under equilibrium plasma conditions (Te = Ti)60. Although further research is required to accurately determine the electron-ion collision frequency in warm dense copper, this work has demonstrated that the XFEL-based time-resolved measurements provide new insights into electron-ion energy relaxation processes and subsequent hydrodynamic expansion induced by high-intensity, femtosecond lasers at both non-relativistic and relativistic intensities62.

Supplementary information

Peer Review File

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-51084-4.

Acknowledgements

The XFEL experiments were performed at the BL2 of SACLA with the approval of the Japan Synchrotron Radiation Research Institute (JASRI) (Proposal No. 2019A8002). This material is based on the work supported by the National Science Foundation under Grant No. 1707357 and 2010502 through the NSF/DOE Partnership in Basic Plasma Science and Engineering. Y. Sentoku and N.I. are supported by JSPS Kakenhi Grants No. JP20H00140, 20K14439, 24H00204, and JP19KK0072. T.Iwasaki, K.K., Y.M., and K.Shigemori are supported by JSPS Kakenhi Grants No.17H02996. Part of the work was performed under the auspices of the US Department of Energy by Lawrence Livermore National Laboratory under contract no. DE-AC52-07NA27344 (A.J.K., T.O., and Y.P.) and funded by the DOE Office of Science Early Career Program under SCW 1265 (A.J.K.) and SCW 1420 (Y.P.). The work of C.B.C., M.F., and S.H.G. was supported by the U.S. Department of Energy (DOE), Office of Science, Fusion Energy Sciences under FWP 100182. C.B.C. was also partially supported by the Natural Sciences and Engineering Research Council of Canada (NSERC). T. Izumi, Y. N., and Y. Sakawa were supported by JSPS Kakenhi Grants No.17H06202 and JSPS Core-to-Core Program B: Asia-Africa Science Plat-forms Grant No. JPJSCCB20190003. S.X.H acknowledges the support of the Department of Energy under award number DE-NA0004144. T.Y. is supported by JSPS Kakenhi Grants No. 19K03788 and 22K03571.

Author contributions

H.S. led the writing the manuscript with the help of the coauthors’ inputs and feedback. H.S., T.Y., T.Iwasaki, K.K., Y.M., T.Izumi, Y.N., K.Shigemori, and Y.Sakawa, performed the experiment and acquired the data. N.H., N.I., Y.Sentoku, S. X. H., and H.S. performed the theoretical and numerical work. T.Y., K.Sueda, and T.T. operated the high-power femtosecond laser. C.B.C and M.F. performed laser cutting of the Cu foils. H.S., S.H.G., A.J.K., Y.P., T.O., and Y.Sentoku developed the original experimental concept. All authors contributed to the work presented here and to the final manuscript.

Peer review

Peer review information

Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. A peer review file is available.

Data availability

Experimental data were generated at the SACLA XFEL facility. The key data underpinning the conclusions of this work are included in the article. All raw data relevant to the findings of this study are available from the corresponding author (H.S.) upon request.

Code availability

The authors declare that the computer code PICLS supporting the findings of this study is fully documented within the paper and its references. Additional inquiries about the codes should be directed to Y. Sentoku.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Wilks SC Kruer WL Tabak M Langdon AB Absorption of ultra-intense laser pulses Phys. Rev. Lett. 1992 69 1383 1386 10.1103/PhysRevLett.69.1383 10047203
Wilks, S. C., Kruer, W. L., Tabak, M. & Langdon, A. B. Absorption of ultra-intense laser pulses. Phys. Rev. Lett. 69, 1383–1386 (1992).10047203 10.1103/PhysRevLett.69.1383
2. Tikhonchuk VT Interaction of a beam of fast electrons with solids Phys. Plasmas 2002 9 1416 10.1063/1.1459061
Tikhonchuk, V. T. Interaction of a beam of fast electrons with solids. Phys. Plasmas 9, 1416 (2002).10.1063/1.1459061
3. Faure J A laser-plasma accelerator producing monoenergetic electron beams Nature 2004 431 541 544 10.1038/nature02963 15457253
Faure, J. et al. A laser-plasma accelerator producing monoenergetic electron beams. Nature 431, 541–544 (2004).15457253 10.1038/nature02963
4. Geddes CGR High-quality electron beams from a laser wakefield accelerator using plasma-channel guiding Nature 2004 431 538 541 10.1038/nature02900 15457252
Geddes, C. G. R. et al. High-quality electron beams from a laser wakefield accelerator using plasma-channel guiding. Nature 431, 538–541 (2004).15457252 10.1038/nature02900
5. Mangles SPD Monoenergetic beams of relativistic electrons from intense laser-plasma interactions Nature 2004 431 535 538 10.1038/nature02939 15457251
Mangles, S. P. D. et al. Monoenergetic beams of relativistic electrons from intense laser-plasma interactions. Nature 431, 535–538 (2004).15457251 10.1038/nature02939
6. Giulietti A Intense γ-ray source in the giant-dipole-resonance range driven by 10-TW laser pulses Phys. Rev. Lett. 2008 101 105002 10.1103/PhysRevLett.101.105002 18851220
Giulietti, A. et al. Intense γ-ray source in the giant-dipole-resonance range driven by 10-TW laser pulses. Phys. Rev. Lett. 101, 105002 (2008).18851220 10.1103/PhysRevLett.101.105002
7. Sarri G Table-top laser-based source of femtosecond, collimated, ultrarelativistic positron beams Phys. Rev. Lett. 2013 110 255002 10.1103/PhysRevLett.110.255002 23829742
Sarri, G. et al. Table-top laser-based source of femtosecond, collimated, ultrarelativistic positron beams. Phys. Rev. Lett. 110, 255002 (2013).23829742 10.1103/PhysRevLett.110.255002
8. Bulanov SV On the problems of relativistic laboratory astrophysics and fundamental physics with super powerful lasers Plasma Phys. Rep. 2015 41 1 51 10.1134/S1063780X15010018
Bulanov, S. V. et al. On the problems of relativistic laboratory astrophysics and fundamental physics with super powerful lasers. Plasma Phys. Rep. 41, 1–51 (2015).10.1134/S1063780X15010018
9. Fleishman GD Revealing the evolution of non-thermal electrons in solar flares using 3D modeling Astrophys. J. 2018 859 17 10.3847/1538-4357/aabae9
Fleishman, G. D. et al. Revealing the evolution of non-thermal electrons in solar flares using 3D modeling. Astrophys. J. 859, 17 (2018).10.3847/1538-4357/aabae9
10. Shang JJS Modeling Plasma via Electron Impact Ionization Aerosp 2018 5 2 10.3390/aerospace5010002
Shang, J. J. S. Modeling Plasma via Electron Impact Ionization. Aerosp 5, 2 (2018).10.3390/aerospace5010002
11. Breizman BN Aleynikov P Hollmann EM Lehnen M Physics of runaway electrons in tokamaks Nucl. Fusion 2019 59 083001 10.1088/1741-4326/ab1822
Breizman, B. N., Aleynikov, P., Hollmann, E. M. & Lehnen, M. Physics of runaway electrons in tokamaks. Nucl. Fusion 59, 083001 (2019).10.1088/1741-4326/ab1822
12. Kidder RE Hot-electron preheat of laser-driven targets Nucl. Fusion 1981 21 145 10.1088/0029-5515/21/2/002
Kidder, R. E. Hot-electron preheat of laser-driven targets. Nucl. Fusion 21, 145 (1981).10.1088/0029-5515/21/2/002
13. Betti R Hurricane OA Inertial-confinement fusion with lasers Nat. Phys. 2016 12 435 448 10.1038/nphys3736
Betti, R. & Hurricane, O. A. Inertial-confinement fusion with lasers. Nat. Phys. 12, 435–448 (2016).10.1038/nphys3736
14. Jiang Z X‐ray spectroscopy of hot solid density plasmas produced by subpicosecond high contrast laser pulses at 1018–1019 W/cm2 Phys 1995 2 1702 1711
Jiang, Z. et al. X‐ray spectroscopy of hot solid density plasmas produced by subpicosecond high contrast laser pulses at 1018–1019 W/cm2. Phys 2, 1702–1711 (1995).
15. Workman J Control of bright picosecond X-ray emission from intense subpicosecond laser-plasma interactions Phys. Rev. Lett. 1995 75 2324 2327 10.1103/PhysRevLett.75.2324 10059275
Workman, J. et al. Control of bright picosecond X-ray emission from intense subpicosecond laser-plasma interactions. Phys. Rev. Lett. 75, 2324–2327 (1995).10059275 10.1103/PhysRevLett.75.2324
16. Saemann A Isochoric heating of solid aluminum by ultrashort laser pulses focused on a tamped target Phys. Rev. Lett. 1999 82 4843 4846 10.1103/PhysRevLett.82.4843
Saemann, A. et al. Isochoric heating of solid aluminum by ultrashort laser pulses focused on a tamped target. Phys. Rev. Lett. 82, 4843–4846 (1999).10.1103/PhysRevLett.82.4843
17. Tabak M Ignition and high gain with ultrapowerful lasers Phys. Plasmas 1994 1 1626 10.1063/1.870664
Tabak, M. et al. Ignition and high gain with ultrapowerful lasers. Phys. Plasmas 1, 1626 (1994).10.1063/1.870664
18. Kodama R Fast heating of ultrahigh-density plasma as a step towards laser fusion ignition Nature 2001 412 798 802 10.1038/35090525 11518960
Kodama, R. et al. Fast heating of ultrahigh-density plasma as a step towards laser fusion ignition. Nature 412, 798–802 (2001).11518960 10.1038/35090525
19. Nuckolls J Wood L Thiessen A Zimmerman G Laser compression of matter to super-high densities: thermonuclear (CTR) applications Nature 1972 239 139 142 10.1038/239139a0
Nuckolls, J., Wood, L., Thiessen, A. & Zimmerman, G. Laser compression of matter to super-high densities: thermonuclear (CTR) applications. Nature 239, 139–142 (1972).10.1038/239139a0
20. Atzeni, S. & Meyer-ter-Vehn, J. The Physics of Inertial Fusion: Beam Plasma Interaction, Hydrodynamics, Hot Dense Matter 1st edn, (Oxford Univ. Press, 2004); 10.1093/acprof:oso/9780198562641.001.0001.
21. Roth M Fast ignition by intense laser-accelerated proton beams Phys. Rev. Lett. 2001 86 436 439 10.1103/PhysRevLett.86.436 11177849
Roth, M. et al. Fast ignition by intense laser-accelerated proton beams. Phys. Rev. Lett. 86, 436–439 (2001).11177849 10.1103/PhysRevLett.86.436
22. Ditmire T Focused energy, a new approach towards inertial fusion energy J. Fusion Energy 2023 42 27 10.1007/s10894-023-00363-x
Ditmire, T. et al. Focused energy, a new approach towards inertial fusion energy. J. Fusion Energy 42, 27 (2023).10.1007/s10894-023-00363-x
23. Fernández JC Progress and prospects of ion-driven fast ignition Nucl. Fusion 2009 49 65004 10.1088/0029-5515/49/6/065004
Fernández, J. C. et al. Progress and prospects of ion-driven fast ignition. Nucl. Fusion 49, 65004 (2009).10.1088/0029-5515/49/6/065004
24. Booth N Laboratory measurements of resistivity in warm dense plasmas relevant to the microphysics of brown dwarfs Nat. Commun. 2015 6 8742 10.1038/ncomms9742 26541650
Booth, N. et al. Laboratory measurements of resistivity in warm dense plasmas relevant to the microphysics of brown dwarfs. Nat. Commun. 6, 8742 (2015).26541650 10.1038/ncomms9742
25. Nazir K X-ray spectroscopic studies of hot, dense iron plasma formed by subpicosecond high intensity KrF laser irradiation Appl. Phys. Lett. 1996 69 3686 3688 10.1063/1.117189
Nazir, K. et al. X-ray spectroscopic studies of hot, dense iron plasma formed by subpicosecond high intensity KrF laser irradiation. Appl. Phys. Lett. 69, 3686–3688 (1996).10.1063/1.117189
26. Davidson SJ Nazir K Rose SJ Smith R Tallents GJ Short-pulse laser opacity measurements J. Quant. Spectrosc. Radiat. Transf. 2000 65 151 160 10.1016/S0022-4073(99)00063-1
Davidson, S. J., Nazir, K., Rose, S. J., Smith, R. & Tallents, G. J. Short-pulse laser opacity measurements. J. Quant. Spectrosc. Radiat. Transf. 65, 151–160 (2000).10.1016/S0022-4073(99)00063-1
27. Malko S Proton stopping measurements at low velocity in warm dense carbon Nat. Commun. 2022 13 2893 10.1038/s41467-022-30472-8 35610200
Malko, S. et al. Proton stopping measurements at low velocity in warm dense carbon. Nat. Commun. 13, 2893 (2022).35610200 10.1038/s41467-022-30472-8
28. Lee RW Finite temperature dense matter studies on next-generation light sources J. Opt. Soc. Am. B 2003 20 770 10.1364/JOSAB.20.000770
Lee, R. W. et al. Finite temperature dense matter studies on next-generation light sources. J. Opt. Soc. Am. B 20, 770 (2003).10.1364/JOSAB.20.000770
29. IFE Basic Research Needs Final Report (2023), https://science.osti.gov/-/media/fes/pdf/workshop-reports/2023/IFE-Basic-Research-Needs-Final-Report.pdf.
30. Hu SX Militzer B Goncharov VN Skupsky S Strong coupling and degeneracy effects in inertial confinement fusion implosions Phys. Rev. Lett. 2010 104 235003 10.1103/PhysRevLett.104.235003 20867248
Hu, S. X., Militzer, B., Goncharov, V. N. & Skupsky, S. Strong coupling and degeneracy effects in inertial confinement fusion implosions. Phys. Rev. Lett. 104, 235003 (2010).20867248 10.1103/PhysRevLett.104.235003
31. Hu SX Impact of first-principles properties of deuterium–tritium on inertial confinement fusion target designs Phys. Plasmas 2015 22 56304 10.1063/1.4917477
Hu, S. X. et al. Impact of first-principles properties of deuterium–tritium on inertial confinement fusion target designs. Phys. Plasmas 22, 56304 (2015).10.1063/1.4917477
32. Gamaly EG Rode AV Physics of ultra-short laser interaction with matter: from phonon excitation to ultimate transformations Prog. Quantum Electron. 2013 37 215 323 10.1016/j.pquantelec.2013.05.001
Gamaly, E. G. & Rode, A. V. Physics of ultra-short laser interaction with matter: from phonon excitation to ultimate transformations. Prog. Quantum Electron. 37, 215–323 (2013).10.1016/j.pquantelec.2013.05.001
33. Glinsky ME Regimes of suprathermal electron transport Phys. Plasmas 1995 2 2796 2806 10.1063/1.871472
Glinsky, M. E. Regimes of suprathermal electron transport. Phys. Plasmas 2, 2796–2806 (1995).10.1063/1.871472
34. Kemp AJ Sentoku Y Sotnikov V Wilks SC Collisional relaxation of superthermal electrons generated by relativistic laser pulses in dense plasma Phys. Rev. Lett. 2006 97 235001 10.1103/PhysRevLett.97.235001 17280208
Kemp, A. J., Sentoku, Y., Sotnikov, V. & Wilks, S. C. Collisional relaxation of superthermal electrons generated by relativistic laser pulses in dense plasma. Phys. Rev. Lett. 97, 235001 (2006).17280208 10.1103/PhysRevLett.97.235001
35. Sawada H Monochromatic 2D K α emission images revealing short-pulse laser isochoric heating mechanism Phys. Rev. Lett. 2019 122 155002 10.1103/PhysRevLett.122.155002 31050520
Sawada, H. et al. Monochromatic 2D K α emission images revealing short-pulse laser isochoric heating mechanism. Phys. Rev. Lett. 122, 155002 (2019).31050520 10.1103/PhysRevLett.122.155002
36. Audebert P Heating of thin foils with a relativistic-intensity short-pulse laser Phys. Rev. Lett. 2002 89 265001 10.1103/PhysRevLett.89.265001 12484827
Audebert, P. et al. Heating of thin foils with a relativistic-intensity short-pulse laser. Phys. Rev. Lett. 89, 265001 (2002).12484827 10.1103/PhysRevLett.89.265001
37. Audebert P Time-resolved plasma spectroscopy of thin foils heated by a relativistic-intensity short-pulse laser Phys. Rev. E 2002 66 066412 10.1103/PhysRevE.66.066412
Audebert, P. et al. Time-resolved plasma spectroscopy of thin foils heated by a relativistic-intensity short-pulse laser. Phys. Rev. E 66, 066412 (2002).10.1103/PhysRevE.66.066412
38. Faenov AY Nonlinear increase of X-ray intensities from thin foils irradiated with a 200 TW femtosecond laser Sci. Rep. 2015 5 13436 10.1038/srep13436 26330230
Faenov, A. Y. et al. Nonlinear increase of X-ray intensities from thin foils irradiated with a 200 TW femtosecond laser. Sci. Rep. 5, 13436 (2015).26330230 10.1038/srep13436
39. Tono K Beamline, experimental stations and photon beam diagnostics for the hard x-ray free electron laser of SACLA N. J. Phys. 2013 15 083035 10.1088/1367-2630/15/8/083035
Tono, K. et al. Beamline, experimental stations and photon beam diagnostics for the hard x-ray free electron laser of SACLA. N. J. Phys. 15, 083035 (2013).10.1088/1367-2630/15/8/083035
40. Yabuuchi T An experimental platform using high-power, high-intensity optical lasers with the hard X-ray free-electron laser at SACLA J. Synchrotron Radiat. 2019 26 585 594 10.1107/S1600577519000882 30855271
Yabuuchi, T. et al. An experimental platform using high-power, high-intensity optical lasers with the hard X-ray free-electron laser at SACLA. J. Synchrotron Radiat. 26, 585–594 (2019).30855271 10.1107/S1600577519000882
41. Sawada H Ultrafast time-resolved 2D imaging of laser-driven fast electron transport in solid density matter using an x-ray free electron laser Rev. Sci. Instrum. 2023 94 033511 10.1063/5.0130953 37012804
Sawada, H. et al. Ultrafast time-resolved 2D imaging of laser-driven fast electron transport in solid density matter using an x-ray free electron laser. Rev. Sci. Instrum. 94, 033511 (2023).37012804 10.1063/5.0130953
42. Dorchies F X-ray absorption K edge as a diagnostic of the electronic temperature in warm dense aluminum Phys. Rev. B 2015 92 085117 10.1103/PhysRevB.92.085117
Dorchies, F. et al. X-ray absorption K edge as a diagnostic of the electronic temperature in warm dense aluminum. Phys. Rev. B 92, 085117 (2015).10.1103/PhysRevB.92.085117
43. Dorchies F Recoules V Non-equilibrium solid-to-plasma transition dynamics using XANES diagnostic Phys. Rep. 2016 657 1 26 10.1016/j.physrep.2016.08.003
Dorchies, F. & Recoules, V. Non-equilibrium solid-to-plasma transition dynamics using XANES diagnostic. Phys. Rep. 657, 1–26 (2016).10.1016/j.physrep.2016.08.003
44. Hansen SB Fluorescence and absorption spectroscopy for warm dense matter studies and ICF plasma diagnostics Phys. Plasmas 2018 25 056301 10.1063/1.5018580
Hansen, S. B. et al. Fluorescence and absorption spectroscopy for warm dense matter studies and ICF plasma diagnostics. Phys. Plasmas 25, 056301 (2018).10.1063/1.5018580
45. Rosmej FB Astapenko VA Khramov ES XFEL and HHG interaction with matter: Effects of ultrashort pulses and random spikes Matter Radiat. Extrem 2021 6 034001 10.1063/5.0046040
Rosmej, F. B., Astapenko, V. A. & Khramov, E. S. XFEL and HHG interaction with matter: Effects of ultrashort pulses and random spikes. Matter Radiat. Extrem 6, 034001 (2021).10.1063/5.0046040
46. Sentoku Y Kemp AJ Numerical methods for particle simulations at extreme densities and temperatures: weighted particles, relativistic collisions and reduced currents J. Comput. Phys. 2008 227 6846 6861 10.1016/j.jcp.2008.03.043
Sentoku, Y. & Kemp, A. J. Numerical methods for particle simulations at extreme densities and temperatures: weighted particles, relativistic collisions and reduced currents. J. Comput. Phys. 227, 6846–6861 (2008).10.1016/j.jcp.2008.03.043
47. Mishra R Leblanc P Sentoku Y Wei MS Beg FN Collisional particle-in-cell modeling for energy transport accompanied by atomic processes in dense plasmas Phys. Plasmas 2013 20 072704 10.1063/1.4812701
Mishra, R., Leblanc, P., Sentoku, Y., Wei, M. S. & Beg, F. N. Collisional particle-in-cell modeling for energy transport accompanied by atomic processes in dense plasmas. Phys. Plasmas 20, 072704 (2013).10.1063/1.4812701
48. Higashi N Iwata N Sano T Mima K Sentoku Y Isochoric heating of solid-density plasmas beyond keV temperature by fast thermal diffusion with relativistic picosecond laser light Phys. Rev. E 2022 105 055202 10.1103/PhysRevE.105.055202 35706231
Higashi, N., Iwata, N., Sano, T., Mima, K. & Sentoku, Y. Isochoric heating of solid-density plasmas beyond keV temperature by fast thermal diffusion with relativistic picosecond laser light. Phys. Rev. E 105, 055202 (2022).35706231 10.1103/PhysRevE.105.055202
49. Hu SX Continuum lowering and fermi-surface rising in strongly coupled and degenerate plasmas Phys. Rev. Lett. 2017 119 65001 10.1103/PhysRevLett.119.065001
Hu, S. X. Continuum lowering and fermi-surface rising in strongly coupled and degenerate plasmas. Phys. Rev. Lett. 119, 65001 (2017).10.1103/PhysRevLett.119.065001
50. Dyer, G. et al. The MEC-U Project at LCLS. in CLEO 2023 JTh2A.72 (Optica Publishing Group, 2023). 10.1364/CLEO_AT.2023.JTh2A.72.
51. Berger M. J. et al. Stopping-power and range tables for electrons, protons and helium ions (Gaithersburg, MD: NIST Physics Laboratory, 1998).
52. Hubbell, J. H. & Seltzer, S. M. Tables of X-ray mass attenuation coefficients and mass energy-absorption coefficients (version 1.4). (National Institute of Standards and Technology, Gaithersburg, MD, 2004) http://physics.nist.gov/xaamdi.
53. Pedrotti, F. L., Pedrotti, L. M. & Pedrotti, L. S. Introduction to Optics. 3rd edn, (Cambridge University Press, Cambridge, P274, 2017).
54. Li J Peng Z Fu Y Diffraction transfer function and its calculation of classic diffraction formula Opt. Commun. 2007 280 243 248 10.1016/j.optcom.2007.08.053
Li, J., Peng, Z. & Fu, Y. Diffraction transfer function and its calculation of classic diffraction formula. Opt. Commun. 280, 243–248 (2007).10.1016/j.optcom.2007.08.053
55. Zhang N Zhu X Yang J Wang X Wang M Time-resolved shadowgraphs of material ejection in intense femtosecond laser ablation of aluminum Phys. Rev. Lett. 2007 99 167602 10.1103/PhysRevLett.99.167602 17995294
Zhang, N., Zhu, X., Yang, J., Wang, X. & Wang, M. Time-resolved shadowgraphs of material ejection in intense femtosecond laser ablation of aluminum. Phys. Rev. Lett. 99, 167602 (2007).17995294 10.1103/PhysRevLett.99.167602
56. Li, J. et al. Ultrafast electron beam imaging of femtosecond laser-induced plasma dynamics. J. Appl. Phys. 107, 083305 (2010).
57. Shepherd, R. et al. Short pulse laser absorption and energy partition at relativistic laser intensities, (LLNL Report UCRL-TR-229579, 2007).
58. Kraus BF Solid-density ion temperature from redshifted and double-peaked stark line shapes Phys. Rev. Lett. 2021 127 205001 10.1103/PhysRevLett.127.205001 34860067
Kraus, B. F. et al. Solid-density ion temperature from redshifted and double-peaked stark line shapes. Phys. Rev. Lett. 127, 205001 (2021).34860067 10.1103/PhysRevLett.127.205001
59. Chung H-K Chen MH Morgan WL Ralchenko Y Lee RW FLYCHK: generalized population kinetics and spectral model for rapid spectroscopic analysis for all elements High. Energy Density Phys. 2005 1 3 12 10.1016/j.hedp.2005.07.001
Chung, H.-K., Chen, M. H., Morgan, W. L., Ralchenko, Y. & Lee, R. W. FLYCHK: generalized population kinetics and spectral model for rapid spectroscopic analysis for all elements. High. Energy Density Phys. 1, 3–12 (2005).10.1016/j.hedp.2005.07.001
60. Eidmann K Meyer-ter-Vehn J Schlegel T Hüller S Hydrodynamic simulation of subpicosecond laser interaction with solid-density matter Phys. Rev. E 2000 62 1202 10.1103/PhysRevE.62.1202
Eidmann, K., Meyer-ter-Vehn, J., Schlegel, T. & Hüller, S. Hydrodynamic simulation of subpicosecond laser interaction with solid-density matter. Phys. Rev. E 62, 1202 (2000).10.1103/PhysRevE.62.1202
61. Atzeni, S. & Meyer-ter-Vehn, J. The physics of inertial fusion: beam plasma interaction, hydrodynamics, hot dense matter (Clarendon, Oxford, 2004).
62. Gamaly EG The physics of ultra-short laser interaction with solids at non-relativistic intensities Phys. Rep. 2011 508 91 243 10.1016/j.physrep.2011.07.002
Gamaly, E. G. The physics of ultra-short laser interaction with solids at non-relativistic intensities. Phys. Rep. 508, 91–243 (2011).10.1016/j.physrep.2011.07.002
