1
|
Wang S, Ma Z, Pan W. Data-driven coarse-grained modeling of polymers in solution with structural and dynamic properties conserved. SOFT MATTER 2020; 16:8330-8344. [PMID: 32785383 DOI: 10.1039/d0sm01019g] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
We present data-driven coarse-grained (CG) modeling for polymers in solution, which conserves the dynamic as well as structural properties of the underlying atomistic system. The CG modeling is built upon the framework of the generalized Langevin equation (GLE). The key is to determine each term in the GLE by directly linking it to atomistic data. In particular, we propose a two-stage Gaussian process-based Bayesian optimization method to infer the non-Markovian memory kernel from the data of the velocity autocorrelation function (VACF). Considering that the long-time behaviors of the VACF and memory kernel for polymer solutions can exhibit hydrodynamic scaling (algebraic decay with time), we further develop an active learning method to determine the emergence of hydrodynamic scaling, which can accelerate the inference process of the memory kernel. The proposed methods do not rely on how the mean force or CG potential in the GLE is constructed. Thus, we also compare two methods for constructing the CG potential: a deep learning method and the iterative Boltzmann inversion method. With the memory kernel and CG potential determined, the GLE is mapped onto an extended Markovian process to circumvent the expensive cost of directly solving the GLE. The accuracy and computational efficiency of the proposed CG modeling are assessed in a model star-polymer solution system at three representative concentrations. By comparing with the reference atomistic simulation results, we demonstrate that the proposed CG modeling can robustly and accurately reproduce the dynamic and structural properties of polymers in solution.
Collapse
Affiliation(s)
- Shu Wang
- Department of Mechanical Engineering, University of Wisconsin-Madison, Madison, WI 53706, USA.
| | - Zhan Ma
- Department of Mechanical Engineering, University of Wisconsin-Madison, Madison, WI 53706, USA.
| | - Wenxiao Pan
- Department of Mechanical Engineering, University of Wisconsin-Madison, Madison, WI 53706, USA.
| |
Collapse
|
2
|
Guarini E, Neumann M, Bellissima S, Colognesi D, Bafile U. Density dependence of the dynamical processes governing the velocity autocorrelation function of a quantum fluid. Phys Rev E 2020; 100:062111. [PMID: 31962502 DOI: 10.1103/physreve.100.062111] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2019] [Indexed: 11/07/2022]
Abstract
We present an exponential mode analysis of the dynamical processes determining the time behavior of the Kubo velocity autocorrelation function (KVAF) of fluid para-H_{2}, as obtained by ring polymer molecular dynamics simulations at various fluid densities. The mechanisms contributing to the decay of the KVAF are thoroughly characterized at a slightly supercritical temperature, in a density interval ranging from the critical point to the fluid-solid transition. We show that the quantum nature of the system does not influence the specific phenomena and decay channels through which a loss in velocity correlation takes place, since these are the same as found in classical fluids. Similarly, a dynamical crossover is observed with increasing density, signaling the onset of a transverse-like dynamics like in classical systems. We also investigate the effect of density on the processes contributing to the most relevant property of a quantum fluid, namely, the large values of the total and zero-point kinetic energy arising through the Heisenberg uncertainty principle.
Collapse
Affiliation(s)
- Eleonora Guarini
- Dipartimento di Fisica e Astronomia, Università degli Studi di Firenze, via G. Sansone 1, I-50019 Sesto Fiorentino, Italy
| | - Martin Neumann
- Fakultät für Physik der Universität Wien, Strudlhofgasse 4, A-1090 Vienna, Austria
| | - Stefano Bellissima
- Consiglio Nazionale delle Ricerche, Istituto di Fisica Applicata "Nello Carrara", via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy
| | - Daniele Colognesi
- Consiglio Nazionale delle Ricerche, Istituto di Fisica Applicata "Nello Carrara", via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy
| | - Ubaldo Bafile
- Consiglio Nazionale delle Ricerche, Istituto di Fisica Applicata "Nello Carrara", via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy
| |
Collapse
|
3
|
Guarini E, Neumann M, Bafile U, Bellissima S, Colognesi D. Dynamical Origin of the Total and Zero-Point Kinetic Energy in a Quantum Fluid. PHYSICAL REVIEW LETTERS 2019; 123:135301. [PMID: 31697543 DOI: 10.1103/physrevlett.123.135301] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/25/2019] [Indexed: 06/10/2023]
Abstract
By applying an exponential mode analysis to ring polymer molecular dynamics simulations of dense fluid parahydrogen, we find that the dynamical processes establishing the time behavior of the Kubo velocity autocorrelation function display the same nature as those already observed in high-density classical fluids. This result permits us to demonstrate that the exponential mode decomposition is a unique tool to identify which dynamical processes lead to one of the most notable properties of quantum fluids: the large value of the mean kinetic energy per particle and the importance of the zero-temperature quantum effects in determining it.
Collapse
Affiliation(s)
- Eleonora Guarini
- Dipartimento di Fisica e Astronomia, Università degli Studi di Firenze, via G. Sansone 1, I-50019 Sesto Fiorentino, Italy
| | - Martin Neumann
- Fakultät für Physik der Universität Wien, Strudlhofgasse 4, A-1090 Wien, Austria
| | - Ubaldo Bafile
- Consiglio Nazionale delle Ricerche, Istituto di Fisica Applicata "Nello Carrara," via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy
| | - Stefano Bellissima
- Consiglio Nazionale delle Ricerche, Istituto di Fisica Applicata "Nello Carrara," via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy
| | - Daniele Colognesi
- Consiglio Nazionale delle Ricerche, Istituto di Fisica Applicata "Nello Carrara," via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy
| |
Collapse
|
4
|
Mizuta K, Ishii Y, Kim K, Matubayasi N. Bridging the gap between molecular dynamics and hydrodynamics in nanoscale Brownian motions. SOFT MATTER 2019; 15:4380-4390. [PMID: 31086871 DOI: 10.1039/c9sm00246d] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Through molecular dynamics simulations, we examined the hydrodynamic behavior of the Brownian motion of fullerene particles based on molecular interactions. The solvation free energy and velocity autocorrelation function (VACF) were calculated by using the Lennard-Jones (LJ) and Weeks-Chandler-Andersen (WCA) potentials for the solute-solvent and solvent-solvent interactions and by changing the size of the fullerene particles. We also measured the diffusion constant of the fullerene particles and the shear viscosity of the host fluid, and then the hydrodynamic radius aHD was quantified from the Stokes-Einstein relation. The aHD value exceeds that of the gyration radius of the fullerene when the solvation free energy exhibits largely negative values using the LJ potential. In contrast, aHD is similar to the size of bare fullerene when the solvation free energy is positive using the WCA potential. Furthermore, the VACF of the fullerene particles is directly comparable with the analytical expressions utilizing the Navier-Stokes equations both in incompressible and compressible forms. A hydrodynamic long-time tail t-3/2 is demonstrated for timescales longer than the kinematic time of the momentum diffusion over the particle size. However, the VACF at shorter timescales deviates from the hydrodynamic description, particularly for smaller fullerene particles and for the LJ potential. This occurs even though the compressible effect is considered when characterizing the decay of the VACF around the sound-propagation timescale over the particle size. These results indicate that the nanoscale Brownian motion is influenced by the solvation structure around the solute particles originating from the molecular interaction.
Collapse
Affiliation(s)
- Keisuke Mizuta
- Division of Chemical Engineering, Department of Materials Engineering Science, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan.
| | | | | | | |
Collapse
|
5
|
Heyes DM, Smith ER, Dini D. Shear stress relaxation and diffusion in simple liquids by molecular dynamics simulations: Analytic expressions and paths to viscosity. J Chem Phys 2019; 150:174504. [DOI: 10.1063/1.5095501] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- D. M. Heyes
- Department of Mechanical Engineering, Imperial College London, Exhibition Road, South Kensington, London SW7 2AZ, United Kingdom
| | - E. R. Smith
- Department of Mechanical Engineering, Imperial College London, Exhibition Road, South Kensington, London SW7 2AZ, United Kingdom
- Department of Mechanical and Aerospace Engineering, Brunel University London, Uxbridge, Middlesex UB8 3PH, United Kingdom
| | - D. Dini
- Department of Mechanical Engineering, Imperial College London, Exhibition Road, South Kensington, London SW7 2AZ, United Kingdom
| |
Collapse
|
6
|
Bellissima S, Neumann M, Bafile U, Colognesi D, Barocchi F, Guarini E. Density and time scaling effects on the velocity autocorrelation function of quantum and classical dense fluid para-hydrogen. J Chem Phys 2019; 150:074502. [DOI: 10.1063/1.5085202] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Affiliation(s)
- S. Bellissima
- Consiglio Nazionale delle Ricerche, Istituto di Fisica Applicata “Nello Carrara,” via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy
| | - M. Neumann
- Fakultät für Physik der Universität Wien, Strudlhofgasse 4, A-1090 Wien, Austria
| | - U. Bafile
- Consiglio Nazionale delle Ricerche, Istituto di Fisica Applicata “Nello Carrara,” via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy
| | - D. Colognesi
- Consiglio Nazionale delle Ricerche, Istituto di Fisica Applicata “Nello Carrara,” via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy
| | - F. Barocchi
- Dipartimento di Fisica e Astronomia, Università degli Studi di Firenze, via G. Sansone 1, I-50019 Sesto Fiorentino, Italy
| | - E. Guarini
- Dipartimento di Fisica e Astronomia, Università degli Studi di Firenze, via G. Sansone 1, I-50019 Sesto Fiorentino, Italy
| |
Collapse
|
7
|
Torres-Carbajal A, Castañeda-Priego R. Friction and diffusion of a nano-colloidal disk in a two-dimensional solvent with a liquid-liquid transition. Phys Chem Chem Phys 2018; 20:6917-6928. [PMID: 29464245 DOI: 10.1039/c7cp08302e] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We report on the friction and diffusion of a single mobile nano-colloidal disk, whose size and mass are one and two orders of magnitude, respectively, greater than the molecules of the host solvent; all particles are restricted to move in a two-dimensional space. Using molecular dynamics simulations, the variation of the transport coefficients as a function of the thermodynamic state of the supporting fluid, in particular, around those states in the neighbourhood of the liquid-liquid phase coexistence, is investigated. The diffusion coefficient is determined through the fit of the mean-square displacement at long times and with the Green-Kubo relationship for the velocity autocorrelation function, whereas the friction coefficient is computed from the correlation of the fluctuating force. From the determination of the transport properties, the applicability of the Stokes-Einstein relation in two dimensions around the second critical point is discussed.
Collapse
Affiliation(s)
- Alexis Torres-Carbajal
- División de Ciencias e Ingenierías, Campus León, Universidad de Guanajuato, Loma del Bosque 103, Lomas del Campestre, 37150 León, Guanajuato, Mexico.
| | | |
Collapse
|
8
|
Han KH, Kim C, Talkner P, Karniadakis GE, Lee EK. Molecular hydrodynamics: Vortex formation and sound wave propagation. J Chem Phys 2018; 148:024506. [PMID: 29331127 DOI: 10.1063/1.5011992] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Abstract
In the present study, quantitative feasibility tests of the hydrodynamic description of a two-dimensional fluid at the molecular level are performed, both with respect to length and time scales. Using high-resolution fluid velocity data obtained from extensive molecular dynamics simulations, we computed the transverse and longitudinal components of the velocity field by the Helmholtz decomposition and compared them with those obtained from the linearized Navier-Stokes (LNS) equations with time-dependent transport coefficients. By investigating the vortex dynamics and the sound wave propagation in terms of these field components, we confirm the validity of the LNS description for times comparable to or larger than several mean collision times. The LNS description still reproduces the transverse velocity field accurately at smaller times, but it fails to predict characteristic patterns of molecular origin visible in the longitudinal velocity field. Based on these observations, we validate the main assumptions of the mode-coupling approach. The assumption that the velocity autocorrelation function can be expressed in terms of the fluid velocity field and the tagged particle distribution is found to be remarkably accurate even for times comparable to or smaller than the mean collision time. This suggests that the hydrodynamic-mode description remains valid down to the molecular scale.
Collapse
Affiliation(s)
- Kyeong Hwan Han
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, South Korea
| | - Changho Kim
- Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Peter Talkner
- Institut für Physik, Universität Augsburg, 86159 Augsburg, Germany
| | - George Em Karniadakis
- Division of Applied Mathematics, Brown University, Providence, Rhode Island 02912, USA
| | - Eok Kyun Lee
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, South Korea
| |
Collapse
|
9
|
Bellissima S, Neumann M, Guarini E, Bafile U, Barocchi F. Density of states and dynamical crossover in a dense fluid revealed by exponential mode analysis of the velocity autocorrelation function. Phys Rev E 2017; 95:012108. [PMID: 28208443 DOI: 10.1103/physreve.95.012108] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2016] [Indexed: 11/07/2022]
Abstract
Extending a preceding study of the velocity autocorrelation function (VAF) in a simulated Lennard-Jones fluid [Phys. Rev. E 92, 042166 (2015)PLEEE81539-375510.1103/PhysRevE.92.042166] to cover higher-density and lower-temperature states, we show that the recently demonstrated multiexponential expansion method allows for a full account and understanding of the basic dynamical processes encompassed by a fundamental quantity as the VAF. In particular, besides obtaining evidence of a persisting long-time tail, we assign specific and unambiguous physical meanings to groups of exponential modes related to the longitudinal and transverse collective dynamics, respectively. We have made this possible by consistently introducing the interpretation of the VAF frequency spectrum as a global density of states in fluids, generalizing a solid-state concept, and by giving to specific spectral components, obtained through the VAF exponential expansion, the corresponding meaning of partial densities of states relative to specific dynamical processes. The clear identification of a high-frequency oscillation of the VAF with the near-top excitation frequency in the dispersion curve of acoustic waves is a neat example of the power of the method. As for the transverse mode contribution, its analysis turns out to be particularly important, because the multiexponential expansion reveals a transition marking the onset of propagating excitations when the density is increased beyond a threshold value. While this finding agrees with the recent literature debating the issue of dynamical crossover boundaries, such as the one identified with the Frenkel line, we can add detailed information on the modes involved in this specific process in the domains of both time and frequency. This will help obtain a still missing full account of transverse dynamics, in both its nonpropagating and propagating aspects which are linked through dynamical transitions depending on both the thermodynamic states and the excitation wave vectors.
Collapse
Affiliation(s)
- S Bellissima
- Dipartimento di Fisica e Astronomia, Università degli Studi di Firenze, via G. Sansone 1, I-50019 Sesto Fiorentino, Italy
| | - M Neumann
- Fakultät für Physik der Universität Wien, Strudlhofgasse 4, A-1090 Wien, Austria
| | - E Guarini
- Dipartimento di Fisica e Astronomia, Università degli Studi di Firenze, via G. Sansone 1, I-50019 Sesto Fiorentino, Italy
| | - U Bafile
- Consiglio Nazionale delle Ricerche, Istituto dei Sistemi Complessi, via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy
| | - F Barocchi
- Dipartimento di Fisica e Astronomia, Università degli Studi di Firenze, via G. Sansone 1, I-50019 Sesto Fiorentino, Italy
| |
Collapse
|
10
|
Kondratyuk ND, Norman GE, Stegailov VV. Self-consistent molecular dynamics calculation of diffusion in higher n-alkanes. J Chem Phys 2016; 145:204504. [PMID: 27908129 DOI: 10.1063/1.4967873] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
Diffusion is one of the key subjects of molecular modeling and simulation studies. However, there is an unresolved lack of consistency between Einstein-Smoluchowski (E-S) and Green-Kubo (G-K) methods for diffusion coefficient calculations in systems of complex molecules. In this paper, we analyze this problem for the case of liquid n-triacontane. The non-conventional long-time tails of the velocity autocorrelation function (VACF) are found for this system. Temperature dependence of the VACF tail decay exponent is defined. The proper inclusion of the long-time tail contributions to the diffusion coefficient calculation results in the consistency between G-K and E-S methods. Having considered the major factors influencing the precision of the diffusion rate calculations in comparison with experimental data (system size effects and force field parameters), we point to hydrogen nuclear quantum effects as, presumably, the last obstacle to fully consistent n-alkane description.
Collapse
Affiliation(s)
- Nikolay D Kondratyuk
- Joint Institute for High Temperatures of the Russian Academy of Sciences, Moscow 125412, Russia
| | - Genri E Norman
- Joint Institute for High Temperatures of the Russian Academy of Sciences, Moscow 125412, Russia
| | - Vladimir V Stegailov
- Joint Institute for High Temperatures of the Russian Academy of Sciences, Moscow 125412, Russia
| |
Collapse
|
11
|
Rees RJ, Spencer MJ. The science and life of Ian K. Snook. MOLECULAR SIMULATION 2015. [DOI: 10.1080/08927022.2016.1091560] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
12
|
Bellissima S, Neumann M, Guarini E, Bafile U, Barocchi F. Time dependence of the velocity autocorrelation function of a fluid: An eigenmode analysis of dynamical processes. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 92:042166. [PMID: 26565227 DOI: 10.1103/physreve.92.042166] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/06/2015] [Indexed: 06/05/2023]
Abstract
The velocity autocorrelation function (VAF), a key quantity in the atomic-scale dynamics of fluids, has been the first paradigmatic example of a long-time tail phenomenon, and much work has been devoted to detecting such long-lasting correlations and understanding their nature. There is, however, much more to the VAF than simply the evidence of this long-time dynamics. A unified description of the VAF from very short to long times, and of the way it changes with varying density, is still missing. Here we show that an approach based on very general principles makes such a study possible and opens the way to a detailed quantitative characterization of the dynamical processes involved at all time scales. From the analysis of molecular dynamics simulations for a slightly supercritical Lennard-Jones fluid at various densities, we are able to evidence the presence of distinct fast and slow decay channels for the velocity correlation on the time scale set by the collision rate. The density evolution of these decay processes is also highlighted. The method presented here is very general, and its application to the VAF can be considered as an important example.
Collapse
Affiliation(s)
- S Bellissima
- Dipartimento di Fisica e Astronomia, Università degli Studi di Firenze, via G. Sansone 1, I-50019 Sesto Fiorentino, Italy
| | - M Neumann
- Fakultät für Physik der Universität Wien, Strudlhofgasse 4, A-1090 Wien, Austria
| | - E Guarini
- Dipartimento di Fisica e Astronomia, Università degli Studi di Firenze, via G. Sansone 1, I-50019 Sesto Fiorentino, Italy
| | - U Bafile
- Consiglio Nazionale delle Ricerche, Istituto dei Sistemi Complessi, via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy
| | - F Barocchi
- Dipartimento di Fisica e Astronomia, Università degli Studi di Firenze, via G. Sansone 1, I-50019 Sesto Fiorentino, Italy
| |
Collapse
|
13
|
Ryltsev RE, Chtchelkatchev NM. Hydrodynamic anomalies in supercritical fluid. J Chem Phys 2014; 141:124509. [DOI: 10.1063/1.4895726] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- R. E. Ryltsev
- Institute of Metallurgy, Ural Division of Russian Academy of Sciences, 620016 Yekaterinburg, Russia
- L.D. Landau Institute for Theoretical Physics, Russian Academy of Sciences, 142432 Chernogolovka, Moscow Region, Russia
| | - N. M. Chtchelkatchev
- Moscow Institute of Physics and Technology, 141700 Moscow, Russia
- Institute for High Pressure Physics, Russian Academy of Sciences, 142190 Troitsk, Russia
- L.D. Landau Institute for Theoretical Physics, Russian Academy of Sciences, 142432 Chernogolovka, Moscow Region, Russia
- Department of Physics and Astronomy, California State University Northridge, Northridge, California 91330, USA
| |
Collapse
|
14
|
Patra TK, Hens A, Singh JK. Vapor-liquid phase coexistence and transport properties of two-dimensional oligomers. J Chem Phys 2012; 137:084701. [DOI: 10.1063/1.4747195] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
15
|
Morishita T, Spencer MJ, Russo SP, Snook IK, Mikami M. Surface reconstruction of ultrathin silicon nanosheets. Chem Phys Lett 2011. [DOI: 10.1016/j.cplett.2011.03.004] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
16
|
Lee SH. Friction between Two Brownian Particles in a Lennard-Jones Solvent: A Molecular Dynamics Simulation Study. B KOREAN CHEM SOC 2010. [DOI: 10.5012/bkcs.2010.31.8.2402] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
17
|
Donkó Z, Goree J, Hartmann P, Liu B. Time-correlation functions and transport coefficients of two-dimensional Yukawa liquids. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2009; 79:026401. [PMID: 19391849 DOI: 10.1103/physreve.79.026401] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/10/2007] [Revised: 10/10/2008] [Indexed: 05/27/2023]
Abstract
The existence of coefficients for diffusion, viscosity, and thermal conductivity is examined for two-dimensional (2D) liquids. Equilibrium molecular dynamics simulations are performed using a Yukawa potential and the long-time behavior of autocorrelation functions is tested. Advances reported here as compared to previous 2D Yukawa liquid simulations include an assessment of the thermal conductivity, using a larger system size to allow meaningful examination of longer times, and development of improved analysis methods. We find that the transport coefficient exists for diffusion at high temperature and viscosity at low temperature, but not in the opposite limits. The thermal conductivity coefficient does not appear to exist at high temperature. Further advances in computing power could improve these assessments by allowing even larger system sizes and longer time series.
Collapse
Affiliation(s)
- Z Donkó
- Research Institute for Solid State Physics and Optics of the Hungarian Academy of Sciences, H-1525 Budapest, P.O. Box 49, Hungary
| | | | | | | |
Collapse
|
18
|
Ott T, Bonitz M, Donkó Z, Hartmann P. Superdiffusion in quasi-two-dimensional Yukawa liquids. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2008; 78:026409. [PMID: 18850948 DOI: 10.1103/physreve.78.026409] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/04/2008] [Indexed: 05/26/2023]
Abstract
The emergence and vanishing of superdiffusion in quasi-two-dimensional Yukawa systems are investigated by molecular dynamics simulations. Using both the asymptotic behavior of the mean-squared displacement of the particles and the long-time tail of the velocity autocorrelation function as indicators of superdiffusion, we confirm the existence of a transition from normal diffusion to superdiffusion in systems changing from a three-dimensional to a two-dimensional character. A connection between superdiffusion and dimensionality is established by the behavior of the projected pair distribution function.
Collapse
Affiliation(s)
- T Ott
- Christian-Albrechts-Universität zu Kiel, Institut für Theoretische Physik und Astrophysik, Leibnizstrasse 15, 24098 Kiel, Germany
| | | | | | | |
Collapse
|
19
|
Ould-Kaddour F, Levesque D. Diffusion of nanoparticles in dense fluids. J Chem Phys 2007; 127:154514. [DOI: 10.1063/1.2794753] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
20
|
Cappelezzo M, Capellari CA, Pezzin SH, Coelho LAF. Stokes-Einstein relation for pure simple fluids. J Chem Phys 2007; 126:224516. [PMID: 17581072 DOI: 10.1063/1.2738063] [Citation(s) in RCA: 63] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The authors employed the equilibrium molecular dynamics technique to calculate the self-diffusion coefficient and the shear viscosity for simple fluids that obey the Lennard-Jones 6-12 potential in order to investigate the validity of the Stokes-Einstein (SE) relation for pure simple fluids. They performed calculations in a broad range of density and temperature in order to test the SE relation. The main goal of this work is to exactly calculate the constant, here denominated by alpha, present in the SE relation. Also, a modified SE relation where a fluid density is raised to a power in the usual expression is compared to the classical expression. According to the authors' simulations slip boundary conditions (alpha=4) can be satisfied in some state points. An intermediate value of alpha=5 was found in some regions of the phase diagram confirming the mode coupling theory. In addition depending on the phase diagram point and the definition of hydrodynamics radius, stick boundary condition (alpha=6) can be reproduced. The authors investigated the role of the hydrodynamic radius in the SE relation using three different definitions. The authors also present calculations for alpha in a hard-sphere system showing that the slip boundary conditions hold at very high density. They discuss possible explanations for their results and the role of the hydrodynamic radius for different definitions in the SE relation.
Collapse
Affiliation(s)
- M Cappelezzo
- Department of Mathematics, State University of Santa Catarina, Campus Universitário Prof. Avelino Marcante s/n, Joinville 89223-100, Santa Catarina, Brazil
| | | | | | | |
Collapse
|
21
|
Williams SR, McGlynn P, Bryant G, Snook IK, van Megen W. Dynamical signatures of freezing: stable fluids, metastable fluids, and crystals. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2006; 74:031204. [PMID: 17025615 DOI: 10.1103/physreve.74.031204] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/01/2006] [Revised: 07/17/2006] [Indexed: 05/12/2023]
Abstract
Mean squared displacements and velocity auto correlation functions are calculated using molecular dynamics for hard spheres under a range of conditions (i) for the equilibrium fluid below freezing; (ii) for the metastable fluid above freezing; and (iii) for the hard sphere crystal, both in the metastable region between freezing and melting, and in the stable region above melting. In addition, simulations are carried out for a metastable Lennard-Jones system. The results confirm recent studies that indicated the disappearance of the classical Alder long-time tail, and show that they apply to systems other than the metastable hard sphere fluid. The implications of these results for our understanding of crystallization and the glass transition are discussed.
Collapse
Affiliation(s)
- Stephen R Williams
- Research School of Chemistry, Australian National University, Canberra, ACT 0200, Australia.
| | | | | | | | | |
Collapse
|
22
|
Dib RFA, Ould-Kaddour F, Levesque D. Long-time behavior of the velocity autocorrelation function at low densities and near the critical point of simple fluids. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2006; 74:011202. [PMID: 16907082 DOI: 10.1103/physreve.74.011202] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/15/2006] [Indexed: 05/11/2023]
Abstract
Numerous theoretical and numerical works have been devoted to the study of the algebraic decrease at large times of the velocity autocorrelation function of particles in a fluid. The derivation of this behavior, the so-called long-time tail, generally based on linearized hydrodynamics, makes no reference to any specific characteristic of the particle interactions. However, in the literature doubts have been expressed about the possibility that by numerical simulations the long-time tail can be observed in the whole fluid phase domain of systems in which the particles interact by soft-core and attractive pair potentials. In this work, extensive and accurate molecular-dynamics simulations establish that the predicted long-time tail of the velocity autocorrelation function exists in a low-density fluid of particles interacting by a soft-repulsive potential and near the liquid-gas critical point of a Lennard-Jones system. These results contribute to the confirmation that the algebraic decay of the velocity autocorrelation function is universal in these fluid systems.
Collapse
Affiliation(s)
- R F A Dib
- Laboratoire de Physique Théorique, Faculté des Sciences, Université de Tlemcen Boîte Postale 119, Tlemcen 13000, Algeria
| | | | | |
Collapse
|
23
|
Meier K, Laesecke A, Kabelac S. Transport coefficients of the Lennard-Jones model fluid. II Self-diffusion. J Chem Phys 2004; 121:9526-35. [PMID: 15538874 DOI: 10.1063/1.1786579] [Citation(s) in RCA: 94] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In an extensive computer simulation study, the transport coefficients of the Lennard-Jones model fluid were determined with high accuracy from equilibrium molecular-dynamics simulations. In the frame of time-correlation function theory, the generalized Einstein relations were employed to evaluate the transport coefficients. This second of a series of four papers presents the results for the self-diffusion coefficient, and discusses and interprets the behavior of this transport coefficient in the fluid region of the phase diagram. The uncertainty of the self-diffusion data is estimated to be 1% in the gas region and 0.5% at high-density liquid states. With the very accurate data, even fine details in the shape of the self-diffusion isotherms are resolved, and the previously little-investigated behavior of the self-diffusion coefficient at low-density gaseous states is analyzed in detail. Finally, aspects of the mass transport mechanisms on the molecular scale are explored by an analysis of the velocity autocorrelation functions.
Collapse
Affiliation(s)
- Karsten Meier
- Institut für Thermodynamik, Helmut-Schmidt-Universität--Universität der Bundeswehr Hamburg, Holstenhofweg 85, D-22043 Hamburg, Germany.
| | | | | |
Collapse
|
24
|
Schmidt JR, Skinner JL. Hydrodynamic boundary conditions, the Stokes–Einstein law, and long-time tails in the Brownian limit. J Chem Phys 2003. [DOI: 10.1063/1.1610442] [Citation(s) in RCA: 107] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
25
|
Snook I, O'Malley B, McPhie M, Daivis P. The approach to the Brownian limit in particulate dispersions. J Mol Liq 2003. [DOI: 10.1016/s0167-7322(02)00158-7] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
26
|
Barnard A, Russo S, Leach G. Nearest neighbour considerations in Stillinger-Weber type potentials for diamond. MOLECULAR SIMULATION 2002. [DOI: 10.1080/0892702021000002476] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|