\textit{\textbf{First-principles}} description of pumped inelastic X-ray scattering: example of K-edge RIXS in graphite

Dit artikel presenteert een *ab initio* raamwerk dat de Bethe-Salpeter-vergelijking en real-time dichtheidsfunctionaaltheorie combineert om tijd-opgeloste resonante inelastische röntgenverstrooiing (RIXS) in optisch gepompte materialen te voorspellen, waarbij de nauwkeurigheid ervan wordt aangetoond door succesvol hoekafhankelijke K-rand RIXS-spectra in grafiet bij verschillende vertragingstijden te modelleren.

De kern

Stel je voor dat je probeert te begrijpen hoe een complexe machine werkt, zoals een vleugel. Je zou naar een liedje kunnen luisteren dat erop wordt gespeeld (dat is als normale spectroscopie), maar Resonant Inelastic X-ray Scattering (RIXS) is alsof je een specifieke toets aanslaat met een hamer gemaakt van licht, luistert naar het geluid dat het maakt, en vervolgens precies analyseert hoe de interne snaren en hamers in reactie daarop trilden. Het vertelt je niet alleen waar de machine van gemaakt is, maar ook hoe de onderdelen bewegen en met elkaar interageren.

Dit artikel introduceert een nieuw, superprecies computerprogramma dat precies voorspelt hoe deze "licht-hamer"-experimenten eruit zullen zien, zelfs wanneer de machine wordt geschud door een tweede, snellere puls van licht (zoals een pomp).

Hier is de uitsplitsing van hun werk met behulp van alledaagse analogieën:

1. Het Probleem: Het Onvoorspelbare Voorspellen

Wetenschappers zijn al lang in staat om "snapshots" van materialen te maken met röntgenstraling. Het voorspellen van hoe die snapshots er precies uitzien — vooral wanneer het materiaal wordt "gepompt" met een laser om het wakker te maken uit zijn slaap — is echter erg moeilijk.

• De Oude Manier: Eerdere computermodellen waren als het kijken naar een menigte mensen en ervan uitgaan dat iedereen stilstaat en alleen handelt. Ze misten hoe mensen (elektronen) daadwerkelijk elkaars hand vasthouden en samen bewegen (een fenomeen dat "excitonische effecten" wordt genoemd).
• De Nieuwe Manier: De auteurs hebben een nieuw framework gebouwd dat fungeert als een hogesnelheid, 3D-filmsimulator. Het kijkt niet alleen naar individuen; het kijkt hoe de hele menigte samen danst, rekening houdend met hoe ze aan elkaar trekken.

2. De Methode: Een Tweestapsdans

De onderzoekers combineerden twee krachtige instrumenten om hun simulatie te creëren:

• Stap 1 (De "Pomp"): Ze gebruikten een tool genaamd RT-TDDFT om te simuleren wat er gebeurt wanneer een laser op het materiaal valt. Stel je voor dat je een zaklamp op een trampoline schijnt; deze tool berekent hoe de trampoline stuiteren en hoe de mensen erop hun gewicht verplaatsen direct nadat het licht de trampoline raakt. Dit geeft hen een "niet-evenwichtige" kaart van waar de elektronen zich bevinden direct na de laserpuls.
• Stap 2 (De "Probe"): Ze gebruikten vervolgens de Bethe-Salpeter Vergelijking (BSE). Denk aan dit als een supernauwkeurig regelboek voor hoe röntgenstraling interageert met die stuiterende trampoline. Het berekent de complexe dans tussen het elektron dat eruit is gekicked en het "gat" (de lege ruimte) dat het heeft achtergelaten.

Door deze te combineren, kunnen ze de "echo" (de verstrooide röntgenstraling) voorspellen voor elke invalshoek van het licht en elke uitgaande hoek van het licht.

3. De Testcase: Grafiet (Het Potloodgrafiet)

Om te bewijzen dat hun methode werkt, testten ze het op grafiet (het spul waar een potlood van gemaakt is).

