Drops on Polymer Brushes: Advances in Thin-Film Modeling of Adaptive Substrates

We briefly review recent advances in the hydrodynamic modeling of the dynamics of droplets on adaptive substrates, in particular, solids that are covered by polymer brushes. Thereby, the focus is on long-wave and full-curvature variants of mesoscopic hydrodynamic models in gradient dynamics form. After introducing the approach for films/drops of nonvolatile simple liquids on a rigid smooth solid substrate, it is first expanded to an arbitrary number of coupled degrees of freedom before considering the specific case of drops of volatile liquids on brush-covered solids. After presenting the model, its usage is illustrated by briefly considering the natural and forced spreading of drops of nonvolatile liquids on a horizontal brush-covered substrate, stick-slip motion of advancing contact lines as well as drops sliding down a brush-covered incline. Finally, volatile liquids are also considered.

Binding potential and wetting behavior of binary liquid mixtures on surfaces

We present a theory for the interfacial wetting phase behavior of binary liquid mixtures on rigid solid substrates, applicable to both miscible and immiscible mixtures. In particular, we calculate the binding potential as a function of the adsorptions, i.e., the excess amounts of each of the two liquids at the substrate. The binding potential fully describes the corresponding interfacial thermodynamics. Our approach is based on classical density functional theory. Binary liquid mixtures can exhibit complex bulk phase behavior, including both liquid-liquid and vapor-liquid phase separation, depending on the nature of the interactions among all the particles of the two different liquids, the temperature, and the chemical potentials. Here we show that the interplay between the bulk phase behavior of the mixture and the properties of the interactions with the substrate gives rise to a wide variety of interfacial phase behaviors, including mixing and demixing situations. We find situations where the final state is a coexistence of up to three different phases. We determine how the liquid density profiles close to the substrate change as the interaction parameters are varied and how these determine the form of the binding potential, which in certain cases can be a multivalued function of the adsorptions. We also present profiles for sessile droplets of both miscible and immiscible binary liquids.

Nonreciprocal Cahn-Hilliard Model Emerges as a Universal Amplitude Equation

Oscillatory behavior is ubiquitous in out-of-equilibrium systems showing spatiotemporal pattern formation. Starting from a linear large-scale oscillatory instability-a conserved-Hopf instability-that naturally occurs in many active systems with two conservation laws, we derive a corresponding amplitude equation. It belongs to a hierarchy of such universal equations for the eight types of instabilities in homogeneous isotropic systems resulting from the combination of three features: large-scale vs small-scale instability, stationary vs oscillatory instability, and instability without and with conservation law(s). The derived universal equation generalizes a phenomenological model of considerable recent interest, namely, the nonreciprocal Cahn-Hilliard model, and may be of a similar relevance for the classification of pattern forming systems as the complex Ginzburg-Landau equation.

Stationary broken parity states in active matter models

We demonstrate that several nonvariational continuum models commonly used to describe active matter as well as other active systems exhibit nongeneric behavior: each model supports asymmetric but stationary localized states even in the absence of pinning at heterogeneities. Moreover, such states only begin to drift following a drift-transcritical bifurcation as the activity increases. Asymmetric stationary states should only exist in variational systems, i.e., in models with gradient structure. In other words, such states are expected in passive systems, but not in active systems where the gradient structure of the model is broken by activity. We identify a "spurious" gradient dynamics structure of these models that is responsible for this nongeneric behavior, and determine the types of additional terms that render the models generic, i.e., with asymmetric states that appear via drift-pitchfork bifurcations and are generically moving. We provide detailed illustrations of our results using numerical continuation of resting and steadily drifting states in both generic and nongeneric cases.

Stick-slip dynamics in the forced wetting of polymer brushes

We study the static and dynamic wetting of adaptive substrates using a mesoscopic hydrodynamic model for a liquid droplet on a solid substrate covered by a polymer brush. First, we show that on the macroscale Young's law still holds for the equilibrium contact angle and that on the mesoscale a Neumann-type law governs the shape of the wetting ridge. Following an analytic and numeric assessment of the static profiles of droplet and wetting ridge, we examine the dynamics of the wetting ridge for a liquid meniscus that is advanced at constant mean speed. In other words, we consider an inverse Landau-Levich case where a brush-covered plate is introduced into (and not drawn from) a liquid bath. We find a characteristic stick-slip motion that emerges when the dynamic contact angle of the stationary moving meniscus decreases with increasing velocity, and relate the onset of slip to Gibbs' inequality and to a cross-over in relevant time scales.

Nonequilibrium configurations of swelling polymer brush layers induced by spreading drops of weakly volatile oil

Polymer brush layers are responsive materials that swell in contact with good solvents and their vapors. We deposit drops of an almost completely wetting volatile oil onto an oleophilic polymer brush layer and follow the response of the system upon simultaneous exposure to both liquid and vapor. Interferometric imaging shows that a halo of partly swollen polymer brush layer forms ahead of the moving contact line. The swelling dynamics of this halo is controlled by a subtle balance of direct imbibition from the drop into the brush layer and vapor phase transport and can lead to very long-lived transient swelling profiles as well as nonequilibrium configurations involving thickness gradients in a stationary state. A gradient dynamics model based on a free energy functional with three coupled fields is developed and numerically solved. It describes experimental observations and reveals how local evaporation and condensation conspire to stabilize the inhomogeneous nonequilibrium stationary swelling profiles. A quantitative comparison of experiments and calculations provides access to the solvent diffusion coefficient within the brush layer. Overall, the results highlight the-presumably generally applicable-crucial role of vapor phase transport in dynamic wetting phenomena involving volatile liquids on swelling functional surfaces.

