Search

Kinetic modeling of laser-ion beam generation from the “break-out afterburner” (BOA) has been modeled for several deuteron-rich solid-density target foils. Modeling the transport of these beams in a beryllium converter shows as much as a fourfold increase in neutron yield over the present state of the art through the use of alternative target materials. Additionally, species-separation dynamics during the BOA can be exploited to control the hardness of the neutron spectra, of interest for, for example, enhancing penetrability in shielded material in active neutron interrogation settings.

Conventional explicit electromagnetic particle-in-cell (PIC) algorithms do not conserve discrete energy exactly. Time-centered fully implicit PIC algorithms can conserve discrete energy exactly, but may introduce large dispersion errors in the light-wave modes. This can lead to intolerable simulation errors where accurate light propagation is needed (e.g. in laser-plasma interactions). In this study, we selectively combine the leap-frog and Crank-Nicolson methods to produce an exactly energy- and charge-conserving relativistic electromagnetic PIC algorithm. Specifically, we employ the leap-frog method for Maxwell's equations, and the Crank-Nicolson method for the particle equations. The semi-implicit formulation still features a timestep CFL, but facilitates exact global energy conservation, exact local charge conservation, and preserves the dispersion properties of the leap-frog method for the light wave. The algorithm employs a new particle pusher designed to maximize efficiency and minimize wall-clock-time impact vs. the explicit alternative. It has been implemented in a code named iVPIC, based on the Los Alamos National Laboratory VPIC code ( https://github.com/losalamos/vpic ). We present numerical results that demonstrate the properties of the scheme with sample test problems: relativistic two-stream instability, Weibel instability, and laser-plasma instabilities.

The engineered macro-pore foam provides a new way to study thermonuclear burn physics by utilizing capsules containing deuterated (D) foam and filling tritium (T) gas in the engineered macro-pores. The implosion of a thermonuclear capsule filled with an engineered macro-pore foam will be complex due to the interaction of a shock wave with the engineered macro-pores. It is our goal to quantify how substantially complex foam structures affect the shape of shock and bulk shock speed. A cylinder-shape shock tube experiment has been designed and performed at the Omega Laser Facility. In order to examine how a foam structure will affect shock propagation, we performed several tests varying (1) engineered macro-pore size, (2) average foam density, and (3) with/without neopentane (C5H12) gas. X-ray radiographic data indicate that shock speed through engineered macro-pore foams depends strongly on average foam density and less on pore size. Experimental shock propagation data helped guide two numerical simulation approaches: (1) a 2D simulation with homogenizing foams rather than explicitly simulating engineered macro-pores and (2) a 2D toroidal-pore approximation adopting a toroidal-tube geometry to model engineered macro-pores.

The origin of the Finite Grid Instability (FGI) is studied by resolving the dynamics in the 1D electrostatic Particle-In-Cell (PIC) model in the spectral domain at the single particle level and at the collective motion level. The spectral fidelity of the PIC model is contrasted with the underlying physical system or the gridless model. The systematic spectral phase and amplitude errors from the charge deposition and field interpolation are quantified for common particle shapes used in the PIC models. It is shown through such analysis and in simulations that the lack of spectral fidelity relative to the physical system due to the existence of aliased spatial modes is the major cause of the FGI in the PIC model.

The MARBLE project is a novel inertial confinement fusion platform for studying the development of atomic mixing and temperature equilibration in inertial confinement fusion implosions and their impact on thermonuclear burn. Experiments involve the laser-driven implosion of capsules filled with deuterated engineered foams whose pores are filled with a gaseous mixture of hydrogen and tritium. By varying the size of the foam pores, we can study the timescale of the development of atomic mix relative to the development of thermal equilibrium between species. In contrast, previous separated reactant experiments have only provided information on the total amount of mix mass. We report on the series of MARBLE experiments [first reported in Haines et al., Nat. Commun. 11, 544 (2020)] performed on the University of Rochester's OMEGA laser facility and detailed and highly resolved three-dimensional radiation-hydrodynamic simulations of the implosions. In both the experimental and simulation results, we observe that the reactants do not achieve thermal equilibrium during the course of the implosion except in atomically mixed regions—i.e., that atomic mixing develops faster than thermal equilibration between species. The results suggest that ion temperature variations in the mixture are at least as important as reactant concentration variations for determining the fusion reaction rates.

Revolver and Double Shell Inertial Confinement Fusion capsule designs hope to achieve a robust volumetric thermonuclear burn via the use of a high-Z pusher shell filled with a cryogenic D–T fuel. Unfortunately, mix of the pusher material into the fuel (gas) may adversely impact the burn performance. Hydrodynamic instability of the metal/gas interface as the mix source is an obvious concern, but 1D effects may also be detrimental. Such effects include plasma diffusion at material interfaces, which has been the subject of numerous theoretical, computational, and experimental investigations. However, other 1D mix mechanisms may exist, which have yet to be thoroughly explored. In particular, plasma kinetic effects may drive the mix when a shock breaks out of the metal/gas interface. Using the state-of-the-art, hybrid (kinetic-ion/fluid electron), multi-ion Vlasov–Fokker–Planck code, iFP, we show that shock-driven kinetic effects can reconfigure the interface and the interfacial width subsequently grows diffusively. Finally, we consider any implications for high-Z pusher designs.

