⚛️ quantum physics

Josephson energy of superconducting junctions: amorphous versus crystalline tunnel barriers

Utilizzando la modellazione NEGF-DFT basata su primi principi, questo studio dimostra che l'energia di Josephson nelle giunzioni superconduttrici è esponenzialmente sensibile alla struttura atomica della barriera di tunnel, rivelando che le barriere di Al $_2$ O $_3$ amorfo esibiscono valori di energia significativamente più elevati e variabili rispetto a quelle cristalline a causa dell'inomogeneità stechiometrica che crea percorsi di tunneling di tipo percolativo.

Autori originali: Wanting Zhang, Aldilene Saraiva-Souza, Félix Beaudoin, Xianghua Kong, Hong Guo, Yu Zhu

Pubblicato 2026-09-10

📖 5 min di lettura🧠 Approfondimento

CC BY 4.0

Autori originali: Wanting Zhang, Aldilene Saraiva-Souza, Félix Beaudoin, Xianghua Kong, Hong Guo, Yu Zhu

Articolo originale sotto licenza CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Questa è una spiegazione generata dall'IA dell'articolo qui sotto. Non è stata scritta né approvata dagli autori. Per precisione tecnica, consulta l'articolo originale. Leggi il disclaimer completo

Nel mondo silenzioso dei circuiti superconduttori, dove l'elettricità scorre senza resistenza, minuscoli componenti chiamati giunzioni Josephson agiscono come gli interruttori e le valvole essenziali. Queste giunzioni sono il cuore dei computer quantistici più avanzati, specificamente di un tipo di qubit noto come transmon. Il comportamento di questi qubit è governato da un valore energetico specifico, che determina quanto velocemente il computer può pensare e quanto chiaramente può distinguere tra diversi stati di informazione. Questa energia dipende interamente dalla facilità con cui gli elettroni possono attraversare per effetto tunnel una barriera di ossido di alluminio incredibilmente sottile, uno strato spesso solo uno o due nanometri. Poiché questa barriera è così sottile, la sua disposizione atomica — se gli atomi siano allineati in una griglia perfetta o disposti in modo disordinato, come in un vetro — ha un impatto enorme su come si muovono gli elettroni. Comprendere esattamente come questo disordine microscopico influenzi l'energia del circuito è fondamentale per costruire macchine quantistiche affidabili, eppure il collegamento tra la struttura atomica caotica della barriera e l'energia macroscopica del dispositivo è rimasto difficile da definire.

Un team di ricercatori ha ora mappato questa connessione costruendo modelli informatici dettagliati di queste giunzioni, confrontando una barriera cristallina perfettamente ordinata con dieci versioni di una barriera amorfa e disordinata. Utilizzando un metodo potente che combina le leggi della meccanica quantistica con la disposizione specifica di ogni atomo, hanno simulato il passaggio degli elettroni attraverso queste strutture. Lo studio si è concentrato su barriere di ossido di alluminio della stessa dimensione, ma con architetture interne differenti. Nel modello cristallino, gli atomi formano un motivo regolare e ripetitivo, mentre nei modelli amorfi, gli atomi sono stati generati da un processo che imita la fusione e il rapido raffreddamento, risultando in un disordine casuale e congelato. I ricercatori hanno poi calcolato la trasmissione degli elettroni attraverso queste barriere per determinare la risultante energia di Josephson, un parametro chiave che stabilisce la frequenza operativa del qubit.

I risultati hanno rivelato una differenza sorprendente tra i due tipi di barriere. Per la singola giunzione cristallina, l'energia calcolata era un valore costante di 0,73 gigahertz. Tuttavia, le dieci giunzioni amorfe hanno mostrato una vasta dispersione di valori. Sebbene l'energia media per il gruppo amorfo fosse di 2,78 gigahertz, i singoli risultati variavano drasticamente, spaziando da un minimo di 0,18 gigahertz a un massimo di 15,1 gigahertz. Ciò significa che i campioni disordinati non hanno prodotto solo un risultato leggermente diverso; hanno prodotto esiti che coprivano quasi due ordini di grandezza. La distribuzione era fortemente asimmetrica, con la maggior parte dei campioni che si raggruppava vicino al valore cristallino, ma con alcuni valori anomali ad altissima energia che alzavano la media. Questa variabilità suggerisce che la disposizione casuale degli atomi nella barriera amorfa crea un paesaggio in cui il percorso per gli elettroni non è uniforme, ma dipende invece dalla specifica disposizione accidentale degli atomi in ogni campione.

Per capire perché ciò accada, i ricercatori hanno esaminato da vicino il movimento degli elettroni attraverso le barriere. Hanno scoperto che il meccanismo di trasporto è l'effetto tunnel quantistico, in cui gli elettroni passano attraverso la barriera anche se non possiedono l'energia necessaria per scavalcarla. Nei campioni amorfi, la composizione chimica dell'ossido non è perfettamente uniforme; alcune regioni sono più ricche di alluminio, mentre altre sono più ricche di ossigeno. Le aree ricche di alluminio agiscono come barriere più basse e facili da attraversare, mentre le aree ricche di ossigeno sono più alte e difficili da superare. In alcuni campioni amorfi, queste regioni a bassa barriera si sono allineate casualmente, creando un percorso continuo, di tipo percolativo, che ha permesso agli elettroni di fluire molto più facilmente rispetto a quanto farebbero attraverso una barriera uniforme. In altri campioni, questi percorsi facili non si sono connessi, costringendo gli elettroni a faticare attraverso barriere più alte. Questo spiega perché alcune giunzioni amorfe avevano valori di energia molto più alti rispetto al riferimento cristallino, mentre altre erano inferiori.

