Excellent Thermoelectric Performance Predicted in Two-Dimensional

May 22, 2017 - Shu-Shen Lyu,. †. Xiao-Ming Wang,. ‡. Yuan-Xiang Fu,. †. Yi Heng,. † and Dong-Chuan Mo*,†. †. School of Chemical Engineerin...
0 downloads 0 Views 2MB Size
Subscriber access provided by UNIV OF ARIZONA

Article

Excellent Thermoelectric Performance Predicted in TwoDimensional Buckled Antimonene: a First-Principles Study Kai-Xuan Chen, Shu-Shen Lyu, Xiao-Ming Wang, Yuan-Xiang Fu, Yi Heng, and Dong-Chuan Mo J. Phys. Chem. C, Just Accepted Manuscript • Publication Date (Web): 22 May 2017 Downloaded from http://pubs.acs.org on May 30, 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.

The Journal of Physical Chemistry C 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 23

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

The Journal of Physical Chemistry

Excellent Thermoelectric Performance Predicted in Two-dimensional Buckled Antimonene: a Firstprinciples Study Kai-Xuan Chena, Shu-Shen Lyua, Xiao-Ming Wangb, Yuan-Xiang Fua, Yi Henga and Dong-Chuan Moa, * a

School of Chemical Engineering and Technology, Sun Yat-sen University, Guangzhou 510275,

China. b

Department of Physics and Astronomy, University of Toledo, 2801 W Bancroft ST, Toledo, OH

43606, USA.

ABSTRACT Nowadays, new emerging two-dimensional (2D) materials have become a hot topic in the field of theoretical physics, material science and nanotechnology engineering due to their high surface area, planar structure and quantum confinement effect. Within two-dimensional framework, we systematically concentrate on the buckled and puckered systems consisting of the VA group elements (denoted as arsenene, antimonene and bismuthene). Among these studied systems, the buckled antimonene harbors a thermoelectric figure of merit (ZT) of 2.15 at room temperature.

ACS Paragon Plus Environment

1

The Journal of Physical Chemistry

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 23

This is probably the highest value that has ever been reported in pristine 2D materials. By simple biaxial strain engineering, the ZT can even get enhanced to 2.9 under 3% tensile strain. The enhancement mainly results from both tuning the electronic structures and reducing the thermal conductance. This work predicts a new promising candidate in thermoelectric devices, based on the fact that buckled antimonene has been lately fabricated and proved to be stable at ambient conditions.

INTRODUCTION Thermoelectric materials, which can be used in the energy conversion between heat and electricity, have attracted a lot of attention due to the global energy crisis. To evaluate the performance of a thermoelectric material, a dimensionless parameter ZT1 is adopted with the definition of

σ S 2T κ

(1)

κ = κ el + κ ph

(2)

ZT =

where σ and S are the electronic conductance and Seebeck coefficient, respectively. κ denotes the thermal conductance which can be contributed both by phonons and electrons, respectively.2 In the commercial field, frequently-used thermoelectric materials have a ZT of around unit, which indicates the low energy conversion efficiency. For many years, the development in high performance thermoelectric materials meets its bottleneck since it is difficult to alter any one parameter without greatly affecting the other transport coefficients.3 Then the research on

ACS Paragon Plus Environment

2

Page 3 of 23

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

The Journal of Physical Chemistry

thermoelectric materials witnessed great progress in 1993, when Hicks and Dresselhaus3-4 pointed out an effective way to significantly enhance the ZT, by nanostructural engineering. Because in nanoscale size, the interrelation among σ, S and κ may be weakened, due to the quantum confinement. In the frontier experimental research of thermoelectric materials, progress has been made. Many works devoted efforts to reducing the thermal conductivity by nanostructural engineering. Some other studies focused on the band engineering or mobility enhancement, from the point view of tuning the electronic properties. For instances, Venkatasubramanian et al.5 reported a ZT of ~2.4 observed in thin-film p-type Bi2Te3/Sb2Te3 superlattice devices by controlling the transport of phonons and electrons. Yan et al.6 fabricated n-type Bi2Te2.7Se0.3 sample (ZT = 1.04 at 380 K) by the hot pressing of ball-milled powder mainly due to the large increase of electrical conductivity. In the work of Guo et al.7, Tl9BiTe6, was postulated to exceed ZT = 1 above 450 K when prepared by zone-melting and reach ZT = 0.86 at 560 K after hot-pressing via the optimization of the carrier concentration. Zhao et al.2 reported a record of ZT ranging from 0.7 to 2.0 at 300 to 773 K, realized in hole-doped tin SnSe crystals enabled by the contribution of multiple electronic valence bands present in SnSe. Nevertheless, it should be noticed that these reported ZT numbers are still far away from the level which promises sufficient energy conversion efficiency in practical thermoelectric devices. In experiments, 2D materials are more likely to fabricate in large length scale and easier to get used in electronic devices.8 Moreover, they offer great flexibility in terms of tuning their electronic properties.9 In recent years, some research based on two-dimensional thermoelectric materials have been carried on, as listed in Table. 1. However, the challenge still remains in seeking high-performance thermoelectric materials in nanoscale, especially 2D materials.

ACS Paragon Plus Environment

3

The Journal of Physical Chemistry

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 23

Therefore, the objective of this work is to design and predict the thermoelectric properties of new proposed 2D materials, which consist of the VA group elements (As, Sb and Bi). Two distinguishable structural forms in this 2D family, denoted as the buckled and puckered systems, are investigated. Table. 1 The room-temperature ZTmax for typical two-dimensional materials. 2D Monolayer

ZTmax at RT

Graphene

0.009410

Graphyne

0.15710

Silicene

0.3611

Germanene

0.4111

Black Phosphorene

>0.612, 1.4411

Puckered Arsenene

0.8513

MoS2

0.7514, 0.5815, 1.3516

MoSe2

0.8814, 1.3916

WS2

0.7214, 1.5216

WSe2

0.9114, 1.8816

Buckled Antimonene

2.15 (This work)

COMPUTATIONAL METHODS Density functional theory (DFT), as implemented in the package Quantum Espresso17, is adopted in our work. Standard solid state pseudopotentials18 (SSSP) within Perdew–Burke– Ernzerhof (PBE) framework are used. A 15×15×1 Monkhorst–Pack k-mesh is employed with

ACS Paragon Plus Environment

4

Page 5 of 23

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

The Journal of Physical Chemistry

energy (charge density) cutoff up to 40 Ry (320 Ry) and a vacuum region of more than 14 Angstroms is used to avoid the periodic image interaction. Convergence thresholds are set to 10-6 Ry/Bohr on force and 10-5 Ry on total energy. Optimized structural parameters are obtained after ion minimization has been fully performed.

For the calculation of electronic structures, ultrasoft pseudopotentials with the Rappe−Rabe− Kaxiras − Joannopoulos (RRKJ) method are adopted with same energy cutoff in ion minimization. The d-orbital electrons are fully taken into consideration. In the process of selfconsistency, a 7×7×1 Monkhorst–Pack k-mesh is employed. The ballistic electron and phonon transports are studied in typical non-equilibrium Green’s function (NEGF) model “LeadConductor-Lead”.19 The interaction of electron–electron, electron–phonon and phonon–phonon are not taken into consideration. Maximally localized Wannier functions (MLWFs) method, as implemented in WanT20, is used to calculated the Bloch states at each k-point in Brillouin zone. To obtain the electron Hamiltonian, 46 Wannier functions are used in buckled systems and 126 functions for the puckered case. The calculation of thermoelectric factors are described as below.14, 21 For ballistic electronic transport, one should firstly calculate the retarded Green’s function of the central conductor as

G r =  ESC − H C − Σ rL − Σ rR 

−1

(3)

where E, HC and SC are the electron energy, Hamiltonian and overlap matrix, respectively. Σr is the self-energy term from the semi-infinite lead, which can be obtained as † † Σ rL = H LC g Lr H LC , Σ rR = H CR g Rr H CR

ACS Paragon Plus Environment

(4)

5

The Journal of Physical Chemistry

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 23

where gr denotes the retarded surface Green’s function from semi-infinite lead. HLC (HCR) represent the Hamiltonian matrix between the left (right) leads and the central conductor. After then, the electronic transmittance matrix T ( E ) can be calculated as Γ β = i ( Σ rβ − Σ aβ ) , β = L, R T ( E ) = Tr ( G r Γ LG a Γ R ) , G a = ( G r )

(5)



(6)

Once the electronic transmittance matrix has been obtained, several thermoelectric factors can be worked out. For convenience, Lorenz functions Ln are introduced as

Ln ( µ , T ) =

 ∂f ( E , µ , T )  2 n µ dE T E × E − × ( ) ( ) −  h∫ ∂E   1

f ( E, µ,T ) = e

E −µ k BT

(7)

(8)

+1

where f(E, µ, T) is the Fermi-Dirac distribution function and µ, T, h are the chemical potential, absolute temperature and Planck constant, respectively. Here, The electronic conductance σ, Seebeck coefficient S and electronic thermal conductance κel can be worked out by Equations (9) -(11), respectively.

σ = q 2 L0

S=

1 L1 × qT L0

ACS Paragon Plus Environment

(9)

(10)

6

Page 7 of 23

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

The Journal of Physical Chemistry

κ el =

1  L2  ×  L2 − 1  T  L0 

(11)

From the definition of ZT in Equations (1)-(2), now there is only one physical parameter left to be mentioned, the phononic thermal conductance term, which is determined from phononic transport. For ballistic phononic transport, the calculation of phonon transmittance is similar with that of the electron transmittance described above. The interatomic force constant (FC) matrix can be obtained firstly from density functional perturbation theory (DFPT). And then the phonon transmittance T (ω ) can be calculated as

2 G r = (ω + iη ) − FC − Σ rL − Σ rR   

−1

T (ω ) = Tr ( G r Γ LG a Γ R )

(12)

(13)

where ω is the phonon frequency. Then the phononic thermal conductance κph can be obtained as

κ ph (T ) =

h ∞ ∂g (ω,T ) T (ω)ω dω ∫ 2π 0 ∂T

1

g (ω , T ) = e

hω k BT

(14)

(15)

−1

where g(ω, T) is the Bose-Einstein distribution function, in which h and kB denote the reduced Planck constant and Boltzmann constant, respectively.

RESULTS AND DISCUSSION

ACS Paragon Plus Environment

7

The Journal of Physical Chemistry

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 8 of 23

A. Structural Information and Phonon Transport. To begin with, we shall introduce the atomistic structures of the 2D family consisting of VA group elements (As, Sb and Bi). As demonstrated in Fig. 1 and Fig. 2, there are two distinguishable structural forms for this 2D family. The first one has a silicene-like crystal structure, denoted as the buckled system, whereas the second one is similar with the black phosphorene22-23, denoted as the puckered system. The buckled systems belong to the D3d group with two non-coplanar atoms included in a hexagonal unit cell. The puckered systems, however, belong to the C2v group with a cubic unit cell, in which four atoms are included in the unit cell. Besides, structural anisotropy in the puckered systems can be observed along armchair and zigzag directions. The detailed crystal parameters of these two forms of 2D materials are tabulated in Table. 2 and Table. 3, respectively, which agree well with the previous research.24-29

Fig. 1 The atomic structures of buckled systems, in which the unit cell is denoted with red solid box. Table. 2 The crystal parameters of the buckled systems.

Structures

a/Å

R/Å

ACS Paragon Plus Environment

θ/°

8

Page 9 of 23

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

The Journal of Physical Chemistry

Buckled Arsenene

3.61

2.51

92.1

Buckled Antimonene

4.12

2.89

91.1

Buckled Bismuthene

4.29

3.03

90.0

Fig. 2 The atomic structures of puckered systems, in which the unit cell is denoted with red solid box. Table. 3 The crystal parameters of the puckered systems.

Structures

a/Å

b/Å

R1 / Å

R2 / Å

θ1 / °

θ2 / °

Puckered Arsenene

4.76

3.68

2.51

2.49

94.6

100.6

Puckered Antimonene

4.75

4.35

2.94

2.86

95.4

102.5

Puckered Bismuthene

4.73

4.51

3.12

3.04

92.5

104.7

It is worth mentioning that for puckered antimonene and bismuthene, there also exist two different kinds of structures, w-Sb/Bi and aw-Sb/Bi.

30

The difference between w-Sb/Bi and aw-

ACS Paragon Plus Environment

9

The Journal of Physical Chemistry

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 10 of 23

Sb/Bi can be found in zigzag views from Fig. 3. It can be distinguished by whether the atoms are coplanar along the up and down surfaces, as shown in Fig. 3(a)-(b). We then simply perform the phonon calculation and plot out the phonon dispersion spectra. Imaginary frequency mode appears in the w-Sb/Bi, indicating the thermodynamic instability. On the contrary, no obvious negative phonon branches can be observed in the aw-Sb/Bi. Consequently, the aw-Sb/Bi is adopted in our study and in the following part, for convenience, aw-Sb/Bi is described as the puckered system.

Fig. 3 The two different atomic structures of puckered systems for (a, b) antimonene and (c, d) bismuthene

ACS Paragon Plus Environment

10

Page 11 of 23

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

The Journal of Physical Chemistry

Fig. 4 The atomistic structures, phonon dispersion and transmittance spectrum of 2 × 2 supercell of buckled (a) arsenene, (b) antimonene, (c) bismuthene and puckered (d) arsenene, (e) antimonene, (f) bismuthene. Fig. 4 shows the atomic structures and the corresponding phonon dispersion spectra of all investigated systems. Reciprocal path along the high symmetry point “G(0, 0, 0)-M(1/2, 0, 0)K(1/3, 1/3, 0)-G(0, 0, 0)” is used for buckled systems and “G(0, 0, 0)-X(1/2, 0, 0)-S(1/2, -1/2, 0)-

ACS Paragon Plus Environment

11

The Journal of Physical Chemistry

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 23

Y(0, -1/2, 0)-G(0, 0, 0)” for puckered systems. Here, a 2×2 supercell is employed instead of the unit cell in order to meet the criteria of nearest-neighbor interaction in the following NEGF part. As can be observed from Fig. 4, the cut-off frequency in dispersion spectra decreases with increasing atomic mass, both in buckled and puckered systems. To be specific, from arsenene to bismuthene, the values decrease from ~300 cm-1 to ~150 cm-1. Also, phononic bandgap can be observed in all these systems. It should be mentioned that only negligible imaginary frequencies can be observed in all these systems, indicating the thermodynamic stability of the proposed structures. This provides a theoretical guarantee for the further study of the electronic and thermoelectric properties. As is known, the low-frequency modes mainly dominate the room-temperature thermal transport properties of 2D materials since high-frequency modes cannot be activated at such temperature.19 As for the puckered systems, due to the anisotropy along the armchair and zigzag directions, the phonon transmittance also exhibits anisotropic phenomenon, as shown in Fig.

4(d)-(f). There is no doubt that the phonon transmittance along the zigzag direction is higher than that along the armchair direction12, 31. Such anisotropy in phonon transmittance would definitely lead to the same trend in phononic thermal conductance.

B. Electronic Structures. In the part of electronic structures, MLWFs are adopted to describe the Bloch states at each k-point in the Brillouin zone. A frozen window of [-2.0 eV, 2.0 eV] is chosen in which the subbands are correctly reproduced compared with the bands plotted from DFT data, as demonstrated in Fig. 5. The band structures of buckled antimonene plotted by MLWFs and DFT fit well within the frozen window. It provides strong evidence for the accurate description of MLWFs method in studying electronic structures of these system. Based on this,

ACS Paragon Plus Environment

12

Page 13 of 23

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

The Journal of Physical Chemistry

Fig. 6 shows the electronic band structures, transmittance and density of states. All investigated systems exhibit semiconductive properties. Wider electronic bandgaps can be observed in the buckled systems than that in the puckered counterparts. From Fig. 6(a)-(c), the bandgaps seem to decrease with increasing atomic mass, with the value of 1.58 eV, 1.25 eV and 0.56 for buckled arsenene, antimonene and bismuthene, respectively. The bandgaps of 0.76 eV, 0.1 eV, 0.1 eV can be found for puckered arsenene, antimonene and bismuthene, respectively, from Fig. 6(d)(f). The electronic transmittance demonstrates weak anisotropy along the armchair and zigzag directions.13 That is to say the structural anisotropy has more impact on the phononic transport than the electronic transport.

Fig. 5 The electronic band structures of buckled antimonene plotted by MLWFs and DFT, respectively.

ACS Paragon Plus Environment

13

The Journal of Physical Chemistry

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 14 of 23

Fig. 6 The electronic band structure and transmittance spectrum of 2×2 supercell of buckled (a) arsenene, (b) antimonene, (c) bismuthene and puckered (d) arsenene, (e) antimonene, (f) bismuthene. Normal DFT calculation generally underestimates the bandgap, compared to that by GW or hybrid functional calculations. However, since the GW or hybrid functional correction won’t change the effective mass greatly, it may not affect the maximal ZT32, but will definitely change the chemical potential at which the ZT is maximized. Since the thermoelectric performance is our main concern in this work and due to the limited resources, the employment of DFT method is reasonable and effective. Moreover, the calculated results agree well with the previous research. For instances, the electronic structures of arsenene are consistent with the research carried out by Kecik et al.33 and Cao et al34 separately. In the work of Aktürk et al.28, the stability and electronic properties of 2D antimonene were studied and the bandgaps for buckled and puckered antimonene were found to be 1.04 eV and 0.16 eV, respectively. Afterwards, buckled

ACS Paragon Plus Environment

14

Page 15 of 23

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

The Journal of Physical Chemistry

(puckered) bismuthene were figured out to harbor a direct (indirect) bandgap of 0.547 eV (0.28 eV) 30. It is worth mentioning that for puckered antimonene and bismuthene, the lowest bandgap doesn’t locate along the high-symmetry path “G-X-S-Y-G”, but at the k-point (denoted as A) of ~(0.04, 0.25, 0) and ~(0.055, 0.255, 0), respectively. These certain k-points are figured out by using a dense mesh in the band structures process.

C. Thermoelectric performance. In the final section, we turn to the study of thermoelectric performance. Fig. 7 shows the thermoelectric factors of buckled and puckered systems at room temperature. The nature of electronic bandgap leads to a low electronic conductance at low chemical potential, at which the Fermi surface locates between the valence and conduction subbands. Even though the Wiedemann–Franz law doesn’t work precisely for semiconductor, it can still be found that the electronic thermal conductance shows similar trend with the electronic conductance. This phenomenon mainly results from their similar dependence on the electronic transmittance. At the same chemical potential, the puckered systems harbor higher electronic conductance than the buckled counterparts due to the narrower bandgaps. According to the work of Johnson et al35 and Fan et al36, the relation between Seebeck coefficient S and bandgap Eg can be described as S ≈ − ( k B / e ) ( Eg / 2k BT + 2 ) , where kB is the Boltzmann constant. Higher Seebeck coefficient can generally be discovered in systems with higher bandgap, as shown in Fig. 7(b). It may qualitatively account for the higher Seebeck values appear in buckled systems rather than the puckered counterparts. As for phononic thermal conductance in Fig. 7(c), obviously the buckled systems harbor higher values than the puckered counterparts, which may be attributed to the suppression in the phonon dispersion spectra in Fig. 4. In particular, the value for puckered systems along armchair direction are smaller than that along zigzag direction, resulting from the strong anisotropy that has been mentioned above.

ACS Paragon Plus Environment

15

The Journal of Physical Chemistry

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 16 of 23

Fig. 7 The (a) electronic conductance σ, (b) Seebeck coefficient S, (c) phonon thermal conductance κph and (d) ZT of the investigated buckled arsenene, antimonene and bismuthene, respectively. The inset in (a) shows the electronic thermal conductance κel. AC and ZZ denote the armchair and zigzag directions, respectively.

ACS Paragon Plus Environment

16

Page 17 of 23

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

The Journal of Physical Chemistry

Fig. 8 (a) The thermoelectric figure of merit of buckled antimonene with biaxial strain deformation. (b) The bandgap and phononic thermal conductance of buckled antimonene with biaxial strain deformation. From buckled to puckered systems, a degeneration of ZT can be observed from Fig. 7(d). This degeneration may mainly due to the decrease of Seebeck coefficient (resulting from the narrower bandgaps). In particular, the buckled antimonene harbors a ZT of 2.15 at room temperature, which is probably the highest value ever reported in pristine 2D materials, as presented in Table.

1. Moreover, in order to further improve the thermoelectric performance, strain engineering is used.37 The biaxial strain is applied to the buckled antimonene at the tensile range of [0%, 8%]. The effect of biaxial strain on the ZT value is significant. As shown in Fig. 8(a), the maximal room-temperature ZT of buckled antimonene get enhanced at low strain but suppressed when strain is larger than 3%. The ZT of 2.9 can be reached at the strain level of 3%, resulting in a 35% enhancement compared to that of the unstrained buckled antimonene. The enhancement mechanism may result from tuning the electronic structures and reducing the phononic thermal conductance by tensile biaxial strain. As shown in Fig. 8(b), the bandgap of buckled antimonene

ACS Paragon Plus Environment

17

The Journal of Physical Chemistry

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 18 of 23

increases at low strain deformation and then decreases when tensile strain is larger than 2%. On the contrary, the phononic thermal conductance get reduced at low strain deformation [0%, 3%] and enhanced at high strain [3%, 8%]. According to the definition of ZT in equation (1), high Seebeck coefficient (generally influenced by electronic bandgaps) and low thermal conductance are both required to achieve a high ZT. In our work, the highest ZT of 2.9 would appears when appropriate strain (3%) is applied to buckled antimonene. Generally in the experimental procedure, some factors may limit the intrinsic thermoelectric properties, such as sample size, atmoistic layer nuber, substrate effect, doping level and so on. So the obtained ZT may be lower than the theorectical value to some extent. But with the advancement of tecnology in nanostructural engineering, the reducement may be weakened. Very recently, the few-layers buckled antimonene has been successfully fabricated and proved to be stable at ambient conditions,38-39 which provide strong experimental support. Samples with desired systems are likely to be prepared in the near future. This would definitely make buckled antimonene a star material in the field of thermoelectric application. Our theoretical prediction highlights an excellent candidate in thermoelectric devices for experimental researchers. It may greatly promote the development of high efficiency thermoelectric devices based on these new proposed 2D materials.

CONCLUSIONS In conclusion, by using first-principles calculations, we study the thermoelectric transport properties of two-dimensional family which consists of VA elements (As, Sb and Bi). Both buckled and puckered forms are investigated. All the systems show semiconductive properties

ACS Paragon Plus Environment

18

Page 19 of 23

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

The Journal of Physical Chemistry

while the buckled systems possess wider bandgap than the puckered counterparts. Our results find out that the buckled systems beat the puckered counterparts in thermoelectric performance and in particular, buckled antimonene exhibits amazingly high thermoelectric performance, with maximal room-temperature ZT reaching 2.15. Imposing biaxial strain deformation, the ZT value can even get enhanced to 2.9, which make buckled antimonene probably the most excellent 2D thermoelectric material ever reported. We believe that by ways of hole/substitution doping or cutting into nanoribbons, the ZT may get improved to a much higher level, which is totally practical in industrial application.

AUTHOR INFORMATION

Corresponding Author *Tel: +86-020-84112151. E-mail: [email protected].

Author Contributions The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript.

Notes The authors declare no competing financial interest. ACKNOWLEDGMENT Financial support from the National Natural Science Foundation of China (Grant No. 51676212) and the Fundamental Research Funds for the Central Universities are gratefully acknowledged.

ACS Paragon Plus Environment

19

The Journal of Physical Chemistry

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 20 of 23

The simulation work is supported by the National Supercomputer Center in Guangzhou and the high-performance grid computing platform of Sun Yat-sen University. ABBREVIATIONS 2D, two-dimensional; ZT, thermoelectric figure of merit; DFT, density functional theory; SSSP, standard solid state pseudopotentials; PBE, Perdew–Burke–Ernzerhof; RRKJ, Rappe−Rabe− Kaxiras−Joannopoulos; NEGF, Non-equilibrium Green’s function; MLWF, maximally localized Wannier function; DFPT, density functional perturbation theory. REFERENCES 1.

Page, A.; Van der Ven, A.; Poudeu, P. F. P.; Uher, C., Origins of phase separation in thermoelectric (Ti, Zr, Hf)NiSn half-Heusler alloys from first principles. J. Mater. Chem. A. 2016, 4, 13949-13956. 2. Zhao, L. D., et al., Ultrahigh power factor and thermoelectric performance in hole-doped single-crystal SnSe. Science 2016, 351, 141-144. 3. Hicks, L. D.; Dresselhaus, M. S., Effect of quantum-well structures on the thermoelectric figure of merit. Phys. Rev. B 1993, 47, 12727-12731. 4. Hicks, L. D.; Dresselhaus, M. S., Thermoelectric figure of merit of a one-dimensional conductor. Phys. Rev. B 1993, 47, 16631-16634. 5. 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. 6. Yan, X.; Poudel, B.; Ma, Y.; Liu, W. S.; Joshi, G.; Wang, H.; Lan, Y.; Wang, D.; Chen, G.; Ren, Z. F., Experimental studies on anisotropic thermoelectric properties and structures of ntype Bi2Te2.7Se0.3. Nano Lett. 2010, 10, 3373-3378. 7. Guo, Q.; Chan, M.; Kuropatwa, B. A.; Kleinke, H., Enhanced thermoelectric properties of variants of Tl9SbTe6 and Tl9BiTe6. Chem. Mater. 2013, 25, 4097-4104. 8. Li, H.; Ying, H.; Chen, X.; Nika, D. L.; Cocemasov, A. I.; Cai, W.; Balandin, A. A.; Chen, S., Thermal conductivity of twisted bilayer graphene. Nanoscale 2014, 6, 13402-13408. 9. Novoselov, K. S.; Mishchenko, A.; Carvalho, A.; Castro Neto, A. H., 2D materials and van der Waals heterostructures. Science 2016, 353, aac9439. 10. Wang, X. M.; Mo, D. C.; Lu, S. S., On the thermoelectric transport properties of graphyne by the first-principles method. J. Chem. Phys. 2013, 138, 204704. 11. Lv, H. Y.; Lu, W. J.; Shao, D. F.; Sun, Y. P., Enhanced thermoelectric performance of phosphorene by strain-induced band convergence. Phys. Rev. B 2014, 90, 085433. 12. Fei, R.; Faghaninia, A.; Soklaski, R.; Yan, J. A.; Lo, C.; Yang, L., Enhanced thermoelectric efficiency via orthogonal electrical and thermal conductances in phosphorene. Nano Lett. 2014, 14, 6393-6399.

ACS Paragon Plus Environment

20

Page 21 of 23

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

The Journal of Physical Chemistry

13. Medrano Sandonas, L.; Teich, D.; Gutierrez, R.; Lorenz, T.; Pecchia, A.; Seifert, G.; Cuniberti, G., Anisotropic thermoelectric response in two-dimensional puckered structures. J. Phys. Chem. C 2016, 120, 18841-18849. 14. Chen, K. X.; Wang, X. M.; Mo, D. C.; Lyu, S. S., Thermoelectric properties of transition metal dichalcogenides: From monolayers to nanotubes. J. Phys. Chem. C 2015, 119, 2670626711. 15. Huang, W.; Da, H.; Liang, G., Thermoelectric performance of MX2 (M = Mo,W; X = S,Se) monolayers. J. Appl. Phys. 2013, 113, 104304. 16. Wickramaratne, D.; Zahid, F.; Lake, R. K., Electronic and thermoelectric properties of fewlayer transition metal dichalcogenides. J. Chem. Phys. 2014, 140, 124710. 17. Giannozzi, P., et al., QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials. J. Phys.: Condens. Matter 2009, 21, 395502. 18. Lejaeghere, K., et al., Reproducibility in density functional theory calculations of solids. Science 2016, 351, aad3000. 19. Jiang, J. W.; Wang, J. S.; Li, B. W., A nonequilibrium Green's function study of thermoelectric properties in single-walled carbon nanotubes. J. Appl. Phys. 2011, 109, 014326. 20. Ferretti, A.; Bonferroni, B.; Calzolari, A.; Nardelli, M. B. WanT code. http://www.wanniertransport.org (accessed Apr 4, 2017). 21. Chen, K. X.; Luo, Z. Y.; Mo, D. C.; Lyu, S. S., WSe2 nanoribbons: new high-performance thermoelectric materials. Phys. Chem. Chem. Phys. 2016, 18, 16337-16344. 22. Liu, T. H.; Chang, C. C., Anisotropic thermal transport in phosphorene: effects of crystal orientation. Nanoscale 2015, 7, 10648-10654. 23. 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, 48544858. 24. Wang, G.; Pandey, R.; Karna, S. P., Atomically thin group v elemental films: theoretical investigations of antimonene allotropes. ACS Appl. Mater. Interfaces 2015, 7, 11490-11496. 25. Zhang, S.; Yan, Z.; Li, Y.; Chen, Z.; Zeng, H., Atomically thin arsenene and antimonene: semimetal-semiconductor and indirect-direct band-gap transitions. Angew. Chem. Int. Ed. Engl. 2015, 54, 3112-3115. 26. Sun, J. T.; Huang, H.; Wong, S. L.; Gao, H. J.; Feng, Y. P.; Wee, A. T., Energy-gap opening in a Bi(110) nanoribbon induced by edge reconstruction. Phys. Rev. Lett. 2012, 109, 246804. 27. Lu, Y., et al., Topological properties determined by atomic buckling in self-assembled ultrathin Bi(110). Nano Lett. 2015, 15, 80-87. 28. Aktürk, O. Ü.; Özçelik, V. O.; Ciraci, S., Single-layer crystalline phases of antimony: Antimonenes. Phys. Rev. B 2015, 91, 235446. 29. Kecik, D.; Durgun, E.; Ciraci, S., Optical properties of single-layer and bilayer arsenene phases. Phys. Rev. B 2016, 94, 205410. 30. Aktürk, E.; Aktürk, O. Ü.; Ciraci, S., Single and bilayer bismuthene: Stability at high temperature and mechanical and electronic properties. Phys. Rev. B 2016, 94, 014115. 31. Zeraati, M.; Vaez Allaei, S. M.; Abdolhosseini Sarsari, I.; Pourfath, M.; Donadio, D., Highly anisotropic thermal conductivity of arsenene: Anab initiostudy. Phys. Rev. B 2016, 93, 085424. 32. Pei, Y.; Wang, H.; Snyder, G. J., Band engineering of thermoelectric materials. Adv. Mater. 2012, 24, 6125-6135.

ACS Paragon Plus Environment

21

The Journal of Physical Chemistry

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 23

33. Kecik, D.; Durgun, E.; Ciraci, S., Stability of single-layer and multilayer arsenene and their mechanical and electronic properties. Phys. Rev. B 2016, 94, 205409. 34. Cao, H. W.; Yu, Z. Y.; Lu, P. F., Electronic properties of monolayer and bilayer arsenene under in-plain biaxial strains. Superlattices Microstruct. 2015, 86, 501-507. 35. Johnson, V. A.; Lark-Horovitz, K., Theory of thermoelectric power in semiconductors with applications to germanium. Phys. Rev. 1953, 92, 226-232. 36. Fan, D. D.; Liu, H. J.; Cheng, L.; Jiang, P. H.; Shi, J.; Tang, X. F., MoS2 nanoribbons as promising thermoelectric materials. Appl. Phys. Lett. 2014, 105, 133113. 37. Yeo, P. S. E.; Sullivan, M. B.; Loh, K. P.; Gan, C. K., First-principles study of the thermoelectric properties of strained graphene nanoribbons. J. Mater. Chem. A. 2013, 1, 10762. 38. Ares, P.; Aguilar-Galindo, F.; Rodriguez-San-Miguel, D.; Aldave, D. A.; Diaz-Tendero, S.; Alcami, M.; Martin, F.; Gomez-Herrero, J.; Zamora, F., Mechanical isolation of highly stable antimonene under ambient conditions. Adv. Mater. 2016, 28, 6332-6336. 39. Ji, J., et al., Two-dimensional antimonene single crystals grown by van der Waals epitaxy. Nat. Commun. 2016, 7, 13352.

ACS Paragon Plus Environment

22

Page 23 of 23

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

The Journal of Physical Chemistry

TOC Graphic

ACS Paragon Plus Environment

23