🔬 condensed matter

Crystal to liquid cross-over for active particles with inverse-square power-law interaction

Este artículo investiga un sistema unidimensional de partículas de carrera y vuelco (run-and-tumble) con repulsión de inverso al cuadrado en una trampa armónica, revelando un distintivo cruce impulsado por la actividad desde un estado de tipo cristal con picos de densidad agudos hacia un perfil de semicírculo de Wigner de tipo líquido, respaldado tanto por cálculos analíticos de covarianza como por simulaciones numéricas.

Autores originales: Saikat Santra, Leo Touzo, Chandan Dasgupta, Abhishek Dhar, Suman Dutta, Anupam Kundu, Pierre Le Doussal, Gregory Schehr, Prashant Singh

Publicado 2026-08-18

📖 6 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Saikat Santra, Leo Touzo, Chandan Dasgupta, Abhishek Dhar, Suman Dutta, Anupam Kundu, Pierre Le Doussal, Gregory Schehr, Prashant Singh

Artículo original bajo licencia CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Esta es una explicación generada por IA del artículo a continuación. No ha sido escrita ni avalada por los autores. Para mayor precisión técnica, consulte el artículo original. Leer descargo de responsabilidad completo

Imagina un mundo donde los objetos diminutos no son empujados por temblores térmicos invisibles, sino por sus propios motores internos. Estas son partículas activas, una clase de materia que se encuentra en todo, desde bacterias que nadan hasta robots autopropulsados. A diferencia de las motas de polvo ordinarias que derivan sin rumbo, estas partículas consumen energía para moverse en una dirección específica durante un tiempo antes de cambiar de rumbo aleatoriamente. Este movimiento constante y autoimpulsado rompe las reglas habituales del equilibrio, creando comportamientos extraños y complejos que los científicos apenas están comenzando a comprender. Cuando muchas de estas partículas activas se agrupan y se confinan, no se mezclan simplemente en una sopa uniforme. En su lugar, pueden organizarse en patrones distintos, alternando entre estructuras similares a sólidos y flujos similares a líquidos de formas que la materia pasiva nunca hace.

Un equipo de investigadores ha mapeado ahora exactamente cómo ocurre esta transformación en un entorno específico y controlado. Estudiaron una línea unidimensional de partículas atrapadas dentro de un campo de fuerza en forma de cuenco, donde cada partícula se repele de sus vecinas con una fuerza que se vuelve más fuerte cuanto más cerca están. Crucialmente, estas partículas fueron modeladas como agentes de "correr y dar tumbos" (run-and-tumble), lo que significa que corren hacia adelante a una velocidad constante y luego dan un vuelco aleatorio para elegir una nueva dirección. Mediante la ejecución de extensas simulaciones por computadora y el desarrollo de nuevas herramientas matemáticas, el equipo descubrió que, a medida que aumenta la velocidad de estas partículas, todo el sistema experimenta una evolución dramática y predecible. Comienza como un cristal rígido, se funde en un líquido suave y, finalmente, se estira hasta convertirse en una amplia nube con forma de campana.

Al principio, cuando las partículas se mueven muy lentamente, se comportan como un cristal congelado. Debido a que la fuerza de repulsión entre ellas es tan fuerte, se bloquean en posiciones precisas y fijas entre sí. Si tomaras una instantánea de dónde se encuentran a lo largo del tiempo, verías picos nítidos y distintos en la densidad, donde cada pico representa una partícula manteniendo su posición. Los investigadores descubrieron que estas posiciones corresponden a los ceros matemáticos de un tipo específico de polinomio, un patrón que se ha conocido durante décadas en la física pasiva pero que aquí emerge del movimiento activo. Mientras la energía interna de las partículas sea baja, estas simplemente se tambalean ligeramente alrededor de estos puntos fijos, manteniendo su estructura ordenada y sólida.

