## Optical investigation of thermoelectric topological crystalline insulator Pb 0 77 . Sn 0 23 . Se

1, ∗ 1 1

Anjan A. Reijnders, Jason Hamilton, Vivian Britto, Jean-Blaise Brubach, 2 Pascale Roy, 2 Quinn D. Gibson, 3 R.J. Cava, 3 and K. S. Burch 4, †

4

1 Department of Physics &amp; Institute for Optical Sciences, University of Toronto, 60 St. George Street, Toronto, ON M5S 1A7, Canada. 2 Synchrotron SOLEIL, L'Orme des Merisiers Saint-Aubin, 91192 Gif-sur-Yvette, France. 3 Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA Department of Physics, Boston College, 140 Commonwealth Ave, Chestnut Hill, Massachusetts 02467, USA

(Dated: August 4, 2014)

Pb 0 77 . Sn 0 23 . Se is a novel alloy of two promising thermoelectric materials PbSe and SnSe that exhibits a temperature dependent band inversion below 300 K. Recent work has shown that this band inversion also coincides with a trivial to nontrivial topological phase transition. To understand how the properties critical to thermoelectric efficiency are affected by the band inversion, we measured the broadband optical response of Pb 0 77 . Sn 0 23 . Se as a function of temperature. We find clear optical evidence of the band inversion at 160 ± 15 K, and use the extended Drude model to accurately determine a T 3 / 2 dependence of the bulk carrier lifetime, associated with electron-acoustic phonon scattering. Due to the high bulk carrier doping level, no discriminating signatures of the topological surface states are found, although their presence cannot be excluded from our data.

## I. INTRODUCTION

Thermoelectrics are a subject of great interest to industry for their many applications as macro and micro-scale Peltier coolers, their waste-heat energy recycling potential (heat-to-electricity conversion) in industrial and consumer products, and their application as highly sensitive compact and robust sensors. 1-5 An academic interest has also persisted as the optimization of thermoelectric efficiency offers an intriguing challenge in which ostensibly mutually exclusive physical properties must be combined in one material for the greatest performance. This is illustrated by the figure of merit ZT = TS σ 2 κ , where S σ , , and κ are the Seebeck coefficient, the DC conductivity, and the thermal conductivity, respectively. For high thermoelectric efficiency, S and σ must therefore be simultaneously enhanced, while κ is deliberately suppressed. This is a nontrivial task, considering the intimate link between the parameters that define ZT . For example, both σ and κ depend on the carrier lifetime τ , while both S and σ are a function of the band dispersion. Hence, the improvement of each individual quantity could result in an overall decrease of ZT , thus making experimental techniques that can disentangle the common limitations of each quantity highly desirable.

optical investigation of strongly correlated systems, 16-20 the extended Drude model offers some benefits previously unexplored in the context of thermoelectrics. Using this approach grants access not only to the temperature dependence of the carrier lifetime, but also to its frequency dependence, shedding light on the various scattering mechanisms that govern σ and κ . Moreover, with both high temperature and energy resolution, 21 broadband (far infrared through visible) optical spectroscopy provides direct access to the free carrier plasma frequency and band mass, both affecting S and σ , as well as spectral weight transfers and band transitions, useful for studying band gaps and band inversions.

Currently, the primary methods of studying thermoelectrics are by DC transport, optics, thermopower experiments and neutron scattering, 6-15 where the standard Drude model is often employed to extract temperature dependent information about carrier lifetimes. A major shortcoming of the Drude model, however, is the assumption that the carrier lifetime is frequency independent. When analyzing optical data, this often results in poor quality Drude fits, and carrier scattering rates with large uncertainties. A solution to this shortcoming is the application of the extended Drude analysis to a sample's optical response. 16 While commonplace in the

To illustrate the versatility of this approach we've performed our measurements on Pb 0 77 . Sn 0 23 . Se, a compound combining the thermoelectric characteristics of its thoroughly studied binary constituents PbSe and SnSe 9-12,22-24 However, unlike its constituents, previous works on Pb 0 77 . Sn 0 23 . Se crystals have shown evidence of a temperature dependent band inversion, although neither the exact transition temperature, nor the temperature dependence of the gap has been identified conclusively. This is largely a result of limited energy and temperature resolution of the experimental probes that were used, with various reports ranging between 80 K and 200 K. 6,8,9,25-27 Hence, broadband optical spectroscopy provides a unique opportunity to study, in a very controlled manner, how carrier dynamics change throughout the inversion. Moreover, many traditional thermoelectric probes, such as transport and Seebeck measurements, rely on band structure models to extract the scattering rate. Such band structure assumptions are severely complicated by the band inversion, further motivating the investigation of Pb 0 77 . Sn 0 23 . Se by optical spectroscopy as features of the band inversion and the carrier lifetime can be measured independently in a single experiment. This discriminating quality becomes even more evident by considering previous

transport studies of the Pb 1 -x Sn Se x system. While Pb 1 -x Sn Se x alloys have been studied for 50 years, 28 only specially prepared low-carrier extrinsic p-type Pb 0 77 . Sn 0 23 . Se was found to show abrupt changes in the DC transport properties across the band inversion. 8 Interestingly, as-grown n-type crystals showed no clear signatures of the band inversion at all. A recent paper by Liang et al. highlighted this peculiarity and was the first to show evidence of the band inversion in n-type Pb 0 77 . Sn 0 23 . Se, employing Nernst and thermopower experiments. 6 Particularly for device applications, a thorough understanding of the electronic properties of both n and p-type crystals is pivotal. 24

The temperature dependent band inversion of Pb 0 77 . Sn 0 23 . Se has also received significant attention from the topological insulator community. It was recently shown that while PbSe is a trivial insulator, SnSe is a topological crystalline insulator (TCI); a state of matter in which crystal point group symmetry protects the persistence of gapless surface states with a spin polarized Dirac dispersion. 25,27,29,30 Interestingly, in Pb 0 77 . Sn 0 23 . Se we have the unique opportunity to tune between topological and trivial phases by temperature, rather than doping, thus negating the effects of disorder. Moreover, since it is challenging to perform ARPES and STM (the primary methods of studying surface states) at various temperatures, temperature dependent broadband spectroscopy makes a great tool for the investigation of surface states in topological (crystalline) insulators. 24,27,31-39 We note that since band gaps are typically tuned by electron-phonon coupling and anharmonicity, the concurrent band inversion and topological phase transition hints at a more subtle process. Hence, Pb 0 77 . Sn 0 23 . Se can reveal new insights into how anharmonicity, electron-phonon coupling and topology all interplay.

