A Monte Carlo simulation study of binary mixtures of charged hard spherocylinders and charged hard spheres

Phase behaviour of mixtures of charged soft disks and spheres

We have investigated the phase behaviour of mixtures of soft disks (Gay-Berne oblate ellipsoids, GB) and soft spheres (Lennard-Jones, LJ) with opposite charge as a model of ionic liquid crystals and colloidal suspensions. We have used constant volume Molecular Dynamics simulations and fixed the stoichiometry of the mixture in order to have electroneutrality; three systems have been selected GB : LJ = 1 : 2, GB : LJ = 1 : 1 and GB : LJ = 2 : 1. For each system we have selected three values of the scaled point charge q* of the GB particles, namely 0.5, 1.0 and 2.0 (and a corresponding negative scaled charge of the LJ particles that depends on the stoichiometric ratio). We have found a very rich mesomorphism with the formation, as a function of the scaled temperature, of the isotropic phase, the discotic nematic phase, the hexagonal columnar phase and crystal phases. While the structure of the high temperature phases was similar in all systems, the hexagonal columnar phases exhibited a highly variable morphology depending on the scaled charge and stoichiometry. On the one hand, GB : LJ = 1 : 2 systems form lamellar structures, akin to smectic phases, with an alternation of layers of disks (exhibiting an hexagonal columnar phase) and layers of LJ particles (in the isotropic phase). On the other hand, for the 2 : 1 stoichiometry we observe the formation of a frustrated hexagonal columnar phase with an alternating tilt direction of the molecular axis. We rationalize these findings based on the structure of the neutral ion pair dominating the behaviour at low temperature and high charge.

Theoretical equations of state for a charged fluid

In this article, we present a molecular thermodynamic study of a system of N particles contained within a volume V and interacting via a hard-core pair potential with an attractive interaction according to the Wolf model for charged systems. This variable-range potential is characterized by three parameters: the repulsive hard-core diameter σ, the energy-well depth ϵ, and the inverse range α; a fourth parameter of the model is a cut-off distance x_c that depends on α according to the relation x_c = 2/α. Two equations of state (EOSs) are presented and derived from thermodynamic perturbation theory and Monte Carlo (MC) simulation data. The first EOS is given by the standard Zwanzig's high-temperature expansion of the Helmholtz free energy, where the first three perturbation terms a₁, a₂, and a₃ were obtained from MC simulations in the canonical ensemble (NVT) and parameterized as functions of α and the reduced density of particles ρ^* = Nσ³/V. The second EOS was obtained from the discrete perturbation theory applied to a discrete representation of the Wolf potential. Results for pressures, internal energies, and isochoric heat capacities are compared to the MC computer simulation data of the Wolf system, including vapor-liquid coexistence curves, for different values of α. Overall, both EOSs give a very good representation of the thermodynamic properties of the Wolf fluid when 0.3 ≤ α ≤ 1.0 and 0.05 ≤ ρ^* ≤ 0.8. Since the Yukawa fluid can reproduce information of screened ionic interactions, we discuss the equivalence between the Wolf and Yukawa fluids in the context of equivalent systems in liquid theory.

A Monte Carlo study of the electrical double layer of a shape-asymmetric electrolyte around a spherical colloid

In this paper, we present a Monte Carlo simulation study on the structure of the electrical double layer around a spherical colloid surrounded by a binary electrolyte composed of spherical and non-spherical ions. Results are provided for the radial distribution functions between the colloid and ions, the orientation correlations between the colloid and non-spherical particles, and the integrated charge. Work is reported mainly for non-spherical particles modeled as spherocylinders, although a particular comparison is made between spherocylindrical particles and dimers. For the conditions investigated here, spherocylinders and dimers produce essentially the same structural information. Additionally, it is shown that spherocylinders mostly orient tangentially to the colloid at its surface; this preferred orientation disappears for larger distances. We also evidence that, near the colloid, the integrated charge attenuates monotonically when the macroparticle is highly charged, whereas for intermediate and low charged states of the colloid, the integrated charge can display charge reversal, overcharging, or both, with magnitudes that are sensitive to the salt concentration and to the localization of charge inside the spherocylinders.

