Phase behavior and interfacial properties of nonadditive mixtures of Onsager rods

Within a second virial theory, we study bulk phase diagrams as well as the free planar isotropic-nematic interface of binary mixtures of nonadditive thin and thick hard rods. For species of the same type, the excluded volume is determined only by the dimensions of the particles, whereas for dissimilar ones it is taken to be larger or smaller than that, giving rise to a nonadditivity that can be positive or negative. We argue that such a nonadditivity can result from modeling of soft interactions as effective hard-core interactions. The nonadditivity enhances or reduces the fractionation at isotropic-nematic (IN) coexistence and may induce or suppress a demixing of the high-density nematic phase into two nematic phases of different composition (N(1) and N(2)), depending on whether the nonadditivity is positive or negative. The interfacial tension between coexisting isotropic and nematic phases shows an increase with increasing fractionation at the IN interface, and complete wetting of the IN(2) interface by the N(1) phase upon approach of the triple-point coexistence. In all explored cases bulk and interfacial properties of the nonadditive mixtures exhibit a striking and quite unexpected similarity with the properties of additive mixtures of different diameter ratio.     

Nematic-nematic demixing in polydisperse thermotropic liquid crystals

We consider the effects of polydispersity on isotropic-nematic phase equilibria in thermotropic liquid crystals, using a Maier-Saupe theory [Z. Naturforsch. A 13A, 564 (1958)] with factorized interactions. A sufficient spread (approximately 50%) in the interaction strengths of the particles leads to phase separation into two or more nematic phases, which can in addition coexist with an isotropic phase. The isotropic-nematic coexistence region widens dramatically as polydispersity is increased, leading to reentrant isotropic-nematic phase separation in some regions of the phase diagram. We show that similar phenomena will occur also for nonfactorized interactions as long as the interaction strength between any two particle species is lower than the mean of the intraspecies interactions.     

Isotropic-nematic transition of hard rods immersed in random sphere matrices

Using replica density functional theory and Monte Carlo computer simulations we investigate a system of annealed hard spherocylinders adsorbed in a matrix of quenched hard spheres. Theoretical predictions for the partition coefficient, defined as the ratio of density of rods in the matrix and that in a reservoir, agree well with simulation results. Theory predicts the isotropic-nematic transition to remain first order upon increasing sphere packing fraction, and to shift towards lower rod densities. This scenario is consistent with our simulation results that clearly show a jump in the nematic order parameter upon increasing the rod density at constant matrix packing fraction, corresponding to the isotropic-nematic transition, even for sphere matrix packing fractions < or approximately equal to 0.3.     

Periodic orientational motions of rigid liquid-crystalline polymers in shear flow

The collective periodic motions of liquid-crystalline polymers in a nematic phase in shear flow have, for the first time, been simulated at the particle level by Brownian dynamics simulations. A wide range of parameter space has been scanned by varying the aspect ratio L/D between 10 and 60 at three different scaled volume fractions Lphi/D and an extensive series of shear rates. The influence of the start configuration of the box on the final motion has also been studied. Depending on these parameters, the motion of the director is either characterized as tumbling, kayaking, log-rolling, wagging, or flow-aligning. The periods of kayaking and wagging motions are given by T=4.2(Lphi/D)gamma(-1) for high aspect ratios. Our simulation results are in agreement with theoretical predictions and recent shear experiments on fd viruses in solution. These calculations of elongated rigid rods have become feasible with a newly developed event-driven Brownian dynamics algorithm.     

Computation of the nonhomogeneous equilibrium states of a rigid-rod solution

