Chemical Tuning of Carrier Type and ... - ACS Publications

Jul 24, 2017 - Jülich−Aachen Research Alliance. (JARA-FIT and JARA-HPC), RWTH Aachen University, 52056 Aachen, Germany. •S Supporting Information...
2 downloads 9 Views 1MB Size
Subscriber access provided by UNIV OF NEWCASTLE

Article

Chemical Tuning of Carrier Type and Concentration in a Homologous Series of Crystalline Chalcogenides Tobias Schäfer, Philipp M. Konze, Jonas D. Huyeng, Volker L. Deringer, Thibault Lesieur, Paul Müller, Markus Morgenstern, Richard Dronskowski, and Matthias Wuttig Chem. Mater., Just Accepted Manuscript • DOI: 10.1021/acs.chemmater.7b01595 • Publication Date (Web): 24 Jul 2017 Downloaded from http://pubs.acs.org on July 24, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Chemistry of Materials is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 11

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chemistry of Materials

Chemical Tuning of Carrier Type and Concentration in a Homologous Series of Crystalline Chalcogenides Tobias Schäfer,† Philipp M. Konze,‡ Jonas D. Huyeng,† Volker L. Deringer,#,‡ Thibault Lesieur,† Paul Müller,‡ Markus Morgenstern,&,§ Richard Dronskowski,*,‡,§ and Matthias Wuttig*,†,§ †Institute

of Physics IA, ‡Institute of Inorganic Chemistry, &Institute of Physics IIB, and §Jülich–Aachen Research Alliance (JARA-FIT and JARA-HPC), RWTH Aachen University, 52056 Aachen, Germany. ABSTRACT: Tellurium-based phase-change materials (PCMs) enable applications from optical and electronic data storage to thermoelectrics and plasmonics, which all demand precise control of electronic properties. These materials contain an unusually large number of vacancies: “stoichiometric” ones that stem from the chemical composition, and “excess” vacancies that act like classical dopants. Here we show how both types of vacancies can be controlled independently in the solid solution Sn(Sb1–xBix)2Te4. We vary x in small steps over the entire compositional range and show that this has a profound effect on the material’s electronic nature: remarkably, we observe a change from p- to n-type conduction at x ≈ 0.7, solely controlled by composition. Our findings lead to a new compositionally (that is, chemically) tunable materials platform that enables precise control of electrical properties.

INTRODUCTION Chalcogenide materials enable a wealth of applications from digital data storage1, 2 to thermoelectrics3 and topological insulators.4-6 In all these, a decisive role is played by the materials’ electrical conductivity, which is proportional to the concentration n of charge carriers and to their mobility μ. Ultimately, independent control of both properties would be highly beneficial for applications. In the field of data storage, rapidly evolving research is concerned with phase-change random access memory (PCRAM), a non-volatile and fast memory technology.1, 2, 7-9 Many phase-change materials (PCMs) are quasibinary alloys of main-group IV and V chalcogenides, such as the widely studied Ge2Sb2Te5 and GeSb2Te4 (often colloquially referred to as “GSTs”). These materials are semiconductors for which one would expect, due to stoichiometry,10, 11 the Fermi energy EF to reside within the band gap. Experimentally, by stark contrast, they have been shown to be holeconducting (p-type) degenerate semiconductors with EF pinned in the valence band. There is a large carrier concentration in the samples, which first-principles computations traced back to a self-doping mechanism.12, 13 Besides, the materials are sensitive to disorder, which strongly influences carrier mobility. If carefully controlled, disorder can be exploited to tune electronic properties14-17—for example, by annealing samples at different temperatures and thus inducing different degrees of ordering. Materials such as the “GST” family with many empty lattice sites (vacancies) are particularly suitable for this. In applications, PCMs are switched between an amorphous phase (high resistance, “zero bit”) and a metastable rocksalt-like crystalline one (low resistance, “one bit”) by heating and quenching cycles. Further, extended heating of the metastable phases leads to structural reordering and results in stable, layered trigonal structures.18 The second aforementioned field where chalcogenides are used is thermoelectric power generation (waste-heat

recovery).3 By construction, thermoelectric devices require both p- and n-type materials, of which the former are generally more abundant. Materials are typically alloyed or nanostructured to lower the thermal conductivity (“electron crystal, phonon glass” concept; ref 19). Again, disorder in the crystalline phases has substantial impact on thermal transport,20, 21 which is relevant for thermoelectric applications. Besides low thermal conductivity, a high Seebeck coefficient S and large electrical conductivity  are mandatory: the carrier concentration (that influences all three quantities) must be chosen to yield a good trade-off. An important difference compared to data-storage applications (switching between amorphous and metastablecrystalline states) is that thermoelectric materials are chiefly used in their thermodynamically stable crystal phase. This is because the application exposes devices and materials to extended periods of high temperature. Finally, in recent years, several telluride compounds have been predicted and proven to develop topologically non-trivial or topologically insulating (TI) surface states.4-6, 22-27 Here, bulk conductivity is regarded a parasitic effect and neither p- nor n-type bulk conductivity is desired, as both self-doping and disorder are unwanted sources of carriers. Hence, for TI-based devices,26 doping approaches are employed to reduce conductivity in the bulk: the aim is different, but once more a careful control of the electronic nature is crucial to enable future applications. In this work, we move from the concept of disorder tuning known in the PCM community to a more general picture including ideas from thermoelectrics and TI research, and we show how chemical control (that is, varying the composition rather than physical parameters) is crucial to enable new properties. We synthesize and study the pseudo-ternary compounds Sn(Sb1–xBix)2Te4; the end members of this pseudo-ternary line (SnSb2Te4 and SnBi2Te4) have both been used in thermoelectrics and shown to develop TI surface states.28-31 By adjusting the chemical composition, we are able to continuously tune the electronic nature

ACS Paragon Plus Environment

Chemistry of Materials

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

of the solid solution and, ultimately, to turn it from a ptype into an n-type conductor. This provides new general insight into chalcogenide materials and helps to realize a chemically tunable materials platform.

Page 2 of 11

by step, and X-ray patterns are recorded at each temperature. SnSb2Te4, SnSbBiTe4, and SnBi2Te4 retain their metastable structure up to 280±5, 270±5, 275±5 °C, respectively, followed by a relatively sharp transition.

