Large-scale quantum-emitter arrays in atomically thin semiconductors

Abstract

Quantum light emitters have been observed in atomically thin layers of transition metal dichalcogenides. However, they are found at random locations within the host material and usually in low densities, hindering experiments aiming to investigate this new class of emitters. Here, we create deterministic arrays of hundreds of quantum emitters in tungsten diselenide and tungsten disulphide monolayers, emitting across a range of wavelengths in the visible spectrum (610-680 nm and 740-820 nm), with a greater spectral stability than their randomly occurring counterparts. This is achieved by depositing monolayers onto silica substrates nanopatterned with arrays of 150-nm-diameter pillars ranging from 60 to 190 nm in height. The nanopillars create localized deformations in the material resulting in the quantum confinement of excitons. Our method may enable the placement of emitters in photonic structures such as optical waveguides in a scalable way, where precise and accurate positioning is paramount.

ARTICLE

Received 24 Jan 2017 | Accepted 24 Feb 2017 | Published 22 May 2017

Carmen Palacios-Berraquero1,*, Dhiren M. Kara1,*, Alejandro R.-P. Montblanch1, Matteo Barbone1,2, Pawel Latawiec3, Duhee Yoon2, Anna K. Ott2, Marko Loncar3, Andrea C. Ferrari2 & Mete Atatre1

Quantum light emitters have been observed in atomically thin layers of transition metal dichalcogenides. However, they are found at random locations within the host material and usually in low densities, hindering experiments aiming to investigate this new class of emitters. Here, we create deterministic arrays of hundreds of quantum emitters in tungsten diselenide and tungsten disulphide monolayers, emitting across a range of wavelengths in the visible spectrum (610680 nm and 740820 nm), with a greater spectral stability than their randomly occurring counterparts. This is achieved by depositing monolayers onto silica substrates nanopatterned with arrays of 150-nm-diameter pillars ranging from 60 to 190 nm in height. The nanopillars create localized deformations in the material resulting in the quantum connement of excitons. Our method may enable the placement of emitters in photonic structures such as optical waveguides in a scalable way, where precise and accurate positioning is paramount.

DOI: 10.1038/ncomms15093 OPEN

Large-scale quantum-emitter arrays in atomically thin semiconductors

1 Cavendish Laboratory, University of Cambridge, JJ Thomson Avenue, Cambridge CB3 0HE, UK. 2 Cambridge Graphene Centre, University of Cambridge, Cambridge CB3 0FA, UK. 3 John A. Paulson School of Engineering and Applied Science, Harvard University, 29 Oxford Street, Cambridge, Massachusetts 02138, USA. * These authors contributed equally to this work. Correspondence and requests for materials should be addressed toA.C.F. (email: mailto:[email protected]

Web End [email protected] ) or to M.A. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 8:15093 | DOI: 10.1038/ncomms15093 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15093

Transition metal dichalcogenides (TMDs) are optically active, semiconducting layered materials (LMs) which can be exfoliated down to monolayers. These are of

particular interest, as they exhibit properties such as an optically accessible valley degree of freedom that is locked to the exciton spin13, strong two-dimensional connement, favouring bound excitonic states4,5 and the opportunity to investigate many-body physics6. Intrinsic properties of LMsatomically precise interfaces, lack of dangling bonds, exibility and the possibility of stacking different LMs into functional heterostructures7,8 make them not only interesting for fundamental physics but also suitable for technological applications9,10. For these reasons, the identication of quantum emitters (QEs) in LMs1116 has generated much excitement in the eld of two-dimensional nanophotonics10,1719 and quantum technologies9,20.

