Phase behaviour of hard board-like particles

Structural and dynamical equilibrium properties of hard board-like particles in parallel confinement

We performed Monte Carlo and dynamic Monte Carlo simulations to model the diffusion of monodispersed suspensions composed of impenetrable cuboidal particles, specifically hard board-like particles (HBPs), in the presence of parallel hard walls. The impact of the walls was investigated by adjusting the size of the simulation box while maintaining constant packing fractions, fixed at η = 0.150, for systems consisting of HBPs with prolate, dual-shaped, and oblate geometries. We observed that increasing the distance between the walls led to the recovery of an isotropic bulk phase, while local particle organization near the walls remained stable. Due to their shape, oblate HBPs exhibit more efficient anchoring at wall surfaces compared to prolate shapes. The formation of nematic-like particle assemblies near the walls, confirmed by theoretical calculations based on density functional theory, significantly influenced local particle dynamics. This effect was particularly pronounced to the extent that a modest portion of cuboids near the walls tended to diffuse exclusively in planes parallel to the confinement, even more efficiently than observed in the bulk regions.

Effect of particle anisotropy on the thermodynamics and kinetics of ordering transitions in hard faceted particles

Monte Carlo simulations were used to study the influence of particle aspect ratio on the kinetics and phase behavior of hard gyrobifastigia (GBF). First, the formation of a highly anisotropic nucleus shape in the isotropic-to-crystal transition in regular GBF is explained by the differences in interfacial free energies of various crystal planes and the nucleus geometry predicted by the Wulff construction. GBF-related shapes with various aspect ratios were then studied, mapping their equations of state, determining phase coexistence conditions via interfacial pinning, and computing nucleation free-energy barriers via umbrella sampling using suitable order parameters. Our simulations reveal a reduction of the kinetic barrier for isotropic-crystal transition upon an increase in aspect ratio, and that for highly oblate and prolate aspect ratios, an intermediate nematic phase is stabilized. Our results and observations also support two conjectures for the formation of the crystalline state from the isotropic phase: that low phase free energies at the ordering phase transition correlate with low transition barriers and that the emergence of a mesophase provides a steppingstone that expedites crystallization.

Active microrheology of colloidal suspensions of hard cuboids

By performing dynamic Monte Carlo simulations, we investigate the microrheology of isotropic suspensions of hard-core colloidal cuboids. In particular, we infer the local viscoelastic behavior of these fluids by studying the dynamics of a probe spherical particle that is incorporated in the host phase and is dragged by an external force. This technique, known as active microrheology, allows one to characterize the microscopic response of soft materials upon application of a constant force, whose intensity spans here three orders of magnitude. By tuning the geometry of cuboids from oblate to prolate as well as the system density, we observe different responses that are quantified by measuring the effective friction perceived by the probe particle. The resulting friction coefficient exhibits a linear regime at forces that are much weaker and larger than the thermal forces, whereas a nonlinear, force-thinning regime is observed at intermediate force intensities.

Hard core lattice gas with third next-nearest neighbor exclusion on triangular lattice: One or two phase transitions?

We obtain the phase diagram of the hard core lattice gas with third nearest neighbor exclusion on the triangular lattice using Monte Carlo simulations that are based on a rejection-free flat histogram algorithm. In a recent paper [Darjani et al., J. Chem. Phys. 151, 104702 (2019)], it was claimed that the lattice gas with third nearest neighbor exclusion undergoes two phase transitions with increasing density with the phase at intermediate densities exhibiting hexatic order with continuously varying exponents. Although a hexatic phase is expected when the exclusion range is large, it has not been seen earlier in hard core lattice gases with short range exclusion. In this paper, by numerically determining the entropies for all densities, we show that there is only a single phase transition in the system between a low-density fluid phase and a high density ordered sublattice phase and that a hexatic phase is absent. The transition is shown to be first order in nature, and the critical parameters are determined accurately.

Rejection-free cluster Wang-Landau algorithm for hard-core lattice gases

