Self-consistent dielectric constant determination for monolayer WSe 2

Frequency-dependent dielectric constant dispersion of monolayer WSe2, ( ) ( ) ( ) 1 2 ε ε ε i ω ω ω = + , was obtained from simultaneously measured transmittance and reflectance spectra. Optical transitions of the trion as well as A-, B-, and C-excitons are clearly resolved in the 2 ε spectrum. A consistent Kramers-Kronig transformation between the 1 ε and 2 ε spectra support the validity of the applied analysis. It is found that the Aand B-exciton splitting in the case of the double-layer WSe2 can be attributed to the spin-orbit coupling, which is larger than that in the monolayer WSe2. In addition, the temperatureinduced evolution of the A-exciton energy and its width are explained by model equations with electron-phonon interactions. © 2018 Optical Society of America under the terms of the OSA Open Access Publishing Agreement


Introduction
In contrast to graphene and boron nitride that can be categorized as a gapless Dirac material and wide gap insulator, respectively [1,2], transition-metal dichalcogenides (TMDCs) such as MoS 2 , MoSe 2 , WS 2 , and WSe 2 , are semiconducting layered materials.Owing to the unique electronic and optical properties that TMDCs possess, they have attracted great interest from many fields [3][4][5].Because of their finite energy gap, TMDCs can outperform graphenebased transistors in terms of the on/off ratio [6].The large spin-orbit splitting in the case of MoS 2 or WSe 2 has led to long sustaining spin properties enabling novel functionalities in the field of valleytronics [7,8].
Because monolayer MoS 2 or WSe 2 are direct gap semiconductors with good radiative recombination efficiency, they can be utilized as light sources of atomic length scale.Due to large exciton binding energy light emission in TMDCs is featured by exciton states even at room temperature [4,9].In the case of bulk crystals, the exciton binding energy can be determined using the following equation: , where * μ is the reduced effective mass of electron and hole quasi-particles and ε is the dielectric constant [10].
However, in case of monolayers, the substantial field intrusion out of the host material leads to a complicated change in the binding energy; it has been shown that a series of exciton sublevels in monolayer MoS 2 do not follow the Rydberg-type energy separation of the hydrogen atom, but they can be explained by field intrusion into the substrate [11].Nevertheless, the dielectric constant and effective mass are still the most deterministic factors for the binding energy of two-dimensional (2D) monolayer excitons.While the effective mass can be obtained by measuring the energy band dispersion using angle-resolved photoemission spectroscopy [12,13], the dielectric constant can be obtained through optical characterizations.
The imaginary part of the dielectric constant, ( ) ε ω , which is proportional to the optical transition probability, provides valuable information about the density of energy states, including the excitons.Though the ellipsometry technique is used to obtain complex dielectric constants of monolayer MoS 2 and WSe 2 [14][15][16], it not only requires dispersion modeling, but also has difficulty in performing microscopic measurements.It has been proposed that the absorbance is proportional to the reflectance contrast in the case of atomically thin layers of graphene [17], In the opaque spectral range, however, the complex dielectric constant cannot be completely determined based on the reflectance or transmittance individually, but values of both quantities are necessary for singular extractions.Several works have reported on the determination of the dielectric constant of monolayer TMDCs based on their reflectance and transmittance [18,19].In particular, Morozov et al. reported on the dielectric constants of monolayer TMDCs obtained using concerted transmittance and reflectance measurements at room temperature [19].
In this study, we report on the frequency-dependent dielectric constants of monolayer WSe 2 .By using the transfer matrix method (TMM), both real and imaginary components of the dielectric constant are obtained from the simultaneously measured reflectance and transmittance spectra.In addition, we study the temperature-induced evolution of the Aexciton and compare the dielectric constant of monolayer with that of the double layer WSe 2 .

