Search

A universal dynamical crossover temperature, Tcr, in glassy liquids, associated with the {\alpha}-\b{eta} bifurcation temperature, TB, has been observed in dielectric spectroscopy and other experiments. Tcr lies significantly above the glass transition temperature. Here, we introduce a new class of glass-forming liquids, binary mixtures of prolate and oblate ellipsoids. This model system exhibits sharp thermodynamic and dynamic anomalies, such as the specific heat jump during heating and a sharp variation in the thermal expansion coefficient around a temperature identified as the glass transition temperature, Tg. The same temperature is obtained from the fit of the calculated relaxation times to the Vogel-Fulcher-Tammann (VFT) form. As the temperature is lowered, the single peak rotational relaxation spectrum splits into two at a temperature TB significantly above the estimated Tg. Similar bifurcation is also observed in the distribution of short-to-intermediate time translational diffusion. Interrogation of the two peaks reveals a lower extent of dynamic heterogeneity in the population of the faster mode. We observe an unexpected appearance of a sharp peak in the product of rotational relaxation time {\tau}2 and diffusion constant D at a temperature Tcr, close to TB, but above the glass transition temperature. Additionally, we coarse-grain the system into cubic boxes, each containing, on average, ~62 particles, to study the average dynamical properties. Clear evidence of large-scale sudden changes in the diffusion coefficient and rotational correlation time signals first-order transitions between low and high-mobility domains.

We show by inherent structure (IS) analysis that the sharp composition dependent breakdown of the Stokes-Einstein relation correlates surprisingly well with an equally sharp non-monotonic variation in the average inherent structure (IS) energy of these mixtures. Further IS analysis reveals the existence of a unique ground state, stabilized by the optimum number of H-bonds at this composition. The surprisingly sharp turnaround behaviour observed in the effective hydrodynamic radius can be traced back to the formation of low energy equilibrium structures at specific compositions., Comment: 16 pages, 8 figures

We study by computer simulations the coupled rotational and translational dynamics of three important linear diatomic molecules, namely, carbon monoxide (CO), nitric oxide (NO) and cyanide ion (CN-) in water. Translational diffusion of these molecules is found to be strongly coupled to their own rotational dynamics which in turn are coupled to similar motions of surrounding water. We examined the validity of hydrodynamic predictions and found them to be largely insufficient, particularly for rotational diffusion. A mode coupling theory approach is developed and applied to understand the complexity of translation-rotation coupling., Comment: 30 pages, 9 figures

We consider a new class of model systems to study systematically the role of molecular shape in the transport properties of dense liquids. Our model is a liquid binary mixture where both the molecules are non-spherical and characterized by a collection of parameters. Although in the real world most of the molecules are non-spherical, only a limited number of theoretical studies exist on the effects of molecular shapes and hardly any have addressed the validity of the hydrodynamic predictions of rotational and translational diffusion of these shapes in liquids. In this work, we study a model liquid consisting of a mixture of prolate and oblate (80:20 mixture) ellipsoidals with interactions governed by a modified Gay-Berne potential for a particular aspect ratio (ratio of the length and diameter of the ellipsoids), at various temperature and pressure conditions. We report calculations of transport properties of this binary mixture by varying temperature over a wide range, at a fixed pressure. We find that for the pressure-density conditions studied, there is no signature of any phase separation, except transitions to the crystalline phase at low temperatures and relatively low pressure (the reason we largely confined our studies to high pressure). We find that for our model binary mixture, both the stick and slip hydrodynamic predictions break down in a major fashion, for both prolates and oblates and particularly so for rotation. Moreover, prolates and oblates themselves display different dynamical features in the mean square displacement and in orientational time correlation functions., Comment: 34 pages, 8 figures

It is commonly believed that melting occurs when mean square displacement (MSD) of a particle of crystalline solid exceeds a threshold value. This is known as the Lindemann criterion, first introduced in the year of 1910 by Lindemann. However, Chakravarty et al. demonstrated that this common wisdom is inadequate because the MSD at melting can be temperature dependent when pressure is also allowed to vary along the coexistence line of the phase diagram [Chakravarty C, Debenedetti P G and Stillinger F H 2007 J. Chem. Phys.126 204508]. We show here by extensive molecular dynamics simulation of both two and three dimensional polydisperse Lennard-Jones solids that particles on the small and large limits of size distribution exhibit substantially different Lindemann ratio at melting. Despite all the dispersion in MSD, melting is found to be first order in both the dimensions at 5-10% dispersity in size., Comment: 22 pages, 6 figures

