The role of attractive interactions in rod–sphere mixtures

Gaussian representation of coarse-grained interactions of liquids: Theory, parametrization, and transferability

Coarse-grained (CG) interactions determined via bottom-up methodologies can faithfully reproduce the structural correlations observed in fine-grained (atomistic resolution) systems, yet they can suffer from limited extensibility due to complex many-body correlations. As part of an ongoing effort to understand and improve the applicability of bottom-up CG models, we propose an alternative approach to address both accuracy and transferability. Our main idea draws from classical perturbation theory to partition the hard sphere repulsive term from effective CG interactions. We then introduce Gaussian basis functions corresponding to the system's characteristic length by linking these Gaussian sub-interactions to the local particle densities at each coordination shell. The remaining perturbative long-range interaction can be treated as a collective solvation interaction, which we show exhibits a Gaussian form derived from integral equation theories. By applying this numerical parametrization protocol to CG liquid systems, our microscopic theory elucidates the emergence of Gaussian interactions in common phenomenological CG models. To facilitate transferability for these reduced descriptions, we further infer equations of state to determine the sub-interaction parameter as a function of the system variables. The reduced models exhibit excellent transferability across the thermodynamic state points. Furthermore, we propose a new strategy to design the cross-interactions between distinct CG sites in liquid mixtures. This involves combining each Gaussian in the proper radial domain, yielding accurate CG potentials of mean force and structural correlations for multi-component systems. Overall, our findings establish a solid foundation for constructing transferable bottom-up CG models of liquids with enhanced extensibility.

Hidden scale invariance in the Gay-Berne model

This paper presents a numerical study of the Gay-Berne liquid crystal model with parameters corresponding to calamitic (rod-shaped) molecules. The focus is on the isotropic and nematic phases at temperatures above unity, where we find strong correlations between the virial and potential-energy thermal fluctuations, reflecting the hidden scale invariance symmetry. This implies the existence of isomorphs, which are curves in the thermodynamic phase diagram of approximately invariant physics. We study numerically one isomorph in the isotropic phase and one in the nematic phase. In both cases, good invariance of the dynamics is demonstrated via data for the mean-square displacement and the reduced-unit time-autocorrelation functions of the velocity, angular velocity, force, torque, and first- and second-order Legendre polynomial orientational order parameters. Deviations from isomorph invariance are observed at short times for the orientational time-autocorrelation functions, which reflects the fact that the moment of inertia is assumed to be constant and thus not isomorph-invariant in reduced units. Structural isomorph invariance is demonstrated from data for the radial distribution functions of the molecules and their orientations. For comparison, all quantities were also simulated along an isochore of similar temperature variation, in which case invariance is not observed. We conclude that the thermodynamic phase diagram of the calamitic Gay-Berne model is essentially one-dimensional in the studied regions as predicted by isomorph theory, a fact that potentially allows for simplifications of future theories and numerical studies.

Field-induced anti-nematic and biaxial ordering in binary mixtures of discotic mesogens and spherical magnetic nanoparticles

Using computer simulations we explore the equilibrium structure and response to external stimuli of complex magnetic hybrids consisting of magnetic particles in discotic liquid crystalline matrices. We show that the anisotropy of the liquid crystalline matrix (either in the nematic or in the columnar phase) promotes the collective orientational ordering of self-assembled magnetic particles. Upon applying an external homogeneous magnetic field in an otherwise isotropic state, the magnetic particles self-assemble into linear-rodlike-chains. At the same time structural changes occur in the matrix. The matrix transforms from an isotropic to a non-conventional anti-nematic state in which the symmetry axis of the discs is, on average, perpendicular to the magnetic field. In addition, a stable biaxial nematic state is found upon applying an external field to an otherwise uniaxial discotic nematic state. These observed morphologies constitute an appealing alternative to binary mixtures of rigid rod-disc system and indicate that non-trivial biaxial ordering can be obtained in the presence of a uniaxial external stimulus.

Dynamics and Pretransitional Effects in C60 Fullerene Nanoparticles and Liquid Crystalline Dodecylcyanobiphenyl (12CB) Hybrid System