In the present work, we measure the temperature dependent broadband reflectance of Pb 0 77 . Sn 0 23 . Se between 5 K - 292 K, and extract its optical conductivity and dielectric function between 6 meV - 6 eV. Subsequent extended Drude analysis reveals the frequency dependent free carrier scattering rate, and is used for a precise determination of the dominant temperature dependent scattering mechanism. We find clear optical evidence of the band inversion temperature at 160 ± 15 K, by studying the free carrier plasma frequency and direct optical band gap transitions. Interestingly, the temperature dependent free carrier scattering rate shows no signatures of the band inversion, and follows T 3 / 2 behaviour in accordance with electron-acoustic phonon scattering. Finally, we find that the Fermi level lies deep in the conduction band (114 ± 22 meV), obscuring the observation of surface states, 31,40-42 and causing the Seebeck coefficient and power factor to be insensitive to the band inversion.

(a)

FIG. 1. (a) The changing slope in the reflectance below the plasma edge (reflectance minimum) indicates a temperature dependent free carrier scattering rate in Pb 0 77 . Sn 0 23 . Se. (b) The optical conductivity reveals a clear separation ( ω c ) between Drude conductivity below 125 meV, and inter/intraband transitions above 125 meV. (c) The zero-crossing of the real part of the dielectric function ( glyph[epsilon1] 1 ), indicative of the screened plasma frequency, can be used to determine glyph[epsilon1] ∞

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0008v1-figure-ee916fd5-85d6-46d5-9f6b-c2707e3aa33e.jpg)

## II. EXPERIMENT

Pb 0 77 . Sn 0 23 . Se samples were grown at Princeton University. 6 Prior to the reflectance measurements, samples were freshly cleaved with a razor blade to produce a lustrous surface. Cleaving by tape was also attempted, but proved unsuccessful in producing high quality surfaces. Temperature dependent reflectance measurements were performed in near-normal geometry between 6 meV - 1.5 eV, using a modified Bruker VERTEX 80v and ARS continuous flow cryostat, described in detail elsewhere 31 (Fig. 1a). Reflectance data below 12 meV were corroborated by repeating the measurements using synchrotron radiation and a custom designed sample chamber with a closed cycle Cryomech cryostat at the AILES Beamline of the SOLEIL synchrotron in France. All infrared data were complemented by partially overlapping room temperature ellipsometry measurements between 0.75 meV -

6 eV.

Using RefFit, a simultaneous variational dielectric fit (VDF) of Reflectance and ellipsometry data yielded all Kramers-Kronig consistent optical constants. All measurements were performed on two different crystals from the same sample batch, showing excellent qualitative resemblance (see inset Fig. 1a for the reflectance of sample 2). Hence, data from only one crystal are discussed.

## III. RESULTS AND ANALYSIS

Fig. 1b shows the real part of the optical conductivity σ 1 ( ω ) of Pb 0 77 . Sn 0 23 . Se in the far and mid infrared. The spectra are dominated by a distinct Drude peak below 125 meV, associated with free carrier conductivity, and separated by a cut-off frequency ω c from the interband transitions above 125 meV.

The simplest approach to understanding the spectral range below ω c is to apply the semi-classical Drude model, in which the complex dielectric function for the free carrier response is given by

$$\hat { \epsilon } ( \omega ) = \epsilon _ { \infty } - \frac { \omega _ { p } ^ { 2 } } { \omega ^ { 2 } + i \frac { \omega } { \tau } } \quad \quad ( 1 ) \quad \$$

where ω p is the bare plasma frequency, τ is the transport lifetime, and glyph[epsilon1] ∞ captures all contributions to the dielectric function other than the Drude conductivity (i.e. interband transitions and polarizability of the static background). The bare plasma frequency is defined as ω 2 p = 4 πne /m 2 b , where n is the carrier density and m b is the band mass.

An immediate consequence of this model is the large, low frequency negative contribution to the real part of the dielectric function (ˆ = glyph[epsilon1] glyph[epsilon1] 1 + iglyph[epsilon1] 2 ) due to the free carrier response. Indeed from eq. 1 it can be seen that in the limit ωτ glyph[greatermuch] 1 the zero-crossing of glyph[epsilon1] 1 corresponds to the screened plasma frequency ˜ ω p , which is related to the bare plasma frequency ω p through ω p = ˜ ω p √ glyph[epsilon1] ∞ . Hence, we expect to observe this zero-crossing in glyph[epsilon1] 1 , which can be seen directly in 1c. The temperature dependence of both ω p and glyph[epsilon1] ∞ are discussed next.

Before analyzing the spectra in detail, it is useful to consider what we expect to observe as a result of the band inversion. For compounds in which the chemical potential is located in the bulk gap, a distinct signature of a band inversion would be the spectral weight redistribution from interband transitions to the Drude conductivity. As the gap closes, the onset frequency of the interband transitions (roughly ω c in Fig. 1c) moves to lower energy, until it merges with the Drude when the gap size reduces to zero. However, since stoichiometric Pb 0 77 . Sn 0 23 . Se is naturally n-type (chemical potential in the conduction band, shown in Fig. 2a), we do not expect to observe such obvious behaviour, and therefore perform a more sophisticated analysis to observe the band inversion. Fig. 2b illustrates how the band inver-

FIG. 2. (a) Band structure cartoon before or after the band inversion (left panel), and during the inversion (right panel), illustrating the allowed (in)direct transitions when the chemical potential is in the bulk conduction band (as is the case in as-grown Pb 0 77 . Sn 0 23 . Se). (b) Signatures of this inversion include a changing free carrier lifetime τ ( T, ω ) and plasma frequency ω p (region I), spectral weight change for indirect gap transitions (region II), and changes in the size of the direct gap (region III). (c) The plasma frequency of Pb 0 77 . Sn 0 23 . Se shows a distinct change in temperature dependence around 160 K, deviating from it's expected quadratic temperature dependence. (d) glyph[epsilon1] ∞ also shots an abrupt change in its temperature dependence at 160 K, further corroborating 160 K as the band inversion temperature.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0008v1-figure-dfd6bd14-0c37-4ce4-aa31-55fe40432e5d.jpg)

sion could be deduced from signatures in σ 1 ( ω ). Region I is dominated by free carrier behaviour, described by eqn. 1. Hence, as the band dispersion transforms across the band inversion (Fig. 2a), changes in both m b and n can result in an unusual temperature dependence of the plasma frequency (where ω 2 p ∝ spectral weight in σ 1 ( ω )). One might also expect the inversion driven changes in the electronic structure to be reflected in the carrier lifetime. This can be seen from the width of the Lorentzian features (centred at 0), while deviations from the Lorentzian lineshape indicate an energy dependent lifetime. Region III can be used to extract the magnitude of the direct gap transitions E dir shown in Fig. 2a. With the chemical potential ( µ ) in the conduction band, the smallest allowed E dir is offset from the conduction band minimum.

