arxiv: v1 [cond-mat.mtrl-sci] 16 Dec 2011

Transcription

1 Coherent Phonon Coupling to Individul loch Sttes in Photoexcited ismuth E. Pplzrou, J. Fure,3, J. Muchin, M. Mrsi,. Tle-Irhimi, I. Reshetnyk,. vn Roekeghem, I. Timrov, N. Vst,. rnud 5, L. Perfetti Lortoire de Physique des Solides, CNRS-UMR 85, Université Pris-Sud, F-95 Orsy, Frnce Lortoire des Solides Irrdiés, Ecole polytechnique-ce/dsm-cnrs UMR 76, F-98 Pliseu, Frnce 3 Lortoire d Optique ppliquée, Ecole polytechnique-enst-cnrs, F-98 Pliseu cedex, Frnce Synchrotron SOLEIL, Sint-uin-P 8, F-99 Gif sur Yvette, Frnce nd 5 Institut de Physique de Rennes (IPR), UMR UR-CNRS 65, F-35 Rennes Cedex, Frnce rxiv:.399v [cond-mt.mtrl-sci] 6 Dec We investigte the temporl evolution of the electronic sttes t the ismuth () surfce y mens of time nd ngle resolved photoelectron spectroscopy. The inding energy of ulk-like nds oscilltes with the frequency of the g phonon mode wheres surfce sttes re insensitive to the coherent displcement of the lttice. strong dependence of the oscilltion mplitude on the electronic wvevector is correctly reproduced y initio clcultions of electron-phonon coupling. esides these oscilltions, ll the electronic sttes lso disply photoinduced shift towrds higher inding energy whose dynmics follows the evolution of the electronic temperture. PCS numers: It is well known tht most elements crystllize in structures with hexgonl or cuic symmetry. One notle exception is ismuth, which insted crystllizes in the 7 rhomohedrl structure. It hs een originlly proposed y Jones nd Peierls tht cuic ismuth would e unstle towrds lttice distortion leding to the douling of the unit cell []. Due to the distorted nture of the groundstte, the tomic positions re very sensitive to the electronic distriution in the conduction nd. This conjecture hs een experimentlly verified y suddenly chnging the occuption numer of the electrons vi the sorption of femtosecond lser pulse [ ]. Due to photoexcittion, lrge coherent phonon modultes the distnce etween the two toms in the unit cell long the () direction. The displcement vector indictes tht the lttice coordintes move towrds hypotheticl phse with single tom per unit cell. In relity the mteril displys non-therml melting efore the symmetric phse cn e ttined [5]. The tomic motion following the photoexcittion of ismuth is nowdys well understood. Time resolved X- ry diffrction nd Density Functionl Theory (DFT) clcultions could ccurtely descrie the mplitude of the oscilltions nd the frequency softening t high excittion densities [6 8]. On the other hnd, the temporl evolution of electronic sttes hs never een directly oserved. Some spects of the electronic dynmics cn e inferred from generl rguments. The genertion mechnism of coherent phonons with lrge mplitude implies tht low energy electronic excittions re coupled to the g mode. It is lso expected tht the electron-phonon coupling of such distorted mteril vries with the electronic wvevector. Nonetheless, the determintion of the coupling mtrix elements for individul loch sttes is still n experimentl chllenge. We pproch this tsk y time resolved photoelectron spectroscopy mesurements. The electronic sttes of the ismuth () surfce re chrcterized y ngle resolved photoemission nd first principles clcultions. We confirm tht the electronic structure displys rich comintion of ulk-like nds, surfce sttes nd surfce resonnces [9, ]. Upon photoexcittion, these electronic sttes disply dynmics tht depends on wvevector nd nd index. The connection of our dt with DFT clcultions provides full chrcteriztion of the coupling etween electrons nd the g mode. In ddition, we oserve tht the electronic sttes do not follow rigidly the motion of the lttice displcement. Insted, the nds lso disply purely electronic shift towrds higher inding energy. The experiments were performed with the FemtoRPES setup, using Ti:Spphire lser tht genertes 35 fs pulses centered t 8 nm with repetition rte of 5 khz. Prt of the em is employed to generte the fourth hrmonic y frequency douling in O crystls (β- O ). The 5 nm proe nd the 8 nm pump re focused on the smple with spot dimeter of µm nd µm, respectively. Their cross-correltion in O crystl hs full width t hlf mximum (FWHM) of 8 fs. n electrosttic spectrometer with energy resolution etter thn mev nd ngulr resolution etter thn.5 degrees nlyses the photoelectrons emitted y the 5 nm em. The () surfce of ismuth hs een otined y sputtering nd nneling cycles of single crystl. ll mesurements hve een performed t the se temperture of 3 K nd t the se pressure of 7 mr. The out-of-equilirium spectr hve een collected with incident pump fluence of.6 mj/cm. Figure shows n intensity mp of the photoelectrons emitted long the Γ-M direction of the rillouin zone. s previously reported, the rek-down of trnsltionl symmetry in the () direction genertes surfce sttes tht intersect the Fermi level nd give rise to Fermi surfce [9]. These evnescent wvefunctions re loclized t the topmost lyers, thus conferring to the surfce of ismuth good metllic properties. lthough, the nd structure supports surfce sttes only ner the Fermi level, some ulk-like sttes nd surfce resonnces re clerly visile t higher inding energy []. Two

