Photoabsorption Tolerance of Intrinsic Point Defects and Oxidation in

Dec 14, 2016 - by employing time-dependent density functional theory, we reveal that there are two types of photoabsorption in BPQDs for both point de...
1 downloads 0 Views 3MB Size
Letter pubs.acs.org/JPCL

Photoabsorption Tolerance of Intrinsic Point Defects and Oxidation in Black Phosphorus Quantum Dots Xianghong Niu, Huabing Shu, Yunhai Li, and Jinlan Wang* Department of Physics, Southeast University, Nanjing, 211189, China S Supporting Information *

ABSTRACT: Black phosphorus quantum dots (BPQDs) exhibit excellent optical and photothermal properties and promising applications in optoelectronics and biomedicine. However, various intrinsic structural defects and oxidation are nearly unavoidable in preparation of BPQDs and how they affect the electronic and optical properties remains unclear. Here, by employing time-dependent density functional theory, we reveal that there are two types of photoabsorption in BPQDs for both point defects and oxidation. A close structure-absorption relation is unraveled: BPQDs are defect-tolerant and show excellent photoabsorption as long as the coordination number (CN) of defective P atoms is 3. By contrast, the unsaturated or oversaturated P atoms with CN ≠ 3 create in-gap-states (IGSs) and completely quench the optical absorption. An effective way to eliminate the IGSs and repair the photoabsorption of defective BPQDs via sufficient hydrogen passivation is further proposed.

S

defective BPQDs with CN = 3 are defect-tolerant and show excellent photoabsorption, while the unsaturated or oversaturated P atoms with CN ≠ 3 create IGSs and greatly damage the optical absorption. More interestingly, unlike MoS228 or BN29 edges suffering from the IGSs, hydrogen passivation can eliminate the IGSs of defective BPQDs and recover the photoabsorption effectively. All calculations were carried out by using the Becke threeparameter Lee−Yang−Parr hybrid (B3LYP)30 exchangecorrelation functional along with the Slater-type double-zata plus polarization (DZP) basis set31 as implemented in the Amsterdam Density Functional program package (ADF2013).32,33 The ground state geometries of BPQDs were first optimized at the DFT level and vibrational frequencies were calculated to verify their dynamical stability. The excitation energies were then calculated using TD-DFT. Ten lowest singlet excited states of BPQDs were taken into consideration. All optimizations were done without any symmetry constraint. Intrinsic Point Defects in BPQDs. Since the as-synthesized BPQDs are rectangular-like20 and BP monolayer is in orthorhombic lattice, we constructed square-like pristine BPQDs with zigzag or armchair edges and with edge length about 2 nm as shown in Figure 1a. The dangling bonds are passivated by hydrogen atoms. Seven different configurations were considered: divacancy, zigzag or armchair reconstructed edge, interstitial with or without H passivation, and zigzag edge or armchair edge with partial H passivation (see Figure 1a (ii−

ingle and few-layer black phosphorus (BP) has attracted worldwide research interest1−4 for its combination advantages of tunable direct bandgap,5,6 moderate on/off ratio,7,8 and sufficiently high carrier mobility,9,10 which are vital for semiconductor devices. Particularly, few-layer BP exhibits high in-plane anisotropy in optical properties,11−13 electrical conductance,14 mechanical properties,15−18 owing to the distinct puckered structure. Besides, BP quantum dots (BPQDs) have also been successfully synthesized through chemical methods19−22 and exhibited many fascinating properties such as UV/vis absorption,20 fluorescence,19,20 nearinfrared photothermal conversion performance, and biocompatibility.21 Size-dependent optical absorption and emission spectra and excellent photoabsorption in larger pristine BPQDs have also been predicted theoretically.23 These intriguing features make BPQDs promising candidates in electronics, optoelectronics, biomedicine, and so on. On the other hand, the presence of defects possibly gives rise to in-gap-states (IGSs), which normally serve as trap states and greatly degrade the performance of devices.24 Therefore, to a large extent, the success of materials in many fields is partially attributed to their tolerance to intrinsic structural defects and environment.24 It has been widely reported that intrinsic point defects and oxidation are inevitable in preparation of BP nanosheets25−27 and BPQDs19−22 using wet chemistry synthesis. However, how they influence the electronic and optical properties of BPQDs remains unclear. In this work, we provide a detailed theoretical study on this issue by using timedependent density functional theory (TD-DFT). Very interestingly, no matter for point defects (PDs) or for oxidization defects (ODs), there always exist two kinds of photoabsorption in defective BPQDs depending on the coordination number (CN) of P atoms. More specifically, the © XXXX American Chemical Society

Received: October 25, 2016 Accepted: December 14, 2016 Published: December 14, 2016 161

DOI: 10.1021/acs.jpclett.6b02486 J. Phys. Chem. Lett. 2017, 8, 161−166

Letter

The Journal of Physical Chemistry Letters

Figure 1. (a) Top and side views of the optimized ground state structures of pristine and defective BPQDs: (i) pristine, (ii) divacancy, (iii) zigzag reconstructed edge, (iv) armchair reconstructed edge, (v) interstitial, (vi) interstitial with H passivation, (vii) zigzag edge with partial H passivation, and (viii) armchair edge with partial H passivation. P and H atoms are colored in orange and white, respectively. The defective P atoms with CN = 3 (CN ≠ 3) are highlighted by magnified green (red) balls. (b) The electronic gaps of pristine and defective BPQDs. (c) Absorption spectrum of pristine and defective BPQDs. A 0.25 eV Lorentzian broadening is employed. For defective BPQDs with the CN ≠ 3, their absorption intensity is totally quenched and thus not shown.

Figure 2. (a) Density of States (DOS) of pristine and intrinsic PD BPQDs. Black lines represent total DOS, red refers to DOS of defective P atoms. The Fermi level is denoted by the purple line. All states are broadened with a Lorentzian width of 0.01 eV. (b) Illustration of the passivation process of unsaturated bond: switching CN from 2 to 3, and the dissociation process of oversaturated bond: switching CN from 4, 5, to 3 after ventilation with hydrogen. The defective P atoms with CN = 3 (CN ≠ 3) are highlighted by magnified green (red) balls. (c) Representative isosurfaces of wave functions of HOMO and LUMO in BPQDs: (i) pristine, (iii) zigzag reconstructed edge, and (iii-H) zigzag reconstructed edge with hydrogen passivation. (d) The absorption spectra for pristine BPQD (i) and zigzag reconstructed edge BPQD with hydrogen passivation (iii-H).

studies.34−37 These optimized PD BPQDs remain planar except that interstitial defect type has slight distortion. Careful

viii)). These defects have been widely studied in BP nanosheets and nanoribbons in earlier experimental or theoretical 162

