Pressure-Driven Enhancement of Topological Insulating State in Tin

Mar 28, 2013 - College of Materials Science and Engineering & State Key Laboratory ... Department of Physics and HiPSEC, University of Nevada, Las Veg...
4 downloads 0 Views 2MB Size
Download PDF
Article pubs.acs.org/JPCC

Pressure-Driven Enhancement of Topological Insulating State in Tin Telluride Dan Zhou,†,‡ Quan Li,*,†,‡ Yanming Ma,† Qiliang Cui,† and Changfeng Chen*,‡ †

College of Materials Science and Engineering & State Key Laboratory of Superhard Materials, Jilin University, Changchun 130012, China ‡ Department of Physics and HiPSEC, University of Nevada, Las Vegas, Nevada 89154, United States

ABSTRACT: Recent discovery of a distinct type of topological insulating state in tin telluride (SnTe) has reignited great interest in this narrow bandgap semiconductor. A pressing task is to understand the stability of this intriguing state of matter and the response of its electronic properties under the influence of external conditions that may alter the underlying fundamental physics. In this work, we examine by first-principles calculations the effect of pressure on the phonon dispersion, electronic band structure, Fermi surface, and charge distribution of SnTe. We show that pressure suppresses a soft optical phonon mode and enhances the stability of the cubic (B1, Fm3̅m) phase of SnTe, which is a key requirement for the observed topological insulating state. Pressure also drives an electronic topological transition that alters the shape and connectivity of the Fermi surface, and it also promotes a charge redistribution that results in an enhanced bonding interaction that strengthens the cubic crystalline structure. Our calculations further demonstrate that the electronic band gap of the cubic phase, another key parameter of the topological insulating state in SnTe, increases in size monotonically with increasing pressure. These results indicate that pressure makes the topological insulating state in the cubic phase of SnTe more stable and robust; pressure is also effective in tuning the electronic band structure, which is expected to have significant impact on a wide range of physical properties crucial to potential applications. insulator phase in SnTe.12 Most recently, the topological crystalline insulator phase and topological phase transition have also been observed in Pb1‑xSnxTe by ARPES measurements.13 Following these exciting discoveries, a pressing task is to explore the behavior of the recently identified topological insulating state and characterize its response to external physical conditions that may alter the underlying fundamental physics. Of particular interest are the possible changes in the stability of its cubic crystalline phase and the nature and size of its electronic band gap, which are essential ingredients for maintaining the topological insulating state in SnTe.10,11 Tin telluride is a direct narrow band gap semiconductor with a gap of 0.18 eV (0.30 eV) at room (4.2 K) temperature.14 It has been extensively studied over the past several decades, and

1. INTRODUCTION Topological insulators are intriguing states of quantum matter characterized by an insulating gap in the bulk and conducting gapless edges or surface states in the boundaries.1−3 Their discovery has generated great interest because of the scientific importance of the observed phenomena and promising potential of these materials for high-temperature spintronics applications.4−9 Search for additional members of this class of materials has been continuing. Recent theoretical work identified tin telluride (SnTe) as a distinct type of topological crystalline insulator, in which the metallic surface states are protected by the mirror symmetry of the crystal, in contrast to the time-reversal symmetry protection in the earlier identified Z2 topological insulators.10,11 Subsequently, angle-resolved photoemission spectra (ARPES) detected surface electronic states consisting of four Dirac cones (even number of so-called band-inversion points) in the first surface Brillouin zone, providing experimental evidence for the topological crystalline © 2013 American Chemical Society

Received: February 25, 2013 Revised: March 26, 2013 Published: March 28, 2013 8437

dx.doi.org/10.1021/jp401928j | J. Phys. Chem. C 2013, 117, 8437−8442

The Journal of Physical Chemistry C

Article