2 M.. M' FIG. : ): Intensity mp of photoelectrons emitted for prllel wvevector long the Γ-M direction. nd nd nd re the surfce resonnce nd the ulk-like stte, respectively. The lue line is ulk nd clculted for perpendiculr wvevector k =.5 Å. ): Intensity mp of photoelectrons emitted t inding energy ɛ =.6 ev. n intense structure originting from nd is visile long the Γ-M direction ut not long Γ-M. nds of different nture re indicted in Fig. y letter nd. Due to the mtrix elements of the photoemission process, nd is clerly visile only when the prllel component of the electronic wvevector is k <.5 Å. nd hs lrger inding energy, disperses more strongly nd intersects nd t the Γ point. The overll greement etween our mesurements nd one step photoemission clcultions is remrkle [].. Kimur et l. [] lso show tht nd hs lrge spin polriztion wheres nd does not. This finding suggests tht is surfce resonnce wheres is ulk nd. In order to confirm it, we performed DFT clcultions of ulk ismuth within the frmework of the locl density pproximtion nd the INIT [] pckge. Since the photoemission process does not conserve the perpendiculr component of the electronic wvevector, the connection etween pek in the photocurrent nd clculted nd is formlly not llowed [3]. Nonetheless, Fig. shows tht the nd clculted long the Γ-M direction for perpendiculr wvevector k =.5 Å mtches resonly well the experimentl intensity distriution.. M Intensity Temperture (K) fs 3 fs C Low 3 High FIG. : ): Photoelectron spectr cquired fs efore nd 3 fs fter the rrivl of the pump em long the Γ-M direction t k =. Å. ): The effective temperture of photoexcited electrons s function of pump proe dely extrcted from the dt (mrks) nd simulted with the model of ref. [8] (red line). C): Intensity mp of photoelectrons emitted t wvevector k =. Å long the Γ-M direction s function of pump proe dely. More insights on the dispersion of nd cn e otined y mpping in the reciprocl spce the photoelectrons t fixed kinetic energy. Figure shows n intensity mp of photoelectrons cquired t different wvevectors for inding energy ɛ = 6 mev. The structure is very intense long Γ-M wheres it vnishes in the direction Γ-M. The experimentl evidence tht nd is invrint for rottion of ut not for rottion of 6 corroortes the ulk chrcter of nd [9]. Photoexcittion y n intense nd ultrfst lser pulse genertes sizle effects on the photoelectron current (see movie in the supplementry mteril). Figure shows the spectrum cquired t k =. Å long the Γ- M direction for two different pump-proe delys. Note the lrge reduction of spectrl weight tking plce 3 fs fter the rrivl of the pump em. The prtil depletion of sttes nd cn e pproximtely descried y Fermi-Dirc distriution with n effective temperture of 9 K. s shown y Fig., the electronic temperture decys s function of pump proe dely with time constnt of 6 ps. The relxtion of the electronic energy density t the surfce is due to the comined effect of electronic diffusion nd het dissiption into the lttice th. We could ccurtely reproduce the temporl evolution of the electronic temperture y solving the coupled differentil equtions in ref. [8] with incident pump fluence of.6 mj/cm, electron therml conductivity k =. W m K nd temperture dependent

