Further details on the phase diagram of hard ellipsoids of revolution

Dense Disordered Jammed Packings of Hard Spherocylinders with a Low Aspect Ratio: A Characterization of Their Structure

This work numerically investigates dense disordered (maximally random) jammed packings of hard spherocylinders of cylinder length L and diameter D by focusing on L/D ∈ [0,2]. It is within this interval that one expects that the packing fraction of these dense disordered jammed packings ϕMRJ hsc attains a maximum. This work confirms the form of the graph ϕMRJ hsc versus L/D: here, comparably to certain previous investigations, it is found that the maximal ϕMRJ hsc = 0.721 ± 0.001 occurs at L/D = 0.45 ± 0.05. Furthermore, this work meticulously characterizes the structure of these dense disordered jammed packings via the special pair-correlation function of the interparticle distance scaled by the contact distance and the ensuing analysis of the statistics of the hard spherocylinders in contact: here, distinctly from all previous investigations, it is found that the dense disordered jammed packings of hard spherocylinders with 0.45 ≲ L/D ≤ 2 are isostatic.

Enhancement in the diffusivity of Brownian spheroids in the presence of spheres

In the present paper, we have extended the simulation technique Brownian cluster dynamics (BCD) to analyze the dynamics of the binary mixture of hard ellipsoids and spheres. The shape dependent diffusional properties have been incorporated into BCD using Perrin's factor and compared with analytical results of a one-component ellipsoidal system. We have investigated pathways to enhance the diffusivity of spheroids in the binary mixture by manipulating the phase behavior of the system through varying the fraction of spheres in the binary mixture. We show that at low volume fraction the spherical particles have a higher diffusion coefficient than the ellipsoids due to the higher friction coefficient. However, at a higher volume fraction, we show that the diffusion coefficient of the ellipsoids increases irrespective of the aspect ratio due to the anisotropic shape.

Equation of state of hard lenses: A combined virial series and simulation approach

We provide highly accurate equation-of-state data for the isotropic phase of hard lenses obtained by means of cluster Monte Carlo simulations. This data is analyzed using a virial approach considering coefficients up to the order eight and Carnahan-Starling type closure relations for the virial series. The comparison with previously investigated systems consisting of hard, oblate ellipsoids of revolution allows insights into the detailed influence of the particle geometry. We propose a generalized Carnahan-Starling approach as a heuristic equation of state for the isotropic phase of hard lenses that in first approximation shows the same dependence on the excess part of the excluded volume as identified for oblate, hard lenses of revolution.

Reexamining equations of state of oblate hard ellipsoids of revolution: Numerical simulation utilizing a cluster Monte Carlo algorithm and comparison to virial theory

We provide highly accurate equation-of-state data determined by means of cluster Monte Carlo simulations for the isotropic phase of oblate hard ellipsoids of revolution. Both equation-of-state data and phase boundaries of the isotropic phase are obtained from relatively large ensembles with typically 1000 particles. The comparison of simulation data with a virial approach gives evidence for the importance of high-order so-far-unknown virial coefficients and therewith many-particle interactions in dense, isotropic systems of anisotropic particles. While a virial approach with a rescaled Carnahan-Starling correction for the unknown, higher-order virial coefficients reproduces the simulation data of moderately anisotropic particles with high accuracy, we suggest for highly anisotropic shapes a simple, heuristic equation of state as a suitable approach.

Calculation of third to eighth virial coefficients of hard lenses and hard, oblate ellipsoids of revolution employing an efficient algorithm

We provide third to eighth virial coefficients of oblate, hard ellipsoids of revolution and hard lenses in dependence on their aspect ratio ν. Employing an algorithm optimized for hard anisotropic shapes, highly accurate data are accessible with comparatively small numerical effort. For both geometries, reduced virial coefficients B[over ̃]_{i}(ν)=B_{i}(ν)/B_{2}^{i-1}(ν) are in first approximation proportional to the inverse excess contribution α^{-1} of their excluded volume. The latter quantity is directly accessible from second virial coefficients and analytically known for convex bodies.

Observation of liquid glass in suspensions of ellipsoidal colloids

