Role of the Electrostatic Interactions in the Stabilization of Ionic Liquid Crystals: Insights from Coarse-Grained MD Simulations of an Imidazolium Model

Coarse-Grained Simulations of Columnar Ionic Liquid Crystals: Comparison with Experiments

We simulate a homologous series of guanidinium-based columnar ionic liquid crystals (ILCs) using coarse-grained molecular dynamics (MD) simulations with the Martini force field. We systematically vary the length of alkyl side chains, ILC-n (n = 8, 12, 16), and compare our results with previous experimental findings. Experimentally, ILC-8 exhibits a narrow mesophase window and weak columnar order, while ILC-12 and ILC-16 display a broad mesophase window and high columnar order. The MD simulations show that ILC-8 forms a percolated structure, whereas the longer chain analogues self-assemble into columns, with columnar assembly becoming more prominent as the side chain length increases, in qualitative agreement with the experiments. Furthermore, the intercolumnar distance increases monotonically with increasing side chain length and decreases with increasing temperature. Finally, we find that the diffusion coefficient and ionic conductivity decrease substantially with increasing chain length, consistent with experimental observations. We attribute this decrease in mobility to the formation of hexagonally ordered columns, which restrict transport more than percolated networks.

Martini on the Rocks: Can a Coarse-Grained Force Field Model Crystals?

Computational chemistry is an important tool in numerous scientific disciplines, including drug discovery and structural biology. Coarse-grained models offer simple representations of molecular systems that enable simulations of large-scale systems. Because there has been an increase in the adoption of such models for simulations of biomolecular systems, critical evaluation is warranted. Here, the stability of the amyloid peptide and organic crystals is evaluated using the Martini 3 coarse-grained force field. The crystals change shape drastically during the simulations. Radial distribution functions show that the distance between backbone beads in β-sheets increases by ∼1 Å, breaking the crystals. The melting points of organic compounds are much too low in the Martini force field. This suggests that Martini 3 lacks the specific interactions needed to accurately simulate peptides or organic crystals without imposing artificial restraints. The problems may be exacerbated by the use of the 12-6 potential, suggesting that a softer potential could improve this model for crystal simulations.

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.

The effect of hydration on the stability of ionic liquid crystals: MD simulations of [C14C1im]Cl and [C14C1im]Cl·H2O

The thermal range of the stability of Ionic Liquid Crystal (ILC) phases of imidazolium ILCs, and the type of the mesophase itself are affected by several molecular structural features, the two prominent ones being the alkyl chain length and the counter-anion. Hydration is also very important: monohydrate samples of 1-alkyl-3-methylimidazolium halides have a higher clearing point and a wider thermal range of the stability of the ionic smectic phase, compared with the analogous anhydrous sample. To understand the reasons, at a microscopic level, for such increased stability due to hydration, we run classical Molecular Dynamics (MD) simulations of a typical ionic liquid crystal, 1-tetradecyl-3-methylimidazolium chloride, and of its monohydrate form. We tested a full-charge non-polarizable force field and a scaled-charge version having the total charge of the ions scaled by a factor of 0.80. Comparison of the structural and dynamic properties with available experimental data reveals that the scaling of the charge by a factor of 0.80 results in a good agreement between simulated and experimental data and it sheds light on the microscopic mechanism responsible for the increased stability of the monohydrated phase. A hydrogen-bond network between water and the chloride anion is established in the ionic layer which increases the stability of the ionic layer; this in turn increases the nano-segregation between the ionic and hydrophobic layers which eventually produce an increased order of the alkylic layer as well.

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.

Molecular-dynamics simulations on the mesophase transition induced by oscillatory shear in imidazolium-based ionic liquid crystals

Molecular dynamics simulations were performed on a 1-dodecyl-3-methylimidazolium hexafluorophosphate ([C₁₂mim][PF₆]) ionic liquid crystal (ILC) with the application of an oscillatory shear. We found that the oscillatory shear can both accelerate and suppress mesophase formation depending on shear amplitude. A small amplitude shear can speed up the mesophase transition dynamics and result in a more ordered mesomorphic structure than that without shear, i.e., an effect of accelerated aging. The mesophase is destabilized when the shear amplitude is large enough, resulting in a smectic A (SmA) to liquid or a smectic B (SmB) to SmA transition, with the mesophase behaviour summarized in an out-of-equilibrium phase diagram. Inside the layer plane a medium-range hexatic order was observed, with the correlation length extending to several nanometres in the shear-induced SmA phase. We rationalize the nonequilibrium mesophase behaviour from the rheology of isotropic liquids, finding a temperature-independent critical relaxation time for the mesophase transition in the translational or rotational dynamics. This finding can be used to predict the mesophase behaviour in the sheared ILCs from the rheology of isotropic liquids.

Ion conformation and orientational order in a dicationic ionic liquid crystal studied by solid-state nuclear magnetic resonance spectroscopy

