Electron-Transport Properties of Few-Layer Black Phosphorus - The

May 11, 2015 - Department of Chemistry and Nebraska Center for Materials and Nanoscience, University of Nebraska-Lincoln, 536 Hamilton Hall, Lincoln, ...
0 downloads 10 Views 476KB Size
Subscriber access provided by UB + Fachbibliothek Chemie | (FU-Bibliothekssystem)

Letter

Electron Transport properties of few-layer black phosphorus Yuehua Xu, Jun Dai, and Xiao Cheng Zeng J. Phys. Chem. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.jpclett.5b00510 • Publication Date (Web): 11 May 2015 Downloaded from http://pubs.acs.org on May 12, 2015

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 Letters 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 19

The Journal of Physical Chemistry Letters

1 2 3

Electron Transport Properties of Few-layer Black Phosphorus 4 5

Yuehua Xu1,2,3 , Jun Dai2 and Xiao Cheng Zeng2,* 7

6

School of Mathematics and Physics, Changzhou University, Changzhou 213164, People’s Republic of China 2. Department of Chemistry and Nebraska Center for Materials and Nanoscience, University of Nebraska-Lincoln, Lincoln, Nebraska 68588, United States 3 Shanghai Key Laboratory of Special Artificial Microstructure Materials and Technology, School of Physics Science and Engineering, Tongji University, Shanghai 200092, China *Email: [email protected]; 1.

18

17

16

15

14

13

12

1

10

9

8

19 20

Abstract 2

21

We perform the first-principles computational study of the effect of number of 24

23

stacking layers and stacking style of the few-layer black phosphorus (BPs) on the 26

25

electronic properties, including transport gap, current-voltage (i-v) relation and 28

27

differential conductance. Our computation is based on the non-equilibrium Green’s 30

29

function approach combined with density functional theory calculations. Specifically, 32

31

we compute electron transport properties of monolayer BP, bilayer BP and trilayer BP, 34

3

as well as bilayer BPs with either AB-, AA- or AC-stacking. We find that the stacking 36

35

number has greater influence to the transport gap than the stacking type. Conversely, 38

37

the stacking type has greater influence to i-v curve and differential conductance than 40

39

to the transport gap. This study offers useful guidance for determining the number of 42

41

stacking layers and the stacking style of few-layer BP sheets in future experimental 4

43

measurements and for potential applications in nanoelectronic devices. 46

45

TOC: 47 48 49 50 51 52 53 54 5 56 57 58 60

59

Keywords: interlayer interaction, two-probe configuration, transmission spectra, i-v curve,

differential conductance 1 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 4

3

The 6

5

recent successful exfoliation of another new two-dimensional atomic-layer

material, namely, layered black phosphorus (BP) or phosphorene from bulk black 8

7

phosphorus1-5 has attracted immediate attention owing to their remarkable electronic 10

9

properties1, 12

1

6-22

. Notably, the monolayer phosphorene possesses a desirable direct

bandgap of 1.5 eV18, 14

13

23

and significant transport anisotropy within the monolayer

plane, as well as the linear dichorism24-26, rendering it a potential candidate for future 16

15

nanoelectronic applications2, 6, 8, 14, 27-29. Few-layer BPs are predicted to have even 18

17

richer electronic properties and higher tunability compared to the monolayer BP23, 25, 20

19

30-34

2

21

such as the high carrier mobility up to 1000 cm2V-1S-1 and layer-dependent direct

bandgaps, ranging from 1.51 eV for monolayer BP to 0.59 eV for five-layer BP. The 24

23

reduction of bandgap in few-layer BPs is largely due to interlayer van der Waals 26

25

interaction24. The predicted bandgap reduction with increasing the number of BP 28

27

layers has been confirmed experimentally by Das et al35, evidenced by the 30

29

layer-dependent transport gap. In turn, determination of the exact layer number of 32

31

few-layer BP based on the measurement of the transport gap can be an effective tool 34

3

for separation of different few-layer BPs. In addition to the stacking number, our 36

35

recent calculation shows that the stacking type can also affect the bandgap of 38

37

few-layer BPs31. For bilayer BPs, three possible stacking types were considered, 40

39

namely AB-, AA- and AC-stacking, whose corresponding bandgaps are 1.04 eV, 0.95 42

eV, and 0.78 eV respectively31. Hence, it is expected that the transport gap of 4

43

41

few-layer BP can be also affected by the stacking type. It would be useful to examine 46

45

whether the transport gap of few-layer BP can be used as an indicator to distinguish 48

47

different stacking type, at least for bilayer BPs. 50

49

Besides the layer-dependent bandgaps, it is also important and timely to study the 52

51

effect of different stacking number and stacking type on electron transport properties, 54

53

such as current-voltage (i-v) curve and differential conductance. These properties are 56

5

closely relevant to the design of practical nanoelectronic devices. Moreover, 58

57

measurement of layer- and stacking-type-dependent i-v curve and differential 60

59

conductance can be also used as supplemental information, besides the measurement 2 ACS Paragon Plus Environment

Page 2 of 19

Page 3 of 19

The Journal of Physical Chemistry Letters

1 2 3

of transport gap, to distinguish the stacking number and stacking types of few-layer 4 5

BPs. 6 7

A sensible modeling of a practical electron-transport device much meet the 9

8

following requirements36-40: (1) Good contact with metallic electrodes; (2) three major 10 1

factors to determine the electronic transparency of contacts, i.e., favorable interface 12 13

geometry and bonding, the electronic density of states, and the potential barrier at the 14 15

interface; and (3) the top contact with electrode (because contacting a single layer 16 17

from side of electrodes is insufficient for good electron injection). We note that Guo 18 19

and coworkers investigated different metals as contact with monolayer BP for future 20 21

