Modeling velocity autocorrelation functions for confined fluids using γ distributions

Using ? distributions to predict self-diffusivities and density of states of fluids confined in carbon nanotubes

The density of states of rare gas atoms confined in carbon nanotubes is analyzed using a recently proposed model based on gamma distributions [Krishnan and Ayappa, J. Chem. Phys., 124 144503 (2006)]. The inputs into the model are the 2nd and 4th frequency moments that are obtained from molecular dynamics simulations. The predicted density of states, velocity autocorrelation functions and self-diffusivities are compared with those obtained from molecular dynamics simulations, for different nanotube loadings and temperatures. All results are reported for argon confined in a (16,16) carbon nanotube. The model predictions are extremely accurate at intermediate reduced densities of rhosigma(3) = 0.3, 0.4, where the majority of the self-diffusivity predictions lie within 10% of the simulation results. Since the frequency moments can be also obtained from Monte Carlo simulations, the study suggests an alternate route to the system dynamics of strongly confined fluids.

Model for dynamics of inhomogeneous and bulk fluids

An accurate model for the density of states (DOS) for strongly inhomogeneous and bulk fluids has been proposed based on gamma distributions. The contribution to the density of states from the collective dynamics is modeled as an incomplete gamma distribution and the high frequency region is obtained from the solution of the memory equation using a sech memory kernel. Using only the frequency moments as input, the model parameters for the collective dynamics are obtained by matching moments of the resulting distribution. The model results in an analytical expression for the self-diffusivity of the fluid. We present results for soft sphere fluids confined in slit-shaped pores as well as bulk soft sphere liquids. Comparisons of the DOS, velocity autocorrelation functions, and memory kernels with molecular dynamics simulations reveal that the model predicts features in the DOS over the entire frequency range and is able to capture changes in the DOS as a function of fluid density and temperature. As a result the predicted VACFs, memory kernels, and self-diffusivities are accurately predicted over a wide range of conditions. Since the frequency moments for bulk liquids can be obtained from pair correlation functions, our method provides a direct route from fluid structure to dynamics. For fluids confined in slit-shaped pores, where the frequency moments are obtained from molecular dynamics simulations, the predicted self-diffusivities capture the resulting oscillations due to variations in the solvation pressure, and in the case of smooth walled pores, the predictions are superior to those obtained using kinetic theory.

Relaxation and Short Time Dynamics of Bulk Liquids and Fluids Confined in Spherical Cavities and Slit Pores

The density of states for bulk and confined fluids have been modeled using a recently proposed gamma distribution (Krishnan, S. H.; Ayappa, K. G. J. Chem. Phys. 2004, 121, 3197). The gamma distribution results in a closed form analytical expression for the velocity autocorrelation function and the relaxation time of the fluid. The two parameters of the gamma distribution are related analytically to the second and fourth frequency moments of the fluid using short time expansions. The predictions by the proposed gamma model are compared with the velocity autocorrelation functions obtained using the theory of instantaneous normal modes (INMs) and from molecular dynamics simulations. The model is applied to a bulk soft sphere liquid and fluids confined in a spherical cavity and slit-shaped pores. The gamma model is able to capture the resulting changes in relaxation time due to changes in density and temperature extremely well for both the bulk liquid and confined inhomogeneous fluid situations. In all cases, the predictions by the gamma model are superior to those obtained from the INM theory. In the case of the fluid confined in a slit pore, the loadings were obtained from a grand canonical Monte Carlo simulation where the pore is equilibrated with a bulk fluid. This is similar to a confinement situation in a surface force apparatus. The predicted relaxation times vs pore widths from the gamma model are seen to accurately capture the oscillations due to formation and disruption of layers within the slit pore.

Ion exchange in reverse micelles

The distribution and dynamics of alkali cations inside Na-AOT reverse micelles have been investigated using Monte Carlo and molecular dynamics simulations. Water is modeled using the extended simple point charge (SPC/E) model. Simulations were carried out for alkali salts of Li+, Na+, K+, and Cs+ placed into the aqueous core of the reverse micelle, for situations corresponding to one and three molecules of added salt. In all cases, we observe that the larger K+ and Cs+ ions exchange with the Na+ counterion; however, the smaller Li+ ion prefers to remains solvated within the core of the reverse micelle. Our study reveals that the oil-water interface of the Na-AOT reverse micelle has the greatest selectivity toward Cs+ followed by K+ and Li+. A model based on enthalpic contributions illustrates that the solvation energies of the different cations in water control the ion-exchange process. The hydration number of the first water shell for Li+ situated in the aqueous core of the reverse micelle with radius R = 14.1 A was similar to that observed at infinite dilution in bulk water.