Search

In polydisperse systems, describing the structure and any structural order parameter (SOP) is not trivial as it varies with the number of species we use to describe the system, M . Depending on the degree of polydispersity, there is an optimum value of M = M0 where we show that the mutual information of the system increases. However, surprisingly the correlation between a recently proposed SOP and the dynamics is highest for M = 1. This effect increases with polydispersity. We find that the SOP at M = 1 is coupled with the particle size, {\sigma}, and this coupling increases with polydispersity and decreases with an increase in M . Careful analysis shows that at lower polydispersity the SOP is a good predictor of the dynamics. However, at higher polydispersity, the dynamics is strongly dependent on {\sigma}. Since the coupling between the SOP and {\sigma} is higher for M = 1 thus, it appears to be a better predictor of the dynamics. We also study the Vibrality an order parameter independent of structural information. Compared to SOP, at high polydispersity we find Vibrality to be a marginally better predictor of the dynamics. However, this high predictive power of Vibrality, which is not there at lower polydispersity, appears to be due to its stronger coupling with {\sigma}. Thus our study suggests that for systems with high polydispersity, the correlation of any order parameter and {\sigma} will affect the correlation between the order parameter and dynamics and need not project a generic predictive power of the order parameter., Comment: 14 pages, 41 figures, 1 table

We present the study of the structure dynamics correlation for systems interacting via attractive Lennard Jones and its repulsive counterpart, the WCA potentials. The structural order parameter (SOP) is the curvature of the microscopic mean field caging potential, termed the softness parameter. When calculated at a particle level, the SOP shows a distribution with particles having different degrees of softness. Although the two systems have similar structures, their average SOP is different. However, this difference is insufficient to explain the well known slowing down of the dynamics in LJ system at low temperatures. The slowing down of the dynamics can be explained in terms of a stronger coupling between the SOP and the dynamics. To understand the origin of this system specific coupling, we study the difference in the microscopic structure between the hard and soft particles by dividing the whole system into two communities. We find that compared to the WCA system for the LJ system, the communities structural differences are more significant and have a much stronger temperature dependence. Thus the study suggests that attractive interaction creates more structurally different communities. Even for the similar value of the average SOP, the difference in the local structure of the two communities is more comprehensive for the LJ system. Our study suggests that this broader difference in the structural communities is probably responsible for stronger coupling between the structure and dynamics. Thus the system specific structure dynamics correlation, which also leads to a faster slowing down in the dynamics, appears to have a structural origin., Comment: 12 pages, 37 figures

In this work, we perform a comparative study of the size dependence of diffusion of charged and neutral solutes in water. The neutral solute in water shows a nonmonotonicity in the size dependence of diffusion. This is usually connected to the well known Levitation effect where it is found that when solute diffuses through the transient solvent cages then for attractive solute-solvent interaction and for a particular size of the solute there is a force balance which leads to the maximum in diffusion. Similar maximum in diffusion of charged solutes has also been observed and connected to Levitation effect. However, earlier studies of ionic diffusion connects this nonmonotonicity to the interplay between hard sphere repulsion and Coulombic attraction. In this work, we show that although the size dependence of both charged and neutral solutes have a nonmonotonicity, there is a stark difference in their behaviour. For charged solute with increase in attraction the maximum shifts to higher solute sizes and has a lower value whereas for neutral solute it remains at the same place and has a higher value. We show by studying the ionic and non-ionic part of the potential that for larger solutes it is the nonionic part which dominates and for smaller solutes the ionic part and the is a transition between them. As the charge on the solute increases, this transition takes place at larger solute sizes which leads to the shift in the diffusivity maxima and reduction of the peak value. We show that although the charged solutes also explore the solvent cage even before we reach the size which Levitates due to Coulombic attraction the diffusion value drops. Thus the origin of diffusivity maxima in charged and neutral solute diffusion is different.