device application. They predicted that Cu (111) surface is the best candidate to form 2 23

an excellent ohmic contact with monolayer BP with a desirable modest binding 25

24

energy per BP unit-cell of 1.30 eV41. We therefore adopt the Cu (111) surface as the 26 27

metal contact to support few-layer BPs in a two-probe configuration for modeling. 28 29

The first-principles computation of transport properties is based on the 31

30

non-equilibrium Green’s function (NEGF) approach combined with the density 32 3

functional theory calculations. The monolayer, AB-stacking bilayer and ABA-stacking 34 35

trilayer BPs, as well as AB-, AA-, and AC-bilayer BPs are chosen as the model 36 37

systems to study the effect of stacking number and stacking type on the transport gap, 38 39

i-v curve, and differential conductance. The zigzag direction of few-layer BPs is 40 41

chosen as the transport direction because the electron mobility in the armchair 42 43

direction is expected to decrease much rapidly due to the strained effect. In the zigzag 45

4

direction, however, the electron mobility can be even enhanced by the strain25. 46 47

Besides the transport properties, our studies suggest possible strategies and 48 49

guiding rules for effective way to identify different stacking number and stacking type 50 51

of few-layer BPs. Specifically, we find that although both stacking number and 52 53

stacking type can affect the transport gap, i-v curve, and differential conductance of 54 5

few-layer BPs, the stacking number incurs greater influence than the stacking type, 56 57

while the stacking type has greater influence to the i-v curve and differential 58 60

59

conductance than to the transport gap, especially for the AB- and AA-stacking bilayer 3 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3

BPs. Hence, the transport gap measured from the transmission spectra is best used to 4 5

identify the stacking number of BPs, but not the stacking type of BPs. 6 7 8 9

Computation of electronic structures and transport properties is carried out using 1

10

the Atomistix Toolkit (ATK) code package based on the non-equilibrium Green’s 12 13

function (NEGF) approach combined with density functional theory (DFT) 15

14

calcualtions42-43. 16 17

To model the two-probe configuration, a supercell is constructed from a slab of 18 19

three atomic layers of Cu with a few-layer BP sheet adsorbed on the Cu (111) surface. 20 21

A vacuum region of at least 20 Å is set in the supercell. Fig. 1 shows the two-probe 2 23

configuration with a monolayer BP. The system is divided into three regions: left 24 25

electrode, right electrode and scattering region. Here, the transport direction c(z) is set 26 27

along the zigzag direction of the BPs only, and the periodic boundary conditions are 28 29

applied in the transverse transport direction, i.e., the system extends infinitely along 30 31

the b(y) direction. The primitive cell of few-layer BP is taken from our previous 3

32

studies30-31. The primitive cell of Cu has lattice constants of b = 4.496 Å and c = 2.596 35

34

Å44. The lattice constants of the electrodes are chosen to be those of (1×1×3) strained 36 37

cell of few-layer BP as b = 4.496 Å and c = 10.382 Å, the same as those for a (1×1×4) 38 39

cell of Cu(111). The lattice-constant matching is reasonable since the mismatch from 40 41

the lattice parameters of Cu (111) surface is only about 1.8% in the b direction and 42 43

4% in the c direction. The effect of strained cell of few-layer BP on the electron 45

4

mobility was discussed in detail in a previous study25. It was shown that when the 46 47

few-layer BP being under a biaxial strain of 6%, the electron mobility along the 49

48

zigzag direction is even higher than unstrained ones25. So it is believed that strained 50 51

few-layer BP still possesses good electron transport properties in the zigzag direction. 52 53

A short channel device is considered with the length of the scattering region in the c(z) 54 5

direction being about 62 Å (see the two inner vertical white lines in Fig. 1 (a)) , 56 57

making the length of few-layer BP channel about 36 Å (see the red lines in Fig. 1(a)). 58 60

59

This length is long enough to avoid the effect of metal-induced gap states for the 4 ACS Paragon Plus Environment

Page 4 of 19

Page 5 of 19

The Journal of Physical Chemistry Letters

1 2 3

few-layer BP (see Fig. 1(b)). 4 5 6 7

(a) 8 9 10 1 12 13 14 15 16 17 18 19 20 21 2 23 24 26

25

(b) 27 28

60

29

PDOS

30

38

37

36

35

34

3

32

PDOS (arb.units)

31 40

20

40

39 0 -1.0

42

41

-0.5

43

0.0 E (eV)

0.5

1.0

4 46

45

Fig. 1 (a) The optimized structure of two-probe configuration with monolayer BP (the 48

47

lower panel is a top view of the two-probe configuration). The light purple balls 50

49

denote phosphorus atoms and the gold balls denote Cu atoms. The length of the 52

51

scattering region (between two inner vertical white lines) in the c(z) direction is about 54

53

62 Å, while the length of the free-standing few-layer BP channel (between two 56

5

vertical red dashed line) is about 38 Å. (b) the Partial density of states (PDOS) of 58

57

phosphorus atoms in the channel. 59 60

5 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3

In the DFT calculation of electronic structures, we use a numerical basis set to 4 5

expand the wave function, i.e., the double zeta plus polarization basis set (DZP). The 6 7

generalized gradient approximation (GGA) is used for the electron exchange and 8 9

correlation and the Troullier-Martins pseudopotentials are used for the atomic cores. 10 1

The density mesh cut-off is chosen 150 Ry to achieve a balance between calculation 12 13

efficiency and accuracy. The k-point grids of 1 × 8 × 100 and 1 × 32 × 100 are used to 14 15

sample the Brillouin zone of the electrodes in the x, y and z direction for the geometry 16 17

relaxation and for computation of total energy, transmission spectra, and i-v curve, 18 19

