We report the observation of optical second harmonic generation (SHG) in single-layer indium selenide (InSe). We measure a second harmonic signal of > 103 cts/s under nonresonant excitation using a home-built confocal microscope and a standard pulsed pico-second laser. We demonstrate that polarization-resolved SHG serves as a fast, non-invasive tool to determine the crystal axes in single-layer InSe and to relate the sharp edges of the flake to the armchair and zigzag edges of the crystal structure. Our experiment determines these angles to an accuracy better than ± 0.2°. Treating the two-dimensional material as a nonlinear polarizable sheet, we determine a second-order sheet polarizability |χsheet(2)|=(17.9±11.0)×1020m2V1 for single-layer InSe, corresponding to an effective nonlinear susceptibility value of |χeff(2)|(223±138)×1012mV1 accounting for the sheet thickness (d ≈ 0.8 nm). We demonstrate that the SHG technique can also be applied to encapsulated samples to probe their crystal orientations. The method is therefore suitable for creating high quality van der Waals heterostructures with control over the crystal directions.

Since the isolation of graphene by mechanical exfoliation,1 leading, inter-alia, to the discovery of the quantum Hall effect in graphene,2 two-dimensional (2D) materials, such as atomic layers of transition metal dichalcogenides (TMDCs) (e.g. MoS2, MoSe2, WS2 and WSe2) and III-IV compounds (e.g. GaS, GaSe) have attracted great attention in materials research on account of their tunable electronic and optical properties.3,4 The ability to combine the 2D layers with hexagonal boron nitride (h-BN) and few-layer graphene offers an opportunity of creating high-performance, 2D opto-electronic devices.5,6

The desire to stack different 2D materials with precise control over the twist angle7–10 creates a need for a fast, non-invasive tool to probe the underlying crystal symmetries and crystallographic orientations. Nonlinear optical techniques, such as second harmonic generation (SHG), provide insight into the properties of surfaces or interfaces,11 particularly among non-centrosymmetric materials. In general, under an incident electric field E(ω) with fundamental angular frequency ω, the second-order nonlinear polarization is determined by a third-rank electric susceptibility tensor

χ(2):Pi(2ω)=ϵ0χijk(2)(2ω;ω,ω)Ej(ω)Ek(ω)
(1)

resulting in the creation of SHG. The second-order nonlinearities of a material can be probed by impinging intense linearly polarized light at angular frequency ω onto its surface and measuring the generated outgoing response at 2ω through an analyzer. By using different polarization combinations and/or by varying the orientation of the interface with respect to the incoming beam, different components of χ(2) can be determined, giving information about the structural symmetry of the material as well as the strength of various nonlinear processes. SHG has been observed in single- and few-layer TMDCs, such as MoS212,13 and WS2,14 and has been used to align the crystal axes of various TMDCs in order to create high-quality heterobilayers with strong interlayer exciton emission.15,16 However, TMDC flakes with an even number of layers exhibit inversion symmetry, resulting in a vanishing second-order nonlinearity (χ(2) = 0) which prevents the observation of SHG for all layer thicknesses.

Among the large family of van der Waals (vdW) crystals, indium selenide (InSe) has emerged as a promising 2D semiconductor due to its highly tunable optical response in the near-infrared to the visible spectrum17 and its high electron mobilities at room and liquid-helium temperatures allowing the quantum Hall effect to be observed.18 Similar to the TMDCs, single-layer InSe belongs to the non-centrosymmetric D3h (6¯m2) point group with only one independent non-zero second-order nonlinear susceptibility tensor element. The susceptibility components satisfy χxxx(2)=χxyy(2)=χyyx(2)=χyxy(2) with x along the armchair direction (Fig. 1 (right)).11 In contrast, the specific stacking order in bulk and few-layer InSe breaks the mirror-plane symmetry characteristic of single-layer InSe, thus maintaining broken inversion symmetry for all layer thicknesses. This provides an opportunity to investigate nonlinearities in thin InSe films layer-by-layer. The dependence of polarization-resolved SHG on the crystallographic axes further opens up an optical means of characterizing the crystal structure and orientation of the thin InSe films. Indeed, observations of SHG in relatively thin InSe sheets (from 9 to 25 nm) with even and odd number of layers have been reported recently.19 However, in the single-layer limit and under non-resonant condition, it remains elusive whether SHG signal can still be detected. In this Letter, we present SHG from encapsulated single- and few-layer InSe with the aim of determining the crystal axes and estimating the effective second-order nonlinear susceptibility of the single-layer.

