Tunable Lattice Constant and Band Gap of Single- and Few-Layer

Mar 22, 2016 - ... (STM), scanning tunneling spectroscopy (STS), and density functional theory (DFT) calculations. ... Optical Materials 2017 64, 356-...
0 downloads 0 Views 1MB Size
Subscriber access provided by UNIV OF CALIFORNIA SAN DIEGO LIBRARIES

Letter

Tunable Lattice Constant and Band Gap of Single- and Few-Layer ZnO Junseok Lee, Dan C. Sorescu, and Xingyi Deng J. Phys. Chem. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.jpclett.6b00432 • Publication Date (Web): 22 Mar 2016 Downloaded from http://pubs.acs.org on March 24, 2016

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

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

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

Tunable Lattice Constant and Band Gap of Singleand Few-Layer ZnO

Junseok Lee,*,†,‡ Dan C. Sorescu,† and Xingyi Deng†,‡ †

National Energy Technology Laboratory, U. S. Department of Energy, Pittsburgh, PA 15236,

United States ‡

AECOM, P.O. Box 618, South Park, PA 15129, United States

*Corresponding author: [email protected]

Abstract Single and few-layer ZnO(0001) (ZnO(nL), n = 1 - 4) grown on Au(111) have been characterized via scanning tunneling microscopy (STM), scanning tunneling spectroscopy (STS), and density functional theory (DFT) calculations. We find that the in-plane lattice constants of the ZnO(nL, n ≤ 3) are expanded compared to that of the bulk wurtzite ZnO(0001). The lattice constant reaches a maximum expansion of 3% in the ZnO(2L) and decreases to the bulk wurtzite ZnO value in the ZnO(4L). The band gap decreases monotonically with increasing number of ZnO layers from 4.48 eV (ZnO(1L)) to 3.42 eV (ZnO(4L)). These results suggest that a transition from a planar to the bulk-like ZnO structure occurs around the thickness of ZnO(4L). The work also demonstrates that the lattice constant and the band gap in ultrathin ZnO can be tuned by controlling the number of layers, providing a basis for further investigation of this material.

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

TOC Graphic

The advent of graphene followed by the discovery of its unique properties1, 2 has initiated the exploration of other two-dimensional (2D) materials such as transition metal dichalcogenides (TMDs).3, 4 The field of 2D materials has become one of the focal points in material science, leading to ample opportunities for creating novel nanostructures with tailored electronic and catalytic properties. In addition to TMDs, other classes of materials such as layered metal oxides can also be classified as 2D materials.5,

6

The diverse applicability of metal oxides in areas

ranging from catalysis to electronics makes them appealing candidates for investigation.7 Among metal oxides, ZnO is one of the most studied materials in its bulk/powder form due to its applications in semiconductor devices8 and its importance in methanol synthesis from syngas mixtures.9 ZnO is a wide band gap semiconductor (Eg = 3.37 eV) with a large exciton binding energy.10 Unlike TMDs, free-standing layers of ultrathin ZnO were difficult to obtain, but recent studies showed that crystalline ZnO nanosheets11 or graphene-like ZnO membranes12 can be prepared, opening up new opportunities for further research. At ambient conditions, the most stable form of ZnO is the wurtzite structure (a schematic is shown in Figure 1a). The polar wurtzite ZnO(0001)-Zn and ZnO(0001)-O surfaces are unstable, 2

ACS Paragon Plus Environment

Page 2 of 22

Page 3 of 22

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

causing reconstructions that stabilize these surfaces.13, 14 However, in the ultrathin ZnO regime, a different stabilization mechanism is at work; where depolarization takes place via the formation of the unreconstructed planar sheets (Figure 1b).15, 16 Structurally, the bonding configuration of the Zn and O atoms in the ZnO layers changes from sp3 tetrahedral in the wurtzite structure to sp2 trigonal coordination in the planar structure. For free-standing ZnO, the critical thickness up to which the planar structure is stable was predicted computationally to be ~8 layers,15, 17 while in the case of ZnO layers supported on a Ag(111) surface, the planar structure was calculated to be stable up to a bilayer.18 Experimentally, it was found that the ZnO layers on Ag(111) and Pd(111) adopt bulk-like structures at a thickness of 4 layers.16, 19 A material’s band gap is another important property which can impact its application in photovoltaics, sensing, and photocatalysis. The band gap of bulk ZnO can be modified via alloying, doping, and defect engineering.8, 20, 21 For ultrathin ZnO, it was calculated that one could tune the band gap via strain engineering in the ZnO/substrate system,17 similar to the results obtained for the monolayer TMDs.22, 23 For stacked layers of ZnO, the band gap was calculated to decrease as the number of layers increased, eventually approaching the bulk value.18, 24, 25 Experimentally, for thicker ZnO films (> 15 ML) grown on Au(111), the band gap was measured to be 2.8 eV.26,