3 3 electron-phonon coupling. eside the lrge chnges of the electronic distriution, we lso oserve strong vritions in the structure of the electronic sttes. s shown y Figure, oth pek nd shift towrds higher inding energies fter the rrivl of the pump pulse. In ddition, the position of pek displys periodic modultions t the frequency of the g mode. Such dynmics rises from the coupling of the coherent phonon to electronic sttes nd cn e directly oserved on the photoelectron intensity mp of Fig. C. The oscilltions re lrge for the ulk nd wheres they fll elow the detection limit for the surfce resonnce nd for ll other surfce sttes. This oservtion leds to two possile conclusions: either the mplitude of the g mode is smller on the topmost ilyer or the surfce sttes hve wek coupling to this mode. The excittion of coherent phonons in semimetls hs een the suject of mny experimentl nd theoreticl works [ 8]. Since ismuth sors photons t 8 nm, the genertion mechnism is minly due to displcive excittion of the chrge density [, 3]. The DFT clcultions cn help us to illustrte this point: figure 3 displys colorscle plot of the electronic density ρ(r, T i ) t T i = 3 K in the ( ) plne. The sudden increse of electronic temperture from 3 K to T f = K chnges the sptil ditriution of electronic density. We show in Fig. 3 tht (ρ(r, T f ) ρ(r, T i ))/ρ(r, T i ) is negtive in the region etween nerest neighours (toms nd ). The trnsfer of electronic chrge out of the onding region increses the repulsion etween the ion cores. s consequence, the toms oscillte round new equilirium position until the energy dissiption into other phonon modes rings the system ck to the initil stte. The lttice motion chnges the inding energy of the loch stte, k y D,k u, where u is the g phonon coordinte nd D,k is the deformtion potentil. Figure 3C shows the inding energy of pek t the se temperture of 3 K s function of pump proe dely. The oscilltions induced y the g mode hve frequency ν =.97 ±.5 THz nd dmping time γ =.6 ±. ps. These vlues compres well with the frequency nd dmping time oserved y trnsient reflectivity in comprle experimentl conditions []. We reproduce the oscilltions of pek y solving the model proposed in Ref. [8] nd using the vlue of the computed deformtion potentil D,k =.9 ev/å for k =. Å. The simulted lttice dynmics ttins mximum displcement vlue u =.6 pm, corresponding to. % of the ody digonl length c =.86 Å. s shown in Fig. 3D, the mesured mplitude of oscilltion is mev for k <.5 Å nd decreses down to mev for k <. Å. t lrger wvevectors the photoelectron intensity of nd is too low for relile extrction of the D,k u vlue. The lrge k -dependence of D,k hs een previously reported for the mplitude mode of TTe 3 [5]. Due to the presence of Chrge Density Wve, the vrition of the interction strength Oscilltion mplitude (mev) C D x %.%.%.3% FIG. 3: ): Colorscle plot of the electronic density t 3 K in the ( ) plne. ): Reltive chnge of the electronic density for n increse of electronic temperture from 3 K to K. C): inding energy of pek for k =. Å long the Γ-M direction s function of pump proe dely (lue mrks) nd temporl evolution otined y initio clcultions (red line) [8]. D): Oscilltion mplitude of pek s function of k (mrks) together with the prediction of DFT theory for lttice displcement u =.6 pm (solid line). E): Clculted energy dispersion of nd long the Γ-M direction when the g phonon mode is ctivted (u ). The clcultions re performed for two toms locted t ±(.39 c + u) long the digonl of the unit cell for u/c =. %,. % nd.3 %. s function of the loch stte hs een explined y the good nesting properties of the Fermi surfce. On the other hnd, our mesurements suggest tht the wvevector dependence of the electron-phonon coupling constitutes generl property of covlent crystls. This spect of the coupling cn e understood y simple tight inding model on liner chin: if the lttice displcement modifies the overlp etween tomic oritls the ind- E

