arxiv: v1 [physics.plasm-ph] 11 Sep 2018

Transcription

1 Identifiction of coupling mechnisms etween ultrintense lser light nd dense plsms rxiv: v1 [physics.plsm-ph] 11 Sep 218 L. Chopineu, 1 A. Lelnc, 1 G. Blclrd, 1, 2 A. Denoeud, 1 M. Thévenet, 2 J-L. Vy, 2 G. Bonnud, 1 Ph. Mrtin, 1 H. Vincenti, 1, nd F. Quéré 1, 1 LIDYL, CEA, CNRS, Université Pris-Scly, CEA Scly, Gif-sur-Yvette, Frnce 2 Lwrence Berkeley Ntionl Lortory, Berkeley, CA 9472, USA (Dted: Septemer 12, 218) The interction of intense lser ems with plsms creted on solid trgets involves rich nonliner physics. Becuse such dense plsms re reflective for lser light, the coupling with the incident em occurs within thin lyer t the interfce etween plsm nd vcuum. One of the min prdigms used to understnd this coupling, known s Brunel mechnism, is expected to e vlid only for very steep plsm surfces. Despite innumerle studies, its vlidity rnge remins uncertin, nd the physics involved for smoother plsm-vcuum interfces is uncler, especilly for ultrhigh lser intensities. We report the first comprehensive experimentl nd numericl study of the lser-plsm coupling mechnisms s function of the plsm interfce steepness, in the reltivistic interction regime. Our results revel cler trnsition from the temporlly-periodic Brunel mechnism to chotic dynmic ssocited to stochstic heting. By reveling the key signtures of these two distinct regimes on experimentl oservles, we provide n importnt lndmrk for the interprettion of future experiments. High-density plsms cn e creted y focusing intense lser pulses on initilly-solid trgets. The interction of such plsms with lser light, investigted for severl decdes, involves rich non-liner physics which is not only of fundmentl interest, ut lso of high relevnce for wide rnge of pplictions over lrge intervl of lser intensities, nd spnning thermonucler fusion [1], lortory strophysics [2], or lser-driven prticle ccelertion [3, 4]. For most if not ll pplictions, depositing lser energy into the plsm is essentil. Due to their high density, lrgely in excess of the so-clled criticl density where the locl electron plsm frequency equls the lser frequency, these plsms however tend to reflect lrge frction of the lser light. The ctul coupling with the incident light field cn only occur either in the undercriticl prt of the density grdient t the plsm-vcuum interfce, where the lser wve propgtes, or within the skin depth of the overcriticl plsm, where the lser wve is evnescent. At such plsm densities, physicists initilly nticipted the min mechnism of energy deposition to e collisionl sorption [5]: electron-ion collisions disrupt the regulr quivering motion of the plsm electrons in the light field, sttisticlly leding to net kinetic energy gin from the lser field. Soon fter the invention of lsers, experimentl studies on the fesiility of lser-driven thermonucler fusion however reveled the importnce of non-collisionl light sorption mechnisms, coming into ply for moderte lser intensities (Iλ W µm 2 /cm 2 ) [6]. For interctions t the surfce of dense plsms, these processes re expected to e most relevnt when the lser em impinges the henri.vincenti@ce.fr fien.quere@ce.fr trget t olique incidence (ngle of incidence θ i ) nd in p-polriztion, such tht n electric field component efficiently drives electron motion long the norml to the trget surfce. Among these so-clled nomlous sorption mechnisms, the resonnt excittion nd susequent dmping of collective electronic plsm wves t the criticl plsm density hs een the first key process to e identified oth theoreticlly nd experimentlly, nd is commonly known s resonnt sorption [7 1]. Since the coupling with lser light occurs t the interfce of the plsm with vcuum, the chrcteristic sptil length L of the plsm density grdient cross this interfce is crucil prmeter. This density grdient is generlly not step-like, due to the unvoidle expnsion of the plsm into vcuum, either during the min lser pulse driving the interction, or even efore this pulse when lser pre-pulses re present (either ccidentlly, or voluntrily introduced). In n influentil pper from 1987, Brunel predicted trnsition from resonnt sorption to new coupling mechnism, tht he ironiclly clled not-so-resonnt, resonnt sorption [11], depending on the vlue of L. He nticipted this mechnism to come into ply when the lser intensity ecomes so strong tht the quivering mplitude of the plsm electrons in the field long the surfce norml gets lrger thn L, such tht the plsm-vcuum interfce cn e modeled s step-like. This simple nd intuitive mechnism, now known s Brunel sorption or vcuum heting, is qulittively nlogous to the intensively studied three-step model of tomic nd moleculr strong-field physics, where n intense lser field (Iλ W µm 2 /cm 2 ) drives the recollision of ionized electrons with their prent ions [12]. Here, in ech opticl cycle, electrons t the trget surfce re drgged out of the plsm into vcuum when the component of the lser E-field norml to the surfce points to the trget (see simultion results in Fig.1).

