Structure and internal dynamics of a side chain liquid crystalline polymer in various phases by molecular dynamics simulations: a step towards coarse graining

Impact of Molecular-level Structural Disruption on Relaxation Dynamics of Polymers with End-on and Side-on Liquid Crystal Moieties

In side-chain liquid crystal polymers (SCLCPs), short side chains are attached on a flexible polymer backbone, and each side chain can have a liquid crystal (LC) group attached at the final bead in either an end-on or a side-on configuration. SCLCPs with random sequences of end-on and side-on LC moieties exhibit nonmonotonic thermal behavior as a function of composition, with some mixed sequences having a lower isotropic to LC phase transition than either purely end-on or side-on configurations. The origin of this nonmonotonic thermal trend lies in the disruption of molecular-level positional ordering and alignment due to the different preferred types of ordering of the different LC attachment types. We compare coarse-grained molecular dynamics (MD) simulations and experiments on SCLCP systems with only one type of LC moiety and demonstrate qualitative agreement in the observed mesophases of end-on and side-on SCLCP systems. Specifically, end-on SCLCPs display a smectic B-like mesophase, with layers of polymer between LC layers, while side-on SCLCPs exhibit a quasi-hexagonal columnar structure of polymer and a nematic surrounding the LC mesophase. Detailed analysis of SCLCP systems with various compositions of these types of LC attachments via MD reveals structural disruption in systems with intermediate compositions. Simulation snapshots and anisotropy ratio measurements show how random SCLCP systems deviate from the expected behavior of prolate or oblate systems in terms of their conformation. This molecular disruption in random SCLCP systems, particularly with a high composition of side-on LC moieties, also significantly impacts the relaxation dynamics. Modifying the composition of the LC type of attachment (molecular structure) is a possible route to tuning both the phase behavior and mechanical response of these systems.

Coarse-grained modeling of polymers with end-on and side-on liquid crystal moieties: Effect of architecture

Mesogens, which are typically stiff rodlike or disklike molecules, are able to self-organize into liquid crystal (LC) phases in a certain temperature range. Such mesogens, or LC groups, can be attached to polymer chains in various configurations including within the backbone (main-chain LC polymers) or at the ends of side-chains attached to the backbone in an end-on or side-on configuration (side-chain LC polymers or SCLCPs), which can display synergistic properties arising from both their LC and polymeric character. At lower temperatures, chain conformations may be significantly altered due to the mesoscale LC ordering; thus, when heated from the LC ordered state through the LC to isotropic phase transition, the chains return from a more stretched to a more random coil conformation. This can cause macroscopic shape changes, which depend significantly on the type of LC attachment and other architectural properties of the polymer. Here, to study the structure-property relationships for SCLCPs with a range of different architectures, we develop a coarse-grained model that includes torsional potentials along with LC interactions of a Gay-Berne form. We create systems of different side-chain lengths, chain stiffnesses, and LC attachment types and track their structural properties as a function of temperature. Our modeled systems indeed form a variety of well-organized mesophase structures at low temperatures, and we predict higher LC-to-isotropic transition temperatures for the end-on side-chain systems than for analogous side-on side-chain systems. Understanding these phase transitions and their dependence on polymer architecture can be useful in designing materials with reversible and controllable deformations.

Nematic liquid crystalline elastomers are aeolotropic materials

Continuum models describing ideal nematic solids are widely used in theoretical studies of liquid crystal elastomers. However, experiments on nematic elastomers show a type of anisotropic response that is not predicted by the ideal models. Therefore, their description requires an additional term coupling elastic and nematic responses, to account for aeolotropic effects. In order to better understand the observed elastic response of liquid crystal elastomers, we analyse theoretically and computationally different stretch and shear deformations. We then compare the elastic moduli in the infinitesimal elastic strain limit obtained from the molecular dynamics simulations with the ones derived theoretically, and show that they are better explained by including nematic order effects within the continuum framework.

Deconstructing the behavior of donor–acceptor copolymers in solution & the melt: the case of PTB7

Molecular dynamics simulations of the donor–acceptor copolymer PTB7 at near experimental scale reveal structure–dynamics correlations in the condensed phase.

Molecular organizations of conical mesogenic fullerenes

We have studied liquid crystal phases formed by fullerenes functionalized with mesogenic groups yielding a cone-shaped molecular structure. We have modelled these shuttlecock-like molecules with a set of Gay-Berne particles grafted with flexible springs to a spherical core and we have studied, using Monte Carlo simulations, their phase organization, also with a view to examining their possible use as candidate organic photovoltaic materials. We have found that, upon cooling from the isotropic phase, the system forms a columnar phase, like in the experimental work of Kato and coworkers [T. Kato et al., Nature, 2002, 419, 702]. However the phase is made of polar stacks extending not more than about ten molecules, which could limit their usefulness in enhancing and directing charge transport for possible photovoltaic applications.

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.

Diffusivity maximum in a reentrant nematic phase

We report molecular dynamics simulations of confined liquid crystals using the Gay–Berne–Kihara model. Upon isobaric cooling, the standard sequence of isotropic–nematic–smectic A phase transitions is found. Upon further cooling a reentrant nematic phase occurs. We investigate the temperature dependence of the self-diffusion coefficient of the fluid in the nematic, smectic and reentrant nematic phases. We find a maximum in diffusivity upon isobaric cooling. Diffusion increases dramatically in the reentrant phase due to the high orientational molecular order. As the temperature is lowered, the diffusion coefficient follows an Arrhenius behavior. The activation energy of the reentrant phase is found in reasonable agreement with the reported experimental data. We discuss how repulsive interactions may be the underlying mechanism that could explain the occurrence of reentrant nematic behavior for polar and non-polar molecules.

A new anisotropic soft-core model for the simulation of liquid crystal mesophases

A new anisotropic soft-core model is presented, which is suitable for the rapid simulation of liquid crystal mesophases. The potential is based on a soft spherocylinder, which can be easily tuned to favor different liquid crystal mesophases. The soft-core nature of the potential makes it suitable for long-time step molecular dynamics or dissipative particle dynamics simulations, particularly as a reference model for mesogens or as an anisotropic solvent for use in combination with atomistic models. Results are presented for two variants of the new potential, which show different mesophase behaviors. Variants of the potential can also be linked together to produce more complicated molecular structures. Here, as an example, results are provided for a model multipedal liquid crystal, which has eight liquid crystalline groups linked to a central core via semiflexible chains. Here, despite the complexity of molecular structure, the model succeeds in showing the spontaneous formation of a liquid crystal phase. The results also demonstrate that there is a very strong coupling between the internal structure of the multipedal mesogen and the molecular order of the phase, with the mesogen spontaneously undergoing major structural rearrangement at the transition to the liquid crystal phase.