Single photon emission has been seen from QEs in tungsten diselenide (WSe2) with both above bandgap1116 and resonant optical excitation21. In addition, QEs in WSe2 and tungsten disulphide (WS2) have been implemented in heterostructures to achieve electrically driven single-photon emission22. However, the origin of QEs in TMDs is still unclear, and has been assigned to both defects1116 and strain gradients21,23,24. Experiments on these QEs have, until now, been reliant on their rare and random occurrence. Deterministic creation of precisely positioned LM QEs in large numbers is important for accelerating the study of these emitters, as well as opening up the prospect for scalability and on-chip applications.

Here, we report a method to create arrays of single-photon emitting QEs in WSe2 and quantum-like emitters in WS2 using a nanopatterned silica substrate. We obtain structures with QE numbers typically in the range of hundreds. The quality of these deterministic QEs surpasses that of their randomly appearing counterparts, with spectral wanderings of B0.1 meV an order of magnitude lower than previous reports1115. Our technique is a crucial rst step towards solving the scalability challenge for LM-based quantum photonic devices.

ResultsNanopatterned substrate preparation and characterization. To create large-scale QE arrays in LMs, we place the active material on patterned structures fabricated on the substrate in order to create spatially localized physical disturbances to the otherwise at LM akes. To this end, we rst pattern arrays of nanopillars of different heights, ranging from 60 to 190 nm, on silica substrates using electron beam lithography. Figure 1a shows a scanning electron microscope image of one such substrate of 130 nm nanopillar height. We place layers of WSe2 and WS2 on the nanopillars as follows. Bulk WSe2 and WS2 crystals are characterized before exfoliation as described in ref. 22. These are then exfoliated on a polydimethylsiloxane layer by micromechanical cleavage7,25. Single-layer (1L) samples are identied rst by optical contrast26, and the selected 1L-WSe2 and 1L-WS2 akes are then placed onto the patterned nanopillar substrate via an all-dry viscoelastic transfer technique due to their higher adhesion to SiO2 (refs 7,19,27), as schematically shown in

Fig. 1b. After exfoliation and transfer, the 1L-WSe2 and 1L-WS2 akes are characterized by Raman spectroscopy28,29, photoluminescence (PL)30 and atomic force microscopy (AFM), conrming the transfer and that the process does not damage the samples (see Supplementary Fig. 1 and Supplementary Note 1 for the corresponding spectra and discussion). Figure 1c is an AFM scan of a 1L-WSe2 ake over a single nanopillar. The bottom panel of Fig. 1c plots the height prole of the 1L-WSe2 ake taken along the dashed pink line. This reveals how the ake (solid pink line) tents over the nanopillar. The blue-shaded area corresponds

to the measured prole of a bare nanopillar. Figure 1d is a dark-eld optical microscopy image of part of a 43,000 mm2 1L-WSe2 ake on a substrate patterned with a 4-mm-spaced nanopillar array with nominal height of 130 nm. The regularly spaced bright spots correspond to nanopillar sites. We see locations providing brighter scattering (two examples are encircled in pink) and others showing fainter intensity (two examples are encircled in blue). By correlating with AFM measurements we nd that the former correspond to locations where the 1L-WSe2 tents over the nanopillars and the latter correspond to locations where the ake is pierced by the nanopillars (see Supplementary Fig. 2). On average, we nd that two-third of the sites are not pierced during the deposition step.

Quantum light from WSe2-based deterministic QEs. Figure 2a is an integrated raster scan map of PL emission at B10 K of six adjacent non-pierced nanopillar sites in the region enclosed by the green dashed line in Fig. 1d. The most prominent feature is the B 10 increase in intensity at the location of every nano-

