Effects of Vacancy Defects on the Electronic and Optical Properties of

Jun 20, 2018 - U.S. Army Research Laboratory, Aberdeen Proving Ground, Adelphi, ... using the density functional theory plus a first-principles-based ...
0 downloads 0 Views 2MB Size
Subscriber access provided by Kaohsiung Medical University

Energy Conversion and Storage; Plasmonics and Optoelectronics

Effects of Vacancy Defects on the Electronic and Optical Properties of Monolayer PbSe Chinedu E Ekuma J. Phys. Chem. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.jpclett.8b01585 • Publication Date (Web): 20 Jun 2018 Downloaded from http://pubs.acs.org on June 21, 2018

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 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 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.

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

Effects of Vacancy Defects on the Electronic and Optical Properties of Monolayer PbSe C. E. Ekuma∗,†,‡ U.S. Army Research Laboratory, Aberdeen Proving Ground, MD 21005-5069 E-mail: [email protected]

∗ To

whom correspondence should be addressed Army Research Laboratory, Aberdeen Proving Ground, MD 21005-5069 ‡ George F. Adams Research Fellow † U.S.

1

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 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 Defect-engineering is promising for tailoring the properties of atomically thin materials. By creating defects via Se vacancies, we study the optoelectronic properties of monolayer PbSe. We obtain the single-particle properties using the density functional theory plus a firstprinciples-based typical medium approximation. The absorption spectra is explored by solving the Bethe-Salpeter equation. Our results reveal that monolayer PbSe is defect sensitive but defect-tolerant. The latter fingerprint is due to the absence of defect-induced in-gap, localized states. Our results predict that Se vacancies are the dominant defect-type in disordered PbSe monolayer. We observe that increasing Se vacancy concentrations δ renormalize the energy bandgap Eg increases from 0.21 eV for the pristine to as high as 0.45 eV at high δ . The high tunability of the optoelectronic properties of monolayer PbSe using defect-engineering makes this material a candidate for exploring flexible electronic with potential technological applications in nanoelectronics.

TOC Graphic

Among the rapidly evolving class of two-dimensional (2D) materials, monolayer lead chalcogenides and their alloys, such as PbSe, are intensively being studied for next-generation nanoelectronics applications. Their unusual and unique properties relative to other semiconductors make them attractive candidate materials for both potential device applications and for exploring fundamental physical phenomena. Some of their unique properties are high mobility, high dielectric 2

ACS Paragon Plus Environment

Page 2 of 20

Page 3 of 20 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 Letters

constant, low-cost, and a narrow fundamental bandgap. Because of these attractive properties, they have been extensively studied for diverse applications such as thermoelectric energy converters and photovoltaics, 1–3 spintronics and optoelectronic devices, 4–6 infrared lasers and diodes. 7,8 Recently, the nanostructure of PbSe and related materials have been explored for creating new materials with novel and tunable properties. 9–12 There is by now a general consensus on the origin of free carriers, e.g., n- (p-) type semiconductor originates from Se (Pb) vacancies in bulk PbSe. 13,14 However, in the case of PbSe monolayer, the nature and stability of the defect states have yet to be systematically studied. Defects inevitably exist in crystals, which could be unintentionally introduced during experimental sample preparation or material characterization. In the case of atomically thin 2D materials, due to their high susceptibility to disorder, the impurity states could dramatically affect their properties for device applications. Defects in materials could be studied with the density functional theory (DFT) 16,17 using the supercell. However, the DFT supercell approach can only describe ordered defect structures. Even more importantly, the supercell approach is prohibitively expensive for exploring low defect concentrations. Also, supercell model is not adequate for calculating the energy of isolated defect states especially “shallow” impurity states in semiconductors 18 and the finite-size effects of the supercell have been reported to adversely diminish the dispersion of the deep-level in-gap states. 19 Depending on the location of the defect-induced states with respect to the bandgap, exciton recombination processes could be radiative or nonradiative. 20–22 For instance, defect-induced localized, deep-level, in-gap states could be detrimental to device performance as they act as local scattering sites, which facilitates nonradiative recombination of photoexcited excitons vis-à-vis, the reduction of the mobility of charge carriers and quantum efficiency. This is one of the major limiting factors of materials’ applications in optoelectronics devices such as field-effect transistors, solar cells, photodetectors, and light emitters. Hence, a proper understanding of the nature of the defect-induced states: deep-level localized states or shallow impurity, resonant levels within the bandgap is of great importance. The former occurs if the defect-induced states lie inside the

3

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 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: (A) The top-view of the atomic structure and (B) the first Brillouin zone of monolayer PbSe along the [001]-direction. (C). The density functional theory band structure (black dashed bands) reproduced by the downfolded Hamiltonian H0 (solid red bands) for the pristine monolayer of PbSe. A direct single-particle bandgap 0.21 eV is predicted at the given X-point of the high symmetry zones of the Brillouin zone. The side plot is the associated density of states of the Pb-p, Se-p states, and the quasiparticle total density of states obtained from the self-consistent quasiparticle Green function and screened Coulomb interaction calculations using VASP. 15 The horizontal black dash line is the Fermi level.

