New Assembly-Free Bulk Layered Inorganic Vertical Heterostructures

Dec 10, 2018 - Mergers and acquisitions in the drug industry last year brought one particularly high price tag: Takeda... POLICY CONCENTRATES ...
0 downloads 0 Views 1MB Size
Subscriber access provided by University of Winnipeg Library

Communication

New Assembly-Free Bulk Layered Inorganic Vertical Heterostructures with Infrared and Optical Bandgaps Evan R Antoniuk, Gowoon Cheon, Aditi Krishnapriyan, Daniel A. Rehn, Yao Zhou, and Evan J. Reed Nano Lett., Just Accepted Manuscript • DOI: 10.1021/acs.nanolett.8b03500 • Publication Date (Web): 10 Dec 2018 Downloaded from http://pubs.acs.org on December 12, 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 28 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

New Assembly-Free Bulk Layered Inorganic Vertical Heterostructures with Infrared and Optical Bandgaps Evan R. Antoniuk†, Gowoon Cheon‡, Aditi Krishnapriyan§, Daniel A. Rehn∥, Yao Zhou§, and Evan J. Reed §* *E-mail: [email protected]. Phone: (650) 723-2971. Fax: (650) 725-4034. †Department

of Chemistry, ‡Department of Applied Physics, §Department of Materials

Science and Engineering, and ∥Department of Mechanical Engineering, Stanford University, Stanford, California 94305, United States

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 28

ABSTRACT

In principle, a nearly endless number of unique van der Waals heterostructures can be created through the vertical stacking of two-dimensional (2D) materials, resulting in unprecedented potential for material design. However, this widely employed synthetic method for generating van der Waals heterostructures is slow, imprecise, and prone to introducing interlayer contaminants when compared with synthesis methods that are scalable to industrially relevant scales. Herein, we study the properties of a new class of layered bulk inorganic materials that has recently been reported, which we call assemblyfree bulk layered inorganic heterostructures, wherein the individual layers are of dissimilar chemical composition, distinguishing them from commonly studied layered materials. We find that these bulk materials exhibit properties similar to vertical heterostructures, but without the complex and unscalable stacking process. Using stateof-the-art computational approaches, we study the electronic properties of livingstonite (HgSb4S8), a naturally occurring mineral that is a bulk lattice-commensurate heterostructure. We find that isolated bilayers of livingstonite have an intralayer HSE-06 band gap of 2.08eV. This is the first report of a naturally occurring van der Waals heterostructure with a calculated band gap in the visible spectrum. We also studied the electronic properties of tetragonal Ti3Bi4O12, Sm2Ti3Bi2O12, orthorhombic Ti3Bi4O12, Nb3Bi5O15, LaTiNbBi2O9 and AgPbBrO and found some of them are potentially well suited for photovoltaic applications. We also provide characterization of the electronic structure of the isolated bilayer and monolayer subcomponents of the bulk

ACS Paragon Plus Environment

2

Page 3 of 28 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

heterostructures. The report of the properties of these materials significantly enhances the library of known van der Waals heterostructures.

KEYWORDS two-dimensional materials; van der Waals heterostructures; density functional theory; layered materials; band gap; heterojunctions TEXT As a result of the ever growing library of two-dimensional (2D) materials, vertical stacks of chemically dissimilar, weakly-bound 2D layered materials known as van der Waals heterostructures have recently garnered significant interest. In particular, these materials allow for control over the chemical composition of the 2D layers and the stacking sequence in the van der Waals heterostructures, which in turn allows for the precise tailoring of the properties of the heterostructure for a specific application1. Furthermore, since the 2D layers are held together through weak van der Waals interactions, these heterostructures can be fabricated without any lattice matching restrictions2. As such, van der Waals heterostructures have been employed in devices including LEDs, photovoltaics, and 2D superconductors.3–5 Owing to their relatively large absorption coefficient, tunable band gap, and ultrathin width, monolayers of 2D transition metal dichalcogenides (TMDs) have widely been explored for use in many photonic applications.6

As a result, recent efforts have

attempted to study photogenerated interlayer excitons in heterostructures composed of

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 28

TMDs including MoS2/WSe2, MoSe2/WSe2, and MoS2/WS2/MoSe2.7–9 The interlayer exciton lifetime in these materials has typically been measured to be on the order of nanoseconds and is substantially lengthened by the spatial separation of the photogenerated electron-hole pairs.9 On the other hand, electrons have been shown to be very quickly and efficiently shuttled between layers in heterostructures. For example, an electron transfer time of 1.5ps has been measured in a trilayer MoS2/WS2/MoSe2 heterostructure.7 Accordingly, these materials have been proposed to be candidates for ultrathin photovoltaic devices.4 However, studies of the charge transfer processes in van der Waals heterostructures have largely been limited to heterostructures composed of monolayers of boron nitride, graphene, or TMDs. Despite the numerous attractive properties of van der Waals heterostructures, their experimental synthesis has faced challenges. One possible synthetic method involves manually stacking individual 2D layers obtained from the mechanical exfoliation of several different bulk layered materials.10

However, this method is susceptible to

contamination as adsorbates attached to the surface of the 2D layers can become trapped in the interface between the two layers. Additionally, the process of manually stacking the layers is arduous and precisely aligning the 2D crystals is challenging. As such, this method faces challenges for large-scale synthesis.11 An alternative method involves the chemical vapor deposition (CVD) growth of layers on top of each other.12 However, due to the restrictive temperature and pressure conditions required for the CVD growth of each 2D crystal, this method has considerable limitations in the variety of heterostructures that it can produce. Although a one-step CVD growth method has been

ACS Paragon Plus Environment

4

Page 5 of 28 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

reported for the synthesis of vertical stacks of MoS2/WS2 heterostructures, this method has limited control of the size of the heterostructures produced.13 Molecular beam epitaxy approaches have also been employed to generate high quality and well oriented vertically stacked heterostructures. Using this approach, atomically thin van der Waals heterostructures of MoSe2/HfSe2, MoTe2/MoS2, SnSe2/WSe2

and hexagonal boron

