Electronic and optical properties of two-dimensional GaN from first

Nocona Sanders, Dylan Bayerl, Guangsha Shi, Kelsey A. Mengle, and Emmanouil Kioupakis*. Department of Materials ...... corrections to DFT/LDA band str...
0 downloads 6 Views 455KB Size
Subscriber access provided by READING UNIV

Communication

Electronic and optical properties of two-dimensional GaN from first principles Nocona Sanders, Dylan Bayerl, Guangsha Shi, Kelsey A. Mengle, and Emmanouil Kioupakis Nano Lett., Just Accepted Manuscript • DOI: 10.1021/acs.nanolett.7b03003 • Publication Date (Web): 25 Oct 2017 Downloaded from http://pubs.acs.org on October 26, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Nano Letters is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 19

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

Nano Letters

Electronic and optical properties of two-dimensional GaN from first principles Nocona Sanders, Dylan Bayerl, Guangsha Shi, Kelsey A. Mengle, and Emmanouil Kioupakis* Department of Materials Science and Engineering, University of Michigan, Ann Arbor 48109, United States Corresponding Author *E-mail: [email protected] *Phone: 734-764-3321

ACS Paragon Plus Environment

1

Nano 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

Page 2 of 19

ABSTRACT: GaN is an important commercial semiconductor for solid-state lighting applications. Atomically thin GaN, a recently synthesized two-dimensional material, is of particular interest because the extreme quantum confinement enables additional control of its light-emitting properties. We performed first-principles calculations based on density functional and many-body perturbation theory to investigate the electronic, optical, and excitonic properties of monolayer and bilayer 2D GaN as a function of strain. Our results demonstrate that light emission from monolayer 2D GaN is blueshifted into the deep ultraviolet range, which is promising for sterilization and water-purification applications. Light emission from bilayer 2D GaN occurs at a similar wavelength to its bulk counterpart, due to the cancellation of the effect of quantum confinement on the optical gap by the quantum-confined Stark shift. Polarized light emission at room temperature is possible via uniaxial in-plane strain, which is desirable for energy-efficient display applications. We compare the electronic and optical properties of freestanding two-dimensional GaN to atomically thin GaN wells embedded within AlN barriers in order to understand how the functional properties are influenced by the presence of barriers. Our results provide microscopic understanding of the electronic and optical characteristics of GaN at the few-layer regime. KEYWORDS: GaN, electronic structure, optical properties, first-principles calculations, twodimensional materials

Gallium nitride (GaN) is an important commercial semiconductor for optoelectronic applications in the visible and near-ultraviolet part of the spectrum.1,2 Many efforts directed at producing deep-ultraviolet light with III-nitrides utilize high-Al content AlGaN alloys grown along the polar c-axis crystallographic direction.3 However, the low light-extraction efficiency

ACS Paragon Plus Environment

2

Page 3 of 19

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

Nano Letters

and the lack of deep-UV-transparent p-type contacts limit the overall efficiency of deep-UV AlGaN LEDs.4 On the other hand, devices employing atomically thin GaN quantum wells and tunnel injection achieve deep-UV light emission5 in the desirable spectral range for sterilization applications.6 These atomically thin wells also exhibit strong excitonic effects7 even at room temperature that may improve the internal quantum efficiency.8 The suppression of the quantumconfined Stark shift7 improves the internal quantum efficiency as well9 and stabilizes the emission wavelength as a function of power.10 In addition to atomically thin quantum wells, freestanding two-dimensional gallium nitride (2D GaN) is also expected to demonstrate desirable properties for optoelectronic applications. Unlike other 2D materials such as graphene or MoS2, bulk GaN crystallizes in the non-layered wurtzite structure11 and is not expected to adopt a 2D form with traditional exfoliation methods. However, 2D GaN has recently been directly synthesized experimentally via graphene encapsulation.12 Subsequent theoretical calculations have confirmed the stability and investigated the properties of 2D GaN, as well as other III-V semiconductors.13-15 Nevertheless, the fundamental electronic and optical properties of freestanding 2D GaN remain unexplored. In this work we investigate the electronic and optical properties of monolayer and bilayer freestanding 2D GaN. We apply first-principles calculations based on density functional (DFT) and many-body perturbation theory to predict accurate electronic band gaps, luminescence energies, and excitonic properties of monolayer and bilayer 2D GaN as a function of strain. Our calculations reveal that the strong inherent polarization field present in 2D GaN amplifies the quantum-confined Stark effect and partially cancels out the effects of confinement on the electronic and optical gaps. As a result, light emission from bilayer 2D GaN occurs at a similar

