Search

We study melting in a two-dimensional system of classical particles with Gaussian-core interactions in disordered environments. The pure system validates the conventional two-step melting with a hexatic phase intervening between the solid and the liquid. This picture is modified in the presence of pinning impurities. A random distribution of pinning centers forces a hexatic-like low temperature phase that transits into a liquid at a single melting temperature $T^{\rm RP}_{\rm m}$. In contrast, pinning centers located at randomly chosen sites of a perfect crystal anchors a solid at low temperatures which undergoes a direct transition to the liquid at $T^{\rm CP}_{\rm m}$. Thus, the two-step melting is lost in either cases of disorder. We discuss the characteristics of melting depending on the nature of the impurities., Comment: 11 pages, 13 figures

Recently, some of us developed a novel model glass-forming liquid with k extra interactions with pseudo neighbours to each liquid particle over and above the regular interactions with its neighbours. Analysis of the structure and dynamics of these systems showed that with an increase in k the systems have more mean-field like properties. This work presents an extensive study of the thermodynamics of the above-mentioned model for several values of k and its correlation with the dynamics. We surprisingly find that the usual thermodynamic integration (TI) method of calculating the entropy provides unphysical results for this model. It predicts the vanishing of configurational entropy at state points at which both the collective and the single-particle dynamics of the system show complete relaxation. We then employ a new method known as the two-phase thermodynamics (2PT) method to calculate the entropy. We find that with an increase in k the difference in the entropy computed using the two methods (2PT and TI) increases. We also find that in the temperature range studied, the entropy calculated via the 2PT method satisfies the Adam-Gibbs (AG) relationship between the relaxation time and the configurational entropy, whereas the entropy calculated via the TI method shows a strong violation of the same. We then apply the 2PT method to calculate the entropy in another system where some fractions of particles are pinned randomly in their equilibrium positions. This system also shows a similar breakdown of the AG relationship as reported earlier. We show that the difference in entropy calculated via the 2PT and TI methods increases with an increase in pinning density. We also find that when the entropy is calculated using the 2PT method, the AG relationship between the dynamics and the entropy holds., Comment: 16 pages, 27 figures, 3 tables

We study various properties of the vibrational normal modes for Coulomb-interacting particles in two-dimensional irregular confinement using numerical simulations. By analyzing the participation ratio and spectral statistics, we characterize the vibrational modes for Coulomb clusters as localized, quasi-localized and delocalized. We also study a novel correlation function to understand the spatial structure of these different kinds of modes and subsequently extract the associated characteristic length scales. We further demonstrate that, at any given temperature, particles exhibiting larger displacement over a time interval comparable to the structural relaxation time, are strongly correlated with the low-frequency quasi-localized modes of the inherent structure corresponding to the initial configuration. Establishing this correlation for Coulomb clusters paves the path to identify the particular feature of the initial configuration that determines the previously observed heterogeneous dynamics of the particles at low temperatures in these systems., Comment: 13 pages, 10 figures

The origins of rapid dynamical slow down in glass forming liquids in the growth of static length scales, possibly associated with identifiable structural ordering, is a much debated issue. Growth of medium range crystalline order (MRCO) has been observed in various model systems to be associated with glassy behaviour. Such observations raise the question about the eventual state reached by a glass former, if allowed to relax for sufficiently long times. Is a slowly growing crystalline order responsible for slow dynamics? Are the molecular mechanisms for glass transition in liquids with and without MRCO the same? If yes, glass formers with MRCO provide a paradigm for understanding glassy behaviour generically. If not, systems with MRCO form a new class of glass forming materials whose molecular mechanism for slow dynamics may be easier to understand in terms of growing crystalline order, and should be approached in that manner, even while they will not provide generic insights. In this study we perform extensive molecular dynamics simulations of a number of glass forming liquids in two dimensions and show that the static and dynamic properties of glasses with MRCO are different from other glass forming liquids with no predominant local order. We also resolve an important issue regarding the so-called Point-to-set method for determining static length scales, and demonstrate it to be a robust, order agnostic, method for determining static correlation lengths in glass formers.

Thermal conductivity of a model glass-forming system in the liquid and glass states is studied using extensive numerical simulations. We show that near the glass transition temperture, where the structural relaxation time becomes very long, the measured thermal conductivity decreases with increasing age. Secondly the thermal conductivity of the disordered solid obtained at low temperatures depends on the cooling rate with which it was prepared, with lower cooling rates leading to lower thermal conductivity. Our analysis links this decrease of the thermal conductivity with increased exploration of lower-energy inherent structures of the underlying potential energy landscape. Further we show that the lowering of conductivity for lower-energy inherent structures is related to the high frequency harmonic modes associated with the inherent structure being less extended., Comment: 10 pages, 14 figures

The Adam-Gibbs relation between relaxation times and the configurational entropy has been tested extensively for glass formers using experimental data and computer simulation results. Although the form of the relation contains no dependence on the spatial dimensionality in the original formulation, subsequent derivations of the Adam-Gibbs relation allow for such a possibility. We test the Adam-Gibbs relation in 2, 3, and 4 spatial dimensions using computer simulations of model glass formers. We find that the relation is valid in 3 and 4 dimensions. But in 2 dimensions, the relation does not hold, and interestingly, no single alternate relation describes the results for the different model systems we study., Comment: Submitted to Phys. Rev. Lett

