Molecular dynamics of liquids modelled by ‘2-Lennard-Jones centres’ pair potentials

Understanding dynamics in coarse-grained models. IV. Connection of fine-grained and coarse-grained dynamics with the Stokes-Einstein and Stokes-Einstein-Debye relations

Applying an excess entropy scaling formalism to the coarse-grained (CG) dynamics of liquids, we discovered that missing rotational motions during the CG process are responsible for artificially accelerated CG dynamics. In the context of the dynamic representability between the fine-grained (FG) and CG dynamics, this work introduces the well-known Stokes-Einstein and Stokes-Einstein-Debye relations to unravel the rotational dynamics underlying FG trajectories, thereby allowing for an indirect evaluation of the effective rotations based only on the translational information at the reduced CG resolution. Since the representability issue in CG modeling limits a direct evaluation of the shear stress appearing in the Stokes-Einstein and Stokes-Einstein-Debye relations, we introduce a translational relaxation time as a proxy to employ these relations, and we demonstrate that these relations hold for the ambient conditions studied in our series of work. Additional theoretical links to our previous work are also established. First, we demonstrate that the effective hard sphere radius determined by the classical perturbation theory can approximate the complex hydrodynamic radius value reasonably well. Furthermore, we present a simple derivation of an excess entropy scaling relationship for viscosity by estimating the elliptical integral of molecules. In turn, since the translational and rotational motions at the FG level are correlated to each other, we conclude that the "entropy-free" CG diffusion only depends on the shape of the reference molecule. Our results and analyses impart an alternative way of recovering the FG diffusion from the CG description by coupling the translational and rotational motions at the hydrodynamic level.

Understanding dynamics in coarse-grained models. III. Roles of rotational motion and translation-rotation coupling in coarse-grained dynamics

This paper series aims to establish a complete correspondence between fine-grained (FG) and coarse-grained (CG) dynamics by way of excess entropy scaling (introduced in Paper I). While Paper II successfully captured translational motions in CG systems using a hard sphere mapping, the absence of rotational motions in single-site CG models introduces differences between FG and CG dynamics. In this third paper, our objective is to faithfully recover atomistic diffusion coefficients from CG dynamics by incorporating rotational dynamics. By extracting FG rotational diffusion, we unravel, for the first time reported to our knowledge, a universality in excess entropy scaling between the rotational and translational diffusion. Once the missing rotational dynamics are integrated into the CG translational dynamics, an effective translation-rotation coupling becomes essential. We propose two different approaches for estimating this coupling parameter: the rough hard sphere theory with acentric factor (temperature-independent) or the rough Lennard-Jones model with CG attractions (temperature-dependent). Altogether, we demonstrate that FG diffusion coefficients can be recovered from CG diffusion coefficients by (1) incorporating "entropy-free" rotational diffusion with translation-rotation coupling and (2) recapturing the missing entropy. Our findings shed light on the fundamental relationship between FG and CG dynamics in molecular fluids.

Molecular insight into the microstructure and microscopic dynamics of pyridinium ionic liquids with different alkyl chains based on temperature response

High temperature is advantageous to the aggregation of the polar regions as well as the nonpolar regions of pyridinium ionic liquids.

Structure and dynamics of liquid CS2: Going from ambient to elevated pressure conditions

Molecular dynamics simulation studies were performed to investigate the structural and dynamic properties of liquid carbon disulfide (CS₂) from ambient to elevated pressure conditions. The results obtained have revealed structural changes at high pressures, which are related to the more dense packing of the molecules inside the first solvation shell. The calculated neutron and X-ray structure factors have been compared with available experimental diffraction data, also revealing the pressure effects on the short-range structure of the liquid. The pressure effects on the translational, reorientational, and residence dynamics are very strong, revealing a significant slowing down when going from ambient pressure to 1.2 GPa. The translational dynamics of the linear CS₂ molecules have been found to be more anisotropic at elevated pressures, where cage effects and librational motions are reflected on the shape of the calculated time correlation functions and their corresponding spectral densities.

Analysis of the transverse and the longitudinal pseudodiffusion of CO2 in sub- and supercritical states: a molecular-dynamics analysis

We have performed molecular-dynamics simulations of CO(2) system along the gas-liquid coexistence curve and on the isochore 94.22 cm(3) mol(-1) (which corresponds to the critical isochore). The calculation has been carried out in order to analyze the diffusion of CO(2) and particularly to figure out how the diffusion coefficient may be decomposed along the molecular axes. This makes it possible to analyze the anisotropy of the diffusion along these axes and to shed light on the microscopic changes which accompany such behavior. This anisotropy is traced back to the effect of the translation-rotation coupling (TRC) along the molecular axes. Along the liquid-gas coexistence curve, the pseudolongitudinal diffusion is found to be more rapid than the transverse one. The opposite trend is found along the isochore 94.22 cm(3) mol(-1). The role of the local structure was explored by calculating intermediate scattering function and the autocorrelation functions for the forces acting along the molecular axes. It is shown that the strength of the TRC effect is correlated to the difference between the relaxation times of the local structure, that of the reorientation along the molecular axes, and that of the translational motion. The analysis of the correlation time and the average mean square force along the longitudinal and transverse directions confirms the anisotropy of the local environment that determines the translational dynamics of a molecule.

Computer Simulation study of polar liquids: static and dynamic properties

The Berne-Pechukas-Kushick potential for linear molecules, together with a central electric point dipole, has been used as the model for polar liquids in a series of molecular dynamics computer simulations. Different simulations used Ewald summation, spherical cut off and spherical cut off with reaction field correction to sum the dipole-dipole interactions, and the relative merits of these methods are discussed. Results for thermodynamic properties and the distribution function h ∆ ( r ) are presented. Results for the autocorrelation functions for dipole orientation and angular velocity are also presented, and discussed in terms of memory functions, particularly with regard to the Steele perturbation theory. The complex, frequency-dependent dielectric constant has been calculated from the memory function of the autocorrelation function of the periodic cube dipole moment, and results are presented in the form of Cole-Cole plots. The ‘dielectric friction’ approach to polar liquids is examined against the evidence of the molecular dynamics data. Attention is drawn to the behaviour of the cross correlations between dipole orientations and their importance to the problem of dielectric relaxation is emphasized.