Molecular dynamics simulation of thermotropic bolaamphiphiles with a swallow-tail lateral chain: formation of cubic network phases

Self-Assembly of Lamellae-in-Lamellae by Double-Tail Cationic Surfactants

The molecular structures of surfactants play a pivotal role in influencing their self-assembly behaviors. In this work, using simulations and experiments, an unconventional hierarchically layered structure in the didodecyldimethylammonium bromide (DDAB)/water binary system: lamellae-in-lamellae is revealed, a new self-assembly structure in surfactant system. This self-assembly structure refers to a lamellar structure with a shorter periodic length (inner lamellae) embedded in a lamellar phase with a longer periodic length (outer lamellae). The normal vectors of these two lamellar regions orient perpendicularly. In addition, it is observed that this lamellar-in-lamellar phase disappears when the two tails of the cationic surfactants become longer. The formation of the lamellar-in-lamellar architecture arises from multiple interacting factors. The key element is that the short tails of the DDAB surfactants enhance hydrophilicity and rigidity, which facilitates the formation of the inner lamellae. Moreover, the lateral monolayer of the inner lamellae provides shielding from the water and prompts the formation of the outer lamellae. These findings indicate that molecular structures and flexibility can profoundly redirect the hierarchical self-assembly behaviors in amphiphilic systems. More broadly, this work presents a new strategy to deliberately program hierarchical nanomaterials by designing specific surfactant molecules to act as tunable scaffolds, reactors, and carriers.

Coarse-Grained Molecular Simulation of Bolapolyphiles with a Multident Lateral Chain: Formation and Structural Analysis of Cubic Network Phases

Bolapolyphiles constitute a versatile class of materials with a demonstrated potential to form a wide variety of complex ordered mesophases. In particular, cubic network phases (like the gyroid, primitive, and diamond phases) have been a target of many studies for their ability to create percolating 3D nanosized channels. In this study, molecular simulations are used to explore the phase behavior of bolapolyphiles containing a rigid rodlike core, associating hydrophilic core ends and a hydrophobic side chain with a multident architecture, i.e., where the branching pattern can vary from bident (two branches) to hexadent (six branches). Upon network phase formation, its skeleton is made up of "nodes" populated by the core ends and "struts" populated by the cores. It is shown that, by varying the side chain length, branching pattern, and attachment point to the core, one can alter the crowding around the cores and hence tune the nodal size and nodal valence (i.e., number of connecting struts) which lead to different types of network morphologies. For example, for a fixed total side chain length, having more branches generates a stronger crowding around the molecular core, driving them to form bundlelike domains with curvier interfaces that result in thinner struts. Also, attaching the lateral chain closer to one core end breaks the symmetry between the environments around the two core ends, leading to networks with bimodal nodal sizes. Importantly, since the characterization of (ordered or partially ordered) network phases is challenging given the potential incompatibilities between the simulation box size with the structure's space group periodic symmetry and the effect of morphological defects, a detailed framework is presented to analyze and fully characterize the unit cell parameters and structure factor of such systems.

Controlling Mirror Symmetry Breaking and Network Formation in Liquid Crystalline Cubic, Isotropic Liquid and Crystalline Phases of Benzil-Based Polycatenars

Spontaneous development of chirality in systems composed of achiral molecules is important for new routes to asymmetric synthesis, chiral superstructures and materials, as well as for the understanding of the mechanisms of emergence of prebiotic chirality. Herein, it is shown that the 4,4'-diphenylbenzil unit is a universal transiently chiral bent building block for the design of multi-chained (polycatenar) rod-like molecules capable of forming a wide variety of helically twisted network structures in the liquid, the liquid crystalline (LC) and the crystalline state. Single polar substituents at the apex of tricatenar molecules support the formation of the achiral (racemic) cubic double network phase with Ia 3 ‾ d symmetry and relatively small twist along the networks. The combination of an alkyl chain with fluorine substitution leads to the homogeneously chiral triple network phase with I23 space group, and in addition, provides a mirror symmetry broken liquid. Replacing F by Cl or Br further increases the twist, leading to a short pitch double gyroid Ia 3 ‾ d phase, which is achiral again. The effects of the structural variations on the network structures, either leading to achiral phases or chiral conglomerates are analyzed.

