Defects in Phosphorene - The Journal of Physical ... - ACS Publications

Aug 18, 2015 - The SV-(5|9) and DV-(5|8|5)-2 defects can introduce unoccupied localized states into phosphorene's fundamental band gap. Specifically, ...
0 downloads 8 Views 13MB Size
Subscriber access provided by UNIV OF CALIFORNIA SAN DIEGO LIBRARIES

Article

Defects in Phosphorene Wei Hu, and Jinlong Yang J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.5b06077 • Publication Date (Web): 18 Aug 2015 Downloaded from http://pubs.acs.org on August 19, 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 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 26

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

Defects in Phosphorene Wei Hu∗,† and Jinlong Yang∗,‡,¶ †Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA ‡Hefei National Laboratory for Physical Sciences at Microscale, University of Science and Technology of China, Hefei, Anhui 230026, China ¶Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China E-mail: [email protected]; [email protected] Phone: +1-5105416034; +86-55163606408

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: Defects are inevitably present in materials and always can affect their properties. Here, We perform first-principles calculations to systematically investigate the stability and electronic structures of ten kinds of point defects in two-dimensional semiconducting phosphorene, including the Stone-Wales (SW-1 and SW-2) defect, single (SV-(5|9) and SV-(55|66)) and double (DV-(5|8|5)-1, DV-(5|8|5)-2, DV-(555|777)-1, DV-(555|777)-2, DV-(555|777)-3 and DV-(4|10|4)) vacancy defects. We find that these defects are all created quite easily in phosphorene with higher areal density compared with graphene and silicene. Most of them are easy to distinguish each other and correlate with their defective atomic structures with simulated scanning tunneling microscopy images at positive bias. The SW, DV-(5|8|5)-1, DV-(555|777) and DV-(4|10|4) defects have little effect on phosphorene’s electronic properties and defective phosphorene monolayers still show semiconducting with similar band gap values to perfect phosphorene. The SV-(5|9) and DV-(5|8|5)-2 defects can introduce unoccupied localized states into phosphorene’s fundamental band gap. Specifically, the SV-(5|9) and SV-(55|66) defects can induce hole doping in phosphorene, and the SV-(5|9) defect can result in local magnetic moments in phosphorene different from all other defects. Keywords: Point Defects, Phosphorene, Electronic Structures, and Density Functional Theory.

Introduction Two-dimensional (2D) atom-thick materials, 1–4 for example, graphene, 5 silicene, 6 hexagonal boron nitride (h-BN) 7 and molybdenum disulphide (MoS2 ), 8 show remarkable properties and potential applications on semiconductors, electronics and composites, and have recently attracted significant interest. In particular, graphene, 5 a 2D sp2 -hybridized carbon elemental monolayer, shows high performance optoelectronic properties. For example, graphene has a high carrier mobility, however, its zero-gap feature limits its wide applications on graphene2 ACS Paragon Plus Environment

Page 2 of 26

Page 3 of 26

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

based nanoelectronics with large-off current and high on-off ratio. The same limitation also exists in silicene, 6 another well-known single silicon monolayer, which has most similar remarkable electronic properties to graphene but with buckled honeycomb structures. 9 Recently, a new 2D phosphorus elemental monolayer, namely phosphorene, 10–20 has been isolated experimentally from bulk black phosphorus through mechanical exfoliation. Compared with graphene and silicene, phosphorene shows more excellent optoelectronic properties and has immediately received outstanding attention. In particular, phosphorene is direct semiconducting, and it has a high hole mobility 10 up to 105 cm2 /V/s and a large drain current modulation 11 up to 105 in electronic devices. Furthermore, phosphorene possesses a thickness-dependent direct bandgap values from 0.3 eV (Bulk) to 1.5 eV (Monolayer). 10 Like graphene, 21 phosphorene also has high performance optical properties. 12 These advantages make phosphorene an ideal candidate for field effect transistors, 11 photovoltaic PN junctions 13 and thin-film solar cells. 15 Therefore, phosphorene can be used as an alternative to graphene, and may lead to new faster semiconductor optoelectronic devices in the future experiments. On the other hand, most properties and applications of materials are always affected by the presence of defects, especially, point defects, 22 which are generated usually by ion or electron irradiations. 23 For 2D materials, defects also can be extremely used to generate innovative devices. 24 Typical point defects in graphene and silicene include Stone-Wales (SW) defect, single and double vacancy (SV and DV) defects. 25–30 Generally, both graphene and silicene have two kinds of SVs (SV-(5|9) and SV-(55|66)) and DVs (DV-(5|8|5) and DV(555|777)), respectively. These defects are inevitably present in graphene and silicene and severely affect their structural and electronic properties, 28,29 thus alter their applications. 31–33 However, little attention has been focused on the defects in phosphorene. 34 In this work, we systematically investigate the stability and electronic structures of typical point defects in semiconducting phosphorene using the density functional theory and ab-initio molecular dynamics calculations. We find that phosphorene has a wide variety of defects

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

