## Abstract

In this publication, details of the calculation of heat conductivity and thermal diffusivity, compressibility, sound velocity and exciton‐polaron dispersion will be shown. The properties of excitons, coupling to phonons, producing thus polarons, but also blocking the phonons as running waves lead to an exciton condensation or exciton liquid. Surprisingly, this exciton liquid is contained in a macroscopic crystal, a solid, nevertheless, which becomes extremely hard due to the exciton liquid and finally exhibits a strange type of superfluid in a two fluid model, where the superfluid phase increases more and more below about 20 K until the whole exciton liquid becomes a superfluid at zero temperature. Never else a superfluid phase has been observed at such high temperatures.

### Keywords

- excition-polarons
- superfluidity
- exciton
- condensation
- heat conductivity
- compressibility
- sound velocity

## 1. Introduction

Excitons are electron‐hole pairs and as such known in many materials, even in Si. Generally, it is not easy to create so many excitons that they can interact with each other and finally can even condense in an exciton liquid. In standard experiments with laser pulses, one can excite in semiconductors electrons from a valence band into a conduction band and then, due to the electron‐hole attraction, the final state of the excited electron drops to somewhat below the bottom of the conduction band. An exciton is thus mobile, but it does not carry an electrical current due to its charge neutrality.

However, if one uses rare‐earth compounds, where the uppermost occupied state is a localized 4f state and the lowest conduction band is a 5d band and one makes an optically induced 4f‐5d transition, the hole in the 4f state will be localized and the electron just below the bottom of the conduction band will also be localized since it binds by Coulomb attraction to its hole, and the whole exciton will stay at the atom where the photoexcitation occurred and it will not be mobile. Then, this excitation energy will decay in typically 10^{−8} sec with the emission of a photon or phonons at the same atom where it is originated, resulting only in a localized excited atom in the lattice.

However, a p^{6}–5d exciton is always possible since the initial state is in a band. Thus, Mitani and Koda [1] found Mott‐Wannier excitons with thermo‐reflectance in Eu chalcogenides with about 4 eV (consult similar TmTe in Figure 1). The Tm mono‐chalcogenides exhibit a metal‐semiconductor transition inasmuch as Tm^{3+}S^{2−} + e is a trivalent metal with one free electron in the 5d conduction band (Figure 1). The occupied 4f^{12} level is about 6.5 eV below the Fermi energy E_{F}, and the empty 4f^{13} level is little above E_{F}. Experimental evidence comes mainly from X-ray photoemission spectroscopy (XPS) and Bremsstrahlen isochromat spectroscopy (BIS) [2]. TmTe on the other hand is a divalent semiconductor Tm^{2+}Te^{2−} with an occupied 4f^{13} level 0.3 eV below the bottom of an empty 5d band (Figure 1). The driving force behind this different character is the crystal field splitting of the 5d band, which depends on the lattice constant due to the different anion radii.

TmSe, on the other hand, with an intermediate anion radius between sulfide and telluride has such a crystal field splitting of the 5d band that the bottom of this band overlaps with the 4f^{13} level. This f‐d hybridization on the one hand leads to some d‐character of the f‐state and as a consequence to a narrow f‐band and on the other hand to some f‐character of the bottom of the d‐band. It has been conventional to describe the new hybridized f‐state as 4f^{13}‐4f^{12}5d, consisting of a quantum mechanically mixed state [3]. This new phenomenon is called intermediate valence, since the valence of rare‐earth ions is defined by the occupation of the f‐state, and thus. TmSe has a valence, between 3^{+} and 2^{+}, in fact 2.85^{+}. This can only be achieved if the 4f‐state is a narrow band, which is partially filled with electrons [3]. But also in the 5d band there are some free electrons, which yield in the visible a coupled plasma resonance of these electrons and are responsible for the copper‐like color of the crystals.

## 2. Material tailoring

One can now make mixtures between the semiconducting TmTe and the metallic intermediate valence TmSe and thus tune the energy gap ΔE between 300 meV and zero (metal) [3, 4]. Experimentally semiconducting TmSe_{1−x}Te_{x} has been created with x = 0.40, 0.55, 0.68, corresponding to ΔE of 40 meV, 110–120 meV, 170 meV. For these compositions the f‐state is so close to the 5d band that some hybridization occurs between the tails of the wave functions. We have the unexpected situation of intermediate valence semiconductors. This in turn means that the originally localized 4f^{13}‐state acquires now some bandwidth in the order of tens of meV.