Perdew−Burke−Ernzerh (PBE) generalized gradient approximation (GGA) functional28 and the Heyd−Scuseria−Ernzerhof (HSE) hybrid functional29,30 with spin−orbit interaction as the exchange-correlation potential were applied to investigate the electronic structures. The projector augmented wave (PAW) potentials31 with 5s25p2 and 5s25p4 electrons as valence were adopted for Sn and Te, respectively. An energy cutoff of 300 eV and 16 × 16 × 16 Monkhorst-Pack k meshes were chosen to ensure that enthalpy calculations were well converged to better than 1 meV/atom. The phonon frequencies were calculated using a direct supercell method (4 × 4 × 4, 128 atoms supercell), which uses the forces obtained by the Hellmann−Feynman theorem.32 Bader charge analysis33 was used to analyze the charge state. Note that the HSE calculation is only used to correct the band energies in assessment of the band gap, but it is not used in the calculation of electronic density of states, Fermi surfaces, phonon, and Bader charge analysis, which are less sensitive to the choice of the functionals.

it has remained a topic of great interest because of its fundamental physics and potential applications in electronic devices. An unusual feature of the electronic structure of SnTe is that the ordering of its conduction and valence band near the Fermi energy is inverted relative to a normal semiconductor like PbTe whose electronic band structure smoothly connects to the atomic limit.10 The so-called “negative band gap” (or band inversion) in SnTe occurs between a valence band maximum near the L points in the Brillouin zone with L6− symmetry and the conduction band minimum with L6+ symmetry, which are opposite to those in PbTe.15−22 A transition from the normal band structure in PbTe to that of SnTe with the nontrivial inverted band gap has been illustrated by examining the band gap evolution in Pb1‑xSnxTe as a function of alloy composition,23 which shows that the band gap of the alloy initially decreases with increasing value of x, leading to its full closure, and then reopens and increases in the inverted direction with further rising x. Moreover, ARPES and ab initio calculations revealed complex Fermi surface structure near the L points, showing topological changes in the constantenergy surface from disconnected pockets, to open tubes, and then to cuboids as the binding energy (or hole-doping) increases.24 The narrow band gap of SnTe is very sensitive to changes of external conditions such as pressure, doping, or temperature. Given the essential role of the inverted band gap in the topological insulating state of SnTe,10,11 it is crucial that we establish an understanding of its behavior under changing external conditions. While the response of the electronic properties of SnTe to changing temperature and hole-doping has been extensively studied,10,14,24−26 the effect of pressure on the structural stability of the cubic phase of SnTe and its electronic properties remains largely unexplored. In this article, we present a systematic first-principles study of the influence of applied pressure on the structural and electronic properties of the topological insulating state that exists in the cubic phase of SnTe. We have calculated the phonon dispersion curve to examine the dynamic stability of the cubic SnTe structure under pressure, and our results indicate that pressure suppresses a soft optical phonon associated with a structural instability, thus strengthening the cubic phase. Our electronic band structure calculations reveal a significant pressure effect on the topological structure of the Fermi surface and its impact on the charge redistribution, which, in turn, strengthens the bonding and structural stability of the cubic phase. The band structure calculations also show that the electronic band gap of SnTe increases in size considerably with applied pressure, making the topological insulating state more robust under pressure. These results show that pressure stabilizes and enhances the topological crystalline insulator state of SnTe and effectively tunes its electronic structure and Fermi surface, which are expected to have significant impact on the transport and optical properties crucial for its applications. In the following sections, we will first provide some details of our computational methods and procedure, and then present the results and discussion on the calculated structural and electronic properties of SnTe under changing pressure. We summarize our main conclusions at the end of the article.