FIG. 1.

Left: Optical image of the InSe flake. The indicated layer thicknesses were determined by optical contrast and verified by atomic force microscopy and photoluminescence measurements. Second harmonic generation measurements were carried out on the parts of the InSe flake which are fully encapsulated in h-BN. The device consists of SiO2 (290 nm)/h-BN (46 nm)/InSe/h-BN (8 nm). Right: Schematic of the InSe crystal structure (top and side view). Purple and orange spheres correspond to indium (In) and selenium (Se) atoms, respectively.

FIG. 1.

Left: Optical image of the InSe flake. The indicated layer thicknesses were determined by optical contrast and verified by atomic force microscopy and photoluminescence measurements. Second harmonic generation measurements were carried out on the parts of the InSe flake which are fully encapsulated in h-BN. The device consists of SiO2 (290 nm)/h-BN (46 nm)/InSe/h-BN (8 nm). Right: Schematic of the InSe crystal structure (top and side view). Purple and orange spheres correspond to indium (In) and selenium (Se) atoms, respectively.

Close modal

Thin unprotected InSe films have optical properties which deteriorate over time, interpreted as a gradual degradation of the crystal due to interaction with oxygen and water in the atmosphere.18,20 Thus, to protect the material from the interaction with the environment, we employed exfoliation and subsequent encapsulation of single- and few-layer InSe in an inert (argon) atmosphere. The resulting InSe structures are stable under ambient conditions even in the single-layer limit.21 Fig. 1 (left) shows an optical image of the sample used to investigate the SHG response of thin InSe flakes mechanically exfoliated from a Bridgman-grown rhombohedral γ-InSe bulk crystal. The exfoliated InSe layers were picked up by a h-BN crystal (8 nm) and deposited onto a large, relatively thick h-BN flake (46 nm) prepared on a Si wafer with 290 nm thick SiO2. InSe layers of different thickness and single-layer steps (≈0.8 nm) were identified by atomic force microscopy and photoluminescence measurements after encapsulation.

SHG spectroscopy was performed at room temperature using a home-built, confocal microscope setup. The optical pulse, centered at a wavelength of 810 nm, was obtained from a Ti:Sapphire laser with 76 MHz repetition rate. All the spectral components of the ∼ 150 fs pulse were retained; dispersion in the optical fiber connecting the laser source to the microscope stretches the pulse to the pico-second domain at the sample (intensity FWHM 36 ± 10 ps). The laser (average power 3.2 mW) was focused to a spot size of about 1.5 μm on the sample by a microscope objective lens (40x, NA = 0.65) at normal incidence and with a fixed linear polarization. The SHG signal was collected by the same objective and directed through a dichroic beamsplitter to a spectrometer equipped with a 300 grooves/mm grating and a nitrogen-cooled silicon charge-coupled device (CCD).

As shown in Fig. 2, an easy-to-measure SHG signal (>103 cts/s) appeared at wavelength 405 nm when the laser beam (wavelength 810 nm) was focussed on encapsulated single- and few-layer InSe. Moreover, the SHG response could be observed for both even and odd number of layers. No observable SHG signal was measured on h-BN. The integrated SHG signal increases by more than two orders of magnitude as the thickness of the InSe flake increases from 1 to 10 layers (inset (a) to Fig. 2). Specifically, the SHG signal depends approximately quadratically on the layer thickness, i.e. ISHG(2ω) ∝ d2, as shown in the inset (a) in Fig. 2 where the integrated SHG intensity ISHG(2ω) of the InSe flakes is plotted versus the layer number N (with d ≈ N × 0.8 nm) on a logarithmic scale. Following,22 the dependence of the SHG intensity on the flake thickness is given by

