Stopping-power enhancement from discrete particle-wake correlations in high-energy-density plasmas.

Three-dimensional (3D) simulations of electron beams propagating in high-energy-density plasmas using the quasistatic Particle-in-Cell (PIC) code QuickPIC demonstrate a significant increase in stopping power when beam electrons mutually interact via their wakes. Each beam electron excites a plasma wave wake of wavelength ∼2πc/ω_{pe}, where c is the speed of light and ω_{pe} is the background plasma frequency. We show that a discrete collection of electrons undergoes a beam-plasma-like instability caused by mutual particle-wake interactions that causes electrons to bunch in the beam, even for beam densities n_{b} for which fluid theory breaks down. This bunching enhances the beam's stopping power, which we call "correlated stopping," and the effect increases with the "correlation number" N_{b}≡n_{b}(c/ω_{pe})^{3}. For example, a beam of monoenergetic 9.7 MeV electrons with N_{b}=1/8, in a cold background plasma with n_{e}=10^{26}cm^{-3} (450 g cm^{-3} DT), has a stopping power of 2.28±0.04 times the single-electron value, which increases to 1220±5 for N_{b}=64. The beam also experiences transverse filamentation, which eventually limits the stopping enhancement.

Neutron Time-of-Flight Measurements of Charged-Particle Energy Loss in Inertial Confinement Fusion Plasmas.

Neutron spectra from secondary ^{3}H(d,n)α reactions produced by an implosion of a deuterium-gas capsule at the National Ignition Facility have been measured with order-of-magnitude improvements in statistics and resolution over past experiments. These new data and their sensitivity to the energy loss of fast tritons emitted from thermal ^{2}H(d,p)^{3}H reactions enable the first statistically significant investigation of charged-particle stopping via the emitted neutron spectrum. Radiation-hydrodynamic simulations, constrained to match a number of observables from the implosion, were used to predict the neutron spectra while employing two different energy loss models. This analysis represents the first test of stopping models under inertial confinement fusion conditions, covering plasma temperatures of k_{B}T≈1-4 keV and particle densities of n≈(12-2)×10^{24} cm^{-3}. Under these conditions, we find significant deviations of our data from a theory employing classical collisions whereas the theory including quantum diffraction agrees with our data.

Kinetic theory molecular dynamics and hot dense matter: theoretical foundations.

Electrons are weakly coupled in hot, dense matter that is created in high-energy-density experiments. They are also mildly quantum mechanical and the ions associated with them are classical and may be strongly coupled. In addition, the dynamical evolution of plasmas under these hot, dense matter conditions involve a variety of transport and energy exchange processes. Quantum kinetic theory is an ideal tool for treating the electrons but it is not adequate for treating the ions. Molecular dynamics is perfectly suited to describe the classical, strongly coupled ions but not the electrons. We develop a method that combines a Wigner kinetic treatment of the electrons with classical molecular dynamics for the ions. We refer to this hybrid method as "kinetic theory molecular dynamics," or KTMD. The purpose of this paper is to derive KTMD from first principles and place it on a firm theoretical foundation. The framework that KTMD provides for simulating plasmas in the hot, dense regime is particularly useful since current computational methods are generally limited by their inability to treat the dynamical quantum evolution of the electronic component. Using the N-body von Neumann equation for the electron-proton plasma, three variations of KTMD are obtained. Each variant is determined by the physical state of the plasma (e.g., collisional versus collisionless). The first variant of KTMD yields a closed set of equations consisting of a mean-field quantum kinetic equation for the electron one-particle distribution function coupled to a classical Liouville equation for the protons. The latter equation includes both proton-proton Coulombic interactions and an effective electron-proton interaction that involves the convolution of the electron density with the electron-proton Coulomb potential. The mean-field approach is then extended to incorporate equilibrium electron-proton correlations through the Singwi-Tosi-Land-Sjolander (STLS) ansatz. This is the second variant of KTMD. The STLS contribution produces an effective electron-proton interaction that involves the electron-proton structure factor, thereby extending the usual mean-field theory to correlated but near equilibrium systems. Finally, a third variant of KTMD is derived. It includes dynamical electrons and their correlations coupled to a MD description for the ions. A set of coupled equations for the one-particle electron Wigner function and the electron-electron and electron-proton correlation functions are coupled to a classical Liouville equation for the protons. This latter variation has both time and momentum dependent correlations.

Molecular dynamics simulations of temperature equilibration in dense hydrogen.

The temperature equilibration rate between electrons and protons in dense hydrogen has been calculated with molecular dynamics simulations for temperatures between 10 and 600eV and densities between 10;{20}cm;{-3}to10;{24}cm;{-3} . Careful attention has been devoted to convergence of the simulations, including the role of semiclassical potentials. We find that for Coulomb logarithms L greater, similar1 , a model by Gericke-Murillo-Schlanges (GMS) [D. O. Gericke, Phys. Rev. E 65, 036418 (2002)] based on a T -matrix method and the approach by Brown-Preston-Singleton [L. S. Brown, Phys. Rep. 410, 237 (2005)] agrees with the simulation data to within the error bars of the simulation. For smaller Coulomb logarithms, the GMS model is consistent with the simulation results. Landau-Spitzer models are consistent with the simulation data for L>4 .