Shock waves are some of the most interesting and prevalent phenomena in astrophysics. They are ubiquitous throughout the universe and their role in the transport of energy into the interstellar medium is fundamental^[1]. Hence, the understanding of the structure of the interstellar medium requires knowledge of the dynamics and evolution of shock waves^[2]. Radiative shocks occur when shocked matter becomes hot enough that radiative transport modifies the shock structure and its dynamics^[3]. In some low density cases, the heated post-shock medium is ionized and emits radiation which leads to radiative cooling. Radiation from the post-shock region can also heat and ionize the unshocked medium ahead of the shock giving rise to a radiative precursor. Radiative shocks are observed in accretion shocks^[4], supernova remnants^[5], bow shocks at the head of stellar jets^[6] and pulsating stars^[7]. Moreover, radiative shocks are also predicted to exhibit thermal cooling instabilities^[8], which occur due to an imbalance between heating and cooling rates, which is a topic of high interest in astrophysics since they could be related to the formation of several astrophysical objects, for example star formation^[2]. For these reasons, the study of radiative shocks is a research area of great interest at present.

The possibility of scaling the magnetohydrodynamics equations between astrophysics and laboratory flows^[9–13] and the development of high-energy density (HED) facilities^[14], such as high-power lasers and fast magnetic pinch generators, have allowed the design of HED laboratory astrophysics experiments. High-power lasers have been commonly used to perform experiments to produce radiative shocks in gases with different approaches: irradiating a pin or a foil within a moderate to high Z background gas^[15–17], directly depositing the laser energy into clustered gases^[18–20] or driving a solid piston into a gas cell or tube^[21–25].

New experiments have been recently performed to investigate the collision and interaction of two counter-propagating piston-driven radiative shocks^[26–29] and their radiative precursors. These experiments open new applications in the area of laboratory astrophysics, such as the study of reflected and transmitted radiative shock and the interaction of a radiative shock with a denser medium in astrophysics scenarios, and were conducted at the Prague Asterix Laser System (PALS)^[30] and Orion high-power laser^[31] facilities. The macroscopic properties of the experiment were simulated with NYM^[32] and PETRA^[33] radiation-hydrodynamics codes. On the other hand, either for the radiation-hydrodynamics simulations or to interpret experimental data, such as the absorption or emission spectra, plasma microscopic properties, such as the composition and radiative properties of the plasma, must be obtained and analyzed since they provide useful and complementary information to radiation-hydrodynamics simulations.

The main objective of this work is to provide a description of microscopic properties of the plasmas generated in the experiments described above. We have therefore performed numerical simulations of the atomic and radiative properties of the post-shock and radiative precursor regions of radiative shocks in xenon. For this purpose, the plasma conditions (mass densities and electron temperatures) used in the microscopic simulations were provided by the radiation-hydrodynamics (rad-hydro) simulations performed with the NYM and PETRA codes. In particular, we have studied the thermodynamic regime of the plasma in both the post-shock and radiative precursor regions, the monochromatic opacities and emissivities, and the charge state distributions (CSDs) for times prior to the shocks collision as after this time the xenon mixes with the piston material (plastic and bromine). An analysis of the specific intensities of the radiation emitted by the plasma of both regions at early times was made, when the interaction between the shocks and the radiative precursors is weak and the spectra are due to a single shock. However, some of the results obtained could also be useful in analyzing the spectra at later times when the radiative precursors interact. The analysis of the intensities may also be helpful for future experiments in assessing whether an experimental capability (such as emission spectroscopy) would be useful. Finally, the radiation-hydrodynamics simulations predicted hydrodynamic instabilities that could be due to the strong cooling of the radiative shocks. For this reason, a theoretical analysis of the possibility of the onset of thermal instabilities in these experiments has been carried out.

This paper is structured as follows. Section 2 is devoted to presenting the theoretical model used in the computation of the microscopic properties. In Section 3, a brief description of the experimental setup and the radiation-hydrodynamics simulations are shown. Section 4 presents the results of the simulations of the radiative and spectroscopic properties of the post-shock and the radiative precursor regions together with the analysis of the thermal instabilities in the cooling layer in the post-shock medium. Finally, conclusions are drawn in Section 5.