Non-reciprocity induces resonances in a two-field Cahn-Hilliard model

We consider a non-reciprocally coupled two-field Cahn-Hilliard system that has been shown to allow for oscillatory behaviour and suppression of coarsening. After introducing the model, we first review the linear stability of steady uniform states and show that all instability thresholds are identical to the ones for a corresponding two-species reaction-diffusion system. Next, we consider a specific interaction of linear modes-a 'Hopf-Turing' resonance-and derive the corresponding amplitude equations using a weakly nonlinear approach. We discuss the weakly nonlinear results and finally compare them with fully nonlinear simulations for a specific conserved amended FitzHugh-Nagumo system. We conclude with a discussion of the limitations of the employed weakly nonlinear approach. This article is part of the theme issue 'New trends in pattern formation and nonlinear dynamics of extended systems'.

From a microscopic inertial active matter model to the Schrödinger equation

Active field theories, such as the paradigmatic model known as 'active model B+', are simple yet very powerful tools for describing phenomena such as motility-induced phase separation. No comparable theory has been derived yet for the underdamped case. In this work, we introduce active model I+, an extension of active model B+ to particles with inertia. The governing equations of active model I+ are systematically derived from the microscopic Langevin equations. We show that, for underdamped active particles, thermodynamic and mechanical definitions of the velocity field no longer coincide and that the density-dependent swimming speed plays the role of an effective viscosity. Moreover, active model I+ contains an analog of the Schrödinger equation in Madelung form as a limiting case, allowing one to find analoga of the quantum-mechanical tunnel effect and of fuzzy dark matter in active fluids. We investigate the active tunnel effect analytically and via numerical continuation.

Exploring bifurcations in Bose-Einstein condensates via phase field crystal models

To facilitate the analysis of pattern formation and the related phase transitions in Bose-Einstein condensates, we present an explicit approximate mapping from the nonlocal Gross-Pitaevskii equation with cubic nonlinearity to a phase field crystal (PFC) model. This approximation is valid close to the superfluid-supersolid phase transition boundary. The simplified PFC model permits the exploration of bifurcations and phase transitions via numerical path continuation employing standard software. While revealing the detailed structure of the bifurcations present in the system, we demonstrate the existence of localized states in the PFC approximation. Finally, we discuss how higher-order nonlinearities change the structure of the bifurcation diagram representing the transitions found in the system.

Symmetry-breaking, motion and bistability of active drops through polarization-surface coupling

Cell crawling crucially depends on the collective dynamics of the acto-myosin cytoskeleton. However, it remains an open question to what extent cell polarization and persistent motion depend on continuous regulatory mechanisms and autonomous physical mechanisms. Experiments on cell fragments and theoretical considerations for active polar liquids have highlighted that physical mechanisms induce motility through splay and bend configurations in a nematic director field. Here, we employ a simple model, derived from basic thermodynamic principles, for active polar free-surface droplets to identify a different mechanism of motility. Namely, active stresses drive drop motion through spatial variations of polarization strength. This robustly induces parity-symmetry breaking and motility even for liquid ridges (2D drops) and adds to splay- and bend-driven pumping in 3D geometries. Intriguingly, then, stable polar moving and axisymmetric resting states may coexist, reminiscent of the interconversion of moving and resting keratocytes by external stimuli. The identified additional motility mode originates from a competition between the elastic bulk energy and the polarity control exerted by the drop surface. As it already breaks parity-symmetry for passive drops, the resulting back-forth asymmetry enables active stresses to effectively pump liquid and drop motion ensues.

Soft wetting with (a)symmetric Shuttleworth effect

The wetting of soft polymer substrates brings in multiple complexities when compared with the wetting on rigid substrates. The contact angle of the liquid is no longer governed by Young’s Law, but is affected by the substrate’s bulk and surface deformations. On top of that, elastic interfaces exhibit a surface energy that depends on how much they are stretched—a feature known as the Shuttleworth effect (or as surface-elasticity). Here, we present two models through which we explore the wetting of drops in the presence of a strong Shuttleworth effect. The first model is macroscopic in character and consistently accounts for large deformations via a neo-Hookean elasticity. The second model is based on a mesoscopic description of wetting, using a reduced description of the substrate’s elasticity. While the second model is more empirical in terms of the elasticity, it enables a gradient dynamics formulation for soft wetting dynamics. We provide a detailed comparison between the equilibrium states predicted by the two models, from which we deduce robust features of soft wetting in the presence of a strong Shuttleworth effect. Specifically, we show that the (a)symmetry of the Shuttleworth effect between the ‘dry’ and ‘wet’ states governs horizontal deformations in the substrate. Our results are discussed in the light of recent experiments on the wettability of stretched substrates.

Gradient-dynamics model for liquid drops on elastic substrates