2 2 x/λ log(n/n c ) t/t p x (m e c) FIG. 1. Temporl dynmics of dense plsm exposed to n ultrintense lser field in the Brunel regime. This grph displys results from Prticle-In-Cell simultion performed for = 2 (I = W/cm 2 for λ = 8 nm), θ i = 55 o, nd density grdient scle length L = λ/1 (with λ the lser wvelength). The gry-scle color mp shows the temporl evolution, during three lser opticl periods T, of the plsm electron density round the trget surfce, while the purple color scle shows the ttosecond light pulses emitted y this plsm (hrmonics 8 to 15). Two representtive trjectories for the prticles forming the expelled electron em (lue color) nd the recolliding electron flux (red color) re lso displyed. The right pnel shows typicl distriution of electrons in the x p x phse spce (x sptil coordinte long the trget norml) t time t/t = Lter in the cycle, when this driving field chnges sign, some of these expelled electrons re pushed ck towrd the prent plsm : s they penetrte into this dense plsm, they escpe the influence of the lser field due to plsm screening, nd propgte llisticlly into the trget (red trjectory in Fig.1). The initil model y Brunel for step-like surfce of perfectly-conducting plsm focused on this returning electron popultion. However, numericl simultions show tht for smoother interfces nd/or higher lser intensities, nother frction of the electrons escpes into vcuum (lue trjectory in Fig.1), typiclly in the form of periodic ttosecond unches. In oth cses, the fst electrons resulting from this suopticl-cycle dynmics crry wy energy cquired from the lser: for convenience, the terms Brunel electrons or Brunel sorption used in this pper will encompss these two popultions. A few yers lter, with the development of high-power femtosecond lsers, Brunel mechnism ppered s n idel toy-model to understnd the interction of dense plsms with these ultrshort pulses. First, their ultrhigh intensities result in lrge electron quivering mplitudes. Second, these pulses re so short tht plsm expnsion during the interction is very limited, potentilly leding to shrp density grdients t the plsm surfce, i.e. smll vlues of L. However, from n experimentl point of view, reching this regime turned out to e much more chllenging thn expected. The key difficulty rose from the unvoidle light pedestl hed of ultrshort lser pulses [13]: this pedestl, if too intense, leds to the premture cretion nd expnsion of the plsm, nd thus to long nd lrgely uncontrolled density grdients t the plsm-vcuum interfce when the min lser pulse hits the trget, mking the Brunel regime inccessile. More generlly, this mjor issue hs considerly complicted the interprettion of most erly experiments on the interction of intense ultrshort lsers with dense plsms. It took more thn n dditionl decde to find methods to efficiently reduce the light pedestl ccompnying ultrshort lser pulses [14 17], nd thus otin temporl contrsts tht t lst mde extremely shrp plsm surfces ccessile nd comptile with ultrhigh lser intensities [18, 19]. Nowdys, Brunel mechnism is most likely t ply in experiments performed on solid trgets with ultrintense lser pulses of suitly high-contrst. Yet, direct experimentl evidence is still elusive, nd its rnge of vlidity is not precisely known so fr. Furthermore, following the historicl development of this topic, the common wisdom tends to e tht when L is incresed, trnsition from Brunel to resonnt sorption should t some point occur [2 22], ut no cler experimentl evidence of this trnsition hs een reported yet. A rod rnge of topicl experiments re now performed worldwide on the interction of ultrintense lser pulses (Iλ 2 > 1 18 W µm 2 /cm 2 ) with dense plsms, driven y pplictions such s lser-driven ion [3, 4, 23, 24] nd electron ccelertion [25 32], or the genertion of intense hrmonics nd/or ttosecond light pulses [22, 33]. Clerly identifying the lser-plsm coupling mechnisms t ply in this interction regime, nd determining the rnge of physicl prmeters where they re relevnt, is essentil for the proper understnding of such experiments. This is wht we chieve in this rticle, y focusing ultrintense femtosecond lser pulses on dense plsm with shrp, controlled nd mesured density grdient scle length L, which we systemticlly vry from L λ to L λ (with λ the lser wvelength). We show how performing nd correlting mesurements of the high-order hrmonics nd reltivistic electrons emerging from the trget provide cler signtures of these couplings mechnisms, nd relte these oservtions to the underlying physics through n dvnced nlysis of 2D nd 3D Prticle-In-Cell (PIC) simultions, solving the coupled Vlsov-Mxwell eqution system. This comprehensive study shows tht Brunel mechnism is indeed the relevnt physicl process for shrp enough plsm-vcuum interfces. As expected, trnsition occurs to different mechnism when the density grdient scle length L is incresed. Mesurements of this trnsition s function of the lser incidence ngle provide confirmtion of Brunel s trnsition criteri sed on the comprison of the electron quivering mplitude with the typicl sptil extent of the interfce. However, we estlish tht in the regime of ultrhigh lser intensities considered here, resonnt sorption plys no significnt role in the regime of lrge L (L λ). The coupling is rther dominted y nother kinetic mechnism, so fr known s stochstic heting, in which collective plsm

3 3 Lnex Lser τ Trget XUV grting MCP θ i ω L ~ λ / 15 fc/mrd u (mrd) c L ~ λ / 1.5 fc/mrd u (mrd) (mrd) ω/ω L 9 (mrd) ω/ω L (mrd) 2 E (MeV) (mrd) (mrd) (mrd) 11 8 E (MeV) FIG. 2. Principle of the experiment nd min experimentl findings. Pnel shows sketch of the experiment. The trget is impinged y controlled prepulse, followed y the min pulse fter n djustle dely τ. Two of the min dignostics re displyed, the lnex screen for the mesurement of the sptil profile of the high-energy electron em, nd the ngulrly-resolved xuv spectrometer. These two dignostics cn e used either seprtely, or simultneously when smll holes re mde in the electron detection ssemly, s shown in the figure. They cn lso e replced y n ngulrly-resolved electron spectrometer. The min experimentl findings for p-polrized lser field re summrized in the two lower pnels nd c: left imges, ngulr emission pttern of reltivistic electrons; ottom right imges, ngulrly-resolved energy spectrum of electrons in the incidence plne ( = ); top right imges, ngulrly-resolved hrmonic spectrum. These illustrte the mjor chnges occurring on these three oservles s the density grdient scle length is incresed from L λ to L λ (θ i = 55 o, = 3.5, τ = 1 ps for pnel, leding to L = λ/15, τ = 1 ps for pnel c, leding to L = λ/1.5). These very contrsted fetures constitute signtures of the different underlying lser-plsm coupling mechnisms. effects ply little role: s initilly suggested in Ref. [34], electrons in the underdense prt of the density grdient gin energy in the interference pttern resulting from the superposition of the incident lser field with the field reflected y the overdense prt of the plsm. It hs een estlished theoreticlly tht t the lser intensities considered here, electron dynmics in such n interference pttern is not integrle, gets chotic, nd cn led to high energy trnsfer from the lser wve to the electron popultion [35]. In this pper, the mplitude of the incident lser field, which determines the intensity on trget, is chrcterized y the dimensionless potentil vector t the pek of the pulse, = ee /mcω = λ[µm](i[w/cm 2 ]/ ) 1/2, with c the speed of light, e the elementry chrge, m the electron mss, ω the lser frequency, nd E the mplitude of the lser electric field. All experiments nd simultions presented here hve een performed with > 1, which corresponds to the interction regime where reltivistic effects ply n importnt role on electron motion. We define n c s the criticl plsm density ssocited to the lser frequency ω (n c = mɛ ω 2 /e 2 = cm 3 for 8 nm lser light). I. DESCRIPTION OF THE EXPERIMENT A sketch of the experiment is presented in the upper prt of Fig.2. A high-power femtosecond lser em is focused on silic trget, which it fully ionizes on thin surfce lyer, thus producing dense plsm (mximum plsm density n cm 3, i.e 4n c for 8 nm light). We use the UHI1 lser t CEA Scly, commercil system delivering 2 to 25 fs pulses with pek power of 1 TW. After correction of its wvefront y n dptive optic system, the em is focused y n off-xis prolic mirror with f-numer of f/6, leding to focl spot of 5 µm dimeter (FWHM in intensity) nd to n estimted pek intensity of W/cm 2 ( 3.5) on trget. By defult, the lser em is p-polrized on trget, ut the polriztion cn e switched to s y inserting thin zeroorder mic hlf wve-plte in the em.