Despite the omnipresence of colloidal suspensions, little is known about the influence of colloid shape on phase transformations, especially in nonequilibrium. To date, real-space imaging results at high concentrations have been limited to systems composed of spherical colloids. In most natural and technical systems, however, particles are nonspherical, and their structural dynamics are determined by translational and rotational degrees of freedom. Using confocal microscopy of fluorescently labeled core-shell particles, we reveal that suspensions of ellipsoidal colloids form an unexpected state of matter, a liquid glass in which rotations are frozen while translations remain fluid. Image analysis unveils hitherto unknown nematic precursors as characteristic structural elements of this state. The mutual obstruction of these ramified clusters prevents liquid crystalline order. Our results give insight into the interplay between local structures and phase transformations. This helps to guide applications such as self-assembly of colloidal superstructures and also gives evidence of the importance of shape on the glass transition in general.

Shape-Anisotropy-Induced Ordered Packings in Cylindrical Confinement

Densest possible packings of identical spheroids in cylindrical confinement have been obtained through Monte Carlo simulations. By varying the shape anisotropy of spheroids and also the cylinder-to-spheroid size ratio, a variety of densest possible crystalline structures have been discovered, including achiral structures with specific orientations of particles and chiral helical structures with rotating orientations of particles. Our findings reveal a transition between confinement-induced chiral ordering and shape-anisotropy-induced orientational ordering and would serve as a guide for the fabrication of crystalline wires using anisotropic particles.

Formation of nematic order in 3D systems of hard colloidal ellipsoids

Suspensions of hard ellipsoidal particles exhibit complex phase behavior as shown by theoretical predictions and simulations of phase diagrams. Here, we report quantitative confocal microscopy experiments of hard prolate colloidal ellipsoids with different aspect ratio a/b. We studied different volume fractions φ of ellipsoids in density and refractive index matched suspensions. Large 3D sample volumes were investigated and the positions as well as the orientations of all ellipsoids were extracted by image analysis routines. By evaluating the translational and orientational order in the system we determined the presence of isotropic and nematic phases. For ellipsoids with a/b = 2.0 we found that isotropic phases form at all φ, while ellipsoids with a/b = 7.0 formed nematic phases at high φ, as expected from theory and simulations. For a/b = 3.5 and a/b = 4.1, however, we observed the absence of long-range orientational order even at φ where nematic phases are expected. We show that local orientational order formed with the emergence of nematic precursors for a/b = 3.5 and short-ranged nematic domains for a/b = 4.1. Our results provide novel insight into the phase behavior and orientational order of ellipsoids with different aspect ratios.

Hard convex lens-shaped particles: Characterization of dense disordered packings

Among the family of hard convex lens-shaped particles (lenses), the one with aspect ratio equal to 2/3 is "optimal" in the sense that the maximally random jammed (MRJ) packings of such lenses achieve the highest packing fraction ϕ_{MRJ}≃0.73 [G. Cinacchi and S. Torquato, Soft Matter 14, 8205 (2018)1744-683X10.1039/C8SM01519H]. This value is only a few percent lower than ϕ_{DKP}=0.76210⋯, the packing fraction of the corresponding densest-known crystalline (degenerate) packings [G. Cinacchi and S. Torquato, J. Chem. Phys. 143, 224506 (2015)JCPSA60021-960610.1063/1.4936938]. By exploiting the appreciably reduced propensity that a system of such optimal lenses has to positionally and orientationally order, disordered packings of them are progressively generated by a Monte Carlo method-based procedure from the dilute equilibrium isotropic fluid phase to the dense nonequilibrium MRJ state. This allows us to closely monitor how the (micro)structure of these packings changes in the process of formation of the MRJ state. The gradual changes undergone by the many structural descriptors calculated here can coherently and consistently be traced back to the gradual increase in contacts between the hard particles until the isostatic mean value of ten contact neighbors per lens is reached at the effectively hyperuniform MRJ state. Compared to the MRJ state of hard spheres, the MRJ state of such optimal lenses is denser (less porous), more disordered, and rattler-free. This set of characteristics makes them good glass formers. It is possible that this conclusion may also hold for other hard convex uniaxial particles with a correspondingly similar aspect ratio, be they oblate or prolate, and that, by using suitable biaxial variants of them, that set of characteristics might further improve.

Potential and limits of a colloid approach to protein solutions