27

However, a systematic measurement of the band gap

modification of ZnO as a function of the number of ZnO layers in the few-layer regime is still lacking in the literature. In this work, we report on the study of the structural and electronic properties of few-layer ZnO systems. (Throughout this work, we refer to ZnO(nL) as the n-layer thick ZnO structure. For example, ZnO(1L) represents the ZnO monolayer.) We have found that the in-plane lattice constants of the ZnO(nL, n ≤ 3) are all expanded with the ZnO(2L) having the largest expansion

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

(3%). We have also found that the band gap of ultrathin ZnO decreases monotonically as the thickness of ZnO increases, reaching a value close to the bulk band gap at ZnO(4L). Briefly, STM/STS experiments were performed in a low-temperature STM system (Omicron) at liquid nitrogen temperature (T = 77.4 K) using several tungsten tips. The bias voltage refers to the sample voltage. To obtain reliable physical dimensions in the STM measurements, calibration of the piezo scanner was carefully performed prior to the measurements of the lattice constants and of the Moiré pattern periodicities (Figure S1). Reactive deposition of Zn in NO2 was employed to grow the ZnO layers as shown in Figures S2 and S3. (See SI for experimental and theoretical details.)

Figure. 1. Schematics of the wurtzite (a) and the planar (b) ZnO structures. (c) STM image of a ZnO(1L) structure on Au(111) (+1.0 V, 50 pA, 100 × 100 nm2). The yellow arrows indicate the close packed direction of Au(111) substrate. Two red ellipses indicate the areas of disturbed

4

ACS Paragon Plus Environment

Page 4 of 22

Page 5 of 22

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

Moiré patterns. The inset shows the apparent height (~2 Å) along the white solid line. At the lower right side of the image, a part of a ZnO(2L) structure is also visible. (d) High resolution STM image (+5 mV, 1.5 nA, 10 × 10 nm2) of the area indicated in the square in (c). Atomically resolved ZnO lattice and the Moiré pattern are clearly visible. The unit cells of ZnO and Moiré pattern are marked in red and green, respectively. (e) A 2D Fourier transform of the STM image in (d) with outer spots for ZnO lattice and inner spots for the Moiré pattern near the center. A straight line shows the alignment of the ZnO lattice and the Moiré pattern. The lattice constant and the Moiré periodicity of the ZnO(1L) from the Fourier transform are found to be 3.29 ± 0.02 Å and 24.02 ± 0.43 Å, respectively.

Figure 1c shows an STM image of a triangular ZnO(1L) on the Au(111) substrate. Previously, small isolated ZnO monolayer structures with apparent heights of 2.0 Å have been observed on Pd(111) and Pt(111), but not on Au(111) or Ag(111).16,

19, 28, 29

The sides of the triangular

ZnO(1L) are aligned along the close-packed direction of the Au(111) substrate as shown in Figure 1c. The apparent height of the ZnO(1L) (~2 Å) agrees well with that reported in previous studies,16, 19, 28 and suggests that the thickness of the ZnO(1L) is considerably reduced compared to the single step height of the bulk ZnO(0001) single crystal (2.6 Å).30 In our previous work,29 we assigned the feature with an apparent height of 3.5 Å to ZnO(1L) but we revise it to be ZnO(2L) based on the new measurements in the current work. We note, however, that the height measured in STM represents an apparent height at a particular bias voltage. Thus, direct comparison with a real physical dimension should be made with caution. A closer inspection of Figure 1c (red ellipses) reveals that the Moiré patterns in some areas are disturbed due to the unlifted herringbone reconstruction (Figure S4).31

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

A high resolution STM image of the ZnO(1L) surface and its 2D Fourier transform image are shown in Figures 1d and 1e, respectively. The STM image shows that the surface of the ZnO(1L) is almost free of defects and/or adsorbates compared to the surfaces of ZnO(0001) or Zn(0001) bulk crystals.14, 30 From the Fourier transform analysis, the lattice constant of the ZnO(1L) is calculated to be 3.29 ± 0.02 Å which is about 1.2% larger than that of the bulk ZnO(0001) lattice constant (3.25 Å). A straight red line drawn through two spots in Figure 1e shows that the ZnO and Au(111) lattices are aligned along the same direction (see Figure S5). Figure 2a shows an STM image of a ZnO structure on Au(111) where up to four layers of ZnO are visible. The height profile along the green solid line is shown in Figure 2b. The average apparent heights of ZnO(nL) (n = 1 - 4) are 2.0 ± 0.1, 3.7 ± 0.1, 5.5 ± 0.1, and 8.0 ± 0.3 Å, respectively (Figure 2c). As the number of ZnO layers increases, the apparent height difference reaches a minimum at the ZnO(2L) (right axis in Figure 2c) suggesting that the top layer of the ZnO(2L) is the thinnest ZnO layer among the four. The apparent height difference between the ZnO(3L) and the ZnO(4L) is almost identical to the single step height of the bulk ZnO structure (dotted line in Figure 2c). This trend resembles the surface x-ray diffraction results found in the ZnO/Ag(111) system.16

