{'fraction_non_alphanumeric': 0.03148902354226681, 'fraction_numerical': 0.020595344525015727, 'mean_word_length': 4.9277723258096175, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 1, 'lorem ipsum': 0, 'www.': 6, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Polarization-dependent optical responses in natural 2D layered mineral teallite

Multi-element layered materials enable the use of stoichiometric variation to engineer their optical responses at subwavelength scale. In this regard, naturally occurring van der Waals minerals allow us to harness a wide range of chemical compositions, crystal structures and lattice symmetries for layered materials under atomically thin limit. Recently, one type of naturally occurring sulfide mineral, ternary teallite has attained significant interest in the context of thermoelectric, optoelectronic, and photovoltaic applications, but understanding of light-matter interactions in such ternary teallite crystals is scarcely available. Herein, polarization-dependent linear and nonlinear optical responses in mechanically exfoliated teallite crystals are investigated including anisotropic Raman modes, wavelength-dependent linear dichroism, optical band gap evolution, and anisotropic third-harmonic generation (THG). Furthermore, the third-order nonlinear susceptibility of teallite crystal is estimated using the thickness-dependent THG emission process. We anticipate that our findings will open the avenue to a better understanding of the tailored light-matter interactions in complex multi-element layered materials and their implications in optical sensors, frequency modulators, integrated photonic circuits, and other nonlinear signal processing applications.


Scientific Reports
| (2021) 11:21895 | https://doi.org/10.1038/s41598-021-01511-z www.nature.com/scientificreports/ fundamental understanding of linear and nonlinear optical responses of teallite crystal yet. Moreover, these insights regarding the crystal structures and other microscopic processes are desirable for performance efficiency enhancement and engineering the application-based optical and optoelectronic responses for the improvisation of nanoscale devices. Motivated with this, herein we explore the polarization-dependent anisotropic linear and nonlinear optical responses of mechanically exfoliated ternary teallite thin crystals. Thin teallite flakes are prepared via the mechanical exfoliation from natural bulk teallite mineral. The prepared teallite flakes are characterized using optical microscope, atomic force microscope (AFM), high-resolution transmission electron microscopy (HRTEM) and energy dispersive X-ray spectroscopy (EDXS) techniques to determine the crystal's information, including flake shape, flake size, layer thickness, surface smoothness, atomic arrangement and chemical composition. Further, the flakes are probed spectroscopically to determine the crystalline axis and optical axis with polarization-resolved Raman spectroscopy and optical absorption spectroscopy. The optical absorption spectra are further analyzed to understand the angle-resolved evolution of optical band gap in teallite semiconductor. In last, anisotropic THG response of teallite crystal with respect to the incident linear polarization of pump beam is investigated and we estimate the third-order nonlinear susceptibility from the thickness-dependent THG emission process. The experimental results are corroborated by the theoretical models. Our results will lead to new insights in better understanding of light-matter interactions in complex multi-element vdW layered materials and advancing many photonics applications such as frequency modulators, integrated photonic circuits, and nonlinear optical signal processing.

Results
Sample preparation and characterization of teallite crystals. Figure 1a shows the simplified view of the teallite crystal structure. As presented in Fig. 1a, teallite forms a puckered bilayer structure similar to that of black phosphorus and herzenbergite SnS with S atoms bonded to Pb or Sn atoms, where the second layer is rotated relative to the first layer by 180° according to the AB stacking order in the unit cell bilayer. It belongs to the orthorhombic crystal system with space group Pnma (No. 62) with the unit cell dimensions of a = 4.09 Å, b = 4.26 Å, c = 11.41 Å and α = β = γ = 90°. The a and b axes correspond to the zigzag and armchair directions, with the c-axis perpendicular to the a-b plane and α , β and γ are the angles between a-, b-and c-axis. Figure 1b shows a picture of a teallite mineral rock, where silvery-gray tightly packed micaceous foliated plates of teallite crystals are presented. The magnified view of a section is shown in Fig. 1c. We prepare the teallite thin flakes of various thicknesses on glass substrate via the mechanical exfoliation of the bulk natural teallite mineral (Monserrat, Oruro, Bolivia) using Nitto tape (SPV 224) and Scotch tape. Figure 1d,e,g,h,j,k,m,n show the optical microscope reflection and transmission images of several prepared teallite flakes with the thicknesses of 9, 24, 36 and 61 nm, respectively. Further, these teallite flakes are scanned using AFM to estimate the thickness and surface roughness. The captured AFM images are shown in Fig. 1f (9 nm), 1i (24 nm), 1l (36 nm), and 1o (61 nm), respectively.
Next, the structural information and chemical composition of teallite crystal are determined using HRTEM and EDXS characterization techniques. Figure 2a shows a HRTEM image of teallite crystal with the determined lattice spacing of 2.02 Å and 2.14 Å along the zigzag direction (a-axis) and the armchair direction (b-axis), with the intersection angle of around 90°. These spacing are found consistent with the (200) and (020) sets of planes for the orthorhombic teallite crystal. Both (110) and ( 110) sets of planes have the lattice spacings of 2.94 Å with the intersection angle of around 86°, and both sets of planes have the angle of around 43° with respect to the (020) planes. Further, the recorded selected area electron diffraction (SAED) pattern is displayed in Fig. 2b, where the spot patterns are captured from the surface normal to the [001] crystal zone axis. The observed main reflection pattern is consistent with the previously reported pure teallite crystal structure 32 . This consistency underlines that the investigated crystal is of high quality and the crystalline nature or purity of the investigated crystal is not affected. Besides, we also observe few weak-intensity reflections in addition to the main reflections. We ascribe that these weak-intensity reflections are due to the slight distortion and modulation in actual teallite structure. Such minor distortion and modulation of the lattice parameters, unit cell dimensions and lattice symmetries in the resultant teallite crystal can be expected due to the geological origin of the naturally grown teallite mineral and the dual phase coexistence of herzenbergite-teallite. Subsequently, the chemical composition of ternary teallite crystal is quantified. Figure 2c shows the collected average EDXS spectrum, while the bright-field TEM image of the scanned region and the corresponding individual elemental mappings are shown in Fig. 2d-g. These recorded mappings emphasize the homogenous distribution of prime constituent elements of lead (Pb), tin (Sn) and sulfur (S). The obtained average elemental composition (as per atomic weight %) is summarized in Table 1, which is used to perform the compositional stoichiometry analysis. The resultant empirical formula is found to be Pb 0.58 Sn 1.58 S 2.00 with the Pb:Sn ratio of 0.367. The obtained formula is not the excellent match with the generic teallite chemical formula PbSnS 2 and the range of the chemical composition of constituent elements are significantly broad. These two factors emphasize the natural occurrence of herzenbergite-teallite series in the probed specimen 37 . It is noteworthy that the obtained formula lies within the variation range in chemical composition between Pb 0.428 Sn 1.572 S 2 and Pb 0.784 Sn 1.216 S 2 according to the chemical analysis of a series of herzenbergite-teallite minerals from Bolivia 38 , where the Pb:Sn ratio is from 0.272 to 0.645. The appearance of nickel (Ni) peak in the EDXS spectrum is ascribed to the TEM grid. Besides, the presence of small amounts of additional elements such as copper (Cu), carbon (C) and oxygen (O) in the recorded EDXS spectrum can be attributed to the mineral geological origins 25,27,37 . In addition, the high contribution of carbon (C) and oxygen (O) can be ascribed to the sample processing for TEM analysis and the underlying carbon support film of TEM grid. We used polymethyl methacrylate (PMMA)-assisted wet transfer method for transferring the mechanically exfoliated teallite flakes from glass substrate to TEM grid, which further leads to the accumulation of carbon-based adsorbates during  39 . Additionally, the recorded Raman modes of the teallite mineral sample are also compared with the Raman modes of the previously reported chemically synthesized pure teallite crystal to further elucidate any significant influence of these doped elements on the intrinsic and optical properties of teallite mineral. We found a good agreement between the observed Raman modes of the teallite mineral and those of the chemically synthesized highly pure teallite crystal. Thus, we conclude that the probed teallite mineral crystal is substantially pure, and the influence of the geological origin induced impurities on the optical responses would be minimal.
Characterization of anisotropic Raman modes and determination of crystal axis in teallite crystals. Next, the teallite flakes are characterized with polarization-resolved Raman spectroscopy to identify the crystal' axes. The flake is illuminated with a 632.8 nm He-Ne laser. The desired input linear polarization is obtained by using a combination of a linear polarizer and a half-wave plate in the excitation path, while an analyzer is engaged in the collection path set to be in the parallel direction with respect to the input beam polarization. Figure 3a shows the recorded Raman spectrum from the 61 nm-thick teallite crystal in parallel polarization configuration, showing a series of intense Raman vibrational modes at 93, 139, 182 and 218 cm −1 . Importantly, there are some similarities between the Raman spectra of teallite and herzenbergite SnS 40,41 . The Raman peaks of 93, 182 and 218 cm −1 are assigned to the A g modes, which are originated from the waving mode, the breathing mode, and the NaCl-type vibrational mode, respectively. While the peak of 139 cm −1 corresponds to the B 3g mode, which is attributed to the NaCl-type vibration along the zigzag direction. To gain further insight of the angular periodicity of these Raman modes, polarization-dependent Raman spectra are studied. Figure 3b shows the color map of the Raman intensity spectrum as per the input linear polarization angle in parallel polarization configuration, where the anisotropic nature of these Raman vibrational modes can be clearly observed. Figure 3c shows the optical microscope reflection image of the 61 nm-thick flake, and the corresponding polar plots of where // denotes the parallel polarization configuration, a, b and e are the amplitudes of Raman tensor elements, and ∅ ba = ∅ b − ∅ a denotes the associated phase difference between a and b elements. Further, we ensure the consistency of the observed angle-resolved Raman modes for teallite crystals with different thicknesses. Figure 3g,k,o show the optical reflection images of the investigated 9, 24 and 36 nm-thick flakes, and the corresponding A g and B 3g modes are also plotted in Fig. 3h-j,l-n,p-r, respectively. The experimental data are denoted as black squares, whereas the theoretical fitting curves are shown with magenta solid lines, indicating a good agreement with each other. It is found that all the A g modes show anisotropic two-lobe patterns with the maxima at 0° and 180°, whereas the B 3g mode shows anisotropic four-lobe pattern with the maxima at 45°, 135°, 225° and 315°. Furthermore, the orientation of A g mode pattern indicates the crystal axis of the probed teallite flake, where θ = 0° direction (along x-axis) and θ = 90° direction (along y-axis) are accredited as the b-axis (armchair direction) and a-axis (zigzag direction) of teallite crystal, respectively.

Polarization-resolved optical absorption and band gap evolution in teallite crystals.
The reduced in-plane lattice symmetry and the prevalent anisotropic Raman vibrational modes of teallite crystal also endorse the presence of linear dichroism. To get further insight, we probe the optical absorption characteristics of the 61 nm flake within the visible frequency regime from 420 to 750 nm using polarization-resolved optical absorption spectroscopy. Figure 4a shows the measured reflectance (R), transmittance (T) and absorbance (A = 1 -R -T) spectra, keeping the input linear polarization fixed along the optical axis of the crystal in the [ 110] lattice   www.nature.com/scientificreports/ direction. Both the reflectance and transmittance spectra show an initial resonance dip around the wavelength of 450 nm and then gradually increase within the range of 450 nm to 750 nm. Also, the reflectance is higher as compared to the transmittance all over the visible frequency range. Accordingly, the absorbance spectrum shows a peak around 450 nm and then decreases up to the wavelength of 750 nm. To further examine the anisotropic features of optical absorption and linear dichroism, angle-resolved absorbance spectra are systematically measured for the 61 nm-thick teallite crystal, as shown in Fig. 4b. The input linear polarization angle is set relative to the b-axis of the crystal. Figure 4c,d show the polar plots of absorbance as per incident linear polarization angle at two different wavelengths of 520 nm and 720 nm, exhibiting anisotropic two-lobe patterns with the maxima at around − 43° and 137°. It is inferred that the teallite crystal anisotropically absorbs photons with the preference along the optical axis in the [ 110] lattice direction, which is around 43° from the armchair direction along the b-axis. The optical axis of teallite crystal is indicated as the x′-axis in Fig. 5a. The experimentally measured data points are theoretically fitted using the equation A(θ) = A x ′ cos 2 (θ) + A y ′ sin 2 (θ) , where A x ′ and A y ′ are the absorbance magnitudes along the optical axis and its perpendicular direction. The absorption anisotropic ratio A x ′ /A y ′ is obtained as 1.08 and 1.29 at 520 nm and 720 nm, respectively, showing that the teallite crystal exhibits strong wavelength-dependent linear dichroism effect. We understand such lattice direction preferentiality for optical absorption in the context of the distorted lattice symmetry due to the structural complexity of teallite crystal compared to the standard herzenbergite SnS crystal. It is still unknown about how the Pb and Sn elements are exactly distributed in ideal PbSnS 2 crystal yet. It has been studied that the Pb element alloying into SnS lattice can induce the sizable structural deformation of the in-plane lattice for PbSnS 2 crystal, where the inclusion of 20% Pb into SnS leads to the lattice parameter a along the zigzag direction increases by 2% and the lattice parameter b along the armchair direction decreases by 1% 31 . Furthermore, there are four distinct structural phases commonly occurring in PbSnS 2 crystal with different layer stacking sequences along the c-axis and different distances between the successive layers. All these structural phases for describing the PbSnS 2 crystal are orthorhombic with the same unit cell dimensions, but exhibit differing lattice symmetry due to different Pb and Sn atom arrangement and ordering 43 . Therefore, the sizable structural deformation of the in-plane lattice and the lattice symmetry distortion of teallite crystal induced by the Pb alloying into SnS will affect the optical responses. In addition, the dual phase coexistence of herzenbergite-teallite in the natural mineral will further play a crucial role in modulating the lattice parameters, unit cell dimensions and the associated lattice symmetry in the resultant teallite crystal 37 . As shown in the EDXS analysis, the Pb:Sn ratio in our sample varies from 0.321 to 0.421. Such broad range variations in Pb:Sn ratio validate the coexistence of the herzenbergite-teallite phase and the imperfect miscibility in the natural mineral 37 . As www.nature.com/scientificreports/ a result, the investigated natural teallite crystals have more complicated distorted unit cell dimensions and lattice symmetries in contrast to the ideal PbSnS 2 crystal structure, which will further affect the optical responses and may lead to the maximum optical absorption along the optical axis in the [ 110] lattice direction, rather than the armchair direction along the b-axis. Nevertheless, further insight of these optical processes in stoichiometric controlled synthesis of complex 2D materials are inevitable to comprehensively explore such intriguing phenomena. Moreover, such stoichiometric variation in these multi-element 2D materials could also strongly modulate the electronic band structure and optical band gap of the material. Hence, as a ternary-element semiconductor, teallite is an interesting specimen to explore the polarization-dependent evolution of optical band gap. The optical band gap of teallite crystal is estimated by the Tauc plot using the relation (αhν) n = C(hν − E g ) 34 , where E g is optical band gap energy, α is the absorption coefficient corresponding to the photon energy hν ( h is the Plank's constant and ν is the frequency of the incident photon), C is a constant and n = 2 , 1/2 , 2/3 indicate an allowed direct transition, an allowed indirect transition, and a forbidden transition, respectively. The value of α can be extracted from the measured reflectance and transmittance spectra shown in Fig. 4a with the flake thickness d using the following equation 44 Figure 4e shows the angle-resolved plot of (αhν) 2 as a function of the photon energy hν for the 61 nm-thick teallite crystal. The analysis of the Tauc plots in between 2.201 eV and 2.594 eV energy region gives a good linear fit for n = 2 , which underlines the allowed direct energy transition and the presence of direct optical band gap www.nature.com/scientificreports/ close to 1.87 eV. This value shows a decent match with the previously reported band gap for PbSnS 2 thin films 32,34 . Interestingly, these extracted values of direct band gap from the Tauc plots show a systematic variation with the incident linear polarization angle. Figure 4f show the polar plot of optical band gap evolution with respect to the incident linear polarization angle. It is observed that the optical band gap of teallite crystal exhibits an anisotropic two-lobe pattern with the maximum value of 1.92 eV along the optical axis in the [ 110] lattice direction and the minimum value of 1.83 eV along its perpendicular direction, which is strongly influenced with the structural symmetry in such complex 2D materials.
Anisotropic THG response and third-order nonlinear susceptibility in teallite crystals. In last, we explore the effect of the low in-plane crystal symmetry of teallite crystal on the nonlinear optical response with polarization-resolved THG emission. The anisotropic THG emission is investigated with a 1560 nm pulse laser source with the spot size of 1.5 µm. Figure 5a shows the transmission microscope image with the green THG emission from the 61 nm-thick teallite flake. The corresponding THG emission spectrum is shown in Fig. 5b, showing a peak at 520 nm which is exactly one third of the pump wavelength. To further reaffirm the THG process, the log-scale plot of THG emission power as a function of the incident pump power is shown in Fig. 5c, in which the cubic power law fit endorses the THG process. Next, the in-plane anisotropic THG emission as a function of the incident pump beam polarization in different teallite flakes are measured. The incident pump polarization is controlled by engaging the combination of a linear polarizer and a half-wave plate in the excitation path. It is worth noting that the anisotropic total THG emission pattern is symmetrical with respect to the optical axis along the x′-axis and its perpendicular direction along the y′-axis. Hence, the x′-and y′-components of THG emission power are filtered out by introducing an analyzer orientated in parallel and perpendicular to the crystal's optical axis. Figure 5d-g show the measured polarization-resolved THG emission power as per the incident pump beam polarization angle for four different flakes with thicknesses of 9, 24, 36, and 61 nm. The THG emission response is highly anisotropic with four-lobe patterns. The primary maxima of THG emission are collected along around − 40° and 140° which are almost along the optical axis in the [ 110] lattice direction of teallite crystal, which is consistent with the case of the measured linear optical absorption pattern, while the secondary maxima are recorded along its perpendicular direction at around 50° and 230°. The experimentally measured data points for x′-component, y′-component, and total THG emission power are shown with red squares, blue circles and black triangles, whereas the theoretical fits are shown with solid lines in the respective colors. We further corroborate the observed anisotropic THG emission with a theoretical nonlinear susceptibility model. Ternary teallite belongs to orthorhombic crystal family, therefore the contracted form of the third-order nonlinear susceptibility tensor can be written as 45,46 where the first term in subscript 1, 2 and 3 denotes x′, y′ and z′ respectively and the second subscript refers to the combination of three components as As per the experimental condition, teallite crystals are probed with linearly polarized electric field of − → E = x ′ (|E| cos θ) + y ′ (|E| sin θ) , where θ is the linear polarization angle relative to the crystal's optical axis. As the polarization resides in the x′-y′ plane only, no significant contribution from the z′-components is accounted. Hence, the nonlinear susceptibility tensor elements containing z′-components can be neglected, and the resultant the electric field components of THG emission and THG intensity components can then be expressed as 12 The theoretical fits from Eqs. (6) and (7) are plotted in the respective colors, showing a good agreement with the measured data. Further, the theoretical fits also enable us to retrieve the relative magnitudes of the nonlinear susceptibility tensor elements χ 11 , χ 18 , χ 22 and χ 29 . The average relative magnitudes for teallite flakes of different thicknesses are found as χ 11 : χ 18 : χ 22 : χ 29 = 1 : 0.037 : 0.86 : 0.029 . Also, the average value of THG anisotropic ratio |χ 11 | 2 /|χ 22 | 2 is extracted as 1.35. Finally, the third-order nonlinear susceptibility χ (3) value of teallite crystal is estimated by measuring the thickness-dependent THG emission. Figure 5h shows the recorded THG emission power as a function of the flake thickness. The average pump power is kept at 1.1 mW with a peak irradiance of 8.65 GW/cm 2 . Interestingly, the THG emission power gradually increases up to 4.3 pW for the 61 nm-thick flake and afterwards strong optical absorption starts playing a crucial role in attenuating the THG  www.nature.com/scientificreports/ signal propagation so that the THG emission exponentially decays for the relatively thicker teallite flakes. In fact, the exponential decay of THG emission enables us to extract the imaginary part of refractive index (k 3 ) at emission wavelength ( 3 = 520 nm) by fitting the thickness-dependent THG emission power equation where d is thickness of flake and A is a constant. Figure 5h shows the corresponding plot with the measured data (black squares) and the fitted curve (green curve) for k 3 = 1.42. We assume the refractive index of teallite is similar to that of herzenbergite, where n 3 = 4.39 at 520 nm and n 1 = 3.61 at 1560 nm along the armchair direction 47 . Thus, taking these values along with other experimental parameters into account ( P (ω) = 1.1 mW, laser pulse width τ = 90 fs, repetition rate f rep = 80 MHz, and spot size W = 1.5 µm at the fundamental wavelength 1 = 1560 nm), the magnitude of χ (3) can be estimated by the following formula 11 where n 1 and n 3 are the real part of refractive index at pump wavelength ( 1 ) and emission wavelength ( 3 ), �k = 6π 1 (n 1 − n 3 ) is the phase mismatch between the fundamental beam and the forward propagating THG emission beam in the transmission optical setup arrangement. Hence, the magnitude of χ (3) for teallite crystal is estimated as 3.49 × 10 -19 m 2 /V 2 , which is within same order of magnitude as the previously reported 2D layered materials.

Discussion
To conclude, we have demonstrated the polarization-dependent anisotropic optical responses from the mechanically exfoliated teallite flakes of various thicknesses. The exfoliated ternary teallite crystals are investigated using HRTEM and EDXS characterization tools to determine the structural and chemical composition. Further, the anisotropic Raman vibrational modes and linear optical properties are probed with polarization-resolved Raman spectroscopy and optical absorption spectroscopy, which determine the armchair direction and the optical axis of teallite crystals, respectively. In addition, teallite crystals are found to exhibit strong modulation of direct optical band gap depending on the incident linear polarization angle. Finally, we investigate anisotropic nonlinear optical properties of teallite crystals with the THG emission measurement and estimate the third-order nonlinear susceptibility value. Our findings also emphasize on the distorted lattice symmetry induced linear and nonlinear optical responses in teallite crystals, due to the element alloying in complex multi-element layered materials. In hindsight, we envisage that these outcomes would lead to a better understanding of the anisotropic optical responses in natural multi-element layered materials and can have implications in optical strain sensors, photodetectors, frequency modulators, encrypted signal processing, and other prototype nonlinear photonic applications.

Methods
Sample preparation. Teallite flakes of various thicknesses are mechanically exfoliated from bulk natural teallite mineral (from Monserrat, Oruro, Bolivia) using Nitto tape (SPV 224) and Scotch tape. The glass substrates with 1 cm × 1 cm are treated with acetone, deionized water and isopropanol followed by ultra-sonication for 20 min to remove the undesirable residues from the surface. These pretreated glass substrates are used to transfer the mechanically exfoliated flakes.
Polarization-resolved Raman spectroscopy. The sample is illuminated with a 632.8 nm He-Ne laser using a 40× objective lens (NA = 0.65) and the back-reflected signal is collected via the same objective lens to a spectrometer (Horiba, iHR 520) with a beam splitter and sets of mirrors. The elastic scattered light is rejected using an edge filter (Semrock, LP02-633RE-25) in the collection path. The desired linear polarization of excitation beam is obtained using a combination of a linear polarizer and a rotating half-wave plate in the illumination path. The collected signal is further passed through a linear polarization analyzer for resolving the parallel polarization component of the Raman spectra.

Polarization-resolved optical absorption spectroscopy.
To record the polarization-resolved optical absorption spectra, a broadband white light source (Thorlabs, SLS201L, 360-2600 nm) is passed through a combination of a linear polarizer and half-wave plate and then focused on the probed sample with a 80× objective lens (NA = 0.5). To measure the reflection spectrum, the back-reflected light is collected from the sample using the same objective lens and routed towards the spectrometer using a beam splitter and sets of mirrors, whereas in the case of the transmission spectrum, the transmitted light through the sample is collected using another 100× objective lens (NA = 0.7) and routed to the spectrometer using mirrors. Further, the reflection and transmission spectra are normalized with the light source spectrum to achieve the reflectance (R) and transmittance (T) spectra. Finally, the absorbance (A) spectrum is obtained by using the relation of A = 1 -R -T. We also estimate the measurement uncertainties in transmittance, reflectance, absorbance, and optical band gap. To perform this, the transmittance (T) and reflectance (R) spectra are recorded multiple times. The average values of these spectra in the respective configuration T avg and R avg are obtained, and the measurement errors of transmittance and reflectance δT and δR are approximated as the difference of maximum and minimum counts at each wavelength. Then the measurement uncertainties in transmittance and reflectance are estimated as δT/T avg and δR/R avg , and the uncertainty in absorbance is obtained as δA/A = δT/T avg + δR/R avg . Further, the uncertainties in transmittance and reflectance are used to estimate the uncertainty in absorption coefficient δα/α by using Eq. (3), and subsequently the uncertainty in optical band gap is estimated as δE g /E g = 2 · δα/α . According to www.nature.com/scientificreports/ our analysis, the measurement uncertainties in transmittance, reflectance and absorbance turn out to be around 0.14% (0.20%), 0.22% (0.31%), and 0.36% (0.51%) at the wavelength of 520 nm (720 nm). The estimated uncertainty in optical band gap measurement is then around 0.66%.

THG emission measurement.
The sample is pumped with a femtosecond laser source at the wavelength of 1560 nm (pulse width 90 fs, repetition rate 80 MHz) using a 40× objective lens (NA = 0.65). The transmitted THG signal is collected through a 100× objective lens (NA = 0.7) and routed towards the imaging camera and the spectrometer (Horiba, iHR 520) for recording the THG image and the corresponding spectra, respectively. To reject the incident pump beam, the collected signal is passed through a shortpass filter in the collection path.