The calculation of plasma radiative and spectroscopic properties requires the use of atomic data, such as energy levels, oscillator strengths, cross-sections and atomic-level populations. The following section briefly describes the models used in this work to calculate these values.

The oscillator strengths, energy levels and photoionization cross-sections were calculated using FAC code^[34], in which a fully relativistic approach based on the Dirac equation was used. The photoionization cross-sections were obtained in the context of the relativistic distorted wave approximation. In this work, the atomic data were calculated by the relativistic detailed configuration accounting (RDCA) approach. The unresolved transition array (UTA)^[35] formalism was used for the bound–bound transitions. Therefore, the transition energies include the UTA shift and the width of each transition is considered for the line profile. Furthermore, a correction to the oscillator strengths due to the configuration interaction within the same nonrelativistic configurations was included.

The selection of the atomic configurations for population kinetics simulations is a key factor but still an open question, overall for complex elements like xenon with a large number of electrons. In this work, we have included configurations with energies up to three times the ionization potential. This choice should be adequate for accurate modeling of thermal plasmas^{[36, 37]}.

According to this criterion, the following configurations were included: (1) ground configuration; (2) single excited configuration from the valence shell, , to shells with ; (3) doubly excited configurations from the valence shell to shells with ; and (4) single excited configurations from shell to shells with .

For high densities, when the plasma reaches local thermodynamic equilibrium (LTE), the population of the ionization stages, , is obtained from the Saha equation given by

 (1)

For arbitrary densities, the atomic-level populations can be determined from the solution of a system of collisional-radiative (CR) rate equations. This set of kinetic equations is given by

 (2)

The set of rates equations in the CR model is coupled to the radiative transfer equation

 (3)

The CR model as well as the Saha–Boltzmann (SB) equations used in this work is implemented in the MIXKIP code^[39], a code developed to calculate plasma atomic-level populations of mono and multicomponent optically thin and thick plasmas in time-dependent and steady-state situations. The atomic processes included in the CR model implemented in MIXKIP are collisional ionization, three-body recombination, spontaneous decay, collisional excitation and deexcitation, radiative recombination, autoionization and electron capture. Furthermore, in order to take into account the effect of external radiation fields in the calculation of atomic-level populations, the radiative driven processes of photoexcitation, photodeexcitation and photoionization are included in the CR model. A more detailed description of the expressions used for the rates of the atomic processes can be found in Ref. [39]. Plasma self-absorption (i.e., opacity effects) is modeled in MIXKIP in an approximate way using the escape factor formalism for the bound–bound opacity^[40].

Once the rate or SB equations (Equations (1) and (2), respectively) are solved, both the plasma average ionization and the plasma CSD, defined as the set of the population densities of the ions present in the plasma for a given condition of density and temperature, can be obtained.

To calculate the radiative properties, we have used the RAPCAL code^[41]. The monochromatic emissivity and opacity ( and , respectively), include bound–bound, bound-free and free–free contributions. In this work, the complete redistribution hypothesis was assumed for the line profile in the bound–bound transitions. This includes natural, Doppler, UTA and electron impact^[42] broadenings. The line shape function is applied with Voigt profiles incorporating all of these broadenings. For the bound-free contribution, we used the photoionization cross-section calculated with the FAC code in the relativistic distorted wave approach. Finally, for the free–free contribution, the Kramers semi-classical expression for the inverse bremsstrahlung^[43] cross-section was used.

The experiments^[27] were conducted on the Orion laser and the setup consisted of gas cells filled with xenon at bar . Plastic discs were attached to opposite octagonal faces, separated by 4 mm, which acted as pistons. Shocks were driven by focusing four laser beams (each beam , , 1 ns pulse duration) onto each piston (laser intensity ). These experiments were diagnosed with point projection X-ray backlighting (XRBL) and laser optical interferometry. XRBL provided information about the post-shock region. For the radiative precursor, laser interferometry was used to measure the free electron density along the probe beam path. The large xenon volume used in the experiments resulted in shocks with quasi-spherical geometry. A collision between the two shocks occurred after followed by the formation of reverse shocks. The XRBL images indicated a shock velocity of around and a reverse shock velocity of around 30 km/s. Figure 1 shows an example of results from rad-hydro simulations and experiments, with XRBL results, shown in Figure 1(c) at 25 ns.

