Hostname: page-component-cd9895bd7-dk4vv Total loading time: 0 Render date: 2024-12-23T20:16:37.281Z Has data issue: false hasContentIssue false

Double layer acceleration by laser radiation

Published online by Cambridge University Press:  14 February 2014

Shalom Eliezer*
Affiliation:
Institute of Nuclear Fusion, Polytechnic University of Madrid, Madrid, Spain Soreq Research Center, Yavne, Israel
Noaz Nissim
Affiliation:
Soreq Research Center, Yavne, Israel
José Maria Martínez Val
Affiliation:
Institute of Nuclear Fusion, Polytechnic University of Madrid, Madrid, Spain
Kunioki Mima
Affiliation:
Institute of Laser Engineering, Osaka University, Osaka, Japan
Heinrich Hora
Affiliation:
Department of Theoretical Physics, University of New South Wales, Sydney, Australia
*
Address correspondence and reprint requests to: S. Eliezer, Institute of Nuclear Fusion, Polytechnic University of Madrid, Madrid, Spain. E-mail: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

It is shown that it is possible to accelerate micro-foils to velocities from 108 cm/s up to relativistic velocities without the disturbance of the Rayleigh-Taylor instability. The acceleration occurs due to the radiation pressure of proper high power lasers. In these systems, the ablation force is negligible relative to the ponderomotive force that dominates the acceleration. The laser irradiances of 1017 W/cm2 < IL < 1021 W/cm2 with a pulse duration of the order of 10 picoseconds can accelerate a micro-foil by the laser radiation pressure to velocities as high as 109 cm/s before breaking by Rayleigh Taylor (RT) instability. Similarly, laser irradiances of IL > 1021 W/cm2 with pulse duration of the order of 10 femtoseconds can accelerate a micro-foil to relativistic velocities without RT breaking. Due to the nature of the accelerating ponderomotive force, in both the relativistic and non-relativistic cases, the structure of the accelerated target contains a double layer (DL) at the interface of the laser-target interaction. The DL acts as a piston during the acceleration process. The influence of the DL surface tension on the RT instability is also analyzed in this paper.

Type
Research Article
Copyright
Copyright © Cambridge University Press 2014 

1. INTRODUCTION

The use of spherical shells in inertial confinement fusion, in comparison with solid spheres, has reduced significantly the peak laser power and energy required to ignite a target of a given mass. The ablation pressure P a that accelerates the shell is a function of the irradiating laser intensity I L which scales approximately like P a ~ I L2/3. For a shell with initial radius R 0 and thickness ΔR 0, the implosion velocity v f scales as v f2 ~P a(R 0R 0). The resulting stagnation pressure is P f ~ P a(R 0R 0)G, where G is a radial convergence factor of the order of 100. The minimum energy of ignition scales as 1/P f2. In order to get ignition, the ablation pressure of about 30 Mbars has to increase to P f ~ 100 Gbars, i.e., by a factor of 3000, implying an aspect ratio R 0R 0 ~ 30. The larger the aspect ratio is the higher the implosion pressure (P f) one gets implying smaller laser energy (~1/ P f2). However, the magnitude of the aspect ratio is limited by Rayleigh Taylor (RT) instability (Rayleigh, Reference Rayleigh1876; Taylor, Reference Taylor1950). The RT instability can cause a disturbance to grow from extremely small amplitudes to large levels that can break the accelerated foil or the shell before it reaches the desired final velocity. Since the beginning of the inertial confinement fusion program the RT instability has been a subject of uncertainties and controversies.

The highest laser-induced pressures, of the order of 1Gbar, have been obtained during the collision of a target with an accelerating foil. This acceleration was achieved by the ablation of laser produced plasma or by X-rays created by laser plasma interaction. The foil absorbs the laser energy, plasma is then obtained and the foil is accelerated like a rocket (Eliezer, Reference Eliezer2002). In this way, the flyer stores kinetic energy from the laser during the acceleration period and delivers it during the collision with a target, in the form of thermal energy. The laser driven flyer can achieve much higher pressures on impact than the shock wave pressure obtained directly by the laser ablation. In 1994, at the Livermore Laboratory, a pressure of about 1Gbar was created by accelerating a foil with soft X-rays from the indirect drive (Cauble et al., Reference Cauble, Phillion, Hoover, Holmes, Kilkenny and Lee1993). Planar foil targets were accelerated to 1000 km/s at the Institute of Laser Engineering at Osaka University (Azechi et al., Reference Azechi, Sakaiya, Watari, Karasik, Saito, Ohtani, Takeda, Hosoda, Shiraga, Nakai, Shigemori, Fujioka, Murakami, Nagatomo, Johzaki, Gardner, Colombant, Bates, Velikovich, Aglitskiy, Weaver, Obenschain, Eliezer, Kodama, Norimatsu, Fujita, Mima and Kan2009). The targets were about 20 µm thick and were made of polystyrene doped with 0.4% (atomic) bromine. The acceleration driver was a 0.35 µm wavelength laser with energy of 1.5 kJ, 2 ns pulse width and an irradiation of 4 × 1014 W/cm2.

