Theoretical study on ambient pressure superconductivity in La$_3$Ni$_2$O$_7$ thin films : structural analysis, model construction, and robustness of $s\pm$-wave pairing

Este estudo teórico demonstra que, embora a estrutura eletrônica de filmes finos de La$_3$Ni$_2$O$_7$ sob pressão ambiente varie com a estrutura cristalina e os detalhes computacionais, a robustez do emparelhamento do tipo $s\pm$ mediado por flutuações de spin de energia finita persiste, embora a redução observada na temperatura crítica em comparação com o volume sob pressão seja melhor explicada por modelos com pequeno hopping intercamada derivados de estruturas de rede experimentais.

O essencial

A Visão Geral: Um Supercondutor à Pressão Ambiente

Imagine um material chamado La₃Ni₂O₇ (um tipo de cristal à base de níquel) que pode conduzir eletricidade com resistência zero (supercondutividade). Cientistas descobriram recentemente que, se você espremer esse material com pressão massiva, ele se torna um supercondutor a cerca de 80 Kelvin (muito frio, mas quente para supercondutores).

Recentemente, pesquisadores encontraram uma maneira de fazer isso acontecer sem espremê-lo. Eles cresceram o material como um filme muito fino sobre um tipo específico de "chão" de cristal (um substrato). O chão era ligeiramente menor que o filme, o que espremia o filme pelos lados, imitando o efeito de alta pressão. Esse filme tornou-se supercondutor a cerca de 40 Kelvin.

A Pergunta: Por que o filme fino funciona a uma temperatura mais baixa (40 K) do que o material maciço espremido (80 K)? E qual é a "receita" exata dentro do material que faz a eletricidade fluir sem resistência?

A Abordagem dos Cientistas: Construindo um Modelo Digital

Os autores deste artigo não apenas adivinharam; eles construíram uma simulação computadorizada detalhada. Pense nisso como um motor de videogame onde eles tentaram recriar a física desse material do zero.

1. O Projeto (Estrutura): Eles olharam para o "projeto" dos átomos. Eles tentaram dois projetos diferentes:
   • O Projeto Teórico: O que os cálculos do computador deles diziam que os átomos deveriam parecer.
   • O Projeto Experimental: O que os cientistas realmente mediram no laboratório recentemente.
2. O Motor (FLEX): Eles usaram um motor matemático complexo chamado FLEX (Aproximação de Troca de Flutuação). Imagine esse motor como um simulador de clima superpreciso. Em vez de prever chuva, ele prevê como os elétrons (as partículas minúsculas que carregam eletricidade) dançam e interagem entre si. Ele leva em conta cada movimento possível que os elétrons podem fazer, não apenas os óbvios.

Principais Descobertas: A "Dança" dos Elétrons

1. O Mistério do "Bolsão γ"

No mundo desses materiais, há uma forma específica da multidão de elétrons chamada bolsão γ. Alguns cientistas pensavam que esse bolsão era essencial para a supercondutividade; outros pensavam que não importava.

• O Veredito do Artigo: Os autores descobriram que se esse "bolsão γ" existe ou desaparece depende inteiramente de qual projeto você usa (teórico vs. experimental) e de como você ajusta a matemática.
• A Analogia: É como olhar para uma multidão através de óculos de cores diferentes. Através de uma lente, você vê um grupo distinto de pessoas (o bolsão); através de outra, eles se misturam.
• O Resultado: Surpreendentemente, não importava. Se o bolsão estava lá ou não, a supercondutividade permanecia forte. A "dança" dos elétrons era robusta o suficiente para lidar com essas mudanças estruturais.

2. A "Cola" que Mantém Tudo Juntos

Como os elétrons se emparelham para conduzir eletricidade? Geralmente, eles precisam de uma "cola".