6

ACS Paragon Plus Environment

Page 6 of 22

Page 7 of 22

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

Figure 2. (a) STM image (+1.0 V, 50 pA, 95 × 95 nm2) of a ZnO structure grown on Au(111) using reactive deposition method at T = 465 K. (For larger ZnO structure images, see STM images in SI.) Up to 4 layers of ZnO are visible in the image. The layer numbers are marked on each corresponding layer. (b) Height profile along the solid green line in (a). The apparent height of the average Au(111) surface plane is set to zero. The apparent height measurements have also been performed at V = -1.0 V but produced similar results. (c) The average apparent height (filled circle, left axis) and the height difference (filled square, right axis) between successive ZnO layers. The dotted line indicates the single step height of the bulk ZnO crystal. The apparent height of the ZnO(1L) itself is used to indicate the height difference between the ZnO(1L) and Au(111) substrate. (d) The lattice constant of ZnO (filled circle, left axis) and the Moiré periodicity (filled square, right axis) as a function of the number of ZnO layers. The axis for the percent lattice expansion with respect to the bulk ZnO lattice (3.25 Å) is also shown in the upper



left part. The Au lattice constant (filled diamond) is derived from ( )   for the





case of 0° rotation angle between the Au substrate and the ZnO overlayer (where dAu(111), dZnO, 7

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

and dMoire are the lattice constants of Au(111), ZnO, and Moiré pattern, respectively).32 The dotted blue line indicates the lattice constant of the wurtzite ZnO bulk and the dotted red line indicates the lattice constant of the ideal Au(111) surface. (e) A schematic of an optimized (9 × 9) ZnO(1L) overlayer on a (16 × 16) three-layer Au slab from DFT calculations. A similar Au supercell model has been used for the ZnO(nL) (n = 2 - 4) structures. The lattice constant of each ZnO(nL) layer was obtained by averaging all the nearest neighbor O-O distances in the top layer excluding the edge area. (f) Comparison of the experimental (black filled squares) and DFT calculated (red filled circles) lattice constants as a function of the number of ZnO layers. The axis on the right indicates the percentage lattice expansion with respect to the bulk ZnO lattice constant.

Figure 2d shows the average in-plane lattice constants of ZnO(nL) (n = 1 – 4) and their Moiré pattern periodicities. The lattice constants of ZnO and the Moiré periodicity of each ZnO layer were obtained by employing a 2D Fourier transform of the corresponding high resolution STM image as demonstrated in Figures 1d and 1e. Interestingly, the lattice constant as a function of the number of the ZnO layers shows a maximum at the ZnO(2L), and then decreases to almost the bulk ZnO value at the ZnO(4L). This suggests that there is an inverse relationship between the in-plane lattice constant and the thickness of the ZnO layer. A similar inverse relationship is found in Figure 2d between the measured periodicity of the Moiré pattern and the lattice constant. This finding is expected from the mathematical relationship between the periodicities of the substrate and the overlayer assuming that the lattice constant of the Au(111) substrate does not change.32 Indeed, the calculated lattice constant of Au(111) underneath each ZnO layer is close to the ideal unreconstructed lattice constant (2.88 Å) of the Au(111) surface. According to the

8

ACS Paragon Plus Environment

Page 8 of 22

Page 9 of 22

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

results in Figures 2c and 2d, the ZnO(2L) is the thinnest among the four layers and has the most expanded in-plane lattice constant compared to the bulk ZnO (~3% expansion). While a previous study suggests a fixed lattice expansion ratio regardless of the number of the ZnO layers,16 our layer-specific measurements based on STM clearly demonstrate that the lateral lattice expansion occurs at varying degrees in these ZnO layers, reflecting changes in the local energy relaxation landscape. In order to model the lateral variations in the lattice constants of the ZnO structures, a slab model containing a (9 × 9) ZnO overlayer on a larger (16 × 16) Au slab was considered. Previous DFT studies of ZnO supported on metallic substrates have employed a so called ‘coincidence’ structure to describe the Moiré patterns and to simplify the modelling process via periodic boundary conditions. For example, a 7/8 coincidence structure between ZnO and the substrate was employed.29 Once a coincidence ratio is selected,18, 29 the lateral lattice constant is forced to be fixed at a certain ratio to the substrate lattice constant due to the periodic boundary condition. Thus, the change in the lateral lattice constant measured in the experiment cannot be reproduced in such periodic models. The (9 × 9) ZnO / (16 × 16) Au slab model in Figure 2e allows relaxation of the ZnO lattice in both lateral and vertical directions. The optimizations have been performed using DFT calculations with a mixed Gaussian and plane-wave formalism as implemented in the Quickstep module of the CP2K program.33, 34 The PBE functional35 corrected to include Grimme-D3 long range dispersion corrections was used in the calculations.36 The calculated lattice constants of ZnO layers are shown in Figure 2f. It should be noted that upon relaxation of the (9 × 9) ZnO / (16 × 16) Au structure, there are some distortions near the edges of ZnO structures but the overall ZnO structure remains planar. In the case of ZnO(4L), a very small corrugation exists with a preference for O-termination. For comparison, the experimental

9

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

Page 10 of 22

results are also plotted in the same figure. The computational results predict that the maximum lattice expansion to take place in the ZnO(3L), whereas the experimental measurements have shown that the ZnO(2L) has the largest lattice constant (~3% expansion with respect to the bulk value). The origin of the discrepancy is not clear but it might arise either from the significantly smaller size of the computational model compared to the real ZnO structure which encompasses multiple Moiré regions, or from other modifications that were not considered in the current computational model such as oxygen vacancies37 and/or the presence of adsorbates such as hydroxyl species. Nevertheless, both experimental and theoretical results suggest that the lateral lattice expansion from the bulk value diminishes with the increasing number of ZnO layers after reaching a maximum. We believe this is due to a delicate balance in the structural stabilization mechanism between the depolarization by interlayer relaxation versus the formation of thermodynamically more stable bulk structures at thicker layers.

10

ACS Paragon Plus Environment

Page 11 of 22

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

Figure 3. (a) (Lower panel) Differential conductance (dI/dV) spectrum of the clean Au(111) surface. The spectrum is shown in logarithmic scale. The surface state peak at V = -0.4 eV is marked with an arrow. (Upper panel) Differential conductance (dI/dV) measurements for each ZnO layer (Lock-in amplifier parameters: Vmod = 20 mV, f = 361 Hz, Vset = 1.0 V, Iset = 50 pA). To better determine the band edges, the logarithm of dI/dV signal is taken.38 The positions of the band edges of each spectrum are indicated with dotted lines. The band gap for ZnO(4L), Eg(4L), is also indicated. The measurements have been performed on the surface of ZnO layers, away from the edges and adsorbate features. (b) Comparison of the experimental (black circle) and DFT calculated band gaps as a function of the number of ZnO layers. Calculated results were obtained at PBE-D3, PBE-D3 + U (U = 8.5 eV), and HSE06 (with 0.25 Hartree-Fock exchange) theoretical levels using the average lattice constants obtained from (9 × 9) ZnO / (16 × 16) Au slab model. The dotted line indicates the band gap of the bulk ZnO. (c) Influence of different lattice constants on the band gap. The band gaps were calculated using PBE-D3 functional for three different sets of in-plane lattice constants corresponding to the optimized bulk ZnO structure, the 7/8 ZnO/Au coincidence structure, and the (9 × 9) ZnO / (16 × 16) Au model used in (b), respectively. (d) Projected density of states (PDOS) of ZnO layers on the Au(111) substrate for 7/8 ZnO/Au coincidence structure. These results were obtained using the PBE-D3 + U (U = 8.5 eV) functional. The onset of the conduction band minimum of each ZnO layer is indicated by a vertical black line. The PDOS curves are displaced vertically to enhance the visibility.

The electronic structures of the ZnO layers have been studied using STS. The differential conductance (dI/dV) spectroscopy in the lower panel in Figure 3a shows the dI/dV spectrum of

11

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

Page 12 of 22

the clean Au(111) area, where the surface state peak is clearly visible near V = -0.4 eV.39 The upper panel in Figure 3a shows dI/dV spectra measured at the surface of each ZnO layer. The spectra were taken at positions away from the edges on relatively large ZnO areas (> 10 × 10 nm2) to avoid the effect of the edge states. We took an average of spectra for each ZnO layer, and determined the onset of the valence band maximum (VBM) and conduction band minimum (CBM) by taking the logarithm of the spectra following previous STS work.38 The band gap of the ZnO layers decreases as the number of layers increases, becoming quite similar to that of the bulk ZnO at ZnO(4L) (Figure 3b).40 The CBM decreases by as much as 1 eV as the number of ZnO layers increases while the VBM remains near V = -2 eV with respect to the Fermi level (EF, V = 0). In the dI/dV spectra of ZnO(1L), there is also a feature at V = ~1.75 eV with an unknown nature at this point. The decreasing trend of the ZnO band gap as a function of the number of layers suggests potential tunability of the band gap of ultrathin ZnO layers. Similar results were found for other TMD systems.3, 41, 42 The variation of the band structure as a function of the number of ZnO layers has also been analyzed using DFT calculations. The band gap calculations were performed using plane-wave basis sets generated using the projector augmented wave (PAW) method of Blöch as implemented in VASP code.43 Firstly, we analyzed the variation of the band gap using PBE-D3, PBE-D3+U (U = 8.5 eV), and HSE06 hybrid functionals. The individual layers were fully optimized while the lateral size of the simulation box was fixed at the average lattice dimension obtained from the calculations of the (9 × 9) ZnO / (16 × 16) Au model. For all ZnO layers, the VBM and the CBM are located at the Γ point in the band structure, indicating direct band gaps (see Figure S6). Figure 3b shows the calculated band gaps along with the experimental results. Among different methods considered, the best agreement with the experimental data was

12

ACS Paragon Plus Environment

Page 13 of 22

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

obtained by hybrid functional HSE06. However, all three functionals predict monotonic decrease of the band gaps as the number of ZnO layers increases. We also analyzed the influence of the in-plane lattice constant change on the band gap of ZnO. Figure 3c shows that there is little change in the predicted band gaps for ZnO(1L) with different lattice constants. Compared to the (9 × 9) ZnO / (16 × 16) Au model, the calculated band gaps are 1.6% and 2.9% lower for the ZnO(4L) with 7/8 coincidence and bulk-like structures, respectively. Thus, the band gap change observed in this work is mainly caused by the stacking of the ZnO layers. Lastly, we considered the effect of the Au substrate on the band gap variation. Relatively low adhesion energies of ZnO(nL, n = 1 - 4) on Au were found in the range of 0.40 – 0.50 eV / ZnO. In addition, a small charge transfer ranging from 0.01 electrons / ZnO (for ZnO(4L)) to 0.02 electrons / ZnO (for ZnO(1L)) was found in the Bader charge analysis. The majority of the transferred charge is found to reside in the first ZnO layer regardless of the thickness of ZnO layer. Overall, both the adhesion energy and Bader charge analyses indicate that the interaction between the ZnO and the Au substrate is weak. We have also analyzed the variation of the projected density of states (PDOS) as shown in Figure 3d. Similar to the trends found for the isolated ZnO layers (Figure 3b), we observe the band gap narrowing as the number of ZnO layers increases.

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

Page 14 of 22

Figure 4. (a) STM image of a ZnO structure on Au(111) (82 × 82 nm2, 1.5 V, 50 pA). Up to 4 layers of ZnO are visible. (b) Normalized dz/dV spectra of ZnO layers. The dz/dV spectrum for ZnO(4L) contains noise due to the changes in the tunneling junction. (c)-(f) Differential conductance (dI/dV) maps of the ZnO structure in (a). The base of the original topographic image in (a) is indicated by a dotted triangle in each figure. (c) ZnO(1L), V = 3.2 V. (d) ZnO(2L), V = 2.5 V. (e) ZnO(3L), V = 2.0 V. (f) ZnO(4L), V = 1.6 V.

14

ACS Paragon Plus Environment

Page 15 of 22

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 layer-specific electronic properties of ultrathin ZnO layers shown above provide a unique opportunity to distinguish among these nanostructures which would otherwise be difficult to be identified solely based on the apparent height measurements. In Figure 4a, we show the differential conductance maps of a ZnO nanostructure at bias voltages that have been identified in the dz/dV spectroscopy measurements.44, 45 In Figure 4b, the normalized dz/dV spectra are shown for each ZnO layer. The peak positions in Figure 4b correspond to the maximum change in the slope of the z(V) spectra (i.e. the energy where the density of states change is maximum). The dz/dV peak position decreases as the number of the ZnO layers increases, following the same trend as the band gap change in Figure 3a. The differential conductance (dI/dV) maps of the ZnO structure are shown in Figures 4c-4f at the peak energy positions identified in dz/dV spectra. For example, only the ZnO(1L) is revealed at V = 3.2 V in the dI/dV map in Figure 4c. Note that there is almost no intensity in the center of the ZnO structure in Figure 4c because the rate of change in the density of states is minimal at this position for ZnO(2L) – ZnO(4L) as seen by the flat spectra at 3.2 eV in Figure 4b. The ZnO(2L) and ZnO(3L) are also clearly revealed in Figures 4d and 4e, respectively. However, in the dI/dV map of the ZnO(4L) shown in Figure 4f, there is a slight contribution from the ZnO(3L). This is caused by the proximity of the energy positions of the dz/dV peaks of the ZnO(3L) and ZnO(4L) as shown in Figure 4b, where the dz/dV peak for the ZnO(3L) extends into the V = 1.6 eV range. Thus, by employing the dI/dV mapping technique, clear identification of each ZnO layer in a complex multilayer ZnO structure is possible, which could be extended to other nanostructure studies. In summary, we find that the in-plane lattice constant of the ultrathin ZnO is expanded in ZnO(nL, n ≤ 3) compared to the lattice constant of the bulk wurtzite structure, and nearly reaches

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

Page 16 of 22

the bulk value in ZnO(4L). The expansion of the lattice reaches a maximum of about 3% at ZnO(2L). The band gap of the ZnO is found to decrease monotonically from 4.48 eV for ZnO(1L) to 3.42 eV for ZnO(4L). These findings are supported by the DFT calculations. The results demonstrate that the lattice constant and the band gap of ultrathin ZnO can be tuned by controlling the number of the layers in the film. Our study provides new insight into the property changes in the few-layer regime of widely used ZnO and establishes a basis for further investigation of this material.

Acknowledgment This technical effort was performed in support of the National Energy Technology Laboratory’s ongoing research under the RES contract DE-FE0004000. This report was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor any agency thereof, nor any of their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. Supporting Information Available: Experimental, theoretical methods, and supporting figures.

16

ACS Paragon Plus Environment

Page 17 of 22

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

References (1) Novoselov, K. S.; Geim, A. K.; Morozov, S. V.; Jiang, D.; Zhang, Y.; Dubonos, S. V.; Grigorieva, I. V.; Firsov, A. A. Electric Field Effect in Atomically Thin Carbon Films. Science 2004, 306, 666-669. (2) Novoselov, K. S.; Falko, V. I.; Colombo, L.; Gellert, P. R.; Schwab, M. G.; Kim, K. A Roadmap for Graphene. Nature 2012, 490, 192-200. (3) Mak, K. F.; Lee, C.; Hone, J.; Shan, J.; Heinz, T. F. Atomically Thin MoS2: A New DirectGap Semiconductor. Phys. Rev. Lett. 2010, 105, 136805. (4) Chhowalla, M.; Shin, H. S.; Eda, G.; Li, L. J.; Loh, K. P.; Zhang, H. The Chemistry of TwoDimensional Layered Transition Metal Dichalcogenide Nanosheets. Nat. Chem. 2013, 5, 263275. (5) Zhang, H. Ultrathin Two-Dimensional Nanomaterials. ACS Nano 2015, 9, 9451-9469. (6) Geim, A. K.; Grigorieva, I. V. Van der Waals Heterostructures. Nature 2013, 499, 419-425. (7) Osada, M.; Sasaki, T. Exfoliated Oxide Nanosheets: New Solution to Nanoelectronics. J. Mater. Chem. 2009, 19, 2503-2511. (8) Janotti, A.; Van de Walle, C. G. Fundamentals of Zinc Oxide as a Semiconductor. Rep. Prog. Phys. 2009, 72, 126501. (9) Behrens, M.; Studt, F.; Kasatkin, I.; Kühl, S.; Hävecker, M.; Abild-Pedersen, F.; Zander, S.; Girgsdies, F.; Kurr, P.; Kniep, B.-L., et al. The Active Site of Methanol Synthesis over Cu/ZnO/Al2O3 Industrial Catalysts. Science 2012, 336, 893-897. (10) Wang, Z. L. Zinc Oxide Nanostructures: Growth, Properties and Applications. J. Phys. Condens. Matter 2004, 16, R829-R858.

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

Page 18 of 22

(11) Demel, J.; Pleštil, J.; Bezdička, P.; Janda, P.; Klementová, M.; Lang, K. Few-Layer ZnO Nanosheets: Preparation, Properties, and Films with Exposed {001} Facets. J. Phys. Chem. C 2011, 115, 24702-24706. (12) Quang, H. T.; Bachmatiuk, A.; Dianat, A.; Ortmann, F.; Zhao, J.; Warner, J. H.; Eckert, J.; Cunniberti, G.; Rümmeli, M. H. In Situ Observations of Free-Standing Graphene-like Mono- and Bilayer ZnO Membranes. ACS Nano 2015, 9, 11408-11413. (13) Dulub, O.; Diebold, U.; Kresse, G. Novel Stabilization Mechanism on Polar Surfaces: ZnO(0001)-Zn. Phys. Rev. Lett. 2003, 90, 016102. (14) Lauritsen, J. V.; Porsgaard, S.; Rasmussen, M. K.; Jensen, M. C. R.; Bechstein, R.; Meinander, K.; Clausen, B. S.; Helveg, S.; Wahl, R.; Kresse, G., et al. Stabilization Principles for Polar Surfaces of ZnO. ACS Nano 2011, 5, 5987-5994. (15) Freeman, C. L.; Claeyssens, F.; Allan, N. L.; Harding, J. H. Graphitic Nanofilms as Precursors to Wurtzite Films: Theory. Phys. Rev. Lett. 2006, 96, 066102. (16) Tusche, C.; Meyerheim, H. L.; Kirschner, J. Observation of Depolarized ZnO(0001) Monolayers: Formation of Unreconstructed Planar Sheets. Phys. Rev. Lett. 2007, 99, 026102. (17) Wu, D.; Lagally, M. G.; Liu, F. Stabilizing Graphitic Thin Films of Wurtzite Materials by Epitaxial Strain. Phys. Rev. Lett. 2011, 107, 236101. (18) Wang, C.; Wang, W.-H.; Lu, F.; Cheng, Y.; Ren, L.; Wang, W.; Liu, H. Pure and Oxidized Ag Substrate Influence on the Phase Transformation and Semiconducting Behaviors of Layered ZnO: A First-Principles Study. J. Phys. Chem. C 2015, 119, 4891-4897. (19) Weirum, G.; Barcaro, G.; Fortunelli, A.; Weber, F.; Schennach, R.; Surnev, S.; Netzer, F. P. Growth and Surface Structure of Zinc Oxide Layers on a Pd(111) Surface. J. Phys. Chem. C 2010, 114, 15432-15439.

18

ACS Paragon Plus Environment

Page 19 of 22

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

(20) Kumar, P.; Singh, J.; Parashar, V.; Singh, K.; Tiwari, R. S.; Srivastava, O. N.; Ramam, K.; Pandey, A. C. Investigations on Structural, Optical and Second Harmonic Generation in Solvothermally Synthesized Pure and Cr-Doped ZnO Nanoparticles. CrystEngComm 2012, 14, 1653-1658. (21) Singh, J.; Kumar, P.; Hui, K. S.; Hui, K. N.; Ramam, K.; Tiwari, R. S.; Srivastava, O. N. Synthesis, Band-Gap Tuning, Structural and Optical Investigations of Mg Doped ZnO Nanowires. CrystEngComm 2012, 14, 5898-5904. (22) He, K.; Poole, C.; Mak, K. F.; Shan, J. Experimental Demonstration of Continuous Electronic Structure Tuning via Strain in Atomically Thin MoS2. Nano Lett. 2013, 13, 29312936. (23) Hui, Y. Y.; Liu, X.; Jie, W.; Chan, N. Y.; Hao, J.; Hsu, Y.-T.; Li, L.-J.; Guo, W.; Lau, S. P. Exceptional Tunability of Band Energy in a Compressively Strained Trilayer MoS2 Sheet. ACS Nano 2013, 7, 7126-7131. (24) Kang, J.; Zhang, Y.; Wen, Y.-H.; Zheng, J.-C.; Zhu, Z.-Z. First-Principles Study on the Structural and Electronic Properties of Ultrathin ZnO Nanofilms. Phys. Lett. A 2010, 374, 10541058. (25) Topsakal, M.; Cahangirov, S.; Bekaroglu, E.; Ciraci, S. First-Principles Study of Zinc Oxide Honeycomb Structures. Phys. Rev. B 2009, 80, 235119. (26) Stavale, F.; Pascua, L.; Nilius, N.; Freund, H.-J. Morphology and Luminescence of ZnO Films Grown on a Au(111) Support. J. Phys. Chem. C 2013, 117, 10552-10557. (27) Stavale, F.; Nilius, N.; Freund, H.-J. STM Luminescence Spectroscopy of Intrinsic Defects in ZnO(0001) Thin Films. J. Phys. Chem. Lett. 2013, 4, 3972-3976.

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

Page 20 of 22

