Partial and complete wetting of droplets of active Brownian particles

We study wetting droplets formed of active Brownian particles in contact with a repulsive potential barrier, in a wedge geometry. Our numerical results demonstrate a transition between partially wet and completely wet states, as a function of the barrier height, analogous to the corresponding surface phase transition in passive fluids. We analyse partially wet configurations characterised by a nonzero contact angle $\theta$ between the droplet surface and the barrier, including the average density profile and its fluctuations. These findings are compared with two equilibrium systems: a Lennard-Jones fluid and a simple contour model for a liquid-vapour interface. We locate the wetting transition where $\cos(\theta)=1$, and the neutral state where $\cos(\theta)=0$. We discuss the implications of these results for possible definitions of surface tensions in active fluids.


I. INTRODUCTION
Active matter systems display a wide range of surprising phenomena in their non-equilibrium steady states.Among the simplest active systems are fluids comprised of self-propelled particles without aligning interactions, which are known as scalar active matter [1][2][3].Many such systems exhibit motility-induced phase separation (MIPS) [4], which resembles equilibrium coexistence of dense and dilute fluid phases.However, the interfacial properties of these active fluids phases differ significantly from their equilibrium counterparts which can lead -for example -to microphase separation and bubbly phases [5].
Given such observations, it is natural to ask about other interfacial properties of active fluids, and their similarities and differences with equilibrium systems.An interesting example occurs when a system undergoing MIPS is placed in contact with a solid or penetrable substrate (or "wall").In this case, one may expect analogues of the rich phenomenology of wetting, as it occurs in equilibrium fluids at liquid-vapor coexistence [6].Processes reminiscent of equilibrium wetting appear to play a crucial role in active systems composed of living cells, soft responsive materials, and embedded energy sources.For example, wetting and dewetting on soft substrates enables tunable adhesion, motility, and shape change of cells [7][8][9][10] and control the motility of bacteria at interfaces [11].
In equilibrium, wetting behaviour can be analysed in several different settings.A famous example is the formation of a liquid droplet on a weakly attractive solid substrate.The contact angle θ of this droplet obeys Young's equation: where γ lv is the liquid-vapor surface tension and similarly γ wv , γ wl are surface tensions between the fluid phases and the wall.Increasing the attraction between the fluid and the wall, the tension γ wl decreases and the droplet spreads out, leading eventually to a wetting transition [12,13] as θ → 0 (specifically, this is the transition from partial to complete wetting, but we term it here "the wetting transition", for simplicity).This transition -and related phenomena such as drying transitions [14,15] may be either first-order or critical, depending on the behaviour of cos θ as the transition is approached.Measurements of contact angle thus permit the characterisation of wetting transitions on planar surfaces [16] as well as related situations such as the filling transitions that occur in capillaries [17].For computational model fluids like Lennard-Jonesium, measurement of contact angles tends to be challenging.However, there are convenient alternative approaches which either exploit the grand canonical ensemble, or a slit geometry with a fluid confined between two walls (and periodic boundaries in the other direction).Droplets do not form in these cases, and one instead focusses on the average density profile ρ(z), as a function of the distance z from the wall.For a grand canonical system with a single wall, one defines the adsorption Γ ≡ ∞ 0 dz (ρ(z) − ρ b ) where ρ b is the bulk density at z → ∞.This quantity is accessible experimentally [18,19] as well as in density functional theories [15,20] and in simulations [20,21] (with the aid of finite-size scaling).Increasing the attraction between the fluid and the wall, a wetting transition occurs when Γ becomes infinite, which may occur either by a smooth divergence (critical wetting) or by a discontinuous jump (first-order wetting), depending on the range of fluid-fluid and wall-fluid interactions.A similar analysis can be performed for the slit geometry in the canonical ensemble, in which case the wetting transition is signalled by a symmetry breaking of the density profile ρ(z) with respect to the two walls [21,22].
In all these equilibrium cases, statistical mechanical theories place strong constraints on the phenomenology.For example, the surface tensions in (1) can be defined unambiguously through gradients of an appropriate free energy, as can the adsorption Γ [15,20].This provides consistency requirements between different ensembles and geometrical settings: studies based on the adsorption and the contact angle both deliver the same results for the locations and properties of surface phase transitions, as long as finite-size effects are controlled.By contrast, active fluids are not ruled by a free energy, and the status of their wetting transitions is much less well-understood.Indeed, there are several different proposals for active generalisations of the liquid-vapour surface tension [23][24][25][26][27][28][29][30][31][32][33].It is not clear a priori whether any suitably generalised version of Young's equation should apply for these systems; if some such generalisation does exist then one may ask which (if any) of the liquid-vapour surface tensions might appear, and what should be used in place of γ wv − γ wl .
In recent work by some of us [22], an approach based on the absorption Γ was used to analyse the wetting properties of a paradigmatic active fluid, comprised of active Brownian particles (ABPs) in d = 2 and d = 3 dimensions.An important difference from equilibrium fluids is that active particles tend to accumulate at walls, even in the absence of attractive interactions [34,35].As a result, an infinitely repulsive 'hard' wall is always wet for these active fluids, in contrast to equilibrium fluids for which a hard wall remains dry [36].However, on replacing a solid wall with a penetrable barrier, behaviour similar to first-order wetting was found in d = 3, for the slit geometry [22].
In this work, we take a complementary approach to Ref. 22, which is to examine the wetting behaviour of droplets in the same 3d system of ABPs.Such studies are numerically challenging due to finite-size effects which appear in the form of large fluctuations of the droplet shape and position.We show that this can be mitigated by confining droplets in a wedge geometry.We identify a discontinuous wetting transition for droplets (where cos θ → 1).The nature and location of this transition are consistent with the slit geometry.We also characterise the situation of neutral wetting (cos θ = 0).We give a critical analysis of Young's equation in this setting.We argue that by defining surface tensions in terms of the probabilities of droplet shape fluctuations, Young's equation holds by definition for the most likely shape, as long as shape fluctuations are controlled by local properties of the interface.In this case, the liquid-vapour surface tension also determines the probability of large-wavelength capillary waves.
Our paper is organised as follows: the ABP model is defined in Sec.II and its wetting behaviour is discussed extensively in Sec.III.Then Sec.IV discusses analogous behaviour in a passive LJ fluid.Sec.V introduces the simple contour model for equilibrium interfacial fluctuations and compares it with the particle-based models.These results are discussed in Sec.VI, which also sum-marises our main conclusions.

II. MODEL AND GEOMETRICAL SETUP
We analyse ABPs as a prototypical model of scalar active matter that displays MIPS in the bulk [37,38].We follow a previous parameterization [39] described in Appendix A, with ABPs interacting via the Weeks-Chandler-Anderson potential of lengthscale σ, and energy scale ε, with coupled translation and rotational diffusion constants D r = 3D t /σ 2 and self-propulsion velocity v 0 defining a Péclet number Pe = v 0 /σD r .
We simulate this model in a fully periodic orthorhombic box of dimensions L x × L y × L z .We consider various box sizes, but we typically elect to work with L z = 20σ, smaller than both L x and L y .In order for a surface phase transition to occur, the system must be at a state point for which bulk liquid-vapor coexistence occurs.Accordingly, we choose model parameters well inside the MIPS region, specifically number density ρ = 0.60 and constant Péclet number Pe = 60 (see [40] for the bulk phase diagram).Applying the lever rule, this corresponds to an approximate liquid fraction of the system f = (ρ − ρ LD )/(ρ HD − ρ LD ) ≈ (0.6 − 0.45)/(1.25 − 0.45) = 0.1875, which results in a cylindrical liquid domain with its axis parallel to the z axis.Such a geometry allows us to monitor the changes in the liquid-vapour interfaces (and, in particular, the contact angle) via its twodimensional projections in the x − y plane.We note that while MIPS is unstable with respect to vapor-crystal phase separation in this model [40,41], crystallisation is not observed on our simulation timescales.
To induce wetting, we introduce a penetrable wall, which we model as a static repulsive barrier [22].To this end, we identify a piecewise surface S and employ an external repulsive potential perpendicular to the surface, which takes the form of a cosine hump: where r is the perpendicular distance from the surface S.
We take d = σ, which corresponds to a short-ranged (ie.thin) potential barrier.Simulations of liquid droplets of ABPs are challenging because of large fluctuations of the droplet shape.This can be mitigated by increasing the system size, but the dilute (vapor) phase of MIPS has a relatively high concentration ρ LD ≈ 0.45, which means that such simulations quickly become expensive, involving very large numbers (hundreds of thousands) of particles.While planar walls are natural for wetting, we have found it helpful to choose S to take a wedge shape; this localises the droplet (reducing fluctuations) and accelerates its nucleation.To achieve this, we choose S to comprise two finite   planes that are joined together at an angle α by a cylindrical section, which gives the wedge a rounded corner.Fig. 1(a-c) illustrate the resulting setup, showing the planar density ρ(x, y) (details of its numerical estimation are given below).The x − y force field generated by the barrier is illustrated in Fig. 1(d).Typical simulations have L x = L y = 100σ and L z = 20σ with the wedge occupying one quadrant of the box.At the density considered, this results in simulations of N = 120 000 particles.Working in three dimensions ensures that the density fluctuations are more controlled than the two-dimensional case and connects to previous evidence for a wetting transition that becomes sharper in the large N limit [22].