E #

in

I

III

IV

Hence, changes in the band dispersion as well as n are likely to affect E dir , resulting in a shift of the interband transition onset in σ 1 ( ω ), as shown in Fig. 2b. Since glyph[epsilon1] ∞ is a measure of the spectral weight in σ 1 above the free carriers, this quantity too should reveal signatures of the band inversion, similar to E dir . Finally, below E dir , region II is sensitive to temperature dependent changes in indirect gap transitions E in , where E in &lt; E dir , as well as impurity scattering responsible for an extended E dir onset (known as an Urbach tail 43 ).

## A. Bulk plasma frequency and glyph[epsilon1] ∞

The spectral weight of the Drude feature (area under σ 1 ) is a direct measure of the bare plasma frequency squared, ω 2 p , and was determined by integrating σ 1 ( ω ) up to a cut-off point ω c using

$$\omega _ { p } ^ { 2 } = \frac { 1 2 0 \Omega } { \pi } \int _ { 0 } ^ { \omega _ { c } } \sigma _ { 1 } ( \omega ^ { \prime } ) \, \mathrm d \omega ^ { \prime } \, \quad \quad ( 2 ) \quad \text{ has} \, \quad \text{a tc}$$

where ω c separates the Drude contribution to σ 1 ( ω ) from the interband transitions (as indicated in Fig. 1b). 20 The bare plasma frequency was then used to extract glyph[epsilon1] ∞ from the zero-crossing in glyph[epsilon1] 1 ( ω ). We note that an alternative method to find glyph[epsilon1] ∞ (which does not rely on the zerocrossing and is described elsewhere 20 ) was also used, and resulted in similar temperature dependent values differing only by a constant offset of a few percent.

The temperature dependence of ω 2 p and glyph[epsilon1] ∞ are shown in Fig. 2c and 2d, respectively. Starting with ω 2 p , we find a reduction as temperature increases, with a distinct discontinuity in the trend around 160 K. In general, a temperature dependent reduction in ω 2 p is an expected result that can be understood as follows. While semiclassically we have ω 2 p ∝ n/m , quantum mechanically the squared plasma frequency calculated using eqn. 2 corresponds to an integration of the density of states ( N ω,T ( )) multiplied by the Fermi function ( f ( ω, T )), so that ω 2 p ( T ) ∝ ∫ N ω,T f ω,T ( ) ( ) d ω . 44,45 Evaluating this integral using the Sommerfeld expansion yields a temperature independent term equal to ω 2 p at 0 K, and a negative T 2 term as a result of thermal smearing of the Fermi function. Hence, the temperature dependence of the squared plasma frequency is expected to follow

$$\omega _ { p } ^ { 2 } ( T ) \approx \omega _ { p } ^ { 2 } ( T = 0 K ) - B T ^ { 2 } \text{ \quad \ \ } ( 3 ) \text{ \ \ } \text{sur} \text{ }$$

where B depends on N ω,T ( ) at the Fermi energy. 44-46 Interestingly, while a linear decline in ω 2 p vs T 2 is expected, Fig. 2c shows two distinct regimes in its temperature dependence below and above 160 ± 15 K. Since B is a function of the density of states, this is a clear signature of an abrupt change in the band dispersion, associated with the band inversion. Indeed, this critical temperature falls within the 80 K - 200 K temperature range over which the band inversion was previously observed. 6,8,9,25,26,47 This also provides a first indication that despite the Fermi level of n-type Pb 0 77 . Sn 0 23 . Se not being in the gap, discontinuities in its optical properties are a useful probe in the exact determination of its band inversion temperature.

Besides clear signatures of the band inversion in ω 2 p , we also expect to see a change in glyph[epsilon1] ∞ at 160 ± 15 K. Since glyph[epsilon1] ∞ is a measure of the spectral weight of all interband transitions, it is sensitive to temperature dependent changes in the band structure. 48 This is confirmed by (Fig. 2d), revealing a clear change in the temperature dependence around 160 K, and will be discussed in greater in section III C.

## B. Surface States

Besides the Seebeck sensitivity to the band inversion at 160 ± 15 K, this critical temperature is also relevant to the topological phase of Pb 0 77 . Sn 0 23 . Se. Previous studies have shown how the band inversion is accompanied by a topological phase transition, where at low temperature topologically protected surface states persist in the bulk band gap. 6,25,26,29,47,49 In our optical data below 160 K (Fig.1a) we see no evidence of Dirac surface states, such as previously observed in topological insulators Bi 2 Te Se, 2 Bi Se , and (Bi,Sb) Te 2 3 2 3 31-36 This is not surprising since we find a bulk chemical potential of 114 ± 22 meV above the conduction band minimum at 5 K (shown later), consistent with previous studies of Pb 1 -x Sn Se with simix lar stoichiometry 6,25,26,47 Such a high bulk µ causes bulk conductance to dominate the optical properties, which can by understood by considering how the measured reflectance relates to the complex dielectric function, as R = (1 | -√ ˆ) glyph[epsilon1] / (1 + √ ˆ) glyph[epsilon1] | 2 . From Fig. 1c and eqn. 1 it is clear that | 1 ± √ ˆ glyph[epsilon1] | glyph[greatermuch] 1 for all frequencies ω &lt; ω ˜ p , resulting in R &gt; 95% over the full far infrared (FIR) range for most temperatures. Such high FIR reflectance even above the topological phase transitions obscures discriminating surface conductance signatures from the bulk conductivity in R ω ( ). Moreover, with the Fermi level so high in the conduction band, mixing of the surface and bulk bands also blurs other distinguishing features such as the bulk vs surface carrier lifetime. However, to ensure we did not overlook surface state signatures in our reflectance data, we used the Boltzmann transport equation to estimate their expected effect on the measured reflectance. Indeed, the surface state signatures are far below our detection sensitivity, as described in the Appendix, and crystals with a lower bulk carrier concentration are required for further optical investigations of the topological crystalline surface state properties of Pb 0 77 . Sn 0 23 . Se.

## C. Interband transitions

To study how the bulk electronic structure is affected by the band inversion, we turn to the band gap transi-