Experimental data are presented showing maximum carbon C6+ ion energies obtained from nm-scaled targets in the relativistic transparent regime for laser intensities between 9 × 1019 and 2 × 1021 W/cm2. When combined with two-dimensional particle-in-cell simulations, these results show a steep linear scaling for carbon ions with the normalized laser amplitude a0 ($a_0 \propto \sqrt ( I)$). The results are in good agreement with a semi-analytic model that allows one to calculate the optimum thickness and the maximum ion energies as functions of a0 and the laser pulse duration τλ for ion acceleration in the relativistic-induced transparency regime. Following our results, ion energies exceeding 100 MeV/amu may be accessible with currently available laser systems.

Laser-plasma interactions in the novel regime of relativistically induced transparency (RIT) have been harnessed to generate intense ion beams efficiently with average energies exceeding 10 MeV/nucleon (>100 MeV for protons) at “table-top” scales in experiments at the LANL Trident Laser. By further optimization of the laser and target, the RIT regime has been extended into a self-organized plasma mode. This mode yields an ion beam with much narrower energy spread while maintaining high ion energy and conversion efficiency. This mode involves self-generation of persistent high magnetic fields (∼104 T, according to particle-in-cell simulations of the experiments) at the rear-side of the plasma. These magnetic fields trap the laser-heated multi-MeV electrons, which generate a high localized electrostatic field (∼0.1 T V/m). After the laser exits the plasma, this electric field acts on a highly structured ion-beam distribution in phase space to reduce the energy spread, thus separating acceleration and energy-spread reduction. Thus, ion beams with narrow energy peaks at up to 18 MeV/nucleon are generated reproducibly with high efficiency (≈5%). The experimental demonstration has been done with 0.12 PW, high-contrast, 0.6 ps Gaussian 1.053 μm laser pulses irradiating planar foils up to 250 nm thick at 2–8 × 1020 W/cm2. These ion beams with co-propagating electrons have been used on Trident for uniform volumetric isochoric heating to generate and study warm-dense matter at high densities. These beam plasmas have been directed also at a thick Ta disk to generate a directed, intense point-like Bremsstrahlung source of photons peaked at ∼2 MeV and used it for point projection radiography of thick high density objects. In addition, prior work on the intense neutron beam driven by an intense deuterium beam generated in the RIT regime has been extended. Neutron spectral control by means of a flexible converter-disk design has been demonstrated, and the neutron beam has been used for point-projection imaging of thick objects. The plans and prospects for further improvements and applications are also discussed.

The nonlinear saturation of crossed-beam energy transfer (CBET) for multispeckled laser beams crossing at arbitrary angles is examined using vector particle-in-cell simulations. CBET is found to saturate on fast (∼10s of picosecond) time scales involving ion trapping and excitation of oblique forward stimulated Raman scattering (FSRS). Ion trapping reduces wave damping and speckle interaction increases wave coherence length, together enhancing energy transfer; ion acoustic wave (IAW) breakup in the direction transverse to the wavenumber increases wave damping and contributes to CBET saturation. The seed beam can become unstable to oblique FSRS, which leads to beam deflection at a large angle and a frequency downshift (by the plasma frequency). FSRS saturates on fast ∼picosecond time scales by electron plasma wave self-focusing, leading to enhanced side-loss hot electrons with energy exceeding 300 keV. This may contribute to fuel preheat but FSRS can be mitigated by the presence of a density gradient. Such growth of FSRS contributes to the saturation of CBET. Scaling simulations show that CBET, as well as FSRS and hot electrons, increases with beam average intensity, beam diameter, and crossing area, but that CBET is limited by the excitation of FSRS and IAW breakups in addition to pump depletion. FSRS deflects the seed beam energy by greater than 40% of the incident beam energy and puts a few percent of the incident beam energy into hot electrons. FSRS limits the efficacy of CBET for symmetry tuning at late stages in the implosion and may account for a large portion of the “missing energy” in implosions that use gas-filled hohlraums.

