We propose a novel model for a glass-forming liquid, which allows us to switch in a continuous manner from a standard three-dimensional liquid to a fully connected mean-field model. This is achieved by introducing k additional particle–particle interactions, which thus augments the effective number of neighbors of each particle. Our computer simulations of this system show that the structure of the liquid does not change with the introduction of these pseudo-neighbors and by means of analytical calculations, and we determine the structural properties related to these additional neighbors. We show that the relaxation dynamics of the system slows down very quickly with the increase in k and that the onset and the mode-coupling temperatures increase. The systems with high values of k follow the mode-coupling theory power law behavior for a larger temperature range compared to the ones with lower values of k. The dynamic susceptibility indicates that the dynamic heterogeneity decreases with the increase in k, whereas the non-Gaussian parameter is independent of it. Thus, we conclude that with the increase in the number of pseudo-neighbors, the system becomes more mean-field-like. By comparing our results with previous studies on mean-field-like systems, we come to the conclusion that the details of how the mean-field limit is approached are important since they can lead to different dynamical behavior in this limit.

We present a study of the dynamics of small solute particles in a solvent medium where the solute is much smaller in size, mimicking the diffusion of small particles in crowded environment. The solute exhibits Fickian diffusion arising from non-Gaussian van Hove correlation function. Our study shows that there are at least two possible origins of this non-Gaussian behaviour. The decoupling of the solute-solvent dynamics and the intermittency in the solute motion, the latter playing a dominant role. In the former scenario when averaged over time long enough to explore different solvent environments the dynamics recovers the Gaussian nature. In case of intermittent dynamics the non-Gaussianity remains even after long averaging and the Gaussian behaviour is obtained at a much longer time. Our study further shows that only for intermediate attractive solute-solvent interaction the dynamics of the solute is intermittent. The intermittency disappears for weaker or stronger attractions.

We present a comparative study of the glass forming ability of binary systems with varying composition, where the systems have similar global crystalline structure (CsCl+fcc). Biased Monte Carlo simulations using umbrella sampling technique show that the free energy cost to create a CsCl nucleus increases as the composition of the smaller particles is decreased. We find that systems with comparatively lower free energy cost to form CsCl nucleus exhibit more pronounced pre-crystalline demixing near the liquid/crystal interface. The structural frustration between the CsCl and fcc crystal demands this demixing. We show that closer to the equimolar mixture, the entropic penalty for demixing is lower and a glass forming system may crystallize when seeded with a nucleus. This entropic penalty as a function of composition shows a non-monotonic behaviour with a maximum at a composition similar to the well known Kob-Anderson (KA) model. Although the KA model shows the maximum entropic penalty and thus maximum frustration against CsCl formation, it also shows a strong tendency towards crystallization into fcc lattice of the larger "A" particles which can be explained from the study of the energetics. Thus for systems closer to the equimolar mixture although it is the requirement of demixing which provides their stability against crystallization, for KA model it is not demixing but slow dynamics and the presence of the "B" particles make it a good glass former. The locally favoured structure around "B" particles is quite similar to the CsCl structure and the incompatibility of CsCl and fcc hinders the fcc structure growth in the KA model. Although the glass forming binary systems studied here are quite similar, differing only in composition, we find that their glass forming ability cannot be attributed to a single phenomenon.

No abstract

In a system of N particles, with continuous size polydispersity, there exists an N(N − 1) number of partial structure factors, making it analytically less tractable. A common practice is to treat the system as an effective one component system, which is known to exhibit an artificial softening of the structure. The aim of this study is to describe the system in terms of M pseudospecies such that we can avoid this artificial softening but, at the same time, have a value of M ≪ N. We use potential energy and pair excess entropy to estimate an optimum number of species, M0. We then define the maximum width of polydispersity, Δσ0, that can be treated as a monodisperse system. We show that M0 depends on the degree and type of polydispersity and also on the nature of the interaction potential, whereas Δσ0 weakly depends on the type of polydispersity but shows a stronger dependence on the type of interaction potential. Systems with a softer interaction potential have a higher tolerance with respect to polydispersity. Interestingly, M0 is independent of system size, making this study more relevant for bigger systems. Our study reveals that even 1% polydispersity cannot be treated as an effective monodisperse system. Thus, while studying the role of polydispersity by using the structure of an effective one component system, care must be taken in decoupling the role of polydispersity from that of the artificial softening of the structure.

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.

Note: This paper is part of the JCP Special Topic on Slow Dynamics.