Molecular dynamics studies of ionically conducting glasses and ionic liquids: wave number dependence of intermediate scattering function

Dynamical heterogeneity is a key feature to characterize both acceleration and slowing down of the dynamics in interacting disordered materials. In the present work, the heterogeneous ion dynamics in both ionically conducting glass and in room temperature ionic liquids are characterized by the combination of the concepts of Lévy distribution and multifractality. Molecular dynamics simulation data of both systems are analyzed to obtain the fractional power law of the k-dependence of the dynamics, which implies the Lévy distribution of length scale. The multifractality of the motion and structures makes the system more complex. Both contributions in the dynamics become separable by using g(k,t) derived from the intermediate scattering function, F(s)(k,t). When the Lévy index obtained from F(s)(k,t) is combined with fractal dimension analysis of random walks and multifractal analysis, all the spatial exponent controlling both fast and slow dynamics are clarified. This analysis is generally applicable to other complex interacting systems and is deemed beneficial for understanding their dynamics.     

"Cooperativity blockage" in the mixed alkali effect as revealed by molecular-dynamics simulations of alkali metasilicate glass

The relaxation dynamics of a complex interacting system can be drastically changed when mixing with another component having different dynamics. In this work, we elucidate the effect of the less mobile guest ions on the dynamics of the more mobile host ions in mixed alkali glasses by molecular-dynamics (MD) simulations. One MD simulation was carried out on lithium metasilicate glass with the guest ions created by freezing some randomly chosen lithium ions at their initial locations at 700 K. A remarkable slowing down of the dynamics of the majority mobile Li ions was observed both in the self-part of the density-density correlation function, Fs(k,t), and in the mean-squared displacements. On the other hand, there is no significant change in the structure. The motion of the Li ions in the unadulterated Li metasilicate glass is dynamically heterogeneous. In the present work, the fast and slow ions were divided into two groups. The number of fast ions, which shows faster dynamics (Levy flight) facilitated by cooperative jumps, decreases considerably when small amount of Li ions are frozen. Consequently there is a large overall reduction of the mobility of the Li ions. The result is also in accordance with the experimental finding in mixed alkali silicate glasses that the most dramatic reduction of ionic conductivity occurs in the dilute foreign alkali limit. Similar suppression of the cooperative jumps is observed in the MD simulation data of mixed alkali system, LiKSiO3. Naturally, the effect found here is appropriately described as "cooperativity blockage." Slowing down of the motion of Li ions also was observed when a small number of oxygen atoms chosen at random were frozen. The effect is smaller than the case of freezing some the Li ions, but it is not negligible. The cooperativity blockage is also implemented by confining the Li metasilicate glass inside two parallel walls formed by freezing Li ions in the same metasilicate glass. Molecular-dynamics simulations were performed on the dynamics of the Li ions in the confined glass. Slowing down of the dynamics is largest near the wall and decreases monotonically with distance away from the wall.     

The mixed alkali effect in ionically conducting glasses revisited: a study by molecular dynamics simulation

When more than two kinds of mobile ions are mixed in ionic conducting glasses and crystals, there is a non-linear decrease of the transport coefficients of either type of ion. This phenomenon is known as the mixed mobile ion effect or Mixed Alkali Effect (MAE), and remains an unsolved problem. We use molecular dynamics simulation to study the complex ion dynamics in ionically conducting glasses including the MAE. In the mixed alkali lithium-potassium silicate glasses and related systems, a distinct part of the van Hove functions reveals that jumps from one kind of site to another are suppressed. Although, consensus for the existence of preferential jump paths for each kind of mobile ions seems to have been reached amongst researchers, the role of network formers and the number of unoccupied ion sites remain controversial in explaining the MAE. In principle, these factors when incorporated into a theory can generate the MAE, but in reality they are not essential for a viable explanation of the ion dynamics and the MAE. Instead, dynamical heterogeneity and "cooperativity blockage" originating from ion-ion interaction and correlation are fundamental for the observed ion dynamics and the MAE. Suppression of long range motion with increased back-correlated motions is shown to be a cause of the large decrease of the diffusivity especially in dilute foreign alkali regions. Support for our conclusion also comes from the fact that these features of ion dynamics are common to other ionic conductors, which have no glassy networks, and yet they all exhibit the MAE.     

Fast and slow dynamics in single and mixed alkali silicate glasses

A molecular dynamics (MD) simulation for single and mixed alkali metasilicates has been performed. The mixed alkali effect (MAE) in the diffusion coefficients is reproduced. The motion of lithium ions in Li₂SiO₃ can be divided into slow (type A) and fast (type B) categories clearly in the glassy state. It turns out that the slow dynamics are caused by localized jump motions with a long tail of the waiting time distribution. On the other hand, fast dynamics of the lithium ions are caused by successive cooperative jump motions and are identified as Lévy flights. The MAE in LiKSiO₃ occurs due to mutual interception of jump paths of both kinds of mobile ions. In such paths, the fast component is considerably suppressed. Namely, the cooperative forward-correlated jump process is sensitive to the blockage of the jump paths. Propagation of the immobilization in the MAE is discussed in relation with the mechanism of the cooperative jump motion. The ratio of the fast and slow dynamics plays an important role to determine the transport properties in the glassy state.     

Thermodynamic scaling of α-relaxation time and viscosity stems from the Johari-Goldstein β-relaxation or the primitive relaxation of the coupling model

By now it is well established that the structural α-relaxation time, τ(α), of non-associated small molecular and polymeric glass-formers obey thermodynamic scaling. In other words, τ(α) is a function Φ of the product variable, ρ(γ)/T, where ρ is the density and T the temperature. The constant γ as well as the function, τ(α) = Φ(ρ(γ)/T), is material dependent. Actually this dependence of τ(α) on ρ(γ)/T originates from the dependence on the same product variable of the Johari-Goldstein β-relaxation time, τ(β), or the primitive relaxation time, τ(0), of the coupling model. To support this assertion, we give evidences from various sources itemized as follows. (1) The invariance of the relation between τ(α) and τ(β) or τ(0) to widely different combinations of pressure and temperature. (2) Experimental dielectric and viscosity data of glass-forming van der Waals liquids and polymer. (3) Molecular dynamics simulations of binary Lennard-Jones (LJ) models, the Lewis-Wahnstr?m model of ortho-terphenyl, 1,4 polybutadiene, a room temperature ionic liquid, 1-ethyl-3-methylimidazolium nitrate, and a molten salt 2Ca(NO(3))(2)·3KNO(3) (CKN). (4) Both diffusivity and structural relaxation time, as well as the breakdown of Stokes-Einstein relation in CKN obey thermodynamic scaling by ρ(γ)/T with the same γ. (5) In polymers, the chain normal mode relaxation time, τ(N), is another function of ρ(γ)/T with the same γ as segmental relaxation time τ(α). (6) While the data of τ(α) from simulations for the full LJ binary mixture obey very well the thermodynamic scaling, it is strongly violated when the LJ interaction potential is truncated beyond typical inter-particle distance, although in both cases the repulsive pair potentials coincide for some distances.