III. WETTING PHENOMENOLOGY FOR ACTIVE BROWNIAN PARTICLES
A. Average local density profiles: contact angles and wetting transition Our analysis of wetting is based on the local density profile, which is obtained by discretizing space in subvolumes of side ℓ = 2σ.We define the instantaneous number of particles in bin j of center b j where I ℓ (x) is 1 if x is within the cube of side ℓ centered at the origin and zero otherwise.The instantaneous number density at site j is defined as and we average along the z dimension to obtain twodimensional density maps ρ(x, y) = L −1 z Lz 0 ρ(x, y, z)dz.Following a relaxation time of 100τ r , the average local density ρ(x, y) is formed by sampling at intervals of τ r , and averaging over a period of 1500τ r .Fig. 1 shows two-dimensional projections of the density map for three different wedge aperture angles α.The density profiles are smooth and -for the value of ε w = 12 studied-exhibit a liquid drop confined within the wedge.From such profiles one can, in principle, extract the vapor-liquid interface and estimate the contact angle θ between the active droplet and the repulsive barrier.For our system, the macroscopic notion of sharp contact between a circular vapor-liquid interface and the barrier is somewhat blurred by the ubiquitous presence of a thin layer of particles all around the wedge.To deal with this, we work with estimates of an apparent contact angle, defined by fitting a circular arc to circular regions of the density profile as described further in Appendix C. If a non-equilibrium analogue of Young's equation applies to these systems (valid in the large-system limit), then θ should be determined only by local properties of the three-phase contact line, so that θ is independent of α.For the three value of α shown in Fig. 1, we find θ ≈ 75 • −85 • : this relatively small range seems consistent with the applicability of Young's equation.Note that for small equilibrium droplets Young's equation is modified by line tension and by Tolman corrections to the Laplace pressure [42][43][44], while in the active case curvature effects on swim pressure of confined ABPs have also been reported [45]; we explicitly neglect these effects in this work.
Henceforth we elect to work at fixed aperture α = π/2.Fig. 2 illustrates how the stationary average density field ρ(x, y) for this value of α varies with the barrier strength ε w .At the largest values of ε w studied, thick liquid layers are present on both the interior and exterior of the wedge and the density profile flattens close to the corner of the wedge and terminates with a rounded shape at the tips.On lowering the barrier height from ε w = 18 to ε w = 15, liquid progressively accumulates in the wedge interior.At ε w ≈ 18 the curvature of the interface between the interior liquid and the vapor changes sign from positive to negative and a recognisable droplet forms within the wedge having an apparent contact angle θ < π/2.The fitting procedure to obtain θ is appropriate only when    the density profile exhibits such a region of negative curvature, ie. for ε w ≲ 18.We assert (and confirm via a comparison with equilibrium wetting in a comparable geometry-see Secs.IV and V B) that the change in sign of the interfacial curvature corresponds to the transition from partial to complete wetting.In other words that the barrier is partially wet, (with θ > 0) for ε w ≲ 18, and is completely wet (ie.θ = 0) for ε w ≳ 18.Within the partially wet regime, most liquid resides in the drop in the wedge interior.However some liquid resides on the exterior wall of the wedge forming a pair of symmetrical "lobes" of liquid-like density.A ε w decreases, progressively more liquid accumulates in the drop whose contact angle θ increases, while the extent of the lobes decreases.In Sec.III D we discuss the finite-size scaling behaviour of the lobes and the liquid drop.At ε w = 10 we find θ ≈ π/2, which -in the context of Young's equation -is interpreted as a neutral point, where the tensions between the barrier and liquid and the barrier and the vapor balance each other.This point separate the partial wetting and partial drying regimes.On decreasing the barrier strength still further, we find that within the timescale of our simulations, the stationary state becomes harder to define: if we initiate the system in a homogeneous density state, the nucleation of the droplet becomes very slow; if we start instead from a pre-formed droplet and instantaneously decrease the barrier strength, the droplet progressively detaches from the barrier.The detachment indicates that the vapor phase is favoured near the barrier in the stationary state.Henceforth, we restrict our analysis to the range of barrier strength for which the liquid is attached to the barrier, where surface physics plays a role.
Our results for the dependence of the measured contact angle on the barrier strength are displayed in Fig. 3.For the weakest accessible barrier strengths at which the droplet is attached to the wedge, we find cos(θ) ≈ 0, complete wetting partial wetting detached FIG. 3. Cosine of the apparent contact angle θ for ABPs in a 90 • wedge for a range of repulsive strengths εw, as calculated from a circular arc fit to the portion of the liquid-vapor interface having negative curvature.The vertical dashed line indicates the transition from partial to complete wetting determined previously via a different approach [22].The hatched area indicates the regime in which the droplet detaches from the barrier.Error bars are bootstrap estimates for three standard deviations.
i.e. the droplet is close to the neutral point.As ε w increases, cos(θ) increases approximately linearly and appears to attain the wetting point cos θ = 1 with a nonzero slope, suggesting a first-order wetting transition around ε w = 18 ± 1.This value of the wetting point accords with our previous analysis of this system in the slit geometry [22].

B. Local density fluctuations
Density fluctuations play a central role in characterising phase transitions, both in passive and active systems.In the case of surface phase transitions studied within a grand canonical framework, the behaviour of the local compressibility (which quantifies the magnitude of the density fluctuations relative to the average density) allows the order of a surface phase transition [15] to be ascertained and provides information on the character of interfaces.
To quantify the density fluctuations in our system, we accumulate the local variance of the density field Var[ρ(x, y)] = ⟨ρ 2 (x, y)⟩ − ⟨ρ(x, y)⟩ 2 on the same scale ℓ = 2σ over which the field is defined.This provides a fine-grained description of the spatial dependence of the density fluctuations.Fig. 2 shows that the fluctuations are greatest at the liquid-vapor interface: both for the droplet in the wedge interior and for the exterior lobes.In particular, the complete wetting regime at large ε w corresponds to small overall fluctuations, with higher values occurring around points of higher curvature.As the repulsive barrier gets weaker (moving from right to left in Fig. 2), the fluctuations at the liquid-vapor interface in the interior increase in magnitude, and become progressively more localised in the vicinity of the contact region between the interface and the barrier.

C. Polarisation field
The orientation of the self-propulsion force of ABPs rotates via a diffusive process.For a bulk system, every particle orientation is equally likely, and thus, the average orientation vector of each particle is zero.When interfaces are formed (as in MIPS), the local stationary orientation can take nonzero values, indicating a local polarisation of the system.
We define the instantaneous polarisation ϕ from the particle orientation n i on a grid of spacing ℓ.For bin j which corresponds to the local average of the orientation field.Its norm is |ϕ| ∈ [0, 1], and it is expected to be negligible in the bulk phases, as the individual orientations diffuse on the unit sphere.Near boundaries and interfaces, however, the net orientation typically follows the density gradient, i.e. normal to the interface.In Fig. 4(a-d), we track the changes to the stationary polarisation that occur as we vary the barrier strength ε w .As expected, the polarisation fields are non-zero only at interfaces and closely follow the density gradient, directed from the dilute towards the denser phase.The principal effect of varying ε w is manifest in the barrier region: For large ε w , the particles at liquid-like densities on both the interior and exterior of the wedge are oriented against the barrier.For weaker barriers (e.g.ε w = 12), the particles inside the droplet change orientation and point towards the droplet interior.This contrasts with the behaviour of exterior particles in the vicinity of the barrier, which are always polarised towards the barrier, independent of ε w .

