Theoretical ab initio Evolution of Satellite Intensity near Threshold for Cu K-shell transitions

Technische Samenvatting: Theoretische Ab Initio Evolutie van Satellietintensiteit nabij de Drempel voor Cu K-shell Transities

Probleemstelling
Theoretische berekeningen van atomaire structuren voor complexe systemen vertrouwen vaak op de abrupte benadering (sudden approximation, SA), die ervan uitgaat dat ionisatie of nucleair verval plaatsvindt op een tijdschaal die veel korter is dan de elektronische relaxatie. Hoewel de SA een uitstekende beschrijving biedt van hoogenergetische processen, is de geldigheid ervan twijfelachtig nabij ionisatiedrempels. In dit laagenergetische regime bezit het uitgaande elektron een lage kinetische energie en bevindt het zich gedurende een relatief lange tijd in de nabijheid van het atomaire systeem, wat de fundamentele aannames van de SA schendt. Als gevolg hiervan nadert het systeem een adiabatische limiet waarbij de elektronische structuur continu evolueert. Eerdere experimentele waarnemingen van dit effect, met name in röntgenabsorptiespectroscopie (XAS), werden beperkt door resolutiebeperkingen. Bovendien, hoewel theoretische modellen zoals die van Thomas en Roy bestaan om de evolutie van schokintensiteit (satellietlijnen) van het adiabatische naar het abrupte regime te voorspellen, is er een behoefte om deze modellen te valideren tegen hoogresolutie emissiedata en om te bepalen of standaard ab initio methoden deze evolutie nauwkeurig kunnen simuleren zonder empirische fittingparameters voor intensiteit.

Methodologie
De auteurs gebruikten een uitgebreide ab initio benadering met behulp van geavanceerde Multiconfiguration Dirac-Fock (MCDF) methoden, specifiek gebruikmakend van de MCDFGME-code. De methodologie omvatte verschillende belangrijke stappen:

1. Atomaire Structuurberekeningen: Alle Cu K-shell transities, inclusief radiatieve en niet-radiatieve snelheden, werden berekend. De atomaire golffuncties werden geoptimaliseerd met de optimized-level methode met volledige relaxatie voor radiatieve snelheden, terwijl niet-radiatieve snelheden gebruikmaakten van continuüm-elektrongolffuncties zonder relaxatie.
2. Bepaling van Schokkansbaarheid (Shake Probability): De kans op shake-up en shake-off processen werd theoretisch berekend. De totale schokkansbaarheid werd afgeleid uit het inwendig product van de initiële en finale toestand radiële golffuncties. De auteurs maakten onderscheid tussen single-electron schokkansbaarheden (met behulp van een binomiale distributiebenadering) en totale schokkansbaarheden, waarbij rekening werd gehouden met monopolaire excitaties.
3. Synthetische Spectrageneratie: Een synthetisch spectrum voor Cu K$\alpha_{1,2}$ transities werd gegenereerd door alle mogelijke vervalpaden te wegen met hun berekende snelheden en statistische gewichten.
4. Modellering van Energie-evolutie: Om de evolutie van lijnintensiteit nabij de drempel te simuleren, ontwikkelden de auteurs een modulatiefunctie. Dit omvatte het convolueren van de natuurlijke Lorentziaanse breedte van geïoniseerde niveaus met een Gaussisch profiel dat de energieverdeling van de excitatiebundel representeert. Een "Effectief Lorentziaans Profiel" werd gebruikt om de asymmetrische aard van de ionisatiedoorsnede nabij de drempel te verklaren.
5. Inbedding van Resonante Excitaties: In het besef dat ionisatiegegevens alleen de experimentele resultaten in het nabij-drempelgebied niet konden verklaren, incorporeerden de auteurs resonante excitaties van Cu(I) en Cu(II) oxidefasen (specifiek $1s \to 3d$ en $1s \to 4p$ transities) geïdentificeerd via X-ray Absorption Near-Edge Structure (XANES) data.
6. Fitting en Analyse: Theoretische resultaten werden vergeleken met experimentele data van Galambosi et al. [8]. Een Bayesiaanse analyse met nested sampling (nested fit) werd toegepast om de meest waarschijnlijke energies verschuivingen en breedtemultipliers voor verschillende vervalprocessen (diagram, shake-off, shake-up en excitatie) te bepalen om de simulatie in lijn te brengen met de experimentele observaties.

