Microscopic dynamics of molecular liquids and glasses: role of orientations and translation-rotation coupling

Diffusion theory of molecular liquids in the energy representation and application to solvation dynamics

The generalized Langevin equation (GLE) formalism is a useful theoretical fundament for analyzing dynamical phenomena rigorously. Despite the systematic formulation of dynamics theories with practical approximations, however, the applicability of GLE-based methods is still limited to simple polyatomic liquids due to the approximate treatment of molecular orientations involved in the static molecular liquid theory. Here, we propose an exact framework of dynamics based on the GLE formalism incorporating the energy representation theory of solution, an alternative static molecular liquid theory. A fundamental idea is the projection of the relative positions and orientations of solvents around a solute onto the solute-solvent interaction, namely the energy coordinate, enabling us to describe the dynamics on a one-dimensional coordinate. Introducing systematic approximations, such as the overdamped limit, leads to the molecular diffusion equation in the energy representation that is described in terms of the distribution function of solvents on the energy coordinate and the diffusion coefficients. The present theory is applied to the solvation dynamics triggered by the photoexcitation of benzonitrile. The long-time behavior of the solvation time correlation function is in good agreement with that obtained by the molecular dynamics simulation.

Reorientational dynamics in highly asymmetric binary low-molecular mixtures-A quantitative comparison of dielectric and NMR spectroscopy results

Previously, we scrutinized the dielectric spectra of a binary glass former made by a low-molecular high-T_g component 2-(m-tertbutylphenyl)-2'-tertbutyl-9,9'-spirobi[9H]fluorene (m-TPTS; T_g = 350 K) and low-T_g tripropyl phosphate (TPP; T_g = 134 K) [Körber et al., Phys. Chem. Chem. Phys. 23, 7200 (2021)]. Here, we analyze nuclear magnetic resonance (NMR) spectra and stimulated echo decays of deuterated m-TPTS-d₄ (²H) and TPP (³¹P) and attempt to understand the dielectric spectra in terms of component specific dynamics. The high-T_g component (α₁) shows relaxation similar to that of neat systems, yet with some broadening upon mixing. This correlates with high-frequency broadening of the dielectric spectra. The low-T_g component (α₂) exhibits highly stretched relaxations and strong dynamic heterogeneities indicated by "two-phase" spectra, reflecting varying fractions of fast and slow liquid-like reorienting molecules. Missing for the high-T_g component, such two-phase spectra are identified down to w_TPP = 0.04, indicating that isotropic reorientation prevails in the rigid high-T_g matrix stretching from close to T_g ^TPP to T_g1 w_TPP. This correlates with low-frequency broadening of the dielectric spectra. Two T_g values are defined: T_g1 (w_TPP) displays a plasticizer effect, whereas T_g2 (w_TPP) passes through a maximum, signaling extreme separation of the component dynamics at low w_TPP. We suggest understanding the latter counter-intuitive feature by referring to a crossover from "single glass" to "double glass" scenario revealed by recent MD simulations. Analyses reveal that a second population of TPP molecules exists, which is associated with the dynamics of the high-T_g component. However, the fractions are lower than suggested by the dielectric spectra. We discuss this discrepancy considering the role of collective dynamics probed by dielectric but not by NMR spectroscopy.

Dynamics theory for molecular liquids based on an interaction site model

Dynamics theories for molecular liquids based on an interaction site model have been developed over the past few decades and proved to be powerful tools to investigate various dynamical phenomena.

Self-consistent generalized Langevin-equation theory for liquids of nonspherically interacting particles

A self-consistent generalized Langevin-equation theory is proposed to describe the self- and collective dynamics of a liquid of linear Brownian particles. The equations of motion for the spherical harmonics projections of the collective and self-intermediate-scattering functions, F_{lm,lm}(k,t) and F_{lm,lm}^{S}(k,t), are derived as a contraction of the description involving the stochastic equations of the corresponding tensorial one-particle density n_{lm}(k,t) and the translational (α=T) and rotational (α=R) current densities j_{lm}^{α}(k,t). Similar to the spherical case, these dynamic equations require as an external input the equilibrium structural properties of the system contained in the projections of the static structure factor, denoted by S_{lm,lm}(k). Complementing these exact equations with simple (Vineyard-like) approximate relations for the collective and the self-memory functions we propose a closed self-consistent set of equations for the dynamic properties involved. In the long-time asymptotic limit, these equations become the so-called bifurcation equations, whose solutions (the nonergodicity parameters) can be written, extending the spherical case, in terms of one translational and one orientational scalar dynamic order parameter, γ_{T} and γ_{R}, which characterize the possible dynamical arrest transitions of the system. As a concrete illustrative application of this theory we determine the dynamic arrest diagram of the dipolar hard-sphere fluid. In qualitative agreement with mode coupling theory, the present self-consistent equations also predict three different regions in the state space spanned by the macroscopic control parameters η (volume fraction) and T* (scaled temperature): a region of fully ergodic states, a region of mixed states, in which the translational degrees of freedom become arrested while the orientational degrees of freedom remain ergodic, and a region of fully nonergodic states.

Dynamical density functional theory for anisotropic colloidal particles

We generalize the formalism of dynamical density functional theory for translational Brownian dynamics toward that of anisotropic colloidal particles which perform both translational and rotational Brownian motion. Using a mean-field approximation for the density functional and a Gaussian-segment model for the rod interaction, the dynamical density functional theory is then applied to a concentrated rod suspension in a confined slab geometry made by two parallel soft walls. The walls are either expanded or compressed and the relaxation behavior is investigated for an equilibrated starting configuration. We find distinctly different orientational ordering during expansion and compression. During expansion we observe preferential parallel ordering of the rods relative to the wall while during compression there is homeotropic ordering perpendicular to the wall. We find a nonexponential relaxation behavior in time. Furthermore, an external field which aligns the rods perpendicular to the walls is turned on or switched off and similar differences in the relaxational dynamics are found. Comparing the theoretical predictions to Brownian dynamics computer simulation data, we find good agreement.