The report shows the strong impact of fullerene C60 nanoparticles on phase transitions and complex dynamics of rod-like liquid crystal dodecylcyanobiphenyl (12CB), within the limit of small concentrations. Studies were carried out using broadband dielectric spectroscopy (BDS) via the analysis of temperature dependences of the dielectric constant, the maximum of the primary loss curve, and relaxation times. They revealed a strong impact of nanoparticles, leading to a ~20% change of dielectric constant even at x = 0.05% of C60 fullerene. The application of the derivative-based and distortion-sensitive analysis showed that pretransitional effects dominate in the isotropic liquid phase up to 65 K above the clearing temperature and in the whole Smectic A mesophase. The impact of nanoparticles on the pretransitional anomaly appearance is notable for the smectic-solid phase transition. The fragility-based analysis of relaxation times revealed the universal pattern of its temperature changes, associated with scaling via the "mixed" ("activated" and "critical") relation. Phase behavior and dynamics of tested systems are discussed within the extended Landau-de Gennes-Ginzburg mesoscopic approach.

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.

The interplay of interfaces, supramolecular assembly, and electronics in organic semiconductors

Organic semiconductors, which include a diverse range of carbon-based small molecules and polymers with interesting optoelectronic properties, offer many advantages over conventional inorganic semiconductors such as silicon and are growing in importance in electronic applications. Although these materials are now the basis of a lucrative industry in electronic displays, many promising applications such as photovoltaics remain largely untapped. One major impediment to more rapid development and widespread adoption of organic semiconductor technologies is that device performance is not easily predicted from the chemical structure of the constituent molecules. Fundamentally, this is because organic semiconductor molecules, unlike inorganic materials, interact by weak non-covalent forces, resulting in significant structural disorder that can strongly impact electronic properties. Nevertheless, directional forces between generally anisotropic organic-semiconductor molecules, combined with translational symmetry breaking at interfaces, can be exploited to control supramolecular order and consequent electronic properties in these materials. This review surveys recent advances in understanding of supramolecular assembly at organic-semiconductor interfaces and its impact on device properties in a number of applications, including transistors, light-emitting diodes, and photovoltaics. Recent progress and challenges in computer simulations of supramolecular assembly and orientational anisotropy at these interfaces is also addressed.

A coarse-grained model of ionic liquid crystals: the effect of stoichiometry on the stability of the ionic nematic phase

We have investigated, by means of molecular dynamics simulations, the phase behaviour of mixtures of charged ellipsoidal Gay-Berne (GB) particles and spherical Lennard-Jones (LJ) particles, as a coarse-grained model of ionic liquid crystals (ILCs). The anisotropic GB particles represent cations usually found in ILCs, for example, pyridinium or bipyridinium salts, while the spherical LJ particles are taken as a model of anions like common halides, hexafluorophosphate and tetrafluoroborate. Here we have focused our attention on the effect of the stoichiometry of the system (that is, the GB : LJ ratio n : m in the salt formula [GB]n[LJ]m) on the stability and thermal range of the ionic liquid crystal phases formed, with special attention to the ionic nematic phase. To isolate the stoichiometry effect, a comparison of four different systems with GB : LJ ratios of 1 : 3, 1 : 2, 1 : 1 and 2 : 1 is made by keeping the packing fraction and the charge of the minor component fixed. Our results suggest a way to improve the stability of the ionic nematic phase by enhancing the anisotropic van der Waals interaction compared to the Coulomb interaction, and by increasing the proportion of anisotropic particles in the mixture.

Understanding Missing Entropy in Coarse-Grained Systems: Addressing Issues of Representability and Transferability

