Optimal preparation and reachable-state constraints in the Mpemba effect

Resumo Técnico: Preparação Ótima e Restrições de Estados Alcançáveis no Efeito Mpemba

Definição do Problema
O efeito Mpemba, caracterizado por um sistema inicialmente mais quente que relaxa mais rapidamente do que um mais frio em direção a um estado final comum, é tipicamente analisado dentro de estruturas cinéticas assumindo condições iniciais não estacionárias selecionadas a priori. No entanto, essa abordagem frequentemente negligencia as restrições físicas impostas pelo protocolo de envelhecimento específico utilizado para preparar esses estados. Na formulação cinética, o efeito depende da memória da história térmica do sistema, codificada em variáveis internas que se desviam do equilíbrio. Embora trabalhos recentes tenham introduzido protocolos (ex: Pontus–Mpemba) que consideram o tempo de preparação, a questão fundamental permanece: quais estados não equilibrados são fisicamente acessíveis através de mecanismos de condução admissíveis, e como essas restrições de acessibilidade limitam a magnitude do efeito Mpemba? Este artigo aborda o problema de determinar o protocolo de envelhecimento ótimo para maximizar o efeito Mpemba dentro de uma estrutura cinética, especificamente para fluidos granulares uniformemente aquecidos, ao considerar explicitamente o conjunto alcançável de estados imposto pelo termostato estocástico.

Metodologia
Os autores modelam um gás granular de $d$ dimensões composto por esferas rígidas com colisões inelásticas (coeficiente de restituição $\alpha$) impulsionado por um termostato estocástico (ruído branco). O estado mesoscópico do sistema é descrito pela temperatura granular $T$ e pelo excesso de curtose $a_2$, que quantifica a não-gaussianidade da função de distribuição de velocidade (VDF). A análise é conduzida dentro da primeira aproximação de Sonine, onde a evolução de $T$ e $a_2$ é governada por um conjunto fechado de equações diferenciais acopladas.

A etapa de preparação é formulada como um problema de controle ótimo:

1. Variáveis de Estado: A temperatura granular $T$ e o excesso de curtose $a_2$.
2. Variável de Controle: A intensidade do termostato estocástico, $\chi(t)$, limitada por $[\chi_{\min}, \chi_{\max}]$.
3. Objetivo: Encontrar o protocolo de controle $\chi(t)$ que extremaniza o excesso de curtose $a_2$ em um tempo final $t_f$. Como a taxa de resfriamento a uma temperatura fixa é diretamente proporcional a $a_2$, extremar $a_2$ corresponde a maximizar ou minimizar a subsequente taxa de relaxação.
4. Técnica de Solução: Os autores aplicam o Princípio do Máximo de Pontryagin (PMP) para derivar as condições necessárias para a otimalidade. Eles analisam o sistema Hamiltoniano, as equações de coestado e as condições de transversalidade para determinar a estrutura do controle ótimo.

Principais Contribções e Resultados

• Estrutura do Protocolo Ótimo: A aplicação do PMP revela que o protocolo de preparação ótimo é sempre um protocolo "one-bang". Isso significa que o controle $\chi(t)$ permanece constante em seu valor máximo ($\chi_{\max}$) ou mínimo ($\chi_{\min}$) durante todo o intervalo de preparação $(0, t_f)$. A escolha específica depende do sinal do excesso de curtose no estado de resfriamento homogêneo (HCS), que é determinado pelo coeficiente de restituição $\alpha$.

  • Para inelasticidade "grande" ($\alpha < \alpha_c \equiv 1/\sqrt{2}$), o controle ótimo para minimizar $a_2$ é $\chi_{\max}$.
  • Para inelasticidade "pequena" ($\alpha > \alpha_c$), o controle ótimo para minimizar $a_2$ é $\chi_{\min}$.
  • A lógica é invertida para maximizar $a_2$.

• Restrições de Estados Alcançáveis: Uma descoberta central é a existência de restrições de estados alcançáveis não triviais. O conjunto de estados não equilibrados acessíveis no plano $(T, a_2)$ é limitado. Especificamente, os valores acessíveis do excesso de curtose estão confinados entre o estado Gaussiano ($a_2 = 0$) e o valor do excesso de curtose no estado de resfriamento homogêneo ($a_2^{\text{HCS}}$).

  • No limite ideal onde $\chi_{\min} = 0$ (resfriamento livre) e $\chi_{\max} \to \infty$ (aquecimento forte instantâneo), os extremos do excesso de curtose são exatamente $0$e$a_2^{\text{HCS}}$.
  • Isso implica que o termostato estocástico não pode gerar não-gaussianidades arbitrariamente grandes; a magnitude da não-gaussianidade é estritamente limitada pelo valor do HCS.

• Validação Numérica: As previsões teóricas foram verificadas utilizando simulações de Monte Carlo de Direção Direta (DSMC). As simulações confirmaram que os extremos do excesso de curtose convergem para os limites teóricos ($0$e$a_2^{\text{HCS}}$) à medida que os limites do controle são relaxados. A concordância entre a teoria do controle ótimo e os dados numéricos é excelente em toda a gama de coeficientes de restituição.

• Implicações para o Efeito Mpemba: O estudo demonstra que a força do efeito Mpemba cinético está intrinsecamente ligada à acessibilidade dos estados não equilibrados subjacentes. Como o conjunto alcançável de $a_2$ é limitado, existe um limite fundamental para a assimetria da taxa de relaxação que pode ser alcançada entre amostras quentes e frias. As "grandes" não-gaussianidades anteriormente hipotetizadas em trabalhos seminais não podem ser geradas apenas por um termostato estocástico.

Significância e Alegações
O artigo afirma estabelecer uma conexão rigorosa entre relaxação anômala, controle ótimo e acessibilidade de estado em sistemas fora do equilíbrio. Ao recategorizar o efeito Mpemba como um problema de preparação de estado ótima sob restrições de alcançabilidade dinâmica, os autores argumentam que o efeito não pode ser discutido independentemente das limitações físicas do protocolo de preparação.

O trabalho destaca que o excesso de curtose atua como uma variável interna crítica que codifica a memória do sistema, de forma semelhante ao seu papel no efeito Kovacs. Os resultados sugerem que o efeito Mpemba máximo atingível não é uma propriedade intrínseca apenas do material, mas é fundamentalmente limitado pela geometria do espaço de estado não equilibrado alcançável definido pelo mecanismo de condução. Os autores observam que, embora o termostato estocástico imponha esses limites específicos, mecanismos de condução alternativos podem produzir conjuntos alcançáveis diferentes, permitindo potencialmente anomalias de relaxação mais fortes, embora isso permaneça um assunto para investigações futuras. O artigo conclui que a engenharia de caminhos de relaxação requer navegar pela geometria específica do espaço de estado acessível.