• O Veredito do Artigo: A cola aqui são as flutuações de spin. Imagine que os elétrons são dançarinos girando. Às vezes, eles oscilam ou flutuam em seu spin. Essas oscilações atuam como uma batida rítmica que ajuda os dançarinos a se emparelhar.
• A Reviravolta: O artigo argumenta que essa "batida" vem de oscilações de alta energia, não apenas dos movimentos lentos e óbvios perto da superfície da multidão de elétrons. Como a cola é baseada nessas oscilações de alta energia, a supercondutividade é muito estável e não se quebra facilmente se a forma da multidão de elétrons mudar ligeiramente.

3. Por que o Filme é Mais Frio (40 K) que o Maciço (80 K)?

Este foi o maior enigma. O filme fino supercondut a metade da temperatura do material maciço pressurizado.

• O Veredito do Artigo: A diferença se resume a um número específico: $|t_\perp|$.
• A Analogia: Pense no material como um prédio de dois andares onde os elétrons podem pular entre os andares.
  • No maciço pressurizado, os andares estão perfeitamente alinhados, e o salto entre eles é forte e fácil (Alto $|t_\perp|$). Isso cria uma pista de dança muito eficiente (80 K).
  • No filme fino, as medições experimentais mostram que os andares estão ligeiramente desalinhados ou o salto é mais fraco (Baixo $|t_\perp|$).
• A Conclusão: Quando os autores usaram o "Projeto Experimental" (que mostrava esse salto mais fraco) em sua simulação, a temperatura supercondutora caiu exatamente para os 40 K observados. Quando usaram o "Projeto Teórico" (que previa um salto mais forte), a temperatura permaneceu alta em 80 K.
• A Lição: A razão pela qual o filme é "mais fraco" é provavelmente porque a estrutura física real do filme tem uma conexão mais fraca entre suas camadas do que a teoria previa.

Resumo em Poucas Palavras

Os cientistas construíram uma simulação de alta tecnologia para entender por que um novo filme supercondutor funciona à pressão ambiente. Eles descobriram que:

1. O mecanismo de emparelhamento é resistente: Os elétrons se emparelham usando "oscilações" de alta energia (flutuações de spin), tornando a supercondutividade muito robusta contra pequenas mudanças na forma do material.
2. O "bolsão γ" não importa: Se uma forma específica de elétron existe ou não não muda o resultado.
3. A queda de temperatura é estrutural: O filme atinge apenas 40 K (em vez de 80 K) porque a distância física real entre as camadas atômicas no filme é ligeiramente diferente do que a teoria previa, tornando o "salto" entre as camadas mais fraco.

O artigo essencialmente diz: "Nós conhecemos a receita para a supercondutividade neste material. A razão pela qual o filme é ligeiramente menos eficiente que o bloco pressurizado é simplesmente porque as camadas do filme não estão tão perfeitamente conectadas quanto pensávamos que estariam."

Resumo técnico

1. Formulação do Problema

A descoberta de supercondutividade de alta temperatura ($T_c \approx 80$ K) no nickelato de bicamada sob pressão $La_3Ni_2O_7$ despertou intenso interesse. Recentemente, supercondutividade em pressão ambiente ($T_c \approx 40$ K) foi observada em filmes finos de $La_3Ni_2O_7$ crescidos em substratos como $SrLaAlO_4$ (SLAO). No entanto, um desafio teórico significativo permanece:

• A Discrepância de $T_c$: A $T_c$ em filmes finos é aproximadamente metade daquela do sistema bulk sob pressão.
• Ambiguidade Estrutural: Há uma grande divergência entre as estruturas cristalinas determinadas pela otimização da Teoria do Funcional da Densidade (DFT) e aquelas determinadas experimentalmente. Especificamente, a estrutura experimental mostra um integral de salto intercamada significativamente menor ($|t_\perp| \approx 0,4$ eV) em comparação com as previsões da DFT ($|t_\perp| \approx 0,6$ eV).
• Controvérsia da Superfície de Fermi: Estudos experimentais e teóricos discordam sobre a existência do "$\gamma$-pocket" (uma superfície de Fermi originada do topo da banda de ligação $d_{3z^2-r^2}$). Sua presença ou ausência é debatida quanto ao seu impacto na simetria de emparelhamento.
• Objetivo: Os autores visam reproduzir teoricamente a supercondutividade em pressão ambiente em filmes finos, determinar a origem da redução de $T_c$ e investigar a robustez da simetria de emparelhamento contra variações estruturais e eletrônicas.