4

ACS Paragon Plus Environment

Page 4 of 20

Page 5 of 20 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 Letters

bandgap. They behave as dominant electron or hole traps as well as recombination centers while the latter is when the defect-induced states are located outside the bandgap. They generally led to a small perturbation to the lattice potential reminiscent of substitutional dopant. Understanding the defect morphology could greatly aid in exploiting the full technological potentials, e.g., topological states, optoelectronic, and other diverse promising applications of monolayer PbSe. Despite the need, no computational data exploring the nature of the defect states in monolayer PbSe have been published so far to our knowledge. In this letter, we use the recently developed first-principles based typical medium dynamical cluster approximation (TMDCA) 23,24 that takes into account the intrinsic, random nature of impurity states in crystals to systemically study the optoelectronic properties of disordered monolayer PbSe. We determine the stability of the defect states in the disordered crystal by calculating the defect formation energy ∆E f , which is the bond energy cost due to the creation of a defect. The calculated ∆E f profiles show that Se vacancies are energetically more favorable to form in monolayer PbSe. Our results reveal that randomly distributed Se vacancies could impact the optoelectronic properties of monolayer PbSe in diverse ways. We show that defect-engineering could be used to tune the electronic properties of monolayer PbSe. Our data show that monolayer PbSe is sensitive to defects, but it is, however, defect-tolerant 25–27 due to the absence of localized, deep-trapped levels within the bandgap. The optoelectronic fingerprint of monolayer PbSe makes it a potential candidate material to explore flexible nanoelectronics and for solar cells applications. To study disorder induced by vacancy in monolayer PbSe, we use Anderson Hamiltonian of noninteracting electrons subjected to random vacancy configurations

H = H0 + ∑ Viσα nαiσ ,

(1)

iασ

where H0 is the pristine, single-particle Hamiltonian, Viσα is a disorder potential, and nαiσ is the number operator; i, α, and σ are site, orbital, and spin indices, respectively. The terms in Eq. 1 denote the single-particle and disorder in that order.

5

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 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 single-particle Hamiltonian is obtained from DFT 16,17 calculations using the linearized augmented planewave (LAPW) method with a modified Becke-Johnson (mBJ) potential, 28 as implemented in WIEN2K. 29 The mBJ potential yields energy bandgap of semiconductors in reasonable agreement with experiment in contrast to standard density functionals and at a fraction of the computational cost of the many-body perturbation theory, e.g., the Green function plus screened Coulomb approximations. 6,28,30,31 We first carried out structural optimization using the PerdewBurke-Ernzerhof 32 exchange-correlation functional and then self-consistent electronic structure calculations of the periodically repeated monolayers of PbSe separated ≈ ∼20 Å vacuum using the optimized lattice constant to obtain the Kohn-Sham eigenvalues and eigenfunctions. The DFT electronic structure calculation used a well-converged basis set with a dense Brillouin zone sampling grid of 12 ×12×1 ~k-points and an Rmin Kmax = 9, where Rmin is the smallest sphere radius in the system and Kmax is the plane wave momentum cutoff. A LAPW sphere radius of 2.50 Bohr was used for Pb and Se, respectively. All calculations were performed relativistically including spinorbit coupling, which could be important for the band structures. 33 Using downfolding technique as implemented in the WANNIER90, 34 we obtain H0 from the converged Kohn-Sham eigenstates by generating a set of Wannier functions. The H0 included the Pb-p and Se-p states, which are the dominant states around the Fermi level. This observation is in agreement with our recent bulk structure calculations. 6 We show in Fig. 1(B) the first Brillouin zone of PbSe [001] monolayer. The 12-bands that make up H0 reproduce the DFT band structure, with a direct single-particle energy bandgap Eg ≈ 0.21 eV at the X-point of the k-space [Fig. 1(C)]. We model the disorder using a bimodal site potential Viσα = Vi ∈ {0,Wb } randomly generated using the probability mass functions P(Vi = W ) = δ and P(Vi = W ) = 0, where 0 (W ) is for the ordered (disordered) site. The vacant site potential energy is set to be far greater than the pristine bandwidth. Taking the Se vacancies as an example, the average Se vacancy concentration δ satisfies the stoichiometry PbSe1−δ of the disordered material. The TMDCA, through a set of self-consistent field (SCF) equations, maps the disordered lattice onto a periodically finite cluster of Nc primitive cells embedded in a typical medium characterized

6

ACS Paragon Plus Environment

Page 6 of 20

Page 7 of 20 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 Letters