D. Finite size effects
The density profiles display the formation of both an interior liquid layer and some accumulation on the exterior of the wedge.These exterior density "lobes" are a consequence of the non-equilibrium accumulation mechanism and its coupling with the barrier strength.In the limit of very repulsive barriers, the liquid covers both the interior and exterior approximately symmetrically, whereas, in the partial wetting regime that occurs for ε w < 18, the liquid is primarily localised as a droplet in the interior.
If we now consider increasing the overall system size at partial wetting, whilst maintaining the overall number density constant, then one expects the volume of liquid in the wedge interior to increase in size accordingly.However, the thickness of the exterior lobes should remain approximately unchanged because they arise from the usual accumulation of ABPs at a barrier, combined with a local balance of particles crossing the barrier: these aspects only depend on the liquid density.
Fig. 5 confirms this: the top panel shows that on increasing the repulsion strength, the thickness of the lobes grows, and eventually saturates as the complete wetting regime is approached.(This thickness is measured as the stationary average of the largest distance between the outer liquid-vapor interface and the exterior of the wedge).The lower panel of Fig. 5 shows that in the partially wet regime ε < 18, the lobe thickness is independent of system volume.Thus, the exterior lobes in the partial wetting regime can be regarded as a finitesize effect such that for a sufficiently large system, the liquid phase is essentially confined in the wedge interior.

E. Flow fields
The internal self-propulsion force renders the system dissipative, promoting flows coupled to the geometric features of the barriers.For example, in the case of ellipsoidal impenetrable barriers, distinctive quadrupolar flow patterns have been reported previously [46].
In order to quantify the flow field, we track the particle displacements ∆r i = r i (t+∆t)−r i (t) with ∆t = 0.005τ r , a time interval that is sufficiently small to allow us to follows the local patterns of motion.We estimate the flow patterns from the α = x, y, z components of the field where êα is the unit vector.We average over the steady state and the z direction to obtain two-dimensional flow maps from which we can extract streamlines.Here the end-points of the wedge correspond to regions of opposing counterflows, reminiscent of the quadrupolar structure of hard objects mentioned above; the flow in the interior of the wedge is from the vapor into the liquid, consistent with the orientation field; particle flow through the wedge occurs only at the wedge corner, where two other counter-flows are formed, and where the flow field rotates in order to be orthogonal to the exterior of the wedge.These overall features persist on decreasing ε w , but the flow patterns become more complex.In particular, the flow field reorients itself around the contact points between the droplet and the barrier, forming two additional local regions of nonzero circulation.Softer barriers also engender interesting changes in the flow patterns inside the liquid droplet.One is a marked increase in the flow through the wedge corner.Another is that while for the largest ε w particles are trapped in the liquid region, with any flow occurring parallel to the barriers (for ex-ample, Fig. 4(h)) in the positive x and y directions, for weaker ε w the particles can traverse the barrier, creating perpendicular flow patterns as seen in Fig. 4(e).These changes are accompanied by a reversal in the net flux direction, with x and y flux components being net negative for small ε w and net positive for large ε w .Forces on asymmetric objects have been previously linked [47] to the flow patterns of active particles, and the changes observed here indicate a change in the direction of the net force on the wedge as we vary the repulsion strength, a prediction that can potentially be tested in experiments.
These results emphasise that the dissipative nonequilibrium flow patterns generated by the soft wedge can be highly non-trivial, with exquisite features that start to develop even well inside the dilute phase.Notwithstanding these fine-grained phenomena, as we show in the next section, coarse-graining over sufficiently large scales reproduces closely the principal features of the density profiles and their fluctuations that occur in comparable equilibrium systems.

IV. COMPARISON WITH EQUILIBRIUM WETTING IN A PASSIVE FLUID
The phenomenology of the active system displays specific spatial density profiles, local density fluctuations and flow patterns.While the non-equilibrium flows have no counterpart in equilibrium, we can compare the stationary density distributions and their fluctuations to a conventional passive system undergoing a wetting transition in an analogous geometry.To do so, we consider the prototypical model for equilibrium wetting in which a liquid droplet of Lennard-Jones particles, in coexistence with the vapor, interacts with an impenetrable wall supplemented by a long ranged attraction, as described by a Lennard-Jones 9-3 interaction potential.Such longranged wall-fluid attractions are necessary for a planar substrate to engender a first-order wetting transition [20].We note that for our active system, the competition between attraction and repulsion that engenders a wetting transition arises from a single repulsive potential barrier alone.This contrasts with passive systems, where the necessary ingredients are separate attractive and repulsive parts of the wall-fluid potential.
To make contact with the behaviour of droplets in our active system, we arrange the wall to form a wedge with aperture angle α = 90 • (see details of the equilibrium model in Appendix B).We control the affinity of particles for the wall by varying the strength of the interaction ε LJ w between the wedge and the fluid.Note, however, that in contrast to the active case where increasing ε w increases the degree of repulsion between the barrier and the fluid particles, here increasing ε LJ w strengthens the wall-fluid attraction.
Fig. 6 shows the evolution of the density profile with increasing attractive strength ε LJ w .The boundary conditions of the simulation box are such that the LJ walls are located on the left and bottom edges while the top and right edges present reflective boundary conditions; the z dimension, orthogonal to the projection plane, is periodic.The qualitative dependence on wall-fluid interaction strength is very similar to that seen for the active case in Fig. 2: a droplet is formed in the corner of the wedge, and gradually spreads out as the attractive strength is increased.Complete wetting occurs at the largest wall strengths where the particles coat the wall (note that the rounding of the density profile at the end of the simulation box is a consequence of the reflective boundary conditions required in this ensemble and is not a genuine point of contact).As in the active case, the changes in the density profiles are accompanied by changes in the local density fluctuations, as shown in Fig. 6.The regions of highest variance (greatest fluctuation) are located at the vapor-liquid interface, and on approaching the neutral point, θ = π/2, the variance exhibits maxima near the point of contact with the wall.
When the attractive wall strength is reduced further, a situation similar to that observed in the active case arises: the liquid drop pinches off and detaches from the wall.Fig. 7 plots cos(θ) as a function of ε LJ w as extracted from circular fits to the interface.The main difference to the active case is that the passive system can access a larger range of cos(θ) before detachment occurs: specifically for the passive system we can stabilize droplets with cos θ ≈ −0.5 (θ ≈ 2π/3), whereas the active system becomes unstable around the neutral point θ = π/2.

V. THEORETICAL MODEL FOR DROPLET SHAPE FLUCTUATIONS
On sufficiently large (hydrodynamic) length scales, it is natural to describe wetting droplets in a continuum formalism, with a sharp interface separating liquid and vapour domains.We discuss here a simple model for the fluctuations of such an interface, which is useful for rationalising the density fluctuations observed in active and passive systems.We describe the interface as a line in two dimensions, in order to model the planar profiles in Fig. 2.

A. Contour Model
For a wedge with opening angle α, we describe the interface in polar co-ordinates, with the origin at the corner of the wedge.The distance from the origin to the inter-   " LJ w = < l a t e x i t s h a 1 _ b a s e 6 4 = " w E F C B q X T l p P z H   face R(ϕ) is parameterised as

