Computer simulations of nematic drops: coupling between drop shape and nematic order.

We perform Monte Carlo computer simulations of nematic drops in equilibrium with their vapor using a Gay-Berne interaction between the rod-like molecules. To generate the drops, we initially perform NPT simulations close to the nematic-vapor coexistence region, allow the system to equilibrate and subsequently induce a sudden volume expansion, followed with NVT simulations. The resultant drops coexist with their vapor and are generally not spherical but elongated, have the rod-like particles tangentially aligned at the surface and an overall nematic orientation along the main axis of the drop. We find that the drop eccentricity increases with increasing molecular elongation, κ. For small κ the nematic texture in the drop is bipolar with two surface defects, or boojums, maximizing their distance along this same axis. For sufficiently high κ, the shape of the drop becomes singular in the vicinity of the defects, and there is a crossover to an almost homogeneous texture; this reflects a transition from a spheroidal to a spindle-like drop.

Liquid-gas separation in colloidal electrolytes.

The liquid-gas transition of an electroneutral mixture of oppositely charged colloids, studied by Monte Carlo simulations, is found in the low-temperature-low-density region. The critical temperature shows a nonmonotonous behavior as a function of the interaction range, kappa(-1), with a maximum at kappasigma approximately 10, implying an island of coexistence in the kappa-rho plane. The system is arranged in such a way that each particle is surrounded by shells of particles with alternating charge. In contrast with the electrolyte primitive model, both neutral and charged clusters are obtained in the vapor phase.

A novel orientation-dependent potential model for prolate mesogens.

An intermolecular potential is introduced for the study of molecular mesogenic fluids. The model combines distinct features of the well-known Gay-Berne and Kihara potentials by incorporating dispersive interactions dependent on the relative pair orientation to a spherocylinder molecular core. Results of a Monte Carlo simulation study focused on the liquid crystal phases exhibited by the model fluid are presented. For the chosen potential parameters, molecular aspect ratio L*=5 and temperatures T*=2, 3, and 5, isotropic, nematic, smectic-A, and hexatic phases are found. The location of the phase boundaries as well as the equation of state of the fluid and further thermodynamical and structural parameters are discussed and contrasted to the Kihara fluid. In comparison to this latter fluid, the model induces the formation of ordered liquid crystalline phases at lower packing fractions and it favors, in particular, the appearance of layered hexatic ordering as a consequence of the greater attractive interaction assigned to the parallel side-to-side molecular pair configurations. The results contribute to the evaluation of the role of specific interaction energies in the mesogenic behavior of prolate molecular liquids in dense environments.

Density functional theory study of the nematic-isotropic transition in an hybrid cell.

We have employed the density functional theory formalism to investigate the nematic-isotropic capillary transitions of a nematogen confined by walls that favor antagonist orientations to the liquid crystal molecules (hybrid cell). We analyze the behavior of the capillary transition as a function of the fluid-substrate interactions and the pore width. In addition to the usual capillary transition between isotropiclike to nematiclike states, we find that this transition can be suppressed when one substrate is wet by the isotropic phase and the other by the nematic phase. Under this condition the system presents interfacelike states which allow us to continuously transform the nematiclike phase to the isotropiclike phase without undergoing a sharp phase transition. Two different mechanisms for the disappearance of the capillary transition are identified. When the director of the nematiclike state is homogeneously planar-anchored with respect to the substrates, the capillary transition ends up in a critical point. This scenario is analogous to the observed in Ising models when confined in slit pores with opposing surface fields which have critical wetting transitions. When the nematiclike state has a linearly distorted director field, the capillary transition continuously transforms in a transition between two nematiclike states.

Parsons-Lee and Monte Carlo study of soft repulsive nematogens.

A general approach based on the Parsons-Lee theory for soft repulsive molecular fluids is employed to investigate the nematogenic behavior of prolate thermotropic liquid crystals over a broad temperature range. The theory is solved for the particular case of the Kihara soft repulsive spherocylinder model, which is mapped into an effective hard core interaction with a temperature-dependent molecular diameter, expected to resemble the average size and shape of the soft molecules at a given temperature. The reduction of the effective molecular diameter with temperature in the Kihara soft repulsive fluid implies implicitly an increase of the elongation of the molecule and induces the stabilization of the nematic phase at smaller effective packing fractions, contrary to what is found for other fluid models. The rationalization of this effect in terms of excluded volume steric arguments is corroborated by the good general agreement between the Parsons-Lee approach and Monte Carlo simulations for the equation of state of the fluid in the vicinity of the isotropic-nematic transition.