The wetting of soft elastic substrates exhibits many features that have no counterpart on rigid surfaces. Modelling the detailed elastocapillary interactions is challenging, and has so far been limited to single contact lines or single drops. Here we propose a reduced long-wave model that captures the main qualitative features of statics and dynamics of soft wetting, but which can be applied to ensembles of droplets. The model has the form of a gradient dynamics on an underlying free energy that reflects capillarity, wettability and compressional elasticity. With the model we first recover the double transition in the equilibrium contact angles that occurs when increasing substrate softness from ideally rigid towards very soft (i.e., liquid). Second, the spreading of single drops of partially and completely wetting liquids is considered showing that known dependencies of the dynamic contact angle on contact line velocity are well reproduced. Finally, we go beyond the single droplet picture and consider the coarsening for a two-drop system as well as for a large ensemble of drops. It is shown that the dominant coarsening mode changes with substrate softness in a nontrivial way.

Suppression of coarsening and emergence of oscillatory behavior in a Cahn-Hilliard model with nonvariational coupling

We investigate a generic two-field Cahn-Hilliard model with variational and nonvariational coupling. It describes, for instance, passive and active ternary mixtures, respectively. Already a linear stability analysis of the homogeneous mixed state shows that activity not only allows for the usual large-scale stationary (Cahn-Hilliard) instability of the well-known passive case but also for small-scale stationary (Turing) and large-scale oscillatory (Hopf) instabilities. In consequence of the Turing instability, activity may completely suppress the usual coarsening dynamics. In a fully nonlinear analysis, we first briefly discuss the passive case before focusing on the active case. Bifurcation diagrams and selected direct time simulations are presented that allow us to establish that nonvariational coupling (i) can partially or completely suppress coarsening and (ii) may lead to the emergence of drifting and oscillatory states. Throughout, we emphasize the relevance of conservation laws and related symmetries for the encountered intricate bifurcation behavior.

Two-dimensional localized states in an active phase-field-crystal model

The active phase-field-crystal (active PFC) model provides a simple microscopic mean field description of crystallization in active systems. It combines the PFC model (or conserved Swift-Hohenberg equation) of colloidal crystallization and aspects of the Toner-Tu theory for self-propelled particles. We employ the active PFC model to study the occurrence of localized and periodic active crystals in two spatial dimensions. Due to the activity, crystalline states can undergo a drift instability and start to travel while keeping their spatial structure. Based on linear stability analyses, time simulations, and numerical continuation of the fully nonlinear states, we present a detailed analysis of the bifurcation structure of resting and traveling states. We explore, for instance, how the slanted homoclinic snaking of steady localized states found for the passive PFC model is modified by activity. Morphological phase diagrams showing the regions of existence of various solution types are presented merging the results from all the analysis tools employed. We also study how activity influences the crystal structure with transitions from hexagons to rhombic and stripe patterns. This in-depth analysis of a simple PFC model for active crystals and swarm formation provides a clear general understanding of the observed multistability and associated hysteresis effects, and identifies thresholds for qualitative changes in behavior.

Efficient calculation of phase coexistence and phase diagrams: application to a binary phase-field-crystal model

We show that one can employ well-established numerical continuation methods to efficiently calculate the phase diagram for thermodynamic systems described by a suitable free energy functional. In particular, this involves the determination of lines of phase coexistence related to first order phase transitions and the continuation of triple points. To illustrate the method we apply it to a binary phase-field-crystal model for the crystallisation of a mixture of two types of particles. The resulting phase diagram is determined for one- and two-dimensional domains. In the former case it is compared to the diagram obtained from a one-mode approximation. The various observed liquid and crystalline phases and their stable and metastable coexistence are discussed as well as the temperature-dependence of the phase diagrams. This includes the (dis)appearance of critical points and triple points. We also relate bifurcation diagrams for finite-size systems to the thermodynamics of phase transitions in the infinite-size limit.

Phase-field-crystal description of active crystallites: Elastic and inelastic collisions

The active Phase-Field-Crystal (aPFC) model combines elements of the Toner-Tu theory for self-propelled particles and the classical Phase-Field-Crystal (PFC) model that describes the transition between liquid and crystalline phases. In the liquid-crystal coexistence region of the PFC model, crystalline clusters exist in the form of localized states that coexist with a homogeneous background. At sufficiently strong activity (related to self-propulsion strength), they start to travel. We employ numerical path continuation and direct time simulations to first investigate the existence regions of different types of localized states in one spatial dimension. The results are summarized in morphological phase diagrams in the parameter plane spanned by activity and mean density. Then we focus on the interaction of traveling localized states, studying their collision behavior. As a result, we distinguish "elastic" and "inelastic" collisions. In the former, localized states recover their properties after a collision, while in the latter, they may completely or partially annihilate, forming resting bound states or various traveling states.

Adaptive stochastic continuation with a modified lifting procedure applied to complex systems

