Ah-SCDFT:A general approach for superconductivity with an-harmonic corrections

Resumen Técnico: Ah-SCDFT – Un enfoque general para la superconductividad con correcciones anharmónicas

Planteamiento del Problema
Los estudios de primeros principios sobre superconductividad dependen frecuentemente de la aproximación armónica, omitiendo los efectos anharmónicos (AHE). Si bien la aproximación armónica es computacionalmente eficiente, no logra capturar fenómenos físicos cruciales en sistemas con grados de libertad blandos, redes con inestabilidad dinámica (indicada por frecuencias de fonones imaginarias) o grandes desplazamientos atómicos a altas temperaturas. La omisión de AHE puede conducir a errores cuantitativos significativos y clasificaciones cualitativas erróneas, tales como predecir incorrectamente la inestabilidad termodinámica, identificar erróneamente ondas de densidad de carga o subestimar sustancialmente la fuerza del acoplamiento electrón-fonón (EPC) y la temperatura crítica de superconductividad ($T_c$). Los métodos existentes, como las fórmulas de McMillan/Allen-Dynes, dependen de parámetros empíricos, mientras que las ecuaciones de Migdal-Eliashberg anisotrópicas estándar y la Teoría del Funcional de la Densidad de la Superconductividad (SCDFT) operan típicamente dentro del marco armónico.

Metodología: Ah-SCDFT
Los autores introducen la "teoría del funcional de la densidad de la superconductividad anharmónica" (ah-SCDFT), un enfoque computacional unificado que integra sistemáticamente las correcciones anharmónicas en la SCDFT estándar. La metodología procede mediante un bucle iterativo autoconsistente:

1. Cálculo Inicial: Se obtiene una geometría de estado fundamental ($S_1$) a 0 K mediante la DFT estándar dentro de la aproximación armónica.
2. Evaluación de la Superconductividad: Las propiedades superconductoras, incluyendo una temperatura crítica inicial ($T_{c1}$), se calculan utilizando formalismos establecidos (ecuaciones de Migdal-Eliashberg anisotrópicas o SCDFT).
3. Renormalización de la Red: Se determina una estructura de equilibrio dependiente de la temperatura para la $T_{c1}$ predicha. Esto se logra minimizando la energía libre de Helmholtz ($F = F_0 + F_{vib}$), donde $F_{vib}$ es la energía libre vibracional calculada utilizando estructuras de equilibrio a temperatura finita obtenidas mediante la aproximación cuasi-armónica (QHA) y el código thermo_pw.
4. Iteración: Los cálculos de superconductividad se repiten utilizando esta geometría renormalizada para obtener una temperatura crítica actualizada ($T_{c2}$).
5. Convergencia: Este proceso itera hasta que la diferencia entre las temperaturas críticas sucesivas satisface $|T_{cn} - T_{cn-1}| < \epsilon$ (por ejemplo, 0.01 K), produciendo la $T_c$ final autoconsistente.
   Los cálculos utilizan el paquete QUANTUM ESPRESSO (QE) y el SUPERCONDUCTING-TOOLKIT, empleando la aproximación de gradiente generalizado PBE y el kernel de intercambio-correlación Sanna2020.

Resultos Clave: Benchmark de MgB2
Los autores validan el enfoque ah-SCDFT utilizando el prototípico superconductor s-wave MgB2 tanto en condiciones ambientales como bajo presión externa.

• Condiciones Ambientales:

  • Aproximación Armónica: Predice una $T_c$ de 36 K, consistente con informes previos pero inferior al valor experimental.
  • Ah-SCDFT: Al incorporar la renormalización de la red anharmónica a 39 K, la $T_c$ predicha aumenta a 39 K, coincidiendo con los datos experimentales.
  • Mecanismo: El aumento se atribuye a que los efectos anharmónicos incrementan la densidad de estados (DOS) electrónica en el nivel de Fermi. Esto aumenta la constante de EPC ($\lambda$) de 0.82 (armónico) a 0.88 (anharmónico). Específicamente, las correcciones anharmónicas suavizan las ramas de fonones asociadas con los átomos de Mg y B, particularmente cerca del punto A, donde las vibraciones en el plano de B exhiben el acoplamiento más fuerte.

• Bajo Presión Externa:

  • El estudio aborda la literatura contradictoria sobre la dependencia de la $T_c$ de MgB2 con la presión (algunos informes sugieren una disminución monotónica, otros una tendencia no monotónica con un punto de inflexión).
  • Hallazgos de Ah-SCDFT: El enfoque reproduce una tendencia no monotónica donde la $T_c$ disminuye inicialmente y luego se recupera, pero desplaza el punto de inflexión hacia presiones significativamente más altas (alrededor de 200 GPa) en comparación con las predicciones anteriores basadas en el modelo armónico (~100 GPa).
  • Efectos de la Presión: A medida que la presión aumenta, la red se endurece, causando que la diferencia entre las predicciones de $T_c$ armónica y anharmónica disminuya. El método también predice inestabilidad estructural por encima de 340 GPa, señalada por frecuencias de fonones imaginarias.

Significancia y Reivindicaciones
El artículo afirma que la ah-SCDFT establece un marco poderoso, eficiente y ampliamente aplicable para el estudio cuantitativo de la superconductividad. Su principal significancia radica en:

1. Precisión: Logra predicciones de alta fidelidad de las propiedades superconductoras (específicamente la $T_c$) al contabilizar explícitamente la renormalización de la red dependiente de la temperatura, resolviendo las discrepancias entre la teoría y el experimento observadas en materiales como el MgB2.
2. Eficiencia: El método incurre solo en un aumento modesto del costo computacional, requiriendo un número limitado de pasos de convergencia para el bucle iterativo, lo que lo hace factible para una aplicación más amplia.
3. Fiabilidad: Proporciona un enfoque consistente para comprender la superconductividad dependiente de la presión y ofrece una herramienta robusta para la predicción y el diseño de nuevos materiales superconductores donde la anharmonicidad juega un papel crítico.