2D rad-hydro simulations were performed using then NYM and PETRA codes. NYM is a Lagrangian code with multigroup implicit Monte Carlo X-ray transport and full laser-interaction physics. This code was, therefore, used to model the laser–piston interaction. The simulations provided by the NYM code were linked to the Eulerian code PETRA which used multigroup X-ray diffusion to study the late time plasma behavior (after ). The collision of the counter-propagating shocks was simulated assuming a fully reflective boundary (for hydrodynamics and radiation) at the center of the diagnostic window (at a distance of 2 mm from the position of the pistons) both for the plasma flow and for the radiation. In general, the simulations accurately reproduced the overall shock dynamics^[33]. As an example, Figures 1(a) and 1(b) show 2D maps of density of matter and electron temperatures of one of the shocks at 16 ns, obtained from the simulations. In Figure 2, we also show the axial profiles of the mass densities and electron temperatures for times analyzed in this work.

Figure 3 shows the axial profiles of the electron temperatures and mass densities of one of the radiative shocks (propagating from left to right) at and 16 ns obtained from 2D rad-hydro simulations. For the latter, the profile of the electron densities is also presented. The shock front is located at the highest temperature and ahead of the shock front a radiative precursor is observed. Simulations suggest that temperatures in the radiative precursor range from to 20 eV at 8 ns. As time progresses, the lowest temperatures in the radiative precursor increase (maybe due to the proximity of the other shock and the overlapping of the two radiative precursors), as shown in Figure 3 right. At 8 ns the whole range of temperatures mentioned above is present. At this early time, we can consider the whole system as a single shock because the interaction between the two shocks is not predominant as the radiative precursors are not yet overlapping.

Figure 3 shows that the characteristic time of the evolution of plasma conditions is of the order of . However, the characteristic time of the dominant atomic processes (electron–ion collisions) in the plasma, , is given by^[44]

 (4)

The CSDs were also calculated assuming LTE, that is, through the SB equations (Equation (1)), obtaining very small differences in the average ionization between LTE and NLTE calculations. These differences could be greater for the CSD since this is a less average quantity than the average ionization. However, the differences in the fractional abundances only begin to emerge for ions at temperatures around 15 eV (although these differences are lower than 5%) and become more noticeable at 18 eV. At this temperature, the LTE model predicts that the most abundant ion is with an average ionization of while the NLTE model predicts the most abundant ion as with an average ionization of . Therefore, the LTE approach is accurate enough for the calculation of the average ionization and CSDs for temperatures lower than 18 eV, with slight differences at higher temperatures.

The range of electron densities obtained from the CRSS simulations in the radiative precursor varies from at 2 eV to at 20 eV, obtaining similar results to the ones provided by the radiation-hydrodynamics simulations, which use SESAME tables for the opacities and equations of state. For example, the electron densities obtained with the atomic kinetic model at and 20 eV were , and , respectively, whereas those obtained with the macroscopic simulation at 8 ns were , and , respectively, with relative differences lower than between both simulations. Since the temperatures used in the CRSS simulations are those provided by the rad-hydro simulations, the agreement between the electron densities obtained indicates that kinetic atomic models implemented in both simulations should provide similar values of the average ionization.

The electron densities of a region of the radiative precursor were experimentally obtained at 18 ns by means of laser interferometry^[27]. The values obtained ranged from to . As the mass density in the precursor does not change significantly from the initial density of the xenon gas fill, the average ionization can be determined. Using the mass density as an input, the CRSS model was used to determine the electron temperatures that provide similar average ionizations and electron densities to the experimental ones. This yielded temperatures from 4.20 to 6.50 eV for the radiative precursor. For that region of the radiative precursor, the hydrodynamic simulations predicted values of the electron densities from to and electron temperatures from to 7.17 eV. These 2D simulations slightly overestimate the values of the electron density with respect to the experimental ones and, therefore, the electron temperatures as well^[27].