Lo studio ha anche affrontato il modo in cui queste scoperte microscopiche si relazionano ai dispositivi reali utilizzati nei laboratori. I modelli informatici rappresentavano minuscole porzioni di materiale, di circa 1,44 per 1,25 nanometri. Una vera giunzione Josephson utilizzata in un computer quantistico è molto più grande, misurando tipicamente 200 per 200 nanometri. Ciò significa che un dispositivo reale contiene più di 20.000 di queste regioni microscopiche. Poiché la giunzione reale è così grande, essa media efficacemente le variazioni estreme osservate nei piccoli modelli. L'energia finale di un dispositivo reale è probabilmente una combinazione di molti diversi percorsi microscopici, che smorza le oscillazioni selvagge viste nelle singole simulazioni. Questo suggerisce che, sebbene il disordine atomico crei una variabilità significativa alla scala più piccola, il dispositivo macroscopico beneficia di un effetto di auto-mediazione che ne stabilizza le prestazioni.

I ricercatori hanno concluso che la variabilità nell'energia di Josephson è una conseguenza diretta dell'inomogeneità stechiometrica — ovvero l'unione irregolare di alluminio e ossigeno — all'interno della barriera amorfa. La presenza di regioni ricche di alluminio e a bassa barriera, che possono connettersi attraverso l'ossido, crea percorsi che potenziano fortemente la conducibilità, portando ai valori anomali elevati. Al contrario, i campioni privi di questi percorsi connessi possono essere ancora meno trasparenti rispetto al cristallo perfetto. Sebbene lo studio si sia basato su simulazioni e che i metodi matematici specifici utilizzati tendano a sottostimare leggermente il gap energetico dell'ossido, le tendenze osservate forniscono una via quantitativa chiara dal microstruttura dell'ossido alla scala energetica del circuito superconduttore. Questo lavoro stabilisce che la struttura atomica casuale della barriera non è solo un dettaglio minore, ma un fattore fondamentale che detta il panorama energetico del computer quantistico, offrendo un nuovo modo di pensare a come questi dispositivi vengano costruiti e a come le loro prestazioni possano essere ottimizzate.

Sintesi Tecnica: Energia di Josephson delle Giunzioni Superconduttive: Barriere Tunnel Amorfhe rispetto a Cristalline

Definizione del Problema
L'energia di Josephson ( $E_J$ ) è un parametro fondamentale che determina la frequenza e l'anarmonicità dei qubit superconduttori transmon. Nelle giunzioni Al/AlO $_x$ /Al, $E_J$ è governata dal tunneling elettronico attraverso una barriera di ossido ultrasottile, esibendo una dipendenza esponenziale dalla struttura atomica della barriera stessa. Mentre i modelli di circuiti macroscopici trattano la giunzione come un elemento concentrato, la barriera reale è strutturalmente e chimicamente non uniforme, contenendo variazioni di spessore, carenze di ossigeno e disomogeneità stechiometriche. Gli esistenti modelli atomistici si sono concentrati principalmente su barriere cristalline o terminazioni d'interfaccia specifiche. Tuttavia, le realistiche barriere di AlO $_x$ sono amorfe, introducendo variazioni stocastiche nella configurazione atomica che una singola struttura non può rappresentare. La sfida consiste nell'stabilire una via quantitativa e atomistica che colleghi la microstruttura delle barriere di ossido amorfe alla scala di energia macroscopica $E_J$ .

Metodologia
Gli autori impiegano un flusso di lavoro di modellazione di dispositivi basato sui primi principi, combinando la Teoria del Funzionale della Densità (DFT) con il formalismo della Funzione di Green fuori dall'equilibrio (NEGF) per simulare il trasporto quantistico attraverso giunzioni Al/Al $_2$ O $_3$ /Al.