A medida que los investigadores aumentaban la velocidad de las partículas, el sistema comenzaba a fundirse. Las partículas empezaban a alejarse más de sus lugares originales y sus movimientos comenzaban a solaparse. Los picos nítidos en el perfil de densidad se suavizaron, fundiéndose hasta que las posiciones distintas de las partículas individuales ya no podían distinguirse. En esta etapa intermedia, el sistema parecía un líquido. Sorprendentemente, la forma general de esta nube líquida coincidía con una famosa curva matemática conocida como el semicírculo de Wigner. Esta forma está típicamente asociada con sistemas pasivos en equilibrio térmico, pero aquí aparecía en un sistema impulsado lejos del equilibrio por el movimiento activo. Los investigadores confirmaron esta transición midiendo cuánto fluctuaban las partículas en comparación con la distancia entre ellas. Cuando estas fluctuaciones se volvieron lo suficientemente grandes como para desdibujar los huecos entre los vecinos, el cristal se había convertido efectivamente en un líquido.

Sin embargo, la historia no terminó ahí. Cuando las partículas fueron impulsadas a moverse aún más rápido, el semicírculo suave comenzó a distorsionarse. El perfil de densidad se estiró, perdiendo su parte superior curva y adoptando una forma de campana ancha que se extendía mucho más que la fase líquida. En este estado altamente activo, las partículas ya no solo se agitaban en su lugar; estaban explorando la trampa con tal vigor que los bordes de la nube quedaron definidos por su velocidad máxima, en lugar del equilibrio de fuerzas. Los investigadores observaron que la densidad en los bordes mismos de esta nube no caía abruptamente, sino que se desvanecía siguiendo una ley de potencia específica, un patrón que parece ser una característica universal de tales sistemas cuando las partículas experimentan una repulsión infinita al contacto.

Para entender por qué ocurrieron estos cambios, el equipo desarrolló un marco teórico basado en la idea de que las partículas estaban realizando pequeñas desviaciones de sus posiciones ideales y congeladas. Al calcular cómo estas desviaciones se correlacionaban entre sí, pudieron predecir la varianza, o la dispersión, de las posiciones de las partículas. Sus cálculos mostraron que la transición del cristal al líquido ocurre cuando la velocidad de la partícula alcanza un cierto umbral, independientemente de cuántas partículas haya en el sistema. La segunda transición, del líquido a la nube con forma de campana, ocurre solo cuando la velocidad es tan alta que escala con el número total de partículas. Este enfoque analítico coincidió perfectamente con sus simulaciones por computadora, confirmando que el cambio del orden al desorden, y luego a un nuevo tipo de desorden activo, está gobernado por reglas simples y predecibles.

Los hallazgos ofrecen una imagen clara de cómo la materia activa se organiza bajo confinamiento. Los investigadores demostraron que, incluso en un sistema impulsado lejos del equilibrio, las partículas pueden establecerse en un estado que imita el comportamiento suave y de tipo líquido de los sistemas pasivos, siempre que la actividad no sea demasiado extrema. Pero una vez que la actividad cruza un límite crítico, el sistema se libera de estas formas familiares, adoptando una nueva forma dictada por el puro impulso de las partículas. Este trabajo no solo clarifica el comportamiento de las partículas activas en una trampa, sino que también sugiere que transiciones similares podrían ser observables en experimentos del mundo real con bacterias o microrobots sintéticos, siempre que puedan ser confinados y sus interacciones controladas. El estudio es un testimonio del poder de combinar la simulación precisa con una teoría rigurosa para revelar el orden oculto dentro de los sistemas complejos y autoimpulsados.

Resumen Técnico: Transición de Cristal a Líquido en Partículas Activas con Interacción de Ley de Potencia Inversa al Cuadrado

