Robust and Pristine Topological Dirac Semimetal Phase in Pressured

Sep 8, 2017 - Very recently, in spite of various efforts in searching for two-dimensional topological Dirac semimetals (2D TDSMs) in phosphorene, ther...
1 downloads 22 Views 1MB Size
Subscriber access provided by UNIVERSITY OF ADELAIDE LIBRARIES

Article

Robust and Pristine Topological Dirac Semimetal Phase in Pressured Two-Dimensional Black Phosphorous Penglai Gong, Bei Deng, Liang-Feng Huang, Liang Hu, Weichao Wang, Dayong Liu, Xingqiang Shi, Zhi Zeng, and Liang-Jian Zou J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.7b08926 • Publication Date (Web): 08 Sep 2017 Downloaded from http://pubs.acs.org on September 9, 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 17

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

Robust and Pristine Topological Dirac Semimetal Phase in Pressured Two-Dimensional Black Phosphorous Peng-Lai Gong,†,‡,⊥ Bei Deng,†,⊥ Liang-Feng Huang,¶ Liang Hu,† Wei-Chao Wang,‡ Da-Yong Liu,∗,§ Xing-Qiang Shi,∗,† Zhi Zeng,§,∥ and Liang-Jian Zou∗,§,∥ †Department of Physics, South University of Science and Technology of China, Shenzhen 518055, China ‡Department of Electronics and Tianjin Key Laboratory of Photo-Electronic Thin Film Device and Technology, Nankai University, Tianjin 300071, China ¶Department of Materials Science and Engineering, Northwestern University, Evanston, Illinois 60208, USA §Key Laboratory of Materials Physics, Institute of Solid State Physics, Chinese Academy of Sciences, P. O. Box 1129, Hefei 230031, China ∥University of Science and Technology of China, Hefei 230026, China ⊥Contributed equally to this work E-mail: [email protected]; [email protected]; [email protected]

1 ACS Paragon Plus Environment

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

Abstract

1

2

Very recently, in spite of various efforts in searching for two dimensional topological

3

Dirac semimetals (2D TDSMs) in phosphorene, there remains a lack of experimentally

4

efficient way to activate such phase transition and the underlying mechanism for the

5

topological phase acquisition is still controversial. Here, from first-principles calcula-

6

tions in combination with a band-sorting technique based on k·p theory, a layer-pressure

7

topological phase diagram is obtained and some of the controversies are clarified. We

8

demonstrate that, compared with tuning by external electric-fields, strain or doping by

9

adsorption, hydrostatic pressure can be an experimentally more feasible way to acti-

10

vate the topological phase transition for 2D TDSM acquisition in phosphorene. More

11

importantly, the resultant TDSM state is a pristine phase possessing a single pair of

12

symmetry-protected Dirac cones right at the Fermi level, in startling contrast to the

13

pressured bulk black phosphorous where only a carrier-mixed Dirac state can be ob-

14

tained. We corroborate that the Dirac points are robust under external perturbation as

15

long as the glide-plane symmetry preserves. Our findings provide a means to realize 2D

16

pristine TDSM in a more achievable manner, which could be crucial in the realization

17

of controllable TDSM states in phosphorene and related 2D materials.

18

Two dimensional (2D) Dirac semimetals, in which Dirac points cross the Fermi level

19

(EF ) being protected by nonsymmorphic crystal symmetries, was first proposed by Young

20

and Kane in 2015. 1 Very recent progress on the strain or electric-field modified phosphorene,

21

as a representative of the rare candidates, has established a link between Dirac cones and

22

the topological nature. 2–4 Given the experimental discovery of three dimensional (3D) topo-

23

logical Dirac semimetals (Na3 Bi, 5 Cd3 As2 6 ) and α-Sn on InSb(111) substrate, 7 it is natural

24

to ask whether the 2D topological Dirac semimetals (TDSMs) can also be realized in an

25

experimentally feasible manner.

26

A simple mechanism for a 2D TDSM phase acquisition has been proposed based on the

27

Stark effect in phosphorene thin films by applying an external electric field. 3,8 However,

28

applying an exceptionally giant electric field on such a system is difficult to realize exper2 ACS Paragon Plus Environment

Page 2 of 17

Page 3 of 17

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

29

imentally (the value of the field required is ∼0.5 V/Å for a four-layer phosphorene 8 ). On

30

the basis of the same mechanism, the experimentally observed Dirac semimetal state, from

31

potassium doping of few-layer phosphorene, is actually an electron-doped TDSM. 9 However,

32

a pristine TDSM in 2D is highly desirable as it can be tuned to be topological insulators

