skip to main content

Title: Monte Carlo simulations of ionization potential depression in dense plasmas

A particle-particle grand canonical Monte Carlo model with Coulomb pair potential interaction was used to simulate modification of ionization potentials by electrostatic microfields. The Barnes-Hut tree algorithm [J. Barnes and P. Hut, Nature 324, 446 (1986)] was used to speed up calculations of electric potential. Atomic levels were approximated to be independent of the microfields as was assumed in the original paper by Ecker and Kröll [Phys. Fluids 6, 62 (1963)]; however, the available levels were limited by the corresponding mean inter-particle distance. The code was tested on hydrogen and dense aluminum plasmas. The amount of depression was up to 50% higher in the Debye-Hückel regime for hydrogen plasmas, in the high density limit, reasonable agreement was found with the Ecker-Kröll model for hydrogen plasmas and with the Stewart-Pyatt model [J. Stewart and K. Pyatt, Jr., Astrophys. J. 144, 1203 (1966)] for aluminum plasmas. Our 3D code is an improvement over the spherically symmetric simplifications of the Ecker-Kröll and Stewart-Pyatt models and is also not limited to high atomic numbers as is the underlying Thomas-Fermi model used in the Stewart-Pyatt model.
Authors:
 [1]
  1. Department of Radiation and Chemical Physics, Institute of Physics ASCR, Na Slovance 2, 182 21 Prague 8 (Czech Republic)
Publication Date:
OSTI Identifier:
22493830
Resource Type:
Journal Article
Resource Relation:
Journal Name: Physics of Plasmas; Journal Volume: 23; Journal Issue: 1; Other Information: (c) 2016 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
70 PLASMA PHYSICS AND FUSION TECHNOLOGY; ALGORITHMS; ALUMINIUM; COMPUTERIZED SIMULATION; COULOMB FIELD; ELECTRIC POTENTIAL; HYDROGEN; MONTE CARLO METHOD; PARTICLE INTERACTIONS; PLASMA; PLASMA DENSITY; SPHERICAL CONFIGURATION; SYMMETRY; THOMAS-FERMI MODEL