due to its low symmetry structure. Furthermore, these defects are all quite easily created in phosphorene with regard to graphene and silicene and most of them are easy distinguish each other and correlate with their defective atomic structures with simulated scanning tunneling microscopy images at positive bias. These defects have different geometric structures and show different stability and electronic structures in phosphorene.

Theoretical Methods and Models We use the density functional theory (DFT) calculations implemented in the VASP package. 35 We choose the generalized gradient approximation of Perdew, Burke, and Ernzerhof (GGA-PBE) 36 owing to its good description of electronic structures of phosphorene, 14 graphene 28 and silicene. 29 All the atomic coordinates are fully relaxed by using the conjugate gradient (CG) algorithm until total energy and atomic forces are converged to 10−5 eV and 0.01 eV/˚ A, respectively. The unit cell lattice parameters of phosphorene is calculated to be a(P) = 4.62 ˚ A and b(P) = 3.30 ˚ A. 14 A large 5 × 7 supercell of phosphorene (140 phosphorus atoms) with a vacuum space about 15 ˚ A in the Z direction is adopted to study the effect of various local defects in phosphorene. All the phosphorus atoms in phosphorene are labeled by numbers as shown in Figure 1. We set the energy cutoff to be 500 eV and adopt a 3 × 3 regular mesh to sample the surface Brillouin zone. Ab-initio molecular dynamics (AIMD) simulations are employed to study the stability of defects in phosphorene. AIMD simulations are performed in a canonical ensemble. The simulations is performed during 2.0 ps with a time step of 2.0 fs at the temperature of 400 K controlled by using the Nose-Hoover method. 37

Results and Discussion We first investigate the geometric properties of defects in phosphorene as shown in Figure 1. Typical point defects, such as SW, SV-(5|9), SV-(55|66), DV-(5|8|5) and DV-(555|777) usu4 ACS Paragon Plus Environment

Page 4 of 26

Page 5 of 26

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

ally generated in graphene 28 and silicene, 29 also can be formed in phosphorene, but which has a wide variety of defects because of lower symmetry structure of phosphorene compared with graphene and silicene. We find that there are two kinds of SW (denoted as SW-1 and SW-2) and DV-(5|8|5) (denoted as DV-(5|8|5)-1 and DV-(5|8|5)-2), three kinds of DV-(555|777) (denoted as DV-(555|777)-1, DV-(555|777)-2 and DV-(555|777)-3) defects in phosphorene as shown in Figure 1. Furthermore, we find a new stable structure of DV defect in phosphorene, namely DV-(4|10|4), as shown in Figure 1(h), which has not been observed in both graphene and silicene. For two kinds of SW defects as shown in Figure 1(a) and (b), different P-P bonds in phosphorene are rotated by about 90 degrees, P67 -P74 and P59 -P67 bonds for SW-1 and SW-2 defects, respectively. Two kinds of SV defects are both produced by removing the P67 atom from phosphorene as shown in Figure 1(d) and (e), but the distances between the P60 and P74 atoms are 3.60 and 2.53 ˚ A in the SV-(5|9) and SV-(55|66) defects, respectively. Thus, P60 and P67 atoms form a P-P bond in SV-(55|66) defect. There are six kinds of DV defects in phosphorene. For two kinds of DV-(5|8|5) defects as shown in Figure 1(f) and (g), different P-P bonds containing two phosphorus atoms are removed from phosphorene, P60 -P67 and P67 -P74 bonds for DV-(5|8|5)-1 and DV-(5|8|5)-2 defects, respectively. The DV(4|10|4) defect is also produced by removing the P67 -P74 bond from phosphorene, but they show different bonding configurations. In the DV-(5|8|5)-2 defect, P59 -P60 and P80 -P81 form two chemical bonds. But in the DV-(4|10|4) defect, P59 -P80 and P60 -P81 form two chemical bonds. Three configurations of DV-(555|777) defects are complex as shown in Figure 1(i), (j) and (k). Different from the DV-(5|8|5) and DV-(4|10|4) defects, they are all produced by removing two non-bonding-together phosphorus atoms, P53 and P59 , P67 and P68 , P60 and P75 for DV-(555|777)-1, DV-(555|777)-2 and DV-(555|777)-3 defects, respectively. To check the stability of defects in phosphorene, the formation energy is defined as