In order to find a measure of the dynamical features of phase separation kinetics during spinodal decomposition of a liquid binary mixture (like water and cyclohexane , water and 2,6 lutidiene or methanol and cyclohexane), we study both the initial fast exponential-like growth (the Cahn-Hilliard regime) and the subsequent cross-over to a much slower, non-exponential long time growth (the so-called scaling regime), by atomistic molecular dynamics (MD) simulation of a structure breaking binary liquid mixture. In particular, we combine our MD simulations with a coarse grained multi scale modelling (CGMSM) capable of capturing both length and time scales of phase separation kinetics within simulation box. The system is quenched from a higher temperature to two lower temperatures well below the coexistence temperature of the phase diagram. We observe a multiscale phase separation dynamics. Initially the growth is exponential up to a regime of 80-100ps having strong dependence over quench depth. Subsequently a cross-over regime appears where the dynamics slows down considerably due to coarsening through a power law phase. For deeper quench power law growth dominates over the initial exponential and the cross-over becomes transient. We find that for the present parameter values for the binary mixture, the initial rapid growth of structure formation is practically over within 200 ps which is followed by slow structural coarsening. When scaled by the respective viscosities, this time translates to 50-200 ns for water-lutidine binary mixture. The last part of dynamics may extend into ms. The dynamics of phase separation is slowest in regions that have equitable distribution of the two species, and one can observe the signatures of "up-hill diffusion" that is a trade mark of spinodal decomposition., Comment: 35 pages, 7 figures (Submitted to PRE)

We use polydispersity in size as a control parameter to explore certain aspects of melting and freezing transitions in a system of Lennard-Jones spheres. Both analytical theory and computer simulations are employed to establish a potentially interesting relationship between observed terminal polydispersity in Lennard-Jones polydisperse spheres and prediction of the same in the integral equation based theoretical analysis of liquid-solid transition. As we increase polydispersity, solid becomes inherently unstable because of the strain built up due to the size disparity. This aspect is studied here by calculating the inherent structure (IS) calculation. With polydispersity at constant volume fraction we find initially a sharp rise of the average IS energy of the crystalline solid until transition polydispersity, followed by a cross over to a weaker dependence of IS energy on polydispersity in the amorphous state. This cross over from FCC to amorphous state predicted by IS analysis agrees remarkably well with the solid-liquid phase diagram (with extension into the metastable phase) generated by non-linear integral equation theories of freezing. Two other well-known criteria of freezing/melting transitions, the Hansen-Verlet rule of freezing and the Lindemann criterion of melting are both shown to be remarkably in good agreement with the above two estimates. Together they seem to indicate a small range of metastability in the liquid-solid transition in polydisperse solids., Comment: 8 pages, 6 figures

We find through computer simulations that the fractional volume change on freezing of polydisperse Lennard-Jones liquid decreases with increasing polydispersity and approaches zero near a terminal polydispersity of 0.11, independent of temperature. The transition however remains first order at terminal polydispersity. Average inherent structure (IS) energy of the crystalline phase increases nearly quadratically with polydispersity indices ({\delta}) and marked by a crossover to nearly constant IS energy in the amorphous phase., Comment: 13 pages, 3 figures

Binary mixtures are known to phase separate via both nucleation and spinodal decomposition, depending on the initial composition and extent of the quench. Here we develop an energy landscape view of phase separation and non-ideality in binary mixtures by exploring their potential energy landscape (PEL) as functions of temperature and composition, by employing a molecular model that promotes structure breaking abilities of the solute-solvent parent binary liquid at low temperatures. PEL that provides the inherent structure (IS) of a system is obtained by removing the kinetic energy (including that of intermolecular vibrations). Broad Distribution of inherent structure energy demonstrates the larger role of entropy in stabilizing the parent liquid of the structure breaking type of binary mixtures. At high temperature, although the parent structure is homogenous, the corresponding inherent structure is found to be always phase separated, with a density pattern that exhibits marked correlation with the energy of inherent structure. Over broad range of intermediate inherent structure energy, bicontinuous phase separation prevails with interpenetrating stripes as signatures of spinodal decomposition. At low inherent structure energy, the structure is largely phase separated with one interface where as at high inherent structure energy, we find nucleation type growth. Interestingly, at low temperature, the average inherent structure energy () exhibits a drop with temperature which signals onset of crystallization in one of the phases while the other remains in the liquid state. The non-ideal composition dependence of viscosity is anti-correlated with average inherent structure energy., Comment: 35 pages, 10 figures (Phys Rev. E submitted)