ISHG(2ω)(lc|χeff(2)|)2sin2πd2lc
(2)

where d is the flake thickness, lc = 2πk ≈ 20 μm is the coherence length and |χeff(2)| is the effective nonlinear susceptibility determined by the sample geometry and the non-zero components of the nonlinear sheet susceptibility tensor χsheet(2) (|χeff(2)|=|χsheet(2)|/d). Thus, in the limit of atomically thin films (i.e. d ≪ lc), the SHG intensity ISHG(2ω) is expected to increase quadratically with the number of layers, in good agreement with our observations (inset (a) to Fig. 2).

FIG. 2.

SHG spectra of single- and few-layer InSe at room temperature. The inset (a) (plotted on a double logarithmic scale) shows the approximately quadratic increase (slope of ∼ 1.8) of the integrated SHG with the number of layers. The error bars in (a) are smaller than the symbol size. The inset (b) shows the SHG spectrum of encapsulated single-layer InSe.

FIG. 2.

SHG spectra of single- and few-layer InSe at room temperature. The inset (a) (plotted on a double logarithmic scale) shows the approximately quadratic increase (slope of ∼ 1.8) of the integrated SHG with the number of layers. The error bars in (a) are smaller than the symbol size. The inset (b) shows the SHG spectrum of encapsulated single-layer InSe.

Close modal

The relatively strong SHG signal facilitates an investigation of the lattice symmetry and crystallographic orientation of the thin flakes. The SHG intensity ISHG(2ω) is strongly dependent on the polarization angle (θθ0) between the laser polarization E(ω) and the armchair direction of the crystal defined in Fig. 3(b). For polarization-resolved SHG, the laser polarization was rotated about the z-axis by a half-wave plate to vary θ, and the SHG signal ISHG(2ω) (SHG polarization parallel to the excitation polarization) was collected using an analyzer located in front of the detector. The polar plot for single-layer InSe in Fig. 3(a) shows a strongly varying 6-fold symmetry of ISHG(2ω)cos23(θθ0). This directly reveals the underlying symmetry and orientation of the single-layer InSe flake. The initial orientation of the sample with respect to the armchair direction of the crystal was determined to be θ0 = 28.1° ± 0.2°. The sharp edges along which the single-layer InSe flake cleaved during exfoliation could therefore be clearly assigned to the armchair and zigzag crystal axes, respectively (Fig. 3(b)).

FIG. 3.

SHG characterization of encapsulated InSe. (a) Polar plot of the “parallel” SHG intensity ISHG(2ω) of single-layer InSe as a function of rotation angle θ. Fitting the angular dependence to ISHG(2ω)cos23(θθ0) (solid purple line), the armchair direction (dark blue arrow) of the sample is determined as θ0 = 28.1° ± 0.2°. This armchair direction is 30° from the zigzag direction (light blue arrow). (b) Optical image of the InSe flakes before encapsulation. The armchair direction forms an angle θ0 with respect to the “laboratory” axis (LAB). The indicated armchair (dark blue arrow) and zigzag (light blue arrow) edges were determined by polarization-resolved SHG performed on the encapsulated sample.

FIG. 3.

SHG characterization of encapsulated InSe. (a) Polar plot of the “parallel” SHG intensity ISHG(2ω) of single-layer InSe as a function of rotation angle θ. Fitting the angular dependence to ISHG(2ω)cos23(θθ0) (solid purple line), the armchair direction (dark blue arrow) of the sample is determined as θ0 = 28.1° ± 0.2°. This armchair direction is 30° from the zigzag direction (light blue arrow). (b) Optical image of the InSe flakes before encapsulation. The armchair direction forms an angle θ0 with respect to the “laboratory” axis (LAB). The indicated armchair (dark blue arrow) and zigzag (light blue arrow) edges were determined by polarization-resolved SHG performed on the encapsulated sample.