• Waarom Grafiet? Het is als een stapel vellen papier. De atomen binnen elk vel zijn stevig aan elkaar geplakt (zoals sterke lijm), maar de vellen zelf zitten slechts losjes aan elkaar vast (zoals een stapel losse papieren). Dit maakt het materiaal zeer "anisotroop", wat betekent dat het heel anders reageert afhankelijk van of je het van de zijkant of van bovenaf bekijkt.
• Het Resultaat: De computersimulatie voorspelde succesvol twee verschillende soorten "noten" die het grafiet zou spelen:
  • $\pi$ (Pi) noten: Deze komen van de elektronen die tussen de vellen bewegen (de losse papieren).
  • $\sigma$ (Sigma) noten: Deze komen van de elektronen die stevig binnen de vellen bewegen (de sterke lijm).
    De simulatie liet zien dat als je licht van de zijkant schijnt, je vooral de "lijm"-noten hoort. Als je het van bovenaf beschijnt, hoor je de "papier"-noten. Dit kwam perfect overeen met de experimenten in de echte wereld.

4. Het "Gepompte" Experiment: De Tafel Schudden

Het meest opwindende deel van het artikel is wat er gebeurt wanneer ze het grafiet "pompen" met een laser voordat ze het met röntgenstraling raken.

• De Analogie: Stel je voor dat het grafiet een kalme vijver is. De laserpomp is als het gooien van een steen in de vijver, wat rimpelingen veroorzaakt. De röntgenstraling is een sonar-ping die wordt verzonden om te zien hoe de rimpelingen het water hebben veranderd.
• De Bevinding: Wanneer het grafiet werd "gepompt", liet de simulatie zien dat de "noten" licht veranderden. Er verschenen nieuwe, zwakke geluiden in het lage energiebereik, en het volume van bestaande geluiden verschoof.
• De Les: De computer voorspelde dat zelfs een korte laserpuls de elektronische "stemming" van het materiaal verandert, waardoor een tijdelijke staat ontstaat die anders is dan de ruststand. De simulatie kwam zo goed overeen met de experimentele gegevens dat het deze subtiele verschuivingen kon zien, wat bewees dat de methode werkt voor "time-resolved" (frame-voor-frame film) studies.

Samenvatting

In eenvoudige termen zegt dit artikel: "We hebben een nieuw, zeer nauwkeurig computermodel gebouwd dat precies kan voorspellen hoe een materiaal zal reageren op röntgenstraling, zelfs wanneer dat materiaal wordt geschud door een laser."

Ze testten dit op grafiet, en de "voorspelling" van de computer kwam perfect overeen met het echte experiment, waarbij correct werd geïdentificeerd hoe de interne structuur van het materiaal (de strakke vellen versus de losse lagen) reageert op licht vanuit verschillende hoeken en op verschillende tijdstippen. Dit geeft wetenschappers een krachtig nieuw hulpmiddel om te begrijpen hoe materialen zich in realtime gedragen, zonder dat ze voor elke aanname dure experimenten hoeven uit te voeren.

Technische samenvatting

Technische Samenvatting: Eerste-principes beschrijving van gepompt inelastisch röntgenverstrooiingspectroscopie

Probleem en Context
Resonante inelastische röntgenverstrooiing (RIXS) is een krachtige spectroscopische techniek voor het onderzoeken van elektronische excitaties met element-, orbitaal- en momentumresolutie. Hoewel theoretische kaders voor RIXS zijn geëvolueerd van onafhankelijke deeltjeksbenaderingen naar veel-deeltjes perturbatietheorie (MBPT) met behulp van de Bethe-Salpeter-vergelijking (BSE) om excitonische effecten te vangen, blijft er een aanzienlijke kloof bestaan in het beschrijven van tijd-opgeloste (gepompt) RIXS. Bestaande theoretische modellen missen vaak het vermogen om tegelijkertijd rekening te houden met:

1. Excitonische effecten in zowel kern- als valentie-excitaties, die cruciaal zijn voor een nauwkeurige spectrale beschrijving.
2. Niet-evenwicht ladingdragerdistributies geïnduceerd door optische pompacties.
3. Willekeurige polarisatiegeometrieën voor inkomende en uitgaande fotonen, die essentieel zijn voor anisotrope materialen.

Tijd-opgeloste RIXS-experimenten (tr-RIXS) komen op, maar blijven beperkt; theoretische modellering vereist een kader dat elektron-gat-correlaties integreert met niet-evenwichtsdynamica.

Methodologie
De auteurs presenteren een ab initio kader dat Real-Time Time-Dependent Density Functional Theory (RT-TDDFT) combineert met de Bethe-Salpeter-vergelijking (BSE) om polarisatie- en tijd-opgeloste RIXS-spectra te voorspellen. De methodologie verloopt als volgt:

1. Theoretisch Formalisme: De aanpak is gebaseerd op de Kramers-Heisenberg-formule voor de dubbel-differentiaal doorsnede (DDCS). Het RIXS-proces wordt beschreven als een tweestapsmechanisme:

   • Absorptie: Een invallend röntgenfoton ($\omega_1$) creëert een kerngat-intermediairstoestand ($|n\rangle$).
   • Emissie: Een valentie-elektron vult het kerngat op, waarbij een foton ($\omega_2$) wordt uitgezonden en een definitieve elektron-gat-paar ($|f\rangle$) achterblijft.
     De DDCS wordt berekend met behulp van transitie-operatoren afgeleid van de BSE-eigenvectoren en eigenwaarden.

2. Integratie van Niet-Evenwichtstoestanden: Om het gepompte geval te modelleren, gebruikt het kader RT-TDDFT (specifiek binnen de exciting-code) om het Kohn-Sham-systeem te propageren onder een optische pomp-laser. Dit levert tijd afhankelijke, niet-evenwicht elektronische bezettingsgetallen op.

3. Workflow:

   • Grondtoestand: Een DFT-berekening vormt het startpunt.
   • BSE-berekeningen: Twee afzonderlijke BSE-berekeningen worden uitgevoerd:
     • Eén voor kernniveau-excitaties (om kern oscillatorsterktes $t^{(1)}$ te verkrijgen).
     • Eén voor valentie-/optische excitaties (om de finale toestandsroutes $t^{(2)}$ te verkrijgen).
   • Niet-evenwichtsextensie: Voor tr-RIXS worden de niet-evenwicht bezettingen uit RT-TDDFT gebruikt als restricties in de BSE-berekeningen.
   • Generalisatie van Polarisatie: De implementatie werd uitgebreid om willekeurige, onafhankelijke polarisatievectoren ($\mathbf{e}_1$ en $\mathbf{e}_2$) voor inkomend en uitgaand licht te verwerken, waarmee de eerdere beperking waarbij $\mathbf{e}_1 = \mathbf{e}_2$ werd doorbroken.

4. Computationele Details: Berekeningen werden uitgevoerd met de all-elektron full-potential package exciting. De studie maakte gebruik van de koolstof K-rand van grafiet, waarbij een GGA-PBE functionaal met van der Waals-correcties werd toegepast. De BSE omvatte 8 bezette en 15 onbezette banden voor optische spectra, en 40 geleidingsbanden voor röntgenspectra.

Belangrijkste Bijdragen

• Verenigd Kader: Het artikel vestigt een volledig parameter-vrij ab initio workflow die RT-TDDFT en BSE combineert om tr-RIXS te simuleren, waarbij zowel veel-deeltjes excitonische effecten als pomp-geïnduceerde niet-evenwichtsdynamica worden gevangen.
• Polarisatiegeneralisatie: De auteurs hebben de RIXS-implementatie gegeneraliseerd om ondersteuning te bieden aan willekeurige polarisatieoriëntaties voor zowel inkomende als verstrooide fotonen, wat onderzoek naar complexe verstrooiingsgeometrieën en gemengde polarisaties mogelijk maakt.
• Implementatie in exciting: De methode is geïmplementeerd binnen de full-potential code exciting, gebruikmakend van de pyBRIXS package voor cross-section berekeningen.