3. RESULTS AND DISCUSSION With the consideration of spin−orbit coupling, the theoretical equilibrium lattice parameters calculated by the standard PBEGGA functional and HSE hybrid functional for SnTe are 6.404 Å and 6.356 Å, which are 1.3% and 0.5% larger than the experimental data of 6.327 Å,34 respectively. To investigate the effect of pressure on the dynamic structural stability of SnTe, we have calculated phonon dispersion relationships at various pressure points. From our recent study of structural evolution of SnTe at high pressure using an integrated approach of angle dispersive synchrotron X-ray diffraction combined with a firstprinciples structural search method, SnTe undergoes pressuredriven phase transitions from the ambient-pressure Fm3̅m structure to three coexisting intermediate phases of Pnma, Cmcm, and a GeS type structure at 5.0 GPa.35 Therefore, the upper limit of pressure in this study was chosen as 5.0 GPa. The obtained results shown in Figure 1 indicate that the cubic (Fm3̅m) structure is unstable with imaginary transverse optical phonon frequencies near the zone center (Γ point). This result is in agreement with the previously reported Raman scattering experiment.25 This structural instability can be understood in terms of the strong coupling between the interband electronic excitation and the transverse optical phonon.36 Because of the softening phonon mode, the Fm3m ̅ structure transforms to a rhombohedral structure with a space group of R3m driven by a small dimerization in the unit cell. Our results show that rising pressure leads to monotonically increasing frequencies of the soft optical phonon mode (Figure 1a). The squared phonon frequency of the optical branch at the Γ point exhibits an almost linearly increasing relationship with pressure (Figure 1b), and the cubic phase is dynamically stable above 0.5 GPa. The current lattice dynamics results demonstrate that rising pressure is able to stabilize the cubic phase of SnTe, which is consistent with the experimental results.25,26 We now turn to the electronic band structure and examine the influence of pressure, especially its effect on the band gap, which is a crucial physical parameter for maintaining the topological insulating state in SnTe and for its applications. We have performed the band structure calculations using the standard PBE-GGA and the HSE hybrid functional. In heavy atoms, relativistic shifts of the electron energy levels lead to a strong spin−orbit coupling. The band structures of SnTe and PbTe have been previously examined by electronic band

2. COMPUTATIONAL DETAILS We have carried out first-principles energetic and electronic band structure calculations using the density functional theory as implemented in the VASP code.27 Both the standard 8438

dx.doi.org/10.1021/jp401928j | J. Phys. Chem. C 2013, 117, 8437−8442

The Journal of Physical Chemistry C

Article

Figure 1. (a) Calculated phonon dispersion relationship of the cubic structure of SnTe at different pressures. The inset shows the first Brillouin zone with selected high-symmetry points indicated. (b) The calculated squared phonon frequency ω2 of the optical branch at the Γ point as a function of pressure.

Figure 2. Band structures of SnTe at (a) 0 GPa with the lattice parameter determined by using the PBE or HSE functionals, (b) no applied pressure, but using the experimental lattice parameter of a = 6.327 Å, and (c) at 5 GPa using the lattice parameter by using the PBE or HSE functionals. (d) The calculated band gap of SnTe as a function of pressure, where results at pressures corresponding to lattice parameters below the experimental value are represented by the dashed lines.

structure calculations with and without the consideration of spin−orbit coupling.15,22,24,37,38 These calculations identify significant effects of the spin−orbit coupling on the band structures for SnTe and PbTe, e.g., splitting in the degenerate bands along high-symmetry directions and modifying the band gaps.15,22,24,37,38 It is therefore necessary to consider the spin− orbit coupling explicitly to describe the electronic properties of SnTe. The calculated results are shown in Figure 2. The results of our calculations using the standard PBE-GGA functional without applied pressure (Figure 2a) are in good agreement with recently reported data.10,23 In SnTe, the L6+ state is mainly derived from Te and the L6− state from Sn. The occupation in the conduction band minimum (L6+ state) and valence band maximum (L6− state) in SnTe is switched (i.e., inverted) near the L points relative to those in PbTe, and this peculiar behavior stems from the different relativistic shifts of the valence states between Pb and Sn, e.g., 2.75 eV for the s state and 0.73 eV for the p state.16 The PBE-GGA calculations yield a small direct band gap of 0.112 eV at the L points with the optimized lattice parameters, indicating its narrow band gap semiconductor nature. This band gap by PEB-GGA is smaller than the experimental value of 0.18 eV measured at room temperature.14 It should be mentioned that the experimental band gap of 0.30 eV14 at the low temperature of 4.2 K corresponds to the rhombohedral phase and, therefore, cannot be directly compared to the calculated result obtained for the cubic phase. As the GGA usually underestimates electronic band gap, we have performed the calculations using a hybrid functional where 25% of the GGA exchange potential is replaced by screened Fock exchange. This type of HSE hybrid functional29,30 usually provides a better description of the electronic band structure, especially the band gap.39−41 Intriguingly, our calculations show that the band gap obtained using the HSE hybrid functional is reduced to 0.004 eV, which is smaller than that using the standard PBE-GGA. This downward shift of the HSE band gap is attributed to the