Close modal

To quantify the nonlinear response of single-layer InSe, we follow the formalism of Ref. 23 where SHG from an ultrathin layer is treated as radiation from a nonlinear, polarizable sheet embedded in a dielectric medium, with boundary effects taken into account by Fresnel transmission coefficients. The second-order susceptibility |χsheet(2)| of InSe can be extracted from measurements of the intensity of the SHG with respect to the driving intensity via

ISHG(2ω)=ω22c3ϵ0tin4(ω)tout2(2ω)|χsheet(2)|2I2(ω)
(3)

where tin(ω) and tout(2ω) account for the local field correction factors due to the dielectric environment determined by transfer matrix methods (“Essential Macleod”). Relating the laser intensity I(ω) and the SHG intensity ISHG(2ω) in Eq. (3) to the experimentally measured time-averaged power values, we obtain |χsheet(2)|=(17.9±11.0)×1020m2V1 for single-layer InSe. To compare to other nonlinear optical bulk materials, we estimate an effective bulk-like nonlinear susceptibility |χeff(2)|=(223±138)×1012mV1 by including the thickness of the 2D material (d ≈ 0.8 nm). This value is similar to the strong second-order optical susceptibility measured for single- and few-layer GaSe.24 

In conclusion, we report an observation of SHG in single-layer InSe under nonresonant excitation, yielding a nonlinear sheet susceptibility with an estimated value of |χeff(2)|223pmV1. Quantitative characterization of the nonlinear response of single- and few-layer InSe reveals the expected quadratic dependence of the SHG signal on the number of layers. The crystalline symmetry was probed by polarization-resolved SHG where the “petal” direction with maximum signal is parallel to the in-plane In-Se or Se-In (armchair) direction. This allowed for a fast and precise (±0.2°) assignment of the sharp edges of the InSe flake to its crystal axes, demonstrating that SHG serves as useful tool for the determination of the orientation of the material’s crystallographic axes. The ability to distinguish different crystallographic axes, even in encapsulated samples, can be exploited to improve the quality of van der Waals heterostructures by stacking various 2D materials with precise twist-angle control.

We thank Immo Söllner, Jan-Philipp Jahn, and Daniel Najer for technical assistance and experiments; Amalia Patanè from the School of Physics and Astronomy, University of Nottingham, for providing the bulk InSe crystals. The work in Basel was financially supported by SNF [Project Number 200020_175748], Swiss Nanoscience Institute (SNI), QCQT PhD School and NCCR QSIT. The work in Manchester was financially supported by the Engineering and Physical Sciences Research Council [Grant Number EP/M012700/1] (EPSRC). R.V.G. acknowledges financial support from the Royal Society Fellowship Scheme. M.H. and D.T. acknowledge support from the EPSRC NowNano Doctoral Training Centre.