4 4 The first key spect of the experiment is tht it ws crried out with high degree of control nd n ccurte knowledge of the plsm density grdient scle length L t the trget surfce, which is prerequisite for the study of the lser-plsm coupling mechnisms. This implies the use of lser pulses of very high temporl contrst, so tht premture cretion of the plsm on trget is voided: this is chieved thnks to doule plsm mirror system plced efore the min experimentl chmer [18], which increses the contrst y out 4 decdes, from 1 9 to 1 13 on 1 ps time scle. This ultrhigh temporl contrst is of prmount importnce for ll experiments presented here. Strting from the very steep density grdient llowed y this ultrhigh temporl contrst, L is then vried in controlled wy, thnks to the introduction of wek prepulse (fluence 1 kj/cm 2 ) t n djustle dely τ efore the min pulse ( τ 15 ps), produced from n edge of the min em using the opticl lyout descried in Ref.[36]. This prepulse is strong enough to ionize the trget nd initite plsm expnsion, t typicl velocity in the 4 to 6 nm/ps rnge. For ll experimentl conditions considered herein, nd in prticulr for ll incidence ngles, we systemticlly mesured L(τ) using the recently-introduced technique of Sptil Domin Interferometry [37] (SDI, see supplementry mteril). The second key spect of the experiment is the comintion of dignostics tht were implemented to study the interction. We concentrted on two types of oservles: the reltivistic electron em emitted y the trget towrds vcuum, nd the em of high-order hrmonics generted round the speculr reflection direction. Two dignostics were used for the electron em. First, lnex screen, plced ehind 13 µm thick luminum foil (to eliminte lser light nd its hrmonics) nd 2 mm thick glss plte (to filter out low-energy electrons), nd which fluorescence ws imged on CCD cmer, provided the sptil profile of the emission of electrons with energies higher thn 1 MeV, t distnce of 1 cm from the trget (left imges in Fig.2-c). Second, we designed new type of mgnetic spectrometer for reltivistic electrons (see supplementry mteril), which provided, for every lser shot, the ngulrly-resolved kinetic energy spectrum of electrons, in the incidence plne (i.e. for = ). This (, E) distriution, with the ngle in the plne of incidence nd E the electron kinetic energy, ws mesured with very lrge ngulr cceptnce of = 5 mrd round the speculr reflection direction ( = ) (ottom right imges in Fig.2-c). The hrmonic em ws chrcterized using n ngulrly-resoled XUV spectrometer [36], with n ngulr cceptnce of 2 mrd round the speculr direction (top right imges in Fig.2-c). The hrmonic spectrum nd the electron em sptil profile were initilly mesured on the sme lser shots, thnks to smll holes in the luminum foil, glss plte nd lnex screen tht let the hrmonic em go through (Fig.2). However, we oserved n excellent shot-to-shot reproduciility of the experimentl results, so tht these multiple dignostics were finlly implemented on different lser shots performed under identicl interction conditions. A simple dditionl dignostic, implemented on seprte lser shots, consisted in mesuring the sptil profile of the lser em reflected y the trget, y inserting frosted glss plte in this em 2 to 3 cm fter the trget, nd mesuring the imge of the lser light scttered y this plte on cmer plced ehind ndpss filter centered on the fundmentl lser frequency. This cn e exploited to determine the plsm reflectivity for the incident lser em. II. EXPERIMENTAL RESULTS The lower prt of Fig.2 summrizes the min finding of the experiment for p polriztion of the incident lser, y presenting the electron em ngulr profiles nd (, E) distriutions, s well s the ngulrly-resolved hrmonic spectr, mesured for two different density grdient scle lengths L, L 1 λ nd L 2 λ: s L is incresed from L 1 to L 2, the mesured signls for ll these oservles rdiclly chnge. Three min differences re oserved: (i) When L λ, the electron emission is predominntly peked t 1 mrd, i.e. close to the direction of lser speculr reflection ( = mrd), with slight shift towrds the trget norml. As L is incresed, it then switches to 2 mrd, direction etween speculr direction nd the tngent to the trget surfce ( = 6 mrd), nd simultneously slightly rodens ngulrly. (ii) Electrons rech energies out twice higher for lrge L (spectrl pek round 1 MeV), with (, E) distriution tht significntly chnges. In the short grdient regime, cler correltion is oserved etween emission ngle nd electron energy E, especilly in the most intense prt of the distriution ( 2 mrd): the electron energy increses s one gets closer to the speculr direction. In contrst, in the long grdient regime, the electron spectrum hrdly vries ngulrly, i.e. no significnt correltion is oserved on this (, E) distriution. (iii) Hrmonic emission is clerly oserved for smll L, ut it drops elow the experimentl detection threshold for lrge L. The detils of the trnsition etween these two regimes re presented in Fig.3, which displys the evolution with L of the electron em ngulr profile in the incidence plne, nd of the hrmonic spectrum. The most importnt point is the quntittive correltion, oserved t short grdients, etween the emission of reltivistic electrons nd the hrmonic signl (see curves in Fig.3c). As L is grdully incresed, the electron signl round = 1 mrd nd the hrmonic signl rech common optimum round L = λ/15, nd then oth quickly decrese. The electron signl on the other side of the speculr direction then grows, ut is not ssocited with ny hrmonic signl. The trnsition etween these two regimes occurs round L λ/5.