ACS Paragon Plus Environment

3

Nano 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

Page 4 of 19

wavelength as its bulk GaN counterpart. The exciton binding energies are much larger than bulk GaN due to increased electron-hole interaction caused by confinement. Monolayer GaN emits in the deep-UV range, while uniaxial in-plane strain allows for linearly polarized light emission at room temperature in both structures. Our results suggest that 2D GaN is a promising material for nonlinear optics, energy-efficient display applications, and germicidal and water-purification processes. Our first-principles calculations are based on DFT and many-body perturbation theory. DFT structural relaxation calculations were performed using the local density approximation (LDA) for the exchange-correlation functional16,17 within Quantum ESPRESSO18 using a planewave basis and norm-conserving pseudopotentials. The 4s, 4p and 3d electrons of Ga were included in the valence. The relaxed in-plane lattice constant is in good agreement with experimental values for bulk GaN (3.19 Å)19, being underestimated by only 0.77%. Hydrogen atoms of integer charge were included to passivate surface dangling bonds. In simulations, it is common to passivate surface dangling bonds of the polar surfaces of compound semiconductors with fictitious hydrogen atoms of partial charge to eliminate charge transfer.20 However, in our simulations we are interested in examining the physical effects of hydrogen passivation, including charge transfer to 2D GaN, and hence we use hydrogen pseudopotentials of integer charge. Furthermore, our choice of pseudopotential is supported by our formation-energy calculations, which confirm that full hydrogen passivation of surface bonds is energetically favorable for the monolayer. A buckled geometry of 2D GaN was investigated, as passivation makes it more stable than the planar variant.12 Phonon calculations with density functional perturbation theory21 and the frozen phonon method22 (supplementary) were used to explore the mechanical stability of the investigated structures (Figure S1). DFT calculations were carried out

ACS Paragon Plus Environment

4

Page 5 of 19

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

Nano Letters

with an 8 x 8 x 1 Monkhorst-Pack mesh. To obtain converged band gap values, an artificial dipole was included in the vacuum space between periodic slabs. The technical details are analyzed in previous work23. The magnitude of the dipole is chosen to cancel out the artificial electric field in the vacuum region between periodic images of the slabs [i.e., the slopes of the corrected electrostatic potential curves become zero away from the slab, Figures 1(c) and 1(d)]. The band structures were calculated with the G0W0 method24 as implemented in the BerkeleyGW package.25 The GW calculations employed semicore norm-conserving pseudopotentials which consider Ga 3s, 3p, and 3d orbitals to be valence states.26 A plane wave cutoff energy of 250 Ry converges semicore pseudopotential DFT eigenvalues to within 20 meV.

To eliminate

interaction between periodic images, a Wigner-Seitz slab truncation was applied to the Coulomb interaction.27 Simulation cell volumes included vacuum space oriented perpendicularly to the slabs such that 99% of charge density was contained within half the cell. Also employed in the GW method were the Hybertsen-Louie generalized plasmon-pole model28 for the dielectric response’s frequency dependence and the static remainder approach29 to speed up summation convergence over unoccupied states. Quasiparticle band structure calculations were performed using a converged screening cutoff energy of 34 Ry and summing over unoccupied states with energy up to 50% of the screening cutoff. Sampling grids included 8 x 8 x 1, 10 x 10 x 1, 12 x 12 x 1, and 16 x 16 x 1 Monkhorst-Pack meshes, and converged values were obtained from extrapolation (Figure S2) in a manner consistent with previous work.30 Quasiparticle band structures were determined by modifying the DFT-level band energies with a scissors-shift operator constructed from GW quasiparticle corrections (Figures S3 and S4). Finally, exciton binding energies were calculated using the Bethe-Salpeter equation (BSE) method. The top three valence bands and the lowest conduction band were included in this calculation.

