Titanium Trisulfide Monolayer as a Potential ... - ACS Publications

Jan 5, 2017 - State Key Laboratory of Digital Manufacturing Equipment and Technology and School of Mechanical Science and Engineering,. Huazhong ...
0 downloads 0 Views 2MB Size
Subscriber access provided by UNIV OF CALIFORNIA SAN DIEGO LIBRARIES

Article

Titanium trisulfide (TiS3) monolayer as a potential thermoelectric material: first-principles based Boltzmann transport study Jie Zhang, Xiao-Lin Liu, Yanwei Wen, Lu Shi, Rong Chen, Huijun Liu, and Bin Shan ACS Appl. Mater. Interfaces, Just Accepted Manuscript • DOI: 10.1021/acsami.6b14134 • Publication Date (Web): 05 Jan 2017 Downloaded from http://pubs.acs.org on January 6, 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.

ACS Applied Materials & Interfaces 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 25

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

ACS Applied Materials & Interfaces

Titanium trisulfide (TiS3) monolayer as a potential thermoelectric material: first-principles based Boltzmann transport study Jie Zhang1, Xiaolin Liu1, Yanwei Wen1, Lu Shi1, Rong Chen2, Huijun Liu3, Bin Shan1,* 1

State Key Laboratory of Material Processing and Die and Mould Technology and

School of Materials Science and Engineering, Huazhong University of Science and Technology, Wuhan 430074, China 2

State Key Laboratory of Digital Manufacturing Equipment and Technology and

School of Mechanical Science and Engineering, Huazhong University of Science and Technology, Wuhan 430074, China 3

Key Laboratory of Artificial Micro- and Nano-Structures of Ministry of Education and School of Physics and Technology, Wuhan University, Wuhan 430072, China

Abstract: Good electronic transport capacity and low lattice thermal conductivity are beneficial for thermoelectric applications. In this study, the potential use as thermoelectric material for the recently synthesized two-dimensional TiS3 monolayer is explored by applying first-principles method combined with Boltzmann transport theory. Our work demonstrates that carrier transport in TiS3 sheet is orientation dependent, caused by the difference of charge density distribution at band edges. Due to a variety of Ti-S bonds with longer lengths, we find TiS3 monolayer show much lower thermal conductivity compared with that of transition metal dichalcogenides such as MoS2. Combined with high power factor along y direction, a considerable n-type ZT value (3.1) can be achieved at moderate carrier concentration, suggesting TiS3 monolayer as a good candidate for thermoelectric applications.

Keywords: transition metal trichalcogenide (TMTC), TiS3 monolayer, carrier mobility, *

thermal

conductivity,

thermoelectric

performance,

first-principles

Author to whom correspondence should be addressed. Electronic mail: [email protected]. 1

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces

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

Page 2 of 25

calculations, Boltzmann transport equation

1. Introduction The global energy crisis and environmental impact of fossil fuels have led to substantial efforts in developing sustainable energy sources. Thermoelectric devices that are capable of converting heat into electricity and vice versa thus have attracted considerable interest from the scientific community. The performance of a thermoelectric material is described by the dimensionless figure of merit:

ZT = S 2σ T / (κ e + κ L ) .

(1)

Where S is the Seebeck coefficient, σ is the electrical conductivity, T is the temperature, and κ e and κ L are electronic and lattice thermal conductivities, respectively. In order to increase the ZT value of a bulk thermoelectric material, one effective strategy is to enhance the power factor (PF= S 2σ ) by modulation doping and band engineering,1,2,3 the other approach is to reduce the lattice thermal conductivity through the formation of nano-precipitates and all-scale hierarchical architectures.4,5,6 In particular, the seminal work of Hicks et al.7,8 and subsequent experiments9,10,11 have shown that low-dimensional structures could maintain high power factor and low thermal conductivity at the same time. As semiconductors with thickness dependent band gap and high carrier mobility, atomic-layer transition metal dichalcogenides (TMDCs) have attracted growing interest for applications in field-effect transistors (FET) and optoelectronic devices. 12 , 13 , 14 , 15 , 16 In addition, these low-dimensional TMDCs have also been explored as potential thermoelectric materials due to the advantage of flexibility, low cost and minimal environmental impact compared with conventional thermoelectric materials such as Bi2Te3 and PbTe. For example, by using lithium intercalation and filtration technique, Wang et al. synthesized the restacked MoS2 thin-film. It is found that the maximum electrical conductivity and Seebeck coefficient could reach 17.5 S/cm and 93.5 µV/K, respectively.17 Yoshida et al. studied the effect of electric field on the thermopower of ultrathin WSe2 single crystal. They figured out that the power 2

ACS Paragon Plus Environment