nitride/graphene have been produced.14–17 Layered hybrid organic-inorganic materials have also been synthesized with solution fabrication methods.18 For example, solution processing methods have been used to provide an easy and low-cost method for the production of hybrid perovskite solar cells, but have challenges in reliably reproducing large-area films and are generally expected to exhibit lower mobilities than entirely inorganic heterostructures.19 An alternative approach for synthesizing inorganic van der Waals heterostructures involves mechanically exfoliating a bulk crystal that is composed of vertical stacks of chemically dissimilar layers. Since this approach does not require any manual stacking of layers, complications related to interlayer contamination and the misalignment of the layers are avoided. Additionally, this process is considerably less labor-intensive than manually stacking the layers. By exfoliating the bulk heterostructures down to single monolayers, this approach may also serve as a simple means of synthesizing for the previously inaccessible 2D materials. To further simplify the synthetic process, van der Waals heterostructures have been created from the mechanical exfoliation of bulk naturally occurring minerals. In particular, the bulk mineral franckeite, whose structure consists of alternating layers of SnS2 and PbS, has been successfully exfoliated to yield a

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

few-layer van der Waals heterostructure.20,21

Page 6 of 28

Furthermore, mechanically exfoliated

franckeite flakes were then implemented into a near-IR photodetector device, demonstrating the potential applications of these materials.20,21 The use of a naturally occurring mineral eliminates the need for any synthesis of artificial materials while also ensuring the air stability of the bulk crystal. Perhaps even more importantly, in the absence of exfoliation, the bulk materials consist of superlattice heterostructures that may be amenable to photonic and other applications involving charge carrier separation or recombination. The interfacial area in these materials scales with the volume of material, rather than the two-thirds power of volume that one might expect for a single interface. Currently, the main limitation of this approach is the small number of known, bulk, layered crystals with dissimilar layers. In earlier work, we implemented a data mining approach to identify weakly bonded layered materials found within the publicly available databases of tens of thousands of experimentally characterized crystals.22 In this earlier work, we identified over 1000 layered materials, 487 materials consisting of weakly bonded one-dimensional or molecular wire subcomponents, and 98 bulk crystals composed of heterostructures of chemically dissimilar 1D and 2D subcomponents from the Materials Project database of mostly experimentally characterized crystals.22,23 The complete list of materials can be found in the Supporting Information in Reference 22. In the present work, we employ Kohn-Sham density functional theory to investigate the feasibility of generating fewlayer van der Waals heterostructures from these newly uncovered bulk crystals. Furthermore, we explore the electronic structure of the heterostructure bilayers and each

ACS Paragon Plus Environment

6

Page 7 of 28 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

of its subcomponent monolayers to provide the first characterization of these few-layer materials. While conventional vertical heterostructures require challenging fabrication techniques that are difficult to scale, we find that these materials exhibit the properties of vertical heterostructures in their bulk form, eliminating the need for complex fabrication techniques of nanostructures altogether. Hence, we refer to these bulk materials as assembly-free vertical heterostructures.

Figure 1. Structures of the layered, bulk, lattice-commensurate heterostructures studied in this work and their corresponding Materials Project ID: (a) Tetragonal Ti3Bi4O12, Mp2333524 (b) Sm2Ti3Bi2O12, Mp-55743825 (c) Orthorhombic Ti3Bi4O12, Mp-62219826 (d) Nb3Bi5O15, Mp-64559227 (e) LaTiNbBi2O9, Mp-63665428 (f) AgPbBrO, Mp-55947029 and (g) HgSb4S8, Mp-54259630. All of the structures shown are DFT relaxed versions of experimentally characterized, bulk structures as given on the Materials Project database. Structures (a), (b), (c), and (e) are all Aurivillius phase perovskites, which are structures consisting of alternating layers of [Bi2O2]2+ and a perovskite-like layer with the general formula [An-1BnO3n+1]2-. Compound (g) is the bulk structure of the naturally occurring mineral, livingstonite. Calculated interlayer distances and their experimental comparisons can be found in the SI.

ACS Paragon Plus Environment

7

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 8 of 28

To identify the most promising heterostructures for study with hybrid DFT calculations, we apply a number of filters to shrink the list of 98 bulk heterostructures. First, we limit our study to only those bulk structures with direct experimental report of their synthesis and structural characterization. Next, we include only heterostructures that are composed of layers of 2D subunits. Finally, we limit this study to materials with a calculated bulk PBE band gap between 1-3eV to increase the likelihood that bilayers of these bulk crystals will also have band gaps in the visible or IR spectrum for photonic applications. Although it has a PBE calculated bulk band gap of 0.537eV in the Materials Project database, we also include AgPbBrO as AgBr/PbO systems have been reported to be studied for use in photodiodes.31 Furthermore, a AgBr/ZnO heterojunction has been reported to exhibit considerable photocatalytic activity for the degradation of rhodamine B, an environmentally harmful and non-biodegradable organic dye used in pigmentation industries.32 Applying these requirements reduces the number of bulk heterostructures down to 7 (Figure 1), which is tractable to be investigated with hybrid DFT calculations. Notably, this list includes livingstonite, a naturally occurring bulk mineral whose structure consists of alternating layers of HgSb2S4 and SbS2.33 Additionally, 4 of these compounds belong to a class of materials known as Aurivillius phase oxides.34–36 These compounds are given by the general formula Bi2An-1BnO3n+3 and are comprised of alternating layers of [Bi2O2]2+ and the perovskite-like layer of [An1BnO3n+1]

2-,

where n is the number of layers.37 These materials have previously been

extensively studied for their ferroelectric38 and oxide ion conductivity39 properties. Additionally, Aurivillius phase perovskites have been explored as visible-light

ACS Paragon Plus Environment

8

Page 9 of 28 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

photocatalysts for the degradation of environmental pollutants.40

Notably, the

Aurivillius phase perovskites studied in this work do not contain Pb, alleviating the environmental concerns typically associated with the processing of lead-based materials including perovskite solar cells.41 To determine the ease at which these bulk heterostructures may be mechanically exfoliated down to few-layer systems, we calculate the interlayer cohesive energy (Equations 1-2) with the PBE exchange-correlation functional (Figure 2).

To better

account for van der Waals interactions, these calculations were performed with the DFTD3 correction method.42 As GGA functionals tend to neglect long-range van der Waals interactions, including this correction resulted in a larger calculated interlayer cohesive energy for all compounds. To validate the use of the D3 correction, the calculated interlayer distances are compared to the experimental values and found to have an average deviation of 0.31Å (see SI).