1.
K. S.
Novoselov
,
A. K.
Geim
,
S. V.
Morozov
,
D.
Jiang
,
Y.
Zhang
,
S. V.
Dubonos
,
I. V.
Grigorieva
, and
A. A.
Firsov
,
Science
306
,
666
(
2004
).
2.
Y.
Zhang
,
Y.-W.
Tan
,
H. L.
Stormer
, and
P.
Kim
,
Nature
438
,
201
(
2005
).
3.
F.
Xia
,
H.
Wang
,
D.
Xiao
,
M.
Dubey
, and
A.
Ramasubramaniam
,
Nat. Photonics
8
,
899
(
2014
).
4.
K. F.
Mak
and
J.
Shan
,
Nat. Photonics
10
,
216
(
2016
).
5.
C. R.
Dean
,
A. F.
Young
,
I.
Meric
,
C.
Lee
,
L.
Wang
,
S.
Sorgenfrei
,
K.
Watanabe
,
T.
Taniguchi
,
P.
Kim
,
K. L.
Shepard
, and
J.
Hone
,
Nat. Nanotechnol.
5
,
722
(
2010
).
6.
A. V.
Kretinin
,
Y.
Cao
,
J. S.
Tu
,
G. L.
Yu
,
R.
Jalil
,
K. S.
Novoselov
,
S. J.
Haigh
,
A.
Gholinia
,
A.
Mishchenko
,
M.
Lozada
,
T.
Georgiou
,
C. R.
Woods
,
F.
Withers
,
P.
Blake
,
G.
Eda
,
A.
Wirsig
,
C.
Hucho
,
K.
Watanabe
,
T.
Taniguchi
,
A. K.
Geim
, and
R. V.
Gorbachev
,
Nano Lett.
14
,
3270
(
2014
).
7.
A.
Mishchenko
,
J. S.
Tu
,
Y.
Cao
,
R. V.
Gorbachev
,
J. R.
Wallbank
,
M. T.
Greenaway
,
V. E.
Morozov
,
S. V.
Morozov
,
M. J.
Zhu
,
S. L.
Wong
,
F.
Withers
,
C. R.
Woods
,
Y.-J.
Kim
,
K.
Watanabe
,
T.
Taniguchi
,
E. E.
Vdovin
,
O.
Makarovsky
,
T. M.
Fromhold
,
V. I.
Fal’ko
,
A. K.
Geim
,
L.
Eaves
, and
K. S.
Novoselov
,
Nat. Nanotechnol.
9
,
808
(
2014
).
8.
J.
Kunstmann
,
F.
Mooshammer
,
P.
Nagler
,
A.
Chaves
,
F.
Stein
,
N.
Paradiso
,
G.
Plechinger
,
C.
Strunk
,
C.
Schüller
,
G.
Seifert
,
D. R.
Reichman
, and
T.
Korn
,
Nat. Phys.
14
,
801
(
2018
).
9.
K.
Liu
,
L.
Zhang
,
T.
Cao
,
C.
Jin
,
D.
Qiu
,
Q.
Zhou
,
A.
Zettl
,
P.
Yang
,
S. G.
Louie
, and
F.
Wang
,
Nat. Commun.
5
,
4966
(
2014
).
10.
A. M.
van der Zande
,
J.
Kunstmann
,
A.
Chernikov
,
D. A.
Chenet
,
Y.
You
,
X.
Zhang
,
P. Y.
Huang
,
T. C.
Berkelbach
,
L.
Wang
,
F.
Zhang
,
M. S.
Hybertsen
,
D. A.
Muller
,
D. R.
Reichman
,
T. F.
Heinz
, and
J. C.
Hone
,
Nano Lett.
14
,
3869
(
2014
).
11.
Y. R.
Shen
,
The principles of nonlinear optics
(
Wiley
,
New York
,
2003
).
12.
N.
Kumar
,
S.
Najmaei
,
Q.
Cui
,
F.
Ceballos
,
P. M.
Ajayan
,
J.
Lou
, and
H.
Zhao
,
Phys. Rev. B
87
,
161403
(
2013
).
13.
L. M.
Malard
,
T. V.
Alencar
,
A. P. M.
Barboza
,
K. F.
Mak
, and
A. M.
de Paula
,
Phys. Rev. B
87
,
201401
(
2013
).
14.
C.
Janisch
,
Y.
Wang
,
D.
Ma
,
N.
Mehta
,
A. L.
Elías
,
N.
Perea-López
,
M.
Terrones
,
V.
Crespi
, and
Z.
Liu
,
Sci. Rep.
4
,
5530
(
2014
).
15.
P.
Rivera
,
J. R.
Schaibley
,
A. M.
Jones
,
J. S.
Ross
,
S.
Wu
,
G.
Aivazian
,
P.
Klement
,
K.
Seyler
,
G.
Clark
,
N. J.
Ghimire
,
J.
Yan
,
D. G.
Mandrus
,
W.
Yao
, and
X.
Xu
,
Nat. Commun.
6
,
6242
(
2015
).
16.
P. K.
Nayak
,
Y.
Horbatenko
,
S.
Ahn
,
G.
Kim
,
J.-U.
Lee
,
K. Y.
Ma
,
A.-R.
Jang
,
H.
Lim
,
D.
Kim
,
S.
Ryu
,
H.
Cheong
,
N.
Park
, and
H. S.
Shin
,
ACS Nano
11
,
4041
(
2017
).
17.
G. W.
Mudd
,
M. R.
Molas
,
X.
Chen
,
V.
Zólyomi
,
K.
Nogajewski
,
Z. R.
Kudrynskyi
,
Z. D.
Kovalyuk
,
G.
Yusa
,
O.
Makarovsky
,
L.
Eaves
,
M.
Potemski
,
V. I.
Fal’ko
, and
A.
Patanè
,
Sci. Rep.
6
,
39619
(
2016
).
18.
D. A.
Bandurin
,
A. V.
Tyurnina
,
G. L.
Yu
,
A.
Mishchenko
,
V.
Zólyomi
,
S. V.
Morozov
,
R. K.
Kumar
,
R. V.
Gorbachev
,
Z. R.
Kudrynskyi
,
S.
Pezzini
,
Z. D.
Kovalyuk
,
U.
Zeitler
,
K. S.
Novoselov
,
A.
Patanè
,
L.
Eaves
,
I. V.
Grigorieva
,
V. I.
Fal’ko
,
A. K.
Geim
, and
Y.
Cao
,
Nat. Nanotechnol.
12
,
223
(
2016
).
19.
S.
Deckoff-Jones
,
J.
Zhang
,
C. E.
Petoukhoff
,
M. K. L.
Man
,
S.
Lei
,
R.
Vajtai
,
P. M.
Ajayan
,
D.
Talbayev
,
J.
Madéo
, and
K. M.
Dani
,
Sci. Rep.
6
,
22620
(
2016
).
20.
O. D.
Pozo-Zamudio
,
S.
Schwarz
,
J.
Klein
,
R. C.
Schofield
,
E. A.
Chekhovich
,
O.
Ceylan
,
E.
Margapoti
,
A. I.
Dmitriev
,
G. V.
Lashkarev
,
D. N.
Borisenko
,
N. N.
Kolesnikov
,
J. J.
Finley
, and
A. I.
Tartakovskii
, arXiv:1506.05619 (
2015
).
21.
Y.
Cao
,
A.
Mishchenko
,
G. L.
Yu
,
E.
Khestanova
,
A. P.
Rooney
,
E.
Prestat
,
A. V.
Kretinin
,
P.
Blake
,
M. B.
Shalom
,
C.
Woods
,
J.
Chapman
,
G.
Balakrishnan
,
I. V.
Grigorieva
,
K. S.
Novoselov
,
B. A.
Piot
,
M.
Potemski
,
K.
Watanabe
,
T.
Taniguchi
,
S. J.
Haigh
,
A. K.
Geim
, and
R. V.
Gorbachev
,
Nano Lett.
15
,
4914
(
2015
).
22.
A.
Yariv
and
P.
Yeh
,
Optical waves in crystals: propagation and control of laser radiation
(
Wiley
,
New York
,
1984
) Chap. 12.
23.
Y. R.
Shen
,
Annu. Rev. Phys. Chem.
40
,
327
(
1989
).
24.
Y.
Tang
,
K. C.
Mandal
,
J. A.
McGuire
, and
C. W.
Lai
,
Phys. Rev. B
94
,
125302
(
2016
).