Page 3 of 25

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

ACS Applied Materials & Interfaces

factor in both the hole and electron side can be optimized and the maximum values are one-order larger than those obtained by changing chemical composition.18 On the theoretical side, with the help of first-principles calculations and semiclassical Boltzmann transport theory, Kumar et al. concluded that the thermoelectric performance of WSe2 is superior to MoSe2 and a peak ZT value of about 0.8 is reached at 1200 K.19 Huang et al. studied the thermoelectric properties of MX2 (M=Mo, W; X=S, Se) monolayers by a ballistic transport approach under linear response regime. They predicted that the ZT of these materials are generally low, among which, p-type MoS2 monolayer shows the highest peak value at room temperature. 20 Although the above works show that TMDCs could be potential thermoelectric materials, the low efficiency is still the main obstacle for their wide application in energy conversion devices. Recently, a new class of 2D layered materials, namely, transition metal trichalcogenide monolayers MX3 (M=Ti, Zr, Hf; X=S, Se, Te), have been successfully prepared.21,22,23,24 Among them, it is found that TiS3 monolayer possesses a direct band gap and much higher carrier mobility,25,26 which are desirable for thermoelectric applications. Moreover, a more complex atomic configuration compared with TMDCs implies that it would own a relatively lower thermal conductivity. For its further applications as thermoelectric material, it is of importance to evaluate the thermoelectric performance of such monolayer. In this work, using first-principles calculations and Boltzmann transport theory, we give a comprehensive study on the electronic and phonon transport properties of TiS3 monolayer. We show that the single-layer structure exhibit a ZT value of 3.1 at room temperature, suggesting that good thermoelectric performance can be achieved in material constituted by light transition metal titanium and sulfur atoms.

2. Methodology Our theoretical calculations are performed within the first-principles plane-wave pseudopotential formulation27,28,29 as implemented in the Vienna Ab-initio Simulation Package (VASP). 30 The generalized gradient approximation (GGA) with the Perdew-Burke-Ernzerhof (PBE)31 exchange-correlation functional is used. The cutoff 3

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces

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

Page 4 of 25

energy for the plane wave basis is set to be 500 eV. In the ionic relaxation and charge density calculation, Monkhorst-Pack k-mesh32 of 13×9×5 and 7×10×1 are employed for bulk and monolayer TiS3, respectively. For the optimization of bulk crystal structure, van der Waals (vdW) interaction is explicitly included. A vacuum thickness larger than 14 Å is adopted along the direction normal to the atomic plane so that the interaction between the neighboring monolayers can be neglected. The convergence threshold is set as 1 meV for total energy and 0.01 eV/Å for atomic force. In order to give accurate electronic transport coefficients, we calculate the band structure of TiS3 system by using the hybrid density functional in the form of Heyd-Scuseria-Ernzerhof (HSE06).33 Within the Boltzmann transport theory in the relaxation time approximation,34 the Seebeck coefficient ( S ) and electrical conductivity ( σ ) are expressed as:

S=

ek B

σ

 ∂f 0 

∫ d ε  − ∂ε Ξ (ε )

ε −µ k BT

,

(2)

 ∂f 0  Ξ ( ε ) ,  ∂ε 

σ = e2 ∫ dε  − where f 0 =

1 e

( ε − µ )/ k BT

+1

(3)

is the equilibrium Fermi function, k B is the Boltzmann

constant, and µ is the chemical potential, which can be obtained by integrating the density of states (DOS) from the desired chemical potential to the Fermi level. In this v v approach, the so-called transport distribution function Ξ(ε ) = ∑ vkv vkvτ kv can be v k

v calculated based on the band structures. Here vkv and τ kv are the group velocity and v relaxation time at state k , respectively. For the calculation of electronic thermal

conductivity, we apply the Wiedemann-Franz law κ e = Lσ T , where L =2.45×10‒8 WΩK‒2 is the Lorenz number. In the phonon calculations, harmonic phonon properties are studied using finite displacement method performed in Phonopy package. 35 The lattice thermal conductivity can be estimated by using iterative self-consistent method for solving the phonon Boltzmann transport equation as implemented in the ShengBTE code.36,37,38 4

ACS Paragon Plus Environment

Page 5 of 25

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

ACS Applied Materials & Interfaces

A 4×5×1 supercell with Γ point and 3×4×1 supercell with 3×3×1 k-mesh are adopted for the calculations of second- and third-order interatomic force constants (IFCs), respectively. In order to get convergent thermal transport parameters, the interactions up to the eighth nearest neighbors are considered for the anharmonic one.