Ef = EPhosphorene − NP ∗ EP

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

Figure 1: (Color online) Geometric structures of (a) perfect and defective phosphorene in the 5 × 7 supercell, including the (b) SW-1, (c) SW-2, (d) SV-(5|9), (e) SV-(55|66), (f) DV(5|8|5)-1, (g) DV-(5|8|5)-2, (h) DV-(4|10|4), (i) DV-(555|777)-1, (j) DV-(555|777)-2 and (k) DV-(555|777)-3 defects. In figure (a), all the phosphorus atoms are labeled by numbers and the red and blue balls denote the phosphorus atoms in two different upper and lower layers of phosphorene, respectively. In other figures, the violet and yellow balls denote unaffected and affected phosphorus atoms, respectively. 6 ACS Paragon Plus Environment

Page 6 of 26

Page 7 of 26

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

where EPhosphorene represents the total energy of defective phosphorene, EP is the energy per phosphorus atom in a perfect phosphorene sheet and NP corresponds to the number of phosphorus atoms in phosphorene. Notice that for perfect phosphorene, Ef (Phosphorene) = 0 eV. Perfect phosphorene shows less stable with a smaller cohesive energy of 3.48 eV/atom compared with graphene (7.90 eV/atom) and silicene (3.96 eV/atom). 29 Our calculated stability of perfect and defective phosphorene, graphene and silicene are summarized in Table 1 and Table 2. Table 1: Calculated stability and electronic properties of perfect and defective phosphorene, including the formation energy Ef (eV ), total magnetic moment µ (µB ) and band gap Eg (eV ). Ef Perfect 0.000 SW-1 1.012 SW-2 1.322 SV-(5|9) 1.626 SV-(55|66) 2.025 DV-(5|8|5)-1 1.906 DV-(5|8|5)-2 3.041 DV-(4|10|4) 2.137 DV-(555|777)-1 2.081 DV-(555|777)-2 2.350 DV-(555|777)-3 2.613

µ 0.000 0.003 0.003 0.980 0.000 0.002 0.002 0.000 0.002 0.002 0.003

Eg 0.905 0.928 0.883 0.190/0.941 Hole doping 0.962 0.559 0.929 0.966 0.962 0.973

Table 2: Comparison results of cohesive energy Ec (eV/atom) of perfect phosphorene, graphene and silicene with corresponding various defects’ formation energy Ef (eV). Phosphorene Reference This work Ec (Perfect) 3.48 Ef (SW) 1.01−1.32 Ef (SV-(5|9)) 1.63 Ef (SV-(55|66)) 2.03 Ef (DV-(5|8|5)) 1.91−3.04 Ef (DV-(555|777)) 2.08−2.61 Ef (DV-(4|10|4)) 2.13

Graphene Ref. 28 7.90 4.50 7.80 − 7.52 6.40 −

7 ACS Paragon Plus Environment

Silicene Ref. 29 3.96 2.09 3.77 3.01 3.70 2.84 −

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

