Nematic order in nanoscopic liquid crystal droplets

Structure and Thermodynamics of Linear, Ring, and Catenane Polymers in Solutions and at Liquid-Liquid Interfaces

In recent decades, advances in the syntheses of mechanically interlocked macromolecules, such as catenanes, have led to much greater interest in the applications of these complexes, from molecular motors and actuators to nanoscale computational memory and nanoswitches. Much remains to be understood, however, regarding how catenated ring compounds behave as a result of the effects of different solvents as well as the effects of solvent/solvent interfaces. In this work, we have investigated, using molecular dynamics simulations, the effects of solvation of poly(ethylene oxide) chains of different topologies─linear, ring, and [2]catenane─in two solvents both considered favorable toward PEO (water, toluene) and at the water/toluene interface. Compared to ring and [2]catenane molecules, the linear PEO chain showed the largest increase in size at the water/toluene interface compared to bulk water or bulk toluene. Perhaps surprisingly, observations indicate that the tendency of all three topologies to extend at the water/toluene interface may have more to do with screening the interaction between the two solvents than with optimizing specific solvent-polymer contacts.

Modeling deformation and chaining of flexible shells in a nematic solvent with finite elements on an adaptive moving mesh

A micrometer-scale elastic shell immersed in a nematic liquid crystal may be deformed by the host if the cost of deformation is comparable to the cost of elastic deformation of the nematic. Moreover, such inclusions interact and form chains due to quadrupolar distortions induced in the host. A continuum theory model using finite elements is developed for this system, using mesh regularization and dynamic refinement to ensure quality of the numerical representation even for large deformations. From this model, we determine the influence of the shell elasticity, nematic elasticity, and anchoring condition on the shape of the shell and hence extract parameter values from an experimental realization. Extending the model to multibody interactions, we predict the alignment angle of the chain with respect to the host nematic as a function of aspect ratio, which is found to be in excellent agreement with experiments.

Nanodrops of Discotic Liquid Crystals: A Monte Carlo Study

We study the morphologies of nematic nanodrops in a vapor of a discotic nematogen by Monte Carlo simulations. The fluid interactions are modeled by a Gay-Berne model with molecular elongations of κ = 0.3 and 0.5 and different values of the energy anisotropy parameter κ' in the range of temperature T in which the nematic coexists with a vapor phase. We considered nanodrops of N = 4000 and 32 000 particles. For κ > κ', we observe that nanodrops are quite spherical (even for N = 4000 nanodrops), with a homogeneous director field for κ = 0.3 and a bipolar nematic configuration with tangential anchoring for κ = 0.5. By increasing the value of κ', nanodrops change from spherical to lens-shaped for κ = 0.3, and for κ = 0.5, spherical nanodrops with homeotropic anchoring and a disclination ring located on its equatorial plane are observed. Although no radial nanodrops are observed, isotropic liquid nanodrops with a paranematic shell and radial texture are observed for temperatures slightly above the vapor-isotropic-nematic triple point when the vapor-isotropic interface is completely wet by the nematic phase.

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.

The dynamics of 4-cyano-4-n-pentylbiphenyl (5CB) mesogen molecules located between graphene layers--MD study

The ultrathin film of mesogenic molecules 5CB confined between graphene walls has been investigated by molecular dynamics (MD) technique. The dynamical observables of 4-cyano-4-n-pentylbiphenyl (5CB) were calculated for several temperatures: the mean square displacement, translational and angular velocity autocorrelation function, and second-rank order parameter.

Rotational behaviour of red blood cells in suspension: a mesoscale simulation study

The nature of blood as a suspension of red blood cells makes computational haemodynamics a demanding task. Our coarse-grained blood model, which builds on a lattice Boltzmann method for soft particle suspensions, enables the study of the collective behaviour of the order of 10(6) cells in suspension. After demonstrating the viscosity measurement in Kolmogorov flow, we focus on the statistical analysis of the cell orientation and rotation in Couette flow. We quantify the average inclination with respect to the flow and the nematic order as a function of shear rate and haematocrit. We further record the distribution of rotation periods around the vorticity direction and find a pronounced peak in the vicinity of the theoretical value for free model cells, even though cell-cell interactions manifest themselves in a substantial width of the distribution.

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.

A large scale molecular dynamics simulation code using the fast multipole algorithm (FMD): performance and application

We present the performance of the fast classical molecular dynamics (MD) code, fast molecular dynamics (FMD), designed for efficient, object-oriented, and scalable large scale simulations, and summarize its application to a liquid crystalline cluster. FMD uses an implementation of the three-dimensional fast multipole method, developed in our group. The fast multipole method offers an efficient way (order O(N)) to handle long range electrostatic interactions, thus, enabling more realistic simulations of large molecular systems. Performance testing was carried out on IBM SP2, SGI Origin 2000, and CRAY T3E massively parallel systems using the MPI massage passing library. The electrostatic forces were tested on models of up to 100,000 randomly placed charges, and on protein and liquid crystalline molecular systems of over 99,000 atoms. Tests on the stability of the method are presented, along with comparisons with direct calculations, the NAMD2 code, and the physical multipole-based cell-multipole method.

Simulation study of the glass transition temperature in poly(methyl methacrylate)

The glass transition in syndiotactic poly (methyl methacrylate) has been studied through atomistic molecular dynamics simulations performed at temperatures in the range from 297 K to 684 K. The mean squared deviations of atoms, monomers, and molecules from their initial positions were analyzed by means of a technique that separates the effects of diffusive motion from the underlying vibrational motion. The diffusive motion shows a novel power-law variation with time, with an exponent that varies continuously from 0.5 below the glass transition temperature T(g) to 1 at high temperatures. The self part of the van Hove correlation functions for both hydrogen atoms and monomers shows structural arrest at the lowest temperature studied. A second peak in the atomic van Hove correlation is attributed to rotation of the CH3 group.