We introduce a rejection-free, flat histogram, cluster algorithm to determine the density of states of hard-core lattice gases. We show that the algorithm is able to efficiently sample low entropy states that are usually difficult to access, even when the excluded volume per particle is large. The algorithm is based on simultaneously evaporating all the particles in a strip and reoccupying these sites with a new appropriately chosen configuration. We implement the algorithm for the particular case of the hard-core lattice gas in which the first k next-nearest neighbors of a particle are excluded from being occupied. It is shown that the algorithm is able to reproduce the known results for k=1,2,3 both on the square and cubic lattices. We also show that, in comparison, the corresponding flat histogram algorithms with either local moves or unbiased cluster moves are less accurate and do not converge as the system size increases.

Tailoring the phase diagram of discotic mesogens

The computational modelling of discotic molecules is a central topic in colloid science that is key for the smart design of a broad range of modern functional materials. This work lays out a versatile interaction model capable of exposing the rich mesogenic behaviour of discotics. A single coarse-grained spheroplatelet core framework is employed to generate a variety of pair interaction anisotropy classes, favouring specific relative orientations of the particles (stacked, side-side, crossed, T-shaped). This paves the way for the systematic tailoring of the discotic liquid phase diagram. Monte Carlo simulations are performed for an ensemble of case studies to illustrate the correlation between the topology of the interaction and the formation of stable nematic, smectic and columnar phases, as well as of less common cubatic, uniaxial and biaxial columnar domains.

Real-space renormalization-group treatment of the Maier-Saupe-Zwanzig model for biaxial nematic structures

The Maier-Saupe-Zwanzig model for the nematic phase transitions in liquid crystals is investigated in a diamond hierarchical lattice. The model takes into account a parameter to describe the biaxiality of the microscopic units. Also, a suitably chosen external field is added to the Hamiltonian to allow the determination of critical parameters associated with the nematic phase transitions. Using the transfer-matrix technique, the free energy and its derivatives are obtained in terms of recursion relations between successive generations of the hierarchical lattice. In addition, a real-space renormalization-group approach is developed to obtain the critical parameters of the same model system. Results of both methods are in excellent agreement. There are indications of two continuous phase transitions. One of them corresponds to a uniaxial-isotropic transition, in the class of universality of the three-state Potts model on the diamond hierarchical lattice. The transition between the biaxial and the uniaxial phases is in the universality class of the Ising model on the same lattice.

Orientational ordering in a fluid of hard kites: A density-functional-theory study

Using density-functional theory we theoretically study the orientational properties of uniform phases of hard kites-two isosceles triangles joined by their common base. Two approximations are used: scaled particle theory and a new approach that better approximates third virial coefficients of two-dimensional hard particles. By varying some of their geometrical parameters, kites can be transformed into squares, rhombuses, triangles, and also very elongated particles, even reaching the hard-needle limit. Thus, a fluid of hard kites, depending on the particle shape, can stabilize isotropic, nematic, tetratic, and triatic phases. Different phase diagrams are calculated, including those of rhombuses, and kites with two of their equal interior angles fixed to 90^{∘}, 60^{∘}, and 75^{∘}. Kites with one of their unequal angles fixed to 72^{∘}, which have been recently studied via Monte Carlo simulations, are also considered. We find that rhombuses and kites with two equal right angles and not too large anisometry stabilize the tetratic phase but the latter stabilize it to a much higher degree. By contrast, kites with two equal interior angles fixed to 60^{∘} stabilize the triatic phase to some extent, although it is very sensitive to changes in particle geometry. Kites with the two equal interior angles fixed to 75^{∘} have a phase diagram with both tetratic and triatic phases, but we show the nonexistence of a particle shape for which both phases are stable at different densities. Finally, the success of the new theory in the description of orientational order in kites is shown by comparing with Monte Carlo simulations for the case where one of the unequal angles is fixed to 72^{∘}. These particles also present a phase diagram with stable tetratic and triatic phases.

Self-assembly of freely-rotating polydisperse cuboids: unveiling the boundaries of the biaxial nematic phase

