## Observation of Shock-Front Separation in Multi-Ion-Species Collisional Plasma Shocks

Tom Byvank, 1, a) Samuel J. Langendorf, 1, b) Carsten Thoma, 2 and Scott C. Hsu 1 1) Los Alamos National Laboratory, Los Alamos, NM 87545 2)

Voss Scientific, Albuquerque, NM 87108

(Dated: 19 March 2020)

We observe shock-front separation and species-dependent shock widths in multi-ion-species collisional plasma shocks, which are produced by obliquely merging plasma jets of a He/Ar mixture (97% He and 3% Ar by initial number density) on the Plasma Liner Experiment [S. C. Hsu et al., IEEE Trans. Plasma Sci. 46 , 1951 (2018)]. Visible plasma emission near the He-I 587.6nm and Ar-II 476.5-514.5-nm lines are simultaneously recorded by splitting a single visible image of the shock into two different fast-framing cameras with different narrow bandpass filters (589 ± 5 nm for observing the He-I line and 500 ± 25 nm for the Ar-II lines). For conditions in these experiments (pre-shock ion and electron densities ≈ 5 × 10 14 cm -3 , ion and electron temperatures of ≈ 2 2 eV, and relative plasma-merging speed of 22 km/s), the observation-. ally inferred magnitude of He/Ar shock-front separation and the shock widths themselves are &lt; 1 cm, which correspond to ∼ 50 post-shock thermal ionion mean free paths. These experimental lengths scales are in reasonable qualitative and quantitative agreement with results from 1D multi-fluid simulations using the chicago code. However, there are differences between the experimentally-inferred and simulation-predicted ionization states and line emission intensities, particularly in the post-shock region. Overall, the experimental and simulation results are consistent with theoretical predictions that the lighter He ions diffuse farther ahead within the overall shock front than the heavier Ar ions.

## I. INTRODUCTION

Supersonic flows generate shocks in astrophysics, aerodynamics, and high-energy-density (HED) plasma experiments. Compared to hydrodynamic shocks in neutral gases, collisional plasma shocks contain ion and electron species, arise due to Coulomb collisions, and are influenced by electromagnetic fields 1-4 . A plasma shock front with multiple ion species contains additional structure compared to a single ion plasma shock. Prior experiments, simulations, and theoretical work explored multi-ion-species effects in the context of inertial confinement fusion (ICF), for which species separation in the fusion fuel potentially leads to neutron yield degradation 5-18 . Interspecies ion separation and velocity separation have been experimentally observed 19-23 . Additional simulation and theoretical research on multiion-species plasmas examined how ion species diffusion causes species separation 24-30 . The present research reports direct observations of the spatial profile of a multi-ion-species shock, showing shock-front separation and species-dependent shock widths in collisional plasma shocks. The experimental length scale results of the shock structure agree with 1D multifluid simulations using the chicago code 31-33 . However, particularly in the post-shock region, there are differences between the experimentally-inferred and simulation-predicted ionization states and line emission intensities. These experimental and simulation results are both consistent with ion species diffusion theory. This experimental data can be used to validate and benchmark numerical simulations of plasma environments with multi-ion-

a) tbyvank@lanl.gov

b) samuel.langendorf@lanl.gov

species collisional plasma shocks, especially in HED, magneto-inertial-fusion (MIF), and ICF experiments.

The organization of this paper is as follows: Sec. II provides background on the ion diffusion theory to which we compare our results, Sec. III describes the experimental setup of colliding plasma jets to generate a shock, Sec. IV presents results from the experimentally observed shock profiles, Sec. V overviews the relevant length scales in our experiments, Sec. VI A discusses the agreement between the multi-fluid simulations and the experiments, Sec. VII gives estimates for the contributions from various diffusion mechanisms to describe our observations, and Sec. VIII highlights our conclusions.

## II. BACKGROUND ON ION DIFFUSION THEORY

Theory and simulations about interspecies ion diffusion predict that lighter ions diffuse farther ahead within a collisional plasma shock (closer to the pre-shock region) than heavier ions 11,12,14-18,24-28,30 . Following Kagan and Tang 24 , in the center-of-mass frame of a multiion-species fluid element, the diffusion flux ¯ f 1 of the lighter species 1 (in our case, He) equals the negative flux ¯ f 2 of the heavier species 2 (in our case, Ar),

$$\bar { f } _ { 1 } = \rho _ { 1 } \bar { v } _ { D 1 } = - \bar { f } _ { 2 } = - \rho _ { 2 } \bar { v } _ { D 2 },$$

with mass densities ρ 1 , ρ 2 , and diffusion velocities ¯ v D 1 , ¯ v D 2 . The diffusion flux is related to the species mass concentration time evolution (continuity equation),

$$\rho \frac { \partial c _ { m 1 } } { \partial t } + \rho \bar { u } \cdot \nabla c _ { m 1 } + \nabla \cdot \bar { f } _ { 1 } = 0,$$

with total mass density ρ , species 1 mass concentration c m 1 = ρ /ρ 1 , and bulk fluid velocity ¯. u The diffusion flux of the lighter species can be written as

$$\stackrel { \sim } { f } _ { 1 } & = - \rho D \nabla c _ { m 1 } - \rho D \left ( \frac { \kappa _ { p } } { P _ { i } } \nabla P _ { i } \right ) \\ & \quad - \rho D \left ( \frac { \kappa _ { T i } } { T _ { i } } \nabla T _ { i } + \frac { \kappa _ { T e } } { T _ { e } } \nabla T _ { e } \right ) \\ & \quad - \rho D \left ( \frac { e \kappa _ { e } } { T _ { i } } \nabla \Phi \right ), \\ \cdot \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,.$$

where the first term is the classical diffusion flux based on the mass concentration gradient, the second term is the barodiffusion flux based on the ion pressure gradient, the third term is the thermodiffusion flux based on the ion and electron temperature gradients, and the fourth term is the electrodiffusion flux based on the electric field (negative gradient of the electric potential). The new variables in Eq. (3) are the classical diffusion coefficient D , barodiffusion ratio κ p , total ion pressure P i , ion thermodiffusion ratio κ Ti , electron thermodiffusion ratio κ Te , ion temperature T i (approximating as equal for each species), electron temperature T e , electrodiffusion ratio κ e , and electric potential Φ.

Within a shock front, the gradients of pressure, temperature, and electric potential point in the direction from the pre-shock region toward the post-shock region . 3 Additionally, for the He/Ar mixture in our experimental shocks, the various diffusion fluxes have the same sign except for the relatively small electron thermodiffusion, as described in Sec. VII. We also assume no initial concentration separation. Therefore, in the center of mass frame, Eqs. (1) and (3) predict that the lighter species diffusion velocity ¯ v D 1 points in the direction from the post-shock region toward the pre-shock region, and the heavier species diffusion velocity ¯ v D 2 points in the opposite direction. In the present research, we directly observe the spatial profile of a plasma shock front containing a He/Ar mixture, offering the opportunity to validate models of multi-ion-species shock evolution based on the theory.

