Probing the onset of hydrodynamization in peripheral p-Pb collisions at $\sqrt{s_{NN}} =$ 5.02 TeV

O essencial

Imagine que você está tentando descobrir a menor poça de água possível que ainda possa agir como um fluido. Se você tem um oceano gigante, ele flui facilmente. Se você tem uma única gota, ela pode apenas ficar parada ou se desintegrar. Mas onde está a linha? Em que tamanho um conjunto de moléculas de água deixa de agir como um fluido e começa a agir como partículas individuais e caóticas?

Este artigo trata de encontrar esse exato "ponto de virada" para o Plasma de Quarks e Glúons (QGP).

O que é o QGP?

Pense no QGP como a "sopa primordial" do universo. É um estado da matéria que existiu apenas frações de segundo após o Big Bang. Nesse estado, os blocos de construção dos átomos (quarks e glúons) estão derretidos juntos e fluem livremente, como um líquido superquente e superdenso.

Normalmente, os cientistas criam essa sopa colidindo dois átomos pesados (como chumbo) a velocidades próximas à da luz. Mas, recentemente, os cientistas notaram algo intrigante: mesmo quando colidem coisas muito menores — como um único próton atingindo um núcleo de chumbo (colisões p-Pb) —, sinais dessa "sopa líquida" aparecem.

A grande questão é: É realmente um líquido, ou é apenas um grupo de partículas quicando de forma caótica?

O Experimento: Colidir Prótons com Chumbo

Os autores deste artigo queriam encontrar o menor tamanho dessa "sopa" que ainda possa ser descrito pelas leis da hidrodinâmica (a matemática usada para descrever líquidos em fluxo).

Eles usaram uma simulação computacional massiva chamada JETSCAPE. Pense nessa simulação como um motor de videogame de alta tecnologia que recria todo o processo de colisão em quatro etapas:

1. A Configuração (TRENTo): Eles montaram o cenário, posicionando os prótons e núcleos de chumbo em suas posições iniciais.
2. O Pré-Jogo (Freestreaming): Antes que o "líquido" se forme, as partículas voam livremente por uma fração minúscula de segundo.
3. O Fluxo (MUSIC): Esta é a parte da hidrodinâmica. A simulação tenta tratar as partículas como um fluido em fluxo.
4. O Após-Jogo (iSS + SMASH): À medida que a sopa esfria, as partículas congelam em prótons reais, píons e outras partículas que os detectores podem observar.

O Teste: Quão "Líquida" é a Sopa?

Para testar se a sopa está realmente se comportando como um fluido, os cientistas observaram algo chamado Fluxo Elíptico.

A Analogia: Imagine dois carros colidindo de frente. Se eles forem perfeitamente redondos e atingirem exatamente o centro, os detritos voam para fora em um círculo. Mas se eles atingirem ligeiramente fora do centro (um golpe de raspão), os detritos voam para fora mais em formato oval (como um futebol americano).

• Se a matéria dentro agir como um fluido perfeito, ela será espremida fortemente nesse formato oval.
• Se a matéria for apenas partículas caóticas quicando, o formato oval será fraco ou inexistente.

Os cientistas executaram sua simulação para colisões "periféricas" (golpes de raspão onde a sobreposição entre o próton e o núcleo de chumbo é pequena). Eles perguntaram: Quão pequena essa sobreposição pode ficar antes que o comportamento fluido se desintegre?

A Reviravolta: O Botão "Tempo de Relaxação"

Em fluidos reais, há um atraso entre o momento em que você empurra o fluido e o momento em que ele responde. Na física, isso é chamado de tempo de relaxação de cisalhamento.

Os autores fizeram uma manobra: eles giraram esse botão de "tempo de relaxação" para configurações extremas.

• Eles perguntaram: "E se o fluido for muito lento para responder? E se for muito rápido?"
• Eles observaram o Fluxo Elíptico (o formato oval) sob essas condições extremas.

A Descoberta: O Ponto de Virada

À medida que simulavam colisões cada vez mais "de raspão" (o que significa que a quantidade de matéria envolvida, ou dN/dy, ficava menor), eles observaram o comportamento fluido.

• O Resultado: Quando a quantidade de matéria caiu para cerca de 7 partículas por unidade de rapidez (dN/dy ≈ 7), o comportamento fluido começou a oscilar e se desintegrar repentinamente.
• A Metáfora: Imagine uma multidão de pessoas tentando se mover como um fluido. Se você tem 100 pessoas, elas fluem suavemente. Se você tem 10, elas ainda podem fluir. Mas se chegar a 7 pessoas, elas começam a bater umas nas outras individualmente, e o "fluxo" suave desaparece.

O artigo conclui que, para colisões próton-chumbo na energia que estudaram, a hidrodinâmica para de funcionar quando o sistema fica menor do que cerca de 7 partículas. Abaixo disso, a "sopa" é pequena demais para agir como um líquido; é apenas um grupo de partículas individuais.

Por Que Isso Importa?

Isso ajuda os cientistas a entender os limites fundamentais da natureza. Isso nos diz que o estado da matéria "líquido" não é mágico; ele tem um requisito mínimo de tamanho. Se o sistema for muito pequeno, as regras da dinâmica dos fluidos não se aplicam mais, e precisamos olhar para as partículas individuais.

Os autores também observaram que seus resultados foram ligeiramente diferentes de seus estudos anteriores sobre colisões maiores (como chumbo-chumbo), provavelmente porque os modelos computacionais usados desta vez eram mais estáveis e lidavam de forma diferente com a fase de "pré-jogo".

Em resumo: Eles encontraram a menor poça de plasma de quarks e glúons que ainda pode ser chamada de "fluido", e acontece que essa poça precisa conter pelo menos cerca de 7 partículas para se manter unida.

