Search tips
Search criteria 


Logo of ncommsLink to Publisher's site
Nat Commun. 2017; 8: 893.
Published online 2017 October 12. doi:  10.1038/s41467-017-00749-4
PMCID: PMC5715017

Ultra-strong nonlinear optical processes and trigonal warping in MoS2 layers


Nonlinear optical processes, such as harmonic generation, are of great interest for various applications, e.g., microscopy, therapy, and frequency conversion. However, high-order harmonic conversion is typically much less efficient than low-order, due to the weak intrinsic response of the higher-order nonlinear processes. Here we report ultra-strong optical nonlinearities in monolayer MoS2 (1L-MoS2): the third harmonic is 30 times stronger than the second, and the fourth is comparable to the second. The third harmonic generation efficiency for 1L-MoS2 is approximately three times higher than that for graphene, which was reported to have a large χ (3). We explain this by calculating the nonlinear response functions of 1L-MoS2 with a continuum-model Hamiltonian and quantum mechanical diagrammatic perturbation theory, highlighting the role of trigonal warping. A similar effect is expected in all other transition-metal dichalcogenides. Our results pave the way for efficient harmonic generation based on layered materials for applications such as microscopy and imaging.


Nonlinear optical phenomena can generate high-energy photons by converting n = 2, 3, 4,… low-energy photons into one high-energy photon. These are usually referred to as second-, third-, and fourth-harmonic generation (SHG, THG, and FHG)1. Due to different selection rules1, 2, harmonic processes are distinct from optically pumped laser phenomena (e.g., optically pumped amplification3), and other typical single-photon processes (e.g., single-photon excited photoluminescence1), in which the energy of the generated photons is smaller than the pump photons. Multiphoton harmonic processes have been widely exploited for various applications (e.g., all-optical signal processing in telecommunications1, 4, medicine5, and data storage6), as well as to study various transitions forbidden under low-energy single-photon excitation5, 6. The physical origin of these processes is the nonlinear polarization induced by an electromagnetic field E. This gives rise to higher harmonic components, the n-th harmonic component amplitude being proportional1 to |E|n. Quantum mechanically, higher-harmonic generation involves the annihilation of n pump photons and generation of a photon with n times the pump energy. Because an n-th order nonlinear optical process requires n photons to be present simultaneously, the probability of higher-order processes is lower than that of lower order1. Thus, higher-order processes are typically weaker and require higher pump intensities7, 8.

Graphene and related materials are at the center of an ever-increasing research effort due to their unique and complementary properties, making them appealing for a wide range of photonic and optoelectronic applications911. Among these, semiconducting transition-metal dichalcogenides (TMDs) are of particular interest due to their direct bandgap when in monolayer (1L) form12, leading to an increase in luminescence by a few orders of magnitude compared with the bulk material12, 13. 1L-MoS2 has a single layer of Mo atoms sandwiched between two layers of S atoms in a trigonal prismatic lattice. Therefore, in contrast to graphene, it is non-centrosymmetric and belongs to the space group D3h1 14. The lack of spatial inversion symmetry makes 1L-MoS2 an interesting material for nonlinear optics, since second-order nonlinear processes are present only in non-centrosymmetric materials1. However, when stacked, MoS2 layers are arranged mirrored with respect to one another14, therefore MoS2 with an even number of layers (EN) is centrosymmetric and belongs to the D3d3 space group14, producing no second harmonic (SH) signal. On the other hand, MoS2 with odd number of layers (ON) is non-centrosymmetric. SHG from 1L-MoS2 was reported by several groups1421.

Here we present combined experimental and theoretical work on nonlinear harmonic generation in 1L and few-layer (FL) MoS2. We report strong THG and FHG from 1L-MoS2. THG is more than one order of magnitude larger than SHG, while FHG has the same magnitude as SHG. This is surprising, since one normally expects the intensity of nonlinear optical processes to decrease with n 1, 2, with the SHG intensity much larger than that in THG and FHG, although even-order processes only exist in non-centrosymmetric materials. Our results show that this expectation is wrong in the case of 1L-MoS2. At sufficiently low photon frequencies (in our experiments the photon energy of the pump is 0.8 eV), SHG only probes the low-energy band structure of 1L-MoS2. This is nearly rotationally invariant2229, but with corrections due to trigonal warping. It is because of these corrections23, 26, 27, fully compatible with the D3h1 space group1, but reducing the full rotational symmetry of the low-energy bands to a three-fold rotational symmetry1, that a finite amplitude of nonlinear harmonic processes can exist at low photon energies in EN-MoS2. The lack of spatial inversion symmetry is a necessary but not sufficient condition for the occurrence of SHG. A purely isotropic band structure gives a vanishing SHG signal3033, despite some terms in the Hamiltonian explicitly breaking inversion symmetry27, 3437. Terms proportional to the σ z Pauli matrix break inversion symmetry. Breaking the continuous rotational symmetry of isotropic models (e.g., by including trigonal warping) is required to obtain a non-zero second-order response in a two-band system. In hexagonal lattices, trigonal warping is a deviation from purely isotropic bands that emerges as one moves away from the corners K and K′ of the Brillouin zone23, 26, 27, 3638. Since the lattice has a honeycomb structure, this distortion displays a three-fold rotational symmetry23, 26, 27, 3638. We demonstrate that the observed THG/SHG intensity ratio can be explained by quantum mechanical calculations based on finite-temperature many-body diagrammatic perturbation theory39 and low-energy continuum-model Hamiltonians that include trigonal warping35. We conclude that, similar to SHG1418, the THG process is sensitive to the number of layers, their symmetry, relative orientation, as well as the elliptical polarization of the excitation light. Similar effects are expected for all other TMDs. This paves the way for the assembly of heterostructures with tailored nonlinear optical properties.