Concerning now the existence of 4f‐5d excitons, we have created a situation where the hole state in the 4f is mobile and the electron is in a 5d state below the bottom of the 5d band. Now we have the possibility of 4f‐5d excitons. Of course the hole in a narrow 4f‐band has a large effective mass, so that the 4f‐5d exciton is a heavy (m_{h} ≈ 50 m_{e}) nearly immobile particle. These excitons without application of external pressure have a low concentration at low temperatures because thermal excitations into the excitonic state are rare.

From Bohr’s formula for the hydrogen radius a_{H} = 0.53 ε Å and from Figure 2 with the reflectivity for ω → 0 = 50%, we obtain . So the static dielectric constant is ε = 34, and the radius of the orbit of the exciton is about 18 Å. This would be a Mott–Wannier exciton. Its binding energy from the optical result is theoretically E

_{B}= 13.6/ε

^{2}= 15 meV, about the same as from the electrical measurement [4], but experimentally from Figure 2 E

_{B}≈ 60 meV. The most complete measurements have been performed on TmSe

_{0.45}Te

_{0.55}with an energy gap ΔE of 110–120 meV, confirmed by far infrared reflectivity (see inset of Figure 2). Similar absorptive peaks as for TmSe

_{0.45}Te

_{0.55}are absolutely missing in other divalent rare‐earth chalcogenides with only localized 4f‐states (Figure 3).

In the fcc rocksalt structure, the 4f^{13}‐band has a maximum at the Γ point of the Brillouin zone and a minimum at the X point. The 5d band dispersion has its minimum at the X point. An optical transition between the maximum of 4f and the minimum of 5d would be an indirect transition and requires maximal the assistance of a Γ–X phonon for k conservation (black curve).

## 3. Creation of excitons

Under hydrostatic pressure, the bottom of the 5d band at X with its exciton level will be lowered with respect to its center of gravity (5dt_{2g}‐5de_{g}) and shown for 5 kbar the exciton level is exactly at the energy of the 4f‐state at Γ (red curve). Now the highest energy electrons in the 4f^{13}‐band can spill without energy loss into the excitonic state at X leaving behind a positive hole. This transition needs the emission or absorption of Γ–X phonons which couple to the excitons. So in fact we are dealing with an exciton‐polaron. With higher pressure, the bottom of the 5d band at X will approach the energy of the 4f^{13}‐state at Γ and the 4f electrons will enter directly the 5d band and perform a first‐order semiconductor–metal transition.

In Figure 4, these transitions can be observed directly with resistivity in the isotherms versus pressure for TmSe_{0.45}Te_{0.55}. We look at first at room temperature (300 K) and find a classical pressure dependence of a resistivity, namely the resistivity of a semiconductor decreases with increasing pressure, because the energy gap ΔE decreases with pressure and bands widen and finally the metallic state is achieved (above 11 kbar). Starting with about 5 kbar and best observed at 5 K, the resistivity now increases by about three orders of magnitude with pressure. This is exactly the pressure where excitons become stable states and electrons from the f‐band, which have been thermally excited into the 5d conduction band, drop into the excitonic state and are no longer available for electric conduction. We have created an excitonic insulator, a term coined by Sir Nevil Mott [7]. With further pressure increase, the resistivity drops again, until now the 4f electrons can enter the 5d band directly, which leads to a first‐order semiconductor–metal transition. Here we want to make a remark of another possibility of the semiconductor-metal transition, namely a Mott transition to an electron-hole plasma or an electron-hole liquid. The experimentally derived exciton concentration is 3.9 x 10^{21}cm^{−3} (see below). This is in fact too high (because of screening effects) for an electron-hole liquid as has been shown by Monnier et al. [8]. There it is calculated that the electron-hole liquid must be less than 10^{20}cm^{−3} excitons. In fact the rare earth nitrides may serve as experimental examples [9].

## 4. Exciton condensation

Since for the exciton creation no energy is needed, their number is enormous. But not all 4f electrons can form excitons, because as electric dipoles and according to the Pauli principle [10, 11] they repel each other. This goes so far that the formation of this incredible high concentration of excitons forces the whole crystal lattice to expand against the applied pressure. We show this in Figure 5 where we measure the lattice constant (Figure 5a) (with strain gauges) and the expansion coefficient (Figure 5b) of the crystal in an isobar at 11.9 kbar. We observe that at about 230 K the lattice expands by 1.6% isostructurally, an enormous amount. The expansion coefficient becomes negative, of course. We even think that the expansion is of first order (dashed–dot line), but the point‐by‐point measurement cannot reproduce this exactly, because we go from the semimetallic state to the excitonic state.

