No Arabic abstract
We analyze multiple new issues concerning activated relaxation in glassy hard sphere fluids and molecular and polymer liquids based on the Elastically Collective Nonlinear Langevin Equation (ECNLE) theory. By invoking a high temperature reference state, a near universality of the apparent dynamic localization length scale is predicted for liquids of widely varying fragility, a result that is relevant to recent simulation studies and quasi-elastic neutron scattering measurements. In contrast, in the same format strongly non-universal behavior is found for the activation barrier that controls long time relaxation. Two measures of cooperativity in ECNLE theory are analyzed. A particle-level total displacement associated with the alpha relaxation event is found to be only of order 1-2 particle diameters and weakly increases with cooling. In contrast, an alternative cooperativity length is defined as the spatial scale required to recover the full barrier and bulk alpha time. This length scale grows strongly with cooling due to the emergence in the deeply supercooled regime of collective long range elastic fluctuations required to allow local hopping. It becomes very large as the laboratory Tg is approached, though is relatively modest at degrees of supercooling accessible with molecular dynamics simulation. The alpha time is found to be exponentially related to this cooperativity length over an enormous number of decades of relaxation time that span the lightly to deeply supercooled regimes. Moreover, the effective barrier height increases almost linearly with the growing cooperativity length scale. An alternative calculation of the collective elastic barrier based on a literal continuum mechanics approach is shown to result in very little change of the theoretical results for bulk properties, but leads to a much smaller and less temperature-sensitive cooperativity length scale.
Theoretical approaches are formulated to investigate the molecular mobility under various cooling rates of amorphous drugs. We describe the structural relaxation of a tagged molecule as a coupled process of cage-scale dynamics and collective molecular rearrangement beyond the first coordination shell. The coupling between local and non-local dynamics behaves distinctly in different substances. Theoretical calculations for the structural relaxation time, glass transition temperature, and dynamic fragility are carried out over twenty-two amorphous drugs and polymers. Numerical results have a quantitatively good accordance with experimental data and the extracted physical quantities using the Vogel-Fulcher-Tammann fit function and machine learning. The machine learning method reveals the linear relation between the glass transition temperature and the melting point, which is a key factor for pharmaceutical solubility. Our predictive approaches are reliable tools for developing drug formulation.
The sluggish and heterogeneous dynamics of glass forming liquids is frequently associated to the transient coexistence of two phases of particles, respectively with an high and low mobility. In the absence of a dynamical order parameter that acquires a transient bimodal shape, these phases are commonly identified empirically, which makes difficult investigating their relation with the structural properties of the system. Here we show that the distribution of single particle diffusivities can be accessed within a Continuous Time Random Walk description of the intermittent motion, and that this distribution acquires a transient bimodal shape in the deeply supercooled regime, thus allowing for a clear identification of the two coexisting phase. In a simple two-dimensional glass forming model, the dynamic phase coexistence is accompanied by a striking structural counterpart: the distribution of the crystalline-like order parameter becomes also bimodal on cooling, with increasing overlap between ordered and immobile particles. This simple structural signature is absent in other models, such as the three-dimesional Kob-Andersen Lennard-Jones mixture, where more sophisticated order parameters might be relevant. In this perspective, the identification of the two dynamical coexisting phases opens the way to deeper investigations of structure-dynamics correlations.
We numerically study the relaxation dynamics of several glass-forming models to their inherent structures, following quenches from equilibrium configurations sampled across a wide range of temperatures. In a mean-field Mari-Kurchan model, we find that relaxation changes from a power-law to an exponential decay below a well-defined temperature, consistent with recent findings in mean-field $p$-spin models. By contrast, for finite-dimensional systems, the relaxation is always algebraic, with a non-trivial universal exponent at high temperatures crossing over to a harmonic value at low temperatures. We demonstrate that this apparent evolution is controlled by a temperature-dependent population of localised excitations. Our work unifies several recent lines of studies aiming at a detailed characterization of the complex potential energy landscape of glass-formers.
It was recently shown that the real part of the frequency-dependent fluidity for several glass-forming liquids of different chemistry conforms to the prediction of the random barrier model (RBM) devised for ac electrical conduction in disordered solids [S. P. Bierwirth textit{et al.}, Phys. Rev. Lett. {bf 119}, 248001 (2017)]. Inspired by these results we introduce a crystallization-resistant modification of the Kob-Andersen binary Lennard-Jones mixture for which the results of extensive graphics-processing unit (GPU)-based molecular-dynamics simulations are presented. We find that the low-temperature mean-square displacement is fitted well by the RBM prediction, which involves no shape parameters. This finding highlights the challenge of explaining why a simple model based on hopping of non-interacting particles in a fixed random energy landscape can reproduce the complex and highly cooperative dynamics of glass-forming liquids.
We theoretically investigate structural relaxation and activated diffusion of glass-forming liquids at different pressures using both the Elastically Collective Nonlinear Langevin Equation (ECNLE) theory and molecular dynamics (MD) simulation. An external pressure restricts local motions of a single molecule within its cage and triggers the slowing down of cooperative mobility. While the ECNLE theory and simulation generally predict a monotonic increase of the glass transition temperature and dynamic fragility with pressure, the simulation indicates a decrease of fragility as pressure above 1000 bar. The structural relaxation time is found to be linearly coupled with the inverse diffusion constant. Remarkably, this coupling is independent of compression. Theoretical calculations agree quantitatively well with simulations and are also consistent with prior works.