FIG. 1. (Left) View 1: Projection of experimental setup, including approximate diagnostic views. Plasma guns mounted on ports of the spherical vacuum chamber launch plasma jets toward the chamber center in the direction shown by the red dashed lines. Due to expansion of the plasma jets, the jet collision takes place well before the chamber center, in the region observed by the interferometry and spectroscopy. The plasma guns are not shown in this View 1 for diagram clarity. (Right) View 2: Corresponding schematic with a 90 ◦ out-of-page rotation of View 1, adding a mounted plasma gun. Not to scale.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00454v3-figure-0e51bbd4-00e9-4b88-8457-f0065ec25731.jpg)

## III. EXPERIMENTAL SETUP

In order to experimentally observe shock-front separation in multi-ion-species collision plasma shocks, we form the shocks by obliquely merging two supersonic plasma jets, and we image the different species within the shock profiles using distinct narrow bandpass filters. In this section, we describe this experimental setup.

We experimentally create multi-ion-species plasma shocks by colliding plasma jets generated from plasma guns in a 2.7-m-diameter spherical vacuum chamber 34-40 . Figure 1 depicts the experimental setup. Within the gun nozzle, a gas puff is pre-ionized into a plasma, and that plasma is accelerated by a 720 kA, 4.5 kV, 5 µ s rise-time pulse. Initially, the gas puff atomic concentration, c a , is 97% He and 3% Ar, corresponding to a mass concentration, c m , of 76.4% He and 23.6% Ar. As the plasma jets propagate toward the chamber center, they expand in vacuum at approximately the sound speed. Well before reaching the chamber center, the two plasma jets merge at a half-angle of 11.6 , with ◦ speeds of 55 ± 5 km/s corresponding to 11 km/s in the direction normal to the shock and, therefore, a relative normal speed between the jets of v rel = 22 km/s. Before the shock formation, individual jets have electron number densities n e ≈ 5 ± 1 × 10 14 cm -3 , electron and ion temperatures T e ≈ T i ≈ 2 2 . +0 2 . -0 5 . eV, and inferred average ionization states ¯ Z of 0 73 . +0 15 . -0 67 . for He and 1 15 . +0 23 . -0 16 . for Ar. These parameters give pre-shock Mach numbers M = v rel / γ [ ( ¯ ZT e + T i ) /m i ] 1 / 2 of 1.7 for He ions, 5.1 for Ar ions, and an average M = 1 8, . where we assume γ = 5 3 and a minimum ion / ¯ = 1. Z The jets have characteristic scale lengths of ∼ 10 cm. Since the interspecies mean free path is much smaller than this overall characteristic length scale (see Sec. V) and the thermal equilibration timescale is much smaller than the characteristic hydrodynamic timescale, the He and Ar should act as a single fluid with T e ≈ T He ≈ T Ar ≈ T i in the absence of gradients, i.e. pre-shock.

The plasma parameters are measured with diagnostics including a photodiode array to infer jet velocity via time-of-flight, time-resolved interferometry for density 41,42 , and emission spectroscopy for temperature and ionization states. The photodiode array consists of two photodiodes per gun spaced 2 cm apart along the jet propagation direction, and jet speed is determined by looking at relative timing of the plasma emission intensity between the two photodiodes. The heterodyne interferometer measures time-resolved lineintegrated electron density by observing relative phase shifts of a laser passing through the plasma compared to a reference path. We take the volumetric number density to be

TABLE I. Pre-shock and (simulation) post-shock parameters.

| Parameter                         | Pre-shock   | Post-shock   |
|-----------------------------------|-------------|--------------|
| Electron Temperature (eV)         | 2.2         | 2.4          |
| Electron Density ( × 10 14 cm - 3 | 5           | 21           |
| He Temperature (eV)               | 2.2         | 4.6          |
| Ar Temperature (eV)               | 2.2         | 5.5          |
| He Average Ionization ¯ Z         | 0.73        | 0.89         |
| Ar Average Ionization ¯ Z         | 1.15        | 1.47         |
| He Atomic Concentration c a 1     | 97%         | 95.8%        |
| Ar Atomic Concentration c a 2     | 3%          | 4.2%         |
| He Mass Concentration c m 1       | 76.4%       | 69.6%        |
| Ar Mass Concentration c m 2       | 23.6%       | 30.4%        |

equal to the line-integrated areal density divided by the path length, and the path length is estimated based on the plasma emission length observed by visible camera images (not shown) taken along a direction perpendicular to the interferometry. The density measurement uncertainty is based upon these estimates of the path length. Spectroscopy data for the optically thin jets and generated shocks are compared with prismspect non-localthermodynamic-equilibrium (non-LTE) atomic-modeling calculations 43 , as used in earlier research 34-40 . Figure 2 displays the experimental temperature and density diagnostic results. As listed in Table I, the pre-shock, experimentally inferred parameters are used to interpret the experimental observations and as input for simulations.

The shock profiles of the different ion species are imaged using a beam splitter to aim the plasma self-emission onto two intensified-charge-coupled-device (ICCD) cameras (10-ns exposure) with narrow bandpass filters 44,45 . Singly ionized Ar-II line emission (near the Ar-II 476.5-514.5-nm lines) is observed with a 500 ± 25-nm filter, and neutral He-I line emission (near the He-I 587.6-nm line) is observed with a 589 ± 5-nm filter 46 . A filter for singly ionized He-II line emission was not used due to the better ICCD camera sensitivity to visible compared to ultraviolet wavelengths and the presence of other stronger lines (Ar-II and/or He-I) near the visible He-II lines. Figure 3 illustrates how the distinct filters are sensitive to the different ion species, based on reference experiments using single-ion-species plasma jet merging. For the 100%-Ar jet merging, the bright region in the Ar-II filtered image correlates to an increase in singly ionized Ar emission. For the 100%-He jet collision, dark bands in the He-I filtered image correlate to a reduction in neutral He emission and a corresponding increase in He ionization. Thus, using these filters, we can separately image the shock profiles of the different ion species.

## IV. EXPERIMENTAL SHOCK PROFILE RESULTS

By imaging the experimental shocks, we can infer spatial separation of the ions within the shock profiles. Our spatially resolved shock profiles also allow us to observe different shock widths for the different species.

The top 4 rows in Fig. 4 show the time evolution of the merging of multi-ion-species plasma jets, observing the plasma self-emission using the Ar-II and He-I narrow bandpass filters. The time sequence is composed of images from different experimental shots, rather than images from multiple times during a single experimental shot. The jets originate from the top left and bottom left in the images, and they merge near the image center. The 0 µ s time is defined as the time when we observe the jets first starting to merge and form a shock structure. Utilizing the interferometry data, we measure shot-to-shot timing variations of &lt; 0 5 . µ s. As the time sequence progresses, the individual jets continue to propagate at 11 km/s inward (toward image horizontal center) and ∼ 54 km/s parallel to the shock front (toward the image right). We use the intensity minima and maxima to infer the shock position. The plasma emission intensity is a function of electron temperature T e and ion density n i , with larger temperatures and larger densities leading to higher ionization states

FIG. 2. (Top 3 plots) Comparison of pre-shock experimental emission spectroscopy (black) with prismspect atomic-modeling calculations (blue). (Row 4) Comparison of experimental timeresolved interferometry electron-density measurements of one individual jet (green), which represents the pre-shock density for a collision, the central density (post-shock) of two jets merging (blue), and the central density (post-shock) in the simulation of two jets merging (red). Times are shifted such that t = 0 is the time at which the jets are first starting to merge. (Bottom) Comparison of post-shock experimental emission spectroscopy (black) with prismspect atomic-modeling calculations (blue).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00454v3-figure-44e31394-8ee7-4df5-9400-69e1a5043711.jpg)

