⚛️ quantum physics

Computing band gaps of periodic materials via sample-based quantum diagonalization

Este artículo introduce un flujo de trabajo de diagonalización cuántica basado en muestreo que predice con precisión las brechas de banda de materiales periódicos como el dióxido de hafnio y el dióxido de circonio mediante la combinación de Hamiltonianos de red autoconsistentes con el muestreo de procesadores cuánticos superconductores, superando a la teoría del funcional de la densidad tradicional y alineándose con los resultados experimentales.

Autores originales: Alan Duriez, Pamela C. Carvalho, Marco Antonio Barroca, Federico Zipoli, Ben Jaderberg, Rodrigo Neumann Barros Ferreira, Kunal Sharma, Antonio Mezzacapo, Benjamin Wunsch, Mathias Steiner

Publicado 2026-07-14

📖 4 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Alan Duriez, Pamela C. Carvalho, Marco Antonio Barroca, Federico Zipoli, Ben Jaderberg, Rodrigo Neumann Barros Ferreira, Kunal Sharma, Antonio Mezzacapo, Benjamin Wunsch, Mathias Steiner

Artículo original bajo licencia CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Esta es una explicación generada por IA del artículo a continuación. No ha sido escrita ni avalada por los autores. Para mayor precisión técnica, consulte el artículo original. Leer descargo de responsabilidad completo

Imagina intentar predecir cómo un material conduce la electricidad o la bloquea, como si intentaras adivinar la altura de un muro entre dos habitaciones. Durante décadas, los científicos han utilizado una receta estándar llamada Teoría del Funcional de la Densidad (DFT) para hacer estas conjeturas. Es como usar un mapa que asume que todo el mundo camina en línea recta e ignora el hecho de que las personas chocan entre sí. Esto funciona bien para algunos materiales, pero cuando llegas a los complicados, como el dióxido de hafnio (HfO2) y el dióxido de zirconio (ZrO2) —materiales utilizados en todo, desde chips de computadora hasta aislantes— el mapa de la "línea recta" falla. Este ignora las interacciones caóticas y accidentadas entre los electrones, lo que conduce a predicciones que son demasiado bajas en comparación con lo que sucede en la vida real.

Entra en escena un nuevo equipo híbrido: una mezcla de supercomputadoras clásicas y un nuevo tipo de computadora cuántica. Los investigadores no solo lanzaron una computadora cuántica al problema; construyeron un puente especial. Primero, utilizaron la DFT tradicional para obtener un boceto aproximado de la estructura electrónica del material. Luego, tradujeron ese boceto en un "Hamiltoniano de red", que es como convertir el material en una cuadrícula de diminutas islas que interactúan entre sí. En esta cuadrícula, añadieron reglas específicas sobre cómo se repelen los electrones (la parte de los "choques") utilizando parámetros que ellos mismos calcularon.

Para resolver el rompecabezas de cómo se comportan estos electrones en esta cuadrícula, utilizaron un método llamado Diagonalización Cuántica Basada en Muestreo (SQD). Piensa en la computadora cuántica no como una calculadora que resuelve toda la ecuación a la vez, sino como un muestreador superrápido. Ejecuta un circuito cuántico (un conjunto específico de instrucciones) millones de veces para recolectar una lista masiva de "instantáneas" o configuraciones de los electrones. En este estudio, ejecutaron estos circuitos en un procesador cuántico superconductor de última generación llamado ibm pittsburgh, que tiene 156 cúbits. Para el HfO2, utilizaron un circuito con 38 cúbits, y para el ZrO2, uno de 46 cúbits.

La magia ocurre cuando toman estas instantáneas cuánticas y las alimentan a una computadora clásica. La computadora clásica observa esta lista curada de las disposiciones electrónicas más importantes y realiza las matemáticas pesadas para encontrar el estado de menor energía. Es como si la computadora cuántica fuera un explorador encontrando los mejores caminos a través de un bosque, y la computadora clásica fuera el cartógrafo dibujando el mapa final basado en esos caminos.

¿Los resultados? El equipo descubrió que esta combinación cuántica-clásica predijo la "brecha de banda" (la pared de energía entre los estados electrónicos ocupados y vacíos) para ambos, HfO2 y ZrO2, con una precisión impresionante. Para el ZrO2, los métodos estándar de DFT fallaron por completo, perdiendo el rango experimental por completo. El método cuántico, sin embargo, aterrizó justo en el punto ideal de lo que los experimentos de laboratorio independientes han medido. De hecho, la diferencia entre su predicción cuántica y un método de referencia clásico altamente preciso (llamado HCI) fue minúscula: solo 1.6 × 10⁻⁸ eV para el HfO2 y 1.3 × 10⁻⁷ eV para el ZrO2.

