Pressure and Size Dependence of Roton Emission and Vortex Creation by Moving Objects in He~II in $T \to 0$ Limit: Generalized Nonlocal Gross-Pitaevskii Model

Este trabajo presenta un estudio numérico mediante un modelo generalizado no local de Gross-Pitaevskii para analizar cómo la presión y el tamaño del obstáculo influyen en las velocidades críticas para la emisión de rotones y la nucleación de vórtices por objetos en movimiento en helio-4 superfluido a temperatura cero, marcando el primer marco teórico que aborda simultáneamente estos mecanismos dependientes de la presión.

Lo esencial

Imagina un mundo donde un líquido puede fluir sin ninguna fricción en absoluto, como un fantasma deslizándose entre tus dedos. Este es el Helio Superfluido (He II), un estado especial de la materia que existe únicamente cuando el helio se enfría a temperaturas cercanas al cero absoluto.

Este artículo es como un laboratorio de simulación de alta tecnología donde los autores someten a prueba este líquido fantasmal. Querían entender qué sucede cuando empujas un objeto pequeño (como una burbuja diminuta o una partícula cargada) a través de este superfluido. Específicamente, querían saber: ¿A qué velocidad puedes empujarlo antes de que el líquido "rompa" su flujo perfecto?

Aquí está la historia de su descubrimiento, desglosada en conceptos simples:

1. Las Dos Maneras en que el Líquido se Rompe

Cuando empujas un objeto a través del helio superfluido, el líquido no simplemente se "arrastra" como el agua. En cambio, reacciona de dos maneras distintas una vez que alcanzas cierto límite de velocidad (llamado velocidad crítica):

• La Explosión de "Rotones": Imagina el líquido como un océano tranquilo. Si empujas el objeto demasiado rápido, no solo haces olas; de repente creas una enjambre de partículas diminutas y energéticas llamadas rotones. Es como si el líquido decidiera de repente fragmentarse en un millón de chispas diminutas y energéticas. Esto ocurre a una velocidad específica.
• El Remolino de "Vórtices": Si lo empujas aún más rápido (o si el objeto es lo suficientemente grande), el líquido comienza a girar. Crea pequeños tornados microscópicos llamados vórtices cuánticos. Estos son como pequeños remolinos que se quedan atrapados en el objeto, arrastrándolo hacia abajo.

El objetivo principal del artículo era determinar exactamente qué tan rápido debes ir para activar las "chispas" (rotones) en comparación con los "remolinos" (vórtices).

2. El Experimento de la Olla a Presión

Los autores no solo observaron el líquido a una sola presión. Simularon lo que sucede a medida que aprietan el helio más y más, desde un vacío (0 bares) hasta el punto donde se convertiría en una roca sólida (aproximadamente 25 bares).

Utilizaron un modelo matemático especial (un "Modelo Generalizado No Local de Gross-Pitaevskii") que actúa como un motor de videojuegos superpreciso. Este motor fue programado para imitar el comportamiento real y complejo de los átomos de helio, incluido el extraño comportamiento de "rotones" que las ecuaciones físicas estándar suelen pasar por alto.

3. El Gran Descubrimiento: Apretar Cambia las Reglas

Aquí está lo que encontraron, usando una analogía simple:

Imagina que intentas correr a través de una multitud de personas.

• A baja presión (multitud suelta): La gente está dispersa. En realidad es bastante difícil iniciar un motín (crear un roton) porque están lejos entre sí. Pero si corres lo suficientemente rápido, podrías tropezar con alguien e iniciar una reacción en cadena de personas cayendo (vórtices).
• A alta presión (multitud apretada): La gente está apretada hombro con hombro. Ahora, es mucho más fácil iniciar un motín (rotones) porque están tan cerca entre sí. Sin embargo, se vuelve más difícil iniciar una reacción en cadena de personas cayendo (vórtices) porque la multitud es tan densa y rígida que resiste girar.

Los Resultados:

• Velocidad de Rotón: A medida que apretaron el helio (aumentaron la presión), la velocidad necesaria para crear esas "chispas" (rotones) disminuyó. No tienes que correr tan rápido para romper el flujo.
• Velocidad de Vórtice: A medida que apretaron el helio, la velocidad necesaria para crear los "remolinos" (vórtices) aumentó. Tienes que correr mucho más rápido para hacer que el líquido gire.

4. El "Punto Dulce" para la Detección

Esto crea una brecha fascinante. A altas presiones, hay un amplio rango de velocidades donde puedes crear las "chispas" (rotones) sin crear los "remolinos" (vórtices).