DOI: 10.1021/acs.jpclett.6b02486 J. Phys. Chem. Lett. 2017, 8, 161−166

Letter

The Journal of Physical Chemistry Letters

zigzag and armchair edge of partial H passivation BPQD, in which P atoms are coordinated by 2, fully H passivation will eliminate the dangling bonds of P atoms and turn them into pristine BPQD. For the interstitial or zigzag edge reconstruction defective BPQD with CN = 4 or 5 P atoms, the excess H atoms are prone to dissociate when ventilated with hydrogen, and the saturated P bonds are thus obtained. Careful orbital analysis shows that the wave function isosurfaces of HOMO and LUMO in pristine BPQD mainly concentrates in the center region, while those of zigzag reconstructed edge BPQD completely localize in the defective region (see Figure 2c). Upon ventilation with hydrogen, the HOMO and LUMO are localized at the center region and the IGSs then disappear, and thereby the quenched absorption spectra are well recovered (see Figure 2d). Now, we have a clear picture of the influence of intrinsic PDs on the electronic structure and optical absorption of BPQDs. There exist two types of photoabsorption depending on the CN of P atoms. Specifically, the saturated P atoms with CN = 3 always show excellent photoabsorption, while the unsaturated or oversaturated P atoms with CN ≠ 3 create in-gap-states (IGSs) and fully eliminate the photoabsorption. However, these IGSs can be easily removed via sufficient hydrogen passivation, and the photoabsorption is then repaired (see Figure 3). Therefore, we can conclude BPQDs have good photoabsorption tolerance to the intrinsic PDs.