The nonhomogeneous equilibrium phase behavior of a solution of rigid rods is analyzed for a periodic one-dimensional system. Stable and unstable equilibrium solutions for the distribution function are computed as extrema of the free energy of the system expressed by the nonhomogeneous generalization of Onsager's [Ann. N.Y. Acad. Sci. 51, 627 (1949)] theory, which models interaction between rods on the scale of a single rod length. Biaxial equilibrium solutions are computed in a periodic system by discretizing the Euler-Lagrange nonlinear integral equation by the finite-element method and using Newton's method to solve the resulting set of nonlinear equations. Stable states for isotropic-nematic coexistence are computed in a periodic system rather than the semi-infinite system used in previous calculations. The density and order parameter profiles evolve monotically from the isotropic phase to the nematic phase. Unstable, nonhomogeneous, equilibrium states are also computed for concentrations of rods that exceed the value for spinodal decomposition. These nonhomogeneous states are characterized by combinations of bend, twist, and splay distortions in physical space and correspond to unstable attractors in the dynamic process of isotropic-nematic spinodal decomposition. For large systems, the nonhomogeneous states develop wide, bulklike nematic regions separated by thin regions with sharp gradients in orientation. The free energy formulation was also used to compute the accurate neutral stability curve; this curve shows the limits of applicability of the low-wave-number approximations frequently used in the study of spinodal decomposition.     

Time-resolved synchrotron X-ray scattering of the crystallization of a soft hexagonal columnar crystal

We have recently shown (Langmuir 2000, 16, 5846) that a soft surfactant hexagonal phase exhibits, above a critical shear rate, a shear-melting transition, where a two-dimensional polycrystalline texture is converted into a liquid of rods aligned along the flow. Moreover, after abrupt cessation of high shear, a monocrystal-like structure is obtained. In this paper, we investigate the kinetics of crystallization from the shear-melted phase by means of time-resolved synchrotron X-ray scattering experiments. After a latency time, the sample is found to continuously evolve from a liquid of rods to a monocrystal of rods. We show that the crystallization results from a coupling between bulk crystallization and interfacial crystallization, due to a specific anchoring of the liquid crystal at the walls of the shear cell. A detailed analysis of the experimental data allows the evaluation of the propagation of the interfacial front and the time evolution of the size of the crystallites. We find that the size scales as t(v) with v approximately equal to 0.3, in fair agreement with recent theories and simulations on grain rotation-induced grain growth of columnar polycrystalline structures.     

Dichroism in dye-doped colloidal liquid crystals

Nematic liquid crystals were obtained in sterically stabilized suspensions of rodlike particles of sepiolite clay, with an average length up to 900 nm and aspect ratio up to 40. In agreement with computer simulations for hard spherocylinders, the isotropic-nematic transition shifted to lower volume fractions with increasing aspect ratio. However, the coexistence gap was broadened noticeably due to particle polydispersity. The sepiolite crystal structure includes channels filled with zeolitic water, which can be replaced by indigo dye molecules. The indigo molecules are constrained inside the zeolitic channels to be aligned along the long axes of the rods. As a result, the colloidal nematic phase showed a marked dichroism, with an order parameter up to 0.5 for magnetically aligned samples, similar to typical values for dye-doped thermotropic liquid crystals.     

Isotropic-nematic interfacial tension of hard and soft rods: application of advanced grand canonical biased-sampling techniques

Coexistence between the isotropic and the nematic phase in suspensions of rods is studied using grand canonical Monte Carlo simulations with a bias on the nematic order parameter. The biasing scheme makes it possible to estimate the interfacial tension gamma(IN) in systems of hard and soft rods. For hard rods with LD=15, we obtain gammaIN approximately 1.4kBT/L2, with L the rod length, D the rod diameter, T the temperature, and kB the Boltzmann constant. This estimate is in good agreement with theoretical predictions, and the order of magnitude is consistent with experiments.     

Stability limits in binary fluids mixtures

The stability limits in binary fluid mixtures are investigated on the basis of the global phase diagram approach employing a model for the attracting hard-sphere fluid. In addition to the diffusion spinodals the mechanical spinodals are included. As a result one finds topologically different types of the diffusion spinodals while only one shape exists for the mechanical spinodals which are present in the region of liquid-vapor equilibria only. The diffusion spinodals represent the underlying properties of the phase behavior. The types of stable phase behavior therefore resemble that of the spinodal behavior. The different shapes of the spinodals can be important for nonequilibrium processes in nature and technology.     