Resumo técnico

Resumo Técnico: Sondando o Início da Hidrodinamização em Colisões p-Pb Periféricas

Declaração do Problema
Embora a hidrodinâmica relativística tenha modelado com sucesso o Plasma de Quarks e Glúons (PQG) em colisões de íons pesados grandes (Pb-Pb, Au-Au), sua aplicabilidade a sistemas de colisão pequenos (próton-próton de alta multiplicidade e próton-núcleo) permanece um tema de intenso debate. A questão central é se as assinaturas de fluxo coletivo observadas nesses sistemas pequenos surgem de efeitos reais de matéria nuclear quente semelhantes ao PQG ou de efeitos de matéria nuclear fria. Especificamente, os autores visam determinar o menor tamanho de sistema (quantificado pela multiplicidade de partículas, $dN/dy$) no qual as descrições hidrodinâmicas se rompem. Trabalhos anteriores dos autores estabeleceram um início de hidrodinamização em $dN/dy \approx 10$ para colisões periféricas de Pb-Pb e Au-Au; este estudo estende essa investigação para o sistema p-Pb menor em $\sqrt{s_{NN}} = 5.02$ TeV.

Metodologia
O estudo utiliza o framework JETSCAPE (versão 3.6.4), uma ferramenta de simulação modular projetada para modelar consistentemente os setores macio e duro de colisões de alta energia. A cadeia de simulação para o setor macio inclui:

1. Condições Iniciais: O modelo TRENTo gera a densidade de energia inicial com base em uma média generalizada das funções de espessura nuclear, incorporando flutuações de núcleons.
2. Pré-equilíbrio: Um módulo de Livre Caminhamento (Freestreaming) baseado na equação de Boltzmann sem colisões preenche a lacuna entre o estado inicial e a fase hidrodinâmica.
3. Hidrodinâmica: O solucionador MUSIC implementa hidrodinâmica viscosa de segunda ordem (formalismo BRSSS). O tempo de relaxação do cisalhamento ($\tau_\pi$), um coeficiente de transporte de segunda ordem que regula modos não hidrodinâmicos, é a variável principal de interesse.
4. Pós-processador de Hádrons: Os módulos iSS (particulação) e SMASH (transporte hadrônico) lidam com o congelamento e o subsequente espalhamento hadrônico.

Para sondar o início da hidrodinamização, os autores variam o tempo de relaxação do cisalhamento ($\tau_\pi$) para valores extremos. Como $\tau_\pi$ está relacionado à razão entre viscosidade de cisalhamento e densidade de entropia ($\eta/s$) por $\tau_\pi(T) = b_\pi \frac{\eta}{s(T)} \frac{1}{T}$, eles empregam duas estratégias distintas de variação:

1. Variando a constante de normalização $b_\pi$ (de 1 a 20) enquanto mantêm o $\eta/s$ dependente da temperatura fixo.
2. Variando os limites de $\eta/s$ dependentes da temperatura (usando o limite inferior KSS e os limites superiores bayesianos) enquanto mantêm $b_\pi$ fixo.

A ruptura do comportamento hidrodinâmico é identificada monitorando o fluxo elíptico ($v_2$) e suas flutuações. Argumentos teóricos sugerem que, se a contribuição dos modos não hidrodinâmicos (regulados por $\tau_\pi$) exceder a dos modos hidrodinâmicos, a descrição de fluido falha, levando a desvios abruptos nas observáveis de fluxo.

Contribuições e Resultados Principais

• Benchmarking: O modelo reproduz com sucesso dados experimentais do ALICE, ATLAS e CMS para espectros de rapidez de partículas carregadas, espectros de momento transversal ($p_T$) e momento transversal médio em função da centralidade (definida por $N_{track}$).
• Análise de Flutuações: Ao simular colisões p-Pb periféricas e variar $\tau_\pi$, os autores observam que as flutuações do fluxo elíptico aumentam à medida que o sistema se move de colisões centrais para periféricas.
• Determinação do Início:
  • Ao variar $b_\pi$, observa-se um aumento abrupto nas flutuações do fluxo elíptico na faixa de centralidade de 75–85% a 70–80%.
  • Ao variar $\eta/s$, as flutuações começam a aumentar abaixo da faixa de centralidade de 70–80%.
  • Essas flutuações correspondem a uma multiplicidade de partículas carregadas de $dN/dy \approx 7$.
• Comparação com Trabalhos Anteriores: Diferentemente do estudo anterior dos autores sobre sistemas maiores (onde o início foi encontrado em $dN/dy \approx 10$), a transição em p-Pb aparece em uma multiplicidade menor. Os autores observam que o desvio abaixo do início inferido é menos abrupto do que em estudos anteriores, possivelmente devido ao uso do solucionador hidrodinâmico MUSIC mais estável em comparação com o ECHO-QGP anteriormente utilizado.

Significado e Afirmações
O artigo afirma fornecer uma estimativa para o menor volume de PQG em colisões p-Pb a 5.02 TeV que pode ser satisfatoriamente modelado por hidrodinâmica de baixa ordem. A conclusão principal é que o comportamento hidrodinâmico se rompe em aproximadamente $dN/dy \approx 7$.

Os autores mantêm um tom modesto quanto à natureza dessa transição, reconhecendo que permanece uma questão em aberto se isso representa um início abrupto ou uma região de transição suave. Eles também destacam que seu estudo isola o setor macio; uma confirmação definitiva do início idealmente exigiria uma investigação simultânea do setor duro (por exemplo, supressão de jatos) para ver como observáveis duros e macios influenciam mutuamente a determinação da hidrodinamização. O trabalho serve como uma aplicação específica da análise de modos não hidrodinâmicos para restringir os limites de validade de modelos hidrodinâmicos em sistemas de colisão pequenos.