⚛️ quantum physics

Josephson energy of superconducting junctions: amorphous versus crystalline tunnel barriers

Utilizando modelagem NEGF-DFT de primeiros princípios, este estudo demonstra que a energia de Josephson em junções supercondutoras é exponencialmente sensível à estrutura atômica da barreira de tunelamento, revelando que barreiras de Al $_2$ O $_3$ amorfas exibem valores de energia significativamente mais altos e variáveis do que as cristalinas devido à inhomogeneidade estequiométrica criando caminhos de tunelamento do tipo percolação.

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

Publicado 2026-09-10

📖 5 min de leitura🧠 Leitura aprofundada

CC BY 4.0

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

Artigo original sob licença CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Esta é uma explicação gerada por IA do artigo abaixo. Não foi escrita nem endossada pelos autores. Para precisão técnica, consulte o artigo original. Ler aviso legal completo

No mundo silencioso dos circuitos supercondutores, onde a eletricidade flui sem resistência, componentes minúsculos chamados junções de Josephson atuam como os interruptores e válvulas essenciais. Essas junções são o coração dos computadores quânticos mais avançados, especificamente um tipo de qubit conhecido como transmon. O comportamento desses qubits é governado por um valor de energia específico, que determina o quão rápido o computador pode pensar e o quão claramente ele pode distinguir entre diferentes estados de informação. Essa energia depende inteiramente de quão facilmente os elétrons podem tunelar através de uma barreira de óxido de alumínio incrivelmente fina, uma camada de apenas um ou dois nanômetros de espessura. Como essa barreira é tão fina, seu arranjo atômico — se os átomos estão alinhados em uma grade perfeita ou desordenados em uma estrutura vítrea e aleatória — tem um impacto massivo em como os elétrons se movem. Compreender exatamente como esse desordem microscópica afeta a energia do circuito é crucial para construir máquinas quânticas confiáveis, mas a conexão entre a estrutura atômica caótica da barreira e a energia macroscópica do dispositivo permaneceu difícil de definir.

Uma equipe de pesquisadores mapeou agora essa conexão ao construir modelos computacionais detalhados dessas junções, comparando uma barreira cristalina perfeitamente ordenada contra dez versões de uma barreira amorfa e desordenada. Usando um método poderoso que combina as leis da mecânica quântica com o arranjo específico de cada átomo, eles simularam como os elétrons viajam através dessas estruturas. O estudo focou em barreiras de óxido de alumínio com a mesma espessura, mas com arquiteturas internas diferentes. No modelo cristalino, os átomos formam um padrão repetitivo e organizado, enquanto nos modelos amorfos, os átomos foram gerados por um processo que imita o derretimento e o resfriamento rápido, resultando em uma desordem aleatória e congelada. Os pesquisadores então calcularam a transmissão de elétrons através dessas barreiras para determinar a energia de Josephson resultante, um parâmetro fundamental que define a frequência de operação do qubit.

Os resultados revelaram uma diferença marcante entre os dois tipos de barreiras. Para a junção de cristal único, a energia calculada foi de 0,73 gigahertz constantes. No entanto, as dez junções amorfas mostraram uma dispersão selvagem de valores. Embora a energia média para o grupo amorfo fosse de 2,78 gigahertz, os resultados individuais variaram dramaticamente, indo de tão baixo quanto 0,18 gigahertz até tão alto quanto 15,1 gigahertz. Isso significa que as amostras desordenadas não produziram apenas um resultado ligeiramente diferente; elas produziram resultados que abrangem quase duas ordens de magnitude. A distribuição foi fortemente enviesada, com a maioria das amostras se agrupando perto do valor cristalino, mas com alguns valores discrepantes (outliers) de energia extremamente alta puxando a média para cima. Essa variabilidade sugere que o arranjo aleatório de átomos na barreira amorfa cria um cenário onde o caminho para os elétrons não é uniforme, mas depende do arranjo específico e acidental de átomos em cada amostra.

Para entender por que isso acontece, os pesquisadores observaram de perto como os elétrons se moviam através das barreiras. Eles descobriram que o mecanismo de transporte é o tunelamento quântico, onde os elétrons passam pela barreira mesmo sem ter energia para escalá-la. Nas amostras amorfas, a composição química do óxido não é perfeitamente uniforme; algumas regiões são mais ricas em alumínio, enquanto outras são mais ricas em oxigênio. As áreas ricas em alumínio atuam como barreiras mais baixas e fáceis de atravessar pelos elétrons, enquanto as áreas ricas em oxigênio são mais altas e difíceis de atravessar. Em algumas das amostras amorfas, essas regiões de baixa barreira coincidiram, criando um caminho contínuo, do tipo percolação, que permitiu que os elétrons fluíssem muito mais facilmente do que passariam por uma barreira uniforme. Em outras amostras, esses caminhos fáceis não se conectaram, forçando os elétrons a lutar através de barreiras mais altas. Isso explica por que algumas junções amorfas tiveram valores de energia muito superiores à referência cristalina, enquanto outras foram inferiores.