and correspondingly more emission. We expect the maximum intensity of the shock to be located very near the locations of maximum temperature and density, but the maximum temperature may not be located precisely at the maximum density position. The shock

FIG. 3. (Top) Plot of NIST strong emission lines (He-I, He-II, Ar-I, Ar-II) 46 and the spectral filter wavelengths. (Bottom 2 Rows) To establish a baseline, time-gated images (10-ns exposure) of single-ion-species plasma jet merging using narrow bandpass filters are shown. The 100%-Ar (left column) and 100%-He (right column) plasma shocks are much more prominently observed using the Ar-II and He-I filters, respectively. Jet merging and shock generation occur near the centers of the images. Red dashed arrows show the direction of individual plasma jet propagation toward the chamber center (as in Fig. 1). In the bottom left image, dots correspond to the intersection locations of the jet collision plane with the spectroscopy (purple) and interferometry (green) lines of sight (as in Fig. 1). The listed counts are over the whole image, for which 4096 counts is fully saturated.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00454v3-figure-e871d6cd-1041-4c47-bf27-e6eeef7529e9.jpg)

generates sharp gradients in density and temperature that will create sharp gradients in intensity. Although the intensity gradient does not map one-to-one with each parameter across the shock, we examine the intensity gradient to indicate the shock gradient location and length scale.

Qualitatively in the top 4 rows in Fig. 4, we see the jets merging to form two intensity peaks originating a finite distance away from the image center. The experimental intensity peaks illustrate an increase in Ar-II and decrease in He-I (corresponding to an increase in He-II) within the shock. Then, near the center of the collision region, the amount of He-I increases and Ar-II decreases. This double-peaked shock profile has been observed in prior work studying single-ion-species plasma jet collisions with a similar plasma gun setup 35,36 . Although this double-peaked feature is not as clearly seen in the shock created from the 100% Ar jet collisions shown in Fig. 3, the purpose of those images in Fig. 3 was to check the effectiveness of our spectral filters, so we did not spend time to experimentally tune the system in order to try to obtain better looking pictures. Compared to the 100% He jet collision in Fig. 3, the 97% He and 3% Ar multi-ion-species jet collision in Fig. 4 shows more emission using the Ar-II filter.

The bottom row in Fig. 4 displays the spatial profile of the shock structure produced by merging multi-ion-species plasma jets, which are directly imaged using the narrow band-

## 97% He 3% Ar Jets

FIG. 4. (Top 4 rows) Multi-ion-species plasma jet merging time evolution images of plasma self-emission using the Ar-II and He-I filters. Times defined such that t = 0 is the time at which the jets are first starting to merge. (Bottom row) Lineouts from the 2 µ s images (in top rows) for the Ar-II (blue line) and He-I (red line) filtered images; peak Ar-II intensity is at intensity of 1, and peak He-II intensity is at intensity of 0. Horizontal bars show the shock widths: the distances between the 50%-90% intensity values. Shaded regions show the standard deviation of distance between peak ion intensities and the zero position.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00454v3-figure-7b73064d-7636-42b2-8207-8a27fcc04977.jpg)

pass filters for Ar-II and He-I line emission. Normal to the shock front, we take lineouts with widths of 20 pixels ( ∼ 1 cm). In the laboratory frame, the shock fronts are moving away from the zero position. Across the emission spatial profile, the lineout intensities are scaled between 0 (lowest emission) and 1 (peak emission). For the Ar-II filtered lineout, 0 corresponds to background (lowest emission of singly ionized Ar), and 1 corresponds to the peak singly ionized Ar emission. For the He-I filtered lineout, 1 corresponds to background (peak emission of neutral He), and 0 corresponds to peak singly ionized He emission (lowest

emission of neutral He). Using these lineouts, we infer the shock-front separation between the peak intensities of Ar and He ions, showing that the He has moved farther ahead of the shock (closer to the pre-shock region) than the Ar. Additionally, we observe the shock width (horizontal bars in the bottom row of Fig. 4 near the -2 cm position), which we approximate by the transition distance between the middle and peak intensity, here taken as 50%-90% of the peak value. The 50% and 90% values are chosen to reduce the influence of varying background levels and spurious pixel intensities on the measurements. These experimental results are obtained using the shock profile at 2 µ s after the beginning of the plasma jet merging. This time is chosen by analyzing an image time sequence and correlating the density increase over time between the experimental time-resolved interferometry and the simulation results; see row 4 of Fig. 2.

Using the shock profile lineouts and plasma parameters, we compare the length scales of our collisional plasma shocks. Taking statistics from 34 experimental shock-front profiles at the 2 µ s time, we find a separation between the shock-front intensity peaks to be 0 68 . ± 0 17 . cm, a He shock width of 0 36 . ± 0 09 . cm (37 ± 9 post-shock He ion-ion mean free paths), and an Ar shock width of 0 52 . ± 0 11 cm (57 . ± 12 post-shock Ar ion-ion mean free paths). Additionally, we find the distance between the peak He ion intensity and the zero position to be 2 01 . ± 0 33 cm, and the distance between the peak Ar ion intensity and the . zero position to be 1 33 . ± 0 17 cm; these distances are shown as the shaded regions in the . bottom row of Fig. 4. Lastly, we find the ratio of the Ar shock width to He shock width to be 1 52 . ± 0 34. . These experimental uncertainty values are the shot-to-shot standard deviations. The experimental shock widths are in reasonable agreement to the theoretical prediction of &gt; 20 post-shock ion-ion mean free paths (for M &lt; 2) 3,17 . Furthermore, similar observations to our measurements of the ratios of shock widths for different species have been reported for neutral gas shocks, for which barodiffusion is present, and the lighter gas species also diffuses farther ahead of the neutral gas shock front compared to the heavier species 4 .

## V. COMPARISON OF SHOCK LENGTH SCALES TO MEAN FREE PATH AND INTERPENETRATION

Here, we compare the quantitative distances measured in the experiments to relevant length scales for the process of shock formation. This comparison allows us to better understand the physical picture.

Table II compares the characteristic length scale results for the multi-ion-species plasma jet merging. Post-shock length scales are calculated based on the following plasma conditions obtained from simulations (Sec. VI A) and listed in Table I: peak electron number density n e = 2 1 . × 10 15 cm -3 , peak T e = 2 4 eV, peak . T i,He = 4 6 eV, peak . T i,Ar = 5 5 eV, . and average ionization states ¯ of 0.89 for He and 1.47 for Ar. Z The post-shock mean free path for an ion (unprimed) colliding with multiple ion species (primed) is 47,48

$$L _ { m f p, i } & = \frac { v _ { t h, i } } { \sum _ { i ^ { \prime } } ( \frac { n _ { i ^ { \prime } } } { n _ { t o t, i } } ) \nu _ { \epsilon } ^ { i, i ^ { \prime } } } \\ &. \quad. \quad. \quad. \quad i \, i ^ { \prime } \, \ e \cdot \quad. \quad. \quad. \quad. \quad.$$