For I L < 1016 W/cm2, the ablation pressure is dominant. for I L ≫  1016 W/cm2, the radiation pressure is the dominant pressure at the solid-vacuum interface and the ablation pressure is negligible. For I L > 1021 W/cm2, it was suggested (Esirkepov et al., Reference Esirkepov, Borghesi, Bulanov, Mourou and Tajima2004; Eliezer Reference Eliezer2012) that relativistic acceleration can be achieved.

We are interested in the classical acceleration of a micro-foil that can be achieved for laser irradiances in the domain 1017 W/cm2 < I L < 1021 W/cm2 with laser pulses of the order of picosecond where the ablation pressure can be neglected relative to the radiation pressure. We also discuss the relativistic acceleration that can be achieved with I L > 1021 W/cm2.

This paper analyzes the foil acceleration by ponderomotive force. In this model, the force acts on the electrons that are accelerated and the ions that follow accordingly. The system of the negative and positive layers is called a double layer (DL). The DL is regions of non-neutral plasma which induce a large potential drop thereby causing very strong electric fields. In our case, we have DL acceleration by the laser irradiance (Clark et al., Reference Clark, Eliezer and Farley1985; Eliezer & Hora, Reference Eliezer and Hora1989; Hora et al., Reference Hora, Lalousis and Eliezer1984).

In Section 2, the classical laser pressure acceleration and the hydrodynamic Rayleigh-Taylor (RT) instability are calculated relevant to laser irradiances in the domain 1017 W/cm2 < I L < 1021 W/cm2. Section 3 analyzes the surface tension contribution to the RT instability. The relativistic hydrodynamic RT instability appropriate for I L ≥ 1021 W/cm2 is summarized in Section 4. The analysis is discussed in the last section.

2. LASER PRESSURE ACCELERATION AND RT INSTABILITY

The radiation pressure P r created by a laser irradiance I L at the solid-vacuum interface is given by

(1)$$P_r=\left({\displaystyle{{I_L } \over c}} \right)\left({1+R} \right).$$

R is the laser reflection from the micro-foil target and c is the light velocity. For the following calculations, we take R = 1 that might be achieved with very short laser pulses τL where τL ~ 10 picoseconds and less. In this case, for a foil where the center of mass is moving in the x direction, one gets from the Newton second law that the acceleration a of the center of mass, in the x-direction, is given by

(2)$$a=\displaystyle{{2I_L } \over {{\rm \rho} cl_x }} \equiv - g\comma \;$$

ρ is the foil density and the x,y,z dimensions of the foil are l x, l y, l z accordingly. In this case, one can analyze the foil fluid motion in a gravitational field described by g. For a fluid motion described by its field velocities ν(x,y,z) without turbulence, ∇× ν = 0, one can define the velocity potential ϕ by the relation ν = ϕ(x,y,z;t). Neglecting the quadratic velocity term (ν2) and assuming a non-compressible fluid one can write the Euler equation in a gravitational field in the following way

(3)$$\displaystyle{{\partial {\rm \phi} } \over {\partial t}}+\displaystyle{P \over {\rm \rho} }+g{\rm x}= 0.$$

Using the relations νx = ∂x/∂t = ∂ϕ/∂x, v y = ∂ϕ/∂y, v z = ∂ϕ/∂z, taking the local pressure P as induced by the surface tension, P = −α(∂2x/∂y 2 + ∂2x/∂z 2), and taking the time derivative of Eq. (3) we get the following equation

(4)$${\rm \rho} g\displaystyle{{\partial {\rm \phi} } \over {\partial x}} + {\rm \rho} \displaystyle{{\partial ^2 {\rm \phi} } \over {\partial t^2 }} - {\rm \alpha} \displaystyle{\partial \over {\partial x}}\left({\displaystyle{{\partial ^2 {\rm \phi} } \over {\partial y^2 }} + \displaystyle{{\partial ^2 {\rm \phi} } \over {\partial z^2 }}} \right)= 0 .$$