Resultaten: Grafiet K-rand RIXS
Het kader werd toegepast op grafiet, een halfmetaal met een sterke elektronische anisotropie tussen in-plane ($\sigma$) en out-of-plane ($\pi$) toestanden.

• Evenwichtsspectra:

  • De berekende absorptiespectra (zowel optisch als kernniveau) reproduceerden experimentele kenmerken, wat de noodzaak van BSE bevestigt om excitonische effecten te vangen, met name voor $\pi \to \sigma^*$ en $\sigma \to \pi^*$ transities.
  • Anisotropie: Evenwicht RIXS-spectra vertoonden een sterke polarisatieafhankelijkheid. In-plane polarisatie ($\mathbf{e} \parallel x$) benadrukte hoogenergetische $\sigma$-afgeleide excitaties ($>8$ eV), terwijl out-of-plane polarisatie ($\mathbf{e} \parallel z$) lage-energie $\pi$-afgeleide excitaties accentueerde.
  • Verstrooiingsgeometrie: Simulaties van een $30^\circ$ invalshoek (overeenkomend met experimentele opstellingen) onthulden een niet-lineaire menging van $\pi$- en $\sigma$-bijdragen, wat interferentie-effecten en versterkte $\pi \to \pi^*$ transities bij lage energieverlies liet zien.

• Gepompt (Niet-Evenwicht) Spectra:

  • Simulaties werden uitgevoerd voor een systeem dat gepompt wordt door een 266 nm laser (45 fs duur) met een vertragingstijd van 150 fs.
  • Spectrale Veranderingen: De gepompt spectra vertoonden sterke gelijkenis met het evenwichtgeval, wat wijst op matige modificaties. Echter, een duidelijke extra lage-energieverlies feature (enkele eV) kwam naar voren voor beide polarisaties, toegeschreven aan modificaties in de elektronische distributie nabij het Fermi-niveau.
  • Gemengde Polarisaties: In opstellingen met gemengde polarisaties vertoonden de niet-evenwicht spectra een herverdeling van de spectrale massa en verbreedde kenmerken vergeleken met het evenwicht.
  • Vergelijking met Experiment: De theoretische resultaten reproduceerden experimentele piekposities en algemene kenmerken goed. Kleine afwijkingen (bijv. piektijden en specifieke verschuivingen) werden toegeschreven aan elektron-fononkoppeling, die niet is opgenomen in het huidige statische kader, en mogelijke zelfabsorptie-effecten in experimenten.

Betekenis en Claims
De auteurs claimen dat hun werk een algemeen en volledig parameter-vrij pad biedt om RIXS-spectra in complexe materialen te berekenen, variërend van metalen tot grote-gap halfgeleiders, voor willekeurige verstrooiingsgeometrieën. Door succesvol de karakteristieke anisotropie van grafiet en de evolutie ervan onder optische pompactie te reproduceren, demonstreert de methode haar vermogen om:

1. Distincte $\pi$- en $\sigma$-afgeleide kenmerken te onderscheiden.
2. De interactie tussen elektronische structuur en polarisatie te vangen.
3. Tijd-afhankelijke fenomenen te modelleren geïnduceerd door optische excitatie.

Het artikel concludeert dat deze aanpak een robuust instrument biedt voor het interpreteren van tr-RIXS experimenten en het begrijpen van niet-evenwichtfenomenen in vaste stoffen, waarbij specifiek de directionele afhankelijkheid van fundamentele excitaties wordt benadrukt. De auteurs blijven bescheiden over discrepanties met het experiment en schrijven deze toe aan bekende beperkingen zoals de uitsluiting van elektron-fononkoppeling en zelfabsorptie-effecten.