Colloidal cuboids have the potential to self-assemble into biaxial liquid crystal phases, which exhibit two independent optical axes. Over the last few decades, several theoretical works have predicted the existence of a wide region of the phase diagram where the biaxial nematic phase would be stable, but imposed rather strong constraints on the particle rotational degrees of freedom. In this work, we use molecular simulation to investigate the impact of size dispersity on the phase behaviour of freely-rotating hard cuboids, here modelled as self-dual-shaped nanoboards. This peculiar anisotropy, exactly in between the oblate and prolate geometry, has been proposed as the most appropriate to promote phase biaxiality. We observe that size dispersity radically changes the phase behaviour of monodisperse systems and leads to the formation of an elusive biaxial nematic phase, being found in a large region of the packing fraction vs. polydispersity phase diagram. Although our results confirm the tendencies reported in past experimental observations on colloidal dispersions of slightly prolate goethite particles, they cannot reproduce the direct isotropic-to-biaxial nematic phase transition observed in these experiments.

Dynamics of hard colloidal cuboids in nematic liquid crystals

We perform dynamic Monte Carlo simulations to investigate the equilibrium dynamics of hard board-like colloidal particles in oblate and prolate nematic liquid crystals. In particular, we characterize the particles' diffusion along the nematic director and perpendicularly to it, and observe a structural relaxation decay that strongly depends on the particle anisotropy. To assess the Gaussianity of their dynamics and eventual occurrence of collective motion, we calculate two- and four-point correlation functions that incorporate the instantaneous values of the diffusion coefficients parallel and perpendicular to the nematic director. Our simulation results highlight the occurrence of Fickian and Gaussian dynamics at short and long times, locate the minimum diffusivity at the self-dual shape, the particle geometry that would preferentially stabilise biaxial nematics, and exclude the existence of dynamically correlated particles.

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.

Fluctuations and phase transitions of uniaxial and biaxial liquid crystals using a theoretically informed Monte Carlo and a Landau free energy density

In this work, we explore fluctuations during phase transitions of uniaxial and biaxial liquid crystals using a phenomenological free energy functional. We rely on a continuum-level description of the liquid crystal ordering with a tensorial parameter and a temperature dependent Landau polynomial expansion of the tensor's invariants. The free energy functional, over a three-dimensional periodic domain, is integrated with a Gaussian quadrature and minimized with a theoretically informed Monte Carlo method. We reconstruct analytical phase diagrams, following Landau and Doi's notations, to verify that the free energy relaxation reaches the global minimum. Importantly, our relaxation method is able to follow the thermodynamic behavior provided by other non-phenomenological approaches; we predict the first order character of the isotropic-nematic transition, and we identify the uniaxial-biaxial transition as second order. Finally, we use a finite-size scaling method, using the nematic susceptibility, to calculate the transition temperatures for 4-Cyano-4'-pentylbiphenyl (5CB) and N-(4-methoxybenzylidene)-4-butylaniline (MBBA). Our results show good agreement with experimental values, thereby validating our minimization method. Our approach is an alternative towards the relaxation of temperature dependent continuum-level free energy functionals, in any geometry, and can incorporate complicated elastic and surface energy densities.

Biaxial nematics of hard cuboids in an external field

By computer simulation, we model the phase behaviour of colloidal suspensions of board-like particles under the effect of an external field and assess the still disputed occurrence of the biaxial nematic (NB) liquid crystal phase. The external field promotes the rearrangement of the initial isotropic (I) or uniaxial nematic (NU) phase and the formation of the NB phase. In particular, very weak field strengths are sufficient to spark a direct I-NB or NU-NB phase transition at the self-dual shape, where prolate and oblate particle geometries fuse into one. By contrast, forming the NB phase at any other geometry requires stronger fields and thus reduces the energy efficiency of the phase transformation. Our simulation results show that self-dual shaped board-like particles with moderate anisotropy are able to form NB liquid crystals under the effect of a surprisingly weak external stimulus and suggest a path to exploit low-energy uniaxial-to-biaxial order switching.

Phase behaviour of colloidal superballs mixed with non-adsorbing polymers