unusual inverted band gap of SnTe. In a normal semiconductor, HSE tends to move the bottom of the conduction band upward while moving the top of the valence band downward, thus further opening up the band gap and correcting the underestimation by the GGA calculations. In SnTe, the bottom of the conduction band and top of the valence band states are switched; the HSE calculations still have the same effect on these states, resulting in the lowering of the bottom of the conduction band and rising of the valence band. This reversed process leads to the significantly reduced band gap as obtained in our calculations. It should be noted, however, that the theoretical lattice parameters obtained in both the GGA and HSE calculations are larger than the experimental value obtained at room temperature. Using the experimental lattice parameter of 6.327 Å, our calculations produced the electronic band gap of 0.173 and 0.094 eV for SnTe using the PBE-GGA and HSE functional, respectively. While the GGA band gap is in better agreement with the experimental value of 0.18 eV, one should exercise caution in drawing conclusions about the accuracy of the GGA and HSE functionals applied to SnTe since the reported calculations do not explicitly consider the temperature effect (i.e., thermal vibration of atomic positions) that may have a significant influence on the band gap.42 It is important that both functionals produce the same trend for the response of the band gap to pressure as we discuss below. Our calculated results presented in Figure 2 show that the inverted band gap of SnTe is sensitive to the change of external pressure. As the pressure increases, both the GGA and HSE calculations predict that the size of the inverted band gap will have a positive pressure dependence. This indicates that the topological crystalline insulator state in SnTe will be enhanced by pressure. Note also that the results using the HSE and PBE8439

dx.doi.org/10.1021/jp401928j | J. Phys. Chem. C 2013, 117, 8437−8442

The Journal of Physical Chemistry C

Article

GGA are in fair agreement above the pressures corresponding to the experimental lattice parameter. With increasing hydrostatic pressure (decreasing volume), the band structures, especially those near the Fermi energy, expand due to the expansion of the bandwidth by pressure. Because of the symmetry requirements, a crossing of the inverted conduction minimum L6+c and valence maximum L6−v is forbidden, and consequently, they will more strongly repel each other under pressure. The repulsion between the states of the same symmetry, i.e., the conduction minimum L6+c and the valence maximum L6+v will increase with increasing pressure, while the conduction band states L6+c and L6−c have different symmetries and can cross at high enough pressure. Our calculations show that the L6+c states at the L points are significantly pushed up under pressure (results at 5 GPa are shown in Figure 2c). Because of the effect of increased repulsion under high pressure, the shape of the conduction band minimum near the L points changes to that of an inverted saddle, and the conduction band maximum is moved slightly away from the L points. Recent ARPES experiments and density functional theory calculations revealed complex Fermi surface evolution in SnTe with increasing hole doping.24 A key finding is an electronic topological transition, also known as Lifshitz transition, as the binding energy increases, which corresponds to a rising level of hole doping. This result has important implications for understanding the transport and optical properties of SnTe, and there is also empirical evidence26 that doping tends to stabilize the cubic phase of SnTe. Here, we examine the evolution of the electronic structure of SnTe near the Fermi energy under pressure. In Figure 3, we plot constant energy surfaces of the cubic (Fm3̅m) SnTe at different energies by shifting the chemical potential downward, which correspond to Fermi surfaces at different hole (p-)doped concentrations, and compare the results taken at 0 and 5 GPa. With increasing downward energy shift (or higher doping concentration), the

