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

Accessing the electronic structure of liquid crystalline semiconductors with bottom-up electronic coarse-graining

Understanding the relationship between multiscale morphology and electronic structure is a grand challenge for semiconducting soft materials. Computational studies aimed at characterizing these relationships require the complex integration of quantum-chemical (QC) calculations, all-atom and coarse-grained (CG) molecular dynamics simulations, and back-mapping approaches. However, these methods pose substantial computational challenges that limit their application to the requisite length scales of soft material morphologies. Here, we demonstrate the bottom-up electronic coarse-graining (ECG) of morphology-dependent electronic structure in the liquid-crystal-forming semiconductor, 2-(4-methoxyphenyl)-7-octyl-benzothienobenzothiophene (BTBT). ECG is applied to construct density functional theory (DFT)-accurate valence band Hamiltonians of the isotropic and smectic liquid crystal (LC) phases using only the CG representation of BTBT. By bypassing the atomistic resolution and its prohibitive computational costs, ECG enables the first calculations of the morphology dependence of the electronic structure of charge carriers across LC phases at the ∼20 nm length scale, with robust statistical sampling. Kinetic Monte Carlo (kMC) simulations reveal a strong morphology dependence on zero-field charge mobility among different LC phases as well as the presence of two-molecule charge carriers that act as traps and hinder charge transport. We leverage these results to further evaluate the feasibility of developing mesoscopic, field-based ECG models in future works. The fully CG approach to electronic property predictions in LC semiconductors opens a new computational direction for designing electronic processes in soft materials at their characteristic length scales.

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.

Emergent spiral vortex of confined biased active particles

Confinement is known to have profound effects on the collective dynamics of many active systems. Here, we investigate a modeled active system in circular confinement consisting of biased active particles, where the direction of active force deviates a biased angle from the principle orientation of the anisotropic interaction. We find that such particles can spontaneously form a spiral vortex with two concentric and counter-rotating regions near the boundary. The emerged vortex can be measured by the vortex order parameter which shows nonmonotonic dependencies on both the biased angle and the strength of the anisotropic interaction. Our work can provide an understanding of such dynamic behaviors and enable different strategies for designing ordered collective behaviors.

Origin and structure of liquid crystalline Blue Phase III

We report here an off-lattice NVT molecular dynamics simulation study of a system of polar chiral ellipsoidal molecules, which spontaneously exhibits Blue Phase III (BPIII), considering coarse-grained attractive-repulsive pair interaction appropriate for anisotropic liquid crystal mesogens. We have observed that suitable selection of chiral and dipolar strengths not only gives rise to thermodynamically stable BPIII but novel Smectic and Bilayered BPIII as well. Further, we have demonstrated that the occurrence of BPIII and its layered counterparts depend crucially on molecular elongation.

Investigation of nematic to smectic phase transition and dynamical properties of strongly confined semiflexible polymers using Langevin dynamics

We investigated the nematic to smectic phase transition for strongly confined semiflexible polymer solutions in slit-like confinements using GPU-accelerated Langevin dynamics. We characterized the phase transitions from the nematic to smectic phases for semi-flexible polymer solutions as the polymer density increased. The dependence for the lyotropic nematic to smectic transition can be collapsed by scaling exponents between 0.2 and 0.3. The smectic C phase is found for all the cases with the polymer orientation director tilted with respect to smectic layer lateral alignment. As the chain rigidity increases, the transition density decreases for systems in which the polymer persistence length (P) to slit height (H) ratios are 1.25, 2.5, 3.75, 5 and 25. We also characterized the polymer dynamics for the isotropic-nematic-smectic transitions. The overall polymer diffusivity decreased steadily as the polymer density increased. We observed anomalous polymer diffusion along the nematic director near the isotropic-nematic transition, similar to previously reported behavior for nematic-forming ellipsoids. Polymer diffusivity decreased sharply by two orders of magnitude upon the nematic-smectic transition.

G.A.M.E.: GPU-accelerated mixture elucidator

GPU acceleration is useful in solving complex chemical information problems. Identifying unknown structures from the mass spectra of natural product mixtures has been a desirable yet unresolved issue in metabolomics. However, this elucidation process has been hampered by complex experimental data and the inability of instruments to completely separate different compounds. Fortunately, with current high-resolution mass spectrometry, one feasible strategy is to define this problem as extending a scaffold database with sidechains of different probabilities to match the high-resolution mass obtained from a high-resolution mass spectrum. By introducing a dynamic programming (DP) algorithm, it is possible to solve this NP-complete problem in pseudo-polynomial time. However, the running time of the DP algorithm grows by orders of magnitude as the number of mass decimal digits increases, thus limiting the boost in structural prediction capabilities. By harnessing the heavily parallel architecture of modern GPUs, we designed a “compute unified device architecture” (CUDA)-based GPU-accelerated mixture elucidator (G.A.M.E.) that considerably improves the performance of the DP, allowing up to five decimal digits for input mass data. As exemplified by four testing datasets with verified constitutions from natural products, G.A.M.E. allows for efficient and automatic structural elucidation of unknown mixtures for practical procedures.Graphical abstract