Swim pressure on walls with curves and corners

The concept of swim pressure quantifies the average force exerted by microswimmers on confining walls in non-equilibrium. Here we explore how the swim pressure depends on the wall curvature and on the presence of sharp corners in the wall. For active Brownian particles at high dilution, we present a coherent framework which describes the force and torque on passive particles of arbitrary shape, in the limit of large particles compared to the persistence length of the swimmer trajectories. The resulting forces can be used to derive, for example, the activity-induced depletion interaction between two disks, as well as to optimize the shape of a tracer particle for high swimming velocity. Our predictions are verifiable in experiments on passive obstacles exposed to a bath of bacteria or artificial microswimmers

Tuning the liquid-liquid transition by modulating the hydrogen bond angular flexibility in a model for water

We propose a simple extension of the well known ST2 model for water [F.H. Stillinger and A. Rahman, J. Chem. Phys. {\bf 60}, 1545 (1974)] that allows for a continuous modification of the hydrogen bond angular flexibility. We show that the bond flexibility affects the relative thermodynamic stability of the liquid and of the hexagonal (or cubic) ice. On increasing flexibility, the liquid-liquid critical point, which in the original ST2 model is located in the no-man's land (i. e. the region where ice is the thermodynamically stable phase) progressively moves to a temperature where the liquid is more stable than ice. Our study definitively proves that the liquid-liquid transition in ST2 is a genuine phenomenon, of high relevance in all tetrahedral network-forming liquids, including water.Comment: Accepted in Phys. Rev. Let

Understanding tetrahedral liquids through patchy colloids

We investigate the structural properties of a simple model for tetrahedral patchy colloids in which the patch width and the patch range can be tuned independently. For wide bond angles, a fully bonded network can be generated by standard Monte Carlo or molecular dynamics simulations of the model, providing a neat method for generating defect-free random tetrahedral networks. This offers the possibility of focusing on the role of the patch angular width on the structure of the fully bonded network. The analysis of the fully bonded configurations as a function of the bonding angle shows how the bonding angle controls the system compressibility, the strength of the pre-peak in the structure factor and ring size distribution. Comparison with models of liquid water and silica allows us to find the best mapping between these continuous potentials and the colloidal one. Building on previous studies focused on the connection between angular range and crystallization, the mapping makes it possible to shed new light on the glass-forming ability of network-forming tetrahedral liquids.Comment: 10 pages, 6 figure

Phase diagram of the ST2 model of water

We evaluate the free energy of the fluid and crystal phases for the ST2 potential [F.H. Stillinger and A. Rahman, J. Chem. Phys. 60, 1545 (1974)] with reaction field corrections for the long-range interactions. We estimate the phase coexistence boundaries in the temperature-pressure plane, as well as the gas-liquid critical point and gas-liquid coexistence conditions. Our study frames the location of the previously identified liquid-liquid critical point relative to the crystalline phase boundaries, and opens the way for exploring crystal nucleation in a model where the metastable liquid-liquid critical point is computationally accessible

Liquid crystals of hard rectangles on flat and cylindrical manifolds

Using the classical density functional theory of freezing and Monte Carlo computer simulations, we explore the liquid-crystalline phase behavior of hard rectangles on flat and cylindrical manifolds. Moreover, we study the effect of a static external field which couples to the rectangles' orientations, aligning them towards a preferred direction. In the flat and field-free case, the bulk phase diagram involves stable isotropic, nematic, tetratic, and smectic phases depending on the aspect ratio and number density of the particles. The external field shifts the transition lines significantly and generates a binematic phase at the expense of the tetratic phase. On a cylindrical manifold, we observe tilted smectic-like order, as obtained by wrapping a smectic layer around a cylinder. We find in general good agreement between our density functional calculations and particle-resolved computer simulations and mention possible setups to verify our predictions in experiments.Comment: 10 pages, 6 figure

Gelling by Heating

We introduce a simple model, a binary mixture of patchy particles, which has been designed to form a gel upon heating. Due to the specific nature of the particle interactions, notably the number and geometry of the patches as well as their interaction energies, the system is a fluid both at high and at low temperatures, whereas at intermediate temperatures the system forms a solid-like disordered open network structure, i.e. a gel. Using molecular dynamics we investigate the static and dynamic properties of this system

Modeling of many-body interactions between elastic spheres through symmetry functions

Simple models for spherical particles with a soft shell have been shown to self-assemble into numerous crystal phases and even quasicrystals. However, most of these models rely on a simple pairwise interaction, which is usually a valid approximation only in the limit of small deformations, i.e. low densities. In this work, we consider a many-body yet simple model for the evaluation of the elastic energy associated with the deformation of a spherical shell. The resulting energy evaluation, however, is relatively expensive for direct use in simulations. We significantly reduce the associated numerical cost by fitting the potential using a set of symmetry functions. We propose a method for selecting a suitable set of symmetry functions that capture the most relevant features of the particle environment in a systematic manner. The fitted interaction potential is then used in Monte Carlo simulations to draw the phase diagram of the system in two dimensions. The system is found to form both a fluid and a hexagonal crystal phase.Comment: 10 pages, 9 figure

Infinite-pressure phase diagram of binary mixtures of (non)additive hard disks

One versatile route to the creation of two-dimensional crystal structures on the nanometer to micrometer scale is the self-assembly of colloidal particles at an interface. Here, we explore the crystal phases that can be expected from the self-assembly of mixtures of spherical particles of two different sizes, which we map to (additive or non-additive) hard-disk mixtures. We map out the infinite-pressure phase diagram for these mixtures, using Floppy Box Monte Carlo simulations to systematically sample candidate crystal structures with up to 12 disks in the unit cell. As a function of the size ratio and number ratio of the two species of particles, we find a rich variety of periodic crystal structures. Additionally, we identify random tiling regions to predict random tiling quasicrystal stability ranges. Increasing non-additivity both gives rise to additional crystal phases and broadens the stability regime for crystal structures involving a large number of large-small contacts, including random tilings. Our results provide useful guidelines for controlling the self-assembly of colloidal particles at interfaces