Phase transitions and topological defects in discotic liquid crystal droplets with planar anchoring: a Monte Carlo simulation study.

In this work we present the results of Monte Carlo (MC) simulations at the isothermal-isobaric ensemble for a discotic liquid crystal (DLC) droplet whose surface promotes edge-on (planar) anchoring. For a given pressure, we simulate an annealing process that enables observation of phase transitions within the spherical droplet. In particular, we report a first order isotropic-nematic transition as well as a nematic-columnar transition at the center of the droplet. We found the appearance of topological defects consisting of two disclination lines with ends at the surface of the sphere. We also observed that both transitions, isotropic-nematic and nematic-columnar, occur at lower temperatures as compared to the unconfined system.

Structural properties and ring defect formation in discotic liquid crystal nanodroplets.

In this work, we performedNpTMonte Carlo simulations of a Gay-Berne discotic liquid crystal confined in a spherical droplet under face-on anchoring and fixed pressure. We find that, in contrast to the unbounded system, a plot of the order parameter as function of temperature does not show a clear evidence of a first-order isotropic-nematic transition. We also find that the impossibility of simultaneously satisfy the uniform director field requirement of a nematic phase with the radial boundary conditions, results in the appearance of a ring disclination line as a stress release mechanism in the interior of the droplet. Under further cooling, a columnar phase appears at the center of the droplet.

Effect of the anchoring strength on the phase behaviour of discotic liquid crystals under face-on confinement.

In this study we have performed molecular dynamics simulations to study a Gay-Berne discotic fluid confined in a slab geometry for a fixed confinement length. Four different anchoring strengths with a homeotropic (face-on) configuration were studied. We found that changing the anchoring strength changes the normal component of the stress tensor, which in turn changes the density of the system's bulk. This phenomenon leads to a shift in the isotropic-nematic transition temperature. We observe that the temperature regions where the nematic phase is present diminishes as the anchoring strength increases. The anchoring strength also affects the nematic-columnar coexistence temperature-region: it spans over more temperatures at higher anchoring strengths.

Influence of anchoring in the phase behaviour of discotic liquid crystals.

Molecular dynamics simulations were performed for a Gay-Berne discotic fluid confined in a slab geometry for two different anchorings: homeotropic (face-on) and planar (edge-on), and for two different confinement lengths. Our results show that the behaviour of the more confined system in the temperature region of the isotropic-nematic transition is critically influenced by the presence of the walls: the growth of the solid-liquid crystal interface spans over the entire width of the cell, and the character of the transition is changed from first order to continuous. For all the confined systems studied, we observe a higher nematic-columnar transition temperature and a smaller nematic phase region in the phase diagram, as compared with the behaviour of the infinite system.

Induced stabilization of columnar phases in binary mixtures of discotic liquid crystals.

Three discotic liquid-crystalline binary mixtures, characterized by their extent of bidispersity in molecular thickness, were investigated with molecular dynamics simulations. Each equimolar mixture contained A-type (thin) and B-type (thick) discogens. The temperature-dependence of the orientational order parameter reveals that A-type liquid samples produce ordered phases more readily, with the (hexagonal) columnar phase being the most structured variant. Moderately and strongly bidisperse mixtures produce globally-segregated samples for temperatures corresponding to ordered phases; the weakly bidisperse mixture displays microheterogeneities. Ordered phases in the B-type liquid are induced partially by the presence of the A-type fluid. In the moderately bidisperse mixture, order is induced through orientational frustration: a mixed prenematic-like phase precedes global segregation to yield nematic and columnar mesophases upon further cooling. In the strongly bidisperse mixture, order is induced less efficiently through a paranematic-like mechanism: a highly-ordered A-type fluid imparts order to B-type discogens found at the interface of a fully-segregated sample. This ordering effect permeates into the disordered B-type domain until nematic and columnar phases emerge upon further cooling. At sufficiently low temperatures, all samples investigated exhibit the (hexagonal) columnar mesophase.

Room temperature ionic liquids: A simple model. Effect of chain length and size of intermolecular potential on critical temperature.