ACS Paragon Plus Environment

5

Nano 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

Page 6 of 19

We performed formation-energy calculations to understand the energetics of formation and H-passivation of 2D GaN. Our results agree with Al Balushi et al. that the hydrogenated buckled structure is the most stable form of 2D GaN.12 Our calculated enthalpy of hydrogenation per H2 molecule is -1.102 eV for the monolayer and -1.551 eV for the bilayer relative to the unpassivated buckled structures. Moreover, the formation enthalpy of H-passivated 2D GaN with respect to bulk GaN and the H2 molecule is -109 meV for monolayer GaN and 82 meV for bilayer GaN. H-passivated monolayer GaN is therefore a thermodynamically stable compound, while the H-passivated GaN bilayer is a metastable structure. The interfacial polarization from opposing gallium and nitrogen-terminated surfaces gives rise to a strong inherent electric field perpendicular to the 2D layer, which works against the effect of quantum confinement and reduces the band gap. The magnitude of the electric field is calculated to be 74 MV/cm in monolayer GaN, and 65 MV/cm in bilayer GaN, as determined by the slope of the plane-averaged electrostatic potential between the H-atom minima (Figure 1c,d). The electric-field values determined from the dipole moment of the structure along the perpendicular axis are of similar order (details in supplementary information). These electricfield values are remarkably larger than the breakdown voltage of bulk GaN (3.3 MV/cm)31, but they do not cause band-gap closure or impact ionization, since the slab thicknesses are sufficiently small such that carriers cannot be accelerated to high energies. On the other hand, these polarization fields introduce a slope to the energy bands and separate electrons from holes, similar to strained nitride quantum wells9. The polarization fields separate the electron and hole wave functions (Figure 1a,b) and lower the gap (quantum confined Stark effect), counteracting the band-gap increase caused by quantum confinement. The lack of inversion symmetry further enables nonlinear optical properties.32

ACS Paragon Plus Environment

6

Page 7 of 19

VBM CBM

(b)

2

|ψ 2 | |ψ |

VBM (a) CBM

1 0.5

Potential Energy (Ry) Potential Energy (Ry)

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

Nano Letters

(c)

0.5

(d)

0 0 -0.5

-0.5

-1

-1

-1.5

-2 0

-1.5

2

4

6

8

c-AxisPosition Position (Å)(Å) c-Axis

10

12

-2 0

2

4

6

8

10

12

14

c-Axis Position (Å) (Å) c-Axis Position

Figure 1. Electron and hole wave functions squared for (a) monolayer and (b) bilayer 2D GaN, and plane-averaged electrostatic potentials for (c) monolayer and (d) bilayer 2D GaN. The polarization fields within each structure spatially separate electrons and holes. The artificial dipole correction in the vacuum region is necessary to cancel out artifacts of the periodic boundary conditions. The distance and slope between peaks of like atoms – indicated by the dashed line in (c) – allow for the calculation of the internal electric field magnitude. The quasiparticle band structures of 2D GaN are shown in Figure 2. Both the valence band maximum (VBM) and conduction band minimum (CBM) occur at Γ for both structures. The VBM is primarily composed of N 2p orbitals while the CBM is a hybridization of Ga 4s, Ga 4p, N 2s, and H 1s states in both structures. The DFT band gap of monolayer GaN is 2.95 eV, which is larger than the DFT gap of bulk GaN (1.79 eV) due to quantum confinement. However,

ACS Paragon Plus Environment

7

Nano Letters

bilayer GaN exhibits a DFT band gap of 1.32 eV that is surprisingly lower than the bulk value. The reason for the gap reduction in bilayer GaN compared to bulk is the Stark effect due to the strong polarization field. However, DFT/LDA significantly underestimates semiconductor band gaps, hence GW corrections are needed to obtain accurate gap values. The GW band gap of monolayer GaN was found to be 6.32 eV, and the band gap of bilayer GaN is 4.00 eV (Table 1). The converged band-gap values are listed in Table 1 and compared to atomically thin GaN wells surrounded by AlN barriers, which were presented in previous work.7 The disparity between the freestanding and the AlN-embedded 2D GaN is due to the polarization fields. As the polarization fields in freestanding 2D GaN are not screened by the surrounding vacuum, the gap is more strongly affected by the quantum-confined Stark effect. While the electric field is stronger in the monolayer than the bilayer structure, the small thickness of the slab inhibits the spatial separation of electrons and holes, hence the Stark effect on the gap is small.