Es importante destacar que el artículo argumenta explícitamente en contra de la idea de que el enfoque estándar de DFT es suficiente para estos materiales, señalando que subestima la brecha de banda porque ignora las fuertes correlaciones electrónicas. También demuestran que, si bien añadir solo un tipo de corrección (interacciones intra-sitio, o "U") ayuda un poco, es la combinación de ambas, las interacciones intra-sitio e inter-sitio (la parte "V", que representa cómo los electrones en diferentes átomos se comunican entre sí), lo que es crucial. Su método superó no solo a la DFT, sino también a otros bancos de pruebas de química cuántica avanzada como el CCSD(T).

Aunque los resultados son un paso significativo hacia adelante, los autores tienen cuidado de enmarcar esto como una simulación exitosa en hardware actual, "pre-tolerante a fallos". No afirmaron haber resuelto todos los problemas de la ciencia de materiales; más bien, demostraron que este flujo de trabajo específico puede producir predicciones útiles y verificables experimentalmente para estos dos dieléctricos específicos. El método es rápido y eficiente en recursos en comparación con otros métodos de alta precisión como la aproximación GW, lo que sugiere que, con mejores procesadores cuánticos en el futuro, este enfoque podría escalarse para abordar materiales aún más complejos. Por ahora, demuestra que incluso con las máquinas cuánticas ruidosas de hoy, podemos empezar a ver la verdadera forma del mundo electrónico.

Resumen Técnico: Computación de Brechas de Banda de Materiales Periódicos mediante la Diagonalización Cuántica Basada en Muestreo

Planteamiento del Problema
La predicción de las propiedades electrónicas de materiales periódicos, particularmente las brechas de banda (band gaps), sigue siendo un desafío significativo en la física del estado sólido computacional. La Teoría del Funcional de la Densidad (DFT) estándar, aunque ampliamente utilizada, depende de aproximaciones de campo medio y de la suposición de electrones no interactuantes. Estas suposiciones a menudo fallan al capturar las correlaciones electrónicas fuertes presentes en los óxidos de metales de transición, lo que conduce a una subestimación sistemática de las brechas de banda en semiconductores y aislantes de banda ancha. Si bien métodos como el modelo de Hubbard Extendido (EH) pueden incorporar efectos de correlación mediante parámetros intra-sitio ( $U$ ) e inter-sitio ( $V$ ), resolver los Hamiltonianos de red fermiónica resultantes para sistemas grandes es computacionalmente prohibitivo para algoritmos clásicos como el Monte Carlo Cuántico (QMC) o el Grupo de Renormalización de Matriz de Densidad (DMRG), que enfrentan limitaciones en su aplicabilidad a sistemas fermiónicos o en su escalabilidad a tres dimensiones.

Metodología
Los autores proponen un flujo de trabajo híbrido cuántico-clásico para simular estados electrónicos y predecir brechas de banda para materiales periódicos. La metodología consta de tres etapas primarias:

Construcción del Hamiltoniano de Red:
- La estructura electrónica se calcula primero mediante DFT (específicamente DFT+U+V) para obtener las densidades del estado fundamental y los parámetros de Hubbard auto-consistentes ( $U$ y $V$ ).
- Estos resultados de DFT se proyectan sobre un conjunto de bases de orbitales atómicos localizados para construir un Hamiltoniano de enlace fuerte (tight-binding).
- Esta proyección genera un Hamiltoniano de red general (Ecuación 1) definido por parámetros de salto ( $t_{pq}$ ) y parámetros de interacción intra-sitio ( $U_p$ ) e inter-sitio ( $V_{pq}$ ) auto-consistentes. El sistema se mapea de tal manera que cada orbital atómico corresponde a un único sitio en la red.
- El Hamiltoniano se expresa en el espacio recíproco pero se investiga específicamente en el punto $\Gamma$ .
Muestreo Cuántico y Diagonalización:
- El Hamiltoniano de red se mapea sobre un Hamiltoniano electrónico generalizado adecuado para el muestreo cuántico, tratando efectivamente el material periódico como un sistema molecular para fines del algoritmo.
- Se emplea el método de Diagonalización Cuántica Basada en Muestreo (SQD). Esto implica:
  - Utilizar un ansatz de Jastrow de Clúster Unitario Local (LUCJ) para aproximar la función de onda del estado fundamental. Los parámetros del ansatz se inicializan utilizando las amplitudes de un cálculo clásico de Coupled-Cluster Singles and Doubles (CCSD).
  - Ejecutar el circuito cuántico parametrizado en un procesador cuántico superconductor (IBM ibm pittsburgh con un QPU Heron r3) para muestrear configuraciones electrónicas (bitstrings).
  - Construir un subespacio de configuración reducido a partir de estas muestras y diagonalizar clásicamente el Hamiltoniano proyectado dentro de este subespacio para obtener las energías del estado fundamental.