pillar. Figure 2b reveals the source of this emission intensity enhancement: spectra taken at each nanopillar location display bright sub-nanometre linewidth emission peaks. Figure 2c demonstrates the single-photon nature of this emission via photon-correlation measurements taken (from left to right) at the rst, third and fourth nanopillar locations. Ten nanometre bandpass lters, indicated by the pink, green and blue highlighted areas in the panels of Fig. 2b, select the spectral windows for the photon-correlation measurements. We obtain g(2)(0) values of0.08680.0645, 0.1700.021 and 0.1820.028, respectively, uncorrected for background emission or detector response. While these surpass those in early reports1115, we expect the quality of the single-photon emission from the QEs to improve under resonant excitation31. Out of the 53 unpierced nanopillar sites in this substrate we nd sub-nm emission peaks in 51 of them, giving B96% yield in QE generation. Their emission wavelength ranges between 730 and 820 nm (see Supplementary Fig. 3 for statistics), equivalent to a redshift distributed between 50 and 280 meV from the unbound exciton emission energy at B1.755 eV (ref. 32), as observed for the naturally occurring

QEs in WSe2 (refs 1115). The ne-structure splitting (200730 meV) and the emission linewidths as narrow as B180 meV (B0.08 nm) are also consistent with previous reports1115 (Supplementary Fig. 3) advocating that these deterministically created QEs are of the same nature as the randomly appearing ones.

Effect of nanopillar height on WSe2-based QEs. To study the effect of nanopillar height, we carry out similar optical measurements of 1L-WSe2 akes deposited on nanopillars of height

B60 and B190 nm. The spectra taken at the 60-nm nanopillars have multiple peaks of B1 nm linewidth on average (see Supplementary Fig. 3 for example spectra). In contrast, Fig. 2d is a representative spectrum taken from the 190 nm nanopillars, displaying a better isolated, single sub-nm emission peak. The inset reveals a 722 meV ne-structure splitting for this

QE. We do not see clear nanopillar height dependence in the emission wavelength and ne-structure splitting (see Supplementary Fig. 3 for statistics). Although we often observe multiple narrow emission lines from nanopillar locations, increasing the nanopillar height reduces the spread in the number of peaks arising at each. We verify this trend in Fig. 2e, a histogram of the probability that a given number of sub-nm emission peaks appear per nanopillar, for the different nanopillar heights (60, 130 and 190 nm in white, blue and purple, respectively). The likelihood of creating a single QE grows as nanopillar height is

2 NATURE COMMUNICATIONS | 8:15093 | DOI: 10.1038/ncomms15093 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15093 ARTICLE

a

d

Non-pierced Pierced

b

1.

2.

PDMS

Si/SiO2

LM

c

Figure 1 | Fabrication and characterization of scalable quantum connement arrays. (a) SEM image of nanopillar substrate, fabricated by electron beam lithography. Black scale bar, 2 mm. (b) Illustration of the fabrication method: (1) mechanical exfoliation of LM on PDMS and all-dry viscoelastic deposition on patterned substrate; and (2) deposited LM on patterned substrate. (c, top) An AFM scan of 1L-WSe2 on a nanopillar. (bottom) The AFM height prole of a bare nanopillar (blue-shaded region) and of the ake deposited over it (pink line), measured along the dashed pink line cut in the top panel. Colour-scale bar represents height in nm and white scale bar 1 mm. (d) Dark eld optical microscopy image (real colour) of 1L-WSe2 on nanopillar substrate (130 nm high, 4 mm separation). The full image corresponds to a 170 mm by 210 mm area. The green box highlights six adjacent nanopillars within the 1L-WSe2 region, measured in Fig. 2. The blue circles indicate two pierced nanopillars, and the pink circles indicate two non-pierced nanopillars. PDMS, polydimethylsiloxane; SEM, scanning electron microscope.

