arxiv: v1 [physics.chem-ph] 4 Feb 2014

Size: px

Start display at page:

Download "arxiv: v1 [physics.chem-ph] 4 Feb 2014"

Ada Russell
5 years ago
Views:

1 Discrete dynamics versus analytic dynamics Søren Toxvaerd DNRF centre Glass and Time, IMFUFA, Department of Sciences, Roskilde University, arxiv: v1 [physics.chem-ph] 4 Feb 2014 Postbox 260, DK-4000 Roskilde, Denmark (Dated: September 19, 2018) Abstract For discrete classical Molecular dynamics (MD) obtained by the Verlet algorithm (VA) with the time increment h there exists a shadow Hamiltonian H with energy Ẽ(h), for which the discrete particle positions lie on the analytic trajectories for H. Here we proof that there, independent of such an analytic analogy, exists an exact hidden energy invariance E for VA dynamics. The fact that the discrete VA dynamics has the same invariances as Newtonian dynamics raises the question, which of the formulations that are correct, or alternatively, the most appropriate formulation of classical dynamics. In this context the relation between the discrete VA dynamics and the (general) discrete dynamics investigated by T. D. Lee [Phys. Lett. 122B, 217 (1983)] is presented and discussed. 1

2 I. INTRODUCTION Molecular Dynamics (MD) generates the time evolution of N classical mechanical particles by discrete time propagation. Almost all the MD are obtained by the Verlet algorithm (VA) [1] where a new position r i (t+δt) of the i th particle with mass m i at time t+δt is obtained from the force f i (t) and the two last discrete positions r i (t+δt) = 2r i (t) r i (t δt)+ δt2 m i f i (t). (1) The algorithm is the central-difference expression for the mass times the acceleration of the particle which equals the force f i, and it appears in the literature with different names (Verlet, leap-frog, velocity Verlet,..) [2]. The algorithm is time reversible and symplectic, and the different reformulations of the algorithm do not change the discrete time evolution and the physics obtained by the VA dynamics. Mathematical investigations [3 5] have proved the existence of a shadow Hamiltonian H [6] for symplectic algorithms. The proof is obtained by an asymptotic expansion, but the series for the shadow Hamiltonian does not converge in the general case. For a review of the asymptotic expansion, its convergence and optimal truncation see [7]. Only the harmonic approximation, E(1) of the first term in this expansion is known explicitly [2, 6]. But inclusion of E(1) in the traditional obtained zero order energy for MD systems with Lennard-Jones (LJ) particles reduces the fluctuation in the energy by a factor of hundred for traditional values of δt [2] and makes it possible to obtain the shadow energy, Ẽ of the analytic dynamics with high precision. The VA algorithm deviates, however, from all other algorithms for classical dynamics by that momenta are not dynamical variables. Furthermore, the discrete VA dynamics for a harmonic oscillator (DDHO), which can be solved exactly, reveals that the DDHO not only has an asymptotic expansion with an underlying analytic shadow Hamiltonian. But the DDHO dynamics also has a (hidden) energy invariance,e which, independent of an existence of an analytic shadow Hamiltonian, is conserved step by step during the discrete time evolution. Below we show that this hidden invariance is a general quality of the discrete VA dynamics, independent of the existence of a shadow Hamiltonian, and that the discrete VA dynamics has the same qualities and conserved invariances as analytic Newtonian dynamics. In order to prove the existence of a hidden energy invariance we must not make use of any analytic tools. This might seem to be a hopeless agenda, but on the other hand 2

