Computational Methods of Wave Propagation for Semiclassical Models of High Harmonic Generation in Bulk Solids

Licht, wanneer het zwak is, gaat door materialen zoals glas of silicium grotendeendeels onveranderd heen, waarbij het slechts buigt of vertraagt. Maar wanneer licht tot extreme intensiteiten wordt gedreven, begint het zich op wilde en complexe manieren te gedragen, waardoor het materiaal gedwongen wordt op een niet-lineaire wijze te reageren. Dit is het domein van de sterkveld-fysica, waar een enkele lichtflits elektronen uit hun atomen kan rukken of ze zo heftig kan doen schudden dat ze nieuwe kleuren licht uitzenden die niet aanwezig waren in de oorspronkelijke bundel. Een van de meest fascinerende uitkomsten van dit proces is hoge harmonische generatie, een fenomeen waarbij een materiaal een foton met lage energie absorbeert en dit opnieuw uitzendt als een foton met een veel hogere energie, waardoor de frequentie van het licht effectief wordt vermenigvuldigd. In gassen is dit proces al decennia goed begrepen en wordt het gebruikt om ultrakorte lichtpulsen te creëren voor het bestuderen van de snelste gebeurtenissen in de natuur. Echter, toen wetenschappers in 2011 voor het eerst ditzelfde effect in kristallen demonstreerden, openden zij een nieuwe deur naar het begrijpen van de elektronische structuur van materialen, maar introduceerden zij ook een nieuwe laag van complexiteit. In een vaste stof interageert het licht niet alleen met een enkel atoom en vertrekt het weer; het reist door een dichte, geordende rooster van miljarden atomen, en het licht zelf verandert terwijl het beweegt, waardoor de puls die het proces aanstuurt, van vorm verandert.

Een team van onderzoekers aan de Princeton University en het Max Planck Institute for Quantum Optics heeft een nieuwe manier ontwikkeld om deze reis te simuleren, door een computermodel te creëren dat bijhoudt hoe intense lichtpulsen door een kristallijn vaste stof reizen, terwijl het tegelijkertijd berekent hoe de elektronen van het kristal reageren. Hun werk richt zich op silicium, een materiaal dat in bijna elk elektronisch apparaat te vinden is, en maakt gebruik van een krachtige sturende laserpuls die slechts 24 femtoseconden lang is, een duur die zo kort is dat licht slechts enkele micrometers aflegt tijdens de gehele gebeurtenis. De onderzoekers bouwden een raamwerk dat de fundamentele vergelijkingen oplost die bepalen hoe elektrische en magnetische velden door de ruimte bewegen, gekoppeld aan vergelijkingen die het kwantummechanische gedrag van elektronen in het kristal beschrijven. Deze aanpak stelt hen in staat om in een virtuele omgeving te zien hoe de lichtpuls evolueert wanneer deze het silicium binnentreedt, hoe deze nieuwe kleuren licht genereert, en hoe die nieuwe kleuren meereizen met de oorspronkelijke puls, waarmee ze ermee interfereren en de vorm ervan veranderen.

De kern van hun ontdekking is dat het pad dat het licht door het materiaal aflegt, enorm veel uitmaakt. Wanneer de onderzoekers de generatie van deze hoogfrequente harmonischen simuleerden zonder rekening te houden met de reis van het licht door het kristal, zagen zij een specifiek patroon van kleuren. Echter, toen zij het licht door het materiaal lieten voortplanten, veranderde het resulterende spectrum drastisch. De intensiteit van de verschillende harmonische kleuren verschoof, de scherpte van de spectrale lijnen verbreedde, en de centrale frequenties van de kleuren bewogen lichtjes. Deze veranderingen waren niet willekeurig; ze waren het directe resultaat van de wijze waarop de lichtpuls zichzelf tijdens de reis van vorm veranderde. De onderzoekers ontdekten dat de puls zichzelf binnen het materiaal kan focussen, waardoor hij in het centrum intenser wordt, wat op zijn beurt meer harmonischen genereert. Ze observeerden ook dat het licht dat reflecteerde aan de voorkant van het kristal er anders uitzag dan het licht dat helemaal door naar de achterkant van het kristal was gereisd. Het getransmitteerde licht, dat vaak gemakkelijker te meten is in experimenten, vertoonde een ander profiel dan het gereflecteerde licht, wat suggereert dat wat wetenschappers in een laboratorium zien, een mengeling is van de intrinsieke elektronische eigenschappen van het materiaal en de effecten van de reis van het licht door de bulk van het materiaal.