respectively. The atomic structures of the junctions are fully optimized by minimizing 20 21

the atomic forces on the atoms to be smaller than 0.05 eV/Å. The separation between 2 23

the two electrodes is also optimized by minimizing the stress in the direction of the 24 25

junctions. 26 27 28 29

Upon DFT optimization, the average distance dz between the topmost layer of Cu 30 31

(111) surface and the bottom layer of few-layer BP is 2.27 Å, 2.26 Å and 2.26 Å for 32 3

the monolayer, AB-stacking bilayer and ABA-staking trilayer, respectively. These 34 35

values are slightly less than 2.31 Å for the monolayer BP obtained by Guo and 37

36

coworkers.41 The shortest bond length dm between a topmost atom on Cu (111) and 38 39

the bottom-layer of BP is 2.34 Å for monolayer BP, 2.35 Å for AB-stacking bilayer 40 41

BP, and 2.35 Å for ABA-stacking trilayer BP, very close to 2.36 Å for the monolayer 43

42

BP, obtained by Guo and workers.41 For the AA- and AC-stacking bilayer BPs, our 4 45

calculated dz is 2.23 Å and 2.26 Å, respectively, while dm is 2.36 Å and 2.35 Å, 46 47

respectively. 48 49

Next, we present results of computed transport gap, zigzag directional i-v curve, 50 51

normalized current (current per layer) and differential conductance in the two-probe 52 53

configuration, i.e., few-layer BP-Cu (111) junctions. First, we examine whether the 54 5

contact between BP and Cu is Omhic. Our calculation shows that the difference in the 56 57

work function b is -0.30 eV between Cu (111) and monolayer BP, which is quite 60

59

58

close to the value -0.36 eV obtained by Guo and coworker41. In addition, we compute 6 ACS Paragon Plus Environment

Page 6 of 19

Page 7 of 19

The Journal of Physical Chemistry Letters

1 3

2

the electronic potential difference V between the lower P atoms and the Cu 4 5

substrate in the contact, which is -0.50 eV for monolayer BP. This result indicates that 6 7

the electrons can be easily injected from Cu to the BP without any energy barrier. 8 9

Hence, the contact between Cu (111) and monolayer BP is a good Ohmic contact. In 10 1

addition, we also check whether the electrons can be easily injected from the 12 13

combined BP-Cu contact to the freestanding part of BP channel (see Fig. 1(a)). 14 15

Computed PDOS of the freestanding part of BP channel can be used to estimate the 16 17

Schottky barrier, which is about 0.36 eV (see Fig. S1). Fig. 2(a) shows the 18 19

transmission spectra of monolayer BP, AB-stacking bilayer BP, and ABA-stacking 20 21

trilayer BP under zero bias. Note that the top layer in bilayer BP and the top two 2 23

layers in trilayer BP are electronically decoupled from the Cu (111) substrate. So their 24 25

contribution to the transmission spectra and to the i-v curves (see below) can be 26 27

deduced from subtraction (see Fig. 2(b)). 28 29

In Fig. 2(a), the transport gap decreases with increasing the layer number. More 30 31

specifically, the transport gap is 1.15 eV for monolayer BP, 0.71 eV for the bilayer BP, 32 3

and 0.63 eV for the trilayer BP. According to the DFT band-structure calculation, the 34 35

bandgap of bare monolayer, bilayer, and trilayer BP (all under a strain due to 36 37

lattice-constant matching) is 1.17 eV, 0.78 eV and 0.67 eV, respectively (see Fig. S2). 38 39

Hence, the computed transport gaps are largely consistent with computed bandgaps of 40 41

free-standing few-layer BP. Notably, these computed transport gaps are nearly the 43

42

same as corresponding experimental values35, i.e., 1.0 eV for monolayer BP, 0.73 eV 4 45

for bilayer BP, and 0.60 eV for trilayer BP. Note also that the reduction of bandgap is 46 47

not linearly correlated with increasing the stacking number, suggesting that the 48 49

transport gap would become less effective to identify the stacking number for 50 51

many-layer BPs. 52 53 54 5 56 57 58 59 60

7 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

Page 8 of 19

1 2 3 4

9

5

8

6

7

7

5

(a)

monolayer BP bilayer BP trilayer BP

4

(b) the 2nd layer of bilayer BP the 2nd and 3rd layers of trilayer BP

6

1

3

5

T(E)

10

T(E)

9

8

4

2

12

3 2

14

1

15

13

0 -1.0

17

16

1

-0.5

0.0 energy (eV)

0.5

1.0

0 -1.0

-0.5

0.0 energy (eV)

0.5

1.0

18 20

19

Fig. 2 (a) Computed zigzag directional transmission spectra of monolayer BP, 2

21

AB-stacking bilayer BP and ABA-stacking trilayer BP under zero bias, and (b) 24

23

deduced contribution to the corresponding transmission spectra from the top layer of 26

25

bilayer BP and top two layers of trilayer BP under zero bias. 27 28 30

29

In Fig. 3, the computed zigzag directional transmission spectra of AB-, AA-, and 32

31

AC-stacking bilayer BPs, as well as the contribution by the top layer to their 34

3

transmission spectra are shown, where the corresponding transport gap is 0.71 eV, 36

35

0.79 eV, and 0.57 eV, respectively. The transport gaps are also close to the computed 38

37

bandgaps of bare AB-, AA, and AC-stacking bilayer BPs (under the same strain), 40

39

which are 0.78 eV, 0.79 eV and 0.57 eV, respectively (see Fig. S3). The AC-stacking 42

41

bilayer gives notably smaller bandgap than the AB- and AA-stacking bilayers, while 4

43

the AA-stacking and AB-stacking bilayer BPs give nearly the same bandgap. So for 46

45

