intense laser - matter interaction provides a route to create energetic particles ( e.g. , electrons , ions , and photons ) using table - top equipment .
clusters , possessing the transparency of gas targets and the high charge density of solid targets , proved to be very efficient absorbers of laser light .
their small size , compared to laser wavelength and skin depth , avoids reflection of the laser beam at the cluster surface as well as the loss of hot electrons into the cold bulk .
in fact , almost 100% absorption of the laser light was reported in experiments with rare gas clusters @xcite .
useful reviews on the subject are refs .
@xcite .    upon irradiation of the rare gas clusters by intense laser light ,
electrons first absorb energy and leave their `` parent '' ions .
this is known as inner ionization , meaning that the electrons are still bound to the cluster but not necessarily to their `` parent '' ions .
the total electric field ( i.e. , laser plus space charge field ) inside the cluster leads to inner ionization up to high charge states not possible with the laser field alone ( ionization ignition @xcite ) . as the laser intensity during the pulse increases , these electrons absorb energy from the laser field and may leave the cluster leading to the positive charging of the cluster known as outer ionization .
thus outer ionization leads to a non - neutral plasma . with the increasing outer ionization
, the restoring force of the ions counteracts ionization ignition so that the latter stops at some point .
the net positive charge left behind finally explodes due to the coulomb repulsion and hydrodynamic pressure , leading to the conversion of electron energy into ion energy . typically mev ions and kev electrons @xcite
are measured in experiments .
it is clear from the described scenario that the understanding of the relevant mechanism(s ) of laser energy absorption leading to the heating of cluster electrons and outer ionization is of great importance for the development of a complete theoretical description .
laser energy absorption by electrons proceeds either through resonances ( linear or nonlinear ) or through non - adiabaticities ( all possible types of collisions ) .
all of these processes lead to dephasing of the current density with respect to the laser field , which , according to poynting s theorem , is a prerequisite for absorption .
collisional absorption via collisions of electrons with ions are of minor importance at near infrared wavelengths @xmath0 nm or greater @xcite whereas it is the dominant absorption mechanism at short wavelengths @xcite , not studied in this paper .
the finite size of the clusters suggests that `` collisions with the cluster boundary '' may be responsible for the energy absorption .
however , this viewpoint is misleading , as will be shown in sec .
[ sec2b ] . during the expansion of the ionic core
, the decreasing charge density @xmath1 leads to the decrease of the mie plasma frequency , @xmath2 ( atomic units @xmath3 are used unless noted otherwise ) .
for very short near infrared laser pulses @xmath4 can not meet the linear resonance ( t)=[lr ] , unless the cluster has sufficiently expanded ( typically after a few hundred femtoseconds ) .
linear resonance , well understood in theory , experiments , and simulations @xcite , is thus ruled out for very short pulses or during the early cycles of a long pulse laser - cluster interaction where ion motion is negligible . in this case ,
nonlinear resonance ( nlr ) , whose origin lies in the anharmonicity of the cluster potential , turns out to be the dominant collisionless absorption mechanism .
in fact , for very short linearly polarized ( lp ) laser pulses , it was clearly shown @xcite that essentially _ all _ electrons that contribute to outer ionization pass through the nlr , which was unequivocally identified as _ the _ collisionless absorption mechanism in the absence of linear resonance .
the eigenfrequency @xmath5 $ ] of a ( laser- ) driven oscillator in an anharmonic potential , being dependent on its excursion amplitude @xmath6 ( or the energy ) , may dynamically meet the nlr = .
[ nlr ] due to the many body nature of the interaction , the identification , the separation , and the interpretation of the absorption mechanisms in molecular dynamics or particle - in - cell ( pic ) simulations are often difficult . recently , a method of identification of the nlr in many - body simulations of rare gas clusters has been proposed @xcite .
the possible importance of nlr was also mentioned or discussed previously @xcite .
the rigid sphere model ( rsm ) @xcite where electrons and ions are modeled by homogeneously charged rigid spheres oscillating against each other is clearly an oversimplification of a real many - particle system such as a cluster .
however , it proves useful for estimating the order of magnitude of the absorbed energy as well as for the calculation of the laser intensity where energy absorption is most efficient , as will be shown in the present work .
moreover , it provides physical understanding and clearly displays nlr @xcite .
the heating of cluster electrons in circularly polarized laser fields has not yet received much attention , at least theoretically .
experiments with rare - gas clusters show almost no effect on the x - ray emission @xcite and ion energy distribution @xcite when laser light of different ellipticity is used .
theoretically , circular polarization is particularly interesting because the above mentioned `` collisions with the cluster boundary '' are strongly suppressed in this case .
hence one may expect energy absorption being less efficient .
nlr , on the other hand , occurs in both cases , and , in fact , the energy absorption turns out to be equally efficient .
moreover , the absorbed energy compares well with the rsm predictions .
the outline of the present paper is as follows : in sec .
[ sec2a ] we briefly review the nlr and the rsm . in sec .
[ sec2b ] the rsm is extended to circular polarization , where nlr is observed as well . in sec .
[ sec2c ] the rsm threshold intensities are calculated .
section [ sec3 ] is devoted to the pic @xcite results for both linear and circular polarization .
finally , we summarize our results in sec .
[ sec4 ] .
throughout this paper we use @xmath7-cycle laser pulses of near infrared wavelength @xmath8 nm and a fixed cluster radius @xmath9 nm unless stated otherwise .
nlr is a robust phenomenon that
qualitatively  is insensitive to cluster and laser parameters .
the ion motion can be neglected in the study of energy absorption in very short laser pulses .
thus the ions just form a static , positively charged background of spherical shape . for not too high laser intensity , the collective motion of the electrons can also be approximated by a homogeneous , rigid sphere of negative charge . in the simplest case the radii of ion and electron sphere
are assumed to be equal . in a more realistic model
the electron cloud expands @xcite .
however , the method of identification of the nlr used in this paper is independent of the degree of electronic expansion .
the center of mass of the electron - ion system is , in good approximation , located at the center of the ion sphere . in an oscillating laser field , the homogeneously
charged electron sphere oscillates back and forth against the positively charged ion sphere .
the equation of motion of the electron center of mass in a lp laser field , polarized along @xmath10 , can be written as + g(r ) = - , [ rsm_eom ] where @xmath11 is the excursion of the electron sphere , normalized to the cluster radius @xmath12 , @xmath13 , @xmath14 is the normalized time in units of the laser period , and the dimensionless electrostatic restoring force @xmath15 is g(r ) = ( ) ^2 \ {    ll r-+ & 0r2 + & r2 +    . .
[ gr ] the first term in the upper line of is the linear force when the displacement of the electron sphere is small , the next two terms are the nonlinear terms which appear due to the partial overlap of the electron cloud with the ion cloud .
the term @xmath16 in the lower line of is the coulomb force between the separated electron sphere and ion sphere .
the quantity on the right hand side of is the normalized driver strength , which is the quiver amplitude @xmath17 of a free electron in the laser field @xmath18 normalized to the cluster radius @xmath12 .    , a @xmath19 cycle laser field @xmath20 , and cluster radius @xmath21  nm . within a narrow field strength interval ( here @xmath22@xmath23 ) the absorbed energy per particle ( solid line ) increases by many orders of magnitude .
the dashed line represents the absorbed energy by a purely harmonic oscillator driven by the same laser field .
[ rsmi],scaledwidth=50.0% ]    earlier work @xcite showed that absorption of laser energy in the rsm is characterized by a threshold driver strength below which absorption is negligible ( harmonic regime ) and above which absorption is almost constant .
figure [ rsmi ] shows this threshold behavior for @xmath24 and a @xmath19-cycle @xmath25-pulse @xmath20 for @xmath26 .
the dashed line is the absorption corresponding to a driven , purely harmonic oscillator + ( ) ^2 = - [ harmonicrsm ] in the laser field @xmath27 .
the laser energy absorbed by a single electron in a @xmath28-cycle laser pulse of period @xmath29 is = - _ 0^n t v(t)e_(t ) dt . [ introdephase ] solving analytically for the velocity @xmath30 and integrating one finds for the absorbed energy per electron @xcite e_0 ^ 2 [ introabseng ] where @xmath31 ^ 2 [ \omegamie^2 - ( 1 - 1/n)^2 \omegalaser^2]^2 \simeq ( \omegamie^2 - \omegalaser^2)^4 $ ] ( for sufficiently large @xmath28 ) was used .
the analytical estimate is plotted in fig .
[ rsmi ] together with the absorbed energy obtained from the numerical solution of .
one sees that the absorbed energy jumps by many orders of magnitude to a higher value after crossing a threshold driver strength . the higher the cluster charge density is , the higher is this threshold driver strength .
the rigid sphere model shows this behavior of efficient absorption above the threshold driver strength at all cluster charge densities independent of the linear resonance condition @xmath32 , contrary to the nanoplasma model @xcite .
since the rigid sphere model does not necessarily requires expansion of the cluster for the efficient absorption of laser energy , it permits us to understand the behavior of energy absorption and the underlying mechanism for very short laser pulses .
equation can be formally rewritten as + ^2 = - .
[ nlho ] equation yields the instantaneous , scaled effective frequency squared ^2= = [ efffrequ ] which passes through unity at the nlr .
the right hand side of is the restoring force divided by the excursion of the electronic cloud , which in the case of harmonic motion would be the square of the characteristic frequency .
( red , upper solid line ) vs laser cycles above the threshold driver strength for @xmath33 .
here we take @xmath34 corresponding to a laser intensity @xmath35 , a @xmath19-cycle @xmath25-pulse @xmath20 of wavelength @xmath8 nm , and a cluster radius @xmath9 nm .
excursion @xmath36 ( green , middle solid line ) and energy of the electron sphere @xmath37 ( blue , lower solid line ) are included in the plot .
outer ionization ( i.e. , @xmath38 ) and occurrence of nlr @xmath39 ^ 2=1 $ ] always coincide ( dashed vertical line ) .
[ rsmii],scaledwidth=50.0% ]    figure [ rsmii ] shows a typical example of the temporal behavior of @xmath39 ^ 2 $ ] above the threshold driver strength for @xmath33 .
since @xmath40 , @xmath39 ^ 2 $ ] starts at this value and drops with increasing driver strength .
it passes through unity at the time indicated by the vertical line , and it is exactly at that time where the electron sphere becomes free ( outer ionization ) .
this incidence is clearly visible from the energy of the electron sphere , which passes through zero and the excursion as well .
outer ionization and occurrence of nlr happens for all driver strengths above the threshold whereas the resonance is never met below the threshold .
since the amplitude of the excursion of the electronic sphere depends upon the driver strength , the excursion amplitude should also be large enough so that the nlr is passed .
the decrease in the effective frequency with the increase of the amplitude of excursion of the electronic sphere in the force field can be understood by analyzing its motion in the corresponding anharmonic potential v(r ) = ^2 r^2 \ {    ll -+ & 0r2 + - & r2 . +    .
[ rsm_pot]the period @xmath29 of oscillation of the electronic sphere in the potential @xmath41 can be approximated by a perturbation series @xcite as @xmath42^k\ , \diff r}{\sqrt{e_\mathrm{tot } - v_0(r)}}.      \label{timeperiod}\ ] ] here , @xmath43 is the total energy of the electronic sphere in the ionic field , @xmath44 is the harmonic oscillator potential and @xmath45 is the perturbation to the potential with @xmath46 and @xmath47 .
the effective frequency is then @xmath48 .
for the excursion @xmath49 and cluster radius @xmath9 nm ( @xmath50a.u . ) we can consider @xmath51 and the approximate potential @xmath52 .
corrections up to @xmath53 yield @xmath54 .
the unperturbed period is @xmath55 and the successive corrections are @xmath56 , @xmath57 , @xmath58 , @xmath59 , @xmath60 , @xmath61 with @xmath62 , @xmath63 , @xmath64 , @xmath65 , @xmath66 and @xmath67 , respectively .
neglecting the higher order term of @xmath15 for @xmath68 , a simpler approximation to the effective frequency can be derived from by replacing the excursion @xmath69 with the excursion amplitude @xmath70 : _ ( ) ( 1 - /16)^1/2 ( 1 - /32 ) .
[ rsm_freq2 ]     vs the excursion amplitude @xmath70 of the electronic sphere for different cluster charge densities @xmath71@xmath72 .
the solid lines are computed from the numerical solution of .
the dashed lines are the corresponding analytical approximations from for @xmath73 up to 6 , and the dashed - dotted lines are from . for the charge density @xmath74 one
expects nlr ( horizontal dashed line ) to occur at the excursion @xmath75 .
[ rsmiii],scaledwidth=50.0% ]    figure [ rsmiii ] shows the effective frequency vs the excursion amplitude @xmath70 of the electronic sphere in the potential for various cluster charge densities @xmath76@xmath72 as calculated from the numerical solution of together with the approximations and .
the effective frequency as calculated from shows good agreement below the excursion @xmath77 and low charge densities ( e.g. , @xmath78 ) . for higher charge densities more corrections [ large number of @xmath73 values in ]
are needed , which are very much cumbersome to calculate .
although fits well with the numerical solution , @xmath79 becomes negative when @xmath80 .
however , the results using agree with the exact ones in fig .  [ rsmiii ] up to @xmath81 for the excursions of interest , i.e. , up to the point where the electronic sphere undergoes nlr . the variation of frequency with the excursion amplitude as shown in fig .
[ rsmiii ] explains why a threshold driver strength is required for an appreciable laser energy absorption as well as for the crossing of the nlr in the fig .
[ rsmii ] : only a driver exceeding a certain threshold field strength will lead to excursions compatible with the nlr condition @xmath82 .
clusters in a circularly polarized ( cp ) laser field received less attention in the literature .
it is not known _ a priori _ how the outer ionization and energy absorption by clusters depend on the laser polarization . in laser - atom interaction
the laser polarization has dramatic effects : since in cp the free electrons do not return to their parent atom all the atomic effects relying on rescattering such as high - order harmonic generation , high - order above - threshold ionization , and nonsequential ionization are strongly suppressed . in the context of clusters ,
the study of the absorption efficiency as a function of the laser polarization can help to discriminate among different absorption mechanisms .
for instance , if laser energy absorption was indeed due to `` collisions with the cluster boundary '' it would be suppressed in cp because the electrons mainly swirl around parallel to the `` cluster boundary '' rather than crossing ( and hence colliding ) with it . however
, as we will show , the absorption of laser energy is largely independent of the laser polarization , thus ruling out `` collisions with the cluster boundary '' as a meaningful absorption mechanism .
let us first extend the rsm to cp . in a cp laser field with electric field components in @xmath10- and @xmath83-direction , the equation of motion for the electronic sphere in the rigid sphere approximation of a cluster can be written as \
{    c _ x + _ y    } + \ {    cr_x + r_y    } = - \ {    c e_^x ( ) + e_^y ( )    } . [
rsm_eomcpl]here , @xmath84 , @xmath85 , @xmath86 , and @xmath87 one identifies @xmath88 and @xmath89 as the two components of the restoring force in .
note that we have divided the electric field components by a factor @xmath90 so that the ponderomotive energy @xmath91 , ( i.e. , the time - averaged quiver energy of a free electron in the laser field ) is the same as in the lp case with the same @xmath92 ( otherwise @xmath93 would be a factor of two higher in the cp case ) .
the square of the effective , time - dependent oscillator frequency in the cp laser field can be written as ^2 = = [ efffrequcpl ] , which has the same form as in the lp case .
( green , top left solid line ) vs laser cycles above the threshold intensity @xmath94 for circular polarization , @xmath33 , and cluster radius @xmath9 nm . here , @xmath95 , corresponding to a laser intensity @xmath96 , @xmath19-cycle @xmath25-pulse with components @xmath97 , @xmath98 , and wavelength @xmath8 nm .
excursion @xmath69 ( blue , middle left solid line ) and energy @xmath37 ( red , bottom left solid line ) of the electron sphere are also plotted . the energy @xmath37 is scaled down by a factor 10 to display within the excursion and the frequency range .
outer ionization ( i.e. , @xmath38 ) and occurrence of nlr @xmath39 ^ 2=1 $ ] always coincide ( dashed vertical ) .
[ rsmiicpl],scaledwidth=50.0% ]    earlier , in fig .
[ rsmii ] , it was shown that nlr and outer ionization in the rsm only occur when a threshold laser intensity is crossed .
the same is true for the occurrence of nlr with cp light .
figure [ rsmiicpl ] shows the temporal behavior of @xmath99 ^ 2 $ ] above the threshold driver strength for a @xmath19-cycle cp laser pulse of wavelength @xmath8 nm . here
the cluster charge density is @xmath72 times over - critical , i.e. , @xmath40 at which @xmath100 starts and drops with increasing driver field during the pulse .
the nlr @xmath101
^ 2=1 $ ] is passed at the time indicated by the vertical line . as in the lp case in fig .
[ rsmii ] , the electron sphere is set free at the time the nlr is passed : the energy of the electron sphere passes through zero , and the excursion sharply increases to a high value .
it is also clear from fig .
[ rsmiicpl ] , that once the electronic sphere is set free , the frequency drops to zero , and the total absorbed energy remains positive .
a zero effective frequency implies an infinite period , i.e , the electron sphere does not return to the ion sphere . the main difference between fig .
[ rsmiicpl ] and fig .
[ rsmii ] is that in the case of the cp laser field the decrease of the effective frequency is smooth ( i.e. , no oscillations ) since the electric field vector rotates but its absolute value remains constant . as a consequence ,
the electron sphere spirals out , staying away from the potential center where @xmath102 . for lp instead , the electron sphere is driven through the origin and hence the effective frequency undergoes oscillations before it drops to the resonance value , as visible in fig .
[ rsmii ] .
nlr is clearly identified in both cases .
nlr occurs above a threshold driver strength . beyond this driver strength the rigid electron sphere gains laser energy which is many order of magnitude higher than below the threshold ( see fig .
[ rsmi ] ) . in an open potential such as
the electron sphere is detached from the ion sphere above the threshold driver strength , i.e. , appreciable energy absorption and outer ionization occur simultaneously .
we would like to remark that in closed potentials nlr occurs as well , as discussed in ref .
@xcite .    the threshold driver strength can be estimated .
the dimensionless potential @xmath103 of the electron sphere in the ionic field can be written as @xmath104 application of a static electric field @xmath105 ( corresponding to the peak field strength of a low - frequency laser field ) , suppresses the potential in one direction by the amount @xmath106 .
the effective potential seen by the electron sphere is @xmath107 with @xmath108 .
the potential barrier vanishes if @xmath109 and @xmath110 , leading to @xmath111 , and the nlr threshold intensity is estimated to be     for linearly ( lp ) and circularly ( cp ) polarized laser fields .
results from the full numerical solution of with lp ( lower solid ) , cp ( top solid ) , the vanishing barrier approximation ( vba - lp , middle solid ) , and the vba corrected for cp ( upper dashed ) are shown ( see text for a discussion ) .
the over - the - barrier approximation ( oba ) ( lower dashed ) underestimates the exact threshold intensities .
[ rsmpert2],scaledwidth=50.0% ]    @xmath112    we call this the vanishing - barrier - approximation ( vba ) .    in atomic ionization the so - called over - the - barrier approximation ( oba ) or
bethe - rule @xcite allows to estimate at which electric field strength a certain atomic charge state dominates .
if one applies the oba to the rsm one obtains the two equations @xmath113 and @xmath114 with @xmath115 the barrier location .
one finds @xmath116 .
this gives the oba threshold intensity of the nlr @xmath117 which underestimates the numerically determined threshold intensity , as we will show now .
figure  [ rsmpert2 ] shows the threshold intensity as a function of the cluster charge density @xmath71@xmath118 ( corresponding to average charge states @xmath119 for xenon ) .
the numerically determined threshold intensities for lp and cp laser light show that when the cluster charge density is low , the nlr occurs almost at the same value of the threshold intensity , irrespective of the polarization .
as the charge density increases the nlr threshold intensity appears to be higher for cp than for lp .
the vba of the threshold intensity is in good agreement with the numerical result for lp whereas the oba underestimates it .
this fact might be related to the recently observed `` enhanced saturation intensities '' in the ionization of finite size systems such as c@xmath120 ( see , e.g. , @xcite and references therein ) , indicating that the latter might neither be a many - electron nor a quantum effect but just due to the finite size of the target .
the difference of the threshold intensities for lp and cp is due to the definition of the cp field ( [ ex ] ) , ( [ ey ] ) where a factor @xmath121 has been introduced in order to render the ponderomotive potential equal for lp and cp .
however , for the threshold intensity it is the electric field ( or the intensity ) that matters , not @xmath122 . for a given @xmath105
the laser intensity is @xmath123 in the lp case but only @xmath124 for cp .
therefore , the upper black , dashed line in fig .
[ rsmpert2 ] shows the vba threshold intensity multiplied by a factor of two , which is in good agreement with the numerical results for the cp laser field at higher charge densities .
so far we have studied the nlr absorption of laser energy in a simplified model system assuming an anharmonic potential generated by the ions in which the homogeneous and rigid electron cloud moves . in reality
, the potential builds up during the interaction with the laser pulse because of ionization .
the delicate interplay of inner ionization , energy absorption by various mechanisms , and outer ionization can be simulated using methods such as pic or molecular dynamics .
previous work @xcite , studying lp short laser pulses , showed that nlr can be clearly identified in such simulations as well .
in this section we present results obtained from three - dimensional pic simulations .
the cluster is exposed to @xmath19-cycle @xmath25-pulses @xmath125 of near infrared wavelength @xmath126 nm , i.e. , the total pulse duration is @xmath127 . since motion of ions
does not play an important role during the entire pulse the ions are assumed fixed , which ensures a well defined , constant mie frequency @xmath128 .
a pic electron has the same charge to mass ratio as a `` real '' electron , that is , @xmath129 in atomic units .
each pic electron moves under the influence of the external laser field and the space charge field @xmath130 due to the potential @xmath131 that is created by all charges ( mapped to the numerical grid ) .
hence the equation of motion of the @xmath132-th pic electron is _
i + ( _ i , t ) = -_(t).[eqpic ] equation is solved self - consistently for all pic electrons . clearly , @xmath133 depends on the position of _ all _ other particles @xmath134 as well . the pic simulation starts with the neutral cluster configuration so that @xmath135 .       vs
laser intensity for charge density @xmath72 times the critical density and a xe@xmath136 cluster ( @xmath137 electrons , mean ionic charge state @xmath138 , and a cluster radius @xmath9 nm ) .
solid lines are pic results , dashed lines are the rsm results for a @xmath19 cycle pulse @xmath125 of wavelength @xmath139  nm .
[ picrsm1],scaledwidth=50.0% ]    figure [ picrsm1 ] shows pic results for the kinetic energy ( ke ) and the total energy ( te ) absorbed per electron in units of @xmath122 , vs the peak laser intensity .
one sees that the absorbed energy per electron is on the order of @xmath122 .
however , the absorbed energy is nonlinear in @xmath122 and displays a maximum before it drops because of the saturation of outer ionization .
the depletion is due to the fact that at a given intensity most ( if not all ) electrons are removed from the cluster ( complete outer ionization ) .
further inner ionization would be required to generate `` fresh '' electrons that could continue to absorb energy .
the maxima in the pic absorption curves are located close to the threshold intensity predicted by the rsm .
the total absorbed energy around the maximum of the pic absorption curves ( rsm as well ) is on the order of @xmath140@xmath141kev which has been also reported in experiments of intense laser clusters interactions @xcite . in pic simulations
no sharp intensity threshold exists since each pic electron sees its own time - dependent field ( space charge field plus the laser field ) .
therefore sharp jumps ( as seen in the rsm , e.g. , fig .  [ rsmi ] ) are absent .
figure [ picrsm2 ] compares the kinetic and potential energies of a single pic electron and the electronic sphere in the rsm .
the energies are plotted vs the excursion @xmath36 ( along the direction of the laser field ) .
the peak laser intensity is @xmath142 w/@xmath143 .
we have checked that there are many pic electrons that leave the cluster similarly to the one presented in fig .
[ picrsm2 ] .
the potential and the kinetic energy gained by the electron sphere in the rsm is divided by the total number of pic electrons so that a quantitative comparison is possible .
the main difference between pic and rsm results is that in the pic simulation all the electrons see initially a zero potential since the potential builds up from the neutral cluster configuration in the course of the interaction with the laser field while in the rsm the electronic sphere oscillates in a prescribed potential .
those pic electrons that stay long inside the cluster experience almost the full ionic potential which then is similar to the rsm potential .     in the laser polarization direction .
the escaping pic electron resembles the equivalent electronic particle in the rsm .
peak laser intensity @xmath142 w/@xmath143 , charge density is @xmath144 .
other parameters as in fig .
[ picrsm1 ] .
[ picrsm2],scaledwidth=50.0% ]    as mentioned earlier , the absorption of energy by a pic electron depends upon the self - consistent potential which develops during the laser pulse due to outer ionization . as a result
different pic electrons move along different trajectories , `` see '' a different potential , and thus are set free at different times .     from the pic simulation vs the excursion @xmath36 in laser polarization direction .
yellow ( solid gray ) lines represent cuts of the potential @xmath145 at time @xmath146 for different @xmath83 and @xmath147 .
the circles ( red ) represent the total energy @xmath148 of individual pic electrons located within the simulation box at that time .
the dashed lines represent the kinetic ( green , dashed gray ) energy ( ke ) and potential ( dashed black ) energy ( pe ) of a pic electron that is outer ionized when the excursion @xmath149 meets the condition @xmath150 at the same time ( see fig .
[ rsmiii ] for the charge density @xmath74 ) .
the peak laser intensity is @xmath142 w/@xmath143 .
other parameters as in fig .
[ picrsm2 ] .
[ picpot],scaledwidth=50.0% ]    figure [ picpot ] shows a snapshot of the collective potential @xmath131 at time @xmath146 for a cut at @xmath147 ( and various @xmath83 throughout the cluster ) .
the most lower curve represents the potential for @xmath151 .
the red circles represent the total energy @xmath152 of individual pic electrons located within the simulation box at that time .
one observes that many pic electrons are accumulated near the bottom - left of the potential well at this time .
these pic electrons with @xmath153 remain bound since the laser field amplitude is already decreasing from the 4th cycle onward .
the kinetic and the total energy of a pic electron which leaves the potential well when its excursion becomes @xmath149 at the same time is also shown .
this excursion approximately satisfies the nlr condition @xmath154 for the charge density @xmath155 , as identified in fig .
[ rsmiii ] with the rsm analysis .
the nlr behavior exhibited by the pic electron in fig .
[ picpot ] is not accidental .
for the sake of an unequivocal and explicit identification of the nlr we now analyze the motion of all individual pic electrons in the same way as it has been done with the motion of the electron sphere in the rsm in sec .
recalling , the equation for the effective , time - dependent oscillator frequency analogous to for the @xmath132-th pic electron reads ^2_,i(t ) = -= [ efffrequpic ] .
earlier we have mentioned that @xmath133 depends on the position of _ all _ other particles @xmath134 and the simulation starts with the charge neutral cluster configuration i.e. , @xmath135 .
hence , a pic electron `` sees '' initially an effective frequency @xmath156 .
the laser field disturbs the charge equilibrium and @xmath157 becomes different from zero .
@xmath157 may be even negative in regions of accumulated electron density ( repulsive potential ) .
as the cluster charges up , @xmath158 quickly increases beyond unity ( where the rsm starts in the first place ) .
the start from @xmath156 , the possibility of negative @xmath157 , and the three - dimensionality are the main differences to the rsm analysis above .    , excursion @xmath159 , and the total energy @xmath152 for the pic electron of fig .  [ picrsm2 ] vs time in laser cycles .
the total energy becomes positive only when the nlr is crossed ( indicated by the vertical dashed line ) .
this result resembles the rsm result in fig .
[ rsmii ] when nlr is met .
the charge density is @xmath74 , the peak laser intensity is @xmath142 w/@xmath143 .
other parameters as in fig .
[ picrsm2 ] .
[ timevsomega ] , scaledwidth=50.0% ]    figure [ timevsomega ] shows the effective frequency squared , the total energy @xmath160 , and the excursion @xmath159 vs time for the pic electron whose energy history is shown in fig .
[ picrsm2 ] .
we define the time when , for a particular electron , @xmath161 becomes @xmath162 as the ionization time of that electron .
data are plotted shortly before the emission of this particle from the cluster potential .
it is clearly visible in fig .
[ timevsomega ] that the pic electron is escaped only when the resonance line @xmath163 is passed .
figure  [ timevsomega ] can be well compared with fig .
[ rsmii ] showing ionization of the rsm via nlr .    in the case of pic simulations
the fulfillment of the nonlinear resonance condition @xmath163 is necessary but not sufficient for ionization .
as the potential builds up , the pic electrons transiently meet the nlr condition , and , in fact , some electrons leave the cluster at that early stage when the potential is still shallow and the laser field is relatively weak .
however , as the potential deepens , pic electrons `` dropping '' below the energy necessary for nlr to occur , behave from then on similar to the rsm and may finally escape only by climbing up in the potential and hitting the nlr @xmath164 . during the 0.4 laser cycles plotted in fig .
[ timevsomega ] the @xmath158-curve displays artificial , short - time scale fluctuations inherent in pic simulations @xcite .
however , we checked that macroscopic observables such as the absorbed energy or the degree of outer ionization are well converged .    ,
( b ) @xmath165 , ( c ) @xmath166 , ( d ) @xmath167 , ( e ) @xmath168 , ( f ) @xmath169 , ( g ) @xmath170 , and ( h ) @xmath171 laser cycles for lp , laser intensity @xmath172 , and @xmath33 .
other parameters as in fig .
[ timevsomega ] .
the radial positions ( in units of @xmath12 ) are color - coded .
electrons become free upon crossing the nlr , i.e. , @xmath173 .
[ picweffvsenglp],scaledwidth=55.0% ]    by following the dynamics of the electrons in the effective frequency vs energy - plane we identify the main pathway to outer ionization and efficient absorption .
figure  [ picweffvsenglp]a  h shows the scaled effective frequencies squared @xmath158 of the individual pic electrons vs their energies @xmath174 they would have if the driver is switched off instantaneously at @xmath175 , @xmath176 , @xmath177 , @xmath178 , @xmath179 , @xmath180 , @xmath181 , and @xmath182 laser cycles , respectively . at the time when , for a particular electron , @xmath183 becomes @xmath162 the ionization occurs for that electron .
the laser intensity is @xmath184 , and the pre - ionized cluster is @xmath72 times over - critical so that @xmath33 .
as is clearly visible in fig .
[ picweffvsenglp ] , each electron reaches positive energy close to the point @xmath173 .
the radial position of each electron is color - coded , indicating that outer ionization occurs at radii around @xmath185 . during the early time of the laser pulse ( fig .
[ picweffvsenglp]a , b ) when many electrons are still inside the cluster , @xmath158 spreads over a wide range , starting from the maximum value @xmath186 down to negative values due to the repulsive force exerted by the compressed electronic cloud .
note that negative values in effective frequency occur mainly at early times where most of the electrons are still inside the cluster .
electrons with positive but very small @xmath183 and @xmath187 represent low energetic electrons removed earlier during the pulse ( see fig .
[ picweffvsenglp]a , b ) .
the occurrence of nlr is less clear for these early leaving electrons . as mentioned above
, these electrons move in a shallow effective potential with @xmath188 when they leave the cluster with ease and with rather low kinetic energy because the laser intensity is still low at the time of their emission .
figures  [ picweffvsenglp]c  f show that most of the electrons escape from the cluster by passing through the channel @xmath173 at radii around @xmath185 .
it is also visible that more and more electrons are driven to positive frequency before they leave the cluster by passing through @xmath173 .
this is so because as more and more electrons are freed , the remaining electrons experience predominantly the force by the ionic background , and they move deep into the potential ( see their negative values in energy ) where they experience the full mie - frequency @xmath189 . in fig .
[ picweffvsenglp]c , d , the few electrons with positive energy but small radii are those driven back to the cluster by the laser field . in figures  [ picweffvsenglp]e , f , electrons are strongly aligned ( no scattered points ) since the laser field is approaching its maximum ( at @xmath190-cycles ) . after the peak of the laser pulse ( figs  [ picweffvsenglp]g ,
h ) the restoring force of the ions on almost all electrons dominates the laser force .
since many of the features of energy absorption and nlr in a lp laser field are common to the case of cp , we only point out the main differences .    , excursion @xmath159 , and the total energy @xmath152 for a pic electron in a cp field vs time in laser cycles .
the total energy becomes positive only when the nlr is crossed ( indicated by vertical , dashed line ) .
this result resembles the rsm result in fig .
[ rsmiicpl ] .
the charge density is @xmath74 , the peak laser intensity is @xmath142 w/@xmath143 .
[ timevsomegacpl ] , scaledwidth=50.0% ]    equation holds in the cp laser field as well .
figure [ timevsomegacpl ] shows the effective frequency squared vs time for one of the pic electrons , together with the total energy @xmath160 and the excursion @xmath159 .
one can see that the pic electron is freed ( i.e. , its total energy becomes positive ) only when the resonance line @xmath163 is passed .    , ( b ) @xmath165 , ( c ) @xmath166 , ( d ) @xmath167 , ( e ) @xmath168 , ( f ) @xmath169 , ( g ) @xmath170 , and ( h ) @xmath171 laser cycles for a laser intensity @xmath172and @xmath33 . other parameters as in fig .
[ picrsm1 ] . the radial positions ( in units of @xmath12 )
are color - coded .
electrons become free upon crossing the nlr , i.e. , @xmath173 .
[ picweffvsengcpl],scaledwidth=55.0% ]    figure  [ picweffvsengcpl ] is the cp analogue of fig .
[ picweffvsenglp ] .
the results are very much similar to the lp case shown in fig .
[ picweffvsenglp ] and the arguments made there apply here as well .
nlr is clearly observed .
the main difference is that the pic electrons are nicer aligned towards the resonance point , even at early times during the laser pulse ( see fig .
[ picweffvsengcpl ] b , c , d ) .
almost no scattered particles are visible because the dynamics mainly consist of swirling around the cluster center rather than oscillating through it .
the number of electrons returning to the cluster is much less so that the recombination and rescattering probability is smaller in the case of cp .
the same is observed in laser - atom interaction experiments , with important consequences for harmonic generation and non - sequential ionization .    for sufficiently high charge density
, the degree of outer ionization ( defined as the ratio of removed electrons to the total number of electrons ) is less compared to the case of lp laser pulses .
the explanation is the same as in sec.[sec2c ] : the effective laser intensity for cp is by a factor of two lower than for lp because of our definition of the cp field ( [ ex ] ) , ( [ ey ] ) .
one may object that , since the denominator in necessarily increases while the numerator decreases for an electron on its way out of the cluster potential , that the passage through a point @xmath191 with @xmath10 _ some _ value @xmath192 is rather the consequence of outer ionization than the mechanism behind it .
however , nlr only occurs at @xmath193 , and the results in fig .
[ picweffvsenglp ] and fig .
[ picweffvsengcpl ] show only little spreading along @xmath158 at @xmath194 . moreover , the fact that _ both _ the single electron energies become positive _ and _ the radii exceed @xmath195 when @xmath196 indicates that nlr is indeed the responsible mechanism behind outer ionization accompanied by efficient absorption of laser energy .     vs
laser intensity for charge densities ( a ) @xmath197 ( linear resonance ) , ( b ) @xmath198 , ( c ) @xmath199 , and ( d ) @xmath200 .
the absorbed energies for lp ( red , lighter gray ) and cp ( blue , darker gray ) using pic ( dashed ) and the rsm ( solid ) are shown .
[ rsmengcpl],title="fig:",scaledwidth=48.0% ] +   vs laser intensity for charge densities ( a ) @xmath197 ( linear resonance ) , ( b ) @xmath198 , ( c ) @xmath199 , and ( d ) @xmath200 .
the absorbed energies for lp ( red , lighter gray ) and cp ( blue , darker gray ) using pic ( dashed ) and the rsm ( solid ) are shown .
[ rsmengcpl],title="fig:",scaledwidth=48.0% ]    finally , fig .
[ rsmengcpl ] shows the average value of the total absorbed energy per electron vs the peak laser intensity for cluster charge densities between @xmath201@xmath72 in cp and lp laser fields .
pic results are compared with the rsm absorption results . the absorbed energy per electron in figs .
[ rsmengcpl]a  d is plotted in units of @xmath202 whereas the same results are shown in figs .
[ rsmengcpl]e  h in units of the ponderomotive energy @xmath93 .
the pic results in figs .
[ rsmengcpl]a  d show that the absorbed energy increases linearly in the log - log representation up to a certain intensity and then tends to saturate due to the saturation of outer ionization .
one sees that the saturation in the pic results occur close to the rsm threshold intensity . since with increasing charge density the restoring force due to the ions increases ,
the saturation of energy absorption in the rsm and pic occur at higher laser intensities as the density increases from @xmath201 to @xmath74 in fig .
[ rsmengcpl]a  d . when outer ionization and the energy absorption saturate with increasing peak laser intensity , the average absorbed energy per electron divided by @xmath93 ( which is proportional to the so - called fractional absorption ) starts decreasing in figs .
[ rsmengcpl]e  f .
figures [ rsmengcpl]a , e show that at linear resonance @xmath197 , the absorbed energy is already high at low values of the laser intensity ( @xmath203 ) , and absorption is very efficient as compared to higher charge densities @xmath198@xmath72 , presented in figs .
[ rsmengcpl]b  d and figs .
[ rsmengcpl]f  h . in fact ,
[ rsmengcpl]e illustrates that the absorbed energy is on the order of @xmath204 ( both in the rsm and in the pic ) before the saturation of outer ionization .
however , one should bear in mind that at too low laser intensities inner ionization would not occur in the first place so that in reality there would be no absorption at all .
energy absorption in cp and lp laser fields at all intensities and all charge densities are almost equally efficient . recalling that with our definition of the cp laser field the ponderomotive potential is equal for lp and cp while the electric field amplitude is not ,
we conclude that for the absorbed energy @xmath122 matters while the nlr threshold is determined by the field strength .
the discrepancy between the pic and rsm results at lower intensities is , again , due to the fact that the rsm provides a very well defined potential , although a more and more shallow one as the density decreases , in which the rigid electron sphere moves . in pic
, the potential is initially zero and builds up in the course of outer ionization . at higher intensities ,
the potentials for the remaining pic electrons are closer to the rsm potential ( see fig .  [ picrsm2 ] ) , explaining the improving agreement between rsm and pic results as the laser intensity increases .
in summary , two different approaches to study collisionless laser energy absorption by clusters , namely ( i ) the rigid sphere model and ( ii ) particle - in - cell simulations , were pursued in this work .
the goal was to identify the dominant mechanism of energy absorption and outer ionization of the cluster electrons in near infrared , short laser pulses where collisional absorption is known to be inefficient .
we showed that the cluster electrons contributing to efficient absorption and outer ionization undergo nonlinear resonance , meaning that the instantaneous frequency of their motion in a time - dependent , anharmonic , effective potential transiently meets the laser frequency .
nonlinear resonance is the only possible absorption mechanism if the laser pulse is too short for the linear resonance to occur ( or during the early cluster dynamics in longer pulses ) and if electron - ion collisions ( inverse bremsstrahlung ) are negligible . in order to prove the occurrence of nonlinear resonance we used a method to analyze the results obtained from particle - in - cell simulations , namely the mapping of the system of electrons and ions that interact through their mean field onto a system of nonlinear oscillators .
the occurrence of nonlinear resonance in the particle - in - cell simulations presented in this paper resembles the nonlinear resonance in the rigid sphere model . for a given cluster charge density
, there is a threshold intensity around which the average electron energy displays a maximum conversion of laser energy .
the threshold intensity can be calculated using the newly introduced vanishing barrier approximation .
the common over - barrier approximation
applicable to atoms  fails in the case of finite - size potentials and underestimates the required laser field strength for ionization .
this fact might be related to the experimentally observed `` enhanced saturation intensity '' in the ionization of finite - size systems such as c@xmath120 .
the efficiency of energy absorption from the laser is almost the same for linear and circular polarization , indicating that `` collisions with the cluster boundary '' as an explanation are misleading .
instead , nonlinear resonance is the main absorption mechanism in both cases . to illustrate this , the rigid sphere model has been extended for the circularly polarized laser pulses in this paper .
t.  ditmire , r.a .
smith , j.w.g .
tisch , and m.h.r .
hutchinson , , 3121 ( 1997 ) .
molecules and clusters in intense laser fields , edited by jan posthumus ( cambridge university press , cambridge , u.k . , 2001 ) .
u.  saalmann and j.m .
rost , j.  phys .
b :  at .
phys .  * 39 * , r39-r77 ( 2006 ) .
krainov and m.b .
smirnov , phys .
370 * , 237 ( 2002 ) .
c.  rose - petruck , k.j .
schafer , k.r .
wilson , and c.p.j .
barty , , 1182 ( 1997 ) .
d.  bauer and a.  macchi , , 033201 ( 2003 ) .
t.  ditmire , j.w.g .
tisch , e.  springate , m.b .
mason , n.  hay , j.  marangos , and m.h.r .
hutchinson , nature ( london ) * 386 * , 54 ( 1997 ) .
v.  kumarappan , m.  krishnamurthy , and d.  mathur , , 043204 ( 2003 ) .
e.  springate , s.a .
aseyev , s.  zamith , and m.j.j .
vrakking , , 053201 ( 2003 ) .
k.  ishikawa and t.  blenski , , 063204 ( 2000 ) .
f.  megi , m.  belkacem , m.a .
bouchene , e.  suraud , and g.  zwicknagel , j.  phys .
phys .  * 36 * , 273 ( 2003 ) .
d.  bauer , j.  phys .
* 37 * , 3085 ( 2004 ) . c.  siedschlag and j.m .
rost , , 043402 ( 2004 ) . t.  ditmire , t.  donnelly , a.m.  rubenchik , r.w .
falcone , and m.d .
perry , , 3379 ( 1996 ) .
t.  dppner , th .
fennel , th .
diederich , j.  tiggesbumker , and k.h .
meiwes - broer , , 013401 ( 2005 ) .
l.  kller , m.  schumacher , j.  khn , s.  teuber , j.  tiggesbumker , and k.h .
meiwes - broer , , 3783 ( 1999 ) .
s.  zamith , t.  martchenko , y.  ni , s.a .
aseyev , h.g .
muller , and m.j.j .
vrakking , , 011201(r ) ( 2004 ) .
i.  last and j.  jortner , , 2215 ( 1999 ) .
ulf saalmann and jan - michael rost , , 223401 ( 2003 ) .
fennel , g.f .
bertsch , and k.h .
meiwes - broer , eur .
j.  d * 29 * , 367 ( 2004 ) . c.  siedschlag and j.m .
rost , , 031401 ( 2005 ) .
t.  martchenko , ch .
siedschlag , s.  zamith , h.  g.  muller , and m.  j.  j.  vrakking , , 053202 ( 2005 ) .
t.  taguchi , t.m .
antonsen , jr . , and h.m .
milchberg , , 205003 ( 2004 ) ; thomas m.  antonsen , jr .
toshihiro taguchi , ayush gupta , john palastro , and howard m.  milchberg , phys .
plasmas * 12 * , 056703 ( 2005 ) .
mulser and m.  kanapathipillai , , 063201 ( 2005 ) ; p.  mulser , m.  kanapathipillai , and d.h.h .
hoffmann , , 103401 ( 2005 ) .
d.  bauer and p.  mulser , j.  phys .
ser .  * 11 * , 169 ( 2005 ) .
korneev , s.v .
popruzhenko , d.f .
zaretsky , w.  becker , laser phys .
* 2 * , 452 ( 2005 ) .
parks , t.e .
cowan , r.b .
stephens , and e.m .
campbell , , 063203 ( 2001 ) .
fomichev , s.v .
popruzhenko , d.f .
zaretsky , and w.  becker , j.  phys .
.  phys .
* 36 * , 3817 ( 2003 ) .
lin , h.  chu , m.  shen , y.  xiao , c.h .
lee , s.  chen , and j.  wang , optics communications * 231 * , 375 ( 2004 ) v.  kumarappan , m.  krishnamurthy , d.  mathur , and l.c .
tribedi , 023203 ( 2001 ) .
s.  ter - avetisyan , m.  schnrer , h.  stiel , u.  vogt , w.  radloff , w.  karpov , w.  sandner , and p.v .
nickles , , 036404 , ( 2001 ) .
m.  krishnamurthy , d.  mathur , and v.  kumarappan , , 033202 ( 2004 ) .
birdsall and a.b .
langdon , _ plasma physics via computer simulation _
( iop , bristol , 1991 ) . g.l .
kotkin and v.g .
serbo , _ collection of problems in classical mechanics _
( pergamon , oxford , 1971 ) .
bethe and e.e .
salpeter , _ quantum mechanics of one- and two - electron atoms _ ( plenum publishing corporation , new york , 1977 ) .
hertel , t.  laarmann , and c.p .
schulz , adv .  at .
* 50 * , 219 ( 2005 ) .
a.  jaro - becker , a.  becker , and f.h.m .
faisal , , 143006 ( 2006 ) .