3 the exact solution for a discrete harmonic oscillator [6] makes no use of analyticity and the exact solution has an energy invariance E which, in the analytic limit is equal to the energy of analytic Newtonian dynamics. II. THE HIDDEN ENERGY INVARIANCE The kinetic energy in analytic dynamics is obtained from the momenta. The positions in Eq. (1) are the only dynamic variables in the discrete VA dynamics, i.e., the momenta, p i, are not. Consequently, an expression for the total momentum of the system requires a choice ofanexpression forthemomentum p i ofthei thparticleintermsofitspositions. The sentences momenta, energy, potential energy, kinetic energy, and work should be given by quotations in discrete VA dynamics to underline the fact, that the N objects in the discrete dynamics only exercise mutual irritations, or forces f i (t n ) at their positions at the discrete time t n = nδt. With the definition of the momenta p i (t n 1,t n ) m i r i (t n ) r i (t n 1 ) δt (2) it follows immediately from the algorithm that the total momentum and angular momentum are conserved for conservative systems with N i f i (t) = 0 [8]. But the momenta and thereby the kinetic energies appear asynchronous with the discrete positions and they are not a function of a single set of the discrete positions. The proof of an invariance, equivalent to the conserved energy in the analytic dynamics is more difficult, but it can be obtained by proving that there exists a hidden energy invariance, E, of the N objects dynamics with the change δe n = 0 by the discrete step which brings the N positions R n with the forces F n at time t n to R n+1 at t n+1. Forsimplicity considern particleswithequalmassesm i = m, andwiththemassincluded in the discrete time increment, i.e. with h δt m. A step with discrete dynamics changes the kineticenergy ofthesystembyδk anditsability, δu, toperforma work, δu = W. Since the momenta and thereby the kinetic energy is given by two sets of positions, a change in kinetic energy is given by three consecutive sets of positions. The proof is obtained by consider two consecutive time steps. A new set of positions, R n+1 is obtained at the n th time step from the two previous sets, R n 1, R n and the forces F n, by which the change in 3

4 the kinetic energy can be defined as δkn 1 ( ) 2 Rn+1 R n 1 ( ) 2 Rn R n 1. (3) 2 h 2 h The definition of the change in the kinetic energy for discrete VA dynamics is consistent with the definition of the momenta, Eq. (2). The forces F n bring the N particles to the positions R n+1 and with a change in the ability, δu n, to perform a work, δu n = W n. For the two steps we define the total change and the discrete dynamics obeys the relation 2δU n = 2W n F n (R n+1 R n 1 ), (4) δu n +δk n = 0. (5) The proof starts by noticing that if one instead of the (NVE) dynamics, obtained by Eq. (1) with a constant time increment h, adjust the (n+1) th time increment h n so W = 0, one obtains a geodesic step (NVU) [9] to the positions R n+1 (U) which differs from R n+1. If is inserted in the Verlet algorithm; Eq. (1) 2W(U) n = F n (R n+1 (U) R n 1 ) = 0 (6) R n+1 (U) = 2R n R n 1 +h 2 nf n, (7) one obtains an expression for h 2 n at the NVU step at time t n [9] h 2 n = 2F n (R n R n 1 ) F 2. (8) n I.e. instead of propagating the system the n th step with the constant time increment h, the increment h n is adjusted to ensure that the system ability to perform a work is unchanged. The NVU step at time t n updates the position to R n+1 (U) and with the geodesic invariance: the constant length of the steps [9] (R n+1 (U) R n ) 2 = (R n R n 1 ) 2, (9) which is obtained by rearranging and squaring Eq. (7) ( (R n+1 (U) R n ) 2 = R n R n 1 2 F ) 2 n (R n R n 1 ) F 2 F n n (10) = (R n R n 1 ) 2, (11) 4

5 i.e. with the change in kinetic energy δk(u) n = 1 ( ) 2 Rn+1 (U) R n 1 ( ) 2 Rn R n 1 = 0. (12) 2 h 2 h So the NVU step to R n+1 (U) obeys δu(u) n = δk(u) n = 0. (13) We are now able to proof the existence of an energy invariance (Eq. (5)) by the VA dynamics, Eq. (1). The proof can e.g. be obtained by deriving the difference between the NVU and the NVE step at t n. The new positions R n+1 and R n+1 (U) are both obtained from R n 1,R n and F n, but with different time increments. With NVE 2δU n = 2W n = F n (R n+1 R n 1 ) = F n (R n+1 R n+1 (U)+R n+1 (U) R n 1 ) = F n (R n+1 R n+1 (U)). (14) The difference R n+1 R n+1 (U) can be obtained from the Verlet algorithm, Eq. (1) and the NVU algorithm, Eqn. (7) and (8), and gives 2δU n = 2W n = (h2 F 2 n +2F n (R n R n 1 )). (15) The change in the kinetic energy is obtained from the NVE algorithm Eq. (1) ( ) 2 ( ) 2 2δKn = Rn+1 R n Rn R n 1 h h = (R n R n 1 ) 2 +h 4 F 2 ( n +2h 2 F n (R n R n 1 ) Rn R n 1 h 2 h = h 2 F 2 n +2F n (R n R n 1 ), (16) and Eq. (5) for the discrete dynamics is obtained from Eqn. (15) and (16). The change in kinetic energy in discrete dynamics must necessarily be obtained from two consecutive steps and the change in the systems ability to perform a work is consistently obtained from the same sets of positions. But, by eliminating R n 1 in Eq. (4), one obtains an expression for the change in the ability per time step, δu n δu n = 1 2 F n (R n+1 R n 1 ) = F n (R n+1 R n )+ h2 2 F2 n. (17) 5 ) 2

