Kyle Robert Cochrane

Direct observation of an abrupt insulator-to-metal transition in dense liquid deuterium

Driving liquid deuterium into metal Quick and powerful compression can force materials to change their properties dramatically. Knudson et al. compressed liquid deuterium to extreme temperatures and pressures using high-energy magnetic pulses at the Sandia Z-machine (see the Perspective by Ackland). Deuterium began to reflect like a mirror during compression, as the electrical conductivity sharply increased. The observed conditions for metallization of deuterium and hydrogen help us to build theoretical models for the universes most abundant element. This a our understanding of the internal layering of gas giant planets such as Jupiter and Saturn. Science, this issue p. 1455; see also p. 1429 Magnetic compression drives an insulator-to-metal transition in dense liquid deuterium. [Also see Perspective by Ackland] Eighty years ago, it was proposed that solid hydrogen would become metallic at sufficiently high density. Despite numerous investigations, this transition has not yet been experimentally observed. More recently, there has been much interest in the analog of this predicted metallic transition in the dense liquid, due to its relevance to planetary science. Here, we show direct observation of an abrupt insulator-to-metal transition in dense liquid deuterium. Experimental determination of the location of this transition provides a much-needed benchmark for theory and may constrain the region of hydrogen-helium immiscibility and the boundary-layer pressure in standard models of the internal structure of gas-giant planets.Eighty years ago, it was proposed that solid hydrogen would become metallic at sufficiently high density. Despite numerous investigations, this transition has not yet been experimentally observed. More recently, there has been much interest in the analog of this predicted metallic transition in the dense liquid, due to its relevance to planetary science. Here, we show direct observation of an abrupt insulator-to-metal transition in dense liquid deuterium. Experimental determination of the location of this transition provides a much-needed benchmark for theory and may constrain the region of hydrogen-helium immiscibility and the boundary-layer pressure in standard models of the internal structure of gas-giant planets.

Magnetically accelerated, ultrahigh velocity flyer plates for shock wave experiments

The intense magnetic field produced by the 20 MA Z accelerator is used as an impulsive pressure source to accelerate metal flyer plates to high velocity for the purpose of performing plate impact, shock wave experiments. This capability has been significantly enhanced by the recently developed pulse shaping capability of Z, which enables tailoring the rise time to peak current for a specific material and drive pressure to avoid shock formation within the flyer plate during acceleration. Consequently, full advantage can be taken of the available current to achieve the maximum possible magnetic drive pressure. In this way, peak magnetic drive pressures up to 490 GPa have been produced, which shocklessly accelerated 850μm aluminum (6061-T6) flyer plates to peak velocities of 34km∕s. We discuss magnetohydrodynamic (MHD) simulations that are used to optimize the magnetic pressure for a given flyer load and to determine the shape of the current rise time that precludes shock formation within the flyer during ac...

Measurements and simulations of the ablation stage of wire arrays with different initial wire sizes

Comparisons of 20mm diameter, 300-wire tungsten arrays with different initial wire sizes were made on the 20MA Sandia Z facility. Radiographic measurements of each wire array, taken at the same point in the current during the wire ablation stage, show systematic differences. A detailed comparison of the radiography and self-emission data with simulations and analytic models suggests that a variation in the mass ablation rate with wire size may be responsible.

Solid liner implosions on Z for producing multi-megabar, shockless compressions

Current pulse shaping techniques, originally developed for planar dynamic material experiments on the Z-machine [M. K. Matzen et al., Phys. Plasmas 12, 055503 (2005)], are adapted to the design of controlled cylindrical liner implosions. By driving these targets with a current pulse shape that prevents shock formation inside the liner, shock heating is avoided along with the corresponding decrease in electrical conductivity ahead of the magnetic diffusion wave penetrating the liner. This results in an imploding liner with a significant amount of its mass in the solid phase and at multi-megabar pressures. Pressures in the solid region of a shaped pulse driven beryllium liner fielded on the Z-machine are inferred to 5.5 Mbar, while simulations suggest implosion velocities greater than 50kms-1. These solid liner experiments are diagnosed with multi-frame monochromatic x-ray backlighting which is used to infer the material density and pressure. This work has led to a new platform on the Z-machine that can be used to perform off-Hugoniot measurements at higher pressures than are accessible through magnetically driven planar geometries.Current pulse shaping techniques, originally developed for planar dynamic material experiments on the Z-machine [M. K. Matzen et al., Phys. Plasmas 12, 055503 (2005)], are adapted to the design of controlled cylindrical liner implosions. By driving these targets with a current pulse shape that prevents shock formation inside the liner, shock heating is avoided along with the corresponding decrease in electrical conductivity ahead of the magnetic diffusion wave penetrating the liner. This results in an imploding liner with a significant amount of its mass in the solid phase and at multi-megabar pressures. Pressures in the solid region of a shaped pulse driven beryllium liner fielded on the Z-machine are inferred to 5.5 Mbar, while simulations suggest implosion velocities greater than 50kms-1. These solid liner experiments are diagnosed with multi-frame monochromatic x-ray backlighting which is used to infer the material density and pressure. This work has led to a new platform on the Z-machine that can be ...

Implosion dynamics and K-shell x-ray generation in large diameter stainless steel wire array Z pinches with various nesting configurations