Dynamic arrest in a liquid of symmetric dumbbells: Reorientational hopping for small molecular elongations

We present extensive equilibrium and out-of-equilibrium molecular-dynamics simulations of a liquid of symmetric dumbbell molecules, for constant packing fraction, as a function of temperature and molecular elongation. We compute diffusion constants as well as odd and even orientational correlators. The notations odd and even refer to the parity of the order l of the corresponding Legendre l polynomial, evaluated for the orientation of the molecular axis relative to its initial position. Rotational degrees of freedom of order l are arrested if, in the long-time limit, the corresponding orientational l correlator does not decay to zero. It is found that for large elongations translational and rotational degrees of freedom freeze at the same temperature. For small elongations only the even rotational degrees of freedom remain coupled to translational motions and arrest at a finite common temperature. On the contrary, the odd rotational degrees of freedom remain ergodic at all investigated temperatures. Hence, in the translationally arrested state, each molecule remains trapped in the cage formed by its neighboring molecules, but is able to perform 180 degrees rotations, which lead to relaxation only for the odd orientational correlators. The temperature dependence of the characteristic time of these residual rotations is well described by an Arrhenius law. Finally, we discuss the evidence in favor of the presence of the type-A transition for the odd rotational degrees of freedom, as predicted by mode-coupling theory for small molecular elongations. This transition is distinct from the type-B transition, associated with the arrest of the translational and even rotational degrees of freedom for small elongations, and with all degrees of freedom for large elongations. Odd orientational correlators are computed for small elongations at very low temperatures in the translationally arrested state. The obtained results suggest that hopping events restore the ergodicity of the odd rotational degrees of freedom at temperatures far below the A transition.

Direct experimental assessment of the strength of orientational correlations in polar liquids

The strength of molecular orientational correlations in polar liquids is assessed by means of comparison of the diffuse scattering patterns of a liquid composed by molecules devoid of permanent electric dipole but having a weak quadrupole moment and those for a liquid composed by permanent molecular dipoles. The extent of orientational correlations within the liquid phases is in both cases assessed by comparison of the liquid radial distributions to those present in the rotator-phase (plastic) crystal phases of both compounds. For such disordered-crystal phases, information concerning orientational correlations is directly derived from the experimental scattering patterns by means of analysis of the diffuse scattering background present beneath the Bragg peaks. The results show that rather than long-ranged, orientational correlations in polar or polarizable liquids are confined within distances comprising the second coordination sphere.

Structural relaxation in supercooled orthoterphenyl

We report molecular-dynamics simulation results performed for a model of molecular liquid orthoterphenyl in supercooled states, which we then compare with both experimental data and mode-coupling-theory (MCT) predictions, aiming at a better understanding of structural relaxation in orthoterphenyl. We pay special attention to the wave number dependence of the collective dynamics. It is shown that the simulation results for the model share many features with experimental data for real system, and that MCT captures the simulation results at the semiquantitative level except for intermediate wave numbers connected to the overall size of the molecule. Theoretical results at the intermediate wave number region are found to be improved by taking into account the spatial correlation of the molecule's geometrical center. This supports the idea that unusual dynamical properties at the intermediate wave numbers, reported previously in simulation studies for the model and discernible in coherent neutron-scattering experimental data, are basically due to the coupling of the rotational motion to the geometrical-center dynamics. However, there still remain qualitative as well as quantitative discrepancies between theoretical prediction and corresponding simulation results at the intermediate wave numbers, which call for further theoretical investigation.

Structural relaxation in a system of dumbbell molecules

The interaction-site-density-fluctuation correlators, the dipole-relaxation functions, and the mean-squared displacements of a system of symmetric dumbbells of fused hard spheres are calculated for two representative elongations of the molecules within the mode-coupling theory for the evolution of glassy dynamics. For large elongations, universal relaxation laws for states near the glass transition are valid for parameters and time intervals similar to the ones found for the hard-sphere system. Rotation-translation coupling leads to an enlarged crossover interval for the mean-squared displacement of the constituent atoms between the end of the von Schweidler regime and the beginning of the diffusion process. For small elongations, the superposition principle for the reorientational alpha process is violated for parameters and time intervals of interest for data analysis, and there is a strong breaking of the coupling of the alpha-relaxation scale for the diffusion process with that for representative density fluctuations and for dipole reorientations.

Idealized glass transitions for a system of dumbbell molecules

The mode-coupling theory for ideal glass transitions in simple systems is generalized to a theory for the glassy dynamics of molecular liquids using the density fluctuations of the sites of the molecule's constituent atoms as the basic structure variables. The theory is applied to calculate the liquid-glass phase diagram and the form factors for the arrested structure of a system of symmetric dumbbells of fused hard spheres. The static structure factors, which enter the equations of motion as input, are calculated as function of the packing fraction phi and the molecule's elongation zeta within the reference-interaction-site-model and Percus-Yevick theories. The critical packing fraction phi(c) for the glass transition is obtained as nonmonotone function of zeta with a maximum near zeta=0.43. A transition line is calculated separating a small-zeta-glass phase with ergodic dipole motion from a large-zeta-glass phase where also the reorientational motion is arrested. The Debye-Waller factors at the transition are found to be somewhat larger for sufficiently elongated systems than those for the simple hard-sphere system, but the wave-number dependence of the glass-form factors is quite similar. The dipole reorientations for zeta> or =0.6 are arrested as strongly as density fluctuations with wave vectors at the position of the first sharp diffraction peak.