Many complex systems occurring in the natural or social sciences or economics are frequently described on a microscopic level, e.g., by lattice- or agent-based models. To analyze the states of such systems and their bifurcation structure on the level of macroscopic observables, one has to rely on equation-free methods like stochastic continuation. Here we investigate how to improve stochastic continuation techniques by adaptively choosing the parameters of the algorithm. This allows one to obtain bifurcation diagrams quite accurately, especially near bifurcation points. We introduce lifting techniques which generate microscopic states with a naturally grown structure, which can be crucial for a reliable evaluation of macroscopic quantities. We show how to calculate fixed points of fluctuating functions by employing suitable linear fits. This procedure offers a simple measure of the statistical error. We demonstrate these improvements by applying the approach in analyses of (i) the Ising model in two dimensions, (ii) an active Ising model, and (iii) a stochastic Swift-Hohenberg model. We conclude by discussing the abilities and remaining problems of the technique.

Thin-film modeling of resting and moving active droplets

We propose a generic model for thin films and shallow drops of a polar active liquid that have a free surface and are in contact with a solid substrate. The model couples evolution equations for the film height and the local polarization in the form of a gradient dynamics supplemented with active stresses and fluxes. A wetting energy for a partially wetting liquid is incorporated allowing for motion of the liquid-solid-gas contact line. This gives a consistent basis for the description of drops of dense bacterial suspensions or compact aggregates of living cells on solid substrates. As example, we analyze the dynamics of two-dimensional active drops (i.e., ridges) and demonstrate how active forces compete with passive surface forces to shape droplets and drive their motion. In our simple two-dimensional scenario we find that defect structures within the polarization profile drastically influence the shape and motility of active droplets. Thus, we can observe a transition from resting to motile droplets via the elimination of defects in the polarization profile. Furthermore, droplet motility is modulated by strong active stresses. Contractile stresses even lead to topological changes, i.e., drop splitting, which is naturally encoded in the evolution equations.

Bifurcations of front motion in passive and active Allen-Cahn-type equations

The well-known cubic Allen-Cahn (AC) equation is a simple gradient dynamics (or variational) model for a nonconserved order parameter field. After revising main literature results for the occurrence of different types of moving fronts, we employ path continuation to determine their bifurcation diagram in dependence of the external field strength or chemical potential. We then employ the same methodology to systematically analyze fronts for more involved AC-type models. In particular, we consider a cubic-quintic variational AC model and two different nonvariational generalizations. We determine and compare the bifurcation diagrams of front solutions in the four considered models.

Collective Cell Migration in Embryogenesis Follows the Laws of Wetting

Collective cell migration is a fundamental process during embryogenesis and its initial occurrence, called epiboly, is an excellent in vivo model to study the physical processes involved in collective cell movements that are key to understanding organ formation, cancer invasion, and wound healing. In zebrafish, epiboly starts with a cluster of cells at one pole of the spherical embryo. These cells are actively spreading in a continuous movement toward its other pole until they fully cover the yolk. Inspired by the physics of wetting, we determine the contact angle between the cells and the yolk during epiboly. By choosing a wetting approach, the relevant scale for this investigation is the tissue level, which is in contrast to other recent work. Similar to the case of a liquid drop on a surface, one observes three interfaces that carry mechanical tension. Assuming that interfacial force balance holds during the quasi-static spreading process, we employ the physics of wetting to predict the temporal change of the contact angle. Although the experimental values vary dramatically, the model allows us to rescale all measured contact-angle dynamics onto a single master curve explaining the collective cell movement. Thus, we describe the fundamental and complex developmental mechanism at the onset of embryogenesis by only three main parameters: the offset tension strength, α, that gives the strength of interfacial tension compared to other force-generating mechanisms; the tension ratio, δ, between the different interfaces; and the rate of tension variation, λ, which determines the timescale of the whole process.

Effects of time-periodic forcing in a Cahn-Hilliard model for Langmuir-Blodgett transfer

The influence of a temporal forcing on the pattern formation in Langmuir-Blodgett transfer is studied employing a generalized Cahn-Hilliard model. The occurring frequency-locking effects allow for controlling the pattern formation process. In the case of one-dimensional (i.e., stripe) patterns one finds various synchronization phenomena such as entrainment between the distance of deposited stripes and the forcing frequency. In two dimensions, the temporal forcing gives rise to the formation of intricate complex patterns such as vertical stripes, oblique stripes, and lattice structures. Remarkably, it is possible to influence the system in the spatial direction perpendicular to the forcing direction leading to synchronization in two spatial dimensions.

Connecting Monotonic and Oscillatory Motions of the Meniscus of a Volatile Polymer Solution to the Transport of Polymer Coils and Deposit Morphology

We study the deposition mechanisms of polymer from a confined  meniscus of volatile liquid. In particular, we investigate the physical processes that are responsible for qualitative changes in the pattern deposition of polymer and the underlying interplay of the state of pattern deposition, motion of the meniscus, and the transport of polymer within the meniscus. As a model system we evaporate a solution of poly(methyl methacrylate) (PMMA) in toluene. Different deposition patterns are observed when varying the molecular mass,  the initial concentration of the solute, and temperature; these  are systematically presented in the form of morphological phase diagrams. The modi of deposition and meniscus motion are correlated. They vary with the ratio between the evaporation-driven convective flux and the diffusive flux of the polymer coils in the solution. In the case of a diffusion-dominated solute transport, the solution monotonically dewets the solid substrate by evaporation, supporting continuous contact line motion and continuous polymer deposition. However, a convection-dominated transport results in an oscillatory ratcheting dewetting-wetting motion of the contact line with more pronounced dewetting phases. The deposition process is then periodic and produces a stripe pattern. The oscillatory motion of the meniscus differs from the well documented stick-slip motion of the meniscus, observed as well, and is attributed to the opposing influences of evaporation and Marangoni stresses, which alternately dominate the deposition process.