calculated energy surfaces evolve from pockets (at Ef −0.3 eV) to interconnected quasicubic tubes (at Ef −0.5 eV) with cube corners at the L points of the first Brillouin Zone, signifying the occurrence a doping-induced electronic topological transition at zero pressure. This result is in agreement with the previous experimental photoemission spectroscopy and ab initio studies24 on the band structure of SnTe. Our calculations reveal that pressure driven evolution of the band structure promotes the electronic topological transition in SnTe. Results in Figure 3 clearly show that at higher pressure (5 GPa) the transition of the constant energy surfaces from disconnected pockets near the L points to interconnected quasicubic tubes occurs at energies closer to the Fermi energy, which corresponds to a lower level of hole doping. This indicates a stronger tendency of the system toward the electronic topological transition under pressure, which accelerates the process at an energy (e.g., Ef −0.3 eV) where doping alone would not have caused the transition. Pressure has an advantage in that it can be applied and released via external control, allowing a reversible modification of the electronic structure and the topology (i.e., connectivity) of the constant-energy surface, which, in turn, can lead to precise and reversible tuning of a wide range of physical properties, such as resistivity, thermoelectric power, Hall coefficient, phonon dispersion, and electron−phonon coupling.43−48 For further insights into the role of pressure in stabilizing the cubic structure of SnTe, we examine the pressure driven charge redistribution in real space associated with the electronic topological transition that occurs in the momentum space. We show in Figure 4 the band decomposed charge density projected onto the (100) plane of the crystal structure. The results indicate a peculiar reverse electron donation from tellurium to tin atoms. We employ the Bader charge analysis, which provides a description of electron transfer, to quantify the amount of charge belonging to each atom at different pressure points. The calculated charge amounts are 6.6310

Figure 3. Evolution of the three-dimensional constant-energy surfaces at different binding energies, corresponding to the Fermi surfaces at different hole(p-) doped levels, at 0 GPa (upper panels) and 5 GPa (lower panels).

Figure 4. Electronic charge density distribution projected onto the (100) plane for the states in an energy window of 0.30 eV below the Fermi energy at 0 GPa (upper panel) and 5 GPa (lower panel). The charge density is measured in electrons per Å3. 8440

dx.doi.org/10.1021/jp401928j | J. Phys. Chem. C 2013, 117, 8437−8442

The Journal of Physical Chemistry C

Article

gap of the cubic phase increases in size with rising pressure. Pressure also induces an electronic topological transition characterized by a change in the connectivity of the constant energy surface in the momentum space, which directly impacts transport and optical properties. These results show that pressure has a strong influence on the fundamental crystal and electronic structure of SnTe, and as a result, it makes the topological insulating state in SnTe more stable and robust. The sensitive pressure tuning of the electronic band gap and Fermi surface topology offers an effective tool to modulate a wide range of physical properties of SnTe that are important to its potential applications.

(3.3690), 6.6231 (3.3769), 6.6156 (3.3844), 6.6087 (3.3913), 6.6021 (3.3979), and 6.5957 (3.4043) electron for tellurium (tin) at 0, 1, 2, 3, 4, and 5 GPa, respectively. This result reveals that the charge states of tellurium and tin atoms change with increasing pressure. The tellurium (tin) atoms lose (gain) electron charge at an average rate of 7.0 × 10−3 e/GPa. It is noted that while SnTe crystallizes naturally in the cubic phase at ambient conditions, it is known to be in a nonstoichiometric state with a deficiency of tin atoms, resulting in an excess of holes in the material.20 Experimental results suggest that the excessive holes in SnTe tend to stabilize its cubic structure.26 From our calculated charge density results, one can see that at ambient conditions the electrons near the valence band edge are mainly contributed from the tellurium atoms (see Figure 4, where valence charge density within 0.3 eV below the Fermi energy is shown). Under pressure, the charges near the valence band maximum form an interconnected quasicubic filament in the momentum space (see Figure 3); correspondingly, in real space, these charges redistribute among the cubic lattice sites, which is expected to impact the bonding interaction and structural stability. Our phonon calculations reveal a persistent hardening of all the phonon branches of the cubic SnTe phase under increasing pressure (see Figure 1a). This clearly indicates a strengthening of the bonding interaction in the cubic structure in response to the pressure driven charge redistribution. Previous studies indicate that a carrier density of about 5 × 1020 cm−3 or 3.2 × 10−2 holes per SnTe formula unit under atmospheric pressure can suppress the transition from cubic to rhombohedral structure at low temperature.26 From our phonon calculations, the pressure to stabilize the cubic (Fm3̅m) phase of SnTe is 0.5 GPa. At this pressure, our charge analysis indicates a loss of about 3.5 × 10−3 e from tellurium per formula unit. This suggests that pressure driven charge density change is very effective in stabilizing the cubic structure of SnTe. It should be noted, however, that the hole doping by tin deficiency is likely distributed locally and randomly in the crystal structure; in contrast, our calculations study charge redistribution that occurs in an ideal crystal lattice. Consequently, it is not meaningful to make a direct comparison of the experimentally determined hole doping by tin deficiency and the theoretically calculated charge transfer, which is a net loss of electrons that can be viewed as an effective hole doping on the tellurium sites. Nevertheless, it is important that the empirical experimental evidence and our calculated results consistently point to the same overall trend of structural (cubic phase) stabilization by charge redistribution in SnTe either introduced by structural defects or driven by applied pressure.