increased. For the 190 nm nanopillars, 50% of all nanopillar sites host a single QE with one emission peak, as indicated by the purple bars. Spectral wandering of the peaks as a function of time also displays a strong dependence on the nanopillar height. To quantify this dependence, we record the maximum range of emission wavelength wandering per QE over tens of seconds. The solid black circles in Fig. 2f correspond to the mean of these values for each group of QEs pertaining to each nanopillar height, for 17 different QEs in total, with the error bars displaying the standard deviation of these distributions. We observe a reduction from a few meV for 60 nm nanopillars to below 0.25 meV (average) for the tallest 190 nm nanopillars (see Supplementary Fig. 4), reaching as low as 0.1 meV. To the best of our knowledge, this is the lowest spectral wandering seen to date in LM QEs1115. Hence, these deterministic QEs are comparable, and even superior, in spectral stability to their randomly appearing counterparts. The dependence of certain QE characteristics on nanopillar height, along with shifts in the delocalized neutral exciton peak (X0) at room temperature30 at the nanopillar locations (see Supplementary Fig. 5), suggest that a localized strain gradient induced by the nanopillars might be playing an

active role in producing QEs, as well as determining their specic optical properties21,23,24.

Deterministic QE creation in WS2 monolayers. The method we present for QE creation is not restricted to a specic LM. We predict a similar effect on different LMs and test this by using 1L-WS2. Despite previous efforts to measure QEs in 1L-WS2, there has only been one previous report of single-photon emission in this material22. Figure 3a shows an integrated PL intensity raster scan map taken at B10 K of a 1L-WS2 on a substrate with 170-nm-high nanopillars square array spaced by 3 mm. The inset shows a true-colour dark-eld optical microscopy image of the same ake, where the red areas (due to uorescence) are 1L-WS2.

Once again, the brighter spots correspond to the unpierced nanopillar locations, as veried by AFM measurements, and show overlap with the bright uorescence spots in the PL intensity image where, similar to WSe2, intensity is increased (here by a factor B4) at every one of the 22 non-pierced nanopillar sites in the ake. Panel 1 of Fig. 3b shows the typical 1L-WS2 emission spectrum at B10 K (ref. 22), measured from a at region of the same ake away from the nanopillars. The X0 and X unbound

NATURE COMMUNICATIONS | 8:15093 | DOI: 10.1038/ncomms15093 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15093

a

WSe2

x x x x x x

0

160 k

1 m

b

PL intensity

750 800 700 750 800 700 750 800 700 750 800 700 750 800 700 750 800Wavelength (nm)

Wavelength (nm)

700

c d

1.0

0.0

(2) ( )

PL intensity

g 0.5

20 0 20 20 0 20 20 0 20

700 750 800

Time delay (ns)

e f

4

3

2

1

0 60 130 190

50

%

25

190 nm 130 nm 60 nm

0 0 1 2 3 4 5 6

wandering (meV)

Number of sub-nm emission peaks per nanopillar

Spectral

Nanopillar height (nm)

Figure 2 | Creation of quantum emitter arrays in 1L-WSe2. (a) Integrated PL intensity raster scan of the region enclosed by the green rectangle in Fig. 1d, taken under 200 nW mm 2, 532 nm laser excitation at 10 K. Green crosses mark the position of the six nanopillars beneath the 1L-WSe2. Colour-scale bar maximum, 160 kcounts s 1. (b) PL spectra taken at each of the corresponding green crosses in a, from left to right respectively, showing the presence of narrow lines at each nanopillar location. (c) Photon correlation measurements corresponding to the ltered spectral regions (10 nm wide) enclosed by the blue, green and pink rectangles, in b, with g(2)(0) 0.0870.065, 0.170.02 and 0.180.03, and rise times of 8.810.80 ns, 6.150.36 ns and

3.080.41 ns, respectively. (d) Spectrum taken from a 1L-WSe2 on a 190 nm nanopillar, showing lower background and a single sub-nm emission peak. Higher-resolution spectrum in the inset reveals the ne-structure splitting of this QE peak. An asymmetry can be seen in the spectrum, which has been previously attributed to a phonon sideband in naturally occurring QEs31. (e) Probability distribution (in %) of the number of emission lines per nanopillar for samples using different nanopillar heights (60, 130 and 190 nm in white, blue and purple, respectively). A trend of higher probability of single QE emission peaks per nanopillar location with increasing height is evident, reaching 50% for 190 nm nanopillars. (f) Increasing nanopillar height also leads to a reduction of spectral wandering. Solid black circles represent the mean value of spectral wandering of several QEs for a given nanopillar height, while the error bars represent the standard deviation of each distribution, both extracted from time-resolved high-resolution spectral measurements (Supplementary Fig. 4). A total number of seven samples was used to collect the statistics necessary for Fig. 2e,f.