Defining the foil in the domain: −l x ≤ x ≤ 0, 0 ≤ y ≤ l y, 0 ≤ z ≤ l z, and assuming the boundary conditions

(5)$$\eqalign{& {\rm v}_x \left({x = - l_x } \right)= \displaystyle{{\partial {\rm \phi} } \over {\partial x}}\left({x = - l_x } \right)= 0 \cr & {\rm v}_y \left({y = 0\comma \; l_y } \right)= \displaystyle{{\partial {\rm \phi} } \over {\partial y}}\left({y = 0} \right)= \displaystyle{{\partial {\rm \phi} } \over {\partial y}}\left({y = l_y } \right)= 0 \comma \cr & {\rm v}_z \left({z = 0\comma \; l_z } \right)= \displaystyle{{\partial {\rm \phi} } \over {\partial z}}\left({z = 0} \right)= \displaystyle{{\partial {\rm \phi} } \over {\partial z}}\left({z = l_z } \right)= 0}$$

we get the following solution of Eq. (4)

(6)$$\eqalign{& {\rm \phi} \left({x\comma \; y\comma \; z\semicolon \; t} \right)= e^{ - i{\rm \omega} t} Ae^{ - kl_x } \left[{e^{k\left({x + l_x } \right)} + e^{ - k\left({x + l_x } \right)} } \right]\;f\left({y\comma \; z} \right)\cr & f\left({y\comma \; z} \right)= \cos \left({\displaystyle{{n_y y} \over {l_y }}} \right)\cos \left({\displaystyle{{n_z z} \over {l_z }}} \right)\semicolon \; {\rm }\left. \matrix{n_y \hfill \cr n_z \hfill} \right\}{\rm = 0\comma \; 1}{\rm \pi} {\rm\comma \; 2}{\rm \pi} {\rm\comma \; 3}{\rm \pi} {\rm\comma \; }...{\rm }}$$

Substituting this solution into Eq. (4) we get the following dispersion relation at the foil back surface x = 0 (the laser irradiates the foil at the front surface, x = −l x, while we are interested in the back surface instability about the surface x = 0).

(7)$${\rm \omega} ^2 = - k\left[{a \cdot {\rm tgh}\left({kl_x } \right)- \displaystyle{{\rm \alpha} \over {\rm \rho} }\left({\displaystyle{{n_y ^2 } \over {l_y ^2 }}+\displaystyle{{n_z ^2 } \over {l_z ^2 }}} \right)} \right].$$

In Eq. (7), we have substituted the acceleration a instead of the gravitational field g (g = −a, a > 0 according to Eq. (2)).

Denoting by ξ the local disturbance at the back surface of the foil we obtain from Eq. (6)

(8)$${\rm \xi} = \vint {\left({\displaystyle{{\partial {\rm \phi} } \over {\partial x}}} \right)_{x = 0} dt = } \;{\rm \xi} _0 e^{ - i{\rm \omega} t}$$

where ξ0 is the initial disturbance at the x = 0 surface and ω is given by Eq. (7). Neglecting the surface tension term (i.e., α = 0), ω is imaginary and one gets the RT instability by using Eq. (2)

(9)$$\eqalign{& \displaystyle{{\rm \xi} \over {{\rm \xi} _0 }} = {\rm exp }\lcub \lsqb {\rm ka} \cdot {\rm tgh}\left({{\rm k}l_{\rm x} } \right)\rsqb ^{1/2} {\rm t}\rcub \equiv \exp \left({\displaystyle{t \over {{\rm \tau} _{RT} }}} \right)\cr & \displaystyle{1 \over {{\rm \tau} _{RT} }} \equiv \left[{\left({\displaystyle{{2{\rm \pi} a} \over L}} \right){\rm tgh}\left({\displaystyle{{2{\rm \pi} l} \over L}} \right)} \right]^{1/2} = \left[{\left({\displaystyle{{4{\rm \pi} I_L } \over {{\rm \rho} clL}}} \right){\rm tgh}\left({\displaystyle{{2{\rm \pi} l} \over L}} \right)} \right]^{1/2} } .$$

In the second term of Eq. (9), the most restrictive constrain on the instability is taken by choosing k = 2π/L, where L is the target dimension orthogonal to the x amplitude. For tgh(kl x) = 1 one gets the standard RT expression (Eliezer, Reference Eliezer2002; Chandrasekhar, Reference Chandrasekhar1961). However in our case L ~ l y = l z ≫ l xl implying tgh(2πl/L) ~ 2πl/L and the following time scale for the RT instability