5 5 fc/mrd u. 1 c.5 S ROM S e1.4 S e2 L/λ.3 ω (mrd) ω/ω L 3 S (.u) 1 FIG. 3. Evolution of the experimentl oservles with the density grdient scle length. The ngulr profile of the reltivistic electron em in the incidence plne (pnel ), nd the emitted hrmonic spectrum (pnel ) re plotted s function of L, for p-polrized lser field. The experimentl prmeters re the sme s in Fig.2. The insets show how these quntities re relted to the imges of Fig.2 (note tht there is no mesured signl elow these insets). Pnel c shows different curves derived from these dtsets: the hrmonic signl integrted from the 2 th to the 25 th order (signl S ROM ), nd the electron signls on the right (5 mrd 15 mrd, signl S e1) nd left sides ( 25 mrd 15 mrd, signl S e2) of the speculr direction, re plotted s function of L. These curves show the trnsition etween the two regimes highlighted in Fig.2, nd revel the quntittive correltion etween the hrmonic signl nd the reltivistic electrons emission for short grdient scle lengths. (mrd) c p-polriztion (mrd) (mrd) fc/mrd d s-polriztion FIG. 4. Effect of the lser polriztion direction on reltivistic electron emission. Pnel nd show the ngulr profiles of the reltivistic electron emission for grdient scle length L 1 = λ/15, respectively for p nd s polriztions of the incident lser. Pnels c nd d show these electron ngulr distriutions, now mesured for longer grdient L 2 = λ/1.5, i.e. eyond the trnsition oserved in Fig.3. Note the different color scles used in pnels -c versus -d. Another importnt difference etween these two interction regimes is the dependence of the oservles on lser polriztion direction, illustrted in Fig.4. In the short grdient regime, the electron nd hrmonic signls re totlly suppressed when the polriztion is switched from p to s. By contrst, for longer grdients, the electron signl is still oserved for s-polriztion, lthough it gets out five times weker fc/mrd 2 These oservtions on the electron nd hrmonic ems clerly point to complete chnge in the coupling mechnism etween the lser field nd the plsm, which we will nlyze in the rest of this rticle. In the following, we will refer to these two interction regimes s the short-grdient nd long-grdient regimes for convenience. This trnsition lso hs consequences on even simpler oservles: Fig.5 (left pnels) displys the sptil intensity profiles of the reflected lser em, mesured on scttering screen in these two distinct coupling regimes. In the short grdient regime, smooth em is oserved, which is lmost unltered compred to the incident lser em: this is the so-clled plsm mirror regime [19], where the plsm cts s usul high-qulity mirror, speculrly reflecting the fundmentl frequency, despite the ultrhigh intensity on trget. By contrst, in the long grdient regime, the em profile is strongly pertured nd strts exhiiting sptil structures tht were not present on the incident em. The term plsm mirror is thus no longer pproprite to this regime, lthough the lser field still intercts with dense -nd hence reflective- plsm. Experimentlly, the sptil profile of the reflected lser em might thus lso e used s n lterntive nd very simple signture of the trnsition in the lser-plsm interction. By sptilly integrting these imges, the vrition of the plsm reflectivity t the fundmentl lser frequency s function of L cn e determined, nd is displyed in the right pnel of Fig.5, for oth p nd s polriztions of the incident lser field.

6 6 (mrd) (mrd) 1.u Re lectivity c p-polriztion s-polriztion L/λ FIG. 5. Reflected fundmentl em nd evolution of the plsm reflectivity. Using scttering screen, the sptil intensity profile of the lser em reflected y the trget cn e mesured. The two imges on the left show typicl em profiles respectively otined in the short (upper imge, L 1 = λ/15) nd long (lower imge, L 2 = λ/1.5) density grdient regimes, for p-polriztion of the incident lser. The lck dshed circles indicte the initil divergence of the top-ht lser em, efore its interction with the trget. From the sptil integrtion of these imges, the reflectivity of the plsm for the fundmentl lser frequency cn e determined, nd is plotted in pnel c s function of L for oth s nd p polriztions (squres nd circles). The lines show the corresponding results of 2D Prticle-In-Cell simultions (see section IV D). All these mesurements hve een repeted over four different experimentl cmpigns on the UHI1 experimentl fcility over the lst four yers, nd ll effects descried ove hve een oserved to e very reproducile nd roust. III. 3D PARTICLE-IN-CELL SIMULATIONS To interpret these experimentl oservtions, we will now turn to Prticle-In-Cell (PIC) simultions of the lser-plsm interction. This requires ensuring tht (i) these PIC simultions re performed in the ctul physicl conditions of the experiment, (ii) they re relile nd properly reproduce the key experimentl findings. To check these two criticl points, we first crry out full 3D simultions of the interction, so tht we cn directly confront the numericl nd experimentl results -especilly the full ngulr pttern of the electron emission, which is only ccessile y 3D simultions. All simultions reported herein hve een performed using the recently-developed WARP+PXR code [38 42]. The specificity of this code is the use of mssivelyprllel high-order spectrl solver for Mxwell s equtions [43], which gretly reduces numericl dispersion of electromgnetic wves s well s numericl noise. This ensures convergence of the simultions for lrger sptil nd temporl mesh steps thn in most other PIC codes, nd thus mkes physiclly-relistic nd relile 3D simultions of the interction with dense plsms computtionlly trctle [43, 44]. Ech 3D simultions reported here required 6.3 millions computtion hours on mssively-prllel mchine [45, 46]. More detiled informtion on the numericl prmeters of these simultions is provided in the supplementry mteril. We performed two 3D simultions, for the sme physicl conditions s in the experiments of Fig.2 nd c, corresponding to fixed lser prmeters ut different density grdients scle lengths, L 1 λ nd L 2 λ. From these simultions, we extrcted the exct sme oservles s those mesured in the experiment: the clculted ngulr profiles nd ngle-energy distriutions of the emitted electron em, s well s the ngulrly-resolved hrmonic spectr, re displyed in Fig.6. Comprison with the lower pnels of Fig.2 shows tht these simultions very well reproduce the two distinct interction regimes oserved experimentlly for ll these oservles. These 3D enchmrk simultions clerly demonstrte oth the reliility of the PIC simultions, s well s our excellent control of the interction conditions in the experiment, which ensures relevnt choice of the physicl prmeters used in the theoreticl study. In the next section, we will further exploit simultions performed with this WARP+PXR code to get detiled insight into the physicl processes underlying these distinct interction regimes. To this end, we will rely on more trctle 2D simultions [47]. Detiled informtion on the numericl prmeters of these simultions is lso provided in the supplementry mteril. IV. PHYSICAL ANALYSIS AND DISCUSSION The strting point for our nlysis of the experimentl results is the joint mesurement of the reltivistic electron nd hrmonic signls, which provides informtion on the temporl structure of the electron emission y the plsm (section IV A). We then discuss the sptil structure of the oserved electron ems (section IV B). In order to understnd the electron heting nd ejection mechnism in the long grdient regime, we turn to detiled nlysis of PIC simultions, presented in section IV C. With the physicl insight provided y this nlysis, we finlly discuss the influence of the lser polriztion (Fig.4) nd the evolution with L of the plsm reflectivity t the fundmentl lser frequency (Fig.5) in section IV D. A. Temporl structure of the electron emission The cler correltion oserved etween the high-energy electron signl nd the hrmonic signl for short grdients (Fig.3c) suggests tht in this regime, the reltivistic electrons re involved in the hrmonic emission. Then, the fct tht highly-contrsted hrmonic com is produced would e n indiction tht this electron emission is periodic in time, eing locked to the driving lser field.