Om dit te bereiken, moesten het team aanzienlijke computationele hindernissen overwinnen. Ze gebruikten een methode die het licht behandelt als een volledige elektromagnetische golf, waarbij zowel de elektrische als de magnetische componenten worden bijgehouden terwijl ze door het kristal naar voren en naar achteren bewegen. Dit onderscheidt zich van oudere methoden die vaak aannemen dat het licht slechts in één richting beweegt of de interactie vereenvoudigen tot een statische respons. Door de volledige vergelijkingen op te lossen, konden zij de subtiele wisselwerking tussen de sturende laser en de nieuw gegenereerde harmonischen vastleggen. Ze modelleerden het siliciumkristal met behulp van gegevens afgeleid van kwantummechanische berekeningen van de atomaire structuur, waarbij ze specif specifiek keken naar hoe de energieniveaus van elektronen veranderen terwijl ze door het rooster van het kristal bewegen. Ze ontdekten dat de elektronen in het silicium voor hun gekozen materiaal en lasercondities op een manier reageerden die accuraat beschreven kon worden door slechts twee hoofdenergiestanden te overwegen: de valentieband waar elektronen normaal gesproken zitten, en de conductieband waar ze vrij kunnen bewegen. Dit vereenvoudigde model stelde hen in staat de complexe simulaties efficiënt uit te voeren terwijl de essentiële fysica behouden bleef.

De onderzoekers verkenden ook hoe het vermogen van het materiaal om licht te buigen, bekend als de brekingsindex, de uitkomst beïnvloedt. Ze ontdekten dat de wijze waarop zij de respons van het materiaal op verschillende kleuren licht modelleerden, cruciaal was. Als zij de absorptie van licht op de specifieke energie inclusieerden waarbij het siliciumkristal van nature energie absorbeert om een elektron van de ene band naar de andere te laten springen, produceerde de simulatie een scherpe afkapwaarde in de gegenereerde harmonischen. Omdat hun model echter al deze specifieke elektronische sprong rekening hield, moesten zij die specifieke absorptiekenmerken uit het brekingsindexmodel verwijderen om te voorkomen dat dit dubbel werd geteld. Wanneer zij dit correct deden, produceerde de simulatie een vloeiend en continu spectrum van harmonischen. Ze testten ook of de trillingen van de kristalatomen, bekend als fononen, een rol speelden. Omdat de lichtpuls zo ongelooflijk snel is, ontdekten zij dat voor de korte afstanden die zij simuleerden, de trage beweging van de atomen de harmonische spectrum niet significant veranderde, wat suggereert dat de elektronische respons te snel plaatsvindt voor de atomen om te reageren.

Een van de meest praktische inzichten van dit werk is de vergelijking tussen twee verschillende manieren om de reis van het licht te simuleren. De onderzoekers vergeleken hun volledige, tijdsgebonden simulatie met een snellere, eenrichtingsbenadering die vaak in het vakgebied wordt gebruikt. Ze ontdekten dat voor zeer dunne plakjes silicium de twee methoden goed overeenkwamen wat betreft de lagere frequentie-harmonischen. Echter, naarm matter de licht verder reisde of naarmatter de frequenties hoger werden, werden de verschillen duidelijk. De eenrichtingsmethode, die negeert dat licht naar achteren reflecteert, slaagde er niet in bepaalde subtiele verschuivingen in het spectrum te vangen die de volledige simulatie wel vastlegde. Dit suggereert dat voor een accurate interpretatie van experimentele gegevens, vooral wanneer men probeert de interne structuur van een materiaal af te leiden uit het licht dat het uitzendt, de volledige, tweerichtingssimulatie noodzakelijk is. De onderzoekers benadrukken dat het licht dat aan de achterkant van het kristal wordt gegenereerd niet simpelweg een som is van wat er aan de voorkant is gegenereerd; de puls is door zijn reis van vorm veranderd, en die vormverandering draagt informatie over de eigenschappen van het materiaal.

De implicaties van dit werk strekken zich uit voorbij alleen silicium. Het raamwerk dat het team heeft gebouwd, is ontworpen om flexibel te zijn, waardoor wetenschappers verschillende materialen en verschillende lasercondities kunnen invoeren om te zien hoe het licht zich gedraagt. Dit is cruciaal voor het gebied van de attoseconde-wetenschap, die tot doel heeft de hoogfrequente lichtpulsen te gebruiken om snapshots te maken van bewegende elektronen binnen atomen en moleculen. Als wetenschappers niet in staat zijn om het signaal gegenereerd door de elektronische structuur van het materiaal te onderscheiden van het signaal dat is veranderd door de reis van het licht door het materiaal, zou hun interpretatie van het experiment gebrekkig kunnen zijn. Door een instrument te bieden dat deze effecten scheidt, hebben de onderzoekers een helderder pad geboden naar het begrijpen van de microscopische wereld. Hun simulaties tonen aan dat de macroscopische reis van de lichtpuls onlosmakelijk verbonden is met de microscopische dans van de elektronen, en dat men om de één echt te begrijpen, de ander moet modelleren. Dit werk beschrijft niet alleen een nieuw fenomeen; het biedt de computationele lens die nodig is om de ware aard van de interactie tussen licht en materie in het extreme te zien, wat de weg vrijmaakt voor preciezere metingen en een dieper begrip van hoe vaste stoffen reageren op de meest intense lichtkrachten.