examination shows that these PD BPQDs can be categorized into two types: The P atoms at the defect region are saturated with CN = 3 (type I), and the P atoms are unsaturated or oversaturated with CN ≠ 3 (type II). More specifically, the divacancy, armchair reconstructed edge and interstitial with H passivation PDs belong to type I (highlighted by magnified green in Figure 1a), while the zigzag reconstructed edge, interstitial, zigzag or armchair edge with partial passivation PDs are type II (highlighted by magnified red in Figure 1a) with CN = 2, 4, 5, respectively. The electronic gap of BPQDs, defined as the energy difference between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), i.e., the HOMO−LUMO gap, is depicted in Figure 1b. A closer examination shows there are two kinds of responses toward the presence of defects: the electronic gap of type I BPQDs (i.e., divacancy, armchair reconstructed edge and interstitial with H passivation) is comparable to that of the pristine BPQDs, while that of type II is reduced dramatically (the zigzag reconstructed edge, interstitial, zigzag edge and armchair edge with partial H passivation). The discrepancy arises from their different bonding nature. For type I BPQDs, as the P atoms are all saturated, the bonding nature therein is similar to that of pristine BPQD. While for type II BPQDs, the defective P atoms are 2, 4, 5-coordinate, the unsaturated or oversaturated P atoms give rise to IGSs in the electronic gap as illustrated in Figure 2a. Note that this is different from the monolayer BN and metal dichalcogenides (MX2), which their IGSs are mostly originated from the chemical disorder introduced by the homoelemental bonds (M−M, X−X, B−B, or N−N).38,39 Regarding of the optical properties, as displayed in Figure 1c, the pristine BPQD exhibits striking photoabsorption in visible light region. The main absorption is dominated by transition from the ground state S0 to the first excited singlet state S1 (see Figure S1). The lowest singlet transition is mostly between the HOMO and LUMO orbitals that are mainly concentrated in the central region of BPQD (see Figure 2c). The presence of PD defects induces significantly different absorption spectra. For the type I BPQDs, the absorption spectrum is greatly broadened, and it is blue-shifted upon the presence of divacancy or interstitial defects with H passivation. The broadened and blue-shifted absorption spectra are derived from the fact that their HOMO and LUMO orbitals are partly localized in the defective region which activates high transitions such as S0 → S2 and S0 → S3 (see Figure S1,2,3). Regarding the case of armchair reconstructed edge defect, it owns similar frontier orbitals to pristine BPQD (see Figure S4), and the S0 → S1 transition remains bright excitation; Nevertheless, the reduced electronic gap leads to the red-shift of the absorption spectrum compared with the pristine BPQDs. For the type II defects, the absorption spectra are completely quenched (see Figure S5). This is because the unsaturated or oversaturated P atoms in type II BPQDs create IGSs in the bandgap and localize the frontier orbitals (see Figure S6−S8), which acts as undesirable blockage for electron excitation.40 Since the unsaturated or oversaturated P atoms in defective BPQDs cause the IGSs and quench the absorption spectra, it is natural to ask whether it is possible to remove these IGSs and recover the absorption spectra via certain treatment? The answer is positive. The IGSs in defective BPQDs can be removed easily by adding additional hydrogen atoms, i.e., further hydrogen passivation. As shown in Figure 2b, for the

Figure 3. Correlation between structure-electronic and gap-absorption in BPQDs with PDs. Type I BPQDs with P CN = 3 show good optical absorption, while the unsaturated or oversaturated P atoms with CN ≠ 3 create IGSs and quench the absorption. Interestingly, type II defects can transform into type I under the circumstance of hydrogen, and the absorption spectrum is well recovered.