O estudo também abordou como essas descobertas microscópicas se relacionam com os dispositivos reais usados em laboratórios. Os modelos computacionais representavam fatias minúsculas de material, com aproximadamente 1,44 por 1,25 nanômetros de tamanho. Uma junção de Josephson real usada em um computador quântico é muito maior, medindo tipicamente 200 por 200 nanômetros. Isso significa que um dispositivo real contém mais de 20.000 dessas regiões microscópicas. Como a junção real é tão grande, ela efetivamente faz a média das variações extremas vistas nos pequenos modelos. A energia final de um dispositivo real é provavelmente uma mistura de muitos caminhos microscópicos diferentes, suavizando as oscilações selvagens vistas nas simulações individuais. Isso sugere que, embora a desordem atômica crie uma variabilidade significativa na menor escala, o dispositivo macroscópico se beneficia de um efeito de autossuavização que estabiliza seu desempenho.

Os pesquisadores concluíram que a variabilidade na energia de Josephson é uma consequência direta da inhomogeneidade estequiométrica — a mistura desigual de alumínio e oxigênio — dentro da barreira amorfa. A presença de regiões ricas em alumínio e de baixa barreira que podem se conectar através do óxido cria caminhos que aumentam fortemente a condutância, levando aos valores discrepantes de alta energia. Por outro lado, amostras sem esses caminhos conectados podem ser ainda menos transparentes do que o cristal perfeito. Embora o estudo tenha dependido de simulações e os métodos matemáticos específicos utilizados tendam a subestimar ligeiramente o gap de energia do óxido, as tendências observadas fornecem uma rota quantitativa clara da estrutura microscópica do óxido para a escala de energia do circuito supercondutor. Este trabalho estabelece que a estrutura atômica aleatória da barreira não é apenas um detalhe menor, mas um fator fundamental que dita o cenário de energia do computador quântico, oferecendo uma nova maneira de pensar sobre como esses dispositivos são construídos e como seu desempenho pode ser otimizado.

Resumo Técnico: Energia de Josephson de Junções Supercondutoras: Barreiras de Tunelamento Amorfas versus Cristalinas

Definição do Problema
A energia de Josephson ( $E_J$ ) é um parâmetro fundamental que determina a frequência e a anharmonicidade de qubits supercondutores do tipo transmon. Em junções de Al/AlO $_x$ /Al, a $E_J$ é governada pelo tunelamento de elétrons através de uma barreira de óxido ultrafina, exibindo uma dependência exponencial da estrutura atômica da barreira. Enquanto os modelos de circuitos macroscópicos tratam a junção como um elemento concentrado, a barreira real é estrutural e quimicamente não uniforme, contendo variações de espessura, deficiências de oxigênio e inhomogeneidades estequiométricas. Os modelos atomísticos existentes focaram amplamente em barreiras cristalinas ou terminações de interface específicas. No entanto, as barreiras de AlO $_x$ realistas são amorfas, introduzindo variações estocásticas na configuração atômica que uma única estrutura não pode representar. O desafio reside em estabelecer uma rota quantitativa e atomística da microestrutura das barreiras de óxido amorfas para a escala de energia macroscópica $E_J$ .

Metodologia
Os autores empregam um fluxo de trabalho de modelagem de dispositivos de primeira ordem combinando a Teoria do Funcional da Densidade (DFT) com o formalismo da Função de Green de Não-Equilíbrio (NEGF) para simular o transporte quântico através de junções Al/Al $_2$ O $_3$ /Al.