Energy (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 19

10

10

8

8

6

(a)

6

4

4

2

2

0

0

-2

-2

-4

-4

-6

K

Γ

(b)

-6

M K

Γ

M

Figure 2. Quasiparticle band structures of (a) monolayer and (b) bilayer 2D GaN. Both materials exhibit a direct gap at Γ. The polarization fields counteract the effects of confinement in the bilayer structure.

ACS Paragon Plus Environment

8

Page 9 of 19

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

Nano Letters

Monolayer

Bilayer

Freestanding

Embedded7

Freestanding

Embedded7

Band Gap (eV)

6.32

4.76-5.44

4.00

4.55-4.69

Exciton Binding Energy (eV)

1.31

0.10-0.21

0.75

0.10-0.17

Optical Gap (eV)

5.01

4.66-5.23

3.25

4.45-4.52

Table 1. Band Gaps, Binding Energies, and Optical Gaps of 2D GaN. The freestanding GaN values are obtained from extrapolation (supplementary). Embedded values are from GaN wells within AlN barriers.7 Excitonic effects are also much stronger in 2D GaN compared to the bulk counterpart. The lowest-exciton binding energy in monolayer GaN is 1.31 eV, whereas in bilayer GaN the binding energy is 0.75 eV. Both values are more than one order of magnitude larger than bulk GaN (0.02 eV)33, due to the increased electron-hole interaction strength caused by quantum confinement. Subtracting the exciton binding energies from the band gap yields a lowest-exciton (i.e., luminescence) energy of 5.01 eV for monolayer GaN, which is comparable to the value of AlN-embedded monolayer-GaN quantum wells. The lowest-exciton energy in bilayer GaN is 3.25 eV, which is slightly smaller than the luminescence energy of bulk GaN. Our results therefore demonstrate that the effects of confinement and Stark effect on the carriers and excitons cancel each other out in bilayer GaN. For thicker structures we anticipate that the Starkeffect contribution dominates, and the luminescence energy is lower than the bulk. However, in these thicker structures the polarization fields are more effectively screened by free carriers, and the quantum-confined Stark effect is suppressed. We further evaluated the radiative lifetimes of excitons in monolayer and bilayer 2D GaN. We first determined the singlet-triplet exchange splitting energies to be 63 meV for

ACS Paragon Plus Environment

9

Nano 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

Page 10 of 19

monolayer GaN and 6 meV for bilayer GaN. The monolayer value is large compared to kBT at room temperature and below, hence excitons in the monolayer predominantly occupy long-lived triplet states with long diffusion lengths. The small exchange splitting for the bilayer implies a large thermal occupation of singlets and shorter radiative lifetimes and diffusion lengths. Following the methodology developed for evaluating radiative exciton lifetimes in twodimensional materials34, the effective radiative lifetime τ  at room temperature is calculated as the thermal average of the radiative lifetimes τ  of exciton states  according to: τ  =

and

     ℏ  

τ  =



  

 



 

') *

∑& '( +, ') *

∑ +, -

,



!" #,

(1)

(2)

where . is the exciton effective mass, /0 is the area of the unit cell, 1 is the square modulus of the BSE exciton transition dipole divided by the number of 2D k-points, and 2 0 is the exciton energy calculated using the BSE method. We assumed the effective mass of electrons and holes to be the same as in bulk GaN (0.2 me and 1.4 me, respectively). The BSE calculation includes the top three valence bands and the lowest conduction band, while the Brillouin zone was sampled with 8 x 8 x 1, 10 x 10 x 1, 12 x 12 x 1, and 16 x 16 x 1 Monkhorst-Pack meshes. Converged values were obtained from extrapolation to an infinitely dense grid. We determined the thermally averaged radiative lifetimes of singlet excitons at 300 K to be 0.6 ns and 4.6 ns for monolayer and bilayer GaN, respectively. These values are comparable to those of other 2D semiconductors34 and approximately one order of magnitude shorter than typical values in InGaN LEDs for typical free-carrier densities (5×1018 cm–3)8. The shorter exciton lifetimes in 2D GaN may be beneficial for improved internal quantum efficiency, as well as telecommunication (Li-Fi) applications.