Coarse-grained (CG) models facilitate efficient simulation of complex systems by integrating out the atomic, or fine-grained (FG), degrees of freedom. Systematically derived CG models from FG simulations often attempt to approximate the CG potential of mean force (PMF), an inherently multidimensional and many-body quantity, using additive pairwise contributions. However, they currently lack fundamental principles that enable their extensible use across different thermodynamic state points, i.e., transferability. In this work, we investigate the explicit energy-entropy decomposition of the CG PMF as a means to construct transferable CG models. In particular, despite its high-dimensional nature, we find for liquid systems that the entropic component to the CG PMF can similarly be represented using additive pairwise contributions, which we show is highly coupled to the CG configurational entropy. This approach formally connects the missing entropy that is lost due to the CG representation, i.e., translational, rotational, and vibrational modes associated with the missing degrees of freedom, to the CG entropy. By design, the present framework imparts transferable CG interactions across different temperatures due to the explicit definition of an additive entropic contribution. Furthermore, we demonstrate that transferability across composition state points, such as between bulk liquids and their mixtures, is also achieved by designing combining rules to approximate cross-interactions from bulk CG PMFs. Using the predicted CG model for liquid mixtures, structural correlations of the fitted CG model were found to corroborate a high-fidelity combining rule. Our findings elucidate the physical nature and compact representation of CG entropy and suggest a new approach for overcoming the transferability problem. We expect that this approach will further extend the current view of CG modeling into predictive multiscale modeling.

Tuning Coulombic interactions to stabilize nematic and smectic ionic liquid crystal phases in mixtures of charged soft ellipsoids and spheres

We have investigated the effect of electrostatic interactions in mixtures of soft ellipsoids and spheres based on the well-known Gay-Berne (GB) and Lennard-Jones (LJ) potential, respectively. These model systems, in their original version, that is without any electrostatic charge, have been thoroughly investigated in the literature both as pure components and mixtures. In particular, mixtures of particles of different shapes, such as spheres and ellipsoids, tend to phase separate because of the excluded volume effects. Common ionic liquid crystals, based on imidazolium or other quaternary ammonium salts, are usually composed of roughly elongated (although flexible) cations and roughly spherical anions, that is, particles with a similar shape such as the GB and LJ models. Therefore, in this work, we present the results of molecular dynamics simulations of mixtures of positively charged GB and negatively charged LJ particles as models of ionic liquid crystals. Interestingly, by modulating the charge of the particles it is possible to stabilize isotropic, nematic, smectic and crystalline ionic phases. The relative weight of Coulomb (a radial, therefore isotropic interaction) and van der Waals (an anisotropic interaction) contributions is a key parameter to tune the stability of various mesophases.

Phase diagram of binary colloidal rod-sphere mixtures from a 3D real-space analysis of sedimentation-diffusion equilibria

Self-assembly of binary particle systems offers many new opportunities for materials science. Here, we studied sedimentation equilibria of silica rods and spheres, using quantitative 3D confocal microscopy. We determined not only pressure, density and order parameter profiles, but also the experimental phase diagram exhibiting a stable binary smectic liquid-crystalline phase (Sm₂). Using computer simulations we confirmed that the Sm₂-phase can be stabilized by entropy alone, which opens up the possibility of combining new materials properties at a wide array of length scales.

Orientational order and translational dynamics of magnetic particle assemblies in liquid crystals

Implementing extensive molecular dynamics simulations we explore the organization of magnetic particle assemblies (clusters) in a uniaxial liquid crystalline matrix comprised of rodlike particles. The magnetic particles are modelled as soft dipolar spheres with diameter significantly smaller than the width of the rods. Depending on the dipolar strength coupling the magnetic particles arrange into head-to-tail configurations forming various types of clusters including rings (closed loops) and chains. In turn, the liquid crystalline matrix induces long range orientational ordering to these structures and promotes their diffusion along the director of the phase. Different translational dynamics are exhibited as the liquid crystalline matrix transforms either from isotropic to nematic or from nematic to smectic state. This is caused due to different collective motion of the magnetic particles into various clusters in the anisotropic environments. Our results offer a physical insight for understanding both the structure and dynamics of magnetic particle assemblies in liquid crystalline matrices.

Probing a self-assembled fd virus membrane with a microtubule

