Source author record

Andrew J. Archer

Andrew J. Archer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

26 published item(s)

preprint2022arXiv

Smectic Layering: Landau theory for a complex-tensor order parameter

Composed of microscopic layers that stack along one direction while maintaining fluid-like positional disorder within layers, smectics are excellent systems for exploring topology, defects and geometric memory in complex confining geometries. However, the coexistence of crystalline-like characteristics in one direction and fluid-like disorder within layers makes lamellar liquid crystals notoriously difficult to model - especially in the presence of defects and large distortions. Nematic properties of smectics can be comprehensively described by the Q-tensor but to capture the features of the smectic layering alone, we develop a phenomenological Landau theory for a complex-tensor order parameter E, which is capable of describing the local degree of lamellar ordering, layer displacement, and orientation of the layers. This theory can account for both parallel and perpendicular elastic contributions. In addition to resolving the potential ambiguities inherent to complex scalar order parameter models, this model reduces to previous employed models of simple smectics, and opens new possibilities for numerical studies on smectics possessing many defects, within complex geometries and under extreme confinement.

preprint2021arXiv

Dynamical density functional theory for the drying and stratification of binary colloidal dispersions

We develop a dynamical density functional theory based model for the drying of colloidal films on planar surfaces. We consider mixtures of two different sizes of hard-sphere colloids. Depending on the solvent evaporation rate and the initial concentrations of the two species, we observe varying degrees of stratification in the final dried films. Our model predicts the various structures described in the literature previously from experiments and computer simulations, in particular the small-on-top stratified films. Our model also includes the influence of adsorption of particles to the interfaces.

preprint2021arXiv

How ice grows from premelting films and water droplets

Close to the triple point, the surface of ice is covered by a thin liquid layer (so-called quasi-liquid layer) which crucially impacts growth and melting rates. Experimental probes cannot observe the growth processes below this layer, and classical models of growth by vapor deposition do not account for the formation of premelting films. Here, we develop a mesoscopic model of liquid-film mediated ice growth, and identify the various resulting growth regimes. At low saturation, freezing proceeds by terrace spreading, but the motion of the buried solid is conveyed through the liquid to the outer liquid-vapor interface. At higher saturations water droplets condense, a large crater forms below, and freezing proceeds undetectably beneath the droplet. Our approach is a general framework that naturally models freezing close to three phase coexistence and provides a first principle theory of ice growth and melting which may prove useful in the geosciences.

preprint2020arXiv

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. 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.

preprint2020arXiv

First order phase transitions and the thermodynamic limit

We consider simple mean field continuum models for first order liquid-liquid demixing and solid-liquid phase transitions and show how the Maxwell construction at phase coexistence emerges on going from finite-size closed systems to the thermodynamic limit. The theories considered are the Cahn-Hilliard model of phase separation, which is also a model for the liquid-gas transition, and the phase field crystal model of the solid-liquid transition. Our results show that states comprising the Maxwell line depend strongly on the mean density with spatially localized structures playing a key role in the approach to the thermodynamic limit.

preprint2020arXiv

Pattern formation in two-dimensional hard-core/soft-shell systems with variable soft shell profiles

Hard-core/soft shell (HCSS) particles have been shown to self-assemble into a remarkably rich variety of structures under compression due to the simple interplay between the hard-core and soft-shoulder length scales in their interactions. Most studies in this area model the soft shell interaction as a square shoulder potential. Although appealing from a theoretical point of view, the potential is physically unrealistic because there is no repulsive force in the soft shell regime, unlike in experimental HCSS systems. To make the model more realistic, here we consider HCSS particles with a range soft shell potential profiles beyond the standard square shoulder form and study the model using both minimum energy calculations and Monte Carlo simulations. We find that by tuning density and the soft shell profile, HCSS particles in the thin shell regime (i.e., shell to core ratio $r_1/r_0 \leq \sqrt{3}$) can form a large range of structures, including hexagons, chains, squares, rhomboids and two distinct zig-zag structures. Furthermore, by tuning the density and $r_1/r_0$, we find that HCSS particles with experimentally realistic linear ramp soft shoulder repulsions can form honeycombs and quasicrystals with 10-fold and 12-fold symmetry. Our study therefore suggests the exciting possibility of fabricating these exotic 2D structures experimentally through colloidal self-assembly.

