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

Dipolar Order Controls Dielectric Response of Glass-Forming Liquids

The dielectric response of liquids reflects both reorientation of single molecular dipoles and collective modes, i.e., dipolar cross-correlations. A recent theory predicts the latter to produce an additional slow peak in the dielectric loss spectrum. Following this idea we argue that in supercooled liquids the high-frequency power law exponent of the dielectric loss β should be correlated with the degree of dipolar order, i.e., the Kirkwood correlation factor g_{K}. This notion is confirmed for 25 supercooled liquids. While our findings support recent theoretical work the results are shown to violate the earlier Kivelson-Madden theory.

Theory and modeling of molecular modes in the NMR relaxation of fluids

Traditional theories of the nuclear magnetic resonance (NMR) autocorrelation function for intra-molecular dipole pairs assume a single-exponential decay, yet the calculated autocorrelation of realistic systems displays a rich, multi-exponential behavior, resulting in anomalous NMR relaxation dispersion (i.e., frequency dependence). We develop an approach to model and interpret the multi-exponential intra-molecular autocorrelation using simple, physical models within a rigorous statistical mechanical development that encompasses both rotational diffusion and translational diffusion in the same framework. We recast the problem of evaluating the autocorrelation in terms of averaging over a diffusion propagator whose evolution is described by a Fokker-Planck equation. The time-independent part admits an eigenfunction expansion, allowing us to write the propagator as a sum over modes. Each mode has a spatial part that depends on the specified eigenfunction and a temporal part that depends on the corresponding eigenvalue (i.e., correlation time) with a simple, exponential decay. The spatial part is a probability distribution of the dipole pair, analogous to the stationary states of a quantum harmonic oscillator. Drawing inspiration from the idea of inherent structures in liquids, we interpret each of the spatial contributions as a specific molecular mode. These modes can be used to model and predict the NMR dipole-dipole relaxation dispersion of fluids by incorporating phenomena on the molecular level. We validate our statistical mechanical description of the distribution in molecular modes with molecular dynamics simulations interpreted without any relaxation models or adjustable parameters: the most important poles in the Padé-Laplace transform of the simulated autocorrelation agree with the eigenvalues predicted by the theory.

Ultrafast Proton Conduction in an Aqueous Electrolyte Confined in Adamantane-like Micropores of a Sulfonated, Aromatic Framework

Sulfonated, cross-linked porous polymers are promising frameworks for aqueous high-performance electrolyte-host systems for electrochemical energy storage and conversion. The systems offer high proton conductivities, excellent chemical and mechanical stabilities, and straightforward water management. However, little is known about mass transport mechanisms in such nanostructured hosts. We report on the synthesis and postsynthetic sulfonation of an aromatic framework (SPAF-2) with a 3D-interconnected nanoporosity and varying sulfonation degrees. Water adsorption produces the system SPAF-2H20. It features proton exchange capacities up to 6 mequiv g-1 and exceptional proton conductivities of about 1 S cm-1. Two contributions are essential for the highly efficient transport. First, the nanometer-sized pores link the charge transport to the diffusion of adsorbed water molecules, which is almost as fast as bulk water. Second, continuous exchange between interface-bound and mobile species enhances the conductivities at elevated temperatures. SPAF-2H20 showcases how to tailor nanostructured electrolyte-host systems with liquid-like conductivities.

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.

NMR Investigation of Water in Salt Crusts: Insights from Experiments and Molecular Simulations

The evaporation of water from bare soil is often accompanied by the formation of a layer of crystallized salt, a process that must be understood in order to address the issue of soil salinization. Here, we use nuclear magnetic relaxation dispersion measurements to better understand the dynamic properties of water within two types of salt crusts: sodium chloride (NaCl) and sodium sulfate (Na₂SO₄). Our experimental results display a stronger dispersion of the relaxation time T₁ with frequency for the case of sodium sulfate as compared to sodium chloride salt crusts. To gain insight into these results, we perform molecular dynamics simulations of salt solutions confined within slit nanopores made of either NaCl or Na₂SO₄. We find a strong dependence of the value of the relaxation time T₁ on pore size and salt concentration. Our simulations reveal the complex interplay between the adsorption of ions at the solid surface, the structure of water near the interface, and the dispersion of T₁ at low frequency, which we attribute to adsorption-desorption events.