7 7 (mrd) 2-2 L = λ / 15 fc/mrd (mrd).u ω/ω L 9 (mrd) E (MeV) (mrd) 2 L = λ / 1.5 fc/mrd u θ -1 2 x (mrd) (mrd) (mrd) ω/ω L (mrd) E (MeV) FIG. 6. 3D PIC simultions of the lser-plsm interction for two different density grdient scle lengths. The physicl conditions of these simultions re mtched to the estimted experimentl conditions of the shots shown in the lower prt of Fig.2: pnel corresponds to density grdient L 1 = λ/15 nd pnel to L 2 = λ/1.5, while ll other physicl prmeters remin the sme. From these simultions, we extrct the sme oservles s those mesured in the experiment: the ngulr profile of the high-energy electron em expelled in vcuum (left imge in ech pnel), the (, E) distriution of these electrons (ottom right imge in ech pnel), nd the ngulrly-resolved hrmonic spectrum (top right imge in ech pnel). The key fetures re the sme s those oserved on the experimentl results (compre with the lower pnels of Fig.2). To support this tenttive interprettion, we consider the PIC simultions of Fig.1 nd Fig.7: they show tht temporl periodicity is indeed chrcteristic of the Brunel mechnism. Electron emission occurs in the form of unches tht re initilly extremely short (in the ttosecond rnge), emitted once every opticl period. For > 1, these electrons rech reltivistic velocities when they escpe the plsm, nd thus induce Doppler effect on the reflected lser field: this results in the genertion of trin of ttosecond light pulses, spced y one lser period, which re clerly oserved on Fig.1 nd Fig.7. This well-identified process is known s the Reltivistic Oscillting Mirror (ROM) effect [19, 48 55]. In simultions, the resulting periodic light emission hs spectrum consisting of com of high-order hrmonics of the lser frequency: this is the origin of the hrmonic signl oserved in our experiment. The electron-hrmonic correltions oserved experimentlly is therefore signture of the periodicity of the dynmics of electrons. As descried in section II, the hrmonic signl is oserved to collpse for longer density grdients L. One possile interprettion for this collpse cn e tht the electron emission ceses to e periodic in time. And indeed, PIC simultions for longer density grdients (Fig.7d) strikingly show tht, in contrst to the Brunel mechnism, electron emission is no longer periodic in this regime. The sence of hrmonic signl in conjunction with the reltivistic electron emission cn thus e considered s signture of the trnsition to new coupling mechnism, ssocited to the very different plsm temporl dynmics oserved in Fig.7. This mechnism will e descried in section IV C. B. Sptil structure of the electron emission We now discuss the sptil properties of the outgoing electron ems, with the support of the simultion results of Fig.7. We show tht in the short grdient regime, this structure is mostly determined y the interction of expelled electrons with the reflected lser field in vcuum (section IV B 1), while in the long grdient regime, it is rther imposed y lrge qusi-sttic surfce fields tht develop in the vicinity of the plsm surfce during the interction (section IV B 2). 1. Short-grdient regime In the cse of short density grdients, the peculir ngulr structure of the electron em hs recently een nlyzed experimentlly nd theoreticlly in Ref. [32]. In this Brunel regime, electrons re expelled from the plsm s very lminr em, with reltivistic velocities initilly qusi-prllel to the direction of speculr reflection (Fig.7). These reltivistic electrons thus co-propgte with the intense reflected lser field, with which they interct in vcuum over distnce of the order of the Ryleigh length. This interction lwys results in the ejection of electrons out of the lser em, nd therefore digs hole in the electron em, centered on the speculr direction, s oserved in our experiment. There re two typicl scenrios for this ejection, depending on the exct initil conditions with which electrons re expelled from the plsm into vcuum. Some electrons explore multiple opticl cycles of the lser field, nd thus oscillte in the field nd get expelled from the