RESULTS AND DISCUSSION Synthesis and Structure of Sn(Sb1–xBix)2Te4. The materials investigated here are heavier homologues of the aforementioned PCMs. One of the parent compounds, SnSb2Te4, is isostructural to GeSb2Te4 and similar in chemical and electronic properties.10, 32-40 From there on, one may substitute the main-group V element, by moving to the heavier homologue, formally written as SnSb2Te4 + 2x Bi → Sn(Sb1–xBix)2Te4 + 2x Sb. In a first step, we synthesized three powder samples, namely, SnSb2Te4 (x = 0), a material close to SnBiSbTe4 (x = 0.57 as determined by EDX), and SnBi2Te4 (x = 1), and investigated their structures. Sputter deposition initially yields rocksalt-like phases in all cases, and powder XRD patterns and refinement results are provided as Supporting Information (Figure S2 and Tables S1–2). As expected, the lattice parameter increases from 6.1906(4) Å (x = 0) via 6.224(5) Å (x = 0.5) to 6.3251(4) Å (x = 1). Upon heating, the samples transform to stable trigonal polymorphs. Here, the in-plane lattice parameter a rises from 4.3158(1) Å via 4.3711(1) Å to 4.4085(1) Å, while the out-of-plane parameter c stays almost constant upon Sb→Bi replacement (to within 0.02 Å). This anisotropy can be explained by first-principles computations, showing that the octahedral distortion around Sn is lower in SnBi2Te4 than in SnSb2Te4 (due to the difference in ionic radii, Figure S1); in consequence, the height of the septuple-layer building block remains largely unchanged. The structures of the end members SnSb2Te4 and SnBi2Te4 agree well with literature data.33, 35, 39, 41, 42 No syntheses or structural data of rocksalt-type SnBi2Te4 and of the intermediate phases have been reported to our knowledge. This transformation from a rocksalt-like metastable to a trigonal stable structure is in accord with previous observations for GeSb2Te4, which suggests similar principles for the occupation of the lattice. Assuming that GeSb2Te4, SnSb2Te4 and SnBi2Te4 are structurally analogous, one can infer the following occupancy of the crystalline unit cell: in the metastable rocksalt-type phase, the anionic-like sublattice is fully occupied with Te atoms, whereas on the interpenetrating second sublattice, Sn, Sb, and Bi, as well as 25% vacancies are statistically distributed.43 A different way to describe the structure is as a dense stacking of 111 planes, of which every second is filled with Te, while those in between are occupied by Sn, Sb, Bi, and vacancies. Upon heating, the vacancies assemble into layers along 111, concomitant with increased cation ordering.14, 44 While this vacancy ordering proceeds continuously during heating, the rocksalt-to-trigonal transition is furthermore characterized by a displacement of the septuple-layers (yielding the long c-axis) and by a narrowing of the vacancy layers. This causes the denser packing in the trigonal form (smaller volume in Figure 1c). This process can be directly tracked in XRD experiments. A capillary filled with non-annealed powder is heated step

Figure 1. (a) Structure models of metastable (left) and stable (right) Sn(Sb1–xBix)2Te4, with exemplary environments of Te shown as close-ups. Sn, Sb, and Bi are dark grey, red, and blue, respectively; Te atoms are light grey, and an exemplary vacancy site on the rocksalt-like lattice is indicated (□). (b) Grazing-incidence XRD for the homologous series. The scans on the left represent the cubic “as-deposited” (that is, non-annealed) state, while the samples characterized on the right have been annealed to 325 °C. The Bi content increases from bottom (red; pure SnSb2Te4) to top (blue; pure SnBi2Te4). (c) Application of Vegard’s Law, comparing experimental volumes per formula unit Sn(Sb1–xBix)2Te4 (additional powder samples as open circles) to computed ones (Methods section).

Summarizing thus far, we synthesized Sn(Sb1–xBix)2Te4 by an alternating co-deposition sputtering process. The samples initially have rocksalt-like structures, and on annealing convert to stable trigonal phases. Our data so far

ACS Paragon Plus Environment

Page 3 of 11

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chemistry of Materials

suggest that Sn(Sb1–xBix)2Te4 behaves similarly to the wellknown GeSb2Te4 in terms of structure. The subsequent discussion will center on electrical data recorded on thin films (~250 nm), rather than bulk samples. Therefore, thin-film samples were synthesized and verified to exhibit similar structure: grazing-incidence XRD patterns for the thin films (Figure 1b) agree well with powder-diffraction results (Figure S2) and evidence that rocksalt-like as well as trigonal structures were successfully produced along the entire line. We note in passing that the lattice parameter of the as-deposited thin-film samples varies slightly compared to that determined for the bulk powder, due to necessary adjustments in the production process; this is much less pronounced for the stable structures. To verify Vegard-like behavior, the experimentally determined cell volumes are shown in Figure 1c. We furthermore performed first-principles computations for structural models of the quasi-ternary series: the results are well in line with the thin-film data (Figure 1c) and discussed in the Supporting Information. In particular, no pronounced structural rearrangements were observed during computations: one can hence conclude complete mixing along the Sn(Sb1–xBix)2Te4 line; this holds for both the rocksalt and the trigonal phases, and for all samples produced in this work. The composition of all samples was verified by EDX, yielding just a minor Te deficit independent of  (Figure S5). This is common for sputter-deposited PCMs.45 How Disorder and Doping Control Electronic Properties. So far, we have shown that Sn(Sb1–xBix)2Te4 is structurally similar to typical PCMs such as GeSb2Te4—what, now, about the electronic nature of both material systems? As reported before, the electronic properties of PCMs are rather intricate and dominated by the distribution of vacancies. In the metastable phase, vacancies randomly occupy the cation lattice, but upon annealing they order into clusters or planes.14 This strongly influences the localization of the electronic wavefunction: ordered vacancy planes do not interrupt extended Bloch states in the crystal, but disordered vacancies do.14 Besides providing scattering centers, they can also generate localized states that mostly concentrate on the band edges. Localized states at the band edges are mostly known from amorphous semiconductors, but crystalline PCMs can feature localized states for at least the first 130 meV from the band edge.17 These localized states at the band edge alter the transport mechanism, from a Drudelike propagation of Bloch waves to a temperatureactivated transport regime. While the former is the typical transport mechanism in metals, providing conductivity by nearly free electrons, the latter relies on carriers localized in the defect states. To participate in transport, those localized carriers either get thermally activated into a delocalized state, or the carrier tunnels from one localized state to another (“hopping”). This is known as Anderson localization and has been observed in PCMs.17, 46 Both in the localized and delocalized transport regime, many PCMs have extremely low mobilities of the order of 0.2 to 25 cm2/Vs, which are more common for amorphous than for crystalline materials.