We find that the SW-1 defect is most easily formed in phosphorene with the smallest formation energy of 1.01 eV among various defects similar to graphene 28 and silicene. 29 For DVs, the DV-(5|8|5) and DV-(555|777) defects are also stable in phosphorene similar to graphene and silicene. Interestingly, we find the most stable DV in phosphorene is DV(4|10|4) with a small formation energy of 2.13 eV, but which can not formed in graphene and silicene. To further study the thermal stability of perfect and defective phosphorene, AIMD simulations in the VASP are performed at 400 K. In the initial state (t = 0.0 ps), perfect and defective phosphorene monolayers are set to optimized geometric structures. We find that the perfect phosphorene and most of defective phosphorene monolayers are stable at 400 K during t = 2.0 ps. The only unstable one is the SV-(55|66) defect, which trends to change into more stable SV-(5|9) defect in phosphorene at 400 K, different from the situation of silicene. 29 At finite temperature T , the areal density NDefect (m−2 ) of defects in 2D materials follows the Arrhenius equation NDefect = NPerfect exp(−Ef /kB T ) where NPerfect is the areal density of atoms in perfect 2D materials, Ef is the formation energy of a defect formed in materials and kB is the Boltzmann constant. For perfect phosphorene, graphene and silicene, their areal densities are NPerfect (Phosphorene) = 2.62 × 1019 m−2 , NPerfect (Graphene) = 3.79 × 1019 m−2 and NPerfect (Silicene) = 1.55 × 1019 m−2 , respectively. The temperature-dependence areal density of the most stable defects in phosphorene, graphene and silicene is calculated and plotted in Figure 2. The results show that these defects have much higher areal density and they are created quite easily in phosphorene compared with graphene and silicene. In order to help recognize these point defects in the future experiments, the scanning tunneling microscopy (STM) images of perfect and defective phosphorene are also simulated at +1.0 and -1.0 V bias as shown in Figure 3. At positive bias (+1.0 V), the STM images 8 ACS Paragon Plus Environment

Page 8 of 26

Page 9 of 26

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 2: (Color online) Areal density of various stable defects (SW, SV-(5|9), SV-(55|66), DV-(5|8|5), DV-(555|777) and DV-(4|10|4)) in (a) phosphorene, (b) graphene and (c) silicene as a function of temperature. of most defects are easy to understand and correlate with their defective atomic structures. But at negative bias (-1.0 V), the STM images of perfect and defective phosphorene are difficult to distinguish each other due to similar height variation for the buckled structures of phosphorene. The different phenomena of STM images of perfect and defective phosphorene can be understood by their projected density of states (PDOS) of a selected phosphorus atom in perfect and defective phosphorene as shown in Figure 4. We find that the valence bands of perfect phosphorene are mainly contributed by P3pz orbitals associated with lone pairs of electrons, and its conduction bands are composed of the sp3 hybridized states of P3s, P3px , P3py and P3pz orbitals (Figure 4(a)). Thus, for perfect phosphorene, there are symmetrical white bright spots with a distance of 3.54 ˚ A in the STM image at -1.0 V bias (valence bands), but symmetrical gray rings in the STM image at +1.0 V bias (conduction bands). For defective phosphorene, the defect-affected phosphorus atoms near the defects all have some occupied states contributing into the valence bands of phosphorene, but the peak strength and positions of these occupied states are mostly overlapped with P3pz states of other unaffected phosphorus atoms. Therefore, perfect and defective phosphorene show similar STM images at negative bias (-1.0 V) and thus they are difficult to distinguish each 9 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

Figure 3: (Color online) Simulated STM images (+1.0 and -1.0 V) of (a) perfect and defective phosphorene with the (b) SW-1, (c) SW-2, (d) SV-(5|9), (e) SV-(55|66), (f) DV-(5|8|5)-1, (g) DV-(5|8|5)-2, (h) DV-(4|10|4), (i) DV-(555|777)-1, (j) DV-(555|777)-2 and (k) DV-(555|777)3 defects.

10 ACS Paragon Plus Environment

Page 10 of 26

Page 11 of 26

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

other at this situation.

