Analysis of microscopic properties of radiative shock experiments performed at the Orion laser facility
1 Introduction
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 Theoretical model
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.
2.1 Atomic data
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 .
2.2 Determination of plasma atomic-level populations
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
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
The set of rates equations in the CR model is coupled to the radiative transfer equation
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 (
2.3 Calculation of radiative properties
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.
3 Description of the experimental setup and radiation-hydrodynamics simulations
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
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
Fig. 1. (a) Simulated mass density and (b) simulated electron temperature at 16 ns. (c) Experimental X-ray backlighting at 25 ns. The dashed lines mark the position of the diagnostic window on the gas-cell targets.
4 Results
4.1 Analysis of the radiative precursor
Figure
Fig. 2. Electron temperature (dashed lines) and mass density profiles of one of the radiative shocks as a function of time and position obtained with the 2D radiative-hydrodynamic simulation.
Fig. 3. Axial electron temperature (orange) and mass density (blue) profiles of one of the radiative shocks at 8 ns and 16 ns, deduced from the 2D radiation-hydrodynamics simulations. An electron density profile is also represented in green at 16 ns.
Figure
Fig. 4. Charge state distribution (CSD) as a function of the electron temperature at the mass density in the radiative precursor ().
The CSDs were also calculated assuming LTE, that is, through the SB equations (Equation (
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
Figure
According to Figure
Fig. 5. Division of layers of the radiative precursor at . Layer 1 is located closest to shock front, and layer 4 furthest.
Fig. 6. Specific intensities of the radiation emitted by different layers in the radiative precursor.
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 (
Fig. 7. Monochromatic opacities of the radiative precursor at four characteristic temperatures (, , and 20 eV).
In Figure
The analysis of the intensities and (Figure
Figure
Figure
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.
4.2 Analysis of the post-shock medium
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
The effect of the increase in mass density on the average ionization is evident when comparing Figures
Fig. 8. (a) Charge state distributions and (b) their monochromatic emissivities for two plasma conditions of the post-shock medium at 8 ns.
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
Figure
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
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
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]
5 Conclusions
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.
[1] 1.Y. B.Zel’dovich and Y. P.Raizer,Physics of Shock Waves and High-Temperature Hydrodynamic Phenomena (Dover Publications Inc., 2002).
[2] 2.C. F.McKee and B. T.Draine,Science252, 397 (1991).10.1126/science.252.5004.397
[3] 3.J. M.Shull and C. F.McKee,Astrophys. J.227, 131 (1979).10.1086/156712
[6] 6.A. C.Raga,G.Mellema,S. J.Arthur,L.Binette,P.Ferruit, and W.Steffen,Rev. Mex. Astron. Astr.35, 123 (1999).
[7] 7.E.Huguet and J. P. J.Lafon,Astron. Astrophys.324, 1046 (1997).
[8] 8.J. H.Hunter, Jr,Astrophys. J.161, 451 (1970).10.1086/150551
[10] 10.D.Ryutov,R. P.Drake, and B.Remington,Astrophys. J. Suppl. Ser.127, 465 (2000).10.1086/313320
[11] 11.J.Castor,Astrophys. Space Sci.307, 207 (2007).10.1007/s10509-006-9236-3
[21] 21.S.Bouquet,C.Stéhle,M.Koenig,J. P.Chieze,A.Benuzzi-Mounaix,D.Batani,S.Leygnac,X.Fleury,H.Merdji,C.Michaut,F.Thais,N.Grandjouan,T.Hall,E.Henry,V.Malka, and J. P. J.Lafon,Phys. Rev. Lett.22, 056501 (2004).
[33] 33.D. L.Youngs,Physica D12, 32 (1984).10.1016/0167-2789(84)90512-8
[34] 34.M. F.Gu,Can. J. Phys.86, 675 (2008).10.1139/p07-197
[36] 36.O.Peyrusse,C.Bauche-Arnoult, and J.Bauche,Phys. Plasmas12, 063302 (2005).10.1063/1.1931109
[38] 38.J. C.Stewart and K. D.Pyatt,Astrophys. J.144, 1203 (1966).10.1086/148714
[41] 41.R.Rodríguez,R.Florido,J. M.Gil,J. G.Rubiano,P.Martel, and E.Mínguez,Laser Part. Beams26, 433 (2008).
[42] 42.M. S.Dimitrijevic and N.Konjevin,Astron. Astrophys.172, 345 (1987).
[43] 43.H. A.Kramers,Philos. Mag.46, 836 (1923).10.1080/14786442308565244
[44] 44.D.Salzmann,Atomic Physics in Hot Plasmas (Oxford University Press, 1998).
[46] 46.http://www.awe.co.uk/app/uploads/2014/10/AWE-Orion-Fact-Sheet-No1.pdf.
[47] 47.G. B.Field,Astrophys. J.142, 531 (1965).10.1086/148317
[48] 48.Yu. A.Shchekinov,Sov. Astron.22, 182 (1978).
Article Outline
R. Rodríguez, G. Espinosa, J. M. Gil, F. Suzuki-Vidal, T. Clayson, C. Stehlé, P. Graham. Analysis of microscopic properties of radiative shock experiments performed at the Orion laser facility[J]. High Power Laser Science and Engineering, 2018, 6(2): 02000e36.