MoS2 flakes are produced by micromechanical cleavage (MC) of bulk MoS2 40 onto Si + 285 nm SiO2. 1L-MoS2 and bilayer (2L-MoS2) flakes are identified by a combination of optical contrast41, 42 and Raman spectroscopy43. Raman spectra are acquired by a Renishaw micro-Raman spectrometer equipped with a 600 line/mm grating and coupled with an Ar+ ion laser at 514.5 nm. Figure 1a shows the MoS2 flakes studied in this work and their Raman signatures. A reference MC graphene sample is also prepared and placed on a similar substrate.

Fig. 1
Optical image and Raman spectra of the MoS2 flakes. a Optical micrograph with single-layer, bilayer, and five-layer areas marked by 1L, 2L, and 5L. b Raman spectra of the same sample

SHG and THG charcterization

Nonlinear optical measurements are carried out with the set-up shown in Fig. 2 44, 45. As excitation source, we use an erbium-doped mode-locked fiber laser with a 50 MHz repetition rate, maximum average power 60 mW, and pulse duration 150 fs, which yields an estimated pulse peak power of ~8 kW46. The laser beam is scanned with a galvo mirror and focused on the sample using a microscope objective. The back-scattered second and third harmonic signals are split into different branches using a dichroic mirror and then detected using photomultiplier tubes (PMTs). For two-channel detection, the light is split into two PMTs using a dichroic mirror with 562 nm cutoff. After the dichroic mirror, the detected wavelength range can be further refined using bandpass filters. The light can also be directed to a spectrometer (OceanOptics QE Pro-FL). The average power on the sample is kept between 10 and 28 mW with a typical measurement time ~5 μs, which prevents sample damage and enables high signal-to-noise-ratio, even with acquisition time per pixel in the μs range.

Fig. 2
Schematic diagram of multiphoton microscope. MLL, linearly polarized mode-locked fiber laser. VA, variable attenuator. QWP, quarter-wave plate used to study the dependence of SHG and THG on the elliptical polarization of the pump light. BP filter, bandpass ...

SHG and THG images of the MoS2 sample are shown in Fig. 3a, b. The SH photon energy is ~1.6 eV, lower than the bandgap of 1L-MoS2 12, 13. This is not unexpected, as harmonic generation can occur when the harmonic energy is below the bandgap1, 47, 48. The SHG signal is generated in 1L-MoS2, while 2L-MoS2 appears dark. As discussed above, the second-order nonlinear response is present in 1L-MoS2, which is non-centrosymmetric. However, when stacked to form 2L-MoS2, MoS2 layers are mirrored14, 15. Therefore, EN-MoS2 is centrosymmetric14, 15, and belongs to the D3d3 space group14, 15, producing no SHG signal. On the other hand, ON-MoS2 flakes are non-centrosymmetric14, 15.

Fig. 3
Multiphoton images of MoS2 flakes. a SHG and b THG map of the flakes in Fig. 1a. c Optical spectrum of the nonlinear signal from 1L-MoS2 with a peak irradiance ~30 GW cm−2

We note that strong THG is detected compared with SHG, even for 1L-MoS2, Fig. 3b. THG was previously reported for a 7L-MoS2 flake18, but here we see it down to 1L-MoS2. Reference 49 followed our work50 and reported THG and SHG from 1L-MoS2, giving effective bulk-like second- and third-order susceptibilities χeff(2) and χeff(3) of 2.9 × 10−11 mV−1 and 2.4 × 10−19 m2V−2, respectively. However, ref. 49 did not provide a detailed explanation of the large THG signal compared to the SHG. Instead it assigned the large THG/SHG ratio to a possible enhancement of THG by the edge of the B exciton. However, refs. 51, 52 demonstrated that SHG is enhanced only when the SHG wavelength overlaps the A or B excitons. A similar behavior is expected for THG. Thus, the explanation in ref. 49 may not be correct. Reference 53 reported high-harmonic (>6th-order) generation in the non-perturbative regime with mid-infrared (IR) excitation (0.3 eV), unlike our THG and FHG results with near-IR excitation (0.8 eV). We do not detect THG from the thickest areas of our flake, with N > 30, as in ref. 18. The output spectrum in Fig. 3c further confirms that we observe both SHG and THG. Peaks for THG and SHG at 520 and 780 nm can be seen, as well as at 390 nm, corresponding to a four-photon process. This is detected only in 1L-MoS2. Its intensity is ~5.5 times lower than SHG, and two orders of magnitude smaller than THG.

