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 estudio teórico demuestra que, si bien la estructura electrónica de las películas delgadas de La$_3$Ni$_2$O$_7$ a presión ambiente varía con la estructura cristalina y los detalles computacionales, la robustez del apareamiento de onda $s\pm$ mediado por fluctuaciones de espín de energía finita persiste, aunque la reducción observada en la temperatura crítica en comparación con el volumen bajo presión se explica mejor mediante modelos con pequeño salto intercapa derivados de estructuras de red experimentales.

Lo esencial

La Gran Imagen: Un Superconductor a Presión Ambiente

Imagina un material llamado La₃Ni₂O₇ (un tipo de cristal a base de níquel) que puede conducir electricidad sin resistencia alguna (superconductividad). Los científicos descubrieron recientemente que si se comprime este material con una presión masiva, se convierte en un superconductor a unos 80 Kelvin (muy frío, pero cálido para los superconductores).

Recientemente, los investigadores encontraron una manera de lograr esto sin comprimirlo. Hicieron crecer el material como una película muy delgada sobre un tipo específico de "suelo" cristalino (un sustrato). El suelo era ligeramente más pequeño que la película, lo que la comprimía desde los lados, imitando el efecto de la alta presión. Esta película se volvió superconductora a alrededor de 40 Kelvin.

La Pregunta: ¿Por qué funciona la película delgada a una temperatura más baja (40 K) que el material masivo comprimido (80 K)? Y, ¿cuál es la "receta" exacta dentro del material que hace que la electricidad fluya sin resistencia?

El Enfoque de los Científicos: Construyendo un Modelo Digital

Los autores de este artículo no solo adivinaron; construyeron una simulación por computadora detallada. Piensa en ello como un motor de videojuegos donde intentaron recrear la física de este material desde cero.

1. El Plano (Estructura): Observaron el "plano" de los átomos. Probaron dos planos diferentes:
   • El Plano Teórico: Lo que sus cálculos por computadora decían que los átomos deberían parecer.
   • El Plano Experimental: Lo que los científicos midieron realmente en el laboratorio recientemente.
2. El Motor (FLEX): Utilizaron un motor matemático complejo llamado FLEX (Aproximación de Intercambio de Fluctuaciones). Imagina este motor como un simulador de clima superpreciso. En lugar de predecir la lluvia, predice cómo bailan e interactúan entre sí los electrones (las diminutas partículas que transportan la electricidad). Tiene en cuenta cada movimiento posible que pueden hacer los electrones, no solo los obvios.

Hallazgos Clave: El "Baile" de los Electrones

1. El Misterio de la "Bolsa-γ"

En el mundo de estos materiales, hay una forma específica de la multitud de electrones llamada la bolsa-γ. Algunos científicos pensaban que esta bolsa era esencial para la superconductividad; otros pensaban que no importaba.

• El Veredicto del Artículo: Los autores descubrieron que si esta "bolsa-γ" existe o desaparece depende enteramente de qué plano uses (teórico vs. experimental) y de cómo ajustes las matemáticas.
• La Analogía: Es como mirar a una multitud a través de gafas de diferentes colores. A través de una lente, ves un grupo distinto de personas (la bolsa); a través de otra, se mezclan.
• El Resultado: Sorprendentemente, no importaba. Ya fuera que la bolsa estuviera presente o no, la superconductividad permaneció fuerte. El "baile" de los electrones era lo suficientemente robusto para manejar estos cambios estructurales.

2. El "Pegamento" que lo Sostiene

¿Cómo se emparejan los electrones para conducir la electricidad? Por lo general, necesitan un "pegamento".

• El Veredicto del Artículo: El pegamento aquí son las fluctuaciones de espín. Imagina que los electrones son bailarines girando. A veces, tambalean o fluctúan en su espín. Estos tambaleos actúan como un ritmo que ayuda a los bailarines a emparejarse.
• El Giro: El artículo argumenta que este "ritmo" proviene de tambaleos de alta energía, no solo de los movimientos lentos y obvios cerca de la superficie de la multitud de electrones. Dado que el pegamento se basa en estos tambaleos de alta energía, la superconductividad es muy estable y no se rompe fácilmente si la forma de la multitud de electrones cambia ligeramente.