tions that are visible above ω c in Fig. 1b. With ARPES measurements showing µ in the conduction band for crystals with identical stoichiometry, 25,26,47 the smallest allowed direct gap transition exceeds the direct bandgap (and thus takes place away from the Brillouin Zone centre), as illustrated in Fig. 2a. This effect is known as the Burstein-Moss shift. 11,50 Assuming a parabolic dispersion for both the conduction and valence band, the energy of this transition can be estimated. By adding µ = ¯ h k 2 2 F / 2 m c to the direct band gap energy E g , and to the energy of the valence band electron with respect to valence band maximum of E v = ¯ h k 2 2 F / 2 m v , the total transition energy at zero energy is given by

$$E _ { d i r } \approx E _ { g } + \mu \left ( 1 + \frac { m _ { c } } { m _ { v } } \right ) \quad \quad ( 4 ) \quad \begin{array} { c } \gtrless \\ \tilde { \Xi } \end{array}$$

The magnitude of E dir can be obtained from a linear fit of the squared absorption coefficient, where α 2 = C ω ( -E dir ), as shown in Fig. 3a. 50 The temperature dependence of E dir is shown in Fig. 3b, where a clear maximum of 315 ± 10 meV is observed around 160 K. A previous study of the doping dependent band inversion in the Pb 1 -x Sn Te system also found a change in sign of x dE dir /dT as the bands invert. 14 Interestingly, for composition Pb 0 44 . Sn 0 56 . Te, in which a temperature dependent inversion takes place, a minimum in E dir is observed at the inversion temperature, in contrast to our maximum. Since the direct gap E g is zero at 160 ± 15 K, eqn. 4 suggests that this maximum is either due to an increase in µ , or a change in m /m c v across the transition. Both scenarios are plausible, and require further study specifically addressing the temperature dependent band mass and µ .

We find a minimum direct transition of 280 meV at 5 K, which is larger than the low temperature direct transition of 200-250 meV shown by previous ARPES results. 25,26,47 This difference is expected considering our larger carrier density compared to previous reports, as we show later in section DC Transport . We also note that the sensitivity to surface states in ARPES conceals bulk signatures, complicating the estimation of bulk transition energies.

We now consider the spectral weight below the direct transition energy E dir . This is often explained as an Urbach tail, which results from disorder in the crystal, and has an exponential energy dependence that changes monotonically with temperature. 11,43 While Fig. 3a does show exponential behaviour below E dir ( ∼ 280 meV), a closer look at α T ( ) at various frequencies in Fig. 3c reveals atypical Urbach behavior. Instead of a monotonic temperature dependent change of the exponential behaviour, a distinct minimum is observed at the 160 K band inversion temperature, while at both lower and higher temperatures the spectral weight in the absorption coefficient increases. Although it is not surprising that transitions below the direct gap are sensitive to the band inversion, the exact details of their temperature dependence are not clear and require future theory/experiment.

FIG. 3. (a) Linear fits of the squared absorption coefficient yield the smallest observable direct band transition energy (xintercept). (b) Direct gap transition extracted from the linear fit. (c) Absorption coefficient (log scale) below E dir , revealing a distinct change in temperature dependence, uncharacteristic for an Urbach tail.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0008v1-figure-024a8324-e23c-4d67-be6b-7bbad78c8e9e.jpg)

Nonetheless, this shows that the absorption/optical conductivity alone could be used to detect the transition.

## D. Extended Drude analysis and scattering rate

To study the bulk carrier lifetime in Pb 0 77 . Sn 0 23 . Se, we now turn to the extended Drude model. From eqn. 1, it can be seen that the free carrier response in the standard Drude model is fully defined by glyph[epsilon1] ∞ , ω p , and τ , where τ is the static transport lifetime resulting from a sum of all scattering processes according to Matthiessen's rule (including impurity, electron-phonon and electron-electron scattering). 51 While the Drude model is frequently used to study scattering mechanisms that can be distinguished by their temperature dependence, the temperature dependent band inversion in Pb 0 77 . Sn 0 23 . Se severely complicates such analysis. Moreover, the standard Drude model does not consider the energy dependence of the lifetime τ ( ω ). Since the carrier lifetime is an important parameter for thermoelectric efficiency 52 we turn to the extended Drude model, in which the lifetime has an energy dependence given by 16

$$\frac { 1 } { \tau ( \omega ) } = - \frac { \omega _ { p } ^ { 2 } } { \omega } \text{Im} \left ( \frac { 1 } { \hat { \epsilon } ( \omega ) - \epsilon _ { \infty } } \right ) \text{ \quad \ \ } ( 5 )$$

Fig. 4 shows the frequency dependent scattering rate obtained using eqn. 5. The black line shows where ω = 1 /τ , below which the quasiparticle transport is considered to

be coherent. 16 To study the nature of the dominant scattering mechanism, the data were fitted (red dotted lines in Fig. 4a) using the empirical relation

$$\frac { 1 } { \tau ( \omega ) } = \frac { 1 } { \tau _ { D C } } + \beta \omega ^ { n } \quad \quad \ \ ( 6 ) \quad \ \underbrace { \widehat { \in } } _ { \mathbb { S } }$$

where 1 /τ DC = αT m + 1 /τ imp . Here, 1 /τ imp accounts for impurity scattering, and α and β are constants that depends on material properties. When n = m = 2, and for certain values of α/β , eqn. 6 resembles the singleparticle scattering rate within Fermi Liquid theory, in which electron-electron scattering dominates the transport properties. 19

Fig. 4b shows 1 /τ DC , obtained from the fits with eqn. 6. As the temperature is raised, T 3 / 2 behavior is observed (dashed red line), associated with electronacoustic phonon scattering in semiconductors. 50,51 This is consistent with most good thermoelectrics, in which electron-acoustic phonon scattering dominates the transport properties. 7,53-55 Perhaps somewhat unexpected, however, is that the band inversion appears to have no bearing on 1 /τ DC . This is contrasted by recent work by L. Wu et al., 36 in which the scattering rate of carriers in (Bi 1 -x In x ) 2 Se 3 was found to diverge going across the topological phase transition as a function of doping.

We note that since we are most sensitive to bulk carriers (with µ high in the conduction band) we don't expect to see such a divergence, but it is worth emphasizing that more resistive Pb 1 -x Sn Se samples (in which the x surface contributions are enhanced) may exhibit similar behaviour. For such samples, the frequency dependent scattering rate may also serve as a novel analytical technique to discriminate surface states from the bulk. Since surface states are protected from backscattering, their lifetimes are expected to exceed the bulk carrier lifetime. This could result in distinct surface and bulk scattering contributions in 1 /τ ( ω ), for which we see no evidence in our samples. From fitting 1 /τ DC , we also find a constant 1 /τ imp = 6 8 . ± 0 3 meV, attributed to impurity scatter-. ing. We note that in metals and degenerate semiconductors such as as-grown Pb 0 77 . Sn 0 23 . Se, the electronic impurity scattering rate 1 /τ σ imp is roughly equal to the thermal impurity scattering lifetime 1 /τ κ imp . 51