There are two different mechanisms for vacancy formation in PCMs. We will, in the following, speak of “stoichiometric” and “excess” vacancies (Figure 2); the former term is used for the empty cation sites arising directly from composition (on the cation sublattice of GeSb2Te4, one quarter of the sites is empty). The formation of these in GeSb2Te4 was previously explained using chemicalbonding analyses:11 vacancy formation lowers EF, depleting previously filled, antibonding regions of the band structure. In addition, there are “excess vacancies”, less abundant in numbers, but it is those vacancies that cause p-type conductivity (approximately one percent of cation sites in GeSb2Te4 and GeTe; Supporting Information). The occurrence of excess vacancies is often referred to as “selfdoping”: the material adjusts the number of carriers without external influence, and a self-doped solid often is insensitive to additional doping. Note that this is also irrespective of the sputtering process, which induces a minor Te deficiency (see above) which would intuitively be expected to lead to Te vacancies; instead, the formation of Ge vacancies is energetically much more favorable, resulting in p-type conductivity.12

Figure 2. Simplified schematic of how disorder and selfdoping control the electrical conductivity in crystalline PCMs. In the metastable phases (left), both stoichiometric and excess vacancies are present; they localize carriers, reducing the mobility  (the stoichiometric vacancies dominate due to higher abundance). In the stable phases (right), vacancies arrange into well-ordered van der Waals gaps, and only the excess vacancies remain. This allows one to separate both “types” of vacancies and to study them independently.

A different perspective on the effect of excess vacancies is that they shift the Fermi energy EF into one band (for ptype materials, the valence band). This renders the materi-

ACS Paragon Plus Environment

Chemistry of Materials

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

al a degenerate semiconductor, which in turn generates carriers. Both views describe the same underlying phenomenon despite focusing on different aspects. The two “types” of vacancies are encoded by separate colors in Figure 2, although it is of course impossible to assign an individual vacancy to one “type”; all implications of the diagram should be in a statistical manner. While only excess vacancies (green) provide carriers, all vacancies may serve as scattering centres if randomly distributed. Due to their overwhelming number, stoichiometric vacancies (purple) govern disorder and localization, and hence mobility in the metastable phase. During annealing, the stoichiometric vacancies order into planes, enhancing mobility. Eventually, upon transition to the stable phase, the vacancies form vdW gaps which renders them irrelevant for mobility tuning. Consequently, there are two approaches to tailor conductivity: either, to control the mobility by changing the ordering17 and/or number45 of stoichiometric vacancies, or, to vary the number of excess vacancies and hence carrier concentration. An approach for the latter is presented in the following. Figure 3 shows the resistivity during annealing of the parent compounds SnSb2Te4 and SnBi2Te4. In this experiment, the same sample has been heated to progressively higher temperatures and cooled down to room temperature after each step. The results for SnSb2Te4 resemble earlier ones for GeSb2Te4.17 With each annealing step, the resistance drops further, and the temperature coefficient of resistivity (TCR = dρ / dT) gradually changes from temperature-activated insulator-like (TCR < 0, indicating Anderson localization) to metal-like behavior (TCR > 0).47 Hall measurements at room temperature of the asdeposited state (point (A) in Figure 3), an intermediate (B), and a highly annealed one (C) reveal that the drop in resistivity (increase in conductivity) is solely due to increased mobility (from 0.1 to 13 cm2 /Vs), while the carrier concentration does not change considerably (Figure S3). This is consistent with previous findings for GeSb2Te4.17 In Figure 2, it would be visualized as a gradual ordering of stoichiometric vacancies. As expected for an Anderson-localized material, similar to GeSb2Te4, SnSb2Te4 shows low mobilities compared to typical crystalline materials.17 During annealing, both conductivity and crystal structure change but in a fundamentally different way. While the change from rocksalt to trigonal structure does not occur below 270 °C, the increase in conductivity progresses throughout the entire annealing process (Figure 3) without any specific correlation with structural changes. Both effects thus go back to separate processes, even though they are both linked to vacancy ordering.

Page 4 of 11

Figure 3. Resistivity of the stoichiometric end members during progressive annealing. The resistivity of SnSb2Te4 (a) decreases upon stepwise annealing at increasingly large temperatures, similar to Ge–Sb–Te PCMs (see text). By contrast, that of SnBi2Te4 (b) stays almost constant.

Compared to SnSb2Te4, the resistivity of SnBi2Te4 behaves very differently (Figure 3b): it shows only a moderate annealing effect, and its resistivity increases during annealing. However, again decoupling carrier concentration and mobility, one can establish an increase of mobility upon annealing (from 1.6 to 6.6 cm2 /Vs; Figure S3). This remains in line with an ordering of stoichiometric vacancies that promotes mobility and eventually leads to a structural transition as well. But here, the increase in mobility is overcompensated by a severe reduction in carrier concentration, from 5 10 cm (A) to 7 10  cm (C) at room temperature. This might stem from the excess vacancies being metastable upon annealing, and hence a synchronous and countervailing tuning of carriers and mobility. Anderson physics and its relief by annealing play a role in both SnSb2Te4 and SnBi2Te4 as indicated by the low mobilities, however the initially high number of carriers in SnBi2Te4 overshadows this effect by pushing the Fermilevel deep into the band. As stated above, the thin-film samples have a thickness of around 250 nm, and in principle finite-size effects can influence electronic property measurements. However, even for our highest-mobility samples, the mean free path stays well below 2 nm, more than two orders of magnitude

ACS Paragon Plus Environment

Page 5 of 11

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chemistry of Materials

below the film thickness. Our samples can therefore be regarded as bulk-like in terms of transport properties. Tuning Carrier Concentration, the Chemical Way. To decouple the influence of both vacancy types, we next aim to vary the number of excess vacancies only. All compounds under study have isovalent IV1V2VI4 compositions, so the concentration of stoichiometric vacancies is constant. As shown in Figure 4a, the resistivity of the asdeposited samples (purple) continuously decreases with increasing Bi content, and the initial mobility of SnBi2Te4 is 16× that of SnSb2Te4 as discussed before. The TCR becomes less insulating-like at the same time; it is therefore obvious that the ordering of stoichiometric vacancies plays an important role.