■

AUTHOR INFORMATION

Corresponding Author

*(Q.L.) Tel: +1-702-895-1714. Fax: +1-702-895-0804. E-mail: [email protected]. (C.C.) Tel: +1-702-895-4230. Fax: +1-702-895-0804. E-mail: [email protected]. Notes

The authors declare no competing financial interest.

■

ACKNOWLEDGMENTS This work was supported by the Department of Energy through Cooperative Agreement DE-FC52-06NA26274 at UNLV and by NSFC (No. 11025418, 91022029, 51202084, 11074089, and 51172087) at JLU.

■

REFERENCES

(1) Hasan, M. Z.; Kane, C. L. Colloquium: Topological Insulators. Rev. Mod. Phys. 2010, 82, 3045. (2) Qi, X. L.; Zhang, S. C. Topological Insulators and Superconductors. Rev. Mod. Phys. 2011, 83, 1057. (3) Moore, J. E. The Birth of Topological Insulators. Nature 2010, 464, 194−198. (4) Hasan, M. Z.; Moore, J. E. Three-Dimensional Topological Insulators. Annu. Rev. Condens. Matter Phys. 2011, 2, 55−78. (5) Xia, Y.; Qian, D.; Hsieh, D.; Wray, L.; Pal, A.; Lin, H.; Bansil, A.; Grauer, D.; Hor, Y. S.; Cava, R. J. Observation of A Large-Gap Topological-Insulator Class with A Single Dirac Cone on the Surface. Nat. Phys. 2009, 5, 398−402. (6) Xu, S. Y.; Xia, Y.; Wray, L. A.; Jia, S.; Meier, F.; Osterwalder, J.; Slomski, B.; Bansil, A.; Lin, H.; Cava, R. J. Topological Phase Transition and Texture Inversion in A Tunable Topological Insulator. Science 2011, 332, 560−564. (7) Schnyder, A. P.; Ryu, S.; Furusaki, A.; Ludwig, A. W. W. Classification of Topological Insulators and Superconductors in Three Spatial Dimensions. Phys. Rev. B 2008, 78, 195125. (8) Mong, R. S. K.; Essin, A. M.; Moore, J. E. Antiferromagnetic Topological Insulators. Phys. Rev. B 2010, 81, 245209. (9) Li, R.; Wang, J.; Qi, X. L.; Zhang, S. C. Dynamical Axion Field in Topological Magnetic Insulators. Nat. Phys. 2010, 6, 284−288. (10) Hsieh, T. H.; Lin, H.; Liu, J.; Duan, W.; Bansil, A.; Fu, L. Topological Crystalline Insulators in the SnTe Material Class. Nat. Commun. 2012, 3, 982. (11) Fu, L. Topological Crystalline Insulators. Phys. Rev. Lett. 2011, 106, 106802. (12) Tanaka, Y.; Ren, Z.; Sato, T.; Nakayama, K.; Souma, S.; Takahashi, T.; Segawa, K.; Ando, Y. Experimental Realization of A Topological Crystalline Insulator in SnTe. Nat. Phys. 2012, 8, 800− 803. (13) Xu, S. Y.; Liu, C.; Alidoust, N.; Neupane, M.; Qian, D.; Belopolski, I.; Denlinger, J. D.; Wang, Y. J.; Lin, H.; Wray, L. A. Observation of A Topological Crystalline Insulator Phase and Topological Phase Transition in Pb1‑xSnxTe. Nat. Commun. 2012, 3, 1192.