Although the electron densities reached in the radiative precursor are not particularly high, the surrounding plasma, modeled by the inclusion of the CL in the population calculation, introduces some differences in the CSDs since the electron temperatures were also low. Thus, for example, at 12 eV, the most relevant ions obtained with the simulation including the CL were , whereas the ones obtained with the simulation for the isolated situation were . This could lead to noticeable changes in the calculation of the spectra. Therefore, the effect of the plasma surrounding must be considered in the atomic kinetic calculations for these ranges of electron temperatures and densities.

In the experiments^{[27, 28]}, optical and X-ray diagnostics were fielded to study the plasma in the transverse direction to the direction of propagation of the radiative shocks. Therefore, radiation emitted from either the shocked material or from the radiative precursor will have to pass through regions of the radiative precursor with different temperatures. This material can absorb some of that radiation, before reaching the detectors. For this reason, the monochromatic opacities of the radiative precursor were studied. The monochromatic opacities in the radiative precursor were calculated for four different characteristic temperatures (4, 10, 16 and 20 eV) and are presented in Figures 4(a) and 4(b).

Figure 4(a) shows that the opacities for temperatures of and 10 eV are quite different. Analysis of the CSDs (Figure 4) shows that at 4 eV there are predominantly four ion species in the plasma, . The CSD at 10 eV shows that the ions species in the plasma are although the abundance of and ions are minor compared to the other ions. Therefore, none of the relevant ions at 4 eV are at 10 eV. This fact explains the significant differences between both spectra. The bound–bound contribution at 4 eV is practically contained in photon energies between 0–30 eV which are associated with single electron transitions from and subshells. For photons with higher energies the main contribution is bound-free. At a temperature of 10 eV an absorption structure in the photon energy range 0–30 eV is present and primarily due to electron transitions in the same shell of ions. Another two significant absorption features are present for a temperature of 10 eV. The first one located around a photon energy of 60 eV is primarily due to electron transitions from to higher shells in ion. For the other two ions present, the transitions involved are those between and subshells. The second absorption feature is around 90 eV. In this case, the most relevant transitions are the ones between the and and between the and subshells. Therefore, from this analysis it is clear that the bound–bound absorption at electron temperatures of and 10 eV is basically restricted to the UV and XUV ranges.

According to Figure 4, the CSDs for electron temperatures of and 20 eV are and ions, respectively. Because of the similarity between them, the absorption structures should be quite similar in the monochromatic opacities at both temperatures, as Figure 4(b) shows. The main differences are in the height of the peaks due to the differences in the ions populations. There are other discrepancies; for example, the peak at around 60–80 eV, due to ion, is present at 16 eV but not at 20 eV. For these two temperatures significant absorption structures in the range 0–30 eV were not obtained on account of the ionization degrees as these temperatures are larger than at and 10 eV. We can also observe that the monochromatic opacities for the two highest temperatures have richer line spectra than the ones for the other two lower temperatures due to a stronger involvement of the shell. The figure shows that the monochromatic opacities have a strong absorption feature in the energy range 70–85 eV. For the ion, the main transitions involved at these energies are the ones from configurations to configurations. However, similar transitions from configuration for ion are shifted to higher photon energies (around 90–100 eV) and, therefore, this ion does not have a relevant contribution to that absorption feature at photon energies around 70–85 eV in the opacity at a temperature 20 eV. Line transitions are detected in the energy ranges 100–160 and 190–250 eV in which transitions to configurations with higher principal quantum number and also from and subshells are involved. Therefore, at the temperatures of and 20 eV there is noticeable absorption of XUV.

