Structural studies of liquid D-glycerol by neutron diffraction

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.

Structures of glass-forming liquids by x-ray scattering: Glycerol, xylitol, and D-sorbitol

Synchrotron x-ray scattering has been used to investigate three liquid polyalcohols of different sizes (glycerol, xylitol, and D-sorbitol) from above the glass transition temperatures T_g to below. We focus on two structural orders: the association of the polar OH groups by hydrogen bonds (HBs) and the packing of the non-polar hydrocarbon groups. We find that the two structural orders evolve very differently, reflecting the different natures of bonding. Upon cooling from 400 K, the O⋯O correlation at 2.8 Å increases significantly in all three systems, indicating more HBs, until kinetic arrests at T_g; the increase is well described by an equilibrium between bonded and non-bonded OH with ΔH = 9.1 kJ/mol and ΔS = 13.4 J/mol/K. When heated above T_g, glycerol loses the fewest HBs per OH for a given temperature rise scaled by T_g, followed by xylitol and by D-sorbitol, in the same order the number of OH groups per molecule increases (3, 5, and 6). The pair correlation functions of all three liquids show exponentially damped density modulations of wavelength 4.5 Å, which are associated with the main scattering peak and with the intermolecular C⋯C correlation. In this respect, glycerol is the most ordered with the most persistent density ripples, followed by D-sorbitol and by xylitol. Heating above T_g causes faster damping of the density ripples with the rate of change being the slowest in xylitol, followed by glycerol and by D-sorbitol. Given the different dynamic fragility of the three liquids (glycerol being the strongest and D-sorbitol being the most fragile), we relate our results to the current theories of the structural origin for the difference. We find that the fragility difference is better understood on the basis of the thermal stability of HB clusters than that of the structure associated with the main scattering peak.

Molecular dynamics simulations vs field-cycling NMR relaxometry: Structural relaxation mechanisms in the glass-former glycerol revisited

We combine field-cycling (FC) relaxometry and molecular dynamics (MD) simulations to study the rotational and translational dynamics associated with the glassy slowdown of glycerol. The ¹H NMR spin-lattice relaxation rates R₁(ω) probed in the FC measurements for different isotope-labelled compounds are computed from the MD trajectories for broad frequency and temperature ranges. We find high correspondence between experiment and simulation. Concerning the rotational motion, we observe that the aliphatic and hydroxyl groups show similar correlation times but different stretching parameters, while the overall reorientation associated with the structural relaxation remains largely isotropic. Additional analysis of the simulation results reveals that transitions between different molecular configurations are slow on the time scale of the structural relaxation at least at sufficiently high temperatures, indicating that glycerol rotates at a rigid entity, but the reorientation is slower for elongated than for compact conformers. The translational contribution to R₁(ω) is well described by the force-free hard sphere model. At sufficiently low frequencies, universal square-root laws provide access to the molecular diffusion coefficients. In both experiment and simulation, the time scales of the rotational and translational motions show an unusually large separation, which is at variance with the Stokes-Einstein-Debye relation. To further explore this effect, we investigate the structure and dynamics on various length scales in the simulations. We observe that a prepeak in the static structure factor S(q), which is related to a local segregation of aliphatic and hydroxyl groups, is accompanied by a peak in the correlation times τ(q) from coherent scattering functions.

Pressure Dependence of the Boson Peak of Glassy Glycerol

The pressure dependence of the boson peak (BP) of glycerol, including its behavior across the liquid-glass transition, has been studied using Raman scattering. A significant increase of the BP frequency was observed with pressure up to 11 GPa at room temperature. The pressure dependence of BP frequency ν_BP is proportional to (1+P/P₀)^1/3, where P and P₀ are the pressure and a constant, respectively, consistent with a soft potential model. The characteristic length of medium range order is close in size to a cyclic trimer of glycerol molecules, as predicted by the medium range order of a BP excitation using molecular dynamics simulations, and the pressure dependence of a characteristic medium range order is nearly constant. The pressure induced structural changes in glycerol can be understood in terms of the shrinkage of voids with cyclic trimers persisting to at least 11 GPa. Pressure dependence of the intermolecular O-H stretching mode indicates that the intermolecular hydrogen bond distances gradually decrease up to the glass transition pressure of ∼5 GPa and become nearly constant in the glassy state, indicating the disappearance of free volume in the dense glass.

Experimental study on the relationship between the frequency-dependent shear viscosity and the intermediate scattering function of representative viscous liquids

The frequency-dependent shear viscosity of two representative viscous liquids, o-terphenyl and glycerin, was experimentally determined at several temperatures and compared with the intermediate scattering functions reported in the literature. A comparison based on mode-coupling theory succeeded in relating the frequency-dependent shear viscosity with the intermediate scattering function at the main peak of the static structure factor. It suggests that the slow relaxation mode of the shear viscosity of both liquids is governed by the density fluctuation at the main peak of the static structure factor, in spite of the differences in the details of their intermolecular interactions.