The self-assembly of highly anisotropic colloidal particles leads to a rich variety of morphologies whose properties are just beginning to be understood. This article uses computer simulations to probe a particle-scale perturbation of a commonly studied colloidal assembly, a monolayer membrane composed of rodlike fd viruses in the presence of a polymer depletant. Motivated by experiments currently in progress, we simulate the interaction between a microtubule and a monolayer membrane as the microtubule "pokes" and penetrates the membrane face-on. Both the viruses and the microtubule are modeled as hard spherocylinders of the same diameter, while the depletant is modeled using ghost spheres. We find that the force exerted on the microtubule by the membrane is zero either when the microtubule is completely outside the membrane or when it has fully penetrated the membrane. The microtubule is initially repelled by the membrane as it begins to penetrate but experiences an attractive force as it penetrates further. We assess the roles played by translational and rotational fluctuations of the viruses and the osmotic pressure of the polymer depletant. We find that rotational fluctuations play a more important role than the translational ones. The dependence on the osmotic pressure of the depletant of the width and height of the repulsive barrier and the depth of the attractive potential well is consistent with the assumed depletion-induced attractive interaction between the microtubule and viruses. We discuss the relevance of these studies to the experimental investigations.

Colloidal Suspensions of Rodlike Nanocrystals and Magnetic Spheres under an External Magnetic Stimulus: Experiment and Molecular Dynamics Simulation

Using experiments and molecular dynamics simulations, we explore magnetic field-induced phase transformations in suspensions of nonmagnetic rodlike and magnetic sphere-shaped particles. We experimentally demonstrate that an external uniform magnetic field causes the formation of small, stable clusters of magnetic particles that, in turn, induce and control the orientational order of the nonmagnetic subphase. Optical birefringence was studied as a function of the magnetic field and the volume fractions of each particle type. Steric transfer of the orientational order was investigated by molecular dynamics (MD) simulations; the results are in qualitative agreement with the experimental observations. By reproducing the general experimental trends, the MD simulation offers a cohesive bottom-up interpretation of the physical behavior of such systems, and it can also be regarded as a guide for further experimental research.

Tunable structures of mixtures of magnetic particles in liquid-crystalline matrices

We investigate the self-organization of a binary mixture of similar sized rods and dipolar soft spheres by means of Monte-Carlo simulations. We model interparticle interactions by employing anisotropic Gay-Berne, dipolar and soft-sphere interactions. In the limit of vanishing magnetic moments we obtain a variety of fully miscible liquid crystalline phases including nematic, smectic and lamellar phases. For the magnetic mixture, we find that the liquid crystalline matrix supports the formation of orientationally ordered ferromagnetic chains. Depending on the relative size of the species the chains align parallel or perpendicular to the director of the rods forming uniaxial or biaxial nematic, smectic and lamellar phases. As an exemplary external perturbation we apply a homogeneous magnetic field causing uniaxial or biaxial ordering to an otherwise isotropic state.

Spontaneous ordering of magnetic particles in liquid crystals: From chains to biaxial lamellae

Using Monte Carlo computer simulations we explore the self-assembly and ordering behavior of a hybrid, soft magnetic system consisting of small magnetic nanospheres in a liquid-crystalline (LC) matrix. Inspired by recent experiments with colloidal rod matrices, we focus on conditions where the sphere and rod diameters are comparable. Already in the absence of a magnetic field, the nematic ordering of the LC can stabilize the formation of magnetic chains along the nematic or smectic director, yielding a state with local (yet no macroscopic) magnetic order. The chains, in turn, increase the overall nematic order, reflecting the complex interplay of the structure formation of the two components. When increasing the sphere diameter, the spontaneous uniaxial ordering is replaced by biaxial lamellar morphologies characterized by alternating layers of rods and magnetic chains oriented perpendicular to the rod's director. These ordering scenarios at zero field suggest a complex response of the resulting hybrid to external stimuli, such as magnetic fields and shear forces.

Shape alloys of nanorods and nanospheres from self-assembly

Mixtures of anisotropic nanocrystals promise a great diversity of superlattices and phase behaviors beyond those of single-component systems. However, obtaining a colloidal shape alloy in which two different shapes are thermodynamically coassembled into a crystalline superlattice has remained a challenge. Here we present a joint experimental-computational investigation of two geometrically ubiquitous nanocrystalline building blocks-nanorods and nanospheres-that overcome their natural entropic tendency toward macroscopic phase separation and coassemble into three intriguing phases over centimeter scales, including an AB2-type binary superlattice. Monte Carlo simulations reveal that, although this shape alloy is entropically stable at high packing fraction, demixing is favored at experimental densities. Simulations with short-ranged attractive interactions demonstrate that the alloy is stabilized by interactions induced by ligand stabilizers and/or depletion effects. An asymmetry in the relative interaction strength between rods and spheres improves the robustness of the self-assembly process.