Oxidation Defects in BPQDs. We constructed seven different configurations in Figure 4a, which can also be categorized into two types according to the bonding nature of P atoms. The type I BPQDs include horizontal oxygen bridge at zigzag edge, vertical oxygen bridge at armchair edge, dangling oxygen at face, interstitial oxygen bridge at surface or interlayer, with P CN = 3. The type II is formed with P CN ≠ 3 due to the unsaturated bond between phosphorus and oxygen for BPQDs with dangling oxygen at zigzag edge or at armchair edge. As shown in Figure 4b, the electronic gap of type II ODs is decreased dramatically because the pz-orbitals of edge P and O atoms form weak unsaturated P−O bonds, which induce the IGSs, as depicted in Figure S9−S11. For type I ODs, the electronic gap of oxygen bridge at zigzag or armchair edge defect is reduced, while the other ODs increase the electronic gap slightly, compared with the pristine BPQD. The discrepancy arises from their different bonding strength. Although the bridge-bonded O atom in zigzag or armchair 163

DOI: 10.1021/acs.jpclett.6b02486 J. Phys. Chem. Lett. 2017, 8, 161−166

Letter

The Journal of Physical Chemistry Letters

Figure 4. (a) Top and side views of optimized OD BPQDs. (i) dangling oxygen at zigzag edge, (ii) dangling oxygen at armchair edge, (iii) horizontal oxygen bridge at zigzag edge, (iv) vertical oxygen bridge at armchair edge, (v) dangling oxygen at face, (vi) interstitial oxygen bridge at face, and (vii) interstitial oxygen bridge at interlayer. P and H atoms are colored in orange and white, respectively. The O atoms forming saturated (unsaturated) bond with P are highlighted by magnified blue (red) ball, corresponding to type I (type II) ODs. (b) The electronic gap of OD BPQDs. (c) Absorption spectra of pristine or type I OD BPQDs. (d) Illustration of the vertical absorption excitation for the OD BPQDs. S0, S1, and Sn denote the ground state, the first singlet excited state and the nth high-lying singlet excited state, respectively. For the sake of clarity, we only marked bright excited states as red lines.

S17, the poor absorption of dangling oxygen at zigzag or armchair edge is well repaired. In summary, we systematically investigate the influence of PDs and ODs on the electronic and optical properties of BPQDs using TD-DFT. Our calculations reveal that BPQDs with PDs or ODs maintain excellent optical absorption features as long as the defective P atoms are saturated with CN = 3, and the optical absorption region is broadened. However, BPQDs containing unsaturated or oversaturated P atoms with CN ≠ 3 will completely destroy the optical absorption due to the introduction of the IGSs. Nevertheless, these IGSs can be eliminated easily by sufficient hydrogen passivation, which repairs the optical absorption of BPQDs. The defect optical tolerance enables BPQDs to serve as promising candidates for electronic and optoelectronic applications.

edge forms saturated bonds with two edge P atoms, two P−O bonds are elongated slightly, i.e., 1.76 and 1.76 Å for zigzag and 1.71 and 1.76 Å for armchair edge, respectively. The elongated P−O bonds induce edge states close to the HOMO for zigzag edge or LUMO for armchair edge defect, leading to the slightly reduced gap (see Figure S9, S12, S13), whereas, for the dangling oxygen atom at surface, the lone pair electron of P atom is attracted by oxygen atom due to large electronegative, forming the strong and polar saturated bond with a bond length of 1.51 Å. In the case of interstitial oxygen bridge at surface and interstitial oxygen bridge at interlayer, the oxygen atom penetrates into the lattice occupying the center position of the P−P bond, forming a stable bridge bond with two P atoms. The last three OD types (the dangling oxygen atom at surface, interstitial oxygen bridge at surface, and interstitial oxygen bridge at interlayer) tend to enlarge the electronic gap and frontier orbitals are delocalized (see Figure S14−S16), which was also observed in 2D phosphorene.41 Similar to the cases of PDs, there are two kinds of photoabsorption response in oxidized BPQDs. As shown in Figure 4c, the type I ODs possess excellent optical absorption in visible light region. The ODs can further tune and broaden photoabsorption region of BPQDs. The horizontal oxygen bridge at zigzag edge defect leads to the red-shift of the absorption spectrum, while the dangling oxygen at surface and interstitial oxygen bridge at surface or interlayer defects cause blue-shift of the absorption region. It is worth noting that the ODs increase the number of bright transitions compared with the only one S0 → S1 bright transition in pristine BPQD as shown in Figure 4d. However, the type II ODs are detrimental to optical absorption due to the existence of IGSs. Similar to the case of PDs, when these ODs are passivated by hydrogen, the pz-orbitals of P and O atoms form saturated P−O bond, which successfully eliminates the IGSs. As displayed in Figure