SHG signals on areas with N = 3, 5, 7 have nearly the same intensity as 1L-MoS2, Fig. 4a. This contrasts ref. 14, where a pump laser at 810 nm was used. We attribute this difference to the fact that photons generated in the second-order nonlinear process in our set-up with a 1560 nm pump have an energy ~1.6 eV (780 nm), below the band gap of 1L-MoS2 12, therefore are not adsorbed, unlike the SHG signal in ref. 14.

Fig. 4
Experimental and theoretical nonlinear optical processes in MoS2. a SHG and THG intensities as functions of N. b Power dependence of SHG and THG in 1L-MoS2. c Experimental and theoretical THG/SHG irradiance ratio as a function of P pump. Different theoretical ...

Second- and third-order nonlinear susceptibilities

Based on the measured SHG and THG intensities, we can estimate the nonlinear susceptibilities χ (2) and χ (3). χ (2) can be calculated from the measured average powers of the fundamental and SH signals as follows54:


where τ is the pulse width, P pump is the average power of the incident fundamental (pump) beam, and P 2ω stands for the generated SH beam power, R is the repetition rate, N a = 0.5 is the numerical aperture, λ 2 = 780 nm is the SH wavelength, τ = τ 2 = 150 fs are the pulse durations at fundamental and SH wavelengths, ϕ=8π01cos-1ρ-ρ1-ρ22ρdρ=3.56 from ref. 54, and n1n2~1.45 are the refractive indexes of the substrate at the wavelengths of the fundamental and SHG, respectively. The effective bulk-like second-order susceptibility of MoS2 (χeff(2)) can be obtained from Eq. (1) with χeff(2)=χs(2)tMoS2, where tMoS2 = 0.65 nm is the 1L-MoS2 thickness10, 24. We obtain the effective second-order susceptibility χeff(2)~2.2 pmV−1 for 1L-MoS2. Reference 49 reported a bulk-like second-order susceptibility 29 pmV−1, which is ~10 times larger than here. However, several other studies reported ~5 pmV−1 for 1560 nm20, 52, 55. Thus, our measured χeff(2) agrees well with earlier values measured with similar excitation wavelength.

The third-order susceptibility χeff(3) can be estimated by comparing the measured THG signal from MoS2 to that of 1L-graphene (SLG):


with t SLG = 0.33 nm the SLG thickness, and THGSLG and THGMoS2 the measured signals from SLG and MoS2, respectively. Our results show that THG from 1L-MoS2 is around three times larger than THGSLG, which indicates that χ (3) of 1L-MoS2 is comparable to that of SLG, in the frequency range of our experiments. Previous reports indicate49, 56 that χSLG(3) is ~10−17–10−19 m2 V−2. Thus, based on Eq. (2), χ (3) of 1L-MoS2 is in the same range. This is remarkable, as SLG is known to have a large χ (3) 5659. Reference 18 reported χ (3) of 7L-MoS2 to be approximately three orders of magnitude smaller than χSLG(3) of ref. 56. χSLG(3) from ref. 56 is much higher than other theoretical58 and experimental49 values. We believe that our measured ratio between 1L-MoS2 and SLG is more accurate, since we measured both materials at the same time under the same conditions.

We note that large discrepancies can be found in earlier reported effective susceptibilities for layered materials (LM). For example, there is a approximately four orders of magnitude difference in χ (3) for SLG (~10−15 m2V−2 in ref. 57; ~10−19 m2V−2 in ref. 49). There is an approximately three orders of magnitude difference in χ (2) reported for 1L-MoS2 at 800 nm (e.g., ~10−7 mV−1 in ref. 15; and ~10−10 mV−1 in ref. 17). Effective susceptibilities are well defined only in three-dimensional materials, since their definition involves a polarization per unit volume1. Therefore, given the large discrepancies in literature, it is better to describe the nonlinear processes in LMs using the ratio between the harmonic signal power and the incident pump power (i.e., harmonic conversion efficiency). In this case, when comparing the efficiencies in our measurements with those in ref. 49, our THG efficiency (~4.76 × 10−10) is ~1.4 times larger than that (~3.38 × 10−10) in ref. 49, while our SHG efficiency (~6.47 × 10−11) is twice that of ref. 49. Since the effective susceptibilities are not well defined for LMs and also depend on the calculation method, we believe that the conversion efficiency is a better figure of merit for LMs.


Our measurements show that the nonlinear response of 1L-MoS2 and SLG are comparable in magnitude, both revealing stronger nonlinear efficiency than three-dimensional nonlinear materials, such as diamond1 and quartz59. This can be explained by considering their effective Hamiltonians27, 3437, 60. The main contribution to THG is paramagnetic. This is described by the square diagram in Supplementary Notes 14 (Supplementary Figs. 16). This paramagnetic contribution is mainly related to the strong inter-band coupling in the effective Hamiltonian, controlled by large velocity scales, vFc300 and v=t0a00.65×c300 for SLG and 1L-MoS2, with c the speed of light. The SLG paramagnetic third-harmonic efficiency (PTHE) is proportional to the square of third-order conductivity. Since39 σyyyy(3)vF2, we get an overall prefactor vF4, which explains the strong nonlinear SLG response. Similarly, for 1L-MoS2, the square diagram contains four paramagnetic current vertices, which gives an overall prefactor v 4, and an integral over the dummy momentum variables, which gives a prefactor 1v2 (see Supplementary Note 3). Therefore, the third-order response function, Πyyyy(3), is proportional to v 2, which implies a scaling of PTHE as v 4. Exciton physics is not considered because our experimental conditions only capture off-resonance transitions.