(28) Liu, B.-H.; McBriarty, M. E.; Bedzyk, M. J.; Shaikhutdinov, S.; Freund, H.-J. Structural Transformations of Zinc Oxide Layers on Pt(111). J. Phys. Chem. C 2014, 118, 28725-28729. (29) Deng, X.; Yao, K.; Sun, K.; Li, W.-X.; Lee, J.; Matranga, C. Growth of Single- and Bilayer ZnO on Au(111) and Interaction with Copper. J. Phys. Chem. C 2013, 117, 11211-11218. (30) Dulub, O.; Boatner, L. A.; Diebold, U. STM Study of the Geometric and Electronic Structure of ZnO(0 0 0 1)-Zn, (0 0 0 1)-O, (1 0 1 0), and (1 1 2 0) Surfaces. Surf. Sci. 2002, 519, 201-217. (31) Cosma, D. A.; Wallbank, J. R.; Cheianov, V.; Fal'ko, V. I. Moire Pattern as a Magnifying Glass for Strain and Dislocations in van der Waals Heterostructures. Faraday Discuss. 2014, 173, 137-143. (32) Hermann, K. Periodic Overlayers and Moire Patterns: Theoretical Studies of Geometric Properties. J. Phys. Condens. Matter 2012, 24, 314210. (33) VandeVondele, J.; Krack, M.; Mohamed, F.; Parrinello, M.; Chassaing, T.; Hutter, J. QUICKSTEP: Fast and Accurate Density Functional Calculations Using a Mixed Gaussian and Plane Waves Approach. Comput. Phys. Commun. 2005, 167, 103-128. (34) Hutter, J.; Iannuzzi, M.; Schiffmann, F.; VandeVondele, J. CP2K: Atomistic Simulations of Condensed Matter Systems. WIREs Comput. Mol. Sci. 2014, 4, 15-25. (35) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865-3868. (36) 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.

20

ACS Paragon Plus Environment

Page 21 of 22

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

(37) Kan, E. J.; Deng, K. M.; Wu, F. Stability of Graphitic-Like Zinc Oxide Layers Under Carriers Doping: A First-Principles Study. Nanoscale 2013, 5, 12111-12114. (38) Ugeda, M. M.; Bradley, A. J.; Shi, S.-F.; da Jornada, F. H.; Zhang, Y.; Qiu, D. Y.; Ruan, W.; Mo, S.-K.; Hussain, Z.; Shen, Z.-X., et al. Giant Bandgap Renormalization and Excitonic Effects in a Monolayer Transition Metal Dichalcogenide Semiconductor. Nat. Mater. 2014, 13, 1091-1095. (39) Davis, L. C.; Everson, M. P.; Jaklevic, R. C.; Shen, W. Theory of the Local Density of Surface States on a Metal: Comparison with Scanning Tunneling Spectroscopy of a Au(111) Surface. Phys. Rev. B 1991, 43, 3821-3830. (40) Mang, A.; Reimann, K.; Rubenacke, S. Band-Gaps, Crystal-Field Splitting, Spin-OrbitCoupling, and Exciton Binding-Energies in ZnO under Hydrostatic-Pressure. Solid State Commun. 1995, 94, 251-254. (41) Bradley, A. J.; Ugeda, M. M.; da Jornada, F. H.; Qiu, D. Y.; Ruan, W.; Zhang, Y.; Wickenburg, S.; Riss, A.; Lu, J.; Mo, S. K., et al. Probing the Role of Interlayer Coupling and Coulomb Interactions on Electronic Structure in Few-Layer MoSe2 Nanostructures. Nano Lett. 2015, 15, 2594-2599. (42) Kim, H. C.; Kim, H.; Lee, J. U.; Lee, H. B.; Choi, D. H.; Lee, J. H.; Lee, W. H.; Jhang, S. H.; Park, B. H.; Cheong, H., et al. Engineering Optical and Electronic Properties of WS2 by Varying the Number of Layers. Acs Nano 2015, 9, 6854-6860. (43) Blöchl, P. E. Projector Augmented-Wave Method. Phys. Rev. B 1994, 50, 17953-17979. (44) Binnig, G.; Frank, K. H.; Fuchs, H.; Garcia, N.; Reihl, B.; Rohrer, H.; Salvan, F.; Williams, A. R. Tunneling Spectroscopy and Inverse Photoemission: Image and Field States. Phys. Rev. Lett. 1985, 55, 991-994.

21

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

(45) Dougherty, D. B.; Maksymovych, P.; Lee, J.; Yates, J. T. Local Spectroscopy of ImagePotential-Derived States: From Single Molecules to Monolayers of Benzene on Cu(111). Phys. Rev. Lett. 2006, 97, 236806.

22

ACS Paragon Plus Environment

Page 22 of 22