preprint2020arXiv

Rounded Layering Transitions on the Surface of Ice

Understanding the wetting properties of premelting films requires knowledge of the film's equation of state, which is not usually available. Here we calculate the disjoining pressure curve of premelting films, and perform a detailed thermodynamic characterization of premelting behavior on ice. Analysis of the density profiles reveals the signature of weak layering phenomena, from one to two and from two to three water molecular layers. However, disjoining pressure curves, which closely follow expectations from a renormalized mean field liquid state theory, show that there are no layering phase transitions in the thermodynamic sense along the sublimation line. Instead, we find that transitions at mean field level are rounded due to capillary wave fluctuations. We see signatures that true first order layering transitions could arise at low temperatures, for pressures between the metastable line of water/vapor coexistence and the sublimation line. The extrapolation of the disjoining pressure curve above water vapor saturation displays a true first order phase transition from a thin to a thick film consistent with experimental observations.

preprint2016arXiv

Crystallisation of soft matter under confinement at interfaces and in wedges

The surface freezing and surface melting transitions exhibited by a model two-dimensional soft matter system is studied. The behaviour when confined within a wedge is also considered. The system consists of particles interacting via a soft purely repulsive pair potential. Density functional theory (DFT) is used to calculate density profiles and thermodynamic quantities. The external potential due to the confining walls is modelled via a hard-wall with an additional repulsive Yukawa potential. The surface phase behaviour depends on the range and strength of this repulsion: When the repulsion strength is weak, the wall promotes freezing at the surface of the wall. The thickness of this frozen layer grows logarithmically as the bulk liquid-solid phase coexistence is approached. Our mean-field DFT predicts that this crystalline layer at the wall must be nucleated (i.e. there is a free energy barrier) and its formation is necessarily a first-order transition, referred to as `prefreezing', by analogy with the prewetting transition. However, in contrast to the latter, prefreezing cannot terminate in a critical point, since the phase transition involves a change in symmetry. If the wall-fluid interaction is sufficiently long ranged and the repulsion is strong enough, surface melting can instead occur. Then the interface between the wall and the bulk crystalline solid becomes wet by the liquid phase as the chemical potential is decreased towards the value at liquid-solid coexistence. It is observed that the finite thickness fluid film at the wall has a broken translational symmetry due to its proximity to the bulk crystal and so the nucleation of the wetting film can be either first-order or continuous. Our mean-field theory predicts that for certain wall potentials there is a premelting critical point analogous to the surface critical point for the prewetting transition. In a wedge...

preprint2016arXiv

Gradient dynamics models for liquid films with soluble surfactant

In this paper we propose equations of motion for the dynamics of liquid films of surfactant suspensions that consist of a general gradient dynamics framework based on an underlying energy functional. This extends the gradient dynamics approach to dissipative non-equilibrium thin film systems with several variables, and casts their dynamic equations into a form that reproduces Onsager's reciprocity relations. We first discuss the general form of gradient dynamics models for an arbitrary number of fields and discuss simple well-known examples with one or two fields. Next, we develop the gradient dynamics (three field) model for a thin liquid film covered by soluble surfactant and discuss how it automatically results in consistent convective (driven by pressure gradients, Marangoni forces and Korteweg stresses), diffusive, adsorption/desorption, and evaporation fluxes. We then show that in the dilute limit, the model reduces to the well-known hydrodynamic form that includes Marangoni fluxes due to a linear equation of state. In this case the energy functional incorporates wetting energy, surface energy of the free interface (constant contribution plus an entropic term) and bulk mixing entropy. Subsequently, as an example, we show how various extensions of the energy functional result in consistent dynamical models that account for nonlinear equations of state, concentration-dependent wettability and surfactant and film bulk decomposition phase transitions. We conclude with a discussion of further possible extensions towards systems with micelles, surfactant adsorption at the solid substrate and bioactive behaviour.

preprint2015arXiv

Liquid drops on a surface: using density functional theory to calculate the binding potential and drop profiles and comparing with results from mesoscopic modelling

The contribution to the free energy for a film of liquid of thickness $h$ on a solid surface, due to the interactions between the solid-liquid and liquid-gas interfaces is given by the binding potential, $g(h)$. The precise form of $g(h)$ determines whether or not the liquid wets the surface. Note that differentiating $g(h)$ gives the Derjaguin or disjoining pressure. We develop a microscopic density functional theory (DFT) based method for calculating $g(h)$, allowing us to relate the form of $g(h)$ to the nature of the molecular interactions in the system. We present results based on using a simple lattice gas model, to demonstrate the procedure. In order to describe the static and dynamic behaviour of non-uniform liquid films and drops on surfaces, a mesoscopic free energy based on $g(h)$ is often used. We calculate such equilibrium film height profiles and also directly calculate using DFT the corresponding density profiles for liquid drops on surfaces. Comparing quantities such as the contact angle and also the shape of the drops, we find good agreement between the two methods. We also study in detail the effect on $g(h)$ of truncating the range of the dispersion forces, both those between the fluid molecules and those between the fluid and wall. We find that truncating can have a significant effect on $g(h)$ and the associated wetting behaviour of the fluid.

preprint2015arXiv

Two-dimensional colloidal fluids exhibiting pattern formation

Fluids with competing short range attraction and long range repulsive interactions between the particles can exhibit a variety of microphase separated structures. We develop a lattice-gas (generalised Ising) model and analyse the phase diagram using Monte Carlo computer simulations and also with density functional theory (DFT). The DFT predictions for the structures formed are in good agreement with the results from the simulations, which occur in the portion of the phase diagram where the theory predicts the uniform fluid to be linearly unstable. However, the mean-field DFT does not correctly describe the transitions between the different morphologies, which the simulations show to be analogous to micelle formation. We determine how the heat capacity varies as the model parameters are changed. There are peaks in the heat capacity at state points where the morphology changes occur. We also map the lattice model onto a continuum DFT that facilitates a simplification of the stability analysis of the uniform fluid.

preprint2014arXiv

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 whilst the cluster centres remain stationary, or (ii) by means of a cluster translation mode, in which the clusters move towards 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 size clusters. As long as the particle hard-core diameter is non-zero, 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 cross-over in the dominance of the two different coarsening modes, that depends on the strength of the attraction between particles, the cluster sizes and the separation distance between clusters.

preprint2013arXiv

An Introduction to Microscopic Theories for Inhomogeneous Liquids: Getting Started with Density Functional Theory

Classical density functional theory (DFT) is a statistical mechanical theory for calculating the density profiles of the molecules in a liquid. It is widely used, for example. to calculate the density distribution of the molecules in the vicinity of a confining wall, the interfacial tension, the wetting behaviour and many other properties of nonuniform liquids. DFT can however be somewhat daunting to students entering the field, because of the many connections to other areas of liquid-state science that are required and used to develop the theories. Here we give an introduction to some of the key ideas, based on a lattice-gas (Ising) model fluid. This builds on knowledge covered in most undergraduate statistical mechanics and thermodynamics courses and so students can quickly get to the stage of calculating density profiles, etc for themselves. We derive a simple DFT for the lattice-gas and present some typical results that can readily be calculated using the theory.

preprint2013arXiv

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.

preprint2013arXiv

Note on the hydrodynamic description of thin nematic films: strong anchoring model

We discuss the long-wave hydrodynamic model for a thin film of nematic liquid crystal in the limit of strong anchoring at the free surface and at the substrate. We rigorously clarify how the elastic energy enters the evolution equation for the film thickness in order to provide a solid basis for further investigation: several conflicting models exist in the literature that predict qualitatively different behaviour. We consolidate the various approaches and show that the long-wave model derived through an asymptotic expansion of the full nemato-hydrodynamic equations with consistent boundary conditions agrees with the model one obtains by employing a thermodynamically motivated gradient dynamics formulation based on an underlying free energy functional. As a result, we find that in the case of strong anchoring the elastic distortion energy is always stabilising. To support the discussion in the main part of the paper, an appendix gives the full derivation of the evolution equation for the film thickness via asymptotic expansion.

preprint2013arXiv

Sedimentation of a two-dimensional colloidal mixture exhibiting liquid-liquid and gas-liquid phase separation: a dynamical density functional theory study

We present dynamical density functional theory results for the time evolution of the density distribution of a sedimenting model two-dimensional binary mixture of colloids. The interplay between the bulk phase behaviour of the mixture, its interfacial properties at the confining walls, and the gravitational field gives rise to a rich variety of equilibrium and non-equilibrium morphologies. In the fluid state, the system exhibits both liquid-liquid and gas-liquid phase separation. As the system sediments, the phase separation significantly affects the dynamics and we explore situations where the final state is a coexistence of up to three different phases. Solving the dynamical equations in two-dimensions, we find that in certain situations the final density profiles of the two species have a symmetry that is different from that of the external potentials, which is perhaps surprising, given the statistical mechanics origin of the theory. The paper concludes with a discussion on this.

preprint2012arXiv

Modelling the formation of structured deposits at receding contact lines of evaporating solutions and suspensions

When a film of a liquid suspension of nanoparticles or a polymer solution is deposited on a surface, it may dewet from the surface and as the solvent evaporates the solute particles/polymer can be deposited on the surface in regular line patterns. In this paper we explore a hydrodynamic model for the process that is based on a long-wave approximation that predicts the deposition of irregular and regular line patterns. This is due to a self-organised pinning-depinning cycle that resembles a stick-slip motion of the contact line. We present a detailed analysis of how the line pattern properties depend on quantities such as the evaporation rate, the solute concentration, the Péclet number, the chemical potential of the ambient vapour, the disjoining pressure, and the intrinsic viscosity. The results are related to several experiments and to depinning transitions in other soft matter systems.

preprint2012arXiv

Thermodynamically consistent description of the hydrodynamics of free surfaces covered by insoluble surfactants of high concentration

In this paper we propose several models that describe the dynamics of liquid films which are covered by a high concentration layer of insoluble surfactant. First, we briefly review the 'classical' hydrodynamic form of the coupled evolution equations for the film height and surfactant concentration that are well established for small concentrations. Then we re-formulate the basic model as a gradient dynamics based on an underlying free energy functional that accounts for wettability and capillarity. Based on this re-formulation in the framework of nonequilibrium thermodynamics, we propose extensions of the basic hydrodynamic model that account for (i) nonlinear equations of state, (ii) surfactant-dependent wettability, (iii) surfactant phase transitions, and (iv) substrate-mediated condensation. In passing, we discuss important differences to most of the models found in the literature.

preprint2011arXiv

Modelling fluids and crystals using a two-component modified phase field crystal model

A modified phase field crystal model in which the free energy may be minimised by an order parameter profile having isolated bumps is investigated. The phase diagram is calculated in one and two dimensions and we locate the regions where modulated and uniform phases are formed and also regions where localised states are formed. We investigate the effectiveness of the phase field crystal model for describing fluids and crystals with defects. We further consider a two component model and elucidate how the structure transforms from hexagonal crystalline ordering to square ordering as the concentration changes. Our conclusion contains a discussion of possible interpretations of the order parameter field.

preprint2011arXiv

Ratcheting of driven attracting colloidal particles: Temporal density oscillations and current multiplicity

We consider the unidirectional particle transport in a suspension of colloidal particles which interact with each other via a pair potential having a hard-core repulsion plus an attractive tail. The colloids are confined within a long narrow channel and are driven along by a DC or an AC external potential. In addition, the walls of the channel interact with the particles via a ratchet-like periodic potential. We use dynamical density functional theory to compute the average particle current. In the case of DC drive, we show that as the attraction strength between the colloids is increased beyond a critical value, the stationary density distribution of the particles loses its stability leading to depinning and a time dependent density profile. Attraction induced symmetry breaking gives rise to the coexistence of stable stationary density profiles with different spatial periods and time-periodic density profiles, each characterized by different values for the particle current.

preprint2010arXiv

Collective shuttling of attracting particles in asymmetric narrow channels

The rectification of a single file of attracting particles subjected to a low frequency ac drive is proposed as a working mechanism for particle shuttling in an asymmetric narrow channel. Increasing the particle attraction results in the file condensing, as signalled by the dramatic enhancement of the net particle current. Magnitude and direction of the current become extremely sensitive to the actual size of the condensate, which can then be made to shuttle between two docking stations, transporting particles in one direction, with an efficiency much larger than conventional diffusive models predict.

preprint2010arXiv

Dynamical model for the formation of patterned deposits at receding contact lines

We describe the formation of deposition patterns that are observed in many different experiments where a three-phase contact line of a volatile nanoparticle suspension or polymer solution recedes. A dynamical model based on a long-wave approximation predicts the deposition of irregular and regular line patterns due to self-organised pinning-depinning cycles corresponding to a stick-slip motion of the contact line. We analyze how the line pattern properties depend on the evaporation rate and solute concentration.

preprint2010arXiv

On the interplay between sedimentation and phase separation phenomena in two-dimensional colloidal fluids

Colloidal particles that are confined to an interface effectively form a two-dimensional fluid. We examine the dynamics of such colloids when they are subject to a constant external force, which drives them in a particular direction over the surface. Such a situation occurs, for example, for colloidal particles that have settled to the bottom of their container, when the container is tilted at an angle, so that they `sediment' to the lower edge of the surface. We focus in particular on the case when there are attractive forces between the colloids which causes them to phase separate into regions of high density and low density and we study the influence of this phase separation on the sedimentation process. We model the colloids as Brownian particles and use both Brownian dynamics computer simulations and dynamical density functional theory (DDFT) to obtain the time evolution of the ensemble average one-body density profiles of the colloids. We consider situations where the external potential varies only in one direction so that the ensemble average density profiles vary only in this direction. We solve the DDFT in one-dimension, by assuming that the density profile only varies in one direction. However, we also solve the DDFT in two-dimensions, allowing the fluid density profile to vary in both the $x$- and $y$-directions. We find that in certain situations the two-dimensional DDFT is clearly superior to its one-dimensional counterpart when compared with the simulations and we discuss this issue.