2. Metodologia

O estudo emprega uma abordagem de múltiplos passos combinando cálculos de primeiros princípios com teoria de perturbação de muitos corpos:

• Otimização Estrutural:

  • Cálculos de primeiros princípios foram realizados usando VASP com o funcional PBEsol.
  • As constantes de rede no plano ($a$) foram fixadas para corresponder a três substratos experimentais: LSAT, LAO e SLAO.
  • A constante de rede fora do plano ($c$) e as posições atômicas foram otimizadas para minimizar a energia total.
  • Duas suposições de simetria foram testadas: a $I4/mmm$ (tetragonal) sugerida experimentalmente e a $Amam$ (ortorrômbica) otimizada teoricamente.
  • Um modelo adicional foi construído usando a estrutura determinada experimentalmente da Ref. [24] (Yue et al.), que apresenta uma distância Ni-Ni distinta.

• Construção do Modelo:

  • Funções de Wannier: Funções de Wannier maximamente localizadas foram extraídas usando WANNIER90, focando nos orbitais Ni-$3d_{x^2-y^2}$ e $3d_{3z^2-r^2}$.
  • Hamiltoniano: Um modelo de ligação forte de bicamada de dois orbitais foi construído.
  • Parâmetros de Interação: Repulsão Coulombiana no sítio ($U$), repulsão interorbital ($U'$), acoplamento de Hund ($J$) e salto de pares ($J'$) foram adotados de cálculos anteriores de RPA restrito ($U \approx 3,6-3,8$ eV).
  • Estrutura de Banda: Cálculos foram realizados com e sem a correção +U (formulação de Dudarev, $U=3$ eV) para avaliar a sensibilidade.

• Análise de Supercondutividade (FLEX):

  • A aproximação de Troca de Flutuações (FLEX) foi aplicada ao modelo de Hubbard de dois orbitais.
  • Este método calcula autoconsistentemente a autoenergia e a interação de emparelhamento mediada por flutuações de spin, levando em conta as dependências completas de momento e frequência da função de Green.
  • A equação de Eliashberg linearizada foi resolvida para obter o autovalor $\lambda$ (um proxy para $T_c$) e a função de gap supercondutor $\Delta(k)$.

3. Contribuições e Resultados Principais

A. Tendências Estruturais e Eletrônicas

• Ângulo de Ligação: À medida que a constante de rede no plano $a$ diminui (aumentando a tensão compressiva), o ângulo de ligação intercamada Ni-O-Ni ($\alpha$) aproxima-se de $180^\circ$, consistente com a transição para a simetria $I4/mmm$. No entanto, mesmo para o substrato SLAO, a otimização DFT não atinge exatamente $180^\circ$ (a simetria $Amam$ persiste ligeiramente), embora a diferença de energia entre $Amam$e$I4/mmm$ seja negligenciável.
• Ocupação Orbital: Reduzir $a$ aumenta a altura do oxigênio apical e o deslocamento de nível ($\Delta E$) entre os orbitais $d_{x^2-y^2}$ e $d_{3z^2-r^2}$. Isso desloca elétrons para o orbital $d_{3z^2-r^2}$, aproximando-o do preenchimento metade (ocupação $\approx 0,95$ para SLAO).
• Integrais de Salto:
  • O salto intercamada $|t_\perp|$ aumenta à medida que $a$ diminui em estruturas otimizadas por DFT (atingindo $\sim 0,6$ eV).
  • Crucialmente, a estrutura determinada experimentalmente produz um $|t_\perp|$ muito menor $\approx 0,45$ eV devido a uma distância Ni-Ni maior.

B. Robustez do Emparelhamento de Onda $s_{\pm}$