Formation of double glass in binary mixtures of anisotropic particles: dynamic heterogeneities in rotations and displacements

We study glass behavior in a mixture of elliptic and circular particles in two dimensions at low temperatures using an orientation-dependent Lennard-Jones potential. The ellipses have a mild aspect ratio (∼1.2) and tend to align at low temperatures, while the circular particles play the role of impurities disturbing the ellipse orientations at a concentration of 20%. These impurities have a size smaller than that of the ellipses and attract them in the homeotropic alignment. As a result, the coordination number around each impurity is mostly 5 or 4 in glassy states. We realize double glass, where both the orientations and the positions are disordered but still hold mesoscopic order. We find a strong heterogeneity in the flip motions of the ellipses, which sensitively depends on the impurity clustering. In our model, a small fraction of the ellipses still undergo flip motions relatively rapidly even at low temperatures. In contrast, the nonflip rotations (with angle changes not close to ±π) are mainly caused by the cooperative configuration changes involving many particles. Then, there arises a long-time heterogeneity in the nonflip rotations closely correlated with the dynamic heterogeneity in displacements.

Self-assembly of 2D membranes from mixtures of hard rods and depleting polymers()

We combine simulations and experiments to elucidate the molecular forces leading to the assembly of two dimensional membrane-like structures composed of a one rod-length thick monolayer of aligned rods from an immiscible suspension of hard rods and depleting polymers. Computer simulations predict that monolayer membranes are thermodynamically stable above a critical rod aspect ratio and below a critical depletion interaction length scale. Outside of these conditions alternative structures such as stacked smectic columns or nematic droplets are thermodynamically stable. These predictions are confirmed by subsequent experiments using a model system of virus rod-like molecules and non-adsorbing polymer. Our work demonstrates that collective molecular protrusion fluctuations alone are sufficient to stabilize membranes composed of homogenous rods with simple excluded volume interactions.

Theoretical calculation of the phase behavior of colloidal membranes

We formulate a density functional theory that describes the phase behavior of hard rods and depleting polymers, as realized in recent experiments on suspensions of fd virus and nonadsorbing polymer. The theory predicts the relative stability of nematic droplets, stacked smectic columns, and a recently discovered phase of isolated monolayers of rods, or colloidal membranes. We find that a minimum rod aspect ratio is required for stability of colloidal membranes and that collective protrusion undulations are the dominant effect that stabilizes this phase. The theoretical predictions are shown to be qualitatively consistent with experimental and computational results.

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 dynamics simulations of nanoparticles in dense isotropic nematogens: the role of matrix-induced long-range repulsive interactions

We have carried out molecular dynamics simulation studies of binary mixtures of spherical nanoparticles (NPs) in a matrix of dense isotropic rod-shaped nematogens, with the size of the nematogen length being similar to that of the NP diameter. NPs at even low concentrations were found to shift the isotropic-nematic (I-N) transition significantly to higher pressure at a given temperature, indicative of long-range perturbation of the nematogenic matrix by the NPs. The NPs were found to be dispersed in the dense isotropic nematogenic matrix over a wide range of NP concentrations due to long-range (compared with the molecular size of the nematogens) repulsion caused by NP-induced local order fluctuations and reduced local orientational correlation in the isotropic nematogenic matrix, in contrast to the phase separation predicted and observed in other studies where the particles were much larger or smaller than the nematogens. Furthermore, since the repulsion observed in the NP-nematogen mixtures is only microscopically long range (on the order of about ten molecular lengths of the nematogens), globally ordered clustering observed in mixtures of colloidal particles in nematic matrices resulting from macroscopically long-range interaction is not observed in our simulations.