Notably, the interlayer cohesive energy value

calculated with Equation 1-2 is the difference between the energy of the heterostructure and the sum of the energies of the two isolated and stress relaxed monolayers. As such, this value includes some component of stress relaxation energy in addition to interlayer van der Waals interactions. Figure 2 also shows an approximate maximum interlayer cohesive energy of potentially exfoliable layered materials of 150 meV/Å2 for comparison.

This upper bound was chosen to correspond with the vdW-optB8843

calculated interlayer cohesive energy of SnSe44, which has previously been reported to form 2D monolayers.45 We find that isolated bilayers of AgPbBrO and HgSb4S8 are weakly bound by interlayer van der Waals interactions with cohesive energies of 5 meV/

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 28

Å2 and 32meV/Å2, respectively, and are likely to be exfoliable. Interestingly, the more strongly bound materials identified here are Aurivillius phase materials which have the general formula of layers of [Bi2O2]2+ and [An-1BnO3n+1]2-. In order to maintain charge neutrality, these materials must contain interlayer charge transfer. This result is further supported by the calculated Bader charges on the atoms at the heterointerface (see SI, Table S4). Since this class of materials is known to have interlayer charge transfer, our choice of screening criteria (PBE calculated bandgap between 1-3 eV) may have disproportionately increased the number of high interlayer binding energy compounds found relative to lower bandgap materials. Neglecting structural relaxation, the calculated bulk cohesive energy of these materials should be twice that of its bilayer cohesive energy per formula unit if the binding is dominated by nearest neighbor van der Waals interactions since the bulk material possesses twice as many bonded interfaces. The two most weakly bonded materials in Figure 2 (HgSb4S8 and AgPbBrO) approximately obey this condition, suggesting van der Waals dominated bonding. However, the oxides exhibit stronger bilayer bonding relative to the bulk, consistent with an ionic charge transfer character. To determine the mechanical stability of the monolayers generated through the exfoliation of these heterostructures in a freely suspended form, we also calculate some of their phonon spectra and elastic constants. The density functional perturbation theory (DFPT) calculated phonon spectrum of the PbO monolayer contains all positive frequencies and its elastic matrix was found to be positive definite, indicating that the PbO monolayer is likely to be mechanically stable as a freely suspended monolayer. Very

ACS Paragon Plus Environment

10

Page 11 of 28 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

recently, atomic sheets of 𝛼-PbO were reported to be synthesized through a sonochemical exfoliation technique and exhibited a direct, visible band gap and higher thermal stability than monolayer MoS2.46 In addition, these heterostructures offer a feasible route to study the properties of charged monolayers such as BiO2+ which have no bulk counterpart due to charge neutrality constraints.

Figure 2. Calculated interlayer cohesive energy per area of all of the heterostructures considered in this study (Equation 1-2). The values shown here were calculated with the PBE functional with the D3 correction method applied to better account for van der Waals interactions between layers. An approximate upper bound for exfoliable materials is depicted for comparison.44 Due to the presence of glide or mirror planes in the monolayers of all heterostructures, the interlayer cohesive energy values are invariant to the stacking sequence of the layers in the heterostructure bilayer and bulk heterostructures. Calculated energies are with respect to charge neutral and zero stress monolayers. Next, to gain some insight into the electronic structure of these materials, we perform band structure calculations on the bulk heterostructures, the bilayers of the

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 28

heterostructures, and each of the separated monolayers (Figure 3). All of the compounds are composed of two non-metallic monolayers since the bulk materials were chosen to be non-metallic. Furthermore, the majority of the compounds studied have PBE band gaps in the visible region, thus suggesting potential solar photonic applications for these materials.

The PBE functional typically considerably underestimates band gaps,

particularly for oxides. In the case of the AgBr and PbO monolayers, the HSE-06 band gaps are 0.79eV and 0.91eV larger than the PBE calculated band gaps, respectively. As such, a similar increase in the band gaps of the other materials may be expected if the HSE-06 functional was to be used. The difference in the band gaps between the bulk heterostructures and their corresponding bilayers in Nb3Bi5O15, Ti3Bi4O12, and Sm2Ti3Bi2O12 also suggests that these materials may be suitable for band-gap engineering through varying the number of vertically stacked layers.

ACS Paragon Plus Environment

12

Page 13 of 28 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

Figure 3. Calculated quasiparticle band gaps of all relaxed bulk heterostructures, the heterostructure bilayers, and their monolayer subcomponents. The band gaps shown here are calculated with the PBE functional except for AgPbBrO and the HgSb4S8 monolayers, which were calculated with HSE-06. Spin-orbit coupling effects were included in the calculation of all PBE band gaps. The chemical composition of each monolayer is shown above each corresponding bar, whereas the bulk chemical composition is shown beneath the bars. Upon ionic relaxation, the BiO2+ monolayer in orthorhombic Ti3Bi4O12 exists in a different structural phase than in the other heterostructures, accounting for its significantly different band gap.

ACS Paragon Plus Environment

13

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 14 of 28

Figure 4. Partial band structure of the a) HSE-06 AgPbBrO heterostructure and its unrelaxed monolayer subcomponents and b) PBE livingstonite heterostructure. The depictions of the set of high symmetry k-points in reciprocal space can be found in Ref. 65. To allow the band structures to be compared, the vacuum level is set to be 0eV for all band diagrams. a) The PbO and AgBr monolayer band structures shown here were calculated with the structure of the monolayers upon cleavage of the heterostructure bilayer without allowing the monolayer structures to be relaxed to illustrate the effect of interlayer interactions on the band structure. The vacuum level of bulk AgPbBrO is taken to be equal to the vacuum level of the isolated AgPbBrO bilayer. The bulk and bilayer AgPbBrO heterostructure and the PbO monolayer all exhibit a direct band gap at Γ of 1.11eV, 0.85eV and 2.98eV respectively. The AgBr monolayer possesses an indirect band gap of 2.85eV with the valence band maximum located at X and the conduction band minimum located at Γ. Notably, the AgPbBrO bilayer valence band edge closely resembles the PbO monolayer valence band edge, whereas the AgPbBrO bilayer conduction band edge closely resembles the AgBr monolayer conduction band edge. This and the relatively small energy dispersion in the bands along the Γ-Z, interlayer direction in the bulk suggests that the electronic interlayer interactions are relatively weak. The calculated electron effective mass in the interlayer and intralayer directions for bulk AgPbBrO is 0.57me and 0.47me, respectively, while the hole effective mass in the interlayer and intralayer directions is 0.76me and 5.75me. b) The livingstonite band structure is calculated along the monoclinic reciprocal space path: Γ-Y-H-C-E-M1-A-X-Γ-