Thermal Stability of π-Conjugated n-Ethylene-Glycol-Terminated Quaterthiophene Oligomers: A Computational and Experimental Study

This work represents a joint computational and experimental study on a series of n-ethylene glycol (PEOn)-terminated quaterthiophene (4T) oligomers for 1 < n < 10 to elucidate their self-assembly behavior into a smectic-like lamellar phase. This study builds on an earlier study for n = 4 that showed that our model predictions were consistent with experimental data on the melting behavior and structure of the lamellar phase, with the latter consisting of crystal-like 4T domains and liquid-like PEO4 domains. The present study aims to understand how the length of the terminal PEOn chains modulates the disordering temperature of the lamellar phase and hence the relative stability of the ordered structure. A simplified bilayer model, where the 4T domains are not explicitly described, is put forward to efficiently estimate the disordering effect of the PEO domains with increasing n; this method is first validated by correctly predicting that layers of alkyl (PE)-capped 4T oligomers (for 1 < n < 10) stay ordered at room temperature. Both 4T-domain implicit and explicit model simulations reveal that the order-disorder temperature decreases with the length of the PEO capping chains, as the associated increase in conformational entropy drives a tendency toward disorder that overtakes the cohesive energy, keeping the ordered packing of the 4T domains.

Following the nucleation pathway from disordered liquid to gyroid mesophase

Mesophases have order intermediate between liquids and crystals and arise in systems with frustration, such as surfactants, block copolymers, and Janus nanoparticles. The gyroid mesophase contains two interpenetrated, nonintersecting chiral networks that give it properties useful for photonics. It is challenging to nucleate a gyroid from the liquid. Elucidating the reaction coordinate for gyroid nucleation could assist in designing additives that facilitate the formation of the mesophase. However, the complexity of the gyroid structure and the extreme weakness of the first-order liquid to gyroid transition make this a challenging quest. Here, we investigate the pathway and transition states for the nucleation of a gyroid from the liquid in molecular simulations with a mesogenic binary mixture. We find that the gyroid nuclei at the transition states have a large degree of positional disorder and are not compact, consistent with the low surface free energy of the liquid-gyroid interface. A combination of bond-order parameters for the minor component is best to describe the passage from liquid to gyroid, among those we consider. The committor analyses, however, show that this best coordinate is not perfect and suggests that accounting for the relative ordering of the two interpenetrated networks in infant nuclei, as well as for signatures of ordering in the major component of the mesophase, would improve the accuracy of the reaction coordinate for gyroid formation and its use to evaluate nucleation barriers. To our knowledge, this study is the first to investigate the reaction coordinate and critical nuclei for the formation of any mesophase from an amorphous phase.

Why Is Gyroid More Difficult to Nucleate from Disordered Liquids than Lamellar and Hexagonal Mesophases?

Block copolymers, surfactants, and biomolecules form lamellar, hexagonal, and gyroid mesophases. Across these systems, the nucleation of lamellar from the disordered liquid is the easiest and the nucleation of gyroid the most challenging. This poses the question of what are the factors that determine the rates of nucleation of the mesophases and whether they are controlled by the complexity of the structures or the thermodynamics of nucleation. Here, we use molecular simulations to investigate the nucleation and thermodynamics of lamellar, hexagonal, and gyroid in a binary mixture of particles that produces the same mesophases as those of surfactants and block copolymers. We demonstrate that a combination of averaged bond-order parameters q̅₂ and q̅₈ identifies and distinguishes the three mesophases. We use these parameters to track the microscopic process of nucleation of each mesophase and investigate the existence of heterogeneous nucleation (cross-nucleation) between mesophases. We estimate the surface tensions of the liquid/mesophase interfaces from nucleation rates using classical nucleation theory and find that they are comparable for the three mesophases with values that are about a third of those expected for liquid-crystal interfaces. The driving forces for nucleation, on the other hand, are quite different and increase in the order gyroid < hexagonal < lamellar at any temperature. We find that the nucleation rates of the mesophases follow the order of their driving forces. We conclude that the difficulty to nucleate the gyroid originates in its lower temperature of melting and extremely low entropy of melting compared to those of the hexagonal and lamellar mesophases.