A model of a room temperature ionic liquid can be represented as an ion attached to an aliphatic chain mixed with a counter ion. The simple model used in this work is based on a short rigid tangent square well chain with an ion, represented by a hard sphere interacting with a Yukawa potential at the head of the chain, mixed with a counter ion represented as well by a hard sphere interacting with a Yukawa potential of the opposite sign. The length of the chain and the depth of the intermolecular forces are investigated in order to understand which of these factors are responsible for the lowering of the critical temperature. It is the large difference between the ionic and the dispersion potentials which explains this lowering of the critical temperature. Calculation of liquid-vapor equilibrium orthobaric curves is used to estimate the critical points of the model. Vapor pressures are used to obtain an estimate of the triple point of the different models in order to calculate the span of temperatures where they remain a liquid. Surface tensions and interfacial thicknesses are also reported.

Fluid-solid coexistence from two-phase simulations: binary colloidal mixtures and square well systems.

Molecular dynamics simulations are performed to clarify the reasons for the disagreement found in a previous publication [G. A. Chapela, F. del Río, and J. Alejandre, J. Chem. Phys. 138(5), 054507 (2013)] regarding the metastability of liquid-vapor coexistence on equimolar charged binary mixtures of fluids interacting with a soft Yukawa potential with κσ = 6. The fluid-solid separation obtained with the two-phase simulation method is found to be in agreement with previous works based on free energy calculations [A. Fortini, A.-P. Hynninen, and M. Dijkstra, J. Chem. Phys. 125, 094502 (2006)] only when the CsCl structure of the solid is used. It is shown that when pressure is increased at constant temperature, the solids are amorphous having different structures, densities, and the diagonal components of the pressure tensor are not equal. A stable low density fluid-solid phase separation is not observed for temperatures above the liquid-vapor critical point. In addition, Monte Carlo and discontinuous molecular dynamics simulations are performed on the square well model of range 1.15σ. A stable fluid-solid transition is observed above the vapor-liquid critical temperature only when the solid has a face centered cubic crystalline structure.

Disclination lines at homogeneous and heterogeneous colloids immersed in a chiral liquid crystal.

In the present work we perform Monte Carlo simulations in the isothermal-isobaric ensemble to study defect topologies formed in a cholesteric liquid crystal due to the presence of a spherical colloidal particle. Topological defects arise because of the competition between anchoring at the colloidal surface and the local director. We consider homogeneous colloids with either local homeotropic or planar anchoring to validate our model by comparison with earlier lattice Boltzmann studies. Furthermore, we perform simulations of a colloid in a twisted nematic cell and discuss the difference between induced and intrinsic chirality on the formation of topological defects. We present a simple geometrical argument capable of describing the complex three-dimensional topology of disclination lines evolving near the surface of the colloid. The presence of a Janus colloid in a cholesteric host fluid reveals a rich variety of defect structures. Using the Frank free energy we analyze these defects quantitatively indicating a preferred orientation of the Janus colloid relative to the cholesteric helix.

Phase equilibria, fluid structure, and diffusivity of a discotic liquid crystal.

Molecular Dynamics simulations were performed for the Gay-Berne discotic fluid parameterized by GB(0.345, 0.2, 1.0, 2.0). The volumetric phase diagram exhibits isotropic (IL), nematic (ND), and two columnar phases characterized by radial distribution functions: the transversal fluid structure varies between a hexagonal columnar (CD) phase (at higher temperatures and pressures) and a rectangular columnar (CO) phase (at lower temperatures and pressures). The slab-wise analysis of fluid dynamics suggests the formation of grain-boundary defects in the CO phase. Longitudinal fluid structure is highly periodic with narrow peaks for the CO phase, suggestive of a near-crystalline (yet diffusive) system, but is only short-ranged for the CD phase. The IL phase does not exhibit anisotropic diffusion. Transversal diffusion is more favorable in the ND phase at all times, but only favorable at short times for the columnar phases. In the columnar phases, a crossover occurs where longitudinal diffusion is favored over transversal diffusion at intermediate-to-long timescales. The anomalous diffusivity is pronounced in both columnar phases, with three identifiable contributions: (a) the rattling of discogens within a transient "interdigitation" cage, (b) the hopping of discogens across columns, and (c) the drifting motion of discogens along the orientation of the director.