Resting and traveling localized states in an active phase-field-crystal model

The conserved Swift-Hohenberg equation (or phase-field-crystal [PFC] model) provides a simple microscopic description of the thermodynamic transition between fluid and crystalline states. Combining it with elements of the Toner-Tu theory for self-propelled particles, Menzel and Löwen [Phys. Rev. Lett. 110, 055702 (2013)PRLTAO0031-900710.1103/PhysRevLett.110.055702] obtained a model for crystallization (swarm formation) in active systems. Here, we study the occurrence of resting and traveling localized states, i.e., crystalline clusters, within the resulting active PFC model. Based on linear stability analyses and numerical continuation of the fully nonlinear states, we present a detailed analysis of the bifurcation structure of periodic and localized, resting and traveling states in a one-dimensional active PFC model. This allows us, for instance, to explore how the slanted homoclinic snaking of steady localized states found for the passive PFC model is amended by activity. A particular focus lies on the onset of motion, where we show that it occurs either through a drift-pitchfork or a drift-transcritical bifurcation. A corresponding general analytical criterion is derived.

Equilibrium Contact Angle and Adsorption Layer Properties with Surfactants

The three-phase contact line of a droplet on a smooth surface can be characterized by the Young equation. It relates the interfacial energies to the macroscopic contact angle θ_e. On the mesoscale, wettability is modeled by a film-height-dependent wetting energy f( h). Macro- and mesoscale descriptions are consistent if γ cos θ_e = γ + f( h_a), where γ and h_a are the liquid-gas interface energy and the thickness of the equilibrium liquid adsorption layer, respectively. Here, we derive a similar consistency condition for the case of a liquid covered by an insoluble surfactant. At equilibrium, the surfactant is spatially inhomogeneously distributed, implying a nontrivial dependence of θ_e on surfactant concentration. We derive macroscopic and mesoscopic descriptions of a contact line at equilibrium and show that they are consistent only if a particular dependence of the wetting energy on the surfactant concentration is imposed. This is illustrated by a simple example of dilute surfactants, for which we show excellent agreement between theory and time-dependent numerical simulations.

Modelling of surfactant-driven front instabilities in spreading bacterial colonies

The spreading of bacterial colonies at solid-air interfaces is determined by the physico-chemical properties of the involved interfaces. The production of surfactant molecules by bacteria is a widespread strategy that allows the colony to efficiently expand over the substrate. On the one hand, surfactant molecules lower the surface tension of the colony, effectively increasing the wettability of the substrate, which facilitates spreading. On the other hand, gradients in the surface concentration of surfactant molecules result in Marangoni flows that drive spreading. These flows may cause an instability of the circular colony shape and the subsequent formation of fingers. In this work, we study the effect of bacterial surfactant production and substrate wettability on colony growth and shape within the framework of a hydrodynamic thin film model. We show that variations in the wettability and surfactant production are sufficient to reproduce four different types of colony growth, which have been described in the literature, namely, arrested and continuous spreading of circular colonies, slightly modulated front lines and the formation of pronounced fingers.

Sliding Drops: Ensemble Statistics from Single Drop Bifurcations

Ensembles of interacting drops that slide down an inclined plate show a dramatically different coarsening behavior as compared to drops on a horizontal plate: As drops of different size slide at different velocities, frequent collisions result in fast coalescence. However, above a certain size individual sliding drops are unstable and break up into smaller drops. Therefore, the long-time dynamics of a large drop ensemble is governed by a balance of merging and splitting. We employ a long-wave film height evolution equation and determine the dynamics of the drop size distribution towards a stationary state from direct numerical simulations on large domains. The main features of the distribution are then related to the bifurcation diagram of individual drops obtained by numerical path continuation. The gained knowledge allows us to develop a Smoluchowski-type statistical model for the ensemble dynamics that well compares to full direct simulations.

Continuous versus Arrested Spreading of Biofilms at Solid-Gas Interfaces: The Role of Surface Forces

We introduce and analyze a model for osmotically spreading bacterial colonies at solid-air interfaces that includes wetting phenomena, i.e., surface forces. The model is based on a hydrodynamic description for liquid suspensions which is supplemented by bioactive processes. We show that surface forces determine whether a biofilm can expand laterally over a substrate and provide experimental evidence for the existence of a transition between continuous and arrested spreading for Bacillus subtilis biofilms. In the case of arrested spreading, the lateral expansion of the biofilm is confined, albeit the colony is biologically active. However, a small reduction in the surface tension of the biofilm is sufficient to induce spreading. The incorporation of surface forces into our hydrodynamic model allows us to capture this transition in biofilm spreading behavior.

Nudged elastic band calculation of the binding potential for liquids at interfaces