3. Results and discussion To confirm the validity of our method, we first study the structure and electronic properties of bulk TiS3. As can be seen from Figure 1(a), TiS3 crystallizes in a layered structure with the space group p21m. The neighboring atomic planes composed by the 1D chains of triangular TiS3 unit are held together by van der Waals forces. Table S1 (Supporting information) lists the structural parameters for bulk TiS3 obtained by using different computational framework. It is found that DFT-D2

39

and

optB86b-vdW40 functionals give the best results.41 Besides, we further find that the band structures based on the lattice constants calculated by the two methods are almost the same. Unless otherwise specified, we will use optB86b-vdW functional to take into account the vdW interactions in the following sections. Figure 2(a) shows the band structure of bulk TiS3 by using HSE06. One can see that TiS3 is an indirect semiconductor, with the valence band maximum (VBM) and conduction band minimum (CBM) located at Γ and Z points, respectively. The band gap obtained is 1.09 eV, which is in good agreement with the previous theoretical and experimental works.25,42,43 The calculations show that our approach is reliable for this system, which is the precondition to accurately predict the transport coefficients. In the following, we will apply the same method to deal with the TiS3 monolayer. Due to weak vdW force between the atomic planes, TiS3 single layer can be obtained by mechanical cleavage or liquid-phase exfoliation.21,25 Its corresponding crystal structure is shown in Figure 1(b) and 1(c). Within the monolayer, the outer and inner sulfur atoms are inequivalent and are defined as S1 and S2, respectively. The optimized lattice constants are a=5.03 Å and b=3.42 Å, which are larger than its analogue MoS2 single-layer (a=b=3.19 Å). The Ti-S1 bond marked by d1 has a length of 2.50 Å, while the bond length of Ti-S2 cross and along the atomic plane direction 5

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces

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

Page 6 of 25

are different, which are represented by d2=2.47 Å and d3=2.67 Å, respectively. Besides, the two outer S1 atoms are covalently bonded with a length of 2.04 Å. Note that there is only a uniform Mo-S bond (2.41 Å) in MoS2. The more complex atomic configuration of TiS3 monolayer suggest that it would exhibit lower lattice thermal conductivity, as will be discussed in the following. Figure 2(b) shows that TiS3 sheet is a direct semiconductor with both VBM and CBM located at Γ point. The calculated band gap is 1.16 eV, which is slightly larger than its bulk counterpart, indicating that there is a finite quantum confinement effect in this system.44,45 In order to understand the origin of the VBM and CBM states, we also plot the orbital-decomposed band structures in Figure 2(b)-(d). One can see that VBM is derived from the hybridization between the p orbitals of S2 and d orbitals of Ti, whereas the CBM is mainly contributed by the Ti d states. As S1 has almost no contribution to the band edge states, it is reasonable to expect that the carrier type and concentration of the system can be modulated by substitution at the S1 sites without significant changes in the shape of band extrema, which is beneficial for thermoelectric applications. Based on the calculated band structure, we estimate the electronic transport coefficients within the Boltzmann transport theory. In this approach, the electrical conductivity is calculated with respect to carrier relaxation time ( τ ), which is obtained by τ =

µ m* e

, where µ and m* are carrier mobility and effective mass,

respectively. For two-dimensional system, µ

can be determined by using

deformation potential (DP) theory on the basis of effective mass approximation:46,47

µ=

2eh3C , 3k BTm*2 E12

where the elastic modulus C is defined as C =

(4)

∂2 E , the deformation S0 ∂(∂l / l0 )2

potential constant E1 = ∆E / (∆l / l ) , in which ∆E is the variation of band edge position (VBM and CBM) with the lattice dilation of ∆l / l . The energies of band extremum are calculated with respect to a common vacuum level. Other parameters in 6

ACS Paragon Plus Environment

Page 7 of 25

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

ACS Applied Materials & Interfaces

the above formulae are total energy E and area S0 of the system. Table 1 lists the calculated mobility and relaxation time of TiS3 monolayer at room temperature. One can see that hole mobility and relaxation time along the x direction are higher than those along the y direction, while it is reversed for the case of electron. To understand this phenomenon, we plot the charge density of band edges in Figure S1 (Supporting information). It can be clearly seen that the charge density of VBM shows spatial continuity along the x direction, indicating a sizable hole mobility along this direction. However, the charge density of CBM has a continuous distribution along y direction, which means electron transport is favorable in this direction. Moreover, we find that y-direction electron mobility has an order of 104 cm2V−1s−1, which is much higher than that of MoS2 monolayer (in the range of 60 ~ 200 cm2V−1s−1)12,48 and has potential use in thermoelectric applications. In Figure 3, we plot the calculated room temperature electronic transport coefficients along x and y directions as a function of carrier concentration ( n ) from 1011 to 1014 cm−2. We see that both the p-type and n-type Seebeck coefficient along two directions are almost identical to each other in the carrier concentration range considered. Besides, the absolute values of Seebeck coefficient increase with decreasing carrier concentration, which coincides with the relation: S ∝ (1/ n )