Sample preparation and methods
Monolayer WSe 2 was synthesized on a SiO 2 (300 nm)/Si substrate using the chemical vapor deposition (CVD) method.The synthesized monolayer film was then transferred onto a quartz plate via the poly(methyl methacrylate) (PMMA) support method in order to perform the transmittance measurement.The flake size of the monolayer was typically ~100 μm 100 μm × , but, it depended on the growth conditions.As indicated by the optical image in Fig. 1(a), in addition to pure monolayers, some double-layers were formed as well.A broadband white light source from a supercontinuum laser (Compact, NKT photonics) was passed through a monochromator, and the filtered beam was then focused using an objective lens onto the monolayer or double layer region of WSe 2 .Both the transmitted (T) and reflected (R) intensities were simultaneously measured using silicon photodiodes.In addition, the q R and q T at adjacent quartz regions were measured, using which, the normalized reflectance ( n R ) and transmittance ( n T ), defined by n R = / q R R and n T = / q T T , respectively, were obtained.The sample was mounted inside an optical cryostat in vacuum and the lock-in technique referenced by an optical chopper was applied to suppress the detection noise.

Results and discussion
Figure 2 shows the n R and n T spectra of monolayer WSe 2 measured at a temperature of 10 K.While n T is less than unity because of the absorption by the WSe 2 layer, n R is greater than one over the entire range.In previous works, it has been demonstrated that the reflectance contrast, defined as ( )/ , is positive and proportional to the absorbance in the case of graphene [17], indicating that n R is more than one in opaque regions.In many studies on monolayer TMDCs, the reflectance was analyzed to represent the optical transition probability or the absorbance.This can be attributed, in part, to the fact that transmittance could not be measured in samples with a silicon substrate.However, as shown in Fig. 2, there is a difference in the spectral shape of the n R and n T spectra.Furthermore, as shown in the inset of Fig. 2, the A-exciton peak near 1.75 eV in n T is slightly differentiated compared with that in n R , which indicates that the reflectance alone cannot provide the precise position and shape of the exciton resonance, especially when the resonance is spectrally narrow.
In an optical system under study, the reflectance or transmittance is dependent on the dielectric constants of the constituent layers.In the opaque spectral region, the dielectric constant cannot be set based on reflectance alone.However, it can be unambiguously determined using both the n R and n T spectra.As previously mentioned, in order to obtain the dielectric constant from the measured n R and n T spectra, we applied the TMM approach [20][21][22].According to the TMM approach, the matrix elements ABCD representing the light propagation from air to quartz through the WSe 2 layer are presented as the multiples of the three sub matrices in Eq. ( 1), which correspond to the electric field transfer at the air/WSe 2 interface, through the WSe 2 layer, and at the WSe 2 /quartz interface, respectively.
Here, ω is the angular frequency and c is the speed of light.In addition, ε and n q are the dielectric constants of WSe 2 and refractive index of quartz, respectively.With the ABCD elements known, the reflectance and transmittance can be obtained using . In this study, based on the measured n R and n T spectra, the TMM approach was used in reverse to obtain the dielectric constant of the WSe 2 layers.According to previous literature, d = 0.65 nm was assumed as the monolayer thickness [23].We prepared a 2D array of ( 1ε , 2 ε ), covering a broad range of real ( 1 ε ) and imaginary ( 2 ε ) components of the dielectric constant, and calculated the ( , ) nc nc R T for each set.Then, for each photon energy, a set of ( 1 ε , 2 ) ε was selected for which the deviation of ( , )  ε ω attains a value of ~21.1 in the zero-frequency limit.We find that the low-energy 1 ε spectrum, which is indicated using a red line in Fig. 3(c), can be modeled by the Sellmeier single oscillator formula of ( ) as a function of photon energy E in the eV unit [26].Here, the energy of the Sellmeier oscillator is differentiated from the compensational oscillator because the optical transitions within our detection range contributes as well in the Kramers-Kronig transformation.
According to Eq. ( 1) the extracted dielectric constant depends on the choice of the monolayer thickness d.As we varied the layer thickness from 0.5 to 0.8 nm it was found that the value of ε d × did not change noticeably.We note that this multiplied value, through the relation of 2 0 σ ωεε d with 0 ε being the vacuum permittivity, corresponds to the 2D optical conductivity, which are widely used in order to explain the optical response of atomically thin monolayers [19,27,28].The ( n R , n T ) measurement as well as the ( 1 ε , 2 ) ε extraction were repeated on the doublelayer WSe 2 region assuming a thickness of 2d .Figure 4 shows a comparison of the ( 1 ε , 2 ) ε spectra of the double-layer WSe 2 with those of the monolayer WSe 2 at 10 K.It can be observed that, for the double-layer WSe 2 , the energies of the A-and B-excitons are 36 meV and 8 meV lower than those of the monolayer WSe 2 , respectively.The energy separation between the A-and B-excitons, which is larger in the case of the double-layer WSe 2 , originates from the spin-orbit splitting of the valence band.It should be noted that the larger A-B splitting in the case of the double-layer over the monolayer was previously explained by the effect of interlayer interactions [29,30].In contrast to the A-and B-excitons, the Cexciton experiences a greater energy shift of ~160 meV as the WSe 2 changes from monolayer to double-layer.There have been some controversial theories about the origin of the Cexcitons-whether they are due to the transition between parallel bands or Van Hove singularities at saddle points in the band structure [31,32].We believe that the layer number dependence of C-exciton energy in combination with the theoretical interpretation might help in addressing this issue.Though the 1 ε spectrum of the double-layer in Fig. 4 is accordingly differentiated from the monolayer, especially near the exciton resonances, the overall shape, including the low energy value of 1 ~26 ε = at 1.5 eV, is similar to the monolayer.As the temperature is increased, the electron-phonon interaction is accelerated and results in spectral broadening of the exciton resonances.We investigated the temperature dependence of excitons by measuring the dielectric constant at different temperatures.As is expected, the evolution of the 2 ε spectrum of the monolayer WSe 2 in Fig. 5 reveals a red shift and broadening of the A-, B-, and C-excitons with an increase in temperature.We plotted the temperature induced variation of A-exciton energy and width in the inset.From data fitting, we find that the energy shift follows the Varshni formula of

