Tuning Coulombic interactions to stabilize nematic and smectic ionic liquid crystal phases in mixtures of charged soft ellipsoids and spheres

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.

Phase coexistence in [C22/C1MIm]+[NO3]- ionic-liquid mixtures and first-order phase transitions from homogeneous liquid to smectic B by varying the cation ratio

We perform molecular dynamics simulations to investigate the transition processes of [C22/C1MIm]+[NO3]- binary mixtures by varying the cation ratio of C22 to C1 at a fixed temperature of 400 K. The cation ratio is tuned by ranging C22 percentage from 0% to 100% with a fixed number of 4096 total simulated ion pairs. Our simulated-annealing results indicate that, at 400 K, pure C1 is a homogeneous liquid whilst pure C22 is an ionic liquid crystal (ILC) of smectic-B (SmB) type. With increasing C22 percentage, the system goes through a first-order phase transition from homogeneous liquid to nano-fragment liquid in the range from 15% to 17.5%, during which some of the individual cationic alkyl side chains locally aggregate to form small bundles "floating" in the polar "solvent" composed of anions and cationic head groups. Although the side chains in each bundle are parallelly aligned, the bundles distribute randomly without a global orientation. As the C22 percentage further increases, another first-order phase transition occurs to bring the system into the SmB ILC phase. Particularly, when the C22 percentage is in the range from 45% to 50%, the SmB phase coexists with the liquid phase containing both individual and bundled alkyl side chains.

Phase behaviour of mixtures of charged soft disks and spheres

We have investigated the phase behaviour of mixtures of soft disks (Gay-Berne oblate ellipsoids, GB) and soft spheres (Lennard-Jones, LJ) with opposite charge as a model of ionic liquid crystals and colloidal suspensions. We have used constant volume Molecular Dynamics simulations and fixed the stoichiometry of the mixture in order to have electroneutrality; three systems have been selected GB : LJ = 1 : 2, GB : LJ = 1 : 1 and GB : LJ = 2 : 1. For each system we have selected three values of the scaled point charge q* of the GB particles, namely 0.5, 1.0 and 2.0 (and a corresponding negative scaled charge of the LJ particles that depends on the stoichiometric ratio). We have found a very rich mesomorphism with the formation, as a function of the scaled temperature, of the isotropic phase, the discotic nematic phase, the hexagonal columnar phase and crystal phases. While the structure of the high temperature phases was similar in all systems, the hexagonal columnar phases exhibited a highly variable morphology depending on the scaled charge and stoichiometry. On the one hand, GB : LJ = 1 : 2 systems form lamellar structures, akin to smectic phases, with an alternation of layers of disks (exhibiting an hexagonal columnar phase) and layers of LJ particles (in the isotropic phase). On the other hand, for the 2 : 1 stoichiometry we observe the formation of a frustrated hexagonal columnar phase with an alternating tilt direction of the molecular axis. We rationalize these findings based on the structure of the neutral ion pair dominating the behaviour at low temperature and high charge.

Towards Nematic Phases in Ionic Liquid Crystals - A Simulation Study

Ionic liquid crystals (ILCs) are soft matter materials with broad liquid crystalline phases and intrinsic electric conductivity. They typically consist of a rod-shaped mesogenic ion and a smaller spherical counter-ion. Their mesomorphic properties can be easily tuned by exchanging the counter ion. ILCs show a strong tendency to form smectic A phases due to the segregation of ionic and the non-ionic molecular segments. Nematic phases are therefore extremely rare in ILCs and the question of why nematic phases are so exceptional in existing ILCs, and how nematic ILCs might be obtained in the future is of vital interest for both the fundamental understanding and the potential applications of ILCs. Here, we present the result of a simulation study, which highlights the crucial role of the location of the ionic charge on the rod-like mesogenic ions in the phase behaviour of ILCs. We find that shifting the charge from the ends towards the centre of the mesogenic ion destabilizes the liquid crystalline state and induces a change from smectic A to nematic phases.

Current Topics in Ionic Liquid Crystals

Ionic liquid crystals (ILCs), that is, ionic liquids exhibiting mesomorphism, liquid crystalline phases, and anisotropic properties, have received intense attention in the past years. Among others, this is due to their special properties arising from the combination of properties stemming from ionic liquids and from liquid crystalline arrangements. Besides interesting fundamental aspects, ILCs have been claimed to have tremendous application potential that again arises from the combination of properties and architectures that are not accessible otherwise, or at least not accessible easily by other strategies. The current review highlights recent developments in ILC research, starting with some key fundamental aspects. Further subjects covered include the synthesis and variations of modern ILCs, including the specific tuning of their mesomorphic behavior. The review concludes with reflections on some applications that may be within reach for ILCs and finally highlights a few key challenges that must be overcome prior and during true commercialization of ILCs.

