1
|
Lan YS, Gu YJ, Li ZG, Li GJ, Liu L, Wang ZQ, Chen QF, Chen XR. Transport properties of a quasisymmetric binary nitrogen-oxygen mixture in the warm dense regime. Phys Rev E 2022; 105:015201. [PMID: 35193253 DOI: 10.1103/physreve.105.015201] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2021] [Accepted: 12/14/2021] [Indexed: 11/07/2022]
Abstract
Transport properties of mixtures in the warm dense matter (WDM) regime play an important role in natural astrophysics. However, a physical understanding of ionic transport properties in quasisymmetric liquid mixtures has remained elusive. Here, we present extensive ab initio molecular dynamics (AIMD) simulations on the ionic diffusion and viscosity of a quasisymmetric binary nitrogen-oxygen (N-O) mixture in a wide warm dense regime of 8-120 kK and 4.5-8.0 g/cm^{3}. Diffusion and viscosity of N-O mixtures with different compositions are obtained by using the Green-Kubo formula. Unlike asymmetric mixtures, the change of proportions in N-O mixtures slightly affects the viscosity and diffusion in the strong-coupling region. Furthermore, the AIMD results are used to build and verify a global pseudo-ion in jellium (PIJ) model for ionic transport calculations. The PIJ model succeeds in reproducing the transport properties of N-O mixtures where ionization has occurred, and provides a promising alternative approach to obtaining comparable results to AIMD simulations with relatively small computational costs. Our current results highlight the characteristic features of the quasisymmetric binary mixtures and demonstrate the applicability of the PIJ model in the WDM regime.
Collapse
Affiliation(s)
- Yang-Shun Lan
- College of Physics, Sichuan University, Chengdu 610064, 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
| | - 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
| | - 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
| | - Guo-Jun Li
- College of Physics, Sichuan University, Chengdu 610064, 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
| | - Lei Liu
- School of Science, Southwest University of Science and Technology, Mianyang 621010, China
| | - Zhao-Qi Wang
- College of Science, Xi'an University of Science and Technology, Xi'an 710054, 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 610064, People's Republic of China
| |
Collapse
|
2
|
Scheiner B, Baalrud SD. Testing thermal conductivity models with equilibrium molecular dynamics simulations of the one-component plasma. Phys Rev E 2019; 100:043206. [PMID: 31770988 DOI: 10.1103/physreve.100.043206] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2019] [Indexed: 11/07/2022]
Abstract
Equilibrium molecular dynamics simulations are used to calculate the thermal conductivity of the one-component plasma via the Green-Kubo formalism over a broad range of Coulomb coupling strength, 0.1≤Γ≤180. These simulations address previous discrepancies between computations using equilibrium versus nonequilibrium methods. Analysis of heat flux autocorrelation functions show that very long (6×10^{5}ω_{p}^{-1}) time series are needed to reduce the noise level to allow ≲2% accuracy. The simulations provide accurate data for Γ≲1. This enables a test of the traditional Landau-Spitzer theory, which is found to agree with the simulations for Γ≲0.3. It also enables tests of theories to address moderate and strong Coulomb coupling. Two are found to provide accurate extensions to the moderate coupling regime of Γ≲10, but none are accurate in the Γ≳10 regime where potential energy transport and coupling between mass flow and stress dominate thermal conduction.
Collapse
Affiliation(s)
- Brett Scheiner
- Department of Physics and Astronomy, University of Iowa, Iowa City, Iowa 52240, USA
| | - Scott D Baalrud
- Department of Physics and Astronomy, University of Iowa, Iowa City, Iowa 52240, USA
| |
Collapse
|
3
|
Liu L, Li ZG, Dai JY, Chen QF, Chen XR. Quantum molecular dynamics study on the proton exchange, ionic structures, and transport properties of warm dense hydrogen-deuterium mixtures. Phys Rev E 2018; 97:063204. [PMID: 30011461 DOI: 10.1103/physreve.97.063204] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2018] [Indexed: 06/08/2023]
Abstract
Comprehensive knowledge of physical properties such as equation of state (EOS), proton exchange, dynamic structures, diffusion coefficients, and viscosities of hydrogen-deuterium mixtures with densities from 0.1 to 5 g/cm^{3} and temperatures from 1 to 50 kK has been presented via quantum molecular dynamics (QMD) simulations. The existing multi-shock experimental EOS provides an important benchmark to evaluate exchange-correlation functionals. The comparison of simulations with experiments indicates that a nonlocal van der Waals density functional (vdW-DF1) produces excellent results. Fraction analysis of molecules using a weighted integral over pair distribution functions was performed. A dissociation diagram together with a boundary where the proton exchange (H_{2}+D_{2}⇌2HD) occurs was generated, which shows evidence that the HD molecules form as the H_{2} and D_{2} molecules are almost 50% dissociated. The mechanism of proton exchange can be interpreted as a process of dissociation followed by recombination. The ionic structures at extreme conditions were analyzed by the effective coordination number model. High-order cluster, circle, and chain structures can be founded in the strongly coupled warm dense regime. The present QMD diffusion coefficient and viscosity can be used to benchmark two analytical one-component plasma (OCP) models: the Coulomb and Yukawa OCP models.
Collapse
Affiliation(s)
- Lei Liu
- Institute of Atomic and Molecular Physics, College of Physical Science and Technology, Sichuan University, Chengdu 610064, People's Republic of China
- National Key Laboratory for Shock Wave and Detonation Physics Research, Institute of Fluid Physics, China 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, China Academy of Engineering Physics, Mianyang 621900, People's Republic of China
| | - Jia-Yu Dai
- 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, China Academy of Engineering Physics, Mianyang 621900, People's Republic of China
| | - Xiang-Rong Chen
- Institute of Atomic and Molecular Physics, College of Physical Science and Technology, Sichuan University, Chengdu 610064, People's Republic of China
| |
Collapse
|
4
|
Reinholz H, Röpke G, Rosmej S, Redmer R. Conductivity of warm dense matter including electron-electron collisions. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 91:043105. [PMID: 25974600 DOI: 10.1103/physreve.91.043105] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/16/2015] [Indexed: 06/04/2023]
Abstract
We present an approach that can resolve the controversy with respect to the role of electron-electron collisions in calculating the dynamic conductivity of dense plasmas. In particular, the dc conductivity is analyzed in the low-density, nondegenerate limit where the Spitzer theory is valid and electron-electron collisions lead to the well-known reduction in comparison to the result considering only electron-ion collisions (Lorentz model). With increasing degeneracy, the contribution of electron-electron collisions to the dc conductivity is decreasing and can be neglected for the liquid metal domain where the Ziman theory is applicable. We give expressions for the effect of electron-electron collisions in calculating the conductivity in the warm dense matter region, i.e., for strongly coupled Coulomb systems at arbitrary degeneracy.
Collapse
Affiliation(s)
- H Reinholz
- Universität Rostock, Institut für Physik, 18051 Rostock, Germany and University of Western Australia School of Physics, WA 6009 Crawley, Australia
| | - G Röpke
- Universität Rostock, Institut für Physik, 18051 Rostock, Germany
| | - S Rosmej
- Universität Rostock, Institut für Physik, 18051 Rostock, Germany
| | - R Redmer
- Universität Rostock, Institut für Physik, 18051 Rostock, Germany
| |
Collapse
|
5
|
Hu SX, Collins LA, Goncharov VN, Boehly TR, Epstein R, McCrory RL, Skupsky S. First-principles opacity table of warm dense deuterium for inertial-confinement-fusion applications. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 90:033111. [PMID: 25314551 DOI: 10.1103/physreve.90.033111] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/25/2014] [Indexed: 06/04/2023]
Abstract
Accurate knowledge of the optical properties of a warm dense deuterium-tritium (DT) mixture is important for reliable design of inertial confinement fusion (ICF) implosions using radiation-hydrodynamics simulations. The opacity of a warm dense DT shell essentially determines how much radiation from hot coronal plasmas can be deposited in the DT fuel of an imploding capsule. Even for the simplest species of hydrogen, the accurate calculation of their opacities remains a challenge in the warm-dense matter regime because strong-coupling and quantum effects play an important role in such plasmas. With quantum-molecular-dynamics (QMD) simulations, we have derived a first-principles opacity table (FPOT) of deuterium (and the DT mixture by mass scaling) for a wide range of densities from ρ(D)=0.5 to 673.518g/cm(3) and temperatures from T=5000K up to the Fermi temperature T(F) for each density. Compared with results from the astrophysics opacity table (AOT) currently used in our hydrocodes, the FPOT of deuterium from our QMD calculations has shown a significant increase in opacity for strongly coupled and degenerate plasma conditions by a factor of 3-100 in the ICF-relevant photon-energy range. As conditions approach those of classical plasma, the opacity from the FPOT converges to the corresponding values of the AOT. By implementing the FPOT of deuterium and the DT mixture into our hydrocodes, we have performed radiation-hydrodynamics simulations for low-adiabat cryogenic DT implosions on the OMEGA laser and for direct-drive-ignition designs for the National Ignition Facility. The simulation results using the FPOT show that the target performance (in terms of neutron yield and energy gain) could vary from ∼10% up to a factor of ∼2 depending on the adiabat of the imploding DT capsule; the lower the adiabat, the more variation is seen in the prediction of target performance when compared to the AOT modeling.
Collapse
Affiliation(s)
- S X Hu
- Laboratory for Laser Energetics, University of Rochester, 250 East River Road, Rochester, New York 14623-1299, USA
| | - L A Collins
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - V N Goncharov
- Laboratory for Laser Energetics, University of Rochester, 250 East River Road, Rochester, New York 14623-1299, USA
| | - T R Boehly
- Laboratory for Laser Energetics, University of Rochester, 250 East River Road, Rochester, New York 14623-1299, USA
| | - R Epstein
- Laboratory for Laser Energetics, University of Rochester, 250 East River Road, Rochester, New York 14623-1299, USA
| | - R L McCrory
- Laboratory for Laser Energetics, University of Rochester, 250 East River Road, Rochester, New York 14623-1299, USA
| | - S Skupsky
- Laboratory for Laser Energetics, University of Rochester, 250 East River Road, Rochester, New York 14623-1299, USA
| |
Collapse
|
6
|
Hu SX, Collins LA, Boehly TR, Kress JD, Goncharov VN, Skupsky S. First-principles thermal conductivity of warm-dense deuterium plasmas for inertial confinement fusion applications. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 89:043105. [PMID: 24827353 DOI: 10.1103/physreve.89.043105] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/11/2013] [Indexed: 06/03/2023]
Abstract
Thermal conductivity (κ) of both the ablator materials and deuterium-tritium (DT) fuel plays an important role in understanding and designing inertial confinement fusion (ICF) implosions. The extensively used Spitzer model for thermal conduction in ideal plasmas breaks down for high-density, low-temperature shells that are compressed by shocks and spherical convergence in imploding targets. A variety of thermal-conductivity models have been proposed for ICF hydrodynamic simulations of such coupled and degenerate plasmas. The accuracy of these κ models for DT plasmas has recently been tested against first-principles calculations using the quantum molecular-dynamics (QMD) method; although mainly for high densities (ρ > 100 g/cm3), large discrepancies in κ have been identified for the peak-compression conditions in ICF. To cover the wide range of density-temperature conditions undergone by ICF imploding fuel shells, we have performed QMD calculations of κ for a variety of deuterium densities of ρ = 1.0 to 673.518 g/cm3, at temperatures varying from T = 5 × 103 K to T = 8 × 106 K. The resulting κQMD of deuterium is fitted with a polynomial function of the coupling and degeneracy parameters Γ and θ, which can then be used in hydrodynamic simulation codes. Compared with the "hybrid" Spitzer-Lee-More model currently adopted in our hydrocode lilac, the hydrosimulations using the fitted κQMD have shown up to ∼20% variations in predicting target performance for different ICF implosions on OMEGA and direct-drive-ignition designs for the National Ignition Facility (NIF). The lower the adiabat of an imploding shell, the more variations in predicting target performance using κQMD. Moreover, the use of κQMD also modifies the shock conditions and the density-temperature profiles of the imploding shell at early implosion stage, which predominantly affects the final target performance. This is in contrast to the previous speculation that κQMD changes mainly the inside ablation process during the hot-spot formation of an ICF implosion.
Collapse
Affiliation(s)
- S X Hu
- Laboratory for Laser Energetics, University of Rochester, Rochester, New York 14623, USA
| | - L A Collins
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - T R Boehly
- Laboratory for Laser Energetics, University of Rochester, Rochester, New York 14623, USA
| | - J D Kress
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - V N Goncharov
- Laboratory for Laser Energetics, University of Rochester, Rochester, New York 14623, USA
| | - S Skupsky
- Laboratory for Laser Energetics, University of Rochester, Rochester, New York 14623, USA
| |
Collapse
|
7
|
Sun H, Kang D, Dai J, Zeng J, Yuan J. Quantum molecular dynamics study on the structures and dc conductivity of warm dense silane. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 89:022128. [PMID: 25353443 DOI: 10.1103/physreve.89.022128] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/14/2013] [Indexed: 06/04/2023]
Abstract
The ionic and electronic structures of warm dense silane at the densities of 1.795, 2.260, 3.382, and 3.844 g/cm(3) have been studied with temperatures from 1000 K to 3 eV using quantum molecular dynamics simulations. At all densities, the structures are melted above 1000 K. The matter states are characterized as polymeric from 1000 to 4000 K and become dense plasma states with further increasing temperature to 1 eV. At two lower densities of 1.795 and 2.260 g/cm(3), silane first dissociates and then becomes the polymeric state via a chain state from the initial crystalline structure. At higher densities, however, no dissociation stage was found. These findings can help us understand how the warm dense matter forms. A rise is found for the direct current electric conductivity at T ∼ 1000 K, indicating the nonmetal-to-metal transition. The conductivity decreases slightly with the increase of temperature, which is due to the more disordered structures at higher temperatures.
Collapse
Affiliation(s)
- Huayang Sun
- Department of Physics, College of Science, National University of Defense Technology, Changsha 410073, People's Republic of China
| | - Dongdong Kang
- Department of Physics, College of Science, National University of Defense Technology, Changsha 410073, People's Republic of China
| | - Jiayu Dai
- Department of Physics, College of Science, National University of Defense Technology, Changsha 410073, People's Republic of China
| | - Jiaolong Zeng
- Department of Physics, College of Science, National University of Defense Technology, Changsha 410073, People's Republic of China
| | - Jianmin Yuan
- Department of Physics, College of Science, National University of Defense Technology, Changsha 410073, People's Republic of China and State Key Laboratory of High Performance Computing, National University of Defense Technology, Changsha 410073, People's Republic of China
| |
Collapse
|