ASSOCIATED CONTENT

* Supporting Information S

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpclett.6b02486. Illustration of the vertical absorption excitation for the pristine and PD BPQDs (Figures S1,5), DOS of BPQDs and OD BPQDs (Figure S9), the frontier orbitals of PD and OD BPQDs (Figures S2−4,6−8,10−16), absorption spectra of repaired hydroxyl termination in armchair or zigzag edge of BPQDs (Figure S17) (PDF)



AUTHOR INFORMATION

ORCID

Xianghong Niu: 0000-0001-6475-8839 Huabing Shu: 0000-0001-9278-185X Yunhai Li: 0000-0003-4692-8114 Jinlan Wang: 0000-0002-4529-874X 164

DOI: 10.1021/acs.jpclett.6b02486 J. Phys. Chem. Lett. 2017, 8, 161−166

Letter

The Journal of Physical Chemistry Letters Notes

(18) Li, Y.; Hu, Z.; Lin, S.; Lai, S. K.; Ji, W.; Lau, S. P. Giant Anisotropic Raman Response of Encapsulated Ultrathin Black Phosphorus by Uniaxial Strain. Adv. Funct. Mater. 2016, DOI: 10.1002/adfm.201600986. (19) Lee, H. U.; Park, S. Y.; Lee, S. C.; Choi, S.; Seo, S.; Kim, H.; Won, J.; Choi, K.; Kang, K. S.; Park, H. G.; et al. Black Phosphorus (BP) Nanodots for Potential Biomedical Applications. Small 2016, 12, 214−219. (20) Zhang, X.; Xie, H.; Liu, Z.; Tan, C.; Luo, Z.; Li, H.; Lin, J.; Sun, L.; Chen, W.; Xu, Z.; et al. Black Phosphorus Quantum Dots. Angew. Chem., Int. Ed. 2015, 54, 3653−3657. (21) Sun, Z.; Xie, H.; Tang, S.; Yu, X. F.; Guo, Z.; Shao, J.; Zhang, H.; Huang, H.; Wang, H.; Chu, P. K. Ultrasmall Black Phosphorus Quantum Dots: Synthesis and Use as Photothermal Agents. Angew. Chem. 2015, 127, 11688−11692. (22) Gao, L. F.; Xu, J. Y.; Zhu, Z. Y.; Hu, C. X.; Zhang, L.; Wang, Q.; Zhang, H. L. Small Molecule-assisted Fabrication of Black Phosphorus Quantum Dots with A Broadband Nonlinear Optical Response. Nanoscale 2016, 8, 15132−15136. (23) Niu, X.; Li, Y.; Shu, H.; Wang, J. Anomalous Size Dependence of Optical Properties in Black Phosphorus Quantum Dots. J. Phys. Chem. Lett. 2016, 7, 370−375. (24) Zherebetskyy, D.; Zhang, Y.; Salmeron, M.; Wang, L. W. Tolerance of Intrinsic Defects in PbS Quantum Dots. J. Phys. Chem. Lett. 2015, 6, 4711−4716. (25) Castellanos-Gomez, A.; Vicarelli, L.; Prada, E.; Island, J. O.; Narasimha-Acharya, K. L.; Blanter, S. I.; Groenendijk, D. J.; Buscema, M.; Steele, G. A.; Alvarez, J. V.; et al. Isolation and Characterization of Few-layer Black Phosphorus. 2D Mater. 2014, 1, 025001. (26) Zhou, Q.; Chen, Q.; Tong, Y.; Wang, J. Light-Induced Ambient Degradation of Few-Layer Black Phosphorus: Mechanism and Protection. Angew. Chem., Int. Ed. 2016, 55, 11437−11441. (27) Yang, B.; Wan, B.; Zhou, Q.; Wang, Y.; Hu, W.; Lv, W.; Chen, Q.; Zeng, Z.; Wen, F.; Xiang, J.; et al. Te-Doped Black Phosphorus Field-Effect Transistors. Adv. Mater. 2016, 28, 9408−9415. (28) Pan, H.; Zhang, Y.-W. Edge-dependent Structural, Electronic and Magnetic Properties of MoS2 Nanoribbons. J. Mater. Chem. 2012, 22, 7280−7290. (29) Zhang, Z.; Guo, W. Energy-gap Modulation of BN Ribbons by Transverse Electric Fields: First-principles Calculations. Phys. Rev. B: Condens. Matter Mater. Phys. 2008, 77, 075403. (30) Stephens, P. J.; Devlin, F. J.; Chabalowski, C. F.; Frisch, M. J. Ab Initio Calculation of Vibrational Absorption and Circular Dichroism Spectra Using Density Functional Force Fields. J. Phys. Chem. 1994, 98, 11623−11627. (31) Van Lenthe, E.; Baerends, E. J. Optimized Slater-Type Basis Sets for the Elements 1−118. J. Comput. Chem. 2003, 24, 1142−1156. (32) te Velde, G.; Bickelhaupt, F. M.; Baerends, E. J.; Fonseca Guerra, C.; Van Gisbergen, S. J. A.; Snijders, J. G.; Ziegler, T. Chemistry with ADF. J. Comput. Chem. 2001, 22, 931−967. (33) Fonseca Guerra, C.; Snijders, J. G.; te Velde, G.; Baerends, E. J. Towards an Order-N DFT Method. Theor. Chem. Acc. 1998, 99, 391− 403. (34) Liu, Y.; Xu, F.; Zhang, Z.; Penev, E. S.; Yakobson, B. I. Twodimensional Mono-elemental Semiconductor with Electronically Inactive Defects: the Case of Phosphorus. Nano Lett. 2014, 14, 6782−6786. (35) Liang, L.; Wang, J.; Lin, W.; Sumpter, B. G.; Meunier, V.; Pan, M. Electronic Bandgap and Edge Reconstruction in Phosphorene Materials. Nano Lett. 2014, 14, 6400−6406. (36) Wang, V.; Kawazoe, Y.; Geng, W. T. Native Point Defects in Few-layer Phosphorene. Phys. Rev. B: Condens. Matter Mater. Phys. 2015, 91, 045433. (37) Zhang, R.; Zhou, X. Y.; Zhang, D.; Lou, W. K.; Zhai, F.; Chang, K. Electronic and Magneto-optical Properties of Monolayer Phosphorene Quantum Dots. 2D Mater. 2015, 2, 045012. (38) Zou, X.; Liu, Y.; Yakobson, B. I. Predicting Dislocations and Grain Boundaries in Two-dimensional Metal-disulfides From the First Principles. Nano Lett. 2013, 13, 253−258.

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work is supported by the NSFC (21525311, 21373045) and SRFDP (20130092110029) in China, the Scientific Research Foundation of Graduate School of Southeast University (YBJJ1620), and Jiangsu Innovation Projects for Graduate Students (KYZZ16_0117) in China. The authors thank the computational resources provided by Southeast University.



