Cross correlations between hard and soft forces: A molecular dynamics study

Rotational dynamics of polyatomic ions in aqueous solutions: From continuum model to mode-coupling theory, aided by computer simulations

Due to the presence of the rotational mode and the distributed surface charges, the dynamical behavior of polyatomic ions in water differs considerably from those of the monatomic ions. However, their fascinating dynamical properties have drawn scant attention. We carry out theoretical and computational studies of a series of well-known polyatomic ions, namely, sulfate, nitrate, and acetate ions. All three ions exhibit different rotational diffusivity, with that of the nitrate ion being considerably larger than the other two. They all defy the hydrodynamic laws of size dependence. Study of the local structure around the ions provides valuable insight into the origin of these differences. We carry out a detailed study of the rotational diffusion of these ions by extensive computer simulation and by using the theoretical approaches of the dielectric friction developed by Fatuzzo-Mason (FM) and Nee-Zwanzig (NZ), and subsequently generalized by Alavi and Waldeck. A critical element of the FM-NZ theory is the decomposition of the total rotational friction, ζ_Rot, into Stokes and dielectric parts. The study shows a dominant role of dielectric friction in the sense that if the ions are made neutral, the nature of diffusion changes and the values become much larger. Our analyses further reveal that the decomposition of total friction into the Stokes and dielectric friction breaks down for sulfate ions but remains semi-quantitatively valid for nitrate and acetate ions. We discuss the relationship between translational and rotational dielectric friction on rigid spherical ions. We develop a self-consistent mode-coupling theory (SC-MCT) formalism that could provide a unified view of rotational friction of polyatomic ions in polar medium. Our SC-MCT shows that the breakdown can be attributed to the change in the microscopic structural features. The mode-coupling theory helps in elucidating the role of coupling between translational and rotational motion of these ions. In fact, these two motions self-consistently determine the value of each other. The reference interaction site model-based MCT suggests an interesting relation between the torque-torque and the force-force time correlation function with the proportionality constant being determined by the geometry and the charge distribution of the polyatomic molecule. We point out several parallelisms between the theories of translational and rotation friction calculations of ions in polar liquids.

SIZE DEPENDENCE OF SOLVATION STRUCTURE AND DYNAMICS OF IONS IN LIQUID N-METHYLACETAMIDE: A MOLECULAR DYNAMICS SIMULATION STUDY

The solvation structure and dynamics of alkali metal (Li+, Na+, K+, Rb+, Cs+) and halide (F-, Cl-, Br-, I-) ions in liquid N -methylacetamide (NMA) are calculated at two different temperatures T = 313 K and 453 K, by using classical molecular dynamics simulations. We have also considered [Formula: see text] and some larger cations such as I +, Me 4 N +, and Et4N+ in this study to investigate the size dependence solvation structure and dynamics of ions in liquid NMA. With the increase of ion size, the self-diffusion coefficients of cations are found to increase and the maximum is observed at Me4N+ , whereas for halide ions the increase of diffusion coefficient with ion size continues up to I- and no maximum is observed. However, the relative increase of the diffusion coefficients of larger ion compared to those of Li+ and F+ are found to be significantly higher at low temperature. Results are very good in agreement with experimental observation.

Molecular dynamics simulation study on the transient response of solvation structure during the translational diffusion of solute

The transient response function of the density profile of the solvent around a solute during the translational diffusion of the solute is formulated based on the generalized Langevin formalism. The resultant theory is applied to both neat Lennard-Jones fluids and cations in liquid water, and the response functions are obtained from the analysis of the molecular dynamics simulations. In the case of the self-diffusion of Lennard-Jones fluids, the responses of the solvation structures are in harmony with conventional pictures based on the mode-coupling theory, that is, the binary collision in the low-density fluids, the backflow effect from medium to high density fluids, and the backscatter effect in the liquids near the triple point. In the case of cations in water, the qualitative behavior is strongly dependent on the size of cations. The pictures similar to simple dense liquids are obtained for the large ion and the neutral molecule, while the solvent waters within the first solvation shell of small ions show an oscillatory response in the short-time region. In particular, the oscillation is remarkably underdumped for lithium ion. The origin of the oscillation is discussed in relation to the theoretical treatment of the translational diffusion of ions in water.