Surface tension of the isotropic-nematic interface

We present the first calculations of the pressure tensor profile in the vicinity of the planar interface between isotropic liquid and nematic liquid crystal, using Onsager's density functional theory and computer simulation. When the liquid crystal director is aligned parallel to the interface, the situation of lowest free energy, there is a large tension on the nematic side of the interface and a small compressive region on the isotropic side. By contrast, for perpendicular alignment, the tension is on the isotropic side. There is excellent agreement between theory and simulation both in the forms of the pressure tensor profiles, and the values of the surface tension.Comment: Minor changes; to appear in Phys. Rev.

Molecular dynamics in arbitrary geometries : parallel evaluation of pair forces

A new algorithm for calculating intermolecular pair forces in molecular dynamics (MD) simulations on a distributed parallel computer is presented. The arbitrary interacting cells algorithm (AICA) is designed to operate on geometrical domains defined by an unstructured, arbitrary polyhedral mesh that has been spatially decomposed into irregular portions for parallelisation. It is intended for nano scale fluid mechanics simulation by MD in complex geometries, and to provide the MD component of a hybrid MD/continuum simulation. The spatial relationship of the cells of the mesh is calculated at the start of the simulation and only the molecules contained in cells that have part of their surface closer than the cut-off radius of the intermolecular pair potential are required to interact. AICA has been implemented in the open source C++ code OpenFOAM, and its accuracy has been indirectly verified against a published MD code. The same system simulated in serial and in parallel on 12 and 32 processors gives the same results. Performance tests show that there is an optimal number of cells in a mesh for maximum speed of calculating intermolecular forces, and that having a large number of empty cells in the mesh does not add a significant computational overhead

Evolution of Non-Equilibrium Profile in Adsorbate Layer under Compressive Strain

We investigate the time evolution of an initial step profile separating a bare substrate region from the rest of the compressively strained adsorbate layer near a commensurate to incommensurate transition. The rate of profile evolution as a function of the mismatch, coverage and the strength of the substrate potential are determined by Brownian molecular dynamics simulations. We find that the results are qualitatively similar to those observed for the Pb/Si(111) system. The anomalously fast time evolution and sharpness of the non-equilibrium profile can be understood through the domain wall creation at the boundary and its subsequent diffusion into the interior of the adsorbate layer.Comment: 6 pages, 7 figures, Tribology Letter

Self-assembly of twisted, multi-sheet aggregates

Hierarchical self-assembly underpins much of the diversity of form and function seen in soft systems, yet the pathways by which they achieve their final form are not always straightforward – intermediate steps, kinetic effects and finite sizes of aggregates all influence the self-assembly pathways of these systems. In this paper, we use molecular dynamics simulations of binary mixtures of spheres and ellipsoidal discs to investigate the self-assembly of anisotropic aggregates with internal structures. Through this, the full aggregation pathways of spontaneously chiral, multi-bilayer and multi-layer assemblies have been tracked and characterised via a semi-qualitative analysis. This includes the unambiguous identification of first-, second- and third-generation hierarchical assemblies within a single simulation. Given the significant challenge of tracking full aggregation pathways in experimental systems, our findings strongly support the notion that molecular simulation has much to contribute to improving our understanding of hierarchical self-assembling systems

Structure, Stresses and Local Dynamics in Glasses

The interrelations between short range structural and elastic aspects in glasses and glass forming liquids pose important and yet unresolved questions. In this paper these relations are analyzed for mono-atomic glasses and stressed liquids with a short range repulsive-attractive pair potentials. Strong variations of the local pressure are found even in a zero temperature glass, whereas the largest values of pressure are the same in both glasses and liquids. The coordination number z(J) and the effective first peak radius depend on the local pressures J's. A linear relation was found between components of site stress tensor and the local elastic constants. A linear relation was also found between the trace of the squares of the local frequencies and the local pressures. Those relations hold for glasses at zero temperature and for liquids. We explain this by a relation between the structure and the potential terms. A structural similarity between liquids and solids is manifested by similar dependencies of the coordination number on the pressures.Comment: 7 pages, 11 figure

A Model for the Elasticity of Compressed Emulsions

We present a new model to describe the unusual elastic properties of compressed emulsions. The response of a single droplet under compression is investigated numerically for different Wigner-Seitz cells. The response is softer than harmonic, and depends on the coordination number of the droplet. Using these results, we propose a new effective inter-droplet potential which is used to determine the elastic response of a monodisperse collection of disordered droplets as a function of volume fraction. Our results are in excellent agreement with recent experiments. This suggests that anharmonicity, together with disorder, are responsible for the quasi-linear increase of $G$ and $\Pi$ observed at $\varphi_c$.Comment: RevTeX with psfig-included figures and a galley macr

External and intrinsic anchoring in nematic liquid crystals: A Monte Carlo study

We present a Monte Carlo study of external surface anchoring in nematic cells with partially disordered solid substrates, as well as of intrinsic anchoring at free nematic interfaces. The simulations are based on the simple hexagonal lattice model with a spatially anisotropic intermolecular potential. We estimate the corresponding extrapolation length $b$ by imposing an elastic deformation in a hybrid cell-like nematic sample. Our estimates for $b$ increase with increasing surface disorder and are essentially temperature--independent. Experimental values of $b$ are approached only when both the coupling of nematic molecules with the substrate and the anisotropy of nematic--nematic interactions are weak.Comment: Revisions primarily in section I

Algorithm for numerical integration of the rigid-body equations of motion

A new algorithm for numerical integration of the rigid-body equations of motion is proposed. The algorithm uses the leapfrog scheme and the quantities involved are angular velocities and orientational variables which can be expressed in terms of either principal axes or quaternions. Due to specific features of the algorithm, orthonormality and unit norms of the orientational variables are integrals of motion, despite an approximate character of the produced trajectories. It is shown that the method presented appears to be the most efficient among all known algorithms of such a kind.Comment: 4 pages, 1 figur

Applying spatial reasoning to topographical data with a grounded geographical ontology

Grounding an ontology upon geographical data has been pro- posed as a method of handling the vagueness in the domain more effectively. In order to do this, we require methods of reasoning about the spatial relations between the regions within the data. This stage can be computationally expensive, as we require information on the location of points in relation to each other. This paper illustrates how using knowledge about regions allows us to reduce the computation required in an efficient and easy to understand manner. Further, we show how this system can be implemented in co-ordination with segmented data to reason abou

Simple Fluids with Complex Phase Behavior

We find that a system of particles interacting through a simple isotropic potential with a softened core is able to exhibit a rich phase behavior including: a liquid-liquid phase transition in the supercooled phase, as has been suggested for water; a gas-liquid-liquid triple point; a freezing line with anomalous reentrant behavior. The essential ingredient leading to these features resides in that the potential investigated gives origin to two effective core radii.Comment: 7 pages including 3 eps figures + 1 jpeg figur