Search

The fusion Z-pinch experiment (FuZE) is a sheared-flow stabilized Z-pinch designed to study the effects of flow stabilization on deuterium plasmas with densities and temperatures high enough to drive nuclear fusion reactions. Results from FuZE show high pinch currents and neutron emission durations thousands of times longer than instability growth times. While these results are consistent with thermonuclear neutron emission, energetically resolved neutron measurements are a stronger constraint on the origin of the fusion production. This stems from the strong anisotropy in energy created in beam-target fusion, compared to the relatively isotropic emission in thermonuclear fusion. In dense Z-pinch plasmas, a potential and undesirable cause of beam-target fusion reactions is the presence of fast-growing, “sausage” instabilities. This work introduces a new method for characterizing beam instabilities by recording individual neutron interactions in plastic scintillator detectors positioned at two different angles around the device chamber. Histograms of the pulse-integral spectra from the two locations are compared using detailed Monte Carlo simulations. These models infer the deuteron beam energy based on differences in the measured neutron spectra at the two angles, thereby discriminating beam-target from thermonuclear production. An analysis of neutron emission profiles from FuZE precludes the presence of deuteron beams with energies greater than 4.65 keV with a statistical uncertainty of 4.15 keV and a systematic uncertainty of 0.53 keV. This analysis demonstrates that axial, beam-target fusion reactions are not the dominant source of neutron emission from FuZE. These data are promising for scaling FuZE up to fusion reactor conditions.

Kinetic-ion, quasineutral, fluid-electron particle-in-cell simulations of interpenetrating carbon–carbon plasma flows in 2D RZ cylindrical geometry are presented. The simulations are initialized with solid density targets that are subsequently irradiated by 1014 W/cm2 intensity lasers using a raytracing package. The ablation, interpenetration, heating, slowing, entrainment, and stagnation of the plasma flows evolve self-consistently within the code. The particle density, velocity phase space, and fits to the velocity distribution functions are used, along with analytical collisional stopping rates, to interpret the dynamics of the flow evolution. Comparisons to multifluid simulations are described and used to highlight ion-kinetic effects in the setup. Synthetic Thomson scattering diagnostic signals are generated using detailed knowledge of the plasma distribution functions. The large scale of the system, 1 × 1 mm for 2 ns, and the detailed dynamics extracted demonstrate that such hybrid codes are powerful tools for the design and evaluation of laboratory-scale high-energy-density plasma physics experiments.Kinetic-ion, quasineutral, fluid-electron particle-in-cell simulations of interpenetrating carbon–carbon plasma flows in 2D RZ cylindrical geometry are presented. The simulations are initialized with solid density targets that are subsequently irradiated by 1014 W/cm2 intensity lasers using a raytracing package. The ablation, interpenetration, heating, slowing, entrainment, and stagnation of the plasma flows evolve self-consistently within the code. The particle density, velocity phase space, and fits to the velocity distribution functions are used, along with analytical collisional stopping rates, to interpret the dynamics of the flow evolution. Comparisons to multifluid simulations are described and used to highlight ion-kinetic effects in the setup. Synthetic Thomson scattering diagnostic signals are generated using detailed knowledge of the plasma distribution functions. The large scale of the system, 1 × 1 mm for 2 ns, and the detailed dynamics extracted demonstrate that such hybrid codes are powerfu...

The first fully kinetic particle-in-cell (PIC) simulations of sheared flow stabilized Z-pinch plasmas show the suppression of the sausage instability by shear, ∂rvz ≠ 0, with flow Mach numbers ≲1, consistent with experimental observations. Experimental investigations of sheared-flow stabilized Z-pinches demonstrated stability for 10 s of microseconds, over 1000 Alfven radial transit times, in quasi steady-state plasmas that are an intermediate between conventional inertial and magnetic confinement systems. The observed stability coincides with the presence of radial shear in axial flow profiles with peak speeds less than Mach 1, and experiments are underway to validate scaling this design to fusion conditions. The experimentally observed stability agrees with models of m = 1 kink mode suppression by sheared flows, but existing models of the m = 0 sausage mode underestimate the efficacy of sheared flow stabilization. These models rely on fluid approximations and find that stabilization requires flows ranging from Mach 1.7 to 4.3, and in some cases, stabilization is not reproduced in the models. This is faster than the measured flows in long-lived plasmas and would necessitate substantial energy convection out of the Z-pinch and the need to drive and sustain supersonic flows in future devices. The MHD models typically used in the literature are invalid in the high-temperature, high-current environments desirable for many Z-pinch applications, and they ignore large Larmor radius effects and viscous dissipation which are known to impact Z-pinch stability. PIC simulations can capture all these effects as well as kinetic instabilities that could influence the performance of high-temperature sheared flow stabilized Z-pinch plasmas. The PIC simulations presented here show the suppression and damping of m = 0 modes by sheared flows ∂rvz = 0.75vA/r0 with flow Mach numbers ≲1. Equivalent stability occurs under plasma conditions ranging from the limits of present-day experimental capabilities to the projected conditions of a sheared flow stabilized Z-pinch reactor.

Dense plasma focus (DPF) devices are conventionally operated with a polarity such that the inner electrode (IE) is the anode. It has been found that interchanging the polarity of the electrodes (i.e., IE as the cathode) can cause an order of magnitude decrease in the neutron yield. This polarity riddle has previously been studied empirically through several experiments and is yet not well understood. We have performed kinetic simulations using the particle-in-cell modeling to investigate the problem. This is the first time that both polarities have been studied with simulations in great detail. In our simulations, we have modeled the entire beam and plasma target formation processes, but we did not consider differences in break-down conditions caused by the two polarities. We have found that when using reverse polarity ions are still accelerated and, in fact, attain similar energy spectra as in the standard polarity case. The difference is that the fields are flipped and thus ions are accelerated in the opposite direction. So, in the reverse polarity case, the majority of the “plasma target” (formed by the imploding plasma) is in the opposite direction of the beam, and thus, the beam hits the IE and produces few neutrons. With a better inner electrode configuration, reverse polarity is able to create a high-quality ion beam as well as a high-density target. Both can be comparable to that generated by standard polarity. Furthermore, we will show that it is easier to add an additional solid catcher target to a DPF device with reverse polarity, potentially enabling it to generate more neutrons than standard polarity.

Collisional kinetic modifications of ion distributions in interpenetrating flows are investigated by irradiating two opposing targets, either CD/CD or CD/CH, on the National Ignition Facility. In the CD/CD case, neutron time-of-flight diagnostics are successfully used to infer the ion temperature, 5–6 keV, and velocity, 500 km/s per flow, of the flows using a multi-fluid approximation of beam-beam nuclear fusion. These values are found to be in agreement with simulations and other diagnostics. However, for CD/CH, the multi-fluid assumption breaks down, as fusion is quasi-thermonuclear in this case and thus more dependent on the details of the ion velocity distribution. Using kinetic-ion, hydrodynamic-electron, and hybrid particle-in-cell modeling, this is found to be partially due to a skewed deviation from a Maxwellian in the ion velocity distribution function resulting from ion-ion collisions. This skew causes a downshift in the mean neutron velocity that partially resolves the observation in the CD/CH case. We note that the discrepancy is not completely resolved via collisional effects alone and may be a signature of collisionless electromagnetic interactions such as the Weibel-filamentation instability.

Many astrophysical systems are effectively “collisionless,” that is, the mean free path for collisions between particles is much longer than the size of the system. The absence of particle collisions does not preclude shock formation, however, as shocks can be the result of plasma instabilities that generate and amplify electromagnetic fields. The magnetic fields required for shock formation may either be initially present, for example, in supernova remnants or young galaxies, or they may be self-generated in systems such as gamma-ray bursts (GRBs). In the case of GRB outflows, the Weibel instability is a candidate mechanism for the generation of sufficiently strong magnetic fields to produce shocks. In experiments on the OMEGA Laser, we have demonstrated a quasi-collisionless system that is optimized for the study of the non-linear phase of Weibel instability growth. Using a proton probe to directly image electromagnetic fields, we measure Weibel-generated magnetic fields that grow in opposing, initially unmagnetized plasma flows. The collisionality of the system is determined from coherent Thomson scattering measurements, and the data are compared to similar measurements of a fully collisionless system. The strong, persistent Weibel growth observed here serves as a diagnostic for exploring large-scale magnetic field amplification and the microphysics present in the collisional-collisionless transition.

In this paper, we detail the diagnostic technique used to infer the spatially resolved electron temperatures and densities in experiments dedicated to investigate the generation of magnetically collimated plasma jets. It is shown that the relative intensities of the resonance transitions in emitting He-like ions can be used to measure the temperature in such recombining plasmas. The intensities of these transitions are sensitive to the plasma density in the range of 1016–1020 cm−3 and to plasma temperature ranges from 10 to 100 eV for ions with a nuclear charge Zn ∼ 10. We show how detailed calculations of the emissivity of F VIII ions allow to determine the parameters of the plasma jets that were created using ELFIE ns laser facility (Ecole Polytechnique, France). The diagnostic and analysis technique detailed here can be applied in a broader context than the one of this study, i.e., to diagnose any recombining plasma containing He-like fluorine ions.

Transport of relativistic electrons in a solid Cu wire target has been modeled with the implicit hybrid particle-in-cell code LSP to investigate the electron energy distribution and energy coupling from the high-intensity, short-pulse laser to electrons entering to the wire. Experiments were performed on the TITAN laser using a 1.5 mm long Cu wire attached to a Au cone tip at the laser intensity of 1 × 1020 W/cm2 which was irradiated into the cone. The simulated Cu Kα wire profile and yields matched the measurements using a two-temperature energy distribution. These modeling results show that the cold component of the energy spectrum can be determined with ±100 keV accuracy from the fit to the initial experimental fall-off of the Kα emission while the simulated profiles were relatively insensitive to the hotter component of the electron distribution (>4 MeV). The slope of measured escaped electrons was used to determine the hotter temperature. Using exponential energy distributions, the laser-to-electron-in-wire coupling efficiencies inferred from the fits decreased from 3.4% to 1.5% as the prepulse energy increases up to 1 J. The comparison of the energy couplings using the exponential and Relativistic Maxwellian distribution functions showed that the energy inferred in the cold component is independent of the type of the distribution function.

A roadmap is suggested and demonstrated experimentally for the production of high-energy (>15 MeV) neutrons using short pulse lasers. Investigation with a 3D Monte Carlo model has been employed to quantify the production of energetic neutrons. Numerical simulations have been performed for three nuclear reactions, d(d,n)3He, 7Li(d,n)8Be, and 7Li(p,n)7Be, driven by monoenergetic ion beams. Quantitative estimates for the driver ion beam energy and number have been made and the neutron spectra and yield in the ion propagation direction have been evaluated for various incident ion energies. In order to generate neutron fluence above a detection limit of 106 neutrons/sr, either ∼1010 protons with energy 20–30 MeV or comparable amount of deuterons with energy 5–10 MeV are required. Experimental verification of the concept with deuterons driven by the Titan laser (peak intensity 2 × 1019 W/cm2, pulse duration of 9 ps, wavelength 1.05 μm, and energy of 360 J) is provided with the generation of neutrons with energy...

The generation of high-energy neutrons using laser-accelerated ions is demonstrated experimentally using the Titan laser with 360 J of laser energy in a 9 ps pulse. In this technique, a short-pulse, high-energy laser accelerates deuterons from a CD2 foil. These are incident on a LiF foil and subsequently create high energy neutrons through the 7Li(d,xn) nuclear reaction (Q = 15 MeV). Radiochromic film and a Thomson parabola ion-spectrometer were used to diagnose the laser accelerated deuterons and protons. Conversion efficiency into protons was 0.5%, an order of magnitude greater than into deuterons. Maximum neutron energy was shown to be angularly dependent with up to 18 MeV neutrons observed in the forward direction using neutron time-of-flight spectrometry. Absolutely calibrated CR-39 detected spectrally integrated neutron fluence of up to 8 × 108 n sr−1 in the forward direction.

Fast ignition requires a precise knowledge of fast electron propagation in a dense hydrogen plasma. In this context, a dedicated HiPER (High Power laser Energy Research) experiment was performed on the VULCAN laser facility where the propagation of relativistic electron beams through cylindrically compressed plastic targets was studied. In this paper, we characterize the plasma parameters such as temperature and density during the compression of cylindrical polyimide shells filled with CH foams at three different initial densities. X-ray and proton radiography were used to measure the cylinder radius at different stages of the compression. By comparing both diagnostics results with 2D hydrodynamic simulations, we could infer densities from 2 to 11 g/cm(3) and temperatures from 30 to 120 eV at maximum compression at the center of targets. According to the initial foam density, kinetic, coupled (sometimes degenerated) plasmas were obtained. The temporal and spatial evolution of the resulting areal densities and electrical conductivities allow for testing electron transport in a wide range of configurations. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3578346]

The relativistic electron transport induced by an ultraintense picosecond laser is experimentally investigated using an x-ray two-dimensional imaging system. Previous studies of the electron beam divergence [R. B. Stephens et al. Phys. Rev. E 69, 066414 (2004), for instance] were based on an x-ray imaging of a fluorescence layer buried at different depths in the target along the propagation axis. This technique required several shots to be able to deduce the divergence of the beam. Other experiments produced single-shot images in a one-dimensional geometry. The present paper describes a new target design producing a single-shot, two-dimensional image of the electrons propagating in the target. Several characteristics of the electron beam are extracted and discussed and Monte Carlo simulations provide a good understanding of the observed beam shape. The proposed design has proven to be efficient, reliable, and promising for further similar studies.

A neutron source for neutron resonance spectroscopy has been developed using high-intensity, short-pulse lasers. This technique will allow robust measurement of interior ion temperature of laser-shocked materials and provide insight into material equation of state. The neutron generation technique uses laser-accelerated protons to create neutrons in LiF through (p,n) reactions. The incident proton beam has been diagnosed using radiochromic film. This distribution is used as the input for a (p,n) neutron prediction code which is validated with experimentally measured neutron yields. The calculation infers a total fluence of 1.8×109 neutrons, which are expected to be sufficient for neutron resonance spectroscopy temperature measurements.

The Bremsstrahlung and K-shell emission from 1×1×1 mm3 planar targets irradiated by a short-pulse 3×1018–8×1019 W/cm2 laser were measured. The Bremsstrahlung was measured using a filter stack spectrometer with spectral discrimination up to 500 keV. K-shell emission was measured using a single photon counting charge coupled device. From Monte Carlo modeling of the target emission, conversion efficiencies into 1–3 MeV electrons of 3%–12%, representing 20%–40% total conversion efficiencies, were inferred for intensities up to 8×1019 W/cm2. Comparisons to scaling laws using synthetic energy spectra generated from the intensity distribution of the focal spot imply slope temperatures less than the ponderomotive potential of the laser. Resistive transport effects may result in potentials of a few hundred kV in the first few tens of microns in the target. This would lead to higher total conversion efficiencies than inferred from Monte Carlo modeling but lower conversion efficiencies into 1–3 MeV electrons.