Molecular Dynamics of Ionic Liquids from Fast-Field Cycling NMR and Molecular Dynamics Simulations

Understanding the connection between the molecular structure of ionic liquids and their properties is of paramount importance for practical applications. However, this connection can only be established if a broad range of physicochemical properties on different length and time scales is already available. Even then, the interpretation of the results often remains ambiguous due to the natural limits of experimental approaches. Here we use fast-field cycling (FFC) to access both translational and rotational dynamics of ionic liquids. These combined with a comprehensive physicochemical characterization and MD simulations provide a toolkit to give insight into the mechanisms of molecular mechanics. The FFC results are consistent with the computer simulation and conventional physicochemical approaches. We show that curling of the side chains around the positively charged cationic core is essential for the properties of ether-functionalized ionic liquids, and we demonstrate that neither geometry nor polarity alone are sufficient to explain the macroscopic properties.

Glassy dynamics in polyalcohols: intermolecular simplicity vs. intramolecular complexity

Using depolarized light scattering, we have recently shown that structural relaxation in a broad range of supercooled liquids follows, to good approximation, a generic line shape with high-frequency power law ω^-1/2. We now continue this study by investigating a systematic series of polyalcohols (PAs), frequently used as model-systems in glass-science, i.a., because the width of their respective dielectric loss spectra varies strongly along the series. Our results reveal that the microscopic origin of the observed relaxation behavior varies significantly between different PAs: while short-chained PAs like glycerol rotate as more or less rigid entities and their light scattering spectra follow the generic shape, long-chained PAs like sorbitol display pronounced intramolecular dynamic contributions on the time scale of structural relaxation, leading to systematic deviations from the generic shape. Based on these findings we discuss an important limitation for observing the generic shape in a supercooled liquid: the dynamics that is probed needs to reflect the intermolecular dynamic heterogeneity, and must not be superimposed by effects of intramolecular dynamic heterogeneity.

A Relation between the Formation of a Hydrogen-Bond Network and a Time-Scale Separation of Translation and Rotation in Molecular Liquids

We study the relation between the translational and rotational motions of liquids, which is anticipated in the framework of the Stokes-Einstein-Debye (SED) treatment. For this purpose, we exploit the fact that ¹H field-cycling nuclear magnetic resonance relaxometry and molecular dynamics simulations provide access to both modes of motion. The experimental and computational findings are fully consistent and show that the time-scale separation between translation and rotation increases from the van der Waals liquid o-terphenyl over ethylene glycol to the hydrogen-bonded liquid glycerol, indicating an increasing degree of breakdown of the SED relation. Furthermore, the simulation results for two ethylene glycol models with different molecular conformations indicate that the translation is more retarded than the rotation when the density of intermolecular hydrogen bonds increases. We conclude that an increasing connectivity of a hydrogen-bond network leads to an increasing time-scale separation and, thus, to a stronger SED violation.

NMR Relaxometry Accessing the Relaxation Spectrum in Molecular Glass Formers