6 In the Newtonian dynamics the existence of a potential energy state function, U(R), is ensured by that the total work done around any closed circuit from R is zero [10]. The VA dynamics is started from two sets of positions R 0 and R 1, and the time increment h. The ability to perform a discrete work U (R) is also a state function and it plays the same role as the potential energy U(R) for analytic Newtonian dynamics. Consider any discrete closed sequence of positions generated with VA dynamics with the time increment h and which starts and ends with the same two configuration R 0 and R 1. The total change in the kinetic energy is i=n δki = (R 2 R 1 ) 2 h i=1 (R 1 R 0) 2 h + (R 3 R 2 ) 2 h (R 2 R 1) 2 h (R 1 R 0) 2 h The start ability is U (R 1 ), and since all the terms in i=n i=1 (δu i +δk i) are zero accordingly to Eq. (5) it implies that The energy invariance = 0. (18) i=n δui = 0. (19) i=1 E = U n (R n)+k n (R n 1,R n ) (20) is given by the start condition for the discrete dynamics and it differs from the energy invariance of Newtonian dynamics. It is a state function, and due to the discrete dynamics itdependsontwoconsecutivesetsofthepositions, R n 1,R n, insteadoftheenergyinvariance in Newtonian dynamics which depends on the positions R(t n ) and the momenta P(t n ) at the same time t n. The two invariances are, however, equal in the analytic limit [10] lim h 0 (δu n +δk n ) = F(t n) δr(t n )+δk(t n )+O(h 2 ) = 0, (21) where the term O(h 2 ) = h2 2 F(t n) 2 is the total deviation from the Newtonian dynamics. The invariance, Eq. (20) does not depend on a convergence of an asymptotic expansion, and it differs also from the shadow energy for the shadow Hamiltonian by that, although the change contains two terms, the expressions for their changes do not make use of a potential, but only of the forces and the discrete positions. It is obtained by noticing that, with a suitable definition of the work and kinetic energy, δun = δk n, and by formulating the requirement that U n(r n ) is a state function. The derivation is a copy of the derivation of the energy invariance for Newtonian dynamics [10]. In Thermodynamics the First law of 6

7 thermodynamics is formulated exactly in the same manner, but as a basic assumption of that the energy function is a state function consisting of two terms which change by work and kinetic energy exchanges, and the present formulations is the corresponding formulation of the energy conservation in dynamics and thermodynamics for discrete VA dynamics. A. Energy conservation in MD with VA dynamics The formulation of energy in discrete VA dynamics by the ability U (R) to perform a discrete work instead of the potential energy U(R) works equally well as the traditional formulation. Molecular Dynamics simulations with VA for N particles are obtained from two consecutive start sets of positions, R 0 and R 1, and these positions define not only the total dynamics evolution, but also the mean value of traditional zero order energy, < E(0) n >, the accurate first order estimate of the shadow energy, Ẽn E(0) n +h 2 E(1) n of theunderlying analyticdynamicsandtheexact energyinvariancee. SincethechangeinE is given in the same manner as the energy conservation by the First law of thermodynamics, we need to define a start ability, U 1 for the discrete dynamics. But the discrete VA dynamics differs in fact neither from the analytic counterpart at this point. In principle we could obtain the ability U 1 at the start of the simulation by determining the discrete work performed by bringing the particles from infinite separations via R 0 to the positions R 1. In the thermodynamics one defines, however, a standard state of energy (enthalpy), and here we will use the potential energy U 1 (R 1 ) at the positions R 1 and the accurate estimate, Ẽ 1 at the start of the dynamics, and obtain and the energy invariance with U 1 = U 1(R 1 ), (22) i=n E (t n+1 ) = Ẽ1 + [δun +δkn] (23) i=1 δu n = F i (R i+1 R i )+ h2 2 F2 i. (24) The energy evolution by MD in double precision arithmetic with VA was determined for two systems. In the first a liquid system of N = 2000 LJ particles at the density ρ = 0.80σ 3 was calibrated at the temperature kt/ǫ = 1. The thermostat [8] was switched off and the 7

