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 artigo apresenta um estudo numérico utilizando um modelo generalizado não local de Gross-Pitaevskii para analisar como a pressão e o tamanho do obstáculo influenciam as velocidades críticas para a emissão de rotons e a nucleação de vórtices por objetos em movimento no hélio-4 superfluido a temperatura zero, marcando a primeira estrutura teórica a abordar simultaneamente esses mecanismos dependentes da pressão.

O essencial

Imagine um mundo onde um líquido pode fluir sem qualquer atrito, como um fantasma deslizando entre os seus dedos. Este é o Hélio Superfluido (He II), um estado especial da matéria que existe apenas quando o hélio é resfriado a temperaturas próximas do zero absoluto.

Este artigo é como um laboratório de simulação de alta tecnologia onde os autores colocam esse líquido fantasmagórico à prova. Eles queriam entender o que acontece quando você empurra um pequeno objeto (como uma bolha minúscula ou uma partícula carregada) através desse superfluido. Especificamente, eles queriam saber: Qual é a velocidade máxima que você pode atingir antes que o líquido "quebre" seu fluxo perfeito?

Aqui está a história de sua descoberta, decomposta em conceitos simples:

1. As Duas Maneiras pelas Quais o Líquido Quebra

Quando você empurra um objeto através do hélio superfluido, o líquido não é simplesmente "arrastado" como a água. Em vez disso, ele reage de duas maneiras distintas assim que você atinge um certo limite de velocidade (chamado de velocidade crítica):

• A Explosão de "Rótons": Imagine o líquido como um oceano calmo. Se você empurrar o objeto muito rápido, não cria apenas ondas; você cria repentinamente um enxame de partículas minúsculas e energéticas chamadas rótons. É como se o líquido decidisse subitamente se fragmentar em um milhão de faíscas minúsculas e energéticas. Isso ocorre em uma velocidade específica.
• O Redemoinho de "Vórtices": Se você empurrá-lo ainda mais rápido (ou se o objeto for grande o suficiente), o líquido começa a girar. Ele cria minúsculos tornados microscópicos chamados vórtices quânticos. Estes são como pequenos redemoinhos que ficam presos ao objeto, arrastando-o para baixo.

O objetivo principal do artigo era descobrir exatamente quão rápido você precisa ir para disparar as "faíscas" (rótons) versus os "redemoinhos" (vórtices).

2. O Experimento da Panela de Pressão

Os autores não observaram o líquido apenas em uma pressão. Eles simularam o que acontece à medida que apertam o hélio cada vez mais, desde o vácuo (0 bar) até o ponto em que ele se transformaria em uma rocha sólida (cerca de 25 bar).

Eles usaram um modelo matemático especial (um "Modelo de Gross-Pitaevskii Generalizado Não Local") que atua como um motor de videogame superpreciso. Este motor foi programado para imitar o comportamento real e complexo dos átomos de hélio, incluindo o estranho comportamento de "róton" que as equações físicas padrão geralmente ignoram.

3. A Grande Descoberta: Apertar Muda as Regras

Aqui está o que eles descobriram, usando uma analogia simples:

Imagine que você está tentando correr através de uma multidão de pessoas.

• Em baixa pressão (multidão solta): As pessoas estão espalhadas. Na verdade, é bastante difícil iniciar um motim (criar um róton) porque elas estão distantes. Mas, se você correr rápido o suficiente, pode tropeçar em alguém e iniciar uma reação em cadeia de pessoas caindo (vórtices).
• Em alta pressão (multidão apertada): As pessoas estão apertadas ombro a ombro. Agora, é muito mais fácil iniciar um motim (rótons) porque elas estão tão próximas. No entanto, torna-se mais difícil iniciar uma reação em cadeia de pessoas caindo (vórtices) porque a multidão é tão densa e rígida que resiste a girar.

Os Resultados:

• Velocidade do Róton: À medida que apertavam o hélio (aumentando a pressão), a velocidade necessária para criar essas "faíscas" (rótons) diminuiu. Você não precisa correr tão rápido para quebrar o fluxo.
• Velocidade do Vórtice: À medida que apertavam o hélio, a velocidade necessária para criar os "redemoinhos" (vórtices) aumentou. Você precisa correr muito mais rápido para fazer o líquido girar.

4. O "Ponto Ideal" para Detecção

Isso cria uma lacuna fascinante. Em altas pressões, há uma ampla faixa de velocidades onde você pode criar as "faíscas" (rótons) sem criar os "redemoinhos" (vórtices).

