1
|
Vorberger J, Chapman DA. Quantum theory for the dynamic structure factor in correlated two-component systems in nonequilibrium: Application to x-ray scattering. Phys Rev E 2018; 97:013203. [PMID: 29448372 DOI: 10.1103/physreve.97.013203] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2017] [Indexed: 06/08/2023]
Abstract
We present a quantum theory for the dynamic structure factors in nonequilibrium, correlated, two-component systems such as plasmas or warm dense matter. The polarization function, which is needed as the input for the calculation of the structure factors, is calculated in nonequilibrium based on a perturbation expansion in the interaction strength. To make our theory applicable for x-ray scattering, a generalized Chihara decomposition for the total electron structure factor in nonequilibrium is derived. Examples are given and the influence of correlations and exchange on the structure and the x-ray-scattering spectrum are discussed for a model nonequilibrium distribution, as often encountered during laser heating of materials, as well as for two-temperature systems.
Collapse
Affiliation(s)
- J Vorberger
- Institute of Radiation Physics, Helmholtz-Zentrum Dresden-Rossendorf e.V., 01328 Dresden, Germany
| | - D A Chapman
- AWE plc, Aldermaston, Reading RG7 4PR, United Kingdom
- Centre for Fusion, Space and Astrophysics, University of Warwick, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
2
|
Dornheim T, Groth S, Sjostrom T, Malone FD, Foulkes WMC, Bonitz M. Ab Initio Quantum Monte Carlo Simulation of the Warm Dense Electron Gas in the Thermodynamic Limit. PHYSICAL REVIEW LETTERS 2016; 117:156403. [PMID: 27768371 DOI: 10.1103/physrevlett.117.156403] [Citation(s) in RCA: 57] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/28/2016] [Indexed: 06/06/2023]
Abstract
We perform ab initio quantum Monte Carlo (QMC) simulations of the warm dense uniform electron gas in the thermodynamic limit. By combining QMC data with the linear response theory, we are able to remove finite-size errors from the potential energy over the substantial parts of the warm dense regime, overcoming the deficiencies of the existing finite-size corrections by Brown et al. [Phys. Rev. Lett. 110, 146405 (2013)]. Extensive new QMC results for up to N=1000 electrons enable us to compute the potential energy V and the exchange-correlation free energy F_{xc} of the macroscopic electron gas with an unprecedented accuracy of |ΔV|/|V|,|ΔF_{xc}|/|F|_{xc}∼10^{-3}. A comparison of our new data to the recent parametrization of F_{xc} by Karasiev et al. [Phys. Rev. Lett. 112, 076403 (2014)] reveals significant deviations to the latter.
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
| | - Travis Sjostrom
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Fionn D Malone
- Department of Physics, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - W M C Foulkes
- Department of Physics, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Michael Bonitz
- Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, D-24098 Kiel, Germany
| |
Collapse
|
3
|
Malone FD, Blunt NS, Brown EW, Lee DKK, Spencer JS, Foulkes WMC, Shepherd JJ. Accurate Exchange-Correlation Energies for the Warm Dense Electron Gas. PHYSICAL REVIEW LETTERS 2016; 117:115701. [PMID: 27661699 DOI: 10.1103/physrevlett.117.115701] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/23/2016] [Indexed: 06/06/2023]
Abstract
The density matrix quantum Monte Carlo (DMQMC) method is used to sample exact-on-average N-body density matrices for uniform electron gas systems of up to 10^{124} matrix elements via a stochastic solution of the Bloch equation. The results of these calculations resolve a current debate over the accuracy of the data used to parametrize finite-temperature density functionals. Exchange-correlation energies calculated using the real-space restricted path-integral formalism and the k-space configuration path-integral formalism disagree by up to ∼10% at certain reduced temperatures T/T_{F}≤0.5 and densities r_{s}≤1. Our calculations confirm the accuracy of the configuration path-integral Monte Carlo results available at high density and bridge the gap to lower densities, providing trustworthy data in the regime typical of planetary interiors and solids subject to laser irradiation. We demonstrate that the DMQMC method can calculate free energies directly and present exact free energies for T/T_{F}≥1 and r_{s}≤2.
Collapse
Affiliation(s)
- Fionn D Malone
- Department of Physics, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - N S Blunt
- University Chemical Laboratory, Cambridge University, Lensfield Road, Cambridge CB2 1EW, United Kingdom
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, Stuttgart 70569, Germany
| | - Ethan W Brown
- Institute for Theoretical Physics, ETH Zürich, Wolfgang Pauli Strasse 27, 8093 Zürich, Switzerland
| | - D K K Lee
- Department of Physics, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - J S Spencer
- Department of Physics, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - W M C Foulkes
- Department of Physics, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - James J Shepherd
- Department of Physics, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| |
Collapse
|
4
|
Schoof T, Groth S, Vorberger J, Bonitz M. Ab Initio Thermodynamic Results for the Degenerate Electron Gas at Finite Temperature. PHYSICAL REVIEW LETTERS 2015; 115:130402. [PMID: 26451539 DOI: 10.1103/physrevlett.115.130402] [Citation(s) in RCA: 50] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/17/2015] [Indexed: 06/05/2023]
Abstract
The uniform electron gas at finite temperature is of key relevance for many applications in dense plasmas, warm dense matter, laser excited solids, and much more. Accurate thermodynamic data for the uniform electron gas are an essential ingredient for many-body theories, in particular, density-functional theory. Recently, first-principles restricted path integral Monte Carlo results became available, which, however, had to be restricted to moderate degeneracy, i.e., low to moderate densities with r_{s}=r[over ¯]/a_{B}≳1. Here we present novel first-principles configuration path integral Monte Carlo results for electrons for r_{s}≤4. We also present quantum statistical data within the e^{4} approximation that are in good agreement with the simulations at small to moderate r_{s}.
Collapse
Affiliation(s)
- T Schoof
- Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, D-24098 Kiel, Germany
| | - S Groth
- Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, D-24098 Kiel, Germany
| | - J Vorberger
- Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, D-24098 Kiel, Germany
| | - M Bonitz
- Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, D-24098 Kiel, Germany
| |
Collapse
|
5
|
Ramazanov TS, Moldabekov ZA, Gabdullin MT. Effective potentials of interactions and thermodynamic properties of a nonideal two-temperature dense plasma. Phys Rev E 2015; 92:023104. [PMID: 26382532 DOI: 10.1103/physreve.92.023104] [Citation(s) in RCA: 45] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2015] [Indexed: 11/07/2022]
Abstract
In this article a dense nonideal, nonisothermal plasma is considered. New effective screened interaction potentials taking into account quantum-mechanical diffraction and symmetry effects have been obtained. The effective potential of the ion-ion interaction in plasmas with a strongly coupled ion subsystem and semiclassical electron subsystem is presented. Based on the obtained effective potentials the analytical expressions for internal energy and the pressure of the considered plasma were obtained.
Collapse
Affiliation(s)
- T S Ramazanov
- Al-Farabi Kazakh National University, Institute of Experimental and Theoretical Physics, Almaty 050040, Kazakhstan
| | - Zh A Moldabekov
- Al-Farabi Kazakh National University, Institute of Experimental and Theoretical Physics, Almaty 050040, Kazakhstan
| | - M T Gabdullin
- Al-Farabi Kazakh National University, Institute of Experimental and Theoretical Physics, Almaty 050040, Kazakhstan
| |
Collapse
|
6
|
Omarbakiyeva YA, Reinholz H, Röpke G. Cluster virial expansion of the equation of state for hydrogen plasma with e-H(2) contributions. Phys Rev E 2015; 91:043103. [PMID: 25974598 DOI: 10.1103/physreve.91.043103] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2014] [Indexed: 11/07/2022]
Abstract
The equation of state of partially ionized hydrogen plasma is considered with special focus on the contribution of the e-H(2) interaction. Traditional semiempirical concepts such as the excluded volume are improved using microscopic approaches to treat the e-H(2) problem. Within a cluster virial expansion, the Beth-Uhlenbeck formula is applied to infer the contribution of bound and scattering states to the temperature-dependent second virial coefficient. The scattering states are calculated using the phase expansion method for the polarization interaction that incorporates experimental data for the e-H(2) scattering cross section. We present results for the scattering phase shifts, differential scattering cross sections, and the second virial coefficient due to the e-H(2) interaction. The influence of this interaction on the composition of the partially ionized hydrogen plasma is confined to the parameter range where both the H(2) and the free-electron components are abundant.
Collapse
Affiliation(s)
- Y A Omarbakiyeva
- Institute of Physics, University of Rostock, D-18051 Rostock, Germany.,International IT University, 050040 Almaty, Kazakhstan
| | - H Reinholz
- Institute of Physics, University of Rostock, D-18051 Rostock, Germany.,The University of Western Australia, Crawley, Western Australia 6009, Australia
| | - G Röpke
- Institute of Physics, University of Rostock, D-18051 Rostock, Germany
| |
Collapse
|
7
|
Fortmann C, Lee HJ, Döppner T, Falcone RW, Kritcher AL, Landen OL, Glenzer SH. Measurement of the adiabatic index in be compressed by counterpropagating shocks. PHYSICAL REVIEW LETTERS 2012; 108:175006. [PMID: 22680877 DOI: 10.1103/physrevlett.108.175006] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/01/2011] [Indexed: 06/01/2023]
Abstract
We report on the first direct measurement of the adiabatic index γ through x-ray Thomson scattering from shock-compressed beryllium. 9 keV x-ray photons probe the bulk properties of matter during the collision of two counterpropagating shocks. This novel experimental technique determines γ by using only the measured mass densities and vanishing particle velocity at the point of shock collision to close the Rankine-Hugoniot equations. We find γ>5/3 at 3× compression, clearly different from ideal gas behavior. At 6× compression, γ shows the convergence to the ideal gas limit, in agreement with linear scaling laws.
Collapse
Affiliation(s)
- C Fortmann
- Lawrence Livermore National Laboratory, Livermore, California 94551, USA
| | | | | | | | | | | | | |
Collapse
|
8
|
Vorberger J, Gericke DO, Bornath T, Schlanges M. Energy and temperature relaxation described by nonequilibrium green's functions. ACTA ACUST UNITED AC 2010. [DOI: 10.1088/1742-6596/220/1/012002] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
9
|
Vorberger J, Gericke DO, Bornath T, Schlanges M. Energy relaxation in dense, strongly coupled two-temperature plasmas. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2010; 81:046404. [PMID: 20481844 DOI: 10.1103/physreve.81.046404] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/25/2009] [Revised: 03/01/2010] [Indexed: 05/29/2023]
Abstract
A quantum kinetic approach for the energy relaxation in strongly coupled plasmas with different electron and ion temperatures is presented. Based on the density operator formalism, we derive a balance equation for the energies of electrons and ions connecting kinetic, correlation, and exchange energies with a quite general expression for the electron-ion energy-transfer rate. The latter is given in terms of the correlation function of density fluctuations which allows for a derivation of increasingly realistic approximation schemes including a coupled-mode expression. The equilibration of the contributions of the total energy including the species temperatures in dense hydrogen and beryllium relevant for inertial confinement fusion is investigated as an example.
Collapse
Affiliation(s)
- J Vorberger
- Department of Physics, University of Warwick, Centre for Fusion, Space and Astrophysics, Coventry CV4 7Al, United Kingdom
| | | | | | | |
Collapse
|
10
|
Fortmann C. Single-particle spectral function for the classical one-component plasma. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2009; 79:016404. [PMID: 19257145 DOI: 10.1103/physreve.79.016404] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/15/2008] [Indexed: 05/27/2023]
Abstract
The spectral function for an electron one-component plasma is calculated self-consistently using the GW;{(0)} approximation for the single-particle self-energy. In this way, correlation effects that go beyond the mean-field description of the plasma are contained, i.e., the collisional damping of single-particle states, the dynamical screening of the interaction, and the appearance of collective plasma modes. Second, a nonperturbative analytic solution for the on-shell GW;{(0)} self-energy as a function of momentum is presented. It reproduces the numerical data for the spectral function with a relative error of less than 10% in the regime where the Debye screening parameter is smaller than the inverse Bohr radius, kappa<1a_{B};{-1} . In the limit of low density, the nonperturbative self-energy behaves as n;{14} , whereas a perturbation expansion leads to the unphysical result of a density-independent self-energy [Fennel and Wilfer, Ann. Phys. (Leipzig) 32, 265 (1974)]. The derived expression will greatly facilitate the calculation of observables in correlated plasmas (transport properties, equation of state) that need the spectral function as an input quantity. This is demonstrated for the shift of the chemical potential, which is computed from the analytical formulas and compared to the GW;{(0)} result. At a plasma temperature of 100eV and densities below 10;{21}cm;{-3} , the two approaches deviate by less than 10% from each other.
Collapse
Affiliation(s)
- C Fortmann
- Institut für Physik, Universität Rostock, 18051 Rostock, Germany.
| |
Collapse
|
11
|
Wünsch K, Hilse P, Schlanges M, Gericke DO. Structure of strongly coupled multicomponent plasmas. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2008; 77:056404. [PMID: 18643173 DOI: 10.1103/physreve.77.056404] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/30/2007] [Revised: 03/02/2008] [Indexed: 05/26/2023]
Abstract
We investigate the short-range structure in strongly coupled fluidlike plasmas using the hypernetted chain approach generalized to multicomponent systems. Good agreement with numerical simulations validates this method for the parameters considered. We found a strong mutual impact on the spatial arrangement for systems with multiple ion species which is most clearly pronounced in the static structure factor. Quantum pseudopotentials were used to mimic diffraction and exchange effects in dense electron-ion systems. We demonstrate that the different kinds of pseudopotentials proposed lead to large differences in both the pair distributions and structure factors. Large discrepancies were also found in the predicted ion feature of the x-ray scattering signal, illustrating the need for comparison with full quantum calculations or experimental verification.
Collapse
Affiliation(s)
- K Wünsch
- Centre for Fusion, Space and Astrophysics, Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom
| | | | | | | |
Collapse
|
12
|
Kraeft WD, Bonitz M. Thermodynamics of a correlated confined plasma II. Mesoscopic system. ACTA ACUST UNITED AC 2006. [DOI: 10.1088/1742-6596/35/1/008] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|