As previously discussed, spectroscopic diagnostics is a useful tool for obtaining information about the plasma conditions. Typically, shell spectrum is commonly used to derive plasma conditions^[45] and due to their more complicated structure, shell spectrum is less useful. Under the conditions analyzed in this work, the ionization degree of xenon is low, both in the radiative precursor and in the post-shock regions, and the ions present in the plasma are those from to , which means that and shells participate in the spectra. These shells contribute considerably more complex spectra than shell, and so only permit estimations of ranges of electron temperatures in the plasma. The method consists in analyzing the presence of the contribution of certain ions to the spectra, taking into account the temperature sensitivity of xenon in this range of plasma conditions. Thus, the specific intensities emitted by the radiative precursor were calculated at different temperatures. Assuming stationary situation for the radiation in Equation (3), the intensities were computed along the beam whose propagation direction is transverse to that of the radiative shocks in the gas cell. The intensity emitted by a region of the radiative precursor must travel through regions with different electron densities and temperatures before the radiation is recorded by the detector. Since the shocks obtained in the experiment have quasi-spherical symmetry (see Figures 1(b) and 1(c)), we have assumed, in order to estimate radiation transport through the radiative precursor, that the electron temperature profiles in the transverse direction are the same as those in the direction of the propagation of the shock represented in Figure 5. This assumption will only be valid until the radiative precursors of the two counter-streaming radiative shocks overlap. After this point, this symmetry is not a true representation of the plasma conditions in the radiative precursor. Therefore, we have restricted ourselves to estimate and analyze the specific intensity due to different regions of the radiative precursor at 8 ns. With that purpose we have divided the radiative precursor into four homogeneous layers, each characterized by an average temperature (see Figure 5), considering the different ranges of temperatures that can be found in this region. Thus, the radiation emitted by a layer must be transported through the layers that lie ahead and the specific intensity of that radiation along the ray is given by

 (5)

In Figure 6, the specific intensities for the three first layers, calculated according to Equation (5), are displayed (black line). These represent the intensity of outgoing radiation from the different layers of the radiative precursor. The specific intensity of the radiation emitted by layer (indicated in Figure 6 with orange line) is given by

 (6)

 (7)

The analysis of the intensities and (Figure 6(a)) provides information about the transmission and emission of Layer 4, since . From the analysis of we can detect that Layer strongly absorbs in the photon energy range 0–40 eV. This absorption is rather important between 15–20 eV as we can see from the purple curve in Figure 6(a). This agrees with the results obtained for the monochromatic opacities (see Figure 7(a)) and with the large thickness of Layer 4. Consequently, the specific intensity for photons with energies between 0–40 eV is mainly due to the radiation emitted by Layer 4. On the other hand, for energies higher than 60 eV the emission of radiation from Layer 4 is negligible with respect to the one from Layer 3 and that part of the spectra is caused by the radiation emitted by Layer 3 after crossing Layer 4, and thus the coincidence between and .

Figure 6(b) shows that self-absorption of Layer 2 is quite significant for photon energies between 60 eV and 70 eV and around 90 eV. This is expected from the important absorption peaks up to in the monochromatic opacity observed in Figure 7(b). The structure of the spectrum for energies up to 70 eV is quite similar to the one obtained for Layer 3, since the radiation emitted in Layer 2 travels through Layers 3 and 4. So, the detection of the structures in the intensity spectra between 0–40 and 50–70 eV could indicate regions of plasma with temperatures lower than 5 eV and around 10 eV, respectively. In the spectrum of Layer 2 new strong features are also detected in the range 70–100 eV and to a lesser extent for energies from 110 to 160 eV. For higher energies the radiation from Layer 3 dominates.

Figure 6(c) reveals the main differences between the emission spectra from Layers 2 and 1 which are, basically: the increase in the structure located in the middle of the two relevant structures at 60–70 and 85–95 eV, respectively, observed in Figure 6(b); a broadening of the peak in that latter range of energy; and also, an increase in the contribution in the range of energies higher than 110 eV. However, the spectra from both layers are quite similar as expected from the analysis of their opacities. In this range of plasma temperatures (16–20 eV) there are not noticeable changes in the CSDs (see Figure 4). Consequently, from a spectroscopic point of view, it is more difficult to distinguish between spectra of different temperatures in the range 16–20 eV than for ranges of lower temperatures.

