No Arabic abstract
Freezing is a fundamental physical phenomenon that has been studied over many decades; yet the role played by surfaces in determining nucleation has remained elusive. Here we report direct computational evidence of surface induced nucleation in supercooled systems with a negative slope of their melting line (dP/dT < 0). This unexpected result is related to the density decrease occurring upon crystallization, and to surface tension facilitating the initial nucleus formation. Our findings support the hypothesis of surface induced crystallization of ice in the atmosphere, and provide insight, at the atomistic level, into nucleation mechanisms of widely used semiconductors.
A theoretical treatment of deeply supercooled liquids is difficult because their properties emerge from spatial inhomogeneities that are self-induced, transient, and nanoscopic. I use computer simulations to analyse self-induced static and dynamic heterogeneity in equilibrium systems approaching the experimental glass transition. I characterise the broad sample-to-sample fluctuations of salient dynamic and thermodynamic properties in elementary mesoscopic systems. Findings regarding local lifetimes and distributions of dynamic heterogeneity are in excellent agreement with recent single molecule studies. Surprisingly broad thermodynamic fluctuations are also found, which correlate well with dynamics fluctuations, thus providing a local test of the thermodynamic origin of slow dynamics.
By means of Raman spectroscopy of liquid microjets we have investigated the crystallization process of supercooled quantum liquid mixtures composed of parahydrogen (pH$_2$) diluted with small amounts of up to 5% of either neon or orthodeuterium (oD$_2$), and of oD$_2$ diluted with either Ne or pH$_2$. We show that the introduction of Ne impurities affects the crystallization kinetics in both the pH$_2$-Ne and oD$_2$-Ne mixtures in terms of a significant reduction of the crystal growth rate, similarly to what found in our previous work on supercooled pH$_2$-oD$_2$ liquid mixtures [M. Kuhnel et {it al.}, Phys. Rev. B textbf{89}, 180506(R) (2014)]. Our experimental results, in combination with path-integral simulations of the supercooled liquid mixtures, suggest in particular a correlation between the measured growth rates and the ratio of the effective particle sizes originating from quantum delocalization effects. We further show that the crystalline structure of the mixture is also affected to a large extent by the presence of the Ne impurities, which likely initiate the freezing process through the formation of Ne crystallites.
Glasses are solid materials whose constituent atoms are arranged in a disordered manner. The transition from a liquid to a glass remains one of the most poorly understood phenomena in condensed matter physics, and still no fully microscopic theory exists that can describe the dynamics of supercooled liquids in a quantitative manner over all relevant time scales. Here we present such a theoretical framework that yields near-quantitative accuracy for the time-dependent correlation functions of a supercooled system over a broad density range. Our approach requires only simple static structural information as input and is based entirely based on first principles. Owing to this first-principles nature, the framework offers a unique platform to study the relation between structure and dynamics in glass-forming matter, and paves the way towards a systematically correctable and ultimately fully quantitative theory of microscopic glassy dynamics.
We investigate the relaxation mechanism of a supercooled tetrahedral liquid at its limit of stability using isothermal isobaric ($NPT$) Monte Carlo (MC) simulations. In similarity with systems which are far from equilibrium but near the onset of jamming [OHern et.al., Phys. Rev. Lett. {bf 93}, 165702 (2004)], we find that the relaxation is characterized by two time-scales: the decay of long-wavelength (slow) fluctuations of potential energy is controlled by the the slope $[partial (G/N)/partial phi]$ of the Gibbs free energy ($G$) at a unique value of per particle potential energy $phi = phi_{mid}$. The short-wavelength (fast) fluctuations are controlled by the bath temperature $T$. The relaxation of the supercooled liquid is initiated with a dynamical crossover after which the potential energy fluctuations are biased towards values progressively lesser than $phi_{mid}$. The dynamical crossover leads to the change of time-scale, i.e., the decay of long-wavelength potential energy fluctuations (intermediate relaxation). Because of the condition [$partial^2 (G/N)/partial phi^2 = 0$] at $phi = phi_{mid}$, the slope $[partial (G/N)/partial phi]$ has a unique value and governs the intermediate stage of relaxation, which ends just after the crossover. In the subsequent stage, there is a relatively rapid crystallization due to lack of long-wavelength fluctuations and the instability at $phi_{mid}$, i.e., the condition that $G$ decreases as configurations with potential energies lower than $phi_{mid}$ are accessed. The dynamical crossover point and the associated change in the time-scale of fluctuations is found to be consistent with the previous studies.
A novel liquid-liquid phase transition has been proposed and investigated in a wide variety of pure substances recently, including water, silica and silicon. From computer simulations using the Stillinger-Weber classical empirical potential, Sastry and Angell [1] demonstrated a first order liquid-liquid transition in supercooled silicon, subsequently supported by experimental and simulation studies. Here, we report evidence for a liquid-liquid critical end point at negative pressures, from computer simulations using the SW potential. Compressibilities exhibit a growing maximum upon lowering temperature below 1500 K and isotherms exhibit density discontinuities below 1120 K, at negative pressure. Below 1120 K, isotherms obtained from constant volume-temperature simulations exhibit non-monotonic, van der Waals-like behavior signaling a first order transition. We identify Tc ~ 1120 +/- 12 K, Pc -0.60 +/- 0.15 GPa as the critical temperature and pressure for the liquid-liquid critical point. The structure of the liquid changes dramatically upon decreasing the temperature and pressure. Diffusivities vary over 4 orders of magnitude, and exhibit anomalous pressure dependence near the critical point. A strong relationship between local geometry quantified by the coordination number, and diffusivity, is seen, suggesting that atomic mobility in both low and high density liquids can usefully be analyzed in terms of defects in the tetrahedral network structure. We have constructed the phase diagram of supercooled silicon. We identify the lines of compressibility, density extrema (maxima and minima) and the spinodal which reveal the interconnection between thermodynamic anomalies and the phase behaviour of the system as suggested in previous works [2-9]