4. CONCLUSIONS In summary, we have performed first-principles electronic and phonon calculations to examine the phase stability and electronic structure evolution in a topological crystalline insulator SnTe under pressure. Our phonon calculations show that pressure suppresses a soft optical phonon mode and stabilizes the cubic (Fm3̅m) phase of SnTe, which is predicted to be dynamically stable above 0.5 GPa in the absence of any carrier (hole) doping. Increasing pressure promotes a peculiar charge transfer in the form of a backdonation of electron from the tellurium to tin sites, which further strengthens the bonding interaction of the cubic structure as evidenced by the hardening of all the phonon branches under pressure. The electronic band structure of SnTe also responds sensitively to applied pressure. The inverted band 8441

dx.doi.org/10.1021/jp401928j | J. Phys. Chem. C 2013, 117, 8437−8442

The Journal of Physical Chemistry C

Article

(14) Esaki, L.; Stiles, P. J. New Type of Negative Resistance in Barrier Tunneling. Phys. Rev. Lett. 1966, 16, 1108−1111. (15) Rabii, S. Energy-Band Structure and Electronic Properties of SnTe. Phys. Rev. 1969, 182, 821−828. (16) Dimmock, J. O.; Melngailis, I.; Strauss, A. J. Band Structure and Laser Action in PbxSn1‑xTe. Phys. Rev. Lett. 1966, 16, 1193−1196. (17) Burke, J. R., Jr.; Allgaier, R. S.; Houston, B. B., Jr.; Babiskin, J.; Siebenmann, P. G. Shubnikov-de Haas Effect in SnTe. Phys. Rev. Lett. 1965, 14, 360−361. (18) Tung, Y. W.; Cohen, M. L. Relativistic Band Structure and Electronic Properties of SnTe, GeTe, and PbTe. Phys. Rev. 1969, 180, 823−826. (19) Allgaier, R. S.; Houston, B. Weak-Field Magnetoresistance and the Valence-Band Structure of SnTe. Phys. Rev. B 1972, 5, 2186−2197. (20) Melvin, J. S.; Hendry, D. C. Self-Consistent Relativistic Energy Bands for Tin Telluride. J. Phys. C 1979, 12, 3003−3012. (21) Tsang, Y. W.; Cohen, M. L. Calculation of the Temperature Dependence of the Energy Gaps in PbTe and SnTe. Phys. Rev. B 1971, 3, 1254−1261. (22) Wei, S.-H.; Zunger, A. Electronic and Structural Anomalies in Lead Chalcogenides. Phys. Rev. B 1997, 55, 13605−13610. (23) Gao, X.; Daw, M. S. Investigation of Band Inversion in (Pb,Sn) Te Alloys Using ab initio Calculations. Phys. Rev. B 2008, 77, 033103. (24) Littlewood, P. B.; Mihaila, B.; Schulze, R. K.; Safarik, D. J.; Gubernatis, J. E.; Bostwick, A.; Rotenberg, E.; Opeil, C. P.; Durakiewicz, T.; Smith, J. L.; et al. Band Structure of SnTe Studied by Photoemission Spectroscopy. Phys. Rev. Lett. 2010, 105, 086404. (25) Sugai, S.; Murase, K.; Katayama, S.; Takaoka, S.; Nishi, S.; Kawamura, H. Carrier Density Dependence of Soft TO-Phonon in SnTe by Raman Scattering. Solid State Commun. 1977, 24, 407−409. (26) Salje, E. K. H.; Safarik, D. J.; Modic, K. A.; Gubernatis, J. E.; Cooley, J. C.; Taylor, R. D.; Mihaila, B.; Saxena, A.; Lookman, T.; Smith, J. L.; et al. Tin Telluride: A Weakly Co-Elastic Metal. Phys. Rev. B 2010, 82, 184112. (27) 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. (28) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865−3868. (29) Heyd, J.; Scuseria, G. E.; Ernzerhof, M. Hybrid Functionals Based on A Screened Coulomb Potential. J. Chem. Phys. 2003, 118, 8207. (30) Heyd, J.; Scuseria, G. E.; Ernzerhof, M. Erratum: “Hybrid Functionals Based on A Screened Coulomb Potential”. J. Chem. Phys. 2006, 124, 219906. (31) Kresse, G.; Joubert, D. From Ultrasoft Pseudopotentials to the Projector Augmented-Wave Method. Phys. Rev. B 1999, 59, 1758− 1775. (32) Alfè, D. PHON: A Program to Calculate Phonons Using the Small Displacement Method. Comput. Phys. Commun. 2009, 180, 2622−2633. (33) Henkelman, G.; Arnaldsson, A.; Jónsson, H. A Fast and Robust Algorithm for Bader Decomposition of Charge Density. Comput. Mater. Sci. 2006, 36, 354−360. (34) Bis, R. F.; Dixon, J. R. Applicability of Vegard’s Law to the PbxSn1‑xTe Alloy System. J. Appl. Phys. 1969, 40, 1918−1921. (35) Zhou, D.; Li, Q.; Ma, Y. M.; Cui, Q. L.; Chen, C. F. Unraveling Convoluted Structural Transitions in SnTe at High Pressure. J. Phys. Chem. C 2013, 117, 5352−5357. (36) Kawamura, H.; Katayama, S.; Takano, S.; Hotta, S. Dielectric Constant and Soft Mode of Pb1‑xSnxTe. Solid State Commun. 1974, 14, 259−261. (37) Zhang, Y.; Ke, X. Z.; Chen, C. F.; Yang, J.; Kent, P. R. C. Thermodynamic Properties of PbTe, PbSe, and PbS: First-Principles Study. Phys. Rev. B 2009, 80, 024304. (38) Rabe, K. M.; Joannopoulos, J. D. Ab Initio Relativistic Pseudopotential Study of the Zero-Temperature Structural Properties of SnTe and PbTe. Phys. Rev. B 1985, 32, 2302−2314.