Coarse-grained simulations of ionic liquid materials: from monomeric ionic liquids to ionic liquid crystals and polymeric ionic liquids

Ionic liquid (IL) materials are promising electrolytes with striking physicochemical properties for energy and environmental applications. Heterogeneous structures and transport quantities of monomeric and polymeric ILs are intrinsically intercorrelated and span multiple spatiotemporal scales, which is more feasible for coarse-grained (CG) simulations than atomistic modelling. Herein we constructed a novel CG model for ethyl-imidazolium tetrafluoroborate ILs with varied cation alkyl chains ranging from C₂ to C₂₀, and the interaction parameters were validated against representative static and dynamic properties that were obtained from atomistic reference simulations and experimental characterizations at relevant thermodynamic states. This CG model was extended to study thermotropic phase behaviors of monomeric ILs and to explore ion association structures and ion transport quantities in polymeric ILs with different architectures. A systematic analysis of structural and dynamical quantities identifies an evolution of liquid morphology from homogeneous to nanosegregated structures and then a smectic mesomorphism via a gradual lengthening of cation alkyl chains, and thereafter a distinct structural transition characterized by a monotonic decrease in orientational and translational order parameters in a sequential heating cascade. Backbone and pendant polymeric ILs exhibit evident anion association structures with cation monomers and polymer chains, and striking intra- and interchain coordinations between cation monomers owing to an intrinsic polymer architecture effect. Such a peculiar ion pairing association leads to a progressive increase in anion intrachain hopping probabilities, and a concomitant decrease in anion interchain hopping events with a gradual lengthening of polymeric ILs. The anion diffusivities in polymeric ILs are intrinsically correlated with ion pairing association lifetimes and ion structural relaxation times via a universal power law correlation D ∼ τ^-1, irrespective of polymer architectures.

The role of counterions in ionic liquid crystals

Previous theoretical studies of calamitic (i.e., rod-like) ionic liquid crystals (ILCs) based on an effective one-species model led to indications of a novel smectic-A phase with a layer spacing being much larger than the length of the mesogenic (i.e., liquid-crystal forming) ions. In order to rule out the possibility that this wide smectic-A phase is merely an artifact caused by the one-species approximation, we investigate an extension that accounts explicitly for cations and anions in ILCs. Our present findings, obtained by grand canonical Monte Carlo simulations, show that the phase transitions between the isotropic and the smectic-A phases of the cation-anion system are in qualitative agreement with the effective one-species model used in the preceding studies. In particular, for ILCs with mesogens (i.e., liquid-crystal forming species) carrying charged sites at their tips, the wide smectic-A phase forms, at low temperatures and within an intermediate density range, in between the isotropic and hexagonal crystal phases. We find that in the ordinary smectic-A phase, the spatial distribution of the counterions of the mesogens is approximately uniform, whereas in the wide smectic-A phase, the small counterions accumulate in between the smectic layers. Due to this phenomenology, the wide smectic-A phase could be interesting for applications, which hinge on the presence of conductivity channels for mobile ions.

Dynamics and Pretransitional Effects in C60 Fullerene Nanoparticles and Liquid Crystalline Dodecylcyanobiphenyl (12CB) Hybrid System

The report shows the strong impact of fullerene C60 nanoparticles on phase transitions and complex dynamics of rod-like liquid crystal dodecylcyanobiphenyl (12CB), within the limit of small concentrations. Studies were carried out using broadband dielectric spectroscopy (BDS) via the analysis of temperature dependences of the dielectric constant, the maximum of the primary loss curve, and relaxation times. They revealed a strong impact of nanoparticles, leading to a ~20% change of dielectric constant even at x = 0.05% of C60 fullerene. The application of the derivative-based and distortion-sensitive analysis showed that pretransitional effects dominate in the isotropic liquid phase up to 65 K above the clearing temperature and in the whole Smectic A mesophase. The impact of nanoparticles on the pretransitional anomaly appearance is notable for the smectic-solid phase transition. The fragility-based analysis of relaxation times revealed the universal pattern of its temperature changes, associated with scaling via the "mixed" ("activated" and "critical") relation. Phase behavior and dynamics of tested systems are discussed within the extended Landau-de Gennes-Ginzburg mesoscopic approach.

Microstructural and Dynamical Heterogeneities in Ionic Liquids

