I. INTRODUCTION
The investigation of states of matter with high energy density (>105 J/cm3) is of great interest across various scientific fields, including astrophysics,1,2 strong-interaction physics,3 plasma physics and thermonuclear fusion,4,5 physics at high voltage and high power,6–8 and particle acceleration techniques.9,10 Creating such extreme states of matter in the laboratory typically involves compression of heated matter or plasma using shock waves induced by various tools such as gas guns, pinching discharges, and high-power lasers.11–13 If it becomes possible to provide a high initial density of matter through the use of solid targets and rather short durations of the heating pulse (τlas ≪ τexp ∼ lplasma/vs, where lplasma is the plasma size and vs is the ion-acoustic velocity), then high-energy-density conditions can be created at least for times not exceeding the plasma expansion time τexp. Note that for times shorter than τexp, the expansion of the plasma is one-dimensional and its density decreases slowly, but for t > τexp, the plasma expansion becomes three-dimensional and its density decreases rapidly. When a laser pulse interacts with a target, the size of the directly heated area dplasma is almost equal to the focal spot and is 1–10 μm. The expansion velocity of the ions is the ion-acoustic velocity, which is not strongly dependent on the plasma temperature and has a typical value of about 107 cm/s. The condition τlas < τexp requires laser pulses with a duration of femtoseconds (τlas < 10–100 fs). To achieve this, femtosecond lasers with relativistic intensity can be used. There has been significant progress in the development of these lasers during the last ten years.14–17
The use of ultra-intense ultrashort pulses enables higher temperatures of the heated electrons and a higher electron density to be achieved.
However, it is possible to increase the density in the absorption region if an ultrashort pulse irradiates targets that have not undergone significant expansion on the front side. In this case, the laser energy couples directly into the solid-density material, and the electron density significantly exceeds the critical density. This enables the generation of a plasma with ultrahigh energy density (>108 J/cm3). However, it requires very high-contrast laser pulses to exclude target expansion during the laser prepulse.14
In this work, we investigate the conditions that allow the generation of an ultrahigh density, ultra-high-temperature plasma with an energy density of more than 108 J/cm3 and a pressure of more than 1 Gbar. We investigate plasma parameters measured by X-ray spectroscopic approaches when targets are driven with a femtosecond laser pulse with a laser intensity on the target of up to 5 × 1021 W/cm2. When the laser contrast is insufficient, but an ultrahigh laser pulse is maintained, the energy concentration in the plasma is still high, owing to the high relativistic critical density, but significantly lower than in the case of solid-density plasma formation. In practice, measuring the laser contrast in the 1012–1015 range is challenging, and it is difficult to determine the contrast level of single shots in advance. Nevertheless, by diagnosing the parameters of the generated plasma, we can deduce the state of the target at the time of the main laser pulse. If the target remains solid, this indicates that the laser contrast is sufficient to minimize the pre-plasma effects. In Ref. 18, it is shown that at a laser intensity of ∼1020 W/cm2, a typical laser contrast of 1010 is required to eliminate pre-plasma effects and that the requirements increase with increasing laser intensity.
II. EXPERIMENT
A. Experimental setup
The experiments were performed at the J-KAREN-P laser facility at the Kansai Photon Science Institute, Japan. The laser is a PW-class hybrid laser system that combines optical parametric chirped pulse amplification (OPCPA) and Ti:Sa chirped pulse amplification technology. In the current experiments, the J-KAREN-P laser operated at a wavelength of 800 nm and delivered high-contrast (up to 1012) pulses with a duration t ∼ 40 fs and an intensity Ilt of up to 5 × 1021 W/cm2 on the target.16,19,20
In our experiments, we systematically investigated the plasma of thin stainless-steel foils under different laser conditions: on-target intensity < 1021 W/cm2 and high contrast > 1010 (referred to as case A); on-target intensity > 1021 W/cm2 and low contrast 105–106 (case B); and on-target intensity > 1021 W/cm2 and high contrast ∼1012 (case C). The laser contrast was measured using a Sequoia cross-correlator (Amplitude Technologies). The parameters of the generated plasma, namely, its electron temperature Te and electron density Ne, were measured using high-resolution X-ray spectroscopy. A sketch of the experimental setup is shown in Fig. 1. The laser beam irradiated a stainless-steel foil (SUS, type AISI 304: Fe 72%, Cr 18%, Ni 10%) with thickness 2 or 5 μm at an angle of 45° to the normal to the target surface. The p-polarized beam was focused by an off-axis parabolic mirror with f/1.3 into a focal spot with diameter d ∼ 2–4 μm.