(39) Heyd, J.; Peralta, J. E.; Scuseria, G. E.; Martin, R. L. Energy Band Gaps and Lattice Parameters Evaluated with the Heyd-ScuseriaErnzerhof Screened Hybrid Functional. J. Chem. Phys. 2005, 123, 174101. (40) Janesko, B. G.; Henderson, T. M.; Scuseria, G. E. Screened Hybrid Density Functionals for Solid-State Chemistry and Physics. Phys. Chem. Chem. Phys. 2008, 11, 443−454. (41) Heyd, J.; Scuseria, G. E. Efficient Hybrid Density Functional Calculations in Solids: Assessment of the Heyd−Scuseria−Ernzerhof Screened Coulomb Hybrid Functional. J. Chem. Phys. 2004, 121, 1187. (42) Zhang, Y.; Ke, X.; Kent, P. R. C.; Yang, J.; Chen, C. Anomalous Lattice Dynamics Near the Ferroelectric Instability in PbTe. Phys. Rev. Lett. 2011, 107, 175503. (43) Fast, L.; Ahuja, R.; Nordström, L.; Wills, J. M.; Johansson, B.; Eriksson, O. Anomaly in c/a Ratio of Zn under Pressure. Phys. Rev. Lett. 1997, 79, 2301−2303. (44) Struzhkin, V. V.; Timofeev, Y. A.; Hemley, R. J.; Mao, H. K. Superconducting Tc and Electron-Phonon Coupling in Nb to 132 GPa: Magnetic Susceptibility at Megabar Pressures. Phys. Rev. Lett. 1997, 79, 4262−4265. (45) Li, Z.; Tse, J. S. Phonon Anomaly in High-Pressure Zn. Phys. Rev. Lett. 2000, 85, 5130−5133. (46) Polian, A.; Gauthier, M.; Souza, S. M.; Trichês, D. M.; Cardoso de Lima, J.; Grandi, T. A. Two-Dimensional Pressure-Induced Electronic Topological Transition in Bi2Te3. Phys. Rev. B 2011, 83, 113106. (47) Li, Q.; Li, Y.; Cui, T.; Wang, Y.; Zhang, L. J.; Xie, Y.; Niu, Y. L.; Ma, Y. M.; Zou, G. T. The Effects of Pressure on the Electronic, Transport and Dynamical Properties of AuX2 (X= Al, Ga and In). J. Phys.: Condens. Matter 2007, 19, 425224. (48) Yelland, E. A.; Barraclough, J. M.; Wang, W.; Kamenev, K. V.; Huxley, A. D. High-Field Superconductivity at An Electronic Topological Transition in URhGe. Nat. Phys. 2011, 7, 890.

8442

dx.doi.org/10.1021/jp401928j | J. Phys. Chem. C 2013, 117, 8437−8442