The wetting behavior of a liquid on solid substrates is governed by the nature of the effective interaction between the liquid-gas and the solid-liquid interfaces, which is described by the binding or wetting potential g(h) which is an excess free energy per unit area that depends on the liquid film height h. Given a microscopic theory for the liquid, to determine g(h), one must calculate the free energy for liquid films of any given value of h, i.e., one needs to create and analyze out-of-equilibrium states, since at equilibrium there is a unique value of h, specified by the temperature and chemical potential of the surrounding gas. Here we introduce a Nudged Elastic Band (NEB) approach to calculate g(h) and illustrate the method by applying it in conjunction with a microscopic lattice density functional theory for the liquid. We also show that the NEB results are identical to those obtained with an established method based on using a fictitious additional potential to stabilize the non-equilibrium states. The advantages of the NEB approach are discussed.

Films, layers, and droplets: The effect of near-wall fluid structure on spreading dynamics

We present a study of the spreading of liquid droplets on a solid substrate at very small scales. We focus on the regime where effective wetting energy (binding potential) and surface tension effects significantly influence steady and spreading droplets. In particular, we focus on strong packing and layering effects in the liquid near the substrate due to underlying density oscillations in the fluid caused by attractive substrate-liquid interactions. We show that such phenomena can be described by a thin-film (or long-wave or lubrication) model including an oscillatory Derjaguin (or disjoining or conjoining) pressure and explore the effects it has on steady droplet shapes and the spreading dynamics of droplets on both an adsorption (or precursor) layer and completely dry substrates. At the molecular scale, commonly used two-term binding potentials with a single preferred minimum controlling the adsorption layer height are inadequate to capture the rich behavior caused by the near-wall layered molecular packing. The adsorption layer is often submonolayer in thickness, i.e., the dynamics along the layer consists of single-particle hopping, leading to a diffusive dynamics, rather than the collective hydrodynamic motion implicit in standard thin-film models. We therefore modify the model in such a way that for thicker films the standard hydrodynamic theory is realized, but for very thin layers a diffusion equation is recovered.

Dip-coating with prestructured substrates: transfer of simple liquids and Langmuir-Blodgett monolayers

When a plate is withdrawn from a liquid bath, either a static meniscus forms in the transition region between the bath and the substrate or a liquid film of finite thickness (a Landau-Levich film) is transferred onto the moving substrate. If the substrate is inhomogeneous, e.g. has a prestructure consisting of stripes of different wettabilities, the meniscus can be deformed or show a complex dynamic behavior. Here we study the free surface shape and dynamics of a dragged meniscus occurring for striped prestructures with two orientations, parallel and perpendicular to the transfer direction. A thin film model is employed that accounts for capillarity through a Laplace pressure and for the spatially varying wettability through a Derjaguin (or disjoining) pressure. Numerical continuation is used to obtain steady free surface profiles and corresponding bifurcation diagrams in the case of substrates with different homogeneous wettabilities. Direct numerical simulations are employed in the case of the various striped prestructures. The final part illustrates the importance of our findings for particular applications that involve complex liquids by modeling a Langmuir-Blodgett transfer experiment. There, one transfers a monolayer of an insoluble surfactant that covers the surface of the bath onto the moving substrate. The resulting pattern formation phenomena can be crucially influenced by the hydrodynamics of the liquid meniscus that itself depends on the prestructure on the substrate. In particular, we show how prestructure stripes parallel to the transfer direction lead to the formation of bent stripes in the surfactant coverage after transfer and present similar experimental results.

Mode instabilities and dynamic patterns in a colony of self-propelled surfactant particles covering a thin liquid layer

We consider a colony of point-like self-propelled surfactant particles (swimmers) without direct interactions that cover a thin liquid layer on a solid support. The particles predominantly swim normal to the free film surface with only a small component parallel to the film surface. The coupled dynamics of the swimmer density and film height profile is captured in a long-wave model allowing for diffusive and convective transport of the swimmers (including rotational diffusion). The dynamics of the film height profile is determined by i) the upward pushing force of the swimmers onto the liquid-gas interface, ii) the solutal Marangoni force due to gradients in the swimmer concentration, and iii) the rotational diffusion of the swimmers together with the in-plane active motion. After reviewing and extending the analysis of the linear stability of the uniform state, we analyse the fully nonlinear dynamic equations and show that point-like swimmers, which only interact via long-wave deformations of the liquid film, self-organise in highly regular (standing, travelling, and modulated waves) and various irregular patterns.

Bimodality and re-entrant behaviour in the hierarchical self-assembly of polymeric nanoparticles

We show that a film of a suspension of polymer grafted nanoparticles on a liquid substrate can be employed to create two-dimensional nanostructures with a remarkable variation in the pattern length scales. The presented experiments also reveal the emergence of concentration-dependent bimodal patterns as well as re-entrant behaviour that involves length scales due to dewetting and compositional instabilities. The experimental observations are explained through a gradient dynamics model consisting of coupled evolution equations for the height of the suspension film and the concentration of polymer. Using a Flory-Huggins free energy functional for the polymer solution, we show in a linear stability analysis that the thin film undergoes dewetting and/or compositional instabilities depending on the concentration of the polymer in the solution. We argue that the formation via 'hierarchical self-assembly' of various functional nanostructures observed in different systems can be explained as resulting from such an interplay of instabilities.

Instabilities of Layers of Deposited Molecules on Chemically Stripe Patterned Substrates: Ridges versus Drops