ACS Paragon Plus Environment

14

Page 15 of 28 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

Z-D-M Z-A D-Y, where Γ-Y is the interlayer direction. The livingstonite heterostructure bilayer and monolayer SbS2 both have an indirect PBE-calculated band gap of 1.12eV and 1.20eV respectively, with the valence band maximum located between Γ and Z and the conduction band minimum located at Γ. Monolayer HgSb2S4 has a direct band gap at Γ of 1.50eV. The calculated electron effective mass in the interlayer and intralayer directions for bulk livingstonite is 0.51me and 0.40me, respectively, while the hole effective mass in the interlayer and intralayer directions is 1.14me and 0.64me. To further characterize the electronic structure of AgPbBrO and HgSb4S8, we perform band alignment calculations to classify the heterojunction as type I (straddling gap), type II (staggered gap), or type III (broken gap). Figure 4 shows the aligned band structure of the AgPbBrO and livingstonite heterostructure bilayers. We approximate the band alignments of weakly bonded heterostructure bilayers determined from their substituent monolayers using the assumption that the weak interlayer interactions will not significantly alter the magnitude of the band gap or the absolute energy of the Fermi level.47 In Figure 4, we can see this approximation is justified for both AgPbBrO and livingstonite bilayers. For AgPbBrO, the bilayer band structure can be well approximated as the superposition of the AgBr and PbO monolayer band structures. In Figure 4b), the livingstonite bilayer valence band and conduction band edge shape and location more closely resembles that of the SbS2 monolayer bands than the HgSb2S4 monolayer bands. Since the position of the VBM and CBM are insensitive to interlayer interactions, livingstonite is probably closest to a type I heterojunction since both the shape and location of the VBM and CBM most closely resemble that of a single monolayer, the SbS2 monolayer. Notably, this is the first report of a naturally occurring bulk heterostructure that can be mechanically exfoliated to yield a type I heterojunction.

Type I

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 28

heterostructures spatially isolate holes and electrons to a single layer which allows for efficient recombination and photon emission. As such, heterostructures of this type have seen use in LED and laser applications.48,49 Figure 5 shows the band alignments of HSE–06 calculations on the monolayers of AgPbBrO and livingstonite. Both HSE-06 and PBE calculations on the relaxed monolayers of livingstonite incorrectly suggest that HgSb4S8 exists as a type II heterojunction with the conduction band edge residing in the SbS2 layer and the valence band edge residing in the HgSb2S4 layer. However, performing the same calculations while constraining the atomic positions and lattice constants of the isolated monolayers to that of the bulk livingstonite heterostructure results in a type I heterojunction shown in Figure 5a). In this case, the relaxation of the monolayers significantly impacts the classification of the heterostructure and results in the misclassification of the livingstonite bilayer.

A

previous report suggests strong optical absorptions in mercury chalcogenide compounds as a result of electronic transitions from the chalcogen s orbitals in the valence band to mercury and chalcogen s orbitals in the conduction band, suggesting strong intralayer absorption.40 Figure 5a) shows AgPbBrO to be a type II heterojunction with the conduction band edge residing in the AgBr layer and the valence band edge residing in the PbO layer. This result agrees with previously proposed band alignments in AgBr/PbO emulsions.31 In the type II band alignment, injected electrons and holes become spatially separated into the conduction band and valence band layers, respectively. Assuming these particles have a sufficiently long lifetime, the electron and hole can then be removed through an

ACS Paragon Plus Environment

16

Page 17 of 28 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

external circuit. As such, heterojunctions of this type allow for charge separation across the two different layers and have been proposed for use in photovoltaics.4

Figure 5. a) HSE-06 calculated band alignment of isolated constituent monolayers of HgSb4S8 (livingstonite) and AgPbBrO constrained to the atomic positions and lattice

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 28

constants of their respective bulk heterostructures. The value located next to each bar represents the energy of the corresponding band edge relative to the vacuum level and the value located next to the arrow is the magnitude of the quasi-particle gap. To align the bands, the band structures of the two monolayers in each heterostructure are shifted such that the calculated vacuum level is 0eV. In agreement with PBE results in Figure 4, these results also suggest that livingstonite and AgPbBrO form type-I and type II heterojunctions, respectively. b) PBE calculated plot of the imaginary part of the KuboGreenwood dielectric function for bulk AgPbBrO, its heterostructure bilayer, and its monolayer substituents. The in-plane component of the dielectric tensor corresponds to either the x or y directions since the dielectric tensor is isotropic in-plane, whereas the out-of-plane component corresponds to the z direction. The bulk and bilayer heterostructure show considerable interlayer absorption at ~0.5eV, consistent with the PBE calculated quasiparticle band gap of 0.39eV. TMDs have shown promise for use in ultrathin photonic devices in part due to their relatively strong light absorption.51 To further investigate the promise of AgPbBrO and livingstonite for use in optoelectronic applications, we determine the effective masses of electrons and holes in the bulk. For AgPbBrO, the HSE-06 calculated electron effective mass in the interlayer and intralayer directions for bulk AgPbBrO is 0.57me and 0.47me, respectively, while the hole effective mass in the interlayer and intralayer directions is 0.76me and 5.75me. Remarkably, the holes in AgPbBrO have a smaller effective mass in the interlayer direction than the intralayer direction. The calculated hole effective masses in AgPbBrO are consistent with previously calculated hole effective masses in bulk PbO of 1.12me in the interlayer direction and 5.91me in the intralayer direction.52 In particular, the large in-plane hole effective mass in PbO has been attributed to the strong localization of valence band electron density on the oxygen atoms.52 Our calculations of the valence band electron density in bulk AgPbBrO (see Supporting Information) also show strong localization of electron density on the oxygen atoms.