Looking at globular proteins with the eyes of a colloid scientist has a long tradition, in fact a significant part of the early colloid literature was focused on protein solutions. However, it has also been recognized that proteins are much more complex than the typical hard sphere-like synthetic model colloids. Proteins are not perfect spheres, their interaction potentials are in general not isotropic, and using theories developed for such particles are thus clearly inadequate in many cases. In this perspective article, we now take a closer look at the field. In particular, we reflect on the fact that modern colloid science has been undergoing a tremendous development, where a multitude of novel systems have been developed in the lab and in silico. During the last decade we have seen a rapidly increasing number of reports on the synthesis of anisotropic, patchy and/or responsive synthetic colloids, that start to resemble their complex biological counterparts. This experimental development is also reflected in a corresponding theoretical and simulation effort. The experimental and theoretical toolbox of colloid science has thus rapidly expanded, and there is obviously an enormous potential for an application of these new concepts to protein solutions, which has already been realized and harvested in recent years. In this perspective article we make an attempt to critically discuss the exploitation of colloid science concepts to better understand protein solutions. We not only consider classical applications such as the attempt to understand and predict solution stability and phase behaviour, but also discuss new challenges related to the dynamics, flow behaviour and liquid-solid transitions found in concentrated or crowded protein solutions. It not only aims to provide an overview on the progress in experimental and theoretical (bio)colloid science, but also discusses current shortcomings in our ability to correctly reproduce and predict the structural and dynamic properties of protein solutions based on such a colloid approach.

Planar Anchoring of C70 Liquid Crystals Using a Covalent Organic Framework Template

The surface-induced anchoring effect is a well-developed technique to control the growth of liquid crystals (LCs). Nevertheless, a defined nanometer-scale template has never been used to induce the anchored growth of LCs with molecular building units. Scanning tunneling microscopy results at the solid/liquid interface reveal that a 2D covalent organic framework (COF-1) can offer an anchoring effect to template C₇₀ molecules into forming several LC mesophases, which cannot be obtained under other conditions. Through comparison with the C₆₀ system, a stepwise breakdown in ordering of C₇₀ LC is observed. The process is described in terms of the effects of molecular anisotropy on the epitaxial growth of molecular crystals. The results suggest that using a surface-confined template to anchor the initial layer of LC molecules can be a modular and potentially broadly applicable approach for organizing molecular mesogens into LCs.

Phase behavior of hard C_{2h}-symmetric particle systems

Using Monte Carlo numerical simulation, this work sketches the phase diagram of systems of certain hard C_{2h}-symmetric particles, formed by gluing two aligned and displaced hard spherocylinders with a cylindrical-length-to-diameter ratio realistically, if viewed not only from the lyotropic colloidal liquid-crystal side but also from the thermotropic low-molecular-mass liquid-crystal side, equal to 5, as a function of the displacement. Several distinctive phases are observed, such as a nonperiodic smectic-B-like phase, a nonperiodic smectic-H-like phase, a smectic-C phase, and a short-layer-spacing uniaxial smectic-A phase but no biaxial nematic phase.

Hard convex lens-shaped particles: metastable, glassy and jammed states

We generate and study dense positionally and/or orientationally disordered, including jammed, monodisperse packings of hard convex lens-shaped particles (lenses). Relatively dense isotropic fluid configurations of lenses of various aspect ratios are slowly compressed via a Monte Carlo method based procedure. Under this compression protocol, while 'flat' lenses form a nematic fluid phase (where particles are positionally disordered but orientationally ordered) and 'globular' lenses form a plastic solid phase (where particles are positionally ordered but orientationally disordered), 'intermediate', neither 'flat' nor 'globular', lenses do not form either mesophase. In general, a crystal solid phase (where particles are both positionally and orientationally ordered) does not spontaneously form during lengthy numerical simulation runs. In correspondence to those volume fractions at which a transition to the crystal solid phase would occur in equilibrium, a 'downturn' is observed in the inverse compressibility factor versus volume fraction curve beyond which this curve behaves essentially linearly. This allows us to estimate the volume fraction at jamming of the dense non-crystalline packings so generated. These packings are nematic for 'flat' lenses and plastic for 'globular' lenses, while they are robustly isotropic for 'intermediate' lenses, as confirmed by the calculation of the τ order metric, among other quantities. The structure factors S(k) of the corresponding jammed states tend to zero as the wavenumber k goes to zero, indicating they are effectively hyperuniform (i.e., their infinite-wavelength density fluctuations are anomalously suppressed). Among all possible lens shapes, 'intermediate' lenses with aspect ratio around 2/3 are special because they are those that reach the highest volume fractions at jamming while being positionally and orientationally disordered and these volume fractions are as high as those reached by nematic jammed states of 'flat' lenses and plastic jammed states of 'globular' lenses. All of their attributes, taken together, make such 'intermediate' lens packings particularly good glass-forming materials.

Phase diagram of hard squares in slit confinement

