1
|
Gao CZ, Zhang CB, Cai Y, Wu Y, Fan ZF, Wang P, Wang JG. Assessment of the electron-proton energy relaxation rates extracted from molecular dynamics simulations in weakly-coupled hydrogen plasmas. Phys Rev E 2023; 107:015203. [PMID: 36797881 DOI: 10.1103/physreve.107.015203] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2022] [Accepted: 12/23/2022] [Indexed: 06/18/2023]
Abstract
Electron-proton energy relaxation rates are assessed using molecular dynamics (MD) simulations in weakly-coupled hydrogen plasmas. To this end, we use various approaches to extract the energy relaxation rate from MD-simulated temperatures, and we find that existing extracting approaches may yield results with a sizable discrepancy larger than the variance between analytical models, which is further verified by well-known case studies. Present results show that two of the extracting approaches can produce identical results, which is attributed to a proper treatment of relaxation evolution. To discriminate the use of various methods, an empirical criterion with respect to initial plasma temperatures is proposed, which can self-consistently explain the cases considered. In addition, for a transient electron-proton plasma, we show that it is possible to extrapolate the Coulomb logarithm from that derived by initial plasma parameters in a single MD calculation, which is reasonably consistent with previous MD data. Our results are helpful to obtain accurate MD-based energy relaxation rates.
Collapse
Affiliation(s)
- Cong-Zhang Gao
- Institute of Applied Physics and Computational Mathematics, Beijing 100088, People's Republic of China
| | - Cun-Bo Zhang
- Institute of Applied Physics and Computational Mathematics, Beijing 100088, People's Republic of China
| | - Ying Cai
- Institute of Applied Physics and Computational Mathematics, Beijing 100088, People's Republic of China
| | - Yong Wu
- Institute of Applied Physics and Computational Mathematics, Beijing 100088, People's Republic of China
| | - Zheng-Feng Fan
- Institute of Applied Physics and Computational Mathematics, Beijing 100088, People's Republic of China
| | - Pei Wang
- Institute of Applied Physics and Computational Mathematics, Beijing 100088, People's Republic of China
| | - Jian-Guo Wang
- Institute of Applied Physics and Computational Mathematics, Beijing 100088, People's Republic of China
| |
Collapse
|
2
|
Adrian PJ, Florido R, Grabowski PE, Mancini R, Bachmann B, Benedict LX, Johnson MG, Kabadi N, Lahmann B, Li CK, Petrasso RD, Rinderknecht HG, Regan SP, Séguin FH, Singleton RL, Sio H, Sutcliffe GD, Whitley HD, Frenje JA. Measurements of ion-electron energy-transfer cross section in high-energy-density plasmas. Phys Rev E 2022; 106:L053201. [PMID: 36559377 DOI: 10.1103/physreve.106.l053201] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2021] [Accepted: 09/02/2022] [Indexed: 06/17/2023]
Abstract
We report on measurements of the ion-electron energy-transfer cross section utilizing low-velocity ion stopping in high-energy-density plasmas at the OMEGA laser facility. These measurements utilize a technique that leverages the close relationship between low-velocity ion stopping and ion-electron equilibration. Shock-driven implosions of capsules filled with D^{3}He gas doped with a trace amount of argon are used to generate densities and temperatures in ranges from 1×10^{23} to 2×10^{24} cm^{-3} and from 1.4 to 2.5 keV, respectively. The energy loss of 1-MeV DD tritons and 3.7-MeV D^{3}He alphas that have velocities lower than the average velocity of the thermal electrons is measured. The energy loss of these ions is used to determine the ion-electron energy-transfer cross section, which is found to be in excellent agreement with quantum-mechanical calculations in the first Born approximation. This result provides an experimental constraint on ion-electron energy transfer in high-energy-density plasmas, which impacts the modeling of alpha heating in inertial confinement fusion implosions, magnetic-field advection in stellar atmospheres, and energy balance in supernova shocks.
Collapse
Affiliation(s)
- P J Adrian
- Plasma Science and Fusion Center, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - R Florido
- iUNAT, Departamento de Física, Universidad de Las Palmas de Gran Canaria, 35017 Las Palmas de Gran Canaria, Spain
| | - P E Grabowski
- Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - R Mancini
- Department of Physics, University of Nevada, Reno, Reno, Nevada 89557, USA
| | - B Bachmann
- Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - L X Benedict
- Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - M Gatu Johnson
- Plasma Science and Fusion Center, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - N Kabadi
- Plasma Science and Fusion Center, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - B Lahmann
- Plasma Science and Fusion Center, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - C K Li
- Plasma Science and Fusion Center, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - R D Petrasso
- Plasma Science and Fusion Center, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - H G Rinderknecht
- Laboratory for Laser Energetics, University of Rochester, Rochester, New York 14623, USA
| | - S P Regan
- Laboratory for Laser Energetics, University of Rochester, Rochester, New York 14623, USA
| | - F H Séguin
- Plasma Science and Fusion Center, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - R L Singleton
- SavantX Research Center, Santa Fe, New Mexico 87501, USA
- School of Mathematics, University of Leeds, Leeds LS2 9JT, United Kingdom
| | - H Sio
- Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - G D Sutcliffe
- Plasma Science and Fusion Center, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - H D Whitley
- Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - J A Frenje
- Plasma Science and Fusion Center, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| |
Collapse
|
3
|
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
|
4
|
Dornheim T, Vorberger J, Bonitz M. Nonlinear Electronic Density Response in Warm Dense Matter. PHYSICAL REVIEW LETTERS 2020; 125:085001. [PMID: 32909774 DOI: 10.1103/physrevlett.125.085001] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/09/2020] [Accepted: 07/20/2020] [Indexed: 06/11/2023]
Abstract
Warm dense matter (WDM)-an extreme state with high temperatures and densities that occurs, e.g., in astrophysical objects-constitutes one of the most active fields in plasma physics and materials science. These conditions can be realized in the lab by shock compression or laser excitation, and the most accurate experimental diagnostics is achieved with lasers and free electron lasers which is theoretically modeled using linear response theory. Here, we present first ab initio path integral Monte Carlo results for the nonlinear density response of correlated electrons in WDM and show that for many situations of experimental relevance nonlinear effects cannot be neglected.
Collapse
Affiliation(s)
- Tobias Dornheim
- Center for Advanced Systems Understanding (CASUS), D-028262 Görlitz, Germany
| | - Jan Vorberger
- Helmholtz-Zentrum Dresden-Rossendorf, Bautzner Landstraße 400, D-01328 Dresden, Germany
| | - Michael Bonitz
- Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, Leibnizstraße 15, D-24098 Kiel, Germany
| |
Collapse
|
5
|
Sprenkle T, Dodson A, McKnight Q, Spencer R, Bergeson S, Diaw A, Murillo MS. Ion friction at small values of the Coulomb logarithm. Phys Rev E 2019; 99:053206. [PMID: 31212549 DOI: 10.1103/physreve.99.053206] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2018] [Indexed: 11/07/2022]
Abstract
Transport properties of high-energy-density plasmas are influenced by the ion collision rate. Traditionally, this rate involves the Coulomb logarithm, lnΛ. Typical values of lnΛ are ≈10-20 in kinetic theories where transport properties are dominated by weak-scattering events caused by long-range forces. The validity of these theories breaks down for strongly coupled plasmas, when lnΛ is of order one. We present measurements and simulations of collision data in strongly coupled plasmas when lnΛ is small. Experiments are carried out in the first dual-species ultracold neutral plasma (UNP), using Ca^{+} and Yb^{+} ions. We find strong collisional coupling between the different ion species in the bulk of the plasma. We simulate the plasma using a two-species fluid code that includes Coulomb logarithms derived from either a screened Coulomb potential or a the potential of mean force. We find generally good agreement between the experimental measurements and the simulations. With some improvements, the mixed Ca^{+} and Yb^{+} dual-species UNP will be a promising platform for testing theoretical expressions for lnΛ and collision cross-sections from kinetic theories through measurements of energy relaxation, stopping power, two-stream instabilities, and the evolution of sculpted distribution functions in an idealized environment in which the initial temperatures, densities, and charge states are accurately known.
Collapse
Affiliation(s)
- Tucker Sprenkle
- Department of Physics and Astronomy, Brigham Young University, Provo, Utah 84602, USA
| | - Adam Dodson
- Department of Physics and Astronomy, Brigham Young University, Provo, Utah 84602, USA
| | - Quinton McKnight
- Department of Physics and Astronomy, Brigham Young University, Provo, Utah 84602, USA
| | - Ross Spencer
- Department of Physics and Astronomy, Brigham Young University, Provo, Utah 84602, USA
| | - Scott Bergeson
- Department of Physics and Astronomy, Brigham Young University, Provo, Utah 84602, USA
| | - Abdourahmane Diaw
- Computational Physics and Methods Group, Los Alamos National Laboratory, Los Alamos, New Mexico 87544, USA
| | - Michael S Murillo
- Department of Computational Mathematics, Science and Engineering, Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|
6
|
Ma Q, Dai J, Kang D, Murillo MS, Hou Y, Zhao Z, Yuan J. Extremely Low Electron-ion Temperature Relaxation Rates in Warm Dense Hydrogen: Interplay between Quantum Electrons and Coupled Ions. PHYSICAL REVIEW LETTERS 2019; 122:015001. [PMID: 31012692 DOI: 10.1103/physrevlett.122.015001] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/02/2017] [Revised: 10/24/2018] [Indexed: 06/09/2023]
Abstract
Theoretical and computational modeling of nonequilibrium processes in warm dense matter represents a significant challenge. The electron-ion relaxation process in warm dense hydrogen is investigated here by nonequilibrium molecular dynamics using the constrained electron force field (CEFF) method. CEFF evolves wave packets that incorporate dynamic quantum diffraction that obviates the Coulomb catastrophe. Predictions from this model reveal temperature relaxation times as much as three times longer than prior molecular dynamics results based on quantum statistical potentials. Through analyses of energy distributions and mean free paths, this result can be traced to delocalization. Finally, an improved GMS [Gericke, Murillo, and Schlanges, Phys. Rev. E 78, 025401 (2008)PRESCM1539-375510.1103/PhysRevE.78.025401] model is proposed, in which the Coulomb logarithms are in good agreement with CEFF results.
Collapse
Affiliation(s)
- Qian Ma
- Department of Physics, National University of Defense Technology, Changsha, Hunan 410073, P. R. China
| | - Jiayu Dai
- Department of Physics, National University of Defense Technology, Changsha, Hunan 410073, P. R. China
| | - Dongdong Kang
- Department of Physics, National University of Defense Technology, Changsha, Hunan 410073, P. R. China
| | - M S Murillo
- Department of Computational Mathematics, Science and Engineering, Michigan State University, East Lansing, Michigan 48824, USA
| | - Yong Hou
- Department of Physics, National University of Defense Technology, Changsha, Hunan 410073, P. R. China
| | - Zengxiu Zhao
- Department of Physics, National University of Defense Technology, Changsha, Hunan 410073, P. R. China
| | - Jianmin Yuan
- Department of Physics, National University of Defense Technology, Changsha, Hunan 410073, P. R. China
- Graduate School of China Academy of Engineering Physics, Beijing 100193, P. R. China
| |
Collapse
|
7
|
Moldabekov ZA, Groth S, Dornheim T, Kählert H, Bonitz M, Ramazanov TS. Structural characteristics of strongly coupled ions in a dense quantum plasma. Phys Rev E 2018; 98:023207. [PMID: 30253556 DOI: 10.1103/physreve.98.023207] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2018] [Indexed: 06/08/2023]
Abstract
The structural properties of strongly coupled ions in dense plasmas with moderately to strongly degenerate electrons are investigated in the framework of the one-component plasma model of ions interacting through a screened pair interaction potential. Special focus is put on the description of the electronic screening in the Singwi-Tosi-Land-Sjölander (STLS) approximation. Different cross-checks and analyses using ion potentials obtained from ground-state quantum Monte Carlo data, the random phase approximation (RPA), and existing analytical models are presented for the computation of the structural properties, such as the pair distribution and the static structure factor, of strongly coupled ions. The results are highly sensitive to the features of the screened pair interaction potential. This effect is particularly visible in the static structure factor. The applicability range of the screened potential computed from STLS is identified in terms of density and temperature of the electrons. It is demonstrated that at r_{s}>1, where r_{s} is the ratio of the mean interelectronic distance to the Bohr radius, electronic correlations beyond RPA have a nonnegligible effect on the structural properties. Additionally, the applicability of the hypernetted chain approximation for the calculation of the structural properties using the screened pair interaction potential is analyzed employing the effective coupling parameter approach.
Collapse
Affiliation(s)
- Zh A Moldabekov
- Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, Leibnizstraße 15, 24098 Kiel, Germany
- Institute for Experimental and Theoretical Physics, Al-Farabi Kazakh National University, 71 Al-Farabi str., 050040 Almaty, Kazakhstan
| | - S Groth
- Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, Leibnizstraße 15, 24098 Kiel, Germany
| | - T Dornheim
- Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, Leibnizstraße 15, 24098 Kiel, Germany
| | - H Kählert
- Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, Leibnizstraße 15, 24098 Kiel, Germany
| | - M Bonitz
- Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, Leibnizstraße 15, 24098 Kiel, Germany
| | - T S Ramazanov
- Institute for Experimental and Theoretical Physics, Al-Farabi Kazakh National University, 71 Al-Farabi str., 050040 Almaty, Kazakhstan
| |
Collapse
|
8
|
Scullard CR, Serna S, Benedict LX, Ellison CL, Graziani FR. Analytic expressions for electron-ion temperature equilibration rates from the Lenard-Balescu equation. Phys Rev E 2018; 97:013205. [PMID: 29448423 DOI: 10.1103/physreve.97.013205] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2017] [Indexed: 06/08/2023]
Abstract
In this work, we elucidate the mathematical structure of the integral that arises when computing the electron-ion temperature equilibration time for a homogeneous weakly coupled plasma from the Lenard-Balescu equation. With some minor approximations, we derive an analytic formula, requiring no input Coulomb logarithm, for the equilibration rate that is valid for moderate electron-ion temperature ratios and arbitrary electron degeneracy. For large temperature ratios, we derive the necessary correction to account for the coupled-mode effect, which can be evaluated very efficiently using ordinary Gaussian quadrature.
Collapse
Affiliation(s)
| | - Susana Serna
- Departament de Matematiques, Universitat Autònoma de Barcelona, 08193 Bellaterra-Barcelona, Spain
| | - Lorin X Benedict
- Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - C Leland Ellison
- Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - Frank R Graziani
- Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| |
Collapse
|
9
|
Groth S, Dornheim T, Bonitz M. Configuration path integral Monte Carlo approach to the static density response of the warm dense electron gas. J Chem Phys 2017; 147:164108. [DOI: 10.1063/1.4999907] [Citation(s) in RCA: 39] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Affiliation(s)
- Simon Groth
- Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, D-24098 Kiel, Germany
| | - Tobias Dornheim
- Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, D-24098 Kiel, Germany
| | - Michael Bonitz
- Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, D-24098 Kiel, Germany
| |
Collapse
|
10
|
Dornheim T, Groth S, Vorberger J, Bonitz M. Permutation-blocking path-integral Monte Carlo approach to the static density response of the warm dense electron gas. Phys Rev E 2017; 96:023203. [PMID: 28950530 DOI: 10.1103/physreve.96.023203] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2017] [Indexed: 06/07/2023]
Abstract
The static density response of the uniform electron gas is of fundamental importance for numerous applications. Here we employ the recently developed ab initio permutation blocking path integral Monte Carlo (PB-PIMC) technique [T. Dornheim et al., New J. Phys. 17, 073017 (2015)10.1088/1367-2630/17/7/073017] to carry out extensive simulations of the harmonically perturbed electron gas at warm dense matter conditions. In particular, we investigate in detail the validity of linear response theory and demonstrate that PB-PIMC allows us to obtain highly accurate results for the static density response function and, thus, the static local field correction. A comparison with dielectric approximations to our new ab initio data reveals the need for an exact treatment of correlations. Finally, we consider a superposition of multiple perturbations and discuss the implications for the calculation of the static response function.
Collapse
Affiliation(s)
- Tobias Dornheim
- Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, D-24098 Kiel, Germany
| | - Simon Groth
- Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, D-24098 Kiel, Germany
| | - Jan Vorberger
- Helmholtz-Zentrum Dresden-Rossendorf, D-01328 Dresden, Germany
| | - Michael Bonitz
- Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, D-24098 Kiel, Germany
| |
Collapse
|