Figure 4: (Color online) Projected density of states (PDOS) of a selected phosphorus atom in (a) perfect and defective phosphorene, including (b) P67 in SW-1, (c) P59 in SW-2, (d) P60 in SV-(5|9) (spin-up plus spin-down states), (e) P74 in SV-(55|66), (f) P59 in DV-(5|8|5)-1, (g) P59 in DV-(5|8|5)-2, (h) P59 in DV-(4|10|4), (i) P60 in DV-(555|777)-1, (j) P60 in DV(555|777)-2 and (k) P68 in DV-(555|777)-3. The Fermi level is marked by green dotted lines and set to zero.

However, defective phosphorene systems show stronger unoccupied states than that of perfect phosphorene, and that is why most defects in phosphorene can be recognized by using the STM images at positive bias. The SV-(55|66) defect in phosphorene is most easily recognized in the STM image, because there is only a white bright spot at +1.0 V bias (Figure 4(e)), which is contributed by the stronger unoccupied states from the defectaffected phosphorus atom P74 (Figure 1(e)). For the SV-(5|9) defect in phosphorene, there are four grey spots at +1.0 V bias (Figure 4(d)), which are contributed by four defect-affected 11 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

phosphorus atoms P59 , P60 , P68 and P74 (Figure 1(d)). The DV-(5|8|5)-2 and DV-(4|10|4) defects in phosphorene are also easily recognized in the STM image at +1.0 V bias due to their high symmetry (Figure 4(g) and (h)). For the DV-(5|8|5)-1 defect in phosphorene, there are a white bright spot and several gray spots at +1.0 V bias (Figure 4(f)). For two kinds of SW defects in phosphorene, they can be distinguished by comparing the atomic structures and STM images at positive bias (Figure 4(b) and (c)). But, three DV-(555|777) defects in phosphorene show similar and disordered STM images at positive bias (Figure 4(i), (j) and (k)) due to their lower symmetry and similar configurations. Further advanced methods need to be introduced to observe them in the experiments. Finally, we check the electronic band structures of perfect and defective phosphorene as shown in Figure 5. Monolayer phosphorene is a direct semiconductor with a band gap of 0.91 eV (Figure 5(a)), agreeing well with previous theoretical studies. 14 We find that the SW, DV-(555|777) and DV-(4|10|4) defects have little effect on phosphorene’s electronic properties, still showing semiconducting with similar band gap values (about 0.9 eV) to perfect phosphorene, different from graphene 28 and silicene. 29 That is because the occupied and unoccupied states contributed by the defect-effaced phosphorus atoms for these defects are respectively close that of other unaffected phosphorus atoms as shown in Figure 4. Defective phosphorene with the DV-(5|8|5) defects are also semiconductors but with different band gaps (0.96 and 0.56 eV respectively for DV-(5|8|5)-1 and DV-(5|8|5)-2). The SV(5|9) and DV-(5|8|5)-2 defects can introduce unoccupied localized states into phosphorene’s fundamental band gap as shown in Figure 5(d) and (g). The introduced states of SV-(5|9) and DV-(5|8|5)-2 defects are contributed by the P60 and P59 atom as shown in Figure 4(d) and (g), respectively. The P60 atom in SV-(5|9) defect shows a dangling bond different from other phosphorus atoms in phosphorene which can form three covalent P-P bonds and a lone pair of electrons via sp3 orbital hybridization of P3s, P3px , P3py and P3pz orbitals. The P59 -P60 atoms in the DV-(5|8|5)-2 defect are bonding very weekly with a longer bond length of 2.54 ˚ A compared with other P-P covalent bonds (2.22 ˚ A).

12 ACS Paragon Plus Environment

Page 12 of 26

Page 13 of 26

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 5: (Color online) Electronic band structures of (a) perfect and defective phosphorene with the (b) SW-1, (c) SW-2, (d) SV-(5|9), (e) SV-(55|66), (f) DV-(5|8|5)-1, (g) DV-(5|8|5)2, (h) DV-(4|10|4), (i) DV-(555|777)-1, (j) DV-(555|777)-2 and (k) DV-(555|777)-3 defects. The Fermi level is marked by green dotted lines and set to zero.

13 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