REFERENCES

(1) Castellanos-Gomez, A. Black Phosphorus: Narrow Gap, Wide Applications. J. Phys. Chem. Lett. 2015, 6, 4280−4291. (2) Reich, E. S. Phosphorene Excites Materials Scientists. Nature 2014, 506, 19. (3) Liu, H.; Du, Y.; Deng, Y.; Ye, P. D. Semiconducting Black Phosphorus: Synthesis, Transport Properties and Electronic Applications. Chem. Soc. Rev. 2015, 44, 2732−2743. (4) Pei, J.; Gai, X.; Yang, J.; Wang, X.; Yu, Z.; Choi, D. Y.; LutherDavies, B.; Lu, Y. Producing Air-stable Monolayers of Phosphorene and Their Defect Engineering. Nat. Commun. 2016, 7, 10450. (5) Tran, V.; Soklaski, R.; Liang, Y.; Yang, L. Layer-controlled Band Gap and Anisotropic Excitons in Few-layer Black Phosphorus. Phys. Rev. B: Condens. Matter Mater. Phys. 2014, 89, 235319. (6) 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. (7) Liu, H.; Neal, A. T.; Zhu, Z.; Luo, Z.; Xu, X.; Tománek, D.; Ye, P. D. Phosphorene: An Unexplored 2D Semiconductor with a High Hole Mobility. ACS Nano 2014, 8, 4033−4041. (8) Xia, F.; Wang, H.; Jia, Y. Rediscovering Black Phosphorus as An Anisotropic Layered Material for Optoelectronics and Electronics. Nat. Commun. 2014, 5, 4458. (9) Qiao, J.; Kong, X.; Hu, Z. X.; Yang, F.; Ji, W. High-mobility Transport Anisotropy and Linear Dichroism in Few-layer Black Phosphorus. Nat. Commun. 2014, 5, 4475. (10) Doganov, R. A.; O’Farrell, E. C.; Koenig, S. P.; Yeo, Y.; Ziletti, A.; Carvalho, A.; Campbell, D. K.; Coker, D. F.; Watanabe, K.; Taniguchi, T.; et al. Transport Properties of Pristine Few-layer Black Phosphorus by Van Der Waals Passivation in An Inert Atmosphere. Nat. Commun. 2015, 6, 6647. (11) Wang, X.; Jones, A. M.; Seyler, K. L.; Tran, V.; Jia, Y.; Zhao, H.; Wang, H.; Yang, L.; Xu, X.; Xia, F. Highly Anisotropic and Robust Excitons in Monolayer Black Phosphorus. Nat. Nanotechnol. 2015, 10, 517−521. (12) Shu, H.; Li, Y.; Niu, X.; Wang, J. The Stacking Dependent Electronic Structure and Optical Properties of Bilayer Black Phosphorus. Phys. Chem. Chem. Phys. 2016, 18, 6085−6091. (13) Guo, Z.; Zhang, H.; Lu, S.; Wang, Z.; Tang, S.; Shao, J.; Sun, Z.; Xie, H.; Wang, H.; Yu, X.-F.; et al. From Black Phosphorus to Phosphorene: Basic Solvent Exfoliation, Evolution of Raman Scattering, and Applications to Ultrafast Photonics. Adv. Funct. Mater. 2015, 25, 6996−7002. (14) Fei, R.; Yang, L. Strain-engineering the Anisotropic Electrical Conductance of Few-layer Black Phosphorus. Nano Lett. 2014, 14, 2884−2889. (15) Qin, G.; Yan, Q. B.; Qin, Z.; Yue, S. Y.; Cui, H. J.; Zheng, Q. R.; Su, G. Hinge-like Structure Induced Unusual Properties of Black Phosphorus and New Strategies to Improve the Thermoelectric Performance. Sci. Rep. 2014, 4, 6946. (16) Jiang, J. W.; Park, H. S. Negative Poisson’s Ratio in Single-layer Black Phosphorus. Nat. Commun. 2014, 5, 4727. (17) Jiang, J. W.; Rabczuk, T.; Park, H. S. A Stillinger-weber Potential for Single-layered Black Phosphorus, and the Importance of Crosspucker Interactions for a Negative Poisson’s Ratio and Edge Stressinduced Bending. Nanoscale 2015, 7, 6059−6068. 165

DOI: 10.1021/acs.jpclett.6b02486 J. Phys. Chem. Lett. 2017, 8, 161−166

Letter

The Journal of Physical Chemistry Letters (39) Liu, Y.; Zou, X.; Yakobson, B. I. Dislocations and Grain Boundaries in Two-Dimensional Boron Nitride. ACS Nano 2012, 6, 7053−7058. (40) Yu, L.; Ruzsinszky, A.; Perdew, J. P. Bending Two-Dimensional Materials To Control Charge Localization and Fermi-Level Shift. Nano Lett. 2016, 16, 2444−2449. (41) Ziletti, A.; Carvalho, A.; Trevisanutto, P. E.; Campbell, D. K.; Coker, D. F.; Castro Neto, A. H. Phosphorene Oxides: Bandgap Engineering of Phosphorene by Oxidation. Phys. Rev. B: Condens. Matter Mater. Phys. 2015, 91, 085407.

166

DOI: 10.1021/acs.jpclett.6b02486 J. Phys. Chem. Lett. 2017, 8, 161−166