Orientational ordering and phase behaviour of binary mixtures of hard spheres and hard spherocylinders

Excluded volume interactions and phase stability in mixtures of hard spheres and hard rods

In this paper we study excluded volume interactions, the free volume fraction available, and the phase behaviour, in mixtures of hard spheres (HS) and hard rods, modeled as spherocylinders. We use free volume theory (FVT) to predict various physical properties and compare to Monte Carlo computer simulations. FVT is used at two levels. We use the original FVT approach in which it is assumed that the correlations of the HS are not affected by the rods. This is compared to a recent, more rigorous, FVT approach which includes excluded volume interactions between the different components at all levels. We find that the novel rigorous FVT approach agrees well with computer simulation results at the level of free volume available, as well as for the phase stability. The FVT predictions show significant quantitative and qualitative deviations with respect to the original FVT approach. The phase transition curves are systematically at higher rod concentrations than previously predicted. Furthermore, the calculations revealed that a certain asphericity is required to induce isostructural fluid-fluid coexistence and the stability region is highly dependent on the size ratio between the rods and the spheres.

Non-monotonic phase behaviour of a mixture containing non-adsorbing particles and polymerising rod-like molecules

HYPOTHESIS

Previous works have shown that many-body interactions induced by dispersants with increasing correlation length will generate a diminishing two-phase region [Soft Matter 14, 6921 (2018)]. We conjecture that the attenuation of the depletion attraction due to many-body interactions is a ubiquitous phenomenon in medium-induced interactions. We propose mixtures of colloidal particles and rod-like polymers as a feasible experimental system for verifying these predictions, since the intra-molecular correlations are not screened in a good solvent for rod-like polymers as they are in flexible polymers. The length of the rods can grow and become the dominant length scale that determines the range of the depletion interactions for the imbedded non-adsorbing particles. Simulations: We study many-body depletion forces induced by polymerizing rod-like polymers on spherical non-adsorbing colloids, using Metropolis Monte Carlo simulations. We also employ a simple mean-field theory to further justify our numerical predictions.

FINDINGS

We demonstrate that the phase diagram displays the same qualitative features that have previously been predicted by many-body theory, for mixtures containing flexible polymers under theta solvent conditions. The contraction of the particle two-phase region that we observe, as the correlation length increases beyond some specific value, could be a signature of the weakening of the depletion caused by many-body effects.

Monte Carlo simulation of cylinders with short-range attractions

Cylindrical or rod-like particles are promising materials for the applications of fillers in nanocomposite materials and additives to control rheological properties of colloidal suspensions. Recent advances in particle synthesis allows for cylinders to be manufactured with short-ranged attractions to study the gelation as a function of packing fraction, aspect ratio and attraction strength. In order to aid in the analysis of small-angle scattering experiments of rod-like particles, computer simulation methods were used to model these particles with specialized Monte Carlo algorithms and tabular superquadric potentials. The attractive interaction between neighboring rods increases with the amount of locally-accessible surface area, thus leading to patchy-like interactions. We characterize the clustering and percolation of cylinders as the attractive interaction increases from the homogenous fluid at relatively low attraction strength, for a variety of aspect ratios and packing fractions. Comparisons with the experimental scattering results are also presented, which are in agreement.

Isotropic-Nematic Transition and Demixing Behavior in Binary Mixtures of Hard Spheres and Hard Spherocylinders Confined in a Disordered Porous Medium: Scaled Particle Theory

We develop the scaled particle theory to describe the thermodynamic properties and orientation ordering of a binary mixture of hard spheres (HS) and hard spherocylinders (HSC) confined in a disordered porous medium. Using this theory, the analytical expressions of the free energy, the pressure, and the chemical potentials of HS and HSC have been derived. The improvement of obtained results is considered by introducing the Carnahan-Starling-like and Parsons-Lee-like corrections. Phase diagrams for the isotropic-nematic transition are calculated from the bifurcation analysis of the integral equation for the orientation singlet distribution function and from the conditions of thermodynamic equilibrium. Both the approaches correctly predict the isotropic-nematic transition at low concentrations of hard spheres. However, the thermodynamic approach provides more accurate results and is able to describe the demixing phenomena in the isotropic and nematic phases. The effects of porous medium on the isotropic-nematic phase transition and demixing behavior in a binary HS/HSC mixture are discussed.

Cholesteric ordering predicted using a coarse-grained polymeric model with helical interactions

The understanding of cholesteric liquid crystals at a molecular level is challenging. Limited insights are available to bridge between molecular structures and macroscopic chiral organization. In the present study, we introduce a novel coarse-grained (CG) molecular model, which is represented by flexible chain particles with helical interactions (FCh), to study the liquid crystalline phase behavior of cholesteric molecules such as double strand DNA and α-helix polypeptides using molecular dynamics (MD) simulations. The isotropic-cholesteric phase transitions of FCh molecules were simulated for varying chain flexibilities. A wall confinement was used to break the periodicity along the cholesteric helix director in order to predict the equilibrium cholesteric pitch. The left-handed cholesteric phase was shown for FCh molecules with right-handed chiral interactions, and a spatially inhomogeneous distribution of the nematic order parameter profile was observed in cholesteric phases. It was found that the chain flexibility plays an important role in determining the macroscopic cholesteric pitch and the structure of the cholesteric liquid crystal phase. The simulations provide insight into the relationship between microscopic molecular characteristics and the macroscopic phase behavior.

Structure and Interfacial Tension of a Hard-Rod Fluid in Planar Confinement

The structural properties and interfacial tension of a fluid of rodlike hard-spherocylinder particles in contact with hard structureless flat walls are studied by means of Monte Carlo simulation. The calculated surface tension between the rod fluid and the substrate is characterized by a nonmonotonic trend as a function of the bulk concentration (density) over the range of isotropic bulk concentrations. As suggested by earlier theoretical studies, a surface-ordering scenario is confirmed by our simulations: the local orientational order close to the wall changes from uniaxial to biaxial nematic when the bulk concentration reaches about 85% of the value at the onset of the isotropic-nematic phase transition. The surface ordering coincides with a wetting transition whereby the hard wall is wetted by a nematic film. Accurate values of the fluid-solid surface tension, the adsorption, and the average particle-wall contact distance are reported (over a broad range of densities into the dense nematic region for the first time), which can serve as a useful benchmark for future theoretical and experimental studies on confined rod fluids. The simulation data are supplemented with predictions from second-virial density functional theory, which are in good qualitative agreement with the simulation results.

Density functional theory for dense nematic liquid crystals with steric interactions

The celebrated work of Onsager on hard particle systems, based on the truncated second order virial expansion, is valid at relatively low volume fractions for large aspect ratio particles. While it predicts the isotropic-nematic phase transition, it does not provide a realistic equation of state in that the pressure remains finite for arbitrarily high densities. In this work, we derive a mean field density functional form of the Helmholtz free energy for nematics with hard core repulsion. In addition to predicting the isotropic-nematic transition, the model provides a more realistic equation of state. The energy landscape is much richer, and the orientational probability distribution function in the nematic phase possesses a unique feature-it vanishes on a nonzero measure set in orientation space.