8 8 c 1 2 By x/λ z/λ d z/λ x/λ t/t 25 FIG. 7. 2D PIC simultions in the two distinct regimes of lser-plsm coupling. These dt re otined from 2D PIC simultions with different density grdients L (λ/15 for nd, nd λ/1.5 for c nd d), while ll other physicl prmeters remin the sme ( = 3.5, θ i = 55 o ). The two upper pnels disply the complete trjectories of selected set of expelled high-energy test electrons (ornge lines), together with the totl y-component of the totl B-field (lue to red color mp) t given time fter the lser-plsm interction. The plsm density profile t the end of the interction is indicted in gry in log scle. The lower pnels shows the temporl evolution of the plsm electron density (gry-scle color mp, in log scle), sptilly-resolved long the norml to the trget surfce, t the center of the focl spot. The emitted ttosecond pulses re superimposed to this density mp in purple. They re clerly visile in pnel, ut re too wek to e oserved in pnel d t/t 3 lser focl volume y the so-clled ponderomotive effect, isotropiclly nd with limited energy gin [56, 57]. They form the ring-shped hlo oserved on the electron em. But most electrons ctully remin round given phse of the reflected field nd rther surf single wvefront of the reflected field, thus escping the lser em lterlly long the lser polriztion direction, nd forming the right pek oserved next to the speculr direction. The side on which this pek forms is determined y the lser phse t which electrons re expelled from the plsm into vcuum: the oservtion of pek on one side only of the ponderomotive hole (etween the speculr direction nd the trget norml) is cler indiction tht electrons re ejected periodiclly once every lser period, when the lser field drgs them out of the plsm, in the form of su-opticl cycle unches. This second set of electrons experiences qusiconstnt E-field from the lser in vcuum until they escpe the focl volume, leding to greter energy gin thn in ponderomotive scttering: this vcuum lser ccelertion (VLA) process ccounts for the oserved symmetry of the (, E) distriution (see Fig.2), where higher energies re oserved on one side of the hole (VLA electrons) thn on the other ( ponderomotive electrons), s well s for the ngle-energy correltions on this distriution [58]. An importnt consequence is tht the Doppler upshift fctor induced y outgoing electrons on the reflected field, which leds to the genertion of highorder hrmonic (ROM mechnism), cnnot e directly deduced from the electron spectr mesured experimentlly, since electrons keep gining energy fter they escped the trget nd emitted high-order hrmonics, nd efore eing detected. For instnce, simultions show tht in the present experiment, the electron Lorentz fctor typiclly vries from γ 2 3 s they re ejected from the plsm nd emit hrmonics, to γ 15 fter their interction with the reflected lser field [32], when they re detected. 2. Long-grdient regime The sptil properties of the electron em oserved in the long-grdient regime descried here hve not een explined in detil yet, to the est of our knowledge. The electron trjectories displyed in Fig.7 show tht the conditions of electron ejection from the plsm re lredy very different in the short nd long grdient regimes. In the second cse, the expelled electron em is no longer lminr, nd rther hs complex velocity distriution. As will e shown in the next section, this feture cn e ttriuted to the chotic chrcter of the electron heting mechnism leding to ejection from the plsm. Furthermore, these 2D PIC simultions show tht qusi-sttic mgnetic field develops t the plsm surfce (see mp of the mgnetic field in Fig.7c, where red re

9 9 indicting high mgnetic field is present close to the plsm surfce). This field grows during the lser-plsm interction, reches n mplitude typiclly of the sme order of mgnitude s the lser mgnetic field, nd then persists even fter the lser pulse hs een reflected y the plsm. A detiled nlysis of electron trjectories in these simultions shows tht this surfce field, which is much lrger thn the one occurring in the short grdient regime (compre mps of the mgnetic field in Fig.7 nd c), strongly deflects the escping electrons towrd the trget surfce. This deflection ccounts for the fct tht the electron ngulr distriutions oserved in this regime re essentilly centered etween the speculr direction nd trget surfce. Test simultions hve een performed to check tht the reflected lser field plys no role in the deflection of electrons fter their ejection from the plsm, in strong contrst to wht is oserved in Brunel regime. Such surfce qusi-sttic fields hve lredy een reported in multiple studies of the interction of intense lsers with dense plsms (see e.g. [59 62]), nd cn e induced y vriety of physicl processes (see e.g. [63 66]). In the present cse, our simultions indicte tht their development cn e ttriuted to the fountin effect descried in [66] where they originte from the plsm cold return current tht compenstes for the lterl chrge ejection from the lser focl volume. This is supported y the fct tht these fields re not oserved t ll in plne wve simultions, where the plsm surfce is homogeneously illuminted y the lser field. In contrst to the short-grdient regime, here the sptil properties of the electron em do not provide much insight into the involved electron heting mechnism. To identify this mechnism, we now turn to more detiled numericl investigtion sed on PIC simultions. C. Electron heting mechnism in the long grdient regime 1. Importnce of the reflected lser field Our nlysis of the electron heting mechnism in the long grdient regime is sed on set of 2D plne wve PIC simultions [47] for three different physicl configurtions. Their key results re summrized in Fig.8, nd shed light on the underlying physicl mechnism. The upper pnels (cse A) correspond to the interction of n ultrintense lser pulse with dense plsm in the long grdient regime, i.e. the sme physicl configurtion s in Fig.7 nd s in our experiment: they disply the temporl evolution of the plsm electron density (pnel ), nd the x p x phse spce distriution of electrons (pnel d) t the time when electron ejection from the plsm is oserved to strt (lue dshed line in pnel ). The middle pnels (cse B) disply the sme quntities, ut now in sitution where the plsm profile hs een truncted for densities n >.4n c cos 2 θ, i.e. keeping only the underdense prt of the plsm, such tht there is hrdly ny reflection of the incident lser y the plsm. This underdense plsm lyer, surrounded y vcuum on oth sides, is irrdited y the sme lser em s efore, ut lso y second em of slightly lower intensity (8%), symmetriclly rriving from the other side of the plsm. The role of this second lser is to emulte the em reflected y the dense prt of the plsm in cse A. The key point here is tht oth the temporl dynmics of the plsm density profile nd the electron phse-spce distriutions look very similr for cses A nd B. By contrst, if the truncted plsm lyer is irrdited y one lser em only (cse C, lower pnels), the plsm dynmics ecomes totlly different. More specificlly, while similr upwrd electron ejections re oserved in cses A (corresponding to the electron signl oserved in our experiment) nd B, this electron emission is strongly reduced in cse C: less electrons re emitted, nd they hve much weker velocities. This toy-model study leds to two importnt conclusions. (i) The comprison of cses A nd B indictes tht in the long grdient regime, the coupling mechnism leding to electron ejection mostly occurs in the underdense prt of the density grdient. (ii) The comprison of cses B nd C indictes tht the overdense prt of the plsm nonetheless plys key role, y producing reflected em tht, when crossing nd interfering with the incident em, strongly modifies the dynmics of electrons in the underdense plsm lyer. As we explin in the next section, the electron heting mechnism coming into ply in such physicl sitution is lredy well-identified in the existing literture. 2. Introduction to stochstic heting Let us strt y considering free electron (i.e. without ny collective plsm effects involved) exposed to two non-colliner ultrintense lser ems. Using quntum description of the field s n ensemle of photons provides simple wy to understnd tht this electron cn gin more energy thn when exposed to single lser em. When single lser em (ssumed to e plne wve) is present, it is well-known tht photon sorption processes re hindered ecuse they do not enle to conserve oth energy nd momentum of the totl system. By contrst, when two non-colliner ems re present, the comined sorption of multiple photons simultneously from oth ems is llowed, ecuse the vilility of photons with different k vectors mkes it possile to conserve oth energy nd momentum of the totl system. In other words, the presence of second em llows for energy sorption y the electron from the lser field, through process tht cn e defined s stimulted multiphoton Compton scttering [67]. For lrge field mplitudes, the lser field cn e treted clssiclly, nd mny previous studies in the literture hve shown tht electron dynmics in these comined

