Shear Thinning and Nonlinear Structural Deformation of Ionic Liquids with Long Alkyl Chains Studied by Molecular Dynamics Simulation

Structural effects of water clusters on viscosity at high shear rates

In this study, we use molecular dynamics simulations of liquid water to investigate how shear thinning affects the viscosity of liquid water by structural changes of the hydrogen bond network. The effect of shear on viscosity can be divided into two parts: shear-induced destruction of the hydrogen bond network and the influence of the water structure on shear viscosity. First, strong shear destroys tetrahedral structures and thus reduces the connectivity of the hydrogen bond network. It is mainly because shear deformation, characterized by compression and expansion axes, respectively, triggers the destruction and formation of hydrogen bonds, resulting in anisotropic effects on water structures. At the same time, shear destroys large clusters and enhances the formation of small ones, resulting in a decrease in average cluster sizes. Second, the change of viscosity obeys a power law relationship with the change of hydrogen bond structures, highlighting a one-to-one correspondence between structure and property. Meanwhile, in order to explain why the structure affects viscosity, we define hydrogen-bond viscosity and find that the cooperative motion of the water structures can promote momentum transfer in the form of aggregations. Hydrogen-bond viscosity accounts for 5%-50% of the total viscosity. Our results elucidate that water structures are the important structural units to explain the change of water properties.

Elucidating the interplay of local and mesoscale ion dynamics and transport properties in aprotic ionic liquids

Ion dynamics and charge transport in 1-methyl-3-octylimidazolium ionic liquids with chloride, bromide, tetrafluoroborate, tricyanomethanide, hexafluorophosphate, triflate, tetrachloroaluminate, bis(trifluoromethylsulfonyl)imide, and heptachlorodialuminate anions are investigated by broadband dielectric spectroscopy, rheology, viscometry, and differential scanning calorimetry. A detailed analysis reveals an anion and temperature-dependent separation of characteristic molecular relaxation rates extracted from various representations of the dielectric spectra. The separation in rates extracted from the electric modulus and conductivity formalisms is interpreted as an experimental signature of significant heterogeneity in the local ion dynamics associated with the structural glass transition, viscosity, and dc ion conductivity. It is further found that the degree of dynamic heterogeneity correlates with the strengths of slow dielectric and mechanical relaxations previously attributed to the dynamics of mesoscale solvophobic aggregates. Increasing local dynamic heterogeneity correlates with an increase in the strength of the slow, aggregate dielectric relaxation and a decrease in the strength of the slow, aggregate mechanical relaxation. Accordingly, increasing local dynamic heterogeneity, brought about by change in temperature and/or cation/anion chemical structure, correlates with an increase in the static dielectric permittivities and a decrease in the contribution of aggregate dynamics to the zero-shear viscosities. The established correlation provides a new ability to distinguish between the influence of mesoscale aggregate shape/morphology versus local and mesoscale ion dynamics on the transport properties of ionic liquids.

Viscoelasticity of Low-Molecular-Weight Polyelectrolytes

Shear-thickening fluids that absorb the impact energy of high-velocity projectiles are of great interest for aerospace and body-armor applications. In such a frame, we investigate transient states of neat and aqueous polyelectrolytes (PE) having low molecular weights and containing poly([2-(methacryloyloxy)ethyl]trimethylammonium) as polycations and poly(acrylamide-co-acrylic acid) as polyanions. We compare results with those of bulk water. We employ nonequilibrium molecular dynamics to simulate oscillatory shear, mainly in the linear viscoelastic regime. We find that neat PE exhibits properties of a viscoelastic solid, whereas water and the aqueous mixture of PE conform to viscoelastic liquids with Maxwellian behavior at low angular frequencies. Terminal relaxation times are ∼0.499 and ∼1.385 ps for water and the aqueous mixture of PE, respectively. At high angular frequencies, storage moduli show anomalous behaviors that correspond to transitions between shear thinning and shear thickening in complex shear viscosities. The change in potential energy with the increase of the angular frequency is mainly driven by intramolecular interactions for neat PE, whereas short-range Coulomb interactions are the major contributions for water and the aqueous mixture of PE. Upon observation of the molecular configurations, only the local polyionic structure in the aqueous mixture of PE shows improvement when increasing the angular frequency, whereas the rest remains barely affected. Thus, the water structure in the aqueous mixture of PE allows the storage of energy elastically through the hydrogen-bond network at large angular frequencies, whereas the mechanical contribution of polyions weakens and fully vanishes at the beginning of shear thinning, explaining the superimposed data with data of bulk water. Our method and findings set the path for future molecular simulations in the nonlinear viscoelastic regime with more complex underlying molecular mechanisms.

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.