excitons are labelled in the gure, while the broad red-shifted emission band arises from weakly localized or defect-related excitons in the 1L-WS2 at low temperatures22, and is present in this material regardless of location. Panels 2 and 3 of Fig. 3b show representative PL emission spectra taken at nanopillars of heights B170 andB190 nm, respectively, where once again sub-nm spectral features arise. We note that we observe ne-structure splitting for WS2 in these QEs, which range from 300 to 810 meV (Supplementary Fig. 6), as represented in the panel insets corresponding to the spectral regions highlighted in red. We also measure the spectrum of several WS2 QEs as a function of time (see Supplementary Fig. 6) and nd all spectral wandering

values below 0.5 meV over 12 min. Figure 3c shows statistics on QE emission wavelength collected for over B80 QEs for 1L-WS2 on 170-nm (white bars) and 190-nm (red bars) nanopillars. The wavelength distribution of the sub-nm emission lines, typically in the 610680 nm region (53300 meV redshift from X0)22, is as narrow as B20 nm for the 190-nm nanopillars. Most nanopillar sites on WS2 show multiple sub-nm lines, suggesting the creation of several QEs at each site for these nanopillar heights. Figure 3d plots a histogram of the number of sub-nm peaks appearing at each nanopillar for both nanopillar heights. The trend is similar to that seen in WSe2, where higher nanopillars lead to a narrower spread in the number of peaks towards a higher likelihood of

4 NATURE COMMUNICATIONS | 8:15093 | DOI: 10.1038/ncomms15093 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15093 ARTICLE

a

b

1.

2.

3.

WS2

2 k

18 k

6 m

PL intensity

600 625 650 675 700

Wavelength (nm)

c

d

50

Count

20

15

10

5

%

25

0 600 625 650 675 700

0 0 1 2 3 4 5 6 7 8

Wavelength of sub-nm emission peaks (nm) Number of sub-nm emission peaks per nanopillar

Figure 3 | Creation of quantum-emitter arrays in 1L-WS2. (a) Integrated PL intensity raster scan of a 1L-WS2 ake deposited on top of a 3 mm spaced, 170-nm-high nanopillar array, taken at 300 nW mm 2, 532 nm laser excitation at 10 K. Colour-scale bar maximum is 18 kcounts s 1. Inset: true-colour DFM image of the same area. The red region corresponds to the WS2 monolayer. (b) PL spectra of 1L-WS2 at 10 K. Panel 1 shows a spectrum taken from a at region away from nanopillars. Red arrows indicate unbound monolayer neutral (X0) and charged (X ) excitons. Panels 2 and 3 show representative spectra of WS2 on 170 and 190 nm nanopillars, respectively. Insets are high-resolution PL spectra of the red-highlighted spectral regions, showing the ne-structure splitting of the peaks. (c) Distribution of the emission wavelengths measured for 1L-WS2 QEs on 170 (black and white) and 190 nm (red)

nanopillars. (d) Distribution of the number of narrow emission lines observed per nanopillar for 1L-WS2 QEs on 170 (black and white) and 190 nm (red) nanopillars. DFM, dark-eld optical microscopy.

creating a single QE at each nanopillar site. We note that we obtain a 95% yield of QE creation in 1L-WS2 on non-pierced nanopillars. Further, B75% of these display two or less sub-nm emission peaks. In contrast, the 60- and the 130-nm-high nanopillars do not result in any QE occurrence (see Supplementary Fig. 7 for examples of these measurements).