10 1 x/λ x/λ log(n/n c ) -2 2 nc cos 2 θ α nc cos 2 θ α nc cos 2 θ c d e f terministic nd involves no stochstic process. This effect is oviously not restricted to isolted free electrons: it cn eqully occur for electrons in n underdense plsms, leding to n energy sorption process where neither collisions nor collective plsm effects ply ny mjor role. Such coupling of the plsm with multiple lser ems hs een studied experimentlly in [72] y exposing n underdense plsm to two lser ems (like in cse B of Fig.8). It is lso known to ply role in electron injection in lser-driven plsm wkefield ccelertors y the colliding pulse scheme [73, 74]. To the est of our knowledge, Ref. [34] ws the first to point out tht this mechnism should lso come into ply when n smll underdense plsm lyer is present t the surfce of dense plsm exposed to single lser em. In this cse, the required second non-colliner lser em results from the reflection of the single input em y the dense plsm. Electrons in the underdense plsm re then exposed to the stnding wve formed in front of the dense plsm y the superposition of the incident nd reflected ems, nd cn gin energy y stochstic heting: this is precisely how we interpret our present experimentl results in the long grdient regime. x/λ 2 3. Numericl evidence for stochstic heting in the long grdient regime t/t p x (m e c) FIG. 8. Set of 2D plne-wve PIC simultions crried out to revel the role of the lser field reflected y the plsm. Ech line corresponds to different physicl cse (see insets in pnels -c). cse A (upper line) corresponds to n overdense plsm with density grdient of scle length L = λ, irrdited y single lser em with n incidence ngle θ i = 55 o nd = 2.5. In cse B, this plsm hs een truncted for n >.4n c cos 2 θ, keeping only the underdense prt of the plsm, which is now irrdited y two lmost identicl lser ems (sme prmeters s in cse A) rriving symmetriclly from oth sides of the plsm lyer. In cse C, this sme plsm lyer is irrdited only y the upper lser em. In these three cses, plots of the temporl evolution of the plsm electron density (pnels -c) nd snpshots of the x p x phse spce distriution of electrons (pnels d-f, tken t the time indicted y the lue line in pnels -c) re displyed. The multilyered phse spce distriutions of pnels d nd e strongly contrst with the smooth regulr distriution of pnel f, nd re typicl of chotic dynmics resulting from repetitive stretching nd folding effect in phse spce. noncolliner fields is not integrle nd gets chotic for high enough lser mplitudes (typiclly.15 for t lest one of the two ems). This results in lrge energy gins, nd this efficient regime of energy sorption y electrons is known s stochstic heting [35, 68 71] -lthough the nme chotic heting would proly e more pproprite here, since the system is perfectly de- With the help of simultions, we now support this interprettion y providing evidence tht stochstic heting indeed occurs in the long grdient regime. To this end, we nlyze the electron phse spce distriutions in cses A, B nd C (pnels d-f in Fig.8) nd their temporl evolution during the lser-plsm interction (see movie M1). When single lser em is present (cse C), electrons re oserved to simply oscillte non-linerly in the lser field, leding to smooth nd regulr phse spce distriution (pnel f nd movie M1). In striking contrst, in cses A nd B, electron dynmics in the stnding wve resulting from the superposition of two non-colliner lser ems is complex: the key point is tht we oserve very strong locl stretching nd folding effect on the phse-spce distriution, round ech node of the stnding wve (see movie M1). Such stretching nd folding effect results in very different trjectories for prticles tht re initilly very close in phse spce: this is known to e one of the most typicl routes to chotic dynmic [75], exemplified in the wellknown horseshoe mp models. These repetitive stretching nd folding eventully result in highly-structured, multilyered phse spce distriution (pnels d-e), where electrons t given sptil position hve complex momentum distriution, typicl of chotic dynmics. The striking contrst etween these highly-structured phse spce distriutions, nd the smooth distriution oserved in cse C (Fig.8f) gin demonstrtes the impct of stochstic heting on electron dynmics in the underdense prt of the plsm. Furthermore, the comprison

11 11 p-polriztion s-polriztion σt σ px Short grdient Long grdient Fundmentl Hrmonics Sttic ields Non-speculr EM sources Prticle kinetic energy FIG. 9. Lypunov exponent for the p x vrile in the long grdient regime, s function of. Simultions hve een performed in the physicl configurtion of cse A of Fig 8. The onset of stochstic heting in the stnding wve formed in the underdense prt of the plsm y the incident nd reflected lser fields is indicted y the fct tht the Lypunov exponent ecomes positive. The clcultion procedure used to deduce this exponent from the results of 2D plne-wve PIC simultions is explined in the supplementry mteril. of these phse spce distriutions with the type of distriution oserved for the Brunel mechnism (Fig.1) is illustrtive of the mjor difference in the dynmics of the system etween the short nd long grdient regimes. The chotic chrcter of the electron dynmics cn e further supported y the clcultion of the Lypunov exponents of plsm electrons (see supplementry mteril), which should e positive in the cse of chotic dynmics. The Lypunov exponent σ px for the p x vrile, otined from 2D plne wve PIC simultions of cse A of Fig.8, re displyed in Fig.9, s function of the incident lser mplitude. This exponent is negtive t low intensity, nd gets positive when.15, thus pointing to chotic dynmics. This threshold in lser intensity is consistent with erly theoreticl investigtions of stochstic heting [35, 68]. D. Effect of the polriztion, plsm reflectivity With the support of the previous physicl nlysis, we finlly discuss the experimentl oservtions on the influence of the lser polriztion direction (Fig.4), nd the evolution of the plsm reflectivity with L, for oth s nd p polriztions (Fig.5). In the short grdient regime, switching the polriztion from p to s is expected to suppress oth electron emission nd hrmonic genertion, s oserved experimentlly (Fig.4-), ecuse the lser E-field component norml to the plsm surfce is the min driving force of Brunel sorption t the lser intensities considered here. In contrst, when stochstic heting is involved, the polriztion direction of the incident lser em is not expected to hve strong influence, since the plsm FIG. 1. Distriution of the initil lser pulse energy fter the interction, in different conditions. These pie chrts summrize how the lser energy is distriuted fter the lser plsm interction, in the short (upper line) nd long (lower line) grdient regimes, for p (left column) nd s (right column) polriztions. Five ctegories hve een numericlly seprted, indicted in the legend: electromgnetic energy in the fundmentl lser frequency in the speculr direction, in the hrmonics of the lser frequency in the speculr direction, in qusi-sttic fields (i.e. with frequency lower thn the lser frequency), in non-sttic fields in the non-speculr direction, nd kinetic energy of plsm prticles. The results re otined from 2D PIC simultions with = 3, θ i = 55 o, nd L 1 = λ/15 for the short grdient regime, L 2 = λ/1.5 in the long grdient regime. surfce only comes into ply y producing reflected wve: this explins why, in the long grdient regime, the electron signl is experimentlly oserved to still persist in s-polriztion (Fig.4c-d). This experimentl oservtion lone mkes strong cse ginst n interprettion of the long grdient regime in terms of resonnt sorption, which should rther e very sensitive to the lser polriztion. In the long grdient regime, sptil degrdtion of the reflected lser em wvefronts is oserved right fter the trget in these simultions (Fig.7c), while the em wvefront is preserved in the short grdient regime (Fig.7). This is qulittively consistent with experimentl oservtions (Fig.5), where the lser em intensity profile fr from the trget is oserved to ecome degrded for long grdients. The chotic chrcter of the electron dynmics, identified in the previous section, ffects the lser em propgtion in the underdense plsm lyer, nd thus provides possile interprettion for this degrdtion of the reflected lser wvefronts. From more quntittive point of view, 2D PIC simultions lso properly reproduce the evolution with L of the plsm reflectivity t the fundmentl lser frequency, for oth s nd p polriztions (Fig.5). For short density grdients, the reflectivity is lower in p ( 5%) thn in s (close to 1%) polriztion, while similr vlues ( 7%) re oserved for oth polriztions in the long grdient regime. This is consistent with the results of Fig. 4 on reltivistic electron emission, nd with our previous interprettion of the inter-

