Interplay of non-local transport and local scattering during electron thermalization and spatial equilibration in laser-excited metals

Cet article emploie une équation de transport de Boltzmann reformulée pour démontrer que, si le transport non local accélère la thermalisation apparente à la surface irradiée en éliminant les porteurs atermaux, il retarde simultanément l'équilibrage complet de l'ensemble du système électronique, révélant une interaction complexe entre transport et diffusion qui varie selon la position et l'énergie.

L'essentiel

Imaginez une piste de danse bondée à l'intérieur d'un bloc de métal. Soudain, une impulsion laser ultra-rapide frappe la surface, comme un DJ qui lâche un beat massif que seuls les danseurs juste devant la scène peuvent entendre. Ces « danseurs » (les électrons) s'excitent, bondissent et bougent frénétiquement, tandis que ceux plus loin dans la salle restent assis calmement.

Cet article traite de ce qui se passe ensuite : comment ces danseurs excités se calment-ils, et comment l'énergie se propage-t-elle dans toute la pièce ?

Les deux forces principales en jeu

Les chercheurs ont découvert que deux phénomènes principaux se produisent simultanément, et qu'ils s'opposent souvent l'un à l'autre :

1. La « Diffusion Locale » (Le choc sur la piste de danse) : Les danseurs excités se cognent entre eux et contre les murs de la pièce (la structure atomique du métal). C'est comme un mosh pit chaotique où tout le monde finit par ralentir pour danser sur un rythme synchronisé et calme. C'est la thermalisation.
2. Le « Transport Non-Local » (La vague de la foule) : Parce que le laser n'a frappé que l'avant, les danseurs à l'avant sont entassés et énergiques, tandis que l'arrière est vide. Naturellement, les danseurs énergiques de l'avant commencent à courir vers l'arrière pour combler l'espace vide. C'est le transport.

La grande surprise : L'effet de « l'illusion »

La découverte la plus intéressante de l'article est un petit tour de passe-passe de la lumière.

Si vous vous tenez à l'entrée (la surface) et que vous observez les danseurs, vous pourriez vous dire : « Wow, ils se sont calmés très vite ! » Les chercheurs ont découvert que le transport donne en fait l'impression que l'avant du métal se refroidit plus rapidement.

Pourquoi ? Parce que les danseurs excités sont littéralement en train de fuir l'entrée pour courir vers l'arrière de la pièce. Ils ne sont pas forcément en train de se calmer à l'avant ; ils sont simplement en train de quitter l'avant. Ainsi, si vous regardez seulement la surface, il semble que le système ait atteint un équilibre paisible très rapidement.

Cependant, l'article soutient que le système entier n'est pas encore réellement calme. Les danseurs courent toujours, essayant de remplir l'arrière de la pièce. La « paix » à l'avant est une illusion causée par le mouvement de la foule qui s'éloigne. Le système complet ne devient véritablement calme qu'une fois que les danseurs se sont répartis uniformément de l'avant vers l'arrière.

L'analogie de la « Fenêtre d'Énergie »

Les chercheurs ont également observé des groupes spécifiques de danseurs en fonction de leur degré d'excitation (leurs niveaux d'énergie).

• Le groupe « Modérément Excité » (Basse Énergie) : Ce sont les danseurs qui sont juste un peu agités. Leur mouvement est principalement contrôlé par la vague de la foule (le transport). Ils se contentent principalement de se déplacer de l'avant encombré vers l'arrière vide.
• Le groupe « Follement Excité » (Haute Énergie) : Ce sont les danseurs qui sautent sur les tables. Leur comportement est principalement contrôlé par les collisions entre eux (la diffusion). Ils perdent leur énergie sauvage en percutant les autres très rapidement, quel que soit l'endroit où ils se trouvent dans la pièce.

L'essentiel

L'article conclut que l'on ne peut pas comprendre ce qui se passe dans un métal frappé par un laser en regardant seulement la surface ou en supposant que tout se passe en un seul point.

• À la surface : On a l'impression que les choses se calment vite parce que les électrons « chauds » s'enfuient.
• À l'intérieur du métal : Le système est en réalité toujours chaotique car ces électrons se propagent, créant un nouveau type de déséquilibre pendant leur trajet vers l'arrière.

Les chercheurs ont construit un nouveau modèle mathématique (comme une simulation ultra-précise de la piste de danse) qui suit à la fois les collisions (la diffusion) et la course (le transport) en même temps. Cela aide les scientifiques à comprendre que, dans les métaux épais, le « refroidissement » ne consiste pas seulement à ralentir ses pieds ; c'est aussi une question de l'endroit où l'on se tient dans la pièce.

Résumé technique

Résumé Technique : Interaction entre le transport non local et la diffusion locale dans les métaux excités par laser

Problématique
L'excitation par laser ultrafast des métaux induit un hors-équilibre électronique tant au niveau local (distribution en énergie) que non local (dans l'espace), particulièrement lorsque l'épaisseur de l'échantillon dépasse la profondeur de pénétration optique. La dynamique subséquente est alors régie par l'interaction entre le transport non local et la diffusion locale (électron-électron et électron-phonon), mais les descriptions microscopiques combinées de ces processus restent rares. Les approches théoriques existantes reposent généralement sur des hypothèses simplificatrices : elles supposent soit un chauffage homogène pour résoudre la diffusion en détail, soit conservent le transport tout en employant des modèles de diffusion approximatifs (par exemple, des approximations par temps de relaxation). Ces approximations entravent la capacité à distinguer les contributions spécifiques de la thermalisation intrinsèque et du transport spatial dans les dynamiques de porteurs mesurées, notamment dans les expériences résolues en énergie comme la spectroscopie de photoémission (PES) résolue en temps.

