Hostname: page-component-78c5997874-t5tsf Total loading time: 0 Render date: 2024-10-31T07:00:37.024Z Has data issue: false hasContentIssue false

Shock Hugoniot Data for Water up to 5 Mbar Obtained with Quartz Standard at High-Energy Laser Facilities

Published online by Cambridge University Press:  01 January 2024

D. Mancelli*
Affiliation:
Universitè de Bordeaux, CNRS, CEA, CELIA, UMR 5107, F-33405 Talence, France Donostia International Physics Center (DIPC), Donostia/San Sebastián, Spain
I. Errea
Affiliation:
Donostia International Physics Center (DIPC), Donostia/San Sebastián, Spain Centro de Física de Materiales (CSIC-UPV/EHU), Donostia/San Sebastián, Spain Fisika Aplikatua Saila, Gipuzkoako Ingeniaritza Eskola, University of the Basque Country (UPV/EHU), Donostia/San Sebastián, Spain
A. Tentori
Affiliation:
Universitè de Bordeaux, CNRS, CEA, CELIA, UMR 5107, F-33405 Talence, France
O. Turianska
Affiliation:
Universitè de Bordeaux, CNRS, CEA, CELIA, UMR 5107, F-33405 Talence, France
H. Larreur
Affiliation:
Universitè de Bordeaux, CNRS, CEA, CELIA, UMR 5107, F-33405 Talence, France
K. Katagiri
Affiliation:
Graduate School of Engineering, Osaka University, Suita, Osaka 565-0871, Japan Institute of Laser Engineering, Osaka University, Suita, Osaka 565-0871, Japan
N. Ozaki
Affiliation:
Graduate School of Engineering, Osaka University, Suita, Osaka 565-0871, Japan Institute of Laser Engineering, Osaka University, Suita, Osaka 565-0871, Japan
N. Kamimura
Affiliation:
Graduate School of Engineering, Osaka University, Suita, Osaka 565-0871, Japan
D. Kamibayashi
Affiliation:
Graduate School of Engineering, Osaka University, Suita, Osaka 565-0871, Japan
K. Ishida
Affiliation:
Graduate School of Engineering, Osaka University, Suita, Osaka 565-0871, Japan
H. Ogura
Affiliation:
Graduate School of Engineering, Osaka University, Suita, Osaka 565-0871, Japan
K. Kawasaki
Affiliation:
Institute of Laser Engineering, Osaka University, Suita, Osaka 565-0871, Japan
Y. Maeda
Affiliation:
Institute of Laser Engineering, Osaka University, Suita, Osaka 565-0871, Japan
Y. Hironaka
Affiliation:
Institute of Laser Engineering, Osaka University, Suita, Osaka 565-0871, Japan
K. Shigemori
Affiliation:
Institute of Laser Engineering, Osaka University, Suita, Osaka 565-0871, Japan
K. Batani
Affiliation:
IPPLM Institute of Plasma Physics and Laser Microfusion, Warsaw, Poland
G. Schaumann
Affiliation:
GSI Helmholtzzentrum für Schwerionenforschung GmbH-Planckstraß1, 64291 Darmstadt, Germany
O. Rosmej
Affiliation:
GSI Helmholtzzentrum für Schwerionenforschung GmbH-Planckstraß1, 64291 Darmstadt, Germany
P. Neumayer
Affiliation:
GSI Helmholtzzentrum für Schwerionenforschung GmbH-Planckstraß1, 64291 Darmstadt, Germany
B. Zielbauer
Affiliation:
GSI Helmholtzzentrum für Schwerionenforschung GmbH-Planckstraß1, 64291 Darmstadt, Germany
A. S. Martynenko
Affiliation:
GSI Helmholtzzentrum für Schwerionenforschung GmbH-Planckstraß1, 64291 Darmstadt, Germany Joint Institute for High Temperatures, RAS, 13 Bd.2 Izhorskaya st., Moscow 125412, Russia
E. D. Filippov
Affiliation:
Joint Institute for High Temperatures, RAS, 13 Bd.2 Izhorskaya st., Moscow 125412, Russia Institute of Applied Physics, RAS, 46 Ulyanov Str., 603950 Nizhny Novgorod, Russia
S. Pikuz
Affiliation:
Joint Institute for High Temperatures, RAS, 13 Bd.2 Izhorskaya st., Moscow 125412, Russia National Research Nuclear University “MEPhI”, 31 Kashirskoe shosse, Moscow 115409, Russia
D. Batani
Affiliation:
Universitè de Bordeaux, CNRS, CEA, CELIA, UMR 5107, F-33405 Talence, France National Research Nuclear University “MEPhI”, 31 Kashirskoe shosse, Moscow 115409, Russia
*
Correspondence should be addressed to D. Mancelli; donaldi.mancelli@u-bordeaux.fr
Rights & Permissions [Opens in a new window]

Abstract