bilayer BPs, it is easier to identify the AC-stacking based on the transport gap, but not 48

47

so to differentiate AB- and AA-stacking. 49 50 51 52 53 54 5 56 57 58 59 60

8 ACS Paragon Plus Environment

Page 9 of 19

The Journal of Physical Chemistry Letters

1 2 3

7

4 5

6

6 7

4

(a)

AB-stacking AA-stacking AC-stacking

3

5

12

1

4

T(E)

10

T(E)

9

8

(b)

2

3 2

13

the 2nd layer of AB-stacking BP the 2nd layer of AA-stacking BP the 2nd layer of AC-stacking BP

1

15

14

1 0 -1.0

17

16

-0.5

0.0 energy (eV)

0.5

1.0

0 -1.0

-0.5

0.0 0.5 energy (eV)

1.0

18 20

19

Fig. 3 (a) Computed zigzag directional transmission spectra of bilayer BP with AA-, 2

21

AB-, and AC-stacking, under zero bias, and (b) the contribution by the top layer to 24

23

their transmission spectra. 26

25

Fig. 4 displays the computed zigzag directional i-v curves (based on the 27

Landauer-Buttiker formula45) for the monolayer BP, AB-stacking bilayer BP, and 30

29

28

ABA-stacking trilayer BP under the bias ranging from 0 to 2 V. It can be seen that the 32

31

amplitude of the current increases with increasing the stacking number. For example, 34

3

under a bias of 2 V the current is 249 nA, 688 nA, and 936 nA for monolayer BP, 36

35

bilayer BP, and trilayer BP, respectively. This behavior can be understood because the 38

37

more stacking number for the few-layer BP, the more transmission channels would 40

39

form. Besides the generic trend of i-v curves for all three cases, some special features 42

41

can still be seen in each i-v curve. Specifically, for monolayer BP, a small decrease in 4

43

current exists between bias 1.7 V and 1.8 V, and for bilayer BP, a similar decrease 46

45

also exists between 1.6 V and 1.7 V. While for trilayer BP, no such decrease in 48

47

current is observed in the bias range from 0 to 2 V. In addition, the deduced 50

49

contribution from the top layer of bilayer BP and top two layer of trilayer BP to the 52

51

current is shown in Fig. 4(b). Fig. 4(c) displays the normalized current (current per 54

53

layer) for monolayer BP, bilayer BP and trilayer BP, respectively. It can be also seen 56

5

that several special features in Fig. 4(b) and Fig. 4(c) also arise in their corresponding 58

57

i-v curve. Hence, one may determine existence of monolayer, bilayer or trilayer BP in 60

59

the sample, based on the amplitude of current as well as the presence and position of 9 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

Page 10 of 19

1 2 3

the decrease of current in i-v curve. 4 5 6 1000

9 10

800

16

15

14

13

Current (nA)

12

1

800

(a)

monolayer BP bilayer BP trilayer BP

600 400

400

200

200

17

(b) the 2nd layer of bilayer BP the 2nd and 3rd layers of trilayer BP

600

Current (nA)

8

7

18 0 0.0

20 21

400

30

29

28

27

26

25

24

Current per layer (nA)

23

2

300

0.5

1.0 Voltage (Volt)

1.5

0 0.0

2.0

1600

(c)

monolayer BP bilayer BP trilayer BP

1200

dI/dV (nS)

19

200

100

0.5

1.0 Voltage (Volt)

1.5

2.0

(d)

monolayer BP bilayer BP trilayer BP

800

400

32

31 0 0.0

35

34

3

0.5

1.0 Voltage (Volt)

1.5

2.0

0 0.0

0.5

1.0 Voltage (Volt)

1.5

1.9

Fig. 4 Computed zigzag directional (a) i-v curves, (c) normalized current and (d) 37

36

differential-conductance curve of the monolayer BP, bilayer BP, and trilayer BP. (b) 39

38

The contribution from the top layer of bilayer BP and top two layers of trilayer BP to 41

40

the corresponding i-v curves. 42 43 45

4

Based on the computed i-v curve (with electron transport being in the zigzag 47

46

direction), the differential conductance di/dv under bias 0 - 1.9 V can be computed 49

48

(see Fig. 4(d)). One can see that the differences among the differential conductance 51

50

are more pronounced compared to the i-v curve. Specifically, for the monolayer BP, a 52 53

valley in the differential conductance can be seen under bias ~1.5 V, whereas for the 54 5

bilayer BP and trilayer BP, a peak, rather than a valley, can be seen within this bias 56 57

region. Moreover, for the monolayer BP, there exists a peak under bias ~1.6 V, while 58 60

59

for the bilayer BP and trilayer BP, there is a valley near this bias region. Hence, monolayer BP can be distinguished from the bilayer and trilayer BPs according to the 10 ACS Paragon Plus Environment

Page 11 of 19

The Journal of Physical Chemistry Letters

1 2 3

characteristic valley in the differential conductance under bias ~1.5 V as well as the 4 5

peak under about 1.6 V in i-v curve. In addition, for the trilayer BP, there is a peak 6 7

under bias ~1.1 V while for the bilayer BP, no such characteristic peak is seen. 8 9

Furthermore, for the trilayer BP, the highest peak is located at bias ~1.8 V, while for 10 1

the bilayer BP, although a similar peak exists at this bias, the peak height is much 12 13

lower than that of the trilayer BP. Hence, the characteristic peaks at 1.1 V and 1.8 V 14 15

in the differential-conductance curve of trilayer BP are important features to 16 17

differentiate it from bilayer BP. 18 19

We have also studied the effects of different stacking type on i-v curve using the 20 21

two-probe configuration. The computed i-v curves for AB-, AA- and AC-bilayer 2 23