This work shows a complete phase diagram of hard squares of side length σ in slit confinement for H < 4.5, H being the wall to wall distance measured in σ units, including the maximal packing fraction limit. The phase diagram exhibits a transition between a single-row parallel 1-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{\square }}$$\end{document}□ and a zigzag 2-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hat{\diamond }$$\end{document}◇ˆ structures for H_c(2) = (2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{{\bf{2}}}$$\end{document}2 − 1) < H < 2, and also another one involving the 1-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{\square }}$$\end{document}□ and 2-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{\square }}$$\end{document}□ structures (two parallel rows) for 2 < H < H_c(3) (H_c(n) = n − 1 + \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{{\bf{2}}{\boldsymbol{n}}-{\bf{1}}}$$\end{document}2n−1/n is the critical wall-to-wall distance for a (n − 1)-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{\square }}$$\end{document}□ to n-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\diamond }$$\end{document}◇ transition and where n-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\diamond }$$\end{document}◇ represents a structure formed by tilted rectangles, each one clustering n stacked squares), and a triple point for H_t \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{\simeq }}$$\end{document}≃ 2.005. In this triple point there coexists the 1-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{\square }}$$\end{document}□, 2-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{\square }}$$\end{document}□, and 2-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hat{\diamond }$$\end{document}◇ˆ structures. For regions H_c(3) < H < H_c(4) and H_c(4) < H < H_c(5), very similar pictures arise. There is a (n − 1)-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{\square }}$$\end{document}□ to a n-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\diamond }$$\end{document}◇ strong transition for H_c(n) < H < n, followed by a softer (n − 1)-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{\square }}$$\end{document}□ to n-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{\square }}$$\end{document}□ transition for n < H < H_c(n + 1). Again, at H \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{\gtrsim }}$$\end{document}≳ n there appears a triple point, involving the (n − 1)-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{\square }}$$\end{document}□, n-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{\square }}$$\end{document}□, and n-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\diamond }$$\end{document}◇ structures. The similarities found for n = 2, 3 and 4 lead us to propose a tentative phase diagram for H_c(n) < H < H_c(n + 1) (n ∈ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{{\mathbb{N}}}}$$\end{document}ℕ, n > 2), where structures (n − 1)-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{\square }}$$\end{document}□, n-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{\square }}$$\end{document}□, and n-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\diamond }$$\end{document}◇ fill the phase diagram. Simulation and Onsager theory results are qualitatively consistent.

Elastic properties of the nematic phase in hard ellipsoids of short aspect ratio

We present a Monte Carlo simulation study of suspensions of hard ellipsoids of revolution. Based on the spatial fluctuations of the orientational order, we have computed the Frank elastic constants for prolate and oblate ellipsoids and compared them to the affine transformation model. The affine transformation model predicts the right order of magnitude of the twist and bend constant but not of the splay constant. In addition, we report the observation of a stable nematic phase at an aspect ratio as low as 2.5.

Dense crystalline packings of ellipsoids

An ellipsoid, the simplest nonspherical shape, has been extensively used as a model for elongated building blocks for a wide spectrum of molecular, colloidal, and granular systems. Yet the densest packing of congruent hard ellipsoids, which is intimately related to the high-density phase of many condensed matter systems, is still an open problem. We discover an unusual family of dense crystalline packings of self-dual ellipsoids (ratios of the semiaxes α:sqrt[α]:1), containing 24 particles with a quasi-square-triangular (SQ-TR) tiling arrangement in the fundamental cell. The associated packing density ϕ exceeds that of the densest known SM2 crystal [ A. Donev et al., Phys. Rev. Lett. 92, 255506 (2004)10.1103/PhysRevLett.92.255506] for aspect ratios α in (1.365, 1.5625), attaining a maximal ϕ≈0.75806... at α=93/64. We show that the SQ-TR phase derived from these dense packings is thermodynamically stable at high densities over the aforementioned α range and report a phase diagram for self-dual ellipsoids. The discovery of the SQ-TR crystal suggests organizing principles for nonspherical particles and self-assembly of colloidal systems.

Critical behavior of hard squares in strong confinement