Inspired by experimental work on colloidal cuboid-polymer dispersions (Rossi et al., Soft Matter, 7, 4139 (2011)) we have theoretically studied the phase behaviour of such mixtures. To that end, free volume theory (FVT) was applied to predict the phase behaviour of mixtures of superballs and non-adsorbing polymer chains in a common solvent. Closed expressions for the thermodynamic properties of a suspension of hard colloidal superballs have been derived, accounting for fluid (F), face-centred cubic (FCC) and simple cubic (SC) phase states. Even though the considered solid phases are approximate, the hard superballs phase diagram semi-quantitatively matches with more evolved methods. The theory developed for the cuboid-polymer mixture reveals a rich phase behaviour, which includes not only isostructural F₁-F₂ coexistence, but also SC₁-SC₂ coexistence, several triple coexistences, and even a quadruple-phase coexistence region (F₁-F₂-SC-FCC). The model proposed offers a tool to asses the stability of cuboid-polymer mixtures in terms of the colloid-to-polymer size ratio.

Hierarchical bounding structures for efficient virial computations: Towards a realistic molecular description of cholesterics

We detail the application of bounding volume hierarchies to accelerate second-virial evaluations for arbitrary complex particles interacting through hard and soft finite-range potentials. This procedure, based on the construction of neighbour lists through the combined use of recursive atom-decomposition techniques and binary overlap search schemes, is shown to scale sub-logarithmically with particle resolution in the case of molecular systems with high aspect ratios. Its implementation within an efficient numerical and theoretical framework based on classical density functional theory enables us to investigate the cholesteric self-assembly of a wide range of experimentally relevant particle models. We illustrate the method through the determination of the cholesteric behavior of hard, structurally resolved twisted cuboids, and report quantitative evidence of the long-predicted phase handedness inversion with increasing particle thread angles near the phenomenological threshold value of 45°. Our results further highlight the complex relationship between microscopic structure and helical twisting power in such model systems, which may be attributed to subtle geometric variations of their chiral excluded-volume manifold.

Hard Competition: Stabilizing the Elusive Biaxial Nematic Phase in Suspensions of Colloidal Particles with Extreme Lengths

We use computer simulations to study the existence and stability of a biaxial nematic N_{b} phase in systems of hard polyhedral cuboids, triangular prisms, and rhombic platelets, characterized by a long (L), medium (M), and short (S) particle axis. For all three shape families, we find stable N_{b} states provided the shape is not only close to the so-called dual shape with M=sqrt[LS] but also sufficiently anisotropic with L/S>9,11,14,23 for rhombi, (two types of) triangular prisms, and cuboids, respectively, corresponding to anisotropies not considered before. Surprisingly, a direct isotropic-N_{b} transition does not occur in these systems due to a destabilization of N_{b} by a smectic (for cuboids and prisms) or a columnar (for platelets) phase at small L/S or by an intervening uniaxial nematic phase at large L/S. Our results are confirmed by a density functional theory provided the third virial coefficient is included and a continuous rather than a discrete (Zwanzig) set of particle orientations is taken into account.

Entropic patchiness drives multi-phase coexistence in discotic colloid-depletant mixtures

Entropy-driven equilibrium phase behaviour of hard particle dispersions can be understood from excluded volume arguments only. While monodisperse hard spheres only exhibit a fluid-solid phase transition, anisotropic hard particles such as rods, discs, cuboids or boards exhibit various multi-phase equilibria. Ordering of such anisotropic particles increases the free volume entropy by reducing the excluded volume between them. The addition of depletants gives rise to an entropic patchiness represented by orientation-dependent attractions resulting in non-trivial phase behaviour. We show that free volume theory is a simple, generic and tractable framework that enables to incorporate these effects and rationalise various experimental findings. Plate-shaped particles constitute the main building blocks of clays, asphaltenes and chromonic liquid crystals that find widespread use in the food, cosmetics and oil industry. We demonstrate that mixtures of platelets and ideal depletants exhibit a strikingly rich phase behaviour containing several types of three-phase coexistence areas and even a quadruple region with four coexisting phases.