화학공학소재연구정보센터
Journal of Chemical Physics, Vol.103, No.6, 2169-2177, 1995
On the Application of Instantaneous Normal-Mode Analysis to Long-Time Dynamics of Liquids
While the applicability of instantaneous normal mode (INM) analysis of liquids to short time dynamics is in principle obvious, its relevance to long time dynamics is not clear. Recent attempts by Keyes and co-workers to apply information obtained from this analysis to self-diffusion in supercooled liquid argon is critically analyzed. By extending the range of frequencies studied we show that both imaginary and real branches of the density of modes are represented better, for large omega, by ln[rho(omega)]similar to omega(2)/T than by ln[rho(omega)]similar to omega(4)/T-2 as advocated by Keyes [J. Chem. Phys. 101, 5081 (1994)]. However, since in the relevant frequency range the two fits almost overlap, the numerical results obtained by Keyes, showing good agreement with the simulation results for self-diffusion in supercooled liquid argon, remain valid even though implications for the frequency dependence of the barrier height distribution change. We also explore other possibilities for extracting information from the INM analysis : (1) The density of "zero force modes," defined as the distribution of normal modes found at the bottom or top of their parabolic potential surfaces, can be computed with no appreciable additional numerical effort. This distribution provides a better representation than the total density of modes for the normal mode distribution at well bottoms and at saddles, however, we find that it makes little difference in quantitativeanalysis. (2) We suggest that the ratio rho(u)(omega)/rho(s)(omega) between the density of modes in the unstable and stable branches provide an estimate for the averaged barrier height distribution for large omega. Using this estimate in a transition state theory calculation of the average hopping time between locally stable liquid configurations and using the resulting time in a calculation of the self-diffusion coefficient yields a very good agreement with results of numerical simulation.