Calculation of the rotational viscosity of a nematic liquid crystal

Structures, thermodynamics and dynamics of topological defects in Gay-Berne nematic liquid crystals

Topological defects are a ubiquitous phenomenon across different physical systems. A better understanding of defects can be helpful in elucidating the physical behaviors of many real materials systems. In nematic liquid crystals, defects exhibit unique optical signatures and can segregate impurities, showing their promise as molecular carriers and nano-reactors. Continuum theory and simulations have been successfully applied to link static and dynamical behaviors of topological defects to the material constants of the underlying nematic. However, further evidence and molecular details are still lacking. Here we perform molecular dynamics simulations of Gay-Berne particles, a model nematic, to examine the molecular structures and dynamics of +1/2 defects in a thin-film nematic. Specifically, we measure the bend-to-splay ratio K₃/K₁ using two independent, indirect measurements, showing good agreement. Next, we study the annihilation event of a pair of ±1/2 defects, of which the trajectories are consistent with experiments and hydrodynamic simulations. We further examine the thermodynamics of defect annihilation in an NVE ensemble, leading us to correctly estimate the elastic modulus by using the energy conservation law. Finally, we explore effects of defect annihilation in regions of nonuniform temperature within these coarse-grained molecular models which cannot be analysed by existing continuum level simulations. We find that +1/2 defects tend to move toward hotter areas and their change of speed in a temperature gradient can be quantitatively understood through a term derived from the temperature dependence of the elastic modulus. As such, our work has provided molecular insights into structures and dynamics of topological defects, presented unique and accessible methods to measure elastic constants by inspecting defects, and proposed an alternative control parameter of defects using temperature gradient.

Molecular Simulation Approaches to the Study of Thermotropic and Lyotropic Liquid Crystals

Over the last decade, the availability of computer time, together with new algorithms capable of exploiting parallel computer architectures, has opened up many possibilities in molecularly modelling liquid crystalline systems. This perspective article points to recent progress in modelling both thermotropic and lyotropic systems. For thermotropic nematics, the advent of improved molecular force fields can provide predictions for nematic clearing temperatures within a 10 K range. Such studies also provide valuable insights into the structure of more complex phases, where molecular organisation may be challenging to probe experimentally. Developments in coarse-grained models for thermotropics are discussed in the context of understanding the complex interplay of molecular packing, microphase separation and local interactions, and in developing methods for the calculation of material properties for thermotropics. We discuss progress towards the calculation of elastic constants, rotational viscosity coefficients, flexoelectric coefficients and helical twisting powers. The article also covers developments in modelling micelles, conventional lyotropic phases, lyotropic phase diagrams, and chromonic liquid crystals. For the latter, atomistic simulations have been particularly productive in clarifying the nature of the self-assembled aggregates in dilute solution. The development of effective coarse-grained models for chromonics is discussed in detail, including models that have demonstrated the formation of the chromonic N and M phases.

Molecular insight into the microstructure and microscopic dynamics of pyridinium ionic liquids with different alkyl chains based on temperature response

High temperature is advantageous to the aggregation of the polar regions as well as the nonpolar regions of pyridinium ionic liquids.

Optimization of the GAFF force field to describe liquid crystal molecules: the path to a dramatic improvement in transition temperature predictions

Systematic optimization of the General Amber Force Field (GAFF) for mesogenic fragments leads to a dramatic improvement in the modelling of liquid crystal clearing points.

Liquid Crystal Properties of the n-Alkyl-cyanobiphenyl Series from Atomistic Simulations with Ab Initio Derived Force Fields

Lengthy molecular dynamics (MD) simulations were performed at constant atmospheric pressure and different temperatures for the series of the 4-n-alkyl-4'-cyanobiphenyls (nCB) with n = 6, 7, and 8. The accurate atomistic force field (Bizzarri, M.; Cacelli, I.; Prampolini, G; Tani, A. J. Phys. Chem. A 2004, 108, 10336), successfully employed to reproduce thermodynamic and transport properties of the 5CB molecule, has here been extended to higher homologues. Nematic and isotropic phases were found for all members of the series, and also, a smectic phase was (tentatively) identified for 8CB at 1 atm and 300 K. Transition temperatures reproduce the experimental values within +/-10 K. Also, structural properties as second and fourth rank orientational order parameters are in good agreement with the corresponding experimental quantities. This means that the well-known odd-even effect, observed for many properties along the nCB series, is well reproduced, despite the narrow range of oscillations, e.g., in clearing temperatures. A detailed analysis of the correlation between molecular properties and odd-even effects is presented.

Molecular orientational and dipolar correlation in the liquid crystal mixture E7: a molecular dynamics simulation study at a fully atomistic level

