Matthieu Marechal

Phase behavior and structure of a new colloidal model system of bowl-shaped particles

We study the phase behavior of bowl-shaped (nano)particles using confocal microscopy and computer simulations. Experimentally, we find the formation of a wormlike fluid phase in which the bowl-shaped particles have a strong tendency to stack on top of each other. However, using free energy calculations in computer simulations, we show that the wormlike phase is out-of-equilibrium and that the columnar phase is thermodynamically stable for sufficiently deep bowls and high densities. In addition, we employ a novel technique based on simulated annealing to predict the crystal structures for shallow bowls. We find four exotic new crystal structures and we determine their region of stability using free energy calculations. We discuss the implications of our results for the development of materials consisting of molecular mesogens or nanoparticles.

Vacancy-stabilized crystalline order in hard cubes

We examine the effect of vacancies on the phase behavior and structure of systems consisting of hard cubes using event-driven molecular dynamics and Monte Carlo simulations. We find a first-order phase transition between a fluid and a simple cubic crystal phase that is stabilized by a surprisingly large number of vacancies, reaching a net vacancy concentration of approximately 6.4% near bulk coexistence. Remarkably, we find that vacancies increase the positional order in the system. Finally, we show that the vacancies are delocalized and therefore hard to detect.

Phase behavior of hard colloidal platelets using free energy calculations

We investigate the phase behavior of a model for colloidal hard platelets and rigid discotic molecules: oblate hard spherocylinders (OHSC). We perform free energy calculations using Monte Carlo simulations to map out the phase diagram as a function of the aspect ratio L∕D of the particles. The phase diagram displays a stable isotropic phase, a nematic liquid crystal phase for L∕D≤0.12, a columnar phase for L∕D≲0.3, a tilted crystal phase for L≲0.45, and an aligned crystal phase for L∕D≳0.45. We compare the results to the known phase diagram of hard cut spheres. Thin cut spheres are almost cylinder-shaped, while the interactions between real discotic mesogens and colloidal platelets are more consistent with the toroidal rims of the OHSC. Since the shapes of the OHSC and the cut spheres are otherwise similar, the phase diagrams of the two types of particles are quite akin. However, the tilted crystal phase for OHSC, which is of a crystal type that is frequently found in experiments on disklike molecules, has not been found for hard cut spheres. Furthermore, although we have found a cubatic phase, it was shown to be definitely unstable, whereas the stability of the cubatic phase of cut spheres is still disputed. Finally, we also show that the phase boundaries differ significantly from those for cut spheres. These are remarkable consequences of a subtle change in particle shape, which show that for a detailed comparison with the phase behavior of experimental particles, the OHSC should be used as a model particle.

Crystal-structure prediction via the Floppy-Box Monte Carlo algorithm: Method and application to hard (non)convex particles

In this paper, we describe the way to set up the floppy-box Monte Carlo (FBMC) method [L. Filion, M. Marechal, B. van Oorschot, D. Pelt, F. Smallenburg, and M. Dijkstra, Phys. Rev. Lett. 103, 188302 (2009)] to predict crystal-structure candidates for colloidal particles. The algorithm is explained in detail to ensure that it can be straightforwardly implemented on the basis of this text. The handling of hard-particle interactions in the FBMC algorithm is given special attention, as (soft) short-range and semi-long-range interactions can be treated in an analogous way. We also discuss two types of algorithms for checking for overlaps between polyhedra, the method of separating axes and a triangular-tessellation based technique. These can be combined with the FBMC method to enable crystal-structure prediction for systems composed of highly shape-anisotropic particles. Moreover, we present the results for the dense crystal structures predicted using the FBMC method for 159 (non)convex faceted particles, on which the findings in [J. de Graaf, R. van Roij, and M. Dijkstra, Phys. Rev. Lett. 107, 155501 (2011)] were based. Finally, we comment on the process of crystal-structure prediction itself and the choices that can be made in these simulations.

Phase behavior and structure of colloidal bowl-shaped particles: simulations.

