Optimal preparation and reachable-state constraints in the Mpemba effect

Resumen Técnico: Preparación Óptima y Restricciones de Estados Alcanzables en el Efecto Mpemba

Planteamiento del Problema
El efecto Mpemba, caracterizado por un sistema inicialmente más caliente que se relaja más rápido que uno más frío hacia un estado final común, se analiza típicamente dentro de marcos cinéticos asumiendo condiciones iniciales no estacionarias seleccionadas a priori. Sin embargo, este enfoque a menudo descuida las restricciones físicas impuestas por el protocolo de envejecimiento específico utilizado para preparar dichos estados. En la formulación cinética, el efecto depende de la memoria de la historia térmica del sistema, codificada en variables internas que se desvían del equilibrio. Si bien trabajos recientes han introducido protocolos (p. ej., Pontus–Mpemba) que consideran el tiempo de preparación, la pregunta fundamental sigue siendo: ¿qué estados fuera del equilibrio son físicamente accesibles mediante mecanismos de excitación admisibles, y cómo estas restricciones de accesibilidad limitan la magnitud del efecto Mpemba? Este artículo aborda el problema de determinar el protocolo de envejecimiento óptimo para maximizar el efecto Mpemba dentro de un marco cinético, específicamente para fluidos granulares calentados uniformemente, considerando explícitamente el conjunto alcanzable de estados impuesto por el termostato estocástico.

Metodología
Los autores modelan un gas granular de $d$ dimensiones de esferas duras con colisiones inelásticas (coeficiente de restitución $\alpha$) impulsado por un termostato estocástico (ruido blanco). El estado mesoscópico del sistema se describe mediante la temperatura granular $T$ y la curtosis excedente $a_2$, que cuantifica la no-gaussianidad de la función de distribución de velocidades (VDF). El análisis se lleva a cabo dentro de la primera aproximación de Sonine, donde la evolución de $T$ y $a_2$ está gobernada por un conjunto cerrado de ecuaciones diferenciales acopladas.

La etapa de preparación se formula como un problema de control óptimo:

1. Variables de Estado: La temperatura granular $T$ y la curtosis excedente $a_2$.
2. Variable de Control: La intensidad del termostato estocástico, $\chi(t)$, acotada por $[\chi_{\min}, \chi_{\max}]$.
3. Objetivo: Encontrar el protocolo de control $\chi(t)$ que extremice la curtosis excedente $a_2$ en un tiempo final $t_f$. Dado que la tasa de enfriamiento a una temperatura fija es directamente proporcional a $a_2$, extremizar $a_2$ corresponde a maximizar o minimizar la subsiguiente tasa de relajación.
4. Técnica de Solución: Los autores aplican el Principio de Máximo de Pontryagin (PMP) para derivar las condiciones necesarias de optimalidad. Analizan el sistema hamiltoniano, las ecuaciones de los coestados y las condiciones de transversalidad para determinar la estructura del control óptimo.

Contribuciones Clave y Resultados

• Estructura del Protocolo Óptimo: La aplicación del PMP revela que el protocolo de preparación óptimo es siempre un protocolo de "un solo golpe" (one-bang). Esto significa que el control $\chi(t)$ permanece constante en su valor máximo ($\chi_{\max}$) o mínimo ($\chi_{\min}$) durante todo el intervalo de preparación $(0, t_f)$. La elección específica depende del signo de la curtosis excedente en el estado de enfriamiento homogéneo (HCS), el cual está determinado por el coeficiente de restitución $\alpha$.

  • Para una inelasticidad "grande" ($\alpha < \alpha_c \equiv 1/\sqrt{2}$), el control óptimo para minimizar $a_2$ es $\chi_{\max}$.
  • Para una inelasticidad "pequeña" ($\alpha > \alpha_c$), el control óptimo para minimizar $a_2$ es $\chi_{\min}$.
  • La lógica se invierte para maximizar $a_2$.

• Restricciones de Estados Alcanzables: Un hallazgo central es la existencia de restricciones de estados alcanzables no triviales. El conjunto de estados fuera del equilibrio accesibles en el plano $(T, a_2)$ está acotado. Específicamente, los valores accesibles de la curtosis excedente están confinados entre el estado gaussiano ($a_2 = 0$) y el valor de la curtosis excedente en el estado de enfriamiento homogéneo ($a_2^{\text{HCS}}$).

  • En el límite ideal donde $\chi_{\min} = 0$ (enfriamiento libre) y $\chi_{\max} \to \infty$ (calentamiento fuerte instantáneo), los extremos de la curtosis excedente son exactamente $0$y$a_2^{\text{HCS}}$.
  • Esto implica que el termostato estocástico no puede generar no-gaussianidades arbitrariamente grandes; la magnitud de la no-gaussianidad está estrictamente limitada por el valor de la HCS.

• Validación Numérica: Las predicciones teóricas fueron verificadas mediante simulaciones de Monte Carlo de Simulación Directa (DSMC). Las simulaciones confirmaron que los extremos de la curtosis excedente convergen a los límites teóricos ($0$y$a_2^{\text{HCS}}$) a medida que se relajan los límites del control. El acuerdo entre la teoría de control óptimo y los datos numéricos es excelente en todo el rango de coeficientes de restitución.

• Implicaciones para el Efecto Mpemba: El estudio demuestra que la fuerza del efecto Mpemba cinético está intrínsecamente ligada a la accesibilidad de los estados fuera del equilibrio subyacentes. Debido a que el conjunto alcanzable de $a_2$ está acotado, existe un límite fundamental a la asimetría de la tasa de relajación que se puede lograr entre muestras calientes y frías. Las grandes no-gaussianidades hipotetizadas previamente en trabajos seminales no pueden ser generadas únicamente por un termostato estocástico.

Significancia y Reivindicaciones
El artículo afirma establecer una conexión rigurosa entre la relajación anómala, el control óptimo y la accesibilidad de estados en sistemas fuera del equilibrio. Al reformular el efecto Mpemba como un problema de preparación de estado óptimo bajo restricciones de alcanzabilidad dinámica, los autores argumentan que el efecto no puede discutirse independientemente de las limitaciones físicas del protocolo de preparación.

El trabajo destaca que la curtosis excedente actúa como una variable interna crítica que codifica la memoria del sistema, de manera similar a su papel en el efecto Kovacs. Los resultados sugieren que el efecto Mpemba máximo alcanzable no es una propiedad intrínseca del material por sí solo, sino que está fundamentalmente limitado por la geometría del espacio de estados fuera del equilibrio alcanzable definido por el mecanismo de excitación. Los autores señalan que, si bien el termostato estocástico impone estos límites específicos, otros mecanismos de excitación podrían producir conjuntos alcanzables diferentes, permitiendo potencialmente efectos de relajación más fuertes, aunque esto sigue siendo objeto de investigación futura. El artículo concluye que la ingeniería de vías de relajación requiere navegar la geometría específica del espacio de estados accesible.