We examine the phase behavior of a quasi-one-dimensional system of hard squares with side-length σ, where the particles are confined between two parallel walls and only nearest-neighbor interactions occur. As in our previous work [Gurin, Varga, and Odriozola, Phys. Rev. E 94, 050603 (2016)]2470-004510.1103/PhysRevE.94.050603, the transfer operator method is used, but here we impose a restricted orientation and position approximation to yield an analytic description of the physical properties. This allows us to study the parallel fluid-like to zigzag solid-like structural transition, where the compressibility and heat capacity peaks sharpen and get higher as H→H_{c}=2sqrt[2]-1≈1.8284 and p→p_{c}=∞. Here H is the width of the channel measured in σ units and p is the pressure. We have found that this structural change becomes critical at the (p_{c},H_{c}) point. The obtained critical exponents belong to the universality class of the one-dimensional Ising model. We believe this behavior holds for the unrestricted orientational and positional case.

Anisotropic magnetic particles in a magnetic field

We characterize the structural properties of magnetic ellipsoidal hematite colloids with an aspect ratio ρ ≈ 2.3 using a combination of small-angle X-ray scattering and computer simulations. The evolution of the phase diagram with packing fraction ϕ and the strength of an applied magnetic field B is described, and the coupling between orientational order of magnetic ellipsoids and the bulk magnetic behavior of their suspension addressed. We establish quantitative structural criteria for the different phase and arrest transitions and map distinct isotropic, polarized non-nematic, and nematic phases over an extended range in the ϕ-B coordinates. We show that upon a rotational arrest of the ellipsoids around ϕ = 0.59, the bulk magnetic behavior of their suspension switches from superparamagnetic to ordered weakly ferromagnetic. If densely packed and arrested, these magnetic particles thus provide persisting remanent magnetization of the suspension. By exploring structural and magnetic properties together, we extend the often used colloid-atom analogy to the case of magnetic spins.

Effect of orientational restriction on monolayers of hard ellipsoids

The effect of out-of-plane orientational freedom on the orientational ordering properties of a monolayer of hard ellipsoids is studied using the Parsons-Lee scaling approach and replica exchange Monte Carlo computer simulation. Prolate and oblate ellipsoids exhibit very different ordering properties, namely, the axes of revolution of prolate particles tend to lean out, while those of oblate ones prefer to lean into the confining plane. The driving mechanism of this is that the particles try to maximize the available free area on the confining surface, which can be achieved by minimizing the cross section areas of the particles with the plane. In the lack of out-of-plane orientational freedom the monolayer of prolate particles is identical to a two-dimensional hard ellipse system, which undergoes an isotropic-nematic ordering transition with increasing density. With gradually switching on the out-of-plane orientational freedom the prolate particles lean out from the confining plane and destabilisation of the in-plane isotropic-nematic phase transition is observed. The system of oblate particles behaves oppositely to that of prolates. It corresponds to a two-dimensional system of hard disks in the lack of out-of-plane freedom, while it behaves similar to that of hard ellipses in the freely rotating case. Solid phases can be realised by lower surface coverage due to the out-of-plane orientation freedom for both oblate and prolate shapes.

Nematic-smectic transition of parallel hard spheroellipsoids

Spheroellipsoids are truncated ellipsoids with spherical end caps. If gradients are assumed to change smoothly at the junction of body and cap, the truncation height z0 determines the geometry uniquely. The resulting model particle has only two shape parameters, namely, the aspect ratio c/a of the basic ellipsoid and the cutoff z0/a. These two parameters can be tuned to yield a continuous transformation between a pure ellipsoid and a spherocylinder. Since parallel hard spherocylinders display a nematic-smectic A phase transition, while ellipsoids do not, the influence of the particle shape on the possibility of a smectic phase may be investigated. A density functional analysis is used to detect the dividing line, in the (c/a, z0/a) plane, between the presence and absence of the N-S transition. Since spheroellipsoids may be useful as generic model particles for anisotropic molecules, we provide a computationally efficient overlap criterion for a pair in a general, non-parallel configuration.

Strong effect of weak charging in suspensions of anisotropic colloids

Suspensions of hard colloidal particles frequently serve as model systems in studies on fundamental aspects of phase transitions. But often colloidal particles that are considered as "hard" are in fact weakly charged. If the colloids are spherical, weak charging has only a weak effect on the structural properties of the suspension, which can be easily corrected for. However, this does not hold for anisotropic particles. We introduce a model for the interaction potential between charged ellipsoids of revolution (spheroids) based on the Derjaguin approximation of Debye-Hückel theory and present a computer simulation study on aspects of the system's structural properties and phase behaviour. In line with previous experimental observations, we find that even a weak surface charge has a strong impact on the correlation functions. A likewise strong impact is seen on the phase behaviour, in particular, we find stable cubatic order in suspensions of oblate ellipsoids.