More interestingly, we find that two SV defects can induce hole doping in phosphorene as shown in Figure 5(d) and (e), which can increase the hole carrier concentration of semiconducting phosphorene. Notice that phosphorene has a high areal density for the SV defects, agreeing well with a high hole mobility observed in phosphorene in recent experiments. 10 Furthermore, among all the point defects in phosphorene, only the SV-(5|9) defect can induce local magnetic moments in phosphorene with a magnetic moment of 0.98 µB due to the presence of dangling bond in P60 atom, and it is ferromagnetic as shown in Figure 6. But, other defective phosphorene monolayers are not magnetic similar to the situation of silicene. 29 The magnetic property of the SV-(5|9) defect makes phosphorene an extremely promising material for spin qubits 38 and quantum spintronics. 39

Figure 6: (Color online) Spin-density isosurfaces for the SV-(5|9) defect in phosphorene. The red and blue color isosurfaces represent spin-up and spin-down states, respectively.

Conclusions We systematically study the stability and electronic structures of defects in semiconducting phosphorene using the density functional theory and ab-initio molecular dynamics calculations. We find that phosphorene has a wide variety of point defects (SW-1, SW-2, SV-(5|9), SV-(55|66), DV-(5|8|5)-1, DV-(5|8|5)-2, DV-(555|777)-1, DV-(555|777)-2, DV-(555|777)-3 14 ACS Paragon Plus Environment

Page 14 of 26

Page 15 of 26

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

and DV-(4|10|4)) due to its low symmetry structure. Furthermore, these defects are all quite easily created in phosphorene compared with graphene and silicene and most of them are easy to distinguish each other and correlate with their defective atomic structures with simulated scanning tunneling microscopy images at positive bias. Defects of different structures shows different stability and electronic structures in phosphorene. The SW, DV-(5|8|5)-1, DV-(555|777) and DV-(4|10|4) defects have little effect on phosphorene’s electronic properties and defective phosphorene monolayers still show semiconducting with similar band gap values to perfect phosphorene. The SV-(5|9) and DV-(5|8|5)-2 defects can introduce unoccupied localized states into phosphorene’s fundamental band gap. Specifically, the SV-(5|9) and SV-(55|66) defects can induce hole doping in phosphorene, and the SV-(5|9) defect can result in local magnetic moments in phosphorene different from all other defects. Our theoretical calculations provide significant insights into the identification of point defects in the further experiments and the understanding of their effects on the electronic properties and potential applications of phosphorene.

Acknowledgement This work is partially supported by the Scientific Discovery through Advanced Computing (SciDAC) Program funded by U.S. Department of Energy, Office of Science, Advanced Scientific Computing Research and Basic Energy Sciences (W. H.). This work is also partially supported by the National Key Basic Research Program (2011CB921404), by NSFC (21421063, 91021004, 21233007), by Chinese Academy of Sciences (CAS) (XDB01020300). We thank the National Energy Research Scientific Computing (NERSC) center, USTCSCC, SCCAS, Tianjin, and Shanghai Supercomputer Centers for the computational resources.

15 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

