Device variability of Josephson junctions induced by interface roughness

Technische Samenvatting: Variabiliteit van Josephson-contacten door Interface-ruwheid

Probleemstelling
Naarmate supergeleidende quantumprocessors opschalen naar de duizenden qubits die vereist zijn voor fouttolerante quantumcomputatie (bijv. Shor's algoritme), is de variabiliteit tussen apparaten een kritieke flessenhals geworden. Supergeleidende qubits, met name transmonen, vertrouwen op Al/AlOx/Al Josephson-contacten die opereren in het tunnelregime. In dit regime hangt de kritische stroom $I_c$ — en daarmee de Josephson-energie $E_J$ — exponentieel af van microscopische details van de barrière. Zelfs kleine variaties in de geometrie van het Al/AlOx-interface, zoals lokale diktefluctuaties of contaminatie, kunnen substantiële verschuivingen in de transitiefrequenties van de qubit en de gate-fidelities induceren. Hoewel uitgebreide experimentele inspanningen gericht zijn op het verbeteren van de reproduceerbaarheid van contacten door procesoptimalisatie en microscopie, ontbrak een kwantitatief theoretisch kader dat microstructuur-wanorde (specifiek interface-ruwheid) verbindt met de energie-statistieken op contactniveau.

Methodologie
De auteurs ontwikkelen een kwantitatief modelleringskader dat de brug slaat tussen de interface-wanorde op atomaire schaal en macroscopische contacttransporteigenschappen. De aanpak integreert drie primaire componenten:

1. Modellering van Interface-ruwheid: Het Al/AlOx-interface wordt gemodelleerd als een Gaussisch willekeurig veld $h(x,y)$ gekenmerkt door twee parameters: de wortelgemiddelde kwadratische (RMS) ruwheidsamplitude $\sigma$ en de transversale correlatielengte $\xi$. Deze parameters zijn afgeleid van bestaande experimentele literatuur (bijv. TEM-studies van korrelgroottes en diktehistogrammen) en nieuwe moleculaire dynamica (MD) simulaties. De MD-simulaties, gebruikmakend van de CHGNet machine learning potentiaal, modelleren de oxidatie van een Al-oppervlak om $\sigma$ onder geïdealiseerde omstandigheden te schatten.
2. Kwantumtransport-benadering: Om de computationele onhandelbaarheid van het oplossen van de Bogoliubov–de Gennes (BdG) vergelijking voor realistische contactafmetingen ($\sim 100 \text{ nm} \times 100 \text{ nm}$) aan te pakken, gebruiken de auteurs twee belangrijke benaderingen:
   • Ambegaokar–Baratoff (AB) Relatie: De kritische stroom wordt berekend via de AB-relatie, die $I_c$ koppelt aan de normale weerstand $R_N$. Dit vermindert de computationele complexiteit door het integreren over het volledige energiespectrum te vermijden en zich te concentreren op de transmissiecoëfficiënt bij de Fermi-energie.
   • Lokale Dikte-benadering: Gezien de Fermi-golflengte ($\lambda_F \approx 3.6$ Å) en de barrière-vervallengte ($\lambda_D \approx 1.9$ Å) veel kleiner zijn dan de correlatielengte van de ruwheid ($\xi \sim 10$ nm), wordt het contact behandeld als een ensemble van onafhankelijke tunnelkanalen. De totale conductantie wordt benaderd als het integraal van lokale conductantie-dichtheden, waarbij de lokale barrière-dikte varieert volgens het ruwheidsprofiel.
3. Numerieke Simulatie: Het model wordt toegepast op een ensemble van 5.000 Josephson-contacten met variërende tunnelbarrière-profielen gegenereerd door het Gaussische willekeurige veldmodel.

Belangrijkste Resultaten
De numerieke simulaties leveren verschillende kwantitatieve inzichten op over het statistische gedrag van $E_J$:

• Log-normale Verdeling: De distributie van de Josephson-energie $E_J$ over het ensemble is niet Gaussisch maar volgt een log-normale verdeling. Dit komt doordat $E_J$ exponentieel afhankelijk is van de barrière-dikte, welke normaal verdeeld is. Bijgevolg is de waarschijnlijkheidsdichtheidsfunctie scheef naar grotere waarden.
• Impact van RMS-ruwheid ($\sigma$): Het vergroten van de RMS-ruwheidsamplitude $\sigma$ leidt tot een significante toename van zowel de gemiddelde waarde $\langle E_J \rangle$ als de variantie $\text{Var}(E_J)$. De toename in het gemiddelde wordt toegeschreven aan de exponentiële gevoeligheid van tunneling, waarbij zeldzame, anomalistisch dunne regio's de totale conductantie domineren.
• Impact van Correlatielengte ($\xi$): Het vergroten van de transversale correlatielengte $\xi$ veroorzaakt een lichte afname in het gemiddelde $\langle E_J \rangle$ maar een sterke toename in de variantie $\text{Var}(E_J)$. De variantie neemt toe omdat een grotere $\xi$ het aantal statistisch onafhankelijke tunnelkanalen vermindert (zelf-averaging), waardoor de fluctuaties tussen monsters worden versterkt. Het gemiddelde neemt licht af door eindige-omvangseffecten; naarmate $\xi$ de afmetingen van het contact nadert, wordt de kans kleiner om de zeldzame, extreem dunne barrièreregio's te realiseren die de tunneling domineren.
• Variabiliteit van Qubit-frequentie: Gebruikmakend van de afgeleide statistieken voor $E_J$ en uitgaande van een typische opladingsenergie $E_C$, voorspelt het model een standaarddeviatie in de qubit-transitiefrequentie $\omega_{01}/2\pi$ van ongeveer 0,364 GHz rond een gemiddelde van 6,145 GHz voor de specifieke geteste parameters.

Betekenis en Claims
Het artikel beweert een "kwantitatief en intuïtief beeld" te bieden van hoe oppervlakte-ruwheid de variabiliteit in de Josephson-energie induceert. Door de microstructuur-parameters ($\sigma$ en $\xi$) direct te koppelen aan de contact-niveau statistieken, biedt dit werk een kader voor het begrijpen van de opbrengst- en uniformiteitsuitdagingen in grootschalige quantumprocessors. De auteurs benadrukken dat hun resultaten "direct relevant zijn voor contactontwerp", wat suggereert dat het controleren van de interface-ruwheidsparameters cruciaal is voor het minimaliseren van de variabiliteit tussen apparaten.

De auteurs behouden een bescheiden standpunt over de reikwijdte van hun werk. Ze merken op dat hun model ervan uitgaat dat $\sigma$ en $\xi$ onafhankelijk en identiek zijn voor beide interfaces in het contact, een aanname die experimentele verificatie vereist onder diverse groeicondities. Bovendien identificeren zij toekomstige richtingen, zoals het verbinden van dit model met de experimentele karakterisering van gefabriceerde contacten en het integreren van hun supergeleidings-solvers met microwave engineering software om het co-design van contacten en circuitcomponenten mogelijk te maken. Het artikel stelt geen nieuwe experimentele protocollen voor of claimt het variabiliteitsprobleem volledig op te lossen, maar vestigt een theoretische baseline voor het kwantificeren van de impact van interface-ruwheid.