2/3 49

.

On the other hand, due to the direction dependence of carrier relaxation time, the electrical conductivity exhibits strong anisotropy. With the increase of carrier concentration, σ shows a reversed behavior with Seebeck coefficient, implying a compromise must be taken to maximize the PF, as indicated in Fig. 3(c). Since the Seebeck coefficient is normally isotropic as mentioned, we can see from the figure that the highest PF appears in the direction where the electrical conductivity is the largest. We next discuss the thermal transport properties of TiS3 monolayer. Figure 4(a) displays the phonon dispersion relation of this structure. It is found that the phonon dispersion is asymmetric along the Γ-X and Γ-Y directions. Based on the formula

7

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces

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

v=

∂w , the acoustic phonon group velocities near Γ point are calculated and listed in ∂q

Table 2. For comparison, we also present corresponding results of MoS2 monolayer. One see that the group velocity is different for that along x and y directions, which means that the thermal transport of TiS3 monolayer will exhibit orientation dependence. Besides, a smaller acoustic phonon velocity compared with that of MoS2 suggest that a lower lattice thermal conductivity in TiS3 monolayer. Figure 4(b) shows the partial phonon DOS of TiS3 monolayer, we find there is a strong hybridization between the Ti and S1 atoms in the frequency region 0 ~ 5 THz. Since the heat is mainly carried by low frequency phonons, we can conclude that heat transport along Ti-S1 bond is the key thermal transfer mechanism in TiS3 monolayer, as will be discussed later. Figure 4(c) plots the calculated lattice thermal conductivity κ L of the TiS3 layer as a function of temperature in the range from 200 to 800 K. For two-dimensional system, the calculated κ L depends on the selection of layer thickness. In our study, we apply the interlayer separation of the bulk TiS3 (8.72 Å). One can find that the lattice thermal conductivity decreases with increasing temperatures, implying the thermal transport is determined by the Umklapp process. It should be noted that the room temperature κ L are 17.9 and 10.2 W/mK along the y and x directions, respectively, which are much lower than that of MoS2 sheet50,51,52 and suggest that TiS3 monolayer could have favorable thermoelectric performance. In order to understand the reason that thermal transport along the Ti-S1 bond dominates in TiS3 monolayer, we have plotted the phonon group velocity and relaxation time of different modes in Figure 5(a) and (b). One can clearly see in the range of 0 ~ 5 THz, the phonon group velocity and relaxation time are much larger than those of high frequency phonons. As a result, the room temperature accumulative thermal conductivity occupies a lot in this frequency region, as shown in Fig. 5(c). It should be noted that the low frequency phonon DOS are mainly contributed by the hybridization between Ti and S1 atoms as mentioned above, which indicates that thermal carry on Ti-S1 bond is the main heat transfer mechanism in TiS3 monolayer. 8

ACS Paragon Plus Environment

Page 8 of 25

Page 9 of 25

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

ACS Applied Materials & Interfaces

Such a behavior is consistent with the fact that the bond length of Ti-S1 is much smaller than that of Ti-S2 along the planar direction. On the other hand, the longer bond for Ti-S1 and Ti-S2 is also the reason for the lower lattice thermal conductivity of TiS3 sheet than its analogue MoS2 monolayer. Figure 6 shows the percentage contribution of different phonon modes to the total lattice thermal conductivity. We can see that along the x direction, three acoustic phonons dominate the thermal transport in TiS3 monolayer. For example, almost 85% of the total thermal conductivity are contributed by acoustic phonons in the temperature range considered. However, it is not the case along the y direction, where the contributions from the optical phonons are quite large. Similar behavior has also been found in phosphorene and can be attributed to the hybridization of low frequency optical phonons with the acoustic phonons (see Fig. 4(a)).53 Note that the thermal conductivity from optical phonons slowly increase with the temperature, which may be due to the fact that more high frequency phonons are stimulated at high temperatures.54 With all the transport coefficients obtained, we can now evaluate the thermoelectric performance of TiS3 monolayer. Figure 7 shows the room temperature ZT values as a function of carrier concentration. For p-type system, the highest ZT of 0.68 can be achieved along the x direction at the concentration of 2.4×1012 cm‒2. However, the optimal ZT is very small along the y direction due to the relatively large lattice thermal conductivity. Moreover, the ZT value of n-type system can reach as high as 3.1 when the carrier concentration is 7.2×1011 cm‒2. Our calculations suggest that TiS3 monolayer can be a promising thermoelectric material by properly adjusting the carrier concentration.

4. Summary We have investigated the electronic, phonon and thermoelectric properties of TiS3 monolayer by first-principles calculations and Boltzmann theory. It is found that the band gap changes little when reducing the dimension of bulk TiS3, which is much different from other two-dimensional materials such as MoS2 and phosphorene. Since the contribution of S1 atom to the band edges is negligible, substitution at S1 site can 9

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces

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

be used to adjust the carrier concentration of the system. We find that the anisotropic electronic transport of TiS3 monolayer is due to the difference of charge density distribution at the band extrema. On the other hand, TiS3 sheet possesses much lower thermal conductivity than its analogue MoS2 single-layer, which is consistent with its relatively complex atomic configuration. At a moderate carrier concentration, TiS3 monolayer possesses superior thermoelectric performance with a maximum room temperature ZT value of 3.1, suggesting that n-type system can be a very promising thermoelectric material.

10

ACS Paragon Plus Environment

Page 10 of 25

Page 11 of 25

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

ACS Applied Materials & Interfaces

Supporting Information Structural parameters of bulk TiS3 obtained by using different computational framework and charge density of VBM and CBM of TiS3 monolayer.

Acknowledgement This work is supported by the National Basic Research Program of China (2013CB934800), National Natural Science Foundation of China (51572097, 51575217), and the Hubei Province Funds for Distinguished Young scientists (2014CFA018 and 2015CFA034). Rong Chen acknowledges the Thousand Young Talents Plan, the Recruitment Program of Global Experts. The calculations are done at the Texas Advanced Computing Center (TACC) at The University of Texas at Austin (http://www.tacc.utexas.edu).

11

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces

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

Page 12 of 25

Table 1. The deformation potential constant E1, elastic modulus C, effective mass m*, carrier mobility µ and relaxation time τ of TiS3 monolayer at room temperature. Direction

carriers

E1 (eV)

C (N/m)

µ (×103

*

m (m0)

2

−1 −1

τ (10−13s)

cm V s ) hole

3.16

84.41

0.30

1.33

2.27

electron

1.08

84.41

1.51

0.45

3.87

hole

−3.99

134.83

1.02

0.12

0.67

electron

0.81

134.83

0.38

20.21

43.66

x

y

Table 2. The group velocity of acoustic phonons (out-of-plane acoustic (ZA), transverse acoustic (TA), longitudinal acoustic (LA)) near Γ point for TiS3 monolayer along x and y directions. The result of MoS2 single-layer is also listed for comparison. direction

ZA (km/s)

TA (km/s)

LA (km/s)

x

0.88

3.01

5.43

y

1.11

2.31

6.16

1.40

3.96

6.47

TiS3 MoS2

12

ACS Paragon Plus Environment

Page 13 of 25

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

ACS Applied Materials & Interfaces

Figure 1. The ball-and-stick model of (a) bulk TiS3, (b) and (c) are side-view and top-view of TiS3 monolayer, respectively. The coordinate axes (x, y, z) and structural parameters (d0, d1, d2, d3) are indicated.

13

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces

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

Figure 2. (a) Calculated band structure of bulk TiS3, and the orbital-decomposed band structure of TiS3 monolayer: (b) the d orbits of Ti (red), (c) the p orbits of S2 (green), (d) the p orbits of S1 (blue). The fermi level is set to 0 eV.

14

ACS Paragon Plus Environment

Page 14 of 25

Page 15 of 25

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

ACS Applied Materials & Interfaces

Figure 3. The room temperature electronic transport coefficients: (a) absolute values of Seebeck coefficient, (b) electrical conductivity, (c) power factor of TiS3 monolayer as a function of carrier concentration from 1011 to 1014 cm−2.

15

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces

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

Figure 4. (a) The phonon dispersion relations and (b) partial phonon DOS of TiS3 monolayer, (c) is the calculated lattice thermal conductivity as a function of temperature along x and y directions.

16

ACS Paragon Plus Environment

Page 16 of 25

Page 17 of 25

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

ACS Applied Materials & Interfaces

Figure 5. The extracted (a) group velocity and (b) relaxation time of different phonon modes for TiS3 monolayer as a function of frequency, (c) is the variation of accumulative lattice thermal conductivity with respect to cutoff frequency at room temperature.

17

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces

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

Figure 6. Percentage contribution of different phonon branches to lattice thermal conductivity of TiS3 monolayer along (a) x and (b) y directions as a function of temperature.

18

ACS Paragon Plus Environment

Page 18 of 25

Page 19 of 25

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

ACS Applied Materials & Interfaces

Figure 7. The calculated (a) p-type and (b) n-type ZT values of TiS3 monolayer as a function of carrier concentration along x and y directions at 300 K.

19

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces

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

References [1] Pei, Y. Z.; Shi, X. Y.; LaLonde, A.; Wang, H.; Chen, L. D.; Snyder, G. J. Convergence of Electronic Bands for High Performance Bulk Thermoelectric. Nature

2011, 473, 66-69. [2] Liu, W.; Tan, X. J.; Yin, K.; Liu, H. J.; Tang, X. F.; Shi, J.; Zhang, Q. J.; Uher, C. Convergence of Conduction Bands as a Means of Enhancing Thermoelectric Performance of n-Type Mg2Si1-xSnx Solid Solutions. Phys. Rev. Lett. 2012, 108, 166601-5. [3] Yu, B.; Zebarjadi, M.; Wang, H.; Lukas, K.; Wang, H. Z.; Wang, D. Z.; Opeil, C.; Dresselhaus, M.; Chen, G.; Ren, Z. F. Enhancement of Thermoelectric Properties by Modulation-Doping in Silicon Germanium Alloy Nanocomposites. Nano Lett. 2012, 12, 2077-2082. [4] Hsu, K. F.; Loo, S.; Guo, F.; Chen, W.; Dyck, J. S.; Uher, C.; Hogan, T.; Polychroniadis, E. K.; Kanatzidis, M. G. Cubic AgPbmSbTe2+m: Bulk Thermoelectric Materials with High Figure of Merit. Science 2004, 303, 818-821. [5] Poudel, B.; Hao, Q.; Ma, Y.; Lan, Y. C.; Minnich, A.; Yu, B.; Yan, X.; Wang, D. Z.; Muto, A.; Vashaee, D.; Chen, X. Y.; Liu, J. M.; Dresselhaus, M. S.; Chen, G.; Ren, Z. F. High-Thermoelectric Performance of Nanostructured Bismuth Antimony Telluride Bulk Alloys. Science 2008, 320, 634-638. [6] Biswas, K.; He, J. Q.; Blum, I. D.; Wu, C. I.; Hogan, T. P.; Seidman, D. N.; Dravid, V. P.; Kanatzidis, M. G. High-performance Bulk Thermoelectrics with All-scale Hierarchical Architectures. Nature 2012, 489, 414-418. [7] Hicks, L. D.; Dresselhaus, M. S. Effect of Quantum-well Structures on the Thermoelectric Figure of Merit. Phys. Rev. B 1993, 47, 12727-12731. [ 8 ] Hicks, L. D.; Dresselhaus, M. S. Thermoelectric Figure of Merit of a One-dimensional Conductor. Phys. Rev. B 1993, 47, 16631-16634. [ 9 ] Venkatasubramanian, R.; Siivola, E.; Colpitts, T.; O’Quinn, B.; Thin-film Thermoelectric Devices with High Room-temperature Figures of Merit. Nature 2001, 413, 597-602. [10] Harman, T. C.; Taylor, P. J.; Walsh, M. P.; LaForge, B. E. Quantum Dot Superlattice Thermoelectric Materials and Devices. Science 2002, 297, 2229-2232. [11] Hochbaum, A. I.; Chen, R. K.; Delgado, R. D.; Liang, W. J.; Garnett, E. C.; Najarian, M.; Majumdar, A.; Yang, P. D. Enhanced Thermoelectric Performance of 20

ACS Paragon Plus Environment

Page 20 of 25

Page 21 of 25

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

ACS Applied Materials & Interfaces

Rough Silicon Nanowires. Nature 2008, 451, 163-167. [12] Radisavljevic, B.; Radenovic, A.; Brivio, J.; Giacometti, V.; Kis, A. Single-layer MoS2 Transistors. Nat. Nanotechnol. 2011, 6, 147-150. [ 13 ] Lopez-Sanchez, O.; Lembke, D.; Kayci, M.; Radenovic, A.; Kis, A. Ultrasensitive Photodetectors based on Monolayer MoS2. Nat. Nanotechnol. 2013, 8, 497-501. [14] Ma, J.; Chen, Y.; Han, Z.; Li, W. Strong Anisotropic Thermal Conductivity of Monolayer WTe2. 2D Mater. 2016, 3, 045010-7. [15] Pawbake, A. S.; Pawar, M. S.; Jadkar, S. R.; Late, D. J. Large Area Chemical Vapor Deposition of Monolayer Transition Metal Dichalcogenides and Their Temperature Dependent Raman Spectroscopy Studies. Nanoscale 2016, 8, 3008-3018. [16] Pawbake, A. S.; Waykar, R. G., Late, D. J.; Jadkar, S. R. Highly Transparent Wafer-Scale Synthesis of Crystalline WS2 Nanoparticle Thin Film for Photodetector and Humidity-Sensing Applications. ACS Appl. Mater. Interfaces 2016, 8, 3359-3365. [17] Wang, T.; Liu, C.; Xu, J.; Zhu, Z.; Liu, E.; Hu, Y.; Li, C.; Jiang, F. Thermoelectric Performance of Restacked MoS2 Nanosheets Thin-film. Nanotechnology 2016, 27, 285703-6. [18] Yoshida, M.; Iizuka, T.; Saito, Y.; Onga, M.; Suzuki, R.; Zhang, Y.; Iwasa, Y.; Shimizu, S. Gate-Optimized Thermoelectric Power Factor in Ultrathin WSe2 Single Crystals. Nano Lett. 2016, 16, 2061-2065. [ 19 ] Kumar, S.; Schwingenschlögl, U. Thermoelectric Response of Bulk and Monolayer MoSe2 and WSe2. Chem. Mater. 2015, 27, 1278-1284. [20] Huang, W.; Da, H.; Liang, G. Thermoelectric Performance of MX2 (M = Mo, W; X = S, Se) Monolayers. J. Appl. Phys. 2013, 113, 104304-7. [21] Island, J. O.; Barawi, M.; Biele, R.; Almázn, A.; Clamagirand, J. M.; Ares, J. R.; Sánchez, C.; Zant, H. S. van der; Álvarez, J. V.; D’Agosta, R.; Ferrer, I. J.; Castellanos-Gomez, A. TiS3 Transistors with Tailored Morphology and Electrical Properties. Adv. Mater. 2015, 27, 2595-2601. [22] Lipatov, A.; Wilson, P. M.; Shekhirev, M.; Teeter, J. D.; Netusil, R.; Sinitskii, A. Few-layered Titanium Trisulfide (TiS3) Field-effect Transistors. Nanoscale 2015, 7, 12291-12296. [23] Island, J. O.; Buscema, M.; Barawi, M.; Clamagirand, J. M.; Ares, J. R.; Sánchez, C.; Ferrer, I. J.; Steele, G. A.; Zant, H. S. van der; Castellanos-Gomez, A. Ultrahigh 21

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces

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

Page 22 of 25

Photoresponse of Few-Layer TiS3 Nanoribbon Transistors. Adv. Opt. Mater. 2014, 2, 641-645. [24] Pawbake, A. S.; Island, J. O.; Flores, E.; Ares, J. R.; Sánchez, C.; Ferrer, I. J.; Jadkar, S. R.; Zant, H. S. J. van der; Castellanos-Gomez, A.; Late, D. J. Temperature-Dependent Raman Spectroscopy of Titanium Trisulfide

(TiS3)

Nanoribbons and Nanosheets. ACS Appl. Mater. Interfaces 2015, 7, 24185-24190. [25] Dai, J.; Zeng, X. C. Titanium Trisulfide Monolayer: Theoretical Prediction of a New Direct-Gap Semiconductor with High and Anisotropic Carrier Mobility. Angew. Chem. Int. Ed. 2015, 54, 7572-7576. [26] Dai, J.; Li, M.; Zeng, X. C. Group IVB Transition Metal Trichalcogenides: a New Class of 2D Layered Materials Beyond Grapheme. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2016, 6, 211-222. [27] Kresse, G.; Hafner, J. Ab. initio Molecular Dynamics for Liquid Metals. Phys. Rev. B 1993, 47, R558-561. [ 28 ] Kresse, G.; Hafner, J. Ab initio Molecular-dynamics Simulation of the Liquid-metal-amorphous-semiconductor Transition in Germanium. Phys. Rev. B 1994, 49, 14251-14269. [29] Kresse, G.; Furthmüller, J. Efficiency of Ab-initio Total Energy Calculations for Metals and Semiconductors Using a Plane-wave Basis Set. Comput. Mater. Sci. 1996, 6, 15-50. [30] 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. [31] Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865-3868. [32] Monkhorst, H. J.; Pack, J. D. Special Points for Brillonin-zone Integrations. Phys. Rev. B 1976, 13, 5188-5192. [33] Heyd, J.; Scuseria, G. E.; Ernzerhof, M. Hybrid Functionals Based on a Screened Coulomb Potential. J. Chem. Phys. 2003, 118, 8207-8215. [34] Madsen, G. K.H.; Singh, D. J. BoltzTraP. A code for Calculating Band-structure Dependent Quantities. Comput. Phys. Commun. 2006, 175, 67-71. [35] Togo, A.; Oba, F.; Tanaka, I. First-principles Calculations of the Ferroelastic Transition between Rutile-type and CaCl2-type SiO2 at High Pressures. Phys. Rev. B

2008, 78, 134106-9. 22

ACS Paragon Plus Environment

Page 23 of 25

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

ACS Applied Materials & Interfaces

[36] Li, W.; Carrete, J.; Katcho, N. A.; Mingo, N. ShengBTE: A Solver of the Boltzmann Transport Equation for Phonons. Comput. Phys. Commun. 2014, 185, 1747-1758. [37] Li, W.; Mingo, N.; Lindsay, L.; Broido, D. A.; Stewart, D. A.; Katcho, N. A. Thermal Conductivity of Diamond Nanowires from First Principles. Phys. Rev. B

2012, 85, 195436-5. [38] Li, W.; Lindsay, L.; Broido, D. A.; Stewart, D. A.; Mingo, N. Thermal Conductivity of Bulk and Nanowire Mg2SixSn1−x Alloys from First Principles. Phys. Rev. B 2012, 86, 174307-8. [39] Grimme, S. Semiempirical GGA-type Density Functional Constructed with a Long-range Dispersion Correction. J. Comput. Chem. 2006, 27, 1787-1799. [40] Klimeš, J.; Bowler, D. R.; Michaelides, A. Van der Waals Density Functionals Applied to Solids. Phys. Rev. B 2011, 83, 195131-13. [41] Furuseth, S.; Brattas, L.; Kjekshus, A. On the Crystal Structures of TiS3, ZrS3, ZrSe3, ZrTe3, HfS3, and HfSe3. Acta Chem. Scand. 1975, 29, 623-621. [42] Wu, J.; Wang, D.; Liu, H.; Lau, W.-M.; Liu, L.-M. An ab initio Study of TiS3: a Promising Electrode Material for Rechargeable Li and Na Ion Batteries. RSC Adv.

2015, 5, 21455-21463. [43] Grimmeiss, H. G.; Rabenau, A.; Hahn, H.; Ness, P. Electrical and Optical Properties of Chalcogenides of the Elements of Group IVB. Z. Elektrochem. 1961, 65, 776-783. [44] Kang, J.; Wang, L.-W. Robust Band Gap of TiS3 Nanofilms. Phys. Chem. Chem. Phys. 2016, 18, 14805-14809. [45] Kang, J.; Sahin, H.; Ozaydin, H. D.; Senger, R. T.; Peeters, F. M. TiS3 Nanoribbons: Width-independent Band Gap and Strain-tunable Electronic Properties. Phys. Rev. B 2015, 92, 075413-8. [46] Bardeen, J.; Shockley, W. Deformation Potentials and Mobilities in Non-Polar Crystals. Phys. Rev. B 1950, 80, 72-80. [47] Xi, J. Y.; Long, M. Q.; Tang, L.; Wang, D.; Shuai, Z. G. First-principles Prediction of Charge Mobility in Carbon and Organic Nanomaterials. Nanoscale 2012, 4, 4348-4369. [48] Cai, Y.; Zhang, G.; Zhang, Y.-W. Polarity-Reversed Robust Carrier Mobility in Monolayer MoS2 Nanoribbons. J. Am. Chem. Soc. 2014, 136, 6269-6275. 23

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces

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

[49] Snyder, G. J.; Toberer, E. S. Complex Thermoelectric Materials. Nat. Mater. 2008,

7, 105-114. [50] Li, W.; Carrete, J.; Mingo, N. Thermal Conductivity and Phonon Linewidths of Monolayer MoS2 from First Principles. Appl. Phys. Lett. 2013, 103, 253103-4. [ 51 ] Gandi, A. N.; Schwingenschlögl, U. Thermal Conductivity of Bulk and Monolayer MoS2. EPL 2016, 113, 36002-6. [ 52 ] Gu, X.; Yang, R. Phonon Transport in Single-layer Transition Metal Dichalcogenides: A First-principles Study. Appl. Phys. Lett. 2014, 105, 131903-5. [53] Zhang, J.; Liu, H. J.; Cheng, L.; Wei, J.; Liang, J. H.; Fan, D. D.; Jiang, P. H.; Shi, J. Thermal Conductivities of Phosphorene Allotropes from First-principles Calculations: A Comparative Study. arXiv, 2016, 1601.07302. [54] Qin, G.; Yan, Q.-B.; Qin, Z.; Yue, S.-Y.; Hu, M.; Su, G. Anisotropic Intrinsic Lattice Thermal Conductivity of Phosphorene from First Principles. Phys. Chem. Chem. Phys. 2015, 17, 4854-4858.

24

ACS Paragon Plus Environment

Page 24 of 25

Page 25 of 25

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

ACS Applied Materials & Interfaces

TOC

25

ACS Paragon Plus Environment