Critical behavior of self-assembled rigid rods on triangular and honeycomb lattices

Using Monte Carlo simulations and finite-size scaling analysis, the critical behavior of self-assembled rigid rods on triangular and honeycomb lattices at intermediate density has been studied. The system is composed of monomers with two attractive (sticky) poles that, by decreasing temperature or increasing density, polymerize reversibly into chains with three allowed directions and, at the same time, undergo a continuous isotropic-nematic (IN) transition. The determination of the critical exponents, along with the behavior of Binder cumulants, indicate that the IN transition belongs to the q=1 Potts universality class.     

Competition between liquid-crystalline ordering and glassy freezing in melts of semiflexible polymers: A monte carlo simulation

We present results of a Monte Carlo simulation of dense melts of semi-flexible polymers using the bond-fluctuation model. The chosen Hamiltonian increases the chain stiffness upon cooling which in turn leads to glass-transition like freezing of the polymer mobility. Employing an efficient simulation algorithm, which is able to equilibrate the simulated systems to lower temperature than the Rouse-type algorithm showing the glassy freezing, we are able to observe an isotropic-nematic phase transition. This transition lies above the glass transition temperature one would extrapolate from the observed freezing behavior.     

Bifurcational analysis of the isotropic-nematic phase transition of rigid rod polymers subjected to biaxial stretching flow

A bifurcational analysis is performed on Doi's equation of nematodynamics that describes the non-equlibrium isotropic-nematic phase transition of rigid rod polymers in the presence of steady biaxial stretching flow. The symmetry of the flow and of the governing order parameter equations are shown to be the source of a rich bifurcation, symmetry breaking, and multistability behavior involving two physically equivalent biaxial nematic phases, one uniaxial nematic phase and one uniaxial paranematic phase. According to the relative intensity of the nematic ordering field and stretching rate, the uniaxial isotropic-biaxial nematic transition may be continuous (2nd order), discontinuous (1st order), or it may exhibit a tricritical non-equilibrium phase transition point. The solutions to the Doi equations of nematodynamics are found to be consistent with those of Khokhlov and Semenov [Macromolecules 15 , 1272 (1982)], which are based on a version of the Onsager theory of isotropic-nematic phase transitions. The present simulations provide a useful guide for orientation control in biaxial stretching flows.     

Applications of polymer rism theory to blends

Polymer RISM (reference interaction site model) theory is a theory of polymer systems in the liquid phase in which account for chemical realism can be made. Results are reported here of phase diagrams (spinodals) for blend systems calculated by means of this theory, using the mean spherical approximation as a closure. The systems investigated are an isotopic blend, a set of homopolymer/copolymer blends, and a model blend containing specific interactions.     

Theoretical analysis of polydispersity in the nematic phase of self-assembled semiflexible chains

In recent simulations of semiflexible equilibrium polymers [X. Lu and J. T. Kindt, J. Chem. Phys. 120, 10328 (2004)] a roughly biexponential distribution of chain lengths was observed in the nematic phase. In this study, we show that a theory representing the nematic phase as an equilibrium mixture of randomly oriented chains below a critical length Lc and chains above Lc obeying a Gaussian orientational distribution reproduces this distribution qualitatively. The agreement between predicted and simulated phase boundaries is improved, with a narrower isotropic-nematic coexistence region, compared with a monoexponential chain distribution. We find, furthermore, that the critical length Lc scales as the inverse of monomer concentration, irrespective of the bond strength of the chains, but that at the phase boundary, the fraction of monomers found in disordered chains peaks at a certain bond strength.     

Glassiness of thermotropic liquid crystals across the isotropic-nematic transition