Méthodologie
Les auteurs présentent un cadre de Boltzmann spatialement résolu qui intègre simultanément la diffusion locale et le transport non local. Le cœur de la méthode est une reformulation de l'équation de transport de Boltzmann dans l'espace des énergies utilisant des intégrales de collision microscopiques complètes. Les caractéristiques méthodologiques clés incluent :

• Description à deux branches : Pour résoudre le mouvement spatial sans résolution complète de l'impulsion, l'impulsion de l'électron est projetée sur la normale à la surface (direction $z$). La dispersion est divisée en deux branches : les électrons se déplaçant vers le volume (bulk) ($v_z > 0$) et les électrons se déplaçant vers la surface ($v_z < 0$).
• Intégrales de collision complètes (FCI) : Le modèle utilise des FCI complètes pour la diffusion électron-électron et électron-phonon dans l'approximation de l'angle aléatoire (random-k). Cela évite les inexactitudes des approximations par temps de relaxation tout en restant plus gérable numériquement que les calculs complets résolus en impulsion.
• Couplage entre diffusion et transport : Le cadre prend explicitement en compte les événements de diffusion qui peuvent changer la direction du mouvement (changement de branche) et traite le transport comme l'évolution spatiale de ces deux branches. Les conditions aux limites aux surfaces avant et arrière sont gérées par la réflexion des électrons entre les branches.
• Configuration de simulation : Le modèle est appliqué à un métal à électrons libres de type aluminium (correspondant aux paramètres de premiers principes pour la vitesse de Fermi et l'absorption) d'une épaisseur de 50 nm. Les simulations couvrent une gamme de fluences laser ($0,05$à$15 \text{ J m}^{-2}$) et comparent une excitation homogène à une excitation de Beer-Lambert inhomogène à densités d'énergie absorbée identiques.

Contributions Clés
La contribution primaire de ce travail est le développement d'un cadre microscopique cohérent qui comble le fossé entre l'équilibrage spatial et la diffusion dans l'espace des énergies. En combinant directement la diffusion et le mouvement des électrons, le modèle couvre implicitement toute la gamme allant du transport balistique au transport diffusif. Cela permet de quantifier directement comment les processus de transport modifient les dynamiques de thermalisation apparentes, une distinction auparavant difficile à isoler en raison de la dépendance aux modèles simplifiés.

Résultats
L'étude apporte plusieurs enseignements critiques sur la dynamique des électrons chauds :

1. Accélération de la thermalisation de surface : Le transport accélère la thermalisation apparente observée à la surface irradiée. En extrayant les porteurs athermiques de la région de surface vers le volume, le transport agit comme un canal de relaxation supplémentaire. Cet effet est plus prononcé à faibles fluences, là où la diffusion intrinsèque électron-électron est plus lente.
2. Retard de l'équilibrage global : Bien que la thermalisation de surface paraisse plus rapide, la redistribution spatiale de l'énergie retarde l'équilibrage complet de l'ensemble du système électronique. La présence de gradients spatiaux entraîne une population persistante de porteurs athermiques à travers la profondeur de l'échantillon pendant des durées plus longues par rapport à une excitation homogène.
3. Dynamique dépendante de la profondeur : Le mécanisme dominant pilotant la dynamique varie considérablement selon la position et l'énergie :
   • Surface avant : La dynamique est pilotée par une compétition entre la diffusion électron-électron et le transport (extraction des porteurs).
   • Surface arrière : La dynamique est dominée par les processus de transport, car l'excitation laser directe y est négligeable. Les électrons athermiques s'y accumulent de manière transitoire en raison de la réflexion.
   • Dépendance en énergie : Près du niveau de Fermi, le transport est le principal moteur des changements de densité spectrale. À des énergies plus élevées (par exemple, 1,0 eV au-dessus du niveau de Fermi), la diffusion électron-électron domine la décroissance, se produisant sur l'échelle de temps de l'impulsion laser.
4. Évolution de la densité spectrale : L'analyse des densités spectrales (imitant les observables de la PES) révèle que l'interprétation des temps de relaxation dépend fortement de la fenêtre d'énergie et de la profondeur sondées. Par exemple, à la surface arrière, les changements de densité spectrale sont presque entièrement pilotés par le transport plutôt que par la diffusion locale.

Signification et Revendications
Les auteurs affirment que leur travail améliore la compréhension du hors-équilibre électronique dans les systèmes optiquement épais et pilotés par laser. Ils soulignent que la thermalisation électronique dans ces systèmes ne peut être interprétée comme un processus purement local. Une conclusion clé est la distinction entre les observables sensibles à la surface et l'état global du système électronique : les mesures de surface peuvent indiquer une thermalisation rapide car le transport extrait les porteurs athermiques du volume sondé, même si le système électronique complet demeure hors équilibre en raison de la redistribution spatiale.

L'article conclut que cette distinction est particulièrement critique pour les mesures ultrafast résolues en énergie, où le canal de relaxation dominant dépend à la fois de la fenêtre d'énergie sondée et de la profondeur de sondage. Le cadre présenté offre une voie microscopique pour connecter le transport, la diffusion et la dynamique des porteurs chauds accessible expérimentalement, offrant une base plus précise pour interpréter les données pertinentes pour les futures applications électroniques impliquant des porteurs chauds.