(10)$${\rm }{\rm \tau} _{RT} = \left({\displaystyle{L \over {2{\rm \pi} }}} \right)\sqrt {\displaystyle{{{\rm \rho} c} \over {2I_L }}} \;{\rm for }\;l{\rm\ll \; L}.$$

Denoting the critical value that the foil breaks by (ξ/ξ0)cr at time t cr, one gets the foil velocity v f before breakdown by RT instability

(11)$$\eqalign{& {\rm v}_f = at_{cr} = a{\rm \tau} _{RT} \ln \left({\displaystyle{{\rm \xi} \over {{\rm \xi} _0 }}} \right)_{cr} \Rightarrow \cr & {\rm v}_f = c\left\{{\left({\displaystyle{{I_L } \over {{\rm \rho} c^3 }}} \right)\left({\displaystyle{L \over {{\rm \pi} l}}} \right)\left[{tgh\left({\displaystyle{{2{\rm \pi} l} \over L}} \right)} \right]^{ - 1} } \right\}^{1/2} ln\left({\displaystyle{{\rm \xi} \over {{\rm \xi} _0 }}} \right)_{cr}. }$$

For I L = 1018 W/cm2, L =10 µm, l = 1 µm, and ρ = 1 g/cm3 we get using Eqs. (2) and (9) the acceleration a = (2/3) × 1019 cm/s2 and the RT time scale τRT =6.54 ps. This interesting result implies that a foil with very small initial disturbances can be accelerated until it breaks at ln(ξ/ξ0)cr ~3 inferring t cr ~20 ps before breaking down by RT instability. In this case, the foil velocity v f is about v f ~ 3aτRT = 1.34 × 108 cm/s.

The final velocity of the accelerated foil before breakdown by Rayleigh Taylor instability according to Eq. (11) is given in Figure 1, where we have chosen ln(ξ/ξ0)cr = 3 and the foil density ρ = 1 g/cm3. The ratio of foil thickness (l) to the transversal dimension (L) are given for (1) 2πl/L= 1, (2) 2πl/L=0.1. As one can see from these figures the laser can accelerate foils to velocities as high as 109 cm/s before breaking by RT instability.

Fig. 1. (Color online) The final velocity of the accelerated foil before breakdown by Rayleigh Taylor instability according to Eq. (11). In this equation ln(ξ/ξ0)cr was chosen to be equal to 3 and the foil density ρ =1 g/cm3. The ratio of foil thickness (l) to the transversal dimension (L) are chosen for (a) 2πl/L = 1, (b) 2πl/L = 0.1.

3. THE SURFACE TENSION AND RAYLEIGH-TAYLOR (RT) INSTABILITY

In this section, we describe the double layer by a capacitor model as given in Figure 2. This capacitor model is controlled by laser irradiances I L where the ponderomotive force dominates the interaction. E x is the electric field inside the double layer, λDL is the distance between the positive and negative DL charges, l is the foil thickness, and δ is the solid density skin depth of the foil. n e and n i are the electron and ion densities accordingly, E x is the electric field solution. The DL is geometrically followed by neutral plasma where the electric field decays within a skin depth δ given by δ = c/ωpe, where ωpe = 5.64 × 104n e = 1.78 × 1016 s−1, implying δ ~ 2 × 10−6 cm.

Fig. 2. (Color online) (a) The capacitor model for laser irradiances I L where the ponderomotive force dominates the interaction. E x is the electric field inside the double layer, λDL is the distance between the positive and negative DL charges, l is the foil thickness and δ is the solid density skin depth of the foil. (b) A schematic figure that our capacitor model is based. n e and n i are the electron and ion densities accordingly, E x is the electric field solution. The DL is geometrically followed by neutral plasma where the electric field decays within a skin depth.

The acceleration of double layers by a high power laser to high velocity was suggested in the literature for a long time. The dispersion relation that takes into account also the double layer (DL) surface tension (Eliezer & Hora, Reference Eliezer and Hora1989; Hora et al., Reference Hora, Min, Eliezer, Lalousis, Pease and Szichman1989) is given in Eq. (7). The double layer motion is RT stable if ω2 > 0 and this can be achieved if the following inequality is satisfied:

(12)$${\rm \alpha} \geq {\rm \rho} aL^2 \cdot {\rm tgh}\left({\displaystyle{{2{\rm \pi} l} \over L}} \right).$$

In Eq. (12), we have taken the most unstable modes as defined in Eq. (7) by (n y = 2π, n z =0) or (n y =0, n z = 2π). The DL surface tension within the context of a capacitor model (see Fig. 2) is given by