A mesoscopic continuum model is employed to analyze the transport mechanisms and structure formation during the redistribution stage of deposition experiments where organic molecules are deposited on a solid substrate with periodic stripe-like wettability patterns. Transversally invariant ridges located on the more wettable stripes are identified as very important transient states and their linear stability is analyzed accompanied by direct numerical simulations of the fully nonlinear evolution equation for two-dimensional substrates. It is found that there exist two different instability modes that lead to different nonlinear evolutions that result (i) at large ridge volume in the formation of bulges that spill from the more wettable stripes onto the less wettable bare substrate and (ii) at small ridge volume in the formation of small droplets located on the more wettable stripes. In addition, the influence of different transport mechanisms during redistribution is investigated focusing on the cases of convective transport with no-slip at the substrate, transport via diffusion in the film bulk and via diffusion at the film surface. In particular, it is shown that the transport process does neither influence the linear stability thresholds nor the sequence of morphologies observed in the time simulation, but only the ratio of the time scales of the different process phases.

A Variational approach to thin film hydrodynamics of binary mixtures

In order to model the dynamics of thin films of mixtures, solutions, and suspensions, a thermodynamically consistent formulation is needed such that various coexisting dissipative processes with cross couplings can be correctly described in the presence of capillarity, wettability, and mixing effects. In the present work, we apply Onsager's variational principle to the formulation of thin film hydrodynamics for binary fluid mixtures. We first derive the dynamic equations in two spatial dimensions, one along the substrate and the other normal to the substrate. Then, using long-wave asymptotics, we derive the thin film equations in one spatial dimension along the substrate. This enables us to establish the connection between the present variational approach and the gradient dynamics formulation for thin films. It is shown that for the mobility matrix in the gradient dynamics description, Onsager's reciprocal symmetry is automatically preserved by the variational derivation. Furthermore, using local hydrodynamic variables, our variational approach is capable of introducing diffusive dissipation beyond the limit of dilute solute. Supplemented with a Flory-Huggins-type mixing free energy, our variational approach leads to a thin film model that treats solvent and solute in a symmetric manner. Our approach can be further generalized to include more complicated free energy and additional dissipative processes.

Nonequilibrium Gibbs' criterion for completely wetting volatile liquids

During the spreading of a liquid over a solid substrate, the contact line can stay pinned at sharp edges until the contact angle exceeds a critical value. At (or sufficiently near) equilibrium, this is known as Gibbs' criterion. Here, we show both experimentally and theoretically that, for completely wetting volatile liquids, there also exists a dynamically-produced contribution to the critical angle for depinning, which increases with the evaporation rate. This suggests that one may introduce a simple modification of the Gibbs' criterion for (de)pinning that accounts for the nonequilibrium effect of evaporation.

Solidification in soft-core fluids: Disordered solids from fast solidification fronts

Using dynamical density functional theory we calculate the speed of solidification fronts advancing into a quenched two-dimensional model fluid of soft-core particles. We find that solidification fronts can advance via two different mechanisms, depending on the depth of the quench. For shallow quenches, the front propagation is via a nonlinear mechanism. For deep quenches, front propagation is governed by a linear mechanism and in this regime we are able to determine the front speed via a marginal stability analysis. We find that the density modulations generated behind the advancing front have a characteristic scale that differs from the wavelength of the density modulation in thermodynamic equilibrium, i.e., the spacing between the crystal planes in an equilibrium crystal. This leads to the subsequent development of disorder in the solids that are formed. In a one-component fluid, the particles are able to rearrange to form a well-ordered crystal, with few defects. However, solidification fronts in a binary mixture exhibiting crystalline phases with square and hexagonal ordering generate solids that are unable to rearrange after the passage of the solidification front and a significant amount of disorder remains in the system.

Stability of liquid films covered by a carpet of self-propelled surfactant particles

We consider a carpet of self-propelled particles at the liquid-gas interface of a liquid film on a solid substrate. The particles exert an excess pressure on the interface and also move along the interface while the swimming direction changes due to rotational diffusion. We study the intricate influence of these self-propelled insoluble surfactants on the stability of the film surface and show that depending on the strength of in-surface rotational diffusion and the absolute value of the in-surface swimming velocity, several characteristic instability modes can occur. In particular, rotational diffusion can either stabilize the film or induce instabilities of different character.

Continuous and discontinuous dynamic unbinding transitions in drawn film flow

When a plate is withdrawn from a liquid bath a coating layer is deposited whose thickness and homogeneity depend on the velocity and the wetting properties of the plate. Using a long-wave mesoscopic hydrodynamic description that incorporates wettability via a Derjaguin (disjoining) pressure we identify four qualitatively different dynamic transitions between microscopic and macroscopic coatings that are out-of-equilibrium equivalents of known equilibrium unbinding transitions. Namely, these are continuous and discontinuous dynamic wetting and emptying transitions. Several of their features have no equivalent at equilibrium.

Collapsed heteroclinic snaking near a heteroclinic chain in dragged meniscus problems

