Analytical Charge Density Profile of Vortex Core in Weak-Coupling Superconductor

Resumen Técnico: Perfil Analítico de la Densidad de Carga del Núcleo de un Vórtice en un Superconductor de Acoplamiento Débil

Planteamiento del Problema
Se sabe desde hace tiempo que la estructura electrónica de los núcleos de vórtices en superconductores de tipo II alberga estados ligados de cuasipartículas discretos (estados de Caroli–de Gennes–Matricon o CdGM). Si bien los cálculos de Bogoliubov–de Gennes (BdG) auto-consistentes acoplados con la ecuación de Poisson han demostrado numéricamente que la densidad de carga dentro de un núcleo de vórtice exhibe oscilaciones de signo alternante con un periodo de $\pi/k_F$, se carecía de una explicación analítica sencilla para este fenómeno. Estudios numéricos previos (p. ej., Machida y Koyama) observaron que la densidad de carga apantallada cambia de signo, pero el mecanismo que impulsa esta reversión de signo —específicamente la interacción entre la distribución de carga de los estados ligados y el apantallamiento electrostático— no había sido derivado desde los primeros principios de forma cerrada. El desafío radica en la complejidad de las ecuaciones de BdG, donde las constantes de normalización son típicamente dependientes del modo y las funciones de onda involucran ecuaciones diferenciales no lineales acopladas, lo que dificulta la sumatoria analítica.

Metodología
El autor emplea un enfoque elemental y explicativo utilizando el modelo más económico que retiene la física esencial del núcleo del vórtice:

1. Sistema Modelo: Un único vórtice lineal en un superconductor de onda-$s$ de acoplamiento débil con un perfil de brecha escalonado ($\Delta(r) = 0$ para $r < \xi$ y $\Delta(r) = \Delta_\infty$ en otro caso). El espectro se toma como el resultado estándar de CdGM, $E_\mu = \mu \Delta^2/E_F$.
2. Construcción de la Función de Onda:
   • Región Interna ($r < \xi$): La brecha desaparece, desacoplando las ecuaciones de BdG. El autor utiliza soluciones exactas de funciones de Bessel ($J_{\mu \mp 1/2}$) en lugar del ansatz WKB/BKJT, el cual es mal definido cerca de la línea del vórtice debido al enrollamiento topológico.
   • Región Externa ($r > \xi$): El envolvente BKJT se linealiza, con coeficientes fijados por las condiciones de acoplamiento en el borde del núcleo.
3. Análisis de Normalización: Un paso crítico consiste en calcular la constante de normalización $C_\mu$. El autor demuestra que, a pesar de la variable ventana clásicamente permitida y las tasas de decaimiento, la constante de normalización es casi independiente del número cuántico de momento angular $\mu$ para la mayor parte de la escalera de CdGM. Esto permite tratar a $C_\mu$ como un factor prefactor constante $\langle C^2 \rangle$ en la suma de modos.
4. Sumatoria de Modos:
   • Dentro del Núcleo: La suma de la densidad de carga se reduce a una suma parcial de Bessel–Parseval. Al aplicar la identidad de completitud para las funciones de Bessel y notar que la restricción topológica excluye el canal $\nu=0$, la suma colapsa en una forma cerrada.
   • Fuera del Núcleo: La suma produce un envolvente decreciente de $1/r$ con un rizo residual.
5. Electrostática: La densidad de carga de los estados ligados sirve como término fuente para la ecuación de Poisson. El autor resuelve la ecuación de Poisson apantallada en el espacio de momentos, tratando por separado el fondo suave y los componentes oscilatorios de $2k_F$.

Contribuciones Clave y Resultados

• Densidad Interna de Forma Cerrada: El autor deriva una expresión analítica para la densidad de carga de los estados ligados dentro del núcleo: $\delta n(r < \xi) \propto 1 - J_0^2(k_F r)$. Este resultado satisface la condición de contorno $n(0)=0$ (debido a la exclusión del canal de Bessel $\nu=0$) y revela que la densidad consiste en un fondo suave más oscilaciones de tipo Friedel con vector de onda $2k_F$.
• Mecanismo de Apantallamiento: El artículo elucida por qué la densidad de carga apantallada alterna en signo. La carga de los estados ligados no integra a cero; por lo tanto, la neutralidad de carga global requiere que los estados extendidos compensen.
  • El apantallamiento de Thomas–Fermi es altamente eficiente para longitudes de onda largas (transferencia de momento baja), neutralizando completamente el componente de fondo suave.
  • Sin embargo, el apantallamiento es ineficiente para transferencias de momento grandes ($q \approx 2k_F$). Consecuentemente, el componente oscilatorio de $2k_F$ de la carga de los estados ligados sobrevive al proceso de apantallamiento.
• Expresión Analítica de la Carga Sobreviviente: La densidad de carga apantallada final viene dada por:
  $$\delta n_{\text{total}}(r) \simeq -\langle C^2 \rangle \frac{4k_F^2}{4k_F^2 + k_{TF}^2} \frac{\sin(2k_F r)}{\pi k_F r}$$
  Esta expresión muestra que la oscilación sobreviviente es una modulación de $2k_F$ de signo alternante. La amplitud se reduce por un factor de $4k_F^2/(4k_F^2 + k_{TF}^2)$, que se sitúa entre $1/2$ y $3/4$ para metales típicos.
• Validación: Los resultados analíticos se verifican frente a soluciones numéricas directas de la ecuación de Poisson utilizando la suma completa de modos de CdGM. Las formas cerradas analíticas concuerdan con los datos numéricos con un error inferior al 1% (el error residual se atribuye a la aproximación de promedio de modo de la constante de normalización).

Significancia y Reivindicaciones
El artículo afirma proporcionar un relato elemental del perfil de densidad de carga con cambio de signo observado en simulaciones numéricas. Su principal significancia radica en:

1. Claridad Mecanística: Identifica que la oscilación no es un efecto delicado, sino una consecuencia directa de la neutralidad de carga y la ineficiencia del apantallamiento de Thomas–Fermi ante transferencias de momento grandes ($2k_F$).
2. Referencia Analítica: Ofrece un referente analítico de forma cerrada (Ecs. 18, 22 y 28) contra el cual se pueden verificar futuros cálculos numéricos auto-consistentes, particularmente en el límite de acoplamiento débil donde el perfil de la brecha es rígido.
3. Distinción de la Dispersión por Impurezas: El autor aclara que, si bien la periodicidad ($2k_F$) es la misma que las oscilaciones de Friedel alrededor de impurezas, el mecanismo es inverso. En la dispersión por impurezas, la fuente es sin estructura y la oscilación surge de la no-analyticidad de la función de respuesta. En el núcleo del vórtice, la fuente misma está centrada en $2k_F$, y la respuesta (apantallamiento) simplemente filtra los componentes suaves.

El autor establece explícitamente que su objetivo es explicativo más que predictivo, y no pretende resolver detalles cuantitativos en el régimen de acoplamiento fuerte o para simetrías de brecha complejas, aunque señala estas como extensiones naturales. El trabajo confirma que los "rizos" en la densidad de carga son rasgos robustos de los núcleos de vórtices en superconductores de acoplamiento débil.