Access the full text.
Sign up today, get DeepDyve free for 14 days.
Chemical exfoliation of MoS via Li-intercalation route has led to many desirable properties and spectacular applications due to the presence of a metastable state in addition to the stable H phase. However, the nature of the specific metastable phase formed, and its basic charge conduction properties have remained controversial. Using spatially resolved Raman spectroscopy (~1 m resolution) and photoelectron spectroscopy (~120 nm resolution), we probe such chemically exfoliated MoS samples in comparison to a mechanically exfoliated H phase sample and confirm that the dominant metastable state formed by this approach is a distorted T’ state with a small semiconducting gap. Investigating two such samples with different extents of Li residues present, we establish that Li ions, not only help to exfoliate MoS into few layer samples, but also contribute to enhancing the relative stability of the metastable state as well as dope the system with electrons, giving rise to a lightly doped small bandgap system with the T’ structure, responsible for its spectacular properties. 1. Introduction Amongst the non-graphene layered materials, transition metal dichalcogenides (TMDs) have been in the focus of the research community for quite a long time because of their fascinating properties and novel applications that emerge upon 2D confinement. Specifically, molybdenum disulphide (MoS ) in this series bears great promises for displaying wide range of electronic [1–5], mechanical [6–9], optical [10–15] and chemical properties [16,17] due to the presence of several polymorphs with distinctly different electronic structures [18–20] including the possibility of exotic properties such as the quantum spin Hall effect . MoS structure is built from two close- 2- 4+ packed planes of S ions, sandwiching the hexagonal arrangement of Mo layer. Depending on the relative positioning of the two sulphur planes with respect to each other, two basic polymorphs of MoS2 are realised, with the H phase having the ABA arrangement and the T phase with the 2- 4+ 2- ABC arrangement of the S -Mo -S trilayers. Thus, the coordination of sulphur ions around the central Mo sites are distinctly different with the H phase having a trigonal prismatic coordination and the T phase an octahedral one. While H is the thermodynamically ground state, octahedral T structure can undergo further distortion to minimize the energy, consistent with the simple molecular orbital picture . Such distortions lead to the formation of various superlattices with 1 2 superlattice with chains of dimers (T’), 2 2 (T’’) and 3 3 (T’’’) periodic units relative to the undistorted 1 1 (T) unit cell [23,24]. MoS , showing semiconducting property with a tuneable band gap of 1.2 – 1.9 eV depending on the number of layers present [25,26] in its H phase can undergo phase transition involving different routes, such as, plasma hot electron transfer [27,28], mechanical strain [29–31], electron-beam irradiation [18,19] and chemical reactions, often at the expense of its semiconducting property. The chemical routes normally involve non-aqueous alkali metal intercalation [32–34] followed by a water exfoliation step, sophisticated electrochemical control [15,35] or expansion of the interlayer distance by hydrothermal synthesis [36,37]. Although theoretical calculations predict that different superlattices of the metastable state show different properties such that T’ and T’’ are small band gap semiconductor [38,39], T’’’ is a ferroelectric insulator , it has been generally assumed that, all the above mentioned processes lead to the formation of the metallic T phase [27,41–43]. However, theoretical calculation predicts T-phase to be not a metastable phase, but an intrinsically unstable state with spontaneous distortions lowering its energy . Direct experimental evidence about the electronic structure of the metastable states has been obtained by probing the states close to Fermi level with photoemission spectroscopy . It was conclusively shown in this paper from both experimental and theoretical considerations that, the lithium intercalation in the MoS sheets leads to the formation of semiconducting T’ state instead of metallic T state. But lithium intercalated MoS has been extensively studied for different application purposes such as hydrogen evolution reaction [44,45] and supercapacitor [46,47] to name a couple, based deeply in the belief of its metallic property. Thus, there is a clear divergence of views along two mutually incompatible scenarios for the metastable phase found in relative abundance in chemically exfoliated MoS . One of these views is based on its many useful device properties, subscribing to a metallic and therefore, the undistorted T phase. The opposing view- point stresses the impossibility of the T phase due to its intrinsic instability, suggesting instead that the metastable phase is in fact a small bandgap, semiconducting T’ phase. The T’ phase is derived from the T phase by a simple dimerization, leading to alternate long and short Mo-Mo distances, as shown in Fig. S1 of the Supplementary Material. In all these discussions, both camps across the dividing line have generally ignored any possible role of Li ions that are used in the process of the chemical exfoliation. We focus our attention on this very aspect to establish that Li ions play a crucial role in this context and our results provide a natural resolution of the apparently conflicting viewpoints described above in terms of a metallic state achieved via Li ion induced doping of electrons into the small bandgap, semiconducting distorted T’ phase. 2. Methodologies Flakes of the stable H phase MoS were obtained by the micromechanical scotch tape exfoliation and were deposited on conducting indium doped tin oxide (ITO) substrates, to avoid the sample from getting charged during photoemission experiments. We characterised several such flakes using AFM and found that most of these have about 4-9 layers of MoS with a typical interlayer separation of about 0.65 nm [48,49](see Fig. S2 a-f in Supplementary Material). The substrate, with the flakes on top of it, are then sealed in a vial at inert atmosphere and treated with 1.6M n- butyl lithium in hexane solution for two days at room temperature. In this step Li MoS is formed. x 2 Without exposing the samples to the ambient, the unreacted n-butyl lithium is then removed from the sample by injecting hexane repeatedly inside the vial and washing it thoroughly. Finally, the water exfoliation step is done by removing the sample from the vial and washed several times with distilled water to remove any excess intercalated Li in the form of lithium hydroxide (LiOH). In order to produce a series of samples having different concentration of Li , we kept everything same in the above-mentioned process except the last water exfoliation step. We washed the samples with distilled water through different numbers of wash-cycles to remove the extra Li from the system, namely 12 times, 6 times, 4 times, 1 time and without any water-wash at all. We refer throughout the text to these samples as 12W, 6W, 4W, 1W and 0W, respectively. In order to ascertain the consequence of Li intercalation and subsequent washing with water in effectively separating the layers of MoS by increasing the van der Waals gap, we investigated several flakes before intercalation, after intercalation and after the water washing (see Fig. S3 a-f in Supplementary Material). These results show that the interlayer gap increases to typically 1.3-1.5 nm and 0.9-1.1 nm with the intercalation and washing, respectively; this markedly enhanced interlayer separation compared to that (0.65 nm) of the H stable phase ensures electronically decoupled single layers of MoS in the chemically treated samples, as reported earlier . The reaction of Li MoS with water is so vigorous that no trace of Li is found in photoemission x 2 spectra after washing the sample just once with distilled water (see Fig. S4 in Supplementary Material), leading to very similar results for the 12W, 6W, 4W and 1W samples, as shown in Fig. S5 in Supplementary Material in terms of Raman spectra of all compounds. Even electron spectroscopic results between the 1W and 12W samples are very similar, while being distinct from the 0W sample, as illustrated in Fig. S6 in Supplementary Material. Therefore, we discuss below results obtained only from 0W and 12W samples, where we have collected data with the highest signal-to-noise ratio. For the sake of comparison with a reference sample, we carried out a complete set of experiments with identical conditions on mechanically exfoliated MoS flakes without the lithium treatment; this sample is referred as H in this paper. Micro-Raman measurements were carried out with LabRAM HR Evolution (HORIBA) confocal microscope-based Raman spectrometer at room temperature with a beam spot of radius 1 m. All Raman spectra were recorded with the excitation laser wavelength and power of 532 nm and 1.2 mW, respectively. Scanning photoelectron microscopy (SPEM) experiments were performed at the ESCA microscopy beamline of Elettra Synchrotron laboratory, with hν = 520 eV, spatial resolution of 120 nm and energy resolution of 180 meV. While detailed spectral features were obtained from specific spots on samples, the spot size is limited by the photon beam size (~120 nm). Composition resolved maps and images of the samples were obtained over a much larger dimension, typically about 40 m x 40 m area, by restoring the photon beam across the sample face with the help of a piezo-driven sample mounting stage. All measurements were carried out -10 with a chamber pressure 10 mbar at ambient temperature. 3. Results and Discussions Raman spectroscopy has proved to be a very powerful technique to understand the vibrational fingerprint of 2D layered materials. In Fig. 1a. Raman spectra, normalised at A peak, obtained 1g from three different samples: (i) H (blue circles), (ii) 12W (red triangles) and (iii) 0W (black dots), are compared. In order to avoid conversion of the metastable phase to the stable one under the excitation source, it was found necessary to record all data using a very low laser power (1.2 mW), explaining the relatively noisy data; however, the peaks and their positions are unambiguously -1 1 identifiable in these spectra. H sample exhibits two peaks at 380 and 404 cm due to E and A 2g 1g responsible for in plane and out of plane opposite vibrations of S and Mo atoms, respectively. For intercalated samples in the lower wavenumber region, four additional peaks, commonly termed J -1 -1 -1 -1 (156 cm ), J (220 cm ), E (283 cm ) and J (330 cm ), are clearly observed. The appearance 2 1g 3  A b. a. 1g 0W 1 12W 2g 𝑴 = 𝜞 ′ 1g 380 400 420 -1 Raman Shift (cm ) E 3 𝑲 1g J 1  2 E 2g 𝜞 100 200 300 400 -1 Raman shift (cm ) Fig. 1. (a) Room temperature Raman Spectra for pristine H and both of the intercalated samples (0W and -1 12W). The inset shows the spectral features in the region of about 370-425 cm in an expanded scale for clarity; (b) Reciprocal space diagram illustrating the relationship between the basal plane projection of the Brillouin zone of H-MoS (solid black lines) and the Brillouin zone of a 2a a superlattice (red dotted 2 0 0 lines) corresponding to the T’ structure. Intensity (a. u.) Intensity (a. u.) of these additional Raman peaks only from chemically exfoliated samples in Fig. 1a is a clear evidence of additional polymorphs of MoS in such cases compared to the pure H phase of the mechanically exfoliated samples. Such results have been extensively reported in the literature and often these new Raman lines have been used to argue the presence of metallic T phase [45,50]. However, this can be easily seen to be erroneous in the following way. The point group describing the symmetry of a single MoS layer of undistorted T-phase is D . Applying arguments from 2 3d group theory, irreducible representation for normal modes of vibrational degrees of freedom can be written as: Γ (1𝑇 ) = 𝐴 + 𝐸 + 2𝐸 + 2𝐴 𝑣𝑖𝑏 1𝑔 𝑔 𝑢 2𝑢 By looking at the basis functions of each irreducible representations from the character table of D point group, A and E modes can be easily recognised as Raman active modes, whereas E 3d 1g g u and A modes are infra-red active. This analysis is also supported by several theoretical study 2u described in the literature . So, only two peaks should appear in Raman spectra of the undistorted T-phase MoS . In Fig. 1a, the presence of four peaks, namely, J , J , E and J peaks, 2 1 2 1g 3 clearly indicates that more than two vibrational normal modes are Raman active for the metastable state present in the chemically exfoliated samples. This is a strong evidence that the metastable state present in this case cannot be the undistorted T-phase. In fact, the frequencies of these four additional peaks in the Raman spectra of chemically exfoliated samples suggests the formation of the distorted T’ phase. The appearance of the additional peaks can be explained in terms of the presence of the superlattice in such distorted structures. Because of wave vector conservation, the peaks in Raman spectra are normally responsible for the normal modes vibrations at the zone centre (k = 0) . But if a superlattice is formed, the Brillouin zone is reduced, and zone boundary points of the undistorted lattice can coincide with the centre of the new zone. A typical example is shown in Fig. 1b, for the formation of 2 1 superlattice of the T’ phase. In this case, the M point of the original Brillouin zone is folded to the centre of the Brillouin zone of the distorted crystal. As a result, apart from point, vibrational modes from M can also emerge as Raman peaks. For the T’-phase irreducible representation for the vibrational normal modes can be written , ( ) Γ 1𝑇 ′ = 6𝐴 + 3𝐵 + 3𝐴 + 6𝐵 𝑣𝑖𝑏 𝑔 𝑔 𝑢 𝑢 where, A and B modes are Raman active. Therefore, group theory predicts the appearance of as g g many as nine peaks for the T’-MoS sample, which is in good agreement with the theoretical calculations . Within the interpretation of the formation of primarily T’ phase as the metastable polymorph in chemically exfoliated MoS , J peak is composed of two different phonon modes, 2 1 namely out of plane opposite vibration of each stripe of the Mo atoms inside the zigzag chain and in plane shearing vibration of one stripe of atoms with respect to another inside a chain. J2 peak is responsible for the phonon modes which tends to restore the original trigonal prismatic structure of H phase by equalising Mo - Mo distance between two zigzag chains. Though theoretically predicted  most intense peak is J , the much larger full width at half maxima (FWHM) has substantially reduced the height of the peak. E peak arises from the anti-symmetric stretching of 1g the adjacent S atom layers with respect to the Mo atom inside a zigzag chain. The slight out of plane phonon mode responsible for J peak tends to break each zigzag chain from its adjacent one in two stripes. It is interesting to note here the intense presence of the additional Raman peaks in the 0W sample and its significant reduction in the 12W sample. Since the only difference between these two samples is washing of the 12W sample with distilled water at room temperature in order to remove + + Li ions, it would appear that the presence of Li ions is essential to give rise to the extensive presence of the metastable state, evidenced by the strong presence of the characteristic, four additional Raman peaks (Fig. 1a) and the removal of the Li ions, as in 12W sample, leads to a sizable conversion of the metastable state to the stable H phase, as evidenced by the reduction of these Raman signals. We note that this is not a thermally-induced gradual conversion of the metastable polymorph to the stable one merely as a function of time, as explicitly checked by us by recording Raman spectra of the 0W sample after the lapse of a period of time equivalent to the time required to obtain 12W sample from a 0W one. Moreover, we also note that the Raman spectra of the 1W sample also show a similar reduction of the characteristic Raman signals. Thus, the presence of Li ions, not only facilitates reaching the metastable phase kinetically, but it also enhances its relative stability thermodynamically. We shall discuss this point further after presenting electron spectroscopic evidences in the following parts. Photoelectron spectra in the Mo 3d region for the two chemically exfoliated samples, namely 0W (black dots) and 12W (red triangles), are presented in Fig. 2a along with the spectrum (blue circles) from the mechanically exfoliated H sample using the highly focused photon beam from an arbitrary point on each sample; it was confirmed that spectral features remain qualitatively the same at different points on several flakes of the same sample. We can easily image the particular flake under investigation by mapping the integrated intensity of the Mo 3d peak between 226.84 5/2 and 230.7 eV over the entire sample area; the resulting Mo 3d intensity maps for flakes of 0W and 12W samples, most investigated for the purpose of this work, are shown in Fig. 2b and 2c. These images clearly show the well-defined presence of the two flakes on the ITO substrate, illustrated by the approximate boundary drawn with the cyan dotted lines. The reference H sample exhibits two narrow peaks at binding energies of 229.4 eV (Mo 3d5/2) and 232.5 eV (Mo 3d ) and a broad, low-intensity feature at ∼226.6 eV due to S 2s contributions (Fig. 3/2 a. b. c. 0W I II Max Max 0W 12W 12W 38.4 X 38.4 m 38.4 X 38.4 m Min Min 234 232 230 228 226 224 Binding energy (eV) e. d. Max Max 0W 12W 38.4 X 38.4 m 38.4 X 38.4 m Min Min Fig. 2. (a) Mo 3d core level spectra for 0W, 12W and H sample; Map of Mo 3d signal intensity over the 5/2 binding energy range of 226.84-230.7 eV in the (b) 0W sample and (c) 12W sample; Map of the intensity ratio II/I over respective energy windows shown in (a) for (d) 0W and (e) 12W sample. Intensity (a. u.) 2a). Upon comparing this spectrum with the spectra taken on the chemically treated samples, one observes two major differences: (i) spectra from the intercalated samples are prominently asymmetric and shifted towards the lower binding energy, (ii) FWHM of these spectra are large compared to the spectrum obtained from the pristine reference H-phase sample. These are clear indications that in case of intercalated samples, the dominant species is chemically shifted with respect to the reference sample, indicating an extensive presence of a metastable polymorph. Larger widths of the spectral features suggest the presence of multiple phases with finite chemical shifts between signals arising from different polymorphs. This is consistent with the observation described in ref. 32, where two distinct phases of MoS2, namely H and T’, were identified as contributing to the spectral features of the chemically exfoliated MoS in contrast to the almost exclusive presence of H- phase signal in mechanically exfoliated MoS . Following the same procedure, Mo 3d peak for the intercalated samples in Fig. 2a can be divided in two different 5/2 energy windows based on dominant contributions of the two different phases. It is clear that the intensity in the region I (228.9 eV – 230.1 eV) is dominated by contributions from the H-phase whereas contribution in region II (227.5 eV - 228.9 eV) is mainly from the metastable T’-phase with hardly any contribution from the stable H phase, as indicated by the spectral features of the pure H phase sample in Fig. 2a. Therefore, the contrast in the intensity ratio II/I, plotted in Fig. 2d and 2e for 0W and 12W, respectively, reveals the relative abundance of the metastable phase in the intercalated samples, with the lighter regions corresponding to the T’ MoS -rich regions. Fig. 2d shows that the T’-phase is reasonably uniformly distributed over the entire 0W sample in a. b. Raw data Raw data Mo 3d 3d Fitted data Mo 3d 5/2 Fitted data 3d 5/2 0W_H 12W_H 0W_T' 12W_T' S 2s S 2s 3d 3/2 3d 3/2 Error Error Background Background S 2s S 2s 234 232 230 228 226 224 234 232 230 228 226 224 Binding energy (eV) Binding energy (eV) Fig. 3. Mo 3d core level spectra (black open circle) and fitting of these spectra with two contributions for (a) 0W and (b) 12W samples. S 2s contribution is also shown. Intensity (a. u.) Intensity (a. u.) abundance. In contrast, Fig. 2e shows an overall markedly reduced abundance of the T’-phase in the 12W sample with some parts of the sample showing even more pronounced absence of the metastable phase than other parts, leading to a patchy distribution of the metastable phase. This is consistent with the conclusions drawn from the Raman studies with the 12W sample showing a significantly reduced average metastable phase contributions. This suggests that washing away of Li ions with water indeed helps to convert metastable T’-phase into the stable H-phase or conversely, retaining Li ions helps to stabilise the metastable state. In order to understand the nature and the origin of the metastable T’ phase, we identify the photoemission pixel with the largest T’ contributions in both samples, enclosed by the solid rectangle, shown in Fig. 2d and 2e and analyse spectra obtained from these pixels in detail. The Mo 3d spectra obtained from the specified spots are shown in Fig. 3a and 3b for 0W and 12W samples, respectively; we note that these spectra are the same ones shown in Fig. 2a overlapping with each other for an easy comparison. Our initial attempt to describe the Mo 3d spectral region as a sum of two chemically-shifted spectral contributions arising from the stable H and the metastable T’ phases did not work at all for the spectral features from the 0W sample, indicating the presence of additional spectral features in this case. In fact, a closer inspection of Fig. 2a, that overlaps the two spectra in Fig. 3a and 3b, reveals additional spectral weights for the 0W sample compared to the 12W one on the higher binding energy side, for example around 230 eV for the Mo 3d and around 233 eV for the Mo 3d regions. A reliable description of the spectral features 5/2 3/2 over the entire energy therefore requires additional components to account for such high binding energy satellite features accompanying the main peaks. Since spectral analysis of the S 2p region, discussed later in this paper, also requires similar satellite features with similar intensity and peak position relative to the main peak as in the Mo 3d region, it appears that there are prominent plasmon loss features accompanying all photoelectron signals from the 0W sample. Therefore, we fitted the spectral features obtained from 0W sample in Fig. 3a in terms of contributions from two phases, H and T’, each consisting of S 2s, Mo 3d and Mo 3d main peaks with each main peak 5/2 3/2 accompanied by plasmon loss satellite replicas with a fixed relative intensity and peak position; details of the spectral fitting have been described in greater detail in the SI. Allowing the relative intensity and the peak position of the plasmon satellite to vary, we obtain a very good description of the entire spectrum as shown in Fig. 3a with 0.23 and 1.3 eV as the relative intensity and the peak position for the plasmon satellite relative to the corresponding main peak, with the fit and the components also shown in Fig. 3a. An identical analysis of the spectrum recorded from the 12W sample yields an equally good fit of the entire spectrum with a relative intensity and peak position of 0.05 and 1.2 eV for the plasmon satellite, as shown in terms of the best fit and corresponding component spectra in Fig. 3b. Thus, there appears to be a considerable decrease in the plasmon satellite contribution in the 12W sample compared to the 0W sample. The lower intensity of the plasmon satellite in the 12W sample compared to the 0W sample suggests that the electron density responsible for this plasmonic signature is in higher concentration for the 0W sample compared to the 12W sample. In fact, we find that the spectral features for the 12W sample can be described reasonably well even without invoking any plasmonic satellite feature consistent with low relative intensity of 0.05 estimated in this satellite for the 12W sample here (see Fig. S7 in Supplementary Material). The resulting components appear with binding energies of Mo 3d peaks of H and T’ components 5/2 at 228.6 and 229.3 eV in 0W sample (Fig. 3a) and at 228.4 and 229.1 eV in 12W sample (Fig. 3b). It is clear that T’ is the dominant phase at the selected pixel for both 0W and 12W samples, with contributions of 75% and 67%, respectively, of the total Mo 3d signal. It is to be noted that this decrease of the T’ contribution in the 12W sample compared to that in the 0W sample refers to the spots identified to have the highest T’ contributions on respective samples. The overall decrease in the T’ contribution in the 12W sample is evident from the map of the relative T’ intensity for the two samples presented in Fig. 2d and 2e. Therefore, a significant amount of T’ a. b. c. Raw data Raw data 0W Fitted data 12W S 2p Fitted data S 2p S 2p H 0W_H 12W_H 0W_T' 12W_T' Defect state Defect state Error Error Background Background 166 164 162 160 166 164 162 160 166 164 162 160 Binding energy (eV) Binding energy (eV) Binding energy (eV) Fig. 4. (a) S 2p core level spectra for pristine H, 0W and 12W samples; Fitting of these spectra with three contributions for (b) 0W and (c) 12W samples. Intensity (a. u.) Intensity (a. u.) Intensity (a. u.) phase appears to convert back to the stable H phase on washing with water and the consequent removal of Li ions. We show the S 2p spectral region from 0W (black dots), 12W (red triangles), and H (blue circles) samples in Fig. 4a. While the S 2p spectrum of the H phase shows the expected spin- orbit split doublet structure, consisting of 2p and 2p signals with sharp features, the S spectra from the 3/2 1/2 Li-treated 0W and 12W samples are rather broad and featureless, indicating the presence of 2- multiple species of S ions in these samples. Substantial additional intensities on the lower binding energy side for the chemically treated samples compared to the H sample suggest the existence of the metastable T’ state as in the case of Mo 3d spectral region (see Fig. 2a). Additional spectral intensities on the higher binding energy side for the treated samples in Fig. 4a, with more enhanced intensity for the 0W sample, is most likely due to the presence of prominent contributions from the plasmonic loss satellites, as also observed in the Mo 3d spectral region. These close similarities prompt us to analyse the S 2p spectral region in a way similar to that for the Mo 3d region for the chemically treated samples, allowing for simultaneous presence of signals from two polymorphs, H and T’, along with plasmonic loss satellites accompanying main peaks with a constant intensity ratio and peak separation. Such an analysis, however, fails to provide a satisfactory description of the S 2p spectral region, as shown in Fig. S8 (in Supplementary Material) suggesting the presence of other species. In order to keep the description as simple as possible, we allowed for an additional 2- species of S to contribute in this spectral region. We find that the description in terms of 2- contributions from three species of S , namely originating from the H, the T’ and the additional S species, is sufficient to provide a consistent and convincing description of the entire spectral regions for both 0W and 12W samples. We show the outcome of this fitting procedure, resulting in the spectral decomposition for the S 2p regions of 0W and 12W samples in Fig. 4b and 4c. The three spectral components are marked as H (blue solid line), T’ (red dashed line), and the additional, third component (green dash-dotted line) termed as the “defect”. While the origin of H and T’ related signals are obvious, being related to patches of the chemically treated samples existing in the corresponding crystallographic forms, the appearance of the third component in the sulphur spectrum can be understood in the following manner. It is clear from our results that there is extensive coexistence of the H and T’ phases in such chemically treated samples. All T-derived phases, such as T, T’, T’’ and T’’’ phases, are structurally related to the H phase by a sliding 2- arrangement of the top S layer to change the ABA arrangement of the hexagonal planes in the H phase to the ABC type of arrangement of the layers in all T derived polymorphs. Extensive coexistence of H and T’ phases would then suggest small patches (< 120 nm of the photon beam size) of these two polymorphic forms, necessitating extensive coexistence of ABA and ABC type 2- arrangements in any given flake of chemically treated MoS . Thus, the top S layer would locally 2- shift between the A type and the C type arrangements with respect to the bottom S layer in nanoscopic domains. The domain boundaries will represent defect lines, necessarily strained to 2- 2- accommodate the local sliding arrangement of the S layer. The S ions along this domain 2- 2- boundaries represent a “defect” S ions with a distinct electronic structure compared to S ions within the H-type or the T’-type polymorphic forms. We note that S 2p region can be fitted reasonably well, considering the plasmons associated with it, in a way very similar to that of the Mo 3d fitting with relative intensities of 0.23 and 0.05 and relative peak positions of 1.5 and 1.2 eV with respect to the corresponding main peaks for the 0W and 12W samples, respectively. These values are almost same to that of the numbers obtained from Mo 3d analysis, indicating the reliability of the fit parameters extracted here and the description of the spectral features in terms of the components. Thus, both Mo 3d region and the S 2p region suggest that the plasmonic satellite contributions decrease with the removal of the Li ions via washing. In order to understand the nature of changes in the electronic structure between 0W and 12W samples and the consequent enhanced stability of the T’ phase in the 0W sample, we carefully compare the components, obtained from the fitting of Mo 3d spectra for both 0W and 12W samples in Fig. 5a by normalizing them at 3d peak; this procedure allows us to readily identify spectral 5/2 shifts between different components as well as between different samples. Here solid and dotted lines refer to 0W and 12W samples, respectively, whereas H and T’ phases have been indicated using red and blue colour. Therefore, solid blue and red lines denote the Mo 3d spectral positions from the T’ and H components of the 0W sample, while the dotted lines represent the corresponding components for the 12W sample. It is most interesting to note that the blue solid line and the dotted line, representing the T’ phase in 0W and 12W samples, respectively, show a shift of 0.2 eV for the dotted line of the 12W sample towards the lower binding energy compared to the blue solid line of the 0W sample. We note that this shift is evident even in the raw a. b. ~ 0.3 eV ~ 0.2 eV 0W_T' 0W_T' 0W_H 0W_H 12W_T' 12W_T' 12W_H 12W_H 166 164 162 160 234 232 230 228 Binding energy (eV) Binding energy (eV) c. d. 0W 12W 0.2 eV 2.5 2.0 1.5 1.0 0.5 0.0 Binding energy (eV) Fig. 5. H and T’ components of (a) Mo 3d core level spectra normalized at Mo 3d peak and (b) S 2p core 5/2 level spectra normalized at S 2p peak for both 0W and 12W samples; (c) Spectra near the Fermi edges 3/2 for pristine H, 0W and 12W sample; (d) Schematic diagram of the shift of the Fermi level due to electron doping of the system, reflecting as a shift in the binding energies of all levels. experimental spectra of 0W and 12W samples, shown in Fig. 2a. We also note that the H components shown by the red solid line and the dotted line in Fig. 5a also show the identical 0.2 eV shift to the lower binding energy for the dotted line corresponding to the H component in the 12W sample. Similar plots are shown in Fig. 5b for S 2p components after normalizing them at the 2p peaks. Once again, we observe similar shifts of 0.3 eV towards lower binding energy side 3/2 for all S 2p components of the 12W sample compared to those of 0W sample. Even the valence band edge of the 12W sample exhibits a 0.2 eV shift to the lower binding energy compared to that of the 0W sample, as illustrated in Fig. 5c. In Fig. 5c, we have also plotted the valence band edge region from the pristine H sample with blue circles, whereas the black dots and red triangles are used to plot the valence band edge of 0W and 12W samples, respectively. All three spectra shown in Fig. 5c exhibit negligible intensity at the Fermi energy (EF), clearly suggesting that none of these samples can be regarded as a metal, for example as would be expected from the undistorted Intensity (a. u.) Intensity (a. u.) Intensity (a. u.) T phase with a large, finite density of states at the E . Fig. 5c also shows that for pristine H sample the intensity is negligible down to 1.0 eV binding energy, indicating a sizeable bandgap which in turn leads to the conclusion that the spectra from the intercalated samples within the range of binding energies smaller than 1.0 eV are contributed only by the metastable, small bandgap T’ phase. The comparison of the spectral features clearly establishes that 0W sample has the leading edge of the valence band shifted towards the higher binding energy compared to that of the washed 12W sample by the energy separation comparable to that of the core level spectra. This consistent spectral shift of about 0.2 eV towards the higher binding energy for all spectral features from the 0W sample cannot be explained by any mechanism other than in terms of some electron doping of the 0W sample, thereby moving up the Fermi energy. Since all spectra are referenced to the Fermi energy in photoelectron spectroscopy, an upward shift of this reference energy, as would happen with electron doping, will appear as a constant shift of the same magnitude towards the higher binding energy for all spectra, as shown schematically in Fig. 5d. The electron doping is most likely dominated by the presence of Li ions, each ion donating one electron to MoS and therefore, determining the pinning of the Fermi energy to some level. With the removal of Li via washing also removes the doped electrons, thereby moving the Fermi energy down, this phenomenon reflecting itself as a shift of all spectroscopic features to a lower binding energy, as observed. We may also relate this electron doping via Li ions to the presence of prominent plasmon features only in the 0W sample. Clearly, the spectral features of the pure H sample show no signature of plasmonic loss satellites. Thoroughly washed 12W sample with Li ions below the detection limit shows an evidence of only a very low intensity plasmonic replica of all spectral features. Thus, it appears that the additionally doped electrons via the presence of Li ions on MoS , residing close to the bottom of the conduction band, contribute to the plasmonic loss features, observed prominently for the 0W sample. It is then easy to visualise that such electrons, doped via Li ions, pinning the Fermi energy near the bottom of the conduction band will also make chemically exfoliated MoS more conducting. This conductivity will depend on the extent of Li ion induced electron doping, but is also aided by the formation of the small bandgap (90 ± 40 meV) metastable T’ phase, such that the thermally excited free charge carriers from the Fermi energy to the conduction band will be significant in number. This scenario also suggests a possible explanation for our observation that the presence of Li ions is essential for the extensive presence of the metastable state, as evidenced by a significant depletion of the metastable state on + removal of Li ions via washing (see Fig. 1a). It is obvious that electron doping will prefer a small bandgap semiconductor like the T’ phase, rather than one like the H phase with a considerably higher energy conduction band. Thus, the presence of Li ions and the consequent electron doping of MoS tilts the energy balance slightly in favour of the doped T’ phase, this advantage being absent in the undoped case. The ready and irreversible conversion of the T’ phase to the stable H phase on various mild perturbations, such as warming the sample at around ~ 98°C , establishes that this phase is indeed a higher energy, or a metastable phase and the H phase is thermodynamically the ground state. So, the ability to reach this alternate state must be kinetically controlled. The fact that mechanical exfoliation does not yield this metastable phase and that the presence of Li ions is necessary to reach this state suggests that Li ions must be contributing to the kinetic pathways to reach the metastable state. After having reached the metastable state, the removal of Li ions by washing leads to the partial conversion of the metastable state to the stable H phase. Since the metastable state was already kinetically stabilised during the Li intercalation, the conversion of the metastable phase to the stable phase must be related to the change in the potential energy landscape with the removal of the Li ions, possibly lowering the barrier locally for the conversion to the stable phase at the same temperature. We also argue that the presence of a doped electron in presence of Li intercalation reduces the energy difference between the stable H phase with its large bandgap and the metastable T’ phase with its small bandgap in favour of the metastable phase, providing a plausible thermodynamic reason for the change in the local potential landscape + + between the Li ion doped state and the Li ion removed states. 4. Conclusion In conclusion, we have shown that chemically exfoliated MoS has an extensive presence of a distorted metastable T’ state that has a small bandgap compared to the stable H phase. Raman spectra clearly establish that this metastable state cannot be the undistorted T form, while the electron spectroscopy directly shows the absence of any significant metallic density of states at the Fermi energy. Investigating samples with different extent of washing to remove Li ions, we show that Li ions dope the small bandgap semiconducting metastable state to make these more conducting than in their pure state, explaining why such chemically exfoliated MoS has often been termed as metallic and this in turn has led to the erroneous conclusion of the undistorted T phase formation. Our results also show that these doped electrons facilitate the formation of the metastable, small bandgap, distorted T’ phase, enhancing their relative stability via doping, this is at the origin of the extensive formation of this metastable state via the chemical exfoliation route. Acknowledgement: The authors thank Smruti Rekha Mahapatra for her help in AFM measurements and G V Sai Manohar for initial Raman characterization (not included in this manuscript). The authors thank the Department of Science and Technology, Government of India, Indo-Italian Program of Cooperation, and International Centre for Theoretical Physics, Trieste, for supporting this research. D.P. acknowledges the Council of Scientific and Industrial Research for a student fellowship. D.D.S. thanks Jamsetji Tata Trust for support. Reference:  T. Li, G. Galli, Electronic properties of MoS nanoparticles, J. Phys. Chem. C. 111 (2007) 16192–16196, https://doi.org/10.1021/jp075424v.  R. Addou, S. McDonnell, D. Barrera, Z. Guo, A. Azcatl, J. Wang, H. Zhu, C.L. Hinkle, M. Quevedo-Lopez, H.N. Alshareef, L. Colombo, J.W.P. Hsu, R.M. Wallace, Impurities and Electronic Property Variations of Natural MoS Crystal Surfaces, ACS Nano. 9 (2015) 9124–9133, https://doi.org/10.1021/acsnano.5b03309.  G. Seifert, H. Terrones, M. Terrones, G. Jungnickel, T. Frauenheim, Structure and electronic properties of MoS nanotubes, Phys. Rev. Lett. 85 (2000) 146–149, https://doi.org/10.1103/PhysRevLett.85.146.  H. Pan, Y.W. Zhang, Tuning the electronic and magnetic properties of MoS2 nanoribbons by strain engineering, J. Phys. Chem. C. 116 (2012) 11752–11757, https://doi.org/10.1021/jp3015782.  A. Steinhoff, M. Rösner, F. Jahnke, T.O. Wehling, C. Gies, Influence of excited carriers on the optical and electronic properties of MoS , Nano Lett. 14 (2014) 3743–3748, https://doi.org/10.1021/nl500595u.  A. Castellanos-Gomez, M. Poot, G.A. Steele, H.S.J. Van Der Zant, N. Agraït, G. Rubio- Bollinger, Elastic properties of freely suspended MoS nanosheets, Adv. Mater. 24 (2012) 772–775, https://doi.org/10.1002/adma.201103965.  C. Ataca, H. Şahin, E. Aktuörk, S. Ciraci, Mechanical and electronic properties of MoS nanoribbons and their defects, J. Phys. Chem. C. 115 (2011) 3934–3941, https://doi.org/10.1021/jp1115146.  J.W. Jiang, H.S. Park, T. Rabczuk, Molecular dynamics simulations of single-layer molybdenum disulphide (MoS ): Stillinger-Weber parametrization, mechanical properties, and thermal conductivity, J. Appl. Phys. 114 (2013), https://doi.org/10.1063/1.4818414.  S. Xiong, G. Cao, Molecular dynamics simulations of mechanical properties of monolayer MoS , Nanotechnology. 26 (2015), https://dx.doi.org/10.1088/0957- 4484/26/18/185705.  J.P. Wilcoxon, P.P. Newcomer, G.A. Samara, Synthesis and optical properties of MoS and isomorphous nanoclusters in the quantum confinement regime, J. Appl. Phys. 81 (1997) 7934–7944, https://doi.org/10.1063/1.365367.  A. Molina-Sánchez, K. Hummer, L. Wirtz, Vibrational and optical properties of MoS : From monolayer to bulk, Surf. Sci. Rep. 70 (2015) 554–586, https://doi.org/10.1016/j.surfrep.2015.10.001.  M. Ye, D. Winslow, D. Zhang, R. Pandey, Y.K. Yap, Recent advancement on the optical properties of two-dimensional molybdenum disulfide (MoS ) thin films, Photonics. 2 (2015) 288–307, https://doi.org/10.3390/photonics2010288.  H. Shi, H. Pan, Y.W. Zhang, B.I. Yakobson, Quasiparticle band structures and optical properties of strained monolayer MoS and WS , Phys. Rev. B. 87 (2013) 1–8, 2 2 https://doi.org/10.1103/PhysRevB.87.155304.  Y. Li, Y. Rao, K.F. Mak, Y. You, S. Wang, C.R. Dean, T.F. Heinz, Probing symmetry properties of few-layer MoS and h-BN by optical second-harmonic generation, Nano Lett. 13 (2013) 3329–3333, https://doi.org/10.1021/nl401561r.  F. Xiong, H. Wang, X. Liu, J. Sun, M. Brongersma, E. Pop, Y. Cui, Li Intercalation in MoS2: In Situ Observation of Its Dynamics and Tuning Optical and Electrical Properties, Nano Lett. 15 (2015) 6777–6784, https://doi.org/10.1021/acs.nanolett.5b02619.  W. Chen, E.J.G. Santos, W. Zhu, E. Kaxiras, Z. Zhang, Tuning the electronic and chemical properties of monolayer MoS adsorbed on transition metal substrates, Nano Lett. 13 (2013) 509–514, https://doi.org/10.1021/nl303909f.  M. Chhowalla, H.S. Shin, G. Eda, L.J. Li, K.P. Loh, H. Zhang, The chemistry of two- dimensional layered transition metal dichalcogenide nanosheets, Nat. Chem. 5 (2013) 263–275, https://doi.org/10.1038/nchem.1589.  S. Kretschmer, H.P. Komsa, P. Bøggild, A. V. Krasheninnikov, Structural Transformations in Two-Dimensional Transition-Metal Dichalcogenide MoS under an Electron Beam: Insights from First-Principles Calculations, J. Phys. Chem. Lett. 8 (2017) 3061–3067, https://doi.org/10.1021/acs.jpclett.7b01177.  Y.C. Lin, D.O. Dumcenco, Y.S. Huang, K. Suenaga, Atomic mechanism of the semiconducting-to-metallic phase transition in single-layered MoS , Nat. Nanotechnol. 9 (2014) 391–396, https://doi.org/10.1038/nnano.2014.64.  X. Zhao, S. Ning, W. Fu, S.J. Pennycook, K.P. Loh, Differentiating Polymorphs in Molybdenum Disulfide via Electron Microscopy, Adv. Mater. 30 (2018) 1–11, https://doi.org/10.1002/adma.201802397.  X. Qian, J. Liu, L. Fu, J. Li, Quantum spin Hall effect in two-dimensional transition metal dichalcogenides, Science 346 (2014) 1344–1348, https://doi.org/10.1126/science.1256815.  I. Song, C. Park, H.C. Choi, Synthesis and properties of molybdenum disulphide: From bulk to atomic layers, RSC Adv. 5 (2015) 7495–7514, https://doi.org/10.1039/C4RA11852A.  T. Hu, R. Li, J. Dong, A new (2 × 1) dimerized structure of monolayer 1T-molybdenum disulfide, studied from first principles calculations, J. Chem. Phys. 139 (2013), https://doi.org/10.1063/1.4827082.  M. Kan, J.Y. Wang, X.W. Li, S.H. Zhang, Y.W. Li, Y. Kawazoe, Q. Sun, P. Jena, Structures and phase transition of a MoS monolayer, J. Phys. Chem. C. 118 (2014) 1515– 1522, https://doi.org/10.1021/jp4076355.  K.F. Mak, C. Lee, J. Hone, J. Shan, T.F. Heinz, Atomically thin MoS : A new direct-gap semiconductor, Phys. Rev. Lett. 105 (2010) 2–5, https://doi.org/10.1103/PhysRevLett.105.136805.  C. Lee, H. Yan, L.E. Brus, T.F. Heinz, J. Hone, S. Ryu, Anomalous lattice vibrations of single- and few-layer MoS , ACS Nano. 4 (2010) 2695–2700, https://doi.org/10.1021/nn1003937.  Y. Kang, S. Najmaei, Z. Liu, Y. Bao, Y. Wang, X. Zhu, N.J. Halas, P. Nordlander, P.M. Ajayan, J. Lou, Z. Fang, Plasmonic Hot Electron Induced Structural Phase Transition in a MoS Monolayer, Adv. Mater. 26 (2014) 6467–6471, https://doi.org/10.1002/adma.201401802.  Y. Kang, Y. Gong, Z. Hu, Z. Li, Z. Qiu, X. Zhu, P.M. Ajayan, Z. Fang, Plasmonic hot electron enhanced MoS photocatalysis in hydrogen evolution, Nanoscale. 7 (2015) 4482– 4488, https://doi.org/10.1039/C4NR07303G.  X. Li, M. Wu, B. Xu, R. Liu, C. Ouyang, Compressive strain induced dynamical stability of monolayer 1T-MX ( M = Mo , W ; X = S , Se ), Mater. Res. Express. 4 (2017) 115018, https://doi.org/10.1088/2053-1591/aa9762.  A. Castellanos-Gomez, R. Roldán, E. Cappelluti, M. Buscema, F. Guinea, H.S.J. Van Der Zant, G.A. Steele, Local strain engineering in atomically thin MoS , Nano Lett. 13 (2013) 5361–5366, https://doi.org/10.1021/nl402875m.  Z. Ding, Q.X. Pei, J.W. Jiang, Y.W. Zhang, Manipulating the thermal conductivity of monolayer MoS via lattice defect and strain engineering, J. Phys. Chem. C. 119 (2015) 16358–16365, https://doi.org/10.1021/acs.jpcc.5b03607.  B. Pal, A. Singh, G. Sharada, P. Mahale, A. Kumar, S. Thirupathaiah, H. Sezen, M. Amati, L. Gregoratti, U. V. Waghmare, D.D. Sarma, Chemically exfoliated MoS layers: Spectroscopic evidence for the semiconducting nature of the dominant trigonal metastable phase, Phys. Rev. B. 96 (2017) 1–7, https://doi.org/10.1103/PhysRevB.96.195426.  K. Leng, Z. Chen, X. Zhao, W. Tang, B. Tian, C.T. Nai, W. Zhou, K.P. Loh, Phase Restructuring in Transition Metal Dichalcogenides for Highly Stable Energy Storage, ACS Nano. 10 (2016) 9208–9215, https://doi.org/10.1021/acsnano.6b05746  G. Eda, H. Yamaguchi, D. Voiry, T. Fujita, M. Chen, M. Chhowalla, Photoluminescence from chemically exfoliated MoS , Nano Lett. 11 (2011) 5111–5116, https://doi.org/10.1021/nl201874w.  M. Azhagurajan, T. Kajita, T. Itoh, Y.G. Kim, K. Itaya, In Situ Visualization of Lithium Ion Intercalation into MoS Single Crystals using Differential Optical Microscopy with Atomic Layer Resolution, J. Am. Chem. Soc. 138 (2016) 3355–3361, https://doi.org/10.1021/jacs.5b11849.  D. Wang, X. Zhang, S. Bao, Z. Zhang, H. Fei, Z. Wu, Phase engineering of a multiphasic 1T/2H MoS catalyst for highly efficient hydrogen evolution, J. Mater. Chem. A. 5 (2017) 2681–2688, https://doi.org/10.1039/C6TA09409K.  Z. Liu, Z. Gao, Y. Liu, M. Xia, R. Wang, N. Li, Heterogeneous Nanostructure Based on 1T-Phase MoS for Enhanced Electrocatalytic Hydrogen Evolution, ACS Appl. Mater. Interfaces. 9 (2017) 25291–25297, https://doi.org/10.1021/acsami.7b05775.  H.L. Zhuang, M.D. Johannes, A.K. Singh, R.G. Hennig, Doping-controlled phase transitions in single-layer MoS , Phys. Rev. B. 96 (2017) 1–8, https://doi.org/10.1103/PhysRevB.96.165305  M. Calandra, Chemically exfoliated single-layer MoS : Stability, lattice dynamics, and catalytic adsorption from first principles, Phys. Rev. B. 88 (2013) 1–6, https://doi.org/10.1103/PhysRevB.88.245428.  S.N. Shirodkar, U. V. Waghmare, Emergence of ferroelectricity at a metal-semiconductor transition in a 1T monolayer of MoS , Phys. Rev. Lett. 112 (2014) 1–5, 2 https://doi.org/10.1103/PhysRevLett.112.157601.  T. Hong, B. Chamlagain, S. Hu, S.M. Weiss, Z. Zhou, Y.Q. Xu, Plasmonic hot electron induced photocurrent response at MoS -metal junctions, ACS Nano. 9 (2015) 5357–5363, https://doi.org/10.1021/acsnano.5b01065.  W. Wang, A. Klots, D. Prasai, Y. Yang, K.I. Bolotin, J. Valentine, Hot Electron-Based Near-Infrared Photodetection Using Bilayer MoS , Nano Lett. 15 (2015) 7440–7444, https://doi.org/10.1021/acs.nanolett.5b02866.  J.S. Kim, J. Kim, J. Zhao, S. Kim, J.H. Lee, Y. Jin, H. Choi, B.H. Moon, J.J. Bae, Y.H. Lee, S.C. Lim, Electrical Transport Properties of Polymorphic MoS , ACS Nano. 10 (2016) 7500–7506, https://doi.org/10.1021/acsnano.6b02267.  D. Voiry, M. Salehi, R. Silva, T. Fujita, M. Chen, T. Asefa, V.B. Shenoy, G. Eda, M. Chhowalla, Conducting MoS nanosheets as catalysts for hydrogen evolution reaction, Nano Lett. 13 (2013) 6222–6227, https://doi.org/10.1021/nl403661s.  M.A. Lukowski, A.S. Daniel, F. Meng, A. Forticaux, L. Li, S. Jin, Enhanced hydrogen evolution catalysis from chemically exfoliated metallic MoS nanosheets, J. Am. Chem. Soc. 135 (2013) 10274–10277, https://doi.org/10.1021/ja404523s.  D. Sarkar, D. Das, S. Das, A. Kumar, S. Patil, K.K. Nanda, D.D. Sarma, A. Shukla, Expanding Interlayer Spacing in MoS for Realizing an Advanced Supercapacitor, ACS Energy Lett. 4 (2019) 1602–1609, https://doi.org/10.1021/acsenergylett.9b00983.  M. Acerce, D. Voiry, M. Chhowalla, Metallic 1T phase MoS nanosheets as supercapacitor electrode materials, Nat. Nanotechnol. 10 (2015) 313–318, https://doi.org/10.1038/nnano.2015.40.  K.D. Rasamani, F. Alimohammadi, Y. Sun, Interlayer-expanded MoS , Mater. Today. 20 (2017) 83–91, https://doi.org/10.1016/j.mattod.2016.10.004.  H.S.S. Ramakrishna Matte, A. Gomathi, A.K. Manna, D.J. Late, R. Datta, S.K. Pati, C.N.R. Rao, MoS and WS analogues of graphene, Angew. Chemie - Int. Ed. 49 (2010) 2 2 4059–4062, https://doi.org/10.1002/anie.201000009.  Q. Ding, F. Meng, C.R. English, M. Cabán-Acevedo, M.J. Shearer, D. Liang, A.S. Daniel, R.J. Hamers, S. Jin, Efficient photoelectrochemical hydrogen generation using heterostructures of Si and chemically exfoliated metallic MoS , J. Am. Chem. Soc. 136 (2014) 8504–8507, https://doi.org/10.1021/ja5025673.  M. Plischke, K.K. Bardhan, R. Leonelli, C. Irwin, Intercalate order-disorder transition in stage-2 Ag TiS , Can. J. Phys. 61 (1983) 397–404, https://doi.org/10.1139/p83-050. x 2  S.J.R. Tan, S. Sarkar, X. Zhao, X. Luo, Y.Z. Luo, S.M. Poh, I. Abdelwahab, W. Zhou, T. Venkatesan, W. Chen, S.Y. Quek, K.P. Loh, Temperature- and Phase-Dependent Phonon Renormalization in 1T’-MoS2, ACS Nano. 12 (2018) 5051–5058, https://doi.org/10.1021/acsnano.8b02649.  S. Ghatak, S. Mukherjee, M. Jain, D.D. Sarma, A. Ghosh, Microscopic origin of low frequency noise in MoS field-effect transistors, APL Mater. 2 (2014), https://doi.org/10.1063/1.4895955.  J. Heising, M.G. Kanatzidis, Exfoliated and restacked MoS and WS : Ionic or neutral 2 2 species? Encapsulation and ordering of hard electropositive cations, J. Am. Chem. Soc. 121 (1999) 11720–11732, https://doi.org/10.1021/ja991644d. Supplementary Material On the origin of metallicity and stability of the metastable phase in chemically exfoliated MoS a a a,+ b b Debasmita Pariari , Rahul Mahavir Varma , Maya N. Nair , Patrick Zeller , Matteo Amati , b c a,* Luca Gregoratti , Karuna Kar Nanda , D. D. Sarma Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru 560012, India Elettra-Sincrotrone Trieste S.C.p.A., SS14, km 163.5 in AREA Science Park, 34149 Basovizza, Trieste, Italy. Materials Research Centre, Indian Institute of Science, Bengaluru 560012, India Corresponding author: email@example.com Present address: CUNY Advanced Science Research Center, 85 St. Nicholas Terrace, New York, NY 10031, USA 1. Structure of undistorted (T) and distorted (T’) octahedral phase: a. b. Fig. S1. Side views of (a) T and (b) T’ phase of MoS where red and yellow balls indicate Mo and S atoms, respectively. Six S atoms are arranged around a Mo atom in a pure and distorted octahedral fashion in T and T’ phases, respectively. For the T phase, Mo-Mo distance is same throughout the chain in contrast to the T’ phase, where dimerization of consecutive Mo atoms leads to the formation of zigzag arrangements of the Mo ions. The side views for both T and T’ phases are shown in Fig. S1a and 2 4+ S1b. These structural changes, representing a Jahn-Teller distortion of the d configuration of Mo from a regular octahedron to a distorted one, changes the ground state electronic structure of MoS from being metallic in T phase to a small bandgap semiconductor in the T’ phase. Moreover, total energy calculations suggest that the distorted phase is relatively more stable compared to the undistorted state; in fact, the undistorted state is found to be unstable and is expected to spontaneously distort. 2. AFM study on the pristine and intercalated flakes: b. c. a. e. f. d. 5.0 3 -2 9 L 8 L 5 L 4.5 2 -3 4.0 -4 3.5 -5 3.0 -1 2.5 -6 -2 2.0 -7 -3 1.5 -4 -8 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.0 0.1 0.2 0.3 Distance (m) Distance (m) Distance (m) Fig. S2. (a), (b), (c) Representative AFM images of the pristine H samples and (d), (e), (f) corresponding height profiles of the flakes represented in (a), (b) and (c) respectively through the cyan lines. Mechanically exfoliated H samples (Fig. S2) were found primarily in the range of 4-9 layers with characteristic layer separation of about 0.65 nm, as reported in the literature [1,2]. In order to illustrate the effect of Li ion intercalation as well as washing these intercalated flakes, we show in Fig. S3 an individual flake of pristine H (also shown in Fig. S2c) sample before intercalation (Fig. S3a), after Li intercalation (Fig. S3b) and after washing with water (Fig. S3c). It clearly shows that the interlayer separation is increased to nearly 1.7 nm on intercalation, which reduces somewhat to about 1.1 nm on substantial removal of Li ions via washing, as has been reported in the literature . These hugely enhanced interlayer separation decouples successive layers, leading to essentially non-interacting single layers of MoS in these chemically treated samples. Height (nm) Height (nm) Height (nm) c. b. a. f. d. e. -6 -8 -2 9 L 9 L 9 L -4 Interlayer Interlayer Interlayer 1 -10 -6 separation separation separation -8 -12 0.65 nm -10 1.7 nm 1.08 nm -1 -14 -12 -2 -14 -16 -3 -16 -4 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 Distance (m) Distance (m) Distance (m) Fig. S3. Representative AFM images of an individual flake (a) before intercalation (H), (b) after Li intercalation (0W) and (c) after washing (12W) with corresponding height profiles represented by the cyan lines shown in (d), (e) and (f) respectively. 3. Evidence for the presence of lithium ions (Li ): In Fig. S4, photoelectron spectra in the Li 1s region for the four chemically exfoliated samples, namely 0W, 1W, 6W and 12W, are presented along with the spectrum from the mechanically b. a. 1W 0W (sample 1) Li 1s 6W Li 1s 0W (sample 2) 12W 60 58 56 54 52 50 48 60 58 56 54 52 50 48 Binding energy (eV) Binding energy (eV) Fig. S4. Li 1s core level spectra for (a) 0W samples and (b) 1W, 6W and 12W and mechanically exfoliated pristine H samples. Height (nm) Intensity (a. u.) Height (nm) Intensity (a. u.) Height (nm) exfoliated H sample. The 0W sample exhibits one peak, responsible for Li 1s, at binding energy of 55.8 eV, shown in Fig. S4a normalized at the peak position. Two representative data, indicated by black and red lines, taken on different samples, are shown here in order to prove the reproducibility of the data. The spectrum in the binding energy window responsible for Li 1s, obtained from other chemically exfoliated samples apart from 0W and pristine H sample, are compared in Fig. S4b. It is quite clear that, only the noise level can be visible in Li 1s binding energy window for these samples, which further proves that the intercalated Li ions react with water so vigorously that, even after a single wash the concentration level of Li goes beyond the detection limit. 4. Raman spectra of the samples with different washing cycles: We carried out preliminary studies for several control groups, namely 1W, 4W and 6W, in addition to 0W and 12W samples. The similarity of all washed samples is ascertained by the Raman spectra, (shown in Fig. S5). 0W 1W 4W 6W 12W 100 200 300 400 500 -1 Raman Shift (cm ) Fig. S5. Comparison of room temperature Raman spectra for the samples with different washing cycles -1 normalized at 405 cm (A peak). 1g Intensity (a. u.) 5. Photoelectron spectra of the samples with different washing cycles: a. b. 0W 0W 3d 5/2 Mo 3d S 2p 1W 1W 12W 12W 3d 3/2 S 2s 234 232 230 228 226 224 165 164 163 162 161 160 159 Binding energy (eV) Binding energy (eV) Fig. S6. Comparison of photoelectron spectra for 0W, 1W and 12W samples in the energy range of (a) Mo 3d and (b) S 2p core level. The similarity of the two limiting washed samples, namely 1W and 12W, are distinct from the 0W sample in terms of their photoelectron spectroscopic results as shown in Fig. S6a and S6b. Intensity (a. u.) Intensity (a. u.) 6. Analysis of Raman spectrum from Group Theory point of view: The presence of extra Raman peaks in intercalated samples can be explained with the help of Group Theory. The point groups, describing the symmetry of single MoS layer of pristine H, T and T’ phases, are D3h, D3d and C2h respectively. For each of the above mentioned point groups, irreducible representations corresponding to the normal modes of vibrations at zone centre (k = 0) along with the Raman and IR active modes are given below: Table 1: Point group symmetry for different phases of MoS along with the comparison between number of Raman active modes predicted from Group Theory and experimental observations. Point Irreducible representation IR Raman No. of No. of prominent group active active Raman Raman peaks modes observed experimentally ′ ′ ′′ ′ ′′ ′ ′ ′′ D ( ) 3 Mostly 2. Presence 3h Γ H = 𝐴 + 2𝐸 + 2𝐴 𝐸 , 𝐴 𝐴 , 𝐸 , 𝐸 𝑣𝑖𝑏 1 2 2 1 ′′ ′′ of very weak 𝐸 + 𝐸 mode i.e. 𝐸 peak 1𝑔 is reported in few studies[4,5]. ( ) D Γ T = 𝐴 + 𝐸 + 2𝐸 𝐸 , 𝐴 𝐴 , 𝐸 2 2 peaks 3d 𝑣𝑖𝑏 1𝑔 𝑔 𝑢 𝑢 2𝑢 1𝑔 𝑔 + 2𝐴 2𝑢 C Γ (T′) = 6𝐴 + 3𝐵 𝐴 , 𝐵 𝐴 , 𝐵 9 6 peaks*[6,7] 2h 𝑣𝑖𝑏 𝑔 𝑔 𝑢 𝑢 𝑔 𝑔 + 3𝐴 + 6𝐵 𝑢 𝑢 * Till now as per our knowledge there is no experimental report of stabilizing pure T’ phase. It is always coexisted with stable H phase. Therefore, it is inevitable to get signature of H phase in Raman spectroscopy, from the samples containing T’ phase. For such samples, three peaks responsible for H phase are accompanied by three more peaks, presumably coming from the presence of T’ phase. 7. Spectral fitting of 12W sample without plasmonic peak: b. a. Raw data Raw data Mo 3d 3d S 2p 5/2 Fitted data Fitted data 12W_2H 12W_2H 12W_1T' 12W_1T' 3d 3/2 Defect state S 2s Error Error Background Background S 2s 234 232 230 228 226 224 166 164 162 160 Binding energy (eV) Binding energy (eV) Fig. S7. Fitting of (a) Mo 3d and (b) S 2p core level spectra for 12W sample without considering the plasmon loss features. Both Mo 3d and S 2p spectra obtained from 12W sample can be fitted with good agreement without considering the plasmon loss features, showing in Fig. S7. This clearly indicates towards almost non-existent nature of the plasmon satellites in 12W sample. Intensity (a. u.) Intensity (a. u.) 8. Fitting the core level spectra: a. Raw data Mo 3d Fitted data 3d 5/2 Mo 3d S 2s Difference Background 3d 3/2 S 2s 234 232 230 228 226 224 Binding energy (eV) b. c. Raw data Raw data S 2p S 2p Fitted data Fitted data 0W_2H 12W_2H 0W_1T' 12W_1T' Error Background Background Error 166 164 162 160 166 164 162 160 Binding energy (eV) Binding energy (eV) Fig. S8. (a) Fitted Mo 3d core level spectra obtained from mechanically exfoliated sample; S 2p core level spectra have been fitted with two components for (b) 0W and (c) 12W samples. As there is only H phase present in mechanically exfoliated sample, the spectra can be decomposed in Mo 3d , Mo 3d and S 2s component, showing in Fig. S8a. The shape of each peak is defined 5/2 3/2 by a Lorentzian function convoluted by a Gaussian function, representing the lifetime and resolution broadening respectively. Throughout our analysis we have used a single Gaussian broadening for all components as the source of the resolution broadening lies in the instrumental artefact. We have imposed the constrain of fixed relative intensity between each spin orbit coupling pair i.e. the intensity ratio of 3d to 3d is 1.5 and similarly the intensity ratio between 2p and 5/2 3/2 3/2 2p is 2, along with the constant energy separation between each of these spin orbit pairs. We 1/2 have used the least number of components required for a reasonable description of the Mo 3d and Intensity (a. u.) Intensity (a. u.) Intensity (a. u.) S 2p spectra obtained from chemically exfoliated samples. As shown in Fig. S8a, only one component is sufficient to fit the spectra, while the Mo 3d spectra for both 0W and 12W samples can be fitted with the help of only two components namely, H and T’ phase. It appears to be impossible to fit S 2p spectra assuming only two components as shown in Fig. S8b and S8c. For both 0W and 12W samples, the relative abundance of the H component is higher compared to that of the T’ component, which is contradicting the results obtained from Mo 3d fitting. The spin orbit coupling components in Mo 3d spectra are clearly distinguishable, indicating a much reliable fit in that case. On the contrary, 2p and 2p peaks are merged together for both the intercalated 3/2 1/2 samples, making the fitting more difficult. In this scenario it is obvious to stick to the relative intensity value of the two components obtained from Mo 3d fitting, which in turn makes the S 2p spectral fitting with two components obsolete. As described in the main text, in order to take care of the extra intensity, appearing in the higher binding energy side of each peaks in Mo 3d spectra of the 0W sample, we have assumed the presence of plasmon satellites. That can be easily done by defining a function, which is the summation of four Lorentzian functions. It is important to note that we have kept the intensity ratio fixed between two spin orbit coupling pairs of both the main and the plasmon peaks but the relative intensity between a main peak and its plasmon has been allowed to vary. Similarly, the peak positions have been varied in such a way that the distance from the main peak to its plasmon for both the spin orbit coupling components remain same. References:  M. Chhowalla, H.S. Shin, G. Eda, L.J. Li, K.P. Loh, H. Zhang, The chemistry of two- dimensional layered transition metal dichalcogenide nanosheets, Nat. Chem. 5 (2013) 263–275. https://doi.org/10.1038/nchem.1589.  D. Wang, X. Zhang, S. Bao, Z. Zhang, H. Fei, Z. Wu, Phase engineering of a multiphasic 1T/2H MoS catalyst for highly efficient hydrogen evolution, J. Mater. Chem. A. 5 (2017) 2681–2688. https://doi.org/10.1039/C6TA09409K.  X. Geng, W. Sun, W. Wu, B. Chen, A. Al-Hilo, M. Benamara, H. Zhu, F. Watanabe, J. Cui, T.P. Chen, Pure and stable metallic phase molybdenum disulfide nanosheets for hydrogen evolution reaction, Nat. Commun. 7 (2016) 1–7. https://doi.org/10.1038/ncomms10672.  S. Jiménez Sandoval, D. Yang, R.F. Frindt, J.C. Irwin, Raman study and lattice dynamics of single molecular layers of MoS , Phys. Rev. B. 44 (1991) 3955–3962. https://doi.org/10.1103/PhysRevB.44.3955.  G. Zhu, J. Liu, Q. Zheng, R. Zhang, D. Li, D. Banerjee, D.G. Cahill, Tuning thermal conductivity in molybdenum disulfide by electrochemical intercalation, Nat. Commun. 7 (2016) 1–9. https://doi.org/10.1038/ncomms13211.  F. Xiong, H. Wang, X. Liu, J. Sun, M. Brongersma, E. Pop, Y. Cui, Li Intercalation in MoS2: In Situ Observation of Its Dynamics and Tuning Optical and Electrical Properties, Nano Lett. 15 (2015) 6777–6784. https://doi.org/10.1021/acs.nanolett.5b02619.  B. Pal, A. Singh, G. Sharada, P. Mahale, A. Kumar, S. Thirupathaiah, H. Sezen, M. Amati, L. Gregoratti, U. V. Waghmare, D.D. Sarma, Chemically exfoliated Mo S2 layers: Spectroscopic evidence for the semiconducting nature of the dominant trigonal metastable phase, Phys. Rev. B. 96 (2017) 1–7. https://doi.org/10.1103/PhysRevB.96.195426.
Condensed Matter – arXiv (Cornell University)
Published: May 2, 2020
Access the full text.
Sign up today, get DeepDyve free for 14 days.