< l a t e x i t s h a 1 _ b a s e 6 4 = " T x G v B + N O t n P h H c 4 c 7 f L h h f h d 3 R 0 = " >
where ϕ is the polar angle with 0 ≤ ϕ ≤ α and the parameter M enforces a short-wavelength cutoff for numerical convenience.Example interfaces are shown in Fig. 8a, they are parameterised by R 0 , B, a 1 , a 2 , . . ., a M .
The following results depend weakly on the cutoff M .We formulate our model in terms of the probability distribution for the droplet shape, encoded by the function R(ϕ).Since we describe the interface as a lineor contour -we refer to this as the contour model.For simplicity, we suppose that the area of the liquid region is fixed at some reference value A 0 , so the probability (or probability density) of a given contour is where 2 dϕ is the area enclosed by the liquid-vapour interface, F[R] is a free-energy-like quantity that controls the shape fluctuations, and Z is a suitable normalisation constant.
By analogy with equilibrium systems we propose that the dominant terms in F on hydrodynamic scales are where ℓ lv is the length of the liquid-vapour interface, d wl is the (total) length of the liquid-wall interface, and γlv , ∆γ w play the role of surface tensions (see below for a more detailed discussion).One sees immediately that d wl = R(0) + R(α) and the interfacial length is In equilibrium, ( 9) is valid on hydrodynamic scales and its parameters are related to the surface tensions between the phases as γlv = γ lv /(k B T ) and ∆γ w = (γ wv − γ wl )/(k B T ): here T is the temperature and k B is Boltzmann's constant. 1In principle other terms might also appear in F but these are constrained by localitythe probability of an interfacial perturbation should not depend on the shape of the interface far away -and by the symmetries of the problem, and the restriction to large length scales.
For active systems, non-local contributions to F cannot be ruled out, but the symmetries of the system (isotropy and translation invariance) suggest that ( 9) is a natural starting point for a hydrodynamic theory.Moreover, for large droplets one sees that ℓ lv and d wl are O(L) and γlv , ∆γ = O(1) which means that shape fluctuations are small, as happens for thermal fluctuations in the thermodynamic limit.As in that case, the result is that the distribution of shapes is dominated by the most likely contour.Applying calculus of variations to minimise F at fixed droplet area A 0 one arrives as Young's equation in the form γlv cos θ = ∆γ w . ( For completeness, a derivation of this result is given in appendix D, which also shows that the curvature of the liquid-vapour interface must be constant everywhere, so that the droplet boundary is an arc of a circle.Analysing small fluctuations about this most-likely shape also recovers the standard theory of capillary fluctuations, with a spectrum proportional to q2 /γ lv .This suggests that the appropriate γ lv that should appear in a Young's equation for active fluids is a capillary surface tension [31], see Sec.VI below for further discussion of this point.

B. Comparison with wetting simulations
To investigate the predictions of this model, we sample contour fluctuations numerically by a straightforward Monte Carlo algorithm as detailed in appendix E. Assuming sharp interfaces, this yields a density field ρ(x, y|t) by taking a value of unity for the liquid and zero for the vapor.These sharp profiles are then averaged to return the equilibrium density profile, and its variance, which are both computed by partitioning the 2d domain into bins of size δx.
Figure 8 illustrates the procedure for a particular choice of ∆γ/γ, connecting collections of paths (Fig. 8(a)) to a density profile and its variance (panels (b) and (c)).For a specific choice of ∆γ/γ the average density profile converges to form a specific contact angle, such that cos θ = ∆γ/γ: in the specific example of panel (b) and (c) we have ∆γ/γ = 0.7 → θ ≈ 45 • .While the density profile closely resembles that produced by the passive system (and the interior of the wedge in the active system), the variance profiles are typically sharper and do not display a marked relative increase of the variance in the regions where the contour contacts the wall.
This deficiency of the model arises from our assumption of perfectly sharp instantaneous interfaces.In reality, the density across the instantaneous interfaces (e.g. in the Lennard-Jones system) evolves smoothly, interpolating between the liquid and the vapor density with a characteristic interface width.(The observed width is also affected by the projection of the three-dimensional system onto a planar density, because of capillary waves along the z-direction.)To account for this width, we perform a convolution of the instantaneous density profiles with a Gaussian kernel of size σ blur .As illustrated in Fig. 9, the convolution promotes the emergence of peaks in local density fluctuations in the vicinity of the contacts with the wedge, similarly to what is observed in both the active and passive particle systems.
It is interesting to quantify the similarities of the fluctuation profiles for the active fluid (Fig. 2) and the contour model as a function of the contact angle (or equivalently in the contour model, ∆γ/γ).To achieve this, we use Monte Carlo to sample contours with a fixed area A. First, we calibrate the spatial discretisation of the contour model on the active model in the partially wet regime ε w = 10, by tuning the binsize δx of the contour model to match the bin location of the contact point in the active system.This gives δx = 0.066 √ A. We fix the discretisation, and by using a blurring scale σ blur = 2δx, we compute the average particle density in each bin.Then, we optimise the parameters of the contour model so that this density profile matches the ABP system at a given state point, by adjusting ∆γ/γ to minimise where ρABP j is the average of the local density ρ j defined as in ( 4), but now scaled to lie between 0 and 1; also ρ contour j is the corresponding (blurred) density for the contour model (which lies between 0 and 1 by definition).
Results are shown in Fig. 10.The contour model generates density profiles which closely match those of the ABPs not only for near-neutral conditions (e.g.10) but also close to the wetting transition ε w = 18.
The best fit profiles are associated with a contact angle cos(θ) = ∆γ/γ which confirms the wetting transition to occur around ε w = 18 (∆γ/γ = 0.95) and accords with direct measurements of θ via circular fits to the interfaces of Fig. 3.In generating each density profile, we take advantage of the low computational cost of the contour model by performing averages over 10 4 independent contours.We also show the local variance of the density.This is not part of the fitting procedure so it can be interpreted as a prediction of the contour model.As expected, the variance is large at the liquid-vapour interface.In the contour model, it also tends to be large where this interface meets the barrier.For the ABPs under partial wetting conditions, the same behaviour is observed, with large fluctuations near points of contact with the barrier.However, close to the transition point ε w = 18, the strongest density fluctuations in the ABPs occur away from these contact points.In fact, the contour model behaves analogously to the passive Lennard-Jones fluid (recall fig.6), with large fluctuations restricted to the contact points.This suggests that the geometrical setup of the active system may be influencing these density fluctuations, via the tips of the wedge, which are not explicitly accounted for either in the LJ system (which has walls terminating at reflective boundary conditions) or the contour model (which assumes infinite walls).Future work should clarify these potential finite-size effects and their coupling with the nonequilibrium fluxes of Fig. 4.
To analyse the density fluctuations in more detail, we characterise the behaviour of the local variance w(x, y) = Var[ρ], along the interface.To identify the interfacial region we take two density thresholds ρ lo , ρ hi and define a top-hat function Π ρ lo ,ρ hi (ρ) that is equal to unity for ρ lo < ρ < ρ hi (which is the interfacial region) and zero otherwise.Then we integrate this local variance over the interfacial region to obtain (the integral runs over the entire system and ρ is the local density).
To determine how much w varies as we move along the interface, we compute a corresponding spatial variance: Var[w] interface = (w(x, y)−w interface ) 2 Π ρ lo ,ρ hi [ρ(x, y)]dxdy.(13) We vary the barrier and the wall strength in the ABP model and the (passive) LJ system respectively, to obtain systems with various contact angles.As different models have different scales for the variance w = Var[ρ(x, y)], we rescale Var[w] interface of the active and passive model by their maximal values.Results are shown in Fig. 11, which also shows corresponding quantities for the contour model, computed as a function of ∆γ/γ.To reduce the noise of our estimates, we repeat the measurements over a number of independent runs, which depend on the complexity of the models (5 for the active system, 20 for the passive system and 100 for the contour model).
The contour model results can be accurately fitted with a polynomial expression (dashed line), and indicate an increasing spatial variance of the fluctuations, as the system moves from complete wetting [cos θ ≈ 1] towards the neutral point [cos θ ≈ 0].The particle-based simulations exhibit more pronounced statistical errors, but follow a comparable trend.
Overall, this figure illustrates that -for a given choice of the contact angle -and notwithstanding the very different specific interaction mechanisms with the wall/barrier, the three model systems display similar characteristic fluctuations of their interfaces, which match quantitatively (up to an overall scale factor).=More generally, Figs. 2, 6 and 10 demonstrate similar phenomenology for ABPs, passive LJ particles and the contour model, indicating that the principal features of the steady state density distribution in the wedge interior are well described by the minimal contour model, within which the most likely droplet shape is determined by a ratio of surface-tension-like quantities.This idea is developed further in the next Section.

