Classification of singular limits for free boundary and singularly perturbed elliptic problems: the Dancer-Yan spikes revisited

Planteamiento del Problema
Este artículo aborda la clasificación de los límites singulares para un problema de frontera libre que surge en la física de plasmas, modelando específicamente los equilibrios de plasma en una sección transversal de un Tokamak. El problema se formula como una ecuación elíptica singularmente perturbada en dimensión $d=2$:
$$\begin{cases} -\Delta v = [v]_+^p & \text{en } \Omega, \\ v = \gamma & \text{en } \partial\Omega, \\ \int_\Omega [v]_+^p = I, \end{cases}$$
donde $\Omega \subset \mathbb{R}^2$ es un dominio acotado $C^{2,\beta}$, $p \in (1, +\infty)$, e $I > 0$ representa una corriente total fija. Los autores investigan el comportamiento asintótico de las soluciones cuando la corriente total $I \to +\infty$. En este límite, la "región de plasma" $\Omega_+ = \{v > 0\}$ se espera que desarrolle una estructura no trivial, potencialmente concentrándose en "picos" (regiones donde $v$ es significativamente mayor que el valor en la frontera).

Si bien el comportamiento de tales problemas en dimensiones superiores ($d \ge 3$) ha sido establecido como una superposición de finitely many "spikes" (específicamente, picos de Dancer-Yan), el caso bidimensional presenta dificultades únicas. En $d=2$, el problema global límite estándar $-\Delta w = [w-1]_+^p$ en $\mathbb{R}^2$ con masa finita no admite soluciones, lo que impide la aplicación directa de argumentos de explosión (blow-up) de dimensiones superiores.

Metodología
Los autores emplean un análisis de explosión refinado adaptado al entorno bidimensional, basándose en el marco de Dancer y Yan pero introduciendo nuevas técnicas de reescalado para manejar la falta de soluciones globales estándar.

1. Transformación de Variables: El problema se reformula utilizando las variables $(\alpha, \psi)$ donde $v = \lambda^{1/(p-1)}(\alpha + \lambda \psi)$ y $\lambda = I^{1/q}$ ($\lambda$ siendo $q$ el conjugado de Hölder de $p$). Esto transforma el problema en una forma singularmente perturbada que involucra un parámetro pequeño $\varepsilon_n \to 0$.
2. Cotas Integrales: El análisis se apoya en dos cotas integrales naturales, denotadas como (H1) y (H2), que controlan la $(p-1)$-masa y la $p$-masa de las soluciones, respectivamente. Estas cotas son cruciales para descartar soluciones de "masa infinita" y asegurar la compacidad.
3. Reescalado Refinado: A diferencia del caso $d \ge 3$ donde los máximos locales permanecen alejados de 1, en $d=2$ los máximos $v_n(x_n)$ convergen a 1. Los autores introducen un reescalado delicado y multiescala involucrando parámetros $s_n$ y $\theta_n$ (derivados de la solución de la ecuación de Emden $\phi$) para capturar la estructura "microscópica" de los picos.
   • Definen una función reescalada $\tilde{v}_n$ y una función normalizada adicional $u_n$ para analizar el perfil cerca de los maximizadores locales.
4. Clasificación de Soluciones Enteras: Un paso técnico clave es la clasificación de soluciones de masa finita para la ecuación plana $-\Delta w = [w]_+^p$ en $\mathbb{R}^2$. Los autores demuestran que tales soluciones son o bien triviales o bien radiales y toman una forma específica que involucra la solución de Emden $\phi$ dentro de un disco y un decaimiento logarítmico fuera de él.
5. Extracción de Múltiples Picos: Para manejar múltiples picos agrupándose en un mismo punto, los autores utilizan un argumento inductivo. Demuestran que dos picos distintos no pueden estar arbitrariamente cerca en relación con sus parámetros de escala y que la presencia de un pico de Tipo I o Tipo II impide la existencia de un pico de "Desvanecimiento" (Fading) en la misma ubicación.
6. Identidad de Pohozaev: Para la estructura global, los autores utilizan la identidad de Pohozaev para derivar una restricción sobre las ubicaciones de los picos, mostrando que estos deben ser puntos críticos de un Hamiltoniano de Kirchhoff-Routh.

Contribuciones Clave y Resultados