ACS Paragon Plus Environment

18

Page 19 of 28 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

We also determine the strength of light interactions with AgPbBrO by calculating the imaginary part of the dielectric function for bulk AgPbBrO, its heterostructure bilayer, and its monolayer substituents (Figure 5b).

The magnitudes of the out of plane

component of the interlayer absorption peaks of the heterostructure are comparable in magnitude to the out of plane intralayer absorption peaks of the AgBr and PbO monolayers. This result as well as the relatively low effective masses provides promising evidence for the efficient generation of spatially separated electrons and holes in bulk AgPbBrO. The calculated electron effective mass in the interlayer and intralayer directions for bulk livingstonite is 0.51me and 0.40me, respectively, while the hole effective mass in the interlayer and intralayer directions is 1.14me and 0.64me. The band offsets in Figure 5a) suggest that the valence band holes must tunnel through the gap of the HgSb2S4 layer to move in the interlayer direction, consistent with the larger effective electron mass for propagation in the interlayer direction. No such restriction is placed on the conduction band electrons.

For comparison, the HSE-06 effective masses of bulk 2H-MoS2 in the

interlayer direction are reported to be 0.49me for electrons and 0.80me for holes.53 The effective masses for livingstonite suggest promise for the use of livingstonite in photonic applications. One of the current state-of-the-art materials used in solar cells is the hybrid inorganic-organic perovskite: methylammonium lead iodide.54 These materials boast low hole and electron effective masses (0.29 me and 0.23me ),55 but have issues associated with stability as well as lead toxicity.

ACS Paragon Plus Environment

19

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 20 of 28

While bulk MoS2 has been reported to exhibit a relatively weak exciton binding of ~84 meV,56 single layer TMDs have been reported to exhibit photogenerated holes and electrons that are tightly bound by Coulombic interactions as excitons. In these single layer systems, the exciton binding energy is typically approximately 0.5eV depending on the local environment and reflects the difference between the quasiparticle band gap and the optical absorption band gap of the material.57 It is possible to speculate that exciton binding energies could also follow a similar trend for livingstonite and AgPbBrO, with bulk binding energies suppressed by dielectric screening and single layer or bilayer structures exhibiting higher potential for strong binding. In summary, we have employed DFT to investigate the electronic and mechanical properties of newly uncovered bulk superlattice van der Waals heterostructures. These materials have been previously reported to be mined or synthesized directly into heterostructure forms, eliminating the need for challenging layer by layer fabrication procedures typically employed to create heterostructures of other layered materials. We find that the interlayer cohesive energies of AgPbBrO and HgSb4S8 are small enough to be suitable for mechanical exfoliation. The band structures of all heterostructure bilayers and their monolayer constituents were calculated. We find that the relatively weak interlayer interactions present in AgPbBrO and livingstonite allows for the electronic band structure of the heterostructure bilayers to be approximated as the superposition of the electronic band structures of their constituent monolayers. Our results suggest that AgPbBrO and livingstonite form type-II and type I heterojunctions, respectively. This is the first report of a naturally-occurring bulk van der Waals heterostructure with a

ACS Paragon Plus Environment

20

Page 21 of 28 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

calculated band gap in the visible spectrum.

We find that both AgPbBrO and

livingstonite exhibit photocarrier effective masses comparable to those reported in bulk MoS2.53 These results suggest promising potential uses for these materials in photocatalytic or photovoltaic applications. All Kohn-Sham density functional theory (KS DFT)58,59 computations were performed with the Vienna ab initio simulation package (VASP 5.4.1)60 with the projector augmented-wave (PAW) method.61 All bilayer heterostructures and 2D monolayers were initially relaxed to an ionic step energy convergence of 10-4 eV using self-consistent, periodic DFT. A vacuum space of ~20Å is included between slabs to ensure there is no interaction between periodic images. A 9x9x1 Monkhorst-Pack k-point mesh is initially used to sample the Brillouin zone with denser meshes used to check for total energy convergence.62 A plane wave basis set is implemented with a kinetic energy cut-off of 520eV for all oxygen containing structures, whereas a cut-off of 340eV is used for HgSb4S8. The exchange-correlation energy is treated with the generalized gradient approximation (GGA) PBE functional for the relaxation of structures.63 Once the relaxed structures have been obtained, the interlayer cohesive energy, Ecoh is then calculated as 𝐸𝑐𝑜ℎ = 𝐸𝑐𝑜ℎ =

((𝐸𝑚𝑜𝑛𝑜𝑙𝑎𝑦𝑒𝑟,1 + 𝐸𝑚𝑜𝑛𝑜𝑙𝑎𝑦𝑒𝑟,2) ― 𝐸𝑏𝑖𝑙𝑎𝑦𝑒𝑟) 𝐴𝑖𝑛𝑡𝑒𝑟𝑓𝑎𝑐𝑒 ((𝐸𝑚𝑜𝑛𝑜𝑙𝑎𝑦𝑒𝑟,1 + 𝐸𝑚𝑜𝑛𝑜𝑙𝑎𝑦𝑒𝑟,2) ―𝐸𝑏𝑢𝑙𝑘) 𝐴𝑖𝑛𝑡𝑒𝑟𝑓𝑎𝑐𝑒

(Equation 1) (Equation 2)

for the bilayer and bulk cohesive energies, respectively. Here, Ebilayer is the free energy of the relaxed heterostructure bilayer, Ebulk is the free energy of the relaxed bulk heterostructure, Emonolayer,1,2 are the free energies of the relaxed individual monolayer constituents of the heterostructure, and Ainterface is the area of the interface between the

ACS Paragon Plus Environment

21

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 22 of 28

two monolayers. Due to the in plane mirror symmetry present in all of the studied monolayers, the relative ordering of the layers in the heterostructure bilayer does not affect the calculated cohesive energy values. The DFT-D3 correction method is employed to provide a more accurate description of interlayer van der Waals interactions.42 Dipole corrections are applied to all slab calculations to correct for the artificial electric field created in the supercell as a result of using periodic boundary conditions.64 The band structures for all heterostructure bilayers and 2D monolayers are calculated along lines of k-points according to the symmetry of the solid.65 The PBE functional is employed for all band structure computations, while the HSE-06 functional is also used to provide a more accurate band structure of AgPbBrO and the monolayers of HgSb4S8.66 Spin-orbit coupling was included in the calculation of all PBE band structures. Band alignment calculations were also performed for HgSb4S8 and AgPbBrO.