Molecular dynamics simulations are reported for the four component nematic liquid crystal mixture E7, which is used commercially. We are able to show the growth of a nematic phase directly from an isotropic liquid over a 100 ns period for an all-atom model, and study orientational and dipole order within the nematic phase. The simulations show that the cyanoterphenyl component of the mixture, 5CT, is more ordered than the three cyanobiphenyl components. The simulations show also that both parallel and anti-parallel dipole correlation take place in E7 but that the strong anti-parallel dipole correlation is localised to particular arrangements of molecules. It is possible to identify two key preferred configurations for molecular pairs in the fluid, which explain the form of the dipole correlation function, g(1)(r).

Atomistic simulations of a thermotropic biaxial liquid crystal

We have performed molecular dynamics simulations of a 2,5-bis-(p-hydroxyphenyl)-1,3,4-oxadiazole mesogen (ODBP-Ph-C(7)) at a fully atomistic level for a range of temperatures within the region that has experimentally been assigned to a biaxial nematic phase. Analysis of the data shows that the simulated nematic phase is biaxial but that the degree of biaxiality is small. The simulations show also the formation of ferroelectric domains in the nematic where the molecular short axis is aligned with the oxadiazole dipoles parallel to each other. Removal of electrostatic interactions leads to destabilization of ferroelectric domains and destabilization of the biaxiality. An additional simulation shows the slow growth of a mesophase directly from the isotropic fluid over a period of approximately 50 ns. This is the first time this has been achieved within the framework of an all-atom model.

Molecular dynamic study of the odd-even effect in some 4-n-alkyl-4'-cyanobiphenyls

Using a united atom model, we performed molecular dynamics (MD) simulations on the 4-n-alkyl-4'-cyanobiphenyl (nCB) liquid crystal homologous series for at a constant reduced temperature in the nematic phase. To evaluate the rotational diffusion coefficient (RDC), the second- and fourth-rank order parameters (OP's) and correlation time were calculated from MD trajectories. We analyzed the rotational viscosity coefficient (RVC) by using the Nemtsov-Zakharov and Fialkowski methods based on statistical-mechanical approaches. The simulated RDC, RVC, and OP's were found to be reasonable when compared with previous theoretical and experimental data. These quantities were also calculated for the rigid core and alkyl chain segments separately for each mesogen, to take a closer look at the molecule tail-chain flexibility. The properties calculated for a rigid core were compared with those for an alkyl chain, and satisfactory results were obtained. Odd-even effects for RDF, RVC, and molecular ordering parameters were calculated for five mesogens.

Molecular dynamics simulations of side chain liquid crystal polymer molecules in isotropic and liquid-crystalline melts

A detailed molecular dynamics simulation study is described for a polysiloxane side chain liquid crystal polymer (SCLCP). The simulations use a coarse-grained model composed of a combination of isotropic and anisotropic interaction sites. On cooling from a fully isotropic polymer melt, we see spontaneous microphase separation into polymer-rich and mesogen-rich regions. Upon application of a small aligning potential during cooling, the structures that form on microphase separation anneal to produce a smectic-A phase in which the polymer backbone is largely confined between the smectic layers. Several independent quenches from the melt are described that vary in the strength of the aligning potential and the degree of cooling. In each quench, defects were found where the backbone chains hop from one backbone-rich region to the next by tunneling through the mesogenic layers. As expected, the number of such defects is found to depend strongly on the rate of cooling. In the vicinity of such a defect, the smectic-A structure of the mesogen-rich layers is disrupted to give nematiclike ordering. Additionally, several extensive annealing runs of approximately 40 ns duration have been carried out at the point of microphase separation. During annealing the polymer backbone is seen to be slowly excluded from the mesogenic layers and lie perpendicular to the smectic-A director. These observations agree with previous assumptions about the structure of a SCLCP and with interpretations of x-ray diffraction and small angle neutron scattering data. The flexible alkyl spacers, which link the backbone to the mesogens, are found to form sublayers around the backbone layer.

Calculation of flexoelectric coefficients for a nematic liquid crystal by atomistic simulation

Equilibrium molecular dynamics calculations have been performed for the liquid crystal molecule n-4-(trans-4-n-pentylcyclohexyl)benzonitrile (PCH5) using a fully atomistic model. Simulation data have been obtained for a series of temperatures in the nematic phase. The simulation data have been used to calculate the flexoelectric coefficients e(s) and e(b) using the linear response formalism of Osipov and Nemtsov [M. A. Osipov and V. B. Nemtsov, Sov. Phys. Crstallogr. 31, 125 (1986)]. The temperature and order parameter dependence of e(s) and e(b) are examined, as are separate contributions from different intermolecular interactions. Values of e(s) and e(b) calculated from simulation are consistent with those found from experiment.