The physics of mixing in plasmas is of fundamental importance to inertial confinement fusion and high energy density laboratory experiments. Two- and three-dimensional (2D and 3D) particle-in-cell simulations with a binary collision model are used to explore kinetic effects arising during the mixing of plasma media. The applicability of the one-dimensional (1D) ambipolarity condition is evaluated in 2D and 3D simulations of a plasma interface with a sinusoidal perturbation. The 1D ambipolarity condition is found to remain valid in 2D and 3D, as electrons and ions flow together required for J = 0. Simulations of perturbed interfaces show that diffusion-induced total pressure imbalance and hydroflows flatten fine interface structures and drive rapid atomic mix. The atomic mix rate from a structured interface is faster than the ∼ t scaling obtained from 1D theory in the small-Knudsen-number limit. Plasma kinetic effects inhibit the growth of the Rayleigh-Taylor instability at small wavelengths and result in a nonmonotonic growth rate scaling with wavenumber k with a maximum at a low k value, much different from Agk (where A is the Atwood number and g is the gravitational constant) as expected in the absence of plasma kinetic effects. Simulations under plasma conditions relevant to MARBLE separated-reactant experiments on Omega and the NIF show kinetic modification of DT fusion reaction rates. With non-Maxwellian distributions and relative drifts between D and T ions, DT reactivity is higher than that inferred from rates using stationary Maxwellian distributions. Reactivity is also found to be reduced in the presence of finite-Knudsen-layer losses.

Recent progress in laser-driven quasi-monoenergetic ion beams enabled the production of uniformly heated warm dense matter. Matter heated rapidly with this technique is under extreme temperatures and pressures, and promptly expands outward. While the expansion speed of an ideal plasma is known to have a square-root dependence on temperature, computer simulations presented here show a linear dependence of expansion speed on initial plasma temperature in the warm dense matter regime. The expansion of uniformly heated 1–100 eV solid density gold foils was modeled with the RAGE radiation-hydrodynamics code, and the average surface expansion speed was found to increase linearly with temperature. The origin of this linear dependence is explained by comparing predictions from the SESAME equation-of-state tables with those from the ideal gas equation-of-state. These simulations offer useful insight into the expansion of warm dense matter and motivate the application of optical shadowgraphy for temperature measurement.

A new class of inertial fusion capsules is presented that combines multishell targets with laser direct drive at low intensity ($2.8\ifmmode\times\else\texttimes\fi{}{10}^{14}\text{ }\text{ }\mathrm{W}/{\mathrm{cm}}^{2}$) to achieve robust ignition. The targets consist of three concentric, heavy, metal shells, enclosing a volume of tens of $\ensuremath{\mu}\mathrm{g}$ of liquid deuterium-tritium fuel. Ignition is designed to occur well ``upstream'' from stagnation, with minimal pusher deceleration to mitigate interface Rayleigh-Taylor growth. Laser intensities below thresholds for laser plasma instability and cross beam energy transfer facilitate high hydrodynamic efficiency ($\ensuremath{\sim}10%$).

When electrons are accelerated to near light-speeds through an overdense plasma by an intense laser beam, the usually opaque plasma becomes optically transparent. High-speed laser experiments provide unprecedented insight into the dynamics of this process.

Experiments demonstrate the ~77× amplification of 0.5 to 3.5-ps pulses of seed light by interaction with Langmuir waves in a low density (1.2 × 1019 cm−3) plasma produced by a 1-ns, 230-J, 1054-nm pump beam with 1.2 × 1014 W/cm2 intensity. The waves are strongly damped (kλD = 0.38, Te = 244 eV) and grow over a ~ 1 mm length, similar to what is experienced by scattered light when it interacts with crossing beams as it exits an ignition target. The amplification reduces when the seed intensity increases above ~1 × 1011 W/cm2, indicating that saturation of the plasma waves on the electron kinetic time scale (

Double shell capsules are predicted to ignite and burn at relatively low temperature (∼3 keV) via volume ignition and are a potential low-convergence path to substantial α-heating and possibly ignition at the National Ignition Facility. Double shells consist of a dense, high-Z pusher, which first shock heats and then performs work due to changes in pressure and volume (PdV work) on deuterium-tritium gas, bringing the entire fuel volume to high pressure thermonuclear conditions near implosion stagnation. The high-Z pusher is accelerated via a shock and subsequent compression of an intervening foam cushion by an ablatively driven low-Z outer shell. A broad capsule design parameter space exists due to the inherent flexibility of potential materials for the outer and inner shells and foam cushion. This is narrowed down by design physics choices and the ability to fabricate and assemble the separate pieces forming a double shell capsule. We describe the key physics for good double shell performance, the trade-offs in various design choices, and the challenges for capsule fabrication. Both 1D and 2D calculations from radiation-hydrodynamic simulations are presented.

The classical transport theory based on Chapman-Enskog methods provides self-consistent approximations for the kinetic flux of mass, heat, and momentum in a fluid limit characterized with a small Knudsen number. The species mass fluxes relative to the center of mass, or “diffusive fluxes,” are expressed as functions of known gradient quantities with kinetic coefficients evaluated using similar analyses for mixtures of gases or plasma components. The sum over species of the diffusive mass fluxes is constrained to be zero in the Lagrange frame, and thus results in a non-zero molar flux leading to a pressure perturbation. At an interface between two species initially in pressure equilibrium, the pressure perturbation driven by the diffusive molar flux induces a center of mass velocity directed from the species of greater atomic mass towards the lighter atomic mass species. As the ratio of the species particle masses increases, this center of mass velocity carries an increasingly greater portion of the mass across the interface and for a particle mass ratio greater than about two, the center of mass velocity carries more mass than the gradient driven diffusion flux. Early time transients across an interface between two species in a 1D plasma regime and initially in equilibrium are compared using three methods; a fluid code with closure in a classical transport approximation, a particle in cell simulation, and an implicit Fokker-Planck solver for the particle distribution functions. The early time transient phenomenology is shown to be similar in each of the computational simulation methods, including a pressure perturbation associated with the stationary “induced” component of the center of mass velocity which decays to pressure equilibrium during diffusion. At early times, the diffusive process generates pressure and velocity waves which propagate outward from the interface and are required to maintain momentum conservation. The energy in the outgoing waves dissipates as heat in viscous regions, and it is hypothesized that these diffusion driven waves may sustain fluctuations in less viscid finite domains after reflections from the boundaries. These fluid dynamic phenomena are similar in gases or plasmas and occur in flow transients with a moderate Knudsen number. The analysis and simulation results show how the kinetic flux, represented in the fluid transport closure, directly modifies the mass averaged flow described with the Euler equations.

We examine the relativistic Buneman instability in systems relevant to high-intensity laser-plasma interactions under conditions of relativistically-induced transparency, as this instability can generate large-amplitude electrostatic waves at low frequencies that are pertinent to ion dynamics in these systems. Ion flows are shown to significantly alter the range of unstable wave numbers and to increase the phase velocities of the unstable modes; we particularly highlight the relativistic effects from both the ion and electron (with transverse motion) populations. These findings are related to the mode structure seen in particle-in-cell simulation results of a short-pulse laser breaking through an initially opaque target with the onset of relativistic transparency. Additionally, driving mechanisms from free energy present in density and velocity gradients are shown to be capable of significantly enhancing the growth rates, and these instabilities furthermore extend the breadth of the unstable wave number range. Lastly, we discuss how the transverse self-generated magnetic fields characteristic of short-pulse interactions can potentially constrain the unstable wave numbers in a non-trivial manner.

We present a particle-in-cell study of linearly polarized laser-ion acceleration systems, in which we use both two-dimensional (2D) and three-dimensional (3D) simulations to characterize the ion acceleration mechanisms in targets which become transparent to the laser pulse during irradiation. First, we perform a target length scan to optimize the peak ion energies in both 2D and 3D, and the predictive capabilities of 2D simulations are discussed. Tracer analysis allows us to isolate the acceleration into stages of target normal sheath acceleration (TNSA), hole boring (HB), and break-out afterburner (BOA) acceleration, which vary in effectiveness based on the simulation parameters. The thinnest targets reveal that enhanced TNSA is responsible for accelerating the most energetic ions, whereas the thickest targets have ions undergoing successive phases of HB and TNSA (in 2D) or BOA and TNSA (in 3D); HB is not observed to be a dominant acceleration mechanism in the 3D simulations. It is in the intermediate optimal regime, both when the laser breaks through the target with appreciable amplitude and when there is enough plasma to form a sustained high density flow, that BOA is most effective and is responsible for the most energetic ions. Eliminating the transverse laser spot size effects by performing a plane wave simulation, we can isolate with greater confidence the underlying physics behind the ion dynamics we observe. Specifically, supplemented by wavelet and FFT analyses, we match the post-transparency BOA acceleration with a wave-particle resonance with a high-amplitude low-frequency electrostatic wave of increasing phase velocity, consistent with that predicted by the Buneman instability.We present a particle-in-cell study of linearly polarized laser-ion acceleration systems, in which we use both two-dimensional (2D) and three-dimensional (3D) simulations to characterize the ion acceleration mechanisms in targets which become transparent to the laser pulse during irradiation. First, we perform a target length scan to optimize the peak ion energies in both 2D and 3D, and the predictive capabilities of 2D simulations are discussed. Tracer analysis allows us to isolate the acceleration into stages of target normal sheath acceleration (TNSA), hole boring (HB), and break-out afterburner (BOA) acceleration, which vary in effectiveness based on the simulation parameters. The thinnest targets reveal that enhanced TNSA is responsible for accelerating the most energetic ions, whereas the thickest targets have ions undergoing successive phases of HB and TNSA (in 2D) or BOA and TNSA (in 3D); HB is not observed to be a dominant acceleration mechanism in the 3D simulations. It is in the intermediate op...

We construct and investigate a set of stochastic differential equations that incorporate the physics of velocity-dependent small-angle Coulomb collisions among the plasma particles in a particle-in-cell simulation. Each particle is scattered stochastically from all the other particles in a simulation cell modeled as one or more Maxwellians. Total energy and momentum are conserved by linear transformation of the velocity increments. In two test simulations the proposed ''particle-moment'' collision algorithm performs well with time steps as large as 10% of the relaxation time - far larger than a particle-pairing collision algorithm, in which pairs of particles are scattered from one another, requires to achieve the same accuracy.

Short pulse laser plasma interaction experiments using diffraction limited beams provide an excellent platform to investigate the fundamental physics of stimulated Raman scattering. Detailed understanding of these laser plasma instabilities impacts the current inertial confinement fusion ignition designs and could potentially impact fast ignition when higher energy lasers are used with longer pulse durations (>1 kJ and >1 ps). Using short laser pulses, experiments can be modeled over the entire interaction time of the laser using particle-in-cell codes to validate our understanding quantitatively. Experiments have been conducted at the Trident laser facility and the Laboratoire pour l'Utilisation des Lasers Intenses (LULI) to investigate stimulated Raman scattering near the threshold of the instability using 527 nm and 1059 nm laser light, respectively, with 1.5–3.0 ps pulses. In both experiments, the interaction beam was focused into pre-ionized helium gas-jet plasma. Measurements of the reflectivity as a function of intensity and kλD were completed at the Trident laser facility, where k is the electron plasma wave number and λD is the plasma Debye length. At LULI, a 300 fs Thomson scattering probe is used to directly measure the density fluctuations of the driven electron plasma and ion acoustic waves. Work is currently underway comparing the results of the experiments with simulations using the VPIC particle-in-cell code. Details of the experimental results are presented in this manuscript.

Los Alamos National Laboratory short pulse experiments have shown using various target cleaning techniques such that heavy ion beams of different charge states can be produced. Furthermore, by controlling the thickness of light ions on the rear of the target, monoenergetic ion pulses can be generated. The spectral shape of the accelerated particles can be controlled to yield a range of distributions, from Maxwellian to ones possessing a monoenergetic peak at high energy. The key lies in understanding and utilizing target surface chemistry. Careful monitoring and control of the surface properties and induction of reactions at different temperatures allows well defined source layers to be formed, which in turn lead to the desired energy spectra in the acceleration process. Theoretical considerations provide understanding of the process of monoenergetic ion production. In addition, numerical modeling has identified a new acceleration mechanism, the laser break-out afterburner that could potentially boost particle energies by up to two orders of magnitude for the same laser parameters. This mechanism may enable application of laser-accelerated ion beams to venues such as compact accelerators, tumor therapy, and ion fast ignition.

With the development of several novel heating sources, scientists can now heat a small sample isochorically above 10,000 K. Although matter at such an extreme state, known as warm dense matter, is commonly found in astrophysics (e.g., in planetary cores) as well as in high energy density physics experiments, its properties are not well understood and are difficult to predict theoretically. This is because the approximations made to describe condensed matter or high-temperature plasmas are invalid in this intermediate regime. A sufficiently large warm dense matter sample that is uniformly heated would be ideal for these studies, but has been unavailable to date. Here we have used a beam of quasi-monoenergetic aluminum ions to heat gold and diamond foils uniformly and isochorically. For the first time, we visualized directly the expanding warm dense gold and diamond with an optical streak camera. Furthermore, we present a new technique to determine the initial temperature of these heated samples from the measured expansion speeds of gold and diamond into vacuum. We anticipate the uniformly heated solid density target will allow for direct quantitative measurements of equation-of-state, conductivity, opacity and stopping power of warm dense matter, benefiting plasma physics, astrophysics and nuclear physics.

In a recent experiment on the Trident laser facility, a laser-driven beam of quasi-monoenergetic aluminum ions was used to heat solid gold and diamond foils isochorically to 5.5 eV and 1.7 eV, respectively. Here theoretical calculations are presented that suggest the gold and diamond were heated uniformly by these laser-driven ion beams. According to calculations and SESAME equation-of-state tables, laser-driven aluminum ion beams achievable on Trident, with a finite energy spread of (delta E)/E ~ 20%, are expected to heat the targets more uniformly than a beam of 140 MeV aluminum ions with zero energy spread. The robustness of the expected heating uniformity relative to the changes in the incident ion energy spectra is evaluated, and expected plasma temperatures of various target materials achievable with the current experimental platform are presented., Comment: 18 pages, 8 figures

The authors use explicit particle-in-cell simulations to model stimulated scattering processes in media with both solitary and multiple laser hot spots. These simulations indicate coupling among hot spots, whereby scattered light, plasma waves, and hot electrons generated in one laser hot spot may propagate to neighboring hot spots, which can be destabilized to enhanced backscatter. A nonlinear statistical model of a stochastic beam exhibiting this coupled behavior is described here. Calibration of the model using particle-in-cell simulations is performed, and a threshold is derived for "detonation" of the beam to high reflectivity.

A new laser-driven ion acceleration mechanism has been identified using particle-in-cell (PIC) simulations. This mechanism allows ion acceleration to GeV energies at vastly reduced laser intensities compared with earlier acceleration schemes. The new mechanism, dubbed “Laser Break-out Afterburner” (BOA), enables the acceleration of carbon ions to greater than 2 GeV energy at a laser intensity of only 1021W/cm2, an intensity that has been realized in existing laser systems. Other techniques for achieving these energies in the literature rely upon intensities of 1024W/cm2or above, i.e., 2–3 orders of magnitude higher than any laser intensity that has been demonstrated to date. Also, the BOA mechanism attains higher energy and efficiency than target normal sheath acceleration (TNSA), where the scaling laws predict carbon energies of 50 MeV/u for identical laser conditions. In the early stages of the BOA, the carbon ions accelerate as a quasi-monoenergetic bunch with median energy higher than that realized recently experimentally.

The interaction of a high intensity, short-pulse laser with a thin target can lead to the generation of a highly collimated beam of fast ions off the rear target surface. These ion beams have the potential to impact inertial confinement fusion applications, including their use in diagnostics and fast-ignition. Recent work by the authors in the modeling of ion acceleration, using both full particle-in-cell and hybrid (particle ions, reduced electron physics) models, is leading to improved understanding of the physics governing ion acceleration.

Quiet direct simulation Monte Carlo (QDSMC) is a new particle simulation technique that is applicable to a broad range of applications where the underlying system dynamics obey Fokker Planck equations. These include hydrodynamics, radiation transport, magnetohydrodynamics, diffusion, and collisional kinetic plasmas. At the beginning of each time step in QDSMC, the weights and abscissas of Gaussian-Hermite quadrature are used to deterministically create particles to sample the random process. At the end of the time step, particles are gathered to the computational mesh to obtain updated distributions of conserved quantities on the mesh and then the particles are destroyed. The creation and destruction of particles allows arbitrary dynamical range to be accessed quiescently with only a small number of particles per computational cell. The application of QDSMC to the simulation of Coulomb collisions is considered in this report, and the method is demonstrated on problems involving the collisional relaxation of non-Maxwellian distributions.

One of the outstanding challenges of solar physics is understanding the mechanism for heating the solar corona, a low-β, low-resistivity plasma driven slowly by the shuffling of “footpoints” (ends of flux tubes) anchored in the photosphere. The corona’s low resistivity means that enormous current densities must be present for Ohmic dissipation to be effective. Two ways that this can occur in a slowly evolving, “quasiequilibrium” medium are by the formation of current sheets or by the excitation of very low-frequency Alfven resonances. Both processes lead to singular current densities in the ideal limit, and both may result from the excitation of very low- (or zero-) frequency shear Alfven resonances. Thus, understanding the conditions under which zero-frequency Alfven resonances can occur may ultimately be important in understanding coronal heating. This paper examines shear Alfven resonances in line-tied magnetic geometries, with special attention being paid to geometries with discontinuities in the fiel...

We model radiation transport by advancing computational photons through phase space with solutions to a set of stochastic differential equations of motion. Random numbers that appear in the equations of motion are sampled with deterministically chosen Gaussian quadrature weights and abscissas. In this way, the advantages of particle Monte Carlo are realized without generating statistical noise. We demonstrate this technique by performing one- and two-dimensional test problems in which gray radiation is energetically coupled to stationary material. Scattering is accomplished with a Fokker–Planck scattering operator. Free streaming, diffusion and Marshak waves are recovered in appropriate limits.

A new approach to particle simulation, called “quiet direct simulation Monte Carlo” (QDSMC), is described that can be applied to many problems of interest, including hydrodynamics, magnetohydrodynamics (MHD), and the modeling of collision plasmas. The essence of QDSMC is the use of carefully chosen weights for the particles (e.g., Gauss–Hermite, for Maxwellian distributions), which are destroyed each time step after the particle information is deposited onto the grid and reconstructed at the beginning of the next time step. The method overcomes the limited dynamical range and statistical noise typically found in particle simulations. In this article QDSMC is applied to hydrodynamics and MHD test problems, and its suitability for modeling semi-collisional plasma dynamics is considered.

A particle-in-cell study of laser-ion acceleration mechanisms in the transparency regime illustrates how two-dimensional (2D) S and P simulations (laser polarization in and out of the simulation plane, respectively) capture different physics characterizing these systems, visible in their entirety often in cost-prohibitive three-dimensional (3D) simulations. The electron momentum anisotropy induced in the target by a laser pulse is dramatically different in the two 2D cases, manifested in differences in target expansion timescales, electric field strengths, and density thresholds for the onset of relativistically induced transparency. In particular, 2D-P simulations exhibit dramatically greater electron heating in the simulation plane, whereas 2D-S ones show a much more isotropic energy distribution, similar to 3D. An ion trajectory analysis allows one to isolate the fields responsible for ion acceleration and to characterize the acceleration regimes in time and space. The artificial longitudinal electron ...

The loss of fuel ions in the Gamow peak and other kinetic effects related to the α particles during ignition, run-away burn, and disassembly stages of an inertial confinement fusion D-T capsule are investigated with a quasi-1D hybrid volume ignition model that includes kinetic ions, fluid electrons, Planckian radiation photons, and a metallic pusher. The fuel ion loss due to the Knudsen effect at the fuel-pusher interface is accounted for by a local-loss model by Molvig et al. [Phys. Rev. Lett. 109, 095001 (2012)] with an albedo model for ions returning from the pusher wall. The tail refilling and relaxation of the fuel ion distribution are captured with a nonlinear Fokker-Planck solver. Alpha heating of the fuel ions is modeled kinetically while simple models for finite alpha range and electron heating are used. This dynamical model is benchmarked with a 3 T hydrodynamic burn model employing similar assumptions. For an energetic pusher (∼40 kJ) that compresses the fuel to an areal density of ∼1.07g/cm2 at ignition, the simulation shows that the Knudsen effect can substantially limit ion temperature rise in runaway burn. While the final yield decreases modestly from kinetic effects of the α particles, large reduction of the fuel reactivity during ignition and runaway burn may require a higher Knudsen loss rate compared to the rise time of the temperatures above ∼25 keV when the broad D-T Gamow peak merges into the bulk Maxwellian distribution.

We measured the average deuterium cluster size within a mixture of deuterium clusters and helium gas by detecting Rayleigh scattering signals. The average cluster size from the gas mixture was comparable to that from a pure deuterium gas when the total backing pressure and temperature of the gas mixture were the same as those of the pure deuterium gas. According to these measurements, the average size of deuterium clusters depends on the total pressure and not the partial pressure of deuterium in the gas mixture. To characterize the cluster source size further, a Faraday cup was used to measure the average kinetic energy of the ions resulting from Coulomb explosion of deuterium clusters upon irradiation by an intense ultrashort pulse. The deuterium ions indeed acquired a similar amount of energy from the mixture target, corroborating our measurements of the average cluster size. As the addition of helium atoms did not reduce the resulting ion kinetic energies, the reported results confirm the utility of using a known cluster source for beam-target-fusion experiments by introducing a secondary target gas.

A fluid model is presented that describes the linear evolution of a thin layer of beam ions in a conducting cylindrical cavity with a narrow rotating beam layer, and a dispersion relation is computed for the low-frequency electromagnetic response of the system. For a broad range of ion ring parameters this configuration is found to be unstable to the generation of Alfvenic fluctuations. Simple estimates are provided of the time needed for these fluctuations to self-trap the ion beam.

Stochastic magnetic fields in a weakly collisional plasma enhance the collisionality of electrons in the medium compared to the classical Spitzer values. The reduction in the thermal conductivity may play an important role in the evolution of galaxy clusters. This work extends prior work by examining the parallel propagation of electrons in plasma media where the magnetic field amplitude sampled along a field line is a stochastic function of the distance. New physics is exhibited when many large-amplitude magnetic mirrors are present; the parallel transport becomes subdiffusive rather than diffusive in this limit. The transition between diffusion and subdiffusion can be obtained from properties of the asymptotic solution to the kinetic equation in the vicinity of a solitary magnetic mirror. This transition has been computed as a function of the power law exponent in the distribution of mirror maxima, and the asymptotic scaling of the mean-squared electron displacement with time has been obtained in these systems as a function of this exponent. Monte Carlo comparisons are in agreement with these results.

The interaction of intense short pulse radiation with thick Al foils is studied as a function of the absorption region density profile and the laser spot shape. Absorption of the light generates relativistic hot electrons near the critical surface. When the density profile is steep with micron scale lengths, and a small spot (8 μm FWHM) the hot electrons are retained near the surface while undergoing strong lateral surface transport through intense thermoelectric magnetic fields. Alternatively, with mild, ∼30 μm, scale length initial profiles the light beam bores a hole in the corona, wraps it with B-field, extinguishes the lateral hot electron flow, and sends a reduced fraction of hot electrons forward in filaments. Finally, broadening the spot to 40 μm and totally flattening it gives strong forward-directed hot electron penetration, which, if achievable, could serve effectively for Fast Ignition and radiography.

Here we present experimental results on laser-driven ion accel- eration from relativistically transparent, overdense plasmas in the break-out afterburner (BOA) regime. Experiments were preformed at the Trident ultra-high contrast laser facility at Los Alamos National Laboratory, and at the Texas Petawatt laser facility, located in the University of Texas at Austin. It is shown that when the target becomes relativistically transparent to the laser, an epoch of dramatic acceleration of ions occurs that lasts until the electron density in the expanding target reduces to the critical density in the non-relativistic limit. For given laser parameters, the optimal target thickness yielding the highest maximum ion energy is one in which this time window for ion acceleration overlaps with the intensity peak of the laser pulse. A simple analytic model of relativistically induced transparency is presented for plasma expansion at the

Summary form only given. The success of Inertial Confinement Fusion (ICF) is to achieve controlled thermonuclear burn in the laboratory which will lead to the commercialization of clean, carbon-free and safe Inertial Fusion Energy (IFE). Both ICF and IFE demand a detailed understanding of the rapidly evolving high energy density plasmas (HEDP) as intense lasers create and nonlinearly modify them. We have developed and tested new design tools for novel ultrafast diagnostics that use nonlinear optical (NLO) techniques to ferret out the complex, nonlinear, kinetic, microscopic dynamics of HEDP. Measuring the slope of the velocity distribution function of a plasma electron or ion species in a velocity sector of interest is one such paramount goal. We accomplish this by (i) adopting the appropriate method of generating a pump laser composed of spike trains of uneven duration and delay (STUD pulses)1, 2, (ii) adopting the appropriate method of detecting and diagnosing the amplified transmission of a stimulated Raman or stimulated Brillouin scattered (SRS or SBS) probe beam, and (iii) utilizing the gain variations of the scattered signal to develop a detailed map of background plasma instabilities. This GeDeCo code is being tested using output from state of the art kinetic simulations3 to emulate the microscopic state of an HED plasma. High-repetition-rate, high-average-power future drivers of IFE will use STUD pulses in order to control undesirable instabilities adaptively.

Experimental data are presented for laser-driven carbon C6+ ion-acceleration, verifying 2D-PIC studies for multi-species targets in the Break-Out Afterburner regime. With Trident's ultra-high contrast at relativistic intensities of 5 × 1020 W/cm2 and nm-scale diamond targets, acceleration of carbon ions has been optimized by using target laser-preheating for removal of surface proton contaminants. Using a high-resolution wide angle spectrometer, carbon C6+ ion energies exceeding 1 GeV or 83 MeV/amu have been measured, which is a 40% increase in maximum ion energy over uncleaned targets. These results are consistent with kinetic plasma modeling and analytic theory.

We present experimental data and PIC simulations on laser driven ion acceleration in the relativistic transparent regime. We measured carbon C6+ ions energies exceeding 1GeV and proton energies exceeding 100MeV. Conversion efficiency of laser light into ions, beam shape and scaling laws are presented for this regime and compared to the Target Normal Sheath acceleration and Radiation pressure acceleration regime.

Mixing at interfaces in dense plasma media is a problem central to inertial confinement fusion and high energy density laboratory experiments. In this work, collisional particle-in-cell simulations are used to explore kinetic effects arising during the mixing of unmagnetized plasma media. Comparisons are made to the results of recent analytical theory in the small Knudsen number limit and while the bulk mixing properties of interfaces are in general agreement, some differences arise. In particular, “super-diffusive” behavior, large diffusion velocity, and large Knudsen number are observed in the low density regions of the species mixing fronts during the early evolution of a sharp interface prior to the transition to a slow diffusive process in the small-Knudsen-number limit predicted by analytical theory. A center-of-mass velocity profile develops as a result of the diffusion process and conservation of momentum.

Here, kinetic effects on inertial confinement fusion have been investigated. In particular, inter-ion-species diffusion and suprathermal ion distribution have been analyzed. The former drives separation of the fuel constituents in the hot reacting core and governs mix at the shell/fuel interface. The latter underlie measurements obtained with nuclear diagnostics, including the fusion yield and inferred ion burn temperatures. Basic mechanisms behind and practical consequences from these effects are discussed.

Two- and three-dimensional particle-in-cell simulations of stimulated Brillouin scattering(SBS) in laser speckle geometry have been analyzed to evaluate the relative importance of competing nonlinear processes in the evolution and saturation of SBS. It is found that ion-trapping-induced wavefront bowing and breakup of ion acoustic waves(IAW) and the associated side-loss of trapped ions dominate electron-trapping-induced IAW wavefront bowing and breakup, as well as the two-ion-wave decay instability over a range of ZTe/Ti conditions and incident laser intensities. In the simulations, the latter instability does not govern the nonlinear saturation of SBS; however, evidence of two-ion-wave decay is seen, appearing as a modulation of the ion acoustic wavefronts. This modulation is periodic in the laser polarization plane, anti-symmetric across the speckle axis, and of a wavenumber matching that of the incident laser pulse. A simple analytic model is provided for how spatial “imprinting” from a high frequency inhomogeneity (in this case, the density modulation from the laser) in an unstable system with continuum eigenmodes can selectively amplify modes with wavenumbers that match that of the inhomogeneity.

A simple energy balance model has been obtained for laser-plasma heating in indirect drive hohlraum plasma that allows rapid temperature scaling and evolution with parameters such as plasma density and composition. Furthermore, this model enables assessment of the effects on plasma temperature of, e.g., adding high-Z dopant to the gas fill or magnetic fields.