We can estimate the maximal number of excitons with the help of Figure 5a, and we observe that the lattice expansion occurs spontaneously when entering the excitonic phase. There must be an energy balance between the lattice energy causing the expansion and the electronic energy of the excitons. The energy balance can be described by the first equation in Figure 6. We take the lattice constant change from Figure 5a to go from 5.93 to 6.03 Å and compute Δl/l and ΔV/V. We choose a pressure of 8.5 kbar and an E_{B} of 70 meV and compute the number of excitons n_{ex} = 3.9 × 10^{21} cm^{−3} (red field). We also can compute the number of Tm ions in the crystal in the fcc structure, and it is n_{Tm} = 1.8 × 10^{22} cm^{−3} (yellow field). In other words, the exciton concentration is about 22% of the atomic density, an enormous amount of excitons. With the exciton orbit of 18 Å, it is quite clear that we have an exciton band or an exciton condensation. Since the exciton couples to a phonon, the condensation is a Bose condensation, not a Bose–Einstein condensation. We can also estimate the Bose condensation temperature shown in Figure 6, where the general accepted formula yields T_{B} = 130 K, the right order of magnitude. The holes of the exciton are in a narrow 4f‐band, and with a pressure change of 5 to 8 kbar (Figure 4), one scans the width of the 4f‐band [12]. The closing rate of the semiconductor rate has been measured to be dΔE/dp = −11 meV/kbar [4], so 3 kbar · 11 meV = 33 meV for the width of the narrow 4f‐band. From this in turn, we use the general estimate that a band width of 1.5 eV yields an effective mass of m_{e} and derive that a band width of 30 meV corresponds to an effective hole mass m_{h} ≈ 50 m_{e}. The excitons are thus heavy bosons.

Here we want to make some remarks about this exciton condensation. Nobody in the world (to the best of our knowledge) has a comparable concentration of excitons which exist as long as we can sustain the pressure and as the liquid Helium lasts, this means for days. We can make all kinds of experiments in this condition, such as electrical conductivity, Hall effect, compressibility, heat conductivity, superfluidity, ultrasound velocity, phonon dispersion and specific heat. Nobody else has these possibilities. But the experiments are very demanding at low temperatures with simultaneous pressure and doing the measurements.

## 5. Phase diagram of semiconductor, excitonic insulator and semimetal

We plot in Figure 7 the coexistence ranges of the intermediate valence semiconductor, the excitonic insulator and the intermediate valence semimetal. We see that the highest temperature for which the excitonic insulator exists is about 260 K and the pressure range is between 7 and 13–14 kbar (pressures applied at room temperature). Experimentally one can only measure isobars in a clamped pressure cell. However, the isobars in Figure 7 are no straight lines, because the pressure applied at room temperature relaxes somewhat at low temperatures. In the inset of Figure 7, we see the Hall effect, which measures the free electron concentration in the 5d band. In the semimetallic state (curve M) at 13 kbar, the electron concentration is about 3 × 10^{21} cm^{−3}. For the excitonic insulator at 8 kbar, the free electron concentration is about 10^{18} cm^{−3} because now the free electrons condense into the excitons and do not contribute anymore to the Hall effect. In fact we observe that the carrier concentration reduces by about three orders of magnitude, the same as has been observed in Figure 4 for the electrical resistivity. The change in resistivity is thus mainly an effect in the carrier concentration though the mobility changes also somewhat [13]. The concentration of the excitons is then 3 × 10^{21} cm^{−3}–10^{18} cm^{−3} = 3 × 10^{21} cm^{−3}, about the same as has been obtained in Figure 4.

We can consider in an analogy a pot with soup. The pot is the hard surrounding of the crystal and inside is a soup of liquid excitons.

In Figure 8, we show in the upper part a proposal from Kohn [14] from 1968 with the three‐phase semiconductor, excitonic insulator and semimetal plotted against the energy gap ΔE with increased pressure going to the left. When ΔE = E_{B}, the excitonic instability starts. In the lower part of the figure, we show the E‐k diagram again for the three phases. It is surprising and satisfying that the foresight of Walter Kohn has practically reached reality by comparing the inset of Figure 8 with the real phases of TmSe_{0.45}Te_{0.55} as shown in Figure 7.

## 6. Isotherm and compressibility

We may ask what the direct evidence for the condensed excitonic state is. Typical for any liquid is its incompressibility. We can, for instance, at 1.5 K, apply an increasing pressure to TmSe_{0.45}Te_{0.55}, and this is shown in Figure 9 [12]. At first, we cool at zero pressure from 300 to 1.5 K and volume and lattice constant decrease. Then, we increase pressure and measure the lattice constant with elastic neutrons through the pressure cell. Of course, lattice constant and volume decrease further, corresponding to a Birch–Mournaghan equation (red curve). This is a very time‐consuming experiment, because for each pressure change the pressure cell had to be heated to room temperature to change to a higher pressure and then cooled down again and adjust the sample in the neutron beam and wait for beam time. Therefore, this experiment has only four points, but at the relevant pressures. As can be seen in Figure 9 when entering the excitonic state, the lattice constant remains unchanged with increasing pressure, which means a compressibility of zero, as shown in Figure 10.

