the rapid development and improvement of ultra - intense laser pulses have opened new areas of physics for fundamental research and have also enabled novel technological applications .
relativistically intense ( @xmath0 w/@xmath1 ) laser pulses readily ionize matter converting it into a plasma .
the interaction with the plasma electrons is the primary channel for the energy transfer from the laser pulse , providing the basis for a wide range of phenomena and applications .
specifically , generation of copious relativistic electrons is the key to x - ray  @xcite and secondary particle sources , such as energetic ions  @xcite , neutrons  @xcite , and positrons  @xcite .
it is therefore critical to understand what controls the generation of relativistic electrons in relevant regimes of laser - plasma interactions .
it is well recognized that the regime in which a relativistically intense laser pulse irradiates a sub - critical plasma is optimal for generating relativistic electrons .
the pulse can propagate through such a plasma , which enables an extended interaction length with the electrons . in experiments aimed at x - ray generation ,
this regime is deliberately achieved by using an expanding gas jet  @xcite .
this regime can also naturally arise in experiments with solid density targets irradiated by a powerful laser pulse due to the presence of a prepulse  @xcite .
the prepulse often delivers a considerable amount of energy , causing the front of the target to expand and form an extended subcritical preplasma prior to the arrival of the main pulse .    the important role played by the interaction of the main pulse with the preplasma in generating
an energetic electron population in experiments with initially solid - density targets was however not immediately recognized .
this is in part due to the fact that in some setups the role of the prepulse can also be detrimental .
it is the case in experiments aimed at ion acceleration , where the pre - expansion at the surface of a thick bulk target caused by the prepulse significantly reduces the effectiveness of ion acceleration  @xcite .
recently , the focus of the ion acceleration research has markedly shifted towards those regimes in which the target becomes transparent to the main relativistically intense pulse and the electron acceleration and heating can be fully utilized  @xcite .
such a regime is achieved by employing an initially ultra - thin target that significantly expands during the pre - pulse and forms a relativistically transparent plasma extending many wavelengths along the direction of the laser beam propagation .
the same regime has also been used to generate energetic electrons that are subsequently converted into a short neutron beam  @xcite .
how energetic electrons are generated strongly depends on the duration of the main laser pulse . in most setups aimed at producing secondary particle sources , the laser pulse irradiates the plasma over a time period that is longer than the characteristic electron response time .
this regime is the opposite to the regime used for the wakefield acceleration  @xcite . as a consequence
, the laser pulse establishes a quasi - steady - state structure in the plasma that slowly evolves on an ion time scale  @xcite . by expelling some of the electrons radially ,
the laser creates a positively charged elongated channel in the plasma with quasi - static transverse and longitudinal electric fields .
new electrons are continuously injected into the channel , typically through the channel opening  @xcite , and then get accelerated and pushed forward by the laser pulse in the presence of the quasi - static fields .
the described regime is often broadly referred to as the direct laser acceleration regime . in terms of applications
, it is critical to know what controls the electron energy gain
. early work on the topic  @xcite indicated that the transverse static electric field can be beneficial for enhancing the electron energy gain beyond what is expected from a single electron irradiated by a plane wave in a vacuum .
recently , there has been a renewed interest in the direct laser acceleration of electrons , as experimental groups shift their focus to regimes of relativistic transparency and also try to optimize the preplasma conditions using multiple laser pulses  @xcite .
recent simulation results have also demonstrated that the direct laser acceleration can be important in the context of the laser wakefield acceleration  @xcite .    in this paper
, we examine the role played by transverse and longitudinal quasi - static electric fields present in a plasma channel in enhancing the electron energy gain from the laser pulse .
we also address the role of electron injection into the laser beam and the limitations imposed by the super - luminosity of the laser field that is induced by the channel .
this paper is based on a body of work performed by us over the last couple of years  @xcite and is designed to serve in part as an overview of some novel aspects of the direct laser acceleration . here
we focus on illustrating the key qualitative concepts and phenomena , while providing references to those publications where one can find more detailed technical analysis .
the nature of the laser - plasma interaction strongly depends on the amplitude of the irradiating laser pulse and on the electron density of the irradiated plasma .
it is convenient to use a dimensionless parameter @xmath2 to quantify the impact of a laser pulse with electric field amplitude @xmath3 and frequency @xmath4 on electron motion . here
@xmath5 is the speed of light and @xmath6 and @xmath7 are the electron charge and mass .
the parameter @xmath8 is often referred to as the normalized laser amplitude .
it is roughly the ratio of the transverse electron momentum induced by the oscillating laser electric field to @xmath9 .
therefore , a laser pulse with a normalized amplitude of @xmath10 would induce relativistic electron motion .
the electron density in the plasma determines whether the laser pulse can propagate into the plasma and accelerate plasma electrons .
the cut - off for a pulse with @xmath11 occurs at a critical density , @xmath12 for which the electron plasma frequency @xmath13 becomes equal to the frequency of the laser pulse . at laser amplitudes @xmath10
, the plasma can become relativistically transparent at electron densities exceeding the critical density @xmath14 .
the adjusted critical density in this case depends on the amplitude of the irradiating laser pulse , because the effect is caused by the relativistic motion of electrons in the strong field of the laser .
the optimal regime for generating copious relativistic electrons is then the regime in which a relativistic amplitude laser pulse @xmath15 irradiates an extended sub - critical plasma @xmath16 . in order to illustrate the key features of the laser - plasma interaction in this regime ,
we have performed a two - dimensional ( 2d ) particle - in - cell ( pic ) simulation whose results are shown in fig .
[ fig:2d_example ] . in this simulation , a uniform sub - critical plasma with initial electron density @xmath17
is irradiated by a laser beam with wavelength @xmath18 @xmath19 m whose amplitude ramps up to @xmath20 and then remains constant .
the laser pulse propagates along the @xmath21-axis and it is linearly polarized , with the laser electric field polarized in the @xmath22-plane .
the ions were kept immobile in this simulation to distinguish more clearly the effect of the long laser pulse .
detailed parameters of the simulation are given in appendix  [ appendix_1 ] .
as the pulse enters the plasma , its ponderomotive force begins to expel some of the electrons out of the laser pulse in the transverse direction producing a channel .
the un - neutralized ion charge generates a counteracting force that , in the example shown in fig .
[ fig:2d_example ] , prevents the channel from becoming fully evacuated .
the laser pulse produces and maintains a steady - state channel if the pulse duration exceeds the characteristic electron response time . the snapshot of the electron density in fig .
[ fig:2d_example ] taken at 1 ps illustrates such a channel .
the positively charged elongated channel generates quasi - static transverse and longitudinal electric fields shown in the lower - left panel of fig .
[ fig:2d_example ] .
the fields have been averaged over ten laser periods and they are normalized to the electric field amplitude @xmath3 of a plane wave with @xmath20 , which is essentially the amplitude of the electric field in the laser pulse .
the time - averaged electric fields are relatively small compared to the amplitude of the oscillating laser electric field that is also present in the channel .
once the steady - state channel structure is established , new electrons are continuously injected into the channel through the opening .
this is particularly clear from the snapshot of the time - averaged current density , whose absolute value is shown in the upper - right panel of fig .
[ fig:2d_example ] .
the injected electrons are accelerated and pushed forward by the laser pulse , producing a steady - state electron current in the channel .
this current generates a quasi - static transverse magnetic field directed in and out of the plane of the simulation ( along the @xmath23-axis ) .
the return current flowing outside of the channel causes for the magnetic field to be localized inside the channel , as shown in the lower - right panel of fig .
[ fig:2d_example ] .
the plotted time - averaged field is normalized to the magnetic field amplitude @xmath24 of a plane wave with @xmath20 .
the quasi - static magnetic field is also relatively weak in this regime , as it is less that 1% of the magnetic field in the laser pulse .    .
color - coded is the electron longitudinal momentum . ]
the presented example illustrates that , in a sub - critical plasma irradiated by a long laser pulse , electron acceleration takes place in a positively charged channel in the presence of extended transverse quasi - static electric and magnetic fields and a localized quasi - static longitudinal electric field .
these fields are relatively small compared to the fields in the laser pulse , not exceeding a few percent .
it is then somewhat unexpected that , for example , the transverse electric field can significantly enhance the electron energy gain  @xcite .
as we show in the following sections , the transverse and longitudinal electric fields can synergistically enhance the electron energy gain , with a well - pronounced threshold , well beyond the energy gain that one would expect in the absence of these fields .
in order to provide the context for the discussion of the electron acceleration in the channel , we briefly review the key features of the electron motion in a vacuum where no static electric or magnetic fields are present . specifically
, we consider a single electron that is initially at rest .
it is irradiated by a plane electromagnetic wave whose normalized amplitude gradually increases from zero to @xmath8 . in what follows
, we refer to this regime as the vacuum regime .    the electron moves according to the following equations : @xmath25 , \label{main_eq:1 } \\ & & \frac{d { \bf{r}}}{d t } = \frac{c}{\gamma } \frac{{\bf{p}}}{m_e c } , \label{main_eq:2}\end{aligned}\ ] ] where @xmath26 is the relativistic factor , @xmath27 and @xmath28 are the electron position and momentum , and @xmath29 is the time . in the regime under consideration ,
@xmath30 and @xmath31 are the electric and magnetic fields of the wave .
it is convenient to express these fields in terms of a normalized vector potential @xmath32 , @xmath33 .
\label{b_wave}\end{aligned}\ ] ] in the case of a plane wave with wave - length @xmath34 propagating along the @xmath21-axis , the vector potential is only a function of a normalized phase @xmath35 without any loss of generality , we assume that the laser electric field is polarized along the @xmath36-axis . in this case , the vector potential @xmath32 only has a @xmath36-component .
we consider a pulse with @xmath37 , where @xmath8 is the maximum normalized amplitude and @xmath38 is a slowly varying envelope that ramps up from zero to unity .
figure [ figure_dla_vac_1 ] shows an electron trajectory for a pulse with @xmath20 obtained by solving eqs .
( [ main_eq:1 ] ) and ( [ main_eq:2 ] ) numerically .
it illustrates the qualitative change that takes place at relativistic wave amplitudes of @xmath10 . at @xmath11 , the dominant force experienced by
the electron is the force from the electric field of the wave . as a result
, the electron oscillates across the laser pulse , while its longitudinal displacement is negligible . at @xmath10 ,
the momentum oscillations induced by the laser electric field become relativistic and , as a consequence of this , the lorentz force becomes important .
this force causes longitudinal electron motion that leads to a trajectory shown in fig .
[ figure_dla_vac_1 ] .    even though the energy from the laser pulse is transferred to the transverse electron oscillations , the lorentz force converts most of this energy into the longitudinal electron motion . indeed , an analytical solution of eqs .
( [ main_eq:1 ] ) and ( [ main_eq:2 ] ) yields  @xcite @xmath39 the electron moves along a parabola in the momentum space , with @xmath40 for @xmath41 .
the maximum @xmath42-factor that the electron can reach is thus given by @xmath43 the change of the longitudinal electron momentum along the electron trajectory is shown in fig .
[ figure_dla_vac_1 ] .
-factor ( lower panel ) as functions of the phase variable @xmath44 in a pulse with @xmath20 .
the color coding is the amplitude of the laser electric field .
the arrows show the direction of the transverse electron velocity along the trajectory . ]
the energy gain is limited by the dephasing between the electron and the wave .
this aspect is illustrated in fig .
[ figure_dla_vac_2 ] , where the electron trajectory and the @xmath42-factor are shown as functions of the phase variable @xmath44 . in the upper panel , the color - coding shows the amplitude of the laser electric field normalized to its maximum amplitude @xmath3 .
the phase of the wave at the electron location , @xmath44 , continuously decreases , because the electron longitudinal velocity is less than the phase velocity @xmath5 .
the electron , shown with a circle , is slipping with respect to the laser wave fronts , thus moving to the left in the upper panel . at @xmath45 , the electron transverse velocity ( shown with arrows )
is anti - parallel to the laser electric field , so the electron is gaining energy from the laser pulse . as the electron continues to slip , the laser electric field eventually changes its sign at @xmath46 and becomes positive . at this point ,
the electron transverse velocity is still positive and the electron starts to lose its energy .
this positive field continues to reduce the transverse electron velocity until it becomes negative at @xmath47 , at which point the electron begins to gain energy again .
the electron experiences energy gain twice every laser cycle , as evident from the plot of the @xmath42-factor in fig .
[ figure_dla_vac_2 ] .
it is worth pointing out that the electron comes momentarily to a complete stop also twice every laser cycle .
the longitudinal distance that the electron travels between the stops to achieve @xmath48 is much longer than the laser wavelength due to the fact that the electron is moving forward and it can be estimated as @xmath49 .
the discussed mechanism of electron acceleration and energy gain in the vacuum regime suggests that , in principle , there are two alternatives for further enhancing the electron energy gain .
one option is to decrease the dephasing between the electron and the wave , allowing the transverse electron velocity to remain antiparallel to the laser electric field for longer in terms of the actual time .
the second option is to change the oscillations of the transverse velocity , so that there is a net energy gain by the electron after a laser cycle . in what follows ,
we show how these scenarios can be realized utilizing relatively weak static longitudinal and transverse electric fields .
as discussed in sec .  [ sec - channel ] , one of the key features of the steady - state channel is the continuous injection of new electrons through the channel opening with a quasi - static longitudinal electric field .
this aspect raises the question of the impact of the quasi - static longitudinal electric field on subsequent electron energy gain during acceleration by the laser pulse along the channel .    in order to pinpoint the effect of the longitudinal electric field , we begin by considering electron motion in a vacuum using the same setup as in sec .  [ sec - vac ] , but with an added narrow region @xmath50 of a uniform longitudinal static electric field @xmath51 .
figure [ figure_dla_ez_1 ] shows numerical solutions of eqs .
( [ main_eq:1 ] ) and ( [ main_eq:2 ] ) for @xmath52 and @xmath53 , where @xmath3 is the electric field amplitude of the laser pulse .
the laser pulse amplitude is again @xmath20 following a gradual initial ramp - up .
the difference between the two examples shown in fig .
[ figure_dla_ez_1 ] is the location of the field region along the electron trajectory . in the first example ( left panel )
, the interaction takes place when the electron momentum is close to its maximum value . the subsequent change in the maximum electron momentum after the interaction is insignificant .
this result is in agreement with the expectation that the localized field we are considering is too weak to transfer significant energy to the electron . in the second example ( right panel ) , the interaction takes place as the electron comes to a stop . in stark contrast to the previous example ,
the maximum electron energy more than doubles during the motion in the wave after the interaction .
that is @xmath53 long .
the only difference between the two cases is the location of the field region . ]
the key role of the longitudinal electric field in these examples is not to directly increase the electron energy , but rather to decrease the electron dephasing with the wave .
we define the dephasing , @xmath54 as the rate at which the wave phase at the electron location , @xmath44 , changes with proper time , @xmath55 , defined as @xmath56 it follows directly from the definition of @xmath44 given by eq .
( [ xi ] ) that @xmath57 the quantity @xmath58 is conserved by eqs .
( [ main_eq:1 ] ) and ( [ main_eq:2 ] ) in the absence of the static electric field , which means that the dephasing defined by eq .
( [ r ] ) is constant before and after the interaction with the longitudinal field .
it can also be shown that during the interaction the static field reduces the dephasing at the rate  @xcite @xmath59 that is independent of the laser pulse amplitude .
the reduced dephasing allows the electron to spend more time being accelerated by the wave and thus gain more energy from the laser pulse . as shown in fig .
[ figure_dla_vac_2 ] of sec .
[ sec - vac ] , the acceleration continues before the phase decreases by @xmath60 .
a reduction in the dephasing then means that the corresponding proper time interval @xmath61 of electron acceleration is indeed longer [ see eq .
( [ r ] ) ] .
the maximum @xmath42-factor that the electron with reduced dephasing achieves  @xcite is @xmath62     fs from a 2d pic simulation with @xmath20 and @xmath63 . ]
not all parts of the electron trajectory are equivalent in terms of the dephasing change induced by the longitudinal field .
the reason for that is the strong dependence of the dephasing on the longitudinal momentum . for simplicity ,
let us assume that that the electric field instantaneously increases the longitudinal momentum by @xmath64 .
prior to the interaction , we have @xmath65 here we used expressions ( [ vac_px ] ) and ( [ vac_py ] ) for the electron momentum in the vacuum regime . assuming that @xmath64 is smaller than @xmath66 and @xmath67 , we find from eq .
( [ r2 ] ) that the dephasing is reduced to @xmath68 evidently , the most favorable parts of the electron trajectory for reducing the dephasing are those parts where the longitudinal momentum is the lowest .
this aspect is at work in the two examples shown in fig .
[ figure_dla_ez_1 ] , where the dephasing drops to @xmath69 on the right and only to @xmath70 on the left .
the resulting energy increase in both cases is in good agreement with eq .
( [ gamma_long ] ) .    , on top of the time - averaged electron density profile .
the lower panel shows the trajectory , with the color - coded dephasing rate @xmath71 , on top of the time - averaged longitudinal electric field . the field and the density
are averaged over ten laser periods at 850 fs . ]    in order to examine how this mechanism is realized in a plasma channel , we have performed a 2d pic simulation in which a long laser pulse with an electric field polarized perpendicular to the plane of the simulation irradiates a plasma slab with @xmath63 .
detailed parameters of the simulation are given in appendix  [ appendix_1 ] .
the normalized laser amplitude is @xmath20 , so the @xmath42-factor in the vacuum regime would be limited by @xmath72 .
however , the electron spectrum in the plasma , as seen in fig .
[ figure_spectrum ] , has a pronounced energetic electron tail with @xmath73 .    a trajectory of one of the energetic electrons is shown in fig .  [ figure_dla_ez_2 ] together with a time - averaged electron density profile and longitudinal electric field .
the electron is injected into the channel near the opening and initially moves against the laser pulse .
the electron typical @xmath42-factor at this stage is comparable to @xmath8 .
the quasi - static negative electric field , that is present in the channel opening , gradually turns the electron around .
as the electron starts to move in the direction of the laser pulse propagation , the regime becomes similar to that considered earlier in this section .
the longitudinal electric field causes the dephasing rate to decrease , as shown in the lower panel of fig .
[ figure_dla_ez_2 ] . similarly to the simplified example , the significant drop in the dephasing occurs as the @xmath42-factor reaches its minimum at @xmath74 fs .
the energy increase takes place as the electron continues its longitudinal motion into the channel after leaving the region with a negative quasi - static electric field .
the reduced dephasing is evident not only from the color - coding in the lower panel of fig .
[ figure_dla_ez_2 ] , but also from the lack of oscillations in @xmath75 in fig .
[ figure_dla_ez_3 ] .
the electron @xmath42-factor peaks at @xmath76 fs , roughly 20 @xmath19 m after the energy increase began .
it exceeds @xmath77 by a factor of 2.3 .
most of the energy comes directly from the laser electric field and not from the longitudinal field , as shown in the upper panel of fig .
[ figure_dla_ez_3 ] .
the laser contribution to the @xmath42-factor is calculated by integrating the work done by @xmath78 and normalizing it to @xmath79 .
-factor and momentum components as functions of time for the electron trajectory shown in fig .
[ figure_dla_ez_2 ] .
the upper panel also shows the integrated contribution to the @xmath42-factor from the work done by the laser electric field . ]
we have so far assumed that the phase velocity in the laser pulse , @xmath80 , is equal to the speed of light , @xmath5 . however , the plasma in the channel and the channel itself cause dispersion , making the phase velocity superluminal , @xmath81 .
the effect of the superluminosity is faster dephasing between the electron and the wave and , as a result , lower energy gain from the wave .
the superluminosity must be taken into account if it leads to a significant change in the dephasing .
in general , the dephasing is determined by by the difference @xmath82 .
we can then conclude that the superluminosity becomes important for @xmath83 , whereas we can set @xmath84 for @xmath85 . using the definitions @xmath86 and @xmath87 , we find that @xmath88 in the case when @xmath84 , the @xmath42-factor is given by eq .
( [ gamma_long ] ) . taking into account this estimate and the expression for @xmath77
, we find that the role of the superluminosity is negligible for  @xcite @xmath89    ultimately , the superluminosity sets a limit on the energy gain from the wave .
it follows from eq .
( [ v_ph_cond_2 ] ) that there is a critical value for @xmath86 , given by @xmath90 for @xmath91 , the dephasing between the electron and the wave is determined primarily by the difference @xmath92 and
the difference between the electron longitudinal velocity and the speed of light is unimportant .
this means that even if the longitudinal electric field lowers @xmath93 well below @xmath94 , the maximum @xmath42-factor that the electron can achieve will be approximately limited by @xmath95 which was estimated from eq .
( [ gamma_long ] ) by setting @xmath96 .     is the location of the center of the window .
the field is normalized to the electric field amplitude @xmath3 of a plane wave with @xmath20 .
the dashed lines represent wavefronts moving with two different superluminal phase velocities @xmath80 . ]    in order to examine the conditions ( [ v_ph_cond_2 ] ) and ( [ v_ph_cond_3 ] ) for the simulation shown in fig
.  [ figure_dla_ez_2 ] , we have plotted the wavefronts in a window moving with the speed of light along the axis of the channel .
figure  [ figure_v_ph ] shows the wave electric field @xmath78 on the axis of the channel in a window that is @xmath97 wide .
the center of the window , located at @xmath98 , slides the entire length of the channel shown in fig .
[ figure_dla_ez_2 ] ( note that @xmath18 @xmath19 m in the simulation ) .
the two dashed lines indicate wavefront slopes for @xmath99 and @xmath100 . inside the channel we then approximately have @xmath99 .
the corresponding critical dephasing is @xmath101 and the maximum @xmath42 that can be attained by accelerating electrons is @xmath102 , which translates into @xmath103 for @xmath20 .    in the light of these estimates ,
a closer look at the electron data shown in figs .
[ figure_dla_ez_2 ] and [ figure_dla_ez_3 ] reveals that the superluminosity likely plays a role in limiting the electron energy gain .
the dephasing rate @xmath93 shown in fig .  [ figure_dla_ez_2 ] falls to @xmath104 along the electron trajectory , which could lead to a factor of ten increase of the @xmath42-factor compared to @xmath77 for @xmath84 .
however , our estimates show that reduction of the dephasing below @xmath101 no longer leads to an enhancement of the @xmath42-factor , capped at @xmath103 , because of the superluminosity .
this is in good agreement with the @xmath42-profile shown in the upper panel of fig .
[ figure_dla_ez_3 ] .
finally , it is worth pointing out that heating of the bulk electrons in the channel to relativistic energies significantly reduces the superluminosity .
the phase velocity of a linear electromagnetic wave propagating through a cold non - relativistic plasma with electron density @xmath105 is @xmath106 .
we then should expect @xmath107 for @xmath63 used in the simulation .
the corresponding wavefront should follow the right dashed line in fig .
[ figure_v_ph ] . the actual @xmath92 is smaller by an order of magnitude , which clearly indicates that the critical density that determines the phase velocity is effectively lowered by an order of magnitude due to relativistic motion induced in the channel by the laser pulse .
this so - called relativistic transparency enables the electron energy gain enhancement by reducing the superluminosity and thus reducing the critical dephasing given by eq .
( [ critical_dephasing ] ) .
inside the channel , the electron acceleration by the laser pulse takes place in the presence of a relatively weak quasi - static transverse electric field . in order to examine the role played by such a field ,
we consider just a single electron irradiated by a plane electromagnetic wave in a fully evacuated cylindrical ion channel shown in fig .
[ figure_cyl_setup ] .    -plane . ]    the electron dynamics in this model is described by eqs .
( [ main_eq:1 ] ) and ( [ main_eq:2 ] ) , where the electric and magnetic fields are given .
these fields are a superposition of the fields of the wave given by eqs .
( [ e_wave ] ) and ( [ b_wave ] ) and the electric field of the channel , @xmath108 where @xmath109 . we have assumed for simplicity that the channel is an infinitely long uniform cylinder that consists of immobile singly - charged ions with density @xmath110 .    it can be shown that , as the electron moves along the channel , the following quantity remains conserved @xmath111 this relation indicates that amplification of electron oscillations across the channel leads to a reduction of @xmath58 .
since @xmath66 is the dominant component of the electron momentum , a significant reduction of @xmath58 implies a significant increase of @xmath66 .
therefore , the relation ( [ q_def ] ) formally points at a direct connection between enhancement of transverse oscillations and a significant increase of the longitudinal electron momentum .
the integral of motion ( [ q_def ] ) sets an upper limit for the amplitude of the transverse oscillations across the channel .
the value of @xmath112 is determined by the initial conditions , so that @xmath113 for an on - axis electron that is initially at rest .
we find from eq .
( [ q_def ] ) that the amplitude of the electron oscillations can not exceed @xmath114 for @xmath113 .
the amplitude of the oscillations approaches @xmath115 as @xmath116 and @xmath117 .
it is important to distinguish transverse electron oscillations in the @xmath36 and @xmath23 directions , because the motion in the @xmath36-direction is driven by the laser electric field , whereas the motion in the @xmath23-direction is only affected by the electric field of the channel .
we first focus on the key aspects of the driven electron oscillations .
we consider an electron that has no initial displacement or momentum in the @xmath23-direction , so that the electron trajectory driven by the laser pulse remains flat , with the electron moving only in the @xmath22-plane .
the equation for the electron motion across the channel can be written in the form that resembles an equation for a driven oscillator  @xcite : @xmath118 where @xmath119 and @xmath55 is the proper time defined by eq .
( [ tau_def ] ) .
the wave amplitude is a function of the phase variable @xmath44 that changes at the rate @xmath120    there are two characteristic frequencies in this case : the frequency of the oscillations induced by the field of the channel electric field and the frequency of the oscillations induced by the field of the laser pulse .
the former is the natural frequency @xmath121 .
the latter is equal to @xmath4 in the vacuum regime , because @xmath122 and , as a result , @xmath123 . according to eq .
( [ q_def ] ) , the relation @xmath124 holds in the presence of the transverse electric field as well while the amplitude of the transverse oscillations is much smaller than @xmath125 .    at low channel densities , such that @xmath126 , the channel electric field has very little impact on the electron motion .
the frequency mismatch means that a resonant interaction that can lead to an amplification of the transverse oscillations is not possible .
the electron motion is essentially the same as in the vacuum regime , with @xmath127 .
the characteristic amplitude of the transverse oscillations , estimated as @xmath128 from eq .
( [ perp_osc ] ) , is indeed much smaller than @xmath115 for @xmath126 and @xmath127 .
-factor ( lower panel ) of an electron irradiated by a plane wave with @xmath129 in a channel whose ion density corresponds to @xmath130 .
the color coding is the amplitude of the laser electric field .
the arrows show the direction of the transverse electron velocity along the trajectory . ]    at higher channel densities for which the two characteristic frequencies become comparable , @xmath131 , the channel electric field becomes capable of significantly changing the phase of the transverse oscillations , which causes their amplitude to grow . the change of phase after significant amplification of the transverse oscillations
has already taken place is illustrated in fig .
[ figure_dla_vac_3 ] , where the electron trajectory and the @xmath42-factor are shown as functions of the phase variable @xmath44 . in this example , the electron is initially at rest on the axis of the channel ( @xmath132 ) whose ion density corresponds to @xmath130 .
the electron is irradiated by a plane wave with @xmath133 and @xmath129 .
the upper panel of fig .
[ figure_dla_vac_3 ] shows that the channel electric field changes the oscillations of the transverse electron velocity , allowing it to remain anti - parallel to the laser electric field over extended segments of the electron trajectory marked with white .
as a result , there is a net energy gain each laser cycle for @xmath134 that leads to a significant increase of the relativistic @xmath42-factor well above what can be achieved in the vacuum regime illustrated in fig .
[ figure_dla_vac_2 ] .
note that @xmath135 along this segment of the electron trajectory .
the example in fig .
[ figure_dla_vac_3 ] confirms that a significant increase of the @xmath42-factor goes hand in hand with the enhancement of the transverse oscillations .
+    in order to determine the conditions for the energy enhancement , we have performed a parameter scan solving eqs .
( [ main_eq:1 ] ) and ( [ main_eq:2 ] ) numerically for a wide range of parameters @xmath8 and @xmath136 that specify the fields acting on the electron . in all the cases ,
the electron is initially at rest on the axis of the channel and the laser pulse amplitude gradually ramps up as @xmath137 \sin(\xi ) ,   & \text{for $ \xi >   \xi_0$;}\\        a_0 \sin(\xi ) ,   & \text{for $ \xi \leq   \xi_0 $ }     \end{cases } \label{a}\ ] ] where @xmath138 and @xmath139 .
the maximum @xmath42-factor , @xmath140 , that the electron attains for different sets of parameters @xmath8 and @xmath136 is shown in fig .
[ invited_pop_figure_8 ] , where @xmath140 is normalized to the maximum @xmath42-factor in the vacuum regime , @xmath77 , for the same laser amplitude @xmath8 .
the energy enhancement has a sharp threshold , as further illustrated by the lineouts in fig .
[ invited_pop_figure_8 ] for three different values of @xmath8 .
the lineouts for @xmath141 and @xmath142 indicate that the threshold and the energy gain above the threshold are determined by a single dimensionless combination @xmath143 the contours of constant @xmath144 plotted in the upper panel of fig .
[ invited_pop_figure_8 ] with dashed lines confirm that this is indeed a general trend for @xmath145 .
the existence of the threshold that depends only on @xmath146 was first discovered in ref .  [ ] and then explored in more detail in ref .
the dimensionless parameter @xmath144 is the ratio of the frequencies of oscillations induced by the field of the channel and by the field of the laser prior to a significant energy enhancement . indeed , taking into account that the electron @xmath42-factor in this regime is essentially @xmath77 , we find that @xmath147 for @xmath41 .
currently , the dependence of @xmath140 on @xmath144 above the threshold remains a robust numerical observation that requires an in - depth theoretical analysis .
-factor attained by injected electrons .
the laser amplitude in all the cases is @xmath20 .
the initial injection phase defines the initial wave amplitude , @xmath148 . ]      in the scan presented in fig .
[ invited_pop_figure_8 ] we have considered the case where initially the electron is already in the channel , but the laser pulse has not yet reached the electron location longitudinally .
on the other hand , in the steady - state channel discussed in sec .  [ sec - channel ] and illustrated in fig .
[ fig:2d_example ] , electrons are injected into the channel from the side with the laser beam already present in the channel .
this raises the question regarding the role of electron injection in determining both the threshold and the energy gain .
we test the role of injection by seeding an electron without any initial momentum into a plane wave with a given initial phase . specifically , the wave amplitude is set to @xmath149 , where @xmath150 is the phase of the laser pulse at the moment of injection defined as @xmath151 .
the electron is placed onto the channel axis at the moment of injection .
figure  [ figure_injection ] shows how the maximum @xmath42-factor attained by the electron depends on the injection phase @xmath150 and the ion density in the channel .
we show only the range of initial phases @xmath152 , because the plot for @xmath153 is identical .
the energy enhancement threshold is present for all the injection phases in fig .
[ figure_injection ] , but the location of the threshold is sensitive to the injection phase .
the threshold is the lowest when initially the electric field of the wave is equal to zero , @xmath154 .
since the wave amplitude is fixed , this means that the injection phase determines the value of the parameter @xmath144 at the threshold .
it should also be pointed that the energy gain increases as the threshold becomes lower .
therefore , we can conclude that proper electron injection into the laser beam can significantly enhance the electron energy gain from the wave , while lowering the corresponding threshold determined by the product @xmath146 .
-plane ( shown in light - blue ) . ]
we have so far deliberately limited the electron motion to the plane of the driven oscillations by placing the electron initially on the axis of the channel . small off - axis displacements in the absence of the laser pulse lead to small free oscillations across the channel since the restoring force is no longer equal zero . in what follows ,
we examine how the out - of - plane oscillations evolve in the presence of a laser pulse with relativistic amplitude @xmath41 .    the equation for the the out - of - plane oscillations can be directly obtain from eq .
( [ perp_osc ] ) by taking into account that the channel is cylindrically symmetric . by setting the right - hand side to zero and replacing @xmath36 with @xmath23 we obtain an equation , @xmath155 that resembles an equation for an oscillator with a natural frequency @xmath121 .
the free oscillations along the @xmath23-axis are coupled to the driven electron motion in the @xmath22-plane via the relativistic @xmath42-factor , @xmath156 .
as long as the amplitude of the free oscillations remains small , the @xmath42-factor is determined predominantly by the driven motion .
taking the @xmath42-factor in the vacuum regime shown in fig .
[ figure_dla_vac_2 ] as an example of the driven motion , we can conclude that the driven motion affects the natural frequency @xmath121 in two ways : it increases the natural frequency roughly by a factor of @xmath8 for @xmath41 and it also strongly modulates the natural frequency . the first aspect has been discussed in detail earlier in this section and it is responsible for the onset of the significant energy enhancement .
the modulation of the natural frequency can make the free oscillations along the @xmath23-axis parametrically unstable .
the @xmath42-factor modulations effectively modulate the restoring force acting on the oscillator described by eq .
( [ free_perp_osc ] ) .
it is then intuitively clear that the oscillations remain stable as long as the modulation frequency is much higher than the natural frequency @xmath121 .
the oscillations become unstable when the two frequencies become comparable .
as already discussed , the driven electron motion below the energy enhancement threshold can be approximated by the solution for the vacuum regime .
we again consider an electron that is irradiated by a plane wave in the channel and assume that the initial displacement is much smaller than @xmath115 . in this case ,
@xmath157 and we then find from the @xmath42 profile shown in fig .
[ figure_dla_vac_2 ] that the modulation frequency is @xmath158 . on the other hand , significantly below the threshold we have @xmath126 , as evident from fig .  [ invited_pop_figure_8 ] . therefore ,
free oscillations should be stable below the threshold .
as we cross the threshold by increasing the ion density , two important changes in the dependence of the @xmath42-factor on @xmath55 take place .
the @xmath42-factor is significantly increased , which greatly increases the natural frequency .
in addition to that , the modulations of the @xmath42-factor become less frequent .
this is because the peaks of enhanced @xmath42 act as modulations and they take multiple laser cycles to develop .
the combination of these two factors changes the relation between the two characteristic frequencies to @xmath159 .
we therefore conclude that the energy enhancement changes the stability of the free oscillations , making them susceptible to the parametric instability .    a parameter scan over the same range of @xmath8 and @xmath136 as in fig .  [ invited_pop_figure_8 ]
reveals that there is a threshold for the amplification of the free transverse oscillations  @xcite .
this threshold matches the energy enhancement threshold ( dashed line in fig .
[ invited_pop_figure_8 ] ) , because the energy enhancement triggers the growth of the electron oscillations along the @xmath23-axis .
the amplitude of the free oscillations grows until it becomes comparable to @xmath115 .
the amplitude of the driven oscillations , i.e. the oscillations along the @xmath36-axis , also has an enhancement threshold that matches the energy enhancement threshold as well .
we find that , in general , the electron energy enhancement is accompanied by a considerable enhancement of the transverse oscillations across the channel , making the electron trajectory considerably three - dimensional with the maximum displacement comparable to @xmath125 .
to illustrate this point , we show in fig .
[ figure_9 ] electron trajectories for @xmath20 and @xmath160 .
these parameters are above the threshold according to fig .
[ invited_pop_figure_8 ] . in the upper panel ,
the electron has no initial out - of - plane displacement , so its trajectory remains flat .
however , an initial displacement of @xmath161 in the lower panel provides a necessary seed for the instability to develop significant out - of - plane oscillations . as a result
, a non - planar electron trajectory quickly develops , accompanying the energy enhancement .
we have examined a regime in which a laser pulse with relativistic intensity irradiates a sub - critical plasma over a time period much longer than the characteristic electron response time .
it is shown that a steady - state channel is formed in the plasma in this case with quasi - static transverse and longitudinal electric fields .
these fields , even though they are much smaller than the electric field in the laser pulse , profoundly alter the electron dynamics .
the longitudinal electric field reduces the longitudinal dephasing between the electron and the wave .
this allows the electron to gain significantly more energy from the wave . the energy gain in this regime
is ultimately limited by the superluminosity of the wave fronts induced by the plasma in the channel , even though the relativistic transparency greatly reduces the role of the plasma .
the transverse electric field alters the oscillations of the transverse electron velocity . as a result
, it can remain anti - parallel to the laser electric field , which leads to significant energy gain by the electron from the wave .
we showed that this process has a sharp threshold determined by a single parameter , @xmath146 .
the threshold can be greatly reduced and the energy gain further enhanced by appropriately injecting electrons into the laser pulse . the threshold for the energy enhancement matches the threshold for the onset of the parametric instability in the direction perpendicular to the plane of the driven electron oscillation .
therefore , the trajectories of electrons with enhanced energies eventually become non - planar if the interaction time and length with the pulse are sufficient for the instability to develop .
the effect of the longitudinal electric field is akin to pre - acceleration of electrons in the longitudinal direction .
analysis of pre - accelerated electrons in a channel  @xcite has showed that pre - acceleration reduces the energy enhancement threshold and increases the electron energy gain .
this indicates that the transverse and longitudinal electric fields can synergistically enhance the electron energy gain .
further in - depth research is required to determine the energy scalings and the key parameters that determine the electron dynamics .
finally , it should be pointed out that even though the presented analysis was motivated by laser - plasma interactions at significantly sub - critical plasma densities , most of the results are much more general and they are not limited just to the regime where @xmath162 .
therefore , the analysis is well suited to make meaningful predictions regarding electron acceleration in near - critical and over - critical plasmas , provided that such plasmas are relativistically transparent to the incoming high - intensity laser pulse  @xcite .
however , wave propagation becomes a crucial aspect in this case and it should be addressed self - consistently  @xcite taking into account transverse laser pulse dimensions and its polarization  @xcite .
ava would like to thank dr .
toma toncian for constructive comments .
this material is based upon work supported by the u.s .
department of energy [ national nuclear security administration ] under award number de - na0002723 .
was also supported by afosr contract no .
fa9550 - 14 - 1 - 0045 , u.s .
department of energy - national nuclear security administration cooperative agreement no .
de - na0002008 , and u.s .
department of energy contract no .
de - fg02- 04er54742 .
v.n.k . and g.s .
were supported by afosr contract no .
fa9550 - 14 - 1 - 0045 and u.s .
department of energy contract nos .
de - sc0007889 and de - sc0010622 .
simulations for this paper were performed using the epoch code ( developed under uk epsrc grants ep / g054940/1 , ep / g055165/1 and ep / g056803/1 ) using hpc resources provided by the texas advanced computing center at the university of texas .
the snapshots shown in figs .
[ fig:2d_example ] and [ figure_dla_ez_2 ] are from two - dimensional particle - in - cell simulations performed using an open - source code epoch  @xcite . in both simulations ,
the initial electron density profile is given by @xmath163 ,   & \text{for $ x <   x_0$;}\\        n_0 ,   & \text{for $ x \geq   x_0 $ } ,     \end{cases } \label{n_sim}\ ] ] where @xmath164 @xmath19 m and @xmath165 @xmath19 m . in fig .
[ fig:2d_example ] , we set @xmath166 , whereas , in fig .
[ figure_dla_ez_2 ] , we set @xmath167
. the electron population is initialized using 10 macro - particles per cell .
the ion density is initially equal to the electron density and the ion density profile is initialized using 4 macro - particles per cell .
the electrons are initially cold , whereas the ions are immobile throughout the simulation .
the size of the domain in the simulation shown in fig .
[ fig:2d_example ] is 180 @xmath19 m ( 12000 cells ) along the @xmath21-axis and 60 @xmath19 m ( 1200 cells ) along the @xmath36-axis .
the size of the domain in the simulation shown in fig .
[ figure_dla_ez_2 ] is 140 @xmath19 m ( 9500 cells ) along the @xmath21-axis and 60 @xmath19 m ( 1200 cells ) along the @xmath36-axis .
the longitudinal resolution is chosen in agreement with the criterion outlined in ref .
[ ] in order to correctly compute the electron energy gain .
we use a laser pulse whose focal plane in the absence of plasma is located at @xmath168 @xmath19 m .
the laser fields have a gaussian profile along the @xmath36-axis , with the electric and magnetic fields in the focal plane given by @xmath169 where @xmath170 @xmath19 m and the temporal profile is @xmath171 ,   & \text{for $ t <   t_0$;}\\        1 ,   & \text{for $ t \geq   t_0 $ } ,     \end{cases } \label{e_sim}\ ] ] with @xmath172 fs and @xmath173 fs . in both simulations , the laser wavelength is 1 @xmath19 m and the peak laser intensity at @xmath174 @xmath19 m is @xmath175 w/@xmath1 , which corresponds to @xmath20 .
the normalizing amplitudes @xmath3 and @xmath24 are maximum electric and magnetic fields in a plane wave with @xmath20 . in the simulation shown in fig .  [ fig:2d_example ] , the laser electric field has @xmath21 and @xmath36 components , whereas the laser magnetic field has only a @xmath23-component . in the simulation shown in fig .  [ figure_dla_ez_2 ]
, the laser electric field has only a @xmath23-component , whereas the laser magnetic field has @xmath21 and @xmath36 components .
10    h .- s .
park , d. m. chambers , h .- k .
chung , r. j. clarke , r. eagleton , e. giraldez , t. goldsack , r. heathcote , n. izumi , m. h. key , j. a. king , j. a. koch , o. l. landen , a. nikroo , p. k. patel , d. f. price , b. a. remington , h. f. robey , r. a. snavely , d. a. steinman , r. b. stephens , c. stoeckl , m. storm , m. tabak , w. theobald , r. p. j. town , j. e. wickersham , and b. b. zhang , phys .
plasmas * 13 * , 056309 ( 2006 ) .
s. kneip , s. r. nagel , c. bellei , n. bourgeois , a. e. dangor , a. gopal , r. heathcote , s. p. d. mangles , j. r. marques , a. maksimchuk , p. m. nilson , k. ta phuoc , s. reed , m. tzoufras , f. s. tsung , l. willingale , w. b. mori , a. rousse , k. krushelnick , and z. najmudin , phys .
* 100 * , 105006 ( 2008 ) .
i. pomerantz , e. mccary , a. meadows , a. arefiev , a. bernstein , c. chester , j. cortez , m. donovan , g. dyer , e. gaul , d. hamilton , d. kuk , a. lestrade , c. wang , t. ditmire , and b. hegelich , phys .
* 113 * , 184801 ( 2014 ) .
m. schollmeier , a. b. sefkow , m. geissel , a. v. arefiev , k. a. flippo , s. a. gaillard , r. p. johnson , m. w. kimmel , d. t. offermann , p. k. rambo , j. schwarz , and t. shimada , phys .
plasmas * 22 * , 043116 ( 2015 ) .
h. chen , a. link , y. sentoku , p. audebert , f. fiuza , a. hazi , r. f. heeter , m. hill , l. hobbs , a. j. kemp , g. e. kemp , s. kerr , d. d. meyerhofer , j. myatt , s. r. nagel , j. park , r. tommasini , and g. j. williams , phys .
plasmas * 22 * , 056705 ( 2015 )
.      j. fuchs , c. a. cecchetti , m. borghesi , t. grismayer , e. dhumieres , p. antici , s. atzeni , p. mora , a. pipahl , l. romagnani , a. schiavi , y. sentoku , t. toncian , p. audebert , and o. willi , phys .
lett . * 99 * , 015002 ( 2007 ) .
h. w. powell , m. king , r. j. gray , d. a. maclellan , b. gonzalez - izquierdo , l. c. stockhausen , g. hicks , n. p. dover , d. r. rusby , d. c. carroll , h. padda , r. torres , s. kar , r. j. clarke , i. o. musgrave , z. najmudin , m. borghesi , d. neely , and p. mckenna , new j. phys . * 17 * , 103033 ( 2015 ) .
d. mangles , b. r. walton , m. tzoufras , z. najmudin , r. j. clarke , a. e. dangor , r. g. evans , s. fritzler , a. gopal , c. hernandez - gomez , w. b. mori , w. rozmus , m. tatarakis , a. g. r. thomas , f. s. tsung , m. s. wei , and k. krushelnick phys .
94 * , 245001 ( 2005 ) .
l. willingale , p. m. nilson , a. g. r. thomas , j. cobble , r. s. craxton , a. maksimchuk , p. a. norreys , t. c. sangster , r. h. h. scott , c. stoeckl , c. zulick , and k. krushelnick , phys .
* 106 * , 105002 ( 2011 ) .
l. willingale , a. g. r. thomas , p. m. nilson , h. chen , j. cobble , r. s. craxton , a. maksimchuk , p. a. norreys , t. c. sangster , r. h. h. scott , c. stoeckl , c. zulick , and k. krushelnick , new j. phys . *
15 * , 025023 ( 2013 ) .
j. peebles , c. mcguffey , c. m. krauland , l. c. jarrott , a. sorokovikova , m. s. wei , j. park , h. chen , h. s. mclean , c. wagner , m. spinks , e. w. gaul , g. dyer , b. m. hegelich , m. martinez , m. donovan , t. ditmire , s. i. krasheninnikov , and f. n. beg , nucl . fusion * 56 * , 016007 ( 2016 ) .
t. d. arber , k. bennett , c. s. brady , a. lawrence - douglas , m. g. ramsay , n. j. sircombe , p. gillies , r. g. evans , h. schmitz , a. r. bell , and c. p. ridgers , plasma phys .
fusion * 57 * , 113001 ( 2015 ) .
t. toncian , c. wang , e. mccary , a. meadows , a.v .
arefiev , j. blakeney , k. serratto , d. kuk , c. chester , r. roycroft , l. gao , h. fu , x. q. yan , j. schreiber , i. pomerantz , a. bernstein , h. quevedo , g. dyer , t. ditmire , and b. m. hegelich , matter and radiation at extremes ( 2016 ) [ http://dx.doi.org/10.1016/j.mre.2015.11.001 ] .
b. gonzalez - izquierdo , r. j. gray , m. king , r. j. dance , r. wilson , j. mccreadie , n. m. h. butler , r. capdessus , s. hawkes , j. s. green , m. borghesi , d. neely , and p. mckenna , nature physics ( 2016 ) [ doi:10.1038/nphys3613 ] .