Spontaneous chemical functionalization via coordination of Au single atoms on monolayer MoS2

See allHide authors and affiliations

Science Advances  02 Dec 2020:
Vol. 6, no. 49, eabc9308
DOI: 10.1126/sciadv.abc9308


Surface functionalization of metallic and semiconducting 2D transition metal dichalcogenides (TMDs) have mostly relied on physi- and chemi-sorption at defect sites, which can diminish the potential applications of the decorated 2D materials, as structural defects can have substantial drawbacks on the electronic and optoelectronic characteristics. Here, we demonstrate a spontaneous defect-free functionalization method consisting of attaching Au single atoms to monolayers of semiconducting MoS2(1H) via S-Au-Cl coordination complexes. This strategy offers an effective and controllable approach for tuning the Fermi level and excitation spectra of MoS2 via p-type doping and enhancing the thermal boundary conductance of monolayer MoS2, thus promoting heat dissipation. The coordination-based method offers an effective and damage-free route of functionalizing TMDs and can be applied to other metals and used in single-atom catalysis, quantum information devices, optoelectronics, and enhanced sensing.


Chemical surface functionalization has been used to tune the electronic, optical, and catalytic properties of two-dimensional (2D) transition metal dichalcogenides (TMDs) (13). Chemical functionalization commonly relies on the presence of lattice defects and physisorption methods to circumvent the chemical inertness of pristine semiconducting TMDs (4). Unfortunately, such routes inevitably modify the surface characteristics as well as the optical, thermal, and transport properties of the atomically thin layers (59). Most of the proposed MoS2 functionalization routes are based on liquid exfoliated flakes, which have reduced and more reactive lateral sizes. In addition, this exfoliation route is not ideal for scaling up processes in the electronics and optoelectronics industries (10). For example, Voiry et al. (6) reported covalent functionalization using chemically exfoliated MoS2, WS2, and MoSe2 with organohalides and its impact on their photoluminescence (PL) emission. Although this functionalization did not entirely rely on defect engineering, it is restricted to TMDs in the 1T phase, and the exfoliation procedure involves harsh chemical treatments, such as butyllithium. Efforts have also been made to functionalize chemical vapor deposition (CVD)–grown TMDs. For example, Ding et al. (2) reported that thiol functionalization of CVD-grown MoS2 led to an increase in the PL intensity due to the passivation of S vacancies by the ligand, thus relying on the presence of inherent defects to achieve functionalization. However, in this case, the thiol does not provide perfect healing to the crystalline structure; thus, the charge mobility of a pristine crystal is not reestablished/recovered. To the best of our knowledge, chemically bonding single atoms, molecules, or clusters to the surface of pristine crystalline TMD monolayers remains a challenge. For metal monochalcogenides, Lei et al. (11) proposed the formation of coordination complexes on the surface InSe and showed that the presence of nonbonding electron pairs can be exploited by applying Lewis acid-base concepts to form coordination bonds with metal ions. As a result, the electronic properties of InSe were tailored and p-type doping was observed through this functionalization approach.

In this work, we developed a route to functionalize monolayers of MoS2 with individual Au atoms via the formation of S-Au-Cl coordination complexes ([Au(MoS2)Clx]) on the TMD surface. The [Au(MoS2)Clx] coordination complexes have been synthesized by taking advantage of the lone pair electrons of the S atoms present in the MoS2 lattice. Unlike substitutional doping or defect passivation, the formation of the metal-MoS2 coordination complexes does not rely on the presence of additional defects on the 1H-MoS2 but induces significant changes in the optical, thermal, and electrical properties. Density functional theory (DFT) calculations confirmed that single AuClx complex would bond to S atoms, rather than randomly absorbed on top of Mo atoms or over the vacant spots of the hexagonal lattice, and scanning transmission electron microscopy (STEM) imaging revealed the presence of Au single atoms fixed on top of S atoms. During the Au functionalization, the Au3+ precursor was spontaneously reduced to Au1+ on the TMD surface before forming the Au─S coordination bonds. This Au─S bond results in an efficient transfer of electrons from the MoS2 valence band with high contribution of S 3p orbitals to the Au 5d and 6s orbitals, thus tuning the Fermi level of MoS2 monolayers. By treating MoS2 with different Au concentrations, we found an effective way to tune the Fermi level of MoS2 via controlled p-type doping, as measured in field-effect transistors (FETs). In addition, Au functionalization leads to high exciton-to-trion ratios in the PL spectra of monolayer MoS2. The stability of the S-Au-Cl coordination sphere leads to the existence of [Au(MoS2)Clx] complexes that results in a major enhancement on the thermal boundary conductance across MoS2 monolayers. This coordination method could also be used to synthesize Ag single atoms on MoS2. The isolation and anchorage of single noble atoms via a solution-phase chemical approach pave the way for large-scale manipulation of single atoms, leading to high-performance catalysis and quantum information applications (12).


Synthesis of [Au(MoS2)Clx] complexes: Single Au atoms in motion