Planteamiento del Problema
Este estudio investiga las propiedades del estado estacionario fuera del equilibrio de un sistema unidimensional compuesto por $N$ partículas de carrera y giro (RTPs, por sus siglas en inglés) confinadas en una trampa armónica. Las partículas interactúan mediante un potencial repulsivo de ley de potencia inversa al cuadrado ( $V(x) \propto 1/x^2$ ), un modelo conocido en la dinámica hamiltoniana como el modelo de Calogero-Moser o el Modelo de Inversa al Cuadrado (ISM). El objetivo principal es comprender cómo la introducción de la actividad (autopropulsión) altera el perfil de densidad y el orden estructural del sistema, examinando específicamente la transición de un estado cristalino a un estado de tipo líquido y, eventualmente, a un régimen fuertemente activo distinto. Este trabajo extiende estudios analíticos previos sobre el movimiento de Dyson-Brown (DBM) activo, que presenta interacciones logarítmicas, al ISM activo para explorar las conexiones entre estos dos sistemas integrables en presencia de un accionamiento fuera del equilibrio.

Metodología
Los autores emplean un enfoque dual que combina simulaciones numéricas extensas y derivaciones analíticas en el límite de ruido débil.

Simulaciones Numéricas: El sistema se evoluciona de acuerdo con la ecuación de movimiento de tipo Langevin para las RTPs con repulsión inversa al cuadrado y un potencial de confinamiento armónico. El perfil de densidad en estado estacionario $\rho(x)$ se calcula promediando las posiciones de las partículas a lo largo del tiempo y de múltiples realizaciones. Los observables clave incluyen:
- El perfil de densidad global y su escalamiento con el tamaño del sistema $N$ .
- El cociente de Lindemann ( $\eta_i$ ), definido como la relación entre la desviación estándar de las posiciones de las partículas y la separación media entre partículas, utilizado para cuantificar la transición cristal-líquido.
- Una medida de distancia ( $\chi$ ), que cuantifica la desviación del perfil de densidad en estado estacionario respecto al perfil de semicírculo de Wigner (WSc), utilizada para identificar la transición a un régimen fuertemente activo.
- La varianza de las brechas entre partículas.
Enfoque Analítico: En el límite de actividad pequeña ( $v_0 \to 0$ ) y tiempo de persistencia grande ( $\gamma \to 0$ ), los autores linealizan las ecuaciones de movimiento alrededor de la configuración de equilibrio (los ceros de los polinomios de Hermite). Utilizan la matriz Hessiana del potencial ISM, denotada como $H_{ISM}$ . Una idea teórica clave utilizada es la relación $H_{ISM} = H_L^2$ , donde $H_L$ es la Hessiana del DBM. Mediante la diagonalización de $H_L$ usando polinomios de Hermite y sus vectores propios asociados, los autores derivan expresiones cerradas para la covarianza y la varianza de los desplazamientos de las partículas y de las brechas entre partículas. Estos resultados se extienden posteriormente a tasas de giro finitas $\gamma$ utilizando técnicas de espacio de Fourier.

Resultos Clave