VI. DISCUSSION AND CONCLUSIONS
We have considered active Brownian particles subject to a finite repulsive external potential as a simple model for surface phase transitions in active fluids.Our findings connect the behaviour of stable droplets to previously identified sharp spontaneous symmetry-breaking transitions from completely to partially wet states in threedimensional systems [22].In order to stabilise active droplets with nonzero contact angles, we utilise a wedgeshaped external potential.The droplets thus formed display distinctive density distributions as well as non-trivial steady currents.
A comparison with the wetting behaviour of a (Lennard-Jones) passive fluid confined by impenetrable, long-range attractive walls shows quantitative correspondences between the characteristic features of the density field and its fluctuations.Increasing the strength of the attractive wall-fluid interaction in the passive system plays an analogous role to increasing the potential barrier height in the active system, promoting a continuous change in curvature in the liquid-vapor interfaces and the promotion of density fluctuations in the vicinity of the contact line that increases as one moves from complete wetting (cos(θ) = 1) through the partial wetting regime to the neutral point (cos(θ) = 0).
Our results for the partial to complete wetting transition in ABPs provide new insights into phase separation phenomena in dry active matter.They suggest that, despite the inherent non-equilibrium mechanisms that engender phase separation and non-trivial flow patterns, the large-scale properties of the active model can be mapped onto an equilibrium one, an assertion that is supported by the finite-size scaling analysis of the density profiles in the partially wet regime.Such a scenario is similar to the hydrodynamic models proposed to rationalise motility-induced phase separation in the bulk [48] and connects to recent attempts to recover capillary-wave-like fluctuations in the nonequilibrium case [27,32,49].It will be interesting to study whether these findings generalize to more complex active matter systems with anisotropic interactions or longranged hydrodynamic effects: in such systems, MIPS may be modified [50] or suppressed [51], hence affecting dramatically one of the key ingredients for surface phase transitions.
The agreement of the contour model with the ABP system supports the correspondence between the interfacial behaviour of MIPS systems and equilibrium fluids.The parameters γlv , ∆γ w appearing in F are similar to (rescaled) surface tensions in equilibrium systems in that they set the probabilities of interfacial fluctuations.Indeed, one may use (8,9) to write the suggestive expression, γlv = − ∂ ∂ℓ lv ln P , (14) analogous to the definition of an equilibrium surface tension as a derivative of the free energy.Similarly ∆γ w = (∂/∂d wl ) ln P .However, it should not be assumed that γlv , ∆γ w are related to mechanical aspects of interfaces, in contrast to equilibrium systems where the surface tension controls anisotropic contributions to the stress tensor [52], as well as the magnitude of capillary fluctuations, and the Laplace pressure.It is a familiar feature of active matter that there is no a priori connection between fluctuations and mechanical forces, because standard fluctuation-dissipation theorems do not hold in these non-equilibrium steady states.In deriving Young's equation, the central (albeit trivial) assumption is that the observed large-scale droplet configuration is the one that maximises the probability, so its shape is naturally determined by probabilistic quantities.
To the extent that the contour model is an accurate description of interfacial fluctuations, the quantity γlv can be deduced from fluctuations of a planar MIPS interface: it is closely related to the capillary surface tensions discussed in [31,32].In equilibrium, ∆γ w is the difference between liquid-wall and vapour-wall surface tensions: these can be computed separately from simulations of the fluid in contact with the wall, under appropriate conditions.In this case, all the quantities that appear in Young's equation are known, and the contact angle θ can be predicted.In the active case, it is not clear how ∆γ w can be estimated without direct simulation of a wetting droplet, so the modified Young's equation (10) cannot predict the value of θ.Instead, one might infer ∆γ w by measuring the contact angle, under the assumption that (8,9) are suitable as a model for droplet-shape fluctuations.Under this assumption, the theory does make nontrivial predictions, for example that the contact angle of a fluid in a wedge should be independent of the aperture angle α, as found (at least approximately) in Fig. 1.Predictions of the most likely shape and its fluctuations would then also available for other wall geometries.However, the existence of steady-state currents in active fluids should be borne in mind since this could mediate nonlocal interaction terms in F, in which case (8,9) would not hold.More detailed tests of this theory would be desirable.
Our measurements of cos(θ) as a function of the wallfluid repulsive interaction strength (fig.3) show the wetting transition in ABPs to be first order.This finding accords with that of our previous work [22] in which the transition was studied using very different methods.Interestingly, however, in equilibrium fluids first order wetting is associated with the presence of long-range wallfluid interactions [20], while in our ABP simulations the repulsive barrier potential Eq. ( 2) is intrinsically very short ranged.It is possible that effective long-ranged interactions may be generated in active matter via manybody correlations and long-ranged flow fields; this is believed to occur at least for the case of a disordered boundary [53].Investigating the matter of what controls the order of the wetting transition in active wetting would be an interesting topic of further work.Other directions suggested by our research include a better understanding of the differences between three and two-dimensional active wetting, and the quantification of finite-size effects on the magnitude of the contact angles.cedures to fit the contact angles from stationary density profiles of the active and passive models: (1) a three-parameter fit with a circular profile; (2) a sevenparameter fit with a Fourier series representing the liquid-vapor interface.
For method (1), we identify the interface between vapor and liquid by thresholding the absolute value of the gradient of the density profile to extract the (x,y) coordinates of points exclusively at the vapor-liquid.We ignore interfacial points close to the walls as they worsen the quality of the result.We then perform a least-square circular fit with parameters R, x c , y x where R is the radius of the circle x c and y c are the coordinates of the centre of the circle.The centre of the circle corresponds to the tip of the wedge only when the contact angle is 90 • .In general, the wedge identifies the sector of a circle with the centre inside or outside of the wedge.When the centre is outside, the contact angle is measured as θ = π/2 − arcsin (|y c |/R), when it is inside, the angle is θ = π/2 + arcsin (|y c |/R).
For method (2), we use the parametrisation of Eq. 7 with M = 5 to fit the interface between the vapor and the liquid.In this case, we include all the points of the interface that are inside the wedge and evaluate the local slope of the interface at a fixed distance d θ = 3σ from the wall in order to approximate the contact angle.We take advantage of the closed form parameterisation (Eq.7) by taking its derivative analytically, which improves the numerical stability of the estimate of the contact angle.
The two methods yield comparable results in the partially wet regime.Method 2 breaks down on the approach to the wetting transition, although it continues to provide a lower bound for the contact angle at high ε w .This is illustrated in Fig. 12.We show that Young's equation can be obtained by minimising the free-energy-like quantity (9), subject to a constraint of fixed enclosed area.We emphasize however than in the active case F is not a thermodynamic free energy, instead it describes the log-probability of shape fluctuations, such that its minimum is the most likely shape.For convenience we consider a planar wall instead of a wedge geometry and define f (x) as the normal distance from the wall to the interface, see Fig. 13.We explain below how results for other geometries can be deduced from this computation.(In order for f to be single-valued our derivation also requires that θ ≤ π/2, this assumption can be also be relaxed.)We assume that the interface meets the wall at x = ±X so f (X) = 0 = f (−X): here X is a variational parameter.The length of the liquidvapour interface is ℓ lv = X −X 1 + f ′2 dx and the length of the liquid-wall interface is 2X.The enclosed area is X −X f dx.The problem can be simplified by a symmetry argument: the optimal droplet shape is symmetric about x = 0 so f ′ (0).We enforce the constraints of enclosed area A 0 and f (X) = 0 by Lagrange multipliers λ, µ.Hence it remains to extremise L[f, X] = 2 X 0 γlv 1 + f ′2 − λf dx−2∆γ w X+λA 0 −µf (X) (D1) (Note that this is a functional of f , also X is a scalar variational parameter and A 0 is the imposed area.) Within the calculus of variations we write f = f + δf and X = X + δX.Substituting this into L, the optimal droplet shape is identified by setting the first variation δL = 0. We find Integration by parts and using f (X) = 0 and f ′ (0) = 0 t e x i t s h a 1 _ b a s e 6 4 = " 4 e c b g / W 4 8 v h o 2 E 2 o 4 0 u Y F G L J u n M = " > A A A C A n i c d V D J S g N B E O 1 x j X G L e h S k M Q g e Z J i 4 J Z 4 M e P G Y g F k g G U J P p y f T p m e h u 0 Y I Q 2 5 6 9 6 q f 4 C 1 4 8 e C P + A U e 9 C P s J A r G 5 U H B 4 7 0 q 3 a F 7 4 9 p 4 M A b G 4 7 h 1 y v i c 2 U A T M J 4 / A L Q r n S M = < / l a t e x i t > ✓ < l a t e x i t s h a 1 _ b a s e 6 4 = " Y p X S I O D U f x w W h r l K 4 h v H a 2 j Z p 2 o = " > A A A C B H i c d V D L S s N A F J 3 4 r P F V d e l m s A i u Q t K n X Y h F N y 4 r 2 A e 0 o U y m k 3 b o T B J m J k I J 3 f o B b n X n 1 p 2 4 d e s n i F / g b z h t F K z o g Q u H c + 7 l 3 H u 9 i F G p b P v d W F h c W l 5 Z z a y Z 6 x u b