MoS2 monolayers were synthesized by CVD on Si/SiO2 substrates using powder precursors (see Materials and Methods for details). Optical microscopy, Raman spectroscopy, and PL spectroscopy were carried out to characterize the MoS2 flakes (see fig. S1). The coordination reaction between the MoS2 monolayers with HAuCl4 was performed by dipping the as-grown MoS2 into an ethanol solution of HAuCl4 exhibiting various concentrations for 10 min, following washing and drying. High-resolution STEM (HR-STEM) imaging using a high-angle annular dark-field (ADF) detector was conducted to study the atomic structure of the Au atoms on MoS2. The ADF intensity changes, depending on the Z-number of atoms (~Z1.6–1.9) (13); thus, the higher Z-number Au atoms stand out in the MoS2 lattice. As shown in Fig. 1A, bright single Au atoms on top of the MoS2 lattice can be observed. These Au atoms remain isolated, rather than aggregated. The largest Au aggregate found by STEM is shown in the inset of Fig. 1A, where four Au atoms occupy neighboring S sites. Note that we consider this as a planar aggregate of Au single atoms, rather than a regular Au cluster, as the atomic spacing corresponds to that of the S in the MoS2 lattice in a planar view, rather than the predicted 1.34 Å interatomic distance between Au atoms in planar clusters (14, 15). Atomic resolution STEM-ADF images confirmed the presence of Au single atoms and their preferred position on the MoS2 monolayer, which is directly on top of the S sites [see Fig. 1 (A and B)]. The DFT-simulated STEM-ADF image (Fig. 1C) was obtained on the basis of the model shown in Fig. 1A (inset) and exhibits an excellent agreement with the experimental image. Single-atom Ag could also be synthesized via AgNO3 ethanol solution displaying similar Ag-MoS2 structures (see fig. S2). An ADF intensity line scan was also performed across the single Au atom shown in Fig. 1 (B and C) (vertical gray line in the figures). In the experimentally acquired ADF intensity line profile (Fig. 1D), besides observing the ADF intensity peak at the Mo site and 2S sites, an intense peak corresponding to the Au-2S site was found matching the simulated ADF intensity at the Au-2S site (Fig. 1E and fig. S3).We believe that Cl atoms are bonded to the Au single ions to complete the Au coordination sphere and balance the charge on the Au atoms in the form of [Au(MoS2)Cl] and [Au(MoS2)Cl3] complexes (see below for details). However, STEM simulations conducted for [Au(MoS2)], [Au(MoS2)Cl], and [Au(MoS2)Cl3] structures (fig. S3) indicate the absence of Cl atoms on our experimental STEM images. We believe that this is due to the electron beam (e-beam) irradiation at 80 keV that may be displacing Cl atoms during imaging. HR-STEM was also used to observe changes in the position of individual Au atoms on the MoS2 surface over time. Figure 1F shows selected frames of a video (see movie S1) showing the movement of a single Au atom on MoS2 under e-beam irradiation. The beam provides energy to the Au atom, leading to a dynamic atom displacement. It can be observed that the Au atom is capable of hopping from one S coordination site to another. A clear trajectory (gray curve) was determined for the atom movement. For the first 100 s of the referred video, the tracked Au atom stayed directly above the same S atom, and it then jumped toward a neighboring S site, occupying that for several more seconds. The displacement distance of each hop is shown in Fig. 1G. The Au atom moves to the nearest S atom most of the time but can cover a distance of 2 nm during one frame (10 s per frame). The displacement distance covered by the Au atom shows that the atom stays almost static when it is on top of S, while it moves much more actively when it is on top of Mo. We could not find a Au atom that would stay above the Mo atom or at any position other than above a S site for several seconds; thus, we conclude that those positions are only transitory and not stable configurations. This observation indicates that stronger interactions between the Au and S atoms make S sites energetically favorable. The movement of single Au atoms, if deterministically manipulated, could be further exploited in future quantum information devices (12).

Fig. 1 Synthesis of Au single atoms on monolayer MoS2 flakes.

(A) HR-STEM images of [Au(MoS2)Clx] featuring the single Au atoms on S sites. Insets: The structure model of Au-MoS2 used for STEM image simulation (top right) and the higher-magnification STEM image of a Au aggregate found on MoS2 (bottom left).(B and D) High-angle ADF (HAADF) images and Z-contrast line scan [along the vertical gray line in (B)] of one Au single atom on MoS2, revealing the position of the Au atom to be directly on top of the S atom. a.u., arbitrary units. (C and E) TEM simulation of one Au atom directly on top of one S atom in MoS2 lattice and Z-contrast line scan of the simulation image. The simulated TEM image and line scan match exactly with experimental data, confirming the atomic structure of Au single atoms directly bonded on S atoms on the MoS2 surface. (F) Time-lapse images of Au single atoms on MoS2. The Au single atoms can move on the MoS2 plane under the e-beam irradiation during STEM characterization. The full video can be viewed in the Supplementary Materials. The gray curves show the trail of one Au atom moving on MoS2. (G) The statistical analysis of the position of one Au atom and displacement distance between each movement.

[Au(MoS2)Clx] complex coordination sphere

To further understand the chemical nature of the functionalization and to characterize the Au oxidation state, x-ray photoelectron spectroscopy (XPS) was carried out. Figure 2A shows the Au 4f core-level spectrum of the [Au(MoS2)Clx] complex (green curve). The spectrum can be well fitted by two sets of doublet peaks (red and blue curves) that can be assigned to Au3+, from the HAuCl4 precursor, and Au1+, indicating spontaneous Au reduction on MoS2. Specifically, the Au3+ signature corresponds to the 4f5/2 (4f7/2) peak at 91.0 (87.4) eV, while the Au1+ signature relates to the 4f5/2 (4f7/2) peak at 88.2 (84.5) eV, which are comparable to other studies reporting Au ions with S-containing ligands (16, 17). Huang et al. (18) have reported that gold can be spontaneously reduced from Au3+ to Au0 nanoparticles in a solution of chemically exfoliated MoS2. In this work, we demonstrate that Au3+ can be reduced into Au1+ in the form of single atoms due to the formation of a coordination bond to the nondefective surface of CVD-grown MoS2 monolayers. On the basis of the standard reduction potentials for the Au3+/Au1+ pair (0.926 V) and MoS2 (−0.090 V), electrons from the MoS2 monolayer can be donated to Au3+ so that Au3+ is reduced to Au1+ (19). In addition, according to Pearson, in hard-soft acid and base theory (20, 21), soft ligands, such as S, stabilize linear Au1+ complexes (22).

Fig. 2 Spectroscopic studies of the [Au(MoS2)Clx] complex.

(A) XPS spectra of Au 4f orbit on Au-MoS2. The green curve is the original spectrum after C1s(248.8 eV) calibration. The blue and red curves are fittings of the original data, which show the presence of Au3+ and Au1+. We believe that the presence of Au3+comes from the unreacted HAuCl4 precursor, whereas the Au1+ from the reduction reaction with MoS2. The purple curve is the sum of the blue and red curves. (B) S 2p orbit of pristine and Au functionalized MoS2. The green curve is the original spectrum after C1s(248.8 eV) calibration. The blue and red curves are fittings that show the S 2p1/2 and 2p3/2 peaks, respectively. The S 2p orbitals shift to higher energy because of the loss of outer electrons to Au─S bonds. (C) PL spectrum of pristine and functionalized MoS2 monolayers. The A exciton is deconvoluted into the trion (X) (blue curve) and exciton (X) (red curve) peaks through Lorentzian functions. After doping, the trion intensity decreased, while the exciton intensity increased, which corresponds well with the p-type doping effect of the Au1+. (D) Exciton-to-trion intensity ratio of pristine and functionalized MoS2 with different precursor Au concentrations.

The formation of Au─S bonds was also confirmed by the S 2p core-level spectra shown in Fig. 2B. In pristine MoS2 monolayers, the S shows a doublet peak corresponding to a 2p1/2 peak at 163.0 eV and a 2p3/2 peak at 161.9 eV. After functionalization, the doublet peak shifted to a higher binding energy (0.8 eV shift), thus confirming the formation of Au─S bonds. The donation of lone pair electrons from S to Au atoms results in a higher electrostatic attraction on the inner shells of S, causing XPS peaks with higher binding energy (17). The shift further confirms the formation of Au─S bonds. Therefore, the formation of Au1+─S bonds contributes to the stabilization, partial reduction, and formation of single Au ions on the MoS2 surface. In addition, the absence of water inhibits the Au1+ disproportionation reaction that results in the formation of Au0 and Au3+, which, in turn, prevents the formation of gold clusters or nanoparticles (22). In addition, XPS data reveal the presence of Cl with a Au:Cl ratio of 1:3, as shown in fig. S5, which could be due to the presence of unreacted Au precursor HAuCl4 and/or from AuClx complexes coordinated on MoS2 surface (see Discussion). Similar peak shifting of S 2p orbitals and additional peak splitting were observed for Ag-MoS2 functionalizations demonstrating the formation of Ag─S bonds (see fig. S4) (23).

To understand the effect of the single AuClx complex functionalization on the optical properties of monolayer MoS2, we carried out PL spectroscopy measurements. As seen in Fig. 2C, PL studies reveal the presence of a pristine MoS2 monolayer characteristic signature, due to its direct band gap at the K(K′) point. The peak shown in Fig. 2C can be deconvoluted, via fitting with two Lorentzians, into a neutral exciton (X; red) at 1.83 eV and a negative trion (X; blue) located at 1.79 eV. It can also be seen that the X trion was markedly quenched, while the neutral exciton X intensity was enhanced after the Au functionalization. The formation of the coordination bonds between the S and Au atoms results in MoS2 electrons being transferred to the Au atoms, thus leading to a decrease in the X trion intensity and an increase in the neutral exciton intensity in MoS2 (3). Therefore, we believe that the recorded changes in the PL after functionalization are caused by the p-type doping of MoS2 via AuClx coordination. The p-type doping level of MoS2 via AuClx coordination can be carefully controlled by exposing the MoS2 monolayers to different concentrations of HAuCl4 solutions. Figure 2D shows the exciton-to-trion intensity ratio as a function of the Au concentration. It can be clearly seen that this ratio increases with the Au concentration, offering an effective approach to sharpen the PL spectra of MoS2. A more detailed discussion on the p-type doping effect is available below (Fig. 3). Raman spectra of monolayer MoS2 before and after the Au functionalization exhibit the characteristic in-plane E′ and out-of-plane A′1 modes, without noticeable frequency nor intensity differences (fig. S1). As shown in table S1, no significant changes can be observed either in the Mo─S bond lengths and the Mo─S─Mo bond angles, for the optimized Au coordination in the 1H phase MoS2.

Fig. 3 Band structure studies of the [Au(MoS2)Clx] complex.

(A) DFT-optimized structure for the [Au(MoS2)Cl3] complex. Chlorine atoms were considered to complete the Au coordination sphere. (B) Charge transfer map of the [Au(MoS2)Cl3] complex showing a higher hole concentration along most of the MoS2 sheet (blue surface), and a pronounced electron concentration on the Mo and S atoms near the Au atom (red surface), commensurate with a p-type doping effect of the Au1+. (C) Electronic band structure of pristine MoS2 and of the [Au(MoS2)Cl3] complex, both calculated without considering the spin-orbit coupling. Electronic transitions are indicated by vertical arrows. In the [Au(MoS2)Cl3] complex, besides the pristine MoS2 excitons A and B, charge transfer transitions are expected between the MoS2 state (blue line) and the Au localized states (red line). (D) Room temperature absorbance spectra of the [Au(MoS2)Clx] complex (dotted line) and of pristine MoS2 (black). A slight blue shift of the A and B excitons and new features associated with charge transfer transitions at energies lower than 1.7 eV are observed, as expected from the electronic band structure of the [Au(MoS2)Cl3] complex shown in (C).

The formation of the Au─S bond was further investigated by DFT calculations, which were performed to identify the most energy-favorable structure for the [Au(MoS2)Clx] complex [shown in Fig. 3 (A and B)]. Cl atoms were added to the Au coordination sphere to balance the charge on the Au atoms and complete its coordination sphere, which is also confirmed by XPS (fig.S5). Thus, three different structures were evaluated, with an isolated Au atom ([Au(MoS2)]), with one Cl atom ([Au(MoS2)Cl]), and with three Cl atoms ([Au(MoS2)Cl3]). After the optimization simulations, three available coordination sites, known as the H (hollow), Mo, and 2S sites, were identified. In the H site, the Au is located at the center of the hexagon formed by the lattice projection on a plane (fig. S6); in the Mo site, the Au sits on top of one Mo atom and is bound to the three adjacent S atoms. Last, the 2S site, where the Au atom is positioned on top of a S atom, was found to be the most stable site by an energy difference of 550 meV (Fig. 3, A and B). Illustration of the H and S sites of Ag-MoS2 complexes can be found in fig. S7 and table S1, where the S site is found to be the most favorable for Ag single atoms. These results are in good agreement with the HR-STEM images depicted in Figs 1B and 2A and fig. S2, where Au and Ag locate preferentially over S atoms. Furthermore, to understand the spontaneous reduction of single AuClx complexes, the charge distribution of the three systems was calculated and analyzed using the Voronoi deformation density method. The Au atom charge in the [Au(MoS2)Cl3] and [Au(MoS2)Cl] optimized geometries were found to be +0.80 and 0.40 eV, respectively, indicating that the charge transfer from the S atoms leads to the spontaneous reduction of Au3+ to Au1+ during the Au─S bond formation. The [Au(MoS2)Cl3] projected charge transfer map shown in Fig. 3B indicates a higher electron concentration in the atoms near the Au─S bond. On the contrary, the MoS2 regions far from the AuCl3 exhibited a higher hole concentration, which can be understood as p-type of doping. For the [Au(MoS2)Cl], the projected charge transfer map shown in fig. S8 indicates a very small hole concentration on MoS2, whereas for Au-MoS2, we found a high electron concentration on the MoS2 and an n-type of doping, showing Cl atoms bound to Au to provide the p-type doping observed experimentally. The [Au(MoS2)Cl] and [Au(MoS2)Cl3] data are in good agreement with XPS, where the Au 4f core-level spectra confirm the Au1+ formation and the Cl 2p spectra indicate the presence of Cl on the surface. In addition, the exciton-to-trion intensity ratio changes observed in the PL spectra (Fig. 2, C and D) show a clear p-type doping of MoS2.

From the calculated electronic band structure, displayed in Fig. 3C, we notice that the MoS2 remained a direct gap semiconductor after functionalization. However, the bands near the Fermi level were drastically affected by the presence of single Au1+ ions coordinated on the MoS2 monolayer, with the corresponding wave functions overlapping with Au orbitals. A localized state is also apparent in the functionalized structure above the Fermi level corresponding to Au d orbitals. The mentioned changes in the band structure are expected to affect the optical properties of MoS2. In this context, optical absorbance measurements conducted in MoS2, before and after functionalization, provide some valuable information that can be correlated to band structure. Figure 3D depicts the average spectra of 10 pristine (black) and [Au(MoS2)Clx] (pink) samples; some of the individual spectra can be found in fig. S9. After Au coordination, there is a small blue shift in the A (1.83 to 1.88 eV) and B (1.98 to 2.02 eV) exciton absorption peaks, as predicted by the DFT calculations (blue arrow in Fig. 3C). The Tauc plots and its representative fitting of absorption edge show that the MoS2 preserved a direct gap after the functionalization (fig. S10) (24). In addition, a new broad absorption band located at a low energy, between 1.5 and 1.7 eV, can be observed in the Au-MoS2 absorption spectra, which can be assigned to charge transfer transitions from localized states in the Au ions to the conduction band, with major contributions from MoS2 electronic states (red arrow in Fig. 3C). Another charge transfer transition can be seen at 1.35 eV, which can be assigned to an electronic transition from states on the valence band showing a mixture of Au and MoS2 to the conduction band (green) arrow in Fig. 3C. It is worth mentioning that the band structure calculated for [Au(MoS2)Cl] (fig. S11) could also explain the electronic transitions observed experimentally. Since both [Au(MoS2)Cl] and [Au(MoS2)Cl3] lead to p-type doping, both complexes could be well formed during the functionalization process; however, the [Au(MoS2)Cl3] has a better fit with the experimental data because of the higher binding energy and p-type doping induced on MoS2.

Fermi level tuning through AuClx coordination

To further investigate the influence of AuClx complex coordination on the electronic properties of monolayer MoS2, we designed and fabricated back-gated FETs (BGFETs) (fig. S12) Figure 4A shows the schematic of the monolayer MoS2 BGFET. Figure 4 (B and C) indicates the transfer characteristics, i.e., drain current (IDS) versus back-gate voltage (VBG) and output characteristics, i.e., drain current (IDS) versus drain bias (VDS), respectively, of a representative pristine MoS2 BGFET. The field-effect mobility (μFE) extracted from the peak transconductance was found to be ~25 cm2/V·s, which is comparable to the mobility values reported for exfoliated single-crystal material (25). This confirms the superior quality of the as-grown monolayer MoS2. The ON state current, which is proportional to the μFE and the charge in the channel, i.e., COX*(VBGVT), where VT is the threshold voltage of the device and COX is the oxide capacitance, was found to be ~120 μA/μm at VDS = 8 V and VBG = 10 V, further confirming the high performance of the monolayer MoS2 BGFET. The subthreshold slope (SS) of the device was found to be 350 mV per decade, which indicates a reasonably clean interface with the back gate despite the wet transfer process carried out on the MoS2 during device fabrication. Figure 4D shows device-to-device variations in pristine MoS2 BGFETs. These same devices were then treated with ethanol solution with various concentration of HAuCl4 such as 10−10, 10−9, and 10−6 M. Figure 4 (E to G) display the transfer characteristics at VDS = 1 V for 10 representative devices after functionalization. Note that MoS2 BGFETs were treated with the lowest concentration of HAuCl4 at first (Fig. 4E), measured, and then subsequently treated with increasing concentrations (Fig. 4, F and G). For higher HAuCl4 concentrations, the threshold voltage becomes increasingly positive, i.e., it shifts toward the right, indicating that the Au coordination acts as a p-dopant. From TEM, XPS, and DFT results, we can then infer that the AuClx complexes coordinate to the S atoms, forming a fixed charge on the channel, which influences the threshold voltage, based on the equation given belowΔVT=ΔQFcox(1)

Fig. 4 Fermi level tuning using gold doping.

(A) Schematic of a BGFET involving coordination of gold on the MoS2 channel. (B) Back-gated transfer characteristics of the MoS2 FET. (C) Output characteristics of the MoS2 FET. (D to G) Drain current (IDS) versus back-gate voltage (VBG) characteristics with a drain bias of 1 V. (D) Pristine, (E) 10−10 M, (F) 10−9 M, and (G) 10−6 M concentrations of HAuCl4 in ethanol.

Here, ΔVT corresponds to the change in threshold voltage, while ΔQF is the change in the fixed charge. As the fixed charge on the semiconducting material changes, the flat band voltage increases, meaning that the Fermi level moves closer to the valence band. Table 1 shows the statistics for threshold voltage, SS, and field-effect mobilities for 10 representative devices corresponding to different concentrations of HAuCl4. The threshold voltage shift in the devices is extracted for the electron branch using the iso-current method for a current value of 100 nA/μm for VD = 1 V. As the concentration of HAuCl4 increases, the coordination of AuClx complexes to the S also increases, which leads to an increased fixed charge on the channel and, ultimately, a shift in the threshold voltage. Figure S13A shows the error bar plot of threshold voltage variations for different doping concentrations. Mobility values and SS from the statistics suggest that neither the ON-state nor the OFF-state device performance are significantly degraded because of the Au doping. It is also observed that mobility and SS values remain almost constant throughout for various Au concentrations. On the basis of the threshold voltage shift, we further computed the number of AuClx complexes coordinated to the MoS2 monolayer using the equation belown=coxΔVTq(2)

Table 1 Threshold voltage, SS, mobility, and ON/OFF ratio statistics of the 10 representative devices for various concentrations of HAuCl4.

View this table:

Here, q is the electronic charge. Figure S13B shows the error bar plot of the number of coordinated AuClx complexes versus concentration for an area of 10 nm by 10 nm. Note that this particular area was chosen to compare the number of Au atoms coordinated to S through Eq. 2 and then compared with the number obtained from HR-STEM. The results obtained for the 10−6 M concentration (ca. eight single complexes per 100 nm2) match perfectly with the number of Au atoms coordinated to S shown in the HR-STEM image in Fig. 1B.

The stability of the single AuClx complex functionalization is further confirmed via FET measurements carried out 1 year after the functionalization process (fig. S15). It was observed that the devices did not degrade over time, ensuring the devices and the functionalization method to be highly reliable. Therefore, we can conclude that our Au functionalization technique clearly results in a precise tuning of the Fermi level positions by adjusting the HAuCl4 solution concentration.

Thermal boundary conductance measurements

To further understand the effects of the single AuClx complex functionalization of MoS2, we carried out thermal boundary conductance measurements via time-domain thermoreflectance (TDTR). Previous measurements via frequency-domain thermoreflectance have demonstrated that the thermal conductance at graphene contacts can be extracted via a multi-modulation frequency approach (26). Recently, Brown et al. (27) have shown that thermal conductance values at the metal/MoSe2 interfaces can be extracted via a TDTR mapping technique. To capture the appropriate conductance at the Al/MoS2/SiO2 interface (see Materials and Methods), we have used the magnitude of the thermoreflectance signal to locate single crystals of MoS2 and subsequently performed full TDTR measurements near the center of the crystals. Figure 5A depicts a representative micrograph where the thermoreflectance overall magnitude was used to find a MoS2 single crystal. The uniformity of the thermoreflectance magnitude in this region suggests that the thermal conductance across the interface is relatively uniform. The TDTR curve and best fit for the conductance, hK, are presented in Fig. 5B for the pristine MoS2 monolayer, with the inset showing the early picosecond acoustic response used to extract the thickness of the Al layer (28, 29). A summary of the hK results obtained for AuClx functionalized MoS2 monolayers at various metal concentrations is presented in Fig. 5C.

Fig. 5 Thermal conductance enhancement from single AuClx complexes.

(A) TDTR magnitude mapping of a pristine single-crystal MoS2 flake. (B) TDTR model and best fit for the conductance at the Al/pristine MoS2/SiO2 interface. The inset shows the picosecond acoustics response at earlier time delays. (C) Results for the thermal boundary conductance at Al/Au-MoS2/SiO2 interfaces.

In general, we observed that the thermal conductance at the Al/MoS2/SiO2 interfaces is commensurate with the Au concentration. In this context, chemical functionalization of graphene via oxygen plasma treatment has been previously reported to increase the conductance by 50 and 100% at these interfaces when Au and Al are chosen as the metal contact, respectively (30, 31). As in these works, we attribute the enhanced conductance at the Al/MoS2/SiO2 interface to the enhanced reactivity of the MoS2 caused by chemical functionalization. As Al is required to be in contact with some fraction of Au ions, the reasons for the enhanced conductance could be attributed to the additional pathway of conduction via the electronic system or through an improvement of the bonding state at the Al/MoS2 interfacial region through functionalization. In addition, functionalization of MoS2 monolayers will inherently change the localized vibrational density of states, thus offering additional pathways of conduction. Various reports have previously examined the thermal conductance at the MoS2/SiO2 interfaces via Raman spectroscopy (32, 33). There, a conductance value of ~15 MW m−2 K−1 was extracted, which is similar to those presented here. However, a direct comparison between the two approaches is not possible, as the conductance at the Al/MoS2/SiO2 interface is inherently measured in our experiment, whereas just the MoS2/SiO2 interface without a metal topcoat is measured in the referred studies (30, 31). The major enhancement on the thermal boundary conductance of monolayer MoS2, obtained upon Au functionalization, can be applied to increase the heat dissipation rate in vertically stacked 2D transistors.


The formation of chemical bonds is crucial for an effective and nondestructive functionalization of 2D TMD materials. In this work, we demonstrated the formation of the [Au(MoS2)Clx] complex, which results in the transfer of electrons from MoS2 valence band formed by Mo 4d and S 3p orbitals into the Au valence orbitals (6s and 5d), as indicated by the high mixing of Mo, S, and Au states on the bands near the Fermi level shown in Fig. 3. In the STEM images shown in Fig. 1, we observe that the distribution of Au single atoms is fixed on top of S sites rather than randomly dispersed. The movie (Supplementary Materials) and trajectory image (Fig. 1G) further confirm that Au atoms prefer to sit on S sites even when moving. The displacement distance analysis also concludes that the Au atoms stay static on S and move immediately to adjacent S when found on top of Mo. These pieces of evidence confirm that the Au atoms have a much higher affinity toward S atoms, indicating chemical interactions between Au and S atoms. STEM simulations conducted for [Au(MoS2)], [Au(MoS2)Cl], and [Au(MoS2)Cl3] structures (fig. S3) indicate the absence of Cl atoms in our experimental STEM images. It is due to the e-beam irradiation at 80 keV, which might have displaced Cl atoms before image acquisition. The [Au(MoS2)Clx] complex is further studied via XPS. The shift of S 2p spectra (Fig. 2B) that indicate electron transfer from S to Au agrees with the formation of a covalent bond through Au coordination with S atoms, where MoS2 acts as ligands donating electrons to Au atoms. XPS data for the Au and Cl core level depicted in fig. S5 show the presence of Cl atoms on the [Au(MoS2)Clx] sample with a Au:Cl ratio of 1:3. We constructed three different structures with isolated Au ion, Au with one Cl, and Au with three Cl on MoS2 via DFT calculations. For [Au(MoS2)Cl] and [Au(MoS2)Cl3] complexes, the electrons are transferred from the MoS2 to Au, promoting the Au reduction and a p-type doping on the MoS2, thus matching our experimental observations. In contrast, for [Au(MoS2)], electrons are transferred from Au atoms to MoS2, leading to n-type doping. When considering 3Cl as additional ligand, the Au─S binding energy (−1.4 eV) is higher than the Au─S binding energy (−0.86 eV) without Cl, which could facilitate the formation of single AuClx complexes and prevent aggregation. Thus, the simulations indicate the need of Cl ligands to complete gold’s coordination sphere and balance the charges.

The Au 4f spectrum also reveals two states of Au present on the MoS2 surface, namely, Au1+ and Au3+. We believe that the presence of Au3+ comes from the unreacted HAuCl4 precursor, while the Au1+ comes from the reduction with MoS2. The standard reduction potential difference between the Au3+/Au1+ pair (AuCl4/AuCl2 = 0.926 V) and MoS2 (−0.090 V) causes the spontaneous reduction of Au3+ to Au1+ by MoS2. In addition and according to Pearson hard-soft acid and base theory, the reduction from Au3+ to Au1+ makes Au atoms softer, thus promoting a strong bonding to S with a more covalent character since S atoms are soft ligands known to stabilize linear Au1+ complexes. Because of the filled 3p S orbitals, MoS2 acts as a pi-donor ligand. This way, when Au1+ ion forms coordination bonds with MoS2, the Au redox potential should decrease considerably, preventing further reduction into its elemental state. For example, [Au(SCN)2] complexes have a redox potential of +0.66 V, whereas [AuCl2] complexes have a much higher redox potential (E0 = +1.15 V), as shown in fig. S14. The absence of elemental Au atoms also prevents the formation of Au clusters and particles since both Au3+ and Au1+ are positively charged ions that show electrostatic repulsion. The repulsion and coordination bonding thus resulted in the stable Au single ions that we see in Fig. 1.

Furthermore, we believe that this coordination-based functionalization strategy is general and could be applied to other metals and TMD systems. Many transition metals such as Fe, Ni, and Pt have empty d orbitals that have been well known for forming coordination complexes with S-based ligands (3436). Since S is considered a soft base in Pearson hard-soft acid and base theory, we would expect soft atoms such as Pt and Ni to form stronger coordination bonds when compared to Fe and Ti. The coordination sphere will also depend on the ligands from the precursors. In addition, since Se and Te are even softer than S, we believe that the interaction between Au and Ag with selenides or tellurides would be stronger. The controlled doping of Se and Te into MoS2 systems may offer a way to control the location of Au single atoms since Au─Se bonds should be more favorable (37).

In this work, we successfully prepared [Au(MoS2)Clx] and Ag-MoS2 coordination complexes where Au and Ag single atoms were bonded to S atoms via coordination bonds. This approach does not rely on defects, such as vacancies, in the MoS2 crystalline layers but does have a significant impact on the optical, electrical, and thermal properties of the functionalized monolayer MoS2. The formation of the coordination complex [Au(MoS2)Clx] led to the transfer of electrons from MoS2 to the noble metal, which introduced p-type doping to the functionalized atomically thin semiconductor. Moreover, the degree of p-type doping can be fine-tuned by varying the Au precursor (HAuCl4) concentrations, thus controlling the Fermi level of MoS2 and the exciton-to-trion relative population. Furthermore, the MoS2 surface functionalization with single noble metal atoms creates midgap states in the electronic structure of MoS2 and yields band displacements near the Fermi level, which effectively modified the absorbance spectra in the visible range and are responsible for shifts in the PL excitonic emissions. The single noble metal atom functionalization has resulted in a major enhancement of the thermal boundary conductance across monolayer MoS2 interfaces, which could be applied to increase the heat dissipation rate of vertically stacked 2D transistors based on this and other semiconducting TMDs. In addition, the synthesis of single atoms introduced in this work could also be exploited in other applications such as single-atom catalysis, quantum information devices, optoelectronics, and enhanced sensing.


Synthesis of monolayer MoS2

Monolayer MoS2 was synthesized by a salt-assisted CVD method, similar to previous publications (38, 39). NaBr (Alfa Aesar, 99%) was ground into a fine powder with a mortar and pestle and then mixed with MoO2 in a 10:1 ratio by weight. For growth of monolayer films, 2 mg of the salt/oxide mixture was placed at the bottom of a ceramic boat and a piece of SiO2 (300 nm)/Si substrate was placed facing down over the mixture, with ca. 2 mm of space between the mixture and substrate. Sulfur powder (100 mg; Alfa Aesar, 95%, 300 mesh) was used as the sulfur source. The growth substrate was placed at the center of a 1-inch-diameter horizontal tube furnace (Lindberg/Blue M), while the sulfur powder was placed upstream outside of the furnace and 30 cm away from the growth substrates. Before the synthesis experiments, the tube was flushed with 400 standard cubic centimeter per minute (sccm) of Ar for 20 min, and then, the flow was reduced to 100 sccm. The furnace was then heated to 800°C for 20 min and then held for 5 min. The sulfur powder was separately heated to 220°C for 5 min and held for 5 min, while the furnace was kept at 800°C.

Synthesis of the [Au(MoS2)Clx] complex

The functionalization was performed by dipping the Si/SiO2 substrate with the CVD MoS2 into an ethanol solution of HAuCl4 with concentrations between 1 × 10−2 to 1 × 10−9 M for 10 min. The functionalized [Au(MoS2)Clx] sample was then immersed in isopropanol (IPA) for a few seconds to remove excess HAuCl4, followed by N2 drying. Last, the sample was kept in vacuum for 10 min before the measurements. To exclude the effect of ethanol or IPA adsorption on the flake surface, pristine samples were also washed with ethanol and IPA, dried with N2, and kept in vacuum for 10 min.

Scanning transmission electron microscopy

STEM imaging was carried out in a FEI Titan3 G2 60/300 operated at 80 kV to reduce irradiation damage. A high-angle ADF (HAADF) detector was used to collect ADF signal. A Gaussian blur filter was applied using the ImageJ software to reduce the noise and enhance the visibility of the detailed structure, but raw images were used for acquiring the line profile of the ADF intensity. STEM-ADF image simulations were conducted using the QSTEM package. Simulation parameters such as acceleration voltage, spherical aberration (C3 and C5), convergence angle, and inner/outer angle for the HAADF detector were set according to experimental conditions. Note that before the TEM imaging, the CVD-grown MoS2 was first transferred to a TEM grid and then functionalized by 1 × 10−6 M HAuCl4 in ethanol solution.

Device fabrication

BGFETs were fabricated using a 50-nm alumina (Al2O3) film as the gate oxide and a stack of Pt/TiN/p++Si as the back gate electrode. First, CVD-grown MoS2 was transferred onto the alumina sample; then, the sample was spin coated with A6 polymethyl methacrylate (PMMA), followed by e-beam lithography to specify the 2-μm channels and then separating them out by sulfur hexafluoride (SF6) etch under 5°C for 30 s. After each step, the sample was rinsed in acetone for 30 min, followed by a rinse in IPA. To define the source and drain contacts, the sample was then spin coated with methyl methacrylate, followed by A3 PMMA. Using e-beam lithography, source and drain contacts were then patterned and further developed by using 1:1 mixture of 4-methyl-2-pentanone and IPA for 60 s. Forty nanometers of Ni and 30 nm of Au were deposited onto the patterns using e-beam evaporation. Lift-off of the evaporated materials was carried out by immersing the sample in acetone for 30 min, followed by an IPA final rinse.

Electrical characterization

Electrical characterization of the fabricated devices and coordinated devices was performed in a Lake Shore CRX-VF probe station under high vacuum condition using a Keysight B1500A parameter analyzer.

Thermal boundary conductance measurements

A nominal 80-nm Al film was deposited on specimens via e-beam evaporation to serve as the metal transducer in our experiment. The 808.5-nm wavelength beam of a 80-MHz Ti:sapphire oscillator was spectrally separated into pump and probe paths. The pump path was amplitude-modulated via an electro-optic modulator at 8.4 MHz, creating a frequency-dependent heating event at the surface of the Al-coated MoS2. The probe was mechanically delayed in time and monitored the thermoreflectance at the sample surface due to the pump-induced heating event. With a 10× objective, the 1/e2 diameters of the pump and probe were 14 and 11 μm, respectively. We fitted the data to the radially symmetric heat diffusion equation for the conductance, hK, of the Al/MoS2/SiO2 interfaces and the thermal conductivity of the underlying thermal oxide. Additional information regarding the specific analyses of TDTR can be found elsewhere (4042). Single-crystal MoS2 flakes were located by mounting our samples on a two-axis stage, which was oriented to minimize changes in the pump/probe radii over the measurement range.

XPS measurements

XPS measurements were conducted in a high-resolution Thermo Scientific XPS with monochromatic Al Kα x-ray source. The binding energies were calibrated with C1s binding energy of 284.8 eV. The analysis of peak fitting was performed on XPSPEAK41 software.

Optical absorbance measurements

Optical absorption measurements were performed using a homemade setup. For excitation, a supercontinuum white light source (NKT Photonics), ranging from 400 to 2400 nm, was coupled to an optical fiber with core diameter of 9 μm and the output was collimated using a 5× Newport objective [0.10 numerical aperture (NA)]. The beam was then reflected using a Thorlabs visible nonpolarizing cube beam splitter (400 to 700 nm) to a Nikon TU Plan Fluor EPI 50× objective (0.8 NA), which focused the beam on the sample at normal incidence. The reflected signal was collected and collimated by the same 50× objective. After transmission through the cube beam splitter, the beam was focused onto the input of a Thorlabs multimode fiber with core diameter of 500 μm using a 10× Newport objective (0.25 NA), connected to a Yokogawa Optical Spectrum Analyzer operation in the 350- to 1700-nm range. In this setup, the position where the spot was focused on the sample was determined by directing the reflected beam to a Thorlabs USB charge-coupled device camera (1024 × 768 resolution) on which an image was formed using a Thorlabs infinity-corrected lens (focus, 200 mm). For imaging, another visible beam splitter was placed before the camera lens and a diffused white light light-emitting diode source was used to illuminate the sample. To obtain absorbance spectra, the bare substrate was used as a reference.

First-principles calculations

Density functional calculations were performed using the SIESTA code (43) with a DZP localized basis set, energy shift of 0.03 eV, and mesh cutoff of 400 rydberg. The exchange-correlation functional was described by the gradient generalized approximation with a Perdew-Burke-Ernzerhof parameterization (44). We used a 4 × 4 diagonal supercell for the electronic band structure and 3 × 2 orthorhombic supercells for real-space charge density plots. Brillouin zone sampling was based on Γ-centered Monkhorst-Pack method with 4 × 4 grid for supercells (45). A vacuum spacing of 20 Å was added to avoid interaction among periodic layers. Structural optimization was performed with a force tolerance of 0.010 eV/Å.


Supplementary material for this article is available at

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: Funding: This work was supported by the Air Force Office of Scientific Research (AFOSR) through grant no. FA9550-18-1-0072. L.S. acknowledges financial support from CNPq (grant no. 408525/2018-5) and high-performance computing facilities from LoboC/NACAD/UFRJ. This work was partly supported by the São Paulo State Foundation (FAPESP; grant nos. 2012/50259-8, 2015/11779-4 and 2017/01817-1), the Brazilian Nanocarbon Institute of Science and Technology (INCT/Nanocarbono), Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), and CAPES-PRINT (Programa Institucional de Internacionalização; grant no. 88887.310281/2018-00). This work has been funded by the Air Force Office of Scientific Research under award no. FA9550-18-1-0352. Author contributions: The manuscript was written through contributions of all the authors. H.L. and D.G. conducted the functionalization experiments and characterizations. D.G. and C.J.S.d.M. preformed the microabsorption measurements. L.S., C.M.M., and D.G. performed the first-principles calculations. A.D. and S.D. helped with device fabrication and tests. E.K. and R.B.N.B. synthesized the pristine MoS2 materials. K.F., T.Z., and F.Z. helped with TEM characterization and simulation. Y.L., A.L.E., and Y.-T.Y. provided guidance and helped analyze the data. R.C.S. helped with XPS measurements and analysis. M.T., S.D., and C.J.S.d.M. supervised the whole work. All authors gave approval to the final version of the manuscript. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.

Stay Connected to Science Advances

Navigate This Article