Generazione della Struttura Atomica:
- Riferimento Cristallino: È stata costruita una supercella commisurata Al(111)– $\alpha$ -Al $_2$ O $_3$ –Al(111) con 420 atomi. L'ossido è stato terminato in modo asimmetrico (Al su un lato, O sull'altro).
- Insieme Amorfo: Dieci barriere indipendenti di Al $_2$ O $_3$ amorfo sono state generate utilizzando un protocollo di dinamica molecolare melt–quench in LAMMPS (utilizzando il campo di forza Matsui). Queste sono state successivamente inserite tra lastre cristalline di Al(111) e rilassate utilizzando il potenziale interatomico di apprendimento automatico CHGNet per gestire le eterogenee interfacce metallo-ossido. Ogni giunzione amorfa contiene 1080 atomi.
- Geometria: I modelli cristallini e amorfi possiedono diverse periodicità laterali e aree di sezione trasversale ( $A_0$ ), che sono esplicitamente considerate nell'analisi del trasporto.
Calcolo del Trasporto Quantistico:
- I calcoli di trasporto sono stati eseguiti utilizzando il pacchetto NanoDCAL con una base di combinazione lineare di orbitali atomici (LCAO) e pseudopotenziali norm-conservanti (livello GGA).
- Il coefficiente di trasmissione dipendente dall'energia $T(E)$ è stato computato per le giunzioni rilassate.
- La conduttanza nello stato normale ( $G_N$ ) è stata derivata dalla trasmissione al livello di Fermi $T(E_F)$ tramite la formula di Landauer, scalata per l'area della cella di simulazione $A_0$ e l'area del dispositivo $A$ (impostata a $200 \times 200$ nm $^2$ per il confronto).
- La corrente critica ( $I_c$ ) e l'energia di Josephson ( $E_J$ ) sono state estratte utilizzando la relazione di Ambegaokar–Baratoff a temperatura zero e la definizione $E_J = \hbar I_c / 2e$ , assumendo un gap superconduttivo $\Delta = 0.18$ meV.
Analisi degli Stati di Scattering:
- Per comprendere il meccanismo di trasporto, gli autori hanno analizzato la probabilità degli stati di scattering al livello di Fermi, cercando specificamente schemi di decadimento evanescente e localizzazione della funzione d'onda nello spazio reale.

Risultati Chiave

Variabilità in $E_J$ : Le dieci giunzioni amorfe esibiscono una vasta distribuzione di energie di Josephson. Il valore medio di $E_J/h$ è 2.78 GHz con una deviazione standard di 4.67 GHz. Al contrario, il riferimento cristallino produce 0.73 GHz.
Distribuzione Asimmetrica: La distribuzione è fortemente asimmetrica. Il valore mediano (0.86 GHz) è vicino al riferimento cristallino, ma i singoli valori spaziano per quasi due ordini di grandezza, da 0.18 GHz (Campione 5) a 15.1 GHz (Campione 6).
Caratteristiche di Trasmissione: Le curve di trasmissione normalizzate per l'area ( $T(E)/A_0$ ) mostrano una struttura di gap isolante simile per tutti i campioni. Tuttavia, la posizione del livello di Fermi rispetto al minimo di trasmissione varia significativamente tra i campioni. Piccoli spostamenti nella posizione del minimo rispetto a $E_F$ causano grandi variazioni in $T(E_F)$ a causa della natura esponenziale del tunneling.
Meccanismo di Trasporto: L'analisi degli stati di scattering conferma che il trasporto avviene tramite tunneling quantistico. La variabilità deriva dall'inomogeneità stechiometrica all'interno dell'ossido amorfo.
- Percorsi di Percolazione: Nei campioni ad alta $E_J$ , regioni ricche di Al (che hanno un'altezza di barriera effettiva inferiore) si connettono attraverso l'ossido, formando percorsi di tunneling di tipo percolativo. Ciò consente agli elettroni di attraversare la barriera più facilmente, aumentando significativamente la conduttanza.
- Localizzazione: Nei campioni a bassa $E_J$ , tali percorsi connessi a bassa barriera sono assenti, e gli stati di scattering rimangono localizzati vicino all'elettrodo incidente, risultando in una trasparenza inferiore rispetto alla giunzione stessa, anche rispetto a quella cristallina.

Significato e Rivendicazioni
Il lavoro stabilisce un collegamento quantitativo tra la struttura atomica microscopica delle barriere di ossido amorfo e l'energia di Josephson macroscopica dei circuiti superconduttori. Gli autori affermano che la grande varianza in $E_J$ osservata nelle giunzioni amorfe non è rumore casuale, ma una conseguenza diretta della formazione stocastica di canali di tunneling di tipo percolativo guidati dalle fluttuazioni locali della stechiometria.

Fondamentalmente, gli autori osservano che, sebbene i singoli modelli atomistici mostrino un'estrema variabilità, una giunzione macroscopica realistica (ad esempio, $200 \times 200$ nm $^2$ ) comprende oltre $2 \times 10^4$ di queste regioni microscopiche. Di conseguenza, l'energia $E_J$ a livello di dispositivo dovrebbe "auto-mediarsi" su queste configurazioni, riflettendo una media dell'insieme piuttosto che gli estremi outlier visti nelle piccole celle di simulazione.

Limitazioni e Direzioni Future
Gli autori riconoscono che il funzionale GGA utilizzato probabilmente sottostima il band gap dell'ossido, portando potenzialmente a valori di trasmissione assoluti leggermente sovrastimati, sebbene gli andamenti tra i campioni dovrebbero rimanere robusti. Suggeriscono che lavori futuri dovrebbero impiegare funzionali GGA+U o ibridi (come HSE) per una maggiore accuratezza quantitativa. Inoltre, propongono che il calcolo di un insieme più ampio di campioni amorfi permetterebbe di ottenere una distribuzione statistica appropriata di $E_J$ da confrontare direttamente con la dispersione sperimentale tra le giunzioni.

Articoli simili