1
|
Salazar R, Cobos C, Jaramillo D, Bayona-Roa C. Comparative analysis of molecular dynamics and method of moments in two-dimensional concentric circular layers. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2024; 36:395603. [PMID: 38917836 DOI: 10.1088/1361-648x/ad5baf] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/29/2024] [Accepted: 06/25/2024] [Indexed: 06/27/2024]
Abstract
In this manuscript, we undertake an examination of a classical plasma deployed on two finite co-planar surfaces: a circular regionΩininto an annular regionΩoutwith a gap in between. It is studied both from the point of view of statistical mechanics and the electrostatics of continua media. We employ a dual perspective: the first one is by using molecular dynamics (MD) simulations to find the system's positional correlation functions and velocity distributions. That by modeling the system as a classical two-dimensional Coulomb plasma of point-like charged particlesq1andq2on the layersΩinandΩoutrespectively with no background density. The second one corresponds to a finite Surface Electrode (SE) composed of planar metallic layers displayed on the regionsΩin,Ωoutat constant voltagesVin,Voutconsidering axial symmetry. The surface charge density is calculated by the Method of Moments (MoM) under the electrostatic approximation. Point-like and differential charges elements interact via a1/r-electric potential in both cases. The thermodynamic averages of the number density, and electric potential due to the plasma depend on the coupling and the charge ratioξ=q1/q2once the geometry of the layers is fixed. On the other hand, the fields due to the SE depend on the layer's geometry and their voltage. In the document, is defined a protocol to properly compare the systems. We show that there are values of the coupling parameter, where the thermodynamic averages computed via MD agree with the results of MoM for attractiveξ=-1and repulsive layersξ = 1.
Collapse
Affiliation(s)
- Robert Salazar
- Dirección de Ingeniería Electrónica, Universidad ECCI, Bogotá, Colombia
| | - Cristian Cobos
- Dirección de Ingeniería Electrónica, Universidad ECCI, Bogotá, Colombia
| | - Diego Jaramillo
- Departamento de Ciencias, Universidad Antonio Nariño, Bogotá, Colombia
| | - Camilo Bayona-Roa
- Centro de Ingeniería Avanzada Investigación y Desarrollo-CIAID, Bogotá, Colombia
| |
Collapse
|
2
|
Filinov AV, Ara J, Tkachenko IM. Dynamic properties and the roton mode attenuation in liquid 3He: ab initio study within the self-consistent method of moments. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2023; 381:20220324. [PMID: 37393939 DOI: 10.1098/rsta.2022.0324] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/29/2022] [Accepted: 02/25/2023] [Indexed: 07/04/2023]
Abstract
The dynamic structure factor and the eigenmodes of density fluctuations in liquid 3He are studied using a novel non-perturbative approach. This new version of the self-consistent method of moments invokes up to nine sum rules and other exact relations, the two-parameter Shannon information entropy maximization procedure, and the ab initio path integral Monte Carlo simulations which provide necessary reliable input information on the system static properties. Detailed analysis is performed of the collective excitations dispersion relations, the modes' decrements and the static structure factor of 3He at the saturated vapour pressure. The results are compared to available experimental data by Albergamo et al. (Albergamo et al. 2007 Phys. Rev. Lett. 99, 205301. (doi:10.1103/PhysRevLett.99.205301)) and Fåk et al. (Fåk et al. 1994 J. Low Temp. Phys. 97, 445-487. (doi:10.1007/BF00754303)). The theory reveals a clear signature of the roton-like feature in the particle-hole segment of the excitation spectrum with a significant reduction of the roton decrement in the wavenumber range [Formula: see text]. The observed roton mode remains a well-defined collective mode even in the particle-hole band, where it is strongly damped. The existence of the roton-like mode in the bulk liquid 3He is confirmed like in other quantum fluids. The phonon branch of the spectrum is in a reasonable agreement with the same experimental data. This article is part of the theme issue 'Dynamic and transient processes in warm dense matter'.
Collapse
Affiliation(s)
- A V Filinov
- Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, Kiel, Germany
| | - J Ara
- Instituto de Tecnología Química, Universitat Politècnica de València, Valencia, Spain
| | - I M Tkachenko
- Departament de Matemàtica Aplicada, Universitat Poliècnica de València, Valencia, Spain
- Al-Farabi Kazakh National University, Almaty, Kazakhstan
| |
Collapse
|
3
|
Fairushin II, Mokshin AV. Collective ion dynamics in Coulomb one-component plasmas within the self-consistent relaxation theory. Phys Rev E 2023; 108:015206. [PMID: 37583226 DOI: 10.1103/physreve.108.015206] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2023] [Accepted: 07/10/2023] [Indexed: 08/17/2023]
Abstract
In this paper, we present the theoretical formalism describing the collective ion dynamics of the nonideal Coulomb classical one-component plasmas on the basis of the self-consistent relaxation theory. The theory is adapted to account for correlations between the frequency relaxation parameters that characterize the three- and four-particle dynamics and the parameters associated with the two-particle dynamics. The dynamic structure factor spectra and dispersion characteristics calculated for a wide range of wave numbers are in agreement with the molecular dynamics simulation data and the results obtained with the theory of the frequency moments. The proposed formalism reproduces all the features inherent to the Coulomb one-component plasmas and requires only knowledge of the coupling parameter and the information about the structure.
Collapse
Affiliation(s)
- Ilnaz I Fairushin
- Department of Computational Physics, Institute of Physics, Kazan Federal University, 420008 Kazan, Russia
| | - Anatolii V Mokshin
- Department of Computational Physics, Institute of Physics, Kazan Federal University, 420008 Kazan, Russia
| |
Collapse
|
4
|
Dhaka A, Subhash PV, Bandyopadhyay P, Sen A. Auto-correlations of microscopic density fluctuations for Yukawa fluids in the generalized hydrodynamics framework with viscoelastic effects. Sci Rep 2022; 12:21883. [PMID: 36536026 PMCID: PMC9763351 DOI: 10.1038/s41598-022-26401-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2022] [Accepted: 12/14/2022] [Indexed: 12/23/2022] Open
Abstract
The present work develops a theoretical procedure for obtaining transport coefficients of Yukawa systems from density fluctuations. The dynamics of Yukawa systems are described in the framework of the generalized hydrodynamic (GH) model that incorporates strong coupling and visco-elastic memory effects by using an exponentially decaying memory function in time. A hydrodynamic matrix for such a system is exactly derived and then used to obtain an analytic expression for the density autocorrelation function (DAF)-a marker of the time dynamics of density fluctuations. The present approach is validated against a DAF obtained from numerical data of Molecular Dynamics (MD) simulations of a dusty plasma system that is a practical example of a Yukawa system. The MD results and analytic expressions derived from the model equations are then used to obtain various transport coefficients and the latter are compared with values available in the literature from other models. The influence of strong coupling and visco-elastic effects on the transport parameters are discussed. Finally, the utility of our calculations for obtaining reliable estimates of transport coefficients from experimentally determined DAF is pointed out.
Collapse
Affiliation(s)
- Ankit Dhaka
- grid.502813.d0000 0004 1796 2986Institute for Plasma Research, Bhat, Gandhinagar, Gujarat 382428 India ,grid.450257.10000 0004 1775 9822Homi Bhabha National Institute, Training School Complex, Anushaktinagar, Mumbai, 400094 India
| | - P. V. Subhash
- grid.450257.10000 0004 1775 9822Homi Bhabha National Institute, Training School Complex, Anushaktinagar, Mumbai, 400094 India ,grid.502813.d0000 0004 1796 2986ITER-India, Institute for Plasma Research, Bhat, Gandhinagar, Gujarat 382428 India
| | - P. Bandyopadhyay
- grid.502813.d0000 0004 1796 2986Institute for Plasma Research, Bhat, Gandhinagar, Gujarat 382428 India ,grid.450257.10000 0004 1775 9822Homi Bhabha National Institute, Training School Complex, Anushaktinagar, Mumbai, 400094 India
| | - A. Sen
- grid.502813.d0000 0004 1796 2986Institute for Plasma Research, Bhat, Gandhinagar, Gujarat 382428 India ,grid.450257.10000 0004 1775 9822Homi Bhabha National Institute, Training School Complex, Anushaktinagar, Mumbai, 400094 India
| |
Collapse
|
5
|
Dynamical Density of Two-Dimensional Dusty Plasmas. IRANIAN JOURNAL OF SCIENCE AND TECHNOLOGY, TRANSACTIONS A: SCIENCE 2022. [DOI: 10.1007/s40995-022-01386-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
6
|
Mokshin AV, Fairushin II, Tkachenko IM. Self-consistent relaxation theory of collective ion dynamics in Yukawa one-component plasmas under intermediate screening regimes. Phys Rev E 2022; 105:025204. [PMID: 35291083 DOI: 10.1103/physreve.105.025204] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2021] [Accepted: 02/01/2022] [Indexed: 06/14/2023]
Abstract
The self-consistent relaxation theory is employed to describe the collective ion dynamics in strongly coupled Yukawa classical one-component plasmas. The theory is applied to equilibrium states corresponding to intermediate screening regimes with appropriate values of the structure and coupling parameters. The information about the structure (the radial distribution function and the static structure factor) and the thermodynamics of the system are sufficient to describe collective dynamics over a wide range of spatial scales, namely, from the extended hydrodynamic to the microscopic dynamics scale. The main experimentally measurable characteristics of the equilibrium collective dynamics of ions-the spectrum of the dynamic structure factor, the dispersion parameters, the speed of sound, and the sound attenuation-are determined within the framework of the theory without using any adjustable parameters. The results demonstrate agreement with molecular dynamics simulations. Thus a direct realization is presented of the key idea of statistical mechanics: for the theoretical description of the collective particle dynamics in equilibrium fluids it is sufficient to know the interparticle interaction potential and the structural characteristics. Comparison with alternative or complementary theoretical approaches is provided.
Collapse
Affiliation(s)
- Anatolii V Mokshin
- Department of Computational Physics, Institute of Physics, Kazan Federal University, 420008 Kazan, Russia
| | - Ilnaz I Fairushin
- Department of Computational Physics, Institute of Physics, Kazan Federal University, 420008 Kazan, Russia
| | - Igor M Tkachenko
- Departament de Matemàtica Aplicada, Universitat Politècnica de València, Camino de Vera s/n, 46022 Valencia, Spain
- Al-Farabi Kazakh National University, al-Farabi Av. 71, 050040 Almaty, Kazakhstan
| |
Collapse
|
7
|
Lucco Castello F, Tolias P. Bridge functions of classical one-component plasmas. Phys Rev E 2022; 105:015208. [PMID: 35193199 DOI: 10.1103/physreve.105.015208] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2021] [Accepted: 01/10/2022] [Indexed: 06/14/2023]
Abstract
In a recent paper, Lucco Castello et al. [arXiv:2107.03537] performed systematic extractions of classical one-component plasma bridge functions from molecular dynamics simulations and provided an accurate parametrization that was incorporated in their isomorph-based empirically modified hypernetted chain approach for Yukawa one-component plasmas. Here the extraction technique and parametrization strategy are described in detail, while the deficiencies of earlier efforts are discussed. The structural and thermodynamic predictions of the updated version of the integral equation theory approach are compared with extensive available simulation results revealing a truly unprecedented level of accuracy in the entire dense liquid region of the Yukawa phase diagram.
Collapse
Affiliation(s)
- F Lucco Castello
- Space and Plasma Physics, Royal Institute of Technology, Stockholm SE-100 44, Sweden
| | - P Tolias
- Space and Plasma Physics, Royal Institute of Technology, Stockholm SE-100 44, Sweden
| |
Collapse
|
8
|
Shahzad A, Manzoor A, Wang W, Mahboob A, Kashif M, He MG. Dynamic Characteristics of Strongly Coupled Nonideal Plasmas. ARABIAN JOURNAL FOR SCIENCE AND ENGINEERING 2021. [DOI: 10.1007/s13369-021-05954-4] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
9
|
Magyar P, Kalman GJ, Hartmann P, Donkó Z. Anomalous behavior of plasma response functions at strong coupling. Phys Rev E 2021; 104:015202. [PMID: 34412281 DOI: 10.1103/physreve.104.015202] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2021] [Accepted: 06/08/2021] [Indexed: 06/13/2023]
Abstract
Using data from equilibrium molecular dynamics computer simulations we have built up a catalog of response functions for the Coulomb one-component plasma over a wide range of Γ coupling values, including the strongly coupled Γ>1 liquid regime. We focus on the domain of negative compressibility (Γ>3), where the proper response displays an acausal behavior, implying a modification of the relation between its real and imaginary parts in the Kramers-Kronig relations. We give a description of the details of this acausal feature, in both the frequency and time domains. We show that the viscoelastic pole of the density response function morphs into an imaginary pole in the upper ω half-plane that is responsible for the anomalous behavior of the response in this coupling range. By examining the plasmon dispersion relation through the dielectric response function, rather than via the peaks of the dynamical structure function, we obtain a more reliable representation for the dispersion. We demonstrate that there is an intimate link between the formation of the roton minimum in the dispersion and the negative compressibility of the system. The feasibility of the extension of our analysis to systems with a short-range interaction is explored.
Collapse
Affiliation(s)
- Péter Magyar
- Institute for Solid State Physics and Optics, Wigner Research Centre for Physics, P.O.B. 49, H-1525 Budapest, Hungary
| | - Gabor J Kalman
- Department of Physics, Boston College, Chestnut Hill, Massachusetts 20467, USA
| | - Péter Hartmann
- Institute for Solid State Physics and Optics, Wigner Research Centre for Physics, P.O.B. 49, H-1525 Budapest, Hungary
| | - Zoltán Donkó
- Institute for Solid State Physics and Optics, Wigner Research Centre for Physics, P.O.B. 49, H-1525 Budapest, Hungary
| |
Collapse
|
10
|
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
|
11
|
Dornheim T, Vorberger J. Finite-size effects in the reconstruction of dynamic properties from ab initio path integral Monte Carlo simulations. Phys Rev E 2020; 102:063301. [PMID: 33466040 DOI: 10.1103/physreve.102.063301] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2020] [Accepted: 11/15/2020] [Indexed: 06/12/2023]
Abstract
We systematically investigate finite-size effects in the dynamic structure factor S(q,ω) of the uniform electron gas obtained via the analytic continuation of ab initio path integral Monte Carlo data for the imaginary-time density-density correlation function F(q,τ). Using the recent scheme by Dornheim et al. [Phys. Rev. Lett. 121, 255001 (2018)PRLTAO0031-900710.1103/PhysRevLett.121.255001], we find that the reconstructed spectra are not afflicted with any finite-size effects for as few as N=14 electrons both at warm dense matter (WDM) conditions and at the margins of the strongly correlated electron liquid regime. Our results further corroborate the high quality of our current description of the dynamic density response of correlated electrons, which is of high importance for many applications in WDM theory and beyond.
Collapse
Affiliation(s)
- Tobias Dornheim
- Center for Advanced Systems Understanding (CASUS), D-02826 Görlitz, Germany
| | - Jan Vorberger
- Helmholtz-Zentrum Dresden-Rossendorf (HZDR), D-01328 Dresden, Germany
| |
Collapse
|
12
|
Pan H, Kalman GJ, Hartmann P, Donkó Z. Strongly coupled Yukawa trilayer liquid: Structure and dynamics. Phys Rev E 2020; 102:043206. [PMID: 33212692 DOI: 10.1103/physreve.102.043206] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2020] [Accepted: 09/09/2020] [Indexed: 11/07/2022]
Abstract
The equilibrium structure and the dispersion relations of collective excitations in trilayer Yukawa systems in the strongly coupled liquid regime are examined. The equilibrium correlations reveal a variety of structures in the liquid phase, reminiscent of the corresponding structures in the solid phase. At small layer separation substitutional disorder becomes the governing feature. Theoretical dispersion relations are obtained by applying the quasilocalized charge approximation (QLCA) formalism, while numerical data are generated by microcanonical molecular dynamics simulations. The dispersions and polarizations of the collective excitations obtained through both of these methods are compared and discussed in detail. We find that the QLCA method is, in general, very satisfactory, but that there are phenomena not covered by the QLCA. In particular, by analyzing the dynamical longitudinal and transverse current fluctuation spectra we discover the existence of a structure not related to the collective mode spectra. This also provides insight into the long-standing problem of the gap frequency discrepancy, observed in strongly coupled layered systems in earlier studies.
Collapse
Affiliation(s)
- Hong Pan
- Department of Physics, Boston College, Chestnut Hill, Massachusetts 02467, USA
| | - Gabor J Kalman
- Department of Physics, Boston College, Chestnut Hill, Massachusetts 02467, USA
| | - Peter Hartmann
- Institute for Solid State Physics and Optics, Wigner Research Centre for Physics, P.O. Box 49, H-1525 Budapest, Hungary.,Center for Astrophysics, Space Physics and Engineering Research (CASPER), Baylor University, 100 Research Pkwy, Waco, Texas 76706, USA
| | - Zoltán Donkó
- Institute for Solid State Physics and Optics, Wigner Research Centre for Physics, P.O. Box 49, H-1525 Budapest, Hungary
| |
Collapse
|
13
|
Wang ZQ, Tang J, Hou Y, Chen QF, Chen XR, Dai JY, Meng XJ, Gu YJ, Liu L, Li GJ, Lan YS, Li ZG. Benchmarking the effective one-component plasma model for warm dense neon and krypton within quantum molecular dynamics simulation. Phys Rev E 2020; 101:023302. [PMID: 32168678 DOI: 10.1103/physreve.101.023302] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2019] [Accepted: 01/15/2020] [Indexed: 11/07/2022]
Abstract
The effective one-component plasma (EOCP) model has provided an efficient approach to obtaining many important thermophysical parameters of hot dense matter [J. Clérouin, et al., Phys. Rev. Lett. 116, 115003 (2016)PRLTAO0031-900710.1103/PhysRevLett.116.115003]. In this paper, we perform extensive quantum molecular dynamics (QMD) simulations to determine the equations of state, ionic structures, and ionic transport properties of neon and krypton within the warm dense matter (WDM) regime where the density (ρ) is up to 12 g/cm^{3} and the temperature (T) is up to 100 kK. The simulated data are then used as a benchmark to explicitly evaluate the EOCP and Yukawa models. It is found that, within present ρ-T regime, the EOCP model can excellently reproduce the diffusion and viscosity coefficients of neon and krypton due to the fact that this model defines a system which nearly reproduces the actual physical states of WDM. Therefore, the EOCP model may be a promising alternative approach to reasonably predicting the transport behaviors of matter in WDM regime at lower QMD computational cost. The evaluation of Yukawa model shows that the consideration of the energy level broadening effect in the average atom model is necessary. Finally, with the help of EOCP model, the Stokes-Einstein relationships about neon and krypton are discussed, and fruitful plasma parameters as well as a practical ρ-T-dependent formula of the effective coupling parameter are obtained. These results not only provide valuable information for future theoretical and experimental studies on dense neon and krypton but also reveal the applicability of the EOCP model and the limitation of the Yukawa model in WDM regime and further support the continuing search for a unified description of ionic transport in dense plasma.
Collapse
Affiliation(s)
- Zhao-Qi Wang
- College of Physics, Sichuan University, Chengdu 610065, People's Republic of China.,National Key Laboratory for Shock Wave and Detonation Physics Research, Institute of Fluid Physics, Chinese Academy of Engineering Physics, Mianyang 621900, People's Republic of China
| | - Jun Tang
- Science and Technology on Surface Physics and Chemistry Laboratory, Mianyang 621908, People's Republic of China
| | - Yong Hou
- Department of Physics, National University of Defense Technology, Changsha 410073, People's Republic of China
| | - Qi-Feng Chen
- National Key Laboratory for Shock Wave and Detonation Physics Research, Institute of Fluid Physics, Chinese Academy of Engineering Physics, Mianyang 621900, People's Republic of China
| | - Xiang-Rong Chen
- College of Physics, Sichuan University, Chengdu 610065, People's Republic of China
| | - Jia-Yu Dai
- Department of Physics, National University of Defense Technology, Changsha 410073, People's Republic of China
| | - Xu-Jun Meng
- Institute of Applied Physics and Computational Mathematics, Beijing 100094, People's Republic of China
| | - Yun-Jun Gu
- National Key Laboratory for Shock Wave and Detonation Physics Research, Institute of Fluid Physics, Chinese Academy of Engineering Physics, Mianyang 621900, People's Republic of China
| | - Lei Liu
- College of Physics, Sichuan University, Chengdu 610065, People's Republic of China.,National Key Laboratory for Shock Wave and Detonation Physics Research, Institute of Fluid Physics, Chinese Academy of Engineering Physics, Mianyang 621900, People's Republic of China
| | - Guo-Jun Li
- College of Physics, Sichuan University, Chengdu 610065, People's Republic of China.,National Key Laboratory for Shock Wave and Detonation Physics Research, Institute of Fluid Physics, Chinese Academy of Engineering Physics, Mianyang 621900, People's Republic of China
| | - Yang-Shun Lan
- College of Physics, Sichuan University, Chengdu 610065, People's Republic of China.,National Key Laboratory for Shock Wave and Detonation Physics Research, Institute of Fluid Physics, Chinese Academy of Engineering Physics, Mianyang 621900, People's Republic of China
| | - Zhi-Guo Li
- National Key Laboratory for Shock Wave and Detonation Physics Research, Institute of Fluid Physics, Chinese Academy of Engineering Physics, Mianyang 621900, People's Republic of China
| |
Collapse
|
14
|
Choi Y, Dharuman G, Murillo MS. High-frequency response of classical strongly coupled plasmas. Phys Rev E 2019; 100:013206. [PMID: 31499843 DOI: 10.1103/physreve.100.013206] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2019] [Indexed: 06/10/2023]
Abstract
The dynamic structure factor (DSF) of the Yukawa system is here obtained with highly converged molecular dynamics (MD) over the entire liquid phase. The data provide a rigorous test of theoretical models of ion-acoustic wave-dispersion relations, the intermediate scattering function, and the high-frequency response. We compare our MD results with seven diverse models, finding good agreement among those that enforce the three basic sum rules for dispersion properties, although one of the models has previously unreported spurious peaks. The MD simulations reveal that at intermediate frequencies ω, the high-frequency response of the DSF follows a power law, going approximately as ω^{-p}, where p>0, and p shows nontrivial dependencies on the wave vector q and the plasma parameters κ and Γ. In contrast, among the seven comparison models, the predicted high-frequency response is found to be independent of {q,κ,Γ}. This high-frequency power suggests a useful fitting form. In addition, these results expose limitations of several models and, moreover, suggest that some approaches are difficult or impossible to extend because of the lack of finite moments. We also find the double-plasmon resonance peak in our MD simulations that none of the theoretical models predicts.
Collapse
Affiliation(s)
- Yongjun Choi
- Institute for Cyber-Enabled Research, Michigan State University, East Lansing, Michigan 48824, USA
| | - Gautham Dharuman
- Computational Mathematics, Science and Engineering, Michigan State University, East Lansing, Michigan 48824, USA
| | - Michael S Murillo
- Computational Mathematics, Science and Engineering, Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|
15
|
Khrapak SA, Khrapak AG, Kryuchkov NP, Yurchenko SO. Onset of transverse (shear) waves in strongly-coupled Yukawa fluids. J Chem Phys 2019; 150:104503. [PMID: 30876343 DOI: 10.1063/1.5088141] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A simple practical approach to describe transverse (shear) waves in strongly-coupled Yukawa fluids is presented. Theoretical dispersion curves, based on hydrodynamic consideration, are shown to compare favorably with existing numerical results for plasma-related systems in the long-wavelength regime. The existence of a minimum wave number below which shear waves cannot propagate and its magnitude are properly accounted in the approach. The relevance of the approach beyond plasma-related Yukawa fluids is demonstrated by using experimental data on transverse excitations in liquid metals Fe, Cu, and Zn, obtained from inelastic x-ray scattering. Some potentially important relations, scalings, and quasi-universalities are discussed. The results should be interesting for a broad community in chemical physics, materials physics, physics of fluids and glassy state, complex (dusty) plasmas, and soft matter.
Collapse
Affiliation(s)
- Sergey A Khrapak
- Institut für Materialphysik im Weltraum, Deutsches Zentrum für Luft- und Raumfahrt (DLR), 82234 Weßling, Germany
| | - Alexey G Khrapak
- Joint Institute for High Temperatures, Russian Academy of Sciences, 125412 Moscow, Russia
| | | | | |
Collapse
|
16
|
Dornheim T, Groth S, Vorberger J, Bonitz M. Ab initio Path Integral Monte Carlo Results for the Dynamic Structure Factor of Correlated Electrons: From the Electron Liquid to Warm Dense Matter. PHYSICAL REVIEW LETTERS 2018; 121:255001. [PMID: 30608805 DOI: 10.1103/physrevlett.121.255001] [Citation(s) in RCA: 42] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/31/2018] [Indexed: 06/09/2023]
Abstract
The accurate description of electrons at extreme density and temperature is of paramount importance for, e.g., the understanding of astrophysical objects and inertial confinement fusion. In this context, the dynamic structure factor S(q,ω) constitutes a key quantity as it is directly measured in x-ray Thomson scattering experiments and governs transport properties like the dynamic conductivity. In this work, we present the first ab initio results for S(q,ω) by carrying out extensive path integral Monte Carlo simulations and developing a new method for the required analytic continuation, which is based on the stochastic sampling of the dynamic local field correction G(q,ω). In addition, we find that the so-called static approximation constitutes a promising opportunity to obtain high-quality data for S(q,ω) over substantial parts of the warm dense matter regime.
Collapse
Affiliation(s)
- T Dornheim
- Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, Leibnizstraße 15, D-24098 Kiel, Germany
| | - S Groth
- Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, Leibnizstraße 15, D-24098 Kiel, Germany
| | - J Vorberger
- Helmholtz-Zentrum Dresden-Rossendorf, D-01328 Dresden, Germany
| | - M Bonitz
- Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, Leibnizstraße 15, D-24098 Kiel, Germany
| |
Collapse
|