No passado, os cientistas lutavam para estudar rótons porque eles eram frequentemente escondidos pelos redemoinhos bagunçados. Os autores sugerem que, ao apertar o hélio até altas pressões, podemos criar um ambiente "limpo" onde os rótons aparecem sozinhos, tornando-os muito mais fáceis de estudar.

5. O Tamanho Importa

O artigo também analisou o tamanho do objeto movendo-se através do líquido.

• Objetos minúsculos (como um único íon): Eles são muito sensíveis. Eles atingem o "limite de róton" primeiro.
• Objetos grandes (como um disco grande): Eles são menos sensíveis aos rótons. Eles tendem a atingir o "limite de vórtice" primeiro, independentemente da pressão.

Resumo

Os autores construíram um microscópio digital para observar o hélio superfluido sob pressão. Eles descobriram que apertar o hélio torna mais fácil criar faíscas de energia (rótons), mas mais difícil criar redemoinhos giratórios (vórtices).

Isso explica por que experimentos no passado observaram comportamentos diferentes em pressões diferentes e sugere que, se quisermos estudar as misteriosas partículas "róton", devemos realizar nossos experimentos sob alta pressão, onde o líquido é mais propenso a revelar seus segredos sem se tornar bagunçado com redemoinhos.

Nota: Os autores admitem que sua simulação foi feita em duas dimensões (uma fatia plana do mundo) porque fazê-lo em 3D completo é computacionalmente pesado demais, mas eles acreditam que a física que descobriram é válida para o mundo real, tridimensional.

Resumo técnico

Resumo Técnico: Dependência da Pressão e do Tamanho na Emissão de Rótons e Criação de Vórtices por Objetos em Movimento no He II no Limite T → 0

Enunciado do Problema
A quebra da superfluidez no Hélio II (He II) é governada por dois mecanismos primários quando um objeto se move através do fluido: a emissão de rótons (excitações elementares) e a nucleação de vórtices quantizados. Embora a dependência da pressão da velocidade crítica para a emissão de rótons ($v_{rot}^{crit}$) e para a criação de vórtices ($v_{v}^{crit}$) tenha sido investigada experimentalmente — particularmente pelo grupo de Lancaster usando íons negativos —, esses fenômenos historicamente careceram de uma descrição teórica unificada. Observações experimentais indicam que, em baixas pressões, a nucleação de vórtices limita a velocidade do íon, whereas, em altas pressões (acima de $\approx 10-12$ bar), a emissão de rótons torna-se o fator limitante. Além disso, a velocidade crítica para a emissão de rótons diminui com a pressão, enquanto a velocidade crítica para a nucleação de vórtices aumenta. Estruturas teóricas existentes, como a equação de Gross-Pitaevskii padrão (GPE), falham em capturar a emissão de rótons porque possuem um espectro de excitação estritamente monótono, incapaz de reproduzir a curva de dispersão de Landau com seu mínimo característico de "róton". Por outro lado, teorias microscópicas frequentemente lutam para descrever ambos os mecanismos dentro de um único arcabouço.

Metodologia
Os autores empregam um modelo de Gross-Pitaevskii generalizado e não local (gnlGPE) para simular a dinâmica do superfluido $^4$He a temperatura zero ($T \to 0$). O modelo é definido por um parâmetro de ordem complexo $\psi(x,t)$ e incorpora um potencial de interação interatômica não local $V(r)$ e um termo de correção de densidade de ordem superior.

1. Calibração do Modelo: Para garantir precisão física em uma faixa de pressões (0 a 24 bar, aproximando-se da pressão de solidificação de $\approx 25$ bar), os autores não utilizam uma forma analítica fixa para o potencial de interação. Em vez disso, invertem a relação de dispersão de Bogoliubov generalizada usando dados experimentais de Godfrin et al. (2021). Isso permite que o modelo reproduza as curvas de dispersão exatas dependentes da pressão, incluindo os ramos de fônon, máxon e róton, bem como o platô de Pitaevskii.
2. Configuração da Simulação: O estudo é conduzido principalmente em uma aproximação bidimensional (2D). Os autores simulam um obstáculo circular (disco) movendo-se através de He II em repouso. O obstáculo é modelado como um potencial externo com perfil gaussiano que esgota completamente a densidade do superfluido.
3. Determinação da Velocidade Crítica:
   • Emissão de Rótons ($v_{rot}^{crit}$): Determinada analisando o diagrama de bifurcação de soluções estacionárias. Os autores utilizam um método de Newton-Raphson acoplado a um método de gradiente bi-conjugado estabilizado para encontrar o número de Mach máximo ($M = U/c$) para o qual soluções estacionárias estáveis existem. Além desse ponto, o sistema torna-se instável e rótons são emitidos.
   • Nucleação de Vórtices ($v_{v}^{crit}$): Determinada via uma abordagem dinâmica. O obstáculo é acelerado até atingir uma velocidade-alvo, e o sistema é evoluído por uma escala de tempo de $\approx 10\xi/c$ (onde $\xi$ é o comprimento de cura) para observar se vórtices quantizados são nucleados.