ACS Paragon Plus Environment

10

Page 11 of 19

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

Nano Letters

Since in-plane strain is a viable method of tuning the electronic properties of atomicallythin materials35,36 we also explore the effect of strain on the luminescence properties of 2D GaN. We applied both compressive and tensile in-plane biaxial strain of up to 5%. Our results (band gap calculated with DFT and adjusted to account for the zero-strain GW and BSE corrections, Figure 3) show that increasing tensile strain first increases the band gap, and subsequently leads to a small decline. On the other hand, compressive strain monotonically reduces the gap by up to 0.3 eV for a 5% strain. To understand these trends, we examined the electron and hole wave functions as a function of strain (Figure S5). Electrons and holes are spatially separated due to the inherent polarization field. While the hole wave function is fully confined inside the slab, the electron wave function noticeably extends beyond the edge of the slab and is only weakly bound within the potential well of the structure. Tensile in-plane strain reduces the thickness of the slab and hence the electron-hole separation, while compressive strain further shifts the electron probability to the surface of the slab, amplifying the Stark effect and the gap reduction (supplementary). The gap decrease at large tensile strains likely results from the interplay between the quantum-confined Stark shift and the gap reduction by the deformation potential. The crystal-field splitting for ±5% biaxial strain spans a range of 1.28-2.80 eV and 0.80-1.63 eV for monolayer GaN and bilayer GaN, respectively.

ACS Paragon Plus Environment

11

Nano 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

Page 12 of 19

Figure 3. Optical gap of monolayer and bilayer 2D GaN as a function of biaxial strain. Compressive strain is a viable method for controlling the band gap of both structures. These results were determined from the DFT gap as a function of strain adjusted to account for the GW and BSE corrections at zero strain. We also investigated the effect of uniaxial compressive and tensile in-plane strain of up to 5% on the luminescence properties. The uniaxial strain breaks the in-plane symmetry and lifts the degeneracy of the top two N 2p valence bands, resulting in optical polarization. The in-plane Poisson ratio is 0.35 for compressive strain and 0.1 for tensile strain along the directions denoted in Figure 4. The band-gap values in Figure 4 are estimated using DFT. For the unstrained monolayer GaN GW corrections increase the crystal-field splitting (i.e., the energy difference

ACS Paragon Plus Environment

12

Page 13 of 19

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

Nano Letters

between the topmost two degenerate valence bands and the third from the top at Γ) from the DFT value of 1.837 eV to 2.360 eV, i.e., an increase of 22%. Since the topmost three valence bands of GaN are of the same N 2p orbital character, we expect that quasiparticle corrections will similarly increase the DFT degeneracy splitting due to strain. Compressive strain results in a larger heavy and light hole splitting than tensile strain, but even a tensile strain of 1% is sufficient to split the bands by more than kBT at room temperature (Figure 4). The resulting emitted light is linearly polarized along the direction of compressive strain. Polarized light emission is also evident from the anisotropic optical absorbance spectra of uniaxially strained 2D GaN (Figure S6). Polarized light emission at room temperature is therefore possible with 2D GaN, which is desirable for energy-efficient display applications. The emission wavelength of the polarized luminescence can be further adjusted in the visible range by alloying with InN to form 2D InGaN.

ACS Paragon Plus Environment

13

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

Valence Band Splitting (eV)

Nano Letters

0.4 0.3 0.2 0.1 0 -0.1 -0.2 -0.3 -0.4 -0.5 -0.6 -5

Page 14 of 19

1 ML 2 ML

1 ML 2 ML -4

-3

-2

-1

0

1

2

3

4

5