Ionic liquids (ILs) are a special category of molten salts solely composed of ions with varied molecular symmetry and charge delocalization. The versatility in combining varied cation-anion moieties and in functionalizing ions with different atoms and molecular groups contributes to their peculiar interactions ranging from weak isotropic associations to strong, specific, and anisotropic forces. A delicate interplay among intra- and intermolecular interactions facilitates the formation of heterogeneous microstructures and liquid morphologies, which further contributes to their striking dynamical properties. Microstructural and dynamical heterogeneities of ILs lead to their multifaceted properties described by an inherent designer feature, which makes ILs important candidates for novel solvents, electrolytes, and functional materials in academia and industrial applications. Due to a massive number of combinations of ion pairs with ion species having distinct molecular structures and IL mixtures containing varied molecular solvents, a comprehensive understanding of their hierarchical structural and dynamical quantities is of great significance for a rational selection of ILs with appropriate properties and thereafter advancing their macroscopic functionalities in applications. In this review, we comprehensively trace recent advances in understanding delicate interplay of strong and weak interactions that underpin their complex phase behaviors with a particular emphasis on understanding heterogeneous microstructures and dynamics of ILs in bulk liquids, in mixtures with cosolvents, and in interfacial regions.

Influence of the ion size on the stability of the smectic phase of ionic liquid crystals

The thermotropic phase behavior of ionic liquids and ionic liquid crystals based on novel N-alkyl-3-methylpyridinium halides, trihalides and dichloroiodates was experimentally studied by polarized optical spectroscopy (POM) and differential scanning calorimetry (DSC) as well as by molecular dynamics (MD) simulation. In the experiments, the existence and thermal range of stability of the smectic phase of these ionic liquid crystals are found to strongly depend on the volume ratio between the cation and anion, that is their relative size. Only compounds with a relatively large volume ratio of the cation to anion, i.e., those with longer cationic alkyl chains and monoatomic halide anions, have a stable smectic A phase. Both melting points and clearing points increase with such a ratio. The MD simulation results qualitatively agree very well with the experimental data and provide molecular details which can explain the experimentally observed phenomena: the stronger van der Waals interactions from the longer alkyl chains and the stronger electrostatic interactions from the smaller anions with a higher charge density increase the stability of both the crystal phase and the smectic phase; this also prevents the ionic layers from easily mixing with the hydrophobic regions, a mechanism that ultimately leads to a nanosegregated isotropic liquid phase.

A coarse-grained model of ionic liquid crystals: the effect of stoichiometry on the stability of the ionic nematic phase

We have investigated, by means of molecular dynamics simulations, the phase behaviour of mixtures of charged ellipsoidal Gay-Berne (GB) particles and spherical Lennard-Jones (LJ) particles, as a coarse-grained model of ionic liquid crystals (ILCs). The anisotropic GB particles represent cations usually found in ILCs, for example, pyridinium or bipyridinium salts, while the spherical LJ particles are taken as a model of anions like common halides, hexafluorophosphate and tetrafluoroborate. Here we have focused our attention on the effect of the stoichiometry of the system (that is, the GB : LJ ratio n : m in the salt formula [GB]n[LJ]m) on the stability and thermal range of the ionic liquid crystal phases formed, with special attention to the ionic nematic phase. To isolate the stoichiometry effect, a comparison of four different systems with GB : LJ ratios of 1 : 3, 1 : 2, 1 : 1 and 2 : 1 is made by keeping the packing fraction and the charge of the minor component fixed. Our results suggest a way to improve the stability of the ionic nematic phase by enhancing the anisotropic van der Waals interaction compared to the Coulomb interaction, and by increasing the proportion of anisotropic particles in the mixture.

A combined LX-NMR and molecular dynamics investigation of the bulk and local structure of ionic liquid crystals

The unique power of NMR spectroscopy in anisotropic media (LX-NMR) as a tool to obtain local and bulk structural information, combined with the effectiveness of molecular dynamics simulations at the atomistic level, shows very attractive potentialities for the study of interesting, even though still poorly understood, materials such as Ionic Liquid Crystals (ILCs). In this work, we focused our attention, in particular, on the orientational ordering of two mesophases: 1-dodecyl-3-methylimidazolium chloride, [C12C1im]Cl, and 1-dodecyl-3-methylimidazolium tetrafluoroborate, [C12C1im][BF4]. Both ILCs were studied by a 2H NMR direct investigation of the molecules forming the phases, suitably deuterated, and by 1H NMR spectroscopy, using the small rigid probe-solutes 1,4-dichlorobenzene (DCB), dissolved in [C12C1im][BF4] and [C12C1im]Cl, and 1,4-dibromobenzene (DBB) dissolved in [C12C1im][BF4], to probe the local, internal structure and organization of the mesophases. The experimental results were then compared with the predictions, by atomistic MD simulations, of the structure of the smectic phase of the two salts, at two selected temperatures, containing a single DCB molecule as a probe. The MD simulations show that the DCB solute is distributed only within the hydrophobic layers of the ILC. Orientational order parameters of the imidazolium cations and of the DCB molecule were obtained and compared with the experiments, showing a general good agreement and allowing a deeper understanding of the microscopic structure of the systems.