(13)$${\rm \alpha} = \left({\displaystyle{{E_{DL} ^2 } \over {8{\rm \pi} }}} \right){\rm \lambda} _{DL} \approx \displaystyle{{I_L {\rm \lambda} _{DL} } \over {c{\rm \varepsilon} _r ^2 }}$$

E DL is the electric field in the DL region; λDL is the distance between the positive and negative DL charges, usually (Eliezer & Hora, Reference Eliezer and Hora1989) of the order of the Debye length λD ≈ 750(T eV/n e)1/2 where T eV and n e are the electron temperature in eV units and the electron density accordingly. εr ≈ (1 − ωp2L2) = (1 − n e/n c) is the dielectric constant in the DL domain. Since the electrons are pushed out of the DL by the ponderomotive force then one can assume that the electron density inside the DL is smaller than the critical density n c, described schematically in Figure 2, implying εr ≈ 1 as a good approximation. Substituting the acceleration from Eq. (2) into Eq. (13) we get the DL stability for

(14)$${\rm \lambda} _{DL} \geq 2\left({\displaystyle{{L^2 } \over l}} \right){\rm tgh}\left({\displaystyle{{2{\rm \pi} l} \over L}} \right).$$

For l/L ≪  1 we can approximate the DL stability by the simple relation

(15)$${\rm \lambda} _{DL} \geq 4{\rm \pi} L .$$

As a typical example we can take λDL of the order of the Debye length λD and for T = 100 eV, n e = 1022 cm−3 we have λDL ~ 7.5 × 10−8 cm ~10 L, a transversal dimension of a nano-foil or rather a nanotube! Therefore it looks that for micro-foils the DL surface tension does not contribute to hydrodynamic stability unless the electron density is few orders of magnitude smaller than 1022 cm−3 and the temperature significantly higher than 100 eV.

4. RELATIVISTIC RAYLEIGH-TAYLOR (RT) INSTABILITY

With the development (Piazza et al., Reference Piazza, Muller, Hatsagortsyan and Keitel2012) of laser irradiances larger than 1021 W/cm2 it became possible to accelerate a micro-foil to relativistic velocities due to the laser pressure P =2I F /c, where I F is the laser irradiance in the accelerated foil rest frame and c is the speed of light and it is assumed that the reflection in the foil rest frame is one. The laser irradiance I L and angular frequency ωL in the laboratory are related to their appropriate values in the foil rest frame (denoted by subscript F) by Doppler equation, based on appropriate Lorenz transformation of the electro-magnetic fields and the definition of I in term of these fields,

(16)$$\eqalign{& I_L = I_F \left({\displaystyle{{{\rm \omega} _L } \over {{\rm \omega} _F }}} \right)^2 = I_F \left({\displaystyle{{1 + {\rm \beta} _f } \over {1 - {\rm \beta} _f }}} \right)\cr & P = \displaystyle{{2I_L } \over c}\left({\displaystyle{{1 - {\rm \beta} _f } \over {1+{\rm \beta} _f }}} \right)}$$

where cβf is the foil velocity in the laboratory frame of reference.

The relativistic Newton law of motion

(17)$$\displaystyle{d \over {dt}}\left[{\left({{\rm \rho} _0 lc} \right)\displaystyle{{{\rm \beta} _f } \over {\sqrt {1 - {\rm \beta} _f ^2 } }}} \right]= \displaystyle{{2I_L } \over c}\left({\displaystyle{{1 - {\rm \beta} _f } \over {1+{\rm \beta} _f }}} \right)\comma \;$$

has been solved (Eliezer et al., Reference Eliezer, Martinez Val and Pinhasi2013)

(18)$$\mathop{\vint}\limits_0^{{\rm \beta} _f} \displaystyle{dx \over \left({1 - x} \right)^{5/2} \left({1 + x} \right)^{1/2}} = \displaystyle{\left(2 - {\rm \beta} _f \right)\sqrt{1 - {\rm \beta} _f ^2} \over 3\left(1 - {\rm \beta}_f \right)^2} - \displaystyle{2 \over 3} = \displaystyle{2It \over {\rm \rho}_0 c^2 l} \equiv \displaystyle{t \over {\rm \tau} } \comma \;$$

for a pulsed laser with constant intensity in order to get micro-foil velocity βf as function of the laser pulse duration t L in units of τ = ρ0c 2l/(2I L), where ρ0 is the micro-foil initial density and l is the foil thickness. Figure 3 describes the micro-foil velocity as a function of the laser pulse duration in units of τ. For example, if ρ0 = 1 g/cm3, l = 0.1 µm, foil cross-section of 10 µm2, then βf = 0.5 for a laser with I Lt L = 4.5 × 108 J/cm2 and energy 45 J. Furthermore, for an irradiance of 1022 W/cm2 one requires a laser with pulse duration of 45 fs.