1L-MoS2 is transparent at this wavelength due to its ~1.9 eV gap12, while SLG absorbs 2.3% of the light61. Therefore, 1L-MoS2 and other TMDs are promising for integration with waveguides or fibers for all-optical nonlinear devices, such as all-optical modulators and signal processing devices, where materials with nonlinear properties are essential11.

The SHG and THG power dependence follows quadratic and cubic trends, Fig. 4b. At our power levels, THG is up to 30 times stronger than SHG. 1L-TMDs have strongly bound excitons that can modify their optical properties6264. The exciton resonances also affect their nonlinear optical responses17, 65, 66. References 51, 55 reported that when the SHG energy is above the A and B excitons, resonance effects are not observed. In our experiments, the energy of 3ω photons is above the A exciton but does not directly overlap with the A or B excitons. Thus, we do not assign the large THG/SHG intensity ratio to an excitonic enhancement, but to the approximate rotational invariance of the 1L-MoS2 band structure at low energies, which is broken by trigonal warping.

SHG is weaker than expected for a non-centrosymmetric material, due to near-isotropic bands contributing to the SHG signal for our low incident photon energies (0.8 eV). Even in the presence of a weak trigonal warping, SHG and THG might be comparable above the threshold for two- and three-photon absorption edges. However, this is not a resonant effect. Resonances only emerge when the laser matches a single level (like an excitonic level) rather than a continuum of states67. In our analysis, SHG would be absent without trigonal warping. But, trigonal warping alone cannot explain the magnitude of the FHG signal compared to SHG and THG.

Figure 4c compares the THG/SHG ratio from experiments and calculations based on the k·p theory35 (see Supplementary Note 1) and finite-temperature diagrammatic perturbation theory39 (see Supplementary Notes 3 and 4). The calculations are a factor 2 smaller than the experiments. Considering the complexity of the nonlinear optical processes and that our calculations ignore high-energy band structure effects29 and many-body renormalizations65, we believe this to be a satisfactory agreement, indicating the importance of trigonal warping in harmonic generation.

FHG generally derives from cascades of lower-order nonlinear multiphoton processes68. With an excitation wavelength of 1560 nm, this could be, e.g., a cascade of two SHG processes, where 780 nm photons are first generated through SHG (ω 1560 nm + ω 1560 nm [implies] ω 780 nm) and then undergo another SHG process (ω 780 nm + ω 780 nm [implies] ω 390 nm). To yield a FHG at 390 nm of the same intensity as SHG at 780 nm in this cascaded process, one would need a conversion efficiency (defined as P 2ω/P pump 1) for the second SHG process (i.e., ω 780 nm + ω 780 nm [implies] ω 390 nm) to be close to unity. However, we observe a conversion efficiency ~10−10 for SHG. Therefore, we conclude that our FHG does not arise from cascaded SHGs. Another possible cascade process is based on THG (ω 1560 nm + ω 1560 nm + ω 1560 nm [implies] ω 520nm) and sum-frequency generation (ω 520 nm + ω 1560 nm [implies] ω 390 nm). We find that THG strongly increases up to N = 5, as for Fig. 4a. Therefore, we expect this cascaded process to have a similar trend with N. However, we only observe FHG in 1L-MoS2. Thus, we also exclude this cascade process, and conclude that this is a direct χ (4) process.

We now consider the dependence of our results on the elliptical polarization of the incident light. We consider an incident laser beam with arbitrary polarization, i.e., E=Eε^± with ε^±=x^cos(θ)±iy^sin(θ). Using the crystal symmetries of 1L-MoS2, we derive (see Supplementary Note 2) the following expressions for the second- and third-order polarizations P (2) and P (3):




Note that θ = 0° corresponds to a linearly polarized laser along the x^ direction, perpendicular to the D3h1 mirror symmetry plane, while θ = 45° corresponds to a circularly polarized laser. From Eq. (3), we expect the intensity of SHG in response to a circularly polarized pump laser to be twice that of a linearly polarized laser. Equation (4) implies vanishing THG in response to a circularly polarized pump laser.

We measure the dependence of SHG and THG on elliptical polarization using a linearly polarized laser and a rotating QWP. Depending on the angle θ between the QWP axes and the laser polarization, the excitation light is linearly (θ =  + m·90°) or circularly (θ = 45° + m·90°) polarized. Figure 5 shows that the experimental data are in agreement with Eqs. (3) and (4). The THG signal is maximum for a linearly polarized excitation laser, while it vanishes for circularly polarized light. SHG is always visible, but its intensity is maximum for circularly polarized light.

Fig. 5
SHG and THG dependence on the pump light polarization. The polar plot angle corresponds to linearly polarized light when θ =  + m·90°, and to circularly polarized pump light when θ ...