Fig. 4c shows the temperature dependence of fitting parameters n and β . Since we were only able to obtain reliable 1 /τ ( ω ) data over a limited range spanning just one decade in energy (as seen in Fig. 4a), we hesitate to attribute much physical meaning to either parameter. Additional measurements at lower frequencies and temperatures beyond the scope of this work, or detailed studies by alternative probes such as high resolution temperature dependent ARPES could offer greater insights. Nevertheless, a considerable discontinuity in both n and β at 160 K, which does not fall within our error bars, does suggest that the frequency dependent scattering rate is tied to the band inversion. Moreover, with n = 2 normally associated with Fermi liquid type electron-electron dominated scattering, 19 this high value of n is quite sur-

(a)

Frequency (cm -1 )

β

FIG. 4. (a) Allometric fits of the frequency dependent scattering rate. The black dashed line corresponds to ω = 1 /τ . (b) A T 3 / 2 fit (red dashed line) of the temperate dependent static scattering rate 1 /τ DC illustrates phonon dominated scattering. (c) While it's difficult to interpret fitting parameters n and β due to the small frequency range over which 1 /τ could be fitted, they both reveal an abrupt change in temperature dependence around the band inversion temperature ∼ 160 K. (d) Comparison of Pb 0 77 . Sn 0 23 . Se resistivity obtained by optical spectroscopy (solid red line), and by standard DC transport methods by Liang et al. 6 (green dotted line), and Dixon et al. 8 (all other data).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0008v1-figure-16f0c325-c382-4e17-ba47-9be6c0fa8cf2.jpg)

prising considering the relatively high temperatures and frequencies up to which this behavior persists, making it an interesting starting point for followup studies.

## E. DC Transport

For comparison with previous transport studies of Pb 0 77 . Sn 0 23 . Se, we calculate its DC conductivity using σ o = ω τ 2 p DC / 60 and find σ o = 6 2 . × 10 3 Ω -1 cm -1 at 5 K. This low temperature value is in good agreement with previous transport results of n-type samples with identical stoichiometry (6 7 . × 10 3 Ω -1 cm -1 by Dixon et al. 8 and 5 5 . × 10 3 Ω -1 cm -1 by Liang et al. 6 ). A comparison of the resistivity with previously published data

is shown in Fig. 4d, where the dashed lines are digitized data from Liang et al. 6 (green) and Dixon et al. 8 (all other colors). We note that all n-type crystals increase in resistance with temperature, typical for metals and degenerate semiconductors. Only the resistivity of a low carrier p-type sample measured by Dixon et al. 8 shows an anomaly between 75 K and 200 K, possibly associated with the band inversion temperature, while none of the n-type curves reveal such signatures. This further illustrates the utility of optical spectroscopy, which can simultaneously determine ω p and 1 /τ to find σ 0 , while also measuring E F (via direct transitions) and the inversion temperature.

A possible explanation for the difference between ntype and p-type resistivity behavior is the asymmetry of the conduction and valance band, which can be tested by comparing their respective band masses. Using the carrier density of N e = 1 7 . × 10 19 cm -3 (from the resistivity curve most similar to our our results), we find a mobility of 2260 cm 2 V -1 s -1 , and a band mass of m b = 0 076m , . e using σ o = 4 πne τ ω /m 2 ( ) b = neµ . This band mass is smaller than the value of 0.1m e found by Dixon et al. for p-type Pb 0 77 . Sn 0 23 . Se, thus confirming the asymmetry between conduction and valence band. Using eqn. 4 we can now estimate the value of µ ( ≈ E F ), using the low temperature direct gap value of E g = 80 ± 20 meV, consistent with previous studies and calculations. 9,25,26 We find a value for µ of 114 ± 22 meV. To further corroborate this value, we derive an expression for E F using the parabolic band approximation. Starting with E F = ¯ h k 2 2 F / 2 m b , and substituting the 3D carrier density n 3 D = 1 (2 π ) 3 4 3 πk 3 F g g s v , we get

$$E _ { F } = \frac { \hbar { ^ { 2 } } } { 2 m _ { b } } \left ( \frac { 6 \pi ^ { 2 } n } { g _ { s } g _ { v } } \right ) ^ { 2 / 3 }$$

where g s = 2 and g v = 4 are the bulk spin and valley degeneracy, respectively. 6,25 Substituting the plasma frequency in SI units into eqn. 7, where ω 2 p = ne /m glyph[epsilon1] 2 b 0 , we find an expression for E F in [eV] as a function of physical constants, ω 2 p in [rad/s] and n in [m -3 ]:

$$E _ { F } = \frac { \hbar { ^ { 2 } } \omega _ { p } ^ { 2 } \epsilon _ { 0 } } { 2 e ^ { 3 } n ^ { 1 / 3 } } \left ( \frac { 6 \pi ^ { 2 } } { g _ { s } g _ { v } } \right ) ^ { 2 / 3 } \quad \quad ( 8 ) \quad \text{He} \\ \text{at}$$

Eqn. 8 yields a value for E F ( ≈ µ ) of 126 meV at 5 K, which is consistent with our value of µ of 114 ± 22 meV obtained from the direct gap transitions (see section. III C), and further confirms the metallic character of naturally doped n-type Pb 0 77 . Sn 0 23 . Se. As explained in Section III B, this further confirms why surface states are not observed in the Reflectance spectra.

## F. Power Factor

Finally, we relate our obtained values for σ o and µ to the Seebeck coefficient S T ( ) and the power factor

(a)

FIG. 5. (a) Seebeck coefficient in the parabolic band assumption, where all changes in E dir are attributed to changes in E F . (b) The Power Factor indicates minimal (if any) sensitivity to the band inversion at 160 ± 15 K, as a result of the high chemical potential.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0008v1-figure-434983fd-d058-4de1-8717-ebafb17e92e2.jpg)

S σ 2 o ( T ). In the parabolic band approximation, the Seebeck coefficient is given by 3 S = 8 π k 2 2 B 3 eh 2 m T b ( π 3 n ) 2 / 3 (where m b is the conduction of valance band mass), which can be expressed as a function of E F (using eqn. 7) as

$$S ( T ) = \left ( \frac { \pi ^ { 2 } k _ { B } } { 3 e } \right ) \left ( \frac { k _ { B } T } { E _ { F } } \right ) \left ( \frac { 1 } { g _ { s } g _ { v } } \right ) ^ { 2 / 3 } \quad ( 9 )$$