Interface structures in ionic liquid crystals

Ionic liquid crystals (ILCs) are anisotropic mesogenic molecules which additionally carry charges. This combination gives rise to a complex interplay of the underlying (anisotropic) contributions to the pair interactions. It promises interesting and distinctive structural and orientational properties to arise in systems of ILCs, combining properties of liquid crystals and ionic liquids. While previous theoretical studies have focused on the phase behavior of ILCs and the structure of the respective bulk phases, in the present study we provide new results, obtained within density functional theory, concerning (planar) free interfaces between an isotropic liquid L and two types of smectic-A phases (SA or SAW). We discuss the structural and orientational properties of these interfaces in terms of the packing fraction profile η(r) and the orientational order parameter profile S2(r) concerning the tilt angle α between the (bulk) smectic layer normal and the interface normal. The asymptotic decay of η(r) and of S2(r) towards their values in the isotropic bulk is discussed, too.

Effect of Central Longitudinal Dipole Interactions on Chiral Liquid-Crystal Phases

Monte Carlo simulations of chiral liquid-crystals, represented by a simple coarse-grained chiral Gay–Berne model, were performed to investigate the effect of central longitudinal dipole interactions on phase behavior. A systematic analysis of the structural properties and phase behavior of both achiral and chiral systems, with dipole interactions, reveals differing effects; strong dipole interactions enhance the formation of layered structures; however, chiral interactions may prevent the formation of such phases under certain conditions. We also observed a short-ranged smectic structure within the cholesteric phases with strong dipole interactions. This constitutes possible evidence of presmectic ordering and/or the existence of chiral line liquid phases, which have previously been observed in X-ray experiments to occur between the smectic twisted grain boundary and cholesteric phases. These results provide a systematic understanding of how the phase behavior of chiral liquid-crystals changes when alterations are made to the strength of dipole interactions.

Molecular dynamics study of mesophase transitions upon annealing of imidazolium-based ionic liquids with long-alkyl chains

Molecular dynamics simulations are performed on a 1-dodecyl-3-methylimidazolium hexafluorophosphate ([C12mim][PF6]) ionic liquid using a united-atom model. The ionic liquid exhibits second step relaxation at temperatures below a crossover point, where the diffusion coefficient shows an Arrhenius to non-Arrhenius transition. Annealing below this crossover temperature makes an isotropic to mesophase transition, where the smectic A (SmA) phase or crystal-like smectic B (SmB) phase forms. Hundreds of nanoseconds are required for completing these transitions. A normal diffusion process is found for anions along the layer-normal and -lateral directions in the SmA phase, but only in the lateral directions in the SmB phase. We find a preserved orientational order for the imidazolium-ring rotational and the alkyl-chain reorientational dynamics in both of the smectic phases.

Charge Transport and Phase Behavior of Imidazolium-Based Ionic Liquid Crystals from Fully Atomistic Simulations

Ionic liquid crystals occupy an intriguing middle ground between room-temperature ionic liquids and mesostructured liquid crystals. Here, we examine a non-polarizable, fully atomistic model of the 1-alkyl-3-methylimidazolium nitrate family using molecular dynamics in the constant pressure-constant temperature ensemble. These materials exhibit a distinct "smectic" liquid phase, characterized by layers formed by the molecules, which separate the ionic and aliphatic moieties. In particular, we discuss the implications this layering may have for electrolyte applications.

Ion pairing in 1-butyl-3-methylpyridinium halide ionic liquids studied using NMR and DFT calculations

Experimental and calculated ¹H, ¹³C and ¹⁵N NMR data of bulk 1-butyl-3-methylimidazolium halides ionic liquids provide key insights on their ion pairing.

Effect of counter ions on the mesogenic ionic N-phenylpyridiniums

Two series of ionic liquid crystals (ILCs) derived from N-phenylpyridinium 1–2 were prepared and their mesomorphic properties were investigated.