Recent theoretical and experimental results have suggested that launching a spherically convergent shock wave at the end of the acceleration phase improves the ignition conditions for inertial confinement fusion (ICF) implosions.[
Chinese Physics B, Volume. 29, Issue 10, (2020)
Hot-electron deposition and implosion mechanisms within electron shock ignition
A hot-electron driven scheme can be more effective than a laser-driven scheme within suitable hot-electron energy and target density. In our one-dimensional (1D) radiation hydrodynamic simulations, 20× pressure enhancement was achieved when the ignitor laser spike was replaced with a 60-keV hot-electron spike in a shock ignition target designed for the National Ignition Facility (NIF), which can lead to greater shell velocity. Higher hot-spot pressure at the deceleration phase was obtained owing to the greater shell velocity. More cold shell material is ablated into the hot spot, and it benefits the increases of the hot-spot pressure. Higher gain and a wider ignition window can be observed in the hot-electron-driven shock ignition.
1. Introduction
Recent theoretical and experimental results have suggested that launching a spherically convergent shock wave at the end of the acceleration phase improves the ignition conditions for inertial confinement fusion (ICF) implosions.[
In a shock-ignition implosion, the main pulse used to assemble the dense core is a conventional low-adiabat laser pulse and the target is typically a thick cryogenic shell.[
An even larger enhancement of the hot-spot pressure can be achieved if the late shock is driven by energetic particles, which is called “electron shock ignition”.[
2. Hot electron energy deposition
Kilo-electron-Volt (keV) energetic electrons can be generated by the nonlinear precess of laser plasma interactions, including the three-wave parametric instabilities (LPI): stimulated Raman scattering (SRS), two plasmon decay (TPD), and filamentation instability, when the laser intensity is around 5 × 1014 W/cm2 to 2 × 1015 W/cm2 for conventional direct- and indirect-drive (ICF).[
The energy deposition of energetic electrons in plasmas has been widely studied.[
The Maxwellian distribution of hot electrons is widely used.[
Figure 1.The normalized Maxwellian distribution of hot-electron temperature of 40 keV, 60 keV, and 100 keV.
The accurate calculations of hot-electron transport and energy deposition are PIC and Monte Carlo simulations. Most Monte Carlo codes use tables for the electron stopping power and transport scattering cross sections obtained for normalized media[
In the Solodov–Betti model, the cross sections for fast electrons scattering off plasma ions and electrons are obtained in the first Born approximation.[
The target used in our work was taken from the 1D hydrocode LILAC target design for NIF shock ignition.[
Figure 2.The energy depositions (solid lines) for the monoenergetic hot electrons with energies of 40 keV, 60 keV, and 100 keV. The density profile is shown as the dashed line.
The integro-differential diffusion equation of the multiple scattering problem in an infinite, homogeneous plasma was studied.[
Figure 1.The normalized Maxwellian distribution of hot-electron temperature of 40 keV, 60 keV, and 100 keV.
A hot-electron transport and energy deposition package is developed according the above model and put into hydrocode DEC2D. Within DT plasma (A = 2.5), the relationship of hot-electron energy deposition to the propagated areal density is presented in Fig. 4. The monoenergic hot electron with an energy of 100 keV is reported to have a stopping range of ∼ 15 mg/cm2 in DT plasma.[
Figure 2.The energy depositions (solid lines) for the monoenergetic hot electrons with energies of 40 keV, 60 keV, and 100 keV. The density profile is shown as the dashed line.
3. Gain characteristics of electron-driven shock ignition
The shock-ignition target on the NIF was designed[
Figure 4.The energy deposition (blue line) and the locally deposited flux (red line) of the 100-keV monoenergetic hot electron in DT plasma.
The hot-electron-driven shock ignition gains versus ignitor launching times are also demonstrated in Fig. 5. A 200-ps hot-electron spike was used in the present study with Maxwellian energy distribution. The hot electron temperature was set at 60 keV according to the recent OMEGA experimental data.[
4. Optimum implosion of electron-driven shock ignition
The ignitor shock strength is one of the most important parameters of shock ignition. Different ignitor shock launching schemes have different ignitor shock properties and can lead to different implosion characteristics. In this work, the laser spike-driven ignitor shocks and the hot-electron spike-driven ignitor shocks are presented in Fig. 6. The laser-driven ignitor shock is chosen as the highest gain ∼ 40 of the LILAC target design, which is shown in Fig. 5 at the laser spike launching time of 9.6 ns. The hot-electron-driven ignitor shock is chosen as the highest gain ∼ 106 of DEC2D simulations, which is also shown in Fig. 5 at the hot-elelctron spike launching time of 10.3 ns. The laser energy is deposited in the corona zone until the critical density nc = 1.1 × 1021 / λ2 cm−3,[
Figure 5.The gains versus the ignitor shock launching times. The black line represents LILAC calculations with the laser-driven shock ignition, and the colored lines are 1D simulations with the hot-electron-driven shock ignition with different hot-electron energies.
For the hot-electron-driven case, the hot electrons are generated in the corona zone and depend on the non-linear lPI processes. For instance, if TPD is dominated, most of the hot electrons are produced in nc/4, which is the same as that on OMEGA direct-drive implosions. If SRS is dominated, the hot electrons are produced in some specific areas of corona, which depends on the plasma parameters and noise levels, and the SRS dominated nonlinear LPI process occurs in the high laser intensity (shock-ignition relevant) experiments, as performed recently on OMEGA.[
The ignitor shock compresses the shell and the center low-density gas region, which is called the “hot spot”. The pressure of the hot spot keeps increasing because of the geometry convergence, the shock transmission, and the mass increase from the ablation of the inner surface of the low-temperature, high-density shell. When the hot-spot pressure equals the shell pressure, the acceleration phase ends, the shell velocity starts to decrease, and the deceleration phase comes out. For the laser-driven shock ignition with a laser launching time of 9.6 ns, the trajectory of the shell peak density and the shell velocity are presented in Fig. 7(a). There is a peak density location jump at around 9.8 ns toward the outside surface, which is caused by the laser-driven shock catching up with the outside shell surface [see Fig. 6(c)], this jump compresses the outside part of the shell into a high density. The laser pulse turns off at 9.8 ns, but the shell velocity keeps increasing from 9.8 ns to 10.6 ns, owing to the compression by the hot-spot pressure to the inside part of the shell. The mass center moves toward the target center. At 10.6 ns, there is a sudden enhancement in the shell velocity from 284 km/s to 294 km/s because the incoming and outcoming shocks collide. During this process a very high density zone (∼ 3× that of the density before collision) can be generated near the inner surface of the shell.
Figure 6.The density and pressure profiles after the ignitor shock launched for 0, 100, and 200 ps. Panels (a)–(c) for the laser spike driven, and panels (d)–(f) for the hot-electron spike driven. The highest gain targets in Fig.
The mass center trajectory and shell velocity of the electron-driven shock ignition with the ignitor shock launching time of 10.3 ns are shown in Fig. 7(b). The stagnation occurs at about 10.8 ns, which is the same for the laser-driven shock ignition in Fig. 7(a). The time duration from the ending of the hot-electron spike to the stagnation is ∼ 300 ps, which is ∼ 1/3 of that in the laser-driven shock ignition. The reason is that the hot-spot pressure is comparable to the shell pressure in the hot-electron-driven shock ignition when the hot-electron spike turns off, and the hot-spot pressure is negligible in comparison with the shell pressure in the laser-driven shock ignition when the laser spike turns off. As a result, it takes less time for the hot-electron-driven shock ignition to achieve an equivalent result between the hot spot and the shell, and it implies less perturbation growth in the deceleration phase, which benefits the implosion performance. At 10.6 ns, a similar sudden enhancement of the shell velocity is observed within the electron shock ignition because of the collision of the incoming and outcoming shocks. However, the shell velocity increases from 385 km/s to 484 km/s, which has more shell velocity enhancement than that of the laser-driven shock ignition (from 284 km/s to 294 km/s). It implies that more intense shock collision occurs and stronger incoming shock can be generated after the collision in the hot-electron-driven shock ignition, which benefits the second piston compression effect of the hot spot. The shell velocity at the end of the laser spike is 250 km/s and the maximum shell velocity (just after the collision) is 294 km/s, and the corresponding shell velocities of the hot-electron-driven shock ignition are 338 km/s and 484 km/s.
In an implosion, the thermonuclear instability is triggered in the central hot spot. Nuclear reactions occur in the hot spot with a cross section, and self-ignition occurs when ρR ≥ 0.3 g/cm2 and Ti ≥ 5 keV in the hot spot. The ignition products includes high-energy α particles and neutrons. The α particles form a burning wave toward the outside low-temperature and high-density shell, and deposits energy in the shell depending on the shell’s areal density and temperature. In simulations, the burn-wave effect is not always considered and the α particles are turned off to investigate the pure hydrodynamic of the implosion. In Fig. 8, the stagnation parameters of the laser-driven shock ignition with a laser spike launching time of 9.6 ns, and the hot-electron-driven shock ignition with a hot-electron spike launching time of 10.2 ns are demonstrated with α particles turned off. The density profiles in Fig. 8(a) show that the shell density of the laser-driven shock ignition at stagnation is greater than that of the hot-electron-driven shock ignition. However, the hot-spot density of the laser-driven shock ignition (∼ 93 g/cm3 in the target center) is less than that of the hot-electron-driven shock ignition (∼ 144 g/cm3 in the target center). The shell thickness of the laser-driven shock ignition is less than that of the hot-electron-driven shock igniton with both of the two cases having almost the same outside shell surface. The hot-spot pressure of the laser-driven shock ignition (378 Gbar in the target center) is much less than that of the hot-electron-driven shock ignition (644 Gbar in the target center). It has been shown above that the hot-electron-driven shock ignition implosion has greater shell velocity compared to that of the laser-driven shock ignition, which implies that greater hot-spot pressure can be achieved within the hot-electron-driven shock ignition, according to the scaling law[
Figure 8.The target density profile (a), pressure profile (b), temperature profile (c), and adiabat profile (d) at stagnation without burn wave for the laser- and hot-electron-driven shock ignitions.
The evolutions of the neutron rate in the laser- and hot-electron-driven shock ignitions with burn wave are shown in Fig. 9(a). The confinement time of the ignition implosion can be written as τ ∼ Rs/cs, where Rs is the hot-spot radius and cs is the sound speed of the burn wave at the peak neutron rate. For the laser-driven shock ignition, τ ∼ 40 ps; for the hot-electron-driven shock ignition, τ ∼ 23 ps. However, the peak neutron rate is greater for the hot-electron-driven shock ignition than that for the laser-driven shock ignition. The peak neutron rate along with the density, pressure, and ion temperature profiles are presented in Figs. 9(b)–9(d). At peak neutron rate, the hot-spot density for the laser-driven shock ignition is 74 g/cm3, and for the hot-electron-driven shock ignition is 165 g/cm3. The hot-spot pressure for the laser-driven shock ignition is 1672 Gbar, and for the hot-electron-driven shock ignition is 5172 Gbar. The ion temperature profiles for the laser- and hot-electron-driven shock ignitions are almost the same in Fig. 9(d). In the electron shock ignition, more cold shell material is ablated into the high temperature hot-spot, which increases the hot-spot density and pressure, and degrades the hot spot temperature.
Figure 9.(a) The neutron rate, (b) target density profile, (c) pressure profile, (d) and ion temperature profile at peak neutron rate with burn wave for the laser- and hot-electron-driven shock ignitions.
5. Conclusions
The hot-electron energy deposition in plasma depends on the hot-electron energy and the material density. In a planar target, for a laser to hot-electron energy conversion efficiency of 25%, greater ablation pressure can be obtained by the hot-electron-driven than that by the laser driven with plasma density greater than 1.6 g/cm3. In our spherical 1D simulations, 20× pressure enhancement is obtained by the hot-electron-driven shock ignition compared to that by the laser-driven shock ignition. The hot-electron-driven scheme can be much more effective than the laser-driven scheme with suitable hot-electron energy and plasma density.
In this work, the Solodov–Betti hot-electron stopping power model was utilized, binary collisions including electron–electron and electron–ion collision and exciting plasma waves effect were considered, the Lewis’ multiple scattering theory was used to calculate the spatial moments of the hot-electron energy deposition. According to the recent OMEGA experiment data, the hot-electron temperature of 60 keV was used in our simulations. More than 2 Gbar of shock was generated in our work within the suitable hot-electron spike launching time. The hot-electron shock ignition can obtain greater shell velocity than that of the laser-driven shock ignition, which implies more hot-spot pressure in the deceleration phase, and more gain in the implosion. The hot-spot temperature profiles of the laser- and hot-electron-driven shock ignition are almost the same, the reason being that greater shell velocity in the hot-electron shock ignition leads to higher hot-spot temperature, and it also ablates more inner surface of the cold shell into the hot spot and decreases the hot-spot temperature. In the simulations more than 10-μg shell mass was ablated into the hot spot in the hot-electron shock-ignition than that in the laser-driven shock ignition, and the ablated mass into the hot spot increased the hot-spot density and pressure and also decreases the hot-spot temperature. On the other side, the hot electron can propagate deeper than the laser, and it leads to a higher adiabat for the shell in the hot-electron driven shock ignition than that of the laser-driven shock ignition, which broadens the shell and degrades the shell areal density. However, the increased shell velocity can surmount the degradation of the shell areal density caused by the increase of the adiabat in the hot-electron-driven shock ignition. The hot-electron-driven shock ignition can be a robust (big ignition window), high gain (> 100) ignition scheme in ICF.
The self-generated electromagnetic field is not taken into considered in these 1D simulations. Three-dimensional (3D) extended-magnetohydrodynamic simulations have proved that the electromagnetic fields lead to significant distortions to the plasmas temperature and density, which degrade the implosion performance.[
[1] H B Cai, K Mima, W M Zhou, T Jozaki, H Nagatomo, A Sunahara, R J Mason. Phys. Rev. Lett., 102(2009).
[2] W Theobald, R Betti, C Stoeckl, K S Anderson, J A Delettrez, V Yu Glebov, V N Goncharov, F J Marshall, D N Maywar, R L McCrory, D D Meyerhofer, P B Radha, T C Sangster, W Seka, D Shvarts, V A Smalyuk, A A Solodov, B Yaakobi, C D Zhou, J A Frenje, C K Li, F H Séguin, R D Petrasso, L J Perkins. Phys. Plasmas, 15(2008).
[3] R Betti, C D Zhou, K S Anderson, J L Perkins, W Theobald, A A Solodov. Phys. Rev. Lett., 98(2007).
[4] S Atzeni, A Marocchino, A Schiavi, G Schurtz. New J. Phys., 15(2013).
[5] W Shang, H Wei, Z Li, R Yi, T Zhu, T Song, C Huang, J Yang. Phys. Plasmas, 20(2013).
[6] M Lafon, X Ribeyre, G Schurtz. Phys. Plasmas, 20(2013).
[7] L J Perkins, R Betti, K N LaFortune, W H Williams. Phys. Rev. Lett., 103(2009).
[8] B Canaud, M Temporal. New J. Phys., 12(2010).
[9] W L Shang, R Betti, S X Hu, K Woo, L Hao, C Ren, A R Christopherson, A Bose, W Theobald. Phys. Rev. Lett., 119(2017).
[10] S Gus’kov, X Ribeyre, M Touati, J L Feugeas, Ph Nicola, V Tikhonchuk. Phys. Rev. Lett., 109(2012).
[11] X Ribeyre, S Gus’kov, J L Feugeas, Ph Nicola, V T Tikhonchuk. Phys. Plasmas, 20(2013).
[12] A R Piriz, S A Piriz, N A Tahir. Phys. Plasmas, 20(2013).
[13] E M Campbell, W J Hogan. Plasma Phys. Control. Fusion, 41, B39(1999).
[14] C D Zhou, R Betti. Phys. Plasmas, 14(2007).
[15] R Nora, R Betti, K S Anderson, A Shvydky, A Bose, K M Woo, A R Christopherson, J A Marozas, T J B Collins, P B Radha, S X Hu, R Epstein, F J Marshall, R L McCrory, T C Sangster, D D Meyerhofer. Phys. Plasmas, 21(2014).
[16] R K Kirkwood, J D Moody, J Kline, E Dewald, S Glenzer, L Divol, P Michel, D Hinkel, R Berger, E Williams, J Milovich, L Yin, H Rose, B MacGowan, O Landen, M Rosen, J Lindl. Plasma Phys. Control. Fusion, 55(2013).
[17] W Seka, D H Edgell, J F Myatt, A V Maximov, R W Short, V N Goncharov, H A Baldis. Phys. Plasmas, 16(2009).
[18] A J Kemp, F Fiuza, A Debayle, T Johzaki, W B Mori, P K Patel, Y Sentoku, L O Silva. Nucl. Fusion, 54(2014).
[19] W Theobald, R Nora, W Seka, M Lafon, K S Anderson, M Hohenberger, F J Marshall, D T Michel, A A Solodov, C Stoeckl, D H Edgell, B Yaakobi, A Casner, C Reverdin, X Ribeyre, A Shvydky, A Vallet, J Peebles, F N Beg, M S Wei, R Betti. Phys. Plasmas, 22(2015).
[20] A A Solodov, R Betti. Phys. Plasmas, 15(2008).
[21] C K Li, R D Petrasso. Phys. Rev. E, 73(2006).
[22] H G Reik, H Risken. Phys. Rev., 124, 777(1961).
[23] J Delettrez, R Epstein, M C Richardson, P A Jaanimagi, B L Henke. Phys. Rev. A, 36, 3926(1987).
[24] K Anderson, R Betti, T A Gardiner. Bull. Am. Phys. Soc., 46, 280(2001).
[25] A Bose, K M Woo, R Nora, R Betti. Phys. Plasmas, 22(2015).
[26] K M Woo, R Betti, A Bose, R Epstein, J A Delettrez, K S Anderson, R Yan, P Y Chang, D Jonathan, M Charissis(2015).
[27] K S Anderson, R Betti, P W McKenty, J B Collins, M Hohenberger, W Theobald, R S Craxton, J A Delettrez, M Lafon, J A Marozas, R Nora, S Skupsky, A Shvydky. Phys. Plasmas, 20(2013).
[28] V N Goncharov, T C Sangster, T R Boehly, S X Hu, I V Igumenshchev, F J Marshall, R L McCrory, D D Meyerhofer, P B Radha, W Seka. Phys. Rev. Lett., 14(2010).
[29] H W Lewis. Phys. Rev., 78, 526(1950).
[30] X Ribeyre, G Schurtz, M Lafon, S Galera, S Weber. Plasma Phys. Control. Fusion, 51(2009).
[31] W Shang, J Yang, Y Dong. Appl. Phys. Lett., 102(2013).
[32] J Lindl. Phys. Plasmas, 2, 3933(1995).
[33] C A Walsh, J P Chittenden, K McGlinchey, N P L Niasse, B D Appelbe. Phys. Rev. Lett., 118(2017).
Get Citation
Copy Citation Text
Wan-Li Shang, Xing-Sen Che, Ao Sun, Hua-Bing Du, Guo-Hong Yang, Min-Xi Wei, Li-Fei Hou, Yi-Meng Yang, Wen-Hai Zhang, Shao-Yong Tu, Feng Wang, Hai-En He, Jia-Min Yang, Shao-En Jiang, Bao-Han Zhang. Hot-electron deposition and implosion mechanisms within electron shock ignition[J]. Chinese Physics B, 2020, 29(10):
Received: Jun. 1, 2020
Accepted: --
Published Online: Apr. 21, 2021
The Author Email: Wan-Li Shang (wanlishang@gmail.com)