Tres Regímenes de Actividad Distintos: El estudio identifica tres regímenes distintos basados en el parámetro de actividad $v_0$ :
- Régimen Débilmente Activo ( $v_0 \ll 1$ ): El sistema se asemeja a un cristal. El perfil de densidad exhibe picos agudos centrados en las posiciones de equilibrio (ceros de Hermite). Las fluctuaciones son pequeñas en comparación con el espaciamiento entre partículas.
- Régimen de Actividad Intermedia ( $v_0 \sim O(1)$ ): A medida que la actividad aumenta, las fluctuaciones de las partículas se vuelven comparables a la distancia entre partículas. Los picos agudos desaparecen y el perfil de densidad se suaviza hacia una forma de semicírculo de Wigner (WSc) en el bulbo. Esta fase "líquida" persiste en un amplio rango de actividad.
- Régimen Fuertemente Activo ( $v_0 \sim O(N)$ ): A una actividad muy alta, el perfil de densidad se desvía significamente del WSc, transformándose en una distribución con forma de campana con un soporte que se extiende hasta $\sim \pm v_0$ . Notablemente, los bordes de esta distribución exhiben una caída de ley de potencia universal $\sim 1/x^3$ .
Criterios de Crossover Cuantitativos:
- La transición de la fase cristalina (con picos) a la fase líquida (suave WSc) ocurre cuando el cociente de Lindemann $\eta_i$ supera un umbral (aprox. 0.5). Los resultados analíticos sugieren que este crossover ocurre cuando $v_0 \sim O(1)$ , independientemente del tamaño del sistema $N$ .
- La transición de la fase líquida (WSc) a la fase fuertemente activa (forma de campana) está marcada por la medida de distancia $\chi$ superando 0.1. Esto ocurre cuando $v_0 \sim O(N)$ .
Predicciones Analíticas para las Fluctuaciones:
- Los autores derivan expresiones exactas para la varianza de las posiciones de las partículas $\langle \delta x_i^2 \rangle$ y la varianza de las brechas entre partículas $g_{i,n}^2$ en el límite de $N$ grande.
- A diferencia del DBM activo, donde las partículas de los bordes exhiben comportamientos de escalamiento diferentes a las partículas del bulbo, el ISM activo muestra que la varianza decae como $1/N$ tanto para las partículas del bulbo como para las de los bordes. Esto se atribuye a la convergencia más rápida de la expansión en serie que involucra el factor $1/k^4$ (derivado de la relación $H_{ISM} = H_L^2$ ) en comparación con el factor $1/k^2$ en el DBM.
- Las predicciones analíticas para la varianza de las brechas cerca del centro de la trampa concuerdan bien con las simulaciones numéricas incluso para $\gamma$ finito y hasta $v_0 \sim O(N)$ .
Comportamiento de los Bordes y Alas:
- En el régimen intermedio, el ISM activo no exhibe las "alas" (colas extendidas) observadas en el DBM activo. Los autores atribuyen esto al hecho de que en el ISM, la fuerza de interacción sobre las partículas de los bordes está dominada por los vecinos más cercanos (escalando como $\sqrt{N}$ ), lo que efectivamente las ancla, mientras que en el DBM, la fuerza está dominada por el bulbo colectivo (escalando como $\sqrt{N}$ pero con coeficientes diferentes), permitiendo mayores fluctuaciones.
- Sin embargo, en el régimen fuertemente activo con un $\gamma$ pequeño, pueden emerger estructuras de tipo "ala" cerca de los bordes debido al dominio de la velocidad de conducción sobre la trampa armónica.

Significancia y Reivindicaciones
El artículo afirma proporcionar una caracterización exhaustiva del ISM activo, estableciendo un paralelo directo con el DBM activo, al tiempo que destaca las diferencias cruciales en el escalamiento de las fluctuaciones y el comportamiento de los bordes.

Universalidad del Semicírculo de Wigner: El estudio refuerza el hallazgo de que el semicírculo de Wigner es un perfil de densidad robusto para partículas activas que interactúan en una trampa armónica, sobreviviendo incluso a altas actividades antes de una segunda transición hacia un perfil de forma de campana.
Tractabilidad Analítica: Al aprovechar la conexión entre las Hesianas del ISM y el DBM, los autores logran derivar fórmulas analíticas para las varianzas de posición y de brechas que siguen siendo precisas a través de un amplio rango de parámetros, incluyendo la fase "líquida".
Diagrama de Fases: La construcción de un diagrama de fases en el plano $(a, 1/\gamma)$ (donde $v_0 = aN$ ) delimita cuantitativamente los límites entre las fases cristalina, líquida y fuertemente activa.
Relevancia Experimental: Los autores sugieren que la transición de fusión cristal-líquido predicha y los perfiles de densidad específicos podrían ser observables experimentalmente en sistemas actuales de partículas atrapadas (por ejemplo, bacterias, partículas de Janus), si se pueden diseñar interacciones de ley de potencia efectivas, aunque señalan que los efectos de interacción aún no han sido plenamente explorados en tales experimentos.

El trabajo concluye que, si bien el ISM activo y el DBM activo comparten características cualitativas (tres regímenes, fase intermedia WSc), la interacción específica de ley de potencia en el ISM conduce a comportamientos cuantitativos distintos en las fluctuaciones y la estabilidad de los bordes, ofreciendo una comprensión más profunda de la mecánica estadística fuera del equilibrio en sistemas integrables.

Más como este