BP-Cu (111) junctions are shown in Fig. 5(a). In general, values of current under bias 24 25

0 - 2 V are close to one another among the three stacking types, although some 26 27

distinctive features can still be discerned. For the AB-stacking, as mentioned above, a 28 29

notable decrease in current exists between bias 1.6 V and 1.7 V, while the current 30 31

exhibits monotonic increase in this bias region for the AA-stacking. For the 32 3

AC-stacking, the maximum current (a peak) exhibits at about 1.9 V while for the 34 35

AB-stacking, the maximum is located under 2.0 V bias. Thus, although the amplitude 36 37

of current is more or less the same for the AB-, AA-, and AC-stacking bilayer BPs, 38 39

the decrease in current between 1.6 V and 1.7 V for AB-stacking and the decrease in 40 41

current between 1.9 V and 2.0 V for AC-stacking can be used as clear indicators to 42 43

differentiate the two stacking types from the AA-stacking type. In addition, the 4 45

deduced contribution from the top layer of AB-, AA-, and AC-stacking bilayer BPs to 46 47

the current is shown in Fig. 5(b). The normalized current is also calculated for the 48 49

AB-, AA-, and AC-stacking bilayer BPs (Fig. 5(c)). It can be seen that several main 50 51

features shown in Fig. 5(b) and 5(c) are also seen in the corresponding i-v curves (Fig. 52 53

5(a)). 54 5 56 57 58 59 60

11 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

Page 12 of 19

1 2 3 800

5 6 600

400

9

Current (nA)

8

7

400

12

1

10

200

13

500

(a)

AB-stacking AA-stacking AC-stacking

Current (nA)

4

the 2nd layer of AB-stacking BP the 2nd layer of AA-stacking BP the 2nd layer of AC-stacking BP

(b)

300 200 100

14 0 0.0

16 17

400

18

28

27

26

25

24

23

2

21

Current per layer (nA)

20

19 300

31

30

1.0 Voltage (Volt)

1.5

0 0.0

2.0

1200

(c)

AB-stacking AA-stacking AC-stacking

0.5

1.0 Voltage (Volt)

1.5

2.0 (d)

AB-stacking AA-stacking AC-stacking

800 200

400

100

0 0.0

29

0.5

dI/dV (nS)

15

0.5

1.0 Voltage (Volt)

1.5

2.0

0 0.0

0.5

1.0 Voltage (Volt)

1.5

1.9

Fig. 5 Computed zigzag directional (a) i-v curve, (c) normalized current (d) 3

32

differential conductance of bilayer BPs with AB-stacking, AA-stacking and 35

34

AC-stacking, (b) the deduced contribution from the top layer of AB-, AA-, and 37

36

AC-stacking bilayer BP to the current. 38 39 41

40

The computed differential-conductance curves for AB-, AA-, and AC-stacking 43

42

bilayer BPs under bias from 0 to 1.9 V is shown in Fig. 5(c). Again, the 45

4

differential-conductance curves exhibit more striking differences compared to the 47

46

corresponding i-v curves. Specifically, for the AA-stacking, a clear valley is seen 49

48

under ~1.5 V bias, whereas for both the AB- and AC-stacking, a peak is seen near this 51

50

bias. Moreover, for both AB- and AC-stacking, a deep valley is seen under ~1.6 V 53

52

and ~1.7 V bias, respectively. Therefore, these valleys under 1.5 V, 1.6 V or 1.7 V 5

54

bias can be used to distinguish the three stacking types from one another. 56 57 58 59 60

12 ACS Paragon Plus Environment

Page 13 of 19

The Journal of Physical Chemistry Letters

1 2 3 0.25

4 5

(a)

0.20

6

0.6

Bias 1.7 V Bias 1.8 V

7 8

0.4

T(E) 1

10 12

T(E)

0.15

9

0.10

14 16 17

0.20

18 19 20

1.8 V -0.5

T(E) 25

24

0.5

1.0

0.0 -1.0

Bias 1.9 V Bias 2.0 V

(c)

0.15

23

0.0

energy (eV)

1.9 V

21 2

0.2

0.10 2.0 V 0.05

26 27 28

0.00

-1.0

30

29

1.6 V

0.1

0.00 -1.0

15

0.3

1.7 V

0.05

13

Bias 1.6 V Bias 1.7 V

(b)

0.5

-0.5

0.0

energy (eV)

0.5

1.0

1.7 V -0.5

0.0

energy (eV)

0.5

1.0

Fig.6 Computed zigzag directional transmission spectra of (a) monolayer BP under bias 1.7 V and 1.8 V; (b) AB-stacking bilayer BP under bias 1.6 V and 1.7 V; and (c) the AC-stacking bilayer BP under 1.9 V and 2.0 V. To highlight the integration of transmission spectra through the bias window, the integral areas are filled with different color.

31 3

32

We now discuss why these characteristic features arise for BPs with different 35

34

layer number or different stacking type. First, for the monolayer BP, the transmission 37

36

spectra under bias 1.7 V and 1.8 V are shown in Fig. 6(a). Clearly, under the bias 1.7 39

38

V, the integration of the transmission spectrum over the bias range [-0.85 V, 0.85 V] 41

40

is greater than that over the range [-0.9 V, 0.9 V]. Thus, the current is expected to 43

42

decrease within the bias range of 1.7 to 1.8 V. Likewise, for the bilayer BP, the 45

4

decrease in current between the bias 1.6 V and 1.7 V can be understood from 47

46

analyzing the corresponding transmission spectrum under bias 1.6 V and 1.7 V (see 49

48

Fig. 6(b)). Clearly, integration of the transmission spectrum over the bias range [-0.8 51

50

V, 0.8 V] is greater than that over the range [-0.85 V, 0.85 V]. In addition, for the 53