with thermal velocity v th,i and energy loss collision rate ν i,i ′ /epsilon1 . Since we are not in the fast or slow limit, we use the full Coulomb collision-rate formula 47 . In this formula, because we add contributions from individual ion species with discrete ionization states, Z , rather than using an average ¯ , Z the collision rate is weighted by the individual species number density fraction n i ′ /n tot,i . We compare the experimentally inferred shock width (as defined above in Sec. IV) to the ion-ion mean free path rather than to the electron-ion mean free path because the density jump and ion temperature jump within a shock front is predicted to occur over a distance of a few ion-ion mean free paths, while the electron temperature jump is predicted to occur over a (longer) distance of a few electron-ion mean free paths 3 . Furthermore, we calculate that the ion collision rates are dominant compared with the neutral collision rates ν n ≈ n πd v n 2 n th,n , where d n ∼ 0 3 nm is the approximate .

TABLE II. Plasma-jet-merging length scales (in centimeters).

| Length scale (cm)                   | He     | ∆    | Ar     |
|-------------------------------------|--------|------|--------|
| He-Ar shock-front separation (exp.) |        | 0.68 |        |
| He-Ar shock-front separation (sim.) |        | 0.50 |        |
| Shock width (exp.)                  | 0.36   |      | 0.52   |
| Shock width (sim.)                  | 0.44   |      | 0.57   |
| Post-shock mean free path (Eq. 4)   | 0.0097 |      | 0.0092 |
| Jet characteristic size             | 10     |      | 10     |
| Pre-shock interpenetration (Eq. 5)  | 0.058  |      | 0.740  |

Van der Waals atomic diameter, and also with the neutral charge exchange collision rates ν CE ≈ n σ n CE th,n v , where σ CE ≈ 3 × 10 -15 cm 2 is the charge exchange cross section 49,50 .

The ion-ion interpenetration length is the pre-shock slowing distance over which ions from an individual jet (test particle, unprimed) will stream through the other jet (field particles, primed) before becoming collisional 36,37,39,40,47,48,51 ,

$$L _ { s, i } & = \frac { v _ { r e l } } { 4 \sum _ { i ^ { \prime } } ( \frac { n _ { i ^ { \prime } } } { n _ { t o t, i } } ) \nu _ { s } ^ { i, i ^ { \prime } } }, \\ \\. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \,. \quad \$$

with relative velocity v rel and slowing-down collision rate ν i,i ′ s . Again, we use the entire collision-rate formula 47 . In our parameter space, the ion-ion slowing dominates over ionelectron slowing, and thus the latter is ignored. Because our jet size is much larger than the interpenetration distance, the merging of the two supersonic plasma jets produces a shock.

## VI. SIMULATION SHOCK PROFILE RESULTS

We now compare the shock profile, length scales, and emission intensities obtained in the experiment to simulation results. The shock profile lengths in the simulations agree with the experimental observations within uncertainty bars. There are disagreements between the emission intensities in the post-processed simulations compared to the experiments, particularly in the post-shock region.

## A. Comparison of Simulation and Experimental Shock Length Scales

We have employed multi-fluid simulations, using the chicago code 31-33 , in which each plasma jet is modeled as a separate ion fluid that can contain different species, to model the plasma collision 35,39,40 , including the interpenetration effects described in Sec. V. chicago is a hybrid particle-in-cell code with the capability to model both ions and electrons as fluid species 31 . The multi-fluid simulation can model interpenetration effects because each jet is treated as a separate fluid, and the jets can exchange momentum and energy via the inter-fluid collision models used in the code (see Eqs. 1 and 2 in Thoma et al., Physics of Plasmas 18, 103507 (2011) 31 ). The code may also be run in a magnetohydrodynamic (MHD) mode in which electron inertia is neglected. In the present research, the multi-ion MHD approach was found to generate identical results as those of lengthier simulations with electron inertia retained. For the 1D simulation presented in this work, the two jets are given an initial density profile, temperature, and velocity based on input from experimental measurements. The simulation grid cell length is 0.06 cm, and the shock profile results are similar to results from a simulation with a 0.015 cm grid cell length. propaceos non-LTE data are used for the ion equations of state and opacities 43 . Collisionality between ions and electrons is determined by a Spitzer model. At initial temperatures below 2 eV, there was no discernible interaction between the colliding jets because the low ionization states lead to large interpenetration lengths. The simulation may need to be better tuned in these

FIG. 5. (Top) Simulated density lineouts for the He and Ar in jets coming from the left (dashed lines) and right (solid lines). Total densities (dotted lines) are the sum of the left and right jets. Ar density values are divided by 0.03 for visualization. (Bottom) Corresponding simulated temperature lineouts. Times are at 2 µ s after jets first start to merge.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00454v3-figure-ae139737-d92d-4161-b4ce-ca65eb670edd.jpg)

low ionization regimes and include neutral collision rates in order to obtain more physically accurate results.

The chicago 1D multi-fluid simulations agree with the experimental results, showing: shock-front separation, He diffusing farther ahead of the shock (closer to the pre-shock region) than Ar, and different shock widths for different ions. Figure 5 displays the simulated densities and temperatures, respectively, for He and Ar in the individual plasma jets coming from the left and right (dashed and solid lines) at 2 µ s after the jets start to merge. We use these individual jet shock profiles to obtain the following quantities. Simulations indicate a separation between the shock-front density peaks to be 0 50 . ± 0 12 cm, a He shock width . of 0 44 . ± 0 12 cm (45 . ± 12 post-shock He ion-ion mean free paths), and an Ar shock width of 0 57 . ± 0 12 cm (62 . ± 13 post-shock Ar ion-ion mean free paths). Additionally, we find the distance between the peak He density and the zero position to be 1 51 . ± 0 06 . cm, and the distance between the peak Ar density and the zero position to be 1 01 . ± 0 06 cm. . We also find the ratio of the Ar shock width to He shock width to be 1 29 . ± 0 45. . The simulation uncertainty values are based upon the 0.06 cm simulation resolution. Except for the distances between the peak intensities and the zero position, these simulation quantities agree with the experimental data within uncertainty bars.

FIG. 6. (Top) Synthetic intensity lineouts of plasma self-emission using the Ar-II (500 ± 25 nm) and He-I (589 ± 5 nm) filters with spect3D 43 . Inputs for these intensities are the simulation total densities and temperatures shown in Fig. 5. Comparing with Fig. 4, the synthetic lineouts do not depict the experimentally observed double-peaked intensity feature. (Bottom) The ionization fractions of neutral He (He I), singly ionized He (He II), singly ionized Ar (Ar II), and doubly ionized Ar (Ar III) corresponding to the above synthetic intensity lineouts. The average ionization values for both He and Ar increase near the center post-shock region, despite the He-I filter intensity remaining large.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00454v3-figure-705ab832-07a7-40f0-b377-629a608282c6.jpg)

## B. Comparison of Simulation and Experimental Intensity Profiles

Through this point in the paper, the experiments and simulations have generally agreed. Because the experimental results in Fig. 4 are emission intensity and the simulation results in Fig. 5 are density and temperature, it is useful to directly compare the experiments with synthetic emission intensities from the simulation. However, in the simulation postshock region, there are significant disagreements to the experiment regarding the plasma state: particularly the ionization state and emission intensity. The simulation post-shock density prediction matches relatively well based on the experimental interferometry (see Fig. 2). The simulation post-shock ion temperature predictions are reasonable based on calculations converting 100% of the ion kinetic energy (at 11 km/s) into thermal energy within the shock, which would generate post-shock T He max ≈ 4 7 eV and . T Ar max ≈ 27 5 eV, . and we expect post-shock temperatures closer to the T He max value due to the relative atomic concentrations of the species. Post-shock ion temperatures were also measured in experiments for single ion species shocks, showing agreement with simulation results 39,40 . Although the simulation post-shock density and temperature values, which influence the intensity results, seem reasonable, the following paragraph will explore differences in the simulation ionization states and synthetic emission intensities compared to the experiments.