Structural heterogeneities at the origin of acoustic and transport anomalies in glycerol glass-former

By means of large scale molecular dynamics simulations, we explore mesoscopic properties of prototypical glycerol glass-former above and below the glass transition. The model used, in excellent agreement with various experimental techniques, permits to carefully study the structure and the vibrational dynamics. We find that a medium range order is present in glycerol glass-former and arises from hydrogen bond network extension. The characteristic size of the structural heterogeneities is related to the anomalous properties of acoustic vibrations (Rayleigh scattering, "mode softening," and Boson Peak) in the glassy state. Finally the characteristic size of these heterogeneities, nearly constant in temperature, is also connected to the cross-over between structural relaxation and diffusion in liquid glycerol.

A combined experimental and computational study of the esterification reaction of glycerol with acetic acid

This work describes theoretical and experimental studies on glycerol esterification to obtain acetins focusing on the obtained isomers. The reaction of glycerol with acetic acid was carried out on Amberlyst 36 wet. Density functional theory calculations on the level of M06-2X functional and 6-311+G(d,p) basis set are carried out and the most stable structures of the reactants and products are located by considering a large number of conformers. The thermodynamics is discussed in terms of the calculated reaction Gibbs free energy. The AIM theory was used to characterize reactants and products. The glycerol esterification with acetic acid is found to be thermodynamically favored, with exothermal property. These agree well with experiments and allow us to explain the relative selectivity of products.

Salty glycerol versus salty water surface organization: bromide and iodide surface propensities

Salty NaBr and NaI glycerol solution interfaces are examined in the OH stretching region using broadband vibrational sum frequency generation (VSFG) spectroscopy. Raman and infrared (IR) spectroscopy are used to further understand the VSFG spectroscopic signature. The VSFG spectra of salty glycerol solutions reveal that bromide and iodide anions perturb the interfacial glycerol organization in a manner similar as that found in aqueous halide salt solutions, thus confirming the presence of bromide and iodide anions at the glycerol surface. Surface tension measurements are consistent with the surface propensity suggested by the VSFG data and also show that the surface excess increases with increasing salt concentration, similar to that of water. In addition, iodide is shown to have more surface prevalence than bromide, as has also been determined from aqueous solutions. These results suggest that glycerol behaves similarly to water with respect to surface activity and solvation of halide anions at its air/liquid interface.

Relationships between Positron Lifetime and Dynamics of Polymers

In this work positron annihilation lifetime spectroscopy (PALS) is employed to study the ortho-Positronium lifetime parameters τ3>, σ3 and I3 in two structurally simple amorphous polymers 1,2-poly(butadiene) and cis-1,4-poly(isoprene). The ortho-Positronium mean lifetime parameter τ3> is compared with dynamics data from broadband dielectric spectroscopy (BDS) and electron spin resonance (ESR) experiments. Coincidences of characteristic temperatures from PALS with BDS and ESR experiments show close relationships between positron lifetime and dynamics of polymers.

The non-Gaussian dynamics of glycerol

We have combined incoherent quasielastic neutron scattering experiments and atomistic molecular simulations to investigate the microscopic dynamics of glycerol moving away from the hydrodynamic limit. We relate changes in the momentum transfer (Q) dependence of the relaxation time to distinct changes of the single-particle dynamics. Going from small to large values of Q, a first crossover at about 0.5 Å(-1) is related to the coupling of the translational diffusion dynamics to the non-Debye structural relaxation, while the second crossover at a Q-value near the main diffraction peak is associated with the Gaussian to non-Gaussian crossover of the short-time molecular dynamics, related to the decaging processes. We offer an unprecedented extension of previous studies on polymeric systems towards the case of the typical low-molecular-weight glass-forming system glycerol.

Positron annihilation response and broadband dielectric spectroscopy: salol

A phenomenological analysis of the ortho-positronium (o-Ps) annihilation from positron annihilation lifetime spectroscopy (PALS) and the dynamics from broadband dielectric spectroscopy (BDS) are reported on a small molecular glass former of intermediate H-bonding and fragility: salol. The dielectric spectra extend over a very broad frequency range of about 2 × 10(-2)-3.5 × 10(11) Hz, providing information on the α-relaxation, the secondary relaxation giving rise to the excess wing, and the shallow high-frequency minimum in the micro- to milli-meter wave range. A number of empirical correlations between the o-Ps lifetime, τ(3)(T), and the various spectral and relaxation features have been observed. Thus, the phenomenological evaluation of the τ(3)(T) dependence of the PALS response of the amorphous sample reveals three characteristic PALS temperatures: T(g)(PALS), T(b1)(L) = 1.15T(g)(PALS) and T(b2)(L) = 1.25T(g)(PALS), which are discussed in relation to similar findings for some typical small molecular vdW- and H-bonded glass formers. A slighter change of the slope at T(b1)(L) appears to be related to the transition from excess wing to the primary α-process-dominated behavior, with the secondary process dominating in the deeply supercooled liquid state below T(b1)(L). The high-temperature plateau effect in the τ (3)(T) plot occurs at T(b2)(L) and agrees with the characteristic Stickel temperature, T(B)(ST), marking a qualitative change of the primary α process, but it does not follow the relation T(b2)(L) < T(α) [τ(3)(T(b2)) < τ(α)]. Both effects at T(b1)(L) and T(b2)(L) correlate with two crossovers in the spectral shape and related non-exponentiality parameter of the structural relaxation, β (KWW). Finally, the application of the two-order parameter (TOP) model to the structural relaxation as represented by the primary α relaxation times from BDS leads to the characteristic TOP temperature, T(m)(c), close to T(b1) from PALS. Within this model the phenomenological interpretation is offered based on changes in the probability of occurrence of solid-like and liquid-like domains to explain the dynamic as well as PALS responses. In summary, all the empirical correlations support further very close connections between the PALS response and the dielectric relaxation behavior in small molecule glass formers.