We study the phase behavior of bowl-shaped particles using computer simulations. These particles were found experimentally to form a metastable wormlike fluid phase in which the bowl-shaped particles have a strong tendency to stack on top of each other [M. Marechal, Nano Lett. 10, 1907 (2010)]. In this work, we show that the transition from the low-density fluid to the wormlike phase has an interesting effect on the equation of state. The simulation results also show that the wormlike fluid phase transforms spontaneously into a columnar phase for bowls that are sufficiently deep. Furthermore, we describe the phase behavior as obtained from free energy calculations employing Monte Carlo simulations. The columnar phase is stable for bowl shapes ranging from infinitely thin bowls to surprisingly shallow bowls. Aside from a large region of stability for the columnar phase, the phase diagram features four novel crystal phases and a region where the stable fluid contains wormlike stacks.

Freezing of parallel hard cubes with rounded edges

The freezing transition in a classical three-dimensional system of rounded hard cubes with fixed, equal orientations is studied by computer simulation and fundamental-measure density functional theory. By switching the rounding parameter s from zero to one, one can smoothly interpolate between cubes with sharp edges and hard spheres. The equilibrium phase diagram of rounded parallel hard cubes is computed as a function of their volume fraction and the rounding parameter s. The second order freezing transition known for oriented cubes at s = 0 is found to be persistent up to s = 0.65. The fluid freezes into a simple-cubic crystal which exhibits a large vacancy concentration. Upon a further increase of s, the continuous freezing is replaced by a first-order transition into either a sheared simple cubic lattice or a deformed face-centered cubic lattice with two possible unit cells: body-centered orthorhombic or base-centered monoclinic. In principle, a system of parallel cubes could be realized in experiments on colloids using advanced synthesis techniques and a combination of external fields.

Stacking in sediments of colloidal hard spheres

We use computer simulations to investigate the crystallization dynamics of sedimenting hard spheres in large systems (hundreds of thousands of particles). We show that slow sedimentation results primarily in face-centered cubic (fcc) stacked crystals, instead of random hexagonal close packed or hexagonal close packed (hcp) crystals. We also find slanted stacking faults, in the fcc regions. However, we attribute the formation of fcc to the free energy difference between fcc and hcp and not to the presence of these slanted stacking faults. Although the free energy difference between hcp and fcc per particle is small (only 10(-3) times the thermal energy), it can become considerable, when multiplied by the number of particles in each domain. The ratio of fcc to hcp obtained from dynamic simulations is in excellent agreement with well-equilibrated Monte Carlo simulations, in which no slanted stacking faults were found. Our results explain a range of experiments on colloids, in which the amount of fcc increases upon lowering the sedimentation rate or decreasing the initial volume fraction.

Packing Confined Hard Spheres Denser with Adaptive Prism Phases

We show that hard spheres confined between two parallel hard plates pack denser with periodic adaptive prismatic structures which are composed of alternating prisms of spheres. The internal structure of the prisms adapts to the slit height which results in close packings for a range of plate separations, just above the distance where three intersecting square layers fit exactly between the plates. The adaptive prism phases are also observed in real-space experiments on confined sterically stabilized colloids and in Monte Carlo simulations at finite pressure.

Inhomogeneous fluids of colloidal hard dumbbells: Fundamental measure theory and Monte Carlo simulations

Recently, a density functional theory for hard particles with shape anisotropy was developed, the extended deconvolution fundamental measure theory (edFMT). We apply edFMT to hard dumbbells, arguably the simplest non-convex shape and readily available experimentally in the form of colloids. We obtain good agreement between edFMT and Monte Carlo simulations for fluids of dumbbells in a slit and for the same system under gravity. This indicates that edFMT can be successfully applied to nearly all colloidal shapes, not just for the convex shapes for which edFMT was originally derived. A theory, such as edFMT, that allows a fast and general way of mapping the phase behavior of anisotropic colloids, can act as a useful guide for the design of colloidal shapes for various applications.

Frustration of the Isotropic-Columnar Phase Transition of Colloidal Hard Platelets by a Transient Cubatic Phase

Using simulations and theory, we show that the cubatic phase is metastable for three model hard platelets. The locally favored structures of perpendicular particle stacks in the fluid prevent the formation of the columnar phase through geometric frustration resulting in vitrification. Also, we find a direct link between structure and dynamic heterogeneities in the cooperative rotation of particle stacks, which is crucial for the devitrification process. Finally, we show that the lifetime of the glassy cubatic phase can be tuned by surprisingly small differences in particle shape.