Figure 1.Experimental setup. A femtosecond laser pulse with wavelength λ = 800 nm is focused at an angle of 45° onto the surface of a steel foil target in a focal spot with a diameter of 2–4 μm to generate a plasma on the front surface of the target. The X rays from the plasma are measured by the X-ray focusing spectrometer (FSSR) positioned at a distance L = 2045 mm from the plasma source at an angle of 8° to the normal to the target surface.
The X-ray spectra of the plasma emission were detected by a focusing spectrometer with spatial resolution (FSSR)21 with a spectral resolution λ/δλ ∼ 3000. The spectrometer was installed at a distance of 2045 mm from the plasma source to provide registration of spectra emitted from the front side of the SUS foil surface at an observation angle ∼8° to the target surface normal.
The spectrometer was equipped with a spherically bent mica crystal with radius of curvature R = 150 mm and lattice spacing 2d = 19.94 Å or an alpha-quartz crystal with 2d = 2.36 Å (orientation 31–40, radius of curvature R = 150 mm). X-ray spectra were detected by an Andor DX-440 X-ray CCD with pixel size 13.5 μm. The CCD matrix was protected against exposure to visible light by two layers of 1 μm-thick polypropylene coated with 0.2 μm Al. In all experiments, the spectrometer was aligned to measure K-shell emission of multiply charged Fe xxv and Fe xxvi ions (lines Heα, and Lyα) and the characteristic lines Kα,β of the Fe i ion in the wavelength range 1.74–1.97 Å. This observation range corresponds to the eighth (m = 8) diffraction order of the mica crystal and the first diffraction order of the α-quartz crystal. To suppress the contribution of plasma emission in the low diffraction orders of the mica, we placed a Mylar (C10H8O4) film in front of the crystal. The thickness of the Mylar film was varied from 6 to 100 μm. As already mentioned, the experiments were carried out for different laser operating regimes that differed primarily in terms of laser contrast and laser intensity on the target surface. The parameters of the experiments are listed in Table I. The measured X-ray spectra emitted by the SUS plasma in the spectral range 1.73–1.97 Å are shown in Fig. 2. The X-ray spectra shown in Figs. 2(a)–2(c) correspond respectively to the experimental conditions of cases A, B, and C described in Table I.

Table 1. Main parameters of three experiments (cases A, B, and C) performed at the J-KAREN-P laser facility.
Table 1. Main parameters of three experiments (cases A, B, and C) performed at the J-KAREN-P laser facility.
Parameter | J-KAREN, case A | J-KAREN-P, case B | J-KAREN-P, case C |
---|
Energy (J) | 7 | 14 | 9 | Pulse duration (fs) | 35 | 35 | 40 | Focal spot diameter (μm) | 6.5 | 4 | 3 | Laser intensity (W/cm2) | 6 × 1020 | 3 × 1021 | 5 × 1021 | Laser contrast | >1010 | 105–106 | >1012 | Target thickness (μm) | 2 | 5 | 5 | FSSR crystal | Alpha-quartz, | Mica, | Alpha-quartz, | | 2d = 2.36 Å | 2d = 19.94 Å | 2d = 2.36 Å | X-ray spectrum | Figure 2(a) | Figure 2(b) | Figure 2(c) |
|

