Technical Comments

Comment on “Insulator-metal transition in dense fluid deuterium”

See allHide authors and affiliations

Science  22 Mar 2019:
Vol. 363, Issue 6433, eaaw0969
DOI: 10.1126/science.aaw0969

Abstract

Celliers et al. (Reports, 17 August 2018, p. 677), in an attempt to reconcile differences in inferred metallization pressures, provide an alternative temperature analysis of the Knudson et al. experiments (Reports, 26 June 2015, p. 1455). We show that this reanalysis implies an anomalously low specific heat for the metallic fluid that is clearly inconsistent with first-principles calculations.

Celliers et al. (1) recently reported observation of the insulator-metal transition in deuterium at a substantially lower pressure (~200 GPa) than that reported by Knudson et al. (2) (~280 GPa). To reconcile this difference in metallization pressure, they suggest a correction to the inferred temperatures reported by Knudson et al. Their reanalysis assumes that (i) both sets of experiments—referred to as the NIF and Z experiments, respectively—sample the same first-order phase boundary; (ii) the NIF and Z pressure states are accurate and reflect the entry and exit of the coexistence region, respectively; and (iii) the boundary exhibits constant latent heat. From these assumptions, and applying a Clausius-Clapeyron analysis, they derive a correction factor ΔT = –0.49T0; that is, the temperatures in the Z experiments decrease by ~600 K for the lowest-T loading path and by nearly 900 K for the highest-T loading path.

Here, we show that the ΔT obtained in their analysis has severe thermodynamic implications, namely an anomalously low specific heat for the metallic liquid. Because the compression path is isentropic, ΔT must be such that the system follows an isentrope, as illustrated by the simplified schematic in Fig. 1. Traversal of the phase boundary isothermally (at S1 and T0) results in an increase in entropy, ΔS, governed by ΔH = T0ΔS. The ΔT (= T0T1) required to account for this ΔS, within the simplified diagram of Fig. 1, can be approximated from the specific heat at constant pressure, CP=TdSdT|P(1)If we assume that CP depends weakly on T and P away from the coexistence boundary,ΔSCPlnT0T1(2)Thus, given CP, one can determine ΔT for a given isentrope. Likewise, given ΔT, one can determine the thermodynamically consistent CP.

Fig. 1

Although ΔH is expected to vary along the phase boundary (going to zero at the critical point), we use ΔH = 2.75 kJ/g for illustration (average of the range given by Celliers et al.). In that case, at T0 = 1400 K, ΔS = 0.48 kB/atom across the phase boundary. In their reanalysis, Celliers et al. concluded that the temperature drops by ~50%, so the logarithm term in Eq. 2 is equal to 0.69. This implies CP = 0.70 kB/atom, much lower than expected on physical grounds [CP = 3.5 kB/atom for liquid lithium metal (3)].

We use tables S-IV and S-V in the supporting information for Pierleoni et al. (4) to extract a realistic value for CP in the metallic fluid, including nuclear quantum effects. A value of 2.6 kB/atom for the metallic fluid is found by differencing the hydrogen results (interpolated to constant P) on the metallic fluid branches away from the phase boundary. Using this value for CP and ΔS = 0.48 kB/atom, one finds T1/T0 = exp(–ΔS/CP) = 0.83 at 1400 K, or ΔT = 240 K, in reasonable agreement with the direct isentrope calculation shown in figure 1 (green curve) of Knudson et al. (2) and much smaller than ΔT = 700 K suggested by Celliers et al.

The above analysis is meant to be illustrative but neglects the real behavior of the isentropes close to the phase boundary. In actuality, an isentrope’s T varies with P and exhibits negative slope close to the transition, consistent with a negative Grüneisen γ. See, for example, the isentropes obtained by thermodynamic integration in figure 1 of Knudson et al. (2). To fully account for the phase boundary slope and the curvature of the isentropes near the boundary, we consider the exact expression for the entropy change along the coexistence line, derived from S(P, T):dSdT)coex=CPT[1CVCPγTBTdPdT)coex](3)where BT is the isothermal bulk modulus, CV is the specific heat at constant volume, and γ is the Grüneisen γ.

For this expression, we again take advantage of the published hydrogen values in the supporting information for Pierleoni et al. (4). Using their figure S-1 for the compressibility, along with their tables S-IV and S-V for 1500 K and 1200 K, respectively, and spanning the coexistence region in the vicinity of 1350 K from the molecular to atomic branches, we find BT = 385 GPa, CV = 5.16 kB/atom (at rs = 1.43), CP = 8.94 kB/atom (interpolating at 189 GPa), and γ = –2.3. In the metallic fluid away from the boundary, we note, in comparison to the 2.6 kB/atom found above, a substantial increase in CP on the boundary, typical of a first-order phase transition. As a consistency check on these values, we consider the general relation between CP and CV, CP = CV(1 + γ2TCV/VBT) = 8.93 kB/atom, and find very good agreement with the extracted value. Note especially that increasing the isentrope curvature through γ2 increases CP. For the slope of the coexistence boundary dP/dTcoex in the vicinity of 1350 K, we use –0.12 GPa/K, which agrees well with the coexistence lines for vdW-DF1 and vdW-DF2, including nuclear quantum effects, around that temperature (2). Regarding the phase boundary slope, we note that the NIF and Z boundary slopes are very similar, and that both are considerably steeper than the vdW-DF1, vdW-DF2, and CEIMC (4) boundaries.

Combining terms and inserting in Eq. 3, we finddSdT)coex=4.05kB/atomT(4)nearly six times larger than implied by the Celliers et al. analysis. For a value of ΔS = 0.48 kB/atom, and integrating, we estimate T1/T0 = 0.89 centered at 1350 K. Applying this to the isentrope that enters the coexistence region at 1416 K in figure 1 of Knudson et al. (2) (green curve), we estimate an exit at 1260 K, in very good agreement with the isentrope obtained there by thermodynamic integration. We conclude from this analysis that the factor of 2 reduction in T along the coexistence boundary suggested by Celliers et al. is not consistent with the equation of state for the metallic liquid, whereas the estimated temperature drops reported in Knudson et al. are consistent. Thus, the argument presented in Celliers et al. is not a viable reconciliation of the two sets of experimental results, and the apparent discrepancy remains.

References and Notes

Acknowledgments: Sandia National Laboratories is a multimission laboratory managed and operated by National Technology and Engineering Solutions of Sandia LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525. This paper describes objective technical results and analysis. Any subjective views or opinions that might be expressed in the paper do not necessarily represent the views of the U.S. Department of Energy or the United States Government. Supported by Deutsche Forschungsgemeinschaft grants SFB 652 and FOR 2440 (R.R.).
View Abstract

Navigate This Article