The band

structures of the two isolated monolayers in each compound were aligned by first calculating the maximum of the Kohn-Sham potential in the middle of the vacuum region for each monolayer. Then, all of the energy bands in each monolayer are shifted such that the vacuum level of the two monolayers are equal. The effective masses of HgSb4S8 and AgPbBrO were determined from the PBE and HSE-06 calculated band structures of each bulk heterostructure. The electron and hole effective masses are calculated by fitting the band structures to Equation 3 and Equation 4, respectively. The calculated effective masses as a function of the number of k-points included can be found in the Supporting Information.

ACS Paragon Plus Environment

22

Page 23 of 28 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

∂2𝐸𝑘 ―1

(Equation 3)

𝑚 ∗ 𝑒 = ℏ2( ∂𝑘2 )

∂2𝐸𝑘 ―1

(Equation 4)

𝑚 ∗ ℎ = ―ℏ2( ∂𝑘2 )

ASSOCIATED CONTENT Supporting Information Charge density plots of the valence band and conduction band of bulk AgPbBrO. Calculated effective masses as a function of number of k-points included. Calculated interlayer distances and their comparison to experimental values for all bulk heterostructures. Calculated density of states of isolated bilayers of AgPbBrO and livingstonite and their monolayer components. Calculated Bader charges on atoms at the heterointerfaces.

AUTHOR INFORMATION Corresponding Author *E-mail: [email protected]. Phone: (650) 723-2971. Fax: (650) 725-4034. ACKNOWLEDGMENT

ACS Paragon Plus Environment

23

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 24 of 28

This work was partially supported by Army Research Office Grant W911NF-15-1-0570, by Office of Naval Research Grant N00014-15-1-2697, by NSF Grant DMR-1455050 and EECS-1436626. REFERENCES (1) (2) (3)

(4) (5) (6) (7) (8)

(9) (10) (11) (12) (13)

Geim, A. K.; Grigorieva, I. V. Van Der Waals Heterostructures. Nature 2013, 499, 419. Liu, Y.; Weiss, N. O.; Duan, X.; Cheng, H.-C.; Huang, Y.; Duan, X. Van Der Waals Heterostructures and Devices. Nat. Rev. Mater. 2016, 1, 16042. Withers, F.; Del Pozo-Zamudio, O.; Mishchenko, A.; Rooney, A. P.; Gholinia, A.; Watanabe, K.; Taniguchi, T.; Haigh, S. J.; Geim, A. K.; Tartakovskii, A. I.; et al. LightEmitting Diodes by Band-Structure Engineering in van Der Waals Heterostructures. Nat. Mater. 2015, 14, 301. Bernardi, M.; Palummo, M.; Grossman, J.C. Extraordinary Sunlight Absorption and One Nanometer Thick Photovoltaics Using Two-Dimensional Monolayer Materials. Nano Lett. 2013, 13, 3664. Qiao, J.-B.; Gong, Y.; Zuo, W.-J.; Wei, Y.-C.; Ma, D.-L.; Yang, H.; Yang, N.; Qiao, K.-Y.; Shi, J.-A.; Gu, L.; et al. One-Step Synthesis of van Der Waals Heterostructures of Graphene and Two-Dimensional Superconducting α-Mo2C. Phys. Rev. B 2017, 95, 201403. Wang, Q. H.; Kalantar-Zadeh, K.; Kis, A.; Coleman, J. N.; Strano, M. S. Electronics and Optoelectronics of Two-Dimensional Transition Metal Dichalcogenides. Nat. Nanotechnol. 2012, 7, 699. Ceballos, F.; Ju, M.-G.; Lane, S. D.; Zeng, X. C.; Zhao, H. Highly Efficient and Anomalous Charge Transfer in van Der Waals Trilayer Semiconductors. Nano Lett. 2017, 17, 1623. Fang, H.; Battaglia, C.; Carraro, C.; Nemsak, S.; Ozdol, B.; Kang, J. S.; Bechtel, H. A.; Desai, S. B.; Kronast, F.; Unal, A. A.; et al. Strong Interlayer Coupling in van Der Waals Heterostructures Built from Single-Layer Chalcogenides. Proc. Natl. Acad. Sci. 2014, 111, 6198. Rivera, P.; Schaibley, J. R.; Jones, A. M.; Ross, J. S.; Wu, S.; Aivazian, G.; Klement, P.; Seyler, K.; Clark, G.; Ghimire, N. J.; et al. Observation of Long-Lived Interlayer Excitons in Monolayer MoSe2–WSe2 Heterostructures. Nat. Commun. 2015, 6, 7242. Novoselov, K. S.; Jiang, D.; Schedin, F.; Booth, T. J.; Khotkevich, V. V.; Morozov, S. V.; Geim, A. K. Two-Dimensional Atomic Crystals. Proc. Natl. Acad. Sci. U. S. A. 2005, 102, 10451. Zhang, Z.; Ji, X.; Shi, J.; Zhou, X.; Zhang, S.; Hou, Y.; Qi, Y.; Fang, Q.; Ji, Q.; Zhang, Y.; et al. Direct Chemical Vapor Deposition Growth and Band-Gap Characterization of MoS2/h-BN van Der Waals Heterostructures on Au Foils. ACS Nano 2017, 11, 4328. Cai, Z.; Liu, B.; Zou, X.; Cheng, H.-M. Chemical Vapor Deposition Growth and Applications of Two-Dimensional Materials and Their Heterostructures. Chem. Rev. 2018, 118, 6091. Woods, J. M.; Jung, Y.; Xie, Y.; Liu, W.; Liu, Y.; Wang, H.; Cha, J. J. One-Step Synthesis of MoS2/WS2 Layered Heterostructures and Catalytic Activity of Defective Transition Metal Dichalcogenide Films. ACS Nano 2016, 10, 2004.

ACS Paragon Plus Environment

24