Uniaxial Strain (%) Figure 4. Energy splitting between the heavy and light holes of monolayer and bilayer GaN as a function of uniaxial strain. Uniaxial strain lifts the band degeneracy, resulting in polarized light emission along the axis of compressive strain. In conclusion, we investigated the electronic and optical properties of 2D GaN as a function of thickness and strain with predictive calculations. Monolayer 2D GaN emits light in the deep-UV range, which is promising for sterilization applications. The long-lived stable triplet excitons of the monolayer may be promising for excitonic applications. Uniaxial in-plane strain results in linearly polarized light emission desirable for display applications. Our results demonstrate that 2D GaN exhibits an array of desirable functional properties and is a synergistic compound between established semiconductors and 2D materials.

ACS Paragon Plus Environment

14

Page 15 of 19

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

Nano Letters

ASSOCIATED CONTENT Supporting Information. The following files are available free of charge. Mechanical stability, details of band gap and exciton binding energy calculations, GW corrections to DFT/LDA band structures, wave functions as a function of strain. (PDF) AUTHOR INFORMATION Corresponding Author *E-mail: [email protected] Notes The authors declare no competing financial interest. ACKNOWLEDGMENT This work was supported by the NSF ECCS-CDS&E program (1607796). Computational resources were provided by the DOE NERSC facility under Contract No. DE-AC0205CH11231. Graphics were generated with VESTA.37 REFERENCES (1) Pimputkar, S.; Speck, J. S.; DenBaars, S. P.; Nakamura, S. Nat. Photon. 2009, 3, 180-182. (2) Kioupakis, E.; Rinke, P.; Janotti, A.; Yan, Q.; Van de Walle, C. G. Energy Conversion: Solid-State Lighting. In Computational Approaches to Energy Materials; Walsh, A., Sokol, A. A., Catlow, C. R. A., Eds.; John Wiley & Sons Ltd.: Oxford, UK, 2013; p 231−259. (3) Taniyasu, Y.; Kasu, M.; & Makimoto, T. Nature 2006, 441, 325-328.

ACS Paragon Plus Environment

15

Nano 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

Page 16 of 19

(4) Kneissl, M.; Kolbe, T.; Chua, C.; Kueller, V.; Lobo, N.; Stellmach, J.; Knauer, A.; Rodriguez, H.; Einfeldt, S.; Yang, Z.; Johnson, N.M.; Weyers, M. Semicond. Sci. Technol. 2010, 26, 014036. (5) Verma, J.; Kandaswamy, P. K.; Protasenko, V.; Verma, A.; Grace Xing, H.; Jena, D. Appl. Phys. Lett. 2013, 102, 041103. (6) Jagger, J. Introduction to research in ultraviolet photobiology. Prentice-Hall. Englewood Cliffs, NJ. 1967. (7) Bayerl, D.; Islam, S. M.; Jones, C. M.; Protasenko, V.; Jena, D.; Kioupakis, E. Appl. Phys. Lett. 2016, 109, 241102. (8) Langer, T.; Chernikov, A.; Kalincev, D.; Gerhard, M.; Bremers, H.; Rossow, U., Koch M.; Hangleiter, A. Appl. Phys. Lett. 2013, 103, 202106. (9) Kioupakis, E.; Yan, Q.; Van de Walle, C. G. Appl. Phys. Lett. 2012, 101, 231107. (10) Speck, J. S.; Chichibu, S. F. MRS Bull. 2009, 34, 304-312. (11) Dingle, R.; Sell, D. D.; Stokowski, S. E.; Ilegems, M. Phys. Rev. B 1971, 4, 1211-1218. (12) Al Balushi, Z. Y.; Wang, K.; Ghosh, R. K.; Vilá, R. A.; Eichfeld, S. M.; Caldwell, J. D.; Qin, X.; Lin, Y.-C.; DeSario, P. A.; Stone, G.; Subramanian, S.; Paul, D. F.; Wallace, R. M.; Datta, S.; Redwing, J. M.; Robinson, J. A. Nat. Mater. 2016, 15, 1166–1171. (13) Lucking, M.C.; Xie, W.; Choe, D. H.; West, D.; Lu, T.M.; Zhang, S.B. 2016 arXiv preprint arXiv:1611.00121.

ACS Paragon Plus Environment

16

Page 17 of 19

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