In order to analyze the effect of the resolution of a hypothetical spectrometer on the calculated spectra, simulations of the theoretical spectra convolved with the resolution of the XUV grating spectrometer of Orion^[46] were performed, which is able to make spectrally resolved measurements of radiation in the range of 1–40 nm with a spectral resolution . Due to this high resolution, and that the xenon ions involved in these experiments have a large number of bound electrons (which produce unresolved structures in the spectra), the effect of the spectrometer on the theoretical spectra was not noticeable.

Rad-hydro simulations indicate that the mass densities in the post-shock region ranged from (with electron densities around ). The maximum electron temperatures are higher than 25 eV at 8 and 12 ns whereas they are lower at later times. A peak electron temperature of was reached at 8 ns and that time was selected for the analysis performed in this section. The CSDs for plasma conditions in the post-shock medium of the maximum and minimum temperatures, respectively, are plotted in Figure 8(a). Although temperatures larger than 20 eV are reached in the post-shock medium due to the increase in mass density (with respect to the radiative precursor) the plasma can be considered to be in LTE.

The effect of the increase in mass density on the average ionization is evident when comparing Figures 3 and 6(a). The most abundant ions at 18 eV in the radiative precursor are , while for the post-shock medium they are (i.e., the same ions that dominate at 10 eV in the radiative precursor) because of the increase in recombination. Another effect of the increase in mass density is a widening of the line transitions in the spectra due to the increase of the collisional broadening which leads to a large line overlapping. This is illustrated in Figure 8(b), where the monochromatic emissivities of the post-shock medium for two plasma conditions are displayed. As the figure shows, there are no detailed lines such as those observed in the monochromatic opacities of Figures 4(a) and 4(b). The main emission features of the shocked medium at 8 ns are for photon energies lower than 120 eV.

The specific intensity of the radiation emitted by the post-shock medium at 8 ns was also analyzed. The width of the shocked medium at this time is relatively small, , and so, the absorption in this region is smaller than seen in the radiative precursor. To calculate the intensity of radiation emitted by the post-shock region this was divided in three layers of of average electron temperatures and 18.22 eV and mass densities of, and , respectively. They correspond to the largest and lowest temperatures in the post-shock medium and also an intermediate one. The outgoing radiation was then transported through the radiative precursor assuming its structure is similar to that shown in Figure 5. This allows the results obtained in the discussion of the intensity of the radiative precursor to be used. Figure 9 shows the radiation emitted by the post-shock medium, denoted as and calculated as in Equation (6), the intensity from the post-shock after passing through Layer 1 of the radiative precursor, (obtained as in Equation (7)). The figure also shows the specific intensity of radiation emitted by both the post-shock medium and the radiative precursor, that would be the outgoing radiation from the plasma, given by Equation (4).

Figure 9 shows that the spectrum of radiation emitted by the post-shock medium, , consists of two peaks at photon energies around and 100 eV that have been broadened due to high mass densities. This has resulted in lines overlapping leading to broadened peaks and a reduction in the depths of the valleys. The absorption of Layer 1 of the radiative precursor (at a temperature of 20 eV) shows peaks at photon energies of around 75 eV and 90–105 eV (see in the Figure), which agrees with the stronger peaks detected in its monochromatic opacity (see Figure 7(b)).

With respect to the total specific intensity, the absorption in the photon energy range 0–60 eV is due to Layers 3 and 4 of the radiative precursor. For photon energies higher than 110 eV the spectrum has a larger contribution of the post-shock medium than that of the radiative precursor and there are not detailed lines. This lack of detailed lines provides information about the range of densities that can be found in the post-shock medium. Similar emission features are observed in the total intensity, in the energy range 60–110 eV, to those found in the intensity of Layer 1 in the radiative precursor, although more broadened. This is due to CSDs in the post-shock medium being quite similar to the radiative precursor layers at and 16 eV. This is because the difference in electron temperatures between the shock and the radiative precursor is not too large and recombination occurs in the post-shock plasma.