<
l a t e x i t s h a 1 _ b a s e 6 4 = " Y p X S I O D U f x w W h r l K 4 h v H a 2 j Z p 2 o = " > A A A C B H i c d V D L S s N A F J 3 4 r P F V d e l m s A i u Q t K n X Y h F N y 4 r 2 A e 0 o U y m k 3 b o T B J m J k I J 3 f o B b n X n 1 p 2 4 d e s n i F / g b z h

10 <
8 b X z B 6 Y g / H y C e 2 p n H w = < / l a t e x i t > 20 l a t e x i t s h a 1 _ b a s e 6 4 = " k h / l a t e x i t > t e x i t s h a 1 _ b a s e 6 4 = " Y p X S I O D U f x w W h r l K 4 h v H a 2 j Z p 2 o = " > A A A C B H i c d V D L S s N A F J 3 4 r P F V d e l m s A i u Q t K n X Y h F N y 4 r 2 A e 0 o U y m k 3 b o T B J m J k I J 3 f o B b n X n 1 p 2 4 d e s n i F / g b z h t F K z o g Q u H c + 7 l 3 H u 9 i F G p b P v d W F h c W l 5 Z z a y Z 6 x u b W 9 v Z n d 2 m D

10 <FIG. 2 .
FIG. 2. Average and variance of the local density in the active system.Upper row: stationary density profiles for phaseseparated active Brownian particles in contact with a finite wedge-shaped potential of repulsive strength εw for increasing values of the barrier strength in the range εw = 10 − 40.The system has periodic boundary conditions.Lower row: stationary profiles of local density fluctuation Var[ρ] for phase-separated active Brownian particles at increasing values of the barrier strength εw.

18 12 <
l a t e x i t s h a 1 _ b a s e 6 4 = " c A u 0 c z

FIG. 4 .
FIG. 4. Polarisation and flow.(a-d) Stationary average of the local particle orientation ϕ for increasing values of the barrier strength, accumulated on a grid of bin size ℓ = 2σ.The color coding indicates the angle of the local orientation; the arrow length corresponds to the average value of the polarization at the grid point, and it is scaled for clarity (leading to some crossings where the polarization is largest).For clarity, the plots depict only the immediate vicinity of the droplet, with grid binning δx = 2σ.Both the dilute phase and the interior of the droplet have negligible average polarization: only the interfaces (between the liquid and the vapor and between the barrier and the fluid phases) display significant contributions orthogonal to the interfaces.(e-h) Streamlines of the particle flux obtained as described in the text for the indicated values of the repulsive barrier strength εw overlayed to the density field.The line thickness scales linearly with the norm of the two-dimensional flux, with the same scaling used for all plots.The streamplots visualise the entire x-y span of the simulations.