References (1) Novoselov, K. S.; Jiang, D.; Schedin, F.; Booth, T. J.; Khotkevich, V. V.; Morozov, S. V.; Geim, A. K. Two-Dimensional Atomic Crystals. Proc. Natl. Acad. Sci. USA 2005, 102, 10451–10453. (2) Mas-Ballest´e, R.; G´ omez-Navarro, C.; G´omez-Herrero, J.; Zamora, F. 2D Materials: to Graphene and Beyond. Nanoscale 2011, 3, 20–30. (3) Osada, M.; Sasaki, T. Two-Dimensional Dielectric Nanosheets: Novel Nanoelectronics from Nanocrystal Building Blocks. Adv. Mater. 2012, 24, 210–228. (4) Xu, M.; Liang, T.; Shi, M.; Chen, H. Graphene-Like Two-Dimensional Materials. Chem. Rev. 2013, 113, 3766–3798. (5) Geim, A. K.; Novoselov, K. S. The Rise of Graphene. Nature Mater. 2007, 6, 183–191. (6) Cahangirov, S.; Topsakal, M.; Akt¨ urk, E.; S¸ahin, H.; Ciraci, S. Two- and One-Dimensional Honeycomb Structures of Silicon and Germanium. Phys. Rev. Lett. 2009, 102, 236804. (7) Watanabe, K.; Taniguchi, T.; Kanda, H. Direct-Bandgap Properties and Evidence for Ultraviolet Lasing of Hexagonal Boron Nitride Single Crystal. Nature Mater. 2004, 3, 404–409. (8) Radisavljevic, B.; Radenovic, A.; Brivio, J.; Giacometti, V.; Kis, A. Single-Layer MoS2 Transistors. Nature Nanotech. 2011, 6, 147–150. (9) Hu, W.; Li, Z.; Yang, J. Structural, Electronic, and Optical Properties of Hybrid Silicene and Graphene Nanocomposite. J. Chem. Phys. 2013, 139, 154704. (10) Liu, H.; Neal, A. T.; Zhu, Z.; Luo, Z.; Xu, X.; Tom´anek, D.; Ye, P. D. Phosphorene: An Unexplored 2D Semiconductor with a High Hole Mobility. ACS Nano 2014, 8, 4033–4041. (11) Li, L.; Yu, Y.; Ye, G.; Ge, Q.; Ou, X.; Wu, H.; Feng, D.; Chen, X.; Zhang, Y. Black Phosphorus Field-Effect Transistors. Nature Nanotech. 2014, 9, 372–377. (12) Qiao, J.; Kong, X.; Hu, Z.-X.; Yang, F.; Ji, W. High-Mobility Transport Anisotropy and Linear Dichroism in Few-Layer Black Phosphorus. Nature Commun. 2014, 5, 4475. (13) Buscema, M.; Groenendijk, D. J.; Steele, G. A.; van der Zant, H. S.; Castellanos-Gomez, A. Photovoltaic Effect in Few-Layer Black Phosphorus PN Junctions Defined by Local Electrostatic Gating. Nature Commun. 2014, 5, 4651.

16 ACS Paragon Plus Environment

Page 16 of 26

Page 17 of 26

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

(14) Guo, H.; Lu, N.; Dai, J.; Wu, X.; Zeng, X. C. Phosphorene Nanoribbons, Phosphorus Nanotubes, and van der Waals Multilayers. J. Phys. Chem. C 2014, 118, 14051–14059. (15) Dai, J.; Zeng, X. C. Bilayer Phosphorene: Effect of Stacking Order on Bandgap and Its Potential Applications in Thin-Film Solar Cells. J. Phys. Chem. Lett. 2014, 5, 1289–1293. (16) Kou, L.; Frauenheim, T.; Chen, C. Phosphorene as a Superior Gas Sensor: Selective Adsorption and Distinct I-V Response. J. Phys. Chem. Lett. 2014, 5, 2675–2681. (17) Kou, L.; Ma, Y.; Smith, S. C.; Chen, C. Anisotropic Ripple Deformation in Phosphorene. J. Phys. Chem. Lett. 2015, 6, 1509–1513. (18) Kou, L.; Chen, C.; Smith, S. C. Phosphorene: Fabrication, Properties, and Applications. J. Phys. Chem. Lett. 2015, 6, 2794–2805. (19) Hu, W.; Wang, T.; Yang, J. Tunable Schottky Contacts in Hybrid Graphene-Phosphorene Nanocomposites. J. Mater. Chem. C 2015, 3, 4756–4761. (20) Hu, W.; Lin, L.; Yang, C. Edge Reconstruction in Armchair Phosphorene Nanoribbons Revealed by Discontinuous Galerkin Density Functional Theory. Phys. Chem. Chem. Phys. 2015, Accepted, DOI: 10.1039/C5CP00333D. (21) Bonaccorso, F.; Sun, Z.; Hasan, T.; Ferrari, A. C. Graphene Photonics and Optoelectronics. Nature Photon. 2010, 4, 611–622. (22) Freysoldt, C.; Grabowski, B.; Hickel, T.; Neugebauer, J.; Kresse, G.; Janotti, A.; de Walle, C. G. V. First-Principles Calculations for Point Defects in Solids. Rev. Mod. Phys. 2014, 86, 253. (23) Komsa, H.-P.; Kotakoski, J.; Kurasch, S.; Lehtinen, O.; Kaiser, U.; Krasheninnikov, A. V. TwoDimensional Transition Metal Dichalcogenides under Electron Irradiation: Defect Production and Doping. Phys. Rev. Lett. 2012, 109, 035503. (24) Zou, X.; Yakobson, B. I. An Open Canvas - 2D Materials with Defects, Disorder, and Functionality. Acc. Chem. Res. 2015, 48, 73–80. (25) Hashimoto, A.; Suenaga, K.; Gloter, A.; Urita, K.; Iijima, S. Direct Evidence for Atomic Defects in Graphene Layers. Nature 2004, 430, 870–873.