Fig. 3. (Color online) Micro-foil velocity βfc as a function of the laser pulse duration in units of τ = ρ0c 2l/(2I L), where ρ0 is the micro-foil initial velocity upon impact and l is the foil thickness.

A crucial question in accelerating a foil to relativistic velocity is its hydrodynamic stability. In particular, the relativistic Rayleigh-Taylor instability was calculated (Pegoraro & Bulanov, Reference Pegoraro and Bulanov2007) as described in Eq. (19).

(19)$$\eqalign{& {\rm \xi} _{NR} = \displaystyle{{\Delta x} \over {x_0 }} = \exp \left({\displaystyle{t \over {{\rm \tau} {}_{NR}}}} \right)\cr & {\rm \xi} _R = \displaystyle{{\Delta x} \over {x_0 }} = \exp \left[{\left({\displaystyle{t \over {{\rm \tau} {}_R}}} \right)^{1/3} } \right]\cr & \displaystyle{{{\rm \tau} _R } \over {{\rm \tau} _{NR} }} = \left({\displaystyle{1 \over {3{\rm \pi} }}} \right)\left({\displaystyle{L \over {l_0 }}} \right)\left({\displaystyle{{I_L } \over {{\rm \rho} _0 c^3 }}} \right)}$$

ξNR and ξR are accordingly the non-relativistic and relativistic development of the instability for an initial disturbance x 0 and L is the target dimension orthogonal to the x amplitude. The non-relativistic time scale τNR we calculate from the second of Eq. (9): τNRRT. We consider the following example: L =10 µm, x 0 =10 nm, ρ0 =1 g/cm3, l = 0.1 µm, I L =1024 W/cm2. Assuming that the foil breaks for ξ = 10 (i.e., Δx ~ l) then the foil breaks at 14.2 fs (τNR = 6.2 fs) for the non-relativistic case while in the relativistic regime the foil is stable during 90.5 fs (τR =39.3 fs). This behavior is understood from the different time dependence of the RT instability in relativistic and non-relativistic cases. Taking into account that during few femtoseconds one can accelerate a micro-foil to relativistic velocities one can achieve stable relativistic acceleration for laser irradiances of the order of 1023 W/cm2.

5. DISCUSSION

In this paper, we analyze the acceleration of a micro-foil by the laser radiation pressure impact. We look for acceleration without laser-plasma ablation. The laser irradiances of interest are in the domain 1017 W/cm2 < I L < 1021 W/cm2 for non-relativistic acceleration and I L ≥ 1021 W/cm2 for the relativistic regime. Due to the nature of the ponderomotive force, in both the relativistic and non-relativistic cases, the structure of the accelerated target contains a double layer as shown in Figure 2. Furthermore, for ultra-high irradiance lasers the RT instability does not apparently disturbs the acceleration as has been derived theoretically in particle in cell simulations (Esirkepov et al., Reference Esirkepov, Borghesi, Bulanov, Mourou and Tajima2004) of a foil and two fluid hydrodynamic simulations of plasma blocks (Lalousis et al., Reference Lalousis, Hora, Eliezer, Martinez Val, Moustaizis, Miley and Mourou2013; Reference Laluosis, Foldes and Hora2012; Hora et al., Reference Hora, Lalousis and Moustaizis2014).

The main result of this research is that the laser radiation pressure can accelerate a micro-foil to very high velocities without breaking it by the RT instabilities. For the non-relativistic domain the laser time duration is of the order of τL~ 10 picoseconds while in the relativistic case the laser pulse is τL~ 10 femtoseconds.

The surface tension is also analyzed using a capacitor model in the context of hydrodynamic stability. It is shown that for solid state density the surface tension cannot stabilize a micro-foil but it might cause stabilization for nano-foils not discussed in this paper (Eliezer et al., Reference Eliezer, Eliaz and Grossman2004; Reference Eliezer, Eliaz, Grossman, Fisher, Gouzman, Henis, Pecker, Horovitz, Fraenkel, Maman and Lereah2005). However, for under dense DL plasmas, the surface tension can be an important stabilizing mechanism.

The non-relativistic velocities that can be achieved are given explicitly in Eq. (11) and described in Figures 1. As one can see from these figures the laser can accelerate foils to velocities as high as 109 cm/s before breaking by RT instability. The relativistic velocities are specified implicitly in Eq. (18) as designated in Figure 3. Therefore, for the state of the art high irradiance lasers the RT instability can be avoided during the acceleration time and relativistic velocities can be achieved.