Since µ is located high in the conduction band, the direct gap transitions are likely only minimally sensitive to changes in the true band gap E g (and thus m c ) as a function of temperature. This is illustrated in the left and right panel in Fig. 2a, showing a similar band dispersion for high µ despite a dramatic change in the band gap. Hence, we can estimate the maximum change in S T ( ) by attributing the temperature dependence of E dir ( T ) (see eqn. 4) entirely to changes in µ T ( ). The result of this estimation, after substituting µ T ( ) for E F in eqn. 9, is plotted in Fig. 5a, while Fig. 5b shows the associated power factor S σ 2 0 . It is clear that neither S nor S σ 2 0 show distinct signatures of the band inversion at 160 ± 15 K. These results suggests that there is room for optimization of the figure of merit ZT in Pb 1 -x Sn Se compounds where x µ is located deep in the conduction band. Further doping to alter the lattice and suppress phonon contributions will likely only affect the electronic response minimally, as illustrated by its insensitivity to the band inversion and the unaffected power factor.

## IV. CONCLUSION

In conclusion, we have performed broadband optical spectroscopy experiments on as-grown n-type Pb 0 77 . Sn 0 23 . Se and extracted its optical constants over a temperature range between 5 K and 293 K. By studying the plasma frequency and band gap transitions, we find clear signatures of the temperature dependent band inversion at 160 ± 15 K, despite the chemical potential being far in the conduction band. We note that since temperature is an easily controllable variable in optical spectroscopy, in principle this could be used to optically determine band inversion temperatures of conductive samples with extremely high precision. While this band inversion is also associated with a topological phase transition, 6,25-27,29,47,49 we do not see any signatures of surface states, and thus don't comment on the topology of the band structure on either side of the band inversion. This is an expected result, however, as optical modelling suggests that the highly conductive bulk carriers dominate the properties of as-grown Pb 0 77 . Sn 0 23 . Se.

By performing extended Drude analysis, we accurately determine the temperature dependent scattering rate, and find that it follows T 3 / 2 behaviour (associated with electron-acoustic phonon interactions) typical for good thermoelectrics, 53 superimposed on a constant impurity scattering background with 1 /τ = 6 8 . ± 0 3 . meV. The frequency dependence of the scattering rate shows no evidence of different scattering channels below and above the band inversion, but data over a greater energy and temperature range are required for more conclusive results. Moreover, neither the DC resistivity, Seebeck coefficient, or the power factor show any evidence of the band inversion. This result confirms earlier transport studies on similar compounds, 9-12,22,23,28 and emphasizes the utility of optical spectroscopy in the study of thermoelectrics and/or band inversions.

For application of Pb 0 77 . Sn 0 23 . Se as either a topological crystalline insulator or thermoelectric, the high electrical conductivity of bulk states is a major prohibitive factor. Moreover, some of the fundamental physics governing the electronic, thermal, and topological properties are difficult to probe due to the high chemical potential in as-grown Pb 0 77 . Sn 0 23 . Se. Hence, future studies on samples in the Pb 1 -x Sn Se series, in which the carx rier density has been reduced, could lead to new insights

∗ reijnder@physics.utoronto.ca

† burchke@bc.edu

1 M. Zebarjadi, K. Esfarjani, M. S. Dresselhaus, Z. F. Ren, and G. Chen, Energy &amp; Environmental Science 5 , 5147 (2012).

2 S. B. Riffat and X. Ma, Applied Thermal Engineering 23 , 913 (2003).

into the interplay of phonons, thermal and electronic bulk lifetimes, and the optical and electronic properties of inversion symmetry protected topological surface states.

A reduction of the chemical potential would also enhance the sensitivity of the Seebeck coefficient to the band inversion, which could potentially result in a dramatic enhancement of the figure of merit ZT around the band inversion. Moreover, alternative stoichiometries also offer an interesting avenue, as previous measurements suggest that the band inversion in Pb 1 -x Sn Se x could be pushed to higher temperatures -relevant to many thermoelectric applications- by increasing the Sn content. 8,9,25 Hence, by optimizing Se doping, an ideal Pb 1 -x Sn Se stoichiometry can be determined, in which x ZT peaks around a desired temperature range. Indeed, such doping studies would also allow for the study of the high temperature interplay between the topologically nontrivial phase and thermoelectric efficiency, which is a direction currently unexplored.

With both thermoelectrics and topological (crystalline) insulators offering numerous appealing applications in industry, while presenting nontrivial physical characterization and optimization challenges, a better understanding of fundamental physics governing the figure of merit ZT and the isolation of surface conductance are both of practical and fundamental interest. The Pb 1 -x Sn Se and Pb x 1 -x Sn Te series offer many exciting x opportunities in this pursuit, particularly through optical investigations, by combining several novel physical phenomena into a single temperature tuneable system.

## V. ACKNOWLEDGEMENTS

We acknowledge Young-June Kim for use of an ellipsometer, as well as Luke Sandilands for very helpful discussions. This work has been partially funded by the Ontario Research Fund, the Natural Sciences and Engineering Research council of Canada, Canada Foundation for Innovation, and the Prins Bernhard Cultuurfonds. The crystal growth at Princeton University was supported by grant N6601-11-1-4110.

We also acknowledge the recent optical and Hall-effect study of Pb 0 77 . Sn 0 23 . Se by Anand et al. , 24 which was performed simultaneously, though independently. This study explores the temperature-driven band inversion of p-type Pb 0 77 . Sn 0 23 . Se, which is found to occur around 100 K.

- 3 G. J. Snyder and E. S. Toberer, Nature Materials 7 , 105 (2008).
- 4 F. J. DiSalvo, Science 285 , 703 (1999).
- 5 J. P. Heremans Font, Acta Physica Polonica A 108 , 609 (2006).
- 6 T. Liang, Q. Gibson, J. Xiong, M. Hirschberger, S. P. Koduvayur, R. J. Cava, and N. P. Ong, Nature Communica-

tions 4 , 2696 (2013).