33

(TIs) or Weyl semimetals by explicit breaking of symmetries. In this regard, in-plane s-

34

train has been proposed as a possible means to induce Dirac cones in monolayer or bilayer

35

phosphorene, 10–12 whereas the critical strain (as large as ∼10% uniaxial strain or 5% biaxial

36

strain based on the DFT-PBE level 11 ) is difficult to be experimentally realized, particularly

37

for the biaxial strain in strong-anisotropic systems such as phosphorene. Furthermore, such

38

large uniaxial/biaxial strains could substantially lower the local symmetry of the system,

39

resulting in large distortion and structural instability due to the accumulative strain energy.

40

Additionally, it remains conflicted in literature about the true phase (i.e., whether the re-

41

sultant phase is a TI or a TDSM) in the strain-induced phase transition. 2,4,10 All the above

42

issues need to be solved or clarified, and examined by a more achievable way.

43

Apart from uniaxial/biaxial strains, the hydrostatic pressures, which to a larger extent

44

preserve the crystal symmetry, and also could be free of the typical epitaxial-mismatch

45

effect, have been proved as a powerful tool to deal with both fundamental and practical

46

issues. Very recently, we have reported that bulk black phosphorous (BP) can convert from

47

a normal insulator (NI) into a 3D Dirac semimetal under a certain hydrostatic pressure,

48

both experimentally and theoretically, but only a carrier-mixed phase is acquired. 13,14 In

49

experiments, hydrostatic pressures can also substantially modify the optical and vibrational

50

properties of 2D systems like monolayer or few-layer MoS2 , 15–20 WS2 , 21 and WSe2 . 22 Starting

51

from this point, the few-layer phosphorene could be feasibly pressurized like other layered

52

2D systems. Therefore, in this Letter, from first-principles methods in conjunction with

53

band-sorting technique based on k·p theory (to solve the above mentioned TI or TDSM

54

conflict), we study the layer-pressure topological phase diagram of few-layer phosphorene.

55

We discovered that hydrostatic pressures can energetically drive the 2D phosphorene from NI

3 ACS Paragon Plus Environment

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

56

to TDSM phase, with the topological-phase-transitions (TPTs) depending on the number of

57

layers. The resultant 2D TDSM state is characterized to be a pristine phase as a consequence

58

of the reduced dimensional effect in the vertical direction. We demonstrate that the Dirac

59

points are robust under external perturbation, including strain, pressure or electric field, as

60

long as the glide-plane symmetry within each sublayer preserves.

61

In order to explore the evolutions of electronic structures of few-layer phosphorene under

62

increasing hydrostatic pressures, we first determine their crystal structures (under each pres-

63

sure) by DFT calculations. We employed the Vienna Ab initio Simulation Package (VASP) 23

64

with the projector augmented wave (PAW) method. 24,25 The interlayer vdW interaction is

65

described by the optB88-vdW functional. 26 A previous study on few-layer MoS2 has pro-

66

posed an extraction method to obtain the 2D structures under pressures. 27 Based on this

67

method, their results reveal that the critical pressure for direct-to-indirect gap transition is

68

∼13 GPa, in qualitative agreement with the conclusion drawn from photoluminescence ex-

69

perimental measurement of ∼16 GPa. 15 The extraction method is a reasonable and practical

70

strategy to mimic the real pressurized procedure for layered systems, because it makes the

71

sample feel the pressure from the other parts of the bulk system, which could serve as a

72

good approximation to the inert pressure-transmitting medium in experiments. Meanwhile,

73

the periodic potential field at the vdW boundary is not such strong that can change the

74

geometric and electronic properties of the 2D BP system, as this is further verified by anoth-

75

er method, namely, the medium method, in which the very realistic pressure-transmitting

76

medium is considered. We find based on a set of comparison test, that basically the two

77

method will yield a same structure for the 2D BP under a given pressure (see Table S2),

78

suggesting that the extraction method is reliable to the cases in this study. As it is compu-

79

tationally much more expensive to consider the medium in each configuration (more than

80

800 atoms in the whole system), especially for the cases with the inclusion of spin-orbital

81

coupling (SOC) effect, here, we chose the extraction method as an appropriate approach to

82

conduct the study. Within the extraction method, in our study the structure of bulk BP was

4 ACS Paragon Plus Environment

Page 4 of 17

Page 5 of 17

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

Figure 1: (a) Side- and (b) top-views of the crystal structure of a 2-layer phosphorene; (c) 2D first Brillouin zone; (d) anisotropic variation of structural parameters of few-layer phosphorene under pressure. 83

