1
|
Lavrinenko Y, Levashov PR, Minakov DV, Morozov IV, Valuev IA. Equilibrium properties of warm dense deuterium calculated by the wave packet molecular dynamics and density functional theory method. Phys Rev E 2021; 104:045304. [PMID: 34781451 DOI: 10.1103/physreve.104.045304] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2021] [Accepted: 09/13/2021] [Indexed: 11/07/2022]
Abstract
A joint simulation method based on the wave packet molecular dynamics and density functional theory (WPMD-DFT) is applied to study warm dense deuterium (nonideal deuterium plasmas). This method was developed recently as an extension of the wave packet molecular dynamics (WPMD) in which the equations of motion are solved simultaneously for classical ions and semiclassical electrons represented as Gaussian wave packets. Compared to the classical molecular dynamics and WPMD simulations, the method of WPMD-DFT provides a more accurate representation of quantum effects such as electron-ion coupling and electron degeneracy. It allows studying nonadiabatic dynamics of electrons and ions in equilibrium and nonequilibrium states while being more accurate and efficient at high densities than WPMD and classical molecular dynamics. In the paper, we discuss particular features of the method such as special boundary conditions and the procedure of isentrope calculation as well as the results obtained by WPMD-DFT for the shock-compressed deuterium. The compression isentrope and principal Hugoniot curves obtained by WPMD-DFT are compared with available experimental data and other simulation approaches to validate the method. It opens up a possibility of further application of the method to study nonequilibrium states and relaxation processes.
Collapse
Affiliation(s)
- Yaroslav Lavrinenko
- Joint Institute for High Temperatures of Russian Academy of Sciences, Moscow 125412, Russia.,Moscow Institute of Physics and Technology, Dolgoprudny 141701, Russia
| | - Pavel R Levashov
- Joint Institute for High Temperatures of Russian Academy of Sciences, Moscow 125412, Russia.,Moscow Institute of Physics and Technology, Dolgoprudny 141701, Russia
| | - Dmitry V Minakov
- Joint Institute for High Temperatures of Russian Academy of Sciences, Moscow 125412, Russia.,Moscow Institute of Physics and Technology, Dolgoprudny 141701, Russia
| | - Igor V Morozov
- Joint Institute for High Temperatures of Russian Academy of Sciences, Moscow 125412, Russia.,Moscow Institute of Physics and Technology, Dolgoprudny 141701, Russia
| | - Ilya A Valuev
- Joint Institute for High Temperatures of Russian Academy of Sciences, Moscow 125412, Russia
| |
Collapse
|
2
|
Arkhipov YV, Ashikbayeva A, Askaruly A, Davletov AE, Dubovtsev DY, Santybayev KS, Syzganbayeva SA, Conde L, Tkachenko IM. Dynamic characteristics of three-dimensional strongly coupled plasmas. Phys Rev E 2020; 102:053215. [PMID: 33327172 DOI: 10.1103/physreve.102.053215] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2020] [Accepted: 10/26/2020] [Indexed: 11/07/2022]
Abstract
The dynamic structure factor and other dynamic characteristics of strongly coupled one-component plasmas have been studied [Yu. V. Arkhipov et al., Phys. Rev. Lett. 119, 045001 (2017)PRLTAO0031-900710.1103/PhysRevLett.119.045001] using the self-consistent version of the method of moments. Within any version of the latter, the system dielectric function satisfies all involved sum rules and other exact relations automatically, and the advantage of this version is that, in addition, the dynamic characteristics (the dynamic structure factor, the dispersion, and decay parameters of the collective modes) are all expressed in terms of the static ones (the static structure factor) without any adjustment to the simulation data. The approach outlined in the aforementioned Letter is justified in detail and applied mainly to the classical Coulomb systems achieving satisfactory agreement with new numerical simulation data. It is shown how the realm of applicability of the method can be extended to partly degenerate and multicomponent systems, even to simple liquids. Some additional theoretical results are presented in the Supplemental Material.
Collapse
Affiliation(s)
- Yu V Arkhipov
- Department of Physics and Technology, IETP, Al-Farabi Kazakh National University, al-Farabi 71, 050040 Almaty, Kazakhstan
| | - A Ashikbayeva
- Department of Physics and Technology, IETP, Al-Farabi Kazakh National University, al-Farabi 71, 050040 Almaty, Kazakhstan
| | - A Askaruly
- Department of Physics and Technology, IETP, Al-Farabi Kazakh National University, al-Farabi 71, 050040 Almaty, Kazakhstan
| | - A E Davletov
- Department of Physics and Technology, IETP, Al-Farabi Kazakh National University, al-Farabi 71, 050040 Almaty, Kazakhstan
| | - D Yu Dubovtsev
- Department of Physics and Technology, IETP, Al-Farabi Kazakh National University, al-Farabi 71, 050040 Almaty, Kazakhstan
| | - Kh S Santybayev
- Department of Physics and Technology, IETP, Al-Farabi Kazakh National University, al-Farabi 71, 050040 Almaty, Kazakhstan
| | - S A Syzganbayeva
- Department of Physics and Technology, IETP, Al-Farabi Kazakh National University, al-Farabi 71, 050040 Almaty, Kazakhstan
| | - L Conde
- Departamento de Física Aplicada a la Ingeniería Aeronáutica, ETSIAE, Universidad Politécnica de Madrid, Plaza del Cardenal Cisneros 3, 28040 Madrid, Spain
| | - I M Tkachenko
- Departament de Matemàtica Aplicada, Universitat Politècnica de València, Camino de Vera s/n, 46022 Valencia, Spain
| |
Collapse
|