Nano Letters

(14) Zhang, H.; Meng, F. S.; Wu, Y. B. Solid State Commun. 2017, 250, 18-22. (15) Prete, M. S.; Mosca Conte, A.; Gori, P.; Bechstedt, F.; Pulci, O. Appl. Phys. Lett. 2017, 110, 012103. (16) Ceperley, D. M.; Alder, B. J. Phys. Rev. Lett. 1980, 45, 566-569. (17) Perdew, J. P.; Zunger, A. Phys. Rev. B 1981, 23, 5048-5079. (18) Giannozzi, P.; Baroni, S.; Bonini, N.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.; Chiarotti, G. L.; Cococcioni, M.; Dabo, I.; Dal Corso, A.; de Gironcoli, S.; Fabris, S.; Fratesi, G.; Gebauer, R.; Gerstmann, U.; Gougoussis, C.; Kokalj, A.; Lazzeri, M.; Martin-Samos, L.; Marzari, N.; Mauri, F.; Mazzarello, R.; Paolini, S.; Pasquarello, A.; Paulatto, L.; Sbraccia, C.; Scandolo, S.; Sclauzero, G.; Seitsonen, A. P.; Smogunov, A.; Umari, P.; Wentzcovitch, R. M. J. Phys. Condens. Matter 2009, 21, 395502. (19) Schulz, H.; Thiemann K.H. Solid State Commun. 1977, 23, 815-819. (20) Huang, X.; Lindgren, E.; Chelikowsky, J.R. Phys. Rev. B 2005, 71, 165328. (21) Baroni, S.; De Gironcoli, S.; Dal Corso, A.; Giannozzi, P. Rev. Mod. Phys. 2001, 73, 515562. (22) Togo, A.; Tanaka, I. Scr. Mater. 2015, 108, 1-5. (23) Bengtsson, L. Phys. Rev. B 1999, 59, 12301-12304. (24) Hybertsen, M. S.; & Louie, S. G. Phys. Rev. B 1986, 34, 5390-5413.

ACS Paragon Plus Environment

17

Nano 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

Page 18 of 19

(25) Deslippe, J.; Samsonidze, G.; Strubbe, D. A.; Jain, M.; Cohen, M. L.; Louie, S. G. Comput. Phys. Commun. 2012, 183, 1269-1289. (26) Tiago, M. L.; Ismail-Beigi, S.; Louie, S. G. Phys. Rev. B 2004, 69, 125212. (27) Ismail-Beigi, S. Phys. Rev. B 2006, 73, 233103. (28) Hybertsen, M.S.; Louie, S.G. Phys. Rev. B 1986, 34, 5390-5413. (29) Deslippe, J.; Samsonidze, G.; Jain, M.; Cohen, M. L.; Louie, S. G. Phys. Rev. B 2013, 87, 165124. (30) Shi, G.; Kioupakis, E. Nano Lett. 2015, 15, 6926-6931. (31) Dobrinsky, A.; Simin, G.; Gaska, R.; Shur, M. ECS Trans. 2013, 58, 129-143. (32) Rondinelli, J. M.; Kioupakis, E. Annu. Rev. Mater. Res. 2015, 45, 491-518. (33) Muth, J. F.; Lee, J. H.; Shmagin, I. K.; Kolbas, R. M.; Casey Jr, H. C.; Keller, B. P.; Mishra, U.K.; DenBaars, S. P. Appl. Phys. Lett. 1997, 71, 2572-2574. (34) Palummo, M.; Bernardi, M.; Grossman, J. C. Nano Lett. 2015, 15, 2794-2800. (35) Huang, Z.; Lü, T. Y.; Wang, H. Q.; Yang, S. W.; Zheng, J. C. Comput. Mater. Sci. 2017, 130, 232-241. (36) Rodin, A. S.; Carvalho, A.; Neto, A. C. Phys. Rev. Lett. 2014, 112, 176801. (37) Momma, K.; Izumi, F. J. Appl. Crystallogr. 2011, 44, 1272-1276.

ACS Paragon Plus Environment

18

Page 19 of 19

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

Nano Letters

For Table of Contents only

ACS Paragon Plus Environment

19