The orientational dynamics of thermotropic liquid crystals across the isotropic-nematic phase transition have traditionally been investigated at long times or low frequencies using frequency domain measurements. The situation has now changed significantly with the recent report of a series of interesting transient optical Kerr effect (OKE) experiments that probed orientational relaxation of a number of calamitic liquid crystals (which consist of rod-like molecules) directly in the time domain, over a wide time window ranging from subpicoseconds to tens of microseconds. The most intriguing revelation is that the decay of the OKE signal at short to intermediate times (from a few tens of picoseconds to several hundred nanoseconds) follows multiple temporal power laws. Another remarkable feature that has emerged from these OKE measurements is the similarity in the orientational relaxation behavior between the isotropic phase of calamitic liquid crystals near the isotropic-nematic transition and supercooled molecular liquids, notwithstanding their largely different macroscopic states. In this article, we present an overview of the understanding that has emerged from recent computational and theoretical studies of calamitic liquid crystals across the isotropic-nematic transition. Topics discussed include (a) single-particle as well as collective orientational dynamics at a short-to-intermediate time window, (b) heterogeneous dynamics in orientational degrees of freedom diagnosed by a non-Gaussian parameter, (c) fragility, and (d) temperature-dependent exploration of underlying energy landscapes as calamitic liquid crystals settle into increasingly ordered mesophases upon cooling from the high-temperature isotropic phase. A comparison of our results with those of supercooled molecular liquids reveals an array of analogous features in these two important classes of soft matter systems. We further find that the onset of growth of the orientational order in the parent nematic phase induces translational order, resulting in smectic-like layers in the potential energy minima of calamitic systems if the parent nematic phase is sandwiched between the high-temperature isotropic phase and the low-temperature smectic phase. We discuss implications of this startling observation. We also discuss recent results on the orientational dynamics of discotic liquid crystals that are found to be rather similar to those of calamitic liquid crystals.     

Molecular alignment of rigid rods in nonrigid spherical pores

We have investigated the orientation ordering of two shish-kebab chains confined by spherically harmonic potentials through Monte Carlo simulations and asymptotic analysis. The rigid rod is modeled as shish-kebab chains consisting of tangent hard spheres aligned in the same axis, and the harmonic potential is chosen to model nonrigid cavities. We first show that the interactions between a rod and the spherically harmonic potential are independent of chain orientation, indicating that the alignment of two confined rods arises from the excluded volume interactions alone. In the strong fields, the order parameter of two confined rods converges to different values, depending on the parity of chain length. From asymptotic order parameters, we find that the rods of odd-number beads rotate more freely even under the limiting strong confinement. However, the two rods of even-number beads are essentially trapped in a configuration of perpendicular alignment through intercalation of their central grooves. We attribute the dependence of the parity of chain length to the different locations of the center-of-mass in a rod for these two cases. Furthermore, we compare the shish-kebab chains with different rod models in the simulations, and utilize these models to explore the effect of the local rod smoothness on molecular alignment. Our findings suggest that increasing local rod smoothness enhances the rotational degree of freedom for confined rods, and the effect of local rod roughness emerges under strong enough applied potentials.     

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.     

Comparative study of temperature dependent orientational relaxation in a model thermotropic liquid crystal and in a model supercooled liquid

Recent optical Kerr effect experiments have revealed a power law decay of the measured signal with a temperature independent exponent at short-to-intermediate times for a number of liquid crystals in the isotropic phase near the isotropic-nematic transition and supercooled molecular liquids above the mode coupling theory critical temperature. In this work, the authors investigate the temperature dependence of short-to-intermediate time orientational relaxation in a model thermotropic liquid crystal across the isotropic-nematic transition and in a binary mixture across the supercooled liquid regime in molecular dynamics simulations. The measure of the experimentally observable optical Kerr effect signal is found to follow a power law decay at short-to-intermediate times for both systems in agreement with recent experiments. In addition, the temperature dependence of the power law exponent is found to be rather weak. As the model liquid crystalline system settles into the nematic phase upon cooling, the decay of the single-particle second-rank orientational time correlation function exhibits a pattern that is similar to what has been observed for supercooled liquids.