We develop a framework for understanding the difference between strong and fragile behavior in the dynamics of glass-forming liquids from the properties of the potential energy landscape. Our approach is based on a master equation description of the activated jump dynamics among the local minima of the potential energy (the so-called inherent structures) that characterize the potential energy landscape of the system. We study the dynamics of a small atomic cluster using this description as well as molecular dynamics simulations and demonstrate the usefulness of our approach for this system. Many of the remarkable features of the complex dynamics of glassy systems emerge from the activated dynamics in the potential energy landscape of the atomic cluster. The dynamics of the system exhibits typical characteristics of a strong supercooled liquid when the system is allowed to explore the full configuration space. This behavior arises because the dynamics is dominated by a few lowest-lying minima of the potential energy and the potential energy barriers between these minima. When the system is constrained to explore only a limited region of the potential energy landscape that excludes the basins of attraction of a few lowest-lying minima, the dynamics is found to exhibit the characteristics of a fragile liquid., Comment: 13 pages, 6 figures

We study the growing time scales and length scales associated with dynamical slow down for a realistic glass former, using computer simulations. We perform finite size scaling to evaluate a length scale associated with dynamical heterogeneity which grows as temperature decreases. However, relaxation times which also grow with decreasing temperature, do not show the same kind of scaling behavior with system size as the dynamical heterogeneity, indicating that relaxation times are not solely determined by the length scale of dynamical heterogeneity. We show that relaxation times are instead determined, for all studied system sizes and temperatures, by configurational entropy, in accordance with the Adam-Gibbs relation, but in disagreement with theoretical expectations based on spin-glass models that configurational entropy is not relevant at temperatures substantially above the critical temperature of mode coupling theory. The temperature dependence of the heterogeneity length scale shows significant deviations from theoretical expectations, and the length scale one may extract from the system size dependence of the configurational entropy has much weaker temperature dependence compared to the heterogeneity length scale at all studied temperatures. Our results provide new insights into the dynamics of glass-forming liquids and pose serious challenges to existing theoretical descriptions.

We have used a mean-field Monte Carlo method to study the zero-temperature synchronous dynamics of a one-pattern model of associative memory with random asymmetric couplings. In the case of symmetric couplings, we find evidence for a transition from a spin-glass-like phase to a ferromagnet-like phase as the acquisition strength of the stored pattern is increased from zero. In the ferromagnetic phase, we find the existence of two types of phase-space structure for $m > 0$ where $m$ is the overlap of the state of the system with the stored pattern: a simple phase-space structure where all initial states with $m > 0$ flow to the attractor corresponding to the stored pattern; and a complex phase-space structure with many attractors with their basins of attraction. The presence of random asymmetry in the couplings results in better retrieval performance of the network by enhancing the size of the basin of attraction of the stored pattern and by making the recall of memory significantly faster.

We report numerical results for the phase diagram in the density-disorder plane of a hard sphere system in the presence of quenched, random, pinning disorder. Local minima of a discretized version of the Ramakrishnan-Yussouff free energy functional are located numerically and their relative stability is studied as a function of the density and the strength of disorder. Regions in the phase diagram corresponding to liquid, glassy and nearly crystalline states are mapped out, and the nature of the transitions is determined. The liquid to glass transition changes from first to second order as the strength of the disorder is increased. For weak disorder, the system undergoes a first order crystallization transition as the density is increased. Beyond a critical value of the disorder strength, this transition is replaced by a continuous glass transition. Our numerical results are compared with those of analytical work on the same system. Implications of our results for the field-temperature phase diagram of type-II superconductors are discussed., Comment: 14 pages, 10 postscript figures (included), submitted to Phys. Rev. E

Properties of the free energy landscape in phase space of a dense hard sphere system characterized by a discretized free energy functional of the Ramakrishnan-Yussouff form are investigated numerically. A considerable number of glassy local minima of the free energy are located and the distribution of an appropriately defined ``overlap'' between minima is calculated. The process of transition from the basin of attraction of a minimum to that of another one is studied using a new ``microcanonical'' Monte Carlo procedure, leading to a determination of the effective height of free energy barriers that separate different glassy minima. The general appearance of the free energy landscape resembles that of a putting green: deep minima separated by a fairly flat structure. The growth of the effective free-energy barriers with increasing density is consistent with the Vogel-Fulcher law, and this growth is primarily driven by an entropic mechanism., Comment: 10 pages, 6 postscript figures, uses iopart.cls and iopart10.clo (included). Invited talk at the ICTP Trieste Conference on "Unifying Concepts in Glass Physics", September 1999. To be published in J. Phys. Cond. Mat

The topography of the free energy landscape in phase space of a dense hard sphere system characterized by a discretized free energy functional of the Ramakrishnan-Yussouff form is investigated numerically using a ``microcanonical'' Monte Carlo procedure. We locate a considerable number of glassy local minima of the free energy and analyze the distributions of the free energy at a minimum and an appropriately defined phase-space ``distance'' between different minima. We find evidence for the existence of pairs of closely related glassy minima (``two-level systems''). We also investigate the way the system makes transitions as it moves from the basin of attraction of a minimum to that of another one after a start under nonequilibrium conditions. This allows us to determine the effective height of free energy barriers that separate a glassy minimum from the others. The dependence of the height of free energy barriers on the density is investigated in detail. The general appearance of the free energy landscape resembles that of a putting green: relatively deep minima separated by a fairly flat structure. We discuss the connection of our results with the Vogel-Fulcher law and relate our observations to other work on the glass transition., Comment: 14 pages plus 10 figs