• Estabilidade de Simetria: O estudo encontra que a simetria de emparelhamento de onda $s_{\pm}$ é robusta. Ela permanece a solução dominante independentemente de:
  • A simetria cristalina específica ($Amam$ vs. $I4/mmm$).
  • A presença ou ausência do $\gamma$-pocket (que depende da correção $U$ e do preenchimento de banda).
  • Os detalhes específicos da topologia da superfície de Fermi.
• Características da Função de Gap:
  • Na representação orbital, a função de gap é dominada pelo emparelhamento intercamada dos orbitais $d_{3z^2-r^2}$.
  • Este gap intercamada é quase independente do momento (emparelhamento intra-célula unitária).
  • A função de gap na representação de banda reflete o peso orbital: é grande onde há caráter $d_{3z^2-r^2}$ e muda de sinal entre as bandas de ligação e antiligação.

C. Origem da $T_c$ Reduzida em Filmes Finos

• Papel de $t_\perp$: A razão principal para a $T_c$ reduzida em filmes finos (comparada ao bulk sob pressão) é identificada como a redução do integral de salto intercamada $|t_\perp|$.
  • Modelos baseados em estruturas otimizadas por DFT (grande $|t_\perp| \approx 0,6$ eV) preveem um alto $\lambda$ semelhante ao bulk sob pressão.
  • Modelos baseados na estrutura determinada experimentalmente (pequeno $|t_\perp| \approx 0,45$ eV) produzem um $\lambda$ significativamente suprimido, consistente com a $T_c$ observada $\approx 40$ K.
• Mecanismo: A supressão de $T_c$ com a redução de $|t_\perp|$ é atribuída à natureza da cola de emparelhamento. O emparelhamento é mediado por flutuações de spin de energia finita (em vez de aninhamento de superfície de Fermi de baixa energia).
  • Flutuações de spin de alta energia (que impulsionam o emparelhamento) são suprimidas quando $|t_\perp|$ é reduzido.
  • Por outro lado, flutuações de baixa energia (associadas ao aninhamento da superfície de Fermi) são aprimoradas pela redução de $|t_\perp|$, mas estas são menos eficazes para o emparelhamento $s_{\pm}$ neste sistema.

4. Significado e Conclusão

• Resolução do Enigma da $T_c$: O artigo fornece uma explicação teórica coerente para a $T_c$ reduzida em filmes finos em pressão ambiente: a estrutura cristalina experimental possui inerentemente um salto intercamada menor ($t_\perp$) do que o previsto pela otimização DFT padrão, o que suprime diretamente a força de emparelhamento mediada por flutuações de spin de energia finita.
• Validação do Mecanismo de Emparelhamento: O estudo reforça a visão que a supercondutividade em $La_3Ni_2O_7$ é impulsionada por flutuações de spin de energia finita mediando um emparelhamento intercamada $d_{3z^2-r^2}$. Este mecanismo é robusto contra variações na topologia da superfície de Fermi (como a presença do $\gamma$-pocket) e detalhes de simetria estrutural.
• Consistência Teórica: Os resultados alinham-se com abordagens não perturbativas (cDMFT, DCA, VMC), confirmando que o sistema pode ser visto como um modelo de Hubbard de bicamada (formado por orbitais $d_{3z^2-r^2}$) embutido em um ambiente multiorbital.
• Direções Futuras: A grande discrepância entre as estruturas cristalinas previstas pela DFT e as observadas experimentalmente (especificamente a distância Ni-Ni) permanece uma questão em aberto. Os autores sugerem que ordens concorrentes (e.g., ondas de densidade) não capturadas no esquema FLEX atual também podem desempenhar um papel na estabilidade e propriedades do material real.

Em resumo, o artigo estabelece que, embora a simetria de emparelhamento ($s_{\pm}$) seja robusta, a magnitude de $T_c$ é altamente sensível ao salto intercamada $t_\perp$, que é ditado pela estrutura atômica precisa do filme fino.