Taking experimental uncertainties into account, we have at least a compressibility just as for diamond. Thus, we can take this experiment as evidence of an excitonic liquid.

## 7. Heat conductivity and superfluidity in the excitonic liquid

We now want to discuss the possibility of superfluidity in the excitonic liquid. Here we resort at first to theory [11, 15]. There is a similarity between pairs of particles: two electrons can condense and produce superconductivity, and an electron–hole pair (exciton) can upon condensation result in superfluidity. A positron pair should also result in superconductivity, but no such experiment is known.

In any case, our exciton condensation may result in superfluidity. What would be the experiment to prove this? In our opinion, this is heat conductivity [16], because it would diverge when the material becomes superfluid, just as in ^{4}He [17].

The experimental arrangement to measure heat conductivity and thermal diffusivity in a pressure cell is described in detail in Ref. [16], but the essence is isobars between 4 and 300 K at various pressures. We show the results of measurements of the heat conductivity λ with isobars at four different pressures, one in the semiconducting range (compare Figure 7) with 7 kbar, one in the semimetallic range at 15 kbar, both outside the excitonic region, and at two pressures 13 and 14 kbar within the excitonic range. Temperature has been measured automatically for each degree. In Figure 11, we collect a few relevant formulae for the heat conductivity λ and the thermal diffusivity a. We see that the heat conductivity depends on the specific heat c_{v} and l_{ph} in direction x, the mean free path for phonon scattering. In short, l_{ph} will increase with decreasing temperature because the density of phonons decreases and we have Umklapp processes involving three phonons. But the specific heat c_{v} definitely will go toward zero for zero temperatures; thus, the heat conductivity outside the excitonic region will display a maximum near 50 K, as well for the semiconducting range (7 kbar) as for the metallic range (15 kbar); this is displayed in Figure 12, and this behavior is quite normal. The difference of the heat conductivity near 300 K for both cases is due to the electronic part of the heat conductivity in the metallic state, and it corresponds roughly to the Wiedemann–Franz relation. This gives confidence to the measurements. We continue with the heat conductivity in the excitonic region at 13 and 14 kbar. We observe an unexpected downward jump in a first‐order transition when entering the excitonic phase. Consulting Figure 7, it is obvious that at different pressures one enters the excitonic phase at different temperatures. At these temperatures and pressures, one enters the insulating excitonic phase mainly from the semimetallic phase, thus with a metal–insulator transition. The downward jumps in the heat conductivity λ reflect the loss of the electronic part of the heat conductivity. The fascinating aspect of the heat conductivity in the excitonic region is the sharp increase of λ below about 20 K, quite in contrast with the λ outside the excitonic region. Since λ follows mainly the specific heat c_{v} and the phonon mean free path l_{ph}, (Figure 11) and c_{v} nevertheless must go to zero for T→0, it is the phonon mean free path which goes faster to infinite than c_{v} toward zero. Finally, it means that the phonon mean free path becomes infinite. When one makes a heat pulse at one end of the crystal, the excited phonon transports its energy without scattering on other phonons to the other side of the crystal, meaning an infinite heat conductivity. This is, however, only possible if the concentration of phonons as running waves is substantially reduced, because most of them couple to the heavy excitons as exciton‐polarons, as we have seen before and thus more or less correspond to local modes.

Unfortunately, the measurements were limited to 4.2 K, because at the time of the measurements one did not realize the implications. In any case 20 K, the onset of the sharp increase of λ with decreasing temperature can be considered as the onset of superfluidity, which, however, is different from the one of ^{4}He, inasmuch as there the onset of superfluidity is a first‐order transition [17]. For our exciton case, we propose a superthermal current in the two‐fluid model, where the superfluid part increases gradually toward zero temperature [17]. A λ‐anomaly in the specific heat as in the first‐order Bose–Einstein transition in ^{4}He is here not to be expected and also not found [17].