While the experiments and simulations quantitatively agree within uncertainty bars regarding the magnitude of the shock-front separation and the He and Ar shock widths, the

simulations predict a single-peaked shock structure near the collision center rather than the double-peaked and wider post-shock region observed in the experiments; see the dotted lines in the top row of Fig. 5 compared to the bottom row of Fig. 4. This difference correlates with the smaller distance between the center zero position and the ion density peaks found in the simulation compared to the experiments. As shown in the top row of Fig. 6, we use spect3D software 43 to obtain synthetic intensity lineouts of the chicago simulation density and temperature parameters. This synthetic intensity lineout does not reproduce the experimentally observed double-peaked shock profile or the decreases in relative intensity near the zero position for either Ar-II or He-I filter. Still, as shown in the bottom row of Fig. 6, the average ionization values for both species are increasing in the post-shock region near the center position of the plasma jet merging, as is expected. A 2D Cartesian simulation of the jet merging process provides qualitatively similar results as the 1D simulation, also depicting a single-peaked shock structure.

These differences in the post-shock region between experiments and simulation will be further assessed in future work. The disagreement seems to occur from the simulation predicting larger amounts of He-I (rather than He-II) in the shock compared to the experiment and also not showing a region of decreased Ar-II near the center of the collision region. Possible sources of this discrepancy include the models of collisionality and thermal conduction in the post-shock region used in the simulations for this plasma parameter regime. Additionally, the simulation EOS and opacity models do not capture time-dependent effects of ionization and/or recombination during the transient process of shock formation, as is present in the experiments. Ultimately, although the experimental and simulated emission intensities disagree in the post-shock region, the separation of ion species within the shock front agrees well between the simulations and experimental observations.

## VII. DOMINANCE OF BARODIFFUSION IN THE PRESENT WORK

The shock-front separation results are consistent among the experiments, simulations, and theoretical predictions of ion species diffusion. In this section, we quantitatively estimate the expected species concentration change within the shock and calculate the relative contributions of the different diffusion mechanisms, showing that barodiffusion dominates in the present research.

As stated earlier in Sec. II: within a shock front, theory predicts that the lighter ion diffusion velocity (in the center of mass frame) points toward the pre-shock region, and the heavier ion diffusion velocity points in the opposite direction. This prediction is consistent with our experimental and simulation results for the He/Ar mixture. The interspecies ion diffusion should change the relative species concentration. Experimentally, we did not directly measure the species concentration along the shock profile, but the simulations show interspecies ion diffusion. Compared to the initial condition of 97% He and 3% Ar atomic concentrations [76.4% He and 23.6% Ar mass concentrations], in the top row in Fig. 5 we find a minimum He atomic concentration in the post-shock region (corresponding to the zero position) of 95.8% [69.6% He mass concentration], and we find a maximum He atomic concentration at the shock front closer to the pre-shock region (corresponding to the ± 2 5 cm position) of 97.7% [80.7% He mass concentration]. .

We now quantitatively explore the effects of different diffusion flux terms on the species concentration. The diffusion terms depend upon the species concentration, masses, number densities, ionization states, and temperatures. The pre-shock and post-shock plasma parameters are listed in Table I. The goal of the following calculations are to serve as order-of-magnitude estimates to show consistency between the theory, experiments, and simulations. First, assuming no initial concentration separation for simplicity, we take ∇ c m 1 = 0 in Eq. 2 and Eq. 3. We also approximate ∂t = 2 µ s. Taking a characteristic length scale as the shock width L SW ≈ 0 4 cm, we estimate . ∇ = 1 /L SW . Therefore, for a quantity Q in Eq. 3, we have ( ∇ Q /Q ) ≈ (1 /L SW )( Q post -shock -Q pre -shock ) /Q pre -shock . For our He/Ar mixture, we then calculate the diffusion coefficients to be D = 600 cm 2 /s, κ p = 1 3, . κ Ti = 0 5, . κ Te ≈ -1, and κ e = 1 3 . 24,28,52,53 . See Appendix A for expressions for

the diffusion coefficients. We also estimate the ambipolar electric field ∇ φ ≈ ∇ T /e e 54,55 . Therefore, the relative contributions for each diffusion mechanism (the quantities in parentheses in Eq. 3) are 20.3 cm -1 for barodiffusion, 1.14 cm -1 for thermodiffusion (1.36 cm -1 for ion thermodiffusion and -0.23 cm -1 for electron thermodiffusion), and 0.30 cm -1 for electrodiffusion. Starting with a He concentration c m 1 = ρ /ρ 1 = 76 4%, in Eq. 1 we ob-. tain a He diffusion speed v D 1 = 1 7 . × 10 -2 cm/ µ s, corresponding in Eq. 2 to a change in concentration of ∂c m 1 = 6 5%. . This estimated calculation result is in agreement with the simulation concentration change from the above paragraph. Furthermore, these estimates show that barodiffusion is the dominant diffusion mechanism in the present research that causes the shock-front separation. This result contrasts with prior ICF experiments, for which the ion thermodiffusion mechanism dominated 19,20 . Additional experimental testing and validation of the theory can be achieved by changing the species concentrations, which influence the contributions of the various diffusion mechanisms.

## VIII. CONCLUSION

Within a collisional multi-ion-species plasma shock front containing 97% He and 3% Ar, we experimentally observe shock-front separation and species-dependent shock widths. Experiments, 1D multi-fluid plasma simulations, and theoretical predictions are all consistent in showing that the lighter He ions diffuse farther ahead within the overall shock front than the heavier Ar ions. The experimental shock profiles of different ion species were directly imaged using narrow bandpass visible wavelength filters. Multi-fluid plasma simulations allowed for reasonably accurate modeling of the plasma jet merging and multi-ion-species effects, though with disagreement on the emission intensities in the post-shock plasma. The fundamental experimental data in the present work can be used to validate models and benchmark numerical simulations of multi-ion-species collisional plasma shocks of relevance to HED, MIF, and ICF experiments. These experimental data are a valuable start to the process of benchmarking codes in this regime by demonstrating where there is agreement and disagreement. Additional work can be performed to obtain species concentration measurements at multiple points across the shock profile during a single experimental shot. These measurements will likely require a more detailed spectroscopic study of the emission line ratios in comparison with the prismspect atomic modeling calculations. More experiments can be conducted in order to acquire data for other times during the shock propagation, which can provide information about time-dependent effects during the process of shock formation including ionization, recombination, and kinetic (non-Maxwellian) effects.

## ACKNOWLEDGMENTS

We acknowledge J. Dunn, K. C. Yates, S. Brockington, A. Case, E. Cruz, F. D. Witherspoon, Y. C. F. Thio, D. Welch, P. Bellan, G. Kagan, and X.-Z. Tang for technical support and/or useful discussions. This work was supported by the Office of Fusion Energy Sciences and the Advanced Research Projects Agency-Energy of the U.S. Department of Energy under contract number DE-AC52-06NA25396.

