Introduction
Upon interaction of lasers with intensity I L > 1018 W/cm2 and pulse duration τL < 20 ps with solid targets the laser radiation pressure may lead to the generation of relativistic shock waves, i.e., the thermodynamic parameters of the equilibrium states are related to the relativistic Rankine–Hugoniot equations. These shock waves may be driven by different mechanisms. One of the models describing these shocks is the laser-piston model. For laser intensities I L > 1018 W/cm2 the laser ponderomotive force pushes electrons ahead, so that the charge separation field forms a double layer. This double layer structure, called a laser piston, drives a shock compression wave moving in the unperturbed over-dense plasma. The double layer (DL) separates the propagation path of the laser pulse from the shocked plasma. The structure of the piston and the relation between its velocity and the laser intensity were described analytically and as well obtained in particle in cell simulations (Esirkepov et al., Reference Esirkepov, Borghesi, Bulanov, Mourou and Tajima2004; Naumova et al., Reference Naumova, Schlegel, Tikhonchuk, Labaune, Sokolov and Mourou2009; Schlegel et al., Reference Schlegel, Naumova, Tikhonchuk, Labaune, Sokolov and Mourou2009; Eliezer et al., Reference Eliezer, Nissim, Raicher and Martinez Val2014, Reference Eliezer, Martinez-Val, Henis, Nissim, Pinhasi, Ravid, Werdiger and Raicher2016; Schmidt and Boine-Frankenheim, Reference Schmidt and Boine-Frankenheim2016). The laser piston as a mechanism of particle acceleration to relativistic velocities was described in papers by Robinson et al., Reference Robinson, Gibbon, Zepf, Kar, Evans and Bellei2009 and Macchi, Reference Macchi2013 and references therein. Two-fluid simulations of laser plasma interaction where the nonlinear ponderomotive force was predominant and was predicted by ultrahigh acceleration of plasma blocks (Hora, Reference Hora2012).
Shock waves may be divided in two types, depending on the dissipation of kinetic energy at the shock front. In collisional shocks the kinetic energy dissipation is provided by Coulomb binary collisions and as a result the shock front is a few mean free path thick. In collisionless shock waves instabilities due to collective behavior of the plasma or reflection of ions from the shock front provide energy dissipation at the front, and the shock front can be several orders of magnitude smaller than the mean free path for binary collisions. In a high-intensity laser plasma interaction ions are usually heated indirectly by the electrons, ions gain thermal energy through electron ion-thermal equilibration or electron-ion instabilities. For short-pulse laser-induced shock waves an interesting question is whether electrons and ions establish their own post-shock temperature. Bulk heating due to laser piston-induced shock waves is considered here. Hot layer formation by sub-picosecond laser irradiation at intensities higher than 1019 W/cm2 has been reported (Theobald et al., Reference Theobald, Akli, Clarke, Delettrez, Freeman, Glenzer, Green, Gregori, Heathcote, Izumi, King, Koch, Kuba, Lancaster, MacKinnon, Key, Mileham, Myatt, Neely, Norreys, Park, Pasley, Patel, Regan, Sawada, Shepherd, Snavely, Stephens, Stoeckl, Storm, Zhang and Sangster2006) and described by a model including refluxing and confinement of supra-thermal in target volume. Shock heating of a thin layer at near solid density to 5 keV by 5·1020 W/cm2 laser radiation pressure driven shock was measured by Akli et al., Reference Akli, Hansen, Kemp, Freeman, Beg, Clark, Chen, Hey, Hatchett, Highbarger, Giraldez, Green, Gregori, Lancaster, Ma, MacKinnon, Norrey, Patel, Shearer, Stephens, Stoeckl, Storm, Theobald, Van Woerkom, Weber and Key2008.
In a recent paper (Eliezer et al., Reference Eliezer, Pinhasi, Martinez-Val, Raicher and Henis2017) we considered the shock wave induced by the laser piston propagating into aluminum at normal density, in particular the heating of the material produced during the shock compression, based on energy dissipation by binary collisions. It was shown that the shocked plasma might have different electron and ion temperatures, depending on the laser intensity and duration. In that paper (Eliezer et al., Reference Eliezer, Pinhasi, Martinez-Val, Raicher and Henis2017) it was assumed that energy was dissipated at a constant rate. The novel aspects of this paper are based on a new energy dissipation rate where the energy deposition rate is temperature dependent. Furthermore, we extend the model to additional materials and analyze the conditions that allow the formation of a collisional shock wave.
Section “Laser piston-induced relativistic shock wave” presents the relativistic Rankine–Hugoniot equations describing the shock wave in the material and the dependence of the shock parameters on the laser intensity. Section “Thermalization in laser-induced shock waves” proposes a model of the plasma heating produced during the laser piston-induced shock wave. Section “Numerical results” presents the numerical results for deuterium, carbon and aluminum targets and Section “Summary” concludes the paper.
Laser piston-induced relativistic shock wave
The shock wave induced by the laser piston is described by the relativistic Rankine–Hugoniot equations, relating the shock pressure P, energy density e, mass density ρ, particle (piston) velocity u p and shock velocity u s. The quantities with subscript zero are the corresponding material parameters before the shock arrival:
where c is the speed of light. The velocities u p and u s are given in the laboratory frame of reference.
The relativistic shock wave of Eq. (1) yields the following non-relativistic well known Hugoniot equations, when the velocities v satisfy v/c ≪ 1, and e = ρc 2 + ρE, where P and ρE are much smaller than ρc 2:
We assume an ideal equation of state (EOS)
Γ = (c p/c v) is the specific heat ratio related to the number of degrees of freedom per particle f by Γ = 1 + (2/f). We assume Γ = (5/3), f = 3. In the non-relativistic limit the above EOS is
The laser piston velocity u p, or βp = (u p/c), as a function of the laser intensity I L can be obtained by solving the relativistic Rankine–Hugoniot equation (1) together with the piston model equation:
where P is the radiation pressure equal to the shock pressure. The relativistic Rankine–Hugoniot equations are solved below for dimensionless variables: compression κ = (ρ/ρ0), dimensionless pressure Π = (P/ρ0c 2) and dimensionless laser intensity ΠL = (I L/ρ0c 3) (Eliezer et al., Reference Eliezer, Nissim, Raicher and Martinez Val2014).
Substituting the ideal EOS [Eq. (3)] with the Rankine–Hugoniot equation (1), we obtain the relativistic Hugoniot equation for dimensionless variables:
In Eqs. (7, 8) κ 0 = (Γ + 1/Γ − 1) is the non-relativistic asymptotic compression in the limit of infinite shock pressure and Π0 = (P 0/ρ0c 2). The solution of the Hugoniot equations for dimensionless variables is material and initial density independent.
We consider here the dimensionless pressure regime:
This regime corresponds to piston velocities in the range of 0.001–0.1c. The compression as a function of the dimensionless pressure, obtained from solving Eqs. (1, 3) is shown in Figure 1a and a zoom displaying the range considered here is shown in Figure 1b. The non-relativistic shock wave regime corresponds to compressions less than four, followed by an intermediate regime where the compression increases very slowly, while the pressure increases by about 10 orders of magnitude and a relativistic regime where there is no upper limit to the compression, unlike for the non-relativistic Hugoniot Eq. (2). The laser intensities generating the pressures considered here are from several times 1017 W/cm2 to 1022 W/cm2, depending on the material initial density. Figure 2 shows the piston/particle and the shock velocity as a function of laser intensity, obtained from the solution of the relativistic Hugoniot equation (1) coupled to the laser piston relation (4), for material initial density ρ0 = 1 g/cm3. Increasing the initial density shifts the laser intensities linearly to higher values for the same piston velocity.
Thermalization in laser-induced shock waves
Due to the shock wave formation in the material, the ions and the electrons behind the shock front are moving with the piston velocity u p. For non-relativistic, and/or nanosecond or longer time duration shock waves the piston work is divided equally to kinetic energy $\int \lpar {\rho u_{\rm p}^2 /2} \rpar dV$ and an increase in internal/thermal energy $\displaystyle \int {3 \over 2}k_{\rm B}\lpar {n_{\rm e}T_{\rm e} + n_{\rm i}T_{\rm i}} \rpar dV $ (Zeldovich and Raizer Reference Zeldovich and Raizer1966). For relativistic and semi-relativistic regimes or shorter laser pulse duration the shocked plasma may have in general different ion and electron shock temperatures, T i and T e (Eliezer et al., Reference Eliezer, Henis, Nissim, Pinhasi and Martinez Val2015). Moreover, for extremely short-laser pulse duration and high-laser intensity, the ions might not have enough time to achieve thermalization among themselves during the shock duration. The electrons reach thermalization at time shorter than the time scales considered here. Therefore, the laser piston work may be partitioned differently between kinetic and internal energy. In this paper we address this issue, estimating the time of temperature equilibration among the ions and the time of temperature equilibration between the electrons and the ions. It is shown that for given laser intensity there is a threshold pulse duration for the formation of a shock wave. In addition, dependent on the laser intensity, the electron and ion temperatures at the end of the pulse duration may be different.
The time-dependent equations for the electron and ion temperatures T e and T i are obtained from the energy conservation of the electrons and the ions:
The electron and ion densities are n e and n i, and k B is the Boltzmann constant. W ie (erg/cm3 s) is the rate of the electron–ion exchange energy density and W B (erg/cm3 s) is the appropriate bremsstrahlung losses of the electron energy density. The energy deposition into the ions $W_{\rm d}^{\rm i} $(erg/cm3) and the electrons $W_{\rm d}^{\rm e} $(erg/cm3) is obtained from the piston work. It is important to point out that the upper limit of $W_{\rm d}^{\rm i} + W_{\rm d}^{\rm e} $ equals 50% of the piston energy density. $\nu _\epsilon ^{\alpha /\beta} $ is the rate of energy loss of test particles α in a background of field particles β. α, β = e, i. These rates are density, temperature, and velocity dependent.
The bremsstrahlung power loss term is:
Here the electron temperature is in eV units and Z av is the average ionization in the material.
The temperature equilibration term between the electrons and the ions is:
The electron ion thermal equilibration time τeq is given by:
where m e and m i are the electron and ion masses and e is the electron charge.
From the fact that 50% of the piston work can contribute to thermal energy and other losses we obtain
The last part of Eq. (15) results from the non-relativistic Hugoniot conservation equations.
The deposited energy W d is calculated after solving the Hugoniot relativistic equations for the pressure and the particle and shock velocities as a function of the laser intensity.
We consider the energy deposition to the electrons and the ions, as test particles and as field particles, separately. The energy deposition is distributed according to the electron and ion masses.
Due to the large ion-electron mass ratio $W_{\rm d}^{\rm i} \gg W_{\rm d}^{\rm e} $.
Assuming that the laser pulse rise time is very small in comparison with the laser pulse duration, at the shock wave front the target particle velocity changes instantaneously from zero to a velocity u p. This is equivalent to a shock wave rise time much shorter than the laser pulse duration. We consider here collisional shock waves, i.e., we assume that the time equilibrations to reach the temperatures T e(t < τL) and T i(t < τL) are obtained by Coulomb binary collisions near the shock wave front. In this model we describe the piston energy deposition of the shocked particles, electrons, and ions into thermal energy using a relaxation rate arising from the interaction of test particles, labeled α, streaming with velocity u p, through a background of field particles, labeled β with a collision frequency of energy deposition (Huba Reference Huba2011) given by:
Here α and β stand for the electrons or ions, α/β denotes kinetic energy transferred from test α to field β particles, m α and m β stands for the electron and ion mass, $\nu _0^{\alpha /\beta} $, the relaxation rate scale and the function ψ are defined by
T β denotes the temperature of the field particles, v is the test particle velocity, v = u p, q α and q β are the charges (q equals the electron charge e for the electrons and Z ave for the ions), and k B is the Boltzmann constant. λαβ = ln Λαβ is the Coulomb logarithm. The transfer rate $\nu _\varepsilon ^{\alpha /\beta} $ is positive for $\varepsilon \gt \varepsilon _\alpha ^* $, and negative for $\varepsilon \lt \varepsilon _\alpha ^* $, where $\varepsilon = (1/2)m_\alpha v_\alpha ^2 $ and $x^* = (m_\beta /m_\alpha )(\varepsilon _\alpha ^* /T_\beta )$ are the solution of ψ ′(x*) = (m α/m β)ψ(x*). Eq. (20) indicates that the energy transfer rates increase linearly with the density and the average ionization and decrease with increasing piston velocity.
The timescale for energy dissipation is related to the collision mean free path. For ion–ion collisions the mean free path is (Huba, Reference Huba2011):
The average ionization Z av is obtained from the calculation of population of the ionization stages as a function of time:
where S z is the ionization coefficient for creating an ion with charge z, R z is the recombination coefficient of an ion with of charge z, $n_{\rm e} = \sum\nolimits_{z = 1}^Z {zn_z} $ is the electron density, and Z is the atomic number.
Analytical expressions for the ionization and recombination rates and values of the ionization potentials of aluminum are given in the Appendix.
The (Z + 1) equations for the ionization states (including the neutral) are solved together with the ion and electron temperatures as a function of time.
The formation of a shock wave requires that the spatial scale of the shocked region is larger than the spatial scale of the shock front width, setting a lower limit on the laser pulse duration τL
where the relaxation time τR is
Equivalent to Eq. (26), the spatial scale of the shocked region must be larger than the ion–ion mean free path:
Therefore, a lower limit on the laser pulse duration is:
Eqs. (10, 11, 24) are numerically integrated up to a time equal to the laser pulse duration.
Numerical results
The following values for the laser piston velocity are considered, 0.001c, 0.01c, and 0.1c. The corresponding shock velocity and compressibility, dimensionless pressure and laser intensity obtained from solving Eqs. 4–8 are given in Table 1.
The laser intensities, pressure and relaxation rate scales $\nu _0^{{\rm i}/{\rm i}} $ and $\nu _0^{{\rm i}/{\rm e}} $ for deuteron, carbon, and aluminum at normal density are given in Tables 2–4. The relaxation rate scales $\nu _0^{{\rm e}/{\rm e}} $ and $\nu _0^{{\rm e}/{\rm i}} $ are not shown, as $W_{\rm d}^{\rm i} \gg W_{\rm d}^{\rm e} $. The relaxation rate scales were calculated assuming that the plasma is fully ionized.
From Eq. (31) and Tables (2–4) it is seen that the formation of a collisional shock wave during the laser pulse duration is not possible for all the materials and the laser piston velocities considered here. From the energy loss rate scales one may approximate the relaxation time τR and obtain an estimate for the lower limit for the laser pulse duration. Table 2 shows that the relaxation time for liquid deuterium at piston velocities of 0.01c and 0.1c is of the order of hundreds of picoseconds and hundreds of nanoseconds, respectively. Therefore, the piston work cannot be transferred into thermal energy of the shocked material on a time scale of few picoseconds. Tables 3 and 4 show that there is no time for piston work deposition in carbon and aluminum at a piston velocity of 0.1c. In other words, the formation of a collisional shock wave is not possible under the above conditions. A lower limit for the laser pulse duration that may allow the formation of a collisional shock wave, resulting from Eq. (31) and estimated with relaxation time scales, is given in Table 5. This estimate is dependent on the material density and atomic number. The temperature-dependent part of the relaxation rates $\nu _\varepsilon ^{{\rm i}/{\rm i}} $ and $\nu _\varepsilon ^{{\rm i}/{\rm e}} $ [Eq. (19)] may shift the value for this lower limit.
The energy loss rates may be increased by the material pre-compression as seen from Eq. (20), and decrease the value of the lower limit for the pulse duration given in Table 5. However, an increase in the initial density will require higher laser intensity.
Solving the electron and the ion temperatures, Eqs. (10, 11) are not dependent on the laser pulse duration. The figures below show the results of the rate equations during a run time of 10 ps.
Figure 3 shows the results of Eqs. (10, 11) for deuterium at liquid normal density at piston velocity u p = 0.001c. The electron and ion temperatures as a function of time are shown in Figure 3a, the energy transfer rates [Eq. (21)] in Figure 3b, and the available energy for deposition W d − n ik BT i − n ek BT e in Figure 3c. It is seen that the energy deposition ends after about 40 fs, when the ion temperature is 0.48 keV and then the ions and the electrons reach equilibrium at about 0.35 keV after about 10 ps. Therefore, an equilibrium shock wave can be generated. For a piston velocity u p = 0.01c, only 13% of the piston work is deposited during a reference time duration of 10 ps, and the ions do not thermalize. Increasing the density to 10 g/cm3 enables thermalization of the ions after 3.6 ps to a temperature T imax = 43.8 keV, while the electrons are heated to T emax = 18.4 keV at this time. Therefore, a non-equilibrium shock wave is generated. Figure 4 shows the electron and ion temperatures and the energy transfer rates for initial density 10 g/cm3 and u p = 0.01c. After the piston work deposition, thermal energy is transferred from the ions to the electrons at a rate of 1/τeq [Eq. (14)].
Figure 5 shows the electron and ion temperatures and the energy transfer rates for carbon at initial density ρ0 = 1 g/cm3 and piston velocity u p = 0.01c. It is seen that a non-equilibrium shock wave is generated. The maximum ion temperature is about 250 keV at 4.6 ps when the piston work deposition ends. From the energy loss rates displayed in Figure 5b, the relaxation time τR is 6.7 ps [Eq. (30)].
Similar calculations for aluminum at normal density show the formation of an equilibrium shock wave at piston velocity u p = 0.001c, i.e. the electrons and the ions reach temperature equilibration during the shock wave, and a non-equilibrium shock wave at u p = 0.01c. Figure 6 shows the electron and ion temperatures, energy deposition rates, and piston work deposition as a function of time for normal density and piston velocity u p = 0.1c. In this case there is no time for collisional shock generation on a time scale of about 10 ps.
We note that to obtain a one-dimensional shock wave, the spot size must be larger than l s the spatial scale of the shocked region:
From Eq. (32) an estimate for the laser pulse energy is $E_{\rm L} = \pi r_{\rm L}^2 \tau _{\rm L}^3 I_{\rm L} \ge (\pi /36)u_{\rm p}^2 \tau _{\rm L}^3 I_{\rm L}$.
Summary
The formation of a collisional shock wave by the light pressure of a short-pulse laser at relativistic intensities and the subsequent heating of a solid layer are considered. The heating model is based on the dissipation of kinetic energy by binary Coulomb collisions at a temperature-dependent rate. We consider laser intensities in the range of 1018 W/cm2 to 1023 W/cm2 for deuterium, carbon, and aluminum targets. It is obtained that the formation of a collisional shock wave is dependent on the laser intensity and the material density. In the framework of the model presented here three cases were distinguished. If the timescale for energy dissipation is shorter than the laser pulse duration a collisional shock wave can be formed, during which the electrons and the ions either may have different temperatures – non-equilibrium shock wave, or may reach equilibrium at equal temperature. In both equilibrium and non-equilibrium shock waves, at early times the electrons are heated by energy transfer from the ions, and eventually the ion temperature overtakes. A third possibility may occur when the timescale for energy dissipation is longer than laser pulse duration, or the ion–ion collision mean free path l ii is larger than the shock width l s. In this case a collisional shock wave is not formed. The collisional shock wave described here may serve as a heating mechanism for a fast ignition scheme. In such a scheme deuterium fuel compressed by nanosecond laser pulses to a density by three orders of magnitude higher than normal liquid density can be heated by the a shock wave generated by the radiation pressure of 1 ps duration laser at several times its intensity of 1023 W/cm2.
Author ORCIDs
Zohar Henis, https://orcid.org/0000-0001-6813-1400.
Appendix: Ionization and recombination rates and ionization potentials
The following analytical expressions for the ionization and recombination rates were used (Eidmann, Reference Eidmann1994):
where I z[eV] is the ionization potential of an ion with charge z, Δz is the number of electrons in the last occupied shell of the ion with charge z.
The recombination rate is generally the contribution of three body and radiative recombination:
Here ξ z+1 is the number of vacancies in the last shell of ion z. In the above ionization and recombination rates T e is in eV units.
The (Z + 1) equations for the ionization states (including the neutral) are solved together with the ion and electron temperatures as a function of time.
The ionization potential in eV of the aluminum ion stages from the NIST Atomic Spectra database was used in the calculations which are: 5.98, 18.82, 28.44, 119.99, 153.82, 190.49, 241.76, 284.64, 330.21, 398.65, 442, 2085, and 2304.14.