CONFINING BOUNDARY CONDITIONS FOR SIMULATION OF ELECTRON-ION PLASMA BY ANTISYMMETRIZED WAVE PACKET MOLECULAR DYNAMICS
Valuev I.A., Morozov I.V.
Joint Institute for High Temperatures of RAS, 125412 Izhorskaya 13 bld. 2,
Moscow, Russia
The method of Wave Packet Molecular dynamics [1] is an approximate quantum method for numerical simulation of many-particle dynamics. In this method single electron wave functions are expanded in a set [2] of floating Gaussian wave packets (WP). A trial many-body wave function is constructed depending on the quantum statistical properties of the simulated ensemble, for example for fermions (electrons) a single Slater determinant antisymmetrized product is usually used for each spin projection. The resulting equations of motion follow from the variational principle [3].
In this work we study a system of electrons in a 3D confinement potential constructed from harmonic walls and a flat floor in each spatial direction. The matrix elements of the proposed potential with Gaussians were obtained analytically. The confined system constructed this way is an excellent testing ground for AWPMD method, because it poses correct quantum boundary conditions; reduces to textbook many-body quantum oscillator in the limit of no interaction; reduces to uniform quantum gas when the floor size is large. The unlimited broadening of the Gaussian wave packets and underestimation of the electron-electron and electron-ion collision frequencies is known to be the major problems of the WPMD method [4] when applied to many-particle systems with homogeneous density, for example plasma systems. We show however that the method is able to excellently describe the thermodynamics of a confined fermionic system. The infinite WP broadening is directly related to the infinite statistical sum in an unconstrained system and does not appear in the confined system. For example, the nearest image boundary conditions guarantee the periodicity of expectation values but do not limit the number of allowed states in the simulated system.
We show that for the confined electron-ion system the AWPMD model does not suffer from the unlimited wave packet broadening problem and the electron wave packet parameters perform finite ergodic motion in the available system phase space. The implementation of the proposed model to simulation of the extended warm dense matter is discussed.
1. A.A. Frost, J. Chem. Phys., 47, 3707 (1967)
2. I. V. Morozov and I. A. Valuev, Contrib. Plasma Phys. 52, 140 (2012)
3. D. Klakow, C. Toepffer, P. G. Reinhard, J. Chem. Phys., 101, 10766 (1994)
4. I.V. Morozov, I. A. Valuev, J. Phys. A, 42, 214044 (2009)
5. H. Feldmeier, J. Schnack, Rev. Modern Phys., 72, 655 (2000)