Page 25 of 28 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) Wofford, J. M.; Nakhaie, S.; Krause, T.; Liu, X.; Ramsteiner, M.; Hanke, M.; Riechert, H.; J. Lopes, J. M. A Hybrid MBE-Based Growth Method for Large-Area Synthesis of Stacked Hexagonal Boron Nitride/Graphene Heterostructures. Sci. Rep. 2017, 7, 43644. (15) Aretouli, K. E.; Tsipas, P.; Tsoutsou, D.; Marquez-Velasco, J.; Xenogiannopoulou, E.; Giamini, S. A.; Vassalou, E.; Kelaidis, N.; Dimoulas, A. Two-Dimensional Semiconductor HfSe2 and MoSe2/HfSe2 van Der Waals Heterostructures by Molecular Beam Epitaxy. Appl. Phys. Lett. 2015, 106, 143105. (16) Diaz, H. C.; Chaghi, R.; Ma, Y.; Batzill, M. Molecular Beam Epitaxy of the van Der Waals Heterostructure MoTe2 on MoS2: Phase, Thermal, and Chemical Stability. 2D Mater. 2015, 2, 044010. (17) Aretouli, K. E.; Tsoutsou, D.; Tsipas, P.; Marquez-Velasco, J.; Aminalragia Giamini, S.; Kelaidis, N.; Psycharis, V.; Dimoulas, A. Epitaxial 2D SnSe2/ 2D WSe2 van Der Waals Heterostructures. ACS Appl. Mater. Interfaces 2016, 8, 23222. (18) Cheng, Z.; Lin, J. Layered Organic–Inorganic Hybrid Perovskites: Structure, Optical Properties, Film Preparation, Patterning and Templating Engineering. CrystEngComm 2010, 12, 2646. (19) Chen, Q.; De Marco, N.; Yang, Y.; Song, T.-B.; Chen, C.-C.; Zhao, H.; Hong, Z.; Zhou, H.; Yang, Y. Under the Spotlight: The Organic–Inorganic Hybrid Halide Perovskite for Optoelectronic Applications. Nano Today 2015, 10, 355. (20) Molina-Mendoza, A. J.; Giovanelli, E.; Paz, W. S.; Niño, M. A.; Island, J. O.; Evangeli, C.; Aballe, L.; Foerster, M.; Zant, H. S. J. van der; Rubio-Bollinger, G.; et al. Franckeite as a Naturally Occurring van Der Waals Heterostructure. Nat. Commun. 2017, 8, 14409. (21) Velicky, M.; Toth, P. S.; Rokowski, A. M.; Rooney, A. P.; Kozikov, A.; Woods, C.R.; Mishchenko, A.; Fumagalli, L.; Yin, J.; Zolyomi, V.; Georgiou, T.; Haigh, S. J.; Novoselov, K. S.; Dryfe, R. A. W. Exfoliation of natural van der Waals heterostructures to a single unit cell thickness. Nat. Commun. 2017, 8, 14410. (22) Cheon, G.; Duerloo, K.-A. N.; Sendek, A. D.; Porter, C.; Chen, Y.; Reed, E. J. Data Mining for New Two- and One-Dimensional Weakly Bonded Solids and Lattice-Commensurate Heterostructures. Nano Lett. 2017, 17, 1915. (23) Jain, A.; Ong, S. P.; Hautier, G.; Chen, W.; Richards, W. D.; Dacek, S.; Cholia, S.; Gunter, D.; Skinner, D.; Ceder, G.; et al. Commentary: The Materials Project: A Materials Genome Approach to Accelerating Materials Innovation. APL Mater. 2013, 1, 011002. (24) Hervoches, C. H.; Lightfoot, P. A Variable-Temperature Powder Neutron Diffraction Study of Ferroelectric Bi4Ti3O12. Chem. Mater. 1999, 11, 3359. (25) Villafuerte-Castrejón, M. E.; Alvarez-Fregoso, O.; Sansores, L. E.; Sànchez-Arjona, A.; Duque, J.; Pomés, R. Crystal Structure Refinement of the Ferroelectric Ceramic Compound Sm2Bi2Ti3O12. Powder Diffr. 1998, 13, 113. (26) Dorrian, J. F.; Newnham, R. E.; Smith, D. K.; Kay, M. I. Crystal Structure of Bi4Ti3O12. Ferroelectrics 1972, 3, 17. (27) Tahara, S.; Shimada, A.; Kumada, N.; Sugahara, Y. Characterization of Bi5Nb3O15 by Refinement of Neutron Diffraction Pattern, Acid Treatment and Reaction of the AcidTreated Product with n-Alkylamines. J. Solid State Chem. 2007, 180, 2517. (28) Nalini, G.; Subbanna, G. N.; Guru Row, T. N. Studies on N=2 Aurivillius Phases: Structure of the Series Bi3−xLaxTiNbO9 (0≤x≤1). Mater. Chem. Phys. 2003, 82, 663. (29) Riebe, H.-J.; Keller, H.-L. AgPbOBr — Ein Neuer Sillén-Typ? Darstellung Und Kristallstruktur. Z. Für Anorg. Allg. Chem. 1991, 597, 151.

ACS Paragon Plus Environment

25

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 26 of 28

