Dynamical and orientational structural crossovers in low-temperature glycerol

Orientational dynamics in supercooled glycerol computed from MD simulations: self and cross contributions

The orientational dynamics of supercooled glycerol is probed using molecular dynamics simulations for temperatures ranging from 323 K to 253 K, through correlation functions of first and second ranks of Legendre polynomials, pertaining respectively to dielectric spectroscopy (DS) and depolarized dynamic light scattering (DDLS). The self, cross, and total correlation functions are compared with relevant experimental data. The computations reveal the low sensitivity of DDLS to cross-correlations, in agreement with what is found in experimental work, and strengthen the idea of directly comparing DS and DDLS data to evaluate the effect of cross-correlations in polar liquids. The analysis of the net static cross-correlations and their spatial decomposition shows that, although cross-correlations extend over nanometric distances, their net magnitude originates, in the case of glycerol, from the first shell of neighbouring molecules. Accessing the angular dependence of the static correlation allows us to get a microscopic understanding of why the rank-1 correlation function is more sensitive to cross-correlation than its rank-2 counterpart.

War and Peace between Electrostatic and van der Walls Forces Regulate Translational and Rotational Diffusion

Diffusion of a Brownian particle or a molecule in a liquid solvent is caused, in the Stokes-Einstein picture, by unbalanced fluctuations of osmotic forces on different sides of the particle. When the particle carries a charge or a higher multipolar moment, this picture is amended by fluctuations of electrostatic forces producing dielectric friction. Dielectric friction slows both the translational and rotational diffusion down. While this picture is well established and is physically sound, standard theories grossly overestimate the magnitude of dielectric friction for small dipolar solutes and larger colloidal particles, such as proteins. Motivated by recent simulation studies, this perspective discusses the interplay between osmotic (van der Waals) and electrostatic forces in promoting molecular and colloidal diffusion. Much can be learned about microscopic friction mechanisms from statistical and dynamical correlations between osmotic and electrostatic forces.

First-principles study of the specific heat of glass at the glass transition with a case study on glycerol

The standard method to determine the transition temperature (T_g) of glasses is the jump in the specific heat,ΔCp. Despite its importance, standard theory for this jump is lacking. The difficulties include lack of proper treatment of the specific heat of liquids, hysteresis, and the timescale issue. The first part of this paper provides a non-empirical method for calculating the specific heat in the glass transition. The method consists of molecular dynamics (MD) simulations based on density-functional theory (DFT) and thermodynamics methods. Calculation of the total energy, which is the heart of DFT, is the most general method for obtaining specific heat for any state of matters. The influence of energy dissipation processes on specific heat is treated by adiabatic MD simulations. The problems of hysteresis and the timescale are alleviated by restricting the scope of calculations to equilibrium states only. The second part of this paper demonstrates the validity and usefulness of the methods by applying to the specific-heat jump of glycerol. By decomposingΔCpinto contributions of the structural, phonon, and thermal expansion energies, an appropriate interpretation for the specific-heat jump has been established: the major contribution toΔCpis the change in the structural energy. From this, a neat energy diagram about the glass transition is obtained. An outcome of this study is verification of the empirical relationship between the fragility and the specific-heat jump. These two quantities scale to the ratiok=Tg/ΔTg, whereΔTgis the width of the transition, through which the two quantities are interrelated.

Hydrogen Bonding and the Structural Properties of Glycerol-Water Mixtures with a Microwave Field: a Molecular Dynamics Study

In a microwave field, the dielectric properties, molecular structures, and hydrogen bonding dynamics of glycerol in its mixtures with water were determined by the molecular dynamics simulation method. The dipole-dipole correlation of glycerol is linked to the field intensity of microwaves. The results show that as the field intensity is increased, even glycerol in the second coordination shell can become correlated with each other. The structures of up to 35 glycerol molecules are observed. More than that, it was observed that lifetimes of glycerol-glycerol hydrogen bonds were prolonged, while the average hydrogen bond number was also increased. Besides, the structures in a strong microwave field mimic the weak C-H⋯O hydrogen bonds seen in high-glycerol concentration mixtures, yet the concentration is lower. These results indicate that with the assistance of the microwave field, glycerol molecules become concentrated and are more likely to establish stable interactions with others. As a consequence, the spherical clusters composed by glycerol molecules in our nanosheet synthesis experiment are easier to form.

Complex Dynamics of Water in Protein Confinement

This paper studies single-molecule and collective dynamics of water confined in protein powders by means of molecular dynamics simulations. The single-particle dynamics show a modest retardation compared to the bulk but become highly stretched in the powder, with the stretching exponent of ≃0.2. The collective dynamics of the total water dipole are affected by intermolecular correlations inside water and by cross-correlations between the water and the protein. The dielectric spectrum of water in the powder has two nearly equal-amplitude peaks: a Debye peak with ≃16 ps relaxation time and a highly stretched peak with the relaxation time of ≃13 ns and a stretching exponent of ≃0.12. The slower relaxation component is not seen in the single-molecule correlation functions and can be assigned to elastic protein motions displacing water in the powder. The loss spectrum of the intermediate scattering function reported by neutron-scattering experiments is also highly stretched, with the high-frequency wing scaling according to a power law. Translational dynamics can become much slower in the powder than in the bulk but are overshadowed by the rotational loss in the overall loss spectrum of neutron scattering.