En el pasado, los científicos lucharon por estudiar los rotones porque a menudo estaban ocultos por los desordenados remolinos. Los autores sugieren que, al apretar el helio a altas presiones, podemos crear un entorno "limpio" donde los rotones aparecen por sí solos, haciéndolos mucho más fáciles de estudiar.

5. El Tamaño Importa

El artículo también examinó el tamaño del objeto que se mueve a través del líquido.

• Objetos diminutos (como un solo ion): Son muy sensibles. Alcanzan primero el "límite de roton".
• Objetos grandes (como un disco grande): Son menos sensibles a los rotones. Tienden a alcanzar primero el "límite de vórtice", independientemente de la presión.

Resumen

Los autores construyeron un microscopio digital para observar el helio superfluido bajo presión. Descubrieron que apretar el helio hace que sea más fácil crear chispas de energía (rotones) pero más difícil crear remolinos giratorios (vórtices).

Esto explica por qué los experimentos del pasado observaron comportamientos diferentes a diferentes presiones y sugiere que, si queremos estudiar las misteriosas partículas "rotones", debemos realizar nuestros experimentos bajo alta presión, donde es más probable que el líquido revele sus secretos sin ensuciarse con remolinos.

Nota: Los autores admiten que su simulación se realizó en dos dimensiones (una rebanada plana del mundo) porque hacerlo en 3D completo es demasiado pesado computacionalmente, pero creen que la física que encontraron se mantiene verdadera para el mundo real, en 3D.

Resumen técnico

Resumen Técnico: Dependencia de la Presión y el Tamaño en la Emisión de Rotones y la Creación de Vórtices por Objetos en Movimiento en He II en el Límite T → 0

Planteamiento del Problema
La ruptura de la superfluidez en Helio II (He II) está gobernada por dos mecanismos principales cuando un objeto se mueve a través del fluido: la emisión de rotones (excitaciones elementales) y la nucleación de vórtices cuantizados. Si bien la dependencia de la presión de la velocidad crítica para la emisión de rotones ($v_{rot}^{crit}$) y la creación de vórtices ($v_{v}^{crit}$) ha sido investigada experimentalmente —particularmente por el grupo de Lancaster utilizando iones negativos—, estos fenómenos han carecido históricamente de una descripción teórica unificada. Las observaciones experimentales indican que a bajas presiones, la nucleación de vórtices limita la velocidad del ion, whereas a altas presiones (por encima de $\approx 10-12$ bar), la emisión de rotones se convierte en el factor limitante. Además, la velocidad crítica para la emisión de rotones disminuye con la presión, mientras que la velocidad crítica para la nucleación de vórtices aumenta. Los marcos teóricos existentes, como la ecuación estándar de Gross-Pitaevskii (GPE), no logran capturar la emisión de rotones porque poseen un espectro de excitación estrictamente monótono, incapaz de reproducir la curva de dispersión de Landau con su característico mínimo "roton". Por el contrario, las teorías microscópicas a menudo luchan por describir ambos mecanismos dentro de un único marco.

Metodología
Los autores emplean un modelo generalizado y no local de Gross-Pitaevskii (gnlGPE) para simular la dinámica del $^4$He superfluido a temperatura cero ($T \to 0$). El modelo se define mediante un parámetro de orden complejo $\psi(x,t)$ e incorpora un potencial de interacción interatómica no local $V(r)$ y un término de corrección de densidad de orden superior.

1. Calibración del Modelo: Para garantizar la precisión física en un rango de presiones (de 0 a 24 bar, acercándose a la presión de solidificación de $\approx 25$ bar), los autores no utilizan una forma analítica fija para el potencial de interacción. En su lugar, invierten la relación de dispersión de Bogoliubov generalizada utilizando datos experimentales de Godfrin et al. (2021). Esto permite que el modelo reproduzca las curvas de dispersión dependientes de la presión exactas, incluyendo las ramas de fonón, maxón y roton, así como la meseta de Pitaevskii.
2. Configuración de la Simulación: El estudio se realiza principalmente en una aproximación bidimensional (2D). Los autores simulan un obstáculo circular (disco) moviéndose a través de He II en reposo. El obstáculo se modela como un potencial externo con un perfil gaussiano que agota completamente la densidad del superfluido.
3. Determinación de la Velocidad Crítica:
   • Emisión de Rotones ($v_{rot}^{crit}$): Determinada analizando el diagrama de bifurcación de soluciones estacionarias. Los autores utilizan un método de Newton-Raphson acoplado con un método de gradiente conjugado bi-estabilizado para encontrar el número de Mach máximo ($M = U/c$) para el cual existen soluciones estacionarias estables. Más allá de este punto, el sistema se vuelve inestable y se emiten rotones.
   • Nucleación de Vórtices ($v_{v}^{crit}$): Determinada mediante un enfoque dinámico. El obstáculo se acelera hasta alcanzar una velocidad objetivo y el sistema se evoluciona durante una escala de tiempo de $\approx 10\xi/c$ (donde $\xi$ es la longitud de curación) para observar si se nuclean vórtices cuantizados.
