statistical mechanics for the phase separation of
play

Statistical mechanics for the phase separation of interacting self - PowerPoint PPT Presentation

Statistical mechanics for the phase separation of interacting self propelled particles J. Barr e, R. Ch etrite, M. Muratori, F. Peruani Laboratoire J.A. Dieudonn e, U. de Nice-Sophia Antipolis. Florence, 05-2014 Discovering the joys of


  1. Statistical mechanics for the phase separation of interacting self propelled particles J. Barr´ e, R. Ch´ etrite, M. Muratori, F. Peruani Laboratoire J.A. Dieudonn´ e, U. de Nice-Sophia Antipolis. Florence, 05-2014

  2. Discovering the joys of research with Stefano 1998 World Cup, Italy vs France. L. Di Biagio misses his penalty kick.

  3. Self propelled particles ◮ Particles with an internal source of free energy that they can convert into systematic movement. ◮ Used to model flocks of animals (from mammals to insects), bacteria, some artificial systems . . . This will be a theoretical talk! ◮ Main question: understand their collective properties ◮ Blooming field; many recent developments (I will not be able to cite all relevant contributions!).

  4. The model system (2D) particle i i ◮ Point particles with an internal angular variable θ i ◮ Move with speed v i along direction θ i + spatial noise ◮ the speed v i may depend on the local density ◮ Particles interact: they tend to align locally

  5. Microscopic equations (2D) • Spatial variables: transport in direction θ i (speed may depend on local density) + noise • Angular variables: interactions promoting local alignment + noise � x i ˙ = v ( n i ) u ( θ i ) + 2 D x σ i ( t ) − γ � � ˙ θ i = sin( θ i − θ j ) + 2 D θ η i ( t ) n i j neighbor of i with n i = local density σ i , η i = gaussian white noises , unit covariance Representative of a class of models with similar large scale properties.

  6. Qualitative behavior • Strong interactions, or ”external field” → local orientation order. Not studied here. • Weak interactions → no local orientation order Large scale dynamics = diffusive. • v depends on the local density ρ → effective diffusion coefficient depends on ρ Possibility of ”motility induced phase separation” (Cates, Tailleur)

  7. Main questions • A macroscopic description? Finite N fluctuations? Stationary measure? Probability distribution of the density? • A very quick review ◮ J. Toner, Y. Tu (1995): phenomenological hydrodynamical equations + noise introduced ”by hand” ◮ E. Bertin, M. Droz, G. Gr´ egoire (2006): write a Boltzmann like equation + expansion close to the phase transition threshold → derivation of Toner-Tu like equations, without noise (many developments from there: Chat´ e et al., Marchetti et al., Ihle . . . ) ◮ Math. literature: P. Degond, S. Motsch (2007); Fokker- Planck like models (locally mean-field); far from the threshold ◮ Keeping finite N fluctuations: J. Tailleur, M. Cates et al. (2008, 2011, 2013): without alignment promoting interactions; Bertin et al. (2013): derive a noise from the microscopic equations for nematics.

  8. Our goals 1. start from microscopic equations 2. derive hydrodynamical equations and noise in a controlled way Noise may have correlations → important to have a microscopic derivation 3. exploit these results to study the dynamical fluctuations of the empirical density (cf Macroscopic Fluctuation Theory). 4. obtain large deviation estimate for the stationary spatial density ρ such as P ( ρ ≈ u ) ≍ e NS [ u ] S = ”entropy”, or ”quasi-potential”. Simple framework: aligning interactions below threshold for local order; density dependent speed ( → clustering possible).

  9. Microscopic equations (simplified), adimensionalized � d ˜ x i 2˜ σ i (˜ = ε ˜ v ( n i ) u ( θ i ) + ε D x � t ) (1) d ˜ t √ d θ i − ˜ γ � 2 η i (˜ = sin( θ i − θ j ) + t ) , (2) d ˜ t n i j neighbor of i with ε = v 0 / ( LD θ ), ˜ D x = D x D θ / v 2 γ = γ/ D θ , ˜ 0 , ˜ t = D θ t . Two important parameters: ˜ D x : ratio spatial diffusion/”active” diffusion ˜ γ : strength of the aligning interaction ε = spatial time scale/angular time scale: small parameter

  10. Strategy Main object of interest: the empirical density ρ ( x , θ, t ) = 1 � δ ( x − x i ( t )) N i Phase space empirical density f ( x , θ, t ) = 1 � δ ( x − x i ( t )) δ ( θ − θ i ( t )) N i 1. Write an equation for f that keeps finite N fluctuations (cf D. Dean 1996): in a sense exact in the large N limit 2. Use the time-scale separation to write an equation for ρ that keeps finite N fluctuations: hoped to be exact in a combined ε → 0 , N → ∞ limit 3. Write a functional Fokker-Planck equation for µ t [ ρ ], the pdf of ρ . 4. Look for a stationary solution of the form µ [ ρ ] ≍ e NS [ ρ ]

  11. A fluctuating non linear Fokker-Planck equation interaction transport � �� � � � � ∂ f γ ∂ � �� � d θ ′ sin( θ − θ ′ ) f ( θ ′ ) ∂ t = − ε ∇ ( v ( ρ ) u ( θ ) f ) + f ρ ∂θ √ 2 D x � √ √ 2 ∂ � � � � + η ( x , θ, t ) f + ε √ ∇ x · � σ ( x , θ, t ) f N ∂θ N � �� � finite N fluctuations ∂ 2 f ∂θ 2 + ε 2 D x ∇ 2 + x f � �� � angular and spatial diffusions Meaning? A dynamical large deviation principle (Dawson 1987). � T P ( f t ≈ g t ) ≍ exp( − NJ [0 , T ] [ g ]) ; J [0 , T ] [ g ] = 1 || ∂ t g − VFP [ g ] || 2 − 1 , g dt 4 0 VFP = nonlinear Vlasov-Fokker-Planck operator= red terms

  12. On the computations √ • Local equilibrium + small deviation (order ε and 1 / N fluctuations) f ( x , θ, t ) = 1 2 πρ ( x , ε α t ) + δ f ( x , θ, t ) • Equation for ρ : slow time scale, depends on δ f √ 2 D x � ∂ρ � � u θ δ f ) + ε 2 D x ∇ 2 ρ + ε � √ ∂ t = − ε ∇ ( v ∇ · ξ ( x , y , t ) (3) N ξ = noise, multiplicative in ρ . • δ f small → obtained by solving a linearized equation • Reintroduce into Eq.(3) → the final equation, a fluctuating PDE for ρ .

  13. Dynamical large deviation principle • Fluctuating PDE for ρ ∂ρ 1 = U [ ρ ]( � x ) + √ ν ( � x , t ) ∂ t N � � 1 v ( ρ ) + D x ∇ 2 ρ U [ ρ ]( � x ) = 2 ∇ · ∇ [ v ( ρ ) ρ ] 1 − ¯ γ 2 � ν ( x , y , t ) ν ( x ′ , y ′ , t ′ ) � x ′ ) δ ( t − t ′ ) = D [ ρ ]( � x ,� • The fluctuating PDE for ρ is a rephrasing of a dynamical large deviation principle ”` a la Dawson” � T P ( ρ ≈ u ) ≍ exp( − NI [0 , T ] [ u ]) with I [0 , T ] [ u ] = 1 || ∂ t u − U [ u ] || 2 − 1 , D 2 0 • This kind of dynamical large deviation principle is the starting point for the macroscopic fluctuation theory (in this case, it is actually trivial . . . )

  14. Yet another formulation: functional Fokker-Planck equation • Ordinary stochastic differential equation for x ∈ R d → PDE (Fokker-Planck) for the pdf of x . • Stochastic PDE for a field ρ → functional equation for µ t [ ρ ], ”pdf” of ρ . drift part � �� � � ∂µ t δ = − d � x x ) ( U [ ρ ]( � x ) µ t ) δρ ( � ∂ t �� � � 1 δ δ x ′ D [ ρ ]( � x ′ ) + d � x d � x ,� x ′ ) µ t δρ ( � δρ ( � 2 N x ) � �� � diffusion part

  15. Results and discussion • When γ < γ c , system effectively at equilibrium (case without aligning interactions: Cates, Tailleur et al. 2007, 2011, 2013) � → computing S is possible, S [ ρ ] = s ( ρ ( x )) d x � � v 2 ( ρ ) + ρ v ( ρ ) v ′ ( ρ ) + 2 D x s ”( ρ ) = − � � 1 − ¯ γ b [ ρ ] b [ ρ ] 2 • → compute S [ ρ ]

  16. Results and discussion • Reasonable to assume v ( ρ ) decreasing → possible phase separation (MIPS = Motility Induced Phase Separation). Role of the interactions? Disorder Homogeneous sp Order D x MIPS 0 0 1 2 3 γ / D θ Sketch of the D x − ˜ γ phase diagram.

  17. Results and discussion Left: entropy s ( ρ ). Right: D x − ρ phase diagram. ◮ Spinodal line very sensitive to the interaction strength (observed in simulations) ◮ Density fluctuations increase when approaching the ordered phase ◮ A strong enough spatial diffusion always prevent phase separation

  18. Conclusion ◮ Nice example where the limiting procedures seem well controlled + a general strategy ◮ Some physical insight in the ”Motility Induced Phase Separation” with aligning interactions ◮ Next step: with a local orientation order → hyperbolic hydrodynamic limit → one cannot expect an effective equilibrium in the same sense ◮ Mathematical theory much less developed in this case . . . (recent works by Mariani, Bertini et al.) Work in progress . . .

Recommend


More recommend