Figure 2.X-ray spectra of stainless-steel plasma measured under the experimental conditions described in Table I: (a) case A; (b) case B; (c) case C.
First, we compare the spectra measured in the experimental cases A and B. Figure 2(a) clearly shows that the intensity of the Fe Heα line is lower, although experiment B was performed with a higher laser pulse energy (14 vs 7 J in experiment A) and intensity (3 × 1021 W/cm2 vs 6 × 1020 W/cm2 in experimental case A. In the spectrum shown in Fig. 2(b), the Heα line of Fe xxv is very poorly separated from the noise, whereas it is well defined in the spectrum in Fig. 2(a).
In case B, the laser contrast is much poorer. This means that the drop in laser contrast to values of 105–106 leads to pre-plasma formation, which drastically (up to two orders of magnitude) decreases the density of the heated plasma and consequently increases the time required to complete the ionization processes by up to two orders. The observed absence of He- and H-like iron lines in the spectrum in Fig. 2(b) makes it impossible to apply X-ray spectroscopic approaches based on the consideration of line emission spectra. Nevertheless, information about the temperature of the bulk of the plasma electrons can be obtained by analyzing the continuous spectrum due to bremsstrahlung, as was done in Refs. 22 and 23. The detailed calculation of the X-ray spectrum measured in case A is presented in Ref. 24 and will not be discussed in detail here. The X-ray spectrum measured in case C is shown in Fig. 2(c). This spectrum is particularly interesting because, as we will see below, it corresponds to the case in which no pre-plasma is formed, and it becomes possible to heat the solid target to high keV temperatures and generate a plasma with an energy density of ∼1 GJ/cm3. In contrast to the spectra in Figs. 2(a) and 2(b), the spectrum in Fig. 2(c) contains intense lines corresponding to the transitions in H- and He-like Fe ions. These can be used to determine the plasma parameters using the following X-ray spectroscopic approaches.
B. Diagnostics approaches
Note that in experiments with ultrahigh laser intensities, the deviation of the contrast level from shot to shot, as well as the deviation of other pulse parameters, can be significant. Therefore, the question arises as to whether it is possible to determine the creation of the pre-plasma with other diagnostic methods.
We propose the use of two X-ray spectroscopic approaches for this purpose. The first approach allows a clear distinction of the conditions for the interaction of the laser pulse with the target when the main pulse is shielded by a pre-plasma or when no pre-plasma occurs, and an ultrahigh energy density plasma is generated. The second approach is based on registering and analyzing specific features that can only appear in the emission spectrum if the plasma has a high density.25,26 These are first, a broadening of the spectral lines due to the Stark effect and, second, the relative intensities of certain transitions in multiply charged ions (see, e.g., Ref. 27).
In our case, the spectral feature in the range of Fe xxv and Fe xxvi ions (lines Heα1 = 1.85 Å, Heα2 = 1.86 Å, and Lyα1 = 1.77 Å, Lyα2 = 1.78 Å) and characteristic lines Kα = 1.94 Å (Kβ = 1.76 Å) of Fe are under consideration. Note that since the measurements of X-ray emission from the femtosecond laser plasma were performed without time resolution, each observed spectrum represents a time-integrated plasma glow consisting of emissions corresponding to the decay of the dense plasma state. As a result, typical spectral features of a dense plasma could be overlapped by spectral features emitted at a later time of the plasma’s existence. Two diagnostic approaches can therefore be used to investigate the interaction between the laser pulse and the dense plasma.
The idea of this first approach is that the time of collisional ionization in the plasma is inversely proportional to the plasma density. Thus, if one knows to some approximation the lifetime of the plasma in the hot state and uses a time-dependent kinetics calculation, one can find maximum ionization states for each ion that correspond to the solid or critical electron plasma density.
If the experimental spectrum shows lines emitted by ions in states that can only be observed in a solid-density plasma, this clearly indicates the absence of a pre-plasma in this experiment. As the ionization time for solid and critical density would differ by two to three orders of magnitude, the inaccuracy in the plasma lifetime estimation may be considered negligible. This approach has already been used in Ref. 24 to investigate the influence of laser pulse contrast on femtosecond laser plasma formation in stainless-steel foils. The method described is suitable for diagnosing pre-plasma formation, although it cannot be used to obtain detailed information about the plasma parameters.
The second diagnostic approach is based on the comparison of experimentally measured spectra with corresponding modeled spectra. It provides more detailed information about the plasma parameters but is more complex and difficult to apply for the following reasons.
In our case, the plasma parameters change rapidly both along and perpendicular to the target surface. To achieve a high laser intensity, the laser beam is usually focused on a focal spot with a diameter of ∼1.5 to 10 μm. Accordingly, the plasma parameters change considerably in the same spatial dimensions. At such dimensions, the spectra of the plasma emission cannot be detected with spatial resolution. The measured spectra are therefore a superposition of X-ray spectra emitted by plasma regions (zones) with very different properties. Since the plasma parameters are continuous functions of the spatial coordinates, it is necessary for an adequate description of the observed spectra to consider a number of plasma zones with different parameters, as has been shown in Refs. 24 and 28–30. It is usually sufficient to consider two to three zones. The concept of plasma zones is illustrated in Fig. 3.

Figure 3.Sketch of femtosecond ultra-intense laser pulse interaction with foil in the plasma zone concept.
The first plasma zone, zone 1 obviously corresponds to the plasma region in which most of the energy of the laser pulse is absorbed. Zone 1 has a transverse size that corresponds to the size of the laser focal spot. If the wavelength of the heating laser pulse is ∼0.8 μm and its intensity is a relativistic intensity Ilas ∼ 1018 to 1021 W/cm2, the existence of a high-intensity laser prepulse leads to the formation of a pre-plasma in the focal spot region. This plasma can reach an electron density around the critical electron density.24,30 However, at ultra-relativistic laser intensities Ilas > 1021 W/cm2, the plasma density in zone 1 becomes much higher, and, under certain experimental conditions,22,31 increases to the value of the relativistic critical density or remains a solid density. Zone 2 is next to zone 1 and is heated by energy transfer processes from zone 1.
Plasma zone 3 is the area of the solid target at a considerable distance from the laser beam axis. The heating of zone 3 is mainly due to hot electrons or intense X-ray photons emitted from the first zone. We therefore assume that the ion density in zones 2 and 3 is equal to the solid density Ni2 = Ni3 = Ni,solid. However, the density in zone 1 could be a solid density Ni1 = Ni,solid if there were no laser prepulse, or, depending on the laser intensity and prepulse, a density corresponding to the relativistic critical density Ne1 = Ncr,rel or the nonrelativistic critical density Ne1 = Ne,cr. This approach suggests that the electron temperature decreases from zone 1 to zone 3, i.e., Te1 > Te2 > Te3.
If the observed spectrum can be described by a sum of emission spectra from zones with different plasma parameters and if the hottest zone has the electron density Ne1, then, by comparing Ne1 with the solid density and the relativistic critical density, it is possible to derive a plasma density in the region where the main laser pulse has mostly been absorbed. It is also possible to directly estimate the energy density of the free plasma particles. The validity of the second approach requires the presence of a sufficient number of spectral transitions in the measured spectra. In addition, a high spectral resolution of the X-ray spectra must be provided. The implementation of the second approach for the spectrum measured in case C is discussed in Sec. III B.
III. RESULTS AND DISCUSSION
A. Calculation of time dependence of formation of He-like and H-like Fe ions
First, we use the first diagnostic approach and consider the time dependences of the formation of He-like Fe xxv and H-like Fe xxvi ions in plasmas with different electron densities and electron temperatures (Fig. 4). Calculations are performed under the assumption of a time-dependent radiational–collisional kinetic model using the computational code PrismSPECT32–34 for the following cases: a solid-density plasma with temperature 3000 eV (shown in black in Fig. 4), relativistic-density plasmas at laser intensities 1 × 1021 W/cm2 (in blue), 3 × 1021 W/cm2 (in red), and 5 × 1021 W/cm2 (in green) with temperatures 1000 eV [Fig. 4(a)] and 2000 eV [Fig. 4(b)] and a critical-density plasma with temperature 3000 eV (in gray). In all cases, the dependences for He-like and H-like Fe ions are indicated by solid and dashed lines, respectively. The calculation results showing the typical plasma lifetime for the formation of He-like Fe xxv and H-like Fe xxvi ions are presented in Table II.

Figure 4.Time dependences of concentrations of He-like Fe xxv (solid lines) and H-like Fe xxvi (dashed lines) ions calculated for electron temperatures of 1000 eV (a) and 2000 eV (b) and critical plasma densities, correspondent to different laser intensities (in different colors). Black curves correspond to the case of solid density (Ni,solid = 8 × 1022 cm−3) and electron temperature Te = 3000 eV. Gray curves correspond to the case of nonrelativistic critical density (Ne,cr ∼ 1 × 1021 cm−3) and electron temperature Te = 3000 eV.

Table 2. Results of the calculation of the time dependence of the formation of He-like Fe xxv and H-like Fe xxvi ions in plasmas with different electron densities and electron temperatures.
Table 2. Results of the calculation of the time dependence of the formation of He-like Fe xxv and H-like Fe xxvi ions in plasmas with different electron densities and electron temperatures.
Plasma density type | Plasma density (cm−3) | Plasma temperature Te (eV) | Time for He-like Fe xxv ion formation (ps) | Time for H-like Fe xxvi ion formation (ps) |
---|
Nonrelativistic critical electron density | Ne,cr ∼ 1 × 1021 | 3000 | 10–20 | 70–100 | Solid ion density | Ni,solid = 8 × 1022 | | 0.01–0.02 | 0.05–0.08 | Relativistic critical electron density for Ilt ∼ 1 × 1021 W/cm2 | Nrel,cr ∼ 3 × 1022 | 2000 | 0.2–0.3 | 1–2 | Relativistic critical electron density for Ilt ∼ 3 × 1021 W/cm2 | Nrel,cr ∼ 5 × 1022 | | 0.1–0.2 | 0.6–1 | Relativistic critical electron density for Ilt ∼ 5 × 1021 W/cm2 | Nrel,cr ∼ 6.5 × 1022 | | 0.07–0.1 | 0.1–0.3 | Relativistic critical electron density for Ilt ∼ 1 × 1021 W/cm2 | Nrel,cr ∼ 3 × 1022 | 1000 | 1–2 | 10–20 | Relativistic critical electron density for Ilt ∼ 3 × 1021 W/cm2 | Nrel,cr ∼ 5 × 1022 | | 0.8–0.1 | 7–10 | Relativistic critical electron density for Ilt ∼ 5 × 1021 W/cm2 | Nrel,cr ∼ 6.5 × 1022 | | 0.2–0.3 | 2–4 |
|
As we can see from Table II and Fig. 4, the Fe ions with main K shell can only be produced at a reasonable plasma lifetime if the plasma density notably exceeds its nonrelativistic critical value Ne,cr ∼ 1021 cm−3. Generation of H-like Fe ions is only possible at laser intensities Ilt ≥ 1021 W/cm2, i.e., if the plasma electron density is not less than 3 × 1022 cm−3 and corresponds to the relativistic critical density. The formation of H-like Fe states becomes more effective in the case of a plasma with a solid density that is heated to a temperature of ∼3000 eV (see Table II).
As can be seen in Fig. 2(a), the laser contrast increases up to 1010, and so the spectral lines of He-like Fe xxv ions can be observed continuously even at a laser intensity Ilt ∼ 6 × 1020 W/cm2, which is lower than in case B. At the same time, the emission of H-like states is almost absent in the spectrum. From the calculations shown in Fig. 4(b), it follows that the plasma density in the focal spot region must have a relativistic critical density to enable such conditions to be achieved within a plasma lifetime of the order of a few picoseconds. This was also shown in Ref. 35.
Let us now look at the spectrum in Fig. 2(c), in which many intense Fe lines of the K shell can be seen. First, we note that the presence of bright emission of H-like Fe xxvi ions in the spectral region of the Lyα line, together with the calculations shown in Fig. 4(b), suggests that the plasma formed under experimental conditions in case C has a bulk electron temperature of the order of 2000 eV and a density higher than the critical density. The observation of intense well-resolved spectra of H- and He-like Fe ions allows us to apply our second approach to determine the plasma parameters.
B. Modeling of x-ray spectrum measured in experimental case C
To model the X-ray spectrum shown in Fig. 2(c), kinetic calculations are performed within the collisional–radiation model, including radiative transfer effects, using the computational code PrismSpect mentioned above. The plasma volume is divided into three plasma zones with significantly different parameters. Each zone is associated with the emission of specific spectral lines.
Let us consider X-ray emission from plasma zone 1. Since we have defined zone 1 as the hottest plasma region, the influence of K-shell emission from this zone is most important and should largely describe the emission of He- and H-like Fe ions. Good agreement between the modeling and the experimental spectrum can best be achieved by varying the values of Te1, Ne1, and the plasma thickness l1. In our considerations, the transverse size of the first plasma zone dplasma1 cannot be larger than the diameter of the laser focal spot. The parameter dplasma1 is limited by dplasma ≤ 2 μm. The longitudinal thickness of the plasma l1 is defined by the expansion of the plasma perpendicular to the surface of the target. Assuming a one-dimensional expansion, the plasma thickness in zone 1 is estimated as l1 ∼ 0.5 μm in the case of Ne1 = Nsolid or l1 ∼ 2 μm in the case of Ne1 = Nrel,cr (for a laser intensity of 5 × 1021 W/cm2). The optical thickness of the plasma is considered using the Biberman–Holstein approximation for the escape factor.
We must emphasize that at an electron density corresponding to the nonrelativistic critical density Ne1 ∼ 1021 cm−3, the spectrum in Fig. 2(c) cannot be described by any values of plasma temperature or plasma size, including the case of an optically thin plasma (l → 0). Figure 5 clearly shows this situation. The modeling in the large Te1 range from 1000 to 5000 eV (the red lines in Fig. 5) does not match the X-ray spectrum measured in experimental case C. Here and in all other figures, the X-ray spectrum measured in case C [see Fig. 2(c)] is shown filled in with gray shading for better illustration and easier differentiation from the modeling results.

Figure 5.Dependence of X-ray spectra emitted from plasma zone 1 with critical electron density (Ne1 = Ne,cr = 1 × 1021 cm−3) on the bulk electron temperature of the plasma Te1; (a) Te1 = 5000 eV; (b) Te1 = 3000 eV; (c) Te1 = 1000 eV. The modeling results are shown by red lines. The X-ray spectrum measured in experimental case C is shown filled in with gray shading.
At first glance, the choice of the relativistic critical density [see Fig. 6(a)] allows the ratio of the relative intensities of the Lyα and Heα resonance lines to be described. However, as in the case of nonrelativistic critical density, the model calculations for the relativistic critical density cannot describe the experimental spectrum in the wavelength range 1.87–1.90 Å, in which the Li-like satellites of the Heα line are located. The situation changes significantly if it is assumed that the plasma generated in zone 1 has a solid density Ni,solid = 8 × 1022 cm−3. For this case, the modeling results for different electron temperatures are shown in Fig. 6(b).

Figure 6.Dependence of X-ray spectra emitted from plasma zone 1 on bulk plasma electron temperature: (a) for electron density Ne = Nrel,cr = 6.5 × 1022 cm−3 (the laser intensity of 5 × 1021 W/cm2) and plasma thickness l1 ∼ 2 μm; (b) for solid density Ni,solid = 8 × 1022 cm−3 and plasma thickness l1 = 0.5 μm. The modeling results are shown by red lines. The X-ray spectrum from Fig. 2(c) is shown filled in with gray shading.
As can be seen in Fig. 6(b), we can describe the experimental spectrum near the Lyα line and the experimental values of the relative intensities of the resonance lines (Lyα1 and Heα1). Good agreement between modeling and measurements in the wide spectral range λ = 1.75–1.87 Å is achieved at a bulk plasma electron temperature Te1 = 1800 eV and a plasma thickness l1 = 0.5 μm.
Since there are no lines associated with transitions in less-ionized Fe ions (Be-, B-, C- and N-like states) for this Te1 in the modeling spectrum, but these lines are present in the measured spectrum, we have considered the emission of these lines as radiation from plasma zone 2, which has a lower temperature. The modeling shows that the contribution of the satellite emission to the resulting intensity in the right wing of the Heα2 line depends strongly on the temperature and density (Figs. 7 and 8). Figure 7 shows the emission spectra of Fe ions calculated at a solid plasma ion density Ni = 8 × 1022 cm−3 for different plasma temperatures. The 0.1% of hot electrons with a temperature Thot = 10 keV were included in the calculations. We will show that the plasma temperature that allows a description of the Li-like Fe ion states is about Te2 ∼ 1000 eV and that a temperature Te2 ∼ 500 eV is required to describe the Be- and B-like states.

Figure 7.Emission spectra of Fe plasma calculated at solid plasma ion density Ni,solid = 8 × 1022 cm−3 for different plasma temperatures Te. The 0.1% of hot electrons with temperature Thot = 10 keV are included in the calculations. The X-ray spectrum from Fig. 2(c) is shown filled in with gray shading.

Figure 8.Emission spectra of Fe plasma calculated as superpositions of emission spectra from zone 1 (Te1 = 1800 eV) with Ne1 = Nsolid and from zone 2 at fixed plasma solid ion density Ni,solid = 8 × 1022 cm−3, for different plasma temperatures, with the inclusion of 0.1% of hot electrons with temperature Thot = 10 keV. The X-ray spectrum from Fig. 2(c) is shown filled in with gray shading.
Let us now fix the ion density Ni,solid = 8 × 1022 cm−3 and consider together the modeling obtained for zone 1 at electron temperature Te1 = 1800 eV [Fig. 6(b), second from top] and the above result for zone 2. Figure 8 shows the sequence of spectra obtained as the sum of the spectra of zones 1 and 2 for different Te2 in the range 500–900 eV. The spectra are considered with equal weight coefficients. Figure 8 shows that the situation is still complicated. By adding a second plasma zone with Te2 ∼ 500–900 eV, either the measured spectrum in the wavelength range 1.89–1.91 Å can be fitted, but with the description of the Li-like lines in the range 1.87–1.89 Å remaining incomplete, or conversely we can obtain good agreement between experiment and modeling for intense Li-like satellites, but with a large discrepancy in the range 1.88–1.91 Å.
Note that in both cases the relative intensity of the spectral lines in the range of 1.85–1.87 Å is significantly higher than the experimental values. The reason for this unsuccessful description of the spectrum is that the value chosen above Te1 = 1800 eV for zone 1, which correctly describes Lyα with the satellites and correctly reflects the ratio of the relative intensities of the resonant Lyα1 and Heα1 lines, was set too low.
This refers to the earlier assumption that the emission of He-like ions is only from the first plasma zone. However, in a solid-density plasma, spectral features associated with transitions in helium-like Fe ions appear to a small extent already at a temperature Te2 = 400 eV [see Fig. 6(b)]. For this reason, the resonance line Heα1 cannot be considered as a spectral feature emitted only by zone 1. Therefore, it is necessary to find a Te1 at which the ratio Lyα1/Heα1 corrects the modeled spectrum near Heα1 and ensures good agreement in the spectral region including Lyα and its satellites. The contribution to the sum of the modeled spectrum from each plasma zone must remain proportional.
It is obvious that choosing a higher Te1 requires a higher value for Te2. For Te1 ∼ 2000 eV, the relative intensities of the intercombination line Heα2 and the Li-like satellites are already lower than the intensity of the resonance line Heα1 [see Fig. 6(b)] and do not correspond to the experimental values. This inevitably leads to the choice of a Te2 that corresponds to an intense excitation of Li-like states and provides a sufficient contribution for zone 2 in the sum modeling the spectrum in the range 1.87–1.90 Å. Since lines of hydrogen-like Fe ions appear in the spectrum at Te2 ≥ 1100 eV, this implies that the effects of zone 2 in the spectral region include Lyα and its satellites. We also note that the choice of electron temperature in zone 1 has an upper limit of ∼3000 eV, since at Te1 > 2500 eV, the contribution of zone 1 to the total spectrum is considerably lower than the contribution of zone 2.
In a few iterations, we gradually increase the plasma temperature in zone 1 from the lowest value Te1 = 1800 eV to Te1 = 3000 eV and adjust the corresponding electron temperature in the range 400 eV ≤ Te2 ≤ 1100 eV. We find that the best agreement between the modeled and measured spectra is achieved in a wide spectral range when considering the electron temperatures Te2 = 2500 eV and Te2 = 900 eV (see Fig. 9).

Figure 9.Emission spectra of iron plasma calculated as superposition of emission spectra from zones 1 and 2 at fixed plasma solid ion density Ni,solid = 8 × 1022 cm−3 and a temperature in zone 1 of Te1 = 2500 eV for different plasma temperatures Te2. The X-ray spectrum from Fig. 2(c) is shown filled in with gray shading.
X-ray emission from the third, coldest, plasma zone 3 is connected with the appearance of the Kα line in the spectrum. This line corresponds to transitions in atoms and ions of iron with an ionization state not higher than 3. These states are generally caused by hot electrons that have heated the target to a temperature of a few tens of eV. The shape of the Kα line and its position in the spectrum can be used to determine the temperature of the plasma electrons in zone 3.35 In case C, this temperature is Te3 ∼ 80 (±30) eV.
If we sum the emission spectra from all three plasma zones, we obtain a good description of the experimentally observed spectrum in the entire observed spectral range. This can be seen in Fig. 10, which shows the X-ray spectrum measured in case C and the model spectrum calculated for the parameters listed in Table III.

Figure 10.Comparison of the results of kinetic modeling obtained under the assumption of X-ray emission from three different plasma zones and the experimentally measured spectrum for a laser intensity on target Ilt = 5 × 1021 W/cm2. (a) Sum of the modeled spectra for the three main plasma zones compared with experimental case C. (b) Contributions of each plasma zone in the sum of the modeled spectra.

Table 3. Plasma parameters used in the calculations to describe the X-ray spectrum measured in experimental case C.
Table 3. Plasma parameters used in the calculations to describe the X-ray spectrum measured in experimental case C.
Parameter | Zone 1 | Zone 2 | Zone 3 |
---|
Ni (cm−3) | 8 × 1022 | 8 × 1022 | 8 × 1022 | Ne (cm−3) | 1.94 × 1024 | 1.8 × 1024 | 5.5 × 1023 | Te (eV) | 2500 (±200) | 900 (±100) | 80 (±30) | Thot (keV) | 0.1%, 10 keV | 0.1%, 10 keV | 0.1%, 10 keV | Electron energy density (J/cm3) | 1.2 × 109 | 0.4 × 109 | 1 × 107 |
|
As can be seen from Fig. 10, to describe the X-ray emission of the plasma generated by irradiating thin SUS foils with femtosecond laser pulses with ultrarelativistic intensities (Ilt ∼ 5 × 1021 W/cm2), it is necessary to consider at least three plasma regions. Of course, the consideration of three zones is a rather rough approximation that only describes the most intense parts of the observed spectrum well, but this is already sufficient for many cases.
If a more precise description of low-intensity spectral features in the 1.88–1.91 Å range is required, we can divide zone 2 into two subzones 2.1 and 2.2 with different electron temperatures and solid densities, as has been done in Ref. 22. For the spectrum in Fig. 2(c), which includes subzone 2.2 with Te2.2 ∼ 400 eV, a more accurate agreement between modeling and experiment can be achieved in the wavelength range 1.88–1.91 Å (see the green curve in Fig. 10). Division of zone 3 into two subzones 3.1 and 3.2 with Te3.1 ∼ 170 eV and Te3.2 ∼ 80 eV also improves the agreement between modeling and experiment in the spectral region near the Kα line.
C. Discussion
In experimental high-energy-density physics, it is important not only to create the plasma in question, but also to measure (or at least estimate) its parameters. In the case of relativistic laser plasmas, the dimensions of the object are very small, owing to the necessity of using high laser intensities and consequently extremely small sizes of the focal laser spot. As a result, the dimensions of the object are smaller than or of the order of the spatial resolution of the X-ray spectrometer, and only spatially integrated emission spectra are available for diagnosis. At the same time, the relativistic plasma is very inhomogeneous despite its small size, and it is usually impossible to describe its emission spectrum, which is caused by transitions in ions of different multiplicity, by single values of density and temperature. A possible solution to this problem is to divide the entire plasma into several zones with different temperatures and densities. To prevent the number of free parameters from becoming too large, it is advisable to include as few zones as possible in the consideration, and the subdivision into zones should be physically justified by the predominance of certain physical processes in each of them.
In this paper, using the spectra of a Fe plasma as an example, we have shown that it is possible to estimate the probable value of the plasma density and electron temperature corresponding to the maximum ionization states for each ion (see Fig. 4 and Table II). We have demonstrated that it is sufficient to consider three plasma zones to describe the emission spectrum of H-, He-, and Li-like Fe ions. The good agreement between experimental measurements and the modeled X-ray spectrum gives confidence in the assumption of different plasma zones, allowing us to determine the electron temperature in each region. In particular, we have observed that the target region heated by the main laser pulse, zone 1, reaches Te1 ∼ 2500 eV, while the peripheral target region around the focal point is heated to Te2 ∼ 900 eV and the region far from the laser axis has a bulk electron temperature Te3 ∼ 80 eV. An intriguing and unexpected observation is the absence of a plasma zone at the critical electron density Ne,cr. This indicates that the laser contrast was high enough to prevent the formation of a pre-plasma and allow direct absorption of the laser energy in the solid-density target. This enabled the generation of record high energy densities of 1.2 GJ/cm3 in the focal spot region and ∼10 MJ/cm3 in the peripheral target region.
IV. CONCLUSION
We have demonstrated the generation of matter with an energy density of ∼1.2 GJ/cm3 (which corresponds to about 12 Gbar pressure) in a laboratory experiment by the interaction of ultrarelativistic high-contrast femtosecond laser pulses with flat SUS foils. A further possibility to increase the energy density will be available with the commissioning of new laser systems with intensities of more than 1023 W/cm2.12,15,36 High-resolution X-ray spectroscopic methods are still effective to study the extreme states of matter produced in new-generation laser systems. However, a deeper understanding of the underlying physical processes requires more complex approaches to the modeling of spectra of multiply charged ions, including nonstationary cases.
ACKNOWLEDGMENTS
Acknowledgment. X-ray measurements, X-ray spectral simulations, and analysis were supported by the Ministry of Science and Higher Education of the Russian Federation (State Assignment No. 075-00270-24-00). Calculations were carried out on the computational resources of the JSCC RAS. The work was carried out within the framework of Program 10 “Experimental laboratory astrophysics and geophysics, NCPM.”