The proposed evidence of superfluidity within the condensed excitonic state necessitates an additional excitation spectrum of other quasiparticles, namely rotons or vortices [17]. λ_{tot} is the sum of individual contributions (Figure 11), and below about 20 K λ_{tot} = λ_{ph} + λ_{ex}. λ_{ph} is the heat conductivity due to uncoupled phonons, which is proportional to T^{3} and can be neglected compared to λ_{ex} at low temperatures. Thus, we obtain for λ_{ex} an Arrhenius law for the increase of the heat conductivity toward zero temperature λ_{ex}_{B}T. This is shown in Figure 13 at 13 kbar. The activation energy or the gap Δ is 1 meV or about 10 K. The application of heat in the heat conductivity experiment can excite quasiparticles, e.g., rotons with gap energy of about 5 K, which is the right order of magnitude. In superfluid ^{4}He, the roton gap is 8.65 K [18].

## 8. Thermal diffusivity

In Figure 14, we display the thermal diffusivity a for the same four pressures as in Figure 12. The thermal diffusivity a = v_{x}l_{x} and thus follows mainly the phonon mean free path l_{x}, with v_{x} being about constant outside the excitonic region, consulting Figure 11. In fact outside the excitonic region with 7 and 15 kbar, it does exactly this, as can be seen in the theoretical curve for l_{x} in Figure 15. For the lowest temperatures, the phonon mean free path becomes the dimension of the crystal and is thus constant. But in the excitonic region, again below about 20 K, the thermal diffusivity increases dramatically. Why then in the excitonic region the dimensions of the crystal do not seem to be important now? Just as in superfluid Helium heat can be transferred not only via phonon‐phonon scattering in a diffuse manner, but ballistically via a highly directional quantum mechanical wave, the second sound. Also above 20 K, there are anomalies, but they can be explained with the velocity of sound v_{x}, which increases now strongly in the excitonic region (see below).

In principle, the two measurements of heat conductivity and thermal diffusivity permit the calculation of the specific heat c_{v} = λ/ρa, with ρ the density (see Figure 11) [16] and we did this in Figure 7 of Ref. [16]. But we never felt very happy with this curve because we divided two point‐by‐point measurements. But the specific heat in the excitonic range is definitely below the one of the specific heat outside this range. But it is also very complex since the density ρ diminishes when entering the excitonic phase, because the crystal expands (see Figure 5). It took us several years before we could make a direct measurement of the specific heat under pressure and below 300 K [6, 19].

## 9. The specific heat

The specific heat c_{v} has been measured for TmSe_{0.45}Te_{0.55} [6, 19] along isobars with 0 kbar, and corresponding to the curves K, N, M in Figure 7, which is shown in Figure 16 [6, 19].

The molar specific heat c_{m} in J/mole K/f.u. has been matched at 300 K to the Dulong‐Petit value of 52 J/mole K/f.u. The specific heat at ambient pressure represents a normal Debye curve (black curve). This curve has been measured by our colleagues at the university of Geneva and ETH Zürich for T > 1 K and T > 0.3 K [20]. A Schottky anomaly due to crystal field splitting of the Tm ions and an exchange splitting due to magnetic order at 0.23 K has been subtracted from the measured curve, and the pure phonon contribution could be plotted as c_{v}/T versus T^{2}; thus, a Debye temperature ϴ of 117 K could be obtained [16]. Curve N is in the semimetallic high‐pressure phase outside the excitonic region, and we find again a normal Debye curve (blue curve), but with a lower Debye temperature than at ambient pressure. This is at first sight surprising since at high pressure a solid becomes harder with a higher Debye temperature, but it has also been observed in reference [16]. A simple explanation can be that with high pressure we change somewhat the degree of valence mixing in the intermediate valence semimetal. Curve M (red curve) starts with about 13 kbar at 300 K in the metallic region, but enters the excitonic region at about 150 K. This occurs with a first‐order transition as we see in Figure 5a. The specific heat should reveal a delta‐function at the phase transition, but experimentally the spike reduces to a Gaussian shape [21]. All curves entering the excitonic region in Figure 7 from the semimetallic region in a first‐order transition (red squares and downward triangles in Figure 7) exhibit the Gaussian anomaly. We now discuss curve K (mauve), which enters the excitonic region in a second‐order transition without a spike in the specific heat. We observe that the specific heat is no longer a Debye curve, but below about 250 K (arrow in Figure 17) one finds a quasi‐linear drop of the specific heat until below about 30 K the specific heat joins the other curves. These measurements reveal a fundamental difference of the thermodynamic phases A and B in Figure 7.

In fact, such a specific heat like curve K with a non‐Debye like curvature has never been seen before. Since the specific heat over a higher and larger temperature interval is entirely given by the phonons (in the absence of magnetic order and special effects like Schottky anomalies), we must conclude a strong renormalization of the phonon spectrum.

As we have stated already several times above the excitons in this indirect semiconductor couple strongly to phonons in a triple particle entity of hole‐electron and phonon as an exciton polaron. But when the phonons couple to the heavy excitons with effective masses of the holes around m_{h} = 50 m_{e} they become more or less localized like a local mode and do no longer contribute significantly to the specific heat. So an essential part of the Debye spectrum of the specific heat is missing.

