All-electron dark matter-electron scattering with random-phase approximation dielectric screening and local field effects

Deze paper introduceert een volledig elektronische framework in de open-source code QCDark2 voor het berekenen van donkere-materie-elektron verstrooiingskansen in kristallijne materialen, waarbij wordt aangetoond dat lokale veld-effecten een cruciale rol spelen bij het nauwkeurig voorspellen van terugstootspectra en gevoeligheidsgrenzen.

De kern

Het Onzichtbare Spook en de Elektronen: Een Reis door het Kristal

Stel je voor dat je probeert een onzichtbaar spook (donkere materie) te vangen. Dit spook is heel licht en beweegt snel, maar het is zo flauw dat het nauwelijks met de wereld om het heen praat. Wetenschappers bouwen daarom supergevoelige valkuilen: kristallen zoals silicium (zoals in je computerchip) of diamant. Als het spook tegen een elektron in het kristal stoot, zou het een kleine vonk moeten geven.

Maar hier zit de twist: kristallen zijn geen lege kamers. Ze zijn vol met elektronen die als een dichte menigte dansen. Als het spook erdoorheen rent, reageert die menigte. Ze schuiven uit de weg, trillen en vormen golven. Dit noemen we dielektrische afscherming.

Deze paper, geschreven door een team van onderzoekers, vertelt ons hoe we die menigte beter moeten begrijpen om het spook echt te kunnen vangen. Hier is de uitleg in simpele taal:

1. Het Probleem: De "Gladde" vs. de "Ruwe" Wereld

Vroeger dachten wetenschappers dat ze de elektronen in een kristal als een gladde, homogene soep konden behandelen. Ze gebruikten wiskundige benaderingen die de details negeerden.

• De analogie: Stel je voor dat je door een drukke stad loopt. De oude manier was alsof je dacht dat de stad een gladde, grijze muur is. Je wist dat je erdoorheen moest, maar je zag de mensen niet.
• De realiteit: In werkelijkheid zijn er straten, pleinen, en mensen die in groepjes staan. Als je hard loopt, bots je tegen specifieke groepen aan. Als je langzaam loopt, reageren ze anders.

De onderzoekers zeggen: "We moeten stoppen met het behandelen van de stad als een muur. We moeten de individuele mensen (elektronen) en hun lokale groepjes zien."

2. De Oplossing: "Lokale Veld-effecten" (LFE)

De paper introduceert een nieuw concept: Lokale Veld-effecten.

• Wat is het? Het is het effect dat ontstaat door de ongelijkheid in de elektronenwolk. Soms zitten elektronen heel dicht bij elkaar (zoals in een drukke kerk), soms ver uit elkaar.
• De analogie: Als je door een drukke menigte loopt, wordt je afgeremd als je in een dichte groep terechtkomt. Als je door een open plek loopt, gaat het sneller. De "lokale veld-effecten" zijn de manier waarop de menigte op die specifieke drukte reageert.
• Waarom is dit belangrijk? De onderzoekers ontdekten dat als je deze lokale effecten negeert, je de kans om het spook te vangen verkeerd berekent.
  • Bij hoge snelheid (het spook rent hard): De lokale effecten remmen de botsing af. Als je dit negeert, denk je dat je meer kans hebt om het spook te vangen dan je eigenlijk hebt. De paper zegt: "We moeten onze verwachtingen verlagen met 20% tot 50%."
  • Bij lage snelheid (het spook loopt langzaam): De lokale effecten zorgen ervoor dat de "golven" in de menigte (plasmonen) breder worden. Dit verandert het patroon van de vonk die je ziet.

3. Twee Soorten Spookjagers

De paper kijkt naar twee soorten donkere materie:

1. Het "Stille" Spook (Halo DM): Dit is de donkere materie die normaal rond de aarde zweeft. Het beweegt langzaam. Voor dit spook is het belangrijk om de hoge snelheid van de botsing goed te begrijpen. Hier bleek dat de lokale effecten de kans op detectie verkleinen.
2. Het "Boosted" Spook (SRDM): Dit zijn spookjes die ergens anders (bijvoorbeeld in de zon) een flinke duw hebben gekregen en nu als een kogel door het kristal vliegen. Voor deze snelle spookjes is het belangrijk om de lage snelheid en de golven in de menigte goed te begrijpen. Hier verandert de lokale effecten het type vonk dat je ziet, waardoor het makkelijker kan zijn om ze te onderscheiden van ruis.

4. De Nieuwe Tool: QCDark2

De onderzoekers hebben een nieuwe computercode gebouwd, genaamd QCDark2.