Use of Parsons-Lee and Onsager theories to predict nematic and demixing behavior in binary mixtures of hard rods and hard spheres

Parsons-Lee and Onsager theories are formulated for the isotropic-nematic transition in a binary mixture of hard rods and hard spheres. Results for the phase coexistence and for the equation of state in both phases for mixtures with different relative sizes and composition are presented. The two theories explain correctly the general behavior observed in experiments and computer simulations for these fluids. In particular, the theory accounts for the destabilization of the nematic phase when spherical or globular macromolecules are added to a system of rodlike colloids, and the entrance of the system into a demixed regime at high volume fractions of the spherical particles. Upon demixing a nematic state rich in rods coexists in equilibrium with an isotropic state much more diluted in the rodlike component. Onsager theory fails on quantitative grounds for aspect ratios of the rodlike molecules smaller than 100, and in the cases where the molar fractions of spheres becomes close to unity. On the contrary, the Parsons-Lee approximation remains accurate down to aspect ratios as small as 5. The spinodal analysis indicates that the isotropic-isotropic and nematic-nematic coexistences become feasible for sufficiently large spheres and long rods, respectively. The latter type of coexistence interferes partially with the isotropic-nematic coexistence regime of interest to the present work. Overall, the study serves to rationalize and control key aspects of the behavior of these binary nematogenic colloidal systems, which can be tuned with an appropriate choice of the relative size and molar fractions of the particles.

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.

A computer simulation study of the formation of liquid crystal nanodroplets from a homogeneous solution

The aggregation of liquid crystal nanodroplets from a homogeneous solution is an important but not well understood step in the preparation of various advanced photonic materials. Here, the authors performed molecular dynamics computer simulations of the formation of liquid crystalline nanodroplets, starting from an isotropic and uniform binary solution of spherical Lennard-Jones (solvent) and elongated ellipsoidal Gay-Berne (solute) rigid particles in low (<10%) concentration. They studied the dynamics of demixing and the mesogen ordering process and characterized the resulting nanodroplets assessing the effect of temperature, composition, and specific solute-solvent interaction on the morphology, structure, and anisotropy. They find that the specific solute-solvent interaction, composition, and temperature can be adjusted to tune the nanodroplet growth and size.

Coarse-grained simulation of amphiphilic self-assembly

The authors present a computer simulation study of amphiphilic self-assembly performed using a computationally efficient single-site model based on Gay-Berne [J. Chem. Phys. 74, 3316 (1981)] and Lennard-Jones particles. Molecular dynamics simulations of these systems show that free self-assembly of micellar, bilayer, and inverse micelle arrangements can be readily achieved for a single model parametrization. This self-assembly is predominantly driven by the anisotropy of the amphiphile-solvent interaction, amphiphile-amphiphile dispersive interactions being found to be of secondary importance. While amphiphile concentration is the main determinant of phase stability, molecular parameters such as head group size and interaction strength also have measurable affects on system properties.

Decoupling phenomena in supercooled liquids: signatures in the energy landscape

A significant deviation from the Debye model of rotational diffusion in the dynamics of orientational degrees of freedom in an equimolar mixture of ellipsoids of revolution and spheres is found to begin at a temperature at which the average inherent structure energy of the system starts falling with drop in temperature. We argue that this onset temperature corresponds to the emergence of the process as a distinct mode of orientational relaxation. Further, we find that the coupling between rotational and translational diffusion breaks down at a still lower temperature where a change occurs in the temperature dependence of the average inherent structure energy.

Anomalous glassy relaxation near the isotropic-nematic phase transition

Dynamical heterogeneity in a system of Gay-Berne ellipsoids near its isotropic-nematic (I-N) transition, and also in an equimolar mixture of Lennard-Jones spheres and Gay-Berne ellipsoids in deeply supercooled regime, is probed by the time evolution of non-Gaussian parameters (NGP). The appearance of a dominant second peak in the rotational NGP near the I-N transition signals the growth of pseudonematic domains. Surprisingly, such a second peak is instead observed in the translational NGP for the glassy binary mixture. Localization of orientational motion near the I-N transition is found to be responsible for the observed anomalous orientational relaxation.