17 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

(26) Amorim, R. G.; Fazzio, A.; Antonelli, A.; Novaes, F. D.; da Silva, A. J. R. Divacancies in Graphene and Carbon Nanotubes. Nano Lett. 2007, 7, 2459–2462. (27) Kou, L.; Tang, C.; Guo, W.; Chen, C. Tunable Magnetism in Strained Graphene with Topological Line Defect. ACS Nano 2011, 5, 1012–1017. (28) Banhart, F.; Kotakoski, J.; Krasheninnikov, A. V. Structural Defects in Graphene. ACS Nano 2011, 5, 26–41. (29) Gaoab, J.; Zhang, J.; Liu, H.; Zhang, Q.; Zhao, J. Structures, Mobilities, Electronic and Magnetic Properties of Point Defects in Silicene. Nanoscale 2013, 5, 9785–9792. (30) Sahin, H.; Sivek, J.; Li, S.; Partoens, B.; Peeters, F. M. Stone-Wales Defects in Silicene: Formation, Stability, and Reactivity of Defect Sites. Phys. Rev. B 2013, 88, 045434. (31) Leenaerts, O.; Partoens, B.; Peeters, F. M. Graphene: A Perfect Nanoballoon. Appl. Phys. Lett. 2008, 93, 193107. (32) Hu, W.; Wu, X.; Li, Z.; Yang, J. Helium Separation via Porous Silicene based Ultimate Membrane. Nanoscale 2013, 5, 9062–9066. (33) Hu, W.; Wu, X.; Li, Z.; Yang, J. Porous Silicene as a Hydrogen Purification Membrane. Phys. Chem. Chem. Phys. 2013, 15, 5753–5757. (34) Liu, Y.; Xu, F.; Zhang, Z.; Penev, E. S.; Yakobson, B. I. Two-Dimensional Mono-Elemental Semiconductor with Electronically Inactive Defects: The Case of Phosphorus. Nano Lett. 2014, 14, 6782–6786. (35) Kresse, G.; Hafner, J. Ab Initio Molecular Dynamics for Liquid Metals. Phys. Rev. B 1993, 47, 558–561. (36) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. (37) Khanna, S. N.; Jena, P. Atomic Clusters: Building Blocks for a Class of Solids. Phys. Rev. B 1995, 51, 13705–13716. (38) Pedersen, T. G.; Flindt, C.; Pedersen, J.; Mortensen, N. A.; Jauho, A.-P.; Pedersen, K. Graphene Antidot Lattices: Designed Defects and Spin Qubits. Phys. Rev. Lett. 2008, 100, 136804. (39) Lundeberg, M. B.; Yang, R.; Renard, J.; Folk, J. A. Defect-Mediated Spin Relaxation and Dephasing in Graphene. Phys. Rev. Lett. 2013, 110, 156601.

18 ACS Paragon Plus Environment

Page 18 of 26

Page 19 of 26

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

19 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

ACS Paragon Plus Environment

Page 20 of 26

Page 21 of 26

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

The Journal of Physical Chemistry

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

ACS Paragon Plus Environment

Page 22 of 26

Page 23 of 26

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

The Journal of Physical Chemistry

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

ACS Paragon Plus Environment

Page 24 of 26

Page 25 of 26

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

The Journal of Physical Chemistry

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

ACS Paragon Plus Environment

Page 26 of 26