4. Parámetros: Las simulaciones cubren presiones $P = 0, 2, 5, 10,$ y $24$ bar. Los tamaños de los obstáculos varían desde $D = 10\xi$ hasta $D = 1000\xi$. Se presta atención específica al tamaño de los iones negativos, que varía con la presión (desde $\approx 34.5$ Å a 0 bar hasta $\approx 21.7$ Å a 24 bar).

Resultados Clave

• Dependencia de la Presión de las Velocidades Críticas: Las simulaciones confirman las tendencias opuestas observadas experimentalmente. A medida que aumenta la presión, el número de Mach crítico para la emisión de rotones ($M_{rot}^{crit}$) disminuye, cayendo de $\approx 0.247$ a 0 bar hasta $\approx 0.128$ a 24 bar. Por el contrario, el número de Mach crítico para la nucleación de vórtices ($M_{v}^{crit}$) aumenta con la presión.
• Dependencia del Tamaño: Para obstáculos pequeños (comparables a los tamaños de iones negativos, $D < 100$ nm), la velocidad crítica está determinada por el umbral más bajo de los dos ($v_{rot}^{crit}$ o $v_{v}^{crit}$). A bajas presiones, $v_{v}^{crit}$ es el factor limitante; a altas presiones, $v_{rot}^{crit}$ es el limitante. A medida que aumenta el tamaño del obstáculo ($D \ge 1000\xi$), las velocidades críticas para todas las presiones colapsan en una sola curva gobernada por la nucleación de vórtices, similar a las predicciones de la GPE estándar.
• La "Brecha": Un hallazgo significativo es el ensanchamiento de la brecha entre $v_{rot}^{crit}$ y $v_{v}^{crit}$ a medida que aumenta la presión. A altas presiones, existe una ventana de velocidad donde se emiten rotones, pero la nucleación de vórtices aún no ha ocurrido. Esto sugiere que la emisión de rotones puede aislarse y observarse experimentalmente a altas presiones sin la complicación inmediata de la generación de vórtices.
• Visualización: Las visualizaciones de densidad confirman que a números de Mach intermedios (por ejemplo, $M=0.28$), se emiten rotones en todas las presiones, mientras que la nucleación de vórtices solo ocurre a números de Mach más altos ($M \approx 0.3$), con el umbral para la nucleación desplazándose a valores más altos a medida que aumenta la presión.

Significado y Afirmaciones
Los autores afirman que este trabajo representa el primer intento de analizar la dependencia de la presión de la emisión de rotones y la nucleación de vórtices, junto con sus velocidades críticas asociadas, dentro de un único marco teórico. Al utilizar la GPE generalizada y no local, el estudio logra cerrar la brecha entre el criterio de Landau (emisión de rotones) y la nucleación hidrodinámica de vórtices.

El artículo enfatiza que, aunque el estudio se restringe a 2D debido a la complejidad computacional de las simulaciones 3D (particularmente en lo que respecta a la determinación de la semilla mínima para la nucleación de vórtices en 3D), los resultados 2D proporcionan una explicación cualitativa y semicuantitativa de la física subyacente que se alinea con las observaciones experimentales. Los autores señalan que sus valores numéricos para las velocidades críticas son comparables con los datos experimentales, validando el enfoque de utilizar relaciones de dispersión derivadas experimentalmente para restringir el potencial de interacción. El trabajo ofrece una explicación motivada físicamente de por qué el carácter del movimiento del ion cambia en la presión crítica $P_{crit} \approx 10-12$ bar, atribuyéndolo al cruce de las curvas dependientes de la presión para $v_{rot}^{crit}$ y $v_{v}^{crit}$.

Los autores mantienen una postura modesta respecto a las implicaciones en 3D, reconociendo que, si bien las simulaciones 2D explican las tendencias generales, la topología específica de la nucleación de vórtices en 3D (por ejemplo, bucles ion-vórtice frente a anillos) sigue siendo una pregunta abierta que requiere investigación futura. También señalan que modelar el ion como una esfera dura perfecta es una simplificación, aunque se espera que las tendencias de presión permanezcan robustas.