preprint2009arXiv

Selectivity in binary fluid mixtures: static and dynamical properties

Selectivity of particles in a region of space can be achieved by applying external potentials to influence the particles in that region. We investigate static and dynamical properties of size selectivity in binary fluid mixtures of two particles sizes. We find that by applying an external potential that is attractive to both kinds of particles, due to crowding effects, this can lead to one species of particles being expelled from that region, whilst the other species is attracted into the region where the potential is applied. This selectivity of one species of particle over the other in a localized region of space depends on the density and composition of the fluid mixture. Applying an external potential that repels both kinds of particles leads to selectivity of the opposite species of particles to the selectivity with attractive potentials. We use equilibrium and dynamical density functional theory to describe and understand the static and dynamical properties of this striking phenomenon. Selectivity by some ion-channels is believed to be due to this effect.

preprint2009arXiv

Solvent mediated interactions between model colloids and interfaces: A microscopic approach

We determine the solvent mediated contribution to the effective potentials for model colloidal or nano- particles dispersed in a binary solvent that exhibits fluid-fluid phase separation. Using a simple density functional theory we calculate the density profiles of both solvent species in the presence of the `colloids', which are treated as external potentials, and determine the solvent mediated (SM) potentials. Specifically, we calculate SM potentials between (i) two colloids, (ii) a colloid and a planar fluid-fluid interface, and (iii) a colloid and a planar wall with an adsorbed wetting film. We consider three different types of colloidal particles: colloid A which prefers the bulk solvent phase rich in species 2, colloid C which prefers the solvent phase rich in species 1, and `neutral' colloid B which has no strong preference for either phase, i.e. the free energies to insert the colloid into either of the coexisting bulk phases are almost equal. When a colloid which has a preference for one of the two solvent phases is inserted into the disfavored phase at statepoints close to coexistence a thick adsorbed `wetting' film of the preferred phase may form around the colloids. The presence of the adsorbed film has a profound influence on the form of the SM potentials.

preprint2006arXiv

Dynamics in inhomogeneous liquids and glasses via the test particle limit

We show that one may view the self and the distinct part of the van Hove dynamic correlation function of a simple fluid as the one-body density distributions of a binary mixture that evolve in time according to dynamical density functional theory. For a test case of soft core Brownian particles the theory yields results for the van Hove function that agree quantitatively with those of our Brownian dynamics computer simulations. At sufficiently high densities the free energy landscape underlying the dynamics exhibits a barrier as a function of the mean particle displacement, shedding new light on the nature of glass formation. For hard spheres confined between parallel planar walls the barrier height oscillates in-phase with the local density, implying that the mobility is maximal between layers, which should be experimentally observable in confined colloidal dispersions.