Effect of atomic number and pressure on plasma pinch properties and characteristic soft x-ray emission in PF1000

In this study, the Lee code is used to compute the characteristics soft x-ray yield (Ysxr) production for nitrogen (N2), oxygen (O2), neon (Ne), and argon (Ar) and bremsstrahlung radiation for hydrogen (H2), deuterium (D2), and helium (He) with pressure variation in PF1000 of 2.5–2.6 MA for D2. In the calculation of characteristic soft x-ray, the corresponding temperature windows of the said gases are set into the code at which they are ionized to their H-like and He-like levels. The focus pinch parameters such as radius ratio (minimum radius of plasma pinch column/anode radius), ion density, specific heat ratio, pinch energy density, self-absorption correction factor, and maximum induced voltage are computed at the optimum pressure of each gas. The obtained pinch plasma temperature range (1.2–2.2) × 106 K of H2, D2, and He is sufficiently high for fully ionized plasmas and the resulting bremsstrahlung radiation (14 J) for He is significantly larger than for H2 (0.26 J) and D2 (0.62 J). The optimum Ysxr of Ne(∼9314 J) at 0.51 Torr with pinch energy density (PED) (∼26 × 108 Jm−3) is found to be the highest whilst for Ar(∼7 J) at 0.019 Torr with (∼1.2 × 108 Jm−3) is the lowest. It is found that the radius ratio (∼0.05) of Ne is 3-fold smaller than that (∼0.16) in Ar. This enhancement of compression in pinch of Ne increases the ion density significantly by a factor of 253 than in Ar gas. Thus, the results show a strong correlation of plasma pinch properties with Ysxr for various gases.


Introduction
The study of pinch plasmas as a bright soft x-ray source is well established. Although it is simple to construct [1,2] the dense plasma focus (DPF) emits x-ray emission efficiently when compared with other pinches including the X-pinch [3] and vacuum spark [4]. Applications in the x-ray spectrum range include spectroscopy [5], microscopy and microlithography [6], for pumping lasers [7], crystallography [8], medical and industrial radiography [9], for back-lighter [10], and for engineering micromachining [11]. A super-hot (∼100s eV) and super-dense (∼10 20 cm −3 ) short-lived plasmas is concentrated in the focus pinch column by electromagnetic acceleration and compression [6]. During the intense pinch compression phase, high levels of ionization are attained with efficient transfer of energy into predominantly radiative recombination and line emission [12] in the soft x-ray wavelength region. It has been shown that the radius of the pinch may be determined by the operation of work-energy balance and magnetic-kinetic pressure balance [13].
The ratio of the minimum radius of plasma pinch column (r min ) to anode radius (a) is, to some extent dependent on temperature through the specific heat ratio (γ). The lower the value of γ, the greater the degree of freedom, hence the greater the compressibility in the pinch plasma column [14]. The lower atomic number gases (e.g. H 2 , D 2 , He) are fully ionized in the focus pinch and the dominant radiations are free-free bremsstrahlung. The high atomic number gases such as Ne to Xe are freely ionizing at the temperatures of the plasma focus pinch and emit predominantly line radiation [15]. The line radiation powers are orders of magnitude greater than Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI. bremsstrahlung power. The analysis [16,17] shows that the line radiation powers are sufficient to cool the high atomic number pinch to the extent of the radiative collapse. On the other hand, the bremsstrahlung in (low atomic number) pinches is hardly sufficient to affect the pinch dynamics.
To optimise the x-ray yield, changes have been made to capacitor bank parameters including bank inductance and discharge current, electrode material, dimension and form [18,19], insulator length and material [18], the gas its operational pressure [10]. Such optimizations are important for applications, especially for specialized material fabrication and modification of surface properties [20]. The Lee code combines the actual electric circuit with the effects of DPF dynamics. Its equations incorporate the thermodynamics and in the pinch phase couple to radiation. The code covers a wide range of gases [21]. The self-absorption of the radiation by the volume of the plasma is included to improve the realism of the code when the plasma becomes optically with sufficient density [15,22].
The code was widely used for the design and understanding of a number of machines including sub-kJ DPF devices [23], kJ machines FNII [24], the UBA [25] hard X-radiation source, the KSU plasma focus [26] and the cascading PF [27]. Information in the output includes trajectories and speeds for the dynamics of the axial and radial phases [28], focus pinch duration and dimensions, pinch properties including densities and temperatures, radiation yields, and characteristics including soft x-ray [28]. The code is used for the systematic optimization of plasma focus devices [28].
In NX2 plasma focus, this code verified the measured values of soft x-ray yield and variation of Ne pressure. The effect of operating pressure on radiation yield for various gases is studied in low-energy DPF device [28,29]. As a megajoule DPF device, x-ray yield investigations in PF1000 are not very frequent since its prime object is to achieve fast neutron yield [30,31]. A numerical experiment was performed for Ne gas with pressure variation in PF1000 [32] to find optimum yield conditions. In this present work, we extend the numerical investigation to study the effect of atomic number and P 0 on radiations and characteristic soft x-ray in the standard PF1000. Radiation-enhanced compression in the focused plasma are studied in lighter (H 2 , D 2 , and He) and heavier (N 2 , O 2 , Ne, and Ar) gases through simulation work. The computed properties including pinch current, pinch duration, ion density, specific heat ration, radius ratio (radius of plasma pinch column/anode radius), pinch column length, energy (from capacitor bank) injected into the pinch plasma (EINP), and (PED) are correlated with radiations and characteristic soft x-ray at the optimum P 0 for each gas in PF1000. Dividing EINP by the volume of the pinch column, PED is computed. The consequent effects of operating pressure and atomic number on radiations and characteristic soft x-ray emission are then discussed. This paper is organized as follows: in the next section, we consider the procedure of numerical experiments to calculate x-rays and radiations configuring the Lee code with the standard parameters of PF1000 and the computed results of characteristics soft x-ray including bremsstrahlung radiations for different gases with pressure variation are shown in section 3. In section 4, we discuss our findings, and finally, section 5 concludes the present work. -´), respectively [32]. The corresponding temperature window for each gas is set so that the PF is operating with the gas ionized to its H-like and He-like levels. Such radiation is referred to be termed the 'characteristic' soft x-ray of the particular gas rather than just the soft x-ray. The gas can also radiate soft x-ray outside the window, but the soft x-ray will not be its characteristic soft x-ray; meaning not from H-like and He-like levels. So, in the RADPF5.15 code, the line radiation yield (Q L ) is computed as the characteristic soft x-ray yield at the corresponding temperature window for each of these gases [32].
The line radiation is computed as follows: Where, we take Y sxr = Q L, by operating in the correct temperature-window according to the gas. The radiation is integrated over the time of the pinch, the soft x-ray energy generated in the pinch depends on: effective ion charge (Z eff ), ion density (N i ) and temperature (T), pinch length (Z max ), and pinch radius, and atomic number (Z). This soft x-ray energy is emitted after reduction by self-absorption of the plasma which has dependence on temperature and density.