Generalized Onsager theory for strongly anisometric patchy colloids

The implications of soft "patchy" interactions on the orientational disorder-order transition of strongly elongated colloidal rods and flat disks is studied within a simple Onsager-van der Waals density functional theory. The theory provides a generic framework for studying the liquid crystal phase behaviour of highly anisometric cylindrical colloids which carry a distinct geometrical pattern of repulsive or attractive soft interactions localized on the particle surface. In this paper, we apply our theory to the case of charged rods and disks for which the local electrostatic interactions can be described by a screened-Coulomb potential. We consider infinitely thin rod like cylinders with a uniform line charge and infinitely thin discotic cylinders with several distinctly different surface charge patterns. Irrespective of the backbone shape the isotropic-nematic phase diagrams of charged colloids feature a generic destabilization of nematic order at low ionic strength, a dramatic narrowing of the biphasic density region, and a reentrant phenomenon upon reducing the electrostatic screening. The low screening regime is characterized by a complete suppression of nematic order in favor of positionally ordered liquid crystal phases.

Monte Carlo computer simulation of sedimentation of charged hard spherocylinders

In this article we present a NVT Monte Carlo computer simulation study of sedimentation of an electroneutral mixture of oppositely charged hard spherocylinders (CHSC) with aspect ratio L/σ = 5, where L and σ are the length and diameter of the cylinder and hemispherical caps, respectively, for each particle. This system is an extension of the restricted primitive model for spherical particles, where L/σ = 0, and it is assumed that the ions are immersed in an structureless solvent, i.e., a continuum with dielectric constant D. The system consisted of N = 2000 particles and the Wolf method was implemented to handle the coulombic interactions of the inhomogeneous system. Results are presented for different values of the strength ratio between the gravitational and electrostatic interactions, Γ = (mgσ)/(e(2)/Dσ), where m is the mass per particle, e is the electron's charge and g is the gravitational acceleration value. A semi-infinite simulation cell was used with dimensions Lx ≈ Ly and Lz = 5Lx, where Lx, Ly, and Lz are the box dimensions in Cartesian coordinates, and the gravitational force acts along the z-direction. Sedimentation effects were studied by looking at every layer formed by the CHSC along the gravitational field. By increasing Γ, particles tend to get more packed at each layer and to arrange in local domains with an orientational ordering along two perpendicular axis, a feature not observed in the uncharged system with the same hard-body geometry. This type of arrangement, known as tetratic phase, has been observed in two-dimensional systems of hard-rectangles and rounded hard-squares. In this way, the coupling of gravitational and electric interactions in the CHSC system induces the arrangement of particles in layers, with the formation of quasi-two dimensional tetratic phases near the surface.

Anomalous columnar order of charged colloidal platelets

Monte Carlo computer simulations are carried out for a model system of like-charged colloidal platelets in the isothermal-isobaric ensemble (NpT). The aim is to elucidate the role of electrostatic interactions on the structure of synthetic clay systems at high particle densities. Short-range repulsions between particles are described by a suitable hard-core model representing a discotic particle. This potential is supplemented with an electrostatic potential based on a Yukawa model for the screened Coulombic potential between infinitely thin disklike macro-ions. The particle aspect-ratio and electrostatic parameters were chosen to mimic an aqueous dispersion of thin, like-charged, rigid colloidal platelets at finite salt concentration. An examination of the fluid phase diagram reveals a marked shift in the isotropic-nematic transition compared to the hard cut-sphere reference system. Several statistical functions, such as the pair correlation function for the center-of-mass coordinates and structure factor, are obtained to characterize the structural organization of the platelets phases. At low salinity and high osmotic pressure we observe anomalous hexagonal columnar structures characterized by interpenetrating columns with a typical intercolumnar distance corresponding to about half of that of a regular columnar phase. Increasing the ionic strength leads to the formation of glassy, disordered structures consisting of compact clusters of platelets stacked into finite-sized columns. These so-called "nematic columnar" structures have been recently observed in systems of charge-stabilized gibbsite platelets. Our findings are corroborated by an analysis of the static structure factor from a simple density functional theory.