Ionic liquids crystals belong to a special class of ionic liquids that exhibit thermotropic liquid-crystalline behavior. Recently, dicationic ionic liquid crystals have been reported with a cation containing two single-charged ions covalently linked by a spacer. In ionic liquid crystals, electrostatic and hydrogen bonding interactions in ionic sublayer and van der Waals interaction in hydrophobic domains are the main forces contributing to the mesophase stabilization and determining the molecular orientational order and conformation. How these properties in dicationic materials are compared to those in conventional monocationic analogs? We address this question using a combination of advanced NMR methods and DFT analysis. Dicationic salt 3,3′-(1,6-hexanediyl)bis(1-dodecylimidazolium)dibromide was studied. Local bond order parameters of flexible alkyl side chains, linker chain, and alignment of rigid polar groups were analyzed. The dynamic spacer effectively “decouples” the motion of two ionic moieties. Hence, local order and alignment in dicationic mesophase were similar to those in analogous single-chain monocationic salts. Bond order parameters in the side chains in the dicationic smectic phase were found consistently lower compared to double-chain monocationic analogs, suggesting decreasing contribution of van der Waals forces. Overall dication reorientation in the smectic phase was characterized by low values of orientational order parameter S. With increased interaction energy in the polar domain the layered structure is stabilized despite less ordered dications. The results emphasized the trends in the orientational order in ionic liquid crystals and contributed to a better understanding of interparticle interactions driving smectic assembly in this and analogous ionic mesogens.

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.

Phase Behaviors of Ionic Liquids Heating from Different Crystal Polymorphs toward the Same Smectic-A Ionic Liquid Crystal by Molecular Dynamics Simulation

Five distinct crystal structures, based on experimental data or constructed manually, of ionic liquid [C14Mim][NO3] were heated in NPT molecular dynamics simulations under the same pressure such that they melted into the liquid crystal (LC) phase and then into the liquid phase. It was found that the more entropy-favored structure had a higher solid-LC transition temperature: Before the transition into the LC, all systems had to go through a metastable state with the side chains almost perpendicular to the polar layers. All those crystals finally melted into the same smectic-A LC structure irrelevant of the initial crystal structure.

Neutron scattering studies on short- and long-range layer structures and related dynamics in imidazolium-based ionic liquids

Alkyl-methyl-imidazolium ionic liquids CnmimX (n: alkyl-carbon number, X: anion) have short-range layer structures consisting of ionic and neutral (alkylchain) domains. To investigate the temperature dependences of the interlayer, interionic group, and inter-alkylchain correlations, we have measured the neutron diffraction (ND) of C16mimPF₆, C9.5mimPF₆, and C8mimPF₆ in the temperature region from 4 K to 470 K. The quasielastic neutron scattering (QENS) of C16mimPF₆ was also measured to study the dynamics of each correlation. C16mimPF₆ shows a first-order transition between the liquid (L) and liquid crystalline (LC) phases at T_c = 394 K. C8mimPF₆ exhibits a glass transition at T_g = 200 K. C9.5mimPF₆, which is a 1:3 mixture between C8mimPF₆ and C10mimPF₆, has both transitions at T_c = 225 K and T_g = 203 K. In the ND experiments, all samples exhibit three peaks corresponding to the correlations mentioned above. The widths of the interlayer peak at ca. 0.2 Å^-1 changed drastically at the L-LC transitions, while the interionic peaks at ca. 1 Å^-1 exhibited a small jump at T_c. The peak position and area of the three peaks did not change much at the transition. The structural changes were minimal at T_g. The QENS experiments demonstrated that the relaxation time of the interlayer motion increased tenfold at T_c, while those of other motions were monotonous in the whole temperature region. The structural and dynamical changes mentioned above are characteristic of the L-LC transition in imidazolium-based ionic liquids.

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.

Surface induced smectic order in ionic liquids - an X-ray reflectivity study of [C22C1im]+[NTf2]. Phys Chem Chem Phys 2018;19:26651-26661. [PMID: 28960006 DOI: 10.1039/c7cp04852a]

Surface induced smectic order was found for the ionic liquid 1-methyl-3-docosylimidazolium bis(trifluoromethlysulfonyl)imide by X-ray reflectivity and grazing incidence scattering experiments. Near the free liquid surface, an ordered structure of alternating layers composed of polar and non-polar moieties is observed. This leads to an oscillatory interfacial profile perpendicular to the liquid surface with a periodicity of 3.7 nm. Small angle X-ray scattering and polarized light microscopy measurements suggest that the observed surface structure is related to fluctuations into a metastable liquid crystalline SmA₂ phase that was found by supercooling the bulk liquid. The observed surface ordering persists up to 157 °C, i.e. more than 88 K above the bulk melting temperature of 68.1 °C. Close to the bulk melting point, we find a thickness of the ordered layer of L = 30 nm. The dependency of L(τ) = Λ ln(τ/τ₁) vs. reduced temperature τ follows a logarithmic growth law. In agreement with theory, the pre-factor Λ is governed by the correlation length of the isotropic bulk phase.

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

We have investigated the effect of electrostatic interactions in mixtures of soft ellipsoids and spheres based on the well-known Gay-Berne (GB) and Lennard-Jones (LJ) potential, respectively. These model systems, in their original version, that is without any electrostatic charge, have been thoroughly investigated in the literature both as pure components and mixtures. In particular, mixtures of particles of different shapes, such as spheres and ellipsoids, tend to phase separate because of the excluded volume effects. Common ionic liquid crystals, based on imidazolium or other quaternary ammonium salts, are usually composed of roughly elongated (although flexible) cations and roughly spherical anions, that is, particles with a similar shape such as the GB and LJ models. Therefore, in this work, we present the results of molecular dynamics simulations of mixtures of positively charged GB and negatively charged LJ particles as models of ionic liquid crystals. Interestingly, by modulating the charge of the particles it is possible to stabilize isotropic, nematic, smectic and crystalline ionic phases. The relative weight of Coulomb (a radial, therefore isotropic interaction) and van der Waals (an anisotropic interaction) contributions is a key parameter to tune the stability of various mesophases.

Highly selective detection of nitrotoluene based on novel lanthanide-containing ionic liquids

Two novel rare-earth ionic liquids demonstrate high selectivity toward nitrotoluene in the presence of other aromatic compounds.