Analytical equation of state with three-body forces: application to noble gases.

We developed an explicit equation of state (EOS) for small non polar molecules by means of an effective two-body potential. The average effect of three-body forces was incorporated as a perturbation, which results in rescaled values for the parameters of the two-body potential. These values replace the original ones in the EOS corresponding to the two-body interaction. We applied this procedure to the heavier noble gases and used a modified Kihara function with an effective Axilrod-Teller-Muto (ATM) term to represent the two- and three-body forces. We also performed molecular dynamics simulations with two- and three-body forces. There was good agreement between predicted, simulated, and experimental thermodynamic properties of neon, argon, krypton, and xenon, up to twice the critical density and up to five times the critical temperature. In order to achieve 1% accuracy of the pressure at liquid densities, the EOS must incorporate the effect of ATM forces. The ATM factor in the rescaled two-body energy is most important at temperatures around and lower than the critical one. Nonetheless, the rescaling of two-body diameter cannot be neglected at liquid-like densities even at high temperature. This methodology can be extended straightforwardly to deal with other two- and three-body potentials. It could also be used for other nonpolar substances where a spherical two-body potential is still a reasonable coarse-grain approximation.

Effect of flexibility on liquid-vapor coexistence and surface properties of tangent linear vibrating square well chains in two and three dimensions.

The effect of flexibility on liquid-vapor and interfacial properties of tangent linear vibrating square well chains is studied. Surface tension, orthobaric densities, vapor pressures, and interfacial thicknesses are reported and analyzed using corresponding states principles. Discontinuous molecular dynamics simulations in two and three dimensions are performed on rigid tangent linear vibrating square well chains of different lengths. In the case of two dimensions, simulation results of completely flexible tangent linear vibrating square well chains are also reported. Properties are calculated for chains of 2-12 monomers. Rigidity is controlled by trapping the first and last monomer in the chain in a vibrating well at half of the distance of the whole chain. Critical property values are reported as obtained from orthobaric densities, surface tensions, and vapor pressures. For the fully flexible chains, the critical temperatures increase with chain length but the effect saturates. In contrast, the critical temperatures increase for the rigid chains until no more critical point is found.

Fractionation in Gay-Berne liquid crystal mixtures.

We present a constant-pressure molecular dynamics simulation study of the phase behavior of binary (50:50) Gay-Berne liquid crystal mixtures consisting of elongated particles with different lengths (LA>LB) and equal diameters. We focus on systems at dense liquid-state conditions. Considering three mixtures characterized by different values of LA(B) and different length ratios q=LB/LA<1, we find complex fluid-fluid phase behavior resulting from the interplay between nematic, smectic-A-type, or smectic-B-type orientational ordering, on the one hand, and demixing into two phases of different composition (fractionation), on the other hand. The driving "forces" of demixing transitions are the temperature and the length ratio. Indeed, in the system characterized by the largest value of q (q=0.86) orientational order occurs already in mixed states, whereas full fractionation is found at q=0.71. The two resulting states are either of type smectic-B-nematic (intermediate temperatures) or smectic-B-smectic-B (low temperatures). In the intermediate case q=0.80 we observe a stepwise ordering and demixing behavior on cooling the system from high temperatures. Moreover, our results show that the stability range of (partially) nematic structures in mixtures of sufficiently small q can be significantly larger than in the pure counterparts, in qualitative agreement with experimental observations.

Phase and interfacial behavior of partially miscible symmetric Lennard-Jones binary mixtures.

We have carried out extensive equilibrium molecular-dynamics simulations to study quantitatively the topology of the temperature versus density phase diagrams and related interfacial phenomena in a partially miscible symmetric Lennard-Jones binary mixture. The topological features are studied as a function of miscibility parameter, alpha = epsilonAB/epsilonAA. Here epsilonAA = epsilonBB and epsilonAB stand for the parameters related to the attractive part of the intermolecular interactions for similar and dissimilar particles, respectively. When the miscibility varies in the range 0 < alpha < 1, a continuous critical line of consolute points Tcons(rho)--critical demixing transition line--appears. This line intersects the liquid-vapor coexistence curve at different positions depending on the values of alpha, yielding mainly three different topologies for the phase diagrams. These results are in qualitative agreement to those found previously for square-well and hard-core Yukawa binary mixtures. The main contributions of the present paper are (i) a quantitative analysis of the phase behavior and (ii) a detailed study of the liquid-liquid interfacial and liquid-vapor surface tensions, as function of temperature and miscibility as well as its relationship to the topological features of the phase diagrams.