52

AC-stacking bilayer, again, the characteristic structures in i-v curve can be understood 5

54

from the corresponding transmission spectrum. The decrease in current between bias 57

56

1.9 V and 2.0 V, the corresponding spectrum under 1.9 V and 2.0 V is shown in Fig. 60

59

58

6(c). Integration of the transmission spectrum over the bias range [-0.95 V, 0.95 V] is greater than that over the range [-1.0 V, 1.0V]. 13 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3

Second, characteristic features in the differential-conductance curves of few-layer 4 5

BPs can be understood through their corresponding i-v curves. For the monolayer BP, 6 7

the valley at ~1.5 V bias stems from the nearly unchanged current at 1.5 V bias (see 8 9

Fig. 4(a)). The peak at ~1.6 V bias can be understood from the trend of rapid increase 10 1

of current between bias 1.6 V and 1.7 V (Fig. 4(a)). For the bilayer BP, the valley at 12 13

~1.6 V bias can be understood from the decreasing current between 1.6 V and 1.7 V 14 15

(Fig. 4(a)). For the trilayer BP, the peaks under ~1.2 V can also be understood from 16 17

the fast increase of current from bias 0.9 V to 1.1 V. Similarly, the peak at ~1.8 V bias 18 19

can be understood from the even faster increase of current than that for bilayer BP. 20 21

Thus the amplitude of peak at 1.8 V bias is greater than that for bilayer BP. 2 23

Lastly, the characteristic features in the differential conductance of AB-, AA- and 24 25

AC-stacking bilayer BPs can be understood from their corresponding i-v curve. For 26 27

the AA-stacking, the valley at ~1.5 V bias is resulted from the slower increase of 28 29

current in this bias region compared to the current in the neighboring region. The 30 31

peaks at 1.5 V in both the AB- and AC-stacking cases are due to the faster increase of 32 3

current in this bias region compared to the current in the neighboring region. 34 35 36 37

We have performed, to our knowledge, the first systematic study of the effect of 38 39

stacking number and stacking type of few-layer BPs on the transport gap, zigzag 40 41

directional i-v curve, and differential conductance. We compute electron transport 42 43

properties of monolayer BP, bilayer BP and trilayer BP, as well as bilayer BPs with 4 45

either AB-, AA- or AC-stacking. We find that the stacking number has greater 46 47

influence to the transport gap than the stacking type. The characteristic features in i-v 48 49

curves of monolayer BP, bilayer BP, and trilayer BP are more pronounced, 50 51

particularly in their relative magnitude of current, than those of AA-, AB-, and 52 53

AC-stacking bilayer BPs. Moreover, the differential conductance deduced from i-v 54 5

curve is much more sensitive to the stacking number than to the stacking type. On the 56 57

other hand, the stacking type has greater influence to i-v curve and differential 58 60

59

conductance than to the transport gap. The difference in the transport gap between 14 ACS Paragon Plus Environment

Page 14 of 19

Page 15 of 19

The Journal of Physical Chemistry Letters

1 2 3

AB- and AA-stacking bilayer BPs is merely 0.08 eV. A notable characteristic feature 4 5

in i-v curve is the decrease in the current between 1.6 V and 1.7 V for AB-stacking. 6 7

The valleys at 1.5 V and 1.6 V bias in the differential-conductance curves of AA- and 8 9

AB-stacking BPs are even more striking. These features can be useful indicators to 10 1

identify the AB-stacking from AA-stacking for bilayer BPs. 12 13

More importantly, our results may offer helpful guidance on how to identify 14 15

stacking number and stacking style of few-layer BP sheets for future experimental 16 17

measurements and applications in nanoelectronic devices. For example, for the 18 19

identification of the stacking number of 1 to 3, the transport gap in the transmission 20 21

spectrum appears to be an effective way, but not so effective for many-layer BPs. To 2 23

identify the AA- and AB-stacking for bilayer BPs, a combination of transport-gap 24 25

data along with zigzag directional i-v curve, and differential-conductance curve is a 26 27

much more effective approach. 28 29 30 31

ACKNOWLEDGMENTS. 32 3

We thank Dr. Anders Blom for valuable discussions. The work was supported by the 34 35

National Natural Science Foundation of China under Grant NO.11404038, the Special 36 37

Foundation for theoretical physics Research Program of China under Grant NO. 38 39

11347149, the Natural Science Foundation of the Jiangsu Higher Education 40 41

Institutions of China under Grant No. 13KJB140001, Jiangsu Government 42 43

Scholarship for Overseas Studies, the High Performance Computing Laboratory of 4 45

Changzhou University. XCZ was supported by the National Science Foundation (NSF) 46 47

through the Nebraska Materials Research Science and Engineering Center (MRSEC) 48 49

(grant No. DMR-1420645). 50 51 53

52

Supporting Information Available: Description of the material included. This material is available free of charge via the Internet 56

5

54

http://pubs.acs.org. 57 58 59 60

15 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

Page 16 of 19

1 2 3

REFERENCES 5

4