- Se utiliza una versión extendida, Ext-SQD, para mejorar la precisión. Este algoritmo elimina las configuraciones de baja probabilidad del conjunto inicial de SQD y aplica operadores de excitación para generar un conjunto de configuración extendido, permitiendo cálculos de energía refinados con espacios de configuración más pequeños.
Cálculo de la Brecha de Banda:
- La brecha de banda directa en el punto $\Gamma$ se deriva de las energías del estado fundamental del Hamiltoniano en tres sectores de simetría: $N_e$ (electrones de valencia), $N_e - 1$ , y $N_e + 1$ .
- La brecha de banda ( $E_g$ ) se calcula mediante la fórmula: $E_g = E[N_e - 1] + E[N_e + 1] - 2E[N_e]$ .

Contribuciones Clave

Representación de Red Generalizada: El artículo introduce un marco general para representar la estructura electrónica de materiales periódicos como un Hamiltoniano de red con parámetros de interacción específicos del material y derivados de forma auto-consistente, cerrando la brecha entre la DFT y la simulación cuántica.
Aplicación de SQD a Sólidos: Los autores adaptan con éxito el método SQD, aplicado previamente a estructuras moleculares, a materiales periódicos mediante el mapeo del modelo de Hubbard Extendido a una forma compatible con las técnicas de muestreo de química cuántica.
Demostración de Hardware: El flujo de trabajo se ejecuta en un procesador cuántico superconductor de última generación (pre-tolerante a fallos, 156 qubits) para simular dióxido de hafnio (HfO $_2$ ) y dióxido de circonio (ZrO $_2$ ).

Resultos

Precisión frente a Referencias Clásicas: El método Ext-SQD demuestra una alta precisión, con errores de energía del estado fundamental respecto a los benchmarks de Interacción de Configuración Completa (HCI) del orden de $10^{-8}$ eV para HfO $_2$ y $10^{-7}$ eV para ZrO $_2$ . El método supera los enfoques estándar de CCSD y CCSD(T) en los espacios de configuración probados.
Predicciones de la Brecha de Banda:
- Para HfO $_2$ , la brecha de banda computada cuánticamente cae dentro del rango de los valores experimentales reportados en la literatura y supera todos los niveles de simulación DFT (incluyendo DFT+U+V) y los resultados de la aproximación GW presentados para comparación.
- Para ZrO $_2$ , las simulaciones de DFT en todos sus niveles fallan en alcanzar el rango de la brecha de banda experimental. En contraste, las brechas de banda computadas cuánticamente logran alcanzar con éxito dicho rango experimental.
Rol de las Interacciones: El estudio destaca que, si bien las interacciones intra-sitio ( $U$ ) causan desplazamientos notables en la DFT, las interacciones inter-sitio ( $V$ ) juegan un papel dominante en la representación de red para capturar con precisión la brecha de banda.
Eficiencia: Se observa que el procedimiento de muestreo cuántico es relativamente rápido, produciendo resultados comparables a la costosa aproximación GW con menos recursos.

Significancia
El artículo afirma que este trabajo demuestra la utilidad de los flujos de trabajo híbridos cuánticos-clásicos en computadoras cuánticas pre-tolerantes a fallos para la ciencia de materiales aplicada. La principal significancia radica en la capacidad de estos métodos para producir predicciones de propiedades útiles y verificables experimentalmente (específicamente brechas de banda) que superan los enfoques de campo medio estándar (DFT) y se alinean con los experimentos de laboratorio. Los autores sugieren que este enfoque efectivamente muestrea las correlaciones electrónicas que son ignoradas por las teorías de campo medio. Aunque el cuello de botella computacional actual sigue siendo el cálculo clásico de la interacción de configuración, los resultados indican que, a medida que la calidad de los procesadores cuánticos mejore, estos flujos de trabajo ofrecerán un camino escalable para simular materiales periódicos complejos con alta precisión.

Más como este