It is a longstanding question whether universality or specificity characterize the molecular dynamics underlying the glass transition of liquids. In particular, there is an ongoing debate to what degree the shape of dynamical susceptibilities is common to various molecular glass formers. Traditionally, results from dielectric spectroscopy and light scattering have dominated the discussion. Here, we show that nuclear magnetic resonance (NMR), primarily field-cycling relaxometry, has evolved into a valuable method, which provides access to both translational and rotational motions, depending on the probe nucleus. A comparison of ¹H NMR results indicates that translation is more retarded with respect to rotation for liquids with fully established hydrogen-bond networks; however, the effect is not related to the slow Debye process of, for example, monohydroxy alcohols. As for the reorientation dynamics, the NMR susceptibilities of the structural (α) relaxation usually resemble those of light scattering, while the dielectric spectra of especially polar liquids have a different broadening, likely due to contributions from cross correlations between different molecules. Moreover, NMR relaxometry confirms that the excess wing on the high-frequency flank of the α-process is a generic relaxation feature of liquids approaching the glass transition. However, the relevance of this feature generally differs between various methods, possibly because of their different sensitivities to small-amplitude motions. As a major advantage, NMR is isotope specific; hence, it enables selective studies on a particular molecular entity or a particular component of a liquid mixture. Exploiting these possibilities, we show that the characteristic Cole-Davidson shape of the α-relaxation is retained in various ionic liquids and salt solutions, but the width parameter may differ for the components. In contrast, the low-frequency flank of the α-relaxation can be notably broadened for liquids in nanoscopic confinements. This effect also occurs in liquid mixtures with a prominent dynamical disparity in their components.

Nuclear Spin Relaxation in Viscous Liquids: Relaxation Stretching of Single-Particle Probes

Spin-lattice relaxation rates R₁(ω,T), probed via high-field and field-cycling nuclear magnetic resonance (NMR), are used to test the validity of frequency-temperature superposition (FTS) for the reorientation dynamics in viscous liquids. For several liquids, FTS is found to apply so that master curves can be generated. The susceptibility spectra are highly similar to those obtained from depolarized light scattering (DLS) and reveal an excess wing. Where FTS works, two approaches are suggested to access the susceptibility: (i) a plot of deuteron R₁(T) vs the spin-spin relaxation rate R₂(T) and (ii) a plot of R₁(T) vs an independently measured reference time τ_ref(T). Using single-frequency scans, (i) allows one to extract the relaxation stretching as well as the NMR coupling constant. Surveying 26 data sets, we find Kohlrausch functions with exponents 0.39 < β_K ≤ 0.67. Plots of the spin-spin relaxation rate R₂─rescaled by the NMR coupling constant─as a function of temperature allow one to test how well site-specific NMR relaxations couple to a given reference process. Upon cooling of flexible molecule liquids, the site-specific dynamics is found to merge, suggesting that near T_g the molecules reorient essentially as a rigid entity. This presents a possible resolution for the much lower stretching parameters reported here at high temperatures that contrast with the ones that were reported to be universal in a recent DLS study close to T_g. Our analysis underlines that deuteron relaxation is a uniquely powerful tool to probe single-particle reorientation.

Molecular and micro-scale heterogeneities in Raman modes of a relaxing polymer glass

We have used Raman spectroscopy to study relaxation dynamics at two different length scales, molecular level and micro-scale in order to probe the presence of cooperative rearranging regions in a polymer glass. Response to slow thermal cycles and fast quench through the glass transition temperature (T_g) is analyzed for film and unprocessed forms of polyvinyl acetate (PVAc). In PVAc film, enhanced disorder and molecular mobility lead to peak broadening by about a factor of 10 compared to unprocessed PVAc. Thermal cycles (10 K min^-1) produce hysteresis in integrated Raman peak intensity (loop areaAINTI).AINTIvalues of film are two orders of magnitude more than unprocessed, indicating more configurational mosaics with higher interfacial energy dissipations. Ageing after 60 K min^-1quench manifests as heterogeneous molecular dynamics of film Raman modes with significant peak-width variations, differentiating high mobility and low mobility modes. Two-dimensional mapping of film Raman modes after quench reveal micro-scale clusters of average size ≈250 molecules having fractal boundaries with fractal dimensiond_f= 1.5, resemblingd_fof percolation clusters below percolation threshold. During thermal cycling and relaxation after a quench, cooperative segmental dynamics with large correlations between skeletal C-C stretch and side branch modes is observed. The observations are analyzed in the context of the random first order transition theory of glasses, which attributes heterogeneous relaxations in glasses to the presence of clusters of variable configurational states.