DiscussionWe presented a simple method for the deterministic creation of scalable arrays of quantum-light emitters embedded in LMs emitting at different regions of the optical spectrum33. The reliability of the technique will accelerate experimental studies of QEs in TMDs, which at present rely on their rather rare and random occurrence1115. In the immediate future, a detailed study is necessary in order to achieve a better understanding of the specic role of nanopillar height and geometry in dening the characteristics of the quantum emission. We expect tunability of the optical emission by varying the shapes of the underlying nanostructures. In this respect, interesting possibilities to realise dynamical circuits using micro-electro-mechanical systems and piezoelectric tuning exist. Heterostructures may enable new routes towards tunnel-coupled quantum devices and the formation of QE molecules. Several approaches are being investigated for the production of wafer-scale samples34,35, which could lead to rapid optimization. While our approach is already compatible with standard silicon processing techniques, it

is nevertheless not restricted to the specic properties of the substrate. In fact, even nanodiamonds of the appropriate dimensions, drop-cast onto silica substrates, are able to create QEs in 1L-WSe2 (see Supplementary Fig. 8 and Supplementary

Note 2). The choice of substrate material will be particularly important, for example, when considering inhomogeneous line-broadening due to charge noise. Further, the exibility in the choice of substrate, in turn, provides an opportunity to create hybrid quantum devices where LM QEs can be coupled to quantum systems in other materials such as spins in diamond and silicon carbide.

Methods

Substrate preparation. The silica nanopillar substrate is fabricated with a high-resolution direct-write lithographic process via spin-on-glass polymer hydrogen silsesquioxane (HSQ)36. First, a wafer with 2 mm thermal oxide is cleaved and then cleaned. HSQ resist (FOx-16, Dow-Corning) is diluted with methyl isobutyl ketone (MIBK) in different ratios and spun onto the substrate, giving variable thickness depending on the dilution. After baking at 90 C for 5 min, the substrate is exposed in an electron beam lithography tool (Elionix F-125) and then developed in a 25% solution of tetramethyl ammonium hydroxide (TMAH) developer and rinsed in methanol. To convert the dened structures into pure SiO2, we apply rapid thermal annealing at 1,000 C in an oxygen atmosphere37, resulting in arrays of sharply-dened sub-100 nm silica nanopillars.

Optical measurements. Room temperature Raman and PL measurements are carried out using a Horiba LabRam HR Evolution microspectrometer equipped with a 100 objective (numerical aperture 0.9) and a spot size o1 mm. The pixel-

NATURE COMMUNICATIONS | 8:15093 | DOI: 10.1038/ncomms15093 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15093

to-pixel spectral resolution for the Raman measurements is B0.5 cm 1. Bragg gratings (BraggGrate) are used to detect the ultralow frequency Raman peaks. The power is kept below 50 mW to prevent heating effects. The excitation wavelength used is 514.5 nm for WSe2 and 457 nm for WS2.

A variable-temperature helium ow cryostat (Oxford Instruments Microstat HiRes2) is used to perform low-temperature PL measurements with a home-built confocal microscope mounted on a three-axis stage (Physik Instrumente M-405DG) with a 5-cm travel range, 200-nm resolution for coarse alignment and a piezo scanning mirror (Physik Instrumente S-334) for high-resolution raster scans. PL is collected using a 1.7-mm working distance objective with a numerical aperture of 0.7 (Nikon S Plan Fluor 60) and detected on a bre-coupled single-

photon-counting module (PerkinElmer: SPCM-AQRH). Photon correlations from a Hanbury Brown and Twiss interferometer are recorded with a time-to-digital converter (quTAU). A double grating spectrometer (Princeton Instruments) is used for acquiring spectra. For PL measurements, the excitation laser (532 nm, Laser Quantum) is suppressed with a long pass lter (550 nm Thorlabs FEL0550).

Data availability. The data that supports the ndings of this study are available from the corresponding author upon request.

References