3. ¿Por qué la Película es más Fría (40 K) que el Material Masivo (80 K)?

Este fue el mayor acertijo. La película delgada es superconductora a la mitad de la temperatura del material masivo presurizado.

• El Veredicto del Artículo: La diferencia se reduce a un número específico: $|t_\perp|$.
• La Analogía: Piensa en el material como un edificio de dos pisos donde los electrones pueden saltar entre los pisos.
  • En el material masivo presurizado, los pisos están perfectamente alineados y el salto entre ellos es fuerte y fácil (Alto $|t_\perp|$). Esto crea una pista de baile muy eficiente (80 K).
  • En la película delgada, las mediciones experimentales muestran que los pisos están ligeramente desalineados o que el salto es más débil (Bajo $|t_\perp|$).
• La Conclusión: Cuando los autores utilizaron el "Plano Experimental" (que mostraba este salto más débil) en su simulación, la temperatura superconductora cayó exactamente a los 40 K observados. Cuando usaron el "Plano Teórico" (que predecía un salto más fuerte), la temperatura se mantuvo alta en 80 K.
• La Lección: La razón por la que la película es "más débil" es probablemente porque la estructura física real de la película tiene una conexión más débil entre sus capas de lo que predijo la teoría.

Resumen en Poca Cosa

Los científicos construyeron una simulación de alta tecnología para entender por qué una nueva película superconductora funciona a presión ambiente. Descubrieron que:

1. El mecanismo de emparejamiento es resistente: Los electrones se emparejan usando "tambaleos" de alta energía (fluctuaciones de espín), lo que hace que la superconductividad sea muy robusta frente a pequeños cambios en la forma del material.
2. La "bolsa-γ" no importa: Si existe o no una forma específica de electrones no cambia el resultado.
3. La caída de temperatura es estructural: La película solo alcanza los 40 K (en lugar de 80 K) porque la distancia física real entre las capas atómicas en la película es ligeramente diferente de lo que predijo la teoría, haciendo que el "salto" entre capas sea más débil.

El artículo esencialmente dice: "Conocemos la receta para la superconductividad en este material. La razón por la que la película es ligeramente menos eficiente que el bloque presurizado es simplemente porque las capas de la película no están tan perfectamente conectadas como pensábamos que estarían".

Resumen técnico

Aquí se presenta un resumen técnico detallado del artículo "Estudio teórico sobre la superconductividad a presión ambiente en películas delgadas de La3Ni2O7: análisis estructural, construcción de modelos y robustez del apareamiento de onda s±".

1. Planteamiento del Problema

El descubrimiento de la superconductividad a alta temperatura ($T_c \approx 80$ K) en el nickelato bicapa presurizado $La_3Ni_2O_7$ despertó un intenso interés. Recientemente, se observó superconductividad a presión ambiente ($T_c \approx 40$ K) en películas delgadas de $La_3Ni_2O_7$ crecidas sobre sustratos como $SrLaAlO_4$ (SLAO). Sin embargo, persiste un desafío teórico significativo:

• La Discrepancia de $T_c$: La $T_c$ en las películas delgadas es aproximadamente la mitad de la del sistema masivo presurizado.
• Ambigüedad Estructural: Existe una gran desviación entre las estructuras cristalinas determinadas por la optimización de la Teoría del Funcional de la Densidad (DFT) y las determinadas experimentalmente. Específicamente, la estructura experimental muestra un integral de salto intercapa significativamente menor ($|t_\perp| \approx 0.4$ eV) en comparación con las predicciones de la DFT ($|t_\perp| \approx 0.6$ eV).
• Controversia de la Superficie de Fermi: Los estudios experimentales y teóricos discrepan sobre la existencia de la "$\gamma$-bolsa" (una superficie de Fermi que se origina en la parte superior de la banda de enlace $d_{3z^2-r^2}$). Su presencia o ausencia se debate en relación con su impacto en la simetría de apareamiento.
• Objetivo: Los autores buscan reproducir teóricamente la superconductividad a presión ambiente en películas delgadas, determinar el origen de la reducción de $T_c$ e investigar la robustez de la simetría de apareamiento frente a variaciones estructurales y electrónicas.