first optimized under a specified hydrostatic pressure, and then these structural parameters

84

were used to construct the n−layer phosphorene, with a vacuum space of 20 Å along the z

85

direction. Based on the constructed 2D structures, the electronic properties is calculated by

86

hybrid functional of Heyd, Scuseria, and Ernzerhof (HSE06) 28 calculations with the inclusion

87

of SOC effect.

88

Few-layer BP is a layered material, in which each layer stacks together by van der Waals

89

(vdW) interactions. Each layer consists of two sublayers and in turn forms a bulking honey-

90

comb structure. The structural model of 2-layer phosphorene, as a representative for a finite

91

n-layer phosphorene (n denotes the number of layers), is shown in Figures 1a,b, where the

92

directions of lattice vectors a and b are along the zigzag (the y axis) and the armchair (the

93

x axis) directions, respectively. The 2D rectangular Brillouin zone of n-layer phosphorene is

94

typified in Figure 1c.

95

With increasing pressure, as shown in Figure 1d, the lattice constant b (armchair direc-

96

tion) and the interlayer spacing between two adjacent phosphorene layers (dint ) are signif-

97

icantly shortened (by 5.6% and 8.8% for b and dint at 4.0 GPa), while a (zigzag direction)

98

and the sublayer distance within one single layer (h) remain almost unchanged. The ther-

99

modynamic and lattice-dynamic stabilities of few-layer phosphorene depend on their atomic

100

and electronic properties, as well as the environment, such as pressure and temperature, etc.

101

In the current work, the stabilities of few-layer phosphorene are in line with that of bulk BP.

5 ACS Paragon Plus Environment

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

(b)

(a)

TDSM NI

Figure 2: Layer-pressure topological phase diagram of n-layer phosphorene. PC and PT are the critical pressures for electronic phase transition and upper limit of thermodynamic stability, respectively (see main text). NI denotes normal insulator (Figure 3a); and TDSM denotes a pure topological Dirac semimetal with Dirac points locating exactly at EF (Figures 3c,d). (b) Three dimensional band structure for the TDSM bands around Γ point, showing a single pair of Dirac cones. 102

The thermodynamic stability of the bulk BP is considered through the enthalpy-pressure re-

103

lationship. We find that the critical pressure of thermodynamic stability PT for bulk BP with

104

orthorhombic phase is 4.6 GPa (see Figure S2), very close to the experimental result (4.7

105

GPa), 29 which indicates a structural phase transition from a A17 (orthorhombic) phase to a

106

A7 (rhombohedral) phase. 30 For few-layer phosphorene, a similar structural phase transition

107

is found in accordance to our calculations because of the same mechanism of the structural

108

phase transition as the bulk BP. 30 Phonon spectra are then calculated to judge whether or

109

not the thermodynamically stable structures (with P ≤ PT ) are lattice-dynamic stable. Our

110

results show that all the pressured structures have no imaginary frequencies below PT . 14

111

The above results ensure that the pressured structures will not experience a structural phase

112

transition below PT ; and this is further confirmed by the medium method, showing that the

113

2D structure is indeed stable and does not favor reconstruction below PT (see Table S2 and

114

Figure S1).

115

Topological phase diagram. The layer-pressure topological phase diagram for few-layer

116

phosphorene (below PT , with consideration of SOC) is displayed in Figure 2, where PC and

117

PT are the critical pressures for electronic and structural phase transitions, respectively. At

118

PC , an electronic phase transition, from NI to TDSM, takes place (see Figure 3b). TDSM

6 ACS Paragon Plus Environment

Page 6 of 17

Page 7 of 17

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

+ -

+-

Figure 3: Band structures of 4-layer phosphorene as a representative for the pressured structures at (a) 0.0 GPa, (b) 4.1 GPa, (c) 4.4 GPa and (d) also 4.4 GPa, but around the Dirac point along the kx and ky directions, respectively. The SOC was also included in the calculation. The parity (even: +, odd: -) of the VBM and CBM at Γ point is labeled in (a) and (c). Dirac point (denoted by green point, labeled with "D") was set as the starting point in the band structure plotting in (d). A band-sorting method is used to sort the bands according to their symmetry (see text). 119

denotes a pristine topological Dirac semimetal phase with only Dirac fermions at EF (see

120

Figure 3c). No TPTs are observed for n = 2 and 3 layers so they are still NIs below PT

121

(dashed red line), this is because of their large band gaps (≥ 0.75 eV) under zero pressure.

122

With increasing number of n, their zero-pressure band gaps gradually decrease, and a NI

123

to TDSM transition happens to occur in the pressure range of PC