1
|
Lad KN, Patel MK, Pratap A. Brownian motion with time-dependent friction and single-particle dynamics in liquids. Phys Rev E 2022; 105:064107. [PMID: 35854483 DOI: 10.1103/physreve.105.064107] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2022] [Accepted: 05/18/2022] [Indexed: 11/07/2022]
Abstract
A microscopic theory of molecular motion in classical monatomic liquids, proposed by Glass and Rice [Phys. Rev. 176, 239 (1968)10.1103/PhysRev.176.239], is revisited and extended to incorporate the dynamic friction in the Brownian description of the atomic diffusion in a mean-time-dependent harmonic force field. A modified, non-Markovian Langevin equation is utilized to derive an equation of motion for the velocity autocorrelation function with time-dependent friction coefficient. Numerical solution of the equation gives an excellent account of the velocity autocorrelation function in Lennard-Jones liquids, liquid alkali, and transition metals over a broad range of density and temperature. Derivation of the equation of motion leads to a self-consistent expression for the time dependence of friction coefficient. Our results demonstrate that the nature of time dependence of the friction coefficient changes dramatically with the liquid density. At low and moderate densities, the dynamic friction decays exponentially whereas it increases exponentially at high liquid densities. Our findings provide an opportunity for a different outlook of the Brownian description of atomic dynamics in liquids.
Collapse
Affiliation(s)
- Kirit N Lad
- Computational Condensed Matter Physics Laboratory, Department of Physics, Sardar Patel University, Vallabh Vidya Nagar-388120, Gujarat, India
| | - Margi K Patel
- Condensed Matter Physics Laboratory, Applied Physics Department, Faculty of Technology & Engineering, The M. S. University of Baroda, Vadodara-390001, Gujarat, India
| | - Arun Pratap
- Condensed Matter Physics Laboratory, Applied Physics Department, Faculty of Technology & Engineering, The M. S. University of Baroda, Vadodara-390001, Gujarat, India
| |
Collapse
|
2
|
Huang BC, Chang CH. Localization-delocalization transition of the instantaneous normal modes of liquid water. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 88:042116. [PMID: 24229125 DOI: 10.1103/physreve.88.042116] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/10/2013] [Revised: 09/06/2013] [Indexed: 06/02/2023]
Abstract
Despite the fact that the localization-delocalization transition (LDT) widely exists in wave systems, quantitative studies on its critical and multifractal properties are mainly focused on solids. In this work, these properties are investigated on the vibrational motions of liquid water. Simulations of up to 18000 molecules on the flexible extended simple point charge water model provide nearly 10(6) instantaneous normal modes. They are shown to undergo an LDT close to the translational transition and exhibit multifractal fluctuations while approaching the LDT. In combination with finite-size scaling, multifractal analysis predicts the critical frequency Im(ω(c))≈-131.6 cm(-1) for unstable modes at room temperature. The estimated critical exponent ν≈1.60 is close to those of other calculated systems in the same Wigner-Dyson class. At the LDT, the fractal spectrum f(α) and the most probable local vibrational intensity α(mc)≈4.04 coincide with those of the Anderson model, which might be additional universal properties of LDT in more general wave systems. The results extend the validity of the multifractal scaling approach beyond Andersonian systems to a Hessian system.
Collapse
Affiliation(s)
- Ban-Chiech Huang
- Institute of Physics, National Chiao Tung University, Hsinchu 300, Taiwan
| | | |
Collapse
|
3
|
Li H, Kobrak MN. Instantaneous Normal Mode Analysis of a Series of Model Molten Salts. Chemphyschem 2012; 13:1934-41. [DOI: 10.1002/cphc.201100727] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2011] [Revised: 10/27/2011] [Indexed: 11/08/2022]
|
4
|
Krishnan SH, Ayappa KG. Modeling velocity autocorrelation functions for confined fluids using γ distributions. J Chem Phys 2004; 121:3197-205. [PMID: 15291631 DOI: 10.1063/1.1768939] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We propose a model for the short-time dynamics of fluids confined in slit-shaped pores. The model has been developed from the observation that the real lobe of the instantaneous normal mode density of states (INM DOS) closely follows a gamma distribution. By proposing that the density of states of the confined fluid can be represented by a gamma distribution, the resulting velocity autocorrelation function (VACF) is constructed such that it is accurate upto the fourth frequency moment. The proposed model results in an analytical expression for the VACF and relaxation times. The VACFs obtained from the model have been compared with the VACFs obtained from molecular dynamic simulations and INM analysis for fluids confined in slit-shaped pores over a wide range of confinement and temperatures. The model is seen to capture the short-time behavior of the VACF extremely accurately and in this region is superior to the predictions of the VACF obtained from the real lobe of the INM DOS. Although the model predicts a zero self-diffusivity, the predicted relaxation times are in better agreement with the molecular dynamics results when compared with those obtained from the INM theory.
Collapse
Affiliation(s)
- S H Krishnan
- Department of Chemical Engineering, Indian Institute of Science, Bangalore 560012
| | | |
Collapse
|
5
|
Cavagna A, Giardina I, Grigera TS. A single saddle model for the -relaxation in supercooled liquids. ACTA ACUST UNITED AC 2003. [DOI: 10.1088/0305-4470/36/43/004] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
6
|
De Santis A, Ercoli A, Rocca D. Relative pair dynamics in simple supercooled liquids: longitudinal contributions. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2003; 68:031205. [PMID: 14524756 DOI: 10.1103/physreve.68.031205] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/12/2002] [Indexed: 05/24/2023]
Abstract
The pair dynamics of simulated argon samples is investigated at the melting (85 K), supercooled (55 K), and quenched (20 K) liquid states, and in the crystal (20 K) state. Tagged pairs, initially lying in a given shell, were divided into incoming and outgoing groups and followed along simulated trajectories. Over them, specific correlation functions deltaB(r(0);t), involving the pair separation vector projected along its initial value r(0) (longitudinal dynamics), have been evaluated. More or less pronounced oscillations are detected according to the temperature of the thermodynamic states (and, obviously, their solid or liquid nature); for each state, they depend on the initial pair distance r(0), too. The oscillations vanish after few picoseconds (fast dynamics) in the case of crystal, whereas in the supercooled liquid they decay towards a plateau, whose height increases with the temperature. It is shown that the power spectrum of deltaB(r(0);t) practically yields the same density of states (DOS) produced by the pair velocity correlation function. The deltaB(r(0);t) functions obtained from the argon crystal at 20 K produce DOS curves dominated by two main frequency contributions, at about 40 and 60 cm(-1) (Einstein and Debye frequency, respectively). Their shape is quite well reproduced by damped harmonic oscillator-like (DHO) functions vibrating at that frequencies. In liquid states, the deltaB(r(0);t) plateau, that forms after the fast DHO dynamics, accounts for the system diffusivity. The relaxation towards the plateau is modeled by an exponential function whose decay time is comparable with the average vibration period. Evidence that the liquid states conserve a certain memory of the vibrational modes of the crystal is obtained. In these states, the DHO functions at the Einstein and Debye frequencies plus an exponential function cannot reproduce the deltaB(r(0);t) shape. A pronounced shoulder, that forms around 0.5 ps, requires the contribution of a third DHO. In the DOS, it yields a band centered below 20 cm(-1) that produces low frequency DOS excess in comparison with the DOS of the crystal. This contribution is present in liquid and supercooled high temperature states and survives near the temperature of the glass transition whereas the diffusion practically vanishes.
Collapse
Affiliation(s)
- Alberto De Santis
- Facoltà di Agraria, DISA and INFM-Viterbo, Università della Tuscia, Via San Camillo De Lellis, 01100 Viterbo, Italy
| | | | | |
Collapse
|
7
|
Agrawal PM, Rice BM, Thompson DL. Molecular dynamics study of the effects of voids and pressure in defect-nucleated melting simulations. J Chem Phys 2003. [DOI: 10.1063/1.1570815] [Citation(s) in RCA: 62] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
8
|
Tanaka H, Mohanty U. On the Debye-Waller factor of hexagonal ice: a computer simulation study. J Am Chem Soc 2002; 124:8085-9. [PMID: 12095352 DOI: 10.1021/ja011927h] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We investigate by molecular dynamics (MD) simulations the temperature dependence of the Debye-Waller (DW) factor of hexagonal ice with 25 different proton-disordered configurations. Each initial configuration is composed of 288 water molecules with no net dipole moment. The intermolecular interaction of water is described by TIP4P potential. Each production run of the simulation is 15 ns or longer. We observe a change in slope of the DW factor around 200 K, which cannot be explained within the framework of either classical or quantum harmonic approximation. Configurations generated by MD simulations are subjected to the steepest descent energy minimization. Analysis of the local energy minimum structures reveals that water molecules above 200 K jump to other lattice sites via some local energy minimum structures which contain some water molecules sitting on the locations other than the lattice sites. As time evolves, these defect molecules move back and forth to the lattice sites yielding defect-free structures. Those motions are responsible for the unusual increase in the DW factor at high temperatures. In making a transition from an energy-minimum structure to another one, a small number of water molecules are involved in a highly cooperative fashion. The larger DW factor at higher temperature arises from jump-like motions of water molecules among these locally stable configurations which may or may not be a family of the proton-disordered ice forms satisfying the "ice rule".
Collapse
Affiliation(s)
- Hideki Tanaka
- Department of Chemistry, Faculty of Science, Okayama University, 3-1-1 Tsushima-naka, Okayama 700-8530, Japan
| | | |
Collapse
|
9
|
Doye JPK, Wales DJ. Saddle points and dynamics of Lennard-Jones clusters, solids, and supercooled liquids. J Chem Phys 2002. [DOI: 10.1063/1.1436470] [Citation(s) in RCA: 158] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023] Open
|
10
|
Chisolm ED, Clements BE, Wallace DC. Mean-atom-trajectory model for the velocity autocorrelation function of monatomic liquids. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2001; 63:031204. [PMID: 11308639 DOI: 10.1103/physreve.63.031204] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/04/2000] [Revised: 10/04/2000] [Indexed: 05/23/2023]
Abstract
We present a model for the motion of an average atom in a liquid or supercooled liquid state and apply it to calculations of the velocity autocorrelation function Z(t) and diffusion coefficient D. The model trajectory consists of oscillations at a distribution of frequencies characteristic of the normal modes of a single potential valley, interspersed with position- and velocity-conserving transits to similar adjacent valleys. The resulting predictions for Z(t) and D agree remarkably well with molecular dynamics simulations of Na at up to almost three times its melting temperature. Two independent processes in the model relax velocity autocorrelations: (a) dephasing due to the presence of many frequency components, which operates at all temperatures but which produces no diffusion, and (b) the transit process, which increases with increasing temperature and which produces diffusion. Because the model provides a single-atom trajectory in real space and time, including transits, it may be used to calculate all single-atom correlation functions.
Collapse
Affiliation(s)
- E D Chisolm
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | | | | |
Collapse
|
11
|
Vardeman CF, Gezelter JD. Comparing Models for Diffusion in Supercooled Liquids: The Eutectic Composition of the Ag−Cu Alloy. J Phys Chem A 2001. [DOI: 10.1021/jp0035784] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Charles F. Vardeman
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556
| | - J. Daniel Gezelter
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556
| |
Collapse
|
12
|
Schwartz SD. Quantum dynamics in condensed phases via extended modes and exact interaction propagator relations. J Chem Phys 2000. [DOI: 10.1063/1.1312280] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
13
|
|
14
|
Li WX, Keyes T. Instantaneous normal mode theory of diffusion and the potential energy landscape: Application to supercooled liquid CS2. J Chem Phys 1999. [DOI: 10.1063/1.479810] [Citation(s) in RCA: 29] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
15
|
Murry RL, Fourkas JT, Li WX, Keyes T. Molecular coordinates for instantaneous normal mode calculations. I. Coordinate dependence. J Chem Phys 1999. [DOI: 10.1063/1.479048] [Citation(s) in RCA: 23] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
16
|
Gezelter JD, Rabani E, Berne BJ. Calculating the hopping rate for diffusion in molecular liquids: CS2. J Chem Phys 1999. [DOI: 10.1063/1.478211] [Citation(s) in RCA: 23] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
17
|
|
18
|
|