2. Metodología

El estudio emplea un enfoque de múltiples pasos que combina cálculos de primeros principios con la teoría de perturbaciones de muchos cuerpos:

• Optimización Estructural:

  • Se realizaron cálculos de primeros principios utilizando VASP con el funcional PBEsol.
  • Las constantes de red en el plano ($a$) se fijaron para coincidir con tres sustratos experimentales: LSAT, LAO y SLAO.
  • La constante de red fuera del plano ($c$) y las posiciones atómicas se optimizaron para minimizar la energía total.
  • Se probaron dos suposiciones de simetría: la $I4/mmm$ (tetragonal) sugerida experimentalmente y la $Amam$ (ortorrómbica) optimizada teóricamente.
  • Se construyó un modelo adicional utilizando la estructura determinada experimentalmente de la Ref. [24] (Yue et al.), que presenta una distancia Ni-Ni distinta.

• Construcción del Modelo:

  • Funciones de Wannier: Se extrajeron funciones de Wannier máximamente localizadas utilizando WANNIER90, centrándose en los orbitales Ni-$3d_{x^2-y^2}$ y $3d_{3z^2-r^2}$.
  • Hamiltoniano: Se construyó un modelo de enlace fuerte bicapa de dos orbitales.
  • Parámetros de Interacción: La repulsión de Coulomb en el sitio ($U$), la repulsión interorbital ($U'$), el acoplamiento de Hund ($J$) y el salto de pares ($J'$) se adoptaron de cálculos previos de RPA restringida ($U \approx 3.6-3.8$ eV).
  • Estructura de Bandas: Se realizaron cálculos con y sin la corrección +U (formulación de Dudarev, $U=3$ eV) para evaluar la sensibilidad.

• Análisis de Superconductividad (FLEX):

  • Se aplicó la aproximación de Intercambio de Fluctuaciones (FLEX) al modelo de Hubbard de dos orbitales.
  • Este método calcula autoconsistentemente la autoenergía y la interacción de apareamiento mediada por fluctuaciones de espín, teniendo en cuenta las dependencias completas del momento y la frecuencia de la función de Green.
  • Se resolvió la ecuación de Eliashberg linealizada para obtener el autovalor $\lambda$ (un proxy para $T_c$) y la función de brecha superconductora $\Delta(k)$.

3. Contribuciones y Resultados Clave

A. Tendencias Estructurales y Electrónicas

• Ángulo de Enlace: A medida que disminuye la constante de red en el plano $a$ (aumentando la tensión compresiva), el ángulo de enlace intercapa Ni-O-Ni ($\alpha$) se acerca a $180^\circ$, consistente con la transición hacia la simetría $I4/mmm$. Sin embargo, incluso para el sustrato SLAO, la optimización DFT no alcanza exactamente $180^\circ$ (la simetría $Amam$ persiste ligeramente), aunque la diferencia de energía entre $Amam$e$I4/mmm$ es negligible.
• Ocupación Orbital: Reducir $a$ aumenta la altura del oxígeno apical y el desplazamiento de nivel ($\Delta E$) entre los orbitales $d_{x^2-y^2}$ y $d_{3z^2-r^2}$. Esto desplaza electrones hacia el orbital $d_{3z^2-r^2}$, acercándolo a la mitad de llenado (ocupación $\approx 0.95$ para SLAO).
• Integrales de Salto:
  • El salto intercapa $|t_\perp|$ aumenta a medida que disminuye $a$ en las estructuras optimizadas por DFT (alcanzando $\sim 0.6$ eV).
  • Crucialmente, la estructura determinada experimentalmente produce un $|t_\perp|$ mucho menor, $\approx 0.45$ eV, debido a una mayor distancia Ni-Ni.

B. Robustez del Apareamiento de Onda $s_{\pm}$

