Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
  • Featured in Physics
  • Editors' Suggestion

Pressure and Phase Equilibria in Interacting Active Brownian Spheres

Alexandre P. Solon, Joakim Stenhammar, Raphael Wittkowski, Mehran Kardar, Yariv Kafri, Michael E. Cates, and Julien Tailleur
Phys. Rev. Lett. 114, 198301 – Published 11 May 2015
Physics logo See Viewpoint: An Equation of State for Active Matter
PDFHTMLExport Citation

Abstract

We derive a microscopic expression for the mechanical pressure P in a system of spherical active Brownian particles at density ρ. Our exact result relates P, defined as the force per unit area on a bounding wall, to bulk correlation functions evaluated far away from the wall. It shows that (i) P(ρ) is a state function, independent of the particle-wall interaction; (ii) interactions contribute two terms to P, one encoding the slow-down that drives motility-induced phase separation, and the other a direct contribution well known for passive systems; and (iii) P is equal in coexisting phases. We discuss the consequences of these results for the motility-induced phase separation of active Brownian particles and show that the densities at coexistence do not satisfy a Maxwell construction on P.

  • Figure
  • Figure
  • Received 18 December 2014

DOI:https://doi.org/10.1103/PhysRevLett.114.198301

© 2015 American Physical Society

Viewpoint

Key Image

An Equation of State for Active Matter

Published 11 May 2015

An equation of state for a gas of self-propelled spheres is a step towards a thermodynamic description of “active” matter, such as bird flocks and tissue.

See more in Physics

Authors & Affiliations

Alexandre P. Solon1, Joakim Stenhammar2, Raphael Wittkowski2, Mehran Kardar3, Yariv Kafri4, Michael E. Cates2, and Julien Tailleur1

  • 1Laboratoire, Matière et Systèmes Complexes, UMR 7057 CNRS/P7, Université Paris Diderot, 75205 Paris Cedex 13, France
  • 2SUPA, School of Physics and Astronomy, University of Edinburgh, Edinburgh EH9 3FD, United Kingdom
  • 3Department of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
  • 4Department of Physics, Technion, Haifa 32000, Israel

Article Text (Subscription Required)

Click to Expand

Supplemental Material (Subscription Required)

Click to Expand

References (Subscription Required)

Click to Expand
Issue

Vol. 114, Iss. 19 — 15 May 2015

Reuse & Permissions
Access Options
CHORUS

Article Available via CHORUS

Download Accepted Manuscript
Author publication services for translation and copyediting assistance advertisement

Authorization Required


×

Images

  • Figure 1
    Figure 1

    Numerical measurements of P0+PI, PS, and PD in single-phase ABP simulations at Péclet number Pe3v0/(Drσ)=40, where σ is the particle diameter. Expressions (12), (14), and (16) for P0+PI and PS show perfect agreement. Also shown are data for Pe=20, unscaled and rescaled by factor 2. This confirms that PS=P0+PI is almost linear in Pe; small deviations arise from the Pe dependence of the correlators. In red is PD for Pe=20,40, with no rescaling. Pe was varied using Dr, at fixed v0 and with Dt=Drσ2/3. Solid lines are fits to piecewise parabolic (PS) and exponential (PD) functions used in the semiempirical equation of state. ρ0 is a near-close-packed density at which v(ρ) vanishes and ρ˜ is the threshold density above which PD>PS. See the Supplemental Material [39] for details.

    Reuse & Permissions
  • Figure 2
    Figure 2

    Simulated coexistence curves (binodals) for ABPs (red) and those calculated via the Maxwell construction (black) on the mechanical pressure P using the semiempirical equation of state for PS and PD fitted from Fig. 1. Dashed lines: predicted high Pe asymptotes for the binodals calculated via f or Pf (lower line) and calculated via P or fP (upper line). Inset: measured binodal pressures and densities (diamonds) fall on the equation-of-state curves but do not match the MC values (horizontal dashed lines). Stars show the P(ρ) relation across the full density range from simulations at Pe=40 and Pe=100. The latter includes two metastable states at low density (high ρ0/ρ) that are yet to phase separate.

    Reuse & Permissions
×

Sign up to receive regular email alerts from Physical Review Letters

Log In

Cancel
×

Search


Article Lookup

Paste a citation or DOI

Enter a citation
×