Belangrijkste Bijdragen

• Validatie van Ab Initio Methoden: De studie demonstreert dat standaard state-of-the-art ab initio methoden, zonder aanpasbare intensiteitsfittingparameters voor ionisatie, succesvol de evolutie van satellietintensiteit nabij ionisatiedrempels kunnen simuleren.
• Kwantificering van Schokkansbaarheden: De auteurs berekenden een totale schokkansbaarheid voor Cu K-shell ionisatie van ongeveer 25,97%, met een uitsplitsing van $\approx$13,08% voor shake-off en $\approx$12,89% voor shake-up.
• Identificatie van Drempel-oorsprong: Het werk identificeert dat de "onder-drempel" satellietintensiteit waargenomen in experimenten niet uitsluitend voortkomt uit schokprocessen, maar significant wordt beïnvloed door resonante excitaties in Cu(I) ($1s \to 4p$) en Cu(II) ($1s \to 3d$) oxidefasen.
• Verfijning van Theoretische Modellen: Het artikel vormt een rigoureuze test van de Thomas en Roy modellen. Het bevestigt de functionele vorm van het Thomas-model voor de adiabatic-to-sudden transitie en wijst op discrepanties in de parameters van het Roy-model wanneer deze worden gefit aan experimentele data, wat suggereert dat de experimentele data fysieke processen bevatten (zoals solid-state effecten) die niet volledig worden gevangen door de eenvoudige atomaire modellen.

Resultaten

• Overeenstemming met Experiment: De gesimuleerde intensiteitsratio's van satellietlijnen tot de totale K$\alpha_{1,2}$ intensiteit vertoonden een goede overeenstemming met experimentele data voor excitatie-energieën boven 9050 eV, waar alle vervalkanalen open zijn. De gesimuleerde ratio bij 10 keV werd gevonden op $0,262 \pm 0,00004$, consistent met de totale schokkansbaarheid berekend in de abrupte benadering.
• Drempelgedrag: Onder 9050 eV faalde de simulatie die enkel op ionisatie gebaseerd was in het reproduceren van de experimentele stijging in satellietintensiteit. Echter, na het opnemen van resonante excitaties van Cu oxiden, slaagde de simulatie erin de experimentele data in het bereik van 8975–9060 eV te reproduceren.
• Model Fitting: Het fitten van het Thomas-model aan de simulatie leverde parameters ($\Delta E$, $r$, $\mu_\infty$) op die in goede overeenstemming waren met zowel de experimentele fit als de theoretische berekeningen. De Roy-model fit vertoonde een redelijke nabijheid tot theoretische waarden, maar significante afwijkingen in bindingsenergieparameters, wat duidt op de beperkingen van het model om de volledige complexiteit van het experimentele spectrum te vatten.
• Spectrale Kenmerken: De inclusie van hogere-orde shake-up transities (bijv. $3s$, $3p$, $2p$ shake-off en $3d$ shake-up) veranderde de theoretische satellietlijn-vorm aanzienlijk vergeleken met eerdere Dirac-Fock berekeningen, wat het belang van gedetailleerde shake-berekeningen voor hoogresolutie metrologie benadrukt.

Significantie
Het artikel claimt dat de primaire significantie ligt in het aantonen dat standaard theoretische berekeningen de gemiddelde evolutie van spectrale intensiteit nabij ionisatiedrempels kunnen voorspellen zonder te vertrouven op empirische intensiteitsfitting. Door succesvol het satellietintensiteit (afkomstig van schokprocessen) te scheiden van directe diagram-intensiteit en de bijdrage van resonante excitaties in oxidefasen te identificeren, biedt het werk een completer fysiek beeld van de nabij-drempel regio. Het werk onderstreept het belang van gedetailleerde shake-up en shake-off berekeningen voor state-of-the-art röntgenspectroscopie en suggereert dat discrepanties tussen theorie en experiment in dit regime vaak toe te schrijven zijn aan solid-state effecten (zoals de aanwezigheid van oxiden) en de specifieke lijnvormen die in data-analyse worden gebruikt, in plaats van een falen van de onderliggende atomaire fysica modellen.