• Estabilidad de Simetría: El estudio encuentra que la simetría de apareamiento de onda $s_{\pm}$ es robusta. Sigue siendo la solución dominante independientemente de:
  • La simetría cristalina específica ($Amam$ vs. $I4/mmm$).
  • La presencia o ausencia de la $\gamma$-bolsa (que depende de la corrección $U$ y del llenado de bandas).
  • Los detalles específicos de la topología de la superficie de Fermi.
• Características de la Función de Brecha:
  • En la representación orbital, la función de brecha está dominada por el apareamiento intercapa de orbitales $d_{3z^2-r^2}$.
  • Esta brecha intercapa es casi independiente del momento (apareamiento intra-celda unitaria).
  • La función de brecha en la representación de bandas refleja el peso orbital: es grande donde existe carácter $d_{3z^2-r^2}$ y cambia de signo entre las bandas de enlace y antienlace.

C. Origen de la $T_c$ Reducida en Películas Delgadas

• Rol de $t_\perp$: La razón principal de la $T_c$ reducida a la mitad en las películas delgadas (en comparación con el masivo presurizado) se identifica como la reducción del integral de salto intercapa $|t_\perp|$.
  • Los modelos basados en estructuras optimizadas por DFT (gran $|t_\perp| \approx 0.6$ eV) predicen un alto $\lambda$ similar al del masivo presurizado.
  • Los modelos basados en la estructura determinada experimentalmente (pequeño $|t_\perp| \approx 0.45$ eV) producen un $\lambda$ significativamente suprimido, consistente con la $T_c$ observada de $\approx 40$ K.
• Mecanismo: La supresión de $T_c$ con la reducción de $|t_\perp|$ se atribuye a la naturaleza del "pegamento" de apareamiento. El apareamiento está mediado por fluctuaciones de espín de energía finita (en lugar del anidamiento de superficie de Fermi de baja energía).
  • Las fluctuaciones de espín de alta energía (que impulsan el apareamiento) se suprimen cuando $|t_\perp|$ se reduce.
  • Por el contrario, las fluctuaciones de baja energía (asociadas al anidamiento de la superficie de Fermi) se ven potenciadas por la reducción de $|t_\perp|$, pero estas son menos efectivas para el apareamiento $s_{\pm}$ en este sistema.

4. Significado y Conclusión

• Resolución del Enigma de $T_c$: El artículo proporciona una explicación teórica coherente para la $T_c$ reducida en películas delgadas a presión ambiente: la estructura cristalina experimental posee intrínsecamente un salto intercapa ($t_\perp$) menor que el predicho por la optimización DFT estándar, lo que suprime directamente la fuerza de apareamiento mediada por fluctuaciones de espín de energía finita.
• Validación del Mecanismo de Apareamiento: El estudio refuerza la visión de que la superconductividad en $La_3Ni_2O_7$ está impulsada por fluctuaciones de espín de energía finita que median un apareamiento intercapa de $d_{3z^2-r^2}$. Este mecanismo es robusto frente a variaciones en la topología de la superficie de Fermi (como la presencia de la $\gamma$-bolsa) y detalles de simetría estructural.
• Consistencia Teórica: Los resultados se alinean con enfoques no perturbativos (cDMFT, DCA, VMC), confirmando que el sistema puede verse como un modelo de Hubbard bicapa (formado por orbitales $d_{3z^2-r^2}$) incrustado en un entorno multiorbital.
• Futuras Direcciones: La gran discrepancia entre las estructuras cristalinas predichas por DFT y las observadas experimentalmente (específicamente la distancia Ni-Ni) sigue siendo una pregunta abierta. Los autores sugieren que órdenes competidores (por ejemplo, ondas de densidad) no capturados en el esquema FLEX actual también podrían desempeñar un papel en la estabilidad y propiedades reales del material.

En resumen, el artículo establece que, si bien la simetría de apareamiento ($s_{\pm}$) es robusta, la magnitud de $T_c$ es altamente sensible al salto intercapa $t_\perp$, el cual está dictado por la estructura atómica precisa de la película delgada.