- 7 O. Delaire, J. Ma, K. Marty, A. F. May, M. A. McGuire, M.-H. Du, D. J. Singh, A. Podlesnyak, G. Ehlers, M. D. Lumsden, and B. C. Sales, Nature Materials 10 , 614 (2011).
- 8 J. R. Dixon and G. F. Hoff, Physical Review B 3 , 4299 (1971).
- 9 A. J. Strauss, Physical Review 157 , 608 (1967).
- 10 M. Baleva, T. Georgiev, and G. Lashkarev, Journal of Physics: Condensed Matter 2 , 2935 (1990).
- 11 Z. M. Gibbs, A. LaLonde, and G. J. Snyder, New Journal of Physics 15 , 075020 (2013).
- 12 C. E. Ekuma, D. J. Singh, J. Moreno, and M. Jarrell, Physical Review B 85 , 085205 (2012).
- 13 J. Zhao, C. D. Malliakas, D. Bugaris, N. Appathurai, V. Karlapati, D. Y. Chung, M. G. Kanatzidis, and U. Chatterjee, arXiv (2014), 1404.1807v1.
- 14 S. O. Ferreira, E. Abramof, P. Motisuke, P. H. O. Rappl, H. Closs, A. Y. Ueta, C. Boschetti, and I. N. Bandeira, Brazilian Journal of Physics 29 , 771 (1999).
- 15 Z. M. Gibbs, H. Kim, H. Wang, R. L. White, F. Drymiotis, M. Kaviany, and G. Jeffrey Snyder, Applied Physics Letters 103 , 2109 (2013).
- 16 M. M. Qazilbash, J. J. Hamlin, R. E. Baumbach, L. Zhang, D. J. Singh, M. B. Maple, and D. N. Basov, Nature Physics 5 , 647 (2009).
- 17 D. N. Basov and T. Timusk, Reviews of Modern Physics 77 , 721 (2005).
- 18 D. Basov, R. Averitt, D. van der Marel, M. Dressel, and K. Haule, Reviews of Modern Physics 83 , 471 (2011).
- 19 U. Nagel, T. Uleksin, T. Room, R. P. S. M. Lobo, P. Lejay, C. C. Homes, J. S. Hall, A. W. Kinross, S. K. Purdy, T. Munsie, T. J. Williams, G. M. Luke, and T. Timusk, Proceedings of the National Academy of Sciences 109 , 19161 (2012).
- 20 J. Hwang, T. Timusk, and G. D. Gu, Journal of Physics: Condensed Matter 19 , 5208 (2007).
- 21 S. I. Mirzaei, D. Stricker, J. N. Hancock, C. Berthod, A. Georges, E. Van Heumen, M. K. Chan, X. Zhao, Y. Li, M. Greven, and D. van der Marel, Proceedings of the National Academy of Sciences 110 , 5774 (2013).
- 22 R. K. Shukla, S. K. Shukla, V. K. Pandey, and P. Awasthi, Physics and Chemistry of Liquids 45 , 169 (2007).
- 23 L.-D. Zhao, S.-H. Lo, Y. Zhang, H. Sun, G. Tan, C. Uher, C. Wolverton, V. P. Dravid, and M. G. Kanatzidis, Nature 508 , 373 (2014).
- 24 N. Anand, S. Buvaev, A. F. Hebard, D. B. Tanner, Z. Chen, Z. Li, K. Choudhary, A. B. Sinnott, G. Gu, and C. Martin, arXiv , 1407.5726 (2014).
- 25 P. Dziawa, B. J. Kowalski, K. Dybko, R. Buczko, A. Szczerbakow, M. Szot, E. glyph[suppress] Lusakowska, T. Balasubramanian, B. M. Wojek, M. H. Berntsen, O. Tjernberg, and T. Story, Nature Materials 11 , 1023 (2012).
- 26 B. M. Wojek, R. Buczko, S. Safaei, P. Dziawa, B. J. Kowalski, M. H. Berntsen, T. Balasubramanian, M. Leandersson, A. Szczerbakow, and P. Kacman, Physical Review B 87 , 115106 (2013).
- 27 X. Xi, X.-G. He, F. Guan, Z. Liu, R. D. Zhong, J. A. Schneeloch, T. S. Liu, G. D. Gu, D. Xu, Z. Chen, X. G. Hong, W. Ku, and G. L. Carr, arXiv , 1406.1726 (2014).
- 28 J. F. Butler, A. R. Calawa, and T. C. Harman, Applied Physics Letters 9 , 427 (1966).
- 29 L. Fu, Physical Review Letters 106 , 106802 (2011).
- 30 S.-Y. Xu, C. Liu, N. Alidoust, M. Neupane, D. Qian, I. Be-
- lopolski, J. D. Denlinger, Y. J. Wang, H. Lin, and L. A. Wray, Nature Communications 3 , 1192 (2012).
- 31 A. Reijnders, Y. Tian, L. J. Sandilands, G. Pohl, I. D. Kivlichan, S. Y. Frank Zhao, S. Jia, M. E. Charles, R. J. Cava, N. Alidoust, S. Xu, M. Neupane, M. Z. Hasan, X. Wang, S. W. Cheong, and K. S. Burch, Physical Review B 89 , 075138 (2014).
- 32 R. Vald´ es Aguilar, A. V. Stier, W. Liu, L. S. Bilbro, D. K. George, N. Bansal, L. Wu, J. Cerne, A. G. Markelz, S. Oh, ˇ and N. P. Armitage, Physical Review Letters 108 , 87403 (2012).
- 33 G. S. Jenkins, D. C. Schmadel, A. B. Sushkov, and H. D. Drew, Physical Review B 87 , 155126 (2013).
- 34 P. Di Pietro, M. Ortolani, O. Limaj, A. Di Gaspare, V. Giliberti, F. giorgianni, M. brahlek, N. bansal, N. koirala, S. Oh, P. Calvani, and S. lupi, Nature 8 , 556 (2013).
- 35 B. C. Chapler, M. Liu, K. Post, H. Stinson, M. Goldflam, A. Richardella, J. Lee, A. Reijnders, K. Burch, N. Samarth, and D. Basov, (unpublished) (2014).
- 36 L. Wu, M. Brahlek, R. V. Aguilar, A. V. Stier, C. M. Morris, Y. Lubashevsky, L. S. Bilbro, N. Bansal, S. Oh, and N. P. Armitage, Nature Physics 9 , 410 (2013).
- 37 Z. Li and J. P. Carbotte, Physical Review B 88 , 45414 (2013).
- 38 X. Xi, C. Ma, Z. Liu, Z. Chen, W. Ku, H. Berger, C. Martin, D. B. Tanner, and G. L. Carr, Physical Review Letters 111 , 155701 (2013).
- 39 M. K. Tran, J. Levallois, P. Lerch, J. Teyssier, A. B. Kuzmenko, G. Aut` es, O. V. Yazyev, A. Ubaldini, E. Giannini, D. van der Marel, and A. Akrap, Physical Review Letters 112 , 047402 (2014).
- 40 P. Di Pietro, F. M. Vitucci, D. Nicoletti, L. Baldassarre, P. Calvani, R. Cava, Y. S. Hor, U. SCHADE, and S. lupi, Physical Review B 86 , 045439 (2012).
- 41 A. D. LaForge, A. Frenzel, B. C. Pursley, T. Lin, X. Liu, J. Shi, and D. N. Basov, Physical Review B 81 (2010).
- 42 A. Akrap, M. Tran, A. Ubaldini, J. Teyssier, E. Giannini, D. van der Marel, P. Lerch, and C. C. Homes, Physical Review B 86 , 235207 (2012).
- 43 C. H. Grein and S. John, Physical Review B 39 , 1140 (1989).
- 44 L. Benfatto, S. Sharapov, N. Andrenacci, and H. Beck, Physical Review B 71 , 104511 (2005).
- 45 F. Carbone, A. Kuzmenko, H. Molegraaf, E. van Heumen, E. Giannini, and D. van der Marel, Physical Review B 74 , 024502 (2006).
- 46 M. Ortolani, P. Calvani, and S. Lupi, Physical Review Letters 94 , 67002 (2005).
- 47 A. Gyenis, I. K. Drozdov, S. Nadj-Perge, O. B. Jeong, J. Seo, I. Pletikosic, T. Valla, G. D. Gu, and A. Yazdani, arXiv: 1306.0043 (2013).
- 48 A. B. Kuzmenko, D. van der Marel, F. Carbone, and F. Marsiglio, New Journal of Physics 9 , 229 (2007).
- 49 Y. Okada, M. Serbyn, H. Lin, D. Walkup, W. Zhou, C. Dhital, M. Neupane, S. Xu, Y. J. Wang, R. Sankar, F. Chou, A. Bansil, M. Z. Hasan, S. D. Wilson, L. Fu, and V. Madhavan, Science , 2823 (2013).
- 50 P. Y. Yu and M. Cardona, Fundamentals of Semiconductors , 4th ed. (Springer, 2010).
- 51 J. Singleton, Band Theory and Electronic Properties of Solids (Oxford University Press, 2010).
- 52 M. Jonson and G. D. Mahan, Physical Review B 21 , 4223 (1980).

- 53 Y. Pei, X. Shi, A. LaLonde, H. Wang, L. Chen, and G. J. Snyder, Nature 473 , 66 (2011).
- 54 H. Wang, Y. Pei, A. D. LaLonde, and G. J. Snyder, Advanced Materials 23 , 1366 (2011).
- 55 Y. I. Ravich, B. A. Efimova, and V. I. Tamarchenko, physica status solidi (b) 43 , 11 (1971).
- 56 S. Adam, E. H. Hwang, V. M. Galitski, and S. Das Sarma, in Proceedings of the National Academy of Sciences (2007) pp. 18392-18397.
- 57 J. Liu, W. Duan, and L. Fu, Physical Review B 88 , 241303 (2013).

(a)

FIG. 6. (a) Reflectance of Pb 0 77 . Sn 0 23 . Se at 194 K, and a model of the 194 K Reflectance that included a surface state layer with ω p, 2 D = 0 34 meV cm . 1 / 2 . (b) Ratios of the temperature dependent reflectance above (black), below (green), and across (blue) the topological phase transition. The dotted red line is the ratio of reflectance spectra from (a), illustrating the much larger effect of temperature on R ω ( ) compared to the introduction of surface states below the topological phase transition.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0008v1-figure-50860171-84b1-4ae2-b20e-66abb4d54699.jpg)

## Appendix A: Expected optical signatures of TCI surface states in Pb 0 77 . Sn 0 23 . Se

To ensure we did not overlook surface state signatures in our reflectance data, we use the Boltzmann transport equation to estimate the expected effect of surface states on the measured reflectance. The 2D surface state plasma frequency can be approximated by ω 2 p, 2 D = ( e /h 2 )2 πk F v F g g s v , where v F and k F are the Fermi velocity and wavevector, and g s and g v are the spin and valley degeneracies. 31,35,56 Ignoring the parabolic correction for the Dirac states, we can use E F = ¯ hk F v F , and express ¯ hω p, 2 D in a convenient form with all SI units

$$\hbar { \omega } _ { p, 2 D } = c \sqrt { e E _ { F } g _ { s } g _ { v } 1 0 ^ { - 5 } } \quad \quad ( A 1 )$$

where c is the speed of light and e is the elementary charge. From previous ARPES measurements, we can estimate a surface Fermi level E ss F around 125 meV. 25 At the risk of overestimating the surface contribution, we can account for sample-to-sample variations by considering an even higher surface state Fermi level of 200 meV. We also note that while g s = 1 and g v = 2 near the Dirac point, additional Dirac states with smaller wavevectors intersect the Fermi surface when the Fermi level crosses the bulk conduction band, thus doubling the surface bands, albeit with different wavevectors. 6,25,49,57 For convenience, we account for these states by setting g v = 4, resulting in a overestimated 2D plasma frequency of 0.34 meV cm 1 / 2 .

Fig. 6a shows the reflectance of Pb 0 77 . Sn 0 23 . Se at 192 K (the topologically trivial regime) as well as a model in which a surface state layer with ¯ hω p, 2 D = 0 34 meV . cm 1 / 2 was added to the 192 K Reflectance. To emphasize the minimal difference, Fig. 6b shows the ratio of both reflectance spectra (red dotted line). The same plot also shows the ratio of reflectance measurements at different temperatures above (black), below (green), and across (blue) the topological phase transition. It is clear that the sharpening bulk plasma edge has a far greater effect on the reflectance, than the surface states. Hence, discriminating surface state signatures are not within our measurement capabilities for samples with such high bulk carrier densities.
Anjan A. Reijnders Jason Hamilton Vivian Britto Jean-Blaise Brubach Pascale Roy Quinn D. Gibson R. J. Cava K. S. Burch
Optical investigation of thermoelectric topological crystalline insulator Pb$_{0.77}$Sn$_{0.23}$Se
Pb$_{0.77}$Sn$_{0.23}$Se is a novel alloy of two promising thermoelectric
materials PbSe and SnSe that exhibits a temperature dependent band inversion
below 300 K. Recent work has shown that this band inversion also coincides with
a trivial to nontrivial topological phase transition. To understand how the
properties critical to thermoelectric efficiency are affected by the band
inversion, we measured the broadband optical response of
Pb$_{0.77}$Sn$_{0.23}$Se as a function of temperature. We find clear optical
evidence of the band inversion at $160\pm15$ K, and use the extended Drude
model to accurately determine a $T^{3/2}$ dependence of the bulk carrier
lifetime, associated with electron-acoustic phonon scattering. Due to the high
bulk carrier doping level, no discriminating signatures of the topological
surface states are found, although their presence cannot be excluded from our
data.