Self-assembly of ellipsoidal particles at fluid-fluid interfaces with an empirical pair potential.

Colloidal particles adsorbed at fluid-fluid interfaces interact via mechanisms that can be specific to the presence of interfaces, for instance, lateral capillary interactions induced by nonspherical particles. Capillary interactions are highly relevant for self-assembly and the formation of surface microstructures, however, these are very challenging to model due to the multibody nature of capillary interactions. This work pursues a direct comparison between our computational modelling approach and experimental results on surface microstructures formed by ellipsoidal particles. We begin by investigating the accuracy of using pairwise interactions to describe the multibody capillary interaction by contrasting exact two- and three-particle interaction energies and we find that the pairwise approximation appears reasonable for the experimentally relevant configurations studied. We then develop an empirical pair potential and use it in Monte-Carlo type simulations to efficiently model the structure formation process for relevant particle properties such as aspect ratio, contact angle and surface coverage, and succeed in reproducing our experimental observations where we spread sterically-stabilised ellipsoidal particles onto an oil-air interface at high surface coverage. At lower surface coverages, we find that the self-assembly process falls into the diffusion-limited colloid aggregation universality class.

Reply to "Comment on 'Temperature-dependent orientational ordering on a spherical surface modeled with a lattice spin model' ".

In Romano's Comment [Phys. Rev. E 91, 046501 (2015)] on Sec. II of our paper [Phys. Rev. E 90, 022502 (2014)], the author suggests that our findings concerning the nature of the ordering transition of our modified Lebwohl-Lasher model with two-dimensional planar rotators on a planar lattice are inconsistent with known mathematical results. We argue in this Reply that our findings are in fact in agreement with previous mathematical and simulation results and that the criticisms raised by Romano have no impact on the results presented in our paper.

Modelling the rheology of anisotropic particles adsorbed on a two-dimensional fluid interface.

We present a general approach based on nonequilibrium thermodynamics for bridging the gap between a well-defined microscopic model and the macroscopic rheology of particle-stabilised interfaces. Our approach is illustrated by starting with a microscopic model of hard ellipsoids confined to a planar surface, which is intended to simply represent a particle-stabilised fluid-fluid interface. More complex microscopic models can be readily handled using the methods outlined in this paper. From the aforementioned microscopic starting point, we obtain the macroscopic, constitutive equations using a combination of systematic coarse-graining, computer experiments and Hamiltonian dynamics. Exemplary numerical solutions of the constitutive equations are given for a variety of experimentally relevant flow situations to explore the rheological behaviour of our model. In particular, we calculate the shear and dilatational moduli of the interface over a wide range of surface coverages, ranging from the dilute isotropic regime, to the concentrated nematic regime.

Nonequilibrium thermodynamics of the shear-transformation-zone model.

The shear-transformation-zone (STZ) model has been applied numerous times to describe the plastic deformation of different types of amorphous systems. We formulate this model within the general equation for nonequilibrium reversible-irreversible coupling (GENERIC) framework, thereby clarifying the thermodynamic structure of the constitutive equations and guaranteeing thermodynamic consistency. We propose natural, physically motivated forms for the building blocks of the GENERIC, which combine to produce a closed set of time evolution equations for the state variables, valid for any choice of free energy. We demonstrate an application of the new GENERIC-based model by choosing a simple form of the free energy. In addition, we present some numerical results and contrast those with the original STZ equations.

Temperature-dependent orientational ordering on a spherical surface modeled with a lattice spin model.

We study the orientational ordering on the surface of a sphere using Monte Carlo and Brownian dynamics simulations of rods interacting with an anisotropic potential. We restrict the orientations to the local tangent plane of the spherical surface and fix the position of each rod to be at a discrete point on the spherical surface. On the surface of a sphere, orientational ordering cannot be perfectly nematic due to the inevitable presence of defects. We find that the ground state of four +1/2 point defects is stable across a broad range of temperatures. We investigate the transition from disordered to ordered phase by decreasing the temperature and find a very smooth transition. We use fluctuations of the local directors to estimate the Frank elastic constant on the surface of a sphere and compare it to the planar case. We observe subdiffusive behavior in the mean square displacement of the defect cores and estimate their diffusion constants.

The Landau free energy of hard ellipses obtained from microscopic simulations.

Systems of two-dimensional hard ellipses of varying aspect ratios and packing fractions are studied by Monte Carlo simulations in the generalised canonical ensemble. From this microscopic model, we extract a coarse-grained macroscopic Landau-de Gennes free energy as a function of packing fraction and orientational order parameter. We separate the free energy into the ideal orientational entropy of non-interacting two-dimensional spins and an excess free energy associated with excluded volume interactions. We further explore the isotropic-nematic phase transition using our empirical expression for the free energy and find that the nature of the phase transition is continuous for the aspect ratios we studied.