Configuring the Lee code for PF1000
The code is configured with the capacitor bank characteristics, tube dimensions and operating voltage and pressure of the standard PF1000 [7]: Bank parameters: Static inductance (L 0 ) = 33.5 nH, Capacitance (C 0 ) = 1332 μF, Stray circuit resistance (r 0 ) = 6.1 mΩ.
Model (mass swept-up and effective current) parameters: for axial phase, ( f m ) = 0.13, ( f c ) = 0.7 and ( f mr ) = 0.35 and ( f cr ) = 0.65 in the radial phase. These model parameters have been obtained by fitting the computed and measured current waveforms at 3.5 Torr D 2 in PF1000 [7]. The experimental current trace for different gases in this device is not available. Though this set of model parameters is expected to vary considerably within a wide range of gases, even in a pressure range for a given gas, these are kept fixed throughout our present series of numerical experiments for different gases [4,15,17].

PF1000 with hydrogen, deuterium, and helium
The effects of atomic number and pressure on pinch plasma properties are tabulated (see table 1). We found the temperature of pinch plasmas as 1.15 × 10 6 K at 2.48 Torr H 2 , 1.36 × 10 6 K at 2.49 Torr D 2 , and 2.20 × 10 6 K at 1.48 Torr He. Since, H 2 , D 2 , and He are fully ionized in the PF pinch and hence the radiations for them are completely bremsstrahlung (Y br ).
The bremsstrahlung radiation from He plasma (14 J) at optimum pressure of 1.48 Torr is considerably higher than from H 2 (0.26 J) and D 2 (0.62 J). This is because of the increased compression and the lower specific heat ratio in He. In He pinch plasma, (see table 1

PF1000 with nitrogen, oxygen, neon, and argon
The computed maximum characteristic soft x-ray yields and pinch properties at optimum pressure for N 2 , O 2 , Ne, and Ar gases are recorded (see table 2). The operating pressure was varied in the range of 0.01 to 2.5 Torr. Through a series of numerical experiments, the obtained optimum pressure for each gas at which the maximum characteristic soft x-ray is presented here.
The computed total discharge currents rise to the peak values of 1978 kA at 6.60 μs, 1906.4 kA at 6.48 μs, 1773.5 kA at 6.23 μs, and 1145.1 kA at 3.29 μs for N 2 , O 2 , Ne, and Ar, respectively (see figure 1). The current peaks come earlier and are lower because the dynamic loads (0.5 × dL/dt = 10 -7 × ln (b/a) × axial speed) added in the circuit due to their axial motion of the current sheath. The peak axial current sheath speeds of N 2 , O 2 , Ne, and Ar are found 8.8 cm μs −1 , 10.0 cm μs −1 , 12.5 cm μs −1 , and 28.2 cm μs −1 , respectively (see table 2). Also, the pinch currents decreases and the dips arise earlier (see figure 1) due to the larger radial dynamic loads with higher magnetic piston speeds for high atomic number gases (see table 2). In our optimization, we considered Q L = Y sxr within the corresponding temperature window for N 2 , O 2 , Ne, and Ar. The maximum characteristic soft x-ray for high atomic number gas is found at lower pressure compared to that of lighter gas (see figure 3). In this study, the maximum characteristic soft x-ray yield (∼9314 J) is obtained Table 1. At optimum pressure, computed maximum radiation yields for H 2 , D 2 , and He with fixed V 0 = 27 kV, f m = 0.13, f c = 0.7, f mr = 0.35, and f cr = 0.65 in PF1000. Whilst for Ar, it requires much higher pinch temperature (∼16.5 × 10 6 K) to reach its optimum characteristic soft x-ray emission regime since need to reach the 16th and 17th ionized level. That's why we need to go to such low pressure (∼0.019 Torr) so as to get to high shock speeds (∼50.5 cm μs −1 ). But such high temperature gives low overall line radiation (∼6.69 J) in Ar, particularly when we need to go to such low  pressures. So, operating PF in Ar to radiate its characteristic soft x-ray is not conducive to much radiation and we do not get radiative collapse, not even radiative cooling.
When we find Ar data over a range of pressures, there is a small range of pressures (∼0.1-0.35 Torr) there is immense line radiation (∼21780 J at 0.14 Torr, r min /a ∼ 0.009, T pinch ∼5 × 10 6 K) (NOT characteristic radiation) and that is the range when radiative effects on dynamics including radiation collapse is observed in Ar. It was found for N 2 , O 2 , and Ne (see table 2), a smaller radius ratio than typical and that is radiation enhancing compression (also to lesser expect γ effects). Especially Ne, we are close to its most compressed state due to radiative enhancement of compression. Figure 4 shows the optimum characteristic soft x-ray yield increases with the high atomic number up to Ne(Z = 10) and then reduces drastically for Ar(Z = 18). The optimum characteristic soft x-ray yield for Ne(∼9314 J) was the highest whilst it was minimum in Ar(∼6.69 J).
Through a series of numerical experiments, we find an inequality relation in radius ratio (k = r min /a): k min (Ar) > k min (N 2 ) > k min (O 2 ) > k min (Ne). It is clear that a reduction in the specific heat ratio leads to a corresponding reduction in the radius ratio. This reduction in k min leads to an increase in the characteristic soft x-ray emission. Within the corresponding temperature window, the γ(Ne) ∼ 1.49, the r min (Ne) ∼ 0.61 cm gives k min ∼ 0.05 at optimum P 0 = 0.51 Torr whereas the γ(Ar) ∼ 1.57, the r min (Ar) of 1.82 cm, k min ∼ 0.16 and at optimum P 0 = 0.019 Torr. The k min (Ne) is a more than 3-fold reduction from k min (Ar) gives an increase in plasma density N i (Ne) = 7.59 × 10 23 m −3 to N i (Ar) = 0.03 × 10 23 m −3 ; an increase in pinch density by 253 times. This enhanced pinch density is conducive to greater production of characteristic soft x-ray yield in Ne from Ar.  In the notation of the code, AB (self-absorption correction factor) = 1 when self-absorption is zero. When AB reaches 1/e, meaning substantial self-absorption of its emission, the radiation mode is considered to transition from volume to surface radiation [17]. We find within the temperature regime, Z eff = 16 consequently AB(Ar) ∼1.0 and in the case of Ne, Z eff = 8.2 consequently AB(Ne) ∼0.98. This higher value of AB(Ar) reduces its characteristic soft x-ray emission from other gases Ne.

Discussions
The energy and intensity of radiations in a DPF device can be controlled by experimental conditions that include operating energy, electric current, electrode geometry, type of gas, and pressure. The lighter gases, H 2 and D 2 are fully ionized at 0.464 × 10 6 K (40 eV) and He is fully ionized at 1.856 × 10 6 K (160 eV) [2]. In this present paper, the pinch plasma temperatures are found as 1.15 × 10 6 K, 1.36 × 10 6 K, and 2.20 × 10 6 K for H 2 , D 2 , and He, respectively. Here, all the obtained temperatures are within the range of fully ionization temperature and hence the radiations for these gases are completely bremsstrahlung. This study presents that the optimum bremsstrahlung radiation increases with high atomic number gases.
In the energy range of 0.7-15 keV, the time-resolved soft x-ray signals (not absolute yield) were measured by Si-PIN detector from PF1000 device operating with pure Ne in the injection of D 2 and without puffing D 2 [33]. It was observed that the soft x-ray emission decreases considerably from this device in operation with Ne without injection of D 2 . Also, the soft x-ray emission from this PF1000 device was studied operating with D 2 and puffing Ne and the intense soft x-ray emission was measured in presence of Ne at the central region of the pinch [34]. However, no quantitative measurement of soft x-ray yield from PF1000 is available for Ne and other gases.
The radiation that come from H-like and He-like ions of an atom is equivalent to K-shell radiation (Y k ) or characteristics soft x-ray. In Z-pinch, K-shell radiation yield follows the scaling law as Y I k peak 4 for low current and Y I k peak 2 for high current [35,36], whilst in PF-pinch, it follows the scaling law as Y I sxr peak

3.2
for low 0.1 to 2.4 MA current [37]. Also, experimentally it is observed that the K-shell radiation from the Z-pinch increases with decreasing pinch radius along with increasing peak current [37]. In our numerical studies using Lee code, we also found a similar trend in dependency of characteristic soft x-ray yield on peak discharge current as well as pinch radius. For example, in Z-pinch with argon gas it was experimentally measured that Y k = 440 J with I peak = 1600 kA and r min = 0.17 cm [38], whilst in our numerical studies we obtained that Y k ∼ 6.69 J with I peak ∼1145 kA and r min ∼ 1.82 cm for argon gas. Here, the peak discharge current of our studies is lower (∼ 455 kA), whilst the pinch radius is very high (∼1.65 cm) compare to the K-shell radiation in Z-pinch [38]. Due to these lower peak current and larger pinch radius, our optimum Y k with argon is too smaller than that from Z-pinch. Again, for neon gas K-shell radiation from Z-pinch plasma has been measured as [39]: Y k = 4000 J with I peak = 1450 kA and r min < 0.2 cm whilst in our numerical studies we obtained as Y k = 9314 J with I peak = 1773.5 kA and r min = 0.61 cm. Here, the peak discharge current of our studies is higher (∼323.5 kA) though the pinch radius is slightly higher (∼ 0.41 cm) compare to those of K-shell radiation in Z-pinch [39]. In addition, experimentally it was observed that the K-shell radiation yield for argon was 800-1000 J whilst for neon it was 4000-5000 J from the same Z-pinch plasma system [40]. That means the H-like and He-like radiations for neon are higher than that of argon which is similar in nature to our numerical studies. For the other two gases (nitrogen and oxygen) the variation trend of soft x-ray yield with peak current and pinch radius like K-shell radiations in Z-pinch. Therefore, we may conclude that our computed results are justified with those from Z-pinch.
By comparing our results of numerical simulation with the research results [4,32], both the characteristics soft x-ray and PED increase with the increase in atomic number (N 2 , O 2 , and Ne) of used gases due to two factors: radius ratio and the effective ion charge while in our study we continue the increase in atomic number (Ar) of gas. We also focus on the effect of induced voltage, absorption correction factor, and ion density on the characteristics soft x-ray as it is found that both Y sxr and PED increase for high atomic number gases (N 2 , O 2 , and Ne) and then sharply reduce for Ar gas because the radius ratio and hence the number of ions inside the pinch of Ar plasma reduces significantly.

Conclusion
In this paper, the effect of atomic number and pressure on pinch properties, characteristic soft x-ray (for N 2 , O 2 , Ne, and Ar) and bremsstrahlung (for H 2 , D 2 , and He) radiations are studied using the modified Lee code in PF1000. The results show the optimum characteristic soft x-ray emission for high atomic number gas is obtained at lower pressure with high pinch temperature. A large enhancement of soft x-ray yield within the temperature window is found when the radius ratio could be reduced to a small value. In these present numerical studies, the emission of the optimum characteristic soft x-ray increases from 3879.93 J (N 2 ) to 9314 J (Ne) then drastically falls to 6.69 J (Ar).
The pinch temperatures of H 2 , D 2 , and He plasmas are found to be sufficiently high for fully ionized plasmas and emit completely bremsstrahlung radiations. The radiation from He plasma (14 J) at optimum pressure of 1.48 Torr is considerably higher than from H 2 (0.26 J) and D 2 (0.62 J).