Crowding effects in binary mixtures of rod-like and spherical particles.

Crowding effects relevant to the phase stability of binary mixtures of rod-like and spherical particles are investigated by means of Monte Carlo simulations in the isobaric NPT ensemble. The two types of particles are represented, respectively, by freely rotating hard spherocylinders of a moderate aspect ratio (L/sigma = 5) and hard spheres of the same diameter sigma. Molar fractions of spheres ranging xHS = 0.00-0.37 are considered with the aim of characterizing the crowding effects on the liquid crystal phases of the hard spherocylinder fluid induced by the spherical component as depleting agent. We find that the addition of the spherical crowder is beneficial for the stabilization of the layers of the rod-like particles characteristic of the smectic phase. On the contrary, the addition of spheres has a negative impact upon the stability of the nematic phase, where the rod-like particles tend to align collectively parallel to each other. Interestingly, the spheres tend to arrange forming rod-like clusters in the nematic phase and lamellar structures in the smectic phase, which is compensated by the entropy gained by the spherocylinder particles in each phase. The main results are in qualitative agreement with recent experimental and theoretical studies and serve to test the prediction of current equations of state for these types of binary mixtures.

Liquid crystal behavior of the Kihara fluid.

The liquid crystal phases of the Kihara fluid have been studied in computer simulations. The work focuses on the isotropic-nematic-smectic-A triple point region, especially relevant for the understanding of the properties and the design of real mesogens with specific phase diagrams. The Kihara interaction resembles more appropriately than other related models, the shape of elongated polymers and biomolecules, and a closer assertion is provided for the role of the configurational entropy and the dispersive interactions in the behavior of such molecules in dense phases or under macromolecular crowding conditions.

Surface and capillary transitions in an associating binary mixture model.

We investigate the phase diagram of a two-component associating fluid mixture in the presence of selectively adsorbing substrates. The mixture is characterized by a bulk phase diagram that displays peculiar features such as closed loops of immiscibility. The presence of the substrates may interfere with the physical mechanism involved in the appearance of these phase diagrams, leading to an enhanced tendency to phase separate below the lower critical solution point. Three different cases are considered: a planar solid surface in contact with a bulk fluid, while the other two represent two models of porous systems, namely, a slit and an array on infinitely long parallel cylinders. We confirm that surface transitions, as well as capillary transitions for a large surface area to volume ratio, are stabilized in the one-phase region. Applicability of our results to experiments reported in the literature is discussed.

Dipolar origin of the gas-liquid coexistence of the hard-core 1:1 electrolyte model.

We present a systematic study of the effect of the ion pairing on the gas-liquid phase transition of hard-core 1:1 electrolyte models. We study a class of dipolar dimer models that depend on a parameter R(c), the maximum separation between the ions that compose the dimer. This parameter can vary from sigma(+/-) that corresponds to the tightly tethered dipolar dimer model to R(c)--> infinity that corresponds to the Stillinger-Lovett description of the free ion system. The coexistence curve and critical point parameters are obtained as a function of R(c) by grand-canonical Monte Carlo techniques. Our results show that this dependence is smooth but nonmonotonic and converges asymptotically towards the free ion case for relatively small values of R(c). This fact allows us to describe the gas-liquid transition in the free ion model as a transition between two dimerized fluid phases. The role of the unpaired ions can be considered as a perturbation of this picture.

Orientational transitions in a nematic liquid crystal confined by competing surfaces.

The effect of confinement on the orientational structure of a nematic liquid crystal model has been investigated by using a version of density-functional theory. We have focused on the case of a nematic confined by opposing flat surfaces, in slab geometry (slit pore), which favor planar molecular alignment (parallel to the surface) and homeotropic alignment (perpendicular to the surface), respectively. The spatial dependence of the tilt angle of the director with respect to the surface normal has been studied, as well as the tensorial order parameter describing the molecular order around the director. For a pore of given width, we find that, for weak surface fields, the alignment of the nematic director is perpendicular to the surface in a region next to the surface favoring homeotropic alignment, and parallel along the rest of the pore, with a sharp interface separating these regions (S phase). For strong surface fields, the director is distorted uniformly, the tilt angle exhibiting a linear dependence on the distance normal to the surface (L phase). Our calculations reveal the existence of a first-order transition between the two director configurations, which is driven by changes in the surface field strength, and also by changes in the pore width. In the latter case the transition occurs, for a given surface field, between the S phase for narrow pores and the L phase for wider pores. A link between the L-S transition and the anchoring transition observed for the semi-infinite case is proposed.