Surface-Induced Ordering on Model Liquid Crystalline Dendrimers

Homogeneous alignment of liquid crystalline dendrimers confined in a slit-pore. A simulation study

In this work we present results from isobaric-isothermal (NPT) Monte Carlo simulation studies of model liquid crystalline dendrimer (LCDr) systems confined in a slit-pore made of two parallel flat walls. The dendrimers are modelled as a collection of spherical and ellipsoidal particles corresponding to the junction points of the dendritic core and to the mesogenic units respectively. Assuming planar uniform (unidirectional) soft anchoring of the mesogenic units on the substrates we investigate the conformational and alignment properties of the LCDr system at different thermodynamic state points. Tractable coarse grained force fields have been used from our previous work. At low pressures the interior of the pore is almost empty, since almost all LCDrs are anchored to the substrates forming two-dimensional smectic-like structures with the mesogens aligned along the aligning direction of the substrates. As the pressure grows the LCDrs occupy the whole pore. However, even at low temperatures, the smectic organization does not transmit in the interior of the pore and is preserved for distances of 2-3 mesogenic diameters from the walls. For this reason, the global orientational order decreases with increasing pressure (density). In the vicinity (2-3 mesogenic diameters) of the pore walls, mesogenic units preserve the smectic structure whose layers are separated by layers of spherical beads. In this region individual LCDrs possess a rod like shape.

GPU-Accelerated Molecular Dynamics Simulation to Study Liquid Crystal Phase Transition Using Coarse-Grained Gay-Berne Anisotropic Potential

Gay-Berne (GB) potential is regarded as an accurate model in the simulation of anisotropic particles, especially for liquid crystal (LC) mesogens. However, its computational complexity leads to an extremely time-consuming process for large systems. Here, we developed a GPU-accelerated molecular dynamics (MD) simulation with coarse-grained GB potential implemented in GALAMOST package to investigate the LC phase transitions for mesogens in small molecules, main-chain or side-chain polymers. For identical mesogens in three different molecules, on cooling from fully isotropic melts, the small molecules form a single-domain smectic-B phase, while the main-chain LC polymers prefer a single-domain nematic phase as a result of connective restraints in neighboring mesogens. The phase transition of side-chain LC polymers undergoes a two-step process: nucleation of nematic islands and formation of multi-domain nematic texture. The particular behavior originates in the fact that the rotational orientation of the mesogenes is hindered by the polymer backbones. Both the global distribution and the local orientation of mesogens are critical for the phase transition of anisotropic particles. Furthermore, compared with the MD simulation in LAMMPS, our GPU-accelerated code is about 4 times faster than the GPU version of LAMMPS and at least 200 times faster than the CPU version of LAMMPS. This study clearly shows that GPU-accelerated MD simulation with GB potential in GALAMOST can efficiently handle systems with anisotropic particles and interactions, and accurately explore phase differences originated from molecular structures.