A liquid film is studied that is deposited onto a flat plate that is inclined at a constant angle to the horizontal and is extracted from a liquid bath at a constant speed. We analyse steady-state solutions of a long-wave evolution equation for the film thickness. Using centre manifold theory, we first obtain an asymptotic expansion of solutions in the bath region. The presence of an additional temperature gradient along the plate that induces a Marangoni shear stress significantly changes these expansions and leads to the presence of logarithmic terms that are absent otherwise. Next, we numerically obtain steady solutions and analyse their behaviour as the plate velocity is changed. We observe that the bifurcation curve exhibits collapsed (or exponential) heteroclinic snaking when the plate inclination angle is above a certain critical value. Otherwise, the bifurcation curve is monotonic. The steady profiles along these curves are characterised by a foot-like structure that is formed close to the meniscus and is preceded by a thin precursor film further up the plate. The length of the foot increases along the bifurcation curve. Finally, we prove with a Shilnikov-type method that the snaking behaviour of the bifurcation curves is caused by the existence of an infinite number of heteroclinic orbits close to a heteroclinic chain that connects in an appropriate three-dimensional phase space the fixed point corresponding to the precursor film with the fixed point corresponding to the foot and then with the fixed point corresponding to the bath.

Patterned deposition at moving contact lines

When a simple or complex liquid recedes from a smooth solid substrate it often leaves a homogeneous or structured deposit behind. In the case of a receding non-volatile pure liquid the deposit might be a liquid film or an arrangement of droplets depending on the receding speed of the meniscus and the wetting properties of the system. For complex liquids with volatile components as, e.g., polymer solutions and particle or surfactant suspensions, the deposit might be a homogeneous or structured layer of solute--with structures ranging from line patterns that can be orthogonal or parallel to the receding contact line via hexagonal or square arrangements of drops to complicated hierarchical structures. We review a number of recent experiments and modelling approaches with a particular focus on mesoscopic hydrodynamic long-wave models. The conclusion highlights open question and speculates about future developments.

Coarsening modes of clusters of aggregating particles

There are two modes by which clusters of aggregating particles can coalesce: The clusters can merge either (i) by the Ostwald ripening process, in which particles diffuse from one cluster to the other while the cluster centers remain stationary, or (ii) by means of a cluster translation mode, in which the clusters move toward each other and join. To understand in detail the interplay between these different modes, we study a model system of hard particles with an additional attraction between them. The particles diffuse along narrow channels with smooth or periodically corrugated walls, so that the system may be treated as one-dimensional. When the attraction between the particles is strong enough, they aggregate to form clusters. The channel potential influences whether clusters can move easily or not through the system and can prevent cluster motion. We use dynamical density functional theory to study the dynamics of the aggregation process, focusing in particular on the coalescence of two equal-sized clusters. As long as the particle hard-core diameter is nonzero, we find that the coalescence process can be halted by a sufficiently strong corrugation potential. The period of the potential determines the size of the final stable clusters. For the case of smooth channel walls, we demonstrate that there is a crossover in the dominance of the two different coarsening modes, which depends on the strength of the attraction between particles, the cluster sizes, and the separation distance between clusters.

Gradient dynamics description for films of mixtures and suspensions: dewetting triggered by coupled film height and concentration fluctuations

A thermodynamically consistent gradient dynamics model for the evolution of thin layers of liquid mixtures, solutions, and suspensions on solid substrates is presented which is based on a film-height- and mean-concentration-dependent free energy functional. It is able to describe a large variety of structuring processes, including coupled dewetting and decomposition processes. As an example, the model is employed to investigate the dewetting of thin films of liquid mixtures and suspensions under the influence of effective long-range van der Waals forces that depend on solute concentration. The occurring fluxes are discussed, and it is shown that spinodal dewetting may be triggered through the coupling of film height and concentration fluctuations. Fully nonlinear calculations provide the time evolution and resulting steady film height and concentration profiles.

Effect of Au nanoparticle spatial distribution on the stability of thin polymer films

The stability of thin poly(methyl-methacrylate) (PMMA) films of low molecular weight on a solid substrate is controlled by the areal coverage of gold nanoparticles (NPs) present at the air-polymer interface. As the polymer becomes liquid the Au NPs are free to diffuse, coalesce, and aggregate while the polymer film can change its morphology through viscous flow. These processes lead at the same time to the formation of a fractal network of Au NPs and to the development of spinodal instabilities of the free surface of the polymer films. For thinner films a single wavelength is observed, while for thicker films two wavelengths compete. With continued heating the aggregation process results in a decrease in coverage, the networks evolve into disordered particle assemblies, while the polymer films flatten again. The disordering occurs first on the smallest scales and coincides (in thicker films) with the disappearance of the smaller wavelength. The subsequent disordering on larger scales causes the films to flatten.

Localized states in the conserved Swift-Hohenberg equation with cubic nonlinearity

The conserved Swift-Hohenberg equation with cubic nonlinearity provides the simplest microscopic description of the thermodynamic transition from a fluid state to a crystalline state. The resulting phase field crystal model describes a variety of spatially localized structures, in addition to different spatially extended periodic structures. The location of these structures in the temperature versus mean order parameter plane is determined using a combination of numerical continuation in one dimension and direct numerical simulation in two and three dimensions. Localized states are found in the region of thermodynamic coexistence between the homogeneous and structured phases, and may lie outside of the binodal for these states. The results are related to the phenomenon of slanted snaking but take the form of standard homoclinic snaking when the mean order parameter is plotted as a function of the chemical potential, and are expected to carry over to related models with a conserved order parameter.