4. Parâmetros: As simulações cobrem pressões $P = 0, 2, 5, 10,$ e $24$ bar. Os tamanhos do obstáculo variam de $D = 10\xi$ a $D = 1000\xi$. Atenção específica é dada ao tamanho dos íons negativos, que variam com a pressão (de $\approx 34,5$ Å a 0 bar até $\approx 21,7$ Å a 24 bar).

Principais Resultados

• Dependência da Pressão das Velocidades Críticas: As simulações confirmam as tendências opostas observadas experimentalmente. À medida que a pressão aumenta, o número de Mach crítico para a emissão de rótons ($M_{rot}^{crit}$) diminui, caindo de $\approx 0,247$ a 0 bar para $\approx 0,128$ a 24 bar. Por outro lado, o número de Mach crítico para a nucleação de vórtices ($M_{v}^{crit}$) aumenta com a pressão.
• Dependência do Tamanho: Para obstáculos pequenos (comparáveis aos tamanhos de íons negativos, $D < 100$ nm), a velocidade crítica é determinada pelo menor dos dois limiares ($v_{rot}^{crit}$ ou $v_{v}^{crit}$). Em baixas pressões, $v_{v}^{crit}$ é o fator limitante; em altas pressões, $v_{rot}^{crit}$ é limitante. À medida que o tamanho do obstáculo aumenta ($D \ge 1000\xi$), as velocidades críticas para todas as pressões colapsam em uma única curva governada pela nucleação de vórtices, semelhante às previsões da GPE padrão.
• A "Fenda": Uma descoberta significativa é o alargamento da lacuna entre $v_{rot}^{crit}$ e $v_{v}^{crit}$ à medida que a pressão aumenta. Em altas pressões, existe uma janela de velocidade onde rótons são emitidos, mas a nucleação de vórtices ainda não ocorreu. Isso sugere que a emissão de rótons pode ser isolada e observada experimentalmente em altas pressões sem a complicação imediata da geração de vórtices.
• Visualização: Visualizações de densidade confirmam que, em números de Mach intermediários (por exemplo, $M=0,28$), rótons são emitidos em todas as pressões, enquanto a nucleação de vórtices ocorre apenas em números de Mach mais altos ($M \approx 0,3$), com o limiar para nucleação deslocando-se para valores mais altos à medida que a pressão aumenta.

Significado e Alegações
Os autores afirmam que este trabalho representa a primeira tentativa de analisar a dependência da pressão da emissão de rótons e da nucleação de vórtices, juntamente com suas velocidades críticas associadas, dentro de um único arcabouço teórico. Ao utilizar a GPE generalizada e não local, o estudo conecta com sucesso a lacuna entre o critério de Landau (emissão de rótons) e a nucleação hidrodinâmica de vórtices.

O artigo enfatiza que, embora o estudo esteja restrito a 2D devido à complexidade computacional de simulações 3D (particularmente no que diz respeito à determinação da semente mínima para nucleação de vórtices em 3D), os resultados 2D fornecem uma explicação qualitativa e semi-quantitativa da física subjacente que se alinha com observações experimentais. Os autores observam que seus valores numéricos para velocidades críticas são comparáveis aos dados experimentais, validando a abordagem de usar relações de dispersão derivadas experimentalmente para restringir o potencial de interação. O trabalho oferece uma explicação motivada fisicamente para a razão pela qual o caráter do movimento do íon muda na pressão crítica $P_{crit} \approx 10-12$ bar, atribuindo-a ao cruzamento das curvas dependentes da pressão para $v_{rot}^{crit}$ e $v_{v}^{crit}$.

Os autores permanecem modestos quanto às implicações em 3D, reconhecendo que, embora as simulações 2D expliquem as tendências gerais, a topologia específica da nucleação de vórtices em 3D (por exemplo, loops íon-vórtice versus anéis) permanece uma questão em aberto que requer investigação futura. Eles também observam que modelar o íon como uma esfera rígida perfeita é uma simplificação, embora as tendências de pressão devam permanecer robustas.