8 energy evolution in the next ten thousand time steps with h = was obtained. Figure Energi Number of steps FIG. 1. Discrete energies of a LJ system with the Verlet algorithm at T,ρ = 1.0, 0.80 and for h= Blue solid line: Traditional energy estimate E(0) n ; green dashed line: Shadow energy Ẽ n E(0) n +h 2 E(1) n ; red dash-dotted line: the energy invariance E. The inset shows Ẽn and E. Excess ability per particle and excess temperature (U * n-u * 1)/N (T * n-t * 1)/N Number of steps FIG. 2. Excess discrete energies end temperatures of a non-equilibrium system of N = 2048 LJ particles at spontaneous melting. Red line: the excess ability per particle (Un U 1 )/N; green long dashed: excess temperature T n T 1 ; blue short dashes: the constant energy invariance E n = E 1. 1 shows the energy evolution E(0) n,ẽn and En for the first hundred time steps. The first order estimate of Ẽn (green dashes) improves the accuracy of the energy determination with a factor of hundred, the energy invariance E is exact (see inset). The constant VA dynamics is obtained from two sets of positions, R 0,R 1 and h and these start values can not contain information about whether the system is in equilibrium 8

9 or not. In order to obtain the evolution of the kinetic energy Kn and the ability U n in a non-equilibrium system with VA dynamics a system was started with two sets of positions which correspond to a non-equilibrium state. The non-equilibrium state was obtained for a system of N = 2048 LJ particles in a fcc solid at kt/ǫ = 1 and density ρ = 1.009σ 3 ( density of coexisting solid at kt/ǫ = 1 [11]) by spontaneously expanding the positions to the density ρ = 0.80σ 3 by a scaling of all the positions R 0,R 1. A LJ systems equilibrium state at the density ρ = 0.80σ 3 is a liquid. The fcc ordered system melted spontaneously, and the conservative systems temperature decreased according to the Second law of Thermodynamics. The change in the temperature at the spontaneous melting is shown with green dashes in Figure 2. The temperature decreased from T1 = 1. within time steps to T 0.7. The differences between Tn and the temperature T n, obtained for the shadow Hamiltonian [8] are of the order 10 7, and they are not visible on the figure. The decrease in the spontaneous temperature at the melting was balanced by a corresponding increase in the ability (red line), and the energy invariance E (blue small dashes) was constant in the conservative system. The two MD simulations (Figure 1 and Figure 2) demonstrate that the traditional and the present (discrete) energy concept work equally well. III. DISCRETE DYNAMICS VERSUS ANALYTIC DYNAMICS The discrete VA dynamics has the same invariances as Newtonian dynamics and it raises the question: Which of these formulations that are correct, or alternatively, the most appropriate formulation of classical dynamics? In this context T. D. Lee in 1983 wrote a paper [12] entitled, Can Time Be a Discrete Dynamical Variable? ; which led to a series of publications by Lee and collaborators on the formulation of fundamental dynamics in terms of difference equations, but with exact invariance under continuous groups of translational and rotational transformations. Quoting Lee [13], he wish to explore an alternative point of view: that physics should be formulated in terms of difference equations and that these difference equations could exhibit all the desirable symmetry properties and conservation laws. Lee s analysis covers not only classical mechanics [12], but also non relativistic quantum mechanics and relativistic quantum field theory [14], and Gauge theory and Lattice Gravity [13]. The discrete dynamics is obtained by treating positions and time, but not momenta, as a discrete dynamical variables, and he obtained a conserved (mean) energy over con- 9

10 secutive time intervals of different lengths. But according to Lee [12] in his formulation of discrete mechanics, there is a fundamental length or time l (in natural units). Given any time interval T = t f t 0, the total number N of discrete points that define the trajectory is given by the integer nearest T/l. The analogy between Lee s formulation of discrete dynamics and VA dynamics is striking. FortheVAdynamics oneuses a unit time increment, h, andthe momenta arenotdynamical variables and they have no impact on the discrete dynamics [2]. The fundamental length and time in quantum electrodynamics are the Planck length l P m and Planck time t P s [15], and they are immensely smaller than the length unit (given by the floating point precision) and time increment used in MD to generate the classical discrete dynamics. But the analogy implies that the discrete VA dynamics obtained by MD is the continuation of the Lee s discrete quantum dynamics for a fundamental length of time t P, as is the analytic classical dynamics of the traditional quantum mechanics, given by the Wigner expansion [16]. The discrete non relativistic quantum mechanics is obtained by Lee using Feynman s path integration formalism, but for discrete positions and a corresponding discrete action, N+1 [ ] (Rn R n 1 ) 2 A D = +(t n t n 1 )V(n) 2(t n t n 1 ) n=1 (25) where R N+1 is the end-positions at time t N+1 and the minimum of A D determines the classical path. The action is a sum over products of time increments and kinetic energies Kn, and Lee has used the symbol V(n), for the average of potential energy in the time intervals [t n 1,t n ]. The momenta for all the paths, given by the discrete nodes R 1,..,R N+1 are obtained from differences, R n R n 1, so the classical VA discrete trajectory is the classical limit path for discrete quantum mechanics with h = t P, as the classical Newtonian trajectory is for the traditional quantum mechanics. There is, however, one important difference between the analytic and the discrete dynamics. The momenta for all the paths in the discrete quantum dynamics are obtained by a difference between discrete sets of positions and they are all asynchronous with the positions. So the Heisenberg uncertainty is a trivial consequence of a discrete quantum electrodynamics with a fundamental length of time t P. Lee motivates his reformulation of the analytic dynamics in the Introduction in [14] by the difficulties of formulating a general unifying theoretical model for dynamics and with 10