## Appendix A: Diffusion Coefficient Calculations

For reference, we list the equations to determine the diffusion coefficients 24,28,52,53 , where subscript '1' refers to the lighter ion species with mass m 1 , and subscript '2' refers to the heavier ion species with mass m 2 . The classical diffusion coefficient, D, is

$$D = \frac { \rho T _ { i } } { A _ { 1 2 } \mu _ { 1 2 } n _ { 1 } \nu _ { 1 2 } } \times \frac { c _ { m 1 } ( 1 - c _ { m 1 } ) } { c _ { m 1 } m _ { 2 } + ( 1 - c _ { m 1 } ) m _ { 1 } } \quad \quad ( \mathbf A 1 )$$

where the new variables are the transport coefficient A 12 ( ∼ 1, here) and reduced mass µ 12 = m m / m 1 2 ( 1 + m 2 ). The collision frequency ν 12 is the energy loss collision rate. Next, the barodiffusion ratio, κ p , is

$$\kappa _ { p } = c _ { m 1 } ( 1 - c _ { m 1 } ) ( m _ { 2 } - m _ { 1 } ) ( \frac { c _ { m 1 } } { m _ { 1 } } + \frac { 1 - c _ { m 1 } } { m _ { 2 } } ) & & ( A 2 )$$

The electrodiffusion ratio, κ e , is

$$\kappa _ { e } = m _ { 1 } m _ { 2 } c _ { m 1 } ( 1 - c _ { m 1 } ) ( \frac { c _ { m 1 } } { m _ { 1 } } + \frac { 1 - c _ { m 1 } } { m _ { 2 } } ) ( \frac { Z _ { 1 } } { m _ { 1 } } - \frac { Z _ { 2 } } { m _ { 2 } } ) \quad \quad ( A 3 )$$

The ion thermodiffusion ratio, κ Ti , is

$$\kappa _ { T i } = m _ { 1 } m _ { 2 } ( \frac { c _ { m 1 } } { m _ { 1 } } + \frac { 1 - c _ { m 1 } } { m _ { 2 } } ) [ \frac { c _ { m 1 } B _ { 1 1 } } { m _ { 1 } } + \frac { ( 1 - c _ { m 1 } ) B _ { 1 2 } } { m _ { 2 } } ] \quad \quad ( A 4 )$$

and the electron thermodiffusion ratio, κ Te , is