Liquid-vapor equilibrium and surface tension of nonconformal molecular fluids.

Molecular dynamics simulations at constant temperature have been performed on the liquid-vapor interface for fluids characterized by a recently introduced three-parameter potential. This potential is a modification of the well-known spherical Kihara interaction and is termed approximate nonconformal (ANC). It has been used successfully to describe many real molecules in the gaseous phase. Besides the usual molecular energy and size, the ANC potential introduces a third parameter s, called softness, to measure the form of the potential profile. Study of these systems shows that their critical and interfacial properties follow very closely those of four selected substances: argon, methane, propane, and hexane. Deviations of the properties predicted from the experimental values are analyzed and their probable causes are determined. The critical properties of ANC fluids and their dependence on s are also obtained via first-order perturbation theory in the form of an augmented van der Waals model. Analysis of the results shows that ANC potential functions can be used as reliable effective interactions for real dense fluids.

Wetting phenomenon in the liquid-vapor phase coexistence of a partially miscible Lennard-Jones binary mixture.

We have carried out extensive equilibrium molecular dynamics simulations to study the structure and the interfacial properties in the liquid-vapor phase coexistence of partially miscible binary Lennard-Jones mixtures. By analyzing the structural properties as a function of the miscibility parameter, alpha, we found that at relatively low temperatures the system separates forming a liquid A-liquid B interface in coexistence with the vapor phase. At higher temperatures and, 0< alpha < or =0.5 , we found a temperature range, T*w (alpha) < or =T*< T*Cons (alpha) , where the liquid phases are wet by the vapor phase. Here, T*w (alpha) represents the wetting transition temperature and T*Cons (alpha) is the consolute temperature of the mixture. However, for 0.5< alpha <1 , no wetting phenomenon occurs. For the particular value, alpha=0.25 , we analyzed quantitatively the T* versus rho* , and P* versus T* phase diagrams and found, T*c approximately 1.25 , and T*Cons approximately 1.25 . We also studied quantitatively, as a function of temperature, the surface tension and the adsorption of molecules at the liquid-liquid interface. It was found that the adsorption shows a jump from a finite negative value up to minus infinity, when the vapor wets the liquid phases, suggesting that the wetting transition is of first order. The calculated phase diagram, together with the wetting phenomenon, strongly suggests the existence of a tricritical point. These results agree well with some experiments carried out in fluid binary mixtures.

Metastable liquid lamellar structures in binary and ternary mixtures of Lennard-Jones fluids.

We have carried out extensive equilibrium molecular dynamics simulations to investigate the liquid-vapor coexistence in partially miscible binary and ternary mixtures of Lennard-Jones fluids. We have studied in detail the time evolution of the density profiles and the interfacial properties in a temperature region of the phase diagram where the condensed phase is demixed. The composition of the mixtures is fixed, 50% for the binary mixture and 33.33% for the ternary mixture. The results of the simulations clearly indicate that in the range of temperatures 78

Phase diagram of symmetric binary fluid mixtures: first-order or second-order demixing.

Binary fluid mixtures of 1:1 concentration can demix in a phase transition of first order or of second order. We analyze the two scenarios in density-concentration space and relate them to the structure of the line at which the demixing coexistence surface cuts the liquid-vapor coexistence surface. These scenarios help us to decide between first and second order for a model of a symmetric Lennard-Jones mixture. An optimized reference hypernetted chain integral equation method is employed for calculating the correlation functions and from there the pressure and chemical potentials. We conclude that demixing of a 1:1 mixture is of first order in the whole range of parameters that we have investigated. We did not find a critical point in the 1:1 concentration plane.