Given that harmonic generation is strongly dependent on the symmetry and stacking of layers and that different 1L-TMDs (e.g., WSe2, MoSe2) all have similar nonlinear response11, 14, 15, 21, one could use heterostructures (e.g., MoS2/WSe2) to engineer SHG and other nonlinear processes for high photon-conversion efficiency for a wide range of applications requiring the generation of higher frequencies. This may lead to the use of LMs and heterostructures for applications utilizing optical nonlinearities (e.g., all-optical devices, frequency combs, high-order harmonic generation, multiphoton microscopy, and therapy etc.).


Determination of MoS2 thickness from SHG and THG signals

SHG and THG for FL-MoS2 (N = 1…7) are studied on the flakes in Fig. 6a. SHG and THG images are shown in Fig. 6b, c. At 1560 nm, the contrast between 1 and 3L areas is small, as well as the contrast between 3, 5, and 7 L regions (Fig. 6b).

Fig. 6
Optical and multiphoton images of few-layer MoS2 flake. a Optical micrograph, b SHG, and c THG images of flake with few-layer areas under 1560 nm excitation

The THG signal increases up to N = 7, Figs. 4a and and6c.6c. On the other hand, the SHG signal (Fig. 6b) is only generated in ON flakes, due to symmetry14. Therefore, areas with intensity between the 3, 5, and 7L regions in Fig. 6c, but dark in SHG, are 4 and 6 L. The dependence of the intensities of THG and SHG on N is plotted in Fig. 4a. The combination of SHG and THG can be used to identify N at least up to 7. The THG signal develops as a function of N. Using Maxwell’s equations for a nonlinear medium with thickness t and considering the slowly varying amplitude approximation1, 69, we obtain:


where I in and I 3ω are the intensities of the incident and THG light, respectively, and χ (3)(−3ω;ω, ω, ω) is the third-order optical susceptibility, nj=1,3=ϵ(1)(jω), with ϵ(1) the TMD linear dielectric function. Δkt is the phase mismatch between the fundamental and third harmonic generated waves.

For Δkt  0, THG adds up quadratically with light propagation length (i.e., t [proportional, variant] N). The signal starts to saturate for N = 6. The possible reasons for subquadratic signal build-up can be either phase mismatch, or absorption13. For THG, Δk = 3k in ± k 3ω, where k in and k 3ω are the wavevectors of the incident and THG signals, respectively, where the plus sign indicates THG generated in the backward direction, while minus identifies forward generated THG. Even for backward generated THG, Δkt  0 for 6L-MoS2 (~4.3 nm70). This rules out phase mismatch as the origin of the signal saturation when N  6. Therefore, we assume that the signal saturation is due to absorption of the third harmonic light.

Diagrammatic nonlinear response theory

To quantify theoretically the strength of nonlinear harmonic generation processes, we generalize the diagrammatic perturbation theory approach39 to the case of TMDs. We combine this technique with a low-energy k·p model Hamiltonian ℋ(k) for 1L-MoS2 35. In such low-energy model, light-matter interactions are treated by employing minimal coupling35, 39, kk + e A(t)/ħ, where A(t) is a time-dependent uniform vector potential. Nonlinear response functions are calculated via the multi-legged Feynman diagrams depicted in Supplementary Figs. 2 and 3.

Data availability

The data that support the findings of this study are available from the corresponding author on request.

Electronic supplementary material


We thank M.J. Huttunen and R.W. Boyd for useful discussions. We acknowledge funding from the Academy of Finland (Nos: 276376, 284548, 295777, 298297, and 304666), TEKES (NP-Nano, OPEC), Royal Academy of Engineering (RAEng) Research Fellowships, Fondazione Istituto Italiano di Tecnologia, the Graphene Flagship, ERC grants Hetero2D, Nokia Foundation, EPSRC Grants EP/K01711X/1, EP/K017144/1, EP/L016087/1, AFOSR COMAS MURI (FA9550-10-1-0558), ONR NECom MURI, CIAN NSF ERC under Grant EEC-0812072, and TRIF Photonics funding from the state of Arizona and the Micronova, Nanofabrication Centre of Aalto University.

Author contributions

Author contributions

A.S., L.K., A.A., S.M., K.K., and Z.S. carried out the multiphoton experiments and data analysis. Samples preparation and characterization were carried out by A.L. and A.C.F. The nonlinear response function calculations were carried out by H.R. and M.P. Sample fabrication, multiphoton characterization, and analysis were coordinated by R.A.N., N.P., H.L., A.C.F., M.P., and Z.S. All authors contributed to the manuscript.


Competing interests

The authors declare no competing financial interests.


Antti Säynätjoki, Lasse Karvonen and Habib Rostami contributed equally to this work.

Electronic supplementary material

Supplementary Information accompanies this paper at 10.1038/s41467-017-00749-4.

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