by a hybridization function. This effective hybridization function is calculated self-consistently, employing the TMDCA condition. 23,35–42 The TMDCA SCF equations were iteratively solved to obtain the cluster density of states ρ c = − π1 ℑGc , which is averaged over a large number of configurations to calculate the wave-vector-resolved, non-self-averaged cluster typical density of states (for details, see Ref. 23 ). To calculate the two-particle Green functions, we used the converged single-particle Green functions obtained from the TMDCA self-consistency calculations to solve the Bethe-Salpeter equations as described in Ref. 24 The TMDCA allows for the study of very small impurity concentrations as low as 0.10% since it uses the same primitive cell of the host material without the need for supercell. While this is advantageous since the computation cost does not scale with the system size as in supercell approaches, it, however, neglects possible structural relaxation processes around the vacancies. To verify the effects of this neglect, we performed relaxed supercell calculations of different sizes: 2×2×1 [Fig. 2(A)] and 4×4×1 (not shown) corresponding to 25% and 6.25% Se vacancy concentrations. The DFT supercell calculations used the mBJ functional as implemented in the Vienna Ab Initio Simulation Package (VASP). 15 We show results obtained with the TMDCA using typical (indigo dashed line) and algebraic (green solid line) averaging methods [Fig. 2(A)]. Both the DFT and TMDCA calculations show single-particle bandgap enhancement. Most importantly, the DFT result shows no defect-induced deep-level states. As the answer was negative since the defect-induced states were all shallow impurity levels, one may expect that the main conclusions of this work, related to PbSe monolayer being defect-tolerant material, do not depend on crystal structure relaxation or spin-orbit interaction. Before proceeding to the discussion of our main results, we will first consider the properties of the pristine PbSe monolayer. The band structure and the associated partial density of states are shown in Fig. 1(C). Our results show that monolayer PbSe is a direct and narrow bandgap semiconductor with a single-particle energy bandgap Eg ≈ 0.21 eV at the X-point of the k-space. The calculated Eg of monolayer PbSe is about 5.0% smaller than the bulk Eg ≈ 0.28 eV. 6 Based on the crystal symmetry, we believe that the smaller Eg of the PbSe monolayer is justified. Aside from the difference in the crystal symmetry of the monolayer (tetragonal symmetry) and the bulk

7

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

(cubic symmetry), the monolayer structure could have a markedly stronger crystal field effects. 43 Also, due to quantum confinement, smaller Eg of the monolayer PbSe is consistent with increased pressure that could decrease the bandgap since the lattice constant of the monolayer PbSe a0 = 4.25 Å [Fig. 1(A)] is smaller than the bulk lattice constant a0 = 6.12 Å. 44

(PbSe)1-δ PbSe1-δ Pb1-δSe

DFT TMDCA

0.4

-1

∆fE (eV)

4

(B)

(A)

D(E) (eV )

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

Page 8 of 20

2

0

0.2

0 E (eV)

0

5

10 15 δ (%)

20

0

Figure 2: (A) A comparison of the density of states obtained using the supercell approach with the TMDCA@DFT for Se vacancy concentration δ = 25%. The DFT calculations used the mBJ potential. The indigo dashed line correspond to the typical density of states. Observe that in all cases, the defect-induced states are all resonant levels with no localized deep-level states within the single-particle bandgap. (B) The vacancy formation energy as a function of Pb, Se, and Pb-Se vacancy concentrations. To gain some insight on the energetics of the vacancy formation processes and the stability of the disordered monolayer PbSe crystal, we calculated the defect formation energy ∆E f = Ev − E0 , where Ev and E0 are the total energies of the vacancy and pristine structure, respectively. We checked the stability of the disordered crystal for the Se-rich, Pb-rich, and when we have vacancies on both the Pb- and Se-sublattices. As shown in Fig. 2(B), ∆E f is positive and increases with 8

ACS Paragon Plus Environment

Page 9 of 20

δ = 2%

δ = 1%

δ = 5%

δ = 10%

3 -1

D(E) (eV )

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 Letters

2 1 0

-2

0

2

-2

0

2 -2 E (eV)

0

2

-2

0

2

Figure 3: TMDCA@DFT results of monolayer PbSe showing the typical density of states at four different Se vacancy concentrations δ . For reference, we also includes the density of states of pristine monolayer PbSe obtained from the mBJ potential. The Fermi level has been shifted to zero at the top of the valence bands. The energy bandgap systematically increased from 0.21 eV for the pristine crystal to 0.45 eV for the disordered crystal at δ = 10%. increasing vacancy concentration. Our data reveal that disorder due to Se vacancies is the most probable to form in monolayer PbSe crystal since the defect formation energy for vacancies on Se-sites are greater than that of Pb-sites and for the case where the Pb- and Se-sublattices are both disordered. Since our calculations show that Se vacancy is the most stable defect structure, unless otherwise stated, the reported results below are for the effects of Se vacancies on the optoelectronic properties of monolayer PbSe. Next, we investigate the effects of Se vacancies on the electronic properties of monolayer PbSe. Figure 3 shows the typical density of states at [low (1 and 2%) and high (5 and 10%)] Se vacancy concentrations δ . As δ increases (Fig. 3), the single-particle bandgap systematically increases as 0.21, 0.24, 0.27, 0.31, and 0.45 for the pristine, 1, 2, 5, and 10% Se vacancy concentrations, respectively. The Fermi level also shifts to higher energy defined with respect to δ = 0, implying n-type doping. The data for Pb vacancies (not shown) show p-doping as expected. In the case where the vacancies are introduced in both Pb and Se sublattices, it was dominant n-doping at low concentrations but evolves to a behavior reminiscent of p − n doping character at high concentra9