1. Clasificación de Límites Singulares (Teorema 1.1):
   Los autores establecen que, bajo las cotas integrales (H1) y (H2), el límite singular de las soluciones cae en una de tres categorías alrededor de cualquier secuencia de maximizadores locales:

   • Desvanecimiento (Vanishing): La solución permanece por debajo del umbral 1 en subconjuntos compactos.
   • Picos de Tipo I: La solución converge a un perfil de "pico de Dancer-Yan". Aquí, la función reescalada converge a la solución entera estándar $w^*$, y la masa está cuantizada.
   • Picos de Tipo II: Un fenómeno nuevo específico de $d=2$. La función reescalada aún converge a $w^*$, pero los parámetros de escala divergen de tal manera que la $p$-masa se anula en el límite, mientras que la $(p-1)$-masa permanece cuantizada.
   • Picos de Desvanecimiento (Fading Spikes): La amplitud del pico decae tan rápidamente que efectivamente desaparece en el límite, con la solución comportándose como una función armónica perturbada por un término pequeño.

   Crucialmente, el artículo demuestra que, a diferencia de las dimensiones superiores, no toda solución es una simple superposición de picos de Dancer-Yan. La presencia de picos de Tipo II y de Desvanecimiento indica una estructura singular más rica en 2D.

2. Estructura Global con Condiciones de Dirichlet (Teorema 1.2):
   Al imponer condiciones de frontera de Dirichlet ($v=0$ en $\partial\Omega$) y una condición de no anulación (NVp) sobre la $p$-masa, los autores proporcionan una clasificación detallada para la solución global, mientras señalan que una descripción completa de la interacción entre diferentes tipos de picos agrupándose en el mismo punto sigue siendo un problema abierto.

   • El conjunto singular $\Sigma$ consiste en un número finito de puntos interiores.
   • La solución es una superposición de picos de Tipo I y Tipo II centrados en estos puntos.
   • La $(p-1)$-masa está cuantizada: es igual a $(N_I + N_{II})I_{p-1}$, donde $N_I$ y $N_{II}$ son los conteos de picos de Tipo I y Tipo II.
   • La $p$-masa no está cuantizada en general; depende de los límites de escala específicos de los picos de Tipo I.
   • La región de plasma $\Omega_{n,+}$ consiste en componentes asintóticamente redondas alrededor de los centros de los picos.
   • Las ubicaciones de los picos de Tipo I son puntos críticos de un Hamiltoniano de Kirchhoff-Routh que involucra la función de Green del dominio.
   • Nota: Aunque los autores logran descartar los picos de "Desvanecimiento" en presencia de picos de Tipo I o Tipo II bajo estas condiciones, declaran explícitamente que una descripción completa de la interacción local cuando múltiples picos de diferentes tipos (Tipo I y Tipo II) se agrupan en el mismo punto es actualmente un problema abierto.

3. Aplicación al Problema Original del Plasma (Teorema 1.3):
   Los autores traducen estos resultados de vuelta a las variables originales $(v, \gamma)$ para el problema del plasma. Muestran que, bajo restricciones asintóticas específicas de los parámetros $\lambda_n$ y $\alpha_n$, las soluciones exhiben las estructuras de picos clasificadas. Proporcionan fórmulas asintóticas precisas para el valor en la frontera $\gamma$ y la corriente total en términos del número y tipo de picos.

Significancia y Reivindicaciones
El artículo afirma proporcionar la primera descripción detallada del límite singular para este problema de frontera libre en dimensión $d=2$. Su significancia principal radica en revelar que el caso bidimensional es fundamentalmente diferente al caso de dimensiones superiores:

• Riqueza de Estructura: La existencia de picos de Tipo II y de Desvanecimiento muestra que el espacio de soluciones es más complejo que una simple unión de perfiles estándar.
• Fenómeno de Cuantización: El trabajo destaca un comportamiento de cuantización distinto donde la $(p-1)$-masa está cuantizada (similar a las ecuaciones de tipo Liouville), pero la $p$-masa no lo está, desafiando la intuición heredada de $d \ge 3$.
• Avance Metodológico: El trabajo introduce un refinamiento necesario en las técnicas de reescalado para manejar el desvanecimiento del valor máximo en 2D, superando el obstáculo de que el problema global límite estándar no posee soluciones.

Los autores declaran explícitamente que una descripción completa de la interacción entre múltiples picos de diferentes tipos agrupándose en el mismo punto sigue siendo un problema abierto, aunque logran descartar ciertas configuraciones (como la coexistencia de picos de Desvanecimiento con picos de Tipo I/II) bajo condiciones de Dirichlet. También conjeturan que los picos de Tipo II podrían no existir en dominios convexos, sugiriendo una dependencia de la geometría del dominio.