Finally, we would like to suggest that the consequences of this paper may be important in discussing the possibility of fast ignition (Basov et al., Reference Basov, Guskov and Feoktistov1992; Tabak et al., Reference Tabak, Hammer, Glinsky, Kruer, Wilks, Woodworth, Campbell, Perry and Mason1994) nuclear fusion triggered by laser or by foil impact induced shockwaves (Betti et al., Reference Betti, Zhou, Anderson, Perkins, Theobald and Sokolov2007; Eliezer & Martinez Val, Reference Eliezer and Martinez Val2011; Murakami et al., Reference Murakami, Nagatomo, Azechi, Ogando and Eliezer2006; Eliezer et al., Reference Eliezer, Murakami and Martinez Val2007). The "fast ignition" idea is to ignite the fuel (e.g., deuterium- tritium, proton-boron11, etc.) with less energy and in order to achieve this goal one has to separate the drivers that compress and ignite the target. First the fuel is compressed then a second driver ignites a small part of the fuel while the created alpha particles heat the rest of the target. The fast ignition problem is that the laser pulse does not penetrate directly into the compressed target; therefore many schemes have been suggested. In this context the comeback of shock waves in inertial fusion energy was published and the ignition by impact was suggested (Eliezer, Reference Eliezer2012; Murakami et al., Reference Murakami, Nagatomo, Azechi, Ogando and Eliezer2006). The RT instability plays a crucial role in achieving ignition (Mima et al., Reference Mima, Murakami, Nakai and Eliezer2009) and therefore the analysis of this paper might be important.

References

REFERENCES