So in Figure 17 we have made a model calculation of the specific heat with the assumption that the optical phonons are bound to the excitons below about 250 K (arrow in Figure 17). The acoustic branches of the phonons are modelled with a Debye and the optical phonons with an Einstein ansatz, respectively [6]. The seeming disappearance of phonons, i.e., the binding of the optical phonons on the excitons is represented with a linear decrease of the density of states below 250 K (inset to Figure 17). The model calculation in Figure 17 represents well the intriguing behavior of curve K in Figure 16. A strong coupling regime for the phonons to excitons prevails, and thus, while cooling, more and more wave‐like phonons become locked onto the excitons, giving no more contribution to the specific heat. This renormalization of the phonon spectrum and the resulting effect on the specific heat has never been seen before, and it is due to the extreme large concentration of exciton–polarons. Regarding now the Debye temperature ϴ of curve K in comparison with curve N, we observe a further reduction of the Debye temperature, i.e., a minimum in the Debye temperature versus pressure in the excitonic region (red curve in Figure 17).

In a quantitative formula, we can express

(Debye) + (Einstein) with ρ_{opt} a temperature‐dependent density of optical phonons (see inset Figure 17). In fact, the model calculation in Figure 17 represents quite well the measured specific heat of curve K in Figure 16. Thus, the acoustic phonons alone exhibit a Dulong‐Petit value c/R of 3 cal/degree and the optical phonons have a temperature‐dependent density, their decrease with temperature representing the increase of excitons–polarons with decreasing temperature. The free optical phonons get lost for, e.g., the thermal conductivity. However, below about 20 K the excitons–polarons take over in the heat conductivity or the thermal diffusivity and with a diverging increase in these entities finally lead to superfluidity.

Curve M in Figure 16 is on the decreasing branch of exciton concentration (see Figure 7) where more and more free electrons in the 5d band are screening the Coulomb interaction between electron and hole. In dissolving the excitons in region B of Figure 7, the electrons from the excitons enhance the 5d electrons, further which leads to a cumulative process and a collective breakdown of the rest excitons in a first‐order transition. Curve M enters the excitonic region at a temperature of about 150 K where the exciton concentrations are already about two orders of magnitude lower than at the maximum.

## 10. Sound velocity

An ultrasound transducer has been glued to one end of the crystal, and with a multiple echo from the other end of the crystal over the known length of the crystal, the sound velocity could be obtained. This is shown in Figure 18 for various pressures. At zero pressure and at 7 and 18 kbar, the sound velocity is about 4000 m/s and there is not much change with pressure. But best seen at 12 kbar, when entering the excitonic phase at 180 K, the sound velocity is enhanced by nearly a factor 2 (see Figure 7). With 10 kbar, we are entering the excitonic phase at 240 K, again with a jump of nearly a factor 2 but near 90 K the pressure loss in the cell was just the size for a reentrant transition to the non‐excitonic phase. This was a unique phenomenon, but supporting the experimental measurements.

As mentioned above in the chapter about the thermal diffusivity a, the upwards jumps in the excitonic region are indeed caused by the jumps in the sound velocity. But not only this, the increase in sound velocity contributes directly to the thermal conductivity above 20 K and is responsible for the bumps in the thermal diffusivity.

The sound velocity is related to the bulk modulus B and its inverse the compressibility. For a cubic material B depends on the elastic moduli c_{ij} as B = 1/3 (c_{11} + 2c_{12}) the elastic moduli instead, depend on the sound velocity as c_{11} = ρv_{L}_{2}_{[100]} and c_{12} = ρ(v_{L}^{2}_{[100]} – 2v_{T2}^{2}_{[110]}). Assuming that in general v_{L} is about 3 times v_{T2} we get the simplified relation B ≈ ρv_{L}^{2}. Thus in the excitonic phase we find a 2 times larger v_{L} and thus a 4 times larger B or a 4 times smaller compressibility. The material gets indeed appreciable harder in the excitonic state.

In Figure 9, we have shown that between 5 and 8 kbar the lattice constant remained practically constant during exciton condensation, meaning that the compressibility is close to zero. Putting a maximal error bar through the points of measurement a bulk modulus B = 20 GPa outside the excitonic region and a bulk modulus B = 70 GPa in the excitonic region could be obtained. From the sound velocity measurement in Figure 18, we calculate a bulk modulus B = 24 GPa outside the excitonic region and one of 100 GPa in the excitonic region. So both types of measurements agree reasonable well and confirm the fact that during exciton condensation the material becomes extremely hard.

