High-order implicit particle-in-cell method for plasma simulations at solid densities

Ion Kinetics and Neutron Generation Associated with Electromagnetic Turbulence in Laboratory-Scale Counterstreaming Plasmas

Electromagnetic turbulence and ion kinetics in counterstreaming plasmas hold great significance in laboratory astrophysics, such as turbulence field amplification and particle energization. Here, we quantitatively demonstrate for the first time how electromagnetic turbulence affects ion kinetics under achievable laboratory conditions (millimeter-scale interpenetrating plasmas with initial velocity of 2000  km/s, density of 4×10^{19}  cm^{-3}, and temperature of 100 eV) utilizing a recently developed high-order implicit particle-in-cell code without scaling transformation. It is found that the electromagnetic turbulence is driven by ion two-stream and filamentation instabilities. For the magnetized scenarios where an applied magnetic field of tens of Tesla is perpendicular to plasma flows, the growth rates of instabilities increase with the strengthening of applied magnetic field, which therefore leads to a significant enhancement of turbulence fields. Under the competition between the stochastic acceleration due to electromagnetic turbulence and collisional thermalization, ion distribution function shows a distinct super-Gaussian shape, and the ion kinetics are manifested in neutron yields and spectra. Our results have well explained the recent unmagnetized experimental observations, and the findings of magnetized scenario can be verified by current astrophysical experiments.

Large-scale kinetic simulations of colliding plasmas within a hohlraum of indirect-drive inertial confinement fusion

The National Ignition Facility has recently achieved successful burning plasma and ignition using the inertial confinement fusion (ICF) approach. However, there are still many fundamental physics phenomena that are not well understood, including the kinetic processes in the hohlraum. Shan et al. [Phys. Rev. Lett. 120, 195001 (2018)0031-900710.1103/PhysRevLett.120.195001] utilized the energy spectra of neutrons to investigate the kinetic colliding plasma in a hohlraum of indirect drive ICF. However, due to the typical large spatial-temporal scales, this experiment could not be well simulated by using available codes at that time. Utilizing our advanced high-order implicit PIC code, LAPINS, we were able to successfully reproduce the experiment on a large scale of both spatial and temporal dimensions, in which the original computational scale was increased by approximately seven to eight orders of magnitude. Not only is the validity of the explanation of the experiment confirmed by our simulations, i.e., the abnormally large width of neutron spectra comes from beam-target nuclear fusions, but also a different physical insight into the source of energetic deuterium ions is provided. The acceleration of deuterium ions can be categorized into two components: one is propelled by a sheath electric field created by the charge separation at the onset, while the other is a result of the reflection of the potential of the shock wave. The robustness of the acceleration mechanism is analyzed with varying initial conditions, e.g., temperatures, drifting velocity, and ion components. This paper might serve as a reference for benchmark simulations of upcoming simulation codes and may be relevant for future research on mixtures and entropy increments at plasma interfaces.

Dynamics effects of tritium reduction on the energy gain of D-T fuel pellet using double cone ignition

In this paper a study of the behavior of Deuterium-Tritium (D-T) plasma nuclear fusion reaction in terms of variations of time and temperature and in the presence of deuterium-tritium sources using double cone ignition is presented. The aim is the determination of the optimum physical conditions with low tritium consumption rate for obtaining the total energy gain with a value of greater than 200.

Double-cone ignition scheme for inertial confinement fusion

While major progress has been made in the research of inertial confinement fusion, significant challenges remain in the pursuit of ignition. To tackle the challenges, we propose a double-cone ignition (DCI) scheme, in which two head-on gold cones are used to confine deuterium-tritium (DT) shells imploded by high-power laser pulses. The scheme is composed of four progressive controllable processes: quasi-isentropic compression, acceleration, head-on collision and fast heating of the compressed fuel. The quasi-isentropic compression is performed inside two head-on cones. At the later stage of the compression, the DT shells in the cones are accelerated to forward velocities of hundreds of km s^-1. The head-on collision of the compressed and accelerated fuels from the cone tips transfer the forward kinetic energy to the thermal energy of the colliding fuel with an increased density. The preheated high-density fuel can keep its status for a period of approximately 200 ps. Within this period, MeV electrons generated by ps heating laser pulses, guided by a ns laser-produced strong magnetic field further heat the fuel efficiently. Our simulations show that the implosion inside the head-on cones can greatly mitigate the energy requirement for compression; the collision can preheat the compressed fuel of approximately 300 g cm^-3 to a temperature above keV. The fuel can then reach an ignition temperature of greater than 5 keV with magnetically assisted heating of MeV electrons generated by the heating laser pulses. Experimental campaigns to demonstrate the scheme have already begun. This article is part of a discussion meeting issue 'Prospects for high gain inertial fusion energy (part 1)'.

Particle-in-cell simulation method for macroscopic degenerate plasmas