Geração de Estrutura Atômica:
- Referência Cristalina: Uma supercélula comensurável de Al(111)– $\alpha$ -Al $_2$ O $_3$ –Al(111) foi construída com 420 átomos. O óxido foi terminado assimetricamente (Al de um lado, O do outro).
- Ensemble Amorfo: Dez barreiras de Al $_2$ O $_3$ amorfas independentes foram geradas usando um protocolo de dinâmica molecular de fusão-resfriamento (melt–quench) no LAMMPS (usando o campo de força de Matsui). Estas foram subsequentemente inseridas entre placas de Al(111) cristalinas e relaxadas usando o potencial interatômico de aprendizado de máquina CHGNet para lidar com as interfaces metal–óxido heterogêneas. Cada junção amorfa contém 1080 átomos.
- Geometria: Os modelos cristalinos e amorfos possuem periodicidades laterais e áreas de seção transversal ( $A_0$ ) diferentes, que são explicitamente contabilizadas na análise de transporte.
Cálculo de Transporte Quântico:
- Os cálculos de transporte foram realizados utilizando o pacote NanoDCAL com uma base de combinação linear de orbitais atômicos (LCAO) e pseudopotenciais de conservação de norma (nível GGA).
- O coeficiente de transmissão dependente da energia $T(E)$ foi computado para as junções relaxadas.
- A condutância no estado normal ( $G_N$ ) foi derivada da transmissão no nível de Fermi $T(E_F)$ via fórmula de Landauer, escalonada pela área da célula de simulação $A_0$ e pela área do dispositivo $A$ (definida como $200 \times 200$ nm $^2$ para comparação).
- A corrente crítica ( $I_c$ ) e a energia de Josephson ( $E_J$ ) foram extraídas usando a relação de Ambegaokar–Baratoff a zero de temperatura e a definição $E_J = \hbar I_c / 2e$ , assumindo um gap supercondutor $\Delta = 0,18$ meV.
Análise de Estado de Espalhamento:
- Para compreender o mecanismo de transporte, os autores analisaram a probabilidade de estados de espalhamento no nível de Fermi, especificamente procurando por padrões de decaimento evanescente e localização da função de onda no espaço real.

Principais Resultados

Variabilidade em $E_J$ : As dez junções amorfas exibem uma ampla distribuição de energias de Josephson. A média de $E_J/h$ é 2,78 GHz com um desvio padrão de 4,67 GHz. Em contraste, a referência cristalina produz 0,73 GHz.
Distribuição Assimétrica: A distribuição é fortemente assimétrica. O valor da mediana (0,86 GHz) é próximo à referência cristalina, mas os valores individuais abrangem quase duas ordens de magnitude, variando de 0,18 GHz (Amostra 5) a 15,1 GHz (Amostra 6).
Características de Transmissão: As curvas de transmissão normalizadas pela área ( $T(E)/A_0$ ) mostram uma estrutura de gap isolante semelhante para todas as amostras. No entanto, a posição do nível de Fermi em relação ao vale de transmissão varia significativamente entre as amostras. Pequenos deslocamentos na posição do vale em relação a $E_F$ resultam em grandes variações em $T(E_F)$ devido à natureza exponencial do tunelamento.
Mecanismo de Transporte: A análise de estado de espalhamento confirma que o transporte ocorre via tunelamento quântico. A variabilidade decorre da inhomogeneidade estequiométrica dentro do óxido amorfo.
- Caminhos de Percolação: Em amostras de alta $E_J$ , regiões ricas em Al (que possuem uma altura de barreira efetiva menor) conectam-se através do óxido, formando caminhos de tunelamento do tipo percolação. Isso permite que os elétrons atravessem a barreira mais facilmente, aumentando significativamente a condutância.
- Localização: Em amostras de baixa $E_J$ , tais caminhos conectados de baixa barreira estão ausentes, e os estados de espalhamento permanecem localizados próximos ao eletrodo incidente, resultando em uma transparência inferior até mesmo à da junção cristalina.

Significância e Alegações
O artigo estabelece um elo quantitativo entre a estrutura atômica microscópica das barreiras de óxido amorfo e a energia de Josephson macroscópica de circuitos supercondutores. Os autores alegam que a grande variância observada em $E_J$ nas junções amorfas não é ruído aleatório, mas uma consequência direta da formação estocástica de canais de tunelamento do tipo percolação impulsionados por flutuações estequiométricas locais.

Crucialmente, os autores observam que, embora modelos atomísticos individuais mostrem variabilidade extrema, uma junção macroscópica realista (por exemplo, $200 \times 200$ nm $^2$ ) engloba mais de $2 \times 10^4$ dessas regiões microscópicas. Consequentemente, espera-se que a $E_J$ do dispositivo realize uma "automedição" (self-average) sobre essas configurações, refletindo uma média de ensemble em vez dos valores extremos vistos nas pequenas células de simulação.

Limitações e Direções Futuras
Os autores reconhecem que o funcional GGA utilizado provavelmente subestima o gap de banda do óxido, o que pode levar a valores de transmissão absoluta ligeiramente superestimados, embora as tendências entre as amostras devam permanecer robustas. Eles sugerem que trabalhos futuros utilizem funcionais GGA+U ou híbridos (por exemplo, HSE) para maior precisão quantitativa. Adicionalmente, propõem que o cálculo de um ensemble maior de amostras amorfas permitiria uma distribuição estatística adequada de $E_J$ para ser comparada diretamente com as dispersões de junção para junção observadas experimentalmente.

Mais como este