Spatial Dependence of Non-Gaussian Diffusion of Nanoparticles in Free-Standing Thin Polymer Films

The effects of defects on the transport mechanisms of lithium ions in organic ionic plastic crystals

Organic ionic plastic crystals (OIPCs) consist of molecular ions of which interactions are strong enough to maintain crystalline order but are weak enough to allow the rotations of the molecular ions at sufficiently high temperatures. When defects such as Schottky vacancies and grain boundaries are introduced into OIPCs, the defects facilitate the transport of dopants such as Li+ ions, for which OIPCs are considered as strong candidates for solid electrolytes. The transport mechanism of dopant ions in OIPCs with defects, however, remains elusive at a molecular level partly because it is hard in experiments to track the dopant ions and control the types of defects systematically. In this work, we perform molecular dynamics simulations for 1,3-dimethylimidazolium hexafluorophosphate ([MMIM][PF6]) OIPCs with Li+ ions doped and show that the transport mechanism of Li+ ions depends on the types and concentrations of defects. A high concentration of Schottky vacancies enhance the overall ion conduction, but decrease the transference number. The transference numbers of Li+ ions in [MMIM][PF6] with grain boundaries are similar to that in [MMIM][PF6] with 0.78 mol% point vacancies. We also find that the transport of ions in OIPCs is strongly heterogeneous and the time scales of the dynamic heterogeneity of the ions are sensitive to the types of defects.

Translation-rotation decoupling of tracers reflects medium-range crystalline order in two-dimensional colloid glasses

The dynamic heterogeneity and the translation-rotation decoupling are the dynamic signatures of glasses and supercooled liquids. Whether and how the dynamic heterogeneity would relate to the local structure of glasses has been a puzzle for decades. In this work we perform molecular dynamics simulations for tracers in both two-dimensional polydisperse colloids (2DPC) and two-dimensional binary colloids (2DBC). In 2DPC glasses, hexatic local structures develop at low enough temperatures and grow quickly along with the dynamic correlation length of the 2DPC, which is well known as the medium-range crystalline order (MRCO). In 2DBC glasses, on the other hand, any explicit local structure has not been reported to grow significantly with the dynamic correlation length at low temperatures. We introduce two different types of tracers into colloidal systems: A diamond tracer that resembles the MRCO of 2DPC glasses and a square tracer that is dissimilar to any local structure of glasses. The translation-rotation decoupling of the diamond tracer in 2DPC glasses is much more significant than that of the square tracer in the same 2DPC glasses. On the other hand, such a tracer shape-dependence of the decoupling is not observed in 2DBC glasses where the local hexatic structure does not develop significantly. We introduce a shape-dependency parameter of the decoupling and find that the shape-dependency parameter grows along with the dynamic correlation length in 2DPC glasses but not in 2DBC glasses. This illustrates that the dynamic heterogeneity and the translation-rotation decoupling of tracers could reveal the local structure that develops in glasses.

Temperature Dependence of Conformational Relaxation of Poly(ethylene oxide) Melts

The time-temperature superposition (TTS) principle, employed extensively for the analysis of polymer dynamics, is based on the assumption that the different normal modes of polymer chains would experience identical temperature dependence. We aim to test the critical assumption for TTS principle by investigating poly(ethylene oxide) (PEO) melts, which have been considered excellent solid polyelectrolytes. In this work, we perform all-atom molecular dynamics simulations up to 300 ns at a range of temperatures for PEO melts. We find from our simulations that the conformations of strands of PEO chains in melts show ideal chain statistics when the strand consists of at least 10 monomers. At the temperature range of T= 400 to 300 K, the mean-square displacements (⟨Δr2(t)⟩) of the centers of mass of chains enter the Fickian regime, i.e., ⟨Δr2(t)⟩∼t1. On the other hand, ⟨Δr2(t)⟩ of the monomers of the chains scales as ⟨Δr2(t)⟩∼t1/2 at intermediate time scales as expected for the Rouse model. We investigate various relaxation modes of the polymer chains and their relaxation times (τn), by calculating for each strand of n monomers. Interestingly, different normal modes of the PEO chains experience identical temperature dependence, thus indicating that the TTS principle would hold for the given temperature range.

Interdomain exchange and the flip-flop of cholesterol in ternary component lipid membranes and their effects on heterogeneous cholesterol diffusion

Cell membranes are heterogeneous with a variety of lipids, cholesterol, and proteins and are composed of domains of different compositions. Such heterogeneous environments make the transport of cholesterol complicated: cholesterol not only diffuses within a particular domain but also travels between domains. Cholesterol also flip-flops between upper and lower leaflets such that cholesterol may reside both within leaflets and in the central region between two leaflets. How the presence of multiple domains and the interdomain exchange of cholesterol would affect the cholesterol transport, however, remains elusive. In this study, therefore, we perform molecular dynamics simulations up to 100μs for ternary component lipid membranes, which consist of saturated lipids (dipalmitoylphosphatidylcholine, DPPC), unsaturated lipids (dilinoleylphosphatidylcholine, DIPC), and cholesterol. The ternary component membranes in our simulations form two domains readily: DPPC and DIPC domains. We find that the diffusion of cholesterol molecules is much more heterogeneous and non-Gaussian than expected for binary component lipid membranes of lipids and cholesterol. The non-Gaussian parameter of the cholesterol molecules is about four times larger in the ternary component lipid membranes than in the binary component lipid membranes. Such non-Gaussian and heterogeneous transport of cholesterol arises from the interplay among the interdomain kinetics, the different diffusivity of cholesterol in different domains, and the flip-flop of cholesterol. This suggests that in cell membranes that consist of various domains and proteins, the cholesterol transport can be very heterogeneous. We also find that the mechanism of the interdomain exchange differs for different domains: cholesterol tends to exit the DIPC domain along the central region of the membrane for the DIPC-to-DPPC transition, while the cholesterol is likely to exit the DPPC domain within the membrane leaflet for the DPPC-to-DIPC transition. Also, the interdomain exchange kinetics of cholesterol for the DPPC-to-DIPC transition is up to 7.9 times slower than the DIPC-to-DPPC transition.

The effects of vacancies and their mobility on the dynamic heterogeneity in 1,3-dimethylimidazolium hexafluorophosphate organic ionic plastic crystals

Organic ionic plastic crystals (OIPCs) are the crystals of electrolytes with a long-range translational order. The rotational modes of ions in OIPCs are, however, activated even in solid phases such that the diffusion of dopants such as lithium ions may be facilitated. OIPCs have been, therefore, considered as good candidates for solid electrolytes. Recent experiments and theoretical studies have suggested that both the translational and the rotational diffusion of ions are quite heterogeneous: the diffusion of some ions are quite fast while other ions of the same kind hardly diffuse, either rotationally or translationally. Such dynamic heterogeneity would be a key to the transport mechanism of dopants in solid state electrolytes. In this work, we investigate the effects of defects on the dynamic heterogeneity of OIPCs. We perform atomistic molecular dynamics simulation of 1,3-dimethylimidazolium hexafluorophosphate ([MMIM][PF6]) with a pair of cation and anion vacancies. At low temperature, vacancies undergo hopping motions toward each other and form a charge-neutral cluster. At high temperature, two vacancies act like a loosely bonded molecule and diffuse together via hopping motions. We find that the translational diffusion of ions is correlated strongly with the vacancy diffusion and becomes heterogeneous when the vacancies hop. The rotation of ions also becomes active when the ions are close to vacancies such that the rotational dynamic heterogeneity strengthens.