For times later than 16 ns the hydrodynamic simulations show a double peak in post-shock medium. Furthermore, the simulations also show the shock front as a rippled layer (see Figures 1(b) and 1(c)). Both phenomena reflect the spatial variations in mass density due to the growth of hydrodynamic instabilities^[28]. Due to the strongly radiative nature of these shocks, this could be thermal instabilities. These types of instabilities are due to an imbalanced heating and cooling rates and are expected to occur when the power of the heating source in the thermal energy equation is not adequate to ensure a steady temperature resulting in cooling becoming the most important process^{[47, 48]}. The criteria for their onset for non-stationary media were established by Shchenikov^[47] (a detailed explanation can be found in that reference) and these criteria were employed in this work.

Thermal instabilities can be classified by comparing the length scale of the initial seeding perturbation and a characteristic scale of the medium, which is the sound crossing length, , given as

 (8)

 (9)

 (10)

For the range of plasma temperatures in the post-shock medium at times later than 18 ns (22–15 eV) the sound crossing lengths were calculated to be . For perturbations whose lengths are comparable to the sound crossing length, the pressure will be constant and the thermal instabilities are classified as isobaric^[48]. This kind of instability produces local density increase that could be the origin of the ones predicted in the hydrodynamic simulations. The criterion for this instability is^[48]

 (11)

In this work we have analyzed the radiative and spectroscopic properties of both the radiative precursor and the post-shock medium of radiative shocks driven into xenon by a piston ablated by the Orion high-power laser. The plasma conditions used for the microscopic simulations were extracted from radiation-hydrodynamics simulations. From the analysis of the thermodynamic regime, we have concluded that the xenon plasma, both in the radiative precursor and the post-shock regions, could be assumed to be in LTE. The CSDs were calculated for plasma conditions achieved in the experiments, yielding a maximum average ionization of . This relatively low ionization implies that xenon ions with a large number of bound electrons are present in the plasma which results in complex spectra. The effect of the surrounding plasma on the ion populations was studied, including the CL in the CRSS simulations. The electron densities achieved in the radiative precursor are not too large, but since the electron temperatures are relatively low, plasma effects should still be considered. These effects are more noticeable in the post-shock medium since the mass density is considerably larger than in the radiative precursor whereas the temperatures are similar in the two regions. Furthermore, the CSDs in the post-shock medium are quite similar to those obtained in the radiative precursor, even though the temperatures in the latter can be lower, due to the increased recombination in the post-shock medium. The monochromatic opacities for several characteristic temperatures found in the radiative precursor were analyzed since the outgoing radiation of the plasma can be transported through this region before reaching a detector. This analysis allowed the plasma absorption to be characterized in terms of temperature, which helps to study the plasma spectra.

The specific intensities of radiation emitted by the radiative precursor and the post-shock medium were also analyzed. For the range of plasma temperatures reached in these experiments, the absorption and emission spectra are in the UV and XUV ranges of photon energies. For the analysis of the intensity, the analysis of the CSDs and the monochromatic opacities and emissivities carried out in this work were very helpful. We have determined that some structures detected in the spectra allow the ranges of plasma temperatures to be predicted, which are associated with ions present in the plasma. However, for the range of electron temperatures 15–20 eV, the CSDs are similar and thus it is difficult to distinguish contributions from different temperatures in that range of temperatures. Therefore, from the spectra we could only conclude that there are regions in the plasma with temperatures in the range 15–20 eV. On the other hand, due to the strong dependence of the first charge states of xenon on temperature, it is easier to distinguish contributions from different temperatures in the range 2–10 eV.

We have also analyzed the intensity of radiation emitted by the post-shock medium. Since the mass density is considerably higher than in the radiative precursor, lines are collisionally broadened, resulting in lines overlapping, and individual lines cannot be identified. This could allow us to estimate the plasma density in the post-shock medium.

Finally, a theoretical analysis of the isobaric thermal instability was made in order to explain the instabilities observed in rad-hydro simulations. Since the criterion of this instability is fulfilled, this could be responsible for the double peaks and ripples observed in the simulations. However, confirmation would require both simulations with radiation turned off and experimental measurements of the length of the seeding perturbations.