Nowadays hydrodynamic equations coupled with external equation of states provided by quantum mechanical calculations is a widely used approach for simulations of macroscopic degenerate plasmas. Although such an approach is proven to be efficient and shows many good features, especially for large scale simulations, it encounters intrinsic challenges when involving kinetic effects. As a complement, here we have invented a fully kinetic numerical approach for macroscopic degenerate plasmas. This approach is based on first principle Boltzmann-Uhling-Uhlenbeck equations coupled with Maxwell's equation, and is eventually achieved via an existing particle-in-cell simulation code named LAPINS. In this approach, degenerate particles obey Fermi-Dirac statistics and nondegenerate particles follow the typical Maxwell-Boltzmann statistics. The equation of motion of both degenerate and nondegenerate particles are governed by long range collective electromagnetic fields and close particle-particle collisions. Especially, Boltzmann-Uhling-Uhlenbeck collisions ensure that evolution of degenerate particles is enforced by the Pauli exclusion principle. The code is applied to several benchmark simulations, including electronic conductivity for aluminium with varying temperatures from 2 eV to 50 eV, thermalization of alpha particles in a cold fuel shell in inertial confinement fusion, and rapid heating of solid sample by short and intense laser pulses.

Observation of a high degree of stopping for laser-accelerated intense proton beams in dense ionized matter

Intense particle beams generated from the interaction of ultrahigh intensity lasers with sample foils provide options in radiography, high-yield neutron sources, high-energy-density-matter generation, and ion fast ignition. An accurate understanding of beam transportation behavior in dense matter is crucial for all these applications. Here we report the experimental evidence on one order of magnitude enhancement of intense laser-accelerated proton beam stopping in dense ionized matter, in comparison with the current-widely used models describing individual ion stopping in matter. Supported by particle-in-cell (PIC) simulations, we attribute the enhancement to the strong decelerating electric field approaching 1 GV/m that can be created by the beam-driven return current. This collective effect plays the dominant role in the stopping of laser-accelerated intense proton beams in dense ionized matter. This finding is essential for the optimum design of ion driven fast ignition and inertial confinement fusion.

A detailed understanding of particle stopping in matter is essential for nuclear fusion and high energy density science. Here, the authors report one order of magnitude enhancement of intense laser-accelerated proton beam stopping in dense ionized matter in comparison with currently used models describing ion stopping in matter.

Uniform warm dense matter formed by direct laser heating in the presence of external magnetic fields

With the recent realization of kilotesla quasistatic magnetic fields, the interaction of a laser with magnetized solids enters an unexplored new regime. In particular, a circularly polarized (CP) laser pulse may propagate in a highly magnetized plasma of any high density without encountering cutoff reflection in the whistler mode. With this, we propose a scheme for producing uniform warm dense matter (WDM) by direct laser heating with a CP laser irradiating onto the target along the magnetic field. It is shown by particle-in-cell simulations, which include advanced ionization dynamics and collision dynamics, moderately intense right-hand CP laser light at 10^{15}W/cm^{2} can propagate in solid aluminum and heat it efficiently to the 100 eV level within picoseconds. By using two laser pulses irradiating from two sides of a thin solid target, uniform heating to WDM can be achieved. This provides a controllable way to create WDM at different temperatures.

Transport of intense particle beams in large-scale plasmas

Transport of particle beams in plasmas is widely employed in fundamental research, industry, and medicine. Due to the high inertia of ion beams, their transport in plasmas is usually assumed to be stable. Here we report the focusing and flapping of intense slab proton beams transporting through large-scale plasmas by using a recently developed kinetic particle-in-cell simulation code. The beam self-focusing effect in the simulation is prominent and agrees well with previous experiments and theories. Moreover, the beam can curve and flap like turbulence as the beam density increases. Simulation and analysis indicate that the self-generated magnetic fields, produced by movement of collisional plasmas, are the dominant driver of such behaviors. By analyzing the spatial growth rate of magnetic energy and energy deposition of injected proton beams, it is found that the focusing and flapping are significantly determined by the injected beam densities and energies. In addition, a remarkable nonlinear beam energy loss is observed. Our research might find application in inertial confinement fusion and also might be of interest to the laboratory astrophysics community.

Particle-in-cell simulation of transport and energy deposition of intense proton beams in solid-state materials

A particle-in-cell (PIC) simulation code is used to investigate the transport and energy deposition of an intense proton beam in solid-state material. This code is able to simulate close particle interactions by using a Monte Carlo binary collision model. Such a model takes into account all related interactions between the incident protons and material particles, e.g., proton-nucleus, proton-bound-electron, and proton-free-electron collisions. This code also includes a Monte Carlo model for the collisional ionization and electron-ion recombination as well as the depression of the ionization potential by shielding of surrounding particles. Moreover, for intense proton beams, in order to include collective electromagnetic effects, significantly speed up the simulation, and simultaneously avoid numerical instabilities, an approach that combines the PIC method with a reduced model of high-density plasma based on Ohm's law is used. Simulation results indicate that the collective electromagnetic effects have a significant influence on the transport and energy deposition of proton beams. The Ohmic electric field would increase the stopping power and leads to a shortened range of proton beams in solid. The magnetic field would localize the energy deposition by collimating proton beams, which would otherwise be deflected by the collisions with nuclei.