INTRODUCTION
The generation of energetic particles and acceleration mechanisms in the ultra-intense laser pulses interaction with plasmas have become an international research focus nowadays (Chen et al., Reference Chen, Unick, Vafaei-Najafabadi, Tsui, Fedosejevs, Naseri, Masson-Laborde and Rozmus2008; Karmakar & Pukhov, Reference Karmakar and Pukhov2007; Schwoerer et al., Reference Schwoerer, Pfotenhauer, Jackel, Amthor, Liesfeld, Ziegler, Sauerbrey, Ledingham and Esirkepov2006; Flippo et al., Reference Flippo, Hegelich, Albright, Yin, Gautier, Letzring, Schollmeier, Schreiber, Schulze and Fernandez2007; Kaluza et al., Reference Kaluza, Schreiber, Santala, Tsakiris, Eidmann, Meyer-ter-Vehn and Witte2004; d'Humires et al., Reference d'Humires, Lefebvre, Gremillet and Malka2005; Wilks et al., Reference Wilks, Langdon, Cowan, Roth, Singh, Hatchett, Key, Pennington, MacKinnon and Snavely2001; Huang et al., Reference Huang, Lan, Duan, Tan, Wang, Shi, Tang and Xi2007, Reference Huang, Bi, Duan, Lan, Wang, Tang and He2008a, Reference Huang, Bi, Duan, Lan, Wang, Tang and He2008b). Some new schemes with novel targets (Limpouch et al., Reference Limpouch, Psikal, Andreev, Platonov and Kawata2008; Schwoerer et al., Reference Schwoerer, Pfotenhauer, Jackel, Amthor, Liesfeld, Ziegler, Sauerbrey, Ledingham and Esirkepov2006; Strangio et al., Reference Strangio, Caruso, Neely, Andreoli, Anzalone, Clarke, Cristofari, Del Prete, Di Giorgio, Murphy, Ricci, Stevens and Tolley2007; Schollmeier et al., Reference Schollmeier, Roth, Blazevic, Brambrink, Cobble, Fernandez, Flippo, Gautier, Habs, Harres, Hegelich, Hesslinga, Hoffmann, Letzring, Nurnberg, Schaumann, Schreiber and Witte2007) to enhance the ion acceleration and improve the quality of ion beams have been proposed. The progress of them can provide fundamental theory for inertial confined fusion (ICF) (Eliezer et al., Reference Eliezer, Murakaml and Val2007) and promote the realization of it effectively. The target normal sheath acceleration (TNSA) (Wilks et al., Reference Wilks, Langdon, Cowan, Roth, Singh, Hatchett, Key, Pennington, MacKinnon and Snavely2001) is a generally accepted mechanism for the ion acceleration. However, the fundamental theory of TNSA is a special self-similar solution for neutral-plasma isothermal expansions into a vacuum (Mora, Reference Mora2003; Huang et al., 2008a, 2008b), which is time-independent in the transformation system. Whatever, the influence of the hot-electron recirculation has not been accounted for in the self-similar solution. Therefore, it is suitable for the target thickness larger than a critical value, L c. Mackinnon et al. (Reference Mackinnon, Sentoku, Patel, Price, Hatchett, Key, Andersen, Snavely and Freeman2002) observed enhancement of proton acceleration by hot-electron recirculation in thin foils whose thickness is less than L c. In addition, Sentoku et al. (Reference Sentoku, Cowan, Kemp and Ruhl2003) predicted an equation to conclude the influence of electron recirculation. However, they did not propose a clear description of electron recirculation. Their physical picture is too simple and not clear because n times of electron recirculations cannot happen all at once and the electron density cannot jump to n times of the initial density. Although Huang et al. (Reference Huang, Lan, Duan, Tan, Wang, Shi, Tang and Xi2007) have proposed the step model to describe the hot-electron recirculation in which the electron density increased step by step, the step model combined the Mora's result and the linear increase of the electron density easily. However, an analytic solution on the bases of the fundamental equations: the equations of continuity and motion of the ions are still needed to reveal the influence of the hot-electron recirculation on maximum ion velocity rationally.
In this paper, a time-dependent solution for neutral-plasma isothermal expansions into a vacuum and a simple model based on the equations of continuity, and motion of the ions are proposed to obtain the time-dependent electron density and maximum ion velocity for the target of arbitrary thickness. The model describes the hot-electron recirculation and the influence of it on the ion acceleration in the laser-foil interactions analytically. The maximum ion velocity given by our model for some target thicknesses have been compared with experiments (Mackinnon et al., Reference Mackinnon, Sentoku, Patel, Price, Hatchett, Key, Andersen, Snavely and Freeman2002; Kaluza et al., Reference Kaluza, Schreiber, Santala, Tsakiris, Eidmann, Meyer-ter-Vehn and Witte2004). Since the decrease of the laser absorption efficiency and the three-dimensional (3D) effect on the electron density have been ignored, the results given by our model are larger than the experimental data for thin targets. In order to understand the differences in some degree, the influence of the decrease of the laser absorption efficiency and the 3D effect on the ion acceleration have been discussed.
TIME-DEPENDENT TARGET NORMAL SHEATH ACCELERATION
For convenience, the physical parameters are normalized as follows:
![\matrix{{\hat t = \omega _{\,pi0} t\comma \quad \hat l=l/\lambda _{D0}\comma \; } \cr{u=v/c_s\comma \quad \hat E=E/E_0\comma \quad \hat \phi=e\phi /k_B T_e }\comma \; \cr{\hat n=n/n_{e0}\comma \quad \hat c=c/c_s\comma \; } \cr}\eqno\lpar 1\rpar](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20160126040956937-0419:S0263034608000748_eqn1.gif?pub-status=live)
where E is the electric field, n (n i or n e) is the density (ion or electron), n e0 is the reference hot-electron density at the rear of the target, φ is the electric potential, is the ion acoustic speed,
is the initial ion plasma frequency,
, c is the light speed and
. Here e is the elemental charge.
When a relativistic laser pulse interacts with a solid target, the laser-produced fast electrons with a unique temperature, , determined by the laser ponderomotive potential are instantly created in front of the target and propagate through the target collisionlessly, and then form a high energy plasma at the rear of the target. The hot-electrons at the rear of the target can be considered to be reflected by sheath field at the ion front (Santos et al., Reference Santos, Amiranoff, Baton, Gremillet, Koenig, Martinolli, Rabec Le Gloahec, Rousseaux, Batani, Bernardinello, Greison and Hall2002; Mora, Reference Mora2003) and come back to the front of the target, because the field there is strongest. No matter, the laser pulse exists or has gone, the hot electrons are assumed to be reflected back and forth between the front of the target and the ion front at the rear of the target. Furthermore, the electron beam is assumed to be in equilibrium with a single temperature k BT e. Let the zero time be when the hot electrons come to the rear of the target for the first time and t 0 = 0. For convenience, the hot-electron speed used is the light speed c. Therefore, the time is −L/c when the laser get to the front target, where L is solid target thickness. The electrons return back from the rear at t = 0 and arrive at the front of the target when t = L/c. For t ≤ t l, assuming the laser intensity is constant, the hot-electron temperature,
, is invariant due to the energy supply from the laser pulse. The plasma expansion is an isothermal expansion. Let
be the time when the hot electrons undergo a circle and reach the ion front again at the first time.
The reference frame used here is . With the transformation, the equations of continuity and motion are obtained easily in the new coordinate system.
The fast-electron density is a function of the parameters: the acceleration time, τ, the target thickness, L, laser intensity, I, laser focus radius, r L, the laser absorption efficiency, η, the incidence angle of the laser pulse, θin, the half-opening angle of fast electrons, θe and amplified spontaneous emission (ASE) duration, . The laser absorption mechanisms in the femtosecond laser-plasma interactions determine the generation rate of the hot electrons, f, and then the laser absorption efficiency, η (Cai et al., Reference Cai, Yu, Zhu, Zheng, Cao and Pei2006). Assuming the time-dependent electron density for thin foils is variable-separating, then:
![\matrix{{n_e \lpar \tau\comma \; \xi \rpar =N\lpar \tau \rpar N_e \lpar \xi \rpar \comma \; N\lpar \tau=\tau _L \rpar =1\comma \; } \cr}\eqno\lpar 2\rpar](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20160126040956937-0419:S0263034608000748_eqn2.gif?pub-status=live)
where is the hot-electron density when
is a constant and N(τ) describes the increase of the maximum electron density due to electron recirculation and the electron generation by the laser-plasma interactions at the front of the target. The physical implication of β1 is the same as that given by Huang et al. (2008) and corresponds to
(where Z and M are the charge number and mass of the reference ion, Z k and M k are the charge number and mass of the concerned ion, respectively) in (Gurevich et al., Reference Gurevich, Anderson and Wilhelmsson1979). For different ions, β1 will be different. Therefore, the acceleration of ions of different masses and charges can be described with our model. If we choose the proton as the reference ion, β1 ≤ 1 for any other ions. Here,
corresponds to the special self-similar solution given by Huang et al. (2008a, 2008b) and Gurevich et al. (Reference Gurevich, Anderson and Wilhelmsson1979) for the impurity ions. β1 = 1 corresponds to the classic self-similar solution given by Mora (Reference Mora2003).
Assuming that f is constant in the interval, t l, the hot-electron total number that propagate through the target when for t l≤t L, and
for t l ≥ t L, where E l is the laser energy. Using Eq. (2) (Kaluza et al., Reference Kaluza, Schreiber, Santala, Tsakiris, Eidmann, Meyer-ter-Vehn and Witte2004), with N t, then n e0 in Eq. (2) can be estimated by:
![n_{e0}={{4.077\eta \lpar L\rpar I_{10^{18} {\rm W/cm^2} } } \over {\lpar \gamma - 1\rpar \lpar 1+\lpar L^{\ast} /r_L \rpar {\rm tan}\lpar \theta _{\rm e} \rpar \rpar ^2 }}\comma \; t_L \leq t_l\comma \; \eqno\lpar 3\rpar](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20160126040956937-0419:S0263034608000748_eqn3.gif?pub-status=live)
where r L is the laser pulse focus radius, is the efficient target thickness, θin is the incidence angle of the laser pulse, and θe ≈ 17° is half-opening angle of the superathermal electrons which was measured by Santos et al. (Reference Santos, Amiranoff, Baton, Gremillet, Koenig, Martinolli, Rabec Le Gloahec, Rousseaux, Batani, Bernardinello, Greison and Hall2002). When r L≫L and
, the angular effect can be neglected. Therefore, the influence of
and electron recirculation become dominated for thin targets.
As usual, assuming that the hot electrons satisfy uniform distribution in the bulk from the front of the target to the ion front and where
is the mean value of f, N(τ) can be roughly described by:
![N\lpar \tau \rpar ={\tau \over {\tau _L }}{{\hat L+\hat l\lpar \tau _L \rpar } \over {\hat l\lpar \tau \rpar +\hat L}}\comma \; \tau \leq \tau _l. \eqno\lpar 4\rpar](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20160126040956937-0419:S0263034608000748_eqn4.gif?pub-status=live)
For τL ≤ τl, the electron density increases continuously until t = t l. With Eq. (4), the maximum value of N(τ) is larger than 1, which can describe hot-electron recirculation (Mackinnon et al., Reference Mackinnon, Sentoku, Patel, Price, Hatchett, Key, Andersen, Snavely and Freeman2002; Sentoku et al., Reference Sentoku, Cowan, Kemp and Ruhl2003; Huang et al., Reference Huang, Lan, Duan, Tan, Wang, Shi, Tang and Xi2007). However, for τL ≥ τl, the maximum value of N(τ) is smaller than 1. Therefore, it is assumed conveniently that the hot-electron recirculation stop synchronously with laser pulse ending. As pointed by Huang et al. (2007), the critical condition for hot-electron recirculation is τL = τl and the critical target thickness is: L c = 0.5ct l. Therefore, hot-electron recirculation should be considered for L ≤ L c.
With the quasi-neutral approximation, . Solving the equations of continuity and motion of ions, the time-dependent ion velocity and potential are
, where
, and
. With this, neutral condition of the plasma is still satisfied. Therefore, the zero-ion-velocity position should be:
or
. For N ≡ 1, the result is the same as the self-similar solution's result (Huang et al., 2008a, 2008b). Therefore the dependence of the ion velocity on the time is:
![u_i=\hbox{ln}\lpar \tau ^{\beta _1 } N\rpar +\beta _1 \tau \partial \,\ln\,N\lpar \tau \rpar /\partial \tau. \eqno\lpar 5\rpar](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20160126040956937-0419:S0263034608000748_eqn5.gif?pub-status=live)
Eq. (5) shows the influence of the time-dependent electron density on the ion acceleration. Now the first part of Eq. (5) is ln(τN) compared with ln(τ) in the classic solution (Mora, Reference Mora2003). For N > 1, the hot-electron recirculation affects the ion acceleration. The second part exists for that means time-dependent. For N ≡ 1, Eq. (5) returns to the self-similar solution (Huang et al., 2008a, 2008b), which is time-independent.
The ion velocity at the ion front is the maximum ion velocity. Beyond the ion front, the ion density is zero and the electron density is still demonstrated by Eq. (2). The field for the ions at , where
is the self-similar variable at the ion front, is shielded by the ions beyond them. Therefore, the ion velocity at the ion front is larger than the ion velocity at
and can not be described by Eq. (5), and should be obtained by another way. There are two methods to obtain the ion-front velocity: a mathematical way used by Huang et al. (2008a, 2008b); a physical way given by Mora (Reference Mora2003). For convenience, the physical way will be used here. Similar with Mora's discussion (Mora, Reference Mora2003), using
, the ion velocity at the ion front is given by:
![u_{i\comma f}=\beta _1 \hbox{ln}\lpar \tau ^2 N\rpar +\beta _1 \tau \partial\, \ln N\lpar \tau \rpar /\partial \tau.\eqno\lpar 6\rpar](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20160126040956937-0419:S0263034608000748_eqn6.gif?pub-status=live)
Eq. (6) shows that the ion velocity at the ion front is larger than that given by Eq. (5) and .
Eqs. (4) and (6) are combined to give: , where
. For
at
,
is the position of the ion front. With
, the beginning time of the ion acceleration is:
, here
. Therefore, the position of the ion front is decided by:
![\eqalign{dF/d\tau &= \lpar \beta _1 F/\lpar F+\beta _1 \tau \rpar \rpar \lsqb \hbox{ln}\lpar \tau ^3 F_L /\tau _L F\rpar +1\rsqb \comma \; \tau \in \lsqb \tau _b\comma \; \tau _{acc} \rsqb .}\eqno\lpar 7\rpar](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20160126040956937-0419:S0263034608000748_eqn7.gif?pub-status=live)
where and
.
However, for or
(L c/L)(L+l f (τL)/(L+l f (τacc)))≈L c/L,
at
, where τacc is the normalized acceleration time and a little longer than the laser pulse duration. Therefore,
, where ηs is the saturate laser absorption efficiency for large target thickness and about 30–50% C is a constant and
is the physical acceleration time. The self-similar ion velocity at the ion front given by Huang et al. (2008a, 2008b) is
, which is suitable for the target thickness larger than L c since the hot-electron recirculation has not been accounted for. Sentoku et al.'s model (2003) shows that:
. However, the key physical relation given by Huang et al.'s (2007) step model is about
.
With above discussions, the dependence of laser absorption efficiency on the target thickness becomes the most important unknown information. With Eq. (3), for and
, the reference electron density
and is proportional to laser absorption efficiency and target thickness. When the laser absorption efficiency tends to zero with
,
and the maximum ion velocity tends to zero too. d'Humires et al. (Reference d'Humires, Lefebvre, Gremillet and Malka2005) have given a curve to describe that, but can it be used for any other target or other parameters of laser pulse? However, the changing law should be similar, the exact relationship have not been obtained. In fact, the dependence of opening angle of electrons on the target thickness and laser parameters is also important for the electron density, and the opening angle of ions to be ascertained, but still a challenge. The angle decides the electron density n e0 as shown in Eq. (3) and increases with target thickness. Therefore, for thin target, the opening angle should be smaller than 17° (Santos et al., Reference Santos, Amiranoff, Baton, Gremillet, Koenig, Martinolli, Rabec Le Gloahec, Rousseaux, Batani, Bernardinello, Greison and Hall2002) and becomes zero when
.
COMPARISON WITH EXPERIMENTS AND DISCUSSION
In fact, to compare our model with the relative experiments exactly, and
should be ascertained first. However, Figure 1 shows a course comparison with the assumptions: laser absorption efficiency and the opening angle of electrons are fixed for different target thicknesses, η = 40%, θe = 17°, and β1 = 1. Other parameters are the same as the compared experiments. With Figure 1a, the results of our model are consistent with the experimental data. In these experiments, Mackinnon et al. (Reference Mackinnon, Sentoku, Patel, Price, Hatchett, Key, Andersen, Snavely and Freeman2002), the contrast ratio of intensity is large enough that the influence of the prepulse can be neglected. Therefore, for L ≳ 1 µm, the laser absorption efficiency sustains a constant. However, if the contrast ratio is not large enough and the influence of the prepulse on η can not be ignored, η should not be assumed to be a constant and the results given by our model are larger than experimental data as shown by Figure 1b. It is concluded that for L ≤ L c, the influence of the hot-electron recirculation on the ion acceleration is evident: maximum ion energy increases quickly with the decrease of the target thickness and larger than the experimental data.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160627145026-84505-mediumThumb-S0263034608000748_fig1g.jpg?pub-status=live)
Fig. 1. (Color online) Comparisons with Kaluza et al.'s and Mackinnon et al.'s experiments, where the unit of I is 1018 W/cm2. Kaluza et al.'s experimental parameters are: λ = 790 nm, r L = 2.5 µm, t l = 150 fs and t ASE = 2.5 ns, I = 1.0 × 1019 W/cm2, I = 1.3 × 1019 W/cm2 and I = 1.5 × 1019 W/cm2 with incidence of 30° onto targets, and laser absorption efficiency is about 40%; Mackinnon et al.'s experimental parameters are: λ = 0.8 µm, r L = 2.5 µm, t l = 100 fs, intensity contrast ratio is 1010: 1 with incidence of 22° onto targets, and laser absorption efficiency is 40%. In our model, the half opening angle of electrons is 17°(Santos et al., Reference Santos, Amiranoff, Baton, Gremillet, Koenig, Martinolli, Rabec Le Gloahec, Rousseaux, Batani, Bernardinello, Greison and Hall2002).
The reasons for larger results given by our model may be:
(1) the parameters of laser absorption efficiency is taken as constants. In fact, the laser absorption efficiency stays constant for the target of thickness large enough and decreases with the decrease of the target thickness for thin foils and tends to zero as shown by Fig. 12 in (d'Humires et al., Reference d'Humires, Lefebvre, Gremillet and Malka2005) with particle-in-cell simulations. In general cases, the analytic function of the dependence of laser absorption efficiency on target thickness, laser intensity, laser pulse duration and ASE duration is still a challenge, although Huang et al. has a half-analytic way to obtain it (Huang et al., Reference Huang, Lan, Duan, Tan, Wang, Shi, Tang and Xi2007). Here, the function can also be concluded with the similar method. It is easy to realized and is not repeated here.
(2) the three dimensional effect is not contained. The electron density does decrease with the plasma expansion as the electrons recirculate and spread laterally with a given opening angle. This dilute effect of the electron density is ignored here. In our sequence articles, it will be considered.
CONCLUSION
In conclusion, the time-dependent isothermal expansion for the target normal sheath proton acceleration is briefly proposed. The influence of hot-electron recirculation on the ion acceleration has been shown by Eq. (6) analytically and by Figure 1 obviously. For L ≤ L c, the hot-electron recirculation enhances the ability of the ion acceleration for thin foils.
ACKNOWLEDGMENTS
This work was supported by the Key Project of Chinese National Programs for Fundamental Research (973 Program) under contract No. 2006CB806004 and the Chinese National Natural Science Foundation under contract No. 10334110.