We propose an atomistic model for correlated particle dynamics in liquids and glasses predicting both slow stretched-exponential relaxation (SER) and fast compressed-exponential relaxation (CER). The model is based on the key concept of elastically i
nteracting local relaxation events. SER is related to slowing down of dynamics of local relaxation events as a result of this interaction, whereas CER is related to the avalanche-like dynamics in the low-temperature glass state. The model predicts temperature dependence of SER and CER seen experimentally and recovers the simple, Debye, exponential decay at high temperature. Finally, we reproduce SER to CER crossover across the glass transition recently observed in metallic glasses.
Determining the microstructure of colloidal suspensions under shear flows has been a challenge for theoretical and computational methods due to the singularly-perturbed boundary-layer nature of the problem. Previous approaches have been limited to th
e case of hard-sphere systems and suffer from various limitations in their applicability. We present a new analytic scheme based on intermediate asymptotics which solves the Smoluchowski diffusion-convection equation including both intermolecular and hydrodynamic interactions. The method is able to recover previous results for the hard-sphere fluid and, for the first time, to predict the radial distribution function (rdf) of attractive fluids such as the Lennard-Jones (LJ) fluid. In particular, a new depletion effect is predicted in the rdf of the LJ fluid under shear. This method can be used for the theoretical modelling and understanding of real fluids subjected to flow, with applications ranging from chemical systems to colloids, rheology, plasmas, and atmospherical science.
It is widely accepted that structural glasses and disordered crystals exhibit anomalies in the their thermal, mechanical and acoustic properties as manifestations of the breakdown of the long-wavelength approximation in a disordered dissipative envir
onment. However, the same type of glassy-like anomalies (i.e. boson peak in the vibrational density of states (VDOS) above the Debye level, peak in the normalized specific heat at $Tsimeq10 K$ etc) have been recently observed also in perfectly ordered crystals, including thermoelectric compounds. Here we present a theory that predicts these surprising effects in perfectly ordered crystals as a result of low-lying (soft) optical phonons. In particular, it is seen that a strong boson peak anomaly (low-energy excess of modes) in the VDOS can be due almost entirely to the presence of low-energy optical phonons, provided that their energy is comparable to that of the acoustic modes at the Brillouin zone boundary. The boson peak is predicted also to occur in the heat capacity at low $T$. In presence of strong damping (which might be due to anharmonicities in the ordered crystals), these optical phonons contribute to the low-$T$ deviation from Debyes $T^{3}$ law, producing a linear-in-$T$ behavior which is typical of glasses, even though no assumptions of disorder whatsoever are made in the model. These findings are relevant for understanding and tuning thermal transport properties of thermoelectric compounds, and possibly for the enhancement of electron-phonon superconductivity.
We compute the dielectric response of glasses starting from a microscopic system-bath Hamiltonian of the Zwanzig-Caldeira-Leggett type and using an ansatz from kinetic theory for the memory function in the resulting Generalized Langevin Equation. The
resulting framework requires the knowledge of the vibrational density of states (DOS) as input, that we take from numerical evaluation of a marginally-stable harmonic disordered lattice, featuring a strong boson peak (excess of soft modes over Debye $simomega_{p}^{2}$ law). The dielectric function calculated based on this ansatz is compared with experimental data for the paradigmatic case of glycerol at $Tlesssim T_{g}$. Good agreement is found for both the reactive (real part) of the response and for the $alpha$-relaxation peak in the imaginary part, with a significant improvement over earlier theoretical approaches, especially in the reactive modulus. On the low-frequency side of the $alpha$-peak, the fitting supports the presence of $sim omega_{p}^{4}$ modes at vanishing eigenfrequency as recently shown in [Phys. Rev. Lett. 117, 035501 (2016)]. $alpha$-wing asymmetry and stretched-exponential behaviour are recovered by our framework, which shows that these features are, to a large extent, caused by the soft boson-peak modes in the DOS.
Amorphous solids or glasses are known to exhibit stretched-exponential decay over broad time intervals in several of their macroscopic observables: intermediate scattering function, dielectric relaxation modulus, time-elastic modulus etc. This behavi
our is prominent especially near the glass transition. In this Letter we show, on the example of dielectric relaxation, that stretched-exponential relaxation is intimately related to the peculiar lattice dynamics of glasses. By reformulating the Lorentz model of dielectric matter in a more general form, we express the dielectric response as a function of the vibrational density of states (DOS) for a random assembly of spherical particles interacting harmonically with their nearest-neighbours. Surprisingly we find that near the glass transition for this system (which coincides with the Maxwell rigidity transition), the dielectric relaxation is perfectly consistent with stretched-exponential behaviour with Kohlrausch exponents $0.56 < beta < 0.65$, which is the range where exponents are measured in most experimental systems. Crucially, the root cause of stretched-exponential relaxation can be traced back to soft modes (boson-peak) in the DOS.
While a significant body of investigations have been focused on the process of protein self-assembly, much less is understood about the reverse process of a filament breaking due to thermal motion into smaller fragments, or depolymerization of subuni
ts from the filament ends. Indirect evidence for actin and amyloid filament fragmentation has been reported, although the phenomenon has never been directly observed either experimentally or in simulations. Here we report the direct observation of filament depolymerization and breakup in a minimal, calibrated model of coarse-grained molecular simulation. We quantify the orders of magnitude by which the depolymerization rate from the filament ends $k_mathrm{off}$ is larger than fragmentation rate $k_{-}$ and establish the law $k_mathrm{off}/k_- = exp [( varepsilon_| - varepsilon_bot) / k_mathrm{B}T ] = exp [0.5 varepsilon / k_mathrm{B}T ]$, which accounts for the topology and energy of bonds holding the filament together. This mechanism and the order-of-magnitude predictions are well supported by direct experimental measurements of depolymerization of insulin amyloid filaments.
Application of shear flow to charge-stabilized aqueous colloidal suspensions is ubiquitous in industrial applications and as a means to achieve controlled field-induced assembly of nanoparticles. Yet, applying shear flow to a charge-stabilized colloi
dal suspension, which is initially monodisperse and in quasi-equilibrium leads to non-trivial clustering phenomena (and sometimes to a gelation transition), dominated by the complex interplay between DLVO interactions and shear flow. The quantitative understanding of these strongly nonequilibrium phenomena is still far from being complete. By taking advantage of a recent shear-induced aggregation rate theory developed in our group, we present here a systematic numerical study, based on the governing master kinetic equation (population-balance) for the shear-induced clustering and breakup of colloids exposed to shear flow. In the presence of sufficiently stable particles, the clustering kinetics is characterized by an initial very slow growth, controlled by repulsion. During this regime, particles are slowly aggregating to form clusters, the reactivity of which increases along with their size growth. When their size reaches a critical threshold, a very rapid, explosive-like growth follows, where shear forces are able to overcome the energy barrier between particles. This stage terminates when a dynamic balance between shear-induced aggregation and cluster breakage is reached. It is also observed that these systems are characterized by a cluster mass distribution that for a long time presents a well-defined bimodality. The model predictions are quantitatively in excellent agreement with available experimental data, showing how the theoretical picture is able to quantitatively account for the underlying nonequilibrum physics.
The atomic theory of elasticity of amorphous solids, based on the nonaffine response formalism, is extended into the nonlinear stress-strain regime by coupling with the underlying irreversible many-body dynamics. The latter is implemented in compact
analytical form using a qualitative method for the many-body Smoluchowski equation. The resulting nonlinear stress-strain (constitutive) relation is very simple, with few fitting parameters, yet contains all the microscopic physics. The theory is successfully tested against experimental data on metallic glasses, and it is able to reproduce the ubiquitous feature of stress-strain overshoot upon varying temperature and shear rate. A clear atomic-level interpretation is provided for the stress overshoot, in terms of the competition between the elastic instability caused by nonaffine deformation of the glassy cage and the stress buildup due to viscous dissipation.
