The Leibfried-Schlömann (LS) equation, a commonly assumed model for the pressure dependence of thermal conductivity , is tested by measurements on compressed H 2 O using a combination of the time-domain thermoreflectance method with the diamond anvil cell technique. The thermal conductivity of ice VII increases by an order of magnitude between 2 and 22 GPa, reaching ≈ 25 W m −1 K −1 . Over a large compression range of ≈4%-33%, the LS equation describes the pressure dependence of of ice VII to better than 20%. 3 With at least 13 polymorphs and the diversity of the hydrogen bond, the behavior of H 2 O under pressure is also a subject of considerable interest in the physics of condensed matter. 4,5 Data on the thermal conductivity of compressed H 2 O are available to a maximum pressure of only 2.4 GPa.
6For dielectrics such as H 2 O ices, oxides, and silicates, thermal conduction is largely controlled by phonon transport. The Leibfried-Schlömann (LS) formula 7,8 is among the most widely used schemes to describe the pressure dependence of . The LS equation is based on a detailed theoretical analysis of phonon transport 8 but it has not been experimentally tested over a range of pressures sufficient to change the Debye frequency, density, and elastic constants of a crystal by large factors. Furthermore, since the LS equation is based on the assumptions that acoustic phonons are the dominant carriers of heat and that the dominant scattering mechanism for acoustic phonons is three-phonon interactions between acoustic modes, its applicability to crystals with multiple atoms per cell has been questioned.9,10 Our previous work 11 showed that the pressure dependence of the cross-plane thermal conductivity of layered muscovite crystal could be adequately described by the LS equation when we assumed that the effective value of the Debye frequency varies as the square root of the cross-plane elastic constant C 33 . Muscovite is highly anisotropic, however, and the applicability of the LS equation in this case is not strictly valid.This study aims to measure the thermal conductivity of H 2 O over a pressure range that was not accessible previously. The data on ice VII, a cubic crystal with a relatively small bulk modulus, allow us to test the LS equation over a large compression ratio. Our method combines the time-domain thermoreflectance (TDTR) method 12 in a diamond anvil cell (DAC) with density functional theory (DFT) calculations of the vibrational density of states (DOS). At room temperature, cubic ice VII (space group P n3m) is stable between 2.1 and ≈60 GPa. 13,14 With its extrapolated zero pressure bulk modulus K 0 = 21.1 ± 0.5 GPa and the pressure derivative K 0 = 4.4 ± 0.1, 15 ice VII is compressed by more than 30% at 22 GPa.Symmetric DACs with 600-μm-culet diamonds and steel gasket were used to compress distilled H 2 O to 22 GPa. Pressure was determined from ruby fluorescence. 16 An 80-nm-thick Al film, coated on a 20-μm-thick sheet of muscovite mica [KAl 2 (Si 3 Al)O 10 (OH) 2 , grade V-1 from SPI Supplies], was lo...