In a recent study we have found that for a large number of systems the configuration entropy at pair level, $S_{c2}$, which is primarily determined by the structural information, vanishes at the mode coupling transition temperature $T_{c}$. Thus it appears that the information of the transition temperature is embedded in the structure of the liquid. In order to investigate this we describe the dynamics of the system at the mean field level and using the concepts of the dynamical density function theory show that the dynamics depends only on the structure of the liquid. Thus this theory is similar in spirit to the microscopic MCT. However unlike microscopic MCT, which predicts a very high transition temperature, the present theory predicts a transition temperature which is similar to $T_{c}$. Thus our study reveals that the information of the mode coupling transition temperature is embedded in the structure of the liquid., Comment: 8 Pages, 3 Figs

As a liquid is progressively supercooled an intriguing weakening of the wavenumber ($q$) dependence of the structural relaxation time $\tau(q)$ in the large q limit is observed both in experiments and simulation studies. Neither the continuous Brownian diffusive dynamics nor the discontinuous activated events can alone explain the anomalous wavenumber dependence. Here we use our recently developed theory that unifies the mode coupling theory (MCT) for continuous dynamics and the random first order transition theory (RFOT) treatment of activated discontinuous motion as a nucleation (instanton) process, to understand the wavenumber dependence of density relaxation. The predicted smooth change in mechanism of relaxation from diffusive to activated, in the crossover regime, is wavevector dependent and is eventually responsible for sub-quadratic, almost linear, $q$ dependence of the relaxation time., Comment: 15 pages and 5 figures, minor changes in the content

In low temperature supercooled liquid, below the ideal mode coupling theory transition temperature, hopping and continuous diffusion are seen to coexist. We present a theory which incorporates interaction between the two processes and shows that hopping can induce continuous diffusion in the otherwise frozen liquid. Several universal features arise from nonlinear interactions between the continuous diffusive dynamics (described here by the mode coupling theory (MCT)) and the activated hopping (described here by the random first order transition theory). We apply the theory to real systems (Salol) to show that the theory correctly predicts the temperature dependence of the non-exponential stretching parameter, $\beta$, and the primary $\alpha$ relaxation timescale, $\tau$. The study explains why, even below the ergodic to non-ergodic transition, the dynamics is well described by MCT. The non-linear coupling between the two dynamical processes modifies the relaxation behavior of the structural relaxation from what would be predicted by a theory with a complete static Gaussian barrier distribution in a manner that may be described as a facilitation effect. Furthermore, the theory explains the observed variation of the stretching exponent $\beta$ with the fragility parameter, $D$., Comment: 8 pages, 4 figures

We present a theoretical analysis of the dynamic structure factor (DSF) of a liquid at and below the mode coupling critical temperature $T_c$, by developing a self-consistent theoretical treatment which includes the contributions both from continuous diffusion, described using general two coupling parameter ($F_{12}$) mode coupling theory (MCT), and from the activated hopping, described using the random first order transition (RFOT) theory, incorporating the effect of dynamical heterogeneity. The theory is valid over the whole temperature plane and shows correct limiting MCT like behavior above $T_{c}$ and goes over to the RFOT theory near the glass transition temperature, $T_{g}$. Between $T_{c}$ and $T_{g}$, the theory predicts that neither the continuous diffusion, described by pure mode coupling theory, nor the hopping motion alone suffices but both contribute to the dynamics while interacting with each other. We show that the interplay between the two contributions conspires to modify the relaxation behavior of the DSF from what would be predicted by a theory with a complete static Gaussian barrier distribution in a manner that may be described as a facilitation effect. Close to $T_c$, coupling between the short time part of MCT dynamics and hopping reduces the stretching given by the F$_{12}$-MCT theory significantly and accelerates structural relaxation. As the temperature is progressively lowered below $T_c$, the equations yield a crossover from MCT dominated regime to the hopping dominated regime. In the combined theory the dynamical heterogeneity is modified because the low barrier components interact with the MCT dynamics to enhance the relaxation rate below $T_c$ and reduces the stretching that would otherwise arise from an input static barrier height distribution., 7 pages, 4 figures