ACS Paragon Plus Environment

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

tions. This could be significant as one could engineer a p − n junction by creating vacancies in both sublattices. Experimentally, this could be achieved by modulating the threshold energy for knock-on damage, which may be markedly different for Pb and Se. Our data show significant redistribution and the enhancement of states around the Fermi level. The defect-induced redistribution and enhancement of the states around the Fermi level without the formation of isolated impurity states 45 have been shown to lead to increased thermoelectric figure of merit, e.g., in PbTe, SnTe, and Bi2 Te3 alloys. 46–51 Observe that for all the Se vacancy concentrations considered herein, aside from the formation of some resonant (virtual bound) states, there are no localized, deep-level defect states within the single-particle bandgap. To properly discern the exact band characters and their contribution to the states in the proximity of the Fermi level, we count the number of states in the band extrema. The character of the states in the proximity of the band extrema [valence band maximum (VBM) and conduction band minimum (CBM)] have been shown to be crucial to whether the defect induced impurity states will be shallow or localized, deep-trapped states. 25–27 As shown in Fig. 4, the single-particle bandgap of monolayer PbSe is predominantly derived from Pb-p and Se-p states. The conduction band minimum is formed by a strong hybridization between Pb-pz and Se-p, which contributed ∼ 64% and 36%, respectively of the total CBM states. The valence band maximum is composed mainly of Se-p, which contributed ∼ 92% of the total VBM states. When disorder due to Se vacancies is introduced in the crystal, while the above picture is generally maintained, the states were renormalization and the lattice symmetry-induced degenerate px+y orbital characters are lifted to form distinct px and py states. The electronic signature of the states that form the single-particle bandgap, namely, the antibonding character of the conduction and bonding character of the valence bands extrema will favor the bonds broken due to Se vacancies to mostly reside inside the bands as resonant states leaving the single-particle bandgap free of defect-induced deep-trapped localized states as schematically shown in Fig. 4 (A). This could be attributed to a symmetry-protected property already observed in some semiconductors and such materials are said to be defect-tolerant. 25–27 This also suggests that

10

ACS Paragon Plus Environment

Page 10 of 20

Page 11 of 20 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 Letters

the defect-induced states are reminiscent of shallow, intrinsic defects with minimal or absence of surface states. 14,26 This picture is different from what is observed in some semiconductors where the defect-induced impurity states appear as deep-trapped levels inside the bandgap [Fig. 4(B)], which could adversely diminish the quantum efficiency of the material for device applications.

Figure 4: Schematics of the orbital character of the valence band maximum and conduction band minimum in a pristine and disordered crystal for (A) monolayer PbSe and (B) a conventional semiconductor of orbital characters α and β , respectively, where α (β ) could be a combination or any of the orbitals: s, p, and d. The absorption spectroscopy provides a detailed and direct information about the electronic signature of materials and has the advantage of being a true “bulk” probe. 6 While there have been a number of optical (experimental and theoretical) studies of both pristine and disordered bulk PbSe, 52–61 there are yet no comprehensive computational studies of the optical properties of disordered monolayer PbSe. Herein, we provide a detailed first-principles based calculations of the absorption spectra of disordered monolayer PbSe. We calculate the absorption spectra, which includes the effects of electron-hole interactions by solving the Bethe-Salpeter equation in the particle-hole channel to obtain the complex dielectric function ε(ω). The two particle calculations used the single-particle Green functions obtained from the converged, self-consistent TMDCA calculations.

11

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

30

20

ω4

Pristine 1% 2% 5% 10%

40

ε2(ω)

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

ω3

Eb ω1

Page 12 of 20

ω2

10

0 0

Egq

2

4 h ω (eV) _

6

8

Figure 5: The imaginary part of the dielectric function for monolayer PbSe as a function of the excitation energy h¯ ω at various Se vacancy concentrations δ . The ω1−4 denote the main exciton peaks. We obtain exciton binding energy Eb = 0.20, 0.24, 0.27, 0.38, and 0.42 eV for the pristine, 1, 2, 5, and 10% vacancy concentrations in that order. The vertical solid black line is the quasiparticle q energy bandgap Eg ≈ 1.51 eV. Increasing the concentration initially leads to the increase and then suppression of the generation of free electron-hole pairs at higher excitation energies. We show in Fig. 5 the imaginary part of the dielectric function ε2 (ω) for various Se vacancy concentrations. The impact of the Se vacancies on the magnitude of ε2 (ω) is an initial increase due to electron-hole delocalization and then a decrease as the vacancy concentration increases. The latter effects could be attributed to spectral weight redistribution. The first peak in the absorption spectra (ω1 ∼ 1.30 eV) is related to the band edge. As Se vacancy concentration increase, ω1 is redshifted, broadened, and decreased systematically to ∼1.0 at δ = 20%. Another prominent structure in the absorption spectra occurs at ω2 with a broad shoulder observed at ω3 . The maximum of the spectra is depicted with ω4 , which occurs at ∼5.4 eV. We note that the absorption spectra of the pristine monolayer and the bulk crystal are qualitatively similar (see, e.g., Refs. 6,55 ). However, key differences are obvious, especially at lower energies where the two resonant features below 2.50 eV (ω1 and ω2 ) is just a shoulder ∼1.90 eV in the bulk structure. 6,55 The differences could be attributed to, among other things, the crystal symmetry (tetragonal versus rock-salt), dimensionality (twodimension versus three-dimension), and the inclusion of the effects of electron-hole interactions in our current data.