On achieving experimental accuracy from molecular dynamics simulations of flexible molecules: aqueous glycerol

The rotational isomeric states (RIS) of glycerol at infinite dilution have been characterized in the aqueous phase via a 1 micros conventional molecular dynamics (MD) simulation, a 40 ns enhanced sampling replica exchange molecular dynamics (REMD) simulation, and a reevaluation of the experimental NMR data. The MD and REMD simulations employed the GLYCAM06/AMBER force field with explicit treatment of solvation. The shorter time scale of the REMD sampling method gave rise to RIS and theoretical scalar 3J(HH) coupling constants that were comparable to those from the much longer traditional MD simulation. The 3J(HH) coupling constants computed from the MD methods were in excellent agreement with those observed experimentally. Despite the agreement between the computed and the experimental J-values, there were variations between the rotamer populations computed directly from the MD data and those derived from the experimental NMR data. The experimentally derived populations were determined utilizing limiting J-values from an analysis of NMR data from substituted ethane molecules and may not be completely appropriate for application in more complex molecules, such as glycerol. Here, new limiting J-values have been derived via a combined MD and quantum mechanical approach and were used to decompose the experimental 3J(HH) coupling constants into population distributions for the glycerol RIS.

Positron-annihilation-lifetime response and broadband dielectric relaxation spectroscopy: diethyl phthalate

We report the results of a combined phenomenological analysis of the data from positron-annihilation-lifetime spectroscopy (PALS) and the relaxation data from broadband dielectric spectroscopy (BDS) on diethyl phthalate (DEP). The ortho-positronium (o-Ps) lifetime, tau3, as a function of temperature over a temperature range from 67 K up to 300 K is compared with the spectral features and the relaxation parameters of the BDS spectra decomposed into the primary alpha and the secondary beta processes in the temperature range from 140 K up to 380 K by using the Williams-Watts scheme. Phenomenological model-free analysis of the tau3-T plot provides the three characteristic PALS temperatures, where the two most pronounced ones at TgPALS=185 K and Tb2=245 K=1.32TgPALS are related to the glass-liquid transition and the onset of a quasiplateau region, respectively. In the case of a weaker bend effect at Tb1=210 K=1.14TgPALS, a number of new coincidences with changes in the dielectric beta process have been found. They concern the changes in width parameter of the distribution function for the beta relaxation time and the activation energy of the betaeff process, a crossover from the Arrhenius to the non-Arrhenius type of temperature dependence as well as with the onset of a short-time tail of the beta relaxation time distribution and finally, with changes in the relaxation strength of the alpha and betaeff processes. All these findings indicate a close connection of the o -Ps annihilation parameters and relaxation characteristics of BDS response for the DEP matrix.

Computational analysis of the potential energy surfaces of glycerol in the gas and aqueous phases: effects of level of theory, basis set, and solvation on strongly intramolecularly hydrogen-bonded systems

The 126 possible conformations of 1,2,3-propanetriol (glycerol) have been studied by ab initio molecular orbital and density functional theory calculations in the gas and aqueous phases at multiple levels of theory and basis sets. The partial potential energy surface for glycerol as well as an analysis of the conformational properties and hydrogen-bonding trends in both phases have been obtained. In the gas phase at the G2(MP2) and CBS-QB3 levels of theory, the important, low-energy conformers are structures 100 and 95. In the aqueous phase at the SM5.42/HF/6-31G* level of theory, the lowest energy conformers are structures 95 and 46. Boltzmann distributions have been determined from these high-level calculations, and good agreement is observed when these distributions are compared to the available experimental data. These calculations indicate that the enthalpic and entropic contributions to the Gibbs free energy are important for an accurate determination of the conformational and energetic preferences of glycerol. Different levels of theory and basis sets were used in order to understand the effects of nonbonded interactions (i.e., intramolecular hydrogen bonding). The efficiency of basis set and level of theory in dealing with the issue of intramolecular hydrogen bonding and reproducing the correct energetic and geometrical trends is discussed, especially with relevance to practical computational methods for larger polyhydroxylated compounds, such as oligosaccharides.