We offer two explanations for this phenomenon: the electron from the exciton enters a 5d‐like orbit, which is much larger than the original 4f orbit it came from, and this in spite of the increasing pressure. Or the excitons, being electric dipoles, repel each other at short distances and large concentrations, creating a counter‐pressure to the applied pressure.

The dominant feature in Figure 18 is the sharp increase by a factor two of the sound velocity when entering the excitonic phase. But also at 300 K in an isotherm taken from Figure 18 with the relevant pressures the sound velocity is changing. Now at 300 K this is shown in Figure 19, but it has nothing to do with excitons and the change is much smaller than the one due to excitons. Here we find a minimum of the sound velocity with increasing pressure, inverse to what we have discussed in the excitonic region. The relation of bulk modulus B with sound velocity is B ≈ ρv_{L}^{2}, and a minimum of sound velocity implies a minimum in the bulk modulus, which is inverse to the compressibility. Thus, the minimum in the sound velocity means a maximum in the compressibility. This can be compared with the compressibility for 300 K in Figure 10, and we obtain a similar curve. At 300 K, the softening of the bulk modulus or a maximum in the compressibility is here due to a change of the degree of valence mixing with pressure. So these completely different experiments (also by different authors [4]) support each other and give again confidence into the experiments.

## 11. Dispersion of exciton‐polarons

We mentioned above regarding Figures 16 and 17 that the Debye temperature in the excitonic region is less than the Debye temperature ϴ = 117 K at 300 K but in Figure 18 we observe that the sound velocity is enhanced in the excitonic region. This seems to be a contradiction since in the Debye model the sound velocity is the slope of a linear phonon dispersion curve where the maximum frequency ω_{max} determines the Debye temperature ϴ. A lower Debye temperature has thus a lower sound velocity and a lower bulk modulus. In order to explain a lower ω_{max} together with a higher sound velocity, we have to leave the simple Debye model for bare longitudinal acoustic (LA) phonons ω _{LA (Γ‐X)} ≈ 14 meV [12], in Figure 20b an exciton with 4f character E_{B} ‐ ^{2}k^{2}/2M, where M is the sum of electron and hole mass of the exciton M = m_{h} + m_{e} ≈ m_{h} with m_{h} ≈ 50 m_{e}.

We can see in Figure 20 that now the dispersion of the exciton‐polaron has indeed simultaneously a steeper slope (larger sound velocity) than the LA phonon and a lower ω_{max} than the simple phonon (smaller Debye temperature).

## 12. Theoretical Models

Since the binding energy E_{B} of the exciton-polaron is with 60-70 meV relatively large also a Frenkel type of exciton-polaron is conceivable. Thus in a theoretical paper [24] it is proposed that the exciton condensation occurs in an extended Falikov-Kimball model [25] where, instead of the original model with localized 4f states a narrow hybridized 4f band is used, which is more realistic in this case. Extensively discussed has also been an effective mass model [26] with large differences between electron and hole mass, just as we proposed above. In a further paper [27] it has been shown, that weakly overlapping Frenkel type excitons can condense. Especially the coupling of excitons with phonons has been discussed in Ref. [28] and the formation of exciton-polarons. Finally in [29] exciton densities and superconductivity (sic) are discussed where for low exciton densities a Bose-Einstein condensate is proposed and for high density a Bardeen-Cooper-Schrieffer condensate should prevail, especially in coupled bilayers. This certainly is not the case in our experiments.

## 13. Conclusion

In this review paper, we treat a special rare‐earth material, TmSe_{0.45}Te_{0.55} which has been tailored so that with moderate pressures (up to 20 kbar) and low temperatures (down to 4 K) an enormous amount of excitons (10^{21} cm^{−3}), about 22% of the atomic density, can be statically obtained. This high concentration of excitons with Bohr orbits of about 18 Å leads to a condensation of excitons, which forms a liquid inside a crystalline surrounding. The existence range of condensed excitons is below 250 K and between 5 and 14 kbar. The condensation is accompanied with a phenomenon of incompressibility and as such with a compressibility near zero. In this condition, the heat conductivity and the thermal diffusivity have been measured in order to investigate a possible superfluidity which has been proposed by Keldysh and Kopaev [11] and Kozlov and Maksimov [15]. Outside the excitonic region, both entities behave quite normal, whereas in the excitonic region the heat conductivity diverges to ever‐increasing values. This can be explained below 20 K within a two‐fluid model, where the superfluid part always increases until at temperature zero the complete condensed excitons become superfluid. Also the thermal diffusivity expands in the excitonic region above the phonon mean free path corresponding to the size of the crystal. This can be explained with the quantum‐mechanical second sound, which is a ballistic transport of heat. The Debye temperature exhibits a minimum in the excitonic region where nevertheless the sound velocity is increasing. These two incompatible measurements can be explained with a strong phonon renormalization in the excitonic region, and the Debye phonon dispersion of LA phonons changes into the dispersion of an exciton–polaron, because every exciton binds to a phonon. This in turn means that the number of free phonons is strongly reduced in the excitonic region so that the specific heat becomes extremely anomalous, far away from a Debye specific heat. In general, it can be said that the anomalous physical properties of condensed excitons are unprecedented.