Figure 4. Changes in electronic properties with chemical composition. (a) Electrical resistivity at 300 K with increasing Bi content for as-deposited cubic (purple) and annealed trigonal (green) samples. (b) Temperature coefficient of resistivity (TCR) for as-deposited and annealed samples at 300 K, respectively. A positive value indicates metal-like behavior (resistivity rising with temperature), while negative values point toward a temperature-activated conduction mechanism.

We are here concerned with compositional, chemical control, and it is therefore desirable to eliminate the influence of vacancy disorder (which overshadows the role of the less abundant excess vacancies). For this, the annealed samples in their stable phase are more promising subjects of study: here, no “stoichiometric vacancies” can hamper the mobility (Figure 2, right). In this case, lower resistivity is seen at both ends of the stoichiometric line than in the mixed region (green in Figure 4a). The end members exhibit a metal-like TCR, while samples at intermediate x show insulator-like resistivity.

Resistivity values alone cannot disentangle the influence of disorder from a change in number of carriers. Therefore, Hall and Seebeck effect measurements are needed for deeper insight (Figure 5). While SnSb2Te4 is holeconducting (p-type) with a high carrier concentration on the order of 1021 cm–3, SnBi2Te4 is an electron conductor (n-type) with comparable carrier density, and hence the majority carrier type changes at an intermediate composition (around x = 0.7). We discuss these three regions in sequence, beginning with the p-type one (0   ≪ 0.7). Starting from p-type SnSb2Te4 and moving to larger x, the number of carriers decreases as does the conductivity. Interpreting this in the framework of Figure 2 suggests a reduction in “excess vacancies” by Bi introduction. Roomtemperature and low-temperature Hall data are nearly superimposable: together with absolute values of 1020−1021 cm−3, this is typical for degenerate semiconductors. This indicates that the Fermi energy of the semiconducting material has been shifted into the valence band. Almost no carriers are generated by increasing the measurement temperature; hence the transport is not temperature-activated, as already indicated by the TCR (Figure 4b). A Drude-like transport model is indeed suitable here. Location of the Fermi Energy. So far, the location of the Fermi energy has been inferred qualitatively, but Seebeck measurements (Figure 5b) can also provide a direct, quantitative determination of EF in a Drude-like regime. In the textbook case of EF in a single isotropic parabolic band, the Seebeck coefficient S can be used to calculate EF using   3  " ||    ! , |#$ |  3 2 where is chosen depending on the scattering mechanism (the most common values chosen are  &½ for acoustic phonons, and  1.5 for ionized impurities). As the temperature dependence of resistivity and Seebeck coefficient in our materials is metal-like (Figure 6), the scattering at room temperature is dominated by phonons, and γ = −½ was chosen accordingly.48 The overall term in parentheses in the above equation is positive, and so the absolute Seebeck coefficient rises with increasing temperature. Starting from SnSb2Te4 and adding Bi, the magnitude of S rises, indicating a reduction in carrier concentration. The temperature dependence of this allows additional insight: in metals (and well-ordered degenerate semiconductors), S increases with temperature. This is visible here as well, justifying the model introduced above to calculate EF. As shown in Figure 5c, EF in SnSb2Te4 is located 0.13 eV below the valence-band edge comparable to GeSb2Te4. On Sb→Bi substitution, it shifts continuously toward the valenceband edge. Starting from the other endpoint, SnBi2Te4, a similar picture evolves: moving toward the p-n-transition, EF approaches the band edge and lowers the number of carriers and the conductivity. In SnBi2Te4, however, EF resides about 0.1 eV in the conduction band, and the majority carriers are n-type (Figure 5). The claims made here for two selected temperatures can be further strengthened by detailed temperature-dependent transport data (Figure 6).

ACS Paragon Plus Environment