(30) Srikrishnan, T.; Nowacki, W. A Redetermination of the Crystal Structure of Livingstonite, HgSb4S8*. Z. Für Krist. 2010, 141, 174. (31) Levy, B. Photochemical Conversion and Storage of Solar Energy; Springer Netherlands: Dordrecht, 1991. (32) Liu, S.; Zheng, M.; Chen, R.; Wang, Z. One-Pot Synthesis of an AgBr/ZnO Hierarchical Structure with Enhanced Photocatalytic Capacity. RSC Adv. 2017, 7, 31230. (33) Niizeki, N.; Buerger, M. J. The Crystal Structure of Livingstonite, HgSb4S8. Z. Für Krist. 1957, 109, 129. (34) Aurivillius, B. Mixed Bismuth Oxides with Layer Lattices. I. The Structure Type of CaNb2Bi2O9. Ark. Kemi. 1949, 1, 463. (35) Aurivillius, B. Mixed Bismuth Oxides with Layer Lattices. II. Structure of Bi4Ti3O12. Ark. Kemi. 1949, 1, 499. (36) Aurivillius, B. Mixed oxides with layer lattices. III. Structure of BaBi4Ti4O15 Ark. Kemi. 1950, 2, 519. (37) Kendall, K. R.; Navas, C.; Thomas, J. K.; zur Loye, H.-C. Recent Developments in Oxide Ion Conductors: Aurivillius Phases. Chem. Mater. 1996, 8, 642. (38) Subbarao, E. C. A Family of Ferroelectric Bismuth Compounds. J. Phys. Chem. Solids 1962, 23, 665. (39) Islam, M. S.; Lazure, S.; Vannier, R.-N.; Nowogrocki, G.; Mairesse, G. J. Structural and computational studies of Bi2WO6 based oxygen ion conductors. Mater. Chem. 1998, 8, 655. (40) Naresh, G.; Mandal, T. K. Excellent Sun-Light-Driven Photocatalytic Activity by Aurivillius Layered Perovskites, Bi₅-XLaxTi₃FeO₁₅ (x = 1, 2). ACS Appl. Mater. Interfaces 2014, 6, 21000. (41) Green, M. A.; Ho-Baillie, A.; Snaith, H. J. The Emergence of Perovskite Solar Cells. Nat. Photonics 2014, 8, 506. (42) Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A Consistent and Accurate Ab Initio Parametrization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. J. Chem. Phys. 2010, 132, 154104. (43) Dion, M.; Rydberg, H.; Schröder, E.; Langreth, D. C.; Lundqvist, B. I. Van Der Waals Density Functional for General Geometries. Phys. Rev. Lett. 2004, 92, 246401. (44) Ashton, M.; Paul, J.; Sinnott, S. B.; Hennig, R. G. Topology-Scaling Identification of Layered Solids and Stable Exfoliated 2D Materials. Phys. Rev. Lett. 2017, 118, 106101. (45) Li, L.; Chen, Z.; Hu, Y.; Wang, X.; Zhang, T.; Chen, W.; Wang, Q. Single-Layer SingleCrystalline SnSe Nanosheets. J. Am. Chem. Soc. 2013, 135, 1213. (46) Kumar, P.; Liu, J.; Ranjan, P.; Hu, Y.; Yamijala, S. S.; Pati, S. K.; Irudayaraj, J.; Cheng, G. J. Alpha Lead Oxide (α-PbO): A New 2D Material with Visible Light Sensitivity. Small 2018, 14, 1703346. (47) Gong, C.; Zhang, H.; Wang, W.; Colombo, L.; Wallace, R. M.; Cho, K. Band Alignment of Two-Dimensional Transition Metal Dichalcogenides: Application in Tunnel Field Effect Transistors. Appl. Phys. Lett. 2013, 103, 053513. (48) Nakamura, S.; Senoh, M.; Iwasa, N.; Nagahama, S. High-Brightness InGaN Blue, Green and Yellow Light-Emitting Diodes with Quantum Well Structures. Jpn. J. Appl. Phys. 1995, 34, L797. (49) Arakawa, Y.; Yariv, A. Quantum Well Lasers–Gain, Spectra, Dynamics. IEEE J. Quantum Electron. 1986, 22, 1887.

ACS Paragon Plus Environment

26

Page 27 of 28 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

(50) Ullah, N.; Murtaza, G.; Khenata, R.; Wong, K. M.; Alahmed, Z. A. Phase Transition, Electronic and Optical Properties of Mercury Chalcogenides under Pressure. Phase Transit. 2014, 87, 571. (51) Lopez-Sanchez, O.; Lembke, D.; Kayci, M.; Radenovic, A.; Kis, A. Ultrasensitive Photodetectors Based on Monolayer MoS2. Nat. Nanotechnol. 2013, 8, 497. (52) Berashevich, J.; Semeniuk, O.; Rowlands, J. A.; Reznik, A. Anisotropy of the Carrier Effective Masses in Bulk α-PbO. EPL Europhys. Lett. 2012, 99, 47005. (53) Peelaers, H.; Van de Walle, C. G. Effects of Strain on Band Structure and Effective Masses in MoS2. Phys. Rev. B 2012, 86, 241401. (54) Momblona, C.; Malinkiewicz, O.; Roldán-Carmona, C.; Soriano, A.; Gil-Escrig, L.; Bandiello, E.; Scheepers, M.; Edri, E.; Bolink, H. J. Efficient Methylammonium Lead Iodide Perovskite Solar Cells with Active Layers from 300 to 900 Nm. APL Mater. 2014, 2, 081504. (55) Giorgi, G.; Fujisawa, J.-I.; Segawa, H.; Yamashita, K. Small Photocarrier Effective Masses Featuring Ambipolar Transport in Methylammonium Lead Iodide Perovskite: A Density Functional Analysis. J. Phys. Chem. Lett. 2013, 4, 4213. (56) Saigal, N.; Sugunakar, V.; Ghosh, S. Exciton Binding Energy in Bulk MoS2: A Reassessment. Appl. Phys. Lett. 2016, 108, 132105. (57) Wang, G.; Chernikov, A.; Glazov, M. M.; Heinz, T. F.; Marie, X.; Amand, T.; Urbaszek, B. Excitons in Atomically Thin Transition Metal Dichalcogenides. Rev. Mod. Phys. 2018, 90, 021001. (58) Hohenberg, P.; Kohn, W. Inhomogeneous Electron Gas. Phys. Rev. 1964, 136, B864. (59) Kohn, W.; Sham, L. J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133. (60) Kresse, G.; Furthmüller, J. Efficient Iterative Schemes for Ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B 1996, 54, 11169. (61) Blöchl, P. E. Projector Augmented-Wave Method. Phys. Rev. B 1994, 50, 17953. (62) Monkhorst, H. J.; Pack, J. D. Special Points for Brillouin-Zone Integrations. Phys. Rev. B 1976, 13, 5188. (63) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865. (64) Makov, G.; Payne, M. C. Periodic Boundary Conditions in Ab Initio Calculations. Phys. Rev. B 1995, 51, 4014. (65) Setyawan, W.; Curtarolo, S. High-Throughput Electronic Band Structure Calculations: Challenges and Tools. Comput. Mater. Sci. 2010, 49, 299. (66) Heyd, J.; Scuseria, G. E.; Ernzerhof, M. J. Chem. Phys. 2003, 118, 8207.

For TOC Only:

ACS Paragon Plus Environment

27

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 28 of 28

ACS Paragon Plus Environment

28