• De analogie: Stel je voor dat je eerder een simpele schets van de stad had (een plattegrond). QCDark2 is een 3D-simulatie met virtual reality, waar je elke persoon in de menigte kunt zien en hoe ze reageren op jou.
• Ze hebben deze tool gebruikt om verschillende kristallen te testen: Silicium, Germanium, Gallium-Arsenide, Siliciumcarbide en Diamant.
• Ze hebben ontdekt dat voor sommige materialen (zoals Germanium) je zelfs naar de "diepe kelders" van het atoom moet kijken (de 3d-orbitalen) om de botsing goed te begrijpen. Eerdere methoden keken daar niet naar.

5. Waarom doet dit er toe?

Als je een jager bent die op een zeldzaam dier jaagt, wil je niet een valk bouwen die te groot is (waardoor je het dier mist) of te klein (waardoor je het dier niet kunt vangen).

• Door de "lokale veld-effecten" mee te nemen, maken de onderzoekers de valk (de detector) preciezer.
• Ze zeggen: "Als je deze details negeert, denk je dat je een spook hebt gezien, terwijl het misschien gewoon ruis was. Of je denkt dat je het niet kunt vinden, terwijl je het juist wel kunt vinden als je de juiste instellingen gebruikt."

Kortom:
Deze paper is een handleiding voor de toekomst van spookjacht. Ze zeggen: "Stop met het benaderen van kristallen als gladde blokken. Kijk naar de ruwe, lokale details van de elektronen. Alleen dan weten we precies waar we moeten zoeken en wat we moeten verwachten als we eindelijk een teken van donkere materie zien."

Het is alsof ze van een wazige foto een haarscherpe foto hebben gemaakt, zodat we eindelijk kunnen zien wat er echt gebeurt in de microscopische wereld van donkere materie.

Technische samenvatting

Probleemstelling

Directe detectie-experimenten voor donkere materie (DM) in vaste stofmaterialen (zoals silicium, germanium, etc.) vereisen nauwkeurige voorspellingen van de verstrooiingskans tussen donkere materie en elektronen. Huidige berekeningen kampen met twee fundamentele tekortkomingen die de nauwkeurigheid beperken:

1. All-electron behandeling: Voor hoge impuls-overdrachten (vele keren de kristalimpuls) zijn benaderingen met pseudopotentialen onvoldoende. Een volledige behandeling van alle elektronen is noodzakelijk, maar computationally intensief.
2. Dielektrische afscherming: De respons van het materiaal wordt bepaald door de dielektrische functie. Eerdere methoden gebruikten vaak analytische benaderingen (zoals de Thomas-Fermi of Lindhard benadering) of berekenden de dielektrische functie binnen de Random Phase Approximation (RPA) maar verwaarloosden lokale veld-effecten (LFEs). LFEs ontstaan door microscopische inhomogeniteiten in de elektronische respons en zijn cruciaal voor zowel lage impuls (plasmon-resonantie) als hoge impuls (meerdere Brillouin-zones).

Geen van de bestaande codes (zoals QEDark, DarkELF, EXCEED-DM, of QCDark1) combineerde tot nu toe tegelijkertijd een all-electron behandeling met een volledige, ab initio RPA-berekening van de dielektrische functie inclusief LFEs over het volledige relevante impuls- en energiedomein. Dit leidt tot onnauwkeurigheden bij het voorspellen van de detectiebereik voor zowel halo-DM (niet-relativistisch) als "boosted" DM (relativistisch, bijvoorbeeld zon-gerelateerd).

Methodologie

De auteurs hebben een nieuw framework ontwikkeld, geïmplementeerd in de open-source code QCDark2, dat de volgende stappen combineert:

• All-electron DFT: De elektronische structuur wordt berekend met behulp van de Python-pakket PySCF, gebruikmakend van gelokaliseerde Gauss-type orbitalen (in plaats van vlakke golven met pseudopotentialen). Dit zorgt voor een nauwkeurige beschrijving van de elektronische toestanden over een breed energiedomein (tot 50 eV, of 150 eV voor diamant).
• RPA Dielektrische Functie: De longitudinale dielektrische functie wordt berekend binnen de Random Phase Approximation (RPA) op basis van Kohn-Sham-golf functies.
• Inclusie van Lokale Veld-effecten (LFEs): De kern van de methode is het oplossen van de volledige RPA-matrix $\epsilon_{\vec{G}\vec{G}'}(\omega, \vec{q})$ in de reciproque ruimte.
  • De matrixinversie van deze matrix (waarbij $\vec{G}$ en $\vec{G}'$ reciproque roostervectoren zijn) middelt over de microscopische inhomogeniteiten en introduceert de LFEs.
  • Omdat het berekenen van alle off-diagonale elementen voor hoge impulsen zeer kostbaar is, gebruiken de auteurs een composiet aanpak: LFEs worden volledig meegenomen tot een impuls-cutoff $q_{LFE}$ (ongeveer 7-8 $\alpha m_e$), en daarboven worden de niet-diagonale elementen verwaarloosd ($\epsilon_{noLFE}$).
• Materialen: De methode is toegepast op Silicium (Si), Germanium (Ge), Gallium Arsenide (GaAs), Silicium Carbide (SiC) en Diamant (C).
• Scattering Rate: De differentiaale verstrooiingsrate wordt berekend voor zowel halo-DM als voor "boosted" DM (zoals zon-gerelateerde DM), waarbij de dynamische structuurfactor $S(\omega, q)$ direct wordt afgeleid van de berekende dielektrische functie via de verliesfunctie $\text{Im}[-1/\epsilon]$.

Belangrijkste Bijdragen

1. QCDark2 Code: De publicatie introduceert QCDark2, de eerste open-source code die een all-electron DFT-benadering combineert met een volledige RPA-dielektrische functie inclusief LFEs voor hoge impulsen.
2. Systematische Vergelijking: Voor het eerst worden RPA-dielektrische functies met LFEs gepresenteerd voor een breed scala aan semiconductoren en isolatoren (Si, Ge, GaAs, SiC, Diamant).
3. Kwantificering van LFEs: Het werk demonstreert kwantitatief hoe LFEs de verstrooiingskansen beïnvloeden in verschillende regimes, wat eerder vaak werd verwaarloosd of slechts benaderd.
4. Validatie: De resultaten worden gevalideerd tegen experimentele data (zoals brekingsindex en EELS) en vergeleken met bestaande codes (DarkELF, QEDark, EXCEED-DM).

Resultaten

• Invloed op de Dynamische Structuurfactor:
  • Lage impuls: LFEs veroorzaken een verbreding van de plasmonpiek in de verliesfunctie, wat beter overeenkomt met experimentele data dan berekeningen zonder LFEs.
  • Hoge impuls: LFEs onderdrukken de dynamische structuurfactor aanzienlijk bij hoge impulsen (boven $\sim 3-4 \alpha m_e$). Voor Silicium wordt de structuurfactor bij lage energieën tot 15% van de waarde zonder LFEs gereduceerd.
• Invloed op Elektronen-Recoil Spectra:
  • Voor halo-DM (niet-relativistisch) leiden LFEs tot een reductie van de voorspelde verstrooiingsrate met 20-50% voor massa's boven een paar MeV, voornamelijk door de onderdrukking bij hoge impuls.
  • Voor boosted DM (bijv. zon-gerelateerd) is het effect complexer: LFEs verbreden de plasmonpiek, wat leidt tot een herschikking van het spectrum naar lagere en hogere energieën (meer elektron-gat paren dan de piek), maar de totale geïntegreerde rate verandert minder drastisch dan bij halo-DM.
• Vergelijking met Andere Afschermingsmodellen:
  • Analytische modellen (zoals Modified Thomas-Fermi) zijn onnauwkeurig voor boosted DM vanwege divergenties bij lage impuls.
  • De Lindhard-benadering doet het redelijk goed voor plasmon-regio's, maar mist materiaalspecifieke kenmerken (zoals interband-matrixelementen en 3d-orbitaal-overgangen in Ge) die door de all-electron RPA-berekening wel worden gevangen.
• Gevoeligheidsprojecties:
  • De projecties tonen aan dat het negeren van LFEs de gevoeligheid van experimenten voor halo-DM overschat.
  • Voor Germanium (Ge) is de rol van all-electron effecten cruciaal: de 3d-orbitalen (die in pseudopotentialen vaak worden "bevroren") dragen significant bij aan de verstrooiing bij hoge energieën, wat de gevoeligheid voor zware mediators verhoogt.

Betekenis en Conclusie

Dit werk markeert een belangrijke stap voorwaarts in de theoretische modellering van donkere materie-detectie in vaste stoffen. Door de combinatie van all-electron behandeling en volledige RPA-afscherming met lokale veld-effecten, biedt QCDark2 de meest accurate beschikbare voorspellingen voor elektronen-recoil spectra.

De belangrijkste implicaties zijn:

1. Correctie van eerdere overschattingen: Experimenten moeten rekening houden met een reductie van 20-50% in hun verwachte signaal voor halo-DM als ze gebruikmaken van deze nauwkeurigere afschermingsmodellen.
2. Unieke signatuur voor Boosted DM: LFEs veranderen de vorm van het spectrum voor boosted DM, wat kan helpen bij het onderscheiden van achtergrondruis.
3. Ressources voor de gemeenschap: De auteurs hebben de vooraf berekende dielektrische functies en de code (QCDark2) openbaar gemaakt, wat een standaard zal worden voor toekomstige studies naar donkere materie in diverse kristalstructuren.

Samenvattend stelt dit papier dat voor een nauwkeurige interpretatie van toekomstige data van experimenten zoals SENSEI, DAMIC-M en SuperCDMS, het onontbeerlijk is om lokale veld-effecten en all-electron effecten te incorporeren in de berekening van de dielektrische respons.