Azechi, H., Sakaiya, T., Watari, T., Karasik, M., Saito, H., Ohtani, K., Takeda, K., Hosoda, H., Shiraga, H., Nakai, M., Shigemori, K., Fujioka, S., Murakami, M., Nagatomo, H., Johzaki, T., Gardner, J., Colombant, D.G., Bates, J.W., Velikovich, A.L., Aglitskiy, Y., Weaver, J., Obenschain, S., Eliezer, S., Kodama, R., Norimatsu, T., Fujita, H., Mima, K. & Kan, H., (2009). Experimental evidence of impact ignition: 100-fold increase of neutron yield by impactor collision. Phys. Rev. Lett. 102, 235002.Google Scholar
Basov, N.G., Guskov, S.Y. & Feoktistov, L.P. (1992). Thermonuclear gain of ICF targets with direct heating of the ignitor. J. Soviet Laser Res. 13, 396.Google Scholar
Betti, R., Zhou, C.D., Anderson, K.S., Perkins, L.J., Theobald, W. & Sokolov, A. A. (2007). Shock ignition of thermonuclear fuel with high areal density. Phys. Rev. Lett. 98, 155001.Google Scholar
Cauble, R., Phillion, D.W., Hoover, T.J., Holmes, N.C., Kilkenny, J.D. & Lee, R.W. (1993). Demonstration of 0.75 Gbar planar shocks in x-ray driven colliding foils. Phys. Rev. Lett. 70, 2102.Google Scholar
Chandrasekhar, S. (1961) Hydrodynamic and Hydromagnetic Stability. Oxford: Oxford University Press.Google Scholar
Clark, P.J., Eliezer, S., Farley, F.J.M., et al. (1985). Laser focus accelerator by relativistic self-focusing and high electric fields in double layers of nonlinear force produced cavitons. AIP Conf. Proc. 130, 380.Google Scholar
Eliezer, S. & Hora, H. (1989). Double layers in laser produced plasmas. Phys. Repts 172, 339.Google Scholar
Eliezer, S. (2002). The Interaction of High Power Lasers with Plasmas. Boca Raton: CRC press.Google Scholar
Eliezer, S., Eliaz, N., Grossman, E. et al. (2004). Synthesis of nanoparticles with femtosecond laser pulses. Phys. Rev. B 69, 144119.Google Scholar
Eliezer, S., Eliaz, N., Grossman, E., Fisher, D., Gouzman, I., Henis, Z., Pecker, S., Horovitz, Y., Fraenkel, M., Maman, S. & Lereah, Y. (2005). Nanoparticles and nanotubes induced by femtosecond lasers. Laser Part. Beams 23, 15.Google Scholar
Eliezer, S., Murakami, M. & Martinez Val, J.M. (2007). Equation of state and optimum compression in inertial fusion energy. Laser Part. Beams 25, 1.Google Scholar
Eliezer, S. & Martinez Val, J.M. (2011). The comeback of shock waves in inertial fusion energy. Laser Part. Beams 29, 175.CrossRefGoogle Scholar
Eliezer, S. (2012). Relativistic acceleration of micro-foils with prospects for fast ignition. Laser Part. Beams 30, 225.Google Scholar
Eliezer, S., Martinez Val, J.M. & Pinhasi, S.V. (2013). Relativistic shock waves in the laboratory. Laser Part. Beams 31, 113.Google Scholar
Esirkepov, T., Borghesi, M., Bulanov, S.V., Mourou, G. & Tajima, T. (2004). Highly efficient relativistic ion generation in the laser-piston regime. Phys. Rev. Lett. 92, 175003.Google Scholar
Hora, H., Lalousis, P. & Eliezer, S. (1984). Analysis of the inverted double layers produced by nonlinear forces in laser produced plasmas. Phys. Rev Lett. 53, 1650.Google Scholar
Hora, H., Min, G., Eliezer, S., Lalousis, P., Pease, R.S. & Szichman, H. (1989). On the surface tension of plasmas. IEEE Trans. Plasma Sci. 17, 283.CrossRefGoogle Scholar
Hora, H., Lalousis, P.& Moustaizis, S. (2014). Picosecond laser pulses for ultrahigh plasma acceleration to ignite fusion flames in solid density fuels. IFSA 2013 Proceedings, Nara, Japan.Google Scholar
Laluosis, P., Foldes, I.B. & Hora, H. (2012). Ultrahigh acceleration of plasma by picosecond terawatt laser pulses for fast ignition of fusion. Laser Part. Beams 30, 233.Google Scholar
Lalousis, P., Hora, H., Eliezer, S., Martinez Val, J.M., Moustaizis, S., Miley, G. H. & Mourou, G. (2013). Shock Mechanisms by ultrahigh laser accelerated plasma blocks in solid density targets for fusion. Phys. Lett. A 377, 885.Google Scholar
Mima, K., Murakami, M., Nakai, S. & Eliezer, S. (2009). Applications of Laser-Plasma Interactions. Boca Raton: CRC Press.Google Scholar
Murakami, M., Nagatomo, H., Azechi, H., Ogando, F. & Eliezer, S. (2006). Innovative ignition scheme for ICF impact fast ignition. Nucl. Fusion 46, 99.Google Scholar
Pegoraro, F. & Bulanov, S.V. (2007). Photon bubbles and ion acceleration in a plasma dominated by the radiation pressure of an electromagnetic pulse. Phys. Rev. Lett. 99, 065002.Google Scholar
Piazza, A.D., Muller, C., Hatsagortsyan, K. Z. & Keitel, C. H. (2012). Extremely high-intensity laser interactions with fundamental quantum systems. Rev. Mod. Phys. 84, 1177.CrossRefGoogle Scholar
Rayleigh, J.W.S. Lord. (1876). On waves. Phil. Mag. 1, 257279.Google Scholar
Tabak, M., Hammer, J., Glinsky, M.E., Kruer, W.L., Wilks, S.C., Woodworth, J., Campbell, E., Perry, M.D. & Mason, R.J. (1994). Ignition and high gain with ultra-powerful lasers. Phys. Plasmas 1, 1626.Google Scholar
Taylor, G. (1950). The instability of liquid surfaces when accelerated in a direction perpendicular to their planes. Proc. Roy. Soc. London A 201, 192.Google Scholar
Figure 0

Fig. 1. (Color online) The final velocity of the accelerated foil before breakdown by Rayleigh Taylor instability according to Eq. (11). In this equation ln(ξ/ξ0)cr was chosen to be equal to 3 and the foil density ρ =1 g/cm3. The ratio of foil thickness (l) to the transversal dimension (L) are chosen for (a) 2πl/L = 1, (b) 2πl/L = 0.1.

Figure 1

Fig. 2. (Color online) (a) The capacitor model for laser irradiances IL where the ponderomotive force dominates the interaction. Ex is the electric field inside the double layer, λDL is the distance between the positive and negative DL charges, l is the foil thickness and δ is the solid density skin depth of the foil. (b) A schematic figure that our capacitor model is based. ne and ni are the electron and ion densities accordingly, Ex is the electric field solution. The DL is geometrically followed by neutral plasma where the electric field decays within a skin depth.

Figure 2

Fig. 3. (Color online) Micro-foil velocity βfc as a function of the laser pulse duration in units of τ = ρ0c2l/(2IL), where ρ0 is the micro-foil initial velocity upon impact and l is the foil thickness.