1. Cao, T. et al. Valley-selective circular dichroism of monolayer molybdenum disulphide. Nat. Commun. 3, 887 (2012).

2. Zeng, H., Dai, J., Yao, W., Xiao, D. & Cui, X. Valley polarization in MoS2 monolayers by optical pumping. Nat. Nanotechnol. 7, 490493 (2012).

3. Mak, K. F., He, K., Shan, J. & Heinz, T. F. Control of valley polarization in monolayer MoS2 by optical helicity. Nat. Nanotechnol. 7, 494498 (2012).

4. Mak, K. F. et al. Tightly bound trions in monolayer MoS2. Nat. Mater. 12, 207211 (2013).

5. He, K. et al. Tightly bound excitons in monolayer WSe2. Phys. Rev. Lett. 113, 26803 (2014).

6. Sidler, M. et al. Fermi polaron-polaritons in charge-tunable atomically thin semiconductors. Nat. Phys. 13, 255261 (2016).

7. Bonaccorso, F. et al. Production and processing of graphene and 2d crystals. Mater. Today 15, 564589 (2012).

8. Akinwande, D., Petrone, N. & Hone, J. Two-dimensional exible nanoelectronics. Nat. Commun. 5, 5678 (2014).

9. Majumdar, A. et al. Hybrid 2D material nanophotonics: a scalable platform for low-power nonlinear and quantum optics. ACS Photonics 2, 11601166 (2015).

10. Xia, F., Wang, H., Xiao, D., Dubey, M. & Ramasubramaniam, A. Two-dimensional material nanophotonics. Nat. Photonics 8, 899907 (2014).

11. Tonndorf, P. et al. Single-photon emission from localized excitons in an atomically thin semiconductor. Optica 2, 347 (2015).

12. Srivastava, A. et al. Optically active quantum dots in monolayer WSe2.

Nat. Nanotechnol. 10, 491496 (2015).13. He, Y.-M. et al. Single quantum emitters in monolayer semiconductors. Nat. Nanotechnol. 10, 497502 (2015).

14. Koperski, M. et al. Single photon emitters in exfoliated WSe2 structures. Nat. Nanotechnol. 10, 503506 (2015).

15. Chakraborty, C., Kinnischtzke, L., Goodfellow, K. M., Beams, R. & Vamivakas, A. N. Voltage-controlled quantum light from an atomically thin semiconductor. Nat. Nanotechnol. 10, 507511 (2015).

16. Tran, T. T., Bray, K., Ford, M. J., Toth, M. & Aharonovich, I. Quantum emission from hexagonal boron nitride monolayers. Nat. Nanotechnol. 11, 3741 (2015).

17. Bonaccorso, F., Sun, Z., Hasan, T. & Ferrari, A. C. Graphene photonics and optoelectronics. Nat. Photonics 4, 611622 (2010).

18. Koppens, F. H. L. et al. Photodetectors based on graphene, other two-dimensional materials and hybrid systems. Nat. Nanotechnol. 9, 780793 (2014).

19. Ferrari, A. C. Science and technology roadmap for graphene, related two-dimensional crystals, and hybrid systems. Nanoscale 7, 45984810 (2015).

20. Aharonovich, I., Englund, D. & Toth, M. Solid-state single-photon emitters. Nat. Photonics 10, 631641 (2016).

21. Branny, A. et al. Discrete quantum dot like emitters in monolayer MoSe2:

spatial mapping, magneto-optics, and charge tuning. Appl. Phys. Lett. 108, 142101 (2016).22. Palacios-Berraquero, C. et al. Atomically thin quantum light-emitting diodes. Nat. Commun. 7, 12978 (2016).

23. Kumar, S., Kaczmarczyk, A. & Gerardot, B. D. Strain-induced spatial and spectral isolation of quantum emitters in mono- and bilayer WSe2. Nano Lett.