12

ACS Paragon Plus Environment

Page 13 of 20 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 Letters

q

To estimate the exciton binding energy, we calculate the quasiparticle energy bandgap Eg ≈ 1.51 eV from the self-consistent quasiparticle Green function and screened Coulomb interaction q

calculations as implemented in VASP. 15 The difference between Eg and the first excitation energy (onset of a prominent feature) in the absorption spectra is the measure of the exciton binding energy Eb . As a function of increasing Se vacancy concentrations, we obtain Eb ≈ 0.20, 0.24, 0.27, 0.38, and 0.42 eV for δ = 0, 1, 2, 5, and 10%, in that order. As expected, our calculated fundamental bandgap is lower than the quasiparticle bandgap due to the exciton binding energy. Interestingly, the calculated exciton binding energies seem high for such a narrow bandgap material. For example, the exciton binding energy in monolayer MoS2 range from 0.2–0.7 eV. 62–67 This is significant since a bound exciton is a fingerprint of a stable material properties. With such Eb , the excitons could, e.g., withstand strong electric fields and high carrier injections. Also, the calculated absorption spectra show relatively large absorption around 1.0 – 2.5 eV, which could be explored for potential solar cell applications. In summary, we have studied the effects of randomly distributed Se vacancies on the electronic and optical properties of monolayer PbSe using a first-principles-based typical medium approach. Our results show that the defect-induced states in disordered monolayer PbSe are predominantly shallow impurity resonant levels, which mostly resides inside the bandgap, a signature of a defecttolerant material. Our calculations further show that the single-particle bandgap of monolayer PbSe could be modulated by Se vacancies since increased δ led to bandgap renormalization. The calculated exciton binding energies further show that monolayer PbSe is stable even at high vacancy concentration. The research was sponsored by the Army Research Laboratory (ARL) and was accomplished under the Cooperative Agreement Number W911NF-11-2-0030 as an ARL Research [George F. Adams] Fellow. The views and conclusions contained in this document are those of the authors and should not be interpreted as representing the official policies, either expressed or implied, of ARL or the U.S. Government. The U.S. Government is authorized to reproduce and distribute reprints for Government purposes notwithstanding any copyright notation herein. This work was supported

13

ACS Paragon Plus Environment

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

in part by a grant of computer time from the DOD High-Performance Computing Modernization Program at the Army Engineer Research and Development Center, Vicksburg, MS.

References (1) Ioffe, A. F. Semiconductor Thermoelements and Thermoelectric and Thermoelectric Cooling; Inforsearch, London, 1957. (2) Wood, C. Materials for Thermoelectric Energy Conversion. Reports on Progress in Physics 1988, 51, 459–539. (3) Rowe, M. CRC Handbook of Thermoelectrics; CRC Press, Boca Raton, 1995. (4) Jin, S.; Wu, H.; Xu, T. Large Rashba splitting in Highly Asymmetric CdTe/PbTe/PbSrTe Quantum Well Structures. Appl. Phys. Lett. 2009, 95, 132105-1–132105-3 . (5) Akimov, B.; Dmitriev, A.; Khohlov, D.; Ryabova, L. Carrier Transport and Non-equilibrium Phenomena in Doped PbTe and Related Materials. Phys. Status Solidi A 1993, 137, 9–55. (6) Ekuma, C. E.; Singh, D. J.; Moreno, J.; Jarrell, M. Optical Properties of PbTe and PbSe. Phys. Rev. B 2012, 85, 085205-1–085205-7 . (7) Preier, H. Recent Advances in Lead-chalcogenide Diode Lasers. Applied physics 1979, 20, 189–206. (8) John, J.; Zogg, H. Infrared P-N-junction Diodes in Epitaxial Narrow Gap PbTe layers on Si Substrates. J. Appl. Phys. 1989, 85, 3364–3367 . (9) Alimoradi, J. M.; Janssen, V. A. E. C.; Evers, W. H.; Tadjine, A.; Delerue, C.; Siebbeles, L. D. A.; van der Zant, H. S. J.; Houtepen, A. J.; Vanmaekelbergh, D. Transport Properties of a Two-Dimensional PbSe Square Superstructure in an Electrolyte-Gated Transistor. Nano Letters 2017, 17, 5238–5243. (10) Zaiats, G.; Shapiro, A.; Yanover, D.; Kauffmann, Y.; Sashchiuk, A.; Lifshitz, E. Optical and Electronic Properties of Nonconcentric PbSe/CdSe Colloidal Quantum Dots. J. Phys. Chem. Lett. 2015, 6, 2444– 2448.

14

ACS Paragon Plus Environment

Page 14 of 20

Page 15 of 20 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 Letters

(11) Groiss, H.; Kaufmann, E.; Springholz, G.; Schwarzl, T.; Hesser, G.; Schäffler, F.; Heiss, W.; Koike, K.; Itakura, T.; Hotei, T.; Yano, M.; Wojtowicz, T. Size Control and Midinfrared Emission of Epitaxial PbTe/CdTe Quantum Dot Precipitates Grown by Molecular Beam Epitaxy. Appl. Phys. Lett. 2007, 91, 222106-1–222106-3. (12) Shen, J.; Jung, Y.; Disa, A. S.; Walker, F. J.; Ahn, C. H.; Cha, J. J. Synthesis of SnTe Nanoplates with 100 and 111 Surfaces. Nano Letters 2014, 14, 4183–4188. (13) Breschi, R.; Camanzi, A.; Fano, V. Defects in PbTe Single Crystals. J. Cryst. Growth 1982, 58, 399– 408 . (14) Bauer, G.; Burkhard, H.; Heinrich, H.; Lopez-Otero, A. Impurity and vacancy states in PbTe. J. Appl. Phys. 1976, 47, 1721–1723. (15) Kresse, G.; Furthmüller, J. Efficiency of Ab-initio Total Energy Calculations for Metals and Semiconductors Using a Plane-wave Basis Set. Comput. Mater. Sci. 1996, 6, 15–50. (16) Hohenberg, P.; Kohn, W. Inhomogeneous Electron Gas. Phys. Rev. 1964, 136, B864–B871. (17) Kohn, W.; Sham, L. J. Self-consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133–A1138 . (18) Pandey, M.; Rasmussen, F. A.; Kuhar, K.; Olsen, T.; Jacobsen, K. W.; Thygesen, K. S. Defect-Tolerant Monolayer Transition Metal Dichalcogenides. Nano Letters 2016, 16, 2234–2239. (19) Schultz, P. A. Theory of Defect Levels and the “Band Gap Problem” in Silicon. Phys. Rev. Lett. 2006, 96, 246401-1–246401-4 . (20) Luque, A.; Martí, A.; Antolín, E.; Tablero, C. Intermediate Bands Versus Levels in Non-radiative Recombination. Physica B 2006, 382, 320–327. (21) Luque, A.; Marti, A.; Stanley, C. Understanding Intermediate-band Solar Cells. Nat Photon 2012, 6, 146–152.

15

ACS Paragon Plus Environment

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

(22) Casavola, M.; van Huis, M. A.; Bals, S.; Lambert, K.; Hens, Z.; Vanmaekelbergh, D. Anisotropic Cation Exchange in PbSe/CdSe Core/Shell Nanocrystals of Different Geometry. Chemistry of Materials 2012, 24, 294–302. (23) Ekuma, C. E.; Dobrosavljevi´c, V.; Gunlycke, D. First-Principles-Based Method for Electron Localization: Application to Monolayer Hexagonal Boron Nitride. Phys. Rev. Lett. 2017, 118, 106404-1– 106404-5 . (24) Ekuma, C.E.; Gunlycke, D. Optical Absorption in Disordered Monolayer Molybdenum Disulfide. Phys. Rev. B 2018, 97, 201414-1–201414-6 (25) Ahmad, S.; Mahanti, S. D.; Hoang, K.; Kanatzidis, M. G. Ab Initio Studies of the Electronic Structure of Defects in PbTe. Phys. Rev. B 2006, 74, 155205-1–155205-13 . (26) Zakutayev, A.; Caskey, C. M.; Fioretti, A. N.; Ginley, D. S.; Vidal, J.; Stevanovic, V.; Tea, E.; Lany, S. Defect Tolerant Semiconductors for Solar Energy Conversion. J. Phys. Chem. Lett. 2014, 5, 1117– 1125. (27) Brandt, R. E.; Stevanovi´c, V.; Ginley, D. S.; Buonassisi, T. Identifying Defect-tolerant Semiconductors with High Minority-carrier Lifetimes: Beyond Hybrid Lead Halide Perovskites. MRS Communications 2015, 5, 265–275. (28) Tran, F.; Blaha, P. Accurate Band Gaps of Semiconductors and Insulators with a Semilocal Exchangecorrelation Potential. Phys. Rev. Lett. 2009, 102, 226401-1–226401-4. (29) Blaha, P.; Schwarz, K.; Madsen, G.; Kvasnicka, D.; Luitz, J. WIEN2K, An Augmented Plane Wave+Local Orbitals Program for Calculating Crystal Structure; K. Schwarz Technical University, Wien, Australia, 2001. (30) Singh, D. J. Electronic Structure Calculations with the Tran-Blaha Modified Becke-Johnson Density Functional. Phys. Rev. B 2010, 82, 205102-1–205102-10. (31) Singh, D. J. Structure and Optical Properties of High Light Output Halide Scintillators. Phys. Rev. B 2010, 82, 155145-1–155145-10.

16

ACS Paragon Plus Environment

Page 16 of 20

Page 17 of 20 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 Letters

(32) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. (33) Wei, S. H.; Zunger, A. Electronic and Structural Anomalies in Lead Chalcogenides. Phys. Rev. B 1997, 55, 13605–13610 . (34) Mostofi, A. A.; Yates, J. R.; Lee, Y.-S.; Souza, I.; Vanderbilt, D.; Marzari, N. wannier90: A Tool for Obtaining Maximally-localised Wannier Functions. Comp. Phys. Commun. 2008, 178, 685–699. (35) Ekuma, C. E.; Terletska, H.; Tam, K.-M.; Meng, Z.-Y.; Moreno, J.; Jarrell, M. Typical Medium Dynamical Cluster Approximation for the Study of Anderson Localization in Three Dimensions. Phys. Rev. B 2014, 89, 081107-1–081107-5. (36) Ekuma, C. E.; Moore, C.; Terletska, H.; Tam, K.-M.; Moreno, J.; Jarrell, M.; Vidhyadhiraja, N. S. Finite-cluster Typical Medium Theory for Disordered Electronic Systems. Phys. Rev. B 2015, 92, 014209-1–014209-18. (37) Ekuma, C. E.; Terletska, H.; Meng, Z. Y.; Moreno, J.; Jarrell, M.; Mahmoudian, S.; Dobrosavljevi´c, V. Effective Cluster Typical Medium Theory for the Diagonal Anderson Disorder Model in One- and Two-dimensions. J. Phys. Conden. Matt. 2014, 26, 274209. (38) Zhang, Y.; Terletska, H.; Moore, C.; Ekuma, C.; Tam, K.-M.; Berlijn, T.; Ku, W.; Moreno, J.; Jarrell, M. Study of Multiband Disordered Systems using the Typical Medium Dynamical Cluster Approximation. Phys. Rev. B 2015, 92, 205111-1–205111-9. (39) Jarrell, M.; Maier, T.; Huscroft, C.; Moukouri, S. Quantum Monte Carlo Algorithm for Nonlocal Corrections to the Dynamical Mean-field Approximation. Phys. Rev. B 2001, 64, 195130-1–19513023. (40) Hettler, M. H.; Mukherjee, M.; Jarrell, M.; Krishnamurthy, H. R. Dynamical Cluster Approximation: Nonlocal Dynamics of Correlated Electron Systems. Phys. Rev. B 2000, 61, 12739–12756. (41) Maier, T.; Jarrell, M.; Pruschke, T.; Hettler, M. H. Quantum Cluster Theories. Rev. Mod. Phys. 2005, 77, 1027–1080.

17

ACS Paragon Plus Environment

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

(42) Ekuma, C. E. Towards the Realization of Systematic, Self-Consistent Typical Medium Theory for Interacting Disordered Systems; LSU Doctoral Dissertations. 2391, 2015. (43) Liu, J.; Qian, X.; Fu, L. Crystal Field Effect Induced Topological Crystalline Insulators in Monolayer IV-VI Semiconductors. Nano Letters 2015, 15, 2657–2661 . (44) Madelung, O.; Rössler, U.; Schulz, M. E. Numerical Data and Functional Relationships in Science and Technology, Landolt-Börnstein, New Series, Group III; Springer, 2006; Vol. 17a, and 22a. (45) Wiendlocha, B. Fermi Surface and Electron Dispersion of PbTe Doped with Resonant Tl Impurity from KKR-CPA Calculations. Phys. Rev. B 2013, 88, 205205-1–205205-7 . (46) Heremans, J. P.; Jovovic, V.; Toberer, E. S.; Saramat, A.; Kurosaki, K.; Charoenphakdee, A.; Yamanaka, S.; Snyder, G. J. Enhancement of Thermoelectric Efficiency in PbTe by Distortion of the Electronic Density of States. Science 2008, 321, 554–557. (47) Jaworski, C. M.; Wiendlocha, B.; Jovovic, V.; Heremans, J. P. Combining Alloy Scattering of Phonons and Resonant Electronic Levels to Reach a High Thermoelectric Figure of Merit in PbTeSe and PbTeS Alloys. Energy Environ. Sci. 2011, 4, 4155–4162. (48) Jaworski, C. M.; Kulbachinskii, V.; Heremans, J. P. Resonant Level Formed by Tin in Bi2 Te3 and the Enhancement of Room-temperature Thermoelectric Power. Phys. Rev. B 2009, 80, 233201-1–233201. (49) Zhang, Q.; Wang, H.; Zhang, Q.; Liu, W.; Yu, B.; Wang, H.; Wang, D.; Ni, G.; Chen, G.; Ren, Z. Effect of Silicon and Sodium on Thermoelectric Properties of Thallium-Doped Lead Telluride-Based Materials. Nano Letters 2012, 12, 2324–2330 . (50) Zhang, Q.; Cao, F.; Liu, W.; Lukas, K.; Yu, B.; Chen, S.; Opeil, C.; Broido, D.; Chen, G.; Ren, Z. Heavy Doping and Band Engineering by Potassium to Improve the Thermoelectric Figure of Merit in p-Type PbTe, PbSe, and PbTe1−y Sey . J. Am. Chem. Soc. 2012, 134, 10031–10038 . (51) Zhang, Q.; Liao, B.; Lan, Y.; Lukas, K.; Liu, W.; Esfarjani, K.; Opeil, C.; Broido, D.; Chen, G.; Ren, Z. High Thermoelectric Performance by Resonant Dopant Indium in Nanostructured SnTe. Proc. Natl. Acad. Sci. U. S. A. 2013, 110, 13261–13266.

18

ACS Paragon Plus Environment

Page 18 of 20

Page 19 of 20 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 Letters

(52) Lawson, W. D. A Method of Growing Single Crystals of Lead Telluride and Lead Selenide. J. Appl. Phys. 1951, 22, 1444–1447 . (53) Dalven, R. in Solid State Physics, H. Ehrenreich, F. Seitz, and D. Turnbell, (Eds.); Academic Press, New York, 1973; Vol. 28. (54) Kohn, S. E.; Yu, P. Y.; Petroff, Y.; Shen, Y. R.; Tsang, Y.; Cohen, M. L. Electronic Band Structure and Optical Properties of PbTe, PbSe, and PbS. Phys. Rev. B 1973, 8, 1477–1488. (55) Suzuki, N.; Sawai, K.; Adachi, S. Optical Properties of PbSe. J. Appl. Phys. 1995, 77, 1249–1255 . (56) Cardona, M.; Greenaway, D. L. Optical Properties and Band Structure of Group IV-VI and Group V Materials. Phys. Rev. 1964, 133, A1685–A1697 . (57) Moss, T. Optical Properties of Semiconductors; Butterworth, London, 1959. (58) Bogacki, F. J.; Sood, A. K.; Yang, C. Y.; Rabii, S.; Fischer, J. E. Thermoreflectance of IV-VI Compounds. Surf. Sci. 1973, 37, 494–507. (59) Bartnik, A. C.; Efros, A. L.; Koh, W.-K.; Murray, C. B.; Wise, F. W. Electronic States and Optical Properties of PbSe Nanorods and Nanowires. Phys. Rev. B 2010, 82, 195313-1–195313-16 . (60) Li, X.-B.; Guo, P.; Zhang, Y.-N.; Peng, R.-F.; Zhang, H.; Liu, L.-M. High Carrier Mobility of Fewlayer PbX (X = S, Se, Te). J. Mater. Chem. C 2015, 3, 6284–6290 . (61) Gao, Y.; Sandeep, C. S. S.; Schins, J. M.; Houtepen, A. J.; Siebbeles, L. D. A. Disorder Strongly Enhances Auger Recombination in Conductive Quantum-dot Solids. Nat. Commun. 2013, 4, 2329 . (62) Radisavljevic, B.; Kis, A. Mobility Engineering and a Metal-Insulator Transition in Monolayer MoS2 . Nat. Mater. 2013, 12, 815–820. (63) Rigosi, A. F.; Hill, H. M.; Rim, K. T.; Flynn, G. W.; Heinz, T. F. Electronic Band Gaps and Exciton Binding Energies in Monolayer MoxW1−x S2 Transition Metal Dichalcogenide Alloys Probed by Scanning Tunneling and Optical Spectroscopy. Phys. Rev. B 2016, 94, 075440-1–075440-6 .

19

ACS Paragon Plus Environment

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

(64) Hill, H. M.; Rigosi, A. F.; Roquelet, C.; Chernikov, A.; Berkelbach, T. C.; Reichman, D. R.; Hybertsen, M. S.; Brus, L. E.; Heinz, T. F. Observation of Excitonic Rydberg States in Monolayer MoS2 and WS2 by Photoluminescence Excitation Spectroscopy. Nano Letters 2015, 15, 2992–2997 . (65) Zhang, C.; Johnson, A.; Hsu, C.-L.; Li, L.-J.; Shih, C.-K. Direct Imaging of Band Profile in Single Layer MoS2 on Graphite: Quasiparticle Energy Gap, Metallic Edge States, and Edge Band Bending. Nano Letters 2014, 14, 2443–2447 . (66) Klots, A. R. et al. Probing Excitonic States in Suspended Two-Dimensional Semiconductors by Photocurrent Spectroscopy. Scientific Reports 2014, 4, 6608 . (67) Yu, Y.; Yu, Y.; Cai, Y.; Li, W.; Gurarslan, A.; Peelaers, H.; Aspnes, D. E.; Van de Walle, C. G.; Nguyen, N. V.; Zhang, Y.-W.; Cao, L. Exciton-dominated Dielectric Function of Atomically Thin MoS2 Films. Scientific Reports 2015, 5, 16996.

20

ACS Paragon Plus Environment

Page 20 of 20