Fig. 4 (
Fig.4(e-h) illustrates the flow patterns occurring at weak (ε w = 12) intermediate (ε w = 18) and strong (ε w = 30) and very strong (ε w = 40) barriers.In all cases, the wedge geometry shapes the structure of the flow field.Consider first the results for the strongest repulsive barrier (ε w = 40).Here the end-points of the wedge correspond to regions of opposing counterflows, reminiscent of the quadrupolar structure of hard objects mentioned above; the flow in the interior of the wedge is from the vapor into the liquid, consistent with the orientation field; particle flow through the wedge occurs only at the wedge corner, where two other counter-flows are formed, and where the flow field rotates in order to be orthogonal to the exterior of the wedge.

FIG. 5 .
FIG.5.Thickness of the 'lobes' ie the adsorbed layer on the exterior of the wedge.The upper panel shows the variation with the barrier strength for a system of size Lx = 100σ.The lower panel shows the scaling with system size for different values of the barrier strength εw both in the complete wet regime (εw = 30) and partially wet regime (εw = 12).
t e x i t s h a 1 _ b a s e 6 4 = " Y p X S I O D U f x w W h r l K 4 h v H a 2 j Z p 2 o = " > A A A C B H i c d V D L S s N A F J 3 4 r P F V d e l m s A i u Q t K n X Y h F N y 4 r 2 A e 0 o U y m k 3 b o T B J m J k I J 3 f o B b n X n 1 p 2 4 d e s n i F / g b z h t F K z o g Q u H c + 7 l 3 H u 9 i F G p b P v d W F h c W l 5 Z z a y Z 6 x u b W 9 v Z n d 2 m D s y g C u H 0 e N i n g m D F x p o g L K j e F e I h E g g r / a K 5 F M w 9 Q Q d D N d G f + T 4 f / k + a e c s p W 6 V L O 1 c 7 A y k y Y B 8 c g C P g g A q o g Q t Q B w 2 A A Q O 3 4 A 7 c G z f G o / F k P K e t C 8 b X z B 6 Y g / H y C e 2 p n H w = < / l a t e x i t > t e x i t s h a 1 _ b a s e 6 4 = " Y p X S I O D U f x w W h r l K 4 h v H a 2 j Z p 2 o = " > A A A C B H i c d V D L S s N A F J 3 4 r P F V d e l m s A i u Q t K n X Y h F N y 4 r 2 A e 0 o U y m k 3 b o T B J m J k I J 3 f o B b n X n 1 p 2 4 d e s n i F / g b z h t F K z o g Q u H c + 7 l 3 H u 9 i F G p b P v d W F h c W l 5 Z z a y Z 6 x u b s y g C u H 0 e N i n g m D F x p o g L K j e F e I h E g g r / a K 5 F M w 9 Q Q d D N d G f + T 4 f / k + a e c s p W 6 V L O 1 c 7 A y k y Y B 8 c g C P g g A q o g Q t Q B w 2 A A Q O 3 4 A 7 c G z f G o / F k P K e t C 8 b X z B 6 Y g / H y C e 2 p n H w = < / l a t e x i t > 20 < l a t e x i t s h a 1 _ b a s e 6 4 = " w EF C B q X T l p P z H E K l W 1 + H Y / G r d W c = " > A A A C E n i c d V D L S g M x F M 3 4 r P U 1 P n Z u Q o v Q 1 T A j 9 u F C K L o R c V H B P q A d S y Z N 2 9 D M g y R T q c P 8 h T s 3 b v U T 3 I k r w R 8 Q 3 f s b p q 2 C 9 X E g c O 4 5 9 3 J v j h M w K q R p v m h T 0 z O z c / O J h e T i 0 v L K q r 6 2 X h F + y D E p Y 5 / 5 v O Y g Q R j 1 S F l S y U g t 4 A S 5 D i N V p 3 c 4 9 K t 9 w g X 1 v T M 5 C I j t o o 5 H 2 x Q j q a S m v t n o I 0 4 C Q Z m q L s 6 j k + M Y 7 j f 1 t G n s F b K 7 2 R z 8 T S z D H C F d T G X e X v N P 1 6 W m / t 5 o + T h 0 i S c x Q 0 L U L T O Q d o S 4 p J i R O N k I B Q k Q 7 q E O q S v q I Z c I O x p d H 8 N t p b R g 2 + f q e R K O 1 O 8 T E X K F G L i O 6 n SR 7 I q f 3 l D 8 y 6 u H s l 2 w I + o F o S Q e H i 9 q h w x K H w 6 j g C 3 K C Z Z s o A j C n K p b I e 4 i j r B U g U 1 s w a 7 D a a c r Y 5 X M 1 / f h / 6 S y Y 1 g 5 I 3 u q I j o A Y y T A F k i B D L B A H h T B E S i B M s D g E t y A W 3 C n X W n 3 2 o P 2 O G 6 d 0 j 5 n N s A E t O c P P a y i v Q = = < / l a t e x i t >

FIG. 6 .FIG. 7 .
FIG.6.Average and variance of the density profile of the equilibrium system.Upper row: density profile for a droplet of Lennard-Jones liquid in coexistence with the vapour in contact with impenetrable, attractive walls of depth εw on the bottom and left sides, and reflective periodic boundary conditions on the top and right sides.Lower row: stationary local density fluctuation profiles for a droplet of Lennard-Jones liquid in coexistence with the vapour in contact with an impenetrable, attractive wall of depth ε LJ w .

1 < 2 <
H 6 P T M p J P Q F 6 w / U W H / m + 3 z 4 P 2 k c W L Z r O Z e o X D 0 D U x T B D t g F + 8 A G F V A F F 6 A G 6 o C A A b g F d + D e u D E e j S f j e d p a M L 5 m t s E M j J d P O + u b B A = = < / l a t e x i t > <l a t e x i t s h a 1 _ b a s e 6 4 = " t Q y 0 o I n q 8 B H c P Q L O s t O y d k R o / P o = " > A A A B / 3 i c d V D L S g M x F M 3 4 q L V a r X b p J l g E F z J k a j v W X c G N y 4 r 2 A e 1 Q M p l M G 5 p 5 k G S E Y e j C D 3 C r n + B O 3 P o p f o F f I Z i 2 C l b 0 w I X D O f d y 7 z 1 u z J l U C L 0 Z K 6 t r 6 7 m N / G Z h a 7 u 4 s 1 v a 2 + / I K B G E t k n E I 9 F z s a S c h b S t m O K 0 F w u K A 5 f T r j u 5 m P n d W y o k i 8 I b l c bU C f A o Z D 4 j W G n p 2 h t a w 1 I F m e d 2 w 0 Y W R C a y b K u B Z u S 0 V q 9 V o W W i O S p N Y 1 D + K O b S 1 r D 0 P v A i k g Q 0 V I R j K f s W i p W T Y a E Y 4 X R a G C S S x p h M 8 I j 2 N Q 1 x Q K W T z U + d w i O t e N C P h K 5 Q w b n 6 c y L D g Z R p 4 O r O A K u x / O 3 N x L + 8 f q L 8 h p O x M E 4 U D c l i k Z 9 w q C I 4 + x t 6 T F C i e K o J J o L p W y E Z Y 4 G J 0 u k s b S G B K 9 h o r K Y 6 m e / 3 4 f + k U z U t 2 6 x f 6 Y h O w A J 5 c A A O w T G w w B l o g k v Q A m 1 A w A j cg w f w a N w Z T 8 a z 8 b J o X T G + Z s p g C c b r J 4 P R m T A = < / l a t e x i t > d l a t e x i t s h a 1 _ b a s e 6 4 = " a P Z j 6 J 7 q F I m L P 7 b q s+ Y N A m v F 1 K 8 = " > A A A B / 3 i c d V D L S s N A F J 3 4 q L V a r X b p Z r A I L i R M E l q 7 L L h x W d E + o A 1 l M p m k Q y c P Z i Z C K F 3 4 A W 7 1 E 9 y J W z / F L / A r B K e t g h U 9 c O F w z r 3 c e 4 + X c i Y V Q m / G 2 v r G Z m G r u F 3 a 2 S 3 v 7 V c O D r s y y Q S h H Z L w R P Q 9 L C l n M e 0 o p j j t p 4 L i y O O 0 5 0 0 u 5 n 7 v l g r J k v h G 5 S l 1 I x z G L G A E K y 1 d + y N 7 V K k h 0 2 n a l l O H y L S R 4 y x I H a G G Y 0 H L R A v U W s a w + l E u 5 O 1 R 5 X 3 o J y S L a K w I x 1 I O L J Q q d 4 q F Y o T T W W m Y S Z p i M s E h H W g a 4 4 h K d 7 o 4 d Q Z P t O L D I B G 6 Y g U X 6 s + J K Y 6 k z C N P d 0 Z Y j e V v b y 7 + 5 Q 0 y F T T d K Y v T T N G Y L B c F G Y c q g f O / o c 8 E J Y r n m m A i m L 4 V k j E W m C i d z s o W E n m C h W M1 0 8 l 8 v w / / J 1 3 b t B p m / U p H d A a W K I I j c A x O g Q X O Q Q t c g j b o A A J C c A 8 e w K N x Z z w Z z 8 b L s n X N + J q p g h U Y r 5 9 0 V Z k m < / l a t e x i t > d l a t e x i t s h a 1 _ b a s e 6 4 = " f 2 G f 4 b J W

FIG. 8 .
FIG. 8. (a) Example paths generated by Monte-Carlo according to the free energy in Eq. 9.The quantities ℓ, d1, d2 are illustrated for a specific contour.(b) Average and (c) variance of the density determined by the statistics of the contours generated at fixed ∆γ/γ = 0.7.For these profiles, no Gaussian convolution is performed.
ε w = ABPs contour model < l a t e x i t s h a 1 _ b a s e 6 4 = " Q y e P e r b I o X q K a U Y 3 Y M U P b j M 1 + S I = " > A A A C G H i c b V D L S g M x F M 3 U q r W + q u 5 0 E 6 y C q z o j + N g I B V 2 4 r G A f 0 J Z y J 8 2 0 o c n M k N w R S i n 4 H X 6 A W / 0 E d + L W n V / g b 5 g + F r b 1 Q H I P 5 9 z L T Y 4 f S 2 H Q d b + d 1 F J 6 e W U 1 s 5 Z d 3 9 j c 2 s 7 t 7 F Z M l G j G y y y S k a 7 5 Y L g U I S + j Q M l r s e a g f M m r f u 9 m 5 F c f u T Y i C h + w H / O m g k 4 o A s E A r d T K 7 T d u u U S g j Q 4 o B a e T Q q + p W 3 B b u b y 9 x 6 C L x J u S f D E d j F F q 5 X 4 a 7 Y g l i o f I J B h T 9 9 w Y m w P Q K J j k w 2 w j M T w G 1 o M O r 1 s a g u K m O R j / Y U i P r d K m Q a T t C Z G O 1 b 8 T A 1 D G 9 J V v O x V g 1 8 x 7 I / E / r 5 5 g c N U c i D B O k I d s s i h I J M W I j g K h b a E 5 Q 9 m 3 B J g W 9 q 2 U d U E D Q x v b z B a m f C 0 6 X R z a Z L z 5 H B Z J 5 a z g X R T O 7 2 1 E R 2 S C D D k g h + S E e O S S F M k d K Z E y Y e S J v J B X 8 u Y 8 O + / O h / M 5 a U 0 5 0 5 k 9 M g P n 6 x f X 5 K H Y < / l a t e x i t > / = 0.0 < l a t e x i t s h a 1 _ b a s e 6 4 = " f 4 w 9 y 3 5 8 z C d D 2 q 6 P V p 1 W W a t k k n M = " > A A A C G X i c b V D L S g N B E J z V q D G + V j 3 m M h g F T 3 F X i I + D E N C D x w j m A U k I v Z P Z Z H B m d 5 n p F U L I w e / w A 7 z q J 3 g T r 5 7 8 A n / D y e N g 1 I K m i 6 p u e q a C R A q D n v f p L C x m l p Z X s q u 5 t f W N z S 1 3 e 6 d m 4 l Q z X m W x j H U j A M O l i H g V B U r e S D Q H F U h e D + 4 u x 3 7 9 n m s j 4 u g W B w l v K + h F I h Q M 0 E o d N 9 + 6 4 h K B t n q g F B x N G 7 2 g X v G 8 1 H E L X t G b g P 4 l / o w U y p l w g k r H / W p 1 Y 5 Y q H i G T Y E z T 9 x J s D 0 G j Y J K P c q 3 U 8 A T Y H f R 4 0 9 I I F D f t 4 e Q T I 3 p g l S 4 N Y 2 0 r Q j p R f 2 4 M Q R k z U I G d V I B 9 8 9 s b i / 9 5 z R T D s / Z Q R E m K P G L T Q 2 E q K c Z 0 n A j t C s 0 Z y o E l w L S w b 6 W s D x o Y 2 t z m r j A V a N H r 4 8 g m 4 / / O 4 S + p H R f 9 k 2 L p x k a 0 T 6 b I k j z Z I 4 f E J 6 e k T K 5 J h V Q J I w / k i T y T F + f R e X X e n P f p 6 I I z 2 9 k l c 3 A + v g F o s q I g < / l a t e x i t > / = 0.95 < l a t e x i t s h a 1 _ b a s e 6 4 = " E d I 3 e h f p L i / z b 9 H N l k / N 8 5 / C W 1 4 6 J f K p 7 e 2 I j 2 y R R Z k i d 7 5 J D 4 5 I y U y T W p k C p h 5 I E 8 k W f y 4 j w 6 r 8 6 b 8 z 4 d X X B m O 7 t k D s 7 H N 2 J L o h w = < / l a t e x i t > / = 0.64 < l a t e x i t s h a 1 _ b a s e 6 4 = " 2 s H q P s p 6 o I G h z W 3 u C l O + F t 0 e j m w y 3 m I O y 6 R 6 X v A u C x f 3 N q J j M k W a 5 M g R O S U e u S I l c k f K p E I Y e S I v 5 J W 8 O c / O u / P h f E 5 H V 5 z Z z i G Z g / P 1 C 2 C t o h s = < / l a t e x i t > / = 0.18 < l a t e x i t s h a 1 _ b a s e 6 4 = " t X 3 T V a E y d e d U Y 6 t l L L H e L c P r F X c = " > A A A C D 3 i c b V C 7 S g N B F J 2 N r x h f G y 0 s b A a j Y B V 2 B R + N G E h j G c E 8 I F n C 7 G Q 2 G T I z u 8 z M J o Q l H + E H 2 O o n p B N b w d L G L x D / w N J J Y m E S D 1 w 4 9 5 x 7 u Z f j R 4 w q 7 T g f V m p p e W V 1 L b 2 e 2 d j c 2 t 6 x s 7 s V F c Y S k z I O W S h r P l K E U U H K m m p G a p E k i P u M V P 1 u c e x X e 0 Q q G o o 7 P Y i I x 1 F b 0 I B i p I 3 U t L O N H p I k U p S Z r g + v X K d p 5 5 y 8 M w F c J O 4 v y V 1 / v 4 3 e 9 4 t f p a b 9 2 r A f T P J k e 6 o e W 8 s / u f V Y x 1 c e g k V U a y J w N N D Q c y g D u E 4 B 9 i i k m D N B o Y g L K n 5 F e I O k g h r k 9 b M F c x 9 S d s d P T T J u P M 5 L J L K a d 4 9 z 5 / d O r n C E Z g i D Q 7 A I T g B L r g A B X A D S q A M M O i D B / A I n q x 7 a 2 Q 9 W y / T 0 Z T 1 u 7 M H Z m C 9 / g D t 7 6 F 8 < / l a t e x i t > " w = 10 < l a t e x i t s h a 1 _ b a s e 6 4 = " M w f W F y p Y K N t v p I w u 7 T p y D o X F 2 d w= " > A A A C D 3 i c b V C 7 S g N B F J 2 N r x h f G y 0 s b A a j Y B V 2 A z 4 a M Z D G M o J 5 Q L K E 2 c k k G T I z u 8 z M J o Q l H + E H 2 O o n p B N b w d L G L x D / w N L J o z C J B y 6 c e 8 6 9 3 M v x Q 0 a V d p x P K 7 G y u r a + k d x M b W 3 v 7 O 7 Z 6 f 2 y C i K J S Q k H L J B V H y n C q C A l T T U j 1 V A S xH 1 G K n 6 3 M P Y r P S I V D c S 9 H o T E 4 6 g t a I t i p I 3 U s N P 1 H p I k V J S Z r g + v 3 V z D z j h Z Z w K 4 T N w Z y d z 8 v I 8 + D g v f x Y b 9 V W 8 G O O J E a M y Q U j X X C b U X I 6 k p Z m S Y q k e K h A h 3 U Z v U D B W I E + X F k 9 e H 8 N Q o T d g K p C m h 4 U T 9 u x E j r t S A + 2 a S I 9 1 R i 9 5 Y / M + r R b p 1 5 c V U h J E m A k 8 P t S I G d Q D H O c A m l Q R r N j A E Y U n N r x B 3 k E R Y m 7 T m r m D u S 9 r u 6 K F J x l 3 M Y Z m U c 1 n 3 I n t + 5 2 T y J 2 C K J D g C x + A M u O A S 5 M E t K I I S w K A P H s E T e L Y e r J H 1 Y r 1 O R x P W b O c A z M F 6 + w X x I a F + < / l a t e x i t > " w = 12 < l a t e x i t s h a 1 _ b a s e 6 4 = " 7 T / W e P F 3 h y P y t h d L 1 Q 9 v 6 e R O O S M = " > A A A C D 3 i c b V C 7 S g N B F J 2 N r x h f G y 0 s b A a j Y B V 2 x V c j B t J Y R j A P S J Y w O 5 l N h s z M L j O z C S H k I / w A W / 2 E d G I r W N r 4 B e I f W D p 5 F C b x w I V z z 7 m X e z l + x R H n 9 8 e s D e G y U B g x C a U p o O F b / b v Q R V 6 r H f T P J k W 6 p e W 8 k / u d V Y x 1 c e X 0 q o l g T g S e H g p h B H c J R D r B B J c G a 9 Q x B W F L z K 8 Q t J B H W J q 2 Z K 5 j 7 k j Z b e m C S c e d z W C S l 0 6 x 7 k T 2 / c z K 5 I z B B E h y A Q 3 A C X H A J c u A W F E A R Y N A F j + A J P F s P 1 t B 6 s V 4 n o w l r u r M H Z m C 9 / Q L 0 U 6 G A < / l a t e x i t > " w = 14 < l a t e x i t s h a 1 _ b a s e 6 4 = " 8 Y m r o + L H o Z 8 B P K P P U 1 x g f o z y e y Q = " > A A A C D 3 i c b V C 7 S g N B F J 2 N r x h f G y 0 s b A a j Y B V 2 B T W N G E h j G c E 8 I F n C 7 G S S D J m Z X W Z m E 8 K S j / A D b P U T 0 o m t Y G n j F 4 h / Y O n k U Z j E A x f O P e d e 7 u X 4 I a N K O 8 6 n 5 j 7 k r Y 7 e m i S c R d z W C b l 8 6 x 7 m b 2 4 c z L 5 E z B F E h y B Y 3 A G X H A F 8 u A W F E E J Y N A H j + A J P F s P 1 s h 6 s V 6 n o w l r t n M A 5 m C 9 / Q L 6 t 6 G E < / l a t e x i t > " w = 18 < l a t e x i t s h a 1 _ b a s e 6 4 = " t O X 7 N K 6 e j T r C T / C z r a s R G / W B v 4 0

FIG. 10 .
FIG. 10.Fitting the density profiles of the active system with the contour model for different values of the barrier strength εw.The upper row of each panel is the local density, while the lower row is the variance.Every panel displays two columns: the left column corresponds to the data from the region of the active systems contained inside the wedge, excluding the barrier region; the right column represents the fitted profiles produced by an ensemble of 10 000 contours.The fit is only performed on the density field, whereas the variance maps are inferred.The numerical values of the barrier strength of the active model εw and the fitted values of ∆γ/γ are also reported.

FIG. 11 .
FIG. 11.Scaled value of Var[w]interface where w = Var[ρ(x, y)] in the interfacial region between the vapor and the liquid, scaled by its maximum.For the contour model we consider, we convolve the density profiles with convolution scale σ blur = 2δx, where δx = 0.13R0.

FIG. 12 .
FIG.12.Comparison of the two independent measures of contact angles in the case of the active droplets.The hatched area indicates the detached regime.

FIG. 13 .
FIG. 13.(a) Sketch of geometry for derivation of Young's equation in Appendix D. (b) Corresponding geometry for the contour model, see Appendix E.