$$\kappa _ { T e } = & m _ { 1 } m _ { 2 } ( \frac { c _ { m 1 } } { m _ { 1 } } + \frac { 1 - c _ { m 1 } } { m _ { 2 } } ) [ \frac { c _ { m 1 } Z _ { 1 } } { m _ { 1 } } + \frac { ( 1 - c _ { m 1 } ) Z _ { 2 } } { m _ { 2 } } ] \\ & \times [ ( 1 - c _ { m 1 } ) B _ { 1 e } - c _ { m 1 } B _ { 2 e } ] \frac { T _ { e } } { T _ { i } } \\ \varepsilon \quad \varepsilon \quad \cdot \varepsilon \quad. \quad \varepsilon \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \end{$$

where B 11 , B 12 , B 1 e , and B 2 e are transport coefficients that can be evaluated numerically. The electron thermodiffusion ratio, κ Te , can also be written as

$$\kappa _ { T e } & = - Z _ { 1 } ^ { 2 } c _ { m 1 } ( 1 - c _ { m 1 } ) ( c _ { m 1 } + \frac { ( 1 - c _ { m 1 } ) m _ { 1 } } { m _ { 2 } } ) \\ & \quad \times ( \frac { m _ { 2 } } { m _ { 1 } } - \frac { Z _ { 2 } ^ { 2 } } { Z _ { 1 } ^ { 2 } } ) \frac { T _ { e } } { T _ { i } } \frac { \beta _ { 0 } } { Z _ { e f f } } \\ \varpi \cdot \cdot \cdot \cdot \end{matrix}$$

where the coefficient β 0 is

$$\mathbb { m } { \Phi } _ { \Phi } = \frac { 3 0 Z _ { e f f } ( 1 1 Z _ { e f f } + 1 5 \sqrt { 2 } ) } { 2 1 7 Z _ { e f f } ^ { 2 } + 6 0 4 \sqrt { 2 } Z _ { e f f } + 2 8 8 }$$

where Z eff is

$$Z _ { e f f } = \frac { n _ { 1 } Z _ { 1 } ^ { 2 } + n _ { 2 } Z _ { 2 } ^ { 2 } } { n _ { 1 } Z _ { 1 } + n _ { 2 } Z _ { 2 } }$$

1 T. Cowling, 'VI. The influence of diffusion on the propagation of shock waves,' The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science 33 , 61-67 (1942).

2

J.

D.

Jukes,

'The structure

of

a

shock wave

in

a

fully ionized

gas,'

Journal of Fluid Mechanics

3

,

275-285 (1957).

3 M. Y. Jaffrin and R. F. Probstein, 'Structure of a Plasma Shock Wave,' Physics of Fluids 7 , 1658 (1964). 4 R. E. Center, 'Measurement of Shock-Wave Structure in Helium-Argon Mixtures,' Physics of Fluids 10 , 1777 (1967).

5 J. R. Rygg, J. A. Frenje, C. K. Li, F. H. S´ eguin, R. D. Petrasso, J. A. Delettrez, V. Y. Glebov, V. N. Goncharov, D. D. Meyerhofer, S. P. Regan, T. C. Sangster, and C. Stoeckl, 'Tests of the hydrodynamic equivalence of direct-drive implosions with different D2 and He3 mixtures,' Physics of Plasmas 13 , 052702 (2006).

6 H. W. Herrmann, J. R. Langenbrunner, J. M. Mack, J. H. Cooley, D. C. Wilson, S. C. Evans, T. J. Sedillo, G. A. Kyrala, S. E. Caldwell, C. S. Young, A. Nobile, J. Wermer, S. Paglieri, A. M. McEvoy, Y. Kim, S. H. Batha, C. J. Horsfield, D. Drew, W. Garbett, M. Rubery, V. Y. Glebov, S. Roberts, and J. A. Frenje, 'Anomalous yield reduction in direct-drive deuterium/tritium implosions due to H3e addition,' Physics of Plasmas 16 , 056312 (2009).

7 D. T. Casey, J. A. Frenje, M. Gatu Johnson, M. J.-E. Manuel, H. G. Rinderknecht, N. Sinenian, F. H. S´ eguin, C. K. Li, R. D. Petrasso, P. B. Radha, J. A. Delettrez, V. Y. Glebov, D. D. Meyerhofer, T. C. Sangster, D. P. McNabb, P. A. Amendt, R. N. Boyd, J. R. Rygg, H. W. Herrmann, Y. H. Kim, and A. D. Bacher, 'Evidence for Stratification of Deuterium-Tritium Fuel in Inertial Confinement Fusion Implosions,' Physical Review Letters 108 , 075002 (2012).

- 8 H. G. Rinderknecht, M. J. Rosenberg, C. K. Li, N. M. Hoffman, G. Kagan, A. B. Zylstra, H. Sio, J. A. Frenje, M. Gatu Johnson, F. H. S´ eguin, R. D. Petrasso, P. Amendt, C. Bellei, S. Wilks, J. Delettrez, V. Y. Glebov, C. Stoeckl, T. C. Sangster, D. D. Meyerhofer, and A. Nikroo, 'Ion Thermal Decoupling and Species Separation in Shock-Driven Implosions,' Physical Review Letters 114 , 025001 (2015).
- 9 D. P. Higginson, J. S. Ross, D. D. Ryutov, F. Fiuza, S. C. Wilks, E. P. Hartouni, R. Hatarik, C. M. Huntington, J. Kilkenny, B. Lahmann, C. K. Li, A. Link, R. D. Petrasso, B. B. Pollock, B. A. Remington, H. G. Rinderknecht, Y. Sakawa, H. Sio, G. F. Swadling, S. Weber, A. B. Zylstra, and H.-S. Park, 'Kinetic effects on neutron generation in moderately collisional interpenetrating plasma flows,' Physics of Plasmas 26 , 012113 (2019).
- 10 H. Sio, J. A. Frenje, A. Le, S. Atzeni, T. J. T. Kwan, M. Gatu Johnson, G. Kagan, C. Stoeckl, C. K. Li, C. E. Parker, C. J. Forrest, V. Glebov, N. V. Kabadi, A. Bose, H. G. Rinderknecht, P. Amendt, D. T. Casey, R. Mancini, W. T. Taitano, B. Keenan, A. N. Simakov, L. Chac´n, S. P. Regan, T. C. o Sangster, E. M. Campbell, F. H. Seguin, and R. D. Petrasso, 'Observations of Multiple Nuclear Reaction Histories and Fuel-Ion Species Dynamics in Shock-Driven Inertial Confinement Fusion Implosions,' Physical Review Letters 122 , 035001 (2019).
- 11 O. Larroche, 'Ion Fokker-Planck simulation of D3 He gas target implosions,' Physics of Plasmas 19 , 122706 (2012).
- 12 C. Bellei, P. A. Amendt, S. C. Wilks, M. G. Haines, D. T. Casey, C. K. Li, R. Petrasso, and D. R. Welch, 'Species separation in inertial confinement fusion fuels,' Physics of Plasmas 20 , 012701 (2013).
- 13 C. Bellei, H. Rinderknecht, A. Zylstra, M. Rosenberg, H. Sio, C. K. Li, R. Petrasso, S. C. Wilks, and P. A. Amendt, 'Species separation and kinetic effects in collisional plasma shocks,' Physics of Plasmas 21 , 056310 (2014).
- 14 A. Inglebert, B. Canaud, and O. Larroche, 'Species separation and modification of neutron diagnostics in inertial-confinement fusion,' EPL (Europhysics Letters) 107 , 65003 (2014).
- 15 E. Vold, R. Rauenzahn, and A. N. Simakov, 'Multi-species plasma transport in 1D direct-drive ICF simulations,' Physics of Plasmas 26 , 032706 (2019).
- 16 P. Amendt, O. L. Landen, H. F. Robey, C. K. Li, and R. D. Petrasso, 'Plasma Barodiffusion in InertialConfinement-Fusion Implosions: Application to Observed Yield Anomalies in Thermonuclear Fuel Mixtures,' Physical Review Letters 105 , 115005 (2010).
- 17 P. Amendt, S. C. Wilks, C. Bellei, C. K. Li, and R. D. Petrasso, 'The potential role of electric fields and plasma barodiffusion on the inertial confinement fusion database,' Physics of Plasmas 18 , 056308 (2011). 18 G. Kagan and X.-Z. Tang, 'Thermo-diffusion in inertially confined plasmas,' Physics Letters A 378 , 1531-1535 (2014).
- 19 S. C. Hsu, T. R. Joshi, P. Hakel, E. L. Vold, M. J. Schmitt, N. M. Hoffman, R. M. Rauenzahn, G. Kagan, X.-Z. Tang, R. C. Mancini, Y. Kim, and H. W. Herrmann, 'Observation of interspecies ion separation in inertial-confinement-fusion implosions,' EPL (Europhysics Letters) 115 , 65001 (2016).
- 20
- T. R. Joshi, P. Hakel, S. C. Hsu, E. L. Vold, M. J. Schmitt, N. M. Hoffman, R. M. Rauenzahn, G. Kagan, X.-Z. Tang, R. C. Mancini, Y. Kim, and H. W. Herrmann, 'Observation and modeling of interspecies ion separation in inertial confinement fusion implosions via imaging x-ray spectroscopy,' Physics of Plasmas 24 , 056305 (2017).
- 21 H. G. Rinderknecht, H.-S. Park, J. S. Ross, P. A. Amendt, S. C. Wilks, J. Katz, N. M. Hoffman, G. Kagan, E. L. Vold, B. D. Keenan, A. N. Simakov, and L. Chac´n, 'Measurements of ion velocity separation and o ionization in multi-species plasma shocks,' Physics of Plasmas 25 , 056312 (2018).
- 22 H. G. Rinderknecht, H.-S. Park, J. S. Ross, P. A. Amendt, D. P. Higginson, S. C. Wilks, D. Haberberger, J. Katz, D. H. Froula, N. M. Hoffman, G. Kagan, B. D. Keenan, and E. L. Vold, 'Highly Resolved Measurements of a Developing Strong Collisional Plasma Shock,' Physical Review Letters 120 , 095001 (2018).
- 23 T. R. Joshi, S. C. Hsu, P. Hakel, N. M. Hoffman, H. Sio, and R. C. Mancini, 'Progress on observations of interspecies ion separation in inertial-confinement-fusion implosions via imaging x-ray spectroscopy,' Phys. Plasmas 26 , 062702 (2019).
- 24 G. Kagan and X.-Z. Tang, 'Electro-diffusion in a plasma with two ion species,' Physics of Plasmas 19 , 082709 (2012).
- 25 C. Bellei and P. A. Amendt, 'Asymptotic separation in multispecies collisional plasma shocks,' Physical Review E 90 , 013101 (2014).
- 26 A. N. Simakov and K. Molvig, 'Hydrodynamic description of an unmagnetized plasma with multiple ion species. I. General formulation,' Physics of Plasmas 23 , 032115 (2016).

27

S. I. Glazyrin, A. S. Kuratov, and V. Y. Bychenkov, 'Separation of ions on the front of a shock wave in

a multicomponent plasma,' JETP Letters

103

,

238-243 (2016).

- 28 A. N. Simakov, B. D. Keenan, W. T. Taitano, and L. Chac´n, 'Plasma ion stratification by weak planar o shocks,' Physics of Plasmas 24 , 092702 (2017).
- 29 B. D. Keenan, A. N. Simakov, L. Chac´n, o and W. T. Taitano, 'Deciphering the kinetic structure of multi-ion plasma shocks,' Physical Review E 96 , 053203 (2017).
- 30 B. D. Keenan, A. N. Simakov, W. T. Taitano, and L. Chac´n, 'Ion species stratification within strong o shocks in two-ion plasmas,' Physics of Plasmas 25 , 032103 (2018).
- 31 C. Thoma, D. R. Welch, R. E. Clark, N. Bruner, J. J. MacFarlane, and I. E. Golovkin, 'Two-fluid electromagnetic simulations of plasma-jet acceleration with detailed equation-of-state,' Physics of Plasmas 18 , 103507 (2011).

- 32 C. Thoma, D. R. Welch, and S. C. Hsu, 'Particle-in-cell simulations of collisionless shock formation via head-on merging of two laboratory supersonic plasma jets,' Physics of Plasmas 20 , 082128 (2013).
- 33 C. Thoma, D. R. Welch, R. E. Clark, D. V. Rose, and I. E. Golovkin, 'Hybrid-PIC modeling of laser-plasma interactions and hot electron generation in gold hohlraum walls,' Physics of Plasmas 24 , 062707 (2017).
- 34 S. C. Hsu, E. C. Merritt, A. L. Moser, T. J. Awe, S. J. E. Brockington, J. S. Davis, C. S. Adams, A. Case, J. T. Cassibry, J. P. Dunn, M. A. Gilmore, A. G. Lynn, S. J. Messer, and F. D. Witherspoon, 'Experimental characterization of railgun-driven supersonic plasma jets motivated by high energy density physics applications,' Physics of Plasmas 19 , 123514 (2012).
- 35 E. C. Merritt, A. L. Moser, S. C. Hsu, J. Loverich, and M. Gilmore, 'Experimental Characterization of the Stagnation Layer between Two Obliquely Merging Supersonic Plasma Jets,' Physical Review Letters 111 , 085003 (2013).
- 36 E. C. Merritt, A. L. Moser, S. C. Hsu, C. S. Adams, J. P. Dunn, A. Miguel Holgado, and M. A. Gilmore, 'Experimental evidence for collisional shock formation via two obliquely merging supersonic plasma jets,' Phys. Plasmas 21 , 055703 (2014).
- 37 A. L. Moser and S. C. Hsu, 'Experimental characterization of a transition from collisionless to collisional interaction between head-on-merging supersonic plasma jetsa),' Physics of Plasmas 22 , 055707 (2015). 38
- S. C. Hsu, S. J. Langendorf, K. C. Yates, J. P. Dunn, S. Brockington, A. Case, E. Cruz, F. D. Witherspoon, M. A. Gilmore, J. T. Cassibry, R. Samulyak, P. Stoltz, K. Schillo, W. Shih, K. Beckwith, and Y. C. F. Thio, 'Experiment to Form and Characterize a Section of a Spherically Imploding Plasma Liner,' IEEE Transactions on Plasma Science 46 , 1951-1961 (2018).
- 39 S. J. Langendorf, K. C. Yates, S. C. Hsu, C. Thoma, and M. Gilmore, 'Experimental Measurements of Ion Heating in Collisional Plasma Shocks and Interpenetrating Supersonic Plasma Flows,' Physical Review Letters 121 , 185001 (2018).
- 40 S. J. Langendorf, K. C. Yates, S. C. Hsu, C. Thoma, and M. Gilmore, 'Experimental study of ion heating in obliquely merging hypersonic plasma jets,' Physics of Plasmas 26 , 082110 (2019).
- 41 E. C. Merritt, A. G. Lynn, M. A. Gilmore, and S. C. Hsu, 'Multi-chord fiber-coupled interferometer with a long coherence length laser,' Review of Scientific Instruments 83 , 033506 (2012).
- 42 E. C. Merritt, A. G. Lynn, M. A. Gilmore, C. Thoma, J. Loverich, and S. C. Hsu, 'Multi-chord fiber-coupled interferometry of supersonic plasma jets (invited),' Review of Scientific Instruments 83 , 10D523 (2012).
- 43 Prism Computational Sciences, 'http://www.prism-cs.com,' .

44

E.

V.

Stenson and

P.

M.

Bellan,

'Dual-Species

Plasmas

Illustrate

MHD

Flows,'

IEEE Transactions on Plasma Science

36

,

1206-1207 (2008).

- 45 E. V. Stenson and P. M. Bellan, 'Magnetically Driven Flows in Arched Plasma Structures,' Physical Review Letters 109 , 075001 (2012).
- 46
- National Institute of Standards and Technology Atomic Spectra Database, 'https://www.nist.gov/pml/atomic-spectra-database,' (2019).
- 47 J. Huba, NRL Plasma Formulary (Naval Research Laboratory, Washington, DC, 2016) p. 31.
- 48 B. A. Trubnikov, 'Particle Interactions in a Fully Ionized Plasma,' Reviews of Plasma Physics 1 , 105 (1965).
- 49 H. Helm, 'The cross section for symmetric charge exchange of He + in He at energies between 0.3 and 8 eV,' Journal of Physics B: Atomic and Molecular Physics 10 , 3683-3697 (1977).
- 50 B. M. Smirnov, 'Tables for Cross Sections of the Resonant Charge Exchange Process,' Physica Scripta 61 , 595-602 (2000).
- 51 S. Messer, A. Case, L. Wu, S. Brockington, and F. D. Witherspoon, 'Nonlinear compressions in merging plasma jets,' Physics of Plasmas 20 , 032306 (2013).
- 52 G. Kagan and X.-Z. Tang, 'Mathematica Code to Calculate Diffusion Ratios,' (2019).
- 53 V. M. Zhdanov, 'Transport Processes in Multicomponent Plasma,' Plasma Physics and Controlled Fusion 44 , 2283-2283 (2002).
- 54 P. A. Amendt, J. L. Milovich, S. C. Wilks, C. K. Li, R. D. Petrasso, and F. H. S´ eguin, 'Electric field and ionization-gradient effects on inertial-confinement-fusion implosions,' Plasma Physics and Controlled Fusion 51 , 124048 (2009).
- 55 X.-Z. Tang, H. L. Berk, Z. Guo, and C. J. McDevitt, 'Reduced Fokker-Planck models for fast particle distribution across a transition layer of disparate plasma temperatures,' Physics of Plasmas 21 , 032707 (2014).
Tom Byvank Samuel J. Langendorf Carsten Thoma Scott C. Hsu
Observation of Shock-Front Separation in Multi-Ion-Species Collisional Plasma Shocks
We observe shock-front separation and species-dependent shock widths in
multi-ion-species collisional plasma shocks, which are produced by obliquely
merging plasma jets of a He/Ar mixture (97% He and 3% Ar by initial number
density) on the Plasma Liner Experiment [S. C. Hsu et al., IEEE Trans. Plasma
Sci. 46, 1951 (2018)]. Visible plasma emission near the He-I 587.6 nm and Ar-II
476.5-514.5 nm lines are simultaneously recorded by splitting a single visible
image of the shock into two different fast-framing cameras with different
narrow bandpass filters (589 +/- 5 nm for observing the He-I line and 500 +/-
25 nm for the Ar-II lines). For conditions in these experiments (pre-shock ion
and electron densities ~5*10^14 cm^-3, ion and electron temperatures of ~2.2
eV, and relative plasma-merging speed of 22 km/s), the observationally inferred
magnitude of He/Ar shock-front separation and the shock widths themselves are <
1 cm, which correspond to ~50 post-shock thermal ion-ion mean free paths. These
experimental lengths scales are in reasonable qualitative and quantitative
agreement with results from 1D multi-fluid simulations using the Chicago code.
However, there are differences between the experimentally-inferred and
simulation-predicted ionization states and line emission intensities,
particularly in the post-shock region. Overall, the experimental and simulation
results are consistent with theoretical predictions that the lighter He ions
diffuse farther ahead within the overall shock front than the heavier Ar ions.