Sandia's Z Machine has been used to magnetically drive shockless compression of materials in a planar configuration to multi-megabar pressure levels, allowing accurate measurements of quasi-isentropic mechanical response at relatively low temperatures in the solid phase. This paper details recent improvements to design and analysis of such experiments, including the use of new data on the mechanical and optical response of lithium fluoride windows. Comparison of windowed and free-surface data on copper to 350 GPa lends confidence to the window correction method. Preliminary results are presented on gold to 500 GPa and platinum to 450 GPa; both appear stiffer than existing models.
Shock compression experiments in the few hundred GPa (multi-Mbar) regime were performed on Lithium Deuteride single crystals. This study utilized the high velocity flyer plate capability of the Sandia Z Machine to perform impact experiments at flyer plate velocities in the range of 17-32 km/s. Measurements included pressure, density, and temperature between ∼190 and 570 GPa along the Principal Hugoniot - the locus of end states achievable through compression by large amplitude shock waves - as well as pressure and density of reshock states up to ∼920 GPa. The experimental measurements are compared with density functional theory calculations, tabular equation of state models, and legacy nuclear driven results that have been reanalyzed using modern equations of state for the shock wave standards used in the experiments.
Single crystal lithium fluoride (LiF), oriented [100], was shock loaded and subsequently shocklessly compressed in two experiments at the Z Machine. Velocimetry measurements were employed to obtain an impactor velocity, shock transit times, and in-situ particle velocities for LiF samples up to ∼1.8 mm thick. A dual thickness Lagrangian analysis was performed on the in-situ velocimetry data to obtain the mechanical response along the loading path of these experiments. An elastic response was observed on one experiment during initial shockless compression from 100 GPa before yielding. The relatively large thickness differences utilized for the dual sample analyses (up to ∼1.8 mm) combined with a relative timing accuracy of ∼0.2 ns resulted in an uncertainty of less than 1% on density and stress at ∼200 GPa peak loading on one experiment and <4% on peak loading at ∼330 GPa for another. The stress-density analyses from these experiments compare favorably with recent equation of state models for LiF.
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.
We have developed the design of Thor: a pulsed power accelerator that delivers a precisely shaped current pulse with a peak value as high as 7 MA to a strip-line load. The peak magnetic pressure achieved within a 1-cm-wide load is as high as 100 GPa. Thor is powered by as many as 288 decoupled and transit-time isolated bricks. Each brick consists of a single switch and two capacitors connected electrically in series. The bricks can be individually triggered to achieve a high degree of current pulse tailoring. Because the accelerator is impedance matched throughout, capacitor energy is delivered to the strip-line load with an efficiency as high as 50%. We used an iterative finite element method (FEM), circuit, and magnetohydrodynamic simulations to develop an optimized accelerator design. When powered by 96 bricks, Thor delivers as much as 4.1 MA to a load, and achieves peak magnetic pressures as high as 65 GPa. When powered by 288 bricks, Thor delivers as much as 6.9 MA to a load, and achieves magnetic pressures as high as 170 GPa. We have developed an algebraic calculational procedure that uses the single brick basis function to determine the brick-triggering sequence necessary to generate a highly tailored current pulse time history for shockless loading of samples. Thor will drive a wide variety of magnetically driven shockless ramp compression, shockless flyer plate, shock-ramp, equation of state, material strength, phase transition, and other advanced material physics experiments.
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.
Aluminum has been used prolifically as an impedance matching standard in the multimegabar regime (1 Mbar = 100 GPa), particularly in nuclear driven, early laser driven, and early magnetically driven flyer plate experiments. The accuracy of these impedance matching measurements depends upon the knowledge of both the Hugoniot and release or reshock response of aluminum. Here, we present the results of several adiabatic release measurements of aluminum from ∼400-1200 GPa states along the principal Hugoniot using full density polymethylpentene (commonly known as TPX), and both ∼190 and ∼110 mg/cc silica aerogel standards. These data were analyzed within the framework of a simple, analytical model that was motivated by a first-principles molecular dynamics investigation into the release response of aluminum, as well as by a survey of the release response determined from several tabular equations of state for aluminum. Combined, this theoretical and experimental study provides a method to perform impedance matching calculations without the need to appeal to any tabular equation of state for aluminum. As an analytical model, this method allows for propagation of all uncertainty, including the random measurement uncertainties and the systematic uncertainties of the Hugoniot and release response of aluminum. This work establishes aluminum for use as a high-precision standard for impedance matching in the multimegabar regime.
Magnetically-driven, planar shockless-compression experiments to multi-megabar pressures were performed on tantalum samples using a stripline target geometry. Free-surface velocity waveforms were measured in 15 cases; nine of these in a dual-sample configuration with two samples of different thicknesses on opposing electrodes, and six in a single-sample configuration with a bare electrode opposite the sample. Details are given on the application of inverse Lagrangian analysis (ILA) to these data, including potential sources of error. The most significant source of systematic error, particularly for single-sample experiments, was found to arise from the pulse-shape dependent free-surface reflected wave interactions with the deviatoric-stress response of tantalum. This could cause local, possibly temporary, unloading of material from a ramp compressed state, and thus multi-value response in wave speed that invalidates the free-surface to in-material velocity mapping step of ILA. By averaging all 15 data sets, a final result for the principal quasi-isentrope of tantalum in stress-strain was obtained to a peak longitudinal stress of 330GPa with conservative uncertainty bounds of ±4.5% in stress. The result agrees well with a tabular equation of state developed at Los Alamos National Laboratory.
S hock compression exper iments in the few hundred GPa (multi - Mabr) regime were performed on Lithium Deuteride (LiD) single crystals . This study utilized the high velocity flyer plate capability of the Sandia Z Machine to perform impact experiments at flyer plate velocities in the range of 17 - 32 km/s. Measurements included pressure, density, and temperature between %7E200 - 600 GPa along the Principal Hugoniot - the locus of end states achievable through compression by large amplitude shock waves - as well as pressure and density of re - shock states up to %7E900 GPa . The experimental measurements are compared with recent density functional theory calculations as well as a new tabular equation of state developed at Los Alamos National Labs.