Conclusion
In this study, by performing TMM analysis of simultaneously measured transmittance and reflectance, frequency-dependent dielectric constant of monolayer WSe 2 was obtained with self-consistency.In addition, thermal shift and broadening of A-exciton resonance acquired from temperature-dependent measurements demonstrates significant exciton relaxation because of electron-phonon interactions.We expect that self-consistent dielectric constant extraction performed in this study will be useful in characterizing other ultrathin layered materials and provide valuable optical data for TMDC-based optoelectronic devices.

Fig. 1 .
Fig. 1.Schematic diagram of the experimental setup for simultaneous reflection and transmission measurements.White light from the supercontinuum laser passes the monochromator and is focused onto WSe 2 layers using an objective lens.Transmitted or reflected beams are simultaneously measured using photodiodes.The optical microscope image of the WSe 2 layers is included in the figure.

Fig. 2 .
Fig. 2. Normalized reflectance ( n R ) and transmittance ( T ) obtained in monolayer WSe 2 at 80 K.The inset shows the comparison of n R and n T spectra around the sharp peak at 1.75 eV.

Fig. 3 . 2 ε 1 ε− 1 ε− 2 ε 2 ε 2 ε 2 ε ω at 3 .
Fig. 3. Real (a) and imaginary (b) spectra of the complex dielectric constant as a function of the photon energy for monolayer WSe 2 at 80 K with the plots in the right highlighting the region near the A-exciton position.The notations A, B, and C in ( ) 2 ε ω indicate positions of A-, B-, and C-excitons, respectively.The Kramers-Kronig transformation of the ( ) 2 ε ω ,

Fig. 4 . 1 ε
Fig. 4. 1 ε and 2 ε spectra of the double-layer WSe 2 at 80 K compared with those of the monolayer WSe 2 plotted with solid and dotted lines, respectively.
. The temperature induced broadening, on the other hand, is well described by the Rudin's relation of

Fig. 5 . 2 ε
Fig. 5. 2 ε spectrum at different temperatures from 10 to 290 K. Inset shows the temperature dependence of the energy and width of A-exciton resonance, where the dotted lines represent the model fitting of the energy and width according to the Varshni formula and Rudin's relation, respectively.

The 2 ε
spectrum exhibits clear optical transitions at A-, B-, and C-excitons as well as the trion resonance.The Kramers-Kronig transformation leads to a dispersion model ( ) region below 1.4 eV.