Nested stainless steel wire array variations were investigated on the 20MA Z machine [R. B. Spielman et al., Phys. Plasmas 5, 2105 (1998)]. In order to reach experimentally observed electron temperatures near 3.8keV and excite the K shell, these ∼6.7keV photon energy x-ray sources must be of large initial diameter (45–80mm) which poses a concern for magnetic Rayleigh–Taylor instability growth. We discuss the implosion dynamics in these large diameter wire arrays, including an analysis of the ablation phase indicating that the prefill material is snowplowed at large radius. Nested array configurations with various mass and radius ratios are compared for instability mitigation and K-shell scaling. Degradation of the K-shell x-ray power and yield was observed for shots that did not have simultaneous implosion of the outer and inner wire arrays. Shots that were designed per this constraint exhibited K-shell yield scaling consistent with the model of J. W. Thornhill et al. [IEEE Trans. Plasma Sci. 34, 2377 (20...

Shock compression of hydrocarbon polymer foam using molecular dynamics

Organic polymers and nanocomposites are increasingly being subjected to extreme environments. Molecular-scale modeling of these materials offers insight into failure mechanisms and response. In previously published work, we used classical molecular dynamics (MD) and density functional theory (DFT) simulations to determine the principal shock Hugoniot for two hydrocarbon polymers, polyethylene (PE) and poly(4-methyl-1-pentene) (PMP). DFT was in excellent agreement with experiment, and one of four classical MD potentials, ReaxFF, was found to be suitable for studies up to 50 GPa. Here, we extend these results to include low-density polymer foams using NEMD techniques. We find good quantitative agreement with both experiment and hydrocode simulations. Further, we have measured local temperatures to investigate the formation of hot spots and polymer dissociation near foam voids.

Wavelength-dependent measurements of optical-fiber transit time, material dispersion, and attenuation

A new, to our knowledge, method for measuring the wavelength dependence of the transit time, material dispersion, and attenuation of an optical fiber is described. We inject light from a 4-ns rise-time pulsed broadband flash lamp into fibers of various lengths and record the transmitted signals with a time-resolved spectrograph. Segments of data spanning a range of approximately 3000 A are recorded from a single flash-lamp pulse. Comparison of data acquired with short and long fibers enables the determination of the transit time and the material dispersion as functions of wavelength dependence for the entire recorded spectrum simultaneously. The wavelength-dependent attenuation is also determined from the signal intensities. The method is demonstrated with experiments using a step-index 200-mum-diameter SiO(2) fiber. The results agree with the transit time determined from the bulk glass refractive index to within ?0.035% for the visible (4000-7200-A) spectrum and 0.12% for the UV (2650-4000-A) spectrum and with the attenuation specified by the fiber manufacturer to within ?10%.

Probing off-Hugoniot states in Ta, Cu, and Al to 1000 GPa compression with magnetically driven liner implosions

We report on a new technique for obtaining off-Hugoniot pressure vs. density data for solid metals compressed to extreme pressure by a magnetically driven liner implosion on the Z-machine (Z) at Sandia National Laboratories. In our experiments, the liner comprises inner and outer metal tubes. The inner tube is composed of a sample material (e.g., Ta and Cu) whose compressed state is to be inferred. The outer tube is composed of Al and serves as the current carrying cathode. Another aluminum liner at much larger radius serves as the anode. A shaped current pulse quasi-isentropically compresses the sample as it implodes. The iterative method used to infer pressure vs. density requires two velocity measurements. Photonic Doppler velocimetry probes measure the implosion velocity of the free (inner) surface of the sample material and the explosion velocity of the anode free (outer) surface. These two velocities are used in conjunction with magnetohydrodynamic simulation and mathematical optimization to obtain the current driving the liner implosion, and to infer pressure and density in the sample through maximum compression. This new equation of state calibration technique is illustrated using a simulated experiment with a Cu sample. Monte Carlo uncertainty quantification of synthetic data establishes convergence criteria for experiments. Results are presented from experiments with Al/Ta, Al/Cu, and Al liners. Symmetric liner implosion with quasi-isentropic compression to peak pressure ∼1000 GPa is achieved in all cases. These experiments exhibit unexpectedly softer behavior above 200 GPa, which we conjecture is related to differences in the actual and modeled properties of aluminum.

Subignition fusion yields generated by fast heating of compressed deuterium–tritium and break-even scaling

A simple model is presented to calculate the fusion yield from the fast heating of compressed deuterium–tritium (DT). The model is applicable when the fusion yield is small enough to neglect self-heating. Since 80% of the fusion yield escapes as 14 MeV neutrons, self-heating is small as long as the fusion yield is less than the deposited energy, i.e., the fusion gain, Q⩽1. We show the model is in good agreement with detailed numerical simulations when this condition is satisfied, as will be the case for fast ignition experiments in the near future. The model is used to calculate the fast heating fusion yields as a function of the important parameters such as the fuel density, deposited energy, and pulse length. The model is also used to obtain the minimum energy, Eq1, and the fuel diameter, d, necessary to obtain Q=1, which is given approximately by the scaling laws Eq1=15.3(ρ/100)−1.5 kJ, and d=16.5(ρ/100)−0.84 μm, where ρ is the DT density. These scaling laws should help the design of Q=1 experiments.

Exploring magnetized liner inertial fusion with a semi-analytic model

In this paper, we explore magnetized liner inertial fusion (MagLIF) [S. A. Slutz et al., Phys. Plasmas 17, 056303 (2010)] using a semi-analytic model [R. D. McBride and S. A. Slutz, Phys. Plasmas 22, 052708 (2015)]. Specifically, we present simulation results from this model that: (a) illustrate the parameter space, energetics, and overall system efficiencies of MagLIF; (b) demonstrate the dependence of radiative loss rates on the radial fraction of the fuel that is preheated; (c) explore some of the recent experimental results of the MagLIF program at Sandia National Laboratories [M. R. Gomez et al., Phys. Rev. Lett. 113, 155003 (2014)]; (d) highlight the experimental challenges presently facing the MagLIF program; and (e) demonstrate how increases to the preheat energy, fuel density, axial magnetic field, and drive current could affect future MagLIF performance.