4 Shift nd Electronic T 3 FIG. : Pek position of the surfce resonnce for prllel wvevector k =. Å s function of pump proe dely (red mrks, left xis), nd electronic temperture (lck line, right xis). ing energy of the loch sttes will chnge differently for different wvevectors. Such k dependence of the electronphonon coupling is essentil for severl prolems of solid stte physics: it ccounts for the presence of Chrge Density Wves tht do not respect the nesting condition [6] nd it llows for the correct estimte of the high trnsition temperture in superconducting Mg [7]. In the cse of ismuth, our dt offer n unique opportunity to test the deformtion potentil extrcted y DFT. We clculte the nd structure of ismuth for severl vlues of the displcement u. Figure 3E shows the dispersion of nd for u etween.% nd.3% of the ody digonl length c =.86 Å. In grement with the mesurements, the inding energy of nd strongly vries ner the Γ point wheres is less sensitive to the phonon coordinte for k =. Å. On the other hnd, Figure 3D shows tht the wvevector dependence of D,k is stronger in the clcultions thn in the experiment. more ccurte simultion of the mesured signl would require the computtion of one-step photoemission spectr []. It is cler in Fig. 3C tht the dynmics of electronic sttes does not depend simply on the lttice coordintes. esides coherent phonon oscilltions, n dditionl effect shifts pek towrds negtive vlues. This photoinduced shift occurs lso in the surfce resonnce nd Temperture (K) surfce sttes. s shown in Fig., the dynmics of pek nerly follows the electronic temperture. Under similr experimentl conditions, the photoinduced increse of inding energy s een oserved lso in the surfce sttes of Gdolinium [8]. n opposite ehviour hs insted een reported in the electronic sttes of mterils with photoinduced phse trnsitions [, 5]. The detiled nlysis of the electronic shift s function of k will e prt of fore coming work. Nonetheless, we notice tht the mgnitude of this component depends on the electronic wvevector nd nd index. wvevector dependent stiffening of the electronic sttes in the onding nd cn e expected ecuse of the reduced ion core screening t elevted electronic temperture. However, the finite temperture DFT clcultions indicte tht the chrge trnsfer of figure 3 ccounts t mximum for mev to the shift oserved in nd. This discrepncy my e due to the fct tht our ulk model of the electronic sttes does not tke into ccount chrge redistriution nd crrier trnsport t the surfce of the smple. lterntively, the dynmicl spect of the electronic screening my require simultions tht go ehind the DFT method [9]. In conclusion, the electronic sttes of photoexcited ismuth oscillte with the frequency of the g mode. y nlyzing the photoelectron spectr collected t different ngles, we reconstruct the deformtion potentil of specific loch sttes. quntittive greement with DFT clcultions confirms tht the photoinduced chrge trnsfer ccounts for the mplitude of the oserved oscilltions. t lst, we mention tht purely electronic effect dds to the temporl evolution of the electronic sttes with the dynmics of the electronic temperture. The FemtoRPES project ws finncilly supported y the RTR Tringle de l Physique, the NR progrm Chires d Excellence (Nr. NR-8-CEXCEC8- -). We lso knowledge the contriutions of DG nd PNNO CCTTONE. Prt of the Density Functionl Theory results hve een otined or reproduced lso y the Quntum Espresso [] pckge. Clcultions were performed using HPC resources from GENCI- CINES (Nr nd project ). [] R. Peierls, More surprises in theoreticl physics (Princeton University Press, New Jersey, 99). [] H. J. Zeiger et l., Phys. Rev. 5, 768 (99). [3] M.F. DeCmp, D.. Reis, P. H. ucksum, nd R. Merlin, Phys. Rev. 6, 93 (). [] M. Hse, K. Mizoguchi, H. Hrim, S. I. Nkshim, nd K. Ski, Phys. Rev. 58, 58 (998). [5] G. Sciini et l., Nture 58, 56 (9). [6] D. M. Fritz et l., Science 35, 633 (7). [7] S. L. Johnson et l., Phys. Rev. Lett., 555 (8). [8] Y. Giret, Gellé, nd rice rnud, Phys. Rev. Lett. 6, 5553 (). [9] Yu. M. Koroteev et l., Phys. Rev. Lett. 93, 63 (). [] C. R. st nd H. Höchst, Phys. Rev. 7, 5 (). []. Kimur et l., Phys. Rev. Lett. 5, 768 (). [] X. Gonze et l., Computer Phys. Commun. 8, 58 (9) [3] E. E. Krsovskii, K. Rossngel,. Fedorov, W. Schttke, L. Kipp, Phys. Rev. Lett. 98, 76 (7). [] L. Perfetti et l., Phys. Rev. Lett. 97, 67 (6). [5] F. Schmitt et l., Science 3, 69 (8). [6] my Y. Liu, Phys. Rev. 79, 55 (9). [7] H. J. Choi,D. Roundy, H. Sun, M. L. Cohen, nd S. G. Louie, Phys. Rev. 66, 53 (). [8] U. ovensiepen, J. Phys.: Condens. Mtter 9, 83

5 5 (7). [9] S. V. Fleev, M. vnschilfgrde,t. Kotni,F. Leonrd, nd M. P. Desjrlis, Phys. Rev. 7, 33 (6). [] P. Ginnozzi et l., J. Phys.: Condens. Mtter, 3955 (9).