1. Boyd, R. W. Nonlinear Optics (Academic Press, 2003).
2. Pavone, F. S. & Campagnola, P. J. Second Harmonic Generation Imaging (Taylor & Francis, 2013).
3. Saleh, B. E. A. & Teich, M. C. Fundamentals of Photonics (Academic Press, 2003).
4. Willner AE, Khaleghi S, Chitgarha MR, Yilmaz OF. All-optical signal processing. J. Lightwave Technol. 2014;32:660–680. doi: 10.1109/JLT.2013.2287219. [Cross Ref]
5. Zipfel WR, Williams RM, Webb WW. Nonlinear magic: multiphoton microscopy in the biosciences. Nat. Biotechnol. 2003;21:1369–1377. doi: 10.1038/nbt899. [PubMed] [Cross Ref]
6. Bhawalkar JD, He GS, Prasad PN. Nonlinear multiphoton processes in organic and polymeric materials. Rep. Prog. Phys. 1996;59:1041. doi: 10.1088/0034-4885/59/9/001. [Cross Ref]
7. Tsang TYF. Optical third-harmonic generation at interfaces. Phys. Rev. A. 1995;52:4116. doi: 10.1103/PhysRevA.52.4116. [PubMed] [Cross Ref]
8. Karvonen L, et al. Investigation of second-and third-harmonic generation in few-layer gallium selenide by multiphoton microscopy. Sci. Rep. 2015;5:10334. doi: 10.1038/srep10334. [PMC free article] [PubMed] [Cross Ref]
9. Bonaccorso F, Sun Z, Hasan T, Ferrari AC. Graphene photonics and optoelectronics. Nat. Photon. 2010;4:611–622. doi: 10.1038/nphoton.2010.186. [Cross Ref]
10. Ferrari AC, et al. Science and technology roadmap for graphene, related two-dimensional crystals, and hybrid systems. Nanoscale. 2015;7:4598–4810. doi: 10.1039/C4NR01600A. [PubMed] [Cross Ref]
11. Sun Z, Martinez A, Wang F. Optical modulators with 2D layered materials. Nat. Photon. 2016;10:227–238. doi: 10.1038/nphoton.2016.15. [Cross Ref]
12. Mak KF, Lee C, Hone J, Shan J, Heinz TF. Atomically thin MoS2: a new direct-gap semiconductor. Phys. Rev. Lett. 2010;105:136805. doi: 10.1103/PhysRevLett.105.136805. [PubMed] [Cross Ref]
13. Eda G, et al. Photoluminescence from chemically exfoliated MoS2. Nano Lett. 2011;11:5111–5116. doi: 10.1021/nl201874w. [PubMed] [Cross Ref]
14. Li Y, et al. Probing symmetry properties of few-layer MoS2 and h-BN by optical second-harmonic generation. Nano Lett. 2013;13:3329–3333. doi: 10.1021/nl401561r. [PubMed] [Cross Ref]
15. Kumar N, et al. Second harmonic microscopy of monolayer MoS2. Phys. Rev. B. 2013;87:161403. doi: 10.1103/PhysRevB.87.161403. [Cross Ref]
16. Wang K, et al. Ultrafast saturable absorption of two-dimensional MoS2 nanosheets. ACS Nano. 2013;7:9260–9267. doi: 10.1021/nn403886t. [PubMed] [Cross Ref]
17. Malard LM, Alencar TV, Barboza APM, Mak KF, de Paula AM. Observation of intense second harmonic generation from MoS2 atomic crystals. Phys. Rev. B. 2013;87:201401. doi: 10.1103/PhysRevB.87.201401. [Cross Ref]
18. Wang R, et al. Third-harmonic generation in ultrathin films of MoS2. ACS Appl. Mater. Interfaces. 2013;6:314–318. doi: 10.1021/am4042542. [PubMed] [Cross Ref]
19. Trolle ML, Seifert G, Pedersen TG. Theory of excitonic second-harmonic generation in monolayer MoS2. Phys. Rev. B. 2014;89:235410. doi: 10.1103/PhysRevB.89.235410. [Cross Ref]
20. Clark D, et al. Strong optical nonlinearity of CVD-grown MoS2 monolayer as probed by wavelength-dependent second-harmonic generation. Phys. Rev. B. 2014;90:121409. doi: 10.1103/PhysRevB.90.121409. [Cross Ref]
21. Seyler KL, et al. Electrical control of second-harmonic generation in a WSe2 monolayer transistor. Nat. Nanotechnol. 2015;10:407–411. doi: 10.1038/nnano.2015.73. [PubMed] [Cross Ref]
22. Kuc A, Zibouche N, Heine T. Influence of quantum confinement on the electronic structure of the transition metal sulfide TS2. Phys. Rev. B. 2011;83:245213. doi: 10.1103/PhysRevB.83.245213. [Cross Ref]
23. Kadantsev ES, Hawrylak P. Electronic structure of a single MoS2 monolayer. Solid State Commun. 2012;152:909–913. doi: 10.1016/j.ssc.2012.02.005. [Cross Ref]
24. Wang QH, Kalantar-Zadeh K, Kis A, Coleman JN, Strano MS. Electronics and optoelectronics of two-dimensional transition metal dichalcogenides. Nat. Nanotechnol. 2012;7:699–712. doi: 10.1038/nnano.2012.193. [PubMed] [Cross Ref]
25. Shi H, Pan H, Zhang Y-W, Yakobson BI. Quasiparticle band structures and optical properties of strained monolayer MoS2 and WS2. Phys. Rev. B. 2013;87:155304. doi: 10.1103/PhysRevB.87.155304. [Cross Ref]
26. Zahid F, Liu L, Zhu Y, Wang J, Guo H. A generic tight-binding model for monolayer, bilayer and bulk MoS2. AIP Adv. 2013;3:052111. doi: 10.1063/1.4804936. [Cross Ref]
27. Kormányos A, et al. Monolayer MoS2: trigonal warping, the γ valley, and spin-orbit coupling effects. Phys. Rev. B. 2013;88:045416. doi: 10.1103/PhysRevB.88.045416. [Cross Ref]
28. Qiu DY, Felipe H, Louie SG. Optical spectrum of MoS2: many-body effects and diversity of exciton states. Phys. Rev. Lett. 2013;111:216805. doi: 10.1103/PhysRevLett.111.216805. [PubMed] [Cross Ref]
29. Gibertini M, Pellegrino FM, Marzari N, Polini M. Spin-resolved optical conductivity of two-dimensional group-VIb transition-metal dichalcogenides. Phys. Rev. B. 2014;90:245411. doi: 10.1103/PhysRevB.90.245411. [Cross Ref]
30. Margulis VA, Muryumin EE, Gaiduk EA. Optical second-harmonic generation from two-dimensional hexagonal crystals with broken space inversion symmetry. J. Phys. Condens. Matter. 2013;25:195302. doi: 10.1088/0953-8984/25/19/195302. [PubMed] [Cross Ref]
31. Wu S, et al. Quantum-enhanced tunable second-order optical nonlinearity in bilayer graphene. Nano Lett. 2012;12:2032–2036. doi: 10.1021/nl300084j. [PubMed] [Cross Ref]
32. Brun SJ, Pedersen TG. Intense and tunable second-harmonic generation in biased bilayer graphene. Phys. Rev. B. 2015;91:205405. doi: 10.1103/PhysRevB.91.205405. [Cross Ref]
33. Hipolito F, Pedersen TG, Pereira VM. Nonlinear photocurrents in two-dimensional systems based on graphene and boron nitride. Phys. Rev. B. 2016;94:045434. doi: 10.1103/PhysRevB.94.045434. [Cross Ref]
34. Xiao D, Liu G-B, Feng W, Xu X, Yao W. Coupled spin and valley physics in monolayers of MoX2 and other group-vi dichalcogenides. Phys. Rev. Lett. 2012;108:196802. doi: 10.1103/PhysRevLett.108.196802. [PubMed] [Cross Ref]
35. Rostami H, Roldán R, Cappelluti E, Asgari R, Guinea F. Theory of strain in single-layer transition metal dichalcogenides. Phys. Rev. B. 2015;92:195402. doi: 10.1103/PhysRevB.92.195402. [Cross Ref]
36. Rostami H, Asgari R, Guinea F. Edge modes in zigzag and armchair ribbons of monolayer MoS2. J. Phys. Condens. Matter. 2016;28:495001. doi: 10.1088/0953-8984/28/49/495001. [PubMed] [Cross Ref]
37. Rostami H, Moghaddam AG, Asgari R. Effective lattice hamiltonian for monolayer MoS2: tailoring electronic structure with perpendicular electric and magnetic fields. Phys. Rev. B. 2013;88:085440. doi: 10.1103/PhysRevB.88.085440. [Cross Ref]
38. Alidoust N, et al. Observation of monolayer valence band spin-orbit effect and induced quantum well states in MoX2. Nat. Commun. 2014;5:4673. doi: 10.1038/ncomms5673. [PubMed] [Cross Ref]
39. Rostami H, Polini M. Theory of third-harmonic generation in graphene: a diagrammatic approach. Phys. Rev. B. 2016;93:161411. doi: 10.1103/PhysRevB.93.161411. [Cross Ref]
40. Bonaccorso F, et al. Production and processing of graphene and 2d crystals. Mater. Today. 2012;15:564–589. doi: 10.1016/S1369-7021(13)70014-2. [Cross Ref]
41. Sundaram, R. et al. Electroluminescence in single layer MoS2. Nano Lett. 13, 1416–1421 (2013). [PubMed]
42. Casiraghi C, et al. Rayleigh imaging of graphene and graphene layers. Nano Lett. 2007;7:2711–7. doi: 10.1021/nl071168m. [PubMed] [Cross Ref]
43. Zhang X, et al. Raman spectroscopy of shear and layer breathing modes in multilayer MoS2. Phys. Rev. B. 2013;87:115413. doi: 10.1103/PhysRevB.87.115413. [Cross Ref]
44. Säynätjoki A, et al. Rapid large-area multiphoton microscopy for characterization of graphene. ACS Nano. 2013;7:8441–8446. doi: 10.1021/nn4042909. [PubMed] [Cross Ref]
45. Kieu K, Jones RJ, Peyghambarian N. Generation of few-cycle pulses from an amplified carbon nanotube mode-locked fiber laser system. IEEE Photon. Technol. Lett. 2010;22:1521–1523. doi: 10.1109/LPT.2010.2063423. [Cross Ref]
46. Kieu K, Jones RJ, Peyghambarian N. High power femtosecond source near 1 micron based on an all-fiber er-doped mode-locked laser. Opt. Express. 2010;18:21350–21355. doi: 10.1364/OE.18.021350. [PubMed] [Cross Ref]
47. Miller RC. Optical second harmonic generation in piezoelectric crystals. Appl. Phys. Lett. 1964;5:17–19. doi: 10.1063/1.1754022. [Cross Ref]
48. Chin AH, Calderón OG, Kono J. Extreme midinfrared nonlinear optics in semiconductors. Phys. Rev. Lett. 2001;86:3292. doi: 10.1103/PhysRevLett.86.3292. [PubMed] [Cross Ref]
49. Woodward RI, et al. Characterization of the second- and third-order nonlinear optical susceptibilities of monolayer MoS2 using multiphoton microscopy. 2D Mater. 2017;4:011006. doi: 10.1088/2053-1583/4/1/011006. [Cross Ref]
50. Säynätjoki, A. et al. Ultra-strong nonlinear optical processes and trigonal warping in MoS2 layers. Preprint at (2016). [PMC free article] [PubMed]
51. Clark D, et al. Near bandgap second-order nonlinear optical characteristics of MoS2 monolayer transferred on transparent substrates. Appl. Phys. Lett. 2015;107:131113. doi: 10.1063/1.4932134. [Cross Ref]
52. Le CT, et al. Impact of selenium doping on resonant second harmonic generation in monolayer MoS2. ACS Photon. 2017;4:38–44. doi: 10.1021/acsphotonics.6b00530. [Cross Ref]
53. Liu H, et al. High-harmonic generation from an atomically thin semiconductor. Nat. Phys. 2017;13:262–265. doi: 10.1038/nphys3946. [Cross Ref]
54. Janisch, C. et al. Extraordinary second harmonic generation in tungsten disulfide monolayers. Sci. Rep. 4, 5530 (2014). [PMC free article] [PubMed]
55. Le CT, et al. Nonlinear optical characteristics of monolayer MoSe2. Ann. Phys. 2016;528:551–559. doi: 10.1002/andp.201600006. [Cross Ref]
56. Kumar N, et al. Third harmonic generation in graphene and few-layer graphite films. Phys. Rev. B. 2013;87:121406. doi: 10.1103/PhysRevB.87.121406. [Cross Ref]
57. Hendry E, Hale PJ, Moger J, Savchenko A, Mikhailov S. Coherent nonlinear optical response of graphene. Phys. Rev. Lett. 2010;105:097401. doi: 10.1103/PhysRevLett.105.097401. [PubMed] [Cross Ref]
58. Cheng J, Vermeulen N, Sipe J. Third order optical nonlinearity of graphene. New J. Phys. 2014;16:053014. doi: 10.1088/1367-2630/16/5/053014. [Cross Ref]
59. Hong S-Y, et al. Optical third-harmonic generation in graphene. Phys. Rev. X. 2013;3:021014.
60. Neto AHC, Guinea F, Peres NM, Novoselov KS, Geim AK. The electronic properties of graphene. Rev. Mod. Phys. 2009;81:109. doi: 10.1103/RevModPhys.81.109. [Cross Ref]
61. Nair RR, et al. Fine structure constant defines visual transparency of graphene. Science. 2008;320:1308–1308. doi: 10.1126/science.1156965. [PubMed] [Cross Ref]
62. Ramasubramaniam A. Large excitonic effects in monolayers of molybdenum and tungsten dichalcogenides. Phys. Rev. B. 2012;86:115409. doi: 10.1103/PhysRevB.86.115409. [Cross Ref]
63. Cheiwchanchamnangij T, Lambrecht WRL. Quasiparticle band structure calculation of monolayer, bilayer, and bulk MoS2. Phys. Rev. B. 2012;85:205302. doi: 10.1103/PhysRevB.85.205302. [Cross Ref]
64. Mak KF, et al. Tightly bound trions in monolayer MoS2. Nat. Mater. 2013;12:207–211. doi: 10.1038/nmat3505. [PubMed] [Cross Ref]
65. Grüning M, Attaccalite C. Second harmonic generation in h-BN and MoS2 monolayers: role of electron-hole interaction. Phys. Rev. B. 2014;89:081102. doi: 10.1103/PhysRevB.89.081102. [Cross Ref]
66. Wang G, et al. Giant enhancement of the optical second-harmonic emission of WSe2 monolayers by laser excitation at exciton resonances. Phys. Rev. Lett. 2015;114:097403. doi: 10.1103/PhysRevLett.114.097403. [PubMed] [Cross Ref]
67. Haug, H. & Koch, S. W. Quantum Theory of the Optical and Electronic Properties of Semiconductors (World Scientific Publishing, 2009).
68. Zhu S-n, Zhu Y-y, Ming N-b. Quasi-phase-matched third-harmonic generation in a quasi-periodic optical superlattice. Science. 1997;278:843–846. doi: 10.1126/science.278.5339.843. [Cross Ref]
69. Butcher, P. N. & Cotter, D. The Elements of Nonlinear Optics (Cambrige University Press, 1990).
70. Radisavljevic B, Radenovic A, Brivio J, Giacometti IV, Kis A. Single-layer MoS2 transistors. Nat. Nanotechnol. 2011;6:147–150. doi: 10.1038/nnano.2010.279. [PubMed] [Cross Ref]

Articles from Nature Communications are provided here courtesy of Nature Publishing Group