((1) Liu, H.; Neal, A. T.; Zhu, Z.; Luo, Z.; Xu, X.; Tománek, D.; Ye, P. D. Phosphorene: An

6

Unexplored 2d Semiconductor with a High Hole Mobility. ACS Nano 2014, 8, 4033-4041.

8

(2) Li, L.; Yu, Y.; Ye, G. J.; Ge, Q.; Ou, X.; Wu, H.; Feng, D.; Chen, X. H.; Zhang, Y. Black

9

7

Phosphorus Field-Effect Transistors. Nat Nano 2014, 9, 372-377.

10

(3) Tromel, M. Atomvolumen, Packungsdichte Der Atome Und Chemische Bindung in

12

Nichtmetallischen Elementen. Acta Crystallographica Section B 2007, 63, 532-536.

13

1

(4) Koenig, S. P.; Doganov, R. A.; Schmidt, H.; Castro Neto, A. H.; Özyilmaz, B. Electric Field Effect

14

in Ultrathin Black Phosphorus. Appl. Phys. Lett. 2014, 104, 103106-1-103106-5.

16

(5) Buscema, M.; Groenendijk, D. J.; Blanter, S. I.; Steele, G. A.; van der Zant, H. S. J.;

17

15

Castellanos-Gomez, A. Fast and Broadband Photoresponse of Few-Layer Black Phosphorus

18

Field-Effect Transistors. Nano Lett. 2014, 14, 3347-3352.

20

(6) Liu, Q.; Zhang, X.; Abdalla, L. B.; Fazzio, A.; Zunger, A. Switching a Normal Insulator into a

21

19

Topological Insulator Via Electric Field with Application to Phosphorene. Nano Lett. 2015, 15,

23

1222-1228.

24

(7) Wu, M.; Fu, H.; Zhou, L.; Yao, K.; Zeng, X. C. Nine New Phosphorene Polymorphs with

25

2

Non-Honeycomb

27

DOI: 10.1021/acs.nanolett.5b01041

28

(8) Cai, Y.; Ke, Q.; Zhang, G.; Zhang, Y.-W. Energetics, Charge Transfer, and Magnetism of Small

29

26

Structures:

A

Much

Extended

Family.

Nano

Lett.

2015,

Molecules Physisorbed on Phosphorene. J. Phys. Chem. C. 2015, 119, 3102-3110.

31

(9) Das, S.; Demarteau, M.; Roelofs, A. Ambipolar Phosphorene Field Effect Transistor. ACS Nano

32

2014, 8, 11730-11738.

3

30

(10) Fei, R.; Faghaninia, A.; Soklaski, R.; Yan, J.-A.; Lo, C.; Yang, L. Enhanced Thermoelectric

35

Efficiency Via Orthogonal Electrical and Thermal Conductances in Phosphorene. Nano Lett. 2014, 14,

36

6393-6399.

37

34

(11) Guan, J.; Zhu, Z.; Tománek, D. Tiling Phosphorene. ACS Nano 2014, 8, 12763-12768.

39

(12) Hashmi, A.; Hong, J. Transition Metal Doped Phosphorene: First-Principles Study. J. Phys. Chem.

40

C. 2015,119,9198-9204

41

38

(13) Hu, T.; Hong, J. First-Principles Study of Metal Adatom Adsorption on Black Phosphorene. J.

43

Phys. Chem. C. 2015, 119, 8199-8207.

4

(14) Kou, L.; Frauenheim, T.; Chen, C. Phosphorene as a Superior Gas Sensor: Selective Adsorption

45

42

and Distinct I–V Response. J.Phys.Chem. Lett. 2014, 5, 2675-2681.

47

(15) Kou, L.; Ma, Y.; Smith, S. C.; Chen, C. Anisotropic Ripple Deformation in Phosphorene.

48

J.Phys.Chem. Lett. 2015, 6, 1509-1513.

49

46

(16) Li, W.; Yang, Y.; Zhang, G.; Zhang, Y.-W. Ultrafast and Directional Diffusion of Lithium in

51

Phosphorene for High-Performance Lithium-Ion Battery. Nano Lett. 2015, 15, 1691-1697.

52

(17) Li, W.; Zhang, G.; Zhang, Y.-W. Electronic Properties of Edge-Hydrogenated Phosphorene

53

50

Nanoribbons: A First-Principles Study. J. Phys. Chem. C. 2014, 118, 22368-22372.

5

(18) Liang, L.; Wang, J.; Lin, W.; Sumpter, B. G.; Meunier, V.; Pan, M. Electronic Bandgap and Edge

56

Reconstruction in Phosphorene Materials. Nano Lett. 2014, 14, 6400-6406.

57

54

58

(19) Srivastava, P.; Hembram, K. P. S. S.; Mizuseki, H.; Lee, K.-R.; Han, S. S.; Kim, S. Tuning the

59

Electronic and Magnetic Properties of Phosphorene by Vacancies and Adatoms. J. Phys. Chem. C. 2015,

60

119, 6530-6538. (20) Sui, X.; Si, C.; Shao, B.; Zou, X.; Wu, J.; Gu, B.-L.; Duan, W. Tunable Magnetism in Transition

16 ACS Paragon Plus Environment

Page 17 of 19

The Journal of Physical Chemistry Letters

1 3

2 Metal Decorated Phosphorene. J. Phys. Chem. C. 2015,119,10059-10063

4

(21) Zhang, R.; Li, B.; Yang, J. A First-Principles Study on Electron Donor and Acceptor Molecules

6

Adsorbed on Phosphorene. J. Phys. Chem. C. 2015, 119, 2871-2878.

7

5

(22) Sa, B.; Li, Y.-L.; Qi, J.; Ahuja, R.; Sun, Z. Strain Engineering for Phosphorene: The Potential

8

Application as a Photocatalyst. J. Phys. Chem. C. 2014, 118, 26560-26568.

10

(23) Tran, V.; Soklaski, R.; Liang, Y.; Yang, L. Layer-Controlled Band Gap and Anisotropic Excitons

1

9

in Few-Layer Black Phosphorus. Phys. Rev. B 2014, 89, 235319-1-235319-6.

12

(24) Qiao, J.; Kong, X.; Hu, Z.-X.; Yang, F.; Ji, W. High-Mobility Transport Anisotropy and Linear

14

Dichroism in Few-Layer Black Phosphorus. Nat Commun 2014, 5, 4475.

15

13

(25) Fei, R.; Yang, L. Strain-Engineering the Anisotropic Electrical Conductance of Few-Layer Black

16

Phosphorus. Nano Lett. 2014, 14, 2884-2889.

18

(26) Rodin, A. S.; Carvalho, A.; Castro Neto, A. H. Strain-Induced Gap Modification in Black

19

17

Phosphorus. Phys. Rev. Lett. 2014, 112, 176801-1-176801-5.

21

(27) Han, X.; Morgan Stewart, H.; Shevlin, S. A.; Catlow, C. R. A.; Guo, Z. X. Strain and Orientation

2

Modulated Bandgaps and Effective Masses of Phosphorene Nanoribbons. Nano Lett. 2014, 14,

23

20

4607-4614.

25

(28) Li, Y.; Yang, S.; Li, J. Modulation of the Electronic Properties of Ultrathin Black Phosphorus by

26

Strain and Electrical Field. J. Phys. Chem. C. 2014, 118, 23970-23976.

27

24

(29) Yao, Q.; Huang, C.; Yuan, Y.; Liu, Y.; Liu, S.; Deng, K.; Kan, E. Theoretical Prediction of

29

Phosphorene and Nanoribbons as Fast-Charging Li Ion Battery Anode Materials. J. Phys. Chem. C.

30

2015, 119, 6923-6928.

31

28

(30) Guo, H.; Lu, N.; Dai, J.; Wu, X.; Zeng, X. C. Phosphorene Nanoribbons, Phosphorus Nanotubes,

3

and Van Der Waals Multilayers. J. Phys. Chem. C. 2014, 118, 14051-14059.

34

(31) Dai, J.; Zeng, X. C. Bilayer Phosphorene: Effect of Stacking Order on Bandgap and Its Potential

35

32

Applications in Thin-Film Solar Cells. J.Phys.Chem. Lett. 2014, 5, 1289-1293.

37

(32) Wang, V.; Kawazoe, Y.; Geng, W. T. Native Point Defects in Few-Layer Phosphorene. Phys. Rev.

38

B 2015, 91, 045433-1-045433-9.

39

36

(33) Guan, J.; Zhu, Z.; Tománek, D. Phase Coexistence and Metal-Insulator Transition in Few-Layer

41

Phosphorene: A Computational Study. Phys. Rev. Lett. 2014, 113, 046804-1-046804-5.

42

(34) Low, T.; Roldán, R.; Wang, H.; Xia, F.; Avouris, P.; Moreno, L. M.; Guinea, F. Plasmons and

43

40

Screening in Monolayer and Multilayer Black Phosphorus. Phys. Rev. Lett. 2014, 113,

45

106802-1-106802-5.

46

(35) Das, S.; Zhang, W.; Demarteau, M.; Hoffmann, A.; Dubey, M.; Roelofs, A. Tunable Transport Gap

47

4

in Phosphorene. Nano Lett. 2014, 14, 5733-5739.

49

(36) Popov,

50

${\Mathrm{Mos}}_{2}$ Monolayers: A Computational Study. Phys. Rev. Lett. 2012, 108,

51

48

I.;

Seifert,

G.;

Tománek,

D.

Designing

Electrical

Contacts

to

156802-1-156802-5.

53

(37) Archambault, C.; Rochefort, A. States Modulation in Graphene Nanoribbons through Metal

54

Contacts. ACS Nano 2013, 7, 5414-5420.

5

52

(38) Liu, H.; Kondo, H.; Ohno, T. Contact Effects of Nickel and Copper on Electron Transport through

57

Graphene. Phys. Rev. B 2012, 86, 155434-1-155434-7.

58

(39) Patrick Wilhite, A. A. V., Jason Tan, Jasper Tan, Toshishige Yamada, Phillip Wang Jeongwon Park

60

59

56

and Cary Y Yang Investigation on Nanotube-Metal Contacts under Different Contact Types. Mater. Lett. 2015, 145, 95-98.

17 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 3

2 (40) Maassen, J.; Ji, W.; Guo, H. First Principles Study of Electronic Transport through a Cu(111)∣

4

Graphene Junction. Appl. Phys. Lett. 2010, 97, 142105-1-142105-3.

6

(41) Gong, K.; Zhang, L.; Ji, W.; Guo, H. Electrical Contacts to Monolayer Black Phosphorus: A

7

5

First-Principles Investigation. Phys. Rev. B 2014, 90, 125441-1-125441-6.

8

(42) Brandbyge, M.; Mozos, J.-L.; Ordejón, P.; Taylor, J.; Stokbro, K. Density-Functional Method for

10

Nonequilibrium Electron Transport. Phys. Rev. B 2002, 65, 165401-1-165401-17.

1

9

(43) José, M. S.; Emilio, A.; Julian, D. G.; Alberto, G.; Javier, J.; Pablo, O.; Daniel, S.-P. The Siesta

12

Method for Ab Initio Order- N Materials Simulation. J. Phys.: Condens. Matter. 2002, 14, 2745-2779.

14

(44) Haynes, W. M. Crc Handbook of Chemistry and Physics. CRC press: Florida, U.S.A, 2013.

15

13

(45) Büttiker, M.; Imry, Y.; Landauer, R.; Pinhas, S. Generalized Many-Channel Conductance Formula

16

with Application to Small Rings. Phys. Rev. B 1985, 31, 6207-6215.

17 18 19 20 21 2 23 24 25 26 27 28 29 30 31 32 3 34 35 36 37 38 39 40 41 42 43 4 45 46 47 48 49 50 51 52 53 54 5 56 57 58 59 60

18 ACS Paragon Plus Environment

Page 18 of 19

Page 19 of 19

The Journal of Physical Chemistry Letters

1 2 3

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

19 ACS Paragon Plus Environment