Device variability of Josephson junctions induced by interface roughness

Resumen Técnico: Variabilidad de Dispositivos de Uniones Josephson Inducida por la Rugosidad de la Interfaz

Planteamiento del Problema
A medida que los procesadores cuánticos superconductores escalan hacia los miles de cúbits requeridos para la computación cuántica tolerante a fallos (por ejemplo, el algoritmo de Shor), la variabilidad de dispositivo a dispositivo ha emergido como un cuello de botella crítico. Los cúbits superconductores, particularmente los transmones, dependen de uniones Josephson de Al/AlOx/Al que operan en el régimen de tunelamiento. En este régimen, la corriente crítica $I_c$ —y, consecuentemente, la energía de Josephson $E_J$— depende exponencialmente de los detalles microscópicos de la barrera. Incluso variaciones menores en la geometría de la interfaz Al/AlOx, tales como fluctuaciones locales en el espesor o contaminación, pueden inducir desplazamientos sustanciales en las frecuencias de transición de los cúbits y en las fidelidades de las puertas. Si bien los esfuerzos experimentales extensos se han centrado en mejorar la reproducibilidad de las uniones mediante la optimización de procesos y la microscopía, se carecía de un marco teórico cuantitativo que conectara el desorden microestructural (específicamente la rugosidad de la interfaz) con las estadísticas de energía a nivel de unión.

Metodología
Los autores desarrollan un marco de modelado cuantitativo que vincula el desorden de la interfaz a escala atómica con las propiedades de transporte macroscópicas de la unión. El enfoque integra tres componentes principales:

1. Modelado de la Rugosidad de la Interfaz: La interfaz Al/AlOx se modela como un campo aleatorio gaussiano $h(x,y)$ caracterizado por dos parámetros: la amplitud de rugosidad raíz media cuadrática (RMS) $\sigma$ y la longitud de correlación transversal $\xi$. Estos parámetros se derivan de la literatura experimental existente (por ejemplo, estudios de TEM de tamaños de grano e histogramas de espesor) y de nuevas simulaciones de dinámica molecular (MD). Las simulaciones de MD, utilizando el potencial de aprendizaje automático CHGNet, modelan la oxidación de una superficie de Al para estimar $\sigma$ bajo condiciones idealizadas.
2. Aproximación de Transporte Cuántico: Para abordar la intratabilidad computacional de resolver la ecuación de Bogoliubov–de Gennes (BdG) para tamaños de unión realistas ($\sim 100 \text{ nm} \times 100 \text{ nm}$), los autores emplean dos aproximaciones clave:
   • Relación de Ambegaokar–Baratoff (AB): La corriente crítica se calcula mediante la relación AB, que vincula $I_c$ con la resistencia normal $R_N$. Esto reduce la complejidad computacional al evitar la integración sobre todo el espectro de energía y centrarse en el coeficiente de transmisión en la energía de Fermi.
   • Aproximación de Espesor Local: Dado que la longitud de onda de Fermi ($\lambda_F \approx 3.6$ Å) y la longitud de decaimiento de la barrera ($\lambda_D \approx 1.9$ Å) son mucho menores que la longitud de correlación de la rugosidad ($\xi \sim 10$ nm), la unión se trata como un conjunto de canales de tunelamiento independientes. La conductancia total se aproxima como la integral de las densidades de conductancia local, donde el espesor local de la barrera varía según el perfil de rugosidad.
3. Simulación Numérica: El modelo se aplica a un conjunto de 5,000 uniones Josephson con perfiles de barrera de túnel variables generados con el modelo de campo aleatorio gaussiano.

Resultados Clave
Las simulaciones numéricas arrojan varias perspectivas cuantitativas sobre el comportamiento estadístico de $E_J$:

• Distribución Log-Normal: La distribución de la energía de Josephson $E_J$ a través del conjunto no es gaussiana, sino que sigue una distribución log-normal. Esto surge porque $E_J$ depende exponencialmente del espesor de la barrera, el cual está distribuido normalmente. Consecuentemente, la función de densidad de probabilidad está sesgada hacia valores mayores.
• Impacto de la Rugosidad RMS ($\sigma$): El aumento de la amplitud de rugosidad RMS $\sigma$ conduce a un incremento significativo tanto en el valor medio $\langle E_J \rangle$ como en la varianza $\text{Var}(E_J)$. El aumento en la media se atribuye a la sensibilidad exponencial del tunelamiento, donde regiones anomalamente delgadas y raras dominan la conductancia total.
• Impacto de la Longitud de Correlación ($\xi$): El aumento de la longitud de correlación transversal $\xi$ causa una ligera disminución en el valor medio $\langle E_J \rangle$, pero un incremento marcado en la varianza $\text{Var}(E_J)$. La varianza aumenta porque un $\xi$ mayor reduce el número de canales de tunelamiento estadísticamente independientes (auto-promediado), potenciando así las fluctuaciones de muestra a muestra. La media disminuye ligeramente debido a efectos de tamaño finito; a medida que $\xi$ se aproxima a las dimensiones de la unión, la probabilidad de realizar las regiones raras de barrera extremadamente delgada que dominan el tunelamiento se reduce.
• Variabilidad de la Frecuencia del Cúbit: Utilizando las estadísticas derivadas para $E_J$ y asumiendo una energía de carga $E_C$ típica, el modelo predice una desviación estándar en la frecuencia de transición del cúbit $\omega_{01}/2\pi$ de aproximadamente 0.364 GHz alrededor de una media de 6.145 GHz para los parámetros específicos probados.

Significancia y Reivindicaciones
El artículo afirma proporcionar una "imagen cuantitativa e intuitiva" de cómo la rugosidad superficial induce la variabilidad en la energía de Josephson. Al conectar los parámetros microestructurales ($\sigma$ y $\xi$) directamente con las estadísticas de la unión, el trabajo ofrece un marco para comprender los desafíos de rendimiento y uniformidad en procesadores cuánticos a gran escala. Los autores enfatizan que sus resultados tienen "relevancia directa para el diseño de uniones", sugiriendo que controlar los parámetros de rugosidad de la interfaz es crucial para minimizar la variabilidad de los dispositivos.

Los autores mantienen una postura modesta respecto al alcance de su trabajo. Señalan que su modelo asume que $\sigma$ y $\xi$ son independientes e idénticos para ambas interfaces en la unión, un supuesto que requiere verificación experimental bajo diversas condiciones de crecimiento. Además, identifican direcciones futuras para conectar este modelo con la caracterización experimental de uniones fabricadas e integrar sus resolvedores de superconductividad con software de ingeniería de microondas para permitir el co-diseño de uniones y elementos de circuito. El artículo no propone nuevos protocolos experimentales ni pretende resolver el problema de la variabilidad por completo, sino que establece una línea base teórica para cuantificar el impacto de la rugosidad de la interfaz.