15, 75677573 (2015).24. Kern, J. et al. Nanoscale positioning of single-photon emitters in atomically thin WSe2. Adv. Mater. 28, 71017105 (2016).

25. Novoselov, K. S. et al. Two-dimensional atomic crystals. Proc. Natl Acad. Sci. USA 102, 1045110453 (2005).

26. Casiraghi, C. et al. Rayleigh imaging of graphene and graphene layers. Nano Lett. 7, 27112717 (2007).

27. Castellanos-Gomez, A. et al. Deterministic transfer of two-dimensional materials by all-dry viscoelastic stamping. 2D Mater. 1, 11002 (2014).

28. Terrones, H. et al. New rst order Raman-active modes in few layered transition metal dichalcogenides. Sci. Rep. 4, 4215 (2014).

29. Zhao, W. et al. Lattice dynamics in mono- and few-layer sheets of WS2 and WSe2. Nanoscale 5, 96779683 (2013).

30. Zhou, B. et al. Evolution of electronic structure in atomically thin sheets of WS2 and WSe2. ACS Nano 7, 791797 (2013).

31. Kumar, S. et al. Resonant laser spectroscopy of localized excitons in monolayer WSe2. Optica 3, 882 (2016).

32. Wang, G. et al. Valley dynamics probed through charged and neutral exciton emission in monolayer WSe2. Phys. Rev. B 90, 075413 (2014).

33. Jiang, M., Kurvits, J. A., Lu, Y., Nurmikko, A. V. & Zia, R. Reusable inorganic templates for electrostatic self-assembly of individual quantum dots, nanodiamonds, and lanthanide-doped nanoparticles. Nano Lett. 15, 50105016 (2015).

34. Kang, K. et al. High-mobility three-atom-thick semiconducting lms with wafer-scale homogeneity. Nature 520, 656660 (2015).

35. Ling, X. et al. Parallel stitching of 2D Materials. Adv. Mater. 28, 23222329 (2016).

36. Grigorescu, A. E. et al. Resists for sub-20-nm electron beam lithography with a focus on HSQ: state of the art. Nanotechnology 20, 292001 (2009).

37. Holzwarth, C. W., Barwicz, T. & Smith, H. I. Optimization of hydrogen silsesquioxane for photonic applications. J. Vac. Sci. Technol. B Microelectron. Nanometer Struct. Process. Meas. Phenom. 25, 2658 (2007).

Acknowledgements

We acknowledge nancial support from the Marie Skodowska-Curie Actions Spin-NANO, Grant No. 676108, EU Graphene Flagship, ERC Grants Hetero2D and PHOENICS, EPSRC Grants EP/K01711X/1, EP/K017144/1, EP/N010345/1, EP/M507799/1, EP/L016087/1, Quantum Technology Hub NQIT EP/M013243/1, the EPSRC Cambridge NanoDTC, Graphene Technology CDT, EP/G037221/1 and the STC Center for Integrated Quantum Materials (NSF Grant No. DMR-1231319). We thank H.S. Knowles andP. Borisova for technical assistance.

Author contributions

M.A., A.C.F. and M.L. managed the project. M.A., D.M.K. and C.P.-B. devised the project. P.L. fabricated the nanopatterned substrates. C.P.-B., A.R.-P.M. and M.B. developed the deposition method and fabricated the samples. C.P.-B., A.R.-P.M., D.M.K., A.K.O., M.B. and D.Y. performed the optical measurements. C.P.-B., A.R.-P.M. and M.B. performed AFM measurements. All authors participated in the results discussions, analysis and the writing of the manuscript.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Palacios-Berraquero, C. et al. Large-scale quantum-emitter arrays in atomically thin semiconductors. Nat. Commun. 8, 15093doi: 10.1038/ncomms15093 (2017).

Publishers note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional afliations.

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

r The Author(s) 2017

6 NATURE COMMUNICATIONS | 8:15093 | DOI: 10.1038/ncomms15093 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications