Anomalous metallic phase in molybdenum disulphide induced via gate-driven organic ion intercalation

Transition metal dichalcogenides exhibit rich phase diagrams dominated by the interplay of superconductivity and charge density waves, which often result in anomalies in the electric transport properties. Here, we employ the ionic gating technique to realize a tunable, non-volatile organic ion intercalation in bulk single crystals of molybdenum disulphide (MoS$_{2}$). We demonstrate that this gate-driven organic ion intercalation induces a strong electron doping in the system without changing the pristine $2H$ crystal symmetry and triggers the emergence of a re-entrant insulator-to-metal transition. We show that the gate-induced metallic state exhibits clear anomalies in the temperature dependence of the resistivity with a natural explanation as signatures of the development of a charge-density wave phase which was previously observed in alkali-intercalated MoS$_{2}$. The relatively large temperature at which the anomalies are observed ($\sim$$150$ K), combined with the absence of any sign of doping-induced superconductivity down to $\sim$$3$ K, suggests that the two phases might be competing with each other to determine the electronic ground state of electron-doped MoS$_2$.

Although in most cases the effects of ionic gating vanish when the gate voltage is removed, it has recently been shown that non-volatile charge doping of macroscopic bulk specimens of layered crystals is attainable via gatedriven intercalation of either hydrogen [61][62][63][64] or organic ions [65]. The latter has been proved to be a very effective tool to tune the electronic ground state of layered FeSe [65] and to boost its superconducting transition temperature [66,67] well beyond that achieved via other means of bulk electron doping. These results make the investigation of the effects of gate-driven organic-ion intercalation in other layered materials highly desirable.
In this work, we employ the ionic gating technique to intercalate organic ions in bulk 2H-MoS 2 single crystals at room temperature. We find that the gate-induced intercalation is partially non-volatile and can be tuned by varying the amount of time during which the fully intercalated samples are kept in the ionic liquid after removal of the gate voltage. Our organic ion-intercalated MoS 2 crystals retain the 2H crystal symmetry and undergo a doping-induced insulatorto-metal transition, with clear anomalies emerging in the temperature-dependence of the resistivity around ∼ 150 K; all samples show no evidence of superconductivity down to ∼ 3 K, but an incipient re-entrant transition to an insulating state is observed at the highest doping levels. We discuss how all the observed features can be naturally ascribed to the development of a doping-induced CDW phase assisted by disorder-induced strong localization, which may be in competition with the otherwise expected superconducting phase.

A. Electric-field-driven ion intercalation
Freshly-cleaved 2H-MoS 2 crystals (SPI supplies; typical size 2 mm×1 mm×50 µm) were electrically contacted via drop-casted spots of silver paste (RS Components) and immersed in a Duran crucible filled with 1-ethyl-3methylimidazolium tetrafluoroborate ionic liquid (EMIM-BF 4 , Sigma Aldrich) together with a platinum (Pt) counter electrode. The setup is sketched in Figure 1a. The gate voltage V G between the Pt electrode and the crystal was applied by an Agilent B2961 power source, at room temperature and in ambient atmosphere. The resistivity ρ was monitored in situ by the standard four-wire method. A constant current I DS ≈ 1 µA, flowing between the drain (D) and source (S) contacts, was supplied by a Keithley 220 current source, and the longitudinal voltage drop V xx between the inner voltage contacts was measured by a HP3457 multimeter. The resistivity was then determined as ρ = V xx I −1 DS t w l −1 , where t and w are the sample thickness and width, and l is the distance between the inner voltage contacts. Common-mode offsets were removed using the current-reversal method. Open-circuit (OC) conditions were obtained by physically disconnecting the gate electrode from the power source. Ex-situ characterizations were performed after extracting the intercalated MoS 2 samples from the cell and rinsing them thoroughly with acetone and ethanol to remove ionic-liquid residues. The samples were stored in standard desiccators either in argon atmosphere or under low vacuum to avoid moisture contamination.

B. Vibrational spectroscopies
Both Raman and Fourier-transform infrared (FT-IR) spectra were measured on freshly-cleaved surfaces in ambient conditions. Raman spectra were acquired using a Renishaw InVia H43662 micro-Raman spectrometer. All spectra were acquired using an excitation wavelength of 514 nm, a laser power < 1 mW focused through a 100X objective, an exposure time of 20 s, and 50 accumulations. FT-IR spectra were measured using a Thermoscientific Nicolet 5700 spectrometer in attenuated total reflectance (ATR) mode (Thermoscientific Smartorbit module) in the range between 900 and 3500 cm −1 . The spectrum of the EMIM-BF 4 ionic liquid was measured by mixing the liquid with potassium bromide. The determination of the peak areas in the bands centered around ∼ 1570 and ∼ 1045 cm −1 in both the liquid and the intercalated MoS 2 samples was performed by fitting the experimental spectra in the two regions to one and four Voigt peak functions respectively, in accordance with the peak assignments.

C. Nano-infrared microscopy
Scanning probe microscopy images were acquired on freshly-cleaved surfaces in ambient conditions by means of a Bruker Anasys nanoIR3-s atomic force microscope (AFM) combined with a multichip tunable quantum cascade laser (QCL; MIRcat-QT Daylight Solutions) covering the IR spectral ranges 1900-1350 cm −1 and 1150-900 cm −1 , using commercial gold-coated silicon tips (Bruker PR-EX-TnIR-A). AFM topography maps were acquired in tapping mode, and the IR absorption maps were simultaneously recorded via heterodyne detection as implemented in the "Tapping AFM-IR" mode of the built-in Analysis Studio software. All maps were acquired with ∼ 60 % of the average QCL laser power (0.5 W), a duty cycle of 15%, and a pulse rate ∼ 300 kHz set as to allow for heterodyne detection of the IR signal amplitude. Absorption spectra were collected with a spectral resolution of 2 cm −1 and 64 co-averages, and normalised to the QCL emission profile at the same laser power. Gwyddion software [68] was employed for image analysis.

D. Kelvin-probe force microscopy
Maps of the surface potential were acquired on freshlycleaved surfaces in ambient conditions by means of a Bruker Innova AFM equipped with a surface potential imaging addon, using standard platinum-iridium-covered silicon tips (Bruker SCM-PIT-V2). Pristine and intercalated samples were mounted simultaneously on the same metallic sample holder and fixed with electrically-conducting adhesive tape. AFM topography maps were acquired in tapping mode during the forward scans, whereas Kelvin-probe force microscopy (KPFM) measurements were performed in lift mode during the backward scan by grounding the sample and biasing the KPFM tip. The 20 × 20 µm 2 maps shown in the following were acquired using a scan rate of 0.3 Hz, a lift height of 50 nm, and a tip AC bias voltage of 3 V. Gwyddion software [68] was employed for image analysis.

E. Resistivity measurements
The measurements of the resistivity as a function of temperature were performed in the high-vacuum chamber of a Cryomech pulse-tube cryocooler with a base temperature of ≈ 2.8 K. The resistivity was determined with the same method employed in the room-temperature gating runs, except that I DS was supplied by an Agilent B2912 sourcemeasure unit, and V xx was measured by an Agilent 34420 nanovoltmeter. In each cooling-warming cycle, only the data recorded during the slower, quasi-static warming to room temperature were considered.

III. RESULTS
The ionic gating technique is widely employed to control the electronic properties of various classes of materials. The principle of the technique is that the material under study constitutes the active channel of an electrochemical transistor, and is thus separated from a gate counterelectrode by an electrolyte such as an ionic liquid. When polarized by the application of a finite V G , the electrolyte allows the mobile ions to accumulate at the sample surface and induce an equal and opposite sheet carrier density in the channel to maintain charge neutrality, building up the electric double layer (EDL) which acts as a nanoscale capacitor with ultra-high capacitance [18][19][20]. The intense electric field in the EDL can then be directly exploited to electrostatically tune this surface charge density as in a conventional field-effect transistor [21][22][23][24][25][26][27][28][29][30][31][32][33], or to activate more complex electrochemical effects including bulk ion intercalation [33][34][35][36][37][38][39][40][41][42][43][44][45]. Here, we make use of the ionic gating setup sketched in Fig. 1a to intercalate positively-charged organic ions in bulk 2H-MoS 2 crystals and induce a semi-permanent electron doping. Figure 1b shows the resistivity ρ of a representative MoS 2 crystal during the gating process at room temperature. At V G = 0, a finite ρ ≈ 18 Ω cm is measured due to the intrinsic n-doping arising from the sulphur vacancies ubiquitous in MoS 2 crystals [69]. As V G is swept from 0 to positive values, ρ remains mostly unaffected as long as V G +1 V (region 1 in Figure 1b) and then starts to decrease, reaching a plateau around ρ ≈ 5.5 Ω cm in the V G range between ∼ +2.75 and ∼ +3.25 V (region 2 in Figure 1b). Further increasing V G to +3.5 V then triggers a sharp decrease in ρ of about 3 orders of magnitude. If the gate voltage is then kept fixed at V G = +3.5 V the resistivity saturates rather quickly ( 10 min, region 3 in Figure 1b) to ρ ≈ 3 mΩ cm and shows no further change (region 4 in Figure 1b). This gate-induced change in ρ is partially volatile: when the potential difference between the Pt counter electrode and the sample is removed by setting the cell in the OC conditions (region 5 in Figure 1b), ρ rapidly increases, saturating to an intermediate value ρ ≈ 150 mΩ cm after about 15 minutes; however, smaller values of ρ can be stabilized by extracting the MoS 2 crystal from the cell before the full OC saturation has taken place. Figure 1c shows the ρ of gated MoS 2 crystals measured ex situ as a function of the release time, i.e. the time span in which each crystal has been left in the ionic liquid in OC conditions. While the lowest-resistivity state (ρ ≈ 3 mΩ cm) observed in region 4 goes lost as soon as the gate circuit is opened, permanent modifications of the intercalated crystals leading to intermediate ρ values ranging from ≈ 25 mΩ cm up to 150 mΩ cm (and monotonically depending on the release time) can thus be obtained. With the exception of the electric transport measurements, all characterizations have been performed on samples at zero (nominal) release time.
Since 2H-MoS 2 is prone to undergo structural transitions upon heavy electron doping [54,56] and intercalation by different ionic species [53,55,[57][58][59][60], the structural phase of both the pristine and the intercalated crystals was assessed via Raman spectroscopy. As shown in Figure 2a, the Raman spectrum of a pristine crystal (solid black line) exhibits the two typical modes of MoS 2 in the 2H crystal structure, the E 1 2g mode at ∼ 384 cm −1 and the A 1g mode at ∼ 409 cm −1 , which correspond to in-plane and out-of-plane vibrations of Mo and S atoms respectively [11,70]. The same two peaks are also observed in the Raman spectrum of intercalated MoS 2 (solid red line), with no sign of the emergence of the peaks at 158, 218 and 334 cm −1 which are associated to the 1T /1T ' MoS 2 polytypes [57][58][59][60]. This indicates that our MoS 2 samples retain their 2H crystal structure upon EMIM intercalation. Furthermore, Lorentzian fitting to the experimental spectra (inset to Figure 2a)   cates that the A 1g mode undergoes a tiny ≈ 0.3 cm −1 redshift in the intercalated sample, whereas the E 1 2g mode remains unchanged. Both features are consistent with the behavior of MoS 2 nanolayers upon gate-induced electron doping [13,71,72], even though a quantitative comparison is prevented by the strong dependence of the effect on the sample thickness.
The successful incorporation of EMIM ions in the MoS 2 lattice was confirmed by FT-IR measurements. As shown in Figure 2b, the infrared spectrum of the pristine MoS 2 crystal (solid black line) is mostly featureless in the investigated frequency ranges, since the sharp infrared peaks typical of the 2H phase emerge at wavenumbers below ∼ 500 cm −1 [73]. Conversely, the spectrum of the intercalated MoS 2 sample (solid red line) exhibits several bands, that are actually typical of organic compounds. Indeed, they fall in the C-H stretching region (between 3200 and 2800 cm −1 ), and in the fingerprint region (between 1800 and 900 cm −1 ). Notably, the same bands are also observed in the spectrum of the EMIM-BF 4 ionic liquid (solid blue line). In accordance with the literature, we assign the structure between 3200 and 3050 cm −1 to the C-H stretching modes involving unsaturated C atoms in the EMIM ring [74][75][76]; the band between 3000 and 2850 cm −1 to the aliphatic C-H modes of the CH 3 group bonded to the EMIM ring, the CH 3 terminal group of the ethyl chain, and the CH 2 group [74][75][76]; the peak around 1572 cm −1 to the EMIMring breathing [74][75][76]; the band around 1450 cm −1 to the (CH 2 )/CCH HCH bending modes [74][75][76]; the peak around 1170 cm −1 to the C-C aliphatic stretching modes [74][75][76]. All these bands are associated purely to the EMIM cation. Conversely, the broad band between 1150 and 900 cm −1 comprises at least four peaks arising both from the cation and the BF 4 anion. Specifically, we assign the peak around 1116 cm −1 to ring HCCH symmetric bending [75,76]; the peak around 1085 cm −1 to C-N stretching [75]; the peak around 1050 cm −1 to a combination of BF 4 asymmetric stretching and C-N symmetric stretching [76]; and the peak around 1021 cm −1 to BF 4 asymmetric stretching [75,76]. As such, only this last peak can be uniquely associated to the presence of the BF 4 anion.
The fact that all the main absorption bands of the EMIM-BF 4 ionic liquid are observed also in the infrared spectrum of the intercalated MoS 2 sample indicates that both the EMIM cation and the BF 4 anion are incorporated in the MoS 2 lattice. However, the relative intensity of the BF 4related bands with respect to that of the EMIM-related ones is much lower in the intercalated MoS 2 sample than in the ionic liquid. Specifically, the peak-area ratio between the BF 4 -stretching mode at ∼ 1021 cm −1 and the EMIMbreathing mode at ∼ 1572 cm −1 in the EMIM-BF 4 ionic liquid (A BF 4 /A EMIM ≈ 12) is about 10 times larger than that in the intercalated MoS 2 sample (A BF 4 /A EMIM ≈ 1.2). Considering that the ionic liquid has a 1:1 EMIM:BF 4 stoichiometry, this allows us to estimate that in the intercalated MoS 2 samples about one BF 4 anion every ∼ 10 EMIM cations is incorporated in the lattice as a contaminant. The fact that the MoS 2 lattice incorporates a much larger amount of cations than of anions confirms a net electron doping in the intercalated samples.
The spatial distribution of the intercalants was investigated by nanoscale infrared spectroscopy [77,78].  phase image (Figure 3b), which indicates that their mechanical properties differ significantly from those of the MoS 2 background. IR absorption spectra were then collected both on one of these "droplets" (point D) and on the exposed surface (point S). As shown in Figure 3b, both spectra display the peaks, specific of the EMIM cation, at about 1572 cm −1 , 1450 cm −1 , and 1170 cm −1 , as well as the broad band between 1150 and 900 cm −1 , that also includes the BF 4 -related peak at ∼ 1021 cm −1 . However, in the spectrum acquired on the droplet (solid blue line) this latter band is much more pronounced than the EMIM-related peaks, as it happens in the FT-IR spectrum of the pure ionic liquid shown in Figure 2b. Conversely, in the spectrum acquired on the exposed MoS 2 surface, the BF 4 -related band is significantly less intense with respect to the EMIM-related peaks: in particular, the spectral weight of the BF 4 -stretching peak at ∼ 1021 cm −1 is strongly suppressed. Combined with the phase contrast, these results strongly suggest that the droplets are due to ionic-liquid residues embedded within the MoS 2 planes (and exposed to the surface after cleaving), and their presence is the main source of BF 4 anions detected by FT-IR. Finally, the spatial distribution of the intercalated EMIM cations was studied by acquiring an IR absorption map in correspondence of the EMIM-breathing mode at ∼ 1572 cm −1 : as shown in Figure 3d, a finite IR absorption at this wavenumber can be observed throughout the entire surface, and it is notably minimum in the areas covered by the ionic-liquid residues; however, a longrange modulation in the IR absorption over a length scale ∼ 10 µm can be distinguished, which indicates that EMIM incorporation in the MoS 2 lattice is inhomogeneous and develops local minima and maxima over a comparable length scale.
The impact of the EMIM intercalation on the electronic structure of MoS 2 was assessed by means of Kelvin-probe force microscopy [79][80][81] on a pristine MoS 2 crystal and on the same intercalated crystal characterized via nano-IR. Figure 4a and b show the AFM topography map and surface potential distribution of the pristine MoS 2 respectively. As expected for the cleaved surface of a layered material, the topography is very flat (S q ≈ 1.5 nm), atomic steps and terraces can be distinguished, and the surface potential is featureless and nearly homogeneous across the entire area. The topography map of intercalated MoS 2 (Figure 4c) again exhibits a much larger background corrugation (S q ≈ 15 nm) and the localized features associated to ionic-liquid contaminants. More importantly, the surface potential of intercalated MoS 2 (Figure 4d  content detected via nanoIR. As clearly shown by the direct comparison of the histograms of the two images reported in Figure 4e, the surface potential of intercalated MoS 2 is also shifted to lower values with respect to that of pristine MoS 2 across the entire area. This reduction in surface potential is a direct proof of electron doping due to the insertion of the EMIM cations in the MoS 2 lattice, since it signifies that in the intercalated MoS 2 the Fermi level has been shifted upwards in energy and closer to the vacuum level. The average shift in Fermi energy ∆E F between the two systems can then be estimated simply as: (1) where e is the elementary charge, 〈S.P.〉 pristine = +12 mV and 〈S.P.〉 intercalated = −163 mV are the average values of the surface potential in pristine and intercalated MoS 2 respectively. From this shift in the Fermi energy the average density of electrons n e doped into the conduction band of MoS 2 can also be determined, since the Fermi level of pristine MoS 2 is known to be pinned to the bottom of the conduction band by the presence of sulphur vacancies in the crystal [69]. In particular, assuming a parabolic three-dimensional dispersion of the density of states of bulk MoS 2 : with a spin degeneracy ν s = 2, a valley degeneracy ν v = 6 [47,82], and an effective mass m * = 0.6m e [82], one obtains n e ∼ 9 × 10 20 cm −3 . If every intercalated EMIM cation were to donate one electron to the MoS 2 conduction band, this would correspond to about one EMIM cation every ∼ 20 MoS 2 formula units, or a nominal stoichiometry of EMIM 0.05 MoS 2 .
Having established the mechanism and amount of electron doping provided by the EMIM intercalation in the MoS 2 crystals, we now turn to consider how this affects their electric transport properties as a function of temperature T . Figure 5a shows the T -dependence of ρ, divided by its value at T = 300 K (shown in Figure 1c), down to T ∼ 3 K and for a series of MoS 2 crystals with increasing release times. The pristine MoS 2 crystal shows the exponentially-increasing ρ with decreasing T typical of insulators where conduction occurs via hopping processes in the localized states in the band tails, as expected [48,69]. The T -dependence of ρ in the intercalated samples is much less steep, and its values at low T are several orders of magnitude lower than that of pristine MoS 2 , consistent with a metallization induced by the electron doping. A closer inspection of the curves for the intercalated samples ( Figure 5b) reveals a complex dependency of ρ on both T and release time. The sample at zero release time exhibits an overall metallic behavior, in the sense that its ρ decreases on going from 300 K to 3 K; however, this T -dependence is not monotonic, and a clear anomaly can be observed in the form of a broad hump appearing between 50 and 150 K. This behavior is best visualized by plotting the first derivative of ρ with respect to T , dρ/d T . As shown in Figure 5c, this anomaly gives rise to  a hump-dip structure in the T -dependence of dρ/d T , with a broad hump at lower T and a sharp dip at higher T , the latter being highlighted by the dashed lines. Additionally, a small resistance upturn is observed for T 20 K, indicating that at low T some degree of carrier localization is present in the sample.
We track the evolution of these two main features -degree of metallicity and resistivity anomaly -on increasing the release time, by defining two distinct figures of merit: for the metallicity, the residual resistivity ratio RRR, defined as usual as ρ(300 K)/ρ(3 K); for the resistivity anomaly, the temperature T * where the dip in dρ/d T is located. Figure 5d shows that, upon increasing the release time (i.e., decreasing the electron doping), T * (violet circles, left scale) decreases in a nearly monotonic fashion from ∼ 155 K at release times below ∼ 2 min, to ∼ 115 K at release times 15 min. At the same time, however, RRR (green diamonds, right scale) exhibits a clearly non-monotonic trend: after an initial increase from ∼ 2.7 to ∼ 5.7 for release times 4 min, it strongly decreases down to ∼ 0.58 at release times 15 min, i.e. when spontaneous deintercalation of the EMIM ions stops. This indicates that the intercalation stage obtained with a release time of 4 min achieves a maximum degree of metallicity, and doping the sample away from it shifts the system towards an insulating behavior.

IV. DISCUSSION
In our EMIM-intercalated MoS 2 crystals, increasing the release time reduces the doping and therefore brings the samples back toward the MoS 2 pristine state of a band insulator. On the other side, reducing the release time increases the doping level and brings the samples closer to an insulating state that, however, must have a different origin, since increasing the concentration of EMIM dopants should in principle provide a larger density of free electrons to the conduction band. Such insulating state must therefore be due to some kind of localization.
A first source of localization could be the randomness caused by intercalation and progressive introduction of disorder in the form of extra scattering centers due to the presence of the ions themselves [36,[83][84][85][86][87]: this type of gate-driven re-entrant transition (from a band insulator, to a metal, and finally to an Anderson insulator) has already been reported in other ion-gated TMDs [88][89][90] and oxides [42]. A second source of localization could instead be associated to the same phenomenon giving rise to the resistivity anomaly, which is often found among TMDs and represents one of the typical signatures of the onset of CDW phases [1]. These consist in a periodic modulation of the charge carrier density accompanied by a distortion of the underlying crystal lattice [91], usually resulting from strong electron-phonon coupling [92] and/or Fermi-surface nesting [93]. CDW phases are very susceptible to charge carrier doping: in most TMD compounds, CDWs are suppressed by increasing doping in favor of the development of superconducting phases [28,34,[94][95][96], but in some specific compounds the CDW order is instead promoted by doping [26,37,52]. Density functional theory calculations have predicted 2H-MoS 2 to fall into the latter category, both in the bulk [53,55] and in the monolayer [54,56] form. Recent experimental evidences have confirmed that, indeed, the onset of a CDW order in MoS 2 bulk crystals and bulk-like flakes is promoted by Li-doping [37] and K-doping [52], as well as by pressure [97]. In this context, T * would mark the critical temperature for the CDW phase transition in EMIMintercalated MoS 2 , placing it in an intermediate range between those of Li x MoS 2 (≈ 230 K) and K x MoS 2 (≈ 85 K).
The onset of CDW order can naturally explain why resistivity anomalies emerge in our EMIM-intercalated MoS 2 crystals, and can also provide an additional possible mechanism for the observed incipient localization at high doping, i.e. the opening of a (possibily partial) gap on the Fermi surface associated to the CDW [52,55]. It is interesting, at this stage, to understand the reasons why superconductivity is not observed in MoS 2 , at least down to 2.8 K, upon EMIM intercalation. The fact is that both Li-doping and K-doping instead do lead to the appearance of a superconducting phase, although with some differences. In K-intercalated MoS 2 , where the CDW order sets in at a low T * , bulk superconductivity occurs with a critical temperature T c ≈ 7 K [51, 52] throughout most of its phase diagram with little change in T c as a function of the K content (between K 0.13 MoS 2 and K 1.5 MoS 2 ) [51]. In Li-intercalated MoS 2 , where T * is high, superconductivity (with T c ≈ 3.7 K) is observed only in samples showing no sign of CDW-related resistivity anomalies [36]. In our EMIM-intercalated crystals, in which the CDW order is ubiquitous, but sets in at intermediate temperatures T * , no sign of superconductivity has been detected. One could object that, maybe, the range of doping contents of our samples is not the right one to observe superconductivity. This is certainly possible, but since the EMIM doping distribution is largely inhomogeneous, the system actually experiences a wide range of local electron doping levels, from very small ones to others, much larger than the high-est average doping content calculated above. This means that, if there were a sufficiently wide range of doping for superconductivity to develop, at least a fraction of the volume should fall in that range. Even if this fraction were not sufficient to ensure superconducting percolation, a reduction in the resistivity should be observed. The fact that this does not happen suggests that the superconducting phase either does not develop at all or, if it does, has a very low T c (smaller than 2.8 K).
This whole picture seems to indicate that the CDW order is competing with superconductivity in MoS 2 , so that the superconducting phase appears with higher T c only when the CDW is suppressed or at least T * is lower. Alternatively, the lack of superconductivity in EMIM-intercalated MoS 2 could be due to the large size of the organic ions, which strongly increases the interlayer spacing: this would make the intercalated system akin to a stack of electronically-decoupled MoS 2 monolayers, where superconductivity occurs with a much lower T c 2 K and is much more easily suppressed by disorder [98,99]. Such decoupling was reported in the case of another TMD compound when intercalated via organic molecules, SnSe 2 , where however the electronic decoupling was shown to increase the superconducting T c [100,101]. Overall, which of the depicted scenarios -or combination thereof -is at the origin of the anomalous metallic transport exhibited by our EMIM-intercalated MoS 2 crystals is left as an open question, making future investigations highly desirable. Electric transport measurements below ∼ 3 K would of course be needed to ascertain whether superconductivity develops at lower temperatures than those achievable in the setup used in the current study. Moreover, lowtemperature structure-sensitive techniques such as Raman spectroscopy, X-ray diffraction, and scanning tunnelling microscopy would be highly beneficial to provide direct evidence for CDW ordering in the system, as well as to track its evolution as a function of temperature and doping.

V. CONCLUSIONS
In summary, we synthesized organic ion-intercalated molybdenum disulphide crystals via the ionic gating method with EMIM-BF 4 ionic liquid. We showed that the intercalation is partially non-volatile and can be tuned by allowing spontaneous de-intercalation to occur when the crystals are left in the electrochemical cell in opencircuit conditions. By combining vibrational spectroscopies, Kelvin-probe force microscopy, and nano-infrared microscopy, we demonstrated that the EMIM ions act as electron donors and are intercalated in the entire sample volume, with their density exhibiting fluctuations over a length scale of the order of ∼ 10 µm. We characterized the temperature dependence of the resistivity in the intercalated crystals, showing how increasing EMIM intercalation drives the system across a re-entrant insulator-to-metal transition, i.e. from the pristine band insulator into an anomalous metallic state, and eventually towards an incipient second insulating state. We further revealed the appear-ance of a doping-induced hump in the resistivity around ∼ 150 K, which we attribute to the development of a chargedensity wave phase on the basis of the established literature on Li-doped and K-doped MoS 2 . A qualitative comparison between these three systems, in particular focused on the relation between the appearance of superconductivity and the onset of the charge-density-wave phase, is suggestive of a competition between these two quantum phases in electron-doped MoS 2 .