Chemistry of Materials

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Having examined the regions adjacent to both endpoints of the compositional line, a discussion of the direct vicinity of the p-to-n-transition ( ' 0.7( should be provided as well. Here, unlike elsewhere, samples do show signatures of temperature activation in their carrier concentrations, as typical for ordinary semiconductors. Additionally, the carrier concentration at x=0.7 follows the trend of the surrounding compositions only at low temperatures. By contrast, at room temperature, the carrier density exceeds that in Bi-richer samples. This might be due to the excitation of carriers into a second band, so that electrons and holes compensate in the Hall signal, which requires a small band gap. In fact, estimates from the Seebeck coefficient point to 60-90 meV (Supporting Information), which agrees with findings in the related materials series (Sb,Bi)2Te3.24, 49-51

Page 6 of 11

ties rise. Although the standard reasoning here is made for band insulators, the predictions likewise hold for Anderson insulators.52, 53 Additionally, the magnitude of the Seebeck coefficient is lower than expected from neighbouring compositions. This coincides with the overestimated Hall carrier concentration and might be an additional sign of a two-channel transport, where both electrons and holes participate (and partly cancel each other) in Hall and Seebeck coefficient. Recalling the low mean free path (< 2 nm) and mobility in these samples, however, puts forward that disorder is relevant in transport as well, as soon as the Fermi-energy approaches the band edge, where localized states are most prominent. Indeed, the data plotted against the inverse temperature (Figure S4) seems to saturate at low temperatures, being unexpected both for bandinsulators and two-carrier transport. Hence, a reasonable model for the Seebeck coefficient should include both disorder (cf. refs 52, 53) and two bands, which is certainly not feasible on the limited available temperature range.

Figure 5. Experimental verification of the p-to-n transition. For the annealed homologous series, measurements are reported of (a) Hall carrier concentration nH, (b) Seebeck coefficient S, and (c) location of the Fermi energy #$ , as calculated from S. In the latter case, negative values indicate that #$ is located in the valence band. Data measured at room temperature (green) and low temperature (grey) are shown.

The data in direct vicinity of the transition shows remarkable features in the Seebeck coefficient as well. The temperature coefficient decreases with temperature, suggesting a transition from the degenerate to the nondegenerate case, or disorder localizing the carriers at reduced temperatures. The Seebeck coefficient of insulators drops with increasing temperature as their carrier densi-

ACS Paragon Plus Environment

Page 7 of 11

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chemistry of Materials

Figure 6. Temperature dependent resistivity, Hall and Seebeck data for the samples shown in Figure 4. (An additional plot of the same data against 1/T is found in Figure S4.)

The data around the transition can be well explained by carriers that can be excited into the conduction band at moderate temperatures already, as the band gap estimated from Seebeck is only around 3 kBT at room temperature. However, findings in some related materials show that the p-to-n-transition is accompanied by a band inversion, which closes and reopens the band gap. While the latter was found in BixSb(1-x),54, 55 no band inversion has been observed in the more closely related Bi2xSb2(1-x)Te3 system.26, 56 A definite distinction of both scenarios cannot be made with electrical data at moderately low temperatures. Hall and Seebeck data at low cryogenic temperatures might help to differentiate between both cases, but precise measurements or simulations of the band structure around the transition are more promising. Both is beyond the scope of this work. Although both chemical substitution (as done here) and the disorder-driven Anderson MIT described previously in PCMs17 decrease the conductivity, these mechanisms are fundamentally different in origin. The Anderson transition is governed by carrier mobility (vacancy ordering) but leaves the absolute number of carriers largely unchanged, whereas we here observe a drastic change in carrier concentration (Figure 5a). In the spirit of Figure 2, the Anderson MIT can be explained by a movement of the mobility edge #) (which separates disorder-localized states from lower-lying extended ones). By contrast, the quantity that shifts here is EF: in SnSb2Te4, it lies well inside the valence band, but Sb→Bi substitution pushes it gradually upward. Around x=0.7, the band gap is crossed, followed by EF continuing its upward movement into the conduction band (Figure 5c). Understanding the mechanisms underlying this behavior will require further, microscopic insight. For example, the calculation of defect formation energies for charged defects has enabled insight into the conduction mechanisms in complex chalcogenides.57 However, this is beyond the scope of this work, owing to the complexity of the fourcomponent system studied here. Additional experimental insight might come from advanced synthesis techniques (e.g., the growth of large single crystals, which would allow to go beyond powder X-ray diffraction), and from local analytical probes such as local energy-dispersive X-ray spectroscopy, or Mößbauer spectroscopy.

CONCLUSIONS The homologous series Sn(Sb1–xBix)2Te4 can be tuned in carrier density by more than one order of magnitude, and exhibits reversal from p- to n-type conductivity, solely by modifying the chemical composition. This is due to a shift of the Fermi edge from deep inside the valence band (SnSb2Te4) toward and through the gap, and finally into the conduction band (SnBi2Te4). While previous studies focused on modifying the mobility edge, we here showed how the Fermi edge can be systematically controlled in crystalline chalcogenides. However, in contrast to known p-n-transitions in main-group IV and V chalcogenides, our

particular system allows for mobility tuning via Anderson localization as well. Our work therefore adds a key aspect to the toolbox for tailoring electronic properties of chalcogenide materials, which promises new applications in thermoelectric or electronic devices—for example, in neuromorphic hardware, where distinct and carefully tailored electronic states are needed to encode digital information beyond simple “ones” and “zeroes”. Finally, as TI surface states have been predicted both for SnBi2Te4 and SnSb2Te4,28-31 a carefully synthesized and highly ordered intermediate species might even become relevant as a bulk-insulating TI.

METHODS Sample Preparation. All samples were prepared by sputter deposition. As substrates for resistivity and Hall effect measurements, simple 2 cm × 2 cm glass cover slips (MenzelGläser) were used, whereas Corning 1737F glass was employed for Seebeck measurements. Contacts for electrical measurements were deposited on the substrates prior to film deposition by thermal evaporation of gold (90 nm) on a 10 nm chromium adhesion layer. Direct-current magnetron sputtering was performed in a customized Von Ardenne LS 320S sputter system. The chamber was evacuated to achieve a base pressure lower than 2×10–6 mbar prior to deposition. Ar was introduced as a process gas with a pressure in the range of 7×10–3 mbar. The target-to-substrate distance was 4.5 cm; samples were revolved over the targets to achieve better homogeneity. SnSb2Te4 and SnBi2Te4 were deposited from stoichiometric targets, while the homologous series was sputtered via codeposition of both by adjusting the sputter power. The films were ≈ 250 nm thick, and all were capped in situ with a 25 nm amorphous ZnS:SiO2 (80:20) cover layer via RF-magnetron sputtering. For powder samples, we deposited layers of several μm thickness on spring steel sheets and delaminated them due to their poor adhesion. Thermal treatment for both thin-film and powder samples was performed under Ar atmosphere, using ramps of 5 K min– 1 and holding times of 30 min. Characterization. Crystal structures of samples were analyzed using X-ray diffraction with three different setups. First, powder samples were measured in the as-deposited and an ex situ annealed state using a STADI P (STOE Darmstadt) powder diffractometer with Cu Kα1 radiation and an image plate detector collecting 140° in 2θ simultaneously. Rietveld refinements were carried out using FullProf.58 Second, powder samples were annealed to higher temperatures in steps of 5 K, while measuring their diffraction pattern for 5 min at each holding temperature. For this, we used a G644 (HUBER Rimsting) powder diffractometer with Cu Kα1 radiation. Finally, a Panalytical X’Pert Pro diffractometer using Cu Kα-radiation was used for thin-film characterization in a grazing-incidence geometry with ω=0.7°. Resistivity measurements were done on square-shaped samples in a custom-built setup employing the Van der Pauw method.59, 60 The Hall carrier density was evaluated with an alternating current (AC) Lock-In technique, whereby the current (at 5 Hz) as well as the applied magnetic field (±0.2 T with 2 Hz) were modulated.61, 62 Two consistent Hall signals were demodulated at 3 Hz and 7 Hz. Cooling with liquid nitrogen allowed for low-temperature measurements. The samples

ACS Paragon Plus Environment

Chemistry of Materials

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

used for Hall measurements were the same as those used for resistivity measurements, so full consistency of the experimental results is ensured. Seebeck measurements were conducted in a custom setup as well, using Peltier stacks for temperature control and an Agilent 34410A DMM for measuring. The temperature differences between both ends of the sample were ±5 K and 0 K, respectively. All electrical measurements were carried out in vacuum. Consistency of all data is ensured since both the Hall and the Seebeck setup can record temperature-dependent two-point resistances. First-Principles Computations. Density-functional theory (DFT) computations were carried out to complement the experiment. To properly treat Te···Te interactions in the layered structures, the “D3” dispersion correction63 to the PBE functional64 was employed, which reproduces the title compounds’ structural parameters very well (Table S3). We used the projector augmented wave (PAW) method65 as implemented in VASP,66-68 a 400 eV cutoff for the plane-wave expansion, and an electronic convergence criterion of 10−6 eV. Structural optimization was performed until residual forces fell below 5×10−3 eV Å−1. Reciprocal space was sampled on Γcentred grids with densities between 0.02 and 0.05 Å−1. To model the quaternary phases, we employed 2×2×1 supercell expansions (84 atoms) for the trigonal structure in a hexagonal cell setup, and 2×2×2 supercells (56 atoms) for their metastable rocksalt-type counterparts. Mixed occupancies in the intermediate phases were modeled using the special quasi-random structure approach as implemented in the Alloy Theoretic Automated Toolkit (ATAT).69, 70

ASSOCIATED CONTENT Supporting Information. Crystal data, refinement results, and powder-diffraction patterns for SnSb2Te4, SnSbBiTe4, and SnBi2Te4 in trigonal and cubic phase. Stoichiometric composition of samples as determined by EDX. Computational results for unit cell volumes, bond angles and density of states. Visualization of carrier mobility and concentration during annealing of SnSb2Te4 and SnBi2Te4. Temperature-dependent resistivity, Hall and Seebeck data of thin-film samples. Estimation of band gap from Seebeck measurements. Power factor for thin-film samples. Supplementary discussion regarding excess vacancies in GeSb2Te4 and GeTe.

AUTHOR INFORMATION Corresponding Authors * E-mail: [email protected] (R.D.); [email protected] (M.W.).

Present Addresses # (V.L.D.) Engineering Laboratory, University of Cambridge, Trumpington Street, Cambridge CB2 1 PZ, United Kingdom

Author Contributions T.S., J.D.H., T.L., and P.M. performed experiments; P.M.K. performed DFT computations. All authors contributed to data analysis and discussions. T.S., P.M.K., V.L.D., and M.W. wrote the manuscript with input from all authors. M.M. and M.W. initialized, and R.D. and M.W. supervised the project.

Notes The authors declare no competing financial interest.

ACKNOWLEDGMENT

Page 8 of 11

This study has been supported by Deutsche Forschungsgemeinschaft (DFG) in the framework of SFB 917 “Nanoswitches” as well as Mo/058/13-2. The research leading to the experimental results has received funding from the European Union Seventh Framework Programme (FP7/2007-2013) under grant agreement no. 340698. Computer time was provided by the Jülich–Aachen Research Alliance (JARA-HPC project 0033).

REFERENCES 1. Wuttig, M.; Yamada, N., Phase-Change Materials for Rewriteable Data Storage. Nat. Mater. 2007, 6, 824–832. 2. Raoux, S.; Wełnic, W.; Ielmini, D., Phase Change Materials and Their Application to Nonvolatile Memories. Chem. Rev. 2010, 110, 240–267. 3. Schneider, M. N.; Rosenthal, T.; Stiewe, C.; Oeckler, O., From phase-change materials to thermoelectrics? Z. Kristallogr. 2010, 225, 463–470. 4. Hasan, M. Z.; Kane, C. L., Colloquium: Topological insulators. Rev. Mod. Phys. 2010, 82, 3045–3067. 5. Qi, X.-L.; Zhang, S.-C., Topological insulators and superconductors. Rev. Mod. Phys. 2011, 83, 1057–1110. 6. Ando, Y., Topological Insulator Materials. J. Phys. Soc. Jpn. 2013, 82, 102001. 7. Matsunaga, T.; Yamada, N.; Kojima, R.; Shamoto, S.; Sato, M.; Tanida, H.; Uruga, T.; Kohara, S.; Takata, M.; Zalden, P., et al., Phase-Change Materials: Vibrational Softening upon Crystallization and Its Impact on Thermal Properties. Adv. Funct. Mater. 2011, 21, 2232–2239. 8. Loke, D.; Lee, T. H.; Wang, W. J.; Shi, L. P.; Zhao, R.; Yeo, Y. C.; Chong, T. C.; Elliott, S. R., Breaking the Speed Limits of PhaseChange Memory. Science 2012, 336, 1566–1569. 9. Deringer, V. L.; Dronskowski, R.; Wuttig, M., Microscopic Complexity in Phase-Change Materials and its Role for Applications. Adv. Funct. Mater. 2015, 25, 6343–6359. 10. Matsunaga, T.; Kojima, R.; Yamada, N.; Kifune, K.; Kubota, Y.; Tabata, Y.; Takata, M., Single Structure Widely Distributed in a GeTe-Sb2Te3 Pseudobinary System: A Rock Salt Structure is Retained by Intrinsically Containing an Enormous Number of Vacancies within its Crystal. Inorg. Chem. 2006, 45, 2235–2241. 11. Wuttig, M.; Lüsebrink, D.; Wamwangi, D.; Wełnic, W.; Gilleßen, M.; Dronskowski, R., The role of vacancies and local distortions in the design of new phase-change materials. Nat. Mater. 2007, 6, 122–128. 12. Edwards, A. H.; Pineda, A. C.; Schultz, P. A.; Martin, M. G.; Thompson, A. P.; Hjalmarson, H. P.; Umrigar, C. J., Electronic structure of intrinsic defects in crystalline germanium telluride. Phys. Rev. B 2006, 73, 045210. 13. Caravati, S.; Bernasconi, M.; Kühne, T. D.; Krack, M.; Parrinello, M., First principles study of crystalline and amorphous Ge2Sb2Te5 and the effects of stoichiometric defects. J. Phys.: Condens. Matter 2009, 21, 499803. 14. Zhang, W.; Thiess, A.; Zalden, P.; Zeller, R.; Dederichs, P. H.; Raty, J.-Y.; Wuttig, M.; Blügel, S.; Mazzarello, R., Role of vacancies in metal–insulator transitions of crystalline phase– change materials. Nat. Mater. 2012, 11, 952–956. 15. Deringer, V. L.; Lumeij, M.; Stoffel, R. P.; Dronskowski, R., Mechanisms of Atomic Motion Through Crystalline GeTe. Chem. Mater. 2013, 25, 2220–2226. 16. Xu, M.; Zhang, W.; Mazzarello, R.; Wuttig, M., Disorder Control in Crystalline GeSb2Te4 Using High Pressure. Adv. Sci. 2015, 2, 1500117. 17. Siegrist, T.; Jost, P.; Volker, H.; Woda, M.; Merkelbach, P.; Schlockermann, C.; Wuttig, M., Disorder–induced localization in crystalline phase–change materials. Nat. Mater. 2011, 10, 202– 208.

ACS Paragon Plus Environment

Page 9 of 11

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chemistry of Materials

18. Schneider, M. N.; Biquard, X.; Stiewe, C.; Schröder, T.; Urban, P.; Oeckler, O., From metastable to stable modifications— in situ Laue diffraction investigation of diffusion processes during the phase transitions of (GeTe)nSb2Te3 (6 < n < 15) crystals. Chem. Commun. 2012, 48, 2192–2194. 19. Slack, G. A., In CRC Handbook of Thermoelectrics, Rowe, D. M., Ed. CRC Press, Boca Raton, FL: 1995; pp 407–408. 20. Siegert, K. S.; Lange, F. R. L.; Sittner, E. R.; Volker, H.; Schlockermann, C.; Siegrist, T.; Wuttig, M., Impact of vacancy ordering on thermal transport in crystalline phase-change materials. Rep. Prog. Phys. 2015, 78, 013001. 21. Kumar, M.; Vora-ud, A.; Seetawan, T.; Han, J. G., Enhancement in Thermoelectric Properties of Cubic Ge2Sb2Te5 Thin Films by Introducing Structural Disorder. Energy Technol. 2016, 4, 375–379. 22. Fu, L.; Kane, C. L., Topological insulators with inversion symmetry. Phys. Rev. B 2007, 76, 045302. 23. Hsieh, D.; Qian, D.; Wray, L.; Xia, Y.; Hor, Y. S.; Cava, R. J.; Hasan, M. Z., A topological Dirac insulator in a quantum spin Hall phase. Nature 2008, 452, 970–974. 24. Chen, Y. L.; Analytis, J. G.; Chu, J. H.; Liu, Z. K.; Mo, S. K.; Qi, X. L.; Zhang, H. J.; Lu, D. H.; Dai, X.; Fang, Z., et al., Experimental Realization of a Three-Dimensional Topological Insulator, Bi2Te3. Science 2009, 325, 178–181. 25. Zhang, H.; Liu, C.-X.; Qi, X.-L.; Dai, X.; Fang, Z.; Zhang, S.C., Topological insulators in Bi2Se3, Bi2Te3 and Sb2Te3 with a single Dirac cone on the surface. Nat. Phys. 2009, 5, 438–442. 26. Zhang, J.; Chang, C.-Z.; Zhang, Z.; Wen, J.; Feng, X.; Li, K.; Liu, M.; He, K.; Wang, L.; Chen, X., et al., Band structure engineering in (Bi1−xSbx)2Te3 ternary topological insulators. Nat. Commun. 2011, 2, 574. 27. Pauly, C.; Liebmann, M.; Giussani, A.; Kellner, J.; Just, S.; Sánchez-Barriga, J.; Rienks, E.; Rader, O.; Calarco, R.; Bihlmayer, G., et al., Evidence for topological band inversion of the phase change material Ge2Sb2Te5. Appl. Phys. Lett. 2013, 103, 243109. 28. Vergniory, M. G.; Menshchikova, T. V.; Silkin, I. V.; Koroteev, Y. M.; Eremeev, S. V.; Chulkov, E. V., Electronic and spin structure of a family of Sn-based ternary topological insulators. Phys. Rev. B 2015, 92, 045134. 29. Menshchikova, T. V.; Eremeev, S. V.; Koroteev, Y. M.; Kuznetsov, V. M.; Chulkov, E. V., Ternary compounds based on binary topological insulators as an efficient way for modifying the Dirac cone. JETP Lett. 2011, 93, 15–20. 30. Eremeev, S. V.; Menshchikova, T. V.; Silkin, I. V.; Vergniory, M. G.; Echenique, P. M.; Chulkov, E. V., Sublattice effect on topological surface states in complex (SnTe)n>1(Bi2Te3)m=1 compounds. Phys. Rev. B 2015, 91, 245145. 31. Menshchikova, T. V.; Eremeev, S. V.; Chulkov, E. V., Electronic structure of SnSb2Te4 and PbSb2Te4 topological insulators. Appl. Surf. Sci. 2013, 267, 1–3. 32. Matsunaga, T.; Yamada, N., Structural investigation of GeSb2Te4: A high-speed phase-change material. Phys. Rev. B 2004, 69, 104111. 33. Talybov, A. G., Electron diffraction study of the structure of SnSb2Te4. Sov. Phys. Crystallogr. 1961, 6, 40–44. 34. Agaev, K. A.; Talybov, A. G., Electron-Diffraction Analysis of the Structure of GeSb2Te4. Sov. Phys. Crystallogr. 1966, 11, 400– 402. 35. Zhukova, T. B.; Zaslavskii, A. I., Crystal structures of the compounds PbBi4Te7, PbBi2Te4, SnBi4Te7, SnBi2Te4, SnSb2Te4, and GeBi4Te7. Sov. Phys. Crystallogr. 1972, 16, 796–800. 36. Mao, Z. L.; Chen, H.; Jung, A. l., The structure and crystallization characteristics of phase change optical disk material Ge1Sb2Te4. J. Appl. Phys. 1995, 78, 2338–2342. 37. Karpinskii, O. G.; Shelimova, L. E.; Kretova, M. A.; Fleurial, J. P., An X-ray study of the mixed-layered compounds of (GeTe)n(Sb2Te3)m homologous series. J. Alloys Compd. 1998, 268, 112–117.

38. Shelimova, L. E.; Karpinskii, O. G.; Kretova, M. A.; Kosyakov, V. I.; Shestakov, V. A.; Zemskov, V. S.; Kuznetsov, F. A., Homologous series of layered tetradymite-like compounds in the Sb-Te and GeTe-Sb2Te3 systems. Inorg. Mater. 2000, 36, 768–775. 39. Oeckler, O.; Schneider, M. N.; Fahrnbauer, F.; Vaughan, G., Atom distribution in SnSb2Te4 by resonant X-ray diffraction. Solid State Sci. 2011, 13, 1157–1161. 40. Volker, H. Disorder and electrical transport in phase change materials. Dissertation, RWTH Aachen University, 2013. 41. Kuznetsova, L. A.; Kuznetsov, V. L.; Rowe, D. M., Thermoelectric properties and crystal structure of ternary compounds in the Ge(Sn,Pb)Te–Bi2Te3 systems. J. Phys. Chem. Solids 2000, 61, 1269–1274. 42. Karpinskii, O. G.; Shelimova, L. E.; Kretova, M. A.; Avilov, E. S.; Zemskov, V. S., X-ray Diffraction Study of Mixed-Layer Compounds in the Pseudobinary System SnTe–Bi2Te3. Inorg. Mater. 2003, 39, 240–246. 43. A numerical prediction of defect concentrations, including anti-site defects and vacancies would be of interest but requires much further computational effort, including at minimum the treatment of charged defects, and likely that of defect complexes. Such endeavors are currently underway but outside the scope of the present work. 44. Da Silva, J. L. F.; Walsh, A.; Lee, H., Insights into the structure of the stable and metastable (GeTe)m(Sb2Te3)n compounds. Phys. Rev. B 2008, 78, 224111. 45. Jost, P.; Volker, H.; Poitz, A.; Poltorak, C.; Zalden, P.; Schäfer, T.; Lange, F. R. L.; Schmidt, R. M.; Hollander, B.; Wirtssohn, M. R., et al., Disorder-Induced Localization in Crystalline Pseudo-Binary GeTe-Sb2Te3 Alloys between Ge3Sb2Te6 and GeTe. Adv. Funct. Mater. 2015, 25, 6399–6406. 46. Volker, H.; Jost, P.; Wuttig, M., Low-Temperature Transport in Crystalline Ge1Sb2Te4. Adv. Funct. Mater. 2015, 25, 6390–6398. 47. We remark that an elaborate analysis would require information about the low-temperature resistivity and extrapolation to 0 K, where the definition of metal and insulator is meaningful. However, here the TCR may serve as an ample criterion for the relevance of Anderson localization. 48. Seeger, K., Semiconductor Physics – An Introduction. 9 ed.; Springer-Verlag Berlin Heidelberg: 2004. 49. Sehr, R.; Testardi, L. R., The optical properties of p-type Bi2Te3–Sb2Te3 alloys between 2–15 microns. J. Phys. Chem. Solids 1962, 23, 1219–1224. 50. Reimann, J.; Güdde, J.; Kuroda, K.; Chulkov, E. V.; Höfer, U., Spectroscopy and dynamics of unoccupied electronic states of the topological insulators Sb2Te3 and Sb2Te2S. Phys. Rev. B 2014, 90, 081106. 51. Michiardi, M.; Aguilera, I.; Bianchi, M.; de Carvalho, V. E.; Ladeira, L. O.; Teixeira, N. G.; Soares, E. A.; Friedrich, C.; Blügel, S.; Hofmann, P., Bulk band structure of Bi2Te3. Phys. Rev. B 2014, 90, 075105. 52. Villagonzalo, C.; Romer, R. A.; Schreiber, M., Thermoelectric transport properties in disordered systems near the Anderson transition. Eur. Phys. J. B 1999, 12, 179–189. 53. Imry, Y.; Amir, A., The Localization Transition at Finite Temperatures: Electric and Thermal Transport. Int. J. Mod. Phys. B 2010, 24, 1789–1810. 54. Teo, J. C. Y.; Fu, L.; Kane, C. L., Surface states and topological invariants in three-dimensional topological insulators: Application to Bi1–xSbx. Phys. Rev. B 2008, 78, 045426. 55. Lenoir, B.; Cassart, M.; Michenaud, J. P.; Scherrer, H.; Scherrer, S., Transport properties of Bi-RICH Bi-Sb alloys. J. Phys. Chem. Solids 1996, 57, 89–99. 56. Weyrich, C.; Drögeler, M.; Kampmeier, J.; Eschbach, M.; Mussler, G.; Merzenich, T.; Stoica, T.; Batov, I. E.; Schubert, J.; Plucinski, L., et al., Growth, characterization, and transport properties of ternary (Bi1−xSbx)2Te3 topological insulator layers. J. Phys.: Condens. Matter 2016, 28, 495501.

ACS Paragon Plus Environment

Chemistry of Materials

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

57. Scanlon, D. O.; King, P. D. C.; Singh, R. P.; de la Torre, A.; McKeown Walker, S.; Balakrishnan, G.; Baumberger, F.; Catlow, C. R. A., Controlling Bulk Conductivity in Topological Insulators: Key Role of Anti-Site Defects. Adv. Mater. 2012, 24, 2154–2158. 58. Rodríguez-Carvajal, J., Recent developments of the program FULLPROF. Commission on powder diffraction (IUCr). Newsletter 2001, 26, 12–19. 59. van der Pauw, L. J., A method of measuring specific resistivity and Hall effect of discs of arbitrary shape. Philips Res. Rep. 1958, 13, 1–9. 60. van der Pauw, L. J., A Method of Measuring the Resistivity and Hall Coefficient on Lamellae and Arbitrary Shape. Philips Techn. Rev. 1958, 20, 220–224. 61. Russell, B. R.; Wahlig, C., A New Method for the Measurement of Hall Coefficients. Rev. Sci. Instrum. 1950, 21, 1028–1029. 62. Rzewuski, H.; Werner, Z., New Double‐Frequency Method for Hall Coefficient Measurements. Rev. Sci. Instrum. 1965, 36, 235–236. 63. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H., A Consistent and Accurate ab initio Parametrization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. J. Chem. Phys. 2010, 132, 154104.

Page 10 of 11

64. Perdew, J. P.; Burke, K.; Ernzerhof, M., Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865-3868. 65. Blöchl, P. E., Projector Augmented-Wave Method. Phys. Rev. B 1994, 50, 17953–17979. 66. Kresse, G.; Furthmüller, J., Efficient Iterative Schemes for Ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B 1996, 54, 11169–11186. 67. Kresse, G.; Furthmüller, J., Efficiency of ab-initio Total Energy Calculations for Metals and Semiconductors Using a PlaneWave Basis Set. Comput. Mater. Sci. 1996, 6, 15–50. 68. Kresse, G.; Joubert, D., From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 1999, 59, 1758-1775. 69. van de Walle, A.; Asta, M.; Ceder, G., The Alloy Theoretic Automated Toolkit: A user guide. CALPHAD: Comput. Coupling Phase Diagrams Thermochem. 2002, 26, 539-553. 70. van de Walle, A.; Tiwary, P.; de Jong, M.; Olmsted, D. L.; Asta, M.; Dick, A.; Shin, D.; Wang, Y.; Chen, L.-Q.; Liu, Z.-K., Efficient stochastic generation of special quasirandom structures. CALPHAD: Comput. Coupling Phase Diagrams Thermochem. 2013, 42, 13-18.

ACS Paragon Plus Environment

Page 11 of 11

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chemistry of Materials

Table of Contents artwork

ACS Paragon Plus Environment

11