In this work, we present experimental results on the behavior of liquid water at megabar pressure. The experiment was performed using the HIPER (High-Intensity Plasma Experimental Research) laser facility, a uniaxial irradiation chamber of GEKKO XII (GXII) at the Institute of Laser Engineering (ILE), and the PHELIX at GSI (GSI Helmholtz Centre for Heavy Ion Research), a single-beam high-power laser facility, to launch a planar shock into solid multilayered water samples. Equation-of-state data of water (H2O are obtained in the pressure range 0.50–4.6 Mbar by tuning the laser-drive parameters. The Hugoniot parameters (pressure, density, etc.) and the shock temperature were simultaneously determined by using VISAR and SOP as diagnostic tools and quartz as the standard material for impedance mismatch experiments. Finally, our experimental results are compared with hydrodynamic simulations tested with different equations of state, showing good compatibility with tabulated SESAME tables for water.

Type
Research Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
Copyright © 2021 D. Mancelli et al.

1. Introduction

Equation of state (EOS) of matter at extreme thermodynamic states is of relevance for several fields in physics, especially in laboratory astrophysics and in inertial confinement fusion (ICF), where precise knowledge of the EOS is required. Water is predicted to be abundant in the outer planets of the Solar System, as well as in many of the recently discovered exoplanets. In particular, it is one of the principal ingredients of the mantles of giant planets such as Uranus, Neptune, and Jupiter [Reference Mazevet, Licari, Chabrier and Potekhin1]. The observation of large and asymmetric magnetic fields in these planets [Reference Ness, Acuña and Behannon2Reference Chen, Na and Curry4] indicated that the mantle is the origin of the field. As the dynamo theory requires the presence of a conductive material, it was suggested that one or all of the main ingredients of the mantle (ammonia, water, and methane, i.e., carbon [Reference Wang, Zhang and Jiang5, Reference Jakubowska, Mancelli and Benocci6]) experience a phase transition to a conducting state. Pioneering theoretical work has been done calculating the properties of the superionic phase of water at planetary conditions [Reference Cavazzoni, Chiarotti, Scandolo, Tosatti, Bernasconi and Parrinello7Reference French, Mattsson, Nettelmann and Redmer9]. Therefore, water at very high pressures has been extensively studied in recent works [Reference Mattsson and Desjarlais8, Reference Kimura, Ozaki and Sano10Reference Nagayama, Mori, Shimada and Nakahara14]. Recently [Reference Millot, Hamel and Rygg15, Reference Millot, Coppari and Rygg16], the superionic phase of water was detected experimentally along the Uranus and Neptune isentropes at about 1.5 Mbar by laser shocking water samples precompressed up to 28 kbar. The new superionic phase is predicted to span the pressure range of 1.5 to 6 Mbars.

Superionic water is a solid system that has high ionic conductivity well below the melting temperature. Whether H2O in planetary interiors is in the superionic or metallic state is of great importance for understanding the source of the planetary magnetic field.

Laser-driven shock compression is a useful technique for generating high pressure and temperature conditions similar to the inner part of Neptune and Uranus. However, measuring P − ρ − T simultaneously is challenging, particularly temperature, as this requires absolute intensity measurements which are not necessary for the determination of the pressure or the density. In this work, we conducted laser-driven shock experiments on H2O samples up to 4.6 Mbar; along with pressure and density, we evaluated the temperature from measured reflectivity and thermal emission of the shocked sample.

The impedance mismatch (IM) method is widely used for the determination of Hugoniot of sample material. Given that one knows Hugoniot of a “standard” material that is used as a reference, the Rankine–Hugoniot (RH) [Reference Zel’dovich and Raizer17] set of equations, expressing the conservation laws, can be used to relate the experimentally measured shock velocity (D s ) in the standard material before the shock reaches the standard-sample interface in the sample after it passes through the interface. In this work, we used z-cut α-quartz (SiO2) as a reference material [Reference Knudson and Desjarlais18, Reference Brygoo, Millot and Loubeyre19] which, at ambient pressure, is transparent to visible light and becomes reflective at pressures above 1 Mbar [Reference Laudernet, Clérouin and Mazevet20]. As diagnostic tools, we used the velocity interferometer system for any reflector (VISAR) [Reference Celliers, Bradley, Collins, Hicks, Boehly and Armstrong21Reference Barker and Schuler25] which allowed us to make a precise measurement of shock wave parameters and characterize the EOS of water. The RH set provides information of P and ρ, but not of T, whereas T is also an important thermodynamic parameter. Indeed, in this work, we also measured the temperature using streaked optical pyrometry (SOP) diagnostic data.

2. Experimental Setup

2.1. GEKKO XII-ILE Laser Facility

The first experiment was carried out on the HIPER (High-Intensity Plasma Experimental Research) laser facility, a uniaxial irradiation chamber of the GEKKO XII (GXII) laser at the Institute of Laser Engineering (ILE), Osaka University. Up to 12 beams of neodymium-doped yttrium aluminum garnet (Nd:YAG) laser, frequency tripled (351 nm), were used in the experiment. The laser pulse temporal profile was approximately square shape in time with full width at half maximum (FWHM) 2.5 ns with typical rise and fall time of 100 ps each. The focal spot diameter was 600 µm flat top. SSD (smoothing by spectral dispersion technique) was applied to smooth out the beams. Kinoform phase plates were also used to achieve uniform irradiation, resulting in good shock front planarity.

Figure 1 shows the experimental setup and target configuration. A typical target assembly consisting of 10 µm parylene (C8H8) will be referred to as CH/200 nm Au/100 µm quartz (with AR@1054&532 both sides) nominally, ∼500 µm water/100 µm quartz (with AR@1054&532 both sides). The laser initially hits the 0.1 µm-thick Al layer, coated over the CH layer, to avoid laser shine through. The CH layer is the actual ablator, and quartz is the pusher layer. The 200 nm gold layer was placed to stop any X-rays from the plasma corona and avoid preheating of water. Water cells were produced at the Technical University of Darmstadt (Germany) target fabrication department. Due to the low Z ablator and the low laser intensity, X-ray radiation is low and characterized by low photon energy. The Au layer is capable to stop the X-ray radiation; thus, preheating of water is negligible.

Figure 1: Experimental setup: sketch of the multilayered target design in planar geometry with its associated thicknesses used in our experiment. Drive laser from top. From bottom, the VISAR and SOP measure the shock velocity and the self-emission of the shocked sample on the rear side.

The primary diagnostics were the VISAR and SOP. Two line-imaging interferometers (VISARs) allow to record time-resolved Doppler shift of the velocity of the fast-moving reflector and also the optical properties such as reflectivity [Reference Celliers, Bradley, Collins, Hicks, Boehly and Armstrong21Reference Yan, Liu and Zhang26]. These VISARs had different velocity-per-fringe (VPF) sensitivities to resolve the 2π-phase jump ambiguities due to the shock velocity jump at material interfaces. The sensitivity of the two VISARs was 7.523 km/s and 4.476 km/s taking into account the refractive index of SiO2. The VISAR probe beam was an injection-seeded Q-switched YAG laser. The pulse duration of the probe was ∼10 ns at the FWHM, and its wavelength was 532 nm. The postprocessing of the VISAR raw data obtained in the experiment allows determining [Reference Neutrino27] the fringe position to 10% of a fringe, while the multiple fringe shifts allow the precision of the shock velocity measurements to be a few percent. Our VISAR analysis showed uncertainties in measured D s of the order of 3%.

To measure the shock temperature, we used SOP [Reference Miller, Boehly and Melchior28Reference He, Jia and Zhang30]. The self-emission of the shocked sample at wavelength 450 nm with 38 nm bandwidth was recorded spatially and temporally resolved, using an absolutely calibrated SOP.

2.2. PHELIX-GSI Laser Facility

The second experiment was conducted at the GSI facility using the PHELIX laser, a flash-lamp-pumped Nd:glass laser utilizing the second harmonic at wavelength 527 nm. The spatial profile was flat top with spot 350 µm FWHM obtained by an appropriate phase plate [Reference Koenig, Faral, Boudenne, Batani, Benuzzi and Bossi31], and the temporal profile was top hat with a duration of 3.5 ns (FWHM). The laser was focused onto a multilayered target with intensity I 5.49 × 1012 W/cm2 to 2.52 × 1013 W/cm2. The ablation pressure was approximately 12 Mbar, was generated in our plastic ablator (parylene with gross chemical formula C8H8), and was estimated from the well-known scaling laws [Reference Mora32Reference Chen, Hu and Zhou34] for given intensities.

(1) P = 8.6 I 10 14 2 / 3 λ 2 / 3 A 2 Z 1 / 3 ,

where A and Z are the atomic mass number and the atomic number of the target material: I is in W/cm2, P in Mbar, and λ in μm. A schematic of the PHELIX experimental setup is shown in Figure 2. Targets were also produced at the Technical University of Darmstadt (Germany) target fabrication department.

Figure 2: Experimental setup: a sketch of the multilayered target design at the PHELIX laser facility in planar geometry with its associated thicknesses used in our experiment. Drive laser from the top. Bottom: the VISAR and SOP diagnostics which measure the shock velocity and the self-emission of the shocked target on the rear side.

The target configuration and the experimental setup are presented in Figure 2. The sample consists of 15 µm C8H8 (parylene)/5 µm epoxy/10 µm Al/7 µm epoxy/100 µm quartz (with AR@1054&527 both sides) nominally, ∼500 µm water/100 µm quartz (with AR@1054&527 both sides). The laser initially hits the 0.1 µm-thick Al layer, coated over the CH layer, to avoid direct laser shine through. Epoxy was used to glue CH/Al and Al/quartz interfaces. The CH layer is the actual ablation layer, and quartz is the pusher layer. The 10 µm Al layer was placed to stop any X-rays from the plasma corona and avoid preheating of the quartz/water layers. The VISAR laser had a wavelength of 660 nm, and the associated sensitivities were 1.285 km/s/f and 4.7 km/s/f, respectively. Also, SOP was looking at the self-emission of the shocked target; however, here, we report only VISAR results and compare with findings on water samples from the GEKKO XII laser facility.

3. Impedance Mismatch Method: Single-Shock Data

The impedance-matching method [Reference Wang, Zhang and Jiang5, Reference Zel’dovich and Raizer17Reference Brygoo, Millot and Loubeyre19, Reference Marshall, Lazicki and Erskine35] was used to estimate the shock state in water after passing the quartz/water interface; an illustrative method is shown in Figure 3. Because of the impedance mismatch at the SiO2/H2O interface, the shock wave produced a transmitted shock into H2O and a reflected rarefaction wave propagating back into quartz. In the rarefaction wave, the shock-compressed quartz undergoes isentropic release until its pressure and particle velocity match those of shocked water. The IM method requires precise knowledge of Hugoniot and release behavior of the standard reference material (quartz in our case) and the Rankine–Hugoniot (RH) jump relations, which are derived from the conservation laws, mass, momentum, and energy to close the system and derive all the remaining thermodynamic parameters:where ρ,Ds,Up,P, and ε denote the density, shock velocity, particle velocity, pressure, and internal energy behind the shock. The initial states are denoted with subscript 0. The third equation plotted in the (P − U) plane gives the so-called Rayleigh line of the material. Figure 3 illustrates such a method, deriving the shock pressure in water and the particle velocity. The measurement of D s Q z gives the shock state P Q z , U p Q z . The crossing of the isentropic release of the quartz standard with the Rayleigh line of the water sample determines the shock state in water P H 2 O , U p H 2 O .

Figure 3: P − Up graphical construction of the impedance mismatch procedure for the highest pressure shot of the quartz standard. The shocked state of H2O is constrained to remain on the incident Rayleigh line (blue dashed line) defined for D s Q z = 21.8 km / s for this case. The solid red and dashed lines are the Hugoniot and the Rayleigh line of (SiO2), respectively, associated with the release curve calculated for D s Q z = 21.8 km / s . The crossing of the Rayleigh line with the release isentrope path (blue square) determines the IM point and provides the experimental value of P H 2 O , U p H 2 O which may or may not lay on the shock polar for water, i.e., may (or may not) be in agreement with the used theoretical EOS.

(2) ρ ρ 0 = D s D s U p ,
(3) P P 0 = ρ 0 D s U p ,
(4) ε ε 0 = 1 2 P + P 0 1 ρ 0 1 ρ .

4. Results and Discussion

Figure 4 presents the raw data obtained with the VISAR for the shot SID-43058. The arrival time of the shock at the Au/SiO2 interface front side defines the time zero (t 0 = 0). Then, while the shock front enters into SiO2, it becomes highly reflective, and its velocity gradually decreases. It is indicated in the raw image with a dotted line as t 0. At t 1 = 7.21 ns, the shock front strikes the H2O layer as shown in Figure 4. The velocity profile and the reflectivity obtained from Figure 4 are shown in Figure 5. The transmitted shock in H2O also becomes reflective. The shock velocity in SiO2, used in the impedance mismatch analysis, was obtained by analysing VISAR images first by using shock chronometry, i.e., the average shock velocity in SiO2 was given by D s Q z = 100 μ m / t 1 t 0 . Then, the analysis is refined by looking at the displacement of the fringe shift and averaging the shock velocities from the two VISAR records taking into account the temporal resolution of the streak camera and the sensitivities of each VISAR. The IM method is applied using the instantaneous value of the shock velocity just before t 1.

Figure 4: Experimentally obtained VISAR images for the shot SID-43058.

Figure 5: Top: evolution of the shock velocity in time for SiO2 and H2O extracted from Figure 4: VISAR-1 (black) and VISAR-2 (yellow). Bottom: the reflectivity vs. time. This shows an increase after the shock enters the SiO2 layer owing to the reverberating wave in SiO2 (at about t ≈ 3-4 ns), and a second increase appears when the shock passes the SiO2-H2O interface at time t 1. Shock velocities are listed in Table 1 (SID-43058).

Table 1: Hugoniot data from the impedance mismatch method with α-quartz as a reference material. The velocity of shock in quartz DQz and water D H 2 O was used in the IM analysis to obtain the particle velocity U p H 2 O , pressure P H 2 O , and density ρ H 2 O on water Hugoniot. The compressibility ρ / ρ 0 H 2 O was calculated by dividing ρ H 2 O by the initial density. Water target was studied at two different high-power laser facilities.

To determine D s H 2 O , we instead could rely only on the VISAR images. The experimental results were supported with numerical simulation utilizing MULTI-1D, which is a one-dimensional radiation hydrodynamic code [Reference Ramis, Schmalz and Meyer-Ter-Vehn36]. Simulations show that, in the present experiment, D s Q z decays. We observe a decay in the case of D s H 2 O . Detailed interpretation of simulations will be shown in Section 6.

Getting Hugoniot points for H2O requires that one knows the release isentrope (RI) of SiO2 from its shocked state. This study uses an EOS model to calculate the release isentrope for each shock state in SiO2 utilizing an in-house Python script. The IM method provides the Hugoniot state of H2O as the locus where the release isentrope of SiO2 crosses the Rayleigh line P = ρ 0 H 2 O D s H 2 O U p of water.

Such a method is shown in Figure 3, with the Hugoniot and the RI of SiO2 (black line and red dotted curve, respectively). Dashed lines denote the Rayleigh lines of SiO2 and water (red and blue), respectively, as well as their associated Hugoniot. The Hugoniot curve is derived from the SESAME database [37] and the QEOS [Reference More, Warren, Young and Zimmerman38], tabulated data 7385 for quartz and 7154 for water. The obtained shocked states of water are compared with other available works and are shown in Figure 6 in the P − ρ plane and in Figure 7 in the P − Up plane.

Figure 6: Water Hugoniot data in the P − ρ plane. Squared data points: our experimental points. The lines correspond to tabulated data from SESAME (SESAME 7150, 7153, and 7154 and P 0 = 10 kbar) [37]. Density functional theory-molecular dynamics (DFT-MD) calculations of Hugoniot of water form French et al. [Reference French, Mattsson, Nettelmann and Redmer9, Reference French, Desjarlais and Redmer39].

Figure 7: Water Hugoniot data in the P − Up plane. Our data showed good agreement with the SESAME models available for water.

5. Temperature Measurement

The shock front emits thermal radiation, and which, in our experiment, was collected by the SOP [Reference Miller, Boehly and Melchior28, Reference Bossi, Hall and Mahdieh44]. The absolute spectral radiance of the shock front was measured at a wavelength of 450 nm with bandwidth 38 nm. Assuming a grey body spectrum, the spectral radiance I (450 nm) is given by Planck’s black body emissivity corrected for the absorptivity of the material (ε = 1 − R) calculated using the reflectivity R of the shock front measured by VISARs since the light is either absorbed or reflected. Then,

(5) I λ , T = A λ 2 h c 2 λ 5 1 e h c / λ k B T 1 .

The temperature was determined using this formula (where h, c, and kB are Planck’s constant, speed of light, and Boltzmann’s constants, respectively). A(λ) is a parameter taking into account the reflectivity, the sensitivity of the streak, the geometry, and the collection optics path.

Solving equation (5) with respect to temperature, we get

(6) T = h c λ k B ln 2 h c 2 A λ / I λ λ 5 + 1 .

SOP calibration has been made in situ using the quartz standard (GEKKO XII) [Reference Guarguaglini, Hernandez and Okuchi42, Reference Sano, Ozaki and Sakaiya45]. For the temperature estimation, we used the mean reflectivity of the two VISAR records which were taken simultaneously for each shot. The results of SOP temperature estimation are shown in Figure 8 where we report two shots at different laser energies with the associated temperature profile in time (Figure 8(a)).

Figure 8: (a) SOP result for shot SN43045 at a laser energy of 1.1 kJ, the highest energy shot obtained in our experiment, and (b) temperature estimation taking into account the mean reflectivity of VISAR-1 and VISAR-2. The initial luminosity (t ≤ 4 ns) is due to scattered laser light.

The temperature results as a function of pressure reaching the water interface are shown in Figure 9 together with results from Lyzenga et al. [Reference Lyzenga, Ahrens, Nellis and Mitchell41], Kimura et al. [Reference Kimura, Ozaki and Sano10], and Guarguaglini et al. [Reference Guarguaglini, Hernandez and Okuchi42] and Hugoniot tabulated curves based on SESAME models (SESAME 7150, 7153, and 7154 and P 0 = 10 kbar) and ab initio Hugoniot from French et al. [Reference French, Desjarlais and Redmer39].

Figure 9: Principal Hugoniot data in the P − T plane obtained in this study from the GEKKO XII-ILE laser facility, compared with Lyzenga et al. [Reference Lyzenga, Ahrens, Nellis and Mitchell41], Kimura et al. [Reference Kimura, Ozaki and Sano10], Guarguaglini et al. [Reference Guarguaglini, Hernandez and Okuchi42], and Batani et al. [Reference Batani, Jakubowska and Benuzzi-Mounaix43]. Also shown are SESAME models for three different tabulated EOS data at slightly different initial densities.

6. Hydro Simulations

1D radiative hydrodynamic simulations were performed with MULTI-1D [Reference Ramis, Schmalz and Meyer-Ter-Vehn36] to comprehend our experimental results. The laser temporal profile of the pulse was flat top in time with a plateau duration of 2.5 ns at FWHM and rise and fall times of 0.1 ns. In the simulation, we utilized the SESAME tables of the following materials which consist of our target SESAME table 7770 for parylene [37], SESAME table 2700 for gold [37], SESAME 7385 for quartz [37], and SESAME tables 7150, 7153, and 7154 for water [37]. Concerning the water EOS, we compared different models of EOS tables coming from the SESAME database [37], from QEOS [Reference More, Warren, Young and Zimmerman38], and from FEOS [Reference Faik, Tauschwitz and Iosilevskiy46], in all cases setting the initial density at ρ 0 = 0.98 (g/cm3) in order to check the validity of the EOS data. The density and pressure map, reproducing the shot SID-43058 obtained from 1D simulation, is shown in Figures 10(a) and 10(b), respectively. Figure 10(c) shows the change of the maximum density postprocessed from Figure 10(a) at interval times, and Figure 10(d) shows the change of pressure along the shock front of the simulated target top (Figure 10(b)), respectively. The use of 1D simulation to interpret the experimental results is confirmed because of the use of appropriate phase plate resulting in a large focal spot (∼600 µm). In addition to this, the justification of the 1D approximation is supported by two experimental observations; such results can be seen in Figure 4 where good planarity of the shock is observed.

Raw data from VISAR images (see Figure 4) indicate that the shock breakout is quite flat both at metal/SiO2 and at SiO2/H2O interfaces. At this point, we must note that 2D effects will result in curvature, initially affecting the edges of the shock front, yet gradually advancing to the center. However, in our experimental results, the curvature of the shock front is not observed suggesting that the 2D effects in hydrodynamics can be neglected. In our case, the velocity of the shock is maintained and decaying quite slowly, and actually (within error bars), the decay obtained by analysing the VISAR images is compatible with the results of 1D simulations. The reduction of shock pressure and velocity during propagation is caused by two phenomena: (a) the release wave from the target front side catching up the traveling shock and (b) 2D effects in shock front propagation. In fact, due to the fair agreement of conducted experiment with 1D simulations, we conclude that (b) is not important. Indeed, a much faster decay of shock pressure and velocity would be expected if 2D effects were important.

Figure 10: Laser drive beam hits the target from the top as indicated by the arrow symbol; the constituent layers are also indicated with the right order (CH/Au/SiO2/water). (a) Hydrodynamic progression of the density map obtained from MULTI-1D reproducing the experimental results of shot SID-43058 with nominal laser intensity 1.4 × 1013 W/cm2; for this particular shot, 3 beams/12 were used for the main drive delivering a total of 252 J on the target. With the above laser parameter, the simulation reproduces well the shock breakout time in the SiO2-H2O inner interface. (b) Variation of the pressure map for the same shot conditions. (c) Time evolution of pressure and density forms MULTI simulation for shot SID-43058. (d) Change of pressure at the shock front up to the maximum simulation time.

7. Conclusion

In summary, we obtained EOS data of water along the principal Hugoniot up to 5 Mbar. Water samples, contained within a multilayered water cell, were dynamically compressed in planar geometry using the high-power laser facilities GEKKO XII (ILE) and PHELIX (GSI). Utilizing quartz as a standard material in both experimental campaigns and the main diagnostics such as VISAR/SOP substantially reduced experimental errors in optimized experimental conditions. The impedance mismatching analysis allowed to verify that P, ρ, (εε 0, and the Hugoniot data are in fair agreement with those predicted by SESAME table 7150, and on the contrary, they show a significant difference concerning the Hugoniot curve calculated using DFT-MD simulations. Also, our experimental outcome showed good agreement with simulations performed with the radiative hydro code MULTI-1D using the SESAME tabulated EOS, the QEOS model, and the FEOS model, a modified version of the QEOS. The agreement with 1D simulations shows indeed that, in our experimental setup, 2D effects in hydrodynamics are negligible, a result which mainly depends on the use of laser focal spots. For a few shots, we could also measure the temperature of shocked material using calibrated SOP diagnostics. Our data confirm previous experimental results and show that, in the pressure range up to 4.6 Mbar, water is in a reflective state.

Data Availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Disclosure

This work has been carried out within the framework of the EUROfusion Consortium. The views and opinions expressed herein do not necessarily reflect those of the European Commission. Our findings were presented at the Satellite Meeting of the 47th European Conference on Plasma Physics.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors would like to acknowledge the support of the laser technical team at GEKKO XII (ILE) and PHELIX (GSI). This work received funding from the Euratom Research and Training Program 2014–2018 and 2019-2020 (Grant agreement no. 633053). The involved teams were operated within the framework of the Enabling Research Project ENR-IFE19.CEA-01, Study of Direct Drive and Shock Ignition for IFE: Theory, Simulations, Experiments, Diagnostics Development. JIHT RAS team members were supported by the Ministry of Science and Higher Education of the Russian Federation (State Assignment no. 075-00460-21-00).

References

Mazevet, S., Licari, A., Chabrier, G., and Potekhin, A. Y., “Ab initio based equation of state of dense water for planetary and exoplanetary modeling,” Astronomy & Astrophysics, vol. 621, p. A128, 2019.10.1051/0004-6361/201833963CrossRefGoogle Scholar
Ness, N. F., Acuña, M. H., Behannon, K. W. et al., “Magnetic fields at Uranus,” Science, vol. 233, no. 4759, pp. 8589, 1986.10.1126/science.233.4759.85CrossRefGoogle ScholarPubMed
Connerney, J. E. P., Acuña, M. H., and Ness, N. F., “The magnetic field of Neptune,” Advances in Space Research, vol. 12, no. 8, pp. 239248, 1992.10.1016/0273-1177(92)90394-DCrossRefGoogle Scholar
Chen, Z., Na, X., Curry, C. B. et al., “Observation of a highly conductive warm dense state of water with ultrafast pump-probe free-electron-laser measurements,” Matter and Radiation at Extremes, vol. 6, no. 5, Article ID 054401, 2021.10.1063/5.0043726CrossRefGoogle Scholar
Wang, P., Zhang, C., Jiang, S. et al., “Density-dependent shock hugoniot of polycrystalline diamond at pressures relevant to ICF,” Matter and Radiation at Extremes, vol. 6, no. 3, Article ID 035902, 2021.10.1063/5.0039062CrossRefGoogle Scholar
Jakubowska, K., Mancelli, D., Benocci, R. et al., “Reflecting laser-driven shocks in diamond in the megabar pressure range,” High Power Laser Science and Engineering, vol. 9, no. e3, 2021.10.1017/hpl.2020.38CrossRefGoogle Scholar
Cavazzoni, C., Chiarotti, G. L., Scandolo, S., Tosatti, E., Bernasconi, M., and Parrinello, M., “Superionic and metallic states of water and ammonia at giant planet conditions,” Science, vol. 283, no. 5398, pp. 4446, 1999.10.1126/science.283.5398.44CrossRefGoogle ScholarPubMed
Mattsson, T. R. and Desjarlais, M. P., “Phase diagram and electrical conductivity of high energy-density water from density functional theory,” Physical Review Letters, vol. 97, Article ID 017801, 2006.10.1103/PhysRevLett.97.017801CrossRefGoogle ScholarPubMed
French, M., Mattsson, T. R., Nettelmann, N., and Redmer, R., “Equation of state and phase diagram of water at ultrahigh pressures as in planetary interiors,” Physical Review B: Condensed Matter, vol. 79, Article ID 054107, 2009.10.1103/PhysRevB.79.054107CrossRefGoogle Scholar
Kimura, T., Ozaki, N., Sano, T. et al., “P-ρ-T measurements of H2O up to 260 GPa under laser-driven shock loading,” The Journal of Chemical Physics, vol. 142, no. 16, Article ID 164504, 2015.10.1063/1.4919052CrossRefGoogle ScholarPubMed
Goldman, N., Reed, E. J., Kuo, I.-F. W., Fried, L. E., Mundy, C. J., and Curioni, A., “Ab initio simulation of the equation of state and kinetics of shocked water,” The Journal of Chemical Physics, vol. 130, no. 12, Article ID 124517, 2009.10.1063/1.3089426CrossRefGoogle ScholarPubMed
Pan, D., Wan, Q., and Galli, G., “The refractive index and electronic gap of water and ice increase with increasing pressure,” Nature Communications, vol. 5, no. 1, p. 3919, 2014.10.1038/ncomms4919CrossRefGoogle ScholarPubMed
Dolan, D. H., Knudson, M. D., Hall, C. A., and Deeney, C., “A metastable limit for compressed liquid water,” Nature Physics, vol. 3, no. 5, pp. 339342, 2007.10.1038/nphys562CrossRefGoogle Scholar
Nagayama, K., Mori, Y., Shimada, K., and Nakahara, M., “Shock hugoniot compression curve for water up to 1 GPa by using a compressed gas gun,” Journal of Applied Physics, vol. 91, no. 1, pp. 476482, 2002.10.1063/1.1421630CrossRefGoogle Scholar
Millot, M., Hamel, S., Rygg, J. R. et al., “Experimental evidence for superionic water ice using shock compression,” Nature Physics, vol. 14, no. 3, pp. 297302, 2018.10.1038/s41567-017-0017-4CrossRefGoogle Scholar
Millot, M., Coppari, F., Rygg, J. R. et al., “Nanosecond X-ray diffraction of shock-compressed superionic water ice,” Nature, vol. 569, no. 7755, pp. 251255, 2019.10.1038/s41586-019-1114-6CrossRefGoogle ScholarPubMed
Zel’dovich, I. A. B. and Raizer, Y. P., “Physics of shock waves and high-temperature hydrodynamic phenomena,” Dover Books on Physics, Dover Publications, Mineola, NY, USA, 2002.Google Scholar
Knudson, M. D. and Desjarlais, M. P., “Adiabatic release measurements in α-quartz between 300 and 1200 GPa: characterization of α-quartz as a shock standard in the multimegabar regime,” Physical Review B: Condensed Matter, vol. 88, Article ID 184107, 2013.10.1103/PhysRevB.88.184107CrossRefGoogle Scholar
Brygoo, S., Millot, M., Loubeyre, P. et al., “Analysis of laser shock experiments on precompressed samples using a quartz reference and application to warm dense hydrogen and helium,” Journal of Applied Physics, vol. 118, no. 19, Article ID 195901, 2015.10.1063/1.4935295CrossRefGoogle Scholar
Laudernet, Y., Clérouin, J., and Mazevet, S., “Ab initiosimulations of the electrical and optical properties of shock-compressed SiO2,” Physical Review B, vol. 70, no. 16, Article ID 165108, 2004.10.1103/PhysRevB.70.165108CrossRefGoogle Scholar
Celliers, P. M., Bradley, D. K., Collins, G. W., Hicks, D. G., Boehly, T. R., and Armstrong, W. J., “Line-imaging velocimeter for shock diagnostics at the omega laser facility,” Review of Scientific Instruments, vol. 75, no. 11, pp. 49164929, 2004.10.1063/1.1807008CrossRefGoogle Scholar
Bloomquist, D. D. and Sheffield, S. A., “Optically recording interferometer for velocity measurements with subnanosecond resolution,” Journal of Applied Physics, vol. 54, no. 4, pp. 17171722, 1983.10.1063/1.332222CrossRefGoogle Scholar
Bolme, C. A. and Ramos, K. J., “Line-imaging velocimetry for observing spatially heterogeneous mechanical and chemical responses in plastic bonded explosives during impact,” Review of Scientific Instruments, vol. 84, no. 8, Article ID 083903, 2013.10.1063/1.4817307CrossRefGoogle ScholarPubMed
Barker, L. M. and Hollenbach, R. E., “Laser interferometer for measuring high velocities of any reflecting surface,” Journal of Applied Physics, vol. 43, no. 11, pp. 46694675, 1972.10.1063/1.1660986CrossRefGoogle Scholar
Barker, L. M. and Schuler, K. W., “Correction to the velocity‐per‐fringe relationship for the VISAR interferometer,” Journal of Applied Physics, vol. 45, no. 8, pp. 36923693, 1974.10.1063/1.1663841CrossRefGoogle Scholar
Yan, Z., Liu, H., Zhang, X. et al., “Dynamics of particles near the surface of a medium under ultra-strong shocks,” Matter and Radiation at Extremes, vol. 6, no. 2, Article ID 026903, 2021.10.1063/5.0030906CrossRefGoogle Scholar
Neutrino, V. T. F. A., “A light, expandable and full featured image analysis tool for research,2011, https://github.com/NeutrinoToolkit/Neutrino.Google Scholar
Miller, J. E., Boehly, T. R., Melchior, A. et al., “Streaked optical pyrometer system for laser-driven shock-wave experiments on omega,” Review of Scientific Instruments, vol. 78, no. 3, Article ID 034903, 2007.10.1063/1.2712189CrossRefGoogle ScholarPubMed
Gregor, M. C., Boni, R., Sorce, A. et al., “Absolute calibration of the omega streaked optical pyrometer for temperature measurements of compressed materials,” Review of Scientific Instruments, vol. 87, no. 11, Article ID 114903, 2016.10.1063/1.4968023CrossRefGoogle ScholarPubMed
He, Z., Jia, G., Zhang, F. et al., “Calibration and verification of streaked optical pyrometer system used for laser-induced shock experiments,” High Power Laser Science and Engineering, vol. 7, Article ID e49, 2019.10.1017/hpl.2019.27CrossRefGoogle Scholar
Koenig, M., Faral, B., Boudenne, J. M., Batani, D., Benuzzi, A., and Bossi, S., “Optical smoothing techniques for shock wave generation in laser-produced plasmas,” Physical Review E, vol. 50, no. 5, pp. R3314ߝR3317, 1994.10.1103/PhysRevE.50.R3314CrossRefGoogle ScholarPubMed
Mora, P., “Theoretical model of absorption of laser light by a plasma,” Physics of Fluids, vol. 25, no. 6, pp. 10511056, 1982.10.1063/1.863837CrossRefGoogle Scholar
Lindl, J., “Development of the indirect‐drive approach to inertial confinement fusion and the target physics basis for ignition and gain,” Physics of Plasmas, vol. 2, no. 11, pp. 39334024, 1995.10.1063/1.871025CrossRefGoogle Scholar
Chen, P., Hu, R., Zhou, H. et al., “Numerical investigation of radiation ablation and acceleration of high-density carbon foils,” Laser and Particle Beams, vol. 38, no. 4, pp. 239243, 2020.10.1017/S0263034620000336CrossRefGoogle Scholar
Marshall, M. C., Lazicki, A. E., Erskine, D. et al., “Developing quartz and molybdenum as impedance-matching standards in the 100-mbar regime,” Physical Review B, vol. 99, no. 17, Article ID 174101, 2019.10.1103/PhysRevB.99.174101CrossRefGoogle Scholar
Ramis, R., Schmalz, R., and Meyer-Ter-Vehn, J., “Multi—a computer code for one-dimensional multigroup radiation hydrodynamics,” Computer Physics Communications, vol. 49, no. 3, pp. 475505, 1988.10.1016/0010-4655(88)90008-2CrossRefGoogle Scholar
T4 Group LANL, “Eos sesame tables developed at the Los Alamos laboratory (sesame report on the los alamos equation-of-state library, report no. lalp-83-4, 1983, and report la-ur-92-3407, 1992, t4 group lanl, Los Alamos)”, Report LA-UR-92-3407, 1992.Google Scholar
More, R. M., Warren, K. H., Young, D. A., and Zimmerman, G. B., “A new quotidian equation of state (QEOS) for hot dense matter,” Physics of Fluids, vol. 31, no. 10, pp. 30593078, 1988.10.1063/1.866963CrossRefGoogle Scholar
French, M., Desjarlais, M. P., and Redmer, R., “Ab initio calculation of thermodynamic potentials and entropies for superionic water,” Physical Review E, vol. 93, Article ID 022140, 2016.10.1103/PhysRevE.93.022140CrossRefGoogle ScholarPubMed
Knudson, M. D., Desjarlais, M. P., Lemke, R. W. et al., “Probing the interiors of the ice giants: shock compression of water to 700 GPa and 3.8 (g/cm3),” Physical Review Letters, vol. 108, Article ID 091102, 2012.10.1103/PhysRevLett.108.091102CrossRefGoogle ScholarPubMed
Lyzenga, G. A., Ahrens, T. J., Nellis, W. J., and Mitchell, A. C., “The temperature of shock‐compressed water,” The Journal of Chemical Physics, vol. 76, no. 12, pp. 62826286, 1982.10.1063/1.443031CrossRefGoogle Scholar
Guarguaglini, M., Hernandez, J.-A., Okuchi, T. et al., “Laser-driven shock compression of “synthetic planetary mixtures” of water, ethanol, and ammonia,” Scientific Reports, vol. 9, no. 1, Article ID 10155, 2019.10.1038/s41598-019-46561-6CrossRefGoogle ScholarPubMed
Batani, D., Jakubowska, K., Benuzzi-Mounaix, A. et al., “Refraction index of shock compressed water in the megabar pressure range,” EPL (Europhysics Letters), vol. 112, no. 3, Article ID 36001, 2015.Google Scholar
Bossi, S., Hall, T. A., Mahdieh, M. et al., “Determination of the color temperature in laser-produced shocks,” Laser and Particle Beams, vol. 15, no. 4, pp. 485493, 1997.10.1017/S0263034600011071CrossRefGoogle Scholar
Sano, T., Ozaki, N., Sakaiya, T. et al., “Laser-shock compression and hugoniot measurements of liquid hydrogen to 55 GPa,” Physical Review B: Condensed Matter, vol. 83, Article ID 054117, 2011.10.1103/PhysRevB.83.054117CrossRefGoogle Scholar
Faik, S., Tauschwitz, A., and Iosilevskiy, I., “The equation of state package FEOS for high energy density matter,” Computer Physics Communications, vol. 227, pp. 117125, 2018.10.1016/j.cpc.2018.01.008CrossRefGoogle Scholar
Figure 0

Figure 1: Experimental setup: sketch of the multilayered target design in planar geometry with its associated thicknesses used in our experiment. Drive laser from top. From bottom, the VISAR and SOP measure the shock velocity and the self-emission of the shocked sample on the rear side.

Figure 1

Figure 2: Experimental setup: a sketch of the multilayered target design at the PHELIX laser facility in planar geometry with its associated thicknesses used in our experiment. Drive laser from the top. Bottom: the VISAR and SOP diagnostics which measure the shock velocity and the self-emission of the shocked target on the rear side.

Figure 2

Figure 3: P − Up graphical construction of the impedance mismatch procedure for the highest pressure shot of the quartz standard. The shocked state of H2O is constrained to remain on the incident Rayleigh line (blue dashed line) defined for DsQz=21.8km/s for this case. The solid red and dashed lines are the Hugoniot and the Rayleigh line of (SiO2), respectively, associated with the release curve calculated for DsQz=21.8km/s. The crossing of the Rayleigh line with the release isentrope path (blue square) determines the IM point and provides the experimental value of PH2O,UpH2O which may or may not lay on the shock polar for water, i.e., may (or may not) be in agreement with the used theoretical EOS.

Figure 3

Figure 4: Experimentally obtained VISAR images for the shot SID-43058.

Figure 4

Figure 5: Top: evolution of the shock velocity in time for SiO2 and H2O extracted from Figure 4: VISAR-1 (black) and VISAR-2 (yellow). Bottom: the reflectivity vs. time. This shows an increase after the shock enters the SiO2 layer owing to the reverberating wave in SiO2 (at about t ≈ 3-4 ns), and a second increase appears when the shock passes the SiO2-H2O interface at time t1. Shock velocities are listed in Table 1 (SID-43058).

Figure 5

Table 1: Hugoniot data from the impedance mismatch method with α-quartz as a reference material. The velocity of shock in quartz DQz and water DH2O was used in the IM analysis to obtain the particle velocity UpH2O, pressure PH2O, and density ρH2O on water Hugoniot. The compressibility ρ/ρ0H2O was calculated by dividing ρH2O by the initial density. Water target was studied at two different high-power laser facilities.

Figure 6

Figure 6: Water Hugoniot data in the P − ρ plane. Squared data points: our experimental points. The lines correspond to tabulated data from SESAME (SESAME 7150, 7153, and 7154 and P0 = 10 kbar) [37]. Density functional theory-molecular dynamics (DFT-MD) calculations of Hugoniot of water form French et al. [9, 39].

Figure 7

Figure 7: Water Hugoniot data in the P − Up plane. Our data showed good agreement with the SESAME models available for water.

Figure 8

Figure 8: (a) SOP result for shot SN43045 at a laser energy of 1.1 kJ, the highest energy shot obtained in our experiment, and (b) temperature estimation taking into account the mean reflectivity of VISAR-1 and VISAR-2. The initial luminosity (t ≤ 4 ns) is due to scattered laser light.

Figure 9

Figure 9: Principal Hugoniot data in the P − T plane obtained in this study from the GEKKO XII-ILE laser facility, compared with Lyzenga et al. [41], Kimura et al. [10], Guarguaglini et al. [42], and Batani et al. [43]. Also shown are SESAME models for three different tabulated EOS data at slightly different initial densities.

Figure 10

Figure 10: Laser drive beam hits the target from the top as indicated by the arrow symbol; the constituent layers are also indicated with the right order (CH/Au/SiO2/water). (a) Hydrodynamic progression of the density map obtained from MULTI-1D reproducing the experimental results of shot SID-43058 with nominal laser intensity 1.4 × 1013 W/cm2; for this particular shot, 3 beams/12 were used for the main drive delivering a total of 252 J on the target. With the above laser parameter, the simulation reproduces well the shock breakout time in the SiO2-H2O inner interface. (b) Variation of the pressure map for the same shot conditions. (c) Time evolution of pressure and density forms MULTI simulation for shot SID-43058. (d) Change of pressure at the shock front up to the maximum simulation time.