TmSe_{0.45}Te_{o.55} is not the only material where these phenomena can be observed. YbO and YbS are similar materials though one will need much larger pressures to close their gaps of about 1 eV [30]. But also Sm_{0.75}La_{0.25}S [31] is a possible candidate for which much lower pressures are needed. So with a good feeling for materials new and exciting effects can be found.

## Acknowledgments

The author wishes to thank Prof. Benno Bucher of HSR, Switzerland, for rechecking the experimental results and fruitful discussions.

## References

- 1.
Mitani T, Koda T. Proceedings of the 12th Int. Conf. Phys. Semiconductors, Stuttgart, Teubner publ. 1974; p. 889 - 2.
Wachter P, Kamba S, Grioni M. Physica B. 1998; 252 :178 - 3.
Wachter P: In: Gschneidner jr. KA, Eyring L, Lander GH, Chopin GR. Editors. Handbook on the Physics and Chemistry of Rare Earths, 19 Lanthanides/Actinides: Physics II. 1991. p. 177 - 4.
Boppart H, Wachter P. Proceedings of Material Research Society. 1984; 22 :341 - 5.
Neuenschwander J, Wachter P. Physical Review B. 1990; 41 :12693 - 6.
Wachter P, Bucher B. Physica B. 2013; 408 :51 - 7.
Mott N. Philosophical Magazine. 1961; 6 :287 - 8.
Monnier R, Rhyner J, Rice TM, Koelling DD. Physical Review B. 1985; 31 :5554 - 9.
Wachter P. Advances in Material Physics and Chemistry. 2015; 5 :96 - 10.
Halperin BI, Rice TM. Reviews in Modern Physics. 1968; 40 :755 - 11.
Keldysh LV, Kopaev AN. Soviet Physics Solid State. 1965; 6 :2219 - 12.
Wachter P. Solid State Communication. 2001; 118 :645 - 13.
Bucher B, Steiner P, Wachter P. Physical Review Letters. 1991; 67 :2717 - 14.
Kohn W. In: Many Body Physics, edited by de Witt C, Balian R. New York: Gordon & Breach; 1968 - 15.
Kozlov AN, Maksimov LA. JETP. 1965; 21 :790 - 16.
Wachter P, Bucher B, Malar J. Physical Review B. 2004; 69 :094502 - 17.
Tilley DR, Tilley J. In: Superfluidity and Superconductivity. edited by Brewer DF. Bristol and New York: Adam Hilger; 1990. P. 35 - 18.
Henshaw DG, Woods ADW. Physical Review. 1961; 121 :1266 - 19.
Bucher B, Park T, Thompson JD, Wachter P. arXiv:0802.3354v2 - 20.
The authors are most grateful to Prof. A. Junod (University of Geneva) and CH. Bergeman ETH Zürich for performing specific heat measurements at ambient pressure - 21.
Schilling A, Fischer RA, Phillips NE, Welp U, Kwok WK, Crabtree GW. Physical Review Letters. 1997; 78 :4833 - 22.
Hodgson JN. In: Optical Absorption and Dispersion in Solids. London: Butler and Tanner; 1970 - 23.
Wachter P, Bucher B. Physica Status Solidi ©. 2006; 3 :18 - 24.
Ihle D, Pfafferot M, Burovski E, Bronold FX, Fehske H. Physical Review B. 2008; 78 :193103 - 25.
Falikov LM, Kimball JC. Physical Review Letters. 1969; 22 :997 - 26.
Bronold FX, Fehske F. Physical Review B. 2006; 74 :165107 - 27.
Zenker B, Ihle D, Bronold FX, Fehske H. Physical Review B. 2012; 85 :121102 - 28.
Zenker B, Fehske H, Beck H. Physical Review B. 2014; 90 :195118 - 29.
Soller H. Journal of Applied Mathematical and Physics. 2015:1218 - 30.
Wachter P. Journal of Alloys and Compounds. 1995; 225 :133 - 31.
Wachter P, Jung A, Pfuner F. Physical Letters A. 2006; 359 :528