11 the Concluding remarks in [13] that (he tries to explore the opposite viewpoint): Difference equations are more fundamental, and differential equations are regarded as approximations. The difference in the energy between the analytic energy and the energy obtained by Eq. (21) for discrete electrodynamics with a unit time increment t P is of the order t 2 P, and it is absolute marginal. The Heisenberg uncertainty between positions and momenta is of the order t P and this uncertainty is an inherent quality of discrete dynamics with a fundamental lengthoftimet P. ThediscreteclassicalVAdynamicsisfundamentallydifferentfromanalytic Newtonian dynamics, but has the same invariances and the dynamics is obtained equally well by both methods. But, on the other hand the traditional quantum mechanics is in all manner fully appropriate and justifies no revision of the formulation, and an eventual revision of the dynamics must be justified by other facts than conservation of the energy by classical Molecular Dynamics simulation with the VA algorithm. ACKNOWLEDGMENTS The author acknowledges useful discussions with Ole J Heilmann and Jeppe C Dyre. The centre for viscous liquid dynamics Glass and Time is sponsored by the Danish National Research Foundation (DNRF) grant No. DNRF61. [1] L. Verlet, Phys. Rev. 159, 98 (1967). [2] S. Toxvaerd, O. J. Heilmann and J. C. Dyre, J. Chem. Phys. 136, (2012). [3] J. M. Sanz-Serna, Acta Numer. 1, 243 (1992). [4] E. Hairer, Ann. Numer. Math. 1, 107 (1994). [5] S. Reich, SIAM J. Numer. Anal. 36, 1549 (1999). [6] S. Toxvaerd, Phys Rev. E, 50, 2271 (1994). (The word shadow Hamiltonian was introduced in this paper, inspired by the terms a slightly perturbed Hamiltonian [H. Yoshida, Phys. Lett. A 150, 262 (1990)] and shadow trajectories [C. Grebogi, S. M. Hammel, J. A: Yorke and T. Saur, Phys. Rev. Lett. 65, 1527 (1990)]). 11

12 [7] E. Hairer, C. Lubich and G. Wanner Geometrical Numerical Integration (Springer Books Archives, 2006) [8] S. Toxvaerd, J. Chem. Phys. 139, (2013). [9] T. S. Ingebrigtsen, S. Toxvaerd, O. J. Heilmann, T. B. Schrøder, and J. C. Dyre, J. Chem. Phys. 135, (2011). [10] H. Goldstein, C. P. Poole, and J. Safko Classical Mechanics Third Edition. (Pearson, 2011). Chapter 1. [11] M. A. Barroso and A. L. Ferreira, J. Chem. Phys. 116, 7145 (2002). [12] T. D. Lee, Phys. Lett. 122 B, 217 (1983). [13] T. D. Lee, J. Stat. Phys. 46, 843 (1987). [14] R. Friedberg and T. D. Lee, Nucl. Phys. B 225 [FS9], 1 (1983). [15] See e.g. L. J. Garay, J. Mod. Phys. A 10, 145 (1995). [16] E. Wigner, Phys. Rev. 40, 749 (1932). 12

arxiv: v1 [cond-mat.other] 5 Oct 2015

arxiv: v1 [cond-mat.other] 5 Oct 2015 a) st@ruc.dk Dynamics of homogeneous nucleation Søren Toxvaerd a) arxiv:1510.01065v1 [cond-mat.other] 5 Oct 2015 DNRF centre Glass and Time, IMFUFA, Department of Sciences, Roskilde University, Postbox