Exploring the interplay of liquid crystal orientation and spherical elastic shell deformation in spatial confinement

The application of liquid crystal technology typically relies on the precise control of molecular orientation at a surface or interface. This control can be achieved through a combination of morphological and chemical methods. Consequently, variations in constrained boundary flexibility can result in a diverse range of phase behaviors. In this study, we delve into the self-assembly of liquid crystals within elastic spatial confinement by using the Gay-Berne model with the aid of molecular dynamics simulations. Our findings reveal that a spherical elastic shell promotes a more regular and orderly alignment of liquid crystals compared to a hard shell. Moreover, during the cooling process, the hard-shell confined system undergoes an isotropic-smectic phase transition. In contrast, the phase behavior within the spherical elastic shell closely mirrors the isotropic-nematic-smectic phase transition observed in bulk systems. This indicates that the orientational arrangement of liquid crystals and the deformations induced by a flexible interface engage in a competitive interplay during the self-assembly process. Importantly, we found that phase behavior could be manipulated by altering the flexibility of the confined boundaries. This insight offers a fresh perspective for the design of innovative materials, particularly in the realm of liquid crystal/polymer composites.

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.

Mesoscopic Modeling of a Highly-Ordered Sanidic Polymer Mesophase and Comparison With Experimental Data

Board-shaped polymers form sanidic mesophases: assemblies of parallel lamellae of stacked polymer backbones separated by disordered side chains. Sanidics vary significantly with respect to polymer order inside their lamellae, making them “stepping stones” toward the crystalline state. Therefore, they are potentially interesting for studying crystallization and technological applications. Building on earlier mesoscopic models of the most disordered sanidics Σ_d, we focus on the other extreme, near-crystalline order, and develop a generic model that captures a highly ordered Σ_r mesophase. Polymers are described by generic hindered-rotation chains. Anisotropic nonbonded potentials, with strengths comparable to the thermal energy, mimic board-like monomer shapes. Lamellae equilibrated with Monte Carlo simulations, for a broad range of model parameters, have intralamellar order typical for Σ_r mesophases: periodically stacked polymers that are mutually registered along their backbones. Our mesophase shows registration on both monomer and chain levels. We calculate scattering patterns and compare with data published for highly ordered sanidic mesophases of two different polymers: polyesters and polypeptoids. Most of the generic structural features that were identified in these experiments are present in our model. However, our mesophase has correlations between chains located in different lamellae and is therefore closer to the crystalline state than the experimental samples.

Self-assembly of freely-rotating polydisperse cuboids: unveiling the boundaries of the biaxial nematic phase

Colloidal cuboids have the potential to self-assemble into biaxial liquid crystal phases, which exhibit two independent optical axes. Over the last few decades, several theoretical works have predicted the existence of a wide region of the phase diagram where the biaxial nematic phase would be stable, but imposed rather strong constraints on the particle rotational degrees of freedom. In this work, we use molecular simulation to investigate the impact of size dispersity on the phase behaviour of freely-rotating hard cuboids, here modelled as self-dual-shaped nanoboards. This peculiar anisotropy, exactly in between the oblate and prolate geometry, has been proposed as the most appropriate to promote phase biaxiality. We observe that size dispersity radically changes the phase behaviour of monodisperse systems and leads to the formation of an elusive biaxial nematic phase, being found in a large region of the packing fraction vs. polydispersity phase diagram. Although our results confirm the tendencies reported in past experimental observations on colloidal dispersions of slightly prolate goethite particles, they cannot reproduce the direct isotropic-to-biaxial nematic phase transition observed in these experiments.

General patchy ellipsoidal particle model for the aggregation behaviors of shape- and/or surface-anisotropic building blocks

We present a general patchy ellipsoidal particle model suitable for conducting dynamics simulations of the aggregation behaviors of various shape- and/or surface-anisotropic colloids, especially patchy ellipsoids with continuously variable shape and tunable patchiness. To achieve higher computational efficiency in dynamics simulations, we employ a multi-GPU acceleration technique based on a domain decomposition algorithm. The validation and performance evaluation of this GPU-assisted model are performed by simulating several typical benchmark systems of non-patchy and patchy ellipsoids. Given the generality and efficiency of our GPU-assisted patchy ellipsoidal particle model, it will provide a highly feasible dynamics simulation framework to investigate the aggregation behaviors of anisotropic soft matter systems comprised of shape- and/or surface-anisotropic building blocks.