12 12 ction: in s-polriztion, for short L, Brunel sorption is suppressed, while for long L stochstic heting is still present. These 2D PIC simultions cn e exploited to get more insight into the redistriution of the initil lser energy fter the lser-plsm interction. The different distriutions otined for oth short nd long grdient regimes, nd for oth p nd s lser polriztions, re displyed s pie chrts in Fig. 1. For short grdients nd p- polriztion, round 25 % of the lser is converted into hrmonics of the lser frequency (mostly low-orders), nd out one third is deposited s kinetic energy of the plsm prticles (mong which the reltivistic electrons oserved in our experiment). When the polriztion is switched to s, these two contriutions get considerly reduced, down to round 5% ech, leding to much higher reflectivity of the fundmentl lser frequency. For long grdients, y contrst, very little energy is converted into hrmonics, regrdless of polriztion. As lredy emphsized, the energy stored into qusi-sttic fields round the plsm surfce significntly increses compred to the short grdient regime. In p-polriztion, the frction of energy going into prticle kinetic energy is only slightly weker thn in the short grdient cse, nd gets reduced y out 5% in s-polriztion. The reflectivity for the fundmentl frequency remins similr for oth polriztions, while more energy goes into scttered light in s-polriztion. E. Conclusion of the physicl nlysis The previous comintion of multiple experimentl oservles nd PIC simultions hs provided strong evidence for the trnsition from Brunel sorption to stochstic heting s the density grdient scle length L is incresed, while no evidence of resonnce sorption hs een oserved in this ultrhigh lser intensity regime. The intuitive physicl insight underlying this trnsition is tht the Brunel mechnism requires shrp interfce, such tht the mplitude of the quivering motion of electrons in the lser field exceeds the length scle of this interfce. On the opposite, stochstic heting cn only occur if electrons re present within the stnding wve interference pttern formed y the incident nd reflected fields: it will therefore e fvored y longer density grdients, for which this interference zone contins lrger numer of prticles provided y the underdense prt of the plsm locted just in front of the lser reflection point. It however remins difficult to predict nlyticlly the vlue L t of the density grdient scle length for which this trnsition occurs. In the next section, we determine this trnsition length experimentlly s function of the lser incidence ngle on trget θ i - n essentil physicl prmeter, since it ffects the point of the density grdient where lser reflection occurs. We lso discuss the effect of this ngle on the properties of the emitted electron ems. V. EFFECT OF THE LASER INCIDENCE ANGLE A. Evolution of the trnsition grdient scle length We hve repeted the previous mesurements, crried out for θ i = 55 o, for five other incidence ngles θ i rnging from 4 o to 65 o. For ech ngle, the density grdient scle length ws systemticlly vried y chnging the prepulse dely, nd ws mesured using the SDI technique. The min outcomes of this experiment re summrized in Fig.11. The left pnels show how the electron em ngulr profile in the incidence plne evolves s function of L, for three different incidence ngles. In ll cses, the sme type of trnsition s reported in Fig.2 nd 3 is oserved, nd it occurs for shorter vlues of L s the incidence ngle is incresed. As efore, we oserve tht this trnsition is correlted with mjor chnges in the electron ngle-energy distriution nd in the hrmonic emission. The experimentl trnsition length L e t deduced from ll these dt is plotted in Fig.11 (lue dots), nd clerly decreses with θ i. To get some qulittive insight into this ngulr dependnce, we consider the strting point of Brunel s model [11], which is tht the physics of the lser-plsm coupling chnges when the quivering mplitude z e of plsm electrons strts exceeding the typicl sptil extent d of the plsm-vcuum interfce. We then define - somewht ritrrily- this sptil extent s d = z s z c, the distnce etween the effective reflective surfce of the plsm, locted t z s (θ i ) such tht n(z s ) = n c cos 2 θ i, nd the fixed point corresponding to the loction of the criticl plsm density, z c such tht n(z c ) = n c (see sketch in Fig.11). d cn e clculted y using the stndrd ssumption of n exponentil density profile t the trget surfce [76]: n(z) = n exp( z/l), with n the mximum plsm density in the trget, reched in z =. This leds to d = 2L ln(cos θ i ). In the simplest possile pproch, the electron quivering mplitude cn e considered to e the one of free electrons in the lser field, which is z e λ/2π 2 /( 2 +1). At sufficiently high intensity, well into the reltivistic regime ( 1), z e λ/2π is qusi-independent of the lser field mplitude. The condition z e = d(θ i ) then leds to simple theoreticl expression for the trnsition grdient scle length (see movie M1): L t λ = 1 4π ln(cos θ i ) (1) This simple nlysis indeed predicts n ngulr dependence of L t on θ i : physiclly, this is ecuse for given plsm density profile t the surfce (i.e. fixed L), the distnce d increses with θ i, due to the wellknown ngulr dependence of the effective criticl den-